Source author record

O. R. Pols

O. R. Pols 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

17works
4topics
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

17 published item(s)

preprint2025arXiv

Understanding post-red giant branch binaries through stable mass transfer

Post-RGB and post-AGB binaries consist of a primary star that has recently evolved off either the RGB or AGB after losing most of its envelope, and a main-sequence companion. They are distinguished by luminosities below and above the RGB tip, respectively. These systems host a stable, dusty circumbinary disc, characterised by a near-infrared excess. Observed Galactic post-AGB and post-RGB binaries have orbital periods and eccentricities inconsistent with binary population synthesis models. Here, we focus on post-RGB binaries, testing whether stable mass transfer can explain their orbital periods by comparing models with the known sample of 38 Galactic post-RGB binaries. We systematically determined luminosities of Galactic post-RGB and post-AGB binaries through SED fitting. We computed evolution models for low- and intermediate-mass binaries with RGB donors at two metallicities using MESA. We selected stable mass transfer models producing primaries with effective temperatures within the observed range. From these models, we find that low-mass post-RGB binaries should follow strict luminosity-orbital period relations. The Galactic post-RGB binaries seem consistent with these relations if their orbits remained eccentric during mass transfer and if the donor filled its Roche lobe at periastron. However, our models are unable to explain the eccentricities themselves. Moreover, post-mass-transfer ages from our models are much longer than predicted dissipation timescales of circumbinary discs. Stable mass transfer seems to explain the orbital periods of Galactic post-RGB binaries. This formation channel can be tested further by obtaining orbits of additional Galactic systems and Magellanic Cloud candidates via long-term radial velocity monitoring. Gaia DR 4 will improve luminosities of Galactic post-RGB binaries, enabling more accurate comparison with luminosity-orbital period relations.

preprint2016arXiv

Face-on accretion onto a protoplanetary disc

Globular clusters (GCs) are known to harbor multiple stellar populations. To explain these observations Bastian et al. suggested a scenario in which a second population is formed by the accretion of enriched material onto the low-mass stars in the initial GC population. The idea is that the low-mass, pre-main sequence stars sweep up gas expelled by the massive stars of the same generation into their protoplanetary disc as they move through the GC core. We perform simulations with 2 different smoothed particle hydrodynamics codes to investigate if a low-mass star surrounded by a protoplanetary disc can accrete the amount of enriched material required in this scenario. We focus on the gas loading rate onto the disc and star as well as on the lifetime of the disc. We find that the gas loading rate is a factor of 2 smaller than the geometric rate, because the effective cross section of the disc is smaller than its surface area. The loading rate is consistent for both codes, irrespective of resolution. The disc gains mass in the high resolution runs, but loses angular momentum on a time scale of 10^4 yrs. Two effects determine the loss of (specific) angular momentum in our simulations: 1) continuous ram pressure stripping and 2) accretion of material with no azimuthal angular momentum. Our study and previous work suggest that the former, dominant process is mainly caused by numerical rather than physical effects, while the latter is not. The latter process causes the disc to become more compact, increasing the surface density profile at smaller radii. The disc size is determined in the first place by the ram pressure when the flow first hits the disc. Further evolution is governed by the decrease in the specific angular momentum of the disc. We conclude that the size and lifetime of the disc are probably not sufficient to accrete the amount of mass required in Bastian et al.'s scenario.

preprint2015arXiv

Carbon-enhanced metal-poor stars: a window on AGB nucleosynthesis and binary evolution. I. Detailed analysis of 15 binary stars with known orbital periods

AGB stars are responsible for producing a variety of elements, including carbon, nitrogen, and the heavy elements produced in the slow neutron-capture process ($s$-elements). There are many uncertainties involved in modelling the evolution and nucleosynthesis of AGB stars, and this is especially the case at low metallicity, where most of the stars with high enough masses to enter the AGB have evolved to become white dwarfs and can no longer be observed. The stellar population in the Galactic halo is of low mass ($\lesssim 0.85M_{\odot}$) and only a few observed stars have evolved beyond the first giant branch. However, we have evidence that low-metallicity AGB stars in binary systems have interacted with their low-mass secondary companions in the past. The aim of this work is to investigate AGB nucleosynthesis at low metallicity by studying the surface abundances of chemically peculiar very metal-poor stars of the halo observed in binary systems. To this end we select a sample of 15 carbon- and $s$-element-enhanced metal-poor (CEMP-$s$) halo stars that are found in binary systems with measured orbital periods. With our model of binary evolution and AGB nucleosynthesis, we determine the binary configuration that best reproduces, at the same time, the observed orbital period and surface abundances of each star of the sample. The observed periods provide tight constraints on our model of wind mass transfer in binary stars, while the comparison with the observed abundances tests our model of AGB nucleosynthesis.

preprint2015arXiv

Carbon-enhanced metal-poor stars: a window on AGB nucleosynthesis and binary evolution. II. Statistical analysis of a sample of 67 CEMP-$s$ stars

Many observed CEMP stars are found in binary systems and show enhanced abundances of $s$-elements. The origin of the chemical abundances of these CEMP-$s$ stars is believed to be accretion in the past of enriched material from a primary star in the AGB phase. We investigate the mechanism of mass transfer and the process of nucleosynthesis in low-metallicity AGB stars by modelling the binary systems in which the observed CEMP-$s$ stars were formed. For this purpose we compare a sample of $67$ CEMP-$s$ stars with a grid of binary stars generated by our binary evolution and nucleosynthesis model. We classify our sample CEMP-$s$ stars in three groups based on the observed abundance of europium. In CEMP$-s/r$ stars the europium-to-iron ratio is more than ten times higher than in the Sun, whereas it is lower than this threshold in CEMP$-s/nr$ stars. No measurement of europium is currently available for CEMP-$s/ur$ stars. On average our models reproduce well the abundances observed in CEMP-$s/nr$ stars, whereas in CEMP-$s/r$ stars and CEMP-$s/ur$ stars the abundances of the light-$s$ elements are systematically overpredicted by our models and in CEMP-$s/r$ stars the abundances of the heavy-$s$ elements are underestimated. In all stars our modelled abundances of sodium overestimate the observations. This discrepancy is reduced only in models that underestimate the abundances of most of the $s$-elements. Furthermore, the abundance of lead is underpredicted in most of our model stars. These results point to the limitations of our AGB nucleosynthesis model, particularly in the predictions of the element-to-element ratios. Finally, in our models CEMP-$s$ stars are typically formed in wide systems with periods above 10000 days, while most of the observed CEMP-$s$ stars are found in relatively close orbits with periods below 5000 days.

preprint2015arXiv

Modelling the observed properties of carbon-enhanced metal-poor stars using binary population synthesis

The stellar population in the Galactic halo is characterised by a large fraction of CEMP stars. Most CEMP stars are enriched in $s$-elements (CEMP-$s$ stars), and some of these are also enriched in $r$-elements (CEMP-$s/r$ stars). One formation scenario proposed for CEMP stars invokes wind mass transfer in the past from a TP-AGB primary star to a less massive companion star which is presently observed. We generate low-metallicity populations of binary stars to reproduce the observed CEMP-star fraction. In addition, we aim to constrain our wind mass-transfer model and investigate under which conditions our synthetic populations reproduce observed abundance distributions. We compare the CEMP fractions and the abundance distributions determined from our synthetic populations with observations. Several physical parameters of the binary stellar population of the halo are uncertain, e.g. the initial mass function, the mass-ratio and orbital-period distributions, and the binary fraction. We vary the assumptions in our model about these parameters, as well as the wind mass-transfer process, and study the consequent variations of our synthetic CEMP population. The CEMP fractions calculated in our synthetic populations vary between 7% and 17%, a range consistent with the CEMP fractions among very metal-poor stars recently derived from the SDSS/SEGUE data sample. The results of our comparison between the modelled and observed abundance distributions are different for CEMP-$s/r$ stars and for CEMP-$s$ stars. For the latter, our simulations qualitatively reproduce the observed distributions of C, Na, Sr, Ba, Eu, and Pb. Contrarily, for CEMP-$s/r$ stars our model cannot reproduce the large abundances of neutron-rich elements such as Ba, Eu, and Pb. This result is consistent with previous studies, and suggests that CEMP-$s/r$ stars experienced a different nucleosynthesis history to CEMP-$s$ stars.

preprint2014arXiv

Surprising detection of an equatorial dust lane on the AGB star IRC+10216

Understanding the formation of planetary nebulae remains elusive because in the preceding asymtotic giant branch (AGB) phase these stars are heavily enshrouded in an optically thick dusty envelope. To further understand the morphology of the circumstellar environments of AGB stars we observe the closest carbon-rich AGB star IRC+10216 in scattered light. When imaged in scattered light at optical wavelengths, IRC+10216 surprisingly shows a narrow equatorial density enhancement, in contrast to the large-scale spherical rings that have been imaged much further out. We use radiative transfer models to interpret this structure in terms of two models: firstly, an equatorial density enhancement, commonly observed in the more evolved post-AGB stars, and secondly, in terms of a dust rings model, where a local enhancement of mass-loss creates a spiral ring as the star rotates. We conclude that both models can be used to reproduce the dark lane in the scattered light images, which is caused by an equatorially density enhancement formed by dense dust rather than a bipolar outflow as previously thought. We are unable to place constraints on the formation of the equatorial density enhancement by a binary system.

preprint2014arXiv

Theoretical uncertainties of the type Ia supernova rate

It is thought that type Ia supernovae (SNe Ia) are explosions of carbon-oxygen white dwarfs (CO WDs). Two main evolutionary channels are proposed for the WD to reach the critical density required for a thermonuclear explosion: the single degenerate scenario (SD), in which a CO WD accretes from a non-degenerate companion, and the double degenerate scenario (DD), in which two CO WDs merge. However, it remains difficult to reproduce the observed SN Ia rate with these two scenarios. With a binary population synthesis code we study the main evolutionary channels that lead to SNe Ia and we calculate the SN Ia rates and the associated delay time distributions. We find that the DD channel is the dominant formation channel for the longest delay times. The SD channel with helium-rich donors is the dominant channel at the shortest delay times. Our standard model rate is a factor five lower than the observed rate in galaxy clusters. We investigate the influence of ill-constrained aspects of single- and binary-star evolution and uncertain initial binary distributions on the rate of type Ia SNe. These distributions, as well as uncertainties in both helium star evolution and common envelope evolution, have the greatest influence on our calculated rates. Inefficient common envelope evolution increases the relative number of SD explosions such that for $α_{\rm ce} = 0.2$ they dominate the SN Ia rate. Our highest rate is a factor three less than the galaxy-cluster SN Ia rate, but compatible with the rate determined in a field-galaxy dominated sample. If we assume unlimited accretion onto WDs, to maximize the number of SD explosions, our rate is compatible with the observed galaxy-cluster rate.

preprint2013arXiv

The elusive origin of Carbon-Enhanced Metal-Poor stars

Carbon-enhanced metal-poor (CEMP) stars count for 9-25% of all the very metal-poor stars of the halo. In at least some CEMP stars the chemical enrichment is believed to be due to wind mass transfer in the past from an AGB donor star on to a low-mass companion. However, binary population synthesis models predict much lower CEMP fractions. As an alternative to the canonical Bondi-Hoyle-Lyttleton (BHL) wind accretion model, recent hydrodynamical simulations suggest an efficient mode of wind mass transfer, called wind Roche-lobe overflow (WRLOF), can reproduce observations of AGB winds in binary systems. We use our population synthesis model to test the consequences of WRLOF on a population of CEMP stars. Compared to previous predictions based on the BHL model we find a modest increase of the fraction of CEMP stars and substantial differences in the distributions of carbon and periods in the population of CEMP stars.

preprint2013arXiv

Wind Roche-lobe overflow: Application to carbon-enhanced metal-poor stars

Carbon-enhanced metal-poor stars (CEMP) are observed as a substantial fraction of the very metal-poor stars in the Galactic halo. Most CEMP stars are also enriched in s-process elements and these are often found in binary systems. This suggests that the carbon enrichment is due to mass transfer in the past from an asymptotic giant branch (AGB) star on to a low-mass companion. Models of binary population synthesis are not able to reproduce the observed fraction of CEMP stars without invoking non-standard nucleosynthesis or a substantial change in the initial mass function. This is interpreted as evidence of missing physical ingredients in the models. Recent hydrodynamical simulations show that efficient wind mass transfer is possible in the case of the slow and dense winds typical of AGB stars through a mechanism called wind Roche-lobe overflow (WRLOF), which lies in between the canonical Bondi-Hoyle-Lyttleton (BHL) accretion and Roche-lobe overflow. WRLOF has an effect on the accretion efficiency of mass transfer and on the angular momentum lost by the binary system. The aim of this work is to understand the overall effect of WRLOF on the population of CEMP stars. To simulate populations of low-metallicity binaries we combined a synthetic nucleosynthesis model with a binary population synthesis code. In this code we implemented the WRLOF mechanism. We used the results of hydrodynamical simulations to model the effect of WRLOF on the accretion efficiency and we took the effect on the angular momentum loss into account by assuming a simple prescription. As a result the number of CEMP stars predicted by our model increases by a factor 1.2-1.8 compared to earlier results that consider the BHL prescription. Moreover, higher enrichments of carbon are produced and the final orbital period distribution is shifted towards shorter periods.

preprint2011arXiv

Binary progenitor models of type IIb supernovae

Massive stars that lose their hydrogen-rich envelope down to a few tenths of a solar mass explode as extended type IIb supernovae, an intriguing subtype that links the hydrogen-rich type II supernovae with the hydrogen-poor type Ib and Ic. The progenitors may be very massive single stars that lose their envelope due to their stellar wind, but mass stripping due to interaction with a companion star in a binary system is currently considered to be the dominant formation channel. We computed an extensive grid of binary models with the Eggleton binary evolution code. The predicted rate from our standard models, which assume conservative mass transfer, is about 6 times smaller than the current rate indicated by observations. It is larger but still comparable to the rate expected from single stars. To recover the observed rate we must generously allow for uncertainties and low accretion efficiencies in combination with limited angular momentum loss from the system. Motivated by the claims of detection and non-detection of companions for a few IIb supernovae, we investigate the properties of the secondary star at the moment of explosion. We identify three cases: (1) the companion appears as a hot O star in about 90% of the cases, (2) the companion becomes an over-luminous B star in about 3% of the cases, or (3) the companion will appear as a K supergiant. The second case, which applies to the well-studied case of SN 1993J and possibly to SN 2001ig, is the least common case and requires that the companion very efficiently accretes the transferred material -in contrast to what is required to recover the overall IIb rate. The fast increasing number of type IIb detections may lead to more solid constraints on model uncertainties, like stellar wind-mass loss and internal mixing.

preprint2010arXiv

Chemically homogeneous evolution in massive binaries

Rotation can have severe consequences for the evolution of massive stars. It is now considered as one of the main parameters, alongside mass and metallicity that determine the final fate of single stars. In massive, fast rotating stars mixing processes induced by rotation may be so efficient that helium produced in the center is mixed throughout the envelope. Such stars evolve almost chemically homogeneously. At low metallicity they remain blue and compact, while they gradually evolve into Wolf-Rayet stars and possibly into progenitors of long gamma-ray bursts. In binaries this type of evolution may occur because of (I) tides in very close binaries, as a result of (II) spin up by mass transfer, as result of (III) a merger of the two stars and (IV) when one of the components in the binary was born with a very high initial rotation rate. As these stars stay compact, the evolutionary channels are very different from what classical binary evolutionary models predict. In this contribution we discuss examples of nearly chemically homogeneous evolution in very close tidally-locked binaries. Even in such very close massive binaries, the stars may remain compact and avoid mass transfer, while Roche lobe overflow and a merger would be inevitable in the classical picture. This type of evolution may provide an alternative path to form tight Wolf-Rayet binaries and massive black hole binaries.

preprint2009arXiv

Massive binaries and the enrichment of the interstellar medium in globular clusters

Abundance anomalies observed in globular cluster stars indicate pollution with material processed by hydrogen burning. Two main sources have been suggested: asymptotic giant branch stars and massive stars rotating near the break-up limit. We discuss the potential of massive binaries as an interesting alternative source of processed material. We discuss observational evidence for mass shedding from interacting binaries. In contrast to the fast, radiatively driven winds of massive stars, this material is typically ejected with low velocity. We expect that it remains inside the potential well of a globular cluster and becomes available for the formation or pollution of a second generation of stars. We estimate that the amount of processed low-velocity material that can be ejected by massive binaries is larger than the contribution of two previously suggested sources combined.

preprint2009arXiv

Massive binaries as the source of abundance anomalies in globular clusters

Abundance anomalies observed in globular cluster stars indicate pollution with material processed by hydrogen burning. Two main sources have been suggested: asymptotic giant branch (AGB) stars and massive stars rotating near the break-up limit (spin stars). We propose massive binaries as an alternative source. We compute the evolution of a 20 Msun star in a close binary and find that it sheds about 10 Msun. The ejecta are enriched in He, N, Na, and Al and depleted in C and O, similar to the abundance patterns observed in gobular cluster stars. In contrast to the fast, radiatively driven winds of massive stars, this material is typically ejected with low velocity. We expect that it remains inside the potential well of a globular cluster and becomes available for the formation or pollution of a second generation of stars. We estimate that the amount of processed low-velocity material ejected by massive binaries is greater than the contribution of AGB stars and spin stars combined, assuming that the majority of massive stars in a proto-globular cluster interact. If we take the possible contribution of intermediate mass stars in binaries into account and assume that the ejecta are diluted with an equal amount of unprocessed material, we find that this scenario can provide enough material to form a second generation of low-mass stars, which is as numerous as the first generation of low-mass stars. In this scenario there is no need to make commonly adopted assumptions, such as preferential loss of the first generation of stars, external pollution of the cluster, or an anomalous initial mass function. [Abridged]

preprint2000arXiv

Comprehensive analytic formulae for stellar evolution as a function of mass and metallicity

We present analytic formulae that approximate the evolution of stars for a wide range of mass and metallicity. Stellar luminosity, radius and core mass are given as a function of age, M and Z, for all phases from the zero-age main-sequence up to, and including, the remnant stages. For the most part we find continuous formulae accurate to within 5% of detailed models. These formulae are useful for purposes such as population synthesis that require very rapid but accurate evaluation of stellar properties, and in particular for use in combination with N-body codes. We describe a mass loss prescription that can be used with these formulae and investigate the resulting stellar remnant distribution.

preprint1995arXiv

Approximate input physics for stellar modelling

We present a simple and efficient, yet reasonably accurate, equation of state, which at the moderately low temperatures and high densities found in the interiors of stars less massive than the Sun is substantially more accurate than its predecessor by Eggleton, Faulkner & Flannery. Along with the most recently available values in tabular form of opacities, neutrino loss rates, and nuclear reaction rates for a selection of the most important reactions, this provides a convenient package of input physics for stellar modelling. We briefly discuss a few results obtained with the updated stellar evolution code.