Search

We model cataclysmic variables (CVs) with solar metallicity donors ($X=0.7,\:Z=0.02$) that evolve to form AM CVn stars through the Evolved CV formation channel using various angular momentum loss mechanisms by magnetic braking ($\mathrm{AML_{MB}}$). We find that the time-scale for $\mathrm{AML_{MB}}$ in our double-dynamo (DD) model is shorter than that of previously used empirical formulae. Owing to the shorter time-scales, a larger parameter space of initial conditions evolves to form AM CVn stars with the DD model than with other models. We perform an analysis of the expected number of AM CVn stars formed through the Evolved CV channel and find about $3$ times as many AM CVn stars as reported before. We evolve these systems in detail with the Cambridge stellar evolution code (STARS) and show that evolved CVs populate a region with orbital period $P_\mathrm{orb}\geq5.5\,\mathrm{hr}$. We evolve our donors beyond their orbital period minimum and find that a significant number become extremely H-exhausted systems. This makes them indistinguishable from systems evolved from the He-star and the White Dwarf (WD) channels in terms of the absence of H in their spectra. We also compare the masses, mass-transfer rates of the donor, and the orbital period with observations. We find that the state of the donor and the absence of H in systems such as YZ LMi and V396 Hya match with our modelled trajectories, while systems such as CR Boo and HP Lib match with our modelled tracks if their actual donor mass lies on the lower-end of the observed mass range., Comment: 15 pages, 15 figures. Accepted for publication in MNRAS

We model helium-rich stars with solar metallicity ($X=0.7,\:Z=0.02$) progenitors that evolve to form AM Canum Venaticorum systems through a helium-star formation channel, with the aim to explain the observed properties of Gaia14aae and ZTFJ1637+49. We show that semi-degenerate, H-exhausted ($X\leq 10^{-5}$), He-rich ($Y\approx0.98$) donors can be formed after a common envelope evolution (CEE) phase if either additional sources of energy are used to eject the common envelope, or a different formalism of CEE is implemented. We follow the evolution of such binary systems after the CEE phase using the Cambridge stellar evolution code, when they consist of a He-star and a white dwarf accretor, and report that the mass, radius, and mass-transfer rate of the donor, the orbital period of the system, and the lack of hydrogen in the spectrum of Gaia14aae and ZTFJ1637+49 match well with our modelled trajectories wherein, after the CEE phase Roche lobe overflow is governed not only by the angular momentum loss (AML) owing to gravitational wave radiation ($\mathrm{AML_{GR}}$) but also an additional AML owing to $\alpha-\Omega$ dynamos in the donor. This additional AML is modelled with our double-dynamo (DD) model of magnetic braking in the donor star. We explain that this additional AML is just a consequence of extending the DD model from canonical cataclysmic variable donors to evolved donors. We show that none of our modelled trajectories match with Gaia14aae or ZTFJ1637+49 if the systems are modelled only with $\mathrm{AML_{GR}}$., Comment: 11 pages, 11 figures. Accepted for publication in the Monthly Notices of the Royal Astronomical Society

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., 20 pages including 9 figures and 5 tables; accepted for publication in The Astrophysical Journal (ApJ)

We find that applying a theoretical wind mass-loss rate from Monte Carlo radiative transfer models for hydrogen-deficient stars results in significantly more leftover hydrogen following stable mass transfer through Roche-lobe overflow than when we use an extrapolation of an empirical fit for Galactic Wolf-Rayet stars, for which a negligible amount of hydrogen remains in a large set of binary stellar evolution computations. These findings have implications for modelling progenitors of Type Ib and Type IIb supernovae. Most importantly, our study stresses the sensitivity of the stellar evolution models to the assumed mass-loss rates and the need to develop a better theoretical understanding of stellar winds.

In a series of recent papers, it has been proposed that high field magnetic white dwarfs are the result of close binary interaction and merging. Population synthesis calculations have shown that the origin of isolated highly magnetic white dwarfs is consistent with the stellar merging hypothesis. In this picture, the observed fields are caused by an alpha-Omega dynamo driven by differential rotation. The strongest fields arise when the differential rotation equals the critical break-up velocity and result from the merging of two stars (one of which has a degenerate core) during common envelope evolution or from the merging of two white dwarfs. We now synthesise a population of binary systems to investigate the hypothesis that the magnetic fields in the magnetic cataclysmic variables also originate during stellar interaction in the common envelope phase. Those systems that emerge from common envelope more tightly bound form the cataclysmic variables with the strongest magnetic fields. We vary the common envelope efficiency parameter and compare the results of our population syntheses with observations of magnetic cataclysmic variables. We find that common envelope interaction can explain the observed characteristics of these magnetic systems if the envelope ejection efficiency is low., Comment: 16 pages, 11 figures, accepted for publication in the Monthly Notices of the Royal Astronomical Society

The origin of magnetic fields in isolated and binary white dwarfs has been investigated in a series of recent papers. One proposal is that magnetic fields are generated through an alpha-omega dynamo during common envelope evolution. Here we present population synthesis calculations showing that this hypothesis is supported by observations of magnetic binaries., Comment: 3 pages, 4 figures, Proceedings of the 21st European Workshop on White Dwarfs held July 23-27, 2018 in Austin, TX, USA

Over half of all observed hot subdwarf B (sdB) stars are found in binaries, and over half of these are found in close configurations with orbital periods of 10$ \,\rm{d}$ or less. In order to estimate the companion masses in these predominantly single-lined systems, tidal locking has frequently been assumed for sdB binaries with periods less than half a day. Observed non-synchronicity of a number of close sdB binaries challenges that assumption and hence provides an ideal testbed for tidal theory. We solve the second-order differential equations for detailed 1D stellar models of sdB stars to obtain the tidal dissipation strength and hence to estimate the tidal synchronization time-scale owing to Zahn's dynamical tide. The results indicate synchronization time-scales longer than the sdB lifetime in all observed cases. Further, we examine the roles of convective overshooting and convective dissipation in the core of sdB stars and find no theoretical framework in which tidally-induced synchronization should occur., Comment: 13 pages, 12 figures, accepted for publication by MNRAS

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., 5 pages, 5 figure, 1 table

The Sch\"onberg-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., Comment: 9 pages, 7 figures, published in MNRAS. Updated to (more closely) match published version

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., Comment: 14 pages, 13 figures. Accepted for publication in MNRAS

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., Comment: Accepted for publication in MNRAS. 15 pages, 13 figures

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., Comment: 19 pages, 8 figures, accepted for publication in MNRAS

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., 9 pages, 7 figures, published in MNRAS. Updated to (more closely) match published version

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}, 12 pages, 8 figures, 1 table, accepted for publication in MNRAS

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., Accepted by MNRAS

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., Comment: Accepted for publication in MNRAS, 23 pages, 17 figures and 11 tables. Abstract was editted to fit within arXiv.org submission requirements

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., Comment: 4 pages, 2 figures, to appear in the ASP conference series proceedings of "Advances in Computational Astrophysics: methods, tools, and outcomes" in Cefalu', Italy, June 13-17, 2011

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., Comment: 14 pages, 11 figures. Accepted for publication in MNRAS

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., Comment: Accepted for publication in MNRAS

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., Comment: 12 pages, 11 figures. Published in MNRAS. Very belatedly updated to (more closely) match published version

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., 6 pages, no figures, accepted by MNRAS

WOS: 000280710300026, dWe 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 breakup. 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. Unlike solar-type chromospherically active stars, with enhanced magnetic activity which leads to angular momentum and mass loss, the gainers in classical Algols have radiative envelopes. We further find that normal radiative tides are far too weak to account for the necessary angular momentum loss. Thus enhanced mass loss in a stellar wind seems to be required to spin-down the gainers in classical Algol systems. 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 M(circle dot) and the other with 3.2 and 2 M(circle dot). 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 breakup 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., Scientific and Technological Research Council of Turkey (TUBITAK)Turkiye Bilimsel ve Teknolojik Arastirma Kurumu (TUBITAK) [109T708]; Churchill College; Science and Technology Facilities CouncilScience & Technology Facilities Council (STFC) [ST/H00243X/1, ST/G00269X/1, ST/J001538/1], AD is grateful to the Institute of Astronomy, Cambridge for hospitality and the Scientific and Technological Research Council of Turkey (TUBITAK) for PhD fellowships. CAT thanks Churchill College for his fellowship and is very grateful for the hospitality of IUCAA, Pune and the University of Canterbury where some work on this paper was carried out. This work has been supported in part by TUBITAK through grant no. 109T708. We thank the referee for his or her helpful comments for improving this paper.

International audience; Due to a program error, the published version of Table 2 included inaccurate entries for the calcium, titanium, and mean alpha-element abundance ratios for the stars in the dwarf spheroidal galaxies from Shetrone et al. (2003). A corrected version is given below. We are grateful to Bruno Letarte for detecting this error. We have also updated the information for additional stars in the Sculptor dwarf spheroidal galaxy from D. Geisler et al. (AJ, 128, 1177 [2004]) following the final publication of that paper.

A long standing problem with asymptotic giant branch (AGB) star models has been their inability to produce the low-luminosity carbon stars in the Large and Small Magellanic Clouds. Dredge-up must begin earlier and extend deeper. We find this for the first time in our models of LMC metallicity. Such features are not found in our models of SMC metallicity. The fully implicit and simultaneous stellar evolution code STARS has been used to calculate the evolution of AGB stars with metallicities of Z=0.008 and Z=0.004, corresponding to the observed metallicities of the Large and Small Magellanic Clouds, respecitively. Third dredge-up occurs in stars of 1Msol and above and carbon stars were found for models between 1Msol and 3Msol. We use the detailed models as input physics for a population synthesis code and generate carbon star luminosity functions. We now find that we are able to reproduce the carbon star luminosity function of the LMC without any manipulation of our models. The SMC carbon star luminosity function still cannot be produced from our detailed models unless the minimum core mass for third dredge-up is reduced by 0.06Msol., 6 pages, 5 figures. Accepted for publication in MNRAS

We study the evolution of stars that may be the progenitors of long-soft gamma-ray burst (GRBs) -- rotating naked helium stars presumed to have lost their envelopes to winds or companions. Our aim is to investigate the formation and development of single and binary systems and from this population evaluate the rates of interesting individual species. Using a rapid binary evolution algorithm, that enables us to model the most complex binary systems and to explore the effect of metallicity on GRB production, we draw the following conclusions. First we find that, if we include an approximate treatment of angular momentum transport by mass loss, the resulting spin rates for single stars become too low to form a centrifugally supported disc that can drive a GRB engine. Second massive stars in binaries result in enough angular momentum -- due to spin-orbit tidal interactions -- to form a centrifugally supported disc and are thus capable of supplying a sufficient number of progenitors. This holds true even if only a small fraction of bursts are visible to a given observer and the GRB rate is several hundred times larger than the observed rate. Third low-metallicity stars aid the formation of a rapidly rotating, massive helium cores at collapse and so their evolution is likely to be affected by the local properties of the ISM. This effect could increase the GRB formation rate by a factor of 5-7 at Z=Z_solar/200. Finally we quantify the effects of mass loss, common-envelope evolution and black-hole formation and show that more stringent constraints to many of these evolution parameters are needed in order to draw quantitative conclusions from population synthesis work., 37 pages, 14 figures, MNRAS in press

We have developed a numerical code with which we study the effects of 2D perturbations on stellar structure. We present new numerical and analytical results on the heating of a main-sequence star in a binary system by its companion., Comment: 8 pages, 2 figures, for publication by PASP for conference, '3D Stellar Evolution Workshop'

We present high quality ground-based VRI images of the site of the Type II-P SN1999em (in NGC1637) taken before explosion, which were extracted from the CFHT archive. We determine a precise position of the SN on these images to an accuracy of 0.17''. The host galaxy is close enough (7.5 +/- 0.5 Mpc) that the bright supergiants are resolved as individual objects, however we show that there is no detection of an object at the SN position before explosion that could be interpreted as the progenitor star. By determining the sensitivity limits of the VRI data, we derive bolometric luminosity limits for the progenitor. Comparing these to standard stellar evolutionary tracks which trace evolution up to the point of core carbon ignition, we initially derive an upper mass limit of approximately 12M_sol. However we present evolutionary calculations that follow 7-12M_sol stars throughout their C-burning lifetime and show that we can restrict the mass of the progenitor even further. Our calculations indicate that progenitors initially of 8-10M_sol, undergoing expected mass loss, can also be excluded because a second dredge up sends them to somewhat higher luminosities than a star of initially 12M_sol. These results limit the progenitor's initial main-sequence mass to a very narrow range of 12 +/- 1 M_sol. We discuss the similarities between the Type II-P SNe 1999em and 1999gi and their progenitor mass limits, and suggest that SN Type II-P originate only in intermediate mass stars of 8-12M_sol, which are in the red supergiant region and that higher mass stars produce the other Type II sub-types. (Abridged)., Replaced with accepted version to appear in ApJ, 30 pages, inc. 6 figures

We investigate the gravonuclear instabilities reported by Bono et al. (1997a,b) during the onset of helium-shell burning at the end of the horizontal-branch (HB) phase. These instabilities are characterized by relaxation oscillations within the helium shell which lead to loops in the evolutionary tracks. We find the occurrence of these instabilities depends critically on how the breathing pulses are suppressed near the end of the HB phase. If they are suppressed by omitting the gravitational energy term in the stellar structure equations, then the helium profile within the core at the end of the HB phase will contain a broad region of varying helium abundance. The helium-burning shell which forms in this region is too thick to be unstable, and gravonuclear instabilities do not occur. If, on the other hand, the breathing pulses are suppressed by prohibiting any increase in the central helium abundance, then the final helium profile can exhibit a large discontinuity at the edge of the helium-exhausted core. The helium shell which forms just exterior to this discontinuity is then much thinner and can be thermally unstable. Even in this case, however, the gravonuclear instabilities disappear as soon as the nuclear burning broadens the helium shell into its characteristic S-shape. We conclude that the gravonuclear instabilities found by Bono et al. are a consequence of the ad hoc procedure used to suppress the breathing pulses., Comment: 6 pages, 5 figures, to appear in "The Galactic Halo: From Globular Clusters to Field Stars", Proceedings of the 35th Liege International Astrophysical Colloquium

We present a method of determining lower limits on the masses of pre-main-sequence (PMS) stars and so constraining the PMS evolutionary tracks. This method uses the red-shifted absorption feature observed in some emission-line profiles of T Tauri stars indicative of infall. The maximum velocity of the accreting material measures the potential energy at the stellar surface, which, combined with an observational determination of the stellar radius, yields the stellar mass. This estimate is a lower limit owing to uncertainties in the geometry and projection effects. Using available data, we show that the computed lower limits can be larger than the masses derived from PMS evolutionary tracks for M < 0.5 M_\odot. Our analysis also supports the notion that accretion streams do not impact near the stellar poles but probably hit the stellar surface at moderate latitudes., Accepted by MNRAS. 7 pages, 2 figures

Dynamical models of star clusters are maturing in the sense that effects other than simple point particle dynamics are taken into account. We summarize the relevance of and prospects for this new generation of N-body models., 6 pages including 3 figures, Latex; to appear in Highlights of Astronomy Vol. 11, Kluwer Academic Publishers, ed. Johannes Andersen

Over the last decade or so, we have been developing the possible existence of highly magnetized white dwarfs with analytical stellar structure models. While the primary aim was to explain the nature of the peculiar overluminous type Ia supernovae, later on, these magnetized stars were found to have even wider ranging implications including those for white dwarf pulsars, soft gamma-ray repeaters and anomalous X-ray pulsars, as well as gravitational radiation. In particular, we have explored in detail the mass-radius relations for these magnetized stars and showed that they can be significantly different from the Chandrasekhar predictions which essentially leads to a new super-Chandrasekhar mass-limit. Recently, using the stellar evolution code STARS, we have successfully modelled their formation and cooling evolution directly from the magnetized main sequence progenitor stars. Here we briefly discuss all these findings and conclude with their current status in the scientific community., 13 pages including 7 figures; based on the talk given in the 16th Marcel Grossmann Meeting (MG16) held during July 5-10, 2021; to appear in the proceeding of MG16

More than two dozen soft gamma-ray repeaters (SGRs) and anomalous X-ray pulsars (AXPs) have been detected so far. These are isolated compact objects. Many of them are either found to be associated with supernova remnants or their surface magnetic fields are directly measured, confirming that they are neutron stars (NSs). However, it has been argued that some SGRs and AXPs are highly magnetized white dwarfs (WDs). Meanwhile, the existence of super-Chandrasekhar WDs has remained to be a puzzle. However, not even a single such massive WD has been observed directly. Moreover, some WD pulsars are detected in electromagnetic surveys and some of their masses are still not confirmed. Here we calculate the signal-to-noise ratio for all these objects, considering different magnetic field configurations and thereby estimate the required time for their detection by various gravitational wave (GW) detectors. For SGRs and AXPs, we show that, if these are NSs, they can hardly be detected by any of the GW detectors, while if they are WDs, Big Bang Observer (BBO), DECi-hertz Interferometer Gravitational wave Observatory (DECIGO) and Advanced Laser Interferometer Antenna (ALIA) would be able to detect them within a few days to a year of integration, depending on the magnetic field strength and its configuration. Similarly, if a super-Chandrasekhar WD has a dominant toroidal field, we show that even Laser Interferometer Space Antenna (LISA) and TianQin would be able to detect it within one year of integration. We also discuss how GWs can confirm the masses of the WD pulsars., 10 pages including 12 figures and 2 tables; Accepted for publication in MNRAS

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., 11 pages including 2 in Appendix, 8 figures including 3 in Appendix (total 9 eps files) and 3 tables (change in ordering figures w.r.to previous version); Accepted for publication in MNRAS

We present a detailed study of KIC 2306740, an eccentric double-lined eclipsing binary system. Kepler satellite data were combined with spectroscopic data obtained with the 4.2 m William Herschel Telescope (WHT). This allowed us to determine precise orbital and physical parameters of this relatively long period (P = 10(d).3) and slightly eccentric (e = 0.3) binary system. The physical parameters have been determined to be M-1 = 1.194 +/- 0.008 M-circle dot, M-2 = 1.078 +/- 0.007 M-circle dot, R-1 = 1.682 +/- 0.004 R-circle dot, R-2 = 1.226 +/- 0.005 R-circle dot, L-1 = 2.8 +/- 0.4 L-circle dot, L-2 = 1.8 +/- 0.2 L-circle dot, and orbital separation a = 26.20 +/- 0.04 R-circle dot through simultaneous solutions of the Kepler light curves and WHT radial velocity data. Binarity effects were extracted from the light curve in order to study intrinsic variations in the residuals. Five significant and more than 100 combination frequencies were detected. We modeled the binary system assuming nonconservative evolution models with the Cambridge STARS (TWIN) code, and we show the evolutionary tracks of the components in the log L-logT plane, the log R-log M plane, and the log P-age plane for both spin and orbital periods together with eccentricity e and log R-1. The model of the nonconservative processes in the code led the system to evolve to the observed system parameters in roughly 5.1 Gyr., Turkish Scientific and Research Council [TUBITAK 117F188]; TUBITAK-BIDEB 2211-C scholarship; TUBITAK-BIDEB 2214-A scholarship; Churchill College; COST (European Cooperation in Science and Technology) Action [CA15117, CA16104], We are very grateful to an anonymous referee for comments and helpful, constructive suggestions, which helped us to improve the paper. The authors gratefully acknowledge the numerous people who have helped the NASA Kepler mission possible. This study was supported by the Turkish Scientific and Research Council (TUBITAK 117F188). D.K. is grateful to the Astronomy Department of the University of Geneva (Geneva Observatory) for the kind hospitality during her visit and gratefully acknowledges the support provided by the TUBITAK-BIDEB 2211-C and 2214-A scholarships. C.A.T. thanks Churchill College for his fellowship. K.Y. and D.K. would like to acknowledge the contribution of COST (European Cooperation in Science and Technology) Action CA15117 and CA16104.

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., Comment: 13 pages, 5 figures, submitted to PASA

The combination of asteroseismologically-measured masses with abundances from detailed analyses of stellar atmospheres challenges our fundamental knowledge of stars and our ability to model them. Ancient red-giant stars in the Galactic thick disc are proving to be most troublesome in this regard. They are older than 5 Gyr, a lifetime corresponding to an initial stellar mass of about $1.2{\mathrm{M}_{\odot}}$. So why do the masses of a sizeable fraction of thick-disc stars exceed $1.3{\mathrm{M}_{\odot}}$, with some as massive as $2.3{\mathrm{M}_{\odot}}$ ? We answer this question by considering duplicity in the thick-disc stellar population using a binary population-nucleosynthesis model. We examine how mass transfer and merging affect the stellar mass distribution and surface abundances of carbon and nitrogen. We show that a few per cent of thick-disc stars can interact in binary star systems and become more massive than $1.3{\mathrm{M}_{\odot}}$. Of these stars, most are single because they are merged binaries. Some stars more massive than $1.3{\mathrm{M}_{\odot}}$ form in binaries by wind mass transfer. We compare our results to a sample of the APOKASC data set and find reasonable agreement except in the number of these thick-disc stars more massive than $1.3{\mathrm{M}_{\odot}}$. This problem is resolved by the use of a logarithmically-flat orbital-period distribution and a large binary fraction., Comment: This is a pre-copyedited, author-produced version of an article accepted for publication in Monthly Notices of the Royal Astronomical Society following peer review. The version of record [Robert Izzard, 'Binary stars in the Galactic thick disc', Monthly Notices of the Royal Astronomical Society, accepted 7th Sep. 2017] is at https://academic.oup.com/mnras/article-lookup/doi/10.1093/mnras/stx2355

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., submitted to MNRAS, comments welcome

Generally the virial theorem provides a relation between various components of energy integrated over a system. This helps us to understand the underlying equilibrium. Based on the virial theorem we can estimate, for example, the maximum allowed magnetic field in a star. Recent studies have proposed the existence of highly magnetized white dwarfs, with masses significantly higher than the Chandrasekhar limit. Surface magnetic fields of such white dwarfs could be more than 10^9 G with the central magnitude several orders higher. These white dwarfs could be significantly smaller in size than their ordinary counterparts (with surface fields restricted to about 10^9 G). In this paper we reformulate the virial theorem for non-rotating, highly magnetized white dwarfs (B-WDs) in which, unlike in previous formulations, the contribution of the magnetic pressure to the magnetohydrostatic balance cannot be neglected. Along with the new equation of magnetohydrostatic equilibrium, we approach the problem by invoking magnetic flux conservation and by varying the internal magnetic field with the matter density as a power law. Either of these choices are supported by previous independent work and neither violates any important physics. They are useful while there is no prior knowledge of field profile within a white dwarf. We then compute the modified gravitational, thermal and magnetic energies and examine how the magnetic pressure influences the properties of such white dwarfs. Based on our results we predict important properties of these B-WDs, which turn out to be independent of our chosen field profiles., Comment: 9 pages incluing 8 figures; accepted for publication in MNRAS

Rotation and magnetism are increasingly recognized as important phenomena in stellar evolution. Surface magnetic fields from a few to $20{,}000\,$G have been observed and models have suggested that magnetohydrodynamic transport of angular momentum and chemical composition could explain the peculiar composition of some stars. Stellar remnants such as white dwarfs have been observed with fields from a few to more than $10^{9}\,$G. We investigate the origin of and the evolution, on thermal and nuclear rather than dynamical time-scales, of an averaged large-scale magnetic field throughout a star's life and its coupling to stellar rotation. Large-scale magnetic fields sustained until late stages of stellar evolution with conservation of magnetic flux could explain the very high fields observed in white dwarfs. We include these effects in the Cambridge stellar evolution code using three time-dependant advection-diffusion equations coupled to the structural and composition equations of stars to model the evolution of angular momentum and the two components of the magnetic field. We present the evolution in various cases for a $3\rm\,M_{\odot}$ star from the beginning to the late stages of its life. Our particular model assumes that turbulent motions, including convection, favour small-scale field at the expense of large-scale field. As a result the large-scale field concentrates in radiative zones of the star and so is exchanged between the core and the envelope of the star as it evolves. The field is sustained until the end of the asymptotic giant branch, when it concentrates in the degenerate core.

We evolve stellar models to study the common envelope (CE) interaction of an early asymptotic giant branch star of initial mass $5\,\rm M_{\odot}$ with a companion star of mass ranging from $0.1$ to $2\,\rm M_{\odot}$. We model the CE as a fast stripping phase in which the primary experiences rapid mass loss and loses about 80 per cent of its mass. The post-CE remnant is then allowed to thermally readjust during a Roche-lobe overflow (RLOF) phase and the final binary system and its orbital period are investigated. We find that the post-CE RLOF phase is long enough to allow nuclear burning to proceed in the helium shell. By the end of this phase, the donor is stripped of both its hydrogen and helium and ends up as carbon-oxygen white dwarf of mass about $0.8\,\rm M_{\odot}$. We study the sensitivity of our results to initial conditions of different companion masses and orbital separations at which the stripping phase begins. We find that the companion mass affects the final binary separation and that helium-shell burning causes the star to refill its Roche lobe leading to post-CE RLOF. Our results show that double mass transfer in such a binary interaction is able to strip the helium and hydrogen layers from the donor star without the need for any special conditions or fine tuning of the binary parameters.

We have simulated a common envelope interaction of a 0.88-M$_{\odot}$, 90-R$_{\odot}$, red giant branch star and a 0.6-M$_{\odot}$, compact companion with the smoothed particle hydrodynamics code, Phantom, from the beginning of the Roche lobe overflow phase to the beginning of the self-regulated inspiral, using three different resolutions. The duration of the Roche lobe overflow phase is resolution dependent and would lengthen with increased resolution beyond the $\sim$20 years observed, while the inspiral phase and the post-common envelope separation are largely independent of resolution. Mass transfer rates through the Lagrangian points drive the orbital evolution during the Roche lobe overflow phase, as predicted analytically. The absolute mass transfer rate is resolution dependent, but always within an order of magnitude of the analytical value. Similarly, the gravitational drag in the simulations is close to the analytical approximation. This gives us confidence that simulations approximate reality. The $L_2$ and $L_3$ outflow observed during Roche lobe overflow remains bound, forming a circumbinary disk that is largely disrupted by the common envelope ejection. However, a longer phase of Roche lobe overflow and weaker common envelope ejection typical of a more stable binary may result in a surviving circumbinary disk. Finally, we examine the density distribution resulting from the interaction for simulations that include or omit the phase of Roche lobe overflow. We conclude that the degree of stability of the Roche lobe phase may modulate the shape of the subsequent planetary nebula, explaining the wide range of post-common envelope planetary nebula shapes observed., Submitted to MNRAS, comments and feedback are welcomed

We review the salient evidence consistent with or predicted by the Hoyle-Wickramasinghe (H-W) thesis of Cometary (Cosmic) Biology. Much of this physical and biological evidence is multifactorial. One particular focus are the recent studies which date the emergence of the complex retroviruses of vertebrate lines at or just before the Cambrian Explosion of ∼500 Ma. Such viruses are known to be plausibly associated with major evolutionary genomic processes. We believe this coincidence is not fortuitous but is consistent with a key prediction of H-W theory whereby major extinction-diversification evolutionary boundaries coincide with virus-bearing cometary-bolide bombardment events. A second focus is the remarkable evolution of intelligent complexity (Cephalopods) culminating in the emergence of the Octopus. A third focus concerns the micro-organism fossil evidence contained within meteorites as well as the detection in the upper atmosphere of apparent incoming life-bearing particles from space. In our view the totality of the multifactorial data and critical analyses assembled by Fred Hoyle, Chandra Wickramasinghe and their many colleagues since the 1960s leads to a very plausible conclusion - life may have been seeded here on Earth by life-bearing comets as soon as conditions on Earth allowed it to flourish (about or just before 4.1 Billion years ago); and living organisms such as space-resistant and space-hardy bacteria, viruses, more complex eukaryotic cells, fertilised ova and seeds have been continuously delivered ever since to Earth so being one important driver of further terrestrial evolution which has resulted in considerable genetic diversity and which has led to the emergence of mankind. ispartof: PROGRESS IN BIOPHYSICS & MOLECULAR BIOLOGY vol:136 pages:3-23 ispartof: location:England status: published

Stars with unusual elemental abundances offer clues about rare astrophysical events or nucleosynthetic pathways. Stars with significantly depleted magnesium and enhanced potassium ([Mg/Fe] < -0.5; [K/Fe] > 1) have to date only been found in the massive globular cluster NGC 2419 and, to a lesser extent, NGC 2808. The origin of this abundance signature remains unknown, as does the reason for its apparent exclusivity to these two globular clusters. Here we present 112 field stars, identified from 454,180 LAMOST giants, that show significantly enhanced [K/Fe] and possibly depleted [Mg/Fe] abundance ratios. Our sample spans a wide range of metallicities (-1.5 < [Fe/H] < 0.3), yet none show abundance ratios of [K/Fe] or [Mg/Fe] that are as extreme as those observed in NGC 2419. If confirmed, the identified sample of stars represents evidence that the nucleosynthetic process producing the anomalous abundances ratios of [K/Fe] and [Mg/Fe] probably occurs at a wide range of metallicities. This would suggest that pollution scenarios that are limited to early epochs (such as Population III supernovae) are an unlikely explanation, although they cannot be ruled out entirely. This sample is expected to help guide modelling attempts to explain the origin of the Mg-K abundance signature.

Wickramasinghe et al. (2014) and Briggs et al. (2015) have proposed that the strong magnetic fields observed in some single white dwarfs (MWDs) are formed by a dynamo driven by differential rotation when two stars, the more massive one with a degenerate core, merge during common envelope (CE) evolution (Ferrario et al., 2015b). We synthesize a population of binaries to investigate if fields in the magnetic cataclysmic variables (MCVs) may also originate during stellar interaction in the CE phase., 2 pages, 1 figure, accepted by Contrib. Astron. Obs. Skalnate Pleso. Proceedings of the conference "Stars with a stable magnetic field: from pre-main sequence to compact remnants", Brno, 2017. Editors: E. Paunzen, M. Netopil, V. Petit

Making use of a new N-body model to describe the evolution of a moderate-size globular cluster we in- vestigate 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.

Diffusion of elements in the atmosphere of a star can drastically affect its surface composition, sometimes leading to unusual mixtures. These chemically peculiar stars can be identified fromthe presence of unusual lines in their spectra. Some hot subdwarf stars show extraordinary abundances of elements such as lead, zirconium and strontium, while the abundance of helium ranges from practically zero to almost 100 percent across the hot subdwarf population. A sequence of extreme horizontal branch star models was generated by producing a number of post-common envelope objects from red giants. The evolution of these subdwarf progenitors was computed with the MESA stellar evolution code from immediately after envelope ejection right up to the ignition of helium in the core. Envelope abundances were calculated at the zero age horizontal branch for models both with and without the presence of diffusion. A small number of simulations also looked at the effects on radiative levitation of these abundances, to test how well diffusion physics is able to reproduce observational data.

A dynamo mechanism driven by differential rotation when stars merge has been proposed to explain the presence of strong fields in certain classes of magnetic stars. In the case of the high field magnetic white dwarfs (HFMWDs), the site of the differential rotation has been variously thought to be the common envelope, the hot outer regions of a merged degenerate core or an accretion disc formed by a tidally disrupted companion that is subsequently accreted by a degenerate core. We have shown previously that the observed incidence of magnetism and the mass distribution in HFMWDs are consistent with the hypothesis that they are the result of merging binaries during common envelope evolution. Here we calculate the magnetic field strengths generated by common envelope interactions for synthetic populations using a simple prescription for the generation of fields and find that the observed magnetic field distribution is also consistent with the stellar merging hypothesis. We use the Kolmogorov-Smirnov test to study the correlation between the calculated and the observed field strengths and find that it is consistent for low envelope ejection efficiency. We also suggest that field generation by the plunging of a giant gaseous planet on to a white dwarf may explain why magnetism among cool white dwarfs (including DZ white dwarfs) is higher than among hot white dwarfs. In this picture a super Jupiter residing in the outer regions of the planetary system of the white dwarf is perturbed into a highly eccentric orbit by a close stellar encounter and is later accreted by the white dwarf., 8 pages, 5 figures, accepted for publication in MNRAS

Convection in the cores of massive stars becomes anisotropic when they rotate. This anisotropy leads to a misalignment of the thermal gradient and the thermal flux, which in turn results in baroclinicity and circulation currents in the upper radiative zone. We show that this induces a much stronger meridional flow in the radiative zone than previously thought. This drives significantly enhanced mixing, though this mixing does not necessarily reach the surface. The extra mixing takes on a similar form to convective overshooting, and is relatively insensitive to the rotation rate above a threshold, and may help explain the large overshoot distances inferred from observations. This has significant consequences for the evolution of these stars by enhancing core-envelope mixing., Comment: 21 pages, 10 figures. Published in MNRAS. Corrected typographic error in text between equations 82 and 83

We present an approach to turbulence closure based on mixing length theory with three-dimensional fluctuations against a two-dimensional background. This model is intended to be rapidly computable for implementation in stellar evolution software and to capture a wide range of relevant phenomena with just a single free parameter, namely the mixing length. We incorporate magnetic, rotational, baroclinic and buoyancy effects exactly within the formalism of linear growth theories with nonlinear decay. We treat differential rotation effects perturbatively in the corotating frame using a novel controlled approximation which matches the time evolution of the reference frame to arbitrary order. We then implement this model in an efficient open source code and discuss the resulting turbulent stresses and transport coefficients. We demonstrate that this model exhibits convective, baroclinic and shear instabilities as well as the magnetorotational instability (MRI). It also exhibits non-linear saturation behaviour, and we use this to extract the asymptotic scaling of various transport coefficients in physically interesting limits., Comment: 17 pages, 20 figures, published in MNRAS. Figures 9-19 and C1 and associated text have been updated relative to MNRAS version to reflect minor changes due to corrections to the numerical code. Typographic errors have been corrected in equations 8, 10 and 11 and Figure 1 caption. Figure 3 legend updated to be clearer, and colors have been coordinated between Figures 1-3