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 present UV and/or optical observations and models of SN 2023ixf, a type II supernova (SN) located in Messier 101 at 6.9 Mpc. Early time (flash) spectroscopy of SN 2023ixf, obtained primarily at Lick Observatory, reveals emission lines of Hi, Hei/ii, Civ, and Niii/iv/vwith a narrow core and broad, symmetric wings arising from the photoionization of dense, close-in circumstellar material (CSM) located around the progenitor star prior to shock breakout. These electron-scattering broadened line profiles persist for ∼8 days with respect to first light, at which time Doppler broadened the features from the fastest SN ejecta form, suggesting a reduction in CSM density atr≳ 1015cm. The early time light curve of SN 2023ixf shows peak absolute magnitudes (e.g.,Mu= −18.6 mag,Mg= −18.4 mag) that are ≳2 mag brighter than typical type II SNe, this photometric boost also being consistent with the shock power supplied from CSM interaction. Comparison of SN 2023ixf to a grid of light-curve and multiepoch spectral models from the non-LTE radiative transfer codeCMFGENand the radiation-hydrodynamics codeHERACLESsuggests dense, solar-metallicity CSM confined tor= (0.5–1) × 1015cm, and a progenitor mass-loss rate ofṀ=102Myr−1. For the assumed progenitor wind velocity ofvw= 50 km s−1, this corresponds to enhanced mass loss (i.e.,superwindphase) during the last ∼3–6 yr before explosion.

     
    more » « less
  2. 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
  3. 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
  4. Abstract

    We compare the core-collapse evolution of a pair of 15.8Mstars with significantly different internal structures, a consequence of the bimodal variability exhibited by massive stars during their late evolutionary stages. The 15.78 and 15.79Mprogenitors have core masses (masses interior to an entropy of 4kBbaryon−1) of 1.47 and 1.78Mand compactness parametersξ1.75of 0.302 and 0.604, respectively. The core-collapse simulations are carried out in 2D to nearly 3 s postbounce and show substantial differences in the times of shock revival and explosion energies. The 15.78Mmodel begins exploding promptly at 120 ms postbounce when a strong density decrement at the Si–Si/O shell interface, not present in the 15.79Mprogenitor, encounters the stalled shock. The 15.79Mmodel takes 100 ms longer to explode but ultimately produces a more powerful explosion. Both the larger mass accretion rate and the more massive core of the 15.79Mmodel during the first 0.8 s postbounce time result in largerνe/ν¯eluminosities and RMS energies along with a flatter and higher-density heating region. The more-energetic explosion of the 15.79Mmodel resulted in the ejection of twice as much56Ni. Most of the ejecta in both models are moderately proton rich, though counterintuitively the highest electron fraction (Ye= 0.61) ejecta in either model are in the less-energetic 15.78Mmodel, while the lowest electron fraction (Ye= 0.45) ejecta in either model are in the 15.79Mmodel.

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