The binary fraction of planetary nebula central stars I. A high-precision, I-band excess search

In an attempt to determine how many planetary nebulae derive from binary interactions, we have started a project to measure their unbiased binary fraction. This number, when compared to the binary fraction of the presumed parent population can give a first handle on the origin of planetary nebulae. By detecting 27 bona fide central stars in the I band we have found that 30% of our sample have an I band excess between one and a few sigmas, possibly denoting companions brighter than M3-4V and with separations smaller than approximately 1000 AU. By accounting for the undetectable companions, we determine a de-biased binary fraction of 67-78% for all companions at all separations. We compare this number to a main sequence binary fraction of (50+/-4)% determined for spectral types F6V-G2V, appropriate if the progenitors of today's PN central star population is indeed the F6V-G2V stars. The error on our estimate could be between 10 and 30%. We conclude that the central star binary fraction may be larger than expected from the putative parent population. Using the more sensitive J band of a subset of 11 central stars, the binary fraction is 54% for companions brighter than approximately M5-6V and with separations smaller than about 900 AU. De-biassing this number we obtain a binary fraction of 100-107%. The two numbers should be the same and the discrepancy is likely due to small number statistics. We also present an accurately vetted compilation of observed main sequence star magnitudes, colours and masses, which can serve as a reference for future studies. We also present synthetic colours of hot stars as a function of temperature (20-170kK) and gravity (log g= 6-8) for Solar and PG1159 compositions.Comment: 22 pages, 6 figures, 12 tables, accepted by MNRA

Testing the binary hypothesis for the formation and shaping of planetary nebulae

There is no quantitative theory to explain why a high 80% of all planetary nebulae are non-spherical. The Binary Hypothesis states that a companion to the progenitor of a central star of planetary nebula is required to shape the nebula and even for a planetary nebula to be formed at all. A way to test this hypothesis is to estimate the binary fraction of central stars of planetary nebulae and to compare it with that of the main sequence population. Preliminary results from photometric variability and the infrared excess techniques indicate that the binary fraction of central stars of planetary nebulae is higher than that of the main sequence, implying that PNe could preferentially form via a binary channel. This article briefly reviews these results and current studies aiming to refine the binary fraction.Comment: EUROWD12 Proceeding

Planetary nebulae : getting closer to an unbiased binary fraction

Why 80% of planetary nebulae are not spherical is not yet understood. The Binary Hypothesis states that a companion to the progenitor of the central star of a planetary nebula is required to shape the nebula and even for a planetary nebula to be formed at all. A way to test this hypothesis is to estimate the binary fraction of central stars of planetary nebula and to compare it with the main sequence population. Preliminary results from photometric variability and infrared excess techniques indicate that the binary fraction of central stars of planetary nebulae is higher than that of the putative main sequence progenitor population, implying that PNe could be preferentially formed via a binary channel. This article briefly reviews these results and future studies aiming to refine the binary fraction.Comment: SF2A 2012 proceeding

Estimating the binary fraction of planetary nebulae central stars

During the past 20 years, the idea that non-spherical planetary nebulae (PN) may need a binary or planetary interaction to be shaped was discussed by various authors. It is now generally agreed that the varied morphologies of PN cannot be fully explained solely by single star evolution. Observationally, more binary central stars of planetary nebulae (CSPN) have been discovered, opening new possibilities to understand the connections between binarity and morphology. So far, \simeq 45 binary CSPN have been detected, most being close systems detected via flux variability. To determine the PN binary fraction, one needs a method to detect wider binaries. We present here recent results obtained with the various techniques described, concentrating on binary infrared excess observations aimed at detecting binaries of any separation.Comment: 2 pages, IAU 283: An Eye To The Future proceeding

Precision Pointing Control System (PPCS) system design and analysis

The precision pointing control system (PPCS) is an integrated system for precision attitude determination and orientation of gimbaled experiment platforms. The PPCS concept configures the system to perform orientation of up to six independent gimbaled experiment platforms to design goal accuracy of 0.001 degrees, and to operate in conjunction with a three-axis stabilized earth-oriented spacecraft in orbits ranging from low altitude (200-2500 n.m., sun synchronous) to 24 hour geosynchronous, with a design goal life of 3 to 5 years. The system comprises two complementary functions: (1) attitude determination where the attitude of a defined set of body-fixed reference axes is determined relative to a known set of reference axes fixed in inertial space; and (2) pointing control where gimbal orientation is controlled, open-loop (without use of payload error/feedback) with respect to a defined set of body-fixed reference axes to produce pointing to a desired target

SPITZER SAGE Observations of Large Magellanic Cloud Planetary Nebulae

We present IRAC and MIPS images and photometry of a sample of previously known planetary nebulae (PNe) from the SAGE survey of the Large Magellanic Cloud (LMC) performed with the Spitzer Space Telescope. Of the 233 known PNe in the survey field, 185 objects were detected in at least two of the IRAC bands, and 161 detected in the MIPS 24 micron images. Color-color and color-magnitude diagrams are presented using several combinations of IRAC, MIPS, and 2MASS magnitudes. The location of an individual PN in the color-color diagrams is seen to depend on the relative contributions of the spectral components which include molecular hydrogen, polycyclic aromatic hydrocarbons (PAHs), infrared forbidden line emission from the ionized gas, warm dust continuum, and emission directly from the central star. The sample of LMC PNe is compared to a number of Galactic PNe and found to not significantly differ in their position in color-color space. We also explore the potential value of IR PNe luminosity functions (LFs) in the LMC. IRAC LFs appear to follow the same functional form as the well-established [O III] LFs although there are several PNe with observed IR magnitudes brighter than the cut-offs in these LFs.Comment: 18 pages, 10 figures, 3 tables, to be published in the Astronomical Journal. Additional online data available at http://www.cfa.harvard.edu/irac/publications

Spectrum and Morphology of the Two Brightest Milagro Sources in the Cygnus Region: MGRO J2019+37 and MGRO J2031+41

The Cygnus region is a very bright and complex portion of the TeV sky, host to unidentified sources and a diffuse excess with respect to conventional cosmic-ray propagation models. Two of the brightest TeV sources, MGRO J2019+37 and MGRO J2031+41, are analyzed using Milagro data with a new technique, and their emission is tested under two different spectral assumptions: a power law and a power law with an exponential cutoff. The new analysis technique is based on an energy estimator that uses the fraction of photomultiplier tubes in the observatory that detect the extensive air shower. The photon spectrum is measured in the range 1 to 200 TeV using the last 3 years of Milagro data (2005-2008), with the detector in its final configuration. MGRO J2019+37 is detected with a significance of 12.3 standard deviations ($\sigma$), and is better fit by a power law with an exponential cutoff than by a simple power law, with a probability $>98$% (F-test). The best-fitting parameters for the power law with exponential cutoff model are a normalization at 10 TeV of $7^{+5}_{-2}\times10^{-10}$ $\mathrm{s^{-1}\: m^{-2}\: TeV^{-1}}$, a spectral index of $2.0^{+0.5}_{-1.0}$ and a cutoff energy of $29^{+50}_{-16}$ TeV. MGRO J2031+41 is detected with a significance of 7.3$\sigma$, with no evidence of a cutoff. The best-fitting parameters for a power law are a normalization of $2.4^{+0.6}_{-0.5}\times10^{-10}$ $\mathrm{s^{-1}\: m^{-2}\: TeV^{-1}}$ and a spectral index of $3.08^{+0.19}_{-0.17}$. The overall flux is subject to an $\sim$30% systematic uncertainty. The systematic uncertainty on the power law indices is $\sim$0.1. A comparison with previous results from TeV J2032+4130, MGRO J2031+41 and MGRO J2019+37 is also presented.Comment: 11 pages, 10 figure

Observation and Spectral Measurements of the Crab Nebula with Milagro

The Crab Nebula was detected with the Milagro experiment at a statistical significance of 17 standard deviations over the lifetime of the experiment. The experiment was sensitive to approximately 100 GeV - 100 TeV gamma ray air showers by observing the particle footprint reaching the ground. The fraction of detectors recording signals from photons at the ground is a suitable proxy for the energy of the primary particle and has been used to measure the photon energy spectrum of the Crab Nebula between ~1 and ~100 TeV. The TeV emission is believed to be caused by inverse-Compton up-scattering scattering of ambient photons by an energetic electron population. The location of a TeV steepening or cutoff in the energy spectrum reveals important details about the underlying electron population. We describe the experiment and the technique for distinguishing gamma-ray events from the much more-abundant hadronic events. We describe the calculation of the significance of the excess from the Crab and how the energy spectrum is fit. The fit is consistent with values measured by IACTs between 1 and 20 TeV. Fixing the spectral index to values that have been measured below 1 TeV by IACT experiments (2.4 to 2.6), the fit to the Milagro data suggests that Crab exhibits a spectral steepening or cutoff between about 20 to 40 TeV.Comment: Submitted to Astrophysical Journa

Design of a split Hopkinson pressure bar with partial lateral confinement

This paper presents the design of a modified split Hopkinson pressure bar (SHPB) where partial lateral con- finement of the specimen is provided by the inertia of a fluid annulus contained in a long steel reservoir. In contrast to unconfined testing, or a constant cell pressure applied before axial loading, lateral restraint is permitted to develop throughout the axial loading: this enables the high-strain-rate shear behaviour of soils to be characterised under conditions which are more representative of buried explosive events. A pressure transducer located in the wall of the reservoir allows lateral stresses to be quantified, and a dispersion-correction technique is used to provide accurate measurements of axial stress and strain. Preliminary numerical modelling is utilised to inform the experimental design, and the capability of the apparatus is demonstrated with specimen results for a dry quartz sand

Instability of LBV-stars against radial oscillations

In this study we consider the nonlinear radial oscillations exciting in LBV--stars with effective temperatures 1.5e4 K <= Teff <= 3e4 K, bolometric luminosities 1.2e6 L_odot <= L <= 1.9e6 L_odot and masses 35.7 M_odot <= M <= 49.1 M_odot. Hydrodynamic computations were carried out with initial conditions obtained from evolutionary sequences of population I stars (X=0.7, Z=0.02) with initial masses from 70M_odot to 90 M_odot. All hydrodynamical models show instability against radial oscillations with amplitude growth time comparable with dynamical time scale of the star. Radial oscillations exist in the form of nonlinear running waves propagating from the boundary of the compact core to the upper boundary of the hydrodynamical model. The velocity amplitude of outer layers is of several hundreds of km/s while the bolometric light amplitude does not exceed 0.2 mag. Stellar oscillations are not driven by the kappa-mechanism and are due to the instability of the gas with adiabatic exponent close to the critical value Gamma_1 = 4/3 due to the large contribution of radiation in the total pressure. The range of the light variation periods (6 day <= P <= 31 day) of hydrodynamical models agrees with periods of microvariability observed in LBV--stars.Comment: 14 pages, 5 figures, submitted to Astronomy Letter