More Like this

1. Electron Heating in 2D Particle-in-cell Simulations of Quasi-perpendicular Low-beta Shocks

   https://doi.org/10.3847/1538-4357/ad1f69  

   Tran, Aaron ; Sironi, Lorenzo ( April 2024 , The Astrophysical Journal)

   We measure the thermal electron energization in 1D and 2D particle-in-cell simulations of quasi-perpendicular, low-beta (β_p= 0.25) collisionless ion–electron shocks with mass ratiom_i/m_e= 200, fast Mach number${\mathcal{M}}_{\mathrm{ms}}=1$–4, and upstream magnetic field angleθ_Bn= 55°–85° from the shock normal$\hat{\mathit{n}}$. It is known that shock electron heating is described by an ambipolar,B-parallel electric potential jump, Δϕ_∥, that scales roughly linearly with the electron temperature jump. Our simulations have$\mathrm{\Delta }{\varphi }_{\parallel }/(0.5m_{\mathrm{i}}{u_{\mathrm{sh}}}^2)\sim 0.1$–0.2 in units of the pre-shock ions’ bulk kinetic energy, in agreement with prior measurements and simulations. Different ways to measureϕ_∥, including the use of de Hoffmann–Teller frame fields, agree to tens-of-percent accuracy. Neglecting off-diagonal electron pressure tensor terms can lead to a systematic underestimate ofϕ_∥in our low-β_pshocks. We further focus on twoθ_Bn= 65° shocks: a${\mathcal{M}}_{\mathrm{s}}=4$(${\mathcal{M}}_{\mathrm{A}}=1.8$) case with a long, 30d_iprecursor of whistler waves along$\hat{\mathit{n}}$, and a${\mathcal{M}}_{\mathrm{s}}=7$(${\mathcal{M}}_{\mathrm{A}}=3.2$) case with a shorter, 5d_iprecursor of whistlers oblique to both$\hat{\mathit{n}}$andB;d_iis the ion skin depth. Within the precursors,ϕ_∥has a secular rise toward the shock along multiple whistler wavelengths and also has localized spikes within magnetic troughs. In a 1D simulation of the${\mathcal{M}}_{\mathrm{s}}=4$,θ_Bn= 65° case,ϕ_∥shows a weak dependence on the electron plasma-to-cyclotron frequency ratioω_pe/Ω_ce, andϕ_∥decreases by a factor of 2 asm_i/m_eis raised to the true proton–electron value of 1836.

    

   more » « less
2. Length and Velocity Scales in Protoplanetary Disk Turbulence

   https://doi.org/10.3847/1538-4357/ad2c89  

   Sengupta, Debanjan ; Cuzzi, Jeffrey N ; Umurhan, Orkan M ; Lyra, Wladimir ( April 2024 , The Astrophysical Journal)

   In the theory of protoplanetary disk turbulence, a widely adopted ansatz, or assumption, is that the turnover frequency of the largest turbulent eddy, Ω_L, is the local Keplerian frequency Ω_K. In terms of the standard dimensionless Shakura–Sunyaevαparameter that quantifies turbulent viscosity or diffusivity, this assumption leads to characteristic length and velocity scales given respectively by$\sqrt{\alpha }H$and$\sqrt{\alpha }c$, in whichHandcare the local gas scale height and sound speed. However, this assumption is not applicable in cases when turbulence is forced numerically or driven by some natural processes such as vertical shear instability. Here, we explore the more general case where Ω_L≥ Ω_Kand show that, under these conditions, the characteristic length and velocity scales are respectively$\sqrt{\alpha /R\prime }H$and$\sqrt{\alpha R\prime }c$, where$R\prime \equiv {\mathrm{\Omega }}_L/{\mathrm{\Omega }}_K$is twice the Rossby number. It follows that$\alpha =\tilde{\alpha }/R\prime$, where$\sqrt{\tilde{\alpha }}c$is the root-mean-square average of the turbulent velocities. Properly allowing for this effect naturally explains the reduced particle scale heights produced in shearing box simulations of particles in forced turbulence, and it may help with interpreting recent edge-on disk observations; more general implications for observations are also presented. For$R\prime >1$, the effective particle Stokes numbers are increased, which has implications for particle collision dynamics and growth, as well as for planetesimal formation.

    

   more » « less
3. Cosmic-Ray Drag and Damping of Compressive Turbulence

   https://doi.org/10.3847/1538-4357/aceef9  

   Bustard, Chad ; Oh, S. Peng ( September 2023 , The Astrophysical Journal)

   While it is well known that cosmic rays (CRs) can gain energy from turbulence via second-order Fermi acceleration, how this energy transfer affects the turbulent cascade remains largely unexplored. Here, we show that damping and steepening of the compressive turbulent power spectrum are expected once the damping time$t_{\mathrm{damp}}\sim \rho v^2/{\dot{E}}_{\mathrm{CR}}\propto E_{\mathrm{CR}}^{-1}$becomes comparable to the turbulent cascade time. Magnetohydrodynamic simulations of stirred compressive turbulence in a gas-CR fluid with diffusive CR transport show clear imprints of CR-induced damping, saturating at${\dot{E}}_{\mathrm{CR}}\sim \widetilde{ϵ}$, where$\widetilde{ϵ}$is the turbulent energy input rate. In that case, almost all of the energy in large-scale motions is absorbed by CRs and does not cascade down to grid scale. Through a Hodge–Helmholtz decomposition, we confirm that purely compressive forcing can generate significant solenoidal motions, and we find preferential CR damping of the compressive component in simulations with diffusion and streaming, rendering small-scale turbulence largely solenoidal, with implications for thermal instability and proposed resonant scattering ofE≳ 300 GeV CRs by fast modes. When CR transport is streaming dominated, CRs also damp large-scale motions, with kinetic energy reduced by up to 1 order of magnitude in realisticE_CR∼E_gscenarios, but turbulence (with a reduced amplitude) still cascades down to small scales with the same power spectrum. Such large-scale damping implies that turbulent velocities obtained from the observed velocity dispersion may significantly underestimate turbulent forcing rates, i.e.,$\widetilde{ϵ}\gg \rho v^3/L$.

    

   more » « less
4. Now You See It, Now You Don’t: Star Formation Truncation Precedes the Loss of Molecular Gas by ∼100 Myr in Massive Poststarburst Galaxies at z ∼ 0.6

   https://doi.org/10.3847/1538-4357/ac3dfa  

   Bezanson, Rachel ; Spilker, Justin S. ; Suess, Katherine A. ; Setton, David J. ; Feldmann, Robert ; Greene, Jenny E. ; Kriek, Mariska ; Narayanan, Desika ; Verrico, Margaret ( February 2022 , The Astrophysical Journal)

   We use ALMA observations of CO(2–1) in 13 massive (M_*≳ 10¹¹M_⊙) poststarburst galaxies atz∼ 0.6 to constrain the molecular gas content in galaxies shortly after they quench their major star-forming episode. The poststarburst galaxies in this study are selected from the Sloan Digital Sky Survey spectroscopic samples (Data Release 14) based on their spectral shapes, as part of the Studying QUenching at Intermediate-z Galaxies: Gas, angu$$\overrightarrow{L}\mathrm{ar}$$momentum, and Evolution ($\mathit{SQuIGG}\vec{L}E$) program. Early results showed that two poststarburst galaxies host large H₂reservoirs despite their low inferred star formation rates (SFRs). Here we expand this analysis to a larger statistical sample of 13 galaxies. Six of the primary targets (45%) are detected, with$M_{{\mathrm{H}}_2}\gtrsim {10}^9$M_⊙. Given their high stellar masses, this mass limit corresponds to an average gas fraction of$〈f_{{\mathrm{H}}_2}\equiv M_{{\mathrm{H}}_2}/M_{*}〉\sim 7\%$or ∼14% using lower stellar masses estimates derived from analytic, exponentially declining star formation histories. The gas fraction correlates with theD_n4000 spectral index, suggesting that the cold gas reservoirs decrease with time since burst, as found in local K+A galaxies. Star formation histories derived from flexible stellar population synthesis modeling support this empirical finding: galaxies that quenched ≲150 Myr prior to observation host detectable CO(2–1) emission, while older poststarburst galaxies are undetected. The large H₂reservoirs and low SFRs in the sample imply that the quenching of star formation precedes the disappearance of the cold gas reservoirs. However, within the following 100–200 Myr, the$\mathit{SQuIGG}\vec{L}E$galaxies require the additional and efficient heating or removal of cold gas to bring their low SFRs in line with standard H₂scaling relations.

    

   more » « less
5. (Nearly) Model-independent Constraints on the Neutral Hydrogen Fraction in the Intergalactic Medium at z ∼ 5–7 Using Dark Pixel Fractions in Lyα and Lyβ Forests

   https://doi.org/10.3847/1538-4357/aca678  

   Jin, Xiangyu ; Yang, Jinyi ; Fan, Xiaohui ; Wang, Feige ; Bañados, Eduardo ; Bian, Fuyan ; Davies, Frederick B. ; Eilers, Anna-Christina ; Farina, Emanuele Paolo ; Hennawi, Joseph F. ; et al ( January 2023 , The Astrophysical Journal)

   Cosmic reionization was the last major phase transition of hydrogen from neutral to highly ionized in the intergalactic medium (IGM). Current observations show that the IGM is significantly neutral atz> 7 and largely ionized byz∼ 5.5. However, most methods to measure the IGM neutral fraction are highly model dependent and are limited to when the volume-averaged neutral fraction of the IGM is either relatively low (${\overline{x}}_{\mathrm{H}\mathrm{I}}\lesssim {10}^{-3}$) or close to unity (${\overline{x}}_{\mathrm{H}\mathrm{I}}\sim 1$). In particular, the neutral fraction evolution of the IGM at the critical redshift range ofz= 6–7 is poorly constrained. We present new constraints on${\overline{x}}_{\mathrm{H}\mathrm{I}}$atz∼ 5.1–6.8 by analyzing deep optical spectra of 53 quasars at 5.73 <z< 7.09. We derive model-independent upper limits on the neutral hydrogen fraction based on the fraction of “dark” pixels identified in the Lyαand Lyβforests, without any assumptions on the IGM model or the intrinsic shape of the quasar continuum. They are the first model-independent constraints on the IGM neutral hydrogen fraction atz∼ 6.2–6.8 using quasar absorption measurements. Our results give upper limits of${\overline{x}}_{\mathrm{H}\mathrm{I}}(z=6.3)<0.79\pm 0.04$(1σ),${\overline{x}}_{\mathrm{H}\mathrm{I}}(z=6.5)<0.87\pm 0.03$(1σ), and${\overline{x}}_{\mathrm{H}\mathrm{I}}(z=6.7)<{0.94}_{-0.09}^{+0.06}$(1σ). The dark pixel fractions atz> 6.1 are consistent with the redshift evolution of the neutral fraction of the IGM derived from Planck 2018.

    

   more » « less