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 $\omega$ Cen (SOPS IV e-94) and showed that this object is most likely not a SySt.

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.

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.

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$~\beta$ 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$~\beta$ 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.

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.

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).

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 $\rho \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.

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.

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.

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.

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.

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.

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, $\Gamma$.The results of our simulations show a correlation between $N_{\rm cv}$ and $\Gamma$ 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.

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.

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.