Search

We present contemporaneous high-angular-resolution millimeter imaging and visible polarimetric imaging of the nearby asymptotic giant branch (AGB) star W Hya to better understand the dynamics and dust formation within a few stellar radii. The star W Hya was observed in two vibrationally excited H2O lines at 268 and 251 GHz with ALMA at a spatial resolution of 16 x 20 mas and at 748 and 820 nm at a resolution of 26 x 27 mas with the VLT/SPHERE-ZIMPOL. ALMA's high spatial resolution allowed us to image strong emission of the vibrationally excited H2O line at 268 GHz (v2 = 2, J_K_a,K_c = 6_5,2 - 7_4,3) over the stellar surface instead of absorption against the continuum, which is expected for thermal excitation. Strong, spotty emission was also detected along and just outside the stellar disk limb at an angular distance of ~40 mas (~1.9 stellar radii), extending to ~60 mas (~2.9 stellar radii). Another H2O line (v2 = 2, J_K_a,K_c = 9_2,8 - 8_3,5) at 251 GHz with a similar upper-level energy was tentatively identified, which shows absorption over the stellar surface. This suggests that the emission over the surface seen in the 268 GHzH2O line is suprathermal or even maser emission. The estimated gas temperature and H2O density are consistent with the radiatively pumped masers. The 268 GHz H2O line reveals global infall at up to ~15 km/s within 2--3 stellar radii, but outflows at up to ~8 km/s are also present. The polarized intensity maps obtained in the visible reveal clumpy dust clouds forming within ~40 mas (~1.9 stellar radii) with a particularly prominent cloud in the SW quadrant and a weaker cloud in the east. The 268 GHz H2O emission overlaps very well with the visible polarized intensity maps, which suggests that both the nonthermal and likely maser H2O emission and the dust originate from dense, cool pockets in the inhomogeneous atmosphere within ~2--3 stellar radii., Comment: 10 pages, 7 figures, published in Astronomy and Astrophysics

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

Our purpose is to study the effect of binary companions located within the first 10 stellar radii from the primary AGB star. In this work, we target the mass-losing carbon star V Hydrae (V Hya), looking for signatures of its companion in the dust forming region of the atmosphere. The star was observed in the L- and N-bands with the VLTI/MATISSE instrument at low spectral resolution. We reconstructed images of V Hya's photosphere and surroundings using the two bands and compared our interferometric observables with VLTI/MIDI and VISIR archival data. To constrain the dust properties, we used DUSTY to model the spectral energy distribution. The star is dominated by dust emission in the L- and N- bands. The VISIR image confirms the presence of a large-scale dusty circumstellar envelope surrounding V Hya. The MATISSE reconstructed images show asymmetric and elongated structures in both infrared bands. In the L-band, we detected an elongated shape of approximately 15 mas, likely to be of photospheric origin. In the N-band, we found a 20 mas extension North-East from the star, and perpendicular to the L-band elongated axis. The position angle and the size of the N-band extension match the prediction of the companion position at MATISSE epoch. By comparing MATISSE N-band with MIDI data, we deduce that the elongation axis in the N-band has rotated since the previous interferometric measurements 13 years ago, supporting the idea that the particle enhancement is related to the dusty clump moving along with the companion. The MATISSE images unveil the presence of a dust enhancement at the companion position, opening new doors for further analysis on the binary interaction with an AGB component., Comment: Accepted for publication in Astronomy & Astrophysics

Eta Carinae underwent the Great Eruption in the 1840s and a Lesser Eruption in the 1890s. Its apparent spectrum, modified by intervening ejecta, the Homunculus and Little Homunculus, continues to evolve but contains information pertaining to events in the 19th century. The LOS spectrum contains narrow absorption velocities, from -122 to -1665 km/s: rungs of a broken ladder caused by shells formed by the interacting winds. Estimated shell origin dates correlate with origin dates of expanding emission structures preceding the Great Eruption. The LOS absorption velocities extend the record post Great Eruption to the Lesser Eruption. We suggest that these shells originated from a binary merger within a triple system. Shells formed not only from periastron passages of the current secondary, but also from ear-like extensions preceding and following the periastron event. Additional models need to be considered., Comment: 18 pages, 5 figures, 1 table

Mass loss at the asymptotic giant branch (AGB) plays an important role not only in the final fates of stars, but also in the chemical evolution of galaxies. Nevertheless, the metallicity effects on AGB mass loss are not yet fully understood. We present spatially resolved observations of an AGB star, V3, in the metal-poor globular cluster 47 Tuc (NGC 104). The AGB star 47 Tuc V3 was observed using the GRAVITY instrument at ESO's Very Large Telescope Interferometer (VLTI) at 2-2.45 micron with a projected baseline length of up to 96 m. The object 47 Tuc V3 has been spatially resolved and stands as the first to attempt to spatially resolve an individual star in a globular cluster. The uniform-disk fit to the observed data results in an angular diameter of ~0.7 mas. Our modeling of the spectral energy distribution and near-infrared interferometric GRAVITY data suggests that the observed data can be explained by an optically thin dust shell with a 0.55 micron optical depth of 0.05-0.25, consisting of metallic iron grains, likely together with effects of the extended atmosphere of the central star. The dust temperature at the inner shell boundary is 500-800 K (corresponding to 23-90 stellar radii), significantly lower than observed in nearby oxygen-rich AGB stars. Radiation pressure on small (< 0.05 micron) iron grains is not sufficient to drive stellar winds. Therefore, iron grains may grow to larger sizes, even in the metal-poor environment. Alternatively, it is possible that the observed iron grain formation is a result of the mass outflow initiated by some other mechanism(s). The sensitivity and angular resolution of VLTI provides a new window onto spatially resolving individual stars in metal-poor globular clusters. This allows us to improve subsequent studies of the metallicity dependence of dust formation and mass loss., Comment: 7 pages, 3 figures, accepted for publication in Astronomy and Astrophysics

From Nov. 2019 to May 2020, the red supergiant star Betelgeuse experienced an unprecedented drop of brightness in the visible domain called the great dimming event. Large atmospheric dust clouds and large photospheric convective features are suspected to be responsible for it. To better understand the dimming event, we used mid-infrared long-baseline spectro-interferometric measurements of Betelgeuse taken with the VLTI/MATISSE instrument before (Dec. 2018), during (Feb. 2020), and after (Dec. 2020) the GDE. We present data in the 3.98 to 4.15\,$\mu$m range to cover SiO spectral features molecules as well as adjacent continuum. We have employed geometrical models, image reconstruction, as well as radiative transfer models to monitor the spatial distribution of SiO over the stellar surface. We find a strongly in-homogeneous spatial distribution of SiO that appears to be looking very different between our observing epochs, indicative of a vigorous activity in the stellar atmosphere. The contrast of our images is small in the pseudo-continuum for all epochs, implying that our MATISSE observations support both cold spot and dust cloud model.

Context. Rocky planets form by the concentration of solid particles in the inner few au regions of planet-forming disks. Their chemical composition reflects the materials in the disk available in the solid phase at the time the planets were forming. Aims. We aim to constrain the structure and dust composition of the inner disk of the young star HD 144432, using an extensive set of infrared interferometric data taken by the Very Large Telescope Interferometer (VLTI), combining PIONIER, GRAVITY, and MATISSE observations. Methods. We introduced a new physical disk model, TGMdust, to image the interferometric data, and to fit the disk structure and dust composition. We also performed equilibrium condensation calculations with GGchem. Results. Our best-fit model has three disk zones with ring-like structures at 0.15, 1.3, and 4.1 au. Assuming that the dark regions in the disk at ~0.9 au and at ~3 au are gaps opened by planets, we estimate the masses of the putative gap-opening planets to be around a Jupiter mass. We find evidence for an optically thin emission ($\tau<0.4$) from the inner two disk zones ($r<4$ au) at $\lambda>3\ \mu$m. Our silicate compositional fits confirm radial mineralogy gradients. To identify the dust component responsible for the infrared continuum emission, we explore two cases for the dust composition, one with a silicate+iron mixture and the other with a silicate+carbon one. We find that the iron-rich model provides a better fit to the spectral energy distribution. Conclusions. We propose that in the warm inner regions ($r<5$ au) of typical planet-forming disks, most if not all carbon is in the gas phase, while iron and iron sulfide grains are major constituents of the solid mixture along with forsterite and enstatite. Our analysis demonstrates the need for detailed studies of the dust in inner disks with new mid-infrared instruments such as MATISSE and JWST/MIRI., Comment: 29 pages, 24 figures

Changes in the flux and spectrum of Eta Carinae since 1900 have been attributed to the evolution of the central binary by some. Others suggest evolution in the occulting ejecta. The brightness jump in the 1940s, which coincided with the appearance of narrow forbidden emission lines, may have been caused by the clearing and ionization of intervening circumstellar ejecta. The brightening changed at a slower pace up through forty years later. Here we continue earlier studies focused on the long-term showing that the forbidden line emission increased in the early 1990s with no noticeable increase in the brightness of the Homunculus. We interpret that the increase in narrow line emission is due to decreased extinction in the LOS from the central binary to the Weigelt clumps. In 2000, the central stellar core increased in brightness at a faster rate without associated changes in the Homunculus. By 2018, hundreds of narrow-line absorptions from singly-ionized metals in our LOS from Eta Carinae disappeared, thought to be caused by increased ionization of metals. These three events (1990, 2000, and 2018) are explained by the dissipation of circumstellar material within the Homunculus close to the binary. Combining these changes with the steadiness of the Homunculus and the primary winds over the past four decades indicates that circumstellar ejecta in our direction have been cleared., Comment: 26 pages, 8 figures, to appear in ApJ

The Very Large Telescope Interferometer (VLTI) has been providing breakthrough images of the dust in the central parsecs of Active Galactic Nuclei (AGN), a key component of the AGN unification scheme and AGN host galaxy interaction. In single IR bands, the images can have multiple interpretations some of which could challenge the unification scheme. This is the case for the archetypal type 2 AGN of NGC1068. The ambiguity is reduced by multi-band temperature maps which are hindered by uncertainty in intra-band alignment. We create a chromatic model capable of simultaneously explaining the VLTI GRAVITY+MATISSE 2$\mu$m-13$\mu$m observations of the AGN in NGC1068. We use a simple disk and wind geometry populated with spherical black body emitters and dust obscuration to create a versatile multi-wavelength model for IR interferometric data of dusty objects. This simple geometry is capable of reproducing the K-N-band VLTI data, explains the complex single band images, and solves the alignment between bands. We find that the resulting geometry is consistent with previous studies. Compared to molecular gas emission, our model wind position angle (PA) of $22^3_2{\deg}$ is close to the mas scale outflowing CO(6-5) PA of ~33{\deg} seen with the ALMA. The equivalent 90{\deg} offset model disk PA is also consistent with the CO(6-5) disk axis of 112{\deg} as well as the mas scale disk axis from CO(2-1), CO(3-2), and HCO$^+$(4-3) of 115$\pm$5{\deg}. Furthermore, the resulting model visually resembles the equivalent achromatic image reconstructions. We conclude that the IR emitting structure surrounding the AGN can indeed be explained by the clumpy disk+wind iteration of the AGN unification scheme. Within the scheme, we find it is best explained as a type 2 and the obscuring dust chemistry is consistent with a mix of olivine silicates and 16$\pm$1% amorphous carbon., Comment: 15 pages, 19 appendix pages, 23 figures, and 2 tables. Accepted for publication in A&A

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

In this paper we present the first-ever $L$- and $M$-band interferometric observations of Circinus, building upon a recent $N$-band analysis. We used these observations to reconstruct images and fit Gaussian models to the $L$ and $M$ bands. Our findings reveal a thin edge-on disk whose width is marginally resolved and is the spectral continuation of the disk imaged in the $N$ band to shorter wavelengths. Additionally, we find a point-like source in the $L$ and $M$ bands that, based on the $LMN$-band spectral energy distribution fit, corresponds to the $N$-band point source. We also demonstrate that there is no trace of direct sightlines to hot dust surfaces in the circumnuclear dust structure of Circinus. By assuming the dust is present, we find that obscuration of A$_{\rm V} \gtrsim 250$ mag is necessary to reproduce the measured fluxes. Hence, the imaged disk could play the role of the obscuring "torus" in the unified scheme of active galactic nuclei. Furthermore, we explored the parameter space of the disk + hyperbolic cone radiative transfer models and identify a simple modification at the base of the cone. Adding a cluster of clumps just above the disk and inside the base of the hyperbolic cone provides a much better match to the observed temperature distribution in the central aperture. This aligns well with the radiation-driven fountain models that have recently emerged. Only the unique combination of sensitivity and spatial resolution of the VLTI allows such models to be scrutinized and constrained in detail. We plan to test the applicability of this detailed dust structure to other MATISSE-observed active galactic nuclei in the future., Comment: Main article: 12 pages, 6 figures. Accepted for publication in A&A. arXiv admin note: text overlap with arXiv:2205.01575

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

Previous STIS long-slit observations of Eta Carinae identified numerous absorption features in both the stellar spectrum, and in the adjacent nebular spectra, along our line-of-sight. The absorption features became temporarily stronger when the ionizing FUV radiation field was reduced by the periastron passage of the secondary star. Subsequently, dissipation of a dusty structure in our LOS has led to a long-term increase in the apparent magnitude of \ec, an increase in the ionizing UV radiation, and the disappearance of absorptions from multiple velocity-separated shells extending across the foreground Homunculus lobe. We use HST/STIS spectro-images, coupled with published infrared and radio observations, to locate this intervening dusty structure. Velocity and spatial information indicate the occulter is ~1000 au in front of Eta Carinae. The Homunculus is a transient structure composed of dusty, partially-ionized ejecta that eventually will disappear due to the relentless rain of ionizing radiation and wind from the current binary system along with dissipation and mixing with the ISM. This evolving complex continues to provide an astrophysical laboratory that changes on human timescales., Comment: 20 pages, 10 figures, submitted to ApJ

$Context.$ The YSO Th 28 possesses a highly collimated jet, which clearly exhibits an asymmetric brightness of its jet lobes at optical and NIR wavelengths. There may be asymmetry in the jet plasma parameters in opposite jet lobes (e.g. electron density, temperature, and outflow velocity). $Aims.$ We examined the Th 28 jet in a 3"x3" where the jet material is collimated and accelerated. Our goal is to map the morphology and determine its physical parameters to determine the physical origin of such asymmetries. $Methods.$ We present $JHK$-spectra of Th 28 obtained with the SINFONI on the (VLT, ESO) in June-July 2015. $Results.$ The [Fe II] emission originates in collimated jet lobes. Two new axial knots are detected at 1" in the blue lobe and 1".2 in the red lobe. The H$_2$ radiation is emitted from an extended region with a radius of $\gtrsim270$ au, which is perpendicular to the jet. The PV diagrams of the bright H$_2$ lines reveal faint H$_2$ emission along both jet lobes as well. The compact and faint H I emission (Pa$\beta$ and Br$\gamma$) comes from two regions, namely from a spherical region around the star and from the jet lobes. The size of the jet launching region is derived as 0".015 ($\sim$3 au at 185 pc), and the initial opening angle of the Th 28 jet is $\sim28^0$, which makes this jet substantially less collimated than most jets from other CTTs. $Conclusions.$ The emission in [Fe II], H$_2$, and H I lines suggests a morphology in which the ionised gas in the disc appears to be disrupted by the jet. The resolved disc-like H$_2$ emission most likely arises in the disc atmosphere from shocks caused by a radial uncollimated wind. The asymmetry of the [Fe II] photocentre shifts with respect to the jet source arises in the immediate vicinity of the driving source of Th28 and suggests that the observed brightness asymmetry is intrinsic as well., Comment: 18 pages, 14 figures

We use the long-slit spectroscopy mode of the VISIR-NEAR experiment to perform diffraction-limited observations of eight nearby Herbig Ae protoplanetary disks. We extract spectra for various locations along the slit with a spectral resolution of R = 300 and perform a compositional fit at each spatial location using spectral templates of silicates and the four PAH bands. This yields the intensity vs. location profiles of each species. Results. We could obtain spatially-resolved intensity profiles of the PAH emission features in the N-band for five objects (AB Aurigae, HD 97048, HD 100546, HD 163296, and HD 169142). We observe two kinds of PAH emission geometry in our sample: centrally-peaked (HD 97048) and ring-like (AB Aurigae, HD 100546, HD 163296, and potentially HD 169142). Comparing the spatial PAH emission profiles with near-infrared scattered light images, we find a strong correlation in the disk sub-structure but a difference in radial intensity decay rate. The PAH emission shows a less steep decline with distance from the star. Finally, we find a correlation between the presence of (sub-) micron-sized silicate grains leading to the depletion of PAH emission within the inner regions of the disks. In this work, we find the following: (1) PAH emission traces the extent of Herbig Ae disks to a considerable radial distance. (2) The correlation between silicate emission within the inner regions of disks and the depletion of PAH emission can result from dust-mixing and PAH coagulation mechanisms and competition over UV photons. (3) For all objects in our sample, PAHs undergo stochastic heating across the entire spatial extent of the disk and are not saturated. (4) The difference in radial intensity decay rates between the PAHs and scattered-light profiles may be attributed to shadowing and dust-settling effects, which affect the scattering grains more than the PAHs.

The binary eta Carinae is the closest example of a very massive star, which may have formed through a merger during its Great Eruption in the mid-nineteenth century. We aimed to confirm and improve the kinematics using a spectroscopic data set taken with the CTIO 1.5 m telescope over the time period of 2008-2020, covering three periastron passages of the highly eccentric orbit. We measure line variability of H-alpha and H-beta, where the radial velocity and orbital kinematics of the primary star were measured from the H-beta emission line using a bisector method. At phases away from periastron, we observed the He II 4686 emission moving opposite the primary star, consistent with a possible Wolf-Rayet companion, although with a seemingly narrow emission line. This could represent the first detection of emission from the companion., Comment: 11 pages, accepted to MNRAS