Gaussian processes (GPs) are commonly used as a model of stochastic variability in astrophysical time series. In particular, GPs are frequently employed to account for correlated stellar variability in planetary transit light curves. The efficient application of GPs to light curves containing thousands to tens of thousands of datapoints has been made possible by recent advances in GP methods, including the celerite method. Here we present an extension of the celerite method to two input dimensions, where, typically, the second dimension is small. This method scales linearly with the total number of datapoints when the noise in each large dimension is proportional to the same celerite kernel and only the amplitude of the correlated noise varies in the second dimension. We demonstrate the application of this method to the problem of measuring precise transit parameters from multiwavelength light curves and show that it has the potential to improve transit parameters measurements by orders of magnitude. Applications of this method include transit spectroscopy and exomoon detection, as well a broader set of astronomical problems.

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 the detection and analysis of the phase curve of KELT-1b at optical wavelengths, analyzing data taken by the Transiting Exoplanet Survey Satellite (TESS). With a mass of ~27 M_J, KELT-1b is a low-mass brown dwarf. Due to the high mass and close proximity of its companion, the host star has a TESS light curve which shows clear ellipsoidal variations. We model the data with a six-component model: secondary eclipse, phase curve accounting for reflected light and thermal emission, Doppler beaming, ellipsoidal variations, stellar activity and the primary transit. We determine the secondary eclipse depth in the TESS bandpass to be 304 +/- 75 parts-per-million (ppm), the most accurate eclipse depth determined so far for KELT-1b. We measure the amplitude of the phase curve to be 128 +/- 27 ppm, with a corresponding eastward offset between the region of maximum brightness and the substellar point of 19.2 +/- 9.6 degrees, in good agreement with Spitzer measurements. We determine day and night brightness temperatures of 3201 +/- 147 K and 1484 +/- 110 K, respectively, slightly higher than those from Spitzer 3.6 and 4.5 micrometer data. A one-dimensional self-consistent atmospheric model can explain the TESS and Spitzer day side brightness temperatures with thermal emission alone and no reflected light. The night side data can be explained by a model with an internal temperature of ~1100 K, which may be related to the inflated radius. The difference between the TESS and Spitzer brightness temperatures can be explained by stronger molecular opacity in the Spitzer bands. On the night side, this opacity is due primarily to CH4 and CO while on the day side it is due to H2-H2 and H2-He collision-induced absorption.

We present the discovery of Kepler-88 d ($P_d = 1403\pm14$ days, $M\mathrm{sin}i_d = 965\pm44\,M_\oplus = 3.04\pm0.14\,M_J$, $e_d = 0.41\pm0.03$) based on six years of radial velocity (RV) follow-up from the W. M. Keck Observatory HIRES spectrograph. Kepler-88 has two previously identified planets. Kepler-88 b (KOI-142.01) transits in the NASA \Kepler\ photometry and has very large transit timing variations. \citet{Nesvorny2013} perfomed a dynamical analysis of the TTVs to uniquely identify the orbital period and mass of the perturbing planet (Kepler-88 c), which was later was confirmed with RVs from the Observatoire de Haute-Provence (OHP, Barros et al. 2014). To fully explore the architecture of this system, we performed photodynamical modeling on the \Kepler\ photometry combined with the RVs from Keck and OHP and stellar parameters from spectroscopy and Gaia. Planet d is not detectable in the photometry, and long-baseline RVs are needed to ascertain its presence. A photodynamical model simultaneously optimized to fit the RVs and \Kepler\ photometry yields the most precise planet masses and orbital properties yet for b and c: $P_b = 10.91647\pm0.00014\,\mathrm{days}$, $M_b=9.5\pm1.2\,M_\oplus$, $P_c=22.2649\pm0.0007\,\mathrm{days}$, and $M_c=214.1\pm5.3\,M_\oplus$. The photodynamical solution also finds that planets b and c have low eccentricites and low mutual inclination, are apsidally anti-aligned, and have conjunctions on the same hemisphere of the star. Continued RV follow-up of systems with small planets will improve our understanding of the link between inner planetary system architectures and giant planets.

One of the main science motivations for the ESA PLAnetary Transit and Oscillations (PLATO) mission is to measure exoplanet transit radii with 3% precision. In addition to flares and starspots, stellar oscillations and granulation will enforce fundamental noise floors for transiting exoplanet radius measurements. We simulate light curves of Earth-sized exoplanets transiting continuum intensity images of the Sun taken by the HMI instrument aboard SDO to investigate the uncertainties introduced on the exoplanet radius measurements by stellar granulation and oscillations. After modeling the solar variability with a Gaussian process, we find that the amplitude of solar oscillations and granulation is of order 100 ppm -- similar to the depth of an Earth transit -- and introduces a fractional uncertainty on the depth of transit of 0.73% assuming four transits are observed over the mission duration. However, when we translate the depth measurement into a radius measurement of the planet, we find a much larger radius uncertainty of 3.6%. This is due to a degeneracy between the transit radius ratio, the limb-darkening, and the impact parameter caused by the inability to constrain the transit impact parameter in the presence of stellar variability. We find that surface brightness inhomogeneity due to photospheric granulation contributes a lower limit of only 2 ppm to the photometry in-transit. The radius uncertainty due to granulation and oscillations, combined with the degeneracy with the transit impact parameter, accounts for a significant fraction of the error budget of the PLATO mission, before detector or observational noise is introduced to the light curve. If it is possible to constrain the impact parameter or to obtain follow-up observations at longer wavelengths where limb-darkening is less significant, this may enable higher precision radius measurements.

The upcoming launch of the James Webb Space Telescope (JWST) combined with the unique features of the TRAPPIST-1 planetary system should enable the young field of exoplanetology to enter into the realm of temperate Earth-sized worlds. Indeed, the proximity of the system (12pc) and the small size (0.12 Rsun) and luminosity (0.05 Lsun) of its host star should make the comparative atmospheric characterization of its seven transiting planets within reach of an ambitious JWST program. Given the limited lifetime of JWST, the ecliptic location of the star that limits its visibility to 100d per year, the large number of observational time required by this study, and the numerous observational and theoretical challenges awaiting it, its full success will critically depend on a large level of coordination between the involved teams and on the support of a large community. In this context, we present here a community initiative aiming to develop a well-defined sequential structure for the study of the system with JWST and to coordinate on every aspect of its preparation and implementation, both on the observational (e.g. study of the instrumental limitations, data analysis techniques, complementary space-based and ground-based observations) and theoretical levels (e.g. model developments and comparison, retrieval techniques, inferences). Depending on the outcome of the first phase of JWST observations of the planets, this initiative could become the seed of a major JWST Legacy Program devoted to the study of TRAPPIST-1.

Transit Timing Variations, or TTVs, can be a very efficient way of constraining masses and eccentricities of multi-planet systems. Recent measurements of the TTVs of TRAPPIST-1 led to an estimate of the masses of the planets, enabling an estimate of their densities. A recent TTV analysis using data obtained in the past two years yields a 34% and 13% increase in mass for TRAPPIST-1b and c, respectively. In most studies to date, a Newtonian N-body model is used to fit the masses of the planets, while sometimes general relativity is accounted for. Using the Posidonius N-body code, we show that in the case of the TRAPPIST-1 system, non-Newtonian effects might be also relevant to correctly model the dynamics of the system and the resulting TTVs. In particular, using standard values of the tidal Love number $k_2$ (accounting for the tidal deformation) and the fluid Love number $k_{2f}$ (accounting for the rotational flattening) leads to differences in the TTVs of TRAPPIST-1b and c similar to the differences caused by general relativity. We also show that relaxing the values of tidal Love number $k_2$ and the fluid Love number $k_{2f}$ can lead to TTVs which differ by as much as a few 10~s on a $3-4$-year timescale, which is a potentially observable level. The high values of the Love numbers needed to reach observable levels for the TTVs could be achieved for planets with a liquid ocean, which, if detected, might then be interpreted as a sign that TRAPPIST-1b and TRAPPIST-1c could have a liquid magma ocean. For TRAPPIST-1 and similar systems, the models to fit the TTVs should potentially account for general relativity, for the tidal deformation of the planets, for the rotational deformation of the planets and, to a lesser extent, for the rotational deformation of the star, which would add up to 7x2+1 = 15 additional free parameters in the case of TRAPPIST-1.

We derive analytic, closed-form solutions for the light curve of a planet transiting a star with a limb darkening profile which is a polynomial function of the stellar elevation, up to arbitrary integer order. We provide improved analytic expressions for the uniform, linear, and quadratic limb-darkened cases, as well as novel expressions for higher order integer powers of limb darkening. The formulae are crafted to be numerically stable over the expected range of usage. We additionally present analytic formulae for the partial derivatives of instantaneous flux with respect to the radius ratio, impact parameter, and limb darkening coefficients. These expressions are rapid to evaluate, and compare quite favorably in speed and accuracy to existing transit light curve codes. We also use these expressions to numerically compute the first partial derivatives of exposure-time averaged transit light curves with respect to all model parameters. An additional application is modeling eclipsing binary or eclipsing multiple star systems in cases where the stars may be treated as spherically symmetric. We provide code which implements these formulae in C++, Python, IDL, and Julia, with tests and examples of usage.

To date a dozen transiting "Tatooines" or circumbinary planets (CBPs) have been discovered, by eye, in the data from the Kepler mission; by contrast, thousands of confirmed circumstellar planets orbiting around single stars have been detected using automated algorithms. Automated detection of CBPs is challenging because their transits are strongly aperiodic with irregular profiles. Here, we describe an efficient and automated technique for detecting circumbinary planets that transit their binary hosts in Kepler light curves. Our method accounts for large transit timing and duration variations (TTVs and TDVs), induced by binary reflex motion, in two ways: 1) We directly correct for large-scale TTVs and TDVs in the light curves by using Keplerian models to approximate binary and CBP orbits; and 2) We allow additional aperiodicities on the corrected light curves by employing the Quasi-periodic Automated Transit Search algorithm (QATS). We demonstrate that our method dramatically improves detection significance using simulated data and two previously identified CBP systems, Kepler-35 and Kepler-64.

We implement a search for exoplanets in campaigns zero through eight (C0-8) of the K2 extension of the Kepler spacecraft. We apply a modified version of the QATS planet search algorithm to K2 light curves produced by the EVEREST pipeline, carrying out the C0-8 search on $1.5 \times 10^5$ target stars with magnitudes in the range of Kp = 9-15. We detect 818 transiting planet candidates, of which 374 were undiscovered by prior searches, with {64,15,5,2,1} in {2,3,4,5,6}-planet multi-planet candidate systems, respectively. Of the new planets detected, 100 orbit M dwarfs, including one which is potentially rocky and in the habitable zone. 154 of our candidates reciprocally transit with our Solar System: they are geometrically aligned to see at least one Solar System planet transit. We find candidates which display transit timing variations and dozens of candidates on both period extremes with single transits or ultra-short periods. We point to evidence that our candidates display similar patterns in frequency and size-period relation as confirmed planets, such as tentative evidence for the radius gap. Confirmation of these planet candidates with follow-up studies will increase the number of K2 planets by up to 50%, and characterization of their host stars will improve statistical studies of planet properties. Our sample includes many planets orbiting bright stars amenable for radial velocity follow-up and future characterization with JWST. We also list the 579 eclipsing binary systems detected as part of this search.

We report on the properties of eclipsing binaries from the Kepler mission with a newly developed photometric modeling code, which uses the light curve, spectral energy distribution of each binary, and stellar evolution models to infer stellar masses without the need for radial velocity measurements. We present solutions and posteriors to orbital and stellar parameters for 728 systems, forming the largest homogeneous catalogue of full Kepler binary parameter estimates to date. Using comparisons to published radial velocity measurements, we demonstrate that the inferred properties (e.g., masses) are reliable for well-detached main-sequence binaries, which make up the majority of our sample. The fidelity of our inferred parameters degrades for a subset of systems not well described by input isochrones, such as short-period binaries that have undergone interactions, or binaries with post-main sequence components. Additionally, we identify 35 new systems which show evidence of eclipse timing variations, perhaps from apsidal motion due to binary tides or tertiary companions. We plan to subsequently use these models to search for and constrain the presence of circumbinary planets in Kepler eclipsing binary systems.

K2-146 is a mid-M dwarf ($M_\star = 0.331 \pm 0.009 M_\odot$; $R_\star = 0.330 \pm 0.010 R_\odot$), observed in Campaigns 5, 16, and 18 of the K2 mission. In Campaign 5 data, a single planet was discovered with an orbital period of $2.6$~days and large transit timing variations due to an unknown perturber. Here we analyze data from Campaigns 16 and 18, detecting the transits of a second planet, c, with an orbital period of $4.0$~days, librating in a 3:2 resonance with planet b. Large, anti-correlated timing variations of both planets exist due to their resonant perturbations. The planets have a mutual inclination of $2.40^\circ\pm0.25^\circ$, which torqued planet c more closely into our line-of-sight. Planet c was grazing in Campaign 5 and thus missed in previous searches; in Campaigns 16 and 18 it is fully transiting, and its transit depth is three times larger. We improve the stellar properties using data from Gaia DR2, and using dynamical fits find that both planets are sub-Neptunes: their masses are $5.77\pm0.18$ and $7.50\pm0.23 M_{\oplus}$ and their radii are $2.04\pm0.06$ and $2.19\pm0.07$ R$_\oplus$, respectively. These mass constraints set the precision record for small exoplanets (a few gas giants have comparable relative precision). These planets lie in the photoevaporation valley when viewed in Radius-Period space, but due to the low-luminosity M-dwarf host star, they lie among the atmosphere-bearing planets when viewed in Radius-Irradiation space. This, along with their densities being 60%-80% that of Earth, suggests that they may both have retained a substantial gaseous envelope.

We present the next generation public exoplanet fitting software, EXOFASTv2. It is capable of fitting an arbitrary number of planets, radial velocity data sets, astrometric data sets, and/or transits observed with any combination of wavelengths. We model the star simultaneously in the fit and provide several state-of-the-art ways to constrain its properties, including taking advantage of the now-ubiquitous all-sky catalog photometry and Gaia parallaxes. EXOFASTv2 can model the star by itself, too. Multi-planet systems are modeled self-consistently with the same underlying stellar mass that defines their semi-major axes through Kepler's law and the planetary period. Transit timing, duration, and depth variations can be modeled with a simple command line option. We explain our methodology and rationale as well as provide an improved version of the core transit model that is both 25\% faster and more accurate. We highlight several potential pitfalls in exoplanet modeling, including the handling of eccentricity in transit-only fits, that the standard exoplanet convention for $\omega$ uses a left-handed coordinate system, contrary to most modern textbooks, how to avoid an important degeneracy when allowing negative companion masses, and a widely unappreciated, potential 10-minute ambiguity in the reported transit times. EXOFASTv2 is available at https://github.com/jdeast/EXOFASTv2 . The code is written in IDL, and includes an executable that can be run freely and legally without an IDL license or any knowledge of the language. Extensive documentation and tutorials are included in the distribution for a variety of example fits. Advanced amateurs and undergrads have successfully performed sophisticated global fits of complex planetary systems with EXOFASTv2. It is therefore a powerful tool for education and outreach as well as the broader professional community.

Context. The Kepler Object of Interest Network (KOINet) is a multi-site network of telescopes around the globe organised for follow-up observations of transiting planet candidate Kepler objects of interest (KOIs) with large transit timing variations (TTVs). The main goal of KOINet is the completion of their TTV curves as the Kepler telescope stopped observing the original Kepler field in 2013. Aims. We ensure a comprehensive characterisation of the investigated systems by analysing Kepler data combined with new ground-based transit data using a photodynamical model. This method is applied to the Kepler-82 system leading to its first dynamic analysis. Methods. In order to provide a coherent description of all observations simultaneously, we combine the numerical integration of the gravitational dynamics of a system over the time span of observations with a transit light curve model. To explore the model parameter space, this photodynamical model is coupled with a Markov chain Monte Carlo algorithm. Results. The Kepler-82b/c system shows sinusoidal TTVs due to their near 2:1 resonance dynamical interaction. An additional chopping effect in the TTVs of Kepler-82c hints to a further planet near the 3:2 or 3:1 resonance. We photodynamically analysed Kepler long- and short-cadence data and three new transit observations obtained by KOINet between 2014 and 2018. Our result reveals a non-transiting outer planet with a mass of $m_f=20.9\pm1.0\;M_\bigoplus$ near the 3:2 resonance to the outermost known planet, Kepler-82c. Furthermore, we determined the densities of planets b and c to the significantly more precise values $\rho_b=0.98_{-0.14}^{+0.10}\;\text{g cm}^{-3}$ and $\rho_c=0.494_{-0.077}^{+0.066}\;\text{g cm}^{-3}$.

We gathered high resolution spectra for an ensemble of 55 bright active and inactive stars using the ARC 3.5 m Telescope Echelle Spectrograph at Apache Point Observatory ($R\approx$31,500). We measured spectroscopic effective temperatures, surface gravities and metallicities for most stars in the sample with SME and MOOG. Our stellar property results are consistent with the photometric effective temperatures from the Gaia DR2 pipeline. We also measured their chromospheric $S$ and $\log R^\prime_\mathrm{HK}$ indices to classify the stars as active or inactive and study the connection between chromospheric activity and starspots. We then attempted to infer the starspot covering fractions on the active stars by modeling their spectra as a linear combination of hot and cool inactive stellar spectral templates. We find that it is critical to use precise colors of the stars to place stringent priors on the plausible spot covering fractions. The inferred spot covering fractions generally increase with the chromospheric activity indicator $\log R^\prime_\mathrm{HK}$, though we are largely insensitive to spot coverages smaller than $f_S \lesssim 20$%. We find a dearth of stars with small $\log R^\prime_\mathrm{HK}$ and significant spot coverages.

Spectroscopic Eclipsing Binaries (SEBs) are fundamental benchmarks in stellar astrophysics and today are observed in breathtaking detail by missions like TESS, Kepler, and APOGEE. We develop a methodology for simultaneous analysis of high precision Kepler light curves and high resolution near-IR spectra from APOGEE and present orbital solutions and evolutionary histories for a subset of SEBs within this overlap. Radial velocities extracted from APOGEE spectra using the Broadening Function technique are combined with Kepler light curves and to yield binary orbital solutions. The Broadening Function approach yields more precise radial velocities than the standard Cross-Correlation Function, which in turn yields more precise orbital parameters and enables the identification of tertiary stars. The orbital periods of these seven SEBs range from 4 to 40 days. Four of the systems (KIC 5285607, KIC 6864859, KIC 6778289, and KIC 4285087) are well-detached binaries. The remaining three systems have apparent tertiary companions, but each exhibits two eclipses along with at least one spectroscopically varying component (KIC 6449358, KIC 6131659, and KIC 6781535). Gaia distances are available for four targets which we use to estimate temperatures of both members of these SEBs. We explore evolutionary histories in H-R diagram space and estimate ages for this subset of our sample. Finally, we consider the implications for the formation pathways of close binary systems via interactions with tertiary companions. Our methodology combined with the era of big data and observation overlap opens up the possibility of discovering and analyzing large numbers of diverse SEBs, including those with high flux ratios and those in triple systems.

KOI-3278 is a self-lensing stellar binary consisting of a white-dwarf secondary orbiting a Sun-like primary star. Kruse and Agol (2014) noticed small periodic brightenings every 88.18 days in the Kepler photometry and interpreted these as the result of microlensing by a white dwarf with about 63$\%$ of the mass of the Sun. We obtained two sets of spectra for the primary that allowed us to derive three sets of spectroscopic estimates for its effective temperature, surface gravity, and metallicity for the first time. We used these values to update the Kruse and Agol (2014) Einsteinian microlensing model, resulting in a revised mass for the white dwarf of $0.539^{+0.022}_{-0.020} \, M_{\odot}$. The spectra also allowed us to determine radial velocities and derive orbital solutions, with good agreement between the two independent data sets. An independent Newtonian dynamical MCMC model of the combined velocities yielded a mass for the white dwarf of $0.5122^{+0.0057}_{-0.0058} \, M_{\odot}$. The nominal uncertainty for the Newtonian mass is about four times better than for the Einsteinian, $\pm 1.1\%$ vs. $\pm 4.1\%$ and the difference between the two mass determinations is $5.2 \%$. We then present a joint Einsteinian microlensing and Newtonian radial velocity model for KOI-3278, which yielded a mass for the white dwarf of $0.5250^{+0.0082}_{-0.0089} \, M_{\odot}$. This joint model does not rely on any white dwarf evolutionary models or assumptions on the white dwarf mass-radius relation. We discuss the benefits of a joint model of self-lensing binaries, and how future studies of these systems can provide insight into the mass-radius relation of white dwarfs.

Kepler-62f is the first exoplanet small enough to plausibly have a rocky composition orbiting within the habitable zone (HZ) discovered by the Kepler Mission. The planet is 1.4 times the size of the Earth and has an orbital period of 267 days. At the time of its discovery, it had the longest period of any small planet in the habitable zone of a multi-planet system. Because of its long period, only four transits were observed during Kepler's interval of observations. It was initially missed by the Kepler pipeline, but the first three transits were identified by an independent search by Eric Agol, and it was identified as a planet candidate in subsequent Kepler catalogs. However in the latest catalog of exoplanets (Thompson et al., 2018), it is labeled as a false positive. Recent exoplanet catalogues have evolved from subjective classification to automatic classifications of planet candidates by algorithms (such as `Robovetter'). While exceptionally useful for producing a uniform catalogue, these algorithms sometimes misclassify planet candidates as a false positive, as is the case of Kepler-62f. In particularly valuable cases, i.e., when a small planet has been found orbiting in the habitable zone (HZ), it is important to conduct comprehensive analyses of the data and classification protocols to provide the best estimate of the true status of the detection. In this paper we conduct such analyses and show that Kepler-62f is a true planet and not a false positive. The table of stellar and planet properties has been updated based on GAIA results.

We describe the circumstances that led to the discovery of Kepler-36b, and the subsequent characterization of its host planetary system. The Kepler-36 system is remarkable for its physical properties: the close separation of the planets, the contrasting densities of the planets despite their proximity, and the short chaotic timescale. Its discovery and characterization was also remarkable for the novelty of the detection technique and for the precise characterization due to the large transit-timing variations caused by the close proximity of the planets, as well as the precise stellar parameters due to asteroseismology. This was the first multi-planet system whose transit data was processed using a fully consistent photometric-dynamical model, using population Markov Chain Monte Carlo techniques to precisely constrain system parameters. Amongst those parameters, the stellar density was found to be consistent with a complementary, concurrent asteroseismic analysis. In a first, the 3D orientation of the planets was constrained from the lack of transit-duration variations. The system yielded insights into the composition and evolution of short-period planet systems. The denser planet appears to have an Earth-like composition, with uncertainties comparable to the highest precision rocky exoplanet measurements, and the planet densities foreshadowed the rocky/gaseous boundary. The formation of this system remains a mystery, but should yield insights into the migration and evolution of compact exoplanet systems.

Of the nine confirmed transiting circumbinary planet systems, only Kepler-47 is known to contain more than one planet. Kepler-47 b (the "inner planet") has an orbital period of 49.5 days and a radius of about $3\,R_{\oplus}$. Kepler-47 c (the "outer planet") has an orbital period of 303.2 days and a radius of about $4.7\,R_{\oplus}$. Here we report the discovery of a third planet, Kepler-47 d (the "middle planet"), which has an orbital period of 187.4 days and a radius of about $7\,R_{\oplus}$. The presence of the middle planet allows us to place much better constraints on the masses of all three planets, where the $1\sigma$ ranges are less than $26\,M_{\oplus}$, between $7-43\,M_{\oplus}$, and between $2-5\,M_{\oplus}$ for the inner, middle, and outer planets, respectively. The middle and outer planets have low bulk densities, with $\rho_{\rm middle} < 0.68$ g cm$^{-3}$ and $\rho_{\rm outer} < 0.26$ g cm$^{-3}$ at the $1\sigma$ level. The two outer planets are "tightly packed," assuming the nominal masses, meaning no other planet could stably orbit between them. All of the orbits have low eccentricities and are nearly coplanar, disfavoring violent scattering scenarios and suggesting gentle migration in the protoplanetary disk.

We report the Transiting Exoplanet Survey Satellite ($TESS$) detection of a multi-planet system orbiting the $V=10.9$ K0 dwarf TOI 125. We find evidence for up to five planets, with varying confidence. Three high signal-to-noise transit signals correspond to sub-Neptune-sized planets ($2.76$, $2.79$, and $2.94\ R_{\oplus}$), and we statistically validate the planetary nature of the two inner planets ($P_b = 4.65$ days, $P_c = 9.15$ days). With only two transits observed, we report the outer object ($P_{.03} = 19.98$ days) as a high signal-to-noise ratio planet candidate. We also detect a candidate transiting super-Earth ($1.4\ R_{\oplus}$) with an orbital period of only $12.7$ hours and a candidate Neptune-sized planet ($4.2\ R_{\oplus}$) with a period of $13.28$ days, both at low signal-to-noise. This system is amenable to mass determination via radial velocities and transit timing variations, and provides an opportunity to study planets of similar size while controlling for age and environment. The ratio of orbital periods between TOI 125 b and c ($P_c/P_b = 1.97$) is slightly smaller than an exact 2:1 commensurability and is atypical of multiple planet systems from $Kepler$, which show a preference for period ratios just $wide$ of first-order period ratios. A dynamical analysis refines the allowed parameter space through stability arguments and suggests that, despite the nearly commensurate periods, the system is unlikely to be in resonance.

We use archival daily spot coverage measurements from Howard et al. (1984) to study the rotational modulation of the Sun as though it were a distant star. A quasi-periodic Gaussian process measures the solar rotation period $P_\mathrm{rot} = 26.3 \pm 0.1$ days, and activity cycle period $P_\mathrm{cyc} = 10.7 \pm 0.3$ years. We attempt to search for evidence of differential rotation in variations of the apparent rotation period throughout the activity cycle and do not detect a clear signal of differential rotation, consistent with the null results of the hare-and-hounds exercise of Aigrain et al. (2015). The full reconstructed solar light curve is available online.

Presently, the Galactic plane receives relatively few observations compared to most of the LSST footprint. While this may address static science, the plane will also represent the highest density of variable Galactic sources. The proper characterization of variability of these sources will benefit greatly from observations at a higher cadence.

LHS 1140 is a nearby mid-M dwarf known to host a temperate rocky super-Earth (LHS 1140 b) on a 24.737-day orbit. Based on photometric observations by MEarth and Spitzer as well as Doppler spectroscopy from HARPS, we report the discovery of an additional transiting rocky companion (LHS 1140 c) with a mass of $1.81\pm0.39~{\rm M_{Earth}}$ and a radius of $1.282\pm0.024~{\rm R_{Earth}}$ on a tighter, 3.77795-day orbit. We also obtain more precise estimates of the mass and radius of LHS 1140 b to be $6.98\pm0.89~{\rm M_{Earth}}$ and $1.727\pm0.032~{\rm R_{Earth}}$. The mean densities of planets b and c are $7.5\pm1.0~\rm{g/cm^3}$ and $4.7\pm1.1~\rm{g/cm^3}$, respectively, both consistent with the Earth's ratio of iron to magnesium silicate. The orbital eccentricities of LHS 1140 b and c are consistent with circular orbits and constrained to be below 0.06 and 0.31, respectively, with 90% confidence. Because the orbits of the two planets are co-planar and because we know from previous analyses of Kepler data that compact systems of small planets orbiting M dwarfs are commonplace, a search for more transiting planets in the LHS 1140 system could be fruitful. LHS 1140 c is one of the few known nearby terrestrial planets whose atmosphere could be studied with the upcoming James Webb Space Telescope.

The TRAPPIST-1 planetary system represents an exceptional opportunity for the atmospheric characterization of temperate terrestrial exoplanets with the upcoming James Webb Space Telescope (JWST). Assessing the potential impact of stellar contamination on the planets' transit transmission spectra is an essential precursor step to this characterization. Planetary transits themselves can be used to scan the stellar photosphere and to constrain its heterogeneity through transit depth variations in time and wavelength. In this context, we present our analysis of 169 transits observed in the optical from space with K2 and from the ground with the SPECULOOS and Liverpool telescopes. Combining our measured transit depths with literature results gathered in the mid/near-IR with Spitzer/IRAC and HST/WFC3, we construct the broadband transmission spectra of the TRAPPIST-1 planets over the 0.8-4.5 $\mu$m spectral range. While planets b, d, and f spectra show some structures at the 200-300ppm level, the four others are globally flat. Even if we cannot discard their instrumental origins, two scenarios seem to be favored by the data: a stellar photosphere dominated by a few high-latitude giant (cold) spots, or, alternatively, by a few small and hot (3500-4000K) faculae. In both cases, the stellar contamination of the transit transmission spectra is expected to be less dramatic than predicted in recent papers. Nevertheless, based on our results, stellar contamination can still be of comparable or greater order than planetary atmospheric signals at certain wavelengths. Understanding and correcting the effects of stellar heterogeneity therefore appears essential to prepare the exploration of TRAPPIST-1's with JWST.

We derive analytic, closed form, numerically stable solutions for the total flux received from a spherical planet, moon or star during an occultation if the specific intensity map of the body is expressed as a sum of spherical harmonics. Our expressions are valid to arbitrary degree and may be computed recursively for speed. The formalism we develop here applies to the computation of stellar transit light curves, planetary secondary eclipse light curves, and planet-planet/planet-moon occultation light curves, as well as thermal (rotational) phase curves. In this paper we also introduce STARRY, an open-source package written in C++ and wrapped in Python that computes these light curves. The algorithm in STARRY is six orders of magnitude faster than direct numerical integration and several orders of magnitude more precise. STARRY also computes analytic derivatives of the light curves with respect to all input parameters for use in gradient-based optimization and inference, such as Hamiltonian Monte Carlo (HMC), allowing users to quickly and efficiently fit observed light curves to infer properties of a celestial body's surface map.

We explore the connection between starspots and plages of three main-sequence stars by studying the chromospheric and photospheric activity over several rotation periods. We present simultaneous photometry and high-resolution ($R\sim 31,500$) spectroscopy of KIC 9652680, a young, superflare-producing G1 star with a rotation period of 1.4 days. Its Kepler light curve shows rotational modulation consistent with a bright hemisphere followed by a relatively dark hemisphere, generating photometric variability with a semi-amplitude of 4%. We find that KIC 9652680 is darkest when its $S$-index of Ca II H & K emission is at its maximum. We interpret this anti-correlation between flux and $S$ to indicate that dark starspots in the photosphere are co-located with the bright plages in the chromosphere, as they are on the Sun. Moving to lower masses and slower rotators, we present K2 observations with simultaneous spectroscopy of EPIC 211928486 (K5V) and EPIC 211966629 (K4V), two active stars in the 650 Myr-old open cluster Praesepe. The K2 photometry reveals that both stars have rotation periods of 11.7 days; while their flux varies by 1 and 2% respectively, their Ca II H & K $S$-indices seem to hold relatively constant as a function of rotational phase. This suggests that extended chromospheric networks of plages are not concentrated into regions of emission centered on the starspots that drive rotational modulation, unlike KIC 9652680. We also note that the Ca II emission of EPIC 211928486 dipped and recovered suddenly over the duration of one rotation, suggesting that the evolution timescale of plages may be of order the rotation period.

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.

We typically measure the radii of transiting exoplanets from the transit depth, which is given by the ratio of cross-sectional areas of the planet and star. However, if a star has dark starspots (or bright regions) distributed throughout the transit chord, the transit depth will be biased towards smaller (larger) values, and thus the inferred planet radius will be smaller (larger) if unaccounted for. We reparameterize the transit light curve to account for "self-contamination" by photospheric inhomogeneities by splitting the parameter $R_p/R_\star$ into two parameters: one for the radius ratio -- which controls the duration of ingress and egress -- and another which measures the possibly contaminated transit depth. We show that this is equivalent to the formulation for contamination by a second star (with positive or negative flux), and that it is sensitive to time-steady inhomogeneity of the stellar photosphere. We use synthetic light curves of spotted stars at high signal-to-noise to show that the radius recovered from measurement of the ingress/egress duration can recover the true radii of planets transiting spotted stars with axisymmetric spot distributions if the limb-darkening parameters are precisely known. We fit time-averaged high signal-to-noise transit light curves from Kepler and Spitzer of ten planets to measure the planet radii and search for evidence of spot distributions. We find that this sample has a range of measured depths and ingress durations which are self-consistent, providing no strong evidence for contamination by spots. However, there is suggestive evidence for occultation of starspots on Kepler-17, and that relatively bright regions are occulted by the planets of Kepler-412 and HD 80606. Future observations with the James Webb Space Telescope may enable this technique to yield accurate planetary radii in the presence of stellar inhomogeneities.

The Kepler Object of Interest Network (KOINet) is a multi-site network of telescopes around the globe organised to follow up transiting planet candidate KOIs with large transit timing variations (TTVs). Its main goal is to complete their TTV curves, as the Kepler telescope no longer observes the original Kepler field. Combining Kepler and new ground-based transit data we improve the modelling of these systems. To this end, we have developed a photodynamical model, and we demonstrate its performance using the Kepler-9 system as an example. Our comprehensive analysis combines the numerical integration of the system's dynamics over the time span of the observations along with the transit light curve model. This model is coupled with a Markov chain Monte Carlo algorithm, allowing the exploration of the model parameter space. Applied to the Kepler-9 long cadence data, short cadence data and 13 new transit observations collected by KOINet between the years 2014 to 2017, our modelling provides well constrained predictions for the next transits and the system's parameters. We have determined the densities of the planets Kepler-9b and 9c to the very precise values of rho_b = 0.439 +/-0.023 g/cm3 and rho_c = 0.322 +/- 0.017 g/cm3. Our analysis reveals that Kepler-9c will stop transiting in about 30 years. This results from strong dynamical interactions between Kepler-9b and 9c, near 2:1 resonance, that leads to a periodic change in inclination. Over the next 30 years the inclination of Kepler-9c (-9b) will decrease (increase) slowly. This should be measurable by a substantial decrease (increase) in the transit duration, in as soon as a few years' time. Observations that contradict this prediction might indicate the presence of additional objects. If this prediction proves true, this behaviour opens up a unique chance to scan the different latitudes of a star.

The Origins Space Telescope (OST) is one of four mission concepts currently being studied by NASA in preparation for the Astrophysics 2020 Decadal Survey. With active cooling (~4 K), OST will be sensitive in mid- to far-IR wavelengths, using imaging and spectroscopy to probe the furthest reaches of our galaxies, trace the path of water through star and planet formation, and place thermochemical constraints on the atmospheres of exoplanets ranging in size from Jupiter to Earth. This contribution to the Exoplanet Science Strategy committee discusses the significant advancements that the OST Mid-Infrared Imager, Spectrometer, and Coronagraph (MISC) instrument can make in studying cool planetary atmospheres. We particularly focus on the atmospheres of transiting rocky planets in the habitable zones of mid-to-late M stars. We discuss how OST thermal infrared observations can significantly enhance our understanding of the temperature structure and molecular abundances of biologically interesting gases on these worlds, including O3, CH4, H2O, and CO2.

The M8V star TRAPPIST-1 hosts seven roughly Earth-sized planets and is a promising target for exoplanet characterization. Kepler/K2 Campaign 12 observations of TRAPPIST-1 in the optical show an apparent rotational modulation with a 3.3 day period, though that rotational signal is not readily detected in the Spitzer light curve at 4.5 $\mu$m. If the rotational modulation is due to starspots, persistent dark spots can be excluded from the lack of photometric variability in the Spitzer light curve. We construct a photometric model for rotational modulation due to photospheric bright spots on TRAPPIST-1 which is consistent with both the Kepler and Spitzer light curves. The maximum-likelihood model with three spots has typical spot sizes of $R_\mathrm{spot}/R_\star \approx 0.004$ at temperature $T_\mathrm{spot} \gtrsim 5300 \pm 200$ K. We also find that large flares are observed more often when the brightest spot is facing the observer, suggesting a correlation between the position of the bright spots and flare events. In addition, these flares may occur preferentially when the spots are increasing in brightness, which suggests that the 3.3 d periodicity may not be a rotational signal, but rather a characteristic timescale of active regions.

Metal pollution in white dwarf atmospheres may be the accreted remnants of planetary objects. After the discovery of disintegrating planetary objects transiting WD 1145+017 (Vanderburg et al. 2015), undergraduates in the University of Washington's Pre- Major in Astronomy Program (Pre-MAP) were inspired to collect photometry of the brightest white dwarfs to hunt for similar transiting objects around other metal-polluted white dwarfs. Prior surveys have yet to make a detection of a transiting planet orbiting a white dwarf, yet white dwarfs are still an attractive target for searches of small, rocky planetary material. Since a typical white dwarf is Earth-sized, transits of Earth-sized planets could have depths >50%, so even low S/N photometry has a chance at discovering transiting material. We identified bright, northern, metal-polluted white dwarfs in the SDSS DR10 white dwarf catalog, and observed five targets with the Astrophysical Research Consortium Small Aperture Telescope (ARCSAT) 0.5-meter telescope at Apache Point Observatory. The ARCSAT photometry had sufficient precision to detect Moon-sized objects or larger at short orbital periods, though no such planets were detected for these targets. We look forward to surveys which may find planets orbiting white dwarfs, such as NASA's TESS, ESA's PLATO, and the Evryscope.

Astrometry from Gaia will measure the positions of stellar photometric centroids to unprecedented precision. We show that the precision of Gaia astrometry is sufficient to detect starspot-induced centroid jitter for nearby stars in the Tycho-Gaia Astrometric Solution (TGAS) sample with magnetic activity similar to the young G-star KIC 7174505 or the active M4 dwarf GJ 1243, but is insufficient to measure centroid jitter for stars with Sun-like spot distributions. We simulate Gaia observations of stars with 10 year activity cycles to search for evidence of activity cycles, and find that Gaia astrometry alone likely can not detect activity cycles for stars in the TGAS sample, even if they have spot distributions like KIC 7174505. We review the activity of the nearby low-mass stars in the TGAS sample for which we anticipate significant detections of spot-induced jitter.

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

Recent Kepler photometry has revealed that about half of white dwarfs (WDs) have periodic, low-level (~ 1e-4 - 1e-3), optical variations. Hubble Space Telescope (HST) ultraviolet spectroscopy has shown that up to about one half of WDs are actively accreting rocky planetary debris, as evidenced by the presence of photospheric metal absorption lines. We have obtained HST ultraviolet spectra of seven WDs that have been monitored for periodic variations, to test the hypothesis that these two phenomena are causally connected, i.e. that the optical periodic modulation is caused by WD rotation coupled with an inhomogeneous surface distribution of accreted metals. We detect photospheric metals in four out of the seven WDs. However, we find no significant correspondence between the existence of optical periodic variability and the detection of photospheric ultraviolet absorption lines. Thus the null hypothesis stands, that the two phenomena are not directly related. Some other source of WD surface inhomogeneity, perhaps related to magnetic field strength, combined with the WD rotation, or alternatively effects due to close binary companions, may be behind the observed optical modulation. We report the marginal detection of molecular hydrogen in WD J1949+4734, only the fourth known WD with detected H_2 lines. We also re-classify J1926+4219 as a carbon-rich He-sdO subdwarf.

During its four years of photometric observations, the Kepler space telescope detected thousands of exoplanets and exoplanet candidates. One of Kepler's greatest heritages has been the confirmation and characterization of hundreds of multi-planet systems via Transit Timing Variations (TTVs). However, there are many interesting candidate systems displaying TTVs on such long time scales that the existing Kepler observations are of insufficient length to confirm and characterize them by means of this technique. To continue with Kepler's unique work we have organized the "Kepler Object of Interest Network" (KOINet). The goals of KOINet are, among others, to complete the TTV curves of systems where Kepler did not cover the interaction timescales well. KOINet has been operational since March, 2014. Here we show some promising first results obtained from analyzing seven primary transits of KOI-0410.01, KOI-0525.01, KOI-0760.01, and KOI-0902.01 in addition to Kepler data, acquired during the first and second observing seasons of KOINet. While carefully choosing the targets we set demanding constraints about timing precision (at least 1 minute) and photometric precision (as good as 1 part per thousand) that were achieved by means of our observing strategies and data analysis techniques. For KOI-0410.01, new transit data revealed a turn-over of its TTVs. We carried out an in-depth study of the system, that is identified in the NASA's Data Validation Report as false positive. Among others, we investigated a gravitationally-bound hierarchical triple star system, and a planet-star system. While the simultaneous transit fitting of ground and space-based data allowed for a planet solution, we could not fully reject the three-star scenario. New data, already scheduled in the upcoming 2018 observing season, will set tighter constraints on the nature of the system.

TRAPPIST-1 hosts seven Earth-sized planets transiting an M8 star. We observed mid-transit times of each of the inner two planets with the Astrophysical Research Consortium (ARC) 3.5 m Telescope at Apache Point Observatory (APO) to help constrain the planet masses with transit timing variations, and we outline a procedure for analyzing transit observations of late-M stars with APO. The transit times of TRAPPIST-1 b and c are $\mathrm{BJD}_{\mathrm{TDB}} = 2457580.87634^{+0.00034}_{-0.00034}$ and $2457558.89477^{+0.00080}_{-0.00085}$, respectively, which will help constrain the planet masses.

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}$.

We present new 3.6 and 4.5 $\mu m$ Spitzer phase curves for the highly irradiated hot Jupiter WASP-33b and the unusually dense Saturn-mass planet HD 149026b. As part of this analysis, we develop a new variant of pixel level decorrelation that is effective at removing intrapixel sensitivity variations for long observations (>10 hours) where the position of the star can vary by a significant fraction of a pixel. Using this algorithm, we measure eclipse depths, phase amplitudes, and phase offsets for both planets at 3.6 $\mu m$ and 4.5 $\mu m$. We use a simple toy model to show that WASP-33b's phase offset, albedo, and heat recirculation efficiency are largely similar to those of other hot Jupiters despite its very high irradiation. On the other hand, our fits for HD 149026b prefer a very high albedo and an unusually high recirculation efficiency. We also compare our results to predictions from general circulation models, and find that while neither are a good match to the data, the discrepancies for HD 149026b are especially large. We speculate that this may be related to its high bulk metallicity, which could lead to enhanced atmospheric opacities and the formation of reflective cloud layers in localized regions of the atmosphere. We then place these two planets in a broader context by exploring relationships between the temperatures, albedos, heat transport efficiencies, and phase offsets of all planets with published thermal phase curves. We find a striking relationship between phase offset and irradiation temperature--the former drops with increasing temperature until around 3400 K, and rises thereafter. Although some aspects of this trend are mirrored in the circulation models, there are notable differences that provide important clues for future modeling efforts.

We explore the occurrence and detectability of planet-planet occultations (PPOs) in exoplanet systems. These are events during which a planet occults the disk of another planet in the same system, imparting a small photometric signal as its thermal or reflected light is blocked. We focus on the planets in TRAPPIST-1, whose orbital planes we show are aligned to within 0.3 degrees at 90% confidence. We present a photodynamical model for predicting and computing PPOs in TRAPPIST-1 and other systems for various assumptions of the planets' atmospheric states. When marginalizing over the uncertainties on all orbital parameters, we find that the rate of PPOs in TRAPPIST-1 is about 1.4 per day. We investigate the prospects for detection of these events with the James Webb Space Telescope, finding that ~10-20 occultations per year of b and c should be above the noise level at 12-15 microns. Joint modeling of several of these PPOs could lead to a robust detection. Alternatively, observations with the proposed Origins Space Telescope should be able to detect individual PPOs at high signal-to-noise. We show how PPOs can be used to break transit timing variation degeneracies, imposing strong constraints on the eccentricities and masses of the planets, as well as to constrain the longitudes of nodes and thus the complete three-dimensional structure of the system. We further show how modeling of these events can be used to reveal a planet's day/night temperature contrast and construct crude surface maps. We make our photodynamical code available on github.

Kepler photometry of the hot Neptune host star HAT-P-11 suggests that its spot latitude distribution is comparable to the Sun's near solar maximum. We search for evidence of an activity cycle in the CaII H & K chromospheric emission $S$-index with archival Keck/HIRES spectra and observations from the echelle spectrograph on the ARC 3.5 m Telescope at APO. The chromospheric emission of HAT-P-11 is consistent with a $\gtrsim 10$ year activity cycle, which plateaued near maximum during the Kepler mission. In the cycle that we observed, the star seemed to spend more time near active maximum than minimum. We compare the $\log R^\prime_{HK}$ normalized chromospheric emission index of HAT-P-11 with other stars. HAT-P-11 has unusually strong chromospheric emission compared to planet-hosting stars of similar effective temperature and rotation period, perhaps due to tides raised by its planet.

Transiting exoplanets in multi-planet systems have non-Keplerian orbits which can cause the times and durations of transits to vary. The theory and observations of transit timing variations (TTV) and transit duration variations (TDV) are reviewed. Since the last review, the Kepler spacecraft has detected several hundred perturbed planets. In a few cases, these data have been used to discover additional planets, similar to the historical discovery of Neptune in our own Solar System. However, the more impactful aspect of TTV and TDV studies has been characterization of planetary systems in which multiple planets transit. After addressing the equations of motion and parameter scalings, the main dynamical mechanisms for TTV and TDV are described, with citations to the observational literature for real examples. We describe parameter constraints, particularly the origin of the mass/eccentricity degeneracy and how it is overcome by the high-frequency component of the signal. On the observational side, derivation of timing precision and introduction to the timing diagram are given. Science results are reviewed, with an emphasis on mass measurements of transiting sub-Neptunes and super-Earths, from which bulk compositions may be inferred.

The growing field of large-scale time domain astronomy requires methods for probabilistic data analysis that are computationally tractable, even with large datasets. Gaussian Processes are a popular class of models used for this purpose but, since the computational cost scales, in general, as the cube of the number of data points, their application has been limited to small datasets. In this paper, we present a novel method for Gaussian Process modeling in one-dimension where the computational requirements scale linearly with the size of the dataset. We demonstrate the method by applying it to simulated and real astronomical time series datasets. These demonstrations are examples of probabilistic inference of stellar rotation periods, asteroseismic oscillation spectra, and transiting planet parameters. The method exploits structure in the problem when the covariance function is expressed as a mixture of complex exponentials, without requiring evenly spaced observations or uniform noise. This form of covariance arises naturally when the process is a mixture of stochastically-driven damped harmonic oscillators -- providing a physical motivation for and interpretation of this choice -- but we also demonstrate that it can be a useful effective model in some other cases. We present a mathematical description of the method and compare it to existing scalable Gaussian Process methods. The method is fast and interpretable, with a range of potential applications within astronomical data analysis and beyond. We provide well-tested and documented open-source implementations of this method in C++, Python, and Julia.

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.

Measuring precise planet masses, densities, and orbital dynamics in individual planetary systems is an important pathway toward understanding planet formation. The WASP-47 system has an unusual architecture that motivates a complex formation theory. The system includes a hot Jupiter ("b") neighbored by interior ("e") and exterior ("d") sub-Neptunes, and a long-period eccentric giant planet ("c"). We simultaneously modeled transit times from the Kepler K2 Mission and 118 radial velocities to determine precise masses, densities, and Keplerian orbital elements of the WASP-47 planets. Combining RVs and TTVs provides a better estimate of the mass of planet d (13.6\pm2.0~M_\odot) than obtained with only RVs ($12.75\pm2.70~M_\odot$) or TTVs ($16.1\pm3.8~M_\odot$). Planets e and d have high densities for their size, consistent with a history of photo-evaporation and/or formation in a volatile-poor environment. Through our RV and TTV analysis, we find that the planetary orbits have eccentricities similar to the solar system planets. The WASP-47 system has three similarities to our own solar system: (1) the planetary orbits are nearly circular and coplanar, (2) the planets are not trapped in mean motion resonances, and (3) the planets have diverse compositions. None of the current single-process exoplanet formation theories adequately reproduce these three characteristics of the WASP-47 system (or our solar system). We propose that WASP-47, like the solar system, formed in two stages: first, the giant planets formed in a gas-rich disk and migrated to their present locations, and second, the high-density sub-Neptunes formed in situ in a gas-poor environment.

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 present an update to the EVEREST K2 pipeline that addresses various limitations in the previous version and improves the photometric precision of the de-trended light curves. We develop a fast regularization scheme for third order pixel level decorrelation (PLD) and adapt the algorithm to include the PLD vectors of neighboring stars to enhance the predictive power of the model and minimize overfitting, particularly for faint stars. We also modify PLD to work for saturated stars and improve its performance on extremely variable stars. On average, EVEREST 2.0 light curves have 10-20% higher photometric precision than those in the previous version, yielding the highest precision light curves at all Kp magnitudes of any publicly available K2 catalog. For most K2 campaigns, we recover the original Kepler precision to at least Kp = 14, and to at least Kp = 15 for campaigns 1, 5, and 6. We also de-trend all short cadence targets observed by K2, obtaining even higher photometric precision for these stars. All light curves for campaigns 0-8 are available online in the EVEREST catalog, which will be continuously updated with future campaigns. EVEREST 2.0 is open source and is coded in a general framework that can be applied to other photometric surveys, including Kepler and the upcoming TESS mission.

We examine the feasibility of detecting auroral emission from the potentially habitable exoplanet Proxima Centauri b. Detection of aurorae would yield an independent confirmation of the planet's existence, constrain the presence and composition of its atmosphere, and determine the planet's eccentricity and inclination, thereby breaking the mass-inclination degeneracy. If Proxima Centauri b is a terrestrial world with an Earth-like atmosphere and magnetic field, we estimate the power at the 5577\AA\ OI auroral line is on the order of 0.1 TW under steady-state stellar wind, or ${\sim} 100 {\times}$ stronger than that on Earth. This corresponds to a planet-star contrast ratio of $10^{-6}-10^{-7}$ in a narrow band about the 5577\AA\ line, although higher contrast ($10^{-4}-10^{-5}$) may be possible during periods of strong magnetospheric disturbance (auroral power $1-10$ TW). We searched the Proxima Centauri b HARPS data for the 5577\AA\ line and for other prominent oxygen and nitrogen lines, but find no signal, indicating that the OI auroral line contrast must be lower than $2\times 10^{-2}$ (with power $\lesssim$ 3,000 TW), consistent with our predictions. We find that observations of 0.1 TW auroral emission lines are likely infeasible with current and planned telescopes. However, future observations with a space-based coronagraphic telescope or a ground-based extremely large telescope (ELT) with a coronagraph could push sensitivity down to terawatt oxygen aurorae (contrast $7\times 10^{-6}$) with exposure times of ${\sim} 1$ day. If a coronagraph design contrast of $10^{-7}$ can be achieved with negligible instrumental noise, a future concept ELT could observe steady-state auroral emission in a few nights.

Extrasolar planets on eccentric short-period orbits provide a laboratory in which to study radiative and tidal interactions between a planet and its host star under extreme forcing conditions. Studying such systems probes how the planet's atmosphere redistributes the time-varying heat flux from its host and how the host star responds to transient tidal distortion. Here, we report the insights into the planet-star interactions in HAT-P-2's eccentric planetary system gained from the analysis of 350 hr of 4.5 micron observations with the Spitzer Space Telescope. The observations show no sign of orbit-to-orbit variability nor of orbital evolution of the eccentric planetary companion, HAT-P-2 b. The extensive coverage allows us to better differentiate instrumental systematics from the transient heating of HAT-P-2 b's 4.5 micron photosphere and yields the detection of stellar pulsations with an amplitude of approximately 40 ppm. These pulsation modes correspond to exact harmonics of the planet's orbital frequency, indicative of a tidal origin. Transient tidal effects can excite pulsation modes in the envelope of a star, but, to date, such pulsations had only been detected in highly eccentric stellar binaries. Current stellar models are unable to reproduce HAT-P-2's pulsations, suggesting that our understanding of the interactions at play in this system is incomplete.

Planet searches using the radial velocity technique show a paucity of companions to solar-type stars within ~5 AU in the mass range of ~10 - 80 M$_{\text{Jup}}$. This deficit, known as the brown dwarf desert, currently has no conclusive explanation. New substellar companions in this region help asses the reality of the desert and provide insight to the formation and evolution of these objects. Here we present 10 new brown dwarf and two low-mass stellar companion candidates around solar-type stars from the Multi-object APO Radial-Velocity Exoplanet Large-Area Survey (MARVELS) of the Sloan Digital Sky Survey III (SDSS-III). These companions were selected from processed MARVELS data using the latest University of Florida Two Dimensional (UF2D) pipeline, which shows significant improvement and reduction of systematic errors over previous pipelines. The 10 brown dwarf companions range in mass from ~13 to 76 M$_{\text{Jup}}$ and have orbital radii of less than 1 AU. The two stellar companions have minimum masses of ~98 and 100 M$_{\text{Jup}}$. The host stars of the MARVELS brown dwarf sample have a mean metallicity of [Fe/H] = 0.03 $\pm$ 0.08 dex. Given our stellar sample we estimate the brown dwarf occurrence rate around solar-type stars with periods less than ~300 days to be ~0.56%.

The Kepler Mission has discovered thousands of exoplanets and revolutionized our understanding of their population. This large, homogeneous catalog of discoveries has enabled rigorous studies of the occurrence rate of exoplanets and planetary systems as a function of their physical properties. However, transit surveys like Kepler are most sensitive to planets with orbital periods much shorter than the orbital periods of Jupiter and Saturn, the most massive planets in our Solar System. To address this deficiency, we perform a fully automated search for long-period exoplanets with only one or two transits in the archival Kepler light curves. When applied to the $\sim 40,000$ brightest Sun-like target stars, this search produces 16 long-period exoplanet candidates. Of these candidates, 6 are novel discoveries and 5 are in systems with inner short-period transiting planets. Since our method involves no human intervention, we empirically characterize the detection efficiency of our search. Based on these results, we measure the average occurrence rate of exoplanets smaller than Jupiter with orbital periods in the range 2-25 years to be $2.0\pm0.7$ planets per Sun-like star.

We report the detections of a giant planet (MARVELS-7b) and a brown dwarf candidate (MARVELS-7c) around the primary star in the close binary system, HD 87646. It is the first close binary system with more than one substellar circum-primary companion discovered to the best of our knowledge. The detection of this giant planet was accomplished using the first multi-object Doppler instrument (KeckET) at the Sloan Digital Sky Survey (SDSS) telescope. Subsequent radial velocity observations using ET at Kitt Peak National Observatory, HRS at HET, the "Classic" spectrograph at the Automatic Spectroscopic Telescope at Fairborn Observatory, and MARVELS from SDSS-III confirmed this giant planet discovery and revealed the existence of a long-period brown dwarf in this binary. HD 87646 is a close binary with a separation of $\sim22$ AU between the two stars, estimated using the Hipparcos catalogue and our newly acquired AO image from PALAO on the 200-inch Hale Telescope at Palomar. The primary star in the binary, HD 87646A, has Teff = 5770$\pm$80K, log(g)=4.1$\pm$0.1 and [Fe/H] = $-0.17\pm0.08$. The derived minimum masses of the two substellar companions of HD 87646A are 12.4$\pm$0.7M$_{\rm Jup}$ and 57.0$\pm3.7$M$_{\rm Jup}$. The periods are 13.481$\pm$0.001 days and 674$\pm$4 days and the measured eccentricities are 0.05$\pm$0.02 and 0.50$\pm$0.02 respectively. Our dynamical simulations show the system is stable if the binary orbit has a large semi-major axis and a low eccentricity, which can be verified with future astrometry observations.

We present EVEREST, an open-source pipeline for removing instrumental noise from K2 light curves. EVEREST employs a variant of pixel level decorrelation (PLD) to remove systematics introduced by the spacecraft's pointing error and a Gaussian process (GP) to capture astrophysical variability. We apply EVEREST to all K2 targets in campaigns 0-7, yielding light curves with precision comparable to that of the original Kepler mission for stars brighter than $K_p \approx 13$, and within a factor of two of the Kepler precision for fainter targets. We perform cross-validation and transit injection and recovery tests to validate the pipeline, and compare our light curves to the other de-trended light curves available for download at the MAST High Level Science Products archive. We find that EVEREST achieves the highest average precision of any of these pipelines for unsaturated K2 stars. The improved precision of these light curves will aid in exoplanet detection and characterization, investigations of stellar variability, asteroseismology, and other photometric studies. The EVEREST pipeline can also easily be applied to future surveys, such as the TESS mission, to correct for instrumental systematics and enable the detection of low signal-to-noise transiting exoplanets. The EVEREST light curves and the source code used to generate them are freely available online.

We present new ways to identify single and multiple moons around extrasolar planets using planetary transit timing variations (TTVs) and transit duration variations (TDVs). For planets with one moon, measurements from successive transits exhibit a hitherto undescribed pattern in the TTV-TDV diagram, originating from the stroboscopic sampling of the planet's orbit around the planet-moon barycenter. This pattern is fully determined and analytically predictable after three consecutive transits. The more measurements become available, the more the TTV-TDV diagram approaches an ellipse. For planets with multiple moons in orbital mean motion resonance (MMR), like the Galilean moons, the pattern is much more complex and addressed numerically in this report. Exomoons in MMR can also form closed, predictable TTV-TDV figures if the drift of the moons' pericenters is sufficiently slow. We find that MMR exomoons produce loops in the TTV-TDV diagram and that the number of these loops is equal to the order of the MMR, or the largest integer in the MMR ratio. We use a Bayesian model and Monte Carlo simulations to test the discoverability of exomoons using TTV-TDV diagrams with current and near-future technology. In a blind test, two of us (BP, DA) successfully retrieved a large moon from simulated TTV-TDV by co-authors MH and RH, which resembled data from a known Kepler planet candidate. Single exomoons with a 10% moon-to-planet mass ratio, like to Pluto-Charon binary, can be detectable in the archival data of the Kepler primary mission. Multi-exomoon systems, however, require either larger telescopes or brighter target stars. Complementary detection methods invoking a moon's own photometric transit or its orbital sampling effect can be used for validation or falsification. A combination of TESS, CHEOPS, and PLATO data would offer a compelling opportunity for exomoon discoveries around bright stars.

We analyze full-orbit phase curve observations of the transiting hot Jupiters WASP-19b and HAT-P-7b at 3.6 and 4.5 $\mu$m obtained using the Spitzer Space Telescope. For WASP-19b, we measure secondary eclipse depths of $0.485\%\pm 0.024\%$ and $0.584\%\pm 0.029\%$ at 3.6 and 4.5 $\mu$m, which are consistent with a single blackbody with effective temperature $2372 \pm 60$ K. The measured 3.6 and 4.5 $\mu$m secondary eclipse depths for HAT-P-7b are $0.156\%\pm 0.009\%$ and $0.190\%\pm 0.006\%$, which are well-described by a single blackbody with effective temperature $2667\pm 57$ K. Comparing the phase curves to the predictions of one-dimensional and three-dimensional atmospheric models, we find that WASP-19b's dayside emission is consistent with a model atmosphere with no dayside thermal inversion and moderately efficient day-night circulation. We also detect an eastward-shifted hotspot, suggesting the presence of a superrotating equatorial jet. In contrast, HAT-P-7b's dayside emission suggests a dayside thermal inversion and relatively inefficient day-night circulation; no hotspot shift is detected. For both planets, these same models do not agree with the measured nightside emission. The discrepancies in the model-data comparisons for WASP-19b might be explained by high-altitude silicate clouds on the nightside and/or high atmospheric metallicity, while the very low 3.6 $\mu$m nightside planetary brightness for HAT-P-7b may be indicative of an enhanced global C/O ratio. We compute Bond albedos of 0 ($<0.08$ at $1\sigma$) and $0.38\pm 0.06$ for WASP-19b and HAT-P-7b, respectively. In the context of other planets with thermal phase curve measurements, we show that WASP-19b and HAT-P-7b fit the general trend of decreasing day-night heat recirculation with increasing irradiation.

We derive the transit timing variations (TTVs) of two planets near a second order mean motion resonance on nearly circular orbits. We show that the TTVs of each planet are given by sinusoids with a frequency of $j n_2-(j-2)n_1$, where $j \ge 3$ is an integer characterizing the resonance and $n_2$ and $n_1$ are the mean motions of the outer and inner planets, respectively. The amplitude of the TTV depends on the mass of the perturbing planet, relative to the mass of the star, and on both the eccentricities and longitudes of pericenter of each planet. The TTVs of the two planets are approximated anti-correlated, with phases of $\phi$ and $\approx \phi+\pi$, where the phase $\phi$ also depends on the eccentricities and longitudes of pericenter. Therefore, the TTVs caused by proximity to a second order mean motion resonance do not in general uniquely determine both planet masses, eccentricities, and pericenters. This is completely analogous to the case of TTVs induced by two planets near a first order mean motion resonance. We explore how other TTV signals, such as the short-period synodic TTV or a first order resonant TTV, in combination with the second order resonant TTV, can break degeneracies. Lastly, we derive approximate formulae for the TTVs of planets near any order eccentricity-type mean motion resonance; this shows that the same basic sinusoidal TTV structure holds for all eccentricity-type resonances. Our general formula reduces to previously derived results near first order mean motion resonances.

Characterization of transiting planets with transit timing variations (TTVs) requires understanding how to translate the observed TTVs into masses and orbital elements of the planets. This can be challenging in multi-planet transiting systems, but fortunately these systems tend to be nearly plane-parallel and low eccentricity. Here we present a novel derivation of analytic formulae for TTVs that are accurate to first order in the planet-star mass ratios and in the orbital eccentricities. These formulae are accurate in proximity to first order resonances, as well as away from resonance, and compare well with more computationally expensive N-body integrations in the low eccentricity, low mass-ratio regime when applied to simulated and to actual multi-transiting Kepler planet systems. We make code available for implementing these formulae.

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 infer dynamical masses in eight multi-planet systems using transit times measured from Kepler's complete dataset, including short-cadence data where available. Of the eighteen dynamical masses that we infer, ten pass multiple tests for robustness. These are in systems; Kepler-26 (KOI-250), Kepler-29 (KOI-738), Kepler-60 (KOI-2086), Kepler-105 (KOI-115), and Kepler-307 (KOI-1576). Kepler-105 c has a density consistent with an Earth-like composition. Strong TTV signals were detected from additional planets, but their inferred masses were sensitive to outliers or consistent solutions could not be found with independently-measured transit times, including planets at; Kepler-49 (KOI-248), Kepler-57 (KOI-1270), Kepler-105 (KOI-115) and Kepler-177 (KOI-523). Strong upper limits on the mass of Kepler-177 c imply an extremely low density ~0.1 g cm$^{-3}$. In most cases, individual orbital eccentricities were poorly constrained due to degeneracies in TTV inversion. For five planet pairs in our sample, strong secular interactions imply a moderate-to-high likelihood of apsidal alignment over a wide range of possible eccentricities. We also find solutions for the three planets known to orbit Kepler-60 in a Laplace-like resonance chain. However, non-librating solutions also match the transit-timing data. For six systems, we calculate more precise stellar parameters than previously known, enabling useful constraints on planetary densities where we have robust mass measurements. Placing these exoplanets on the mass-radius diagram, we find a wide range of densities is observed among sub-Neptune mass planets and that the range in observed densities is anti-correlated with incident flux.

As lower-mass stars often host multiple rocky planets, gravitational interactions among planets can have significant effects on climate and habitability over long timescales. Here we explore a specific case, Kepler-62f, a potentially habitable planet in a five-planet system with a K2V host star. N-body integrations reveal the stable range of initial eccentricities for Kepler-62f is $0.00 \leqslant e \leqslant 0.32$, absent the effect of additional, undetected planets. We simulate the tidal evolution of Kepler-62f in this range and find that, for certain assumptions, the planet can be locked in a synchronous rotation state. Simulations using LMD Generic GCM indicate that with 3 bars of CO$_2$ in its atmosphere, Kepler-62f would only be warm enough for surface liquid water at the upper limit of this eccentricity range, providing it has a high planetary obliquity (between 60$^\circ$ and 90$^\circ$). A climate similar to modern-day Earth is possible for the entire range of stable eccentricities if atmospheric CO$_2$ is increased to 5-bar levels. In a low-CO$_2$ case, simulations with CCSM4 and LMD Generic GCM indicate that increases in planetary obliquity and orbital eccentricity coupled with an orbital configuration that places the summer solstice at or near pericenter permit regions of the planet with above-freezing surface temperatures. This may melt ice sheets formed during colder seasons. If Kepler-62f is synchronously rotating and has an ocean, CO$_2$ levels above 3 bars would be required to distribute enough heat to the night side of the planet to avoid atmospheric freeze-out and permit a large enough region of open water at the planet's substellar point to remain stable. Overall, we find multiple plausible combinations of orbital and atmospheric properties that permit surface liquid water on Kepler-62f.

We report the discovery of the sixth known eclipsing double white dwarf (WD) system, SDSS J1152+0248, with a 2.3968 +/- 0.0003 h orbital period, in data from the Kepler Mission's K2 continuation. Analysing and modelling the K2 data together with ground-based fast photometry, spectroscopy, and radial-velocity measurements, we determine that the primary is a DA-type WD with mass M1 = 0.47 +/- 0.11 Msun, radius R1 = 0.0197 +/- 0.0035 Rsun, and cooling age t1 = 52 +/- 36 Myr. No lines are detected, to within our sensitivity, from the secondary WD, but it is likely also of type DA. Its central surface brightness, as measured from the secondary eclipse, is 0.31 of the primary's surface brightness. Its mass, radius, and cooling age, respectively, are M2 = 0.44 +/- 0.09 Msun, R2 = 0.0223 +0.0064 -0.0050 Rsun, and t2 = 230 +/- 100 Myr. SDSS J1152+0248 is a near twin of the double-lined eclipsing WD system CSS 41177.

We examine the light and color evolution of the T Tauri binary KH 15D through photometry obtained at wavelengths between 0.55 and 8.0 $\mu$m. The data were collected with ANDICAM on the 1.3 m SMARTS telescope at Cerro-Tololo Inter-American Observatory and with IRAC on the Spitzer Space Telescope. We show that the system's circumbinary ring, which acts as a screen that covers and uncovers different portions of the binary orbit as the ring precesses, has reached an orientation where the brighter component (star B) fully or nearly fully emerges during each orbital cycle. The fainter component (star A) remains fully occulted by the screen at all phases. The leading and trailing edges of the screen move across the sky at the same rate of $\sim$15 meters per second, consistent with expectation for a ring with a radius and width of $\sim$4 AU and a precession period of $\sim$6500 years. Light and color variations continue to indicate that the screen is sharp edged and opaque at \emph{VRIJH} wavelengths. However, we find an increasing transparency of the ring edge at 2.2, 3.6, and 4.5 $\mu$m. Reddening seen at the beginning of the eclipse that occurred during the CSI 2264 campaign particularly suggests selective extinction by a population of large dust grains. Meanwhile, the gradual bluing observed while star B is setting is indicative of forward scattering effects at the edge of the ring. The SED of the system at its bright phase shows no evidence of infrared excess emission that can be attributed to radiation from the ring or other dust component out to 8 microns.

Determining which small exoplanets have stony-iron compositions is necessary for quantifying the occurrence of such planets and for understanding the physics of planet formation. Kepler-10 hosts the stony-iron world Kepler-10b (K10b), and also contains what has been reported to be the largest solid silicate-ice planet, Kepler-10c (K10c). Using 220 radial velocities (RVs), including 72 precise RVs from Keck-HIRES of which 20 are new from 2014-2015, and 17 quarters of Kepler photometry, we obtain the most complete picture of the Kepler-10 system to date. We find that K10b (Rp=1.47 Re) has mass 3.72$\pm$0.42 Me and density 6.46$\pm$0.73 g/cc. Modeling the interior of K10b as an iron core overlaid with a silicate mantle, we find that the iron core constitutes 0.17$\pm$0.11 of the planet mass. For K10c (Rp=2.35 Re) we measure Mp=13.98$\pm$1.79 Me and $\rho$=5.94$\pm$0.76 g/cc, significantly lower than the mass computed in Dumusque et al. (2014, 17.2$\pm$1.9 Me). Internal compositional modeling reveals that at least $10\%$ of the radius of Kepler-10c is a volatile envelope composed of hydrogen-helium ($0.2\%$ of the mass, $16\%$ of the radius) or super-ionic water ($28\%$ of the mass, $29\%$ of the radius). Analysis of only HIRES data yields a higher mass for K10b and a lower mass for K10c than does analysis of the HARPS-N data alone, with the mass estimates for K10c formally inconsistent by 3$\sigma$. Splitting the RVs from each instrument leads to inconsistent measurements for the mass of planet c in each data set. This suggests that time-correlated noise is present and that the uncertainties in the planet masses (especially K10c) exceed our formal estimates. Transit timing variations (TTVs) of K10c indicate the likely presence of a third planet in the system, KOI-72.X. The TTVs and RVs are consistent with KOI-72.X having an orbital period of 24, 71, or 101 days, and a mass from 1-7 Me.

This is a joint summary of the reports from the three Astrophysics Program Analysis Groups (PAGs) in response to the "Planning for the 2020 Decadal Survey" charge given by the Astrophysics Division Director Paul Hertz. This joint executive summary contains points of consensus across all three PAGs. Additional findings specific to the individual PAGs are reported separately in the individual reports. The PAGs concur that all four large mission concepts identified in the white paper as candidates for maturation prior to the 2020 Decadal Survey should be studied in detail. These include the Far-IR Surveyor, the Habitable-Exoplanet Imaging Mission, the UV/Optical/IR Surveyor, and the X-ray Surveyor. This finding is predicated upon assumptions outlined in the white paper and subsequent charge, namely that 1) major development of future large flagship missions under consideration are to follow the implementation phases of JWST and WFIRST; 2) NASA will partner with the European Space Agency on its L3 Gravitational Wave Surveyor; 3) the Inflation Probe be classified as a probe-class mission to be developed according to the 2010 Decadal Survey report. If these key assumptions were to change, this PAG finding would need to be re-evaluated. The PAGs find that there is strong community support for the second phase of this activity - maturation of the four proposed mission concepts via Science and Technology Definition Teams (STDTs). The PAGs find that there is strong consensus that all of the STDTs contain broad and interdisciplinary representation of the science community. Finally, the PAGs find that there is community support for a line of Probe-class missions within the Astrophysics mission portfolio (condensed).

The primary Kepler Mission provided nearly continuous monitoring of ~200,000 objects with unprecedented photometric precision. We present the final catalog of eclipsing binary systems within the 105 square degree Kepler field of view. This release incorporates the full extent of the data from the primary mission (Q0-Q17 Data Release). As a result, new systems have been added, additional false positives have been removed, ephemerides and principal parameters have been recomputed, classifications have been revised to rely on analytical models, and eclipse timing variations have been computed for each system. We identify several classes of systems including those that exhibit tertiary eclipse events, systems that show clear evidence of additional bodies, heartbeat systems, systems with changing eclipse depths, and systems exhibiting only one eclipse event over the duration of the mission. We have updated the period and galactic latitude distribution diagrams and included a catalog completeness evaluation. The total number of identified eclipsing and ellipsoidal binary systems in the Kepler field of view has increased to 2878, 1.3% of all observed Kepler targets. An online version of this catalog with downloadable content and visualization tools is maintained at http://keplerEBs.villanova.edu.

We utilize multi-dimensional simulations of varying equatorial jet strength to predict wavelength dependent variations in the eclipse times of gas-giant planets. A displaced hot-spot introduces an asymmetry in the secondary eclipse light curve that manifests itself as a measured offset in the timing of the center of eclipse. A multi-wavelength observation of secondary eclipse, one probing the timing of barycentric eclipse at short wavelengths and another probing at longer wavelengths, will reveal the longitudinal displacement of the hot-spot and break the degeneracy between this effect and that associated with the asymmetry due to an eccentric orbit. The effect of time offsets was first explored in the IRAC wavebands by Williams et. al (2006). Here we improve upon their methodology, extend to a broad ranges of wavelengths, and demonstrate our technique on a series of multi-dimensional radiative-hydrodynamical simulations of HD 209458b with varying equatorial jet strength and hot-spot displacement. Simulations with the largest hot-spot displacement result in timing offsets of up to 100 seconds in the infrared. Though we utilize a particular radiative hydrodynamical model to demonstrate this effect, the technique is model independent. This technique should allow a much larger survey of hot-spot displacements with JWST then currently accessible with time-intensive phase curves, hopefully shedding light on the physical mechanisms associated with thermal energy advection in irradiated gas-giants.

Direct imaging of extrasolar planets with future space-based coronagraphic telescopes may provide a means of detecting companion moons at wavelengths where the moon outshines the planet. We propose a detection strategy based on the positional variation of the center of light with wavelength, "spectroastrometry." This new application of this technique could be used to detect an exomoon, to determine the exomoon's orbit and the mass of the host exoplanet, and to disentangle of the spectra of the planet and moon. We consider two model systems, for which we discuss the requirements for detection of exomoons around nearby stars. We simulate the characterization of an Earth-Moon analog system with spectroastrometry, showing that the orbit, the planet mass, and the spectra of both bodies can be recovered. To enable the detection and characterization of exomoons we recommend that coronagraphic telescopes should extend in wavelength coverage to 3 micron, and should be designed with spectroastrometric requirements in mind.

We present full-orbit phase curve observations of the eccentric ($e\sim 0.08$) transiting hot Jupiter WASP-14b obtained in the 3.6 and 4.5 $\mu$m bands using the \textit{Spitzer Space Telescope}. We use two different methods for removing the intrapixel sensitivity effect and compare their efficacy in decoupling the instrumental noise. Our measured secondary eclipse depths of $0.1882\%\pm 0.0048\%$ and $0.2247\%\pm 0.0086\%$ at 3.6 and 4.5 $\mu$m, respectively, are both consistent with a blackbody temperature of $2402\pm 35$ K. We place a $2\sigma$ upper limit on the nightside flux at 3.6 $\mu$m and find it to be $9\%\pm 1\%$ of the dayside flux, corresponding to a brightness temperature of 1079 K. At 4.5 $\mu$m, the minimum planet flux is $30\%\pm 5\%$ of the maximum flux, corresponding to a brightness temperature of $1380\pm 65$ K. We compare our measured phase curves to the predictions of one-dimensional radiative transfer and three-dimensional general circulation models. We find that WASP-14b's measured dayside emission is consistent with a model atmosphere with equilibrium chemistry and a moderate temperature inversion. These same models tend to over-predict the nightside emission at 3.6 $\mu$m, while under-predicting the nightside emission at 4.5 $\mu$m. We propose that this discrepancy might be explained by an enhanced global C/O ratio. In addition, we find that the phase curves of WASP-14b ($7.8 M_{\mathrm{Jup}}$) are consistent with a much lower albedo than those of other Jovian mass planets with thermal phase curve measurements, suggesting that it may be emitting detectable heat from the deep atmosphere or interior processes.

In this work we present Spitzer 3.6 and 4.5 micron secondary eclipse observations of five new cool (<1200 K) transiting gas giant planets: HAT-P-19b, WASP-6b, WASP-10b, WASP-39b, and WASP-67b. We compare our measured eclipse depths to the predictions of a suite of atmosphere models and to eclipse depths for planets with previously published observations in order to constrain the temperature- and mass-dependent properties of gas giant planet atmospheres. We find that the dayside emission spectra of planets less massive than Jupiter require models with efficient circulation of energy to the night side and/or increased albedos, while those with masses greater than that of Jupiter are consistently best-matched by models with inefficient circulation and low albedos. At these relatively low temperatures we expect the atmospheric methane to CO ratio to vary as a function of metallicity, and we therefore use our observations of these planets to constrain their atmospheric metallicities. We find that the most massive planets have dayside emission spectra that are best-matched by solar metallicity atmosphere models, but we are not able to place strong constraints on metallicities of the smaller planets in our sample. Interestingly, we find that the ratio of the 3.6 and 4.5 micron brightness temperatures for these cool transiting planets is independent of planet temperature, and instead exhibits a tentative correlation with planet mass. If this trend can be confirmed, it would suggest that the shape of these planets' emission spectra depends primarily on their masses, consistent with the hypothesis that lower-mass planets are more likely to have metal-rich atmospheres.

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.

The Kepler mission has yielded a large number of planet candidates from among the Kepler Objects of Interest (KOIs), but spectroscopic follow-up of these relatively faint stars is a serious bottleneck in confirming and characterizing these systems. We present motivation and survey design for an ongoing project with the SDSS-III multiplexed APOGEE near-infrared spectrograph to monitor hundreds of KOI host stars. We report some of our first results using representative targets from our sample, which include current planet candidates that we find to be false positives, as well as candidates listed as false positives that we do not find to be spectroscopic binaries. With this survey, KOI hosts are observed over ~20 epochs at a radial velocity precision of 100-200 m/s. These observations can easily identify a majority of false positives caused by physically-associated stellar or substellar binaries, and in many cases, fully characterize their orbits. We demonstrate that APOGEE is capable of achieving RV precision at the 100-200 m/s level over long time baselines, and that APOGEE's multiplexing capability makes it substantially more efficient at identifying false positives due to binaries than other single-object spectrographs working to confirm KOIs as planets. These APOGEE RVs enable ancillary science projects, such as studies of fundamental stellar astrophysics or intrinsically rare substellar companions. The coadded APOGEE spectra can be used to derive stellar properties (T_eff, log(g)) and chemical abundances of over a dozen elements to probe correlations of planet properties with individual elemental abundances.

Gravitational interactions between planets in transiting exoplanetary systems lead to variations in the times of transit that are diagnostic of the planetary masses and the dynamical state of the system. Here we show that synodic "chopping" contributions to these transit timing variations (TTVs) can be used to uniquely measure the masses of planets without full dynamical analyses involving direct integration of the equations of motion. We present simple analytic formulae for the chopping signal, which are valid (generally <10% error) for modest eccentricities e <~ 0.1. Importantly, these formulae primarily depend on the mass of the perturbing planet, and therefore the chopping signal can be used to break the mass/free-eccentricity degeneracy which can appear for systems near first order mean motion resonances. Using a harmonic analysis, we apply these TTV formulae to a number of Kepler systems which had been previously analyzed with full dynamical analyses. We show that when chopping is measured, the masses of both planets can be determined uniquely, in agreement with previous results, but without the need for numerical orbit integrations. This demonstrates how mass measurements from TTVs may primarily arise from an observable chopping signal. The formula for chopping can also be used to predict the number of transits and timing precision required for future observations, such as those made by TESS or PLATO, in order to infer planetary masses through analysis of TTVs.

We report the discovery of one newly confirmed planet ($P=66.06$ days, $R_{\rm{P}}=2.68\pm0.17R_\oplus$) and mass determinations of two previously validated Kepler planets, Kepler-289 b ($P=34.55$ days, $R_{\rm{P}}=2.15\pm0.10R_\oplus$) and Kepler-289-c ($P=125.85$ days, $R_{\rm{P}}=11.59\pm0.10R_\oplus$), through their transit timing variations (TTVs). We also exclude the possibility that these three planets reside in a $1:2:4$ Laplace resonance. The outer planet has very deep ($\sim1.3%$), high signal-to-noise transits, which puts extremely tight constraints on its host star's stellar properties via Kepler's Third Law. The star PH3 is a young ($\sim1$ Gyr as determined by isochrones and gyrochronology), Sun-like star with $M_*=1.08\pm0.02M_\odot$, $R_*=1.00\pm0.02R_\odot$, and $T_{\rm{eff}}=5990\pm38$ K. The middle planet's large TTV amplitude ($\sim5$ hours) resulted either in non-detections or inaccurate detections in previous searches. A strong chopping signal, a shorter period sinusoid in the TTVs, allows us to break the mass-eccentricity degeneracy and uniquely determine the masses of the inner, middle, and outer planets to be $M=7.3\pm6.8M_\oplus$, $4.0\pm0.9M_\oplus$, and $M=132\pm17M_\oplus$, which we designate PH3 b, c, and d, respectively. Furthermore, the middle planet, PH3 c, has a relatively low density, $\rho=1.2\pm0.3$ g/cm$^3$ for a planet of its mass, requiring a substantial H/He atmosphere of $2.1^{+0.8}_{-0.3}%$ by mass, and joins a growing population of low-mass, low-density planets.

We report secondary eclipse photometry of the hot Jupiter XO-3b in the 4.5~$\mu$m band taken with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. We measure individual eclipse depths and center of eclipse times for a total of twelve secondary eclipses. We fit these data simultaneously with two transits observed in the same band in order to obtain a global best-fit secondary eclipse depth of $0.1580\pm 0.0036\%$ and a center of eclipse phase of $0.67004\pm 0.00013 $. We assess the relative magnitude of variations in the dayside brightness of the planet by measuring the size of the residuals during ingress and egress from fitting the combined eclipse light curve with a uniform disk model and place an upper limit of 0.05$\%$. The new secondary eclipse observations extend the total baseline from one and a half years to nearly three years, allowing us to place an upper limit on the periastron precession rate of $2.9\times 10^{-3}$ degrees/day the tightest constraint to date on the periastron precession rate of a hot Jupiter. We use the new transit observations to calculate improved estimates for the system properties, including an updated orbital ephemeris. We also use the large number of secondary eclipses to obtain the most stringent limits to date on the orbit-to-orbit variability of an eccentric hot Jupiter and demonstrate the consistency of multiple-epoch Spitzer observations.

The hot Jupiter HD 209458b is particularly amenable to detailed study as it is among the brightest transiting exoplanet systems currently known (V-mag = 7.65; K-mag = 6.308) and has a large planet-to-star contrast ratio. HD 209458b is predicted to be in synchronous rotation about its host star with a hot spot that is shifted eastward of the substellar point by superrotating equatorial winds. Here we present the first full-orbit observations of HD 209458b, in which its 4.5 $\mu$m emission was recorded with $Spitzer$/IRAC. Our study revises the previous 4.5 $\mu$m measurement of HD 209458b's secondary eclipse emission downward by $\sim$35% to $0.1391%^{+0.0072%}_{-0.0069%}$, changing our interpretation of the properties of its dayside atmosphere. We find that the hot spot on the planet's dayside is shifted eastward of the substellar point by $40.9^{\circ}\pm{6.0^{\circ}}$, in agreement with circulation models predicting equatorial superrotation. HD 209458b's dayside (T$_{bright}$ = 1499 $\pm$ 15 K) and nightside (T$_{bright}$ = 972 $\pm$ 44 K) emission indicates a day-to-night brightness temperature contrast smaller than that observed for more highly irradiated exoplanets, suggesting that the day-to-night temperature contrast may be partially a function of the incident stellar radiation. The observed phase curve shape deviates modestly from global circulation model predictions potentially due to disequilibrium chemistry or deficiencies in the current hot CH$_{4}$ line lists used in these models. Observations of the phase curve at additional wavelengths are needed in order to determine the possible presence and spatial extent of a dayside temperature inversion, as well as to improve our overall understanding of this planet's atmospheric circulation.

We report on the orbital architectures of Kepler systems having multiple planet candidates identified in the analysis of data from the first six quarters of Kepler data and reported by Batalha et al. (2013). These data show 899 transiting planet candidates in 365 multiple-planet systems and provide a powerful means to study the statistical properties of planetary systems. Using a generic mass-radius relationship, we find that only two pairs of planets in these candidate systems (out of 761 pairs total) appear to be on Hill-unstable orbits, indicating ~96% of the candidate planetary systems are correctly interpreted as true systems. We find that planet pairs show little statistical preference to be near mean-motion resonances. We identify an asymmetry in the distribution of period ratios near first-order resonances (e.g., 2:1, 3:2), with an excess of planet pairs lying wide of resonance and relatively few lying narrow of resonance. Finally, based upon the transit duration ratios of adjacent planets in each system, we find that the interior planet tends to have a smaller transit impact parameter than the exterior planet does. This finding suggests that the mode of the mutual inclinations of planetary orbital planes is in the range 1.0-2.2 degrees, for the packed systems of small planets probed by these observations.

Over 40% of Sun-like stars are bound in binary or multistar systems. Stellar remnants in edge-on binary systems can gravitationally magnify their companions, as predicted 40 years ago. By using data from the Kepler spacecraft, we report the detection of such a "self-lensing" system, in which a 5-hour pulse of 0.1% amplitude occurs every orbital period. The white dwarf stellar remnant and its Sun-like companion orbit one another every 88.18 days, a long period for a white dwarf-eclipsing binary. By modeling the pulse as gravitational magnification (microlensing) along with Kepler's laws and stellar models, we constrain the mass of the white dwarf to be ~63% of the mass of our Sun. Further study of this system, and any others discovered like it, will help to constrain the physics of white dwarfs and binary star evolution.

Kepler-13Ab (= KOI-13.01) is a unique transiting hot Jupiter. It is one of very few known short-period planets orbiting a hot A-type star, making it one of the hottest planets currently known. The availability of Kepler data allows us to measure the planet's occultation (secondary eclipse) and phase curve in the optical, which we combine with occultations observed by warm Spitzer at 4.5 mic and 3.6 mic and a ground-based occultation observation in the Ks band (2.1 mic). We derive a day-side hemisphere temperature of 2,750 +- 160 K as the effective temperature of a black body showing the same occultation depths. Comparing the occultation depths with one-dimensional planetary atmosphere models suggests the presence of an atmospheric temperature inversion. Our analysis shows evidence for a relatively high geometric albedo, Ag= 0.33 +0.04 -0.06. While measured with a simplistic method, a high Ag is supported also by the fact that the one-dimensional atmosphere models underestimate the occultation depth in the optical. We use stellar spectra to determine the dilution, in the four wide bands where occultation was measured, due to the visual stellar binary companion 1.15 +- 0.05" away. The revised stellar parameters measured using these spectra are combined with other measurements leading to revised planetary mass and radius estimates of Mp = 4.94 - 8.09 Mjup and Rp = 1.406 +- 0.038 Rjup. Finally, we measure a Kepler mid-occultation time that is 34.0 +- 6.9 s earlier than expected based on the mid-transit time and the delay due to light travel time, and discuss possible scenarios.

We report the detection and characterization of two short-period, Neptune-sized planets around the active host star Kepler-210. The host star's parameters derived from those planets are (a) mutually inconsistent and (b) do not conform to the expected host star parameters. We furthermore report the detection of transit timing variations (TTVs) in the O-C diagrams for both planets. We explore various scenarios that explain and resolve those discrepancies. A simple scenario consistent with all data appears to be one that attributes substantial eccentricities to the inner short-period planets and that interprets the TTVs as due to the action of another, somewhat longer period planet. To substantiate our suggestions, we present the results of N-body simulations that modeled the TTVs and that checked the stability of the Kepler-210 system.

Transit timing variations (TTVs) have proven to be a powerful technique for confirming Kepler planet candidates, for detecting non-transiting planets, and for constraining the masses and orbital elements of multi-planet systems. These TTV applications often require the numerical integration of orbits for computation of transit times (as well as impact parameters and durations); frequently tens of millions to billions of simulations are required when running statistical analyses of the planetary system properties. We have created a fast code for transit timing computation, TTVFast, which uses a symplectic integrator with a Keplerian interpolator for the calculation of transit times (Nesvorny et al. 2013). The speed comes at the expense of accuracy in the calculated times, but the accuracy lost is largely unnecessary, as transit times do not need to be calculated to accuracies significantly smaller than the measurement uncertainties on the times. The time step can be tuned to give sufficient precision for any particular system. We find a speed-up of at least an order of magnitude relative to dynamical integrations with high precision using a Bulirsch-Stoer integrator.

The Kepler mission has discovered over 2500 exoplanet candidates in the first two years of spacecraft data, with approximately 40% of them in candidate multi-planet systems. The high rate of multiplicity combined with the low rate of identified false-positives indicates that the multiplanet systems contain very few false-positive signals due to other systems not gravitationally bound to the target star (Lissauer, J. J., et al., 2012, ApJ 750, 131). False positives in the multi- planet systems are identified and removed, leaving behind a residual population of candidate multi-planet transiting systems expected to have a false-positive rate less than 1%. We present a sample of 340 planetary systems that contain 851 planets that are validated to substantially better than the 99% confidence level; the vast majority of these have not been previously verified as planets. We expect ~2 unidentified false-positives making our sample of planet very reliable. We present fundamental planetary properties of our sample based on a comprehensive analysis of Kepler light curves and ground-based spectroscopy and high-resolution imaging. Since we do not require spectroscopy or high-resolution imaging for validation, some of our derived parameters for a planetary system may be systematically incorrect due to dilution from light due to additional stars in the photometric aperture. None the less, our result nearly doubles the number of verified exoplanets.

We extend the statistical analysis of Lissauer et al. (2012, ApJ 750, 112), which demonstrates that the overwhelming majority of Kepler candidate multiple transiting systems (multis) represent true transiting planets, and develop therefrom a procedure to validate large numbers of planet candidates in multis as bona fide exoplanets. We show that this statistical framework correctly estimates the abundance of false positives already identified around Kepler targets with multiple sets of transit-like signatures based on their abundance around targets with single sets of transit-like signatures. We estimate the number of multis that represent split systems of one or more planets orbiting each component of a binary star system. We use the high reliability rate for multis to validate more than one dozen particularly interesting multi-planet systems are validated in a companion paper by Rowe et al. (2014, ApJ, this issue). We note that few very short period (P < 1.6 days) planets orbit within multiple transiting planet systems and discuss possible reasons for their absence. There also appears to be a shortage of planets with periods exceeding a few months in multis.

We report on the masses, sizes, and orbits of the planets orbiting 22 Kepler stars. There are 49 planet candidates around these stars, including 42 detected through transits and 7 revealed by precise Doppler measurements of the host stars. Based on an analysis of the Kepler brightness measurements, along with high-resolution imaging and spectroscopy, Doppler spectroscopy, and (for 11 stars) asteroseismology, we establish low false-positive probabilities for all of the transiting planets (41 of 42 have a false-positive probability under 1%), and we constrain their sizes and masses. Most of the transiting planets are smaller than 3X the size of Earth. For 16 planets, the Doppler signal was securely detected, providing a direct measurement of the planet's mass. For the other 26 planets we provide either marginal mass measurements or upper limits to their masses and densities; in many cases we can rule out a rocky composition. We identify 6 planets with densities above 5 g/cc, suggesting a mostly rocky interior for them. Indeed, the only planets that are compatible with a purely rocky composition are smaller than ~2 R_earth. Larger planets evidently contain a larger fraction of low-density material (H, He, and H2O).

The past three decades have seen prodigious advances in astronomy and astrophysics. Beginning with the exploration of our solar system and continuing through the pioneering Explorers and Great Observatories of today, NASA missions have made essential contributions to these advances. This roadmap presents a science-driven 30-year vision for the future of NASA Astrophysics that builds on these achievements to address some of our most ancient and fundamental questions: Are we alone? How did we get here? How does the universe work? The search for the answers constitutes the Enduring Quests of this roadmap. Building on the priorities identified in New Worlds, New Horizons, we envision future science investigations laid out in three Eras, with each representing roughly ten years of mission development in a given field. The immediate Near-Term Era covers ongoing NASA-led activities and planned missions. This will be followed by the missions of the Formative Era, which will build on the preceding technological developments and scientific discoveries, with remarkable capabilities that will enable breakthroughs across the landscape of astrophysics. These will then lay the foundations for the Daring Visions of the Visionary Era: missions and explorations that will take us deep into unchartered scientific and technological terrain. The roadmap outlined herein will require the vision and wherewithal to undertake highly ambitious programs over the next 30 years. The discoveries that emerge will inspire generations of citizen scientists young and old, and inspire all of humanity for decades to come.

The well-known bluer-when-brighter trend observed in quasar variability is a signature of the complex processes in the accretion disk, and can be a probe of the quasar variability mechanism. Using a sample of 604 variable quasars with repeat spectra in SDSS-I/II, we construct difference spectra to investigate the physical causes of this bluer-when-brighter trend. The continuum of our composite difference spectrum is well-fit by a power-law, with a spectral index in excellent agreement with previous results. We measure the spectral variability relative to the underlying spectra of the quasars, which is independent of any extinction, and compare to model predictions. We show that our SDSS spectral variability results cannot be produced by global accretion rate fluctuations in a thin disk alone. However, we find that a simple model of a inhomogeneous disk with localized temperature fluctuations will produce power-law spectral variability over optical wavelengths. We show that the inhomogeneous disk will provide good fits to our observed spectral variability if the disk has large temperature fluctuations in many independently varying zones, in excellent agreement with independent constraints from quasar microlensing disk sizes, their strong UV spectral continuum, and single-band variability amplitudes. Our results provide an independent constraint on quasar variability models, and add to the mounting evidence that quasar accretion disks have large localized temperature fluctuations.

Unlike hot Jupiters or other gas giants, super-Earths are expected to have a wide variety of compositions, ranging from terrestrial bodies like our own to more gaseous planets like Neptune. Observations of transiting systems, which allow us to directly measure planet masses and radii and constrain atmospheric properties, are key to understanding the compositional diversity of the planets in this mass range. Although Kepler has discovered hundreds of transiting super-Earth candidates over the past four years, the majority of these planets orbit stars that are too far away and too faint to allow for detailed atmospheric characterization and reliable mass estimates. Ground-based transit surveys focus on much brighter stars, but most lack the sensitivity to detect planets in this size range. One way to get around the difficulty of finding these smaller planets in transit is to start by choosing targets that are already known to contain super-Earth sized bodies detected using the radial velocity technique. Here we present results from a Spitzer program to observe six of the most favorable RV detected super-Earth systems, including HD 1461, HD 7924, HD 156668, HIP 57274, and GJ 876. We find no evidence for transits in any of their 4.5 micron flux light curves, and place limits on the allowed transit depths and corresponding planet radii that rule out even the most dense and iron-rich compositions for these objects. We also observed HD 97658, but the observation window was based on a possible ground-based transit detection (Henry et al. 2011) that was later ruled out; thus the window did not include the predicted time for the transit detection recently made by MOST (Dragomir et al. 2013).

We report secondary eclipse photometry of two hot Jupiters, WASP-48b and HAT-P-23b, at 3.6 and 4.5 um taken with the InfraRed Array Camera aboard the Spitzer Space Telescope during the warm Spitzer mission and in the H and Ks bands with the Wide Field IR Camera at the Palomar 200-inch Hale Telescope. WASP-48b and HAT-P-23b are Jupiter-mass and twice Jupiter-mass objects orbiting an old, slightly evolved F star and an early G dwarf star, respectively. In the H, Ks, 3.6 um, and 4.5 um bands, respectively, we measure secondary eclipse depths of 0.047% +/- 0.016%, 0.109% +/- 0.027%, 0.176% +/- 0.013%, and 0.214% +/- 0.020% for WASP-48b. In the Ks, 3.6 um, and 4.5 um bands, respectively, we measure secondary eclipse depths of 0.234% +/- 0.046%, 0.248% +/- 0.019%, and 0.309% +/- 0.026% for HAT-P-23b. For WASP-48b and HAT-P-23b, respectively, we measure delays of 2.6 +/- 3.9 minutes and 4.0 +/- 2.4 minutes relative to the predicted times of secondary eclipse for circular orbits, placing 2-sigma upper limits on |e cos(omega)| of 0.0053 and 0.0080, both of which are consistent with circular orbits. The dayside emission spectra of these planets are well-described by blackbodies with effective temperatures of 2158 +/- 100 K (WASP-48b) and 2154 +/- 90 K (HAT-P-23b), corresponding to moderate recirculation in the zero albedo case. Our measured eclipse depths are also consistent with one-dimensional radiative transfer models featuring varying degrees of recirculation and weak thermal inversions or no inversions at all. We discuss how the absence of strong temperature inversions on these planets may be related to the activity levels and metallicities of their host stars.

Since 2005, the Pre-Major in Astronomy Program (Pre-MAP) at the University of Washington (UW) Department of Astronomy has made a concentrated effort to recruit and retain underrepresented undergraduates in science, technology, engineering and mathematics (STEM). This paper evaluates Pre-MAP in the context of the larger UW student population using data compiled by the University's student database. We evaluate the Pre-MAP program in terms of our goals of recruiting a more diverse population than the University and in terms of a higher fraction of students successfully completing degrees. We find that Pre-MAP serves a higher percentage of underrepresented minorities and equal percentages of women compared to entering freshmen classes at UW. Additionally, Pre-MAP has a higher percentage of degree completion with higher average GPA's and similar time to completion when compared to UW as a whole and other STEM majors, particularly with students that place into lower-level math courses (such as basic algebra or pre-calculus).

Variations in the timing of transiting exoplanets provide a powerful tool detecting additional planets in the system. Thus, the aim of this paper is to discuss the plausibility of transit timing variations on the Qatar-1 system by means of primary transit light curves analysis. Furthermore, we provide an interpretation of the timing variation. We observed Qatar-1 between March 2011 and October 2012 using the 1.2 m OLT telescope in Germany and the 0.6 m PTST telescope in Spain. We present 26 primary transits of the hot Jupiter Qatar-1b. In total, our light curves cover a baseline of 18 months. We report on indications for possible long-term transit timing variations (TTVs). Assuming that these TTVs are true, we present two different scenarios that could explain them. Our reported $\sim$ 190 days TTV signal can be reproduced by either a weak perturber in resonance with Qatar-1b, or by a massive body in the brown dwarf regime. More observations and radial velocity monitoring are required to better constrain the perturber's characteristics. We also refine the ephemeris of Qatar-1b, which we find to be \mbox{$T_0 = 2456157.42204 \pm 0.0001$ \bjdtdb} and \mbox{$P = 1.4200246 \pm 0.0000007$ days}, and improve the system orbital parameters.

A large fraction of white dwarfs (WDs) may host planets in their habitable zones. These planets may provide our best chance to detect bio-markers on a transiting exoplanet, thanks to the diminished contrast ratio between the Earth-sized WD and its Earth-sized planets. The JWST is capable of obtaining the first spectroscopic measurements of such planets, yet there are no known planets around WDs. Here we propose to take advantage of the unique capability of the Kepler spacecraft in the 2-Wheels mode to perform a transit survey that is capable of identifying the first planets in the habitable zone of a WD. We propose to obtain Kepler time-series photometry of 10,000 WDs in the SDSS imaging area to search for planets in the habitable zone. Thanks to the large field of view of Kepler, for the first time in history, a large number of WDs can be observed at the same time, which is essential for discovering transits. Our proposed survey requires a total of 200 days of observing time, and will find up to 100 planets in the WD habitable zone. This survey will maintain Kepler's spirit of searching for habitable Earths, but near new hosts. With few-day observations and minute-cadences per field, it will also open up a completely unexplored discovery space. In addition to planets, this survey is sensitive to pulsating WDs, as well as eclipsing short period stellar and substellar companions. These have important implications for constraining the double WD merger rate and their contribution to Type Ia supernovae and the gravitational wave foreground. Given the relatively low number density of our targets, this program can be combined with other projects that would benefit from high cadence and many-fields observations with Kepler, e.g. a transit survey of a magnitude-limited, complete sample of nearby M dwarfs or asteroseismology of variable stars (e.g. RR Lyrae) in the same fields.

Exoplanetary transmission spectroscopy in the near-infrared using Hubble/NICMOS is currently ambiguous because different observational groups claim different results from the same data, depending on their analysis methodologies. Spatial scanning with Hubble/WFC3 provides an opportunity to resolve this ambiguity. We here report WFC3 spectroscopy of the giant planets HD209458b and XO-1b in transit, using spatial scanning mode for maximum photon-collecting efficiency. We introduce an analysis technique that derives the exoplanetary transmission spectrum without the necessity of explicitly decorrelating instrumental effects, and achieves nearly photon-limited precision even at the high flux levels collected in spatial scan mode. Our errors are within 6-percent (XO-1) and 26-percent (HD209458b) of the photon-limit at a spectral resolving power of 70, and are better than 0.01-percent per spectral channel. Both planets exhibit water absorption of approximately 200 ppm at the water peak near 1.38 microns. Our result for XO-1b contradicts the much larger absorption derived from NICMOS spectroscopy. The weak water absorption we measure for HD209458b is reminiscent of the weakness of sodium absorption in the first transmission spectroscopy of an exoplanet atmosphere by Charbonneau et al. (2002). Model atmospheres having uniformly-distributed extra opacity of 0.012 cm^2 per gram account approximately for both our water measurement and the sodium absorption in this planet. Our results for HD209458b support the picture advocated by Pont et al. (2013) in which weak molecular absorptions are superposed on a transmission spectrum that is dominated by continuous opacity due to haze and/or dust. However, the extra opacity needed for HD209458b is grayer than for HD189733b, with a weaker Rayleigh component.

We report the discovery of a candidate brown dwarf or a very low mass stellar companion (MARVELS-5b) to the star HIP 67526 from the Multi-object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The radial velocity curve for this object contains 31 epochs spread over 2.5 years. Our Keplerian fit using a Markov Chain Monte Carlo approach, reveals that the companion has an orbital period of $90.2695^{+0.0188}_{-0.0187}$ days, an eccentricity of $0.4375 \pm 0.0040$ and a semi-amplitude of $2948.14^{+16.65}_{-16.55}$ m s$^{-1}$. Using additional high-resolution spectroscopy, we find the host star has an effective temperature $T_{\rm{eff}}=6004 \pm 34$ K, a surface gravity $\log g$ [cgs] $=4.55 \pm 0.17$ and a metallicity [Fe/H] $=+0.04 \pm 0.06$. The stellar mass and radius determined through the empirical relationship of Torres et al. (2010), yields 1.10$\pm$0.09 $M_{\sun}$ and 0.92$\pm$0.19 $R_{\sun}$. The minimum mass of MARVELS-5b is $65.0 \pm 2.9 M_{Jup}$, indicating that it is likely to be either a brown dwarf or a very low mass star, thus occupying a relatively sparsely-populated region of the mass function of companions to solar-type stars. The distance to this system is 101$\pm$10 pc from the astrometric measurements of Hipparcos. No stellar tertiary is detected in the high-contrast images taken by either FastCam lucky imaging or Keck adaptive optics imaging, ruling out any star with mass greater than 0.2$M_{\sun}$ at a separation larger than 40 AU.

We present photometry of the extrasolar planet WASP-5b in the 3.6 and 4.5 micron bands taken with the Spitzer Space Telescope's Infrared Array Camera as part of the extended warm mission. By examining the depth of the planet's secondary eclipse at these two wavelengths, we can place joint constraints on the planet's atmospheric pressure-temperature profile and chemistry. We measure secondary eclipse depths of 0.197% +/- 0.028% and 0.227% +/- 0.025% in the 3.6 micron and 4.5 micron bands, respectively. Our observations are best matched by models showing a hot dayside and, depending on our choice of model, a weak thermal inversion or no inversion at all. We measure a mean offset from the predicted center of eclipse of 0.078 +/- 0.032 hours, translating to ecos(omega) = 0.0031 +/- 0.0013 and consistent with a circular orbit. We see no evidence for any eclipse timing variations comparable to those reported in a previous transit study.

We present 11.2 micron observations of the gravitationally lensed, radio-loud z_s=2.64 quasar MG0414+0534, obtained using the Michelle camera on Gemini North. We find a flux ratio anomaly of A2/A1= 0.93 +/- 0.02 for the quasar images A1 and A2. When combined with the 11.7 micron measurements from Minezaki et al. (2009), the A2/A1 flux ratio is nearly 5-sigma from the expected ratio for a model based on the two visible lens galaxies. The mid-IR flux ratio anomaly can be explained by a satellite (substructure), 0.3" Northeast of image A2, as can the detailed VLBI structures of the jet produced by the quasar. When we combine the mid-IR flux ratios with high-resolution VLBI measurements, we find a best-fit mass between 10^(6.2) and 10^(7.5) M_sol inside the Einstein radius for a satellite substructure modeled as a singular isothermal sphere at the redshift of the main lens (z_l=0.96). We are unable to set an interesting limit on the mass to light ratio due to its proximity to the quasar image A2. While the observations used here were technically difficult, surveys of flux anomalies in gravitational lenses with the James Webb Space Telescope will be simple, fast, and should well constrain the abundance of substructure in dark matter haloes.

Following Ford et al. (2011, 2012) and Steffen et al. (2012) we derived the transit timing of 1960 Kepler KOIs using the pre-search data conditioning (PDC) light curves of the first twelve quarters of the Kepler data. For 721 KOIs with large enough SNRs, we obtained also the duration and depth of each transit. The results are presented as a catalog for the community to use. We derived a few statistics of our results that could be used to indicate significant variations. Including systems found by previous works, we have found 130 KOIs that showed highly significant TTVs, and 13 that had short-period TTV modulations with small amplitudes. We consider two effects that could cause apparent periodic TTV - the finite sampling of the observations and the interference with the stellar activity, stellar spots in particular. We briefly discuss some statistical aspects of our detected TTVs. We show that the TTV period is correlated with the orbital period of the planet and with the TTV amplitude.

The Kepler-11 planetary system contains six transiting planets ranging in size from 1.8 to 4.2 times the radius of Earth. Five of these planets orbit in a tightly-packed configuration with periods between 10 and 47 days. We perform a dynamical analysis of the system based upon transit timing variations observed in more than three years of \ik photometric data. Stellar parameters are derived using a combination of spectral classification and constraints on the star's density derived from transit profiles together with planetary eccentricity vectors provided by our dynamical study. Combining masses of the planets relative to the star from our dynamical study and radii of the planets relative to the star from transit depths together with deduced stellar properties yields measurements of the radii of all six planets, masses of the five inner planets, and an upper bound to the mass of the outermost planet, whose orbital period is 118 days. We find mass-radius combinations for all six planets that imply that substantial fractions of their volumes are occupied by constituents that are less dense than rock. The Kepler-11 system contains the lowest mass exoplanets for which both mass and radius have been measured.

We report the discovery of a highly eccentric, double-lined spectroscopic binary star system (TYC 3010-1494-1), comprising two solar-type stars that we had initially identified as a single star with a brown dwarf companion. At the moderate resolving power of the MARVELS spectrograph and the spectrographs used for subsequent radial-velocity (RV) measurements (R ~ <30,000), this particular stellar binary mimics a single-lined binary with an RV signal that would be induced by a brown dwarf companion (Msin(i)~50 M_Jup) to a solar-type primary. At least three properties of this system allow it to masquerade as a single star with a very low-mass companion: its large eccentricity (e~0.8), its relatively long period (P~238 days), and the approximately perpendicular orientation of the semi-major axis with respect to the line of sight (omega~189 degrees). As a result of these properties, for ~95% of the orbit the two sets of stellar spectral lines are completely blended, and the RV measurements based on centroiding on the apparently single-lined spectrum is very well fit by an orbit solution indicative of a brown dwarf companion on a more circular orbit (e~0.3). Only during the ~5% of the orbit near periastron passage does the true, double-lined nature and large RV amplitude of ~15 km/s reveal itself. The discovery of this binary system is an important lesson for RV surveys searching for substellar companions; at a given resolution and observing cadence, a survey will be susceptible to these kinds of astrophysical false positives for a range of orbital parameters. Finally, for surveys like MARVELS that lack the resolution for a useful line bisector analysis, it is imperative to monitor the peak of the cross-correlation function for suspicious changes in width or shape, so that such false positives can be flagged during the candidate vetting process.

Since the discovery of the first exoplanet we have known that other planetary systems can look quite unlike our own. However, until recently we have only been able to probe the upper range of the planet size distribution. The high precision of the Kepler space telescope has allowed us to detect planets that are the size of Earth and somewhat smaller, but no previous planets have been found that are smaller than those we see in our own Solar System. Here we report the discovery of a planet significantly smaller than Mercury. This tiny planet is the innermost of three planets that orbit the Sun-like host star, which we have designated Kepler-37. Owing to its extremely small size, similar to that of Earth's Moon, and highly irradiated surface, Kepler-37b is probably a rocky planet with no atmosphere or water, similar to Mercury.

We present Warm Spitzer/IRAC secondary eclipse time series photometry of three short-period transiting exoplanets, HAT-P-3b, HAT-P-4b and HAT-P-12b, in both the available 3.6 and 4.5 micron bands. HAT-P-3b and HAT-P-4b are Jupiter-mass, objects orbiting an early K and an early G dwarf stars, respectively. For HAT-P-3b we find eclipse depths of 0.112%+0.015%-0.030% (3.6 micron) and 0.094%+0.016%-0.009% (4.5 micron). The HAT-P-4b values are 0.142%+0.014%-0.016% (3.6 micron) and 0.122%+0.012%-0.014% (4.5micron). The two planets' photometry is consistent with inefficient heat redistribution from their day to night sides (and low albedos), but it is inconclusive about possible temperature inversions in their atmospheres. HAT-P-12b is a Saturn-mass planet and is one of the coolest planets ever observed during secondary eclipse, along with hot Neptune GJ 436b and hot Saturn WASP-29b. We are able to place 3$\sigma$ upper limits on the secondary eclipse depth of HAT-P-12b in both wavelengths: < 0.042% (3.6 micron) and <0.085% (4.5 micron). We discuss these results in the context of the {\it Spitzer} secondary eclipse measurements of GJ 436b and WASP-29b. It is possible that we do not detect the eclipses of HAT-P-12b due to high eccentricity, but find that weak planetary emission in these wavelengths is a more likely explanation. We place 3 sigma upper limits on the |e cos(omega)| quantity (where e is eccentricity and omega is the argument of periapsis) for HAT-P-3b (<0.0081) and HAT-P-4b (<0.0042), based on the secondary eclipse timings.

We present the detection of five planets -- Kepler-62b, c, d, e, and f -- of size 1.31, 0.54, 1.95, 1.61 and 1.41 Earth radii, orbiting a K2V star at periods of 5.7, 12.4, 18.2, 122.4 and 267.3 days, respectively. The outermost planets (Kepler-62e & -62f) are super-Earth-size (1.25 < planet radius/earth radius < 2.0) planets in the habitable zone (HZ) of their host star, receiving 1.2 +- 0.2 and 0.41 +- 0.05 times the solar flux at Earth's orbit. Theoretical models of Kepler-62e and -62f for a stellar age of ~7 Gyr suggest that both planets could be solid: either with a rocky composition or composed of mostly solid water in their bulk.

The Apache Point Survey of Transit Lightcurves of Exoplanets (APOSTLE) observed 10 transits of XO-2b over a period of three years. We present measurements which confirm previous estimates of system parameters like the normalized semi-major axis (a/R_{*}), stellar density (\rho_{*}), impact parameter (b) and orbital inclination (i_{orb}). Our errors on system parameters like a/R_{*} and \rho_{*} have improved by ~40% compared to previous best ground-based measurements. Our study of the transit times show no evidence for transit timing variations and we are able to rule out co-planar companions with masses \ge 0.20 \mearth\ in low order mean motion resonance with XO-2b. We also explored the stability of the XO-2 system given various orbital configurations of a hypothetical planet near the 2:1 mean motion resonance. We find that a wide range of orbits (including Earth-mass perturbers) are both dynamically stable and produce observable TTVs. We find that up to 51% of our stable simulations show TTVs that are smaller than the typical transit timing errors (~20 sec) measured for XO-2b, and hence remain undetectable.

We describe the discovery of a likely brown dwarf (BD) companion with a minimum mass of 31.7 +/- 2.0 M_Jup to GSC 03546-01452 from the MARVELS radial velocity survey, which we designate as MARVELS-6b. For reasonable priors, our analysis gives a probability of 72% that MARVELS-6b has a mass below the hydrogen-burning limit of 0.072 M_Sun, and thus it is a high-confidence BD companion. It has a moderately long orbital period of 47.8929 +0.0063/-0.0062 days with a low eccentricty of 0.1442 +0.0078/-0.0073, and a semi-amplitude of 1644 +12/-13 m/s. Moderate resolution spectroscopy of the host star has determined the following parameters: T_eff = 5598 +/- 63, log g = 4.44 +/- 0.17, and [Fe/H] = +0.40 +/- 0.09. Based upon these measurements, GSC 03546-01452 has a probable mass and radius of M_star = 1.11 +/- 0.11 M_Sun and R_star = 1.06 +/- 0.23 R_Sun with an age consistent with less than ~6 Gyr at a distance of 219 +/- 21 pc from the Sun. Although MARVELS-6b is not observed to transit, we cannot definitively rule out a transiting configuration based on our observations. There is a visual companion detected with Lucky Imaging at 7.7 arcsec from the host star, but our analysis shows that it is not bound to this system. The minimum mass of MARVELS-6b exists at the minimum of the mass functions for both stars and planets, making this a rare object even compared to other BDs.

We present the first secondary eclipse and phase curve observations for the highly eccentric hot Jupiter HAT-P-2b in the 3.6, 4.5, 5.8, and 8.0 \mu m bands of the Spitzer Space Telescope. The 3.6 and 4.5 \mu m data sets span an entire orbital period of HAT-P-2b, making them the longest continuous phase curve observations obtained to date and the first full-orbit observations of a planet with an eccentricity exceeding 0.2. We present an improved non-parametric method for removing the intrapixel sensitivity variations in Spitzer data at 3.6 and 4.5 \mu m that robustly maps position-dependent flux variations. We find that the peak in planetary flux occurs at 4.39+/-0.28, 5.84+/-0.39, and 4.68+/-0.37 hours after periapse passage with corresponding maxima in the planet/star flux ratio of 0.1138%+/-0.0089%, 0.1162%+/-0.0080%, and 0.1888%+/-0.0072% in the 3.6, 4.5, and 8.0 \mu m bands respectively. We compare our measured secondary eclipse depths to the predictions from a one-dimensional radiative transfer model, which suggests the possible presence of a transient day side inversion in HAT-P-2b's atmosphere near periapse. We also derive improved estimates for the system parameters, including its mass, radius, and orbital ephemeris. Our simultaneous fit to the transit, secondary eclipse, and radial velocity data allows us to determine the eccentricity and argument of periapse of HAT-P-2b's orbit with a greater precision than has been achieved for any other eccentric extrasolar planet. We also find evidence for a long-term linear trend in the radial velocity data. This trend suggests the presence of another substellar companion in the HAT-P-2 system, which could have caused HAT-P-2b to migrate inward to its present-day orbit via the Kozai mechanism.

We present EXOFAST, a fast, robust suite of routines written in IDL which is designed to fit exoplanetary transits and radial velocity variations simultaneously or separately, and characterize the parameter uncertainties and covariances with a Differential Evolution Markov Chain Monte Carlo method. We describe how our code incorporates both data sets to simultaneously derive stellar parameters along with the transit and RV parameters, resulting in more self-consistent results on an example fit of the discovery data of HAT-P-3b that is well-mixed in under five minutes on a standard desktop computer. We describe in detail how our code works and outline ways in which the code can be extended to include additional effects or generalized for the characterization of other data sets -- including non-planetary data sets. We discuss the pros and cons of several common ways to parameterize eccentricity, highlight a subtle mistake in the implementation of MCMC that could bias the inferred eccentricity of intrinsically circular orbits to significantly non-zero results, discuss a problem with IDL's built-in random number generator in its application to large MCMC fits, and derive a method to analytically fit the linear and quadratic limb darkening coefficients of a planetary transit. Finally, we explain how we achieved improved accuracy and over a factor of 100 improvement in the execution time of the transit model calculation. Our entire source code, along with an easy-to-use online interface for several basic features of our transit and radial velocity fitting, are available online at http://astroutils.astronomy.ohio-state.edu/exofast .

Planets embedded within dust disks may drive the formation of large scale clumpy dust structures by trapping dust into resonant orbits. Detection and subsequent modeling of the dust structures would help constrain the mass and orbit of the planet and the disk architecture, give clues to the history of the planetary system, and provide a statistical estimate of disk asymmetry for future exoEarth-imaging missions. Here we present the first search for these resonant structures in the inner regions of planetary systems by analyzing the light curves of hot Jupiter planetary candidates identified by the Kepler mission. We detect only one candidate disk structure associated with KOI 838.01 at the 3-sigma confidence level, but subsequent radial velocity measurements reveal that KOI 838.01 is a grazing eclipsing binary and the candidate disk structure is a false positive. Using our null result, we place an upper limit on the frequency of dense exozodi structures created by hot Jupiters. We find that at the 90% confidence level, less than 21% of Kepler hot Jupiters create resonant dust clumps that lead and trail the planet by ~90 degrees with optical depths >~5*10^-6, which corresponds to the resonant structure expected for a lone hot Jupiter perturbing a dynamically cold dust disk 50 times as dense as the zodiacal cloud.

We present transit observations of the WASP-2 exoplanet system by the Apache Point Survey of Transit Lightcurves of Exoplanets (APOSTLE) program. Model fitting to these data allows us to improve measurements of the hot-Jupiter exoplanet WASP-2b and its orbital parameters by a factor of ~2 over prior studies; we do not find evidence for transit depth variations. We do find reduced chi^2 values greater than 1.0 in the observed minus computed transit times. A sinusoidal fit to the residuals yields a timing semi-amplitude of 32 seconds and a period of 389 days. However, random rearrangements of the data provide similar quality fits, and we cannot with certainty ascribe the timing variations to mutual exoplanet interactions. This inconclusive result is consistent with the lack of incontrovertible transit timing variations (TTVs) observed in other hot-Jupiter systems. This outcome emphasizes that unique recognition of TTVs requires dense sampling of the libration cycle (e.g. continuous observations from space-based platforms). However, even in systems observed with the Kepler spacecraft, there is a noted lack of transiting companions and TTVs in hot-Jupiter systems. This result is more meaningful, and indicates that hot-Jupiter systems, while they are easily observable from the ground, do not appear to be currently configured in a manner favorable to the detection of TTVs. The future of ground-based TTV studies may reside in resolving secular trends, and/or implementation at extreme quality observing sites to minimize atmospheric red noise.

The Apache Point Survey of Transit Lightcurves of Exoplanets (APOSTLE) observed eleven transits of TrES-3b over two years in order to constrain system parameters and look for transit timing and depth variations. We describe an updated analysis protocol for APOSTLE data, including the reduction pipeline, transit model and Markov Chain Monte Carlo analyzer. Our estimates of the system parameters for TrESb are consistent with previous estimates to within the 2\sigma\ confidence level. We improved the errors (by 10--30%) on system parameters like the orbital inclination ($i_{\text{orb}}$), impact parameter (b) and stellar density (\rho$_{\star}$) compared to previous measurements. The near-grazing nature of the system, and incomplete sampling of some transits, limited our ability to place reliable uncertainties on individual transit depths and hence we do not report strong evidence for variability. Our analysis of the transit timing data show no evidence for transit timing variations and our timing measurements are able to rule out Super-Earth and Gas Giant companions in low order mean motion resonance with TrES-3b.

We present an eccentric, short-period brown dwarf candidate orbiting the active, slightly evolved subgiant star TYC 2087-00255-1, which has effective temperature T_eff = 5903+/-42 K, surface gravity log (g) = 4.07+/-0.16 (cgs), and metallicity [Fe/H] = -0.23+/-0.07. This candidate was discovered using data from the first two years of the Multi-object APO Radial Velocity Exoplanets Large-area Survey (MARVELS), which is part of the third phase of Sloan Digital Sky Survey. From our 38 radial velocity measurements spread over a two-year time baseline, we derive a Keplerian orbital fit with semi-amplitude K=3.571+/-0.041 km/s, period P=9.0090+/-0.0004 days, and eccentricity e=0.226+/-0.011. Adopting a mass of 1.16+/-0.11 Msun for the subgiant host star, we infer that the companion has a minimum mass of 40.0+/-2.5 M_Jup. Assuming an edge-on orbit, the semimajor axis is 0.090+/-0.003 AU. The host star is photometrically variable at the \sim1% level with a period of \sim13.16+/-0.01 days, indicating that the host star spin and companion orbit are not synchronized. Through adaptive optics imaging we also found a point source 643+/-10 mas away from TYC 2087-00255-1, which would have a mass of 0.13 Msun if it is physically associated with TYC 2087-00255-1 and has the same age. Future proper motion observation should be able to resolve if this tertiary object is physically associated with TYC 2087-00255-1 and make TYC 2087-00255-1 a triple body system. Core Ca II H and K line emission indicate that the host is chromospherically active, at a level that is consistent with the inferred spin period and measured v_{rot}*sin i, but unusual for a subgiant of this T_eff. This activity could be explained by ongoing tidal spin-up of the host star by the companion.

We present a detailed three-dimensional radiative-hydrodynamical simulation of the well known irradiated exoplanet HD189733b. Our model solves the fully compressible Navier-Stokes equations coupled to wavelength-dependent radiative transfer throughout the entire planetary envelope. We provide detailed comparisons between the extensive observations of this system and predictions calculated directly from the numerical models. The atmospheric dynamics is characterized by supersonic winds that fairly efficiently advect energy from the dayside to the nightside. Super-rotating equatorial jets form for a wide range of pressures from 10^-5 to 10 bars while counter rotating jets form at higher latitudes. Calculated transit spectrum agree well with the data from the infrared to the UV including the strong Rayleigh scattering seen at short wavelength, though we slightly under-predict the observations at wavelengths shorter then ~0.6 microns. Our predicted emission spectrum agrees remarkably well at 5.8 and 8 microns, but slightly over-predicts the emission at 3.6 and 4.5 microns when compared to the latest analysis by Knutson et. al (2012). Our simulated IRAC phasecurves agree fairly well with the amplitudes of variations, shape, and phases of minimum and maximum flux. However, we over-predict the peak amplitude at 3.6 and 4.5 microns, and slightly under-predict the location of the phasecurve maximum and minimum. These simulations include, for the first time in a multi-dimensional simulation, a strong Rayleigh scattering component to the absorption opacity, necessary to explain observations in the optical and UV. The agreement between our models and observations suggest that including the effects of condensates in simulations as the dominant form of opacity will be very important in future models.

We present a new algorithm for detecting transiting extrasolar planets in time-series photometry. The Quasiperiodic Automated Transit Search (QATS) algorithm relaxes the usual assumption of strictly periodic transits by permitting a variable, but bounded, interval between successive transits. We show that this method is capable of detecting transiting planets with significant transit timing variations (TTVs) without any loss of significance -- "smearing" -- as would be incurred with traditional algorithms; however, this is at the cost of an slightly-increased stochastic background. The approximate times of transit are standard products of the QATS search. Despite the increased flexibility, we show that QATS has a run-time complexity that is comparable to traditional search codes and is comparably easy to implement. QATS is applicable to data having a nearly uninterrupted, uniform cadence and is therefore well-suited to the modern class of space-based transit searches (e.g., Kepler, CoRoT). Applications of QATS include transiting planets in dynamically active multi-planet systems and transiting planets in stellar binary systems.

We derive the first secondary eclipse map of an exoplanet, HD 189733b, based on Spitzer IRAC 8 micron data. We develop two complementary techniques for deriving the two dimensional planet intensity: regularized slice mapping and spherical harmonic mapping. Both techniques give similar derived intensity maps for the infrared day-side flux of the planet, while the spherical harmonic method can be extended to include phase variation data which better constrain the map. The longitudinal offset of the day-side hot spot is consistent with that found in prior studies, strengthening the claim of super-rotating winds, and eliminating the possibility of phase variations being caused by stellar variability. The latitude of the hot-spot is within 10.1 deg (68% confidence) of the planet's equator, confirming the predictions of general circulation models for hot Jupiters and indicative of a small planet obliquity.

We report the detection of Kepler-47, a system consisting of two planets orbiting around an eclipsing pair of stars. The inner and outer planets have radii 3.0 and 4.6 times that of the Earth, respectively. The binary star consists of a Sun-like star and a companion roughly one-third its size, orbiting each other every 7.45 days. With an orbital period of 49.5 days, eighteen transits of the inner planet have been observed, allowing a detailed characterization of its orbit and those of the stars. The outer planet's orbital period is 303.2 days, and although the planet is not Earth-like, it resides within the classical "habitable zone", where liquid water could exist on an Earth-like planet. With its two known planets, Kepler-47 establishes that close binary stars can host complete planetary systems.

We discuss the discovery and characterization of the circumbinary planet Kepler-38b. The stellar binary is single-lined, with a period of 18.8 days, and consists of a moderately evolved main-sequence star (M_A = 0.949 +/- 0.059 solar masses and R_A = 1.757 +/- 0.034 solar radii) paired with a low-mass star (M_B = 0.249 +/- 0.010 solar masses and R_B = 0.2724 +/- 0.0053 solar radii) in a mildly eccentric (e=0.103) orbit. A total of eight transits due to a circumbinary planet crossing the primary star were identified in the Kepler light curve (using Kepler Quarters 1 through 11), from which a planetary period of 105.595 +/- 0.053 days can be established. A photometric dynamical model fit to the radial velocity curve and Kepler light curve yields a planetary radius of 4.35 +/- 0.11 Earth radii, or 1.12 +/- 0.03 Neptune radii. Since the planet is not sufficiently massive to observably alter the orbit of the binary from Keplerian motion, we can only place an upper limit on the mass of the planet of 122 Earth masses (7.11 Neptune masses or 0.384 Jupiter masses) at 95% confidence. This upper limit should decrease as more Kepler data become available.

We confirm 27 planets in 13 planetary systems by showing the existence of statistically significant anti-correlated transit timing variations (TTVs), which demonstrates that the planet candidates are in the same system, and long-term dynamical stability, which places limits on the masses of the candidates---showing that they are planetary. %This overall method of planet confirmation was first applied to \kepler systems 23 through 32. All of these newly confirmed planetary systems have orbital periods that place them near first-order mean motion resonances (MMRs), including 6 systems near the 2:1 MMR, 5 near 3:2, and one each near 4:3, 5:4, and 6:5. In addition, several unconfirmed planet candidates exist in some systems (that cannot be confirmed with this method at this time). A few of these candidates would also be near first order MMRs with either the confirmed planets or with other candidates. One system of particular interest, Kepler-56 (KOI-1241), is a pair of planets orbiting a 12th magnitude, giant star with radius over three times that of the Sun and effective temperature of 4900 K---among the largest stars known to host a transiting exoplanetary system.

We present new, full-orbit observations of the infrared phase variations of the canonical hot Jupiter HD 189733b obtained in the 3.6 and 4.5 micron bands using the Spitzer Space Telescope. When combined with previous phase curve observations at 8.0 and 24 micron, these data allow us to characterize the exoplanet's emission spectrum as a function of planetary longitude. We utilize improved methods for removing the effects of intrapixel sensitivity variations and accounting for the presence of time-correlated noise in our data. We measure a phase curve amplitude of 0.1242% +/- 0.0061% in the 3.6 micron band and 0.0982% +/- 0.0089% in the 4.5 micron band. We find that the times of minimum and maximum flux occur several hours earlier than predicted for an atmosphere in radiative equilibrium, consistent with the eastward advection of gas by an equatorial super-rotating jet. The locations of the flux minima in our new data differ from our previous observations at 8 micron, and we present new evidence indicating that the flux minimum observed in the 8 micron is likely caused by an over-shooting effect in the 8 micron array. We obtain improved estimates for HD 189733b's dayside planet-star flux ratio of 0.1466% +/- 0.0040% at 3.6 micron and 0.1787% +/- 0.0038% at 4.5 micron; these are the most accurate secondary eclipse depths obtained to date for an extrasolar planet. We compare our new dayside and nightside spectra for HD 189733b to the predictions of models from Burrows et al. (2008) and Showman et al. (2009). We find that HD 189733b's 4.5 micron nightside flux is 3.3 sigma smaller than predicted by the Showman et al. models, which assume that the chemistry is in local thermal equilibrium. We conclude that this discrepancy is best-explained by vertical mixing, which should lead to an excess of CO and correspondingly enhanced 4.5 micron absorption in this region. [abridged]

We report the discovery via radial velocity of a short-period (P = 2.430420 \pm 0.000006 days) companion to the F-type main sequence star TYC 2930-00872-1. A long-term trend in the radial velocities indicates the presence of a tertiary stellar companion with $P > 2000$ days. High-resolution spectroscopy of the host star yields T_eff = 6427 +/- 33 K, log(g) = 4.52 +/- 0.14, and [Fe/H]=-0.04 +/- 0.05. These parameters, combined with the broad-band spectral energy distribution and parallax, allow us to infer a mass and radius of the host star of M_1=1.21 +/- 0.08 M_\odot and R_1=1.09_{-0.13}^{+0.15} R_\odot. We are able to exclude transits of the inner companion with high confidence. The host star's spectrum exhibits clear Ca H and K core emission indicating stellar activity, but a lack of photometric variability and small v*sin(I) suggest the primary's spin axis is oriented in a pole-on configuration. The rotational period of the primary from an activity-rotation relation matches the orbital period of the inner companion to within 1.5 \sigma, suggesting they are tidally locked. If the inner companion's orbital angular momentum vector is aligned with the stellar spin axis, as expected through tidal evolution, then it has a stellar mass of M_2 ~ 0.3-0.4 M_\odot. Direct imaging limits the existence of stellar companions to projected separations < 30 AU. No set of spectral lines and no significant flux contribution to the spectral energy distribution from either companion are detected, which places individual upper mass limits of M < 1.0 M_\odot, provided they are not stellar remnants. If the tertiary is not a stellar remnant, then it likely has a mass of ~0.5-0.6 M_\odot, and its orbit is likely significantly inclined from that of the secondary, suggesting that the Kozai-Lidov mechanism may have driven the dynamical evolution of this system.

We report on the long-term dynamical evolution of the two-planet Kepler-36 system, which we studied through numerical integrations of initial conditions that are consistent with observations of the system. The orbits are chaotic with a Lyapunov time of only ~10 years. The chaos is a consequence of a particular set of orbital resonances, with the inner planet orbiting 34 times for every 29 orbits of the outer planet. The rapidity of the chaos is due to the interaction of the 29:34 resonance with the nearby first order 6:7 resonance, in contrast to the usual case in which secular terms in the Hamiltonian play a dominant role. Only one contiguous region of phase space, accounting for ~4.5% of the sample of initial conditions studied, corresponds to planetary orbits that do not show large scale orbital instabilities on the timescale of our integrations (~200 million years). The long-lived subset of the allowed initial conditions are those that satisfy the Hill stability criterion by the largest margin. Any successful theory for the formation of this system will need to account for why its current state is so close to unstable regions of phase space.

In the Solar system the planets' compositions vary with orbital distance, with rocky planets in close orbits and lower-density gas giants in wider orbits. The detection of close-in giant planets around other stars was the first clue that this pattern is not universal, and that planets' orbits can change substantially after their formation. Here we report another violation of the orbit-composition pattern: two planets orbiting the same star with orbital distances differing by only 10%, and densities differing by a factor of 8. One planet is likely a rocky `super-Earth', whereas the other is more akin to Neptune. These planets are thirty times more closely spaced--and have a larger density contrast--than any adjacent pair of planets in the Solar system.

TYC 4110-01037-1 has a low-mass stellar companion, whose small mass ratio and short orbital period are atypical amongst solar-like (Teff ~< 6000 K) binary systems. Our analysis of TYC 4110-01037-1 reveals it to be a moderately aged (~<5 Gyr) solar-like star having a mass of 1.07 +/- 0.08 MSun and radius of 0.99 +/- 0.18 RSun. We analyze 32 radial velocity measurements from the SDSS-III MARVELS survey as well as 6 supporting radial velocity measurements from the SARG spectrograph on the 3.6m TNG telescope obtained over a period of ~2 years. The best Keplerian orbital fit parameters were found to have a period of 78.994 +/- 0.012 days, an eccentricity of 0.1095 +/- 0.0023, and a semi-amplitude of 4199 +/- 11 m/s. We determine the minimum companion mass (if sin i = 1) to be 97.7 +/- 5.8 MJup. The system's companion to host star mass ratio, >0.087 +/- 0.003, places it at the lowest end of observed values for short period stellar companions to solar-like (Teff ~< 6000 K) stars. One possible way to create such a system would be if a triple-component stellar multiple broke up into a short period, low q binary during the cluster dispersal phase of its lifetime. A candidate tertiary body has been identified in the system via single-epoch, high contrast imagery. If this object is confirmed to be co-moving, we estimate it would be a dM4 star. We present these results in the context of our larger-scale effort to constrain the statistics of low mass stellar and brown dwarf companions to FGK-type stars via the MARVELS survey.

Ongoing millimeter VLBI observations with the Event Horizon Telescope allow unprecedented study of the innermost portion of black hole accretion flows. Interpreting the observations requires relativistic, time-dependent physical modeling. We discuss the comparison of radiative transfer calculations from general relativistic MHD simulations of Sagittarius A* and M87 with current and future mm-VLBI observations. This comparison allows estimates of the viewing geometry and physical conditions of the Sgr A* accretion flow. The viewing geometry for M87 is already constrained from observations of its large-scale jet, but, unlike Sgr A*, there is no consensus for its millimeter emission geometry or electron population. Despite this uncertainty, as long as the emission region is compact, robust predictions for the size of its jet launching region can be made. For both sources, the black hole shadow may be detected with future observations including ALMA and/or the LMT, which would constitute the first direct evidence for a black hole event horizon.

The supermassive black hole candidate at the center of M87 drives an ultra-relativistic jet visible on kiloparsec scales, and its large mass and relative proximity allow for event horizon scale imaging with very long baseline interferometry at millimeter wavelengths (mm-VLBI). Recently, relativistic magneto-hydrodynamic (MHD) simulations of black hole accretion flows have proven capable of launching magnetically-dominated jets. We construct time-dependent disc/jet models of the innermost portion of the M87 nucleus by performing relativistic radiative transfer calculations from one such simulation. We identify two types of models, jet-dominated or disc/jet, that can explain the spectral properties of M87, and use them to make predictions for current and future mm-VLBI observations. The Gaussian source size for the favored sky orientation and inclination from observations of the large-scale jet is 33-44 microarcseconds (~4-6 Schwarzschild radii) on current mm-VLBI telescopes, very similar to existing observations of Sgr A*. The black hole shadow, direct evidence of an event horizon, should be visible in future measurements using baselines between Hawaii and Mexico. Both models exhibit variability at millimeter wavelengths with factor of ~2 amplitudes on year timescales. For the low inclination of M87, the counter-jet dominates the event horizon scale millimeter wavelength emission from the jet-forming region.

We have analyzed Warm Spitzer/IRAC observations of the secondary eclipses of three planets, XO-4b, HAT-P-6b and HAT-P-8b. We measure secondary eclipse amplitudes at 3.6{\mu}m and 4.5{\mu}m for each target. XO-4b exhibits a stronger eclipse depth at 4.5{\mu}m than at 3.6{\mu}m, which is consistent with the presence of a temperature inversion. HAT-P-8b shows a stronger eclipse amplitude at 3.6{\mu}m, and is best-described by models without a temperature inversion. The eclipse depths of HAT-P-6b can be fitted with models with a small or no temperature inversion. We consider our results in the context of a postulated relationship between stellar activity and temperature inversions and a relationship between irradiation level and planet dayside temperature, as discussed by Knutson et al. (2010) and Cowan & Agol (2011), respectively. Our results are consistent with these hypotheses, but do not significantly strengthen them. To measure accurate secondary eclipse central phases, we require accurate ephemerides. We obtain primary transit observations and supplement them with publicly available observations to update the orbital ephemerides of the three planets. Based on the secondary eclipse timing, we set upper boundaries for e cos(\omega) for HAT-P-6b, HAT-P-8b and XO-4b and find that the values are consistent with circular orbits.

We present theoretical wavelength-dependent transit light curves for the giant planet HD209458b based on a number of state of the art 3D radiative hydrodynamical models. By varying the kinematic viscosity in the model we calculate observable signatures associated with the emergence of a super-rotating circumplanetary jet that strengthens with decreased viscosity. We obtain excellent agreement between our mid-transit transit spectra and existing data from Hubble and Spitzer, finding the best fit for intermediate values of viscosity. We further exploit dynamically driven differences between eastern and western hemispheres to extract the spectral signal imparted by a circumplanetary jet. We predict that: (i) the transit depth should decrease as the jet becomes stronger; (ii) the measured transit times should show timing offsets of up to 6 seconds at wavelengths with higher opacity, which increases with jet strength; (iii) wavelength-dependent differences between ingress and egress spectra increase with jet strength; (iv) the color-dependent transit shape should vary more strongly with wavelength for stronger jets. These techniques and trends should be valid for other hot Jupiters as well. Observations of transit timing offsets may be accessible with current instrumentation, though the other predictions may require the capabilities JWST and other future missions. Hydrodynamical models utilized solve the 3D Navier-Stokes equations together with decoupled thermal and radiative energy equations and wavelength dependent stellar heating.

Building on the legacy of the Sloan Digital Sky Survey (SDSS-I and II), SDSS-III is a program of four spectroscopic surveys on three scientific themes: dark energy and cosmological parameters, the history and structure of the Milky Way, and the population of giant planets around other stars. In keeping with SDSS tradition, SDSS-III will provide regular public releases of all its data, beginning with SDSS DR8 (which occurred in Jan 2011). This paper presents an overview of the four SDSS-III surveys. BOSS will measure redshifts of 1.5 million massive galaxies and Lya forest spectra of 150,000 quasars, using the BAO feature of large scale structure to obtain percent-level determinations of the distance scale and Hubble expansion rate at z<0.7 and at z~2.5. SEGUE-2, which is now completed, measured medium-resolution (R=1800) optical spectra of 118,000 stars in a variety of target categories, probing chemical evolution, stellar kinematics and substructure, and the mass profile of the dark matter halo from the solar neighborhood to distances of 100 kpc. APOGEE will obtain high-resolution (R~30,000), high signal-to-noise (S/N>100 per resolution element), H-band (1.51-1.70 micron) spectra of 10^5 evolved, late-type stars, measuring separate abundances for ~15 elements per star and creating the first high-precision spectroscopic survey of all Galactic stellar populations (bulge, bar, disks, halo) with a uniform set of stellar tracers and spectral diagnostics. MARVELS will monitor radial velocities of more than 8000 FGK stars with the sensitivity and cadence (10-40 m/s, ~24 visits per star) needed to detect giant planets with periods up to two years, providing an unprecedented data set for understanding the formation and dynamical evolution of giant planet systems. (Abridged)

We analyze 26 archival Kepler transits of the exo-Neptune HAT-P-11b, supplemented by ground-based transits observed in the blue (B-band) and near-IR (J-band). Both the planet and host star are smaller than previously believed; our analysis yields Rp=4.31 +/-0.06 Earth-radii, and Rs = 0.683 +/-0.009 solar radii, both about 3-sigma smaller than the discovery values. Our ground-based transit data at wavelengths bracketing the Kepler bandpass serve to check the wavelength dependence of stellar limb darkening, and the J-band transit provides a precise and independent constraint on the transit duration. Both the limb darkening and transit duration from our ground-based data are consistent with the new Kepler values for the system parameters. Our smaller radius for the planet implies that its gaseous envelope can be less extensive than previously believed, being very similar to the H-He envelope of GJ436b and Kepler-4b. HAT-P-11 is an active star, and signatures of star spot crossings are ubiquitous in the Kepler transit data. We develop and apply a methodology to correct the planetary radius for the presence of both crossed and uncrossed star spots. Star spot crossings are concentrated at phases -0.002 and +0.006. This is consistent with inferences from Rossiter-McLaughlin measurements that the planet transits nearly perpendicular to the stellar equator. We identify the dominant phases of star spot crossings with active latitudes on the star, and we infer that the stellar rotational pole is inclined at about 12 +/-5 degrees to the plane of the sky. We point out that precise transit measurements over long durations could in principle allow us to construct a stellar Butterfly diagram, to probe the cyclic evolution of magnetic activity on this active K-dwarf star.

In this paper we describe a uniform analysis of eight transits and eleven secondary eclipses of the extrasolar planet GJ 436b obtained in the 3.6, 4.5, and 8.0 micron bands using the IRAC instrument on the Spitzer Space Telescope between UT 2007 June 29 and UT 2009 Feb 4. We find that the best-fit transit depths for visits in the same bandpass can vary by as much as 8% of the total (4.7 sigma significance) from one epoch to the next. Although we cannot entirely rule out residual detector effects or a time-varying, high-altitude cloud layer in the planet's atmosphere as the cause of these variations, we consider the occultation of active regions on the star in a subset of the transit observations to be the most likely explanation. We reconcile the presence of magnetically active regions with the lack of significant visible or infrared flux variations from the star by proposing that the star's spin axis is tilted with respect to our line of sight, and that the planet's orbit is therefore likely to be misaligned. These observations serve to illustrate the challenges associated with transmission spectroscopy of planets orbiting late-type stars; we expect that other systems, such as GJ 1214, may display comparably variable transit depths. Our measured 8 micron secondary eclipse depths are consistent with a constant value, and we place a 1 sigma upper limit of 17% on changes in the planet's dayside flux in this band. Averaging over the eleven visits gives us an improved estimate of 0.0452% +/- 0.0027% for the secondary eclipse depth. We combine timing information from our observations with previously published data to produce a refined orbital ephemeris, and determine that the best-fit transit and eclipse times are consistent with a constant orbital period. [ABRIDGED]

To date the search for habitable Earth-like planets has primarily focused on nuclear burning stars. I propose that this search should be expanded to cool white dwarf stars that have expended their nuclear fuel. I define the continuously habitable zone of white dwarfs, and show that it extends from ~0.005 to 0.02 AU for white dwarfs with masses from 0.4 to 0.9 solar masses, temperatures less than 10,000 K, and habitable durations of at least 3 Gyr. As they are similar in size to Earth, white dwarfs may be deeply eclipsed by terrestrial planets that orbit edge-on, which can easily be detected with ground-based telescopes. If planets can migrate inward or reform near white dwarfs, I show that a global robotic telescope network could carry out a transit survey of nearby white dwarfs placing interesting constraints on the presence of habitable Earths. If planets were detected, I show that the survey would favor detection of planets similar to Earth: similar size, temperature, rotation period, and host star temperatures similar to the Sun. The Large Synoptic Survey Telescope could place even tighter constraints on the frequency of habitable Earths around white dwarfs. The confirmation and characterization of these planets might be carried out with large ground and space telescopes.

We aim to use microlensing taking place in the lensed quasar Q2237+0305 to study the structure of the broad line region and measure the size of the region emitting the CIV and CIII] lines. Methods: Based on 39 spectrophotometric monitoring data points obtained between Oct. 2004 and Dec. 2007, we derived lightcurves for the CIV and CIII] emission lines. We used three different techniques to analyse the microlensing signal. Different components of the lines (narrow, broad and very broad) are identified and studied. We built a library of simulated microlensing lightcurves that reproduce the signal observed in the continuum and in the lines provided only the source size is changed. A Bayesian analysis scheme is then developed to derive the size of the various components of the BLR. Results: 1. The half-light radius of the region emitting the CIV line is found to be R_CIV ~ 66^{+110}_{-46} lt-days = 0.06$^{+0.09}_{-0.04}$ pc = 1.7$^{+2.8}_{-1.1}$\,10$^{17}$ cm (at 68.3% CI). Similar values are obtained for CIII]. Relative sizes of the carbon-line and V-band continuum emitting-regions are also derived with median values of R(line)/R(cont) in the range 4 to 29, depending of the FWHM of the line component. 2. The size of the CIV emitting region agrees with the Radius-Luminosity relationship derived from reverberation mapping. Using the virial theorem we derive the mass of the black hole in Q2237+0305 to be M_BH ~ 10^{8.3+/-0.3} M_sun. 3. We find that the CIV and CIII] lines are produced in at least 2 spatially distinct regions, the most compact one giving rise to the broadest component of the line. The broad and narrow line profiles are slightly different for CIV and CIII]. 4. Our analysis suggests a different structure for the CIV and FeII+III emitting regions, with the latter produced in the inner part of the BLR or in a less extended emitting region than CIV.

We have obtained the first time-resolved, disc-integrated observations of Earth's poles with the Deep Impact spacecraft as part of the EPOXI Mission of Opportunity. These data mimic what we will see when we point next-generation space telescopes at nearby exoplanets. We use principal component analysis (PCA) and rotational lightcurve inversion to characterize color inhomogeneities and map their spatial distribution from these unusual vantage points, as a complement to the equatorial views presented in Cowan et al. (2009). We also perform the same PCA on a suite of simulated rotational multi-band lightcurves from NASA's Virtual Planetary Laboratory 3D spectral Earth model. This numerical experiment allows us to understand what sorts of surface features PCA can robustly identify. We find that the EPOXI polar observations have similar broadband colors as the equatorial Earth, but with 20-30% greater apparent albedo. This is because the polar observations are most sensitive to mid-latitudes, which tend to be more cloudy than the equatorial latitudes emphasized by the original EPOXI Earth observations. The cloudiness of the mid-latitudes also manifests itself in the form of increased variability at short wavelengths in the polar observations, and as a dominant gray eigencolor in the south polar observation. We construct a simple reflectance model for a snowball Earth. By construction, our model has a higher Bond albedo than the modern Earth; its surface albedo is so high that Rayleigh scattering does not noticeably affect its spectrum. The rotational color variations occur at short wavelengths due to the large contrast between glacier ice and bare land in those wavebands. Thus we find that both the broadband colors and diurnal color variations of such a planet would be easily distinguishable from the modern-day Earth, regardless of viewing angle.

We present an IDL graphical user interface-driven software package designed for the analysis of extrasolar planet transit light curves. The Transit Analysis Package (TAP) software uses Markov Chain Monte Carlo (MCMC) techniques to fit light curves using the analytic model of Mandel and Agol (2002). The package incorporates a wavelet based likelihood function developed by Carter and Winn (2009) which allows the MCMC to assess parameter uncertainties more robustly than classic chi-squared methods by parameterizing uncorrelated "white" and correlated "red" noise. The software is able to simultaneously analyze multiple transits observed in different conditions (instrument, filter, weather, etc). The graphical interface allows for the simple execution and interpretation of Bayesian MCMC analysis tailored to a user's specific data set and has been thoroughly tested on ground-based and Kepler photometry. AutoKep provides a similar GUI for the preparation of Kepler MAST archive data for analysis by TAP or any other analysis software. This paper describes the software release and provides instructions for its use.

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

Theoretical arguments and observations suggest that the atmospheres of Brown Dwarfs and planets are very dynamic on chemical and on physical time scales. The modelling of such substellar atmospheres has, hence, been much more demanding than initially anticipated. This Splinter (http://star-www.st-and.ac.uk/~ch80/CS16/MultiDSplinter_CS16.html) has combined new developments in atmosphere modelling, with novel observational techniques, and new challenges arising from planetary and space weather observations.

Active galactic nuclei (AGN) have been observed to vary stochastically with 10-20 rms amplitudes over a range of optical wavelengths where the emission arises in an accretion disk. Since the accretion disk is unlikely to vary coherently, local fluctuations may be significantly larger than the global rms variability. We investigate toy models of quasar accretion disks consisting of a number of regions, n, whose temperatures vary independently with an amplitude of \sigma_T in dex. Models with large fluctuations (\sigma_T=0.35-0.50) in 100-1000 independently fluctuating zones for every factor of two in radius can explain the observed discrepancy between thin accretion disk sizes inferred from microlensing events and optical luminosity while matching the observed optical variability. For the same range of \sigma_T, inhomogeneous disk spectra provide excellent fits to the HST quasar composite without invoking global Compton scattering atmospheres to explain the high levels of observed UV emission. Simulated microlensing light curves for the Einstein cross from our time-varying toy models are well fit using a time-steady power-law temperature disk, and produce magnification light curves that are consistent with current microlensing observations. Deviations due to the inhomogeneous, time-dependent disk structure should occur above the 1% level in the light curves, detectable in future microlensing observations with millimag sensitivity.

We present three transits of GJ 1214b, observed as part of the Apache Point Observatory Survey of Transit Lightcurves of Exoplanets (APOSTLE). We used APOSTLE r-band lightcurves in conjunction with previously gathered data of GJ 1214b to re-derive system parameters. By using parameters such as transit duration and ingress/egress length we are able to reduce the degeneracies between parameters in the fitted transit model, which is a preferred condition for Markov Chain Monte Carlo techniques typically used to quantify uncertainties in measured parameters. The joint analysis of this multi-wavelength dataset confirms earlier estimates of system parameters including planetary orbital period, the planet-to-star radius ratio and stellar density. We fit the photometric spectralenergy distribution of GJ 1214 to derive stellar luminosity, which we then use to derive its absolute mass and radius. From these derived stellar properties and previously published radial velocity data we were able to refine estimates of the absolute parameters for the planet GJ 1214b. Transit times derived from our study show no evidence for strong transit timing variations. Some lightcurves we present show features that we believe are due to stellar activity. During the first night we observed a rise in the out-of-eclipse flux of GJ 1214 with a characteristic fast-rise exponential decay shape commonly associated with stellar flares. On the second night we observed a minor brightening during transit, which we believe might have been caused by the planet obscuring a star-spot on the stellar disk.

We present a new short-period brown dwarf candidate around the star TYC 1240-00945-1. This candidate was discovered in the first year of the Multi-object APO Radial Velocity Exoplanets Large-area Survey (MARVELS), which is part of the third phase of the Sloan Digital Sky Survey (SDSS-III), and we designate the brown dwarf as MARVELS-1b. MARVELS uses the technique of dispersed fixed-delay interferometery to simultaneously obtain radial velocity measurements for 60 objects per field using a single, custom-built instrument that is fiber fed from the SDSS 2.5-m telescope. From our 20 radial velocity measurements spread over a ~370 d time baseline, we derive a Keplerian orbital fit with semi-amplitude K=2.533+/-0.025 km/s, period P=5.8953+/-0.0004 d, and eccentricity consistent with circular. Independent follow-up radial velocity data confirm the orbit. Adopting a mass of 1.37+/-0.11 M_Sun for the slightly evolved F9 host star, we infer that the companion has a minimum mass of 28.0+/-1.5 M_Jup, a semimajor axis 0.071+/-0.002 AU assuming an edge-on orbit, and is probably tidally synchronized. We find no evidence for coherent instrinsic variability of the host star at the period of the companion at levels greater than a few millimagnitudes. The companion has an a priori transit probability of ~14%. Although we find no evidence for transits, we cannot definitively rule them out for companion radii ~<1 R_Jup.

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

We measure secondary eclipses of the hot giant exoplanets CoRoT-1 at 3.6 and 4.5 microns, and CoRoT-2 at 3.6 microns, both using Warm Spitzer. We find that the Warm Spitzer mission is working very well for exoplanet science. For consistency of our analysis we also re-analyze archival cryogenic Spitzer data for secondary eclipses of CoRoT-2 at 4.5 and 8 microns. We compare the total data for both planets, including optical eclipse measurements by the CoRoT mission, and ground-based eclipse measurements at 2 microns, to existing models. Both planets exhibit stronger eclipses at 4.5 than at 3.6 microns, which is often indicative of an atmospheric temperature inversion. The spectrum of CoRoT-1 is best reproduced by a 2460K blackbody, due either to a high altitude layer that strongly absorbs stellar irradiance, or an isothermal region in the planetary atmosphere. The spectrum of CoRoT-2 is unusual because the 8 micron contrast is anomalously low. Non-inverted atmospheres could potentially produce the CoRoT-2 spectrum if the planet exhibits line emission from CO at 4.5 microns, caused by tidal-induced mass loss. However, the viability of that hypothesis is questionable because the emitting region cannot be more than about 30-percent larger than the planet's transit radius, based on the ingress and egress times at eclipse. An alternative possibility to account for the spectrum of CoRoT-2 is an additional opacity source that acts strongly at wavelengths less than 5 microns, heating the upper atmosphere while allowing the deeper atmosphere seen at 8 microns to remain cooler. We obtain a similar result as Gillon et al. for the phase of the secondary eclipse of CoRoT-2, implying an eccentric orbit with e*cos(omega)=-0.0030 +/- 0.0004.

We present a semi-analytic model atmosphere for close-in exoplanets that captures the essential physics of phase curves: orbital and viewing geometry, advection, and re-radiation. We calibrate the model with the well-characterized transiting planet, HD 189733b, then compute light curves for seven of the most eccentric transiting planets. We present phase variations for a variety of different radiative times and wind speeds. In the limit of instant re-radiation, the light curve morphology is entirely dictated by the planet's eccentricity and argument of pericenter: the light curve maximum leads or trails the eclipse depending on whether the planet is receding from or approaching the star at superior conjunction, respectively. For a planet with non-zero radiative timescales, the phase peak occurs early for super- rotating winds, and late for sub-rotating winds. We find that for a circular orbit, the timing of the phase variation maximum with respect to superior conjunction indicates the direction of the dominant winds, but cannot break the degeneracy between wind speed and radiative time. For circular planets the phase minimum occurs half an orbit away from the phase maximum -despite the fact that the coolest longitudes are always near the dawn terminator- and therefore does not convey any additional information. In general, increasing the advective frequency or the radiative time has the effect of reducing the peak-to-trough amplitude of phase variations, but there are interesting exceptions to these trends. Lastly, eccentric planets with orbital periods significantly longer than their radiative time exhibit "ringing" whereby the hot spot generated at periastron rotates in and out of view. The existence of ringing makes it possible to directly measure the wind speed (the frequency of the ringing) and the radiative time constant (the damping of the ringing).

Transit light curves for stellar continua have only one minimum and a "U" shape. By contrast, transit curves for optically thin chromospheric emission lines can have a "W" shape because of stellar limb-brightening. We calculate light curves for an optically thin shell of emission and fit these models to time-resolved observations of Si IV absorption by the planet HD209458b. We find that the best fit Si IV absorption model has R_p,SIV/R_*= 0.34+0.07-0.12, similar to the Roche lobe of the planet. While the large radius is only at the limit of statistical significance, we develop formulae applicable to transits of all optically thin chromospheric emission lines.

We present observations of seven transits and seven eclipses of the transiting planet system HD 189733 taken with Spitzer IRAC at 8 microns. We use a new correction for the detector ramp variation with a double-exponential function. Our main findings are: (1) an upper limit on the variability of the day-side planet flux of 2.7% (68% confidence); (2) the most precise set of transit times measured for a transiting planet, with an average accuracy of 3 seconds; (3) a lack of transit-timing variations, excluding the presence of second planets in this system above 20% of the mass of Mars in low-order mean-motion resonance at 95% confidence; (4) a confirmation of the planet's phase variation, finding the night side is 64% as bright as the day side, as well as an upper limit on the night-side variability of 17% (68% confidence); (5) a better correction for stellar variability at 8 micron causing the phase function to peak 3.5 hrs before secondary eclipse, confirming that the advection and radiation timescales are comparable at the 8 micron photosphere; (6) variation in the depth of transit, which possibly implies variations in the surface brightness of the portion of the star occulted by the planet, posing a fundamental limit on non-simultaneous multi-wavelength transit absorption measurements of planet atmospheres; (7) a measurement of the infrared limb-darkening of the star, in agreement with stellar atmosphere models; (8) an offset in the times of secondary eclipse of 69 sec, which is mostly accounted for by a 31 sec light travel time delay and 33 sec delay due to the shift of ingress and egress by the planet hot spot; this confirms that the phase variation is due to an offset hot spot on the planet; (9) a retraction of the claimed eccentricity of this system due to the offset of secondary eclipse; and (10) high precision measurements of the parameters of this system.

Recent high resolution observations of the Galactic center black hole allow for direct comparison with accretion disk simulations. We compare two-temperature synchrotron emission models from three dimensional, general relativistic magnetohydrodynamic simulations to millimeter observations of Sgr A*. Fits to very long baseline interferometry and spectral index measurements disfavor the monochromatic face-on black hole shadow models from our previous work. Inclination angles \le 20 degrees are ruled out to 3 \sigma. We estimate the inclination and position angles of the black hole, as well as the electron temperature of the accretion flow and the accretion rate, to be i=50+35-15 degrees, \xi=-23+97-22 degrees, T_e=(5.4 +/- 3.0)x10^10 K and Mdot=(5+15-2)x10^-9 M_sun / yr respectively, with 90% confidence. The black hole shadow is unobscured in all best fit models, and may be detected by observations on baselines between Chile and California, Arizona or Mexico at 1.3mm or .87mm either through direct sampling of the visibility amplitude or using closure phase information. Millimeter flaring behavior consistent with the observations is present in all viable models, and is caused by magnetic turbulence in the inner radii of the accretion flow. The variability at optically thin frequencies is strongly correlated with that in the accretion rate. The simulations provide a universal picture of the 1.3mm emission region as a small region near the midplane in the inner radii of the accretion flow, which is roughly isothermal and has \nu/\nu_c ~ 1-20, where \nu_c is the critical frequency for thermal synchrotron emission.

CoRoT, the pioneer space-based transit search, steadily provides thousands of high-precision light curves with continuous time sampling over periods of up to 5 months. The transits of a planet perturbed by an additional object are not strictly periodic. By studying the transit timing variations (TTVs), additional objects can be detected in the system. A transit timing analysis of CoRoT-1b is carried out to constrain the existence of additional planets in the system. We used data obtained by an improved version of the CoRoT data pipeline (version 2.0). Individual transits were fitted to determine the mid-transit times, and we analyzed the derived $O-C$ diagram. N-body integrations were used to place limits on secondary planets. No periodic timing variations with a period shorter than the observational window (55 days) are found. The presence of an Earth-mass Trojan is not likely. A planet of mass greater than $\sim 1$ Earth mass can be ruled out by the present data if the object is in a 2:1 (exterior) mean motion resonance with CoRoT-1b. Considering initially circular orbits: (i) super-Earths (less than 10 Earth-masses) are excluded for periods less than about 3.5 days, (ii) Saturn-like planets can be ruled out for periods less than about 5 days, (iii) Jupiter-like planets should have a minimum orbital period of about 6.5 days.

The recent VLBI observation of the Galactic center black hole candidate Sgr A* at 1.3mm shows source structure on event-horizon scales. This detection enables a direct comparison of the emission region with models of the accretion flow onto the black hole. We present the first results from time-dependent radiative transfer of general relativistic MHD simulation data, and compare simulated synchrotron images at black hole spin a=0.9 with the VLBI measurements. After tuning the accretion rate to match the millimeter flux, we find excellent agreement between predicted and observed visibilities, even when viewed face-on (i < 30 degrees). VLBI measurements on 2000-3000km baselines should constrain the inclination. The data constrain the accretion rate to be (1.0-2.3)x10^-9 M_sun / yr with 99% confidence, consistent with but independent of prior estimates derived from spectroscopic and polarimetric measurements. Finally, we compute light curves, which show that magnetic turbulence can directly produce flaring events with .5 hour rise times, 2-3.5 hour durations and 40-50% flux modulation, in agreement with observations of Sgr A* at millimeter wavelengths.

We monitor the star HD 149026 and its Saturn-mass planet at 8.0 micron over slightly more than half an orbit using the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. We find an increase of 0.0227% +/- 0.0066% (3.4 sigma significance) in the combined planet-star flux during this interval. The minimum flux from the planet is 45% +/- 19% of the maximum planet flux, corresponding to a difference in brightness temperature of 480 +/- 140 K between the two hemispheres. We derive a new secondary eclipse depth of 0.0411% +/- 0.0076% in this band, corresponding to a dayside brightness temperature of 1440 +/- 150 K. Our new secondary eclipse depth is half that of a previous measurement (3.0 sigma difference) in this same bandpass by Harrington et al. (2007). We re-fit the Harrington et al. (2007) data and obtain a comparably good fit with a smaller eclipse depth that is consistent with our new value. In contrast to earlier claims, our new eclipse depth suggests that this planet's dayside emission spectrum is relatively cool, with an 8 micron brightness temperature that is less than the maximum planet-wide equilibrium temperature. We measure the interval between the transit and secondary eclipse and find that that the secondary eclipse occurs 20.9 +7.2 / -6.5 minutes earlier (2.9 sigma) than predicted for a circular orbit, a marginally significant result. This corresponds to e*cos(omega) = -0.0079 +0.0027 / -0.0025 where e is the planet's orbital eccentricity and omega is the argument of pericenter.

We present the first resolved mid-IR (11 micron) observations of the four-image quasar lens H1413+117 using the Michelle camera on Gemini North. All previous observations (optical, near-IR, and radio) of this lens show a "flux anomaly," where the image flux ratios cannot be explained by a simple, central lens galaxy. We attempt to reproduce the mid-IR flux ratios, which are insensitive to extinction and microlensing, by modeling the main lens as a singular isothermal ellipsoid. This model fails to reproduce the flux ratios. However, we can explain the flux ratios simply by adding to the model a nearby galaxy detected in the H-band by HST/NICMOS-NIC2. This perturbing galaxy lies 4.0" from the main lens and it has a critical radius of 0.63" +/- 0.02" which is similar to that of the main lens, as expected from their similar H-band fluxes. More remarkably, this galaxy is not required to obtain a good fit to the system astrometry, so this represents the first clear detection of an object through its effect on the image fluxes of a gravitational lens. This is a parallel to the detections of visible satellites from astrometric anomalies, and provides a proof of the concept of searching for substructure in galaxies using anomalous flux ratios.

We use gravitational microlensing of the four images of the z=0.658 quasar RXJ1131-1231 to measure the sizes of the optical and X-ray emission regions of the quasar. The (face-on) scale length of the optical disk at rest frame 400 nm is 1.3 10^15cm, while the half-light radius of the rest frame 0.3-17 keV X-ray emission is 2.3 10^14cm. The formal uncertainties are factors of 1.6 and 2.0, respectively. With the exception of the lower limit on the X-ray size, the results are very stable against any changes in the priors used in the analysis. Based on the Hbeta line-width, we estimate that the black hole mass is ~10^8 Msun, which corresponds to a gravitational radius of r_g~2 10^13 cm. Thus, the X-ray emission is emerging on scales of ~10r_g and the 400 nm emission on scales of ~70 r_g. A standard thin disk of this size should be significantly brighter than observed. Possible solutions are to have a flatter temperature profile or to scatter a large fraction of the optical flux on larger scales after it is emitted. While our calculations were not optimized to constrain the dark matter fraction in the lens galaxy, dark matter dominated models are favored. With well-sampled optical and X-ray light curves over a broad range of frequencies there will be no difficulty in extending our analysis to completely map the structure of the accretion disk as a function of wavelength.

A long standing goal in astrophysics is to directly observe the immediate environment of a black hole with angular resolution comparable to the event horizon. Realizing this goal would open a new window on the study of General Relativity in the strong field regime, accretion and outflow processes at the edge of a black hole, the existence of an event horizon, and fundamental black hole physics (e.g., spin). Steady long-term progress on improving the capability of Very Long Baseline Interferometry (VLBI) at short wavelengths has now made it extremely likely that this goal will be achieved within the next decade. The most compelling evidence for this is the recent observation by 1.3mm VLBI of Schwarzschild radius scale structure in SgrA*, the compact source of radio, submm, NIR and xrays at the center of the Milky Way. SgrA* is thought to mark the position of a ~4 million solar mass black hole, and because of its proximity and estimated mass presents the largest apparent event horizon size of any black hole candidate in the Universe. Over the next decade, existing and planned mm/submm facilities will be combined into a high sensitivity, high angular resolution "Event Horizon Telescope" that will bring us as close to the edge of black hole as we will come for decades. This white paper describes the science case for mm/submm VLBI observations of both SgrA* and M87 (a radio loud AGN of a much more luminous class that SgrA*). We emphasize that while there is development and procurement involved, the technical path forward is clear, and the recent successful observations have removed much of the risk that would normally be associated with such an ambitious project.

[Abridged] To simulate the kinds of observations that will eventually be obtained for exoplanets, the Deep Impact spacecraft obtained light curves of Earth at seven wavebands spanning 300-1000 nm as part of the EPOXI mission of opportunity. In this paper we analyze disc-integrated light curves, treating Earth as if it were an exoplanet, to determine if we can detect the presence of oceans and continents. We present two observations each spanning one day, taken at gibbous phases. The rotation of the planet leads to diurnal albedo variations of 15-30%, with the largest relative changes occuring at the reddest wavelengths. To characterize these variations in an unbiased manner we carry out a principal component analysis of the multi-band light curves; this analysis reveals that 98% of the diurnal color changes of Earth are due to only 2 dominant eigencolors. We use the time-variations of these two eigencolors to construct longitudinal maps of the Earth, treating it as a non-uniform Lambert sphere. We find that the spectral and spatial distributions of the eigencolors correspond to cloud-free continents and oceans; this despite the fact that our observations were taken on days with typical cloud cover. We also find that the near-infrared wavebands are particularly useful in distinguishing between land and water. Based on this experiment we conclude that it should be possible to infer the existence of water oceans on exoplanets with time-resolved broadband observations taken by a large space-based coronagraphic telescope.

The four-image gravitationally lensed quasar QSO 2237+0305 is microlensed by stars in the lens galaxy. The amplitude of microlensing variability can be used to infer the relative size of the quasar as a function of wavelength; this provides a test of quasar models. Toward this end, we present Spitzer Space Telescope Infrared Spectrograph and Infrared Array Camera (IRAC) observations of QSO 2237+0305, finding the following. (1) The infrared (IR) spectral energy distribution (SED) is similar to that of other bright radio-quiet quasars, contrary to an earlier claim. (2) A dusty torus model with a small opening angle fits the overall shape of the IR SED well, but the quantitative agreement is poor due to an offset in wavelength of the silicate feature. (3) The flux ratios of the four lensed images can be derived from the IRAC data despite being unresolved. We find that the near-IR fluxes are increasingly affected by microlensing toward shorter wavelengths. (4) The wavelength dependence of the IRAC flux ratios is consistent with the standard quasar model in which an accretion disk and a dusty torus both contribute near 1 micron in the rest frame. This is also consistent with recent IR spectropolarimetry of nearby quasars.

Relativistic radiative transfer problems require the calculation of photon trajectories in curved spacetime. We present a novel technique for rapid and accurate calculation of null geodesics in the Kerr metric. The equations of motion from the Hamilton-Jacobi equation are reduced directly to Carlson's elliptic integrals, simplifying algebraic manipulations and allowing all coordinates to be computed semi-analytically for the first time. We discuss the method, its implementation in a freely available FORTRAN code, and its application to toy problems from the literature.

Gliese 436 is an M dwarf with a mass of 0.45 Msun and hosts the extrasolar planet GL 436b [3, 6, 7, 2], which is currently the least massive transiting planet with a mass of ~23.17 Mearth [10], and the only planet known to transit an M dwarf. GL 436b represents the first transiting detection of the class of extrasolar planets known as "Hot Neptunes" that have masses within a few times that of Neptune's mass (~17 Mearth) and orbital semimajor axis <0.1 AU about the host star. Unlike most other known transiting extrasolar planets, GL 436b has a high eccentricity (e~0.16). This brings to light a new parameter space for habitability zones of extrasolar planets with host star masses much smaller than typical stars of roughly a solar mass. This unique system is an ideal candidate for orbital perturbation and transit-time variation (TTV) studies to detect smaller, possibly Earth-mass planets in the system. In April 2008 we began a long-term intensive campaign to obtain complete high-precision light curves using the Apache Point Observatory's 3.5-meter telescope, NMSU's 1-meter telescope (located at APO), and Sommers Bausch Observatory's 24" telescope. These light curves are being analyzed together, along with amateur and other professional astronomer observations. Results of our analysis are discussed. Continued measurements over the next few years are needed to determine if additional planets reside in the system, and to study the impact of other manifestations on the light curves, such as star spots and active regions.

We present Spitzer 8 micron transit observations of the extrasolar planet HD 149026b. At this wavelength, transit light curves are weakly affected by stellar limb-darkening, allowing for a simpler and more accurate determination of planetary parameters. We measure a planet-star radius ratio of Rp/Rs = 0.05158 +/- 0.00077, and in combination with ground-based data and independent constraints on the stellar mass and radius, we derive an orbital inclination of i = 85.4 +0.9/-0.8 degrees and a planet radius of Rp = 0.755 +/- 0.040 R_jup. These measurements further support models in which the planet is greatly enriched in heavy elements.

We present the continuation of our long-term spectroscopic monitoring of the gravitationally lensed quasar QSO 2237+0305. We investigate the chromatic variations observed in the UV/optical continuum of both quasar images A and B, and compare them with numerical simulations to infer the energy profile of the quasar accretion disk. Our procedure combines the microlensing ray-shooting technique with Bayesian analysis, and derives probability distributions for the source sizes as a function of wavelength. We find that the effective caustic crossing timescale is 4.0+/-1.0 months. Using a robust prior on the effective transverse velocity, we find that the source responsible for the UV/optical continuum has an energy profile well reproduced by a power-law R lambda^{zeta} with zeta=1.2+/-0.3, where R is the source size responsible for the emission at wavelength lambda. This is the first accurate, model-independent determination of the energy profile of a quasar accretion disk on such small scales.

We present the main results of the first long-term spectrophotometric monitoring of the ``Einstein cross'' Q2237+0305 and of the single-epoch spectra of the lensed quasar J1131-1231. From October 2004 to December 2006, we find that two prominent microlensing events affect images A & B in Q2237+0305 while images C & D remain grossly unaffected by microlensing on a time scale of a few months. Microlensing in A & B goes with chromatic variations of the quasar continuum. We observe stronger micro-amplification in the blue than in the red part of the spectrum, as expected for continuum emission arising from a standard accretion disk. Microlensing induced variations of the CIII] emission are observed both in the integrated line intensity and profile. Finally, we also find that images C & D are about 0.1-0.3 mag redder than images A & B. The spectra of images A-B-C in J1131-1231 reveal that, in April 2003, microlensing was at work in images A and C. We find that microlensing de-amplifies the continuum emission and the Broad Line Region (BLR) in these images. Contrary to the case of Q2237+0305, we do not find evidence for chromatic microlensing of the continuum emission. On the other hand, we observe that the Balmer and MgII broad line profiles are deformed by microlensing. These deformations imply an anti-correlation between the width of the emission line and the size of the corresponding emitting region. Finally, the differential microlensing of the FeII emission suggests that the bulk of FeII is emitted in the outer parts of the BLR while another fraction of FeII is produced in a compact region.

We present new Spitzer observations of the phase variation of the hot Jupiter HD 189733b in the MIPS 24 micron bandpass, spanning the same part of the planet's orbit as our previous observations in the IRAC 8 micron bandpass (Knutson et al. 2007). We find that the minimum hemisphere-averaged flux from the planet in this bandpass is 76 +/- 3% of the maximum flux; this corresponds to minimum and maximum hemisphere-averaged brightness temperatures of 984 +/- 48 K and 1220 +/- 47 K, respectively. The planet reaches its maximum flux at an orbital phase of 0.396 +/- 0.022, corresponding to a hot region shifted 20-30 degrees east of the substellar point. Because tidally locked hot Jupiters would have enormous day-night temperature differences in the absence of winds, the small amplitude of the observed phase variation indicates that the planet's atmosphere efficiently transports thermal energy from the day side to the night side at the 24 micron photosphere, leading to modest day-night temperature differences. The similarities between the 8 and 24 micron phase curves for HD 189733b lead us to conclude that the circulation on this planet behaves in a fundamentally similar fashion across the range of pressures sensed by these two wavelengths. One-dimensional radiative transfer models indicate that the 8 micron band should probe pressures 2-3 times greater than at 24 micron, although the uncertain methane abundance complicates the interpretation. If these two bandpasses do probe different pressures, it would indicate that the temperature varies only weakly between the two sensed depths, and hence that the atmosphere is not convective at these altitudes. (abridged)

We present limits on transit timing variations and secondary eclipse depth variations at 8 microns with the Spitzer Space Telescope IRAC camera. Due to the weak limb darkening in the infrared and uninterrupted observing, Spitzer provides the highest accuracy transit times for this bright system, in principle providing sensitivity to secondary planets of Mars mass in resonant orbits. Finally, the transit data provides tighter constraints on the wavelength- dependent atmospheric absorption by the planet.

We present Spitzer 8 micron transit observations of the extrasolar planet system HD 149026. At this wavelength, transit light curves are weakly affected by stellar limb-darkening, allowing for a simpler and more accurate determination of planetary parameters. We measure a planet-star radius ratio of R_p/R_s = 0.05158 +/- 0.00077, and in combination with ground-based data and independent constraints on the stellar mass and radius, we derive an orbital inclination of i = 85.4 +0.9/-0.8 deg. and a planet radius of 0.755 +/- 0.040 Jupiter radii. These measurements further support models in which the planet is greatly enriched in heavy elements.

One of the most exciting results of the Spitzer era has been the ability to construct longitudinal brightness maps from the infrared phase variations of hot Jupiters. We presented the first such map in Knutson et al. (2007), described the mapping theory and some important consequences in Cowan & Agol (2008) and presented the first multi waveband map in Knutson et al. (2008). In these proceedings, we begin by putting these maps in historical context, then briefly describe the mapping formalism. We then summarize the differences between the complementary N-Slice and Sinusoidal models and end with some of the more important and surprising lessons to be learned from a careful analytic study of the mapping problem.

We present the results of the first long-term (2.2 years) spectroscopic monitoring of a gravitationally lensed quasar, namely the Einstein Cross Q2237+0305. We spatially deconvolve deep VLT/FORS1 spectra to accurately separate the spectrum of the lensing galaxy from the spectra of the quasar images. Accurate cross-calibration of the observations at 31 epochs from October 2004 to December 2006 is carried out using foreground stars observed simultaneously with the quasar. The quasar spectra are further decomposed into a continuum component and several broad emission lines. We find prominent microlensing events in the quasar images A and B, while images C and D are almost quiescent on a timescale of a few months. The strongest variations are observed in the continuum, and their amplitude is larger in the blue than in the red, consistent with microlensing of an accretion disk. Variations in the intensity and profile of the broad emission lines are also reported, most prominently in the wings of the CIII] and in the center of the CIV emission lines. During a strong microlensing episode observed in quasar image A, the broad component of the CIII] is more magnified than the narrow component. In addition, the emission lines with higher ionization potentials are more magnified than the lines with lower ionization potentials, consistent with the stratification of the broad line region (BLR) infered from reverberation mapping observations.

The periodic eclipses of the pre-main-sequence binary, KH 15D, have been explained by a circumbinary dust ring inclined to the orbital plane, which causes occultations of the stars as they pass behind the ring edge. We compute the extinction and forward scattering of light by the edge of the dust ring to explain (1) the gradual slope directly preceding total eclipse, (2) the gradual decline at the end of ingress, and (3) the slight rise in flux at mid-eclipse. The size of the forward scattering halo indicates that the dust grains have a radius of a ~ 6 (D/3 AU) microns, where D is the distance of the edge of the ring from the system barycenter. This dust size estimate agrees well with estimates of the dust grain size from polarimetry, adding to the evidence that the ring lies at several AU. Finally, the ratio of the fluxes inside and outside eclipse independently indicates that the ring lies at a few astronomical units.

We describe how to generate a longitudinal brightness map for a tidally locked exoplanet from its phase function light curve. We operate under a number of simplifying assumptions, neglecting limb darkening/brightening, star spots, detector ramps, as well as time-variability over a single planetary rotation. We develop the transformation from a planetary brightness map to a phase function light curve and simplify the expression for the case of an edge-on system. We introduce two models--composed of longitudinal slices of uniform brightness, and sinusoidally varying maps, respectively--which greatly simplify the transformation from map to light curve. We discuss numerical approaches to extracting a longitudinal map from a phase function light curve, explaining how to estimate the uncertainty in a computed map and how to choose an appropriate number of fit parameters. We demonstrate these techniques on a simulated map and discuss the uses and limitations of longitudinal maps. The sinusoidal model provides a better fit to the planet's underlying brightness map, although the slice model is more appropriate for light curves which only span a fraction of the planet's orbit. Regardless of which model is used, we find that there is a maximum of ~5 free parameters which can be meaningfully fit based on a full phase function light curve, due to the insensitivity of the latter to certain modes of the map.

We report on 2MASS J01542930+0053266, a faint eclipsing system composed of two M dwarfs. The variability of this system was originally discovered during a pilot study of the 2MASS Calibration Point Source Working Database. Additional photometry from the Sloan Digital Sky Survey yields an 8-passband lightcurve, from which we derive an orbital period of 2.6390157 +/- 0.0000016 days. Spectroscopic followup confirms our photometric classification of the system, which is likely composed of M0 and M1 dwarfs. Radial velocity measurements allow us to derive the masses (M_1 = 0.66 +/- 0.03 M_sun; M_2 = 0.62 +/- 0.03 M_sun) and radii (R_1 = 0.64 +/- 0.08 R_sun; R_2 = 0.61 +/- 0.09 R_sun) of the components, which are consistent with empirical mass-radius relationships for low-mass stars in binary systems. We perform Monte Carlo simulations of the lightcurves which allow us to uncover complicated degeneracies between the system parameters. Both stars show evidence of H-alpha emission, something not common in early-type M dwarfs. This suggests that binarity may influence the magnetic activity properties of low-mass stars; activity in the binary may persist long after the dynamos in their isolated counterparts have decayed, yielding a new potential foreground of flaring activity for next generation variability surveys.

We present the results of the first long-term (2.2 years) spectroscopic monitoring of a gravitationally lensed quasar, namely the Einstein Cross Q2237+0305. The goal of this paper is to present the observational facts to be compared in follow-up papers with theoretical models to constrain the inner structure of the source quasar. We spatially deconvolve deep VLT/FORS1 spectra to accurately separate the spectrum of the lensing galaxy from the spectra of the quasar images. Accurate cross-calibration of the 58 observations at 31-epoch from October 2004 to December 2006 is carried out with non-variable foreground stars observed simultaneously with the quasar. The quasar spectra are further decomposed into a continuum component and several broad emission lines to infer the variations of these spectral components. We find prominent microlensing events in the quasar images A and B, while images C and D are almost quiescent on a timescale of a few months. The strongest variations are observed in the continuum of image A. Their amplitude is larger in the blue (0.7 mag) than in the red (0.5 mag), consistent with microlensing of an accretion disk. Variations in the intensity and profile of the broad emission lines are also reported, most prominently in the wings of the CIII] and center of the CIV emission lines. During a strong microlensing episode observed in June 2006 in quasar image A, the broad component of the CIII] is more highly magnified than the narrow component. In addition, the emission lines with higher ionization potentials are more magnified than the lines with lower ionization potentials, consistent with the results obtained with reverberation mapping. Finally, we find that the V-band differential extinction by the lens, between the quasar images, is in the range 0.1-0.3 mag.

Precision infrared photometry from Spitzer has enabled the first direct studies of light from extrasolar planets, via observations at secondary eclipse in transiting systems. Current Spitzer results include the first longitudinal temperature map of an extrasolar planet, and the first spectra of their atmospheres. Spitzer has also measured a temperature and precise radius for the first transiting Neptune-sized exoplanet, and is beginning to make precise transit timing measurements to infer the existence of unseen low mass planets. The lack of stellar limb darkening in the infrared facilitates precise radius and transit timing measurements of transiting planets. Warm Spitzer will be capable of a precise radius measurement for Earth-sized planets transiting nearby M-dwarfs, thereby constraining their bulk composition. It will continue to measure thermal emission at secondary eclipse for transiting hot Jupiters, and be able to distinguish between planets having broad band emission versus absorption spectra. It will also be able to measure the orbital phase variation of thermal emission for close-in planets, even non-transiting planets, and these measurements will be of special interest for planets in eccentric orbits. Warm Spitzer will be a significant complement to Kepler, particularly as regards transit timing in the Kepler field. In addition to studying close-in planets, Warm Spitzer will have significant application in sensitive imaging searches for young planets at relatively large angular separations from their parent stars.

The search for habitable planets like Earth around other stars fulfils an ancient imperative to understand our origins and place in the cosmos. The past decade has seen the discovery of hundreds of planets, but nearly all are gas giants like Jupiter and Saturn. Recent advances in instrumentation and new missions are extending searches to planets the size of the Earth, but closer to their host stars. There are several possible ways such planets could form, and future observations will soon test those theories. Many of these planets we discover may be quite unlike Earth in their surface temperature and composition, but their study will nonetheless inform us about the process of planet formation and the frequency of Earth-like planets around other stars.

We present results from Spitzer Space Telescope observations of the mid-infrared phase variations of three short-period extrasolar planetary systems: HD 209458, HD 179949 and 51 Peg. We gathered IRAC images in multiple wavebands at eight phases of each planet's orbit. We find the uncertainty in relative photometry from one epoch to the next to be significantly larger than the photon counting error at 3.6 micron and 4.5 micron. We are able to place 2-sigma upper limits of only 2% on the phase variations at these wavelengths. At 8 micron the epoch-to-epoch systematic uncertainty is comparable to the photon counting noise and we detect a phase function for HD 179949 which is in phase with the planet's orbit and with a relative peak-to-trough amplitude of 0.00141(33). Assuming that HD 179949b has a radius R_J < R_p < 1.2R_J, it must recirculate less than 21% of incident stellar energy to its night side at the 1-sigma level (where 50% signifies full recirculation). If the planet has a small Bond albedo, it must have a mass less than 2.4 M_J (1-sigma). We do not detect phase variations for the other two systems but we do place the following 2-sigma upper limits: 0.0007 for 51 Peg, and 0.0015 for HD 209458. Due to its edge-on configuration, the upper limit for HD 209458 translates, with appropriate assumptions about Bond albedo, into a lower limit on the recirculation occuring in the planet's atmosphere. HD 209458b must recirculate at least 32% of incident stellar energy to its night side, at the 1-sigma level, which is consistent with other constraints on recirculation from the depth of secondary eclipse depth at 8 micron and the low optical albedo. These data indicate that different Hot Jupiter planets may experience different recirculation efficiencies.

"Hot Jupiter" extrasolar planets are expected to be tidally locked because they are close (<0.05 astronomical units, where 1 AU is the average Sun-Earth distance) to their parent stars, resulting in permanent daysides and nightsides. By observing systems where the planet and star periodically eclipse each other, several groups have been able to estimate the temperatures of the daysides of these planets. A key question is whether the atmosphere is able to transport the energy incident upon the dayside to the nightside, which will determine the temperature at different points on the planet's surface. Here we report observations of HD 189733, the closest of these eclipsing planetary systems, over half an orbital period, from which we can construct a 'map' of the distribution of temperatures. We detected the increase in brightness as the dayside of the planet rotated into view. We estimate a minimum brightness temperature of 973 +/- 33 K and a maximum brightness temperature of 1212 +/- 11 K at a wavelength of 8 microns, indicating that energy from the irradiated dayside is efficiently redistributed throughout the atmosphere, in contrast to a recent claim for another hot Jupiter. Our data indicate that the peak hemisphere-integrated brightness occurs 16$\pm$6 degrees before opposition, corresponding to a hot spot shifted east of the substellar point. The secondary eclipse (when the planet moves behind the star) occurs 120 +/- 24 s later than predicted, which may indicate a slightly eccentric orbit.

In the coming decades, research in extrasolar planets aims to advance two goals: 1) detecting and characterizing low-mass planets increasingly similar to the Earth, and 2) improving our understanding of planet formation. We present a new planet detection method that is capable of making large advances towards both of these objectives and describe a modest network of telescopes that is able to make the requisite observations. In a system where a known planet transits its host star, a second planet in that system will cause the time between transits to vary. These transit timing variations can be used to infer the orbital elements and mass of the perturbing planet even if it has a mass that is smaller than the mass of the Earth. This detection technique complements other techniques because it is most sensitive in mean-motion resonances where, due to degeneracies, other techniques have reduced sensitivity. Small ground-based observatories have already exceeded the photometric precision necessary to detect sub-Earth mass planets. However, TTV planet searches are currently limited by the relatively small number of high-precision transit data and insufficient observing time on existing telescopes. These issues will be compounded as the number of known transiting planets suitable for TTV study will increase substantially in the near future. A relatively modest investment in a ground-based network of small ($\sim 0.5 {\rm m}$ telescopes could provide the needed coverage and so dramatically increase the effectiveness of transit timing observations.

We present results from Chandra and XMM-Newton observations of the low-ionization broad absorption line (LoBAL) quasar H 1413+117. Our spatial and spectral analysis of a recent deep Chandra observation confirms a microlensing event in a previous Chandra observation performed about 5 years earlier. We present constraints on the structure of the accretion flow in H 1413+117 based on the time-scale of this microlensing event. Our analysis of the combined spectrum of all the images indicates the presence of two emission peaks at rest-frame energies of 5.35 keV and 6.32 keV detected at the > 98% and > 99% confidence levels, respectively. The double peaked Fe emission line is fit well with an accretion-disk line model, however, the best-fitting model parameters are neither well constrained nor unique. Additional observations are required to constrain the model parameters better and to confirm the relativistic interpretation of the double peaked Fe Kalpha line. Another possible interpretation of the Fe emission is fluorescent Fe emission from the back-side of the wind. The spectra of images C and D show significant high-energy broad absorption features that extend up to rest-frame energies of 9 keV and 15 keV respectively. We propose that a likely cause of these differences is significant variability of the outflow on time-scales that are shorter than the time-delays between the images. The Chandra observation of H 1413+117 has made possible for the first time the detection of the inner regions of the accretion disk and/or wind and the high ionization component of the outflowing wind of a LoBAL quasar.

In a transiting planetary system, the presence of a second planet will cause the time interval between transits to vary. These transit timing variations (TTV) are particularly large near mean-motion resonances and can be used to infer the orbital elements of planets with masses that are too small to detect by any other means. I present the results of a study of simulated data where I show the potential that this planet detection technique has to detect and characterize secondary planets in transiting systems. These results have important ramifications for planetary transit searches since each transiting system presents an opportunity for additional discoveries through a TTV analysis. I present such an analysis for 13 transits of the HD 209458 system that were observed with the Hubble Space Telescope. This analysis indicates that a putative companion in a low-order, mean-motion resonance can be no larger than the mass of the Earth and constitutes, to date, the most sensitive probe for extrasolar planets that orbit main sequence stars. The presence or absence of small planets in low-order, mean-motion resonances has implications for theories of the formation and evolution of planetary systems. Since TTV is most sensitive in these regimes, it should prove a valuable tool not only for the detection of additional planets in transiting systems, but also as a way to determine the dominant mechanisms of planet formation and the evolution of planetary systems.

I derive analytic scalings for coronagraphic imaging searches for extrasolar planets. I compute the efficiency of detecting planets about any given star, and from this compute dimensionless distribution functions for the detected planets as a function of planet-star distance and distance to the host stars. I find the following for blind planet surveys: (1) the optimum wavelength is between 4000-5000 Angstroms for Earth-like planets and 4200-5800 Angstroms for Jovian planets; (2) between 21-32% of the number of planets per decade of radius can be detected with an optimized survey; (3) target stars should be ranked from greatest to least by their luminosity divided by distance to the sixth or eighth power, depending on the dominant source of noise for the survey; (4) surveys targeting all main sequence stars will detect ~3 times as many planets as surveys only targeting G-type stars; and (5) stellar populations with different metallicities should have exposure times that vary with the cube of the metallicity. I apply these results to the current suite of proposed coronagraphic satellite telescopes, of which TPF-C is the most powerful, but a much smaller telescope, TOPS, may have a significant chance of detecting Earth-sized planets due to its small inner working angle and high throughput. The most significant uncertainty in these results is the noise contribution of Exo-zodiacal light. These results can be applied to designing coronagraphs, comparing proposed telescope designs, optimizing the observing strategies, determining the properties of detected planet populations, and selecting target stars.

We present a combined analysis of all publicly available, visible HST observations of transits of the planet HD 209458b. We derive the times of transit, planet radius, inclination, period, and ephemeris. The transit times are then used to constrain the existence of secondary planets in the system. We show that planets near an Earth mass can be ruled out in low-order mean-motion resonance, while planets less than an Earth mass are ruled out in interior, 2:1 resonance. We also present a combined analysis of the transit times and 68 high precision radial velocity measurements of the system. These results are compared to theoretical predictions for the constraints that can be placed on secondary planets.

The galaxy image pair CSL-1 has been a leading candidate for a cosmic string lens. High quality imaging data from the Hubble Space Telescope presented here show that it is not a lens but a pair of galaxies. The galaxies show different orientations of their principal axes, not consistent with any lens model. We present a new direct test of the straight-string lens model, using a displaced difference of the image from itself to exclude CSL-1 at high confidence.

The presence of a second planet in a known, transiting-planet system will cause the time between transits to vary. These variations can be used to constrain the orbital elements and mass of the perturbing planet. We analyse the set of transit times of the TrES-1 system given in Charbonneau et al. (2005). We find no convincing evidence for a second planet in the TrES-1 system from that data. By further analysis, we constrain the mass that a perturbing planet could have as a function of the semi-major axis ratio of the two planets and the eccentricity of the perturbing planet. Near low-order, mean-motion resonances (within ~1% fractional deviation), we find that a secondary planet must generally have a mass comparable to or less than the mass of the Earth--showing that this data is the first to have sensitivity to sub Earth-mass planets. We compare the sensitivity of this technique to the mass of the perturbing planet with future, high-precision radial velocity measurements.

AM CVn systems are a rare (about a dozen previously known) class of cataclysmic variables, arguably encompassing the shortest orbital periods (down to about 10 minutes) of any known binaries. Both binary components are thought to be degenerate (or partially so), likely with mass-transfer from a helium-rich donor onto a white dwarf, driven by gravitational radiation. Although rare, AM CVn systems are of high interest as possible SN Ia progenitors, and because they are predicted to be common sources of gravity waves in upcoming experiments such as LISA. We have identified four new AM CVn candidates from the Sloan Digital Sky Survey (SDSS) spectral database. All four show hallmark spectroscopic characteristics of the AM CVn class: each is devoid of hydrogen features, and instead shows a spectrum dominated by helium. All four show double-peaked emission, indicative of helium-dominated accretion disks. Limited time-series CCD photometric follow-on data have been obtained for three of the new candidates from the ARC 3.5m; most notably, a 28.3 minute binary period with sharp, deep eclipses is discovered in one case, SDSS J0926+3624. This is the first confirmed eclipsing AM CVn, and our data allow initial estimates of binary parameters for this ultracompact system. The four new SDSS objects also provide a substantial expansion of the currently critically-small sample of AM CVn systems.

The transits of a planet on a Keplerian orbit occur at time intervals exactly equal to the period of the orbit. If a second planet is introduced the orbit is not Keplerian and the transits are no longer exactly periodic. We compute the magnitude of these variations in the timing of the transits, dt. We investigate analytically several limiting cases: (i) interior perturbing planets with much smaller periods; (ii) exterior perturbing planets on eccentric orbits with much larger periods; (iii) both planets on circular orbits with arbitrary period ratio but not in resonance; and (iv) planets on initially circular orbits locked in resonance. Case (iv) is perhaps the most interesting case since some systems are known to be in resonances and the perturbations are the largest. As long as the perturber is more massive than the transiting planet, the timing variations would be of order of the period regardless of the perturber mass! For lighter perturbers, we show that the timing variations are smaller than the period by the perturber to transiting planet mass ratio. An earth mass planet in 2:1 resonance with a 3-day period transiting planet (e.g. HD 209458b) would cause timing variations of order 3 minutes, which would be accumulated over a year. These are easily detectable with current ground-based measurements. For the case of both planets on eccentric orbits, we compute numerically the transit timing variations for several cases of known multiplanet systems assuming they were edge-on. Transit timing measurements may be used to constrain the masses and radii of the planetary system and, when combined with radial velocity measurements, to break the degeneracy between mass and radius of the host star. (abstract truncated)

We present new results uncovered by a re-analysis of a Chandra observation of the gravitationally lensed, low-ionization broad absorption line (LoBAL) quasar H 1413+117. Previous analyses of the same Chandra observation led to the detection of a strong, redshifted Fe Kalpha line from the combined spectrum of all images. We show that the redshifted Fe Kalpha line is only significant in the brighter image A. The X-ray flux fraction of image A is larger by a factor of 1.55 +/- 0.17 than the optical R-band flux fraction, indicating that image A is significantly enhanced in the X-ray band. We also find that the Fe Kalpha line and the continuum are enhanced by different factors. A microlensing event could explain both the energy-dependent magnification and the significant detection of Fe Kalpha line emission in the spectrum of image A only. In the context of this interpretation we provide constraints on the spatial extent of the inferred scattered continuum and reprocessed Fe Kalpha line emission regions in a LoBAL quasar.

KH 15D is a strongly variable T Tauri star in the young star cluster NGC 2264 that shows a decrease in flux of 3.5 magnitudes lasting for 18 days and repeating every 48 days. The eclipsing material is likely due to orbiting dust or rocky bodies in a partial ring or warped disk that periodically occults the star. We measured the polarized spectrum in and out of eclipse at the Keck and Palomar observatories. Outside of the eclipse, the star exhibited low polarization consistent with zero. During eclipse, the polarization increased dramatically to ~2% across the optical spectrum, while the spectrum had the same continuum shape as outside of eclipse and exhibited emission lines of much larger equivalent width, as previously seen. From the data, we conclude that (a) the scattering region is uneclipsed; (b) the scattering is nearly achromatic; (c) the star is likely completely eclipsed so that the flux during eclipse is entirely due to scattered light, a conclusion also argued for by the shape of the ingress and egress. We argue that the scattering is not due to electrons, but may be due to large dust grains of size ~10 micron, similar to the interplanetary grains which scatter the zodiacal light. We construct a warped-disk model with an extended dusty atmosphere which reproduces the main features of the lightcurve, namely (a) a gradual decrease before ingress due to extinction in the atmosphere (similar for egress); (b) a sharper decrease within ingress due to the optically-thick base of the atmosphere; (c) a polarized flux during eclipse which is 0.1% of the total flux outside of eclipse, which requires no fine-tuning of the model. (abridged)

We make predictions for the rate of discovery of eclipsing white dwarf-main sequence (WD-MS) binaries in terrestrial-planet transit searches, taking the planned Kepler and Eddington missions as examples. We use a population synthesis model to characterize the Galactic WD-MS population, and we find that, despite increased noise due to stellar variability compared with the typical planetary case, discovery of >100 non-accreting, eclipsing WD-MS systems is likely using Kepler and Eddington, with periods of 2-20 days and transit amplitudes of |delta m|~0.0003-0.00003 magnitudes. Follow-up observations of these systems could accurately test the theoretical white dwarf mass-radius relation or theories of binary star evolution.

We prove a gravitational lensing theorem: the magnification of a source of uniform brightness by a foreground spherical lens is mu =1+pi(2R_E^2-R_L^2)/A, where A is the area of the source and R_E and R_L are the Einstein radius and size of the lens projected into the source plane; this provides an accurate approximation to the exact magnification for R_L^2,R_E^2 << A. Remarkably, this result is independent of the shape of the source or position of the lens (except near the edges). We show that this formula can be generalized to include limb-darkening of a circular source by simply inserting the surface-brightness at the position of the foreground object (divided by the average surface-brightness of the star). We also show that similar formulae apply for a point-mass lens contained in a shear field and mass sheet, and for an ensemble of point masses as long as the Einstein radii are much smaller than the source size. This theorem may be used to compute transit or microlensing lightcurves for which the foreground star or planet has a size and Einstein radius much smaller than the background star.

We present the observations of the gravitationally lensed system QSO 2237+0305 performed with the Chandra X-ray Observatory on 2000 Sept. 6, and on 2001 Dec. 8 for 30.3 ks and 9.5 ks, respectively. Imaging analysis resolves the four X-ray images of the Einstein Cross. A possible fifth image is detected; however, this detection less certain. Fits to the combined spectrum of all images of the Einstein Cross assuming a simple power law with Galactic and intervening absorption at the lensing galaxy yield a photon index of 1.90(+0.05,-0.05). For the first observation, this spectral model yields a 0.4-8.0 keV X-ray flux of 4.6e-13 erg cm-2 s-1 and a 0.4-8.0 keV lensed luminosity of 1.0e46 erg s-1. The source exhibits variability both over long and short time scales. The X-ray flux has dropped by 20% between the two observations, and the K-S test showed that image A is variable at the 97% confidence level within the first observation. Furthermore, a possible time-delay of 2.7(+0.5,-0.9) hours between images A and B with image A leading is detected in the first observation. The X-ray flux ratios of the images are consistent with the optical flux ratios which are affected by microlensing suggesting that the X-ray emission is also microlensed. A comparison between our measured column densities and those inferred from extinction measurements suggests a higher dust-to-gas ratio in the lensing galaxy than the average value of our Galaxy. Finally, we report the detection at the 99.99% confidence level of a broad emission feature near the redshifted energy of the Fe K\alpha line in only the spectrum of image A. The rest frame energy, width, and equivalent width of this feature are E = 5.7(+0.2,-0.3) keV, sigma = 0.87(+0.30,-0.15) keV, and EW = 1200(+300,-200) eV, respectively.

Occultation and microlensing are different limits of the same phenomena of one body passing in front of another body. We derive a general exact analytic expression which describes both microlensing and occultation in the case of spherical bodies with a source of uniform brightness and a non-relativistic foreground body. We also compute numerically the case of a source with quadratic limb-darkening. In the limit that the gravitational deflection angle is comparable to the angular size of the foreground body, both microlensing and occultation occur as the objects align. Such events may be used to constrain the size ratio of the lens and source stars, the limb-darkening coefficients of the source star, and the surface gravity of the lens star (if the lens and source distances are known). Application of these results to microlensing during transits in binaries and giant-star microlensing are discussed. These results unify the microlensing and occultation limits and should be useful for rapid model fitting of microlensing, eclipse, and "microccultation" events.

We present exact analytic formulae for the eclipse of a star described by quadratic or nonlinear limb darkening. In the limit that the planet radius is less than a tenth of the stellar radius, we show that the exact lightcurve can be well approximated by assuming the region of the star blocked by the planet has constant surface brightness. We apply these results to the HST observations of HD 209458, showing that the ratio of the planetary to stellar radii is 0.1207+-0.0003. These formulae give a fast and accurate means of computing lightcurves using limb-darkening coefficients from model atmospheres which should aid in the detection, simulation, and parameter fitting of planetary transits.

The MACHO and OGLE collaborations have argued that the three longest duration bulge microlensing events are likely caused by nearby black holes, given the small velocities measured with microlensing parallax and nondetection of the lenses. However, these events may be due to lensing by more numerous lower mass stars at greater distances. We find a posteriori probabilities of 76%, 16%, and 4% that the three longest events are black holes, assuming a Salpeter initial mass function (IMF) and a 40 M cutoff for neutron star progenitors; the numbers depend strongly on the assumed mass function but favor a black hole for the longest event for most standard IMFs. The longest events (>600 days) have an a priori 26% probability of being black holes for a standard mass function. We propose a new technique for measuring the lens mass function using the mass distribution of long events measured with the Advanced Camera for Surveys on the Hubble Space Telescope, the Very Large Telescope Interferometer, the Space Interferometry Mission, or the Global Astrometric Interferometer for Astrophysics.

Published parametric models of the Einstein Cross gravitational lens demonstrate that the image geometry can be reproduced by families of models. In particular, the slope of the mass-profile for the lens galaxy is unconstrained. However, recent models predict a dependence of image flux ratios on the slope of the mass profile. We use this dependence to constrain the mass profile by calculating the likelihood of the slope using published mid-IR flux ratios (including microlensing variability). We find that the galaxy is likely to be flatter than isothermal, and therefore that the mass-to-light ratio is decreasing in the inner kpc.

We use published mid-IR and V-band flux ratios for images A and B of Q2237+0305 to demonstrate that the size of the mid-IR emission region has a scale comparable to or larger than the microlens Einstein Radius (ER, ~10^17 cm for solar mass stars). Q2237+0305 has been monitored extensively in the R and V-bands for ~15 years. The variability record shows significant microlensing variability of the optical emission region, and has been used by several studies to demonstrate that the optical emission region is much smaller than the ER for solar-mass objects. For the majority of the monitoring history, the optical flux ratios have differed significantly from those predicted by macro-models. In contrast, recent observations in mid-IR show flux ratios similar to those measured in the radio, and to predictions of some lens models, implying that the mid-IR flux is emitted from a region that is at least 2 orders of magnitude larger than the optical emission region. We have calculated the likeli-hood of the observed mid-IR flux ratio as a function of mid-IR source size given the observed V-band flux ratio. The expected flux ratio for a source having dimensions of ~1 ER is a sensitive function of the macro model adopted. However we find that the probability of source size given the observed flux ratios is primarily sensitive to the ratio of the macro-model magnifications. The majority of published macro models for Q2237+0305 yield a flux ratio for images B and A of 0.8 - 1.1. By combining probabilities from the ratios A/B and C/D we infer that the diameter of a circular IR emission region is >1ER with >95% confidence. For microlensing by low-mass stars, this source size limit rules out non-thermal processes such as synchrotron as mechanisms for mid-IR emission.

We present results from monitoring of the distant (z = 2.64), gravitationally lensed quasar MG J0414+0534 with the Chandra X-ray Observatory. An Fe Ka line at 6.49 +/- 0.09 keV (rest-frame) with an equivalent width of ~ 190eV consistent with fluorescence from a cold medium is detected at the 99 percent confidence level in the spectrum of the brightest image A. During the last two observations of our monitoring program we detected a five-fold increase of the equivalent width of a narrow Fe Ka line in the spectrum of image B but not in the brighter image A whereas image C is too faint to resolve the line. The continuum emission component of image B did not follow the sudden enhancement of the iron line in the last two observations. We propose that the sudden increase in the iron line strength from ~ 190eV to 900eV only in image B can be explained with a caustic crossing due to microlensing that selectively enhances a strip of the line emission region of the accretion disk. The non-enhancement of the continuum emission in the spectrum of image B suggests that the continuum emission region is concentrated closer to the center of the accretion disk than the iron line emission region and the magnification caustic has not reached close enough to the former region to amplify it. A model of a caustic crossing event predicts discontinuities in the light-curve of the magnification and provides an upper limit of ~ 5 x 10^(-4) pc on the outer radius of the Fe Ka emission region. The non-detection of any relativistic or Doppler shifts of the iron line in the spectrum of image B implies that the magnification caustic for the last two observations was located at a radius greater than ~ 100 gravitational radii.

Galactic stellar-population-synthesis models, chemical-enrichment models, and long-duration Bulge microlensing events indicate about N_tot=(0.1-1) billion stellar-mass black holes reside in our Galaxy. We study X-ray emission from accretion from the interstellar medium on to isolated black holes. Although isolated black holes may be fewer in number than neutron stars, N_NS ~ 1 billion, their higher masses, <M>~9 M_sun, and smaller space velocities, v ~ 40 km/s, result in Bondi-Hoyle accretion rates ~4.E3 times higher than for neutron stars. We estimate that ~1.E4 N_9 isolated black holes within the Milky Way should accrete at > 1.E15 g/s, comparable to accretion rates inferred for black-hole X-ray binaries given a total number of black holes N_tot= N_9 1.E9. If black holes accrete at efficiencies only ~1.E-4(N_NS/N_tot)^0.8 of the neutron-star accretion efficiency, a comparable number of each may be detectable. We make predictions for the number of isolated accreting black holes in our Galaxy which can be detected with X-ray surveys as a function of efficiency, concluding that all-sky surveys at a depth of F = f 1.E-15 erg/cm^2/s/dex can find N(>F)~ 1.E4 N_9 (e5/f)^1.2 isolated accreting black holes for a velocity dispersion of 40 km/s and X-ray accretion efficiency = e5 1.E-5. Deeper surveys of the Galactic plane with Chandra and XMM-Newton may find tens of these objects per year, depending on the efficiency. We argue that a minimum mass can be derived for microlensing black-hole candidates if they are detected in the X-ray.

We fit general relativistic, geometrically thin accretion disk models with non-LTE atmospheres to near simultaneous multiwavelength data of 3C~273, extending from the optical to the far ultraviolet. Our model fits show no flux discontinuity associated with a hydrogen Lyman edge, but they do exhibit a spectral break which qualitatively resembles that seen in the data. This break arises from relativistic smearing of Lyman emission edges which are produced locally at tens of gravitational radii in the disk. We discuss the possible effects of metal line blanketing on the model spectra, as well as the substantial Comptonization required to explain the observed soft X-ray excess. Our best fit accretion disk model underpredicts the near ultraviolet emission in this source, and also has an optical spectrum which is too red. We discuss some of the remaining physical uncertainties, and suggest in particular that an extension of our models to the slim disk regime and/or including nonzero magnetic torques across the innermost stable circular orbit may help resolve these discrepancies.

Recent developments in the spectropolarimetric study of compact objects, specifically black holes (stellar and massive) and neutron stars are reviewed. The lectures are organized around five topics: disks, jets, outflows, neutron stars and black holes. They emphasize physical mechanisms and are intended to bridge the gap between the fundamentals of polarimetry and the phenomenology of observed cosmic sources of polarized radiation, as covered by the other lecturers. There has been considerable recent progress in spectropolarimetry from radio through optical frequencies and this is producing some unique diagnostics of the physical conditions around compact objects. It is argued that there is a great need to develop a correspondingly sensitive polarimetric capability at ultraviolet through gamma-ray energies.

We extend our models of the vertical structure and emergent radiation field of accretion disks around supermassive black holes described in previous papers of this series. Our models now include both a self-consistent treatment of Compton scattering and the effects of continuum opacities of the most important metal species (C, N, O, Ne, Mg, Si, S, Ar, Ca, Fe, Ni). With these new effects incorporated, we compute the predicted spectrum from black holes accreting at nearly the Eddington luminosity (L/L_Edd = 0.3) and central masses of 10^6, 10^7, and 10^8 M_sun. We also consider two values of the Shakura-Sunyaev alpha parameter, 0.1 and 0.01. Although it has little effect when M > 10^8 M_sun, Comptonization grows in importance as the central mass decreases and the central temperature rises. It generally produces an increase in temperature with height in the uppermost layers of hot atmospheres. Compared to models with coherent electron scattering, Comptonized models have enhanced EUV/soft X-ray emission, but they also have a more sharply declining spectrum at very high frequencies. Comptonization also smears the hydrogen and the He II Lyman edges. The effects of metals on the overall spectral energy distribution are smaller than the effects of Comptonization for these parameters. Compared to pure hydrogen-helium models, models with metal continuum opacities have reduced flux in the high frequency tail, except at the highest frequencies, where the flux is very low. Metal photoionization edges are not present in the overall disk-integrated model spectra. In addition to our new grid of models, we also present a simple analytic prescription for the vertical temperature structure of the disk in the presence of Comptonization, and show under what conditions a hot outer layer (a corona) is formed.

The standard equilibrium for radiation-dominated accretion disks has long been known to be viscously, thermally, and convectively unstable, but the nonlinear development of these instabilities---hence the actual state of such disks---has not yet been identified. By performing local two-dimensional hydrodynamic simulations of disks, we demonstrate that convective motions can release heat sufficiently rapidly as to substantially alter the vertical structure of the disk. If the dissipation rate within a vertical column is proportional to its mass, the disk settles into a new configuration thinner by a factor of two than the standard radiation-supported equilibrium. If, on the other hand, the vertically-integrated dissipation rate is proportional to the vertically-integrated total pressure, the disk is subject to the well-known thermal instability. Convection, however, biases the development of this instability toward collapse. The end result of such a collapse is a gas pressure-dominated equilibrium at the original column density.

Recent spectropolarimetry of QSOs has revealed a surprising rise in polarization in the Lyman continuum of several objects. We discuss several recent attempts to interpret this feature, including the role of Lyman limit absorption in PG 1222+228. We present new theoretical results involving electron scattering in a hot corona or wind above an accretion disk, and polarization resulting from the relativistic returning radiation. These mechanisms can lead to potentially observable polarization at short wavelengths, but neither has quantitative success in fitting the observed Lyman continuum polarization rises. A renewed capability for ultraviolet spectropolarimetry from space is urgently needed to provide further clues to the nature of this phenomenon.

Using the Long Wavelength Spectrometer on Keck I, we have imaged the gravitationally lensed radio quiet quasi-stellar object (QSO) 2237+0305 at 8.9 and 11.7 microns for the first time. The mid-infrared flux ratios are inconsistent with the optical flux ratios, but agree with the radio flux ratios and with some published gravitational lens models. These flux ratios indicate that the infrared emission is not affected by microlensing, which rules out the synchrotron emission model. The infrared emission is likely produced by hot dust extended on a length scale of more than 0.03 parsecs. The spectral energy distribution further implies a narrow range of dust temperatures, suggesting that the dust may be located in a shell extending between ~1 to 3 parsecs from the nucleus, and intercepting about half of the QSO luminosity.

The central regions of the gravitationally lensed quasar Q2237+0305 can be indirectly resolved on nano-arcsecond scales if viewed spectrophotometricly during a microlensing high magnification event (HME). Q2237+0305 is currently being monitored from the ground (eg. OGLE collaboration, Apache Point Observatory), with the goal, among others, of triggering ground and spacecraft based target of opportunity (TOO) observations of an HME. In this work we investigate the rate of change (trigger) in image brightness that signals an imminent HME and importantly, the separation between the trigger and the event peak. In addition, we produce colour dependent model light-curves by combining high-resolution microlensing simulations with a realistic model for a thermal accretion disc source. We make hypothetical target of opportunity spectroscopic observations using our determination of the appropriate trigger as a guide. We find that if the source spectrum varies with source radius, a 3 observation TOO program should be able to observe a microlensing change in the continuum slope following a light-curve trigger with a success rate of >80%.

The recent detection of polarized radiation from Sgr A* requires a non-thermal electron distribution for the emitting plasma. The Faraday rotation measure must be small, placing strong limits on the density and magnetic field strength. We show that these constraints rule out advection-dominated accretion flow models. We construct a simple two-component model which can reproduce both the radio to mm spectrum and the polarization. This model predicts that the polarization should rise to nearly 100% at shorter wavelengths. The first component, possibly a black-hole powered jet, is compact, low density, and self-absorbed near 1 mm with ordered magnetic field, relativistic Alfven speed, and a non-thermal electron distribution. The second component is poorly constrained, but may be a convection-dominated accretion flow with dM/dt~10^-9 M_Sun/yr, in which feedback from accretion onto the black hole suppresses the accretion rate at large radii. The black hole shadow should be detectable with sub-mm VLBI.

In recent years, the evidence for the existence of an ultra-compact concentration of dark mass associated with the radio source Sgr A* in the Galactic Center has become very strong. However, an unambiguous proof that this object is indeed a black hole is still lacking. A defining characteristic of a black hole is the event horizon. To a distant observer, the event horizon casts a relatively large ``shadow'' with an apparent diameter of ~10 gravitational radii due to bending of light by the black hole, nearly independent of the black hole spin or orientation. The predicted size (~30 micro-arcseconds) of this shadow for Sgr A* approaches the resolution of current radio-interferometers. If the black hole is maximally spinning and viewed edge-on, then the shadow will be offset by ~8 micro-arcseconds from the center of mass, and will be slightly flattened on one side. Taking into account scatter-broadening of the image in the interstellar medium and the finite achievable telescope resolution, we show that the shadow of Sgr A* may be observable with very long-baseline interferometry at sub-millimeter wavelengths, assuming that the accretion flow is optically thin in this region of the spectrum. Hence, there exists a realistic expectation of imaging the event horizon of a black hole within the next few years.

We have constructed a grid of non-LTE disk models for a wide range of black hole mass and mass accretion rate, for several values of viscosity parameter alpha, and for two extreme values of the black hole spin: the maximum-rotation Kerr black hole, and the Schwarzschild (non-rotating) black hole. Our procedure calculates self-consistently the vertical structure of all disk annuli together with the radiation field, without any approximations imposed on the optical thickness of the disk, and without any ad hoc approximations to the behavior of the radiation intensity. The total spectrum of a disk is computed by summing the spectra of the individual annuli, taking into account the general relativistic transfer function. The grid covers nine values of the black hole mass between M = 1/8 and 32 billion solar masses with a two-fold increase of mass for each subsequent value; and eleven values of the mass accretion rate, each a power of 2 times 1 solar mass/year. The highest value of the accretion rate corresponds to 0.3 Eddington. We show the vertical structure of individual annuli within the set of accretion disk models, along with their local emergent flux, and discuss the internal physical self-consistency of the models. We then present the full disk-integrated spectra, and discuss a number of observationally interesting properties of the models, such as optical/ultraviolet colors, the behavior of the hydrogen Lyman limit region, polarization, and number of ionizing photons. Our calculations are far from definitive in terms of the input physics, but generally we find that our models exhibit rather red optical/UV colors. Flux discontinuities in the region of the hydrogen Lyman limit are only present in cool, low luminosity models, while hotter models exhibit blueshifted changes in spectral slope.

Magnetic connections to the plunging region can exert stresses on the inner edge of an accretion disk around a black hole. We recompute the relativistic corrections to the thin-disk dynamics equations when these stresses take the form of a time-steady torque on the inner edge of the disk. The additional dissipation associated with these stresses is concentrated relatively close outside the marginally stable orbit, scaling as r to the -7/2 at large radius. As a result of these additional stresses: spin-up of the central black hole is retarded; the maximum spin-equilibrium accretion efficiency is 36%, and occurs at a/M=0.94; the disk spectrum is extended toward higher frequencies; line profiles (such as Fe K-alpha) are broadened if the line emissivity scales with local flux; limb-brightening, especially at the higher frequencies, is enhanced; and the returning radiation fraction is substantially increased, up to 58%. This last effect creates possible explanations for both synchronized continuum fluctuations in AGN, and polarization rises shortward of the Lyman edge in quasars. We show that no matter what additional stresses occur, when a/M < 0.36, the second law of black hole dynamics sets an absolute upper bound on the accretion efficiency.

We show how analysis of a quasar high-magnification microlensing event may be used to construct a map of the frequency-dependent surface brightness of the quasar accretion disk. The same procedure also allows determination of the disk inclination angle, the black hole mass (modulo the caustic velocity), and possibly the black hole spin. This method depends on the validity of one assumption: that the optical and ultraviolet continuum of the quasar is produced on the surface of an azimuthally symmetric, flat equatorial disk, whose gas follows prograde circular orbits in a Kerr spacetime (and plunges inside the marginally stable orbit). Given this assumption, we advocate using a variant of first-order linear regularization to invert multi-frequency microlensing lightcurves to obtain the disk surface brightness as a function of radius and frequency. The other parameters can be found by minimizing chi-square in a fashion consistent with the regularized solution for the surface brightness. We present simulations for a disk model appropriate to the Einstein Cross quasar, an object uniquely well-suited to this approach. These simulations confirm that the surface brightness can be reconstructed quite well near its peak, and that there are no systematic errors in determining the other model parameters. We also discuss the observational requirements for successful implementation of this technique.

MHD turbulence is generally believed to have two important functions in accretion disks: it transports angular momentum outward, and the energy in its shortest wavelength modes is dissipated into the heat that the disks radiate. In this paper we examine a pair of mechanisms which may play an important role in regulating the amplitude and spectrum of this turbulence: photon diffusion and viscosity. We demonstrate that in radiation pressure-dominated disks, photon damping of compressive MHD waves is so rapid that it likely dominates all other dissipation mechanisms.

We have calculated the spectrum and polarization of a standard thin accretion disk with parameters appropriate for a bright quasar. This model improves upon previous work by including ultraviolet metal line opacities, assumed for now to be in LTE. Though not yet fully self-consistent, our calculations demonstrate that metal lines can change the spectral slope, reduce the polarization, and reduce the Lyman edge feature in accretion disk spectra. Some observational differences between quasar spectra and accretion disk models might be reconciled with the inclusion of metal lines.

Equipartition magnetic fields can dramatically affect the polarization of radiation emerging from accretion disk atmospheres in active galactic nuclei. We extend our previous work on this subject by exploring the interaction between Faraday rotation and absorption opacity in local, plane-parallel atmospheres with parameters appropriate for accretion discs. Faraday rotation in pure scattering atmospheres acts to depolarize the radiation field by rotating the polarization planes of photons after last scattering. Absorption opacity in an unmagnetized atmosphere can increase or decrease the polarization compared to the pure scattering case, depending on the thermal source function gradient. Combining both Faraday rotation and absorption opacity, we find the following results. If absorption opacity is much larger than scattering opacity throughout the atmosphere, then Faraday rotation generally has only a small effect on the emerging polarization because of the small electron column density along a photon mean free path. However, if the absorption opacity is not too large and it acts alone to increase the polarization, then the effects of Faraday rotation can be enhanced over those in a pure scattering atmosphere. Finally, while Faraday rotation often depolarizes the radiation field, it can in some cases increase the polarization when the thermal source function does not rise too steeply with optical depth. We confirm the correctness of the Silant'ev (1979) analytic calculation of the high magnetic field limit of the pure scattering atmosphere, which we incorrectly disputed in our previous paper.