We present constraints on Ks-band emission from one of the nearest short hard gamma-ray bursts, GRB 160821B, at z = 0.16, at three epochs. We detect a red relativistic afterglow from the jetted ...emission in the first epoch but do not detect any excess kilonova emission in the second two epochs. We compare upper limits obtained with Keck I/MOSFIRE to multi-dimensional radiative transfer models of kilonovae, that employ composition-dependent nuclear heating and LTE opacities of heavy elements. We discuss eight models that combine toroidal dynamical ejecta and two types of wind and one model with dynamical ejecta only. We also discuss simple, empirical scaling laws of predicted emission as a function of ejecta mass and ejecta velocity. Our limits for GRB 160821B constrain the ejecta mass to be lower than 0.03 M for velocities greater than 0.1 c. At the distance sensitivity range of advanced LIGO, similar ground-based observations would be sufficiently sensitive to the full range of predicted model emission including models with only dynamical ejecta. The color evolution of these models shows that I-K color spans 7-16 mag, which suggests that even relatively shallow infrared searches for kilonovae could be as constraining as optical searches.
Abstract
We report our Spitzer Space Telescope observations and detections of the binary neutron star merger GW170817. At 4.5 μm, GW170817 is detected at 21.9 mag AB at +43 days and 23.9 mag AB at ...+74 days after merger. At 3.6 μm, GW170817 is not detected to a limit of 23.2 mag AB at +43 days and 23.1 mag AB at +74 days. Our detections constitute the latest and reddest constraints on the kilonova/macronova emission and composition of heavy elements. The 4.5 μm luminosity at this late phase cannot be explained by elements exclusively from the first abundance peak of the r-process. Moreover, the steep decline in the Spitzer band, with a power-law index of 3.4 ± 0.2, can be explained by a few of the heaviest isotopes with half-life around 14 d dominating the luminosity (e.g. 140Ba, 143Pr, 147Nd, 156Eu, 191Os, 223Ra, 225Ra, 233Pa, 234Th) or a model with lower deposition efficiency. This data offers evidence that the heaviest elements in the second and third r-process abundance peak were indeed synthesized. Our conclusion is verified by both analytics and network simulations and robust despite intricacies and uncertainties in the nuclear physics. Future observations with Spitzer and James Webb Space Telescope will further illuminate the relative abundance of the synthesized heavy elements.
We present a dust spectral energy distribution (SED) and binary stellar population analysis revisiting the dust production rates (DPRs) in the winds of carbon-rich Wolf-Rayet (WC) binaries and their ...impact on galactic dust budgets. DustEM SED models of 19 Galactic WC "dustars" reveal DPRs of M yr−1 and carbon dust condensation fractions, χC, between 0.002% and 40%. A large (0.1-1.0 m) dust grain size composition is favored for efficient dustars where χC 1%. Results for dustars with known orbital periods verify a power-law relation between χC, orbital period, WC mass-loss rate, and wind velocity consistent with predictions from theoretical models of dust formation in colliding-wind binaries. We incorporated dust production into Binary Population and Spectral Synthesis (BPASS) models to analyze dust production rates from WC dustars, asymptotic giant branch stars (AGBs), red supergiants (RSGs), and core-collapse supernovae (SNe). BPASS models assuming constant star formation (SF) and a coeval 106 M stellar population were performed at low, Large Magellanic Cloud (LMC)-like, and solar metallicities (Z = 0.001, 0.008, and 0.020). Both constant SF and coeval models indicate that SNe are net dust destroyers at all metallicities. Constant SF models at LMC-like metallicities show that AGB stars slightly outproduce WC binaries and RSGs by factors of 2-3, whereas at solar metallicities WC binaries are the dominant source of dust for ∼60 Myr until the onset of AGBs, which match the dust input of WC binaries. Coeval population models show that, for "bursty" SF, AGB stars dominate dust production at late times (t 70 Myr).
Abstract
We present Spitzer/InfraRed Array Camera observations of dust formation from six extragalactic carbon-rich Wolf-Rayet (WC) binary candidates in low-metallicity (
Z
≲ 0.65
Z
⊙
) environments ...using multiepoch mid-infrared (IR) imaging data from the SPitzer InfraRed Intensive Transients Survey (SPIRITS). Optical follow-up spectroscopy of SPIRITS 16ln, 19q, 16df, 18hb, and 14apu reveals emission features from C
iv
λ
5801−12 and/or the C
iii–iv
λ
4650/He
ii
λ
4686 blend that are consistent with early-type WC stars. We identify SPIRITS 16ln as the variable mid-IR counterpart of the recently discovered colliding-wind WC4 + O binary candidate, N604-WRXc, located in the subsolar metallicity NGC 604 H
ii
region in M33. We interpret the mid-IR variability from SPIRITS 16ln as a dust-formation episode in an eccentric colliding-wind WC binary. SPIRITS 19q, 16df, 14apu, and 18hb exhibit absolute 3.6 magnitudes exceeding that of one of the most IR-luminous dust-forming WC systems known, WR 104 (
M
3.6
≲ −12.3). An analysis of dust formation in the mid-IR outburst from SPIRITS 19q reveals a high dust production rate of
M
⊙
yr
−1
, which may therefore exceed that of the most efficient dust-forming WC systems known. We demonstrate that efficient dust formation is feasible from early-type WC binaries in the theoretical framework of colliding-wind binary dust formation if the systems host an O-type companion with high mass-loss rates (
M
⊙
yr
−1
). This efficient dust formation from early-type WC binaries highlights their potential role as significant sources of dust in low-metallicity environments.
Abstract
The Wolf–Rayet (WR) binary system WR140 is a close (0.9–16.7 mas; ref.
1
) binary star consisting of an O5 primary and WC7 companion
2
and is known as the archetype of episodic ...dust-producing WRs. Dust in WR binaries is known to form in a confined stream originating from the collision of the two stellar winds, with orbital motion of the binary sculpting the large-scale dust structure into arcs as dust is swept radially outwards. It is understood that sensitive conditions required for dust production in WR140 are only met around periastron when the two stars are sufficiently close
2–4
. Here we present multiepoch imagery of the circumstellar dust shell of WR140. We constructed geometric models that closely trace the expansion of the intricately structured dust plume, showing that complex effects induced by orbital modulation may result in a ‘Goldilocks zone’ for dust production. We find that the expansion of the dust plume cannot be reproduced under the assumption of a simple uniform-speed outflow, finding instead the dust to be accelerating. This constitutes a direct kinematic record of dust motion under acceleration by radiation pressure and further highlights the complexity of the physical conditions in colliding-wind binaries.
Abstract
We analyze pre-explosion near- and mid-infrared (IR) imaging of the site of SN 2023ixf in the nearby spiral galaxy M101 and characterize the candidate progenitor star. The star displays ...compelling evidence of variability with a possible period of ≈1000 days and an amplitude of Δ
m
≈ 0.6 mag in extensive monitoring with the Spitzer Space Telescope since 2004, likely indicative of radial pulsations. Variability consistent with this period is also seen in the near-IR
J
and
K
s
bands between 2010 and 2023, up to just 10 days before the explosion. Beyond the periodic variability, we do not find evidence for any IR-bright pre-supernova outbursts in this time period. The IR brightness (
M
K
s
=
−
10.7
mag) and color (
J
−
K
s
= 1.6 mag) of the star suggest a luminous and dusty red supergiant. Modeling of the phase-averaged spectral energy distribution (SED) yields constraints on the stellar temperature (
T
eff
=
3500
−
1400
+
800
K) and luminosity (
log
L
/
L
⊙
=
5.1
±
0.2
). This places the candidate among the most luminous Type II supernova progenitors with direct imaging constraints, with the caveat that many of these rely only on optical measurements. Comparison with stellar evolution models gives an initial mass of
M
init
= 17 ± 4
M
⊙
. We estimate the pre-supernova mass-loss rate of the star between 3 and 19 yr before explosion from the SED modeling at
M
̇
≈
3
×
10
−
5
to 3 × 10
−4
M
⊙
yr
−1
for an assumed wind velocity of
v
w
= 10 km s
−1
, perhaps pointing to enhanced mass loss in a pulsation-driven wind.
We present a mid-infrared (IR) sample study of nearby ultraluminous X-ray sources (ULXs) using multiepoch observations with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. ...Spitzer/IRAC observations taken after 2014 were obtained as part of the Spitzer Infrared Intensive Transients Survey. Our sample includes 96 ULXs located within 10 Mpc. Of the 96 ULXs, 12 have candidate counterparts consistent with absolute mid-IR magnitudes of supergiants, and 16 counterparts exceeded the mid-IR brightness of single supergiants and are thus more consistent with star clusters or non-ULX background active galactic nuclei. The supergiant candidate counterparts exhibit a bimodal color distribution in a Spitzer/IRAC color-magnitude diagram, where "red" and "'blue" ULXs fall in IRAC colors 3.6 - 4.5 ∼ 0.7 and 3.6 - 4.5 ∼ 0.0, respectively. The mid-IR colors and absolute magnitudes of four "red" and five "blue" ULXs are consistent with those of supergiant Be (sgBe) and red supergiant (RSG) stars, respectively. Although "blue," RSG-like mid-IR ULX counterparts likely host RSG mass donors; we propose that "red" counterparts are ULXs exhibiting the "Be phenomenon" rather than hosts of sgBe mass donors. We show that the mid-IR excess from the "red" ULXs is likely due to thermal emission from circumstellar or circumbinary dust. Using dust as a probe for total mass, we estimate mass-loss rates of M ˙ ∼ 1 × 10 − 4 M yr−1 in dust-forming outflows of red ULXs. Based on the transient mid-IR behavior and its relatively flat spectral index, = −0.19 0.1, we suggest that the mid-IR emission from Holmberg IX X-1 originates from a variable jet.
Palomar Gattini-IR is a new wide-field, near-infrared (NIR) robotic time domain survey operating at Palomar Observatory. Using a 30 cm telescope mounted with a H2RG detector, Gattini-IR achieves a ...field of view (FOV) of 25 sq. deg. with a pixel scale of 8 7 in J-band. Here, we describe the system design, survey operations, data processing system and on-sky performance of Palomar Gattini-IR. As a part of the nominal survey, Gattini-IR scans 7500 square degrees of the sky every night to a median 5 depth of 15.7 AB mag outside the Galactic plane. The survey covers 15,000 square degrees of the sky visible from Palomar with a median cadence of 2 days. A real-time data processing system produces stacked science images from dithered raw images taken on sky, together with point-spread function (PSF)-fit source catalogs and transient candidates identified from subtractions within a median delay of 4 hr from the time of observation. The calibrated data products achieve an astrometric accuracy (rms) of 0 7 with respect to Gaia DR2 for sources with signal-to-noise ratio > 10, and better than 0 35 for sources brighter than 12 Vega mag. The photometric accuracy (rms) achieved in the PSF-fit source catalogs is better than 3% for sources brighter than 12 Vega mag and fainter than the saturation magnitude of 8.5 Vega mag, as calibrated against the Two Micron All Sky Survey catalog. The detection efficiency of transient candidates injected into the images is better than 90% for sources brighter than the 5 limiting magnitude. The photometric recovery precision of injected sources is 3% for sources brighter than 13 mag, and the astrometric recovery rms is 0 9. Reference images generated by stacking several field visits achieve depths of 16.5 AB mag over 60% of the sky, while it is limited by confusion in the Galactic plane. With a FOV 40× larger than any other existing NIR imaging instrument, Gattini-IR is probing the reddest and dustiest transients in the local universe such as dust obscured supernovae in nearby galaxies, novae behind large columns of extinction within the galaxy, reddened microlensing events in the Galactic plane and variability from cool and dust obscured stars. We present results from transients and variables identified since the start of the commissioning period.
Luminous red novae (LRNe) are astrophysical transients associated with the partial ejection of a binary system’s common envelope shortly before its merger. Here we present the results of our ...photometric and spectroscopic follow-up campaign of AT 2018bwo (DLT 18x), a LRN discovered in NGC 45, and investigate its progenitor system using binary stellar-evolution models. The transient reached a peak magnitude of
M
r
= −10.97 ± 0.11 and maintained this brightness during its optical plateau of
t
p
= 41 ± 5 days. During this phase, it showed a rather stable photospheric temperature of ∼3300 K and a luminosity of ∼10
40
erg s
−1
. Although the luminosity and duration of AT 2018bwo is comparable to the LRNe V838 Mon and M31-2015LRN, its photosphere at early times appears larger and cooler, likely due to an extended mass-loss episode before the merger. Toward the end of the plateau, optical spectra showed a reddened continuum with strong molecular absorption bands. The IR spectrum at +103 days after discovery was comparable to that of a M8.5 II type star, analogous to an extended AGB star. The reprocessed emission by the cooling dust was also detected in the mid-infrared bands ∼1.5 years after the outburst. Archival
Spitzer
and
Hubble
Space Telescope data taken 10−14 yrs before the transient event suggest a progenitor star with
T
prog
∼ 6500 K,
R
prog
∼ 100
R
⊙
, and
L
prog
= 2 × 10
4
L
⊙
, and an upper limit for optically thin warm (1000 K) dust mass of
M
d
< 10
−6
M
⊙
. Using stellar binary-evolution models, we determined the properties of binary systems consistent with the progenitor parameter space. For AT 2018bwo, we infer a primary mass of 12–16
M
⊙
, which is 9–45% larger than the ∼11
M
⊙
obtained using single-star evolution models. The system, consistent with a yellow-supergiant primary, was likely in a stable mass-transfer regime with −2.4 ≤ log(
Ṁ
/
M
⊙
yr
−1
) ≤ −1.2 a decade before the main instability occurred. During the dynamical merger, the system would have ejected 0.15–0.5
M
⊙
with a velocity of ∼500 km s
−1
.
Abstract
WR 137 (HD 192641) is a binary system consisting of a carbon-rich Wolf–Rayet (W-R) star and an Oe companion star in a 13 yr orbit. Near periastron, the winds of the two stars collide and ...form carbonaceous dust. We obtained three mid-infrared grism spectra of the system with SOFIA and FORCAST during the last year of SOFIA’s operations in 2021 July, 2021 February, and 2022 May (Cycle 9). Within these spectra, we have identified several wind lines from He
i
, He
ii
, C
iii
, and C
iv
that are emitted from the W-R wind as well as a weak emission feature around 6.3–6.4
μ
m that may have shifted its peak flux from 6.29 to 6.41
μ
m through this time period. The weak feature grew as the continuum dust emission grew while the W-R emission appeared to decline due to lower contrast with the continuum. Furthermore, we observe that the peak of the feature shifts to redder wavelengths during the observations. We compare this feature to the unidentified infrared feature and other emission lines identified in dusty carbon-rich W-R (WC) binaries. For WR 137, we speculate that mixing of the winds in the system with the Oe star’s disk is important for starting the dust formation and that it is less important as dust formation continues. Previous infrared photometry shows “minieruptions” of dust production, which could then be explained with variations of the Oe star disk.