Far-infrared astronomy has advanced rapidly since its inception in the late 1950's, driven by a maturing technology base and an expanding community of researchers. This advancement has shown that observations at far-infrared wavelengths are important in nearly all areas of astrophysics, from the search for habitable planets and the origin of life, to the earliest stages of galaxy assembly in the first few hundred million years of cosmic history. The combination of a still developing portfolio of technologies, particularly in the field of detectors, and a widening ensemble of platforms within which these technologies can be deployed, means that far-infrared astronomy holds the potential for paradigm-shifting advances over the next decade. In this review, we examine current and future far-infrared observing platforms, including ground-based, sub-orbital, and space-based facilities, and discuss the technology development pathways that will enable and enhance these platforms to best address the challenges facing far-infrared astronomy in the 21st century.

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

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.

Theoretical models and spectroscopic observations of newborn stars suggest that protoplantary disks have an inner "wall" at a distance set by the disk interaction with the star. Around T Tauri stars, the size of this disk hole is expected to be on a 0.1-AU scale that is unresolved by current adaptive optics imaging, though some model-dependent constraints have been obtained by near-infrared interferometry. Here we report the first measurement of the inner disk wall around a solar-mass young stellar object, YLW 16B in the {\rho} Ophiuchi star-forming region, by detecting the light travel time of the variable radiation from the stellar surface to the disk. Consistent time lags were detected on two nights, when the time series in H (1.6 {\mu}m) and K (2.2 {\mu}m) bands were synchronized while the 4.5 {\mu}m emission lagged by 74.5 +/- 3.2 seconds. Considering the nearly edge-on geometry of the disk, the inner rim should be 0.084 AU from the protostar on average, with an error of order 0.01 AU. This size is likely larger than the range of magnetospheric truncations, and consistent with an optically and geometrically thick disk front at the dust sublimation radius at ~1500 K. The widths of the cross-correlation functions between the data in different wavebands place possible new constraints on the geometry of the disk.

We provide CoRoT and Spitzer light curves, as well as broad-band multi-wavelength photometry and high resolution, multi- and single-epoch spectroscopy for 17 classical T Tauris in NGC 2264 whose CoRoT light curves (LCs) exemplify the "stochastic" LC class as defined in Cody et al. (2014). The most probable physical mechanism to explain the optical variability in this LC class is time-dependent mass accretion onto the stellar photosphere, producing transient hot spots. As evidence in favor of this hypothesis, multi-epoch high resolution spectra for a subset of these stars shows that their veiling levels also vary in time and that this veiling variability is consistent in both amplitude and timescale with the optical LC morphology. Furthermore, the veiling variability is well-correlated with the strength of the HeI 6678A emission line, a feature predicted by models to arise in accretion shocks on or near the stellar photosphere. Stars with accretion burst LC morphology (Stauffer et al. 2014) are also attributed to variable mass accretion. Both the stochastic and accretion burst LCs can be explained by a simple model of randomly occurring flux bursts, with the stochastic LC class having a higher frequency of lower amplitude events. Based on their UV excesses, veiling, and mean Ha equivalent widths, members of the stochastic LC class have only moderate time-averaged mass accretion rates. The most common feature of their Ha profiles is for them to exhibit blue-shifted absorption features, most likely originating in a disk wind. The lack of periodic signatures in the LCs suggests that little of the variability is due to long-lived hot spots rotating into or out of our line of sight; instead, the primary driver of the observed photometric variability is likely to be instabilities in the inner disk that lead to variable mass accretion.

To move one step forward toward a Galactic distribution of planets, we present the first planet sensitivity analysis for microlensing events with simultaneous observations from space and the ground. We present this analysis for two such events, OGLE-2014-BLG-0939 and OGLE-2014-BLG-0124, which both show substantial planet sensitivity even though neither of them reached high magnification. This suggests that an ensemble of low to moderate magnification events can also yield significant planet sensitivity and therefore probability to detect planets. The implications of our results to the ongoing and future space-based microlensing experiments to measure the Galactic distribution of planets are discussed.

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.

Space-based microlens parallax measurements are a powerful tool for understanding planet populations, especially their distribution throughout the Galaxy. However, if space-based observations of the microlensing events must be specifically targeted, it is crucial that microlensing events enter the parallax sample without reference to the known presence or absence of planets. Hence, it is vital to define objective criteria for selecting events where possible and to carefully consider and minimize the selection biases where not possible so that the final sample represents a controlled experiment. We present objective criteria for initiating observations and determining their cadence for a subset of events, and we define procedures for isolating subjective decision making from information about detected planets for the remainder of events. We also define procedures to resolve conflicts between subjective and objective selections. These procedures maximize planet sensitivity of the sample as a whole by allowing for planet detections even if they occur before satellite observations for objectively-selected events and by helping to trigger fruitful follow-up observations for subjectively-chosen events. This paper represents our public commitment to these procedures, which is a necessary component of enforcing objectivity on the experimental protocol.

We identify nine young stellar objects (YSOs) in the NGC 2264 star-forming region with optical {\em CoRoT} light curves exhibiting short-duration, shallow, periodic flux dips. All of these stars have infrared (IR) excesses that are consistent with their having inner disk walls near the Keplerian co-rotation radius. The repeating photometric dips have FWHM generally less than one day, depths almost always less than 15%, and periods (3<P<11 days) consistent with dust near the Keplerian co-rotation period. The flux dips vary considerably in their depth from epoch to epoch, but usually persist for several weeks and, in two cases, were present in data collected on successive years. For several of these stars, we also measure the photospheric rotation period and find that the rotation and dip periods are the same, as predicted by standard "disk-locking" models. We attribute these flux dips to clumps of material in or near the inner disk wall, passing through our line of sight to the stellar photosphere. In some cases, these dips are also present in simultaneous {\em Spitzer} IRAC light curves at 3.6 and 4.5 microns. We characterize the properties of these dips, and compare the stars with light curves exhibiting this behavior to other classes of YSO in NGC 2264. A number of physical mechanisms could locally increase the dust scale height near the inner disk wall, and we discuss several of those mechanisms; the most plausible mechanisms are either a disk warp due to interaction with the stellar magnetic field or dust entrained in funnel-flow accretion columns arising near the inner disk wall.

This article summarizes a workshop held on March, 2014, on the potential of the James Webb Space Telescope (JWST) to revolutionize our knowledge of the physical properties of exoplanets through transit observations. JWST's unique combination of high sensitivity and broad wavelength coverage will enable the accurate measurement of transits with high signal-to-noise. Most importantly, JWST spectroscopy will investigate planetary atmospheres to determine atomic and molecular compositions, to probe vertical and horizontal structure, and to follow dynamical evolution, i.e. exoplanet weather. JWST will sample a diverse population of planets of varying masses and densities in a wide variety of environments characterized by a range of host star masses and metallicities, orbital semi-major axes and eccentricities. A broad program of exoplanet science could use a substantial fraction of the overall JWST mission.

We present the Coordinated Synoptic Investigation of NGC 2264, a continuous 30-day multi-wavelength photometric monitoring campaign on more than 1000 young cluster members using 16 telescopes. The unprecedented combination of multi-wavelength, high-precision, high-cadence, and long-duration data opens a new window into the time domain behavior of young stellar objects. Here we provide an overview of the observations, focusing on results from Spitzer and CoRoT. The highlight of this work is detailed analysis of 162 classical T Tauri stars for which we can probe optical and mid-infrared flux variations to 1% amplitudes and sub-hour timescales. We present a morphological variability census and then use metrics of periodicity, stochasticity, and symmetry to statistically separate the light curves into seven distinct classes, which we suggest represent different physical processes and geometric effects. We provide distributions of the characteristic timescales and amplitudes, and assess the fractional representation within each class. The largest category (>20%) are optical "dippers" having discrete fading events lasting ~1-5 days. The degree of correlation between the optical and infrared light curves is positive but weak; notably, the independently assigned optical and infrared morphology classes tend to be different for the same object. Assessment of flux variation behavior with respect to (circum)stellar properties reveals correlations of variability parameters with H$\alpha$ emission and with effective temperature. Overall, our results point to multiple origins of young star variability, including circumstellar obscuration events, hot spots on the star and/or disk, accretion bursts, and rapid structural changes in the inner disk.

We present Submillimeter Array (SMA) $\lambda =$ 0.88 and 1.3 mm broad band observations, and the Jansky Very Large Array (VLA) observations in $\rm{NH_3}$ $(J,K) = (1,1)$ up to $(5,5)$, $\rm{H_2O}$ and $\rm{CH_3OH}$ maser lines toward the two most massive molecular clumps in infrared dark cloud (IRDC) G11.11-0.12. Sensitive high-resolution images reveal hierarchical fragmentation in dense molecular gas from the $\sim 1$ pc clump scale down to $\sim 0.01$ pc condensation scale. At each scale, the mass of the fragments is orders of magnitude larger than the Jeans mass. This is common to all four IRDC clumps we studied, suggesting that turbulence plays an important role in the early stages of clustered star formation. Masers, shock heated $\rm{NH_3}$ gas, and outflows indicate intense ongoing star formation in some cores while no such signatures are found in others. Furthermore, chemical differentiation may reflect the difference in evolutionary stages among these star formation seeds. We find $\rm{NH_3}$ ortho/para ratios of $1.1\pm0.4$, $2.0\pm0.4$, and $3.0\pm0.7$ associated with three outflows, and the ratio tends to increase along the outflows downstream. Our combined SMA and VLA observations of several IRDC clumps present the most in depth view so far of the early stages prior to the hot core phase, revealing snapshots of physical and chemical properties at various stages along an apparent evolutionary sequence.

Based on more than four weeks of continuous high cadence photometric monitoring of several hundred members of the young cluster NGC 2264 with two space telescopes, NASA's Spitzer and the CNES CoRoT (Convection, Rotation, and planetary Transits), we provide high quality, multi-wavelength light curves for young stellar objects (YSOs) whose optical variability is dominated by short duration flux bursts, which we infer are due to enhanced mass accretion rates. These light curves show many brief -- several hour to one day -- brightenings at optical and near-infrared (IR) wavelengths with amplitudes generally in the range 5-50% of the quiescent value. Typically, a dozen or more of these bursts occur in a thirty day period. We demonstrate that stars exhibiting this type of variability have large ultraviolet (UV) excesses and dominate the portion of the u-g vs. g-r color-color diagram with the largest UV excesses. These stars also have large Halpha equivalent widths, and either centrally peaked, lumpy Halpha emission profiles or profiles with blue-shifted absorption dips associated with disk or stellar winds. Light curves of this type have been predicted for stars whose accretion is dominated by Rayleigh-Taylor instabilities at the boundary between their magnetosphere and inner circumstellar disk, or where magneto-rotational instabilities modulate the accretion rate from the inner disk. Amongst the stars with the largest UV excesses or largest Halpha equivalent widths, light curves with this type of variability greatly outnumber light curves with relatively smooth sinusoidal variations associated with long-lived hot spots. We provide quantitative statistics for the average duration and strength of the accretion bursts and for the fraction of the accretion luminosity associated with these bursts.

This paper presents a study of the rate and efficiency of star formation in the NGC 6334 star forming region. We obtained observations at J, H, and Ks taken with the NOAO Extremely Wide-Field Infrared Imager (NEWFIRM) and combined them with observations taken with the Infrared Array Camera (IRAC) camera on the Spitzer Space Telescope at wavelengths {\lambda} = 3.6, 4.5, 5.8, and 8.0 {\mu}m. We also analyzed previous observations taken at 24 {\mu}m using the Spitzer MIPS camera as part of the MIPSGAL survey. We have produced a point source catalog with >700,000 entries. We have identified 2283 YSO candidates, 375 Class I YSOs and 1908 Class II YSOs using a combination of existing IRAC-based color classification schemes that we have extended and validated to the near-IR for use with warm Spitzer data. We have identified multiple new sites of ongoing star formation activity along filamentary structures extending tens of parsecs beyond the central molecular ridge of NGC 6334. By mapping the extinction we derived an estimate for the gas mass, 2.2 x 10^5 M_sun. The heavy concentration of protostars along the dense filamentary structures indicates that NGC 6334 may be undergoing a "mini-starburst" event with {\Sigma}{SFR}>8.2 M_sun Myr^-1 pc^-2 and SFE>0.10. We have used these estimates to place NGC 6334 in the Kennicutt-Schmidt diagram to help bridge the gap between observations of local low-mass star forming regions and star formation in other galaxies.

Multi-epoch Spitzer Space Telescope 24 micron data is utilized from the MIPSGAL and Taurus Legacy surveys to detect asteroids based on their relative motion. These infrared detections are matched to known asteroids and rotationally averaged diameters and albedos are derived using the Near Earth Asteroid Model (NEATM) in conjunction with Monte Carlo simulations for 1835 asteroids ranging in size from 0.2 to 143.6 km. A small subsample of these objects was also detected by IRAS or MSX and the single wavelength albedo and diameter fits derived from this data are within 5% of the IRAS and/or MSX derived albedos and diameters demonstrating the robustness of our technique. The mean geometric albedo of the small main belt asteroids in this sample is p_V = 0.138 with a sample standard deviation of 0.105. The albedo distribution of this sample is far more diverse than the IRAS or MSX samples. The cumulative size-frequency distribution of asteroids in the main belt at small diameters is directly derived. Completeness limits of the optical and infrared surveys are discussed.

We investigate the environment of the very low-luminosity object L1521F IRS using data from the Taurus Spitzer Legacy Survey. The MIPS 160 micron image shows both extended emission from the Taurus cloud as well as emission from multiple cold cores over a 1 X 2 deg region. Analysis shows that the cloud dust temperature is 14.2 +- 0.4 K and the extinction ratio is A_160/A_K = 0.010 +- 0.001 up to A_V ~ 4 mag. We find kappa_160 = 0.23 +- 0.046 cm^2/g for the specific opacity of the gas-dust mixture. Therefore, for dust in the Taurus cloud we find the 160 um opacity is significantly higher than that measured for the diffuse ISM, but not too different from dense cores, even at modest extinction values. Furthermore, the 160 um image shows features that do not appear in the IRAS 100 um image. We identify six regions as cold cores, i.e. colder than 14.2 K, all of which have counterparts in extinction maps or C18O maps. We compare the effects of L1521F IRS on its natal core and find there is no evidence for dust heating at 160 or 100 um by the embedded source. From the infrared luminosity L_TIR = 0.024 Lo we find L_bol = 0.034 - 0.046L_o, thus confirming the source's low-luminosity. Comparison of L1521F IRS with theoretical simulations for the very early phases of star formation appears to rule out the first core collapse phase. The evolutionary state appears similar to or younger than the class 0 phase, and the estimated mass is likely to be substellar.

We present several corrections for point source photometry to be applied to data from the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. These corrections are necessary because of characteristics of the IRAC arrays and optics and the way the instrument is calibrated in-flight. When these corrections are applied, it is possible to achieve a ~2% relative photometric accuracy for sources of adequate signal to noise in an IRAC image.

Water is predicted to be among, if not the most abundant molecular species after hydrogen in the atmospheres of close-in extrasolar giant planets (hot-Jupiters) Several attempts have been made to detect water on an exoplanet, but have failed to find compelling evidence for it or led to claims that should be taken with caution. Here we report an analysis of recent observations of the hot-Jupiter HD189733b taken during the transit, where the planet passed in front of its parent star. We find that absorption by water vapour is the most likely cause of the wavelength-dependent variations in the effective radius of the planet at the infrared wavelengths 3.6, 5.8 and 8 microns. The larger effective radius observed at visible wavelengths may be due to either star variability or the presence of clouds/hazes. We explain the most recent thermal infrared observations of the planet during secondary transit behind the star, reporting a non-detection of water on HD189733b, as being a consequence of the nearly isothermal vertical profile of the planet.s atmosphere. Our results show that water is detectable on extrasolar planets using the primary transit technique and that the infrared should be a better wavelength region than the visible, for such searches.

The Infrared Array Camera (IRAC) on the Spitzer Space Telescope is absolutely calibrated by comparing photometry on a set of A stars near the north ecliptic pole to predictions based on ground-based observations and a stellar atmosphere model. The brightness of point sources is calibrated to an accuracy of 3%, relative to models for A star stellar atmospheres, for observations performed and analyzed in the same manner as the calibration stars. This includes corrections for location of the star in the array and the location if the centroid within the peak pixel. Long-term stability of the IRAC photometry was measured by monitoring the brightness of A dwarfs and K giants (near the north ecliptic pole) observed several times per month; the photometry is stable to 1.5% (rms) over a year. Intermediate-time-scale stability of the IRAC photometry was measured by monitoring at least one secondary calibrator (near the ecliptic plane) every 12 hr while IRAC is in nominal operations; the intermediate-term photometry is stable with a 1% dispersion (rms). One of the secondary calibrators was found to have significantly time-variable (5%) mid-infrared emission, with period (7.4 days) matching the optical light curve; it is possibly a Cepheid variable.

The optically-dark globule IC 1396A is revealed using Spitzer images at 3.6, 4.5, 5.8, 8, and 24 microns to be infrared-bright and to contain a set of previously unknown protostars. The mid-infrared colors of the 24 microns detected sources indicate several very young (Class I or 0) protostars and a dozen Class II stars. Three of the new sources (IC 1396A: gamma, delta, and epsilon) emit over 90% of their bolometric luminosities at wavelengths greater than 3 microns, and they are located within ~0.02 pc of the ionization front at the edge of the globule. Many of the sources have spectra that are still rising at 24 microns. The two previously-known young stars LkHa 349 a and c are both detected, with component c harboring a massive disk and component a being bare. Of order 5% of the mass of material in the globule is presently in the form of protostars in the 10^5 to 10^6 yr age range. This high star formation rate was likely triggered by radiation from a nearby O star.