Search

Aims. We present a new sub-grid model, HYACINTH – HYdrogen And Carbon chemistry in the INTerstellar medium in Hydro simulations – for computing the non-equilibrium abundances of H2 and its carbon-based tracers, namely CO, C, and C+, in cosmological simulations of galaxy formation. Methods. The model accounts for the unresolved density structure in simulations using a variable probability distribution function of sub-grid densities and a temperature-density relation. Included is a simplified chemical network that has been tailored for hydrogen and carbon chemistry within molecular clouds and easily integrated into large-scale simulations with minimal computational overhead. As an example, we applied HYACINTH to a simulated galaxy at redshift z ~ 2.5 in post-processing and compared the resulting abundances with observations. Results. The chemical predictions from HYACINTH are in reasonable agreement with high-resolution molecular-cloud simulations at different metallicities. By post-processing a galaxy simulation with HYACINTH, we reproduced the H I − H2 transition as a function of the hydrogen column density NH for both Milky-Way-like and Large-Magellanic-Cloud-like conditions. We also matched the NCO versus NH2 relation inferred from absorption measurements towards Milky-Way molecular clouds, although most of our post-processed regions occupy the same region as (optically) dark molecular clouds in the NCO – NH2 plane. Column density maps reveal that CO is concentrated in the peaks of the H2 distribution, while atomic carbon more broadly traces the bulk of H2 in our post-processed galaxy. Based on both the column density maps and the surface density profiles oŕ the different gas species in the post-processed galaxy, we find that C+ maintains a substantially high surŕace density out to ~10 kpc as opposed to other components that exhibit a higher central concentration. This is similar to the extended [C II] emission ŕound in some recent observations at high redshifts. [ABSTRACT FROM AUTHOR]

The interstellar medium is threaded by a hierarchy of filaments from large scales (∼100 pc) to small scales (∼0.1 pc). The masses and lengths of these nested structures may reveal important constraints for cloud formation and evolution, but it is difficult to investigate from an evolutionary perspective using single observations. In this work, we extract simulated molecular clouds from the 'Cloud Factory' galactic-scale ISM suite in combination with 3D Monte Carlo radiative transfer code polaris to investigate how filamentary structure evolves over time. We produce synthetic dust continuum observations in three regions with a series of snapshots and use the filfinder algorithm to identify filaments in the dust derived column density maps. When the synthetic filaments mass and length are plotted on an mass–length (M – L) plot, we see a scaling relation of L ∝ M 0.45 similar to that seen in observations, and find that the filaments are thermally supercritical. Projection effects systematically affect the masses and lengths measured for the filaments, and are particularly severe in crowded regions. In the filament M – L diagram we identify three main evolutionary mechanisms: accretion, segmentation, and dispersal. In particular we find that the filaments typically evolve from smaller to larger masses in the observational M – L plane, indicating the dominant role of accretion in filament evolution. Moreover, we find a potential correlation between line mass and filament growth rate. Once filaments are actively star forming they then segment into smaller sections, or are dispersed by internal or external forces. [ABSTRACT FROM AUTHOR]

Supernovae (SNe) inject ∼1051 erg in the interstellar medium, thereby shocking and heating the gas. A substantial fraction of this energy is later lost via radiative cooling. We present a post-processing module for the flash code to calculate the cooling radiation from shock-heated gas using collisional excitation data from mappings v. When applying this tool to a simulated SN remnant (SNR), we find that most energy is emitted in the EUV. However, optical emission lines ([O iii ], [N ii ], [S ii ], H α, H β) are usually best observable. Our shock detection scheme shows that [S ii ] and [N ii ] emissions arise from the thin shell surrounding the SNR, while [O iii ], H  |$\rm \alpha$|⁠ , and H  |$\rm \beta$| originate from the volume-filling hot gas inside the SNR bubble. We find that the optical emission lines are affected by the SNR's complex structure and its projection on to the plane of the sky because the escaping line luminosity can be reduced by 10–80 per cent due to absorption along the line of sight. Additionally, the subtraction of contaminating background radiation is required for the correct classification of an SNR on the oxygen or sulphur BPT diagrams. The electron temperature and density obtained from our synthetic observations match well with the simulation but are very sensitive to the assumed metallicity. [ABSTRACT FROM AUTHOR]

We present magnetohydrodynamic (MHD) simulations of the star-forming multiphase interstellar medium (ISM) in stratified galactic patches with gas surface densities Σgas = 10, 30, 50, and 100 |$\mathrm{M_\odot \, pc^{-2}}$|⁠. The silcc project simulation framework accounts for non-equilibrium thermal and chemical processes in the warm and cold ISM. The sink-based star formation and feedback model includes stellar winds, hydrogen-ionizing UV radiation, core-collapse supernovae, and cosmic ray (CR) injection and diffusion. The simulations follow the observed relation between Σgas and the star formation rate surface density ΣSFR. CRs qualitatively change the outflow phase structure. Without CRs, the outflows transition from a two-phase (warm and hot at 1 kpc) to a single-phase (hot at 2 kpc) structure. With CRs, the outflow always has three phases (cold, warm, and hot), dominated in mass by the warm phase. The impact of CRs on mass loading decreases for higher Σgas and the mass loading factors of the CR-supported outflows are of order unity independent of ΣSFR. Similar to observations, vertical velocity dispersions of the warm ionized medium (WIM) and the cold neutral medium (CNM) correlate with the star formation rate as |$\sigma _\mathrm{z} \propto \Sigma _\mathrm{SFR}^a$|⁠ , with a ∼ 0.20. In the absence of stellar feedback, we find no correlation. The velocity dispersion of the WIM is a factor ∼2.2 higher than that of the CNM, in agreement with local observations. For |$\Sigma _\mathrm{SFR} \gtrsim 1.5 \times 10^{-2}\, \mathrm{M}_\odot \, \mathrm{yr}^{-1}\, \mathrm{kpc}^{-2}$| the WIM motions become supersonic. [ABSTRACT FROM AUTHOR]

Modern graphics processing units come with dedicated hardware to perform ray/triangle intersections and bounding volume hierarchy traversal. While the primary use case for this hardware is photorealistic 3-D computer graphics, with careful algorithm design scientists can also use this special-purpose hardware to accelerate general-purpose computations, such as point containment queries. This article explains the principles behind these techniques and their application to vector field visualization of large simulation data using particle tracing. [ABSTRACT FROM AUTHOR]

The treatment of radiative transfer with multiple radiation sources is a critical challenge in simulations of star formation and the interstellar medium (ISM). In this paper, we present the novel TreeRay method for solving general radiative transfer problems, based on reverse ray-tracing combined with tree-based accelerated integration. We implement TreeRay in the adaptive mesh refinement code flash , as a module of the tree solver developed by Wünsch et al. However, the method itself is independent of the host code and can be implemented in any grid-based or particle-based hydrodynamics code. A key advantage of TreeRay is that its computational cost is independent of the number of sources, making it suitable for simulations with many point sources (e.g. massive star clusters) as well as simulations where diffuse emission is important. A very efficient communication and tree-walk strategy enable TreeRay to achieve almost ideal parallel scalings. TreeRay can easily be extended with sub-modules to treat radiative transfer at different wavelengths and to implement related physical processes. Here, we focus on ionizing radiation and use the on-the-spot approximation to test the method and its parameters. The ability to set the tree solver time-step independently enables the speedy calculation of radiative transfer in a multiphase ISM, where the hydrodynamic time-step is typically limited by the sound speed of the hot gas produced in stellar wind bubbles or supernova remnants. We show that complicated simulations of star clusters with feedback from multiple massive stars become feasible with TreeRay. [ABSTRACT FROM AUTHOR]

We characterize in detail the two ∼0.3 pc long filamentary structures found within the subsonic region of Barnard 5. We use combined Robert C. Byrd Green Bank Telescope and Very Large Array observations of the molecular lines NH3(1,1) and (2,2) at a resolution of 1800 au, as well as James Clerk Maxwell Telescope continuum observations at 850 and 450 μm at a resolution of 4400 and 3000 au, respectively. We find that both filaments are highly supercritical with a mean mass per unit length, M/L, of ∼80 M⊙ pc−1 after background subtraction, with local increases reaching values of ∼150 M⊙ pc−1. This would require a magnetic field strength of ∼500 μG to be stable against radial collapse. We extract equidistant cuts perpendicular to the spine of the filament and fit a modified Plummer profile as well as a Gaussian to each of the cuts. The filament widths (deconvolved FWHM) range between 6500 and 7000 au (∼0.03 pc) along the filaments. This equals ∼twice the radius of the flat inner region. We find an anticorrelation between the central density and this flattening radius, suggestive of contraction. Further, we also find a strong correlation between the power-law exponent at large radii and the flattening radius. We note that the measurements of these three parameters fall in a plane and derive their empirical relation. Our high-resolution observations provide direct constraints on the distribution of the dense gas within supercritical filaments showing pre- and protostellar activity. [ABSTRACT FROM AUTHOR]

The degree to which the formation and evolution of clouds and filaments in the interstellar medium is regulated by magnetic fields remains an open question. Yet the fundamental properties of the fields (strength and 3D morphology) are not readily observable. We investigate the potential for recovering magnetic field information from dust polarization, the Zeeman effect, and the Faraday rotation measure (RM) in a SILCC-Zoom magnetohydrodynamic (MHD) filament simulation. The object is analysed at the onset of star formation and it is characterized by a line-mass of about |$\mathrm{\left(M/L\right) \sim 63\ \mathrm{M}_{\odot }\ pc^{-1}}$| out to a radius of |$1\,$|  pc and a kinked 3D magnetic field morphology. We generate synthetic observations via polaris radiative transfer (RT) post-processing and compare with an analytical model of helical or kinked field morphology to help interpreting the inferred observational signatures. We show that the tracer signals originate close to the filament spine. We find regions along the filament where the angular dependence with the line of sight (LOS) is the dominant factor and dust polarization may trace the underlying kinked magnetic field morphology. We also find that reversals in the recovered magnetic field direction are not unambiguously associated to any particular morphology. Other physical parameters, such as density or temperature, are relevant and sometimes dominant compared to the magnetic field structure in modulating the observed signal. We demonstrate that the Zeeman effect and the RM recover the line-of-sight magnetic field strength to within a factor 2.1–3.4. We conclude that the magnetic field morphology may not be unambiguously determined in low-mass systems by observations of dust polarization, Zeeman effect, or RM , whereas the field strengths can be reliably recovered. [ABSTRACT FROM AUTHOR]

We present a new approach to extract the power-law part of a density/column-density probability density function (ρ-pdf/ N -pdf) in star-forming clouds. This approach is based on the mathematical method bPlfit of Virkar & Clauset (2014, Annals of Applied Statistics, 8, 89) and it assesses the power-law part of an arbitrary distribution, without any assumptions about the other parts of this distribution. The slope and deviation point are derived as averaged values as the number of bins is varied. Neither parameter is sensitive to spikes and other local features of the tail. This adapted bPlfit method is applied to two different sets of data from numerical simulations of star-forming clouds at scales 0.5 and 500 pc, and it displays ρ-pdf and N -pdf evolution in agreement with a number of numerical and theoretical studies. Applied to Herschel data on the regions Aquila and Rosette, the method extracts pronounced power-law tails, consistent with those seen in simulations of evolved clouds. [ABSTRACT FROM AUTHOR]

Sources of X-rays such as active galactic nuclei and X-ray binaries are often variable by orders of magnitude in luminosity over time-scales of years. During and after these flares the surrounding gas is out of chemical and thermal equilibrium. We introduce a new implementation of X-ray radiative transfer coupled to a time-dependent chemical network for use in 3D magnetohydrodynamical simulations. A static fractal molecular cloud is irradiated with X-rays of different intensity, and the chemical and thermal evolution of the cloud are studied. For a simulated |$10^5\, \mathrm{M}_\odot$| fractal cloud, an X-ray flux <0.01 erg cm−2 s−1 allows the cloud to remain molecular, whereas most of the CO and H2 are destroyed for a flux of ≥1 erg cm−2 s−1. The effects of an X-ray flare, which suddenly increases the X-ray flux by 105×, are then studied. A cloud exposed to a bright flare has 99 per cent of its CO destroyed in 10–20 yr, whereas it takes >103 yr for 99 per cent of the H2 to be destroyed. CO is primarily destroyed by locally generated far-UV emission from collisions between non-thermal electrons and H2; He+ only becomes an important destruction agent when the CO abundance is already very small. After the flare is over, CO re-forms and approaches its equilibrium abundance after 103–105 yr. This implies that molecular clouds close to Sgr A⋆ in the Galactic Centre may still be out of chemical equilibrium, and we predict the existence of clouds near flaring X-ray sources in which CO has been mostly destroyed but H is fully molecular. [ABSTRACT FROM AUTHOR]

We use magnetohydrodynamical simulations of converging warm neutral medium flows to analyse the formation and global evolution of magnetized and turbulent molecular clouds subject to supernova feedback from massive stars. We show that supernova feedback alone fails to disrupt entire, gravitationally bound, molecular clouds, but is able to disperse smallsized (~10 pc) regions on time-scales of less than 1 Myr. Efficient radiative cooling of the supernova remnant as well as strong compression of the surrounding gas result in non-persistent energy and momentum input from the supernovae. However, if the time between subsequent supernovae is short and they are clustered, large hot bubbles form that disperse larger regions of the parental cloud. On longer time-scales, supernova feedback increases the amount of gas with moderate temperatures (T ≈ 300-3000 K). Despite its inability to disrupt molecular clouds, supernova feedback leaves a strong imprint on the star formation process. We find an overall reduction of the star formation efficiency by a factor of 2 and of the star formation rate by roughly factors of 2-4. [ABSTRACT FROM AUTHOR]

Stars are born in turbulent, magnetized filamentary molecular clouds, typically as members of star clusters. Several remarkable technical advances enable observations of magnetic structure and field strengths across many physical scales, from galactic scales on which giant molecular clouds (GMCs) are assembled, down to the surfaces of magnetized accreting young stars. These are shedding new light on the role of magnetic fields in star formation. Magnetic fields affect the gravitational fragmentation and formation of filamentary molecular clouds, the formation and fragmentation of magnetized disks, and finally to the shedding of excess angular momentum in jets and outflows from both the disks and young stars. Magnetic fields play a particularly important role in angular momentum transport on all of these scales. Numerical simulations have provided an important tool for tracking the complex process of the collapse and evolution of protostellar gas since several competing physical processes are at play - turbulence, gravity, MHD, and radiation fields. This paper focuses on the role of magnetic fields in three crucial regimes of star formation: the formation of star clusters emphasizing fragmentation, disk formation and the origin of early jets and outflows, to processes that control the spin evolution of young stars. [ABSTRACT FROM PUBLISHER]