Multi-messenger astrophysics is a fast-growing, interdisciplinary field that combines data, which vary in volume and speed of data processing, from many different instruments that probe the Universe using different cosmic messengers: electromagnetic waves, cosmic rays, gravitational waves and neutrinos. In this Expert Recommendation, we review the key challenges of real-time observations of gravitational wave sources and their electromagnetic and astroparticle counterparts, and make a number of recommendations to maximize their potential for scientific discovery. These recommendations refer to the design of scalable and computationally efficient machine learning algorithms; the cyber-infrastructure to numerically simulate astrophysical sources, and to process and interpret multi-messenger astrophysics data; the management of gravitational wave detections to trigger real-time alerts for electromagnetic and astroparticle follow-ups; a vision to harness future developments of machine learning and cyber-infrastructure resources to cope with the big-data requirements; and the need to build a community of experts to realize the goals of multi-messenger astrophysics.

The field of astronomy has arrived at a turning point in terms of size and complexity of both datasets and scientific collaboration. Commensurately, algorithms and statistical models have begun to adapt --- e.g., via the onset of artificial intelligence --- which itself presents new challenges and opportunities for growth. This white paper aims to offer guidance and ideas for how we can evolve our technical and collaborative frameworks to promote efficient algorithmic development and take advantage of opportunities for scientific discovery in the petabyte era. We discuss challenges for discovery in large and complex data sets; challenges and requirements for the next stage of development of statistical methodologies and algorithmic tool sets; how we might change our paradigms of collaboration and education; and the ethical implications of scientists' contributions to widely applicable algorithms and computational modeling. We start with six distinct recommendations that are supported by the commentary following them. This white paper is related to a larger corpus of effort that has taken place within and around the Petabytes to Science Workshops (https://petabytestoscience.github.io/).

Software is a critical part of modern research, and yet there are insufficient mechanisms in the scholarly ecosystem to acknowledge, cite, and measure the impact of research software. The majority of academic fields rely on a one-dimensional credit model whereby academic articles (and their associated citations) are the dominant factor in the success of a researcher's career. In the petabyte era of astronomical science, citing software and measuring its impact enables academia to retain and reward researchers that make significant software contributions. These highly skilled researchers must be retained to maximize the scientific return from petabyte-scale datasets. Evolving beyond the one-dimensional credit model requires overcoming several key challenges, including the current scholarly ecosystem and scientific culture issues. This white paper will present these challenges and suggest practical solutions for elevating the role of software as a product of the research enterprise.

We search for transiting exoplanets around the star $\beta$ Pictoris using high resolution spectroscopy and Doppler imaging that removes the need for standard star observations. These data were obtained on the VLT with UVES during the course of an observing campaign throughout 2017 that monitored the Hill sphere transit of the exoplanet $\beta$ Pictoris b. We utilize line profile tomography as a method for the discovery of transiting exoplanets. By measuring the exoplanet distortion of the stellar line profile, we remove the need for reference star measurements. We demonstrate the method with white noise simulations, and then look at the case of $\beta$ Pictoris, which is a $\delta$ Scuti pulsator. We describe a method to remove the stellar pulsations and perform a search for any transiting exoplanets in the resultant data set. We inject fake planet transits with varying orbital periods and planet radii into the spectra and determine the recovery fraction. In the photon noise limited case we can recover planets down to a Neptune radius with an $\sim$80% success rate, using an 8 m telescope with a $R\sim 100,000$ spectrograph and 20 minutes of observations per night. The pulsations of $\beta$ Pictoris limit our sensitivity to Jupiter-sized planets, but a pulsation removal algorithm improves this limit to Saturn-sized planets. We present two planet candidates, but argue that their signals are most likely caused by other phenomena. We have demonstrated a method for searching for transiting exoplanets that (i) does not require ancillary calibration observations, (ii) can work on any star whose rotational broadening can be resolved with a high spectral dispersion spectrograph and (iii) provides the lowest limits so far on the radii of transiting Jupiter-sized exoplanets around $\beta$ Pictoris with orbital periods from 15 days to 200 days with >50% coverage.

This report provides an overview of recent work that harnesses the Big Data Revolution and Large Scale Computing to address grand computational challenges in Multi-Messenger Astrophysics, with a particular emphasis on real-time discovery campaigns. Acknowledging the transdisciplinary nature of Multi-Messenger Astrophysics, this document has been prepared by members of the physics, astronomy, computer science, data science, software and cyberinfrastructure communities who attended the NSF-, DOE- and NVIDIA-funded "Deep Learning for Multi-Messenger Astrophysics: Real-time Discovery at Scale" workshop, hosted at the National Center for Supercomputing Applications, October 17-19, 2018. Highlights of this report include unanimous agreement that it is critical to accelerate the development and deployment of novel, signal-processing algorithms that use the synergy between artificial intelligence (AI) and high performance computing to maximize the potential for scientific discovery with Multi-Messenger Astrophysics. We discuss key aspects to realize this endeavor, namely (i) the design and exploitation of scalable and computationally efficient AI algorithms for Multi-Messenger Astrophysics; (ii) cyberinfrastructure requirements to numerically simulate astrophysical sources, and to process and interpret Multi-Messenger Astrophysics data; (iii) management of gravitational wave detections and triggers to enable electromagnetic and astro-particle follow-ups; (iv) a vision to harness future developments of machine and deep learning and cyberinfrastructure resources to cope with the scale of discovery in the Big Data Era; (v) and the need to build a community that brings domain experts together with data scientists on equal footing to maximize and accelerate discovery in the nascent field of Multi-Messenger Astrophysics.

astroquery is a collection of tools for requesting data from databases hosted on remote servers with interfaces exposed on the internet, including those with web pages but without formal application program interfaces (APIs). These tools are built on the Python requests package, which is used to make HTTP requests, and astropy, which provides most of the data parsing functionality. astroquery modules generally attempt to replicate the web page interface provided by a given service as closely as possible, making the transition from browser-based to command-line interaction easy. astroquery has received significant contributions from throughout the astronomical community, including several significant contributions from telescope archives. astroquery enables the creation of fully reproducible workflows from data acquisition through publication. This paper describes the philosophy, basic structure, and development model of the astroquery package. The complete documentation for astroquery can be found at http://astroquery.readthedocs.io/.

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 inter-operable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy project.

We present estimates of the basic physical properties -- including size and albedo -- of the extreme Centaur 2013 AZ60. These properties have been derived from optical and thermal infrared measurements. Our optical measurements revealed a likely full period of ~9.4 h with a shallow amplitude of 4.5%. By combining optical brightness information and thermal emission data, we are able to derive a diameter of 62.3 +/- 5.3 km and a geometric albedo of 2.9% -- corresponding to an extremely dark surface. Additionally, our finding of ~> 50 Jm^{-2}K^{-1}s^{-1/2} for the thermal inertia is also noticeably for objects in such a distance. The results of dynamical simulations yield an unstable orbit, with a 50% probability that the target will be ejected from the Solar System within 700,000 years. The current orbit of this object as well as its instability could imply a pristine cometary surface. This possibility is in agreement with the observed low geometric albedo and red photometric colour indices for the object, which are a good match for the surface of a dormant comet -- as would be expected for a long-period cometary body approaching perihelion. Despite the fact it was approaching ever closer to the Sun, however, the object exhibited star-like profiles in each of our observations, lacking any sign of cometary activity. By the albedo, 2013 AZ60 is a candidate for the darkest body among the known TNOs.

We report the discovery of WTS-2 b, an unusually close-in 1.02-day hot Jupiter (Mp=1.12MJ, Rp=1.363RJ) orbiting a K2V star, which has a possible gravitationally-bound M-dwarf companion at 0.6 arcsec separation contributing ~20 percent of the total flux in the observed J-band light curve. The planet is only 1.5 times the separation from its host star at which it would be destroyed by Roche lobe overflow, and has a predicted remaining lifetime of just ~40 Myr, assuming a tidal dissipation quality factor of Q'*=10^6. Q'* is a key factor in determining how frictional processes within a host star affect the orbital evolution of its companion giant planets, but it is currently poorly constrained by observations. We calculate that the orbital decay of WTS-2 b would correspond to a shift in its transit arrival time of T_shift~17 seconds after 15 years assuming Q'*=10^6. A shift less than this would place a direct observational constraint on the lower limit of Q'* in this system. We also report a correction to the previously published expected T_shift for WASP-18 b, finding that T_shift=356 seconds after 10 years for Q'*=10^6, which is much larger than the estimated 28 seconds quoted in WASP-18 b discovery paper. We attempted to constrain Q'* via a study of the entire population of known transiting hot Jupiters, but our results were inconclusive, requiring a more detailed treatment of transit survey sensitivities at long periods. We conclude that the most informative and straight-forward constraints on Q'* will be obtained by direct observational measurements of the shift in transit arrival times in individual hot Jupiter systems. We show that this is achievable across the mass spectrum of exoplanet host stars within a decade, and will directly probe the effects of stellar interior structure on tidal dissipation.

The WFCAM Transit Survey (WTS) is a near-infrared transit survey running on the United Kingdom Infrared Telescope (UKIRT), designed to discover planets around M dwarfs. The WTS acts as a poor-seeing backup programme for the telescope, and represents the first dedicated wide-field near-infrared transit survey. In this paper we describe the observing strategy of the WTS and the processing of the data to generate lightcurves. We describe the basic properties of our photometric data, and measure our sensitivity based on 950 observations. We show that the photometry reaches a precision of ~4mmag for the brightest unsaturated stars in lightcurves spanning almost 3 years. Optical (SDSS griz) and near-infrared (UKIRT ZYJHK) photometry is used to classify the target sample of 4600 M dwarfs with J magnitudes in the range 11-17. Most have spectral-types in the range M0-M2. We conduct Monte Carlo transit injection and detection simulations for short period (<10 day) Jupiter- and Neptune-sized planets to characterize the sensitivity of the survey. We investigate the recovery rate as a function of period and magnitude for 4 hypothetical star-planet cases: M0-2+Jupiter, M2-4+Jupiter, M0-2+Neptune, M2-4+Neptune. We find that the WTS lightcurves are very sensitive to the presence of Jupiter-sized short-period transiting planets around M dwarfs. Hot Neptunes produce a much weaker signal and suffer a correspondingly smaller recovery fraction. Neptunes can only be reliably recovered with the correct period around the rather small sample (~100) of the latest M dwarfs (M4-M9) in the WTS. The non-detection of a hot-Jupiter around an M dwarf by the WFCAM Transit Survey allows us to place an upper limit of 1.7-2.0 per cent (at 95 per cent confidence) on the planet occurrence rate.

We present an analysis of the photometric variability of M dwarfs in the WFCAM Transit Survey. Although periodic lightcurve variability in low mass stars is generally dominated by photospheric star spot activity, M dwarf variability in the J band has not been as thoroughly investigated as at visible wavelengths. Spectral type estimates for a sample of over 200,000 objects are made using spectral type-colour relations, and over 9600 dwarfs (J<17) with spectral types later than K7 were found. The light curves of the late-type sample are searched for periodicity using a Lomb-Scargle periodogram analysis. A total of 68 periodic variable M dwarfs are found in the sample with periods ranging from 0.16 days to 90.33 days, with amplitudes in the range of ~0.009 to ~0.115 in the J band. We simulate active M dwarfs with a range of latitude-independent spot coverages and estimate a periodically variable fractions of 1-3 per cent for stars where spots cover more than 10 per cent of the star's surface. Our simulated spot distributions indicate that operating in the J band, where spot contrast ratios are minimised, enables variability in only the most active of stars to be detected. These findings affirm the benefits of using the $J$ band for planetary transit searches compared to visible bands. We also serendipitously find a \Delta J>0.2 mag flaring event from an M4V star in our sample.

We report on the discovery of HAT-P-11b, the smallest radius transiting extrasolar planet (TEP) discovered from the ground, and the first hot Neptune discovered to date by transit searches. HAT-P-11b orbits the bright (V=9.587) and metal rich ([Fe=H] = +0.31 +/- 0.05) K4 dwarf star GSC 03561-02092 with P = 4.8878162 +/- 0.0000071 days and produces a transit signal with depth of 4.2 mmag. We present a global analysis of the available photometric and radial-velocity data that result in stellar and planetary parameters, with simultaneous treatment of systematic variations. The planet, like its near-twin GJ 436b, is somewhat larger than Neptune (17Mearth, 3.8Rearth) both in mass Mp = 0.081 +/- 0.009 MJ (25.8 +/- 2.9 Mearth) and radius Rp = 0.422 +/- 0.014 RJ (4.73 +/- 0.16 Rearth). HAT-P-11b orbits in an eccentric orbit with e = 0.198 +/- 0.046 and omega = 355.2 +/- 17.3, causing a reflex motion of its parent star with amplitude 11.6 +/- 1.2 m/s, a challenging detection due to the high level of chromospheric activity of the parent star. Our ephemeris for the transit events is Tc = 2454605.89132 +/- 0.00032 (BJD), with duration 0.0957 +/- 0.0012 d, and secondary eclipse epoch of 2454608.96 +/- 0.15 d (BJD). The basic stellar parameters of the host star are M* = 0.809+0.020-0.027 Msun, R* = 0.752 +/- 0.021 Rsun and Teff = 4780 +/- 50 K. Importantly, HAT-P-11 will lie on one of the detectors of the forthcoming Kepler mission. We discuss an interesting constraint on the eccentricity of the system by the transit light curve and stellar parameters. We also present a blend analysis, that for the first time treats the case of a blended transiting hot Jupiter mimicing a transiting hot Neptune, and proves that HAT-P-11b is not such a blend.

We present refined parameters for the extrasolar planetary system HAT-P-2 (also known as HD 147506), based on new radial velocity and photometric data. HAT-P-2b is a transiting extrasolar planet that exhibits an eccentric orbit. We present a detailed analysis of the planetary and stellar parameters, yielding consistent results for the mass and radius of the star, better constraints on the orbital eccentricity, and refined planetary parameters. The improved parameters for the host star are M_star = 1.36 +/- 0.04 M_sun and R_star = 1.64 +/- 0.08 R_sun, while the planet has a mass of M_p = 9.09 +/- 0.24 M_Jup and radius of R_p = 1.16 +/- 0.08 R_Jup. The refined transit epoch and period for the planet are E = 2,454,387.49375 +/- 0.00074 (BJD) and P = 5.6334729 +/- 0.0000061 (days), and the orbital eccentricity and argument of periastron are e = 0.5171 +/- 0.0033 and omega = 185.22 +/- 0.95 degrees. These orbital elements allow us to predict the timings of secondary eclipses with a reasonable accuracy of ~15 minutes. We also discuss the effects of this significant eccentricity including the characterization of the asymmetry in the transit light curve. Simple formulae are presented for the above, and these, in turn, can be used to constrain the orbital eccentricity using purely photometric data. These will be particularly useful for very high precision, space-borne observations of transiting planets.

We present HATNet observations of XO-5b, confirming its planetary nature based on evidence beyond that described in the announcement of Burke et al. (2008), namely, the lack of significant correlation between spectral bisector variations and orbital phase. In addition, using extensive spectroscopic measurements spanning multiple seasons, we investigate the relatively large scatter in the spectral line bisectors. We also examine possible blended stellar configurations (hierarchical triples, chance alignments) that can mimic the planet signals, and we are able to show that none are consistent with the sum of all the data. The analysis of the S activity index shows no significant stellar activity. Our results for the planet parameters are consistent with values in Burke et al. (2008), and we refine both the stellar and planetary parameters using our data. XO-5b orbits a slightly evolved, late G type star with mass M_s = 0.88 +/- 0.03, radius R_s = 1.08 +/- 0.04, and metallicity close to solar. The planetary mass and radius are M_p = 1.059 +/- 0.028 M_Jup and R_p = 1.109 +/- 0.050 R_Jup, respectively, corresponding to a mean density of 0.96 -0.11 +0.14 g/cm^3. The ephemeris for the orbit is P = 4.187757 +/- 0.000011, E= 2454552.67168 +/- 0.00029 (BJD) with transit duration of 0.1307 +/- 0.0013 d. By measuring four individual transit centers, we found no signs for transit timing variations. The planet XO-5b is notable for its anomalously high Safronov number, and has a high surface gravity when compared to other transiting exoplanets with similar period.

We report the identification of 32 transiting-planet candidates in HATNet field G205. We describe the procedures that we have used to follow up these candidates with spectroscopic and photometric observations, and we present a status report on our interpretation of the 28 candidates for which we have follow-up observations. Eight are eclipsing binaries with orbital solutions whose periods are consistent with their photometric ephemerides; two of these spectroscopic orbits are singled-lined and six are double-lined. For one of the candidates, a nearby but fainter eclipsing binary proved to be the source for the HATNet light curve, due to blending in the HATNet images. Four of the candidates were found to be rotating more rapidly than vsini = 50 km/s and were not pursued further. Thirteen of the candidates showed no significant velocity variation at the level of 0.5 to 1.0 km/s . Seven of these were eventually withdrawn as photometric false alarms based on an independent reanalysis using more sophisticated tools. Of the remaining six, one was put aside because a close visual companion proved to be a spectroscopic binary, and two were not followed up because the host stars were judged to be too large. Two of the remaining candidates are members of a visual binary, one of which was previously confirmed as the first HATNet transiting planet, HAT-P-1b. In this paper we confirm that the last of this set of candidates is also a a transiting planet, which we designate HAT-P-8b, with mass Mp = 1.52 +/- 0.18/0.16 Mjup, radius Rp = 1.50 +/- 0.08/0.06 Rjup, and photometric period P = 3.076320 +/- 0.000004 days. HAT-P-8b has an inflated radius for its mass, and a large mass for its period. The host star is a solar-metallicity F dwarf, with mass M* = 1.28 +/- 0.04 Msun and Rp = 1.58 +/- 0.08/0.06 Rsun.