BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Christopher A. Tout

Christopher A. Tout appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRastro-ph.HEastro-ph.EPastro-phastro-ph.GAphysics.flu-dynastro-ph.COphysics.gen-ph

Catalog footprint

What is connected

45works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

45 published item(s)

preprint2024arXiv

One model to rule them all: magnetic braking from CVs to low-mass stars

We present the results of the study of cataclysmic variables (CVs) and AM Canum Venaticorum (AM CVn) stars with our double dynamo (DD) formalism of angular momentum loss (AML) by magnetic braking (MB). We show that (1) our MB model reproduces the period gap ($2\lesssim P_\mathrm{orb}/\,\mathrm{hr}\lesssim3$) and the period minimum spike ($P_\mathrm{orb}\approx 80\, \mathrm{min}$) in CV distribution, (2) evolved CVs, where the donor star commences Roche lobe overflow (RLOF) close to or just beyond the end of the main-sequence, populate the region in and beyond the period gap, and are more likely to be detected at $P_\mathrm{orb}\geq 5.5 \,\mathrm{hr}$. This contaminates the mass-radius fit of long-period CV donors. We show that (3) several evolved CVs become AM CVn stars with $10\lesssim P_\mathrm{orb}/\,\mathrm{min}\lesssim 65$. Their evolution, driven by $\mathrm{AML_{MB}}$ and AML by gravitational radiation (GR, $\mathrm{AML_{GR}}$), leaves them extremely H-exhausted to the point of being indistinguishable from AM CVn stars formed via the He-star and the White Dwarf (WD) channels in terms of the absence of H in their spectra. We further show that (4) owing to the presence of a significant radiative region, intermediate-mass giants/sub-giants, which are progenitors of AM CVn stars formed through the He-star channel, may undergo common envelope evolution that does not behave classically, (5) several AM CVn systems with extremely bloated donors, such as Gaia14aae, ZTFJ1637+49 and SRGeJ045359.9+622444 do not match any modelled trajectories if these systems are modelled only with $\mathrm{AML_{GR}}$, (6) the uncertainties in MB greatly affect modelling results. This, in turn, affects our efforts to distinguish between different AM CVn formation channels and their relative importance. Finally, we find that (7) a similar MB prescription also explains the spin-down of single, low-mass stars.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A unified model for the evolution of cataclysmic variables

We give an updated version of the analytical equation of state used in the Cambridge stellar evolution code (STARS) as a free to use open-source package that we have used to model cool white dwarfs down to temperatures $\log_{10}(T_\mathrm{eff}/\mathrm{K})\:=\;3$. With this update in the STARS code we model the secular evolution of cataclysmic variable (CV) stars using a double dynamo model wherein there is an interplay between two $α-Ω$ dynamos, one in the convective envelope and the other at the boundary of a slowly rotating shrinking radiative core and the growing convective envelope. We confirm that this model provides a physical formalism for the interrupted magnetic braking paradigm. In addition, our model also provides a mechanism for extra angular momentum loss below the period gap. We construct the relative probability distribution of orbital periods $P_\mathrm{orb}$ using the {mass} distribution of white dwarfs in cataclysmic variables and find that our model excellently reproduces the period gap and the observed period minimum spike in CV distribution. We also compare the evolutionary trajectories from our model with those of other empirical models and find agreement between the two. We also report good agreement between our modelled systems and observational data.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Equilibrium Tide: An Updated Prescription for Population Synthesis Codes

We present an updated prescription for the equilibrium tides suitable for population synthesis codes. A grid of 1D evolutionary models was created and the viscous time-scale was calculated for each detailed model. A metallicity dependent power-law relation was fitted to both the convective cores and convective envelopes of the models. The prescription was implemented into the population synthesis code BSE and predicts an 16.5% reduction in the overall number of merges, with those involving main-sequence stars most affected. The new prescription also reduces the overall supernova rate by 3.6% with individual channels being differently affected. The single degenerate Ia supernova occurrence is reduced by 12.8%. The merging of two Carbon Oxygen white dwarfs to cause a Ia supernova occurs 16% less frequently. The number of sub-synchronously rotating stars in close binaries is substantially increased with our prescription, as is the number of non-circularized systems at the start of common-envelope evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Effect of field dissipation and cooling on the mass-radius relation of strongly magnetised white dwarfs

We investigate the luminosity suppression and its effect on the mass-radius relation as well as cooling evolution of highly magnetised white dwarfs. Based on the effect of magnetic field relative to gravitational energy, we suitably modify our treatment of the radiative opacity, magnetostatic equilibrium and degenerate core equation of state to obtain the structural properties of these stars. Although the Chandrasekhar mass limit is retained in the absence of magnetic field and irrespective of the luminosity, strong central fields of about $10^{14}\, {\rm G}$ can yield super-Chandrasekhar white dwarfs with masses up to $1.9\, M_{\odot}$. Smaller white dwarfs tend to remain super-Chandrasekhar for sufficiently strong central fields even when their luminosity is significantly suppressed to $10^{-16}\ L_{\odot}$. Owing to the cooling evolution and simultaneous field decay over $10\ {\rm Gyr}$, the limiting masses of small magnetised white dwarfs can fall to $1.5\ M_{\odot}$ over time. However the majority of these systems still remain practically hidden throughout their cooling evolution because of their high fields and correspondingly low luminosities. Utilising the stellar evolution code $\textit{STARS}$, we obtain close agreement with the analytical mass limit estimates and this suggests that our analytical formalism is physically motivated. Our results argue that super-Chandrasekhar white dwarfs born due to strong field effects may not remain so for long. This explains their apparent scarcity in addition to making them hard to detect because of their suppressed luminosities.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Convective Differential Rotation in Stars and Planets I: Theory

We derive the scaling of differential rotation in both slowly- and rapidly-rotating convection zones using order of magnitude methods. Our calculations apply across stars and fluid planets and all rotation rates, as well as to both magnetized and purely hydrodynamic systems. We find shear $|R\nablaΩ|$ of order the angular frequency $Ω$ for slowly-rotating systems with $Ω\ll |N|$, where $N$ is the \brvs\ frequency, and find that it declines as a power-law in $Ω$ for rapidly-rotating systems with $Ω\gg |N|$. We further calculate the meridional circulation rate and baroclinicity and examine the magnetic field strength in the rapidly rotating limit. Our results are in general agreement with simulations and observations and we perform a detailed comparison with those in a companion paper.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Convective Differential Rotation in Stars and Planets II: Observational and Numerical Tests

Differential rotation is central to a great many mysteries in stars and planets. In Part I we predicted the order of magnitude and scaling of the differential rotation in both hydrodynamic and magnetohydrodynamic convection zones. Our results apply to both slowly- and rapidly-rotating systems, and provide a general picture of differential rotation in stars and fluid planets. We further calculated the scalings of the meridional circulation, entropy gradient and baroclinicity. In this companion paper we compare these predictions with a variety of observations and numerical simulations. With a few exceptions we find that these are consistent in both the slowly-rotating and rapidly-rotating limits. Our results help to localize core-envelope shear in red~giant stars, suggest a rotation-dependent frequency shift in the internal gravity waves of massive stars and potentially explain observed deviations from von Zeipel's gravity darkening in late-type stars.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Suppression of luminosity and mass-radius relation of highly magnetized white dwarfs

We explore the luminosity L of magnetized white dwarfs and its effect on the mass-radius relation. We self-consistently obtain the interface between the electron degenerate gas dominated inner core and the outer ideal gas surface layer or envelope by incorporating both the components of gas throughout the model white dwarf. This is obtained by solving the set of magnetostatic equilibrium, photon diffusion and mass conservation equations in the Newtonian framework, for different sets of luminosity and magnetic field. We appropriately use magnetic opacity, instead of Kramer's opacity, wherever required. We show that the Chandrasekhar-limit is retained, even at high luminosity upto about 10^{-2} solar luminosity but without magnetic field, if the temperature is set constant inside the interface. However there is an increased mass for large-radius white dwarfs, an effect of photon diffusion. Nevertheless, in the presence of strong magnetic fields, with central strength of about 10^{14} G, super-Chandrasekhar white dwarfs, with masses of about 1.9 solar mass, are obtained even when the temperature inside the interface is kept constant. Most interestingly, small-radius magnetic white dwarfs remain super-Chandrasekhar even if their luminosity decreases to as low as about 10^{-20} solar luminosity. However, their large-radius counterparts in the same mass-radius relation merge with Chandrasekhar's result at low L. Hence, we argue for the possibility of highly magnetized, low luminous super-Chandrasekhar mass white dwarfs which, owing to their faintness, can be practically hidden.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Unresolved stellar companions with Gaia DR2 astrometry

For stars with unresolved companions, motions of the centre of light and that of mass decouple, causing a single-source astrometric model to perform poorly. We show that such stars can be easily detected with the reduced chi2 statistic, or RUWE, provided as part of Gaia DR2. We convert RUWE into the amplitude of the image centroid wobble, which, if scaled by the source distance, is proportional to the physical separation between companions (for periods up to several years). We test this idea on a sample of known spectroscopic binaries and demonstrate that the amplitude of the centroid perturbation scales with the binary period and the mass ratio as expected. We apply this technique to the Gaia DR2 data and show how the binary fraction evolves across the Hertzsprung--Russell diagram. The observed incidence of unresolved companions is high for massive young stars and drops steadily with stellar mass, reaching its lowest levels for white dwarfs. We highlight the elevated binary fraction for the nearby Blue Stragglers and Blue Horizontal Branch stars. We also illustrate how unresolved hierarchical triples inflate the relative velocity signal in wide binaries. Finally, we point out a hint of evidence for the existence of additional companions to the hosts of extrasolar hot jupiters.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Cosmic Biology in Perspective

A series of astronomical observations obtained over the period 1986 to 2018 supports the idea that life is a cosmic rather than a purely terrestrial or planetary phenomenon. These include (1) the detection of biologically relevant molecules in interstellar clouds and in comets, (2) mid-infrared spectra of interstellar grains and the dust from comets, (3) a diverse set of data from comets including the Rosetta mission showing consistency with biology and (4) the frequency of Earth-like or habitable planets in the Galaxy. We argue that the conjunction of all the available data suggests the operation of cometary biology and interstellar panspermia rather than the much weaker hypothesis of comets being only the source of the chemical building blocks of life. We conclude with specific predictions on the properties expected of extra-terrestrial life if it is discovered on Enceladus, Europa or beyond. A radically different biochemistry elsewhere can be considered as a falsification of the theory of interstellar panspermia.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A dynamical gravitational wave source in a dense cluster

Making use of a new N-body model to describe the evolution of a moderate-size globular cluster we investigate the characteristics of the population of black holes within such a cluster. This model reaches core-collapse and achieves a peak central density typical of the dense globular clusters of the Milky Way. Within this high-density environment we see direct confirmation of the merging of two stellar remnant black-holes in a dynamically-formed binary, a gravitational wave source. We describe how the formation, evolution and ultimate ejection/destruction of binary systems containing black holes impacts the evolution of the cluster core. Also, through comparison with previous models of lower density, we show that the period distribution of black hole binaries formed through dynamical interactions in this high-density model favours the production of gravitational wave sources. We confirm that the number of black holes remaining in a star cluster at late times and the characteristics of the binary black hole population depend on the nature of the star cluster, critically on the number density of stars and by extension the relaxation timescale.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The proper motion of HV2112: A TŻO candidate in the SMC

The candidate Thorne-Żytkow object (TŻO), HV2112, is becoming a well-studied if enigmatic object. A key point of its candidacy as a TŻO is whether or not it resides in the Small Magellanic Cloud (SMC). HV2112 has detections in a series of photometric catalogues which have resulted in contradictory estimates of its proper motion and, therefore, its membership within the SMC. This letter seeks to resolve the issue of the SMC membership of HV2112 through a reanalysis of extant photometric data. We also demonstrate the difficulties and downfalls inherent in considering a range of catalogue proper motions. We conclude that the proper motion, and associated ancillary radial velocity, positional and photometric properties, are fully consistent with HV2112 being within the SMC and thus it remains a candidate TŻO.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Core radii and common-envelope evolution

Many classes of objects and events are thought to form in binary star systems after a phase in which a core and companion spiral to smaller separation inside a common envelope (CE).Such a phase can end with the merging of the two stars or with the ejection of the envelope to leave a surviving binary system.The outcome is usually predicted by calculating the separation to which the stars must spiral to eject the envelope, assuming that the ratio of the core--envelope binding energy to the change in orbital energy is equal to a constant efficiency factor $α$. If either object would overfill its Roche lobe at this end-of-CE separation, then the stars are assumed to merge. It is unclear what critical radius should be compared to the end-of-CE Roche lobe for stars which have developed cores before the start of a CE phase. After improving the core radius formulae in the widely used BSE rapid evolution code, we compare the properties of populations in which the critical radius is chosen to be the pre-CE core radius or the post-CE stripped remnant radius. Our improvements to the core radius formulae and the uncertainty in the critical radius significantly affect the rates of merging in CE phases of most types. We find the types of systems for which these changes are most important.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Formation of redbacks via accretion induced collapse

We examine the growing class of binary millisecond pulsars known as redbacks. In these systems the pulsar's companion has a mass between 0.1 and about 0.5 solar masses in an orbital period of less than 1.5 days. All show extended radio eclipses associated with circumbinary material. They do not lie on the period-companion mass relation expected from the canonical intermediate-mass X-ray binary evolution in which the companion filled its Roche lobe as a red giant and has now lost its envelope and cooled as a white dwarf. The redbacks lie closer to, but usually at higher period than, the period-companion mass relation followed by cataclysmic variables and low-mass X-ray binaries. In order to turn on as a pulsar mass accretion on to a neutron star must be sufficiently weak, considerably weaker than expected in systems with low-mass main-sequence companions driven together by magnetic braking or gravitational radiation. If a neutron star is formed by accretion induced collapse of a white dwarf as it approaches the Chandrasekhar limit some baryonic mass is abruptly lost to its binding energy so that its effective gravitational mass falls. We propose that redbacks form when accretion induced collapse of a white dwarf takes place during cataclysmic variable binary evolution because the loss of gravitational mass makes the orbit expand suddenly so that the companion no longer fills its Roche lobe. Once activated, the pulsar can ablate its companion and so further expand the orbit and also account for the extended eclipses in the radio emission of the pulsar that are characteristic of these systems. The whole period-companion mass space occupied by the redbacks can be populated in this way.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Merging Binary Stars and the magnetic white dwarfs

A magnetic dynamo driven by differential rotation generated when stars merge can explain strong fields in certain classes of magnetic stars, including the high field magnetic white dwarfs (HFMWDs). In their case the site of the differential rotation has been variously proposed to be within a common envelope, the massive hot outer regions of a merged degenerate core or an accretion disc formed by a tidally disrupted companion that is subsequently incorporated into a degenerate core. We synthesize a population of binary systems to investigate the stellar merging hypothesis for observed single HFMWDs. Our calculations provide mass distribution and the fractions of white dwarfs that merge during a common envelope phase or as double degenerate systems in a post common envelope phase. We vary the common envelope efficiency parameter alpha and compare with observations. We find that this hypothesis can explain both the observed incidence of magnetism and the mass distribution of HFMWDs for a wide range of alpha. In this model, the majority of the HFMWDs are of the Carbon Oxygen type and merge within a common envelope. Less than about a quarter of a per cent of HFMWDs originate from double degenerate stars that merge after common envelope evolution and these populate the high-mass tail of the HFMWD mass distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Pre-mainsequence stellar evolution in N-body models

We provide a set of analytic fits to the radii of pre-mainsequence stars in the mass range 0.1 < M/Msun < 8.0. We incorporate the formulae in N-body cluster models for evolution from the beginning of pre-main sequence. In models with 1,000 stars and high initial cluster densities, pre-mainsequence evolution causes roughly twice the number of collisions between stars than in similar models with evolution begun only from the zero-age main sequence. The collisions are often all part of a runaway sequence that creates one relatively massive star.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

An Exoplanet's Response to Anisotropic Stellar Mass-Loss During Birth and Death

The birth and death of planets may be affected by mass outflows from their parent stars during the T-Tauri or post-main-sequence phases of stellar evolution. These outflows are often modelled to be isotropic, but this assumption is not realistic for fast rotators, bipolar jets and supernovae. Here we derive the general equations of motion for the time evolution of a single planet, brown dwarf, comet or asteroid perturbed by anisotropic mass loss in terms of a complete set of planetary orbital elements, the ejecta velocity, and the parent star's co-latitude and longitude. We restrict our application of these equations to 1) rapidly rotating giant stars, and 2) arbitrarily-directed jet outflows. We conclude that the isotropic mass-loss assumption can safely be used to model planetary motion during giant branch phases of stellar evolution within distances of hundreds of au. In fact, latitudinal mass loss variations anisotropically affect planetary motion only if the mass loss is asymmetric about the stellar equator. Also, we demonstrate how constant-velocity, asymmetric bipolar outflows in young systems incite orbital inclination changes. Consequently, this phenomenon readily tilts exoplanetary orbits external to a nascent disc on the order of degrees.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Planetary nebulae after common-envelope phases initiated by low-mass red giants

It is likely that at least some planetary nebulae are composed of matter which was ejected from a binary star system during common-envelope (CE) evolution. For these planetary nebulae the ionizing component is the hot and luminous remnant of a giant which had its envelope ejected by a companion in the process of spiralling-in to its current short-period orbit. A large fraction of CE phases which end with ejection of the envelope are thought to be initiated by low-mass red giants, giants with inert, degenerate helium cores. We discuss the possible end-of-CE structures of such stars and their subsequent evolution to investigate for which structures planetary nebulae are formed. We assume that a planetary nebula forms if the remnant reaches an effective temperature greater than 30 kK within 10^4 yr of ejecting its envelope. We assume that the composition profile is unchanged during the CE phase so that possible remnant structures are parametrized by the end-of-CE core mass, envelope mass and entropy profile. We find that planetary nebulae are expected in post-CE systems with core masses greater than about 0.3 solar masses if remnants end the CE phase in thermal equilibrium. We show that whether the remnant undergoes a pre-white dwarf plateau phase depends on the prescribed end-of-CE envelope mass. Thus, observing a young post-CE system would constrain the end-of CE envelope mass and post-CE evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Great Escape III: Placing post-main-sequence evolution of planetary and binary systems in a Galactic context

Our improving understanding of the life cycle of planetary systems prompts investigations of the role of the Galactic environment before, during and after Asymptotic Giant Branch (AGB) stellar evolution. Here, we investigate the interplay between stellar mass loss, Galactic tidal perturbations, and stellar flybys for evolving stars which host one planet, smaller body or stellar binary companion and reside in the Milky Way's bulge or disc. We find that the potential evolutionary pathways from a main sequence (MS) to a white dwarf (WD) planetary system are a strong function of Galactocentric distance only with respect to the prevalence of stellar flybys. Planetary ejection and collision with the parent star should be more common towards the bulge. At a given location anywhere in the Galaxy, if the mass loss is adiabatic, then the secondary is likely to avoid close flybys during AGB evolution, and cannot eventually escape the resulting WD because of Galactic tides alone. Partly because AGB mass loss will shrink a planetary system's Hill ellipsoid axes by about 20 to 40 per cent, Oort clouds orbiting WDs are likely to be more depleted and dynamically excited than on the MS.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Most Magnetic Stars

Observations of magnetic A, B and O stars show that the poloidal magnetic flux per unit mass has an upper bound of 10^-6.5 G cm^2/g. A similar upper bound is found for magnetic white dwarfs even though the highest magnetic field strengths at their surfaces are much larger. For magnetic A and B stars there also appears to be a well defined lower bound below which the incidence of magnetism declines rapidly. According to recent hypotheses, both groups of stars may result from merging stars and owe their strong magnetism to fields generated by a dynamo mechanism as they merge. We postulate a simple dynamo that generates magnetic field from differential rotation. The growth of magnetic fields is limited by the requirement that the poloidal field stabilizes the toroidal and vice versa. While magnetic torques dissipate the differential rotation, toroidal field is generated from poloidal by an Omega dynamo. We further suppose that mechanisms that lead to the decay of toroidal field lead to the generation of poloidal. Both poloidal and toroidal fields reach a stable configuration which is independent of the size of small initial seed fields but proportional to the initial differential rotation. We pose the hypothesis that strongly magnetic stars form from the merging of two stellar objects. The highest fields are generated when the merge introduces differential rotation that amounts to critical break up velocity within the condensed object. Calibration of a simplistic dynamo model with the observed maximum flux per unit mass for main-sequence stars and white dwarfs indicates that about 1.5x10^-4 of the decaying toroidal flux must appear as poloidal. The highest fields in single white dwarfs are generated when two degenerate cores merge inside a common envelope or when two white dwarfs merge by gravitational-radiation angular momentum loss.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Nature of Millisecond Pulsars with Helium White Dwarf Companions

We examine the growing data set of binary millisecond pulsars that are thought to have a helium white dwarf companion. These systems are believed to form when a low- to intermediate-mass companion to a neutron star fills its Roche lobe between central hydrogen exhaustion and core helium ignition. We confirm that our own stellar models reproduce a well-defined period-companion mass relation irrespective of the details of the mass transfer process. With magnetic braking this relation extends to periods of less than 1d for a 1Msun giant donor. With this and the measured binary mass functions we calculate the orbital inclination of each system for a given pulsar mass. We expect these inclinations to be randomly oriented in space. If the masses of the pulsars were typically 1.35Msun then there would appear to be a distinct dearth of high-inclination systems. However if the pulsar masses are more typically 1.55 to 1.65Msun then the distribution of inclinations is indeed indistinguishable from random. If it were as much as 1.75Msun then there would appear to be an excess of high-inclination systems. Thus with the available data we can argue that the neutron star masses in binary millisecond pulsars recycled by mass transfer from a red giant typically lie around 1.6Msun and that there is no preferred inclination at which these systems are observed. Hence there is reason to believe that pulsar beams are either sufficiently broad or show no preferred direction relative to the pulsar's spin axis which is aligned with the binary orbit. This is contrary to some previous claims, based on a subset of the data available today, that there might be a tendency for the pulsar beams to be perpendicular to their spin.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Quasi-stars, giants and the Schönberg-Chandrasekhar limit

The Schönberg-Chandrasekhar (SC) limit is a well-established result in the understanding of stellar evolution. It provides an estimate of the point at which an evolved isothermal core embedded in an extended envelope begins to contract. We investigate contours of constant fractional mass in terms of homology invariant variables U and V and find that the SC limit exists because the isothermal core solution does not intersect all the contours for an envelope with polytropic index 3. We find that this analysis also applies to similar limits in the literature including the inner mass limit for polytropic models of quasi-stars. Consequently, any core solution that does not intersect all the fractional mass contours exhibits an associated limit and we identify several relevant cases where this is so. We show that a composite polytrope is at a fractional core mass limit when its core solution touches but does not cross the contour of the corresponding fractional core mass. We apply this test to realistic models of helium stars and find that stars typically expand when their cores are near a mass limit. Furthermore, it appears that stars that evolve into giants have always first exceeded an SC-like limit.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Stellar evolution of massive stars with a radiative alpha-omega dynamo

Models of rotationally-driven dynamos in stellar radiative zones have suggested that magnetohydrodynamic transport of angular momentum and chemical composition can dominate over the otherwise purely hydrodynamic processes. A proper consideration of the interaction between rotation and magnetic fields is therefore essential. Previous studies have focused on a magnetic model where the magnetic field strength is derived as a function of the stellar structure and angular momentum distribution. We have adapted our one-dimensional stellar rotation code, RoSE, to model the poloidal and toroidal magnetic field strengths with a pair of time-dependent advection-diffusion equations coupled to the equations for the evolution of the angular momentum distribution and stellar structure. This produces a much more complete, though still reasonably simple, model for the magnetic field evolution. Our model reproduces well observed surface nitrogen enrichment of massive stars in the Large Magellanic Cloud. In particular it reproduces a population of slowly-rotating nitrogen-enriched stars that cannot be explained by rotational mixing alone alongside the traditional rotationlly-enriched stars. The model further predicts a strong mass-dependency for the dynamo-driven field. Above a threshold mass, the strength of the magnetic dynamo decreases abruptly and so we predict that more massive stars are much less likely to support a dynamo-driven field than less massive stars.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Great Escape II: Exoplanet Ejection from Dying Multiple Star Systems

Extrasolar planets and belts of debris orbiting post-main-sequence single stars may become unbound as the evolving star loses mass. In multiple star systems, the presence or co-evolution of the additional stars can significantly complicate the prospects for orbital excitation and escape. Here, we investigate the dynamical consequences of multi-phasic, nonlinear mass loss and establish a criterion for a system of any stellar multiplicity to retain a planet whose orbit surrounds all of the parent stars. For single stars which become white dwarfs, this criterion can be combined with the Chandrasekhar Limit to establish the maximum allowable mass loss rate for planet retention. We then apply the criterion to circumbinary planets in evolving binary systems over the entire stellar mass phase space. Through about 10^5 stellar evolutionary track realizations, we characterize planetary ejection prospects as a function of binary separation, stellar mass and metallicity. This investigation reveals that planets residing at just a few tens of AU from a central concentration of stars are susceptible to escape in a wide variety of multiple systems. Further, planets are significantly more susceptible to ejection from multiple star systems than from single star systems for a given system mass. For system masses greater than about 2 Solar masses, multiple star systems represent the greater source of free-floating planets.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The structure and evolution of quasi-stars

The existence of bright quasars at high redshifts implies that supermassive black holes were able to form in the early Universe. Though a number of mechanisms to achieve this have been proposed, none yet stands out. A recent suggestion is the formation of quasi-stars, initially stellar-mass black holes accreting from hydrostatic giant-like envelopes of gas, formed from the monolithic collapse of pre-galactic gas clouds. In this work, we modify the Cambridge STARS stellar evolution package to construct detailed models of the evolution of these objects. We find that, in all of our models, the black hole inside the envelope is able to reach slightly more than one-tenth of the total mass of the system before hydrostatic equilibrium breaks down. This breakdown occurs after a few million years of evolution. We show that the mechanism which causes the hydrostatic evolution to end is present in polytropic models. We also show that the solutions are highly sensitive to the size of the inner boundary radius and that no physical solutions exist if the inner boundary is chosen to be less than about 0.3 of the Bondi radius.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Towards a unified model of stellar rotation II: Model-dependent characteristics of stellar populations

Rotation has a number of important effects on the evolution of stars. Apart from structural changes because of the centrifugal force, turbulent mixing and meridional circulation caused by rotation can dramatically affect a star's chemical evolution. This leads to changes in the surface temperature and luminosity as well as modifying its lifetime. Observationally rotation decreases the surface gravity, causes enhanced mass loss and leads to surface abundance anomalies of various chemical isotopes. The replication of these physical effects with simple stellar evolution models is very difficult and has resulted in the use of numerous different formulations to describe the physics. Using stellar evolution calculations based on several physical models we discuss the features of the resulting simulated stellar populations which can help to distinguish between the models.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Interacting compact binaries: modeling mass transfer in eccentric systems

We discuss mass transfer in eccentric binaries containing a white dwarf and a neutron star (WD--NS binaries). We show that such binaries are produced from field binaries following a series of mass transfer episodes that allow the white dwarf to form before the neutron star. We predict the orbital properties of binaries similar to the observed WD--NS binary J1141+6545, and show that they will undergo episodic mass transfer from the white dwarf to the neutron star. Furthermore, we describe oil-on-water, a two-phase SPH formalism that we have developed in order to model mass transfer in such binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Runaway stars as progenitors of supernovae and gamma-ray bursts

When a core collapse supernova occurs in a binary system, the surviving star as well as the compact remnant emerging from the SN, may reach a substantial space velocity. With binary population synthesis modelling at solar and one fifth of solar metallicity, we predict the velocities of such runaway stars or binaries. We compile predictions for runaway OB stars, red supergiants and Wolf-Rayet stars. For those stars or binaries which undergo a second stellar explosion we compute their further evolution and the distance travelled until a Type II or Type Ibc SN or a long or short gamma-ray burst occurs. We find our predicted population of OB runaway stars broadly matches the observed population of stars but, to match the fastest observed WR runaway stars, we require that black holes receive an asymmetric kick upon formation. We find that at solar metallicity Type Ic SN progenitors travel shorter distances than the progenitors of other SN types because they are typically more massive and thus have shorter lifetimes. Those of Type IIP SN can fly farthest about 48 pc on average at solar metallicity. In considering the consequences of assuming that the progenitors of long GRBs are spun-up secondary stars that experience quasi-homogeneous evolution, we find that such evolution has a dramatic effect on the population of runaway WR stars and that some 30 per cent of GRBs could occur a hundred parsecs or more from their initial positions. We also consider mergers of double compact object binaries consisting of neutron stars and/or black holes. We find the most common type of visible mergers are neutron star--black hole mergers that are roughly ten times more common than neutron star--neutron star mergers. We also find that there may be a population of low-velocity neutron stars that are ejected from a binary rather than by their own natal kick.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Slowing down atomic diffusion in subdwarf B stars: mass loss or turbulence?

Subdwarf B stars show chemical peculiarities that cannot be explained by diffusion theory alone. Both mass loss and turbulence have been invoked to slow down atomic diffusion in order to match observed abundances. The fact that some sdB stars show pulsations gives upper limits on the amount of mass loss and turbulent mixing allowed. Consequently, non-adiabatic asteroseismology has the potential to decide which process is responsible for the abundance anomalies. We compute for the first time seismic properties of sdB models with atomic diffusion included consistently during the stellar evolution. The diffusion equations with radiative forces are solved for H, He, C, N, O, Ne, Mg, Fe and Ni. We examine the effects of various mass-loss rates and mixed surface masses on the abundances and mode stability. It is shown that the mass-loss rates needed to simulate the observed He abundances (10^{-14}<=Mdot [Msun/yr]<=10^{-13}) are not consistent with observed pulsations. We find that for pulsations to be driven the rates should be Mdot<=10^{-15} Msun/yr. On the other hand, weak turbulent mixing of the outer 10^{-6} Msun can explain the He abundance anomalies while still allowing pulsations to be driven. The origin of the turbulence remains unknown but the presence of pulsations gives tight constraints on the underlying turbulence model.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Spindown of massive rotating stars

Models of rapidly rotating massive stars at low metallicities show significantly different evolution and higher metal yields compared to non-rotating stars. We estimate the spin-down time-scale of rapid rotating non-convective stars supporting an alpha-Omega dynamo. The magnetic dynamo gives rise to mass loss in a magnetically controlled stellar wind and hence stellar spin down owing to loss of angular momentum. The dynamo is maintained by strong horizontal rotation-driven turbulence which dominates over the Parker instability. We calculate the spin-down time-scale and find that it could be relatively short, a small fraction of the main-sequence lifetime. The spin-down time-scale decreases dramatically for higher surface rotations suggesting that rapid rotators may only exhibit such high surface velocities for a short time, only a small fraction of their main-sequence lifetime.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Tidal Warping of Be Star Decretion Discs

Rapidly rotating Be stars are observed as shell stars when the decretion disc is viewed edge on. Transitions between the two implies that the discs may be warped and precessing. Type II X-ray outbursts are thought to occur when the warped disc interacts with the fast stellar wind. We suggest that tides from a misaligned companion neutron star can cause the observed effects. We make numerical models of a Be star decretion disc in which the spin of the Be star is misaligned with the orbital axis of a neutron star companion. Tidal torques from the neutron star truncate the disc at a radius small enough that the neutron star orbit does not intersect the disc unless the eccentricity or misalignment is very large. A magnetic torque from the Be star that is largest at the equator, where the rotation is fastest, is approximated by an inner boundary condition. There are large oscillations in the mass and inclination of the disc as it moves towards a steady state. These large variations may explain the observed changes from Be star to Be shell star and vice-versa and also the Type II X-ray outbursts. We find the tidal timescale on which the disc warps, precesses and reaches a steady state to be around a year up to a few hundred years. If present, the oscillations in mass and disc inclination occur on a fraction of this timescale depending on the orbital parameters of the binary. The timescales associated with the tidal torque for observed Be star binaries suggest that these effects are important in all but the longest period binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Tidally enhanced stellar wind: a way to make the symbiotic channel to Type Ia supernova viable

In the symbiotic (or WD+RG) channel of the single degenerate scenario for type Ia supernovae (SNe Ia) the explosions occur a relatively long time after star formation. The birthrate from this channel would be too low to account for all observed SNe Ia were it not for some mechanism to enhance the rate of accretion on to the white dwarf. A tidally enhanced stellar wind, of the type which has been postulated to explain many phenomena related to giant star evolution in binary systems, can do this. Compared to mass-stripping this model extends the space of SNe Ia progenitors to longer orbital periods and hence increases the birthrate to about $0.0069\,{\rm yr^{-1}}$ for the symbiotic channel. Two symbiotic stars, T CrB and RS Oph, considered to be the most likely progenitors of SNe Ia through the symbiotic channel, are well inside the period--companion mass space predicted by our models.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Towards a unified model of stellar rotation

The effects of rapid rotation on stellar evolution can be profound. We are now beginning to gather enough data to allow a realistic comparison between different physical models. Two key tests for any theory of stellar rotation are first whether it can match observations of the enrichment of nitrogen, and potentially other elements, in clusters containing rapid rotators and secondly whether it can reproduce the observed broadening of the main sequence in the Hertzsprung-Russel diagram. Models of stellar rotation have been steadily increasing in number and complexity over the past two decades but the lack of data makes it difficult to determine whether such additions actually give a closer reflection of reality. One of the most poorly explored features of stellar rotation models is the treatment of angular momentum transport within convective zones. If we treat the core as having uniform specific angular momentum the angular momentum distribution in the star, for a given surface rotation, is dramatically different from what it is when we assume the star rotates as a solid body. The uniform specific angular momentum also generates strong shears which can drive additional transport of chemical elements close to the boundary of a convection zone. A comparison of different models and their reproduction of observable properties with otherwise identical input physics is essential to properly distinguish between them. We compare detailed grids of stellar evolution tracks of intermediate and high-mass stars produced using several models for rotation generated with our new stellar rotation code.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A Common Envelope Binary Star Origin of Long Gamma-ray Bursts

The stellar origin of gamma-ray bursts can be explained by the rapid release of energy in a highly collimated, extremely relativistic jet. This in turn appears to require a rapidly spinning highly magnetised stellar core that collapses into a magnetic neutron star or a black hole within a relatively massive envelope. They appear to be associated with type Ib/c supernovae but, with a birthrate of around 10^{-6}-10^{-5} per year per galaxy, they are considerably rarer than such supernovae in general. To satisfy all these requirements we hypothesize a binary star model that ends with the merging of an oxygen neon white dwarf with the carbon-oxygen core of a naked helium star during a common envelope phase of evolution. The rapid spin and high magnetic field are natural consequences of such a merging. The evolution that leads to these progenitors is convoluted and so naturally occurs only very rarely. To test the hypothesis we evolve a population of progenitors and find that the rate is as required. At low metallicity we calculate that a similar fraction of stars evolve to this point and so would expect the gamma-ray burst rate to correlate with the star formation rate in any galaxy. This too is consistent with observations. These progenitors, being of intermediate mass, differ radically from the usually postulated high-mass stars. Thus we can reconcile observations that the bursts occur close to but not within massive star associations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Does GD 356 have a Terrestrial Planetary Companion?

GD 356 is unique among magnetic white dwarfs because it shows Zeeman-split Balmer lines in pure emission. The lines originate from a region of nearly uniform field strength (delta B/B is approximately 0.1) that covers 10 per cent of the stellar surface in which there is a temperature inversion. The energy source that heats the photosphere remains a mystery but it is likely to be associated with the presence of a companion. Based on current models we use archival Spitzer IRAC observations to place a new and stringent upper limit of 12 Jupiter masses for the mass of such a companion. In the light of this result and the recent discovery of a 115 min photometric period for GD 356, we exclude previous models that invoke accretion and revisit the unipolar inductor model that has been proposed for this system. In this model a highly conducting planet with a metallic core orbits the magnetic white dwarf and, as it cuts through field lines, a current is set flowing between the two bodies. This current dissipates in the photosphere of the white dwarf and causes a temperature inversion. Such a planet is unlikely to have survived the RGB/AGB phases of evolution so we argue that it may have formed from the circumstellar disc of a disrupted He or CO core during a rare merger of two white dwarfs. GD 356 would then be a white dwarf counterpart of the millisecond binary pulsar PSR 1257+12 which is known to host a planetary system.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Spin Angular Momentum Evolution of the Long Period Algols

We consider the spin angular momentum evolution of the accreting components of Algol-type binary stars. In wider Algols the accretion is through a disc so that the accreted material can transfer enough angular momentum to the gainer that material at its equator should be spinning at break-up. We demonstrate that even a small amount of mass transfer, much less than required to produce today's mass ratios, transfers enough angular momentum to spin the gainer up to this critical rotation velocity. However the accretors in these systems have spins typically between 10 and $40\,$per cent of the critical rate. So some mechanism for angular momentum loss from the gainers is required. We consider generation of magnetic fields in the radiative atmospheres in a differentially rotating star and the possibility of angular momentum loss driven by strong stellar winds in the intermediate mass stars, such as the primaries of the Algols. Differential rotation, induced by the accretion itself, may produce such winds which carry away enough angular momentum to reduce their rotational velocities to the today's observed values. We apply this model to two systems with initial periods of 5\,d, one with initial masses 5 and $3\,\rm{M}_{\odot}$ and the other with 3.2 and $2\,\rm{M}_{\odot}$. Our calculations show that, if the mass outflow rate in the stellar wind is about $10\,$per cent of the accretion rate and the dipole magnetic field is stronger than about $1\,$kG, the spin rate of the gainer is reduced to below break-up velocity even in the fast phase of mass transfer. Larger mass loss is needed for smaller magnetic fields. The slow rotation of the gainers in the classical Algol systems is explained by a balance between the spin-up by mass accretion and spin-down by a stellar wind linked to a magnetic field.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Supernova Kicks and Misaligned Be Star Binaries

Be stars are rapidly spinning B stars surrounded by an outflowing disc of gas in Keplerian rotation. Be star/X-ray binary systems contain a Be star and a neutron star. They are found to have non-zero eccentricities and there is evidence that some systems have a misalignment between the spin axis of the star and the spin axis of the binary orbit. The eccentricities in these systems are thought to be caused by a kick to the neutron star during the supernova that formed it. Such kicks would also give rise to misalignments. In this paper we investigate the extent to which the same kick distribution can give rise to both the observed eccentricity distribution and the observed misalignments. We find that a Maxwellian distribution of velocity kicks with a low velocity dispersion, $σ_k \approx 15\rm km s^{-1}$, is consistent with the observed eccentricity distribution but is hard to reconcile with the observed misalignments, typically $i \ge 25^\circ$. Alternatively a higher velocity kick distribution, $σ_k = 265 \rm km s^{-1}$, is consistent with the observed misalignments but not with the observed eccentricities, unless post-supernova circularisation of the binary orbits has taken place. We discuss briefly how this might be achieved.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Supernova Kicks and Misaligned Microquasars

The low-mass X-ray binary microquasar GRO J1655-40 is observed to have a misalignment between the jets and the binary orbital plane. Since the current black hole spin axis is likely to be parallel to the jets, this implies a misalignment between the spin axis of the black hole and the binary orbital plane. It is likely the black holes formed with an asymmetric supernova which caused the orbital axis to misalign with the spin of the stars. We ask whether the null hypothesis that the supernova explosion did not affect the spin axis of the black hole can be ruled out by what can be deduced about the properties of the explosion from the known system parameters. We find that this null hypothesis cannot be disproved but we find that the most likely requirements to form the system include a small natal black hole kick (of a few tens of km/s) and a relatively wide pre-supernova binary. In such cases the observed close binary system could have formed by tidal circularisation without a common envelope phase.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The chemical composition of donors in AM CVn stars and ultra-compact X-ray binaries: observational tests of their formation

We study the formation of ultra-compact binaries (AM CVn stars and ultra-compact X-ray binaries) with emphasis on the surface chemical abundances of the donors in these systems. Hydrogen is not convincingly detected in the spectra of these systems. Three different proposed formation scenarios involve different donor stars, white dwarfs, helium stars or evolved main-sequence stars. Using detailed evolutionary calculations we show that the abundances of helium WD donors and evolved main-sequence stars are close to equilibrium CNO-processed material, and the detailed abundances correlate with the core temperature and thus mass of the MS progenitors. Evolved MS donors typically have traces of H left. For hybrid or CO white dwarf donors, the carbon and oxygen abundances depend on the temperature of the helium burning and thus on the helium core mass of the progenitors. For helium star donors in addition to their mass, the abundances depend strongly on the amount of helium burnt before mass transfer starts and can range from unprocessed and thus almost equal to CNO-processed matter, to strongly processed and thus C/O rich and N-deficient. We briefly discuss the relative frequency of these cases for helium star donors, based on population synthesis results. Finally we give diagnostics for applying our results to observed systems and find that the most important test is the N/C ratio, which can indicate the formation scenario as well as, in some cases, the mass of the progenitor of the donor. In addition, if observed, the N/O, O/He and O/C ratios can distinguish between helium star and WD donors. Applied to the known systems we find evidence for WD donors in the AM CVn systems GP Com, CE 315 and SDSS J0804+16 and evidence for hybrid WD or very evolved helium star donors in the UCXBs 4U 1626-67 and 4U 0614+09. [Abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The evolution of low-metallicity asymptotic giant branch stars and the formation of carbon-enhanced metal-poor stars

We investigate the behaviour of asymptotic giant branch (AGB) stars between metallicities Z = 10-4 and Z = 10-8 . We determine which stars undergo an episode of flash-driven mixing, where protons are ingested into the intershell convection zone, as they enter the thermally pulsing AGB phase and which undergo third dredge-up. We find that flash-driven mixing does not occur above a metallicity of Z = 10-5 for any mass of star and that stars above 2 M do not experience this phenomenon at any metallicity. We find carbon ingestion (CI), the mixing of carbon into the tail of hydrogen burning region, occurs in the mass range 2 M to around 4 M . We suggest that CI may be a weak version of the flash-driven mechanism. We also investigate the effects of convective overshooting on the behaviour of these objects. Our models struggle to explain the frequency of CEMP stars that have both significant carbon and nitrogen enhancement. Carbon can be enhanced through flash-driven mixing, CI or just third dredge up. Nitrogen can be enhanced through hot bottom burning and the occurrence of hot dredge-up also converts carbon into nitrogen. The C/N ratio may be a good indicator of the mass of the primary AGB stars.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Shape of an Accretion Disc in a Misaligned Black Hole Binary

We model the overall shape of an accretion disc in a semi-detached binary system in which mass is transfered on to a spinning black hole the spin axis of which is misaligned with the orbital rotation axis. We assume the disc is in a steady state. Its outer regions are subject to differential precession caused by tidal torques of the companion star. These tend to align the outer parts of the disc with the orbital plane. Its inner regions are subject to differential precession caused by the Lense-Thirring effect. These tend to align the inner parts of the disc with the spin of the black hole. We give full numerical solutions for the shape of the disc for some particular disc parameters. We then show how an analytic approximation to these solutions can be obtained for the case when the disc surface density varies as a power law with radius. These analytic solutions for the shape of the disc are reasonably accurate even for large misalignments and can be simply applied for general disc parameters. They are particularly useful when the numerical solutions would be slow.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Critical Mass Transfer in Double-Degenerate Type Ia Supernovae

Doubly-degenerate binary systems consisting of two white dwarfs both composed of carbon and oxygen and close enough that mass is transferred from the less massive to the more massive are possible progenitors of type Ia supernovae. If the mass transfer rate is slow enough that the accreting white dwarf can reach 1.38 solar masses then it can ignite carbon degenerately at its centre. This can lead to a thermonuclear runaway and thence a supernova explosion. However if the accretion rate is too high the outer layers of the white dwarf heat up too much and carbon ignites there non-degenerately. A series of mild carbon flashes can then propagate inwards and convert the carbon to neon relatively gently. There is no thermonuclear runaway and no supernova. We have determined the critical rate for accretion and find it to be generally about two-fifths of the Eddington rate. In the case of actual bianry stars the mass transfer rate falls off as mass is transferred so we find that this rate can be very much exceeded at the onset of mass transfer. We determine combinations of white dwarf masses that can lead to central degenerate carbon ignition.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Direct N-body Modelling of Stellar Populations: Blue Stragglers in M67

We present a state-of-the-art N-body code which includes a detailed treatment of stellar and binary evolution as well as the cluster dynamics. This code is ideal for investigating all aspects relating to the evolution of star clusters and their stellar populations. It is applicable to open and globular clusters of any age. We use the N-body code to model the blue straggler population of the old open cluster M67. Preliminary calculations with our binary population synthesis code show that binary evolution alone cannot explain the observed numbers or properties of the blue stragglers. On the other hand, our N-body model of M67 generates the required number of blue stragglers and provides formation paths for all the various types found in M67. This demonstrates the effectiveness of the cluster environment in modifying the nature of the stars it contains and highlights the importance of combining dynamics with stellar evolution. We also perform a series of N = 10000 simulations in order to quantify the rate of escape of stars from a cluster subject to the Galactic tidal field.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

The Binary Second Sequence in Cluster Colour--Magnitude Diagrams

We show how the second sequence seen lying above the main sequence in cluster colour magnitude diagrams results from binaries with a large range of mass ratios and not just from those with equal masses. We conclude that the presence of a densely populated second sequence, with only sparse filling in between it and the single star main sequence, does not necessarily imply that binary mass ratios are close to unity.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

The Theoretical Mass--Magnitude Relation of Low-Mass Stars and its Metallicity Dependence

We investigate the dependence of theoretically generated mass - (absolute magnitude) relations on stellar models. Using up to date physics we compute models in the mass range 0.1 < m < 1M_sun. We compare the solar-metallicity models with our older models, with recent models computed by others, and also with an empirical mass - (absolute magnitude) relation that best fits the observed data. At a given mass below 0.6M_sun the effective temperatures differ substantially from model to model. However taken individually each set of models is in good agreement with observations in the mass - luminosity plane. A minimum in the derivative dm/dM_V at M_V = 11.5, which is due to H_2 formation and establishment of a fully convective stellar interior, is present in all photometric bands, for all models. This minimum leads to a maximum in the stellar luminosity function for Galactic disk stars at M_V = 11.5, M_bol = 9.8. Stellar models should locate this maximum in the stellar luminosity function at the same magnitude as observations. Models which incorporate the most realistic theoretical atmospheres and the most recent equation of state and opacities can satisfy this constraint. These models are also in best agreement with the most recent luminosity - (effective temperature) and mass-luminosity data. Each set of our models of a given metallicity (with 0.2 > [Fe/H] > -2.3) shows a maximum in -dm/dM_bol, which moves to brighter bolometric magnitudes with decreasing metallicity. The change in location of the maximum, as a function of [Fe/H], follows the location of structure in luminosity functions for stellar populations with different metal abundances. This structure seen in all observed stellar populations can be accounted for by the mass--luminosity relation.

0 citations0 reviewsNo code linkedOpen access