We investigate the transport of spectrally resolved cosmic ray (CR) protons with kinetic energies between $1-100$ GeV within the dynamic, multiphase interstellar medium (ISM), using a two-moment CR fluid solver applied to a TIGRESS MHD simulation with conditions similar to the solar neighborhood. Our CR transport prescription incorporates space- and momentum-dependent CR scattering coefficients $\sigma=\kappa^{-1}$, computed from the local balance between streaming-driven Alfvèn wave growth and damping processes. We find that advection combines with momentum-dependent diffusion to produce a CR distribution function $f(p)\propto~p^{-\gamma}$ with $\gamma\approx4.6$ that agrees with observations, steepened from an injected power law slope $\gamma_\mathrm{inj}=4.3$. The CR pressure is uniform in the highly diffusive, mostly neutral midplane region, but decreases exponentially in the ionized extraplanar region where scattering is efficient. To interpret these numerical results, we develop a two-zone analytic model that captures and links the two (physically and spatially) distinct regimes of CR transport in the multiphase, dynamic ISM. At low momenta, CR transport is dominated by gas advection, while at high momenta, both advection and diffusion contribute. At high momentum, the analytic prediction for the spectral slope approaches $\gamma=(4/3)\gamma_\mathrm{inj}-1$, and the predicted scaling of grammage with momentum is $X\propto p^{1-\gamma_\mathrm{inj}/3}$, consistent with the simulations. These results support a physical picture in which CRs are confined within the neutral midplane by the surrounding ionized gas, with their escape regulated by both the CR scattering rate in the ionized extraplanar gas and the velocity and Alfvén speed of that gas, at effective speed $v_\mathrm{c,eff}\approx(1/2)[\kappa_\parallel~d(v+v_\mathrm{A,i})/dz]^{1/2}$.
Radiative cooling can drive dynamics in multi-phase gas. A dramatic example is hydrodynamic `shattering', the violent, pressure-driven fragmentation of a cooling cloud which falls drastically out of pressure balance with its surroundings. We run MHD simulations to understand how shattering is influenced by magnetic fields. In MHD, clouds do not `shatter' chaotically. Instead, after initial fragmentation, both hot and cold phases coherently `stream' in long-lived, field-aligned, self-sustaining gas flows, at high speed ($\sim 100 \, {\rm km \, s^{-1}}$). MHD thermal instability also produces such flows. They are due to the anisotropic nature of MHD pressure support, which only operates perpendicular to B-fields. Thus, even when $P_{\rm B} + P_{\rm gas} \approx$const, pressure balance only holds perpendicular to B-fields. Field-aligned gas pressure variations are unopposed, and results in gas velocities $v \sim (2 \Delta P/\rho)^{1/2}$ from Bernoulli's principle. Strikingly, gas in adjacent flux tubes $\textit{counter-stream}$ in opposite directions. We show this arises from a cooling-induced, MHD version of the thin shell instability. Magnetic tension is important both in enabling corrugational instability and modifying its non-linear evolution. Even in high $\beta$ hot gas, streaming can arise, since magnetic pressure support grows as gas cools and compresses. Thermal conduction increases the sizes and velocities of streaming cloudlets, but does not qualitatively modify dynamics. These results are relevant to the counter-streaming gas flows observed in solar coronal rain, as well as multi-phase gas cooling and condensation in the ISM, CGM and ICM.
We model the transport and spectral evolution of 1-100 GeV cosmic ray (CR) electrons (CREs) in TIGRESS MHD simulations of the magnetized, multiphase interstellar medium. We post-process a kpc-sized galactic disk patch representative of the solar neighborhood using a two-moment method for CR transport that includes advection, streaming, and diffusion. The diffusion coefficient is set by balancing wave growth via the CR streaming instability against wave damping (nonlinear Landau and ion-neutral collisions), depending on local gas and CR properties. Implemented energy loss mechanisms include synchrotron, inverse Compton, ionization, and bremsstrahlung. We evaluate CRE losses by different mechanisms as a function of energy and distance from the midplane, and compare loss timescales to transport and diffusion timescales. This comparison shows that CRE spectral steepening above p = 1 GeV/c is due to a combination of energy-dependent transport and losses. Our evolved CRE spectra are consistent with direct observations in the solar neighborhood, with a spectral index that steepens from an injected value of -2.3 to an energy dependent value between -2.7 and -3.3. We also show that the steepening is independent of the injection spectrum. Finally, we present potential applications of our models, including to the production of synthetic synchrotron emission. Our simulations demonstrate that the CRE spectral slope can be accurately recovered from pairs of radio observations in the range 1.5-45 GHz.
The simulation cost for cosmological simulation-based inference can be decreased by combining simulation sets of varying fidelity. We propose an approach to such multi-fidelity inference based on feature matching and knowledge distillation. Our method results in improved posterior quality, particularly for small simulation budgets and difficult inference problems.
One important question in active galactic nucleus (AGN) is how gas is brought down to the galaxy center. Both internal secular evolution (torque induced by non-axisymmetric galactic structures such as bars) and external processes (e.g. mergers or interactions) are expected to redistribute the angular momentum (AM) and transport gas inward. However, it is still under debate whether these processes can significantly affect AGN activities. Here we for the first time report that AGN fraction increases with the difference of kinematic position angles ($\Delta PA\equiv|PA_{\mathrm{gas}}-PA_{\mathrm{star}}|$) between ionized gas ($PA_{\mathrm{gas}}$) and stellar disks ($PA_{\mathrm{star}}$) in blue and green galaxies, meanwhile this fraction remains roughly constant for red galaxies. Also the high luminosity AGN fraction increases with $\Delta PA$ while the low luminosity AGN fraction is independent with $\Delta PA$. These observational results support a scenario in which the interaction between accreted and pre-existing gas provides the AM loss mechanism, thereby the gas inflow fuels the central BH activities, and the AM loss efficiency is positively correlated with the $\Delta PA$.
SRGA J144459.2$-$604207 is a newly confirmed accreting millisecond X-ray pulsar and type I X-ray burster. We present the broadband X-ray timing and spectral behaviors of SRGA J144459.2$-$604207 during its 2024 outburst. The data were collected from NICER, Einstein Probe, IXPE, Insight-HXMT, NuSTAR and INTEGRAL observations. X-ray pulsations have been detected for the 1.5$-$90 keV energy range throughout the `ON' phase of the outburst from MJD $\sim 60355-60385$. We refined the orbital and spin ephemerides assuming a circular orbit, and found that the pulsar was in a spin-up state during MJD $\sim$ 60361--60377 showing a significant spin-up rate $\dot{\nu}$ of $(3.15\pm 0.36)\times10^{-13}~{\rm Hz~s^{-1}}$. Around MJD $\sim 60377$ a swing was detected in the spin evolution accompanied by significantly enhanced pulsed emission. We studied the pulse profile morphology during the X-ray bursts as observed by Insight-HXMT, IXPE and NuSTAR. During the bursts, pulsations were detected across the 2$-$60 keV with shapes broadly consistent with those observed for the persistent emission. We found, however, that the `burst' pulse profiles exhibit significant phase offsets relative to the pre- and post-burst profiles. These offsets systematically decrease with increasing energy, $\Delta \phi\approx0.15$, 0.11 and 0.02 for IXPE, Insight-HXMT ME and HE in 2$-$8, 5$-$30 and 20$-$60 keV, respectively, and $\Delta \phi\approx0.21$, 0.10 and 0.07 for NuSTAR in 3$-$10, 20$-$35 and 35$-$60 keV, respectively, compared to the pre- and post-burst profiles. We performed a joint spectral analysis of quasi-simultaneous NICER, NuSTAR, and Insight-HXMT data for two epochs. The resulting spectra from both observations were consistent and well-described by an absorbed thermal Comptonization model, nthcomp, plus relativistic reflection, relxillCp.
this https URL ). Comments and endorsements welcomed via the link provided or at sukrit@arizona.edu. Comments/Endorsements received by August 15 will be incorporated into the published version (PASP conference proceedings)
this https URL ). Comments and endorsements welcomed via the link provided or at sukrit@arizona.edu. Comments/Endorsements received by August 15 will be incorporated into the published version (PASP conference proceedings)
this https URL ). Comments are welcome!