Source author record

T. Khouri

T. Khouri appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

22works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2025arXiv

ATOMIUM: Continuum emission and evidence of dust enhancement from binary motion

Low- and intermediate-mass stars on the asymptotic giant branch (AGB) account for a significant portion of the dust and chemical enrichment in their host galaxy. Here we present ALMA observations of the continuum emission at 1.24 mm around a sample of 17 stars from the ATOMIUM survey. From our analysis of the stellar contributions to the continuum flux, we find that the semi-regular variables all have smaller physical radii and fainter monochromatic luminosities than the Mira variables. Comparing these properties with pulsation periods, we find a positive trend between stellar radius and period only for the Mira variables with periods above 300 days and a positive trend between the period and the monochromatic luminosity only for the red supergiants and the most extreme AGB stars with periods above 500 days. We find that the continuum emission at 1.2 mm can be classified into four groups. "Featureless" continuum emission is confined to the (unresolved) regions close to the star for five stars in our sample, relatively uniform extended flux is seen for four stars, tentative elongated features are seen for three stars, and the remaining five stars have unique or unusual morphological features in their continuum maps. These features can be explained by binary companions to 10 out of the 14 AGB stars in our sample. Based on our results we conclude that there are two modes of dust formation: well established pulsation-enhanced dust formation and our newly proposed companion-enhanced dust formation. If the companion is located close to the AGB star, in the wind acceleration region, then additional dust formed in the wake of the companion can increase the mass lost through the dust driven wind. This explains the different dust morphologies seen around our stars and partly accounts for a large scatter in literature mass-loss rates, especially among semiregular stars with small pulsation periods.

preprint2025arXiv

X TrA through the eyes of MATISSE: More evidence of clumpy molecular layers around C-type asymptotic giant branch stars

Aims. The goal of this study is to further the understanding of the wind formation mechanism in asymptotic giant branch (AGB) stars through the analysis of the close environment (within a few stellar radii) of the carbon star X TrA. Methods. X TrA was observed for the first time with the Mid-Infrared SpectroScopic Experiment instrument (MATISSE) in the L and N bands in low spectral resolution mode (R=30), and its close surroundings were mapped in specific wavelength ranges corresponding to specific molecules ($C_2H_2$ and HCN, at 3.1 and 3.8 $μ$m) and dust (amorphous carbon and, for example, Sic at 11.3 $μ$m), via image reconstruction techniques. Results. The angular diameter of the star ranges from 10 mas in the L band pseudo-continuum (3.5 $μ$m) to 20 mas at 3.1 and 11.3 $μ$m. The reconstructed images show some mild elongated features (along the east-west direction) and asymmetric protrusions, which are most evident around 3.1 $μ$m. Imaging results highlight the clumpy nature of the circumstellar environment, starting from the photospheric region up to more distant layers. Conclusions. The angular diameters found for X TrA in the image data are in agreement with previous photospheric diameter estimates (following VLTI/MIDI 8-13 $μ$m observations), and their wavelength dependence is similar to values found for other carbon stars observed with MATISSE (R Scl and V Hya). The 3.1 $μ$m images presented here show highly asymmetric features, another case of a C-rich star with irregular morphologies close to the stellar disk; this supports the notion that the $C_2H_2+HCN$ abundance distribution usually originates from a clumpy layer around carbon stars.

preprint2022arXiv

First detection of AlF line emission towards M-type AGB stars

The nucleosynthesis production of fluorine (F) is still a matter of debate. Asymptotic giant branch (AGB) stars are one of the main candidates for F production. However, their contribution to the total F budget is not fully known due to the lack of observations. In this paper, we report the detection of AlF line emission, one of the two main carriers of F in the gas-phase in the outflow of evolved stars, towards five nearby oxygen-rich AGB stars, $o$ Ceti, R Leo, IK Tau, R Dor, and W Hya. From spatially resolved observations, we estimated the AlF emitting region with a radius $\sim11R_{\star}$ for $o$ Ceti and $\sim9R_{\star}$ for R Leo. From population diagram analysis, we report the AlF column densities of $\sim 5.8\times10^{15}$ cm$^{-2}$ and $\sim 3\times10^{15}$ cm$^{-2}$ for $o$ Ceti and R Leo, respectively, within these regions. For $o$ Ceti, we used the C$^{18}$O ($v=0$, $J=3-2$) observations to estimate the H$_2$ column density of the emitting region. We found a fractional abundance of $f_{\rm AlF/H_2}\sim(2.5\pm1.7)\times10^{-8}$. This gives a lower limit on the F budget in $o$ Ceti and is compatible with the solar F budget $f_{\rm F/H_2}=(5\pm2)\times10^{-8}$. For R Leo, a fractional abundance $f_{\rm AlF/H_2}=(1.2\pm0.5)\times10^{-8}$ is estimated. For other sources, we cannot precisely determine the emitting region based on the available data. Assuming an emitting region with a radius of $\sim 11R_{\star}$ and the rotational temperatures derived for $o$ Ceti and R Leo, we crudely approximated the AlF column density to be $\sim(1.2-1.5)\times10^{15}$ cm$^{-2}$ in W Hya, $\sim(2.5-3.0)\times10^{14}$ cm$^{-2}$ in R Dor, and $\sim(0.6-1.0)\times10^{16}$ cm$^{-2}$ in IK Tau. These result in fractional abundances within a range of $f_{\rm AlF/H_2}\sim(0.1-4)\times10^{-8}$ in W Hya, R Dor, and IK Tau.

preprint2021arXiv

ATOMIUM: ALMA tracing the origins of molecules in dust forming oxygen rich M-type stars: Motivation, sample, calibration, and initial results

This overview paper presents ATOMIUM, a Large Programme in Cycle 6 with the Atacama Large Millimeter-submillimeter Array (ALMA). The goal of ATOMIUM is to understand the dynamics and the gas phase and dust formation chemistry in the winds of evolved asymptotic giant branch (AGB) and red supergiant (RSG) stars. A more general aim is to identify chemical processes applicable to other astrophysical environments. 17 oxygen-rich AGB and RSG stars spanning a range in (circum)stellar parameters and evolutionary phases were observed in a homogeneous observing strategy allowing for an unambiguous comparison. Data were obtained between 213.83 and 269.71 GHz at high (0.025-0.050 arcsec), medium (0.13-0.24 arcsec), and low (about 1 arcsec) angular resolution. The sensitivity per 1.3 km/s channel was 1.5-5 mJy/beam. 13 molecules were designated as primary molecules in the survey: CO, SiO, AlO, AlOH, TiO, TiO2, HCN, SO, SO2, SiS, CS, H2O, and NaCl. The scientific motivation, survey design, sample properties, data reduction, and an overview of the data products are described; and we highlight one scientific result - the wind kinematics of the ATOMIUM sources. The ATOMIUM sources often have a slow wind acceleration, and a fraction of the gas reaches a velocity which can be up to a factor of two times larger than previously reported terminal velocities assuming isotropic expansion, and the wind kinematic profiles establish that the radial velocity described by the momentum equation for a spherical wind structure cannot capture the complexity of the velocity field. In 15 sources, some molecular transitions other than 12CO v=0 J=2-1 reach a higher outflow velocity, with a spatial emission zone that is often greater than 30 stellar radii, but much less than the extent of CO. Binary interaction with a (sub)stellar companion might (partly) explain the non-monotonic behaviour of the projected velocity field.

preprint2020arXiv

CO and HCN isotopologue ratios in the outflows of AGB stars

Isotopologue line intensity ratios of circumstellar molecules have been widely used to trace the photospheric elemental isotopic ratios of evolved stars. However, depending on the molecular species and the physical conditions of the environment, the circumstellar isotopologue ratio may deviate considerably from the stellar atmospheric value. In this paper, we aim to examine how the CO and HCN abundance ratios vary radially due to chemical reactions in the outflows of AGB stars and the effect of excitation and optical depth on the resulting line intensity ratios. We find that the circumstellar 12CO/13CO can deviate from its atmospheric value by up to 25-94% and 6-60% for C- and O-type CSEs, respectively. We show that variations of the intensity of the ISRF and the gas kinetic temperature can significantly influence the CO isotopologue ratio in the outer CSEs. On the contrary, the H12CN/H13CN ratio is stable for all tested mass-loss rates. The RT modeling shows that the integrated line intensity ratio of CO of different rotational transitions varies significantly for stars with intermediate mass-loss rates due to combined chemical and excitation effects. In contrast, the excitation conditions for the both HCN isotopologues are the same. We demonstrate the importance of using the isotopologue abundance profiles from chemical models as inputs to RT models in the interpretation of isotopologue observations. Previous studies of CO isotopologue ratios are based on multi-transition data for individual sources and it is difficult to estimate the errors in the reported values due to assumptions that are not entirely correct according to this study. If anything, previous studies may have overestimated the circumstellar 12CO/13CO abundance ratio. The use of the HCN as a tracer of C isotope ratios is affected by fewer complicating problems, provided one accounts corrections for high optical depths.

preprint2020arXiv

DEATHSTAR: Nearby AGB stars with the Atacama Compact Array I. CO envelope sizes and asymmetries: A new hope for accurate mass-loss-rate estimates

This is the first publication of the DEATHSTAR project. The goal of the project is to reduce the uncertainties of observational estimates of mass-loss rates from Asymptotic Giant Branch (AGB) stars. Line emission from 12CO J=2-1 and 3-2 were mapped using the ACA. In this initial analysis, the emission distribution was fit to a Gaussian distribution in the uv-plane. Detailed radiative transfer analysis will be presented in the future. The axes of the best-fit Gaussian at the line center of the 12CO J=2-1 emission gives a first indication of the size of the emitting region. Furthermore, the fitting results, such as the major and minor axis, center position, and the goodness of fit across both lines, constrain the symmetry of the emission distribution. We find that the CO envelope sizes are, in general, larger for C-type than for M-type AGB stars, which is expected if the CO/H2 ratio is larger in C-type stars. Furthermore, a relation between the 12CO J=2-1 size and circumstellar density is shown that, while in broad agreement with photodissociation calculations, reveals large scatter and systematic differences between the stellar types. The majority of the sources have CO envelopes that are consistent with a spherically symmetric, smooth outflow. For about a third of the sources, indications of strong asymmetries are found. This is consistent with previous interferometric investigations of northern sources. Smaller scale asymmetries are found in a larger fraction of sources. These results for CO envelope radii and shapes can be used to constrain detailed radiative transfer modeling of the same stars so as to determine mass-loss rates that are independent of photodissociation models. For a large fraction of the sources, observations at higher spatial resolution will be necessary to further investigate the complex circumstellar dynamics revealed by our ACA observations.

preprint2020arXiv

Determination of the Fermi surface and field-induced quasi-particle tunneling around the Dirac nodal-loop in ZrSiS

Unambiguous and complete determination of the Fermi surface is a primary step in understanding the electronic properties of topical metals and semi-metals, but only in a relatively few cases has this goal been realized. In this work, we present a systematic high-field quantum oscillation study up to 35 T on ZrSiS, a textbook example of a nodal-line semimetal with only linearly dispersive bands crossing the Fermi energy. The topology of the Fermi surface is determined with unprecedented precision and all pockets are identified by comparing the measured angle dependence of the quantum oscillations to density functional theory calculations. Comparison of the Shubnikov-de Haas and de Haas-van Alphen oscillations at low temperatures and analysis of the respective Dingle plots reveal the presence of significantly enhanced scattering on the electron pocket. Above a threshold field that is aligned along the c-axis of the crystal, the specific cage-like Fermi surface of ZrSiS allows for electron-hole tunneling to occur across finite gaps in momentum space leading to quantum oscillations with a complex frequency spectrum. Additional high-frequency quantum oscillations signify magnetic breakdown orbits that encircle the entire Dirac nodal loop. We suggest that the persistence of quantum oscillations in the resistivity to high temperatures is caused by Stark interference between orbits of nearly equal masses.

preprint2020arXiv

Inner dusty envelope of the AGB stars W~Hydrae, SW\,Virginis, and R~Crateris using SPHERE/ZIMPOL

Stars with initial masses between $\sim0.8$ and 8~$M_\odot$ present copious mass loss during the asymptotic giant branch (AGB) at the end of their lives. The accepted mass-loss mechanism requires radiation pressure acting on dust grains that form in the extended AGB stellar atmospheres. The details of this process are not yet well understood, however. Using the extreme-adaptive-optics imager and polarimeter SPHERE/ZIMPOL, we observed light polarised by grains around W\,Hya, SW\,Vir, and R\,Crt, which have mass-loss rates between 10$^{-7}$ and 10$^{-6}~M_\odot~{\rm yr^{-1}}$. We find the distribution of dust to be asymmetric around the three targets. A biconical morphology is seen for R Crt, with a position angle that is very similar to those inferred from interferometric observations of maser emission and of mid-infrared continuum emission. The cause of the biconical outflow cannot be directly inferred from the ZIMPOL data. The dust grains polarise light more efficiently at 0.65~$μ$m for R\,Crt and SW\,Vir and at 0.82~$μ$m for W\,Hya. This indicates that at the time of the observations, the grains around SW\,Vir and R\,Crt had sizes $< 0.1~μ$m, while those around W\,Hya were larger, with sizes $\gtrsim 0.1~μ$m. The asymmetric distribution of dust around R\,Crt makes the interpretation more uncertain for this star, however. We find that polarised light is produced already from within the visible photosphere of W~Hya, which we reproduce using models with an inner dust shell that is optically thick to scattering. The radial profile of the polarised light observed around W\,Hya reveal a steep dust density profile. We find the wind-acceleration region of W\,Hya to extend to at least $\sim 7~R_\star$, in agreement with theoretical predictions of acceleration up to $\sim 12~R_\star$.

preprint2019arXiv

Unusual phase boundary of the magnetic-field-tuned valence transition in CeOs$_4$Sb$_{12}$

The phase diagram of the filled skutterudite CeOs$_4$Sb$_{12}$ has been mapped in fields $H$ of up to 60 T and temperatures $T$ down to 0.5 K using resistivity, magnetostriction, and MHz conductivity. The valence transition separating the semimetallic, low-$H$, low-$T$, $\cal{L}$ phase from the metallic high-$H$, high-$T$ $\cal{H}$ phase exhibits a very unusual, wedge-shaped phase boundary, with a non-monotonic gradient alternating between positive and negative. This is quite different from the text-book "elliptical" phase boundary usually followed by valence transitions. Analysis of Shubnikov-de Haas oscillations within the $\cal{H}$ phase reveals an effective mass that increases as $H$ drops toward the $\cal{H-L}$ phase boundary, suggesting proximity to a quantum-critical point. The associated magnetic fluctuations may be responsible for the anomalous $H,T$ dependence of the valence transition at high $H$, whereas the low$-H$, high$-T$ portion of the phase boundary may rather be associated with the proximity of CeOs$_4$Sb$_{12}$ to a topological semimetal phase induced by uniaxial stress.

preprint2016arXiv

ALMA observations of the vibrationally-excited rotational CO transition $v=1, J=3-2$ towards five AGB stars

We report the serendipitous detection with ALMA of the vibrationally-excited pure-rotational CO transition $v=1, J=3-2$ towards five asymptotic giant branch (AGB) stars, $o$ Cet, R Aqr, R Scl, W Aql, and $π^1$ Gru. The observed lines are formed in the poorly-understood region located between the stellar surface and the region where the wind starts, the so-called warm molecular layer. We successfully reproduce the observed lines profiles using a simple model. We constrain the extents, densities, and kinematics of the region where the lines are produced. R Aqr and R Scl show inverse P-Cygni line profiles which indicate infall of material onto the stars. The line profiles of $o$ Cet and R Scl show variability. The serendipitous detection towards these five sources shows that vibrationally-excited rotational lines can be observed towards a large number of nearby AGB stars using ALMA. This opens a new possibility for the study of the innermost regions of AGB circumstellar envelopes.

preprint2016arXiv

Anisotropic and strong negative magneto-resistance in the three-dimensional topological insulator Bi2Se3

We report on high-field angle-dependent magneto-transport measurements on epitaxial thin films of Bi2Se3, a three-dimensional topological insulator. At low temperature, we observe quantum oscillations that demonstrate the simultaneous presence of bulk and surface carriers. The magneto- resistance of Bi2Se3 is found to be highly anisotropic. In the presence of a parallel electric and magnetic field, we observe a strong negative longitudinal magneto-resistance that has been consid- ered as a smoking-gun for the presence of chiral fermions in a certain class of semi-metals due to the so-called axial anomaly. Its observation in a three-dimensional topological insulator implies that the axial anomaly may be in fact a far more generic phenomenon than originally thought.

preprint2016arXiv

Constraints on the H2O formation mechanism in the wind of carbon-rich AGB stars

Context. The recent detection of warm H$_2$O vapor emission from the outflows of carbon-rich asymptotic giant branch (AGB) stars challenges the current understanding of circumstellar chemistry. Two mechanisms have been invoked to explain warm H$_2$O vapor formation. In the first, periodic shocks passing through the medium immediately above the stellar surface lead to H$_2$O formation. In the second, penetration of ultraviolet interstellar radiation through a clumpy circumstellar medium leads to the formation of H$_2$O molecules in the intermediate wind. Aims. We aim to determine the properties of H$_2$O emission for a sample of 18 carbon-rich AGB stars and subsequently constrain which of the above mechanisms provides the most likely warm H$_2$O formation pathway. Methods, Results, and Conclusions. See paper.

preprint2016arXiv

High-temperature quantum Hall effect in finite gapped HgTe quantum wells

We report on the observation of the quantum Hall effect at high temperatures in HgTe quantum wells with a finite band gap and a thickness below and above the critical thickness $d_\textnormal{c}$ that separates a conventional semiconductor from a two-dimensional topological insulator. At high carrier concentrations we observe a quantized Hall conductivity up to 60\,K with energy gaps between Landau Levels of the order of 25\,meV, in good agreement with the Landau Level spectrum obtained from $\mathbf{k\cdot p}$-calculations. Using the scaling approach for the plateau-plateau transition at $ν=2\rightarrow 1$, we find the scaling coefficient $κ=0.45 \pm 0.04$ to be consistent with the universality of scaling theory and we do not find signs of increased electron-phonon interaction to alter the scaling even at these elevated temperatures. Comparing the high temperature limit of the quantized Hall resistance in HgTe quantum wells with a finite band gap with room temperature experiment in graphene, we find the energy gaps at the break-down of the quantization to exceed the thermal energy by the same order of magnitude.

preprint2016arXiv

Linear magnetoresistance in a quasi-free two dimensional electron gas in an ultra-high mobility GaAs quantum well

We report a magnetotransport study of an ultra-high mobility ($\barμ\approx 25\times 10^6$\,cm$^2$\,V$^{-1}$\,s$^{-1}$) $n$-type GaAs quantum well up to 33 T. A strong linear magnetoresistance (LMR) of the order of 10$^5$ % is observed in a wide temperature range between 0.3 K and 60 K. The simplicity of our material system with a single sub-band occupation and free electron dispersion rules out most complicated mechanisms that could give rise to the observed LMR. At low temperature, quantum oscillations are superimposed onto the LMR. Both, the featureless LMR at high $T$ and the quantum oscillations at low $T$ follow the empirical resistance rule which states that the longitudinal conductance is directly related to the derivative of the transversal (Hall) conductance multiplied by the magnetic field and a constant factor $α$ that remains unchanged over the entire temperature range. Only at low temperatures, small deviations from this resistance rule are observed beyond $ν=1$ that likely originate from a different transport mechanism for the composite fermions.

preprint2016arXiv

Study of the inner dust envelope and stellar photosphere of the AGB star R Doradus using SPHERE/ZIMPOL

We use high-angular-resolution images obtained with SPHERE/ZIMPOL to study the photosphere, the warm molecular layer, and the inner wind of the close-by oxygen-rich AGB star R Doradus. We present observations in filters V, cntH$α$, and cnt820 and investigate the surface brightness distribution of the star and of the polarised light produced in the inner envelope. Thanks to second-epoch observations in cntH$α$, we are able to see variability on the stellar photosphere. We find that in the first epoch the surface brightness of R Dor is asymmetric in V and cntH$α$, the filters where molecular opacity is stronger, while in cnt820 the surface brightness is closer to being axisymmetric. The second-epoch observations in cntH$α$ show that the morphology of R Dor changes completely in a timespan of 48 days to a more axisymmetric and compact configuration. The polarised intensity is asymmetric in all epochs and varies by between a factor of 2.3 and 3.7 with azimuth for the different images. We fit the radial profile of the polarised intensity using a spherically symmetric model and a parametric description of the dust density profile, $ρ(r)=ρ_\circ r^{-n}$. On average, we find exponents of $- 4.5 \pm 0.5$ that correspond to a much steeper density profile than that of a wind expanding at constant velocity. The dust densities we derive imply an upper limit for the dust-to-gas ratio of $\sim 2\times10^{-4}$ at 5.0 $R_\star$. Given the uncertainties in observations and models, this value is consistent with the minimum values required by wind-driving models for the onset of a wind, of $\sim 3.3\times10^{-4}$. However, if the steep density profile we find extends to larger distances from the star, the dust-to-gas ratio will quickly become too small for the wind of R Dor to be driven by the grains that produce the scattered light.

preprint2015arXiv

Dusty wind of W Hya. Multi-wavelength modelling of the present-day and recent mass-loss

Low- and intermediate-mass stars go through a period of intense mass-loss at the end of their lives in a phase known as the asymptotic giant branch (AGB). During the AGB a significant fraction of their initial mass is expelled in a stellar wind. This process controls the final stages of their evolution and contributes to the chemical evolution of galaxies. However, the wind-driving mechanism of AGB stars is not yet well understood, especially so for oxygen-rich sources. Characterizing both the present-day mass-loss and wind structure and the evolution of the mass-loss rate of such stars is paramount to advancing our understanding of this processes. We modelled the dust envelope of W Hya using an advanced radiative transfer code. The dust model was analysed in the light of a previously calculated gas-phase wind model and compared to measurements available in the literature, such as infrared spectra, infrared images, and optical scattered light fractions. We find that the dust spectrum of W Hya can partly be explained by a gravitationally bound dust shell that probably is responsible for most of the amorphous Al$_2$O$_3$ emission. The composition of the large ($\sim$\,0.3\,$μ$m) grains needed to explain the scattered light cannot be constrained, but probably is dominated by silicates. Silicate emission in the thermal infrared was found to originate from beyond 40 AU from the star and we find that they need to have substantial near-infrared opacities to be visible at such large distances. The increase in near-infrared opacity of the dust at these distances roughly coincides with a sudden increase in expansion velocity as deduced from the gas-phase CO lines. Finally, the recent mass loss of W Hya is confirmed to be highly variable and we identify a strong peak in the mass-loss rate that occurred about 3500 years ago and lasted for a few hundred years.

preprint2015arXiv

H2O vapor excitation in dusty AGB envelopes. A PACS view of OH 127.8+0.0

Context: AGB stars lose a large percentage of their mass in a dust-driven wind. This creates a circumstellar envelope, which can be studied through thermal dust emission and molecular emission lines. In the case of high mass-loss rates, this study is complicated by the high optical depths and the intricate coupling between gas and dust radiative transfer characteristics. An important aspect of the physics of gas-dust interactions is the strong influence of dust on the excitation of several molecules, including H2O. Aims: The dust and gas content of the envelope surrounding the high mass-loss rate OH/IR star OH 127.8+0.0, as traced by Herschel observations, is studied, with a focus on the H2O content and the dust-to-gas ratio. We report detecting a large number of H2O vapor emission lines up to J = 9 in the Herschel data, for which we present the measured line strengths. Methods: The treatments of both gas and dust species are combined using two numerical radiative transfer codes. The method is illustrated for both low and high mass-loss-rate sources. Specifically, we discuss different ways of assessing the dust-to-gas ratio: 1) from the dust thermal emission spectrum and the CO molecular gas line strengths; 2) from the momentum transfer from dust to gas and the measured gas terminal velocity; and 3) from the determination of the required amount of dust to reproduce H2O lines for a given H2O vapor abundance. These three diagnostics probe different zones of the outflow, for the first time allowing an investigation of a possible radial dependence of the dust-to-gas ratio. Results: ... Continued in article.

preprint2014arXiv

ALMA Observations of Anisotropic Dust Mass-loss in the Inner Circumstellar Environment of the Red Supergiant VY Canis Majoris

The processes leading to dust formation and the subsequent role it plays in driving mass loss in cool evolved stars is an area of intense study. Here we present high resolution ALMA Science Verification data of the continuum emission around the highly evolved oxygen-rich red supergiant VY CMa. These data enable us to study the dust in its inner circumstellar environment at a spatial resolution of 129 mas at 321 GHz and 59 mas at 658 GHz, thus allowing us to trace dust on spatial scales down to 11 R$_{\star}$ (71 AU). Two prominent dust components are detected and resolved. The brightest dust component, C, is located 334 mas (61 R$_{\star}$) South East of the star and has a dust mass of at least $2.5\times 10^{-4}$ M$_{\odot}$. It has a dust emissivity spectral index of $β=-0.1$ at its peak, implying that it is optically thick at these frequencies with a cool core of $T_{d}\lesssim 100$ K. Interestingly, not a single molecule in the ALMA data has emission close to the peak of this massive dust clump. The other main dust component, VY, is located at the position of the star and contains a total dust mass of $4.0 \times 10^{-5} $M$_{\odot}$. It also contains a weaker dust feature extending over $60$ R$_{\star}$ to the North with the total component having a typical dust emissivity spectral index of $β=0.7$. We find that at least $17\%$ of the dust mass around VY CMa is located in clumps ejected within a more quiescent roughly spherical stellar wind, with a quiescent dust mass loss rate of $5 \times 10^{-6}$ M$_{\odot} $yr$^{-1}$. The anisotropic morphology of the dust indicates a continuous, directed mass loss over a few decades, suggesting that this mass loss cannot be driven by large convection cells alone.

preprint2014arXiv

The wind of W Hya as seen by Herschel - I. The CO envelope

Asymptotic giant branch (AGB) stars lose their envelopes by means of a stellar wind whose driving mechanism is not understood well. Characterizing the composition and thermal and dynamical structure of the outflow provides constraints that are essential for understanding AGB evolution, including the rate of mass loss and isotopic ratios. We modeled the CO emission from the wind of the low mass-loss rate oxygen-rich AGB star W Hya using data obtained by the HIFI, PACS, and SPIRE instruments onboard the Herschel Space Observatory and ground-based telescopes. $^{12}$CO and $^{13}$CO lines are used to constrain the intrinsic $^{12}$C/$^{13}$C ratio from resolved HIFI lines. The acceleration of the outflow up to about 5.5 km/s is quite slow and can be represented by a beta-type velocity law with index 5. Beyond this point, acceleration up the terminal velocity of 7 km/s is faster. Using the J=10-9, 9-8, and 6-5 transitions, we find an intrinsic $^{12}$C/$^{13}$C ratio of $18\pm10$ for W Hya, where the error bar is mostly due to uncertainties in the $^{12}$CO abundance and the stellar flux around 4.6 $μ$m. To match the low-excitation CO lines, these molecules need to be photo-dissociated at about 500 stellar radii. The radial dust emission intensity profile measured by PACS images at 70 $μ$m shows substantially stronger emission than our model predicts beyond 20 arcsec. The initial slow acceleration of the wind implies inefficient wind driving in the lower part of the envelope. The final injection of momentum in the wind might be the result of an increase in the opacity thanks to the late condensation of dust species. The derived intrinsic isotopologue ratio for W Hya is consistent with values set by the first dredge-up and suggestive of an initial mass of 2 M$_\odot$ or more. However, the uncertainty in the main-sequence mass derived based on this isotopologic ratio is large.

preprint2014arXiv

The wind of W Hya as seen by Herschel. II. The molecular envelope of W Hya

The evolution of low- and intermediate-mass stars on the asymptotic giant branch (AGB) is mainly controlled by the rate at which these stars lose mass in a stellar wind. Understanding the driving mechanism and strength of the stellar winds of AGB stars and the processes enriching their surfaces with products of nucleosynthesis are paramount to constraining AGB evolution and predicting the chemical evolution of galaxies. In a previous paper we have constrained the structure of the outflowing envelope of W Hya using spectral lines of the $^{12}$CO molecule. Here we broaden this study by modelling an extensive set of H$_{2}$O and $^{28}$SiO lines observed by the three instruments on board Herschel using a state-of-the-art molecular excitation and radiative transfer code. The oxygen isotopic ratios and the $^{28}$SiO abundance profile can be connected to the initial stellar mass and to crucial aspects of dust formation at the base of the stellar wind, respectively. The modelling of H$_{2}$O and $^{28}$SiO confirms the properties of the envelope model of W Hya derived from $^{12}$CO lines. We find an H$_2$O ortho-to-para ratio of 2.5\,$^{+2.5}_{-1.0}$, consistent with what is expected for an AGB wind. The O$^{16}$/O$^{17}$ ratio indicates that W Hya has an initial mass of about 1.5 M$_\odot$. Although the ortho- and para-H$_{2}$O lines observed by HIFI appear to trace gas of slightly different physical properties, a turbulence velocity of $0.7\pm0.1$ km s$^{-1}$ fits the HIFI lines of both spin isomers and those of $^{28}$SiO well. The ortho- and para-H$_2^{16}$O and $^{28}$SiO abundances relative to H$_{2}$ are $(6^{+3}_{-2}) \times 10^{-4}$, $(3^{+2}_{-1}) \times 10^{-4}$, and $(3.3\pm 0.8)\times 10^{-5}$, respectively. Assuming a solar silicon-to-carbon ratio, the $^{28}$SiO line emission model is consistent with about one-third of the silicon atoms being locked up in dust particles.

preprint2011arXiv

Herschel/HIFI observations of O-rich AGB stars : molecular inventory

Spectra, taken with the heterodyne instrument, HIFI, aboard the Herschel Space Observatory, of O-rich asymptotic giant branch (AGB) stars which form part of the guaranteed time key program HIFISTARS are presented. The aim of this program is to study the dynamical structure, mass-loss driving mechanism, and chemistry of the outflows from AGB stars as a function of chemical composition and initial mass. We used the HIFI instrument to observe nine AGB stars, mainly in the H2O and high rotational CO lines We investigate the correlation between line luminosity, line ratio and mass-loss rate, line width and excitation energy. A total of nine different molecules, along with some of their isotopologues have been identified, covering a wide range of excitation temperature. Maser emission is detected in both the ortho- and para-H2O molecules. The line luminosities of ground state lines of ortho- and para-H2O, the high-J CO and NH3 lines show a clear correlation with mass-loss rate. The line ratios of H2O and NH3 relative to CO J=6-5 correlate with the mass-loss rate while ratios of higher CO lines to the 6-5 is independent of it. In most cases, the expansion velocity derived from the observed line width of highly excited transitions formed relatively close to the stellar photosphere is lower than that of lower excitation transitions, formed farther out, pointing to an accelerated outflow. In some objects, the vibrationally excited H2O and SiO which probe the acceleration zone suggests the wind reaches its terminal velocity already in the innermost part of the envelope, i.e., the acceleration is rapid. Interestingly, for R Dor we find indications of a deceleration of the outflow in the region where the material has already escaped from the star.

preprint2010arXiv

A high resolution line survey of IRC+10216 with Herschel. First results: Detection of warm silicon dicarbide SiC2

We present the first results of a high-spectral-resolution survey of the carbon-rich evolved star IRC+10216 that was carried out with the HIFI spectrometer onboard Herschel. This survey covers all HIFI bands, with a spectral range from 488 to 1901GHz. In this letter we focus on the band-1b spectrum, in a spectral range 554.5-636.5GHz, where we identified 130 spectral features with intensities above 0.03 K and a signal-to-noise ratio >5. Detected lines arise from HCN, SiO, SiS, CS, CO, metal-bearing species and, surprisingly, silicon dicarbide (SiC2). We identified 55 SiC2 transitions involving energy levels between 300 and 900 K. By analysing these rotational lines, we conclude that SiC2 is produced in the inner dust formation zone, with an abundance of ~2x10^-7 relative to molecular hydrogen. These SiC2 lines have been observed for the first time in space and have been used to derive an SiC2 rotational temperature of ~204 K and a source-averaged column density of ~6.4x10^15 cm^-2. Furthermore, the high quality of the HIFI data set was used to improve the spectroscopic rotational constants of SiC2.