Ross P. Church

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, varying the internal inclination variation of the orbits. 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 5 degrees. Our analysis shows that the Kepler data contain a population of planets larger than four Earth radii in single-transit systems that cannot 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 by planet-planet collisions, and less frequently by ejections or collisions with the star, yielding transits in agreement with the large planets in the Kepler single-transit systems. The resulting population of massive double-planet systems nevertheless cannot 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.

Red giant stellar collisions in the Galactic Centre

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. We show that the missing objects are low-mass giants and asymptotic giant branch stars in the range 1-3 M-circle dot. 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 ROLE OF THERMOHALINE MIXING IN INTERMEDIATE- AND LOW-METALLICITY GLOBULAR CLUSTERS

It is now widely accepted that globular cluster red giant branch (RGB) 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 imprints 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 RGBs 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.

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 3He. We present the evolution of [C/Fe], [N/Fe] and 12C/13C for three models: a carbon-normal metal-poor star, and two stars that have accreted material from a 1.5 M⊙ AGB companion, one having received 0.01 M⊙ of material and the other having received 0.1 M⊙. 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.

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 4000 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. (Less)

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.

On the numerical treatment and dependence of thermohaline mixing in red giants

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

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 characterize 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 per cent 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 per cent 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.

Hubble Space Telescope observations of the host galaxies and environments of calcium-rich supernovae

Calcium-rich supernovae (SNe) represent a significant challenge for our understanding of the fates of stellar systems. They are less luminous than other 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 (~1/3 are offset >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 five members of the class three exhibiting large offsets and two coincident with the disc 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.

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.