Researcher profile

S. Peng Oh

S. Peng Oh contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
4topics
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)

preprint2023arXiv

Turbulent Reacceleration of Streaming Cosmic Rays

Subsonic, compressive turbulence transfers energy to cosmic rays (CRs), a process known as non-resonant reacceleration. It is often invoked to explain observed ratios of primary to secondary CRs at $\sim \rm GeV$ energies, assuming wholly diffusive CR transport. However, such estimates ignore the impact of CR self-confinement and streaming. We study these issues in stirring box magnetohydrodynamic (MHD) simulations using Athena++, with field-aligned diffusive and streaming CR transport. For diffusion only, we find CR reacceleration rates in good agreement with analytic predictions. When streaming is included, reacceleration rates depend on plasma $β$. Due to streaming-modified phase shifts between CR and gas variables, they are slower than canonical reacceleration rates in low-$β$ environments like the interstellar medium (ISM) but remain unchanged in high-$β$ environments like the intracluster medium (ICM). We also quantify the streaming energy loss rate in our simulations. For sub-Alfvénic turbulence, it is resolution-dependent (hence unconverged in large scale simulations) and heavily suppressed -- by an order of magnitude -- compared to the isotropic loss rate $v_{A} \cdot \nabla P_{\rm CR} / P_{\rm CR} \sim v_{A}/L_{0}$, due to misalignment between the mean field and isotropic CR gradients. Counterintuitively, and unlike acceleration efficiencies, CR losses are almost independent of magnetic field strength over $β\sim 1-100$ and are, therefore, not the primary factor behind lower acceleration rates when streaming is included. While this paper is primarily concerned with how turbulence affects CRs, in a follow-up paper (Bustard and Oh, in prep), we consider how CRs affect turbulence by diverting energy from the MHD cascade, altering the pathway to gas heating and steepening the turbulent power spectrum.

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

Turbulent Heating in a Stratified Medium

There is considerable evidence for widespread subsonic turbulence in galaxy clusters, most notably from {\it Hitomi}. Turbulence is often invoked to offset radiative losses in cluster cores, both by direct dissipation and by enabling turbulent heat diffusion. However, in a stratified medium, buoyancy forces oppose radial motions, making turbulence anisotropic. This can be quantified via the Froude number ${\rm Fr}$, which decreases inward in clusters as stratification increases. We exploit analogies with MHD turbulence to show that wave-turbulence interactions increase cascade times and reduces dissipation rates $ε\propto {\rm Fr}$. Equivalently, for a given energy injection/dissipation rate $ε$, turbulent velocities $u$ must be higher compared to Kolmogorov scalings. High resolution hydrodynamic simulations show excellent agreement with the $ε\propto {\rm Fr}$ scaling, which sets in for ${\rm Fr} < 0.1$. We also compare previously predicted scalings for the turbulent diffusion coefficient $D \propto {\rm Fr}^2$ and find excellent agreement, for ${\rm Fr} < 1$. However, we find a different normalization, corresponding to stronger diffusive suppression by more than an order of magnitude. Our results imply that turbulent diffusion is more heavily suppressed by stratification, over a much wider radial range, than turbulent dissipation. Thus, the latter potentially dominates. Furthermore, this shift implies significantly higher turbulent velocities required to offset cooling, compared to previous models. These results are potentially relevant to turbulent metal diffusion (which is likewise suppressed), and to planetary atmospheres.

preprint2021arXiv

The Cosmic Ray Staircase: the Outcome of the Cosmic Ray Acoustic Instability

Recently, cosmic rays (CRs) have emerged as a leading candidate for driving galactic winds. Small-scale processes can dramatically affect global wind properties. We run two-moment simulations of CR streaming to study how sound waves are driven unstable by phase-shifted CR forces and CR heating. We verify linear theory growth rates. As the sound waves grow non-linear, they steepen into a quasi-periodic series of propagating shocks; the density jumps at shocks create CR bottlenecks. The depth of a propagating bottleneck depends on both the density jump and its velocity; ΔP_c is smaller for rapidly moving bottlenecks. A series of bottlenecks creates a CR staircase structure, which can be understood from a convex hull construction. The system reaches a steady state between growth of new perturbations, and stair mergers. CRs are decoupled at plateaus, but exert intense forces and heating at stair jumps. The absence of CR heating at plateaus leads to cooling, strong gas pressure gradients and further shocks. If bottlenecks are stationary, they can drastically modify global flows; if their propagation times are comparable to dynamical times, their effects on global momentum and energy transfer are modest. The CR acoustic instability is likely relevant in thermal interfaces between cold and hot gas, as well as galactic winds. Similar to increased opacity in radiative flows, the build-up of CR pressure due to bottlenecks can significantly increase mass outflow rates, by up to an order of magnitude. It seeds unusual forms of thermal instability, and the shocks could have distinct observational signatures.

preprint2020arXiv

Fluid Simulations of Cosmic Ray Modified Shocks

We consider cosmic ray (CR) modified shocks with both streaming and diffusion in the two-fluid description. Previously, numerical codes were unable to incorporate streaming in this demanding regime, and have never been compared against analytic solutions. First, we find a new analytic solution highly discrepant in acceleration efficiency from the standard solution. It arises from bi-directional streaming of CRs away from the subshock, similar to a Zeldovich spike in radiative shocks. Since fewer CRs diffuse back upstream, this results in a much lower acceleration efficiency, typically $\sim 10\%$ as opposed to $\sim 50\%$ found in previous analytic work. At Mach number $\gtrsim 10$, the new solution bifurcates into 3 branches, with efficient, intermediate and inefficient CR acceleration. Our two-moment code (Jiang & Oh 2018) accurately recovers these solutions across the entire parameter space probed, with no ad hoc closure relations. For generic initial conditions, the inefficient branch is the most robust and preferred solution. The intermediate branch is unstable, while the efficient branch appears only when the inefficient branch is not allowed (for CR dominated or high plasma $β$ shocks). CR modified shocks have very long equilibration times ($\sim 1000$ diffusion time) required to develop the precursor, which must be resolved by $\gtrsim 10$ cells for convergence. Non-equilibrium effects, poor resolution and obliquity of the magnetic field all reduce CR acceleration efficiency. Shocks in galaxy scale simulations will generally contribute little to CR acceleration without a subgrid prescription.

preprint2011arXiv

On Lyman-limit Systems and the Evolution of the Intergalactic Ionizing Background

We study the properties of self-shielding intergalactic absorption systems and their implications for the ionizing background. We find that cosmological simulations post-processed with detailed radiative transfer calculations generally are able to reproduce the observed abundance of Lyman-limit systems, and we highlight possible discrepancies between the observations and simulations. This comparison tests cosmological simulations at overdensities of ~100. Furthermore, we show that the properties of Lyman-limit systems in these simulations, in simple semi-analytic arguments, and as suggested by recent observations indicate that a small change in the ionizing emissivity of the sources would have resulted in a much larger change in the amplitude of the intergalactic HI-ionizing background (with this scaling strengthening with increasing redshift). This strong scaling could explain the rapid evolution in the Lyman-alpha forest transmission observed at z = 6. Our calculations agree with the suggestion of simpler models that the comoving ionizing emissivity was constant or even increasing from z = 3 to 6. Our calculations also provide a more rigorous estimate than in previous studies for the clumping factor of intergalactic gas after reionization, which we estimate was 2 - 3 at z = 6.

preprint2009arXiv

Accretion onto Seed Black Holes in the First Galaxies

The validity of the hypothesis that the massive black holes in high redshift quasars grew from stellar-sized &#34;seeds&#34; is contingent on a seed&#39;s ability to double its mass every few ten million years. This requires that the seed accrete at approximately the Eddington-limited rate. In the specific case of radiatively efficient quasiradial accretion in a metal-poor protogalactic medium, for which the Bondi accretion rate is often prescribed in cosmological simulations of massive black hole formation, we examine the effects of the radiation emitted near the black hole&#39;s event horizon on the structure of the surrounding gas flow. We find that the radiation pressure from photoionization significantly reduces the steady-state accretion rate and renders the quasiradial accretion flow unsteady and inefficient. The time-averaged accretion rates are a small fraction of the Eddington-limited accretion rate for Thomson scattering. The pressure of Ly-alpha photons trapped near the HII region surrounding the black hole may further attenuate the inflow. These results suggest that an alternative to quasiradial, radiatively efficient Bondi-like accretion should be sought to explain the rapid growth of quasar-progenitor seed black holes.

preprint2006arXiv

Model Atmospheres for Irradiated Giant Stars: Implications for the Galactic Center

Irradiation of a stellar atmosphere by an external source (e.g. an AGN) changes its structure and therefore its spectrum. Using a state-of-the-art stellar atmosphere code, we calculate the infrared spectra of such irradiated and transformed stars. We show that the original spectrum of the star, which is dominated by molecular bands, changes dramatically when irradiated even by a low-luminosity AGN ($L_{\rm X} = 10^{33}$ erg s$^{-1}$), becoming dominated by atomic lines in absorption. We study the changes in the spectrum of low-mass carbon- and oxygen-rich giant stars as they are irradiated by a modest AGN, similar to the one at the Galactic center (GC). The resulting spectra are similar to those of the faintest S-cluster stars observed in the GC. The spectrum of a star irradiated by a much brighter AGN, like that powered by a tidally disrupted star, is very different from that of any star currently observed near the GC. For the first time we have discovered that the structure of the atmosphere of an irradiated giant changes dramatically and induces a double inversion layer. We show that irradiation at the current level can explain the observed trend of CO band intensities decreasing as a function of increasing proximity to Sg $A^{*}$. This may indicate that (contrary to previous claims) there is no paucity of old giants in the GC, which coexist simultaneously with young massive stars.

preprint2004arXiv

Quasar Feedback: the Missing Link in Structure Formation

We consider the impact of quasar outflows on structure formation. Such outflows are potentially more important than galactic winds, which appear insufficient to produce the level of preheating inferred from X-ray observations of galaxy clusters. Using a simple analytical model for the distribution of quasars with redshift, coupled with a one-dimensional Sedov-Taylor model for outflows, we are able to make robust statements about their impact of on structure formation. As large regions of the IGM are heated above a critical entropy of approximately 100 keV cm^2, cooling become impossible within them, regardless of changes in density. On quasar scales, this has the effect of inhibiting further formation, resulting in the observed fall-off in their number densities below z = 2. On galaxy scales, quasar feedback fixes the turn-over scale in the galaxy luminosity function (L_*) as the nonlinear scale at the redshift of strong feedback. The galaxy luminosity function then remains largely fixed after this epoch, consistent with recent observations and in contrast to the strong evolution predicted in more standard galaxy-formation models. Finally, strong quasar feedback explains why the intracluster medium is observed to have been pre-heated to entropy levels just above S_crit, the minimum excess that would not have been erased by cooling. The presence of such outflows is completely consistent with the observed properties of the Lyman-alpha forest at z ~ 2, but is expected to have a substantial and detectable impact on Compton distortions observed in the microwave background and the multiphase properties of the &#34;warm-hot&#34; (z=0) circumgalactic medium.