Consider two regions situated so far from each other that they are not causally connected. Thus the regions take generally speaking independent configurations. Therefore a so-called topological defect, analogous to a dislocation in ferromagnetic crystals, appears on the boundary between the regions. A simplest type of such a defect is analogous to a point dislocation. In typical Grand Unification Theories (GUT) such objet represents the magnetic monopole. It behaves as a particle with mass
$$
m_{mon}\approx \frac{M_X}{\alpha_U}.
$$
Here ${M_X}\approx 10^{16}\,GeV$ is the GUT energy scale and $\alpha_U\approx 1/40$ is the effective coupling constant. The magnetic monopoles could be created in the hot Universe at the phase transition connected to spontaneous symmetry breaking: when temperature of the Universe falls lower than $T_c\approx E_{GUT}\approx 10^{16}\,GeV,$ the Higgs field presented by $X$ and $Y$ massive bosons acquires non-zero vacuum mean.
Due to their giant mass, the monopoles become nonrelativistic component of the energy density in the Universe right after their creation. The monopole density is expected to roughly equal to unity in each isolated region. Size of such region is determined by the distance passed by light during the time period $t_c$ from the Big Bang to the phase transition. This distance simply equals to the particle horizon at time $t_c.$ If the universe was dominated by radiation up to that time, then $a\sim t^{1/2}$ and therefore the particle horizon equals to $a(t)\int_0^{t_c}\frac{dt'}{a(t')}=2t_c$. Then the predicted monopole number density reads
$$
n_{mon}\approx \frac{1}{\left( 2t_c \right)^3}
$$
The time period ${t_c}$ can be estimated from the relation
$$
t_c=\frac 14\sqrt{\frac{45}{\pi^3g^*}} M_{Pl}T_{c}^{-2}\approx 10^{-39}\, sec.
$$
As the monopoles are non-relativistic particles, then their contribution into the energy density is
$$
\rho_{mon}=n_{mon}m_{mon}\approx \frac{1}{\left( 2t_c \right)^3}\frac{M_X}{\alpha_U}
\approx 2\times 10^{57}\,GeV^4.
$$
Compare this value to the photon concentration at the same time
$$
\rho_{\gamma }=\frac{\pi^2}{15}T_c^4\approx 2\times 10^{63}\,GeV^4.
$$