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

Christoph Federrath

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

ResearcherUnclaimed source record
astro-ph.GAastro-ph.SRphysics.flu-dynastro-ph.IMastro-ph.COphysics.comp-phphysics.plasm-phastro-ph.HEphysics.data-anDistributed, Parallel, and Cluster Computinghep-phmath-phmath.MPnlin.CDPerformancephysics.ao-ph

Catalog footprint

What is connected

94works
16topics
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

94 published item(s)

preprint2026arXiv

Turbulence driving in a star-forming Milky-Way-type galaxy

The life-cycle, structure, and dynamics of the interstellar medium (ISM) is regulated by turbulence. Complex physical processes, including supernova (SN) explosions, shear, and gravitational collapse, drive and maintain turbulence, but it is still an open question what turbulence driving mode is primarily excited by these different mechanisms. The turbulence driving parameter, b, can be used to quantify the ratio of solenoidal to compressive modes in the acceleration field that drives the turbulence. Compressive driving is characterised by b ~ 1, while purely solenoidal driving gives b ~ 0.3. To quantify the turbulence in the galactic ISM, we investigate the time evolution of b, as well as the turbulent Mach number, and plasma beta (thermal-to-magnetic pressure ratio), and its correlation with star formation in the magnetised warm neutral medium (WNM) of the TIGRESS shearing-box simulations of a kpc-sized patch of a Milky-Way-like galaxy, over a 100 Myr time period (~ half an orbital time). In this simulation the turbulence is driven by a combination of shear, gravitational collapse, and star formation feedback in the form of radiation and SNe. We find that the turbulence driving parameter fluctuates in time between b ~ 0.4 and b ~ 1. We find a time-dependent correlation of b with star formation activity, such that high star formation rates follow about one turbulent turnover time (~ 10 Myr) after phases of highly compressive driving (b > 0.5). About 20 Myr after the peak in star formation, type-B SN feedback drives up the WNM fraction and turbulent Mach numbers, and reduces plasma beta and the driving to b ~ 0.4-0.5.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

What Sets the Star Formation Rate of Molecular Clouds? The Density Distribution as a Fingerprint of Compression and Expansion Rates

We use a suite of 3D simulations of star-forming molecular clouds, with and without stellar feedback, magnetic fields, and driven turbulence, to study the compression and expansion rates of the gas as functions of density. We show that, around the mean density, supersonic turbulence promotes rough equilibrium between the amounts of compressing and expanding gas, consistent with continuous gas cycling between high and low density states. We find that the inclusion of protostellar jets produces rapidly expanding and compressing low-density gas. We find that the gas mass flux peaks at the transition between the lognormal and power-law forms of the density probability distribution function (PDF). This is consistent with the transition density tracking the post-shock density, which promotes an enhancement of mass at this density (i.e., shock compression and filament formation). At high densities, the gas dynamics are dominated by self-gravity: the compression rate in all of our runs matches the rate of the run with only gravity, suggesting that processes other than self-gravity have little effect at these densities. The net gas mass flux becomes constant at a density below the sink formation threshold, where it equals the star formation rate. The density at which the net gas mass flux equals the star formation rate is one order of magnitude lower than our sink threshold density, corresponds to the formation of the second power-law tail in the density PDF, and sets the overall star formation rates of these simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Energy balance and Alfvén Mach numbers in compressible magnetohydrodynamic turbulence with a large-scale magnetic field

Energy equipartition is a powerful theoretical tool for understanding astrophysical plasmas. It is invoked, for example, to measure magnetic fields in the interstellar medium (ISM), as evidence for small-scale turbulent dynamo action, and, in general, to estimate the energy budget of star-forming molecular clouds. In this study we motivate and explore the role of the volume-averaged root-mean-squared (rms) magnetic coupling term between the turbulent, $δ\mathbf{B}$ and large-scale, $\mathbf{B}_0$ fields, $\left< (δ\mathbf{B}\cdot\mathbf{B}_0)^{2} \right>^{1/2}_{\mathcal{V}}$. By considering the second moments of the energy balance equations we show that the rms coupling term is in energy equipartition with the volume-averaged turbulent kinetic energy for turbulence with a sub-Alfvénic large-scale field. Under the assumption of exact energy equipartition between these terms, we derive relations for the magnetic and coupling term fluctuations, which provide excellent, parameter-free agreement with time-averaged data from 280 numerical simulations of compressible MHD turbulence. Furthermore, we explore the relation between the turbulent, mean-field and total Alfvén Mach numbers, and demonstrate that sub-Alfvénic turbulence can only be developed through a strong, large-scale magnetic field, which supports an extremely super-Alfvénic turbulent magnetic field. This means that the magnetic field fluctuations are significantly subdominant to the velocity fluctuations in the sub-Alfvénic large-scale field regime. Throughout our study, we broadly discuss the implications for observations of magnetic fields and understanding the dynamics in the magnetised ISM.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Fundamental scales in the kinematic phase of the turbulent dynamo

The turbulent dynamo is a powerful mechanism that converts turbulent kinetic energy to magnetic energy. A key question regarding the magnetic field amplification by turbulence, is, on what scale, $k_{\rm p}$, do magnetic fields become most concentrated? There has been some disagreement about whether $k_{\rm p}$ is controlled by the viscous scale, $k_ν$ (where turbulent kinetic energy dissipates), or the resistive scale, $k_η$ (where magnetic fields dissipate). Here we use direct numerical simulations of magnetohydrodynamic turbulence to measure characteristic scales in the kinematic phase of the turbulent dynamo. We run $104$-simulations with hydrodynamic Reynolds numbers of $10 \leq {\rm Re} \leq 3600$, and magnetic Reynolds numbers of $270 \leq {\rm Rm} \leq 4000$, to explore the dependence of $k_{\rm p}$ on $k_ν$ and $k_η$. Using physically motivated models for the kinetic and magnetic energy spectra, we measure $k_ν$, $k_η$ and $k_{\rm p}$, making sure that the obtained scales are numerically converged. We determine the overall dissipation scale relations $k_ν= (0.025^{+0.005}_{-0.006})\, k_{\rm turb}\, {\rm Re}^{3/4}$ and $k_η= (0.88^{+0.21}_{-0.23})\, k_ν\, {\rm Pm}^{1/2}$, where $k_{\rm turb}$ is the turbulence driving wavenumber and ${\rm Pm}={\rm Rm}/{\rm Re}$ is the magnetic Prandtl number. We demonstrate that the principle dependence of $k_{\rm p}$ is on $k_η$. For plasmas where ${\rm Re} \gtrsim 100$, we find that $k_{\rm p} = (1.2_{-0.2}^{+0.2})\, k_η$, with the proportionality constant related to the power-law `Kazantsev' exponent of the magnetic power spectrum. Throughout this study, we find a dichotomy in the fundamental properties of the dynamo where ${\rm Re} > 100$, compared to ${\rm Re} < 100$. We report a minimum critical hydrodynamic Reynolds number, ${\rm Re}_{\rm crit} = 100$ for bonafide turbulent dynamo action.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Ion Alfvén velocity fluctuations and implications for the diffusion of streaming cosmic rays

The interstellar medium (ISM) of star-forming galaxies is magnetized and turbulent. Cosmic rays (CRs) propagate through it, and those with energies from $\sim\,\rm{GeV} - \rm{TeV}$ are likely subject to the streaming instability, whereby the wave damping processes balances excitation of resonant ionic Alfvén waves by the CRs, reaching an equilibrium in which the propagation speed of the CRs is very close to the local ion Alfvén velocity. The transport of streaming CRs is therefore sensitive to ionic Alfvén velocity fluctuations. In this paper we systematically study these fluctuations using a large ensemble of compressible MHD turbulence simulations. We show that for sub-Alfvénic turbulence, as applies for a strongly magnetized ISM, the ionic Alfvén velocity probability density function (PDF) is determined solely by the density fluctuations from shocked gas forming parallel to the magnetic field, and we develop analytical models for the ionic Alfvén velocity PDF up to second moments. For super-Alfvénic turbulence, magnetic and density fluctuations are correlated in complex ways, and these correlations as well as contributions from the magnetic fluctuations sets the ionic Alfvén velocity PDF. We discuss the implications of these findings for underlying "macroscopic" diffusion mechanisms in CRs undergoing the streaming instability, including modeling the macroscopic diffusion coefficient for the parallel transport in sub-Alfvénic plasmas. We also describe how, for highly-magnetized turbulent gas, the gas density PDF, and hence column density PDF, can be used to access information about ionic Alfvén velocity structure from observations of the magnetized ISM.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Multiphase turbulence in galactic halos: effect of the driving

Supernova explosions, active galactic nuclei jets, galaxy--galaxy interactions and cluster mergers can drive turbulence in the circumgalactic medium (CGM) and in the intracluster medium (ICM). However, the exact nature of turbulence forced by these sources and its impact on the different statistical properties of the CGM/ICM and their global thermodynamics is still unclear. To investigate the effects of different types of forcing, we conduct high resolution ($1008^3$ resolution elements) idealised hydrodynamic simulations with purely solenoidal (divergence-free) forcing, purely compressive (curl-free) forcing, and natural mixture forcing (equal fractions of the two components). The simulations also include radiative cooling. We study the impact of the three different forcing modes (sol, comp, mix) on the morphology of the gas, its temperature and density distributions, sources and sinks of enstrophy, i.e., solenoidal motions, as well as the kinematics of hot ($\sim10^7~\mathrm{K}$) X-ray emitting and cold ($\sim10^4~\mathrm{K}$) H$α$ emitting gas. We find that compressive forcing leads to stronger variations in density and temperature of the gas as compared to solenoidal forcing. The cold phase gas forms large-scale filamentary structures for compressive forcing and misty, small-scale clouds for solenoidal forcing. The cold phase gas has stronger large-scale velocities for compressive forcing. The natural mixture forcing shows kinematics and gas distributions intermediate between the two extremes, the cold-phase gas occurs as both large-scale filaments and small-scale misty clouds.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The driving mode of shock-driven turbulence

Turbulence in the interstellar medium (ISM) is crucial in the process of star formation. Shocks produced by supernova explosions, jets, radiation from massive stars, or galactic spiral-arm dynamics are amongst the most common drivers of turbulence in the ISM. However, it is not fully understood how shocks drive turbulence, in particular whether shock driving is a more solenoidal(rotational, divergence-free) or a more compressive (potential, curl-free) mode of driving turbulence. The mode of turbulence driving has profound consequences for star formation, with compressive driving producing three times larger density dispersion, and an order of magnitude higher star formation rate than solenoidal driving. Here, we use hydrodynamical simulations of a shock inducing turbulent motions in a structured, multi-phase medium. This is done in the context of a laser-induced shock, propagating into a foam material, in preparation for an experiment to be performed at the National Ignition Facility (NIF). Specifically, we analyse the density and velocity distributions in the shocked turbulent medium, and measure the turbulence driving parameter $b=(σ^{2 Γ}_{ρ/\langle ρ\rangle}-1)^{1/2} (1-σ_{ρ\langle ρ\rangle}^{-2})^{-1/2}\mathcal{M}^{-1}Γ^{-1/2}$ with the density dispersion $σ_{ρ/ \langle ρ\rangle}$, the turbulent Mach number $\mathcal{M}$, and the polytropic exponent $Γ$. Purely solenoidal and purely compressive driving correspond to $b \sim 1/3$ and $b \sim 1$, respectively. Using simulations in which a shock is driven into a multi-phase medium with structures of different sizes and $Γ< 1$, we find $b \sim 1$ for all cases, showing that shock-driven turbulence is consistent with strongly compressive driving.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Velocity Statistics of Turbulent Clouds in the Presence of Gravity, Magnetic fields, Radiation, and Outflow Feedback

The interaction of turbulence, magnetic fields, self-gravity, and stellar feedback within molecular clouds is crucial for understanding star formation. We study the effects of self-gravity and outflow feedback on the properties of the turbulent velocity via the structure function over length scales from $\sim$ 0.01 pc to 2 pc. We analyze a series of three-dimensional, magnetohydrodynamical (MHD) simulations of star cluster formation. We find outflow feedback can change the scaling of velocity fluctuations but still roughly being in between Kolmogorov and Burgers turbulence. We observe that self-gravity and protostellar outflows increase the velocity fluctuations over all length scales. Outflows can amplify the velocity fluctuations by up to a factor of $\sim$7 on scales $\sim$ 0.01 - 0.2 pc and drive turbulence up to a scale of $\sim$ 1 pc. The amplified velocity fluctuations provide more support against gravity and enhance fragmentation on small scales. The self-gravity's effect is more significant on smaller dense clumps and it increases the fraction of the compressive velocity component up to a scale of $\sim$ 0.2 pc. However, outflow feedback drives both solenoidal and compressive modes, but it induces a higher fraction of solenoidal modes relative to compressive modes. Thus, with outflows, the dense core ends up with a slightly higher fraction of solenoidal modes. We find that the compressible fraction is fairly constant with about 1/3 on scales $\sim$ 0.1 - 0.2 pc. The combined effect of enhanced velocity dispersion and reduced compressive fraction contributes to a reduction in the star formation rate.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Theoretical and Observational Evidence for Coriolis Effects in Coronal Magnetic Fields Via Direct Current Driven Flaring Events

All stars produce explosive surface events such as flares and coronal mass ejections. These events are driven by the release of energy stored in coronal magnetic fields, generated by the stellar dynamo. However, it remains unclear if the energy deposition in the magnetic fields is driven by direct or alternating currents. Recently, we presented observational measurements of the flare intensity distributions for a sample of $\sim10^5$ stars across the main sequence observed by $\textit{TESS}$, all of which exhibited power-law distributions similar to those observed in the Sun, albeit with varying slopes. Here we investigate the mechanisms required to produce such a distribution of flaring events via direct current energy deposition, in which coronal magnetic fields braid, reconnect, and produce flares. We adopt a topological model for this process which produces a power-law distribution of energetic flaring events. We expand this model to include the Coriolis effect, which we demonstrate produces a shallower distribution of flare energies in stars that rotate more rapidly (corresponding to a weaker decline in occurrence rates toward increasing flare energies). We present tentative evidence for the predicted rotation-power-law index correlation in the observations. We advocate for future observations of stellar flares that would improve our measurements of the power-law exponents, and yield key insights into the underlying dynamo mechanisms that underpin the self-similar flare intensity distributions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Turbulence Generation by Shock Interaction with a Highly Non-Uniform Medium

An initially planar shock wave propagating into a medium of non-uniform density will be perturbed, leading to the generation of post-shock velocity perturbations. Using numerical simulations we study this phenomenon in the case of highly-non-uniform density (order-unity normalized variance, $σ_ρ/\overlineρ \sim 1$) and strong shocks (shock Mach numbers $\overline{M}_s \gtrsim 10$). This leads to a highly disrupted shock and a turbulent post-shock flow. We simulate this interaction for a range of shock drives and initial density configurations meant to mimic those which might be presently achieved in experiments. Theoretical considerations lead to scaling relations, which are found to reasonably predict the post-shock turbulence properties. The turbulent velocity dispersion and turbulent Mach number are found to depend on the pre-shock density dispersion and shock speed in a manner consistent with the linear Richtymer-Meshkov instability prediction. We also show a dependence of the turbulence generation on the scale of density perturbations. The post-shock pressure and density, which can be substantially reduced relative to the unperturbed case, are found to be reasonably predicted by a simplified analysis that treats the extended shock transition region as a single normal shock.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Turbulent dynamo in the two-phase interstellar medium

Magnetic fields are a dynamically important component of the turbulent interstellar medium (ISM) of star-forming galaxies. These magnetic fields are due to a dynamo action, which is a process of converting turbulent kinetic energy to magnetic energy. A dynamo that acts at scales less than the turbulent driving scale is known as the turbulent dynamo. The ISM is a multiphase medium and observations suggest that the properties of magnetic fields differ with the phase. Here, we aim to study how the properties of the turbulent dynamo depend on the phase. We simulate the non-isothermal turbulent dynamo in a two-phase medium (most previous work assumes an isothermal gas). We show that the warm phase ($T\ge10^3~{\rm K}$) is transsonic and the cold phase ($T<10^3~{\rm K}$) is supersonic. We find that the growth rate of magnetic fields in the exponentially growing stage is similar in both phases. We compute the terms responsible for amplification and destruction of vorticity and show that in both phases vorticity is amplified due to turbulent motions, further amplified by the baroclinic term in the warm phase, and destroyed by the term for viscous interactions in the presence of logarithmic density gradients in the cold phase. We find that the final ratio of magnetic to turbulent kinetic energy is lower in the cold phase due to a stronger Lorentz force. We show that the non-isothermal turbulent dynamo is significantly different from its isothermal counterpart and this demonstrates the need for studying the turbulent dynamo in a multiphase medium.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

VETTAM: A scheme for radiation hydrodynamics with adaptive mesh refinement using the variable Eddington tensor method

We present Variable Eddington Tensor-closed Transport on Adaptive Meshes (\texttt{VETTAM}), a new algorithm to solve the equations of radiation hydrodynamics (RHD) with support for adaptive mesh refinement (AMR) in a frequency-integrated, two-moment formulation. The method is based on a non-local Variable Eddington Tensor (VET) closure computed with a hybrid characteristics scheme for ray tracing. We use a Godunov method for the hyperbolic transport of radiation with an implicit backwards-Euler temporal update to avoid the explicit timestep constraint imposed by the light-crossing time, and a fixed-point Picard iteration scheme to handle the nonlinear gas-radiation exchange term, with the two implicit update stages jointly iterated to convergence. We also develop a modified wave-speed correction method for AMR, which we find to be crucial for obtaining accurate results in the diffusion regime. We demonstrate the robustness of our scheme with a suite of pure radiation and RHD tests, and show that it successfully captures the streaming, static diffusion, and dynamic diffusion regimes and the spatial transitions between them, casts sharp shadows, and yields accurate results for rates of momentum and energy exchange between radiation and gas. A comparison between different closures for the radiation moment equations, with the Eddington approximation (0th-moment closure) and the $M_1$ approximation (1st-moment closure), demonstrates the advantages of the VET method (2nd-moment closure) over the simpler closure schemes. \texttt{VETTAM} has been coupled to the AMR \texttt{FLASH} (magneto-)hydrodynamics code and we summarize by reporting performance features and bottlenecks of our implementation.

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

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Magnetic fields in the Milky Way from pulsar observations: effect of the correlation between thermal electrons and magnetic fields

Pulsars can act as an excellent probe of the Milky Way magnetic field. The average strength of the Galactic magnetic field component parallel to the line of sight can be estimated as $\langle B_\parallel \rangle = 1.232 \, \text{RM}/\text{DM}$, where $\text{RM}$ and $\text{DM}$ are the rotation and dispersion measure of the pulsar. However, this assumes that the thermal electron density and magnetic field of the interstellar medium are uncorrelated. Using numerical simulations and observations, we test the validity of this assumption. Based on magnetohydrodynamical simulations of driven turbulence, we show that the correlation between the thermal electron density and the small-scale magnetic field increases with increasing Mach number of the turbulence. We find that the assumption of uncorrelated thermal electron density and magnetic fields is valid only for subsonic and transsonic flows, but for supersonic turbulence, the field strength can be severely overestimated by using $1.232 \, \text{RM}/\text{DM}$. We then correlate existing pulsar observations from the Australia Telescope National Facility with regions of enhanced thermal electron density and magnetic fields probed by ${^{12} \mathrm {CO}}$ data of molecular clouds, magnetic fields from the Zeeman splitting of the 21 cm line, neutral hydrogen column density, and H$α$ observations. Using these observational data, we show that the thermal electron density and magnetic fields are largely uncorrelated over kpc scales. Thus, we conclude that the relation $\langle B_\parallel \rangle = 1.232 \, \text{RM}/\text{DM}$ provides a good estimate of the magnetic field on Galactic scales but might break down on sub - kpc scales.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Non-ideal MHD simulations of subcritical prestellar cores with non-equilibrium chemistry

Non-ideal magnetohydrodynamic (MHD) effects are thought to be gravity's closest ally in overcoming the support of magnetic fields and in forming stars. Here, we modify the publicly available version of the adaptive mesh refinement code FLASH (Fryxell et al. 2000; Dubey et al. 2008) to include a detailed treatment of non-ideal MHD and study such effects in collapsing prestellar cores. We implement two very extended non-equilibrium chemical networks, the largest of which is comprised of $\sim$ 300 species and includes a detailed description of deuterium chemistry. The ambipolar-diffusion, Ohmic and Hall resistivities are then self-consistently calculated from the abundances of charged species. We present a series of 2-dimensional axisymmetric simulations where we vary the chemical model, cosmic-ray ionization rate, and grain distribution. We benchmark our implementation against ideal MHD simulations and previously-published results. We show that, at high densities ($n_{\rm{H_2}}>~10^6~\rm{cm^{-3}}$), the ion that carries most of the perpendicular and parallel conductivities is not $\rm{H_3^+}$ as was previously thought, but is instead $\rm{D_3^+}$.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Reconstructing three-dimensional densities from two-dimensional observations of molecular gas

Star formation has long been known to be an inefficient process, in the sense that only a small fraction $ε_{\rm ff}$ of the mass of any given gas cloud is converted to stars per cloud free-fall time. However, developing a successful theory of star formation will require measurements of both the mean value of $ε_{\rm ff}$ and its scatter from one molecular cloud to another. Because $ε_{\rm ff}$ is measured relative to the free-fall time, such measurements require accurate determinations of cloud volume densities. Efforts to measure the volume density from two-dimensional projected data, however, have thus far relied on treating molecular clouds as simple uniform spheres, while their real shapes are likely filamentary and their density distributions far from uniform. The resulting uncertainty in the true volume density is likely one of the major sources of error in observational estimates of $ε_{\rm ff}$. In this paper, we use a suite of simulations of turbulent, magnetized, radiative, self-gravitating star-forming clouds to examine whether it is possible to obtain more accurate volume density estimates and thereby reduce this error. We create mock observations from simulations, and show that current analysis methods relying on the spherical assumption likely yield ~ 0.26 dex underestimations and ~ 0.51 dex errors in volume density estimates, corresponding to a ~ 0.13 dex overestimation and a ~ 0.25 dex scatter in $ε_{\rm ff}$, comparable to the scatter in observed cloud samples. We build a predictive model that uses information accessible in two-dimensional measurements -- most significantly the Gini coefficient of the surface density distribution -- to estimate volume density with ~ 0.3 dex less scatter. We test our method on a recent observation of the Ophiuchus cloud, and show that it successfully reduces the $ε_{\rm ff}$ scatter.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The Effects of Magnetic Fields and Outflow Feedback on the Shape and Evolution of the Density PDF in Turbulent Star-Forming Clouds

Using a suite of 3D hydrodynamical simulations of star-forming molecular clouds, we investigate how the density probability distribution function (PDF) changes when including gravity, turbulence, magnetic fields, and protostellar outflows and heating. We find that the density PDF is not lognormal when outflows and self-gravity are considered. Self-gravity produces a power-law tail at high densities and the inclusion of stellar feedback from protostellar outflows and heating produces significant time-varying deviations from a lognormal distribution at the low densities. The simulation with outflows has an excess of diffuse gas compared to the simulations without outflows, exhibits increased average sonic Mach number, and maintains a slower star formation rate over the entire duration of the run. We study the mass transfer between the diffuse gas in the lognormal peak of the PDF, the collapsing gas in the power-law tail, and the stars. We find that the mass fraction in the power-law tail is constant, such that the stars form out of the power-law gas at the same rate at which the gas from the lognormal part replenishes the power-law. We find that turbulence does not provide significant support in the dense gas associated with the power-law tail. When including outflows and magnetic fields in addition to driven turbulence, the rate of mass transfer from the lognormal to the power-law, and then to the stars, becomes significantly slower, resulting in slower star formation rates and longer depletion times.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The physics of gas phase metallicity gradients in galaxies

We present a new model for the evolution of gas phase metallicity gradients in galaxies from first principles. We show that metallicity gradients depend on four ratios that collectively describe the metal equilibration timescale, production, transport, consumption, and loss. Our model finds that most galaxy metallicity gradients are in equilibrium at all redshifts. When normalized by metal diffusion, metallicity gradients are governed by the competition between radial advection, metal production, and accretion of metal-poor gas from the cosmic web. The model naturally explains the varying gradients measured in local spirals, local dwarfs, and high-redshift star-forming galaxies. We use the model to study the cosmic evolution of gradients across redshift, showing that the gradient in Milky Way-like galaxies has steepened over time, in good agreement with both observations and simulations. We also predict the evolution of metallicity gradients with redshift in galaxy samples constructed using both matched stellar masses and matched abundances. Our model shows that massive galaxies transition from the advection-dominated to the accretion-dominated regime from high to low redshifts, which mirrors the transition from gravity-driven to star formation feedback-driven turbulence. Lastly, we show that gradients in local ultraluminous infrared galaxies (major mergers) and inverted gradients seen both in the local and high-redshift galaxies may not be in equilibrium. In subsequent papers in this series, we show that the model also explains the observed relationship between galaxy mass and metallicity gradients, and between metallicity gradients and galaxy kinematics.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Honing and proofing Astrophysical codes on the road to Exascale. Experiences from code modernization on many-core systems

The complexity of modern and upcoming computing architectures poses severe challenges for code developers and application specialists, and forces them to expose the highest possible degree of parallelism, in order to make the best use of the available hardware. The Intel$^{(R)}$ Xeon Phi$^{(TM)}$ of second generation (code-named Knights Landing, henceforth KNL) is the latest many-core system, which implements several interesting hardware features like for example a large number of cores per node (up to 72), the 512 bits-wide vector registers and the high-bandwidth memory. The unique features of KNL make this platform a powerful testbed for modern HPC applications. The performance of codes on KNL is therefore a useful proxy of their readiness for future architectures. In this work we describe the lessons learnt during the optimisation of the widely used codes for computational astrophysics P-Gadget-3, Flash and Echo. Moreover, we present results for the visualisation and analysis tools VisIt and yt. These examples show that modern architectures benefit from code optimisation at different levels, even more than traditional multi-core systems. However, the level of modernisation of typical community codes still needs improvements, for them to fully utilise resources of novel architectures.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Implementation of stellar heating feedback in simulations of star cluster formation: effects on the initial mass function

Explaining the initial mass function (IMF) of stars is a long-standing problem in astrophysics. The number of complex mechanisms involved in the process of star cluster formation, such as turbulence, magnetic fields and stellar feedback, make understanding and modeling the IMF a challenging task. In this paper, we aim to assert the importance of stellar heating feedback in the star cluster formation process and its effect on the shape of the IMF. We use an analytical sub-grid model to implement the radiative feedback in fully three-dimensional magnetohydrodynamical (MHD) simulations of star cluster formation, with the ultimate objective of obtaining numerical convergence on the IMF. We compare a set of MHD adaptive-mesh-refinement (AMR) simulations with three different implementations of the heating of the gas: 1) a polytropic equation of state (EOS), 2) a spherically symmetric stellar heating feedback, and 3) our newly developed polar heating model that takes into account the geometry of the accretion disc and the resulting shielding of stellar radiation by dust. For each of the three heating models, we analyse the distribution of stellar masses formed in ten molecular cloud simulations with different realizations of the turbulence to obtain a statistically representative IMF. We conclude that stellar heating feedback has a profound influence on the number of stars formed and plays a crucial role in controlling the IMF. We find that the simulations with the polar heating model achieve the best convergence on the observed IMF.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Magnetic field fluctuations in anisotropic, supersonic turbulence

The rich structure that we observe in molecular clouds is due to the interplay between strong magnetic fields and supersonic (turbulent) velocity fluctuations. The velocity fluctuations interact with the magnetic field, causing it too to fluctuate. Using numerical simulations, we explore the nature of such magnetic field fluctuations, $\vec{δB}$, over a wide range of turbulent Mach numbers, $\mathcal{M} = 2 - 20$ (i.e., from weak to strong compressibility), and Alfvén Mach numbers, $\mathcal{M}_{\text{A}0} = 0.1 - 100$ (i.e., from strong to weak magnetic mean fields, $B_0$). We derive a compressible quasi-static fluctuation model from the magnetohydrodynamical (MHD) equations and show that velocity gradients parallel to the mean magnetic field give rise to compressible modes in sub-Alfvénic flows, which prevents the flow from becoming two-dimensional, as is the case in incompressible MHD turbulence. We then generalise an analytical model for the magnitude of the magnetic fluctuations to include $\mathcal{M}$, and find $|\vec{δB}| = δB = c_s\sqrt{πρ_0}\mathcal{M}\mathcal{M}_{\text{A}0}$, where $c_s$ is the sound speed and $ρ_0$ is the mean density of gas. This new relation fits well in the strong $B$-field regime. We go on to study the anisotropy between the perpendicular ($ B_{\perp}$) and parallel ($ B_{\parallel}$) fluctuations and the mean-normalised fluctuations, which we find follow universal scaling relations, invariant of $\mathcal{M}$. We provide a detailed analysis of the morphology for the $δB_{\perp}$ and $δB_{\parallel}$ probability density functions and find that eddies aligned with $B_0$ cause parallel fluctuations that reduce $B_{\parallel}$ in the most anisotropic simulations. We discuss broadly the implications of our fluctuation models for magnetised gases in the interstellar medium.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Magnetic fields in elliptical galaxies: an observational probe of the fluctuation dynamo action

Fluctuation dynamos are thought to play an essential role in magnetized galaxy evolution, saturating within $\sim0.01~$Gyr and thus potentially acting as seeds for large-scale dynamos. However, unambiguous observational confirmation of the fluctuation dynamo action in a galactic environment is still missing. This is because, in spiral galaxies, it is difficult to differentiate between small-scale magnetic fields generated by a fluctuation dynamo and those due to the tangling of the large-scale field. We propose that observations of magnetic fields in elliptical galaxies would directly probe the fluctuation dynamo action. This is motivated by the fact that in ellipticals, due to their lack of significant rotation, the conventional large-scale dynamo is absent and the fluctuation dynamo is responsible for controlling the magnetic field strength and structure. By considering turbulence injected by Type Ia supernova explosions and possible magnetic field amplification by cooling flows, we estimate expected magnetic field strengths of $0.2~-~1 \,μ{\rm G}$ in quiescent elliptical galaxies. We use a semi-analytic model of galaxy formation to estimate the distribution and redshift evolution of field strengths, tentatively finding a decrease in magnetic field strength with decreasing redshift. We analyse a sample of radio sources that exhibit the Laing-Garrington (LG) effect (radio polarization asymmetry in jets) and infer magnetic field strengths between $0.14~-~1.33 \,μ{\rm G}$ for a uniform thermal electron density and between $1.36~-~6.21\,μ{\rm G}$ for the thermal electron density following the King profile. We examine observational techniques for measuring the magnetic field saturation state in elliptical galaxies, focusing on Faraday RM grids, the LG effect, synchrotron emission, and gravitational lensing, finding appealing prospects for future empirical analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Molecular cloud formation by compression of magnetized turbulent gas subjected to radiative cooling

Complex turbulent motions of magnetized gas are ubiquitous in the interstellar medium. The source of this turbulence, however, is still poorly understood. Previous work suggests that compression caused by supernova shockwaves, gravity, or cloud collisions, may drive the turbulence to some extent. In this work, we present three-dimensional (3D) magnetohydrodynamic (MHD) simulations of contraction in turbulent, magnetized clouds from the warm neutral medium (WNM) of the ISM to the formation of cold dense molecular clouds, including radiative heating and cooling. We study different contraction rates and find that observed molecular cloud properties, such as the temperature, density, Mach number, and magnetic field strength, and their respective scaling relations, are best reproduced when the contraction rate equals the turbulent turnover rate. In contrast, if the contraction rate is significantly larger (smaller) than the turnover rate, the compression drives too much (too little) turbulence, producing unrealistic cloud properties. We find that the density probability distribution function evolves from a double log-normal representing the two-phase ISM, to a skewed, single log-normal in the dense, cold phase. For purely hydrodynamical simulations, we find that the effective driving parameter of contracting cloud turbulence is natural to mildly compressive (\mbox{$b\sim0.4$--$0.5$}), while for MHD turbulence, we find \mbox{$b\sim0.3$--$0.4$}, i.e., solenoidal to naturally mixed. Overall, the physical properties of the simulated clouds that contract at a rate equal to the turbulent turnover rate, indicate that large-scale contraction may explain the origin and evolution of turbulence in the ISM.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the turbulence driving mode of expanding HII regions

We investigate the turbulence driving mode of ionizing radiation from massive stars on the surrounding interstellar medium (ISM). We run hydrodynamical simulations of a turbulent cloud impinged by a plane-parallel ionization front. We find that the ionizing radiation forms pillars of neutral gas reminiscent of those seen in observations. We quantify the driving mode of the turbulence in the neutral gas by calculating the driving parameter $b$, which is characterised by the relation $σ_s^2 = \ln({1+b^2\mathcal{M}^2})$ between the variance of the logarithmic density contrast $σ_s^2$ (where $s = \ln({ρ/ρ_0})$ with the gas density $ρ$ and its average $ρ_0$), and the turbulent Mach number $\mathcal{M}$. Previous works have shown that $b\sim1/3$ indicates solenoidal (divergence-free) driving and $b\sim1$ indicates compressive (curl-free) driving, with $b\sim1$ producing up to ten times higher star formation rates than $b\sim1/3$. The time variation of $b$ in our study allows us to infer that ionizing radiation is inherently a compressive turbulence driving source, with a time-averaged $b\sim 0.76 \pm 0.08$. We also investigate the value of $b$ of the pillars, where star formation is expected to occur, and find that the pillars are characterised by a natural mixture of both solenoidal and compressive turbulent modes ($b\sim0.4$) when they form, and later evolve into a more compressive turbulent state with $b\sim0.5$--$0.6$. A virial parameter analysis of the pillar regions supports this conclusion. This indicates that ionizing radiation from massive stars may be able to trigger star formation by producing predominately compressive turbulent gas in the pillars.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantifying the effects of spatial resolution and noise on galaxy metallicity gradients

Metallicity gradients are important diagnostics of galaxy evolution, because they record the history of events such as mergers, gas inflow and star-formation. However, the accuracy with which gradients can be measured is limited by spatial resolution and noise, and hence measurements need to be corrected for such effects. We use high resolution (~20 pc) simulation of a face-on Milky Way mass galaxy, coupled with photoionisation models, to produce a suite of synthetic high resolution integral field spectroscopy (IFS) datacubes. We then degrade the datacubes, with a range of realistic models for spatial resolution (2 to 16 beams per galaxy scale length) and noise, to investigate and quantify how well the input metallicity gradient can be recovered as a function of resolution and signal-to-noise ratio (SNR) with the intention to compare with modern IFS surveys like MaNGA and SAMI. Given appropriate propagation of uncertainties and pruning of low SNR pixels, we show that a resolution of 3-4 telescope beams per galaxy scale length is sufficient to recover the gradient to ~10-20% uncertainty. The uncertainty escalates to ~60% for lower resolution. Inclusion of the low SNR pixels causes the uncertainty in the inferred gradient to deteriorate. Our results can potentially inform future IFS surveys regarding the resolution and SNR required to achieve a desired accuracy in metallicity gradient measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Spinning Bar and a Star-formation Inefficient Repertoire: Turbulence in Hickson Compact Group NGC7674

The physics regulating star formation (SF) in Hickson Compact Groups (HCG) has thus far been difficult to describe, due to their unique kinematic properties. In this study we expand upon previous works to devise a more physically meaningful SF relation able to better encompass the physics of these unique systems. We combine CO(1--0) data from the Combined Array from Research in Millimeter Astronomy (CARMA) to trace the column density of molecular gas $Σ_\mathrm{gas}$ and deep H$α$ imaging taken on the Southern Astrophysical Research (SOAR) Telescope tracing $Σ_\mathrm{SFR}$ to investigate star formation efficiency across face-on HCG, NGC7674. We find a lack of universality in star formation, with two distinct sequences present in the $Σ_\mathrm{gas}-Σ_\mathrm{SFR}$ plane; one for inside and one for outside the nucleus. We devise a SF relation based on the multi-freefall nature of gas and the critical density, which itself is dependent on the virial parameter $α_{\mathrm{vir}}$, the ratio of turbulent to gravitational energy. We find that our modified SF relation fits the data and describes the physics of this system well with the introduction of a virial parameter of about 5--10 across the galaxy. This $α_{\mathrm{vir}}$ leads to an order-of-magnitude reduction in SFR compared to $α_{\mathrm{vir}}\approx 1$ systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Star formation in cloud cores -- simulations and observations of dense molecular cores and the formation of solar mass stars

Star formation is inefficient. Recent advances in numerical simulations and theoretical models of molecular clouds show that the combined effects of interstellar turbulence, magnetic fields and stellar feedback can explain the low efficiency of star formation. The star formation rate is highly sensitive to the driving mode of the turbulence. Solenoidal driving may be more important in the Central Molecular Zone, compared to more compressive driving agents in spiral-am clouds. Both theoretical and observational efforts are underway to determine the dominant driving mode of turbulence in different Galactic environments. New observations with ALMA, combined with other instruments such as CARMA, JCMT and the SMA begin to reveal the magnetic field structure of dense cores and protostellar disks, showing highly complex field geometries with ordered and turbulent field components. Such complex magnetic fields can give rise to a range of stellar masses and jet/outflow efficiencies in dense cores and protostellar accretion disks.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The importance of magnetic fields for the initial mass function of the first stars

Magnetic fields play an important role for the formation of stars in both local and high-redshift galaxies. Recent studies of dynamo amplification in the first dark matter haloes suggest that significant magnetic fields were likely present during the formation of the first stars in the Universe at redshifts of 15 and above. In this work, we study how these magnetic fields potentially impact the initial mass function (IMF) of the first stars. We perform 200 high-resolution, three-dimensional (3D), magneto-hydrodynamic (MHD) simulations of the collapse of primordial clouds with different initial turbulent magnetic field strengths as predicted from turbulent dynamo theory in the early Universe, forming more than 1100 first stars in total. We detect a strong statistical signature of suppressed fragmentation in the presence of strong magnetic fields, leading to a dramatic reduction in the number of first stars with masses low enough that they might be expected to survive to the present day. Additionally, strong fields shift the transition point where stars go from being mostly single to mostly multiple to higher masses. However, irrespective of the field strength, individual simulations are highly chaotic, show different levels of fragmentation and clustering, and the outcome depends on the exact realisation of the turbulence in the primordial clouds. While these are still idealised simulations that do not start from cosmological initial conditions, our work shows that magnetic fields play a key role for the primordial IMF, potentially even more so than for the present-day IMF.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Launching of Cold Clouds by Galaxy Outflows III: The Influence of Magnetic Fields

Motivated by observations of outflowing galaxies, we investigate the combined impact of magnetic fields and radiative cooling on the evolution of cold clouds embedded in a hot wind. We perform a collection of three-dimensional adaptive mesh refinement, magnetohydrodynamical simulations that span two resolutions, and include fields that are aligned and transverse to the oncoming, super-Alfvénic material. Aligned fields have little impact on the overall lifetime of the clouds over the non-magnetized case, although they do increase the mixing between the wind and cloud material by a factor of $\approx 3.$ Transverse fields lead to magnetic draping, which isolates the clouds, but they also squeeze material in the direction perpendicular to the field lines, which leads to rapid mass loss. A resolution study suggests that the magnetized simulations have somewhat better convergence properties than non-magnetized simulations, and that a resolution of 64 zones per cloud radius is sufficient to accurately describe these interactions. We conclude that the combined effects of radiative cooling and magnetic fields are dependent on field orientation, but are unlikely to enhance cloud lifetimes beyond the effect of radiative cooling alone.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Predicted Properties of Helium-Enriched Globular Cluster Progenitors at High Redshift

Globular cluster progenitors may have been detected by \textit{HST}, and are predicted to be observable with \textit{JWST} and ground-based extremely-large telescopes with adaptive optics. This has the potential to elucidate the issue of globular cluster formation and the origins of significantly helium-enriched subpopulations, a problem in Galactic astronomy with no satisfactory theoretical solution. Given this context, we use model stellar tracks and isochrones to investigate the predicted observational properties of helium-enriched stellar populations in globular cluster progenitors. We find that, relative to helium-normal populations, helium-enriched ($ΔY=+0.12$) stellar populations similar to those inferred in the most massive globular clusters, are expected, modulo some rapid fluctuations in the first $\sim$30 Myr, to be brighter and redder in the rest frame. At fixed age, stellar mass, and metallicity, a helium-enriched population is predicted to converge to being $\sim$0.40 mag brighter at $λ\approx 2.0\, μm$, and to be 0.30 mag redder in the \textit{JWST}-NIRCam colour $(F070W-F200W)$, and to actually be fainter for $λ\lesssim 0.50 \, μm$. Separately, we find that the time-integrated shift in ionizing radiation is a negligible $\sim 5\%$, though we show that the Lyman-$α$ escape fraction could end up higher for helium-enriched stars.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The SAMI Galaxy Survey: Gas velocity dispersions in low-$z$ star-forming galaxies and the drivers of turbulence

We infer the intrinsic ionised gas kinematics for 383 star-forming galaxies across a range of integrated star-formation rates (SFR $\in [10^{-3}, 10^2]$ M$_\odot$ yr$^{-1}$) at $z \lesssim 0.1$ using a consistent 3D forward-modelling technique. The total sample is a combination of galaxies from the SAMI Galaxy Survey and DYNAMO survey. For typical low-$z$ galaxies taken from the SAMI Galaxy Survey, we find the vertical velocity dispersion ($σ_{v, z}$) to be positively correlated with measures of star-formation rate, stellar mass, HI gas mass, and rotational velocity. The greatest correlation is with star-formation rate surface density ($Σ_\text{SFR}$). Using the total sample, we find $σ_{v, z}$ increases slowly as a function of integrated star-formation rate in the range SFR $\in$ [$10^{-3}$, 1] M$_\odot$ yr$^{-1}$ from $17\pm3$ km s$^{-1}$ to $24\pm5$ km s$^{-1}$ followed by a steeper increase up to $σ_{v, z}$ $\sim 80$ km s$^{-1}$ for SFR $\gtrsim 1$ M$_\odot$ yr$^{-1}$. This is consistent with recent theoretical models that suggest a $σ_{v, z}$ floor driven by star-formation feedback processes with an upturn in $σ_{v, z}$ at higher SFR driven by gravitational transport of gas through the disc.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The sonic scale revealed by the world's largest supersonic turbulence simulation

Understanding the physics of turbulence is crucial for many applications, including weather, industry, and astrophysics. In the interstellar medium (ISM), supersonic turbulence plays a crucial role in controlling the gas density and velocity structure, and ultimately the birth of stars. Here we present a simulation of interstellar turbulence with a grid resolution of 10048^3 cells that allows us to determine the position and width of the sonic scale (l_s) - the transition from supersonic to subsonic turbulence. The simulation simultaneously resolves the supersonic and subsonic cascade, v(l) ~ l^p, where we measure p_sup = 0.49 +/- 0.01 and p_sub = 0.39 +/- 0.02, respectively. We find that l_s agrees with the relation l_s / L = phi_s Mach^(-1/p_sup), where Mach is the three-dimensional Mach number, and L is either the driving scale of turbulence or the diameter of a molecular cloud. If L is the driving scale, we measure phi_s = 0.42 (+0.12) (-0.09), primarily because of the separation between the driving scale and the start of the supersonic cascade. For a supersonic cascade extending beyond the cloud scale, we get phi_s = 0.91 (+0.25) (-0.20). In both cases, phi_s < 1, because we find that the supersonic cascade transitions smoothly to the subsonic cascade over a factor of 3 in scale, instead of a sharp transition. Our measurements provide quantitative input for turbulence-regulated models of filament structure and star formation in molecular clouds.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Turbulence in stratified atmospheres: implications for the intracluster medium

The gas motions in the intracluster medium (ICM) are governed by stratified turbulence. Stratified turbulence is fundamentally different from Kolmogorov (isotropic, homogeneous) turbulence; kinetic energy not only cascades from large to small scales, but it is also converted into buoyancy potential energy. To understand the density and velocity fluctuations in the ICM, we conduct high-resolution ($1024^2\times 1536$ grid points) hydrodynamical simulations of subsonic turbulence (with rms Mach number $\mathcal{M}\approx 0.25$) and different levels of stratification, quantified by the Richardson number $\mathrm{Ri}$, from $\mathrm{Ri}=0$ (no stratification) to $\mathrm{Ri}=13$ (strong stratification). We quantify the density, pressure and velocity fields for varying stratification because observational studies often use surface brightness fluctuations to infer the turbulent gas velocities of the ICM. We find that the standard deviation of the logarithmic density fluctuations ($σ_s$), where $s=\ln(ρ/\left<ρ(z)\right>)$, increases with $\mathrm{Ri}$. For weakly stratified subsonic turbulence ($\mathrm{Ri}\lesssim10$, $\mathcal{M}<1$), we derive a new $σ_s$--$\mathcal{M}$--$\mathrm{Ri}$ relation, $σ_s^2=\ln(1+b^2\mathcal{M}^4+0.09\mathcal{M}^2\mathrm{Ri}H_P/H_S)$, where $b=1/3$--$1$ is the turbulence driving parameter, and $H_P$ and $H_S$ are the pressure and entropy scale heights respectively. We further find that the power spectrum of density fluctuations, $P(ρ_k/\left<ρ\right>)$, increases in magnitude with increasing $\mathrm{Ri}$, whereas the velocity power spectrum is invariant. Thus, the ratio between density and velocity power spectra strongly depends on $\mathrm{Ri}$. Pressure fluctuations, on the other hand, are independent of stratification and only depend on $\mathcal{M}$.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Filaments and striations: anisotropies in observed, supersonic, highly-magnetised turbulent clouds

Stars form in highly-magnetised, supersonic turbulent molecular clouds. Many of the tools and models that we use to carry out star formation studies rely upon the assumption of cloud isotropy. However, structures like high-density filaments in the presence of magnetic fields, and magnetosonic striations introduce anisotropies into the cloud. In this study we use the two-dimensional (2D) power spectrum to perform a systematic analysis of the anisotropies in the column density for a range of Alfvén Mach numbers ($\mathcal{M}_A=0.1$--$10$) and turbulent Mach numbers ($\mathcal{M}=2$--$20$), with 20 high-resolution, three-dimensional (3D) turbulent magnetohydrodynamic simulations. We find that for cases with a strong magnetic guide field, corresponding to $\mathcal{M}_A<1$, and $\mathcal{M}\lesssim 4$, the anisotropy in the column density is dominated by thin striations aligned with the magnetic field, while for $\mathcal{M}\gtrsim 4$ the anisotropy is significantly changed by high-density filaments that form perpendicular to the magnetic guide field. Indeed, the strength of the magnetic field controls the degree of anisotropy and whether or not any anisotropy is present, but it is the turbulent motions controlled by $\mathcal{M}$ that determine which kind of anisotropy dominates the morphology of a cloud.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A comparison between grid and particle methods on the small-scale dynamo in magnetised supersonic turbulence

We perform a comparison between the smoothed particle magnetohydrodynamics (SPMHD) code, Phantom, and the Eulerian grid-based code, Flash, on the small-scale turbulent dynamo in driven, Mach 10 turbulence. We show, for the first time, that the exponential growth and saturation of an initially weak magnetic field via the small-scale dynamo can be successfully reproduced with SPMHD. The two codes agree on the behaviour of the magnetic energy spectra, the saturation level of magnetic energy, and the distribution of magnetic field strengths during the growth and saturation phases. The main difference is that the dynamo growth rate, and its dependence on resolution, differs between the codes, caused by differences in the numerical dissipation and shock capturing schemes leading to differences in the effective Prandtl number in Phantom and Flash.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Extreme Scale-out SuperMUC Phase 2 - lessons learned

In spring 2015, the Leibniz Supercomputing Centre (Leibniz-Rechenzentrum, LRZ), installed their new Peta-Scale System SuperMUC Phase2. Selected users were invited for a 28 day extreme scale-out block operation during which they were allowed to use the full system for their applications. The following projects participated in the extreme scale-out workshop: BQCD (Quantum Physics), SeisSol (Geophysics, Seismics), GPI-2/GASPI (Toolkit for HPC), Seven-League Hydro (Astrophysics), ILBDC (Lattice Boltzmann CFD), Iphigenie (Molecular Dynamic), FLASH (Astrophysics), GADGET (Cosmological Dynamics), PSC (Plasma Physics), waLBerla (Lattice Boltzmann CFD), Musubi (Lattice Boltzmann CFD), Vertex3D (Stellar Astrophysics), CIAO (Combustion CFD), and LS1-Mardyn (Material Science). The projects were allowed to use the machine exclusively during the 28 day period, which corresponds to a total of 63.4 million core-hours, of which 43.8 million core-hours were used by the applications, resulting in a utilization of 69%. The top 3 users were using 15.2, 6.4, and 4.7 million core-hours, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Magnetic field amplification in turbulent astrophysical plasmas

Magnetic fields play an important role in astrophysical accretion discs, and in the interstellar and intergalactic medium. They drive jets, suppress fragmentation in star-forming clouds and can have a significant impact on the accretion rate of stars. However, the exact amplification mechanisms of cosmic magnetic fields remain relatively poorly understood. Here I start by reviewing recent advances in the numerical and theoretical modelling of the 'turbulent dynamo', which may explain the origin of galactic and inter-galactic magnetic fields. While dynamo action was previously investigated in great detail for incompressible plasmas, I here place particular emphasis on highly compressible astrophysical plasmas, which are characterised by strong density fluctuations and shocks, such as the interstellar medium. I find that dynamo action works not only in subsonic plasmas, but also in highly supersonic, compressible plasmas, as well as for low and high magnetic Prandtl numbers. I further present new numerical simulations from which I determine the growth of the turbulent (un-ordered) magnetic field component ($B_\mathrm{turb}$) in the presence of weak and strong guide fields ($B_0$). I vary $B_0$ over 5 orders of magnitude and find that the dependence of $B_\mathrm{turb}$ on $B_0$ is relatively weak, and can be explained with a simple theoretical model in which the turbulence provides the energy to amplify $B_\mathrm{turb}$. Finally, I discuss some important implications of magnetic fields for the structure of accretion discs, the launching of jets, and the star formation rate of interstellar clouds.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the universality of interstellar filaments: theory meets simulations and observations

Filaments are ubiquitous in the universe. Recent observations have revealed that stars and star clusters form preferentially along dense filaments. Understanding the formation and properties of filaments is therefore a crucial step in understanding star formation. Here we perform three-dimensional high-resolution magnetohydrodynamical simulations that follow the evolution of molecular clouds and the formation of filaments and stars. We apply a filament detection algorithm and compare simulations with different combinations of physical ingredients: gravity, turbulence, magnetic fields and jet/outflow feedback. We find that gravity-only simulations produce significantly narrower filament profiles than observed, while simulations that include turbulence produce realistic filament properties. For these turbulence simulations, we find a remarkably universal filament width of 0.10 +/- 0.02 pc, which is independent of the star formation history of the clouds. We derive a theoretical model that provides a physical explanation for this characteristic filament width, based on the sonic scale (lambda_sonic) of molecular cloud turbulence. Our derivation provides lambda_sonic as a function of the cloud diameter L, the velocity dispersion sigma_v, the gas sound speed c_s, and the ratio of thermal to magnetic pressure, plasma beta. For typical cloud conditions in the Milky Way spiral arms, we find lambda_sonic = 0.04-0.16 pc, in excellent agreement with the filament width of 0.05-0.15 pc from observations. Consistent with the theoretical model assumptions, we find that the velocity dispersion inside the filaments is subsonic and supersonic outside. We further explain the observed p=2 scaling of the filament density profile, rho ~ r^(-p) with the collision of two planar shocks forming a filament at their intersection.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

SHOCKFIND - An algorithm to identify magnetohydrodynamic shock waves in turbulent clouds

The formation of stars occurs in the dense molecular cloud phase of the interstellar medium. Observations and numerical simulations of molecular clouds have shown that supersonic magnetised turbulence plays a key role for the formation of stars. Simulations have also shown that a large fraction of the turbulent energy dissipates in shock waves. The three families of MHD shocks --- fast, intermediate and slow --- distinctly compress and heat up the molecular gas, and so provide an important probe of the physical conditions within a turbulent cloud. Here we introduce the publicly available algorithm, SHOCKFIND, to extract and characterise the mixture of shock families in MHD turbulence. The algorithm is applied to a 3-dimensional simulation of a magnetised turbulent molecular cloud, and we find that both fast and slow MHD shocks are present in the simulation. We give the first prediction of the mixture of turbulence-driven MHD shock families in this molecular cloud, and present their distinct distributions of sonic and Alfvenic Mach numbers. Using subgrid one-dimensional models of MHD shocks we estimate that ~0.03 % of the volume of a typical molecular cloud in the Milky Way will be shock heated above 50 K, at any time during the lifetime of the cloud. We discuss the impact of this shock heating on the dynamical evolution of molecular clouds.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The reliability of observational measurements of column density probability distribution functions

Probability distribution functions (PDFs) of column densities are an established tool to characterize the evolutionary state of interstellar clouds. Using simulations, we show to what degree their determination is affected by noise, line-of-sight contamination, field selection, and the incomplete sampling in interferometric measurements. We solve the integrals that describe the convolution of a cloud PDF with contaminating sources and study the impact of missing information on the measured column density PDF. The effect of observational noise can be easily estimated and corrected for if the root mean square (rms) of the noise is known. For $σ_{noise}$ values below 40% of the typical cloud column density, $N_{peak}$, this involves almost no degradation of the accuracy of the PDF parameters. For higher noise levels and narrow cloud PDFs the width of the PDF becomes increasingly uncertain. A contamination by turbulent foreground or background clouds can be removed as a constant shield if the PDF of the contamination peaks at a lower column or is narrower than that of the observed cloud. Uncertainties in the definition of the cloud boundaries mainly affect the low-column density part of the PDF and the mean density. As long as more than 50% of a cloud are covered, the impact on the PDF parameters is negligible. In contrast, the incomplete sampling of the uv plane in interferometric observations leads to uncorrectable distortions of the PDF of the produced maps. An extension of ALMA's capabilities would allow us to recover the high-column density tail of the PDF but we found no way to measure the intermediate and low column density part of the underlying cloud PDF in interferometric observations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The small-scale turbulent dynamo in smoothed particle magnetohydrodynamics

Supersonic turbulence is believed to be at the heart of star formation. We have performed smoothed particle magnetohydrodynamics (SPMHD) simulations of the small-scale dynamo amplification of magnetic fields in supersonic turbulence. The calculations use isothermal gas driven at rms velocity of Mach 10 so that conditions are representative of star-forming molecular clouds in the Milky Way. The growth of magnetic energy is followed for 10 orders in magnitude until it reaches saturation, a few percent of the kinetic energy. The results of our dynamo calculations are compared with results from grid-based methods, finding excellent agreement on their statistics and their qualitative behaviour. The simulations utilise the latest algorithmic developments we have developed, in particular, a new divergence cleaning approach to maintain the solenoidal constraint on the magnetic field and a method to reduce the numerical dissipation of the magnetic shock capturing scheme. We demonstrate that our divergence cleaning method may be used to achieve $\nabla \cdot {\bf B}=0$ to machine precision, albeit at significant computational expense.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The world's largest turbulence simulations

Understanding turbulence is critical for a wide range of terrestrial and astrophysical applications. Here we present first results of the world's highest-resolution simulation of turbulence ever done. The current simulation has a grid resolution of 10048^3 points and was performed on 65536 compute cores on SuperMUC at the Leibniz Supercomputing Centre (LRZ). We present a scaling test of our modified version of the FLASH code, which updates the hydrodynamical equations in less than 3 micro seconds per cell per time step. A first look at the column density structure of the 10048^3 simulation is presented and a detailed analysis is provided in a forthcoming paper.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A universal, turbulence-regulated star formation law: from Milky Way clouds to high-redshift disk and starburst galaxies

Whilst the star formation rate (SFR) of molecular clouds and galaxies is key in understanding galaxy evolution, the physical processes which determine the SFR remain unclear. This uncertainty about the underlying physics has resulted in various different star formation laws, all having substantial intrinsic scatter. Extending upon previous works that define the column density of star formation (Sigma_SFR) by the gas column density (Sigma_gas), we develop a new universal star formation (SF) law based on the multi-freefall prescription of gas. This new SF law relies predominantly on the probability density function (PDF) and on the sonic Mach number of the turbulence in the star-forming clouds. By doing so we derive a relation where the star formation rate (SFR) correlates with the molecular gas mass per multi-freefall time, whereas previous models had used the average, single-freefall time. We define a new quantity called maximum (multi-freefall) gas consumption rate (MGCR) and show that the actual SFR is only about 0.4% of this maximum possible SFR, confirming the observed low efficiency of star formation. We show that placing observations in this new framework (Sigma_SFR vs. MGCR) yields a significantly improved correlation with 3-4 times reduced scatter compared to previous SF laws and a goodness-of-fit parameter R^2=0.97. By inverting our new relationship, we provide sonic Mach number predictions for kpc-scale observations of Local Group galaxies as well as unresolved observations of local and high-redshift disk and starburst galaxies that do not have independent, reliable estimates for the turbulent cloud Mach number.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Filament formation in wind-cloud interactions. I. Spherical clouds in uniform magnetic fields

Filamentary structures are ubiquitous in the interstellar medium, yet their formation, internal structure, and longevity have not been studied in detail. We report the results from a comprehensive numerical study that investigates the characteristics, formation, and evolution of filaments arising from magnetohydrodynamic interactions between supersonic winds and dense clouds. Here we improve on previous simulations by utilising sharper density contrasts and higher numerical resolutions. By following multiple density tracers, we find that material in the envelopes of the clouds is removed and deposited downstream to form filamentary tails, while the cores of the clouds serve as footpoints and late-stage outer layers of these tails. Aspect ratios >12, subsonic velocity dispersions ~0.1-0.3 of the wind sound speed, and magnetic field amplifications ~100 are found to be characteristic of these filaments. We also report the effects of different magnetic field strengths and orientations. The magnetic field strength regulates vorticity production: sinuous filamentary towers arise in non-magnetic environments, while strong magnetic fields inhibit small-scale perturbations at boundary layers making tails less turbulent. Magnetic field components aligned with the direction of the flow favour the formation of pressure-confined flux ropes inside the tails, whilst transverse components tend to form current sheets. Softening the equation of state to nearly isothermal leads to suppression of dynamical instabilities and further collimation of the tail. Towards the final stages of the evolution, we find that small cloudlets and distorted filaments survive the break-up of the clouds and become entrained in the winds, reaching velocities ~0.1 of the wind speed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Inefficient star formation through turbulence, magnetic fields and feedback

Star formation is inefficient. Only a few percent of the available gas in molecular clouds forms stars, leading to the observed low star formation rate (SFR). The same holds when averaged over many molecular clouds, such that the SFR of whole galaxies is again surprisingly low. Indeed, considering the low temperatures, molecular clouds should be highly gravitationally unstable and collapse on their global mean freefall timescale. And yet, they are observed to live about 10-100 times longer, i.e., the SFR per freefall time (SFR_ff) is only a few percent. Thus, other physical mechanisms must counteract the quick global collapse. Turbulence, magnetic fields and stellar feedback have been proposed as regulating agents, but it is still unclear which of these processes is the most important and what their relative contributions are. Here we run high-resolution simulations including gravity, turbulence, magnetic fields, and jet/outflow feedback. We confirm that clouds collapse on a mean freefall time, if only gravity is considered, producing stars at an unrealistic rate. In contrast, if turbulence, magnetic fields, and feedback are included step-by-step, the SFR is reduced by a factor of 2-3 with each additional physical ingredient. When they all act in concert, we find a constant SFR_ff = 0.04, currently the closest match to observations, but still about a factor of 2-4 higher than the average. A detailed comparison with other simulations and with observations leads us to conclude that only models with turbulence producing large virial parameters, and including magnetic fields and feedback can produce realistic SFRs.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Nonlinear closures for scale separation in supersonic magnetohydrodynamic turbulence

Turbulence in compressible plasma plays a key role in many areas of astrophysics and engineering. The extreme plasma parameters in these environments, e.g. high Reynolds numbers, supersonic and super-Alfvenic flows, however, make direct numerical simulations computationally intractable even for the simplest treatment -- magnetohydrodynamics (MHD). To overcome this problem one can use subgrid-scale (SGS) closures -- models for the influence of unresolved, subgrid-scales on the resolved ones. In this work we propose and validate a set of constant coefficient closures for the resolved, compressible, ideal MHD equations. The subgrid-scale energies are modeled by Smagorinsky-like equilibrium closures. The turbulent stresses and the electromotive force (EMF) are described by expressions that are nonlinear in terms of large scale velocity and magnetic field gradients. To verify the closures we conduct a priori tests over 137 simulation snapshots from two different codes with varying ratios of thermal to magnetic pressure ($β_\mathrm{p} = 0.25, 1, 2.5, 5, 25$) and sonic Mach numbers ($M_s = 2, 2.5, 4$). Furthermore, we make a comparison to traditional, phenomenological eddy-viscosity and $α-β-γ$ closures. We find only mediocre performance of the kinetic eddy-viscosity and $α-β-γ$ closures, and that the magnetic eddy-viscosity closure is poorly correlated with the simulation data. Moreover, three of five coefficients of the traditional closures exhibit a significant spread in values. In contrast, our new closures demonstrate consistently high correlation and constant coefficient values over time and and over the wide range of parameters tested. Important aspects in compressible MHD turbulence such as the bi-directional energy cascade, turbulent magnetic pressure and proper alignment of the EMF are well described by our new closures.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Saturation of the Turbulent Dynamo

The origin of strong magnetic fields in the Universe can be explained by amplifying weak seed fields via turbulent motions on small spatial scales and subsequently transporting the magnetic energy to larger scales. This process is known as the turbulent dynamo and depends on the properties of turbulence, i.e. on the hydrodynamical Reynolds number and the compressibility of the gas, and on the magnetic diffusivity. While we know the growth rate the magnetic energy in the linear regime, the saturation level, i.e. the ratio of magnetic energy to turbulent kinetic energy that can be reached, is not known from analytical calculations. In this paper we present the first scale-dependent saturation model based on an effective turbulent resistivity which is determined by the turnover timescale of turbulent eddies and the magnetic energy density. The magnetic resistivity increases compared to the Spitzer value and the effective scale on which the magnetic energy spectrum is at its maximum moves to larger spatial scales. This process ends when the peak reaches a characteristic wavenumber k* which is determined by the critical magnetic Reynolds number. The saturation level of the dynamo also depends on the type of turbulence and differs for the limits of large and small magnetic Prandtl numbers Pm. With our model we find saturation levels between 43.8% and 1.3% for Pm>>1 and between 2.43% and 0.135% for Pm<<1, where the higher values refer to incompressible turbulence and the lower ones to highly compressible turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The density structure and star formation rate of non-isothermal polytropic turbulence

The interstellar medium of galaxies is governed by supersonic turbulence, which likely controls the star formation rate (SFR) and the initial mass function (IMF). Interstellar turbulence is non-universal, with a wide range of Mach numbers, magnetic fields strengths, and driving mechanisms. Although some of these parameters were explored, most previous works assumed that the gas is isothermal. However, we know that cold molecular clouds form out of the warm atomic medium, with the gas passing through chemical and thermodynamic phases that are not isothermal. Here we determine the role of temperature variations by modelling non-isothermal turbulence with a polytropic equation of state (EOS), where pressure and temperature are functions of gas density, P~rho^Gamma, T~rho^(Gamma-1). We use grid resolutions of 2048^3 cells and compare polytropic exponents Gamma=0.7 (soft EOS), Gamma=1 (isothermal EOS), and Gamma=5/3 (stiff EOS). We find a complex network of non-isothermal filaments with more small-scale fragmentation occurring for Gamma<1, while Gamma>1 smoothes out density contrasts. The density probability distribution function (PDF) is significantly affected by temperature variations, with a power-law tail developing at low densities for Gamma>1. In contrast, the PDF becomes closer to a lognormal distribution for Gamma<=1. We derive and test a new density variance - Mach number relation that takes Gamma into account. This new relation is relevant for theoretical models of the SFR and IMF, because it determines the dense gas mass fraction of a cloud, from which stars form. We derive the SFR as a function of Gamma and find that it decreases by a factor of ~5 from Gamma=0.7 to Gamma=5/3.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The density variance - Mach number relation in isothermal and non-isothermal adiabatic turbulence

The density variance - Mach number relation of the turbulent interstellar medium is relevant for theoretical models of the star formation rate, efficiency, and the initial mass function of stars. Here we use high-resolution hydrodynamical simulations with grid resolutions of up to 1024^3 cells to model compressible turbulence in a regime similar to the observed interstellar medium. We use Fyris Alpha, a shock-capturing code employing a high-order Godunov scheme to track large density variations induced by shocks. We investigate the robustness of the standard relation between the logarithmic density variance (sigma_s^2) and the sonic Mach number (M) of isothermal interstellar turbulence, in the non-isothermal regime. Specifically, we test ideal gases with diatomic molecular (gamma = 7/5) and monatomic (gamma = 5/3) adiabatic indices. A periodic cube of gas is stirred with purely solenoidal forcing at low wavenumbers, leading to a fully-developed turbulent medium. We find that as the gas heats in adiabatic compressions, it evolves along the relationship in the density variance - Mach number plane, but deviates significantly from the standard expression for isothermal gases. Our main result is a new density variance - Mach number relation that takes the adiabatic index into account: sigma_s^2 = ln {1+b^2*M^[(5*gamma+1)/3]} and provides good fits for b*M <= 1. A theoretical model based on the Rankine-Hugoniot shock jump conditions is derived, sigma_s^2 = ln {1+(gamma+1)*b^2*M^2/[(gamma-1)*b^2*M^2+2]}, and provides good fits also for b*M > 1. We conclude that this new relation for adiabatic turbulence may introduce important corrections to the standard relation, if the gas is not isothermal.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Collective outflow from a small multiple stellar system

The formation of high-mass stars is usually accompanied by powerful protostellar outflows. Such high-mass outflows are not simply scaled-up versions of their lower-mass counterparts, since observations suggest that the collimation degree degrades with stellar mass. Theoretically, the origins of massive outflows remain open to question because radiative feedback and fragmentation of the accretion flow around the most massive stars, with M > 15 M_Sun, may impede the driving of magnetic disk winds. We here present a three-dimensional simulation of the early stages of core fragmentation and massive star formation that includes a subgrid-scale model for protostellar outflows. We find that stars that form in a common accretion flow tend to have aligned outflow axes, so that the individual jets of multiple stars can combine to form a collective outflow. We compare our simulation to observations with synthetic H_2 and CO observations and find that the morphology and kinematics of such a collective outflow resembles some observed massive outflows, such as Cepheus A and DR 21. We finally compare physical quantities derived from simulated observations of our models to the actual values in the models to examine the reliability of standard methods for deriving physical quantities, demonstrating that those methods indeed recover the actual values to within a factor of 2-3.

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

Principal Component Analysis of Molecular Clouds: Can CO reveal the dynamics?

We use Principal Component Analysis (PCA) to study the gas dynamics in numerical simulations of typical MCs. Our simulations account for the non-isothermal nature of the gas and include a simplified treatment of the time-dependent gas chemistry. We model the CO line emission in a post-processing step using a 3D radiative transfer code. We consider mean number densities n_0 = 30, 100, 300 cm^{-3} that span the range of values typical for MCs in the solar neighbourhood and investigate the slope α_{PCA} of the pseudo structure function computed by PCA for several components: the total density, H2 density, 12CO density, 12CO J = 1 -> 0 intensity and 13CO J = 1 -> 0 intensity. We estimate power-law indices α_{PCA} for different chemical species that range from 0.5 to 0.9, in good agreement with observations, and demonstrate that optical depth effects can influence the PCA. We show that when the PCA succeeds, the combination of chemical inhomogeneity and radiative transfer effects can influence the observed PCA slopes by as much as ~ +/- 0.1. The method can fail if the CO distribution is very intermittent, e.g. in low-density clouds where CO is confined to small fragments.

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Unfolding the Laws of Star Formation: The Density Distribution of Molecular Clouds

The formation of stars shapes the structure and evolution of entire galaxies. The rate and efficiency of this process are affected substantially by the density structure of the individual molecular clouds in which stars form. The most fundamental measure of this structure is the probability density function of volume densities (rho-PDF), which determines the star formation rates predicted with analytical models. This function has remained unconstrained by observations. We have developed an approach to quantify rho-PDFs and establish their relation to star formation. The rho-PDFs instigate a density threshold of star formation and allow us to quantify the star formation efficiency above it. The rho-PDFs provide new constraints for star formation theories and correctly predict several key properties of the star-forming interstellar medium.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A measurement of the turbulence-driven density distribution in a non-star-forming molecular cloud

Molecular clouds are supersonically turbulent. This turbulence governs the initial mass function and the star formation rate. In order to understand the details of star formation, it is therefore essential to understand the properties of turbulence, in particular the probability distribution of density in turbulent clouds. We present formaldehyde volume density measurements of a non-star-forming cloud along the line of sight towards W49A. We use these measurements in conjunction with total mass estimates from 13CO to infer the shape of the density probability distribution function. This method is complementary to measurements of turbulence via the column density distribution and should be applicable to any molecular cloud with detected CO. We show that turbulence in this cloud is probably compressively driven, with a compressive-to-total Mach number ratio $b = \mathcal{M}_C/\mathcal{M}>0.4$. We measure the standard deviation of the density distribution, constraining it to the range $1.5 < σ_s < 1.9$ assuming that the density is lognormally distributed. This measurement represents an essential input into star formation laws. The method of averaging over different excitation conditions to produce a model of emission from a turbulent cloud is generally applicable to optically thin line observations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Connection between dense gas mass fraction, turbulence driving, and star formation efficiency of molecular clouds

We examine the physical parameters that affect the accumulation of gas in molecular clouds to high column densities where the formation of stars takes place. In particular, we analyze the dense gas mass fraction (DGMF) in a set of self-gravitating, isothermal, magnetohydrodynamic turbulence simulations including sink particles to model star formation. We find that the simulations predict close to exponential DGMFs over the column density range N(H2) = 3-25 x 10^{21} cm^{-2} that can be easily probed via, e.g., dust extinction measurements. The exponential slopes correlate with the type of turbulence driving and also with the star formation efficiency. They are almost uncorrelated with the sonic Mach number and magnetic-field strength. The slopes at early stages of cloud evolution are steeper than at the later stages. A comparison of these predictions with observations shows that only simulations with relatively non-compressive driving (b ~< 0.4) agree with the DGMFs of nearby molecular clouds. Massive infrared dark clouds can show DGMFs that are in agreement with more compressive driving. The DGMFs of molecular clouds can be significantly affected by how compressive the turbulence is on average. Variations in the level of compression can cause scatter to the DGMF slopes, and some variation is indeed necessary to explain the spread of the observed DGMF slopes. The observed DGMF slopes can also be affected by the clouds' star formation activities and statistical cloud-to-cloud variations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the universality of supersonic turbulence

Compressible turbulence shapes the structure of the interstellar medium of our Galaxy and likely plays an important role also during structure formation in the early Universe. The density PDF and the power spectrum of such compressible, supersonic turbulence are the key ingredients for theories of star formation. However, both the PDF and the spectrum are still a matter of debate, because theoretical predictions are limited and simulations of supersonic turbulence require enormous resolutions to capture the inertial-range scaling. To advance our limited knowledge of compressible turbulence, we here present and analyse the world's largest simulations of supersonic turbulence. We compare hydrodynamic models with numerical resolutions of 256^3-4096^3 mesh points and with two distinct driving mechanisms, solenoidal (divergence-free) driving and compressive (curl-free) driving. We find convergence of the density PDF, with compressive driving exhibiting a much wider and more intermittent density distribution than solenoidal driving. Analysing the power spectrum of the turbulence, we find a pure velocity scaling close to Burgers turbulence with P(v) k^(-2) for both driving modes in our hydrodynamical simulations with Mach = 17. The spectrum of the density-weighted velocity rho^(1/3)v, however, does not provide the previously suggested universal scaling for supersonic turbulence. We find that the power spectrum P(rho^(1/3)v) scales with wavenumber as k^(-1.74) for solenoidal driving, close to incompressible Kolmogorov turbulence, k^(-5/3), but is significantly steeper with k^(-2.10) for compressive driving. We show that this is consistent with a recent theoretical model for compressible turbulence that predicts P(rho^(1/3)v) k^(-19/9) in the presence of a strong div(v) component as is produced by compressive driving and remains remarkably constant throughout the supersonic turbulent cascade.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Smoothed Particle Magnetohydrodynamics Simulations of Protostellar Jets and Turbulent Dynamos

We presents results from Smoothed Particle Magnetohydrodynamics simulations of collapsing molecular cloud cores, and dynamo amplification of the magnetic field in the presence of Mach 10 magnetised turbulence. Our star formation simulations have produced, for the first time ever, highly collimated magnetised protostellar jets from the first hydrostatic core phase. Up to 40% of the initial core mass may be ejected through this outflow. The primary difficulty in performing these simulations is maintaining the divergence free constraint of the magnetic field, and to address this issue, we have developed a new divergence cleaning method which has allowed us to stably follow the evolution of these protostellar jets for long periods. The simulations performed of supersonic MHD turbulence are able to exponentially amplify magnetic energy by up to 10 orders of magnitude via turbulent dynamo. To reduce numerical dissipation, a new shock detection algorithm is utilised which is able to track magnetic shocks throughout a large range of magnetic field strengths.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The origin of physical variations in the star formation law

Observations of external galaxies and of local star-forming clouds in the Milky Way have suggested a variety of star formation laws, i.e., simple direct relations between the column density of star formation (Sigma_SFR: the amount of gas forming stars per unit area and time) and the column density of available gas (Sigma_gas). Extending previous studies, we show that these different, sometimes contradictory relations for Milky Way clouds, nearby galaxies, and high-redshift discs and starbursts can be combined in one universal star formation law in which Sigma_SFR is about 1% of the local gas collapse rate, Sigma_gas/t_ff, but a significant scatter remains in this relation. Using computer simulations and theoretical models, we find that the observed scatter may be primarily controlled by physical variations in the Mach number of the turbulence and by differences in the star formation efficiency. Secondary variations can be induced by changes in the virial parameter, turbulent driving and magnetic field. The predictions of our models are testable with observations that constrain both the Mach number and the star formation efficiency in Milky Way clouds, external disc and starburst galaxies at low and high redshift. We also find that reduced telescope resolution does not strongly affect such measurements when Sigma_SFR is plotted against Sigma_gas/t_ff.

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Star Formation Rate of Molecular Clouds

We review recent advances in the analytical and numerical modeling of the star formation rate in molecular clouds and discuss the available observational constraints. We focus on molecular clouds as the fundamental star formation sites, rather than on the larger-scale processes that form the clouds and set their properties. Molecular clouds are shaped into a complex filamentary structure by supersonic turbulence, with only a small fraction of the cloud mass channeled into collapsing protostars over a free-fall time of the system. In recent years, the physics of supersonic turbulence has been widely explored with computer simulations, leading to statistical models of this fragmentation process, and to the prediction of the star formation rate as a function of fundamental physical parameters of molecular clouds, such as the virial parameter, the rms Mach number, the compressive fraction of the turbulence driver, and the ratio of gas to magnetic pressure. Infrared space telescopes, as well as ground-based observatories have provided unprecedented probes of the filamentary structure of molecular clouds and the location of forming stars within them.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A new density variance - Mach number relation for subsonic and supersonic, isothermal turbulence

The probability density function (PDF) of the gas density in subsonic and supersonic, isothermal, driven turbulence is analyzed with a systematic set of hydrodynamical grid simulations with resolutions up to 1024^3 cells. We performed a series of numerical experiments with root mean square (r.m.s.) Mach number M ranging from the nearly incompressible, subsonic (M=0.1) to the highly compressible, supersonic (M=15) regime. We study the influence of two extreme cases for the driving mechanism by applying a purely solenoidal (divergence-free) and a purely compressive (curl-free) forcing field to drive the turbulence. We find that our measurements fit the linear relation between the r.m.s. Mach number and the standard deviation of the density distribution in a wide range of Mach numbers, where the proportionality constant depends on the type of the forcing. In addition, we propose a new linear relation between the standard deviation of the density distribution and the standard deviation of the velocity in compressible modes, i.e. the compressible component of the r.m.s. Mach number. In this relation the influence of the forcing is significantly reduced, suggesting a linear relation between the standard deviation of the density distribution and the standard deviation of the velocity in compressible modes, independent of the forcing, ranging from the subsonic to the supersonic regime.

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

On the Star Formation Efficiency of Turbulent Magnetized Clouds

We study the star formation efficiency (SFE) in simulations and observations of turbulent, magnetized, molecular clouds. We find that the probability density functions (PDFs) of the density and the column density in our simulations with solenoidal, mixed, and compressive forcing of the turbulence, sonic Mach numbers of 3-50, and magnetic fields in the super- to the trans-Alfvenic regime, all develop power-law tails of flattening slope with increasing SFE. The high-density tails of the PDFs are consistent with equivalent radial density profiles, rho ~ r^(-kappa) with kappa ~ 1.5-2.5, in agreement with observations. Studying velocity-size scalings, we find that all the simulations are consistent with the observed v ~ l^(1/2) scaling of supersonic turbulence, and seem to approach Kolmogorov turbulence with v ~ l^(1/3) below the sonic scale. The velocity-size scaling is, however, largely independent of the SFE. In contrast, the density-size and column density-size scalings are highly sensitive to star formation. We find that the power-law slope alpha of the density power spectrum, P(rho,k) ~ k^alpha, or equivalently the Delta-variance spectrum of column density, DV(Sigma,l) ~ l^(-alpha), switches sign from alpha < 0 for SFE ~ 0 to alpha > 0 when star formation proceeds (SFE > 0). We provide a relation to compute the SFE from a measurement of alpha. Studying the literature, we find values ranging from alpha = -1.6 to +1.6 in observations covering scales from the large-scale atomic medium, over cold molecular clouds, down to dense star-forming cores. From those alpha values, we infer SFEs and find good agreement with independent measurements based on young stellar object (YSO) counts, where available. Our SFE-alpha relation provides an independent estimate of the SFE based on the column density map of a cloud alone, without requiring a priori knowledge of star-formation activity or YSO counts.

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

Star Formation in the First Galaxies I: Collapse Delayed by Lyman-Werner Radiation

We investigate the process of metal-free star formation in the first galaxies with a high-resolution cosmological simulation. We consider the cosmologically motivated scenario in which a strong molecule-destroying Lyman-Werner (LW) background inhibits effective cooling in low-mass haloes, delaying star formation until the collapse or more massive haloes. Only when molecular hydrogen (H2) can self-shield from LW radiation, which requires a halo capable of cooling by atomic line emission, will star formation be possible. To follow the formation of multiple gravitationally bound objects, at high gas densities we introduce sink particles which accrete gas directly from the computational grid. We find that in a 1 Mpc^3 (comoving) box, runaway collapse first occurs in a 3x10^7 M_sun dark matter halo at z~12 assuming a background intensity of J21=100. Due to a runaway increase in the H2 abundance and cooling rate, a self-shielding, supersonically turbulent core develops abruptly with ~10^4 M_sun in cold gas available for star formation. We analyze the formation of this self-shielding core, the character of turbulence, and the prospects for star formation. Due to a lack of fragmentation on scales we resolve, we argue that LW-delayed metal-free star formation in atomic cooling haloes is very similar to star formation in primordial minihaloes, although in making this conclusion we ignore internal stellar feedback. Finally, we briefly discuss the detectability of metal-free stellar clusters with the James Webb Space Telescope.

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 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
preprint2012arXiv

The Star Formation Rate of Turbulent Magnetized Clouds: Comparing Theory, Simulations, and Observations

We derive and compare six theoretical models for the star formation rate (SFR) - the Krumholz & McKee (KM), Padoan & Nordlund (PN), and Hennebelle & Chabrier (HC) models, and three multi-freefall versions of these, suggested by HC - all based on integrals over the log-normal distribution of turbulent gas. We extend all theories to include magnetic fields, and show that the SFR depends on four basic parameters: (1) virial parameter alpha_vir; (2) sonic Mach number M; (3) turbulent forcing parameter b, which is a measure for the fraction of energy driven in compressive modes; and (4) plasma beta=2(M_A/M)^2 with the Alfven Mach number M_A. We compare all six theories with MHD simulations, covering cloud masses of 300 to 4x10^6 solar masses and Mach numbers M = 3 to 50 and M_A = 1 to infinity, with solenoidal (b=1/3), mixed (b=0.4) and compressive turbulent (b=1) forcings. We find that the SFR increases by a factor of four between M=5 and 50 for compressive forcing and alpha_vir~1. Comparing forcing parameters, we see that the SFR is more than 10x higher with compressive than solenoidal forcing for Mach 10 simulations. The SFR and fragmentation are both reduced by a factor of two in strongly magnetized, trans-Alfvenic turbulence compared to hydrodynamic turbulence. All simulations are fit simultaneously by the multi-freefall KM and multi-freefall PN theories within a factor of two over two orders of magnitude in SFR. The simulated SFRs cover the range and correlation of SFR column density with gas column density observed in Galactic clouds, and agree well for star formation efficiencies SFE = 1% to 10% and local efficiencies epsilon = 0.3 to 0.7 due to feedback. We conclude that the SFR is primarily controlled by interstellar turbulence, with a secondary effect coming from magnetic fields.

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

A robust numerical scheme for highly compressible magnetohydrodynamics: Nonlinear stability, implementation and tests

The ideal MHD equations are a central model in astrophysics, and their solution relies upon stable numerical schemes. We present an implementation of a new method, which possesses excellent stability properties. Numerical tests demonstrate that the theoretical stability properties are valid in practice with negligible compromises to accuracy. The result is a highly robust scheme with state-of-the-art efficiency. The scheme's robustness is due to entropy stability, positivity and properly discretised Powell terms. The implementation takes the form of a modification of the MHD module in the FLASH code, an adaptive mesh refinement code. We compare the new scheme with the standard FLASH implementation for MHD. Results show comparable accuracy to standard FLASH with the Roe solver, but highly improved efficiency and stability, particularly for high Mach number flows and low plasma beta. The tests include 1D shock tubes, 2D instabilities and highly supersonic, 3D turbulence. We consider turbulent flows with RMS sonic Mach numbers up to 10, typical of gas flows in the interstellar medium. We investigate both strong initial magnetic fields and magnetic field amplification by the turbulent dynamo from extremely high plasma beta. The energy spectra show a reasonable decrease in dissipation with grid refinement, and at a resolution of 512^3 grid cells we identify a narrow inertial range with the expected power-law scaling. The turbulent dynamo exhibits exponential growth of magnetic pressure, with the growth rate twice as high from solenoidal forcing than from compressive forcing. Two versions of the new scheme are presented, using relaxation-based 3-wave and 5-wave approximate Riemann solvers, respectively. The 5-wave solver is more accurate in some cases, and its computational cost is close to the 3-wave solver.

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

MHD simulations of a supernova-driven ISM and the warm ionized medium using a positivity preserving ideal MHD scheme

We present new 3D magnetohydrodynamic (MHD) simulations of a supernova-driven, stratified interstellar medium. These simulations were run using the Waagan (2009) positivity preserving scheme for ideal MHD implemented in the Flash code. The scheme is stable even for the Mach numbers approaching 100 found in this problem. We have previously shown that the density distribution arising from hydrodynamical versions of these simulations creates low-density pathways through which Lyman continuum photons can travel to heights |z| > 1 kpc. This naturally produces the warm ionized medium through photoionization due primarily to O stars near the plane. However, our earlier models reproduce the peak but not the width of the observed emission measure distribution. Here, we examine whether inclusion of magnetic fields and a greater vertical extent to the simulation domain produce a gas distribution that better matches the observations. We further study the change of magnetic energy over time in our models, showing that it appears to reach a steady state after a few hundred megayears, presumably supported by a turbulent dynamo driven by the supernova explosions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Modeling H2 formation in the turbulent ISM: Solenoidal versus compressive turbulent forcing

We present results from high-resolution three-dimensional simulations of the turbulent interstellar medium that study the influence of the nature of the turbulence on the formation of molecular hydrogen. We have examined both solenoidal (divergence-free) and compressive (curl-free) turbulent driving, and show that compressive driving leads to faster H2 formation, owing to the higher peak densities produced in the gas. The difference in the H2 formation rate can be as much as an order of magnitude at early times, but declines at later times as the highest density regions become fully molecular and stop contributing to the total H2 formation rate. We have also used our results to test a simple prescription suggested by Gnedin et al. (2009) for modeling the influence of unresolved density fluctuations on the H2 formation rate in large-scale simulations of the ISM. We find that this approach works well when the H2 fraction is small, but breaks down once the highest density gas becomes fully molecular.

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

Statistical properties of supersonic turbulence in the Lagrangian and Eulerian frameworks

We present a systematic study of the influence of different forcing types on the statistical properties of supersonic, isothermal turbulence in both the Lagrangian and Eulerian frameworks. We analyse a series of high-resolution, hydrodynamical grid simulations with Lagrangian tracer particles and examine the effects of solenoidal (divergence-free) and compressive (curl-free) forcing on structure functions, their scaling exponents, and the probability density functions of the gas density and velocity increments. Compressively driven simulations show a significantly larger density contrast, a more intermittent behaviour, and larger fractal dimension of the most dissipative structures at the same root mean square Mach number. We show that the absolute values of Lagrangian and Eulerian structure functions of all orders in the integral range are only a function of the root mean square Mach number, but independent of the forcing. With the assumption of a Gaussian distribution for the probability density function of the velocity increments on large scales, we derive a model that describes this behaviour.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Turbulence Spectrum of Molecular Clouds in the Galactic Ring Survey: A Density-Dependent PCA Calibration

Turbulence plays a major role in the formation and evolution of molecular clouds. The problem is that turbulent velocities are convolved with the density of an observed region. To correct for this convolution, we investigate the relation between the turbulence spectrum of model clouds, and the statistics of their synthetic observations obtained from Principal Component Analysis (PCA). We apply PCA to spectral maps generated from simulated density and velocity fields, obtained from hydrodynamic simulations of supersonic turbulence, and from fractional Brownian motion fields with varying velocity, density spectra, and density dispersion. We examine the dependence of the slope of the PCA structure function, alpha_PCA, on intermittency, on the turbulence velocity (beta_v) and density (beta_n) spectral indexes, and on density dispersion. We find that PCA is insensitive to beta_n and to the log-density dispersion sigma_s, provided sigma_s < 2. For sigma_s>2, alpha_PCA increases with sigma_s due to the intermittent sampling of the velocity field by the density field. The PCA calibration also depends on intermittency. We derive a PCA calibration based on fBms with sigma_s<2 and apply it to 367 CO spectral maps of molecular clouds in the Galactic Ring Survey. The average slope of the PCA structure function, <alpha_PCA>=0.62\pm0.2, is consistent with the hydrodynamic simulations and leads to a turbulence velocity exponent <beta_v>=2.06\pm0.6 for a non-intermittent, low density dispersion flow. Accounting for intermittency and density dispersion, the coincidence between the PCA slope of the GRS clouds and the hydrodynamic simulations suggests beta_v~1.9, consistent with both Burgers and compressible intermittent turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A comparison between grid and particle methods on the statistics of driven, supersonic, isothermal turbulence

We compare the statistics of driven, supersonic turbulence at high Mach number using FLASH a widely used Eulerian grid-based code and PHANTOM, a Lagrangian SPH code at resolutions of up to 512^3 in both grid cells and SPH particles. We find excellent agreement between codes on the basic statistical properties: a slope of k^-1.95 in the velocity power spectrum for hydrodynamic, Mach 10 turbulence, evidence in both codes for a Kolmogorov-like slope of k^-5/3 in the variable rho^1/3 v as suggested by Kritsuk et al. and a log-normal PDF with a width that scales with Mach number and proportionality constant b=0.33-0.5 in the density variance-Mach number relation. The measured structure function slopes are not converged in either code at 512^3 elements. We find that, for measuring volumetric statistics such as the power spectrum slope and structure function scaling, SPH and grid codes give roughly comparable results when the number of SPH particles is approximately equal to the number of grid cells. In particular, to accurately measure the power spectrum slope in the inertial range, in the absence of subgrid models, requires at least 512^3 computational elements in either code. On the other hand the SPH code was found to be better at resolving dense structures, giving max. densities at a resolution of 128^3 particles that were similar those resolved in the grid code at 512^3 cells, reflected also in the high density tail of the PDF. We find SPH to be more dissipative at comparable numbers of computational elements in statistics of the velocity field, but correspondingly less dissipative than the grid code in the statistics of density weighted quantities such as rho^1/3 v. For SPH simulations of high Mach number turbulence we find it important to use sufficient non-linear beta-viscosity to prevent particle interpenetration in shocks (we require beta=4 instead of the default beta=2).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A method for reconstructing the PDF of a 3D turbulent density field from 2D observations

We introduce a method for calculating the probability density function (PDF) of a turbulent density field in three dimensions using only information contained in the projected two-dimensional column density field. We test the method by applying it to numerical simulations of hydrodynamic and magnetohydrodynamic turbulence in molecular clouds. To a good approximation, the PDF of log(normalised column density) is a compressed, shifted version of the PDF of log(normalised density). The degree of compression can be determined observationally from the column density power spectrum, under the assumption of statistical isotropy of the turbulence.

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

Numerical and semi-analytic core mass distributions in supersonic isothermal turbulence

We investigate the influence of the turbulence forcing on the mass distributions of gravitationally unstable cores by postprocessing data from simulations of non-selfgravitating isothermal supersonic turbulence with varying resolution. In one set of simulations solenoidal forcing is applied, while the second set uses purely compressive forcing to excite turbulent motions. From the resulting density field, we compute the mass distribution of gravitationally unstable cores by means of a clump-finding algorithm. Using the time-averaged probability density functions of the mass density, semi-analytic mass distributions are calculated from analytical theories. We apply stability criteria that are based on the Bonnor-Ebert mass resulting from the thermal pressure and from the sum of thermal and turbulent pressure. Although there are uncertainties in the application of the clump-finding algorithm, we find systematic differences in the mass distributions obtained from solenoidal and compressive forcing. Compressive forcing produces a shallower slope in the high-mass power-law regime compared to solenoidal forcing. The mass distributions also depend on the Jeans length resulting from the choice of the mass in the computational box, which is freely scalable for non-selfgravitating isothermal turbulence. Provided that all cores are numerically resolved and most cores are small compared to the length scale of the forcing, the normalised core mass distributions are found to be close to the semi-analytic models. Especially for the high-mass tails, the Hennebelle-Chabrier theory implies that the additional support due to turbulent pressure is important.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Smoothed Particle Hydrodynamics: Turbulence and MHD

In this paper we discuss recent applications of the Smoothed Particle Hydrodynamics (SPH) method to the simulation of supersonic turbulence in the interstellar medium, as well as giving an update on recent algorithmic developments in solving the equations of magnetohydrodynamics (MHD) in SPH. Using high resolution calculations (up to 134 million particles), we find excellent agreement with grid-based results on a range of measures including the power spectrum slope in both the velocity field and the density-weighted velocity rho^(1/3) v, the latter showing a Kolmogorov-like k^-5/3 scaling as proposed by Kritsuk et al. (2007). We also find good agreement on the statistics of the Probability Distribution Function (PDF) and structure functions, independently confirming the scaling found by Schmidt, Federrath & Klessen (2008). On Smoothed Particle Magnetohydrodynamics (SPMHD) we have recently wasted a great deal of time and effort investigating the vector potential as an alternative to the Euler potentials formulation, in the end concluding that using the vector potential has even more severe problems than the standard (B-field based) SPMHD approach.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The density variance -- Mach number relation in supersonic, isothermal turbulence

We examine the relation between the density variance and the mean-square Mach number in supersonic, isothermal turbulence, assumed in several recent analytic models of the star formation process. From a series of calculations of supersonic, hydrodynamic turbulence driven using purely solenoidal Fourier modes, we find that the `standard' relationship between the variance in the log of density and the Mach number squared, i.e., sigma^2_(ln rho/rhobar)=ln (1+b^2 M^2), with b = 1/3 is a good fit to the numerical results in the supersonic regime up to at least Mach 20, similar to previous determinations at lower Mach numbers. While direct measurements of the variance in linear density are found to be severely underestimated by finite resolution effects, it is possible to infer the linear density variance via the assumption of log-normality in the Probability Distribution Function. The inferred relationship with Mach number, consistent with sigma_(rho/rhobar) ~ b M with b=1/3, is, however, significantly shallower than observational determinations of the relationship in the Taurus Molecular Cloud and IC5146 (both consistent with b~ 0.5), implying that additional physics such as gravity is important in these clouds and/or that turbulent driving in the ISM contains a significant compressive component. Magnetic fields are not found to change this picture significantly, in general reducing the measured variances and thus worsening the discrepancy with observations.

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