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.

In hierarchical models of structure formation, the first galaxies form in low-mass dark matter potential wells, probing the behavior of dark matter on kiloparsec (kpc) scales. Even though these objects are below the detection threshold of current telescopes, future missions will open an observational window into this emergent world. In this Letter we investigate how the first galaxies are assembled in a `fuzzy' dark matter (FDM) cosmology where dark matter is an ultralight $\sim 10^{-22}$~eV boson and the primordial stars are expected to form along dense dark matter filaments. Using a first-of-its-kind cosmological hydrodynamical simulation, we explore the interplay between baryonic physics and unique wavelike features inherent to FDM. In our simulation, the dark matter filaments show coherent interference patterns on the boson de Broglie scale and develop cylindrical soliton-like cores which are unstable under gravity and collapse into kpc-scale spherical solitons. Features of the dark matter distribution are largely unaffected by the baryonic feedback. On the contrary, the distributions of gas and stars, which do form along the entire filament, exhibit central cores imprinted by dark matter -- a smoking gun signature of FDM.

We present a suite of six fully cosmological, three-dimensional simulations of the collapse of an atomic cooling halo in the early Universe. We use the moving-mesh code arepo with an improved primordial chemistry network to evolve the hydrodynamical and chemical equations. The addition of a strong Lyman-Werner background suppresses molecular hydrogen cooling and permits the gas to evolve nearly isothermally at a temperature of about 8000 K. Strong gravitational torques effectively remove angular momentum and lead to the central collapse of gas, forming a supermassive protostar at the center of the halo. We model the protostar using two methods: sink particles that grow through mergers with other sink particles, and a stiff equation of state that leads to the formation of an adiabatic core. We impose threshold densities of $10^8$, $10^{10}$, and $10^{12}\,\text{cm}^{-3}$ for the sink particle formation and the onset of the stiff equation of state to study the late, intermediate, and early stages in the evolution of the protostar, respectively. We follow its growth from masses $\simeq 10\,\text{M}_\odot$ to $\simeq 10^5\,\text{M}_\odot$, with an average accretion rate of $\langle\dot{M}_\star\rangle \simeq 2\,\text{M}_\odot\,\text{yr}^{-1}$ for sink particles, and $\simeq 0.8 - 1.4\,\text{M}_\odot\,\text{yr}^{-1}$ for the adiabatic cores. At the end of the simulations, the HII region generated by radiation from the central object has long detached from the protostellar photosphere, but the ionizing radiation remains trapped in the inner host halo, and has thus not yet escaped into the intergalactic medium. Fully coupled, radiation-hydrodynamics simulations hold the key for further progress.

The Schr\"odinger-Poisson equations describe the behavior of a superfluid Bose-Einstein condensate under self-gravity with a 3D wave function. As $\hbar/m\to 0$, $m$ being the boson mass, the equations have been postulated to approximate the collisionless Vlasov-Poisson equations also known as the collisionless Boltzmann-Poisson equations. The latter describe collisionless matter with a 6D classical distribution function. We investigate the nature of this correspondence with a suite of numerical test problems in 1D, 2D, and 3D along with analytic treatments when possible. We demonstrate that, while the density field of the superfluid always shows order unity oscillations as $\hbar/m\to 0$ due to interference and the uncertainty principle, the potential field converges to the classical answer as $(\hbar/m)^{2}$. Thus, any dynamics coupled to the superfluid potential is expected to recover the classical collisionless limit as $\hbar/m\to 0$. The quantum superfluid is able to capture rich phenomena such as multiple phase-sheets, shell-crossings, and warm distributions. Additionally, the quantum pressure tensor acts as a regularizer of caustics and singularities in classical solutions. This suggests the exciting prospect of using the Schr\"odinger-Poisson equations as a low-memory method for approximating the high-dimensional evolution of the Vlasov-Poisson equations. As a particular example we consider dark matter composed of ultra-light axions, which in the classical limit ($\hbar/m\to 0$) is expected to manifest itself as collisionless cold dark matter.

We present a model for the evolution of supermassive protostars from their formation at $M_\star \simeq 0.1\,\text{M}_\odot$ until their growth to $M_\star \simeq 10^5\,\text{M}_\odot$. To calculate the initial properties of the object in the optically thick regime we follow two approaches: based on idealized thermodynamic considerations, and on a more detailed one-zone model. Both methods derive a similar value of $n_{\rm F} \simeq 2 \times 10^{17} \,\text{cm}^{-3}$ for the density of the object when opacity becomes important, i.e. the opacity limit. The subsequent evolution of the growing protostar is determined by the accretion of gas onto the object and can be described by a mass-radius relation of the form $R_\star \propto M_\star^{1/3}$ during the early stages, and of the form $R_\star \propto M_\star^{1/2}$ when internal luminosity becomes important. For the case of a supermassive protostar, this implies that the radius of the star grows from $R_\star \simeq 0.65 \,{\rm AU}$ to $R_\star \simeq 250 \,{\rm AU}$ during its evolution. Finally, we use this model to construct a sub-grid recipe for accreting sink particles in numerical simulations. A prime ingredient thereof is a physically motivated prescription for the accretion radius and the effective temperature of the growing protostar embedded inside it. From the latter, we can conclude that photo-ionization feedback can be neglected until very late in the assembly process of the supermassive object.

We perform a post-processing radiative feedback analysis on a 3D ab initio cosmological simulation of an atomic cooling halo under the direct collapse black hole (DCBH) scenario. We maintain the spatial resolution of the simulation by incorporating native ray-tracing on unstructured mesh data, including Monte Carlo Lyman-alpha (Ly{\alpha}) radiative transfer. DCBHs are born in gas-rich, metal-poor environments with the possibility of Compton-thick conditions, $N_H \gtrsim 10^{24} {\rm cm}^{-2}$. Therefore, the surrounding gas is capable of experiencing the full impact of the bottled-up radiation pressure. In particular, we find that multiple scattering of Ly{\alpha} photons provides an important source of mechanical feedback after the gas in the sub-parsec region becomes partially ionized, avoiding the bottleneck of destruction via the two-photon emission mechanism. We provide detailed discussion of the simulation environment, expansion of the ionization front, emission and escape of Ly{\alpha} radiation, and Compton scattering. A sink particle prescription allows us to extract approximate limits on the post-formation evolution of the radiative feedback. Fully coupled Ly{\alpha} radiation hydrodynamics will be crucial to consider in future DCBH simulations.

We study the star formation process at galactic scales and the role of rotation through numerical simulations of spiral and starburst galaxies using the adaptive mesh refinement code Enzo. We focus on the study of three integrated star formation laws found in the literature: the Kennicutt-Schmidt (KS) and Silk-Elmegreen (SE) laws, and the dimensionally homogeneous equation proposed by Escala (2015) $\Sigma_{\rm SFR} \propto \sqrt{G/L}\Sigma_{\rm gas}^{1.5}$. We show that using the last we take into account the effects of the integration along the line of sight and find a unique regime of star formation for both types of galaxies, suppressing the observed bi-modality of the KS law. We find that the efficiencies displayed by our simulations are anti-correlated with the angular velocity of the disk $\Omega$ for the three laws studied in this work. Finally, we show that the dimensionless efficiency of star formation is well represented by an exponentially decreasing function of $-1.9\Omega t_{\rm ff}^{\rm ini}$, where $t_{\rm ff}^{\rm ini}$ is the initial free-fall time. This leads to a unique galactic star formation relation which reduces the scatter of the bi-modal KS, SE, and Escala (2015) relations by 43%, 43%, and 35% respectively.

We present the highest-resolution three-dimensional simulation to date of the collapse of an atomic cooling halo in the early Universe. We use the moving-mesh code arepo with the primordial chemistry module introduced in Greif (2014), which evolves the chemical and thermal rate equations for over more than 20 orders of magnitude in density. Molecular hydrogen cooling is suppressed by a strong Lyman-Werner background, which facilitates the near-isothermal collapse of the gas at a temperature of about $10^4\,$K. Once the central gas cloud becomes optically thick to continuum emission, it settles into a Keplerian disc around the primary protostar. The initial mass of the protostar is about $0.1\,{\rm M}_\odot$, which is an order of magnitude higher than in minihaloes that cool via molecular hydrogen. The high accretion rate and efficient cooling of the gas catalyse the fragmentation of the disc into a small protostellar system with 5-10 members. After about 12 yr, strong gravitational interactions disrupt the disc and temporarily eject the primary protostar from the centre of the cloud. By the end of the simulation, a secondary clump has collapsed at a distance of $\simeq 150\,$au from the primary clump. If this clump undergoes a similar evolution as the first, the central gas cloud may evolve into a wide binary system. High accretion rates of both the primary and secondary clumps suggest that fragmentation is not a significant barrier for forming at least one massive black hole seed.

We use the Adaptive Mesh Refinement code Enzo to model the interstellar medium in isolated local disk galaxies. The simulation includes a treatment for star formation and stellar feedback. We get a highly supersonic turbulent disk, which is fragmented at multiple scales and characterized by a multi-phase interstellar medium. We show that a Kennicutt-Schmidt (KS) relation only holds when averaging over large scales. However, values of star formation rates and gas surface densities lie close in the plot for any averaging size. This suggests an intrinsic relation between stars and gas at cell-size scales, which dominates over the global dynamical evolution. To investigate this effect, we develop a method to simulate the creation of stars based on the density field from the snapshots, without running the code again. We also investigate how the star formation law is affected by the characteristic star formation timescale, the density threshold and the efficiency considered in the recipe. We find that the slope of the law might vary from ~1.4 for a free-fall timescale, to ~1.0 for a constant depletion timescale. We further demonstrate that a power-law is recovered just by assuming that the mass of the new stars is a fraction of the mass of the cell $m_\star=\epsilon\rho_{\rm gas}\Delta x^3$, with no other physical criteria required. We show that both efficiency and density threshold do not affect the slope, but the right combination of them can adjust the normalization of the relation, which in turn could explain a possible bi-modality in the law.

We study the gravitational stability of gaseous streams in the complex environment of a galaxy merger, because mergers are known to be places of ongoing massive cluster formation and bursts of star formation. We find an analytic stability parameter for case of gaseous streams orbiting around the merger remnant. We test our stability criteria using hydrodynamical simulations of galaxy mergers, obtaining satisfactory results. We find that our criteria successfully predicts the streams that will be gravitationally unstable to fragment into clumps.