Search

In the context of Fuzzy Dark Matter (FDM) we study the core formation in the presence of an Ideal Gas (IG). Our analysis is based on the solution of the Schr\"odinger-Poisson-Euler system of equations that drives the evolution of FDM together with a compressible IG, both coupled through the gravitational potential they produce. Starting from random initial conditions for both FDM and IG, with dominant FDM, we study the evolution of the system until it forms a nearly relaxed, virialized and close to hydrostatic equilibrium core, surrounded by an envelope of the two components. We find that the core corresponds to Newtonian Fermion-Boson Stars (FBS). If the IG is used to model luminous matter, our results indicate that FBS behave as attractor core solutions of structure formation of FDM along with baryonic matter., Comment: 6 pages, 5 figures

The post-inflationary Universe can pass through a long epoch of effective matter-dominated expansion. This era may allow for both the parametric amplification of initial fluctuations and the gravitational collapse of inflaton perturbations. We perform first-of-their-kind high-resolution simulations that span the resonant phase and the subsequent gravitational collapse of the inflaton field by seguing from a full Klein-Gordon treatment of resonance to a computationally efficient Schr\"odinger-Poisson description that accurately captures the gravitational dynamics when most quanta are nonrelativistic. We consider a representative example in which resonance generates $\mathcal{O}(10^{-1})$ overdensities and gravitational collapse follows promptly as resonance ends. We observe the formation of solitonic cores inside inflaton halos and complex gravitational dynamics on scales of $10^{-27}\,\mathrm{m}$, greatly extending the possible scope of nonlinear post-inflationary gravitational dynamics., Comment: 11 pages, 8 figures, matches version accepted in Phys Rev D

Following inflation, the Universe may pass through an early matter-dominated phase supported by the oscillating inflaton condensate. Initially small fluctuations in the condensate grow gravitationally on subhorizon scales and can collapse to form nonlinear ``inflaton halos''. Their formation and subsequent tidal interactions will source gravitational waves, resulting in a stochastic background in the present Universe. We extend N-body simulations that model the growth and interaction of collapsed structures to compute the resulting gravitational wave emission. The spectrum of this radiation is well-matched by semi-analytical estimates based on the collapse of inflaton halos and their tidal evolution. We use this semi-analytic formalism to infer the spectrum for scenarios where the early matter-dominated phase gives way to a thermalized universe at temperatures as low as $100\,\mathrm{MeV}$ and we discuss the possible experimental opportunities created by this signal in inflationary models in which thermalization takes place long after inflation has completed., Comment: 9 pages, 4 figures. Published version, Editors' suggestion

Primordial magnetic fields (PMFs) are possible candidates for explaining the observed magnetic fields in galaxy clusters. Two competing scenarios of primordial magnetogenesis have been discussed in the literature: inflationary and phase-transitional. We study the amplification of both large- and small-scale correlated magnetic fields, corresponding to inflation- and phase transition-generated PMFs, in a massive galaxy cluster. We employ high-resolution magnetohydrodynamic cosmological zoom-in simulations to resolve the turbulent motions in the intracluster medium. We find that the turbulent amplification is more efficient for the large-scale inflationary models, while the phase transition-generated seed fields show moderate growth. The differences between the models are imprinted on the spectral characteristics of the field (such as the amplitude and the shape of the magnetic power spectrum) and therefore on the final correlation length. We find a one order of magnitude difference between the final strengths of the inflation- and phase transition-generated magnetic fields, and a factor of 1.5 difference between their final coherence scales. Thus, the final configuration of the magnetic field retains information about the PMF generation scenarios. Our findings have implications for future extragalactic Faraday rotation surveys with the possibility of distinguishing between different magnetogenesis scenarios., Comment: 21 pages, 13 figures, comments welcome

The Peak-Patch algorithm is used to identify the densest minicluster seeds in the initial axion density field simulated from string decay. The fate of these dense seeds is found by tracking the subsequent gravitational collapse in cosmological $N$-body simulations. We find that miniclusters at late times are well described by NFW profiles, although for around 80\% of simulated miniclusters a single power-law density profile of $r^{-2.9}$ is an equally good fit due to the unresolved scale radius. Under the assumption that all miniclusters with an unresolved scale radius are described by a power-law plus axion star density profile, we identify a significant number of miniclusters that might be dense enough to give rise to gravitational microlensing if the axion mass is $0.2 \,\mathrm{meV}\lesssim m_a \lesssim 3\,\mathrm{meV}$. Higher resolution simulations resolving the inner structure and axion star formation are necessary to explore this possibility further., Comment: 10+8 pages, 15 figures, v2 further comments on density profile fits, matches version accepted in Phys Rev D

The Universe may pass through an effectively matter-dominated epoch between inflation and Big Bang Nucleosynthesis during which gravitationally bound structures can form on subhorizon scales. In particular, the inflaton field can collapse into inflaton halos, forming "large scale" structure in the very early universe. We combine N-body simulations with high-resolution zoom-in regions in which the non-relativistic Schr\"odinger-Poisson equations are used to resolve the detailed, wave-like structure of inflaton halos. Solitonic cores form inside them, matching structure formation simulations with axion-like particles in the late-time universe. We denote these objects \textit{inflaton stars}, by analogy with boson stars. Based on a semi-analytic formalism we compute their overall mass distribution which shows that some regions will reach overdensities of $10^{15}$ if the early matter-dominated epoch lasts for 20 $e$-folds. The radii of the most massive inflaton stars can shrink below the Schwarzschild radius, suggesting that they could form primordial black holes prior to thermalization., Comment: 9 pages, 6 figures, accepted by PRD

Fuzzy dark matter (FDM) made of ultra-light bosonic particles is a viable alternative to cold dark matter (CDM) with clearly distinguishable small-scale features in collapsed structures. On large scales, it behaves gravitationally like CDM deviating only by a cut-off in the initial power spectrum and can be studied using N-body methods. In contrast, wave interference effects near the de Broglie scale result in new phenomena unique to FDM. Interfering modes in filaments and halos yield a stochastically oscillating granular structure which condenses into solitonic cores during halo formation. Investigating these highly non-linear wave phenomena requires the spatially resolved numerical integration of the Schr\"odinger equation. In previous papers we introduced a hybrid zoom-in scheme that combines N-body methods to model the large-scale gravitational potential around and the mass accretion onto pre-selected halos with simulations of the Schr\"odinger-Poisson equation to capture wave-like effects inside these halos. In this work, we present a new, substantially improved reconstruction method for the wave function inside of previously collapsed structures. We demonstrate its capabilities with a deep zoom-in simulation of a well-studied sub-$L_\ast$-sized galactic halo from cosmological intitial conditions. With a particle mass of $m = 2.5\times 10^{-22}\,$eV and halo mass $M_{\text{vir}}=1.7\times 10^{11}\,M_{\odot}$ in a ($60$h${^{-1}}$ comoving Mpc)${}^{3}$ cosmological box, it reaches an effective resolution of 20 comoving pc. This pushes the values of $m$ and $M$ accessible to simulations significantly closer to those relevant for studying galaxy evolution in the allowed range of FDM masses., Comment: 5 pages, 5 figures

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

This work studies the formation and growth of boson stars and their surrounding miniclusters by gravitational condensation using non-linear dynamical numerical methods. Fully dynamical attractive and repulsive self-interactions are also considered for the first time. In the case of pure gravity, we numerically prove that the growth of boson stars inside halos slows down and saturates as has been previously conjectured, and detail its conditions. Self-interactions are included using the Gross-Pitaevskii-Poisson equations. We find that in the case of strong attractive self-interactions the boson stars can become unstable and collapse, in agreement with previous stationary computations. At even stronger coupling, the condensate fragments. Repulsive self-interactions, as expected, promote boson star formation, and lead to solutions with larger radii., Comment: 15 pages, 21 figures

The distinctive effects of fuzzy dark matter are most visible at non-linear galactic scales. We present the first simulations of mixed fuzzy and cold dark matter, obtained with an extended version of the Nyx code. Fuzzy (or ultralight, or axion-like) dark matter dynamics are governed by the comoving Schr\"odinger-Poisson equation. This is evolved with a pseudospectral algorithm on the root grid, and with finite differencing at up to six levels of adaptive refinement. Cold dark matter is evolved with the existing N-body implementation in Nyx. We present the first investigations of spherical collapse in mixed dark matter models, focusing on radial density profiles, velocity spectra and soliton formation in collapsed halos. We find that the effective granule masses decrease in proportion to the fraction of fuzzy dark matter which quadratically suppresses soliton growth, and that a central soliton only forms if the fuzzy dark matter fraction is greater than 10\%. The Nyx framework supports baryonic physics and key astrophysical processes such as star formation. Consequently, AxioNyx will enable increasingly realistic studies of fuzzy dark matter astrophysics.

Axion-like particle (ALP) dark matter shows distinctive behavior on scales where wavelike effects dominate over self-gravity. Ultralight axions are candidates for fuzzy dark matter (FDM) whose de Broglie wavelength in virialized halos reaches scales of kiloparsecs. Important features of FDM scenarios are the formation of solitonic halo cores, suppressed small-scale perturbations, and enhanced gravitational relaxation. More massive ALPs, including the QCD axion, behave like CDM on galactic scales but may be clumped into axion miniclusters if they were produced after inflation. Just as FDM halos, axion miniclusters may host the formation of coherent bound objects (axion stars) by Bose-Einstein condensation. This article presents a selection of topics in this field that are currently under active investigation., Comment: To be published in Progress in Particle and Nuclear Physics 2020

We study the gravitational collapse of axion dark matter fluctuations in the post-inflationary scenario, so-called axion miniclusters, with N-body simulations. Largely confirming theoretical expectations, overdensities begin to collapse in the radiation-dominated epoch and form an early distribution of miniclusters with masses up to $10^{-12}\,M_\odot$. After matter-radiation equality, ongoing mergers give rise to a steep power-law distribution of minicluster halo masses. The density profiles of well-resolved halos are NFW-like to good approximation. The fraction of axion dark matter in these bound structures is $\sim 0.75$ at redshift $z=100$., Comment: 8 pages, 6 figures, accepted by PRL

We present zoom-in simulations of fuzzy dark matter (FDM) halos including baryons and star-formation with sufficient resolution to follow the formation and evolution of central solitons. We find that their properties are determined by the local dark matter velocity dispersion in the combined dark matter-baryon gravitational potential. This motivates a simple prescription to estimate the radial density profiles of FDM cores in the presence of baryons. As cores become more massive and compact if baryons are included, galactic rotation curve measurements are likely harder to reconcile with FDM., Comment: 6 pages, 6 figures, matches published version

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

We study the formation and the subsequent mass growth of axion stars inside axion miniclusters. Numerically solving the Schr\"odinger-Poisson equations with realistic initial conditions we find that the axion stars exhibit similar properties to solitonic cores in ultralight bosonic dark matter halos in terms of their radial density profiles and large-amplitude oscillations. A merger of two axion stars confirms a previously found empirical law for the mass of the merged axion star. Monitoring the axion star masses over time, we observe a mass growth consistent with the mass increase of Bose stars in the kinetic regime reported by Levkov et al., confirming that the mass evolution of axion stars can be understood in terms of wave condensation. Based on this result, we predict a saturation of mass growth in relation to the minicluster mass consistent with the core-halo mass relation previously found for ultra-light bosonic dark matter halos., Comment: 7 pages, 5 figures, accepted by PRD

The fuzzy dark matter (FDM) model treats DM as a bosonic field with astrophysically large de Broglie wavelength. A striking feature of this model is $\mathcal{O}(1)$ fluctuations in the dark matter density on time scales which are shorter than the gravitational timescale. Including for the first time the effect of core oscillations, we demonstrate how such fluctuations lead to heating of star clusters, and thus an increase in their size over time. From the survival of the old star cluster in Eridanus II we infer $m_a\gtrsim 0.6\rightarrow 1\times 10^{-19}\text{ eV}$ within modelling uncertainty if FDM is to compose all of the DM, and derive constraints on the FDM fraction at lower masses. The subhalo mass function in the Milky Way implies $m_a\gtrsim 0.8\times 10^{-21}\text{ eV}$ to successfully form Eridanus II. The window between $10^{-21}\text{ eV}\lesssim m_a\lesssim 10^{-20}\text{ eV}$ is affected by narrow band resonances, and the limited applicability of the diffusion approximation. Some of this window may be consistent with observations of Eridanus II and more detailed investigations are required., Comment: 6 pages, 2 figures; v2 comments and clarifications added, matches version accepted for publication in PRL

Axions are increasingly favoured as a candidate particle for the dark matter in galaxies, since they satisfy the observational requirements for cold dark matter and are theoretically well motivated. Fluctuations in the axion field give rise to stable localised overdensities known as axion stars, which, for the most massive, compact cases, are potential neutron star mimickers. In principle, there are no fundamental arguments against the multi-messenger observations of GW170817/GRB170817A/AT2017gfo arising from the merger of a neutron star with a neutron star mimicker, rather than from a binary neutron star. To constrain this possibility and better understand the astrophysical signatures of a neutron star--axion star (NSAS) merger, we present in this work a detailed example case of a NSAS merger based on full 3D numerical relativity simulations, and give an overview of the many potential observables - ranging from gravitational waves, to optical and near-infrared electromagnetic signals, radio flares, fast radio bursts, gamma ray bursts, and neutrino emission. We discuss the individual channels and estimate to which distances current and future observatories might be able to detect such a NSAS merger. Such signals could constrain the unknown axion mass and its couplings to standard baryonic matter,thus enhancing our understanding of the dark matter sector of the Universe.

Axions are a potential dark matter candidate, which may condense and form self gravitating compact objects, called axion stars (ASs). In this work, we study for the first time head-on collisions of relativistic ASs with black holes (BHs) and neutron stars (NSs). In the case of BH-AS mergers we find that, in general, the largest scalar clouds are produced by mergers of low compactness ASs and spinning BHs. Although in most of the cases which we study the majority of the mass is absorbed by the BH within a short time after the merger, in favourable cases the remaining cloud surrounding the final BH remnant can be as large as 30% of the initial axion star mass, with a bosonic cloud mass of order ~ 0.1 M_BH and peak energy density comparable to that obtained in a superradiant build up. This provides a dynamical mechanism for the formation of long lived scalar hair, which could lead to observable signals in cases where the axion interacts with baryonic matter around the BH, or where it forms the seed of a future superradiant build up in highly spinning cases. Considering NS-AS collisions we find two possible final states (i) a BH surrounded by a (small) scalar cloud, or (ii) a stable NS enveloped in an axion cloud of roughly the same mass as the initial AS. Whilst for low mass ASs the NS is only mildly perturbed by the collision, a larger mass AS gives rise to a massive ejection of baryonic mass from the system, purely due to gravitational effects. Therefore, even in the absence of a direct axion coupling to baryonic matter, NS-AS collisions could give rise to electromagnetic observables in addition to their gravitational wave signatures., Comment: 13 pages, 14 figures

We simulate the formation and evolution of ultralight bosonic dark matter halos from cosmological initial conditions. Using zoom-in techniques we are able to resolve the detailed interior structure of the halos. We observe the formation of solitonic cores and confirm the core-halo mass relation previously found by Schive et al. The cores exhibit strong quasi-normal oscillations that remain largely undamped on evolutionary timescales. On the other hand, no conclusive growth of the core mass by condensation or relaxation can be detected. In the incoherent halo surrounding the cores, the scalar field density profiles and velocity distributions show no significant deviation from collisionless N-body simulations on scales larger than the coherence length. Our results are consistent with the core properties being determined mainly by the coherence length at the time of virialization, whereas the Schr\"odinger-Vlasov correspondence explains the halo properties when averaged on scales greater than the coherence length., Comment: 8 pages, 8 figures, accepted by PRD

A recent article by Galtier and Nazarenko [1] proposed that weakly nonlinear gravitational waves could result in a turbulent cascade, with energy flowing from high to low frequency modes or vice versa. This is an interesting proposition for early universe cosmology because it could suggest some "natural" initial conditions for the gravitational background. In this paper we use the ADM formalism to show that, given some simple and, arguably, natural assumptions, such initial conditions lead to expansion (or collapse) of the spacetime on a timescale much faster than that of the turbulent cascade, meaning that the cascade is unlikely to have sufficient time to develop under general conditions. We suggest possible ways in which the expansion could be mitigated to give the cascade time to develop., Comment: 4 pages

We study tidal stripping of fuzzy dark matter (FDM) subhalo cores using simulations of the Schr\"{o}dinger-Poisson equations and analyze the dynamics of tidal disruption, highlighting the differences with standard cold dark matter. Mass loss outside of the tidal radius forces the core to relax into a less compact configuration, lowering the tidal radius. As the characteristic radius of a solitonic core scales inversely with its mass, tidal stripping results in a runaway effect and rapid tidal disruption of the core once its central density drops below $4.5$ times the average density of the host within the orbital radius. Additionally, we find that the core is deformed into a tidally locked ellipsoid with increasing eccentricities until it is completely disrupted. Using the core mass loss rate, we compute the minimum mass of cores that can survive several orbits for different FDM particle masses and compare it with observed masses of satellite galaxies in the Milky Way., Comment: 11 pages, 12 figures, accepted for publication in PRD

Lyman-$\alpha$ emitters (LAEs) are a promising probe of the large-scale structure at high redshift, $z\gtrsim 2$. In particular, the Hobby-Eberly Telescope Dark Energy Experiment aims at observing LAEs at 1.9 $

In the context of structure formation with ultralight axion dark matter, we offer an alternative explanation for the mass relation of solitonic cores and their host halos observed in numerical simulations. Our argument is based entirely on the mass gain that occurs during major mergers of binary cores and largely independent of the initial core-halo mass relation assigned to hosts that have just collapsed. We find a relation between the halo mass $M_h$ and corresponding core mass $M_c$, $M_c\propto M_h^{2\beta-1}$, where $(1-\beta)$ is the core mass loss fraction. Following the evolution of core masses in stochastic merger trees, we find empirical evidence for our model. Our results are useful for statistically modeling the effects of dark matter cores on the properties of galaxies and their substructures in axion dark matter cosmologies., Comment: 6 pages, 5 figures, accepted for publication in PRD

We derive the halo mass function (HMF) for fuzzy dark matter (FDM) by solving the excursion set problem explicitly with a mass-dependent barrier function, which has not been done before. We find that compared to the naive approach of the Sheth--Tormen HMF for FDM, our approach has a higher cut off mass and the cut off mass changes less strongly with redshifts. Using merger trees constructed with a modified version of the Lacey & Cole formalism that accounts for suppressed small scale power and the scale-dependent growth of FDM haloes and the semi-analytic GALACTICUS code, we study the statistics of halo substructure including the effects from dynamical friction and tidal stripping. We find that if the dark matter is a mixture of cold dark matter (CDM) and FDM, there will be a suppression on the halo substructure on small scales which may be able to solve the Missing Satellites Problem faced by the pure CDM model. The suppression becomes stronger with increasing FDM fraction or decreasing FDM mass. Thus, it may be used to constrain the FDM model., Comment: 12 pages, 18 figures, accepted for publication in MNRAS

We study cosmological structure formation with ultralight axion dark matter, or "fuzzy dark matter (FDM), using a particle-mesh scheme to account for the quantum pressure arising in the Madelung formulation of the Schr\"odinger-Poisson equations. Subpercent-level energy conservation and correct linear behavior are demonstrated. Whereas the code gives rise to the same core-halo profiles as direct simulations of the Schr\"odinger equation, it does not reproduce the detailed interference patterns. In cosmological simulations with FDM initial conditions, we find a maximum relative difference of O($10\%$) in the power spectrum near the quantum Jeans length compared to using a standard N-body code with identical initial conditions. This shows that the effect of quantum pressure during nonlinear structure formation cannot be neglected for precision constraints on a dark matter component consisting of ultralight axions., Comment: 7 pages, 8 figures, this version was accepted in PRD

Using three-dimensional simulations, we study the dynamics and final structure of merging solitonic cores predicted to form in ultra-light axion dark matter halos. The classical, Newtonian equations of motion of a self-gravitating scalar field are described by the Schr\"odinger-Poisson equations. We investigate mergers of ground state (boson star) configurations with varying mass ratios, relative phases, orbital angular momenta and initial separation with the primary goal to understand the mass loss of the emerging core by gravitational cooling. Previous results showing that the final density profiles have solitonic cores and NFW-like tails are confirmed. In binary mergers, the final core mass does not depend on initial phase difference or angular momentum and only depends on mass ratio, total initial mass, and total energy of the system. For non-zero angular momenta, the otherwise spherical cores become rotating ellipsoids. The results for mergers of multiple cores are qualitatively identical., Comment: 10 pages, 22 figures

The gas in galaxy clusters is heated by shock compression through accretion (outer shocks) and mergers (inner shocks). These processes additionally produce turbulence. To analyse the relation between the thermal and turbulent energies of the gas under the influence of non-adiabatic processes, we performed numerical simulations of cosmic structure formation in a box of 152 Mpc comoving size with radiative cooling, UV background, and a subgrid scale model for numerically unresolved turbulence. By smoothing the gas velocities with an adaptive Kalman filter, we are able to estimate bulk flows toward cluster cores. This enables us to infer the velocity dispersion associated with the turbulent fluctuation relative to the bulk flow. For halos with masses above $10^{13}\,M_\odot$, we find that the turbulent velocity dispersions averaged over the warm-hot intergalactic medium (WHIM) and the intracluster medium (ICM) are approximately given by powers of the mean gas temperatures with exponents around 0.5, corresponding to a roughly linear relation between turbulent and thermal energies and transonic Mach numbers. However, turbulence is only weakly correlated with the halo mass. Since the power-law relation is stiffer for the WHIM, the turbulent Mach number tends to increase with the mean temperature of the WHIM. This can be attributed to enhanced turbulence production relative to dissipation in particularly hot and turbulent clusters., Comment: 21 pages, 19 figures, accepted for publication by MNRAS

We propose an anisotropic generalisation of the line correlation function (ALCF) to separate and quantify phase information in the large-scale structure of galaxies. The line correlation function probes the strictly non-linear regime of structure formation and since phase information drops out of the power spectrum, the line correlation function provides a complementary tool to commonly used techniques based on two-point statistics. Furthermore, it is independent of linear bias as well as the Gaussian variance on the modulus of the density field and thus may also prove to be advantageous compared to the bispectrum or similar higher-order statistics for certain cases. For future applications it is vital, though, to be able to account for observational effects that cause anisotropies in the distribution of galaxies. Based on a number of numerical studies, we find that our ALCF is well suited to accomplish this task and we demonstrate how the Alcock-Paczynski effect and kinematical redshift-space distortions can in principle be measured via the ALCF., Comment: 14 pages

We present a subgrid-scale model for the Multi-phase Interstellar medium, Star formation, and Turbulence (MIST) and explore its behaviour in high-resolution large-eddy simulations of isolated disc galaxies. MIST follows the evolution of a clumpy cold and a diffuse warm component of the gas within a volume element which exchange mass and energy via various cooling, heating and mixing processes. The star formation rate is dynamically computed from the state of the gas in the cold phase. An important feature of MIST is the treatment of unresolved turbulence in the two phases and its interaction with star formation and feedback by supernovae. This makes MIST a particularly suitable model for the interstellar medium in galaxy simulations. We carried out a suite of simulations varying fundamental parameters of our feedback implementation. Several observational properties of galactic star formation are reproduced in our simulations, such as an average star formation efficiency ~1%, a typical velocity dispersion around ~10 km/s in star-forming regions, and an almost linear relationship between the column densities of star formation and dense molecular gas., Comment: 21 pages, 20 figures, accepted for publication by MNRAS

We investigate the radiative transfer of Lyman alpha photons through simplified anisotropic gas distributions, which represent physically motivated extensions of the popular 'shell-models'. Our study is motivated by the notion that (i) shell models do not always reproduce observed Lyman alpha spectral line profiles, (ii) (typical) shell models do not allow for the escape of ionizing photons, and (iii) the observation & expectation that winds are more complex, anisotropic phenomena. We examine the influence of inclination on the Lyman alpha spectra, relative fluxes and escape fractions. We find the flux to be enhanced/suppressed by factors up to a few depending on the parameter range of the models, corresponding to a boost in equivalent width of the same amplitude if we neglect dust. In general, lower mean optical depths tend to reduce the impact of anisotropies as is expected. We find a correlation between an observed peak in the -- occasionally triple-peaked -- spectrum at the systemic velocity and the existence of a low optical depth cavity along the line of sight. This can be of importance in the search for ionizing photons leaking from high- redshift galaxies since these photons will also be able to escape through the cavity., Comment: 11 pages, accepted for publication in A&A

We investigate the correlations between the observed fraction of Lyman alpha emission from star-forming galaxies and the large scale structure by post-processing snapshots of a large, high resolution hydrodynamical simulation with a Lyman alpha radiative transfer code at redshifts z=4, 3, 2. We find correlations of the observed fraction with density, density gradient along the line of sight, velocity and velocity gradient along the line of sight, all within the same order of magnitude (tens of percent). Additionally, a correlation with the angular momentum of the dark matter halo is detected. In contrast to a previous study, we find no significant deformation of the 2-point correlation function due to selection effects from radiative transfer in the IGM within the limited statistics of the simulation volume., Comment: 13 pages; accepted for publication in A&A

We derive oscillatory signals in correlation functions in two-field open inflation by means of the in-in formalism; such signatures are caused by resonances between oscillations in the tunnelling field and fluctuations in the inflaton during the curvature dominated, intermediate and subsequent inflationary regime. While amplitudes are model-dependent, we find distinct oscillations in the power and bi-spectrum that can act as a direct probe of the curvature dominated phase and thus, indirectly, strengthen the claim of the string landscape if they were observed. We comment on the prospects of detecting these tell-tale signs in current experiments, which is challenging, but not impossible. At the technical level, we pay special attention to the applicability conditions for truncation fluctuations to the light (inflaton) field and derive upper limits on the oscillation amplitude of the heavy field. A violation of these bounds requires a multi-field analysis at the perturbed level., Comment: 42 pages, 9 figures

We explore the amplification of magnetic fields in the high-redshift Universe. For this purpose, we perform high-resolution cosmological simulations following the formation of primordial halos with \sim10^7 M_solar, revealing the presence of turbulent structures and complex morphologies at resolutions of at least 32 cells per Jeans length. Employing a turbulence subgrid-scale model, we quantify the amount of unresolved turbulence and show that the resulting turbulent viscosity has a significant impact on the gas morphology, suppressing the formation of low-mass clumps. We further demonstrate that such turbulence implies the efficient amplification of magnetic fields via the small-scale dynamo. We discuss the properties of the dynamo in the kinematic and non-linear regime, and explore the resulting magnetic field amplification during primordial star formation. We show that field strengths of \sim10^{-5} G can be expected at number densities of \sim5 cm^{-3}., Comment: 7 pages, 6 figures, 1 table. Proceedings article for the 2012 Fall Meeting of the Astronomische Gesellschaft (AG 2012) in Hamburg, submitted for the yearbook series "Reviews in Modern Astronomy", volume 25, of the Astronomische Gesellschaft

In a landscape of compactifications with different numbers of macroscopic dimensions, it is possible that our universe has nucleated from a vacuum where some of our four large dimensions were compact while other, now compact, directions were macroscopic. From our perspective, this shapeshifting can be perceived as an anisotropic background spacetime. As an example, we present a model where our universe emerged from a tunneling event which involves the decompactification of two dimensions compactified on the two-sphere. In this case, our universe is of the Kantowski-Sachs type and therefore homogeneous and anisotropic. We study the deviations from statistical isotropy of the Cosmic Microwave Background induced by the anisotropic curvature, with particular attention to the anomalies. The model predicts a quadrupolar power asymmetry with the same sign and acoustic oscillations as found by WMAP. The amplitude of the effect is however too small given the current estimated bound on anisotropic curvature derived from the quadrupole., Comment: 17 pages, 4 figures. v3: extended discussion of CMB anomalies; version accepted by Phys. Rev. D

Motivated by cosmological first-order phase transitions we examine the nucleation and evolution of vacuum bubbles in non-vacuum environments. Non-standard backgrounds can be relevant in the context of rapid tunneling processes on the landscape. Utilising complex time methods, we show that tunneling rates can be notably modified in the case of dynamical FRW backgrounds. We give a classification of the importance of the effect in terms of the relevant dynamical time scales. For both the bubble nucleation and evolution analysis we make use of the thin-wall approximation. From the classical bubble evolution on homogeneous matter backgrounds via the junction method, we find that the inflation of vacuum bubbles is very sensitive to the presence of ambient matter and quantify this statement. We also employ inhomogeneous matter models (LTB) and models that undergo a rapid phase transition (FRW) as a background and discuss in which cases potentially observable imprints on the bubble trajectory can remain., Comment: 8 pages, 4 figures; proceedings (PoS) of the International Workshop on Cosmic Structure and Evolution, Bielefeld, Germany, 23-25 September 2009

In this letter we discuss cosmological first order phase transitions with de Sitter bubbles nucleating on (inhomogeneous) matter backgrounds. The de Sitter bubble can be a toy model for an inflationary phase of universes like our own. Using the thin wall approximation and the Israel junction method we trace the classical evolution of the formed bubbles within a compound model. We first address homogeneous ambient space (FRW model) and already find that bubbles nucleated in a dust dominated background cannot expand. For an inhomogeneous dust background (LTB model) we describe cases with at least initially expanding bubbles. Yet, an ensuing passage of the bubble wall through ambient curvature inhomogeneities remains unnoticed for observers inside the bubble. Notable effects also for interior observers are found in the case of a rapid background phase transition in a FRW model., Comment: 4 pages, 2 figures, proceedings (Annalen der Physik) of the Grassmannian Conference in Fundamental Cosmology (Grasscosmofun'09) Szczecin, Poland, 14-19 September 2009

In the context of bubble universes produced by a first-order phase transition with large nucleation rates compared to the inverse dynamical time scale of the parent bubble, we extend the usual analysis to non-vacuum backgrounds. In particular, we provide semi-analytic and numerical results for the modified nucleation rate in FLRW backgrounds, as well as a parameter study of bubble walls propagating into inhomogeneous (LTB) or FLRW spacetimes, both in the thin-wall approximation. We show that in our model, matter in the background often prevents bubbles from successful expansion and forces them to collapse. For cases where they do expand, we give arguments why the effects on the interior spacetime are small for a wide range of reasonable parameters and discuss the limitations of the employed approximations., Comment: 29 pages, 8 figures, typos corrected, matches published version

Colliding flows are a commonly used scenario for the formation of molecular clouds in numerical simulations. Due to the thermal instability of the warm neutral medium, turbulence is produced by cooling. We carry out a two-dimensional numerical study of such colliding flows in order to test whether statistical properties inferred from adaptive mesh refinement (AMR) simulations are robust with respect to the applied refinement criteria. We compare probability density functions of various quantities as well as the clump statistics and fractal dimension of the density fields in AMR simulations to a static-grid simulation. The static grid with 2048^2 cells matches the resolution of the most refined subgrids in the AMR simulations. The density statistics is reproduced fairly well by AMR. Refinement criteria based on the cooling time or the turbulence intensity appear to be superior to the standard technique of refinement by overdensity. Nevertheless, substantial differences in the flow structure become apparent. In general, it is difficult to separate numerical effects from genuine physical processes in AMR simulations., Comment: 6 pages, 6 figures, submitted to A&A

We performed numerical simulations of supersonic isothermal turbulence driven by mostly compressive large-scale forcing, using both a static grid and adaptive mesh refinement with an effective resolution N=768^3. After a transient phase dominated by shocks, turbulence evolves into a steady state with an RMS Mach number about 2.5, in which cloud-like structures of over-dense gas are surrounded by highly rarefied gas. The index of the turbulence energy spectrum function beta = 2.0 in the shock-dominated phase. As the flow approaches statistical equilibrium, the spectrum flattens, with beta = 1.9. For the scaling exponent of the root mean square velocity fluctuation, we obtain gamma = 0.43 from the velocity structure functions of second order. These results are well within the range of observed scaling properties for the velocity dispersion in molecular clouds. Calculating structure functions of order p=1,...,5, we find for all scaling exponents significant deviations from the Kolmogorov-Burgers model proposed by Boldyrev. Our results are very well described by a general log-Poisson model with a higher degree of intermittency, which implies an influence of the forcing on the scaling properties. Contrary to previous numerical results for isothermal turbulence, we obtain a skewed probability density function of the mass density fluctuations that is not consistent with log-normal statistics and entails a substantially higher fraction of mass in the density peaks than implied by the Padoan-Nordlund relation between the variance of the density fluctuations and the Mach number. In conclusion, it seems necessary to account for the production mechanism of turbulence in the ISM., Comment: 20 pages, 18 figures, submitted to A&A

The sensitivity of inflationary spectra to initial conditions is addressed in the context of a phenomenological model that breaks Lorentz invariance by dissipative effects above some threshold energy $\Lambda$. These effects are obtained dynamically by coupling the fluctuation modes to extra degrees of freedom which are unobservable below $\Lambda$. Because of the strong dissipative effects in the early propagation, only the state of the extra degrees of freedom is relevant for the power spectrum. If this state is the ground state, and if $\Lambda$ is much larger than the Hubble scale $H$, the standard spectrum is recovered. Using analytical and numerical methods, we calculate the modifications for a large class of dissipative models. For all of these, we show that the leading modification (in an expansion in $H/\Lambda$) is linear in the decay rate evaluated at horizon exit, and that high frequency superimposed oscillations are not generated. The modification is negative when the decay rate decreases slower than the cube of $H$, which means that there is a loss of power on the largest scales., Comment: 12 pages, 6 figures. Minor changes after refereeing and editing

We reconsider trans-Planckian corrections to inflationary spectra by taking into account a physical effect which has been overlooked and which could have important consequences. We assume that the short length scale characterizing the new physics is endowed with a finite width, the origin of which could be found in quantum gravity. As a result, the leading corrections responsible for superimposed osillations in the CMB temperature anisotropies are generically damped by the blurring of the UV scale. To determine the observational ramifications of this damping, we compare it to that which effectively occurs when computing the angular power spectrum of temperature anisotropies. The former gives an overall change of the oscillation amplitudes whereas the latter depends on the angular scale. Therefore, in principle they could be distinguished. In any case, the observation of superimposed oscillations would place tight constraint on the variance of the UV cutoff.

Inflation can act as a space-time microscope for Planck or string scale effects, leaving potentially observable traces in the primordial perturbation spectrum. I discuss two frameworks that were used recently to study this phenomenon: nonlinear dispersion and short distance uncertainty., Comment: 9 pages. Talk presented at COSMO-01, Rovaniemi, Finland, August 29 -- September 4, 2001

We investigate the quantum evolution of large black holes that nucleate spontaneously in de Sitter space. By numerical computation in the s-wave and one-loop approximations, we verify claims that such black holes can initially "anti-evaporate" instead of shrink. We show, however, that this is a transitory effect. It is followed by an evaporating phase, which we are able to trace until the black holes are small enough to be treated as Schwarzschild. Under generic perturbations, the nucleated geometry is shown to decay into a ring of de Sitter regions connected by evaporating black holes. This confirms that de Sitter space is globally unstable and fragments into disconnected daughter universes., Comment: 10 pages, 8 figures, to appear in PRD