Every cosmology lecturer these days is confronted with teaching the modern cosmological standard model $\Lambda$CDM, and there are many approaches to do this. However, the danger is imminent that it is presented to students as something set into stone, merely to be accepted as a fact based on the plenty of evidences we have. This is even more critical, given that the standard model of cosmology confronts us with entities not yet fully understood, namely a cosmological constant $\Lambda$ and cold dark matter. In this article, we report on an advanced cosmology course exercise, conducted in computer lab, which was conceived as a means to have students experience first-hand why the $\Lambda$CDM model has become so prevalent in the interpretation of modern cosmological data. To this end, we focused on the cosmic microwave background radiation (CMB) and calculated theoretical temperature and matter power spectra, using the modern Boltzmann code CLASS. By comparing and analyzing the outcome for $\Lambda$CDM, as well as for three other exotic cosmological models, the students were able to grasp the impact of cosmological parameters on CMB observables, and also to understand some of the complicated CMB physics in a direct way. Our chosen examples are not exhaustive and can be easily modified or expanded, so we express the hope that this article will serve as a valuable resource for interested students and lecturers.

We consider an alternative cold dark matter candidate, ultralight bosons ($m>10^{-22}$eV) described by a complex scalar field (SFDM) with global U(1) symmetry, with comoving particle number density conserved after particle production during standard reheating. We allow for repulsive self-interaction. In a Lambda-SFDM universe, SFDM starts relativistic, evolving from stiff (w=1) to radiation-like (w=1/3), becoming nonrelativistic (w=0) at late times. Thus, a stiff-SFDM-dominated era precedes the familiar radiation-dominated era. SFDM particle mass $m$ and quartic self-interaction strength \lambda, are therefore constrained by cosmological observables, N_{eff}, the effective number of neutrino species during BBN, and z_{eq}, the matter-radiation equality redshift. Since the stochastic gravitational wave background (SGWB) from inflation is amplified during the stiff-SFDM-dominated era, it can also contribute a radiationlike component large enough to affect these observables. Remarkably, this amplification makes this SGWB detectable by current GW experiments, e.g., aLIGO/Virgo and LISA, for Lambda-SFDM models satisfying cosmological constraints, for a range of reheat temperatures T_{re} and currently allowed values of tensor-to-scalar ratio $r$. For given r and $\lambda/(mc^2)^2$, the marginally-allowed Lambda-SFDM model for each T_{re} has the smallest m that satisfies cosmological constraints. For example, for marginally-allowed models with r=0.01 and $\lambda/(mc^2)^2=10^{-18}$eV$^{-1}$cm$^3$, null detection by the aLIGO O1 run excludes 8.75*10^3<T_{re} (GeV)<1.7*10^5 at 95% confidence, demonstrating that GW experiments already place a new kind of cosmological constraint on SFDM. A wider parameter range should be accessible to aLIGO/Virgo O5, with potential to detect this signature of Lambda-SFDM. For this same illustrative family, 3-sigma detection is predicted for 600<T_{re} (GeV)<10^7.

We study the impact on electroweak baryogenesis from a swifter cosmological expansion induced by dark matter. We detail the experimental bounds that one can place on models that realize it, and we investigate the modifications of these bounds that result from a non-standard cosmological history. The modifications can be sizeable if the expansion rate of the Universe increases by several orders of magnitude. We illustrate the impact through the example of scalar field dark matter, which can alter the cosmological history enough to enable a strong-enough first-order phase transition in the Standard Model when it is supplemented by a dimension six operator directly modifying the Higgs boson potential. We show that due to the modified cosmological history, electroweak baryogenesis can be realized, while keeping deviations of the triple Higgs coupling below HL-LHC sensitivies. The required scale of new physics to effectuate a strong-enough first order phase transition can change by as much as twenty percent as the expansion rate increases by six orders of magnitude.

Dark Stars are stellar objects made (almost entirely) of hydrogen and helium, but powered by the heat from Dark Matter annihilation, rather than by fusion. They are in hydrostatic and thermal equilibrium, but with an unusual power source. Weakly Interacting Massive Particles (WIMPs), among the best candidates for dark matter, can be their own antimatter and can annihilate inside the star, thereby providing a heat source. Although dark matter constitutes only $\lesssim 0.1\%$ of the stellar mass, this amount is sufficient to power the star for millions to billions of years. Thus, the first phase of stellar evolution in the history of the Universe may have been dark stars. We review how dark stars come into existence, how they grow as long as dark matter fuel persists, and their stellar structure and evolution. The studies were done in two different ways, first assuming polytropic interiors and more recently using the MESA stellar evolution code; the basic results are the same. Dark stars are giant, puffy ($\sim$ 10 AU) and cool (surface temperatures $\sim$10,000 K) objects. We follow the evolution of dark stars from their inception at $\sim 1 M_\odot$ as they accrete mass from their surroundings to become supermassive stars, some even reaching masses $> 10^6 M_\odot$ and luminosities $>10^{10} L_\odot$, making them detectable with the upcoming James Webb Space Telescope. Once the dark matter runs out and the dark star dies, it may collapse to a black hole; thus dark stars may provide seeds for the supermassive black holes observed throughout the Universe and at early times. Other sites for dark star formation may exist in the Universe today in regions of high dark matter density such as the centers of galaxies. The current review briefly discusses dark stars existing today, but focuses on the early generation of dark stars.

The cosmic evolution of a dark matter model which behaves relativistically in the early Universe is explored. Dark matter is described as a complex scalar field, whose earliest evolution is characterized by a stiff equation of state ($p \simeq \rho$). In this phase, it is the dominant component in the Universe. We present constraints from Big Bang nucleosynthesis and primordial gravity waves from inflation. Also, we study how the associated enhanced expansion rate due to the stiff phase might facilitate a first-order electroweak symmetry breaking phase transition, in light of the recently measured value of the Higgs boson mass.

We use the stellar evolution code MESA to study dark stars. Dark stars (DSs), which are powered by dark matter (DM) self-annihilation rather than by nuclear fusion, may be the first stars to form in the Universe. We compute stellar models for accreting DSs with masses up to 10^6 M_{sun}. The heating due to DM annihilation is self-consistently included, assuming extended adiabatic contraction of DM within the minihalos in which DSs form. We find remarkably good overall agreement with previous models, which assumed polytropic interiors. There are some differences in the details, with positive implications for observability. We found that, in the mass range of 10^4 -10^5 M_{sun}, our DSs are hotter by a factor of 1.5 than those in Freese et al.(2010), are smaller in radius by a factor of 0.6, denser by a factor of 3 - 4, and more luminous by a factor of 2. Our models also confirm previous results, according to which supermassive DSs are very well approximated by (n=3)-polytropes. We also perform a first study of dark star pulsations. Our DS models have pulsation modes with timescales ranging from less than a day to more than two years in their rest frames, at z ~ 15, depending on DM particle mass and overtone number. Such pulsations may someday be used to identify bright, cool objects uniquely as DSs; if properly calibrated, they might, in principle, also supply novel standard candles for cosmological studies.

(Abridged) Despite the success of Cold Dark Matter (CDM) in explaining a wide range of observations, the microscopic nature of dark matter is still unknown. Attempts to detect WIMPs, the most commonly studied form, have not yet succeeded. Some predictions are even in apparent conflict with observations (e.g. cuspy-cored halos and `missing satellite problem'). Therefore, it is important to consider different forms of CDM. We focus on dark matter comprised of ultralight bosons that form a Bose-Einstein Condensate (BEC), described by a complex scalar field, and solve the Klein-Gordon and Einstein field equations to evolve the Friedmann-Robertson-Walker universe. We find, in addition to phases of radiation-, matter- and \Lambda-domination, an earlier phase of scalar-field-domination. Moreover, while WIMP CDM is non-relativistic at all times after it decouples, BEC scalar field dark matter (SFDM) is relativistic at early times, evolving from stiff to radiation-like, before it becomes non-relativistic (CDM-like) at late times. The timing of transitions between these phases yields fundamental constraints on SFDM model parameters, particle mass m and self-interaction coupling strength \lambda. We derive the range of particle parameters required to match observations of the evolving background universe, including the CMB and abundances of light elements produced by BBN, characterized by N_eff, the effective number of neutrino species, and the epoch of matter-radiation equality z_eq. This yields m >= 2.4*10^{-21} eV/c^2 and 9.5*10^{-19} eV^{-1}cm^3 <= \lambda/(mc^2)^2 <= 4*10^{-17} eV^{-1}cm^3. Our model accommodates current observations in which N_eff at BBN is higher than at z_eq, as probed by the CMB, otherwise unexplained by WIMP CDM. SFDM without self-interaction (`Fuzzy Dark Matter') cannot satisfy current BBN constraints within 68% confidence and is therefore disfavored.

If cosmological cold dark matter (CDM) consists of light enough bosonic particles that their phase-space density exceeds unity, they will comprise a Bose-Einstein condensate (BEC). The nature of this BEC-CDM as a quantum fluid may then distinguish it dynamically from the standard form of CDM involving a collisionless gas of non-relativistic particles that interact purely gravitationally. We summarize some of the dynamical properties of BEC-CDM that may lead to observable signatures in galactic halos and present some of the bounds on particle mass and self-interaction coupling strength that result from a comparison with observed galaxies.

The nature of the cosmological dark matter remains elusive. Recent studies have advocated the possibility that dark matter could be composed of ultra-light, self-interacting bosons, forming a Bose-Einstein condensate in the very early Universe. We consider models which are charged under a global U(1)-symmetry such that the dark matter number is conserved. It can then be described as a classical complex scalar field which evolves in an expanding Universe. We present a brief review on the bounds on the model parameters from cosmological and galactic observations, along with the properties of galactic halos which result from such a dark matter candidate.

(Abridged) Extensions of the standard model of particle physics predict very light bosons, ranging from about 10^{-5} eV for the QCD axion to 10^{-33} eV for ultra-light particles, which could be the cold dark matter (CDM) in the Universe. If so, their phase-space density must be high enough to form a Bose-Einstein condensate (BEC). The fluid-like nature of BEC-CDM dynamics differs from that of standard collisionless CDM (sCDM), so observations of galactic haloes may distinguish them. sCDM has problems with galaxy observations on small scales, which BEC-CDM may overcome for a large range of particle mass m and self-interaction strength g. For quantum-coherence on galactic scales of radius R and mass M, either the de-Broglie wavelength lambda_deB <~ R, requiring m >~ m_H \cong 10^{-25}(R/100 kpc)^{-1/2}(M/10^{12} M_solar)^{-1/2} eV, or else lambda_deB << R but self-interaction balances gravity, requiring m >> m_H and g >> g_H \cong 2 x 10^{-64} (R/100 kpc)(M/10^{12} M_solar)^{-1} eV cm^3. Here we study the largely-neglected effects of angular momentum. Spin parameters lambda \cong 0.05 are expected from tidal-torquing by large-scale structure, just as for sCDM. Since lab BECs develop quantum vortices if rotated rapidly enough, we ask if this angular momentum is sufficient to form vortices in BEC haloes, affecting their structure with potentially observable consequences. The minimum angular momentum for this, L_{QM} = $\hbar M/m$, requires m >= 9.5 m_H for lambda = 0.05, close to the particle mass required to influence structure on galactic scales. We study the equilibrium of self-gravitating, rotating BEC haloes which satisfy the Gross-Pitaevskii-Poisson equations, to calculate if and when vortices are energetically favoured. Vortices form as long as self-interaction is strong enough, which includes a large part of the range of m and g of interest for BEC-CDM haloes.

If cold dark matter elementary particles form a Bose-Einstein condensate, their superfluidity may distinguish them from other forms of cold dark matter, including creation of quantum vortices. We demonstrate here that such vortices are favoured in strongly-coupled condensates, while this is not the case for axions, which are generally presumed to form a Bose-Einstein condensate but are effectively non-interacting.

In the past two decades, it has been established by high-resolution observations of early-type galaxies that their nuclear surface brightness and corresponding stellar mass densities are characterized by cusps. In this paper, we present a new spherical analytical model family describing mild cuspy centres. We study isotropic and anisotropic models of Osipkov-Merritt type. It is shown that the associated distribution functions and intrinsic velocity dispersions can be represented analytically in a unified way in terms of hypergeometric series, allowing thus a straightforward comparison of these important global quantities for galaxies having underlying mass densities which may differ significantly in their degree of central cuspiness or radial falloff.