Researcher profile

P. J. Davis

P. J. Davis contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
8works
0followers
1topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

8 published item(s)

preprint2014arXiv

Binary evolution using the theory of osculating orbits: conservative Algol evolution

Our aim is to calculate the evolution of Algol binaries within the framework of the osculating orbital theory, which considers the perturbing forces acting on the orbit of each star arising from mass exchange via Roche lobe overflow (RLOF). The scheme is compared to results calculated from a `classical' prescription. Using our stellar binary evolution code BINSTAR, we calculate the orbital evolution of Algol binaries undergoing case A and case B mass transfer, by applying the osculating scheme. The velocities of the ejected and accreted material are evaluated by solving the restricted three-body equations of motion, within the ballistic approximation. This allows us to determine the change of linear momentum of each star, and the gravitational force applied by the mass transfer stream. Torques applied on the stellar spins by tides and mass transfer are also considered. Using the osculating formalism gives shorter post-mass transfer orbital periods typically by a factor of 4 compared to the classical scheme, owing to the gravitational force applied onto the stars by the mass transfer stream. Additionally, during the rapid phase of mass exchange, the donor star is spun down on a timescale shorter than the tidal synchronization timescale, leading to sub-synchronous rotation. Consequently, between 15 and 20 per cent of the material leaving the inner-Lagrangian point is accreted back onto the donor (so-called `self-accretion'), further enhancing orbital shrinkage. Self-accretion, and the sink of orbital angular momentum which mass transfer provides, may potentially lead to more contact binaries. Even though Algols are mainly considered, the osculating prescription is applicable to all types of interacting binaries, including those with eccentric orbits.

preprint2014arXiv

The formation of long-period eccentric binaries with a helium white dwarf

The recent discovery of long-period eccentric binaries hosting a He-WD or a sdB star has been challenging binary-star modelling. Based on accurate determinations of the stellar and orbital parameters for IP Eri, a K0 + He-WD system, we propose an evolutionary path that is able to explain the observational properties of this system and, in particular, to account for its high eccentricity (0.25). Our scenario invokes an enhanced-wind mass loss on the first red giant branch (RGB) in order to avoid mass transfer by Roche-lobe overflow, where tides systematically circularize the orbit. We explore how the evolution of the orbital parameters depends on the initial conditions and show that eccentricity can be preserved and even increased if the initial separation is large enough. The low spin velocity of the K0 giant implies that accretion of angular momentum from a (tidally-enhanced) RGB wind should not be efficient.

preprint2013arXiv

Critically-rotating accretors and non-conservative evolution in Algols

During the mass-transfer phase in Algol systems, a large amount of mass and angular momentum are accreted by the gainer star which can be accelerated up to its critical Keplerian velocity. The fate of the gainer once it reaches this critical value is unclear. We investigate the orbital and stellar spin evolution in semi-detached binary systems, specifically for systems with rapidly rotating accretors. Our aim is to better distinguish between the different spin-down mechanisms proposed which can consistently explain the slow rotation observed in Algols' final states and assess the degree of non-conservatism due to the formation of a hotspot. We use our state-of-the-art binary evolution code, Binstar, which incorporates a detailed treatment of the orbital and stellar spin, including all torques due to mass transfer, the interactions between a star and its accretion disc, tidal effects and magnetic braking. We also present a new prescription for mass loss due to the formation of a hotspot based on energy conservation. The coupling between the star and the disc via the boundary layer prevents the gainer from exceeding the critical rotation. Magnetic-field effects, although operating, are not the dominant spin-down mechanism for sensible field strengths. Spin down owing to tides is 2-4 orders of magnitudes too weak to compensate the spinning-up torque due to mass accretion. Moreover, we find that the final separation strongly depends on the spin-down mechanism. The formation of a hotspot leads to a large event of mass loss during the rapid phase of mass transfer. The degree of conservatism strongly depends on the opacity of the impacted material. A statistical study as well as new observational constraints are needed to find the optimal set of parameters (magnetic-field strength, hotspot geometry,...) to reproduce Algol evolutions.

preprint2013arXiv

Mass transfer in eccentric binary systems using the binary evolution code BINSTAR

We present the first calculations of mass transfer via RLOF for a binary system with a significant eccentricity using our new binary stellar evolution code. The study focuses on a 1.50+1.40 Msun main sequence binary with an eccentricity of 0.25, and an orbital period of about 0.7 d. The reaction of the stellar components due to mass transfer is analyzed, and the evolution of mass transfer during the periastron passage is compared to recent smooth particle hydrodynamics (SPH) simulations. The impact of asynchronism and non-zero eccentricity on the Roche lobe radius, and the effects of tidal and rotational deformation on the stars' structures, are also investigated. Calculations were performed using the state-of-the-art binary evolution code BINSTAR, which calculates simultaneously the structure of the two stars and the evolution of the orbital parameters. The evolution of the mass transfer rate during an orbit has a Gaussian-like shape, with a maximum at periastron, in qualitative agreement with SPH simulations. The Roche lobe radius is modified by the donor star's spin and the orbital eccentricity. This has a significant impact on both the duration and the rate of mass transfer. We find that below some critical rotation rate, mass transfer never occurs, while above some threshold, mass is transferred over the entire orbit. Tidal and rotational deformation of the donor star causes it to become over-sized, enhancing the mass transfer rate further by about a factor of ten, leading to non-conservative mass transfer. The modulation of mass transfer rate with orbital phase produces short-term variability in the surface luminosity and radius of each star. The longer-term behaviour shows, in accordance with studies of circular systems with radiative stars, that the donor becomes ever small and under-luminous, while the converse is the case for the accretor.

preprint2011arXiv

Is the common envelope ejection efficiency a function of the binary parameters?

We reconstruct the common envelope (CE) phase for the current sample of observed white dwarf-main sequence post-common envelope binaries (PCEBs). We apply multi-regression analysis in order to investigate whether correlations exist between the CE ejection efficiencies, alpha_CE, inferred from the sample, and the binary parameters: white dwarf mass, secondary mass, orbital period at the point the CE commences, or the orbital period immediately after the CE phase. We do this with and without consideration for the internal energy of the progenitor primary giants&#39; envelope. Our fits should pave the first steps towards an observationally motivated recipe for calculating alpha_CE using the binary parameters at the start of the CE phase, which will be useful for population synthesis calculations or models of compact binary evolution. If we do consider the internal energy of the giants&#39; envelope, we find a statistically significant correlation between alpha_CE and the white dwarf mass. If we do not, a correlation is found between alpha_CE and the orbital period at the point the CE phase commences. Furthermore, if the internal energy of the progenitor primary envelope is taken into account, then the CE ejection efficiencies are within the canonical range 0<alpha_CE<=1, although PCEBs with brown dwarf secondaries still require alpha_CE>=1.

preprint2011arXiv

Post common envelope binaries from SDSS. XII: The orbital period distribution

The complexity of the common envelope phase and of magnetic stellar wind braking currently limits our understanding of close binary evolution. Because of their intrinsically simple structure, observational population studies of white dwarf plus main sequence (WDMS) binaries hold the potential to test theoretical models and constrain their parameters. The Sloan Digital Sky Survey (SDSS) has provided a large and homogeneously selected sample of WDMS binaries, which we are characterising in terms of orbital and stellar parameters. We have obtained radial velocity information for 385 WDMS binaries from follow-up spectroscopy, and for an additional 861 systems from the SDSS sub-spectra. Radial velocity variations identify 191 of these WDMS binaries as post common envelope binaries (PCEBs). Orbital periods of 58 PCEBs were subsequently measured, predominantly from time-resolved spectroscopy, bringing the total number of SDSS PCEBs with orbital parameters to 79. Observational biases inherent to this PCEB sample were evaluated through extensive Monte Carlo simulations. We find that 21-24% of all SDSS WDMS binaries have undergone common envelope evolution, which is in good agreement with published binary population models and high-resolution HST imaging of WDMS binaries unresolved from the ground. The bias corrected orbital period distribution of PCEBs ranges from 1.9 h to 4.3 d and follows approximately a normal distribution in log(Porb), peaking at ~10.3 h. There is no observational evidence for a significant population of PCEBs with periods in the range of days to weeks. The large and homogeneous sample of SDSS WDMS binaries provides the means to test fundamental predictions of binary population models, and hence to observationally constrain the evolution of all close compact binaries.

preprint2009arXiv

A comprehensive population synthesis study of post-common envelope binaries

We apply population synthesis techniques to calculate the present day population of post-common envelope binaries (PCEBs) for a range of theoretical models describing the common envelope (CE) phase. Adopting the canonical energy budget approach we consider models where the ejection efficiency, $α_{\rmn{CE}}$ is either a constant, or a function of the secondary mass. We obtain the envelope binding energy from detailed stellar models of the progenitor primary, with and without the thermal and ionization energy, but we also test a commonly used analytical scaling. We also employ the alternative angular momentum budget approach, known as the $γ$-algorithm. We find that a constant, global value of $α_{\rmn{CE}} \ga 0.1$ can adequately account for the observed population of PCEBs with late spectral-type secondaries. However, this prescription fails to reproduce IK Pegasi, which has a secondary with spectral type A8. We can account for IK Pegasi if we include thermal and ionization energy of the giant&#39;s envelope, or if we use the $γ$-algorithm. However, the $γ$-algorithm predicts local space densities that are 1 to 2 orders of magnitude greater than estimates from observations. In contrast, the canonical energy budget prescription with an initial mass ratio distribution that favours unequal initial mass ratios gives a local space density which is in good agreement with observations, and best reproduces the observed distribution of PCEBs. Finally, all models fail to reproduce the sharp decline for orbital periods, $P_{\rmn{orb}} \ga 1$ d in the orbital period distribution of observed PCEBs, even if we take into account selection effects against systems with long orbital periods and early spectral-type secondaries.