Circumstellar Dust Research Articles

Radiation-driven acceleration in the expanding WR140 dust shell

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 companion2 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 close2–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.

Nested dust shells around the Wolf–Rayet binary WR 140 observed with JWST

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.

Modeling Extinction and Reddening Effects by Circumstellar Dust in the Betelgeuse Envelope in the Presence of Radiative Torque Disruption

Circumstellar dust forms and evolves within the envelope of evolved stars, including asymptotic giant branch (AGB) and red supergiant (RSG) stars. The extinction of stellar light by circumstellar dust is vital for interpreting RSG/AGB observations and determining high-mass RSG progenitors of core-collapse supernovae. However, circumstellar dust properties are not well understood. Modern understanding of dust evolution suggests that intense stellar radiation can radically change dust properties across the circumstellar envelope through the RAdiative Torque Disruption (RAT-D) mechanism. In this paper, we study the impacts of RAT-D on the grain size distribution (GSD) of circumstellar dust and model its effects on photometric observations of α Orionis (Betelgeuse). Due to the RAT-D effects, large grains formed in the dust formation zone are disrupted into smaller species of size a < 0.5 μm. Using the GSD constrained by the RAT-D effects, we model the visual extinction of background stars and Betelgeuse. We find that the extinction decreases at near-UV, optical, and IR wavelengths while increasing at far-UV wavelengths. The resulting flux potentially reproduces the observation from the near-UV to the near-IR range. Our results can be used to explain dust extinction and photometric observations of other RSG/AGB stars.

Optical Time-series Photometry of the Symbiotic Nova V1835 Aquilae

We present time-series CCD photometry in the BVRI passbands of the recently identified symbiotic nova V1835 Aquilae (NSV 11749) over an interval of 5.1 yr with 7–14 day cadence, observed during its quiescence. We find slow light variations with a range of ∼0.9 mag in V and ∼0.3 mag in I. Analysis of these data show strong periodicity at 419 ± 10 days, which we interpret to be the system’s orbital period. A dip in the otherwise-sinusoidal phased light curve suggests a weak ellipsoidal effect due to tidal distortion of the giant star, which in turn opens the possibility that V1835 Aql transfers some of its mass to the hot component via Roche lobe overflow rather than via a stellar wind. We also find evidence that V1835 Aql is an S-type symbiotic star, relatively free of circumstellar dust, and include it among the nuclear burning group of symbiotics. Finally, we provide photometry, periods, and light curve classifications for 22 variable stars in the field around V1835 Aql, about half of which are newly identified.

Mira variables in the Milky Way’s nuclear stellar disc: discovery and classification

ABSTRACT The properties of the Milky Way’s nuclear stellar disc give crucial information on the epoch of bar formation. Mira variables are promising bright candidates to study the nuclear stellar disc, and through their period–age relation dissect its star formation history. We report on a sample of 1782 Mira variable candidates across the central $3\times 3\, \mathrm{deg}^2$ of the Galaxy using the multi-epoch infrared VISTA Variables in Via Lactea (VVV) survey. We describe the algorithms employed to select candidate variable stars and then model their light curves using periodogram and Gaussian process methods. By combining with WISE, 2MASS, and other archival photometry, we model the multiband light curves to refine the periods and inspect the amplitude variation between different photometric bands. The infrared brightness of the Mira variables means many are too bright and missed by VVV. However, our sample follows a well-defined selection function as expected from artificial star tests. The multiband photometry is modelled using stellar models with circumstellar dust that characterize the mass-loss rates. We demonstrate how ≳90 per cent of our sample is consistent with O-rich chemistry. Comparison to period–luminosity relations demonstrates that the bulk of the short period stars are situated at the Galactic Centre distance. Many of the longer period variables are very dusty, falling significantly under the O-rich Magellanic Cloud and solar neighbourhood period–luminosity relations and exhibit high mass-loss rates of $\sim 2.5\times 10^{-5}M_\odot \, \mathrm{yr}^{-1}$. The period distribution appears consistent with the nuclear stellar disc forming $\gtrsim 8\, \mathrm{Gyr}$ ago, although it is not possible to disentangle the relative contributions of the nuclear stellar disc and the contaminating bulge.

An ALMA 1.3 millimeter Search for Debris Disks around Solar-type Stars in the Pleiades

Millimeter emission from debris disks around stars of different ages provides constraints on the collisional evolution of planetesimals. We present ALMA 1.3 millimeter observations of a sample of 76 Solar-type stars in the ∼115 Myr old Pleiades star cluster. These ALMA observations complement previous infrared observations of this sample by providing sensitivity to emission from circumstellar dust at lower temperatures, corresponding to debris at radii comparable to the Kuiper Belt and beyond. The observations obtain a beam size of 1.″5 (200 au) and a median rms noise of 54 μJy beam−1, which corresponds to a fractional luminosity L dust/L star ∼ 10−4 for 40 K dust for a typical star in the sample. The ALMA images show no significant detections of the targeted stars. We interpret these limits in the context of a steady-state collisional cascade model for debris disk evolution that provides a good description of observations of the field population near the Sun but is not well-calibrated on younger populations. The ALMA nondetections of the Pleiades systems are compatible with the disk flux predictions of this model. We find no high fractional luminosity outliers from these ALMA data that could be associated with enhanced collisions resulting from activity not accounted for by steady-state evolution. However, we note that two systems (H ii 1132 and HD 22680) show 24 μm excess much higher than the predictions of this model, perhaps due to unusually high dust production from dynamical events involving planets.

Infrared Excess of a Large OB Star Sample

The infrared (IR) excess from OB stars is commonly considered to be a contribution from ionized stellar wind or circumstellar dust. With the newly published Large Sky Area Multi-Object Fiber Spectroscopy Telescope (LAMOST)-OB catalog and Galactic O-Star Spectroscopic Survey data, this work steps further on understanding the IR excess of OB stars. Based on a forward-modeling approach comparing the spectral slope of observational spectral energy distributions and photospheric models, 1147 stars are found to have IR excess out of 7818 stars with good-quality photometric data. After removing the objects in the sightline of dark clouds, 532 (∼7%) B-type stars and 118 (∼23%) O-type stars are identified to be true OB stars with circumstellar IR excess emission. The ionized stellar wind model and the circumstellar dust model are adopted to explain the IR excess, and Bayes factors are computed to quantitatively compare the two. It is shown that the IR excess can be accounted for by the stellar wind for about 65% cases, of which 33% by free–free emission and 32% by synchrotron radiation. Other 30% sources could have and 4% should have a dust component or other mechanisms to explain the sharp increase in flux at λ > 10 μm. The parameters of the dust model indicate a large-scale circumstellar halo structure, which implies the origin of the dust from the birthplace of the OB stars. A statistical study suggests that the proportion with IR excess in OB stars increases with the stellar effective temperature and luminosity, and that there is no systematic change in the mechanism for IR emission with stellar parameters.

On the provenance of GEMS, a quarter century post discovery

The provenance of GEMS (glass with embedded metal and sulfides) in cometary type interplanetary dust particles is investigated using analytical scanning transmission electron microscopy and secondary ion mass spectrometry. We review the current state of knowledge and closely examine the densities, elemental compositions and distributions, iron oxidation states, and isotopic compositions of a subset of GEMS in chondritic porous interplanetary dust. We find that GEMS are underdense with estimated densities that are 35–65% of compositionally equivalent crystalline aggregates. GEMS low densities result in a lower contribution to the bulk compositions of IDPs than has been assumed based on their volume fraction. We also find that element/Si ratios, assumed to be primary (indigenous), are instead perturbed by contamination and secondary alteration, including pulse heating during atmospheric entry. Fe in pyrrhotite inclusions was oxidized and Mg, S, Ca, and Fe were depleted relative to lithophile Al and Si, resulting in reduction in element/Si ratios. Because they trap outgassing elements, Fe-rich oxide rims that formed on the surface of GEMS are serendipitous “witness plates” to the changes in composition that accompany atmospheric entry. As a result of alteration, GEMS elemental compositions cannot reliably inform about their provenance. Except for highly anomalous oxygen isotope ratios measured in some large GEMS grains that indicate a contribution from circumstellar dust, oxygen isotope compositions are generally poor indicators of provenance. Prior work indicates that most GEMS fall close to the terrestrial oxygen isotope composition, which, however, does not exclude a presolar interstellar origin. Nitrogen isotopic compositions are more diagnostic. Elevated 15N/14N ratios indicate that GEMS accreted in conjunction with formation of organic matter by ion–molecule reactions in a cold (<50 K) presolar environment like the extreme outer nebula or interstellar medium. Considering all observations, we conclude that GEMS are most likely processed interstellar silicates.

Faint Debris Disk Peering through the Superflare Light Echo

We present the detectability of strong mid-infrared (MIR) light echoes from faint debris disks illuminated by bright superflares of M-dwarf stars. Circumstellar dust grains around an M-dwarf star are simultaneously heated by superflare radiation. One can thus expect their re-emission in the MIR wavelength regime. According to our model calculations for the Proxima Centauri system, the nearest M-dwarf star system, thermal emission echoes from an inner (r < 1 au) debris disk with a total mass down to that of the solar system’s zodiacal dust are expected to emerge at wavelengths longer than ∼10 μm with a strength comparable to or greater than a white-light superflare. Also, observable echoes from inner- (r ≲ 0.5 au) debris disks irradiated by energetic ( ≳1033.5 erg) superflares of nearby (D < 3 pc) M dwarfs are expected. Our simulation results indicate that superflare monitoring using high-speed optical instruments like OASES and its prompt follow-up using ground-based MIR instruments, such as TAO/MIMIZUKU, can detect these MIR light echoes from debris disks around solar neighborhood flare stars.

Inefficient Growth of SiOx Grains: Implications for Circumstellar Outflows

To explain observations of abundant circumstellar dust and high stellar wind velocity, most models simply postulate the efficient nucleation and growth of silicate dust particles. Here, we report measurement of the SiO–(SiOx)n grain sticking coefficient in a microgravity sounding rocket experiment, indicating very inefficient (0.005–0.016) grain formation from the vapor. Application of this measurement to radiative-driven winds in oxygen-rich asymptotic giant branch stars indicates that the initial grain condensate population should consist of very tiny dust particles in very large numbers. Aggregation of this dust population will produce low-dimension fractal aggregates that should couple well to the stellar radiation field and efficiently drive stellar mass loss.

The Extreme Scarcity of Dust-enshrouded Red Supergiants: Consequences for Producing Stripped Stars via Winds

Quiescent mass loss during the red supergiant (RSG) phase has been shown to be far lower than prescriptions typically employed in single-star evolutionary models. Importantly, RSG winds are too weak to drive the production of Wolf-Rayet (WR) stars and stripped-envelope supernovae (SE-SNe) at initial masses of roughly 20–40 M ⊙. If single stars are to make WR stars and SE-SNe, this shifts the burden of mass loss to rare dust-enshrouded RSGs (DE-RSGs), objects claimed to represent a short-lived, high-mass-loss phase. Here, we take a fresh look at the purported DE-RSGs. By modeling the mid-IR excesses of the full sample of RSGs in the Large Magellanic Cloud, we find that only one RSG has both a high mass-loss rate ( ≥ 10−4 M ⊙ yr−1) and a high optical circumstellar dust extinction (7.92 mag). This RSG is WOH G64, and it is the only one of the 14 originally proposed DE-RSGs that is actually dust enshrouded. The rest appear to be either normal RSGs without strong IR excess, or lower-mass asymptotic giant branch stars. Only one additional object in the full catalog of RSGs (not previously identified as a DE-RSG) shows strong mid-IR excess. We conclude that if DE-RSGs do represent a pre-SN phase of enhanced in single stars, it is extremely short-lived, only capable of removing ≤2 M ⊙ of material. This rules out the single-star post-RSG pathway for the production of WR stars, luminous blue variables, and SE-SNe. Single-star models should not employ -prescriptions based on these extreme objects for any significant fraction of the RSG phase.

The Effects of Circumstellar Dust Scattering on the Light Curves and Polarizations of Type Ia Supernovae* *Supported by the National Natural Science Foundation of China.

Observational signatures of the circumstellar material (CSM) around Type Ia supernovae (SNe Ia) provide a unique perspective on their progenitor systems. The pre-supernova evolution of the SN progenitors may naturally eject CSM in most of the popular scenarios of SN Ia explosions. In this study, we investigate the influence of dust scattering on the light curves and polarizations of SNe Ia. A Monte Carlo method is constructed to numerically solve the process of radiative transfer through the CSM. Three types of geometric distributions of the CSM are considered: spherical shell, axisymmetric disk, and axisymmetric shell. We show that both the distance of the dust from the SN and the geometric distribution of the dust affect the light curve and color evolutions of SN. We found that the geometric location of the hypothetical circumstellar dust may not be reliably constrained based on photometric data alone, even for the best observed cases such as SN 2006X and SN 2014J, due to the degeneracy of CSM parameters. Our model results show that a time sequence of broadband polarimetry with appropriate time coverage from a month to about one year after explosion can provide unambiguous limits on the presence of circumstellar dust around SNe Ia.

The Great Dimming of Betelgeuse seen by the Himawari-8 meteorological satellite

Betelgeuse, one of the most studied red supergiant stars1,2, dimmed in the optical range by ~1.2 mag between late 2019 and early 2020, reaching a historical minimum3–5 called ‘the Great Dimming’. Thanks to enormous observational effort to date, two hypotheses remain that can explain the Dimming1: a decrease in the effective temperature6,7 and an enhancement of the extinction caused by newly produced circumstellar dust8,9. However, the lack of multiwavelength monitoring observations, especially in the mid-infrared, where emission from circumstellar dust can be detected, has prevented us from closely examining these hypotheses. Here we present 4.5 yr, 16-band photometry of Betelgeuse between 2017 and 2021 in the 0.45–13.5 μm wavelength range making use of images taken by the Himawari-810 geostationary meteorological satellite. By examining the optical and near-infrared light curves, we show that both a decreased effective temperature and increased dust extinction may have contributed by almost equal amounts to the Great Dimming. Moreover, using the mid-infrared light curves, we find that the enhanced circumstellar extinction actually contributed to the Dimming. Thus, the Dimming event of Betelgeuse provides us with an opportunity to examine the mechanism responsible for the mass loss of red supergiants, which affects the fate of massive stars as supernovae11.

Planetesimals around stars with TESS (PAST) – II. An M dwarf ‘dipper’ star with a long-lived disc in the TESS continuous viewing zone

ABSTRACT Studies of T Tauri discs inform planet formation theory; observations of variability due to occultation by circumstellar dust are a useful probe of unresolved, planet-forming inner discs, especially around faint M dwarf stars. We report observations of 2M0632, an M dwarf member of the Carina young moving group that was observed by Transiting Exoplanet Survey Satellite over two 1-yr intervals. The combined light curve contains &amp;gt;300 dimming events, each lasting a few hours, and as deep as 40 per cent (0.55 magnitudes). These stochastic events are correlated with a distinct, stable 1.86-d periodic signal that could be stellar rotation. Concurrent ground-based, multiband photometry show reddening consistent with interstellar medium-like dust. The star’s excess emission in the infrared and emission lines in optical and infrared spectra reveal a T Tauri-like accretion disc around the star. We confirm membership of 2M0632 in the Carina group by a Bayesian analysis of its Galactic space motion and position. We combine stellar evolution models with Gaia photometry and constraints on Teff, luminosity, and the absence of detectable lithium in the photosphere to constrain the age of the group and 2M0632 to 40–60 Myr, consistent with earlier estimates. 2M0632 joins a handful of long-lived discs which challenge the canon that disc lifetimes are ≲10 Myr. All known examples surround M dwarfs, suggesting that lower X-ray/ultraviolet irradiation and slower photoevaporation by these stars can dramatically affect disc evolution. The multiplanet systems spawned by long-lived discs probably experienced significant orbital damping and migration into close-in, resonant orbits, and perhaps represented by the TRAPPIST-1 system.

THA 15−31: Discovery with VLT/X-shooter and Swift/UVOT of a new symbiotic star of the accreting-only variety

We report the discovery and characterization of a new symbiotic star of the accreting-only variety, which we observed in the optical/near-infrared (NIR) with VLT/X-shooter and in the X-rays/ultraviolet with Swift/UVOT+XRT. The new symbiotic star, THA 15−31, was previously described as a pre-main sequence star belonging to the Lupus 3 association. Our observations, ancillary data, and Gaia EDR3 parallax indicate that THA 15−31 is a symbiotic star composed of an M6III red giant and an accreting companion, is subject to EB − V = 0.38 reddening, and is located at a distance of ∼12 kpc and at 1.8 kpc above the Galactic plane in the outskirts of the Bulge. The luminosity of the accreting companion is ∼100 L⊙, placing THA 15−31 among the symbiotic stars accreting at a high rate (2.5 × 10−8 M⊙ yr−1 if the accretion is occurring on a white dwarf of 1 M⊙). The observed emission lines originate primarily from H I, He I, and Fe II, with no He II or other high-excitation lines observed; a sharp central absorption superimposed on the Balmer emission lines is observed, while all other lines have a simple Gaussian-like profile. The emission from the companion dominates over the M6III red giant at U and B-band wavelengths, and is consistent with an origin primarily in an optically thick accretion disk. No significant photometric variability is observed at optical or NIR wavelengths, suggesting either a face-on orbital orientation and/or that the red giant is far from Roche-lobe filling conditions. The profile of emission lines supports a low orbital inclination if they form primarily in the accretion disk. An excess emission is present in AllWISE W3 (12 μm) and W4 (22 μm) data, radiating a luminosity ≥35 L⊙, consistent with thermal emission from optically thin circumstellar dust.

Scattered Radiation of Protoplanetary Disks

Scattered radiation of circumstellar (CS) dust plays an important role in the physics of young stars. Its observational manifestations are various but more often they are connected with the appearance of intrinsic polarization in young stars and their CS disks. In our brief review we consider two classes of astrophysical objects in which the participation of scattered radiation is key for understanding their nature. First of all, these are irregular variables (UX Ori type stars). The modern idea of their nature and the mechanism of their variability has been formed thanks to synchronous observations of their linear polarization and brightness. The second class of objects is the CS disks themselves. Their detailed investigation became possible due to observations in polarized light using a coronographic technique and large telescopes.

Quasi-periodic Dimming of the ∼130 Myr-old Debris-Disk Hosting Star HD 240779 is not Persistent

Abstract Quasi-periodic (1.5 days) dimming (by circumstellar dust) of the 135 Myr-old AB Doradus moving group member HD 240779 was detected in photometry by TESS in late 2018. Similar observations two years later show no such signal, and ground-based photometry indicate that the signal was absent in late 2019. This suggests that the source of the dust did not survive long after 2018, e.g., it was a disrupted planetesimal, or that dust production by the body is episodic, analogous to the “evaporating” planets detected by Kepler.

Young Exoplanet Transit Initiative follow-up observations of the T Tauri star CVSO 30 with transit-like dips

ABSTRACT The T Tauri star CVSO 30, also known as PTFO 8-8695, was studied intensively with ground-based telescopes as well as with satellites over the last decade. It showed a variable light curve with additional repeating planetary transit-like dips every ∼10.8 h. However, these dimming events changed in depth and duration since their discovery and from autumn 2018 on, they were not even present or near the predicted observing times. As reason for the detected dips and their changes within the complex light curve, e.g. a disintegrating planet, a circumstellar dust clump, stellar spots, possible multiplicity, and orbiting clouds at a Keplerian co-rotating radius were discussed and are still under debate. In this paper, we present additional optical monitoring of CVSO 30 with the meter class telescopes of the Young Exoplanet Transit Initiative in Asia and Europe over the last 7 yr and characterize CVSO 30 with the new Early Data Release 3 of the European Space Agency-Gaia mission. As a result, we describe the evolution of the dimming events in the optical wavelength range since 2014 and present explanatory approaches for the observed variabilities. We conclude that orbiting clouds of gas at a Keplerian co-rotating radius are the most promising scenario to explain most changes in CVSO 30’s light curve.

SOFIA/FORCAST Monitoring of the Dust Emission from R Aqr: Start of the Eclipse

Abstract We present mid-infrared (mid-IR) spectra from our continued monitoring of R Aquarii, the nearest symbiotic Mira, using the Stratospheric Observatory for Infrared Astronomy (SOFIA). New photometric and spectroscipic data were obtained with the Faint Object infraRed CAmera for the SOFIA Telescope in 2018 and 2019 after the system had started its “eclipse,” during which it became two magnitudes fainter in the visual. The mid-IR flux, in particular the 10 μm silicate feature, has strengthened compared with the previous cycles. Radiative transfer models for the circumstellar dust emission were calculated for the new spectra, and recalculated for those previously obtained using more appropriate values of the near-IR magnitudes to constrain the properties of the asymptotic giant branch spectra heating the dust. The modeling shows that the luminosity dependence on pulsation phase is not affected by the onset of the eclipse, and that the increase in the mid-IR flux is due to a higher dust density. The models also confirm our earlier results that micron-size grains are present, and that no changes in the grain composition are required to explain the variations in the spectra.

A Fresh Look at AGB Stars in Galactic Open Clusters with Gaia: Impact on Stellar Models and the Initial–Final Mass Relation

Abstract Benefiting from the Gaia second and early third releases of photometric and astrometric data, we examine the population of asymptotic giant branch (AGB) stars that appear in the fields of intermediate-age and young open star clusters. We identify 49 AGB star candidates, brighter than the tip of the red giant branch, with a good to high cluster membership probability. Among them, we find 19 TP-AGB stars with known spectral type: 4 M stars, 3 MS/S stars, and 12 C stars. By combining observations, stellar models, and radiative transfer calculations that include the effect of circumstellar dust, we characterize each star in terms of initial mass, luminosity, mass-loss rate, core mass, period, and mode of pulsation. The information collected helps us shed light on the TP-AGB evolution at solar-like metallicity, placing constraints on the third dredge-up process, the initial masses of carbon stars, stellar winds, and the initial–final mass relation (IFMR). In particular, we find that two bright carbon stars, MSB 75 and BM IV 90, members of the clusters NGC 7789 and NGC 2660 (with similar ages of ≃ 1.2–1.6 Gyr and initial masses 2.1 ≳ M i /M ⊙ ≳ 1.9), have unusually high core masses, M c ≈ 0.67–0.7 M ⊙. These results support the findings of a recent work (Marigo et al. 2020) that identified a kink in the IFMR, which interrupts its monotonic trend just at the same initial masses. Finally, we investigate two competing scenarios to explain the M c data: the role of stellar winds in single-star evolution, and binary interactions through the blue straggler channel.