Showing papers in "The Astronomical Journal in 2018"

The Astropy Project: Building an Open-science Project and Status of the v2.0 Core Package

Abstract: The Astropy Project supports and fosters the development of open-source and openly developed Python packages that provide commonly needed functionality to the astronomical community. A key element of the Astropy Project is the core package astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package, as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of interoperable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy Project.

Binary Companions of Evolved Stars in APOGEE DR14: Search Method and Catalog of ∼5000 Companions

Abstract: Multi-epoch radial velocity measurements of stars can be used to identify stellar, sub-stellar, and planetary-mass companions. Even a small number of observation epochs can be informative about companions, though there can be multiple qualitatively different orbital solutions that fit the data. We have custom-built a Monte Carlo sampler (The Joker) that delivers reliable (and often highly multi-modal) posterior samplings for companion orbital parameters given sparse radial-velocity data. Here we use The Joker to perform a search for companions to 96,231 red-giant stars observed in the APOGEE survey (DR14) with $\\geq 3$ spectroscopic epochs. We select stars with probable companions by making a cut on our posterior belief about the amplitude of the stellar radial-velocity variation induced by the orbit. We provide (1) a catalog of 320 companions for which the stellar companion properties can be confidently determined, (2) a catalog of 4,898 stars that likely have companions, but would require more observations to uniquely determine the orbital properties, and (3) posterior samplings for the full orbital parameters for all stars in the parent sample. We show the characteristics of systems with confidently determined companion properties and highlight interesting systems with candidate compact object companions.

Estimating Distance from Parallaxes. IV. Distances to 1.33 Billion Stars in Gaia Data Release 2

Abstract: For the vast majority of stars in the second Gaia data release, reliable distances cannot be obtained by inverting the parallax. A correct inference procedure must instead be used to account for the nonlinearity of the transformation and the asymmetry of the resulting probability distribution. Here, we infer distances to essentially all 1.33 billion stars with parallaxes published in the second Gaia data release. This is done using a weak distance prior that varies smoothly as a function of Galactic longitude and latitude according to a Galaxy model. The irreducible uncertainty in the distance estimate is characterized by the lower and upper bounds of an asymmetric confidence interval. Although more precise distances can be estimated for a subset of the stars using additional data (such as photometry), our goal is to provide purely geometric distance estimates, independent of assumptions about the physical properties of, or interstellar extinction toward, individual stars. We analyze the characteristics of the catalog and validate it using clusters. The catalog can be queried using ADQL at http://gaia.ari.uni-heidelberg.de/tap.html (which also hosts the Gaia catalog) and downloaded from http://www.mpia.de/~calj/gdr2_distances.html.

The TESS Input Catalog and Candidate Target List

Abstract: An erratum for this article has been published in 2018 AJ 156 183, https://doi.org/10.3847/1538-3881/aade86

The California-Kepler Survey. VII. Precise Planet Radii Leveraging Gaia DR2 Reveal the Stellar Mass Dependence of the Planet Radius Gap

Abstract: The distribution of planet sizes encodes details of planet formation and evolution. We present the most precise planet size distribution to date based on Gaia parallaxes, Kepler photometry, and spectroscopic temperatures from the California-Kepler Survey. Previously, we measured stellar radii to 11% precision using high-resolution spectroscopy; by adding Gaia astrometry, the errors are now 3%. Planet radius measurements are, in turn, improved to 5% precision. With a catalog of ~1000 planets with precise properties, we probed in fine detail the gap in the planet size distribution that separates two classes of small planets, rocky super-Earths and gas-dominated sub-Neptunes. Our previous study and others suggested that the gap may be observationally under-resolved and inherently flat-bottomed, with a band of forbidden planet sizes. Analysis based on our new catalog refutes this; the gap is partially filled in. Two other important factors that sculpt the distribution are a planet's orbital distance and its host-star mass, both of which are related to a planet's X-ray/UV irradiation history. For lower-mass stars, the bimodal planet distribution shifts to smaller sizes, consistent with smaller stars producing smaller planet cores. Details of the size distribution including the extent of the "sub-Neptune desert" and the width and slope of the gap support the view that photoevaporation of low-density atmospheres is the dominant evolutionary determinant of the planet size distribution.

APOGEE Data Releases 13 and 14: Data and Analysis

Abstract: Alfred P. Sloan Foundation; U.S. Department of Energy Office of Science; Center for High-Performance Computing at the University of Utah; Brazilian Participation Group; Carnegie Institution for Science; Carnegie Mellon University; Chilean Participation Group; French Participation Group; Harvard-Smithsonian Center for Astrophysics; Instituto de Astrofisica de Canarias; Johns Hopkins University; Kavli Institute for the Physics and Mathematics of the Universe (IPMU)/University of Tokyo; Lawrence Berkeley National Laboratory; Leibniz Institut fur Astrophysik Potsdam (AIP); Max-Planck-Institut fur Astronomie (MPIA Heidelberg); Max-Planck-Institut fur Astrophysik (MPA Garching); Max-Planck-Institut fur Extra-terrestrische Physik (MPE); National Astronomical Observatories of China; New Mexico State University; New York University; University of Notre Dame; Observatorio Nacional/MCTI; Ohio State University; Pennsylvania State University; Shanghai Astronomical Observatory; United Kingdom Participation Group; Universidad Nacional Autonoma de Mexico; University of Arizona; University of Colorado Boulder; University of Oxford; University of Portsmouth; University of Utah; University of Virginia; University of Washington; University of Wisconsin; Vanderbilt University; Yale University; National Science Foundation [AST-1109178]; NSF [AST-1616636, AST-1211673]; Premium Postdoctoral Research Program of the Hungarian Academy of Sciences; Hungarian NKFI Grants of the Hungarian National Research, Development and Innovation Office [K-119517]; Birgit and Hellmuth Hertz Foundation (via the Royal Physiographic Society of Lund); Crafoord Foundation; Stiftelsen Olle Engkvist Byggmastare

Identifying Exoplanets with Deep Learning: A Five-planet Resonant Chain around Kepler-80 and an Eighth Planet around Kepler-90

Abstract: NASA's Kepler Space Telescope was designed to determine the frequency of Earth-sized planets orbiting Sun-like stars, but these planets are on the very edge of the mission's detection sensitivity. Accurately determining the occurrence rate of these planets will require automatically and accurately assessing the likelihood that individual candidates are indeed planets, even at low signal-to-noise ratios. We present a method for classifying potential planet signals using deep learning, a class of machine learning algorithms that have recently become state-of-the-art in a wide variety of tasks. We train a deep convolutional neural network to predict whether a given signal is a transiting exoplanet or a false positive caused by astrophysical or instrumental phenomena. Our model is highly effective at ranking individual candidates by the likelihood that they are indeed planets: 98.8% of the time it ranks plausible planet signals higher than false positive signals in our test set. We apply our model to a new set of candidate signals that we identified in a search of known Kepler multi-planet systems. We statistically validate two new planets that are identified with high confidence by our model. One of these planets is part of a five-planet resonant chain around Kepler-80, with an orbital period closely matching the prediction by three-body Laplace relations. The other planet orbits Kepler-90, a star which was previously known to host seven transiting planets. Our discovery of an eighth planet brings Kepler-90 into a tie with our Sun as the star known to host the most planets.

The California-Kepler Survey. IV. Metal-rich Stars Host a Greater Diversity of Planets

Abstract: Probing the connection between a star's metallicity and the presence and properties of any associated planets offers an observational link between conditions during the epoch of planet formation and mature planetary systems. We explore this connection by analyzing the metallicities of Kepler target stars and the subset of stars found to host transiting planets. After correcting for survey incompleteness, we measure planet occurrence: the number of planets per 100 stars with a given metallicity M. Planet occurrence correlates with metallicity for some, but not all, planet sizes and orbital periods. For warm super-Earths having P = 10–100 days and R_P = 1.0–1.7 R⊕, planet occurrence is nearly constant over metallicities spanning −0.4 to +0.4 dex. We find 20 warm super-Earths per 100 stars, regardless of metallicity. In contrast, the occurrence of warm sub-Neptunes (R_P= 1.7–4.0 R⊕) doubles over that same metallicity interval, from 20 to 40 planets per 100 stars. We model the distribution of planets as df ∝ 10^(βM) dM, where β characterizes the strength of any metallicity correlation. This correlation steepens with decreasing orbital period and increasing planet size. For warm super-Earths β = -0.3^(+0.2)_(-0.2), while for hot Jupiters β = +3.4^(+0.9)_(-0.8). High metallicities in protoplanetary disks may increase the mass of the largest rocky cores or the speed at which they are assembled, enhancing the production of planets larger than 1.7 R⊕. The association between high metallicity and short-period planets may reflect disk density profiles that facilitate the inward migration of solids or higher rates of planet–planet scattering.

The California-Kepler Survey. V. Peas in a Pod: Planets in a Kepler Multi-planet System Are Similar in Size and Regularly Spaced*

Abstract: We have established precise planet radii, semimajor axes, incident stellar fluxes, and stellar masses for 909 planets in 355 multi-planet systems discovered by Kepler. In this sample, we find that planets within a single multi-planet system have correlated sizes: each planet is more likely to be the size of its neighbor than a size drawn at random from the distribution of observed planet sizes. In systems with three or more planets, the planets tend to have a regular spacing: the orbital period ratios of adjacent pairs of planets are correlated. Furthermore, the orbital period ratios are smaller in systems with smaller planets, suggesting that the patterns in planet sizes and spacing are linked through formation and/or subsequent orbital dynamics. Yet, we find that essentially no planets have orbital period ratios smaller than 1.2, regardless of planet size. Using empirical mass–radius relationships, we estimate the mutual Hill separations of planet pairs. We find that 93% of the planet pairs are at least 10 mutual Hill radii apart, and that a spacing of ~20 mutual Hill radii is most common. We also find that when comparing planet sizes, the outer planet is larger in 65% ± 0.4% of cases, and the typical ratio of the outer to inner planet size is positively correlated with the temperature difference between the planets. This could be the result of photo-evaporation.

Global Climate and Atmospheric Composition of the Ultra-hot Jupiter WASP-103b from HST and Spitzer Phase Curve Observations

Abstract: We present thermal phase curve measurements for the hot Jupiter WASP-103b observed with Hubble/WFC3 and Spitzer/IRAC. The phase curves have large amplitudes and negligible hotspot offsets, indicative of poor heat redistribution to the nightside. We fit the phase variation with a range of climate maps and find that a spherical harmonics model generally provides the best fit. The phase-resolved spectra are consistent with blackbodies in the WFC3 bandpass, with brightness temperatures ranging from 1880 ± 40 K on the nightside to 2930 ± 40 K on the dayside. The dayside spectrum has a significantly higher brightness temperature in the Spitzer bands, likely due to CO emission and a thermal inversion. The inversion is not present on the nightside. We retrieved the atmospheric composition and found that it is moderately metal-enriched ($[{\rm{M}}/{\rm{H}}]={23}_{-13}^{+29}\times \mathrm{solar}$) and the carbon-to-oxygen ratio is below 0.9 at 3σ confidence. In contrast to cooler hot Jupiters, we do not detect spectral features from water, which we attribute to partial H2O dissociation. We compare the phase curves to 3D general circulation models and find that magnetic drag effects are needed to match the data. We also compare the WASP-103b spectra to brown dwarfs and young, directly imaged companions. We find that these objects have significantly larger water features, indicating that surface gravity and irradiation environment play an important role in shaping the spectra of hot Jupiters. These results highlight the 3D structure of exoplanet atmospheres and illustrate the importance of phase curve observations for understanding their complex chemistry and physics.

The APOGEE-2 Survey of the Orion Star-forming Complex. II. Six-dimensional Structure

Abstract: NSF [AST-1449476]; Research Corporation via a Time Domain Astrophysics Scialog award [24217]; UNAM-DGAPA-PAPIIT, Mexico [IN103017]; FONDECYT REGULAR project [1170476]; CONICYT PAI Concurso Nacional de Insercion en la Academia [PAI79160052]; Ministry for the Economy, Development and Tourism, Programa Iniciativa Cientica Milenio grant [IC120009]; Spanish Ministry of Economy and Competitiveness (MINECO) [AYA-2017-88254-P]; Alfred P. Sloan Foundation; U.S. Department of Energy Office of Science; Center for High-Performance Computing at the University of Utah; Brazilian Participation Group; Carnegie Institution for Science, Carnegie Mellon University; Chilean Participation Group; French Participation Group; Harvard-Smithsonian Center for Astrophysics; Instituto de Astrofisica de Canarias; Johns Hopkins University; Kavli Institute for the Physics and Mathematics of the Universe (IPMU)/University of Tokyo; Lawrence Berkeley National Laboratory; Leibniz Institut fur Astrophysik Potsdam (AIP); Max-Planck-Institut fur Astronomie (MPIA Heidelberg); Max-Planck- Institut fur Astrophysik (MPA Garching); Max-Planck-Institut fur Extraterrestrische Physik (MPE); National Astronomical Observatories of China; New Mexico State University; New York University; University of Notre Dame; Observatario Nacional/MCTI; Ohio State University; Pennsylvania State University; Shanghai Astronomical Observatory; United Kingdom Participation Group; Universidad Nacional Autonoma de Mexico; University of Arizona; University of Colorado Boulder; University of Oxford; University of Portsmouth; University of Utah; University of Virginia; University of Washington; University of Wisconsin; Vanderbilt University; Yale University

A Population Study of Gaseous Exoplanets

Abstract: We present here the analysis of 30 gaseous extrasolar planets, with temperatures between 600 and 2400 K and radii between 0.35 and 1.9 R Jup. The quality of the HST/WFC3 spatially scanned data combined with our specialized analysis tools allow us to study the largest and most self-consistent sample of exoplanetary transmission spectra to date and examine the collective behavior of warm and hot gaseous planets rather than isolated case studies. We define a new metric, the Atmospheric Detectability Index (ADI) to evaluate the statistical significance of an atmospheric detection and find statistically significant atmospheres in around 16 planets out of the 30 analyzed. For most of the Jupiters in our sample, we find the detectability of their atmospheres to be dependent on the planetary radius but not on the planetary mass. This indicates that planetary gravity plays a secondary role in the state of gaseous planetary atmospheres. We detect the presence of water vapour in all of the statistically detectable atmospheres, and we cannot rule out its presence in the atmospheres of the others. In addition, TiO and/or VO signatures are detected with 4σ confidence in WASP-76 b, and they are most likely present in WASP-121 b. We find no correlation between expected signal-to-noise and atmospheric detectability for most targets. This has important implications for future large-scale surveys.

Nbodykit: an open-source, massively parallel toolkit for large-scale structure

Abstract: Author(s): Hand, N; Feng, Y; Beutler, F; Li, Y; Modi, C; Seljak, U; Slepian, Z | Abstract: We present nbodykit, an open-source, massively parallel Python toolkit for analyzing large-scale structure (LSS) data. Using Python bindings of the Message Passing Interface, we provide parallel implementations of many commonly used algorithms in LSS. nbodykit is both an interactive and scalable piece of scientific software, performing well in a supercomputing environment while still taking advantage of the interactive tools provided by the Python ecosystem. Existing functionality includes estimators of the power spectrum, two- and three-point correlation functions, a friends-of-friends grouping algorithm, mock catalog creation via the halo occupation distribution technique, and approximate N-body simulations via the FastPM scheme. The package also provides a set of distributed data containers, insulated from the algorithms themselves, that enables nbodykit to provide a unified treatment of both simulation and observational data sets. nbodykit can be easily deployed in a high-performance computing environment, overcoming some of the traditional difficulties of using Python on supercomputers. We provide performance benchmarks illustrating the scalability of the software. The modular, component-based approach of nbodykit allows researchers to easily build complex applications using its tools. The package is extensively documented at http://nbodykit.readthedocs.io, which also includes an interactive set of example recipes for new users to explore. As open-source software, we hope nbodykit provides a common framework for the community to use and develop in confronting the analysis challenges of future LSS surveys.

The Exoplanet Population Observation Simulator. I. The Inner Edges of Planetary Systems

Abstract: National Aeronautics and Space Administration [NNX15AD94G]; NASAs Science Mission Directorate

An Update to the EVEREST K2 Pipeline: Short Cadence, Saturated Stars, and Kepler-like Photometry Down to Kp = 15

Abstract: We present an update to the EVEREST K2 pipeline that addresses various limitations in the previous version and improves the photometric precision of the light curves. We develop a fast regularization scheme for pixel-level decorrelation (PLD) and adapt the algorithm to include the PLD vectors of other 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 variable stars, although some high-frequency variables may still suffer from overfitting. On average, EVEREST 2.0 light curves have 10–20% higher photometric precision than those in version 1, yielding the highest-precision light curves at all magnitudes of any publicly available K2 catalog. For most K2 campaigns, we recover the original Kepler precision to at least = 14, and to at least = 15 for campaigns 1, 5, 6, and 13. We also detrend most short-cadence targets observed by K2, obtaining even higher photometric precision for these stars. Like all aggressive, flexible models, EVEREST is prone to overfitting, and may cause a decrease in transit depths by ~10%; we urge users to mask signals of interest using our open-source software, which we show removes this bias. Light curves for campaigns 0–8 and 10–13 are available online in the EVEREST catalog, which will be updated with future campaigns. EVEREST 2.0 is open source and is coded in a framework that can be adapted to other photometric surveys, including Kepler and the upcoming TESS mission.

Rotation of Low-mass Stars in Upper Scorpius and ρ Ophiuchus with K2

Abstract: We present an analysis of K2 light curves (LCs) for candidate members of the young Upper Sco (USco) association (~8 Myr) and the neighboring ρ Oph embedded cluster (~1 Myr). We establish ~1300 stars as probable members, ~80% of which are periodic. The phased LCs have a variety of shapes which can be attributed to physical causes ranging from stellar pulsation and stellar rotation to disk-related phenomena. We identify and discuss a number of observed behaviors. The periods are ~0.2–30 days with a peak near 2 days and the rapid period end nearing breakup velocity. M stars in the young USco region rotate systematically faster than GK stars, a pattern also present in K2 data for the older Pleiades and Praesepe systems. At higher masses (types FGK), the well-defined period–color relationship for slowly rotating stars seen in the Pleiades and Praesepe systems is not yet present in USco. Circumstellar disks are present predominantly among the more slowly rotating M stars in USco, with few disks in the subday rotators. However, M dwarfs with disks rotate faster on average than FGK systems with disks. For four of these disked M dwarfs, we provide direct evidence for disk locking based on the K2 LC morphologies. Our preliminary analysis shows a relatively mass-independent spin-up by a factor of ~3.5 between USco and the Pleiades, then mass-dependent spin-down between Pleiades and Praesepe.

275 Candidates and 149 Validated Planets Orbiting Bright Stars in K2 Campaigns 0-10

Abstract: Since 2014, NASA's K2 mission has observed large portions of the ecliptic plane in search of transiting planets and has detected hundreds of planet candidates. With observations planned until at least early 2018, K2 will continue to identify more planet candidates. We present here 275 planet candidates observed during Campaigns 0–10 of the K2 mission that are orbiting stars brighter than 13 mag (in Kepler band) and for which we have obtained high-resolution spectra (R = 44,000). These candidates are analyzed using the vespa package in order to calculate their false-positive probabilities (FPP). We find that 149 candidates are validated with an FPP lower than 0.1%, 39 of which were previously only candidates and 56 of which were previously undetected. The processes of data reduction, candidate identification, and statistical validation are described, and the demographics of the candidates and newly validated planets are explored. We show tentative evidence of a gap in the planet radius distribution of our candidate sample. Comparing our sample to the Kepler candidate sample investigated by Fulton et al., we conclude that more planets are required to quantitatively confirm the gap with K2 candidates or validated planets. This work, in addition to increasing the population of validated K2 planets by nearly 50% and providing new targets for follow-up observations, will also serve as a framework for validating candidates from upcoming K2 campaigns and the Transiting Exoplanet Survey Satellite, expected to launch in 2018.

Characterizing Earth Analogs in Reflected Light: Atmospheric Retrieval Studies for Future Space Telescopes

Abstract: Space-based high contrast imaging mission concepts for studying rocky exoplanets in reflected light are currently under community study. We develop an inverse modeling framework to estimate the science return of such missions given different instrument design considerations. By combining an exoplanet albedo model, an instrument noise model, and an ensemble Markov chain Monte Carlo sampler, we explore retrievals of atmospheric and planetary properties for Earth twins as a function of signal-to-noise ratio (SNR) and resolution ($R$). Our forward model includes Rayleigh scattering, single-layer water clouds with patchy coverage, and pressure-dependent absorption due to water vapor, oxygen, and ozone. We simulate data at $R = 70$ and $R = 140$ from 0.4-1.0 $\mu$m with SNR $ = 5, 10, 15, 20$ at 550 nm (i.e., for HabEx/LUVOIR-type instruments). At these same SNR, we simulate data for WFIRST paired with a starshade, which includes two photometric points between 0.48-0.6 $\mu$m and $R = 50$ spectroscopy from 0.6-0.97 $\mu$m. Given our noise model for WFIRST-type detectors, we find that weak detections of water vapor, ozone, and oxygen can be achieved with observations with at least $R = 70$ / SNR$\ = 15$, or $R = 140$ / SNR$\ = 10$ for improved detections. Meaningful constraints are only achieved with $R = 140$ / SNR$\ = 20$ data. The WFIRST data offer limited diagnostic information, needing at least SNR = 20 to weakly detect gases. Most scenarios place limits on planetary radius, but cannot constrain surface gravity and, thus, planetary mass.

Korea Microlensing Telescope Network Microlensing Events from 2015: Event-finding Algorithm, Vetting, and Photometry

Abstract: We present microlensing events in the 2015 Korea Microlensing Telescope Network (KMTNet) data and our procedure for identifying these events. In particular, candidates were detected with a novel "completed event" microlensing event-finder algorithm. The algorithm works by making linear fits to a (t0,teff,u0) grid of point-lens microlensing models. This approach is rendered computationally efficient by restricting u0 to just two values (0 and 1), which we show is quite adequate. The implementation presented here is specifically tailored to the commission-year character of the 2015 data, but the algorithm is quite general and has already been applied to a completely different (non-KMTNet) data set. We outline expected improvements for 2016 and future KMTNet data. The light curves of the 660 "clear microlensing" and 182 "possible microlensing" events that were found in 2015 are presented along with our policy for their public release.

The Stellar Membership of the Taurus Star-forming Region*

Abstract: The high-precision astrometry from the second data release of the Gaia mission has made it possible to greatly improve the census of members of nearby clusters and associations. I have applied the Gaia data to the Taurus star-forming region, refining the sample of known members and identifying candidates for undiscovered members. The resulting samples of members and candidates provide the best constraints to date on the distribution of ages and the initial mass function (IMF) in Taurus. Several studies over the last 30 years have proposed the existence of a population of older stars (>=10 Myr) that is associated with the Taurus clouds. The data from Gaia demonstrate that such a population does not exist. Meanwhile, previous IMF estimates for small fields surrounding the Taurus aggregates have exhibited a surplus of K7-M0 stars (0.7-0.8 Msun) relative to star-forming clusters like IC 348 and the Orion Nebula Cluster. However, that difference disappears when the new census of the entire region is considered, which should be complete for spectral types earlier than M6-M7 at A_J<1. Thus, there is little variation in the stellar IMF across the 3-4 orders of magnitude in stellar density that are present in nearby star-forming regions. Finally, I note that the proper motions of two previously known members, KPNO 15 and 2MASS J04355209+2255039, indicate that they may have been ejected from the same location within the L1536 cloud ~7200 years ago.

APOGEE Data Releases 13 and 14: Stellar Parameter and Abundance Comparisons With Independent Analyses

Abstract: Data from the SDSS-IV/Apache Point Observatory Galactic Evolution Experiment (APOGEE-2) have been released as part of SDSS Data Releases 13 (DR13) and 14 (DR14). These include high-resolution H-band spectra, radial velocities, and derived stellar parameters and abundances. DR13, released in 2016 August, contained APOGEE data for roughly 150,000 stars, and DR14, released in 2017 August, added about 110,000 more. Stellar parameters and abundances have been derived with an automated pipeline, the APOGEE Stellar Parameter and Chemical Abundance Pipeline (ASPCAP). We evaluate the performance of this pipeline by comparing the derived stellar parameters and abundances to those inferred from optical spectra and analysis for several hundred stars. For most elements-C, Na, Mg, Al, Si, S, Ca, Cr, Mn, Ni-the DR14 ASPCAP analyses have systematic differences with the comparisons samples of less than 0.05 dex (median), and random differences of less than 0.15 dex (standard deviation). These differences are a combination of the uncertainties in both the comparison samples as well as the ASPCAP analysis. Compared to the references, magnesium is the most accurate alpha-element derived by ASPCAP, and shows a very clear thin/thick disk separation, while nickel is the most accurate iron-peak element (besides iron itself). (Less)

Phase Curves of WASP-33b and HD 149026b and a New Correlation between Phase Curve Offset and Irradiation Temperature

Abstract: We present new 3.6 and 4.5 μ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 hr) 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 and 4.5 μ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. We also compare our results to predictions from general circulation models, and we find that while neither planet matches the models well, 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.

A First Catalog of Variable Stars Measured by the Asteroid Terrestrial-impact Last Alert System (ATLAS)

Abstract: The Asteroid Terrestrial-impact Last Alert System (ATLAS) carries out its primary planetary defense mission by surveying about 13000 deg^2 at least four times per night. The resulting data set is useful for the discovery of variable stars to a magnitude limit fainter than r~18, with amplitudes down to 0.01 mag for bright objects. Here we present a Data Release One catalog of variable stars based on analyzing 142 million stars measured at least 100 times in the first two years of ATLAS operations. Using a Lomb-Scargle periodogram and other variability metrics, we identify 4.7 million candidate variables which we analyze in detail. Through Space Telescope Science Institute, we publicly release lightcurves for all of them, together with a vector of 169 classification features for each star. We do this at the level of unconfirmed candidate variables in order to provide the community with a large set of homogeneously analyzed photometry and avoid pre-judging which types of objects others may find most interesting. We use machine learning to classify the candidates into fifteen different broad categories based on lightcurve morphology. About 10% (430,000 stars) pass extensive tests designed to screen out spurious variability detections: we label these as `probable' variables. Of these, 230,000 receive specific classifications as eclipsing binaries, pulsating, Mira-type, or sinusoidal variables: these are the `classified' variables. New discoveries among the probable variables number more than 300,000, while 150,000 of the classified variables are new, including about 10,000 pulsating variables, 2,000 Mira stars, and 70,000 eclipsing binaries.

An Optical Transmission Spectrum for the Ultra-hot Jupiter WASP-121b Measured with the Hubble Space Telescope

Abstract: We present an atmospheric transmission spectrum for the ultra-hot Jupiter WASP-121b, measured using the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope. Across the 0.47–1 μm wavelength range, the data imply an atmospheric opacity comparable to—and in some spectroscopic channels exceeding—that previously measured at near-infrared wavelengths (1.15–1.65 μm). Wavelength-dependent variations in the opacity rule out a gray cloud deck at a confidence level of 3.7σ and may instead be explained by VO spectral bands. We find a cloud-free model assuming chemical equilibrium for a temperature of 1500 K and a metal enrichment of 10–30× solar matches these data well. Using a free-chemistry retrieval analysis, we estimate a VO abundance of -6.6^{+0.2}/{-0.3} dex. We find no evidence for TiO and place a 3σ upper limit of −7.9 dex on its abundance, suggesting TiO may have condensed from the gas phase at the day–night limb. The opacity rises steeply at the shortest wavelengths, increasing by approximately five pressure scale heights from 0.47 to 0.3 μm n wavelength. If this feature is caused by Rayleigh scattering due to uniformly distributed aerosols, it would imply an unphysically high temperature of 6810 ± 1530 K. One alternative explanation for the short-wavelength rise is absorption due to SH (mercapto radical), which has been predicted as an important product of non-equilibrium chemistry in hot Jupiter atmospheres. Irrespective of the identity of the NUV absorber, it likely captures a significant amount of incident stellar radiation at low pressures, thus playing a significant role in the overall energy budget, thermal structure, and circulation of the atmosphere.

Measuring the Recoverability of Close Binaries in Gaia DR2 with the Robo-AO Kepler Survey

Abstract: We use the Robo-AO survey of Kepler planetary candidate host stars, the largest adaptive optics survey yet performed, to measure the recovery rate of close stellar binaries in Gaia DR2. We find that Gaia recovers binaries down to 1" at magnitude contrasts as large as 6; closer systems are not resolved, regardless of secondary brightness. Gaia DR2 binary detection does not have a strong dependence on the orientation of the stellar pairs. We find 177 nearby stars to Kepler planetary candidate host stars in Gaia DR2 that were not detected in the Robo-AO survey, almost all of which are faint (G>20); the remainder were largely targets observed by Robo-AO in poor conditions. If the primary star is the host, the impact on the radii estimates of planet candidates in these systems is likely minimal; many of these faint stars, however, could be faint eclipsing binaries that are the source of a false positive planetary transit signal. With Robo-AO and Gaia combined, we find that 18.7% of Kepler planet candidate hosts have nearby stars within 4". We also find 36 nearby stars in Gaia DR2 around 35 planetary candidate host stars detected with K2. The nearby star fraction rate for K2 planetary candidates is significantly lower than that for the primary Kepler mission. The binary recovery rate of Gaia will improve initial radius estimates of future TESS planet candidates significantly, however ground-based high-resolution follow-up observations are still needed for precise characterization and confirmation. The sensitivity of Gaia to closely separated binaries is expected to improve in later data releases.

A Catalog of Cool Dwarf Targets for the Transiting Exoplanet Survey Satellite

Abstract: P.S.M. acknowledges support from the NASA Exoplanet Research Program (XRP) under Grant No. NNX15AG08G issued through the Science Mission Directorate. B.R.-A. acknowledges the support from CONICYT PAI/Concurso Nacional Insercion en la Academia, Convocatoria 2015 79150050. This work was performed in part under contract with the Jet Propulsion Laboratory (JPL) funded by NASA through the Sagan Fellowship Program executed by the NASA Exoplanet Science Institute. This research has made use of the NASA Exoplanet Archive, which is operated by the California Institute of Technology, under contract with the National Aeronautics and Space Administration under the Exoplanet Exploration Program (Akeson et al. 2013).This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC,. https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement. (NNX15AG08G - NASA Exoplanet Research Program (XRP); 79150050 - CONICYT PAI/Concurso Nacional Insercion en la Academia, Convocatoria; NASA through the Sagan Fellowship Program; European Space Agency (ESA))

The Super Earth–Cold Jupiter Relations

Abstract: We report relations between inner ($ 1$ au) giant planets (mass $>0.3~M_{\rm J}$, or cold Jupiters) around Sun-like stars, based on data from both ground-based radial velocity (RV) observations and the Kepler mission. We find that cold Jupiters appear three times more often around hosts of super Earths than they do around field stars. Given the prevalence of the super Earth systems, their cold Jupiters can account for nearly all cold Jupiters. In other words, cold Jupiters are almost certainly ($\sim90\%$) companied by super Earths. A few corollaries follow: (1) around metal-rich ([Fe/H]$>0.1$) stars, the fraction of super Earths with cold Jupiters can rise to $60\%$ or higher; (2) the inner architecture can be strongly impacted by the outer giant and we report some observational evidence for this; (3) planetary systems like our own, with cold Jupiters but no super Earths, should be rare ($\sim1\%$). The strong correlation between super Earths and cold Jupiters establish that super Earths and cold Jupiters do not compete for solid material, rather, they share similar origins, with the cold Jupiter formation requiring a somewhat more stringent condition. Lastly, we propose a few immediate observational tests of our results, using ground-based RV observations and ongoing/planned space missions.

New Rotation Period Measurements for M Dwarfs in the Southern Hemisphere: An Abundance of Slowly Rotating, Fully Convective Stars

Abstract: Stellar rotation periods are valuable both for constraining models of angular momentum loss and for under- standing how magnetic features impact inferences of exoplanet parameters. Building on our previous work in the northern hemisphere, we have used long-term, ground-based photometric monitoring from the MEarth Observatory to measure 234 rotation periods for nearby, southern hemisphere M dwarfs. Notable examples include the exoplanet hosts GJ 1132, LHS 1140, and Proxima Centauri. We find excellent agreement between our data and K2 photometry for the overlapping subset. Amongst the sample of stars with the highest quality datasets, we recover periods in 66%; as the length of the dataset increases, our recovery rate approaches 100%. The longest rotation periods we detect are around 140 days, which we suggest represent the periods that are reached when M dwarfs are as old as the local thick disk (about 9 Gyr).

Revisiting the Phase Curves of WASP-43b: Confronting Re-analyzed Spitzer Data with Cloudy Atmospheres

Abstract: Recently acquired Hubble and Spitzer phase curves of the short-period hot Jupiter WASP-43b make it an ideal target for confronting theory with data. On the observational front, we re-analyze the 3.6 and 4.5 μm Spitzer phase curves and demonstrate that our improved analysis better removes residual red noise due to intra-pixel sensitivity, which leads to greater fluxes emanating from the nightside of WASP-43b, thus reducing the tension between theory and data. On the theoretical front, we construct cloud-free and cloudy atmospheres of WASP-43b using our Global Circulation Model (GCM), THOR, which solves the non-hydrostatic Euler equations (compared to GCMs that typically solve the hydrostatic primitive equations). The cloud-free atmosphere produces a reasonable fit to the dayside emission spectrum. The multi-phase emission spectra constrain the cloud deck to be confined to the nightside and have a finite cloud-top pressure. The multi-wavelength phase curves are naturally consistent with our cloudy atmospheres, except for the 4.5 μm phase curve, which requires the presence of enhanced carbon dioxide in the atmosphere of WASP-43b. Multi-phase emission spectra at higher spectral resolution, as may be obtained using the James Webb Space Telescope, and a reflected-light phase curve at visible wavelengths would further constrain the properties of clouds in WASP-43b.

An HST/WFC3 thermal emission spectrum of the hot Jupiter HAT-P-7b

Abstract: Secondary eclipse observations of several of the hottest hot Jupiters show featureless, blackbody-like spectra or molecular emission features, which are consistent with thermal inversions being present in those atmospheres. Theory predicts a transition between warmer atmospheres with thermal inversions and cooler atmospheres without inversions, but the exact transition point is unknown. In order to further investigate this issue, we observed two secondary eclipses of the hot Jupiter HAT-P-7b with the Hubble Space Telescope (HST) WFC3 instrument and combined these data with previous Spitzer and Kepler secondary eclipse observations. The HST and Spitzer data can be well fit by a blackbody with T = 2692 ± 14 K, and the Kepler data point constrains the geometric albedo to A g = 0.077 ± 0.006. We modeled these data with a three-dimensional (3D) GCM and one-dimensional (1D) self-consistent forward models. The 1D models indicate that the atmosphere has a thermal inversion, weak heat redistribution, and water dissociation that limits the range of pressures probed. This result suggests that WFC3 observations of HAT-P-7b and possibly some other ultra-hot Jupiters appear blackbody-like because they probe a region near the tropopause where the atmospheric temperature changes slowly with pressure. Additionally, the 1D models constrain the atmospheric metallicity ($[{\rm{M}}/{\rm{H}}]=-{0.87}_{-0.34}^{+0.38}$) and the carbon-to-oxygen ratio (C/O ≺ 1 at 99% confidence). The solar composition 3D GCM matches the Spitzer data but generally underpredicts the flux in the WFC3 bandpass and cannot reproduce its featureless shape. This discrepancy could be explained by high atmospheric drag or nightside clouds and may be better understood through further observation with the James Webb Space Telescope.