Astrophysical observations indicate that there is roughly five times more dark matter in the Universe than ordinary baryonic matter, with an even larger amount of the Universe's energy content due to dark energy. So far, the microscopic properties of these dark components have remained shrouded in mystery. In addition, even the five percent of ordinary matter in our Universe has yet to be understood, since the Standard Model of particle physics lacks any consistent explanation for the predominance of matter over antimatter. Inspired by these central problems of modern physics, we present here a direct search for interactions of antimatter with dark matter, and place direct constraints on the interaction of ultra-light axion-like particles $-$ one of the dark-matter candidates $-$ and antiprotons. If antiprotons exhibit a stronger coupling to these dark-matter particles than protons, such a CPT-odd coupling could provide a link between dark matter and the baryon asymmetry in the Universe. We analyse spin-flip resonance data acquired with a single antiproton in a Penning trap [Smorra et al., Nature 550, 371 (2017)] in the frequency domain to search for spin-precession effects from ultra-light axions with a characteristic frequency governed by the mass of the underlying particle. Our analysis constrains the axion-antiproton interaction parameter $f_a/C_{\overline{p}}$ to values greater than $0.1$ to $0.6$ GeV in the mass range from $2 \times 10^{-23}$ to $4 \times 10^{-17}\,$eV/$c^2$, improving over astrophysical antiproton bounds by up to five orders of magnitude. In addition, we derive limits on six combinations of previously unconstrained Lorentz-violating and CPT-violating terms of the non-minimal Standard Model Extension.

Atmospheric retrieval determines the properties of an atmosphere based on its measured spectrum. The low signal-to-noise ratio of exoplanet observations require a Bayesian approach to determine posterior probability distributions of each model parameter, given observed spectra. This inference is computationally expensive, as it requires many executions of a costly radiative transfer (RT) simulation for each set of sampled model parameters. Machine learning (ML) has recently been shown to provide a significant reduction in runtime for retrievals, mainly by training inverse ML models that predict parameter distributions, given observed spectra, albeit with reduced posterior accuracy. Here we present a novel approach to retrieval by training a forward ML surrogate model that predicts spectra given model parameters, providing a fast approximate RT simulation that can be used in a conventional Bayesian retrieval framework without significant loss of accuracy. We demonstrate our method on the emission spectrum of HD 189733 b and find Bhattacharyya coefficients of 97.74 -- 99.74% between our 1D marginalized posterior distributions and those of the Bayesian Atmospheric Radiative Transfer (BART) code. Our retrieval method is ~20x faster than BART when run on an Intel i7-4770 central processing unit (CPU). Neural-network computation using an NVIDIA Titan Xp graphics processing unit is ~600x faster than BART on that CPU.

Fundamental coding and software development skills are increasingly necessary for success in nearly every aspect of astronomical and astrophysical research as large surveys and high resolution simulations become the norm. However, professional training in these skills is inaccessible or impractical for many members of our community. Students and professionals alike have been expected to acquire these skills on their own, apart from formal classroom curriculum or on-the-job training. Despite the recognized importance of these skills, there is little opportunity to develop them - even for interested researchers. To ensure a workforce capable of taking advantage of the computational resources and the large volumes of data coming in the next decade, we must identify and support ways to make software development training widely accessible to community members, regardless of affiliation or career level. To develop and sustain a technology capable astronomical and astrophysical workforce, we recommend that agencies make funding and other resources available in order to encourage, support and, in some cases, require progress on necessary training, infrastructure and policies. In this white paper, we focus on recommendations for how funding agencies can lead in the promotion of activities to support the astronomy and astrophysical workforce in the 2020s.

As astronomy moves into the era of large-scale time-domain surveys, we are seeing a flood of new transient and variable sources which will reach biblical proportions with the advent of LSST. A key strategic challenge for astronomy in this era is the lack of suitable spectroscopic followup facilities. In response to this need, we have developed the PolyOculus approach for producing large-area-equivalent telescopes by using fiber optics to link modules of multiple semi-autonomous, small, inexpensive, commercial-off-the-shelf telescopes. Crucially, this scalable design has construction costs which are $>10x$ lower than equivalent traditional large-area telescopes. In addition, PolyOculus is inherently highly automated and well-suited for remote operations. Development of this technology will enable the expansion of major research efforts in the LSST era to a host of smaller universities and colleges, including primarily-undergraduate institutions, for budgets consistent with their educational expenditures on similar facilities. We propose to develop and deploy a 1.6-m prototype demonstrator at the Mt. Laguna Observatory in California, followed by a full-scale 5-meter-class PolyOculus facility for linkage to existing and upcoming time-domain surveys.

We propose an experiment, the Cosmic Accelerometer, designed to yield velocity precision of $\leq 1$ cm/s with measurement stability over years to decades. The first-phase Cosmic Accelerometer, which is at the scale of the Astro2020 Small programs, will be ideal for precision radial velocity measurements of terrestrial exoplanets in the Habitable Zone of Sun-like stars. At the same time, this experiment will serve as the technical pathfinder and facility core for a second-phase larger facility at the Medium scale, which can provide a significant detection of cosmological redshift drift on a 6-year timescale. This larger facility will naturally provide further detection/study of Earth twin planet systems as part of its external calibration process. This experiment is fundamentally enabled by a novel low-cost telescope technology called PolyOculus, which harnesses recent advances in commercial off the shelf equipment (telescopes, CCD cameras, and control computers) combined with a novel optical architecture to produce telescope collecting areas equivalent to standard telescopes with large mirror diameters. Combining a PolyOculus array with an actively-stabilized high-precision radial velocity spectrograph provides a unique facility with novel calibration features to achieve the performance requirements for the Cosmic Accelerometer.

Software is a critical part of modern research, and yet there are insufficient mechanisms in the scholarly ecosystem to acknowledge, cite, and measure the impact of research software. The majority of academic fields rely on a one-dimensional credit model whereby academic articles (and their associated citations) are the dominant factor in the success of a researcher's career. In the petabyte era of astronomical science, citing software and measuring its impact enables academia to retain and reward researchers that make significant software contributions. These highly skilled researchers must be retained to maximize the scientific return from petabyte-scale datasets. Evolving beyond the one-dimensional credit model requires overcoming several key challenges, including the current scholarly ecosystem and scientific culture issues. This white paper will present these challenges and suggest practical solutions for elevating the role of software as a product of the research enterprise.

Jupiter's banded structure undergoes strong temporal variations, changing the visible and infrared appearance of the belts and zones in a complex and turbulent way due to physical processes that are not yet understood. In this study we use ground-based 5-$\mu$m infrared data captured between 1984 and 2018 by 8 different instruments mounted on the Infrared Telescope Facility in Hawai'i and on the Very Large Telescope in Chile to analyze and characterize the long-term variability of Jupiter's cloud-forming region at the 1-4 bar pressure level. The data show a large temporal variability mainly at the equatorial and tropical latitudes, with a smaller temporal variability at mid-latitudes. We also compare the 5-$\mu$m-bright and -dark regions with the locations of the visible zones and belts and we find that these regions are not always co-located, specially in the southern hemisphere. We also present Lomb-Scargle and Wavelet Transform analyzes in order to look for possible periodicities of the brightness changes that could help us understand their origin and predict future events. We see that some of these variations occur periodically in time intervals of 4-8 years. The reasons of these time intervals are not understood and we explore potential connections to both convective processes in the deeper weather layer and dynamical processes in the upper troposphere and stratosphere. Finally we perform a Principal Component analysis to reveal a clear anticorrelation on the 5-$\mu$m brightness changes between the North Equatorial Belt and the South Equatorial Belt, suggesting a possible connection between the changes in these belts.

We report Spitzer Space Telescope observations during predicted transits of the exoplanet Proxima Centauri b. As the nearest terrestrial habitable-zone planet we will ever discover, any potential transit of Proxima b would place strong constraints on its radius, bulk density, and atmosphere. Subsequent transmission spectroscopy and secondary-eclipse measurements could then probe the atmospheric chemistry, physical processes, and orbit, including a search for biosignatures. However, our photometric results rule out planetary transits at the 200~ppm level at 4.5$~{\mu}m$, yielding a 3$\sigma$ upper radius limit of 0.4~$R_\rm{\oplus}$ (Earth radii). Previous claims of possible transits from optical ground- and space-based photometry were likely correlated noise in the data from Proxima Centauri's frequent flaring. Follow-up observations should focus on planetary radio emission, phase curves, and direct imaging. Our study indicates dramatically reduced stellar activity at near-to-mid infrared wavelengths, compared to the optical. Proxima b is an ideal target for space-based infrared telescopes, if their instruments can be configured to handle Proxima's brightness.

The James Webb Space Telescope (JWST) presents the opportunity to transform our understanding of planets and the origins of life by revealing the atmospheric compositions, structures, and dynamics of transiting exoplanets in unprecedented detail. However, the high-precision, time-series observations required for such investigations have unique technical challenges, and prior experience with other facilities indicates that there will be a steep learning curve when JWST becomes operational. In this paper we describe the science objectives and detailed plans of the Transiting Exoplanet Community Early Release Science (ERS) Program, which is a recently approved program for JWST observations early in Cycle 1. The goal of this project, for which the obtained data will have no exclusive access period, is to accelerate the acquisition and diffusion of technical expertise for transiting exoplanet observations with JWST, while also providing a compelling set of representative datasets that will enable immediate scientific breakthroughs. The Transiting Exoplanet Community ERS Program will exercise the time-series modes of all four JWST instruments that have been identified as the consensus highest priorities, observe the full suite of transiting planet characterization geometries (transits, eclipses, and phase curves), and target planets with host stars that span an illustrative range of brightnesses. The observations in this program were defined through an inclusive and transparent process that had participation from JWST instrument experts and international leaders in transiting exoplanet studies. Community engagement in the project will be centered on a two-phase Data Challenge that culminates with the delivery of planetary spectra, time-series instrument performance reports, and open-source data analysis toolkits in time to inform the agenda for Cycle 2 of the JWST mission.

The globular cluster (GC) system of our Galaxy contains four planetary nebulae (PNe): K 648 (or Ps 1) in M15, IRAS 18333-2357 in M22, JaFu 1 in Pal 6, and JaFu 2 in NGC 6441. Because single-star evolution at the low stellar mass of present-epoch GCs was considered incapable of producing visible PNe, their origin presented a puzzle. We imaged the PN JaFu 1 with the Hubble Space Telescope (HST) to obtain photometry of its central star (CS) and high-resolution morphological information. We imaged IRAS 18333-2357 with better depth and resolution, and we analyzed its archival HST spectra to constrain its CS temperature and luminosity. All PNe in Galactic GCs now have high-quality HST data, allowing us to improve CS mass estimates. We find reasonably consistent masses between 0.53 and 0.58 Msun for all four objects, though estimates vary when adopting different stellar evolutionary calculations. The CS mass of IRAS 18333-2357, though, depends strongly on its temperature, which remains elusive due to reddening uncertainties. For all four objects, we consider their CS and nebular masses, their morphologies, and other incongruities to assess the likelihood that these objects formed from binary stars. Although generally limited by uncertainties (~0.02 Msun) in post-AGB tracks and core mass vs. luminosity relations, the high-mass CS in K 648 indicates a binary origin. The CS of JaFu 1 exhibits compact bright [O III] and Halpha emission, like EGB 6, suggesting a binary companion or disk. Evidence is weaker for a binary origin of JaFu 2.

We present {\em Spitzer} secondary-eclipse observations of the hot Jupiter HAT-P-13 b in the 3.6 {\micron} and 4.5 {\micron} bands. HAT-P-13 b inhabits a two-planet system with a configuration that enables constraints on the planet's second Love number, \math{k\sb{2}}, from precise eccentricity measurements, which in turn constrains models of the planet's interior structure. We exploit the direct measurements of \math{e \cos \omega} from our secondary-eclipse data and combine them with previously published radial velocity data to generate a refined model of the planet's orbit and thus an improved estimate on the possible interval for \math{k\sb{2}}. We report eclipse phases of \math{0.49154 \pm 0.00080} and \math{0.49711 \pm 0.00083} and corresponding \math{e \cos \omega} estimates of \math{-0.0136 \pm 0.0013} and \math{-0.0048 \pm 0.0013}. Under the assumptions of previous work, our estimate of \math{k\sb{2}} of 0.81 {\pm} 0.10 is consistent with the lower extremes of possible core masses found by previous models, including models with no solid core. This anomalous result challenges both interior models and the dynamical assumptions that enable them, including the essential assumption of apsidal alignment. We also report eclipse depths of 0.081\% {\pm} 0.008\% in the 3.6 {\micron} channel and 0.088 \% {\pm} 0.028 \% in the 4.5 {\micron} channel. These photometric results are non-uniquely consistent with solar-abundance composition without any thermal inversion.

Time-correlated noise is a significant source of uncertainty when modeling exoplanet light-curve data. A correct assessment of correlated noise is fundamental to determine the true statistical significance of our findings. Here we review three of the most widely used correlated-noise estimators in the exoplanet field, the time-averaging, residual-permutation, and wavelet-likelihood methods. We argue that the residual-permutation method is unsound in estimating the uncertainty of parameter estimates. We thus recommend to refrain from this method altogether. We characterize the behavior of the time averaging's rms-vs.-bin-size curves at bin sizes similar to the total observation duration, which may lead to underestimated uncertainties. For the wavelet-likelihood method, we note errors in the published equations and provide a list of corrections. We further assess the performance of these techniques by injecting and retrieving eclipse signals into synthetic and real Spitzer light curves, analyzing the results in terms of the relative-accuracy and coverage-fraction statistics. Both the time-averaging and wavelet-likelihood methods significantly improve the estimate of the eclipse depth over a white-noise analysis (a Markov-chain Monte Carlo exploration assuming uncorrelated noise). However, the corrections are not perfect, when retrieving the eclipse depth from Spitzer datasets, these methods covered the true (injected) depth within the 68\% credible region in only $\sim$45--65\% of the trials. Lastly, we present our open-source model-fitting tool, Multi-Core Markov-Chain Monte Carlo ({MC$^3$}). This package uses Bayesian statistics to estimate the best-fitting values and the credible regions for the parameters for a (user-provided) model. {MC$^3$} is a Python/C code, available at https://github.com/pcubillos/MCcubed.

The James Webb Space Telescope will revolutionize transiting exoplanet atmospheric science due to its capability for continuous, long-duration observations and its larger collecting area, spectral coverage, and spectral resolution compared to existing space-based facilities. However, it is unclear precisely how well JWST will perform and which of its myriad instruments and observing modes will be best suited for transiting exoplanet studies. In this article, we describe a prefatory JWST Early Release Science (ERS) program that focuses on testing specific observing modes to quickly give the community the data and experience it needs to plan more efficient and successful future transiting exoplanet characterization programs. We propose a multi-pronged approach wherein one aspect of the program focuses on observing transits of a single target with all of the recommended observing modes to identify and understand potential systematics, compare transmission spectra at overlapping and neighboring wavelength regions, confirm throughputs, and determine overall performances. In our search for transiting exoplanets that are well suited to achieving these goals, we identify 12 objects (dubbed "community targets") that meet our defined criteria. Currently, the most favorable target is WASP-62b because of its large predicted signal size, relatively bright host star, and location in JWST's continuous viewing zone. Since most of the community targets do not have well-characterized atmospheres, we recommend initiating preparatory observing programs to determine the presence of obscuring clouds/hazes within their atmospheres. Measurable spectroscopic features are needed to establish the optimal resolution and wavelength regions for exoplanet characterization. Other initiatives from our proposed ERS program include testing the instrument brightness limits and performing phase-curve observations.(Abridged)

We present the results of a meeting on numerical simulations of ionized nebulae held at the University of Kentucky in conjunction with the celebration of the 70th birthdays of Profs. Donald Osterbrock and Michael Seaton.

We present an open-source Thermochemical Equilibrium Abundances (TEA) code that calculates the abundances of gaseous molecular species. The code is based on the methodology of White et al. (1958) and Eriksson (1971). It applies Gibbs free-energy minimization using an iterative, Lagrangian optimization scheme. Given elemental abundances, TEA calculates molecular abundances for a particular temperature and pressure or a list of temperature-pressure pairs. We tested the code against the method of Burrows & Sharp (1999), the free thermochemical equilibrium code CEA (Chemical Equilibrium with Applications), and the example given by White et al. (1958). Using their thermodynamic data, TEA reproduces their final abundances, but with higher precision. We also applied the TEA abundance calculations to models of several hot-Jupiter exoplanets, producing expected results. TEA is written in Python in a modular format. There is a start guide, a user manual, and a code document in addition to this theory paper. TEA is available under a reproducible-research, open-source license via https://github.com/dzesmin/TEA.

[Abridged] We have only been able to comprehensively characterize the atmospheres of a handful of transiting planets, because most orbit faint stars. TESS will discover transiting planets orbiting the brightest stars, enabling, in principle, an atmospheric survey of 10^2 to 10^3 bright hot Jupiters and warm sub-Neptunes. Uniform observations of such a statistically significant sample would provide leverage to understand---and learn from---the diversity of short-period planets. We argue that the best way to maximize the scientific returns of TESS is with a follow-up space mission consisting of a ~1 m telescope with an optical--NIR spectrograph: it could measure molecular absorption for non-terrestrial planets, as well as eclipses and phase variations for the hottest jovians. Such a mission could observe up to 10^3 transits per year, thus enabling it to survey a large fraction of the bright (J<11) TESS planets. JWST could be used to perform detailed atmospheric characterization of the most interesting transiting targets (transit, eclipse, and---when possible---phase-resolved spectroscopy). TESS is also expected to discover a few temperate terrestrial planets transiting nearby M-Dwarfs. Characterizing these worlds will be time-intensive: JWST will need months to provide tantalizing constraints on the presence of an atmosphere, planetary rotational state, clouds, and greenhouse gases. Future flagship missions should be designed to provide better constraints on the habitability of M-Dwarf temperate terrestrial planets.

With an equilibrium temperature of 1200 K, TrES-1 is one of the coolest hot Jupiters observed by {\Spitzer}. It was also the first planet discovered by any transit survey and one of the first exoplanets from which thermal emission was directly observed. We analyzed all {\Spitzer} eclipse and transit data for TrES-1 and obtained its eclipse depths and brightness temperatures in the 3.6 {\micron} (0.083 % {\pm} 0.024 %, 1270 {\pm} 110 K), 4.5 {\micron} (0.094 % {\pm} 0.024 %, 1126 {\pm} 90 K), 5.8 {\micron} (0.162 % {\pm} 0.042 %, 1205 {\pm} 130 K), 8.0 {\micron} (0.213 % {\pm} 0.042 %, 1190 {\pm} 130 K), and 16 {\micron} (0.33 % {\pm} 0.12 %, 1270 {\pm} 310 K) bands. The eclipse depths can be explained, within 1$\sigma$ errors, by a standard atmospheric model with solar abundance composition in chemical equilibrium, with or without a thermal inversion. The combined analysis of the transit, eclipse, and radial-velocity ephemerides gives an eccentricity $e = 0.033^{+0.015}_{-0.031}$, consistent with a circular orbit. Since TrES-1's eclipses have low signal-to-noise ratios, we implemented optimal photometry and differential-evolution Markov-chain Monte Carlo (MCMC) algorithms in our Photometry for Orbits, Eclipses, and Transits (POET) pipeline. Benefits include higher photometric precision and \sim10 times faster MCMC convergence, with better exploration of the phase space and no manual parameter tuning.

WASP-12b was the first planet reported to have a carbon-to-oxygen ratio (C/O) greater than one in its dayside atmosphere. However, recent work to further characterize its atmosphere and confirm its composition has led to incompatible measurements and divergent conclusions. Additionally, the recent discovery of stellar binary companions ~1" from WASP-12 further complicates the analyses and subsequent interpretations. We present a uniform analysis of all available Hubble and Spitzer Space Telescope secondary-eclipse data, including previously-unpublished Spitzer measurements at 3.6 and 4.5 microns. The primary controversy in the literature has centered on the value and interpretation of the eclipse depth at 4.5 microns. Our new measurements and analyses confirm the shallow eclipse depth in this channel, as first reported by Campo and collaborators and used by Madhusudhan and collaborators to infer a carbon-rich composition. To explain WASP-12b's observed dayside emission spectrum, we implemented several recent retrieval approaches. We find that when we exclude absorption due to C2H2 and HCN, which are not universally considered in the literature, our models require implausibly large atmospheric CO2 abundances, regardless of the C/O. By including C2H2 and HCN in our models, we find that a physically-plausible carbon-rich solution achieves the best fit to the available photometric and spectroscopic data. In comparison, the best-fit oxygen-rich models have abundances that are inconsistent with the chemical equilibrium expectations for hydrogen-dominated atmospheres and are 670 times less probable. Our best-fit solution is also 7.3*10^{6} times more probable than an isothermal blackbody model.

WASP-43b is one of the closest-orbiting hot Jupiters, with a semimajor axis of a = 0.01526 +/- 0.00018 AU and a period of only 0.81 days. However, it orbits one of the coolest stars with a hot Jupiter (Tstar = 4520 +/- 120 K), giving the planet a modest equilibrium temperature of Teq = 1440 +/- 40 K, assuming zero Bond albedo and uniform planetary energy redistribution. The eclipse depths and brightness temperatures from our jointly fit model are 0.347% +/- 0.013% and 1670 +/- 23 K at 3.6 {\mu}m and 0.382% +/- 0.015% and 1514 +/- 25 K at 4.5 {\mu}m. The eclipse timings improved the estimate of the orbital period, P, by a factor of three (P = 0.81347436 +/- 1.4*10-7 days) and put an upper limit on the eccentricity (e = 0.010+0.010 -0.007). We use our Spitzer eclipse depths along with four previously reported ground-based photometric observations in the near-infrared to constrain the atmospheric properties of WASP-43b. The data rule out a strong thermal inversion in the dayside atmosphere of WASP-43b. Model atmospheres with no thermal inversions and fiducial oxygen-rich compositions are able to explain all the available data. However, a wide range of metallicities and C/O ratios can explain the data. The data suggest low day-night energy redistribution in the planet, consistent with previous studies, with a nominal upper limit of about 35% for the fraction of energy incident on the dayside that is redistributed to the nightside.

Exoplanet WASP-14b is a highly irradiated, transiting hot Jupiter. Joshi et al. calculate an equilibrium temperature Teq of 1866 K for zero albedo and reemission from the entire planet, a mass of 7.3 +/- 0.5 Jupiter masses and a radius of 1.28 +/- 0.08 Jupiter radii. Its mean density of 4.6 g/cm3 is one of the highest known for planets with periods less than 3 days. We obtained three secondary eclipse light curves with the Spitzer Space Telescope. The eclipse depths from the best jointly fit model are $0.224\%$ +/- $0.018\%$ at 4.5 {\mu}m and $0.181\%$ +/- $0.022\%$ at 8.0 {\mu}m. The corresponding brightness temperatures are 2212 +/- 94 K and 1590 +/- 116 K. A slight ambiguity between systematic models suggests a conservative 3.6 {\mu}m eclipse depth of $0.19\%$ +/- $0.01\%$ and brightness temperature of 2242 +/- 55 K. Although extremely irradiated, WASP-14b does not show any distinct evidence of a thermal inversion. In addition, the present data nominally favor models with day night energy redistribution less than $~30\%$. The current data are generally consistent with oxygen-rich as well as carbon-rich compositions, although an oxygen-rich composition provides a marginally better fit. We confirm a significant eccentricity of e = 0.087 +/- 0.002 and refine other orbital parameters.

Transit-spectroscopy of exoplanets is one of the key observational techniques to characterize the extrasolar planet and its atmosphere. The observational challenges of these measurements require dedicated instrumentation and only the space environment allows an undisturbed access to earth-like atmospheric features such as water or carbon-dioxide. Therefore, several exoplanet-specific space missions are currently being studied. One of them is EChO, the Exoplanet Characterization Observatory, which is part of ESA's Cosmic Vision 2015-2025 program, and which is one of four candidates for the M3 launch slot in 2024. In this paper we present the results of our assessment study of the EChO spectrometer, the only science instrument onboard this spacecraft. The instrument is a multi-channel all-reflective dispersive spectrometer, covering the wavelength range from 400 nm to 16 microns simultaneously with a moderately low spectral resolution. We illustrate how the key technical challenge of the EChO mission - the high photometric stability - influences the choice of spectrometer concept and drives fundamentally the instrument design. First performance evaluations underline the fitness of the elaborated design solution for the needs of the EChO mission.

WASP-8b has 2.18 times Jupiter's mass and is on an eccentric ($e=0.31$) 8.16-day orbit. With a time-averaged equilibrium temperature of 948 K, it is one of the least-irradiated hot Jupiters observed with the Spitzer Space Telescope. We have analyzed six photometric light curves of WASP-8b during secondary eclipse observed in the 3.6, 4.5, and 8.0 {\microns} Infrared Array Camera bands. The eclipse depths are $0.113\pm 0.018$%, $0.069\pm 0.007$%, and $0.093\pm 0.023$%, respectively, giving respective brightness temperatures of 1552, 1131, and 938 K. We characterized the atmospheric thermal profile and composition of the planet using a line-by-line radiative transfer code and a Markov-chain Monte Carlo sampler. The data indicated no thermal inversion, independently of any assumption about chemical composition. We noted an anomalously high 3.6-{\microns} brightness temperature (1552 K); by modeling the eccentricity-caused thermal variation, we found that this temperature is plausible for radiative time scales less than $\sim 10^{2}$ hours. However, as no model spectra fit all three data points well, the temperature discrepancy remains as an open question.

We present new warm Spitzer occultation photometry of WASP-26 at 3.6{\mu}m and 4.5{\mu}m along with new transit photometry taken in the g,r and i bands. We report the first detection of the occultation of WASP-26b, with occultation depths at 3.6{\mu}m and 4.5{\mu}m of 0.00126 +/- 0.00013 and 0.00149 +/- 0.00016 corresponding to brightness temperatures of 1825+/-80K and 1725+/-89K, respectively. We find that the eccentricity of the orbit is consistent with a circular orbit at the 1{\sigma} level with a 3{\sigma} upper limit of e < 0.04. According to the activity-inversion relation of Knutson et al. (2010), WASP-26b is predicted to host a thermal inversion. The brightness temperatures deduced from the eclipse depths are consistent with an isothermal atmosphere, although it is within the uncertainties that the planet may host a weak thermal inversion. The data are equally well fit by atmospheric models with or without a thermal inversion. We find that variation in activity of solar-like stars does not change enough over the time-scales of months or years to change the interpretation of the Knutson et al. (2010) activity-inversion relation, provided that the measured activity level is averaged over several nights. Further data are required to fully constrain the thermal structure of the atmosphere because the planet lies very close to the boundary between atmospheres with and without a thermal inversion.

We report detection of thermal emission from the exoplanet WASP-19b at 3.6, 4.5, 5.8 and 8.0 micron. We used the InfraRed Array Camera on the Spitzer Space Telescope to observe two occultations of WASP-19b by its host star. We combine our new detections with previous measurements of WASP-19b's emission at 1.6 and 2.09 micron to construct a spectral energy distribution of the planet's dayside atmosphere. By comparing this with model-atmosphere spectra, we find that the dayside atmosphere of WASP-19b lacks a strong temperature inversion. As WASP-19 is an active star (log RHK = -4.50 +/- 0.03), this finding supports the hypothesis of Knutson, Howard & Isaacson (2010) that inversions are suppressed in hot Jupiters orbiting active stars. The available data are unable to differentiate between a carbon-rich and an oxygen-rich atmosphere.

We present new lightcurves of the massive hot Jupiter system WASP-18 obtained with the Spitzer spacecraft covering the entire orbit at 3.6 micron and 4.5 micron. These lightcurves are used to measure the amplitude, shape and phase of the thermal phase effect for WASP-18b. We find that our results for the thermal phase effect are limited to an accuracy of about 0.01% by systematic noise sources of unknown origin. At this level of accuracy we find that the thermal phase effect has a peak-to-peak amplitude approximately equal to the secondary eclipse depth, has a sinusoidal shape and that the maximum brightness occurs at the same phase as mid-occultation to within about 5 degrees at 3.6 micron and to within about 10 degrees at 4.5 micron. The shape and amplitude of the thermal phase curve imply very low levels of heat redistribution within the atmosphere of the planet. We also perform a separate analysis to determine the system geometry by fitting a lightcurve model to the data covering the occultation and the transit. The secondary eclipse depths we measure at 3.6 micron and 4.5 micron are in good agreement with previous measurements and imply a very low albedo for WASP-18b. The parameters of the system (masses, radii, etc.) derived from our analysis are in also good agreement with those from previous studies, but with improved precision. We use new high-resolution imaging and published limits on the rate of change of the mean radial velocity to check for the presence of any faint companion stars that may affect our results. We find that there is unlikely to be any significant contribution to the flux at Spitzer wavelengths from a stellar companion to WASP-18. We find that there is no evidence for variations in the times of eclipse from a linear ephemeris greater than about 100 seconds over 3 years.

Aims. We observe occultations of WASP-24b to measure brightness temperatures and to determine whether or not its atmosphere exhibits a thermal inversion (stratosphere). Methods. We observed occultations of WASP-24b at 3.6 and 4.5 {\mu}m using the Spitzer Space Telescope. It has been suggested that there is a correlation between stellar activity and the presence of inversions, so we analysed existing HARPS spectra in order to calculate log R'HK for WASP-24 and thus determine whether or not the star is chromospherically active. We also observed a transit of WASP-24b in the Str\"{o}mgren u and y bands, with the CAHA 2.2-m telescope. Results. We measure occultation depths of 0.159 \pm 0.013 per cent at 3.6 {\mu}m and 0.202 \pm 0.018 per cent at 4.5 {\mu}m. The corresponding planetary brightness temperatures are 1974 \pm 71 K and 1944 \pm 85 K respectively. Atmosphere models with and without a thermal inversion fit the data equally well; we are unable to constrain the presence of an inversion without additional occultation measurements in the near-IR. We find log R'HK = -4.98 \pm 0.12, indicating that WASP-24 is not a chromospherically active star. Our global analysis of new and previously-published data has refined the system parameters, and we find no evidence that the orbit of WASP-24b is non-circular. Conclusions. These results emphasise the importance of complementing Spitzer measurements with observations at shorter wavelengths to gain a full understanding of hot Jupiter atmospheres.

We report the detection of UCF-1.01, a strong exoplanet candidate with a radius 0.66 +/- 0.04 times that of Earth (R_{\oplus}). This sub-Earth-sized planet transits the nearby M-dwarf star GJ 436 with a period of 1.365862 +/- 8x10^{-6} days. We also report evidence of a 0.65 +/- 0.06 R_{\oplus} exoplanet candidate (labeled UCF-1.02) orbiting the same star with an undetermined period. Using the Spitzer Space Telescope, we measure the dimming of light as the planets pass in front of their parent star to assess their sizes and orbital parameters. If confirmed, UCF-1.01 and UCF-1.02 would be called GJ 436c and GJ 436d, respectively, and would be part of the first multiple-transiting-planet system outside of the Kepler field. Assuming Earth-like densities of 5.515 g/cm^3, we predict both candidates to have similar masses (~0.28 Earth-masses, M_{\oplus}, 2.6 Mars-masses) and surface gravities of ~0.65 g (where g is the gravity on Earth). UCF-1.01's equilibrium temperature (T_{eq}, where emitted and absorbed radiation balance for an equivalent blackbody) is 860 K, making the planet unlikely to harbor life as on Earth. Its weak gravitational field and close proximity to its host star imply that UCF-1.01 is unlikely to have retained its original atmosphere; however, a transient atmosphere is possible if recent impacts or tidal heating were to supply volatiles to the surface. We also present additional observations of GJ 436b during secondary eclipse. The 3.6-micron light curve shows indications of stellar activity, making a reliable secondary eclipse measurement impossible. A second non-detection at 4.5 microns supports our previous work in which we find a methane-deficient and carbon monoxide-rich dayside atmosphere.

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.

The dayside of HD 149026b is near the edge of detectability by the Spitzer Space Telescope. We report on eleven secondary-eclipse events at 3.6, 4.5, 3 x 5.8, 4 x 8.0, and 2 x 16 microns plus three primary-transit events at 8.0 microns. The eclipse depths from jointly-fit models at each wavelength are 0.040 +/- 0.003% at 3.6 microns, 0.034 +/- 0.006% at 4.5 microns, 0.044 +/- 0.010% at 5.8 microns, 0.052 +/- 0.006% at 8.0 microns, and 0.085 +/- 0.032% at 16 microns. Multiple observations at the longer wavelengths improved eclipse-depth signal-to-noise ratios by up to a factor of two and improved estimates of the planet-to-star radius ratio (Rp/Rs = 0.0518 +/- 0.0006). We also identify no significant deviations from a circular orbit and, using this model, report an improved period of 2.8758916 +/- 0.0000014 days. Chemical-equilibrium models find no indication of a temperature inversion in the dayside atmosphere of HD 149026b. Our best-fit model favors large amounts of CO and CO2, moderate heat redistribution (f=0.5), and a strongly enhanced metallicity. These analyses use BiLinearly-Interpolated Subpixel Sensitivity (BLISS) mapping, a new technique to model two position-dependent systematics (intrapixel variability and pixelation) by mapping the pixel surface at high resolution. BLISS mapping outperforms previous methods in both speed and goodness of fit. We also present an orthogonalization technique for linearly-correlated parameters that accelerates the convergence of Markov chains that employ the Metropolis random walk sampler. The electronic supplement contains light-curve files and supplementary figures.

We have investigated the 2009 July impact event on Jupiter using the ZEUS-MP 2 three-dimensional hydrodynamics code. We studied the impact itself and the following plume development. Eight impactors were considered: 0.5 km and 1 km porous (\rho = 1.760 g cm^{-3}) and non-porous (\rho = 2.700 g cm^{-3}) basalt impactors, and 0.5 km and 1 km porous (\rho = 0.600 g cm^{-3}) and non-porous \rho = 0.917 g cm^{-3}) ice impactors. The simulations consisted of these bolides colliding with Jupiter at an incident angle of \theta = 69 degrees from the vertical and with an impact velocity of v = 61.4 km s^{-1}. Our simulations show the development of relatively larger, faster plumes created after impacts involving 1 km diameter bodies. Comparing simulations of the 2009 event with simulations of the Shoemaker-Levy 9 events reveals a difference in plume development, with the higher incident angle of the 2009 impact leading to a shallower terminal depth and a smaller and slower plume. We also studied the amount of dynamical chaos present in the simulations conducted at the 2009 incident angle. Compared to the chaos of the SL9 simulations, where \theta is approximately 45 degrees, we find no significant difference in chaos at the higher 2009 incident angle.

We report on the detection of a transit of the super-Earth 55 Cnc e with warm Spitzer in IRAC's 4.5-micron band. Our MCMC analysis includes an extensive modeling of the systematic effects affecting warm Spitzer photometry, and yields a transit depth of 410 +- 63 ppm, which translates to a planetary radius of 2.08 +- 0.16 R_Earth as measured in IRAC 4.5-micron channel. A planetary mass of 7.81 +- 0.58 M_Earth is derived from an extensive set of radial-velocity data, yielding a mean planetary density of 4.8 +- 1.3 g cm-3. Thanks to the brightness of its host star (V = 6, K = 4), 55 Cnc e is a unique target for the thorough characterization of a super-Earth orbiting around a solar-type star.

The transiting exoplanet WASP-18b was discovered in 2008 by the Wide Angle Search for Planets (WASP) project. The Spitzer Exoplanet Target of Opportunity Program observed secondary eclipses of WASP-18b using Spitzer's Infrared Array Camera (IRAC) in the 3.6 micron and 5.8 micron bands on 2008 December 20, and in the 4.5 micron and 8.0 micron bands on 2008 December 24. We report eclipse depths of 0.30 +/- 0.02%, 0.39 +/- 0.02%, 0.37 +/- 0.03%, 0.41 +/- 0.02%, and brightness temperatures of 3100 +/- 90, 3310 +/- 130, 3080 +/- 140 and 3120 +/- 110 K in order of increasing wavelength. WASP-18b is one of the hottest planets yet discovered - as hot as an M-class star. The planet's pressure-temperature profile most likely features a thermal inversion. The observations also require WASP-18b to have near-zero albedo and almost no redistribution of energy from the day-side to the night side of the planet.

We report the detection of thermal emission at 4.5 and 8 micron from the planet WASP-17b. We used Spitzer to measure the system brightness at each wavelength during two occultations of the planet by its host star. By combining the resulting light curves with existing transit light curves and radial velocity measurements in a simultaneous analysis, we find the radius of WASP-17b to be 2.0 Rjup, which is 0.2 Rjup larger than any other known planet and 0.7 Rjup larger than predicted by the standard cooling theory of irradiated gas giant planets. We find the retrograde orbit of WASP-17b to be slightly eccentric, with 0.0012 < e < 0.070 (3 sigma). Such a low eccentricity suggests that, under current models, tidal heating alone could not have bloated the planet to its current size, so the radius of WASP-17b is currently unexplained. From the measured planet-star flux-density ratios we infer 4.5 and 8 micron brightness temperatures of 1881 +/- 50 K and 1580 +/- 150 K, respectively, consistent with a low-albedo planet that efficiently redistributes heat from its day side to its night side.

The nearby extrasolar planet GJ 436b--which has been labelled as a 'hot Neptune'--reveals itself by the dimming of light as it crosses in front of and behind its parent star as seen from Earth. Respectively known as the primary transit and secondary eclipse, the former constrains the planet's radius and mass, and the latter constrains the planet's temperature and, with measurements at multiple wavelengths, its atmospheric composition. Previous work using transmission spectroscopy failed to detect the 1.4-\mu m water vapour band, leaving the planet's atmospheric composition poorly constrained. Here we report the detection of planetary thermal emission from the dayside of GJ 436b at multiple infrared wavelengths during the secondary eclipse. The best-fit compositional models contain a high CO abundance and a substantial methane (CH4) deficiency relative to thermochemical equilibrium models for the predicted hydrogen-dominated atmosphere. Moreover, we report the presence of some H2O and traces of CO2. Because CH4 is expected to be the dominant carbon-bearing species, disequilibrium processes such as vertical mixing and polymerization of methane into substances such as ethylene may be required to explain the hot Neptune's small CH4-to-CO ratio, which is at least 10^5 times smaller than predicted.

We have studied plume formation after a Jovian comet impact using the ZEUS-MP 2 hydrodynamics code. The three-dimensional models followed objects with 500, 750, and 1000 meter diameters. Our simulations show the development of a fast, upward-moving component of the plume in the wake of the impacting comet that "pinches off" from the bulk of the cometary material ~50 km below the 1 bar pressure level, ~100 km above the depth of greatest mass and energy deposition. The fast-moving component contains about twice the mass of the initial comet, but consists almost entirely (>99.9%) of Jovian atmosphere rather than cometary material. The ejecta rise mainly along the impact trajectory, but an additional vertical velocity component due to buoyancy establishes itself within seconds of impact, leading to an asymmetry in the ejecta with respect to the entry trajectory. The mass of the upward-moving component follows a velocity distribution M(>v) approximately proportional to v^-1.4 (v^-1.6 for the 750 m and 500 m cases) in the velocity range 0.1 < v < 10 km/s.

We observed two secondary eclipses of the exoplanet WASP-12b using the Infrared Array Camera on the Spitzer Space Telescope. The close proximity of WASP-12b to its G-type star results in extreme tidal forces capable of inducing apsidal precession with a period as short as a few decades. This precession would be measurable if the orbit had a significant eccentricity, leading to an estimate of the tidal Love number and an assessment of the degree of central concentration in the planetary interior. An initial ground-based secondary eclipse phase reported by Lopez-Morales et al. (0.510 +/- 0.002) implied eccentricity at the 4.5 sigma level. The spectroscopic orbit of Hebb et al. has eccentricity 0.049 +/- 0.015, a 3 sigma result, implying an eclipse phase of 0.509 +/- 0.007. However, there is a well documented tendency of spectroscopic data to overestimate small eccentricities. Our eclipse phases are 0.5010 +/- 0.0006 (3.6 and 5.8 microns) and 0.5006 +/- 0.0007 (4.5 and 8.0 microns). An unlikely orbital precession scenario invoking an alignment of the orbit during the Spitzer observations could have explained this apparent discrepancy, but the final eclipse phase of Lopez-Morales et al. (0.510 -0.006 / +0.007) is consistent with a circular orbit at better than 2 sigma. An orbit fit to all the available transit, eclipse, and radial-velocity data indicates precession at <1 sigma; a non-precessing solution fits better. We also comment on analysis and reporting for Spitzer exoplanet data in light of recent re-analyses.

The carbon-to-oxygen ratio (C/O) in a planet provides critical information about its primordial origins and subsequent evolution. A primordial C/O greater than 0.8 causes a carbide-dominated interior as opposed to a silicate-dominated composition as found on Earth, and the atmospheres can also differ from those in the Solar System. The solar C/O is 0.54. Here we report an analysis of dayside multi-wavelength photometry of the transiting hot-Jupiter WASP-12b that reveals C/O >= 1 in its atmosphere. The atmosphere is abundant in CO. It is depleted in water vapor and enhanced in methane by over two orders of magnitude each compared to a solar-abundance chemical-equilibrium model at the expected temperatures. We also find that the extremely irradiated atmosphere (T > 2,500 K) of WASP-12b lacks a prominent thermal inversion, or a stratosphere, and has very efficient day- night energy circulation. The absence of a strong thermal inversion is in stark contrast to theoretical predictions for the most highly irradiated hot-Jupiter atmospheres.

We report the detection of 24 micron variations from the planet-hosting upsilon Andromedae system consistent with the orbital periodicity of the system's innermost planet, upsilon And b. We find a peak-to-valley phase curve amplitude of 0.00130 times the mean system flux. Using a simple model with two hemispheres of constant surface brightness and assuming a planetary radius of 1.3 Jupiter radii gives a planetary temperature contrast of >900 K and an orbital inclination of >28 degrees. We further report the largest phase offset yet observed for an extrasolar planet: the flux maximum occurs ~80 degrees before phase 0.5. Such a large phase offset is difficult to reconcile with most current atmospheric circulation models. We improve on earlier observations of this system in several important ways: (1) observations of a flux calibrator star demonstrate the MIPS detector is stable to 10^-4 on long timescales, (2) we note that the background light varies systematically due to spacecraft operations, precluding use of this background as a flux calibrator (stellar flux measured above the background is not similarly affected), and (3) we calibrate for flux variability correlated with motion of the star on the MIPS detector. A reanalysis of our earlier observations of this system is consistent with our new result.

On 1998 November 14, Saturn and its rings occulted the star GSC 0622-00345. The occultation latitude was 55.5 degrees S. This paper analyzes the 2.3 {\mu}m light curve derived by Harrington & French. A fixed-baseline isothermal fit to the light curve has a temperature of 140 +/- 3 K, assuming a mean molecular mass of 2.35 AMU. The thermal profile obtained by numerical inversion is valid between 1 and 60 {\mu}bar. The vertical temperature gradient is >0.2 K/km more stable than the adiabatic lapse rate, but it still shows the alternating-rounded-spiked features seen in many temperature gradient profiles from other atmospheric occultations and usually attributed to breaking gravity (buoyancy) waves. We conduct a wavelet analysis of the thermal profile, and show that, even with our low level of noise, scintillation due to turbulence in Earth's atmosphere can produce large temperature swings in light-curve inversions. Spurious periodic features in the "reliable" region of a wavelet amplitude spectrum can exceed 0.3 K in our data. We also show that gravity-wave model fits to noisy isothermal light curves can lead to convincing wave "detections". We provide new significance tests for localized wavelet amplitudes, wave model fits, and global power spectra of inverted occultation light curves by assessing the effects of pre- and post-occultation noise on these parameters. Based on these tests, we detect several significant ridges and isolated peaks in wavelet amplitude, to which we fit a gravity wave model. We also strongly detect the global power spectrum of thermal fluctuations in Saturn's atmosphere, which resembles the "universal" (modified Desaubies) curve associated with saturated spectra of propagating gravity waves on Earth and Jupiter.

On 1998 November 14, Saturn and its rings occulted the star GSC 0622-00345. We observed atmospheric immersion with NSFCAM at the National Aeronautics and Space Administration's Infrared Telescope Facility on Mauna Kea, Hawaii. Immersion occurred at 55.5\circ S planetocentric latitude. A 2.3 {\mu}m, methane-band filter suppressed reflected sunlight. Atmospheric emersion and ring data were not successfully obtained. We describe our observation, light-curve production, and timing techniques, including improvements in aperture positioning, removal of telluric scintillation effects, and timing. Many of these techniques are known within the occultation community, but have not been described in the reviewed literature. We present a light curve whose signal-to-noise ratio per scale height is 267, among the best ground-based signals yet achieved, despite a disadvantage of up to 8 mag in the stellar flux compared to prior work.

We present a comparative study of the thermal emission of the transiting exoplanets WASP-1b and WASP-2b using the Spitzer Space Telescope. The two planets have very similar masses but suffer different levels of irradiation and are predicted to fall either side of a sharp transition between planets with and without hot stratospheres. WASP-1b is one of the most highly irradiated planets studied to date. We measure planet/star contrast ratios in all four of the IRAC bands for both planets (3.6-8.0um), and our results indicate the presence of a strong temperature inversion in the atmosphere of WASP-1b, particularly apparent at 8um, and no inversion in WASP-2b. In both cases the measured eclipse depths favor models in which incident energy is not redistributed efficiently from the day side to the night side of the planet. We fit the Spitzer light curves simultaneously with the best available radial velocity curves and transit photometry in order to provide updated measurements of system parameters. We do not find significant eccentricity in the orbit of either planet, suggesting that the inflated radius of WASP-1b is unlikely to be the result of tidal heating. Finally, by plotting ratios of secondary eclipse depths at 8um and 4.5um against irradiation for all available planets, we find evidence for a sharp transition in the emission spectra of hot Jupiters at an irradiation level of 2 x 10^9 erg/s/cm^2. We suggest this transition may be due to the presence of TiO in the upper atmospheres of the most strongly irradiated hot Jupiters.

We present here the results of our observations of TrES-2 using the Infrared Array Camera on Spitzer. We monitored this transiting system during two secondary eclipses, when the planetary emission is blocked by the star. The resulting decrease in flux is 0.127%+-0.021%, 0.230%+-0.024%, 0.199%+-0.054%, and 0.359%+-0.060%, at 3.6 microns, 4.5 microns, 5.8 microns, and 8.0 microns, respectively. We show that three of these flux contrasts are well fit by a black body spectrum with T_{eff}=1500 K, as well as by a more detailed model spectrum of a planetary atmosphere. The observed planet-to-star flux ratios in all four IRAC channels can be explained by models with and without a thermal inversion in the atmosphere of TrES-2, although with different atmospheric chemistry. Based on the assumption of thermochemical equilibrium, the chemical composition of the inversion model seems more plausible, making it a more favorable scenario. TrES-2 also falls in the category of highly irradiated planets which have been theoretically predicted to exhibit thermal inversions. However, more observations at infrared and visible wavelengths would be needed to confirm a thermal inversion in this system. Furthermore, we find that the times of the secondary eclipses are consistent with previously published times of transit and the expectation from a circular orbit. This implies that TrES-2 most likely has a circular orbit, and thus does not obtain additional thermal energy from tidal dissipation of a non-zero orbital eccentricity, a proposed explanation for the large radius of this planet.

We report measurements of the thermal emission of the young and massive planet CoRoT-2b at 4.5 and 8 microns with the Spitzer Infrared Array Camera (IRAC). Our measured occultation depths are 0.510 +- 0.042 % and 0.41 +- 0.11 % at 4.5 and 8 microns, respectively. In addition to the CoRoT optical measurements, these planet/star flux ratios indicate a poor heat distribution to the night side of the planet and are in better agreement with an atmosphere free of temperature inversion layer. Still, the presence of such an inversion is not definitely ruled out by the observations and a larger wavelength coverage is required to remove the current ambiguity. Our global analysis of CoRoT, Spitzer and ground-based data confirms the large mass and size of the planet with slightly revised values (Mp = 3.47 +- 0.22 Mjup, Rp = 1.466 +- 0.044 Rjup). We find a small but significant offset in the timing of the occultation when compared to a purely circular orbital solution, leading to e cos(omega) = -0.00291 +- 0.00063 where e is the orbital eccentricity and omega is the argument of periastron. Constraining the age of the system to be at most of a few hundreds of Myr and assuming that the non-zero orbital eccentricity is not due to a third undetected body, we model the coupled orbital-tidal evolution of the system with various tidal Q values, core sizes and initial orbital parameters. For log(Q_s') = 5 - 6, our modelling is able to explain the large radius of CoRoT-2b if log(Q_p') <= 5.5 through a transient tidal circularization and corresponding planet tidal heating event. Under this model, the planet will reach its Roche limit within 20 Myr at most.

Over 300 extrasolar planets (exoplanets) have been detected orbiting nearby stars. We now hope to conduct a census of all planets around nearby stars and to characterize their atmospheres and surfaces with spectroscopy. Rocky planets within their star's habitable zones have the highest priority, as these have the potential to harbor life. Our science goal is to find and characterize all nearby exoplanets; this requires that we measure the mass, orbit, and spectroscopic signature of each one at visible and infrared wavelengths. The techniques for doing this are at hand today. Within the decade we could answer long-standing questions about the evolution and nature of other planetary systems, and we could search for clues as to whether life exists elsewhere in our galactic neighborhood.

We report Spitzer/IRAC photometry of the transiting giant exoplanet HAT-P-1b during its secondary eclipse. This planet lies near the postulated boundary between the pM and pL-class of hot Jupiters, and is important as a test of models for temperature inversions in hot Jupiter atmospheres. We derive eclipse depths for HAT-P-1b, in units of the stellar flux, that are: 0.080% +/- 0.008%,[3.6um], 0.135% +/- 0.022%,[4.5um],0.203% +/- 0.031%,[5.8um], and $0.238% +/- 0.040%,[8.0um]. These values are best fit using an atmosphere with a modest temperature inversion, intermediate between the archetype inverted atmosphere (HD209458b) and a model without an inversion. The observations also suggest that this planet is radiating a large fraction of the available stellar irradiance on its dayside, with little available for redistribution by circulation. This planet has sometimes been speculated to be inflated by tidal dissipation, based on its large radius in discovery observations, and on a non-zero orbital eccentricity allowed by the radial velocity data. The timing of the secondary eclipse is very sensitive to orbital eccentricity, and we find that the central phase of the eclipse is 0.4999 +/- 0.0005. The difference between the expected and observed phase indicates that the orbit is close to circular, with a 3-sigma limit of |e cosw| < 0.002.

We report the results of infrared (8 micron) transit and secondary eclipse photometry of the hot Neptune exoplanet, GJ436b using Spitzer. The nearly photon-limited precision of these data allow us to measure an improved radius for the planet, and to detect the secondary eclipse. The transit (centered at HJD = 2454280.78149 +/-0.00016) shows the flat-bottomed shape typical of infrared transits, and it precisely defines the planet-to-star radius ratio (0.0839 +/-0.0005), independent of the stellar properties. However, we obtain the planetary radius, as well as the stellar mass and radius, by fitting to the transit curve simultaneously with an empirical mass-radius relation for M-dwarfs (M=R). We find Rs=Ms=0.47 +/-0.02 in solar units, and Rp=27,600 +/-1170 km (4.33 +/-0.18 Earth radii). This radius significantly exceeds the radius of a naked ocean planet, and requires a gasesous hydrogen-helium envelope. The secondary eclipse occurs at phase 0.587 +/-0.005, proving a significant orbital eccentricity (e=0.15 +/-0.012). The amplitude of the eclipse (5.7 +/-0.8e-4) indicates a brightness temperature for the planet of T=712 +/-36K. If this is indicative of the planet's physical temperature, it suggests the occurrence of tidal heating in the planet. An uncharacterized second planet likely provides ongoing gravitational perturbations, to maintain GJ436b's orbit eccentricity over long time scales.

We report infrared photometry of the extrasolar planet HD 209458b during the time of secondary eclipse (planet passing behind the star). Observations were acquired during two secondary eclipses at the NASA Infrared Telescope Facility (IRTF) in September 2003. We used a circular variable filter (1.5-percent bandpass) centered at 3.8 microns to isolate the predicted flux peak of the planet at this wavelength. Residual telluric absorption and instrument variations were removed by offsetting the telescope to nearby bright comparison stars at a high temporal cadence. Our results give a secondary eclipse depth of 0.0013 +/- 0.0011, not yet sufficient precision to detect the eclipse, whose expected depth is approximately 0.002 - 0.003. We here elucidate the current observational limitations to this technique, and discuss the approach needed to achieve detections of hot Jupiter secondary eclipses at 3.8 microns from the ground.

Of the over 200 known extrasolar planets, 14 exhibit transits in front of their parent stars as seen from Earth. Spectroscopic observations of the transiting planets can probe the physical conditions of their atmospheres. One such technique can be used to derive the planetary spectrum by subtracting the stellar spectrum measured during eclipse (planet hidden behind star) from the combined-light spectrum measured outside eclipse (star + planet). Although several attempts have been made from Earth-based observatories, no spectrum has yet been measured for any of the established extrasolar planets. Here we report a measurement of the infrared spectrum (7.5--13.2 micron) of the transiting extrasolar planet HD209458b. Our observations reveal a hot thermal continuum for the planetary spectrum, with approximately constant ratio to the stellar flux over this wavelength range. Superposed on this continuum is a broad emission peak centered near 9.65 micron that we attribute to emission by silicate clouds. We also find a narrow, unidentified emission feature at 7.78 micron. Models of these ``hot Jupiter'' planets predict a flux peak near 10 micron, where thermal emission from the deep atmosphere emerges relatively unimpeded by water absorption, but models dominated by water fit the observed spectrum poorly.

The star upsilon Andromeda is orbited by three known planets, the innermost of which has an orbital period of 4.617 days and a mass at least 0.69 that of Jupiter. This planet is close enough to its host star that the radiation it absorbs overwhelms its internal heat losses. Here we present the 24 micron light curve of this system, obtained with the Spitzer Space Telescope. It shows a clear variation in phase with the orbital motion of the innermost planet. This is the first demonstration that such planets possess distinct hot substellar (day) and cold antistellar (night) faces.

We have measured the infrared transit of the extrasolar planet HD209458b using the Spitzer Space Telescope. We observed two primary eclipse events (one partial and one complete transit) using the 24 micron array of the Multiband Imaging Photometer for Spitzer (MIPS). We analyzed a total of 2392 individual images (10-second integrations) of the planetary system, recorded before, during, and after transit. We perform optimal photometry on the images and use the local zodiacal light as a short-term flux reference. At this long wavelength, the transit curve has a simple box-like shape, allowing robust solutions for the stellar and planetary radii independent of stellar limb darkening, which is negligible at 24 microns. We derive a stellar radius of R$_*$ = 1.06 $\pm$ 0.07 R$_\sun$, a planetary radius of R$_p$ = 1.26 $\pm$ 0.08 R$_J$, and a stellar mass of 1.17 M$_\sun$. Within the errors, our results agree with the measurements at visible wavelengths. The 24-micron radius of the planet therefore does not differ significantly compared to the visible result. We point out the potential for deriving extrasolar transiting planet radii to high accuracy using transit photometry at slightly shorter IR wavelengths where greater photometric precision is possible.

We have developed a new method that uses wavelet analysis to remove interference fringe patterns from images. This method is particularly useful for flat fields in the common case where fringes vary between the calibration and object data. We analyze the efficacy of this method by creating fake flats with fictitious fringes and removing the fringes. We find that the method removes 90% of the fringe pattern if its amplitude is equal to the random noise level and 60% if the fringe amplitude is $\approx 1/10$ of the noise level. We also present examples using real flat field frames. A routine written in the Interactive Data Language (IDL) that implements this algorithm is available from the authors and as an attachment to this paper.

We have run high-resolution, three-dimensional, hydrodynamic simulations of the impact of comet Shoemaker-Levy 9 into the atmosphere of Jupiter. We find that the energy deposition profile is largely similar to the previous two-dimensional calculations of Mac Low and Zahnle (1994), though perhaps somewhat broader in the range of height over which the energy is deposited. As with similar calculations for impacts into the Venusian atmosphere, there is considerable sensitivity in the results to small changes in the initial conditions, indicating dynamical chaos. We calculated the median depth of energy deposition (the height z at which 50% of the bolide's energy has been released) per run. The mean value among runs is approx 70 km below the 1-bar level, for a 1-km diameter impactor of porous ice of density rho=0.6 g cm^{-3}. The standard deviation among these runs is 14 km. We find little evidence of a trend in these results with the resolution of the calculations (up to 57 cells across the impactor radius, or 8.8-m resolution), suggesting that resolutions as low as 16 grid cells across the radius of the bolide may yield good results for this particular quantity.

We report detection of strong infrared thermal emission from the nearby (d=19 pc) transiting extrasolar planet HD189733b, by measuring the flux decrement during its prominent secondary eclipse. A 6-hour photometric sequence using Spitzer's infrared spectrograph in peak-up imaging mode at 16-microns shows the secondary eclipse depth to be 0.551 +/-0.030%, with accuracy limited by instrumental baseline uncertainties, but with 32-sigma precision (0.017%) on the detection. The 16-micron brightness temperature of this planet (1117+/-42K) is very similar to the Spitzer detections of TrES-1 and HD209458b, but the observed planetary flux (660 micro-Jy) is an order of magnitude greater. This large signal will allow a detailed characterization of this planet in the infrared. Our photometry has sufficient signal-to-noise (~400 per point) to motivate a search for structure in the ingress/egress portions of the eclipse curve, caused by putative thermal structure on the disk of the planet. We show that by binning our 6-second sampling down to 6-minute resolution, we detect the modulation in the intensity derivative during ingress/egress due to the overall shape of the planet, but our sensitivity is not yet sufficient to distinguish between realistic models of the temperature distribution across the planet's disk. We point out the potential for extending Spitzer secondary eclipse detections down to the regime of transiting hot Neptunes, if such systems are discovered among nearby lower main sequence stars.

Recent Spitzer infrared (IR) observations of two transiting hot Jupiters during their secondary eclipses have provided the first direct detection of planets orbiting other stars (Charbonneau et al. 2005; Deming et al. 2005). We here elaborate on some aspects of our detection of HD209458b at 24 microns, and we compare to the detection of TrES-1 by Charbonneau et al. Spitzer will eventually determine the IR spectral energy distribution of these and similar hot Jupiters, opening the new field of comparative exoplanetology. For now, we have only three Spitzer data points, augmented by upper limits from the ground. We here interpret the available measurements from a purely observational perspective, and we point out that a blackbody spectrum having T=1100K can account for all current IR measurements, within the errors. This will surely not remain true for long, since ongoing Spitzer observations will be very sensitive to the IR characteristics of hot Jupiters.

A class of extrasolar giant planets - the so-called `hot Jupiters' - orbit within 0.05 AU of their primary stars. These planets should be hot and so emit detectable infrared radiation. The planet HD 209458b is an ideal candidate for the detection and characterization of this infrared light because it is eclipsed by the star. This planet has an anomalously large radius (1.35 times that of Jupiter), which may be the result of ongoing tidal dissipation, but this explanation requires a non-zero orbital eccentricity (~0.03), maintained by interaction with a hypothetical second planet. Here we report detection of infrared (24 micron) radiation from HD 209458b, by observing the decrement in flux during secondary eclipse, when the planet passes behind the star. The planet's 24 micron flux is 55 +/- 10 micro-Jy (1 sigma), with a brightness temperature of 1130 +/- 150 Kelvins, confirming the predicted heating by stellar irradiation. The secondary eclipse occurs at the midpoint between transits of the planet in front of the star (to within +/- 7 min, 1 sigma), which means that a dynamically significant orbital eccentricity is unlikely.

We have revisited the search for carbon monoxide absorption features in transmission during the transit of the extrasolar planet HD 209458b. We acquired 1077 high resolution spectra at 2 microns using NIRSPEC on Keck II during three transits. Our sensitivity is sufficient to test the degree of CO absorption in the first overtone bands during transit, based on plausible models of the planetary atmosphere. We compare to theoretical tangent geometry absorption spectra, computed by adding height-invariant ad hoc temperature pertubations to the model atmosphere of Sudarsky et al., and by treating cloud height as an adjustable parameter. We do not detect CO absorption. Our analysis indicates a weakening similar to the case of sodium, suggesting that a general masking mechanism is at work in the planetary atmosphere. If this masking is provided by high clouds, our analysis defines the maximum cloud top pressure (i.e., minimum height) as a function of the model atmospheric temperature. For the relatively hot model used by Charbonneau et al. to interpret their sodium detection, our CO limit requires cloud tops at or above 3.3 mbar, and these clouds must be opaque at a wavelength of 2 microns. High clouds comprised of submicron-sized particles are already present in some models, but may not provide sufficient opacity to account for our CO result. Cooler model atmospheres, having smaller atmospheric scale heights and lower CO mixing ratios, may alleviate this problem to some extent. However, even models 500K cooler that the Sudarsky et al. model require clouds above the 100 mbar level to be consistent with our observations. Our null result therefore requires that clouds exist at an observable level in the atmosphere of HD 209458b, unless this planet is dramatically colder than current belief.

We report infrared spectroscopy of the unusual eclipsing pre-main sequence object KH 15D, obtained using NIRSPEC on Keck II. During eclipse, observations using low spectral resolution (1000) reveal the presence of prominent molecular hydrogen emission in 5 lines near 2 microns. The relative line strengths are consistent with thermal excitation at T = 2800 +/-300K. Observations out of eclipse, at both low and high spectral resolution (20,000), show reduced contrast with the stellar continuum. The change in contrast for the strongest line, 1-0 S(1), is consistent with an approximately constant emission line superposed on a variable stellar continuum. Emission in the 1-0 S(1) line is observed to extend by more than 4 arc-sec both east and west of the stellar point spread function (more than 3000 AU). Observed at high spectral resolution, the velocity and intensity structure of the 1-0 S(1) profile are both asymmetric. East of the stellar PSF (by 1.1 - 2.3 arc-sec) the emission is blueshifted (-63 km/sec), and has significantly greater intensity than the marginally redshifted component (+2 km/sec, consistent with zero) which dominates west of the stellar PSF. The spatial extent of the emission, and the excitation temperature, suggest shock-excitation of ambient gas by a bipolar outflow from the star and/or disk. However, it is difficult to account for the observed radial velocity unless the outflow axis is inclined significantly to the plane of the sky.

We search for an infrared signature of the transiting extrasolar planet HD 209458b during secondary eclipse. Our method, which we call `occultation spectroscopy,' searches for the disappearance and reappearance of weak spectral features due to the exoplanet as it passes behind the star and later reappears. We argue that at the longest infrared wavelengths, this technique becomes preferable to conventional `transit spectroscopy'. We observed the system in the wing of the strong nu-3 band of methane near 3.6 microns during two secondary eclipses, using the VLT/ISAAC spectrometer at a spectral resolution of 3300. Our analysis, which utilizes a model template spectrum, achieves sufficient precision to expect detection of the spectral structure predicted by an irradiated, low-opacity (cloudless), low-albedo, thermochemical equilibrium model for the exoplanet atmosphere. However, our observations show no evidence for the presence of this spectrum from the exoplanet, with the statistical significance of the non-detection depending on the timing of the secondary eclipse, which depends on the assumed value for the orbital eccentricity. Our results reject certain specific models of the atmosphere of HD 209458b as inconsistent with our observations at the 3-sigma level, given assumptions about the stellar and planetary parameters.

We present new HST/WFPC2 imagery and STIS long-slit spectroscopy of the planetary nebula NGC 7009. The primary goal was to obtain high spatial resolution of the intrinsic line ratio [O III] 4364/5008 and thereby evaluate the electron temperature (T_e) and the fractional mean-square T_e variation (t_A^2) across the nebula. The WFPC2 T_e map is rather uniform; almost all values are between 9000 - 11,000 K, with the higher T_e's closely coinciding with the inner He^++-zone. The results indicate very small values - <~ 0.01 - for t_A^2 throughout. Our STIS data allow an even more direct determination of T_e and t_A^2, albeit for a much smaller area than with WFPC2. We present results from binning the data along the slit into tiles that are 0.5'' square (matching the slit width). The average [O III] temperature using 45 tiles (excluding the central star and STIS fiducial bars) is 10,139 K; t_A^2 is 0.0035. The measurements of T_e reported here are an average along each line of sight. Therefore, despite finding remarkably low t_A^2, we cannot completely rule out temperature fluctuations along the line of sight as the cause of the large abundance discrepancy between heavy element abundances inferred from collisionally excited emission lines compared to those derived from recombination lines.

The WFPC2 camera aboard the HST was used to obtain images of the planetary nebula BD+30 3639 at two epochs separated by 5.663 years. The expansion of the nebula in the H-alpha and [N II] bands has been measured using several methods. Detailed expansion maps for both emission lines were constructed from nearly 200 almost independent features. There is good agreement between the (independent) H-alpha and [N II] proper motions. The central velocity split is measured from the STIS echelle spectra of the C II] 2326A multiplet to be +/-36.3 km/s at a position angle of 99 degrees. The angular displacement along this slit position was measured to be 4.25 mas/yr at 2".47 from the center. We constructed a tilted ellipsoidal shell model by fitting the radio brightness variation of the 5 and 15 GHz VLA observations, and making use of the ground-based echelle spectra from Bryce & Mellema (1999), to estimate distance. Our model has an axial ratio of 1.56, is inclined to the line of sight by 9.7 degrees, and exhibits an expansion in the plane of the sky which is 2/3 that in the radial direction, leading to a distance of 1.2 kpc. Not all the kinematic data fits this simple model, so the distance must still be regarded as uncertain. Based on the recent model atmosphere of Crowther et al. (2002), a distance of 1.2 kpc implies a stellar luminosity of 4250 L_sun. The kinematic age of the nebula varies somewhat from region to region; a good average value is 800 years, while the expansion along the position of the echelle slit gives about 600 years.

Spectroscopy in the infrared provides a means to assess important properties of the plasma in gaseous nebulae. We present some of our own work that illustrates the need for interactions between the themes of this conference - astronomical data, atomic data, and plasma simulations. We undertook Infrared Space Observatory (ISO) observations with the intent of better understanding the effects of density variations in nebulae, particularly planetary nebulae (PNs), by determining average electron densities from the flux ratios of several fine-structure, IR emission lines. Instead, we are able to ascertain only minor density information because of several instances of the observed line flux ratios being out of range of the theoretical predictions using current atomic data. In these cases, the ISO data cannot presently be used to derive electron density, but rather provide direction for needed improvements in the atomic collision strengths. We have detected an unidentified (uid) strong emission line in an ISO/SWS spectrum of the Orion Nebula. The line has a rest wavelength 2.89350$\pm$0.00003 $\mu$m. A long-slit UKIRT observation confirms the presence of this line and shows that the emission is spatially extended and appears to be coincident with the brightest part of the ionized region. We do not detect the uid line in our SWS02 spectra of any of the several bright PNs which we observed for a comparable time. The need for basic atomic data, in this case wavelengths to aid species identification, is paramount for future progress. We look toward the future with a brief synopsis of upcoming or planned IR missions.

We model the Comet Shoemaker-Levy 9 - Jupiter impact plumes to calculate synthetic plume views, atmospheric infall fluxes, and debris patterns. Our plume is a swarm of ballistic particles with one of several mass-velocity distributions (MVD). The swarm is ejected instantaneously and uniformly into a cone from its apex. Upon falling to the ejection altitude, particles slide with horizontal deceleration following one of several schemes. The model ignores hydrodynamic and Coriolis effects. We adjust plume tilt, opening angle, and minimum velocity, and choose MVD and sliding schemes, to create impact patterns that match observations. Our best match uses the power-law MVD from the numerical impact model of Zahnle and Mac Low, with velocity cutoffs at 4.5 and 11.8 km/sec, cone opening angle of 75 degrees, cone tilt of 30 degrees from vertical, and a sliding constant deceleration of 1.74 m/sec^2. A mathematically-derived feature of Zahnle and Mac Low's published cumulative MVD is a thin shell of mass at the maximum velocity, corresponding to the former atmospheric shock front. This vanguard contains 22% of the mass and 45% of the energy of the plume, and accounts for several previously-unexplained observations, including the large, expanding ring seen at 3.2 microns by McGregor et al. and the ``third precursors'' and ``flare'' seen near 300 and 1000 sec, respectively, in the infrared lightcurves. We present synthetic views of the plumes in flight and after landing and derive infall fluxes of mass, energy, and vertical momentum as a function of time and position on the surface. These fluxes initialize a radiative-hydrodynamic atmosphere model (Paper II of this series) that calculates the thermal and dynamical response of the atmosphere and produces synthetic lightcurves.

We model the plume "splashback" phase of the SL9 collisions with Jupiter using the ZEUS-3D hydrodynamic code. We modified the Zeus code to include gray radiative transport, and we present validation tests. We couple the infalling mass and momentum fluxes of SL9 plume material (from paper I) to a jovian atmospheric model. A strong and complex shock structure results. The modeled shock temperatures agree well with observations, and the structure and evolution of the modeled shocks account for the appearance of high excitation molecular line emission after the peak of the continuum light curve. The splashback region cools by radial expansion as well as by radiation. The morphology of our synthetic continuum light curves agree with observations over a broad wavelength range (0.9 to 12 microns). A feature of our ballistic plume is a shell of mass at the highest velocities, which we term the "vanguard". Portions of the vanguard ejected on shallow trajectories produce a lateral shock front, whose initial expansion accounts for the "third precursors" seen in the 2-micron light curves of the larger impacts, and for hot methane emission at early times. Continued propagation of this lateral shock approximately reproduces the radii, propagation speed, and centroid positions of the large rings observed at 3-4 microns by McGregor et al. The portion of the vanguard ejected closer to the vertical falls back with high z-component velocities just after maximum light, producing CO emission and the "flare" seen at 0.9 microns. The model also produces secondary maxima ("bounces") whose amplitudes and periods are in agreement with observations.

The region surrounding the Llano de Chajnantor, a high altitude plateau in the Atacama Desert in northern Chile, has caught the attention of the astronomical community for its potential as an observatory site. Combining high elevation and extremely low atmospheric water content, the Llano has been chosen as the future site of the Atacama Large Millimeter Array. We have initiated a campaign to investigate the astronomical potential of the region in the optical/infrared. Here, we report on an aspect of our campaign aimed at establishing a seeing benchmark to be used as a reference for future activities in the region. After a brief description of the region and its climate, we describe the results of an astronomical seeing campaign, carried out with a Differential Image Motion Monitor that operates at 0.5 micron wavelength. The seeing at the Llano level of 5000 m, measured over 7 nights in May 1998, yielded a median FWHM of 1.1". However, the seeing decreased to 0.7" at a modest 100 m gain above the plateau (Cerro Chico), as measured over 38 nights spread between July 1998 and October 2000. Neither of these represents the best seeing expected in the region; the set of measurements provides a reference base for simultaneous dual runs at Cerro Chico and at other sites of interest in the region, currently underway. A comparison between simultaneous measurements at Cerro Chico and Cerro Paranal indicates that seeing at Cerro Chico is about 12% better than at Paranal. The percentage of optically photometric nights in the Chajnantor region is about 60%, while that of nights useful for astronomical work is near 80%.

Spectroscopic observations of a proto-planetary nebula He 3-1475 with the Space Telescope Imaging Spectrograph (STIS) reveal the kinematics of its high (1200 km/s) velocity jets. The jets are formed at a large (0.15 pc) distance from its central star by collimation of an asymmetric stellar wind in a pair of conical shocks seen in Wide Field Planetary Camera (WFPC2) images. The jets consist of several pairs of knots symmetrically distributed with respect to the central star, with most knots exhibiting a head-tail morphology. Large (up to 650 km/s) radial velocity gradients are seen within the knots on subarcsec scales, with velocities decreasing from the knot heads toward their trailing tails. These large velocity gradients are a sign of efficient deceleration of jets by a much slower bipolar outflow. The inclination angle of the bipolar outflow is equal to 40 degrees, based on Doppler shifts of the scattered stellar Halpha line. Its velocity is equal to 140 km/s at a distance of 0.23 pc from the star, and increases monotonically with the radial distance from the star. A comparison of new WFPC2 [N II] 6584 images with older WFPC2 images reveals expansion of the jets. The measured jet proper motions in combination with their radial velocities imply that He 3-1475 is a Galactic Bulge star at a distance of 8 kpc, located 800 pc above the Galactic plane. Its very high luminosity (25,000 solar luminosities) implies that He 3-1475 must be significantly more massive than a typical AGB star within the Galactic Bulge, perhaps because of a past mass transfer and/or a merger event in an interacting binary system.

The jets seen in NGC 6543 motivated us to search for similar features in other PNe. Our HST snapshot program looked for jets using the [N II] and H-alpha narrow-band filters. Although spectacular jets are found in proto-PNe, true jets seem rare among mature PNe. In the later group, IC 4593 is the best case in our sample. We distinguish between jets and a number of interesting ``jet-like'' features, e.g., ``cometary'' structures with a dense globule at the end facing the central star.