{'fraction_non_alphanumeric': 0.03320005711837784, 'fraction_numerical': 0.03374744157265934, 'mean_word_length': 4.572006365203554, 'pattern_counts': {'":': 0, '<': 9, '<?xml version=': 0, '>': 3, 'https://': 0, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
High-resolution Imaging of Transiting Extrasolar Planetary systems (HITEP). II. Lucky Imaging results from 2015 and 2016

The formation and dynamical history of hot Jupiters is currently debated, with wide stellar binaries having been suggested as a potential formation pathway. Additionally, contaminating light from both binary companions and unassociated stars can significantly bias the results of planet characterisation studies, but can be corrected for if the properties of the contaminating star are known. We search for binary companions to known transiting exoplanet host stars, in order to determine the multiplicity properties of hot Jupiter host stars. We also characterise unassociated stars along the line of sight, allowing photometric and spectroscopic observations of the planetary system to be corrected for contaminating light. We analyse lucky imaging observations of 97 Southern hemisphere exoplanet host stars, using the Two Colour Instrument on the Danish 1.54m telescope. For each detected companion star, we determine flux ratios relative to the planet host star in two passbands, and measure the relative position of the companion. The probability of each companion being physically associated was determined using our two-colour photometry. A catalogue of close companion stars is presented, including flux ratios, position measurements, and estimated companion star temperature. For companions that are potential binary companions, we review archival and catalogue data for further evidence. For WASP-77AB and WASP-85AB, we combine our data with historical measurements to determine the binary orbits, showing them to be moderately eccentric and inclined to the line of sight and planetary orbital axis. Combining our survey with the similar Friends of Hot Jupiters survey, we conclude that known hot Jupiter host stars show a deficit of high mass stellar companions compared to the field star population; however, this may be a result of the biases in detection and target selection by ground-based surveys.


Introduction
In the last two decades, many exoplanets have been discovered in orbital configurations that provide challenges for planet formation theory. One class of planets which has generated much interest are the hot Jupiters, a group of gas giants living extremely close to their host stars, with orbital periods of less than 10 days, with some examples having periods below a day. At the time of their discovery, planet formation theory suggested that gas giants form far from their host stars, where temperatures in the protoplanetary disc were low enough for volatiles -such as water and methane -to condense, permitting these planets to accrete sufficient mass to hold onto a gaseous envelope (e.g. Pollack et al. 1996;Boss 1995). Theoretical models aiming to explain the formation of these planets generally suggest that the volves the exchange of angular momentum between the binary and planetary orbits (Wu & Murray 2003;Fabrycky & Tremaine 2007;Naoz et al. 2012).
If wide binaries are indeed responsible for hot Jupiter formation, it would therefore be expected that many hot Jupiters would be found in binary systems. It is also of interest to determine the stellar environments in which hot Jupiters do not exist. For example, it is expected that binary companions that are too close will inhibit planet formation (Jang-Condell 2015), with some observational evidence supporting this theory (Kraus et al. 2016). It has also been suggested that wide binary companions, 1000au or so from the planet host star, could also have an influence on the planet at later stages if the binary orbit is modified by the galactic tide and close stellar encounters (Kaib et al. 2013). High resolution imaging allows the detection and confirmation of wide stellar binaries, through the use of multi-colour photometry and common proper motion analysis.
High resolution imaging is also able to discover physically unassociated stars in the vicinity of suspected transiting exoplanet host stars. Transit surveys, especially those operating from the ground, suffer from a high rate of false positives, with grazing, blended, or high mass ratio eclipsing binaries being able to mimic the signals of a planet (Brown 2003;Günther et al. 2017). However, space-based missions are also liable to false positives, due to faint targets or small photometric/spectroscopic signals making ground-based follow-up observations difficult, as was recently shown in the case of three claimed transiting planets from the K2 mission (Cabrera et al. 2017).
When the planetary signal is real, it is still important to consider the effect of other stars. Contaminating light can significantly alter the derived bulk properties of a planet (Buchhave et al. 2011;Evans et al. 2016b). Additionally, a nearby star dissimilar to the planet host star can alter or mimic various features that are expected to be seen in planetary atmospheres, with cool, red stars producing slopes towards the blue that resemble Rayleigh scattering (Southworth & Evans 2016).
In this paper, we present the continuation of our lucky imaging survey of transiting exoplanet host stars in the Southern hemisphere. In Evans et al. (2016a), hereafter Paper I, we presented lucky imaging observations of 101 planet host stars, detecting 51 companions located within 5 arcseconds of a target star. The current paper covers lucky imaging observations of 97 targets, including re-observations of several interesting objects, with two colour photometry allowing us to assess the likelihood of each detected companion star being bound.

Observations
Lucky imaging observations of 97 targets were carried out using the Two Colour Instrument (TCI) on the Danish 1.54m telescope, located at La Silla, Chile. Following the observing strategy in Paper I, our targets are all publicly announced transiting exoplanet host stars, taken from the TEPCat 1 database of wellstudied transiting extrasolar planets (Southworth 2011). The observations reported in this work were taken during the periods April-September 2015 and April-September 2016. Targets were generally observed once in each year, with new systems being added to our target lists as they were published. Targets for which no companions were detected in our 2014 and 2015 observations, or for which companions were determined to be physically unassociated, were not re-observed in 2016. A summary of targets observed as part of this work is given in Table 1, with further 1 http://www.astro.keele.ac.uk/jkt/tepcat/ Table 1. Summary of transiting exoplanet systems considered in this paper. Systems with potential or confirmed binary companions are shown in bold (see Section 4). The TCI is a lucky imager designed for simultaneous twocolour photometry, using Electron Multiplying CCD (EMCCD) detectors, and is described in detail by Skottfelt et al. (2015b). Light is split between two cameras by a dichroic, with the 'red' camera receiving light redwards of 655nm, and the 'visual' camera receiving light between 466nm and 655nm, with a second dichroic sending light bluewards of 466nm to a focus system. No filters are fitted to the system, and so the passband for each camera is defined only by the dichroics and the sensitivity of the EMCCD chips. We refer to the passband of the visual and red cameras as v TCI and r TCI respectively. The designed passband of the red camera is similar to the combination of the SDSS i and z filters, or a wider version of the Cousins I filter. The visual camera has a central wavelength that is located between the SDSS g and r filters, similar to the Johnson V filter's central wavelength, but has a significantly different response to any of the mentioned filters. Both detectors consist of a 512×512 pixel array with a pixel scale of approximately ∼0.09 arcsec/pixel, giving a field of view of approximately 45 × 45 .

Systems studied in this work
The target systems were selected by brightness (9 ≤ V ≤ 15), with exposure times chosen to match this brightness, ranging from 60s to 900s. For most targets, the default electron multiplication gain of 300 e − /photon was used, but targets brighter than V = 10.5 were observed with a lower gain of 100 e − /photon, in order to keep counts within the range of the analogue to digital converter. Observations were performed at an airmass of 1.5 or less, and in seeing conditions of 0.6 or better. Targets were observed simultaneously with the visual and red cameras, with the exception of the observation of WASP-104 on the night of 2015-

Data reduction and analysis
In general, our reduction process follows the outline given in Paper I . No major changes have been made to the reduction pipeline, as described in Harpsøe et al. (2012) and Skottfelt et al. (2015b), which automatically reduces data from the TCI. The process includes bias and flat field corrections, subtraction of cosmic rays, tip-tilt correction, and measurement of the quality of the lucky imaging exposures. The ranked exposures are then divided into ten cuts, with the exposures in each cut being combined into a single stacked image. The percentage boundaries between each cut are: 1,2,5,10,20,50,90,98,99,100. Therefore, the first cut contains the top 1% of images ranked by quality, the second cut contains the next 1% of images, and so on.
We apply the correction for smearing described in Paper I, which removes apparent smears across the image that appear when bright stars are present, thought to be caused by charge transfer inefficiency. We then run our star detection algorithm, again described in Paper I, which remains unchanged. In brief, we take the first seven cuts in quality (covering the best 90% of the exposures), and convolve each cut with a Gaussian of standard deviation 4 pixels (FWHM 11.7 pixels). Each convolved image is then subtracted from the original, in order to filter out high frequency noise in the image. We then divide each image into a series of annuli centred on the target star, with each annulus having a width of 0.5 pixels. For each annulus, a sigmaclipped mean and standard deviation are calculated, and any pixels more than 1.7σ from the annulus mean are flagged as a candidate detection, this value having been chosen as the best compromise between false positives and missed detections. The detections from each of the seven cuts are then combined, with any groups of adjacent pixels marked as detections being grouped together into a single result.
Each detection is then verified by eye using an interactive program, with the main source of spurious detections being asymmetries of the target star's PSF. Additionally, in crowded fields it is possible for several overlapping stars to be combined by the pipeline into a single detection, or for the identified position of a faint star to be offset slightly due to an overlapping bright star; the user is able to correct the identifications in such cases.
The star detection algorithm was run only on the data from the red camera, as this passband is more sensitive to low mass stellar companions. Additionally, the visual images generally have a lower angular resolution; firstly, is not possible to independently focus the red and visual cameras, with the instrument focus being set based on the red camera images. As a result; secondly, due to the turbulent properties of the atmosphere, the effects of seeing are worse at shorter wavelengths (Fried 1966).

Detection limits
For a given pixel on an image, we determined the mean and standard deviation of the counts inside a box of n × n pixels, centred on the pixel in question, with n being the FWHM of stars on the image. The number of counts corresponding to a 5σ excess were then calculated, giving the faintest detectable star. This process was repeated for each pixel on the image, with pixels then being binned by separation from the target star, and the median detection limit for each bin calculated. Limits were converted from A region that is entirely black indicates that a companion would not be detected in any observation, whilst a region that is white indicates that a companion would be detected in all observations. All detection curves remain essentially flat beyond 10 . counts to magnitudes using the counts on the central star within the same n×n box. This was repeated for each of the seven lucky imaging cuts, from which a final detection limit was compiled for each observation. Figs. 1 and 2 illustrate the detection limits for all targets in 2015 and 2016 respectively. Additionally, the detection limits for each target are tabulated in Tables 2 and 3.

Photometric calibration
We adopt the set of theoretical colour indices presented in Paper I, which give the expected magnitude difference between v TCI and r TCI , for main sequence stars and brown dwarfs between 2300K and 12000K. For stars below 7000K, the expected luminosity in both the visual and red was also presented, our calibration for radius (and hence luminosity) not being valid above Table 2. Observations from 2015 analysed in this work. The "1% FWHM" column lists the FWHM of the target star in the best 1% of lucky imaging exposures. The "Detection Limit" columns indicate the calculated detection limits in r TCI , determined using the method discussed in Section 3.  Table 3. Observations from 2016 analysed in this work. The "1% FWHM" column lists the FWHM of the target star in the best 1% of lucky imaging exposures. The "Detection Limit" columns indicate the calculated detection limits in r TCI , determined using the method discussed in Section 3.1. The full version of table is available in electronic form at the CDS. this temperature. These two calibrations allowed us to estimate the temperature of a candidate companion star, and then to determine whether its brightness is consistent with being at the same distance as the target star, or instead that it is too bright or too faint compared to the target, indicating that it is closer or further away.
With the much more extensive two-colour dataset presented in this work, it became apparent that our previous temperatureradius calibration was not sufficiently accurate below ∼ 3500K, with known physical companions being fainter than expected. Comparison with the Baraffe et al. (2015) stellar models for an age of 1Gyr shows that our radius calibration agrees to within a few percent down to 3500K, but predicts much larger radii than the models at cooler temperatures. We therefore use the Baraffe et al. (2015) model radii below 4000K, where the two sets of radii begin to deviate.
We adopt the previously calculated atmospheric extinction of 0.08 ± 0.03 mag/airmass for the (v − r) TCI colour index. In Paper I, we calibrated the predicted values of (v − r) TCI to the measured values. We present a new zero point offset based on our 2015 and 2016 data of −0.405±0.011 magnitudes -i.e. a star with a predicted colour of 0.000 would be measured as 0.405.
Our analysis in Paper I showed a potential temperaturedependent colour offset, with hot standard stars having a slightly larger offset of ∼ 0.50 magnitudes, whilst cooler stars clustered around an offset of ∼ 0.44 magnitudes. The sample of target stars in this paper cover 3270 − 6600K, and we find no correlation between temperature and colour offset. We find a scatter of 0.07mag in our colour measurements, similar to the variation from 0.44 to 0.50 magnitudes. We attribute the scatter to atmospheric transparency variations, errors in the temperatures from which the expected colours were derived, and intrinsic physical variations due to effects such as metallicity.

Host star distances
To convert a companion's angular separation from the host stars to projected separation, and hence begin to consider the orbital properties of the companion, we must know the distance to the host star. Many, but not all, hot Jupiter host stars have photometric or spectroscopic distance estimates published; however, the methodology used to derive these distances is varied, and in a number of cases, is simply not stated. With the release of the Tycho-Gaia astrometric solutions (TGAS, Michalik et al. 2015), approximately half of the TEP hosts in surveys such as WASP and HAT now have measured parallaxes, but the uncertainties on the derived distances become significant beyond 300pc.
We therefore choose to derive our own distance estimates for our targets in a homogeneous fashion, using the K-band surface brightness-effective temperature relation presented in Kervella et al. (2004), as previously done in Paper I. Stellar radii and effective temperatures were taken from the TEPCat database, using the version of the database from 2016-05-01 2 . K band magnitudes were taken from 2MASS K s magnitudes (Skrutskie et al. 2006), using the relation K = K s + 0.044 (Carpenter 2001).
To correct for interstellar extinction, we derived an initial distance estimate with no correction for extinction, from which we then obtained a value for E (B − V) from the three-dimensional dust map of Green et al. (2015), which provides measurements at several distances along a given line of sight. The extinction in K s was then calculated, using A(K s ) = 0.306 E (B − V) (Yuan et al. 2013), and the extinction-corrected value of K s used to calculate a new distance estimate. This process was iterated until the change in calculated distance was less than 1% between iterations.
The E (B − V) data provided by Green et al. (2015) is tabulated in steps of 0.5 in distance modulus, between which we interpolated linearly to get our values of E (B − V). For stars closer than the nearest reddening measurement provided by the map, we assumed E (B − V) = 0 at zero distance, and again interpolated linearly out to the nearest measurement. The map does not cover the entire sky, with a region towards the Southern pole currently being unmapped -being away from the galactic plane, this is an area of relatively low reddening, and we therefore assumed extinction was negligible for targets in this area.  Fig. 3. A comparison between previously published distances d Lit and distances derived from the TGAS parallaxes d Gaia . The figure shows all published HAT and WASP systems where both values are available, with each system plotted using the log of the mean of the two distance values, and the fractional difference between the two values. The green data point at top left illustrates the median error bar, with individual error bars excluded for clarity. The solid and dashed lines indicate the mean difference and RMS scatter between the two sets of data; the dotted line shows where points would fall if there were no systematic offset. We find that the exoplanet literature distances are systematically 4.4 ± 1.5% lower than the TGAS distances, with an RMS scatter of 14% between the two sets of data.
We compare both our distance values and those in the exoplanet literature for all HAT and WASP survey planets to those derived from the TGAS dataset in Figures 3 and 4. We adopt the distances and uncertainties calculated for the TGAS dataset by Astraatmadja & Bailer-Jones (2016) using their anisotropic 'Milky Way' prior on the distribution of stars, which shows the best agreement with Cepheid distances within 2kpc. We find that the literature distances are systematically lower than those found by TGAS, with a mean ratio of d Lit /d TGAS = 0.956±0.015 (2.9σ) with an RMS scatter of 14%. Our K-band derived distances show a negligible bias and a smaller scatter, with d Ks /d TGAS = 0.997 ± 0.018 (0.2σ) and an RMS scatter of 9%. Our derived distances, and the values used to calculate them, can be found in Table 4.

Interstellar reddening
As part of our analysis of companions, temperatures are derived from our photometric data, using the colour-temperature calibrations in Paper I. To correct for interstellar reddening, which would bias our results towards lower temperatures, we use the E (B − V) values determined in the previous section, and listed in Table 4. Between approximately 350nm and 1000nm, the relationship between interstellar extinction (measured in magnitudes) and wavelength is almost linear (e.g. Fitzpatrick 1999). We therefore assumed that the TCI colour excess E (v − r) TCI is linearly related to the standard Johnson E (B − V) colour excess by some constant X, We used the effective temperatures from TEPCat to calculate the expected intrinsic colours for our observed target stars. The colours for each target were then measured from our data, cor-  Fig. 3, comparing distances derived using the K-band surface brightness-effective temperature relation d Ks to those derived from the TGAS parallaxes d Gaia . We find that our derived distances are in good agreement with the TGAS-derived distances, with our values being only 0.3 ± 1.7% higher on average, with an RMS scatter of 9% between the two sets of data. Obs. minus calc. colour Fig. 5. The offset between the expected and measured colours of our target stars after correction for atmospheric effects and interstellar reddening, plotted against the expected E (B − V) for each target. Note that the x axis is logarithmic, as the majority of our targets have very low reddening values. All targets agree to within ±0.15 dex, corresponding to agreement within ±300K.
rected for the zero-point offset and atmospheric extinction, and the colour excess (difference between observed and calculated) was derived. A least squares fit was then performed to derive a value for the reddening factor X, using the measured colour excesses and the value of E (B − V) that was derived when calculating the distance to each star. We calculated a value of X = 0.53 ± 0.06, and the agreement between expected and measured colours after applying this correction is shown in Fig. 5.

Companion colour analysis
From our two-colour photometry and the colour indices presented in Paper I, we are able to estimate the temperature of detected companion stars. By comparing this to the literature Table 4. Assumed properties of our target stars, including distances estimated using the K-band surface brightness-effective temperature relation, and the values used to calculate the distances. Where no entry is given in the E (B − V) column, no extinction estimate was available, and it was assumed that the effect of extinction was negligible. The full version of temperature of the host star and assuming both are main sequence stars, we then calculate the expected flux ratio in both v TCI and r TCI if the target and companion form a physical binary and hence are at the same distance, using our temperature-radius calibration. Background stars will then appear fainter than expected, whilst foreground stars will be too bright, allowing us to assign a low probability of physical association to stars in a single epoch, without long-term proper motion monitoring. However, further evidence is required to reach definite conclusions. Chance alignments can result in stars with similar distances to appear close together, especially in cases where the target star is more distant or located in the plane of the galaxy. In addition, companions may not be main sequence stars -for example, the K2-02 (HIP 116454) system includes a white dwarf companion , which would be classified as a distant background object from its blue colour.
No further observations of standard stars have occurred since Paper I, and the available data do not allow us to meaningfully constrain atmospheric extinction, which could bias our colourbased temperatures. To account for this, we assume that the colours of the target stars should match the colour predicted by the published spectroscopic temperatures listed in Table 4, and use the offset between the measured and expected colours to determine the atmospheric extinction for each observation.

Astrometric calibration
To determine the detector scales and orientations for both the visual and red cameras, we compared TCI observations of the cores of globular clusters (GCs) to reference data, mainly derived from Hubble Space Telescope (HST) observations. The GC observations were obtained as part of a separate variability survey (Skottfelt et al. 2015a;Figuera Jaimes et al. 2016), and consist of a large number of observations at varying hour angles, dates, and atmospheric conditions for each cluster.
Star positions were extracted using the DAOPHOT program (Stetson 1987), using the best 1% of exposures. Reference stellar astrometry was taken from a number of sources, listed in Table 5. To account for pointing offsets between different observations, we cross-correlated each image with a reference image of the same cluster that was taken in good conditions and was well centred on the cluster.
We transformed the reference star positions to the measured positions using a four-parameter fit: the detector plate scale, in mas/px; the detector rotation eastwards from North; and the (x, y) coordinates of the cluster centre for each image. For each iteration in the fitting process, the reference star positions were transformed and each measured star was paired with the nearest reference star. The data were fitted using the least-squares method, with the distances between each pair of stars being minimized. We assumed that the effect of proper motions was negligible, and that any effect would average out due to the large number of reference stars.
A correction was made for atmospheric dispersion in the TCI data, which modifies the apparent position of stars, following Eq. 12 of Gubler & Tytler (1998) and utilising their slaRefco routine to determine the required coefficients. Dispersion varies slightly with meteorological conditions, but archival meteorological data for the La Silla observatory during much of 2015 is missing 3 . We therefore assumed typical atmospheric conditions of T = 10 C, P = 770 mbar, and a relative humidity of 10%. We assumed a temperature lapse rate of 0.0065K m −1 (Gubler & Tytler 1998), an observatory altitude of 2340 m and latitude of −29.264 • , and a wavelength of 800 nm for the observed light.
We inspected our astrometric results for correlations with a number of variables, and noted a variation with hour angle of up to 0.05 • , shown in Fig. 6. The Danish 1.54m telescope is equatorially mounted, and so should not suffer from field rotation. However, if the polar axis is misaligned by some small angle θ, this will result in field rotation that is most extreme towards the poles, and will vary sinusoidally with hour angle; both of these effects are seen. This effect can be modelled as the field rotation seen by an alt-az telescope situated at a latitude of 90 − θ, and we fitted a polar axis offset of θ = 0.017 ± 0.003 • . We corrected the orientation for each observation based on the telescope coordinates stored in the image FITS headers, and the final derived astrometric solutions for each cluster are given in Table 5.
For both cameras, the derived detector scales from the individual clusters show good agreement. However, systematic offsets on the order of 0.1 • were found in the calculated detector orientation between different clusters, significantly larger than the RMS scatter of ∼ 0.01 • for observations of the same cluster -the offset between the individual measurements for NGC 5139 and NGC 5286 are clearly visible in Fig. 7. The offsets in orientation are common between the cameras, with the relative angle between the two detectors being constant at 2.54 ± 0.02 • .
In some cases, multiple sets of high quality astrometry were available for the same cluster, which should result in the same astrometric solution if the reference data have been correctly calibrated. For NGC 6121, the data in Libralato et al. (2014) are matched to the UCAC4 reference frame, but agree well with the ACS Survey data. The two sets of data for NGC 104 use independently determined reference frames based on the HST guide stars, but also agree on detector orientation -however, the detector scales are significantly different, and the scale we derive from McLaughlin et al. (2006) is noticeably larger than for any other data set. For NGC 5139, Bellini et al. (2010) and Anderson & van der Marel (2010) adopted the HST-based reference frame of Anderson et al. (2008), and the results from all three sets of data reflect this common reference frame; the positions in Bellini et al. (2009b)    reference frame (Bellini et al. 2009a), and are apparently orientated by ∼ 0.1 • relative to the ACS Survey reference frame. Rotation offsets between datasets have been noted previously for HST-derived data, with Anderson & van der Marel (2010) noting a difference of 0.1 • between the ACS Survey and 2MASS for NGC 5139. The HST has been shown to match commanded roll angles to far better than 0.1 • (van der Marel et al. 2007), and if the orientation differences do result from the reference data, we consider the most likely source of errors to be the astrometric calibration itself. For most of our sets of reference data, astrometry has been placed on an absolute reference frame using the HST guide stars (Anderson et al. 2008). In normal operation, two guide stars are used by HST, with the astrometric uncertainties dominated by the accuracy of the Guide Star Catalog II (GSCII), which has an RMS scatter of ∼ 0.3 in stellar positions (Koekemoer et al. 2006). We created a simple model of FGS star position measurements, based on the FGS parameters given in Nelan et al. (2016) and the assumption that pairs of GSCII stars would be uniformly distributed across the available FGS detector space, with normally-distributed positional errors of σ = 0.3 on each star. From this model, we find an r.m.s. error of 0.025 • in the measured orientation of HST derived from two guide stars, with the range of offsets reaching up to 0.1 • .
We conclude that the systematic offsets we see are due to the limited accuracy of our HST reference data, and that these are the limiting factor in determining our detector orientation and scale. Our final astrometric results are therefore a scale of 88.91± 0.05mas/px and orientation of 1.30 ± 0.06 • eastwards of North for the red camera, and 89.07 ± 0.05mas/px and −1.25 ± 0.06 • eastwards of North for the visual camera. The fit results from Bellini et al. (2010) and Anderson & van der Marel (2010) were not used when calculating these final values, as they share the same astrometric calibration as the ACS Survey, and therefore are not independent estimates of our detector orientation.

Dichroic leak at 410nm
Light entering the TCI is split between the red and visual arms by a dichroic mirror with generally simple transmission/reflection properties, transmitting light redwards of 655nm to the red camera, and reflecting bluer light towards the visual camera. The transmission properties are less well behaved between 350 and 430nm, with the dichroic transmitting some light towards the red camera, the main transmission feature being located at 400nm. It was previously assumed that the only effect of the leak was a small contamination of any red camera photometry, which was included in the colour calibrations presented in Paper I.
As the TCI does not have an atmospheric dispersion corrector (ADC), the apparent position of a star will change with wavelength. Given a wavelength difference of approximately 400nm between the blue leak and the bulk of the light reaching the red camera, it is possible for a star to appear twice on an image. At a zenith distance of 40 • (airmass 1.3) and typical atmospheric  conditions, 4 the secondary image formed by the blue light will be separated from the main image of the star by an arcsecond, and hence is easily resolved with the TCI. The secondary image will appear at a position angle corresponding to the parallactic angle q (when measured from main image to secondary image, q + 180 • ). Fig. 8 shows two observations of the visual binary WASP-85, with both components having visible secondary images that vary in position with the observing conditions. Atmospheric dispersion also causes the main PSF of the target to be elongated, reducing our sensitivity to companions that fall along the parallactic angle.

Lucky Imaging Results
All detected companions are listed in Tables 9 and 10, covering our 2015 and 2016 observations respectively. Where targets were observed in 2014, we retain the original companion numbering system for ease of comparison. Newly detected companions are assigned numbers in order of increasing separation. In the following short paragraphs, we consider in detail companions with colours and magnitudes indicating that they are likely to be bound to the target host stars, using archival photometry and astrometry where possible. Also mentioned are companions that are particularly bright or close to the planet host stars, due to the increased fraction of contaminating light that these stars produce. In Table 6, we provide a summary for companions whose distances are consistent to within 2σ of the relevant planet host star, excluding companions where there is other available evidence (e.g. discrepant proper motions). Due to stellar crowding and high reddening resulting in many stars having apparently consistent distances, this table does not include the majority of companions to CoRoT-04, CoRoT-22, and HAT-P-45.
In considering individual companions of interest, archival data were sourced from a variety of catalogues. Where available, we used ground-based proper motions from the NOMAD (Zacharias et al. 2004), PPMXL (Roeser et al. 2010), UCAC4 (Zacharias et al. 2013), and URAT1 (Zacharias et al. 2015) catalogues to assess whether pairs of stars showed common proper motion, but we note that many of our detected companions are too faint to have reliable proper motions, if any at all. We also used archival position measurements to compare to our own, with 2MASS (Cutri et al. 2003)  The 2MASS catalogue also includes JHK photometry for some companions, although we were rarely able to derive properties for companions from 2MASS alone, as dwarf stars below ∼ 0.7M show little variation in their JHK colours, and the faintness of most detected companions resulted in large photometric uncertainties. More rarely, resolved photometry was available in the optical, mainly from the SDSS DR9 (Ahn et al. 2012) and APASS (Henden et al. 2016) catalogues. Finally, some of our targets are included in the Washington Double Star catalog (WDS, Mason et al. 2001), which includes historical astrometry and photometry, with records extending back to the 19th century in a few cases.
CoRoT-2. Companion 1 is a bright star near CoRoT-2 that was noted by Alonso et al. (2008), with several lucky imaging observations since Wöllert & Brandner 2015, and Paper I). Additionally, Schröter et al. (2011) characterised the companion spectroscopically, measuring a temperature between 3900K and 4100K (95% confidence interval), and finding its radial velocity to be consistent with CoRoT-2. We measure a somewhat lower temperature of 3660 ± 90K for the companion, and find no other companions consistent with being physically associated.
CoRoT-22. Of the companions observed in 2016, it is unsurprising that a significant portion have distances that are consistent with CoRoT-22, given its relatively large distance of nearly 600pc (Moutou et al. 2014) and location in the galactic plane. We found no companions closer than a star located 3.3 to the North, which was previously imaged by Moutou et al. (2014) using VLT/NACO. One notable companion, number 8, is located 8.4 to the Southwest, 3.9 mag fainter in r TCI , with an estimated distance that differs by only 0.2σ from that of CoRoT-22, perhaps indicating that it is a bound companion in a wide orbit.
However, due to the high probability of chance alignments, further evidence is required before any definite conclusions can be made.
CoRoT-28. From our 2015 observations, we detected a number of faint background stars within the CoRoT photometric mask that were not identified by Cabrera et al. (2015). The brightest companion falling within the photometric mask, companion 3, is 5.9 mag fainter in r TCI at a separation of 5.9 , and all other stars are at least 6.4 magnitudes fainter in the same filter. Companion 37 is the brightest star detected in our observations, at a separation of 22.5 , is the only companion with a colour consistent with being bound; however, neither the NOMAD, PP-MXL, or UCAC4 catalogues support common proper motion, and so this star is likely unassociated with CoRoT-28.
HAT-P-27. A close companion located 0.65 from HAT-P-27 was announced by  and also observed by Ngo et al. (2016). However, our observations do not sufficiently resolve these two stars due to the companion's faintness, with the companion being visible as a small distortion in the 2015 data, and obscured by the PSF elongation caused by atmospheric dispersion in the 2016 data.
HAT-P-30. Both companion stars observed in 2014 were reobserved in 2015 and 2016. Companion 1 has previously been analysed by Enoch et al. (2011), Adams et al. (2013) Ginski et al. (2013), and Ngo et al. (2015). Common proper motion was preferred by the analysis in Paper I and that of Ngo et al. (2015), and our new separation measurements are consistent with this conclusion. We derive a colour temperature of 3800 ± 130K, and the measured flux ratios in r TCI and v TCI correspond to those expected of a bound companion, and so we conclude that these two stars are physically associated.
Companion 2 is included in the WDS as a 'C' component in the system. In paper I, we noted that 2MASS J and H magnitudes for the companion made it inconsistent with being bound. Our two colour 2015 data are in agreement with this conclusion, with a colour temperature of 5700 ± 500K, the star being much too faint to be a bound companion at this temperature. In addition, the URAT1 catalogue indicates that the two stars do not exhibit common proper motion.
HAT-P-35. All four companion stars found in Paper I were reobserved in 2015 and 2016. Companion 1 is the closest at 0.93 , and has also been observed by  and Ngo et al. (2016), with the latter work confirming common proper motion, but this companion is not sufficiently resolved in either of our datasets from 2015 or 2016 for reliable astrometry. Two colour photometry indicates a temperature in the range 3700 − 4200K, somewhat higher than the values of 3525 ± 76K and 3563 ± 70K derived by Ngo et al. (2016).
For companion 2 at 9 , we derive a temperature of 3560 ± 70K for this star, with its measured flux ratios being consistent with a bound companion to HAT-P-35. Resolved SDSS and 2MASS photometry are available for this companion, with the BT-Settl model isochrones (Allard et al. 2012), indicating a M = 0.5M star (T eff = 3680K) at the same distance as HAT-P-35, using the Bakos et al. (2012) age and distance estimates for the system. Combining 2MASS and Gaia astrometry with our own, we find no significant trends in separation or position angle, though we note that HAT-P-35's proper motion is low. The NOMAD catalogue gives proper motions for the star that are consistent at about 1σ, whilst the URAT1 catalogue clearly indicates common proper motion. It therefore seems likely that this companion is also physically associated with HAT-P-35, making the system a relatively unusual example of a hierarchical triple with three resolved components.
HAT-P-41. In Paper I, we reported the detection of a new companion at 1.0 with a PA of 190 • , the magnitude difference being measured as ∆r TCI = 4.4. However, despite the good seeing in our 2015 and 2016 observations, we do not re-detect this companion; similarly, the companion was not detected in  or , despite our reported magnitude being above their 5σ detection limit.
Following the discovery that the TCI can produce secondary images of the target star, discussed in Sect. 3.5, we calculated the expected properties of a secondary image for our 2014 observations 5 . The predicted separation of 0.8 and PA of 195 • correspond closely with our 2014 detection, and we therefore conclude that this was a spurious detection.
Companion 1 at 3.6 was previously observed by , , and Ngo et al. (2016), with the latter finding the star to be consistent with a bound companion based on its colours, although their proper motion analysis was inconclusive. With a temperature of 4470 ± 120 K, we find the companion's distance is consistent with HAT-P-41 at the 1σ level.
HAT-P-45. Due to its position near the galactic plane, a large number of faint, red stars were observed near HAT-P-45 in 2015, many of which have consistent distances. It is highly likely that many, if not all of these stars are unassociated with HAT-P-45, and further evidence would be required to determine whether any are likely to be bound.
HATS-1. In 2015 we re-observed our previously detected companion 1, located at a separation of 11.5 and 8mag fainter in r TCI . The companion has a position measured in Gaia DR1, with no significant trend in separation or position angle being seen. With a colour temperature of 3130 ± 100K, the companion appears slightly fainter than expected of a star at the same distance as HATS-1 at the 2σ level.
We also detected a new companion located 6 to the South of HATS-1, and 10mag fainter in r TCI . This star falls almost directly beneath one of the diffraction spikes originating from HATS-1, and as a result our photometry is badly contaminated; whilst we do measure a very blue colour for this companion (v TCI − r TCI = −0.5 ± 0.3), we consider this unreliable.
Both of these companions are separate from the very wide common proper motion companion identified by Mugrauer et al. (2014).
HATS-2. The close companion at 1.1 found in 2014 was re-observed in 2015 and 2016. The UCAC4, NOMAD, and PP-MXL catalogues all agree that HATS-2 is moving predominantly in RA at a rate of approximately −45mas/yr, which would lead to increasing separation and position angle given the current position of the companion. Our astrometric measurements do indicate positive trends, but these are not statistically significant (0.3σ in separation, 0.6σ in position angle). Analysis of the colour of the companion indicates that it is much too faint to be physically associated, its flux ratios being 4σ too faint to match our measured temperatures of 4300 ± 300K (2015 observation) and 3900 ± 220K (2016).
HATS-3. We reobserved the close companion at 3.6 in 2016, deriving a temperature of 4000±300K, with the star therefore being 4σ too faint to be physically associated.
HATS-11. Located in a fairly crowded field, we detected 28 companion stars in an observation in 2016. Companion 1 is lo-5 Weather conditions at time of observing were retrieved via ESO's La Silla Meteorology Query Form, available at http://archive.eso. org/wdb/wdb/asm/meteo_lasilla/form Article number, page 9 of 27 A&A proofs: manuscript no. hitep2_accepted_20170921 cated at only 1.4 , and is 4.6 mag fainter in r TCI . This star is too faint at the 2σ level, with a temperature of 3920 ± 160K.
Companions 2, 3, and 7 have colours and flux ratios that give distances that are consistent with HATS-11 to within 1σ. Given HATS-11's distance of 900pc and a location near the galactic plane, the probability of chance alignments of stars with calculated distances that agree within our uncertainties is relatively high.
HATS-14. Of the five stars detected in 2016, only companion 1 at 7.8 has a colour and magnitude consistent with being a bound object at less than 2σ. Companion 3 has a diffuse, extended profile, and is likely a background galaxy.
HATS-26. Two objects were detected in our 2016 observations. Companion 1 has a thin, extended profile, which we identify as an edge-on spiral galaxy. When treated by our reduction pipeline as a star, the centre of this galaxy is measured to be ∼ 5 from HATS-26 and 7 mag fainter, with a relatively blue colour corresponding to a black-body temperature of 6700K. Visual inspection shows that the galaxy's orientation and extended profile result in emission as close as 2 to HATS-26, and with the planet being a tempting target for transmission spectroscopy (Espinoza et al. 2016), it is worth noting the presence of such a close background object.
HATS-27. Companion 5 has similar brightness to HATS-27 at a separation of 21.9 , with colours marginally consistent with being a bound object, the companion appearing to be nearer than HATS-27 at 2σ. Being well separated and of similar brightness to HATS-27, the star is resolved in various catalogues. The companion's colours in the APASS9 and 2MASS catalogues indicate a K3/4 star, with HATS-27 being F5/6; the expected difference in brightness would be on the order of 3 mag in V, much larger than the APASS9 difference of 0.90 ± 0.06 mag.
The NOMAD, PPMXL, and UCAC4 catalogues disagree significantly with each other regarding the proper motions of the two stars, but none indicate that the pair show common proper motion. With no evidence favouring the bound scenario, we conclude that this companion is an unassociated foreground star.
K2-02 (HIP 116454/EPIC 60021410). In their analysis of this system, Vanderburg et al. (2015) noted a faint companion star, with archival imaging clearly showing that the two stars have common proper motion. SDSS photometry indicated that, whilst 7 magnitudes fainter, the companion is bluer than HIP 116454, leading Vanderburg et al. (2015) to characterise the faint star as a white dwarf. Our two-colour imaging supports this hypothesis, with our derived colour temperature of 7300 ± 200K being in agreement with the SDSS-derived value of 7500±200K.
K2-21 (EPIC 206011691). Petigura et al. (2015) presented Keck NIRC2 observations in the K cont filter as part of their characterisation of the K2-21 system, finding no companion stars with ∆K < 7.5 at separations beyond 0.5 . We detect a single faint companion at 2.99 ± 0.07 , 8.8 ± 0.5 mag fainter in r TCI and undetected in v TCI . Whilst the sensitivity curve and image presented in Petigura et al. (2015) cover only the inner two arcseconds, the observations did cover the position of our newly discovered companion, with no star being detectable at this position (D. Ciardi, priv. comm.); this non-detection in K cont places limits on the companion's properties. Using the Baraffe et al. (2015) isochrones, K2-21's temperature of 4043 ± 375K (Petigura et al. 2015), and assuming our r TCI magnitude difference to be comparable to I c , we estimate K2-21 to have I c − K = 1.8mag, and hence the companion I c − K < 2.5. This corresponds to a main sequence star with T eff 3000K, which would appear much brighter if bound to K2-21; we therefore conclude that the companion is an unassociated background star. K2-24 (EPIC 203771098). Similarly to K2-21, a Keck NIRC2 observation with no detected companion stars was presented in Petigura et al. (2016), using the Br-γ filter. We detect several companions, including a faint star at 3.82 ± 0.09 . With ∆r TCI = 8.66 ± 0.18 and ∆v TCI = 8.88 ± 0.16, the companion has a relatively similar colour to K2-24; assuming the same brightness difference in Br-γ, the companion would have been near the detection limit of the observations of Petigura et al. (2016). With a colour temperature of 5400±600K, the star is clearly much too faint to be a bound main sequence companion.
K2-27 (EPIC 201546283). Montet et al. (2015) detected a companion star to K2-27 in a PHARO observation taken on 2015 February 4. The companion was measured to be at a separation of 2.98±0.05 with a position angle of 177.7 ± 0.4 • , 3.72 ± 0.06 mag fainter in K s (Montet 2017, priv. comm.). We reobserved the companion in 2016 and measured differential magnitudes of ∆v TCI = 6.97 ± 0.13 mag and ∆r TCI = 5.38 ± 0.06 mag, corresponding to a temperature of 3150 ± 80K and a photometric distance consistent with a bound companion to K2-27. Using the Baraffe et al. (2015) isochrones and our measured temperature, we find that the K s magnitude difference also supports this hypothesis.
Although not close enough to be included in our observations, we note the presence of a nearby bright star visible on DSS images, EPIC 201546361, at a separation of 32 to the West. The PPMXL and UCAC4 catalogues show the two stars to be moving southwards, with their proper motions consistent well within 1σ. The NOMAD catalogue gives a similar direction of motion, but with a larger 2σ discrepancy between the individual motions. The URAT1 catalogue also supports common proper motion, but disagrees significantly with all other catalogues as to the direction of motion, showing the pair of stars to be moving to the Southeast; it is not clear why these values are so discrepant.
As the latter pair of stars are both bright and are wellseparated from each other, photometry is available in a number of literature surveys. To determine temperatures and distances for the two stars, we fitted the ATLAS9 model atmospheres (Castelli & Kurucz 2004) to the literature photometry the two stars. We used the GALEX near-UV magnitude (Bianchi et al. 2012); APASS B, V, r, and i magnitudes (Henden et al. 2016); 2MASS J, H, and K s magnitudes (Skrutskie et al. 2006); and WISE W1 and W2 magnitudes (Cutri et al. 2012) in our fit. We fixed log g = 4.5 for the fits, and zero reddening was assumed. The effect of the faint companion to K2-27 was neglected, the flux contribution being only 6% in K s and less at shorter wavelengths.
For K2-27, we determined a temperature of 5250 ± 130K, consistent with the literature value of 5320 ± 70K , and a distance of 236 ± 19pc. For EPIC 201546361 we derive a temperature of 5880 ± 140K and distance of 220 ± 15pc. As the distances to the two stars are consistent with one another, as well as common proper motion being favoured by all catalogues we inspected, we conclude that it is likely that these two bright stars are physically associated. This hypothesis will be easily tested with future Gaia data releases, which will provide accurate parallaxes and proper motions for these two bright stars.
K2-31 (EPIC 204129699). We detected a faint star located 8.4 from K2-31 in our 2016 observations, 7.8 mag fainter in r TCI . The companion is included in both 2MASS and Gaia DR1, and there is a highly significant trend of increasing separation over time, leading us to conclude that this companion is likely to be a background star.
K2-44 (EPIC 201295312). A companion star located 8 arcsec. from K2-44 is included in numerous catalogues, and was observed by Crossfield et al. (2016) using NIRC2 on Keck and PHARO on the Palomar 200 inch on unspecified dates, with separations and relative magnitudes in K being reported. We observed this companion in 2016, with its colours and flux ratios indicating that it is likely a foreground object. Combining our astrometry with measurements from 2MASS and Gaia DR1, there are significant trends of increasing separation (8.3σ) and position angle (3.7σ), indicating that the two stars are unassociated. This conclusion is further supported by the URAT1 catalogue, which gives proper motions for the two stars that are inconsistent at the 4σ level.
KELT-15. Of the nearby stars detected in 2016, only companion 1 is consistent with being bound to KELT-15. Located 6.13 away at a position angle of 283 • , we determine a temperature of 3210 ± 70K which is consistent with a bound object at 1σ, the companion being somewhat brighter than expected of a bound companion.
WASP-8. The faint, close companion to WASP-8 was noted by Queloz et al. (2010), with measurements of the pair being recorded in the WDS as far back as 1930. The companion was detected in Paper I, and also imaged by Ngo et al. (2015), with all observations indicating a bound companion. We derive a temperature of 3550±80K for the companion, somewhat cooler than our Paper I temperature of 3740 ± 100K, but in good agreement with the 3380 − 3670K range quoted by Ngo et al. (2015).
WASP-26. Smalley et al. (2010) noted a bright visual companion to WASP-26 at a separation of 15 , and using archival photometry found the star to have a distance consistent with WASP-26. Furthermore, no relative motion was found when comparing the positions of the two stars between the 1950s Palomar Observatory Sky Survey (POSS) and the 2MASS survey, indicating common proper motion. Our two-colour photometry supports the conclusion that the two stars have a common distance, with our temperature of 4510 ± 210K being in agreement with the Smalley et al. (2010) determination of 4600 ± 120K. We find no evidence of a trend in separation or position angle, and also note that the PPMXL catalogue gives proper motions for the two stars that agree at the 2σ level.
WASP-36. Four companion stars to WASP-36 were noted by Smith et al. (2012), who concluded that none of the stars were physically associated based on proper motion catalogues. However, we find that the measurements in ground-based NOMAD, PPMXL, and UCAC4 catalogues suffer from significant confusion between the sources, with multiple spurious entries and no consistent proper motion measurements for any of the companions, and no proper motion determination at all for the closest companion. The 2MASS catalogue contains resolved JHK measurements for the companions, which are all consistent with temperatures between 2000K and 4000K -more accurate determination is not possible due to the large photometric uncertainties, and the small change in JHK colours with temperature in this regime.
Two-colour imaging indicates that companion 4 is likely to be a foreground object, whilst companion 2 appears to be a background star; the remaining companions, 1 and 3, both have colours that are consistent with being physically associated with WASP-36 within 1σ. Whilst companion 1 has also been observed by  and Ngo et al. (2016), it is not yet possible to draw any conclusions from the common proper motion analysis.
WASP-45. We reobserved the previously detected companion star in 2016, and find that its two-colour photometry is consistent with a bound object, with a temperature of 3020 ± 90K. Astrometry for the pair exists in Gaia DR1, with WASP-45 also having a measured proper motion based on the Tycho-Gaia Astrometric Solution (TGAS), moving at 70mas/yr at a position angle of 131 • . If the companion were a stationary background star, a change of 0.14 would be expected between our 2014 and 2016 observations -this is not the case, with no significant trend present in either separation or position angle. We therefore conclude that these two stars form a physical binary.
WASP-49. A close companion star at 2.2 was discovered by Lendl et al. (2016), and also detected in our 2014 observations. Lendl et al. (2016) also obtained a spectrum of the companion relative to WASP-49, finding it to be redder than WASP-49, but not characterising it further. We derive a colour temperature of 3230 ± 100K for the companion, with the measured flux ratios being slightly inconsistent with a bound object at 1σ, the companion being brighter than expected. Gaia DR1 resolves the two stars, including TGAS proper motions for WASP-49, but with relatively high uncertainties on the position of the companion. There is a slight decreasing trend in separation and increasing trend in position angle, as expected of a stationary star given WASP-49's motion, but these trends are not statistically significant; the situation for this companion is therefore inconclusive. From two-colour photometry, we conclude that all other companion stars, including that at 9 reported by Lendl et al. (2012), are background objects.
WASP-54. We re-observed a companion star at 5.7 in both 2015 and 2016. We determine a temperature of 3000±110K, and find that it is consistent with being a bound object. No significant trend is seen in separation or proper motion, but more measurements are required to definitively confirm common proper motion.
WASP-55. Companion 1 at 4.3 was re-observed in 2015 and 2016, with two-colour photometry from both years indicating a temperature of 3340 ± 90K. Our measured flux ratios are consistent with a physically associated companion at 1σ, the companion being slightly brighter than expected. WASP-55's proper motion is too low to confirm common proper motion from the available data, but no significant trends are seen in our astrometric measurements.
WASP-70 and WASP-77. The WASP-70 and WASP-77 systems both have bright, nearby companions with similar properties to the planet host stars. As in Paper I, we support the conclusion that WASP-70AB and WASP-77AB are binary systems. From our two-colour photometry of both B components, we measure temperatures that are slightly lower than the literature values. For WASP-70B, we derive a temperature of 4510±230K, compared to the previous measurement of 4900 ± 200K by Anderson et al. (2014b); similarly, we measure 4570 ± 240K for WASP-77B, compared to a previous value of 4700 ± 100K .
WASP-85. The stars in the WASP-85AB system were characterised by Brown (2015), measuring a spectroscopic temperature of 5250 ± 90K for the secondary star, in good agreement with the weighted mean of our 2015 and 2016 measurements of 5340 ± 340K. The companion has since been observed by , Tokovinin et al. (2016), andSchmitt et al. (2016), although we note that the latter two papers are separate reductions of the same data. Combining these measurements with historical data recorded in the WDS, orbital motion is clearly seen, with a 15 • change in position angle since 1880, although no clear trend is visible in separation. The historic data are not sufficiently precise to allow us to place any meaningful constraints on the binary orbit, and further high-precision astrometry is required to fully determine the orbit. WASP-87. Anderson et al. (2014a) noted a bright star 8.2 to the Southeast of WASP-87, and suggested that the two stars may form a bound WASP-87AB system, based on resolved 2MASS photometry indicating a temperature of 5700 ± 150K for the secondary. It was also stated that the two stars had "similar" proper motions in the UCAC4 catalogue, but we note that the differences are significantly larger than the quoted uncertainties. Additionally, there are significant differences between the proper motions quoted by the UCAC4, PPMXL (primary star only), and NOMAD catalogues, likely due to confusion between sources in a crowded field.
Based on our observations in 2015 and 2016, we find that WASP-87B appears slightly too faint compared to WASP-87A based on its colour temperature. However, our analysis does not account for the fact that WASP-87A is somewhat evolved, and hence appears brighter than predicted by our temperature-radius relationship -we therefore conclude that these two stars are consistent with a bound pair.
WASP-94. WASP-94 is another system with a bright companion, with the addition of a non-transiting planet orbiting WASP-94B. Our measured temperature of 6400 ± 500K for the B component is in agreement with the Neveu-VanMalle et al.
(2014) determination of 6040 ± 90K. Teske et al. (2016) measured elemental abundances of the two stars, claiming that significant differences were seen, despite both stars having similar stellar properties and orbiting giant planets. Under the assumption that the differences were caused by planet formation, the authors hypothesised that additional stellar companions may have affected the planet formation process. Based on our sensitivity limits given in Table 3 and the temperatures of WASP-94 A and B (6198K and 6112K respectively, citealt2016ApJ...819...19T) we exclude bound companion stars hoter than 3800K beyond 0.8 , 3300K beyond 1.5 , and 2900K beyond 2.5 around both stars. Using the mass-temperature relations presented in Paper I, these correspond to companion masses of 0.59M , 0.41M , and 0.23M respectively.
WASP-104. A companion star at 6.8 was detected in both 2015 and 2016, although data were only obtained in r TCI in 2015 due to a technical fault. The companion is found to be slightly too bright in our 2016 photometry, but the v TCI measurement is significantly contaminated by a diffraction spike. Resolved photometry for the companion is available in 2MASS and SDSS, with the colour-spectral type relations presented by Hawley et al. (2002) indicating a spectral type of ∼M6 for the companion. Using the BT-Settl model isochrones (Allard et al. 2012), we find that the SDSS riz photometry is consistent with the companion being a bound companion. Further supporting this hypothesis, the URAT1 catalogue indicates common proper motion, although with large uncertainties for the motion of both stars.
WASP-108. Located in a crowded field, we identify two stars with colours consistent with being bound to WASP-108. Companion 10 is only 2 mag fainter in r TCI , but the UCAC4, NO-MAD, and PPMXL catalogues all indicate that this companion does not have common proper motion with WASP-108, and so we conclude that it is an unassociated star.
For companion 4 at 8.8 , we measure a temperature of 2950 ± 70K in 2016, and find its flux ratios to be consistent with a bound object. However, a chance alignment is relatively likely given the crowding of the region and the faintness of this star, and further evidence is required to confirm whether these two stars are physically associated.
WASP-110. A bright companion to WASP-110 was noted by Anderson et al. (2014a), who measured a separation of 4.589 ± 0.016 and a magnitude difference of 2.872±0.012 mag in I c . We measure a similar magnitude difference in r TCI of 2.284 ± 0.006 mag, and a smaller difference in v TCI of 2.298 ± 0.004 mag, indicating that the companion is hotter than WASP-110 with T eff 8400K. Whilst not resolved in most catalogues, visual inspection of multicolour DSS and 2MASS images supports the hypothesis that the two stars have similar colours.
Various catalogues indicate that WASP-110 has a high proper motion, but with a significant scatter in the actual values, likely due to confusion between the two unresolved stars. The UCAC4 catalogue does resolve the pair, giving significantly different proper motions for the two, indicating that they are not physically associated. Our astrometry, combined with Gaia DR1 positions and the separation measured by Anderson et al. (2014a), indicate trends in both separation (increasing) and position angle (decreasing), in accordance with the UCAC4 proper motions. We therefore conclude that these two stars are merely a chance alignment.
WASP-111. We detected a bright companion 5.0 from WASP-111 in both 2015 and 2016, with its colour and flux ratios being consistent with a bound companion with T eff = 3710 ± 60K. 2MASS JHK photometry and Tycho/YB6 B photometry recorded in NOMAD also support this conclusion. The NOMAD catalogue indicates that the two stars have common proper motion, albeit with large uncertainties on the companion's motion. Our astrometry, combined with measurements from the 2MASS and Gaia DR1 catalogues, show no significant trends in separation or position angle since 2000. We conclude that these two stars form a bound system. WASP-121. Companion 1 was detected in 2016, at a separation of 7.6 . Its red colour is consistent with a bound companion, being 7.5 mag fainter in r TCI and 9.6 mag fainter in v TCI . This star is included in the 2MASS catalogue, in which it is ∼ 5 mag fainter in JHK, further indicating that it is a cool object, but the 2MASS photometry is too uncertain to classify the object further. Comparing the 2MASS positions to our own astrometry, we find no large change in position angle or separation, but further astrometric observations are required to confirm whether common proper motion is present.
WASP-123. We detected a companion at 4.8 in 2016, and derived a temperature of 3240 ± 80K, the companion being marginally consistent with a bound companion. Archival 2MASS astrometry indicates that no significant change in separation or position angle has occurred, but the 2MASS position is too uncertain to confirm common proper motion.
WASP-129. We find companion 2, located 5.9 to the West of WASP-129, to be consistent with a bound companion, being 5.6 mag fainter in r TCI and a photometric temperature of 3360 ± 90K. 2MASS astrometry is consistent with no change in the relative position of the companion, but the proper motion of WASP-129 is too low to determine if the two stars show common proper motion.
WASP-133. A very close companion was detected in 2016, separated from the target star by only 2.4 , ∼ 6.7 mag fainter. From its colour, we determine a temperature of 5000 ± 400K, indicating that it is only ∼ 1000K cooler than WASP-133. As a result, we conclude that this companion is likely to be a distant background object in a chance alignment.

Orbits of binary companions
For a handful of wide binaries hosting transiting hot Jupiters, historical astrometry of the stellar binary is present in the Washington Double Star Catalog (WDS). For WASP-77AB and WASP-85AB, orbital motion is clearly present, but the data cover only a small fraction of the orbit. This results in loosely constrained (and often multi-modal) orbital parameters that are difficult for Markov Chain Monte Carlo algorithms to explore. Instead, we adopt the Bayesian rejection-sampling algorithm detailed by Blunt et al. (2017), which has been shown to converge significantly faster in such cases. For each system, the algorithm was run until 1,000,000 orbital fits had been accepted.
The WDS does not record uncertainties for the measured astrometry, and only for a few of the most recent results are uncertainties quoted in the original literature. Where we were unable to find any measurements, we assumed uncertainties of 0.1 in separation and 1 • in position angle, these values being chosen by visual inspection of the scatter between measurements.
For both systems, radial velocity measurements were obtained for both of the stellar components as part of the planet follow-up programs, in order to confirm which of the stars exhibited radial velocity variations and hence hosts a planet. Compared to the long orbital period of the binary, these measurements provide only a snapshot of the radial velocity of the system, but do allow us to break the degeneracy in the longitude of ascending node, Ω, and hence constrain the full 3 dimensional motion of the system. However, without significant phase coverage it is not possible to include as a free parameter the systematic shifts in measured radial velocity caused by gravitational redshift and convective blueshift, as has been done previously in the case of α Cen (Pourbaix & Boffin 2016;Kervella et al. 2017). These effects vary on a line-by-line basis, currently require full 3D hydrodynamical modelling to match Solar observations, and also depend on both instrumental effects and the reduction method used (Allende Prieto et al. 2013); given the variety of instruments and reduction pipelines used to derive the available radial velocities, we have instead opted to account for these effects with an increased uncertainty of 0.5km/s. We show the fitted orbits on the RA-Dec plane and show the mean, median, etc. of the resulting parameter distributions in the following two subsections. Additional tables containing the astrometric data used, figures showing the predicted separation/position angle, and figures showing the correlations between the parameters can be found in Appendix A.

WASP-77
WASP-77AB was observed numerous times between 1903 and 1933, with the data indicating a change in position angle of ∼ 2 • and an increase in separation of ∼ 0.3 . In addition, radial velocity measurements of the individual A and B components were presented by Maxted et al. (2013), from which we derive a radial velocity difference of 1.1 ± 0.5km/s at a mean BJD of 2455847. We constrain the total mass of the system using the stellar parameters from Maxted et al. (2013) assuming the errors on each component's mass are independent, giving M Tot = 1.71± 0.07M .
Our resulting parameter distributions are summarised in Table 7. A selection of 100 orbits drawn from our accepted fits are plotted in Fig. 9. We find that the binary orbit has a semi-major axis of 461 +200 −140 au, is nearly edge-on (i = 77 +5 −7 • ), and is moderately eccentric (e = 0.51 +0.26 −0.22 ). The constraint provided by the radial velocity measurements is sufficient to break the degeneracy in Ω.

WASP-85
WASP-85AB has been regularly observed since 1881. Historic measurements show a large scatter in separation, but a clear trend is seen in position angle, from approximately 113 • to modern measurements near 100 • . Brown (2015) obtained radial velocity measurements of the two individual stars using HARPS in good seeing, from which we derive a radial velocity difference of −0.5 ± 0.5km/s at a mean BJD of 2456024. We constrain the total mass of the system using the stellar parameters from Brown (2015) assuming the errors on each component's mass are independent, giving M Tot = 1.92 ± 0.10M . Our resulting parameter distributions are summarised in Table 8. A selection of 100 orbits drawn from our accepted fits are plotted in Fig. 10. We derive a = 148 +52 −23 au, e = 0.43 +0.13 −0.25 and i = 140 +16 −12 • . However, we note that these parameters are strongly correlated, with smaller semi-major axes corresponding to more eccentric orbits that are closer to face on, and that the distributions are clearly non-gaussian, particularly for a and e. The radial velocity constraint allows us to prefer one of the orbital orientations, but is not sufficient to reject the other.
When considering dynamical interactions between the planet and the outer star, one of the most important parameters is the relative inclination of the two orbits. For both the transiting planet and resolved binary, the orbital inclinations i pl and i bin relative to the line of sight is easily determined, and as mentioned previously, the longitude of ascending node of a visual binary orbit Ω bin can be determined if radial velocity measurements are available. However, Ω pl is entirely unconstrained, as the orientation of a transiting planet's orbit with respect to the plane of the sky is unknown unless the planet is directly imaged, which is not possible for hot Jupiters. Therefore, in order to derive the distribution of the relative inclination between the planetary and binary orbits from our data, we assume a uniform prior of 0 • ≤ Ω pl < 360 • .

Application to planet migration theory
The proposed eccentric Lidov-Kozai mechanism of hot Jupiter migration invokes a scenario in which an inclined wide binary companion forces an initially 'cold' planet into a highly eccentric orbit, which then shrinks through tidal friction during close encounters with the planet host star (see Naoz 2016 for a review of the topic). One well-known prediction of this theory is that the migration will induce a spin-orbit misalignment between the planet and its host star, by changing the orientation of the planet's orbital plane. Theoretical work has also predicted the distribution of inclination angles between inner (planetary) and outer (binary) orbital planes, with numerical studies indicating that the distribution peaks for inclination angles around 40 − 60 • and 120 − 140 • (Fabrycky & Tremaine 2007;Wu et al. 2007;Naoz et al. 2012). For WASP-77, the projected inclination between the two orbital planes is 22±16 • , and for the WASP-85 system it is 60±17 • , adopting the planetary orbital inclinations from Maxted et al. (2013) and Močnik et al. (2016) respectively. Determining the true inclination between the orbits is complicated by our lack of knowledge regarding the true orientation of the planetary orbit, as it is impossible to determine Ω p without astrometric measurements of the orbit or direct imaging of the planet. We estimate the true inclination angle using a Monte Carlo simulation, choosing 10 6 values for Ω p from a uniform distribution (0 • ≤ Ω p < 360 • ), with our resulting 1-σ confidence limits being 36 • < i < 144 • for WASP-77 and 59 • < i < 121 • for WASP-85.
We briefly note that many analyses of the Lidov-Kozai effect assume a circular orbit for the outer (stellar) binary. However, it has been shown that significantly different behaviour can result from eccentric outer orbits (Naoz et al. 2011(Naoz et al. , 2013 -in particular, the Lidov-Kozai effect is able to occur in many more orbital configurations for eccentric systems than for circular. Our fitted binary orbits do indeed have non-negligible eccentricity, which is in line with results from studies of field stars, which show that many wide stellar binaries are eccentric (Raghavan et al. 2010)

Discussion
We find good agreement with other high-resolution imaging studies that also observed our targets, with Paper I discussing a small number of faint or close-in companions that were discovered by others, but fell below our detection limits. Since Paper I, we also reported on the discovery of a close-in companion to WASP-20 at a separation of 0.26 , which fell below the sensitivity limit of our lucky imaging observations (Evans et al. 2016b). For the targets newly added to our survey since Paper I, no directly imaged companions went undetected in our data. However, this is at least partly due to the fact that no other survey has yet covered many of these these recently announced TEP systems, with the recent work of Ngo et al. (2016) mainly covering systems announced before 2013.
A small number of new targets in our survey have direct imaging observations reported in an individual basis, with some AO or lucky imaging observations being performed during the follow-up of KELT and HAT-South targets, whilst K2 targets are regularly inspected with adaptive optics imaging. Our observations have detected a handful of very faint companions to K2 targets that were at or below the quoted sensitivity limits of the AO observations, with the results otherwise being comparable; comments on individual systems are given in the previous section.

The lack of bright companions to hot Jupiters
In the reduction of our data, it was noted that among the companions that were likely to be bound, there is a deficit of bright stars compared to faint ones: of the 45 such companions listed in Table 6, none are brighter than the planet host star, only two are less than one mag fainter (4 ± 3%), and only five are less than two mag fainter (11 ± 5%). Across our range of target stars, a binary companion 2 mag fainter in r TCI approximately corresponds to a binary mass ratio q of 0.75, and so 11 ± 5% of hot Jupiter binaries detected have q > 0.75. This contrasts with the results of recent surveys of field stars, where the distribution of mass ratios (and hence brightness ratios) shows no preference towards low values, with the results of Raghavan et al. (2010) finding that 31 ± 5% of field binaries have q > 0.75. If this deficit of bright, high-mass companions does indeed exist in the population of hot Jupiter hosting systems, it would indicate that the formation of these close-in planets is inhibited in such binaries systems. This may be due to the suppression of planet formation in protoplanetary discs (e.g. Jang-Condell 2015), or instead due to their influence on planetary migration.
This deficit was also seen in the Friends of Hot Jupiters survey (FoHJ, Ngo et al. 2015Ngo et al. , 2016, with all but 2 of their 35 candidate binary companions being at least 2.5 mag fainter in the K band, but the opposite result was found by Wang et al. (2015) in a sample of hot Jupiters from the Kepler mission, with six out of their seven detected companion stars being within 2.5 mag of the planet host star. Wang et al. (2015) hypothesised that their high fraction of bright companions could be a result of the migration mechanism of hot Jupiters, with higher mass stellar companions being more easily able to migrate planets inwards via the Kozai-Lidov mechanism. The lack of bright companions to hot Jupiters discovered by ground-based planet surveys was attributed to selection biases, with the detection of such planets being inhibited by photometric dilution, combined with cautious D. F. Evans et al.: High-resolution Imaging of Transiting Extrasolar Planetary systems (HITEP).
follow-up strategies that avoid systems showing any potential evidence of binary nature, such as multiple sets of spectral lines (Triaud et al. 2017). However, due to the qualitative nature of target selection, there are no detailed assessments of the completeness of ground-based planet surveys, and hence we cannot assess the likely number of hot Jupiters in wide stellar binaries that have been missed by such surveys.
On the other hand, even if the majority of hot Jupiters do live in high mass binaries as suggested by Wang et al. (2015), there is still a notable population of hot Jupiters for which no stellar companion has been identified. The FoHJ survey found that 49% 6 of systems in their sample have no binary stellar companion . This suggests that some hot Jupiters formed in stellar binaries that have since been disrupted, leaving only a single star, or alternatively that multiple hot Jupiter formation mechanisms operate, as has been proposed by Nelson et al. (2017).

Companion detection using Gaia
Future detection and confirmation of well-separated binary companions will be made significantly easier with Gaia data. Companions of moderate brightness (G 18) at typical distances of a few hundred parsecs are expected to have parallaxes measured to better than 10% accuracy (de Bruijne et al. 2014), sufficient to rule out the vast majority of background stars present in our images. Proper motions for even the faintest companions are expected to have errors less than 1 mas/yr, based on the expected parallax-proper motion error relations (Gaia Collaboration et al. 2016), sufficient to show common proper motion for all but the slowest moving planet host stars.
Gaia's ability to directly resolve close companions falls below separations of a couple of arcseconds, especially for faint companions; however, the expected limits are similar to the sensitivity achieved by our ground-based lucky imaging survey. Gaia will reliably resolve equal-brightness components at separations greater than 0.6 (95% completeness), but successful detections of companions 6-7 magnitudes fainter will only occur in approximately 5% of cases at the same separation (de Bruijne et al. 2015).
Astrometric detections of orbital motion by Gaia will allow some unresolved stellar companions to be detected, with the limits from exoplanet detection studies suggesting that pairs of solar-mass stars could be detected out to orbits of a few tens of au (Sozzetti 2010) for systems within 200pc. Combined with resolved detections of companions, Gaia should provide good coverage for bright companions to the closest transiting exoplanet systems, where the detection limits of the two methods overlap. However, the detection of faint, close-in binary companions will remain possible only with techniques such as adaptive optics and interferometry -for example, none of the companions detected by Kraus et al. (2016) using the non-redundant aperture masking technique are expected to be directly resolved by Gaia, and the majority are thought to be in long period orbits that will not produce a detectable astrometric signal.

Conclusions
We present lucky imaging observations of 97 transiting planet systems. We analyse these data to create a catalogue of stars within the vicinity of our targets, providing relative astrometry and photometry for each companion. Where two colour photometry was possible, we estimate the stellar temperature of each detected star, and hence determine whether it is consistent with being a bound main sequence companion to the planetary system. For viable candidates, we review evidence available in the literature and public catalogues to either support or refute the 'bound companion' hypothesis. Combining our data with historical astrometry, we fit the orbits of the wide stellar binaries WASP-77AB and WASP-85AB, and find both binaries to be moderately eccentric with orbital planes that are at least slightly inclined relative to the planetary orbital plane.
We consider the mass ratio distribution of our plausible binary star systems, in addition to the sample of confirmed binaries from Ngo et al. (2015Ngo et al. ( , 2016, and find that our targets very rarely have high mass stellar companions (q > 0.75). Due to the complex selection biases present in ground-based planet surveys, it is difficult to determine if our target sample is a true representation of all hot Jupiters. However, the significant proportion of hot Jupiters with no detected stellar companion suggests that either formation via the binary Lidov-Kozai mechanism does not occur in all hot Jupiter systems, or instead that many stellar binaries were catastrophically disrupted at an earlier stage, leaving the hot Jupiter in a non-multiple stellar system. A&A proofs: manuscript no. hitep2_accepted_20170921 Notes. Companions with distances consistent at 2σ or better are included, except in the case of significant evidence against physical association (e.g. discrepant proper motions). No companions were found to have discrepant distance measurements from literature photometry or spectroscopy. We also exclude most companions to CoRoT-04, CoRoT-22, HAT-P-45, and HATS-11 from this table, due to a high fraction of background stars being identified with consistent distances.  Notes. * The eccentricity distribution is bimodal, with moderate eccentricity and high eccentricity solutions (e 0.5 and e 0.95). Table 8. Characteristics of the probability distributions for the orbital parameters of WASP-85AB. The 16th, 50th, and 84th percentiles are denoted by P 16 , P 50 and P 84 . As the values of the longitude of the ascending node (Ω) and argument of periapsis (ω) are not fully constrained by the radial velocity measurements, we present these values both over the full range (0-360 • ) and the restricted range usually quoted for astrometric orbits (0-180 • ). Notes. * The probability distribution for ω is multimodal and is highly correlated with other parameters Article number, page 19 of 27 A&A proofs: manuscript no. hitep2_accepted_20170921 Table 9. Astrometry, photometry, and derived values for our 2015 observations. The column headed σ dist gives the probability that the companion is at the same distance as the associated exoplanet host star, which is calculated by comparing the observed flux ratio to the flux ratio expected if the two stars were bound, based on their temperatures. Our flux ratio prediction is based on the assumption that both stars are on the main sequence. The full version of   Table 10. Astrometry, photometry, and derived values for our 2016 observations. The column headed σ dist gives the probability that the companion is at the same distance as the associated exoplanet host star, which is calculated by comparing the observed flux ratio to the flux ratio expected if the two stars were bound, based on their temperatures. Our flux ratio prediction is based on the assumption that both stars are on the main sequence. The full version of Pos. Angle (deg.)   . Density distributions of the parameters of the 1,000,000 accepted fits to the WASP-85AB binary orbit, based on the model described in Sect. 5. The consecutive contours encompass 68.3%, 95.5%, 99.7% and 99.9% of the samples. The symbols used are defined as: a -semi-major axis; e -eccentricity; i -inclination; ω -longitude of the ascending node; Ω -argument of periapsis; M Tot -total mass of the system. Due to a lack of a strong radial velocity constraint on ω and Ω, we have limited these two values to the range 0 ≤ x < 180.