Source author record

Ross P. Church

Ross P. Church 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

20works
6topics
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

20 published item(s)

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.

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.

preprint2019arXiv

Super-Earth ingestion can explain the anomalously high metal abundances of M67 Y2235

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

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

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.

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.

preprint2014arXiv

HV2112, a Thorne-Zytkow Object or a Super Asymptotic Giant Branch Star

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

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.

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

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.

preprint2012arXiv

The Role of Thermohaline Mixing in Intermediate- and Low-Metallicity Globular Clusters

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

preprint2011arXiv

Coordinates and 2MASS and OGLE identifications for all stars in Arp's 1965 finding chart for Baade's Window

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

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.

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

Thermohaline Mixing and its Role in the Evolution of Carbon and Nitrogen Abundances in Globular Cluster Red Giants: The Test Case of Messier 3

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

preprint2010arXiv

White-Dwarf Kicks and Implications for Barium Stars

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

preprint2009arXiv

The depletion of carbon by extra mixing in metal-poor giants

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

preprint2009arXiv

Why do low-mass stars become red giants?

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