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.

The main-belt comets: The Pan-STARRS1 perspective

We analyze a set of 760 475 observations of 333 026 unique main-belt objects obtained by the Pan-STARRS1(PS1) survey telescope between 2012 May 20 and 2013 November 9, a period during which PS1 discoveredtwo main-belt comets, P/2012 T1 (PANSTARRS) and P/2013 R3 (Catalina-PANSTARRS). PS1 cometdetection procedures currently consist of the comparison of the point spread functions (PSFs) of movingobjects to those of reference stars, and the flagging of objects that show anomalously large radial PSFwidths for human evaluation and possible observational follow-up. Based on the number of missed discoveryopportunities among comets discovered by other observers, we estimate an upper limit comet discoveryefficiency rate of 70% for PS1. Additional analyses that could improve comet discovery yields infuture surveys include linear PSF analysis, modeling of trailed stellar PSFs for comparison to trailed movingobject PSFs, searches for azimuthally localized activity, comparison of point-source-optimized photometryto extended-source-optimized photometry, searches for photometric excesses in objects withknown absolute magnitudes, and crowd-sourcing. Analysis of the discovery statistics of the PS1 surveyindicates an expected fraction of 59 MBCs per 106 outer main-belt asteroids (corresponding to a totalexpected population of 140 MBCs among the outer main-belt asteroid population with absolute magnitudesof 12 < HV < 19:5), and a 95% confidence upper limit of 96 MBCs per 106 outer main-belt asteroids(corresponding to a total of 230 MBCs), assuming a detection efficiency of 50%. We note howeverthat significantly more sensitive future surveys (particularly those utilizing larger aperture telescopes)could detect many more MBCs than estimated here. Examination of the orbital element distribution ofall known MBCs reveals an excess of high eccentricities (0:1 < e < 0:3) relative to the background asteroidpopulation. Theoretical calculations show that, given these eccentricities, the sublimation rate for atypical MBC is orders of magnitude larger at perihelion than at aphelion, providing a plausible physicalexplanation for the observed behavior of MBCs peaking in observed activity strength near perihelion.These results indicate that the overall rate of mantle growth should be slow, consistent with observationalevidence that MBC activity can be sustained over multiple orbit passages.

SN 2012ec: mass of the progenitor from PESSTO follow-up of the photospheric phase

We present the results of a photometric and spectroscopic monitoring campaign of SN 2012ec, which exploded in the spiral galaxy NGC 1084, during the photospheric phase. The photometric light curve exhibits a plateau with luminosity L = 0.9 x 10(42) erg s(-1) and duration similar to 90 d, which is somewhat shorter than standard Type II-P supernovae (SNe). We estimate the nickel mass M(Ni-56) = 0.040 +/- 0.015 M-circle dot from the luminosity at the beginning of the radioactive tail of the light curve. The explosion parameters of SN 2012ec were estimated from the comparison of the bolometric light curve and the observed temperature and velocity evolution of the ejecta with predictions from hydrodynamical models. We derived an envelope mass of 12.6 M-circle dot, an initial progenitor radius of 1.6 x 10(13) cm and an explosion energy of 1.2 foe. These estimates agree with an independent study of the progenitor star identified in pre-explosion images, for which an initial mass of M = 14-22 M-circle dot was determined. We have applied the same analysis to two other Type II-P SNe (SNe 2012aw and 2012A), and carried out a comparison with the properties of SN 2012ec derived in this paper. We find a reasonable agreement between the masses of the progenitors obtained from pre-explosion images and masses derived from hydrodynamical models. We estimate the distance to SN 2012ec with the standardized candle method (SCM) and compare it with other estimates based on other primary and secondary indicators. SNe 2012A, 2012aw and 2012ec all follow the standard relations for the SCM for the use of Type II-P SNe as distance indicators.

Supersolar Ni/Fe production in the Type IIP SN 2012ec

SN 2012ec is a Type IIP supernova (SN) with a progenitor detection and comprehensive photospheric phase observational coverage. Here, we present Very Large Telescope and Public ESO Spectroscopic Survey of Transient Objects observations of this SN in the nebular phase. We model the nebular [O I] lambda lambda 6300, 6364 lines and find their strength to suggest a progenitor main-sequence mass of 13-15 M-circle dot. SN2012ec is unique among hydrogen-rich SNe in showing a distinct line of stable nickel [Ni II] lambda 7378. This line is produced by Ni-58, a nuclear burning ash whose abundance is a sensitive tracer of explosive burning conditions. Using spectral synthesis modelling, we use the relative strengths of [Ni II] lambda 7378 and [Fe II] lambda 7155 (the progenitor of which is Ni-56) to derive a Ni/Fe production ratio of 0.20 +/- 0.07 (by mass), which is a factor 3.4 +/- 1.2 times the solar value. High production of stable nickel is confirmed by a strong [Ni II] 1.939 mu m line. This is the third reported case of a core-collapse SN producing a Ni/Fe ratio far above the solar value, which has implications for core-collapse explosion theory and galactic chemical evolution models.

Pre-Main-Sequence Accretion in the Low Metallicity Galactic Star-forming Region Sh 2-284

We present optical spectra of pre-main-sequence (PMS) candidates around the Ha region taken with the Southern African Large Telescope in the low metallicity (Z) Galactic region Sh 2-284, which includes the open cluster Dolidze 25 with an atypical low metallicity of Z similar to 1/5 Z(circle dot). It has been suggested on the basis of both theory and observations that PMS mass-accretion rates, (M) over dot(acc), are a function of Z. We present the first sample of spectroscopic estimates of mass-accretion rates for PMS stars in any low-Z star-forming region. Our data set was enlarged with literature data of H alpha emission in intermediate-resolution R-band spectroscopy. Our total sample includes 24 objects spanning a mass range between 1 and 2 M-circle dot and with a median age of approximately 3.5 Myr. The vast majority (21 out of 24) show evidence for a circumstellar disk on the basis of Two Micron All Sky Survey and Spitzer infrared photometry. We find (M) over dot(acc) in the 1-2 M-circle dot interval to depend quasi-quadratically on stellarmass, with (M) over dot(acc) proportional to M-*(2.4 +/- 0.35), and inversely with stellar age, with (M) over dot(acc) proportional to t(*)(-0.7 +/- 0.4). Furthermore, we compare our spectroscopic (M) over dot(acc) measurements with solar Z Galactic PMS stars in the same mass range, but, surprisingly find no evidence for a systematic change in (M) over dot(acc) with Z. We show that literature accretion-rate studies are influenced by detection limits, and we suggest that (M) over dot(acc) may be controlled by factors other than Z(*), M-*, and age.

Detecting the progenitors of core collapse supernovae

The masses and the evolutionary states of the progenitors of core-collapse supernovae are not well constrained by direct observations. Stellar evolution theory generally predicts that massive stars with initial masses less than about 30M_sol should undergo core-collapse when they are cool M-type supergiants. However the only two detections of a SN progenitor before explosion are SN1987A and SN1993J, and neither of these was an M-type supergiant. Attempting to identify the progenitors of supernovae is a difficult task, as precisely predicting the time of explosion of a massive star is impossible for obvious reasons. There are several different types of supernovae which have different spectral and photometric evolution, and how exactly these are related to the evolutionary states of the progenitor stars is not currently known. I will describe a novel project which may allow the direct identification of core-collapse supernovae progenitors on pre-explosion images of resolved, nearby galaxies. This project is now possible with the excellent image archives maintained by several facilities and will be enhanced by the new initiatives to create Virtual Observatories, the earliest of which ASTROVIRTEL is already producing results.

The Decaying Long-period Oscillation of a Stellar Megaflare

We analyze and interpret the oscillatory signal in the decay phase of the U-band light curve of a stellar megaflare observed on 2009 January 16 on the dM4.5e star YZ CMi. The oscillation is well approximated by an exponentially decaying harmonic function. The period of the oscillation is found to be 32 minutes, the decay time about 46 minutes, and the relative amplitude 15%. As this observational signature is typical of the longitudinal oscillations observed in solar flares at extreme ultraviolet and radio wavelengths, associated with standing slow magnetoacoustic waves, we suggest that this megaflare may be of a similar nature. In this scenario, macroscopic variations of the plasma parameters in the oscillations modulate the ejection of non-thermal electrons. The phase speed of the longitudinal (slow magnetoacoustic) waves in the flaring loop or arcade, the tube speed, of about 230 km s-1 would require a loop length of about 200 Mm. Other mechanisms, such as standing kink oscillations, are also considered.

Transit timing variations in WASP-10b induced by stellar activity

The hot-JupiterWASP-10bwas reported by Maciejewski et al. to showtransit timing variations (TTVs) with an amplitude of ~3.5 min. These authors proposed that the observed TTVs were caused by a 0.1M_Jup perturbing companion with an orbital period of ~5.23 d, and hence, close to the outer 5:3 mean-motion resonance with WASP-10b. To test this scenario, we present eight new transit light curves of WASP-10b obtained with the Faulkes Telescope North and the Liverpool Telescope. The new light curves, together with 22 previously published ones, were modelled with a Markov Chain Monte Carlo transit fitting code. Transit depth differences reported forWASP-10b are thought to be due to starspot-induced brightness modulation of the host star. Assuming the star is brighter at the activity minimum, we favour a small planetary radius. We find Rp = 1.039^+0.043 _-0.049R_Jup in agreement with Johnson et al. and Maciejewski et al. Recent studies find no evidence for a significant eccentricity in this system. We present consistent system parameters for a circular orbit and refine the orbital ephemeris ofWASP-10b. Our homogeneously derived transit times do not support the previous claimed TTV signal, which was strongly dependent on two previously published transits that have been incorrectly normalized. Nevertheless, a linear ephemeris is not a statistically good fit to the transit times of WASP-10b. We show that the observed transit time variations are due to spot occultation features or systematics. We discuss and exemplify the effects of occultation spot features in the measured transit times and show that despite spot occultation during egress and ingress being difficult to distinguish in the transit light curves, they have a significant effect in the measured transit times. We conclude that if we account for spot features, the transit times of WASP-10b are consistent with a linear ephemeris with the exception of one transit (epoch 143) which is a partial transit. Therefore, there is currently no evidence for the existence of a companion to WASP-10b. Our results support the lack of TTVs of hot-Jupiters reported for the Kepler sample. 

Large- and small-scale structure of the intermediate- and high-velocity clouds towards the LMC and SMC.

We employ Ca II K and Na I D interstellar absorption-line spectroscopy of early-type stars in the Large and Small Magellanic Clouds (LMC, SMC) to investigate the large- and small-scale structure in foreground intermediate- and high-velocity clouds (I/HVCs). Data include FLAMES-GIRAFFE Ca II K observations of 403 stars in four open clusters, plus FEROS or UVES spectra of 156 stars in the LMC and SMC. The FLAMES observations are amongst the most extensive probes to date of Ca II structures on ∼20 arcsec scales in Magellanic I/HVCs. From the FLAMES data within a 0 ∘.∘.∘.5 field of view, the Ca II K equivalent width in the I/HVC components towards three clusters varies by factors of ≥10. There are no detections of molecular gas in absorption at intermediate or high velocities, although molecular absorption is present at LMC and Galactic velocities towards some sightlines. The FEROS/UVES data show Ca II K I/HVC absorption in ∼60 per cent of sightlines. The range in the Ca II/Na I ratio in I/HVCs is from –0.45 to +1.5 dex, similar to previous measurements for I/HVCs. In 10 sightlines we find Ca II/O I ratios in I/HVC gas ranging from 0.2 to 1.5 dex below the solar value, indicating either dust or ionization effects. In nine sightlines I/HVC gas is detected in both H I and Ca II at similar velocities, implying that the two elements form part of the same structure.