Apparent Motion of the Circumstellar Envelope of CQ Tau in Scattered Light

The study of spiral structures in protoplanetary disks is of great importance for understanding the processes in the disks, including planet formation. Bright spiral arms were detected in the disk of young star CQ Tau by Uyama et al. in the H and L bands. The spiral arms are located inside the gap in millimeter-sized dust, discovered earlier using Atacama Large Millimeter/submillimeter Array observations. To explain the gap, Ubeira Gabellini et al. proposed the existence of a planet with the semimajor axis of 20 au. We obtained multi-epoch observations of a spiral feature in the circumstellar envelope of CQ Tau in the I c band using a novel technique of differential speckle polarimetry. The observations covering a period from 2015 to 2021 allow us to estimate the pattern speed of the spiral: −0.°2 ± 1.°1 yr−1 (68% credible interval; positive value indicates counterclockwise rotation), assuming a face-on orientation of the disk. This speed is significantly smaller than expected for a companion-induced spiral, if the perturbing body has a semimajor axis of 20 au. We emphasize that the morphology of the spiral structure is likely to be strongly affected by shadows of a misaligned inner disk detected by Eisner et al.

A study of many-year spectral variability of the Ae Herbig star HD 36112

We present the results of long-term high-dispersion spectral observations (R = 20000) of the Ae Herbig star HD 36112 in the regions of the Ha emission line and the NaI D resonance doublet lines. They show that parameters of the Ha emission line demonstrate complicated variability on several time scales: 1) variability from night to night caused by inhomogeneity of the circumstellar envelope; 2) variability on a time scale of about 1200d characterized by a variation of the equivalent width, intensity, and other emission parameters; 3) variability on a time scale of more than 4000d observed as a many-year trend in variations of parameters of the Ha emission line. We associate these results with variability of physical and kinematic conditions in the inner regions of the accretion disk and wind. The most probable mechanism of this variability is a process of planet formation in the circumstellar disk.

SOFIA upGREAT/FIFI-LS Emission-line Observations of Betelgeuse during the Great Dimming of 2019/2020

We report NASA-DLR SOFIA upGREAT circumstellar [O i] 63.2 μm and [C ii] 157.7 μm emission profiles and FIFI-LS [O i] 63.2 μm, [O i] 145.5 μm, and [C ii] 157.7 μm fluxes obtained shortly after Betelgeuse’s 2019/2020 Great Dimming event. Haas et al. noted a potential correlation between the [O i] 63.2 μm flux and V magnitude based on three Kuiper Airborne Observatory observations made with the CGS and FIFI instruments. The FIFI observation was obtained when V ≃ 0.88 and revealed a 3σ non-detection at a quarter of the previous CGS flux measurement made when V ≃ 0.35. A potential explanation could be a change in dust-gas drag heating by circumstellar silicates caused by variations in the photospheric radiation field. SOFIA observations provide a unique test of this correlation because the V-band brightness went to its lowest value on record, V ≃ 1.61, with the SOFIA observations being made when V FIFI−LS ≃ 1.51 and V upGREAT ≃ 1.36. The upGREAT spectra show a [O i] 63.2 μm flux larger than previous space observatory measurements obtained when V ≃ 0.58. The profile is consistent with formation in the slower, more turbulent inner S1 outflow, while the [C ii] 157.7 μm profile is consistent with formation farther out in the faster S2 outflow. Modeling of dust-gas drag heating, combined with 25 yr of Wing three-filter and V photometry, reveals that it is unlikely that the S1 circumstellar envelope and [O i] 63.2 μm fluxes are dominated by the dust-gas drag heating and that another heating source is also active. The [O i] 63.2 μm profile is hard to reconcile with existing outflow velocity models.

Maser emission from the CO envelope of the asymptotic giant branch star W Hydrae

Context. Observation of CO emission around asymptotic giant branch (AGB) stars is the primary method to determine gas mass-loss rates. While radiative transfer models have shown that molecular levels of CO can become mildly inverted, causing maser emission, CO maser emission has yet to be confirmed observationally. Aims. High-resolution observations of the CO emission around AGB stars now have the brightness temperature sensitivity to detect possible weak CO maser emission. Methods. We used high angular resolution observations taken with the Atacama Large Millimeter/submillimeter Array (ALMA) to observe the small-scale structure of CO J = 3−2 emission around the oxygen-rich AGB star W Hya. Results. We find CO maser emission amplifying the stellar continuum with an optical depth τ ≈−0.55. The maser predominantly amplifies the limb of the star because CO J = 3−2 absorption from the extended stellar atmosphere is strongest towards the centre of the star. Conclusions. The CO maser velocity corresponds to a previously observed variable component of high-frequency H2O masers and with the OH maser that was identified as the amplified stellar image. This implies that the maser originates beyond the acceleration region and constrains the velocity profile since we find the population inversion primarily in the inner circumstellar envelope. We find that inversion can be explained by the radiation field at 4.6 μm and that the existence of CO maser emission is consistent with the estimated mass-loss rates for W Hya. However, the pumping mechanism requires a complex interplay between absorption and emission lines in the extended atmosphere. Excess from dust in the circumstellar envelope of W Hya is not sufficient to contribute significantly to the required radiation field at 4.6 μm. The interplay between molecular lines that cause the pumping can be constrained by future multi-level CO observations.

SN 2018bsz: significant dust formation in a nearby superluminous supernova

We investigate the thermal emission and extinction from dust associated with the nearby superluminous supernova (SLSN) 2018bsz. Our dataset has daily cadence and simultaneous optical and near-infrared coverage up to ~ 100 days, together with late time (+1.7 yr) MIR observations. At 230 days after light curve peak the SN is not detected in the optical, but shows a surprisingly strong near-infrared excess, with r - J > 3 mag and r - Ks > 5 mag. The time evolution of the infrared light curve enables us to investigate if the mid-infrared emission is from newly formed dust inside the SN ejecta, from a pre-existing circumstellar envelope, or interstellar material heated by the radiation from the SN. We find the latter two scenarios can be ruled out, and a scenario where new dust is forming in the SN ejecta at epochs > 200 days can self-consistently reproduce the evolution of the SN flux. We can fit the spectral energy distribution well at +230d with 5x10-4 solar mass of carbon dust, increasing over the following several hundred days to 10-2 solar mass by +535d. SN 2018bsz is the first SLSN showing evidence for dust formation within the SN ejecta, and appears to form ten times more dust than normal core-collapse SNe at similar epochs. Together with their preference for low mass, low metallicity host galaxies, we suggest that SLSNe may be a significant contributor to dust formation in the early Universe.

AGB Stars and Their Circumstellar Envelopes. I. the VULCAN Code

The interplay between AGB interiors and their outermost layers, where molecules and dust form, is a problem of high complexity. As a consequence, physical processes like mass loss, which depend on the chemistry of the circumstellar envelope, are often oversimplified. The best candidates to drive mass-loss in AGB stars are dust grains, which trap the outgoing radiation and drag the surrounding gas. Grains build up, however, is far from being completely understood. Our aim is to model both the physics and the chemistry of the cool expanding layers around AGB stars in order to characterize the on-going chemistry, from atoms to dust grains. This has been our rationale to develop ab initio VULCAN, a FORTRAN hydro code able to follow the propagation of shocks in the circumstellar envelopes of AGB stars. The version presented in this paper adopts a perfect gas law and a very simplified treatment of the radiative transfer effects and dust nucleation. In this paper, we present preliminary results obtained with our code.

How to disentangle geometry and mass-loss rate from AGB-star spectral energy distributions

Context. High-angular-resolution observations of asymptotic giant branch (AGB) stars often reveal non-spherical morphologies for the gas and dust envelopes. Aims. We aim to make a pilot study to quantify the impact of different geometries (spherically symmetric, spiral-shaped, and disc-shaped) of the dust component of AGB envelopes on spectral energy distributions (SEDs), mass estimates, and subsequent mass-loss rate (MLR) estimates. We also estimate the error made on the MLR if the SED is fitted by an inappropriate geometrical model. Methods. We use the three-dimensional Monte-Carlo-based radiative-transfer code RADMC-3D to simulate emission from dusty envelopes with different geometries (but fixed spatial extension). We compare these predictions with each other, and with the SED of the AGB star EP Aqr that we use as a benchmark since its envelope is disc-like and known to harbour spiral arms, as seen in CO. Results. The SEDs involving the most massive envelopes are those for which the different geometries have the largest impact, primarily on the silicate features at 10 and 18 μm. These different shapes originate from large differences in optical depths. Massive spirals and discs appear akin to black bodies. Optically thick edge-on spirals and discs (with dust masses of 10−4 and 10−5 M⊙) exhibit black-body SEDs that appear cooler than those from face-on structures and spheres of the same mass, while optically thick face-on distributions appear as warmer emission. We find that our more realistic models, combined spherical and spiral distributions, are 0.1 to 0.5 times less massive than spheres with similar SEDs. More extreme, less realistic scenarios give that spirals and discs are 0.01 to 0.05 times less massive than corresponding spheres. This means that adopting the wrong geometry for an AGB circumstellar envelope may result in a MLR that is incorrect by as much as one to two orders of magnitude when derived from SED fitting.

Interstellar nitrile anions: Detection of C3N− and C5N− in TMC-1

We report on the first detection of C3N− and C5N− towards the cold dark core TMC-1 in the Taurus region, using the Yebes 40 m telescope. The observed C3N/C3N− and C5N/C5N− abundance ratios are ∼140 and ∼2, respectively; that is similar to those found in the circumstellar envelope of the carbon-rich star IRC +10216. Although the formation mechanisms for the neutrals are different in interstellar (ion-neutral reactions) and circumstellar clouds (photodissociation and radical-neutral reactions), the similarity of the C3N/C3N− and C5N/C5N− abundance ratios strongly suggests a common chemical path for the formation of these anions in interstellar and circumstellar clouds. We discuss the role of radiative electronic attachment, reactions between N atoms and carbon chain anions Cn−, and that of H− reactions with HC3N and HC5N as possible routes to form CnN−. The detection of C5N− in TMC-1 gives strong support for assigning to this anion the lines found in IRC +10216, as it excludes the possibility of a metal-bearing species, or a vibrationally excited state. New sets of rotational parameters have been derived from the observed frequencies in TMC-1 and IRC +10216 for C5N− and the neutral radical C5N.