# Non-topological soliton

In quantum field theory, a non-topological soliton (NTS) is a soliton field configuration possessing, contrary to a topological one, a conserved Noether charge and stable against transformation into usual particles of this field for the following reason. For fixed charge Q, the mass sum of Q free particles exceeds the energy (mass) of the NTS so that the latter is energetically favorable to exist.
The interior region of an NTS is occupied by vacuum different from the ambient vacuum. The vacuums are separated by the surface of the NTS representing a domain wall configuration (topological defect), which also appears in field theories with broken discrete symmetry. Infinite domain walls contradict cosmology, but the surface of an NTS is closed and finite, so its existence would not be contradictory. If the topological domain wall is closed, it shrinks because of wall tension; however, due to the structure of the NTS surface, it does not shrink since the decrease of the NTS volume would increase its energy.

## Introduction

Quantum field theory has been developed to predict the scattering probability of elementary particles. However, in the mid 1970s it was found out that this theory predicts one more class of stable compact objects: non-topological solitons (NTS). The NTS represents an unusual coherent state of matter, called also bulk matter. Models were suggested for the NTS to exist in forms of stars, quasars, the dark matter and nuclear matter.
A NTS configuration is the lowest energy solution of classical equations of motion possessing a spherical symmetry. Such a solution has been found for a rich variety of field Lagrangians. One can associate the conserved charge with global, local, Abelian and non-Abelian symmetry. It appears to be possible that the NTS configuration exists with both bosons as well as with fermions. In different models either one and the same field carries the charge and binds the NTS, or there are two different fields: charge carrier and binding field.
The spatial size of the NTS configuration may be elementary small or astronomically large, depending on the model fields and constants. The NTS size could increase with its energy until the gravitation complicates its behavior and finally causes the collapse. In some models, the NTS charge is bounded by the stability (or metastability) condition.

## Simple examples
### One field

For a complex scalar field with the U(1) invariant Lagrange density
    ${\mathcal {L}}=|\partial _{\mu }\Phi |^{2}-U(|\Phi |)\,$
the NTS is a ball with radius R filled with the field $\Phi =(\phi _{0}/{\sqrt {2}})e^{i\omega t}$. Here $\phi _{0}$ is a constant inside the ball except for a thin surface coat where it sharply drops to the global U(1) symmetrical minimum of $U(|\Phi |)$. The value $\phi _{0}$ is adjusted so that it minimises the energy of the configuration
    $E={\Big{\frac {\phi _{0}^{2}\omega ^{2}}{2}}+U({\frac {\phi _{0}}{\sqrt {2}}}){\Big )}{\frac {4}{3}}\pi R^{3}.\,\,\,\,\,\,\,\,\,\,\,\,\,(1)$
Since the U(1) symmetry gives the conserved current $j_{\mu }=-i(\Phi ^{*}\partial _{\mu }\Phi -\partial _{\mu }\Phi ^{*}\Phi ),\,$
the ball possesses the conserved charge
    $Q=\int j_{0}d^{3}x=\omega \varphi _{0}^{2}{\frac {4}{3}}\pi R^{3}.$
The minimization of the energy (1) with R gives
    $E=Q{\sqrt {\frac {2U(\varphi _{0}/{\sqrt {2}})}{\phi _{0}^{2}}}}.\,\,\,\,\,\,\,\,\,\,\,\,(2)$
The charge conservation allows the decay of the ball into Q particles exactly. This decay is energetically unprofitable if the sum mass Qm exceed the energy (2). Therefore, for the NTS existence it is necessary to have
    ${\frac {2U(\varphi _{0}/{\sqrt {2}})}{\varphi _{0}^{2}}}\equiv {\frac {U(|\Phi |)}{|\Phi |^{2}}}{\Bigg |}_{\min }<m^{2}.$
The thin wall approximation, which was used above, allows to omit the gradient term $E_{\text{surface}}$ in the expression for energy (1), since $E_{\text{surface}}\sim U(\varphi _{0}/{\sqrt {2}})R^{2}\Delta R\ll U(\varphi _{0}/{\sqrt {2}})R^{3}$. This approximation is valid for $Q\gg 1$ and is justified by the exact solution of the motion equation.

### Two fields

The NTS configuration for a couple of interacting scalar fields is sketched here. The Lagrange density
    ${\mathcal {L}}=|\partial _{\mu }\Phi |^{2}+{\frac {1}{2}}(\partial _{\mu }\sigma )^{2}-g|\Phi |^{2}\sigma ^{2}-{\frac {\lambda }{4}}(\sigma ^{2}-\sigma _{0}^{2})^{2}$
is invariant under U(1) transformation of the complex scalar field $\Phi .$ Let this field depends on time and coordinate simply as $\Phi ({\vec {r}},t)=(\phi (r)/{\sqrt {2}})e^{i\omega t}$. It carries the conserved charge $Q=4\pi \omega \int \limits _{0}^{\infty }\phi ^{2}(r)r^{2}dr$. In order to check that the energy of the configuration is smaller than Qm, one should either to calculate this energy numerically or to use the variational method. For trial functions $\phi (r)=\varphi _{0}(\sin \omega r/r)$ and $\sigma (r)=0$ for r < R,
    $\phi (r)=0{\text{ and }}\sigma (r)=\sigma _{0}{\Big1-\exp(-(r-R)/l){\Big )}{\text{ for }}r>R,\,$
the energy in the large Q limit is approximately equal to $E\approx {\frac {\pi Q}{R}}+{\frac {\pi }{3}}\lambda \sigma _{0}^{4}R^{3}$.
The minimization with R gives the upper estimation $E_{up}={\frac {4}{3}}\pi \lambda ^{1/4}Q^{3/4}\sigma _{0}$
for the energy of the exact solution of motion equations $\omega ^{2}\phi +{\vec {\partial ^{2}}}\phi -g\sigma ^{2}\phi =0$ and ${\vec {\partial ^{2}}}\sigma -g\phi ^{2}\sigma -\lambda \sigma (\sigma ^{2}-\sigma _{0}^{2})=0$.
It is indeed smaller than $Qg\sigma _{0}$ for Q exceeding the crucial charge
    $Q_{\min }={{\Big{\frac {4\pi }{3g}}{\Big )}}^{4}\lambda .$

### Fermion plus scalar

If instead of boson, fermions carry the conserved charge, an NTS also exists. At this time one could take
    ${\mathcal {L}}=\sum _{k=1}^{N}\left({\frac {i}{2}}{\overleftrightarrow {{\overline {\Psi }}_{k}\gamma ^{\mu }\partial _{\mu }\Psi _{k}}}-(m+g\sigma ){\overline {\Psi }}_{k}\Psi _{k}\right)-U(\sigma )+{\frac {1}{2}}(\partial _{\mu }\sigma )^{2}.\,\,\,\,\,(3)$
N is the number of fermion species in the theory. Q can't exceed N due to the Pauli exclusive principle if the fermions are in the coherent state. This time the NTS energy E is bound by
    $E_{\text{up}}={\frac {4}{3}}\pi {\sqrt {2}}U^{1/4}\left(-{\frac {m}{g}}\right)N^{3/4},$
See Friedberg/Lee.

## Stability
### Classical stability

The condition $E(Q)<Qm$ only allows to assert the NTS stability against a decay into free particles. The equation of motion gives $E(Q)$ only on a classical level. At least two things should be taken into account: (i) the decay into smaller pieces (fission) and (ii) the quantum correction for $E(Q)$.
The condition of stability against the fission looks as follows:
    ${\frac {d^{2}E(Q)}{dQ^{2}}}<0.\,\,\,\,\,\,\,\,\,\,\,\,\,(4)$
It signifies that $E(Q_{1})+E(Q_{2})>E(Q_{1}+Q_{2})$. This condition is satisfied for the NTS in examples 2.2 and 2.3. The NTS in example 2.1, called also Q-ball, is stable against the fission as well, even though the energy (2) does not satisfy (4): one has to recollect the omitted gradient surface energy and to add it to the Q-ball energy (1). Perturbatively, $E_{\text{surface}}\propto R^{2}(Q)\propto Q^{2/3}$. Thus
    ${\frac {d^{2}(E_{\text{surface}}+Q{\sqrt {2U(\phi _{0}/{\sqrt {2}})/\phi _{0}^{2}}})}{dQ^{2}}}<0.$
Another job $E_{\text{surface}}$ does, is to set $Q_{\min }$ for the thin-wall description of Q-ball: for small Q the surface becomes thicker, $E_{\text{surface}}$ grows and kills the energy gain $Q(m-{\sqrt {U(|\Phi |)/|\Phi |^{2}}}{\Big |}_{\min })$. However the formalism for the thick-wall approximation has been developed by Kusenko who says that for small Q, NTS also exists.

### Quantum correction

As for quantum correction, it also diminishes the binding energy per charge $m-E(Q)/Q$ for small NTS, making them unstable. The small NTS are especially important for the fermion case, since it is naturally to expect rather small number of fermions species N in (3), and consequently, Q. For Q=2 the quantum correction decreases the binding energy by 23%. For Q=1 a calculation based on the path integral method has been carried out by Baacke. The quantum energy has been derived as a time derivative of the one-loop fermion effective action
    $S_{\text{eff}}=-i\ln \det {\Big{\frac {i\gamma ^{\mu }{\partial }_{\mu }-g\sigma (x)}{i\gamma ^{\mu }{\partial }_{\mu }-g\sigma _{0}}}{\Big )}.$
This calculation gives the loop energy of the order of binding energy. In order to find the quantum correction following the canonical method of quantization, one has to solve the Schrödinger equation for the Hamiltonian built with quantum expansion of field functions. For the boson field NTS it reads
    $\Phi =1/{\sqrt {2}}(\phi _{cl}({\vec {r}}-{\vec {R(t)}})+\sum _{n}q_{n}(t)\beta ({\vec {r}}))e^{i\zeta (t)},\sigma =\sigma _{cl}({\vec {r}}-{\vec {R(t)}})+\sum _{n}q_{n}(t)\alpha ({\vec {r}}).$
Here $\phi _{cl}\,$ and $\sigma _{cl}$ are the solutions of the classical equation of motion, ${\vec {R(t)}}$ represents motion of the mass center, $\zeta (t)$ is the over-all phase, $q_{n},n=5,6,\dots ,\infty$ are the vibration coordinates, by analogy with the oscillator decomposition of photon field
    $A^{\mu }=\sum _{k}a_{k}^{\mu }(t)e^{i{\vec {k}}{\vec {r}}}+c.c.$
For this calculation the smallness of four-interaction constant is essential, since the Hamiltonian is taken in the lowest order of that constant. The quantum decreasing of the binding energy increases the minimal charge $Q_{\min }$ making the NTS metastable between old and new values of this charge.
NTSs in some models become unstable as Q exceeds some stable charge $Q_{\max }:E(Q>Q_{\max })>Qm$. For example, NTS with fermions carrying a gauge charge has $E(Q)\propto (C_{1}Q^{4/3}+C_{2}Q^{2})^{2/3}$ exceeding Qm for Q large enough as well as for small one. Besides, the gauged NTS probably is unstable against a classical decay without conservation of its charge due to complicated vacuum structure of the theory. Generally, the NTS charge is limited by the gravitational collapse: $E(Q)/M_{Pl}^{2}R(Q)<1$.

### Particle emission

If one adds to the Q-ball Lagrange density an interaction with massless fermion $\Psi$
    ${\mathcal {L}}_{\Psi }+{\mathcal {L}}_{\Psi \Phi }=\Psi ^{+}(i\partial _{0}+i{\vec {\partial }}{\vec {\sigma }})\Psi -ig\Phi \Psi ^{+}\sigma _{2}\Psi ^{*}+ig\Phi ^{*}\Psi ^{T}\sigma _{2}\Psi$
which is also U(1) invariant assuming the global charge for boson twice as for fermion, Q-ball once created begins to emit its charge with $\Psi$-pairs, predominantly from its surface. The evaporation rate per unit area $dN/dtdS<{\BigU(|\Phi |)/|\Phi |^{2}{\Big |}_{\min }{\Big )}^{3/2}/192\pi ^{2}$.
The ball of trapped right-handed Majorana neutrinos in $SU(2)_{L}SU(2)_{R}$ symmetric electroweak theory loses its charge (the number of trapped particles) through the neutrino-antineutrino annihilation by emitting photons from the whole volume.
The third example for a NTS metastable due to particle emission is the gauged non-Abelian NTS. The massive (outside the NTS) member of fermionic multiplet decays into a massless one and a gauged boson also massless in the NTS. Then the massless fermion carries away the charge since it does not interact at all with the Higgs field.
Three last examples represent a class for NTS metastable due to emission of particles which do not participate in the NTS construction. One more similar example: because of the Dirac mass term $m({\overline {\Psi }}_{R}\Psi _{L}+{\overline {\Psi }}_{L}\Psi _{R})$ , right-handed neutrinos convert to left-handed ones. That happens at the surface of neutrino ball mentioned above. Left-handed neutrinos are very heavy inside the ball and they are massless outside it. So they go away carrying the energy and diminishing the number of particles inside. This "leakage" appears to be much slower than the annihilation onto photons.

## Soliton-stars
### Q-star

As the charge Q grows and E(Q) the order of $M_{Pl}^{2}R(Q)$, the gravitation becomes important for NTS. A proper name for such an object is a star. A boson-field Q-star looks like a big Q-ball. The way gravity changes E(Q) dependence is sketched here. It is the gravity what makes $d^{2}E(Q)/dQ^{2}<0$ for Q-star — stabilize it against the fission.
Q-star with fermions has been described by Bahcall/Selipsky. Similar the NTS of Friedberg & Lee, the fermion field carrying a global conserved charge, interacts with a real scalar field.
    ${\mathcal {L}}={\frac {i}{2}}{\overleftrightarrow {{\overline {\Psi }}\gamma ^{\mu }\partial _{\mu }\Psi }}-m(\sigma ){\overline {\Psi }}\Psi -U(\sigma )+{\frac {1}{2}}(\partial _{\mu }\sigma )^{2}.$
The $\sigma$ inside Q-star moves from a global maximum of the potential changing the mass of fermions and making them bound.
But this time Q is not the number of different fermion species but it is the large number of one and the same kind particles in the Fermi gas state. Then for the fermion field description one has to use $k_{F}(r)$ instead of $\Psi (r)$ and the condition of pressure equilibrium instead of the Dirac equation for $\Psi (r)$. Another unknown function is the scalar field $\sigma (r)$ profile which obeys the following motion equation : ${\vec {\partial }}^{2}\sigma -(dm(\sigma )/d\sigma )\langle {\overline {\Psi }}\Psi \rangle -dU(\sigma )/d\sigma =0$. Here$\langle {\overline {\Psi }}\Psi \rangle$ is the scalar density of fermions, averaged on statistical ensemble:
    $\langle {\overline {\Psi }}\Psi \rangle =\int {\frac {m}{(k^{2}+m^{2})^{1/2}}}{\frac {2d^{3}k}{(2\pi )^{3}}},\quad 0<k<k_{F}.$
Fermi energy of the fermion gas $\varepsilon _{F}=(k_{F}^{2}+m^{2})^{1/2}$.
Neglecting the derivatives of $\sigma (r)$ for large Q, that equation together with the pressure equilibrium equation $P_{f}-U=0$, constitute a simple system which gives $k_{F}$ and $\sigma$ inside the NTS. They are constant since we have neglected the derivatives. The fermion pressure
    $P_{f}=\int {\frac {k^{2}}{3(k^{2}+m^{2})^{1/2}}}{\frac {2d^{3}k}{(2\pi )^{3}}},\quad 0<k<k_{F}.$
For example, if $m(\sigma )=g\sigma$ and $U(\sigma )=\lambda /4(\sigma ^{2}-\sigma _{0}^{2})^{2}$, then $\sigma =0$ and $k_{F}=\varepsilon _{F}=(3\pi ^{2}\lambda )^{1/4}\sigma _{0}$. That means fermions appear to be massless in the NTS. Then the full fermion energy $E_{f}=3P_{f}$. For an NTS with the volume $V$ and the charge $Q=V(k_{F}^{3}/3\pi ^{2})$, its energy is proportional to the charge: $E(Q)=(U+E_{f})V=\varepsilon _{F}Q$.
The described above fermion Q-star has been considered as a model for neutron star in the effective hadron field theory.

### Soliton star

If the scalar field potential $U(\sigma )$ has two degenerate or almost degenerate minima, one of them have to be the real (true) minimum in which we happen to leave. Inside NTS $\sigma$ occupies another one. In such a model non-zero vacuum energy appears only at the NTS surface, not in its volume. This allows for the NTS to be very big without falling in gravitational collapse.
That is the case in the left-right symmetric electroweak theory. For a scale of symmetry breaking about 1 TeV, $\nu$-ball of trapped right-handed massless neutrino might have the mass (energy) about 10$^8$ solar masses and was considered as a possible model for quasar.
For the degenerate potential $U(\sigma )=\mu ^{2}\sigma ^{2}(1-\sigma /\sigma _{0})^{2}/2$ both boson and fermion soliton stars were investigated.
A complex scalar field could alone form the state of gravitational equilibrium possessing the astronomically large conserved number of particles. Such objects are called minisoliton stars because of their microscopic size.

## Non-topological soliton with standard fields

Could a system of the Higgs field and some fermion field of the Standard model be in the state of Friedberg & Lee NTS ? That is more possible for a heavy fermion field: for a such one the energy gain would be the most because it does lose its large mass in the NTS interior, were the Yukawa term $g\sigma {\overline {\Psi }}\Psi$ vanishes due to $\sigma \simeq 0$. The more so if the vacuum energy in the NTS interior $U(0)=\lambda \sigma _{0}^{4}/4$ is large, that would mean the large Higgs mass $m_{H}=\sigma _{0}{\sqrt {2\lambda }}$. The large fermion mass implies strong Yukawa coupling $g$.
Calculation shows that the NTS solution is energetically favored over a plane wave (free particle) only if $g>2.3$ for even very small $m_{H}$. For $m_{H}$ =350 GeV (this is the point were $\lambda =1$ for experimentally known $\sigma _{0}=$250 GeV) the coupling $g$ must be more than five.
The next question is whether or not multi-fermion NTS like a fermion Q-star is stable in the Standard model. If we restrict ourself by one fermion species, then the NTS has god the gauge charge. One can estimate the energy of gauged NTS as follows:
    $E=(4\pi )^{1/3}(3/4)^{5/3}Q^{4/3}/R+\beta e^{2}Q^{2}/R+(4\pi R^{3}/3)\lambda \sigma _{0}^{4}/4.$
Here $R$ and $Q$ are its radius and charge, the first term is the kinetic energy of the fermi-gas, the second is the Coulomb energy, $\beta$ takes into account the charge distribution inside the NTS and the latest one gives the volume vacuum energy. Minimization with $R$ gives the NTS energy as a function of its charge:
    $E(Q)=4\sigma _{0}(\pi \lambda )^{1/4}/3\left((4\pi )^{1/3}(3/4)^{5/3}Q^{4/3}+\beta e^{2}Q^{2}\right)^{3/4}.$
An NTS is stable if $E(Q)$ is smaller than the sum of masses for $Q$ particles at infinite distance each from other. That is case for some $Q$, but such a $E(Q)$ dependence allows the fission for any $Q$.
Why could not quarks be bound in a hadron like in NTS. Friedberg and Lee investigated such a possibility. They assumed quarks getting huge masses from their interaction with a scalar field $\sigma$. Thus free quarks are heavy and escape from detection. The NTS built with quarks and $\sigma$ fields demonstrate static properties of hadrons with 15% accuracy. That model demands SU(3) symmetry additional to the color one in order to preserve the later unbroken so that QCD gluons get large masses by SU(3) symmetry breaking outside hadrons and also avoid detection.
Nuclei have been considered as NTS's in the effective theory of strong interaction which is easier to deal with than QCD.

## Solitonogenesis
### Trapped particles

The way NTS's could be born by depends on whether or not the Universe carries a net charge. If it does not then NTS could be formed from random fluctuations of the charge. Those fluctuations grow up, disturb the vacuum and create NTS configurations.
If the net charge is present, i.e. charge asymmetry exists with a parameter $\eta =(|n_{\Phi }-n_{\overline {\Phi }}|)/(n_{\Phi }+n_{\overline {\Phi }})$, NTS could be simply born as the space became divided onto finite regions of true and false vacuum during the phase transition in the early Universe. Those occupied by the NTS (false) vacuum are almost ready NTSs. The scenario of the region formation depends on the phase transition order.
If the first order phase transition occurs, then nucleating bubbles of true vacuum grow and percolate, shrinking regions filled with the false vacuum. The later are preferable for charged particles to live in due to their smaller masses, so those regions become NTSs.
In case of the second order phase transition as temperature drops below the crucial value $T_{c}$ the space consist of interconnecting regions of both false and true vacua with characteristic size $\xi \propto T_{c}^{-1}$. This interconnection "freezes out" as its rate becomes smaller than the expansion rate of the Universe at Ginzburg temperature $T_{G}$, then the regions of two vacua percolate.
But if the false vacuum energy is large enough, $\Lambda >0.8U_{M}$ on the plot, the false vacuum forms finite clusters (NTS's) surrounded by the percolated true vacuum. The trapped charge stabilizes clusters against collapse.
In the second scenario of the NTS formation the number of born $Q$-charged NTS's per unit volume is simply the number density of clusters holding $Q$ particles. Their number density is given by $n(r)=br^{-3/2}\exp {(-cr)}/V_{\xi }$, here b and c are constants of the order of unit, $r\simeq (L/2\xi )^{3}$ is the number of correlation volumes $V_{\xi }$ in a cluster of size $L$. The number of particle in a cluster is $Q(r)\simeq rV_{\xi }\eta n_{Q}$, here $n_{Q}$ is the charge density in the universe at Ginzburg temperature. Thus big clusters are born very rarely and if the minimum stable charge $Q_{\min }$ is present, then overwhelming majority of born NTS carries $Q_{\min }$.
For the following Lagrange density with biased discrete symmetry
    ${\mathcal {L}}=|\partial _{\mu }\Phi |^{2}+(\partial _{\mu }\sigma )^{2}/2-\lambda _{1}(\sigma ^{2}-\sigma _{0}^{2})^{2}/8-\lambda _{2}(\sigma -\sigma _{0})^{3}\sigma _{0}/3-h|\Phi |^{2}(\sigma -\sigma _{0})^{2}-g|\Phi |^{4}-\Lambda$
with
    $\lambda _{1}/\lambda _{2}=0.15,\,$
    $\Lambda =0.6\lambda _{1}\sigma _{0}^{4}$ and $\xi =1/\lambda _{1}T_{G},$
it appears to be $Q_{\min }=18\lambda _{1}/h^{2}$ and
    $n(Q_{\min })\propto (\eta /\lambda _{1}Q_{\min })^{3/2}\exp {(-cQ_{\min }\lambda _{1}^{3}/\eta )}.$

### Field condensate

The net charge could be also placed in the complex scalar field condensate $\langle \Phi \rangle$ instead of free particles. This condensate could consist of spatially homogeneous $\Phi =f(t)\exp {(i\theta (t))}/{\sqrt {2}}$ and $|\Phi |^{2}=f^{2}(t)/2$ provides its potential to be at minimum as the universe cools down and the temperature correction changes the form of the potential. Such a model was treated to explain the baryon asymmetry.
If the field potential allows Q-ball to exist, then they could be born from this condensate as the charge volume density $j_{0}=f^{2}\partial _{0}\theta (t)$ drops in course of the universe expansion and becomes equal to Q-balls charge density. As follows from the equation of motion for $\theta (t)$, this density $j_{0}$ changes with the expansion as the minus third power of scale factor $a(t)$ for the expanding space-time with the differential length element $ds^{2}=dt^{2}-a^{2}(t)(dx^{2}+dy^{2}+dz^{2})$.
Breaking the condensate onto Q-balls appears to be favorable over further dilution of the homogeneous charge density by expansion. The total charge in a comoving volume $Q=\int j_{0}a^{3}d^{3}x$ stays fixed of course.
The condensation of $\Phi$ could occur at high temperature of the universe, due to the negative temperature correction to its mass: $m(T)=m^{2}-\lambda ^{2}/2$ which provides with minimum its potential $U(|\Phi |)=m^{2}|\Phi |^{2}-\lambda |\Phi |^{4}/2+\gamma |\Phi |^{6}$. Here the last term is induced by the interaction $h\chi ^{2}|\Phi |^{2}$ with additional field $\chi$ that has to be introduced in order to satisfy the Q-ball existence condition $U(|\Phi |)/|\Phi |^{2}{\Big |}_{\min }<m^{2}$. At the temperature relevant to relevant Q-balls formation $\chi$ appears only through virtual process (loops) because it is heavy. An alternative way to satisfy the Q=ball existence condition is to appeal to the non-Abelian symmetry.

## Further evolution

Once formed, the NTSs undergo complicated evolution, losing and acquiring the charge by interaction with each other and surrounding particles. Depending on theory parameters, they could either disappear at all or get statistical equilibrium and "freeze out" at some temperature of the universe, or be born "frozen out" if their interaction is slower than expansion rate at $T_{G}$. In the first and the second cases, their up-to-date abundance (if any) has nothing to do with that at the moment of formation.
Since an NTS is a composite object, it has to demonstrate properties different from those of a single particle, e.g. evaporation emission, excitation levels, scattering form-factor. Cosmic observations of such phenomena could provide the unique information about the physics beyond the ability of accelerators.