skip to main content


Title: SNART: Synchrotron Self Absorption Model For Non-relativistic Transients
Abstract

We presentSNART(DeMarchi & Finstad 2023), software for radio synchrotron self-absorption analysis, a generalized version of the model introduced in Chevalier, Chevalier & Fransson.SNART ingests flux density and frequency observations, fits them jointly, and outputs system physical parameters: magnetic fieldB, shock radiusR, post shock energyU, electron number densityne, circumstellar densityρCSM, and mass loss degenerate with wind velocityM˙/vwind.SNART is written in Python and is publicly available via GitHub. The repository hosts an in-depth derivation of the model and a detailed description of parameter definitions in the literature.SNART is a generalized treatment of synchrotron self absorption that leaves the choice of values forp(the power-law index of the electron distribution),q(shock front acceleration),θ(electron pitch angle, commonlyπ/2), andf(the “filling factor,” often 0.5) explicit.

 
more » « less
NSF-PAR ID:
10408747
Author(s) / Creator(s):
; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
Research Notes of the AAS
Volume:
7
Issue:
4
ISSN:
2515-5172
Format(s):
Medium: X Size: Article No. 77
Size(s):
Article No. 77
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    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
  2. 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
  3. Abstract

    We compare 500 pc scale, resolved observations of ionized and molecular gas for thez∼ 0.02 starbursting disk galaxy IRAS08339+6517, using measurements from KCWI and NOEMA. We explore the relationship of the star-formation-driven ionized gas outflows with colocated galaxy properties. We find a roughly linear relationship between the outflow mass flux (Σ̇out) and star formation rate surface density (ΣSFR),Σ̇outΣSFR1.06±0.10, and a strong correlation betweenΣ̇outand the gas depletion time, such thatΣ̇outtdep1.1±0.06. Moreover, we find these outflows are so-calledbreakoutoutflows, according to the relationship between the gas fraction and disk kinematics. Assuming that ionized outflow mass scales with total outflow mass, our observations suggest that the regions of highest ΣSFRin IRAS08 are removing more gas via the outflow than through the conversion of gas into stars. Our results are consistent with a picture in which the outflow limits the ability of a region of a disk to maintain short depletion times. Our results underline the need for resolved observations of outflows in more galaxies.

     
    more » « less
  4. Abstract

    State transitions in black hole X-ray binaries are likely caused by gas evaporation from a thin accretion disk into a hot corona. We present a height-integrated version of this process, which is suitable for analytical and numerical studies. With radiusrscaled to Schwarzschild units and coronal mass accretion rateṁcto Eddington units, the results of the model are independent of black hole mass. State transitions should thus be similar in X-ray binaries and an active galactic nucleus. The corona solution consists of two power-law segments separated at a break radiusrb∼ 103(α/0.3)−2, whereαis the viscosity parameter. Gas evaporates from the disk to the corona forr>rb, and condenses back forr<rb. Atrb,ṁcreaches its maximum,ṁc,max0.02(α/0.3)3. If atrrbthe thin disk accretes withṁd<ṁc,max, then the disk evaporates fully before reachingrb, giving the hard state. Otherwise, the disk survives at all radii, giving the thermal state. While the basic model considers only bremsstrahlung cooling and viscous heating, we also discuss a more realistic model that includes Compton cooling and direct coronal heating by energy transport from the disk. Solutions are again independent of black hole mass, andrbremains unchanged. This model predicts strong coronal winds forr>rb, and aT∼ 5 × 108K Compton-cooled corona forr<rb. Two-temperature effects are ignored, but may be important at small radii.

     
    more » « less
  5. 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