K2-19 hosts a planetary system composed of two outer planets, b and c, with size of $7.0\pm 0.2~R_\oplus$ and $4.1\pm0.2~R_\oplus$ , and an inner planet, d, with a radius of $1.11\pm 0.05 R_\oplus$. A recent analysis of Transit-Timing Variations (TTVs) suggested b and c are close to but not in 3:2 mean motion resonance (MMR) because the classical resonant angles circulate. Such an architecture challenges our understanding of planet formation. Indeed, planet migration through the protoplanetary disk should lead to a capture into the MMR. Here, we show that the planets are in fact, locked into the 3:2 resonance despite circulation of the conventional resonant angles and aligned periapses. However, we show that such an orbital configuration cannot be maintained for more than a few hundred million years due to the tidal dissipation experienced by planet d. The tidal dissipation remains efficient because of a secular forcing of the innermost planet eccentricity by planets b and c. While the observations strongly rule out an orbital solution where the three planets are on close to circular orbits, it remains possible that a fourth planet is affecting the TTVs such that the four planet system is consistent with the tidal constraints.

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.

The dynamical stability of tightly packed exoplanetary systems remains poorly understood. While for a two-planet system a sharp stability boundary exists, numerical simulations of three and more planet systems show that they can experience instability on timescales up to billions of years. Moreover, an exponential trend between the planet orbital separation measured in units of Hill radii and the survival time has been reported. While these findings have been observed in numerous numerical simulations, little is known of the actual mechanism leading to instability. Contrary to a constant diffusion process, planetary systems seem to remain dynamically quiescent for most of their lifetime before a very short unstable phase. In this work, we show how the slow chaotic diffusion due to the overlap of three-body resonances dominates the timescale leading to the instability for initially coplanar and circular orbits. While the last instability phase is related to scattering due to two-planet mean motion resonances (MMR), for circular orbits the two-planets MMR are too far separated to destabilize systems initially away from them. We develop an analytical model to generalize the empirical trend obtained for equal mass and equally-spaced planets to general systems. We obtain an analytical estimate of the survival time consistent with simulations over four orders of magnitude for the planet to star mass ratio $\epsilon$, and 6 to 8 orders of magnitude for the instability time. We also confirm that measuring the orbital spacing in terms of Hill radii is not adapted and that the right spacing unit scales as $\epsilon^{1/4}$. We predict that beyond a certain spacing, the three-planet resonances are not overlapped, which results in an increase of the survival time. We finally discuss the extension of our result to more general systems, containing more planets on initially non circular orbits.

Single-binary scattering may lead to an exchange where the single object captures a component of the binary, forming a new binary. This has been well studied in encounters between a star--planet pair and a single star. Here we explore the application of the exchange mechanism to a planet--satellite pair and another planet in the gravitational potential of a central star. As a case study, we focus on encounters between a satellite-bearing object and Neptune. We investigate whether Neptune can capture satellites from that object and if the captured satellites have orbits analogous to the Neptunian moons Triton and Nereid. Using $N$-body simulations, we study the capture probability at different encounter distances. Post-capture, we use a simple analytical argument to estimate how the captured orbits evolve under collisional and tidal effects. We find that the average capture probability reaches $\sim$$10\%$ if Neptune penetrates the donor planet's satellite system. Most moons grabbed by Neptune acquire highly eccentric orbits. Post-capture, around half of those captured, especially those on tight orbits, can be circularised, either by tides only or by collisions+tides, turning into Triton-like objects. Captures further out, on the other hand, stay on wide and eccentric orbits like that of Nereid. Both moon types can be captured in the same encounter and they have wide distributions in orbital inclination. Therefore, Triton naturally has a $\sim$50\% chance of being retrograde. A similar process potentially applies to an exoplanetary system, and our model predicts that exomoons can jump from one planet to another during planetary scattering. Specifically, there should be two distinct populations of captured moons: one on close-in circular orbits and the other on far-out eccentric orbits. The two populations may have highly inclined prograde or retrograde orbits.

Planetary systems formed in clusters may be subject to stellar encounter flybys. Here we create a diverse range of representative planetary systems with different orbital scales and planets' masses and examine encounters between them in a typical open cluster. We first explore the close-in multi-super earth systems $\lesssim0.1$ au. They are resistant to flybys in that only ones inside a few au can destabilise a planet or break the resonance between such planets. But these systems may capture giant planets onto wide orbits from the intruding star during distant flybys. If so, the original close-in small planets' orbits may be tilted together through Kozai--Lidov mechanism, forming a ``cold'' system that is significantly inclined against the equator of the central host. Moving to the intermediately-placed planets around solar-like stars, we find that the planets' mass gradient governs the systems' long-term evolution post-encounter: more massive planets have better chances to survive. Also, a system's angular momentum deficit, a quantity describing how eccentric/inclined the orbits are, measured immediately after the encounter, closely relates to the longevity of the systems -- whether or not and when the systems turn unstable in the ensuing evolution millions of years post-encounter. We compare the orbits of the surviving planets in the unstable systems through (1) the immediate consequence of the stellar fly or (2) internal interplanetary scattering long post-encounter and find that those for the former are systematically colder. Finally, we show that massive wide-orbit multi-planet systems like that of HR 8799 can be easily disrupted and encounters at a few hundreds of au suffice.

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.

The forthcoming Laser Interferometer Space Antenna (LISA) will probe the population of coalescing massive black hole (MBH) binaries up to the onset of structure formation. Here we simulate the galactic-scale pairing of $\sim10^6 M_\odot$ MBHs in a typical, non-clumpy main-sequence galaxy embedded in a cosmological environment at $z = 7-6$. In order to increase our statistical sample, we adopt a strategy that allows us to follow the evolution of six secondary MBHs concomitantly. We find that the magnitude of the dynamical-friction induced torques is significantly smaller than that of the large-scale, stochastic gravitational torques arising from the perturbed and morphologically evolving galactic disc, suggesting that the standard dynamical friction treatment is inadequate for realistic galaxies at high redshift. The dynamical evolution of MBHs is very stochastic, and a variation in the initial orbital phase can lead to a drastically different time-scale for the inspiral. Most remarkably, the development of a galactic bar in the host system either significantly accelerates the inspiral by dragging a secondary MBH into the centre, or ultimately hinders the orbital decay by scattering the MBH in the galaxy outskirts. The latter occurs more rarely, suggesting that galactic bars overall promote MBH inspiral and binary coalescence. The orbital decay time can be an order of magnitude shorter than what would be predicted relying on dynamical friction alone. The stochasticity, and the important role of global torques, have crucial implications for the rates of MBH coalescences in the early Universe: both have to be accounted for when making predictions for the upcoming LISA observatory.

We report the detection of a transiting super-Earth-sized planet (R=1.39+-0.09 Rearth) in a 1.4-day orbit around L 168-9 (TOI-134),a bright M1V dwarf (V=11, K=7.1) located at 25.15+-0.02 pc. The host star was observed in the first sector of the Transiting Exoplanet Survey Satellite (TESS) mission and, for confirmation and planet mass measurement, was followed up with ground-based photometry, seeing-limited and high-resolution imaging, and precise radial velocity (PRV) observations using the HARPS and PFS spectrographs. Combining the TESS data and PRV observations, we find the mass of L168-9 b to be 4.60+-0.56 Mearth, and thus the bulk density to be 1.74+0.44-0.33 times larger than that of the Earth. The orbital eccentricity is smaller than 0.21 (95% confidence). This planet is a Level One Candidate for the TESS Mission's scientific objective - to measure the masses of 50 small planets - and is one of the most observationally accessible terrestrial planets for future atmospheric characterization.

We investigate the possibility of finding Earth-like planets in the habitable zone of 34 nearby FGK-dwarfs, each known to host one giant planet exterior to their habitable zone detected by RV. First we simulate the dynamics of the planetary systems in their present day configurations and determine the fraction of stable planetary orbits within their habitable zones. Then, we postulate that the eccentricity of the giant planet is a result of an instability in their past during which one or more other planets were ejected from the system. We simulate these scenarios and investigate whether planets orbiting in the habitable zone survive the instability. Explicitly we determine the fraction of test particles, originally found in the habitable zone, which remain in the habitable zone today. We label this fraction the resilient habitability of a system. We find that for most systems the probability of planets existing [or surviving] on stable orbits in the habitable zone becomes significantly smaller when we include a phase of instability in their history. We present a list of candidate systems with high resilient habitability for future observations. These are: HD 95872, HD 154345, HD 102843, HD 25015, GJ 328, HD 6718 and HD 150706. The known planets in the last two systems have large observational uncertainties on their eccentricities, which propagate into large uncertainties on their resilient habitability. Further observational constraints of these two eccentriciti

The atmospheres of close-in planets are strongly influenced by mass loss driven by the high-energy (X-ray and extreme ultraviolet, EUV) irradiation of the host star, particularly during the early stages of evolution. We recently developed a framework to exploit this connection and enable us to recover the past evolution of the stellar high-energy emission from the present-day properties of its planets, if the latter retains some remnants of their primordial hydrogen-dominated atmospheres. Furthermore, the framework can also provide constraints on planetary initial atmospheric mass fractions. The constraints on the output parameters improve when more planets can be simultaneously analysed. This makes the Kepler-11 system, which hosts six planets with bulk densities between 0.66 and 2.45g cm^{-3}, an ideal target. Our results indicate that the star has likely evolved as a slow rotator (slower than 85\% of the stars with similar masses), corresponding to a high-energy emission at 150 Myr of between 1-10 times that of the current Sun. We also constrain the initial atmospheric mass fractions for the planets, obtaining a lower limit of 4.1% for planet c, a range of 3.7-5.3% for planet d, a range of 11.1-14% for planet e, a range of 1-15.6% for planet f, and a range of 4.7-8.7% for planet g assuming a disc dispersal time of 1 Myr. For planet b, the range remains poorly constrained. Our framework also suggests slightly higher masses for planets b, c, and f than have been suggested based on transit timing variation measurements. We coupled our results with published planet atmosphere accretion models to obtain a temperature (at 0.25 AU, the location of planet f) and dispersal time of the protoplanetary disc of 550 K and 1 Myr, although these results may be affected by inconsistencies in the adopted system parameters.

Statistical analyses from exoplanet surveys around low-mass stars indicate that super-Earth and Neptune-mass planets are more frequent than gas giants around such stars, in agreement with core accretion theory of planet formation. Using precise radial velocities derived from visual and near-infrared spectra, we report the discovery of a giant planet with a minimum mass of 0.46 Jupiter masses in an eccentric 204-day orbit around the very low-mass star GJ 3512. Dynamical models show that the high eccentricity of the orbit is most likely explained from planet-planet interactions. The reported planetary system challenges current formation theories and puts stringent constraints on the accretion and migration rates of planet formation and evolution models, indicating that disc instability may be more efficient in forming planets than previously thought.

Most giant exoplanets discovered by radial velocity surveys have much higher eccentricities than those in the solar system. The planet--planet scattering mechanism has been shown to match the broad eccentricity distribution, but the highest-eccentricity planets are often attributed to Kozai-Lidov oscillations induced by a stellar companion. Here we investigate whether the highly eccentric exoplanet population can be produced entirely by scattering. We ran 500 N-body simulations of closely packed giant-planet systems that became unstable under their own mutual perturbations. We find that the surviving bound planets can have eccentricities up to $e > 0.99$, with a maximum of 0.999017 in our simulations. This suggests that there is no maximum eccentricity that can be produced by planet--planet scattering. Importantly, we find that extreme eccentricities are not extremely rare; the eccentricity distribution for all giant exoplanets with $e > 0.3$ is consistent with all planets concerned being generated by scattering. Our results show that the discovery of planets with extremely high eccentricities does not necessarily signal the action of the Kozai-Lidov mechanism.

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.

Wide field surveys for transiting planets are well suited to searching diverse stellar populations, enabling a better understanding of the link between the properties of planets and their parent stars. We report the discovery of HAT-P-69b (TOI 625.01) and HAT-P-70b (TOI 624.01), two new hot Jupiters around A stars from the HATNet survey which have also been observed by the Transiting Exoplanet Survey Satellite (TESS). HAT-P-69b has a mass of 3.58 +0.58/-0.58 MJup and a radius of 1.676 +0.051/-0.033 RJup, residing in a prograde 4.79-day orbit. HAT-P-70b has a radius of 1.87 +0.15/-0.10 RJup and a mass constraint of < 6.78 (3 sigma) MJup, and resides in a retrograde 2.74-day orbit. We use the confirmation of these planets around relatively massive stars as an opportunity to explore the occurrence rate of hot Jupiters as a function of stellar mass. We define a sample of 47,126 main-sequence stars brighter than Tmag=10 that yields 31 giant planet candidates, including 18 confirmed planets, 3 candidates, and 10 false positives. We find a net hot Jupiter occurrence rate of 0.41+/-0.10 % within this sample, consistent with the rate measured by Kepler for FGK stars. When divided into stellar mass bins, we find the occurrence rate to be 0.71+/-0.31% for G stars, 0.43+/-0.15% for F stars, and 0.26+/-0.11% for A stars. Thus, at this point, we cannot discern any statistically significant trend in the occurrence of hot Jupiters with stellar mass.

TESS is revolutionising the search for planets orbiting bright and nearby stars. In sectors 3 and 4, TESS observed TOI-402 (TIC-120896927), a bright V=9.1 K1 dwarf also known as HD 15337, and found two transiting signals with period of 4.76 and 17.18 days and radius of 1.90 and 2.21\,\Rearth. This star was observed as part of the radial-velocity search for planets using the HARPS spectrometer, and 85 precise radial-velocity measurements were obtained over a period of 14 years. In this paper, we analyse the HARPS radial-velocity measurements in hand to confirm the planetary nature of these two signals. By reanalysing TESS photometry and host star parameters using EXOFASTv2, we find that TOI-402.01 and TOI-402.02 have periods of 4.75642$\pm$0.00021 and 17.1784$\pm$0.0016 days and radii of 1.70$\pm$0.06 and 2.52$\pm$0.11\,\Rearth\,(precision 3.6 and 4.2\%), respectively. By analysing the HARPS radial-velocity measurements, we find that those planets are both super-Earths with masses of 7.20$\pm$0.81 and 8.79$\pm$1.67\,\Mearth\,(precision 11.3 and 19.0\%), and small eccentricities compatible with zero at 2$\sigma$. Although having rather similar masses, the radius of these two planets is really different, putting them on different sides of the radius gap. With stellar irradiation 160 times more important than Earth for TOI-402.01 and only 29 times more for TOI-402.02, it is likely that photo-evaporation is at the origin of this radius difference. Those two planets, being in the same system and therefore being in the same irradiation environment are therefore extremely important to perform comparative exoplanetology across the evaporation valley and thus bring constraints on the mechanisms responsible for the radius gap.

Stars formed in clusters can encounter other stars at close distances. In typical open clusters in the Solar neighbourhood containing hundreds or thousands of member stars, ten to twenty per cent of Solar-mass member stars are expected to encounter another star at distances closer than 100 au. These close encounters strongly perturb the planetary systems, directly causing ejection of planets or their capture by the intruding star, as well as exciting the orbits. Using extensive $N$-body simulations, we study such fly-by encounters between two Solar System analogues, each with four giant planets from Jupiter to Neptune. We quantify the rates of loss and capture immediately after the encounter, e.g., the Neptune analogue is lost in one in four encounters within 100 au, and captured by the flying-by star in one in twelve encounters. We then perform long-term (up to 1 Gyr) simulations investigating the ensuing post-encounter evolution. We show that large numbers of planets are removed from systems due to planet--planet interactions and that captured planets further enhance the system instability. While encounters can initially leave a planetary system containing more planets by inserting additional ones, the long-term instability causes a net reduction in planet number. A captured planet ends up on a retrograde orbit in half of the runs in which it survives for 1 Gyr; also, a planet bound to its original host star but flipped during the encounter may survive. Thus, encounters between planetary systems are a channel to create counter-rotating planets, This would happen in around 1% of systems, and such planets are potentially detectable through astrometry or direct imaging.

Young exoplanets are snapshots of the planetary evolution process. Planets that orbit stars in young associations are particularly important because the age of the planetary system is well constrained. We present the discovery of a transiting planet larger than Neptune but smaller than Saturn in the 45 Myr Tucana-Horologium young moving group. The host star is a visual binary, and our follow-up observations demonstrate that the planet orbits the G6V primary component, DS Tuc A (HD 222259A, TIC 410214986). We first identified transits using photometry from the Transiting Exoplanet Survey Satellite (TESS; alerted as TOI 200.01). We validated the planet and improved the stellar parameters using a suite of new and archival data, including spectra from SOAR/Goodman, SALT/HRS and LCO/NRES; transit photometry from Spitzer; and deep adaptive optics imaging from Gemini/GPI. No additional stellar or planetary signals are seen in the data. We measured the planetary parameters by simultaneously modeling the photometry with a transit model and a Gaussian process to account for stellar variability. We determined that the planetary radius is $5.70\pm0.17$ Earth radii and that the orbital period is 8.1 days. The inclination angles of the host star's spin axis, the planet's orbital axis, and the visual binary's orbital axis are aligned within 15 degrees to within the uncertainties of the relevant data. DS Tuc Ab is bright enough (V=8.5) for detailed characterization using radial velocities and transmission spectroscopy.

Our aim is to investigate the nature and formation of brown dwarfs by adding a new well-characterised object to the small sample of less than 20 transiting brown dwarfs. One brown dwarf candidate was found by the KESPRINT consortium when searching for exoplanets in the K2 space mission Campaign 16 field. We combined the K2 photometric data with a series of multi-colour photometric observations, imaging and radial velocity measurements to rule out false positive scenarios and to determine the fundamental properties of the system. We report the discovery and characterisation of a transiting brown dwarf in a 5.17 day eccentric orbit around the slightly evolved F7V star EPIC 212036875. We find a stellar mass of 1.15+/-0.08 M$_\odot$, a stellar radius of 1.41+/-0.05 R$_\odot$, and an age of 5.1+/-0.9 Gyr. The mass and radius of the companion brown dwarf are 51+/-2 MJ and 0.83+/-0.03 RJ, respectively, corresponding to a mean density of 108+15-13 g cm-3. EPIC 212036875 b is a rare object that resides in the brown dwarf desert. In the mass-density diagram for planets, brown dwarfs and stars, we find that all giant planets and brown dwarfs follow the same trend from ~0.3 MJ to the turn-over to hydrogen burning stars at ~73 MJ. EPIC 212036875 b falls close to the theoretical model for mature H/He dominated objects in this diagram as determined by interior structure models, as well as the empirical fit. We argue that EPIC 212036875 b formed via gravitational disc instabilities in the outer part of the disc, followed by a quick migration. Orbital tidal circularisation may have started early in its history for a brief period when the brown dwarf's radius was larger. The lack of spin-orbit synchronisation points to a weak stellar dissipation parameter which implies a circularisation timescale of >23 Gyr, or suggests an interaction between the magnetic and tidal forces of the star and the brown dwarf.

We report the discovery of a super-Earth and a sub-Neptune transiting the star HD 15337 (TOI-402, TIC 120896927), a bright (V=9) K1 dwarf observed by the Transiting Exoplanet Survey Satellite (TESS) in Sectors 3 and 4. We combine the TESS photometry with archival HARPS spectra to confirm the planetary nature of the transit signals and derive the masses of the two transiting planets. With an orbital period of 4.8 days, a mass of 7.51(+1.09)(-1.01) M_Earth, and a radius of 1.64+/-0.06 R_Earth, HD 15337b joins the growing group of short-period super-Earths known to have a rocky terrestrial composition. The sub-Neptune HD 15337c has an orbital period of 17.2 days, a mass of 8.11(+1.82)(-1.69) M_Earth, and a radius of 2.39+/-0.12 R_Earth, suggesting that the planet might be surrounded by a thick atmospheric envelope. The two planets have similar masses and lie on opposite sides of the radius gap, and are thus an excellent testbed for planet formation and evolution theories. Assuming that HD 15337c hosts a hydrogen-dominated envelope, we employ a recently developed planet atmospheric evolution algorithm in a Bayesian framework to estimate the history of the high-energy (extreme ultraviolet and X-ray) emission of the host star. We find that at an age of 150 Myr, the star possessed on average between 3.7 and 127 times the high-energy luminosity of the current Sun.

Many white dwarf stars show signs of having accreted smaller bodies, implying that they may host planetary systems. A small number of these systems contain gaseous debris discs, visible through emission lines. We report a stable 123.4min periodic variation in the strength and shape of the CaII emission line profiles originating from the debris disc around the white dwarf SDSSJ122859.93+104032.9. We interpret this short-period signal as the signature of a solid body held together by its internal strength.

Giant molecular clouds (GMCs) are believed to affect the biospheres of planets as their host star passes through them. We simulate the trajectories of stars and GMCs in the Galaxy and determine how often stars pass through GMCs. We find a strong decreasing dependence with Galactocentric radius, and with the velocity perpendicular to the Galactic plane, $V_z$. The \textit{XY}-component of the kinematic heating of stars was shown to not affect the GMC hit rate, unlike the \textit{Z}-dependence ($V_z$) implies that stars hit fewer GMCs as they age. GMCs are locations of star formation, therefore we also determine how often stars pass near supernovae. For the supernovae the decrease with $V_z$ is steeper as how fast the star passes through the GMC determines the probability of a supernova encounter. We then integrate a set of Sun-like trajectories to see the implications for the Sun. We find that the Sun hits $1.6\pm1.3$ GMCs per Gyr which results in $1.5\pm1.1$ or (with correction for clustering) $0.8\pm 0.6$ supernova closer than 10 pc per Gyr. The different supernova frequencies are from whether one considers multiple supernova per GMC crossing (few Myr) as separate events. We then discuss the effect of the GMC hits on the Oort cloud, and the Earth's climate due to accretion, we also discuss the records of distant supernova. Finally, we determine Galactic Habitable Zone using our model. For the thin disk we find it to lie between 5.8-8.7 kpc and for the thick disk to lie between 4.5-7.7 kpc.

The globular clusters of large spiral galaxies can be divided into two populations: one which formed in-situ and one which comprises clusters tidally stripped away from other galaxies. In this paper we investigate the contribution to the outer globular cluster population in the M31 galaxy through donation of clusters from dwarf galaxies. We test this numerically by comparing the contribution of globular clusters from simulated encounters to the observed M31 globular cluster population. To constrain our simulations, we specifically investigate the outermost globular cluster in the M31 system, MGC1. The remote location of MGC1 favours the idea of it being captured, however, the cluster is devoid of features associated with tidal interactions. Hence we separate simulations where tidal features are present and where they are hidden. We find that our simulated encounters can place clusters on MGC1-like orbits. In addition, we find that tidal stripping of clusters from dwarf galaxies leaves them on orbits having a range of separations, broadly matching those observed in M31. We find that the specific energies of globular clusters captured by M31 closely matches those of the incoming host dwarf galaxies. Furthermore, in our simulations we find an equal number of accreted clusters on co-rotating and counter-rotating orbits within M31 and use this to infer the fraction of clusters that has been accreted. We find that even close in roughly 50% of the clusters are accreted, whilst this figure increases to over 80% further out.

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.

The first neutron star-neutron star merger was detected in by LIGO/Virgo in a galaxy in which the majority of star formation was taking place a long time ago (11 Gyr). LIGO/Virgo estimated that local cosmic NS-NS merger rate is 110-3840 Gpc^-3 yr^-1. Only some extreme evolutionary models can generate NS-NS merger rates in old host galaxies consistent with the LIGO/Virgo estimate. However, we show that these models generate rates exceeding empirical Galactic NS-NS merger rates based on the large population of Milky Way NS-NS binaries. Typically, current evolutionary models produce NS-NS merger rates that are consistent with the Milky Way empirical rates. However, these models generate local cosmic NS-NS merger rate in old host galaxies that are below the LIGO/Virgo estimate. The reason behind this tension is the predicted delay time distribution between star formation and NS-NS mergers that favors short delays. Evolutionary models produce a generic steep power-law NS-NS delay time distribution. This limits NS-NS merger rates in old host galaxies. However, we show that such distribution is consistent with observations of Galactic NS-NS binaries; 50% of which show very long merger times (much longer than Hubble time). Once model distributions are convolved with continuous prolonged (10 Gyr) star formation in the Galactic disk, then 20-70% (depending on a model) of the predicted NS-NS population has very long current Galactic merger times. Although NS-NS binaries are formed predominantly with short delay times, many of short delay time systems merge and do not make it to the present, while long delay time systems survive and contribute to the current Galactic NS-NS population. This study highlights the tension between the current evolutionary predictions and the observation of the first NS-NS merger in an old host galaxy.

We report on the confirmation of a transiting giant planet around the relatively hot (Teff = 6801 $\pm$ 56 K) star HD2685, whose transit signal was detected in Sector 1 data of the TESS mission. We confirmed the planetary nature of the transit signal by using Doppler velocimetric measurements with CHIRON, CORALIE and FEROS, as well as photometric data with CHAT and LCOGT. From the photometry and radial velocities joint analysis, we derived the following parameters for HD2685 $b$: $P$=4.12692$\pm$0.00004 days, M$_P$=1.18 $\pm$ 0.09 $M_J$ and $R_P$=1.44 $\pm$ 0.01 $R_J$. This system is a typical example of an inflated transiting Hot-Jupiter in a circular orbit. Given the host star apparent visual magnitude ($V$ = 9.6 mag), this is one of the brightest known stars hosting a transiting Hot-Jupiter, and a good example of the upcoming systems that will be detected by TESS during the two-year primary mission. This is also an excellent target for future ground and space based atmospheric characterization as well as a good candidate for measuring the projected spin-orbit misalignment angle via the Rossiter-McLaughlin effect.

Data from the newly-commissioned \textit{Transiting Exoplanet Survey Satellite} (TESS) has revealed a "hot Earth" around LHS 3844, an M dwarf located 15 pc away. The planet has a radius of $1.32\pm 0.02$ $R_\oplus$ and orbits the star every 11 hours. Although the existence of an atmosphere around such a strongly irradiated planet is questionable, the star is bright enough ($I=11.9$, $K=9.1$) for this possibility to be investigated with transit and occultation spectroscopy. The star's brightness and the planet's short period will also facilitate the measurement of the planet's mass through Doppler spectroscopy.

Context. Gaia DR2 offers unparalleled precision on stars' parallaxes and proper motions. This allows the prediction of microlensing events for which the lens stars (and any planets they possess) are nearby and may be well studied and characterised. Aims. We identify a number of potential microlensing events that will occur before the year 2035.5, 20 years from the Gaia DR2 reference epoch. Methods. We query Gaia DR2 for potential lenses within 100 pc, extract parallaxes and proper motions of the lenses and background sources, and identify potential lensing events. We estimate the lens masses from Priam effective temperatures, and use these to calculate peak magnifications and the size of the Einstein radii relative to the lens stars' habitable zones. Results. We identify 7 future events with a probability > 10% of an alignment within one Einstein radius. Of particular interest is DR2 5918299904067162240 (WISE J175839.20-583931.6), magnitude G = 14.9, which will lens a G = 13.9 background star in early 2030, with a median 23% net magnification. Other pairs are typically fainter, hampering characterisation of the lens (if the lens is faint) or the ability to accurately measure the magnification (if the source is much fainter than the lens). Of timely interest is DR2 4116504399886241792 (2MASS J17392440-2327071), which will lens a background star in July 2020, albeit with weak net magnification (0.03%). Median magnifications for the other 5 high-probability events range from 0.3% to 5.3%. The Einstein radii for these lenses are 1-10 times the radius of the habitable zone, allowing these lensing events to pick out cold planets around the ice line, and filling a gap between transit and current microlensing detections of planets around very low-mass stars. Conclusions. We provide a catalogue of the predicted events to aid future characterisation efforts... [abridged]

Planet-planet collisions are a common outcome of instability in systems of transiting planets close to the star, as well as occurring during in-situ formation of such planets from embryos. Previous N-body studies of instability amongst transiting planets have assumed that collisions result in perfect merging. Here, we explore the effects of implementing a more realistic collision prescription on the outcomes of instability and in-situ formation at orbital radii of a few tenths of an au. There is a strong effect on the outcome of the growth of planetary embryos, so long as the debris thrown off in collisions is rapidly removed from the system (which happens by collisional processing to dust, and then removal by radiation forces) and embryos are small ($< 0.1\mathrm{\,M}_\oplus$). If this is the case, then systems form fewer detectable ($\gtrsim1\mathrm{\,M}_\oplus$) planets than systems evolved under the assumption of perfect merging in collisions. This provides some contribution to the "Kepler Dichotomy": the observed over-abundance of single-planet systems. The effects of changing the collision prescription on unstable mature systems of super-Earths are less pronounced. Perfect mergers only account for a minority of collision outcomes in such systems, but most collisions resulting in mass loss are grazing impacts in which only a few per cent. of mass is lost. As a result, there is little impact on the final masses and multiplicities of the systems after instability when compared to systems evolved under the assumption that collisions always result in perfect merging.

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.

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

The star EPIC 210894022 has been identified from a light curve acquired through the K2 space mission as possibly orbited by a transiting planet. Our aim is to confirm the planetary nature of the object and derive its fundamental parameters. We combine the K2 photometry with reconnaissance spectroscopy and radial velocity (RV) measurements obtained using three separate telescope and spectrograph combinations. The spectroscopic synthesis package SME has been used to derive the stellar photospheric parameters that were used as input to various stellar evolutionary tracks in order to derive the parameters of the system. The planetary transit was also validated to occur on the assumed host star through adaptive imaging and statistical analysis. The star is found to be located in the background of the Hyades cluster at a distance at least 4 times further away from Earth than the cluster itself. The spectrum and the space velocities of EPIC 210894022 strongly suggest it to be a member of the thick disk population. We find that the star is a metal poor ([Fe/H]=-0.53+/-0.05 dex) and alpha-rich somewhat evolved solar-like object of spectral type G3 with Teff=5730+/-50 K, logg=4.15+/-0.1 (cgs), radius of 1.3+/-0.1 R_Sun, and mass of 0.88+/-0.02 M_Sun. The RV detection together with the imaging confirms with a high level of significance that the transit signature is caused by a super-Earth orbiting the star EPIC 210894022. We measure a mass of 8.6+/-3.9 M_Earth and a radius of 1.9+/-0.2 R_Earth. A second more massive object with a period longer than about 120 days is indicated by a long term linear acceleration. With an age of > 10 Gyrs this system is one of the oldest where planets is hitherto detected. Further studies of this planetary system is important since it contains information about the planetary formation process during a very early epoch of the history of our Galaxy.

The distribution of exoplanet masses is not primordial. After the initial stage of planet formation is complete, the gravitational interactions between planets can lead to the physical collision of two planets, or the ejection of one or more planets from the system. When this occurs, the remaining planets are typically left in more eccentric orbits. Here we use present-day eccentricities of the observed exoplanet population to reconstruct the initial mass function of exoplanets before the onset of dynamical instability. We developed a Bayesian framework that combines data from N-body simulations with present-day observations to compute a probability distribution for the planets that were ejected or collided in the past. Integrating across the exoplanet population, we obtained an estimate of the initial mass function of exoplanets. We find that the ejected planets are primarily sub-Saturn type planets. While the present-day distribution appears to be bimodal, with peaks around $\sim 1 M_{\rm J}$ and $\sim 20 M_\oplus$, this bimodality does not seem to be primordial. Instead, planets around $\sim 60 M_\oplus$ appear to be preferentially removed by dynamical instabilities. Attempts to reproduce exoplanet populations using population synthesis codes should be mindful of the fact that the present population has been depleted of intermediate-mass planets. Future work should explore how the system architecture and multiplicity might alter our results.

Recent years have seen growing interest in the streaming instability as a candidate mechanism to produce planetesimals. However, these investigations have been limited to small-scale simulations. We now present the results of a global protoplanetary disk evolution model that incorporates planetesimal formation by the streaming instability, along with viscous accretion, photoevaporation by EUV, FUV, and X-ray photons, dust evolution, the water ice line, and stratified turbulence. Our simulations produce massive (60-130 $M_\oplus$) planetesimal belts beyond 100 au and up to $\sim 20 M_\oplus$ of planetesimals in the middle regions (3-100 au). Our most comprehensive model forms 8 $M_\oplus$ of planetesimals inside 3 au, where they can give rise to terrestrial planets. The planetesimal mass formed in the inner disk depends critically on the timing of the formation of an inner cavity in the disk by high-energy photons. Our results show that the combination of photoevaporation and the streaming instability are efficient at converting the solid component of protoplanetary disks into planetesimals. Our model, however, does not form enough early planetesimals in the inner and middle regions of the disk to give rise to giant planets and super-Earths with gaseous envelopes. Additional processes such as particle pileups and mass loss driven by MHD winds may be needed to drive the formation of early planetesimal generations in the planet forming regions of protoplanetary disks.

We study how close-in systems such as those detected by Kepler are affected by the dynamics of bodies in the outer system. We consider two scenarios: outer systems of giant planets potentially unstable to planet--planet scattering, and wide binaries that may be capable of driving Kozai or other secular variations of outer planets' eccentricities. Dynamical excitation of planets in the outer system reduces the multiplicity of Kepler-detectable planets in the inner system in $\sim20-25\%$ of our systems. Accounting for the occurrence rates of wide-orbit planets and binary stars, $\approx18\%$ of close-in systems could be destabilised by their outer companions in this way. This provides some contribution to the apparent excess of systems with a single transiting planet compared to multiple, however, it only contributes at most $25\%$ of the excess. The effects of the outer dynamics can generate systems similar to Kepler-56 (two coplanar planets significantly misaligned with the host star) and Kepler-108 (two significantly non-coplanar planets in a binary). We also identify three pathways to the formation of eccentric warm Jupiters resulting from the interaction between outer and inner systems: direct inelastic collision between an eccentric outer and an inner planet, secular eccentricity oscillations that may "freeze out" when scattering resolves in the outer system; and scattering in the inner system followed by "uplift", where inner planets are removed by interaction with the outer planets. In these scenarios, the formation of eccentric warm Jupiters is a signature of a past history of violent dynamics among massive planets beyond $\sim1$ au.

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.

Many observed giant planets lie on eccentric orbits. Such orbits could be the result of strong scatterings with other giant planets. The same dynamical instability that produces these scatterings may also cause habitable planets in interior orbits to become ejected, destroyed, or be transported out of the habitable zone. We say that a habitable planet has resilient habitability if it is able to avoid ejections and collisions and its orbit remains inside the habitable zone. Here we model the orbital evolution of rocky planets in planetary systems where giant planets become dynamically unstable. We measure the resilience of habitable planets as a function of the observed, present-day masses and orbits of the giant planets. We find that the survival rate of habitable planets depends strongly on the giant planet architecture. Equal-mass planetary systems are far more destructive than systems with giant planets of unequal masses. We also establish a link with observation; we find that giant planets with present-day eccentricities higher than 0.4 almost never have a habitable interior planet. For a giant planet with an present-day eccentricity of 0.2 and semimajor axis of 5 AU orbiting a Sun-like star, 50% of the orbits in the habitable zone are resilient to the instability. As semimajor axis increases and eccentricity decreases, a higher fraction of habitable planets survive and remain habitable. However, if the habitable planet has rocky siblings, there is a significant risk of rocky planet collisions that would sterilize the planet.

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.

We investigate the prospects for the capture of the proposed Planet 9 from other stars in the Sun's birth cluster. Any capture scenario must satisfy three conditions: the encounter must be more distant than ~150 au to avoid perturbing the Kuiper belt; the other star must have a wide-orbit planet (a>~100au); the planet must be captured onto an appropriate orbit to sculpt the orbital distribution of wide-orbit Solar System bodies. Here we use N-body simulations to show that these criteria may be simultaneously satisfied. In a few percent of slow close encounters in a cluster, bodies are captured onto heliocentric, Planet 9-like orbits. During the ~100 Myr cluster phase, many stars are likely to host planets on highly-eccentric orbits with apastron distances beyond 100 au if Neptune-sized planets are common and susceptible to planet--planet scattering. While the existence of Planet 9 remains unproven, we consider capture from one of the Sun's young brethren a plausible route to explain such an object's orbit. Capture appears to predict a large population of Trans-Neptunian Objects (TNOs) whose orbits are aligned with the captured planet, and we propose that different formation mechanisms will be distinguishable based on their imprint on the distribution of TNOs.

The size distribution of asteroids in the solar system suggests that they formed top-down, with 100-1000 km bodies forming from the gravitational collapse of dense clumps of small solid particles. We investigate the conditions under which solid particles can form dense clumps in a protoplanetary disc. We used a hydrodynamic code to model the solid-gas interaction in disc. We found that particles down to millimeter size can form dense clumps, but only in regions where solids make $\sim$ 8% of the local surface density. More generally, we mapped the range of particle sizes and concentrations that is consistent with the formation of particle clumps.

We explore the role of dynamics in shaping planetary system multiplicities, focussing on two particular problems. (1) We propose that the lack of close-in super-Earths in hot Jupiter systems is a signature of the migration history of the hot Jupiters and helps to discriminate between different mechanisms of migration. We present N-body simulations of dynamical migration scenarios where proto-hot Jupiters are excited to high eccentricities prior to tidal circularisation and orbital decay. We show that in this scenario, the eccentric giant planet typically destroys planets in the inner system, in agreement with the observed lack of close super-Earth companions to hot Jupiters. (2) We explore the role of the dynamics of outer systems in affecting the multiplicities of close-in systems such as those discovered by Kepler. We consider specifically the effects of planet--planet scattering and Kozai perturbations on an exterior giant planet on the architecture of the inner system, and evaluate the ability of such scenarios to reduce the inner system's multiplicity and contribute to the observed excess of single Kepler planets.

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 present deep 15.7-GHz observations made with the Arcminute Microkelvin Imager Large Array in two fields previously observed as part of the Tenth Cambridge (10C) survey. These observations allow the source counts to be calculated down to 0.1 mJy, a factor of five deeper than achieved by the 10C survey. The new source counts are consistent with the extrapolated fit to the 10C source count, and display no evidence for either steepening or flattening of the counts. There is thus no evidence for the emergence of a significant new population of sources (e.g. starforming) at 15.7 GHz flux densities above 0.1 mJy, the flux density level at which we expect starforming galaxies to begin to contribute. Comparisons with the de Zotti et al. model and the SKADS Simulated Sky show that they both underestimate the observed number of sources by a factor of two at this flux density level. We suggest that this is due to the flat-spectrum cores of radio galaxies contributing more significantly to the counts than predicted by the models.

Hot Jupiters are giant planets on orbits a few hundredths of an AU. They do not share their system with low-mass close-in planets, despite these latter being exceedingly common. Two migration channels for hot Jupiters have been proposed: through a protoplanetary gas disc or by tidal circularisation of highly-eccentric planets. We show that highly-eccentric giant planets that will become hot Jupiters clear out any low-mass inner planets in the system, explaining the observed lack of such companions to hot Jupiters. A less common outcome of the interaction is that the giant planet is ejected by the inner planets. Furthermore, the interaction can implant giant planets on moderately-high eccentricities at semimajor axes $<1$ AU, a region otherwise hard to populate. Our work supports the hypothesis that most hot Jupiters reached their current orbits following a phase of high eccentricity, possibly excited by other planetary or stellar companions.

The size distribution of asteroids and Kuiper belt objects in the solar system is difficult to reconcile with a bottom-up formation scenario due to the observed scarcity of objects smaller than $\sim$100 km in size. Instead, planetesimals appear to form top-down, with large $100-1000$ km bodies forming from the rapid gravitational collapse of dense clumps of small solid particles. In this paper we investigate the conditions under which solid particles can form dense clumps in a protoplanetary disk. We use a hydrodynamic code to model the interaction between solid particles and the gas inside a shearing box inside the disk, considering particle sizes from sub-millimeter-sized chondrules to meter-sized rocks. We find that particles down to millimeter sizes can form dense particle clouds through the run-away convergence of radial drift known as the streaming instability. We make a map of the range of conditions (strength of turbulence, particle mass-loading, disk mass, and distance to the star) which are prone to producing dense particle clumps. Finally, we estimate the distribution of collision speeds between mm-sized particles. We calculate the rate of sticking collisions and obtain a robust upper limit on the particle growth timescale of $\sim$$10^5$ years. This means that mm-sized chondrule aggregates can grow on a timescale much smaller than the disk accretion timescale ($\sim$$10^6 - 10^7$ years). Our results suggest a pathway from the mm-sized grains found in primitive meteorites to fully formed asteroids. We speculate that asteroids may form from a positive feedback loop in which coagualation leads to particle clumping driven by the streaming instability. This clumping, in turn reduces collision speeds and enhances coagulation.} Future simulations should model coagulation and the streaming instability together to explore this feedback loop further.

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.

The solar system started to form about 4.56 Gyr ago and despite the long intervening time span, there still exist several clues about its formation. The three major sources for this information are meteorites, the present solar system structure and the planet-forming systems around young stars. In this introduction we give an overview of the current understanding of the solar system formation from all these different research fields. This includes the question of the lifetime of the solar protoplanetary disc, the different stages of planet formation, their duration, and their relative importance. We consider whether meteorite evidence and observations of protoplanetary discs point in the same direction. This will tell us whether our solar system had a typical formation history or an exceptional one. There are also many indications that the solar system formed as part of a star cluster. Here we examine the types of cluster the Sun could have formed in, especially whether its stellar density was at any stage high enough to influence the properties of today's solar system. The likelihood of identifying siblings of the Sun is discussed. Finally, the possible dynamical evolution of the solar system since its formation and its future are considered.

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.

In this chapter we consider two formation channels for blue straggler stars: 1) the merger of two single stars via a collision, and 2) those produced via mass transfer within a binary. We review how computer simulations show that stellar collisions are likely to lead to relatively little mass loss and are thus effective in producing a young population of more-massive stars. The number of blue straggler stars produced by collisions will tend to increase with cluster mass. We review how the current population of blue straggler stars produced from primordial binaries decreases with increasing cluster mass. This is because exchange encounters with third, single stars in the most massive clusters tend to reduce the fraction of binaries containing a primary close to the current turn-off mass. Rather, their primaries tend to be somewhat more massive and have evolved off the main sequence, filling their Roche lobes in the past, often converting their secondaries into blue straggler stars (but more than 1 Gyr or so ago and thus they are no longer visible today as blue straggler stars).

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.

The innermost parsec around Sgr A* has been found to play host to two discs or streamers of O and W-R stars. They are misaligned by an angle approaching 90 degrees. That the stars are approximately coeval indicates that they formed in the same event rather than independently. We have performed SPH simulations of the infall of a single prolate cloud towards a massive black hole. As the cloud is disrupted, the large spread in angular momentum can, if conditions allow, lead to the creation of misaligned gas discs. In turn, stars may form within those discs. We are now investigating the origins of these clouds in the Galactic Centre (GC) region.

This chapter concerns the long-term dynamical evolution of planetary systems from both theoretical and observational perspectives. We begin by discussing the planet-planet interactions that take place within our own Solar System. We then describe such interactions in more tightly-packed planetary systems. As planet-planet interactions build up, some systems become dynamically unstable, leading to strong encounters and ultimately either ejections or collisions of planets. After discussing the basic physical processes involved, we consider how these interactions apply to extrasolar planetary systems and explore the constraints provided by observed systems. The presence of a residual planetesimal disc can lead to planetary migration and hence cause instabilities induced by resonance crossing; however, such discs can also stabilise planetary systems. The crowded birth environment of a planetary system can have a significant impact: close encounters and binary companions can act to destabilise systems, or sculpt their properties. In the case of binaries, the Kozai mechanism can place planets on extremely eccentric orbits which may later circularise to produce hot Jupiters.

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.

We follow the near radial infall of a prolate cloud onto a 4 x 10^6 Msun supermassive black hole in the Galactic Centre using smoothed particle hydrodynamics (SPH). We show that a prolate cloud oriented perpendicular to its orbital plane naturally produces a spread in angular momenta in the gas which can translate into misaligned discs as is seen in the young stars orbiting Sagittarius A*. A turbulent or otherwise highly structured cloud is necessary to avoid cancelling too much angular momentum through shocks at closest approach. Our standard model of a 2 x 10^4 Msun gas cloud brought about the formation of a disc within 0.3 pc from the black hole and a larger, misaligned streamer at 0.5 pc. A total of 1.5 x 10^4 Msun of gas formed these structures. Our exploration of the simulation parameter space showed that when star formation occurred, it resulted in top-heavy IMFs with stars on eccentric orbits with semi-major axes 0.02 to 0.3 pc and inclinations following the gas discs and streamers. We suggest that the single event of an infalling prolate cloud can explain the occurrence of multiple misaligned discs of young stars.

The AMI Galactic Plane Survey (AMIGPS) is a large area survey of the outer Galactic plane to provide arcminute resolution images at milli-Jansky sensitivity in the centimetre-wave band. Here we present the first data release of the survey, consisting of 868 deg^2 of the Galactic plane, covering the area 76 deg \lessapprox l \lessapprox 170 deg between latitudes of |b| \lessapprox 5 deg, at a central frequency of 15.75 GHz (1.9 cm). We describe in detail the drift scan observations which have been used to construct the maps, including the techniques used for observing, mapping and source extraction, and summarise the properties of the finalized datasets. These observations constitute the most sensitive Galactic plane survey of large extent at centimetre-wave frequencies greater than 1.4 GHz.

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.

Massive black holes have been discovered in all closely examined galaxies with high velocity dispersion. The case is not as clear for lower-dispersion systems such as low-mass galaxies and globular clusters. Here we suggest that above a critical velocity dispersion of roughly 40 km/s, massive central black holes will form in relaxed stellar systems at any cosmic epoch. This is because above this dispersion primordial binaries cannot support the system against deep core collapse. If, as previous simulations show, the black holes formed in the cluster settle to produce a dense subcluster, then given the extremely high densities reached during core collapse the holes will merge with each other. For low velocity dispersions and hence low cluster escape speeds, mergers will typically kick out all or all but one of the holes due to three-body kicks or the asymmetric emission of gravitational radiation. If one hole remains, it will tidally disrupt stars at a high rate. If none remain, one is formed after runaway collisions between stars, then it tidally disrupts stars at a high rate. The accretion rate after disruption is many orders of magnitude above Eddington. If, as several studies suggest, the hole can accept matter at that rate because the generated radiation is trapped and advected, then it will grow quickly and form a massive central black hole.

We investigate whether stellar-mass black holes have to receive natal kicks in order to explain the observed distribution of low-mass X-ray binaries containing black holes within our Galaxy. Such binaries are the product of binary evolution, where the massive primary has exploded forming a stellar-mass black hole, probably after a common envelope phase where the system contracted down to separations of order 10-30 Rsun. We perform population synthesis calculations of these binaries, applying both kicks due to supernova mass-loss and natal kicks to the newly-formed black hole. We then integrate the trajectories of the binary systems within the Galactic potential. We find that natal kicks are in fact necessary to reach the large distances above the Galactic plane achieved by some binaries. Further, we find that the distribution of natal kicks would seem to be similar to that of neutron stars, rather than one where the kick velocities are reduced by the ratio of black hole to neutron-star mass (i.e. where the kicks have the same momentum). This result is somewhat surprising; in many pictures of stellar-mass black-hole formation, one might have expected black holes to receive kicks having the same momentum (rather than the same speed) as those given to neutron stars.

We present 16-GHz AMI SZ observations of 19 clusters with L_X >7x10^37 W (h50=1) selected from the LoCuS survey (0.142<z<0.295) and of A1758b, in the FoV of A1758a. We detect 17 clusters with 5-23sigma peak surface brightnesses. Cluster parameters are obtained using a Bayesian cluster analysis. We fit isothermal beta-models to our data and assume the clusters are virialized (with all the kinetic energy in gas internal energy). Our gas temperature, T_AMI, is derived from AMI SZ data, not from X-ray spectroscopy. Cluster parameters internal to r500 are derived assuming HSE. We find: (i) Different gNFW parameterizations yield significantly different parameter degeneracies. (ii) For h70 = 1, we find the virial radius r200 to be typically 1.6+/-0.1 Mpc and the total mass M_T(r200) typically to be 2.0-2.5xM_T(r500).(iii) Where we have found M_T X-ray (X) and weak-lensing (WL) values in the literature, there is good agreement between WL and AMI estimates (with M_{T,AMI}/M_{T,WL} =1.2^{+0.2}_{-0.3} and =1.0+/-0.1 for r500 and r200, respectively). In comparison, most Suzaku/Chandra estimates are higher than for AMI (with M_{T,X}/M_{T,AMI}=1.7+/-0.2 within r500), particularly for the stronger mergers.(iv) Comparison of T_AMI to T_X sheds light on high X-ray masses: even at large r, T_X can substantially exceed T_AMI in mergers. The use of these higher T_X values will give higher X-ray masses. We stress that large-r T_SZ and T_X data are scarce and must be increased. (v) Despite the paucity of data, there is an indication of a relation between merger activity and SZ ellipticity. (vi) At small radius (but away from any cooling flow) the SZ signal (and T_AMI) is less sensitive to ICM disturbance than the X-ray signal (and T_X) and, even at high r, mergers affect n^2-weighted X-ray data more than n-weighted SZ, implying significant shocking or clumping or both occur even in the outer parts of mergers.

We present 16-GHz observations using the Arcminute Microkelvin Imager (AMI) of 11 clusters with 7 x 10^{37}W < L_X < 11 x 10^{37}W (h_{50}=1.0) selected from the Local Cluster Substructure Survey (LoCuSS) and compare them to X-ray data. We use a fast, Bayesian cluster analysis to explore the high-dimensional parameter space of the cluster-plus-sources model and obtain robust cluster parameter estimates in the presence of radio point sources, receiver noise and primordial CMB anisotropy. Our analysis fits a spherical, isothermal beta-model to our data and assumes the cluster follows the theoretical mass-temperature relation. Large-scale cluster parameters internal to r_{500} are derived under the assumption of hydrostatic equilibrium. Posterior distributions for the large-scale parameters of 8 of our clusters are given; SZ effects towards Abell 1704 and Zw0857.9+2107 were not detected and our spherical beta-profile was found to be an inadequate fit to the decrement on our map for Abell 2409.

We present an interesting Sunyaev-Zel'dovich (SZ) detection in the first of the Arcminute Microkelvin Imager (AMI) 'blind', degree-square fields to have been observed down to our target sensitivity of 100{\mu}Jy/beam. In follow-up deep pointed observations the SZ effect is detected with a maximum peak decrement greater than 8 \times the thermal noise. No corresponding emission is visible in the ROSAT all-sky X-ray survey and no cluster is evident in the Palomar all-sky optical survey. Compared with existing SZ images of distant clusters, the extent is large (\approx 10') and complex; our analysis favours a model containing two clusters rather than a single cluster. Our Bayesian analysis is currently limited to modelling each cluster with an ellipsoidal or spherical beta-model, which do not do justice to this decrement. Fitting an ellipsoid to the deeper candidate we find the following. (a) Assuming that the Evrard et al. (2002) approximation to Press & Schechter (1974) correctly gives the number density of clusters as a function of mass and redshift, then, in the search area, the formal Bayesian probability ratio of the AMI detection of this cluster is 7.9 \times 10^4:1; alternatively assuming Jenkins et al. (2001) as the true prior, the formal Bayesian probability ratio of detection is 2.1 \times 10^5:1. (b) The cluster mass is MT,200 = 5.5+1.2\times 10^14h-1M\odot. (c) Abandoning a physical model with num- -1.3 70 ber density prior and instead simply modelling the SZ decrement using a phenomenological {\beta}-model of temperature decrement as a function of angular distance, we find a central SZ temperature decrement of -295+36 {\mu}K - this allows for CMB primary anisotropies, receiver -15 noise and radio sources. We are unsure if the cluster system we observe is a merging system or two separate clusters.

We present 16 GHz (1.9 cm) deep radio continuum observations made with the Arcminute Microkelvin Imager (AMI) of a sample of low-mass young stars driving jets. We combine these new data with archival information from an extensive literature search to examine spectral energy distributions (SEDs) for each source and calculate both the radio and sub-mm spectral indices in two different scenarios: (1) fixing the dust temperature (Td) according to evolutionary class; (2) allowing Td to vary. We use the results of this analysis to place constraints on the physical mechanisms responsible for the radio emission. From AMI data alone, as well as from model fitting to the full SED in both scenarios, we find that 80 per cent of the objects in this sample have spectral indices consistent with free-free emission. We find an average spectral index in both Td scenarios consistent with free-free emission. We examine correlations of the radio luminosity with bolometric luminosity, envelope mass, and outflow force and find that these data are consistent with the strong correlation with envelope mass seen in lower luminosity samples. We examine the errors associated with determining the radio luminosity and find that the dominant source of error is the uncertainty on the opacity index, beta. We examine the SEDs for variability in these young objects, and find evidence for possible radio flare events in the histories of L1551 IRS 5 and Serpens SMM 1.

We present a structural analysis of halo star clusters in M31 based on deep Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) imaging. The clusters in our sample span a range in galactocentric projected distance from 13 to 100 kpc and thus reside in rather remote environments. Ten of the clusters are classical globulars, while four are from the Huxor et al. (2005, 2008) population of extended, old clusters. For most clusters, contamination by M31 halo stars is slight, and so the profiles can be mapped reliably to large radial distances from their centres. We find that the extended clusters are well fit by analytic King (1962) profiles with ~20 parsec core radii and ~100 parsec photometric tidal radii, or by Sersic profiles of index ~1 (i.e. approximately exponential). Most of the classical globulars also have large photometric tidal radii in the range 50-100 parsec, however the King profile is a less good fit in some cases, particularly at small radii. We find 60% of the classical globular clusters exhibit cuspy cores which are reasonably well described by Sersic profiles of index ~2-6. Our analysis also reinforces the finding that luminous classical globulars, with half-light radii <10 parsec, are present out to radii of at least 100 kpc in M31, which is in contrast to the situation in the Milky Way where such clusters (other than the unusual object NGC 2419) are absent beyond 40 kpc.

We present an analysis of observations made with the Arcminute Microkelvin Imager (AMI) and the Canada-France-Hawaii Telescope (CFHT) of six galaxy clusters in a redshift range of 0.16--0.41. The cluster gas is modelled using the Sunyaev--Zel'dovich (SZ) data provided by AMI, while the total mass is modelled using the lensing data from the CFHT. In this paper, we: i) find very good agreement between SZ measurements (assuming large-scale virialisation and a gas-fraction prior) and lensing measurements of the total cluster masses out to r_200; ii) perform the first multiple-component weak-lensing analysis of A115; iii) confirm the unusual separation between the gas and mass components in A1914; iv) jointly analyse the SZ and lensing data for the relaxed cluster A611, confirming our use of a simulation-derived mass-temperature relation for parameterizing measurements of the SZ effect.

Most Sunyaev--Zel'dovich (SZ) and X-ray analyses of galaxy clusters try to constrain the cluster total mass and/or gas mass using parameterised models and assumptions of spherical symmetry and hydrostatic equilibrium. By numerically exploring the probability distributions of the cluster parameters given the simulated interferometric SZ data in the context of Bayesian methods, and assuming a beta-model for the electron number density we investigate the capability of this model and analysis to return the simulated cluster input quantities via three rameterisations. In parameterisation I we assume that the T is an input parameter. We find that parameterisation I can hardly constrain the cluster parameters. We then investigate parameterisations II and III in which fg(r200) replaces temperature as a main variable. In parameterisation II we relate M_T(r200) and T assuming hydrostatic equilibrium. We find that parameterisation II can constrain the cluster physical parameters but the temperature estimate is biased low. In parameterisation III, the virial theorem replaces the hydrostatic equilibrium assumption. We find that parameterisation III results in unbiased estimates of the cluster properties. We generate a second simulated cluster using a generalised NFW (GNFW) pressure profile and analyse it with an entropy based model to take into account the temperature gradient in our analysis and improve the cluster gas density distribution. This model also constrains the cluster physical parameters and the results show a radial decline in the gas temperature as expected. The mean cluster total mass estimates are also within 1 sigma from the simulated cluster true values. However, we find that for at least interferometric SZ analysis in practice at the present time, there is no differences in the AMI visibilities between the two models. This may of course change as the instruments improve.

We present deep 1.8 cm (16 GHz) radio continuum imaging of seven young stellar objects in the Taurus molecular cloud. These objects have previously been extensively studied in the sub-mm to NIR range and their SEDs modelled to provide reliable physical and geometrical parametres.We use this new data to constrain the properties of the long-wavelength tail of the greybody spectrum, which is expected to be dominated by emission from large dust grains in the protostellar disk. We find spectra consistent with the opacity indices expected for such a population, with an average opacity index of beta = 0.26+/-0.22 indicating grain growth within the disks. We use spectra fitted jointly to radio and sub-mm data to separate the contributions from thermal dust and radio emission at 1.8 cm and derive disk masses directly from the cm-wave dust contribution. We find that disk masses derived from these flux densities under assumptions consistent with the literature are systematically higher than those calculated from sub-mm data, and meet the criteria for giant planet formation in a number of cases.

The Planck Early Release Compact Source Catalogue includes 26 sources with no obvious matches in other radio catalogues (of primarily extragalactic sources). Here we present observations made with the Arcminute Microkelvin Imager Small Array (AMI SA) at 15.75 GHz of the eight of the unmatched sources at declination > +10 degrees. Of the eight, four are detected and are associated with known objects. The other four are not detected with the AMI SA, and are thought to be spurious.

The fraction of binary stars is an important ingredient to interpret globular cluster dynamical evolution and their stellar population. We investigate the properties of main-sequence binaries measured in a uniform photometric sample of 59 Galactic globular clusters that were observed by HST WFC/ACS as a part of the Globular Cluster Treasury project. We measured the fraction of binaries and the distribution of mass-ratio as a function of radial location within the cluster, from the central core to beyond the half-mass radius. We studied the radial distribution of binary stars, and the distribution of stellar mass ratios. We investigated monovariate relations between the fraction of binaries and the main parameters of their host clusters. We found that in nearly all the clusters, the total fraction of binaries is significantly smaller than the fraction of binaries in the field, with a few exceptions only. Binary stars are significantly more centrally concentrated than single MS stars in most of the clusters studied in this paper. The distribution of the mass ratio is generally flat (for mass-ratio parameter q>0.5). We found a significant anti-correlation between the binary fraction in a cluster and its absolute luminosity (mass). Some, less significant correlation with the collisional parameter, the central stellar density, and the central velocity dispersion are present. There is no statistically significant relation between the binary fraction and other cluster parameters. We confirm the correlation between the binary fraction and the fraction of blue stragglers in the cluster.

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.

We present deep radio continuum observations of the cores identified as deeply embedded young stellar objects in the Serpens molecular cloud by the Spitzer c2d programme at a wavelength of 1.8cm with the Arcminute Microkelvin Imager Large Array (AMI-LA). These observations have a resolution of ~30arcsec and an average sensitivity of 19microJy/beam. The targets are predominantly Class I sources, and we find the detection rate for Class I objects in this sample to be low (18%) compared to that of Class 0 objects (67%), consistent with previous works. For detected objects we examine correlations of radio luminosity with bolometric luminosity and envelope mass and find that these data support correlations found by previous samples, but do not show any indiction of the evolutionary divide hinted at by similar data from the Perseus molecular cloud when comparing radio luminosity with envelope mass. We conclude that envelope mass provides a better indicator for radio luminosity than bolometric luminosity, based on the distribution of deviations from the two correlations. Combining these new data with archival 3.6cm flux densities we also examine the spectral indices of these objects and find an average spectral index of 0.53+/-1.14, consistent with the canonical value for a partially optically thick spherical or collimated stellar wind. However, we caution that possible inter-epoch variability limits the usefulness of this value, and such variability is supported by our identification of a possible flare in the radio history of Serpens SMM 1.

Black holes exceeding a billion solar masses have been detected at redshifts greater than six. The rapid formation of these objects may suggest a massive early seed or a period of growth faster than Eddington. Here we suggest a new mechanism along these lines. We propose that in the process of hierarchical structure assembly, dense star clusters can be contracted on dynamical time scales due to the nearly free-fall inflow of self-gravitating gas with a mass comparable to or larger than that of the clusters. This increases the velocity dispersion to the point that the few remaining hard binaries can no longer effectively heat the cluster, and the cluster goes into a period of homologous core collapse. The cluster core can then reach a central density high enough for fast mergers of stellar-mass black holes and hence the rapid production of a black hole seed that could be $10^5 M_\odot$ or larger.

A large population of fragile, wide (> 1000 AU) binary systems exists in the Galactic field and halo. These wide binary stars cannot be primordial because of the high stellar density in star forming regions, while formation by capture in the Galactic field is highly improbable. We propose that these binary systems were formed during the dissolution phase of star clusters (see Kouwenhoven et al. 2010, for details). Stars escaping from a dissolving star cluster can have very similar velocities, which can lead to the formation of a wide binary systems. We carry out N-body simulations to test this hypothesis. The results indicate that this mechanism explains the origin of wide binary systems in the Galaxy. The resulting wide binary fraction and semi-major axis distribution depend on the initial conditions of the dissolving star cluster, while the distributions in eccentricity and mass ratio are universal. Finally, since most stars are formed in (relatively tight) primordial binaries, we predict that a large fraction of the wide "binary stars" are in fact higher-order multiple systems.

We present follow-up observations of two galaxy clusters detected blindly via the Sunyaev-Zel'dovich (SZ) effect and released in the Planck Early Release Compact Source Catalogue. We use the Arcminute Microkelvin Imager, a dual-array 14-18 GHz radio interferometer. After radio source subtraction, we find a SZ decrement of integrated flux density -1.08+/-0.10 mJy toward PLCKESZ G121.11+57.01, and improve the position measurement of the cluster, finding the centre to be RA 12 59 36.4, Dec +60 04 46.8, to an accuracy of 20 arcseconds. The region of PLCKESZ G115.71+17.52 contains strong extended emission, so we are unable to confirm the presence of this cluster via the SZ effect.

We present deep radio continuum observations of the cores identified as deeply embedded young stellar objects in the Perseus molecular cloud by the Spitzer c2d programme at a wavelength of 1.8 cm with the Arcminute Microkelvin Imager Large Array (AMI-LA). We detect 72% of Class 0 objects from this sample and 31% of Class I objects. No starless cores are detected. We use the flux densities measured from these data to improve constraints on the correlations between radio luminosity and bolometric luminosity, infrared luminosity and, where measured, outflow force. We discuss the differing behaviour of these objects as a function of protostellar class and investigate the differences in radio emission as a function of core mass. Two of four possible very low luminosity objects (VeLLOs) are detected at 1.8 cm.

We have observed an area of approximatley 27 deg^2 to an rms noise level of less than 0.2 mJy at 15.7 GHz, using the Arcminute Microkelvin Imager Large Array. These observations constitute the most sensitive radio-source survey of any extent (greater than approximately 0.2 deg^2) above 1.4 GHz. This paper presents the techniques employed for observing, mapping and source extraction. We have used a systematic procedure for extracting information and producing source catalogues, from maps with varying noise and uv-coverage. We have performed simulations to test our mapping and source-extraction procedures, and developed methods for identifying extended, overlapping and spurious sources in noisy images. In an accompanying paper, AMI Consortium: Davies et al. 2010, the first results from the 10C survey, including the deep 15.7-GHz source count, are presented.

The first results from the Tenth Cambridge (10C) Survey of Radio Sources, carried out using the AMI Large Array (LA) at an observing frequency of 15.7 GHz, are presented. The survey fields cover an area of approximately 27 sq. degrees to a flux-density completeness of 1 mJy. Results for some deeper areas, covering approximately 12 sq. degrees, wholly contained within the total areas and complete to 0.5 mJy, are also presented. The completeness for both areas is estimated to be at least 93 per cent. The source catalogue contains 1897 entries and is available at www.mrao.cam.ac.uk/surveys/10C. It has been combined with that of the 9C Survey to calculate the 15.7-GHz source counts. A broken power law is found to provide a good parameterisation of the differential count between 0.5 mJy and 1 Jy. The measured count has been compared to that predicted by de Zotti et al. (2005). The model displays good agreement with the data at the highest flux densities but under-predicts the integrated count between 0.5 mJy and 1 Jy by about 30 per cent. Entries from the source catalogue have been matched to those contained in the catalogues of NVSS and FIRST (both of which have observing frequencies of 1.4 GHz). This matching provides evidence for a shift in the typical 1.4-to-15.7-GHz spectral index of the 15.7-GHz-selected source population with decreasing flux density towards sub-mJy levels - the spectra tend to become less steep. Automated methods for detecting extended sources have been applied to the data; approximately 5 per cent of the sources are found to be extended relative to the LA synthesised beam of approximately 30 arcsec. Investigations using higher-resolution data showed that most of the genuinely extended sources at 16 GHz are classical doubles, although some nearby galaxies and twin-jet sources were also identified.

We present 13.9-18.2 GHz observations of the Sunyaev-Zel'dovich (SZ) effect towards Abell 2146 using the Arcminute Microkelvin Imager (AMI). The cluster is detected with a peak SNR ratio of 13 sigma in the radio source subtracted map. Comparison of the SZ and X-ray images suggests that they both have extended regions which lie approximately perpendicular to one another, with their emission peaks significantly displaced. These features indicate non-uniformities in the distributions of the gas temperature and pressure, indicative of a cluster merger. We use a Bayesian cluster analysis to explore the high-dimensional parameter space of the cluster-plus-sources model to obtain cluster parameter estimates in the presence of radio point sources, receiver noise and primordial CMB anisotropy; the probability of SZ + CMB primordial structure + radio sources + receiver noise to CMB + radio sources + receiver noise is 3 x 10^{6}:1. We compare the results from three different cluster models. Our preferred model exploits the observation that the gas fractions do not appear to vary greatly between clusters. Given the relative masses of the two merging systems in Abell 2146, the mean gas temperature can be deduced from the virial theorem (assuming all of the kinetic energy is in the form of internal gas energy) without being affected significantly by the merger event, provided the primary cluster was virialized before the merger. In this model we fit a simple spherical isothermal beta-model, despite the inadequacy of this model for a merging system like Abell 2146, and assume the cluster follows the mass-temperature relation of a virialized, singular, isothermal sphere. We note that this model avoids inferring large-scale cluster parameters internal to r_200 under the widely used assumption of hydrostatic equilibrium. We find that at r_200 M_T= 4.1 \pm 0.5 x 10^{14} h^{-1}M_sun and T=4.5 \pm 0.5 keV.

We present observations from the Small Array of the Arcminute Microkelvin Imager (AMI) of eight high X-ray luminosity galaxy cluster systems selected from the Local Cluster Substructure Survey (LoCuSS) sample.We detect the Sunyaev-Zel'dovich (SZ) effect in seven of these clusters. With the assumptions that galaxy clusters are isothermal, have a density profile described by a spherical b -model and obey the theoretical M-T relation, we are able to derive cluster parameters at r200 from our SZ data. With the additional assumption of hydrostatic equilibrium we are able to derive parameters at r500. We present posterior probability distributions for cluster parameters such as mass, radius and temperature (TSZ, MT). Combining our sample with that of AMI Consortium: Rodr'\iguez-Gonz\'alvez et al. (2011) and using large-radius X-ray temperature estimates (TX) from Chandra and Suzaku observations, we find that there is reasonable correspondence between TX and TSZ,MT values at low TX, but that for clusters with TX above around 6keV the correspondence breaks down with TX exceeding TSZ, MT; we stress that this finding is based on just ten clusters.

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 summarize results of a search for X-ray-emitting binary stars in the massive globular cluster Omega Centauri (NGC 5139) using Chandra and HST. ACIS-I imaging reveals 180 X-ray sources, of which we estimate that 45-70 are associated with the cluster. We present 40 identifications, most of which we have obtained using ACS/WFC imaging with HST that covers the central 10'x10' of the cluster. Roughly half of the optical IDs are accreting binary stars, including 9 very faint blue stars that we suggest are cataclysmic variables near the period limit. Another quarter comprise a variety of different systems all likely to contain coronally active stars. The remaining 9 X-ray-bright stars are an intriguing group that appears redward of the red giant branch, with several lying along the anomalous RGB. Future spectroscopic observations should reveal whether these stars are in fact related to the anomalous RGB, or whether they instead represent a large group of "sub-subgiants" such as have been seen in smaller numbers in other globular and open clusters.

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.

Most of the observed extrasolar planets are found on tight and often eccentric orbits. The high eccentricities are not easily explained by planet-formation models, which predict that planets should be on rather circular orbits. Here we explore whether fly-bys involving planetary systems with properties similar to those of the gas giants in the solar system, can produce planets with properties similar to the observed planets. Using numerical simulations, we show that fly-bys can cause the immediate ejection of planets, and sometimes also lead to the capture of one or more planets by the intruder. More common, however, is that fly-bys only perturb the orbits of planets, sometimes leaving the system in an unstable state. Over time-scales of a few million to several hundred million years after the fly-by, this perturbation can trigger planet-planet scatterings, leading to the ejection of one or more planets. For example, in the case of the four gas giants of the solar system, the fraction of systems from which at least one planet is ejected more than doubles in 10^8 years after the fly-by. The remaining planets are often left on more eccentric orbits, similar to the eccentricities of the observed extrasolar planets. We combine our results of how fly-bys effect solar-system-like planetary systems, with the rate at which encounters in young stellar clusters occur. For example, we measure the effects of fly-bys on the four gas giants in the solar system. We find, that for such systems, between 5 and 15 per cent suffer ejections of planets in 10^8 years after fly-bys in typical open clusters. Thus, encounters in young stellar clusters can significantly alter the properties of any planets orbiting stars in clusters. As a large fraction of stars which populate the solar neighbourhood form in stellar clusters, encounters can significantly affect the properties of the observed extrasolar planets.

We perform deep 1.8 cm radio continuum imaging towards thirteen protostellar regions selected from the Spitzer c2d small clouds and cores programme at high resolution (25") in order to detect and quantify the cm-wave emission from deeply embedded young protostars. Within these regions we detect fifteen compact radio sources which we identify as radio protostars including two probable new detections. The sample is in general of low bolometric luminosity and contains several of the newly detected VeLLO sources. We determine the 1.8 cm radio luminosity to bolometric luminosity correlation, L_rad -L_bol, for the sample and discuss the nature of the radio emission in terms of the available sources of ionized gas. We also investigate the L_rad-L_IR correlation and suggest that radio flux density may be used as a proxy for the internal luminosity of low luminosity protostars.

We present observations using the Small Array of the Arcminute Microkelvin Imager (AMI; 14-18 GHz) of four Abell and three MACS clusters spanning 0.171-0.686 in redshift. We detect Sunyaev-Zel'dovich (SZ) signals in five of these without any attempt at source subtraction, although strong source contamination is present. With radio-source measurements from high-resolution observations, and under the assumptions of spherical $\beta$-model, isothermality and hydrostatic equilibrium, a Bayesian analysis of the data in the visibility plane detects extended SZ decrements in all seven clusters over and above receiver noise, radio sources and primary CMB imprints. Bayesian evidence ratios range from 10^{11}:1 to 10^{43}:1 for six of the clusters and 3000:1 for one with substantially less data than the others. We present posterior probability distributions for, e.g., total mass and gas fraction averaged over radii internal to which the mean overdensity is 1000, 500 and 200, r_200 being the virial radius. Reaching r_200 involves some extrapolation for the nearer clusters but not for the more-distant ones. We find that our estimates of gas fraction are low (compared with most in the literature) and decrease with increasing radius. These results appear to be consistent with the notion that gas temperature in fact falls with distance (away from near the cluster centre) out to the virial radius.

We report new cm-wave measurements at five frequencies between 15 and 18GHz of the continuum emission from the reportedly anomalous "region 4" of the nearby galaxy NGC6946. We find that the emission in this frequency range is significantly in excess of that measured at 8.5GHz, but has a spectrum from 15-18GHz consistent with optically thin free-free emission from a compact HII region. In combination with previously published data we fit four emission models containing different continuum components using the Bayesian spectrum analysis package radiospec. These fits show that, in combination with data at other frequencies, a model with a spinning dust component is slightly preferred to those that possess better-established emission mechanisms.

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.

Over the past few decades, numerous wide (>1000 au) binaries in the Galactic field and halo have been discovered. Their existence cannot be explained by the process of star formation or by dynamical interactions in the field, and their origin has long been a mystery. We explain the origin of these wide binaries by formation during the dissolution phase of young star clusters: an initially unbound pair of stars may form a binary when their distance in phase-space is small. Using N-body simulations, we find that the resulting wide binary fraction in the semi-major axis range 1000 au - 0.1 pc for individual clusters is 1-30%, depending on the initial conditions. The existence of numerous wide binaries in the field is consistent with observational evidence that most clusters start out with a large degree of substructure. The wide binary fraction decreases strongly with increasing cluster mass, and the semi-major axis of the newly formed binaries is determined by the initial cluster size. The resulting eccentricity distribution is thermal, and the mass ratio distribution is consistent with gravitationally-focused random pairing. As a large fraction of the stars form in primordial binaries, we predict that a large number of the observed 'wide binaries' are in fact triple or quadruple systems. By integrating over the initial cluster mass distribution, we predict a binary fraction of a few per cent in the semi-major axis range 1000 au - 0.1 pc in the Galactic field, which is smaller than the observed wide binary fraction. However, this discrepancy may be solved when we consider a broad range of cluster morphologies.

We present 25 arcsecond resolution radio images of five Lynds Dark Nebulae (L675, L944, L1103, L1111 & L1246) at 16 GHz made with the Arcminute Microkelvin Imager (AMI) Large Array. These objects were previously observed with the AMI Small Array to have an excess of emission at microwave frequencies relative to lower frequency radio data. In L675 we find a flat spectrum compact radio counterpart to the 850 micron emission seen with SCUBA and suggest that it is cm-wave emission from a previously unknown deeply embedded young protostar. In the case of L1246 the cm-wave emission is spatially correlated with 8 micron emission seen with Spitzer. Since the MIR emission is present only in Spitzer band 4 we suggest that it arises from a population of PAH molecules, which also give rise to the cm-wave emission through spinning dust emission.

The majority of stars in the Galactic field and halo are part of binary or multiple systems. A significant fraction of these systems have orbital separations in excess of thousands of astronomical units, and systems wider than a parsec have been identified in the Galactic halo. These binary systems cannot have formed through the 'normal' star-formation process, nor by capture processes in the Galactic field. We propose that these wide systems were formed during the dissolution phase of young star clusters. We test this hypothesis using N-body simulations of evolving star clusters and find wide binary fractions of 1-30%, depending on initial conditions. Moreover, given that most stars form as part of a binary system, our theory predicts that a large fraction of the known wide 'binaries' are, in fact, multiple systems.

Small angular scale (high l) studies of cosmic microwave background anisotropies require accurate knowledge of the statistical properties of extragalactic sources at cm-mm wavelengths. We have used a 30 GHz dual-beam receiver (OCRA-p) on the Torun 32-m telescope to measure the flux densities of 121 sources in VSA fields selected at 15 GHz with the Ryle Telescope. We have detected 57 sources above a limiting flux density of 5mJy, of which 31 sources have a flux density greater than 10mJy, which is our effective completeness limit. From these measurements we derive a surface density of sources above 10mJy at 30 GHz of 2.0+/-0.4 per square degree. This is consistent with the surface density obtained by Mason et al. (2009) who observed a large sample of sources selected at a much lower frequency (1.4 GHz). We have also investigated the dependence of the spectral index distribution on flux density by comparing our results with those for sources above 1 Jy selected from the WMAP 22 GHz catalogue. We conclude that the proportion of steep spectrum sources increases with decreasing flux density, qualitatively consistent with the predictions of deZotti et al. (2005). We find no evidence for an unexpected population of sources whose spectra rise towards high frequencies, which would affect our ability to interpret current high resolution CMB observations at 30 GHz and above.

Over the last decades, numerous wide (>1000 AU) binaries have been discovered in the Galactic field and halo. The origin of these wide binaries cannot be explained by star formation or by dynamical interactions in the Galactic field. We explain their existence by wide binary formation during the dissolution phase of young star clusters. In this scenario, two single stars that leave the dissolving cluster at the same time, in the same direction, and with similar velocities, form a new, very wide binary. Using N-body simulations we study how frequently this occurs, and how the orbital parameters of such binaries depend on the properties of the cluster from which they originate. The resulting wide binary fraction for individual star clusters is 1-30%, depending on the initial conditions. As most stars form as part of a binary or multiple system, we predict that a large fraction of these wide binaries are in fact wide triple and quadruple systems.

Observations at 14.2 to 17.9 GHz made with the AMI Small Array towards fourteen Lynds Dark Nebulae with a resolution of 2' are reported. These sources are selected from the SCUBA observations of Visser et al. (2001) as small angular diameter clouds well matched to the synthesized beam of the AMI Small Array. Comparison of the AMI observations with radio observations at lower frequencies with matched uv-plane coverage is made, in order to search for any anomalous excess emission which can be attributed to spinning dust. Possible emission from spinning dust is identified as a source within a 2' radius of the Scuba position of the Lynds dark nebula, exhibiting an excess with respect to lower frequency radio emission. We find five sources which show a possible spinning dust component in their spectra. These sources have rising spectral indices in the frequency range 14.2--17.9 GHz. Of these five one has already been reported, L1111, we report one new definite detection, L675, and three new probable detections (L944, L1103 and L1246). The relative certainty of these detections is assessed on the basis of three criteria: the extent of the emission, the coincidence of the emission with the Scuba position and the likelihood of alternative explanations for the excess. Extended microwave emission makes the likelihood of the anomalous emission arising as a consequence of a radio counterpart to a protostar or a proto-planetary disk unlikely. We use a 2' radius in order to be consistent with the IRAS identifications of dark nebulae (Parker 1988), and our third criterion is used in the case of L1103 where a high flux density at 850 microns relative to the FIR data suggests a more complicated emission spectrum.

In our first paper we described three regions of the 9C survey of radio sources with the Ryle telescope at 15.2 GHz, constituting a total area of 520 deg^2 to a completeness limit of approximately 25 mJy. Here we report on a series of deeper regions, amounting to an area of 115 deg^2 complete to approximately 10 mJy and of 29 deg^2 complete to approximately 5.5 mJy. We have investigated the source counts and the distributions of the 1.4 to 15.2 GHz spectral index (alpha) for these deeper samples. The whole catalogue of 643 sources is available online. Down to our lower limit of 5.5 mJy we detect no evidence for any change in the differential source count from the earlier fitted count above 25 mJy. We have matched both our new and earlier catalogues with the NRAO VLA Sky Survey (NVSS) catalogue at 1.4 GHz. For samples of sources selected at 15.2 GHz, in three flux density ranges, we detect a significant shift in the median value of alpha; samples with higher flux densities have higher proportions of sources with flat and rising spectra. We suggest that this observed shift is consistent with a model containing two distinct source populations having differently sloped source counts. Samples selected at 1.4 GHz contain significantly smaller proportions of sources with flat and rising spectra. Also, in our area complete to approximately 10 mJy, we find 5 sources between 10 to 15 mJy, amounting to 4.3 per cent of sources in this range, with no counterpart in the NVSS catalogue. These results illustrate the problems inherent in using a low frequency catalogue to characterise the source population at a much higher frequency and emphasise the value of our blind 15.2 GHz survey.

Using the Arcminute Microkelvin Imager (AMI) at 16 GHz and the Very Small Array (VSA) at 33 GHz to make follow-up observations of sources in the New Extragalactic WMAP Point Source Catalogue, we have investigated the flux-density variability in a complete sample of 97 sources over timescales of a few months to approximately 1.5 years. We find that 53 per cent of the 93 sources, for which we have multiple observations, are variable, at the 99 per cent confidence level, above the flux density calibration uncertainties of approximately 4 per cent at 16 GHz; the fraction of sources having varied by more than 20 per cent is 15 per cent at 16 GHz and 20 per cent at 33 GHz. Not only is this common occurrence of variability at high frequency of interest for source physics, but strategies for coping with source contamination in CMB work must take this variability into account. There is no strong evidence of a correlation between variability and flux density for the sample as a whole. Using a maximum-likelihood method, we calculate the correlation in the variability at the two frequencies in a subset of sources classified as variable from both the AMI and VSA data and find the Pearson product-moment correlation coefficient to be very high (0.955 +/- 0.034). We also find the degree of variability at 16 GHz (0.202 +/- 0.028) to be very similar to that at 33 GHz (0.224 +/- 0.039). Finally, we have investigated the relationship between variability and spectral index, alpha_13.9^33.75 (where S proportional to nu^-alpha), and find a significant difference in the spectral indices of the variable sources (-0.06 +/- 0.05) and non-variable sources (0.13 +/- 0.04).

We present follow-up observations of 97 point sources from the Wilkinson Microwave Anisotropy Probe (WMAP) 3-year data, contained within the New Extragalactic WMAP Point Source (NEWPS) catalogue between declinations of -4 and +60 degrees; the sources form a flux-density-limited sample complete to 1.1 Jy (approximately 5 sigma) at 33 GHz. Our observations were made at 16 GHz using the Arcminute Microkelvin Imager (AMI) and at 33 GHz with the Very Small Array (VSA). 94 of the sources have reliable, simultaneous -- typically a few minutes apart -- observations with both telescopes. The spectra between 13.9 and 33.75 GHz are very different from those of bright sources at low frequency: 44 per cent have rising spectra (alpha < 0.0), where flux density is proportional to frequency^-alpha, and 93 per cent have spectra with alpha < 0.5; the median spectral index is 0.04. For the brighter sources, the agreement between VSA and WMAP 33-GHz flux densities averaged over sources is very good. However, for the fainter sources, the VSA tends to measure lower values for the flux densities than WMAP. We suggest that the main cause of this effect is Eddington bias arising from variability.

We report on the discovery of three new pulsars in the first blind survey of the north Galactic plane (45 degrees < l < 135 degrees; |b| < 1 degrees with the Giant Meterwave Radio telescope (GMRT) at an intermediate frequency of 610 MHz. The survey covered 106 square degrees with a sensitivity of roughly 1 mJy to long-period pulsars (pulsars with period longer than 1 s). The three new pulsars have periods of 318, 933, and 1056 ms. Their timing parameters and flux densities, obtained in follow up observations with the Lovell Telescope at Jodrell Bank and the GMRT, are presented. We also report on pulse nulling behaviour in one of the newly discovered pulsars, PSR J2208+5500.

We present observations between 1.4 and 18 GHz confirming that G64.5+0.9 is new Galactic shell supernova remnant, using the Very Large Array and the Arcminute Microkelvin Imager. The remnant is a shell ~8 arcmin in diameter with a spectral index of alpha = 0.47 +/- 0.03 (with alpha defined such that flux density S varies with frequency nu as S proportional to nu to the power of -alpha). There is also emission near the centre of the shell, ~1 arcmin in extent, with a spectral index of alpha = 0.81 +/- 0.02. We do not find any evidence for spectral breaks for either source within our frequency range. The nature of the central object is unclear and requires further investigation, but we argue that is most unlikely to be extragalactic. It is difficult to avoid the conclusion that it is associated with the shell, although its spectrum is very unlike that of known pulsar wind nebulae.

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.

We analyze a ~70 ksec Chandra ACIS-I exposure of the globular cluster Omega Centauri (NGC 5139). The ~17 amin x 17 amin field of view fully encompasses three core radii and almost twice the half-mass radius. We detect 180 sources to a limiting flux of ~4.3x10^-16 erg/cm^2/s (Lx = 1.2x10^30 erg/s at 4.9 kpc). After accounting for the number of active galactic nuclei and possible foreground stars, we estimate that 45-70 of the sources are cluster members. Four of the X-ray sources have previously been identified as compact accreting binaries in the cluster--three cataclysmic variables (CVs) and one quiescent neutron star. Correlating the Chandra positions with known variable stars yields eight matches, of which five are probable cluster members that are likely to be binary stars with active coronae. Extrapolating these optical identifications to the remaining unidentified X-ray source population, we estimate that 20-35 of the sources are CVs and a similar number are active binaries. This likely represents most of the CVs in the cluster, but only a small fraction of all the active binaries. We place a 2-sigma upper limit of Lx < 3x10^30 erg/s on the integrated luminosity of any additional faint, unresolved population of sources in the core. We explore the significance of these findings in the context of primordial vs. dynamical channels for CV formation. The number of CVs per unit mass in Omega Cen is at least 2-3 times lower than in the field, suggesting that primordial binaries that would otherwise lead to CVs are being destroyed in the cluster environment.

We present observations between 14.2 and 17.9 GHz of 12 reported supernova remnants (SNRs) made with the Arcminute Microkelvin Imager Small Array (AMI SA). In conjunction with data from the literature at lower radio frequencies, we determine spectra of these objects. For well-studied SNRs (Cas A, Tycho's SNR, 3C58 and the Crab Nebula), the results are in good agreement with spectra based on previous results. For the less well-studied remnants the AMI SA observations provide higher-frequency radio observations than previously available, and better constrain their radio spectra. The AMI SA results confirm a spectral turnover at ~11 GHz for the filled-centre remnant G74.9+1.2. We also see a possible steepening of the spectrum of the filled-centre remnant G54.1+0.3 within the AMI SA frequency band compared with lower frequencies. We confirm that G84.9+0.5, which had previously been identified as a SNR, is rather an HII region and has a flat radio spectrum.

We present observations of the Lynds' dark nebula LDN 1111 made at microwave frequencies between 14.6 and 17.2 GHz with the Arcminute Microkelvin Imager (AMI). We find emission in this frequency band in excess of a thermal free--free spectrum extrapolated from data at 1.4 GHz with matched uv-coverage. This excess is > 15 sigma above the predicted emission. We fit the measured spectrum using the spinning dust model of Drain & Lazarian (1998a) and find the best fitting model parameters agree well with those derived from Scuba data for this object by Visser et al. (2001).

Most observed extrasolar planets have masses similar to, but orbits very different from, the gas giants of our solar system. Many are much closer to their parent stars than would have been expected and their orbits are often rather eccentric. We show that some of these planets might have formed in systems much like our solar system, i.e. in systems where the gas giants were originally on orbits with a semi-major axis of several au, but where the masses of the gas giants were all rather similar. If such a system is perturbed by another star, strong planet-planet interactions follow, causing the ejection of several planets while leaving those remaining on much tighter and more eccentric orbits. The eccentricity distribution of these perturbed systems is very similar to that of the observed extrasolar planets with semi-major axis between 1 and 6 au.

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.

The Arcminute Microkelvin Imager is a pair of interferometer arrays operating with six frequency channels spanning 13.9-18.2 GHz, with very high sensitivity to angular scales 30''-10'. The telescope is aimed principally at Sunyaev-Zel'dovich imaging of clusters of galaxies. We discuss the design of the telescope and describe and explain its electronic and mechanical systems.

We present the results from realistic N-body modelling of massive star clusters in the Magellanic Clouds. We have computed eight simulations with N ~ 10^5 particles; six of these were evolved for at least a Hubble time. The aim of this modelling is to examine the possibility of large-scale core expansion in massive star clusters and search for a viable dynamical origin for the radius-age trend observed for such objects in the Magellanic Clouds. We identify two physical processes which can lead to significant and prolonged cluster core expansion: mass-loss due to rapid stellar evolution in a primordially mass segregated cluster, and heating due to a retained population of stellar-mass black holes. These two processes operate over different time-scales - the former occurs only at early times and cannot drive core expansion for longer than a few hundred Myr, while the latter typically does not begin until several hundred Myr have passed but can result in core expansion lasting for many Gyr. We investigate the behaviour of these expansion mechanisms in clusters with varying degrees of primordial mass segregation and in clusters with varying black hole retention fractions. In combination, the two processes can lead to a wide variety of evolutionary paths on the radius-age plane, which fully cover the observed cluster distribution and hence define a dynamical origin for the radius-age trend in the Magellanic Clouds. We discuss the implications of core expansion for various aspects of globular cluster research, as well as the possibility of observationally inferring the presence of a population of stellar-mass black holes in a cluster.

We present observations between 14.2 and 17.9 GHz of sixteen Galactic HII regions made with the Arcminute Microkelvin Imager (AMI). In conjunction with data from the literature at lower radio frequencies we investigate the possibility of a spinning dust component in the spectra of these objects. We conclude that there is no significant evidence for spinning dust towards these sources and measure an average spectral index of 0.15+/-0.07 between 1.4 and 17.9 GHz for the sample.

We report on the discovery of three new pulsars in the first blind survey of the north Galactic plane (45 < l < 135 ; |b| < 1) with the Giant Meterwave Radio telescope (GMRT) at an intermediate frequency of 610 MHz. The timing parameters, obtained in follow up observations with the Lovell Telescope at Jodrell Bank Observatory and the GMRT, are presented.

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.

We present the results of realistic N-body modelling of massive star clusters in the Magellanic Clouds, aimed at investigating a dynamical origin for the radius-age trend observed in these systems. We find that stellar-mass black holes, formed in the supernova explosions of the most massive cluster stars, can constitute a dynamically important population. If a significant number of black holes are retained (here we assume complete retention), these objects rapidly form a dense core where interactions are common, resulting in the scattering of black holes into the cluster halo, and the ejection of black holes from the cluster. These two processes heat the stellar component, resulting in prolonged core expansion of a magnitude matching the observations. Significant core evolution is also observed in Magellanic Cloud clusters at early times. We find that this does not result from the action of black holes, but can be reproduced by the effects of mass-loss due to rapid stellar evolution in a primordially mass segregated cluster.

Pinpointing the progenitors of long duration gamma-ray bursts (LGRBs) remains an extremely important question, although it is now clear that at least a fraction of LGRBs originate in the core collapse of massive stars in type Ic supernovae, the pathways to the production of these stars, and their initial masses, remain uncertain. Rotation is thought to be vital in the creation of LGRBs, and it is likely that black hole creation is also necessary. We suggest that these two constraints can be met if the GRB progenitors are very massive stars (>20 solar masses) and are formed in tight binary systems. Using simple models we compare the predictions of this scenario with observations and find that the location of GRBs on their host galaxies are suggestive of main-sequence masses in excess of 20 solar masses, while 50% of the known compact binary systems may have been sufficiently close to have had the necessary rotation rates for GRB creation. Thus, massive stars in compact binaries are a likely channel for at least some fraction of LGRBs.

Recent comparative observations of long duration gamma-ray bursts (LGRBs) and core collapse supernovae (cc SN) host galaxies demonstrate that these two, highly energetic transient events are distributed very differently upon their hosts. LGRBs are much more concentrated on their host galaxy light than cc SN. Here we explore the suggestion that this differing distribution reflects different progenitor masses for LGRBs and cc SN. Using a simple model we show that, assuming cc SN arise from stars with main sequence masses $>$8 M$_{\odot}$, GRBs are likely to arise from stars with initial masses $>$ 20 M$_{\odot}$. This difference can naturally be explained by the requirement that stars which create a LGRB must also create a black hole.

In this letter we consider the evolution of a planetary system around a star inside a wide binary. We simulate numerically the evolution of the planetary orbits for both co-planar and highly-inclined systems. We find that the Kozai mechanism operates in the latter case. This produces a highly eccentric outer planet whose orbit crosses those of some of the inner planets. Strong planet-planet interactions then follow resulting in the ejection of one or more planets. We note that planetary systems resembling our solar system, formed around single stars in stellar clusters may exchange into binaries and thus will be vulnerable to planet stripping. This process will reduce the number of solar-system like planetary systems, and may produce at least some of the observed extra-solar planets.

We consider gamma--ray bursts produced by the merger of a massive white dwarf with a neutron star. We show that these are likely to produce long--duration GRBs, in some cases definitely without an accompanying supernova, as observed recently. This class of burst would have a strong correlation with star formation, and occur close to the host galaxy. However rare members of the class need not be near star--forming regions, and could have any type of host galaxy. Thus a long--duration burst far from any star--forming region would also be a signature of this class. Estimates based on the existence of a known progenitor suggest that our proposed class may be an important contributor to the observed GRB rate.

Recent progress on the nature of short duration Gamma Ray Bursts (GRBs) has shown that a fraction of them originate in the local universe. These systems may well be the result of giant flares from Soft Gamma Repeaters (SGRs) believed to be magnetars (neutron stars with extremely large magnetic fields >= 10^14 G). If these magnetars are formed via the core collapse of massive stars, then it would be expected that the bursts should originate from predominantly young stellar populations. However, correlating the positions of BATSE short bursts with structure in the local universe reveals a correlation with all galaxy types, including those with little or no ongoing star formation. This is a natural outcome if, in addition to magnetars forming via the core collapse of massive stars, they also form via Accretion Induced Collapse following the merger of two white dwarfs, one of which is magnetic. We investigate this possibility and find that the rate of magnetar production via WD-WD mergers in the Milky Way is comparable to the rate of production via core collapse. However, while the rate of magnetar production by core collapse is proportional to the star formation rate, the rate of production via WD-WD mergers (which have long lifetimes) is proportional to the stellar mass density, which is concentrated in early-type systems. Therefore magnetars produced via WD-WD mergers may produce SGR giant flares which can be identified with early type galaxies. We also comment on the possibility that this mechanism could produce a fraction of the observed short duration GRB population at low redshift.

We present a review of some of the current major challenges in stellar cluster research, including young clusters, globular clusters, and galactic nuclei. Topics considered include: primordial mass segregation and runaway mergers, expulsion of gas from clusters, the production of stellar exotica seen in some clusters (eg blue stragglers and extreme horizontal--branch stars), binary populations within clusters, the black--hole population within stellar clusters, the final parsec problem, stellar dynamics around a massive black hole, and stellar collisions. The Modest Questions posed here are the outcome of discussions which took place at the Modest-6A workshop held in Lund, Sweden, in December, 2005. Modest-6A was organised as part of the activities of the Modest Collaboration (see www.manybody.org for further details)

Aims: This paper presents a new homogeneous catalogue of blue straggler stars (BSS) in Galactic open clusters. Methods: Photometric data for 216 clusters were collected from the literature and 2782 BSS candidates were extracted from 76 of them. Results: We found that the anticorrelation of BSS frequency vs. total magnitude identified in similar studies conducted on Galactic globular clusters extends to the open cluster regime: clusters with smaller total magnitude tend to have higher BSS frequencies. Moreover, a clear correlation between the BSS frequency (obtained normalising the total number of BSS either to the total cluster mass or, for the older clusters, to the total number of clump stars) and the age of the clusters was found. A simple model is developed here to try to explain this last and new result. The model allows us to ascertain the important effect played by mass loss in the evolution of open clusters.

Cosmological long-duration gamma-ray bursts (LGRBs) are thought to originate from the core collapse to black holes of stripped massive stars. Those with sufficient rotation form a centrifugally-supported torus whose collapse powers the GRB. We investigate the role of tidal locking within a tight binary as a source of the necessary angular momentum. We find that the binary orbit must be no wider than a few solar radii for a torus to form upon core collapse. Comparing this criterion to the observed population of binaries containing two compact objects suggests that rotation may have been important in the formation of up to 50% of the observed systems. As these systems created a neutron star and not a black hole they presumably did not produce highly luminous GRBs. We suggest instead that they make the subset of GRBs in the relatively local universe which have much lower luminosity.

There has been recent speculation (Davies & King 2005) that the cores of intermediate-mass stars stripped of their envelopes by tidal interaction with the supermassive black hole in the Galactic centre could form a population observationally similar to the so-called Sgr A* cluster or `S' stars, which have close eccentric orbits around the hole. We model the evolution of such stars, and show that the more luminous end of the population may indeed appear similar to young B stars within the observational limits of the Galactic Centre region. Whether some or all of these cluster stars can be accounted for in this manner depends strongly on the assumed IMF of the loss cone stars and the scattering rate. If most of the observed stars are in fact scattered from the Galactic Centre inner cusp region itself then the population of ~20 to current observational limits may be reproduced. However, this only works if the local relaxation time is small and relies on the cusp stars themselves being young, i.e. it is dependent on some star formation being possible in the central few parsecs. Conversely, we obtain a possible constraint on the tidal stripping rate of `normal'-IMF stars if there are not to be red stars visible in the Sgr A* cluster.

We study close encounters involving massive main sequence stars and the evolution of the exotic products of these encounters as common--envelope systems or possible hypernova progenitors. We show that parabolic encounters between low-- and high--mass stars and between two high--mass stars with small periastrons result in mergers on timescales of a few tens of stellar freefall times (a few tens of hours). We show that such mergers of unevolved low--mass stars with evolved high--mass stars result in little mass loss ($\sim0.01$ M$_{\odot}$) and can deliver sufficient fresh hydrogen to the core of the collision product to allow the collision product to burn for several million years. We find that grazing encounters enter a common--envelope phase which may expel the envelope of the merger product. The deposition of energy in the envelopes of our merger products causes them to swell by factors of $\sim100$. If these remnants exist in very densely-populated environments ($n\gtrsim10^{7}$ pc$^{-3}$), they will suffer further collisions which may drive off their envelopes, leaving behind hard binaries. We show that the products of collisions have cores rotating sufficiently rapidly to make them candidate hypernova/gamma--ray burst progenitors and that $\sim0.1%$ of massive stars may suffer collisions, sufficient for such events to contribute significantly to the observed rates of hypernovae and gamma--ray bursts.

Recent progress on the nature of short duration gamma-ray bursts has shown that a fraction of them originate in the local universe. These systems may well be the result of giant flares from soft gamma-repeaters (highly magnetized neutron stars commonly known as magnetars). However, if these neutron stars are formed via the core collapse of massive stars then it would be expected that the bursts should originate from predominantly young stellar populations, while correlating the positions of BATSE short bursts with structure in the local universe reveals a correlation with all galaxy types, including those with little or no ongoing star formation. This is a natural outcome if, in addition to magnetars formed via the core collapse of massive stars they also form via Accretion Induced Collapse following the merger of two white dwarfs, one of which is magnetic. We investigate this possibility and find that the rate of magnetar production via WD-WD mergers in the Milky Way is comparable to the rate of production via core-collapse. However, while the rate of production of magnetars by core collapse is proportional to the star formation rate, the rate of production via WD-WD mergers (which have long lifetimes) is proportional to the stellar mass density, which is concentrated in early-type systems. Therefore magnetars produced via WD-WD mergers may produce SGR giant flares which can be identified with early type galaxies. We also comment on the possibility that this mechanism could produce a fraction of the observed short duration burst population at higher redshift.

Two short (<2 s) g-ray bursts (GRBs) have recently been localized and fading afterglow counterparts detected. The combination of these two results left unclear the nature of the host galaxies of the bursts, because one was a star-forming dwarf, while the other was probably an elliptical galaxy. Here we report the X-ray localization of a short burst (GRB 050724) with unusual g-ray and X-ray properties. The X-ray afterglow lies off the centre of an elliptical galaxy at a redshift of z=0.258, coincident with the position determined by ground-based optical and radio observations. The low level of star formation typical for elliptical galaxies makes it unlikely that the burst originated n a supernova explosion. A supernova origin was also ruled out for GRB 050709, even though that burst took place in a galaxy with current star formation. The isotropic energy for the short bursts is 2-3 orders of magnitude lower han that for the long bursts. Our results therefore suggest that an alternative source of bursts -- the of binary systems of neutron stars or a neutron star-black hole pair -- are the progenitors of short bursts.

Attempts of Magain et al (2005) to detect the host galaxy of the bright QSO HE0450--2958 have not been successful. We suggest that the supermassive black hole powering the QSO was ejected from the observed ULIRG at the same redshift and at 1.5 arcsec distance. Ejection could have either be caused by recoil due to gravitational wave emission from a coalescing binary of supermassive black holes or the gravitational slingshot of three or more supermassive black holes in the ongoing merger of galaxies which triggered the starburst activity in the ULIRG. We discuss implications for the possible hierarchical build-up of supermassive black holes from intermediate and/or stellar mass black holes, and for the detection of coalescing supermassive binary black holes by LISA.

We consider the origin of the so-called S stars orbiting the supermassive black hole at the very center of the Galaxy. These are usually assumed to be massive main-sequence stars. We argue instead that they are the remnants of low-to-intermediate mass red giants which have been scattered on to near-radial orbits and tidally stripped as they approach the central black hole. Such stars retain only low-mass envelopes and thus have high effective temperatures. Our picture simultaneously explains why S stars have tightly-bound orbits, and the observed depletion of red giants in the very center of the Galaxy.

We investigate the structure and evolution of blue stragglers stars which were formed from direct stellar collisions between main sequence stars in globular clusters. In particular, we look at the rotational evolution of the products of off-axis collisions. As found in previous work, such blue stragglers initially have too high an angular momentum to contract down to the main sequence. We consider angular momentum loss through either disc locking or locking to an outflowing wind and show that both methods allow the merged object to shed sufficient angular momentum to contract down to the main sequence.

The globular cluster M15 was observed on three occasions with the High Resolution Camera on board Chandra in 2001 in order to investigate the X-ray source population in the cluster centre. After subtraction of the two bright central sources, four faint sources were identified within 50 arcsec of the core. One of these sources is probably the planetary nebula, K648, making this the first positive detection of X-rays from a planetary nebula inside a globular cluster. Another two are identified with UV variables (one previously known), which we suggest are cataclysmic variables (CVs). The nature of the fourth source is more difficult to ascertain, and we discuss whether it is possibly a quiescent soft X-ray transient (qSXT) or also a CV.

We present a simple, semi--analytical description for the final stages of mergers of black hole (BH) -- neutron star (NS) systems. Such systems are of much interest as gravitational wave sources and gamma--ray burst progenitors. Numerical studies show that in general the neutron star is not disrupted at the first phase of mass transfer. Instead, what remains of the neutron star is left on a wider, eccentric, orbit. We consider the evolution of such systems as they lose angular momentum via gravitational radiation and come into contact for further phases of mass transfer. During each mass transfer event the neutron star mass is reduced until a critical value where mass loss leads to a rapid increase in the stellar radius. At this point Roche lobe overflow shreds what remains of the neutron star, most of the mass forming a disc around the black hole. Such a disc may be massive enough to power a gamma--ray burst. The mass of the neutron star at the time of disruption (and therefore the disc mass) is largely independent of the initial masses of the black hole and neutron star, indicating that BH--NS star mergers may be standard candles.

We discuss the main properties of the Galactic globular cluster (GC) blue straggler stars (BSS), as inferred from our new catalog containing nearly 3000 BSS. The catalog has been extracted from the photometrically homogeneous V vs. (B-V) color-magnitude diagrams (CMD) of 56 GCs, based on WFPC2 images of their central cores. In our analysis we used consistent relative distances based on the same photometry and calibration. The number of BSS has been normalized to obtain relative frequencies (F_{BSS}) and specific densities (N_S) using different stellar populations extracted from the CMD. The cluster F_{BSS} is significantly smaller than the relative frequency of field BSS. We find a significant anti-correlation between the BSS relative frequency in a cluster and its total absolute luminosity (mass). There is no statistically significant trend between the BSS frequency and the expected collision rate. F_{BSS} does not depend on other cluster parameters, apart from a mild dependence on the central density. PCC clusters act like normal clusters as far as the BSS frequency is concerned. We also show that the BSS luminosity function for the most luminous clusters is significantly different, with a brighter peak and extending to brighter luminosities than in the less luminous clusters. These results imply that the efficiency of BSS production mechanisms and their relative importance vary with the cluster mass.

Recent HST observations of a large sample of globular clusters reveal that every cluster contains between 40 and 400 blue stragglers. The population does not correlate with either stellar collision rate (as would be expected if all blue stragglers were formed via collisions) or total mass (as would be expected if all blue stragglers were formed via the unhindered evolution of a subset of the stellar population). In this paper, we support the idea that blue stragglers are made through both channels. The number produced via collisions tends to increase with cluster mass. In this paper we show how the current population produced from primordial binaries decreases with increasing cluster mass; exchange encounters with third, single, stars in the most massive clusters tend to reduce the fraction of binaries containing a primary close to the current turn-off mass. Rather their primaries tend to be somewhat more massive (~1-3 M_sun) and have evolved off the main sequence, filling their Roche lobes in the past, often converting their secondaries into blue stragglers (but more than 1 Gyr or so ago and thus they are no longer visible as blue stragglers). We show that this decline in the primordial blue straggler population is likely to be offset by the increase in the number of blue stragglers produced via collisions. The predicted total blue straggler population is therefore relatively independent of cluster mass, thus matching the observed population. This result does not depend on any particular assumed blue straggler lifetime.

Recent three dimensional, high-resolution simulations of neutron star coalescences are analysed to assess whether short gamma-ray bursts (GRBs) could originate from such encounters. The two most popular modes of energy extraction -- namely the annihilation of $\nu\bar{\nu}$ and magnetohydrodynamic processes- are explored in order to investigate their viability in launching GRBs. We find that $\nu\bar{\nu}$ annihilation can provide the necessary stresses to drive a highly relativistic expansion. But unless the outflow is beamed into less than one percent of the solid angle this mechanism may fail to explain the apparent isotropized energies implied for short GRBs at cosmological distances. We argue that the energetic, neutrino-driven wind that accompanies the merger event will have enough pressure to provide adequate collimation to the $\nu\bar{\nu}$-annihilation-driven jet, thereby comfortably satisfying constraints on event rate and apparent luminosity. We also assess magnetic mechanisms to transform the available energy into a GRB. If the central object does not collapse immediately into a black hole it will be convective and it is expected to act as an effective large scale dynamo amplifying the seed magnetic fields to a few times $10^{17}$ G within a small fraction of a second. The associated spindown time scale is 0.2 s, coinciding with the typical duration of a short GRB. The efficiencies of the various assessed magnetic processes are high enough to produce isotropized luminosities in excess of $10^{52}$ ergs/s even without beaming.

The incidence of evaporating 'black widow' pulsars (BWPs) among all millisecond pulsars (MSPs) is far higher in globular clusters than in the field. This implies a special formation mechanism for them in clusters. Cluster MSPs in wide binaries with WD companions exchange them for turnoff-mass stars. These new companions eventually overflow their Roche lobes because of encounters and tides. The millisecond pulsars eject the overflowing gas from the binary, giving mass loss on the binary evolution timescale. The systems are only observable as BWPs at epochs where this evolution is slow, making the mass loss transparent and the lifetime long. This explains why observed BWPs have low-mass companions. We suggest that at least some field BWPs were ejected from globular clusters or entered the field population when the cluster itself was disrupted.

We investigate the observed depletion of red giants in the cores of post-core-collapse globular clusters. In particular, the evolutionary scenario we consider is a binary consisting of two low-mass stars which undergoes two common envelope phases. The first common envelope phase occurs when the primary is a red giant resulting in a helium white dwarf and main sequence star in a detached binary. The second common envelope phase occurs shortly after the secondary becomes a red giant. During the second common envelope phase the degenerate helium cores merge resulting in a core mass greater than the helium burning limit and the formation of a horizontal branch star. We show that this evolutionary route is enhanced in post-core-collapse clusters by stellar encounters. These encounters increase the population of binary secondaries which would have evolved onto the red giant branch in the recent past.

The binaries PSR J1141-6545 and PSR B2303+46 each appear to contain a white dwarf which formed before the neutron star. We describe an evolutionary pathway to produce these two systems. In this scenario, the primary transfers its envelope onto the secondary which is then the more massive of the two stars, and indeed sufficiently massive later to produce a neutron star via a supernova. The core of the primary produces a massive white dwarf which enters into a common envelope with the core of the secondary when the latter evolves off the main sequence. During the common envelope phase, the white dwarf and the core of the secondary spiral together as the envelope is ejected. The evolutionary history of PSR J1141-6545 and PSR B2303+46 differ after this phase. In the case of PSR J1141--6545, the secondary (now a helium star) evolves into contact transferring its envelope onto the white dwarf. We propose that the vast majority of this material is in fact ejected from the system. The remains of the secondary then explode as a supernova producing a neutron star. Generally the white dwarf and neutron star will remain bound in tight, often eccentric, systems resembling PSR J1141-6545. These systems will spiral in and merge on a relatively short timescale and may make a significant contribution to the population of gamma ray burst progenitors. In PSR B2303+46, the helium-star secondary and white dwarf never come into contact. Rather the helium star loses its envelope via a wind, which increases the binary separation slightly. Only a small fraction of such systems will remain bound when the neutron star is formed (as the systems are wider). Those systems which are broken up will produce a population of high-velocity white dwarfs and neutron stars.

We suggest that the collapsing core of a massive rotating star may fragment to produce two or more compact objects. Their coalescence under gravitational radiation gives the resulting black hole or neutron star a significant kick velocity, which may explain those observed in pulsars. A gamma-ray burst can result only when this kick is small. Thus only a small fraction of core-collapse supernovae produce gamma-ray bursts. The burst may be delayed significantly (hours - days) after the supernova, as suggested by recent observations. If this picture is correct, core-collapse supernovae should be significant sources of gravitational radiation with a chirp signal similar to a coalescing binary.

We present the results of 3D, high-resolution calculations of neutron star mergers. Using Newtonian gravity, but adding the forces from the emission of gravitational waves, we solve the hydrodynamics equations using the smoothed particle hydrodynamics method. The microscopic properties of the neutron star matter are described in terms of a new nuclear equation of state (EOS) based on the tables of Shen et al. (1998). The effects of neutrino emission of all flavours are accounted for. The new EOS is noticeably stiffer than the Lattimer-Swesty-EOS used in previous investigations. We find in general somewhat lower temperatures than previous investigations. This allows a substantial fraction of heavy nuclei to be present even in the inner parts of the torus that forms around the central object and yields important contributions to the neutrino opacity. The central object shows a rapid differential rotation with periods of $\sim$ 2ms, which we suspect to stabilize the central object on (at least) the simulation time scale against gravitational collapse. The differential rotation could wind up initial magnetic fields to $\sim 10^17$ G.

We present the results of multiple simulations of open clusters, modelling the dynamics of a population of brown dwarf members. We consider the effects of a large range of primordial binary populations, including the possibilities of having brown dwarf members contained within a binary system. We also examine the effects of various cluster diameters and masses. Our examination of a population of wide binary systems containing brown dwarfs, reveals evidence for exchange reactions whereby the brown dwarf is ejected from the system and replaced by a heavier main-sequence star. We find that there exists the possibility of hiding a large fraction of the brown dwarfs contained within the primordial binary population. We conclude that it is probable that the majority of brown dwarfs are contained within primordial binary systems which then hides a large proportion of them from detection.

We investigate whether the recently-observed population of high-velocity white dwarfs can be derived from a population of binaries residing initially within the thin disk of the Galaxy. In particular we consider binaries where the primary is sufficiently massive to explode as a type II supernova. A large fraction of such binaries are broken up when the primary then explodes as a supernova owing to the combined effects of the mass loss from the primary and the kick received by the neutron star on its formation. For binaries where the primary evolves to fill its Roche lobe, mass transfer from the primary leads to the onset of a common envelope phase during which the secondary and the core of the primary spiral together as the envelope is ejected. Such binaries are the progenitors of X-ray binaries if they are not broken up when the primary explodes. For those systems which are broken up, a large number of the secondaries receive kick velocities ~100 - 200 km/s and subsequently evolve into white dwarfs. We compute trajectories within the Galactic potential for this population of stars and relate the birthrate of these stars over the entire Galaxy to those seen locally with high velocities relative to the LSR. We show that for a reasonable set of assumptions concerning the Galactic supernova rate and the binary population, our model produces a local number density of high-velocity white dwarfs compatible with that inferred from observations. We therefore propose that a population of white dwarfs originating in the thin disk may make a significant contribution to the observed population of high-velocity white dwarfs.

We performed high-resolution simulations of two stellar collisions relevant for stars in globular clusters. We considered one head-on collision and one off-axis collision between two 0.6 M_sun main sequence stars. We show that a resolution of about 100 000 particles is sufficient for most studies of the structure and evolution of blue stragglers. We demonstrate conclusively that collision products between main-sequence stars in globular clusters do not have surface convection zones larger than 0.004 M_sun after the collision, nor do they develop convection zones during the `pre-main-sequence' thermal relaxation phase of their post-collision evolution. Therefore, any mechanism which requires a surface convection zone (i.e. chemical mixing or angular momentum loss via a magnetic wind) cannot operate in these stars. We show that no disk of material surrounding the collision product is produced in off-axis collisions. The lack of both a convection zone and a disk proves a continuing problem for the angular momentum evolution of blue stragglers in globular clusters.

The importance of stellar encounters in producing stellar exotica in dense stellar clusters is reviewed. We discuss how collisions between main-sequence stars may be responsible for the production of blue stragglers in globular clusters. We also discuss the possible pathways to the production of X-ray binaries, cataclysmic variables, and millisecond pulsars in globular clusters. Neutron stars in globular clusters are likely to exchange into binaries containing moderate-mass main-sequence stars, replacing the lower-mass components of the original systems. These binaries will become intermediate-mass X-ray binaries (IMXBs), once the moderate-mass star evolves off the main-sequence, as mass is transferred onto the neutron star possibly spinning it up in the process. Such systems may be responsible for the population of millisecond pulsars (MSPs) that has been observed in globular clusters. Additionally, the period of mass-transfer (and thus X-ray visibility) in the vast majority of such systems will have occurred 5 - 10 Gyr ago thus explaining the observed relative paucity of X-ray binaries today, given the large MSP population.

We investigate models for the class of ultraluminous non-nuclear X-ray sources (ULXs) seen in a number of galaxies and probably associated with star-forming regions. Models where the X-ray emission is assumed to be isotropic run into several difficulties. In particular formation of sufficient numbers of the required ultramassive black-hole X-ray binaries is problematic, and the likely transient behaviour of the resulting systems is not in good accord with observation. The assumption of mild X-ray beaming suggests instead that ULXs may represent a shortlived but extremely common stage in the evolution of a wide class of X-ray binaries. The best candidate for this is the phase of thermal-timescale mass transfer inevitable in many intermediate and high-mass X-ray binaries. This in turn suggests a link with the Galactic microquasars. The short lifetimes of high-mass X-ray binaries would explain the association of ULXs with episodes of star formation. These considerations still allow the possibility that individual ULXs may contain extremely massive black holes.

We consider the survivability of planetary systems in the globular cluster 47 Tucanae. We compute the cross sections for the breakup of planetary systems via encounters with single stars and binaries. We also compute the cross sections to leave planets on eccentric orbits. We find that wider planetary systems (d > 0.3 AU) are likely to be broken up in the central regions of 47 Tucanae (within the half-mass radius of the cluster). However tighter systems, and those in less-dense regions may survive. Tight systems will certainly survive in less-dense clusters where subsequent surveys should be conducted.

We introduce an adaptation of the well known Tree+SPH numerical scheme to Post Newtonian (PN) hydrodynamics and gravity. Our code solves the (0+1+2.5)PN equations. These equations include Newtonian hydrodynamics and gravity (0PN), the first order relativistic corrections to those (1PN) and the lowest order gravitational radiation terms (2.5PN). We test various aspects of our code using analytically solvable test problems. We then proceed to study the 1PN effects on binary neutron star coalescence by comparing calculations with and without the 1PN terms. We find that the effect of the 1PN terms is rather small. The largest effect arises with a stiff equation of state for which the maximum rest mass density increases. This could induce black hole formation. The gravitational wave luminosity is also affected.

We investigate how the formation and evolution of extrasolar planetary systems can be affected by stellar encounters that occur in the crowded conditions of a stellar cluster. Using plausible estimates of cluster evolution, we show how planet formation may be supressed in globular clusters while planets wider than 0.1 AU that do form in such environments can be ejected from their stellar system. Less crowded systems such as open clusters have a much reduced effect on any planetary system. Planet formation is unaffected in open clusters and only the wider planetary systems will be disrupted during the cluster's lifetime. The potential for free-floating planets in these environments is also discussed.

We show that the usual picture of supersoft X-ray binary evolution as driven by conservative thermal-timescale mass transfer cannot explain the short orbital periods of RX J0537.7-7034 (3.5 hr) and 1E 0035.4-7230 (4.1 hr). Non-conservative evolution may produce such periods, but requires very significant mass loss, and is highly constrained.

Results of 3D, Newtonian hydrodynamic calculations of the last stages of the inspiral and the final coalescence of neutron star binary systems are described. The focus is on sligthly asymmetric systems with the asymmetry stemming either from different spins or different masses of the binary components. Implications for nucleosynthesis and gamma ray bursts are discussed. Test calculations for the used viscosity scheme are provided in the appendix.

Cygnus X-2 appears to be the descendant of an intermediate-mass X-ray binary (IMXB). Using Mazzitelli's (1989) stellar code we compute detailed evolutionary sequences for the system and find that its prehistory is sensitive to stellar input parameters, in particular the amount of core overshooting during the main-sequence phase. With standard assumptions for convective overshooting a case B mass transfer starting with a 3.5 M_sun donor star is the most likely evolutionary solution for Cygnus X-2. This makes the currently observed state rather short-lived, of order 3 Myr, and requires a formation rate > 1e-7 - 1e-6 per yr of such systems in the Galaxy. Our calculations show that neutron star IMXBs with initially more massive donors (> 4 M_sun) encounter a delayed dynamical instability; they are unlikely to survive this rapid mass transfer phase. We determine limits for the age and initial parameters of Cygnus X-2 and calculate possible dynamical orbits of the system in a realistic Galactic potential, given its observed radial velocity. We find trajectories which are consistent with a progenitor binary on a circular orbit in the Galactic plane inside the solar circle that received a kick velocity < 200 km/s at the birth of the neutron star. The simulations suggests that about 7% of IMXBs receiving an arbitrary kick velocity from a standard kick velocity spectrum would end up in an orbit similar to Cygnus X-2, while about 10% of them reach yet larger Galactocentric distances.

We simulate collisions involving red-giant stars in the centre of our galaxy. Such encounters may explain the observed paucity of highly luminous red giants within the central 0.2pc. The masses of the missing stars are likely to be in the range 2 to 8 solar masses. Recent models of the galactic centre cluster's density and velocity distributions are used to calculate two-body collision rates. In particular we use stellar-evolution models to calculate the number of collisions a star will have during different evolutionary phases. We find that the number of two-body collisions per star is $\lo 1$ in the central 0.1 to 0.2 pc, depending strongly on the galactocentric radius. Using a 3D numerical hydrodynamics code (SPH) we simulate encounters involving cluster stars of various masses with 2 and 8 solar-mass red giants. The instantaneous mass loss in such collisions is rarely enough to destroy either giant. A fraction of the collisions do, however, lead to the formation of common envelope systems where the impactor and giant's core are enshrouded by the envelope of the giant. Such systems may evolve to expel the envelope, leaving a tight binary; the original giant is destroyed. The fraction of collisions that produce common envelope systems is sensitive to the local velocity dispersion and hence galactocentric radius. Using our collision-rate calculations we compute the time-scales for a giant star to suffer such a collision within the galactic centre. These time-scales are >10^{9-10}years and so are longer than the lifetimes of stars more-massive than 2 solar masses. Thus the observed paucity of luminous giants is unlikely to be due to the formation of common envelope systems as a result of two-body encounters involving giant stars.

The merger of compact binaries, especially black holes and neutron stars, is frequently invoked to explain gamma-ray bursts (GRB's). In this paper, we present three dimensional hydrodynamical simulations of the relatively neglected mergers of white dwarfs and black holes. During the merger, the white dwarf is tidally disrupted and sheared into an accretion disk. Nuclear reactions are followed and the energy release is negligible. Peak accretion rates are ~0.05 Msun/s (less for lower mass white dwarfs) lasting for approximately a minute. Many of the disk parameters can be explained by a simple analytic model which we derive and compare to our simulations. This model can be used to predict accretion rates for white dwarf and black hole (or neutron star) masses which are not simulated in this paper. Although the mergers studied here create disks with larger radii, and longer accretion times than those from the merger of double neutron stars, a larger fraction of the merging star's mass becomes part of the disk. Thus the merger of a white dwarf and a black hole could produce a long duration GRB. The event rate of these mergers may be as high as 1/Myr per galaxy.

We present the results of 3D Newtonian SPH simulations of the merger of a neutron star binary. The microscopic properties of matter are described by the physical equation of state of Lattimer and Swesty (LS-EOS). To test for the robustness of our results we check the sensitivity to the approximations of our model as well as to the binary system parameters. The main and new result is that for the realistic LS-EOS, depending on the initial spin, between 4e-3 and 4e-2 solar masses of material become unbound. If, as suggested, large parts of this matter consist of r-process nuclei, neutron star mergers could account for the whole observed r-process material in the Galaxy.

Recent calculations place questions on the ability of supernovae to produce r-process nuclei in the correct amounts. We present results from 3D Newtonian SPH calculations of mergers of equal mass neutron star binaries. We find the amounts of ejected mass to be very promising for an explanation of the observed r-process abundances. Preliminary abundance calculations show that practically all the material is subject to r-process conditions. The calculated abundance patterns reproduce the basic features of the solar r-process abundances very well. Especially the peak around A= 195 is easily reproduced without any tuning of the initial entropies. Thus, eventually all the observed r-process material could be explained by mass ejection during neutron star mergers.

We have identified a candidate 1-2 x 10^5 year old luminous white dwarf in NGC 1818, a young star cluster in the Large Magellanic Cloud. This discovery strongly constrains the boundary mass M_c at which stars stop forming neutron stars and start forming white dwarfs, to M_c > 7.6 Msun.

We explore the role that red giants might play in the central regions of Active Galactic Nuclei. Due to their large radii and the low binding energy of the stellar envelope, giants are vulnerable to envelope stripping from collisions with the accretion disk. Using hydrodynamic simulations we show that such collisions will typically deposit a substantial fraction of the envelope mass into the disk on each passage. Repeated encounters will then lead to the complete destruction of the star save for the dense core. We estimate the rate of fuel supply by this mechanism using simple models for the AGN disk and central stellar cluster. If the central stellar density is of order 10^7 solar masses per cubic pc, then stripping of giants could account for the activity of typical AGN provided that the accretion disk extends beyond 0.1 pc. For AGN with smaller disks, or clusters of lower central density, giant stripping could be an important source of gas enriched via stellar nucleosynthesis.

We report on an extensive series of numerical experiments of binary--binary scattering, analysing the cross--section for close approach during interactions for a range of hard binary parameters of interest in globular cluster cores. We consider the implied rate for tidal interactions for different globular clusters and compare our results with previous, complementary estimates of stellar collision rates in globular clusters. We find that the collision rate for binary--binary encounters dominates in low density clusters if the binary fraction in the cluster is larger than $0.2$ for wide main--sequence binaries. In dense clusters binary--single interactions dominate the collision rate and the core binary fraction must be $\ltorder 0.1$ per decade in semi--major axis or too many collisions take place compared to observations. The rates are consistent if binaries with semi--major axes $\sim 100 AU$ are overabundant in low density clusters or if breakup and ejection substantially lowers the binary fraction in denser clusters. Given reasonable assumptions about fractions of binaries in the cores of low density clusters such as NGC~5053, we cannot account for all the observed blue stragglers by stellar collisions during binary encounters, suggesting a substantial fraction may be due to coalescence of tight primordial binaries.

Using King-Mitchie Models, we compute encounter rates between the various stellar species in the globular clusters $\omega$ Cen, and 47 Tuc. We also compute event rates for encounters between single stars and a population of primordial binaries. Using these rates, and what we have learnt from hydrodynamical simulations of encounters performed earlier, we compute the production rates of objects such as low-mass X-ray binaries (LMXBs), smothered neutron stars and blue stragglers (massive main-sequence stars). If 10\% of the stars are contained in primordial binaries, the production rate of interesting objects from encounters involving these binaries is as large as that from encounters between single stars. For example, encounters involving binaries produce a significant number of blue stragglers in both globular cluster models. The number of smothered neutron stars may exceed the number of low-mass X-ray binaries (LMXBs) by a factor of 5-20, which may help explain why millisecond pulsars are observed to outnumber LMXBs in globular clusters.

In studying encounters between binaries and single stars, one is interested in three classes of events: exchanges of stars, hardening of the original binary by a third star, and the production of merged objects. We present a means for computing cross sections for these three outcomes for an arbitrary binary and single star as might be found in the core of a globular cluster. The cross sections for a number of binaries in various stellar populations are then computed. We consider multiple encounters and the ultimate fate of a population of binaries fed into the cores of different globular cluster models. We see that the presence of only a relatively small number of binaries (containing 10\% of the stars) will boost the production rate of astrophysically-interesting objects by a factor of at least a few over the rates expected from encounters between single stars. In particular, the ratio of smothered neutron stars to low-mass X-ray binaries (LMXBs) may be greatly increased, possibly explaining, in part, the excess of millisecond pulsars compared to LMXBs.

Recent surveys of the blue straggler (BS) population in the galactic globular cluster M3 (NGC 5272) give the first complete characterization of the number density of BSs as a function of radius over an entire globular cluster. The BSs in M3 are overabundant at large radii, and significantly underabundant at intermediate radii. Here we present the result of a simulation of the dynamical evolution of a population of BSs in a multi--mass model of M3. Assuming the BSs were formed in the core through binary interactions (Hut \&\ Verbunt 1983, Leonard 1989, Sigurdsson \&\ Phinney 1993, Hut \etal 1992, Davies, Benz and Hills 1994), and given some very general assumptions about the recoil that occurs during stellar mergers in interacting binaries; we find an excellent fit to the observed radial distribution of BSs, suggesting strongly that most of the BSs in M3 were formed through binary collisions in the core.

We present 3D Newtonian simulations of the coalescence of two neutron stars, using a Smoothed Particle Hydrodynamics (SPH) code. We begin the simulations with the two stars in a hard, circular binary, and have them spiral together as angular momentum is lost through gravitational radiation at the rate predicted by modeling the system as two point masses. We model the neutron stars as hard polytropes ($\gamma=2.4$) of equal mass, and investigate the effect of the initial spin of the two stars on the coalescence. The process of coalescence, from initial contact to the formation of an axially symmetric object, takes only a few orbital periods. Some of the material from the two neutron stars is shed, forming a thick disk around the central, coalesced object. The mass of this disk is dependent on the initial neutron star spins; higher spin rates resulting in greater mass loss, and thus more massive disks. For spin rates that are most likely to be applicable to real systems, the central coalesced object has a mass of $2.4M_\odot$, which is tantalizingly close to the maximum mass allowed by any neutron star equation of state for an object that is supported in part by rotation. Using a realistic nuclear equation of state, we estimate the temperature of the material after the coalescence. We find that the central object is at a temperature of $\sim 10$MeV, whilst the disk is heated by shocks to a temperature of 2-4MeV.