Primordial black holes formed through the collapse of cosmological density fluctuations have been hypothesised as contributors to the dark matter content of the Universe. At the same time, their mergers could contribute to the recently observed population of gravitational-wave sources. We investigate the scenario in which primordial black holes form binaries at late times in the Universe. Specifically, we re-examine the mergers of primordial black holes in small clusters of ~30 objects in the absence of initial binaries. Binaries form dynamically through Newtonian gravitational interactions. These binaries act as heat sources for the cluster, increasing the cluster's velocity dispersion, which inhibits direct mergers through gravitational-wave two-body captures. Meanwhile, three-body encounters of tight binaries are too rare to tighten binaries sufficiently to allow them to merge through gravitational-wave emission. We conclude that in the absence of initial binaries, merger rates of primordial black holes in the Bird et al. (2016) initial cluster configuration are at least an order of magnitude lower than previously suggested, which makes gravitational-wave detections of such sources improbable.

Supermassive black holes (SMBHs) are found in most galactic nuclei. A significant fraction of these nuclei also contain a nuclear stellar cluster (NSC) surrounding the SMBH. In this paper, we consider the idea that the NSC forms first, from the merger of several stellar clusters that may contain intermediate-mass black holes (IMBHs). These IMBHs can subsequently grow in the NSC and form an SMBH. We carry out $N$-body simulations of the simultaneous merger of three stellar clusters to form an NSC. We investigate the outcome of the simulated runs containing zero, one, two and three IMBHs. We find that IMBHs can efficiently sink to the centre of the merged cluster. If multiple merging clusters contain an IMBH, we find that an IMBH binary is likely to form and subsequently merge by gravitational wave emission. We show that these mergers are catalyzed by dynamical interactions with surrounding stars, which systematically harden the binary and increase its orbital eccentricity. The seed SMBH will be ejected from the NSC by the recoil kick produced when two IMBHs merge, if their mass ratio $q\gtrsim 0.15$. We find that the other channel, namely, dynamical interactions between an IMBH binary and a third single IMBH, is rarely effective in removing the seed. If the seed is ejected then no SMBH will form in the NSC. This is a natural pathway to explain those galactic nuclei that contain an NSC but apparently lack an SMBH, such as M33. However, if an IMBH is retained then it can seed the growth of an SMBH through gas accretion and tidal disruption of stars.

ESA recently called for new "Science Ideas" to be investigated in terms of feasibility and technological developments -- for technologies not yet sufficiently mature. These ideas may in the future become candidates for M or L class missions within the ESA Science Program. With the launch of Gaia in December 2013, Europe entered a new era of space astrometry following in the footsteps of the very successful Hipparcos mission from the early 1990s. Gaia is the successor to Hipparcos, both of which operated in optical wavelengths, and Gaia is two orders of magnitude more accurate in the five astrometric parameters and is surveying four orders of magnitude more stars in a vast volume of the Milky Way. The combination of the Hipparcos/Tycho-2 catalogues with the first early Gaia data release will give improved proper motions over a long ~25 year baseline. The final Gaia solution will also establish a new optical reference frame by means of quasars, by linking the optical counterparts of radio (VLBI) sources defining the orientation of the reference frame, and by using the zero proper motion of quasars to determine a non-rotating frame. A weakness of Gaia is that it only operates at optical wavelengths. However, much of the Galactic centre and the spiral arm regions, important for certain studies, are obscured by interstellar extinction and this makes it difficult for Gaia to deeply probe. Traditionally, this problem is overcome by switching to the infra-red but this was not possible with Gaia's CCDs. Additionally, to scan the entire sky and make global absolute parallax measurements the spacecraft must have a constant rotation and this requires that the CCDs operate in TDI mode, increasing their complexity.

Binary stellar systems form a large fraction of the Galaxy's stars. They are useful as laboratories for studying the physical processes taking place within stars, and must be correctly taken into account when observations of stars are used to study the structure and evolution of the Galaxy. We present a sample of 12760 well-characterised double-lined spectroscopic binaries that are appropriate for statistical studies of the binary populations. They were detected as SB2s using a t-distributed stochastic neighbour embedding (t-SNE) classification and a cross-correlation analysis of GALAH spectra. This sample consists mostly of dwarfs, with a significant fraction of evolved stars and several dozen members of the giant branch. To compute parameters of the primary and secondary star ($T_{\rm eff[1,2]}$, $\log g_{[1,2]}$, [Fe/H], $V_{r[1,2]}$, $v_{\rm mic[1,2]}$, $v_{\rm broad[1,2]}$, $R_{[1,2]}$, and $E(B-V)$), we used a Bayesian approach that includes a parallax prior from Gaia DR2, spectra from GALAH, and apparent magnitudes from APASS, Gaia DR2, 2MASS, and WISE. The derived stellar properties and their distributions show trends that are expected for a population of close binaries (a $<$ 10 AU) with mass ratios $0.5 \leq q \leq 1$. The derived metallicity of these binary stars is statistically lower than that of single dwarf stars from the same magnitude-limited sample.

We investigate the hypothesis that ingestion of a terrestrial or super-Earth planet could cause the anomalously high metal abundances seen in a turn-off star in the open cluster M67, when compared to other turn-off stars in the same cluster. We show that the mass in convective envelope of the star is likely only $3.45\,\times 10^{-3}\,{\rm M}_\odot$, and hence $5.2\,{\rm M}_\oplus$ of rock is required to obtain the observed 0.128 dex metal enhancement. Rocky planets dissolve entirely in the convective envelope if they enter it with sufficiently tangential orbits: we find that the critical condition for dissolution is that the planet's radial speed must be less than 40% of its total velocity at the stellar surface; or, equivalently, the impact parameter must be greater than about 0.9. We model the delivery of rocky planets to the stellar surface both by planet-planet scattering in a realistic multi-planet system, and by Lidov-Kozai cycles driven by a more massive planetary or stellar companion. In both cases almost all planets that are ingested arrive at the star on grazing orbits and hence will dissolve in the surface convection zone. We conclude that super-Earth ingestion is a good explanation for the metal enhancement in M67 Y2235, and that a high-resolution spectroscopic survey of stellar abundances around the turn-off and main sequence of M67 has the potential to constrain the frequency of late-time dynamical instability in planetary systems.

Supermassive black holes are found in most galactic nuclei. A large fraction of these nuclei also contain a nuclear stellar cluster surrounding the black hole. Here we consider the idea that the nuclear stellar cluster formed first and that the supermassive black hole grew later. In particular we consider the merger of three stellar clusters to form a nuclear stellar cluster, where some of these clusters contain a single intermediate-mass black hole (IMBH). In the cases where multiple clusters contain IMBHs, we discuss whether the black holes are likely to merge and whether such mergers are likely to result in the ejection of the merged black hole from the nuclear stellar cluster. In some cases, no supermassive black hole will form as any merger product is not retained. This is a natural pathway to explain those galactic nuclei that contain a nuclear stellar cluster but apparently lack a supermassive black hole; M33 being a nearby example. Alternatively, if an IMBH merger product is retained within the nuclear stellar cluster, it may subsequently grow, e.g. via the tidal disruption of stars, to form a supermassive black hole.

We investigate the possibility that the Sun could have been born in M67 by carrying out $N$-body simulations of M67-like clusters in a time-varying Galactic environment, and following the galactic orbits of stars that escape from them. We find that model clusters that occupy similar orbits to M67 today can be divided up into three groups. Hot clusters are born with a high initial $z$-velocity, depleted clusters are born on cold orbits but are destroyed by GMC encounters in the Galactic disc, and scattered clusters are born on cold orbits and survive with more than 1000 stars at an age of 4.6 Gyr. We find that all cluster models in all three cluster groups have stellar escapers that are kinematicaly similar to the Sun. Hot clusters having the lowest such fraction $f_{\odot} = 0.06$ %, whilst depleted clusters have the highest fraction, $f_{\odot} = 6.61$ %. We calculate that clusters that are destroyed in the Galactic disc have a specific frequency of escapers that end up on solar-like orbits that is $\sim$ 2 times that of escapers from clusters that survive their journey.

The 4MOST instrument is a multi-object spectrograph that will address Galactic and extragalactic science cases simultaneously by observing targets from a large number of different surveys within each science exposure. This parallel mode of operation and the survey nature of 4MOST require some distinct 4MOST-specific operational features within the overall operations model of ESO. The main feature is that the 4MOST Consortium will deliver, not only the instrument, but also contractual services to the user community, which is why 4MOST is also described as a facility. This white paper concentrates on information particularly useful to answering the forthcoming Call for Letters of Intent.

The signatures of the formation and evolution of a galaxy are imprinted in its stars. Their velocities, ages, and chemical compositions present major constraints on models of galaxy formation, and on various processes such as the gas inflows and outflows, the accretion of cold gas, radial migration, and the variability of star formation activity. Understanding the evolution of the Milky Way requires large observational datasets of stars via which these quantities can be determined accurately. This is the science driver of the 4MOST MIlky way Disc And BuLgE High-Resolution (4MIDABLE-HR) survey: to obtain high-resolution spectra at $R \sim 20\,000$ and to provide detailed elemental abundances for large samples of stars in the Galactic disc and bulge. High data quality will allow us to provide accurate spectroscopic diagnostics of two million stellar spectra: precise radial velocities; rotation; abundances of many elements, including those that are currently only accessible in the optical, such as Li, s-, and r-process; and multi-epoch spectra for a sub-sample of stars. Synergies with complementary missions like Gaia and TESS will provide masses, stellar ages and multiplicity, forming a multi-dimensional dataset that will allow us to explore and constrain the origin and structure of the Milky Way.

The mechanisms of the formation and evolution of the Milky Way are encoded in the orbits, chemistry and ages of its stars. With the 4MOST MIlky way Disk And BuLgE Low-Resolution Survey (4MIDABLE-LR) we aim to study kinematic and chemical substructures in the Milky Way disc and bulge region with samples of unprecedented size out to larger distances and greater precision than conceivable with Gaia alone or any other ongoing or planned survey. Gaia gives us the unique opportunity for target selection based almost entirely on parallax and magnitude range, hence increasing the efficiency in sampling larger Milky Way volumes with well-defined and effective selection functions. Our main goal is to provide a detailed chrono-chemo-kinematical extended map of our Galaxy and the largest Gaia follow-up down to $G = 19$ magnitudes (Vega). The complex nature of the disc components (for example, large target densities and highly structured extinction distribution in the Milky Way bulge and disc area), prompted us to develop a survey strategy with five main sub-surveys that are tailored to answer the still open questions about the assembly and evolution of our Galaxy, while taking full advantage of the Gaia data.

We will study the formation history of the Milky Way, and the earliest phases of its chemical enrichment, with a sample of more than 1.5 million stars at high galactic latitude. Elemental abundances of up to 20 elements with a precision of better than 0.2 dex will be derived for these stars. The sample will include members of kinematically coherent substructures, which we will associate with their possible birthplaces by means of their abundance signatures and kinematics, allowing us to test models of galaxy formation. Our target catalogue is also expected to contain 30,000 stars at a metallicity of less than one hundredth that of the Sun. This sample will therefore be almost a factor of 100 larger than currently existing samples of metal-poor stars for which precise elemental abundances are available (determined from high-resolution spectroscopy), enabling us to study the early chemical evolution of the Milky Way in unprecedented detail.

The current status of and motivation for the 4MOST Survey Strategy, as developed by the Consortium science team, are presented here. Key elements of the strategy are described, such as sky coverage, number of visits and total exposure times in different parts of the sky, and how to deal with different observing conditions. The task of organising the strategy is not simple, with many different surveys that have vastly different target brightnesses and densities, sample completeness levels, and signal-to-noise requirements. We introduce here a number of concepts that we will use to ensure all surveys are optimised. Astronomers who are planning to submit a Participating Survey proposal are strongly encouraged to read this article and any relevant 4MOST Survey articles in this issue of The Messenger such that they can optimally complement and benefit from the planned surveys of the 4MOST Consortium.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

We analyse N-body simulations of star-forming regions to investigate the effects of external far and extreme ultra-violet photoevaporation from massive stars on protoplanetary discs. By varying the initial conditions of simulated star-forming regions, such as the spatial distribution, net bulk motion (virial ratio), and density, we investigate which parameters most affect the rate at which discs are dispersed due to external photoevaporation. We find that disc dispersal due to external photoevaporation is faster in highly substructured star-forming regions than in smooth and centrally concentrated regions. Sub-virial star-forming regions undergoing collapse also show higher rates of disc dispersal than regions that are in virial equilibrium or are expanding. In moderately dense ($\sim$100 M$_{\odot}$ pc$^{-3}$) regions, half of all protoplanetary discs with radii $\geq$ 100 AU are photoevaporated within 1 Myr, three times faster than is currently suggested by observational studies. Discs in lower-density star-forming regions ($\sim$10 M$_{\odot}$ pc$^{-3}$) survive for longer, but half are still dispersed on short timescales ($\sim$2 Myr). This demonstrates that the initial conditions of the star forming regions will greatly impact the evolution and lifetime of protoplanetary discs. These results also imply that either gas giant planet formation is extremely rapid and occurs before the gas component of discs is evaporated, or gas giants only form in low-density star-forming regions where no massive stars are present to photoevaporate gas from protoplanetary discs.

A deep understanding of the Milky Way galaxy, its formation and evolution requires observations of huge numbers of stars. Stellar photometry, therefore, provides an economical method to obtain intrinsic stellar parameters. With the addition of distance information - deriving reliable ages from photometry is a possibility. We have developed a Bayesian method that generates 2D probability maps of a star's age and metallicity from photometry and parallax using isochrones. Our synthetic tests show that including a near-UV passband enables us to break the degeneracy between a star's age and metallicity for certain evolutionary stages. It is possible to find well-constrained ages and metallicities for turn-off and sub-giant stars with colours including a U band and a parallax with uncertainty less than ~20%. Metallicities alone are possible for the main sequence and giant branch. We find good agreement with the literature when we apply our method to the Gaia benchmark stars, particularly for turn-off and young stars. Further tests on the old open cluster NGC 188, however, reveal significant limitations in the stellar isochrones. The ages derived for the cluster stars vary with evolutionary stage, such that turn-off ages disagree with those on the sub-giant branch, and metallicities vary significantly throughout. Furthermore, the parameters vary appreciably depending on which colour combinations are used in the derivation. We identify the causes of these mismatches and show that improvements are needed in the modelling of giant branch stars and in the creation and calibration of synthetic near-UV photometry. Our results warn against applying isochrone fitting indiscriminately. In particular, the uncertainty on the stellar models should be quantitatively taken into account. Further efforts to improve the models will result in significant advancements in our ability to study the Galaxy.

In the era of large-scale surveys of stars in the Milky Way, stellar ages are crucial for studying the evolution of the Galaxy. But determining ages of field stars is notoriously difficult; therefore, we attempt to determine benchmark ages for the extensively studied Gaia benchmark stars which can be used for validation purposes. By searching the literature for age estimates from different methods and deriving new ages based on Bayesian isochrone fitting, we are able to put reliable limits on the ages of 16 out of the 33 benchmark stars. The giants with well-defined ages are all young, and an expansion of the sample to include older giants with asteroseismic ages would be beneficial. Some of the stars have surface parameters inconsistent with isochrones younger than 16 Gyr. Including $\alpha$-enhancement in the models when relevant resolves some of these cases, but others clearly highlight discrepancies between the models and observations. We test the impact of atomic diffusion on the age estimates by fitting to the actual surface metallicity of the models instead of the initial value and find that the effect is negligible except for a single turn-off star. Finally, we show that our ability to determine isochrone-based ages for large spectroscopic surveys largely mirrors our ability to determine ages for these benchmark stars, except for stars with $\log g \gtrsim 4.4$ dex since their location in the HR diagram is almost age insensitive. Hence, isochrone fitting does not constrain their ages given the typical uncertainties of spectroscopic stellar parameters.

Extremely metal-poor (EMP) stars are an integral piece in the puzzle that is the early Universe, and although anomolous subclasses of EMP stars such as carbon-enhanced metal-poor (CEMP) stars are well-studied, they make up less than half of all EMP stars with [Fe/H] $\sim -3.0$. The amount of carbon depletion occurring on the red giant branch (carbon offset) is used to determine the evolutionary status of EMP stars, and this offset will differ between CEMP and normal EMP stars. The depletion mechanism employed in stellar models (from which carbon offfsets are derived) is very important, however the only widely available carbon offsets in the literature are derived from stellar models using a thermohaline mixing mechanism that cannot simultaneously match carbon and lithium abundances to observations for a single diffusion coeffcient. Our stellar evolution models utilise a modified thermohaline mixing model that can match carbon and lithium in the metal-poor globular cluster NGC 6397. We compare our models to the bulk of the EMP star sample at [Fe/H] $= -3$ and show that our modified models follow the trend of the observations and deplete less carbon compared to the standard thermohaline mixing theory. We conclude that stellar models that employ the standard thermohaline mixing formalism overestimate carbon offsets and hence CEMP star frequencies, particularly at metallicities where carbon-normal stars dominate the EMP star population.

The Gaia DR2 proper motion places HV 2112 in the Small Magellanic Cloud (SMC), consistent with it being a Thorne-\.Zytkow object.

Many field stars reside in binaries, and the analysis and interpretation of photometric and spectroscopic surveys must take this into account. We have developed a model to predict how binaries influence the scientific results inferred from large spectroscopic surveys. Based on the rapid binary evolution code BSE, it allows us to model a representative population of binaries and generate synthetic survey observations. We describe this model in detail, and apply it to the radial velocity variation of subgiant and giant stars in the Galactic disc, as observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE), part of the Sloan Digital Sky Survey III. APOGEE provides an excellent data set for testing our binary models since a large fraction of the stars have been observed multiple times. By comparing our model to the APOGEE observations we constrain the initial binary fraction of solar-metallicity stars in the sample to be $f_{\rm b,0}=0.35\pm0.01$, in line with the solar neighbourhood. We find that the binary fraction is higher at lower metallicities, consistent with other observational studies. Our model matches the shape of the high-velocity scatter in APOGEE, which suggests that most velocity variability above 0.5 km/s comes from binaries. Our exploration of binary initial properties shows that APOGEE is mostly sensitive to binaries with periods between 3 and 3000 years, and is largely insensitive to the detailed properties of the population. We can, however, rule out a population where the mass of the lower-mass star is drawn from the IMF independently of its companion.

The luminous X-ray binary 47 Tuc X9 shows radio and X-ray emission consistent with a stellar-mass black hole accreting from a carbon-oxygen white dwarf. Its location, in the core of the massive globular cluster 47 Tuc, hints at a dynamical origin. We assess the stability of mass transfer from a carbon-oxygen white dwarf onto compact objects of various masses, and conclude that for mass transfer to proceed stably the accretor must, in fact, be a black hole. Such systems can form dynamically by the collision of a stellar-mass black hole with a giant star. Tidal dissipation of energy in the giant's envelope leads to a bound binary with a pericentre separation less than the radius of the giant. An episode of common-envelope evolution follows, which ejects the giant's envelope. We find that the most likely target is a horizontal-branch star, and that a realistic quantity of subsequent dynamical hardening is required for the resulting binary to merge via gravitational wave emission. Observing one binary like 47 Tuc X9 in the Milky Way globular cluster system is consistent with the expected formation rate. The observed 6.8-day periodicity in the X-ray emission may be driven by eccentricity induced in the UCXB's orbit by a perturbing companion.

We present our models of the effect of binaries on high-resolution spectroscopic surveys. We want to determine how many binary stars will be observed, whether unresolved binaries will contaminate measurements of chemical abundances, and how we can use spectroscopic surveys to better constrain the population of binary stars in the Galaxy. Using a rapid binary-evolution algorithm that enables modelling of the most complex binary systems we generate a series of large binary populations in the Galactic disc and evaluate the results. As a first application we use our model to study the binary fraction in APOGEE giants. We find tentative evidence for a change in binary fraction with metallicity.

We perform hydrodynamic simulations of mass transfer in binaries that contain a white dwarf and a neutron star (WD-NS binaries), and measure the specific angular momentum of material lost from the binary in disc winds. By incorporating our results within a long-term evolution model we measure the long-term stability of mass transfer in these binaries. We find that only binaries containing helium white dwarfs with masses less than a critical mass of $M_{\rm WD,crit}=0.2\,M_\odot$ undergo stable mass transfer and evolve into ultra-compact $X$-ray binaries. Systems with higher-mass white dwarfs experience unstable mass transfer, which leads to tidal disruption of the white dwarf. Our low critical mass compared to the standard jet-only model of mass loss arises from the efficient removal of angular momentum in the mechanical disc winds which develop at highly super-Eddington mass-transfer rates. We find that the eccentricities expected for WD-NS binaries when they come into contact do not affect the loss of angular momentum, and can only affect the long-term evolution if they change on shorter timescales than the mass-transfer rate. Our results are broadly consistent with the observed numbers of both ultra-compact $X$-ray binaries and radio pulsars with white dwarf companions. The observed calcium-rich gap transients are consistent with the merger rate of unstable systems with higher-mass white dwarfs.

HR 8799 is a star accompanied by four massive planets on wide orbits. The observed planetary configuration has been shown to be unstable on a timescale much shorter than the estimated age of the system (~ 30 Myr) unless the planets are locked into mean motion resonances. This condition is characterised by small-amplitude libration of one or more resonant angles that stabilise the system by preventing close encounters. We simulate planetary systems similar to the HR 8799 planetary system, exploring the parameter space in separation between the orbits, planetary masses and distance from the Sun to the star. We find systems that look like HR 8799 and remain stable for longer than the estimated age of HR 8799. None of our systems are forced into resonances. We find, with nominal masses and in a narrow range of orbit separations, that 5 of 100 systems match the observations and lifetime. Considering a broad range of orbit separations, we find 12 of 900 similar systems. The systems survive significantly longer because of their slightly increased initial orbit separations compared to assuming circular orbits from the observed positions. A small increase in separation leads to a significant increase in survival time. The low eccentricity the orbits develop from gravitational interaction is enough for the planets to match the observations. With lower masses, but still comfortably within the estimated planet mass uncertainty, we find 18 of 100 matching and long-lived systems in a narrow orbital separation range. In the broad separation range, we find 82 of 900 matching systems. Our results imply that the planets in the HR 8799 system do not have to be in strong mean motion resonances.

Could the velocity spread, increasing with time, in the Galactic disk be explained as a result of gravitational interactions of stars with giant molecular clouds (GMCs) and spiral arms? Do the old open clusters high above the Galactic plane provide clues to this question? We explore the effects on stellar orbits of scattering by inhomogeneities in the Galactic potential due to GMCs, spiral arms and the Galactic bar, and whether high-altitude clusters could have formed in orbits closer to the Galactic plane and later been scattered. Simulations of test-particle motions are performed in a realistic Galactic potential. The effects of the internal structure of GMCs are explored. The destruction of clusters in GMC collisions is treated in detail with N-body simulations of the clusters. The observed velocity dispersions of stars as a function of time are well reproduced. The GMC structure is found to be significant, but adequate models produce considerable scattering effects. The fraction of simulated massive old open clusters, scattered into orbits with |z| > 400 pc, is typically 0:5%, in agreement with the observed number of high-altitude clusters and consistent with the present formation rate of massive open clusters. The heating of the thin Galactic disk is well explained by gravitational scattering by GMCs and spiral arms, if the local correlation between the GMC mass and the corresponding voids in the gas is not very strong. Our results suggest that the high-altitude metal-rich clusters were formed in orbits close to the Galactic plane and later scattered to higher orbits. It is possible, though not very probable, that the Sun formed in such a cluster before scattering occurred.

Calcium-rich supernovae represent a significant challenge for our understanding of the fates of stellar systems. They are less luminous than other supernova (SN) types and they evolve more rapidly to reveal nebular spectra dominated by strong calcium lines with weak or absent signatures of other intermediate- and iron-group elements, which are seen in other SNe. Strikingly, their explosion sites also mark them out as distinct from other SN types. Their galactocentric offset distribution is strongly skewed to very large offsets (around one third are offset greater than 20 kpc), meaning they do not trace the stellar light of their hosts. Many of the suggestions to explain this extreme offset distribution have invoked the necessity for unusual formation sites such as globular clusters or dwarf satellite galaxies, which are therefore difficult to detect. Building on previous work attempting to detect host systems of nearby Ca-rich SNe, we here present Hubble Space Telescope imaging of 5 members of the class - 3 exhibiting large offsets and 2 coincident with the disk of their hosts. We find no underlying sources at the explosion sites of any of our sample. Combining with previous work, the lack of a host system now appears to be a ubiquitous feature amongst Ca-rich SNe. In this case the offset distribution is most readily explained as a signature of high-velocity progenitor systems that have travelled significant distances before exploding.

We conduct a series of comparisons between spectroscopic and photometric observations of globular clusters and stellar models to examine their predictive power. Data from medium-to-high resolution spectroscopic surveys of lithium allow us to investigate first dredge-up and extra mixing in two clusters well separated in metallicity. Abundances at first dredge-up are satisfactorily reproduced but there is preliminary evidence to suggest that the models overestimate the luminosity at which the surface composition first changes in the lowest-metallicity system. Our models also begin extra mixing at luminosities that are too high, demonstrating a significant discrepancy with observations at low metallicity. We model the abundance changes during extra mixing as a thermohaline process and determine that the usual diffusive form of this mechanism cannot simultaneously reproduce both the carbon and lithium observations. Hubble Space Telescope photometry provides turnoff and bump magnitudes in a large number of globular clusters and offers the opportunity to better test stellar modelling as function of metallicity. We directly compare the predicted main-sequence turn-off and bump magnitudes as well as the distance-independent parameter $\Delta M_V ~^{\rm{MSTO}}_{\rm{bump}}$. We require 15 Gyr isochrones to match the main-sequence turn-off magnitude in some clusters and cannot match the bump in low-metallicity systems. Changes to the distance modulus, metallicity scale and bolometric corrections may impact on the direct comparisons but $\Delta M_V ~^{\rm{MSTO}}_{\rm{bump}}$, which is also underestimated from the models, can only be improved through changes to the input physics. Overshooting at the base of the convective envelope with an efficiency that is metallicity dependent is required to reproduce the empirically determined value of $\Delta M_V ~^{\rm{MSTO}}_{\rm{bump}}$.

Instabilities in planetary systems can result in the ejection of planets from their host system, resulting in free-floating planets (FFPs). If this occurs in a star cluster, the FFP may remain bound to the star cluster for some time and interact with the other cluster members until it is ejected. Here, we use $N$-body simulations to characterise close star-planet and planet-planet encounters and the dynamical fate of the FFP population in star clusters containing $500-2000$ single or binary star members. We find that FFPs ejected from their planetary system at low velocities typically leave the star cluster 40% earlier than their host stars, and experience tens of close ($<1000$ AU) encounters with other stars and planets before they escape. The fraction of FFPs that experiences a close encounter depends on both the stellar density and the initial velocity distribution of the FFPs. Approximately half of the close encounters occur within the first 30 Myr, and only 10% occur after 100 Myr. The periastron velocity distribution for all encounters is well-described by a modified Maxwell-Bolzmann distribution, and the periastron distance distribution is linear over almost the entire range of distances considered, and flattens off for very close encounters due to strong gravitational focusing. Close encounters with FFPs can perturb existing planetary systems and their debris structures, and they can result in re-capture of FFPs. In addition, these FFP populations may be observed in young star clusters in imaging surveys; a comparison between observations and dynamical predictions may provide clues to the early phases of stellar and planetary dynamics in star clusters.

In recent years much interest has been shown in the process of thermohaline mixing in red giants. In low and intermediate mass stars this mechanism first activates at the position of the bump in the luminosity function, and has been identified as a likely candidate for driving the slow mixing inferred to occur in these stars. One particularly important consequence of this process, which is driven by a molecular weight inversion, is the destruction of lithium. We show that the degree of lithium destruction, or in some cases production, is extremely sensitive to the numerical details of the stellar models. Within the standard 1D diffusion approximation to thermohaline mixing, we find that dfferent evolution codes, with their default numerical schemes, can produce lithium abundances that differ from one another by many orders of magnitude. This disagreement is worse for faster mixing. We perform experiments with four independent stellar evolution codes, and derive conditions for the spatial and temporal resolution required for a converged numerical solution. The results are extremely sensitive to the timesteps used. We find that predicted lithium abundances published in the literature until now should be treated with caution.

We present deep VLT and HST observations of the nearest examples of calcium-rich 'gap' transients -- rapidly evolving transient events, with a luminosity intermediate between novae and supernovae. These sources are frequently found at large galactocentric offsets, and their progenitors remain mysterious. Our observations find no convincing underlying quiescent sources coincident with the locations of these transients, allowing us to rule out a number of potential progenitor systems. The presence of surviving massive-star binary companions (or other cluster members) are ruled out, providing an independent rejection of a massive star origin for these events. Dwarf satellite galaxies are disfavoured unless one invokes as yet unknown conditions that would be extremely favourable for their production in the lowest mass systems. Our limits also probe the majority of the globular cluster luminosity function, ruling out the presence of an underlying globular cluster population at high significance, and thus the possibility that they are created via dynamical interactions in dense globular cluster cores. Given the lack of underlying systems, previous progenitor suggestions have difficulty reproducing the remote locations of these transients, even when considering solely halo-borne progenitors. Our preferred scenario is that calcium-rich transients are high velocity, kicked systems, exploding at large distances from their natal site. Coupled with a long-lived progenitor system post-kick, this naturally explains the lack of association these transients have with their host stellar light, and the extreme host-offsets exhibited. Neutron star -- white dwarf mergers may be a promising progenitor system in this scenario.

The very bright red star HV2112 in the Small Magellanic Cloud could be a massive Thorne-Zytkow Object, a supergiant-like star with a degenerate neutron core. With its luminosity of over $10^5\,\rm L_\odot$, it could also be a super asymptotic giant branch star, a star with an oxygen/neon core supported by electron degeneracy and undergoing thermal pulses with third dredge up. Both TZOs and SAGB stars are expected to be rare. Abundances of heavy elements in HV2112's atmosphere, as observed to date, do not allow us to distinguish between the two possibilities based on the latest models. Molybdenum and rubidium can be enhanced by both the irp-process in a TZO or by the s-process in SAGB stars. Lithium can be generated by hot bottom burning at the base of the convective envelope in either. HV2112's enhanced calcium could thus be the key determinant. A SAGB star is not able to synthesise its own calcium but it may be possible to produce this in the final stages of the process that forms a TZO, when the degenerate electron core of a giant star is tidally disrupted by a neutron star. Hence our calculations indicate that HV2112 is most likely a genuine TZO.

PLATO 2.0 has recently been selected for ESA's M3 launch opportunity (2022/24). Providing accurate key planet parameters (radius, mass, density and age) in statistical numbers, it addresses fundamental questions such as: How do planetary systems form and evolve? Are there other systems with planets like ours, including potentially habitable planets? The PLATO 2.0 instrument consists of 34 small aperture telescopes (32 with 25 sec readout cadence and 2 with 2.5 sec candence) providing a wide field-of-view (2232 deg2) and a large photometric magnitude range (4-16 mag). It focusses on bright (4-11 mag) stars in wide fields to detect and characterize planets down to Earth-size by photometric transits, whose masses can then be determined by ground-based radial-velocity follow-up measurements. Asteroseismology will be performed for these bright stars to obtain highly accurate stellar parameters, including masses and ages. The combination of bright targets and asteroseismology results in high accuracy for the bulk planet parameters: 2%, 4-10% and 10% for planet radii, masses and ages, respectively. The planned baseline observing strategy includes two long pointings (2-3 years) to detect and bulk characterize planets reaching into the habitable zone (HZ) of solar-like stars and an additional step-and-stare phase to cover in total about 50% of the sky. PLATO 2.0 will observe up to 1,000,000 stars and detect and characterize hundreds of small planets, and thousands of planets in the Neptune to gas giant regime out to the HZ. It will therefore provide the first large-scale catalogue of bulk characterized planets with accurate radii, masses, mean densities and ages. This catalogue will include terrestrial planets at intermediate orbital distances, where surface temperatures are moderate. Coverage of this parameter range with statistical numbers of bulk characterized planets is unique to PLATO 2.0.

We use N-body simulations of star cluster evolution to explore the hypothesis that short-lived radioactive isotopes found in meteorites, such as 26-Al, were delivered to the Sun's protoplanetary disc from a supernova at the epoch of Solar System formation. We cover a range of star cluster formation parameter space and model both clusters with primordial substructure, and those with smooth profiles. We also adopt different initial virial ratios - from cool, collapsing clusters to warm, expanding associations. In each cluster we place the same stellar population; the clusters each have 2100 stars, and contain one massive 25M_Sun star which is expected to explode as a supernova at about 6.6Myr. We determine the number of Solar (G)-type stars that are within 0.1 - 0.3pc of the 25M_Sun star at the time of the supernova, which is the distance required to enrich the protoplanetary disc with the 26-Al abundances found in meteorites. We then determine how many of these G-dwarfs are unperturbed `singletons'; stars which are never in close binaries, nor suffer sub-100au encounters, and which also do not suffer strong dynamical perturbations. The evolution of a suite of twenty initially identical clusters is highly stochastic, with the supernova enriching over 10 G-dwarfs in some clusters, and none at all in others. Typically only ~25 per cent of clusters contain enriched, unperturbed singletons, and usually only 1 - 2 per cluster (from a total of 96 G-dwarfs in each cluster). The initial conditions for star formation do not strongly affect the results, although a higher fraction of supervirial (expanding) clusters would contain enriched G-dwarfs if the supernova occurred earlier than 6.6Myr. If we sum together simulations with identical initial conditions, then ~1 per cent of all G-dwarfs in our simulations are enriched, unperturbed singletons.

The planet candidates discovered by the Kepler mission provide a rich sample to constrain the architectures and relative inclinations of planetary systems within approximately 0.5 AU of their host stars. We use the triple-transit systems from the Kepler 16-months data as templates for physical triple-planet systems and perform synthetic transit observations. We find that all the Kepler triple-transit and double-transit systems can be produced from the triple-planet templates, given a low mutual inclination of around five degrees. Our analysis shows that the Kepler data contains a population of planets larger than four Earth radii in single-transit systems that can not arise from the triple-planet templates. We explore the hypothesis that high-mass counterparts of the triple-transit systems underwent dynamical instability to produce a population of massive double-planet systems of moderately high mutual inclination. We perform N-body simulations of mass-boosted triple-planet systems and observe how the systems heat up and lose planets, most frequently by planet-planet collisions, yielding transits in agreement with the large planets in the Kepler single-transit systems. The resulting population of massive double-planet systems can nevertheless not explain the additional excess of low-mass planets among the observed single-transit systems and the lack of gas-giant planets in double-transit and triple-transit systems. Planetary instability of systems of triple gas-giant planets can be behind part of the dichotomy between systems hosting one or more small planets and those hosting a single giant planet. The main part of the dichotomy, however, is more likely to have arisen already during planet formation when the formation, migration or scattering of a massive planet, triggered above a threshold metallicity, suppressed the formation of other planets in sub-AU orbits.

We consider a popular model for long-duration gamma-ray bursts, in which the progenitor star, a stripped helium core, is spun up by tidal interactions with a black- hole companion in a compact binary. We perform population synthesis calculations to produce a representative sample of such binaries, and model the effect that the companion has on material that falls back on to the newly-formed black hole. Taking the results of hydrodynamic models of black-hole formation by fallback as our starting point, we show that the companion has two main effects on the fallback process. First, a break forms in the accretion curve at around 10 000 s. Secondly, subsequent to the break, we expect to see a flare of total energy around 0.1 foe. We predict that the break time is set largely by the semi-major axis of the binary at the time of explosion, and that this correlates negatively with the flare energy. Although comparison with observations is non-trivial, we show that our predicted break times are comparable to those found in the X-ray light curves of canonical long-duration gamma-ray bursts. Similarly, the flare properties that we predict are consistent with the late-time flares observed in a sub-sample of bursts.

It is now widely accepted that globular cluster red giant branch stars owe their strange abundance patterns to a combination of pollution from progenitor stars and in situ extra mixing. In this hybrid theory a first generation of stars imprint abundance patterns into the gas from which a second generation forms. The hybrid theory suggests that extra mixing is operating in both populations and we use the variation of [C/Fe] with luminosity to examine how efficient this mixing is. We investigate the observed red giant branches of M3, M13, M92, M15 and NGC 5466 as a means to test a theory of thermohaline mixing. The second parameter pair M3 and M13 are of intermediate metallicity and our models are able to account for the evolution of carbon along the RGB in both clusters. Although, in order to fit the most carbon-depleted main-sequence stars in M13 we require a model whose initial [C/Fe] abundance leads to a carbon abundance lower than is observed. Furthermore our results suggest that stars in M13 formed with some primary nitrogen (higher C+N+O than stars in M3). In the metal-poor regime only NGC 5466 can be tentatively explained by thermohaline mixing operating in multiple populations. We find thermohaline mixing unable to model the depletion of [C/Fe] with magnitude in M92 and M15. It appears as if extra mixing is occurring before the luminosity function bump in these clusters. To reconcile the data with the models would require first dredge-up to be deeper than found in extant models.

We compare the observed spatial offsets of short gamma-ray bursts from their host galaxies with their predicted distributions, assuming that they originate in double neutron star binaries that form from field stars. We find that, for the majority of bursts, this model is sufficient to explain the observed offsets, although there is a trend towards larger offsets than predicted. One burst, GRB 060502B, has an offset that is clearly anomalous. We discuss possible reasons for the large offsets, including host galaxy misidentification, and suggest that some of the largest-offset bursts may originate in the merger of double neutron star binaries that form dynamically in the cores of globular clusters.

We discuss mass transfer in eccentric binaries containing a white dwarf and a neutron star (WD--NS binaries). We show that such binaries are produced from field binaries following a series of mass transfer episodes that allow the white dwarf to form before the neutron star. We predict the orbital properties of binaries similar to the observed WD--NS binary J1141+6545, and show that they will undergo episodic mass transfer from the white dwarf to the neutron star. Furthermore, we describe oil-on-water, a two-phase SPH formalism that we have developed in order to model mass transfer in such binaries.

Aims: We seek to provide 2MASS and OGLE identifications and coordinates for all stars in the finding chart published by Arp\,(1965). This chart covers the low extinction area around NGC 6522, also known as Baade's window, at coordinates (l,b)=(1.02,-3.92). Methods: A cross correlation, using numerical techniques, was performed between a scan of the original finding chart from Arp (1965) and 2MASS and OGLE-II images and stellar coordinates. Results: We provide coordinates for all stars in Arp's finding chart and 2MASS and OGLE identifications wherever possible. Two identifications in quadrant II do not appear in the original finding chart.

We present the observed offsets of short-duration gamma-ray bursts (SGRBs) from their putative host galaxies and compare them to the expected distributions of merging compact object binaries, given the observed properties of the hosts. We find that for all but one burst in our sample the offsets are consistent with this model. For the case of bursts with massive elliptical host galaxies, the circular velocities of the hosts' haloes exceed the natal velocities of almost all our compact object binaries. Hence the extents of the predicted offset distributions for elliptical galaxies are determined largely by their spatial extents. In contrast, for spiral hosts the galactic rotation velocities are smaller than typical binary natal velocities and the predicted burst offset distributions are more extended than the galaxies. One SGRB, 060502B, apparently has a large radial offset that is inconsistent with an origin in a merging galactic compact binary. Although it is plausible that the host of GRB 060502B is mis-identified, our results show that the large offset is compatible with a scenario where at least a few per cent of SGRBs are created by the merger of compact binaries that form dynamically in globular clusters.

We review the observational evidence for extra mixing in stars on the red giant branch (RGB) and discuss why thermohaline mixing is a strong candidate mechanism. We recall the simple phenomenological description of thermohaline mixing, and aspects of mixing in stars in general. We use observations of M3 to constrain the form of the thermohaline diffusion coefficient and any associated free parameters. This is done by matching [C/Fe] and [N/Fe] along the RGB of M3. After taking into account a presumed initial primordial bimodality of [C/Fe] in the CN-weak and CN-strong stars our thermohaline mixing models can explain the full spread of [C/Fe]. Thermohaline mixing can produce a significant change in [N/Fe] as a function of absolute magnitude on the RGB for initially CN-weak stars, but not for initially CN-strong stars, which have so much nitrogen to begin with that any extra mixing does not significantly affect the surface nitrogen composition.

To date, there have been several detections of high-mass black hole binaries in both the Milky Way and other galaxies. For some of these, the spin parameter of the black hole has been estimated. As many of these systems are quite tight, a suggested origin of the spin is angular momentum imparted by the synchronous rotation of the black hole progenitor with its binary companion. Using Cygnus X-1, the best studied high-mass black hole binary, we investigate this possibility. We find that such an origin of the spin is not likely, and our results point rather to the spin being the result of processes during the collapse.

The formation mechanism of the barium stars is thought to be well understood. Barium-rich material, lost in a stellar wind from a thermally-pulsing asymptotic-giant branch star in a binary system, is accreted by its companion main-sequence star. Now, many millions of years later, the primary is an unseen white dwarf and the secondary has itself evolved into a giant which displays absorption lines of barium in its spectrum and is what we call a barium star. A similar wind-accretion mechanism is also thought to form the low-metallicity CH and carbon-enhanced metal-poor stars. Qualitatively the picture seems clear but quantitatively it is decidedly murky: several key outstanding problems remain which challenge our basic understanding of binary-star physics. Barium stars with orbital periods less than about 4,000 days should -- according to theory -- be in circular orbits because of tidal dissipation, yet they are often observed to be eccentric. Only one barium-star period longer than 10^4 days has been published although such stars are predicted to exist in large numbers. In this paper we attempt to shed light on these problems. First, we consider the impact of kicking the white dwarf at its birth, a notion which is supported by independent evidence from studies of globular clusters. Second, we increase the amount of orbital angular momentum loss during wind mass transfer, which shrinks barium-star binaries to the required period range. We conclude with a discussion of possible physical mechanisms and implications of a kick, such as the break up of wide barium-star binaries and the limits imposed on our models by observations.

We provide a brief review of thermohaline physics and why it is a candidate extra mixing mechanism during the red giant branch (RGB). We discuss how thermohaline mixing (also called $\delta$ $\mu$ mixing) during the RGB due to helium-3 burning, is more complicated than the operation of thermohaline mixing in other stellar contexts (such as following accretion from a binary companion). We try to use observations of carbon depletion in globular clusters to help constrain the formalism and the diffusion coefficient or mixing velocity that should be used in stellar models. We are able to match the spread of carbon depletion for metal poor field giants but are unable to do so for cluster giants, which may show evidence of mixing prior to even the first dredge-up event.

We consider whether stellar collisions can explain the observed depletion of red giants in the Galactic center. We model the stellar population with two different IMFs: 1) the Miller-Scalo and 2) a much flatter IMF. In the former case, low-mass main-sequence stars dominate the population, and collisions are unable to remove red giants out to 0.4 pc although brighter red giants much closer in may be depleted via collisions with stellar-mass black holes. For a much flatter IMF, the stellar population is dominated by compact remnants (i.e. black holes, white dwarfs and neutron stars). The most common collisions are then those between main-sequence stars and compact remnants. Such encounters are likely to destroy the main-sequence stars and thus prevent their evolution into red giants. In this way, the red-giant population could be depleted out to 0.4 pc matching observations. If this is the case, it implies the Galactic center contains a much larger population of stellar-mass black holes than would be expected from a regular IMF. This may in turn have implications for the formation and growth of the central supermassive black hole.

We present MONTAGE, a post-processing nucleosynthesis code that combines a traditional network for isotopes lighter than calcium with a rapid algorithm for calculating the s-process nucleosynthesis of the heavier isotopes. The separation of those parts of the network where only neutron-capture and beta-decay reactions are significant provides a substantial advantage in computational efficiency. We present the yields for a complete set of s-process isotopes for a 3 Mo, Z = 0.02 stellar model, as a demonstration of the utility of the approach. Future work will include a large grid of models suitable for use in calculations of Galactic chemical evolution.

There is an apparent dichotomy between the metal-poor ([Fe/H] <= -2) yet carbon-normal giants and their carbon-rich counterparts. The former undergo significant depletion of carbon on the red giant branch after they have undergone first dredge-up, whereas the latter do not appear to experience significant depletion. We investigate this in the context that the extra mixing occurs via the thermohaline instability that arises due to the burning of helium-3. We present the evolution of [C/Fe], [N/Fe] and carbon-12/carbon-13 for three models: a carbon-normal metal-poor star, and two stars that have accreted material from a 1.5 solar mass AGB companion, one having received 0.01 solar masses of material and the other having received 0.1 solar masses. We find the behaviour of the carbon-normal metal-poor stars is well reproduced by this mechanism. In addition, our models also show that the efficiency of carbon-depletion is significantly reduced in carbon-rich stars. This extra-mixing mechanism is able to reproduce the observed properties of both carbon-normal and carbon-rich stars.

To measure the onset of mass transfer in eccentric binaries we have developed a two-phase SPH technique. Mass transfer is important in the evolution of close binaries, and a key issue is to determine the separation at which mass transfer begins. The circular case is well understood and can be treated through the use of the Roche formalism. To treat the eccentric case we use a newly-developed two phase system. The body of the donor star is made up from high-mass "water" particles, whilst the atmosphere is modelled with low-mass "oil" particles. Both sets of particles take part fully in SPH interactions. To test the technique we model circular mass-transfer binaries containing a 0.6 Msun donor star and a 1 Msun white dwarf; such binaries are thought to form cataclysmic variable (CV) systems. We find that we can reproduce a reasonable CV mass-transfer rate, and that our extended atmosphere gives a separation that is too large by aproximately 16%, although its pressure scale height is considerably exaggerated. We use the technique to measure the semi-major axis required for the onset of mass transfer in binaries with a mass ratio of q=0.6 and a range of eccentricities. Comparing to the value obtained by considering the instantaneous Roche lobe at pericentre we find that the radius of the star required for mass transfer to begin decreases systematically with increasing eccentricity.

An N-body code containing live stellar evolution through combination of the software packages NBODY6 and STARS is presented. Operational details of the two codes are outlined and the changes that have been made to combine them discussed. We have computed the evolution of clusters of 10 000 stars using the combined code and we compare the results with those obtained using NBODY6 and the synthetic stellar evolution code SSE. We find that, providing the physics package within STARS is set up correctly to match the parameters of the models used to construct SSE, the results are very similar. This provides a good indication that the new code is working well. We also demonstrate how this physics can be changed simply in the new code with convective overshooting as an example. Similar changes in SSE would require considerable reworking of the model fits. We conclude by outlining proposed future development of the code to include more complete models of single stars and binary star systems.

We revisit the problem of why stars become red giants. We modify the physics of a standard stellar evolution code in order to determine what does and what does not contribute to a star becoming a red giant. In particular, we have run tests to try to separate the effects of changes in the mean molecular weight and in the energy generation. The implications for why stars become red giants are discussed. We find that while a change in the mean molecular weight is necessary (but not sufficient) for a 1 solar mass star to become a red giant, this is not the case in a star of 5 solar masses. It therefore seems that there may be more than one way to make a giant.

We show that collisions with stellar--mass black holes can partially explain the absence of bright giant stars in the Galactic Centre, first noted by Genzel et al, 1996. We show that the missing objects are low--mass giants and AGB stars in the range 1-3 M$_{\odot}$. Using detailed stellar evolution calculations, we find that to prevent these objects from evolving to become visible in the depleted K bands, we require that they suffer collisions on the red giant branch, and we calculate the fractional envelope mass losses required. Using a combination of Smoothed Particle Hydrodynamic calculations, restricted three--body analysis and Monte Carlo simulations, we compute the expected collision rates between giants and black holes, and between giants and main--sequence stars in the Galactic Centre. We show that collisions can plausibly explain the missing giants in the $10.5<K<12$ band. However, depleting the brighter ($K<10.5$) objects out to the required radius would require a large population of black hole impactors which would in turn deplete the $10.5<K<12$ giants in a region much larger than is observed. We conclude that collisions with stellar--mass black holes cannot account for the depletion of the very brightest giants, and we use our results to place limits on the population of stellar--mass black holes in the Galactic Centre.

Massive black holes (MBHs) with a mass below ~ 1e7 Msun are likely to reside at the centre of dense stellar nuclei shaped by 2-body relaxation, close interactions with the MBH and direct collisions. In this contribution, we stress the role of mass segregation of stellar-mass black holes into the innermost tenths of a parsec and point to the importance of hydrodynamical collisions between stars. At the Galactic centre, collisions must affect giant stars and some of the S-stars.

We have computed the galactic trajectories of twelve hypervelocity stars (HVSs) under the assumption that they originated in the Galactic Centre. We show that eight of these twelve stars are bound to the Galaxy. We consider the subsequent trajectories of the bound stars to compute their characteristic orbital period, which is 2 Gyr. All eight bound stars are moving away from the centre of the Galaxy, which implies that the stars' lifetimes are less than 2 Gyr. We thus infer that the observed HVSs are massive main sequence stars, rather than blue horizontal branch stars. The observations suggest that blue HVSs are ejected from the Galactic Centre roughly every 15 Myr. This is consistent with the observed population of blue stars in extremely tight orbits round the central super-massive black hole (SMBH), the so-called S-stars, if we assume that the HVSs are produced by the breakup of binaries. One of the stars in such a binary is ejected at high velocities to form a HVS; the other remains bound to the SMBH as an S-star. We further show that the one high-velocity system observed to be moving towards the Galactic Centre, SDSS J172226.55+594155.9, could not have originated in the Galactic Centre; rather, we identify it as a halo object.

Most stars are formed in a cluster or association, where the number density of stars can be high. This means that a large fraction of initially-single stars will undergo close encounters with other stars and/or exchange into binaries. We describe how such close encounters and exchange encounters can affect the properties of a planetary system around a single star. We define a singleton as a single star which has never suffered close encounters with other stars or spent time within a binary system. It may be that planetary systems similar to our own solar system can only survive around singletons. Close encounters or the presence of a stellar companion will perturb the planetary system, often leaving planets on tighter and more eccentric orbits. Thus planetary systems which initially resembled our own solar system may later more closely resemble some of the observed exoplanet systems.

The stars that populate the solar neighbourhood were formed in stellar clusters. Through N-body simulations of these clusters, we measure the rate of close encounters between stars. By monitoring the interaction histories of each star, we investigate the singleton fraction in the solar neighbourhood. A singleton is a star which formed as a single star, has never experienced any close encounters with other stars or binaries, or undergone an exchange encounter with a binary. We find that, of the stars which formed as single stars, a significant fraction are not singletons once the clusters have dispersed. If some of these stars had planetary systems, with properties similar to those of the solar system, the planets orbits may have been perturbed by the effects of close encounters with other stars or the effects of a companion star within a binary. Such perturbations can lead to strong planet-planet interactions which eject several planets, leaving the remaining planets on eccentric orbits. Some of the single stars exchange into binaries. Most of these binaries are broken up via subsequent interactions within the cluster, but some remain intact beyond the lifetime of the cluster. The properties of these binaries are similar to those of the observed binary systems containing extra-solar planets. Thus, dynamical processes in young stellar clusters will alter significantly any population of solar-system-like planetary systems. In addition, beginning with a population of planetary systems exactly resembling the solar system around single stars, dynamical encounters in young stellar clusters may produce at least some of the extra-solar planetary systems observed in the solar neighbourhood.