Source author record

M. B. Davies

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

29works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

29 published item(s)

preprint2025arXiv

Improved characterization of the TOI-2141 system: a dense sub-Neptune with non-transiting inner and outer companions

We aim to refine the fundamental parameters of the TOI-2141 planetary system, which includes a transiting sub-Neptune orbiting a Sun-like star in a relatively long orbit of 18.26 days, by combining new photometric and spectroscopic observations. We analyze new space-based photometry from TESS and CHEOPS as well as 61 radial velocity measurements from HARPS-N. We perform individual and joint photometric and RV analyses using several modeling tools within a Bayesian model comparison framework. We refine the radius and mass of the transiting planet TOI-2141 b to 3.15 $\pm$ 0.04 $R_\oplus$ and 20.1 $\pm$ 1.6 $M_\oplus$, respectively, five and two times more precise than the previously reported values. Our radial velocity analysis reveals two additional non-transiting companions with orbital periods of 5.46 and 60.45 days. Despite the innermost planet's high geometric transit probability, we find no evidence for transits in the photometric data. The bulk properties of TOI-2141 b suggest a significant volatile envelope atop an Earth-like core, with modeling indicating a hydrogen-rich atmosphere that may have experienced mild photoevaporation over the system's history. Planets b and c must exhibit a modest mutual inclination of at least 2.4 degrees.

preprint2022arXiv

55 Cancri e's occultation captured with CHEOPS

Past occultation and phase-curve observations of the ultra-short period super-Earth 55 Cnc e obtained at visible and infrared wavelengths have been challenging to reconcile with a planetary reflection and emission model. In this study, we analyse a set of 41 occultations obtained over a two-year timespan with the CHEOPS satellite. We report the detection of 55 Cnc e's occultation with an average depth of $12\pm3$ ppm. We derive a corresponding 2-$σ$ upper limit on the geometric albedo of $A_g < 0.55$ once decontaminated from the thermal emission measured by Spitzer at 4.5$μ$m. CHEOPS's photometric performance enables, for the first time, the detection of individual occultations of this super-Earth in the visible and identifies short-timescale photometric corrugations likely induced by stellar granulation. We also find a clear 47.3-day sinusoidal pattern in the time-dependent occultation depths that we are unable to relate to stellar noise, nor instrumental systematics, but whose planetary origin could be tested with upcoming JWST occultation observations of this iconic super-Earth.

preprint2022arXiv

A stellar occultation by the transneptunian object (50000) Quaoar observed by CHEOPS

Stellar occultation is a powerful technique that allows the determination of some physical parameters of the occulting object. The result depends on the photometric accuracy, the temporal resolution, and the number of chords obtained. Space telescopes can achieve high photometric accuracy as they are not affected by atmospheric scintillation. Using ESA's CHEOPS space telescope, we observed a stellar occultation by the Transneptunian object (50000) Quaoar. We compare the obtained chord with previous occultations by this object and determine its astrometry with sub-milliarcsecond precision. Also, we determine upper limits to the presence of a global methane atmosphere on the occulting body. We predicted and observed a stellar occultation by Quaoar using the CHEOPS space telescope. We measured the occultation light curve from this data-set and determined the dis- and re-appearance of the star behind the occulting body. Furthermore, a ground-based telescope in Australia was used to constrain Quaoar's limb. Combined with results from previous works, these measurements allow us to obtain a precise position of Quaoar at the occultation time. We present results obtained from the first stellar occultation by a Transneptunian object (TNO) using space telescope orbiting Earth. It was the occultation by Quaoar observed on 2020 June 11. We used the CHEOPS light curve to obtain a surface pressure upper limit of 85 nbar for the detection of a global methane atmosphere. Also, combining this observation with a ground-based observation we fit Quaoar's limb to determine its astrometric position with an uncertainty below 1.0 mas. This observation is a first of its kind, and it shall be considered as a proof of concept of stellar occultation observations of Transneptunian objects with space telescopes orbiting Earth. Moreover, it shows significant prospects for the James Webb Space Telescope.

preprint2022arXiv

Analysis of Early Science observations with the CHaracterising ExOPlanets Satellite (CHEOPS) using pycheops

CHEOPS(CHaracterising ExOPlanet Satellite) is an ESA S-class mission that observes bright stars at high cadence from low-Earth orbit. The main aim of the mission is to characterize exoplanets that transit nearby stars using ultrahigh precision photometry. Here we report the analysis of transits observed by CHEOPS during its Early Science observing programme for four well-known exoplanets: GJ436b, HD106315b, HD97658b and GJ1132b. The analysis is done using pycheops, an open-source software package we have developed to easily and efficiently analyse CHEOPS light curve data using state-of-the-art techniques that are fully described herein. We show that the precision of the transit parameters measured using CHEOPS is comparable to that from larger space telescopes such as Spitzer Space Telescope and Kepler. We use the updated planet parameters from our analysis to derive new constraints on the internal structure of these four exoplanets.

preprint2022arXiv

CHEOPS geometric albedo of the hot Jupiter HD 209458b

We report the detection of the secondary eclipse of the hot Jupiter HD 209458b in optical/visible light using the CHEOPS space telescope. Our measurement of 20.4 +/- 3.3 ppm translates into a geometric albedo of A_g = 0.096 +/- 0.016. The previously estimated dayside temperature of about 1500 K implies that our geometric albedo measurement consists predominantly of reflected starlight and is largely uncontaminated by thermal emission. This makes the present result one of the most robust measurements of A_g for any exoplanet. Our calculations of the bandpass-integrated geometric albedo demonstrate that the measured value of A_g is consistent with a cloud-free atmosphere, where starlight is reflected via Rayleigh scattering by hydrogen molecules, and the water and sodium abundances are consistent with stellar metallicity. We predict that the bandpass-integrated TESS geometric albedo is too faint to detect and that a phase curve of HD 209458b observed by CHEOPS would have a distinct shape associated with Rayleigh scattering if the atmosphere is indeed cloud free.

preprint2022arXiv

Fundamental effective temperature measurements for eclipsing binary stars -- III. SPIRou near-infrared spectroscopy and CHEOPS photometry of the benchmark G0V star EBLM J0113+31

EBLM J0113+31 is moderately bright (V=10.1), metal-poor ([Fe/H]$\approx-0.3$) G0V star with a much fainter M dwarf companion on a wide, eccentric orbit (=14.3 d). We have used near-infrared spectroscopy obtained with the SPIRou spectrograph to measure the semi-amplitude of the M dwarf's spectroscopic orbit, and high-precision photometry of the eclipse and transit from the CHEOPS and TESS space missions to measure the geometry of this binary system. From the combined analysis of these data together with previously published observations we obtain the following model-independent masses and radii: $M_1 = 1.029 \pm 0.025 M_{\odot}$, $M_2 = 0.197 \pm 0.003 M_{\odot}$, $R_1 = 1.417 \pm 0.014 R_{\odot}$, $R_2 = 0.215 \pm 0.002 R_{\odot}$. Using $R_1$ and the parallax from Gaia EDR3 we find that this star's angular diameter is $θ= 0.0745 \pm 0.0007$ mas. The apparent bolometric flux of the G0V star corrected for both extinction and the contribution from the M dwarf ($<0.2$ per cent) is ${\mathcal F}_{\oplus,0} = (2.62\pm 0.05)\times10^{-9}$ erg.cm$^{-2}$.s$^{-1}$. Hence, this G0V star has an effective temperature $T_{\rm eff,1} = 6124{\rm\,K} \pm 40{\rm \,K\,(rnd.)} \pm 10 {\rm \,K\,(sys.)}$. EBLM J0113+31 is an ideal benchmark star that can be used for "end-to-end" tests of the stellar parameters measured by large-scale spectroscopic surveys, or stellar parameters derived from asteroseismology with PLATO. The techniques developed here can be applied to many other eclipsing binaries in order to create a network of such benchmark stars.

preprint2022arXiv

Investigating the architecture and internal structure of the TOI-561 system planets with CHEOPS, HARPS-N and TESS

We present a precise characterization of the TOI-561 planetary system obtained by combining previously published data with TESS and CHEOPS photometry, and a new set of $62$ HARPS-N radial velocities (RVs). Our joint analysis confirms the presence of four transiting planets, namely TOI-561 b ($P = 0.45$ d, $R = 1.42$ R$_\oplus$, $M = 2.0$ M$_\oplus$), c ($P = 10.78$ d, $R = 2.91$ R$_\oplus$, $M = 5.4$ M$_\oplus$), d ($P = 25.7$ d, $R = 2.82$ R$_\oplus$, $M = 13.2$ M$_\oplus$) and e ($P = 77$ d, $R = 2.55$ R$_\oplus$, $M = 12.6$ M$_\oplus$). Moreover, we identify an additional, long-period signal ($>450$ d) in the RVs, which could be due to either an external planetary companion or to stellar magnetic activity. The precise masses and radii obtained for the four planets allowed us to conduct interior structure and atmospheric escape modelling. TOI-561 b is confirmed to be the lowest density ($ρ_{\rm b} = 3.8 \pm 0.5$ g cm$^{-3}$) ultra-short period (USP) planet known to date, and the low metallicity of the host star makes it consistent with the general bulk density-stellar metallicity trend. According to our interior structure modelling, planet b has basically no gas envelope, and it could host a certain amount of water. In contrast, TOI-561 c, d, and e likely retained an H/He envelope, in addition to a possibly large water layer. The inferred planetary compositions suggest different atmospheric evolutionary paths, with planets b and c having experienced significant gas loss, and planets d and e showing an atmospheric content consistent with the original one. The uniqueness of the USP planet, the presence of the long-period planet TOI-561 e, and the complex architecture make this system an appealing target for follow-up studies.

preprint2022arXiv

Spi-OPS: Spitzer and CHEOPS confirm the near-polar orbit of MASCARA-1 b and reveal a hint of dayside reflection

The light curves of tidally locked hot Jupiters transiting fast-rotating, early-type stars are a rich source of information about both the planet and star, with full-phase coverage enabling a detailed atmospheric characterisation of the planet. Although it is possible to determine the true spin-orbit angle $Ψ$, a notoriously difficult parameter to measure, from any transit asymmetry resulting from gravity darkening induced by the stellar rotation, the correlations that exist between the transit parameters have led to large disagreements in published values of $Ψ$ for some systems. We aimed to study these phenomena in the light curves of the ultra-hot Jupiter MASCARA-1 b. We obtained optical CHEOPS transit and occultation light curves of MASCARA-1 b, and analysed them jointly with a Spitzer/IRAC 4.5 $μ$m full-phase curve. When fitting the CHEOPS and Spitzer transits together, the degeneracies are greatly diminished and return results consistent with previously published Doppler tomography. Placing priors informed by the tomography achieves even better precision, allowing a determination of $Ψ=72.1^{+2.5}_{-2.4}$ deg. From the occultations and phase variations, we derived dayside and nightside temperatures of $3062^{+66}_{-68}$ K and $1720\pm330$ K, respectively. In addition, we could separately derive geometric albedo $A_g=0.171^{+0.066}_{-0.068}$ and spherical albedo $A_s=0.266^{+0.097}_{-0.100}$ from the CHEOPS data, and Bond albedo $A_B=0.057^{+0.083}_{-0.101}$ from the Spitzer phase curve. Where possible, priors informed by Doppler tomography should be used when fitting transits of fast-rotating stars, though multi-colour photometry may also unlock an accurate measurement of $Ψ$. Our approach to modelling the phase variations at different wavelengths provides a template for how to separate thermal emission from reflected light in spectrally resolved JWST phase curves.

preprint2022arXiv

The atmosphere and architecture of WASP-189 b probed by its CHEOPS phase curve

Gas giants orbiting close to hot and massive early-type stars can reach dayside temperatures that are comparable to those of the coldest stars. These "ultra-hot Jupiters" have atmospheres made of ions and atomic species from molecular dissociation and feature strong day-to-night temperature gradients. Photometric observations at different orbital phases provide insights on the planet atmospheric properties. We analyse the photometric observations of WASP-189 acquired with the instrument CHEOPS to derive constraints on the system architecture and the planetary atmosphere. We implement a light curve model suited for asymmetric transit shape caused by the gravity-darkened photosphere of the fast-rotating host star. We also model the reflective and thermal components of the planetary flux, the effect of stellar oblateness and light-travel time on transit-eclipse timings, the stellar activity and CHEOPS systematics. From the asymmetric transit, we measure the size of the ultra-hot Jupiter WASP-189 b, $R_p=1.600^{+0.017}_{-0.016}\,R_J$, with a precision of 1%, and the true orbital obliquity of the planetary system $Ψ_p=89.6\pm1.2°$ (polar orbit). We detect no significant hotspot offset from the phase curve and obtain an eclipse depth $δ_\text{ecl}=96.5^{+4.5}_{-5.0}\,\text{ppm}$, from which we derive an upper limit on the geometric albedo: $A_g<0.48$. We also find that the eclipse depth can only be explained by thermal emission alone in the case of extremely inefficient energy redistribution. Finally, we attribute the photometric variability to the stellar rotation, either through superficial inhomogeneities or resonance couplings between the convective core and the radiative envelope.

preprint2022arXiv

Transit timing variations of AU Microscopii b and c

Here we report large-amplitude transit timing variations (TTVs) for AU\,Microcopii b and c as detected in combined TESS (2018, 2020) and CHEOPS (2020, 2021) transit observations. AU Mic is a young planetary system with a debris disk and two transiting warm Neptunes. A TTV on the order of several minutes was previously reported for AU Mic b, which was suggested to be an outcome of mutual perturbations between the planets in the system. In 2021, we observed AU Mic b (five transits) and c (three transits) with the CHEOPS space telescope to follow-up the TTV of AU Mic b and possibly detect a TTV for AU Mic c. When analyzing TESS and CHEOPS 2020--2021 measurements together, we find that a prominent TTV emerges with a full span of ~23 minutes between the two TTV extrema. Assuming that the period change results from a periodic process - such as mutual perturbations - we demonstrate that the times of transits in the summer of 2022 are expected to be 30--85 minutes later than predicted by the available linear ephemeris.

preprint2021arXiv

CHEOPS observations of the HD 108236 planetary system: A fifth planet, improved ephemerides, and planetary radii

The detection of a super-Earth and three mini-Neptunes transiting the bright ($V$ = 9.2 mag) star HD 108236 (also known as TOI-1233) was recently reported on the basis of TESS and ground-based light curves. We perform a first characterisation of the HD 108236 planetary system through high-precision CHEOPS photometry and improve the transit ephemerides and system parameters. We characterise the host star through spectroscopic analysis and derive the radius with the infrared flux method. We constrain the stellar mass and age by combining the results obtained from two sets of stellar evolutionary tracks. We analyse the available TESS light curves and one CHEOPS transit light curve for each known planet in the system. We find that HD 108236 is a Sun-like star with $R_{\star}=0.877\pm0.008 R_{\odot}$, $M_{\star}=0.869^{+0.050}_{-0.048} M_{\odot}$, and an age of $6.7_{-5.1}^{+4.0}$ Gyr. We report the serendipitous detection of an additional planet, HD 108236 f, in one of the CHEOPS light curves. For this planet, the combined analysis of the TESS and CHEOPS light curves leads to a tentative orbital period of about 29.5 days. From the light curve analysis, we obtain radii of $1.615\pm0.051$, $2.071\pm0.052$, $2.539_{-0.065}^{+0.062}$, $3.083\pm0.052$, and $2.017_{-0.057}^{+0.052}$ $R_{\oplus}$ for planets HD 108236 b to HD 108236 f, respectively. These values are in agreement with previous TESS-based estimates, but with an improved precision of about a factor of two. We perform a stability analysis of the system, concluding that the planetary orbits most likely have eccentricities smaller than 0.1. We also employ a planetary atmospheric evolution framework to constrain the masses of the five planets, concluding that HD 108236 b and HD 108236 c should have an Earth-like density, while the outer planets should host a low mean molecular weight envelope.

preprint2021arXiv

Six transiting planets and a chain of Laplace resonances in TOI-178

Determining the architecture of multi-planetary systems is one of the cornerstones of understanding planet formation and evolution. Resonant systems are especially important as the fragility of their orbital configuration ensures that no significant scattering or collisional event has taken place since the earliest formation phase when the parent protoplanetary disc was still present. In this context, TOI-178 has been the subject of particular attention since the first TESS observations hinted at a 2:3:3 resonant chain. Here we report the results of observations from CHEOPS, ESPRESSO, NGTS, and SPECULOOS with the aim of deciphering the peculiar orbital architecture of the system. We show that TOI-178 harbours at least six planets in the super-Earth to mini-Neptune regimes, with radii ranging from 1.152(-0.070/+0.073) to 2.87(-0.13/+0.14) Earth radii and periods of 1.91, 3.24, 6.56, 9.96, 15.23, and 20.71 days. All planets but the innermost one form a 2:4:6:9:12 chain of Laplace resonances, and the planetary densities show important variations from planet to planet, jumping from 1.02(+0.28/-0.23) to 0.177(+0.055/-0.061) times the Earth's density between planets c and d. Using Bayesian interior structure retrieval models, we show that the amount of gas in the planets does not vary in a monotonous way, contrary to what one would expect from simple formation and evolution models and unlike other known systems in a chain of Laplace resonances. The brightness of TOI-178 allows for a precise characterisation of its orbital architecture as well as of the physical nature of the six presently known transiting planets it harbours. The peculiar orbital configuration and the diversity in average density among the planets in the system will enable the study of interior planetary structures and atmospheric evolution, providing important clues on the formation of super-Earths and mini-Neptunes.

preprint2019arXiv

The Kepler-11 system: evolution of the stellar high-energy emission and {initial planetary} atmospheric mass fractions

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

preprint2016arXiv

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

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.

preprint2015arXiv

The formation of the solar system

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

preprint2014arXiv

The PLATO 2.0 Mission

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.

preprint2012arXiv

The structure of star clusters in the outer halo of M31

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

preprint2011arXiv

The ACS Survey of Galactic Globular Clusters. XII. Photometric Binaries along the Main-Sequence

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

preprint2010arXiv

The formation of very wide binaries during the star cluster dissolution phase

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

preprint2010arXiv

The Impact of Stellar Collisions in the Galactic Center

We consider whether stellar collisions can explain the observed depletion of red giants in the Galactic center. We model the stellar population with two different IMFs: 1) the Miller-Scalo and 2) a much flatter IMF. In the former case, low-mass main-sequence stars dominate the population, and collisions are unable to remove red giants out to 0.4 pc although brighter red giants much closer in may be depleted via collisions with stellar-mass black holes. For a much flatter IMF, the stellar population is dominated by compact remnants (i.e. black holes, white dwarfs and neutron stars). The most common collisions are then those between main-sequence stars and compact remnants. Such encounters are likely to destroy the main-sequence stars and thus prevent their evolution into red giants. In this way, the red-giant population could be depleted out to 0.4 pc matching observations. If this is the case, it implies the Galactic center contains a much larger population of stellar-mass black holes than would be expected from a regular IMF. This may in turn have implications for the formation and growth of the central supermassive black hole.

preprint2009arXiv

The origin of very wide binary systems

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

preprint2005arXiv

An origin for short g-ray bursts unassociated with current star formation

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

preprint2002arXiv

High Resolution Calculations of Merging Neutron Stars I: Model Description and Hydrodynamic Evolution

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

preprint2001arXiv

Post Newtonian SPH

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

preprint1998arXiv

Coalescing Neutron Stars: A Solution to the R-Process Problem?

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

preprint1998arXiv

Mass ejection in neutron star mergers

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

preprint1996arXiv

Close Approach during Hard Binary--Binary Scattering

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

preprint1995arXiv

A Stellar Audit: The Computation of Encounter Rates for 47 Tucanae and omega Centauri

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

preprint1995arXiv

The Binary Zoo: The Calculation of Production Rates of Binaries Through 2+1 Encounters in Globular Clusters

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