Abstract We present a new suite of numerical simulations of the star-forming interstellar medium (ISM) in galactic disks using the TIGRESS-NCR framework. Distinctive aspects of our simulation suite are (1) sophisticated and comprehensive numerical treatments of essential physical processes including magnetohydrodynamics, self-gravity, and galactic differential rotation, as well as photochemistry, cooling, and heating coupled with direct ray-tracing UV radiation transfer and resolved supernova feedback and (2) wide parameter coverage including the variation in metallicity over , gas surface density Σgas∼ 5–150M⊙pc−2, and stellar surface density Σstar∼ 1–50M⊙pc−2. The range of emergent star formation rate surface density is ΣSFR∼ 10−4–0.5M⊙kpc−2yr−1, and ISM total midplane pressure isPtot/kB= 103–106cm−3K, withPtotequal to the ISM weight . For given Σgasand Σstar, we find . We provide an interpretation based on the pressure-regulated feedback-modulated (PRFM) star formation theory. The total midplane pressure consists of thermal, turbulent, and magnetic stresses. We characterize feedback modulation in terms of the yield ϒ, defined as the ratio of each stress to ΣSFR. The thermal feedback yield varies sensitively with both weight and metallicity as , while the combined turbulent and magnetic feedback yield shows weaker dependence . The reduction in ΣSFRat low metallicity is due mainly to enhanced thermal feedback yield, resulting from reduced attenuation of UV radiation. With the metallicity-dependent calibrations we provide, PRFM theory can be used for a new subgrid star formation prescription in cosmological simulations where the ISM is unresolved.
more »
« less
Planck Dust Polarization Power Spectra Are Consistent with Strongly Supersonic Turbulence
Abstract The polarization of the cosmic microwave background is rich in information but obscured by foreground emission from the Milky Way’s interstellar medium (ISM). To uncover relationships between the underlying turbulent ISM and the foreground power spectra, we simulated a suite of driven, magnetized, turbulent models of the ISM, varying the fluid properties via the sonic Mach number, , and magnetic (Alfvén) Mach number, . We measure the power spectra of density (ρ), velocity (v), magnetic field (H), total projected intensity (T), parity-even polarization (E), and parity-odd polarization (B). We find that the slopes of all six quantities increase with . Most increase with , while the magnetic field spectrum steepens with . By comparing spectral slopes ofEandBto those measured by Planck, we infer typical values of and for the ISM. As the fluid velocity increases, , the ratio of BB power to EE power increases to approach a constant value near the Planck-observed value of ∼0.5, regardless of the magnetic field strength. We also examine correlation coefficients between projected quantities, and find thatrTE≈ 0.3, in agreement with Planck, for appropriate combinations of and . Finally, we consider parity-violating correlationsrTBandrEB.
more »
« less
- Award ID(s):
- 2009870
- PAR ID:
- 10535095
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 972
- Issue:
- 1
- ISSN:
- 0004-637X
- Format(s):
- Medium: X Size: Article No. 26
- Size(s):
- Article No. 26
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract We present a survey of 1D kinetic particle-in-cell simulations of quasi-parallel nonrelativistic shocks to identify the environments favorable for electron acceleration. We explore an unprecedented range of shock speedsvsh≈ 0.067–0.267c, Alfvén Mach numbers , sonic Mach numbers , as well as the proton-to-electron mass ratiosmi/me= 16–1836. We find that high Alfvén Mach number shocks can channel a large fraction of their kinetic energy into nonthermal particles, self-sustaining magnetic turbulence and acceleration to larger and larger energies. The fraction of injected particles is ≲0.5% for electrons and ≈1% for protons, and the corresponding energy efficiencies are ≲2% and ≈10%, respectively. The extent of the nonthermal tail is sensitive to the Alfvén Mach number; when , the nonthermal electron distribution exhibits minimal growth beyond the average momentum of the downstream thermal protons, independently of the proton-to-electron mass ratio. Acceleration is slow for shocks with low sonic Mach numbers, yet nonthermal electrons still achieve momenta exceeding the downstream thermal proton momentum when the shock Alfvén Mach number is large enough. We provide simulation-based parameterizations of the transition from thermal to nonthermal distribution in the downstream (found at a momentum around ), as well as the ratio of nonthermal electron to proton number density. The results are applicable to many different environments and are important for modeling shock-powered nonthermal radiation.more » « less
-
Abstract The abundance of faint dwarf galaxies is determined by the underlying population of low-mass dark matter (DM) halos and the efficiency of galaxy formation in these systems. Here, we quantify potential galaxy formation and DM constraints from future dwarf satellite galaxy surveys. We generate satellite populations using a suite of Milky Way (MW)–mass cosmological zoom-in simulations and an empirical galaxy–halo connection model, and assess sensitivity to galaxy formation and DM signals when marginalizing over galaxy–halo connection uncertainties. We find that a survey of all satellites around one MW-mass host can constrain a galaxy formation cutoff at peak virial masses of at the 1σlevel; however, a tail toward low prevents a 2σmeasurement. In this scenario, combining hosts with differing bright satellite abundances significantly reduces uncertainties on at the 1σlevel, but the 2σtail toward low persists. We project that observations of one (two) complete satellite populations can constrain warm DM models withmWDM≈ 10 keV (20 keV). Subhalo mass function (SHMF) suppression can be constrained to ≈70%, 60%, and 50% that in cold dark matter (CDM) at peak virial masses of 108, 109, and 1010M⊙, respectively; SHMF enhancement constraints are weaker (≈20, 4, and 2 times that in CDM, respectively) due to galaxy–halo connection degeneracies. These results motivate searches for faint dwarf galaxies beyond the MW and indicate that ongoing missions like Euclid and upcoming facilities including the Vera C. Rubin Observatory and Nancy Grace Roman Space Telescope will probe new galaxy formation and DM physics.more » « less
-
Abstract The Kruskal–Szekeres coordinate construction for the Schwarzschild spacetime could be interpreted simply as a squeezing of thet-line into a single point, at the event horizon . Starting from this perspective, we extend the Kruskal charting to spacetimes with two horizons, in particular the Reissner–Nordström manifold, . We develop a new method to construct Kruskal-like coordinates through casting the metric in new null coordinates, and find two algebraically distinct ways to chart , referred to as classes: type-I and type-II within this work. We pedagogically illustrate our method by crafting two compact, conformal, and global coordinate systems labeled and as an example for each class respectively, and plot the corresponding Penrose diagrams. In both coordinates, the metric differentiability can be promoted to in a straightforward way. Finally, the conformal metric factor can be written explicitly in terms of thetandrfunctions for both types of charts. We also argued that the chart recently reported in Soltani (2023 arXiv:2307.11026) could be viewed as another example for the type-II classification, similar to .more » « less
-
Abstract We consider an -system, that is, a compact Riemannian manifold with boundary, endowed with a magnetic field and a potential. On simple -systems, we study the -ray transform in order to obtain new boundary rigidity results for -systems. We show that there is an explicit relation between the -ray transform and the magnetic one, which allows us to apply results from Dairbekovet al(2007Adv. Math.216535–609) to our case. Regarding rigidity, we show that there exists a generic set of simple -systems, which is open and dense, such that any two -systems close to an element in it and having the same boundary action function, must bek-gauge equivalent.more » « less