Search

Over the last few decades, there has been considerable interest in the violation of the sacred "Chandrasekhar" mass limit of white dwarfs (WDs). Peculiar over-luminous type Ia supernovae (such as SNLS-03D3bb) lend observational support to the idea that these super-Chandrasekhar WDs exist. Our group, for more than a decade, has been actively working on the theoretical possibility of these objects through the presence of the star's magnetic field. The magnetic field greatly contributes to the existence of these massive WDs, both through classical and quantum effects. In this work, we explore super-Chandrasekhar WDs, formed via evolution from a main sequence star, as a result of the classical effects of the star's magnetic field. We obtain super-Chandrasekhar WDs and new mass limit(s), depending on the magnetic field geometry. We explore the full evolution and stability of these objects from the main sequence stage through the one-dimensional stellar evolution code STARS. In order to do so, we have appropriately modified the given codes by introducing magnetic effect and cooling. Our simulation confirms that massive WDs are possible in the presence of a magnetic field satisfying underlying stability., Comment: 12 pages, 4 figures; To be published in Astrophysics and Space Science Proceedings, titled "The Relativistic Universe: From Classical to Quantum, Proceedings of the International Symposium on Recent Developments in Relativistic Astrophysics", Gangtok, December 11-13, 2023: to felicitate Prof. Banibrata Mukhopadhyay on his 50th Birth Anniversary", Editors: S Ghosh & A R Rao, Springer Nature

Binary stars are responsible for many unusual astrophysical phenomena, including some important explosive cosmic events. The stability criteria for rapid mass transfer and common-envelope evolution are fundamental to binary star evolution. They determine the mass, mass ratio, and orbital distribution of systems such as X-ray binaries and merging gravitational-wave sources. We use our adiabatic mass-loss model to systematically survey metal-poor and solar-metallicity donor thresholds for dynamical timescale mass transfer. The critical mass ratios qad are systematically explored, and the impact of metallicity and nonconservative mass transfer are studied. For metal-poor radiative-envelope donors, qad are smaller than those for solar-metallicity stars at the same evolutionary stage. However, qad do the opposite for convective-envelope donors. Nonconservative mass transfer significantly decreases qad for massive donors. This is because it matters how conservative mass transfer is during the thermal timescale phase immediately preceding a delayed dynamical mass transfer. We apply our theoretical predictions to observed high-mass X-ray binaries that have overfilled their Roche lobes and find a good agreement with their mass ratios. Our results can be applied to study individual binary objects or large samples of binary objects with binary population synthesis codes., Comment: Submitted to ApJ. Comments are welcome

We present the observational mass-radius (M-R) relation for a sample of 47 magnetized white dwarfs (WDs) with the magnetic field strength (B) ranging from 1 to 773 MG, identified from the SDSS data release 7 (DR7). We derive their effective temperature, surface gravity (log g), luminosity, radius, and mass. While atmospheric parameters are derived using a Virtual Observatory Spectral Energy Distribution Analyzer (VOSA), the mass is derived using their location in the HR diagram in comparison with the evolutionary tracks of different masses. We implement this mass measurement instead of a more traditional method of deriving masses from log g, which is unreliable as is based on SED and generates errors from other physical parameters involved. The main disadvantage of this method is that we need to assume a core composition of WDs. As it is complicated to identify the exact composition of these WDs from low-resolution spectra, we use tracks for the masses 0.2 to 0.4 solar mass assuming a He-core, 0.5 to 1.0 solar mass assuming CO core, and above solar mass assuming O-Ne-Mg core. We compare the observed M-R relation with those predicted by the finite temperature model by considering different B, which are well in agreement considering their relatively low surface fields, less than or of the order of 10^9 G. Currently, there is no direct observational detection of magnetized WDs with B > 10^9 G. We propose that our model can be further extrapolated to higher B, which may indicate the existence of super-Chandrasekhar mass (M > 1.4 solar mass) WDs at higher B., Comment: 9 latex pages including 10 figures (11 pdf figures); accepted for publication in MNRAS

We extend our model of magnetic braking (MB), driven by an $\alpha-\Omega$ dynamo mechanism, from fully convective M-dwarfs (FCMDs) to explain the surface and internal spin $P_\mathrm{spin}$ evolution of partly convective dwarfs (PCDs) starting from the disc-dispersal stage to the main-sequence turnoff. In our model, the spin of the core is governed by shear at the core-envelope boundary while the spin of the envelope is governed by MB and shear. We show that (1) the most massive FCMDs experience a stronger spin-down than PCDs and less massive FCMDs, (2) the stalled spin-down and enhanced activity of K-dwarfs and the pileup of G-dwarfs older than a few Gyr are stellar-structure- and MB-dependent, and weakly dependent on core-envelope coupling effects, (3) our expression of the core-envelope convergence time-scale $\tau_\mathrm{converge}(M_\ast,\,P_\mathrm{spin})$ between a few 10 to 100~Myr strongly depends on stellar structure but weakly on MB strength and shear, such that fast and massive rotators achieve corotation earlier, (4) our estimates of the surface magnetic fields are in general agreement with observations and our wind mass loss evolution explains the weak winds from the solar analog $\pi^1$ UMa and (5) with our model the massive young Sun hypothesis as a solution to the faint young Sun problem can likely be ruled out, because the maximum mass lost by winds from our Sun with our model is about an order of magnitude smaller than required to solve the problem., Comment: Accepted for publication in MNRAS. Minor changes from the original pre-print

Period bouncers are cataclysmic variables (CVs) that have evolved past their orbital period minimum. The strong disagreement between theory and observations of the relative fraction of period bouncers is a severe shortcoming in the understanding of CV evolution. We test the implications of the hypothesis that magnetic braking (MB), which is suggested to be an additional angular momentum loss (AML) mechanism for CVs below the period gap ($P_\mathrm{orb}\lesssim 120$ min), weakens around their period minimum. We compute the evolution of CV donors below the period gap using the MESA code, assuming that the evolution of the system is driven by AML by gravitational wave radiation (GWR) and MB. We parametrize the MB strength as $\mathrm{AML_{MB}}=\kappa\mathrm{AML_{GWR}}$. We compute two qualitatively different sets of models, one where $\kappa$ is a constant and the other where $\kappa$ depends on stellar parameters. We find that two crucial effects drive the latter set of models. (1) A decrease in $\kappa$ as CVs approach the period minimum stalls their evolution so that they spend a long time in the observed period minimum spike ($80\lesssim P_\mathrm{orb}/\,\mathrm{min}\lesssim 86$). Here, they become difficult to distinguish from pre-bounce systems in the spike. (2) A strong decrease in the mass-transfer rate makes them virtually undetectable as they evolve further. So, the CV stalls around the period minimum and then `disappears'. This reduces the number of detectable bouncers. Physical processes, such as dynamo action, white dwarf magnetism, and dead zones, may cause such a weakening of MB at short orbital periods. The weakening magnetic braking formalism provides a possible solution to the problem of the lack of period bouncers in CV observational surveys., Comment: Accepted for publication in A&A Letters. Significant revision from the original version

Accreting main-sequence stars expand significantly when the mass accretion timescale is much shorter than their thermal timescales. This occurs during mass transfer from an evolved giant star onto a main-sequence companion in a binary system, and is an important phase in the formation of compact binaries including X-ray binaries, cataclysmic variables, and gravitational-wave sources. In this study, we compute 1D stellar models of main-sequence accretors with different initial masses and accretion rates. The calculations are used to derive semi-analytical approximations to the maximum expansion radius. We assume that mass transfer remains fully conservative as long as the inflated accretor fits within its Roche lobe, leading stars to behave like hamsters, stuffing excess material behind their expanding cheeks. We suggest a physically motivated prescription for the mass growth of such "hamstars", which can be used to determine mass-transfer efficiency in rapid binary population synthesis models. With this prescription, we estimate that progenitors of high-mass X-ray binaries and gravitational-wave sources may have experienced highly non-conservative mass transfer. In contrast, for low-mass accretors, the accretion timescale can exceed the thermal timescale by a larger factor without causing significant radial expansion., Comment: 12 pages, 6 figures. Accepted for publication in ApJ Letters

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., Comment: 20 pages, 12 figures. Refereed manuscript accepted for publication in PoS-SISSA. Based on an invited talk at the conference titled 'The Golden Age of Cataclysmic Variables and Related Objects VI' held from 4 to 9 September 2023 in Palermo, Italy

The common envelope evolution (CEE) is vital in forming short orbital period compact binaries. It covers many objects, such as double compact merging binaries, type Ia supernovae progenitors, binary pulsars, and X-ray binaries. Knowledge about the common envelope (CE) eject efficiency still needs to be improved, though progress has been made recently. Short orbital period hot subdwarf B star plus white dwarf binaries are the most straightforward samples to constrain CEE physics. We apply the known orbital period-white dwarf relation to constrain the sdB progenitor of seven sdB+WD binaries with a known inclination angle. The average value of the CE efficiency parameter is 0.32, which is consistent with previous studies. However, the CE efficiency might not be a constant but is a function of the initial mass ratio based on well-constrained sdB progenitor mass and evolutionary stage. Our results can be used as physical inputs for binary population synthesis simulations on related objects. A similar method can also be applied to study other short orbital period WD binaries., Comment: 16 pages, 12 figures, 2 tables, accepted for publication in The Astrophysical Journal

We extend a magnetic braking (MB) model, which has been used earlier to address the evolution of cataclysmic variables, to address the spin period $P_\mathrm{spin}$ evolution of fully convective M dwarf (FCMD) stars. The MB mechanism is an $\alpha-\Omega$ dynamo, which leads to stellar winds that carry away angular momentum. We model our MB torque such that the FCMDs experience a MB torque, approximately scaling as $P_\mathrm{spin}^{-1}$ at shorter periods, before transitioning into a Skumanich-type MB torque, scaling as $P_\mathrm{spin}^{-3}$. We also implement a parametrized reduction in the wind mass loss owing to the entrapment of winds in dead zones. We choose a set of initial conditions and vary the two free parameters in our model to find a good match of our spin trajectories with open clusters containing FCMDs such as NGC2547, Pleiades, NGC2516 and Praesepe. We find that our model can explain the long spin periods of field stars and that a spread in spin distribution persists till over 3 Gyr. An advantage of our model is in relating physically motivated estimations of the magnetic field strength and stellar wind to properties of the stellar dynamo, which other models often remain agnostic about. We track the spin dependence of the wind mass losses, Alfv\'en radii and surface magnetic fields and find good agreement with observations. We discuss the implications of our results on the effect of the host FCMD on any orbiting exoplanets and our plans to extend this model to explain solar-like stars in the future., Comment: 16 pages, 11 figures. Accepted for publications in MNRAS

We use the rapid binary stellar evolution code BINARY_C to estimate the rate of merging neutron stars with numerous combinations of envelope ejection efficiency and natal kick dispersion. We find a peak in the local rate of merging neutron stars around $\alpha \approx 0.3$$-$$0.4$, depending on the metallicity, where $\alpha$ is the efficiency of utilising orbital energy to unbind the envelope. The peak height decreases with increasing electron-capture supernova kick dispersion $\sigma_\mathrm{ECSN}$. We explain the peak as a competition between the total number of systems that survive the common-envelope phase increasing with $\alpha$ and their separation, which increases with $\alpha$ as well. Increasing $\alpha$ reduces the fraction of systems that merge within a time shorter than the age of the Universe and results in different mass distributions for merging and non-merging double neutron stars. This offers a possible explanation for the discrepancy between the Galactic double neutron star mass distribution and the observed massive merging neutron star event GW190425. Within the $\alpha$$-$$\sigma_\mathrm{ECSN}$ parameter space that we investigate, the rate of merging neutron stars spans several orders of magnitude up to more than $1\times 10^{3} \, \mathrm{Gpc}^{-3}\,\mathrm{yr}^{-1}$ and can be higher than the observed upper limit or lower than the observed lower limit inferred thus far from merging neutron stars detected by gravitational waves. Our results stress the importance of common-envelope physics for the quantitative prediction and interpretation of merging binary neutron star events in this new age of gravitational wave astronomy., Comment: 16 pages, 10 figures, this is an author-produced version of an article accepted for publication in MNRAS following peer review. The version of record is now available online at https://doi.org/10.1093/mnras/stad971 . (N.b. the tables are better formatted in the arXiv article.) Our data are available online at https://zenodo.org/record/7811486#.ZDJt6I5BzJU

The stability criteria of rapid mass transfer and common envelope evolution are fundamental in binary star evolution. They determine the mass, mass ratio and orbital distribution of many important systems, such as X-ray binaries, Type Ia supernovae and merging gravitational wave sources. We use our adiabatic mass-loss model to systematically survey the intermediate-mass stars' thresholds for dynamical-timescale mass transfer. The impact of metallicity on the stellar responses and critical mass ratios is explored. Both tables ($Z=0.001$) and fitting formula ($Z=0.001$ and $Z=0.02$) of critical mass ratios of intermediate-mass stars are provided. An application of our results to intermediate-mass X-ray binaries (IMXBs) is discussed. We find that the predicted upper limit to mass ratios, as a function of orbital period, is consistent with the observed IMXBs that undergo thermal or nuclear timescale mass transfer. According to the observed peak X-ray luminosity $L_\mathrm{X}$, we predict the range of $L_\mathrm{X}$ for IMXBs as a function of the donor mass and the mass transfer timescale., Comment: 16 pages, 19 figures, 3 tables, accepted for publication in The Astrophysical Journal

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 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., Comment: 22 pages, 17 Figures, Accepted for publication 13/05/22

Common envelope evolution (CEE) physics plays a fundamental role in the formation of binary systems, such as mergering stellar gravitational wave sources, pulsar binaries and type Ia supernovae. A precisely constrained CEE has become more important in the age of large surveys and gravitational wave detectors. We use an adiabatic mass loss model to explore how the total energy of the donor changes as a function of the remnant mass. This provides a more self-consistent way to calculate the binding energy of the donor. For comparison, we also calculate the binding energy through integrating the total energy from the core to the surface. The outcome of CEE is constrained by total energy conservation at the point at which both component's radii shrink back within their Roche lobes. We apply our results to 142 hot subdwarf binaries. For shorter orbital period sdBs, the binding energy is highly consistent. For longer orbital period sdBs in our samples, the binding energy can differ by up to a factor of 2. The CE efficiency parameter $\beta_\mathrm{CE}$ becomes smaller than $\alpha_\mathrm{CE}$ for the final orbital period $\log_{10} P_{\mathrm{orb}}/\mathrm{d} > -0.5$. We also find the mass ratios $\log_{10} q$ and CE efficiency parameters $\log_{10} \alpha_{\mathrm{CE}}$ and $\log_{10} \beta_{\mathrm{CE}}$ linearly correlate in sdBs, similarly to De Marco et al. (2010) for post-AGB binaries., Comment: 18 pages, 17 figures, 3 tables, Accepted for publication in ApJ, minor updates for the final proof

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

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., Comment: 10 pages including 12 figures and 2 tables; Accepted for publication in MNRAS

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

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

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\Omega|$ of order the angular frequency $\Omega$ for slowly-rotating systems with $\Omega \ll |N|$, where $N$ is the \brvs\ frequency, and find that it declines as a power-law in $\Omega$ for rapidly-rotating systems with $\Omega \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., Comment: 26 pages, 3 figures. Accepted in MNRAS. Comments welcome

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., Comment: 27 pages. 15 figures. Accepted in MNRAS. Comments welcome

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., Comment: 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

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

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., Comment: Published in MNRAS

Asteroseismological analysis of NY Vir suggests that at least the outer 55 per cent of the star (in radius) rotates as a solid body and is tidally synchronized to the orbit. Detailed calculation of tidal dissipation rates in NY Vir fails to account for this synchronization. Recent observations of He core burning stars suggest that the extent of the convective core may be substantially larger than that predicted with theoretical models. We conduct a parametric investigation of sdB models generated with the Cambridge STARS code to artificially extend the radial extent of the convective core. These models with extended cores still fail to account for the synchronization. Tidal synchronization may be achievable with a non-MLT treatment of convection., Comment: 7 Pages, 7 Figures, 1 Table, Submitted 21/12/2018, Accepted 20/02/2019 MNRAS

Theoretical models of stellar evolution predict that most of the lithium inside a star is destroyed as the star becomes a red giant. However, observations reveal that about 1% of red giants are peculiarly rich in lithium, often exceeding the amount in the interstellar medium or predicted from the Big Bang. With only about 150 lithium-rich giants discovered in the past four decades, and no distinguishing properties other than lithium enhancement, the origin of lithium-rich giant stars is one of the oldest problems in stellar astrophysics. Here we report the discovery of 2,330 low-mass (1 to 3$\,M_\odot$) lithium-rich giant stars, which we argue are consistent with internal lithium production that is driven by tidal spin-up by a binary companion. Our sample reveals that most lithium-rich giants have helium-burning cores ($80^{+7}_{-6}\%$), and that the frequency of lithium-rich giants rises with increasing stellar metallicity. We find that while planet accretion may explain some lithium-rich giants, it cannot account for the majority that have helium-burning cores. We rule out most other proposed explanations as the primary mechanism for lithium-rich giants, including all stages related to single star evolution. Our analysis shows that giants remain lithium-rich for only about two million years. A prediction from this lithium depletion timescale is that most lithium-rich giants with a helium-burning core have a binary companion., Comment: ApJ accepted. Catalog of 2,330 lithium-rich giants available at https://zenodo.org/record/3515777

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

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

We argue that an increased temperature in star-forming clouds alters the stellar initial mass function to be more bottom-light than in the Milky Way. At redshifts $z \gtrsim 6$, heating from the cosmic microwave background radiation produces this effect in all galaxies, and it is also present at lower redshifts in galaxies with very high star formation rates (SFRs). A failure to account for it means that at present, photometric template fitting likely overestimates stellar masses and star formation rates for the highest-redshift and highest-SFR galaxies. In addition this may resolve several outstanding problems in the chemical evolution of galactic halos., Comment: 9 pages, 5 figures. Published in MNRAS. Added further references

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., Comment: Submitted to MNRAS, comments and feedback are welcomed

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., Comment: Accepted to MNRAS 10 pages, 8 figures

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

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

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.

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

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., Comment: 8 pages, 5 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., Comment: 13 pages, 16 figures, accepted by MNRAS

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

Diffusion of elements in the atmosphere and envelope of a star can drastically alter its surface composition, leading to extreme chemical peculiarities. We consider the case of hot subdwarfs, where surface helium abundances range from practically zero to almost 100 percent. Since hot subdwarfs can form via a number of different evolution channels, a key question concerns how the formation mechanism is connected to the present surface chemistry. A sequence of extreme horizontal branch star models was generated by producing post-common envelope stars from red giants. Evolution was computed with MESA from envelope ejection up to core-helium ignition. Surface abundances were calculated at the zero-age horizontal branch for models with and without diffusion. A number of simulations also included radiative levitation. The goal was to study surface chemistry during evolution from cool giant to hot subdwarf and determine when the characteristic subdwarf surface is established. Only stars leaving the giant branch close to core-helium ignition become hydrogen-rich subdwarfs at the zero-age horizontal branch. Diffusion, including radiative levitation, depletes the initial surface helium in all cases. All subdwarf models rapidly become more depleted than observations allow. Surface abundances of other elements follow observed trends in general, but not in detail. Additional physics is required., Comment: 11 pages, 8 figures, 2 tables, accepted for publication in MNRAS

The most frequently seen binary companions to millisecond pulsars (MSPs) are helium white dwarfs (He WDs). The standard rejuvenation mechanism, in which a low- to intermediate-mass companion to a neutron star fills its Roche lobe between central hydrogen exhaustion and core helium ignition, is the most plausible formation mechanism. We have investigated whether the observed population can realistically be formed via this mechanism. We used the Cambridge STARS code to make models of Case B RLOF with Reimers' mass loss from the donor. We find that the range of initial orbital periods required to produce the currently observed range of orbital periods of MSPs is extremely narrow. To reduce this fine tuning, we introduce a companion enhanced wind (CEW) that strips the donor of its envelope more quickly so that systems can detach at shorter periods. Our models indicate that the fine tuning can be significantly reduced if a CEW is active. Because significant mass is lost owing to a CEW we expect some binary pulsars to accrete less than the 0.1Msun needed to spin them up to millisecond periods. This can account for mildly recycled pulsars present along the entire Mc-Porb relation. Systems with Pspin > 30ms are consistent with this but too few of these mildly recycled pulsars have yet been observed to make a significant comparison., Comment: 10 pages, 9 figures

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

We study the interaction between stellar irradiation and tidal heating in gaseous planets with short orbital periods. The intentionally simplified atmospheric model we employ makes the problem analytically tractable and permits the derivation of useful scaling relations. We show that many tidal models provide thermal feedback, producing interior radiative zones and leading to enhanced g-mode dissipation with a wide spectrum of resonances. These resonances are dynamically tuned by the thermal feedback, and so represent a novel form of thermomechanical feedback, coupling vibrational modes to the very slow thermal evolution of the planet. We then show that stellar irradiation allows the heat produced by these modes to be trapped at depth with high efficiency, leading to entropy increase in the central convective region, as well as expansion of the planet's radius sufficient to match observed swelling. We find that thermally driven winds play an essential role in this process by making the thermal structure of the atmosphere spherically symmetric within a few scale heights of the photosphere. We characterise the relationship between the swelling factor, the orbital period and the host star and determine the timescale for swelling. We show that these g-modes suffice to produce bloating on the order of the radius of the planet over $\mathrm{Gyr}$ timescales when combined with significant insolation and we provide analytic relations for the relative magnitudes of tidal heating and insolation., Comment: 15 pages, 5 figures, accepted to MNRAS

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., Comment: 10 pages, 4 figures, submitted to PASA May 17 2016, resubmitted June 30 2016

The candidate Thorne-\.Zytkow object (T\.ZO), HV2112, is becoming a well-studied if enigmatic object. A key point of its candidacy as a T\.ZO 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\.ZO., Comment: 5 pages, 5 figure, 1 table

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., Comment: 12 Pages, 6 figures, accepted for publication in the MNRAS

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., Comment: 12 pages, 7 figures

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 $\alpha$. 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., Comment: Published in MNRAS. 12 pages, 8 figures. Minor changes to match published version

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., Comment: 5 pages, submitted to MNRAS letters. Comments welcome