Black holes of a billion solar masses are observed in the infant Universe a few hundred million years after the big bang. The direct collapse of protogalactic gas clouds in primordial haloes with ...T... ≥ 10... K provides the most promising way to assemble massive black holes. In this study, we aim to determine the characteristic mass scale of seed black holes and the time evolution of the accretion rates resulting from the direct collapse model. We explore the formation of supermassive black holes via cosmological large eddy simulations (LES) by employing sink particles and following their evolution for 20 000 yr after the formation of the first sink. As the resulting protostars were shown to have cool atmospheres in the presence of strong accretion, we assume here that UV feedback is negligible during this calculation. We confirm this result in a comparison run without sinks. Our findings show that black hole seeds with characteristic mass of 10... M... are formed in the presence of strong Lyman-Werner flux which leads to an isothermal collapse. The characteristic mass is about two times higher in LES compared to the implicit large eddy simulations. The accretion rates increase with time and reach a maximum value of 10 M... yr... after 104 yr. Our results show that the direct collapse model is clearly feasible as it provides the expected mass of the seed black holes. (ProQuest: ... denotes formulae/symbols omitted.)
We investigate the radiative transfer of Lyα 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 Lyα spectral line profiles; (ii) (typical) shell models do not allow for the escape of ionizing photons; and (iii) the observation and expectation that winds are more complex, anisotropic phenomena. We examine the influence of inclination on the Lyα 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.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
Context. Lyman-α emitters (LAEs) are a promising probe of the large-scale structure at high redshift, z ≳ 2. In particular, the Hobby-Eberly Telescope Dark Energy Experiment aims at observing LAEs at ...1.9 < z < 3.5 to measure the baryon acoustic oscillation (BAO) scale and the redshift-space distortion (RSD). However, it has been pointed out that the complicated radiative transfer (RT) of the resonant Lyman-α emission line generates an anisotropic selection bias in the LAE clustering on large scales, s ≳ 10 Mpc. This effect could potentially induce a systematic error in the BAO and RSD measurements. Also, there exists a recent claim to have observational evidence of the effect in the Lyman-α intensity map, albeit statistically insignificant. Aims. We aim at quantifying the impact of the Lyman-α RT on the large-scale galaxy clustering in detail. For this purpose, we study the correlations between the large-scale environment and the ratio of an apparent Lyman-α luminosity to an intrinsic one, which we call the “observed fraction”, at 2 < z < 6. Methods. We apply our Lyman-α RT code by post-processing the full Illustris simulations. We simply assume that the intrinsic luminosity of the Lyman-α emission is proportional to the star formation rate of galaxies in Illustris, yielding a sufficiently large sample of LAEs to measure the anisotropic selection bias. Results. We find little correlation between large-scale environment and the observed fraction induced by the RT, and hence a smaller anisotropic selection bias than has previously been claimed. We argue that the anisotropy was overestimated in previous work due to insufficient spatial resolution; it is important to keep the resolution such that it resolves the high-density region down to the scale of the interstellar medium, that is, ~1 physical kpc. We also find that the correlation can be further enhanced by assumptions in modeling intrinsic Lyman-α emission.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
Baryonic streaming motions produced prior to the epoch of recombination became supersonic during the cosmic dark ages. Various studies suggest that such streaming velocities change the halo ...statistics and also influence the formation of Population III stars. In this study, we aim to explore the impact of streaming velocities on the formation of supermassive black holes at z>10 via the direct collapse scenario. To accomplish this goal, we perform cosmological large eddy simulations for two haloes of a few times 107M⊙ with initial streaming velocities of 3, 6 and 9 km s−1. These massive primordial haloes illuminated by the strong Lyman–Werner flux are the potential cradles for the formation of direct collapse seed black holes. To study the evolution for longer times, we employ sink particles and track the accretion for 10 000 years. Our findings show that higher streaming velocities increase the circular velocities from about 14 to 16 km s−1. They also delay the collapse of haloes for a few million years, but do not have any significant impact on the halo properties such as turbulent energy, radial velocity, density and accretion rates. Sink particles of about ∼105M⊙ are formed at the end of our simulations and no clear distribution of sink masses is observed in the presence of streaming motions. It is further found that the impact of streaming velocities is less severe in massive haloes compared to the minihaloes as reported in the previous studies.
The development of turbulent gas flows in the intra-cluster medium and in the core of a galaxy cluster is studied by means of adaptive mesh refinement (AMR) cosmological simulations. A series of six ...runs was performed, employing identical simulation parameters but different criteria for triggering the mesh refinement. In particular, two different AMR strategies were followed, based on the regional variability of control variables of the flow and on the overdensity of subclumps, respectively. We show that both approaches, albeit with different results, are useful to get an improved resolution of the turbulent flow in the ICM. The vorticity is used as a diagnostic for turbulence, showing that the turbulent flow is not highly volume filling but has a large area-covering factor, in agreement with previous theoretical expectations. The measured turbulent velocity in the cluster core is larger than 200 km s−1, and the level of turbulent pressure contribution to the cluster hydrostatic equilibrium is increased by using the improved AMR criteria.
The injection and evolution of turbulence in the intergalactic medium is studied by means of mesh-based hydrodynamical simulations, including a subgrid-scale (SGS) model for small-scale unresolved ...turbulence. The simulations show that the production of turbulence has a different redshift dependence in the intracluster medium (ICM) and the warm-hot intergalactic medium (WHIM). We show that the turbulence in the ICM is produced chiefly by merger-induced shear flows, whereas the production in the WHIM is dominated by shock interactions. Secondly, the effect of dynamical pressure support on the gravitational contraction has been studied. This turbulent support is stronger in the WHIM gas at baryon overdensities 1 ≲δ≲ 100 and less relevant for the ICM. Although the relative mass fraction of the gas with large vorticity is considerable (52 per cent in the ICM), we find that for only about 10 per cent in mass this is dynamically relevant, namely not associated with an equally large thermal pressure support. According to this result, a significant non-thermal pressure support counteracting the gravitational contraction is a localized characteristic in the cosmic flow, rather than a widespread feature.
We investigate turbulence generated by cosmological structure formation by means of large eddy simulations using adaptive mesh refinement. In contrast to the widely used implicit large eddy ...simulations, which resolve a limited range of length-scales and treat the effect of turbulent velocity fluctuations below the grid scale solely by numerical dissipation, we apply a subgrid-scale model for the numerically unresolved fraction of the turbulence energy. For simulations with adaptive mesh refinement, we utilize a new methodology that allows us to adjust the scale-dependent energy variables in such a way that the sum of resolved and unresolved energies is globally conserved. We test our approach in simulations of randomly forced turbulence, a gravitationally bound cloud in a wind, and the Santa Barbara cluster. To treat inhomogeneous turbulence, we introduce an adaptive Kalman filtering technique that separates turbulent velocity fluctuations on resolved length-scales from the non-turbulent bulk flow. From the magnitude of the fluctuating component and the subgrid-scale turbulence energy, a total turbulent velocity dispersion of several 100 km s−1 is obtained for the Santa Barbara cluster, while the low-density gas outside the accretion shocks is nearly devoid of turbulence. The energy flux through the turbulent cascade and the dissipation rate predicted by the subgrid-scale model correspond to dynamical time-scales around 5 Gyr, independent of numerical resolution.
Hot and turbulent gas in clusters Schmidt, W; Engels, J. F; Niemeyer, J. C ...
Monthly notices of the Royal Astronomical Society,
06/2016, Volume:
459, Issue:
1
Journal Article
Peer reviewed
Open access
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 towards cluster cores. This enables us to infer the velocity dispersion associated with the turbulent fluctuation relative to the bulk flow. For haloes with masses above 1013 M⊙, 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.
Abstract
Major mergers are considered to be a significant source of turbulence in clusters. We performed a numerical simulation of a major merger event using nested-grid initial conditions, adaptive ...mesh refinement, radiative cooling of primordial gas and a homogeneous ultraviolet background. By calculating the microscopic viscosity on the basis of various theoretical assumptions and estimating the Kolmogorov length from the turbulent dissipation rate computed with a subgrid-scale model, we are able to demonstrate that most of the warm–hot intergalactic mediums can sustain a fully turbulent state only if the magnetic suppression of the viscosity is considerable. Accepting this as premise, it turns out that ratios of turbulent and thermal quantities change only little in the course of the merger. This confirms the tight correlations between the mean thermal and non-thermal energy content for large samples of clusters in earlier studies, which can be interpreted as second self-similarity on top of the self-similarity for different halo masses. Another long-standing question is how and to which extent turbulence contributes to the support of the gas against gravity. From a global perspective, the ratio of turbulent and thermal pressures is significant for the clusters in our simulation. On the other hand, a local measure is provided by the compression rate, i.e. the growth rate of the divergence of the flow. Particularly for the intracluster medium, we find that the dominant contribution against gravity comes from thermal pressure, while compressible turbulence effectively counteracts the support. For this reason, it appears to be too simplistic to consider turbulence merely as an effective enhancement of thermal energy.
We present a one-equation subgrid scale model that evolves the turbulence energy corresponding to unresolved velocity fluctuations in large eddy simulations. The model is derived in the context of ...the Germano consistent decomposition of the hydrodynamical equations. The eddy-viscosity closure for the rate of energy transfer from resolved toward subgrid scales is localised by means of a dynamical procedure for the computation of the closure parameter. Therefore, the subgrid scale model applies to arbitrary flow geometry and evolution. For the treatment of microscopic viscous dissipation a semi-statistical approach is used, and the gradient-diffusion hypothesis is adopted for turbulent transport. A priori tests of the localised eddy-viscosity closure and the gradient-diffusion closure are made by analysing data from direct numerical simulations. As an a posteriori testing case, the large eddy simulation of thermonuclear combustion in forced isotropic turbulence is discussed. We intend the formulation of the subgrid scale model in this paper as a basis for more advanced applications in numerical simulations of complex astrophysical phenomena involving turbulence.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK