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 a selection of methods for automatically constructing an optimal kernel model for difference image analysis which require very few external parameters to control the kernel design. Each method consists of two components; namely, a kernel design algorithm to generate a set of candidate kernel models, and a model selection criterion to select the simplest kernel model from the candidate models that provides a sufficiently good fit to the target image. We restricted our attention to the case of solving for a spatially-invariant convolution kernel composed of delta basis functions, and we considered 19 different kernel solution methods including six employing kernel regularisation. We tested these kernel solution methods by performing a comprehensive set of image simulations and investigating how their performance in terms of model error, fit quality, and photometric accuracy depends on the properties of the reference and target images. We find that the irregular kernel design algorithm employing unregularised delta basis functions, combined with either the Akaike or Takeuchi information criterion, is the best kernel solution method in terms of photometric accuracy. Our results are validated by tests performed on two independent sets of real data. Finally, we provide some important recommendations for software implementations of difference image analysis.

We present the Signal Detection using Random-Forest Algorithm (SIDRA). SIDRA is a detection and classification algorithm based on the Machine Learning technique (Random Forest). The goal of this paper is to show the power of SIDRA for quick and accurate signal detection and classification. We first diagnose the power of the method with simulated light curves and try it on a subset of the Kepler space mission catalogue. We use five classes of simulated light curves (CONSTANT, TRANSIT, VARIABLE, MLENS and EB for constant light curves, transiting exoplanet, variable, microlensing events and eclipsing binaries, respectively) to analyse the power of the method. The algorithm uses four features in order to classify the light curves. The training sample contains 5000 light curves (1000 from each class) and 50000 random light curves for testing. The total SIDRA success ratio is $\geq 90\%$. Furthermore, the success ratio reaches 95 - 100$\%$ for the CONSTANT, VARIABLE, EB, and MLENS classes and 92$\%$ for the TRANSIT class with a decision probability of 60$\%$. Because the TRANSIT class is the one which fails the most, we run a simultaneous fit using SIDRA and a Box Least Square (BLS) based algorithm for searching for transiting exoplanets. As a result, our algorithm detects 7.5$\%$ more planets than a classic BLS algorithm, with better results for lower signal-to-noise light curves. SIDRA succeeds to catch 98$\%$ of the planet candidates in the Kepler sample and fails for 7$\%$ of the false alarms subset. SIDRA promises to be useful for developing a detection algorithm and/or classifier for large photometric surveys such as TESS and PLATO exoplanet future space missions.

Context: Understanding the source of systematic errors in photometry is essential for their calibration. Aims: We investigate how photometry performed on difference images can be influenced by errors in the photometric scale factor. Methods: We explore the equations for difference image analysis (DIA) and we derive an expression describing how errors in the difference flux, the photometric scale factor and the reference flux are propagated to the object photometry. Results: We find that the error in the photometric scale factor is important, and while a few studies have shown that it can be at a significant level, it is currently neglected by the vast majority of photometric surveys employing DIA. Conclusions: Minimising the error in the photometric scale factor, or compensating for it in a post-calibration model, is crucial for reducing the systematic errors in DIA photometry.

We present an analysis of the Qatar-1 and TrES-5 transiting exoplanetary systems, which contain Jupiter-like planets on short-period orbits around K-dwarf stars. Our data comprise a total of 20 transit light curves obtained using five medium-class telescopes, operated using the defocussing technique. The average precision we reach in all our data is $RMS_{Q} = 1.1$ mmag for Qatar-1 ($V = 12.8$) and $RMS_{T} = 1.0$ mmag for TrES-5 ($V = 13.7$). We use these data to refine the orbital ephemeris, photometric parameters, and measured physical properties of the two systems. One transit event for each object was observed simultaneously in three passbands ($gri$) using the BUSCA imager. The QES survey light curve of Qatar-1 has a clear sinusoidal variation on a period of $P_{\star} = 23.697 \pm 0.123$\,d, implying significant starspot activity. We searched for starspot crossing events in our light curves, but did not find clear evidence in any of the new datasets. The planet in the Qatar-1 system did not transit the active latitudes on the surfaces of its host star. Under the assumption that $P_{\star}$ corresponds to the rotation period of Qatar-1\,A, the rotational velocity of this star is very close to the $v \sin i_\star$ value found from observations of the Rossiter-McLaughlin effect. The low projected orbital obliquity found in this system thus implies a low absolute orbital obliquity, which is also a necessary condition for the transit chord of the planet to avoid active latitudes on the stellar surface.

We present three newly discovered sub-Jupiter mass planets from the SuperWASP survey: WASP-54b is a heavily bloated planet of mass 0.636$^{+0.025}_{-0.024}$ \mj and radius 1.653$^{+0.090}_{-0.083}$ \rj. It orbits a F9 star, evolving off the main sequence, every 3.69 days. Our MCMC fit of the system yields a slightly eccentric orbit ($e=0.067^{+0.033}_{-0.025}$) for WASP-54b. We investigated further the veracity of our detection of the eccentric orbit for WASP-54b, and we find that it could be real. However, given the brightness of WASP-54 V=10.42 magnitudes, we encourage observations of a secondary eclipse to draw robust conclusions on both the orbital eccentricity and the thermal structure of the planet. WASP-56b and WASP-57b have masses of 0.571$^{+0.034}_{-0.035}$ \mj and $0.672^{+0.049}_{-0.046}$ \mj, respectively; and radii of $1.092^{+0.035}_{-0.033}$ \rj for WASP-56b and $0.916^{+0.017}_{-0.014}$ \rj for WASP-57b. They orbit main sequence stars of spectral type G6 every 4.67 and 2.84 days, respectively. WASP-56b and WASP-57b show no radius anomaly and a high density possibly implying a large core of heavy elements; possibly as high as $\sim$50 M$_{\oplus}$ in the case of WASP-57b. However, the composition of the deep interior of exoplanets remain still undetermined. Thus, more exoplanet discoveries such as the ones presented in this paper, are needed to understand and constrain giant planets' physical properties.

We have studied over 1600 Am stars at a photometric precision of 1 mmag with SuperWASP photometric data. Contrary to previous belief, we find that around 200 Am stars are pulsating delta Sct and gamma Dor stars, with low amplitudes that have been missed in previous, less extensive studies. While the amplitudes are generally low, the presence of pulsation in Am stars places a strong constraint on atmospheric convection, and may require the pulsation to be laminar. While some pulsating Am stars have been previously found to be delta Sct stars, the vast majority of Am stars known to pulsate are presented in this paper. They will form the basis of future statistical studies of pulsation in the presence of atomic diffusion.

We present light curves and periods of 53 candidates for short period eclipsing binary stars identified by SuperWASP. These include 48 newly identified objects with periods <2x10^4 seconds (~0.23d), as well as the shortest period binary known with main sequence components (GSC2314-0530 = 1SWASP J022050.85+332047.6) and four other previously known W UMa stars (although the previously reported periods for two of these four are shown to be incorrect). The period distribution of main sequence contact binaries shows a sharp cut-off at a lower limit of around 0.22d, but until now, very few systems were known close to this limit. These new candidates will therefore be important for understanding the evolution of low mass stars and to allow investigation of the cause of the period cut-off.

The WASP (Wide Angle Search for Planets) project is an exoplanet transit survey that has been automatically taking wide field images since 2004. Two instruments, one in La Palma and the other in South Africa, continually monitor the night sky, building up light curves of millions of unique objects. These light curves are used to search for the characteristics of exoplanetary transits. This first public data release (DR1) of the WASP archive makes available all the light curve data and images from 2004 up to 2008 in both the Northern and Southern hemispheres. A web interface (www.wasp.le.ac.uk/public/) to the data allows easy access over the Internet. The data set contains 3 631 972 raw images and 17 970 937 light curves. In total the light curves have 119 930 299 362 data points available between them.

We report the discovery from WASP-South of a new Jupiter-like extrasolar planet, WASP-16b, which transits its solar analog host star every 3.12 days. Analysis of the transit photometry and radial velocity spectroscopic data leads to a planet with R_p = 1.008+-0.071 R_Jup and M_p =0.855+-0.059 M_Jup, orbiting a host star with R_* = 0.946+-0.054 R_sun and M_* = 1.022+-0.101 M_sun. Comparison of the high resolution stellar spectrum with synthetic spectra and stellar evolution models indicates the host star is a near-solar metallicity ([Fe/H]= 0.01+-0.10) solar analog (Teff = 5700+-150 K, log g= 4.5+-0.2) of intermediate age (Tau = 2.3+5.8-2.2 Gyr).

We present the results of a photometric survey of rotation rates in the Coma Berenices (Melotte 111) open cluster, using data obtained as part of the SuperWASP exoplanetary transit-search programme. The goal of the Coma survey was to measure precise rotation periods for main-sequence F, G and K dwarfs in this intermediate-age (~600 Myr) cluster, and to determine the extent to which magnetic braking has caused the stellar spin periods to converge. We find a tight, almost linear relationship between rotation period and J-K colour with a root-mean square scatter of only 2 percent. The relation is similar to that seen among F, G and K stars in the Hyades. Such strong convergence can only be explained if angular momentum is not at present being transferred from a reservoir in the deep stellar interiors to the surface layers. We conclude that the coupling timescale for angular momentum transport from a rapidly-spinning radiative core to the outer convective zone must be substantially shorter than the cluster age, and that from the age of Coma onward, stars rotate effectively as solid bodies. The existence of a tight relationship between stellar mass and rotation period at a given age supports the use of stellar rotation period as an age indicator in F, G and K stars of Hyades age and older. We demonstrate that individual stellar ages can be determined within the Coma population with an internal precision of order 9 percent (RMS), using a standard magnetic braking law in which rotation period increases with the square root of stellar age. We find that a slight modification to the magnetic-braking power law, P proportional to t^0.56, yields rotational and asteroseismological ages in good agreement for the Sun and other stars of solar age for which p-mode studies and photometric rotation periods have been published.

We report the discovery of WASP-13b, a low-mass $ M_p = 0.46 ^{+ 0.06}_{- 0.05} M_J$ transiting exoplanet with an orbital period of $4.35298 \pm 0.00004$ days. The transit has a depth of 9 mmag, and although our follow-up photometry does not allow us to constrain the impact parameter well ($0 < b < 0.46$), with radius in the range $R_p \sim 1.06 - 1.21 R_J$ the location of WASP-13b in the mass-radius plane is nevertheless consistent with H/He-dominated, irradiated, low core mass and core-free theoretical models. The G1V host star is similar to the Sun in mass (M$_{*} = 1.03^{+0.11}_ {- 0.09} M_{\odot}$) and metallicity ([M/H]=$0.0\pm0.2$), but is possibly older ($8.5^{+ 5.5}_{- 4.9}$ Gyr).

We report the discovery of a low-density exoplanet transiting an 11th magnitude star in the Southern hemisphere. WASP-15b, which orbits its host star with a period P=3.7520656+-0.0000028d has a mass M_p=0.542+-0.050M_J and radius R_p=1.428+-0.077R_J, and is therefore the one of least dense transiting exoplanets so far discovered (rho_p=0.247+-0.035g cm^-3). An analysis of the spectrum of the host star shows it to be of spectral type around F5, with an effective temperature T_eff=6300+-100K and [Fe/H]=-0.17+-0.11.

We report the discovery of WASP-6b, an inflated sub-Jupiter mass planet transiting every 3.3610060 +0.0000022-0.0000035 days a mildly metal-poor solar-type star of magnitude V=11.9. A combined analysis of the WASP photometry, high-precision followup transit photometry and radial velocities yield a planetary mass M_p = 0.503 +0.019-0.038 M_jup and radius R_p = 1.224 +0.051-0.052 R_jup, resulting in a density rho_p = 0.27 +-0.05 rho_jup. The mass and radius for the host star are M_s = 0.88 +0.05-0.08 M_sun and R_s = 0.870 +0.025-0.036 R_sun. The non-zero orbital eccentricity e = 0.054 +0.018-0.015 that we measure suggests that the planet underwent a massive tidal heating ~1 Gyr ago that could have contributed to its inflated radius. High-precision radial velocities obtained during a transit allow us to measure a sky-projected angle between the stellar spin and orbital axis Beta = 11 +14-18 deg. In addition to similar published measurements, this result favors a dominant migration mechanism based on tidal interactions with a protoplanetary disk.

We report the discovery of WASP-10b, a new transiting extrasolar planet (ESP) discovered by the WASP Consortium and confirmed using NOT FIES and SOPHIE radial velocity data. A 3.09 day period, 29 mmag transit depth, and 2.36 hour duration are derived for WASP-10b using WASP and high precision photometric observations. Simultaneous fitting to the photometric and radial velocity data using a Markov-chain Monte Carlo procedure leads to a planet radius of 1.28R_J, a mass of 2.96M_J and eccentricity of ~0.06. WASP-10b is one of the more massive transiting ESPs, and we compare its characteristics to the current sample of transiting ESP, where there is currently little information for masses greater than ~2M_J and non-zero eccentricities. WASP-10's host star, GSC 2752-00114 (USNO-B1.0 1214-0586164) is among the fainter stars in the WASP sample, with V=12.7 and a spectral type of K5. This result shows promise for future late-type dwarf star surveys.

We report the discovery of a 7.3 Mjup exoplanet WASP-14b, one of the most massive transiting exoplanets observed to date. The planet orbits the tenth-magnitude F5V star USNO-B1 11118-0262485 with a period of 2.243752 days and orbital eccentricity e = 0.09. A simultaneous fit of the transit light curve and radial velocity measurements yields a planetary mass of 7.3+/-0.5 Mjup and a radius of 1.28+/-0.08 Rjup. This leads to a mean density of about 4.6 g/cm^3 making it densest transiting exoplanets yet found at an orbital period less than 3 days. We estimate this system to be at a distance of 160+/-20 pc. Spectral analysis of the host star reveals a temperature of 6475+/-100 K, log g = 4.07 cm/s^2 and vsin i = 4.9+/-1.0 km/s, and also a high lithium abundance, log N(Li} = 2.84+/-0.05. The stellar density, effective temperature and rotation rate suggest an age for the system of about 0.5-1.0 Gyr.

We report the discovery of WASP-4b, a large transiting gas-giant planet with an orbital period of 1.34 days. This is the first planet to be discovered by the SuperWASP-South observatory and CORALIE collaboration and the first planet orbiting a star brighter than 16th magnitude to be discovered in the Southern hemisphere. A simultaneous fit to high-quality lightcurves and precision radial-velocity measurements leads to a planetary mass of 1.22 +/- 0.1 MJup and a planetary radius of 1.42 +/- 0.08 RJup. The host star is USNO-B1.0 0479-0948995, a G7V star of visual magnitude 12.5. As a result of the short orbital period, the predicted surface temperature of the planet is 1776 K, making it an ideal candidate for detections of the secondary eclipse at infrared wavelengths.

We report the discovery of WASP-3b, the third transiting exoplanet to be discovered by the WASP and SOPHIE collaboration. WASP-3b transits its host star USNO-B1.0 1256-0285133 every 1.846834+-0.000002 days. Our high precision radial-velocity measurements present a variation with amplitude characteristic of a planetary-mass companion and in-phase with the light-curve. Adaptive optics imaging shows no evidence for nearby stellar companions, and line-bisector analysis excludes faint, unresolved binarity and stellar activity as the cause of the radial-velocity variations. We make a preliminary spectroscopic analysis of the host star finding it to have Teff = 6400+-100 K and log g = 4.25+-0.05 which suggests it is most likely an unevolved main sequence star of spectral type F7-8V. Our simultaneous modelling of the transit photometry and reflex motion of the host leads us to derive a mass of 1.76 +0.08 -0.14 M_J and radius 1.31 +0.07-0.14 R_J for WASP-3b. The proximity and relative temperature of the host star suggests that WASP-3b is one of the hottest exoplanets known, and thus has the potential to place stringent constraints on exoplanet atmospheric models.

Transiting extrasolar planets constitute only a small fraction of the range of stellar systems found to display periodic, shallow dimmings in wide-field surveys employing small-aperture camera arrays. Here we present an efficient selection strategy for follow-up observations, derived from analysis of the light curves of a sample of 67 SuperWASP targets that passed the selection tests we used in earlier papers, but which have subsequently been identified either as planet hosts or as astrophysical false positives. We determine the system parameters using Markov-chain Monte Carlo analysis of the SuperWASP light curves. We use a constrained optimisation of chi-squared combined with a Bayesian prior based on the main-sequence mass and radius expected from the 2MASS J-H colour. The Bayesian nature of the analysis allows us to quantify both the departure of the host star from the main-sequence mass-radius relation and the probability that the companion radius is less than 1.5 Jupiter radii. When augmented by direct light curve analyses that detect binaries with unequal primary and secondary eclipses, and objects with aperture blends that are resolved by SuperWASP, we find that only 13 of the original 67 stars, including the three known planets in the sample, would qualify for follow-up. This suggests that planet discovery "hit rates" better than one-in-five should be achievable. In addition, the stellar binaries that qualify are likely to have astrophysically interesting stellar or sub-stellar secondaries.

We present optical lightcurves of 428 periodic variable stars coincident with ROSAT X-ray sources, detected using the first run of the SuperWASP photometric survey. Only 68 of these were previously recognised as periodic variables. A further 30 of these objects are previously known pre-main sequence stars, for which we detect a modulation period for the first time. Amongst the newly identified periodic variables, many appear to be close eclipsing binaries, their X-ray emission is presumably the result of RS CVn type behaviour. Others are probably BY Dra stars, pre-main sequence stars and other rapid rotators displaying enhanced coronal activity. A number of previously catalogued pulsating variables (RR Lyr stars and Cepheids) coincident with X-ray sources are also seen, but we show that these are likely to be misclassifications. We identify four objects which are probable low mass eclipsing binary stars, based on their very red colour and light curve morphology.

A recent study demonstrated that there is significant covariance structure in the noise on data from ground-based photometric surveys designed to detect transiting extrasolar planets. Such correlation in the noise has often been overlooked, especially when predicting the number of planets a particular survey is likely to find. Indeed, the shortfall in the number of transiting extrasolar planets discovered by such surveys seems to be explained by co-variance in the noise. We analyse SuperWASP (Wide Angle Search for Planets) data and determine that there is a significant amount of correlated systematic noise present. After modelling the potential planet catch, we conclude that this noise places a significant limit on the number of planets that SuperWASP is likely to detect; and that the best way to boost the signal-to-noise ratio and limit the impact of co-variant noise is to increase the number of observed transits for each candidate transiting planet.

We present a model of the stellar populations in the fields observed by one of the SuperWASP-N cameras in the 2004 observing season. We use the Besancon Galactic model to define the range of stellar types and metallicities present, and populate these objects with transiting extra-solar planets using the metallicity relation of Fischer & Valenti (2005). We investigate the ability of SuperWASP to detect these planets in the presence of realistic levels of correlated systematic noise (`red noise'). We find that the number of planets that transit with a signal-to-noise ratio of 10 or more increases linearly with the number of nights of observations. Based on a simulation of detection rates across 20 fields observed by one camera, we predict that a total of 18.6 \pm 8.0 planets should be detectable from the SuperWASP-N 2004 data alone. The best way to limit the impact of co-variant noise and increase the number of detectable planets is to boost the signal-to-noise ratio, by increasing the number of observed transits for each candidate transiting planet. This requires the observing baseline to be increased, by spending a second observing season monitoring the same fields.

We have detected low-amplitude radial-velocity variations in two stars, USNO-B1.0 1219-0005465 (GSC 02265-00107 = WASP-1) and USNO-B1.0 0964-0543604 (GSC 00522-01199 = WASP-2). Both stars were identified as being likely host stars of transiting exoplanets in the 2004 SuperWASP wide-field transit survey. Using the newly-commissioned radial-velocity spectrograph SOPHIE at the Observatoire de Haute-Provence, we found that both objects exhibit reflex orbital radial-velocity variations with amplitudes characteristic of planetary-mass companions and in-phase with the photometric orbits. Line-bisector studies rule out faint blended binaries as the cause of either the radial-velocity variations or the transits. We perform preliminary spectral analyses of the host stars, which together with their radial-velocity variations and fits to the transit light curves, yield estimates of the planetary masses and radii. WASP-1b and WASP-2b have orbital periods of 2.52 and 2.15 days respectively. Given mass estimates for their F7V and K1V primaries we derive planet masses 0.80 to 0.98 and 0.81 to 0.95 times that of Jupiter respectively. WASP-1b appears to have an inflated radius of at least 1.33 R_Jup, whereas WASP-2b has a radius in the range 0.65 to 1.26 R_Jup.