Researcher profile

Suoqing Ji

Suoqing Ji contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

9 published item(s)

preprint2022arXiv

Cosmic-Ray Driven Outflows to Mpc Scales from $L_{\ast}$ Galaxies

We study the effects of cosmic rays (CRs) on outflows from star-forming galaxies in the circum and inter-galactic medium (CGM/IGM), in high-resolution, fully-cosmological FIRE-2 simulations (accounting for mechanical and radiative stellar feedback, magnetic fields, anisotropic conduction/viscosity/CR diffusion and streaming, and CR losses). We showed previously that massive ($M_{\rm halo}\gtrsim 10^{11}\,M_{\odot}$), low-redshift ($z\lesssim 1-2$) halos can have CR pressure dominate over thermal CGM pressure and balance gravity, giving rise to a cooler CGM with an equilibrium density profile. This dramatically alters outflows. Absent CRs, high gas thermal pressure in massive halos "traps" galactic outflows near the disk, so they recycle. With CRs injected in supernovae as modeled here, the low-pressure halo allows "escape" and CR pressure gradients continuously accelerate this material well into the IGM in "fast" outflows, while lower-density gas at large radii is accelerated in-situ into "slow" outflows that extend to $>$Mpc scales. CGM/IGM outflow morphologies are radically altered: they become mostly volume-filling (with inflow in a thin mid-plane layer) and coherently biconical from the disk to $>$Mpc. The CR-driven outflows are primarily cool ($T\sim10^{5}\,$K) and low-velocity. All of these effects weaken and eventually vanish at lower halo masses ($\lesssim 10^{11}\,M_{\odot}$) or higher redshifts ($z\gtrsim 1-2$), reflecting the ratio of CR to thermal+gravitational pressure in the outer halo. We present a simple analytic model which explains all of the above phenomena.

preprint2022arXiv

Effects of Different Cosmic Ray Transport Models on Galaxy Formation

Cosmic rays (CRs) with ~GeV energies can contribute significantly to the energy and pressure budget in the interstellar, circumgalactic, and intergalactic medium (ISM, CGM, IGM). Recent cosmological simulations have begun to explore these effects, but almost all studies have been restricted to simplified models with constant CR diffusivity and/or streaming speeds. Physical models of CR propagation/scattering via extrinsic turbulence and self-excited waves predict transport coefficients which are complicated functions of local plasma properties. In a companion paper, we consider a wide range of observational constraints to identify proposed physically-motivated cosmic-ray propagation scalings which satisfy both detailed Milky Way (MW) and extra-galactic $γ$-ray constraints. Here, we compare the effects of these models relative to simpler 'diffusion+streaming' models on galaxy and CGM properties at dwarf through MW mass scales. The physical models predict large local variations in CR diffusivity, with median diffusivity increasing with galacto-centric radii and decreasing with galaxy mass and redshift. These effects lead to a more rapid dropoff of CR energy density in the CGM (compared to simpler models), in turn producing weaker effects of CRs on galaxy star formation rates (SFRs), CGM absorption profiles and galactic outflows. The predictions of the more physical CR models tend to lie 'in between' models which ignore CRs entirely and models which treat CRs with constant diffusivity.

preprint2022arXiv

Numerical Study of Cosmic Ray Confinement through Dust Resonant Drag Instabilities

We investigate the possibility of cosmic ray (CR) confinement by charged dust grains through resonant drag instabilities (RDIs). We perform magnetohydrodynamic particle-in-cell simulations of magnetized gas mixed with charged dust and cosmic rays, with the gyro-radii of dust and GeV CRs on $\sim\mathrm{AU}$ scales fully resolved. As a first study, we focus on one type of RDI wherein charged grains drift super-Alfv{é}nically, with Lorentz forces strongly dominating over drag forces. Dust grains are unstable to the RDIs and form concentrated columns and sheets, whose scale grows until saturating at the simulation box size. Initially perfectly-streaming CRs are strongly scattered by RDI-excited Alfv{é}n waves, with the growth rate of the CR perpendicular velocity components equaling the growth rate of magnetic field perturbations. These rates are well-predicted by analytic linear theory. CRs finally become isotropized and drift at least at $\sim v_\mathrm{A}$ by unidirectional Alfvén waves excited by the RDIs, with a uniform distribution of the pitch angle cosine $μ$ and a flat profile of the CR pitch angle diffusion coefficient $D_{μμ}$ around $μ= 0$, without the "$90$ degree pitch angle problem." With CR feedback on the gas included, $D_{μμ}$ decreases by a factor of a few, indicating a lower CR scattering rate, because the backreaction on the RDI from the CR pressure adds extra wave damping, leading to lower quasi-steady-state scattering rates. Our study demonstrates that the dust-induced CR confinement can be very important under certain conditions, e.g., the dusty circumgalactic medium around quasars or superluminous galaxies.

preprint2022arXiv

Survival and mass growth of cold gas in a turbulent, multiphase medium

Astrophysical gases are commonly multiphase and highly turbulent. In this work, we investigate the survival and growth of cold gas in such a turbulent, multi-phase medium using three-dimensional hydrodynamical simulations. Similar to previous work simulating coherent flow (winds), we find that cold gas survives if the cooling time of the mixed gas is shorter than the Kelvin-Helmholtz time of the cold gas clump (with some weak additional Mach number dependence). However, there are important differences. Near the survival threshold, the long-term evolution is highly stochastic, and subject to the existence of sufficiently large clumps. In a turbulent flow, the cold gas continuously fragments, enhancing its surface area. This leads to exponential mass growth, with a growth time given by the geometric mean of the cooling and the mixing time. The fragmentation process leads to a large number of small droplets which follow a scale-free $\mathrm{d} N/\mathrm{d} m \propto m^{-2}$ mass distribution, and dominate the area covering fraction. Thus, whilst survival depends on the presence of large `clouds', these in turn produce a `fog' of smaller droplets tightly coupled to the hot phase which are probed by absorption line spectroscopy. We show with the aid of Monte-Carlo simulations that the simulated mass distribution emerges naturally due to the proportional mass growth and the coagulation of droplets. We discuss the implications of our results for convergence criteria of larger scale simulations and observations of the circumgalactic medium.

preprint2022arXiv

Testing Physical Models for Cosmic Ray Transport Coefficients on Galactic Scales: Self-Confinement and Extrinsic Turbulence at GeV Energies

The microphysics of ~GeV cosmic ray (CR) transport on galactic scales remain deeply uncertain, with almost all studies adopting simple prescriptions (e.g. constant-diffusivity). We explore different physically-motivated, anisotropic, dynamical CR transport scalings in high-resolution cosmological FIRE simulations of dwarf and ~$L_{\ast}$ galaxies where scattering rates vary with local plasma properties motivated by extrinsic turbulence (ET) or self-confinement (SC) scenarios, with varying assumptions about e.g. turbulent power spectra on un-resolved scales, Alfven-wave damping, etc. We self-consistently predict observables including $γ$-rays ($L_γ$), grammage, residence times, and CR energy densities to constrain the models. We demonstrate many non-linear dynamical effects (not captured in simpler models) tend to enhance confinement. For example, in multi-phase media, even allowing arbitrary fast transport in neutral gas does not substantially reduce CR residence times (or $L_γ$), as transport is rate-limited by the ionized WIM and 'inner CGM' gaseous halo ($10^{4}-10^{6}$ K gas within 10-30 kpc), and $L_γ$ can be dominated by trapping in small 'patches.' Most physical ET models contribute negligible scattering of ~1-10 GeV CRs, but it is crucial to account for anisotropy and damping (especially of fast modes) or else scattering rates would violate observations. We show that the most widely-assumed scalings for SC models produce excessive confinement by factors >100 in the WIM and inner CGM, where turbulent and Landau damping dominate. This suggests either a breakdown of quasi-linear theory used to derive the CR transport parameters in SC, or that other novel damping mechanisms dominate in intermediate-density ionized gas.

preprint2022arXiv

The impact of cosmic rays on dynamical balance and disk-halo interaction in Lstar disk galaxies

Cosmic rays (CRs) are an important component in the interstellar medium (ISM), but their effect on the dynamics of the disk-halo interface (< 10 kpc from the disk) is still unclear. We study the influence of CRs on the gas above the disk with high-resolution FIRE-2 cosmological simulations of late-type Lstar galaxies at redshift around zero. We compare runs with and without CR feedback (with constant anisotropic diffusion around 3e29 cm^2/s and streaming). Our simulations capture the relevant disk halo interactions, including outflows, inflows, and galactic fountains. Extra-planar gas in all of the runs satisfies dynamical balance, where total pressure balances the weight of the overlying gas. While the kinetic pressure from non-uniform motion (>1-kpc scale) dominates in the midplane, thermal and bulk pressures (or CR pressure if included) take over at large heights. We find that with CR feedback, (1) the warm (1e4K) gas is slowly accelerated by CRs; (2) the hot (> 5e5K) gas scale height is suppressed; (3) the warm-hot (2e4-5e5K) medium becomes the most volume-filling phase in the disk-halo interface. We develop a novel conceptual model of the near-disk gas dynamics in low-redshift Lstar galaxies: With CRs, the disk-halo interface is filled with CR-driven warm winds and hot super-bubbles that are propagating into the CGM with a small fraction falling back to the disk. Without CRs, most outflows from hot superbubbles are trapped by the existing hot halo and gravity, so typically they form galactic fountains.

preprint2020arXiv

But What About... Cosmic Rays, Magnetic Fields, Conduction, & Viscosity in Galaxy Formation

We present a suite of high-resolution cosmological simulations, using the FIRE-2 feedback physics together with explicit treatment of magnetic fields, anisotropic conduction and viscosity, and cosmic rays (CRs) injected by supernovae (including anisotropic diffusion, streaming, adiabatic, hadronic and Coulomb losses). We survey systems from ultra-faint dwarf ($M_{\ast}\sim 10^{4}\,M_{\odot}$, $M_{\rm halo}\sim 10^{9}\,M_{\odot}$) through Milky Way masses, systematically vary CR parameters (e.g. the diffusion coefficient $κ$ and streaming velocity), and study an ensemble of galaxy properties (masses, star formation histories, mass profiles, phase structure, morphologies). We confirm previous conclusions that magnetic fields, conduction, and viscosity on resolved ($\gtrsim 1\,$pc) scales have small effects on bulk galaxy properties. CRs have relatively weak effects on all galaxy properties studied in dwarfs ($M_{\ast} \ll 10^{10}\,M_{\odot}$, $M_{\rm halo} \lesssim 10^{11}\,M_{\odot}$), or at high redshifts ($z\gtrsim 1-2$), for any physically-reasonable parameters. However at higher masses ($M_{\rm halo} \gtrsim 10^{11}\,M_{\odot}$) and $z\lesssim 1-2$, CRs can suppress star formation by factors $\sim 2-4$, given relatively high effective diffusion coefficients $κ\gtrsim 3\times10^{29}\,{\rm cm^{2}\,s^{-1}}$. At lower $κ$, CRs take too long to escape dense star-forming gas and lose energy to hadronic collisions, producing negligible effects on galaxies and violating empirical constraints from $γ$-ray emission. But around $κ\sim 3\times10^{29}\,{\rm cm^{2}\,s^{-1}}$, CRs escape the galaxy and build up a CR-pressure-dominated halo which supports dense, cool ($T\ll 10^{6}$ K) gas that would otherwise rain onto the galaxy. CR heating (from collisional and streaming losses) is never dominant.

preprint2020arXiv

Project AMIGA: The Circumgalactic Medium of Andromeda

Project AMIGA (Absorption Maps In the Gas of Andromeda) is a large ultraviolet Hubble Space Telescope program, which has assembled a sample of 43 QSOs that pierce the circumgalactic medium (CGM) of Andromeda (M31) from R=25 to 569 kpc (25 of them probing gas from 25 kpc to about the virial radius-Rvir = 300 kpc-of M31). Our large sample provides an unparalleled look at the physical conditions and distribution of metals in the CGM of a single galaxy using ions that probe a wide range of gas phases (Si II, Si III, Si IV, C II, C IV, and O VI, the latter being from the Far Ultraviolet Spectroscopic Explorer). We find that Si III and O VI have near unity covering factor maintained all the way out to 1.2Rvir and 1.9Rvir, respectively. We show that Si III is the dominant ion over Si II and Si IV at any R. While we do not find that the properties of the CGM of M31 depend strongly on the azimuth, we show that they change remarkably around 0.3-0.5Rvir, conveying that the inner regions of the CGM of M31 are more dynamic and have more complicated multi-phase gas-structures than at R>0.5Rvir. We estimate the metal mass of the CGM within Rvir as probed by Si II, Si III, and Si IV is 2x10^7 Msun and by O VI is >8x10^7 Msun, while the baryon mass of the 10^4-10^5.5 K gas is ~4x10^10 (Z/0.3 Zsun)^(-1) Msun within Rvir. We show that different zoom-in cosmological simulations of L* galaxies better reproduce the column density profile of O VI with R than Si III or the other studied ions. We find that observations of the M31 CGM and zoom-in simulations of L* galaxies have both lower ions showing higher column density dispersion and dependence on R than higher ions, indicating that the higher ionization structures are larger and/or more broadly distributed.

preprint2020arXiv

Properties of the Circumgalactic Medium in Cosmic Ray-Dominated Galaxy Halos

We investigate the impact of cosmic rays (CRs) on the circumgalactic medium (CGM) in FIRE-2 simulations, for ultra-faint dwarf through Milky Way (MW)-mass halos hosting star-forming (SF) galaxies. Our CR treatment includes injection by supernovae, anisotropic streaming and diffusion along magnetic field lines, collisional and streaming losses, with constant parallel diffusivity $κ\sim3\times10^{29}\,\mathrm{cm^2\ s^{-1}}$ chosen to match $γ$-ray observations. With this, CRs become more important at larger halo masses and lower redshifts, and dominate the pressure in the CGM in MW-mass halos at $z\lesssim 1-2$. The gas in these &#34;CR-dominated&#34; halos differs significantly from runs without CRs: the gas is primarily cool (a few $\sim10^{4}\,$K), and the cool phase is volume-filling and has a thermal pressure below that needed for virial or local thermal pressure balance. Ionization of the &#34;low&#34; and &#34;mid&#34; ions in this diffuse cool gas is dominated by photo-ionization, with O VI columns $\gtrsim 10^{14.5}\,\mathrm{cm^{-2}}$ at distances $\gtrsim 150\,\mathrm{kpc}$. CR and thermal gas pressure are locally anti-correlated, maintaining total pressure balance, and the CGM gas density profile is determined by the balance of CR pressure gradients and gravity. Neglecting CRs, the same halos are primarily warm/hot ($T\gtrsim 10^{5}\,$K) with thermal pressure balancing gravity, collisional ionization dominates, O VI columns are lower and Ne VIII higher, and the cool phase is confined to dense filaments in local thermal pressure equilibrium with the hot phase.