Showing 1–38 of 38 results for author: Belloni, D

Double white dwarf binary population in MOCCA star clusters -- Comparisons with observations of close and wide binaries

Authors: Lucas Hellström, Mirosław Giersz, Arkadiusz Hypki, Diogo Belloni, Abbas Askar, Grzegorz Wiktorowicz

Abstract: There could be a significant population of double white dwarf binaries (DWDs) inside globular clusters (GCs), however, these are often too faint to be individually observed. We have utilized a large number GC models evolved with the Monte Carlo Cluster Simulator (MOCCA) code, to create a large statistical dataset of DWDs. These models include multiple-stellar populations, resulting in two distinct… ▽ More There could be a significant population of double white dwarf binaries (DWDs) inside globular clusters (GCs), however, these are often too faint to be individually observed. We have utilized a large number GC models evolved with the Monte Carlo Cluster Simulator (MOCCA) code, to create a large statistical dataset of DWDs. These models include multiple-stellar populations, resulting in two distinct initial populations: one dense and another less dense. Due to the lower density of one population, a large number of objects escape during the early GC evolution, leading to a high mass-loss rate. In this dataset we have analysed three main groups of DWDs, namely in-cluster binaries, escaped binaries, and isolated evolution of primordial binaries. We compared the properties of these groups to observations of close and wide binaries. We find that the number of escaping DWDs is significantly larger than the number of in-cluster binaries and those that form via the isolated evolution of all promiridial binaries in our GC models. This suggests that dynamics play an important role in the formation of DWDs. For close binaries, we found a good agreement in the separations of escaped binaries and isolated binaries, but in-cluster binaries showed slight differences. We could not reproduce the observed extremely low mass WDs due to the limitations of our stellar and binary evolution prescriptions. For wide binaries, we also found a good agreement in the separations and masses, after accounting for observational selection effects. We conclude that, even though the current observational samples of DWDs are extremely biased and incomplete, our results compare reasonably well with observations. △ Less

Submitted 9 May, 2024; v1 submitted 7 May, 2024; originally announced May 2024.

Comments: 15 pages, 9 figures

Reevaluating the Origin of Detectable Cataclysmic Variables in Globular Clusters: Testing the Importance of Dynamics

Authors: Liliana Rivera Sandoval, Diogo Belloni, Miriam Ramos Arevalo

Abstract: Based on the current detectable cataclysmic variable (CV) population in Galactic globular clusters (GCs), we show that there is not a clear relation between the number of sources per unit of mass and the stellar encounter rate, the cluster mass, or the cluster central density. If any, only in the case of core-collapsed GCs could there be an anticorrelation with the stellar encounter rate. Our find… ▽ More Based on the current detectable cataclysmic variable (CV) population in Galactic globular clusters (GCs), we show that there is not a clear relation between the number of sources per unit of mass and the stellar encounter rate, the cluster mass, or the cluster central density. If any, only in the case of core-collapsed GCs could there be an anticorrelation with the stellar encounter rate. Our findings contrast with previous studies where clear positive correlations were identified. Our results suggest that correlations between faint X-ray sources, from which often conclusions for the CV population are drawn, and the GC parameters considered here, are likely influenced by other type of X-ray sources, including other types of compact binaries, which have X-ray luminosities similar to CVs. The findings presented here also suggest that the role of primordial systems is more important than previously believed and that dynamical formation has less influence in the current detectable CV population. The long-standing paradigm that GCs are efficient factories of CVs formed via dynamical interactions does not seem to be supported by current observations. △ Less

Submitted 4 April, 2024; originally announced April 2024.

Comments: Published on ApJ Letters on March 27, 2024

Journal ref: The Astrophysical Journal Letters, Volume 964, Number 2, 2024

The formation of the magnetic symbiotic star FN Sgr

Authors: Diogo Belloni, Joanna Mikołajewska, Matthias R. Schreiber

Abstract: To shed light on the origin of magnetic symbiotic stars, we investigated the system FN Sgr in detail. We searched for a reasonable formation pathway to explain its stellar and binary parameters including the magnetic field of the accreting white dwarf. We used the MESA code to carry out pre-CE and post-CE binary evolution and determined the outcome of CE evolution assuming the energy formalism. Fo… ▽ More To shed light on the origin of magnetic symbiotic stars, we investigated the system FN Sgr in detail. We searched for a reasonable formation pathway to explain its stellar and binary parameters including the magnetic field of the accreting white dwarf. We used the MESA code to carry out pre-CE and post-CE binary evolution and determined the outcome of CE evolution assuming the energy formalism. For the origin and evolution of the white dwarf magnetic field, we adopted the crystallization scenario. We found that FN Sgr can be explained as follows. First, a non-magnetic white dwarf is formed through CE evolution. Later, during post-CE evolution, the white dwarf starts to crystallize and a weak magnetic field is generated. After a few hundred Myr, the magnetic field penetrates the white dwarf surface and becomes detectable. Meanwhile, its companion evolves and becomes an evolved red giant. Subsequently, the white dwarf accretes part of the angular momentum from the red giant stellar winds. As a result, the white dwarf spin period decreases and its magnetic field reaches super-equipartition, getting amplified due to a rotation- and crystallization-driven dynamo. The binary then evolves into a symbiotic star, with a magnetic white dwarf accreting from an evolved red giant through atmospheric Roche-lobe overflow. We conclude that the rotation- and crystallization-driven dynamo scenario, or any age-dependent scenario, can explain the origin of magnetic symbiotic stars reasonably well. This adds another piece to the pile of evidence supporting this scenario. If our formation channel is correct, our findings suggest that white dwarfs in most symbiotic stars formed through CE evolution might be magnetic, provided that the red giant has spent >3 Gyr as a main-sequence star. △ Less

Submitted 19 April, 2024; v1 submitted 13 February, 2024; originally announced February 2024.

Comments: Accepted for publication in A&A

Journal ref: A&A 686, A226 (2024)

The cataclysmic variable orbital period gap: More evident than ever

Authors: Matthias R. Schreiber, Diogo Belloni, Axel D. Schwope

Abstract: Recently, large and homogeneous samples of cataclysmic variables (CVs) identified by the Sloan Digital Sky Survey (SDSS) were published. In these samples, the famous orbital period gap, which is a dearth of systems in the orbital period range ~2-3 hr and the defining feature of most evolutionary models for CVs, has been claimed not to be clearly present. If true, this finding would completely chan… ▽ More Recently, large and homogeneous samples of cataclysmic variables (CVs) identified by the Sloan Digital Sky Survey (SDSS) were published. In these samples, the famous orbital period gap, which is a dearth of systems in the orbital period range ~2-3 hr and the defining feature of most evolutionary models for CVs, has been claimed not to be clearly present. If true, this finding would completely change our picture of CV evolution. In this Letter we focus on potential differences with respect to the orbital period gap between CVs in which the magnetic field of the white dwarf is strong enough to connect with that of the donor star, so-called polars, and non-polar CVs as the white dwarf magnetic field in polars has been predicted to reduce the strength of angular momentum loss through magnetic braking. We separated the SDSS I-IV sample of CVs into polars and non-polar systems and performed statistical tests to evaluate whether the period distributions are bimodal as predicted by the standard model for CV evolution or not. We confirm the existence of a period gap in the SDSS I-IV sample of non-polar CVs with >98 per cent confidence. The boundaries of the orbital period gap are 147 and 191 minutes, with the lower boundary being different to previously published values (129 min). The orbital period distribution of polars from SDSS I-IV is clearly different and does not show a similar period gap. The SDSS samples as well as previous samples of CVs are consistent with the standard theory of CV evolution. Magnetic braking does indeed seem get disrupted around the fully convective boundary, which causes a detached phase during CV evolution. In polars, the white dwarf magnetic field reduces the strength of magnetic braking and consequently the orbital period distribution of polars does not display an equally profound and extended period gap as non-polars. △ Less

Submitted 3 February, 2024; originally announced February 2024.

Comments: Accepted for publication in Astronomy & Astrophysics Letters

Formation of long-period post-common-envelope binaries II. Explaining the self-lensing binary KOI 3278

Authors: Diogo Belloni, Matthias R. Schreiber, Monica Zorotovic

Abstract: The vast majority of close binaries containing a compact object form through common-envelope (CE) evolution. Despite this importance, we struggle to even understand the energy budget of CE evolution. For decades, observed long-period post-CE binaries have been interpreted as evidence for additional energies to contribute during CE evolution. We have recently shown that this argument is based on si… ▽ More The vast majority of close binaries containing a compact object form through common-envelope (CE) evolution. Despite this importance, we struggle to even understand the energy budget of CE evolution. For decades, observed long-period post-CE binaries have been interpreted as evidence for additional energies to contribute during CE evolution. We have recently shown that this argument is based on simplified assumptions for all long-period post-CE binaries containing massive white dwarfs. The only remaining post-CE binary star that has been claimed to require contributions from additional energy sources to understand its formation is KOI 3278. Here we address in detail the potential evolutionary history of KOI 3278. In particular, we investigated whether extra energy sources, such as recombination energy, are indeed required to explain its existence. We used the 1D stellar evolution code MESA to carry out binary evolution simulations and searched for potential formation pathways for KOI 3278 that are able to explain its observed properties. We found that KOI 3278 can be explained if the white dwarf progenitor filled its Roche lobe during a helium shell flash. In this case, the orbital period of KOI 3278 can be reproduced if the CE binding energy is calculated taking into account gravitational energy and thermodynamic internal energy. While the CE evolution that led to the formation of KOI 3278 must have been efficient, that is, most of the available orbital energy must have been used to unbind the CE, recombination energy is not required. We conclude that currently not a single observed post-CE binary requires to assume energy sources other than gravitational and thermodynamic energy to contribute to CE evolution. KOI 3278, however, remains an intriguing post-CE binary as, unlike its siblings, understanding its existence requires highly efficient CE ejection. △ Less

Submitted 8 April, 2024; v1 submitted 30 January, 2024; originally announced January 2024.

Comments: Accepted for publication in A&A

Journal ref: A&A 687, A12 (2024)

Formation of long-period post-common-envelope binaries I. No extra energy is needed to explain oxygen-neon white dwarfs paired with AFGK-type main-sequence stars

Authors: Diogo Belloni, Monica Zorotovic, Matthias R. Schreiber, Steven G. Parsons, Maxwell Moe, James A. Garbutt

Abstract: In this first in a series of papers related to long-period post-common-envelope (CE) binaries, we investigated whether extra energy is required or not to explain the currently known post-CE binaries with sufficiently long orbital periods consisting of oxygen-neon white dwarfs with AFGK-type main-sequence star companions. We carried out binary population simulations with the BSE code and searched f… ▽ More In this first in a series of papers related to long-period post-common-envelope (CE) binaries, we investigated whether extra energy is required or not to explain the currently known post-CE binaries with sufficiently long orbital periods consisting of oxygen-neon white dwarfs with AFGK-type main-sequence star companions. We carried out binary population simulations with the BSE code and searched for their formation pathways. Unlike what has been claimed for a long time, we show that all such post-CE binaries can be explained by assuming inefficient CE evolution, which is consistent with results achieved for the remaining post-CE binaries. There is therefore no need for an extra energy source. We also found that for CE efficiency close to 100%, post-CE binaries hosting oxygen-neon white dwarfs with orbital periods as long as a thousand days can be explained. For all known systems we found formation pathways consisting of CE evolution triggered when a highly evolved (i.e. the envelope mass being comparable to the core mass) thermally-pulsing asymptotic giant branch star fills its Roche lobe at an orbital period of several thousand days. Due to the sufficiently low envelope mass and sufficiently long orbital period, the resulting post-CE orbital period can easily be several tens of days. We conclude that the known post-CE binaries with oxygen-neon white dwarfs and AFGK-type main-sequence stars can be explained without invoking any energy source other than orbital and thermal energy. Our results strengthen the idea that the most common formation pathway of the overall population of post-CE binaries hosting white dwarfs is through inefficient CE evolution. △ Less

Submitted 18 March, 2024; v1 submitted 15 January, 2024; originally announced January 2024.

Comments: Accepted for publication in A&A

Journal ref: A&A 686, A61 (2024)

A Chandra Survey of Milky Way Globular Clusters -- IV. Periodic X-ray sources

Authors: Tong Bao, Zhiyuan Li, Zhongqun Cheng, Diogo Belloni

Abstract: We present a systematic search for periodic X-ray sources in 10 Galactic globular clusters (GCs) utilizing deep archival Chandra observations. By applying the Gregory-Loredo algorithm, we detect 28 periodic signals among 27 independent X-ray sources in 6 GCs, which include 21 newly discovered ones in the X-ray band. The remaining 4 GCs exhibit no periodic X-ray sources, mainly due to a relatively… ▽ More We present a systematic search for periodic X-ray sources in 10 Galactic globular clusters (GCs) utilizing deep archival Chandra observations. By applying the Gregory-Loredo algorithm, we detect 28 periodic signals among 27 independent X-ray sources in 6 GCs, which include 21 newly discovered ones in the X-ray band. The remaining 4 GCs exhibit no periodic X-ray sources, mainly due to a relatively lower sensitivity of the data. Through analysis of their X-ray timing and spectral properties, complemented with available optical and ultraviolet information, we identify 21 of these periodic sources as cataclysmic variables (CVs). Combining with 11 periodic CVs in 47 Tuc similarly identified in the X-ray band, we compile the most comprehensive sample to date of GC CVs with a probable orbital period. The scarcity of old, short-period CVs in GCs compared to the Galactic inner bulge and solar neighborhood, can be attributed to both a selection effect favouring younger, dynamically-formed systems and the hindrance of CV formation through primordial binary evolution by stellar dynamical interactions common to the GC environment. Additionally, we identify a significant fraction of the GC CVs, most with an orbital period below or within the CV period gap, as probable magnetic CVs, but in the meantime there is a deficiency of luminous intermediate polars in the GC sample compared to the solar neighborhood. △ Less

Submitted 24 November, 2023; originally announced November 2023.

Comments: 20 pages, 6 figures; accepted for publication in MNRAS

Evidence for saturated and disrupted magnetic braking from samples of detached close binaries with M and K dwarfs

Authors: Diogo Belloni, Matthias R. Schreiber, Maxwell Moe, Kareem El-Badry, Ken J. Shen

Abstract: Context. Recent observations of close detached eclipsing M and K dwarf binaries have provided substantial support for magnetic saturation when stars rotate sufficiently fast, leading to a magnetic braking (MB) torque proportional to the spin of the star. Aims. We investigated here how strong MB torques need to be to reproduce the observationally-inferred relative numbers of white dwarf plus M dw… ▽ More Context. Recent observations of close detached eclipsing M and K dwarf binaries have provided substantial support for magnetic saturation when stars rotate sufficiently fast, leading to a magnetic braking (MB) torque proportional to the spin of the star. Aims. We investigated here how strong MB torques need to be to reproduce the observationally-inferred relative numbers of white dwarf plus M dwarf post-common-envelope binaries under the assumption of magnetic saturation. Methods. We carried out binary population simulations with the BSE code adopting empirically-derived inter-correlated main-sequence binary distributions as initial binary populations and compared the simulation outcomes with observations. Results. We found that the dearth of extreme mass ratio binaries in the inter-correlated initial distributions is key to reproduce the large fraction of post-common-envelope binaries hosting low-mass M dwarfs (${\sim0.1-0.2}$ M$_\odot$). In addition, orbital angular momentum loss rates due to MB should be high for M dwarfs with radiative cores and orders of magnitude smaller for fully convective stars to explain the observed dramatic change of the fraction of short-period binaries at the fully convective boundary. Conclusions. We conclude that saturated but disrupted, that is, dropping drastically at the fully convective boundary, MB can explain the observations of both close main-sequence binaries containing M and K dwarfs and post-common-envelope binaries. Whether a similar prescription can explain the spin down rates of single stars and of binaries containing more massive stars needs to be tested. △ Less

Submitted 7 November, 2023; originally announced November 2023.

Comments: Accepted for publication in A&A

Period bouncers as detached magnetic cataclysmic variables

Authors: Matthias R. Schreiber, Diogo Belloni, Jan van Roestel

Abstract: The general prediction that more than half of all CVs have evolved past the period minimum is in strong disagreement with observational surveys, which show that the relative number of these objects is just a few per cent. Here, we investigate whether a large number of post-period minimum CVs could detach because of the appearance of a strong white dwarf magnetic field potentially generated by a ro… ▽ More The general prediction that more than half of all CVs have evolved past the period minimum is in strong disagreement with observational surveys, which show that the relative number of these objects is just a few per cent. Here, we investigate whether a large number of post-period minimum CVs could detach because of the appearance of a strong white dwarf magnetic field potentially generated by a rotation- and crystallization-driven dynamo. We used the MESA code to calculate evolutionary tracks of CVs incorporating the spin evolution and cooling as well as compressional heating of the white dwarf. If the conditions for the dynamo were met, we assumed that the emerging magnetic field of the white dwarf connects to that of the companion star and incorporated the corresponding synchronization torque, which transfers spin angular momentum to the orbit. We find that for CVs with donor masses exceeding 0.04 Msun, magnetic fields are generated mostly if the white dwarfs start to crystallize before the onset of mass transfer. It is possible that a few white dwarf magnetic fields are generated in the period gap. For the remaining CVs, the conditions for the dynamo to work are met beyond the period minimum, when the accretion rate decreased significantly. Synchronization torques cause these systems to detach for several Gyrs even if the magnetic field strength of the white dwarf is just one MG. If the rotation- and crystallization-driven dynamo - which is currently the only mechanism that can explain several observational facts related to magnetism in CVs and their progenitors - or a similar temperature-dependent mechanism is responsible for the generation of magnetic field in white dwarfs, most CVs that have evolved beyond the period minimum must detach for several Gyrs at some point. This reduces the predicted number of semi-detached period bouncers by up to 60-80 per cent. △ Less

Submitted 26 October, 2023; originally announced October 2023.

Comments: A&A letters, in press

Reversing the verdict: Cataclysmic variables could be the dominant progenitors of AM CVn binaries after all

Authors: Diogo Belloni, Matthias R. Schreiber

Abstract: Context. AM CVn binaries are potential progenitors of thermonuclear supernovae and strong sources of persistent gravitational wave radiation. For a long time, it has been believed that these systems cannot descend from cataclysmic variables (CVs), at least not in large numbers, because the initial conditions need to be fine-tuned and, even worse, the resulting surface hydrogen abundance would be h… ▽ More Context. AM CVn binaries are potential progenitors of thermonuclear supernovae and strong sources of persistent gravitational wave radiation. For a long time, it has been believed that these systems cannot descend from cataclysmic variables (CVs), at least not in large numbers, because the initial conditions need to be fine-tuned and, even worse, the resulting surface hydrogen abundance would be high enough to be detected which contradicts a defining feature of AM CVn binaries. Aims. Here we show that both claimed weaknesses of the CV formation channel for AM CVn binaries are model-dependent and rely on poorly constrained assumptions for magnetic braking. Methods. We performed binary evolution simulations with the MESA code for different combinations of post-common-envelope white dwarf and companion masses as well as orbital periods assuming the CARB model for strong magnetic braking. Results. We found that AM CVn binaries with extremely-low surface hydrogen abundances are one natural outcome of CV evolution if the donor star has developed a non-negligible helium core prior to the onset of mass transfer. In this case, after hydrogen envelope exhaustion during CV evolution, the donor becomes degenerate and its surface hydrogen abundance substantially drops and becomes undetectable. Our simulations also show that the CV formation channel is able to explain the observed AM CVn binaries with very low mass and bloated donor stars (Gaia14aae and ZTF J1637+49). Conclusions. CVs with evolved donors are likely the progenitors of at least a fraction of AM CVn binaries. △ Less

Submitted 21 August, 2023; v1 submitted 30 May, 2023; originally announced May 2023.

Comments: Accepted for publication in A&A

Journal ref: A&A 678, A34 (2023)

Formation and Evolution of Accreting Compact Objects

Authors: Diogo Belloni, Matthias R. Schreiber

Abstract: Accreting compact objects are crucial to understand several important astrophysical phenomena such as Type Ia supernovae, gravitational waves, or X-ray and $γがんま$-ray bursts. In addition, they are natural laboratories to infer fundamental properties of stars, to investigate high-energy phenomena and accretion processes, to test theories of stellar and binary evolution, to explore interactions between… ▽ More Accreting compact objects are crucial to understand several important astrophysical phenomena such as Type Ia supernovae, gravitational waves, or X-ray and $γがんま$-ray bursts. In addition, they are natural laboratories to infer fundamental properties of stars, to investigate high-energy phenomena and accretion processes, to test theories of stellar and binary evolution, to explore interactions between high-density plasma and very strong magnetic fields, to examine the interplay between binary evolution and dynamical interactions (in the case they belong to dense star clusters), and they can even be used as a probe for the assembling process of galaxies over cosmic time-scales. Despite the fundamental importance of accreting compact objects for astrophysics and recent progress with the comprehension of these fascinating objects, we still do not fully understand how they form and evolve. In this chapter, we will review the current theoretical status of our knowledge on these objects, and will discuss standing problems and potential solutions to them. △ Less

Submitted 22 June, 2023; v1 submitted 15 March, 2023; originally announced March 2023.

Comments: Invited chapter for the Handbook of X-ray and Gamma-ray Astrophysics (Editors: Cosimo Bambi, Andrea Santangelo; Publisher: Springer Singapore)

A candidate magnetic helium core white dwarf in the globular cluster NGC 6397

Authors: Manuel Pichardo Marcano, Liliana E. Rivera Sandoval, Thomas J. Maccarone, Rene D. Rohrmann, Craig O. Heinke, Diogo Belloni, Leandro G. Althaus, Arash Bahramian

Abstract: We report a peculiar variable blue star in the globular cluster NGC 6397, using Hubble Space Telescope optical imaging. Its position in the colour-magnitude diagrams, and its spectrum, are consistent with this star being a helium core white dwarf (He WD) in a binary system. The optical light curve shows a periodicity at 18.5 hours. We argue that this periodicity is due to the rotation of the WD an… ▽ More We report a peculiar variable blue star in the globular cluster NGC 6397, using Hubble Space Telescope optical imaging. Its position in the colour-magnitude diagrams, and its spectrum, are consistent with this star being a helium core white dwarf (He WD) in a binary system. The optical light curve shows a periodicity at 18.5 hours. We argue that this periodicity is due to the rotation of the WD and possibly due to magnetic spots on the surface of the WD. This would make this object the first candidate magnetic He WD in any globular cluster (GC), and the first candidate magnetic WD in a detached binary system in any GC and one of the few He WDs with a known rotation period and of magnetic nature. Another possibility is that this system is a He WD in a binary system with another WD or another degenerate object, which would make this object one of the few candidate non-accreting double degenerate binaries in any GC. △ Less

Submitted 7 March, 2023; originally announced March 2023.

Comments: Accepted for publication in MNRAS

MOCCA: Dynamics and evolution of binary stars of multiple stellar populations in tidally filling and underfilling globular star clusters

Authors: Arkadiusz Hypki, Mirek Giersz, Jongsuk Hong, Agostino Leveque, Abbas Askar, Diogo Belloni, Magdalena Otulakowska-Hypka

Abstract: We present an upgraded version of the \MOCCA code for the study of dynamical evolution of globular clusters (GCs) and its first application to the study of evolution of multiple stellar populations. We explore initial conditions spanning different structural parameters for the first (FG) and second generation of stars (SG) and we analyze their effect on the binary dynamics and survival. Here, we f… ▽ More We present an upgraded version of the \MOCCA code for the study of dynamical evolution of globular clusters (GCs) and its first application to the study of evolution of multiple stellar populations. We explore initial conditions spanning different structural parameters for the first (FG) and second generation of stars (SG) and we analyze their effect on the binary dynamics and survival. Here, we focus on the number ratio of FG and SG binaries, its spatial variation, and the way their abundances are affected by various cluster initial properties. We find that present-day SG stars are more abundant in clusters that were initially tidally filling. Conversely, FG stars stay more abundant in clusters that were initially tidally underfilling. We find that the ratio between binary fractions is not affected by the way we calculate these fractions (e.g. only main-sequence binaries (MS) or observational binaries, i.e. MS stars $> 0.4 M_{\odot}$ mass ratios $> 0.5$). This implies that the MS stars themselves are a very good proxy for probing entire populations of FG and SG. We also discuss how it relates to the observations of Milky Way GCs. We show that \MOCCA models are able to reproduce the observed range of SG fractions for Milky Way GCs for which we know these fractions. We show how the SG fractions depend on the initial conditions and provide some constraints for the initial conditions to have more numerous FG or SG stars at the Hubble time. △ Less

Submitted 3 October, 2022; v1 submitted 11 May, 2022; originally announced May 2022.

Comments: Acceped to MNRAS

Magnetic dynamos in white dwarfs -III: explaining the occurrence of strong magnetic fields in close double white dwarfs

Authors: Matthias. R. Schreiber, Diogo Belloni, Monica Zorotovic, Sarai Zapata, Boris T. Gänsicke, Steven G. Parsons

Abstract: The origin of strong magnetic fields in white dwarfs has been a puzzle for decades. Recently, a dynamo mechanism operating in rapidly rotating and crystallizing white dwarfs has been suggested to explain the occurrence rates of strong magnetic fields in white dwarfs with close low-mass main sequence star companions. Here we investigate whether the same mechanism may produce strong magnetic fields… ▽ More The origin of strong magnetic fields in white dwarfs has been a puzzle for decades. Recently, a dynamo mechanism operating in rapidly rotating and crystallizing white dwarfs has been suggested to explain the occurrence rates of strong magnetic fields in white dwarfs with close low-mass main sequence star companions. Here we investigate whether the same mechanism may produce strong magnetic fields in close double white dwarfs. The only known strongly magnetic white dwarf that is part of a close double white dwarf system, the magnetic component of NLTT 12758, is rapidly rotating and likely crystallizing and therefore the proposed dynamo mechanism represents an excellent scenario for the origin of its magnetic field. Presenting a revised formation scenario for NLTT 12758, we find a natural explanation for the rapid rotation of the magnetic component. We furthermore show that it is not surprising that strong magnetic fields have not been detected in all other known double white dwarfs. We therefore conclude that the incidence of magnetic fields in close double white dwarfs supports the idea that a rotation and crystallization driven dynamo plays a major role in the generation of strong magnetic fields in white dwarfs. △ Less

Submitted 16 April, 2022; originally announced April 2022.

Comments: accepted for publication in MNRAS

The White Dwarf Binary Pathways Survey VII: Evidence for a bi-modal distribution of post mass transfer systems?

Authors: F. Lagos, M. R. Schreiber, S. G. Parsons, O. Toloza, B. T. Gänsicke, M. S. Hernandez, L. Schmidtobreick, D. Belloni

Abstract: Binary systems consisting of a white dwarf and a main-sequence companion with orbital periods up to $\approx 100$ d are often thought to be formed through common envelope evolution which is still poorly understood. To provide new observational constraints on the physical processes involved in the formation of these objects, we are conducting a large-scale survey of close binaries consisting of a w… ▽ More Binary systems consisting of a white dwarf and a main-sequence companion with orbital periods up to $\approx 100$ d are often thought to be formed through common envelope evolution which is still poorly understood. To provide new observational constraints on the physical processes involved in the formation of these objects, we are conducting a large-scale survey of close binaries consisting of a white dwarf and an A to K-type companion. Here we present three systems with eccentric orbits and orbital periods between $\approx10-42$ d discovered by our survey. Based on HST spectroscopy and high angular resolution images obtained with SPHERE-IRDIS, we find that two of these systems are most likely triple systems while the remaining one could be either a binary or a hierarchical triple but none of them is a post common envelope binary (PCEB). The discovery of these systems shows that our survey is capable to detect systems with orbital periods of the order of weeks, but all six PCEBs we have previously discovered have periods below 2.5 d. We suggest that the fact that all of the systems we identify with periods of the order of weeks are not PCEBs indicates a transition between two different mechanisms responsible for the formation of very close ($\lesssim 10$ d) and somewhat wider WD+AFGK binaries: common envelope evolution and non-conservative stable mass transfer. △ Less

Submitted 16 March, 2022; originally announced March 2022.

Comments: 12 pages, 6 figures, 4 tables. Accepted for publication in MNRAS

MSC Class: 85-02

Constraining the Evolution of Cataclysmic Variables via the Masses and Accretion Rates of their Underlying White Dwarfs

Authors: A. F. Pala, B. T. Gänsicke, D. Belloni, S. G. Parsons, T. R. Marsh, M. R. Schreiber, E. Breedt, C. Knigge, E. M. Sion, P. Szkody, D. Townsley, L. Bildsten, D. Boyd, M. J. Cook, D. De Martino, P. Godon, S. Kafka, V. Kouprianov, K. S. Long, B. Monard, G. Myers, P. Nelson, D. Nogami, A. Oksanen, R. Pickard , et al. (6 additional authors not shown)

Abstract: We report on the masses ($M_\mathrm{WD}$), effective temperatures ($T_\mathrm{eff}$) and secular mean accretion rates ($\langle \dot{M} \rangle$) of 43 cataclysmic variable (CV) white dwarfs, 42 of which were obtained from the combined analysis of their $\mathit{Hubble~Space~Telescope}$ ultraviolet data with the parallaxes provided by the Early Third Data Release of the $\mathit{Gaia}$ space missi… ▽ More We report on the masses ($M_\mathrm{WD}$), effective temperatures ($T_\mathrm{eff}$) and secular mean accretion rates ($\langle \dot{M} \rangle$) of 43 cataclysmic variable (CV) white dwarfs, 42 of which were obtained from the combined analysis of their $\mathit{Hubble~Space~Telescope}$ ultraviolet data with the parallaxes provided by the Early Third Data Release of the $\mathit{Gaia}$ space mission, and one from the white dwarf gravitational redshift. Our results double the number of CV white dwarfs with an accurate mass measurement, bringing the total census to 89 systems. From the study of the mass distribution, we derive $\langle M_\mathrm{WD} \rangle = 0.81^{+0.16}_{-0.20}\,\mathrm{M_\odot}$, in perfect agreement with previous results, and find no evidence of any evolution of the mass with orbital period. Moreover, we identify five systems with $M_\mathrm{WD} < 0.5\mathrm{M_\odot}$, which are most likely representative of helium-core white dwarfs, showing that these CVs are present in the overall population. We reveal the presence of an anti-correlation between the average accretion rates and the white dwarf masses for the systems below the $2-3\,$h period gap. Since $\langle \dot{M} \rangle$ reflects the rate of system angular momentum loss, this correlation suggests the presence of an additional mechanism of angular momentum loss that is more efficient at low white dwarf masses. This is the fundamental concept of the recently proposed empirical prescription of consequential angular momentum loss (eCAML) and our results provide observational support for it, although we also highlight how its current recipe needs to be refined to better reproduce the observed scatter in $T_\mathrm{eff}$ and $\langle \dot{M} \rangle$, and the presence of helium-core white dwarfs. △ Less

Submitted 26 November, 2021; originally announced November 2021.

Comments: Accepted for publication in MNRAS. 22 pages of main body (6 tables and 14 figures) and 54 pages of appendices. Appendix B includes the best-fitting parameters and models to the HST data

Breaking the degeneracy in magnetic cataclysmic variable X-ray spectral modeling using X-ray light curves

Authors: Diogo Belloni, Claudia V. Rodrigues, Matthias R. Schreiber, Manuel Castro, Joaquim E. R. Costa, Takayuki Hayashi, Isabel J. Lima, Gerardo J. M. Luna, Murilo Martins, Alexandre S. Oliveira, Steven G. Parsons, Karleyne M. G. Silva, Paulo E. Stecchini, Teresa J. Stuchi, Monica Zorotovic

Abstract: We present an analysis of mock X-ray spectra and light curves of magnetic cataclysmic variables using an upgraded version of the 3D CYCLOPS code. This 3D representation of the accretion flow allows us to properly model total and partial occultation of the post-shock region by the white dwarf as well as the modulation of the X-ray light curves due to the phase-dependent extinction of the pre-shock… ▽ More We present an analysis of mock X-ray spectra and light curves of magnetic cataclysmic variables using an upgraded version of the 3D CYCLOPS code. This 3D representation of the accretion flow allows us to properly model total and partial occultation of the post-shock region by the white dwarf as well as the modulation of the X-ray light curves due to the phase-dependent extinction of the pre-shock region. We carried out detailed post-shock region modeling in a four-dimensional parameter space by varying the white dwarf mass and magnetic field strength as well as the magnetosphere radius and the specific accretion rate. To calculate the post-shock region temperature and density profiles, we assumed equipartition between ions and electrons, took into account the white dwarf gravitational potential, the finite size of the magnetosphere and a dipole-like magnetic field geometry, and considered cooling by both bremsstrahlung and cyclotron radiative processes. By investigating the impact of the parameters on the resulting X-ray continuum spectra, we show that there is an inevitable degeneracy in the four-dimensional parameter space investigated here, which compromises X-ray continuum spectral fitting strategies and can lead to incorrect parameter estimates. However, the inclusion of X-ray light curves in different energy ranges can break this degeneracy, and it therefore remains, in principle, possible to use X-ray data to derive fundamental parameters of magnetic cataclysmic variables, which represents an essential step toward understanding their formation and evolution. △ Less

Submitted 22 July, 2021; originally announced July 2021.

Comments: 42 pages, accepted for publication in The Astrophysical Journal Supplement Series

Magnetic dynamos in white dwarfs -- II. Relating magnetism and pollution

Authors: Matthias R. Schreiber, Diogo Belloni, Boris T. Gaensicke, Steven G. Parsons

Abstract: We investigate whether the recently suggested rotation and crystallization driven dynamo can explain the apparent increase of magnetism in old metal polluted white dwarfs. We find that the effective temperature distribution of polluted magnetic white dwarfs is in agreement with most/all of them having a crystallizing core and increased rotational velocities are expected due to accretion of planeta… ▽ More We investigate whether the recently suggested rotation and crystallization driven dynamo can explain the apparent increase of magnetism in old metal polluted white dwarfs. We find that the effective temperature distribution of polluted magnetic white dwarfs is in agreement with most/all of them having a crystallizing core and increased rotational velocities are expected due to accretion of planetary material which is evidenced by the metal absorption lines. We conclude that a rotation and crystallization driven dynamo offers not only an explanation for the different occurrence rates of strongly magnetic white dwarfs in close binaries, but also for the high incidence of weaker magnetic fields in old metal polluted white dwarfs. △ Less

Submitted 17 June, 2021; originally announced June 2021.

Comments: Accepted for publication in MNRAS Letters

Magnetic dynamos in white dwarfs -- I. Explaining the dearth of bright intermediate polars in globular clusters

Authors: Diogo Belloni, Matthias R. Schreiber, Maurizio Salaris, Thomas J. Maccarone, Monica Zorotovic

Abstract: Recently, Bahramian et al. investigated a large sample of globular clusters (GCs) and found that bright intermediate polars (IPs) are a factor of 10 less frequent in GCs than in the Galactic field. We theoretically investigate here this discrepancy based on GC numerical simulations. We found that, due to disruptive dynamical interaction, there is on average a reduction of only half of bright IP pr… ▽ More Recently, Bahramian et al. investigated a large sample of globular clusters (GCs) and found that bright intermediate polars (IPs) are a factor of 10 less frequent in GCs than in the Galactic field. We theoretically investigate here this discrepancy based on GC numerical simulations. We found that, due to disruptive dynamical interaction, there is on average a reduction of only half of bright IP progenitors, which is clearly not enough to explain the observed deficiency. However, if the rotation- and crystallization-driven dynamo scenario recently proposed by Schreiber et al. is incorporated in the simulations, the observed rareness of bright IPs in GCs can be reproduced. This is because bright cataclysmic variables in GCs are typically very old systems ($\gtrsim$ 10 Gyr), with white dwarfs that almost fully crystallized before mass transfer started, which does not allow strong magnetic fields to be generated. The observed mass density of bright IPs in GCs can be recovered if around one third of the bright cataclysmic variables dynamically formed through mergers have magnetic field strengths similar to those of IPs. We conclude that the observed paucity of bright IPs in GCs is a natural consequence of the newly proposed rotation- and crystallization-driven dynamo scenario. △ Less

Submitted 24 May, 2021; originally announced May 2021.

Comments: letter accepted for publication by MNRAS

Preparing the next gravitational million-body simulations: Evolution of single and binary stars in Nbody6++GPU, MOCCA and McLuster

Authors: A. W. H. Kamlah, A. Leveque, R. Spurzem, M. Arca Sedda, A. Askar, S. Banerjee, P. Berczik, M. Giersz, J. Hurley, D. Belloni, L. Kühmichel, L. Wang

Abstract: We present the implementation of updated stellar evolution recipes in the codes \texttt{Nbody6++GPU, MOCCA} and \texttt{McLuster}. We test them through numerical simulations of star clusters containing $1.1\times 10^5$ stars (with $2.0\times 10^4$ in primordial hard binaries) performing high-resolution direct $N$-body (\texttt{Nbody6++GPU}) and Monte-Carlo (\texttt{MOCCA}) simulations to an age of… ▽ More We present the implementation of updated stellar evolution recipes in the codes \texttt{Nbody6++GPU, MOCCA} and \texttt{McLuster}. We test them through numerical simulations of star clusters containing $1.1\times 10^5$ stars (with $2.0\times 10^4$ in primordial hard binaries) performing high-resolution direct $N$-body (\texttt{Nbody6++GPU}) and Monte-Carlo (\texttt{MOCCA}) simulations to an age of 10~Gyr. We compare models implementing either delayed or core-collapse supernovae mechanisms, a different mass ratio distribution for binaries, and white dwarf natal kicks enabled/disabled. Compared to \texttt{Nbody6++GPU}, the \texttt{MOCCA} models appear to be denser, with a larger scatter in the remnant masses, and a lower binary fraction on average. The \texttt{MOCCA} models produce more black holes (BHs) and helium white dwarfs (WDs), whilst \texttt{Nbody6++GPU} models are characterised by a much larger amount of WD-WD binaries. The remnant kick velocity and escape speed distributions are similar for the BHs and neutron stars (NSs), and some NSs formed via electron-capture supernovae, accretion-induced collapse or merger-induced collapse escape the cluster in all simulations. The escape speed distributions for the WDs, on the other hand, are very dissimilar. We categorise the stellar evolution recipes available in \texttt{Nbody6++GPU}, \texttt{MOCCA} and \texttt{Mcluster} into four levels: the one implemented in previous \texttt{Nbody6++GPU} and \texttt{MOCCA} versions (\texttt{level A}), state-of-the-art prescriptions (\texttt{level B}), some in a testing phase (\texttt{level C}), and those that will be added in future versions of our codes. △ Less

Submitted 21 December, 2021; v1 submitted 17 May, 2021; originally announced May 2021.

Comments: 31 pages, 15 figures, 5 tables, Accepted by MNRAS

The origin and evolution of magnetic white dwarfs in close binary stars

Authors: Matthias R. Schreiber, Diogo Belloni, Boris T. Gaensicke, Steven G. Parsons, Monica Zorotovic

Abstract: The origin of magnetic fields in white dwarfs remains a fundamental unresolved problem in stellar astrophysics. In particular, the very different fractions of strongly (exceeding 1 MG) magnetic white dwarfs in evolutionarily linked populations of close white dwarf binary stars cannot be reproduced by any scenario suggested so far. Strongly magnetic white dwarfs are absent among detached white dwar… ▽ More The origin of magnetic fields in white dwarfs remains a fundamental unresolved problem in stellar astrophysics. In particular, the very different fractions of strongly (exceeding 1 MG) magnetic white dwarfs in evolutionarily linked populations of close white dwarf binary stars cannot be reproduced by any scenario suggested so far. Strongly magnetic white dwarfs are absent among detached white dwarf binary stars that are younger than approximately 1 Gyr. In contrast, in semi-detached cataclysmic variables in which the white dwarf accretes from a low-mass star companion, more than one third host a strongly magnetic white dwarf. Here we present binary star evolutionary models that include the spin evolution of accreting white dwarfs and crystallization of their cores, as well as magnetic field interactions between both stars. We show that a crystallization- and rotation-driven dynamo similar to those working in planets and low-mass stars can generate strong magnetic fields in the white dwarfs in cataclysmic variables which explains their large fraction among the observed population. When the magnetic field generated in the white dwarfs connects with that of the secondary stars, synchronization torques and reduced angular momentum loss cause the binary to detach for a relatively short period of time. The few known strongly magnetic white dwarfs in detached binaries, including AR Sco, are in this detached phase. △ Less

Submitted 29 April, 2021; originally announced April 2021.

Journal ref: Published by Nature Astronomy (2021)

Search for magnetic accretion in SW Sextantis systems

Authors: I. J. Lima, C. V. Rodrigues, C. E. Ferreira Lopes, P. Szkody, F. J. Jablonski, A. S. Oliveira, K. M. G. Silva, D. Belloni, M. S. Palhares, S. Shugarov, R. Baptista, L. A. Almeida

Abstract: SW Sextantis systems are nova-like cataclysmic variables that have unusual spectroscopic properties, which are thought to be caused by an accretion geometry having part of the mass flux trajectory out of the orbital plane. Accretion onto a magnetic white dwarf is one of the proposed scenarios for these systems. To verify this possibility, we analysed photometric and polarimetric time-series data f… ▽ More SW Sextantis systems are nova-like cataclysmic variables that have unusual spectroscopic properties, which are thought to be caused by an accretion geometry having part of the mass flux trajectory out of the orbital plane. Accretion onto a magnetic white dwarf is one of the proposed scenarios for these systems. To verify this possibility, we analysed photometric and polarimetric time-series data for a sample of six SW Sex stars. We report possible modulated circular polarization in BO Cet, SW Sex, and UU Aqr with periods of 11.1, 41.2 and 25.7 min, respectively, and less significant periodicities for V380 Oph at 22 min and V442 Oph at 19.4 min. We confirm previous results that LS Peg shows variable circular polarization. However, we determine a period of 18.8 min, which is different from the earlier reported value. We interpret these periods as the spin periods of the white dwarfs. Our polarimetric results indicate that 15% of the SW Sex systems have direct evidence of magnetic accretion. We also discuss SW Sex objects within the perspective of being magnetic systems, considering the latest findings about cataclysmic variables demography, formation and evolution. △ Less

Submitted 2 March, 2021; originally announced March 2021.

Comments: 29 pages, 17 figures

Properties of Cataclysmic Variables in Globular Clusters

Authors: Diogo Belloni, Liliana E. Rivera Sandoval

Abstract: The study of star clusters plays an important role in our understanding of the Universe since these systems are natural laboratories for testing theories of stellar dynamics and evolution. Particularly, globular clusters (GCs) are one of the most important objects for studying the formation and the physical nature of exotic systems which in turn provide basic information and tools that can help us… ▽ More The study of star clusters plays an important role in our understanding of the Universe since these systems are natural laboratories for testing theories of stellar dynamics and evolution. Particularly, globular clusters (GCs) are one of the most important objects for studying the formation and the physical nature of exotic systems which in turn provide basic information and tools that can help us to understand the formation and evolution processes of star clusters themselves, galaxies and, in general, the young Universe. Among the most interesting objects in GCs are the cataclysmic variables (CVs), which are interacting binaries harboring a white dwarf accreting from a low-mass companion. Since GC densities are sufficiently high that dynamical encounters involving binaries should be common, CV progenitors are expected to be affected by dynamics in some way in the early stages. In this article we review the formation channels and the influence of dynamics on the CV population in GCs. In particular, we review recent progress in numerical simulations. Furthermore, we discuss observational properties of CVs in GCs and the techniques used to identify and study them. We focus the discussion on the multi-wavelength observations carried out with HST and Chandra on the best-studied GCs NGC 6397, NGC 6752, 47 Tucanae and Omega Centauri; on the recent spectroscopic findings with MUSE, and on updates regarding the correlation between the number of faint X-ray sources and the cluster stellar encounter rate. Finally, we discuss some observational prospects that might potentially help future investigations. △ Less

Submitted 28 August, 2020; originally announced August 2020.

Comments: Invited review at the Palermo workshop on "The Golden Age of Cataclysmic Variables and Related Objects"; Refereed manuscript accepted for publication by PoS-SISSA

On the absence of symbiotic stars in globular clusters

Authors: Diogo Belloni, Joanna Mikołajewska, Krystian Iłkiewicz, Matthias R. Schreiber, Mirek Giersz, Liliana E. Rivera Sandoval, Claudia V. Rodrigues

Abstract: Even though plenty of symbiotic stars (SySts) have been found in the Galactic field and nearby galaxies, not a single one has ever been confirmed in a Galactic globular cluster (GC). We investigate the lack of such systems in GCs for the first time by analysing 144 GC models evolved with the MOCCA code, which have different initial properties and are roughly representative of the Galactic GC popul… ▽ More Even though plenty of symbiotic stars (SySts) have been found in the Galactic field and nearby galaxies, not a single one has ever been confirmed in a Galactic globular cluster (GC). We investigate the lack of such systems in GCs for the first time by analysing 144 GC models evolved with the MOCCA code, which have different initial properties and are roughly representative of the Galactic GC population. We focus here on SySts formed through the wind-accretion channel, which can be consistently modelled in binary population synthesis codes. We found that the orbital periods of the majority of such SySts are sufficiently long (${\gtrsim10^3}$ d) so that, for very dense GC models, dynamical interactions play an important role in destroying their progenitors before the present day (${\sim11-12}$ Gyr). In less dense GC models, some SySts are still predicted to exist. However, these systems tend to be located far from the central parts (${\gtrsim70}$ per cent are far beyond the half-light radius) and are sufficiently rare (${\lesssim1}$ per GC per Myr), which makes their identification rather difficult in observational campaigns. We propose that future searches for SySts in GCs should be performed in the outskirts of nearby low-density GCs with sufficiently long half-mass relaxation times and relatively large Galactocentric distances. Finally, we obtained spectra of the candidate proposed in $ωおめが$ Cen (SOPS IV e-94) and showed that this object is most likely not a SySt. △ Less

Submitted 30 June, 2020; v1 submitted 11 April, 2020; originally announced April 2020.

Comments: accepted by MNRAS, 12 pages

Are white dwarf magnetic fields in close binaries generated during common-envelope evolution?

Authors: Diogo Belloni, Matthias R. Schreiber

Abstract: Understanding the origin of the magnetic fields in white dwarfs (WDs) has been a puzzle for decades. A scenario that has gained considerable attention in the past years assumes that such magnetic fields are generated through a dynamo process during common-envelope evolution. We performed binary population models using an up-to-date version of the BSE code to confront the predictions of this model… ▽ More Understanding the origin of the magnetic fields in white dwarfs (WDs) has been a puzzle for decades. A scenario that has gained considerable attention in the past years assumes that such magnetic fields are generated through a dynamo process during common-envelope evolution. We performed binary population models using an up-to-date version of the BSE code to confront the predictions of this model with observational results. We found that this hypothesis can only explain the observed distribution of WD magnetic fields in polars and pre-polars and the low-temperature WDs in pre-polars if it is re-scaled to fit the observational data. Furthermore, in its present version, the model fails to explain the absence of young close detached WD+M-dwarf binaries harbouring hot magnetic WDs and predicts that the overwhelming majority of WDs in close binaries should be strongly magnetic, which is also in serious conflict with the observations. We conclude that either the common-envelope dynamo scenario needs to be substantially revised or a different mechanism is responsible for the generation of strong WD magnetic fields in close binaries. △ Less

Submitted 23 December, 2019; v1 submitted 18 October, 2019; originally announced October 2019.

Comments: 7 pages; accepted for publication in MNRAS

Evidence for reduced magnetic braking in polars from binary population models

Authors: Diogo Belloni, Matthias R. Schreiber, Anna F. Pala, Boris T. Gänsicke, Mónica Zorotovic, Claudia V. Rodrigues

Abstract: We present the first population synthesis of synchronous magnetic cataclysmic variables, called polars, taking into account the effect of the white dwarf (WD) magnetic field on angular momentum loss. We implemented the reduced magnetic braking (MB) model proposed by Li, Wu & Wickramasinghe into the Binary Stellar Evolution (BSE) code recently calibrated for cataclysmic variable (CV) evolution. We… ▽ More We present the first population synthesis of synchronous magnetic cataclysmic variables, called polars, taking into account the effect of the white dwarf (WD) magnetic field on angular momentum loss. We implemented the reduced magnetic braking (MB) model proposed by Li, Wu & Wickramasinghe into the Binary Stellar Evolution (BSE) code recently calibrated for cataclysmic variable (CV) evolution. We then compared separately our predictions for polars and non-magnetic CVs with a large and homogeneous sample of observed CVs from the Sloan Digital Sky Survey. We found that the predicted orbital period distributions and space densities agree with the observations if period bouncers are excluded. For polars, we also find agreement between predicted and observed mass transfer rates, while the mass transfer rates of non-magnetic CVs with periods ${\gtrsim3}$ hr drastically disagree with those derived from observations. Our results provide strong evidence that the reduced MB model for the evolution of highly magnetized accreting WDs can explain the observed properties of polars. The remaining main issues in our understanding of CV evolution are the origin of the large number of highly magnetic WDs, the large scatter of the observed mass transfer rates for non-magnetic systems with periods ${\gtrsim3}$ hr, and the absence of period bouncers in observed samples. △ Less

Submitted 23 December, 2019; v1 submitted 14 October, 2019; originally announced October 2019.

Comments: 15 pages, accepted for publication in MNRAS

Optical observations and cyclops post-shock region modelling of the polar V348 Pav

Authors: Alexandre S. Oliveira, Claudia V. Rodrigues, Matheus S. Palhares, Marcos P. Diaz, Diogo Belloni, Karleyne M. G. Silva

Abstract: Post-shock regions (PSR) of polar cataclysmic variables produce most of their luminosity and give rise to high circular polarization in optical wavelengths and strong variability on the white dwarf rotation period, which are distinctive features of these systems. To investigate the polar candidate V348 Pav, we obtained a comprehensive observational set including photometric, polarimetric, and spec… ▽ More Post-shock regions (PSR) of polar cataclysmic variables produce most of their luminosity and give rise to high circular polarization in optical wavelengths and strong variability on the white dwarf rotation period, which are distinctive features of these systems. To investigate the polar candidate V348 Pav, we obtained a comprehensive observational set including photometric, polarimetric, and spectroscopic data, which was used to constrain the post-shock properties of the system. The object presents high circular polarization ($\sim30$ per cent) and high He II 4686A to H$~βべーた$ line ratio, confirming it is a polar. From both radial velocities and light curves, we determined an orbital period of 79.98~min, close to the orbital period minimum of CVs. The H$~βべーた$ radial velocity curve has a semi-amplitude of $141.4\pm1.5$ km s$^{-1}$. Doppler tomography showed that most of the spectral line emission in this system is originated in the region of the companion star facing the WD, possibly irradiated by the emission related to the PSR. We modelled the PSR using the Cyclops code. The PSR density and temperature profiles, obtained by a proper solution of the hydro-thermodynamic equations, were used in a 3D radiative transfer solution that takes into account the system geometry. We could reproduce the V348 Pav \textit{B, V, R}, and \textit{I} photometric and polarimetric data using a model with a WD magnetic field of $\sim28$ MG, a WD mass of $\sim0.85$~M$_{\odot}$ and a low ($\sim25^{\circ}$) orbital inclination. These values for the WD mass and orbital inclination are consistent with the measured radial velocities. △ Less

Submitted 2 September, 2019; originally announced September 2019.

Are most Cataclysmic Variables in Globular Clusters dynamically formed?

Authors: Diogo Belloni, Mirek Giersz, Liliana E. Rivera Sandoval, Abbas Askar, Pawel Ciecielag

Abstract: We have been investigating populations of cataclysmic variables (CVs) in a set of more than 300 globular cluster (GC) models evolved with the MOCCA code. One of the main questions we have intended to answer is whether most CVs in GCs are dynamically formed or not. Contrary to what has been argued for a long time, we found that dynamical destruction of primordial CV progenitors is much stronger in… ▽ More We have been investigating populations of cataclysmic variables (CVs) in a set of more than 300 globular cluster (GC) models evolved with the MOCCA code. One of the main questions we have intended to answer is whether most CVs in GCs are dynamically formed or not. Contrary to what has been argued for a long time, we found that dynamical destruction of primordial CV progenitors is much stronger in GCs than dynamical formation of CVs. In particular, we found that, on average, the detectable CV population is predominantly composed of CVs formed via a typical common envelope phase (> 70 per cent). However, core-collapsed models tend to have higher fractions of bright CVs than non-core-collapsed ones, which suggests then that the formation of CVs is indeed slightly favoured through strong dynamical interactions in core-collapsed GCs, due to the high stellar densities in their cores. △ Less

Submitted 15 July, 2019; originally announced July 2019.

Comments: 4 pages, 1 figures, contribution to the proceedings of the IAUえーゆー Symposium No. 351, 2019, "Star Clusters: from the Milky Way to the Early Universe"

MOCCA-SURVEY database I. Accreting white dwarf binary systems in globular clusters -- IV. cataclysmic variables -- properties of bright and faint populations

Authors: Diogo Belloni, Mirek Giersz, Liliana E. Rivera Sandoval, Abbas Askar, Paweł Ciecieląg

Abstract: We investigate here populations of cataclysmic variables (CVs) in a set of 288 globular cluster (GC) models evolved with the MOCCA code. This is by far the largest sample of GC models ever analysed with respect to CVs. Contrary to what has been argued for a long time, we found that dynamical destruction of primordial CV progenitors is much stronger in GCs than dynamical formation of CVs, and that… ▽ More We investigate here populations of cataclysmic variables (CVs) in a set of 288 globular cluster (GC) models evolved with the MOCCA code. This is by far the largest sample of GC models ever analysed with respect to CVs. Contrary to what has been argued for a long time, we found that dynamical destruction of primordial CV progenitors is much stronger in GCs than dynamical formation of CVs, and that dynamically formed CVs and CVs formed under no/weak influence of dynamics have similar white dwarf mass distributions. In addition, we found that, on average, the detectable CV population is predominantly composed of CVs formed via typical common envelope phase (CEP) ($\gtrsim70$ per cent), that only $\approx2-4$ per cent of all CVs in a GC is likely to be detectable, and that core-collapsed models tend to have higher fractions of bright CVs than non-core-collapsed ones. We also consistently show, for the first time, that the properties of bright and faint CVs can be understood by means of the pre-CV and CV formation rates, their properties at their formation times and cluster half-mass relaxation times. Finally, we show that models following the initial binary population proposed by Kroupa and set with low CEP efficiency better reproduce the observed amount of CVs and CV candidates in NGC 6397, NGC 6752 and 47 Tuc. To progress with comparisons, the essential next step is to properly characterize the candidates as CVs (e.g. by obtaining orbital periods and mass ratios). △ Less

Submitted 12 November, 2018; originally announced November 2018.

Comments: 18 pages, 13 figures; accepted for publication in MNRAS

No cataclysmic variables missing: higher merger rate brings into agreement observed and predicted space densities

Authors: Diogo Belloni, Matthias R. Schreiber, Mónica Zorotovic, Krystian Iłkiewicz, Jarrod R. Hurley, Mirek Giersz, Felipe Lagos

Abstract: The predicted and observed space density of cataclysmic variables (CVs) have been for a long time discrepant by at least an order of magnitude. The standard model of CV evolution predicts that the vast majority of CVs should be period bouncers, whose space density has been recently measured to be $ρろー\lesssim 2 \times 10^{-5}$ pc$^{-3}$. We performed population synthesis of CVs using an updated vers… ▽ More The predicted and observed space density of cataclysmic variables (CVs) have been for a long time discrepant by at least an order of magnitude. The standard model of CV evolution predicts that the vast majority of CVs should be period bouncers, whose space density has been recently measured to be $ρろー\lesssim 2 \times 10^{-5}$ pc$^{-3}$. We performed population synthesis of CVs using an updated version of the Binary Stellar Evolution (BSE) code for single and binary star evolution. We find that the recently suggested empirical prescription of consequential angular momentum loss (CAML) brings into agreement predicted and observed space densities of CVs and period bouncers. To progress with our understanding of CV evolution it is crucial to understand the physical mechanism behind empirical CAML. Our changes to the BSE code are also provided in details, which will allow the community to accurately model mass transfer in interacting binaries in which degenerate objects accrete from low-mass main-sequence donor stars. △ Less

Submitted 14 June, 2018; originally announced June 2018.

Comments: 11 pages, 4 figures, 1 table; accepted for publication in MNRAS

The origin of the first neutron star -- neutron star merger

Authors: K. Belczynski, A. Askar, M. Arca-Sedda, M. Chruslinska, M. Donnari, M. Giersz, M. Benacquista, R. Spurzem, D. Jin, G. Wiktorowicz, D. Belloni

Abstract: The first neutron star-neutron star (NS-NS) merger was discovered on August 17, 2017 through gravitational waves (GW170817) and followed with electromagnetic observations. This merger was detected in an old elliptical galaxy with no recent star formation. We perform a suite of numerical calculations to understand the formation mechanism of this merger. We probe three leading formation mechanisms o… ▽ More The first neutron star-neutron star (NS-NS) merger was discovered on August 17, 2017 through gravitational waves (GW170817) and followed with electromagnetic observations. This merger was detected in an old elliptical galaxy with no recent star formation. We perform a suite of numerical calculations to understand the formation mechanism of this merger. We probe three leading formation mechanisms of double compact objects: classical isolated binary star evolution, dynamical evolution in globular clusters and nuclear cluster formation to test whether they are likely to produce NS-NS mergers in old host galaxies. Our simulations with optimistic assumptions show current NS-NS merger rates at the level of 10^-2 yr^-1 from binary stars, 5 x 10^-5 yr^-1 from globular clusters and 10^-5 yr^-1 from nuclear clusters for all local elliptical galaxies (within 100 Mpc^3). These models are thus in tension with the detection of GW170817 with an observed rate 1.5 yr^-1 (per 100 Mpc^3; LIGO/Virgo estimate). Our results imply that either (i) the detection of GW170817 by LIGO/Virgo at their current sensitivity in an elliptical galaxy is a statistical coincidence; or that (ii) physics in at least one of our three models is incomplete in the context of the evolution of stars that can form NS-NS mergers; or that (iii) another very efficient (unknown) formation channel with a long delay time between star formation and merger is at play. △ Less

Submitted 28 March, 2018; v1 submitted 2 December, 2017; originally announced December 2017.

Comments: A&A: accepted

Dynamical equivalence, the origin of the Galactic field stellar and binary population, and the initial radius--mass relation of embedded clusters

Authors: Diogo Belloni, Pavel Kroupa, Helio J. Rocha-Pinto, Mirek Giersz

Abstract: In order to allow a better understanding of the origin of Galactic field populations, dynamical equivalence of stellar-dynamical systems has been postulated by Kroupa and Belloni et al. to allow mapping of solutions of the initial conditions of embedded clusters such that they yield, after a period of dynamical processing, the Galactic field population. Dynamically equivalent systems are defined t… ▽ More In order to allow a better understanding of the origin of Galactic field populations, dynamical equivalence of stellar-dynamical systems has been postulated by Kroupa and Belloni et al. to allow mapping of solutions of the initial conditions of embedded clusters such that they yield, after a period of dynamical processing, the Galactic field population. Dynamically equivalent systems are defined to initially and finally have the same distribution functions of periods, mass ratios and eccentricities of binary stars. Here we search for dynamically equivalent clusters using the {\sc mocca} code. The simulations confirm that dynamically equivalent solutions indeed exist. The result is that the solution space is next to identical to the radius--mass relation of Marks \& Kroupa, $\left( r_h/{\rm pc} \right)= 0.1^{+0.07}_{-0.04}\, \left( M_{\rm ecl}/{\rm M}_\odot \right)^{0.13\pm0.04}$. This relation is in good agreement with the observed density of molecular cloud clumps. According to the solutions, the time-scale to reach dynamical equivalence is about 0.5~Myr which is, interestingly, consistent with the lifetime of ultra-compact HII regions and the time-scale needed for gas expulsion to be active in observed very young clusters as based on their dynamical modelling. △ Less

Submitted 2 December, 2017; v1 submitted 21 November, 2017; originally announced November 2017.

Comments: 6 pages, 2 figures; accepted for publication in MNRAS, minor changes made to the text after proofs

On the initial binary population for star cluster simulations

Authors: Diogo Belloni, Abbas Askar, Mirek Giersz, Pavel Kroupa, Helio J. Rocha-Pinto

Abstract: Colour-magnitude diagrams (CMDs) are powerful tools that might be used to infer stellar properties in globular clusters (GCs), for example, the binary fraction and their mass ratio ($q$) distribution. In the past few years, observations have revealed that q distributions of GC main-sequence (MS) binaries are generally flat, and a distribution characterized by a strong increase towards q $\approx$… ▽ More Colour-magnitude diagrams (CMDs) are powerful tools that might be used to infer stellar properties in globular clusters (GCs), for example, the binary fraction and their mass ratio ($q$) distribution. In the past few years, observations have revealed that q distributions of GC main-sequence (MS) binaries are generally flat, and a distribution characterized by a strong increase towards q $\approx$ 1 is not typical in GCs. In numerical simulations of GC evolution with the initial binary population (IBP) described by Kroupa, synthetic CMD colour distributions exhibit a peak associated with binaries that have q $\approx$ 1. While the Kroupa IBP reproduces binary properties in star-forming regions, clusters and the Galactic field, the peak in the q distribution towards q $\approx$ 1 observed for GC simulations is not consistent with distributions derived from observations. The objective of this paper is to refine and further improve the physical formulation of pre-main-sequence eigenevolution proposed by Kroupa in order to achieve CMD colour distributions of simulated GC models similar to those observed in real GCs, and to get a similarly good agreement with binary properties for late-type binaries in the Galactic field. We present in this paper a modified Kroupa IBP, in which early-type stars follow observational distributions, and late-type stars are generated according to slightly modified pre-main-sequence eigenevolution prescriptions. Our modifications not only lead to a qualitatively good agreement with respect to long-term observations of late-type binaries in the Galactic field, but also resolve the above-mentioned problem related to binary distributions in GC models. △ Less

Submitted 13 July, 2017; originally announced July 2017.

Comments: 18 pages, 10 figures, 4 tables; accepted for publication in MNRAS

MOCCA-SURVEY database I. Accreting white dwarf binary systems in globular clusters -- III. Cataclysmic variables -- Implications of model assumptions

Authors: Diogo Belloni, Mónica Zorotovic, Matthias R. Schreiber, Nathan W. C. Leigh, Mirek Giersz, Abbas Askar

Abstract: In this third of a series of papers related to cataclysmic variables (CVs) and related objects, we analyse the population of CVs in a set of 12 globular cluster models evolved with the MOCCA Monte Carlo code, for two initial binary populations (IBPs), two choices of common-envelope phase (CEP) parameters, and three different models for the evolution of CVs and the treatment of angular momentum los… ▽ More In this third of a series of papers related to cataclysmic variables (CVs) and related objects, we analyse the population of CVs in a set of 12 globular cluster models evolved with the MOCCA Monte Carlo code, for two initial binary populations (IBPs), two choices of common-envelope phase (CEP) parameters, and three different models for the evolution of CVs and the treatment of angular momentum loss. When more realistic models and parameters are considered, we find that present-day cluster CV duty cycles are extremely-low ($\lesssim 0.1$ per cent) which makes their detection during outbursts rather difficult. Additionally, the IBP plays a significant role in shaping the CV population properties, and models that follow the Kroupa IBP are less affected by enhanced angular momentum loss. We also predict from our simulations that CVs formed dynamically in the past few Gyr (massive CVs) correspond to bright CVs (as expected), and that faint CVs formed several Gyr ago (dynamically or not) represent the overwhelming majority. Regarding the CV formation rate, we rule out the notion that it is similar irrespective of the cluster properties. Finally, we discuss the differences in the present-day CV properties related to the IBPs, the initial cluster conditions, the CEP parameters, formation channels, the CV evolution models, and the angular momentum loss treatments. △ Less

Submitted 6 March, 2017; originally announced March 2017.

Comments: 18 pages, 7 figures, 2 tables; accepted for publication in MNRAS

Dynamical Formation of Cataclysmic Variables in Globular Clusters

Authors: Jongsuk Hong, Enrico Vesperini, Diogo Belloni, Mirek Giersz

Abstract: The formation and evolution of X-ray sources in globular clusters is likely to be affected by the cluster internal dynamics and the stellar interactions in the cluster dense environment.Several observational studies have revealed a correlation between the number of X-ray sources and the stellar encounter rate and provided evidence of the role of dynamics in the formation of X-ray binaries. We have… ▽ More The formation and evolution of X-ray sources in globular clusters is likely to be affected by the cluster internal dynamics and the stellar interactions in the cluster dense environment.Several observational studies have revealed a correlation between the number of X-ray sources and the stellar encounter rate and provided evidence of the role of dynamics in the formation of X-ray binaries. We have performed a survey of Monte-Carlo simulations aimed at exploring the connection between the dynamics and formation of cataclysmic variables (CVs) and the origin of the observed correlation between the number of these objects, $N_{\rm cv}$, and the stellar encounter rate, $Γがんま$.The results of our simulations show a correlation between $N_{\rm cv}$ and $Γがんま$ as found in observational data, illustrate the essential role played by dynamics, and shed light on the dynamical history behind this correlation. CVs in our simulations are more centrally concentrated than single stars with masses close to those of turn-off stars, although this trend is stronger for CVs formed from primordial binaries undergoing exchange encounters, which include a population of more massive CVs absent in the group of CVs formed from binaries not suffering any component exchange. △ Less

Submitted 2 November, 2016; originally announced November 2016.

Comments: 6 pages, 5 figures. Accepted for publication in MNRAS

MOCCA-SURVEY database I. Accreting white dwarf binary systems in globular clusters -- II. Cataclysmic variables -- progenitors and population at birth

Authors: Diogo Belloni, Mirek Giersz, Helio J. Rocha-Pinto, Nathan Leigh, Abbas Askar

Abstract: This is the second in a series of papers associated with cataclysmic variables (CVs) and related objects, formed in a suite of simulations for globular cluster evolution performed with the MOCCA Monte Carlo code. We study the properties of our simulated CV populations throughout the entire cluster evolution. We find that dynamics extends the range of binary CV progenitor properties, causing CV for… ▽ More This is the second in a series of papers associated with cataclysmic variables (CVs) and related objects, formed in a suite of simulations for globular cluster evolution performed with the MOCCA Monte Carlo code. We study the properties of our simulated CV populations throughout the entire cluster evolution. We find that dynamics extends the range of binary CV progenitor properties, causing CV formation from binary progenitors that would otherwise not become CVs. The CV formation rate in our simulations can be separated into two regimes: an initial burst ($\lesssim$ 1 Gyr) connected with the formation of the most massive white dwarfs, followed by a nearly constant formation rate. This result holds for all models regardless of the adopted initial conditions, even when most CVs form dynamically. Given the cluster age-dependence of CV properties, we argue that direct comparisons to observed Galactic field CVs could be misleading, since cluster CVs can be up to 4 times older than their field counterparts. Our results also illustrate that, due mainly to unstable mass transfer, some CVs that form in our simulations are destroyed before the present-day. Finally, some field CVs might have originated from globular clusters, as found in our simulations, although the fraction of such escapers should be small relative to the entire Galactic field CV population. △ Less

Submitted 2 December, 2016; v1 submitted 8 September, 2016; originally announced September 2016.

Comments: 19 pages, 7 figures, 2 tables; minor changes made to the text, accepted for publication in MNRAS

Cataclysmic variables in Globular clusters: First results on the analysis of the MOCCA simulations database

Authors: Diogo Belloni, Mirek Giersz, Abbas Askar, Arkadiusz Hypki

Abstract: In this first investigation of the MOCCA database with respect to cataclysmic variables, we found that for models with Kroupa initial distributions, considering the standard value of the efficiency of the common-envelope phase adopted in BSE, no single cataclysmic variable was formed only via binary stellar evolution, i. e., in order to form them, strong dynamical interactions have to take place.… ▽ More In this first investigation of the MOCCA database with respect to cataclysmic variables, we found that for models with Kroupa initial distributions, considering the standard value of the efficiency of the common-envelope phase adopted in BSE, no single cataclysmic variable was formed only via binary stellar evolution, i. e., in order to form them, strong dynamical interactions have to take place. Our results also indicate that the population of cataclysmic variables in globular clusters are, mainly, in the last stage of their evolution and observational selection effects can change drastically the expected number and properties of observed cataclysmic variables. △ Less

Submitted 26 July, 2016; originally announced July 2016.

Comments: 4 pages, 3 figures. Presented at the MODEST 16/Cosmic Lab conference in Bologna, Italy, April 18-22 2016. To be pusblished in Mem. S. A. It. Conference Series

MOCCA-SURVEY database I. Accreting white dwarf binary systems in globular clusters -- I. cataclysmic variables -- present-day population

Authors: Diogo Belloni, Mirek Giersz, Abbas Askar, Nathan Leigh, Arkadiusz Hypki

Abstract: In this paper, which is the first in a series of papers associated with cataclysmic variables and related objects, we introduce the CATUABA code, a numerical machinery written for analysis of the MOCCA simulations, and show some first results by investigating the present-day population of cataclysmic variables in globular clusters. Emphasis was given on their properties and the observational selec… ▽ More In this paper, which is the first in a series of papers associated with cataclysmic variables and related objects, we introduce the CATUABA code, a numerical machinery written for analysis of the MOCCA simulations, and show some first results by investigating the present-day population of cataclysmic variables in globular clusters. Emphasis was given on their properties and the observational selection effects when observing and detecting them. In this work we analysed in this work six models, including three with Kroupa distributions of the initial binaries. We found that for models with Kroupa initial distributions, considering the standard value of the efficiency of the common envelope phase adopted in BSE, no single cataclysmic variable was formed only via binary stellar evolution, i. e., in order to form them, strong dynamical interactions have to take place. We show and explain why this is inconsistent with observational and theoretical results. Our results indicate that the population of cataclysmic variables in globular clusters is, mainly, in the last stage of their evolution and observational selection effects can drastically change the expected number of observed cataclysmic variables. We show that the probability of observing them during the outbursts is extremely small and conclude that the best way of looking for cataclysmic variables in globular clusters is by searching for variabilities during quiescence, instead of during outbursts. For that, one would need a very deep observation which could reach magnitudes $\gtrsim$ 27 mag. Finally, we argue that cataclysmic variables in globular clusters are not necessarily magnetic. △ Less

Submitted 11 August, 2016; v1 submitted 26 July, 2016; originally announced July 2016.

Comments: 21 pages, 5 figures, 2 tables; minor changes made to the text, accepted for publication in MNRAS