With more than 1000 hours of observation from Feb 2016 to Oct 2019, the Spitzer Exploration Program Red Worlds (ID: 13067, 13175 and 14223) exclusively targeted TRAPPIST-1, a nearby (12pc) ultracool dwarf star orbited by seven transiting Earth-sized planets, all well-suited for a detailed atmospheric characterization with the upcoming JWST. In this paper, we present the global results of the project. We analyzed 88 new transits and combined them with 100 previously analyzed transits, for a total of 188 transits observed at 3.6 or 4.5 $\mu$m. We also analyzed 29 occultations (secondary eclipses) of planet b and eight occultations of planet c observed at 4.5 $\mu$m to constrain the brightness temperatures of their daysides. We identify several orphan transit-like structures in our Spitzer photometry, but all of them are of low significance. We do not confirm any new transiting planets. We estimate for TRAPPIST-1 transit depth measurements mean noise floors of $\sim$35 and 25 ppm in channels 1 and 2 of Spitzer/IRAC, respectively. most of this noise floor is of instrumental origins and due to the large inter-pixel inhomogeneity of IRAC InSb arrays, and that the much better interpixel homogeneity of JWST instruments should result in noise floors as low as 10ppm, which is low enough to enable the atmospheric characterization of the planets by transit transmission spectroscopy. We construct updated broadband transmission spectra for all seven planets which show consistent transit depths between the two Spitzer channels. We identify and model five distinct high energy flares in the whole dataset, and discuss our results in the context of habitability. Finally, we fail to detect occultation signals of planets b and c at 4.5 $\mu$m, and can only set 3$\sigma$ upper limits on their dayside brightness temperatures (611K for b 586K for c).

We present a new method employing machine learning techniques for measuring astrophysical features by correcting systematics in IRAC high precision photometry using Random Forests. The main systematic in IRAC light curve data is position changes due to unavoidable telescope motions coupled with an intrapixel response function. We aim to use the large amount of publicly available calibration data for the single pixel used for this type of work (the sweet spot pixel) to make a fast, easy to use, accurate correction to science data. This correction on calibration data has the advantage of using an independent dataset instead of using the science data on itself, which has the disadvantage of including astrophysical variations. After focusing on feature engineering and hyperparameter optimization, we show that a boosted random forest model can reduce the data such that we measure the median of ten archival eclipse observations of XO-3b to be 1459 +- 200 parts per million. This is a comparable depth to the average of those in the literature done by seven different methods, however the spread in measurements is 30-100% larger than those literature values, depending on the reduction method. We also caution others attempting similar methods to check their results with the fiducial dataset of XO-3b as we were also able to find models providing initially great scores on their internal test datasets but whose results significantly underestimated the eclipse depth of that planet.

We present preliminary trigonometric parallaxes of 184 late-T and Y dwarfs using observations from Spitzer (143), USNO (18), NTT (14), and UKIRT (9). To complete the 20-pc census of $\ge$T6 dwarfs, we combine these measurements with previously published trigonometric parallaxes for an additional 44 objects and spectrophotometric distance estimates for another 7. For these 235 objects, we estimate temperatures, sift into five 150K-wide $T_{\rm eff}$ bins covering the range 300-1050K, determine the completeness limit for each, and compute space densities. To anchor the high-mass end of the brown dwarf mass spectrum, we compile a list of early- to mid-L dwarfs within 20 pc. We run simulations using various functional forms of the mass function passed through two different sets of evolutionary code to compute predicted distributions in $T_{\rm eff}$. The best fit of these predictions to our L, T, and Y observations is a simple power-law model with $\alpha \approx 0.6$ (where $dN/dM \propto M^{-\alpha}$), meaning that the slope of the field substellar mass function is in rough agreement with that found for brown dwarfs in nearby star forming regions and young clusters. Furthermore, we find that published versions of the log-normal form do not predict the steady rise seen in the space densities from 1050K to 350K. We also find that the low-mass cutoff to formation, if one exists, is lower than $\sim$5 $M_{Jup}$, which corroborates findings in young, nearby moving groups and implies that extremely low-mass objects have been forming over the lifetime of the Milky Way.

Y dwarfs provide a unique opportunity to study free-floating objects with masses $<$30 M$_{Jup}$ and atmospheric temperatures approaching those of known Jupiter-like exoplanets. Obtaining distances to these objects is an essential step towards characterizing their absolute physical properties. Using Spitzer/IRAC [4.5] images taken over baselines of $\sim$2-7 years, we measure astrometric distances for 22 late-T and early Y dwarfs, including updated parallaxes for 18 objects and new parallax measurements for 4 objects. These parallaxes will make it possible to explore the physical parameter space occupied by the coldest brown dwarfs. We also present the discovery of 6 new late-T dwarfs, updated spectra of two T dwarfs, and the reclassification of a new Y dwarf, WISE J033605.04$-$014351.0, based on Keck/NIRSPEC $J$-band spectroscopy. Assuming that effective temperatures are inversely proportional to absolute magnitude, we examine trends in the evolution of the spectral energy distributions of brown dwarfs with decreasing effective temperature. Surprisingly, the Y dwarf class encompasses a large range in absolute magnitude in the near- to mid-infrared photometric bandpasses, demonstrating a larger range of effective temperatures than previously assumed. This sample will be ideal for obtaining mid-infrared spectra with the James Webb Space Telescope because their known distances will make it easier to measure absolute physical properties.

Sagittarius A* (Sgr A*) is the variable radio, near-infrared (NIR), and X-ray source associated with accretion onto the Galactic center black hole. We present an analysis of the most comprehensive NIR variability dataset of Sgr A* to date: eight 24-hour epochs of continuous monitoring of Sgr A* at 4.5 $\mu$m with the IRAC instrument on the Spitzer Space Telescope, 93 epochs of 2.18 $\mu$m data from Naos Conica at the Very Large Telescope, and 30 epochs of 2.12 $\mu$m data from the NIRC2 camera at the Keck Observatory, in total 94,929 measurements. A new approximate Bayesian computation method for fitting the first-order structure function extracts information beyond current Fast Fourier Transformation (FFT) methods of power spectral density (PSD) estimation. With a combined fit of the data of all three observatories, the characteristic coherence timescale of Sgr A* is $\tau_{b} = 243^{+82}_{-57}$ minutes ($90\%$ credible interval). The PSD has no detectable features on timescales down to 8.5 minutes ($95\%$ credible level), which is the ISCO orbital frequency for a dimensionless spin parameter $a = 0.92$. One light curve measured simultaneously at 2.12 and 4.5 $\mu$m during a low flux-density phase gave a spectral index $\alpha_s = 1.6 \pm 0.1$ ($F_\nu \propto \nu^{-\alpha_s}$). This value implies that the Sgr A* NIR color becomes bluer during higher flux-density phases. The probability densities of flux densities of the combined datasets are best fit by log-normal distributions. Based on these distributions, the Sgr A* spectral energy distribution is consistent with synchrotron radiation from a non-thermal electron population from below 20 GHz through the NIR.

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.

We apply the transit light curve self-contamination technique of Morris et al. (2018) to search for the effect of stellar activity on the transits of the ultracool dwarf TRAPPIST-1 with 2018 Spitzer photometry. The self-contamination method fits the transit light curves of planets orbiting spotted stars, allowing the host star to be a source of contaminating positive or negative flux which influences the transit depths but not the ingress/egress durations. We find that none of the planets show statistically significant evidence for self-contamination by bright or dark regions of the stellar photosphere. However, we show that small-scale magnetic activity, analogous in size to the smallest sunspots, could still be lurking in the transit photometry undetected.

Sgr A*, the supermassive black hole (SMBH) at the center of our Milky Way Galaxy, is known to be a variable source of X-ray, near-infrared (NIR), and submillimeter (submm) radiation and therefore a prime candidate to study the electromagnetic radiation generated by mass accretion flow onto a black hole and/or a related jet. Disentangling the power source and emission mechanisms of this variability is a central challenge to our understanding of accretion flows around SMBHs. Simultaneous multiwavelength observations of the flux variations and their time correlations can play an important role in obtaining a better understanding of possible emission mechanisms and their origin. This paper presents observations of two flares that both apparently violate the previously established patterns in the relative timing of submm/NIR/X-ray flares from Sgr A*. One of these events provides the first evidence of coeval structure between NIR and submm flux increases, while the second event is the first example of the sequence of submm/X-ray/NIR flux increases all occurring within ~1 hr. Each of these two events appears to upend assumptions that have been the basis of some analytic models of flaring in Sgr A*. However, it cannot be ruled out that these events, even though unusual, were just coincidental. These observations demonstrate that we do not fully understand the origin of the multiwavelength variability of Sgr A*, and show that there is a continued and important need for long-term, coordinated, and precise multiwavelength observations of Sgr A* to characterize the full range of variability behavior.

Context. The TRAPPIST-1 system hosts seven Earth-sized, temperate exoplanets orbiting an ultra-cool dwarf star. As such, it represents a remarkable setting to study the formation and evolution of terrestrial planets that formed in the same protoplanetary disk. While the sizes of the TRAPPIST-1 planets are all known to better than 5% precision, their densities have significant uncertainties (between 28% and 95%) because of poor constraints on the planet's masses. Aims.The goal of this paper is to improve our knowledge of the TRAPPIST-1 planetary masses and densities using transit-timing variations (TTV). The complexity of the TTV inversion problem is known to be particularly acute in multi-planetary systems (convergence issues, degeneracies and size of the parameter space), especially for resonant chain systems such as TRAPPIST-1. Methods. To overcome these challenges, we have used a novel method that employs a genetic algorithm coupled to a full N-body integrator that we applied to a set of 284 individual transit timings. This approach enables us to efficiently explore the parameter space and to derive reliable masses and densities from TTVs for all seven planets. Results. Our new masses result in a five- to eight-fold improvement on the planetary density uncertainties, with precisions ranging from 5% to 12%. These updated values provide new insights into the bulk structure of the TRAPPIST-1 planets. We find that TRAPPIST-1\,c and e likely have largely rocky interiors, while planets b, d, f, g, and h require envelopes of volatiles in the form of thick atmospheres, oceans, or ice, in most cases with water mass fractions less than 5%.

The recently detected TRAPPIST-1 planetary system, with its seven planets transiting a nearby ultracool dwarf star, offers the first opportunity to perform comparative exoplanetology of temperate Earth-sized worlds. To further advance our understanding of these planets' compositions, energy budgets, and dynamics, we are carrying out an intensive photometric monitoring campaign of their transits with the $\textit{Spitzer Space Telescope}$. In this context, we present 60 new transits of the TRAPPIST-1 planets observed with $\textit{Spitzer}$/IRAC in February and March 2017. We combine these observations with previously published $\textit{Spitzer}$ transit photometry and perform a global analysis of the resulting extensive dataset. This analysis refines the transit parameters and provides revised values for the planets' physical parameters, notably their radii, using updated properties for the star. As part of our study, we also measure precise transit timings that will be used in a companion paper to refine the planets' masses and compositions using the transit timing variations method. TRAPPIST-1 shows a very low level of low-frequency variability in the IRAC 4.5-$\mu$m band, with a photometric RMS of only 0.11$\%$ at a 123-s cadence. We do not detect any evidence of a (quasi-)periodic signal related to stellar rotation. We also analyze the transit light curves individually, to search for possible variations in the transit parameters of each planet due to stellar variability, and find that the $\textit{Spitzer}$ transits of the planets are mostly immune to the effects of stellar variations. These results are encouraging for forthcoming transmission spectroscopy observations of the TRAPPIST-1 planets with the $\textit{James Webb Space Telescope}$.

Although the final observations of the Spitzer Warm Mission are currently scheduled for March 2019, it can continue operations through the end of the decade with no loss of photometric precision. As we will show, there is a strong science case for extending the current Warm Mission to December 2020. Spitzer has already made major impacts in the fields of exoplanets (including microlensing events), characterizing near Earth objects, enhancing our knowledge of nearby stars and brown dwarfs, understanding the properties and structure of our Milky Way galaxy, and deep wide-field extragalactic surveys to study galaxy birth and evolution. By extending Spitzer through 2020, it can continue to make ground-breaking discoveries in those fields, and provide crucial support to the NASA flagship missions JWST and WFIRST, as well as the upcoming TESS mission, and it will complement ground-based observations by LSST and the new large telescopes of the next decade. This scientific program addresses NASA's Science Mission Directive's objectives in astrophysics, which include discovering how the universe works, exploring how it began and evolved, and searching for life on planets around other stars.

The TRAPPIST-1 system is the first transiting planet system found orbiting an ultra-cool dwarf star. At least seven planets similar to Earth in radius and in mass were previously found to transit this host star. Subsequently, TRAPPIST-1 was observed as part of the K2 mission and, with these new data, we report the measurement of an 18.77 d orbital period for the outermost planet, TRAPPIST-1h, which was unconstrained until now. This value matches our theoretical expectations based on Laplace relations and places TRAPPIST-1h as the seventh member of a complex chain, with three-body resonances linking every member. We find that TRAPPIST-1h has a radius of 0.727 Earth radii and an equilibrium temperature of 173 K. We have also measured the rotational period of the star at 3.3 d and detected a number of flares consistent with a low-activity, middle-aged, late M dwarf.

One focus of modern astronomy is to detect temperate terrestrial exoplanets well-suited for atmospheric characterisation. A milestone was recently achieved with the detection of three Earth-sized planets transiting (i.e. passing in front of) a star just 8% the mass of the Sun 12 parsecs away. Indeed, the transiting configuration of these planets with the Jupiter-like size of their host star - named TRAPPIST-1 - makes possible in-depth studies of their atmospheric properties with current and future astronomical facilities. Here we report the results of an intensive photometric monitoring campaign of that star from the ground and with the Spitzer Space Telescope. Our observations reveal that at least seven planets with sizes and masses similar to the Earth revolve around TRAPPIST-1. The six inner planets form a near-resonant chain such that their orbital periods (1.51, 2.42, 4.04, 6.06, 9.21, 12.35 days) are near ratios of small integers. This architecture suggests that the planets formed farther from the star and migrated inward. The seven planets have equilibrium temperatures low enough to make possible liquid water on their surfaces.

We measure the 4.5 $\mu$m thermal emission of five transiting hot Jupiters, WASP-13b, WASP-15b, WASP-16b, WASP-62b and HAT-P-22b using channel 2 of the Infrared Array Camera (IRAC) on the {\sl Spitzer Space Telescope}. Significant intrapixel sensitivity variations in Spitzer IRAC data require careful correction in order to achieve precision on the order of several hundred parts per million (ppm) for the measurement of exoplanet secondary eclipses. We determine eclipse depths by first correcting the raw data using three independent data reduction methods. The Pixel Gain Map (PMAP), Nearest Neighbors (NNBR), and Pixel Level Decorrelation (PLD) each correct for the intrapixel sensitivity effect in Spitzer photometric time-series observations. The results from each methodology are compared against each other to establish if they reach a statistically equivalent result in every case and to evaluate their ability to minimize uncertainty in the measurement. We find that all three methods produce reliable results. For every planet examined here NNBR and PLD produce results that are in statistical agreement. However, the PMAP method appears to produce results in slight disagreement in cases where the stellar centroid is not kept consistently on the most well characterized area of the detector. We evaluate the ability of each method to reduce the scatter in the residuals as well as in the correlated noise in the corrected data. The NNBR and PLD methods consistently minimize both white and red noise levels and should be considered reliable and consistent. The planets in this study span equilibrium temperatures from 1100 to 2000~K and have brightness temperatures that require either high albedo or efficient recirculation. However, it is possible that other processes such as clouds or disequilibrium chemistry may also be responsible for producing these brightness temperatures.

We have used the Spitzer Space Telescope in February 2016 to obtain high cadence, high signal-to-noise, 17-hour duration light curves of Neptune at 3.6 and 4.5 $\mu$m. The light curve duration was chosen to correspond to the rotation period of Neptune. Both light curves are slowly varying with time, with full amplitudes of 1.1 mag at 3.6 $\mu$m and 0.6 mag at 4.5 $\mu$m. We have also extracted sparsely sampled 18-hour light curves of Neptune at W1 (3.4 $\mu$m) and W2 (4.6 $\mu$m) from the WISE/NEOWISE archive at six epochs in 2010-2015. These light curves all show similar shapes and amplitudes compared to the Spitzer light curves but with considerable variation from epoch to epoch. These amplitudes are much larger than those observed with Kepler/K2 in the visible (amplitude $\sim$0.02 mag) or at 845 nm with the Hubble Space Telescope in 2015 and at 763 nm in 2016 (amplitude $\sim$ 0.2 mag). We interpret the Spitzer and WISE light curves as arising entirely from reflected solar photons, from higher levels in Neptune's atmosphere than for K2. Methane gas is the dominant opacity source in Neptune's atmosphere, and methane absorption bands are present in the HST 763, and 845 nm, WISE W1, and Spitzer 3.6 $\mu$m filters.

We examine the repeatability, reliability, and accuracy of differential exoplanet eclipse depth measurements made using the InfraRed Array Camera (IRAC) on the Spitzer Space Telescope during the post-cryogenic mission. We have re-analyzed an existing 4.5 {\mu}m data set, consisting of 10 observations of the XO-3b system during secondary eclipse, using seven different techniques for removing correlated noise. We find that, on average, for a given technique, the eclipse depth estimate is repeatable from epoch to epoch to within 156 parts per million (ppm). Most techniques derive eclipse depths that do not vary by more than a factor 3 of the photon noise limit. All methods but one accurately assess their own errors: for these methods, the individual measurement uncertainties are comparable to the scatter in eclipse depths over the 10 epoch sample. To assess the accuracy of the techniques as well as to clarify the difference between instrumental and other sources of measurement error, we have also analyzed a simulated data set of 10 visits to XO-3b, for which the eclipse depth is known. We find that three of the methods (BLISS mapping, Pixel Level Decorrelation, and Independent Component Analysis) obtain results that are within three times the photon limit of the true eclipse depth. When averaged over the 10 epoch ensemble, 5 out of 7 techniques come within 60 ppm of the true value. Spitzer exoplanet data, if obtained following current best practices and reduced using methods such as those described here, can measure repeatable and accurate single eclipse depths, with close to photon-limited results.

Motivated by a high Spitzer IRAC oversubscription rate, we present a new technique of randomly and sparsely sampling phase curves of hot Jupiters. Snapshot phase curves are enabled by technical advances in precision pointing as well as careful characterization of a portion of the central pixel on the array. This method allows for observations which are a factor of roughly two more efficient than full phase curve observations, and are furthermore easier to insert into the Spitzer observing schedule. We present our pilot study from this program using the exoplanet WASP-14b. Data of this system were taken both as a sparsely sampled phase curve as well as a staring mode phase curve. Both datasets as well as snapshot style observations of a calibration star are used to validate this technique. By fitting our WASP-14b phase snapshot dataset, we successfully recover physical parameters for the transit and eclipse depths as well as amplitude and maximum and minimum of the phase curve shape of this slightly eccentric hot Jupiter. We place a limit on the potential phase to phase variation of these parameters since our data are taken over many phases over the course of a year. We see no evidence for eclipse depth variations compared to other published WASP-14b eclipse depths over a 3.5 year baseline.

HAT-P-20b is a giant exoplanet that orbits a metal-rich star. The planet itself has a high total density, suggesting that it may also have a high metallicity in its atmosphere. We analyze two eclipses of the planet in each of the 3.6- and 4.5 micron bands of Warm Spitzer. These data exhibit intra-pixel detector sensitivity fluctuations that were resistant to traditional decorrelation methods. We have developed a simple, powerful, and radically different method to correct the intra-pixel effect for Warm Spitzer data, which we call pixel-level decorrelation (PLD). PLD corrects the intra-pixel effect very effectively, but without explicitly using - or even measuring - the fluctuations in the apparent position of the stellar image. We illustrate and validate PLD using synthetic and real data, and comparing the results to previous analyses. PLD can significantly reduce or eliminate red noise in Spitzer secondary eclipse photometry, even for eclipses that have proven to be intractable using other methods. Our successful PLD analysis of four HAT-P-20b eclipses shows a best-fit blackbody temperature of 1134 +/-29K, indicating inefficient longitudinal transfer of heat, but lacking evidence for strong molecular absorption. We find sufficient evidence for variability in the 4.5 micron band that the eclipses should be monitored at that wavelength by Spitzer, and this planet should be a high priority for JWST spectroscopy. All four eclipses occur about 35 minutes after orbital phase 0.5, indicating a slightly eccentric orbit. A joint fit of the eclipse and transit times with extant RV data yields e(cos{omega}) = 0.01352 (+0.00054, -0.00057), and establishes the small eccentricity of the orbit to high statistical confidence. Given the existence of a bound stellar companion, HAT-P-20b is another excellent candidate for orbital evolution via Kozai migration or other three-body mechanism.

We present the first detection from the Spitzer Space Telescope of 4.5 micron variability from Sgr A*, the emitting source associated with the Milky Way's central black hole. The >23 hour continuous light curve was obtained with the IRAC instrument in 2013 December. The result characterizes the variability of Sgr A* prior to the closest approach of the G2 object, a putative infalling gas cloud that orbits close to Sgr A*. The high stellar density at the location of Sgr A* produces a background of ~250 mJy at 4.5 microns in each pixel with a large pixel-to-pixel gradient, but the light curve for the highly variable Sgr A* source was successfully measured by modeling and removing the variations due to pointing wobble. The observed flux densities range from the noise level of ~0.7 mJy rms in a 6.4-s measurement to ~10 mJy. Emission was seen above the noise level ~34% of the time. The light curve characteristics, including the flux density distribution and structure function, are consistent with those previously derived at shorter infrared wavelengths. We see no evidence in the light curve for activity attributable to the G2 interaction at the observing epoch, ~100 days before the expected G2 periapsis passage. The IRAC light curve is more than a factor of two longer than any previous infrared observation, improving constraints on the timescale of the break in the power spectral distribution of Sgr A* flux densities. The data favor the longer of the two previously published values for the timescale.

Herschel FIR continuum data obtained as part of the Hi-GAL survey have been used, together with the GLIMPSE 8 $\mu$m and MIPSGAL 24 $\mu$m data, to attempt the first 3D-decomposition of dust emission associated with atomic, molecular and ionized gas at 15 arcmin angular resolution. Our initial test case is a 2$\times$2 square degrees region centred on (l,b)=(30$^{\circ}$,0$^{\circ}$), a direction that encompasses the origin point of the Scutum-Crux Arm at the tip of the Galactic Bar. Coupling the IR maps with velocity maps specific for different gas phases (HI 21cm, $^{12}$CO and $^{13}$CO, and RRLs), we estimate the properties of dust blended with each of the gas components and at different Galactocentric distances along the LOS. A statistical Pearson's coefficients analysis is used to study the correlation between the column densities and the intensity of the IR emission. This analysis provides evidence that the 2$\times$2 square degree field under consideration is characterized by the presence of a gas component not accounted for by the standard tracers, possibly associated with warm H$_{2}$ and cold HI. We demonstrate that the IR radiation in the range 8 $\mu$m $<$ $\lambda$ $<$ 500 $\mu$m is systematically dominated by emission originating within the Scutum-Crux Arm. By applying an inversion method, we recover the dust emissivities associated with atomic, molecular and ionized gas. Using the DustEM model we obtain an indication for PAHs depletion in the diffuse ionized gas. However, the main goal of this work is to discuss the impact of the missing column density associated with the dark gas component on the accurate evaluation of the dust properties, and to shed light on the limitations of the inversion method approach when this is applied to a small section of the Galactic Plane and when the working resolution allows sufficient de-blending of the gas components along the LOS.

The dominant non-instrumental background source for space-based infrared observatories is the zo- diacal light. We present Spitzer Infrared Array Camera (IRAC) measurements of the zodiacal light at 3.6, 4.5, 5.8, and 8.0 {\mu}m, taken as part of the instrument calibrations. We measure the changing surface brightness levels in approximately weekly IRAC observations near the north ecliptic pole (NEP) over the period of roughly 8.5 years. This long time baseline is crucial for measuring the annual sinusoidal variation in the signal levels due to the tilt of the dust disk with respect to the ecliptic, which is the true signal of the zodiacal light. This is compared to both Cosmic Background Explorer Diffuse Infrared Background Experiment (COBE DIRBE) data and a zodiacal light model based thereon. Our data show a few percent discrepancy from the Kelsall et al. (1998) model including a potential warping of the interplanetary dust disk and a previously detected overdensity in the dust cloud directly behind the Earth in its orbit. Accurate knowledge of the zodiacal light is important for both extragalactic and Galactic astronomy including measurements of the cosmic infrared background, absolute measures of extended sources, and comparison to extrasolar interplanetary dust models. IRAC data can be used to further inform and test future zodiacal light models.

Using the Infrared Spectrograph on board the Spitzer Space Telescope, we have detected emission in the S(0), S(1), and S(2) pure-rotational (v=0-0) transitions of molecular hydrogen (H2) towards 6 positions in two translucent high Galactic latitude clouds, DCld 300.2-16.9 and LDN 1780. The detection of these lines raises important questions regarding the physical conditions inside low-extinction clouds that are far from ultraviolet radiation sources. The ratio between the S(2) flux and the flux from PAHs at 7.9 microns averages 0.007 for these 6 positions. This is a factor of about 4 higher than the same ratio measured towards the central regions of non-active Galaxies in the Spitzer Infrared Nearby Galaxies Survey (SINGS). Thus the environment of these translucent clouds is more efficient at producing rotationally excited H2 per PAH-exciting photon than the disks of entire galaxies. Excitation analysis finds that the S(1) and S(2) emitting regions are warm (T >300K), but comprise no more than 2% of the gas mass. We find that UV photons cannot be the sole source of excitation in these regions and suggest mechanical heating via shocks or turbulent dissipation as the dominant cause of the emission. The clouds are located on the outskirts of the Scorpius-Centaurus OB association and may be dissipating recent bursts of mechanical energy input from supernova explosions. We suggest that pockets of warm gas in diffuse or translucent clouds, integrated over the disks of galaxies, may represent a major source of all non-active galaxy H2 emission.

We present results of mid-infrared spectroscopic mapping observations of six star-forming regions in the Small Magellanic Cloud from the Spitzer Spectroscopic Survey of the SMC (S4MC). We detect the mid-IR emission from polycyclic aromatic hydrocarbons (PAHs) in all of the mapped regions, greatly increasing the range of environments where PAHs have been spectroscopically detected in the SMC. We investigate the variations of the mid-IR bands in each region and compare our results to studies of the PAH bands in the SINGS sample and in a sample of low-metallicity starburst galaxies. PAH emission in the SMC is characterized by low ratios of the 6-9 micron features relative to the 11.3 micron feature and weak 8.6 and 17.0 micron features. Interpreting these band ratios in the light of laboratory and theoretical studies, we find that PAHs in the SMC tend to be smaller and less ionized than those in higher metallicity galaxies. Based on studies of PAH destruction, we argue that a size distribution shifted towards smaller PAHs cannot be the result of processing in the interstellar medium, but instead reflects differences in the formation of PAHs at low metallicity. Finally, we discuss the implications of our observations for our understanding of the PAH life-cycle in low-metallicity galaxies---namely that the observed deficit of PAHs may be a consequence of PAHs forming with smaller average sizes and therefore being more susceptible to destruction under typical interstellar medium conditions.

The fullerene C$_{60}$ has four infrared-active vibrational transitions at 7.0, 8.5, 17.4 and 18.9 $\mu$m. We have previously observed emission features at 17.4 and 18.9 $\mu$m in the reflection nebula NGC 7023 and demonstrated spatial correlations suggestive of a common origin. We now confirm our earlier identification of these features with C$_{60}$ by detecting a third emission feature at 7.04 $\pm$ 0.05 $\mu$m in NGC 7023. We also report the detection of these three C$_{60}$ features in the reflection nebula NGC 2023. Our spectroscopic mapping of NGC 7023 shows that the 18.9 $\mu$m C$_{60}$ feature peaks on the central star and that the 16.4 $\mu$m emission feature due to polycyclic aromatic hydrocarbons peaks between the star and a nearby photodissociation front. The observed features in NGC 7023 are consistent with emission from UV-excited gas-phase C$_{60}$. We find that 0.1-0.6% of interstellar carbon is in C$_{60}$; this abundance is consistent with those from previous upper limits and possible fullerene detections in the interstellar medium. This is the first firm detection of neutral C$_{60}$ in the interstellar medium.

We investigate the density structure of the interstellar medium using new high-resolution maps of the 8 micron, 24 micron, and 70 micron surface brightness towards a molecular cloud in the Gum Nebula, made as part of the Spitzer Space Telescope Galactic First Look Survey. The maps are correlated with 100 micron images measured with IRAS. At 24 and 70 micron, the spatial power spectrum of surface brightness follows a power law with spectral index -3.5. At 24 micron, the power law behavior is remarkably consistent from the 0.2 degree size of our maps down to the 5 arcsecond spatial resolution. Thus, the structure of the 24 micron emission is self-similar even at milliparsec scales. The combined power spectrum produced from Spitzer 24 micron and IRAS 25 micron images is consistent with a change in the power law exponent from -2.6 to -3.5. The decrease may be due to the transition from a two-dimensional to three-dimensional structure. Under this hypothesis, we estimate the thickness of the emitting medium to be 0.3 pc.

The (2P3/2 -> 2P1/2) transition of singly--ionized carbon, [CII], is the primary coolant of diffuse interstellar gas. We describe observations of [CII] emission towards nine high Galactic latitude translucent molecular clouds, made with the long wavelength spectrometer on board the Infrared Space Observatory. To understand the role of dust grains in processing the interstellar radiation field (ISRF) and heating the gas, we compare the [CII] integrated intensity with the far-infrared (FIR) integrated surface brightness for the 101 sampled lines of sight. We find that [CII] is linearly correlated with FIR, and the average ratio is equal to that measured with the COBE satellite for all high-latitude Milky Way gas. There is a significant decrease that was not detected with COBE in [CII] emissivity at high values of FIR. Our sample splits naturally into two populations depending on the 60um/100um surface brightness ratio, or color: ``warm'' positions and ``cold'' positions. A transition from sources with warm to those with cold 60/100 colors coincides approximately with the transition from constant to decreasing [CII] emissivity. We model the [CII] and far-infrared emission under conditions of thermal equilibrium, using the simplifying assumptions that, in all regions heated by the ISRF, the most important source of gas heating is the photoelectric effect on grains and the most important source of gas cooling is [CII] emission. The model matches the data well. There are no statistically significant differences in the derived values of the ISRF intensity and the photoelectric heating efficiency for warm and cold sources. The observed variations in the [CII] emissivity and the 60/100 colors can be understood entirely in terms of the attenuation and softening of the ISRF by translucent clouds, not changes in dust properties.

AST/RO, a 1.7 m diameter telescope for astronomy and aeronomy studies at wavelengths between 200 and 2000 microns, was installed at the South Pole during the 1994-1995 Austral summer. The telescope operates continuously through the Austral winter, and is being used primarily for spectroscopic studies of neutral atomic carbon and carbon monoxide in the interstellar medium of the Milky Way and the Magellanic Clouds. The South Pole environment is unique among observatory sites for unusually low wind speeds, low absolute humidity, and the consistent clarity of the submillimeter sky. Four heterodyne receivers, an array receiver, three acousto-optical spectrometers, and an array spectrometer are installed. A Fabry-Perot spectrometer using a bolometric array and a Terahertz receiver are in development. Telescope pointing, focus, and calibration methods as well as the unique working environment and logistical requirements of the South Pole are described.

The rotational transitions of carbon monoxide (CO) are the primary means of investigating the density and velocity structure of the molecular interstellar medium. Here we study the lowest four rotational transitions of CO towards high-latitude translucent molecular clouds (HLCs). We report new observations of the J = (4-3), (2-1), and (1-0) transitions of CO towards eight high-latitude clouds. The new observations are combined with data from the literature to show that the emission from all observed CO transitions is linearly correlated. This implies that the excitation conditions which lead to emission in these transitions are uniform throughout the clouds. Observed 13CO/12CO (1-0) integrated intensity ratios are generally much greater than the expected abundance ratio of the two species, indicating that the regions which emit 12CO (1-0) radiation are optically thick. We develop a statistical method to compare the observed line ratios with models of CO excitation and radiative transfer. This enables us to determine the most likely portion of the physical parameter space which is compatible with the observations. The model enables us to rule out CO gas temperatures greater than 30K since the most likely high-temperature configurations are 1 pc-sized structures aligned along the line of sight. The most probable solution is a high density and low temperature (HDLT) solution. The CO cell size is approximately 0.01 pc (2000 AU). These cells are thus tiny fragments within the 100 times larger CO-emitting extent of a typical high-latitude cloud. We discuss the physical implications of HDLT cells, and we suggest ways to test for their existence.

We present a simple, semianalytical approach to compute the line emission of C+, C and CO in Photodissociation Regions of varying metallicity that depends on very few parameters and naturally incorporates the clumpiness of the interstellar medium into the calculations. We use it to study the dependence of the I_[C II]/I_CO and I_[C I]/I_CO line ratios on the metallicity. We show how to modify it to include the effects of density and radiation field, and we have compared it with observational data. We find that the model explains the observed trend of enhanced [C II]/CO line ratio with decreasing metallicity as the natural result of the augmented fraction of photodissociated gas in a clump. We show that enhanced [C II]/CO ratios can be produced by lowering the upper limit of the clump size distribution, as may happen in 30 Doradus. We also find that the available data favors a [C I]/CO intensity ratio essentially independent of metallicity, albeit the paucity of observations does not exclude other possibilities. This is difficult to understand if most of the C is produced by UV photons as is the case with C+. Finally, we study the prediction of the model for the trend of the X_CO factor with metallicity. Comparison with previous observational studies yields good agreement.