(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

A key legacy of the recently launched TESS mission will be to provide the astronomical community with many of the best transiting exoplanet targets for atmospheric characterization. However, time is of the essence to take full advantage of this opportunity. JWST, although delayed, will still complete its nominal five year mission on a timeline that motivates rapid identification, confirmation, and mass measurement of the top atmospheric characterization targets from TESS. Beyond JWST, future dedicated missions for atmospheric studies such as ARIEL require the discovery and confirmation of several hundred additional sub-Jovian size planets (R_p < 10 R_Earth) orbiting bright stars, beyond those known today, to ensure a successful statistical census of exoplanet atmospheres. Ground-based ELTs will also contribute to surveying the atmospheres of the transiting planets discovered by TESS. Here we present a set of two straightforward analytic metrics, quantifying the expected signal-to-noise in transmission and thermal emission spectroscopy for a given planet, that will allow the top atmospheric characterization targets to be readily identified among the TESS planet candidates. Targets that meet our proposed threshold values for these metrics would be encouraged for rapid follow-up and confirmation via radial velocity mass measurements. Based on the catalog of simulated TESS detections by Sullivan et al. (2015), we determine appropriate cutoff values of the metrics, such that the TESS mission will ultimately yield a sample of $\sim300$ high-quality atmospheric characterization targets across a range of planet size bins, extending down to Earth-size, potentially habitable worlds.

At a time when ALMA produces spectacular high resolution images of gas and dust in circumstellar disks, the next observational frontier in our understanding of planet formation and the chemistry of planet-forming material may be found in the mid- to far-infrared wavelength range. A large, actively cooled far-infrared telescope in space will offer enormous spectroscopic sensitivity improvements of 3-4 orders of magnitude, making it possible to uniquely survey certain fundamental properties of planet formation. Specifically, the Origins Space Telescope (OST), a NASA flagship concept to be submitted to the 2020 decadal survey, will provide a platform that allows complete surveys of warm and cold water around young stars of all masses and across all evolutionary stages, and to measure their total planet-forming gas mass using the ground-state line of HD. While this white paper is formulated in the context of the NASA Origins Space Telescope concept, it can be applied in general to inform any future space-based, cold far-infrared observatory.

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

The Large Synoptic Survey Telescope (LSST) is uniquely equipped to search for Solar System bodies due to its unprecedented combination of depth and wide field coverage. Over a ten-year period starting in 2022, LSST will generate the largest catalog of Solar System objects to date. The main goal of the LSST Solar System Science Collaboration (SSSC) is to facilitate the efforts of the planetary community to study the planets and small body populations residing within our Solar System using LSST data. To prepare for future survey cadence decisions and ensure that interesting and novel Solar System science is achievable with LSST, the SSSC has identified and prioritized key Solar System research areas for investigation with LSST in this roadmap. The ranked science priorities highlighted in this living document will inform LSST survey cadence decisions and aid in identifying software tools and pipelines needed to be developed by the planetary community as added value products and resources before the planned start of LSST science operations.

We review the evidence for buried ice in the asteroid belt; specifically the questions around the so-called Main Belt Comets (MBCs). We summarise the evidence for water throughout the Solar System, and describe the various methods for detecting it, including remote sensing from ultraviolet to radio wavelengths. We review progress in the first decade of study of MBCs, including observations, modelling of ice survival, and discussion on their origins. We then look at which methods will likely be most effective for further progress, including the key challenge of direct detection of (escaping) water in these bodies.

The Large Synoptic Survey Telescope is designed to provide an unprecedented optical imaging dataset that will support investigations of our Solar System, Galaxy and Universe, across half the sky and over ten years of repeated observation. However, exactly how the LSST observations will be taken (the observing strategy or "cadence") is not yet finalized. In this dynamically-evolving community white paper, we explore how the detailed performance of the anticipated science investigations is expected to depend on small changes to the LSST observing strategy. Using realistic simulations of the LSST schedule and observation properties, we design and compute diagnostic metrics and Figures of Merit that provide quantitative evaluations of different observing strategies, analyzing their impact on a wide range of proposed science projects. This is work in progress: we are using this white paper to communicate to each other the relative merits of the observing strategy choices that could be made, in an effort to maximize the scientific value of the survey. The investigation of some science cases leads to suggestions for new strategies that could be simulated and potentially adopted. Notably, we find motivation for exploring departures from a spatially uniform annual tiling of the sky: focusing instead on different parts of the survey area in different years in a "rolling cadence" is likely to have significant benefits for a number of time domain and moving object astronomy projects. The communal assembly of a suite of quantified and homogeneously coded metrics is the vital first step towards an automated, systematic, science-based assessment of any given cadence simulation, that will enable the scheduling of the LSST to be as well-informed as possible.

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

We report the discovery of three moderately high-mass transiting hot Jupiters from the HATSouth survey: HATS-22b, HATS-23b and HATS-24b. These planets add to the numbers of known planets in the ~2MJ regime. HATS-22b is a 2.74+/-0.11 MJ mass and 0.953+0.048/-0.029 RJ radius planet orbiting a V = 13.455 +/- 0.040 sub-silar mass (M_star = 0.759+/-0.019 M_sun; R_star = 0.759+/-0.019 R_sun) K-dwarf host star on an eccentric (e = 0.079 +/- 0.026) orbit. This planet's high planet-to-stellar mass ratio is further evidence that migration mechanisms for hot Jupiters may rely on exciting orbital eccentricities that bring planets closer to their parent stars followed by tidal circularisation. HATS-23b is a 1.478 +/- 0.080 MJ mass and 1.69 +/- 0.24 RJ radius planet on a grazing orbit around a V = 13.901 +/- 0.010 G-dwarf with properties very similar to those of the Sun (M_star = 1.115 +/- 0.054 M_sun; R_star = 1.145 +/- 0.070 R_sun). HATS-24b orbits a moderately bright V = 12.830 +/- 0.010 F-dwarf star (M_star = 1.218 +/- 0.036 M_sun; R_star = 1.194+0.066/-0.041 R_sun). This planet has a mass of 2.39 +0.21/-0.12 MJ and an inflated radius of 1.516 +0.085/-0.065 RJ.

Around the time of its perihelion passage the observability of 67P/Churyumov-Gerasimenko from Earth was limited to very short windows each morning from any given site, due to the low solar elongation of the comet. The peak in the comet's activity was therefore difficult to observe with conventionally scheduled telescopes, but was possible where service/queue scheduled mode was possible, and with robotic telescopes. We describe the robotic observations that allowed us to measure the total activity of the comet around perihelion, via photometry (dust) and spectroscopy (gas), and compare these results with the measurements at this time by Rosetta's instruments. The peak of activity occurred approximately two weeks after perihelion. The total brightness (dust) largely followed the predictions from Snodgrass et al. 2013, with no significant change in total activity levels from previous apparitions. The CN gas production rate matched previous orbits near perihelion, but appeared to be relatively low later in the year.

We present the discovery of the transiting exoplanets HAT-P-65b and HAT-P-66b, with orbital periods of 2.6055 d and 2.9721 d, masses of $0.527 \pm 0.083$ M$_{J}$ and $0.783 \pm 0.057$ M$_{J}$ and inflated radii of $1.89 \pm 0.13$ R$_{J}$ and $1.59^{+0.16}_{-0.10}$ R$_{J}$, respectively. They orbit moderately bright ($V=13.145 \pm 0.029$, and $V=12.993 \pm 0.052$) stars of mass $1.212 \pm 0.050$ M$_{\odot}$ and $1.255^{+0.107}_{-0.054}$ M$_{\odot}$. The stars are at the main sequence turnoff. While it is well known that the radii of close-in giant planets are correlated with their equilibrium temperatures, whether or not the radii of planets increase in time as their hosts evolve and become more luminous is an open question. Looking at the broader sample of well-characterized close-in transiting giant planets, we find that there is a statistically significant correlation between planetary radii and the fractional ages of their host stars, with a false alarm probability of only 0.0041%. We find that the correlation between the radii of planets and the fractional ages of their hosts is fully explained by the known correlation between planetary radii and their present day equilibrium temperatures, however if the zero-age main sequence equilibrium temperature is used in place of the present day equilibrium temperature then a correlation with age must also be included to explain the planetary radii. This suggests that, after contracting during the pre-main-sequence, close-in giant planets are re-inflated over time due to the increasing level of irradiation received from their host stars. Prior theoretical work indicates that such a dynamic response to irradiation requires a significant fraction of the incident energy to be deposited deep within the planetary interiors.

We report the discovery by the HATSouth network of HATS-18 b: a 1.980 +/- 0.077 Mj, 1.337 +0.102 -0.049 Rj planet in a 0.8378 day orbit, around a solar analog star (mass 1.037 +/- 0.047 Msun, and radius 1.020 +0.057 -0.031 Rsun) with V=14.067 +/- 0.040 mag. The high planet mass, combined with its short orbital period, implies strong tidal coupling between the planetary orbit and the star. In fact, given its inferred age, HATS-18 shows evidence of significant tidal spin up, which together with WASP-19 (a very similar system) allows us to constrain the tidal quality factor for Sun-like stars to be in the range 6.5 <= lg(Q*/k_2) <= 7 even after allowing for extremely pessimistic model uncertainties. In addition, the HATS-18 system is among the best systems (and often the best system) for testing a multitude of star--planet interactions, be they gravitational, magnetic or radiative, as well as planet formation and migration theories.

This report records and discusses the Third Workshop on Sustainable Software for Science: Practice and Experiences (WSSSPE3). The report includes a description of the keynote presentation of the workshop, which served as an overview of sustainable scientific software. It also summarizes a set of lightning talks in which speakers highlighted to-the-point lessons and challenges pertaining to sustaining scientific software. The final and main contribution of the report is a summary of the discussions, future steps, and future organization for a set of self-organized working groups on topics including developing pathways to funding scientific software; constructing useful common metrics for crediting software stakeholders; identifying principles for sustainable software engineering design; reaching out to research software organizations around the world; and building communities for software sustainability. For each group, we include a point of contact and a landing page that can be used by those who want to join that group's future activities. The main challenge left by the workshop is to see if the groups will execute these activities that they have scheduled, and how the WSSSPE community can encourage this to happen.

This technical report records and discusses the Second Workshop on Sustainable Software for Science: Practice and Experiences (WSSSPE2). The report includes a description of the alternative, experimental submission and review process, two workshop keynote presentations, a series of lightning talks, a discussion on sustainability, and five discussions from the topic areas of exploring sustainability; software development experiences; credit & incentives; reproducibility & reuse & sharing; and code testing & code review. For each topic, the report includes a list of tangible actions that were proposed and that would lead to potential change. The workshop recognized that reliance on scientific software is pervasive in all areas of world-leading research today. The workshop participants then proceeded to explore different perspectives on the concept of sustainability. Key enablers and barriers of sustainable scientific software were identified from their experiences. In addition, recommendations with new requirements such as software credit files and software prize frameworks were outlined for improving practices in sustainable software engineering. There was also broad consensus that formal training in software development or engineering was rare among the practitioners. Significant strides need to be made in building a sense of community via training in software and technical practices, on increasing their size and scope, and on better integrating them directly into graduate education programs. Finally, journals can define and publish policies to improve reproducibility, whereas reviewers can insist that authors provide sufficient information and access to data and software to allow them reproduce the results in the paper. Hence a list of criteria is compiled for journals to provide to reviewers so as to make it easier to review software submitted for publication as a "Software Paper."

The Flexible Image Transport System (FITS) standard has been a great boon to astronomy, allowing observatories, scientists and the public to exchange astronomical information easily. The FITS standard is, however, showing its age. Developed in the late 1970s the FITS authors made a number of implementation choices for the format that, while common at the time, are now seen to limit its utility with modern data. The authors of the FITS standard could not appreciate the challenges which we would be facing today in astronomical computing. Difficulties we now face include, but are not limited to, having to address the need to handle an expanded range of specialized data product types (data models), being more conducive to the networked exchange and storage of data, handling very large datasets and the need to capture significantly more complex metadata and data relationships. There are members of the community today who find some (or all) of these limitations unworkable, and have decided to move ahead with storing data in other formats. This reaction should be taken as a wakeup call to the FITS community to make changes in the FITS standard, or to see its usage fall. In this paper we detail some selected important problems which exist within the FITS standard today. It is not our intention to prescribe specific remedies to these issues; rather, we hope to call attention of the FITS and greater astronomical computing communities to these issues in the hopes that it will spur action to address them.

The Flexible Image Transport System (FITS) standard has been a great boon to astronomy, allowing observatories, scientists and the public to exchange astronomical information easily. The FITS standard, however, is showing its age. Developed in the late 1970s, the FITS authors made a number of implementation choices that, while common at the time, are now seen to limit its utility with modern data. The authors of the FITS standard could not anticipate the challenges which we are facing today in astronomical computing. Difficulties we now face include, but are not limited to, addressing the need to handle an expanded range of specialized data product types (data models), being more conducive to the networked exchange and storage of data, handling very large datasets, and capturing significantly more complex metadata and data relationships. There are members of the community today who find some or all of these limitations unworkable, and have decided to move ahead with storing data in other formats. If this fragmentation continues, we risk abandoning the advantages of broad interoperability, and ready archivability, that the FITS format provides for astronomy. In this paper we detail some selected important problems which exist within the FITS standard today. These problems may provide insight into deeper underlying issues which reside in the format and we provide a discussion of some lessons learned. It is not our intention here to prescribe specific remedies to these issues; rather, it is to call attention of the FITS and greater astronomical computing communities to these problems in the hope that it will spur action to address them.

By using the dense coverage of the extrasolar planet survey project HATNet, we Fourier analyze 381 high-probability members of the nearby open cluster Praesepe (Beehive/M44/NGC 2632). In addition to the detection of 10 variables (of \delta Scuti and other types), we identify 180 rotational variables (including the two known planet hosts). This sample increases the number of known rotational variables in this cluster for spectral classes earlier than M by more than a factor of three. These stars closely follow a color/magnitude -- period relation from early F to late K stars. We approximate this relation by polynomials for an easier reference to the rotational characteristics in different colors. The total (peak-to-peak) amplitudes of the large majority (94%) of these variables span the range of 0.005 to 0.04 mag. The periods cover a range from 2.5 to 15 days. These data strongly confirm that Praesepe and the Hyades have the same gyrochronological ages. Regarding the two planet hosts, Pr0211 (the one with the shorter orbital period) has a rotational period that is ~2 days shorter than the one expected from the main rotational pattern in this cluster. This, together with other examples discussed in the paper, may hint that star-planet interaction via tidal dissipation can be significant in some cases in the rotational evolution of stars hosting Hot Jupiters.