Lucas Hellström

Current gravitational wave detectors are sensitive to coalescing black holes and neutron stars. However, double white dwarfs (DWDs) have long been recognized as promising sources of gravitational waves, and upcoming detectors such as LISA will be able to observe these systems in abundance. Double white dwarfs are expected to be the dominant gravitational wave (GW) sources in parts of the LISA frequency range, making it crucial to understand their formation for future detections. The Milky Way contains many white dwarfs (WDs) in both the field and star clusters, promising a rich population of DWDs for LISA. However, the large number of sources may make it difficult to resolve individual binaries, and DWDs in the field and clusters often have similar properties, complicating the identification of their origins from GW signals alone. In this work, we focus on eccentric and tight DWDs, which cannot form in the field, but require dynamical interactions in dense clusters to increase their eccentricity after circularization through mass transfer phases and common-envelope evolution during binary evolution. These binaries may also form in three- and four-body dynamical interactions, whereby a DWD binary may directly form with high eccentricity and low separation. Our results show that we should expect eccentric and tight DWDs in clusters that can provide a meaningful GW signal; however, the number is low, with an upper limit of 10-15 in the MW. These can be used to independently obtain the distances of their host clusters.

Using the MOCCA code, we study the evolution of globular clusters (GCs) with multiple stellar populations. For this purpose, the MOCCA code has been significantly extended to take into account the formation of an enriched population of stars from re-accreted gas with a time delay after the formation of the pristine population of stars. The possibility of cluster migration in the host galaxy and the fact that the pristine population can be described by a model not in virial equilibrium are also taken into account. Gas re-accretion and cluster migration have a decisive impact on the observational parameters of clusters and the ratio of the number of objects between the pristine and enriched populations. The obtained results, together with observational data, suggest a speculative refinement of the AGB scenario that makes it possible to explain some observational data, such as the ratio of the pristine to the enriched populations, the observational fact that for some GCs the pristine population is more concentrated than the enriched one, and possibly a correlation between the ratio of the number of enriched stars to the total number of stars and the mass of the cluster. In this scenario, it is important to take into account the environment in which the cluster lives, the conditions in the galaxy when it formed, and the fact that a significant part of the GCs associated with the Galaxy come from dwarf galaxies that merged with the Milky Way. The initial conditions of GCs in our simulations differ from the widely used typical models, as they require GCs to fill the Roche lobe rather than being highly concentrated within it, imposing strong constraints on their formation locations within the galaxy.

Current gravitational wave detectors are sensitive to coalescing black holes and neutron stars. However, double white dwarfs (DWDs) have long been recognized as promising sources of gravitational waves, and upcoming detectors like LISA will be able to observe these systems in abundance. DWDs are expected to be the dominant gravitational wave (GW) sources in parts of the LISA frequency range, making it crucial to understand their formation for future detections. The Milky Way contains many white dwarfs (WDs) in both the field and star clusters, promising a rich population of DWDs for LISA. However, the large number of sources may make it difficult to resolve individual binaries, and DWDs in the field and clusters often have similar properties, complicating the identification of their origins from GW signals alone. In this work, we focus on eccentric and tight DWDs, which cannot form in the field, but require dynamical interactions in dense clusters to increase their eccentricity after circularization through mass transfer phases and common-envelope evolution during binary evolution. These binaries may also form in three- and four-body dynamical interactions where a DWD binary may directly form with high eccentricity and low separation. Our results show that we should expect eccentric and tight DWDs in clusters that can provide meaningful GW signal, however, the number is low; with an upper limit of 10-15 in the MW. These can be used to independently obtain distances of their host cluster.

Recent findings show that, in some Milky Way globular clusters (GCs), pristine red giant branch (RGB) stars are more centrally concentrated than enriched ones. This contradicts most multiple-population formation scenarios, which predict that the enriched population (2P) should initially be more concentrated than the pristine population (1P). We analyzed a MOCCA GC model that exhibits a higher spatial concentration of 1P RGB stars than 2P RGB stars at 13 Gyr. The MOCCA models assume the asymptotic giant branch (AGB) scenario, in which 2P stars are initially more concentrated than 1P stars. Our results indicate that the observed spatial distributions of multiple populations, and possibly their kinematics, are significantly shaped by dynamical interactions. These interactions preferentially eject 2P RGB progenitors from the central regions, leading to a transient overconcentration of 1P RGB stars at late times. This effect is particularly relevant for GCs with present-day masses of a few 10⁵ M_⊙, which have retained only about 10–20% of their initial mass. Such clusters may appear dynamically young due to heating from a black hole subsystem, even if they have undergone significant mass loss and dynamical evolution. Additionally, the relatively small number of RGB stars in these clusters suggests that interpreting the spatial distributions of multiple populations solely from RGB stars may lead to biased conclusions about the overall distribution of 1P and 2P. The apparent overconcentration of the 1P relative to the 2P is likely a transient effect driven by the preferential removal of 2P RGB progenitors via strong dynamical encounters. MOCCA models of multiple stellar populations based on the AGB scenario may explain anomalous features observed in some Galactic GCs, such as NGC 3201 and NGC 6101.

We explored the evolution of various properties of multiple-population globular clusters (GCs) for a broad range of initial conditions. We simulated over 200 GC models using the MOCCA Monte Carlo code and find that the present-day properties of core and half-light radii and the ratio of the number of second-generation (SG) stars to the total number of stars (N_SG/N_TOT) of these models cover the observed values of these quantities for Milky Way GCs. Starting with a relatively small value of the SG fraction (N_SG/N_TOT ~ 0.25) and a SG system concentrated in the inner regions of the cluster, we find, in agreement with previous studies, that systems in which the first-generation (FG) is initially tidally filling or slightly tidally underfilling best reproduce the observed ratios of N_SG/N_TOT and have values of the core and half-light radii typical of those of many Galactic globular clusters. Models in which the FG is initially tidally underfilling retain values of N_SG /N_TOT close to their initial values. These simulations expand previous investigations and serve to further constrain the viable range of initial parameters and better understand their influence on present-day GC properties. The results of this investigation also provide the basis for our future survey aimed at building specific models to reproduce the observed trends (or lack thereof) between the properties of multiple stellar populations and other cluster properties.

There could be a significant population of double white dwarf binaries (DWDs) inside globular clusters (GCs); however, these binaries 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 the other 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 analyzed three main groups of DWDs, namely in-cluster binaries, escaped binaries, and binaries formed from the 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 primordial 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. Even though the current observational samples of DWDs are extremely biased and incomplete, we conclude that our results compare reasonably well with observations.

In the cores of dense stellar clusters, close gravitational encounters between binary and single stars can frequently occur. Using the tsunami code, we computed the outcome of a large number of binary–single interactions involving two black holes (BHs) and a star to check how the inclusion of orbital energy losses due to tidal dissipation can change the outcome of these chaotic interactions. Each interaction was first simulated without any dissipative processes and then we systematically added orbital energy losses due to gravitational wave emission [using post-Newtonian (PN) corrections] and dynamical tides and recomputed the interactions. We find that the inclusion of tides increases the number of BH–star mergers by up to 75 per cent; however, it does not affect the number of BH–BH mergers. These results highlight the importance of including orbital energy dissipation due to dynamical tides during few-body encounters and evolution of close binary systems within stellar cluster simulations. Consistent with previous studies, we find that the inclusion of PN terms increases the number of BH–BH mergers during binary–single encounters. However, BH–star mergers are largely unaffected by the inclusion of these terms.

In the cores of dense stellar clusters, close gravitational encounters between binary and single stars can frequently occur. Using the Tsunami code, we computed the outcome of a large number of binary-single interactions involving two black holes (BHs) and a star to check how the inclusion of orbital energy losses due to tidal dissipation can change the outcome of these chaotic interactions. Each interaction was first simulated without any dissipative processes and then we systematically added orbital energy losses due to gravitational wave emission (using post-Newtonian (PN) corrections) and dynamical tides and recomputed the interactions. We find that the inclusion of tides increases the number of BH-star mergers by up to 75 per cent but it does not affect the number of BH-BH mergers. These results highlight the importance of including orbital energy dissipation due to dynamical tides during few-body encounters and evolution of close binary systems within stellar cluster simulations. Consistent with previous studies, we find that the inclusion of PN terms increases the number of BH-BH mergers during binary-single encounters. However, BH-star mergers are largely unaffected by the inclusion of these terms.