Showing 1–50 of 63 results for author: Church, R

Gaia dr3 reveals the complex dynamical evolution within star clusters

Authors: C. Viscasillas Vázquez, L. Magrini, N. Miret-Roig, N. J. Wright, J. Alves, L. Spina, R. P. Church, G. Tautvaišienė, S. Randich

Abstract: Context: Star clusters, composed of stars born from the same molecular cloud, serve as invaluable natural laboratories for understanding the fundamental processes governing stellar formation and evolution. Aims: This study aims to investigate correlations between the Mean Interdistance ($\bar{D_{i}}$), Mean Closest Interdistance ($\bar{D_{c}}$) and Median Weighted Central Interdistance (… ▽ More Context: Star clusters, composed of stars born from the same molecular cloud, serve as invaluable natural laboratories for understanding the fundamental processes governing stellar formation and evolution. Aims: This study aims to investigate correlations between the Mean Interdistance ($\bar{D_{i}}$), Mean Closest Interdistance ($\bar{D_{c}}$) and Median Weighted Central Interdistance ($\bar{D_{cc}}$) with the age of star clusters, examining their evolutionary trends and assessing the robustness of these quantities as possible age indicators. Methods: We selected a sample of open clusters in the solar region and with a representative number of members (e.g. well populated and without outliers). The interdistances are derived from the spatial distribution of member stars within a cluster. Their evolution over time allows us to use them as an age indicators for star clusters. Results: Our investigation reveals a high-significant correlation between the interdistances and cluster age. Considering the full sample of clusters between 7 and 9 kpc, the relationship is very broad. This is due to uncertainties in parallax, which increase with increasing distance. In particular, we must limit the sample to a maximum distance from the Sun of about 200 pc to avoid artificial effects on cluster shape and on the spatial distribution of their stars along the line of sight. Conclusions: By conservatively restraining the distance to a maximum of $\sim$200 pc, we have established a relationship between the interdistances and the age of the clusters. In our sample, the relationship is mainly driven by the internal expansion of the clusters, and is marginally affected by external perturbative effects. Such relation might enhance our comprehension of cluster dynamics and might be used to derive cluster dynamical ages. △ Less

Submitted 30 June, 2024; originally announced July 2024.

Comments: This paper was accepted for publication in Astronomy & Astrophysics on June 28, 2024. The document includes 9 main figures and 4 additional figures in the appendix. Additionally, it contains one main table and an extract of another table in the appendix, available at the CDS

The PLATO Mission

Authors: Heike Rauer, Conny Aerts, Juan Cabrera, Magali Deleuil, Anders Erikson, Laurent Gizon, Mariejo Goupil, Ana Heras, Jose Lorenzo-Alvarez, Filippo Marliani, Cesar Martin-Garcia, J. Miguel Mas-Hesse, Laurence O'Rourke, Hugh Osborn, Isabella Pagano, Giampaolo Piotto, Don Pollacco, Roberto Ragazzoni, Gavin Ramsay, Stéphane Udry, Thierry Appourchaux, Willy Benz, Alexis Brandeker, Manuel Güdel, Eduardo Janot-Pacheco , et al. (801 additional authors not shown)

Abstract: PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observati… ▽ More PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases. △ Less

Submitted 8 June, 2024; originally announced June 2024.

A Wide Metallicity Range for Gyr-old Stars in the Nuclear Star Cluster

Authors: Brian Thorsbro, Rebecca Forsberg, Georges Kordopatis, Alessandra Mastrobuono-Battisti, Ross P. Church, R. Michael Rich, Nils Ryde, Mathias Schultheis, Shogo Nishiyama

Abstract: We report metallicities for three $\sim$Gyr-old stars in the Milky Way nuclear star cluster (NSC) using high-resolution near-infrared spectroscopy. We derive effective temperatures from a calibration with Sc line strength, which yields results in good agreement with other methods, and metallicities from spectral fits to Fe I lines. Our derived metallicities range from -1.2 < [Fe/H] < +0.5, a span… ▽ More We report metallicities for three $\sim$Gyr-old stars in the Milky Way nuclear star cluster (NSC) using high-resolution near-infrared spectroscopy. We derive effective temperatures from a calibration with Sc line strength, which yields results in good agreement with other methods, and metallicities from spectral fits to Fe I lines. Our derived metallicities range from -1.2 < [Fe/H] < +0.5, a span of 1.7 dex. In addition we use isochrone projection to obtain masses of 1.6 to 4.3 M$_\odot$, and ages assuming single-star evolution. The oldest of these stars is 1.5 Gyr while the youngest and most metal-rich is only 100 Myr. The wide range in metallicity poses interesting questions concerning the chemical evolution and enrichment of the NSC and adds to the evidence for the presence of a young, metal-rich population in the NSC. We suggest that the candidate intermediate-age, metal-poor ([Fe/H] = -1.2) star may be best explained as a blue straggler from an underlying old population. △ Less

Submitted 30 April, 2024; originally announced April 2024.

Comments: 6 pages, 4 figures

Journal ref: ApJL 958 L18 (2023)

The effects of Galactic model uncertainties on LISA observations of double neutron stars

Authors: Anatole Storck, Ross Church

Abstract: Observations of binaries containing pairs of neutron stars using the upcoming space-based gravitational wave observatory, LISA, have the potential to improve our understanding of neutron star physics and binary evolution. In this work we assess the effect of changing the model of the Milky Way's kinematics and star formation history on predictions of the population of double neutron stars that wil… ▽ More Observations of binaries containing pairs of neutron stars using the upcoming space-based gravitational wave observatory, LISA, have the potential to improve our understanding of neutron star physics and binary evolution. In this work we assess the effect of changing the model of the Milky Way's kinematics and star formation history on predictions of the population of double neutron stars that will be detected and resolved by LISA. We conclude that the spatial distribution of these binaries is insensitive to the choice of galactic models, compared to the stochastic variation induced by the small sample size. In particular, the time-consuming computation of the binaries' Galactic orbits is not necessary. The distributions of eccentricity and gravitational-wave frequency are, however, affected by the choice of star-formation history. Binaries with eccentricities e>0.1, which can be measured by LISA observations, are mostly younger than 100 Myr. We caution that comparisons between different predictions for LISA observations need to use consistent star formation histories, and that the Galactic star formation history should be taken into account in the analysis of the observations themselves. The lack of strong dependency on Galactic models means that LISA detection of double neutron star binaries may provide a relatively clean probe of massive binary star evolution. △ Less

Submitted 27 February, 2023; v1 submitted 3 May, 2022; originally announced May 2022.

Comments: 10 pages, 10 Figures, 4 Tables, accepted for publication in MNRAS

Astrophysics with the Laser Interferometer Space Antenna

Authors: Pau Amaro Seoane, Jeff Andrews, Manuel Arca Sedda, Abbas Askar, Quentin Baghi, Razvan Balasov, Imre Bartos, Simone S. Bavera, Jillian Bellovary, Christopher P. L. Berry, Emanuele Berti, Stefano Bianchi, Laura Blecha, Stephane Blondin, Tamara Bogdanović, Samuel Boissier, Matteo Bonetti, Silvia Bonoli, Elisa Bortolas, Katelyn Breivik, Pedro R. Capelo, Laurentiu Caramete, Federico Cattorini, Maria Charisi, Sylvain Chaty , et al. (134 additional authors not shown)

Abstract: The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery… ▽ More The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe. △ Less

Submitted 25 May, 2023; v1 submitted 11 March, 2022; originally announced March 2022.

Journal ref: Living Reviews in Relativity, Volume 26, Article number: 2 (2023)

Influence of tidal dissipation on outcomes of binary-single encounters between stars and black holes in stellar clusters

Authors: Lucas Hellström, Abbas Askar, Alessandro A. Trani, Mirek Giersz, Ross P. Church, Johan Samsing

Abstract: In the cores of dense stellar clusters, close gravitational encounters between binary and single stars can frequently occur. Using the Tsunami code, we computed the outcome of a large number of binary-single interactions involving two black holes (BHs) and a star to check how the inclusion of orbital energy losses due to tidal dissipation can change the outcome of these chaotic interactions. Each… ▽ More In the cores of dense stellar clusters, close gravitational encounters between binary and single stars can frequently occur. Using the Tsunami code, we computed the outcome of a large number of binary-single interactions involving two black holes (BHs) and a star to check how the inclusion of orbital energy losses due to tidal dissipation can change the outcome of these chaotic interactions. Each interaction was first simulated without any dissipative processes and then we systematically added orbital energy losses due to gravitational wave emission (using post-Newtonian (PN) corrections) and dynamical tides and recomputed the interactions. We find that the inclusion of tides increases the number of BH-star mergers by up to 75 per cent but it does not affect the number of BH-BH mergers. These results highlight the importance of including orbital energy dissipation due to dynamical tides during few-body encounters and evolution of close binary systems within stellar cluster simulations. Consistent with previous studies, we find that the inclusion of PN terms increases the number of BH-BH mergers during binary-single encounters. However, BH-star mergers are largely unaffected by the inclusion of these terms. △ Less

Submitted 29 September, 2022; v1 submitted 28 February, 2022; originally announced March 2022.

Comments: 14 pages, 6 figures, 12 tables, accepted for publication in MNRAS

Formation of supermassive black holes in galactic nuclei II: retention and growth of seed intermediate-mass black holes

Authors: Abbas Askar, Melvyn B. Davies, Ross P. Church

Abstract: In many galactic nuclei, a nuclear stellar cluster (NSC) co-exists with a supermassive black hole (SMBH). In this work, we explore the idea that the NSC forms before the SMBH through the merger of several stellar clusters that may contain intermediate-mass black holes (IMBHs). These IMBHs can subsequently grow by mergers and accretion to form an SMBH. To check the observable consequences of this p… ▽ More In many galactic nuclei, a nuclear stellar cluster (NSC) co-exists with a supermassive black hole (SMBH). In this work, we explore the idea that the NSC forms before the SMBH through the merger of several stellar clusters that may contain intermediate-mass black holes (IMBHs). These IMBHs can subsequently grow by mergers and accretion to form an SMBH. To check the observable consequences of this proposed SMBH seeding mechanism, we created an observationally motivated mock population of galaxies, in which NSCs are constructed by aggregating stellar clusters that may or may not contain IMBHs. We model the growth of IMBHs in the NSCs through gravitational wave (GW) mergers with other IMBHs and gas accretion. In the case of GW mergers, the merged BH can either be retained or ejected depending on the GW recoil kick it receives. The likelihood of retaining the merged BH increases if we consider growth of IMBHs in the NSC through gas accretion. We find that nucleated lower-mass galaxies ($\rm M_{\star} \lesssim 10^{9} \ M_{\odot}$; e.g. M33) have an SMBH seed occupation fraction of about 0.3 to 0.5. This occupation fraction increases with galaxy stellar mass and for more massive galaxies ($\rm 10^{9} \ M_{\odot} \lesssim \rm M_{\star} \lesssim 10^{11} \ M_{\odot}$), it is between 0.5 and 0.8, depending on how BH growth is modelled. These occupation fractions are consistent with observational constraints. Furthermore, allowing for BH growth also allows us to reproduce the observed diversity in the mass range of SMBHs in the $\rm M_{\rm NSC} - M_{\rm BH}$ plane. △ Less

Submitted 20 December, 2021; v1 submitted 22 July, 2021; originally announced July 2021.

Comments: 17 pages, 20 figures, 2 tables, accepted for publication in MNRAS

Transients from ONe White-Dwarf - Neutron-Star/Black-Hole Mergers

Authors: Alexey Bobrick, Yossef Zenati, Hagai B. Perets, Melvyn B. Davies, Ross Church

Abstract: We conduct the first 3D hydrodynamic simulations of oxygen-neon white dwarf-neutron star/black hole mergers (ONe WD-NS/BH mergers). Such mergers constitute a significant fraction, and may even dominate, the inspiral rates of all WD-NS binaries. We post-process our simulations to obtain the nuclear evolution of these systems and couple the results to a supernova spectral synthesis code to obtain th… ▽ More We conduct the first 3D hydrodynamic simulations of oxygen-neon white dwarf-neutron star/black hole mergers (ONe WD-NS/BH mergers). Such mergers constitute a significant fraction, and may even dominate, the inspiral rates of all WD-NS binaries. We post-process our simulations to obtain the nuclear evolution of these systems and couple the results to a supernova spectral synthesis code to obtain the first light curves and spectra for these transients. We find that the amount of $^{56}$Ni synthesised in these mergers grows as a strong function of the WD mass, reaching typically $0.05$ and up to $0.1\,{\rm M}_\odot$ per merger. Photodisintegration leads to similar amounts of $^4$He and about a ten times smaller amount of $^1$H. The nuclear yields from these mergers, in particular those of $^{55}$Mn, may contribute significantly to Galactic chemical evolution. The transients expected from ONe WD-NS mergers are dominantly red/infrared, evolve on month-long timescales and reach bolometric magnitudes of up to -16.5. The current surveys must have already detected these transients or are, alternatively, putting strong constraints on merger scenarios. The properties of the expected transients from WD-NS mergers best agree with faint type Iax supernovae. The Vera Rubin Observatory (LSST) will be detecting up to thousands of merging ONe WD-NS systems per year. We simulate a subset of our models with 2D axisymmetric FLASH code to investigate why they have been challenging for previous studies. We find that the likely main challenge has been effectively modelling the nuclear statistical equilibrium regime in such mergers. △ Less

Submitted 13 January, 2024; v1 submitted 7 April, 2021; originally announced April 2021.

Comments: Accepted version, comments welcome

Close stellar encounters at the Galactic Centre I: The effect on the observed stellar populations

Authors: Alessandra Mastrobuono-Battisti, Ross P. Church, Melvyn B. Davies

Abstract: We model the effects of collisions and close encounters on the stellar populations observed in the Milky Way nuclear stellar cluster (NSC). Our analysis is based on $N$-body simulations in which the NSC forms by accretion of massive stellar clusters around a supermassive black hole. We attach stellar populations to our $N$-body particles and follow the evolution of their stars, and the rate of col… ▽ More We model the effects of collisions and close encounters on the stellar populations observed in the Milky Way nuclear stellar cluster (NSC). Our analysis is based on $N$-body simulations in which the NSC forms by accretion of massive stellar clusters around a supermassive black hole. We attach stellar populations to our $N$-body particles and follow the evolution of their stars, and the rate of collisions and close encounters. The most common encounters are collisions between pairs of main-sequence stars, which lead to mergers: destructive collisions between main-sequence stars and compact objects are rare. We find that the effects of collisions on the stellar populations are small for three reasons. First, our models possess a core which limits the maximum stellar density. Secondly, the velocity dispersion in the NSC is similar to the surface escape velocities of the stars, which minimises the collision rate. Finally, whilst collisions between main-sequence stars destroy bright giants by accelerating their evolution, they also create them by accelerating the evolution of lower-mass stars. These two effects approximately cancel out. We also investigate whether the G2 cloud could be a fuzzball: a compact stellar core which has accreted a tenuous envelope in a close encounter with a red giant. We conclude that fuzzballs with cores below $2\,M_\odot$ have thermal times-scales too short to reproduce G2. A fuzzball with a black-hole core could reproduce the surface properties of G2 but the production rate of such objects in our model is low. △ Less

Submitted 18 May, 2021; v1 submitted 1 April, 2021; originally announced April 2021.

Comments: 17 pages, 15 figures, 4 tables. Accepted for publication in MNRAS

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

Authors: Abbas Askar, Melvyn B. Davies, Ross P. Church

Abstract: 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. W… ▽ More 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. △ Less

Submitted 12 January, 2021; v1 submitted 8 June, 2020; originally announced June 2020.

Comments: 19 pages, 17 figures, 6 tables, accepted for publication in MNRAS

The GALAH survey: Multiple stars and our Galaxy. I. A comprehensive method for deriving properties of FGK binary stars

Authors: G. Traven, S. Feltzing, T. Merle, M. Van der Swaelmen, K. Čotar, R. Church, T. Zwitter, Y. -S. Ting, C. Sahlholdt, M. Asplund, J. Bland-Hawthorn, G. De Silva, K. Freeman, S. Martell, S. Sharma, D. Zucker, S. Buder, A. Casey, V. D'Orazi, J. Kos, G. Lewis, J. Lin, K. Lind, J. Simpson, D. Stello , et al. (2 additional authors not shown)

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

Submitted 30 April, 2020; originally announced May 2020.

Comments: Accepted for publication in A&A

MSC Class: 85-XX

Journal ref: A&A 638, A145 (2020)

Merger rates in primordial black hole clusters without initial binaries

Authors: Valeriya Korol, Ilya Mandel, M. Coleman Miller, Ross P. Church, Melvyn B. Davies

Abstract: 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. Specific… ▽ More 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. △ Less

Submitted 29 June, 2020; v1 submitted 8 November, 2019; originally announced November 2019.

Comments: 8 pages, 4 figures

Journal ref: MNRAS 496, 994-1000 (2020)

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

Authors: Ross P. Church, Alexander J. Mustill, Fan Liu

Abstract: 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 requ… ▽ More 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. △ Less

Submitted 19 August, 2019; originally announced August 2019.

Comments: 12 pages, 10 figures, submitted to MNRAS

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

Authors: Melvyn B. Davies, Abbas Askar, Ross P. Church

Abstract: 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 c… ▽ More 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. △ Less

Submitted 31 July, 2019; originally announced July 2019.

Comments: 4 pages, 3 figures, 1 table. Contribution to the proceedings of the IAU Symposium No. 351: "Star Clusters: from the Milky Way to the Early Universe"

Stellar escapers from M67 can reach solar-like Galactic orbits

Authors: T. G. Jørgensen, R. P. Church

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

Submitted 23 May, 2019; originally announced May 2019.

4MOST Consortium Survey 4: Milky Way Disc and Bulge High-Resolution Survey (4MIDABLE-HR)

Authors: T. Bensby, M. Bergemann, J. Rybizki, B. Lemasle, L. Howes, M. Kovalev, O. Agertz, M. Asplund, P. Barklem, C. Battistini, L. Casagrande, C. Chiappini, R. Church, S. Feltzing, D. Ford, O. Gerhard, I. Kushniruk, G. Kordopatis, K. Lind, I. Minchev, P. McMillan, H. -W. Rix, N. Ryde, G. Traven

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

Submitted 1 April, 2019; v1 submitted 6 March, 2019; originally announced March 2019.

Comments: Part of the 4MOST issue of The Messenger, published in preparation of 4MOST Community Workshop, see http://www.eso.org/sci/meetings/2019/4MOST.html

Journal ref: The Messenger, 2019, 175, 35

4MOST Consortium Survey 3: Milky Way Disc and Bulge Low-Resolution Survey (4MIDABLE-LR)

Authors: C. Chiappini, I. Minchev, E. Starkenburg, F. Anders, N. Gentile Fusillo, O. Gerhard, G. Guiglion, A. Khalatyan, G. Kordopatis, B. Lemasle, G. Matijevic, A. B. de Andrade Queiroz, A. Schwope, M. Steinmetz, J. Storm, G. Traven, P. -E. Tremblay, M. Valentini, R. Andrae, A. Arentsen, M. Asplund, T. Bensby, M. Bergemann, L. Casagrande, R. Church , et al. (13 additional authors not shown)

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

Submitted 1 April, 2019; v1 submitted 6 March, 2019; originally announced March 2019.

Comments: Part of the 4MOST issue of The Messenger, published in preparation of 4MOST Community Workshop, see http://www.eso.org/sci/meetings/2019/4MOST.html

Journal ref: The Messenger, 2019, 175, 30

4MOST Consortium Survey 2: The Milky Way Halo High-Resolution Survey

Authors: N. Christlieb, C. Battistini, P. Bonifacio, E. Caffau, H. -G. Ludwig, M. Asplund, P. Barklem, M. Bergemann, R. Church, S. Feltzing, D. Ford, E. K. Grebel, C. J. Hansen, A. Helmi, G. Kordopatis, M. Kovalev, A. Korn, K. Lind, A. Quirrenbach, J. Rybizki, Á. Skúladóttir, E. Starkenburg

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

Submitted 1 April, 2019; v1 submitted 6 March, 2019; originally announced March 2019.

Comments: Part of the 4MOST issue of The Messenger, published in preparation of 4MOST Community Workshop, see http://www.eso.org/sci/meetings/2019/4MOST.html

Journal ref: The Messenger, 2019, 175, 26

4MOST Survey Strategy Plan

Authors: G. Guiglion, C. Battistini, C. P. M. Bell, T. Bensby, T. Boller, C. Chiappini, J. Comparat, N. Christlieb, R. Church, M. -R. L. Cioni, L. Davies, T. Dwelly, R. S. de Jong, S. Feltzing, A. Gueguen, L. Howes, M. Irwin, I. Kushniruk, M. I Lam, J. Liske, R. McMahon, A. Merloni, P. Norberg, A. S. G. Robotham, O. Schnurr , et al. (8 additional authors not shown)

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

Submitted 1 April, 2019; v1 submitted 6 March, 2019; originally announced March 2019.

Comments: Part of the 4MOST issue of The Messenger, published in preparation of 4MOST Community Workshop, see http://www.eso.org/sci/meetings/2019/4MOST.html

Journal ref: The Messenger, 2019, 175, 17

4MOST Scientific Operations

Authors: C. J. Walcher, M. Banerji, C. Battistini, C. P. M. Bell, O. Bellido-Tirado, T. Bensby, J. M. Bestenlehner, T. Boller, J. Brynnel, A. Casey, C. Chiappini, N. Christlieb, R. Church, M. -R. L. Cioni, S. Croom, J. Comparat, L. J. M. Davies, R. S. de Jong, T. Dwelly, H. Enke, S. Feltzing, D. Feuillet, M. Fouesneau, D. Ford, S. Frey , et al. (43 additional authors not shown)

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

Submitted 1 April, 2019; v1 submitted 6 March, 2019; originally announced March 2019.

Comments: Part of the 4MOST issue of The Messenger, published in preparation of 4MOST Community Workshop, see http://www.eso.org/sci/meetings/2019/4MOST.html

Journal ref: The Messenger, 2019, 175, 12

4MOST: Project overview and information for the First Call for Proposals

Authors: R. S. de Jong, O. Agertz, A. Agudo Berbel, J. Aird, D. A. Alexander, A. Amarsi, F. Anders, R. Andrae, B. Ansarinejad, W. Ansorge, P. Antilogus, H. Anwand-Heerwart, A. Arentsen, A. Arnadottir, M. Asplund, M. Auger, N. Azais, D. Baade, G. Baker, S. Baker, E. Balbinot, I. K. Baldry, M. Banerji, S. Barden, P. Barklem , et al. (313 additional authors not shown)

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

Submitted 1 April, 2019; v1 submitted 6 March, 2019; originally announced March 2019.

Comments: Part of the 4MOST issue of The Messenger, published in preparation of 4MOST Community Workshop, see http://www.eso.org/sci/meetings/2019/4MOST.html

Journal ref: The Messenger, 2019, 175, 3

Rapid destruction of protoplanetary discs due to externalphotoevaporation in star-forming regions

Authors: Rhana B. Nicholson, Richard J. Parker, Ross P. Church, Melvyn B. Davies, Niamh M. Fearon, Sam R. J. Walton

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

Submitted 28 February, 2019; originally announced February 2019.

Comments: 14 pages, 5 figures

Benchmark ages for the Gaia benchmark stars

Authors: Christian L. Sahlholdt, Sofia Feltzing, Lennart Lindegren, Ross P. Church

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

Submitted 5 October, 2018; originally announced October 2018.

Comments: Accepted in MNRAS. 69 pages (18 for main text, 11 for appendix, and 40 for extra figures)

Thermohaline mixing in extremely metal-poor stars

Authors: Kate Henkel, Amanda I. Karakas, Andrew R. Casey, Ross P. Church, John C. Lattanzio

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

Submitted 6 August, 2018; originally announced August 2018.

Comments: 8 pages, 4 figures, 1 table

Journal ref: The Astrophysical Journal Letters, 2018, Volume 863, Number 1

Gaia DR2 confirms that candidate Thorne-Żytkow object HV 2112 is in the SMC

Authors: Paul J. McMillan, Ross P. Church

Abstract: The Gaia DR2 proper motion places HV 2112 in the Small Magellanic Cloud (SMC), consistent with it being a Thorne-Żytkow object. The Gaia DR2 proper motion places HV 2112 in the Small Magellanic Cloud (SMC), consistent with it being a Thorne-Żytkow object. △ Less

Submitted 26 April, 2018; originally announced April 2018.

Comments: 2 pages, submitted to RNAAS

Stellar multiplicity in high-resolution spectroscopic surveys. I. Application to APOGEE subgiants and giants

Authors: Edita Stonkutė, Ross P. Church, Sofia Feltzing, Jennifer A. Johnson

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

Submitted 24 April, 2018; originally announced April 2018.

Comments: 13 pages, 10 figures, submitted to MNRAS

Estimating stellar ages and metallicities from parallaxes and broadband photometry - successes and shortcomings

Authors: Louise M. Howes, Lennart Lindegren, Sofia Feltzing, Ross P. Church, Thomas Bensby

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

Submitted 16 November, 2018; v1 submitted 23 April, 2018; originally announced April 2018.

Comments: Accepted for publication in A&A. 22 pages, 19 figures (+ 6 pages of appendix). Abstract abridged for the arXiv

Journal ref: A&A 622, A27 (2019)

Formation constraints indicate a black-hole accretor in 47 Tuc X9

Authors: Ross P. Church, Jay Strader, Melvyn B. Davies, Alexey Bobrick

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

Submitted 2 January, 2018; originally announced January 2018.

Comments: 6 pages, 3 figures, published in ApJL

Journal ref: ApJL 851 (2017) 4

Stellar multiplicity in the Milky Way Galaxy

Authors: Edita Stonkute, Ross P. Church, Sofia Feltzing, Jennifer A. Johnson

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

Submitted 21 August, 2017; originally announced August 2017.

Comments: 2 pages, 1 figure, to appear in the proceedings of the IAU Symposium 334 "Rediscovering our Galaxy", Potsdam, 10-14 July 2017, eds. C. Chiappini, I. Minchev, E. Starkenburg, M. Valentini

Mass transfer in white dwarf-neutron star binaries

Authors: Alexey Bobrick, Melvyn B. Davies, Ross P. Church

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

Submitted 8 February, 2017; originally announced February 2017.

Comments: Accepted to MNRAS. For videos from the runs, please see http://www.astro.lu.se/~alexey/animations.html

GaiaNIR: Combining optical and Near-Infra-Red (NIR) capabilities with Time-Delay-Integration (TDI) sensors for a future Gaia-like mission

Authors: David Hobbs, Erik Høg, Alcione Mora, Cian Crowley, Paul McMillan, Piero Ranalli, Ulrike Heiter, Carme Jordi, Nigel Hambly, Ross Church, Brown Anthony, Paolo Tanga, Laurent Chemin, Jordi Portell, Fran Jiménez-Esteban, Sergei Klioner, Francois Mignard, Johan Fynbo, Łukasz Wyrzykowski, Krzysztof Rybicki, Richard I. Anderson, Alberto Cellino, Claus Fabricius, Michael Davidson, Lennart Lindegren

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

Submitted 22 May, 2020; v1 submitted 23 September, 2016; originally announced September 2016.

Comments: 27 Pages

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

Authors: Ylva Gotberg, Melvyn B. Davies, Alexander Mustill, Anders Johansen, Ross P. Church

Abstract: 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 prevent… ▽ More 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. △ Less

Submitted 24 June, 2016; originally announced June 2016.

Comments: Accepted for publication in A&A

Journal ref: A&A 592, A147 (2016)

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

Authors: Bengt Gustafsson, Ross Church, Melvyn B. Davies, Hans Rickman

Abstract: 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 an… ▽ More 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. △ Less

Submitted 23 May, 2016; v1 submitted 10 May, 2016; originally announced May 2016.

Comments: 19 pages, 15 figures

Report number: NORDITA-2016-36

Journal ref: A&A 593, A85 (2016)

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

Authors: J. D. Lyman, A. J. Levan, P. A. James, C. R. Angus, R. P. Church, M. B. Davies, N. R. Tanvir

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

Submitted 25 February, 2016; originally announced February 2016.

Comments: 11 pages, accepted for publication in MNRAS

Diagnostics of Stellar Modelling from Spectroscopy and Photometry of Globular Clusters

Authors: George Angelou, Valentina D'Orazi, Thomas Constantino, Ross Church, Richard Stancliffe, John Lattanzio

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

Submitted 7 April, 2015; originally announced April 2015.

Comments: 20 pages, 11 Figures, 4 Tables, Accepted for publication in MNRAS

Close encounters involving free-floating planets in star clusters

Authors: Long Wang, M. B. N. Kouwenhoven, Xiaochen Zheng, Ross P. Church, Melvyn B. Davies

Abstract: 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 t… ▽ More 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. △ Less

Submitted 10 March, 2015; originally announced March 2015.

Comments: Accepted for publication in MNRAS; 18 pages, 12 figures

Journal ref: MNRAS 449, 3543-3558 (2015)

On the Numerical Treatment and Dependence of Thermohaline Mixing in Red Giants

Authors: John Lattanzio, Lionel Siess, Ross Church, George Angelou, Richard Stancliffe, Carolyn Doherty, Thomas Stephen, Simon Campbell

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

Submitted 2 November, 2014; v1 submitted 23 October, 2014; originally announced October 2014.

Comments: 19 pages, 18 figures, accepted for publication by MNRAS; revised version has corrected Fig 17

The progenitors of calcium-rich transients are not formed in situ

Authors: Joe Lyman, Andrew Levan, Ross Church, Melvyn Davies, Nial Tanvir

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

Submitted 6 August, 2014; originally announced August 2014.

Comments: Accepted for publication in MNRAS, 11 pages, 9 figures

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

Authors: Christopher A. Tout, Anna N. Zytkow, Ross P. Church, Herbert H. B. Lau

Abstract: 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 st… ▽ More 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. △ Less

Submitted 23 June, 2014; originally announced June 2014.

Comments: 5 pages, submitted to MNRAS letters. Comments welcome

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

Authors: Richard J. Parker, Ross P. Church, Melvyn B. Davies, Michael R. Meyer

Abstract: 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… ▽ More 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. △ Less

Submitted 11 October, 2013; originally announced October 2013.

Comments: 14 pages, 5 figures, accepted for publication in MNRAS

The PLATO 2.0 Mission

Authors: H. Rauer, C. Catala, C. Aerts, T. Appourchaux, W. Benz, A. Brandeker, J. Christensen-Dalsgaard, M. Deleuil, L. Gizon, M. -J. Goupil, M. Güdel, E. Janot-Pacheco, M. Mas-Hesse, I. Pagano, G. Piotto, D. Pollacco, N. C. Santos, A. Smith, J. -C., Suárez, R. Szabó, S. Udry, V. Adibekyan, Y. Alibert, J. -M. Almenara , et al. (137 additional authors not shown)

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

Submitted 4 March, 2014; v1 submitted 2 October, 2013; originally announced October 2013.

Comments: 63 pages, 17 figures, submitted to Experimental Astronomy (ExA)

Can planetary instability explain the Kepler dichotomy?

Authors: Anders Johansen, Melvyn B. Davies, Ross P. Church, Viktor Holmelin

Abstract: 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-tran… ▽ More 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. △ Less

Submitted 18 August, 2012; v1 submitted 28 June, 2012; originally announced June 2012.

Comments: Accepted for publication in ApJ

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

Authors: Ross P. Church, Chunglee Kim, Andrew J. Levan, Melvyn B. Davies

Abstract: 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… ▽ More 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. △ Less

Submitted 3 July, 2012; v1 submitted 2 May, 2012; originally announced May 2012.

Comments: 7 pages, 6 figures. Accepted by MNRAS

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

Authors: George C. Angelou, Richard J. Stancliffe, Ross P. Church, John C. Lattanzio, Graeme H. Smith

Abstract: 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 u… ▽ More 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. △ Less

Submitted 13 February, 2012; originally announced February 2012.

Comments: 13 Pages, 3 figures. Accepted for publication in the Astrophysical Journal

Interacting compact binaries: modeling mass transfer in eccentric systems

Authors: Ross P. Church, Melvyn B. Davies, Alexey Bobrick, Christopher A. Tout

Abstract: 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 under… ▽ More 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. △ Less

Submitted 19 October, 2011; originally announced October 2011.

Comments: 4 pages, 2 figures, to appear in the ASP conference series proceedings of "Advances in Computational Astrophysics: methods, tools, and outcomes" in Cefalu', Italy, June 13-17, 2011

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

Authors: Ross P. Church, Andrew J. Levan, Melvyn B. Davies, Nial Tanvir

Abstract: 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 0605… ▽ More 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. △ Less

Submitted 19 October, 2011; originally announced October 2011.

Comments: 4 pages, 4 figures. To appear in the proceedings of "GRBs as probes: from the progenitor's environment to the high redshift Universe", Como, Italy, May 2011

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

Authors: Ross P. Church, Jennifer A. Johnson, Sofia Feltzing

Abstract: 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 2MAS… ▽ More 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. △ Less

Submitted 28 March, 2011; originally announced March 2011.

Comments: 30 pages, accepted by A&A as a Research Note

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

Authors: Ross P. Church, Andrew J. Levan, Melvyn B. Davies, Nial Tanvir

Abstract: 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 circu… ▽ More 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. △ Less

Submitted 5 January, 2011; originally announced January 2011.

Comments: 12 pages, 13 figures. Accepted for publication in MNRAS

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

Authors: George C. Angelou, Ross P. Church, Richard J. Stancliffe, John C. Lattanzio, Graeme H. Smith

Abstract: 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 paramet… ▽ More 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. △ Less

Submitted 9 December, 2010; originally announced December 2010.

Comments: 33 Pages, 10 Figures. Accepted for publication in The Astrophysical Journal

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

Authors: Magnus Axelsson, Ross P. Church, Melvyn B. Davies, Andrew J. Levan, Felix Ryde

Abstract: 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, t… ▽ More 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. △ Less

Submitted 19 November, 2010; originally announced November 2010.

Comments: 5 pages, 2 tables, 1 figure. Accepted for publication in MNRAS