Showing 1–50 of 237 results for author: Federrath, C

Wide-binary eccentricity distribution in young star clusters: dependence on the binary separation and mass

Authors: S. S. Mathew, S. Xu, C. Federrath, Y. Hu, A. Seta

Abstract: We study the wide-binary eccentricity ($e$) distribution in young star clusters and the role of turbulence in setting the form of the $e$ distribution using magnetohydrodynamical (MHD) simulations of star cluster formation. The simulations incorporate gravity, turbulence, magnetic fields, protostellar heating, and jets/outflows. We find that (1) simulations that employ purely compressive turbulenc… ▽ More We study the wide-binary eccentricity ($e$) distribution in young star clusters and the role of turbulence in setting the form of the $e$ distribution using magnetohydrodynamical (MHD) simulations of star cluster formation. The simulations incorporate gravity, turbulence, magnetic fields, protostellar heating, and jets/outflows. We find that (1) simulations that employ purely compressive turbulence driving produce binaries with a superthermal $e$ distribution ($α>1$ in $p(e) \propto e^α$), while simulations with purely solenoidal driving or natural mixture of driving modes produce subthermal/thermal distributions ($α\leq$ 1), (2) the $e$ distribution over the full range of binary separations in our simulations is set at the early stages of the star cluster formation process, (3) while binaries (separation of $r_{\mathrm{pair}} \leq 1000\, \mathrm{AU}$) have subthermal to thermal $e$ distributions ($α\sim 0.8$), wide binaries ($r_{\mathrm{pair}} > 1000\, \mathrm{AU}$) have a superthermal distribution ($α\sim 1.8$), and (4) low-mass binary systems (system masses of $M_{\mathrm{sys}} \leq 0.8\, \mathrm{M_\odot}$) have a highly superthermal distribution ($α\sim 2.4$), whereas high-mass systems ($M_{\mathrm{sys}} > 0.8\, \mathrm{M_\odot}$) exhibit a subthermal/thermal distribution ($α\sim 0.8$). The binary eccentricity distribution is often modelled as a thermal distribution. However, our results suggest that the $e$ distribution depends on the range of separation of the sampled binaries, which agrees with the findings from recent Gaia observations. We conclude that the dependence of the $e$ distribution on the binary separation and mass is linked to the binary formation mechanism governed by the turbulent properties of the parent cloud. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: 15 pages, 9 figures, 1 table (MNRAS submitted, first referee report received)

On the limitations of H alpha luminosity as a star formation tracer in spatially resolved observations

Authors: Zipeng Hu, Benjamin D. Wibking, Mark R. Krumholz, Christoph Federrath

Abstract: This study examines the limitations of H$α$ luminosity as a tracer of star formation rates (SFR) in spatially resolved observations. We carry out high-resolution simulations of a Milky Way-like galaxy including both supernova and photoionization feedback, and from these we generate synthetic H$α$ emission maps that we compare to maps of the true distribution of young stellar objects (YSOs) on scal… ▽ More This study examines the limitations of H$α$ luminosity as a tracer of star formation rates (SFR) in spatially resolved observations. We carry out high-resolution simulations of a Milky Way-like galaxy including both supernova and photoionization feedback, and from these we generate synthetic H$α$ emission maps that we compare to maps of the true distribution of young stellar objects (YSOs) on scales from whole-galaxy to individual molecular clouds ($\lesssim 100$ pc). Our results reveal significant spatial mismatches between H$α$ and true YSO maps on sub-100 pc scales, primarily due to ionizing photon leakage, with a secondary contribution from young stars drifting away from their parent molecular clouds. On small scales these effect contribute significantly to the observed anti-correlation between gas and star formation, such that there is noticeably less anti-correlation if we replace an H$α$-based star formation map with a YSO-based one; this in turn implies that previous studies have underestimated the time it takes for young stars to disperse their parent molecular clouds. However, these effects are limited in dense regions with hydrogen columns $N_\mathrm{H} > 3 \times 10^{21}$ cm$^{-2}$, where the H$α$- and YSO-based SFR maps show better agreement. Based on this finding we propose a calibration model that can precisely measure the SFR of large molecular clouds (mean radius > 100 pc) with a combination of H$α$ and CO observations, which provides a foundation for future study of star formation processes in extragalactic molecular clouds. △ Less

Submitted 16 June, 2024; originally announced June 2024.

Nexus: A framework for simulations of idealised galaxies

Authors: Thor Tepper-Garcia, Joss Bland-Hawthorn, Eugene Vasiliev, Oscar Agertz, Romain Teyssier, Christoph Federrath

Abstract: Motivated by the need for realistic, dynamically self-consistent, evolving galaxy models that avoid the inherent complexity of full, and zoom-in, cosmological simulations, we have developed Nexus, an integral, flexible framework to create synthetic galaxies made of both collisionless and gaseous components. Nexus leverages the power of publicly available, tried-and-tested packages: i) the stellar-… ▽ More Motivated by the need for realistic, dynamically self-consistent, evolving galaxy models that avoid the inherent complexity of full, and zoom-in, cosmological simulations, we have developed Nexus, an integral, flexible framework to create synthetic galaxies made of both collisionless and gaseous components. Nexus leverages the power of publicly available, tried-and-tested packages: i) the stellar-dynamics, action-based library AGAMA; and ii) the Adaptive Mesh Refinement, N-body/hydrodynamical code Ramses, modified to meet our needs. In addition, we make use of a proprietary module to account for realistic galaxy formation (sub-grid) physics, including star formation, stellar feedback, and chemical enrichment. As a framework to perform controlled experiments of idealised galaxies, Nexus' basic functionality consists in the generation of bespoke initial conditions (ICs) for any desired galaxy model, which are advanced in time to simulate the system's evolution. The fully self-consistent ICs are generated with a distribution-function based approach, as implemented in the galaxy modelling module of AGAMA -- up to now restricted to collisionless components, extended in this work to treat two types of gaseous configurations: (i) hot halos; and (ii) gas disks. For the first time, we are able to construct equilibrium models with disc gas fractions in the range $0 < f_{\rm gas} < 1$, needed to model high-redshift galaxies. The framework is ideally suited to the study of galactic ecology, specifically how stars and gas work together over billions of years. As a validation of our framework, we reproduce - and improve upon - several isolated galaxy model setups reported in earlier studies. Finally, we showcase Nexus by presenting an interesting type of `nested bar' galaxy class. Future upgrades of Nexus will include magneto-hydrodynamics and highly energetic particle (`cosmic ray') heating. △ Less

Submitted 1 June, 2024; originally announced June 2024.

Comments: To be submitted to MNRAS; we welcome general comments and pointers to relevant, missing references

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.

Comments: 22 pages. 10 figures. Currently the world's largest MHD turbulence simulation on a grid of 10,080^3 (Re and Rm over a million). Comments welcome, email james.beattie@princeton.edu

Probing the role of self-gravity in clouds impacted by AGN-driven winds

Authors: Ankush Mandal, Dipanjan Mukherjee, Christoph Federrath, Geoffrey V. Bicknell, Nicole P. H. Nesvadba, Andrea Mignone

Abstract: The impact of winds and jet-inflated bubbles driven by active galactic nuclei (AGN) are believed to significantly affect the host galaxy's interstellar medium (ISM) and regulate star formation. To explore this scenario, we perform a suite of hydrodynamic simulations to model the interaction between turbulent star-forming clouds and highly pressurised AGN-driven outflows, focusing on the effects of… ▽ More The impact of winds and jet-inflated bubbles driven by active galactic nuclei (AGN) are believed to significantly affect the host galaxy's interstellar medium (ISM) and regulate star formation. To explore this scenario, we perform a suite of hydrodynamic simulations to model the interaction between turbulent star-forming clouds and highly pressurised AGN-driven outflows, focusing on the effects of self-gravity. Our results demonstrate that the cloudlets fragmented by the wind can become gravitationally bound, significantly increasing their survival time. While external pressurisation leads to a global collapse of the clouds in cases of weaker winds ($10^{42}-10^{43}~{\rm erg~s^{-1}}$), higher-power winds ($10^{44}-10^{45}~{\rm erg~s^{-1}}$) disperse the gas and cause localised collapse of the cloudlets. We also demonstrate that a kinetic energy-dominated wind is more efficient in accelerating and dispersing the gas than a thermal wind with the same power. The interaction can give rise to multi-phase outflows with velocities ranging from a few 100 to several 1000~${\rm km\,s^{-1}}$. The mass outflow rates are tightly correlated with the wind power, which we explain by an ablation-based mass-loss model. Moreover, the velocity dispersion and the virial parameter of the cloud material can increase by up to one order of magnitude through the effect of the wind. Even though the wind can suppress or quench star formation for about 1 Myr during the initial interaction, a substantial number of gravitationally bound dense cloudlets manage to shield themselves from the wind's influence and subsequently undergo rapid gravitational collapse, leading to an enhanced star formation rate (SFR). △ Less

Submitted 16 May, 2024; originally announced May 2024.

Comments: Accepted for publication in MNRAS. Accepted 2024 May 16. Received 2024 May 16; in original form 2024 January 18

Turbulence statistics of HI clouds entrained in the Milky Way's nuclear wind

Authors: Isabella A. Gerrard, Karlie A. Noon, Christoph Federrath, Enrico M. Di Teodoro, Antoine Marchal, N. M. McClure-Griffiths

Abstract: The interstellar medium (ISM) is ubiquitously turbulent across many physically distinct environments within the Galaxy. Turbulence is key in controlling the structure and dynamics of the ISM, regulating star formation, and transporting metals within the Galaxy. We present the first observational measurements of turbulence in neutral hydrogen entrained in the hot nuclear wind of the Milky Way. Usin… ▽ More The interstellar medium (ISM) is ubiquitously turbulent across many physically distinct environments within the Galaxy. Turbulence is key in controlling the structure and dynamics of the ISM, regulating star formation, and transporting metals within the Galaxy. We present the first observational measurements of turbulence in neutral hydrogen entrained in the hot nuclear wind of the Milky Way. Using recent MeerKAT observations of two extra-planar HI clouds above (gal. lat.$\,\sim7.0^{\circ}$) and below (gal. lat.$\,\sim-3.9^{\circ}$) the Galactic disc, we analyse centroid velocity and column density maps to estimate the velocity dispersion ($σ_{v,\mathrm{3D}}$), the turbulent sonic Mach number ($\mathcal{M}$), the volume density dispersion ($σ_{ρ/ρ_0}$), and the turbulence driving parameter ($b$). We also present a new prescription for estimating the spatial temperature variations of HI in the presence of related molecular gas. We measure these turbulence quantities on the global scale of each cloud, but also spatially map their variation across the plane-of-sky extent of each cloud by using a roving kernel method. We find that the two clouds share very similar characteristics of their internal turbulence, despite their varying latitudes. Both clouds are in the sub-to-trans-sonic Mach regime, and have primarily compressively-driven ($b\sim1$) turbulence. Given that there is no known active star-formation present in these clouds, this may be indicative of the way the cloud-wind interaction injects energy into the entrained atomic material on parsec scales. △ Less

Submitted 28 April, 2024; originally announced April 2024.

Comments: Accepted to MNRAS, 15 pages, 13 figures

Protostellar spin-up and fast rotator formation through binary star formation

Authors: Rajika L. Kuruwita, Christoph Federrath, Marina Kounkel

Abstract: (Edited) Many fast rotator stars (rotation periods of < 2 days) are found in unresolved binaries with separations of tens of au. This correlation leads to the question of whether the formation of binary stars inherently produces fast rotators. We aim to understand whether the formation of companions plays a role in spinning up stars. We use magneto-hydrodynamical simulations to study the formation… ▽ More (Edited) Many fast rotator stars (rotation periods of < 2 days) are found in unresolved binaries with separations of tens of au. This correlation leads to the question of whether the formation of binary stars inherently produces fast rotators. We aim to understand whether the formation of companions plays a role in spinning up stars. We use magneto-hydrodynamical simulations to study the formation of multiple star systems from turbulent and non-turbulent protostellar cores. We track the angular momentum accreted by individual star and inner disc systems by using a sink particle technique. We run a resolution study to extrapolate protostellar properties. We find in all simulations that the primary star can experience spin-up events that are correlated with the formation of companions. The primary star can spin up by up to 84% of its pre-fragmentation angular momentum and by up to 18% of its pre-fragmentation mass-specific angular momentum. The mechanism for the spin-up is gravitational disc instabilities in the circumstellar disc around the primary star, leading to the accretion of material with high specific angular momentum. The simulations that experience the strongest disc instabilities fragment to form companions. Simulations with weaker spin-up events experience disc instabilities triggered by a companion flyby, and the disc instability in these cases does not produce further fragments. We conclude that the primary star in multiple star systems may end up with a higher spin than single stars. This is because gravitational instabilities in the circumstellar disc around the primary star can trigger a spin-up event. In the strongest spin-up events, the instability is likely to cause disc fragmentation and the formation of companions. This companion formation coupled with shorter disc lifetimes, because the companion truncates the circumstellar disc, can help produce fast rotators. △ Less

Submitted 22 April, 2024; originally announced April 2024.

Comments: 14 pages, 11 Figures. Submitted to A&A

Magnetic field morphology and evolution in the Central Molecular Zone and its effect on gas dynamics

Authors: R. G. Tress, M. C. Sormani, P. Girichidis, S. C. O. Glover, R. S. Klessen, R. J. Smith, E. Sobacchi, L. Armillotta, A. T. Barnes, C. Battersby, K. R. J. Bogue, N. Brucy, L. Colzi, C. Federrath, P. García, A. Ginsburg, J. Göller, H P. Hatchfield, C. Henkel, P. Hennebelle, J. D. Henshaw, M. Hirschmann, Y. Hu, J. Kauffmann, J. M. D. Kruijssen , et al. (12 additional authors not shown)

Abstract: The interstellar medium in the Milky Way's Central Molecular Zone (CMZ) is known to be strongly magnetised, but its large-scale morphology and impact on the gas dynamics are not well understood. We explore the impact and properties of magnetic fields in the CMZ using three-dimensional non-self gravitating magnetohydrodynamical simulations of gas flow in an external Milky Way barred potential. We f… ▽ More The interstellar medium in the Milky Way's Central Molecular Zone (CMZ) is known to be strongly magnetised, but its large-scale morphology and impact on the gas dynamics are not well understood. We explore the impact and properties of magnetic fields in the CMZ using three-dimensional non-self gravitating magnetohydrodynamical simulations of gas flow in an external Milky Way barred potential. We find that: (1) The magnetic field is conveniently decomposed into a regular time-averaged component and an irregular turbulent component. The regular component aligns well with the velocity vectors of the gas everywhere, including within the bar lanes. (2) The field geometry transitions from parallel to the Galactic plane near $z=0$ to poloidal away from the plane. (3) The magneto-rotational instability (MRI) causes an in-plane inflow of matter from the CMZ gas ring towards the central few parsecs of $0.01-0.1$ M$_\odot$ yr$^{-1}$ that is absent in the unmagnetised simulations. However, the magnetic fields have no significant effect on the larger-scale bar-driven inflow that brings the gas from the Galactic disc into the CMZ. (4) A combination of bar inflow and MRI-driven turbulence can sustain a turbulent vertical velocity dispersion of $σ_z \simeq 5$ km s$^{-1}$ on scales of $20$ pc in the CMZ ring. The MRI alone sustains a velocity dispersion of $σ_z \simeq 3$ km s$^{-1}$. Both these numbers are lower than the observed velocity dispersion of gas in the CMZ, suggesting that other processes such as stellar feedback are necessary to explain the observations. (5) Dynamo action driven by differential rotation and the MRI amplifies the magnetic fields in the CMZ ring until they saturate at a value that scales with the average local density as $B \simeq 102 (n/10^3 cm^{-3})^{0.33}$ $μ$G. Finally, we discuss the implications of our results within the observational context in the CMZ. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Comments: submitted to A&A, comments welcome

Turbulent gas-rich disks at high redshift: bars & bulges in a radial shear flow

Authors: Joss Bland-Hawthorn, Thor Tepper-Garcia, Oscar Agertz, Christoph Federrath

Abstract: Recent observations of high-redshift galaxies ($z \lesssim 7$) reveal that a substantial fraction have turbulent, gas-rich disks with well-ordered rotation and elevated levels of star formation. In some instances, disks show evidence of spiral arms, with bar-like structures. These remarkable observations have encouraged us to explore a new class of dynamically self-consistent models using our hydr… ▽ More Recent observations of high-redshift galaxies ($z \lesssim 7$) reveal that a substantial fraction have turbulent, gas-rich disks with well-ordered rotation and elevated levels of star formation. In some instances, disks show evidence of spiral arms, with bar-like structures. These remarkable observations have encouraged us to explore a new class of dynamically self-consistent models using our hydrodynamic N-body simulation framework that mimic a plausible progenitor of the Milky Way at high redshift. We explore disk gas fractions of $f_{\rm gas} = 0, 20, 40, 60, 80, 100\%$ and track the creation of stars and metals. The high gas surface densities encourage vigorous star formation, which in turn couples with the gas to drive turbulence. We explore three distinct histories: (i) there is no ongoing accretion and the gas is used up by the star formation; (ii) the star-forming gas is replenished by cooling in the hot halo gas; (iii) in a companion paper, we revisit the models in the presence of a strong perturbing force. At low $f_{\rm disk}$ ($<0.3$), where $f_{\rm disk}$ is the mass fraction of stars relative to dark matter within 2.2 $R_{\rm disk}$, a bar does not form in a stellar disk; this remains true even when gas dominates the inner disk potential. For a dominant baryon disk ($f_{\rm disk} \gtrsim 0.5$) at all gas fractions, the turbulent gas forms a strong "radial shear flow" that leads to an intermittent star-forming bar within about 500 Myr; turbulent gas speeds up the formation of bars compared to gas-free models. For $f_{\rm gas} \lesssim 60\%$, all bars survive, but for higher gas fractions, the bar devolves into a central bulge after 1 Gyr. The star-forming bars are reminiscent of recent discoveries in high-redshift ALMA observations of gaseous disks. △ Less

Submitted 1 May, 2024; v1 submitted 8 February, 2024; originally announced February 2024.

Comments: 40 pages, 21 figures, ApJ, in press (April 2024) - simulation movies at http://www.physics.usyd.edu.au/turbo_disks

Amplification and saturation of turbulent magnetic field in collapsing primordial gas clouds

Authors: Sho Higashi, Hajime Susa, Christoph Federrath, Gen Chiaki

Abstract: Recent numerical studies suggest that magnetic fields play an important role in primordial star formation in the early universe. However, the detailed evolution of the magnetic field in the collapse phase still has uncertainties because of the complicated physics associated with turbulence in a collapsing magnetized system. Here, we perform a suite of numerical MHD simulations that follow the coll… ▽ More Recent numerical studies suggest that magnetic fields play an important role in primordial star formation in the early universe. However, the detailed evolution of the magnetic field in the collapse phase still has uncertainties because of the complicated physics associated with turbulence in a collapsing magnetized system. Here, we perform a suite of numerical MHD simulations that follow the collapse of magnetized, turbulent primordial gas clouds to investigate the evolution of the magnetic field associated with the turbulence, assuming a polytropic equation of state with exponent $γ_{\rm eff}$ and with various numerical resolutions. In addition, we generalize the analytic theory of magnetic field growth/saturation so that it can deal with various exponents $γ_{\rm eff}$ and turbulence energy spectra. We find that the numerical results are well reproduced by the theory for various $γ_{\rm eff}$ through the collapse phase during the formation of the first stars. The magnetic field is eventually amplified by a factor of $10^{12}$ -- $10^{15}$ due to kinematic and non-linear turbulent dynamo effects and reaches 3% -- 100% of the equipartition level, depending on $γ_{\rm eff}$. We also find that the transition between the kinematic and non-linear stages can be analytically estimated. These results indicate that the strong magnetic field accompanied by supersonic turbulence is a general property and suggest that it can play a crucial role in the formation of the first stars. △ Less

Submitted 18 January, 2024; originally announced January 2024.

Comments: 15 pages, 8 figures, 1 tables, accepted for publication in ApJ

Critical magnetic Reynolds number of the turbulent dynamo in collisionless plasmas

Authors: Radhika Achikanath Chirakkara, Amit Seta, Christoph Federrath, Matthew W. Kunz

Abstract: The intracluster medium of galaxy clusters is an extremely hot and diffuse, nearly collisionless plasma, which hosts dynamically important magnetic fields of $\sim μ{\rm G}$ strength. Seed magnetic fields of much weaker strength of astrophysical or primordial origin can be present in the intracluster medium. In collisional plasmas, which can be approximated in the magneto-hydrodynamical (MHD) limi… ▽ More The intracluster medium of galaxy clusters is an extremely hot and diffuse, nearly collisionless plasma, which hosts dynamically important magnetic fields of $\sim μ{\rm G}$ strength. Seed magnetic fields of much weaker strength of astrophysical or primordial origin can be present in the intracluster medium. In collisional plasmas, which can be approximated in the magneto-hydrodynamical (MHD) limit, the turbulent dynamo mechanism can amplify weak seed fields to strong dynamical levels efficiently by converting turbulent kinetic energy into magnetic energy. However, the viability of this mechanism in weakly collisional or completely collisionless plasma is much less understood. In this study, we explore the properties of the collisionless turbulent dynamo by using three-dimensional hybrid-kinetic particle-in-cell simulations. We explore the properties of the collisionless turbulent dynamo in the kinematic regime for different values of the magnetic Reynolds number, ${\rm Rm}$, initial magnetic-to-kinetic energy ratio, $(E_{\rm{mag}}/E_{\rm{kin}})_{\rm{i}}$, and initial Larmor ratio, $(r_{\rm Larmor}/L_{\rm box})_{\rm i}$, i.e., the ratio of the Larmor radius to the size of the turbulent system. We find that in the `un-magnetised' regime, $(r_{\rm Larmor}/L_{\rm box})_{\rm i} > 1$, the critical magnetic Reynolds number for the dynamo action ${\rm Rm}_{\rm crit} \approx 107 \pm 3$. In the `magnetised' regime, $(r_{\rm Larmor}/L_{\rm box})_{\rm i} \lesssim 1$, we find a marginally higher ${\rm Rm}_{\rm crit} = 124 \pm 8$. We find that the growth rate of the magnetic energy does not depend on the strength of the seed magnetic field when the initial magnetisation is fixed. We also study the distribution and evolution of the pressure anisotropy in the collisionless plasma and compare our results with the MHD turbulent dynamo. △ Less

Submitted 16 January, 2024; originally announced January 2024.

Comments: 17 pages, 16 figures

Journal ref: Monthly Notices of the Royal Astronomical Society (2023)

A Galactic Eclipse: The Small Magellanic Cloud is Forming Stars in Two, Superimposed Systems

Authors: Claire E. Murray, Sten Hasselquist, Joshua E. G. Peek, Christina Willecke Lindberg, Andres Almeida, Yumi Choi, Jessica E. M. Craig, Helga Denes, John M. Dickey, Enrico M. Di Teodoro, Christoph Federrath, Isabella A. Gerrard, Steven J. Gibson, Denis Leahy, Min-Young Lee, Callum Lynn, Yik Ki Ma, Antoine Marchal, N. M. McClure-Griffiths, David Nidever, Hiep Nguyen, Nickolas M. Pingel, Elizabeth Tarantino, Lucero Uscanga, Jacco Th. van Loon

Abstract: The structure and dynamics of the star-forming disk of the Small Magellanic Cloud (SMC) have long confounded us. The SMC is widely used as a prototype for galactic physics at low metallicity, and yet we fundamentally lack an understanding of the structure of its interstellar medium (ISM). In this work, we present a new model for the SMC by comparing the kinematics of young, massive stars with the… ▽ More The structure and dynamics of the star-forming disk of the Small Magellanic Cloud (SMC) have long confounded us. The SMC is widely used as a prototype for galactic physics at low metallicity, and yet we fundamentally lack an understanding of the structure of its interstellar medium (ISM). In this work, we present a new model for the SMC by comparing the kinematics of young, massive stars with the structure of the ISM traced by high-resolution observations of neutral atomic hydrogen (HI) from the Galactic Australian Square Kilometer Array Pathfinder survey (GASKAP-HI). Specifically, we identify thousands of young, massive stars with precise radial velocity constraints from the Gaia and APOGEE surveys and match these stars to the ISM structures in which they likely formed. By comparing the average dust extinction towards these stars, we find evidence that the SMC is composed of two structures with distinct stellar and gaseous chemical compositions. We construct a simple model that successfully reproduces the observations and shows that the ISM of the SMC is arranged into two, superimposed, star-forming systems with similar gas mass separated by ~5 kpc along the line of sight. △ Less

Submitted 12 December, 2023; originally announced December 2023.

Comments: ApJ accepted. 20 pages, 18 figures

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.

Comments: 28 pages, 18 figures, submitted to MNRAS. Comments welcome. Re-uploaded on 11 Dec with compressed large figures to reduce file size

The role of magnetic fields in disc galaxies: spiral arm instability

Authors: Raghav Arora, Christoph Federrath, Robi Banerjee, Bastian Körtgen

Abstract: Context. Regularly-spaced, star-forming regions along the spiral arms of nearby galaxies provide insight into the early stages and initial conditions of star formation. The regular separation of these star-forming regions suggests spiral arm instability as their origin. Aims. We explore the effects of magnetic fields on the spiral arm instability. Methods. We use three-dimensional global magne… ▽ More Context. Regularly-spaced, star-forming regions along the spiral arms of nearby galaxies provide insight into the early stages and initial conditions of star formation. The regular separation of these star-forming regions suggests spiral arm instability as their origin. Aims. We explore the effects of magnetic fields on the spiral arm instability. Methods. We use three-dimensional global magnetohydrodynamical simulations of isolated spiral galaxies, comparing three different initial plasma $β$ values (ratios of thermal to magnetic pressure) of $β=\infty$, $50$, and $10$. We perform Fourier analysis to calculate the separation of the over-dense regions formed from the spiral instability. We then compare the separations with observations. Results. We find that the spiral arms in the hydro case ($β= \infty$) are unstable, with the fragments initially connected by gas streams, reminiscent of Kelvin-Helmholtz instability. In the $β= 50$ case, the spiral arms fragment, but the fragments separate earlier and tend to be elongated in the direction perpendicular to the spiral arms. However, in the $β= 10$ run the arms are stabilised against fragmentation by magnetic pressure. The spiral arms in the unstable cases fragment into regularly-spaced, over-dense regions. We determine their separation to be $\sim 0.5$ kpc in the hydro and $\sim 0.65$ kpc in the $β= 50$ case, both in agreement with the observations of nearby galaxies. We find a smaller median characteristic wavelength of the over-densities to be $0.73^{+0.31}_{-0.36}$ kpc in the hydro case, compared to $0.98^{+0.49}_{-0.46}$ kpc in the $β= 50$ case. Moreover, we find a higher growth rate of the over-densities in the $β= 50$ run compared to the hydro run. We observe magnetic hills and valleys along the fragmented arms in the $β= 50$ run, which is characteristic of the Parker instability. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: 14 pages, 10 figures, submitted to Astronomy & Astrophysics

Numerical viscosity and resistivity in MHD turbulence simulations

Authors: Lakshmi Malvadi Shivakumar, Christoph Federrath

Abstract: To ensure that magnetohydrodynamical (MHD) turbulence simulations accurately reflect the physics, it is critical to understand numerical dissipation. Here we determine the hydrodynamic and magnetic Reynolds number (Re and Rm) as a function of linear grid resolution N, in MHD simulations with purely numerical viscosity and resistivity (implicit large eddy simulations; ILES). We quantify the numeric… ▽ More To ensure that magnetohydrodynamical (MHD) turbulence simulations accurately reflect the physics, it is critical to understand numerical dissipation. Here we determine the hydrodynamic and magnetic Reynolds number (Re and Rm) as a function of linear grid resolution N, in MHD simulations with purely numerical viscosity and resistivity (implicit large eddy simulations; ILES). We quantify the numerical viscosity in the subsonic and supersonic regime, via simulations with sonic Mach numbers of Mach=0.1 and Mach=10, respectively. We find Re=(N/N_Re)^p_Re, with p_Re=[1.2,1.4] and N_Re=[0.8,1.7] for Mach=0.1, and p_Re=[1.5,2.0] and N_Re=[0.8,4.4] for Mach=10, and we find Rm=(N/N_Rm)^p_Rm, with p_Rm=[1.3,1.5] and N_Rm=[1.1,2.3] for Mach=0.1, and p_Rm=[1.2,1.6] and N_Rm=[0.1,0.7] for Mach=10. The resulting magnetic Prandtl number (Pm=Rm/Re) is consistent with a constant value of Pm=1.3+/-1.1 for Mach=0.1, and Pm=2.0+/-1.4 for Mach=10. We compare our results with an independent study in the subsonic regime and find excellent agreement in p_Re and p_Rm, and agreement within a factor of ~2 for N_Re and N_Rm (due to differences in the codes and solvers). We compare these results to the target Re and Rm set in direct numerical simulations (DNS, i.e., using explicit viscosity and resistivity) from the literature. This comparison and our ILES relations can be used to determine whether a target Re and Rm can be achieved in a DNS for a given N. We conclude that for the explicit (physical) dissipation to dominate over the numerical dissipation, the target Reynolds numbers must be set lower than the corresponding numerical values. △ Less

Submitted 17 November, 2023; originally announced November 2023.

Comments: 19 pages, 9 figures, submitted to MNRAS, comments welcome

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.

Comments: 25 pages, 15 figures, submitted to MNRAS, json-file w/ data

A new method for spatially resolving the turbulence driving mixture in the ISM with application to the Small Magellanic Cloud

Authors: Isabella A. Gerrard, Christoph Federrath, Nickolas M. Pingel, Naomi M. McClure-Griffiths, Antoine Marchal, Gilles Joncas, Susan E. Clark, Snežana Stanimirović, Min-Young Lee, Jacco Th. van Loon, John Dickey, Helga Dénes, Yik Ki Ma, James Dempsey, Callum Lynn

Abstract: Turbulence plays a crucial role in shaping the structure of the interstellar medium. The ratio of the three-dimensional density contrast ($σ_{ρ/ρ_0}$) to the turbulent sonic Mach number ($\mathcal{M}$) of an isothermal, compressible gas describes the ratio of solenoidal to compressive modes in the turbulent acceleration field of the gas, and is parameterised by the turbulence driving parameter:… ▽ More Turbulence plays a crucial role in shaping the structure of the interstellar medium. The ratio of the three-dimensional density contrast ($σ_{ρ/ρ_0}$) to the turbulent sonic Mach number ($\mathcal{M}$) of an isothermal, compressible gas describes the ratio of solenoidal to compressive modes in the turbulent acceleration field of the gas, and is parameterised by the turbulence driving parameter: $b=σ_{ρ/ρ_0}/\mathcal{M}$. The turbulence driving parameter ranges from $b=1/3$ (purely solenoidal) to $b=1$ (purely compressive), with $b=0.38$ characterising the natural mixture (1/3~compressive, 2/3~solenoidal) of the two driving modes. Here we present a new method for recovering $σ_{ρ/ρ_0}$, $\mathcal{M}$, and $b$, from observations on galactic scales, using a roving kernel to produce maps of these quantities from column density and centroid velocity maps. We apply our method to high-resolution HI emission observations of the Small Magellanic Cloud (SMC) from the GASKAP-HI survey. We find that the turbulence driving parameter varies between $b\sim 0.3$ and $b\sim 1.0$ within the main body of the SMC, but the median value converges to $b\sim0.51$, suggesting that the turbulence is overall driven more compressively ($b>0.38$). We observe no correlation between the $b$ parameter and HI or H$α$ intensity, indicating that compressive driving of HI turbulence cannot be determined solely by observing HI or H$α$ emission density, and that velocity information must also be considered. Further investigation is required to link our findings to potential driving mechanisms such as star-formation feedback, gravitational collapse, or cloud-cloud collisions. △ Less

Submitted 19 September, 2023; originally announced September 2023.

Comments: 20 pages, 16 figures, accepted to MNRAS

Multiphase Neutral Interstellar Medium: Analyzing Simulation with H I 21cm Observational Data Analysis Techniques

Authors: Soumyadeep Bhattacharjee, Nirupam Roy, Prateek Sharma, Amit Seta, Christoph Federrath

Abstract: Several different methods are regularly used to infer the properties of the neutral interstellar medium (ISM) using atomic hydrogen (H I) 21cm absorption and emission spectra. In this work, we study various techniques used for inferring ISM gas phase properties, namely the correlation between brightness temperature and optical depth $(T_B(v)$, $τ(v))$ at each channel velocity ($v$), and decomposit… ▽ More Several different methods are regularly used to infer the properties of the neutral interstellar medium (ISM) using atomic hydrogen (H I) 21cm absorption and emission spectra. In this work, we study various techniques used for inferring ISM gas phase properties, namely the correlation between brightness temperature and optical depth $(T_B(v)$, $τ(v))$ at each channel velocity ($v$), and decomposition into Gaussian components, by creating mock spectra from a 3D magnetohydrodynamic simulation of a two-phase, turbulent ISM. We propose a physically motivated model to explain the $T_B(v)-τ(v)$ distribution and relate the model parameters to properties like warm gas spin temperature and cold cloud length scales. Two methods based on Gaussian decomposition -- using only absorption spectra and both absorption and emission spectra -- are used to infer the column density distribution as a function of temperature. In observations, such analysis reveals the puzzle of large amounts (significantly higher than in simulations) of gas with temperature in the thermally unstable range of $\sim$200 K to $\sim$2000 K and a lack of the expected bimodal (two-phase) temperature distribution. We show that, in simulation, both methods are able to recover the actual gas distribution in the simulation till temperatures $\lesssim2500$~K (and the two-phase distribution in general) reasonably well. We find our results to be robust to a range of effects such as noise, varying emission beam size, and simulation resolution. This shows that the observational inferences are unlikely to be artifacts, thus highlighting a tension between observations and simulations. We discuss possible reasons for this tension and ways to resolve it. △ Less

Submitted 23 November, 2023; v1 submitted 10 September, 2023; originally announced September 2023.

Comments: 22 pages (including appendixes), 16 figures, 3 tables, Accepted for publication in MNRAS

Structure of iso-density sets in supersonic isothermal turbulence

Authors: Fabien Thiesset, Christoph Federrath

Abstract: Context: The gas density structure of the cold molecular phase of the interstellar medium is the main controller of star formation. Aims: A theoretical framework is proposed to describe the structural content of the density field in isothermal supersonic turbulence. Methods: It makes use of correlation and structure functions of the phase indicator field defined for different iso-density values. T… ▽ More Context: The gas density structure of the cold molecular phase of the interstellar medium is the main controller of star formation. Aims: A theoretical framework is proposed to describe the structural content of the density field in isothermal supersonic turbulence. Methods: It makes use of correlation and structure functions of the phase indicator field defined for different iso-density values. The relations between these two-point statistics and the geometrical features of iso-density sets such as the volume fraction, the surface density, the curvature, and fractal characteristics are provided. An exact scale-by-scale budget equation is further derived revealing the role of the turbulent cascade and dilation on the structural evolution of the density field. Although applicable to many flow situations, this tool is here first invoked for characterising supersonic isothermal turbulence, using data from the currently best-resolved numerical simulation. Results: We show that iso-density sets are surface fractals rather than mass fractals, with dimensions that markedly differ between dilute, neutral, and dense regions. The surface-size relation is established for different iso-density values. We further find that the turbulent cascade of iso-density sets is directed from large towards smaller scales, in agreement with the classical picture that turbulence acts to concentrate more surface into smaller volumes. Intriguingly, there is no range of scales that complies with a constant transfer rate in the cascade, challenging our fundamental understanding of interstellar turbulence. Finally, we recast the virial theorem in a new formulation drawing an explicit relation between the aforementioned geometrical measures and the dynamics of iso-density sets. △ Less

Submitted 14 August, 2023; originally announced August 2023.

Comments: 15 pages, 14 figures

Journal ref: Astronomy and Astrophysics, 676, A12 (2023)

The Rapid ASKAP Continuum Survey III: Spectra and Polarisation In Cutouts of Extragalactic Sources (SPICE-RACS) First Data Release

Authors: Alec J. M. Thomson, David McConnell, Emil Lenc, Timothy J Galvin, Lawrence Rudnick, George Heald, Catherine L. Hale, Stefan W. Duchesne, Craig S. Anderson, Ettore Carretti, Christoph Federrath, B. M. Gaensler, Lisa Harvey-Smith, Marijke Haverkorn, Aidan W. Hotan, Yik Ki Ma, Tara Murphy, N. M. McClure-Griffith, Vanessa A. Moss, Shane P. O'Sullivan, Wasim Raja, Amit Seta, Cameron L. Van Eck, Jennifer L. West, Matthew T. Whiting , et al. (1 additional authors not shown)

Abstract: The Australian SKA Pathfinder (ASKAP) radio telescope has carried out a survey of the entire Southern Sky at 887.5MHz. The wide area, high angular resolution, and broad bandwidth provided by the low-band Rapid ASKAP Continuum Survey (RACS-low) allow the production of a next-generation rotation measure (RM) grid across the entire Southern Sky. Here we introduce this project as Spectral and Polarisa… ▽ More The Australian SKA Pathfinder (ASKAP) radio telescope has carried out a survey of the entire Southern Sky at 887.5MHz. The wide area, high angular resolution, and broad bandwidth provided by the low-band Rapid ASKAP Continuum Survey (RACS-low) allow the production of a next-generation rotation measure (RM) grid across the entire Southern Sky. Here we introduce this project as Spectral and Polarisation in Cutouts of Extragalactic sources from RACS (SPICE-RACS). In our first data release, we image 30 RACS-low fields in Stokes $I$, $Q$, $U$ at 25'' angular resolution, across 744 to 1032MHz with 1MHz spectral resolution. Using a bespoke, highly parallelised, software pipeline we are able to rapidly process wide-area spectro-polarimetric ASKAP observations. Notably, we use 'postage stamp' cutouts to assess the polarisation properties of \ncomponents\ radio components detected in total intensity. We find that our Stokes $Q$ and $U$ images have an rms noise of ~80$μ$Jy/PSF, and our correction for instrumental polarisation leakage allows us to characterise components with >1% polarisation fraction over most of the field of view. We produce a broadband polarised radio component catalogue that contains \nrms\ RM measurements over an area of ~1300deg^2 with an average error in RM of 1.6+1.1-1.0rad/m^2, and an average linear polarisation fraction 3.4+3.0-1.6%. We determine this subset of components using the conditions that the polarised signal-to-noise ratio is $>8$, the polarisation fraction is above our estimated polarised leakage, and the Stokes $I$ spectrum has a reliable model. Our catalogue provides an areal density of $4\pm2$ RMs/deg^2; an increase of $\sim4$ times over the previous state-of-the-art (Taylor et al. 2009). Meaning that, having used just 3% of the RACS-low sky area, we have produced the 3rd largest RM catalogue to date. This catalogue has broad applications for studying... △ Less

Submitted 14 July, 2023; originally announced July 2023.

Comments: 42 pages, 24 figures, 6 tables. Accepted for publications in PASA

Measuring the hot ICM velocity structure function using XMM-Newton observations

Authors: Efrain Gatuzz, R. Mohapatra, C. Federrath, J. S. Sanders, A. Liu, S. A. Walker, C. Pinto

Abstract: It has been shown that the gas velocities within the intracluster medium (ICM) can be measured by applying novel XMM-Newton EPIC-pn energy scale calibration, which uses instrumental Cu Ka as reference for the line emission. Using this technique, we have measured the velocity distribution of the ICM for clusters involving AGN feedback and sloshing of the plasma within the gravitational well (Virgo… ▽ More It has been shown that the gas velocities within the intracluster medium (ICM) can be measured by applying novel XMM-Newton EPIC-pn energy scale calibration, which uses instrumental Cu Ka as reference for the line emission. Using this technique, we have measured the velocity distribution of the ICM for clusters involving AGN feedback and sloshing of the plasma within the gravitational well (Virgo and Centaurus) and a relaxed one (Ophiuchus). We present a detailed study of the kinematics of the hot ICM for these systems. First, we compute the velocity probability distribution functions (PDFs) from the velocity maps. We find that for all sources the PDF follows a normal distribution, with a hint for a multimodal distribution in the case of Ophiuchus. Then, we compute the velocity structure function (VSF) for all sources in order to study the variation with scale as well as the nature of turbulence in the ICM. We measure a turbulence driving scale of 10-20 kpc for the Virgo cluster, while the Ophiuchus cluster VSF reflects the absence of strong interaction between the ICM and a powerful Active Galactic Nucleus (AGN) at such spatial scales. For the former, we compute a dissipation time larger than the jet activity cycle, thus indicating that a more efficient heating process than turbulence is required to reach equilibrium. This is the first time that the VSF of the hot ICM has been computed using direct velocity measurements from X-ray astronomical observations. △ Less

Submitted 5 July, 2023; originally announced July 2023.

Comments: 10 pages, 9 Figures

Fitting Probability Distribution Functions in Turbulent Star-Forming Molecular Clouds

Authors: Avery Kiihne, Sabrina M. Appel, Blakesley Burkhart, Vadim A. Semenov, Christoph Federrath

Abstract: We use a suite of 3D simulations of star-forming molecular clouds, with and without stellar feedback and magnetic fields, to investigate the effectiveness of different fitting methods for volume and column density probability distribution functions (PDFs). The first method fits a piecewise lognormal and power-law (PL) function to recover PDF parameters such as the PL slope and transition density.… ▽ More We use a suite of 3D simulations of star-forming molecular clouds, with and without stellar feedback and magnetic fields, to investigate the effectiveness of different fitting methods for volume and column density probability distribution functions (PDFs). The first method fits a piecewise lognormal and power-law (PL) function to recover PDF parameters such as the PL slope and transition density. The second method fits a polynomial spline function and examines the first and second derivatives of the spline to determine the PL slope and the functional transition density. We demonstrate that fitting a spline allows us to directly determine if the data has multiple PL slopes. The first PL (set by the transition between lognormal and PL function) can also be visualized in the derivatives directly. In general, the two methods produce fits that agree reasonably well for volume density but vary for column density, likely due to the increased statistical noise in column density maps as compared to volume density. We test a well-known conversion for estimating volume density PL slopes from column density slopes and find that the spline method produces a better match (\c{hi}2 of 2.38 vs \c{hi}2 of 5.92), albeit with a significant scatter. Ultimately, we recommend the use of both fitting methods on column density data to mitigate the effects of noise. △ Less

Submitted 18 May, 2023; originally announced May 2023.

Comments: 14 pages, 8 figures, Submitted to ApJ. Comments welcome

Kinematics of Galactic Centre clouds shaped by shear-seeded solenoidal turbulence

Authors: Maya A. Petkova, J. M. Diederik Kruijssen, Jonathan D. Henshaw, Steven N. Longmore, Simon C. O. Glover, Mattia C. Sormani, Lucia Armillotta, Ashley T. Barnes, Ralf S. Klessen, Francisco Nogueras-Lara, Robin G. Tress, Jairo Armijos-Abendaño, Laura Colzi, Christoph Federrath, Pablo García, Adam Ginsburg, Christian Henkel, Sergio Martín, Denise Riquelme, Víctor M. Rivilla

Abstract: The Central Molecular Zone (CMZ; the central ~ 500 pc of the Galaxy) is a kinematically unusual environment relative to the Galactic disc, with high velocity dispersions and a steep size-linewidth relation of the molecular clouds. In addition, the CMZ region has a significantly lower star formation rate (SFR) than expected by its large amount of dense gas. An important factor in explaining the low… ▽ More The Central Molecular Zone (CMZ; the central ~ 500 pc of the Galaxy) is a kinematically unusual environment relative to the Galactic disc, with high velocity dispersions and a steep size-linewidth relation of the molecular clouds. In addition, the CMZ region has a significantly lower star formation rate (SFR) than expected by its large amount of dense gas. An important factor in explaining the low SFR is the turbulent state of the star-forming gas, which seems to be dominated by rotational modes. However, the turbulence driving mechanism remains unclear. In this work, we investigate how the Galactic gravitational potential affects the turbulence in CMZ clouds. We focus on the CMZ cloud G0.253+0.016 (`the Brick'), which is very quiescent and unlikely to be kinematically dominated by stellar feedback. We demonstrate that several kinematic properties of the Brick arise naturally in a cloud-scale hydrodynamics simulation that takes into account the Galactic gravitational potential. These properties include the line-of-sight velocity distribution, the steepened size-linewidth relation, and the predominantly solenoidal nature of the turbulence. Within the simulation, these properties result from the Galactic shear in combination with the cloud's gravitational collapse. This is a strong indication that the Galactic gravitational potential plays a crucial role in shaping the CMZ gas kinematics, and is a major contributor to suppressing the SFR by inducing predominantly solenoidal turbulent modes. △ Less

Submitted 1 August, 2023; v1 submitted 21 April, 2023; originally announced April 2023.

Comments: 7 pages, 8 figures; accepted to MNRAS (July 24th 2023)

Spatiotemporal dynamics of transonic shock-wave/turbulent-boundary-layer interactions in an overexpanded planar nozzle

Authors: Justin Kin Jun Hew, Emanuele Martelli, Mahdi Davoodianidalik, Rod W. Boswell, Christoph Federrath, Matthew Shadwell

Abstract: We perform a combined numerical and experimental study to investigate the transonic shock-wave/turbulent-boundary-layer interactions (STBLI) in a shock-induced separated subscale planar nozzle with fully-expanded Mach number,$M_j = 1.05$ and jet Reynolds number $Re \sim 10^5$. The nozzle configuration is tested via time-resolved schlieren visualisation. While numerous studies have been conducted o… ▽ More We perform a combined numerical and experimental study to investigate the transonic shock-wave/turbulent-boundary-layer interactions (STBLI) in a shock-induced separated subscale planar nozzle with fully-expanded Mach number,$M_j = 1.05$ and jet Reynolds number $Re \sim 10^5$. The nozzle configuration is tested via time-resolved schlieren visualisation. While numerous studies have been conducted on the high Reynolds number separated flowfields, little is known on the weak shock wave unsteadiness present in low nozzle pressure ratio (NPR) transonic nozzles. Therefore, numerical simulations are carried out with high resolution three-dimensional delayed detached eddy simulation (DDES), to study the spatiotemporal dynamics of wall pressure signals and unsteady shock interactions. The transient statistics considered include spectral Fourier and wavelet-based analysis and dynamic mode decomposition (DMD). The spectral analyses reveal energetic low frequency modes corresponding to the staging behaviour of shock unsteadiness, and high frequencies linked to the characteristics of the Kelvin-Helmholtz instabilities in the downstream turbulent mixing layer. The mechanisms for the low frequency unsteadiness is educed through modal decomposition and spectral analysis, wherein it is found that the downstream perturbations within the separation bubble play a major role in not only closing the aeroacoustic feedback loop, but allowing the continual evolution and sustainment of low frequency unsteadiness. An analysis via the vortex sheet method is also carried out to characterise the screech production, by assuming an upstream propagating guided jet mode △ Less

Submitted 13 April, 2023; originally announced April 2023.

Comments: 19 pages, 17 figures, submitted to PoF

Can we observe the ion-neutral drift velocity in prestellar cores?

Authors: Aris Tritsis, Shantanu Basu, Christoph Federrath

Abstract: Given the low ionization fraction of molecular clouds, ambipolar diffusion is thought to be an integral process in star formation. However, chemical and radiative-transfer effects, observational challenges, and the fact that the ion-neutral drift velocity is inherently very small render a definite detection of ambipolar diffusion extremely non-trivial. Here, we study the ion-neutral drift velocity… ▽ More Given the low ionization fraction of molecular clouds, ambipolar diffusion is thought to be an integral process in star formation. However, chemical and radiative-transfer effects, observational challenges, and the fact that the ion-neutral drift velocity is inherently very small render a definite detection of ambipolar diffusion extremely non-trivial. Here, we study the ion-neutral drift velocity in a suite of chemodynamical, non-ideal magnetohydrodynamic (MHD), two-dimensional axisymmetric simulations of prestellar cores where we alter the temperature, cosmic-ray ionization rate, visual extinction, mass-to-flux ratio, and chemical evolution. Subsequently, we perform a number of non-local thermodynamic equilibrium (non-LTE) radiative-transfer calculations considering various idealized and non-idealized scenarios in order to assess which factor (chemistry, radiative transfer and/or observational difficulties) is the most challenging to overcome in our efforts to detect the ion-neutral drift velocity. We find that temperature has a significant effect in the amplitude of the drift velocity with the coldest modelled cores (T = 6 K) exhibiting drift velocities comparable to the sound speed. Against expectations, we find that in idealized scenarios (where two species are perfectly chemically co-evolving) the drift velocity ``survives" radiative-transfer effects and can in principle be observed. However, we find that observational challenges and chemical effects can significantly hinder our view of the ion-neutral drift velocity. Finally, we propose that $\rm{HCN}$ and $\rm{HCNH^+}$, being chemically co-evolving, could be used in future observational studies aiming to measure the ion-neutral drift velocity. △ Less

Submitted 16 March, 2023; originally announced March 2023.

Comments: 14 pages, 11 figures. Accepted for publication in MNRAS

Compressible turbulence in the interstellar medium: New insights from a high-resolution supersonic turbulence simulation

Authors: Renaud Ferrand, Sébastien Galtier, Fouad Sahraoui, Christoph Federrath

Abstract: The role of supersonic turbulence in structuring the interstellar medium (ISM) remains an unsettled question. Here, this problem is investigated using a newexact law of compressible isothermal hydrodynamic turbulence, which involves two-point correlations in physical space. The new law is shown to have a compact expression that contains a single flux term reminiscent of the incompressible case and… ▽ More The role of supersonic turbulence in structuring the interstellar medium (ISM) remains an unsettled question. Here, this problem is investigated using a newexact law of compressible isothermal hydrodynamic turbulence, which involves two-point correlations in physical space. The new law is shown to have a compact expression that contains a single flux term reminiscent of the incompressible case and a source term with a simple expression whose sign is given by the divergence of the velocity. The law is then used to investigate the properties of such a turbulence at integral Mach number $4$ produced by a massive numerical simulation with a grid resolution of $10,048^3$ points. The flux (resp. source) term was found to have positive (resp. negative) contribution to the total energy cascade rate, which is interpreted as a direct cascade amplified by compression, while their sum is constant in the inertial range. Using a local (in space) analysis it is shown that the source is mainly driven by filamentary structures in which the flux is negligible. Taking positive defined correlations reveals the existence of different turbulent regimes separated by the sonic scale, which determines the scale over which the non-negligible source modifies the scaling of the flux. Our study provides new insight into the dynamics and structures of supersonic interstellar turbulence. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Comments: 8 pages, 6 Figures

Journal ref: The Astrophysical Journal, Volume 904, Number 2, Page 160, 2020

Multiphase condensation in cluster halos: interplay of cooling, buoyancy and mixing

Authors: Rajsekhar Mohapatra, Prateek Sharma, Christoph Federrath, Eliot Quataert

Abstract: Gas in the central regions of cool-core clusters and other massive halos has a short cooling time ($\lesssim1~\mathrm{Gyr}$). Theoretical models predict that this gas is susceptible to multiphase condensation, in which cold gas is expected to condense out of the hot phase if the ratio of the thermal instability growth time scale ($t_{\mathrm{ti}}$) to the free-fall time ($t_{\mathrm{ff}}$) is… ▽ More Gas in the central regions of cool-core clusters and other massive halos has a short cooling time ($\lesssim1~\mathrm{Gyr}$). Theoretical models predict that this gas is susceptible to multiphase condensation, in which cold gas is expected to condense out of the hot phase if the ratio of the thermal instability growth time scale ($t_{\mathrm{ti}}$) to the free-fall time ($t_{\mathrm{ff}}$) is $t_{\mathrm{ti}}/t_{\mathrm{ff}}\lesssim10$. The turbulent mixing time $t_{\mathrm{mix}}$ is another important time scale: if $t_{\mathrm{mix}}$ is short enough, the fluctuations are mixed before they can cool. In this study, we perform high-resolution ($512^2\times768$--$1024^2\times1536$ resolution elements) hydrodynamic simulations of turbulence in a stratified medium, including radiative cooling of the gas. We explore the parameter space of $t_{\mathrm{ti}}/t_{\mathrm{ff}}$ and $t_{\mathrm{ti}}/t_{\mathrm{mix}}$ relevant to galaxy and cluster halos. We also study the effect of the steepness of the entropy profile, the strength of turbulent forcing and the nature of turbulent forcing (natural mixture vs. compressive modes) on multiphase gas condensation. We find that larger values of $t_{\mathrm{ti}}/t_{\mathrm{ff}}$ or $t_{\mathrm{ti}}/t_{\mathrm{mix}}$ generally imply stability against multiphase gas condensation, whereas larger density fluctuations (e.g., due to compressible turbulence) promote multiphase gas condensation. We propose a new criterion $\min(t_{\mathrm{ti}}/\min(t_{\mathrm{mix}},t_\mathrm{ff}))\lesssim c_2\times\exp(c_1σ_s)$ for when the halo becomes multiphase, where $σ_s$ denotes the amplitude of logarithmic density fluctuations and $c_1\simeq6$, $c_2\simeq1.8$ from an empirical fit to our results. △ Less

Submitted 29 August, 2023; v1 submitted 18 February, 2023; originally announced February 2023.

Comments: MNRAS-accepted version, 16+2 (appendix) pages, 11 figures, simulation movies available at this playlist on youtube: https://youtube.com/playlist?list=PLuaNgQ1v_KMZlkKXdB7hcaQ7-hb0hmY7G

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

Authors: Sabrina M. Appel, Blakesley Burkhart, Vadim A. Semenov, Christoph Federrath, Anna L. Rosen, Jonathan C. Tan

Abstract: 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 cycl… ▽ More 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. △ Less

Submitted 18 January, 2023; originally announced January 2023.

Comments: 22 pages, 13 figures. Submitted to ApJ. Comments welcome

Lagrangian statistics of a shock-driven turbulent dynamo in decaying turbulence

Authors: Justin Kin Jun Hew, Christoph Federrath

Abstract: Small-scale fluctuating magnetic fields of order $n$G are observed in supernova shocks and galaxy clusters, where its amplification is likely caused by the Biermann battery mechanism. However, these fields cannot be amplified further without the turbulent dynamo, which generates magnetic energy through the stretch-twist-fold (STF) mechanism. Thus, we present here novel three-dimensional magnetohyd… ▽ More Small-scale fluctuating magnetic fields of order $n$G are observed in supernova shocks and galaxy clusters, where its amplification is likely caused by the Biermann battery mechanism. However, these fields cannot be amplified further without the turbulent dynamo, which generates magnetic energy through the stretch-twist-fold (STF) mechanism. Thus, we present here novel three-dimensional magnetohydrodynamic (MHD) simulations of a laser-driven shock propagating into a stratified, multiphase medium, to investigate the post-shock turbulent magnetic field amplification via the turbulent dynamo. The configuration used here is currently being tested in the shock tunnel at the National Ignition Facility (NIF). In order to probe the statistical properties of the post-shock turbulent region, we use $384 \times 512 \times 384$ tracers to track its evolution through the Lagrangian framework, thus providing a high-fidelity analysis of the shocked medium. Our simulations indicate that the growth of the magnetic field, which accompanies the near-Saffman kinetic energy decay ($E_{\textrm{kin}} \propto t^{-1.15})$ without turbulence driving, exhibits slightly different characteristics as compared to periodic box simulations. Seemingly no distinct phases exist in its evolution, because the shock passage and time to observe the magnetic field amplification during the turbulence decay are very short ($\sim\!0.3$ of a turbulent turnover time). Yet, the growth rate is still consistent with those expected for compressive (curl-free) turbulence driving in subsonic, compressible turbulence. Phenomenological understanding of the dynamics of the magnetic and velocity fields are also elucidated via Lagrangian frequency spectra, which are consistent with the expected inertial range scalings in the Eulerian-Lagrangian bridge. △ Less

Submitted 16 February, 2023; v1 submitted 15 January, 2023; originally announced January 2023.

Comments: 15 pages, 19 figures. Accepted in MNRAS. Comments are welcome

Analytical and Numerical Studies of the Non-uniformity induced Type II Asymmetric Cap Shock Mach Reflection in Over-expanded Supersonic Jets

Authors: Justin Kin Jun Hew, Rod W. Boswell, Christoph Federrath, Rajesh Gopalapillai, Vinoth Paramanantham

Abstract: A combined analytical and numerical study is conducted to investigate the asymmetric cap-shock non-uniform Mach Reflection (csMR) phenomenon outside of an over-expanded supersonic jet. Prior analytical works have only considered the wedge-induced steady symmetric and asymmetric Mach reflection (MR) configurations, as well as symmetric MR in an open jet. However, there is another structure occurrin… ▽ More A combined analytical and numerical study is conducted to investigate the asymmetric cap-shock non-uniform Mach Reflection (csMR) phenomenon outside of an over-expanded supersonic jet. Prior analytical works have only considered the wedge-induced steady symmetric and asymmetric Mach reflection (MR) configurations, as well as symmetric MR in an open jet. However, there is another structure occurring in nozzle flow fields known as a non-uniformity induced cap-shock pattern. We derive a new analytical model to predict the wave structure of the asymmetric csMR in the absence of internal shocks by extending on a prior symmetric Mach reflection model. Different from the wedge flow case, flow non-uniformity is incorporated by assuming different upstream Mach numbers in both upper and lower domains, where the asymmetry is predicted through averaged flowfields and slipstream inclination angles. The numerical approach utilises an Euler solver for comparisons to the developed theory. It is found that the model adequately predicts the shock structure obtained from numerical simulations, and can be utilised for various sets of parameters to capture the overall direct Mach reflection o[DiMR+DiMR] configuration. The von Neumann criterion is also well captured by the new analytic model, along with the Mach stem profile and shock curvatures. Based on the analytical and numerical observations, a hypothesis is also made regarding the stability of MR structures within an over-expanded jet. △ Less

Submitted 13 February, 2023; v1 submitted 10 October, 2022; originally announced October 2022.

Comments: 33 pages, 19 figures

Outflows Driven by Direct and Reprocessed Radiation Pressure in Massive Star Clusters

Authors: Shyam H. Menon, Christoph Federrath, Mark R. Krumholz

Abstract: We use three-dimensional radiation hydrodynamic (RHD) simulations to study the formation of massive star clusters under the combined effects of direct ultraviolet (UV) and dust-reprocessed infrared (IR) radiation pressure. We explore a broad range of mass surface density $Σ\sim 10^2$-$10^5 \, \mathrm{M}_{\odot} \, \mathrm{pc}^{-2}$, spanning values typical of weakly star-forming galaxies to extrem… ▽ More We use three-dimensional radiation hydrodynamic (RHD) simulations to study the formation of massive star clusters under the combined effects of direct ultraviolet (UV) and dust-reprocessed infrared (IR) radiation pressure. We explore a broad range of mass surface density $Σ\sim 10^2$-$10^5 \, \mathrm{M}_{\odot} \, \mathrm{pc}^{-2}$, spanning values typical of weakly star-forming galaxies to extreme systems such as clouds forming super-star clusters, where radiation pressure is expected to be the dominant feedback mechanism. We find that star formation can only be regulated by radiation pressure for $Σ\lesssim 10^3 \, \mathrm{M}_{\odot} \, \mathrm{pc}^{-2}$, but that clouds with $Σ\lesssim 10^5 \, \mathrm{M}_{\odot} \, \mathrm{pc}^{-2}$ become super-Eddington once high star formation efficiencies ($\sim 80 \%$) are reached, and therefore launch the remaining gas in a steady outflow. These outflows achieve mass-weighted radial velocities of $\sim 15$ - $30 \,\mathrm{km} \, \mathrm{s}^{-1}$, which is $\sim 0.5$ - $2.0$ times the cloud escape speed. This suggests that radiation pressure is a strong candidate to explain recently observed molecular outflows found in young super-star clusters in nearby starburst galaxies. We quantify the relative importance of UV and IR radiation pressure in different regimes, and deduce that both are equally important for $Σ\sim 10^3 \, \mathrm{M}_{\odot} \, \mathrm{pc}^{-2}$, whereas clouds with higher (lower) density are increasingly dominated by the IR (UV) component. Comparison with control runs without either the UV or IR bands suggests that the outflows are primarily driven by the impulse provided by the UV component, while IR radiation has the effect of rendering a larger fraction of gas super-Eddington, and thereby increasing the outflow mass flux by a factor of $\sim 2$. △ Less

Submitted 20 March, 2023; v1 submitted 6 October, 2022; originally announced October 2022.

Comments: 15 pages, 11 figures. MNRAS accepted. v2: Minor changes in text made to address referee comments

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.

Comments: 16 pages. 11 figures. Accepted to MNRAS. Part I focuses on the journey towards saturation. Part II focuses on the decay into saturation

The role of the turbulence driving mode for the Initial Mass Function

Authors: Sajay Sunny Mathew, Christoph Federrath, Amit Seta

Abstract: Turbulence is a critical ingredient for star formation, yet its role for the initial mass function (IMF) is not fully understood. Here we perform magnetohydrodynamical (MHD) simulations of star cluster formation including gravity, turbulence, magnetic fields, stellar heating and outflow feedback to study the influence of the mode of turbulence driving on IMF. We find that simulations that employ p… ▽ More Turbulence is a critical ingredient for star formation, yet its role for the initial mass function (IMF) is not fully understood. Here we perform magnetohydrodynamical (MHD) simulations of star cluster formation including gravity, turbulence, magnetic fields, stellar heating and outflow feedback to study the influence of the mode of turbulence driving on IMF. We find that simulations that employ purely compressive turbulence driving (COMP) produce a higher fraction of low-mass stars as compared to simulations that use purely solenoidal driving (SOL). The characteristic (median) mass of the sink particle (protostellar) distribution for COMP is shifted to lower masses by a factor of ~ 1.5 compared to SOL. Our simulation IMFs capture the important features of the observed IMF form. We find that turbulence-regulated theories of the IMF match our simulation IMFs reasonably well in the high-mass and low-mass range, but underestimate the number of very low-mass stars, which form towards the later stages of our simulations and stop accreting due to dynamical interactions. Our simulations show that for both COMP and SOL, the multiplicity fraction is an increasing function of the primary mass, although the multiplicity fraction in COMP is higher than that of SOL for any primary mass range. We find that binary mass ratio distribution is independent of the turbulence driving mode. The average specific angular momentum of the sink particles in SOL is a factor of 2 higher than that for COMP. Overall, we conclude that the turbulence driving mode plays a significant role in shaping the IMF. △ Less

Submitted 25 November, 2022; v1 submitted 18 August, 2022; originally announced August 2022.

Comments: 26 pages, 16 figures, 4 tables, accepted for publication in MNRAS, minor changes, part of methodology is same as in Mathew & Federrath 2021. arXiv admin note: text overlap with arXiv:2106.06521

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.

Comments: Astronomy and Astrophysics in presss

Journal ref: A&A 666, A165 (2022)

Infrared Radiation Feedback Does Not Regulate Star Cluster Formation

Authors: Shyam H. Menon, Christoph Federrath, Mark R. Krumholz

Abstract: We present 3D radiation-hydrodynamical (RHD) simulations of star cluster formation and evolution in massive, self-gravitating clouds, whose dust columns are optically thick to infrared (IR) photons. We use \texttt{VETTAM} -- a recently developed, novel RHD algorithm, which uses the Variable Eddington Tensor (VET) closure -- to model the IR radiation transport through the cloud. We also use realist… ▽ More We present 3D radiation-hydrodynamical (RHD) simulations of star cluster formation and evolution in massive, self-gravitating clouds, whose dust columns are optically thick to infrared (IR) photons. We use \texttt{VETTAM} -- a recently developed, novel RHD algorithm, which uses the Variable Eddington Tensor (VET) closure -- to model the IR radiation transport through the cloud. We also use realistic temperature ($T$) dependent IR opacities ($κ$) in our simulations, improving upon earlier works in this area, which used either constant IR opacities or simplified power laws ($κ\propto T^2$). We investigate the impact of the radiation pressure of these IR photons on the star formation efficiency (SFE) of the cloud, and its potential to drive dusty winds. We find that IR radiation pressure is unable to regulate star formation or prevent accretion onto the star clusters, even for very high gas surface densities ($Σ> 10^5 M_{\odot} \, \mathrm{pc}^{-2}$), contrary to recent semi-analytic predictions and simulation results using simplified treatments of the dust opacity. We find that the commonly adopted simplifications of $κ\propto T^2$ or constant $κ$ for the IR dust opacities leads to this discrepancy, as those approximations overestimate the radiation force. By contrast, with realistic opacities that take into account the micro-physics of the dust, we find that the impact of IR radiation pressure on star formation is very mild, even at significantly high dust-to-gas ratios ($\sim 3$ times solar), suggesting that it is unlikely to be an important feedback mechanism in controlling star formation in the ISM. △ Less

Submitted 22 September, 2022; v1 submitted 28 June, 2022; originally announced June 2022.

Comments: 28 pages, 19 figures. MNRAS accepted. v2: Minor changes made to address referee comments

Rotation measure structure functions with higher-order stencils as a probe of small-scale magnetic fluctuations and its application to the Small and Large Magellanic Clouds

Authors: Amit Seta, Christoph Federrath, Jack D. Livingston, N. M. McClure-Griffiths

Abstract: Magnetic fields and turbulence are important components of the interstellar medium (ISM) of star-forming galaxies. It is challenging to measure the properties of the small-scale ISM magnetic fields (magnetic fields at scales smaller than the turbulence driving scale). Using numerical simulations, we demonstrate how the second-order rotation measure (RM, which depends on thermal electron density,… ▽ More Magnetic fields and turbulence are important components of the interstellar medium (ISM) of star-forming galaxies. It is challenging to measure the properties of the small-scale ISM magnetic fields (magnetic fields at scales smaller than the turbulence driving scale). Using numerical simulations, we demonstrate how the second-order rotation measure (RM, which depends on thermal electron density, $n_{\rm e}$, and magnetic field, $b$) structure function can probe the properties of small-scale $b$. We then apply our results to observations of the Small and Large Magellanic Clouds (SMC and LMC). First, using Gaussian random $b$, we show that the characteristic scale where the RM structure function flattens is approximately equal to the correlation length of $b$. We also show that computing the RM structure function with a higher-order stencil (more than the commonly-used two-point stencil) is necessary to accurately estimate the slope of the structure function. Then, using Gaussian random $b$ and lognormal $n_{\rm e}$ with known power spectra, we derive an empirical relationship between the slope of the power spectrum of $b$, $n_{\rm e}$, and RM. We apply these results to the SMC and LMC and estimate the following properties of small-scale $b$: correlation length ($160~\pm 21~{\rm pc}$ for the SMC and $87~\pm~17~{\rm pc}$ for the LMC), strength ($14~\pm 2~μ{\rm G}$ for the SMC and $15~\pm 3~μ{\rm G}$ for the LMC), and slope of the magnetic power spectrum ($-1.3~\pm~0.4$ for the SMC and $-1.6~\pm~0.1$ for the LMC). We also find that $n_{\rm e}$ is practically constant over the estimated $b$ correlation scales. △ Less

Submitted 13 October, 2022; v1 submitted 28 June, 2022; originally announced June 2022.

Comments: 27 pages (including 6 appendices), 29 figures (16 in the main text and 13 in the appendices), and 4 tables (2 in the main text and 2 in the appendices); MNRAS accepted

Environmental variation of the low-mass IMF

Authors: Tabassum S Tanvir, Mark R Krumholz, Christoph Federrath

Abstract: We use a series of magnetohydrodynamic simulations including both radiative and protostellar outflow feedback to study the environmental variation of the initial mass function. The simulations represent a carefully-controlled experiment whereby we keep all dimensionless parameters of the flow constant except for those related to feedback. We show that radiation feedback suppresses the formation of… ▽ More We use a series of magnetohydrodynamic simulations including both radiative and protostellar outflow feedback to study the environmental variation of the initial mass function. The simulations represent a carefully-controlled experiment whereby we keep all dimensionless parameters of the flow constant except for those related to feedback. We show that radiation feedback suppresses the formation of lower mass objects more effectively as the surface density increases, but this only partially compensates for the decreasing Jeans mass in denser environments. Similarly, we find that protostellar outflows are more effective at suppressing the formation of massive stars in higher surface density environments. The combined effect of these two trends is towards an IMF with a lower characteristic mass and a narrower overall mass range in high surface density environments. We discuss the implications of these findings for the interpretation of observational evidence of IMF variation in early-type galaxies. △ Less

Submitted 14 September, 2022; v1 submitted 10 June, 2022; originally announced June 2022.

Comments: 15 pages, 13 figures, accepted by MNRAS

Multiphase turbulence in galactic halos: effect of the driving

Authors: Rajsekhar Mohapatra, Christoph Federrath, Prateek Sharma

Abstract: 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 effe… ▽ More 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. △ Less

Submitted 7 June, 2022; originally announced June 2022.

Comments: 21 pages, 12 figures, accepted for publication in MNRAS. Simulation movies are available at this link: https://www.youtube.com/watch?v=qBsJti2R0HY&list=PLuaNgQ1v_KMaovGyz-7jbFha-HvoXpG2P

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

Authors: Seth Davidovits, Christoph Federrath, Romain Teyssier, Kumar S. Raman, David C. Collins, Sabrina R. Nagel

Abstract: 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… ▽ More 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. △ Less

Submitted 31 May, 2022; originally announced May 2022.

Comments: 15 pages, 10 figures

The driving mode of shock-driven turbulence

Authors: Saee Dhawalikar, Christoph Federrath, Seth Davidovits, Romain Teyssier, Sabrina R. Nagel, Bruce A. Remington, David C. Collins

Abstract: 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, diverg… ▽ More 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. △ Less

Submitted 28 May, 2022; originally announced May 2022.

Comments: 18 pages, 15 figures

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.

Comments: Accepted in MNRAS

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.

Comments: 15 pages, 10 figures; accepted for publication in MNRAS, 2022 April 4

MSC Class: 85-10

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.

Comments: Invited submission accepted for publication in Frontiers in Astronomy and Space Sciences: Multi-scale Magnetic Field Measurements in the Multi-Phase Interstellar Medium; 46 pages; 12 Figures, 2 tables

A radiation hydrodynamics scheme on adaptive meshes using the Variable Eddington Tensor (VET) closure

Authors: Shyam H. Menon, Christoph Federrath, Mark R. Krumholz, Rolf Kuiper, Benjamin D. Wibking, Manuel Jung

Abstract: We present a new algorithm to solve the equations of radiation hydrodynamics (RHD) in a frequency-integrated, two-moment formulation. Novel features of the algorithm include i) the adoption of a non-local Variable Eddington Tensor (VET) closure for the radiation moment equations, computed with a ray-tracing method, ii) support for adaptive mesh refinement (AMR), iii) use of a time-implicit Godunov… ▽ More We present a new algorithm to solve the equations of radiation hydrodynamics (RHD) in a frequency-integrated, two-moment formulation. Novel features of the algorithm include i) the adoption of a non-local Variable Eddington Tensor (VET) closure for the radiation moment equations, computed with a ray-tracing method, ii) support for adaptive mesh refinement (AMR), iii) use of a time-implicit Godunov method for the hyperbolic transport of radiation, and iv) a fixed-point Picard iteration scheme to accurately handle the stiff nonlinear gas-radiation energy exchange. Tests demonstrate that our scheme works correctly, yields accurate rates of energy and momentum transfer between gas and radiation, and obtains the correct radiation field distribution even in situations where more commonly used -- but less accurate -- closure relations like the Flux-limited Diffusion and Moment-1 approximations fail. Our scheme presents an important step towards performing RHD simulations with increasing spatial and directional accuracy, effectively improving their predictive capabilities. △ Less

Submitted 22 March, 2022; originally announced March 2022.

Comments: 7 pages, 2 figures, to appear in the proceedings of IAUS362: The predictive power of computational astrophysics as discovery tool (D. Bisikalo, T. Hanawa, C. Boily, J. Stone, eds.). arXiv admin note: substantial text overlap with arXiv:2202.08778

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

Authors: Yue Hu, Christoph Federrath, Siyao Xu, Sajay Sunny Mathew

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

Submitted 21 April, 2022; v1 submitted 2 March, 2022; originally announced March 2022.

Comments: 11 pages, 9 figures, accepted for publication in MNRAS

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.

Comments: Accepted in MNRAS. 20 pages. 16 figures. 1 table

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

Authors: Shyam H. Menon, Christoph Federrath, Mark R. Krumholz, Rolf Kuiper, Benjamin D. Wibking, Manuel Jung

Abstract: 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.… ▽ More 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. △ Less

Submitted 17 February, 2022; originally announced February 2022.

Comments: 23 pages, 15 figures. Accepted for publication in MNRAS

Turbulent dynamo in the two-phase interstellar medium

Authors: Amit Seta, Christoph Federrath

Abstract: 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… ▽ More 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. △ Less

Submitted 17 May, 2022; v1 submitted 16 February, 2022; originally announced February 2022.

Comments: 18 pages (including 3 appendices), 17 figures (12 in the main text and 5 in the appendices), and 3 tables; accepted in MNRAS

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.

Comments: 19 pages, 11 figures, accepted for publication in ApJ

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

Authors: Aris Tritsis, Christoph Federrath, Karen Willacy, Konstantinos Tassis

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

Submitted 21 December, 2021; originally announced December 2021.

Comments: 18 pages, 14 figures. Accepted for publication in MNRAS