Complex organic molecules in protoplanetary disks

Context. Protoplanetary disks are vital objects in star and planet formation, possessing all the material, gas and dust, which may form a planetary system orbiting the new star. Small, simple molecules have traditionally been detected in protoplanetary disks; however, in the ALMA era, we expect the molecular inventory of protoplanetary disks to significantly increase.

Aims. We investigate the synthesis of complex organic molecules (COMs) in protoplanetary disks to put constraints on the achievable chemical complexity and to predict species and transitions which may be observable with ALMA.

Methods. We have coupled a 2D steady-state physical model of a protoplanetary disk around a typical T Tauri star with a large gas-grain chemical network including COMs. We compare the resulting column densities with those derived from observations and perform ray-tracing calculations to predict line spectra. We compare the synthesised line intensities with current observations and determine those COMs which may be observable in nearby objects. We also compare the predicted grain-surface abundances with those derived from cometary comae observations.

Results. We find COMs are efficiently formed in the disk midplane via grain-surface chemical reactions, reaching peak grain-surface fractional abundances similar to 10(-6)-10(-4) that of the H nuclei number density. COMs formed on grain surfaces are returned to the gas phase via non-thermal desorption; however, gas-phase species reach lower fractional abundances than their grain-surface equivalents, similar to 10(-12)-10(-7). Including the irradiation of grain mantle material helps build further complexity in the ice through the replenishment of grain-surface radicals which take part in further grain-surface reactions. There is reasonable agreement with several line transitions of H2CO observed towards T Tauri star-disk systems. There is poor agreement with HC3(N) lines observed towards LkCa 15 and GO Tau and we discuss possible explanations for these discrepancies. The synthesised line intensities for CH3OH are consistent with upper limits determined towards all sources. Our models suggest CH3OH should be readily observable in nearby protoplanetary disks with ALMA; however, detection of more complex species may prove challenging, even with ALMA "Full Science" capabilities. Our grain-surface abundances are consistent with those derived from cometary comae observations providing additional evidence for the hypothesis that comets (and other planetesimals) formed via the coagulation of icy grains in the Sun's natal disk.

Three-dimensional pure deflagration models with nucleosynthesis and synthetic observables for type ia supernovae

We investigate whether pure deflagration models ofChandrasekhar-mass carbon-oxygen white dwarf stars can account for one or more subclass of the observed population of Type Ia supernova (SN Ia) explosions. We compute a set of 3D full-star hydrodynamic explosion models, in which the deflagration strength is parametrized using the multispot ignition approach. For each model, we calculate detailed nucleosynthesis yields in a post-processing step with a 384 nuclide nuclear network. We also compute synthetic observables with our 3D Monte Carlo radiative transfer code for comparison with observations. For weak and intermediate deflagration strengths (energy release E {less-than or approximate} 1.1 × 10 erg), we find that the explosion leaves behind a bound remnant enriched with 3 to 10 per cent (by mass) of deflagration ashes. However, we do not obtain the large kick velocities recently reported in the literature. We find that weak deflagrations with E ~ 0.5 × 10 erg fit well both the light curves and spectra of 2002cx-like SNe Ia, and models with even lower explosion energies could explain some of the fainter members of this subclass. By comparing our synthetic observables with the properties of SNe Ia, we can exclude the brightest, most vigorously ignited models as candidates for any observed class of SN Ia: their B-V colours deviate significantly from both normal and 2002cx-like SNe Ia and they are too bright to be candidates for other subclasses.

The spotty donor star in the X-ray transient Cen X-4

We accurately determine the fundamental system parameters of the neutron star X-ray transient Cen X-4 solely using phase-resolved high-resolution UV-Visual Echelle Spectrograph spectroscopy. We first determine the radial-velocity curve of the secondary star and then model the shape of the phase-resolved absorption line profiles using an X-ray binary model. The model computes the exact rotationally broadened, phase-resolved spectrum and does not depend on assumptions about the rotation profile, limb-darkening coefficients and the effects of contamination from an accretion disc. We determine the secondary star-to-neutron star binary mass ratio to be 0.1755 ± 0.0025, which is an order of magnitude more accurate than previous estimates. We also constrain the inclination angle to be 32^{+8}_{-2} degrees. Combining these values with the results of the radial-velocity study gives a neutron star mass of 1.94^{+0.37}_{-0.85}M⊙ consistent with previous estimates. Finally, we perform the first Roche tomography reconstruction of the secondary star in an X-ray binary. The tomogram reveals surface inhomogeneities that are due to the presence of cool starspots. A large cool polar spot, similar to that seen in Doppler images of rapidly rotating isolated stars, is present on the Northern hemisphere of the K7 secondary star and we estimate that ~4 percent of the total surface area of the donor star is covered with spots.This evidence for starspots supports the idea that magnetic braking plays an important role in the evolution of low-mass X-ray binaries.

A window on exoplanet dynamical histories: Rossiter-McLaughlin observations of WASP-13b and WASP-32b

We present Rossiter-McLaughlin observations of WASP-13b and WASP-32b and determine the sky-projected angle between the normal of the planetary orbit and the stellar rotation axis (λ). WASP-13b and WASP-32b both have prograde orbits and are consistent with alignment with measured sky-projected angles of λ =8°^{+13}_{-12} and λ =-2°^{+17}_{-19}, respectively. Both WASP-13 and WASP-32 have Teff < 6250 K, and therefore, these systems support the general trend that aligned planetary systems are preferentially found orbiting cool host stars. A Lomb-Scargle periodogram analysis was carried out on archival SuperWASP data for both systems. A statistically significant stellar rotation period detection (above 99.9 per cent confidence) was identified for the WASP-32 system with Prot =11.6 ± 1.0 days. This rotation period is in agreement with the predicted stellar rotation period calculated from the stellar radius, R*, and vsin i if a stellar inclination of i* =90° is assumed. With the determined rotation period, the true 3D angle between the stellar rotation axis and the planetary orbit, ψ, was found to be ψ = 11° ± 14°. We conclude with a discussion on the alignment of systems around cool host stars with Teff < 6150 K by calculating the tidal dissipation time-scale. We find that systems with short tidal dissipation time-scales are preferentially aligned and systems with long tidal dissipation time-scales have a broad range of obliquities.

Electron impact excitation of Mg VIII: Collision strengths, transition probabilities and theoretical EUV and soft X-ray line intensities for Mg VIII

Context: Mg VIII emission lines are observed in a range of astronomical objects such as the Sun, other cool stars and in the coronal line region of Seyfert galaxies. Under coronal conditions Mg VIII emits strongly in the extreme ultraviolet (EUV) and soft X-ray spectral regions which makes it an ideal ion for plasma diagnostics.

Aims. Two theoretical atomic models, consisting of 125 fine structure levels, are developed for the Mg VIII ion. The 125 levels arise from the 2s(2)2p, 2s(2)p2, 2p(3), 2s(2)3s, 2s(2)3p, 2s(2)3d, 2s2p3s, 2s2p3p, 2s2p3d, 2p(2)3s, 2p(2)3p and 2p(2)3d configurations. Electron impact excitation collision strengths and radiative transition probabilities are calculated for both Mg VIII models, compared with existing data, and the best model selected to generate a set of theoretical emission line intensities. The EUV lines, covering 312-790 angstrom, are compared with existing solar spectra (SERTS-89 and SUMER), while the soft X-ray transitions (69-97 angstrom) are examined for potential density diagnostic line ratios and also compared with the limited available solar and stellar observational data.

Methods. The R-matrix codes Breit-Pauli RMATRXI and RMATRXII are utilised, along with the PSTGF code, to calculate the collision strengths for two Mg VIII models. Collision strengths are averaged over a Maxwellian distribution to produce the corresponding effective collision strengths for use in astrophysical applications. Transition probabilities are also calculated using the CIV3 atomic structure code. The best data are then incorporated into the modelling code CLOUDY and line intensities generated for a range of electron temperatures and densities appropriate to solar and stellar coronal plasmas.

Results. The present effective collision strengths are compared with two previous calculations. Good levels of agreement are found with the most recent, but there are large differences with the other for forbidden transitions. The resulting line intensities compare favourably with the observed values from the SERTS-89 and SUMER spectra. Theoretical soft X-ray emission lines are presented and several density diagnostic line ratios examined, which are in reasonable agreement with the limited observational data available.

The Kepler-10 Planetary System Revisited by HARPS-N: A Hot Rocky World and a Solid Neptune-Mass Planet

Kepler-10b was the first rocky planet detected by the Kepler satellite and confirmed with radial velocity follow-up observations from Keck-HIRES. The mass of the planet was measured with a precision of around 30%, which was
insufficient to constrain models of its internal structure and composition in detail. In addition to Kepler-10b, a second planet transiting the same star with a period of 45 days was statistically validated, but the radial velocities were only
good enough to set an upper limit of 20 M⊕ for the mass of Kepler-10c. To improve the precision on the mass for planet b, the HARPS-N Collaboration decided to observe Kepler-10 intensively with the HARPS-N spectrograph
on the Telescopio Nazionale Galileo on La Palma. In total, 148 high-quality radial-velocity measurements were obtained over two observing seasons. These new data allow us to improve the precision of the mass determination for Kepler-10b to 15%. With a mass of 3.33 ± 0.49 M⊕ and an updated radius of 1.47+0.03 −0.02 R⊕, Kepler-10b has a density of 5.8 ± 0.8 g cm−3, very close to the value predicted by models with the same internal structure and composition as the Earth. We were also able to determine a mass for the 45-day period planet Kepler-10c, with an even better precision of 11%. With a mass of 17.2 ± 1.9 M⊕ and radius of 2.35+0.09 −0.04 R⊕, Kepler-10c has a density of 7.1 ± 1.0 g cm−3. Kepler-10c appears to be the first strong evidence of a class of more massive solid planets with longer orbital periods

Roche tomography of cataclysmic variables - VI. Differential rotation of AE Aqr - not tidally locked!

We present Roche tomograms of the K4V secondary star in the cataclysmic variable AE Aqr, reconstructed from two data sets taken 9 d apart, and measure the differential rotation of the stellar surface. The tomograms show many large, cool starspots, including a large high-latitude spot and a prominent appendage down the trailing hemisphere. We find two distinct bands of spots around 22° and 43° latitude, and estimate a spot coverage of 15.4-17 per cent on the Northern hemisphere. Assuming a solar-like differential rotation law, the differential rotation of AE Aqr was measured using two different techniques. The first method yields an equator-pole lap time of 269 d and the second yields a lap time of 262 d. This shows that the star is not fully tidally locked, as was previously assumed for CVs, but has a co-rotation latitude of ˜40°. We discuss the implications that these observations have on stellar dynamo theory, as well as the impact that spot traversal across the L1 point may have on accretion rates in CVs as well as some of their other observed properties. The entropy landscape technique was applied to determine the system parameters of AE Aqr. For the two independent data sets, we find M1 = 1.20 and 1.17 M⊙, M2 = 0.81 and 0.78 M⊙, and orbital inclinations of 50° to 51° at optimal systemic velocities of γ = -64.7 and -62.9 km s-1.

Understanding Astrophysical Noise from Stellar Surface Magneto-Convection

To obtain cm/s precision, stellar surface magneto-convection must be disentangled from observed radial velocities (RVs). In order to understand and remove the convective signature, we create Sun-as-a-star model observations based on a 3D magnetohydrodynamic solar simulation. From these Sun-as-a-star model observations, we find several line characteristics are correlated with the induced RV shifts. The aim of this campaign is to feed directly into future high precision RV studies, such as the search for habitable, rocky worlds, with forthcoming spectrographs such as ESPRESSO.

Herschel key program HERITAGE:A far-infrared source catalog for the Magellanic Clouds

Observations from the HERschel Inventory of the Agents of Galaxy Evolution (HERITAGE ) have been used to identify dusty populations of sources in the Large and Small Magellanic Clouds (LMC and SMC). We conducted the study using the HERITAGE catalogs of point sources available from the Herschel Science Center from both the Photodetector Array Camera and Spectrometer (PACS; 100 and 160 μm) and Spectral and Photometric Imaging Receiver (SPIRE; 250, 350, and 500 μm) cameras. These catalogs are matched to each other to create a Herschel band-merged catalog and then further matched to archival Spitzer IRAC and MIPS catalogs from the Spitzer Surveying the Agents of Galaxy Evolution (SAGE) and SAGE-SMC surveys to create single mid- to far-infrared (far-IR) point source catalogs that span the wavelength range from 3.6 to 500 μm. There are 35,322 unique sources in the LMC and 7503 in the SMC. To be bright in the FIR, a source must be very dusty, and so the sources in the HERITAGE catalogs represent the dustiest populations of sources. The brightest HERITAGE sources are dominated by young stellar objects (YSOs), and the dimmest by background galaxies. We identify the sources most likely to be background galaxies by first considering their morphology (distant galaxies are point-like at the resolution of Herschel) and then comparing the flux distribution to that of the Herschel Astrophysical Terahertz Large Area Survey (ATLAS ) survey of galaxies. We find a total of 9745 background galaxy candidates in the LMC HERITAGE images and 5111 in the SMC images, in agreement with the number predicted by extrapolating from the ATLAS flux distribution. The majority of the Magellanic Cloud-residing sources are either very young, embedded forming stars or dusty clumps of the interstellar medium. Using the presence of 24 μm emission as a tracer of star formation, we identify 3518 YSO candidates in the LMC and 663 in the SMC. There are far fewer far-IR bright YSOs in the SMC than the LMC due to both the SMC's smaller size and its lower dust content. The YSO candidate lists may be contaminated at low flux levels by background galaxies, and so we differentiate between sources with a high ("probable") and moderate ("possible ") likelihood of being a YSO. There are 2493/425 probable YSO candidates in the LMC/SMC. Approximately 73% of the Herschel YSO candidates are newly identified in the LMC, and 35% in the SMC. We further identify a small population of dusty objects in the late stages of stellar evolution including extreme and post-asymptotic giant branch, planetary nebulae, and supernova remnants. These populations are identified by matching the HERITAGE catalogs to lists of previously identified objects in the literature. Approximately half of the LMC sources and one quarter of the SMC sources are too faint to obtain accurate ample FIR photometry and are unclassified.

Superluminous supernovae from PESSTO

We present optical spectra and light curves for three hydrogen-poor superluminous supernovae followed by the Public ESO Spectroscopic Survey of Transient Objects (PESSTO). Time series spectroscopy from a fewdays aftermaximum light to 100 d later shows them to be fairly typical of this class, with spectra dominated by Ca II, MgII, FeII, and Si II, which evolve slowly over most of the post-peak photospheric phase. We determine bolometric light curves and apply simple fitting tools, based on the diffusion of energy input by magnetar spin-down, Ni-56 decay, and collision of the ejecta with an opaque circumstellar shell. We investigate how the heterogeneous light curves of our sample (combined with others from the literature) can help to constrain the possible mechanisms behind these events. We have followed these events to beyond 100-200 d after peak, to disentangle host galaxy light from fading supernova flux and to differentiate between the models, which predict diverse behaviour at this phase. Models powered by radioactivity require unrealistic parameters to reproduce the observed light curves, as found by previous studies. Both magnetar heating and circumstellar interaction still appear to be viable candidates. A large diversity is emerging in observed tail-phase luminosities, with magnetar models failing in some cases to predict the rapid drop in flux. This would suggest either that magnetars are not responsible, or that the X-ray flux from the magnetar wind is not fully trapped. The light curve of one object shows a distinct rebrightening at around 100 d after maximum light. We argue that this could result either from multiple shells of circumstellar material, or from a magnetar ionization front breaking out of the ejecta.

A Bayesian method for detecting stellar flares

We present a Bayesian-odds-ratio-based algorithm for detecting stellar flares in light-curve data. We assume flares are described by a model in which there is a rapid rise with a half-Gaussian profile, followed by an exponential decay. Our signal model also contains a polynomial background model required to fit underlying light-curve variations in the data, which could otherwise partially mimic a flare. We characterize the false alarm probability and efficiency of this method under the assumption that any unmodelled noise in the data is Gaussian, and compare it with a simpler thresholding method based on that used in Walkowicz et al. We find our method has a significant increase in detection efficiency for low signal-to-noise ratio (S/N) flares. For a conservative false alarm probability our method can detect 95 per cent of flares with S/N less than 20, as compared to S/N of 25 for the simpler method. We also test how well the assumption of Gaussian noise holds by applying the method to a selection of 'quiet' Kepler stars. As an example we have applied our method to a selection of stars in Kepler Quarter 1 data. The method finds 687 flaring stars with a total of 1873 flares after vetos have been applied. For these flares we have made preliminary characterizations of their durations and and S/N.

The Well-Aligned Orbit of Wasp-84b: Evidence for Disk Migration of a Hot Jupiter

We report the sky-projected orbital obliquity (spin–orbit angle) of WASP-84 b, a 0.69MJup planet in an 8.52 day orbit around a G9V/K0V star, to be λ = −0.3 ± 1.7°. We obtain a true obliquity of ψ = 17.3 ± 7.7° from a measurement of the inclination of the stellar spin axis with respect to the sky plane. Due to the young age and the weak tidal forcing of the system, we suggest that the orbit of WASP-84b is unlikely to have both realigned and circularized from the misaligned and/or eccentric orbit likely to have arisen from high-eccentricity migration. Therefore we conclude that the planet probably migrated via interaction with the protoplanetary disk. This would make it the first “hot Jupiter” (P d < 10 ) to have been shown to have migrated via this pathway. Further, we argue that the distribution of obliquities for planets orbiting cool stars (Teff < 6250 K) suggests that high-eccentricity migration is an important pathway for the formation of short-orbit, giant planets.

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 mu m data and improved imaging quality at 100 and 160 mu m compared to previous observations in 2010. Our Herschel spectra show only weak CO line emission, and provide an upper limit for the 63 mu m [O-I] line flux, eliminating the possibility that line contaminations distort the previously estimated dustmass. The far-infrared spectral energy distribution (SED) is well fitted by thermal emission from cold dust. The newly measured 70 mu m flux constrains the dust temperature, limiting it to nearly a single temperature. The far-infrared emission can be fitted by 0.5 +/- 0.1M(circle dot) 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.3M(circle dot) of amorphous carbon and 0.5M(circle dot) of silicates, totalling 0.8M(circle dot) 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.

Characterizing K2 Planet Discoveries: A super-Earth transiting the bright K-dwarf HIP 116454

We report the first planet discovery from the two-wheeled Kepler (K2) mission: HIP 116454 b. The host star HIP 116454 is a bright (V = 10.1, K = 8.0) K1 dwarf with high proper motion and a parallax-based distance of 55.2 +/- 5.4 pc. Based on high-resolution optical spectroscopy, we find that the host star is metal-poor with [Fe/H]= -0.16 +/- 0.08 and has a radius R-star = 0.716 +/- 0.024 R-circle dot and mass M-star = 0.775 +/- 0.027M(circle dot). The star was observed by the Kepler spacecraft during its Two-Wheeled Concept Engineering Test in 2014 February. During the 9 days of observations, K2 observed a single transit event. Using a new K2 photometric analysis technique, we are able to correct small telescope drifts and recover the observed transit at high confidence, corresponding to a planetary radius of R-p = 2.53 +/- 0.18 R-circle plus. Radial velocity observations with the HARPS-N spectrograph reveal a 11.82 +/- 1.33 M-circle plus planet in a 9.1 day orbit, consistent with the transit depth, duration, and ephemeris. Follow-up photometric measurements from the MOST satellite confirm the transit observed in the K2 photometry and provide a refined ephemeris, making HIP 116454 b amenable for future follow-up observations of this latest addition to the growing population of transiting super-Earths around nearby, bright stars.

Possible binary progenitors for the Type Ib supernova iPTF13bvn

Cao et al. reported a possible progenitor detection for the Type Ib supernovae iPTF13bvn for the first time. We find that the progenitor is in fact brighter than the magnitudes previously reported by approximately 0.7-0.2 mag with a larger error in the bluer filters. We compare our new magnitudes to our large set of binary evolution models and find that many binary models with initial masses in the range of 10-20M(circle dot) match this new photometry and other constraints suggested from analysing the supernova. In addition, these lower mass stars retain more helium at the end of the model evolution indicating that they are likely to be observed as Type Ib supernovae rather than their more massive, Wolf-Rayet counter parts. We are able to rule out typical Wolf-Rayet models as the progenitor because their ejecta masses are too high and they do not fit the observed SED unless they have a massive companion which is the observed source at the supernova location. Therefore only late-time observations of the location will truly confirm if the progenitor was a helium giant and not a Wolf-Rayet star.