skip to main content

Title: Microphysics of Relativistic Collisionless Electron-ion-positron Shocks

We perform particle-in-cell simulations to elucidate the microphysics of relativistic weakly magnetized shocks loaded with electron-positron pairs. Various external magnetizationsσ≲ 10−4and pair-loading factorsZ±≲ 10 are studied, whereZ±is the number of loaded electrons and positrons per ion. We find the following: (1) The shock becomes mediated by the ion Larmor gyration in the mean field whenσexceeds a critical valueσLthat decreases withZ±. AtσσLthe shock is mediated by particle scattering in the self-generated microturbulent fields, the strength and scale of which decrease withZ±, leading to lowerσL. (2) The energy fraction carried by the post-shock pairs is robustly in the range between 20% and 50% of the upstream ion energy. The mean energy per post-shock electron scales asE¯eZ±+11. (3) Pair loading suppresses nonthermal ion acceleration at magnetizations as low asσ≈ 5 × 10−6. The ions then become essentially thermal with mean energyE¯i, while electrons form a nonthermal tail, extending fromEZ±+11E¯itoE¯i. Whenσ= 0, particle acceleration is enhanced by the formation of intense magnetic cavities that populate the precursor during the late stages of shock evolution. Here, the maximum energy of the nonthermal ions and electrons keeps growing over the duration of the simulation. Alongside the simulations, we develop theoretical estimates consistent with the numerical results. Our findings have important implications for models of early gamma-ray burst afterglows.

more » « less
Award ID(s):
Author(s) / Creator(s):
; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Medium: X Size: Article No. 74
["Article No. 74"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    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 ratiomi/me= 200, fast Mach numberMms=1–4, and upstream magnetic field angleθBn= 55°–85° from the shock normalnˆ. 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Δϕ/(0.5miush2)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: aMs=4(MA=1.8) case with a long, 30diprecursor of whistler waves alongnˆ, and aMs=7(MA=3.2) case with a shorter, 5diprecursor of whistlers oblique to bothnˆandB;diis 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 theMs=4,θBn= 65° case,ϕshows a weak dependence on the electron plasma-to-cyclotron frequency ratioωpece, andϕdecreases by a factor of 2 asmi/meis raised to the true proton–electron value of 1836.

    more » « less
  2. Abstract

    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 (x¯HI103) or close to unity (x¯HI1). In particular, the neutral fraction evolution of the IGM at the critical redshift range ofz= 6–7 is poorly constrained. We present new constraints onx¯HIatz∼ 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 ofx¯HI(z=6.3)<0.79±0.04(1σ),x¯HI(z=6.5)<0.87±0.03(1σ), andx¯HI(z=6.7)<0.940.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
  3. Abstract

    We use ALMA observations of CO(2–1) in 13 massive (M*≳ 1011M) 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, anguLarmomentum, and Evolution (SQuIGGLE) program. Early results showed that two poststarburst galaxies host large H2reservoirs 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, withMH2109M. Given their high stellar masses, this mass limit corresponds to an average gas fraction offH2MH2/M*7%or ∼14% using lower stellar masses estimates derived from analytic, exponentially declining star formation histories. The gas fraction correlates with theDn4000 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 H2reservoirs 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, theSQuIGGLEgalaxies require the additional and efficient heating or removal of cold gas to bring their low SFRs in line with standard H2scaling relations.

    more » « less
  4. Abstract

    We measure the CO-to-H2conversion factor (αCO) in 37 galaxies at 2 kpc resolution, using the dust surface density inferred from far-infrared emission as a tracer of the gas surface density and assuming a constant dust-to-metal ratio. In total, we have ∼790 and ∼610 independent measurements ofαCOfor CO (2–1) and (1–0), respectively. The mean values forαCO (2–1)andαCO (1–0)are9.35.4+4.6and4.22.0+1.9Mpc2(Kkms1)1, respectively. The CO-intensity-weighted mean is 5.69 forαCO (2–1)and 3.33 forαCO (1–0). We examine howαCOscales with several physical quantities, e.g., the star formation rate (SFR), stellar mass, and dust-mass-weighted average interstellar radiation field strength (U¯). Among them,U¯, ΣSFR, and the integrated CO intensity (WCO) have the strongest anticorrelation with spatially resolvedαCO. We provide linear regression results toαCOfor all quantities tested. At galaxy-integrated scales, we observe significant correlations betweenαCOandWCO, metallicity,U¯, and ΣSFR. We also find thatαCOin each galaxy decreases with the stellar mass surface density (Σ) in high-surface-density regions (Σ≥ 100Mpc−2), following the power-law relationsαCO(21)Σ0.5andαCO(10)Σ0.2. The power-law index is insensitive to the assumed dust-to-metal ratio. We interpret the decrease inαCOwith increasing Σas a result of higher velocity dispersion compared to isolated, self-gravitating clouds due to the additional gravitational force from stellar sources, which leads to the reduction inαCO. The decrease inαCOat high Σis important for accurately assessing molecular gas content and star formation efficiency in the centers of galaxies, which bridge “Milky Way–like” to “starburst-like” conversion factors.

    more » « less
  5. Abstract

    The repeating fast radio burst FRB 20190520B is an anomaly of the FRB population thanks to its high dispersion measure (DM = 1205 pc cm−3) despite its low redshift ofzfrb= 0.241. This excess has been attributed to a large host contribution of DMhost≈ 900 pc cm−3, far larger than any other known FRB. In this paper, we describe spectroscopic observations of the FRB 20190520B field obtained as part of the FLIMFLAM survey, which yielded 701 galaxy redshifts in the field. We find multiple foreground galaxy groups and clusters, for which we then estimated halo masses by comparing their richness with numerical simulations. We discover two separateMhalo> 1014Mgalaxy clusters atz= 0.1867 and 0.2170 that are directly intersected by the FRB sight line within their characteristic halo radiusr200. Subtracting off their estimated DM contributions, as well that of the diffuse intergalactic medium, we estimate a host contribution ofDMhost=430220+140or280170+140pccm3(observed frame), depending on whether we assume that the halo gas extends tor200or 2 ×r200. This significantly smaller DMhost—no longer the largest known value—is now consistent with Hαemission measures of the host galaxy without invoking unusually high gas temperatures. Combined with the observed FRB scattering timescale, we estimate the turbulent fluctuation and geometric amplification factor of the scattering layer to beF˜G4.511(pc2km)1/3, suggesting that most of the gas is close to the FRB host. This result illustrates the importance of incorporating foreground data for FRB analyses both for understanding the nature of FRBs and to realize their potential as a cosmological probe.

    more » « less