How to process radio occultation data: 1. From time series of frequency residuals to vertical profiles of atmospheric and ionospheric properties

Expertise in processing radio occultation observations, which provide vertical profiles of atmospheric and ionospheric properties from measurements of the frequency of radio signals, is not widespread amongst the planetary science community. In order to increase the population of radio occultation processing experts, which will have positive consequences for this field, here we provide detailed instructions for one critical aspect of radio occultation data processing: how to obtain a series of bending angles as a function of the ray impact parameter from a time series of frequency residuals. As developed, this tool is valid only for one-way, single frequency occultations at spherically symmetric objects, and is thus not immediately applicable to either two-way occultations, such as those of Mars Express, or occultations at oblate objects, such as Jupiter or Saturn. This tool is demonstrated successfully on frequency residuals from a Mars Global Surveyor occultation at Mars, and the resultant set of bending angles and impact parameters are used to obtain vertical profiles of ionospheric electron density, neutral atmospheric number density, mass density, pressure, and temperature via the usual Abel transform. The root-mean-square difference between electron densities in the ionospheric profile derived herein and archived electron densities is 7×10[superscript 8] m[superscript −3]. At the lowest altitudes, temperatures in the neutral atmospheric profile derived herein differ from archived neutral temperatures by less than 0.1 K. Software programs that implement these procedures accompany this paper and may be used to extract scientifically useful data products from lower-level data sets

A Correlation Between Stellar Activity and Hot Jupiter Emission Spectra

We present evidence for a correlation between the observed properties of hot Jupiter emission spectra and the activity levels of the host stars measured using Ca II H & K emission lines. We find that planets with dayside emission spectra that are well-described by standard 1D atmosphere models with water in absorption (HD 189733, TrES-1, TrES-3, WASP-4) orbit chromospherically active stars, while planets with emission spectra that are consistent with the presence of a strong high-altitude temperature inversion and water in emission orbit quieter stars. We estimate that active G and K stars have Lyman alpha fluxes that are typically a factor of 4-7 times higher than quiet stars with analogous spectral types, and propose that the increased UV flux received by planets orbiting active stars destroys the compounds responsible for the formation of the observed temperature inversions. In this paper we also derive a model-independent method for differentiating between these two atmosphere types using the secondary eclipse depths measured in the 3.6 and 4.5 micron bands on the Spitzer Space Telescope, and argue that the observed correlation is independent of the inverted/non-inverted paradigm for classifying hot Jupiter atmospheres.Comment: 9 pages, 5 figures, accepted for publication in ApJ. The updated paper includes spectra for ten additional systems and a new section discussing the connection between chromospheric activity and UV flu

Mass and luminosity evolution of young stellar objects

A model of protostar mass and luminosity evolution in clusters gives new estimates of cluster age, protostar birthrate, accretion rate and mean accretion time. The model assumes constant protostar birthrate, core-clump accretion, and equally likely accretion stopping. Its parameters are set to reproduce the initial mass function, and to match protostar luminosity distributions in nearby star-forming regions. It obtains cluster ages and birthrates from the observed numbers of protostars and pre-main sequence (PMS) stars, and from the modal value of the protostar luminosity. In 31 embedded clusters and complexes the global cluster age is 1-3 Myr, matching available estimates based on optical spectroscopy and evolutionary tracks. This method of age estimation is simpler than optical spectroscopy, and is more useful for young embedded clusters where optical spectrocopy is not possible. In the youngest clusters, the protostar fraction decreases outward from the densest gas, indicating that the local star-forming age increases outward from a few 0.1 Myr in small protostar-dominated zones to a few Myr in large PMS-dominated zones.Comment: To appear in The Astrophysical Journal, Part

The Statistics of Albedo and Heat Recirculation on Hot Exoplanets

[Abridged] If both the day-side and night-side effective temperatures of a planet can be measured, it is possible to estimate its Bond albedo, 0<A_B<1, as well as its day-night heat redistribution efficiency, 0<epsilon<1. We attempt a statistical analysis of the albedo and redistribution efficiency for 24 transiting exoplanets that have at least one published secondary eclipse. For each planet, we show how to calculate a sub-stellar equilibrium temperature, T_0, and associated uncertainty. We then use a simple model-independent technique to estimate a planet's effective temperature from planet/star flux ratios. We use thermal secondary eclipse measurements -those obtained at lambda>0.8 micron- to estimate day-side effective temperatures, T_d, and thermal phase variations -when available- to estimate night-side effective temperature. We strongly rule out the "null hypothesis" of a single A_B and epsilon for all 24 planets. If we allow each planet to have different parameters, we find that low Bond albedos are favored (A_B<0.35 at 1 sigma confidence), which is an independent confirmation of the low albedos inferred from non-detection of reflected light. Our sample exhibits a wide variety of redistribution efficiencies. When normalized by T_0, the day-side effective temperatures of the 24 planets describe a uni-modal distribution. The dimensionless quantity T_d/T_0 exhibits no trend with the presence or absence of stratospheric inversions. There is also no clear trend between T_d/T_0 and T_0. That said, the 6 planets with the greatest sub-stellar equilibrium temperatures (T>2400 K) have low epsilon, as opposed to the 18 cooler planets, which show a variety of recirculation efficiencies. This hints that the very hottest transiting giant planets are qualitatively different from the merely hot Jupiters.Comment: 12 pages, 7 figures, ApJ accepted. Substantial Changes: more planets and data, refined analysis, different conclusion

Infrared Eclipses of the Strongly Irradiated Planet WASP-33b, and Oscillations of its Host Star

We observe two secondary eclipses of the strongly irradiated transiting planet WASP-33b in the Ks band, and one secondary eclipse each at 3.6- and 4.5 microns using Warm Spitzer. This planet orbits an A5V delta-Scuti star that is known to exhibit low amplitude non-radial p-mode oscillations at about 0.1-percent semi-amplitude. We detect stellar oscillations in all of our infrared eclipse data, and also in one night of observations at J-band out of eclipse. The oscillation amplitude, in all infrared bands except Ks, is about the same as in the optical. However, the stellar oscillations in Ks band have about twice the amplitude as seen in the optical, possibly because the Brackett-gamma line falls in this bandpass. We use our best-fit values for the eclipse depth, as well as the 0.9 micron eclipse observed by Smith et al., to explore possible states of the exoplanetary atmosphere, based on the method of Madhusudhan and Seager. On this basis we find two possible states for the atmospheric structure of WASP-33b. One possibility is a non-inverted temperature structure in spite of the strong irradiance, but this model requires an enhanced carbon abundance (C/O>1). The alternative model has solar composition, but an inverted temperature structure. Spectroscopy of the planet at secondary eclipse, using a spectral resolution that can resolve the water vapor band structure, should be able to break the degeneracy between these very different possible states of the exoplanetary atmosphere. However, both of those model atmospheres absorb nearly all of the stellar irradiance with minimal longitudinal re-distribution of energy, strengthening the hypothesis of Cowan et al. that the most strongly irradiated planets circulate energy poorly. Our measurement of the central phase of the eclipse yields e*cos(omega)=0.0003 +/-0.00013, which we regard as being consistent with a circular orbit.Comment: 23 pages, 9 figures, 3 tables, accepted for the Astrophysical Journa

The Balloon-Borne Large Aperture Submillimeter Telescope (BLAST) 2005: A 10 deg^2 Survey of Star Formation in Cygnus X

We present Cygnus X in a new multi-wavelength perspective based on an unbiased BLAST survey at 250, 350, and 500 micron, combined with rich datasets for this well-studied region. Our primary goal is to investigate the early stages of high mass star formation. We have detected 184 compact sources in various stages of evolution across all three BLAST bands. From their well-constrained spectral energy distributions, we obtain the physical properties mass, surface density, bolometric luminosity, and dust temperature. Some of the bright sources reaching 40 K contain well-known compact H II regions. We relate these to other sources at earlier stages of evolution via the energetics as deduced from their position in the luminosity-mass (L-M) diagram. The BLAST spectral coverage, near the peak of the spectral energy distribution of the dust, reveals fainter sources too cool (~ 10 K) to be seen by earlier shorter-wavelength surveys like IRAS. We detect thermal emission from infrared dark clouds and investigate the phenomenon of cold ``starless cores" more generally. Spitzer images of these cold sources often show stellar nurseries, but these potential sites for massive star formation are ``starless" in the sense that to date there is no massive protostar in a vigorous accretion phase. We discuss evolution in the context of the L-M diagram. Theory raises some interesting possibilities: some cold massive compact sources might never form a cluster containing massive stars; and clusters with massive stars might not have an identifiable compact cold massive precursor.Comment: 42 pages, 31 Figures, 6 table

Secondary Eclipse Photometry of WASP-4b with Warm Spitzer

We present photometry of the giant extrasolar planet WASP-4b at 3.6 and 4.5 micron taken with the Infrared Array Camera on board the Spitzer Space Telescope as part of Spitzer's extended warm mission. We find secondary eclipse depths of 0.319+/-0.031% and 0.343+/-0.027% for the 3.6 and 4.5 micron bands, respectively and show model emission spectra and pressure-temperature profiles for the planetary atmosphere. These eclipse depths are well fit by model emission spectra with water and other molecules in absorption, similar to those used for TrES-3 and HD 189733b. Depending on our choice of model, these results indicate that this planet has either a weak dayside temperature inversion or no inversion at all. The absence of a strong thermal inversion on this highly irradiated planet is contrary to the idea that highly irradiated planets are expected to have inversions, perhaps due the presence of an unknown absorber in the upper atmosphere. This result might be explained by the modestly enhanced activity level of WASP-4b's G7V host star, which could increase the amount of UV flux received by the planet, therefore reducing the abundance of the unknown stratospheric absorber in the planetary atmosphere as suggested in Knutson et al. (2010). We also find no evidence for an offset in the timing of the secondary eclipse and place a 2 sigma upper limit on |ecos(omega)| of 0.0024, which constrains the range of tidal heating models that could explain this planet's inflated radius.Comment: 8 pages, 7 figures (some in color), accepted for publication in Ap

STARSPOTS AND SPIN-ORBIT ALIGNMENT IN THE WASP-4 EXOPLANETARY SYSTEM

We present photometry of four transits of the exoplanet WASP-4b, each with a precision of approximately 500 ppm and a time sampling of 40-60 s. We have used the data to refine the estimates of the system parameters and ephemerides. During two of the transits we observed a short-lived, low-amplitude anomaly that we interpret as the occultation of a starspot by the planet. We also found evidence for a pair of similar anomalies in previously published photometry. The recurrence of these anomalies suggests that the stellar rotation axis is nearly aligned with the orbital axis, or else the starspot would not have remained on the transit chord. By analyzing the timings of the anomalies we find the sky-projected stellar obliquity to be lambda = –1[superscrip +14 over subscript –12] degrees. This result is consistent with (and more constraining than) a recent observation of the Rossiter-McLaughlin effect. It suggests that the planet migration mechanism preserved the initially low obliquity, or else that tidal evolution has realigned the system. Future applications of this method using data from the CoRoT and Kepler missions will allow spin-orbit alignment to be probed for many other exoplanets.United States. National Aeronautics and Space Administration (NASA Origins program award NNX09AD36G)United States. National Aeronautics and Space Administration (NASA Origins program award NNX09AB33G)Massachusetts Institute of Technology (Class of 1942)Caja de Ahorros y Pensiones de BarcelonaSpace Telescope Science Institute (Hubble Fellowship grant HF-51267.01-A

SPICES: Spectro-Polarimetric Imaging and Characterization of Exoplanetary Systems

SPICES (Spectro-Polarimetric Imaging and Characterization of Exoplanetary Systems) is a five-year M-class mission proposed to ESA Cosmic Vision. Its purpose is to image and characterize long-period extrasolar planets and circumstellar disks in the visible (450 - 900 nm) at a spectral resolution of about 40 using both spectroscopy and polarimetry. By 2020/22, present and near-term instruments will have found several tens of planets that SPICES will be able to observe and study in detail. Equipped with a 1.5 m telescope, SPICES can preferentially access exoplanets located at several AUs (0.5-10 AU) from nearby stars ($<$25 pc) with masses ranging from a few Jupiter masses to Super Earths ($\sim$2 Earth radii, $\sim$10 M$_{\oplus}$) as well as circumstellar disks as faint as a few times the zodiacal light in the Solar System

Overview of the massive young star-forming complex study in infrared and X-ray (MYStIX) project

The Massive Young Star-Forming Complex Study in Infrared and X-ray (MYStIX) seeks to characterize 20 OB-dominated young clusters and their environs at distances d ≤ 4 kpc using imaging detectors on the Chandra X-ray Observatory, Spitzer Space Telescope, and the United Kingdom InfraRed Telescope. The observational goals are to construct catalogs of star-forming complex stellar members with well-defined criteria and maps of nebular gas (particularly of hot X-ray-emitting plasma) and dust. A catalog of MYStIX Probable Complex Members with several hundred OB stars and 31,784 low-mass pre-main sequence stars is assembled. This sample and related data products will be used to seek new empirical constraints on theoretical models of cluster formation and dynamics, mass segregation, OB star formation, star formation triggering on the periphery of H II regions, and the survivability of protoplanetary disks in H II regions. This paper gives an introduction and overview of the project, covering the data analysis methodology and application to two star-forming regions: NGC 2264 and the Trifid Nebula