Search

Classical Wolf-Rayet stars are the descendants of massive OB stars that have lost their hydrogen envelopes and are burning helium in their cores prior to exploding as type Ib/c supernovae. The mechanisms for losing their hydrogen envelopes are either through binary interactions or through strong stellar winds potentially coupled with episodic mass-loss. Amongst the bright classical WR stars, the binary system WR\,137 (HD\,192641; WC7d + O9e) is the subject of this paper. This binary is known to have a 13-year period and produces dust near periastron. Here we report on interferometry with the CHARA Array collected over a decade of time and providing the first visual orbit for the system. We combine these astrometric measurements with archival radial velocities to measure masses of the stars of $M_{\rm WR} = 9.5\pm3.4 M_\odot$ and $M_{\rm O} = 17.3\pm 1.9 M_\odot$ when we use the most recent \textit{Gaia} distance. These results are then compared to predicted dust distribution using these orbital elements, which match the observed imaging from \textit{JWST} as discussed recently by Lau et al. Furthermore, we compare the system to the BPASS models, finding that the WR star likely formed through stellar winds and not through binary interactions. However, the companion O star did likely accrete some material from the WR's mass-loss to provide the rotation seen today that drives its status as an Oe star., Comment: Accepted to ApJ

Several long-period binaries with a carbon-rich Wolf-Rayet star and an O star produce dust in their wind collisions. In eccentric binaries, this is seen most strongly near periastron passage. The exact conditions leading to dust creation require orbital properties to be determined, which is difficult owing to their long periods. Recently, the binary system WR 125 (WC7+O9III) began a dust creation episode seen through an infrared outburst first detected by NEOWISE-R, which was the first outburst detected since 1991. We present new near- and mid-infrared photometry, which we use to show consistency between the two outbursts and derive an orbital period of 28.12$^{+0.10}_{-0.05}$ yr. We use a long time-series of optical spectra to place the first constraints on its orbital elements, on the assumption that this system will produce dust near periastron. The orbit has a mild eccentricity of 0.29$\pm$0.12 and is only derived for the Wolf-Rayet component, as the O star's radial velocities have noise that is likely larger than the expected semi-amplitude of the orbit. We also present SOFIA/FORCAST grism spectroscopy to examine the infrared spectral energy distribution (SED) of the dust during this outburst, comparing its properties to other WCd binaries, deriving a dust temperature of 580 K in 2021. This collection of observations will allow us to plan future observations of this system and place the system in the context of dust-creating Wolf-Rayet binaries., Comment: accepted to ApJ

We present infrared aperture masking interferometry (AMI) observations of newly formed dust from the colliding winds of the massive binary system Wolf-Rayet (WR) 137 with JWST using the Near Infrared Imager and Slitless Spectrograph (NIRISS). NIRISS AMI observations of WR 137 and a point-spread-function calibrator star, HD~228337, were taken using the F380M and F480M filters in 2022 July and August as part of the Director's Discretionary Early Release Science (DD-ERS) program 1349. Interferometric observables (squared visibilities and closure phases) from the WR 137 "interferogram" were extracted and calibrated using three independent software tools: ImPlaneIA, AMICAL, and SAMpip. The analysis of the calibrated observables yielded consistent values except for slightly discrepant closure phases measured by ImPlaneIA. Based on all three sets of calibrated observables, images were reconstructed using three independent software tools: BSMEM, IRBis, and SQUEEZE. All reconstructed image combinations generated consistent images in both F380M and F480M filters. The reconstructed images of WR 137 reveal a bright central core with a $\sim300$ mas linear filament extending to the northwest. A geometric colliding-wind model with dust production constrained to the orbital plane of the binary system and enhanced as the system approaches periapsis provided a general agreement with the interferometric observables and reconstructed images. Based on a colliding-wind dust condensation analysis, we suggest that dust formation within the orbital plane of WR 137 is induced by enhanced equatorial mass-loss from the rapidly rotating O9 companion star, whose axis of rotation is aligned with that of the orbit., Comment: 18 pages, 8 figures, Accepted for publication in ApJ. Updated plotting error in Fig. 2

WR 137 (HD 192641) is a binary system consisting of a carbon-rich Wolf-Rayet star and an Oe companion star in a 13-year 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 July 2021, February 2021, and May 2022 (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 Wolf-Rayet wind as well as a weak emission feature around 6.3-6.4 $\mu$m that may have shifted its peak flux from 6.29 to 6.41$\mu$m through this time period. The weak feature grew as the continuum dust emission grew while the WR 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 UIR feature and other emission lines identified in dusty 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 "mini-eruptions" of dust production which could then be explained with variations of the Oe star disk., Comment: 13 pages, accepted to ApJ

Wolf-Rayet (WR) 140 is the archetypal periodic dust-forming colliding-wind binary that hosts a carbon-rich WR (WC) star and an O-star companion with an orbital period of 7.93 years and an orbital eccentricity of 0.9. Throughout the past several decades, multiple dust-formation episodes from WR 140 have been observed that are linked to the binary orbit and occur near the time of periastron passage. Given its predictable dust-formation episodes, WR 140 presents an ideal astrophysical laboratory for investigating the formation and evolution of dust in the hostile environment around a massive binary system. In this paper, we present near- and mid-infrared (IR) spectroscopic and imaging observations of WR 140 with Subaru/SCExAO+CHARIS, Keck/NIRC2+PyWFS, and Subaru/COMICS taken between 2020 June and Sept that resolve the circumstellar dust emission linked to its most recent dust-formation episode in 2016 Dec. Our spectral energy distribution (SED) analysis of WR 140's resolved circumstellar dust emission reveals the presence of a hot ($T_\mathrm{d}\sim1000$ K) near-IR dust component that is co-spatial with the previously known and cooler ($T_\mathrm{d}\sim500$ K) mid-IR dust component composed of $300-500$ {\AA}-sized dust grains. We attribute the hot near-IR dust emission to the presence of nano-sized ("nanodust") grains and suggest they were formed from grain-grain collisions or the rotational disruption of the larger grain size population by radiative torques in the strong radiation field from the central binary. Lastly, we speculate on the astrophysical implications of nanodust formation around colliding-wind WC binaries, which may present an early source of carbonaceous nanodust in the interstellar medium., Comment: 21 pages, 8 figures, Accepted for publication in ApJ

Massive colliding-wind binaries that host a Wolf-Rayet (WR) star present a potentially important source of dust and chemical enrichment in the interstellar medium (ISM). However, the chemical composition and survival of dust formed from such systems is not well understood. The carbon-rich WR (WC) binary WR~140 presents an ideal astrophysical laboratory for investigating these questions given its well-defined orbital period and predictable dust-formation episodes every 7.93 years around periastron passage. We present observations from our Early Release Science program (ERS1349) with the James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) Medium-Resolution Spectrometer (MRS) and Imager that reveal the spectral and spatial signatures of nested circumstellar dust shells around WR~140. MIRI MRS spectroscopy of the second dust shell and Imager detections of over 17 shells formed throughout the past $\gtrsim130$ years confirm the survival of carbonaceous dust grains from WR~140 that are likely carriers of "unidentified infrared" (UIR)-band features at 6.4 and 7.7 $\mu$m. The observations indicate that dust-forming WC binaries can enrich the ISM with organic compounds and carbonaceous dust., Comment: Published in Nature Astronomy on Oct 12, 2022; 21 pages, 5 figures, 2 tables

We present the detection of a broad 8 $\mu$m feature in newly formed dust around the carbon-rich Wolf-Rayet (WC) binary WR 125 from N-band low-resolution (NL; R$\sim$250) spectroscopy between 7.3-13.6 $\mu$m and N-band (11.7 $\mu$m) and Q-band (18.8 $\mu$m) imaging with Subaru/COMICS in 2019 October. WR 125 is a colliding wind binary (${\rm WC7+O9}$) that exhibited renewed dust formation starting in 2018, $\sim$28 years after its first dust formation episode had been observed. We also compare our infrared photometry with historical observations and revise the dust-formation period of WR 125 to 28.1 years. Archival infrared spectra of five dusty WC stars, WR 48a, WR 98a, WR 104, WR 112 and WR 118, obtained with ISO/SWS are reanalyzed and compared with the WR 125 spectrum to search for a similar feature. We analyze the dusty WC spectra using two different extinction curves to investigate the impact of interstellar extinction correction on the presence and/or properties of the 8 $\mu$m feature. All of the dusty WC spectra dereddened with the two different extinction curves show a broad feature around 8 $\mu$m (FWHM$\sim$1-2 $\mu$m). We suggest that these 8 $\mu$m features seen in the dusty WC spectra are related to the Class C unidentified infrared (UIR) features., Comment: 17 pages, 7 figures, 4 tables, accepted to ApJ

Mid-infrared photometry of the Wolf-Rayet star HD 38030 in the Large Magellanic Cloud from the NEOWISE-R mission show it to have undergone a dust-formation episode in 2018 and the dust to have cooled in 2019-20. New spectroscopy with the MagE spectrograph on the Magellan I Baade Telescope in 2019 and 2020 show absorption lines attributable to a companion of type near O9.7III-IV. We found a significant shift in the radial velocity of the C IV 5801-12 blend compared with the RVs measured in 1984 and 1993. The results combine to suggest that HD 38030 is a colliding-wind binary having short-lived dust formation episodes, like the Galactic systems WR 140 and WR 19, but at intervals in excess of 20 yr., Comment: 8 pages, 6 figures, accepted for publication in the MNRAS

We present optical and infrared (IR) light curves of NaSt1, also known as Wolf-Rayet (WR) 122, with observations from Palomar Gattini-IR (PGIR), the Zwicky Transient Facility (ZTF), the Katzman Automatic Imaging Telescope (KAIT), the Asteroid Terrestrial-impact Last Alert System (ATLAS), and the All-Sky Automated Survey for Supernovae (ASAS-SN). We identify a $P=309.7\pm0.7$ d photometric period from the optical and IR light curves that reveal periodic, sinusoidal variability between 2014 July and 2021 July. We also present historical IR light curves taken between 1983 July and 1989 May, which show variability consistent with the period of the present-day light curves. In the past, NaSt1 was brighter in the $J$ band with larger variability amplitudes than the present-day PGIR values, suggesting that NaSt1 exhibits variability on longer ($\gtrsim$ decade) timescales. Sinusoidal fits to the recent optical and IR light curves show that the amplitude of NaSt1's variability differs at various wavelengths and also reveal significant phase offsets of $17.0\pm2.5$ d between the ZTF $r$ and PGIR $J$ light curves. We interpret the $310$ d photometric period from NaSt1 as the orbital period of an enshrouded massive binary. We suggest that the photometric variability of NaSt1 may arise from variations in the line-of-sight optical depth toward circumstellar optical/IR emitting regions throughout its orbit due to colliding-wind dust formation. We speculate that past mass transfer in NaSt1 may have been triggered by Roche-lobe overflow (RLOF) during an eruptive phase of an Ofpe/WN9 star. Lastly, we argue that NaSt1 is no longer undergoing RLOF mass transfer., Comment: 20 Pages, 7 Figures, 5 Tables, 1 Table Attachment, Accepted to ApJ

We present spectroscopy of the P~Cygni profile of the 1.083-micron He I line in the WC7 + O5 colliding-wind binary (CWB) WR 140 (HD 193793), observed in 2008, before its periastron passage in 2009, and in 2016-17, spanning the subsequent periastron passage. Both absorption and emission components showed strong variations. The variation of the absorption component as the O5 star was occulted by the wind-collision region (WCR) sets a tight constraint on its geometry. While the sightline to the O5 star traversed the WCR, the strength and breadth of the absorption component varied significantly on time-scales of days. An emission sub-peak was observed on all our profiles. The variation of its radial velocity with orbital phase was shown to be consistent with formation in the WCR as it swung round the stars in their orbit. Modelling the profile gave a measure of the extent of the sub-peak forming region. In the phase range 0.93-0.99, the flux in the sub-peak increased steadily, approximately inversely proportionally to the stellar separation, indicating that the shocked gas in the WCR where the line was formed was adiabatic. After periastron, the sub-peak flux was anomalously strong and varied rapidly, suggesting formation in clumps down-stream in the WCR. For most of the time, its flux exceeded the 2-10-keV X-ray emission, showing it to be a significant coolant of the shocked wind., Comment: 19 pages, 14 figures, accepted for publication in the MNRAS

We present Spitzer/IRAC observations of dust formation from six extragalactic carbon-rich Wolf-Rayet (WC) binary candidates in low-metallicity (Z $\lesssim0.65$ Z$_\odot$) environments using multi-epoch 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 $\lambda5801\text{-}12$~and/or the C III-IV $\lambda4650$ He II $\lambda4686$~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 sub-solar 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 one of most IR-luminous dust-forming WC systems known, WR~104 (M$_\mathrm{[3.6]}\lesssim-12.3$). An analysis of dust formation in the mid-IR outburst from SPIRITS~19q reveals a high dust production rate of $\dot{M}_d\gtrsim2\times10^{-6}$ M$_\odot$ 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 a high mass-loss rate ($\dot{M}\gtrsim1.6\times10^{-6}$ M$_\odot$ yr$^{-1}$). This efficient dust-formation from early-type WC binaries highlights their potential role as significant sources of dust in low-metallicity environments., Comment: 28 pages, 10 Figures, and 5 Tables. Accepted to ApJ on 2021 Jan 4 (Replacement of version submitted to arXiv on 2020 Nov 19)

WR 112 is a dust-forming carbon-rich Wolf-Rayet (WC) binary with a dusty circumstellar nebula that exhibits a complex asymmetric morphology, which traces the orbital motion and dust formation in the colliding winds of the central binary. Unraveling the complicated circumstellar dust emission around WR 112 therefore provides an opportunity to understand the dust formation process in colliding-wind WC binaries. In this work, we present a multi-epoch analysis of the circumstellar dust around WR 112 using seven high spatial resolution (FWHM $\sim0.3-0.4''$) N-band ($\lambda \sim12$ $\mu$m) imaging observations spanning almost 20 years and includes newly obtained images from Subaru/COMICS in Oct 2019. In contrast to previous interpretations of a face-on spiral morphology, we observe clear evidence of proper motion of the circumstellar dust around WR 112 consistent with a nearly edge-on spiral with a $\theta_s=55^\circ$ half-opening angle and a $\sim20$-yr period. The revised near edge-on geometry of WR 112 reconciles previous observations of highly variable non-thermal radio emission that was inconsistent with a face-on geometry. We estimate a revised distance to WR 112 of $d = 3.39^{+0.89}_{-0.84}$ kpc based on the observed dust expansion rate and a spectroscopically derived WC terminal wind velocity of $v_\infty= 1230\pm260$ km s$^{-1}$. With the newly derived WR 112 parameters we fit optically-thin dust spectral energy distribution models and determine a dust production rate of $\dot{M}_d=2.7^{+1.0}_{-1.3}\times10^{-6}$ M$_\odot$ yr$^{-1}$, which demonstrates that WR 112 is one of the most prolific dust-making WC systems known., Comment: 17 pages, 9 figures, 1 animated gif, accepted for publication in ApJ

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 $\dot{M}_d\sim10^{-10}-10^{-6}$ M$_\odot$ yr$^{-1}$ and carbon dust condensation fractions, $\chi_C$, between $0.002 - 40\%$. A large ($0.1 - 1.0$ $\mu$m) dust grain size composition is favored for efficient dustars where $\chi_C\gtrsim1\%$. Results for dustars with known orbital periods verify a power-law relation between $\chi_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 co-eval $10^6$ M$_\odot$ 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 co-eval models show 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 metallicites WC binaries are the dominant source of dust for $\sim60$ Myr until the onset of AGBs, which match the dust input of WC binaries. Co-eval population models show that for "bursty" SF, AGB stars dominate dust production at late times ($t\gtrsim 70$ Myr)., Comment: 34 pages, 12 figures, 7 tables, Accepted for publication in ApJ

Photometry at 3.4 and 4.6 micron of 128 Population~I WC type Wolf-Rayet stars in the Galaxy and 12 in the Large Magellanic Cloud (LMC) observed in the WISE NEOWISE-R survey was searched for evidence of circumstellar dust emission and its variation. Infrared spectral energy distributions (SEDs) were assembled, making use of archival r, i, Z and Y photometry to determine reddening and stellar wind levels for the WC stars found in recent IR surveys and lacking optical photometry. From their SEDs, ten apparently non-variable stars were newly identified as dust makers, including three, WR 102-22, WR 110-10 and WR 124-10, having subtype earlier than WC8-9, the first such stars to show this phenomenon. The 11 stars found to show variable dust emission include six new episodic dust-makers, WR 47c, WR 75-11, WR 91-1, WR 122-14 and WR 125-1 in the Galaxy and HD 38030 in the LMC. Of previously known dust makers, NEOWISE-R photometry of WR 19 captured its rise to maximum in 2018, confirming the 10.1-y period, that of WR 125 the beginning of a new episode of dust formation suggesting a period near 28.3~y., while that of HD 36402 covered almost a whole period and forced revision of it to 5.1~y., Comment: accepted to MNRAS

We report on interferometric observations with the CHARA Array of two classical Wolf-Rayet stars in suspected binary systems, namely WR 137 and WR 138. In both cases, we resolve the component stars to be separated by a few milliarcseconds. The data were collected in the H-band, and provide a measure of the fractional flux for both stars in each system. We find that the WR star is the dominant H-band light source in both systems ($f_{\rm WR, 137} = 0.59\pm0.04$; $f_{\rm WR, 138} = 0.67\pm0.01$), which is confirmed through both comparisons with estimated fundamental parameters for WR stars and O dwarfs, as well as through spectral modeling of each system. Our spectral modeling also provides fundamental parameters for the stars and winds in these systems. The results on WR 138 provide evidence that it is a binary system which may have gone through a previous mass-transfer episode to create the WR star. The separation and position of the stars in the WR 137 system together with previous results from the IOTA interferometer provides evidence that the binary is seen nearly edge-on. The possible edge-on orbit of WR 137 aligns well with the dust production site imaged by the Hubble Space Telescope during a previous periastron passage, showing that the dust production may be concentrated in the orbital plane., Comment: 11 pages, 4 tables, 7 figures, accepted to MNRAS

Suzaku observations of the Wolf-Rayet binary WR 140 (WC7pd+O5.5fc) were made at four different times around periastron passage in 2009 January. The spectra changed in shape and flux with the phase. As periastron approached, the column density of the low-energy absorption increased, which indicates that the emission from the wind-wind collision plasma was absorbed by the dense W-R wind. The spectra can be mostly fitted with two different components: a warm component with kT=0.3--0.6 keV and a dominant hot component with kT~3 keV. The emission measure of the dominant, hot component is not inversely proportional to the distance between the two stars. This can be explained by the O star wind colliding before it has reached its terminal velocity, leading to a reduction in its wind momentum flux. At phases closer to periastron, we discovered a cool plasma component in a recombining phase, which is less absorbed. This component may be a relic of the wind-wind collision plasma, which was cooled down by radiation, and may represent a transitional stage in dust formation., Comment: 16 pages, 16 figures, accepted for publication in PASJ

We present infrared photometry of the WC8 Wolf-Rayet system WR 48a observed with telescopes at ESO, the SAAO and the AAT between 1982 and 2011 which show a slow decline in dust emission from the previously reported outburst in 1978--79 until about 1997, when significant dust emission was still evident. This was followed by a slow rise, accelerating to reach and overtake the first (1978) photometry, demonstrating that the outburst observed in 1978--79 was not an isolated event, but that they recur at intervals of 32+ years. This suggests that WR 48a is a long-period dust maker and colliding-wind binary (CWB). The locus of WR 48a in the (H-L), K colour-magnitude diagram implies that the rate of dust formation fell between 1979 and about 1997 and then increased steadily until 2011. Superimposed on the long-term variation are secondary (`mini') eruptions in (at least) 1990, 1994, 1997, 1999 and 2004, characteristic of relatively brief episodes of additional dust formation. Spectra show evidence for an Oe or Be companion to the WC8 star, supporting the suggestion that WR 48a is a binary system and indicating a system luminosity consistent with the association of WR 48a and the young star clusters Danks 1 and Danks 2. The range of dust formation suggests that these stars are in an elliptical orbit having e ~ 0.6. The size of the orbit implied by the minimum period, together with the WC wind velocity and likely mass-loss rate, implies that the post-shock WC wind is adiabatic throughout the orbit -- at odds with the observed dust formation. A similar conflict is observed in the `pinwheel' dust-maker WR 112., Comment: 15 pages, 8 figures, accepted for publication in MNRAS

A survey for H2 1-0 S(1) emission in 16 Herbig-Haro (HH) objects and three exciting stars for HH objects is reported. Eleven HH objects which show low-excitation optical spectra exhibit H2 emission. One object (HH 43) is more than twice as bright as any previously reported HH object. In addition, spectra in the range 1.6-2.55 microns are reported for HH 43 and HH 120, and a 2.0-2.55 micron spectrum is presented for HH 26. The spectra yield estimates of the H2 density and temperature ranges in these objects. The role of ultraviolet H2 emission-line fluorescence in HH 43 with respect to cascading among excited vibrational states of the ground electronic state is discussed. Models which may account for the combined ultraviolet, optical, and near-IR spectra of HHs are briefly analyzed.