ALMA hints at the presence of two companions in the disk around HD 100546

HD 100546 is a well-studied Herbig Be star-disk system that likely hosts a close-in companion with compelling observational evidence for an embedded protoplanet at 68 AU. We present ALMA observations of the HD 100546 disk which resolve, for the first time, the gas and dust structure at (sub)mm wavelengths. The CO emission (at 345.795 GHz) originates from an extensive molecular disk (390 AU in radius) whereas the continuum emission is more compact (230 AU in radius) suggesting radial drift of the mm-sized grains. The CO emission is similar in extent to scattered light images indicating well-mixed gas and um-sized grains in the disk atmosphere. Assuming an azimuthally-symmetric disk, the continuum visibilities at long baselines (> 100 klambda) are reproduced by a compact ring with a width of 21 AU centered at 26 AU. An outer component is required to fit the short baselines: assuming a flat brightness distribution, the best-fit model is a ring with a width of 75 AU centered at 190 AU. The influence of a companion and protoplanet on the dust evolution is investigated. The companion at 10 AU facilitates the accumulation of mm-sized grains within a compact ring, ~20-30 AU, by ~10 Myr. The injection of a protoplanet at 1 Myr hastens the ring formation (~1.2 Myr) and also triggers the development of an outer ring (~100-200 AU). These observations provide additional evidence for the presence of a close-in companion and hint at dynamical clearing by a protoplanet at 68 AU.

Synthetic light curves and spectra for three-dimensional delayed-detonation models of type ia supernovae

In a companion paper, Seitenzahl et al. have presented a set of three-dimensional delayed detonation models for thermonuclear explosions of near-Chandrasekhar-mass white dwarfs (WDs). Here,we present multidimensional radiative transfer simulations that provide synthetic light curves and spectra for those models. The model sequence explores both changes in the strength of the deflagration phase (which is controlled by the ignition configuration in our models) and the WD central density. In agreement with previous studies, we find that the strength of the deflagration significantly affects the explosion and the observables. Variations in the central density also have an influence on both brightness and colour, but overall it is a secondary parameter in our set of models. In many respects, the models yield a good match to the observed properties of normal Type Ia supernovae (SNe Ia): peak brightness, rise/decline time-scales and synthetic spectra are all in reasonable agreement. There are, however, several differences. In particular, the models are systematically too red around maximum light, manifest spectral line velocities that are a little too high and yield I-band light curves that do not match observations. Although some of these discrepancies may simply relate to approximations made in the modelling, some pose real challenges to the models. If viewed as a complete sequence, our models do not reproduce the observed light-curve width- luminosity relation (WLR) of SNe Ia: all our models show rather similar B-band decline rates, irrespective of peak brightness. This suggests that simple variations in the strength of the deflagration phase in Chandrasekhar-mass deflagration-to-detonation models do not readily explain the observed diversity of normal SNe Ia. This may imply that some other parameter within the Chandrasekhar-mass paradigm is key to the WLR, or that a substantial fraction of normal SNe Ia arise from an alternative explosion scenario.

Search for the Return of Activity in Active Asteroid 176P/LINEAR

We present the results of a search for the reactivation of active asteroid 176P/LINEAR during its 2011 perihelion passage using deep optical observations obtained before, during, and after that perihelion passage. Deep composite images of 176P constructed from data obtained between 2011 June and 2011 December show no visible signs of activity, while photometric measurements of the object during this period also show no significant brightness enhancements similar to that observed for 176P between 2005 November and 2005 December when it was previously observed to be active. An azimuthal search for dust emission likewise reveals no evidence for directed emission (i.e., a tail, as was previously observed for 176P), while a one-dimensional surface brightness profile analysis shows no indication of a spherically symmetric coma at any time in 2011. We conclude that 176P did not in fact exhibit activity in 2011, at least not on the level on which it exhibited activity in 2005, and suggest that this could be due to the devolatization or mantling of the active site responsible for its activity in 2005.

Main-belt Comet P/2012 T1 (PANSTARRS)

We present initial results from observations and numerical analyses aimed at characterizing the main-belt comet P/2012 T1 (PANSTARRS). Optical monitoring observations were made between 2012 October and 2013 February using the University of Hawaii 2.2 m telescope, the Keck I telescope, the Baade and Clay Magellan telescopes, Faulkes Telescope South, the Perkins Telescope at Lowell Observatory, and the Southern Astrophysical Research Telescope. The object's intrinsic brightness approximately doubles from the time of its discovery in early October until mid-November and then decreases by ~60% between late December and early February, similar to photometric behavior exhibited by several other main-belt comets and unlike that exhibited by disrupted asteroid (596) Scheila. We also used Keck to conduct spectroscopic searches for CN emission as well as absorption at 0.7 μm that could indicate the presence of hydrated minerals, finding an upper limit CN production rate of Q CN <1.5 × 1023 mol s-1, from which we infer a water production rate of Q_H_2O100 Myr and is unlikely to be a recently implanted interloper from the outer solar system, while a search for potential asteroid family associations reveals that it is dynamically linked to the ~155 Myr old Lixiaohua asteroid family. Some of the data presented herein were obtained at the W. M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California, and the National Aeronautics and Space Administration, and made possible by the generous financial support of the W. M. Keck Foundation, the Magellan Telescopes located at Las Campanas Observatory, Chile, and the Southern Astrophysical Research (SOAR) telescope, which is a joint project of the Ministério da Ciência, Tecnologia, e Inovação (MCTI) da República Federativa do Brasil, the U.S. National Optical Astronomy Observatory (NOAO), the University of North Carolina at Chapel Hill (UNC), and Michigan State University (MSU).

Observational Constraints on the Catastrophic Disruption Rate of Small Main Belt Asteroids

We have calculated 90% confidence limits on the steady-state rate of catastrophic disruptions of main belt asteroids in terms of the absolute magnitude at which one catastrophic disruption occurs per year  as a function of the post-disruption increase in brightness (Δm) and subsequent brightness decay rate (τ  ). The confidence limits were calculated using the brightest unknown main belt asteroid (V=18.5) detected with the Pan-STARRS1 (Pan-STARRS1) telescope. We measured the Pan-STARRS1’s catastrophic disruption detection efficiency over a 453-day interval using the Pan-STARRS moving object processing system (MOPS) and a simple model for the catastrophic disruption event’s photometric behavior in a small aperture centered on the catastrophic disruption event. We then calculated the  contours in the ranges from  and  encompassing measured values from known cratering and disruption events and our model’s predictions. Our simplistic catastrophic disruption model suggests that  and  which would imply that H_0≳28—strongly inconsistent withH_{0,B2005=23.26±0.02} predicted by Bottke et al. (Bottke, W.F., Durda, D.D., Nesvorný, D., Jedicke, R., Morbidelli, A., Vokrouhlický, D., Levison, H.F. [2005]. Icarus, 179, 63–94.) using purely collisional models. However, if we assume that H_0=H0,B2005 our results constrain , inconsistent with our simplistic impact-generated catastrophic disruption model. We postulate that the solution to the discrepancy is that >99% of main belt catastrophic disruptions in the size range to which this study was sensitive (∼100 m) are not impact-generated, but are instead due to fainter rotational breakups, of which the recent discoveries of disrupted asteroids P/2013 P5 and P/2013 R3 are probable examples. We estimate that current and upcoming asteroid surveys may discover up to 10 catastrophic disruptions/year brighter than V=18.5.

Sublimation-Driven Activity in Main-Belt Comet 313p/Gibbs

We present an observational and dynamical study of newly discovered main-belt comet 313P/Gibbs. We find that the object is clearly active both in observations obtained in 2014 and in precovery observations obtained in 2003 by the Sloan Digital Sky Survey, strongly suggestingthat its activity is sublimation-driven. This conclusion is supported by a photometric analysis showing an increase in the total brightness of the comet over the 2014 observing period, and dust modeling resultsshowing that the dust emission persists over at least three months during both active periods, where we find start dates for emission nolater than 2003 July 24 ± 10 for the 2003 active period and 2014 July 28 ± 10 for the 2014 active period. From serendipitous observations by the Subaru Telescope in 2004 when the object was apparently inactive, we estimate that the nucleus has an absolute R-band magnitude of HR = 17.1 ± 0.3, corresponding to aneffective nucleus radius of re ∼ 1.00 ± 0.15 km.The object’s faintness at that time means we cannot rule out the presence of activity, and so this computed radius should be consideredan upper limit. We find that 313P’s orbit is intrinsically chaotic, having a Lyapunov time of Tl = 12,000 yr and beinglocated near two three-body mean-motion resonances with Jupiter andSaturn, 11J-1S-5A and 10J+12S-7A, yet appears stable over >50 Myr in an apparent example of stable chaos. We furthermore find that 313P is the second main-belt comet, after P/2012 T1 (PANSTARRS), to belong tothe ∼155 Myr old Lixiaohua asteroid family.

Absolute Magnitudes of Pan-STARRS PS1 Asteroids

Absolute magnitude (H) of an asteroid is a fundamental parameter describing the size and the apparent brightness of the body. Because of its surface shape, properties and changing illumination, the brightness changes with the geometry and is described by the phase function governed by the slope parameter (G). Although many years have been spent on detailed observations of individual asteroids to provide H and G, vast majority of minor planets have H based on assumed G and due to the input photometry from multiple sources the errors of these values are unknown. We compute H of ~ 180 000 and G of few thousands asteroids observed with the Pan-STARRS PS1 telescope in well defined photometric systems. The mean photometric error is 0.04 mag. Because on average there are only 7 detections per asteroid in our sample, we employed a Monte Carlo (MC) technique to generate clones simulating all possible rotation periods, amplitudes and colors of detected asteroids. Known asteroid colors were taken from the SDSS database. We used debiased spin and amplitude distributions dependent on size, spectral class distributions of asteroids dependent on semi-major axis and starting values of G from previous works. H and G (G12 respectively) were derived by phase functions by Bowell et al. (1989) and Muinonen et al. (2010). We confirmed that there is a positive systematic offset between H based on PS1 asteroids and Minor Planet Center database up to -0.3 mag peaking at 14. Similar offset was first mentioned in the analysis of SDSS asteroids and was believed to be solved by weighting and normalizing magnitudes by observatory codes. MC shows that there is only a negligible difference between Bowell's and Muinonen's solution of H. However, Muinonen's phase function provides smaller errors on H. We also derived G and G12 for thousands of asteroids. For known spectral classes, slope parameters agree with the previous work in general, however, the standard deviation of G in our sample is twice as larger, most likely due to sparse phase curve sampling. In the near future we plan to complete the H and G determination for all PS1 asteroids (500,000) and publish H and G values online. This work was supported by NASA grant No. NNX12AR65G.

Solar phase functions of cometary nuclei

Relatively few measurements of the solar phase function of cometary nuclei exist, despite the importance of this parameter in determining accurate sizes and its use in modeling surface properties. We make use of robotic telescopes and servicemode observing to monitor cometary nuclei over months at a time, combining intensive observations at a single epoch with regular short light-curve segments to efficiently account for brightness changes due to both nucleus rotation and changing solar phase angle. We present our latest results on comets 8P/Tuttle, 14P/Wolf, 67P/Churyumov- Gerasimenko and 110P/Hartley 3.

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. 

SN 2009ib: a Type II-P supernova with an unusually long plateau

We present optical and near-infrared photometry and spectroscopy of SN 2009ib, a Type II-P supernova in NGC 1559. This object has moderate brightness, similar to those of the intermediate-luminosity SNe 2008in and 2009N. Its plateau phase is unusually long, lasting for about 130 d after explosion. The spectra are similar to those of the subluminous SN 2002gd, with moderate expansion velocities. We estimate the Ni-56 mass produced as 0.046 +/- A 0.015 M-aS (TM). We determine the distance to SN 2009ib using both the expanding photosphere method (EPM) and the standard candle method. We also apply EPM to SN 1986L, a Type II-P SN that exploded in the same galaxy. Combining the results of different methods, we conclude the distance to NGC 1559 as D = 19.8 +/- A 3.0 Mpc. We examine archival, pre-explosion images of the field taken with the Hubble Space Telescope, and find a faint source at the position of the SN, which has a yellow colour [(V - I)(0) = 0.85 mag]. Assuming it is a single star, we estimate its initial mass as M-ZAMS = 20 M-aS (TM). We also examine the possibility, that instead of the yellow source the progenitor of SN 2009ib is a red supergiant star too faint to be detected. In this case, we estimate the upper limit for the initial zero-age main sequence (ZAMS) mass of the progenitor to be similar to 14-17 M-aS (TM). In addition, we infer the physical properties of the progenitor at the explosion via hydrodynamical modelling of the observables, and estimate the total energy as similar to 0.55 x 10(51) erg, the pre-explosion radius as similar to 400 R-aS (TM), and the ejected envelope mass as similar to 15 M-aS (TM), which implies that the mass of the progenitor before explosion was similar to 16.5-17 M-aS (TM).

Zooming In on the Progenitors of Superluminous Supernovae With the HST

We present Hubble Space Telescope (HST) rest-frame ultraviolet imaging of the host galaxies of 16 hydrogen-poor superluminous supernovae (SLSNe), including 11 events from the Pan-STARRS Medium Deep Survey. Taking advantage of the superb angular resolution of HST, we characterize the galaxies' morphological properties, sizes, and star formation rate (SFR) densities. We determine the supernova (SN) locations within the host galaxies through precise astrometric matching and measure physical and host-normalized offsets as well as the SN positions within the cumulative distribution of UV light pixel brightness. We find that the host galaxies of H-poor SLSNe are irregular, compact dwarf galaxies, with a median half-light radius of just 0.9 kpc. The UV-derived SFR densities are high ([Sigma(SFR)] similar or equal to 0.1M(circle dot) yr(-1) kpc(-1)), suggesting that SLSNe form in overdense environments. Their locations trace the UV light of their host galaxies, with a distribution intermediate between that of long-duration gamma-ray bursts (LGRBs; which are strongly clustered on the brightest regions of their hosts) and a uniform distribution (characteristic of normal core-collapse SNe), though cannot be statistically distinguished from either with the current sample size. Taken together, this strengthens the picture that SLSN progenitors require different conditions than those of ordinary core-collapse SNe to form and that they explode in broadly similar galaxies as do LGRBs. If the tendency for SLSNe to be less clustered on the brightest regions than are LGRBs is confirmed by a larger sample, this would indicate a different, potentially lower-mass progenitor for SLSNe than LRGBs.

Ultraluminous Supernovae as a New Probe of the Interstellar Medium in Distant Galaxies

We present the Pan-STARRS1 discovery and light curves, and follow-up MMT and Gemini spectroscopy of an ultraluminous supernova (ULSN; dubbed PS1-11bam) at a redshift of z = 1.566 with a peak brightness of M UV ≈ -22.3 mag. PS1-11bam is one of the highest redshift spectroscopically confirmed SNe known to date. The spectrum exhibits broad absorption features typical of previous ULSNe (e.g., C II, Si III), and strong and narrow Mg II and Fe II absorption lines from the interstellar medium (ISM) of the host galaxy, confirmed by an [O II]λ3727 emission line at the same redshift. The equivalent widths of the Fe II λ2600 and Mg II λ2803 lines are in the top quartile of the quasar intervening absorption system distribution, but are weaker than those of gamma-ray burst intrinsic absorbers (i.e., GRB host galaxies). We also detect the host galaxy in pre-explosion Pan-STARRS1 data and find that its UV spectral energy distribution is best fit with a young stellar population age of τ* ≈ 15-45 Myr and a stellar mass of M * ≈ (1.1-2.6) × 109 M ⊙ (for Z = 0.05-1 Z ⊙). The star formation rate inferred from the UV continuum and [O II]λ3727 emission line is ≈10 M ⊙ yr-1, higher than in previous ULSN hosts. PS1-11bam provides the first direct demonstration that ULSNe can serve as probes of the ISM in distant galaxies. The depth and red sensitivity of PS1 are uniquely suited to finding such events at cosmologically interesting redshifts (z ~ 1-2); the future combination of LSST and 30 m class telescopes promises to extend this technique to z ~ 4.

SN2014J gamma rays from the 56Ni decay chain

Context. The detection and measurement of gamma-ray lines from the decaychain of 56Ni provides unique information about the explosionin supernovae. SN2014J at 3.3 Mpc is a sufficiently-nearby supernova oftype Ia so that such measurements have been feasible with the gamma-rayspectrometer SPI on ESA's INTEGRAL gamma-ray observatory.

Aims:The 56Ni freshly produced in the supernova is understood topower the optical light curve, because it emits gamma rays upon itsradioactive decay first to 56Co and then to 56Fe.Gamma-ray lines from 56Co decay are expected to becomedirectly visible through the white dwarf material several weeks afterthe explosion, as they progressively penetrate the overlying material ofthe supernova envelope, which is diluted as it expands. The lines areexpected to be Doppler-shifted or broadened from the kinematics of the56Ni ejecta. We aim to exploit high-resolution gamma-rayspectroscopy with the SPI spectrometer on INTEGRAL toward constrainingthe 56Ni distribution and kinematics in this supernova.

Methods: We use the observations with the SPI spectrometer onINTEGRAL, together with an improved instrumental background method.

Results: We detect the two main lines from 56Co decay at847 and 1238 keV, which are significantly Doppler-broadened, and atintensities (3.65 ± 1.21) × 10-4 and (2.27± 0.69) × 10-4 ph cm-2s-1, respectively, at their brightness maximum. We measuretheir rise toward a maximum after about 60-100 days and a declinethereafter. The intensity ratio of the two lines is found to beconsistent with expectations from 56Co decay (0.62 ±0.28 at brightness maximum, the expected ratio is 0.68). We find thatthe broad lines seen in the late, gamma-ray transparent phase are notrepresentative of the early gamma-ray emission, and notice instead thatthe emission spectrum is complex and irregular until the supernova isfully transparent to gamma rays, with progressive uncovering of the bulkof 56Ni. We infer that the explosion morphology is notspherically symmetric, both in the distribution of 56Ni andin the unburnt material which occults the 56Co emission.After we compare light curves from different plausible models, theresulting 56Ni mass is determined to be 0.49 ± 0.09M⊙.

A comparative study of Type II-P and II-L supernova rise times as exemplified by the case of LSQ13cuw

We report on our findings based on the analysis of observations of the Type II-L supernova LSQ13cuw within the framework of currently accepted physical predictions of core-collapse supernova explosions. LSQ13cuw was discovered within a day of explosion, hitherto unprecedented for Type II-L supernovae. This motivated a comparative study of Type II-P and II-L supernovae with relatively well-constrained explosion epochs and rise times to maximum (optical) light. From our sample of twenty such events, we find evidence of a positive correlation between the duration of the rise and the peak brightness. On average, SNe II-L tend to have brighter peak magnitudes and longer rise times than SNe II-P. However, this difference is clearest only at the extreme ends of the rise time versus peak brightness relation. Using two different analytical models, we performed a parameter study to investigate the physical parameters that control the rise time behaviour. In general, the models qualitatively reproduce aspects of the observed trends. We find that the brightness of the optical peak increases for larger progenitor radii and explosion energies, and decreases for larger masses. The dependence of the rise time on mass and explosion energy is smaller than the dependence on the progenitor radius. We find no evidence that the progenitors of SNe II-L have significantly smaller radii than those of SNe II-P.

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).