More Like this

1. The Origin of Power-law Spectra in Relativistic Magnetic Reconnection

   https://doi.org/10.3847/2041-8213/acfe7c  

   Zhang, Hao ; Sironi, Lorenzo ; Giannios, Dimitrios ; Petropoulou, Maria ( October 2023 , The Astrophysical Journal Letters)

   Magnetic reconnection is often invoked as a source of high-energy particles, and in relativistic astrophysical systems it is regarded as a prime candidate for powering fast and bright flares. We present a novel analytical model—supported and benchmarked with large-scale three-dimensional kinetic particle-in-cell simulations in electron–positron plasmas—that elucidates the physics governing the generation of power-law energy spectra in relativistic reconnection. Particles with Lorentz factorγ≳ 3σ(here,σis the magnetization) gain most of their energy in the inflow region, while meandering between the two sides of the reconnection layer. Their acceleration time is$t_{\mathrm{acc}}\sim \gamma {\eta }_{\mathrm{rec}}^{-1}{\omega }_{\mathrm{c}}^{-1}\simeq 20\gamma {\omega }_{\mathrm{c}}^{-1}$, whereη_rec≃ 0.06 is the inflow speed in units of the speed of light andω_c=eB₀/mcis the gyrofrequency in the upstream magnetic field. They leave the region of active energization aftert_esc, when they get captured by one of the outflowing flux ropes of reconnected plasma. We directly measuret_escin our simulations and find thatt_esc∼t_accforσ≳ few. This leads to a universal (i.e.,σ-independent) power-law spectrum${\mathit{dN}}_{\mathrm{free}}/d\gamma \propto {\gamma }^{-1}$for the particles undergoing active acceleration, and$\mathit{dN}/d\gamma \propto {\gamma }^{-2}$for the overall particle population. Our results help to shed light on the ubiquitous presence of power-law particle and photon spectra in astrophysical nonthermal sources.

    

   more » « less
2. Microphysics of Relativistic Collisionless Electron-ion-positron Shocks

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

   Grošelj, Daniel ; Sironi, Lorenzo ; Beloborodov, Andrei M. ( July 2022 , The Astrophysical Journal)

   We perform particle-in-cell simulations to elucidate the microphysics of relativistic weakly magnetized shocks loaded with electron-positron pairs. Various external magnetizationsσ≲ 10⁻⁴and 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 as${\overline{E}}_{\mathrm{e}}\propto {\left(Z_{\pm }+1\right)}^{-1}$. (3) Pair loading suppresses nonthermal ion acceleration at magnetizations as low asσ≈ 5 × 10⁻⁶. The ions then become essentially thermal with mean energy${\overline{E}}_{\mathrm{i}}$, while electrons form a nonthermal tail, extending from$E\sim {\left(Z_{\pm }+1\right)}^{-1}{\overline{E}}_{\mathrm{i}}$to${\overline{E}}_{\mathrm{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
3. Tidal Disruption Event Demographics with the Zwicky Transient Facility: Volumetric Rates, Luminosity Function, and Implications for the Local Black Hole Mass Function

   https://doi.org/10.3847/2041-8213/acf216  

   Yao, Yuhan ; Ravi, Vikram ; Gezari, Suvi ; van Velzen, Sjoert ; Lu, Wenbin ; Schulze, Steve ; Somalwar, Jean J. ; Kulkarni, S. R. ; Hammerstein, Erica ; Nicholl, Matt ; et al ( September 2023 , The Astrophysical Journal Letters)

   We conduct a systematic tidal disruption event (TDE) demographics analysis using the largest sample of optically selected TDEs. A flux-limited, spectroscopically complete sample of 33 TDEs is constructed using the Zwicky Transient Facility over 3 yr (from 2018 October to 2021 September). We infer the black hole (BH) mass (M_BH) with host galaxy scaling relations, showing that the sampleM_BHranges from 10^5.1M_⊙to 10^8.2M_⊙. We developed a survey efficiency corrected maximum volume method to infer the rates. The rest-frameg-band luminosity function can be well described by a broken power law of$\varphi (L_g)\propto {\left[{\left(L_g/L_{\mathrm{bk}}\right)}^{0.3}+{\left(L_g/L_{\mathrm{bk}}\right)}^{2.6}\right]}^{-1}$, withL_bk= 10^43.1erg s⁻¹. In the BH mass regime of 10^5.3≲ (M_BH/M_⊙) ≲ 10^7.3, the TDE mass function follows$\varphi (M_{\mathrm{BH}})\propto M_{\mathrm{BH}}^{-0.25}$, which favors a flat local BH mass function (${\mathit{dn}}_{\mathrm{BH}}/d\log M_{\mathrm{BH}}\approx \mathrm{constant}$). We confirm the significant rate suppression at the high-mass end (M_BH≳ 10^7.5M_⊙), which is consistent with theoretical predictions considering direct capture of hydrogen-burning stars by the event horizon. At a host galaxy mass ofM_gal∼ 10¹⁰M_⊙, the average optical TDE rate is ≈3.2 × 10⁻⁵galaxy⁻¹yr⁻¹. We constrain the optical TDE rate to be [3.7, 7.4, and 1.6] × 10⁻⁵galaxy⁻¹yr⁻¹in galaxies with red, green, and blue colors.

    

   more » « less
4. Cooling Envelope Model for Tidal Disruption Events

   https://doi.org/10.3847/2041-8213/ac90ba  

   Metzger, Brian D. ( September 2022 , The Astrophysical Journal Letters)

   We present a toy model for the thermal optical/UV/X-ray emission from tidal disruption events (TDEs). Motivated by recent hydrodynamical simulations, we assume that the debris streams promptly and rapidly circularize (on the orbital period of the most tightly bound debris), generating a hot quasi-spherical pressure-supported envelope of radiusR_v∼ 10¹⁴cm (photosphere radius ∼10¹⁵cm) surrounding the supermassive black hole (SMBH). As the envelope cools radiatively, it undergoes Kelvin–Helmholtz contractionR_v∝t⁻¹, its temperature risingT_eff∝t^1/2while its total luminosity remains roughly constant; the optical luminosity decays as$\nu L_{\nu }\propto R_v^2{T}_{\mathrm{eff}}\propto t^{-3/2}$. Despite this similarity to the mass fallback rate${\dot{M}}_{\mathrm{fb}}\propto t^{-5/3}$, envelope heating from fallback accretion is subdominant compared to the envelope cooling luminosity except near optical peak (where they are comparable). Envelope contraction can be delayed by energy injection from accretion from the inner envelope onto the SMBH in a regulated manner, leading to a late-time flattening of the optical/X-ray light curves, similar to those observed in some TDEs. Eventually, as the envelope contracts to near the circularization radius, the SMBH accretion rate rises to its maximum, in tandem with the decreasing optical luminosity. This cooling-induced (rather than circularization-induced) delay of up to several hundred days may account for the delayed onset of thermal X-rays, late-time radio flares, and high-energy neutrino generation, observed in some TDEs. We compare the model predictions to recent TDE light-curve correlation studies, finding both agreement and points of tension.

    

   more » « less
5. The Pan-STARRS1 z > 5.6 Quasar Survey. III. The z ≈ 6 Quasar Luminosity Function

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

   Schindler, Jan-Torge ; Bañados, Eduardo ; Connor, Thomas ; Decarli, Roberto ; Fan, Xiaohui ; Farina, Emanuele Paolo ; Mazzucchelli, Chiara ; Nanni, Riccardo ; Rix, Hans-Walter ; Stern, Daniel ; et al ( January 2023 , The Astrophysical Journal)

   We present thez≈ 6 type-1 quasar luminosity function (QLF), based on the Pan-STARRS1 (PS1) quasar survey. The PS1 sample includes 125 quasars atz≈ 5.7–6.2, with −28 ≲M₁₄₅₀≲ −25. With the addition of 48 fainter quasars from the SHELLQs survey, we evaluate thez≈ 6 QLF over −28 ≲M₁₄₅₀≲ −22. Adopting a double power law with an exponential evolution of the quasar density (Φ(z) ∝ 10^k(z−6);k= −0.7), we use a maximum likelihood method to model our data. We find a break magnitude of$M^{*}=-{26.38}_{-0.60}^{+0.79}\mathrm{mag}$, a faint-end slope of$\alpha =-{1.70}_{-0.19}^{+0.29}$, and a steep bright-end slope of$\beta =-{3.84}_{-1.21}^{+0.63}$. Based on our new QLF model, we determine the quasar comoving spatial density atz≈ 6 to be$n(M_{1450}<-26)={1.16}_{-0.12}^{+0.13}{\mathrm{cGpc}}^{-3}$. In comparison with the literature, we find the quasar density to evolve with a constant value ofk≈ −0.7, fromz≈ 7 toz≈ 4. Additionally, we derive an ionizing emissivity of${ϵ}_{912}(z=6)={7.23}_{-1.02}^{+1.65}\times {10}^{22}\mathrm{erg}{\mathrm{s}}^{-1}{\mathrm{Hz}}^{-1}{\mathrm{cMpc}}^{-3}$, based on the QLF measurement. Given standard assumptions, and the recent measurement of the mean free path by Becker et al. atz≈ 6, we calculate an Hiphotoionizing rate of Γ_{H I}(z= 6) ≈ 6 × 10⁻¹⁶s⁻¹, strongly disfavoring a dominant role of quasars in hydrogen reionization.

    

   more » « less