Chemistry and Distribution of Daughter Species in the Circumstellar Envelopes of O-Rich AGB Stars

Context. Thanks to the advent of Herschel and ALMA, new high-quality observations of molecules present in the circumstellar envelopes of asymptotic giant branch (AGB) stars are being reported that reveal large differences from the existing chemical models. New molecular data and more comprehensive models of the chemistry in circumstellar envelopes are now available.
Aims: The aims are to determine and study the important formation and destruction pathways in the envelopes of O-rich AGB stars and to provide more reliable predictions of abundances, column densities, and radial distributions for potentially detectable species with physical conditions applicable to the envelope surrounding IK Tau.
Methods: We use a large gas-phase chemical model of an AGB envelope including the effects of CO and N2 self-shielding in a spherical geometry and a newly compiled list of inner-circumstellar envelope parent species derived from detailed modeling and observations. We trace the dominant chemistry in the expanding envelope and investigate the chemistry as a probe for the physics of the AGB phase by studying variations of abundances with mass-loss rates and expansion velocities.
Results: We find a pattern of daughter molecules forming from the photodissociation products of parent species with contributions from ion-neutral abstraction and dissociative recombination. The chemistry in the outer zones differs from that in traditional PDRs in that photoionization of daughter species plays a significant role. With the proper treatment of self-shielding, the N → N2 and C+→ CO transitions are shifted outward by factors of 7 and 2, respectively, compared with earlier models. An upper limit on the abundance of CH4 as a parent species of (≲2.5 × 10-6 with respect to H2) is found for IK Tau, and several potentially observable molecules with relatively simple chemical links to other parent species are determined. The assumed stellar mass-loss rate, in particular, has an impact on the calculated abundances of cations and the peak-abundance radius of both cations and neutrals: as the mass-loss rate increases, the peak abundance of cations generally decreases and the peak-abundance radius of all species moves outwards. The effects of varying the envelope expansion velocity and cosmic-ray ionization rate are not as significant.

The VLT-FLAMES Tarantula Survey XXII. Multiplicity properties of the B-type stars

We investigate the multiplicity properties of 408 B-type stars observed in the 30 Doradus region of the Large Magellanic Cloud with multi-epoch spectroscopy from the VLT-FLAMES Tarantula Survey (VFTS). We use a cross-correlation method to estimate relative radial velocities from the helium and metal absorption lines for each of our targets. Objects with significant radial-velocity variations (and with an amplitude larger than 16 km s^-1) are classified as spectroscopic binaries. We find an observed spectroscopic binary fraction (defined by periods of <10^3.5 d and mass ratios >0.1) for the B-type stars, f<inf>B</inf>(obs) = 0.25 ± 0.02, which appears constant across the field of view, except for the two older clusters (Hodge 301 and SL 639). These two clusters have significantly lower binary fractions of 0.08 ± 0.08 and 0.10 ± 0.09, respectively. Using synthetic populations and a model of our observed epochs and their potential biases, we constrain the intrinsic multiplicity properties of the dwarf and giant (i.e. relatively unevolved) B-type stars in 30 Dor. We obtain a present-day binary fraction f<inf>B</inf>(true) = 0.58 ± 0.11, with a flat period distribution. Within the uncertainties, the multiplicity properties of the B-type stars agree with those for the O stars in 30 Dor from the VFTS.

On Variable Selection of a Nonlinear Non-parametric System with a Limited Data Set: A Stepwise Algorithm

his paper considers a problem of identification for a high dimensional nonlinear non-parametric system when only a limited data set is available. The algorithms are proposed for this purpose which exploit the relationship between the input variables and the output and further the inter-dependence of input variables so that the importance of the input variables can be established. A key to these algorithms is the non-parametric two stage input selection algorithm.

Self-consistent modelling of line-driven hot-star winds with Monte Carlo radiation hydrodynamics

Radiative pressure exerted by line interactions is a prominent driver of outflows in astrophysical systems, being at work in the outflows emerging from hot stars or from the accretion discs of cataclysmic variables, massive young stars and active galactic nuclei. In this work, a new radiation hydrodynamical approach to model line-driven hot-star winds is presented. By coupling a Monte Carlo radiative transfer scheme with a finite volume fluid dynamical method, line-driven mass outflows may be modelled self-consistently, benefiting from the advantages of Monte Carlo techniques in treating multiline effects, such as multiple scatterings, and in dealing with arbitrary multidimensional configurations. In this work, we introduce our approach in detail by highlighting the key numerical techniques and verifying their operation in a number of simplified applications, specifically in a series of self-consistent, one-dimensional, Sobolev-type, hot-star wind calculations. The utility and accuracy of our approach are demonstrated by comparing the obtained results with the predictions of various formulations of the so-called CAK theory and by confronting the calculations with modern sophisticated techniques of predicting the wind structure. Using these calculations, we also point out some useful diagnostic capabilities our approach provides. Finally, we discuss some of the current limitations of our method, some possible extensions and potential future applications.

Coronal properties of planet-bearing stars

Context. Do extrasolar planets affect the activity of their host stars? Indications for chromospheric activity enhancement have been found for a handful of targets, but in the X-ray regime, conclusive observational evidence is still missing. Aims: We want to establish a sound observational basis to confirm or reject major effects of Star-Planet Interactions (SPI) in stellar X-ray emissions. Methods: We therefore conduct a statistical analysis of stellar X-ray activity of all known planet-bearing stars within 30 pc distance for dependencies on planetary parameters such as mass and semimajor axis. Results: In our sample, there are no significant correlations of X-ray luminosity or the activity indicator L_X/L_bol with planetary parameters which cannot be explained by selection effects. Conclusions: Coronal SPI seems to be a phenomenon which might only manifest itself as a strong effect for a few individual targets, but not to have a major effect on planet-bearing stars in general.

CSI 2264: Simultaneous Optical and Infrared Light Curves of Young Disk-bearing Stars in NGC 2264 with CoRoT and Spitzer—Evidence for Multiple Origins of Variability

We present the Coordinated Synoptic Investigation of NGC 2264, a continuous 30 day multi-wavelength photometric monitoring campaign on more than 1000 young cluster members using 16 telescopes. The unprecedented combination of multi-wavelength, high-precision, high-cadence, and long-duration data opens a new window into the time domain behavior of young stellar objects. Here we provide an overview of the observations, focusing on results from Spitzer and CoRoT. The highlight of this work is detailed analysis of 162 classical T Tauri stars for which we can probe optical and mid-infrared flux variations to 1% amplitudes and sub-hour timescales. We present a morphological variability census and then use metrics of periodicity, stochasticity, and symmetry to statistically separate the light curves into seven distinct classes, which we suggest represent different physical processes and geometric effects. We provide distributions of the characteristic timescales and amplitudes and assess the fractional representation within each class. The largest category (>20%) are optical "dippers" with discrete fading events lasting ~1-5 days. The degree of correlation between the optical and infrared light curves is positive but weak; notably, the independently assigned optical and infrared morphology classes tend to be different for the same object. Assessment of flux variation behavior with respect to (circum)stellar properties reveals correlations of variability parameters with Hα emission and with effective temperature. Overall, our results point to multiple origins of young star variability, including circumstellar obscuration events, hot spots on the star and/or disk, accretion bursts, and rapid structural changes in the inner disk. Based on data from the Spitzer and CoRoT missions. The CoRoT space mission was developed and is operated by the French space agency CNES, with participation of ESA's RSSD and Science Programmes, Austria, Belgium, Brazil, Germany, and Spain.

Indications for an influence of hot Jupiters on the rotation and activity of their host stars

Context. The magnetic activity of planet-hosting stars is an importantfactor for estimating the atmospheric stability of close-in exoplanetsand the age of their host stars. It has long been speculated thatclose-in exoplanets can influence the stellar activity level. However,testing for tidal or magnetic interaction effects in samples ofplanet-hosting stars is difficult because stellar activity hindersexoplanet detection, so that stellar samples with detected exoplanetsshow a bias toward low activity for small exoplanets.

Aims: Weaim to test whether exoplanets in close orbits influence the stellarrotation and magnetic activity of their host stars.

Methods: Wedeveloped a novel approach to test for systematic activity-enhancementsin planet-hosting stars. We use wide (several 100 AU) binary systems inwhich one of the stellar components is known to have an exoplanet, whilethe second stellar component does not have a detected planet andtherefore acts as a negative control. We use the stellar coronal X-rayemission as an observational proxy for magnetic activity and analyzeobservations performed with Chandra and XMM-Newton.

Results: Wefind that in two systems for which strong tidal interaction can beexpected the planet-hosting primary displays a much higher magneticactivity level than the planet-free secondary. In three systems forwhich weaker tidal interaction can be expected the activity levels ofthe two stellar components agree with each other.

Conclusions:Our observations indicate that the presence of Hot Jupiters may inhibitthe spin-down of host stars with thick outer convective layers. Possiblecauses for this effect include a transfer of angular momentum from theplanetary orbit to the stellar rotation through tidal interaction, ordifferences during the early evolution of the system, where the hoststar may decouple from the protoplanetary disk early because of a gapopened by the forming Hot Jupiter.

Young Stellar Object VARiability (YSOVAR): Long Timescale Variations in the Mid-infrared

The YSOVAR (Young Stellar Object VARiability) Spitzer Space Telescope observing program obtained the first extensive mid-infrared (3.6 and 4.5 μm) time series photometry of the Orion Nebula Cluster plus smaller footprints in 11 other star-forming cores (AFGL 490, NGC 1333, Mon R2, GGD 12-15, NGC 2264, L1688, Serpens Main, Serpens South, IRAS 20050+2720, IC 1396A, and Ceph C). There are ~29,000 unique objects with light curves in either or both IRAC channels in the YSOVAR data set. We present the data collection and reduction for the Spitzer and ancillary data, and define the "standard sample" on which we calculate statistics, consisting of fast cadence data, with epochs roughly twice per day for ~40 days. We also define a "standard sample of members" consisting of all the IR-selected members and X-ray-selected members. We characterize the standard sample in terms of other properties, such as spectral energy distribution shape. We use three mechanisms to identify variables in the fast cadence data—the Stetson index, a χ2 fit to a flat light curve, and significant periodicity. We also identified variables on the longest timescales possible of six to seven years by comparing measurements taken early in the Spitzer mission with the mean from our YSOVAR campaign. The fraction of members in each cluster that are variable on these longest timescales is a function of the ratio of Class I/total members in each cluster, such that clusters with a higher fraction of Class I objects also have a higher fraction of long-term variables. For objects with a YSOVAR-determined period and a [3.6]-[8] color, we find that a star with a longer period is more likely than those with shorter periods to have an IR excess. We do not find any evidence for variability that causes [3.6]-[4.5] excesses to appear or vanish within our data set; out of members and field objects combined, at most 0.02% may have transient IR excesses.

Magnetospheric Structure and Atmospheric Joule Heating of Habitable Planets Orbiting M-dwarf Stars

We study the magnetospheric structure and the ionospheric Joule Heating of planets orbiting M-dwarf stars in the habitable zone using a set of magnetohydrodynamic models. The stellar wind solution is used to drive a model for the planetary magnetosphere, which is coupled with a model for the planetary ionosphere. Our simulations reveal that the space environment around close-in habitable planets is extreme, and the stellar wind plasma conditions change from sub- to super-Alfvénic along the planetary orbit. As a result, the magnetospheric structure changes dramatically with a bow shock forming in the super-Alfvénic sectors, while no bow shock forms in the sub-Alfvénic sectors. The planets reside most of the time in the sub-Alfvénic sectors with poor atmospheric protection. A significant amount of Joule Heating is provided at the top of the atmosphere as a result of the intense stellar wind. For the steady-state solution, the heating is about 0.1%-3% of the total incoming stellar irradiation, and it is enhanced by 50% for the time-dependent case. The significant Joule Heating obtained here should be considered in models for the atmospheres of habitable planets in terms of the thickness of the atmosphere, the top-side temperature and density, the boundary conditions for the atmospheric pressure, and particle radiation and transport. Here we assume constant ionospheric Pedersen conductance similar to that of the Earth. The conductance could be greater due to the intense EUV radiation leading to smaller heating rates. We plan to quantify the ionospheric conductance in future study.

YSOVAR: Mid-infrared Variability in the Star-forming Region Lynds 1688

The emission from young stellar objects (YSOs) in the mid-infrared (mid-IR) is dominated by the inner rim of their circumstellar disks. We present IR data from the Young Stellar Object VARiability (YSOVAR) survey of ~800 objects in the direction of the Lynds 1688 (L1688) star-forming region over four visibility windows spanning 1.6 yr using the Spitzer Space Telescope in its warm mission phase. Among all light curves, 57 sources are cluster members identified based on their spectral energy distribution and X-ray emission. Almost all cluster members show significant variability. The amplitude of the variability is larger in more embedded YSOs. Ten out of 57 cluster members have periodic variations in the light curves with periods typically between three and seven days, but even for those sources, significant variability in addition to the periodic signal can be seen. No period is stable over 1.6 yr. Nonperiodic light curves often still show a preferred timescale of variability that is longer for more embedded sources. About half of all sources exhibit redder colors in a fainter state. This is compatible with time-variable absorption toward the YSO. The other half becomes bluer when fainter. These colors can only be explained with significant changes in the structure of the inner disk. No relation between mid-IR variability and stellar effective temperature or X-ray spectrum is found.

Explaining the Coexistence of Large-scale and Small-scale Magnetic Fields in Fully Convective Stars

Despite the lack of a shear-rich tachocline region, low-mass fully convective (FC) stars are capable of generating strong magnetic fields, indicating that a dynamo mechanism fundamentally different from the solar dynamo is at work in these objects. We present a self-consistent three-dimensional model of magnetic field generation in low-mass FC stars. The model utilizes the anelastic magnetohydrodynamic equations to simulate compressible convection in a rotating sphere. A distributed dynamo working in the model spontaneously produces a dipole-dominated surface magnetic field of the observed strength. The interaction of this field with the turbulent convection in outer layers shreds it, producing small-scale fields that carry most of the magnetic flux. The Zeeman–Doppler-Imaging technique applied to synthetic spectropolarimetric data based on our model recovers most of the large-scale field. Our model simultaneously reproduces the morphology and magnitude of the large-scale field as well as the magnitude of the small-scale field observed on low-mass FC stars.

YSOVAR: Mid-infrared Variability Among YSOs in the Star Formation Region GGD12-15

We present an IR-monitoring survey with the Spitzer Space Telescope of the star-forming region GGD 12-15. More than 1000 objects were monitored, including about 350 objects within the central 5′, which is found to be especially dense in cluster members. The monitoring took place over 38 days and is part of the Young Stellar Object VARiability project. The region was also the subject of a contemporaneous 67 ks Chandra observation. The field includes 119 previously identified pre-main sequence star candidates. X-rays are detected from 164 objects, 90 of which are identified with cluster members. Overall, we find that about half the objects in the central 5′ are young stellar objects (YSOs) based on a combination of their spectral energy distribution, IR variability, and X-ray emission. Most of the stars with IR excess relative to a photosphere show large amplitude (>0.1 mag) mid-infrared (mid-IR) variability. There are 39 periodic sources, and all but one is found to be a cluster member. Almost half of the periodic sources do not show IR excesses. Overall, more than 85% of the Class I, flat spectrum, and Class II sources are found to vary. The amplitude of the variability is larger in more embedded YSOs. Most of the Class I/ II objects exhibit redder colors in a fainter state, which is compatible with time-variable extinction. A few become bluer when fainter, which can be explained with significant changes in the structure of the inner disk. A search for changes in the IR due to X-ray events is carried out, but the low number of flares prevented an analysis of the direct impact of X-ray flares on the IR light curves. However, we find that X-ray detected Class II sources have longer timescales for change in the MIR than a similar set of non-X-ray detected Class IIs.

YSOVAR: Mid-infrared Variability in NGC 1333

As part of the Young Stellar Object VARiability (YSOVAR) program, wemonitored NGC 1333 for ∼35 days at 3.6 and 4.5 μm using theSpitzer Space Telescope. We report here on the mid-infrared variabilityof the point sources in the ∼10‧ × ∼20‧ areacentered on 03:29:06, +31:19:30 (J2000). Out of 701 light curves ineither channel, we find 78 variables over the YSOVAR campaign. Abouthalf of the members are variable. The variable fraction for the mostembedded spectral energy distributions (SEDs) (Class I, flat) is higherthan that for less embedded SEDs (Class II), which is in turn higherthan the star-like SEDs (Class III). A few objects have amplitudes(10–90th percentile brightness) in [3.6] or [4.5] > 0.2 mag; amore typical amplitude is 0.1–0.15 mag. The largest color changeis >0.2 mag. There are 24 periodic objects, with 40% of them beingflat SED class. This may mean that the periodic signal is primarily fromthe disk, not the photosphere, in those cases. We find 9 variableslikely to be “dippers,” where texture in the disk occultsthe central star, and 11 likely to be “bursters,” whereaccretion instabilities create brightness bursts. There are 39 objectsthat have significant trends in [3.6]–[4.5] color over thecampaign, about evenly divided between redder-when-fainter (consistentwith extinction variations) and bluer-when-fainter. About a third of the17 Class 0 and/or jet-driving sources from the literature are variableover the YSOVAR campaign, and a larger fraction (∼half) are variablebetween the YSOVAR campaign and the cryogenic-era Spitzer observations(6–7 years), perhaps because it takes time for the envelope torespond to changes in the central source. The NGC 1333 brown dwarfs donot stand out from the stellar light curves in any way except there is amuch larger fraction of periodic objects (∼60% of variable browndwarfs are periodic, compared to ∼30% of the variables overall).

An age-activity calibration for old low-mass stars

Low-mass stars are highly interesting targets: we are able to detect planets in their habitable zones, and upcoming searches for biomarkers in exoplanet atmospheres will focus on low-mass star systems due to their ubiquity and proximity. We aim to develop an age-activity calibration for old low-mass stars, using wide binary systems consisting of an M or K dwarf and a white dwarf. The age of the system is determined by the WD cooling time plus its progenitor lifetime, yielding reliable ages in the regime >1 Gyr. For an exploratory sample of 7 systems where we have already derived ages, we propose to perform Chandra ACIS-S observations to determine the X-ray luminosities of the M dwarfs and correlate their stellar activity with age. We ask for a total observing time of 110 ks.