Source author record

Melvyn B. Davies

Melvyn B. Davies appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

40works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

40 published item(s)

preprint2022arXiv

A pair of Sub-Neptunes transiting the bright K-dwarf TOI-1064 characterised with CHEOPS

We report the discovery and characterisation of a pair of sub-Neptunes transiting the bright K-dwarf TOI-1064 (TIC 79748331), initially detected in TESS photometry. To characterise the system, we performed and retrieved CHEOPS, TESS, and ground-based photometry, HARPS high-resolution spectroscopy, and Gemini speckle imaging. We characterise the host star and determine $T_{\rm eff, \star}=4734\pm67$ K, $R_{\star}=0.726\pm0.007$ $R_{\odot}$, and $M_{\star}=0.748\pm0.032$ $M_{\odot}$. We present a novel detrending method based on PSF shape-change modelling and demonstrate its suitability to correct flux variations in CHEOPS data. We confirm the planetary nature of both bodies and find that TOI-1064 b has an orbital period of $P_{\rm b}=6.44387\pm0.00003$ d, a radius of $R_{\rm b}=2.59\pm0.04$ $R_{\oplus}$, and a mass of $M_{\rm b}=13.5_{-1.8}^{+1.7}$ $M_{\oplus}$, whilst TOI-1064 c has an orbital period of $P_{\rm c}=12.22657^{+0.00005}_{-0.00004}$ d, a radius of $R_{\rm c}=2.65\pm0.04$ $R_{\oplus}$, and a 3$σ$ upper mass limit of 8.5 ${\rm M_{\oplus}}$. From the high-precision photometry we obtain radius uncertainties of $\sim$1.6%, allowing us to conduct internal structure and atmospheric escape modelling. TOI-1064 b is one of the densest, well-characterised sub-Neptunes, with a tenuous atmosphere that can be explained by the loss of a primordial envelope following migration through the protoplanetary disc. It is likely that TOI-1064 c has an extended atmosphere due to the tentative low density, however further RVs are needed to confirm this scenario and the similar radii, different masses nature of this system. The high-precision data and modelling of TOI-1064 b are important for planets in this region of mass-radius space, and it allows us to identify a trend in bulk density-stellar metallicity for massive sub-Neptunes that may hint at the formation of this population of planets.

preprint2022arXiv

Uncovering the true periods of the young sub-Neptunes orbiting TOI-2076

Context: TOI-2076 is a transiting three-planet system of sub-Neptunes orbiting a bright (G = 8.9 mag), young ($340\pm80$ Myr) K-type star. Although a validated planetary system, the orbits of the two outer planets were unconstrained as only two non-consecutive transits were seen in TESS photometry. This left 11 and 7 possible period aliases for each. Aims: To reveal the true orbits of these two long-period planets, precise photometry targeted on the highest-probability period aliases is required. Long-term monitoring of transits in multi-planet systems can also help constrain planetary masses through TTV measurements. Methods: We used the MonoTools package to determine which aliases to follow, and then performed space-based and ground-based photometric follow-up of TOI-2076 c and d with CHEOPS, SAINT-EX, and LCO telescopes. Results: CHEOPS observations revealed a clear detection for TOI-2076 c at $P=21.01538^{+0.00084}_{-0.00074}$ d, and allowed us to rule out three of the most likely period aliases for TOI-2076 d. Ground-based photometry further enabled us to rule out remaining aliases and confirm the $P=35.12537\pm0.00067$ d alias. These observations also improved the radius precision of all three sub-Neptunes to $2.518\pm0.036$, $3.497\pm0.043$, and $3.232\pm0.063$ $R_\oplus$. Our observations also revealed a clear anti-correlated TTV signal between planets b and c likely caused by their proximity to the 2:1 resonance, while planets c and d appear close to a 5:3 period commensurability, although model degeneracy meant we were unable to retrieve robust TTV masses. Their inflated radii, likely due to extended H-He atmospheres, combined with low insolation makes all three planets excellent candidates for future comparative transmission spectroscopy with JWST.

preprint2021arXiv

Formation of super-massive black holes in galactic nuclei I: delivering seed intermediate-mass black holes in massive stellar clusters

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, and investigate the outcome of 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$. 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.

preprint2021arXiv

Hot Jupiters, cold kinematics: High phase space densities of host stars reflect an age bias

Context. The birth environments of planetary systems are thought to influence planet formation and orbital evolution, through external photoevaporation and stellar flybys. Recent work has claimed observational support for this, in the form of a correlation between the properties of planetary systems and the local Galactic phase space density of the host star. In particular, Hot Jupiters are found overwhelmingly around stars in regions of high phase space density, which may reflect a formation environment with high stellar density. Aims. We instead investigate whether the high phase space density may have a galactic kinematic origin: Hot Jupiter hosts may be biased towards being young and therefore kinematically cold, because tidal inspiral leads to the destruction of the planets on Gyr timescales, and the velocity dispersion of stars in the Galaxy increases on similar timescales. Methods. We use 6D positions and kinematics from Gaia for the Hot Jupiter hosts and their neighbours, and construct distributions of the phase space density. We investigate correlations between the stars' local phase space density and peculiar velocity. Results. We find a strong anticorrelation between the phase space density and the host star's peculiar velocity with respect to the Local Standard of Rest. Therefore, most stars in "high-density" regions are kinematically cold, which may be caused by the aforementioned bias towards detecting Hot Jupiters around young stars before the planets' tidal destruction. Conclusions. We do not find evidence in the data for Hot Jupiter hosts preferentially being in phase space overdensities compared to other stars of similar kinematics, nor therefore for their originating in birth environments of high stellar density.

preprint2021arXiv

Metal pollution of the solar white dwarf by solar system small bodies

White dwarfs (WDs) often show metal lines in their spectra, indicating accretion of asteroidal material. Our Sun is to become a WD in several Gyr. Here, we examine how the solar WD accretes from the three major small body populations: the main belt asteroids (MBAs), Jovian trojan asteroids (JTAs), and trans-Neptunian objects (TNOs). Owing to the solar mass loss during the giant branch, 40\% of the JTAs are lost but the vast majority of MBAs and TNOs survive. During the WD phase, objects from all three populations are sporadically scattered onto the WD, implying ongoing accretion. For young cooling ages $\lesssim 100$ Myr, accretion of MBAs predominates; our predicted accretion rate $\sim10^6$ g/s falls short of observations by two orders of magnitude. On Gyr timescales, thanks to the consumption of the TNOs that kicks in $\gtrsim 100$ Myr, the rate oscillates around $10^6-10^7$ g/s until several Gyr and drops to $\sim10^5$ g/s at 10 Gyr. Our solar WD accretion rate from 1 Gyr and beyond agrees well with those of the extrasolar WDs. We show that for the solar WD, the accretion source region evolves in an inside-out pattern. Moreover, in a realistic small body population with individual sizes covering a wide range as WD pollutants, the accretion is dictated by the largest objects. As a consequence, the accretion rate is lower by an order of magnitude than that from a population of bodies of a uniform size and the same total mass and shows greater scatter.

preprint2020arXiv

Capture of satellites during planetary encounters A case study of the Neptunian moons Triton and Nereid

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.

preprint2020arXiv

Effects of capturing a wide-orbit planet on planetary systems: system stability and Habitable Zone bombardment rates

A large fraction of stars are formed in dense clusters. In the cluster, close encounters between stars at distances less than 100 au are common. It has been shown that during close encounters planets can transfer between stars. Such captured planets will be on different orbits compared to planets formed in the system, often on very wide, eccentric and inclined orbits. We examine how these captured planets affect Kuiper-belt like asteroid belts in their new systems, and how this affects habitable planets in the system. We show that these captured planets can destabilize the asteroid belt, and we show that the fraction of the asteroid that make it past the giant planets into the system to impact the habitable planet is independent of the captured planets orbital plane, whereas the fraction of the asteroids that are removed and the rate at which they are removed depend strongly on the captured planets pericentre and inclination. We then examine all possible outcomes of planet capture and find that when a Jupiter-mass planet is captured it will in 40\% of cases destabilize the planets in the system, in 40\% of cases deplete the asteroid belt in a few Myr, i.e. not posing much risk to life on terrestrial planets which would be expected to develop later. In the final 20\% of cases the result will be a flux of impactors 5-10 times greater than that on Earth that can persist for several Gyr, quite detrimental to the development of life on the planet.

preprint2020arXiv

Fly-by encounters between two planetary systems II: Exploring the interactions of diverse planetary system architectures

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.

preprint2020arXiv

Global torques and stochasticity as the drivers of massive black hole pairing in the young Universe

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.

preprint2020arXiv

Making massive stars in the Galactic Centre via accretion onto low-mass stars within an accretion disc

The origin of the population of very massive stars observed within $\sim 0.4$ pc of the supermassive black hole in the Galactic Centre is a mystery. Tidal forces from the black hole would likely inhibit {\it in situ} star formation whilst the youth of the massive stars would seem to exclude formation elsewhere followed by transportation (somehow) into the Galactic centre. Here we consider a third way to produce these massive stars from the lower-mass stars contained in the nuclear stellar cluster which surrounds the supermassive black hole. A passing gas cloud can be tidally shredded by the supermassive black hole forming an accretion disc around the black hole. Stars embedded within this accretion disc will accrete gas from the disc via Bondi-Hoyle accretion, where the accretion rate onto a star, $\dot{M}_\star \propto M_\star^2$. This super-exponential growth of accretion can lead to a steep increase in stellar masses, reaching the required 40-50 M$_\odot$ in some cases. The mass growth rate depends sensitively on the stellar orbital eccentricities and their inclinations. The evolution of the orbital inclinations and/or their eccentricities as stars are trapped by the disc, and their orbits are circularised, will increase the number of massive stars produced. Thus accretion onto low-mass stars can lead to a top heavy stellar mass function in the Galactic Centre and other galactic nuclei. The massive stars produced will pollute the environment via supernova explosions and potentially produce compact binaries whose mergers may be detectable by the LIGO-VIRGO gravitational waves observatories.

preprint2020arXiv

Merger rates in primordial black hole clusters without initial binaries

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.

preprint2020arXiv

Resonance in the K2-19 system is at odds with its high reported eccentricities

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.

preprint2019arXiv

Resilient habitability of nearby exoplanet systems

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

preprint2019arXiv

The Ecology of the Galactic Centre: Nuclear Stellar Clusters and Supermassive Black Holes

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.

preprint2016arXiv

Gravitational scattering of stars and clusters and the heating of the Galactic disk

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.

preprint2016arXiv

How to form asteroids from mm-sized grains

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.

preprint2016arXiv

Is there an exoplanet in the Solar System?

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.

preprint2016arXiv

Long-term stability of the HR 8799 planetary system without resonant lock

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.

preprint2016arXiv

Survival of habitable planets in unstable planetary systems

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.

preprint2015arXiv

Close encounters involving free-floating planets in star clusters

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.

preprint2015arXiv

How to form planetesimals from mm-sized chondrules and chondrule aggregates

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.

preprint2015arXiv

The destruction of inner planetary systems during high-eccentricity migration of gas giants

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.

preprint2014arXiv

Formation Channels for Blue Straggler Stars

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).

preprint2013arXiv

Misaligned streamers around a galactic centre black hole from a single cloud's infall

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.

preprint2013arXiv

Supernova enrichment and dynamical histories of solar-type stars in clusters

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.

preprint2013arXiv

The Long-Term Dynamical Evolution of Planetary Systems

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.

preprint2012arXiv

An upper limit to the velocity dispersion of relaxed stellar systems without massive black holes

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.

preprint2012arXiv

Can planetary instability explain the Kepler dichotomy?

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.

preprint2012arXiv

Investigating stellar-mass black hole kicks

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.

preprint2012arXiv

The properties of long gamma-ray bursts in massive compact binaries

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.

preprint2011arXiv

Implications for the origin of short gamma-ray bursts from their observed positions around their host galaxies

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.

preprint2011arXiv

Interacting compact binaries: modeling mass transfer in eccentric systems

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.

preprint2011arXiv

Short gamma-ray bursts: evidence for an origin in globular clusters?

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.

preprint2011arXiv

SMBH Formation via Gas Accretion in Nuclear Stellar Clusters

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.

preprint2011arXiv

The origin of very wide binary stars

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.

preprint2010arXiv

A Deep Multiwavelength View of Binaries in Omega Centauri

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.

preprint2010arXiv

On the origin of black hole spin in high-mass black hole binaries: Cygnus X-1

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.

preprint2010arXiv

The effects of fly-bys on planetary systems

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.

preprint2009arXiv

A Chandra Study of the Galactic Globular Cluster Omega Centauri

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.

preprint2003arXiv

High Resolution Calculations of Merging Neutron Stars III: Gamma-Ray Bursts

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 $ν\barν$ and magnetohydrodynamic processes- are explored in order to investigate their viability in launching GRBs. We find that $ν\barν$ 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 $ν\barν$-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.