Search

We review dark energy models which can present non-negligible fluctuations on scales smaller than Hubble radius. Both linear and nonlinear evolutions of dark energy fluctuations are discussed. The linear evolution has a well-established framework, based on linear perturbation theory in General Relativity, and is well studied and implemented in numerical codes. We highlight the main results from linear theory to explain how dark energy perturbations become important on the scales of interest for structure formation. Next, we review some attempts to understand the impact of clustering dark energy models in the nonlinear regime, usually based on generalizations of the Spherical Collapse Model. We critically discuss the proposed generalizations of the Spherical Collapse Model that can treat clustering dark energy models and their shortcomings. Proposed implementations of clustering dark energy models in halo mass functions are reviewed. We also discuss some recent numerical simulations capable of treating dark energy fluctuations. Finally, we make an overview of the observational predictions based on these models., Comment: Accepted version of a review published in Universe journal, 34 pages, 7 figures

We study an expanding two-fluid model of non-relativistic dark matter and radiation which are allowed to interact during a certain time span and to establish an approximate thermal equilibrium. Such interaction which generates an effective bulk viscous pressure at background level is expected to be relevant for times around the transition from radiation to matter dominance. We quantify the magnitude of this pressure for dark matter particles masses within the range $1 {\rm eV} \lesssim m_{\chi} \lesssim 10 {\rm eV}$ around the matter-radiation equality epoch (i.e., redshift $z_{\rm eq}\sim 3400$) and demonstrate that the existence of a transient bulk viscosity has consequences which may be relevant for addressing current tensions of the standard cosmological model: i) the additional (negative) pressure contribution modifies the expansion rate around $z_{\rm eq}$, yielding a larger $H_0$ value and ii) large-scale structure formation is impacted by suppressing the amplitude of matter overdensity growth via a new viscous friction-term contribution to the M\'esz\'aros effect. As a result, both the $H_0$ and the $S_8$ tensions of the current standard cosmological model are significantly alleviated., Comment: 8 pages/5 figures. V.2: version accepted for publication in PRD

A range of cosmological observations demonstrate an accelerated expansion of the Universe, and the most likely explanation of this phenomenon is a cosmological constant. Given the importance of understanding the underlying physics, it is relevant to investigate alternative models. This article uses numerical simulations to test the consistency of one such alternative model. Specifically, this model has no cosmological constant, instead the dark matter particles have an extra force proportional to velocity squared, somewhat reminiscent of the magnetic force in electrodynamics. The constant strength of the force is the only free parameter. Since bottom-up structure formation creates cosmological structures whose internal velocity dispersions increase in time, this model may mimic the temporal evolution of the effect from a cosmological constant. It is shown that models with force linearly proportional to internal velocites, or models proportional to velocity to power three or more cannot mimic the accelerated expansion induced by a cosmological constant. However, models proportional to velocity squared are still consistent with the temporal evolution of a Universe with a cosmological model., Comment: 5 pages, 3 figures, ApJ format

Both simulation and observational data have shown that the spin and shape of dark matter halos are correlated with their nearby large-scale environment. As structure formation on different scales is strongly coupled, it is trick to disentangle the formation of halo with the large-scale environment, making it difficult to infer which is the driving force for the correlation between halo spin/shape with the large-scale structure. In this paper, we use N-body simulation to produce twin Universes that share the same initial conditions on small scales but different on large scales. This is achieved by changing the random seeds for the phase of those k modes smaller than a given scale in the initial conditions. In this way, we are able to disentangle the formation of halo and large-scale structure, making it possible to investigate how halo spin and shape correspond to the change of environment on large scales. We identify matching halo pairs in the twin simulations as those sharing the maximum number of identical particles within each other. Using these matched halo pairs, we study the cross match of halo spin and their correlation with the large-scale structure. It is found that when the large-scale environment changes (eigenvector) between the twin simulations, the halo spin has to rotate accordingly, although not significantly, to maintain the universal correlation seen in each simulation. Our results suggest that the large-scale structure is the main factor to drive the correlation between halo properties and their environment., Comment: 12 pages, 8 figures, accepted for publication in ApJ

An ultra-light bosonic particle of mass around $10^{-22}\,\mathrm{eV}/c^2$ is of special interest as a dark matter candidate, as it both has particle physics motivations, and may give rise to notable differences in the structures on highly non-linear scales due to the manifestation of quantum-physical wave effects on macroscopic scales, which could address a number of contentious small-scale tensions in the standard cosmological model, $\Lambda$CDM. Using a spectral technique, we here discuss simulations of such fuzzy dark matter (FDM), including the full non-linear wave dynamics, with a comparatively large dynamic range and for larger box sizes than considered previously. While the impact of suppressed small-scale power in the initial conditions associated with FDM has been studied before, the characteristic FDM dynamics are often neglected; in our simulations, we instead show the impact of the full non-linear dynamics on physical observables. We focus on the evolution of the matter power spectrum, give first results for the FDM halo mass function directly based on full FDM simulations, and discuss the computational challenges associated with the FDM equations. FDM shows a pronounced suppression of power on small scales relative to cold dark matter (CDM), which can be understood as a damping effect due to 'quantum pressure'. In certain regimes, however, the FDM power can exceed that of CDM, which may be interpreted as a reflection of order-unity density fluctuations occurring in FDM. In the halo mass function, FDM shows a significant abundance reduction below a characteristic mass scale only. This could in principle alleviate the need to invoke very strong feedback processes in small galaxies to reconcile $\Lambda$CDM with the observed galaxy luminosity function, but detailed studies that also include baryons will be needed to ultimately judge the viability of FDM., Comment: 16 pages, 14 figures; revised version submitted to MNRAS

Sterile neutrinos ($\nu_s$) that mix with active neutrinos ($\nu_a$) are interesting dark matter candidates with a rich cosmological and astrophysical phenomenology. In their simplest incarnation, their production is severely constrained by a combination of structure formation observations and X-ray searches. We show that if active neutrinos couple to an oscillating condensate of a very light $L_{\mu}-L_{\tau}$ gauge field, resonant $\nu_a$-$\nu_s$ oscillations can occur in the early universe, consistent with $\nu_s$ constituting all of the dark matter, while respecting X-ray constraints on $\nu_s\to\nu_a\gamma$ decays. Interesting deviations from standard solar and atmospheric neutrino oscillations can persist to the present., Comment: 15 pages, 8 figures. v2: some clarifications and references added, small errors corrected, main conclusions unchanged

In the next decade, cosmological surveys will have the statistical power to detect the absolute neutrino mass scale. N-body simulations of large-scale structure formation play a central role in interpreting data from such surveys. Yet these simulations are Newtonian in nature. We provide a quantitative study of the limitations to treating neutrinos, implemented as N-body particles, in N-body codes, focusing on the error introduced by neglecting special relativistic effects. Special relativistic effects are potentially important due to the large thermal velocities of neutrino particles in the simulation box. We derive a self-consistent theory of linear perturbations in Newtonian and non-relativistic neutrinos and use this to demonstrate that N-body simulations overestimate the neutrino free-streaming scale, and cause errors in the matter power spectrum that depend on the initial redshift of the simulations. For $z_{i} \lesssim 100$, and neutrino masses within the currently allowed range, this error is $\lesssim 0.5\%$, though represents an up to $\sim 10\%$ correction to the shape of the neutrino-induced suppression to the cold dark matter power spectrum. We argue that the simulations accurately model non-linear clustering of neutrinos so that the error is confined to linear scales., Comment: 15 pages, 7 figures. Matching prd accepted version

The dark matter (DM) can consist of the primordial black holes (PBHs) in addition to the conventional weakly interacting massive particles (WIMPs). The Poisson fluctuations of the PBH number density produce the isocurvature perturbations which can dominate the matter power spectrum at small scales and enhance the early structure formation. We study how the WIMP annihilation from those early formed structures can affect the CMB (in particular the E-mode polarization anisotropies and $y$-type spectral distortions) and global 21cm signals. Our studies would be of particular interest for the light (sub-GeV) WIMP scenarios which have been less explored compared with the mixed DM scenarios consisting of PBHs and heavy ($\gtrsim 1$ GeV) WIMPs. For instance, for the self-annihilating DM mass $m_{\chi}=1$ MeV and the thermally averaged annihilation cross section $\langle \sigma v \rangle \sim 10^{-30} \rm cm^3/s$, the latest Planck CMB data requires the PBH fraction with respect to the whole DM to be at most ${\cal O}(10^{-3})$ for the sub-solar mass PBHs and an even tighter bound (by a factor $\sim 5$) can be obtained from the global 21-cm measurements., Comment: 14 pages, 10 figures

Morphologies of adsorbed molecular films are of interest in a wide range of applications. To study the epitaxial growth of these systems in computer simulations requires access to long time and length scales, and one typically resorts to kinetic Monte Carlo (KMC) simulations. However, KMC simulations require as input transition rates and their dependence on external parameters (such as temperature). Experimental data allow only limited and indirect access to these rates, and models are often oversimplified. Here, we follow a bottom-up approach and aim at systematically constructing all relevant rates for an example system that has shown interesting properties in experiments, buckminsterfullerene on a calcium fluoride substrate. We develop classical force fields (both atomistic and coarse-grained) and perform molecular dynamics simulations of the elementary transitions in order to derive explicit expressions for the transition rates with a minimal number of free parameters.

In the hierarchical scenario of structure formation, galaxy clusters are the ultimate virialised products in mass and time. Hot baryons in the intracluster medium (ICM) and cold baryons in galaxies inhabit a dark matter dominated halo. Internal processes, accretion, and mergers can perturb the equilibrium, which is established only at later times. However, the cosmic time when thermalisation is effective is still to be assessed. Here we show that massive clusters in the observed universe attained an advanced thermal equilibrium $\sim~1.8~\text{Gyr}$ ago, at redshift $z =0.14\pm0.06$, when the universe was $11.7\pm0.7~\text{Gyr}$ old. Hot gas is mostly thermalised after the time when cosmic densities of matter and dark energy match. We find in a statistically nearly complete and homogeneous sample of 120 clusters from the {\it Planck} Early Sunyaev-Zel'dovich (ESZ) sample that the kinetic energy traced by the galaxy velocity dispersion is a faithful probe of the gravitational energy since a look back time of at least $\sim5.4~\text{Gyr}$, whereas the efficiency of hot gas in converting kinetic to thermal energy, as measured through X-ray observations in the core-excised area within $r_{500}$, steadily increases with time. The evolution is detected at the $\sim 98$ per cent probability level. Our results demonstrate that halo mass accretion history plays a larger role for cluster thermal equilibrium than radiative physics. The evolution of hot gas is strictly connected to the cosmic structure formation., Comment: 10 pages, 3 figures, in press on MNRAS

We derive structure formation limits on dark matter (DM) composed of keV-scale axion-like particles (ALPs), produced via freeze-in through the interactions with photons and Standard Model (SM) fermions. We employ Lyman-alpha (Ly-{\alpha}) forest data sets as well as the observed number of Milky Way (MW) subhalos. We compare results obtained using Maxwell-Boltzmann and quantum statistics for describing the SM bath. It should be emphasized that the presence of logarithmic divergences complicates the calculation of the production rate, which can not be parameterized with a simple power law behaviour. The obtained results, in combination with X-ray bounds, exclude the possibility for a photophilic "frozen-in" ALP DM with mass below $\sim 19\,\mathrm{keV}$. For the photophobic ALP scenario, in which DM couples primarily to SM fermions, the ALP DM distribution function is peaked at somewhat lower momentum and hence for such realization we find weaker limits on DM mass. Future facilities, such as the upcoming Vera C. Rubin observatory, will provide measurements with which the current bounds can be significantly improved to $\sim 80\,\mathrm{keV}$., Comment: 18 pages, 5 figures, minor updates in the discussion, main results unaltered

Quasi-N-body simulations, such as FastPM, provide a fast way to simulate cosmological structure formation, but have yet to adequately include the effects of massive neutrinos. We present a method to include neutrino particles in FastPM, enabling computation of the CDM and total matter power spectra to percent-level accuracy in the non-linear regime. The CDM-neutrino cross-power can also be computed at a sufficient accuracy to constrain cosmological observables. To avoid the shot noise that typically plagues neutrino particle simulations, we employ a quasi-random algorithm to sample the relevant Fermi-Dirac distribution when setting the initial neutrino thermal velocities. We additionally develop an effective distribution function to describe a set of non-degenerate neutrinos as a single particle to speed up non-degenerate simulations. The simulation is accurate for the full range of physical interest, $M_\nu \lesssim 0.6$eV, and applicable to redshifts $z\lesssim2$. Such accuracy can be achieved by initializing particles with the two-fluid approximation transfer functions (using the REPS package). Convergence can be reached in $\sim 25$ steps, with a starting redshift of $z=99$. Probing progressively smaller scales only requires an increase in the number of CDM particles being simulated, while the number of neutrino particles can remain fixed at a value less than or similar to the number of CDM particles. In turn, the percentage increase in runtime-per-step due to neutrino particles is between $\sim 5-20\%$ for runs with $1024^3$ CDM particles, and decreases as the number of CDM particles is increased. The code has been made publicly available, providing an invaluable resource to produce fast predictions for cosmological surveys and studying reconstruction., Comment: 20 pages, 8 figures, 1 table. Code available at https://github.com/fastpm/fastpm

We study the highly magnified arc SGAS J122651.3+215220 caused by a star-forming galaxy at $z_s=2.93$ crossing the lensing caustic cast by the galaxy cluster SDSS J1226+2152 ($z_l=0.43$), using Hubble Space Telescope observations. We report in the arc several asymmetric surface brightness features whose angular separations are a fraction of an arcsecond from the lensing critical curve and appear to be highly but unequally magnified image pairs of underlying compact sources, with one brightest pair having clear asymmetry consistently across four filters. One explanation of unequal magnification is microlensing by intracluster stars, which induces independent flux variations in the images of individual or groups of source stars in the lensed galaxy. For a second possibility, intracluster dark matter subhalos invisible to telescopes effectively perturb lensing magnifications near the critical curve and give rise to persistently unequal image pairs. Our modeling suggests, at least for the most prominent identified image pair, that the microlensing hypothesis is in tension with the absence of notable asymmetry variation over a six-year baseline, while subhalos of $\sim 10^6$--$10^8\,M_\odot$ anticipated from structure formation with Cold Dark Matter typically produce stationary and sizable asymmetries. We judge that observations at additional times and more precise lens models are necessary to stringently constrain temporal variability and robustly distinguish between the two explanations. The arc under this study is a scheduled target of a Director's Discretionary Early Release Science program of the James Webb Space Telescope, which will provide deep images and a high-resolution view with integral field spectroscopy., Comment: New version accepted by MNRAS; 18 pages including references and appendices, 13 figures and 4 tables; major revision of Sec. 3.2 and Figure 4 presenting improved data analysis; original conclusion strengthened

A strongly interacting massive particle (SIMP) is an interesting candidate for dark matter (DM) because its self-interaction cross section can be naturally strong enough to address the astrophysical problem of small-scale structure formation. A simple model was proposed by assuming a monopole condensation, where composite SIMP comes from a "strongly interacting" U(1)$_{\rm d}$ gauge theory. In the original model, the DM relic abundance is determined by the $3\to2$ annihilation process via the Wess-Zumino-Witten term. In this letter, we discuss that the DM relic abundance is naturally determined also by a semi-annihilation process via a kinetic mixing between the hypercharge gauge boson and the dark U(1)$_{\rm d}$ gauge boson (dark photon). The dark photon can be discovered by LDMX-style missing momentum experiments in the near future., Comment: 8 pages, 2 figures; v2: figure updated, minor comments added

The Epoch of Reionization (EoR) depends on the complex astrophysics governing the birth and evolution of the first galaxies and structures in the intergalactic medium. EoR models rely on cosmic microwave background (CMB) observations, and in particular the large-scale E-mode polarization power spectra (EE PS), to help constrain their highly uncertain parameters. However, rather than directly forward-modelling the EE PS, most EoR models are constrained using a summary statistic -- the Thompson scattering optical depth, $\tau_e$. Compressing CMB observations to $\tau_e$ requires adopting a basis set for the EoR history. The common choice is the unphysical, redshift-symmetric hyperbolic tangent (Tanh) function, which differs in shape from physical EoR models based on hierarchical structure formation. Combining public EoR and CMB codes, 21cmFAST and CLASS, here we quantify how inference using the $\tau_e$ summary statistic impacts the resulting constraints on galaxy properties and EoR histories. Using the last Planck 2018 data release, we show that the marginalized constraints on the EoR history are more sensitive to the choice of the basis set (Tanh vs physical model) than to the CMB likelihood statistic ($\tau_e$ vs PS). For example, EoR histories implied by the growth of structure show a small tail of partial reionization extending to higher redshifts. However, biases in inference using $\tau_e$ are negligible for the Planck 2018 data. Using EoR constraints from high-redshift observations including the quasar dark fraction, galaxy UV luminosity functions and CMB EE PS, our physical model recovers $\tau_e=0.0569^{+0.0081}_{-0.0066}$., Comment: 11 pages, 5 figures, Updated to match the published version. Added Gamma evolution; Only minor changes

The early Universe may have passed through an extended period of matter-dominated expansion following inflation and prior to the onset of radiation domination. Sub-horizon density perturbations grow gravitationally during such an epoch, collapsing into bound structures if it lasts long enough. The strong analogy between this phase and structure formation in the present-day universe allows the use of N-body simulations and approximate methods for halo formation to model the fragmentation of the inflaton condensate into inflaton halos. For a simple model we find these halos have masses of up to $20\,\mathrm{kg}$ and radii of the order of $10^{-20}\,\mathrm{m}$, roughly $10^{-24}$ seconds after the Big Bang. We find that the N-body halo mass function matches predictions of the mass-Peak Patch method and the Press-Schechter formalism within the expected range of scales. A long matter-dominated phase would imply that reheating and thermalization occurs in a universe with large variations in density, potentially modifying the dynamics of this process. In addition, large overdensities can source gravitational waves and may lead to the formation of primordial black holes., Comment: 12 pages, 8 figures, accepted by PRD

Using $N$-body simulations of cosmological large-scale structure formation, for the first time, we show that the anisotropic primordial non-Gaussianity (PNG) causes a scale-dependent modification, given by $1/k^2$ at small $k$ limit, in the three-dimensional power spectra of halo shapes (intrinsic alignments), whilst the conventional power spectrum of halo number density field remains unaffected. We discuss that wide-area imaging and spectrocopic surveys observing the same region of the sky allow us to constrain the quadrupole PNG coefficient $f_{\rm NL}^{s=2}$ at a precision comparable with or better than that of the cosmic microwave background., Comment: 9 pages, 4 figures; an accepted version for publication in PRD

We study the future prospects of the 21cm forest observations on the axion-like dark matter when the spontaneous breaking of the global Peccei-Quinn (PQ) symmetry occurs after the inflation. The large isocurvature perturbations of order unity sourced from axion-like particles can result in the enhancement of minihalo formation, and the subsequent hierarchical structure formation can affect the minihalo abundance whose masses can exceed ${\cal O}(10^4) M_{\odot}$ relevant for the 21cm forest observations. We show that the 21cm forest observations are capable of probing the axion-like particle mass in the range $10^{-18}\lesssim m_a \lesssim 10^{-12}$ eV for the temperature independent axion mass. For the temperature dependent axion mass, the zero temperature axion mass scale for which the 21cm forest measurements can be affected is extended further to as big as of order $10^{-6}$ eV., Comment: 11 pages, 6 figures, published in PRD

We study structure formation in a set of cosmological simulations to uncover the scales in the initial density field that gave rise to the formation of present-day structures. Our simulations share a common primordial power spectrum (here Lambda-CDM), but the introduction of hierarchical variations of the phase information allows us to systematically study the scales that determine the formation of structure at later times. We consider the variance in z=0 statistics such as the matter power spectrum and halo mass function. We also define a criterion for the existence of individual haloes across simulations, and determine what scales in the initial density field contain sufficient information for the non-linear formation of unique haloes. We study how the characteristics of individual haloes such as the mass and concentration, as well as the position and velocity, are affected by variations on different scales, and give scaling relations for haloes of different mass. Finally, we use the example of a cluster-mass halo to show how our hierarchical parametrisation of the initial density field can be used to create variants of particular objects. With properties such as mass, concentration, kinematics and substructure of haloes set on distinct and well-determined scales, and its unique ability to introduce variations localised in real space, our method is a powerful tool to study structure formation in cosmological simulations., Comment: 18 pages, 22 figures, submitted to MNRAS

The nature of dark matter sets the timeline for the formation of first collapsed halos and thus affects the sources of reionization. Here, we consider two different models of dark matter: cold dark matter (CDM) and thermal warm dark matter (WDM), and study how they impact the epoch of reionization (EoR) and its 21-cm observables. Using a suite of simulations, we find that in the WDM scenarios, the structure formation on small scales gets suppressed resulting in a smaller number of low mass dark matter halos compared to the CDM scenario. Assuming that the efficiency of sources in producing ionizing photons remain the same, this leads to a lower number of total ionizing photons produced at any given cosmic time and thus in a delay in the reionization process. We also find visual differences in the neutral hydrogen (HI) topology and in 21-cm maps in case of the WDM compared to the CDM. However, differences in the 21-cm power spectra, at the same neutral fraction, are found to be small. Thus, we focus on the non-Gaussianity in the EoR 21-cm signal, quantified through its bispectrum. We find that the 21-cm bispectra (driven by the HI topology) are significantly different in WDM models when compared with CDM, even for same mass averaged neutral fractions. This establishes that the 21-cm bispectrum is a unique and promising way to differentiate between different dark matter models, and can be used to constrain the nature of the dark matter in the future EoR observations., Comment: 13 pages, 13 figures, accepted for publication in MNRAS

We study the statistical mechanics of binary systems under gravitational interaction of the Modified Newtonian Dynamics (MOND) in three-dimensional space. Considering the binary systems, in the microcanonical and canonical ensembles, we show that in the microcanonical systems, unlike the Newtonian gravity, there is a sharp phase transition, with a high-temperature homogeneous phase and a low temperature clumped binary one. Defining an order parameter in the canonical systems, we find a smoother phase transition and identify the corresponding critical temperature in terms of physical parameters of the binary system., Comment: 14 pages, 5 figures

The forthcoming Laser Interferometer Space Antenna (LISA) will probe the population of coalescing massive black hole (MBH) binaries up to the onset of structure formation. Here we simulate the galactic-scale pairing of $\sim10^6 M_\odot$ MBHs in a typical, non-clumpy main-sequence galaxy embedded in a cosmological environment at $z = 7-6$. In order to increase our statistical sample, we adopt a strategy that allows us to follow the evolution of six secondary MBHs concomitantly. We find that the magnitude of the dynamical-friction induced torques is significantly smaller than that of the large-scale, stochastic gravitational torques arising from the perturbed and morphologically evolving galactic disc, suggesting that the standard dynamical friction treatment is inadequate for realistic galaxies at high redshift. The dynamical evolution of MBHs is very stochastic, and a variation in the initial orbital phase can lead to a drastically different time-scale for the inspiral. Most remarkably, the development of a galactic bar in the host system either significantly accelerates the inspiral by dragging a secondary MBH into the centre, or ultimately hinders the orbital decay by scattering the MBH in the galaxy outskirts. The latter occurs more rarely, suggesting that galactic bars overall promote MBH inspiral and binary coalescence. The orbital decay time can be an order of magnitude shorter than what would be predicted relying on dynamical friction alone. The stochasticity, and the important role of global torques, have crucial implications for the rates of MBH coalescences in the early Universe: both have to be accounted for when making predictions for the upcoming LISA observatory., Comment: Accepted for publication in MNRAS; 15 pages, 10 Figures, 2 Tables

If the symmetry breaking inducing the axion occurs after the inflation, the large axion isocurvature perturbations can arise due to a different axion amplitude in each causally disconnected patch. This causes the enhancement of the small-scale density fluctuations which can significantly affect the evolution of structure formation. The epoch of the small halo formation becomes earlier and we estimate the abundance of those minihalos which can host the neutral hydrogen atoms to result in the 21cm fluctuation signals. We find that the future radio telescopes, such as the SKA, can put the axion mass bound of order $m_a \gtrsim 10^{-13}$ eV for the simple temperature-independent axion mass model, and the bound can be extended to of order $m_a \gtrsim 10^{-8}$eV for a temperature-dependent axion mass., Comment: 15 pages, 10 figures

The recent discovery of diffuse dwarf galaxies that are deficient in dark matter appears to challenge the current paradigm of structure formation in our Universe. We describe the numerical experiments to determine if the so-called dark matter deficient galaxies (DMDGs) could be produced when two gas-rich, dwarf-sized galaxies collide with a high relative velocity of $\sim 300\,{\rm kms^{-1}}$. Using idealized high-resolution simulations with both mesh-based and particle-based gravito-hydrodynamics codes, we find that DMDGs can form as high-velocity galaxy collisions separate dark matter from the warm disk gas which subsequently is compressed by shock and tidal interaction to form stars. Then using a large simulated universe IllustrisTNG, we discover a number of high-velocity galaxy collision events in which DMDGs are expected to form. However, we did not find evidence that these types of collisions actually produced DMDGs in the TNG100-1 run. We argue that the resolution of the numerical experiment is critical to realize the "collision-induced" DMDG formation scenario. Our results demonstrate one of many routes in which galaxies could form with unconventional dark matter fractions., Comment: 18 pages including 13 figures. Accepted for publication in ApJ. For associated movies and information, see https://ejshin0.github.io/page/category/visualizations.html and https://physics.snu.ac.kr/cosmo/index/research.html#dmdg

Perturbations from inflation evolve into large scale structure of the late universe, and encode abundant cosmic structure formation physics. We allow freedom in the primordial power spectrum, rather than assuming a power law scale dependence, to study its impact on cosmological parameter determination. Combining various generations of cosmic microwave background (CMB) data and galaxy redshift survey data, we investigate the constraints on reconstruction of the primordial curvature perturbation power spectrum and the late time cosmology, especially the sum of neutrino masses. We quantify how each successive generation, in CMB and galaxy surveys, provides significant improvements, often by factors of several. By using CMB polarization information over a broad range of angular scales, and galaxy redshift data in many bins of redshift, one can allow inflationary freedom and still constrain parameters comparably to assuming power law dependence. The primordial power spectrum can be reconstructed at the subpercent level in a dozen wavenumber bins, while simultaneously fitting the sum of neutrino masses to 14 meV., Comment: 19 pages, 12 figures, 2 tables, v2 accepted for publication

The mass function of dark matter halos is one of the most fundamental statistics in structure formation. Many theoretical models (such as Press-Schechter theory) are based on the notion that it could be universal, meaning independent of redshift and cosmology, when expressed in the appropriate variables. However, simulations exhibit persistent non-universalities in the mass functions of the virial mass and other commonly used spherical overdensity definitions. We systematically study the universality of mass functions over a wide range of mass definitions, for the first time including the recently proposed splashback mass, Msp. We confirm that, in LambdaCDM cosmologies, all mass definitions exhibit varying levels of non-universality that increase with peak height and reach between 20% and 500% at the highest masses we can test. Mvir, M200m, and Msp exhibit similar levels of non-universality. There are, however, two regimes where the splashback mass functions are significantly more universal. First, they are universal to 10% at z, 15 pages, 7 figures. Small changes to match the published version

Predicting the merger timescale ($\tau_{\rm merge}$) of merging dark matter halos, based on their orbital parameters and the structural properties of their hosts, is a fundamental problem in gravitational dynamics that has important consequences for our understanding of cosmological structure formation and galaxy formation. Previous models predicting $\tau_{\rm merge}$ have shown varying degrees of success when compared to the results of cosmological $N$-body simulations. We build on this previous work and propose a new model for $\tau_{\rm merge}$ that draws on insights derived from these simulations. We find that published predictions can provide reasonable estimates for $\tau_{\rm merge}$ based on orbital properties at infall, but tend to underpredict $\tau_{\rm merge}$ inside the host virial radius ($R_{200}$) because tidal stripping is neglected, and overpredict it outside $R_{200}$ because the host mass is underestimated. Furthermore, we find that models that account for orbital angular momentum via the circular radius $R_{\rm circ}$ underpredict (overpredict) $\tau_{\rm merge}$ for bound (unbound) systems. By fitting for the dependence of $\tau_{\rm merge}$ on various orbital and host halo properties,we derive an improved model for $\tau_{\rm merge}$ that can be applied to a merging halo at any point in its orbit. Finally, we discuss briefly the implications of our new model for $\tau_{\rm merge}$ for semi-analytical galaxy formation modelling., Comment: 12 pages, 11 figures, accepted for publication in MNRAS

The simplest model of mirror sector dark matter maintains exact mirror symmetry, but has a baryon abundance $\Omega_{b'} = \beta \Omega_b$ and a suppressed temperature $T' = x T$ in the mirror sector; hence it depends only on two parameters, $\beta,x$. For sufficiently small $x$, early cosmological observables may not constrain mirror baryons from constituting all of the dark matter despite their strong self-interactions, depending on the unknown details of structure formation in the hidden sector. Here we close this loophole by simulating mirror structure formation, mapping out the allowed regions of parameter space using cosmological and astronomical data. We find that the Milky Way disk surface density and bulge mass constrain $\Omega_{b'}\lesssim 0.3 \Omega_{b}$ at the highest $T'$ allowed by BBN and CMB ($T'=0.5 T$), or $\Omega_{b'}\lesssim 0.8 \Omega_{b}$ at lower values of $T'$. We also briefly discuss the realization of the necessary temperature asymmetry between the SM and the mirror sector in our model with unbroken mirror symmetry., Comment: 30 pages, 8 figures

Cold dark matter (DM) models for structure formation predict that DM subhalos are present in the Galaxy. In the standard paradigm of DM as weakly interacting massive particle, subhalos are expected to shine in gamma rays and to provide a signal detectable with current instruments, notably with the Large Area Telescope (LAT) aboard the Fermi~satellite. This is the main motivation behind searches for DM signals towards dwarf spheroidal galaxies and unidentified Fermi-LAT sources. A significant angular extension detected from unassociated sources located at relatively high latitudes is considered a "smoking gun" signature for identifying DM subhalos. In the present work, we systematically explore, by means of state-of-the-art models of cold DM halos in the Galaxy, the detectability of extended subhalos with Fermi-LAT. We simulate a DM signal exploring different assumptions of subhalos distribution in the Galaxy and DM profile, and reconstruct its flux through a realistic Fermi-LAT analysis pipeline. In the most optimistic case, we show that a detection of extended DM subhalos can be made for annihilation cross sections higher than $3 \times 10^{-26}$ cm$^3$/s (for a 100 GeV DM mass), still compatible with existing gamma-ray constraints, and that, in this case, the preference for extension of the source (vs point-like hypothesis) is significant. For fainter signals, instead, halos not only do not show significant extension, but they are not even detectable significantly as point-like sources., Comment: 12 pages, 11 figures. Comments as welcome!

To generate initial conditions for cosmological $N$-body simulations, one needs to prepare a uniform distribution of simulation particles, so-called the pre-initial condition (pre-IC). The standard method to construct the pre-IC is to place the particles on the lattice grids evenly spaced in the three-dimensional spatial coordinates. However, even after the initial displacement of each particle according to cosmological perturbations, the particle distribution remains to display an artificial anisotropy. Such an artifact causes systematic effects in simulations at later time until the evolved particle distribution sufficiently erases the initial anisotropy. In this paper, we study the impacts of the pre-IC on the anisotropic separate universe simulation, where the effect of large-scale tidal field on structure formation is taken into account using the anisotropic expansion in a local background (simulation volume). To quantify the impacts, we compare the simulations employing the standard grid pre-IC and the glass one, where the latter is supposed to suppress the initial anisotropy. We show that the artificial features in the grid pre-IC simulations are seen until $z\sim 9$, while the glass pre-IC simulations appear to be stable and accurate over the range of scales we study. From these results we find that a coupling of the large-scale tidal field with matter clustering is enhanced compared to the leading-order prediction of perturbation theory in the quasi non-linear regime in the redshift range $5\lesssim z\lesssim 15$, indicating the importance of tidal field on structure formation at such high redshifts, e.g. during the epoch of reionization., Comment: 11 pages, 8 figures, 1 table, accepted for publication in MNRAS, added a subsection and Appendix

Axion cold dark matter from standard misalignment typically requires a decay constant $f_a~\gtrsim~10^{11}$ GeV. Kinetic misalignment and parametric resonance easily allow lower values of $f_a$ when the radial Peccei-Quinn (PQ) symmetry breaking field takes large initial values. Here, we consider the effects of inflation on kinetic misalignment and parametric resonance. We assume that the initial PQ field value is determined by quantum fluctuations, and is set by the Hubble parameter during inflation, $H_I$, and the PQ field mass. PQ field oscillations begin before or after the completion of reheating after inflation at a temperature $T_R$. We determine the range of $f_a$ and the inflationary parameters $(H_I, T_R)$ consistent with axion dark matter for a quartic potential for the PQ field. We find that $4\times 10^8$ GeV $< f_a < 10^{11}$ GeV can consistently produce axion dark matter. A significant portion of the allowed parameter space predicts rare kaon decays, $K_L \rightarrow (��^0 + \rm{missing \; energy})$, and/or suppression of structure formation on small scales., 40 pages, 6 figures

We study how statistical properties of supermassive black holes depend on the frequency and conditions for massive seed formation in cosmological simulations of structure formation. We develop a novel method to recalculate detailed growth histories and merger trees of black holes within the framework of the Illustris simulation for several seed formation models, including a physically motivated model where black hole seeds only form in progenitor galaxies that conform to the conditions for direct collapse black hole formation. While all seed models considered here are in a broad agreement with present observational constraints on black hole populations from optical, UV and X-ray studies, we find that they lead to widely different black hole number densities and halo occupation fractions which are currently observationally unconstrained. In terms of future electromagnetic spectrum observations, the faint-end quasar luminosity function and the low-mass-end black hole-host galaxy scaling relations are very sensitive to the specific massive seed prescription. Specifically, the direct collapse model exhibits a seeding efficiency which decreases rapidly with cosmic time and produces much fewer black holes in low mass galaxies, in contrast to the original Illustris simulation. We further find that the total black hole merger rate varies by more than one order of magnitude for different seed models, with the redshift evolution of the chirp mass changing as well. Supermassive black hole merger detections with LISA and International Pulsar Timing Array may hence provide the most direct means of constraining massive black hole seed formation in the early Universe., Comment: 20 pages, 17 figures; Published in MNRAS

The first paper in this series showed that QCD axion dark matter, as a highly-correlated Bose fluid, contains extra-classical physics on cosmological scales. The nature of the derived extra-classical physics is exchange-correlation interactions induced by the constraints of symmetric particle exchange and inter-axion correlations from self-gravitation. The paper also showed that the impact of extra-classical physics on early structure formation is marginal, as the exchange-correlation interaction is inherently non-linear. This paper continues the study of axion structure formation into the non-linear regime by now considering the case of full collapse and virialization. The N-body method is chosen to solve the collapse problem, and its algorithms are derived for a condensed Bose fluid. Simulations of isolated gravitational collapse are performed for both Bose and cold dark matter fluids using a prototype N-body code. Unique Bose structures are found to survive even the most violent collapses. Bose post-collapse features include dynamical changes to global structures, creation of new broad substructures, violations of classical binding energy conditions, and new fine structures. Effective models of the novel structures are constructed and possibilities for their discovery are discussed., Comment: 29 pages, 27 figures

The mass-concentration relation of dark matter halos reflects the assembly history of objects in hierarchical structure formation scenarios, and depends on fundamental quantities in cosmology such as the slope of the primordial matter power-spectrum. This relation is unconstrained by observations on sub-galactic scales. We derive the first measurement of the mass-concentration relation using the image positions and flux ratios from eleven quadruple-image strong gravitational lenses (quads) in the mass range $10^{6} - 10^{10} M_{\odot}$, assuming cold dark matter. Our analysis framework includes both subhalos and line of sight halos, marginalizes over nuisance parameters describing the lens macromodel, accounts for finite source effects on lensing observables, and simultaneously constrains the normalization and logarithmic slope of the mass-concentration relation, and the normalization of the subhalo mass function. At $z=0$, we constrain the concentration of $10^{8} M_{\odot}$ halos $c=12_{-5}^{+6}$ at $68 \%$ CI, and $c=12_{-9}^{+15}$ at $95 \%$ CI. For a $10^{7} M_{\odot}$ halo, we obtain $68 \%$ ($95 \%$) constraints $c=15_{-8}^{+9}$ ($c=15_{-11}^{+18}$), while for $10^{9} M_{\odot}$ halos $c=10_{-4}^{+7}$ ($c=10_{-7}^{+14}$). These results are consistent with the theoretical predictions from mass-concentration relations in the literature, and establish strong lensing by galaxies as a powerful probe of halo concentrations on sub-galactic scales across cosmological distance., Comment: fix typos, edit wording

Nowadays cosmological inflation is the most accepted mechanism to explain the primordial seeds that led to the structure formation observed in the Universe. Current observations are in well agreement to initial adiabatic conditions, which imply that single-scalar-field inflation may be enough to describe the early Universe. However, there are several scenarios where the existence of more than a single field could be relevant during this period, for instance, the situation where the so-called spectator is present. Within the spectator scenario we can find the possibility that an ultra-light scalar field dark matter candidate could coexist with the inflaton. In this work we study this possibility where the additional scalar field could be massive or self-interacting. We use isocurvature observations to constrain the free parameters of the model., Comment: 11 pages, 4 figures

Structure formation in turbulence is effectively an instability of "plasma" formed by fluctuations serving as particles. These "particles" are quantumlike; namely, their wavelengths are non-negligible compared to the sizes of background coherent structures. The corresponding "kinetic equation" describes the Wigner matrix of the turbulent field, and the coherent structures serve as collective fields. This formalism is usually applied to manifestly quantumlike or scalar waves. Here, we extend it to compressible Navier--Stokes turbulence, where the fluctuation Hamiltonian is a five-dimensional matrix operator and diverse modulational modes are present. As an example, we calculate these modes for a sinusoidal shear flow and find two modulational instabilities. One of them is specific to supersonic flows, and the other one is a Kelvin--Helmholtz-type instability that is a generalization of the known zonostrophic instability. This work serves as a stepping stone toward improving the understanding of magnetohydrodynamic turbulence, which can be approached similarly., Comment: 11 pages, 4 figures

We intend to understand cosmological structure formation within the framework of superfluid models of dark matter with finite temperatures. Of particular interest is the evolution of small-scale structures where the pressure and superfluid properties of the dark matter fluid are prominent. We compare the growth of structures in these models with the standard cold dark matter paradigm and non-superfluid dark matter. The equations for superfluid hydrodynamics were computed numerically in an expanding $\Lambda$CDM background with spherical symmetry; the effect of various superfluid fractions, temperatures, interactions, and masses on the collapse of structures was taken into consideration. We derived the linear perturbation of the superfluid equations, giving further insights into the dynamics of the superfluid collapse. We found that while a conventional dark matter fluid with self-interactions and finite temperatures experiences a suppression in the growth of structures on smaller scales, as expected due to the presence of pressure terms, a superfluid can collapse much more efficiently than was naively expected due to its ability to suppress the growth of entropy perturbations and thus gradients in the thermal pressure. We also found that the cores of the dark matter halos initially become more superfluid during the collapse, but eventually reach a point where the superfluid fraction falls sharply. The formation of superfluid dark matter halos surrounded by a normal fluid dark matter background is therefore disfavored by the present work., Comment: 10 pages, 9 figures. Accepted by A&A

A very large dynamic range with simultaneous capture of both large and small scales in the simulations of cosmic structures is required for correct modelling of many cosmological phenomena, particularly at high redshift. This is not always available, or when it is, it makes such simulations very expensive. We present a novel subgrid method for modelling low-mass ($10^5\, {\rm M}_\odot \le M_{\rm halo}\le 10^9\, {\rm M}_\odot$) haloes, which are otherwise unresolved in large-volume cosmological simulations limited in numerical resolution. In addition to the deterministic halo bias that captures the average property, we model its stochasticity that is correlated in time. We find that the instantaneous binned distribution of the number of haloes is well approximated by a lognormal distribution, with overall amplitude modulated by this ‘temporal correlation bias’. The robustness of our new scheme is tested against various statistical measures, and we find that temporally correlated stochasticity generates mock halo data that is significantly more reliable than that from temporally uncorrelated stochasticity. Our method can be applied for simulating processes that depend on both the small- and large-scale structures, especially for those that are sensitive to the evolution history of structure formation such as the process of cosmic reionization. As a sample application, we generate a mock distribution of medium-mass (108 ≤ M/M⊙ ≤ 109) haloes inside a 500 Mpc $\, h^{-1}$, 3003 grid simulation box. This mock halo catalogue bears a reasonable statistical agreement with a halo catalogue from numerically resolved haloes in a smaller box, and therefore will allow a very self-consistent sets of cosmic reionization simulations in a box large enough to generate statistically reliable data.

Yielding transition in isotropic soft materials under superposition of orthogonal deformation fields is known to follow von Mises criterion. However, in anisotropic soft materials von Mises criterion fails owing to preferred directions associated with the system. In this work we study a model anisotropic yield stress system: electrorheological (ER) fluids that show structure formation in the direction of electric field. We subject the ER fluids to superposition of orthogonal stress fields that leads to different yield stress values. We obtain a yielding state diagram by plotting normalized rotational shear stress against normalized radial shear stress corresponding to yield point for a given electric field. Remarkably, the state diagram validates the Hill yielding criterion, which is a general yielding criterion for materials having anisotropy along three orthogonal directions, originally developed for metallic systems. Validation of Hill criterion suggests the universality of its application to anisotropic systems including conventional anisotropic soft materials having yield stress., Comment: 6 pages, 7 figures

In a broad class of scenarios, inflation is followed by an extended era of matter-dominated expansion during which the inflaton condensate is nonrelativistic on subhorizon scales. During this phase density perturbations grow to the point of nonlinearity and collapse into bound structures. This epoch strongly resembles structure formation with ultra-light axion-like particles. This parallel permits us to adapt results from studies of cosmological structure formation to describe the nonlinear dynamics of this post-inflationary epoch. We show that the inflaton condensate fragments into "inflaton clusters", analogues of axion dark matter halos in present-day cosmology. Moreover, solitonic objects or "inflaton stars" can form inside these clusters, leading to density contrasts as large as $10^6$ in the post-inflationary universe., Comment: 11 pages; 3 figures; reformatted, references added, minor changes

Three general cases of dynamical interacting dark energy models ($\mathcal{D}$-class) are investigated in the context of Brans-Dicke cosmology. Some of important cosmological quantities are calculated for every cases as a function of redshift parameter. The most important part of this paper deals with fitting models with two different expansion history: (SNIa+BAO$_A$+$Omh^{2}$ and SNIa+BAO$_A$+H(z)) and with two different sets of data for Hubble parameter. This provides a remarkable feature to could analytically see the effects of each analyzes and each data sets on final results. The best fitted values of parameters according to these analyzes and data points, $\chi_{tot}^2/dof$, AIC and BIC are reported. By these diagnostic tools we found that some of these models have no chance against $\Lambda$CDM, even without need to study the structure formation, and could be ruled out. While some (e.g. $\mathcal{BD-D}C2$ and $\mathcal{BD-D}A^*$) render the best fit quality,i.e. the value of AIC and BIC and figures show that they fit perfectly with overall data and reveals a strong evidence in favor of these two models against $\Lambda$CDM., Comment: Accepted for publication in Astrophysics and Space Science journal. It contains 32 pages, 6 tables and 10 figures. arXiv admin note: substantial text overlap with arXiv:1609.03816

The rate at which the universe is expanding today is a fundamental parameter in cosmology which governs our understanding of structure formation and dark energy. However, current measurements of the Hubble constant, $H_0$, show a significant tension ($\sim 4-6\sigma$) between early- and late-Universe observations. There are ongoing efforts to check the diverse observational results and also to investigate possible theoretical ways to resolve the tension~-- which could point to radical extensions of the standard model. Here we demonstrate the potential of next-generation spectroscopic galaxy surveys to shed light on the Hubble constant tension. Surveys such as those with Euclid and the Square Kilometre Array (SKA) are expected to reach sub-percent precision on Baryon Acoustic Oscillation (BAO) measurements of the Hubble parameter, with a combined redshift coverage of $0.1, Comment: 6 pages, 4 figures, 1 table. JCAP accepted

We study the cosmological properties of a dynamical of dark energy (DE) component determined by a Steep Equation of State (SEoS) $w(z)=w_0+w_i\frac{(z/z_T)^q}{1+(z/z_T)^q}$. The SEoS has a transition at $z_T$ between two pivotal values ($w_i, w_0$) which can be taken as an early time and present day values of $w$ and the steepness is given by $q$. We describe the impact of this dynamical DE at background and perturbative level. The steepness of the transition has a better cosmological fit than a conventional CPL model with $w=w_0+w_a(1-a)$. Furthermore, we analyze the impact of steepness of the transition in the growth of matter perturbations and structure formation. This is manifest in the linear matter power spectrum, $P(k)$, the logarithmic growth function, $f\sigma_8(z)$, and the differential mass function $dn/d\log M(z=0)$. The differences in these last three quantities is at a percent-level using the same cosmological baseline parameters in our SEoS and a $\Lambda CDM$ model. However, we find an increase in the power spectrum, producing a bump at $k\approx k_T$ with $k_T\equiv a_TH(a_T)$ the mode associated to the time of the steep transition ($a_T = 1/(1+z_T)$). Different dynamics of DE lead to a different amount of DM at present time which has an impact in Power Spectrum and accordingly in structure formation., Comment: 10 pages, 7 figures

The majority of the ordinary matter in the local Universe has been heated by strong structure formation shocks and resides in a largely unexplored hot, diffuse, X-ray emitting plasma that permeates the halos of galaxies, galaxy groups and clusters, and the cosmic web. We propose a next-generation "Cosmic Web Explorer" that will permit a complete and exhaustive understanding of these unseen baryons. This will be the first mission capable to reach the accretion shocks located several times farther than the virial radii of galaxy clusters, and reveal the out-of-equilibrium parts of the intra-cluster medium which are live witnesses to the physics of cosmic accretion. It will also enable a view of the thermodynamics, kinematics, and chemical composition of the circumgalactic medium in galaxies with masses similar to the Milky Way, at the same level of detail that $Athena$ will unravel for the virialized regions of massive galaxy clusters, delivering a transformative understanding of the evolution of those galaxies in which most of the stars and metals in the Universe were formed. Finally, the proposed X-ray satellite will connect the dots of the large-scale structure by mapping, at high spectral resolution, as much as 100% of the diffuse gas hotter than $10^6$ K that fills the filaments of the cosmic web at low redshifts, down to an over-density of 1, both in emission and in absorption against the ubiquitous cosmic X-ray background, surveying at least 1600 square degrees over 5 years in orbit. This requires a large effective area (~10 m$^2$ at 1 keV) over a large field of view ($\sim1$ deg$^2$), a megapixel cryogenic microcalorimeter array providing integral field spectroscopy with a resolving power $E/\Delta E$ = 2000 at 0.6 keV and a spatial resolution of 5 arcsec in the soft X-ray band, and a low and stable instrumental background ensuring high sensitivity to faint, extended emission., Comment: White paper submitted in response to ESA's Voyage 2050 Call. Accepted for publication in Experimental Astronomy

We study the nonlinear regime of large scale structure formation considering a dynamical dark energy (DE) component determined by a Steep Equation of State parametrization (SEoS) $w(z)=w_0+w_i\frac{(z/z_T)^q}{1+(z/z_T)^q}$. In order to perform the model exploration at low computational cost, we modified the public code L-PICOLA. We incorporate the DE model by means of the first and second-order matter perturbations in the Lagrangian frame and the expansion parameter. We analyze deviations of SEoS models with respect to $\Lambda$CDM in the non-linear matter power spectrum ($P_k$), the halo mass function (HMF), and the two-point correlation function (2PCF). On quantifying the nature of steep (SEoS-I) and smooth transitions in DE field (CPL-lim), no signature of steep transition is observed, rather found the overall impact of DE behaviors in $P_k$ at level of $\sim 2-3\%$ and $\sim 3-4\%$ differences w.r.t $\Lambda$CDM at $z=0$ respectively. HMF shows the possibility to distinguish between the models at the high mass ends. The best-fitted model assuming only background and linear perturbations dubbed as SEoS-II largely deviates from $\Lambda$CDM and current observations on studying the nonlinear growth. This large deviation in SEoS-II also quantified the combined effect of the dynamical DE and the larger amount of matter contained, $\Omega_{m0}$ and $H_{0}$ accordingly. 2PCF results are relatively robust with $\sim 1-2 \%$ deviation for SEoS-I and CPL-lim and a significant deviation for SEoS-II throughout $r$ from $\Lambda$CDM. Finally, we conclude that the search for viable DE models (like the SEoS) must include non-linear growth constraints., Comment: 23 pages, 8 figues, 2 tables

We investigate the dynamics of axion-like particle (ALP) dark matter where the field range is enlarged by a monodromy. The monodromy potential allows sufficient production of dark matter also at larger couplings to the Standard Model particles. The potential typically features a number of "wiggles" that lead to a rapid growth of fluctuations. Using classical-statistical field theory simulations we go beyond the linear regime and treat the system in the non-linear and even non-perturbative regime. For sufficiently strong wiggles the initially homogeneous field is completely converted into fluctuations. The fluctuations correspond to dark matter particles with a non-vanishing velocity and we consider the corresponding restrictions from structure formation as well as the effects on today's dark matter density. Since all the dark matter is made up from these strong fluctuations, the dark matter density features large, $\mathcal{O}(1)$ fluctuations at scales $\lesssim 10^{6}\,{\rm km}\sqrt{{\rm eV}/m_a}$., Comment: 26 pages, 8 figures; v2: extended discussion about the field ending up in the true vs false minimum, extended discussion about the field becoming inhomogeneous vs staying homogeneous, Fig. 8 updated, clarifications and references added, published version

The origin of supermassive black holes (SMBHs) that inhabit the centers of massive galaxies is largely unconstrained. Remnants from supermassive stars (SMSs) with masses around 10,000 solar masses provide the ideal seed candidates, known as direct collapse black holes. However, their very existence and formation environment in the early Universe are still under debate, with their supposed rarity further exacerbating the problem of modeling their ab-initio formation. SMS models have shown that rapid collapse, with an infall rate above a critical value, in metal-free haloes is a requirement for the formation of a proto-stellar core which will then form an SMS. Using a radiation hydrodynamics simulation of early galaxy formation, we show the natural emergence of metal-free haloes both massive enough, and with sufficiently high infall rates, to form an SMS. We find that haloes that are exposed to both a Lyman-Werner intensity of J_LW ~ 3 J_21 and that undergo at least one period of rapid growth early in their evolution are ideal cradles for SMS formation. This rapid growth induces substantial dynamical heating, amplifying the existing Lyman-Werner suppression originating from a group of young galaxies 20 kiloparsecs away. Our results strongly indicate that structure formation dynamics, rather than a critical Lyman-Werner (LW) flux, may be the main driver of massive black hole formation in the early Universe. We find that massive black hole seeds may be much more common in overdense regions of the early Universe than previously considered with a comoving number density up to 10^-3 Mpc^-3., Comment: To appear in Nature on 23 Jan 2019 for Advanced Online Publication. Authors' latest version. 13 pages in Nature print format, 13 figures, 1 table

OB associations are the prime star forming sites in galaxies. However the detailed formation process of such stellar systems still remains a mystery. In this context, identifying the presence of substructures may help tracing the footprints of their formation process. Here, we present a kinematic study of the two massive OB associations Cygnus OB2 and Carina OB1 using the precise astrometry from the Gaia Data Release 2 and radial velocities. From the parallaxes of stars, these OB associations are confirmed to be genuine stellar systems. Both Cygnus OB2 and Carina OB1 are composed of a few dense clusters and a halo which have different kinematic properties: the clusters occupy regions of 5-8 parsecs in diameter and display small dispersions in proper motion, while the halos spread over tens of parsecs with a 2-3 times larger dispersions in proper motion. This is reminiscent of the so-called "line width-size" relation of molecular clouds related to turbulence. Considering that the kinematics and structural features were inherited from those of their natal clouds would then imply that the formation of OB associations may result from structure formation driven by supersonic turbulence, rather than from the dynamical evolution of individual embedded clusters., Comment: 20 pages, 12 figures, 4 tables, accepted for publication in MNRAS

The emission in the hyperfine structure 21 cm line of atomic hydrogen arising in the halos with masses $\sim10^6-10^{10}$ M$_\odot$ from the Dark Ages in the models with Warm Dark Matter (WDM) is analysed. The halos are assumed to be formed from Gaussian density peaks of cosmological density perturbations at $10\lesssim z\lesssim50$. Semi-analytical modelling of the formation of individual spherical halos in multi-component models shows that gas in them has the kinetic temperature in the range of $60-800$ K under adiabatic compression of the collapsing halo, and the temperature of each halo depends on the time of virialization. It is shown that inelastic collisions between neutral hydrogen atoms are the dominant excitation mechanism for hyperfine structure levels, which pulls the spin temperature closer to the kinetic temperature. The brightness temperature of individual halos is in the range of 1-10 K, depending on the mass of the halo and its virialization redshift, and increasing as these two increase. The apparent angular radii of such halos are in the range 0.06-1.2 arcseconds, their surface number density decreasing exponentially from a few per arcmin$^2$ for the lowest mass and redshift to nearly zero for higher values. Assuming a 1 MHz observation bandwidth the surface number density of the halo at various redshifts is evaluated as well as beam-averaged differential antenna temperatures and fluxes of hydrogen emission from halos of different masses. The beam-averaged signal strongly depends on the cut-off scale in the mass function of dark ages halos that may be caused by free-streaming of WDM particles. The finding is compared with the upper limits on the amplitude of the power spectrum of the hydrogen 21-cm line fluctuations derived from the recent observation data obtained with MWA and LOFAR., Comment: 18 pages, 10 figures, 7 tables; accepted for publication in Physics of the Dark Universe

Bose-Einstein Condensate Dark Matter (BECDM; also known as Fuzzy Dark Matter) is motivated by fundamental physics and has recently received significant attention as a serious alternative to the established Cold Dark Matter (CDM) model. We perform cosmological simulations of BECDM gravitationally coupled to baryons and investigate structure formation at high redshifts ($z \gtrsim 5$) for a boson mass $m=2.5\cdot 10^{-22}~{\rm eV}$, exploring the dynamical effects of its wavelike nature on the cosmic web and the formation of first galaxies. Our BECDM simulations are directly compared to CDM as well as to simulations where the dynamical quantum potential is ignored and only the initial suppression of the power spectrum is considered -- a Warm Dark Matter-like ("WDM") model often used as a proxy for BECDM. Our simulations confirm that "WDM" is a good approximation to BECDM on large cosmological scales even in the presence of the baryonic feedback. Similarities also exist on small scales, with primordial star formation happening both in isolated haloes and continuously along cosmic filaments; the latter effect is not present in CDM. Global star formation and metal enrichment in these first galaxies are delayed in BECDM/"WDM" compared to the CDM case: in BECDM/"WDM" first stars form at $z\sim 13$/$13.5$ while in CDM star formation starts at $z\sim 35$. The signature of BECDM interference, not present in "WDM", is seen in the evolved dark matter power spectrum: although the small scale structure is initially suppressed, power on kpc scales is added at lower redshifts. Our simulations lay the groundwork for realistic simulations of galaxy formation in BECDM., Comment: 19 pages, 12 figures, submitted to MNRAS