BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Leen Decin

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

ResearcherUnclaimed source record
astro-ph.SRastro-ph.EPastro-ph.GAastro-ph.IM

Catalog footprint

What is connected

24works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

24 published item(s)

preprint2022arXiv

Grid of pseudo-2D chemistry models for tidally locked exoplanets -- II. The role of photochemistry

Photochemistry is expected to change the chemical composition of the upper atmospheres of irradiated exoplanets through the dissociation of species, such as methane and ammonia, and the association of others, such as hydrogen cyanide. Although primarily the high altitude day side should be affected by photochemistry, it is still unclear how dynamical processes transport photochemical species throughout the atmosphere, and how these chemical disequilibrium effects scale with different parameters. In this work we investigate the influence of photochemistry in a two-dimensional context, by synthesizing a grid of photochemical models across a large range of temperatures. We find that photochemistry can strongly change the atmospheric composition, even up to depths of several bar in cool exoplanets. We further identify a sweet spot for the photochemical production of hydrogen cyanide and acetylene, two important haze precursors, between effective temperatures of 800 and 1400 K. The night sides of most cool planets (effective temperature < 1800 K) are shown to host photochemistry products, transported from the day side by horizontal advection. Synthetic transmission spectra are only marginally affected by photochemistry, but we suggest that observational studies probing higher altitudes, such as high-resolution spectroscopy, take photochemistry into account.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Grid of Pseudo-2D Chemistry Models for Tidally-Locked Exoplanets. I. The Role of Vertical and Horizontal Mixing

The atmospheres of synchronously rotating exoplanets are intrinsically three-dimensional, and fast vertical and horizontal winds are expected to mix the atmosphere, driving the chemical composition out of equilibrium. Due to the longer computation times associated with multi-dimensional forward models, horizontal mixing has only been investigated for a few case studies. In this paper, we aim to generalize the impact of horizontal and vertical mixing on the chemistry of exoplanet atmospheres over a large parameter space. We do this by applying a sequence of post-processed forward models for a large grid of synchronously rotating gaseous exoplanets, where we vary the effective temperature (between 400 K and 2600 K), surface gravity, and rotation rate. We find that there is a dichotomy in the horizontal homogeneity of the chemical abundances. Planets with effective temperatures below 1400 K tend to have horizontally homogeneous, vertically quenched chemical compositions, while planets hotter than 1400 K exhibit large compositional day-night differences for molecules such as methane. Furthermore, we find that the planet's rotation rate impacts the planetary climate, and thus also the molecular abundances and transmission spectrum. By employing a hierarchical modelling approach, we assess the relative importance of disequilibrium chemistry on the exoplanet transmission spectrum, and conclude that the temperature has the most profound impact. Temperature differences are also the main cause of limb asymmetries, which we estimate could be observable with the James Webb Space Telescope. This work highlights the value of applying a consistent modelling setup to a broad parameter space in exploratory theoretical research.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Bottom-up dust nucleation theory in oxygen-rich evolved stars I. Aluminium oxide clusters

Aluminum oxide (alumina, Al$_{2}$O$_{3}$) is a promising candidate as a primary dust condensate in the atmospheres of oxygen-rich evolved stars. Therefore, alumina \textit{seed} particles might trigger the onset of stellar dust formation and of stellar mass loss in the wind. However, the formation of alumina dust grains is not well understood.} {To shed light on the initial steps of cosmic dust formation (i.e. nucleation) in oxygen-rich environments by a quantum-chemical bottom-up approach.} {Starting with an elemental gas-phase composition, we construct a detailed chemical-kinetic network describing the formation and destruction of aluminium-bearing molecules and dust-forming (Al$_{2}$O$_{3}$)$_n$ clusters up to the size of dimers ($n$=2) coagulating to tetramers ($n=$4). Intermediary species include the prevalent gas-phase molecules AlO and AlOH, and Al$_x$O$_y$ clusters with $x=$1$-$5, $y=$1$-$6. The resulting extensive network is applied to two model stars, representing a semi-regular variable and a Mira-type, and to different circumstellar gas trajectories including a non-pulsating outflow and a pulsating model. The growth of larger-sized (Al$_{2}$O$_{3}$)$_n$ clusters with $n=$4$-$10 is described by the temperature-dependent Gibbs free energies of the most favourable structures (i.e. the global minima clusters) as derived from global optimisation techniques and calculated by density functional theory. We provide energies, bond characteristics, electrostatic properties and vibrational spectra of the clusters as a function of size $n$ and compare these to corundum corresponding to the crystalline bulk limit ($n\rightarrow \infty$). {The circumstellar aluminium gas-phase chemistry in oxygen-rich giants is primarily controlled by AlOH and AlO, which are tightly coupled by the reactions AlO+H$_2$, AlO+H$_2$O, and their reverse. Models of ...}

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Atomium: A high-resolution view on the highly asymmetric wind of the AGB star Pi1 Gruis. I. First detection of a new companion and its effect on the inner wind

The nebular circumstellar environments of cool evolved stars are known to harbour a rich morphological complexity of gaseous structures on different length scales. A large part of these density structures are thought to be brought about by the interaction of the stellar wind with a close companion. The S-type asymptotic giant branch star Pi1 Gruis, which has a known companion at ~440 au and is thought to harbour a second, closer-by (<10 au) companion, was observed with the Atacama Large Millimeter/submillimeter Array as part of the ATOMIUM Large programme. In this work, the brightest CO, SiO, and HCN molecular line transitions are analysed. The continuum map shows two maxima, separated by 0.04'' (6 au). The CO data unambiguously reveal that Pi1 Gru's circumstellar environment harbours an inclined, radially outflowing, equatorial density enhancement. It contains a spiral structure at an angle of 38+/-3 deg with the line-of-sight. The HCN emission in the inner wind reveals a clockwise spiral, with a dynamical crossing time of the spiral arms consistent with a companion at a distance of 0.04'' from the asymptotic giant branch star, which is in agreement with the position of the secondary continuum peak. The inner wind dynamics imply a large acceleration region, consistent with a beta-law power of ~6. The CO emission suggests that the spiral is approximately Archimedean within 5'', beyond which this trend breaks down as the succession of the spiral arms becomes less periodic. The SiO emission at scales smaller than 0.5'' exhibits signatures of gas in rotation, which is found to fit the expected behaviour of gas in the wind-companion interaction zone. An investigation of SiO maser emission reveals what could be a stream of gas accelerating from the surface of the AGB star to the companion. Using these dynamics, we have tentatively derived an upper limit on the companion mass to be ~1.1 Msol.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Equatorial anti-rotating day side wind flow in WASP-43b elicited by deep wind jets?

We present WASP-43b climate simulations with deep wind jets (down to 700~bar) that are linked to retrograde (westward) flow at the equatorial day side for $p<0.1$~bar. Retrograde flow inhibits efficient eastward heat transport and naturally explains the small hotspot shift and large day-night-side gradient of WASP-43b ($P_{\text{orb}}=P_{\text{rot}}=0.8135$~days) observed with Spitzer. We find that deep wind jets are mainly associated with very fast rotations ($P_{\text{rot}}=P_{\text{orb}}\leq 1.5$~days) which correspond to the Rhines length smaller than $2$ planetary radii. We also diagnose wave activity that likely gives rise to deviations from superrotation. Further, we show that we can achieve full steady state in our climate simulations by imposing a deep forcing regime for $p>10$~bar: convergence time scale $τ_{\text{conv}}=10^6-10^8$~s to a common adiabat, as well as linear drag at depth ($p\geq 200$~bar), which mimics to first order magnetic drag. Lower boundary stability and the deep forcing assumptions were also tested with climate simulations for HD~209458b ($P_{\text{orb}}=P_{\text{rot}}=3.5$~days). HD~209458b simulations always show shallow wind jets (never deeper than 100~bar) and unperturbed superrotation. If we impose a fast rotation ($P_{\text{orb}}=P_{\text{rot}}=0.8135$~days), also the HD~209458b-like simulation shows equatorial retrograde flow at the day side. We conclude that the placement of the lower boundary at $p=200$~bar is justified for slow rotators like HD~209458b, but we suggest that it has to be placed deeper for fast-rotating, dense hot Jupiters ($P_{\text{orb}}\leq 1.5$~days) like WASP-43b. Our study highlights that the deep atmosphere may have a strong influence on the observable atmospheric flow in some hot Jupiters.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Indications for very high metallicity and absence of methane for the eccentric exo-Saturn WASP-117b

We investigate the atmospheric composition of the long period ($P_{\rm orb}=$ 10 days), eccentric exo-Saturn WASP-117b. WASP-117b could be in atmospheric temperature and chemistry similar to WASP-107b. In mass and radius WASP-117b is similar to WASP-39b, which allows a comparative study of these planets. We analyze a near-infrared transmission spectrum of WASP-117b taken with Hubble Space Telescope/WFC3 G141, which was reduced with two independent pipelines. High resolution measurements were taken with VLT/ESPRESSO in the optical. We report the robust ($3σ$) detection of a water spectral feature. Using a 1D atmosphere model with isothermal temperature, uniform cloud deck and equilibrium chemistry, the Bayesian evidence of a retrieval analysis of the transmission spectrum indicates a preference for a high atmospheric metallicity ${\rm [Fe/H]}=2.58^{+0.26}_{-0.37}$ and clear skies. The data are also consistent with a lower-metallicity composition ${\rm [Fe/H]}<1.75$ and a cloud deck between $10^{-2.2} - 10^{-5.1}$ bar, but with weaker Bayesian preference. We retrieve a low CH$_4$ abundance of $<10^{-4}$ volume fraction within $1 σ$ and $<2\cdot 10^{-1}$ volume fraction within $3 σ$. We cannot constrain the equilibrium temperature between theoretically imposed limits of 700 and 1000~K. Further observations are needed to confirm quenching of CH$_4$ with $K_{zz}\geq 10^8$~cm$^2$/s. We report indications of Na and K in the VLT/ESPRESSO high resolution spectrum with substantial Bayesian evidence in combination with HST data.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Magritte, a modern software library for 3D radiative transfer: I. Non-LTE atomic and molecular line modelling

Radiative transfer is a key component in almost all astrophysical and cosmological simulations. We present Magritte: a modern open-source software library for 3D radiative transfer. It uses a deterministic ray-tracer and formal solver, i.e. it computes the radiation field by tracing rays through the model and solving the radiative transfer equation in its second-order form along a fixed set of rays originating from each point. Magritte can handle structured and unstructured input meshes, as well as smoothed-particle hydrodynamics (SPH) particle data. In this first paper, we describe the numerical implementation, semi-analytic tests and cross-code benchmarks for the non-LTE line radiative transfer module of Magritte. This module uses the radiative transfer solver to self-consistently determine the populations of the quantised energy levels of atoms and molecules using an accelerated Lambda iteration (ALI) scheme. We compare Magritte with the established radiative transfer solvers Ratran (1D) and Lime (3D) on the van Zadelhoff benchmark and present a first application to a simple Keplerian disc model. Comparing with Lime, we conclude that Magritte produces more accurate and more precise results, especially at high optical depth, and that it is faster.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Measurements of Low Temperature Rate Coefficients for the Reaction of CH with CH$_{2}$O and Application to Dark Cloud and AGB Stellar Wind Models

Rate coefficients have been measured for the reaction of CH radicals with formaldehyde, CH$_{2}$O, over the temperature range 31 - 133 K using a pulsed Laval nozzle apparatus combined with pulsed laser photolysis and laser induced fluorescence spectroscopy. The rate coefficients are very large and display a distinct decrease with decreasing temperature below 70 K, although classical collision rate theory fails to reproduce this temperature dependence. The measured rate coefficients have been parameterized and used as input for astrochemical models for both dark cloud and AGB stellar outflow scenarios. The models predict a distinct change (up to a factor of two) in the abundance of ketene, H$_{2}$CCO, which is the major expected molecular product of the CH + CH$_{2}$O reaction.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Influence of Stellar Flares on the Chemical Composition of Exoplanets and Spectra

More than 3000 exoplanets have been detected so far, and more and more spectroscopic observations of exoplanets are performed. Future instruments are eagerly awaited as they will be able to provide spectroscopic data with a greater accuracy and sensitivity than what is currently available. An important aspect to consider is temporal stellar atmospheric disturbances that can influence the planetary composition, and hence spectra, and potentially can lead to incorrect assumptions about the steady-state atmospheric composition of the planet. We focus on perturbations that come from the host star in the form of flare events that significantly increase the photon flux impingement on the exoplanet atmosphere. In some cases, and particularly for M stars, this sudden increase may last for several hours. We aim at answering the question to what extent a stellar flare is able to modify the chemical composition of the planetary atmosphere and, therefore influence the resulting spectra. We use a 1D thermo-photochemical model to study the neutral atmospheric composition of two hypothetic planets located around the star AD Leo. This active star has already been observed during a flare. We use the spectroscopic data from this flare event to simulate the evolution of the chemical composition of the atmospheres of the two hypothetic planets. We compute synthetic spectra to evaluate the implications for observations. The increase of the incoming photon flux affects the chemical abundances of some important species down to altitudes associated with an atmospheric pressure of 1 bar, that can lead to variations in planetary spectra if performed during transit.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Connecting the dots II: Phase changes in the climate dynamics of tidally locked terrestrial exoplanets

We investigate 3D atmosphere dynamics for tidally locked terrestrial planets with an Earth-like atmosphere and irradiation for different rotation periods ($P_{rot}=1-100$ days) and planet sizes ($R_P=1-2 R_{Earth}$) with unprecedented fine detail. We could precisely identify three climate state transition regions that are associated with phase transitions in standing tropical and extra tropical Rossby waves. We confirm that the climate on fast rotating planets may assume multiple states ($P_{rot}\leq 12$ days for $R_P=2 R_{Earth}$). Our study is, however, the first to identify the type of planetary wave associated with different climate states: The first state is dominated by standing tropical Rossby waves with fast equatorial superrotation. The second state is dominated by standing extra tropical Rossby waves with high latitude westerly jets with slower wind speeds. For very fast rotations ($P_{rot}\leq 5$~days for $R_P=2 R_{Earth}$), we find another climate state transition, where the standing tropical and extra tropical Rossby wave can both fit on the planet. Thus, a third state with a mixture of the two planetary waves becomes possible that exhibits three jets. Different climate states may be observable, because the upper atmosphere's hot spot is eastward shifted with respect to the substellar point in the first state, westward shifted in the second state and the third state shows a longitudinal 'smearing' of the spot across the substellar point. We show, furthermore, that the largest fast rotating planet in our study exhibits atmosphere features known from hot Jupiters like fast equatorial superrotation and a temperature chevron in the upper atmosphere.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

New chemical scheme for studying carbon-rich exoplanet atmospheres

Atmospheres with a high C/O ratio are expected to contain an important quantity of hydrocarbons, including heavy molecules (with more than 2 carbon atoms). To study correctly these C-rich atmospheres, a chemical scheme adapted to this composition is necessary. We have implemented a chemical scheme that can describe the kinetics of species with up to 6 carbon atoms. This chemical scheme has been developed with specialists of combustion and validated through experiments on a wide range of T and P. This chemical network is available on the online database KIDA. We have created a grid of 12 models to explore different thermal profiles and C/O ratios. For each of them, we have compared the chemical composition determined with a C0-C2 chemical scheme (species with up to 2 carbon atoms) and with the C0-C6 scheme. We found no difference in the results obtained with the two schemes when photolyses are not included in the model, whatever the temperature of the atmosphere. In contrast, when there is photochemistry, differences can appear in the upper atmosphere. These differences are found for all the tested PT profiles in the case that the C/O ratio is above 1. When the C/O ratio of the atmosphere is solar, differences are only found at temperatures lower than 1000K. The differences linked to the use of different chemical schemes do not have important influence on the synthetic spectra. However, we have confirmed that C2H2 and HCN as possible tracers of warm C-rich atmospheres. The use of this new chemical scheme is mandatory to model atmospheres with a high C/O ratio and, in particular, if one is interested in studying in details the photochemistry. If one is just interested in the synthetic spectra, the use of a smaller scheme may be sufficient.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Simplified models of stellar wind anatomy for interpreting high-resolution data: Analytical approach to embedded spiral geometries

Recent high-resolution observations have shown stellar winds to harbour complexities which strongly deviate from spherical symmetry, generally assumed as standard wind model. One such morphology is the archimedean spiral, generally believed to be formed by binary interactions, which has been directly observed in multiple sources. We seek to investigate the manifestation in the observables of spiral structures embedded in the spherical outflows of cool stars. We aim to provide an intuitive bedrock with which upcoming ALMA data can be compared and interpreted. By means of an extended parameter study, we model rotational CO emission from the stellar outflow of asymptotic giant branch stars. To this end, we develop a simplified analytical parametrised description of a 3D spiral structure. This model is embedded into a spherical wind, and fed into the 3D radiative transfer code LIME, which produces 3D intensity maps throughout velocity space. Subsequently, we investigate the spectral signature of rotational transitions of CO of the models, as well as the spatial aspect of this emission by means of wide-slit PV diagrams. Additionally, the potential for misinterpretation of the 3D data in a 1D context is quantified. Finally, we simulate ALMA observations to explore the impact of interefrometric noise and artifacts on the emission signatures. The spectral signatures of the CO rotational transition v=0 J=3-2 are very efficient at concealing the dual nature of the outflow. Only a select few parameter combinations allow for the spectral lines to disclose the presence of the spiral structure. The inability to disentangle the spiral from the spherical signal can result in an incorrect interpretation in a 1D context. Consequently, erroneous mass loss rates would be calculated...

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The EChO science case

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A stubbornly large mass of cold dust in the ejecta of Supernova 1987A

We present new Herschel photometric and spectroscopic observations of Supernova 1987A, carried out in 2012. Our dedicated photometric measurements provide new 70 micron data and improved imaging quality at 100 and 160 micron compared to previous observations in 2010. Our Herschel spectra show only weak CO line emission, and provide an upper limit for the 63 micron [O I] line flux, eliminating the possibility that line contaminations distort the previously estimated dust mass. The far-infrared spectral energy distribution (SED) is well fitted by thermal emission from cold dust. The newly measured 70 micron flux constrains the dust temperature, limiting it to nearly a single temperature. The far-infrared emission can be fitted by 0.5+-0.1 Msun of amorphous carbon, about a factor of two larger than the current nucleosynthetic mass prediction for carbon. The observation of SiO molecules at early and late phases suggests that silicates may also have formed and we could fit the SED with a combination of 0.3 Msun of amorphous carbon and 0.5 Msun of silicates, totalling 0.8 Msun of dust. Our analysis thus supports the presence of a large dust reservoir in the ejecta of SN 1987A. The inferred dust mass suggests that supernovae can be an important source of dust in the interstellar medium, from local to high-redshift galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Connecting the dots: A versatile model for the atmospheres of tidally locked Super-Earths

Radiative equilibrium temperatures are calculated for the troposphere of a tidally locked Super-Earth based on a simple greenhouse model, using Solar System data as a guideline. These temperatures provide in combination with a Newtonian relaxation scheme thermal forcing for a 3D atmosphere model using the dynamical core of the Massachusetts Institute of Technology global circulation model (MITgcm). Our model is of the same conceptional simplicity than the model of Held & Suarez (1994) and is thus computationally fast. Furthermore, because of the coherent, general derivation of radiative equilibrium temperatures, our model is easily adaptable for different planets and atmospheric scenarios. As a case study relevant for Super-Earths, we investigate a Gl581g-like planet with Earth-like atmosphere and irradiation and present results for two representative rotation periods of $P_{rot}$ = 10 days and $P_{rot}$ = 36.5 days. Our results provide proof of concept and highlight interesting dynamical features for the rotating regime 3 < $P_{rot}$ < 100 days, which was shown by Edson et al. (2011) to be an intermediate regime between equatorial superrotation and divergence. We confirm that the $P_{rot}$ = 10 days case is more dominated by equatorial superrotation dynamics than the $P_{rot}$ = 36.5 days case, which shows diminishing influence of standing Rossby- Kelvin waves and increasing influence of divergence at the top of the atmosphere. We argue that this dynamical regime change relates to the increase in Rossby deformation radius, in agreement with previous studies. However, we also pay attention to other features that are not or only in partial agreement with other studies, like, e.g., the number of circulation cells and their strength, the role and extent of thermal inversion layers, and the details of heat transport

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Influence of different parameters on the chemical composition of warm neptunes

We developed a 1D photo-thermochemical model to study the atmosphere of warm exoplanets. The chemical scheme used in this model is completely new in planetology and has been constructed in collaboration with specialists of combustion. It has been validated as a whole through experiments on a large range of temperature (300 - 2500 K) and pressure (1 mbar - 100 bar), allowing to study a wide variety of exoplanets. We have used this chemical model to study the atmosphere of two warm Neptunes, GJ3470b and GJ436b, and the influence of different parameters (vertical mixing, metallicity, temperature, . . . ) on their chemical composition. We present here the results obtained in these studies.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Using numerical models of bow shocks to investigate the circumstellar medium of massive stars

Many massive stars travel through the interstellar medium at supersonic speeds. As a result they form bow shocks at the interface between the stellar wind. We use numerical hydrodynamics to reproduce such bow shocks numerically, creating models that can be compared to observations. In this paper we discuss the influence of two physical phenomena, interstellar magnetic fields and the presence of interstellar dust grains on the observable shape of the bow shocks of massive stars. We find that the interstellar magnetic field, though too weak to restrict the general shape of the bow shock, reduces the size of the instabilities that would otherwise be observed in the bow shock of a red supergiant. The interstellar dust grains, due to their inertia can penetrate deep into the bow shock structure of a main sequence O-supergiant, crossing over from the ISM into the stellar wind. Therefore, the dust distribution may not always reflect the morphology of the gas. This is an important consideration for infrared observations, which are dominated by dust emission. Our models clearly show, that the bow shocks of massive stars are useful diagnostic tools that can used to investigate the properties of both the stellar wind as well as the interstellar medium.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Can the magnetic field in the Orion arm inhibit the growth of instabilities in the bow shock of Betelgeuse?

Many evolved stars travel through space at supersonic velocities, which leads to the formation of bow shocks ahead of the star where the stellar wind collides with the interstellar medium (ISM). Herschel observations of the bow shock of $α$-Orionis show that the shock is almost free of instabilities, despite being, at least in theory, subject to both Kelvin-Helmholtz and Rayleigh-Taylor instabilities. A possible explanation for the lack of instabilities lies in the presence of an interstellar magnetic field. We wish to investigate whether the magnetic field of the interstellar medium (ISM) in the Orion arm can inhibit the growth of instabilities in the bow shock of $α$-Orionis. We used the code MPI-AMRVAC to make magneto-hydrodynamic simulations of a circumstellar bow shock, using the wind parameters derived for $α$-Orionis and interstellar magnetic field strengths of $B\,=\,1.4,\, 3.0$, and $5.0\, μ$G, which fall within the boundaries of the observed magnetic field strength in the Orion arm of the Milky Way. Our results show that even a relatively weak magnetic field in the interstellar medium can suppress the growth of Rayleigh-Taylor and Kelvin-Helmholtz instabilities, which occur along the contact discontinuity between the shocked wind and the shocked ISM. The presence of even a weak magnetic field in the ISM effectively inhibits the growth of instabilities in the bow shock. This may explain the absence of such instabilities in the Herschel observations of $α$-Orionis.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Herschel/HIFI search for H2-17O and H2-18O in IRC+10216: constraints on models for the origin of water vapor

We report the results of a sensitive search for the minor isotopologues of water, H2-17O and H2-18O, toward the carbon-rich AGB star IRC+10216 (a.k.a. CW Leonis) using the HIFI instrument on the Herschel Space Observatory. This search was motivated by the fact that any detection of isotopic enhancement in the H2-17O and H2-18O abundances would have strongly implicated CO photodissociation as the source of the atomic oxygen needed to produce water in a carbon-rich circumstellar envelope. Our observations place an upper limit of 1/470 on the H2-17O/H2-16O abundance ratio. Given the isotopic 17O/16O ratio of 1/840 inferred previously for the photosphere of IRC+10216, this result places an upper limit of a factor 1.8 on the extent of any isotope-selective enhancement of H2-17O in the circumstellar material, and provides an important constraint on any model that invokes CO photodissociation as the source of O for H2O production. In the context of the clumpy photodissociation model proposed previously for the origin of water in IRC+10216, our limit implies that 12C-16O (not 13C-16O or SiO) must be the dominant source of 16O for H2O production, and that the effects of self-shielding can only have reduced the 12C-16O photodissociation rate by at most a factor ~ 2.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Late Stages of Stellar Evolution - Herschel's contributions

Cool objects glow in the infrared. The gas and solid-state species that escape the stellar gravitational attraction of evolved late-type stars in the form of a stellar wind are cool, with temperatures typically $\la$1500\,K, and can be ideally studied in the infrared. These stellar winds create huge extended circumstellar envelopes with extents approaching $10^{19}$\,cm. In these envelopes, a complex kinematical, thermodynamical and chemical interplay determines the global and local structural parameters. Unraveling the wind acceleration mechanisms and deriving the complicated structure of the envelopes is important to understand the late stages of evolution of ~97% of stars in galaxies as our own Milky Way. That way, we can also assess the significant chemical enrichment of the interstellar medium by the mass loss of these evolved stars. The Herschel Space Observatory is uniquely placed to study evolved stars thanks to the excellent capabilities of the three infrared and sub-millimeter instruments on board: PACS, SPIRE and HIFI. In this review, I give an overview of a few important results obtained during the first two years of Herschel observations in the field of evolved low and intermediate mass stars, and I will show how the Herschel observations can solve some historical questions on these late stages of stellar evolution, but also add some new ones.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Computing the dust distribution in the bowshock of a fast moving, evolved star

We study the hydrodynamical behavior occurring in the turbulent interaction zone of a fast moving red supergiant star, where the circumstellar and interstellar material collide. In this wind-interstellar medium collision, the familiar bow shock, contact discontinuity, and wind termination shock morphology forms, with localized instability development. Our model includes a detailed treatment of dust grains in the stellar wind, and takes into account the drag forces between dust and gas. The dust is treated as pressureless gas components binned per grainsize, for which we use ten representative grainsize bins. Our simulations allow to deduce how dust grains of varying sizes become distributed throughout the circumstellar medium. We show that smaller dust grains (radius <0.045 micro-meters) tend to be strongly bound to the gas and therefore follow the gas density distribution closely, with intricate finestructure due to essentially hydrodynamical instabilities at the wind-related contact discontinuity. Larger grains which are more resistant to drag forces are shown to have their own unique dust distribution, with progressive deviations from the gas morphology. Specifically, small dust grains stay entirely within the zone bound by shocked wind material. The large grains are capable of leaving the shocked wind layer, and can penetrate into the shocked or even unshocked interstellar medium. Depending on how the number of dust grains varies with grainsize, this should leave a clear imprint in infrared observations of bowshocks of red supergiants and other evolved stars.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dust distribution in circumstellar shells

We present numerical simulations of the hydrodynamical interactions that produce circumstellar shells. These simulations include several scenarios, such as wind-wind interaction and wind-ISM collisions. In our calculations we have taken into account the presence of dust in the stellar wind. Our results show that, while small dust grains tend to be strongly coupled to the gas, large dust grains are only weakly coupled. As a result, the distribution of the large dust grains is not representative of the gas distribution. Combining these results with observations may give us a new way of validating hydrodynamical models of the circumstellar medium.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Circumstellar molecular composition of the oxygen-rich AGB star IK Tau: I. Observations and LTE chemical abundance analysis

The aim of this paper is to study the molecular composition in the circumstellar envelope around the oxygen-rich star IK Tau. We observed IK Tau in several (sub)millimeter bands using the APEX telescope during three observing periods. To determine the spatial distribution of the $\mathrm{^{12}CO(3-2)}$ emission, mapping observations were performed. To constrain the physical conditions in the circumstellar envelope, multiple rotational CO emission lines were modeled using a non local thermodynamic equilibrium radiative transfer code. The rotational temperatures and the abundances of the other molecules were obtained assuming local thermodynamic equilibrium. An oxygen-rich Asymptotic Giant Branch star has been surveyed in the submillimeter wavelength range. Thirty four transitions of twelve molecular species, including maser lines, were detected. The kinetic temperature of the envelope was determined and the molecular abundance fractions of the molecules were estimated. The deduced molecular abundances were compared with observations and modeling from the literature and agree within a factor of 10, except for SO$_2$, which is found to be almost a factor 100 stronger than predicted by chemical models. From this study, we found that IK Tau is a good laboratory to study the conditions in circumstellar envelopes around oxygen-rich stars with (sub)millimeter-wavelength molecular lines. We could also expect from this study that the molecules in the circumstellar envelope can be explained more faithful by non-LTE analysis with lower and higher transition lines than by simple LTE analysis with only lower transition lines. In particular, the observed CO line profiles could be well reproduced by a simple expanding envelope model with a power law structure.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Herschel/HIFI observations of IRC+10216: water vapor in the inner envelope of a carbon-rich AGB star

We report the results of observations of ten rotational transitions of water vapor toward the carbon-rich AGB (asymptotic giant branch) star IRC+10216 (CW Leonis), carried out with Herschel's HIFI instrument. Each transition was securely detected by means of observations using the dual beam switch mode of HIFI. The measured line ratios imply that water vapor is present in the inner outflow at small distances (few x 1.E+14 cm) from the star, confirming recent results reported by Decin et al. from observations with Herschel's PACS and SPIRE instruments. This finding definitively rules out the hypothesis that the observed water results from the vaporization of small icy objects in circular orbits. The origin of water within the dense C-rich envelope of IRC+10216 remains poorly understood. We derive upper limits on the H2-17O/H2-16O and H2-18O/H2-16O isotopic abundance ratios of ~ 5.E-3 (3 sigma), providing additional constraints on models for the origin of the water vapor in IRC+10216.

0 citations0 reviewsNo code linkedOpen access