-
Magnetized compressible turbulence with a fluctuation dynamo and Reynolds numbers over a million
Authors:
James R. Beattie,
Christoph Federrath,
Ralf S. Klessen,
Salvatore Cielo,
Amitava Bhattacharjee
Abstract:
Supersonic magnetohydrodynamic (MHD) turbulence is a ubiquitous state for many astrophysical plasmas. However, even the basic statistics for this type of turbulence remains uncertain. We present results from supersonic MHD turbulence simulations at unparalleled resolutions, with plasma Reynolds numbers of over a million. In the kinetic energy spectrum we find a break between the scales that are do…
▽ More
Supersonic magnetohydrodynamic (MHD) turbulence is a ubiquitous state for many astrophysical plasmas. However, even the basic statistics for this type of turbulence remains uncertain. We present results from supersonic MHD turbulence simulations at unparalleled resolutions, with plasma Reynolds numbers of over a million. In the kinetic energy spectrum we find a break between the scales that are dominated by kinetic energy, with spectral index $-2$, and those that become strongly magnetized, with spectral index $-3/2$. By analyzing the Helmholtz decomposed kinetic energy spectrum, we find that the compressible modes are not passively mixed through the cascade of the incompressible modes. At high magnetic Reynolds number, above $10^5$, we find a power law in the magnetic energy spectrum with spectral index $-9/5$. On the strongly magnetized, subsonic scales the plasma tends to self-organize into locally relaxed regions, where there is strong alignment between the current density, magnetic field, velocity field and vorticity field, depleting both the nonlinearities and magnetic terms in the MHD equations, which we attribute to plasma relaxation on scales where the magnetic fluctuations evolve on shorter timescales than the velocity fluctuations. This process constrains the cascade to inhomogenous, volume-poor, fractal surfaces between relaxed regions, which has significant repercussions for understanding the nature of magnetized turbulence in astrophysical plasmas and the saturation of the fluctuation dynamo.
△ Less
Submitted 29 May, 2024; v1 submitted 26 May, 2024;
originally announced May 2024.
-
Taking control of compressible modes: bulk viscosity and the turbulent dynamo
Authors:
James R. Beattie,
Christoph Federrath,
Neco Kriel,
Justin Kin Jun Hew,
Amitava Bhattacharjee
Abstract:
Many polyatomic astrophysical plasmas are compressible and out of chemical and thermal equilibrium, and yet, due to Stokes' hypothesis, a means to carefully control the decay of compressible modes in these systems has largely been neglected. This is especially important for small-scale, turbulent dynamo processes, which are known to be sensitive to the effects of compression. To control the viscou…
▽ More
Many polyatomic astrophysical plasmas are compressible and out of chemical and thermal equilibrium, and yet, due to Stokes' hypothesis, a means to carefully control the decay of compressible modes in these systems has largely been neglected. This is especially important for small-scale, turbulent dynamo processes, which are known to be sensitive to the effects of compression. To control the viscous properties of the compressible modes, we perform supersonic, visco-resistive dynamo simulations with additional bulk viscosity $ν_{\rm bulk}$, deriving a new $ν_{\rm bulk}$ Reynolds number $\rm{Re}_{\rm bulk}$, and viscous Prandtl number $\rm{P}ν\equiv \rm{Re}_{\rm bulk} / \rm{Re}_{\rm shear}$, where $\rm{Re}_{\rm shear}$ is the shear viscosity Reynolds number. For $10^{-3} \leq \rm{P}ν\leq \infty$, we explore a broad range of statistics critical to the dynamo problem, including the integral and spectral energy ratios, growth rates, and the magnetic $E_{\rm mag}(k)$ and kinetic $E_{\rm kin}(k)$ energy spectrum. We derive a general framework for decomposing $E_{\rm mag}$ growth rates into incompressible and compressible terms via orthogonal tensor decompositions of $\nabla\otimes\mathbf{v}$, where $\mathbf{v}$ is the fluid velocity. We find that compressible modes play a dual role, growing and decaying $E_{\rm mag}$, and that field-line stretching is the main driver of growth, even in supersonic dynamos. In the absence of $ν_{\rm bulk}$, compressible modes pile up on small-scales, creating an apparent spectral bottleneck, which disappears for $\rm{P}ν\approx 1$. As $\rm{P}ν$ decreases, compressible modes are dissipated at increasingly larger scales, in turn suppressing incompressible modes through a coupling between viscosity operators. We emphasise the importance of further understanding the role of $ν_{\rm bulk}$ in compressible astrophysical plasmas.
△ Less
Submitted 11 December, 2023; v1 submitted 6 December, 2023;
originally announced December 2023.
-
Fundamental MHD scales -- II: the kinematic phase of the supersonic small-scale dynamo
Authors:
Neco Kriel,
James R. Beattie,
Christoph Federrath,
Mark R. Krumholz,
Justin Kin Jun Hew
Abstract:
The small-scale dynamo (SSD) amplifies weak magnetic fields exponentially fast via kinetic motions. While there exist well-established theories for SSDs in incompressible flows, many astrophysical SSDs operate in supersonic turbulence. To understand the impact of compressibility on amplified magnetic fields, we perform an extensive set of visco-resistive SSD simulations, covering a wide range of s…
▽ More
The small-scale dynamo (SSD) amplifies weak magnetic fields exponentially fast via kinetic motions. While there exist well-established theories for SSDs in incompressible flows, many astrophysical SSDs operate in supersonic turbulence. To understand the impact of compressibility on amplified magnetic fields, we perform an extensive set of visco-resistive SSD simulations, covering a wide range of sonic Mach number $\mathcal{M}$, hydrodynamic Reynolds number Re, and magnetic Prandtl number Pm. We develop robust methods for measuring kinetic and magnetic energy dissipation scales $\ell_ν$ and $\ell_η$, as well as the scale at which magnetic fields are strongest $\ell_p$ during the kinematic phase of these simulations. We show that $\ell_ν/\ell_η\sim$ Pm$^{1/2}$ is a universal feature in the kinematic phase of Pm $\geq 1$ SSDs, regardless of $\mathcal{M}$ or Re, and we confirm earlier predictions that SSDs operating in incompressible plasmas (either $\mathcal{M} \leq 1$ or Re $<$ Re$_{\rm crit} \approx 100$) concentrate magnetic energy at the smallest scales allowed by magnetic dissipation, $\ell_p \sim \ell_η$, and produce fields organised with field strength and field-line curvature inversely correlated. However, we show that these predictions fail for compressible SSDs ($\mathcal{M} > 1$ and Re $>$ Re$_{\rm crit}$), where shocks concentrate magnetic energy in large-scale, over-dense, coherent structures, with size $\ell_p \sim (\ell_{\rm turb} / \ell_{\rm shock})^{1/3} \ell_η\gg \ell_η$, where $\ell_{\rm shock} \sim \mathcal{M}^2 / [$Re $ (\mathcal{M} - 1)^2]$ is shock width, and $\ell_{\rm turb}$ is the turbulent outer scale; magnetic field-line curvature becomes almost independent of the field strength. We discuss the implications for galaxy mergers and for cosmic-ray transport models in the interstellar medium that are sensitive to field-line curvature statistics.
△ Less
Submitted 25 October, 2023;
originally announced October 2023.
-
Growth or Decay -- I: universality of the turbulent dynamo saturation
Authors:
James R. Beattie,
Christoph Federrath,
Neco Kriel,
Philip Mocz,
Amit Seta
Abstract:
The turbulent small-scale dynamo (SSD) is likely to be responsible for the magnetisation of the interstellar medium (ISM) that we observe in the Universe today. The SSD efficiently converts kinetic energy $E_{\rm kin}$ into magnetic energy $E_{\rm mag} $, and is often used to explain how an initially weak magnetic field with $E_{\rm mag} \ll E_{\rm kin}$ is amplified, and then maintained at a leve…
▽ More
The turbulent small-scale dynamo (SSD) is likely to be responsible for the magnetisation of the interstellar medium (ISM) that we observe in the Universe today. The SSD efficiently converts kinetic energy $E_{\rm kin}$ into magnetic energy $E_{\rm mag} $, and is often used to explain how an initially weak magnetic field with $E_{\rm mag} \ll E_{\rm kin}$ is amplified, and then maintained at a level $E_{\rm mag} \lesssim E_{\rm kin}$. Usually, this process is studied by initialising a weak seed magnetic field and letting the turbulence grow it to saturation. However, in this Part I of the Growth or Decay series, using three-dimensional, visco-resistive magnetohydrodynamical turbulence simulations up to magnetic Reynolds numbers of 2000, we show that the same final state in the integral quantities, energy spectra, and characteristic scales of the magnetic field can also be achieved if initially $E_{\rm mag} \sim E_{\rm kin}$ or even if initially $E_{\rm mag} \gg E_{\rm kin}$. This suggests that the final saturated state of the turbulent dynamo is set by the turbulence and the material properties of the plasma, independent of the initial structure or amplitude of the magnetic field. We discuss the implications this has for the maintenance of magnetic fields in turbulent plasmas and future studies exploring the dynamo saturation.
△ Less
Submitted 15 June, 2023; v1 submitted 21 September, 2022;
originally announced September 2022.
-
Understanding star formation in molecular clouds IV. Column density PDFs from quiescent to massive molecular clouds
Authors:
N. Schneider,
V. Ossenkopf-Okada,
S. Clarke,
R. S. Klessen,
S. Kabanovic,
T. Veltchev,
S. Bontemps,
S. Dib,
T. Csengeri,
C. Federrath,
J. Di Francesco,
F. Motte,
Ph. Andre,
D. Arzoumanian,
J. R. Beattie,
L. Bonne,
P. Didelon,
D. Elia,
V. Koenyves,
A. Kritsuk,
B. Ladjelate,
Ph. Myers,
S. Pezzuto,
J. F. Robitaille,
A. Roy
, et al. (4 additional authors not shown)
Abstract:
We present N-PDFs of 29 Galactic regions obtained from Herschel imaging at high angular resolution, covering diffuse and quiescent clouds, and those showing low-, intermediate-, and high-mass star formation (SF), and characterize the cloud structure using the Delta-variance tool. The N-PDFs are double-log-normal at low column densities, and display one or two power law tails (PLTs) at higher colum…
▽ More
We present N-PDFs of 29 Galactic regions obtained from Herschel imaging at high angular resolution, covering diffuse and quiescent clouds, and those showing low-, intermediate-, and high-mass star formation (SF), and characterize the cloud structure using the Delta-variance tool. The N-PDFs are double-log-normal at low column densities, and display one or two power law tails (PLTs) at higher column densities. For diffuse, quiescent, and low-mass SF clouds, we propose that the two log-normals arise from the atomic and molecular phase, respectively. For massive clouds, we suggest that the first log-normal is built up by turbulently mixed H2 and the second one by compressed (via stellar feedback) molecular gas. Nearly all clouds have two PLTs with slopes consistent with self-gravity, where the second one can be flatter or steeper than the first one. A flatter PLT could be caused by stellar feedback or other physical processes that slow down collapse and reduce the flow of mass toward higher densities. The steeper slope could arise if the magnetic field is oriented perpendicular to the LOS column density distribution. The first deviation point (DP), where the N-PDF turns from log-normal into a PLT, shows a clustering around values of a visual extinction of AV (DP1) around 2-5. The second DP, which defines the break between the two PLTs, varies strongly. Using the Delta-variance, we observe that the AV value, where the slope changes between the first and second PLT, increases with the characteristic size scale in the variance spectrum. We conclude that at low column densities, atomic and molecular gas is turbulently mixed, while at high column densities, the gas is fully molecular and dominated by self-gravity. The best fitting model N-PDFs of molecular clouds is thus one with log-normal low column density distributions, followed by one or two PLTs.
△ Less
Submitted 29 July, 2022;
originally announced July 2022.
-
Turbulent diffusion of streaming cosmic rays in compressible, partially ionised plasma
Authors:
Matt L. Sampson,
James R. Beattie,
Mark R. Krumholz,
Roland M. Crocker,
Christoph Federrath,
Amit Seta
Abstract:
Cosmic rays (CRs) are a dynamically important component of the interstellar medium (ISM) of galaxies. The $\sim$GeV CRs that carry most CR energy and pressure are likely confined by self-generated turbulence, leading them to stream along magnetic field lines at the ion Alfvén speed. However, the consequences of self-confinement for CR propagation on galaxy scales remain highly uncertain. In this p…
▽ More
Cosmic rays (CRs) are a dynamically important component of the interstellar medium (ISM) of galaxies. The $\sim$GeV CRs that carry most CR energy and pressure are likely confined by self-generated turbulence, leading them to stream along magnetic field lines at the ion Alfvén speed. However, the consequences of self-confinement for CR propagation on galaxy scales remain highly uncertain. In this paper, we use a large ensemble of magnetohydrodynamical turbulence simulations to quantify how the basic parameters describing ISM turbulence -- the sonic Mach number, $\mathcal{M}$ (plasma compressibility), Alfvén Mach number, $\mathcal{M}_{A0}$ (strength of the large-scale field with respect to the turbulence), and ionisation fraction by mass, $χ$ -- affect the transport of streaming CRs. We show that the large-scale transport of CRs whose small-scale motion consists of streaming along field lines is well described as a combination of streaming along the mean field and superdiffusion both along (parallel to) and across (perpendicular to) it; $\mathcal{M}_{A0}$ drives the level of anisotropy between parallel and perpendicular diffusion and $χ$ modulates the magnitude of the diffusion coefficients, while in our choice of units, $\mathcal{M}$ is unimportant except in the sub-Alfvénic ($\mathcal{M}_{A0} \lesssim 0.5$) regime. Our finding that superdiffusion is ubiquitous potentially explains the apparent discrepancy between CR diffusion coefficients inferred from measurements close to individual sources compared to those measured on larger, Galactic scales. Finally, we present empirical fits for the diffusion coefficients as a function of plasma parameters that may be used as sub-grid recipes for global interstellar medium, galaxy or cosmological simulations.
△ Less
Submitted 3 November, 2022; v1 submitted 17 May, 2022;
originally announced May 2022.
-
Fundamental scales in the kinematic phase of the turbulent dynamo
Authors:
Neco Kriel,
James R. Beattie,
Amit Seta,
Christoph Federrath
Abstract:
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),…
▽ More
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.
△ Less
Submitted 6 April, 2022; v1 submitted 2 April, 2022;
originally announced April 2022.
-
Ion Alfvén velocity fluctuations and implications for the diffusion of streaming cosmic rays
Authors:
James R. Beattie,
Mark R. Krumholz,
Christoph Federrath,
Matt Sampson,
Roland M. Crocker
Abstract:
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 ver…
▽ More
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.
△ Less
Submitted 7 August, 2022; v1 submitted 25 March, 2022;
originally announced March 2022.
-
Energy balance and Alfvén Mach numbers in compressible magnetohydrodynamic turbulence with a large-scale magnetic field
Authors:
James R. Beattie,
Mark R. Krumholz,
Raphael Skalidis,
Christoph Federrath,
Amit Seta,
Roland M. Crocker,
Philip Mocz,
Neco Kriel
Abstract:
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-squa…
▽ More
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.
△ Less
Submitted 24 July, 2022; v1 submitted 25 February, 2022;
originally announced February 2022.
-
Theoretical and Observational Evidence for Coriolis Effects in Coronal Magnetic Fields Via Direct Current Driven Flaring Events
Authors:
Darryl Z. Seligman,
Leslie A. Rogers,
Adina D. Feinstein,
Mark R. Krumholz,
James R. Beattie,
Christoph Federrath,
Fred C. Adams,
Marco Fatuzzo,
Maximilian N. Günther
Abstract:
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 dis…
▽ More
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.
△ Less
Submitted 5 March, 2022; v1 submitted 10 January, 2022;
originally announced January 2022.
-
Why take the square root? An assessment of interstellar magnetic field strength estimation methods
Authors:
R. Skalidis,
J. Sternberg,
J. R. Beattie,
V. Pavlidou,
K. Tassis
Abstract:
The magnetic field strength in interstellar clouds can be estimated indirectly by using the spread of dust polarization angles ($δθ$). The method developed by Davis 1951 and by Chandrasekhar and Fermi 1953 (DCF) assumes that incompressible magnetohydrodynamic (MHD) fluctuations induce the observed dispersion of polarization angles, deriving $B\propto 1/δθ$ (or, $δθ\propto M_{A}$, in terms of the A…
▽ More
The magnetic field strength in interstellar clouds can be estimated indirectly by using the spread of dust polarization angles ($δθ$). The method developed by Davis 1951 and by Chandrasekhar and Fermi 1953 (DCF) assumes that incompressible magnetohydrodynamic (MHD) fluctuations induce the observed dispersion of polarization angles, deriving $B\propto 1/δθ$ (or, $δθ\propto M_{A}$, in terms of the Alfvénic Mach number). However, observations show that the interstellar medium (ISM) is highly compressible. Recently, Skalidis & Tassis 2021 (ST) relaxed the incompressibility assumption and derived instead $B\propto 1/\sqrt{δθ}$ ($δθ\propto M_{A}^2$). We explored what the correct scaling is in compressible and magnetized turbulence with numerical simulations. We used 26 magnetized, ideal-MHD numerical simulations with different types of forcing. The range of $M_{A}$ and sonic Mach numbers $M_{s}$ explored are $0.1 \leq M_{A} \leq 2.0$ and $0.5 \leq M_{s} \leq 20$. We created synthetic polarization maps and tested the assumptions and accuracy of the two methods. The synthetic data have a remarkable consistency with the $δθ\propto M_{A}^{2}$ scaling, which is inferred by ST, while the DCF scaling fails to follow the data. The ST method shows an accuracy better than $50\%$ over the entire range of $M_{A}$ explored; DCF performs adequately only in the range of $M_{A}$ for which it has been optimized through the use of a "fudge factor". For low $M_{A}$, DCF is inaccurate by factors of tens. The assumptions of the ST method reflect better the physical reality in clouds with compressible and magnetized turbulence, and for this reason the method provides a much better estimate of the magnetic field strength over the DCF method.
△ Less
Submitted 22 September, 2021;
originally announced September 2021.
-
The density distribution and physical origins of intermittency in supersonic, highly magnetised turbulence with diverse modes of driving
Authors:
James R. Beattie,
Philip Mocz,
Christoph Federrath,
Ralf S. Klessen
Abstract:
The probability density function (PDF) of the logarithmic density contrast, $s=\ln (ρ/ρ_0)$, with gas density $ρ$ and mean density $ρ_0$, for hydrodynamical supersonic turbulence is well-known to have significant non-Gaussian (intermittent) features that monotonically increase with the turbulent Mach number, $\mathcal{M}$. By studying the mass- and volume-weighted $s$-PDF for an ensemble of 36 sub…
▽ More
The probability density function (PDF) of the logarithmic density contrast, $s=\ln (ρ/ρ_0)$, with gas density $ρ$ and mean density $ρ_0$, for hydrodynamical supersonic turbulence is well-known to have significant non-Gaussian (intermittent) features that monotonically increase with the turbulent Mach number, $\mathcal{M}$. By studying the mass- and volume-weighted $s$-PDF for an ensemble of 36 sub-to-trans-Alfvenic mean-field, supersonic, isothermal turbulence simulations with different modes of driving, relevant to molecular gas in the cool interstellar medium, we show that a more intricate picture emerges for the non-Gaussian nature of $s$. Using four independent measures of the non-Gaussian components, we find hydrodynamical-like structure in the highly magnetised plasma for $\mathcal{M} \lesssim 4$. However, for $\mathcal{M} \gtrsim 4$, the non-Gaussian signatures disappear, leaving approximately Gaussian $s$-statistics -- exactly the opposite of hydrodynamical turbulence in the high-$\mathcal{M}$ limit. We also find that the non-Gaussian components of the PDF increase monotonically with more compressive driving modes. To understand the $\mathcal{M} \lesssim 4$ non-Gaussian features we use one-dimensional (1D) pencil beams to explore the dynamics along and across the large-scale magnetic field, $\mathbf{B}_0$. We discuss kinetic, density and magnetic field fluctuations from the pencil beams, and identify physical sources of non-Gaussian components to the PDF as single, strong shocks coupled to fast magnetosonic compressions that form along $\mathbf{B}_0$. We discuss the Gaussianisation of the $\mathcal{M} \gtrsim 4$ $s$-fields through the lens of two phenomenologies: the self-similarity of the $s$-field and homogenisation of the dynamical timescales between the over- and under-dense regions in the compressible gas.
△ Less
Submitted 15 October, 2022; v1 submitted 21 September, 2021;
originally announced September 2021.
-
First extragalactic measurement of the turbulence driving parameter: ALMA observations of the star-forming region N159E in the Large Magellanic Cloud
Authors:
Piyush Sharda,
Shyam H. Menon,
Christoph Federrath,
Mark R. Krumholz,
James R. Beattie,
Katherine E. Jameson,
Kazuki Tokuda,
Blakesley Burkhart,
Roland M. Crocker,
Charles J. Law,
Amit Seta,
Terrance J. Gaetz,
Nickolas M. Pingel,
Ivo R. Seitenzahl,
Hidetoshi Sano,
Yasuo Fukui
Abstract:
Studying the driving modes of turbulence is important for characterizing the impact of turbulence in various astrophysical environments. The driving mode of turbulence is parameterized by $b$, which relates the width of the gas density PDF to the turbulent Mach number; $b\approx 1/3$, $1$, and $0.4$ correspond to driving that is solenoidal, compressive, and a natural mixture of the two, respective…
▽ More
Studying the driving modes of turbulence is important for characterizing the impact of turbulence in various astrophysical environments. The driving mode of turbulence is parameterized by $b$, which relates the width of the gas density PDF to the turbulent Mach number; $b\approx 1/3$, $1$, and $0.4$ correspond to driving that is solenoidal, compressive, and a natural mixture of the two, respectively. In this work, we use high-resolution (sub-pc) ALMA $^{12}$CO ($J$ = $2-1$), $^{13}$CO ($J$ = $2-1$), and C$^{18}$O ($J$ = $2-1$) observations of filamentary molecular clouds in the star-forming region N159E (the Papillon Nebula) in the Large Magellanic Cloud (LMC) to provide the first measurement of turbulence driving parameter in an extragalactic region. We use a non-local thermodynamic equilibrium (NLTE) analysis of the CO isotopologues to construct a gas density PDF, which we find to be largely log-normal in shape with some intermittent features indicating deviations from lognormality. We find that the width of the log-normal part of the density PDF is comparable to the supersonic turbulent Mach number, resulting in $b \approx 0.9$. This implies that the driving mode of turbulence in N159E is primarily compressive. We speculate that the compressive turbulence could have been powered by gravo-turbulent fragmentation of the molecular gas, or due to compression powered by H I flows that led to the development of the molecular filaments observed by ALMA in the region. Our analysis can be easily applied to study the nature of turbulence driving in resolved star-forming regions in the local as well as the high-redshift Universe.
△ Less
Submitted 19 October, 2021; v1 submitted 8 September, 2021;
originally announced September 2021.
-
A multi-shock model for the density variance of anisotropic, highly-magnetised, supersonic turbulence
Authors:
James R. Beattie,
Philip Mocz,
Christoph Federrath,
Ralf S. Klessen
Abstract:
Shocks form the basis of our understanding for the density and velocity statistics of supersonic turbulent flows, such as those found in the cool interstellar medium (ISM). The variance of the density field, $σ^2_{ρ/ρ_0}$, is of particular interest for molecular clouds (MCs), the birthplaces of stars in the Universe. The density variance may be used to infer underlying physical processes in an MC,…
▽ More
Shocks form the basis of our understanding for the density and velocity statistics of supersonic turbulent flows, such as those found in the cool interstellar medium (ISM). The variance of the density field, $σ^2_{ρ/ρ_0}$, is of particular interest for molecular clouds (MCs), the birthplaces of stars in the Universe. The density variance may be used to infer underlying physical processes in an MC, and parameterises the star formation (SF) rate of a cloud. However, models for $σ^2_{ρ/ρ_0}$ all share a common feature -- the variance is assumed to be isotropic. This assumption does not hold when a trans/sub-Alfvénic mean magnetic field, $\vec{B}_0$, is present in the cloud, which observations suggest is relevant for some MCs. We develop an anisotropic model for $σ_{ρ/ρ_0}^2$, using contributions from hydrodynamical and fast magnetosonic shocks that propagate orthogonal to each other. Our model predicts an upper bound for $σ_{ρ/ρ_0}^2$ in the high Mach number $(\mathcal{M})$ limit as small-scale density fluctuations become suppressed by the strong $\vec{B}_0$. The model reduces to the isotropic $σ_{ρ/ρ_0}^2-\mathcal{M}$ relation in the hydrodynamical limit. To validate our model, we calculate $σ_{ρ/ρ_0}^2$ from 12~high-resolution, three-dimensional, supersonic, sub-Alfvénic magnetohydrodynamical (MHD) turbulence simulations and find good agreement with our theory. We discuss how the two MHD shocks may be the bimodally oriented over-densities observed in some MCs and the implications for SF theory in the presence of a sub-Alfvénic $\vec{B}_0$. By creating an anisotropic, supersonic density fluctuation model, this study paves the way for SF theory in the highly anisotropic regime of interstellar turbulence.
△ Less
Submitted 8 April, 2021; v1 submitted 31 January, 2021;
originally announced February 2021.
-
The sonic scale revealed by the world's largest supersonic turbulence simulation
Authors:
Christoph Federrath,
Ralf S. Klessen,
Luigi Iapichino,
James R. Beattie
Abstract:
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 t…
▽ More
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.
△ Less
Submitted 12 November, 2020;
originally announced November 2020.
-
Magnetic field fluctuations in anisotropic, supersonic turbulence
Authors:
James R. Beattie,
Christoph Federrath,
Amit Seta
Abstract:
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,…
▽ More
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.
△ Less
Submitted 28 July, 2020; v1 submitted 27 July, 2020;
originally announced July 2020.
-
Filaments and striations: anisotropies in observed, supersonic, highly-magnetised turbulent clouds
Authors:
James R. Beattie,
Christoph Federrath
Abstract:
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 spec…
▽ More
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.
△ Less
Submitted 29 November, 2019;
originally announced November 2019.
-
The relation between the turbulent Mach number and observed fractal dimensions of turbulent clouds
Authors:
James R. Beattie,
Christoph Federrath,
Ralf S. Klessen,
Nicola Schneider
Abstract:
Supersonic turbulence is a key player in controlling the structure and star formation potential of molecular clouds (MCs). The three-dimensional (3D) turbulent Mach number, $\mathcal{M}$, allows us to predict the rate of star formation. However, determining Mach numbers in observations is challenging because it requires accurate measurements of the velocity dispersion. Moreover, observations are l…
▽ More
Supersonic turbulence is a key player in controlling the structure and star formation potential of molecular clouds (MCs). The three-dimensional (3D) turbulent Mach number, $\mathcal{M}$, allows us to predict the rate of star formation. However, determining Mach numbers in observations is challenging because it requires accurate measurements of the velocity dispersion. Moreover, observations are limited to two-dimensional (2D) projections of the MCs and velocity information can usually only be obtained for the line-of-sight component. Here we present a new method that allows us to estimate $\mathcal{M}$ from the 2D column density, $Σ$, by analysing the fractal dimension, $\mathcal{D}$. We do this by computing $\mathcal{D}$ for six simulations, ranging between $1$ and $100$ in $\mathcal{M}$. From this data we are able to construct an empirical relation, $\log\mathcal{M}(\mathcal{D}) = ξ_1(\text{erfc}^{-1} [(\mathcal{D}-\mathcal{D}_{\text{min}})/Ω] + ξ_2),$ where $\text{erfc}^{-1}$ is the inverse complimentary error function, $\mathcal{D}_{\text{min}} = 1.55 \pm 0.13$ is the minimum fractal dimension of $Σ$, $Ω= 0.22 \pm 0.07$, $ξ_1 = 0.9 \pm 0.1$ and $ξ_2 = 0.2 \pm 0.2$. We test the accuracy of this new relation on column density maps from $Herschel$ observations of two quiescent subregions in the Polaris Flare MC, `saxophone' and `quiet'. We measure $\mathcal{M} \sim 10$ and $\mathcal{M} \sim 2$ for the subregions, respectively, which is similar to previous estimates based on measuring the velocity dispersion from molecular line data. These results show that this new empirical relation can provide useful estimates of the cloud kinematics, solely based upon the geometry from the column density of the cloud.
△ Less
Submitted 2 July, 2019;
originally announced July 2019.
-
The relation between the true and observed fractal dimensions of turbulent clouds
Authors:
James R. Beattie,
Christoph Federrath,
Ralf S. Klessen
Abstract:
Observations of interstellar gas clouds are typically limited to two-dimensional (2D) projections of the intrinsically three-dimensional (3D) structure of the clouds. In this study, we present a novel method for relating the 2D projected fractal dimension ($\mathcal{D}_{\text{p}}$) to the 3D fractal dimension ($\mathcal{D}_{\text{3D}}$) of turbulent clouds. We do this by computing the fractal dime…
▽ More
Observations of interstellar gas clouds are typically limited to two-dimensional (2D) projections of the intrinsically three-dimensional (3D) structure of the clouds. In this study, we present a novel method for relating the 2D projected fractal dimension ($\mathcal{D}_{\text{p}}$) to the 3D fractal dimension ($\mathcal{D}_{\text{3D}}$) of turbulent clouds. We do this by computing the fractal dimension of clouds over two orders of magnitude in turbulent Mach number $(\mathcal{M} = 1-100)$, corresponding to seven orders of magnitude in spatial scales within the clouds. This provides us with the data to create a new empirical relation between $\mathcal{D}_{\text{p}}$ and $\mathcal{D}_{\text{3D}}$. The proposed relation is $\mathcal{D}_{\text{3D}}(\mathcal{D}_{\text{p}}) = Ω_1 erfc ( ξ_1 erfc^{-1}[ (\mathcal{D}_{\text{p}} - \mathcal{D}_{\text{p,min}})/Ω_2 ] + ξ_2 ) + \mathcal{D}_{\text{3D,min}}$, where the minimum 3D fractal dimension, $\mathcal{D}_{\text{3D,min}} = 2.06 \pm 0.35$, the minimum projected fractal dimension, $\mathcal{D}_{\text{p,min}} = 1.55 \pm 0.13$, $Ω_1 = 0.47 \pm 0.18$, $Ω_2 = 0.22 \pm 0.07$, $ξ_1 = 0.80 \pm 0.18$ and $ξ_2 = 0.26 \pm 0.19$. The minimum 3D fractal dimension, $\mathcal{D}_{\text{3D,min}} = 2.06 \pm 0.35$, indicates that in the high $\mathcal{M}$ limit the 3D clouds are dominated by planar shocks. The relation between $\mathcal{D}_{\text{p}}$ and $\mathcal{D}_{\text{3D}}$ of molecular clouds may be a useful tool for those who are seeking to understand the 3D structures of molecular clouds, purely based upon 2D projected data and shows promise for relating the physics of the turbulent clouds to the fractal dimension.
△ Less
Submitted 20 May, 2019; v1 submitted 13 May, 2019;
originally announced May 2019.