Researcher profile

Dominik R. G. Schleicher

Dominik R. G. Schleicher contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
39works
0followers
12topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

39 published item(s)

preprint2025arXiv

Rapid formation of a very massive star >50000 $M_\odot$ and subsequently an IMBH from runaway collisions. Direct N-body and Monte Carlo simulations of dense star clusters

Context. We present simulations of a massive young star cluster using \textsc{Nbody6++GPU} and \textsc{MOCCA}. The cluster is initially more compact than previously published models, with one million stars, a total mass of $5.86 \times 10^5~\mathrm{M}_{\odot}$, and a half-mass radius of $0.1~\mathrm{pc}$. Aims. We analyse the formation and growth of a very massive star (VMS) through successive stellar collisions and investigate the subsequent formation of an intermediate-mass black hole (IMBH) in the core of a dense star cluster. Methods. We use both direct \textit{N}-body and Monte Carlo simulations, incorporating updated stellar evolution prescriptions (SSE/BSE) tailored to massive stars and VMSs. These include revised treatments of stellar radii, rejuvenation, and mass loss during collisions. While the prescriptions represent reasonable extrapolations into the VMS regime, the internal structure and thermal state of VMSs formed through stellar collisions remain uncertain, and future work may require further refinement. Results. We find that runaway stellar collisions in the cluster core produce a VMS exceeding $5 \times 10^4~\mathrm{M}_{\odot}$ within 5 Myr, which subsequently collapses into an IMBH. Conclusions. Our model suggests that dense stellar environments may enable the formation of very massive stars and massive black hole seeds through runaway stellar collisions. These results provide a potential pathway for early black hole growth in star clusters and offer theoretical context for interpreting recent JWST observations of young, compact clusters at high redshift.

preprint2022arXiv

Origin of eclipsing time variations: Contributions of different modes of the dynamo-generated magnetic field

The possibility to detect circumbinary planets and to study stellar magnetic fields through eclipsing time variations (ETVs) in binary stars has sparked an increase of interest in this area of research. We revisit the connection between stellar magnetic fields and the gravitational quadrupole moment $Q_{xx}$ and compare different dynamo-generated ETV models with our simulations. We present magnetohydrodynamical simulations of solar mass stars with rotation periods of 8.3, 1.2, and 0.8 days and perform a detailed analysis of the magnetic and quadrupole moment using spherical harmonic decomposition. The extrema of $Q_{xx}$ are associated with changes in the magnetic field structure. This is evident in the simulation with a rotation period of 1.2 days. Its magnetic field has a more complex behavior than in the other models, as the large-scale nonaxisymmetric field dominates throughout the simulation and the axisymmetric component is predominantly hemispheric. This triggers variations in the density field that follow the magnetic field asymmetry with respect to the equator, affecting the $zz$ component of the inertia tensor, and thus modulating $Q_{xx}$. The magnetic field of the two other runs are less variable in time and more symmetric with respect to the equator, such that the variations in the density are weaker, and therefore only small variations in $Q_{xx}$ are seen. If interpreted via the classical Applegate mechanism (tidal locking), the quadrupole moment variations obtained in the current simulations are about two orders of magnitude below those deduced from observations of post-common-envelope binaries. However, if no tidal locking is assumed, our results are compatible with the observed ETVs.

preprint2022arXiv

Simulations of fully convective M dwarfs: dynamo action with varying magnetic Prandtl numbers

M dwarfs are low-mass main-sequence stars, the most numerous type of stars in the solar neighbourhood, which are known to have significant magnetic activity. The aim of this work is to explore the dynamo solutions and magnetic fields of fully convective M dwarfs with varying magnetic Prandtl numbers ${\rm Pr_M}$, and a rotation period (${\rm P_{rot}}$) of 43 days. ${\rm Pr_M}$ is known to play an important role in the dynamo action; dynamos for low-$\rm Pr_M$ and large-$\rm Pr_M$ have very different properties. We performed three-dimensional magnetohydrodynamical (MHD) numerical simulations with the ``star-in-a-box'' model using stellar parameters for an M5 dwarf with $0.21 M_{\odot}$. We found that the dynamo solutions are sensitive to ${\rm Pr_M}$. The simulations at this rotation period present periodic cycles of the large-scale magnetic field up to ${\rm Pr_M} \leq 2$; for higher values the cycles disappear and irregular solutions to arise. Our results are consistent with previous studies and suggest that the dynamos operating in fully convective stars behave similarly as those in partially convective stars.

preprint2022arXiv

Synthetic observations using POLARIS: an application to simulations of massive prestellar cores

Young massive stars are usually found embedded in dense and massive molecular clumps which are known for being highly obscured and distant. During their formation process, the degree of deuteration can be used as a potential indicator of the very early formation stages. This is particularly effective when employing the abundance of H$_2$D$^+$. However, its low abundances and large distances make detections in massive sources hard to achieve. We present an application of the radiative transfer code POLARIS, with the goal to test the observability of the ortho-H$_2$D$^+$ transition $1_{10}$-$1_{11}$ (372.42 GHz) using simulations of high-mass collapsing cores that include deuteration chemistry. We analyzed an early and a late stage of the collapse of a 60 M$_{\odot}$ core, testing different source distances. For all cases, we generated synthetic single-dish and interferometric observations and studied the differences in both techniques. The column densities we derive are comparable to values reported for similar sources. These estimates depend on the extent over which they are averaged, and sources with compact emission they can be highly affected by beam dilution. Combined ALMA-ACA observations improve in signal-to-noise ratio and lead to better column density estimates as compared to ALMA alone. We confirm the feasibility to study ortho-H$_2$D$^+$ emission up to distances of 7 kpc. We provide a proof-of-concept of our framework for synthetic observations and highlight its importance when comparing numerical simulations with real observations. This work also proves how relevant it is to combine single-dish and interferometric measurements to derive appropriate source column densities.

preprint2021arXiv

Magnetic field amplification in accretion discs around the first stars: implications for the primordial IMF

Magnetic fields play an important role in the dynamics of present-day molecular clouds. Recent work has shown that magnetic fields are equally important for primordial clouds, which form the first stars in the Universe. While the primordial magnetic field strength on cosmic scales is largely unconstrained, theoretical models strongly suggest that a weak seed field existed in the early Universe. We study how the amplification of such a weak field can influence the evolution of accretion discs around first stars, and thus affect the primordial initial mass function (IMF). We perform a suite of 3D ideal magneto-hydrodynamic (MHD) simulations with different initial field strengths and numerical resolutions. We find that, in simulations with sufficient spatial resolution to resolve the Jeans scale during the collapse, even initially weak magnetic fields grow exponentially to become dynamically important due to both the so-called 'small-scale turbulent dynamo' and the 'large-scale mean-field dynamo'. Capturing the small-scale dynamo action depends primarily on how well we resolve the Jeans length, while capturing the large-scale dynamo depends on the Jeans resolution as well as the maximum absolute resolution. Provided enough resolution, we find that fragmentation does not depend strongly on the initial field strength, because even weak fields grow to become strong. However, fragmentation in runs with magnetic fields differs significantly from those without magnetic fields. We conclude that the development of dynamically strong magnetic fields during the formation of the first stars is likely inevitable, and that these fields had a significant impact on the primordial IMF.

preprint2020arXiv

A new proxy to estimate the cosmic-ray ionisation rate in dense cores

Cosmic rays are a global source of ionisation, and the ionisation fraction represents a fundamental parameter in the interstellar medium. Ions couple to magnetic fields, affect the chemistry, and the dynamics of star-forming regions as well as planetary atmospheres. However, the cosmic-ray ionisation rate represents one of the bottlenecks for astrochemical models, and its determination is one of the most puzzling problems in astrophysics. While for diffuse clouds reasonable values have been provided from H$_3^+$ observations, for dense clouds, due to the lack of rotational transitions, this is not possible, and estimates are strongly biased by the employed model. We present here an analytical expression, obtained from first principles, to estimate the cosmic-ray ionisation rate from observational quantities. The theoretical predictions are validated with high-resolution three-dimensional numerical simulations and applied to the well known core L1544; we obtained an estimate of $ζ_2 \sim 2-3 \times 10^{-17}$ s$^{-1}$. Our results and the analytical formulae provided represent the first model-independent, robust tool to probe the cosmic-ray ionisation rate in the densest part of star-forming regions (on spatial scales of $R \leq 0.05$ pc). An error analysis is presented to give statistical relevance to our study.

preprint2020arXiv

Synthetic observations of deuterated molecules in massive prestellar cores

Young massive stars are usually found embedded in dense massive molecular clumps and are known for being highly obscured and distant. During their formation process, deuteration is regarded as a potentially good indicator of the very early formation stages. In this work, we test the observability of the ground-state transition of ortho-H$_2$D$^+$ $J_{\rm {K_a, K_c}} = 1_{10}$-$1_{11} $ by performing interferometric and single-dish synthetic observations using magneto-hydrodynamic simulations of high-mass collapsing molecular cores, including deuteration chemistry. We studied different evolutionary times and source distances (from 1 to 7 kpc) to estimate the information loss when comparing the column densities inferred from the synthetic observations to the column densities in the model. We mimicked single-dish observations considering an APEX-like beam and interferometric observations using CASA and assuming the most compact configuration for the ALMA antennas. We found that, for centrally concentrated density distributions, the column densities are underestimated by about 51% in the case of high-resolution ALMA observations ($\leqslant$1") and up to 90% for APEX observations (17"). Interferometers retrieve values closer to the real ones, however, their finite spatial sampling results in the loss of contribution from large-scale structures due to the lack of short baselines. We conclude that, the emission of o-H$_2$D$^+$ in distant massive dense cores is faint and would require from $\sim$1 to $\sim$7 hours of observation at distances of 1 and 7 kpc, respectively, to achieve a 14$σ$ detection in the best case scenario. Additionally, the column densities derived from such observations will certainly be affected by beam dilution in the case of single-dishes and spatial filtering in the case of interferometers.

preprint2019arXiv

Magneto-hydrodynamical origin of eclipsing time variations in post-common-envelope binaries for solar mass secondaries

Eclipsing time variations have been observed for a wide range of binary systems, including post-common-envelope binaries. A frequently proposed explanation, apart from the possibility of having a third body, is the effect of magnetic activity, which may alter the internal structure of the secondary star, particularly its quadrupole moment, and thereby cause quasi-periodic oscillations. Here we present two compressible non-ideal magneto-hydrodynamical (MHD) simulations of the magnetic dynamo in a solar mass star, one of them with three times the solar rotation rate ("slow rotator"), the other one with twenty times the solar rotation rate ("rapid rotator"), to account for the high rotational velocities in close binary systems. For the slow rotator, we find that both the magnetic field and the stellar quadrupole moment change in a quasi-periodic manner, leading to O-C (observed - corrected times of the eclipse) variations of ~0.025 s. For the rapid rotator, the behavior of the magnetic field as well as the quadrupole moment changes become considerably more complex, due to the less coherent dynamo solution. The resulting O-C variations are of the order 0.13 s. The observed system V471~Tau shows two modes of eclipsing time variations, with amplitudes of 151 s and 20 s, respectively. However, the current simulations may not capture all relevant effects due to the neglect of the centrifugal force and self-gravity. Considering the model limitations and that the rotation of V471 Tau is still a factor of 2.5 faster than our rapid rotator, it may be conceivable to reach the observed magnitudes.

preprint2015arXiv

Modeling gravitational instabilities in self-gravitating protoplanetary disks with adaptive mesh refinement techniques

The astonishing diversity in the observed planetary population requires theoretical efforts and advances in planet formation theories. Numerical approaches provide a method to tackle the weaknesses of current planet formation models and are an important tool to close gaps in poorly constrained areas. We present a global disk setup to model the first stages of giant planet formation via gravitational instabilities (GI) in 3D with the block-structured adaptive mesh refinement (AMR) hydrodynamics code ENZO. With this setup, we explore the impact of AMR techniques on the fragmentation and clumping due to large-scale instabilities using different AMR configurations. Additionally, we seek to derive general resolution criteria for global simulations of self-gravitating disks of variable extent. We run a grid of simulations with varying AMR settings, including runs with a static grid for comparison, and study the effects of varying the disk radius. Adopting a marginally stable disk profile (Q_init=1), we validate the numerical robustness of our model for different spatial extensions, from compact to larger, extended disks (R_disk = 10, 100 and 300 AU, M_disk ~ 0.05 M_Sun, M_star = 0.646 M_Sun). By combining our findings from the resolution and parameter studies we find a lower limit of the resolution to be able to resolve GI induced fragmentation features and distinct, turbulence inducing clumps. Irrespective of the physical extension of the disk, topologically disconnected clump features are only resolved if the fragmentation-active zone of the disk is resolved with at least 100 cells, which holds as a minimum requirement for all global disk setups. Our simulations illustrate the capabilities of AMR-based modeling techniques for planet formation simulations and underline the importance of balanced refinement settings to reproduce fragmenting structures.

preprint2015arXiv

The chemical evolution of self-gravitating primordial disks

Numerical simulations show the formation of self-gravitating primordial disks during the assembly of the first structures in the Universe, in particular during the formation of Pop.~III and supermassive stars. Their subsequent evolution is expected to be crucial to determine the mass scale of the first cosmological objects, which depends on the temperature of the gas and the dominant cooling mechanism. Here, we derive a one-zone framework to explore the chemical evolution of such disks and show that viscous heating leads to the collisional dissociation of an initially molecular gas. The effect is relevant on scales of 10 AU (1000 AU) for a central mass of 10 M_sun (10^4 M_sun) at an accretion rate of 10^{-1} M_sun yr^{-1}, and provides a substantial heat input to stabilize the disk. If the gas is initially atomic, it remains atomic during the further evolution, and the effect of viscous heating is less significant. The additional thermal support is particularly relevant for the formation of very massive objects, such as the progenitors of the first supermassive black holes. The stabilizing impact of viscous heating thus alleviates the need for a strong radiation background as a means of keeping the gas atomic.

preprint2014arXiv

Low-metallicity star formation: Relative impact of metals and magnetic fields

Low-metallicity star formation poses a central problem of cosmology, as it determines the characteristic mass scale and distribution for the first and second generations of stars forming in our Universe. Here, we present a comprehensive investigation assessing the relative impact of metals and magnetic fields, which may both be present during low-metallicity star formation. We show that the presence of magnetic fields generated via the small-scale dynamo stabilises the protostellar disc and provides some degree of support against fragmentation. In the absence of magnetic fields, the fragmentation timescale in our model decreases by a factor of ~10 at the transition from Z=0 to Z>0, with subsequently only a weak dependence on metallicity. Similarly, the accretion timescale of the cluster is set by the large-scale dynamics rather than the local thermodynamics. In the presence of magnetic fields, the primordial disc can become completely stable, therefore forming only one central fragment. At Z>0, the number of fragments is somewhat reduced in the presence of magnetic fields, though the shape of the mass spectrum is not strongly affected in the limits of the statistical uncertainties. The fragmentation timescale, however, increases by roughly a factor of 3 in the presence of magnetic fields. Indeed, our results indicate comparable fragmentation timescales in primordial runs without magnetic fields and Z>0 runs with magnetic fields.

preprint2014arXiv

Planet formation from the ejecta of common envelopes

The close binary system NN Serpentis must have gone through a common envelope phase before the formation of its white dwarf. During this phase, a substantial amount of mass was lost from the envelope. The recently detected orbits of circumbinary planets are likely inconsistent with planet formation before the mass loss.We explore whether new planets may have formed from the ejecta of the common envelope and derive the expected planetary mass as a function of radius.We employed the Kashi & Soker model to estimate the amount of mass that is retained during the ejection event and inferred the properties of the resulting disk from the conservation of mass and angular momentum. The resulting planetary masses were estimated from models with and without radiative feedback. We show that the observed planetary masses can be reproduced for appropriate model parameters. Photoheating can stabilize the disks in the interior, potentially explaining the observed planetary orbits on scales of a few AU. We compare the expected mass scale of planets for 11 additional systems with observational results and find hints of two populations, one consistent with planet formation from the ejecta of common envelopes and the other a separate population that may have formed earlier. The formation of the observed planets from the ejecta of common envelopes seems feasible. The model proposed here can be tested through refined observations of additional post-common envelope systems. While it appears observationally challenging to distinguish between the accretion on pre-existing planets and their growth from new fragments, it may be possible to further constrain the properties of the protoplanetary disk through additional observations of current planetary candidates and post-common envelope binary systems.

preprint2014arXiv

The Turbulent Dynamo in Highly Compressible Supersonic Plasmas

The turbulent dynamo may explain the origin of cosmic magnetism. While the exponential amplification of magnetic fields has been studied for incompressible gases, little is known about dynamo action in highly-compressible, supersonic plasmas, such as the interstellar medium of galaxies and the early Universe. Here we perform the first quantitative comparison of theoretical models of the dynamo growth rate and saturation level with three-dimensional magnetohydrodynamical simulations of supersonic turbulence with grid resolutions of up to 1024^3 cells. We obtain numerical convergence and find that dynamo action occurs for both low and high magnetic Prandtl numbers Pm = nu/eta = 0.1-10 (the ratio of viscous to magnetic dissipation), which had so far only been seen for Pm >= 1 in supersonic turbulence. We measure the critical magnetic Reynolds number, Rm_crit = 129 (+43, -31), showing that the compressible dynamo is almost as efficient as in incompressible gas. Considering the physical conditions of the present and early Universe, we conclude that magnetic fields need to be taken into account during structure formation from the early to the present cosmic ages, because they suppress gas fragmentation and drive powerful jets and outflows, both greatly affecting the initial mass function of stars.

preprint2014arXiv

X-Ray Emission from Star-Forming Galaxies - Signatures of Cosmic Rays and Magnetic Fields

The evolution of magnetic fields in galaxies is still an open problem in astrophysics. In nearby galaxies the far-infrared-radio correlation indicates the coupling between magnetic fields and star formation. The correlation arises from the synchrotron emission of cosmic ray electrons traveling through the interstellar magnetic fields. However, with an increase of the interstellar radiation field (ISRF), inverse Compton scattering becomes the dominant energy loss mechanism of cosmic ray electrons with a typical emission frequency in the X-ray regime. The ISRF depends on the one hand on the star formation rate and becomes stronger in starburst galaxies, and on the other hand increases with redshift due to the evolution of the cosmic microwave background. With a model for the star formation rate of galaxies, the ISRF, and the cosmic ray spectrum, we can calculate the expected X-ray luminosity resulting from the inverse Compton emission. Except for galaxies with an active galactic nucleus the main additional contribution to the X-ray luminosity comes from X-ray binaries. We estimate this contribution with an analytical model as well as with an observational relation, and compare it to the pure inverse Compton luminosity. Using data from the Chandra Deep Field Survey and far-infrared observations from ALMA we then determine upper limits for the cosmic ray energy. Assuming that the magnetic energy in a galaxy is in equipartition with the energy density of the cosmic rays, we obtain upper limits for the magnetic field strength. Our results suggest that the mean magnetic energy of young galaxies is similar to the one in local galaxies. This points toward an early generation of galactic magnetic fields, which is in agreement with current dynamo evolution models.

preprint2013arXiv

A new interpretation of the far-infrared - radio correlation and the expected breakdown at high redshift

(Abrigded) Observations of galaxies up to z 2 show a tight correlation between far-infrared and radio continuum emission. We explain the far-infrared - radio continuum correlation by relating star formation and magnetic field strength in terms of turbulent magnetic field amplification, where turbulence is injected by supernova explosions from massive stars. We calculate the expected amount of turbulence in galaxies based on their star formation rates, and infer the expected magnetic field strength due to turbulent dynamo amplification. We estimate the timescales for cosmic ray energy losses via synchrotron emission, inverse Compton scattering, ionization and bremsstrahlung emission, probing up to which redshift strong synchrotron emission can be maintained. We find that the correlation between star formation rate and magnetic field strength in the local Universe can be understood as a result of turbulent magnetic field amplification. If the typical gas density in the interstellar medium increases at high z, we expect an increase of the magnetic field strength and the radio emission, as indicated by current observations. Such an increase would imply a modification of the far-infrared - radio correlation. We expect a breakdown when inverse Compton losses start dominating over synchrotron emission. For a given star formation surface density, we calculate the redshift where the breakdown occurs, yielding z (Sigma_SFR/0.0045 M_solar kpc^{-2} yr^{-1})^{1/(6-alpha/2)}. In this relation, the parameter αdescribes the evolution of the characteristic ISM density in galaxies as (1+z)^α. Both the possible raise of the radio emission at high redshift and the final breakdown of the far-infrared -- radio correlation at a critical redshift will be probed by the Square Kilometre Array (SKA) and its pathfinders, while the typical ISM density in galaxies will be probed with ALMA.

preprint2013arXiv

Magnetic Field Amplification in Young Galaxies

The Universe at present is highly magnetized, with fields of the order of a few 10^-5 G and coherence lengths larger than 10 kpc in typical galaxies like the Milky Way. We propose that the magnetic field was amplified to this values already during the formation and the early evolution of the galaxies. Turbulence in young galaxies is driven by accretion as well as by supernova (SN) explosions of the first generation of stars. The small-scale dynamo can convert the turbulent kinetic energy into magnetic energy and amplify very weak primordial magnetic seed fields on short timescales. The amplification takes place in two phases: in the kinematic phase the magnetic field grows exponentially, with the largest growth on the smallest non-resistive scale. In the following non-linear phase the magnetic energy is shifted towards larger scales until the dynamo saturates on the turbulent forcing scale. To describe the amplification of the magnetic field quantitatively we model the microphysics in the interstellar medium (ISM) of young galaxies and determine the growth rate of the small-scale dynamo. We estimate the resulting saturation field strengths and dynamo timescales for two turbulent forcing mechanisms: accretion-driven turbulence and SN-driven turbulence. We compare them to the field strength that is reached, when only stellar magnetic fields are distributed by SN explosions. We find that the small-scale dynamo is much more efficient in magnetizing the ISM of young galaxies. In the case of accretion-driven turbulence a magnetic field strength of the order of 10^-6 G is reached after a time of 24-270 Myr, while in SN-driven turbulence the dynamo saturates at field strengths of typically 10^-5 G after only 4-15 Myr. This is considerably shorter than the Hubble time. Our work can help to understand why present-day galaxies are highly magnetized.

preprint2013arXiv

Massive black hole factories: Supermassive and quasi-star formation in primordial halos

Supermassive stars and quasi-stars (massive stars with a central black hole) are both considered as potential progenitors for the formation of supermassive black holes. They are expected to form from rapidly accreting protostars in massive primordial halos. We explore how long rapidly accreting protostars remain on the Hayashi track, implying large protostellar radii and weak accretion luminosity feedback. We assess the potential role of energy production in the nuclear core, and determine what regulates the evolution of such protostars into quasi-stars or supermassive stars. We follow the contraction of characteristic mass scales in rapidly accreting protostars, and infer the timescales for them to reach nuclear densities. We compare the characteristic timescales for nuclear burning with those for which the extended protostellar envelope can be maintained. We find that the extended envelope can be maintained up to protostellar masses of 3.6x10^8 \dot{m}^3 solar, where \dot{m} denotes the accretion rate in solar masses per year. We expect the nuclear core to exhaust its hydrogen content in 7x10^6 yrs. If accretion rates \dot{m}>>0.14 can still be maintained at this point, a black hole may form within the accreting envelope, leading to a quasi-star. Alternatively, the accreting object will gravitationally contract to become a main-sequence supermassive star. Due to the limited gas reservoir in dark matter halos with 10^7 solar masses, the accretion rate onto the central object may drop at late times, implying the formation of supermassive stars as the typical outcome of direct collapse. However, if high accretion rates are maintained, a quasi-star with an interior black hole may form.

preprint2013arXiv

Star formation and accretion in the circumnuclear disks of active galaxies

We explore the evolution of supermassive black holes (SMBH) centered in a circumnuclear disk (CND) as a function of the mass supply from the host galaxy and considering different star formation laws, which may give rise to a self-regulation via the injection of supernova-driven turbulence. A system of equations describing star formation, black hole accretion and angular momentum transport was solved for an axisymmetric disk in which the gravitational potential includes contributions from the black hole, the disk and the hosting galaxy. Our model extends the framework provided by Kawakatu et al. (2008) by separately considering the inner and outer part of the disk, and by introducing a potentially non-linear dependence of the star formation rate on the gas surface density and the turbulent velocity. The star formation recipes are calibrated using observational data for NGC 1097, while the accretion model is based on turbulent viscosity as a source of angular momentum transport in a thin viscous accretion disk. We find that current data provide no strong constraint on the star formation recipe, and can in particular not distinguish between models entirely regulated by the surface density, and models including a dependence on the turbulent velocity. The evolution of the black hole mass, on the other hand, strongly depends on the applied star formation law, as well as the mass supply from the host galaxy. We suggest to explore the star formation process in local AGN with high-resolution ALMA observations to break the degeneracy between different star formation models.

preprint2013arXiv

The small-scale dynamo: Breaking universality at high Mach numbers

(Abridged) The small-scale dynamo may play a substantial role in magnetizing the Universe under a large range of conditions, including subsonic turbulence at low Mach numbers, highly supersonic turbulence at high Mach numbers and a large range of magnetic Prandtl numbers Pm, i.e. the ratio of kinetic viscosity to magnetic resistivity. Low Mach numbers may in particular lead to the well-known, incompressible Kolmogorov turbulence, while for high Mach numbers, we are in the highly compressible regime, thus close to Burgers turbulence. In this study, we explore whether in this large range of conditions, a universal behavior can be expected. Our starting point are previous investigations in the kinematic regime. Here, analytic studies based on the Kazantsev model have shown that the behavior of the dynamo depends significantly on Pm and the type of turbulence, and numerical simulations indicate a strong dependence of the growth rate on the Mach number of the flow. Once the magnetic field saturates on the current amplification scale, backreactions occur and the growth is shifted to the next-larger scale. We employ a Fokker-Planck model to calculate the magnetic field amplification during the non-linear regime, and find a resulting power-law growth that depends on the type of turbulence invoked. For Kolmogorov turbulence, we confirm previous results suggesting a linear growth of magnetic energy. For more general turbulent spectra, where the turbulent velocity v_t scales with the characteristic length scale as u_\ell\propto \ell^{\vartheta}, we find that the magnetic energy grows as (t/T_{ed})^{2\vartheta/(1-\vartheta)}, with t the time-coordinate and T_{ed} the eddy-turnover time on the forcing scale of turbulence. For Burgers turbulence, \vartheta=1/2, a quadratic rather than linear growth may thus be expected, and a larger timescale until saturation is reached.

preprint2012arXiv

Generation of strong magnetic fields via the small-scale dynamo during the formation of the first stars

Here we summarize our recent results of high-resolution computer simulations on the turbulent amplification of weak magnetic seed fields showing that such fields will be exponentially amplified also during the gravitational collapse reminiscent to the situation during primordial star formation. The exponential magnetic field amplification is driven by the turbulent small-scale dynamo that can be only observed in computer simulations if the turbulent motions in the central core are sufficiently resolved. We find that the Jeans length, which determines the central core region, has to be resolved by at least 30 grid cells to capture the dynamo activity. We conclude from our studies that strong magnetic fields will be unavoidably created already during the formation of the first stars in the Universe, potentially influencing their evolution and mass distribution.

preprint2012arXiv

Magnetic field amplification during gravitational collapse - Influence of initial conditions on dynamo evolution and saturation

We study the influence of initial conditions on the magnetic field amplification during the collapse of a magnetised gas cloud. We focus on the dependence of the growth and saturation level of the dynamo generated field on the turbulent properties of the collapsing cloud. In particular, we explore the effect of varying the initial strength and injection scale of turbulence and the initial uniform rotation of the collapsing magnetised cloud. In order to follow the evolution of the magnetic field in both the kinematic and the nonlinear regime, we choose an initial field strength of $\simeq 1\,\mkG$ with the magnetic to kinetic energy ratio, $E_{\rm m}/E_{\rm k} \sim 10^{-4}$. Both gravitational compression and the small-scale dynamo initially amplify the magnetic field. Further into the evolution, the dynamo-generated magnetic field saturates but the total magnetic field continues to grow because of compression. The saturation of the small-scale dynamo is marked by a change in the slope of $B/ρ^{2/3}$ and by a shift in the peak of the magnetic energy spectrum from small scales to larger scales. For the range of initial Mach numbers explored in this study, the dynamo growth rate increases as the Mach number increases from $v_{\rm rms}/c_{\rm s}\sim 0.2$ to 0.4 and then starts decreasing from $v_{\rm rms}/c_{\rm s}\sim 1.0$. We obtain saturation values of $E_{\rm m}/E_{\rm k} = 0.2 - 0.3$ for these runs. Simulations with different initial injection scales of turbulence also show saturation at similar levels. For runs with different initial rotation of the cloud, the magnetic energy saturates at $E_{\rm m}/E_{\rm k}\sim 0.2 - 0.4$ of the equipartition value. (Abridged)

preprint2012arXiv

Magnetic fields during high redshift structure formation

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}.

preprint2012arXiv

Small-Scale Dynamo Action in Primordial Halos

The first galaxies form due to gravitational collapse of primordial halos. During this collapse, weak magnetic seed fields get amplified exponentially by the small-scale dynamo - a process converting kinetic energy from turbulence into magnetic energy. We use the Kazantsev theory, which describes the small-scale dynamo analytically, to study magnetic field amplification for different turbulent velocity correlation functions. For incompressible turbulence (Kolmogorov turbulence), we find that the growth rate is proportional to the square root of the hydrodynamic Reynolds number, Re^(1/2). In the case of highly compressible turbulence (Burgers turbulence) the growth rate increases proportional to Re^(1/3). With a detailed chemical network we are able to follow the chemical evolution and determine the kinetic and magnetic viscosities (due to Ohmic and ambipolar diffusion) during the collapse of the halo. This way, we can calculate the growth rate of the small-scale dynamo quantitatively and predict the evolution of the small-scale magnetic field. As the magnetic energy is transported to larger scales on the local eddy-timescale, we obtain an estimate for the magnetic field on the Jeans scale. Even there, we find that equipartition with the kinetic energy is reached on small timescales. Dynamically relevant field structures can thus be expected already during the formation of the first objects in the Universe.

preprint2012arXiv

The impact of thermodynamics on gravitational collapse: filament formation and magnetic field amplification

Stars form by the gravitational collapse of interstellar gas. The thermodynamic response of the gas can be characterized by an effective equation of state. It determines how gas heats up or cools as it gets compressed, and hence plays a key role in regulating the process of stellar birth on virtually all scales, ranging from individual star clusters up to the galaxy as a whole. We present a systematic study of the impact of thermodynamics on gravitational collapse in the context of high-redshift star formation, but argue that our findings are also relevant for present-day star formation in molecular clouds. We consider a polytropic equation of state, P = k rho^Gamma, with both sub-isothermal exponents Gamma < 1 and super-isothermal exponents Gamma > 1. We find significant differences between these two cases. For Gamma > 1, pressure gradients slow down the contraction and lead to the formation of a virialized, turbulent core. Weak magnetic fields are strongly tangled and efficiently amplified via the small-scale turbulent dynamo on timescales corresponding to the eddy-turnover time at the viscous scale. For Gamma < 1, on the other hand, pressure support is not sufficient for the formation of such a core. Gravitational contraction proceeds much more rapidly and the flow develops very strong shocks, creating a network of intersecting sheets and extended filaments. The resulting magnetic field lines are very coherent and exhibit a considerable degree of order. Nevertheless, even under these conditions we still find exponential growth of the magnetic energy density in the kinematic regime.

preprint2012arXiv

Turbulent magnetic field amplification from the smallest to the largest magnetic Prandtl numbers

The small-scale dynamo provides a highly efficient mechanism for the conversion of turbulent into magnetic energy. In astrophysical environments, such turbulence often occurs at high Mach numbers, implying steep slopes in the turbulent spectra. It is thus a central question whether the small-scale dynamo can amplify magnetic fields in the interstellar or intergalactic media, where such Mach numbers occur. To address this long-standing issue, we employ the Kazantsev model for turbulent magnetic field amplification, systematically exploring the effect of different turbulent slopes, as expected for Kolmogorov, Burgers, the Larson laws and results derived from numerical simulations. With the framework employed here, we give the first solution encompassing the complete range of magnetic Prandtl numbers, including Pm << 1, Pm ~ 1 and Pm >> 1. We derive scaling laws of the growth rate as a function of hydrodynamic and magnetic Reynolds number for Pm << 1 and Pm >> 1 for all types of turbulence. A central result concerns the regime of Pm ~ 1, where the magnetic field amplification rate increases rapidly as a function of Pm. This phenomenon occurs for all types of turbulence we explored. We further find that the dynamo growth rate can be decreased by a few orders of magnitude for turbulence spectra steeper than Kolmogorov. We calculate the critical magnetic Reynolds number Rm_c for magnetic field amplification, which is highest for the Burgers case. As expected, our calculation shows a linear behaviour of the amplification rate close to the threshold proportional to Rm-Rm_c. Based on the Kazantsev model, we therefore expect the existence of the small-scale dynamo for any given value of Pm, as long as the magnetic Reynolds number is above the critical threshold.

preprint2012arXiv

WIMP DM and first stars: suppression of fragmentation in primordial star formation

We present the first 3D simulations to include the effects of dark matter annihilation feedback during the collapse of primordial mini-halos. We begin our simulations from cosmological initial conditions and account for dark matter annihilation in our treatment of the chemical and thermal evolution of the gas. The dark matter is modelled using an analytical density profile that responds to changes in the peak gas density. We find that the gas can collapse to high densities despite the additional energy input from the dark matter. No objects supported purely by dark matter annihilation heating are formed in our simulations. However, we find that the dark matter annihilation heating has a large effect on the evolution of the gas following the formation of the first protostar. Previous simulations without dark matter annihilation found that protostellar discs around Population III stars rapidly fragmented, forming multiple protostars that underwent mergers or ejections. When dark matter annihilation is included, however, these discs become stable to radii of 1000 AU or more. In the cases where fragmentation does occur, it is a wide binary that is formed.

preprint2011arXiv

A new Jeans resolution criterion for (M)HD simulations of self-gravitating gas: Application to magnetic field amplification by gravity-driven turbulence

Cosmic structure formation is characterized by the complex interplay between gravity, turbulence, and magnetic fields. The processes by which gravitational energy is converted into turbulent and magnetic energies, however, remain poorly understood. Here, we show with high-resolution, adaptive-mesh simulations that MHD turbulence is efficiently driven by extracting energy from the gravitational potential during the collapse of a dense gas cloud. Compressible motions generated during the contraction are converted into solenoidal, turbulent motions, leading to a natural energy ratio of E_sol/E_tot of approximately 2/3. We find that the energy injection scale of gravity-driven turbulence is close to the local Jeans scale. If small seeds of the magnetic field are present, they are amplified exponentially fast via the small-scale dynamo process. The magnetic field grows most efficiently on the smallest scales, for which the stretching, twisting, and folding of field lines, and the turbulent vortices are sufficiently resolved. We find that this scale corresponds to about 30 grid cells in the simulations. We thus suggest a new minimum resolution criterion of 30 cells per Jeans length in (magneto)hydrodynamical simulations of self-gravitating gas, in order to resolve turbulence on the Jeans scale, and to capture minimum dynamo amplification of the magnetic field. Due to numerical diffusion, however, any existing simulation today can at best provide lower limits on the physical growth rates. We conclude that a small, initial magnetic field can grow to dynamically important strength on time scales significantly shorter than the free-fall time of the cloud.

preprint2011arXiv

Lyman alpha emission from the first galaxies: Implications of UV backgrounds and the formation of molecules

The Lyman alpha line is a robust tracer of high redshift galaxies. We present estimates of Lyman alpha emission from a protogalactic halo illuminated by UV background radiation fields with various intensities. For this purpose, we performed cosmological hydrodynamics simulations with the adaptive mesh refinement code FLASH, including a detailed network for primordial chemistry,comprising the formation of primordial molecules, a multi-level model for the hydrogen atom as well as the photo-ionization and photo-dissociation processes in a UV background. We find that the presence of a background radiation field J_21 excites the emission of Lyman alpha photons, increasing the Lyman alpha luminosity up to two orders of magnitude. For a halo of \sim 10^10 M_sun, we find that a maximum flux of 5 \times 10^-15 erg cm^-2 s^-1 is obtained for J21 \times f_esc = 0.1, where f_esc is the escape fraction of the ionizing radiation. Depending on the environmental conditions, the flux may vary by three orders of magnitude. For J_21 \times f_esc > 0.1 the Lyman alpha luminosity decreases as the atomic hydrogen abundance becomes rather small. The fluxes derived here can be probed using Subaru and the upcoming James Webb Space Telescope. The emission of Lyman alpha photons is extended and comes from the envelope of the halo rather than its core. In the center of the halo, line trapping becomes effective above columns of 10^22 cm^-2 and suppresses the emission of Lyman alpha. In addition, cooling by primordial molecules may decrease the gas temperature in the central region, which further reduces Lyman alpha emission. In the central core, H_2 is photo-dissociated for a background flux of J_21 \geq 1000. For weaker radiation fields, i.e. J_21 < 0.1, H_2 and HD cooling are particularly strong in the center of the halo, leading to gas temperatures as low as \sim 100 K.

preprint2011arXiv

Mach Number Dependence of Turbulent Magnetic Field Amplification: Solenoidal versus Compressive Flows

We study the growth rate and saturation level of the turbulent dynamo in magnetohydrodynamical simulations of turbulence, driven with solenoidal (divergence-free) or compressive (curl-free) forcing. For models with Mach numbers ranging from 0.02 to 20, we find significantly different magnetic field geometries, amplification rates, and saturation levels, decreasing strongly at the transition from subsonic to supersonic flows, due to the development of shocks. Both extreme types of turbulent forcing drive the dynamo, but solenoidal forcing is more efficient, because it produces more vorticity.

preprint2011arXiv

Magnetic fields during primordial star formation

Recent FERMI observations provide a lower limit of 10^{-15} G for the magnetic field strength in the intergalactic medium (IGM). This is consistent with theoretical expectations based on the Biermann battery effect, which predicts such IGM fields already at redshifts z~10. During gravitational collapse, such magnetic fields can be amplified by compression and by turbulence, giving rise to the small-scale dynamo. On scales below the Jeans length, the eddy turnover timescale is much shorter than the free-fall timescale, so that saturation can be reached during collapse. This scenario has been tested and confirmed with magneto-hydrodynamical simulations following the collapse of a turbulent, weakly magnetized cloud. Based on a spectral analysis, we confirm that turbulence is injected on the Jeans scale. For the power spectrum of the magnetic field, we obtain the Kazantsev slope which is characteristic for the small-scale dynamo. A calculation of the critical length scales for ambipolar diffusion and Ohmic dissipation shows that these scales are always small enough to allow significant amplification of the magnetic field by small-scale eddies. We discuss potential implications for the protostellar accretion disk, with particular focus on the magneto-rotational instability, which may change the morphology of the disk and reduce the accretion rate by a factor of a few.

preprint2011arXiv

Magnetic fields in the Large-Scale Structure of the Universe

Magnetic fields appear to be ubiquitous in astrophysical environments. Their existence in the intracluster medium is established through observations of synchrotron emission and Faraday rotation. On the other hand, the nature of magnetic fields outside of clusters, where observations are scarce and controversial, remains largely unknown. In this chapter, we review recent developments in our understanding of the nature and origin of intergalactic magnetic fields, and in particular, intercluster fields. A plausible scenario for the origin of galactic and intergalactic magnetic fields is for seed fields, created in the early universe, to be amplified by turbulent flows induced during the formation of the large scale structure. We present several mechanisms for the generation of seed fields both before and after recombination. We then discuss the evolution and role of magnetic fields during the formation of the first starts. We describe the turbulent amplification of seed fields during the formation of large scale structure and the nature of the magnetic fields that arise. Finally, we discuss implications of intergalactic magnetic fields.

preprint2011arXiv

Primordial magnetic field constraints from the end of reionization

Primordial magnetic fields generated in the early universe are subject of considerable investigation, and observational limits on their strength are required to constrain the theory. Due to their impact on the reionization process, the strength of primordial fields can be limited using the latest data on reionization and the observed UV-luminosity function of high-redshift galaxies. Given the steep faint-end slope of the luminosity function, faint galaxies contribute substantial ionizing photons, and the low-luminosity cutoff has an impact on the total budget thereof. Magnetic pressure from primordial fields affects such cutoff by preventing collapse in halos with mass below 10^{10} M_solar (B_0 / 3 nG)^3, with B_0 the co-moving field strength. In this letter, the implications of these effects are consistently incorporated in a simplified model for reionization, and the uncertainties due to the cosmological parameters, the reionization parameters and the observed UV luminosity function are addressed. We show that the observed ionization degree at z\sim7 leads to the strongest upper limit of B_0\lsim 2-3nG. Stronger limits could follow from measurements of high ionization degree at z>7.

preprint2010arXiv

Black hole formation in primordial galaxies: chemical and radiative conditions

In massive primordial galaxies, the gas may directly collapse and form a single central massive object if cooling is suppressed. Line cooling by molecular hydrogen can be suppressed in the presence of a strong soft-ultraviolet radiation field, but the role played by other cooling mechanisms is less clear. In optically thin gas, Lyman-Alpha cooling can be very effective, maintaining the gas temperature below 10^4 K over many orders of magnitude in density. However, the large neutral hydrogen column densities present in primordial galaxies render them highly optically thick to Lyman-Alpha photons. In this letter, we examine in detail the effects of the trapping of these Lyman-Alpha photons on the thermal and chemical evolution of the gas. We show that despite the high optical depth in the Lyman series lines, cooling is not strongly suppressed, and proceeds via other atomic hydrogen transitions, in particular the 2s-1s and the 3-2 transitions. At densities larger than 10^9 cm^{-3}, collisional dissociation of molecular hydrogen becomes the dominant cooling process and decreases the gas temperature to about 5000 K. The gas temperature evolves with density as $T \propto ρ^{γ_{\rm eff} - 1}$, with $γ_{\rm eff} = 0.97-0.98$. The evolution is thus very close to isothermal, and so fragmentation is possible, but unlikely to occur during the initial collapse. However, after the formation of a massive central object, we expect that later-infalling, higher angular momentum material will form an accretion disk that may be unstable to fragmentation, which may give rise to star formation with a top-heavy IMF.

preprint2010arXiv

Probing high-redshift quasars with ALMA. I. Expected observables and potential number of sources

(abridged) We explore how ALMA observations can probe high-redshift galaxies in unprecedented detail. We discuss the main observables that are excited by the large-scale starburst, and formulate expectations for the chemistry and the fluxes in the center of active galaxies, where chemistry may be driven by the absorption of X-ray photons. We show that such X-ray dominated regions (XDRs) should be large enough to be resolved with ALMA, and predict the expected amount of emission in CO and various fine-structure lines. We discuss how such XDRs can be distinguished from a strong starburst on the same spatial scales based on the CO line SED. Our models are compared to known sources like NGC 1068 and APM 08279. We also analyze the properties of the z=6.42 quasar SDSS J114816.64+525150.3, and find that the observed emission in CO, [CII] and [CI] requires a dense warm and a low-density cold gas component. We estimate the expected number of sources at redshifts higher than 6, finding that one could expect one black hole with $10^6$ solar masses per arcmin$^2$.

preprint2010arXiv

Small-scale dynamo action during the formation of the first stars and galaxies. I. The ideal MHD limit

We explore the amplification of magnetic seed fields during the formation of the first stars and galaxies. During gravitational collapse, turbulence is created from accretion shocks, which may act to amplify weak magnetic fields in the protostellar cloud. Numerical simulations showed that such turbulence is sub-sonic in the first star-forming minihalos, and highly supersonic in the first galaxies with virial temperatures larger than 10^4 K. We investigate the magnetic field amplification during the collapse both for Kolmogorov and Burgers-type turbulence with a semi-analytic model that incorporates the effects of gravitational compression and small-scale dynamo amplification. We find that the magnetic field may be substantially amplified before the formation of a disk. On scales of 1/10 of the Jeans length, saturation occurs after ~10^8 yr. Although the saturation behaviour of the small-scale dynamo is still somewhat uncertain, we expect a saturation field strength of the order ~10^{-7} n^{0.5} G in the first star-forming halos, with n the number density in cgs units. In the first galaxies with higher turbulent velocities, the magnetic field strength may be increased by an order of magnitude, and saturation may occur after 10^6 to 10^7 yr. In the Kolmogorov case, the magnetic field strength on the integral scale (i.e. the scale with most magnetic power) is higher due to the characteristic power-law indices, but the difference is less than a factor of 2 in the saturated phase. Our results thus indicate that the precise scaling of the turbulent velocity with length scale is of minor importance. They further imply that magnetic fields will be significantly enhanced before the formation of a protostellar disk, where they may change the fragmentation properties of the gas and the accretion rate.

preprint2010arXiv

The formation of supermassive black holes in the first galaxies

We discuss the formation of supermassive black holes in the early universe, and how to probe their subsequent evolution with the upcoming mm/sub-mm telescope ALMA. We first focus on the chemical and radiative conditions for black hole formation, in particular considering radiation trapping and molecular dissociation effects. We then turn our attention towards the magnetic properties in the halos where the first black holes form, and show that the presence of turbulence may lead to a magnetic dynamo, which could support the black hole formation process by providing an efficient means of transporting the angular momentum. We finally focus on observable properties of high-redshift black holes with respect to ALMA, and discuss how to distinguish between chemistry driven by the starburst and chemistry driven by X-rays from the black hole.

preprint2010arXiv

The generation of strong magnetic fields during the formation of the first stars

Cosmological hydrodynamical simulations of primordial star formation suggest that the gas within the first star-forming halos is turbulent. This has strong implications on the subsequent evolution, in particular on the generation of magnetic fields. Using high-resolution numerical simulations, we show that in the presence of turbulence, weak seed magnetic fields are exponentially amplified by the small-scale dynamo during the formation of the first stars. We conclude that strong magnetic fields are generated during the birth of the first stars in the universe, potentially modifying the mass distribution of these stars and influencing the subsequent cosmic evolution. We find that the presence of the small-scale turbulent dynamo can only be identified in numerical simulations in which the turbulent motions in the central core are resolved with at least 32 grid cells.

preprint2008arXiv

Influence of primordial magnetic fields on 21 cm emission

Magnetic fields in the early universe can significantly alter the thermal evolution and the ionization history during the dark ages. This is reflected in the 21 cm line of atomic hydrogen, which is coupled to the gas temperature through collisions at high redshifts, and through the Wouthuysen-Field effect at low redshifts. We present a semi-analytic model for star formation and the build-up of a Lyman alpha background in the presence of magnetic fields, and calculate the evolution of the mean 21 cm brightness temperature and its frequency gradient as a function of redshift. We further discuss the evolution of linear fluctuations in temperature and ionization in the presence of magnetic fields and calculate the effect on the 21 cm power spectrum. At high redshifts, the signal is increased compared to the non-magnetic case due to the additional heat input into the IGM from ambipolar diffusion and the decay of MHD turbulence. At lower redshifts, the formation of luminous objects and the build-up of a Lyman alpha background can be delayed by a redshift interval of 10 due to the strong increase of the filtering mass scale in the presence of magnetic fields. This tends to decrease the 21 cm signal compared to the zero-field case. In summary, we find that 21 cm observations may become a promising tool to constrain primordial magnetic fields.