ALMA-resolved salt emission traces the chemical footprint and inner wind morphology of VY Canis Majoris

(abreviated) We aim to study the inner-wind structure (R<250 Rstar) of the well-known red supergiant VY CMa. We analyse high spatial resolution (~0".24x0".13) ALMA Science Verification (SV) data in band 7 in which four thermal emission lines of gaseous sodium chloride (NaCl) are present at high signal-to-noise ratio. For the first time, the NaCl emission in the inner wind region of VY CMa is spatially resolved. The ALMA observations reveal the contribution of up to four different spatial regions. The NaCl emission pattern is different compared to the dust continuum and TiO2 emission already analysed from the ALMA SV data. The emission can be reconciled with an axisymmetric geometry, where the lower density polar/rotation axis has a position angle of ~50 degrees measured from north to east. However, this picture can not capture the full morphological diversity, and discrete mass ejection events need to be invoked to explain localized higher-density regions. The velocity traced by the gaseous NaCl line profiles is significantly lower than the average wind terminal velocity, and much slower than some of the fastest mass ejections, signalling a wide range of characteristic speeds for the mass loss. Gaseous NaCl is detected far beyond the main dust condensation region. Realising the refractory nature of this metal halide, this hints at a chemical process preventing all NaCl from condensing onto dust grains. We show that in the case of the ratio of the surface binding temperature to the grain temperature being ~50, only some 10% of NaCl remains in gaseous form, while for lower values of this ratio thermal desorption efficiently evaporates NaCl. Photodesorption by stellar photons seems not to be a viable explanation for the detection of gaseous NaCl at 220 Rstar from the central star, and instead, we propose shock-induced sputtering driven by localized mass ejection events as alternative.

SN 2011A: A low-luminosity interacting transient with a double plateau and strong sodium absorpton

We present optical photometry and spectroscopy of the optical transient SN 2011A. Our data span 140 days after discovery including BVRI u′g′r′i′z′ photometry and 11 epochs of optical spectroscopy. Originally classified as a type IIn supernova (SN IIn) due to the presence of narrow Hα emission, this object shows exceptional characteristics. First, the light curve shows a double plateau, a property only observed before in the impostor SN 1997bs. Second, SN 2011A has a very low luminosity (M_V=-15.72), placing it between normal luminous SNe IIn and SN impostors. Third, SN 2011A shows low velocity and high equivalent width absorption close to the sodium doublet, which increases with time and is most likely of circumstellar origin. This evolution is also accompanied by a change in line profile; when the absorption becomes stronger, a P Cygni profile appears. We discuss SN 2011A in the context of interacting SNe IIn and SN impostors, which appears to confirm the uniqueness of this transient. While we favor an impostor origin for SN 2011A, we highlight the difficulty in differentiating between terminal and non-terminal interacting transients.

The Kepler-454 system: a small, not-rocky inner planet, a Jovian world, and a distant companion

Kepler-454 (KOI-273) is a relatively bright (V = 11.69 mag), Sun-like star that hosts a transiting planet candidate in a 10.6 day orbit. From spectroscopy, we estimate the stellar temperature to be 5687 ± 50 K, its metallicity to be [m/H] = 0.32 ± 0.08, and the projected rotational velocity to be v sin i <2.4 km s^-1. We combine these values with a study of the asteroseismic frequencies from short cadence Kepler data to estimate the stellar mass to be , the radius to be 1.066 ± 0.012 R_o, and the age to be Gyr. We estimate the radius of the 10.6 day planet as 2.37 ± 0.13 R_⊕. Using 63 radial velocity observations obtained with the HARPS-N spectrograph on the Telescopio Nazionale Galileo and 36 observations made with the HIRES spectrograph at the Keck Observatory, we measure the mass of this planet to be 6.8 ± 1.4 M_⊕. We also detect two additional non-transiting companions, a planet with a minimum mass of 4.46 ± 0.12 M_J in a nearly circular 524 day orbit and a massive companion with a period >10 years and mass >12.1 M_J. The 12 exoplanets with radii ⊕ and precise mass measurements appear to fall into two populations, with those ⊕ following an Earth-like composition curve and larger planets requiring a significant fraction of volatiles. With a density of 2.76 ± 0.73 g cm^-3, Kepler-454b lies near the mass transition between these two populations and requires the presence of volatiles and/or H/He gas.

The VLT-FLAMES Tarantula Survey: XIX. B-type supergiants: Atmospheric parameters and nitrogen abundances to investigate the role of binarity and the width of the main sequence

Context: Model atmosphere analyses have been previously undertaken for both Galactic and extragalactic B-type supergiants. By contrast, little attention has been given to a comparison of the properties of single supergiants and those that are members of multiple systems. 

Aims: Atmospheric parameters and nitrogen abundances have been estimated for all the B-type supergiants identified in the VLT-FLAMES Tarantula survey. These include both single targets and binary candidates. The results have been analysed to investigate the role of binarity in the evolutionary history of supergiants. 

Methods: tlusty non-local thermodynamic equilibrium (LTE) model atmosphere calculations have been used to determine atmospheric parameters and nitrogen abundances for 34 single and 18 binary supergiants. Effective temperatures were deduced using the silicon balance technique, complemented by the helium ionisation in the hotter spectra. Surface gravities were estimated using Balmer line profiles and microturbulent velocities deduced using the silicon spectrum. Nitrogen abundances or upper limits were estimated from the Nii spectrum. The effects of a flux contribution from an unseen secondary were considered for the binary sample. Results. We present the first systematic study of the incidence of binarity for a sample of B-type supergiants across the theoretical terminal age main sequence (TAMS). To account for the distribution of effective temperatures of the B-type supergiants it may be necessary to extend the TAMS to lower temperatures. This is also consistent with the derived distribution of mass discrepancies, projected rotational velocities and nitrogen abundances, provided that stars cooler than this temperature are post-red supergiant objects. For all the supergiants in the Tarantula and in a previous FLAMES survey, the majority have small projected rotational velocities. The distribution peaks at about 50 km s^-1 with 65% in the range 30 km s^-1 ≤ ν_e sin i ≤ 60 km s^-1. About ten per cent have larger ve sin i (≥100 km s^-1), but surprisingly these show little or no nitrogen enhancement. All the cooler supergiants have low projected rotational velocities of ≤70 km s^-1 and high nitrogen abundance estimates, implying that either bi-stability braking or evolution on a blue loop may be important. Additionally, there is a lack of cooler binaries, possibly reflecting the small sample sizes. Single-star evolutionary models, which include rotation, can account for all of the nitrogen enhancement in both the single and binary samples. The detailed distribution of nitrogen abundances in the single and binary samples may be different, possibly reflecting differences in their evolutionary history. 

Conclusions: The first comparative study of single and binary B-type supergiants has revealed that the main sequence may be significantly wider than previously assumed, extending to T_eff = 20 000 K. Some marginal differences in single and binary atmospheric parameters and abundances have been identified, possibly implying non-standard evolution for some of the sample. This sample as a whole has implications for several aspects of our understanding of the evolutionary status of blue supergiants.

The VLT-FLAMES Tarantula Survey: XVIII. Classifications and radial velocities of the B-type stars

We present spectral classifications for 438 B-type stars observed as part of the VLT-FLAMES Tarantula Survey (VFTS) in the 30 Doradus region of the Large Magellanic Cloud. Radial velocities are provided for 307 apparently single stars, and for 99 targets with radial-velocity variations which are consistent with them being spectroscopic binaries. We investigate the spatial distribution of the radial velocities across the 30 Dor region, and use the results to identify candidate runaway stars. Excluding potential runaways and members of two older clusters in the survey region (SL 639 and Hodge 301), we determine a systemic velocity for 30 Dor of 271.6 ± 12.2 km s^-1 from 273 presumed single stars. Employing a 3σ criterion we identify nine candidate runaway stars (2.9% of the single stars with radial-velocity estimates). The projected rotational velocities of the candidate runaways appear to be significantly different to those of the full B-type sample, with a strong preference for either large (≥345 km s^-1) or small (≤65 km s^-1) rotational velocities. Of the candidate runaways, VFTS 358 (classified B0.5: V) has the largest differential radial velocity (-106.9 ± 16.2 km s^-1), and a preliminary atmospheric analysis finds a significantly enriched nitrogen abundance of 12 + log (N/H) ≳ 8.5. Combined with a large rotational velocity (υ_e sin i = 345 ± 22 km s^-1), this is suggestive of past binary interaction for this star.

Rotational velocities of single and binary O-type stars in the Tarantula Nebula

Rotation is a key parameter in the evolution of massive stars, affecting their evolution, chemical yields, ionizing photon budget, and final fate. We determined the projected rotational velocity, υ e sin i, of ~330 O-type objects, i.e. ~210 spectroscopic single stars and ~110 primaries in binary systems, in the Tarantula nebula or 30 Doradus (30 Dor) region. The observations were taken using VLT/FLAMES and constitute the largest homogeneous dataset of multi-epoch spectroscopy of O-type stars currently available. The most distinctive feature of the υ e sin i distributions of the presumed-single stars and primaries in 30 Dor is a low-velocity peak at around 100 km s^-1. Stellar winds are not expected to have spun-down the bulk of the stars significantly since their arrival on the main sequence and therefore the peak in the single star sample is likely to represent the outcome of the formation process. Whereas the spin distribution of presumed-single stars shows a well developed tail of stars rotating more rapidly than 300 km s^-1, the sample of primaries does not feature such a high-velocity tail. The tail of the presumed-single star distribution is attributed for the most part - and could potentially be completely due - to spun-up binary products that appear as single stars or that have merged. This would be consistent with the lack of such post-interaction products in the binary sample, that is expected to be dominated by pre-interaction systems. The peak in this distribution is broader and is shifted toward somewhat higher spin rates compared to the distribution of presumed-single stars. Systems displaying large radial velocity variations, typical for short period systems, appear mostly responsible for these differences.

Non-radial instabilities and progenitor asphericities in core-collapse supernovae

Since core-collapse supernova simulations still struggle to produce robust neutrino-driven explosions in 3D, it has been proposed that asphericities caused by convection in the progenitor might facilitate shock revival by boosting the activity of non-radial hydrodynamic instabilities in the post-shock region. We investigate this scenario in depth using 42 relativistic 2D simulations with multigroup neutrino transport to examine the effects of velocity and density perturbations in the progenitor for different perturbation geometries that obey fundamental physical constraints (like the anelastic condition). As a framework for analysing our results, we introduce semi-empirical scaling laws relating neutrino heating, average turbulent velocities in the gain region, and the shock deformation in the saturation limit of non-radial instabilities. The squared turbulent Mach number, 〈Ma2〉, reflects the violence of aspherical motions in the gain layer, and explosive runaway occurs for 〈Ma2〉 ≳ 0.3, corresponding to a reduction of the critical neutrino luminosity by ∼25∼25 per cent compared to 1D. In the light of this theory, progenitor asphericities aid shock revival mainly by creating anisotropic mass flux on to the shock: differential infall efficiently converts velocity perturbations in the progenitor into density perturbations δρ/ρ at the shock of the order of the initial convective Mach number Maprog. The anisotropic mass flux and ram pressure deform the shock and thereby amplify post-shock turbulence. Large-scale (ℓ = 2, ℓ = 1) modes prove most conducive to shock revival, whereas small-scale perturbations require unrealistically high convective Mach numbers. Initial density perturbations in the progenitor are only of the order of Ma2progMaprog2 and therefore play a subdominant role.

The large- and small-scale Ca IIK structure of the Milky Way from observations of Galactic and Magellanic sightlines

Aims. The large and small-scale (pc) structure of the Galactic interstellar medium can be investigated by utilising spectra of early-type stellar probes of known distances in the same region of the sky. This paper determines the variation in line strength of Ca ii at 3933.661 Å as a function of probe separation for a large sample of stars, including a number of sightlines in the Magellanic Clouds. 

Methods. FLAMES-GIRAFFE data taken with the Very Large Telescope towards early-type stars in 3 Galactic and 4 Magellanic open clusters in Ca ii are used to obtain the velocity, equivalent width, column density, and line width of interstellar Galactic calcium for a total of 657 stars, of which 443 are Magellanic Cloud sightlines. In each cluster there are between 43 and 111 stars observed. Additionally, FEROS and UVES Ca ii K and Na i D spectra of 21 Galactic and 154 Magellanic early-type stars are presented and combined with data from the literature to study the calcium column density - parallax relationship. 

Results. For the four Magellanic clusters studied with FLAMES, the strength of the Galactic interstellar Ca ii K equivalent width on transverse scales from ∼0.05-9 pc is found to vary by factors of ∼1.8-3.0, corresponding to column density variations of ∼0.3-0.5 dex in the optically-thin approximation. Using FLAMES, FEROS, and UVES archive spectra, the minimum and maximum reduced equivalent widths for Milky Way gas are found to lie in the range ∼35-125 mÅ and ∼30-160 mÅ for Ca ii K and Na i D, respectively. The range is consistent with a previously published simple model of the interstellar medium consisting of spherical cloudlets of filling factor ∼0.3, although other geometries are not ruled out. Finally, the derived functional form for parallax (π) and Ca ii column density (N_CaII) is found to be π(mas) = 1 / (2.39 × 10^-13 × N_CaII (cm^-2) + 0.11). Our derived parallax is ∼25 per cent lower than predicted by Megier et al. (2009, A&A, 507, 833) at a distance of ∼100 pc and ∼15 percent lower at a distance of ∼200 pc, reflecting inhomogeneity in the Ca ii distribution in the different sightlines studied.

Discovery of WASP-113b and WASP-114b, two inflated hot-Jupiters with contrasting densities

We present the discovery and characterisation of the exoplanets WASP-113b and WASP-114b by the WASP survey, SOPHIE and CORALIE. The planetary nature of the systems was established by performing follow-up photometric and spectroscopic observations. The follow-up data were combined with the WASP-photometry and analysed with an MCMC code to obtain system parameters. The host stars WASP-113 and WASP-114 are very similar. They are both early G-type stars with an effective temperature of ~5900K, [Fe/H] ~0.12 and T_{eff} ~4.1 dex. However, WASP-113 is older than WASP-114. Although the planetary companions have similar radii, WASP-114b is almost 4 times heavier than WASP-113b. WASP-113b has a mass of 0.48 M_{Jup} and an orbital period of ~4.5 days; WASP-114b has a mass of 1.77 M_{Jup} and an orbital period of ~1.5 days. Both planets have inflated radii, in particular WASP-113 with a radius anomaly of Re=0.35. The high scale height of WASP-113b (~950 km ) makes it a good target for follow-up atmospheric observations.

The Outer Solar Systems Origins Survey. I. Design and First-quarter Discoveries

We report the discovery, tracking, and detection circumstances for 85 trans-Neptunian objects (TNOs) from the first 42 deg2 of the Outer Solar System Origins Survey. This ongoing r-band solar system survey uses the 0.9 deg2 field of view MegaPrime camera on the 3.6 m Canada–France–Hawaii Telescope. Our orbital elements for these TNOs are precise to a fractional semimajor axis uncertainty <0.1%. We achieve this precision in just two oppositions, as compared to the normal three to five oppositions, via a dense observing cadence and innovative astrometric technique. These discoveries are free of ephemeris bias, a first for large trans-Neptunian surveys. We also provide the necessary information to enable models of TNO orbital distributions to be tested against our TNO sample. We confirm the existence of a cold "kernel" of objects within the main cold classical Kuiper Belt and infer the existence of an extension of the "stirred" cold classical Kuiper Belt to at least several au beyond the 2:1 mean motion resonance with Neptune. We find that the population model of Petit et al. remains a plausible representation of the Kuiper Belt. The full survey, to be completed in 2017, will provide an exquisitely characterized sample of important resonant TNO populations, ideal for testing models of giant planet migration during the early history of the solar system.

Molecular Hydrogen Emission from Protoplanetary Disks

We have modeled self-consistently the density and temperature profiles of gas and dust in protoplanetary disks, taking into account irradiation from a central star. Making use of this physical structure, we have calculated the level populations of molecular hydrogen and the line emission from the disks. As a result, we can reproduce the observed strong line spectra of molecular hydrogen from protoplanetary disks, both in the ultraviolet (UV) and the near-infrared, but only if the central star has a strong UV excess radiation.

Molecular Distributions in the Inner Regions of Protostellar Disks

The distributions of molecules in the inner regions of a protostellar disk are presented. These were calculated using an uncoupled chemical/dynamical model, with a numerical integration of the vertical disk structure. A comparison between models with and without the effects of X-ray ionisation is made, and molecules are identified which are good tracers of the ionisation level in this part of the disk, notably CN and C_2H. In the region considered in this paper (r