Many aspects of the evolution of stars, and in particular the evolution of binary stars, remain beyond our ability to model them in detail. Instead, we rely on observations to guide our often phenomenological models and pin down uncertain model parameters. To do this statistically requires population synthesis. Populations of stars modelled on computers are compared to populations of stars observed with our best telescopes. The closest match between observations and models provides insight into unknown model parameters and hence the underlying astrophysics. In this brief review, we describe the impact that modern big-data surveys will have on population synthesis, the large parameter space problem that is rife for the application of modern data science algorithms, and some examples of how population synthesis is relevant to modern astrophysics.

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

The combination of asteroseismologically-measured masses with abundances from detailed analyses of stellar atmospheres challenges our fundamental knowledge of stars and our ability to model them. Ancient red-giant stars in the Galactic thick disc are proving to be most troublesome in this regard. They are older than 5 Gyr, a lifetime corresponding to an initial stellar mass of about $1.2{\mathrm{M}_{\odot}}$. So why do the masses of a sizeable fraction of thick-disc stars exceed $1.3{\mathrm{M}_{\odot}}$, with some as massive as $2.3{\mathrm{M}_{\odot}}$ ? We answer this question by considering duplicity in the thick-disc stellar population using a binary population-nucleosynthesis model. We examine how mass transfer and merging affect the stellar mass distribution and surface abundances of carbon and nitrogen. We show that a few per cent of thick-disc stars can interact in binary star systems and become more massive than $1.3{\mathrm{M}_{\odot}}$. Of these stars, most are single because they are merged binaries. Some stars more massive than $1.3{\mathrm{M}_{\odot}}$ form in binaries by wind mass transfer. We compare our results to a sample of the APOKASC data set and find reasonable agreement except in the number of these thick-disc stars more massive than $1.3{\mathrm{M}_{\odot}}$. This problem is resolved by the use of a logarithmically-flat orbital-period distribution and a large binary fraction.

Context. Classical novae are thermonuclear explosions that take place in the envelopes of accreting white dwarfs in stellar binary systems. The material transferred onto the white dwarf piles up under degenerate conditions, driving a thermonuclear runaway. In those outbursts, about 10-7 - 10-3 Msun, enriched in CNO and, sometimes, other intermediate-mass elements (e.g., Ne, Na, Mg, or Al, for ONe novae) are ejected into the interstellar medium. The large concentrations of metals spectroscopically inferred in the nova ejecta reveal that the (solar-like) material transferred from the secondary mixes with the outermost layers of the underlying white dwarf. Aims. Most theoretical models of nova outbursts reported to date yield, on average, outflows characterized by O > C, from which only oxidized condensates (e.g, O-rich grains) would be expected, in principle. Methods. To specifically address whether CO novae can actually produce C-rich dust, six different hydrodynamic nova models have been evolved, from accretion to the expansion and ejection stages, with different choices for the composition of the substrate with which the solar-like accreted material mixes. Updated chemical profiles inside the H-exhausted core have been used, based on stellar evolution calculations for a progenitor of 8 Msun through H and He-burning phases. Results. We show that these profiles lead to C-rich ejecta after the nova outburst. This extends the possible contribution of novae to the inventory of presolar grains identified in meteorites, particularly in a number of carbonaceous phases (i.e., nanodiamonds, silicon carbides and graphites).

A grid of evolutionary sequences of stars in the mass range $1.2$-$7$ M$_{\odot}$, with solar-like initial composition is presented. We focus on this mass range in order to estimate the masses and calculate the CNO surface abundances of a sample of observed red giants. The stellar models are calculated from the zero-age main sequence till the early asymptotic giant branch (AGB) phase. Stars of M $\leqslant$ $2.2$M$_{\odot}$ are evolved through the core helium flash. In this work, an approach is adopted that improves the mass determination of an observed sample of 21 RGB and early AGB stars. This approach is based on comparing the observationally derived effective temperatures and absolute magnitudes with the calculated values based on our evolutionary tracks in the Hertzsprung-Russell diagram. A more reliable determination of the stellar masses is achieved by using evolutionary tracks extended to the range of observation. In addition, the predicted CNO surface abundances are compared to the observationally inferred values in order to show how far standard evolutionary calculation can be used to interpret available observations and to illustrate the role of convective mixing. We find that extra mixing beyond the convective boundary determined by the Schwarzschild criterion is needed to explain the observational oxygen isotopic ratios in low mass stars. The effect of recent determinations of proton capture reactions and their uncertainties on the $^{16}$O$/^{17}$O and $^{14} $N$/^{15}$N ratios is also shown. It is found that the $^{14}$N$($ p$,\gamma)^{15}$O reaction is important for predicting the $^{14}$N$/^{15}$N ratio in red giants.

Stars in the mass range ~8 - 12 $M_{\odot }$ are the most numerous massive stars. This mass range is critical because it may lead to supernova (SN) explosion, so it is important for the production of heavy elements and the chemical evolution of the galaxy. We investigate the critical transition mass ($M_{up}$), which is the minimum initial stellar mass that attains the conditions for hydrostatic carbon burning. Stars of masses < $M_{up}$ evolve to the Asymptotic Giant Branch and then develop CO White Dwarfs, while stars of masses $\geqslant $ $M_{up}$ ignite carbon in a partially degenerate CO core and form electron degenerate ONe cores. These stars evolve to the Super AGB (SAGB) phase and either become progenitors of ONe White Dwarfs or eventually explode as electron-capture SN (EC-SN). We study the sensitivity of $M_{up}$ to the C-burning reaction rate and to the treatment of convective mixing. In particular, we show the effect of a recent determination of the $^{12}C+$ $^{12}C$ fusion rate, as well as the extension of the convective core during hydrogen and helium burning on $M_{up}$ in solar metallicity stars. We choose the 9$M_{\odot }$ model to show the detailed characteristics of the evolution with the new C-burning rate.

A systematic investigation on the third dredge up in a 3M$_{\odot}$, solar metallicity AGB star will be presented. The model evolves from the main sequence up to the Asymptotic Giant Branch (AGB). Intermediate mass stars are important because they contribute significantly via the slow neutron capture nucleosynthesis. The aim of this work is to gain insight on the behaviour of the AGB star during thermal pulsation. This investigation is based on an extended numerical simulation of the evolutionary phases and full, consistent AGB model calculations. In particular, the convective structure during pulsation will be studied, giving particular emphasis to the analysis of the stability of the Schwarzschild boundary that will eventually determine the occurrence of Third Dredge Up (hereafter referred to as TDUP). We provide a brief description of our updated evolutionary code and focus primarily on the obtaining the TDUP after 14 thermal pulses. We elaborate on the non-standard treatment of convection known as "overshooting" and its efficiency in achieving the TDUP event associated with the thermal pulsations (helium shell flashes). However, our knowledge of the physics behind the TDUP event, and the theoretical description of the overshooting mechanism suffer from embarrassing gaps and are, therefore, not properly treated in stellar modelling. This is not surprising, particularly with the oversimplified descriptions of important mixing mechanisms like convection in stellar models. However, since such event drives the $^{13}$C nucleosynthesis of interest for the following s-process activation, we present a full assessment of this process and its implications as calculations will further proceed in different directions.

We investigate the effects of a modification of the $^{14}$N(p,{\gamma})$^{15}$O reaction rate, as suggested by recent evaluations, on the formation and extension of the blue loops encountered during the evolution of the stars in the mass range 5M$_{\odot}$ to 12M$_{\odot}$. We show that the blue loops of stars in the mass range 5M$_{\odot}$ to 8M$_{\odot}$, that is the range of super ABG stars, are severely affected by a modification of the important $^{14}$N(p,{\gamma})$^{15}$O reaction rate. We also show that the blue loops can be restored if envelope overshooting is included, which is necessary to explain the observations of the Cepheid stars.

We present the evolution of a 6M$_{\odot}$ star, of solar-like initial metallicity, and investigate the effects of key nuclear reaction rates, as well as the treatment of the convective mixing on its evolution along the Cepheid instability strip. In particular, we study the effect of recent estimates of the $^{14}$N(p,{\gamma})$^{15}$O reaction on the formation and extension of the blue loop during core helium burning. We also investigate the effects induced on this blue loop by the adoption of non-standard convective mixing prescriptions, as well as the implications of modifying the Mixing Length Theory.