BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Robi Banerjee

Robi Banerjee appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRastro-ph.GAastro-ph.COphysics.flu-dynastro-ph.IMphysics.comp-phastro-phastro-ph.EPphysics.plasm-phastro-ph.HEhep-phnlin.CDphysics.data-an

Catalog footprint

What is connected

50works
13topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

50 published item(s)

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Efficient highly-subsonic turbulent dynamo and growth of primordial magnetic fields

We present the first study on the amplification of magnetic fields by the turbulent dynamo in the highly subsonic regime, with Mach numbers ranging from $10^{-3}$ to $0.4$. We find that for the lower Mach numbers the saturation efficiency of the dynamo, $(E_{\mathrm{mag}}/E_{\mathrm{kin}})_{\mathrm{sat}}$, increases as the Mach number decreases. Even in the case when injection of energy is purely through longitudinal forcing modes, $(E_{\mathrm{mag}}/E_{\mathrm{kin}})_{\mathrm{sat}}$ $\gtrsim 10^{-2}$ at a Mach number of $10^{-3}$. We apply our results to magnetic field amplification in the early Universe and predict that a turbulent dynamo can amplify primordial magnetic fields to $\gtrsim$ $10^{-16}$ Gauss on scales up to 0.1 pc and $\gtrsim$ $10^{-13}$ Gauss on scales up to 100 pc. This produces fields compatible with lower limits of the intergalactic magnetic field inferred from blazar $γ$-ray observations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the coexistence of the streaming instability and the vertical shear instability in protoplanetary disks

The streaming instability is a leading candidate mechanism to explain the formation of planetesimals. Yet, the role of this instability in the driving of turbulence in protoplanetary disks, given its fundamental nature as a linear hydrodynamical instability, has so far not been investigated in detail. We study the turbulence that is induced by the streaming instability as well as its interaction with the vertical shear instability. For this purpose, we employ the FLASH Code to conduct two-dimensional axisymmetric global disk simulations spanning radii from $1$ au to $100$ au, including the mutual drag between gas and dust as well as the radial and vertical stellar gravity. If the streaming instability and the vertical shear instability start their growth at the same time, we find the turbulence in the dust mid-plane layer to be primarily driven by the streaming instability. It gives rise to vertical gas motions with a Mach number of up to ${\sim}10^{-2}$. The dust scale height is set in a self-regulatory manner to about $1\%$ of the gas scale height. In contrast, if the vertical shear instability is allowed to saturate before the dust is introduced into our simulations, then it continues to be the main source of the turbulence in the dust layer. The vertical shear instability induces turbulence with a Mach number of ${\sim}10^{-1}$ and thus impedes dust sedimentation. Nonetheless, we find the vertical shear instability and the streaming instability in combination to lead to radial dust concentration in long-lived accumulations which are significantly denser than those formed by the streaming instability alone. Thus, the vertical shear instability may promote planetesimal formation by creating weak overdensities that act as seeds for the streaming instability.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Simulating the Formation of Massive Protostars: I. Radiative Feedback and Accretion Disks

We present radiation hydrodynamic simulations of collapsing protostellar cores with initial masses of 30, 100, and 200 M$_{\odot}$. We follow their gravitational collapse and the formation of a massive protostar and protostellar accretion disk. We employ a new hybrid radiative feedback method blending raytracing techniques with flux-limited diffusion for a more accurate treatment of the temperature and radiative force. In each case, the disk that forms becomes Toomre-unstable and develops spiral arms. This occurs between 0.35 and 0.55 freefall times and is accompanied by an increase in the accretion rate by a factor of 2-10. Although the disk becomes unstable, no other stars are formed. In the case of our 100 and 200 M$_{\odot}$ simulation, the star becomes highly super-Eddington and begins to drive bipolar outflow cavities that expand outwards. These radiatively-driven bubbles appear stable, and appear to be channeling gas back onto the protostellar accretion disk. Accretion proceeds strongly through the disk. After 81.4 kyr of evolution, our 30 M$_{\odot}$ simulation shows a star with a mass of 5.48 M$_{\odot}$ and a disk of mass 3.3 M$_{\odot}$, while our 100 M$_{\odot}$ simulation forms a 28.8 M$_{\odot}$ mass star with a 15.8 M$_{\odot}$ disk over the course of 41.6 kyr, and our 200 M$_{\odot}$ simulation forms a 43.7 M$_{\odot}$ star with an 18 M$_{\odot}$ disk in 21.9 kyr. In the absence of magnetic fields or other forms of feedback, the masses of the stars in our simulation do not appear limited by their own luminosities.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Supernova Feedback in Molecular Clouds: Global Evolution and Dynamics

We use magnetohydrodynamical simulations of converging warm neutral medium flows to analyse the formation and global evolution of magnetised and turbulent molecular clouds subject to supernova feedback from massive stars. We show that supernova feedback alone fails to disrupt entire, gravitationally bound, molecular clouds, but is able to disperse small--sized (~10 pc) regions on timescales of less than 1 Myr. Efficient radiative cooling of the supernova remnant as well as strong compression of the surrounding gas result in non-persistent energy and momentum input from the supernovae. However, if the time between subsequent supernovae is short and they are clustered, large hot bubbles form that disperse larger regions of the parental cloud. On longer timescales, supernova feedback increases the amount of gas with moderate temperatures (T~300-3000 K). Despite its inability to disrupt molecular clouds, supernova feedback leaves a strong imprint on the star formation process. We find an overall reduction of the star formation efficiency by a factor of 2 and of the star formation rate by roughly factors of 2-4.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A GPU accelerated Barnes-Hut Tree Code for FLASH4

We present a GPU accelerated CUDA-C implementation of the Barnes Hut (BH) tree code for calculating the gravitational potential on octree adaptive meshes. The tree code algorithm is implemented within the FLASH4 adaptive mesh refinement (AMR) code framework and therefore fully MPI parallel. We describe the algorithm and present test results that demonstrate its accuracy and performance in comparison to the algorithms available in the current FLASH4 version. We use a MacLaurin spheroid to test the accuracy of our new implementation and use spherical, collapsing cloud cores with effective AMR to carry out performance tests also in comparison with previous gravity solvers. Depending on the setup and the GPU/CPU ratio, we find a speedup for the gravity unit of at least a factor of 3 and up to 60 in comparison to the gravity solvers implemented in the FLASH4 code. We find an overall speedup factor for full simulations of at least factor 1.6 up to a factor of 10

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Extragalactic magnetic fields unlikely generated at the electroweak phase transition

In this letter we show that magnetic fields generated at the electroweak phase transition are most likely too weak to explain the void magnetic fields apparently observed today unless they have considerable helicity. We show that, in the simplest estimates, the helicity naturally produced in conjunction with the baryon asymmetry is too small to explain observations, which require a helicity fraction at least of order $10^{-14}-10^{-10}$ depending on the void fields constraint used. Therefore new mechanisms to generate primordial helicity are required if magnetic fields generated during the electroweak phase transition should explain the extragalactic fields.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Formation of star-forming clouds from the magnetised, diffuse interstellar medium

Molecular clouds, the birthplaces of stars in galaxies, form dynamically from the diffuse atomic gas of the interstellar medium (ISM). The ISM is also threaded by magnetic fields which have a large impact on its dynamics. In particular, star forming regions must be magnetically supercrit- ical in order to accomodate gas clumps which can collapse under their own weight. Based on a parameter study of three dimensional magneto-hydrodyamical (MHD) simulations, we show that the long-standing problem of how such supercritical regions are generated is still an open issue.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Impact Of Magnetic Fields On Molecular Cloud Formation & Evolution

We use magnetohydrodynamical simulations of converging flows to investigate the process of molecular cloud formation and evolution out of the magnetised ISM. Here, we investigate whether the observed subcritical HI clouds can become supercritical and hence allow the formation of stars within them. To do so, we vary the turbulent Mach number of the flows, as well as the initial magnetic field strength. We show that dense cores are able to build up under all conditions, but that star formation in these cores is either heavily delayed or completely suppressed if the initial field strength is B>3 microGauss. To probe the effect of magnetic diffusion, we introduce a tilting angle between the flows and the uniform background magnetic field, which mimics non--ideal MHD effects. Even with highly inclined flows, the formed cores are devoid of star formation, because no magnetically supercritical regions are build up. Hence we conclude, that the problem of how supercritical cloud cores are generated still persists.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Planet formation in post-common-envelope binaries

To understand the evolution of planetary systems, it is important to investigate planets in highly evolved stellar systems, and to explore the implications of their observed properties with respect to potential formation scenarios. Observations suggest the presence of giant planets in post-common-envelope binaries (PCEBs). A particularly well-studied system with planetary masses of 1.7 M_J and 7.0 M_J is NN Ser. We show here that a pure first-generation scenario where the planets form before the common envelope (CE) phase and the orbits evolve due to the changes in the gravitational potential is inconsistent with the current data. We propose a second-generation scenario where the planets are formed from the material that is ejected during the CE, which may naturally explain the observed planetary masses. In addition, hybrid scenarios where the planets form before the CE and evolve due to the accretion of the ejected gas appear as a realistic possibility.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Radiation Hydrodynamics using Characteristics on Adaptive Decomposed Domains for Massively Parallel Star Formation Simulations

We present an algorithm for solving the radiative transfer problem on massively parallel computers using adaptive mesh refinement and domain decomposition. The solver is based on the method of characteristics which requires an adaptive raytracer that integrates the equation of radiative transfer. The radiation field is split into local and global components which are handled separately to overcome the non-locality problem. The solver is implemented in the framework of the magneto-hydrodynamics code FLASH and is coupled by an operator splitting step. The goal is the study of radiation in the context of star formation simulations with a focus on early disc formation and evolution. This requires a proper treatment of radiation physics that covers both the optically thin as well as the optically thick regimes and the transition region in particular. We successfully show the accuracy and feasibility of our method in a series of standard radiative transfer problems and two 3D collapse simulations resembling the early stages of protostar and disc formation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A general hybrid radiation transport scheme for star formation simulations on an adaptive grid

Radiation feedback plays a crucial role in the process of star formation. In order to simulate the thermodynamic evolution of disks, filaments, and the molecular gas surrounding clusters of young stars, we require an efficient and accurate method for solving the radiation transfer problem. We describe the implementation of a hybrid radiation transport scheme in the adaptive grid-based FLASH general magnetohydrodynamics code. The hybrid scheme splits the radiative transport problem into a raytracing step and a diffusion step. The raytracer captures the first absorption event, as stars irradiate their environments, while the evolution of the diffuse component of the radiation field is handled by a flux-limited diffusion (FLD) solver. We demonstrate the accuracy of our method through a variety of benchmark tests including the irradiation of a static disk, subcritical and supercritical radiative shocks, and thermal energy equilibration. We also demonstrate the capability of our method for casting shadows and calculating gas and dust temperatures in the presence of multiple stellar sources. Our method enables radiation-hydrodynamic studies of young stellar objects, protostellar disks, and clustered star formation in magnetized, filamentary environments.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Magnetic field amplification by the small-scale dynamo in the early Universe

In this paper we show that the Universe is already strongly magnetized at very early epochs during cosmic evolution. Our calculations are based on the efficient amplification of weak magnetic seed fields, which are unavoidably present in the early Universe, by the turbulent small-scale dynamo. We identify two mechanisms for the generation of turbulence in the radiation dominated epoch where velocity fluctuations are produced by the primordial density perturbation and by possible first-order phase transitions at the electroweak or QCD scales. We show that all the necessities for the small-scale dynamo to work are fulfilled. Hence, this mechanism, operating due to primordial density perturbations, guarantees fields with comoving field strength $B_0\sim10^{-6}\varepsilon^{1/2}$ nG on scales up to $λ_c\sim0.1$ pc, where $\varepsilon$ is the saturation efficiency. The amplification of magnetic seed fields could be even larger if there are first-order phase transitions in the early Universe. Where, on scales up to $λ_c\sim100$ pc, the comoving field strength due to this mechanism will be $B_0\sim10^{-3}\varepsilon^{1/2}$ nG at the present time. Such fields, albeit on small scales, can play an important role in structure formation and could provide an explanation to the apparently observed magnetic fields in the voids of the large-scale structure.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Modeling jet and outflow feedback during star cluster formation

Powerful jets and outflows are launched from the protostellar disks around newborn stars. These outflows carry enough mass and momentum to transform the structure of their parent molecular cloud and to potentially control star formation itself. Despite their importance, we have not been able to fully quantify the impact of jets and outflows during the formation of a star cluster. The main problem lies in limited computing power. We would have to resolve the magnetic jet-launching mechanism close to the protostar and at the same time follow the evolution of a parsec-size cloud for a million years. Current computer power and codes fall orders of magnitude short of achieving this. In order to overcome this problem, we implement a subgrid-scale (SGS) model for launching jets and outflows, which demonstrably converges and reproduces the mass, linear and angular momentum transfer, and the speed of real jets, with ~ 1000 times lower resolution than would be required without SGS model. We apply the new SGS model to turbulent, magnetized star cluster formation and show that jets and outflows (1) eject about 1/4 of their parent molecular clump in high-speed jets, quickly reaching distances of more than a parsec, (2) reduce the star formation rate by about a factor of two, and (3) lead to the formation of ~ 1.5 times as many stars compared to the no-outflow case. Most importantly, we find that jets and outflows reduce the average star mass by a factor of ~ 3 and may thus be essential for understanding the characteristic mass of the stellar initial mass function.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Star formation efficiency in turbulent clouds

Here we present a simple, but nevertheless, instructive model for the star formation efficiency in turbulent molecular clouds. The model is based on the assumption of log-normal density distribution which reflects the turbulent nature of the interstellar medium (ISM). Together with the number count of cloud cores, which follows a Salpeter-like core mass function (CMF), and the minimum mass for the collapse of individual cloud cores, given by the local Jeans mass, we are able to derive the SFE for clouds as a function of their Jeans masses. We find a very generic power-law, SFE \propto (M_cloud/M_J)^{-0.26} and a maximum SFE_max \sim 1/3 for the Salpeter case. This result is independent of the turbulent Mach number but fairly sensitive to variations of the CMF.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Earliest Stages of Star and Planet Formation: Core Collapse, and the Formation of Disks and Outflows

(Abridged) In this review we focus on the observations and theory of the formation of early disks and outflows, and their connections with the first phases of planet formation. Large rotationally supported circumstellar disks, although common around more evolved young stellar objects, are rarely detected during the earliest, "Class 0" phase; however, a few excellent candidates have been discovered recently around both low and high mass protostars. In this early phase, prominent outflows are ubiquitously observed; they are expected to be associated with at least small magnetized disks. Disk formation - once thought to be a simple consequence of the conservation of angular momentum during hydrodynamic core collapse - is far more subtle in magnetized gas. In this case, the rotation can be strongly magnetically braked. Indeed, both analytic arguments and numerical simulations have shown that disk formation is suppressed in the strict ideal magnetohydrodynamic (MHD) limit for the observed level of core magnetization. We review what is known about this "magnetic braking catastrophe", possible ways to resolve it, and the current status of early disk observations. Outflows are also intimately linked to disk formation; they are a natural product of magnetic fields and rotation and are important signposts of star formation. We review new developments on early outflow generation since PPV. The properties of early disks and outflows are a key component of planet formation in its early stages and we review these major connections.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Cloud formation in colliding flows: influence of the choice of cooling function

We study the influence of the choice of cooling function on the formation of molecular clouds in high-resolution three-dimensional simulations of converging flows. We directly compare the results obtained using the simple, parametrized cooling function introduced by Koyama & Inutsuka (2002) and used by a number of converging flow studies with the results of the detailed calculation of the non-equilibrium chemistry and thermal balance of the gas. We find that a number of the cloud properties, such as the mass and volume filling fractions of cold gas, are relatively insensitive to the choice of cooling function. On the other hand, the cloud morphology and the large-scale velocity distribution of the gas do strongly depend on the cooling function. We show that the differences that we see can largely be explained by differences in the way that Lyman-alpha cooling is treated in the two complementary approaches, and that a proper non-equilibrium treatment of the ionisation and recombination of the gas is necessary in order to model the high-temperature cooling correctly. We also investigate the properties of the dense clumps formed within the cloud. In agreement with previous models, we find that the majority of these clumps are not self-gravitating, suggesting that some form of large-scale collapse of the cloud may be required in order to produce gravitationally unstable clumps and hence stars. Overall, the physical properties of the dense clumps are similar in both simulations, suggesting that they do not depend strongly on the choice of cooling function. However, we do find a systematic difference of around 10K in the mean temperatures of the clumps produced by the two models.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Morphologies of protostellar outflows: An ALMA view

The formation of stars is usually accompanied by the launching of protostellar outflows. Observations with the Atacama Large Millimetre/sub-millimetre Array (ALMA) will soon revolutionalise our understanding of the morphologies and kinematics of these objects. In this paper, we present synthetic ALMA observations of protostellar outflows based on numerical magnetohydrodynamic collapse simulations. We find significant velocity gradients in our outflow models and a very prominent helical structure within the outflows. We speculate that the disk wind found in the ALMA Science Verification Data of HD 163296 presents a first instance of such an observation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Are molecular outflows around high-mass stars driven by ionization feedback?

The formation of massive stars exceeding 10 solar masses usually results in large-scale molecular outflows. Numerical simulations, including ionization, of the formation of such stars show evidence for ionization-driven molecular outflows. We here examine whether the outflows seen in these models reproduce the observations. We compute synthetic ALMA and CARMA maps of CO emission lines of the outflows, and compare their signatures to existing single-dish and interferometric data. We find that the ionization-driven models can only reproduce weak outflows around high-mass star-forming regions. We argue that expanding H II regions probably do not represent the dominant mechanism for driving observed outflows. We suggest instead that observed outflows are driven by the collective action of the outflows from the many lower-mass stars that inevitably form around young massive stars in a cluster.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Magnetic Field Amplification by Small-Scale Dynamo Action: Dependence on Turbulence Models and Reynolds and Prandtl Numbers

The small-scale dynamo is a process by which turbulent kinetic energy is converted into magnetic energy, and thus is expected to depend crucially on the nature of turbulence. In this work, we present a model for the small-scale dynamo that takes into account the slope of the turbulent velocity spectrum v(l) ~ l^theta, where l and v(l) are the size of a turbulent fluctuation and the typical velocity on that scale. The time evolution of the fluctuation component of the magnetic field, i.e., the small-scale field, is described by the Kazantsev equation. We solve this linear differential equation for its eigenvalues with the quantum-mechanical WKB-approximation. The validity of this method is estimated as a function of the magnetic Prandtl number Pm. We calculate the minimal magnetic Reynolds number for dynamo action, Rm_crit, using our model of the turbulent velocity correlation function. For Kolmogorov turbulence (theta=1/3), we find that the critical magnetic Reynolds number is approximately 110 and for Burgers turbulence (theta=1/2) approximately 2700. Furthermore, we derive that the growth rate of the small-scale magnetic field for a general type of turbulence is Gamma ~ Re^((1-theta)/(1+theta)) in the limit of infinite magnetic Prandtl numbers. For decreasing magnetic Prandtl number (down to Pm approximately larger than 10), the growth rate of the small-scale dynamo decreases. The details of this drop depend on the WKB-approximation, which becomes invalid for a magnetic Prandtl number of about unity.

0 citations0 reviewsNo code linkedOpen access
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)

0 citations0 reviewsNo code linkedOpen access
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}.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The long-term evolution of warped, magnetised discs, and precessing outflows in collapsing pre-stellar cores

(abridged) The nature of early Class 0/I protostellar discs is not clearly understood. Early protostellar discs are needed to drive molecular outflows and jets observed in star forming regions, but there has been some debate to how they form. From a theoretical perspective, the consequences of disc and outflow generation are crucial to understanding the very nature of how stars are assembled. We have performed 3D ideal magnetohydrodynamic (MHD) simulations of collapsing Bonnor-Ebert spheres, employing sink particles with a radius of 3.2 AU alongside an AMR grid and using a cooling function to model radiative cooling of the gas. This has allowed us to explore 2-8x10^4 yr further into the evolution of an early Class 0 disc-outflow system than previous simulations. Our outflow is precessing, kinked, turbulent, contains internal shocks and has a scale of 0.1 pc end-to-end. We form a rotationally dominated disc with a radius of 100 AU embedded inside a transient, unstable, flattened, rotating core extending out to 2000 AU. The larger flattened structure launches a low speed wind (v_r < 1.5 km/s) dominated by B_phi, while the inner disc launches a centrifugally driven jet dominated by B_p with speeds up to 20$ km/s. From the inner disk, the value of dM/dt_out/dM/dt_in ~ 0.1, wheras in the outer core dM/dt_out/dM/dt_in ~ 1.0. The inner disc becomes unstable to a warping instability due to the magnetic structure of the outflow and warps to 30 deg with respect to the z-axis by the end of the simulation. The envelope is cleared out and is less massive than the disc. We measure star formation efficiencies of eta_core=0.63 (and growing), higher than theoretical predictions. This indicates that outflows are not as efficient at expelling envelope mass as some current models estimate. We discuss the relevance of our disc misalignment concerning the formation of mis-aligned hot Jupiters.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Small-Scale Dynamo and Non-Ideal MHD in Primordial Star Formation

We study the amplification of magnetic fields during the formation of primordial halos. The turbulence generated by gravitational infall motions during the formation of the first stars and galaxies can amplify magnetic fields very efficiently and on short timescales up to dynamically significant values. Using the Kazantsev theory, which describes the so-called small-scale dynamo - a magnetohydrodynamical process converting kinetic energy from turbulence into magnetic energy - we can then calculate the growth rate of the small-scale magnetic field. Our calculations are based on a detailed chemical network and we include non-ideal magnetohydrodynamical effects such as ambipolar diffusion and Ohmic dissipation. We follow the evolution of the magnetic field up to larger scales until saturation occurs on the Jeans scale. Assuming a weak magnetic seed field generated by the Biermann battery process, both Burgers and Kolmogorov turbulence lead to saturation within a rather small density range. Such fields are likely to become relevant after the formation of a protostellar disk and, thus, could influence the formation of the first stars and galaxies in the Universe.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Comparing Numerical Methods for Isothermal Magnetized Supersonic Turbulence

We employ simulations of supersonic super-Alfvenic turbulence decay as a benchmark test problem to assess and compare the performance of nine astrophysical MHD methods actively used to model star formation. The set of nine codes includes: ENZO, FLASH, KT-MHD, LL-MHD, PLUTO, PPML, RAMSES, STAGGER, and ZEUS. We present a comprehensive set of statistical measures designed to quantify the effects of numerical dissipation in these MHD solvers. We compare power spectra for basic fields to determine the effective spectral bandwidth of the methods and rank them based on their relative effective Reynolds numbers. We also compare numerical dissipation for solenoidal and dilatational velocity components to check for possible impacts of the numerics on small-scale density statistics. Finally, we discuss convergence of various characteristics for the turbulence decay test and impacts of various components of numerical schemes on the accuracy of solutions. We show that the best performing codes employ a consistently high order of accuracy for spatial reconstruction of the evolved fields, transverse gradient interpolation, conservation law update step, and Lorentz force computation. The best results are achieved with divergence-free evolution of the magnetic field using the constrained transport method, and using little to no explicit artificial viscosity. Codes which fall short in one or more of these areas are still useful, but they must compensate higher numerical dissipation with higher numerical resolution. This paper is the largest, most comprehensive MHD code comparison on an application-like test problem to date. We hope this work will help developers improve their numerical algorithms while helping users to make informed choices in picking optimal applications for their specific astrophysical problems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Importance of the Initial Conditions for Star Formation - I. Cloud Evolution and Morphology

We present a detailed parameter study of collapsing turbulent cloud cores, varying the initial density profile and the initial turbulent velocity field. We systematically investigate the influence of different initial conditions on the star formation process, mainly focusing on the fragmentation, the number of formed stars, and the resulting mass distributions. Our study compares four different density profiles (uniform, Bonnor-Ebert type, $ρ\propto r^{-1.5}$, and $ρ\propto r^{-2}$), combined with six different supersonic turbulent velocity fields (compressive, mixed, and solenoidal, initialised with two different random seeds each) in three-dimensional simulations using the adaptive-mesh refinement, hydrodynamics code FLASH. The simulations show that density profiles with flat cores produce hundreds of low-mass stars, either distributed throughout the entire cloud or found in subclusters, depending on the initial turbulence. Concentrated density profiles always lead to the formation of one high-mass star in the centre of the cloud and, if at all, low-mass stars surrounding the central one. In uniform and Bonnor-Ebert type density distributions, compressive initial turbulence leads to local collapse about 25% earlier than solenoidal turbulence. However, central collapse in the steep power-law profiles is too fast for the turbulence to have any significant influence. We conclude that (I) the initial density profile and turbulence mainly determine the cloud evolution and the formation of clusters, (II) the initial mass function (IMF) is not universal for all setups, and (III) that massive stars are much less likely to form in flat density distributions. The IMFs obtained in the uniform and Bonnor-Ebert type density profiles are more consistent with the observed IMF, but shifted to lower masses.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Importance of the Initial Conditions for Star Formation - II. Fragmentation Induced Starvation and Accretion Shielding

We investigate the impact of different initial conditions for the initial density profile and the initial turbulence on the formation process of protostellar clusters. We study the collapse of dense molecula r cloud cores with three-dimensional adaptive mesh refinement simulations. We focus our discussion on the distribution of the gas among the protostellar objects in the turbulent dynamical cores. Despite the large variations in the initial configurations and the resulting gas and cluster morphology we find that all stellar clusters follow a very similar gas accretion behaviour. Once secondary protostars begin to form, the central region of a cluster is efficiently shielded from further accretion. Hence, objects located close to the centre are starved of material, as indicated by a strong decrease of the central acc retion rate. This Fragmentation Induced Starvation occurs not only in rotationally supported discs and filaments, but also in more spherically symmetric clusters with complex chaotic motions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Importance of the Initial Conditions for Star Formation - III: Statistical Properties of Embedded Protostellar Clusters

We investigate the formation of protostellar clusters during the collapse of dense molecular cloud cores with a focus on the evolution of potential and kinetic energy, the degree of substructure, and the early phase of mass segregation. Our study is based on a series of hydrodynamic simulations of dense cores, where we vary the initial density profile and the initial turbulent velocity. In the three-dimensional adaptive mesh refinement simulations, we follow the dynamical formation of filaments and protostars until a star formation efficiency of 20%. Despite the different initial configurations, the global ensemble of all protostars in a setup shows a similar energy evolution and forms sub-virial clusters with an energy ratio $E_\mathrm{kin}/|E_\mathrm{pot}|\sim0.2$. Concentrating on the innermost central region, the clusters show a roughly virialised energy balance. However, the region of virial balance only covers the innermost $\sim10-30%$ of all the protostars. In all simulations with multiple protostars, the total kinetic energy of the protostars is higher than the kinetic energy of the gas cloud, although the protostars only contain 20% of the total mass. The clusters vary significantly in size, mass, and number of protostars, and show different degrees of substructure and mass segregation. Flat density profiles and compressive turbulent modes produce more subclusters then centrally concentrated profiles and solenoidal turbulence. We find that dynamical relaxation and hence dynamical mass segregation is very efficient in all cases from the very beginning of the nascent cluster, i.e., during a phase when protostars are constantly forming and accreting.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Molecular Cloud Evolution IV: Magnetic Fields, Ambipolar Diffusion, and the Star Formation Efficiency

We investigate the formation and evolution of giant molecular clouds (GMCs) by the collision of convergent warm neutral medium (WNM) streams in the interstellar medium, in the presence of magnetic fields and ambipolar diffusion (AD), focusing on the evolution of the star formation rate (SFR) and efficiency (SFE), as well as of the mass-to-magnetic-flux ratio (M2FR) in the forming clouds. We find that: 1) Clouds formed by supercritical inflow streams proceed directly to collapse, while clouds formed by subcritical streams first contract and then re-expand, oscillating on the scale of tens of Myr. 2) Our suite of simulations with initial magnetic field strength of 2, 3, and 4 $μ\G$ show that only supercritical or marginal critical streams lead to reasonable star forming rates. 3) The GMC's M2FR is a generally increasing function of time, whose growth rate depends on the details of how mass is added to the GMC from the WNM. 4) The M2FR is a highly fluctuating function of position in the clouds. 5) In our simulations, the SFE approaches stationarity, because mass is added to the GMC at a similar rate at which it converts mass to stars. In such an approximately stationary regime, the SFE provides a proxy of the supercritical mass fraction in the cloud. 6) We observe the occurrence of buoyancy of the low-M2FR regions within the gravitationally-contracting GMCs, so that the latter naturally segregate into a high-density, high-M2FR "core" and a low-density, low-M2FR "envelope", without the intervention of AD. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Radiative Feedback in Massive Star and Cluster Formation

Understanding the origin of high-mass stars is central to modern astrophysics. We shed light on this problem with simulations using a novel, adaptive-mesh, ray-tracing algorithm. These simulations consistently follow the gravitational collapse of a massive molecular cloud core, the subsequent build-up and fragmentation of the accretion disk surrounding the nascent star, and, for the first time, the interaction between its intense UV radiation field and the infalling material. We show that ionization feedback can neither stop protostellar mass growth nor suppress fragmentation. We discuss the effects of feedback by ionizing and non-ionizing radiation on the evolution of the stellar cluster. The accretion is not limited by radiative feedback but by the formation of low-mass companions in a process we call "fragmentation-induced starvation". This behavior consistently reproduces the observed relation between the most massive star and the total mass of stars in a cluster. We show that magnetic fields reduce the star formation rate and lead to the formation of more massive stars.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Statistical analysis of the mass-to-flux ratio in turbulent cores: effects of magnetic field reversals and dynamo amplification

We study the mass-to-flux ratio (M/Φ) of clumps and cores in simulations of supersonic, magnetohydrodynamical turbulence for different initial magnetic field strengths. We investigate whether the (M/Φ)-ratio of core and envelope, R = (M/Φ)_{core}/(M/Φ)_{envelope} can be used to distinguish between theories of ambipolar diffusion and turbulence-regulated star formation. We analyse R for different Lines-of-Sight (LoS) in various sub-cubes of our simulation box. We find that, 1) the average and median values of |R| for different times and initial magnetic field strengths are typically greater, but close to unity, 2) the average and median values of |R| saturate at average values of |R| ~ 1 for smaller magnetic fields, 3) values of |R| < 1 for small magnetic fields in the envelope are caused by field reversals when turbulence twists the field lines such that field components in different directions average out. Finally, we propose two mechanisms for generating values |R| ~< 1 for the weak and strong magnetic field limit in the context of a turbulent model. First, in the weak field limit, the small-scale turbulent dynamo leads to a significantly increased flux in the core and we find |R| ~< 1. Second, in the strong field limit, field reversals in the envelope also lead to values |R| ~< 1. These reversals are less likely to occur in the core region where the velocity field is more coherent and the internal velocity dispersion is typically subsonic.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Interplay of Magnetic Fields, Fragmentation and Ionization Feedback in High-Mass Star Formation

Massive stars disproportionately influence their surroundings. How they form has only started to become clear recently through radiation gas dynamical simulations. However, until now, no simulation has simultaneously included both magnetic fields and ionizing radiation. Here we present the results from the first radiation-magnetohydrodynamical (RMHD) simulation including ionization feedback, comparing an RMHD model of a 1000 M_sol rotating cloud to earlier radiation gas dynamical models with the same initial density and velocity distributions. We find that despite starting with a strongly supercritical mass to flux ratio, the magnetic field has three effects. First, the field offers locally support against gravitational collapse in the accretion flow, substantially reducing the amount of secondary fragmentation in comparison to the gas dynamical case. Second, the field drains angular momentum from the collapsing gas, further increasing the amount of material available for accretion by the central, massive, protostar, and thus increasing its final mass by about 50% from the purely gas dynamical case. Third, the field is wound up by the rotation of the flow, driving a tower flow. However, this flow never achieves the strength seen in low-mass star formation simulations for two reasons: gravitational fragmentation disrupts the circular flow in the central regions where the protostars form, and the expanding H II regions tend to further disrupt the field geometry. Therefore, outflows driven by ionization heating look likely to be more dynamically important in regions of massive star formation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Time Variability in Simulated Ultracompact and Hypercompact HII Regions

Ultracompact and hypercompact HII regions appear when a star with a mass larger than about 15 solar masses starts to ionize its own environment. Recent observations of time variability in these objects are one of the pieces of evidence that suggest that at least some of them harbor stars that are still accreting from an infalling neutral accretion flow that becomes ionized in its innermost part. We present an analysis of the properties of the HII regions formed in the 3D radiation-hydrodynamic simulations presented by Peters et al. as a function of time. Flickering of the HII regions is a natural outcome of this model. The radio-continuum fluxes of the simulated HII regions, as well as their flux and size variations are in agreement with the available observations. From the simulations, we estimate that a small but non-negligible fraction (~ 10 %) of observed HII regions should have detectable flux variations (larger than 10 %) on timescales of ~ 10 years, with positive variations being more likely to happen than negative variations. A novel result of these simulations is that negative flux changes do happen, in contrast to the simple expectation of ever growing HII regions. We also explore the temporal correlations between properties that are directly observed (flux and size) and other quantities like density and ionization rates.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

H II regions: Witnesses to massive star formation

We describe the first three-dimensional simulation of the gravitational collapse of a massive, rotating molecular cloud that includes heating by both non-ionizing and ionizing radiation. We find that as the first protostars gain sufficient mass to ionize the accretion flow, their H II regions are initially gravitationally trapped, but soon begin to rapidly fluctuate between trapped and extended states, in agreement with observations. Over time, the same ultracompact H II region can expand anisotropically, contract again, and take on any of the observed morphological classes. In their extended phases, expanding H II regions drive bipolar neutral outflows characteristic of high-mass star formation. The total lifetime of H II regions is given by the global accretion timescale, rather than their short internal sound-crossing time. The pressure of the hot, ionized gas does not terminate accretion. Instead the final stellar mass is set by fragmentation-induced starvation. Local gravitational instabilities in the accretion flow lead to the build-up of a small cluster of stars, all with relatively high masses due to heating from accretion radiation. These companions subsequently compete with the initial high-mass star for the same common gas reservoir and limit its mass growth. Our findings show that the most significant differences between the formation of low-mass and high-mass stars are all explained as the result of rapid accretion within a dense, gravitationally unstable, ionized flow.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Limiting Accretion onto Massive Stars by Fragmentation-Induced Starvation

Massive stars influence their surroundings through radiation, winds, and supernova explosions far out of proportion to their small numbers. However, the physical processes that initiate and govern the birth of massive stars remain poorly understood. Two widely discussed models are monolithic collapse of molecular cloud cores and competitive accretion. To learn more about massive star formation, we perform simulations of the collapse of rotating, massive, cloud cores including radiative heating by both non-ionizing and ionizing radiation using the FLASH adaptive mesh refinement code. These simulations show fragmentation from gravitational instability in the enormously dense accretion flows required to build up massive stars. Secondary stars form rapidly in these flows and accrete mass that would have otherwise been consumed by the massive star in the center, in a process that we term fragmentation-induced starvation. This explains why massive stars are usually found as members of high-order stellar systems that themselves belong to large clusters containing stars of all masses. The radiative heating does not prevent fragmentation, but does lead to a higher Jeans mass, resulting in fewer and more massive stars than would form without the heating. This mechanism reproduces the observed relation between the total stellar mass in the cluster and the mass of the largest star. It predicts strong clumping and filamentary structure in the center of collapsing cores, as has recently been observed. We speculate that a similar mechanism will act during primordial star formation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Modeling Collapse and Accretion in Turbulent Gas Clouds: Implementation and Comparison of Sink Particles in AMR and SPH

We implemented sink particles in the adaptive mesh refinement (AMR) hydrodynamics code FLASH. Sink particles are created in regions of local gravitational collapse, and their trajectories and accretion can be followed over many dynamical times. We perform a series of tests including the time integration of circular and elliptical orbits, the collapse of a Bonnor-Ebert sphere and a rotating, fragmenting cloud core. We compare the collapse of a highly unstable singular isothermal sphere to the theory by Shu (1977), and show that the sink particle accretion rate is in excellent agreement with the theoretical prediction. To model eccentric orbits and close encounters of sink particles accurately, we show that a very small timestep is often required, for which we implemented subcycling of the N-body system. We emphasize that a sole density threshold for sink particle creation is insufficient in supersonic flows, if the density threshold is below the opacity limit. In that case, the density can exceed the threshold in strong shocks that do not necessarily lead to local collapse. Additional checks for bound state, gravitational potential minimum, Jeans instability and converging flows are absolutely necessary for a meaningful creation of sink particles. We apply our new sink particle module for FLASH to the formation of a stellar cluster, and compare to a smoothed particle hydrodynamics (SPH) code with sink particles. Our comparison shows encouraging agreement of gas properties, indicated by column density distributions and radial profiles, and of sink particle formation times and positions. We find excellent agreement in the number of sink particles formed, and in their accretion and mass distributions.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Understanding Spatial and Spectral Morphologies of Ultracompact H II Regions

The spatial morphology, spectral characteristics, and time variability of ultracompact H II regions provide strong constraints on the process of massive star formation. We have performed simulations of the gravitational collapse of rotating molecular cloud cores, including treatments of the propagation of ionizing and non-ionizing radiation. We here present synthetic radio continuum observations of H II regions from our collapse simulations, to investigate how well they agree with observation, and what we can learn about how massive star formation proceeds. We find that intermittent shielding by dense filaments in the gravitationally unstable accretion flow around the massive star leads to highly variable H II regions that do not grow monotonically, but rather flicker, growing and shrinking repeatedly. This behavior appears able to resolve the well-known lifetime problem. We find that multiple ionizing sources generally form, resulting in groups of ultracompact H II regions, consistent with observations. We confirm that our model reproduces the qualitative H II region morphologies found in surveys, with generally consistent relative frequencies. We also find that simulated spectral energy distributions (SEDs) from our model are consistent with the range of observed H II region SEDs, including both regions showing a normal transition from optically thick to optically thin emission, and those with intermediate spectral slopes. In our models, anomalous slopes are solely produced by inhomogeneities in the H II region, with no contribution from dust emission at millimeter or submillimeter wavelengths. We conclude that many observed characteristics of ultracompact H II regions appear consistent with massive star formation in fast, gravitationally unstable, accretion flows.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Clump morphology and evolution in MHD simulations of molecular cloud formation

Abridged: We study the properties of clumps formed in three-dimensional weakly magnetized magneto-hydrodynamic simulations of converging flows in the thermally bistable, warm neutral medium (WNM). We find that: (1) Similarly to the situation in the classical two-phase medium, cold, dense clumps form through dynamically-triggered thermal instability in the compressed layer between the convergent flows, and are often characterised by a sharp density jump at their boundaries though not always. (2) However, the clumps are bounded by phase-transition fronts rather than by contact discontinuities, and thus they grow in size and mass mainly by accretion of WNM material through their boundaries. (3) The clump boundaries generally consist of thin layers of thermally unstable gas, but these layers are often widened by the turbulence, and penetrate deep into the clumps. (4) The clumps are approximately in both ram and thermal pressure balance with their surroundings, a condition which causes their internal Mach numbers to be comparable to the bulk Mach number of the colliding WNM flows. (5) The clumps typically have mean temperatures 20 < T < 50 K, corresponding to the wide range of densities they contain (20 < n < 5000 pcc) under a nearly-isothermal equation of state. (6) The turbulent ram pressure fluctuations of the WNM induce density fluctuations that then serve as seeds for local gravitational collapse within the clumps. (7) The velocity and magnetic fields tend to be aligned with each other within the clumps, although both are significantly fluctuating, suggesting that the velocity tends to stretch and align the magnetic field with it. (8) The typical mean field strength in the clumps is a few times larger than that in the WNM. (9) The magnetic field strength has a mean value of B ~ 6 mu G ...

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Massive star formation via high accretion rates and early disk-driven outflows

We present an investigation of massive star formation that results from the gravitational collapse of massive, magnetized molecular cloud cores. We investigate this by means of highly resolved, numerical simulations of initial magnetized Bonnor-Ebert-Spheres that undergo collapse and cooling. By comparing three different cases - an isothermal collapse, a collapse with radiative cooling, and a magnetized collapse - we show that massive stars assemble quickly with mass accretion rates exceeding 10^-3 Msol/yr. We confirm that the mass accretion during the collapsing phase is much more efficient than predicted by selfsimilar collapse solutions, i.e. dM/dt ~ c^3/G. We find that during protostellar assembly the mass accretion reaches 20 - 100 c^3/G. Furthermore, we determined the self-consistent structure of bipolar outflows that are produced in our three dimensional magnetized collapse simulations. These outflows produce cavities out of which radiation pressure can be released, thereby reducing the limitations on the final mass of massive stars formed by gravitational collapse. Moreover, we argue that the extraction of angular momentum by disk-threaded magnetic fields and/or by the appearance of bars with spiral arms significantly enhance the mass accretion rate, thereby helping the massive protostar to assemble more quickly.

0 citations0 reviewsNo code linkedOpen access