Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interstellar medium”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Chemical abundances in Seyfert galaxies – Ⅸ. Helium abundance estimates

For the first time, the helium abundance relative to hydrogen (He/H), which relied on direct measurements of the electron temperature, has been derived in the narrow line regions (NLRs) from a local sample of Seyfert 2 nuclei. In view of this, optical emission line intensities [3000 < λ(Å) < 7000] of 65 local Seyfert 2 nuclei (z < 0.2), taken from Sloan Digital Sky Survey Data Release 15 and additional compilation from the literature, were considered. Here we used photoionization model grid to derive an Ionization Correction Factor (ICF) for the neutral helium. The application of this ICF indicates that the NLRs of Seyfert 2 present a neutral helium fraction of ~50 percent in relation to the total helium abundance. We find that Seyfert 2 nuclei present helium abundance ranging from 0.60 to 2.50 times the solar value, while ~85 percent of the sample present oversolar abundance values. The derived (He/H)–(O/H) abundance relation from the Seyfert 2 is stepper than that of star-forming regions (SFs) and this difference could be due to excess of helium injected into the interstellar medium by the winds of Wolf–Rayet stars. From a regression to zero metallicity, by using Seyfert 2 estimates combined with SFs estimates, we obtained a primordial helium mass fraction Y p = 0.2441 ± 0.0037, a value in good agreement with the one inferred from the temperature fluctuations of the cosmic microwave background by the Planck Collaboration, i.e. $Y_{p}^{Planck}$ =0.2471 ± 0.0003.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Energy balance and Alfvén Mach numbers in compressible magnetohydrodynamic turbulence with a large-scale magnetic field

Energy equipartition is a powerful theoretical tool for understanding astrophysical plasmas. It is invoked, for example, to measure magnetic fields in the interstellar medium (ISM), as evidence for small-scale turbulent dynamo action, and, in general, to estimate the energy budget of star-forming molecular clouds. Here, in this study, we motivate and explore the role of the volume-averaged root-mean-squared (rms) magnetic coupling term between the turbulent, δB , and large-scale, B 0 , fields, $\langle$(δB·B 0 ) 2 $\rangle$$_{V}^{1/2}$. 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 magnetohydrodynamic (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 magnetized ISM.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A GALAH view of the chemical homogeneity and ages of stellar strings identified in Gaia

ABSTRACT The advent of Gaia has led to the discovery of nearly 300 elongated stellar associations (called ‘strings’) spanning hundreds of parsecs in length and mere tens of parsecs in width. These newfound populations present an excellent laboratory for studying the assembly process of the Milky Way thin disc. In this work, we use data from GALAH DR3 to investigate the chemical distributions and ages of 18 newfound stellar populations, 10 of which are strings and 8 of which are compact in morphology. We estimate the intrinsic abundance dispersions in [X/H] of each population and compare them with those of both their local fields and the open cluster (OC) M 67. We find that all but one of these groups are more chemically homogeneous than their local fields. Furthermore, half of the strings, namely Theias 139, 169, 216, 303, and 309, have intrinsic [X/H] dispersions that range between 0.01 and 0.07 dex in most elements, equivalent to those of many OCs. These results provide important new observational constraints on star formation and the chemical homogeneity of the local interstellar medium (ISM). We investigate each population’s Li and chemical clock abundances (e.g. [Sc/Ba], [Ca/Ba], [Ti/Ba], and [Mg/Y]) and find that the ages suggested by chemistry generally support the isochronal ages in all but six structures. This work highlights the unique advantages that chemistry holds in the study of kinematically related stellar groups.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Using host galaxy spectroscopy to explore systematics in the standardization of Type Ia supernovae

We use stacked spectra of the host galaxies of photometrically identified Type Ia supernovae (SNe Ia) from the Dark Energy Survey (DES) to search for correlations between Hubble diagram residuals and the spectral properties of the host galaxies. Utilizing full spectrum fitting techniques on stacked spectra binned by Hubble residual, we find no evidence for trends between Hubble residuals and properties of the host galaxies that rely on spectral absorption features (<1.3σ), such as stellar population age, metallicity, and mass-to-light ratio. However, we find significant trends between the Hubble residuals and the strengths of [O ii] (4.4σ) and the Balmer emission lines (3σ). These trends are weaker than the well-known trend between Hubble residuals and host galaxy stellar mass (7.2σ) that is derived from broad-band photometry. After light-curve corrections, we see fainter SNe Ia residing in galaxies with larger line strengths. We also find a trend (3σ) between Hubble residual and the Balmer decrement (a measure of reddening by dust) using H β and H γ. The trend, quantified by correlation coefficients, is slightly more significant in the redder SNe Ia, suggesting that the bluer SNe Ia are relatively unaffected by dust in the interstellar medium of the host and that dust contributes to current Hubble diagram scatter impacting the measurement of cosmological parameters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The density distribution and physical origins of intermittency in supersonic, highly magnetized turbulence with diverse modes of driving

The probability density function (PDF) of the logarithmic density contrast, s = ln (ρ/ρ 0 ), with gas density ρ and mean density ρ 0 , for hydrodynamical supersonic turbulence is well known to have significant non-Gaussian (intermittent) features that monotonically increase with the turbulent Mach number, M. By studying the mass- and volume-weighted s-PDF for an ensemble of 36 sub-to-trans-Alfvenic mean-field, supersonic, isothermal turbulence simulations with different modes of driving, relevant to molecular gas in the cool interstellar medium, we show that a more intricate picture emerges for the non-Gaussian nature of s. Using four independent measures of the non-Gaussian components, we find hydrodynamical-like structure in the highly magnetized plasma for M≲4. However, for M≳4, the non-Gaussian signatures disappear, leaving approximately Gaussian s-statistics – exactly the opposite of hydrodynamical turbulence in the high-M limit. We also find that the non-Gaussian components of the PDF increase monotonically with more compressive driving modes. To understand the M≲4 non-Gaussian features, we use one-dimensional pencil beams to explore the dynamics along and across the large-scale magnetic field, B 0 . Here we discuss kinetic, density, and magnetic field fluctuations from the pencil beams, and identify physical sources of non-Gaussian components to the PDF as single, strong shocks coupled to fast magnetosonic compressions that form along B 0 . We discuss the Gaussianization of the M≳4 s-fields through the lens of two phenomenologies: the self-similarity of the s-field and homogenization of the dynamical time-scales between the over- and underdense regions in the compressible gas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

CHANG-ES XXV: H i imaging of nearby edge-on galaxies – Data Release 4

ABSTRACT We present the ${\rm H}\, {\small I}$ distribution of galaxies from the Continuum Haloes in Nearby Galaxies – an EVLA Survey (CHANG-ES). Though the observational mode was not optimized for detecting ${\rm H}\, {\small I}$, we successfully produce ${\rm H}\, {\small I}$ cubes for 19 galaxies. The moment-0 maps from this work are available on CHANG-ES data release website (i.e. https://www.queensu.ca/changes). Our sample is dominated by star-forming, ${\rm H}\, {\small I}$-rich galaxies at distances from 6.27 to 34.1 Mpc. ${\rm H}\, {\small I}$ interferometric images on two of these galaxies (NGC 5792 and UGC 10288) are presented here for the first time, while 12 of our remaining sample galaxies now have better ${\rm H}\, {\small I}$ spatial resolutions and/or sensitivities of intensity maps than those in existing publications. We characterize the average scale heights of the ${\rm H}\, {\small I}$ distributions for a subset of most inclined galaxies (inclination > 80 deg), and compare them to the radio continuum intensity scale heights, which have been derived in a similar way. The two types of scale heights are well correlated, with similar dependence on disc radial extension and star formation rate surface density but different dependence on mass surface density. This result indicates that the vertical distribution of the two components may be governed by similar fundamental physics but with subtle differences.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

State-to-state rate coefficients for HCS+ in rotationally inelastic collisions with H2 at low temperatures

ABSTRACT HCS+ ions have been detected in several regions of the interstellar medium (ISM), but an accurate determination of the chemical-physical conditions in the molecular clouds where this molecule is observed requires detailed knowledge of the collisional rate coefficients with the most common colliders in those environments. In this work, we study the dynamics of rotationally inelastic collisions of HCS+ + H2 at low temperature, and report, for the first time, a set of rate coefficients for this system. We used a recently developed potential energy surface for the HCS+–H2 van der Waals complex and computed state-to-state rotational rate coefficients for the lower rotational states of HCS+ in collision with both para- and ortho-H2, analysing the influence of the computed rate coefficients on the determination of critical densities. Additionally, the computed rate coefficients are compared with those obtained by scaling the ones from HCS+ in collision with He (an approximation that is sometimes used when data is lacking), and large differences are found. Furthermore, the approximation of using the rates for the HCO+ + H2 collision as a rough approximation for those of the HCS+ + H2 system is also evaluated. Finally, the complete set of de-excitation rate coefficients for the lowest 30 rotational states of HCS+ by collision with H2 is reported from 5 to 100 K.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

SILCC – VII. Gas kinematics and multiphase outflows of the simulated ISM at high gas surface densities

ABSTRACT 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.

Astronomy & Astrophysics↗
DOE OSTI2023

How do supernova remnants cool? – I. Morphology, optical emission lines, and shocks

ABSTRACT 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.

Astronomy & Astrophysics↗
DOE OSTI2023

Growth or Decay – I: universality of the turbulent dynamo saturation

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 kin into magnetic energy E mag , and is often used to explain how an initially weak magnetic field with E mag $\ll$ E kin is amplified, and then maintained at a level E mag ≲ E 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, here 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 mag ~ E kin or even if initially E mag $\gg$ E 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.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A resolved study of the inner regions of nearby galaxies with an excess of young massive stars: missing link in the AGN–starburst connection?

ABSTRACT We have selected galaxies with very high levels of H α emission (EQW(H α) >700 Å.) in their central regions from the final data release of the MaNGA survey. Our study focuses on 14 very well-resolved nearby galaxies with stellar masses in the range 9.5 < log M*/(M⊙) < 11.5. We investigate a variety of procedures for selecting galaxy regions that are likely to harbour excess populations of young massive stars, finding that selection in the 2D space of extinction-corrected H α EQW and [S iii]/[S ii] line ratio produces the best results. By comparing stacked spectra covering these regions with stacked spectra covering normal starburst regions with 100 Å<EQW(H α) <200 Å, we obtain the following main results: (1) Clear signatures of excess Wolf–Rayet stars are found in half of the H α excess regions, (2) galaxy regions containing excess Wolf–Rayet stars are more often associated with the presence of high-ionization emission lines characteristic of accreting black holes. Excess [Ne iii] is detected in 4 out of 8 of the WR regions and there are tentative [Fe x] detections in 2 galaxies. (3) Regions of the galaxy with excess Wolf–Rayet stars are located where the interstellar medium has larger ionized gas turbulent velocities and higher neutral gas overdensities. We make a first attempt to constrain changes in the high-mass end of the stellar initial mass function (IMF) using the HR-pyPopStar evolutionary population synthesis models that include high-wavelength-resolution theoretical atmosphere libraries for Wolf–Rayet stars.

Kauffmann, Guinevere↗
DOE OSTI2023

Tree-based solvers for adaptive mesh refinement code FLASH – III: a novel scheme for radiation pressure on dust and gas and radiative transfer from diffuse sources

Radiation is an important contributor to the energetics of the interstellar medium, yet its transport is difficult to solve numerically. We present a novel approach towards solving radiative transfer of diffuse sources via backwards ray tracing. Here, we focus on the radiative transfer of infrared radiation and the radiation pressure on dust. The new module, TREERAY/RADPRESSURE, is an extension to the novel radiative transfer method TREERAY implemented in the grid-based Magneto-Hydrodynamics code FLASH. In TREERAY/RADPRESSURE, every cell and every star particle is a source of infrared radiation. We also describe how gas, dust, and radiation are coupled via a chemical network. This allows us to compute the local dust temperature in thermal equilibrium, leading to a significantly improvement over the classical grey approximation. In several tests, we demonstrate that the scheme produces the correct radiative intensities as well as the correct momentum input by radiation pressure. Subsequently, we apply our new scheme to model massive star formation from a collapsing, turbulent core of 150 M ⊙ . We include the effects of both, ionizing and infrared radiation on the dynamics of the core. We find that the newborn massive star prevents fragmentation in its proximity due to radiative heating. Over time, dust and radiation temperature equalize, while the gas temperature can be either warmer due to shock heating or colder due to insufficient dust–gas coupling. Compared to gravity, the effects of radiation pressure are insignificant for the stellar mass on the simulated time-scale in this work.

(ISM:) HII regions↗
DOE OSTI2023

Chemical abundances in Seyfert galaxies – X. Sulphur abundance estimates

For the first time, the sulphur abundance relative to hydrogen (S/H) in the narrow-line regions of a sample of Seyfert 2 nuclei (Sy 2s) has been derived via direct estimation of the electron temperature. Narrow emission-line intensities from the Sloan Digital Sky Survey (SDSS) Data Release 17 (DR17) [in the wavelength range 3000 < λ(Å) < 9100] and from the literature for a sample of 45 nearby (z < 0.08) Sy 2s were considered. Our direct estimates indicate that Sy 2s have similar temperatures in the gas region where most of the $\rm S\rm ^{+}$ ions are located in comparison with that of star-forming regions (SFs). However, Sy 2s present higher temperature values (${\sim} 10\, 000$ K) in the region where most of the $\rm S\rm ^{2+}$ ions are located relative to that of SFs. We derive the total sulphur abundance in the range of $6.2 \: \lesssim 12+\log (\rm S/H) \: \lesssim \: 7.5$, corresponding to 0.1–1.8 times the solar value. These sulphur abundance values are lower by ∼0.4 dex than those derived in SFs with similar metallicity, indicating a distinct chemical enrichment of the interstellar medium (ISM) for these object classes. The sulphur abundance relative to oxygen (S/O) values for our Sy 2 sample present an abrupt (∼0.5 dex) decrease with increasing oxygen abundance relative to hydrogen (O/H) for the high-metallicity regime [$\rm 12+\log (O/H) \: \gtrsim 8.7)$], what is not seen for the SFs. However, when our Sy 2 estimates are combined with those from a large sample of SFs, we did not find any dependence between S/O and O/H.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The rate of extreme coronal line emitting galaxies in the Sloan Digital Sky Survey and their relation to tidal disruption events

High-ionization iron coronal lines (CLs) are a rare phenomenon observed in galaxy and quasi-stellar object spectra that are thought to be created by high-energy emission from active galactic nuclei and certain types of transients. In cases known as extreme coronal line emitting galaxies (ECLEs), these CLs are strong and fade away on a time-scale of years. The most likely progenitors of these variable CLs are tidal disruption events (TDEs), which produce sufficient high-energy emission to create and sustain the CLs over these time-scales. To test the possible connection between ECLEs and TDEs, we present the most complete variable ECLE rate calculation to date and compare the results to TDE rates from the literature. To achieve this, we search for ECLEs in the Sloan Digital Sky Survey (SDSS). We detect sufficiently strong CLs in 16 galaxies, more than doubling the number previously found in SDSS. Using follow-up spectra from the Dark Energy Spectroscopic Instrument and Gemini Multi-Object Spectrograph, Wide-field Infrared Survey Explorer mid-infrared observations, and Liverpool Telescope optical photometry, we find that none of the nine new ECLEs evolve in a manner consistent with that of the five previously discovered variable ECLEs. Using this sample of five variable ECLEs, we calculate the galaxy-normalized rate of variable ECLEs in SDSS to be $R_\mathrm{G}=3.6~^{+2.6}_{-1.8}~(\mathrm{statistical})~^{+5.1}_{-0.0}~(\mathrm{systematic})\times 10^{-6}~\mathrm{galaxy}^{-1}~\mathrm{yr}^{-1}$. The mass-normalized rate is $R_\mathrm{M}=3.1~^{+2.3}_{-1.5}~(\mathrm{statistical})~^{+4.4}_{-0.0}~(\mathrm{systematic})\times 10^{-17}~\mathrm{M_\odot ^{-1}}~\mathrm{yr}^{-1}$ and the volumetric rate is $R_\mathrm{V}=7~^{+20}_{-5}~(\mathrm{statistical})~^{+10}_{-0.0}~(\mathrm{systematic})\times 10^{-9}~\mathrm{Mpc}^{-3}~\mathrm{yr}^{-1}$. Our rates are one to two orders of magnitude lower than TDE rates from the literature, which suggests that only 10–40 per cent of all TDEs produce variable ECLEs. Additional uncertainties in the rates arising from the structure of the interstellar medium have yet to be included.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

SDSS-IV MaNGA: calibration of astrophysical line-widths in the Hα region using HexPak observations

ABSTRACT We have re-observed $\rm \sim$40 low-inclination, star-forming galaxies from the MaNGA survey (σ ∼ 65 km s−1) at ∼6.5 times higher spectral resolution (σ ∼ 10 km s−1) using the HexPak integral field unit on the WIYN 3.5-m telescope. The aim of these observations is to calibrate MaNGA’s instrumental resolution and to characterize turbulence in the warm interstellar medium and ionized galactic outflows. Here we report the results for the Hα region observations as they pertain to the calibration of MaNGA’s spectral resolution. Remarkably, we find that the previously reported MaNGA line-spread-function (LSF) Gaussian width is systematically underestimated by only 1 per cent. The LSF increase modestly reduces the characteristic dispersion of H ii regions-dominated spectra sampled at 1–2 kpc spatial scales from 23 to 20 km s−1 in our sample, or a 25 per cent decrease in the random-motion kinetic energy. This commensurately lowers the dispersion zeropoint in the relation between line-width and star-formation rate surface-density in galaxies sampled on the same spatial scale. This modest zero-point shift does not appear to alter the power-law slope in the relation between line-width and star-formation rate surface-density. We also show that adopting a scheme whereby corrected line-widths are computed as the square root of the median of the difference in the squared measured line width and the squared LSF Gaussian avoids biases and allows for lower signal-to-noise data to be used reliably.

Chattopadhyay, Sabyasachi↗
DOE OSTI2024

Detecting rare neutral atomic-carbon absorbers with a deep neural network

ABSTRACT C i absorbers play an important role as indicators for exploring the presence of cold gas in the interstellar medium of galaxies. However, the current data base of C i absorbers is very limited due to their weak absorption feature and rarity. Here, we report results from a search of C i λλ1560, 1656 absorption lines using Mg ii absorbers as signposts with modified deep learning algorithms, which provides a very quick way to search for weak C i absorber candidates. A total of 107 C i absorbers were detected, which nearly doubles the size of previously known samples. In addition, we found 17 C i absorbers to be associated with 2175 Å dust absorbers (2DAs), i.e. about 16 per cent C i absorbers are associated with 2DAs. Comparing the average dust depletion patterns of C i absorbers with those of damped Lyman α absorbers (DLAs), Mg ii absorbers, Ca ii absorbers, and 2175 Å dust absorbers (2DAs) shows that C i absorbers generally have environments with more dust than DLAs, Mg ii, and Ca ii absorbers, but similar to dust in 2DAs. Similarity between the dust depletion pattern of C i absorbers to that of the warm disc in the Milky Way indicates that C i absorption clouds are possibly associated with disc components in distant galaxies. Therefore, C i absorbers are confirmed to be excellent probes to trace cold gas and dust in the Universe.

Ge, Jian↗
DOE OSTI2025

Taking control of compressible modes: bulk viscosity and the turbulent dynamo

Many polyatomic astrophysical plasmas are compressible and out of chemical and thermal equilibrium, introducing a bulk viscosity into the plasma via the internal degrees of freedom of the molecular composition, directly impacting the decay of compressible modes, $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$. This is especially important for small-scale, turbulent dynamo processes in the interstellar medium (ISM), which are known to be sensitive to the effects of compression. To control the viscous properties of $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$, we perform trans-sonic, visco-resistive dynamo simulations with additional bulk viscosity $\nu _{\text{bulk}}$, deriving a new $\nu _{\text{bulk}}$ Reynolds number $\text{Re}_{\text{bulk}}$, and viscous Prandtl number $\text{P}\nu \equiv \text{Re}_{\text{bulk}}/ \text{Re}_{\text{shear}}$, where $\text{Re}_{\text{shear}}$ is the shear viscosity Reynolds number. We derive a framework for decomposing $E_{\rm mag}$ growth rates into incompressible and compressible terms via orthogonal tensor decompositions of $\boldsymbol {\nabla }\otimes \mathrm{{\boldsymbol {\mathit {v}}}}$, where $\mathrm{{\boldsymbol {\mathit {v}}}}$ is the fluid velocity. We find that $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$ play a dual role, growing and decaying $E_{\rm mag}$, and that field-line stretching is the main driver of growth, even in compressible dynamos. In the absence of $\nu _{\text{bulk}}$ ($\text{P}\nu \rightarrow \infty$), $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$ pile up on small-scales, creating a spectral bottleneck, which disappears for $\text{P}\nu \approx 1$. As $\text{P}\nu$ decreases, $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$ are dissipated at increasingly larger scales, in turn suppressing incompressible modes through a coupling between high-k modes. We emphasize the importance of further understanding the role of $\nu _{\text{bulk}}$ in compressible astrophysical plasmas, which we estimate could be as strong as the shear viscosity in the cold ISM, and highlight that compressible direct numerical simulations without bulk viscosity have unresolved compressible mode dissipation scales.

MHD↗
DOE OSTI2025

Rotational excitation and de-excitation of magnesium mono-sulphide (MgS) by He collisions

ABSTRACT Magnesium mono-sulphide (MgS) plays a crucial role in astrochemical processes, particularly in the interstellar medium (ISM), where metal-sulphur chemistry influences molecular formation and evolution. This work presents a detailed study of the rotational excitation and de-excitation of MgS through collisions with helium (He) atoms, the second most abundant species in the ISM after hydrogen, which includes both atomic (H) and molecular forms (H2). The focus on MgS–He collisions arises from He's high abundance, chemical inertness, and simpler electronic structure, which make it well suited for quantum scattering calculations. These characteristics establish He as an ideal candidate for initial studies, providing fundamental data for future investigations involving H2. The study uses quantum scattering methods to calculate the collisional rate coefficients over a broad temperature range. These rates are critical for interpreting observational data on MgS and predicting its abundance in space. The interaction potential between MgS and He is calculated using the rigid rotor approximation and the Jacobi coordinate system, employing the CCSD(T)-F12a/aug-cc-pVTZ method for accurate two-dimensional potential energy surface. The study explores the anisotropic nature of the MgS–He interaction, which favours odd Δj rotational transitions at low collision energies. The inelastic cross-sections for rotational transitions involving up to 16 rotational levels of MgS were computed up to 1000 cm−1, enabling the calculation of rate coefficients up to 150 K for Δj = ±1, ±2, and ± 3 rotational transitions. The results show that Δj = 1 transitions dominate at low temperatures, while Δj = 2 transitions become more significant at higher temperatures. This study provides valuable data for interpreting future astrophysical observations of MgS. The findings also propose new rotational transitions for MgS detection in space, enhancing our ability to track and study this molecule in various cosmic environments.

Hendaoui, Hamza (ORCID:0000000218641872)↗