MOSFiT: Modular Open Source Fitter for Transients Guillochon, James; Nicholl, Matt; Villar, V. Ashley ...
The Astrophysical journal. Supplement series,
05/2018, Letnik:
236, Številka:
1
Journal Article
Recenzirano
Odprti dostop
Much of the progress made in time-domain astronomy is accomplished by relating observational multiwavelength time-series data to models derived from our understanding of physical laws. This goal is ...typically accomplished by dividing the task in two: collecting data (observing), and constructing models to represent that data (theorizing). Owing to the natural tendency for specialization, a disconnect can develop between the best available theories and the best available data, potentially delaying advances in our understanding new classes of transients. We introduce MOSFiT: the Modular Open Source Fitter for Transients, a Python-based package that downloads transient data sets from open online catalogs (e.g., the Open Supernova Catalog), generates Monte Carlo ensembles of semi-analytical light-curve fits to those data sets and their associated Bayesian parameter posteriors, and optionally delivers the fitting results back to those same catalogs to make them available to the rest of the community. MOSFiT is designed to help bridge the gap between observations and theory in time-domain astronomy; in addition to making the application of existing models and creation of new models as simple as possible, MOSFiT yields statistically robust predictions for transient characteristics, with a standard output format that includes all the setup information necessary to reproduce a given result. As large-scale surveys such as that conducted with the Large Synoptic Survey Telescope (LSST), discover entirely new classes of transients, tools such as MOSFiT will be critical for enabling rapid comparison of models against data in statistically consistent, reproducible, and scientifically beneficial ways.
We present detailed multiwavelength observations of GRB 161219B at z = 0.1475, spanning the radio to X-ray regimes, and the first Atacama Large Millimeter/submillimeter Array (ALMA) light curve of a ...γ-ray burst (GRB) afterglow. The centimeter- and millimeter-band observations before 8.5 days require emission in excess of that produced by the afterglow forward shock (FS). These data are consistent with radiation from a refreshed reverse shock (RS) produced by the injection of energy into the FS, signatures of which are also present in the X-ray and optical light curves. We infer a constant-density circumburst environment with an extremely low density, , and show that this is a characteristic of all strong RS detections to date. The Karl G. Lansky Very Large Array (VLA) observations exhibit unexpected rapid variability on roughly minute timescales, indicative of strong interstellar scintillation. The X-ray, ALMA, and VLA observations together constrain the jet break time, days, yielding a wide jet opening angle of , implying beaming-corrected γ-ray and kinetic energies of erg and erg, respectively. Comparing the RS and FS emission, we show that the ejecta are only weakly magnetized, with relative magnetization, , compared to the FS. These direct, multifrequency measurements of a refreshed RS spanning the optical to radio bands highlight the impact of radio and millimeter data in probing the production and nature of GRB jets.
Abstract
Massive, rapidly spinning magnetar remnants produced as a result of binary neutron-star (BNS) mergers may deposit a fraction of their energy into the surrounding kilonova ejecta, powering a ...synchrotron radio signal from the interaction of the ejecta with the circumburst medium. Here, we present 6.0 GHz Very Large Array (VLA) observations of nine, low-redshift short gamma-ray bursts (GRBs;
z
< 0.5) on rest-frame timescales of ≈2.4–13.9 yr following the bursts. We place 3
σ
limits on radio continuum emission of
F
ν
≲ 6–20
μ
Jy at the burst positions, or
L
ν
≲ (0.6–8.3) × 10
28
erg s
−1
Hz
−1
. Comparing these limits with new light-curve modeling that properly incorporates relativistic effects, we obtain limits on the energy deposited into the ejecta of
E
ej
≲ (0.6–6.7) × 10
52
erg (
erg) for an ejecta mass of 0.03
M
⊙
(0.1
M
⊙
). We present a uniform reanalysis of 27 short GRBs with 5.5–6.0 GHz observations, and find that ≳50% of short GRBs did not form stable magnetar remnants in their mergers. Assuming short GRBs are produced by BNS mergers drawn from the Galactic BNS population plus an additional component of high-mass GW194025-like mergers in a fraction
f
GW190425
of cases, we place constraints on the maximum mass of a nonrotating neutron star (NS; Tolman–Oppenheimer–Volkoff mass;
M
TOV
), finding
for
f
GW190425
= 0.4; this limit increases for larger values of
f
GW190425
. The detection (or lack thereof) of radio remnants in untargeted surveys such as the VLA Sky Survey (VLASS) could provide more stringent constraints on the fraction of mergers that produce stable remnants. If
radio remnants are discovered in VLASS, this suggests that short GRBs are a biased population of BNS mergers in terms of the stability of the remnants they produce.
A REVERSE SHOCK IN GRB 160509A Laskar, Tanmoy; Alexander, Kate D.; Berger, Edo ...
The Astrophysical journal,
12/2016, Letnik:
833, Številka:
1
Journal Article
Recenzirano
Odprti dostop
ABSTRACT We present the second multi-frequency radio detection of a reverse shock in a γ-ray burst. By combining our extensive radio observations of the Fermi-Large Area Telescope γ-ray burst 160509A ...at z = 1.17 up to 20 days after the burst with Swift X-ray observations and ground-based optical and near-infrared data, we show that the afterglow emission comprises distinct reverse shock and forward shock contributions: the reverse shock emission dominates in the radio band at 10 days, while the forward shock emission dominates in the X-ray, optical, and near-infrared bands. Through multi-wavelength modeling, we determine a circumburst density of , supporting our previous suggestion that a low-density circumburst environment is conducive to the production of long-lasting reverse shock radiation in the radio band. We infer the presence of a large excess X-ray absorption column, NH 1.5 × 1022 , and a high rest-frame optical extinction, AV 3.4 mag. We identify a jet break in the X-ray light curve at , and thus derive a jet opening angle of , yielding a beaming-corrected kinetic energy and radiated γ-ray energy of erg and erg (1-104 keV, rest frame), respectively. Consistency arguments connecting the forward shocks and reverse shocks suggest a deceleration time of s T90, a Lorentz factor of , and a reverse-shock-to-forward-shock fractional magnetic energy density ratio of . Our study highlights the power of rapid-response radio observations in the study of the properties and dynamics of γ-ray burst ejecta.
Abstract
We present 1.3 mm (230 GHz) observations of the recent and nearby Type II supernova, SN 2023ixf, obtained with the Submillimeter Array (SMA) at 2.6–18.6 days after explosion. The ...observations were obtained as part the SMA Large Program, POETS (Pursuit of Extragalactic Transients with the SMA). We do not detect any emission at the location of SN 2023ixf, with the deepest limits of
L
ν
(230 GHz) ≲ 8.6 × 10
25
erg s
−1
Hz
−1
at 2.7 and 7.7 days, and
L
ν
(230 GHz) ≲ 3.4 × 10
25
erg s
−1
Hz
−1
at 18.6 days. These limits are about a factor of 2 times dimmer than the millimeter emission from SN 2011dh (IIb), about 1 order of magnitude dimmer compared to SN 1993J (IIb) and SN 2018ivc (IIL), and about 30 times dimmer than the most luminous nonrelativistic SNe in the millimeter band (Type IIb/Ib/Ic). Using these limits in the context of analytical models that include synchrotron self-absorption and free–free absorption, we place constraints on the proximate circumstellar medium around the progenitor star, to a scale of ∼2 × 10
15
cm, excluding the range
M
̇
∼
few
×
10
−
6
−
10
−
2
M
⊙
yr
−1
(for a wind velocity,
v
w
= 115 km s
−1
, and ejecta velocity,
v
ej
∼ (1 − 2) × 10
4
km s
−1
). These results are consistent with an inference of the mass-loss rate based on optical spectroscopy (∼2 × 10
−2
M
⊙
yr
−1
for
v
w
= 115 km s
−1
), but are in tension with the inference from hard X-rays (∼7 × 10
−4
M
⊙
yr
−1
for
v
w
= 115 km s
−1
). This tension may be alleviated by a nonhomogeneous and confined CSM, consistent with results from high-resolution optical spectroscopy.
We present detailed multifrequency, multiepoch radio observations of GRB 140304A at z = 5.283 from 1 to 86 GHz and from 0.45 to 89 days. The radio and millimeter data exhibit unusual multiple ...spectral components, which cannot be simply explained by standard forward and reverse shock scenarios. Through detailed multiwavelength analysis spanning radio to X-rays, we constrain the forward shock parameters to Ek,iso 4.9 × 1054 erg, A * 2.6 × 10−2, ϵ e 2.5 × 10−2, ϵ B 5.9 × 10−2, p 2.6, and θ jet 1 1, yielding a beaming-corrected γ-ray and kinetic energy, E γ 2.3 × 1049 erg and E K 9.5 × 1050 erg, respectively. We model the excess radio emission as due to a combination of a late-time reverse shock (RS) launched by a shell collision, which also produces a rebrightening in the X-rays at 0.26 days, and either a standard RS or diffractive interstellar scintillation (ISS). Under the standard RS interpretation, we invoke consistency arguments between the forward and reverse shocks to derive a deceleration time, tdec 100 s, the ejecta Lorentz factor, Γ(tdec) 300, and a low RS magnetization, RB 0.6. Our observations highlight both the power of radio observations in capturing RS emission and thus constraining the properties of GRB ejecta and central engines and the challenge presented by ISS in conclusively identifying RS emission in GRB radio afterglows.
We report the white light flare rates for 10 ultracool dwarfs using Kepler K2 short-cadence data. Among our sample stars, two have spectral type M6, three are M7, three are M8, and two are L0. Most ...of our targets are old low-mass stars. We identify a total of 283 flares in all of the stars in our sample, with Kepler energies in the range log EKp ∼ (29-33.5) erg. Using the maximum-likelihood method of line fitting, we find that the flare frequency distribution (FFD) for each star in our sample follows a power law with slope − in the range −(1.3-2.0). We find that cooler objects tend to have shallower slopes. For some of our targets, the FFD follows either a broken power law, or a power law with an exponential cutoff. For the L0 dwarf 2MASS J12321827-0951502, we find a very shallow slope (− = −1.3) in the Kepler energy range (0.82-130) × 1030 erg: this L0 dwarf has flare rates which are comparable to those of high-energy flares in stars of earlier spectral types. In addition, we report photometry of two superflares: one on the L0 dwarf 2MASS J12321827-0951502 and another on the M7 dwarf 2MASS J08352366+1029318. In the case of 2MASS J12321827-0951502, we report a flare brightening by a factor of ∼144 relative to the quiescent photospheric level. Likewise, for 2MASS J08352366+1029318, we report a flare brightening by a factor of ∼60 relative to the quiescent photospheric level. These two superflares have bolometric (ultraviolet/optical/infrared) energies 3.6 × 1033 erg and 8.9 × 1033 erg respectively, while the full width half maximum timescales are very short, ∼2 min. We find that the M8 star TRAPPIST-1 is more active than the M8.5 dwarf 2M03264453+1919309, but less active than another M8 dwarf (2M12215066-0843197).
A measurement of the wind speed on a brown dwarf Allers, Katelyn N; Vos, Johanna M; Biller, Beth A ...
Science (American Association for the Advancement of Science),
04/2020, Letnik:
368, Številka:
6487
Journal Article
Recenzirano
Odprti dostop
Zonal (latitudinal) winds dominate the bulk flow of planetary atmospheres. For gas giant planets such as Jupiter, the motion of clouds can be compared with radio emissions from the magnetosphere, ...which is connected to the planet's interior, to determine the wind speed. In principle, this technique can be applied to brown dwarfs and/or directly imaged exoplanets if periods can be determined for both the infrared and radio emissions. We apply this method to measure the wind speeds on the brown dwarf 2MASS J10475385+2124234. The difference between the radio period of 1.751 to 1.765 hours and infrared period of 1.741 ± 0.007 hours implies a strong wind (+650 ± 310 meters per second) proceeding eastward. This could be due to atmospheric jet streams and/or low frictional drag at the bottom of the atmosphere.
Abstract
We present a detailed compilation and analysis of the X-ray phase space of low- to intermediate-redshift (0 ≤
z
≤ 1) transients that consolidates observed light curves (and theory where ...necessary) for a large variety of classes of transient/variable phenomena in the 0.3–10 keV energy band. We include gamma-ray burst afterglows, supernovae, supernova shock breakouts and shocks interacting with the environment, tidal disruption events and active galactic nuclei, fast blue optical transients, cataclysmic variables, magnetar flares/outbursts and fast radio bursts, cool stellar flares, X-ray binary outbursts, and ultraluminous X-ray sources. Our overarching goal is to offer a comprehensive resource for the examination of these ephemeral events, extending the X-ray duration–luminosity phase space (DLPS) to show luminosity evolution. We use existing observations (both targeted and serendipitous) to characterize the behavior of various transient/variable populations. Contextualizing transient signals in the larger DLPS serves two primary purposes: to identify areas of interest (i.e., regions in the parameter space where one would expect detections, but in which observations have historically been lacking), and to provide initial qualitative guidance in classifying newly discovered transient signals. We find that while the most luminous (largely extragalactic) and least luminous (largely Galactic) part of the phase space is well populated at
t
> 0.1 days, intermediate-luminosity phenomena (
L
X
= 10
34
–10
42
erg s
−1
) represent a gap in the phase space. We thus identify
L
X
= 10
34
–10
42
erg s
−1
and
t
= 10
−4
to 0.1 days as a key discovery phase space in transient X-ray astronomy.