Search

Advancing our understanding of astrophysical turbulence is bottlenecked by the limited resolution of numerical simulations that may not fully sample scales in the inertial range. Machine-learning (ML) techniques have demonstrated promise in upscaling resolution in both image analysis and numerical simulations (i.e., superresolution). Here we employ and further develop a physics-constrained convolutional neural network ML model called "MeshFreeFlowNet" (MFFN) for superresolution studies of turbulent systems. The model is trained on both the simulation images and the evaluated partial differential equations (PDEs), making it sensitive to the underlying physics of a particular fluid system. We develop a framework for 2D turbulent Rayleigh–Bénard convection generated with the Dedalus code by modifying the MFFN architecture to include the full set of simulation PDEs and the boundary conditions. Our training set includes fully developed turbulence sampling Rayleigh numbers (Ra) of Ra = 106–1010. We evaluate the success of the learned simulations by comparing the power spectra of the direct Dedalus simulation to the predicted model output and compare both ground-truth and predicted power spectral inertial range scalings to theoretical predictions. We find that the updated network performs well at all Ra studied here in recovering large-scale information, including the inertial range slopes. The superresolution prediction is overly dissipative at smaller scales than that of the inertial range in all cases, but the smaller scales are better recovered in more turbulent than laminar regimes. This is likely because more turbulent systems have a rich variety of structures at many length scales compared to laminar flows. [ABSTRACT FROM AUTHOR]

We employ the Feedback In Realistic Environments (FIRE-2) physics model to study how the properties of giant molecular clouds (GMCs) evolve during galaxy mergers. We conduct a pixel-by-pixel analysis of molecular gas properties in both the simulated control galaxies and galaxy major mergers. The simulated GMC pixels in the control galaxies follow a similar trend in a diagram of velocity dispersion (σ v ) versus gas surface density (Σmol) to the one observed in local spiral galaxies in the Physics at High Angular resolution in Nearby GalaxieS (PHANGS) survey. For GMC pixels in simulated mergers, we see a significant increase of a factor of 5–10 in both Σmol and σ v , which puts these pixels above the trend of PHANGS galaxies in the σ v versus Σmol diagram. This deviation may indicate that GMCs in the simulated mergers are much less gravitationally bound compared with simulated control galaxies with virial parameters (α vir) reaching 10–100. Furthermore, we find that the increase in α vir happens at the same time as the increase in global star formation rate, which suggests that stellar feedback is responsible for dispersing the gas. We also find that the gas depletion time is significantly lower for high- α vir GMCs during a starburst event. This is in contrast to the simple physical picture that low- α vir GMCs are easier to collapse and form stars on shorter depletion times. This might suggest that some other physical mechanisms besides self-gravity are helping the GMCs in starbursting mergers collapse and form stars. [ABSTRACT FROM AUTHOR]

The driving of turbulence in galaxies is deeply connected with the physics of feedback, star formation, outflows, accretion, and radial transport in disks. The velocity dispersion of gas in galaxies therefore offers a promising observational window into these processes. However, the relative importance of each of these mechanisms remains controversial. In this work we revisit the possibility that turbulence on galactic scales is driven by the direct impact of accreting gaseous material on the disk. We measure this effect in a disk-like star-forming galaxy in IllustrisTNG, using the high-resolution cosmological magnetohydrodynamical simulation TNG50. We employ Lagrangian tracer particles with a high time cadence of only a few million years to identify accretion and other events. The energies of particles are measured by stacking the events in bins of time around the event. The average effect of each event is measured by fitting explicit models for the kinetic and turbulent energies as a function of time. These measurements are corroborated by cross-correlating the turbulent energy with other time series and searching for signals of causality, i.e., asymmetries across zero time lag. We find that accretion contributes to the large-scale turbulent kinetic energy even if it does not dominate in this ∼5 × 109 M ⊙ stellar mass galaxy. Extrapolating this finding to a range of galaxy masses, we find that there are regimes where energy from direct accretion may dominate the turbulent energy budget, particularly in disk outskirts, galaxies less massive than the Milky Way, and at redshift ∼2. [ABSTRACT FROM AUTHOR]

A supersonic relative velocity between dark matter (DM) and baryons (the stream velocity) at the time of recombination induces the formation of low-mass objects with anomalous properties in the early universe. We widen the scope of the "Supersonic Project" paper series to include objects we term Dark Matter + Gas Halos Offset by Streaming (DM GHOSts)—diffuse, DM-enriched structures formed because of a physical offset between the centers of mass of DM and baryonic overdensities. We present an updated numerical investigation of DM GHOSts and Supersonically Induced Gas Objects (SIGOs), including the effects of molecular cooling, in high-resolution hydrodynamic simulations using the AREPO code. Supplemented by an analytical understanding of their ellipsoidal gravitational potentials, we study the population-level properties of these objects, characterizing their morphology, spin, radial mass, and velocity distributions in comparison to classical structures in non-streaming regions. The stream velocity causes deviations from sphericity in both the gas and DM components and lends greater rotational support to the gas. Low-mass (≲105.5 M ⊙) objects in regions of streaming demonstrate core-like rotation and mass profiles. Anomalies in the rotation and morphology of DM GHOSts could represent an early universe analog to observed ultra-faint dwarf galaxies with variations in DM content and unusual rotation curves. [ABSTRACT FROM AUTHOR]

We present an analytic model for clustered supernovae (SNe) feedback in galaxy disks, incorporating the dynamical evolution of superbubbles formed from spatially overlapping SNe remnants. We propose two realistic outcomes for the evolution of superbubbles in galactic disks: (1) the expansion velocity of the shock front falls below the turbulent velocity dispersion of the interstellar medium in the galaxy disk, whereupon the superbubble stalls and fragments, depositing its momentum entirely within the galaxy disk; or (2) the superbubble grows in size to reach the gas scale height, breaking out of the galaxy disk and driving galactic outflows/fountains. In either case, we find that superbubble breakup/breakout almost always occurs before the last Type II SN (≲40 Myr) in the recently formed star cluster, assuming a standard high-end initial mass function slope, and scalings between stellar lifetimes and masses. The threshold between these two cases implies a break in the effective strength of feedback in driving turbulence within galaxies, and a resulting change in the scalings of, for example, star formation rates with gas surface density (the Kennicuttâ€"Schmidt relation) and the star formation efficiency in galaxy disks. [ABSTRACT FROM AUTHOR]

We present the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project. CAMELS is a suite of 4233 cosmological simulations of volume each: 2184 state-of-the-art (magneto)hydrodynamic simulations run with the AREPO and GIZMO codes, employing the same baryonic subgrid physics as the IllustrisTNG and SIMBA simulations, and 2049 N-body simulations. The goal of the CAMELS project is to provide theory predictions for different observables as a function of cosmology and astrophysics, and it is the largest suite of cosmological (magneto)hydrodynamic simulations designed to train machine-learning algorithms. CAMELS contains thousands of different cosmological and astrophysical models by way of varying Ωm, σ8, and four parameters controlling stellar and active galactic nucleus feedback, following the evolution of more than 100 billion particles and fluid elements over a combined volume of. We describe the simulations in detail and characterize the large range of conditions represented in terms of the matter power spectrum, cosmic star formation rate density, galaxy stellar mass function, halo baryon fractions, and several galaxy scaling relations. We show that the IllustrisTNG and SIMBA suites produce roughly similar distributions of galaxy properties over the full parameter space but significantly different halo baryon fractions and baryonic effects on the matter power spectrum. This emphasizes the need for marginalizing over baryonic effects to extract the maximum amount of information from cosmological surveys. We illustrate the unique potential of CAMELS using several machine-learning applications, including nonlinear interpolation, parameter estimation, symbolic regression, data generation with Generative Adversarial Networks, dimensionality reduction, and anomaly detection. [ABSTRACT FROM AUTHOR]

The complex interplay of magnetohydrodynamics, gravity, and supersonic turbulence in the interstellar medium (ISM) introduces a non-Gaussian structure that can complicate a comparison between theory and observation. In this paper, we show that the wavelet scattering transform (WST), in combination with linear discriminant analysis (LDA), is sensitive to non-Gaussian structure in 2D ISM dust maps. WST-LDA classifies magnetohydrodynamic (MHD) turbulence simulations with up to a 97% true positive rate in our testbed of 8 simulations with varying sonic and Alfvénic Mach numbers. We present a side-by-side comparison with two other methods for non-Gaussian characterization, the reduced wavelet scattering transform (RWST) and the three-point correlation function (3PCF). We also demonstrate the 3D-WST-LDA, and apply it to the classification of density fields in position–position–velocity (PPV) space, where density correlations can be studied using velocity coherence as a proxy. WST-LDA is robust to common observational artifacts, such as striping and missing data, while also being sensitive enough to extract the net magnetic field direction for sub-Alfvénic turbulent density fields. We include a brief analysis of the effect of point-spread functions and image pixelization on 2D-WST-LDA applied to density fields, which informs the future goal of applying WST-LDA to 2D or 3D all-sky dust maps to extract hydrodynamic parameters of interest. [ABSTRACT FROM AUTHOR]

Supersonically induced gas objects (SIGOs) primarily form in the early universe, outside of dark matter halos due to the presence of a relative stream velocity between baryons and dark matter. These structures may be the progenitors of globular clusters. Since SIGOs are made out of pristine gas, we investigate the effect of atomic cooling on their properties. We run a suite of simulations by using the moving-mesh code arepo , with and without baryon dark matter relative velocity and with and without the effects of atomic cooling. We show that SIGO's density, temperature, and prolateness are determined by gravitational interactions rather than cooling. The cold gas fraction in SIGOs is much higher than that of dark matter halos. Specifically, we show that SIGO's characteristically low temperature and extremely high gas density forges a nurturing ground for the earliest star formation sites. [ABSTRACT FROM AUTHOR]

Semianalytic models (SAMs) are a promising means of tracking the physical processes associated with galaxy formation, but many of their approximations have not been rigorously tested. As part of the Simulating Multiscale Astrophysics to Understand Galaxies project, we compare predictions from the FIRE-2 hydrodynamical "zoom-in" simulations to those from the Santa Cruz SAM run on the same halo merger trees, with an emphasis on the global mass flow cycle. Our study includes 13 halos spanning low-mass dwarfs (Mvir ∼ 1010M⊙ at z = 0), intermediate-mass dwarfs (Mvir ∼ 1011M⊙), and Milky Way–mass galaxies (Mvir ∼ 1012M⊙). The SAM and FIRE-2 predictions agree relatively well with each other in terms of stellar and interstellar medium mass but differ dramatically on circumgalactic medium mass (the SAM is lower than FIRE-2 by ∼3 orders of magnitude for dwarfs). Strikingly, the SAM predicts higher gas accretion rates for dwarfs compared to FIRE-2 by factors of ∼10–100, and this is compensated for with higher mass outflow rates in the SAM. We argue that the most severe model discrepancies are caused by the lack of preventative stellar feedback and the assumptions for halo gas cooling and recycling in the SAM. As a first step toward resolving these model tensions, we present a simple yet promising new preventative stellar feedback model in which the energy carried by supernova-driven winds is allowed to heat some fraction of gas outside of halos to at least the virial temperature such that accretion is suppressed. [ABSTRACT FROM AUTHOR]

Radiative shock waves in the Cygnus Loop and other supernova remnants show different morphologies in [O iii ] and Hα emission. We use HST spectra and narrowband images to study the development of turbulence in the cooling region behind a shock on the west limb of the Cygnus Loop. We refine our earlier estimates of shock parameters that were based upon ground-based spectra, including ram pressure, vorticity, and magnetic field strength. We apply several techniques, including Fourier power spectra and the Rolling Hough Transform, to quantify the shape of the rippled shock front as viewed in different emission lines. We assess the relative importance of thermal instabilities, the thin shell instability, upstream density variations, and upstream magnetic field variations in producing the observed structure. [ABSTRACT FROM AUTHOR]

We have obtained a contiguous set of long-slit spectra of a shock wave in the Cygnus Loop to investigate its structure, which is far from the morphology predicted by 1D models. Proper motions from Hubble Space Telescope images combined with the known distance to the Cygnus Loop provide an accurate shock speed. Earlier analyses of shock spectra estimated the shock speed, postshock density, temperature, and elemental abundances. In this paper we determine several more shock parameters: a more accurate shock speed, ram pressure, density, compression ratio, dust destruction efficiency, magnetic field strength, and vorticity in the cooling region. From the derived shock properties we estimate the emissivities of synchrotron emission in the radio and pion decay emission in the gamma-rays. Both are consistent with the observations if we assume simple adiabatic compression of ambient cosmic-rays as in the van der Laan mechanism. We also find that, although the morphology is far from that predicted by 1D models and the line ratios vary dramatically from point to point, the average spectrum is matched reasonably well by 1D shock models with the shock speed derived from the measured proper motion. A subsequent paper will analyze the development of turbulence in the cooling zone behind the shock. [ABSTRACT FROM AUTHOR]

We investigate turbulent properties of the non-star-forming, translucent molecular cloud MBM 16 by applying the statistical technique of a two-dimensional spatial power spectrum (SPS) on the neutral hydrogen (H I) observations obtained by the Galactic Arecibo L-Band Feed Array H I survey. The SPS, calculated over the range of spatial scales from 0.1 to 17 pc, is well represented with a single power-law function, with a slope ranging from –3.3 to –3.7 and being consistent over the velocity range of MBM 16 for a fixed velocity channel thickness. However, the slope varies significantly with the velocity slice thickness, suggesting that both velocity and density contribute to H I intensity fluctuations. By using this variation, we estimate the slope of three-dimensional density fluctuations in MBM 16 to be –3.7 ± 0.2. This is significantly steeper than what has been found for H I in the Milky Way plane, the Small Magellanic Cloud, or the Magellanic Bridge, suggesting that interstellar turbulence in MBM 16 is driven on scales >17 pc and that the lack of stellar feedback could be responsible for the steep power spectrum. [ABSTRACT FROM AUTHOR]

The Fourier power spectrum is one of the most widely used statistical tools to analyze the nature of magnetohydrodynamic (MHD) turbulence in the interstellar medium. Lazarian & Pogosyan predicted that the spectral slope should saturate to -3 for an optically thick medium and many observations exist in support of their prediction. However, there have not been any numerical studies to date for testing these results. We analyze the spatial power spectrum of MHD simulations with a wide range of sonic and Alfvénic Mach numbers, which include radiative transfer effects of the 13CO transition. We numerically confirm the predictions of Lazarian & Pogosyan that the spectral slope of line intensity maps of an optically thick medium saturates to -3. Furthermore, for very optically thin supersonic CO gas, where the density or CO abundance values are too low to excite emission in all but the densest shock compressed gas, we find that the spectral slope is shallower than expected from the column density. Finally,we find that mixed optically thin/thick COgas,which has average optical depths on the order of unity, shows mixed behavior: for super-Alfvénic turbulence, the integrated intensity power spectral slopes generally follow the same trend with sonic Mach number as the true column density power spectrum slopes. However, for sub-Alfvénic turbulence the spectral slopes are steeper with values near -3 which are similar to the very optically thick regime. [ABSTRACT FROM AUTHOR]

We study the effects of radiative transfer on the probability distribution functions (PDFs) of simulations of magnetohydrodynamic turbulence in the widely studied 13CO 2-1 transition. We find that the integrated intensity maps generally follow a log-normal distribution, with the cases that have τ ≈ 1 best matching the PDF of the column density. We fit a two-dimensional variance-sonic Mach number relationship to our logarithmic PDFs of the form σ2 ln(Σ/Σ0) = A × ln(1 + b2M2 s ) and find that, for parameter b = 1/3, parameter A depends on the radiative transfer environment. We also explore the variance, skewness, and kurtosis of the linear PDFs finding that higher moments reflect both higher sonic Mach number and lower optical depth. Finally, we apply the Tsallis incremental PDF function and find that the fit parameters depend on both Mach numbers, but also are sensitive to the radiative transfer parameter space, with the τ ≈ 1 case best fitting the incremental PDF of the true column density. We conclude that, for PDFs of low optical depth cases, part of the gas is always subthermally excited so that the spread of the line intensities exceeds the spread of the underlying column densities and hence the PDFs do not reflect the true column density. Similarly, PDFs of optically thick cases are dominated by the velocity dispersion and therefore do not represent the true column density PDF. Thus, in the case of molecules like carbon monoxide, the dynamic range of intensities, structures observed, and, consequently, the observable PDFs are less determined by turbulence and more often determined by radiative transfer effects. [ABSTRACT FROM AUTHOR]

Magnetohydrodynamic turbulence is able to create hierarchical structures in the interstellar medium (ISM) that are correlated on a wide range of scales via the energy cascade. We use hierarchical tree diagrams known as dendrograms to characterize structures in synthetic position-position-velocity (PPV) emission cubes of isothermal magnetohydrodynamic turbulence. We show that the structures and degree of hierarchy observed in PPV space are related to the presence of self-gravity and the global sonic and Alfvénic Mach numbers. Simulations with higher Alfvénic Mach number, self-gravity and supersonic flows display enhanced hierarchical structure. We observe a strong dependency on the sonic and Alfvénic Mach numbers and self-gravity when we apply the statistical moments (i.e., mean, variance, skewness, kurtosis) to the leaf and node distribution of the dendrogram. Simulations with self-gravity, larger magnetic field and higher sonic Mach number have dendrogram distributions with higher statistical moments. Application of the dendrogram to three-dimensional density cubes, also known as position-position-position (PPP) cubes, reveals that the dominant emission contours in PPP and PPV are related for supersonic gas but not for subsonic. We also explore the effects of smoothing, thermal broadening, and velocity resolution on the dendrograms in order to make our study more applicable to observational data. These results all point to hierarchical tree diagrams as being a promising additional tool for studying ISM turbulence and star forming regions for obtaining information on the degree of self-gravity, the Mach numbers and the complicated relationship between PPV and PPP data. [ABSTRACT FROM AUTHOR]

Faraday rotation of linearly polarized radio signals provides a very sensitive probe of fluctuations in the interstellar magnetic field and ionized gas density resulting from magnetohydrodynamic (MHD) turbulence. We used a set of statistical tools to analyze images of the spatial gradient of linearly polarized radio emission (|∇P|) for both observational data from a test image of the Southern Galactic Plane Survey (SGPS) and isothermal three-dimensional simulations of MHD turbulence. Visually, in both observations and simulations, a complex network of filamentary structures is seen. Our analysis shows that the filaments in |∇P| can be produced both by interacting shocks and random fluctuations characterizing the non-differentiable field of MHD turbulence. The latter dominates for subsonic turbulence, while the former is only present in supersonic turbulence. We show that supersonic and subsonic turbulence exhibit different distributions as well as different morphologies in the maps of |∇P|. Particularly, filaments produced by shocks show a characteristic "double jump" profile at the sites of shock fronts resulting from delta function-like increases in the density and/or magnetic field, while those produced by subsonic turbulence show a single jump profile. In order to quantitatively characterize these differences, we use the topology tool known as the genus curve as well as the probability distribution function moments of the image distribution. We find that higher values for the moments correspond to cases of |∇P| with larger sonic Mach numbers. The genus analysis of the supersonic simulations of |∇P| reveals a "swiss cheese" topology, while the subsonic cases have characteristics of a "clump" topology. Based on the analysis of the genus and the higher order moments, the SGPS test region data have a distribution and morphology that match subsonic- to transonic-type turbulence, which confirms what is now expected for the warm ionized medium. [ABSTRACT FROM AUTHOR]

In an effort to characterize the Mach numbers of interstellar medium (ISM) magnetohydrodynamic (MHD) turbulence, we study the probability distribution functions (PDFs) of spatial increments of density, velocity, and magnetic field for 14 ideal isothermal MHD simulations at a resolution of 512³. In particular, we fit the PDFs using the Tsallis function and study the dependency of the fit parameters on the compressibility and magnetization of the gas. We find that the Tsallis function fits PDFs of MHD turbulence well, with fit parameters showing sensitivities to the sonic and Alfvén Mach numbers. For three-dimensional density, column density, and Position-Position-Velocity data, we find that the amplitude and width of the PDFs show a dependency on the sonic Mach number. We also find that the width of the PDF is sensitive to the global Alfvénic Mach number especially in cases where the sonic number is high. These dependencies are also found for mock observational cases, where cloud-like boundary conditions, smoothing, and noise are introduced. The ability of Tsallis statistics to characterize the sonic and Alfvénic Mach numbers of simulated ISM turbulence points to it being a useful tool in the analysis of the observed ISM, especially when used simultaneously with other statistical techniques. [ABSTRACT FROM AUTHOR]

The dynamics of the interstellar medium (ISM) are strongly affected by turbulence, which shows increased anisotropy in the presence of a magnetic field. We expand upon the Esquivel & Lazarian method to estimate the Alfvén Mach number using the structure function anisotropy in velocity centroid data from Position-Position-Velocity maps. We utilize three-dimensional magnetohydrodynamic simulations of fully developed turbulence, with a large range of sonic and Alfvénic Mach numbers, to produce synthetic observations of velocity centroids with observational characteristics such as thermal broadening, cloud boundaries, noise, and radiative transfer effects of carbon monoxide. In addition, we investigate how the resulting anisotropy-Alfvén Mach number dependency found in Esquivel & Lazarian might change when taking the second moment of the Position-Position-Velocity cube or when using different expressions to calculate the velocity centroids. We find that the degree of anisotropy is related primarily to the magnetic field strength (i.e., Alfvén Mach number) and the line-of-sight orientation, with a secondary effect on sonic Mach number. If the line of sight is parallel to up to ≈45 deg off of the mean field direction, the velocity centroid anisotropy is not prominent enough to distinguish different Alfvénic regimes. The observed anisotropy is not strongly affected by including radiative transfer, although future studies should include additional tests for opacity effects. These results open up the possibility of studying the magnetic nature of the ISM using statistical methods in addition to existing observational techniques. [ABSTRACT FROM AUTHOR]