Enceladus is believed to have a saltwater global ocean with a mean depth of at least 30~km, heated from below at the ocean-core interface and cooled at the top, where the ocean loses heat to the icy lithosphere above. This scenario suggests an important role for vertical convection to influence the interior properties and circulation of Enceladus' ocean. Additionally, the ice shell that encompasses the ocean has dramatic meridional thickness variations that, in steady state, must be sustained against processes acting to remove these ice thickness anomalies. One mechanism that would maintain variations in the ice shell thickness involves spatially-separated regions of freezing and melting at the ocean-ice interface. Here, we use an idealized, dynamical ocean model forced by an observationally-guided density forcing at the ocean-ice interface to argue that Enceladus' interior ocean should support a meridional overturning circulation. This circulation establishes an interior density structure that is more complex than in studies that have focused only on convection, including a shallow freshwater lens in the polar regions. Spatially-separated sites of ice formation and melt enable Enceladus to sustain a significant vertical and horizontal stratification, which impacts interior heat transport, and is critical for understanding the relationship between a global ocean and the planetary energy budget. The presence of low salinity layers near the polar ocean-ice interface also influences whether samples measured from the plumes are representative of the global ocean.

Stellar activity is the major roadblock on the path to finding true Earth-analogue planets with the Doppler technique. Thus, identifying new indicators that better trace magnetic activity (i.e. faculae and spots) is crucial to aid in disentangling these signals from that of a planet's Doppler wobble. In this work, we investigate activity related features as seen in disk-integrated spectra from the HARPS-N solar telescope. We divide high-activity spectral echelle orders by low-activity master templates (as defined using both log R'HK and images from the Solar Dynamics Observatory, SDO), creating "relative spectra". With resolved images of the surface of the Sun (via SDO), the faculae and spot filling factors can be calculated, giving a measure of activity independent of, and in addition to, log R'HK. We find pseudo-emission (and pseudo-absorption) features in the relative spectra that are similar to those reported in our previous work on alpha Cen B. In alpha Cen B, the features are shown to correlate better to changes in faculae filling factor than spot filling factor. In this work we more confidently identify changes in faculae coverage of the visible hemisphere of the Sun as the source of features produced in the relative spectra. Finally, we produce trailed spectra to observe the RV component of the features, which show that the features move in a redward direction as one would expect when tracking active regions rotating on the surface of a star.

We report the discovery, by the Next Generation Transit Survey (NGTS), of two hot-Jupiters NGTS-8b and NGTS-9b. These orbit a V = 13.68 K0V star (Teff = 5241 +/- 50 K) with a period of 2.49970 days, and a V = 12.80 F8V star (Teff = 6330 +/- 130 K) in 4.43527 days, respectively. The transits were independently verified by follow-up photometric observations with the SAAO 1.0-m and Euler telescopes, and we report on the planetary parameters using HARPS, FEROS and CORALIE radial velocities. NGTS-8b has a mass, 0.93 +0.04 -0.03 MJ and a radius, 1.09 +/- 0.03 RJ similar to Jupiter, resulting in a density of 0.89 +0.08 -0.07 g cm-3. This is in contrast to NGTS-9b, which has a mass of 2.90 +/- 0.17 MJ and a radius of 1.07 +/- 0.06 RJ , resulting in a much greater density of 2.93 +0.53 -0.49 g cm-3. Statistically, the planetary parameters put both objects in the regime where they would be expected to exhibit larger than predicted radii. However, we find that their radii are in agreement with predictions by theoretical non-inflated models.

Using the first 18 months of the Pan-STARRS 1 survey we have identified 33 candidate high-amplitude objects for follow-up observations and carried out observations of 22 asteroids. 4 of the observed objects were found to have observed amplitude $A_{obs}\geq 1.0$ mag. We find that these high amplitude objects are most simply explained by single rubble pile objects with some density-dependent internal strength, allowing them to resist mass shedding even at their highly elongated shapes. 3 further objects although below the cut-off for 'high-amplitude' had a combination of elongation and rotation period which also may require internal cohesive strength, depending on the density of the body. We find that none of the 'high-amplitude asteroids' identified here require any unusual cohesive strengths to resist rotational fission. 3 asteroids were sufficiently observed to allow for shape and spin pole models to be determined through light curve inversion. 45864 was determined to have retrograde rotation with spin pole axes $\lambda=218\pm 10^{\circ}, \beta=-82\pm 5^{\circ}$ and asteroid 206167 was found to have best fit spin pole axes $\lambda= 57 \pm 5^{\circ}$, $\beta=-67 \pm 5^{\circ}$. An additional object not initially measured with $A_{obs}>1.0$ mag, 49257, was determined to have a shape model which does suggest a high-amplitude object. Its spin pole axes were best fit for values $\lambda=112\pm 6^{\circ}, \beta=6\pm 5^{\circ}$. In the course of this project to date no large super-fast rotators ($P_{rot} < 2.2$ h) have been identified.

We report the discovery of NGTS-2b, an inflated hot-Jupiter transiting a bright F5V star (2MASS J14202949-3112074; $T_{\rm eff}$=$6478^{+94}_{-89}$ K), discovered as part of the Next Generation Transit Survey (NGTS). The planet is in a P=4.51 day orbit with mass $0.74^{+0.13}_{-0.12}$ M$_{J}$, radius $1.595^{+0.047}_{-0.045}$ R$_{J}$ and density $0.226^{+0.040}_{-0.038}$ g cm$^{-3}$; therefore one of the lowest density exoplanets currently known. With a relatively deep 1.0% transit around a bright V=10.96 host star, NGTS-2b is a prime target for probing giant planet composition via atmospheric transmission spectroscopy. The rapid rotation ($vsin$i=$15.2\pm0.8$ km s$^{-1}$) also makes this system an excellent candidate for Rossiter-McLaughlin follow-up observations, to measure the sky-projected stellar obliquity. NGTS-2b was confirmed without the need for follow-up photometry, due to the high precision of the NGTS photometry.

We have discovered a new, near-equal mass, eclipsing M dwarf binary from the Next Generation Transit Survey. This system is only one of 3 field age ($>$ 1 Gyr), late M dwarf eclipsing binaries known, and has a period of 1.74774 days, similar to that of CM~Dra and KOI126. Modelling of the eclipses and radial velocities shows that the component masses are $M_{\rm pri}$=0.17391$^{+0.00153}_{0.00099}$ $M_{\odot}$, $M_{\rm sec}$=0.17418$^{+0.00193}_{-0.00059}$ $M_{\odot}$; radii are $R_{\rm pri}$=0.2045$^{+0.0038}_{-0.0058}$ $R_{\odot}$, $R_{\rm sec}$=0.2168$^{+0.0047}_{-0.0048}$ $R_{\odot}$. The effective temperatures are $T_{\rm pri} = 2995\,^{+85}_{-105}$ K and $T_{\rm sec} = 2997\,^{+66}_{-101}$ K, consistent with M5 dwarfs and broadly consistent with main sequence models. This pair represents a valuable addition which can be used to constrain the mass-radius relation at the low mass end of the stellar sequence.

We present high cadence detections of two superflares from a bright G8 star (V = 11.56) with the Next Generation Transit Survey (NGTS). We improve upon previous superflare detections by resolving the flare rise and peak, allowing us to fit a solar flare inspired model without the need for arbitrary break points between rise and decay. Our data also enables us to identify substructure in the flares. From changing starspot modulation in the NGTS data we detect a stellar rotation period of 59 hours, along with evidence for differential rotation. We combine this rotation period with the observed \textit{ROSAT} X-ray flux to determine that the star's X-ray activity is saturated. We calculate the flare bolometric energies as $5.4^{+0.8}_{-0.7}\times10^{34}$ and $2.6^{+0.4}_{-0.3}\times10^{34}$ erg and compare our detections with G star superflares detected in the \textit{Kepler} survey. We find our main flare to be one of the largest amplitude superflares detected from a bright G star. With energies more than 100 times greater than the Carrington event, our flare detections demonstrate the role that ground-based instruments such as NGTS can have in assessing the habitability of Earth-like exoplanets, particularly in the era of \textit{PLATO}.

Previous analyses of the spectra of OB-type stars in the Magellanic Clouds have identified targets with low projected rotational velocities and relatively high nitrogen abundances. The VFTS obtained spectroscopy for B-type 434 stars. We have estimated atmospheric parameters and nitrogen abundances using TLUSTY model atmospheres for 54 B-type targets that appear to be single, have projected rotational velocities <80 \kms and were not classified as supergiants. In addition, nitrogen abundances for 34 similar stars from a previous survey have been re-evaluated. Approximately 75-80% of the targets have nitrogen enhancements of less than 0.3 dex, consistent with them having experienced only small amounts of mixing. However, stars with low projected rotational velocities <40 \kms and significant nitrogen enrichments are found in both our samples and simulations imply that these cannot all be rapidly rotating objects observed near pole-on. Adopting an enhancement threshold of 0.6 dex, we observed five stars, yet stellar evolution models with rotation predict only 1.25$\pm$1.11 for random stellar viewing inclinations. The excess of such objects is estimated to be 20-30% of all stars with current rotational velocities of less than 40 \kms and 2-4% of the total single B-type sample. These estimates constitute lower limits for stars that appear inconsistent with current grids of stellar evolutionary models. Including targets with smaller nitrogen enhancements of greater than 0.2 dex implies larger percentages of targets that are inconsistent with current evolutionary models, viz. 70% of the stars with rotational velocities less than 40 \kms and 6-8% of the total single stellar population. We consider possible explanations of which the most promising would appear to be breaking due to magnetic fields or stellar mergers with subsequent magnetic braking.

One-opposition near-Earth asteroids (NEAs) are growing in number, and they must be recovered to prevent loss and mismatch risk, and to improve their orbits, as they are likely to be too faint for detection in shallow surveys at future apparitions. We aimed to recover more than half of the one-opposition NEAs recommended for observations by the Minor Planet Center (MPC) using the Isaac Newton Telescope (INT) in soft-override mode and some fractions of available D-nights. During about 130 hours in total between 2013 and 2016, we targeted 368 NEAs, among which 56 potentially hazardous asteroids (PHAs), observing 437 INT Wide Field Camera (WFC) fields and recovering 280 NEAs (76% of all targets). Engaging a core team of about ten students and amateurs, we used the THELI, Astrometrica, and the Find_Orb software to identify all moving objects using the blink and track-and-stack method for the faintest targets and plotting the positional uncertainty ellipse from NEODyS. Most targets and recovered objects had apparent magnitudes centered around V~22.8 mag, with some becoming as faint as V~24 mag. One hundred and three objects (representing 28% of all targets) were recovered by EURONEAR alone by Aug 2017. Orbital arcs were prolonged typically from a few weeks to a few years; our oldest recoveries reach 16 years. The O-C residuals for our 1,854 NEA astrometric positions show that most measurements cluster closely around the origin. In addition to the recovered NEAs, 22,000 positions of about 3,500 known minor planets and another 10,000 observations of about 1,500 unknown objects (mostly main-belt objects) were promptly reported to the MPC by our team. Four new NEAs were discovered serendipitously in the analyzed fields, increasing the counting to nine NEAs discovered by the EURONEAR in 2014 and 2015.

We present the discovery of NGTS-1b, a hot-Jupiter transiting an early M-dwarf host ($T_{eff}=3916^{+71}_{-63}~K$) in a P=2.674d orbit discovered as part of the Next Generation Transit Survey (NGTS). The planet has a mass of $0.812^{+0.066}_{-0.075}~M_{J}$, making it the most massive planet ever discovered transiting an M-dwarf. The radius of the planet is $1.33^{+0.61}_{-0.33}~R_{J}$. Since the transit is grazing, we determine this radius by modelling the data and placing a prior on the density from the population of known gas giant planets. NGTS-1b is the third transiting giant planet found around an M-dwarf, reinforcing the notion that close-in gas giants can form and migrate similar to the known population of hot Jupiters around solar type stars. The host star shows no signs of activity, and the kinematics hint at the star being from the thick disk population. With a deep (2.5%) transit around a $K=11.9$ host, NGTS-1b will be a strong candidate to probe giant planet composition around M-dwarfs via JWST transmission spectroscopy.

The Next Generation Transit Survey (NGTS), operating in Paranal since 2016, is a wide-field survey to detect Neptunes and super-Earths transiting bright stars, which are suitable for precise radial velocity follow-up and characterisation. Thereby, its sub-mmag photometric precision and ability to identify false positives are crucial. Particularly, variable background objects blended in the photometric aperture frequently mimic Neptune-sized transits and are costly in follow-up time. These objects can best be identified with the centroiding technique: if the photometric flux is lost off-centre during an eclipse, the flux centroid shifts towards the centre of the target star. Although this method has successfully been employed by the Kepler mission, it has previously not been implemented from the ground. We present a fully-automated centroid vetting algorithm developed for NGTS, enabled by our high-precision auto-guiding. Our method allows detecting centroid shifts with an average precision of 0.75 milli-pixel, and down to 0.25 milli-pixel for specific targets, for a pixel size of 4.97 arcsec. The algorithm is now part of the NGTS candidate vetting pipeline and automatically employed for all detected signals. Further, we develop a joint Bayesian fitting model for all photometric and centroid data, allowing to disentangle which object (target or background) is causing the signal, and what its astrophysical parameters are. We demonstrate our method on two NGTS objects of interest. These achievements make NGTS the first ground-based wide-field transit survey ever to successfully apply the centroiding technique for automated candidate vetting, enabling the production of a robust candidate list before follow-up.

We present a summary of the campaign of remote observations that supported the European Space Agency's Rosetta mission. Telescopes across the globe (and in space) followed comet 67P/Churyumov-Gerasimenko from before Rosetta's arrival until nearly the end of mission in September 2016. These provided essential data for mission planning, large-scale context information for the coma and tails beyond the spacecraft, and a way to directly compare 67P with other comets. The observations revealed 67P to be a relatively `well behaved' comet, typical of Jupiter family comets and with activity patterns that repeat from orbit-to-orbit. Comparison between this large collection of telescopic observations and the in situ results from Rosetta will allow us to better understand comet coma chemistry and structure. This work is just beginning as the mission ends -- in this paper we present a summary of the ground-based observations and early results, and point to many questions that will be addressed in future studies.

A spectroscopic analysis has been undertaken for the B-type multiple systems (excluding those with supergiant primaries) in the VLT-FLAMES Tarantula Survey (VFTS). Projected rotational velocities, $v$sin$i$, for the primaries have been estimated using a Fourier Transform technique and confirmed by fitting rotationally broadened profiles. A subset of 33 systems with $v$sin$i$$\leq$ 80 km/s have been analysed using a TLUSTY grid of model atmospheres to estimate stellar parameters and surface abundances for the primaries. The effects of a potential flux contribution from an unseen secondary have also been considered. For 20 targets it was possible to reliably estimate their effective temperatures (T$_{eff}$) but for the other 13 objects it was only possible to provide a constraint of 20,000$\leq$T$_{eff}$$\leq$26,000 K -- the other parameters estimated for these targets will be consequently less reliable. The estimated stellar properties are compared with evolutionary models and are generally consistent with their membership of 30 Doradus, while the nature of the secondaries of 3 SB2 system is discussed. A comparison with a sample of single stars with $v$sin$i$$\leq$ 80 km/s obtained from the VFTS and analysed with the same techniques implies that the atmospheric parameters and nitrogen abundances of the two samples are similar. However, the binary sample may have a lack of primaries with significant nitrogen enhancements, which would be consistent with them having low rotational velocities and having effectively evolved as single stars without significant rotational mixing. This result, which may be actually a consequence of the limitations of the pathfinder investigation presented in this paper, should be considered as a motivation for spectroscopic abundance analysis of large samples of binary stars, with high quality observational data.

A detailed knowledge of stellar activity is crucial for understanding stellar dynamos, as well as pushing exoplanet radial-velocity detection limits towards Earth analogue confirmation. We directly compare archival HARPS spectra taken at the minimum in $\alpha$ Cen B's activity cycle to a high-activity state when clear rotational modulation of $\log{R'_{HK}}$ is visible. Relative to the inactive spectra, we find a large number of narrow pseudo-emission features in the active spectra with strengths that are rotationally modulated. These features most likely originate from plage, spots, or a combination of both. They also display radial velocity variations of $\sim$300 m s$^{-1}$ - consistent with an active region rotating across the stellar surface. Furthermore, we see evidence that some of the lines originating from the `active immaculate' photosphere appear broader relative to the `inactive immaculate' case. This may be due to enhanced contributions of e.g. magnetic bright points to these lines, which then causes additional line broadening. More detailed analysis may enable measurements of plage and spot coverage using single spectra in the future.

We study the spatial correlations between the H$\alpha$ emission and different types of massive stars in two local galaxies, the Large Magellanic Cloud (LMC) and Messier 33. We compare these to correlations derived for core-collapse supernovae (CCSNe) in the literature to connect CCSNe of different types with the initial masses of their progenitors and to test the validity of progenitor mass estimates which use the pixel statistics method. We obtain samples of evolved massive stars in both galaxies from catalogues with good spatial coverage and/or completeness, and combine them with coordinates of main-sequence stars in the LMC from the SIMBAD database. We calculate the spatial correlation of stars of different classes and spectral types with H$\alpha$ emission. We also investigate the effects of distance, noise and positional errors on the pixel statistics method. A higher correlation with H$\alpha$ emission is found to correspond to a shorter stellar lifespan, and we conclude that the method can be used as an indicator of the ages, and therefore initial masses, of SN progenitors. We find that the spatial distributions of type II-P SNe and red supergiants of appropriate initial mass ($\gtrsim$9 $M_{\odot}$) are consistent with each other. We also find the distributions of type Ic SNe and WN stars with initial masses $\gtrsim$20 $M_{\odot}$ consistent, while supergiants with initial masses around 15 $M_{\odot}$ are a better match for type IIb and II-L SNe. The type Ib distribution corresponds to the same stellar types as type II-P, which suggests an origin in interacting binaries. On the other hand, we find that luminous blue variable stars show a much stronger correlation with H$\alpha$ emission than do type IIn SNe.

We present the results of large scale, three-dimensional magneto-hydrodynamics simulations of disc-winds for different initial magnetic field configurations. The jets are followed from the source to 90 AU scale, which covers several pixels of HST images of nearby protostellar jets. Our simulations show that jets are heated along their length by many shocks. We compute the emission lines that are produced, and find excellent agreement with observations. The jet width is found to be between 20 and 30 AU while the maximum velocities perpendicular to the jet is found to be up to above 100 km/s. The initially less open magnetic field configuration simulations results in a wider, two-component jet; a cylindrically shaped outer jet surrounding a narrow and much faster, inner jet. These simulations preserve the underlying Keplerian rotation profile of the inner jet to large distances from the source. However, for the initially most open magnetic field configuration the kink mode creates a narrow corkscrew-like jet without a clear Keplerian rotation profile and even regions where we observe rotation opposite to the disc (counter-rotating). The RW Aur jet is narrow, indicating that the disc field in that case is very open meaning the jet can contain a counter-rotating component that we suggests explains why observations of rotation in this jet has given confusing results. Thus magnetized disc winds from underlying Keplerian discs can develop rotation profiles far down the jet that are not Keplerian.

This white paper offers cautionary observations about the planning and development of new, large radio astronomy instruments. Complexity is a strong cost driver so every effort should be made to assign differing science requirements to different instruments and probably different sites. The appeal of shared resources is generally not realized in practice and can often be counterproductive. Instrument optimization is much more difficult with longer lists of requirements, and the development process is longer and less efficient. More complex instruments are necessarily further behind the technology state of the art because of longer development times. Including technology R&D in the construction phase of projects is a growing trend that leads to higher risks, cost overruns, schedule delays, and project de-scoping. There are no technology breakthroughs just over the horizon that will suddenly bring down the cost of collecting area. Advances come largely through careful attention to detail in the adoption of new technology provided by industry and the commercial market. Radio astronomy instrumentation has a very bright future, but a vigorous long-term R&D program not tied directly to specific projects needs to be restored, fostered, and preserved.

Gravitational lensing calculation using a direct inverse ray-shooting approach is a computationally expensive way to determine magnification maps, caustic patterns, and light-curves (e.g. as a function of source profile and size). However, as an easily parallelisable calculation, gravitational ray-shooting can be accelerated using programmable graphics processing units (GPUs). We present our implementation of inverse ray-shooting for the NVIDIA G80 generation of graphics processors using the NVIDIA Compute Unified Device Architecture (CUDA) software development kit. We also extend our code to multiple-GPU systems, including a 4-GPU NVIDIA S1070 Tesla unit. We achieve sustained processing performance of 182 Gflop/s on a single GPU, and 1.28 Tflop/s using the Tesla unit. We demonstrate that billion-lens microlensing simulations can be run on a single computer with a Tesla unit in timescales of order a day without the use of a hierarchical tree code.

The Atacama Large Millimeter/submillimeter Array (ALMA) is an international radio telescope under construction in the Atacama Desert of northern Chile. ALMA is situated on a dry site at 5000 m elevation, allowing excellent atmospheric transmission over the instrument wavelength range of 0.3 to 10 mm. ALMA will consist of two arrays of high-precision antennas. One, of up to 64 12-m diameter antennas, is reconfigurable in multiple patterns ranging in size from 150 meters up to ~15 km. A second array is comprised of a set of four 12-m and twelve 7-m antennas operating in one of two closely packed configurations ~50 m in diameter. The instrument will provide both interferometric and total-power astronomical information on atomic, molecular and ionized gas and dust in the solar system, our Galaxy, and the nearby to high-redshift universe. In this paper we outline the scientific drivers, technical challenges and planned progress of ALMA.