BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Christoph Federrath

Christoph Federrath contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.GAastro-ph.SRphysics.flu-dynastro-ph.COastro-ph.IMphysics.plasm-phastro-ph.HEphysics.comp-phDistributed, Parallel, and Cluster Computinghep-phPerformancephysics.ao-ph

Trust snapshot

Quick read

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

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

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

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

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

Building this graph slice

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

Published work

35 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&#39; 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 &#34;macroscopic&#34; 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&#39;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 &#39;small-scale turbulent dynamo&#39; and the &#39;large-scale mean-field dynamo&#39;. 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&#39;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&#39;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