In this paper, we consider triple systems composed of main-sequence (MS) stars, and their internal evolution due to stellar and binary evolution. Our focus is on triples that produce white dwarfs (WDs), where Roche lobe overflow of an evolving tertiary triggers accretion onto the inner binary via a circumbinary disk (CBD) driving it toward a mass ratio of unity. We present a combination of analytic- and population synthesis-based calculations performed using the \texttt{SeBa} code to constrain the expected frequency of such systems, given a realistic initial population of MS triples, and provide the predicted distributions of orbital periods. We identify the parameter space for triples that can accommodate a CBD, to inform future numerical simulations of suitable initial conditions. We find that $\lesssim$ 10\% of all MS triples should be able to accommodate a CBD around the inner binary, and compute lower limits for the production rates. This scenario broadly predicts mergers of near equal-mass binaries, producing blue stragglers (BSs), Type Ia supernovae, gamma ray bursts and gravitational wave-induced mergers, along with the presence of an outer WD tertiary companion. We compare our predicted distributions to a sample of field BS binaries, and argue that our proposed mechanism explains the observed range of orbital periods. Finally, the mechanism considered here could produce hypervelocity MS stars, WDs and even millisecond pulsars with masses close to the Chandrasekhar mass limit, and be used to constrain the maximum remnant masses at the time of any supernova explosion.

We study chaos and L\'evy flights in the general gravitational three-body problem. We introduce new metrics to characterize the time evolution and final lifetime distributions, namely Scramble Density $\mathcal{S}$ and the LF index $\mathcal{L}$, that are derived from the Agekyan-Anosova maps and homology radius $R_{\mathcal{H}}$. Based on these metrics, we develop detailed procedures to isolate the ergodic interactions and L\'evy flight interactions. This enables us to study the three-body lifetime distribution in more detail by decomposing it into the individual distributions from the different kinds of interactions. We observe that ergodic interactions follow an exponential decay distribution similar to that of radioactive decay. Meanwhile, L\'evy flight interactions follow a power-law distribution. L\'evy flights in fact dominate the tail of the general three-body lifetime distribution, providing conclusive evidence for the speculated connection between power-law tails and L\'evy flight interactions. We propose a new physically-motivated model for the lifetime distribution of three-body systems and discuss how it can be used to extract information about the underlying ergodic and L\'evy flight interactions. We discuss mass ejection probabilities in three-body systems in the ergodic limit and compare it to previous ergodic formalisms. We introduce a novel mechanism for a three-body relaxation process and discuss its relevance in general three-body systems.

Runaway stellar collisions in dense star clusters are invoked to explain the presence of very massive stars or blue stragglers in the center of those systems. This process has also been explored for the first star clusters in the Universe and shown to yield stars that may collapse at some points into an intermediate mass black hole. Although the early evolution of star clusters requires the explicit modeling of the gas out of which the stars form, these calculations would be extremely time-consuming and often the effects of the gas can be accurately treated by including a background potential to account for the extra gravitational force. We apply this approximation to model the early evolution of the first dense star clusters formed in the Universe by performing $N$-body simulations, our goal is to understand how the additional gravitational force affects the growth of a very massive star through stellar mergers in the central parts of the star cluster. Our results show that the background potential increases the velocities of the stars, causing an overall delay in the evolution of the clusters and in the runaway growth of a massive star at the center. The population of binary stars is lower due to the increased kinetic energy of the stars, initially reducing the number of stellar collisions, and we show that relaxation processes are also affected. Despite these effects, the external potential enhances the mass of the merger product by a factor $\sim$2 if the collisions are maintained for long times.

We study the stellar kinematic properties and spatial distribution of the association Cyg OB3 using precise astrometric data from Gaia DR2. All known O- and B-type stars in Cyg OB3 region with positions, parallaxes and proper motions available are included, comprising a total of 41 stars. The majority of stars are found to be concentrated at a heliocentric distance of 2.0 +/- 0.3 kpc. The mean peculiar velocity of the sample after removing Galactic rotation and solar motion is ~22 km/s, dominated by the velocity component towards the Galactic center. The relative position and velocity of the black hole X-ray binary Cyg X-1 with respect to the association suggest that Cyg OB3 is most likely its parent association. The peculiar kinematic properties of some of the stars are revealed and are suggestive of past stellar encounters. The sample includes a previously known runaway star HD 227018, and its high peculiar velocity of ~50 km/s is confirmed with Gaia. We estimated the velocities of stars relative to the association and the star HD 225577 exhibits peculiar velocity smaller than its velocity relative to the association. The star has lower value of proper motion than the rest of the sample. The results suggest a slowly expanding nature of the association, which is supported by the small relative speeds <20 km/s with respect to the association for a majority of the sample stars.

Mergers of binaries comprised of compact objects can give rise to explosive transient events, heralding the birth of exotic objects which cannot be formed through single star evolution. Using a large number of direct N-body simulations, we explore the possibility that a white dwarf (WD) is dynamically driven to tidal disruption by a stellar-mass black hole (BH) as a consequence of the joint effects of gravitational wave (GW) emission and Lidov-Kozai oscillations imposed by the tidal field of a outer tertiary companion orbiting the inner BH-WD binary. We explore the sensitivity of our results to the distributions of natal kick velocities imparted to the BH and WD upon formation, adiabatic mass loss, semi-major axes and eccentricities of the triples, and stellar mass ratios. We find rates of WD-TDEs in the range $1.2\times 10^{-3}-1.4$ Gpc$^{-3}$ yr$^{-1}$ for $z\leq 0.1$, rarer than stellar TDEs in triples by a factor of $\sim 3$--$30$. The uncertainty in the TDE rates may be greatly reduced in the future using gravitational wave (GW) observations of Galactic binaries and triples with LISA. WD-TDEs may give rise to high energy X-ray or gamma-ray transients of duration similar to long gamma-ray bursts but lacking the signatures of a core-collapse supernova, while being accompanied by a supernova-like optical transient which lasts for only days. WD--BH and WD--NS binaries will also emit GWs in the LISA band before the TDE. The discovery and identification of triple-induced WD-TDE events by future time domain surveys and/or GWs could enable the study of the demographics of BHs in nearby galaxies.

The mass, rate, and spins of merging stellar-mass binary black holes (BBHs) detected by aLIGO and Advanced Virgo provide challenges to traditional BBH formation and merger scenarios. An active galactic nucleus (AGN) disk provides a promising additional merger channel, because of the powerful influence of the gas that drives orbital evolution, makes encounters dissipative, and leads to migration. Previous work showed that stellar mass black holes (sBHs) in an AGN disk migrate to regions of the disk, known as migration traps, where positive and negative gas torques cancel out, leading to frequent BBH formation. Here we build on that work by simulating the evolution of additional sBHs that enter the inner disk by either migration or inclination reduction. We also examine whether the BBHs formed in our models have retrograde or prograde orbits around their centers of mass with respect to the disk, determining the orientation, relative to the disk, of the spin of the merged BBHs. We find that orbiters entering the inner disk form BBHs, that are commonly asymmetric in mass, with sBHs on resonant orbits near the migration trap. When these black holes reach 80 solar masses, they form BBHs with sBHs in the migration trap, which over 10 Myr reach about 1000 solar masses. We also find that 68% of the BBHs in our simulation orbit in the retrograde direction, which implies our merger channel could easily produce the low values of the dimensionless aligned spin, chi_eff, observed in a majority of BBH mergers detected thus far.

We present a novel method to constrain the past collisional evolution of observed globular cluster (GC) systems, in particular their mass functions. We apply our method to a pair of galaxies hypothesized to have recently undergone an episode of violent relaxation due to a strong galaxy-galaxy interaction, namely NGC 1052-DF2 and NGC 1052-DF4. We begin by exploring the observational evidence for a collisional origin for these two recently discovered ultra-diffuse galaxies observed in the NGC 1052 group, posited in the literature to be dark matter (DM)-free. We compute the timescales for infall to the central nucleus due to dynamical friction (DF) for the GCs in these galaxies, using the shortest of these times to constrain how long ago a galaxy-galaxy interaction could have occurred. We go on to quantify the initial GC numbers and densities needed for significant collisional evolution to occur within the allotted times, and show that, if the hypothesis of a previous galaxy-galaxy interaction is correct, a paucity of low-mass GCs should be revealed by deeper observational surveys. If any are found, they should be more spatially extended than the currently observed GC population. Finally, we apply our method to these galaxies, in order to illustrate its efficacy in constraining their dynamical evolution. Our results motivate more complete observations of the GC luminosity functions in these galaxies, in addition to future studies aimed at combining the method presented here with a suite of numerical simulations in order to further constrain the origins of the curious GC populations in these (and other) galaxies.

The discovery of planetary systems outside of the solar system has challenged some of the tenets of planetary formation. Among the difficult-to-explain observations, are systems with a giant planet orbiting a very-low mass star, such as the recently discovered GJ~3512b planetary system, where a Jupiter-like planet orbits an $M$-star in a tight and eccentric orbit. Systems such as this one are not predicted by the core accretion theory of planet formation. Here we suggest a novel mechanism, in which the giant planet is born around a more typical Sun-like star ($M_{*,1}$), but is subsequently exchanged during a dynamical interaction with a flyby low-mass star ($M_{*,2}$). We perform state-of-the-art $N$-body simulations with $M_{*,1}=1M_\odot$ and $M_{*,2}=0.1M_\odot$ to study the statistical outcomes of this interaction, and show that exchanges result in high eccentricities for the new orbit around the low-mass star, while about half of the outcomes result in tighter orbits than the planet had around its birth star. We numerically compute the cross section for planet exchange, and show that an upper limit for the probability per planetary system to have undergone such an event is $\Gamma\sim 4.4(M_{\rm c}/100M_\odot)^{-2}(a_{\rm p}/{\rm AU}) (\sigma/1\,{\rm km}\,{\rm s}^{-1})^{5}$Gyr$^{-1}$, where $a_{\rm p}$ is the planet semi-major axis around the birth star, $\sigma$ the velocity dispersion of the star cluster, and $M_{\rm c}$ the total mass of the star cluster. Hence these planet exchanges could be relatively common for stars born in open clusters and groups, should already be observed in the exoplanet database, and provide new avenues to create unexpected planetary architectures.

The discovery of Exoplanetary Systems has challenged some of the theories of planet formation, which assume unperturbed evolution of the host star and its planets. However, in star clusters the interactions with flyby stars and binaries may be relatively common during the lifetime of a planetary system. Here, via high-resolution $N$-body simulations of star-planet systems perturbed by interlopers (stars and binaries), we explore the reconfiguration to the planetary system due to the encounters. In particular, via an exploration focused on the strong scattering regime, we derive the fraction of encounters which result in planet ejections, planet transfers and collisions by the interloper star/binary, as a function of the characteristics of the environment (density, velocity dispersion), and for different masses of the flyby star/binary. We find that binary interlopers can significantly increase the cross section of planet ejections and collisions, while they only slightly change the cross section for planet transfers. Therefore, in environments with high binary fractions, floating planets are expected to be relatively common, while in environments with low binary fractions, where the cross sections of planet ejection and transfer are comparable, the rate of planet exchanges between two stars will be comparable to the rate of production of free-floating planets.

We study the dynamical evolution of eccentric massive black hole binaries (MBHBs) interacting with unbound stars by means of an extensive set of three body scattering experiments. Compared to previous studies, we extend the investigation down to a MBHB mass ratio of $q=m_2/m_1=10^{-4}$, where $m_1$ and $m_2$ are the masses of the primary and secondary hole respectively. Contrary to a simple extrapolation from higher mass ratios, we find that for $q\lesssim 10^{-3}$ the eccentricity growth rate becomes negative, i.e., the binary {\it circularises} as it shrinks. This behaviour is due to the subset of interacting stars captured in metastable counter-rotating orbits; those stars tend to extract angular momentum from the binary, promoting eccentricity growth for $q>10^{-3}$, but tend to inject angular momentum into the binary driving it towards circularisation for $q<10^{-3}$. The physical origin of this behaviour requires a detailed study of the orbits of this subset of stars and is currently under investigation. Our findings might have important consequences for intermediate MBHs (IMBHs) inspiralling onto MBHs (e.g. a putative $10^3\rm M_{\odot}$ black hole inspiralling onto SgrA$^*$).

The dynamical production of low-mass X-ray binaries and brighter cataclysmic variables (CVs) in dense globular clusters is well-established. We investigate how the X-ray emissivity of fainter X-ray binaries (principally CVs and coronally active binaries) varies between different environments. We compile calculations (largely from the literature) of the X-ray emissivity of old stellar populations, including open and globular clusters and several galaxies. We investigate three literature claims of unusual X-ray sources in low-density stellar populations. We show that a suggested quiescent neutron star in the open cluster NGC 6819 is a foreground M dwarf. We show that the suggested diffuse X-ray emission from an old nova shell in the globular cluster NGC 6366 is actually a background galaxy cluster. And we show that a suggested population of quiescent X-ray binaries in the Sculptor Dwarf Galaxy is mostly (perhaps entirely) background galaxies. We find that above densities of $10^4$ M$_{\odot}$/pc$^3$, the X-ray emissivity of globular clusters increases, due to dynamical production of X-ray emitting systems. Below this density, globular clusters have lower X-ray emissivity than the other populations, and we do not see a strong dependence of X-ray emissivity due to density effects. We find significant correlations between X-ray emissivity and binary fraction, metallicity, and density. Sampling these fits via bootstrap techniques gives less significant correlations, but confirms the effect of metallicity on low-density populations, and that of density on the full globular cluster sample.

The merger rate of stellar-mass black hole binaries (sBHBs) inferred by the Advanced Laser Interferometer Gravitational-Wave Observatory (LIGO) suggests the need for an efficient source of sBHB formation. Active galactic nucleus (AGN) disks are a promising location for the formation of these sBHBs, as well as binaries of other compact objects, because of powerful torques exerted by the gas disk. These gas torques cause orbiting compact objects to migrate towards regions in the disk where inward and outward torques cancel, known as migration traps. We simulate the migration of stellar mass black holes in an example of a model AGN disk, using an augmented N-body code that includes analytic approximations to migration torques, stochastic gravitational forces exerted by turbulent density fluctuations in the disk, and inclination and eccentricity dampening produced by passages through the gas disk, in addition to the standard gravitational forces between objects. We find that sBHBs form rapidly in our model disk as stellar-mass black holes migrate towards the migration trap. These sBHBs are likely to subsequently merge on short time-scales. The process continues, leading to the build-up of a population of over-massive stellar-mass black holes. The formation of sBHBs in AGN disks could contribute significantly to the sBHB merger rate inferred by LIGO.

We investigate the role of the Keplerian tidal field generated by a supermassive black hole (SMBH) on the three-body dynamics of stellar mass black holes. We consider two scenarios occurring close to the SMBH: the breakup of unstable triples and three-body encounters between a binary and a single. These two cases correspond to the hard and soft binary cases, respectively. The tidal field affects the breakup of triples by tidally limiting the system, so that the triples break earlier with lower breakup velocity, leaving behind slightly larger binaries (relative to the isolated case). The breakup direction becomes anisotropic and tends to follow the shape of the Hill region of the triple, favouring breakups in the radial direction. Furthermore, the tidal field can torque the system, leading to angular momentum exchanges between the triple and its orbit about the SMBH. This process changes the properties of the final binary, depending on the initial angular momentum of the triple. Finally, the tidal field also affects binary-single encounters: binaries tend to become both harder and more eccentric with respect to encounters that occur in isolation. Consequently, single-binary scattering in a deep Keplerian potential produces binaries with shorter gravitational wave merger timescales.

The three-body problem is arguably the oldest open question in astrophysics, and has resisted a general analytic solution for centuries. Various implementations of perturbation theory provide solutions in portions of parameter space, but only where hierarchies of masses or separations exist. Numerical integrations show that bound, non-hierarchical triples of Newtonian point particles will almost always disintegrate into a single escaping star and a stable, bound binary, but the chaotic nature of the three-body problem prevents the derivation of tractable analytic formulae deterministically mapping initial conditions to final outcomes. However, chaos also motivates the assumption of ergodicity, suggesting that the distribution of outcomes is uniform across the accessible phase volume. Here, we use the ergodic hypothesis to derive a complete statistical solution to the non-hierarchical three-body problem, one which provides closed-form distributions of outcomes (e.g. binary orbital elements) given the conserved integrals of motion. We compare our outcome distributions to large ensembles of numerical three-body integrations, and find good agreement, so long as we restrict ourselves to "resonant" encounters (the ~50% of scatterings that undergo chaotic evolution). In analyzing our scattering experiments, we identify "scrambles" (periods in time where no pairwise binaries exist) as the key dynamical state that ergodicizes a non-hierarchical triple. The generally super-thermal distributions of survivor binary eccentricity that we predict have notable applications to many astrophysical scenarios. For example, non-hierarchical triples produced dynamically in globular clusters are a primary formation channel for black hole mergers, but the rates and properties of the resulting gravitational waves depend on the distribution of post-disintegration eccentricities.

Recent studies show that the majority of blue straggler stars (BSSs) in old open clusters are formed through mass transfer from an evolved star onto a main-sequence companion, resulting in a BSS and white dwarf (WD) in a binary system. We present constraints on the mass transfer histories of two BSS-WD binaries in the open cluster NGC 188, using WD temperatures and surface gravities measured with HST COS far-ultraviolet spectroscopy. Adopting a Gaia-based cluster distance of $1847\pm107$ pc, we determine that one system, WOCS 4540, formed through Case C mass transfer resulting in a CO-core white dwarf with $T_{\text{eff}}=17000^{+140}_{-200}$ K and a $\log g=7.80^{+0.06}_{-0.06}$, corresponding to a mass of $0.53^{+0.03}_{-0.03}$ $M_{\odot}$ and a cooling age of $105^{+6}_{-5}$ Myr. The other system, WOCS 5379, formed through Case B mass transfer resulting in a He-core white dwarf with $T_{\text{eff}}=15500^{+170}_{-150}$ K and a $\log g=7.50^{+0.06}_{-0.05}$, corresponding to a mass of $0.42^{+0.02}_{-0.02}$ $M_{\odot}$ and an age of $250^{+20}_{-20}$ Myr. The WD parameters are consistent across four different cluster distance assumptions. We determine possible progenitor binary systems with a grid of accretion models using MESA, and investigate whether these systems would lead to stable or unstable mass transfer. WOCS 4540 likely resulted from stable mass transfer during periastron passage in an eccentric binary system, while WOCS 5379 challenges our current understanding of the expected outcomes for mass transfer from red giant branch stars. Both systems are examples of the value in using detailed analyses to fine-tune our physical understanding of binary evolutionary processes.

Stars passing too close to a black hole can produce tidal disruption events (TDEs), when the tidal force across the star exceeds the gravitational force that binds it. TDEs have usually been discussed in relation to massive black holes that reside in the centres of galaxies or lurk in star clusters. In this paper, we investigate the possibility that triple stars hosting a stellar black hole (SBH) may be sources of TDEs. We start from a triple system made up of three main sequence (MS) stars and model the supernova (SN) kick event that led to the production of an inner binary comprised of a SBH. We evolve these triples in isolation with a high precision $N$-body code and study their TDEs as a result of Kozai-Lidov oscillations. We explore a variety of distributions of natal kicks imparted during the SN event, various maximum initial separations for the triples, and different distributions of eccentricities. We show that the main parameter that governs the properties of the SBH-MS binaries which produce a TDE in triples is the mean velocity of the natal kick distribution. Smaller $\sigma$'s lead to larger inner and outer semi-major axes of the systems that undergo a TDE, smaller SBH masses, and longer timescales. We find that the fraction of systems that produce a TDE is roughly independent of the initial conditions, while estimate a TDE rate of $7.3\times 10^{-5}-4.1 \ \mathrm{yr}^{-1}$, depending on the prescriptions adopted for the SBH natal kicks. This rate is almost comparable to the expected TDE rate for massive black holes.

We investigate the old open cluster M67 using ultraviolet photometric data of Ultra-Violet Imaging Telescope in multi-filter far-UV bands. M67, well known for the presence of several blue straggler stars (BSS), has been put to detailed tests to understand their formation pathways. Currently, there are three accepted formation channels: mass transfer due to Roche-lobe overflow in binary systems, stellar mergers either due to dynamical collisions or through coalescence of close binaries. So far, there had not been any confirmed detection of a white dwarf (WD) companion to any of the BSSs in this cluster. Here, we present the detection of WD companions to 5 bright BSSs in M67. The multiwavelength spectral energy distributions covering 0.12 -11.5 $\mu$m range, were found to require binary spectral fits for 5 BSSs, consisting of a cool (BSS) and a hot companion. The parameters (Luminosity, Temperature, Radius and Mass) of the hot companions suggest them to be WDs with mass in the range 0.2 - 0.35 M$_{\odot}$ with T$_{eff}$ $\sim$ 11000 - 24000 K.

The old open cluster M67, populated with blue straggler stars (BSSs), is a well known test bed to study the BSS formation pathways. Here, we report the first direct detection of a white dwarf (WD) companion to a BSS in M67, using far-UV images from the \textit{Ultra Violet Imaging telescope} ({\it UVIT}) on {\it ASTROSAT}. Near-simultaneous observations in three far-UV bands combined with {\it GALEX}, {\it IUE}, ground and space based photometric data covering 0.14 -11.5 $\mu$m range for WOCS1007 were found to require a binary fit to its spectral energy distribution (SED), consisting of a BSS and a hot companion. On the other hand, a single spectral fit was found to be satisfactory for the SEDs of two other BSSs, WOCS1006 and WOCS2011, with the latter showing a deficient far-UV flux. The hot companion of WOCS1007 is found to have a T$_{eff}$ $\sim$ 13250-13750K and a radius of 0.09$\pm$0.01 R$_\odot$. A comparison with WD models suggests it to be a low mass WD ($\sim$ 0.18$M_\odot$), in agreement with the kinematic mass from the literature. As a low mass WD ($<$ 0.4$M_\odot$) necessitates formation through mass transfer (MT) in close binaries, WOCS1007 with a known period of 4.2 days along with its fast rotation, is likely to be formed by a case A or case B binary evolution.

In this paper, we continue our study on the evolution of black holes (BHs) that receive velocity kicks at the origin of their host star cluster potential. We now focus on BHs in rotating clusters that receive a range of kick velocities in different directions with respect to the rotation axis. We perform N-body simulations to calculate the trajectories of the kicked BHs and develop an analytic framework to study their motion as a function of the host cluster and the kick itself. Our simulations indicate that for a BH that is kicked outside of the cluster's core, as its orbit decays in a rotating cluster the BH will quickly gain angular momentum as it interacts with stars with high rotational frequencies. Once the BH decays to the point where its orbital frequency equals that of local stars, its orbit will be circular and dynamical friction becomes ineffective since local stars will have low relative velocities. After circularization, the BH's orbit decays on a longer timescale than if the host cluster was not rotating. Hence BHs in rotating clusters will have longer orbital decay times. The timescale for orbit circularization depends strongly on the cluster's rotation rate and the initial kick velocity, with kicked BHs in slowly rotating clusters being able to decay into the core before circularization occurs. The implication of the circularization phase is that the probability of a BH undergoing a tidal capture event increases, possibly aiding in the formation of binaries and high-mass BHs.

Nuclear star clusters that surround supermassive black holes (SMBHs) in galactic nuclei are thought to contain large numbers of black holes (BHs) and neutron stars (NSs), a fraction of which form binaries and could merge by Kozai-Lidov oscillations (KL). Triple compact objects are likely to be present, given what is known about the multiplicity of massive stars, whose life ends either as a NS or a BH. In this paper, we present a new possible scenario for merging BHs and NSs in galactic nuclei. We study the evolution of a triple black hole (BH) or neutron star (NS) system orbiting an SMBH in a galactic nucleus by means of direct high-precision $N$-body simulations, including Post-Newtonian terms. We find that the four-body dynamical interactions can increase the KL angle window for mergers compared to the binary case and make BH and NS binaries merge on shorter timescales. We show that the merger fraction can be up to $\sim 5$--$8$ times higher for triples than for binaries. Therefore, even if the triple fraction is only $\sim 10\%$--$20\%$ of the binary fraction, they could contribute to the merger events observed by LIGO/VIRGO in comparable numbers.

We combine sophisticated high precision scattering experiments, together with results from the Millenium-II simulation, to compute the cosmic merger rate of bound compact object (CO) binaries dynamically interacting with supermassive black hole binaries (SMBHBs). We consider binaries composed of white dwarfs (WDs), neutron stars (NSs) and black holes (BHs). The overall merger rates for WD-WD, NS-NS, BH-BH, BH-NS binaries and EBBH (eccentric binaries of black holes) from redshift $\sim 5$ are found to be $4.32\times 10^3\,{\rm yr}^{-1}$($5.93\times10^2\,{\rm yr}^{-1}$ for Type Ia SNe), $82.7\,{\rm yr}^{-1}$, $96.3\,{\rm yr}^{-1}$, $13.1\,{\rm yr}^{-1}$ and $148\,{\rm yr}^{-1}$ , respectively, for a nominal CO binary fraction in the Galactic centre of 0.1. We calculate the distance ($R$) distribution of the merger sites with respect to the host galaxies of the binaries. The distribution shows a wide range of distances up to $\sim$Mpc; this tail is produced by escaped hyper-velocity CO binaries. Due to the differences in the matter density of the surrounding environment, merger events with different $R$ are expected to display significantly different signatures in their EM counterparts. In particular, merger events (and especially NS-NS) producing a relativistic jet but occurring in the intergalactic medium will have very weak afterglow radiation relative to their prompt emission. These events, which we call 'off-center', can only be produced from a close encounter between CO binaries and SMBHBs; hence the detection of such merger events would indicate the existence of nearby SMBHBs, and in particular with high mass ratio, produced in the aftermath of a major galaxy merger.

Nuclear star clusters surrounding supermassive black holes (SMBHs) in galactic nuclei contain large numbers of stars, black holes (BHs) and neutron stars (NSs), a fraction of which are likely to form binaries. These binaries were suggested to form a triple system with the SMBH, which acts as a perturber and may enhance BH and NS mergers via the Lidov-Kozai mechanism. We follow-up previous studies, but for the first time perform an extensive statistical study of BH-BH, NS-NS and BH-NS binary mergers by means of direct high-precision regularized $N$-body simulations, including Post-Newtonian (PN) terms up to order PN2.5. We consider different SMBH masses, slopes for the BH mass function, binary semi-major axis and eccentricity distributions, and different spatial distributions for the binaries. We find that the merger rates are a decreasing function of the SMBH mass and are in the ranges $\sim 0.17$-$0.52 \ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$, $\sim 0.06$-$0.10 \ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$ and $\sim 0.04$-$0.16 \ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$ for BH-BH, BH-NS and NS-NS binaries, respectively. However, the rate estimate from this channel remains highly uncertain and depends on the specific assumptions regarding the star-formation history in galactic nuclei and the supply rate of compact objects. We find that $\sim 10\%$--$20\%$ of the mergers enter the LIGO band with eccentricities $\gtrsim 0.1$. We also compare our results to the secular approximation, and show that $N$-body simulations generally predict a larger number of mergers. Finally, these events can also be observable via their electromagnetic counterparts, thus making these compact object mergers especially valuable for cosmological and astrophysical purposes.

In this paper we consider a scenario where the currently observed hypervelocity stars in our Galaxy have been ejected from the Galactic center as a result of dynamical interactions with an intermediate-mass black hole (IMBH) orbiting the central supermassive black hole (SMBH). By performing 3-body scattering experiments, we calculate the distribution of the ejected stars' velocities given various parameters of the IMBH-SMBH binary: IMBH mass, semimajor axis and eccentricity. We also calculate the rates of change of the BH binary orbital elements due to those stellar ejections. One of our new findings is that the ejection rate depends (although mildly) on the rotation of the stellar nucleus (its total angular momentum). We also compare the ejection velocity distribution with that produced by the Hills mechanism (stellar binary disruption) and find that the latter produces faster stars on average. Also, the IMBH mechanism produces an ejection velocity distribution which is flattened towards the BH binary plane while the Hills mechanism produces a spherically symmetric one. The results of this paper will allow us in the future to model the ejection of stars by an evolving BH binary and compare both models with \textit{Gaia} observations, for a wide variety of environments (galactic nuclei, globular clusters, the Large Magellanic Clouds, etc.).

Prior to the detection of black holes (BHs) via the gravitational waves (GWs) they generate at merger, the presence of BHs was inferred in X-ray binaries, mostly via dynamical measurements, with masses in the range between $\sim 5-20~M_\odot$. The LIGO discovery of the first BHs via GWs was surprising in that the two BHs that merged had masses of $35.6^{+4.8}_{-3.0}$ and $30.6^{+3.0}_{-4.4}\,M_\odot$, which are both above the range inferred from X-ray binaries. With 20 BH detections from the O1/O2 runs, the distribution of masses remains generally higher than the X-ray inferred one, while the effective spins are generally lower, suggesting that, at least in part, the GW-detected population might be of dynamical origin rather than produced by the common evolution of field binaries. Here we perform high-resolution N-body simulations of a cluster of isolated BHs with a range of initial mass spectra and upper mass cut-offs, and study the resulting binary mass spectrum resulting from the dynamical interactions. Our clusters have properties similar to those of the massive remnants in an OB association $\sim 10 \, \mathrm{Myr}$ after formation. We perform a likelihood analysis for each of our dynamically-formed binary population against the data from the O1 and O2 LIGO/Virgo runs. We find that an initial mass spectrum $M_{\rm BH}\propto M^{-2.35}$ with an upper mass cutoff $M_{\rm max}\sim 50M_\odot$ is favored by the data, together with a slight preference for a merger rate that increases with redshift.

We report the discovery of a hot white dwarf (WD) companion to a blue straggler star (BSS) in the globular cluster (GC) NGC 5466, based on observations from the Ultra-Violet Imaging Telescope (UVIT) on board AstroSat. The Spectral Energy Distribution (SED) of the Far-UV detected BSS NH 84 was constructed by combining the flux measurements from 4 filters of UVIT, with GALEX, GAIA and other ground-based observations. The SED of NH 84 reveals the presence of a hot companion to the BSS. The temperature and radius of the BSS (T$_{\mathrm{eff}} = 8000^{+1000}_{-250}$ K, R/R$_\odot = 1.44 \pm 0.05$) derived from Gemini spectra and SED fitting using Kurucz atmospheric models are consistent with each other. The temperature and radius of the hotter companion of NH 84 (T$_{\mathrm{eff}} = 32,000 \pm 2000$ K, R/R$_\odot = 0.021 \pm 0.007$) derived by fitting Koester WD models to the SED suggest that it is likely to be a hot WD. The radial velocity derived from the spectra along with the proper motion from GAIA DR2 confirms NH 84 to be a kinematic member of the cluster. This is the second detection of a BSS-WD candidate in a GC, and the first in the outskirts of a low density GC. The location of this BSS in NGC 5466 along with its dynamical age supports the mass-transfer pathway for BSS formation in low density environments.

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.

About a century ago, Jeans (1919) discovered that if binary stars reach a state approximating energy equipartition, for example through many dynamical encounters that exchange energy, their eccentricity distribution can be described by : dN/de = 2e. This is referred to as the thermal eccentricity distribution, and has been widely used for initial conditions in theoretical investigations of binary stars. However, observations suggest that the eccentricity distributions of most observed binaries, and particularly those with masses < 5 Msun, are flatter than thermal and follow more closely to a uniform distribution. Nonetheless, it is often argued that dynamical interactions in a star cluster would quickly thermalize the binaries, which could justify imposing a thermal eccentricity distribution at birth for all binaries. In this paper we investigate the validity of this assumption. We develop our own rapid semi-analytic model for binary evolution in star clusters, and also compare with detailed N-body and Monte Carlo star cluster models. We show that, for nearly all binaries, dynamical encounters fail to convert an initially uniform eccentricity distribution to thermal within a star cluster's lifetime. Thus, if a thermal eccentricity distribution is observed, it is likely imprinted upon formation rather than through subsequent long-term dynamical processing. Theoretical investigations that initialize all binaries with a thermal distribution will make incorrect predictions for the evolution of the binary population. Such models may overpredict the merger rate for binaries with modest orbital separations by a factor of about two.

We propose a formation mechanism for twin blue stragglers (BSs) in compact binaries that involves mass transfer from an evolved outer tertiary companion on to the inner binary via a circumbinary disk. We apply this scenario to the observed double BS system Binary 7782 in the old open cluster NGC 188, and show that its observed properties are naturally reproduced within the context of the proposed model. We predict the following properties for twin BSs: (1) For the outer tertiary orbit, the initial orbital period should lie between 220 days $\lesssim$ P$_{\rm out}$ $\lesssim$ 1100 days, assuming initial masses for the inner binary components of $m_{\rm 1} = 1.1$ M$_{\odot}$ and $m_{\rm 2} =$ 0.9 M$_{\odot}$ and an outer tertiary mass of $m_{\rm 3} = 1.4$ M$_{\odot}$. After Roche-lobe overflow, the outer star turns into a white dwarf (WD) of mass 0.43 to 0.54\,\MSun. There is a correlation between the mass of this WD and the outer orbital period: more massive WDs will be on wider orbits. (3) The rotational axes of both BSs will be aligned with each other and the orbital plane of the outer tertiary WD. (4) The BSs will have roughly equal masses, independent of their initial masses (since the lower mass star accretes the most). The dominant accretor should, therefore, be more enriched by the accreted material. Hence, one of the BSs will appear to be more enriched by either He, C and O or by s-process elements, if the donor started Roche lobe overflow on, respectively, the red giant or asymptotic giant branch. (5) Relative to old clusters, twin BSs in close binaries formed from the proposed mechanism should be more frequent in the Galactic field and open clusters with ages $\lesssim$ 4-6 Gyr, since then the donor will have a radiative envelope. (6) The orbit of the binary BS will have a small semi-major axis (typically $\aplt 0.3$\,au) and be close to circular ($e \aplt 0.2$).

We present a semi-analytic model for self-consistently evolving a population of globular clusters (GCs) in a given host galaxy across cosmic time. We compute the fraction of GCs still hosting intermediate-mass black holes (IMBHs) at a given redshift in early and late type galaxies of different masses and sizes, and the corresponding rate of tidal disruption events (TDEs), both main-sequence (MS) and white dwarf (WD) stars. We find that the integrated TDE rate for the entire GC population can exceed the corresponding rate in a given galactic nucleus and that $\sim 90$% of the TDEs reside in GCs within a maximum radius of $\sim 2-15$ kpc from the host galaxy's center. This suggests that observational efforts designed to identify TDEs should not confine themselves to galactic nuclei alone, but should also consider the outer galactic halo where massive old GCs hosting IMBHs would reside. Indeed, such off-centre TDEs as predicted here may already have been observed. MS TDE rates are more common than WD TDE rates by a factor 30 (100) at $z\leq 0.5$ ($z=2$). We also calculate the rate of IMBH-SBH mergers across cosmic time, finding that the typical IMRI rate at low redshift is of the order of $\sim 0.5-3$ Gpc$^{-3}$ yr$^{-1}$, which becomes as high as $\sim 100$ Gpc$^{-3}$ yr$^{-1}$ near the peak of GC formation. Advanced LIGO combined with VIRGO, KAGRA, ET and LISA will be able to observe the bottom-end and top-end of the IMBH population, respectively.

For a sample of 38 Galactic globular clusters (GCs), we confront the observed distributions of blue straggler (BS) proper motions and masses (derived from isochrone fitting) from the BS catalog of Simunovic & Puzia with theoretical predictions for each of the two main competing BS formation mechanisms. These are mass transfer from an evolved donor on to a main-sequence (MS) star in a close binary system, and direct collisions involving MS stars during binary encounters. We use the \texttt{FEWBODY} code to perform simulations of single-binary and binary-binary interactions. This provides collisional velocity and mass distributions for comparison to the observed distributions. Most clusters are consistent with BSs derived from a dynamically relaxed population, supportive of the binary mass-transfer scenario. In a few clusters, including all the post-core collapse clusters in our sample, the collisional velocities provide the best fit.

Hypervelocity binaries (HVBs) can be produced by a close encounter between stellar binaries and supermassive black hole (SMBH)-intermediate massive black hole (IMBH) binaries. We perform scattering experiments between tight stellar binaries and an SMBH-IMBH binary (within the observationally allowed parameters), using high-precision N-body integrations. Simulation results show that only a few HVBs can be produced in the Galactic center if it harbors a $4\times10^3M_\odot$ IMBH. However, we find that tight HVBs can be efficiently produced in prograde coplanar orbits with respect to the SMBH-IMBH binary. Therefore, the trajectories of the HVBs provide a promising way to constrain the IMBH orbital plane, if one exists. More importantly, we find that the ejected tight binaries are eccentric. For other galaxies hosting equal-mass SMBH binaries for which the HVB ejection rate is much higher relative to SMBH-IMBH binaries, these ejected compact binaries constitute promising extragalactic merger sources that could be detected by gravitational wave detectors.

We present a drag force model for evolving chaotic few-body interactions with the inclusion of orbital energy losses, such as tidal dissipation and gravitational wave (GW) emission. The main effect from such losses is the formation of two-body captures, that for compact objects result in GW mergers, and for stars lead to either compact binaries, mergers or disruptions. Studying the inclusion of energy loss terms in few-body interactions is therefore likely to be important for modeling and understanding the variety of transients that soon will be observed by current and upcoming surveys. However, including especially tides in few-body codes has been shown to be technically difficult and computationally heavy, which has lead to very few systematic tidal studies. In this paper we derive a drag force term that can be used to model the effects from tidal, as well as other, energy losses in few-body interactions, if the two-body orbit averaged energy loss is known a priori. This drag force model is very fast to evolve, and gives results in agreement with other approaches, including the impulsive and affine tide approximations.

Stars passing too close to a super massive black hole (SMBH) can produce tidal disruption events (TDEs). Since the resulting stellar debris can produce an electromagnetic flare, TDEs are believed to probe the presence of single SMBHs in galactic nuclei, which otherwise remain dark. In this paper, we show how stars orbiting an IMBH secondary are perturbed by an SMBH primary. We find that the evolution of the stellar orbits are severely affected by the primary SMBH due to secular effects and stars orbiting with high inclinations with respect to the SMBH-IMBH orbital plane end their lives as TDEs due to Kozai-Lidov oscillations, hence illuminating the secondary SMBH/IMBH. Above a critical SMBH mass of $\approx 1.15 \times 10^8$ M$_{\odot}$, no TDE can occur for typical stars in an old stellar population since the Schwarzschild radius exceeds the tidal disruption radius. Consequently, any TDEs due to such massive SMBHs will remain dark. It follows that no TDEs should be observed in galaxies with bulges more massive than $\approx 4.15\times 10^{10}$ M$_{\odot}$, unless a lower-mass secondary SMBH or IMBH is also present. The secular mechanism for producing TDEs considered here therefore offers a useful probe of SMBH-SMBH/IMBH binarity in the most massive galaxies. We further show that the TDE rate can be $\approx 10^{-4}-10^{-3}$ yr$^{-1}$, and that most TDEs occur on $\approx 0.5$ Myr. Finally, we show that stars may be ejected with velocities up to thousands of km s$^{-1}$, which could contribute to the observed population of Galactic hypervelocity stars.

In this paper, we study the secular dynamical evolution of binaries composed of intermediate-mass and stellar-mass black holes (IMBHs and SBHs, respectively) in orbit about a central super-massive black hole (SMBH) in galactic nuclei. Such BH triplets could form via the inspiral of globular clusters toward galactic nuclei due to dynamical friction, or even major/minor galaxy mergers. We perform, for reasonable initial conditions that we justify, sophisticated $N$-body simulations that include both regularization and Post-Newtonian corrections. We find that mass segregation combined with Kozai-Lidov oscillations induced by the primary SMBH can effectively merge IMBH-SBH binaries on time-scales much shorter than gravitational wave emission alone. Moreover, the rate of such extreme mass ratio inspirals could be high ($\sim 1\ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$) in the local Universe, but these are expected to be associated with recent GC infall or major/minor mergers, making the observational signatures of such events (e.g., tidal debris) good diagnostics for searching for SMBH-IMBH-SBH mergers. A small fraction could also be associated with tidal disruption events by the IMBH-SBH during inspiral.

The Kepler mission has detected a number of transiting circumbinary planets (CBPs). Although currently not detected, exomoons could be orbiting some of these CBPs, and they might be suitable for harboring life. A necessary condition for the existence of such exomoons is their long-term dynamical stability. Here, we investigate the stability of exomoons around the Kepler CBPs using numerical $N$-body integrations. We determine regions of stability and obtain stability maps in the (a_m,i_pm) plane, where a_m is the initial exolunar semimajor axis with respect to the CBP, and i_pm is the initial inclination of the orbit of the exomoon around the planet with respect to the orbit of the planet around the stellar binary. Ignoring any dependence on i_pm, for most Kepler CBPs the stability regions are well described by the location of the 1:1 mean motion commensurability of the binary orbit with the orbit of the moon around the CBP. This is related to a destabilizing effect of the binary compared to the case if the binary were replaced by a single body, and which is borne out by corresponding 3-body integrations. For high inclinations, the evolution is dominated by Lidov-Kozai oscillations, which can bring moons in dynamically stable orbits to close proximity within the CBP, triggering strong interactions such as tidal evolution, tidal disruption, or direct collisions. This suggests that there is a dearth of highly-inclined exomoons around the Kepler CBPs, whereas coplanar exomoons are dynamically allowed.

In this paper, the fourth in the series, we continue our study of combinatorics in chaotic Newtonian dynamics. We focus once again on the chaotic four-body problem in Newtonian gravity assuming finite-sized particles, and interactions that produce direct collisions between any two particles. Our long-term goal is to predict the probability of a given collision event occurring over the course of an interaction, as a function of the numbers and properties of the particles. In previous papers, we varied the number of interacting particles, as well as the distributions of particle radii and masses. Here, we refine the methods developed in these preceding studies, and arrive at a final and robust methodology that can be used to study collisional dynamics in a variety of astrophysical contexts, ranging from stars in star clusters, galaxies in galaxy groups and clusters and even the collisional growth of planetesimals in protoplanetary disks. We further present and refine the concept of a Collision Rate Diagram (CRD), the primary tool we use to quantify the relative rates for different collision scenarios to occur. The agreement between our final theoretical model and the results of numerical scattering simulations is excellent.

We present a half-life formalism for describing the disruption of gravitationally-bound few-body systems, with a focus on binary-binary scattering. For negative total encounter energies, the four-body problem has three possible decay products in the point particle limit. For each decay product and a given set of initial conditions, we obtain directly from numerical scattering simulations the half-life for the distribution of disruption times. As in radioactive decay, the half-lives should provide a direct prediction for the relative fractions of each decay product. We test this prediction with simulated data and find good agreement with our hypothesis. We briefly discuss applications of this feature of the gravitational four-body problem to populations of black holes in globular clusters. This paper, the second in the series, builds on extending the remarkable similarity between gravitational chaos at the macroscopic scale and radioactive decay at the microscopic scale to larger-N systems.

Black hole mergers detectable with LIGO can occur in active galactic nucleus (AGN) disks. Here we parameterize the merger rates, the mass spectrum and the spin spectrum of black holes (BH) in AGN disks. The predicted merger rate spans $\sim 10^{-4}-10^{4} \rm{Gpc}^{-1} \rm{yr}^{-1}$, so upper limits from LIGO ($<212\rm{Gpc}^{-1}\rm{yr}^{-1}$) already constrain it. The predicted mass spectrum has the form of a broken power-law consisting of a pre-existing BH powerlaw mass spectrum and a harder powerlaw mass spectrum resulting from mergers. The predicted spin spectrum is multi-peaked with the evolution of retrograde spin BH in the gas disk playing a key role. We outline the large uncertainties in each of these LIGO observables for this channel and we discuss ways in which they can be constrained in the future.

We present a formalism for constructing schematic diagrams to depict chaotic three-body interactions in Newtonian gravity. This is done by decomposing each interaction in to a series of discrete transformations in energy- and angular momentum-space. Each time a transformation is applied, the system changes state as the particles re-distribute their energy and angular momenta. These diagrams have the virtue of containing all of the quantitative information needed to fully characterize most bound or unbound interactions through time and space, including the total duration of the interaction, the initial and final stable states in addition to every intervening temporary meta-stable state. As shown via an illustrative example for the bound case, prolonged excursions of one of the particles, which by far dominates the computational cost of the simulations, are reduced to a single discrete transformation in energy- and angular momentum-space, thereby potentially mitigating any computational expense. We further generalize our formalism to sequences of (unbound) three-body interactions, as occur in dense stellar environments during binary hardening. Finally, we provide a method for dynamically evolving entire populations of binaries via three-body scattering interactions, using a purely analytic formalism. In principle, the techniques presented here are adaptable to other three-body problems that conserve energy and angular momentum.

The evolution of main sequence binaries that reside in the Galactic Centre can be heavily influenced by the central supermassive black hole (SMBH). Due to these perturbative effects, the stellar binaries in dense environments are likely to experience mergers, collisions or ejections through secular and/or non-secular interactions. More direct interactions with the central SMBH are thought to produce hypervelocity stars (HVSs) and tidal disruption events (TDEs). In this paper, we use N-body simulations to study the dynamics of stellar binaries orbiting a central SMBH primary with an outer SMBH secondary orbiting this inner triple. The effects of the secondary SMBH on the event rates of HVSs, TDEs and stellar mergers are investigated, as a function of the SMBH-SMBH binary mass ratio. Our numerical experiments reveal that, relative to the isolated SMBH case, the TDE and HVS rates are enhanced for, respectively, the smallest and largest mass ratio SMBH-SMBH binaries. This suggests that the observed event rates of TDEs and HVSs have the potential to serve as a diagnostic of the mass ratio of a central SMBH-SMBH binary. The presence of a secondary SMBH also allows for the creation of hypervelocity binaries. Observations of these systems could thus constrain the presence of a secondary SMBH in the Galactic Centre.

We assess the contribution of dynamical hardening by direct three-body scattering interactions to the rate of stellar-mass black hole binary (BHB) mergers in galactic nuclei. We derive an analytic model for the single-binary encounter rate in a nucleus with spherical and disk components hosting a super-massive black hole (SMBH). We determine the total number of encounters $N_{\rm GW}$ needed to harden a BHB to the point that inspiral due to gravitational wave emission occurs before the next three-body scattering event. This is done independently for both the spherical and disk components. Using a Monte Carlo approach, we refine our calculations for $N_{\rm GW}$ to include gravitational wave emission between scattering events. For astrophysically plausible models we find that typically $N_{\rm GW} \lesssim$ 10. We find two separate regimes for the efficient dynamical hardening of BHBs: (1) spherical star clusters with high central densities, low velocity dispersions and no significant Keplerian component; and (2) migration traps in disks around SMBHs lacking any significant spherical stellar component in the vicinity of the migration trap, which is expected due to effective orbital inclination reduction of any spherical population by the disk. We also find a weak correlation between the ratio of the second-order velocity moment to velocity dispersion in galactic nuclei and the rate of BHB mergers, where this ratio is a proxy for the ratio between the rotation- and dispersion-supported components. Because disks enforce planar interactions that are efficient in hardening BHBs, particularly in migration traps, they have high merger rates that can contribute significantly to the rate of BHB mergers detected by the advanced Laser Interferometer Gravitational-Wave Observatory.

We consider how dynamical friction acts on black holes that receive a velocity kick while located at the center of a gravitational potential, analogous to a star cluster, due to either a natal kick or the anisotropic emission of gravitational waves during a black hole-black hole merger. Our investigation specifically focuses on how well various Chandrasekhar-based dynamical friction models can predict the orbital decay of kicked black holes with $m_{bh} \lesssim 100 M_\odot$ due to an inhomogeneous background stellar field. In general, the orbital evolution of a kicked black hole follows that of a damped oscillator where two-body encounters and dynamical friction serve as sources of damping. However, we find models for approximating the effects of dynamical friction do not accurately predict the amount of energy lost by the black hole if the initial kick velocity $v_{k}$ is greater than the stellar velocity dispersion $\sigma$. For all kick velocities, we also find that two-body encounters with nearby stars can cause the energy evolution of a kicked BH to stray significantly from standard dynamical friction theory as encounters can sometimes lead to an energy gain. For larger kick velocities, we find the orbital decay of a black hole departs from classical theory completely as the black hole's orbital amplitude decays linearly with time as opposed to exponentially. Therefore, we have developed a linear decay formalism which scales linearly with black hole mass and $\frac{v_{k}}{\sigma}$ in order to account for the variations in the local gravitational potential.

In this paper, the third in the series, we continue our study of combinatorics in chaotic Newtonian dynamics. We study the chaotic four-body problem in Newtonian gravity assuming finite-sized particles, and we focus on interactions that produce direct collisions between any two stars. Our long-term goal is to construct an equation that gives the probability of a given collision event occurring over the course of the interaction, as a function of the total encounter energy and angular momentum as well as the numbers and properties of the particles. In previous papers, we varied the number of interacting particles and the distribution of particle radii, for all equal mass particles. Here, we focus on the effects of different combinations of particle masses. We develop an analytic formalism for calculating the time-scales for different collision scenarios to occur. Our analytic time-scales reproduce the simulated time-scales when gravitational focusing is included. We present a method for calculating the relative rates for different types of collisions to occur, assuming two different limits for the particle orbits; radial and tangential. These limits yield relative collision probabilities that bracket the probabilities we obtain directly from numerical scattering experiments, and are designed to reveal important information about the (time-averaged) trajectories of the particles as a function of the interaction parameters. Finally, we present a Collision Rate Diagram (CRD), which directly compares the predictions of our analytic rates to the simulations and quantifies the quality of the agreement. The CRD will facilitate refining our analytic collision rates in future work, as we expand in to the remaining parameter space.

The mass evaporation rate of globular clusters evolving in a strong Galactic tidal field is derived through the analysis of large, multi-mass $N$-body simulations. For comparison, we also study the same evaporation rates using MOCCA Monte Carlo models for globular cluster evolution. Our results show that the mass evaporation rate is a dynamical value; that is, far from a constant single number found in earlier analytical work and commonly used in the literature. Moreover, the evaporation rate derived with these simulations is higher than values previously published. These models also show that the value of the mass evaporation rate depends on the strength of the tidal field. We give an analytical estimate of the mass evaporation rate as a function of time and galactocentric distance xi(R_{GC},t). Upon extrapolating this formula to smaller R_{GC} values, our results provide tentative evidence for a very high xi value at small R_{GC}. Our results suggest that the corresponding mass loss in the inner Galactic potential could be high and it should be accounted for when star clusters pass within it. This has direct relevance to nuclear cluster formation/growth via the infall of globular clusters through dynamical friction. As an illustrative example, we estimate how the evaporation rate increases for a ~10^5 solar masses globular cluster that decays through dynamical friction into the galactic centre. We discuss the findings of this work in relation to the formation of nuclear star clusters by inspiraling globular clusters.

We study the formation of runaway stars due to binary-binary (2+2) interactions in young star-forming clusters and/or associations. This is done using a combination of analytic methods and numerical simulations of 2+2 scattering interactions, both in isolation and in a homogeneous background potential. We focus on interactions that produce two single stars and a binary, and study the outcomes as a function of the depth of the background potential, within a range typical of cluster cores. As reference parameters for the observational properties, we use those observed for the system of runaway stars AE Aur and $\mu$ Col and binary $\iota$ Ori. We find that the outcome fractions have no appreciable dependence on the depth of the potential, and neither do the velocities of the ejected single stars. However, as the potential gets deeper and a larger fraction of binaries remain trapped, two binary populations emerge, with the escaped component having higher speeds and shorter semi-major axes than the trapped one. Additionally, we find that the relative angles between the ejected products are generally large. In particular, the angle between the ejected fastest star and the escaped binary is typically $\gtrsim 120-135^{\circ}$, with a peak at around $160^{\circ}$. However, as the potential gets deeper, the angle distribution becomes broader. Finally, we discuss the implications of our results for the interpretation of the properties of the runaway stars AE Aur and $\mu$ Col.

Sub-subgiants are a new class of stars that are optically redder than normal main-sequence stars and fainter than normal subgiant stars. Sub-subgiants, and the possibly related red stragglers (which fall to the red of the giant branch), occupy a region of the color-magnitude diagram that is predicted to be devoid of stars by standard stellar evolution theory. In previous papers we presented the observed demographics of these sources and defined possible theoretical formation channels through isolated binary evolution, the rapid stripping of a subgiant's envelope, and stellar collisions. Sub-subgiants offer key tests for single- and binary-star evolution and stellar collision models. In this paper, we synthesize these findings to discuss the formation frequencies through each of these channels. The empirical data, our analytic formation rate calculations, and analyses of sub-subgiants in a large grid of Monte Carlo globular cluster models suggest that the binary evolution channels may be the most prevalent, though all channels appear to be viable routes to sub-subgiant creation (especially in higher-mass globular clusters). Multiple formation channels may operate simultaneously to produce the observed sub-subgiant population. Finally, many of these formation pathways can produce stars in both the sub-subgiant and red straggler (and blue straggler) regions of the color-magnitude diagram, in some cases as different stages along the same evolutionary sequence.

In this paper, we study the 3-body products (two single stars and a binary) of binary-binary (2+2) scattering interactions. This is done using a combination of analytic methods and numerical simulations of 2+2 scattering interactions, both in isolation and in a homogeneous background potential. We derive analytically a simple formula relating the angle between the velocity vectors of the two ejected single stars and the orbital separation of the remaining binary. We compare our analytic formulation to numerical scattering simulations, and illustrate that the agreement is excellent, both in isolation and in a homogeneous background potential. Our results are ideally suited for application to the GAIA database, which is expected to identify many hundred runaway stars. The analytic relation presented here has the potential to identify runaway stars formed dynamically with high confidence. Finally, by applying our method to the runaways AE Aur and $\mu$ Col, we illustrate that it can be used to constrain the history of the background potential, which was denser than the presently observed density in the case of the Trapezium cluster.

Sub-subgiants are stars observed to be redder than normal main-sequence stars and fainter than normal subgiant (and giant) stars in an optical color-magnitude diagram. The red straggler stars, which lie redward of the red giant branch, may be related and are often grouped together with the sub-subgiants in the literature. These stars defy our standard theory of single-star evolution, and are important tests for binary evolution and stellar collision models. In total, we identify 65 sub-subgiants and red stragglers in 16 open and globular star clusters from the literature; 50 of these, including 43 sub-subgiants, pass our strict membership selection criteria (though the remaining sources may also be cluster members). In addition to their unique location on the color-magnitude diagram, we find that at least 58% (25/43) of sub-subgiants in this sample are X-ray sources with typical 0.5-2.5 keV luminosities of order 10^30 - 10^31 erg/s. Their X-ray luminosities and optical-to-X-ray flux ratios are similar to those of RS CVn active binaries. At least 65% (28/43) of the sub-subgiants in our sample are variables, 21 of which are known to be radial-velocity binaries. Typical variability periods are <15 days. At least 33% (14/43) of the sub-subgiants are H-alpha emitters. These observational demographics provide strong evidence that binarity is important for sub-subgiant formation. Finally, we find that the number of sub-subgiants per unit mass increases toward lower-mass clusters, such that the open clusters in our sample have the highest specific frequencies of sub-subgiants.

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.

We perform a large suite of $N=4$ numerical scattering experiments between two identical binaries consisting of identical point particles in a (continuous) background potential. For investigative purposes, albeit without loss of generality, we assume that the potential corresponds to a uniform (natal or star-forming) gas medium. We explore a range of constant gas densities, from $n=10~ {\rm cm}^{-3}$ to $10^{5}~ {\rm cm}^{-3}$. These densities are relevant for various astrophysical environments, including molecular clouds (i.e., star-forming regions) and denser, fragmented cores within these clouds. Our primary goal is to characterize the effects of the background potential on the subsequent stellar dynamics. We consider the outcome probabilities as well as the properties of any binaries formed during the binary-binary encounters, such as the distributions of binary binding energies and eccentricities. We also present the final velocity distributions of the ejected single stars. The background potential has two important effects on the stellar dynamics: 1) The potential acts to reset the zero-point of the total system energy, which in turn affects the types and properties of the products of the encounter; 2) For higher $n$ and weakly bound systems (i.e., large semimajor axes), the stellar dynamics are significantly affected when stars become trapped in the potential, oscillating around the system centre of mass (CM). This, in turn, increases the number of scattering events between stars (single, binary or triple) near the CM and makes it harder for single stars to escape to infinity. This ultimately leads to the preferential ionization of triples and wide binaries and the survival of compact binaries, with the single stars escaping at very high ejection velocities.

We consider a novel mechanism to account for the observed distance-redshift relation. This is done by presenting a toy model for the large-scale matter distribution in a static Universe. Our model mainly concerns particles with masses far below those in the Standard Model of Particle Physics. The model is founded on three main assumptions: (1) a mass spectrum dN$_{\rm i}$/dm$_{\rm i}$ $=$ $\beta$m$_{\rm i}^{-\alpha}$ (where $\alpha$ and $\beta$ are both positive constants) for low-mass particles with m$_{\rm i}$ $\ll$ 10$^{-22}$ eV $\ll$ M$_{\rm P}$, where M$_{\rm P}$ is the Planck mass; (2) a particle mass-wavelength relation of the form $\lambda_{\rm i} =$ $\hbar$/$\delta_{\rm i}$m$_{\rm i}$c, where $\delta_{\rm i} =$ $\eta$m$_{\rm i}^{\gamma}$ and $\eta$ and $\gamma$ are both constants; and (3) For such low-mass particles, locality can only be defined on large spatial scales, comparable to or exceeding the particle wavelengths. We use our model to derive the cosmological redshift characteristic of the Standard Model of Cosmology, which becomes a gravitational redshift in our model. We compare the results of our model to empirical data and show that, in order to reproduce the sub-linear form of the observed distance-redshift relation, our model requires $\alpha$ + $\gamma$ $<$ 0. We further place our toy model in the context of the Friedmann Universe via a superposition of Einstein Universes, each with its own scale factor a$_{\rm i}$. Given the overwhelming evidence supporting an expanding Universe, we then address possible modifications to our base model that would be required to account for the available empirical constraints, including the addition of some initial expansion. Finally, we consider potentially observable distinctions between the cosmological redshift and our proposed mechanism to account for the observed distance-redshift relation.

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.

In this paper, we study the chaotic four-body problem in Newtonian gravity. Assuming point particles and total encounter energies $\le$ 0, the problem has three possible outcomes. We describe each outcome as a series of discrete transformations in energy space, using the diagrams first presented in Leigh \& Geller (2012; see the Appendix). Furthermore, we develop a formalism for calculating probabilities for these outcomes to occur, expressed using the density of escape configurations per unit energy, and based on the Monaghan description originally developed for the three-body problem. We compare this analytic formalism to results from a series of binary-binary encounters with identical point particles, simulated using the \texttt{FEWBODY} code. Each of our three encounter outcomes produces a unique velocity distribution for the escaping star(s). Thus, these distributions can potentially be used to constrain the origins of dynamically-formed populations, via a direct comparison between the predicted and observed velocity distributions. Finally, we show that, for encounters that form stable triples, the simulated single star escape velocity distributions are the same as for the three-body problem. This is also the case for the other two encounter outcomes, but only at low virial ratios. This suggests that single and binary stars processed via single-binary and binary-binary encounters in dense star clusters should have a unique velocity distribution relative to the underlying Maxwellian distribution (provided the relaxation time is sufficiently long), which can be calculated analytically.

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.

We consider the origins of enigmatic stellar populations in four Local Group galactic nuclei, specifically the Milky Way, M31, M32 and M33. These are centrally concentrated blue stars, found in three out of the four nuclear star clusters (NSCs) considered here. Their origins are unknown, but could include blue straggler (BS) stars, extended horizontal branch stars and young recently formed stars. Here, we calculate order-of-magnitude estimates for various collision rates, as a function of the host NSC environment and distance from the cluster centre. These rates are sufficiently high that BSs, formed via collisions between main sequence (MS) stars, could contribute non-negligibly ($\sim$ 1-10% in mass) to every surface brightness profile, with the exception of the Milky Way. Stellar evolution models show that the envelopes of red giant branch (RGB) stars must be nearly completely stripped to significantly affect their photometric appearance, which requires multiple collisions. Hence, the collision rates for individual RGB stars are only sufficiently high in the inner $\lesssim$ 0.1 pc of M31 and M32 for RGB destruction to occur. Collisions between white dwarfs and MS stars, which should ablate the stars, could offer a steady and significant supply of gas in every NSC in our sample. The gas could either fragment to form new stars, or accrete onto old MS stars already present. Thus, collisional processes could contribute significantly to the observed blue excesses in M31 and M33; future studies should be aimed at better constraining theoretical predictions to compliment existing and future observational data.

We discuss a new scenario for the formation of intermediate mass black holes in dense star clusters. In this scenario, intermediate mass black holes are formed as a result of dynamical interactions of hard binaries containing a stellar mass black hole, with other stars and binaries. We discuss the necessary conditions to initiate the process of intermediate mass black hole formation and the influence of an intermediate mass black hole on the host global globular cluster properties. We discuss two scenarios for intermediate mass black hole formation. The SLOW and FAST scenarios. They occur later or earlier in the cluster evolution and require smaller or extremely large central densities, respectively. In our simulations, the formation of intermediate mass black holes is highly stochastic. In general, higher formation probabilities follow from larger cluster concentrations (i.e. central densities). We further discuss possible observational signatures of the presence of intermediate mass black holes in globular clusters that follow from our simulations. These include the spatial and kinematic structure of the host cluster, possible radio, X-ray and gravitational wave emissions due to dynamical collisions or mass-transfer and the creation of hypervelocity main sequence escapers during strong dynamical interactions between binaries and an intermediate mass black hole. All simulations discussed in this paper were performed with the MOCCA Monte Carlo code. MOCCA accurately follows most of the important physical processes that occur during the dynamical evolution of star clusters but, as with other dynamical codes, it approximates the dissipative processes connected with stellar collisions and binary mergers.

We very briefly discuss proposed in the literature possible scenarios for intermediate mass black holes (IMBH) formation in globular clusters. We also discuss the results of the MOCCA simulations of about 2000 models (BigSurvey) regarding the distribution of events connected with electromagnetic and gravitational radiations, namely: mass transfer on IMBH, collisions and mergers with IMBH and mergers with IMBH due to gravitational radiation. The rates of these events are very small, so their observation is very improbable.

What is the difference between a long-lived unstable (or quasi-stable) multiple star system and a bona fide star cluster? In this paper, we present a possible framework to address this question, by studying the distributions of disruption times for chaotic gravitational encounters as a function of the number of interacting particles. To this end, we perform a series of numerical scattering experiments with the \texttt{FEWBODY} code, to calculate the distributions of disruption times as a function of both the particle number N and the virial coefficient k. The subsequent distributions are fit with a physically-motivated function, consisting of an initial exponential decay followed by a very slowly decreasing tail at long encounter times due to long-lived quasi-stable encounters. We find three primary features characteristic of the calculated distributions of disruption times. These are: (1) the system half-life increases with increasing particle number, (2) the fraction of long-lived quasi-stable encounters increases with increasing particle number and (3) both the system half-life and the fraction of quasi-stable encounters increase with decreasing virial coefficient. We discuss the significance of our results for collisional dynamics, and consider the extrapolation of our results to larger-N systems. We suggest that this could potentially offer a clear and unambiguous distinction between star clusters and (unstable or quasi-stable) multiple star systems. Although we are limited by very small-number statistics, our results tentatively suggest that (for our assumptions) this transition occurs at a critical particle number of order 100.

Galactic nuclei typically host either a Nuclear Star Cluster (NSC, prevalent in galaxies with masses $\lesssim 10^{10}M_\odot$) or a Massive Black Hole (MBH, common in galaxies with masses $\gtrsim 10^{12}M_\odot$). In the intermediate mass range, some nuclei host both a NSC and a MBH. In this paper, we explore scaling relations between NSC mass (${\cal M}_{\rm NSC}$) and host galaxy total stellar mass (${\cal M}_{\star,\rm gal}$) using a large sample of NSCs in late- and early-type galaxies, including a number of NSCs harboring a MBH. Such scaling relations reflect the underlying physical mechanisms driving the formation and (co)evolution of these central massive objects. We find $\sim\!1.5\sigma$ significant differences between NSCs in late- and early-type galaxies in the slopes and offsets of the relations $r_{\rm eff,NSC}$--${\cal M}_{\rm NSC}$, $r_{\rm eff, NSC}$--${\cal M}_{\star,\rm gal}$ and ${\cal M}_{\rm NSC}$--${\cal M}_{\star,\rm gal}$, in the sense that $i)$ NSCs in late-types are more compact at fixed ${\cal M}_{\rm NSC}$ and ${\cal M}_{\star,\rm gal}$; and $ii)$ the ${\cal M}_{\rm NSC}$--${\cal M}_{\star,\rm gal}$ relation is shallower for NSCs in late-types than in early-types, similar to the ${\cal M}_{\rm BH}$--${\cal M}_{\star,\rm bulge}$ relation. We discuss these results in the context of the (possibly ongoing) evolution of NSCs, depending on host galaxy type. For NSCs with a MBH, we illustrate the possible influence of a MBH on its host NSC, by considering the ratio between the radius of the MBH sphere of influence and $r_{\rm eff, NSC}$. NSCs harbouring a sufficiently massive black hole are likely to exhibit surface brightness profile deviating from a typical King profile.

Binary mass transfer is at the forefront of some of the most exciting puzzles of modern astrophysics, including Type Ia supernovae, gamma-ray bursts, and the formation of most observed exotic stellar populations. Typically, the evolution is assumed to proceed in isolation, even in dense stellar environments such as star clusters. In this paper, we test the validity of this assumption via the analysis of a large grid of binary evolution models simulated with the SeBa code. For every binary, we calculate analytically the mean time until another single or binary star comes within the mean separation of the mass-transferring binary, and compare this time-scale to the mean time for stable mass transfer to occur. We then derive the probability for each respective binary to experience a direct dynamical interruption. The resulting probability distribution can be integrated to give an estimate for the fraction of binaries undergoing mass transfer that are expected to be disrupted as a function of the host cluster properties. We find that for lower-mass clusters ($\lesssim 10^4$ M$_{\odot}$), on the order of a few to a few tens of percent of binaries undergoing mass-transfer are expected to be interrupted by an interloping single, or more often binary, star, over the course of the cluster lifetime, whereas in more massive globular clusters we expect $\ll$ 1% to be interrupted. Furthermore, using numerical scattering experiments performed with the FEWBODY code, we show that the probability of interruption increases if perturbative fly-bys are considered as well, by a factor $\sim 2$.

We present results of a Hubble Space Telescope far-ultraviolet (FUV) survey searching for white dwarf (WD) companions to blue straggler stars (BSSs) in open cluster NGC 188. The majority of NGC 188 BSSs (15 of 21) are single-lined binaries with properties suggestive of mass-transfer formation via Roche lobe overflow, specifically through an asymptotic giant branch star transferring mass to a main sequence secondary, yielding a BSS binary with a WD companion. In NGC 188, a BSS formed by this mechanism within the past 400 Myr will have a WD companion hot and luminous enough to be directly detected as a FUV photometric excess with HST. Comparing expected BSS FUV emission to observed photometry reveals four BSSs with WD companions above 12,000 K (younger than 250 Myr) and three WD companions with temperatures between 11,000-12,000 K. These BSS+WD binaries all formed through recent mass transfer. The location of the young BSSs in an optical color-magnitude diagram (CMD) indicates that distance from the zero-age main sequence does not necessarily correlate with BSS age. There is no clear CMD separation between mass transfer-formed BSSs and those likely formed through other mechanisms, such as collisions. The seven detected WD companions place a lower limit on the mass-transfer formation frequency of 33%. We consider other possible formation mechanisms by comparing properties of the BSS population to theoretical predictions. We conclude that 14 BSS binaries likely formed from mass transfer, resulting in an inferred mass-transfer formation frequency of approximately 67%.

We use N-body simulations to model the 12 Gyr evolution of a suite of star clusters with identical initial stellar mass functions over a range of initial cluster masses, sizes, and orbits. Our models reproduce the distribution of present-day global stellar mass functions that is observed in the Milky Way globular cluster population. We find that the slope of a star cluster's stellar mass function is strongly correlated with the fraction of mass that the cluster has lost, independent of the cluster's initial mass, and nearly independent of its orbit and initial size. Thus, the mass function - initial mass relation can be used to determine a Galactic cluster's initial total stellar mass, if the initial stellar mass function is known. We apply the mass function - initial mass relation presented here to determine the initial stellar masses of 33 Galactic globular clusters, assuming an universal Kroupa initial mass function. Our study suggests that globular clusters had initial masses that were on average a factor of 4.5 times larger than their present day mass, with three clusters showing evidence for being 10 times more massive at birth.

Strong encounters between single stars and binaries play a pivotal role in the evolution of star clusters. Such encounters can also dramatically modify the orbital parameters of binaries, exchange partners in and out of binaries, and are a primary contributor to the rate of physical stellar collisions in star clusters. Often, these encounters are studied under the approximation that they happen quickly enough and within a small enough volume to be considered isolated from the rest of the cluster. In this paper, we study the validity of this assumption through the analysis of a large grid of single - binary and binary - binary scattering experiments. For each encounter we evaluate the encounter duration, and compare this with the expected time until another single or binary star will join the encounter. We find that for lower-mass clusters, similar to typical open clusters in our Galaxy, the percent of encounters that will be "interrupted" by an interloping star or binary may be 20-40% (or higher) in the core, though for typical globular clusters we expect <1% of encounters to be interrupted. Thus, the assumption that strong encounters occur in relative isolation breaks down for certain clusters. Instead, many strong encounters develop into more complex "mini-clusters", which must be accounted for in studying, for example, the internal dynamics of star clusters, and the physical stellar collision rate.

The separation distribution for M-dwarf binaries in the ASTRALUX survey is narrower and peaking at smaller separations than the distribution for solar-type binaries. This is often interpreted to mean that M-dwarfs constitute a continuous transition from brown dwarfs (BDs) to stars. Here a prediction for the M-dwarf separation distribution is presented, using a dynamical population synthesis (DPS) model in which "star-like" binaries with late-type primaries ($\lesssim1.5 M_{\rm sun}$) follow universal initial distribution functions and are dynamically processed in their birth embedded clusters. A separate "BD-like" population has both its own distribution functions for binaries and initial mass function (IMF), which overlaps in mass with the IMF for stars. Combining these two formation modes results in a peak on top of a wider separation distribution for late M-dwarfs consistent with the late ASTRALUX sample. The DPS separation distribution for early M-dwarfs shows no such peak and is in agreement with the M-dwarfs in Multiples (MinMS) data. We note that the latter survey is potentially in tension with the early ASTRALUX data. Concluding, the ASTRALUX and MinMS data are unable to unambiguously distinguish whether or not BDs are a continuous extension of the stellar IMF. Future observational efforts are needed to fully answer this interesting question. The DPS model predicts that binaries outside the sensitivity range of the ASTRALUX survey remain to be detected. For application to future data, we present a means to observationally measure the overlap of the putative BD-like branch and the stellar branch. We discuss the meaning of universal star formation and distribution functions.

In a virialized stellar system, the mean-square velocity is a direct tracer of the energy per unit mass of the system. Here, we exploit this to estimate and compare root-mean-square velocities for a large sample of nuclear star clusters and their host (late- or early-type) galaxies. Traditional observables, such as the radial surface brightness and second-order velocity moment profiles, are subject to short-term variations due to individual episodes of matter infall and/or star formation. The total mass, energy and angular momentum, on the other hand, are approximately conserved. Thus, the total energy and angular momentum more directly probe the formation of galaxies and their nuclear star clusters, by offering access to more fundamental properties of the nuclear cluster-galaxy system than traditional observables. We find that there is a strong correlation, in fact a near equality, between the root-mean-square velocity of a nuclear star cluster and that of its host. Thus, the energy per unit mass of a nuclear star cluster is always comparable to that of its host galaxy. We interpret this as evidence that nuclear star clusters do not form independently of their host galaxies, but rather that their formation and subsequent evolution are coupled. We discuss how our results can potentially be used to offer a clear and observationally testable prediction to distinguish between the different nuclear star cluster formation scenarios, and/or quantify their relative contributions.

We develop a formalism for calculating probabilities for the outcomes of stellar dynamical interactions, based on results from $N$-body scattering experiments. We focus here on encounters involving up to six particles and calculate probabilities for direct stellar collisions; however our method is in principle valid for larger particle numbers. Our method relies on the binomial theorem, and is applicable to encounters involving any combination of particle radii. We further demonstrate that our base model is valid to within a few percent for any combination of particle masses, provided the minimum mass ratio is within a factor of a few from unity. This method is particularly suitable for models of collisional systems involving large numbers of stars, such as globular clusters, old open clusters and galactic nuclei, where small subsets of stars may regularly have very close encounters, and the direct integration of all such encounters is computationally expensive. Variations of our method may also be used to treat other encounter outcomes, such as ejections and exchanges.

We will discuss the evolution of star clusters with an large initial binary fraction, up to 95%. The initial binary population is chosen to follow the invariant orbital-parameter distributions suggested by Kroupa (1995). The Monte Carlo MOCCA simulations of star cluster evolution are compared to the observations of Milone et al. (2012) for photometric binaries. It is demonstrated that the observed dependence on cluster mass of both the binary fraction and the ratio of the binary fractions inside and outside of the half mass radius are well recovered by the MOCCA simulations. This is due to a rapid decrease in the initial binary fraction due to the strong density-dependent destruction of wide binaries described by Marks, Kroupa & Oh (2011). We also discuss a new scenario for the formation of intermediate mass black holes in dense star clusters. In this scenario, intermediate mass black holes are formed as a result of dynamical interactions of hard binaries containing a stellar mass black hole, with other stars and binaries. We will discuss the necessary conditions to initiate the process of intermediate mass black hole formation and the dependence of its mass accretion rate on the global cluster properties.

We consider the effects that accretion from the interstellar medium onto the particles of an N-body system has on the rate of two-body relaxation. To this end, we derive an accretion-modified relaxation time by adapting Spitzer's two-component model to include the damping effects of accretion. We consider several different mass-dependencies and efficiency factors for the accretion rate, as well as different mass ratios for the two components of the model. The net effect of accretion is to accelerate mass segregation by increasing the average mass $\bar{m}$, since the relaxation time is inversely proportional to $\bar{m}$. Under the assumption that the accretion rate increases with the accretor mass, there are two additional effects that accelerate mass segregation. First, accretion acts to increase the range of any initial mass spectrum, quickly driving the heaviest members to even higher masses. Second, accretion acts to reduce the velocities of the accretors due to conservation of momentum, and it is the heaviest members that are affected the most. Using our two-component model, we quantify these effects as a function of the accretion rate, the total cluster mass, and the component masses. We conclude by discussing the implications of our results for the dynamical evolution of primordial globular clusters, primarily in the context of black holes formed from the most massive stellar progenitors.

In this paper, we discuss the origin of the observed correlation between cluster concentration c and present-day mass function (PDMF) slope {\alpha} reported by De Marchi, Paresce & Pulone. This relation can either be reproduced from universal initial conditions combined with some dynamical mechanism(s) that alter(s) the cluster structure and mass function over time, or it must arise early on in the cluster lifetime, such as during the gas-embedded phase of cluster formation. Using a combination of Monte Carlo and N-body models for globular cluster evolution performed with the MOCCA and NBODY6 codes, respectively, we explore a number of dynamical mechanisms that could affect the observed relation. For the range of initial conditions considered here, our results are consistent with an universal initial binary fraction \approx 10% (which does not, however, preclude 100%) and an universal initial stellar mass function resembling the standard Kroupa distribution. Most of the dispersion observed in the c-{\alpha} relation can be attributed to two-body relaxation and Galactic tides. However, dynamical processes alone could not have reproduced the dispersion in concentration, and we require at least some correlation between the initial concentration and the total cluster mass. We argue that the origin of this trend could be connected to the gas-embedded phase of cluster evolution.

(abridged) In this paper, we constrain the properties of primordial binary populations in Galactic globular clusters using the MOCCA Monte Carlo code for cluster evolution. Our results are compared to the observations of Milone et al. (2012) using the photometric binary populations as proxies for the true underlying distributions, in order to test the hypothesis that the data are consistent with an universal initial binary fraction near unity and the binary orbital parameter distributions of Kroupa (1995). With the exception of a few possible outliers, we find that the data are to first-order consistent with the universality hypothesis. Specifically, the present-day binary fractions inside the half-mass radius r$_{\rm h}$ can be reproduced assuming either high initial binary fractions near unity with a dominant soft binary component as in the Kroupa distribution combined with high initial densities (10$^4$-10$^6$ M$_{\odot}$ pc$^{-3}$), or low initial binary fractions ($\sim$ 5-10%) with a dominant hard binary component combined with moderate initial densities near their present-day values (10$^2$-10$^3$ M$_{\odot}$ pc$^{-3}$). This apparent degeneracy can be broken using the binary fractions outside r$_{\rm h}$- only high initial binary fractions with a significant soft component combined with high initial densities can contribute to reproducing the observed anti-correlation between the binary fractions outside r$_{\rm h}$ and the total cluster mass. We further illustrate using the simulated present-day binary orbital parameter distributions and the technique introduced in Leigh et al. (2012) that the relative fractions of hard and soft binaries can be used to further constrain the initial cluster density and mass-density relation. Our results favour an initial mass-density relation of the form r$_{\rm h} \propto$ M$_{\rm clus}^{\alpha}$ with $\alpha <$ 1/3.

In this paper, we address the question: What is the probability of stellar-mass black hole (BH) binaries co-existing in a globular cluster with an intermediate-mass black hole (IMBH)? Our results suggest that the detection of one or more BH binaries can strongly constrain the presence of an IMBH in most Galactic globular clusters. More specifically, the detection of one or more BH binaries could strongly indicate against the presence of an IMBH more massive than $\gtrsim 10^3$ M$_{\rm \odot}$ in roughly 80\% of the clusters in our sample. To illustrate this, we use a combination of N-body simulations and analytic methods to weigh the rate of formation of BH binaries against their ejection and/or disruption rate via strong gravitational interactions with the central (most) massive BH. The eventual fate of a sub-population of stellar-mass BHs (with or without binary companions) is for all BHs to be ejected from the cluster by the central IMBH, leaving only the most massive stellar-mass BH behind to form a close binary with the IMBH. During each phase of evolution, we discuss the rate of inspiral of the central BH-BH pair as a function of both the properties of the binary and its host cluster.

We use N-body simulations to explore the influence of orbital eccentricity on the dynamical evolution of star clusters. Specifically we compare the mass loss rate, velocity dispersion, relaxation time, and the mass function of star clusters on circular and eccentric orbits. For a given perigalactic distance, increasing orbital eccentricity slows the dynamical evolution of a cluster due to a weaker mean tidal field. However, we find that perigalactic passes and tidal heating due to an eccentric orbit can partially compensate for the decreased mean tidal field by energizing stars to higher velocities and stripping additional stars from the cluster, accelerating the relaxation process. We find that the corresponding circular orbit which best describes the evolution of a cluster on an eccentric orbit is much less than its semi-major axis or time averaged galactocentric distance. Since clusters spend the majority of their lifetimes near apogalacticon, the properties of clusters which appear very dynamically evolved for a given galactocentric distance can be explained by an eccentric orbit. Additionally we find that the evolution of the slope of the mass function within the core radius is roughly orbit-independent, so it could place additional constraints on the initial mass and initial size of globular clusters with solved orbits. We use our results to demonstrate how the orbit of Milky Way globular clusters can be constrained given standard observable parameters like galactocentric distance and the slope of the mass function. We then place constraints on the unsolved orbits of NGC 1261,NGC 6352, NGC 6496, and NGC 6304 based on their positions and mass functions.

Populations of multiple stars inside clustered regions are known to change through dynamical interactions. The efficiency of binary disruption is thought to be determined by stellar density. King and collaborators recently investigated the multiplicity properties in young star forming regions and in the Galactic field. They concluded that stellar density-dependent modification of a universal initial binary population (the standard or null hypothesis model) cannot explain the observations. We re-visit their results, analyzing the data within the framework of different model assumptions, namely non-universality without dynamical modification and universality with dynamics. We illustrate that the standard model does account for all known populations if regions were significantly denser in the past. Some of the effects of using present-day cluster properties as proxies for their past values are emphasized and that the degeneracy between age and density of a star forming region can not be omitted when interpreting multiplicity data. A new analysis of the Corona Australis region is performed within the standard model. It is found that this region is likely as unevolved as Taurus and an initial density of $\approx190M_{\odot}\;pc^{-3}$ is required to produce the presently observed binary population, which is close to its present-day density.

In this paper, we consider how gas damping affects the dynamical evolution of gas-embedded star clusters. Using a simple three-component (i.e. one gas and two stellar components) model, we compare the rates of mass segregation due to two-body relaxation, accretion from the interstellar medium, and gas dynamical friction in both the supersonic and subsonic regimes. Using observational data in the literature, we apply our analytic predictions to two different astrophysical environments, namely galactic nuclei and young open star clusters. Our analytic results are then tested using numerical simulations performed with the NBSymple code, modified by an additional deceleration term to model the damping effects of the gas. The results of our simulations are in reasonable agreement with our analytic predictions, and demonstrate that gas damping can significantly accelerate the rate of mass segregation. A stable state of approximate energy equilibrium cannot be achieved in our model if gas damping is present, even if Spitzer's Criterion is satisfied. This instability drives the continued dynamical decoupling and subsequent ejection (and/or collisions) of the more massive population. Unlike two-body relaxation, gas damping causes overall cluster contraction, reducing both the core and half-mass radii. If the cluster is mass segregated (and/or the gas density is highest at the cluster centre), the latter contracts faster than the former, accelerating the rate of core collapse.

Several possible formation pathways for blue straggler stars have been developed recently, but no one pathway has yet been observationally confirmed for a specific blue straggler. Here we report the first findings from a Hubble Space Telescope ACS/SBC far-UV photometric program to search for white dwarf companions to blue straggler stars. We find three hot and young white dwarf companions to blue straggler stars in the 7-Gyr open cluster NGC 188, indicating that mass transfer in these systems ended less than 300 Myr ago. These companions are direct and secure observational evidence that these blue straggler stars were formed through mass transfer in binary stars. Their existence in a well-studied cluster environment allows for observational constraints of both the current binary system and the progenitor binary system, mapping the entire mass transfer history.

Over the last few decades, observational surveys have revealed that high-order multiple-star systems (e.g. triples, quadruples, etc.), and triples in particular, are common in our Galaxy. In this paper, we consider the dynamical significance of this transformation in our understanding of stellar multiplicity. Using empirically constrained binary and triple fractions in those star clusters for which these values are available in the literature, we compare analytic rates for encounters involving single, binary, and triple stars. Our results show that, \textit{even for relatively low triple fractions, dynamical interactions involving triples occur roughly as often as encounters involving either single or binary stars alone, particularly in low-mass star clusters}. More generally, using empirically-derived multiple star catalogues for the young star-forming association Taurus-Auriga and the Galactic field, we show that the data are consistent with the gravitationally-focused cross section for encounters increasing with increasing multiplicity. Consequently, triple stars, and even higher-order multiples, could be more important than previously realized for a number of astrophysical phenomena, including the formation and destruction of compact binaries and various types of stellar exotica, and the dynamical evolution of star clusters.

(abridged) We consider the effect of stellar remnants on the interstellar medium of a massive star cluster following the initial burst of star formation. We argue that accretion onto stellar-mass black holes (BHs) is an effective mechanism for rapid gas depletion in clusters of all masses, as long as they contain progenitor stars more massive than \gtrsim 50\msun. This scenario is attractive for the progenitor systems of present-day massive globular clusters (GCs) which likely had masses M \gtrsim 10^7\msun. In such clusters, supernovae and stellar winds cannot provide a plausible explanation for the sudden removal of the primordial gas reservoir that is required to explain their complex chemical enrichment history. In order to consider different regimes in the gas accretion rate onto stellar-mass BHs, we consider both the Bondi-Hoyle and Eddington approximations. For either model, our results show that the gas can be significantly depleted within only a few tens of Myrs. This process will affect the distribution of BH masses, and may accelerate the dynamical decoupling of the BH population and, ultimately, their dynamical ejection. Moreover, the timescales for gas depletion are sufficiently short that the accreting BHs could significantly affect the chemistry of subsequent star formation episodes. The gas depletion times and final mass in BHs are sensitive to the assumed model for the accretion rate, and to the initial mass of the most massive BH which, in turn, is determined by the upper mass cut-off of the stellar IMF. Our results imply that the remnant accretion history can have an important bearing on the observed present-day cluster mass-to-light ratio. In particular, we show that an increase of the upper mass cut-off with decreasing metallicity could contribute to the observed anti-correlation between the mass-to-light ratio and the metallicity of GCs.

(abridged) We use newly available empirical binary fractions for globular clusters to carry out a direct test of the binary evolution hypothesis, and of collisional channels that involve binary stars. More specifically, using the previously reported correlation between blue straggler numbers and core mass as a benchmark, we test for correlations with the number of binary stars, as well as with the rates of single-single, single-binary, and binary-binary encounters. Surprisingly, we find that the simple correlation with core mass remains by far the strongest predictor of blue straggler population size, even in our joint models. This is despite the fact that the binary fractions themselves strongly anti-correlate with core mass, just as expected in the binary evolution model. At first sight, these results do not fit neatly with either binary evolution or collisional models in their simplest forms. Arguably the simplest and most intriguing possibility to explain this unexpected result is that observational errors on the core binary fractions are larger than the true intrinsic dispersion associated with their dependence on core mass. In the context of the binary evolution model, this would explain why the combination of binary fraction and core mass is a poorer predictor of blue straggler numbers than core mass alone. It would also imply that core mass is a remarkably clean predictor of core binary fractions. This would be of considerable importance for the dynamical evolution of globular clusters, and provides an important benchmark for models attempting to understand their present-day properties.

(Abridged) We combine published photometry for the nuclear star clusters (NSCs) and stellar spheroids of 51 low-mass, early-type galaxies in the Virgo cluster with empirical mass-to-light ratios, in order to complement previous studies that explore the dependence of NSC masses (M_{NSC}) on various properties of their host galaxies. We confirm a roughly linear relationship between M_{NSC} and luminous host spheroid mass (M_{Sph}), albeit with considerable scatter (0.57 dex). We estimate velocity dispersions from the virial theorem, assuming all galaxies in our sample share a common DM fraction and are dynamically relaxed. We then find that M_{NSC} \sim \sigma^{2.73\pm 0.29}, with a slightly reduced scatter of 0.54 dex. This confirms recent results that the shape of the M_{CMO} - \sigma relation is different for NSCs and super-massive black holes (SMBHs). We discuss this result in the context of the generalized idea of "central massive objects" (CMOs). In order to assess which physical parameters drive the observed NSC masses, we also carry out a joint multi-variate power-law fit to the data. In this, we allow M_{NSC} to depend on M_{Sph} and R_{Sph} (and hence implicitly on \sigma), as well as on the size of the globular cluster reservoir. When considered together, the dependences on M_{Sph} and R_{Sph} are roughly consistent with the virial theorem, and hence M_{NSC} \propto \sigma^2. However, the only statistically significant correlation is a linear scaling between M_{NSC} and M_{Sph}. We compare M_{NSC} with predictions for two popular models for NSC formation, namely i) globular cluster infall due to dynamical friction, and ii) in-situ formation during the early phases of galaxy formation that is regulated via momentum feedback from stellar winds and/or supernovae. Neither model can directly predict the observations, and we discuss possible interpretations of our findings.

In this paper, we study small-N gravitational dynamics involving up to six objects. We perform a large suite of numerical scattering experiments involving single, binary, and triple stars. This is done using the FEWBODY numerical scattering code, which we have upgraded to treat encounters involving triple stars. We focus on outcomes that result in direct physical collisions between stars, within the low angular momentum and high absolute orbital energy regime. The dependence of the collision probability on the number of objects involved in the interaction, N, is found for fixed total energy and angular momentum. Our results are consistent with a collision probability that increases approximately as N^2. Interestingly, this is also what is expected from the mean free path approximation in the limit of very large N. A more thorough exploration of parameter space will be required in future studies to fully explore this potentially intriguing connection. This study is meant as a first step in an on-going effort to extend our understanding of small-N collisional dynamics beyond the three- and four-body problems and into the realm of larger-N.

We present a new technique to quantify cluster-to-cluster variations in the observed present-day stellar mass functions of a large sample of star clusters. Our method quantifies these differences as a function of both the stellar mass and the total cluster mass, and offers the advantage that it is insensitive to the precise functional form of the mass function. We applied our technique to data taken from the ACS Survey for Globular Clusters, from which we obtained completeness-corrected stellar mass functions in the mass range 0.25-0.75 M$_{\odot}$ for a sample of 27 clusters. The results of our observational analysis were then compared to Monte Carlo simulations for globular cluster evolution spanning a range of initial mass functions, total numbers of stars, concentrations, and virial radii. We show that the present-day mass functions of the clusters in our sample can be reproduced by assuming an universal initial mass function for all clusters, and that the cluster-to-cluster differences are consistent with what is expected from two-body relaxation. A more complete exploration of the initial cluster conditions will be needed in future studies to better constrain the precise functional form of the initial mass function. This study is a first step toward using our technique to constrain the dynamical histories of a large sample of old Galactic star clusters and, by extension, star formation in the early Universe.

(abridged) This thesis presents theoretical and statistical techniques broadly related to systems of dynamically-interacting particles composed of several different types of populations. They are applied to observations of dense star clusters (SCs) in order to study gravitational interactions between stars. We present a new analytic method of quantifying the frequency of encounters involving single, binary and triple stars. With this technique, we have shown that dynamical encounters involving triple stars occur commonly in at least some SCs, and that they are likely to be an important dynamical channel for stellar mergers to occur. We have also used our techniques to analyze observational data for a large sample of SCs taken from the ACS Survey for Globular Clusters. The results of this analysis are as follows: (1) We have compiled a homogeneous catalogue of stellar populations for every cluster in our sample, including main-sequence (MS), red giant branch, horizontal branch and blue straggler (BS) stars. (2) With this catalogue, we have quantified the effects of the cluster dynamics in determining the relative sizes and spatial distributions of these stellar populations. (3) These results are particularly interesting for BSs since they provide compelling evidence that they are descended from binary stars. (4) Our analysis of the MS populations is consistent with a remarkably universal initial stellar mass function in old massive SCs in the Milky Way. This is a new result with important implications for our understanding of star formation in the early Universe and, more generally, the history of our Galaxy. Finally, we describe how our techniques are ideally suited for application to a number of other outstanding puzzles of modern astrophysics, including chemical reactions in the interstellar medium and mergers between galaxies in galaxy clusters.

We present an analytic model for blue straggler formation in globular clusters. We assume that blue stragglers are formed only through stellar collisions and binary star evolution, and compare our predictions to observed blue straggler numbers taken from the catalogue of Leigh, Sills & Knigge (2011). We can summarize our key results as follows: (1) Binary star evolution consistently dominates blue straggler production in all our best-fitting models. (2) In order to account for the observed sub-linear dependence of blue straggler numbers on the core masses (Knigge, Leigh & Sills 2009), the core binary fraction must be inversely proportional to the total cluster luminosity and should always exceed at least a few percent. (3) In at least some clusters, blue straggler formation must be enhanced by dynamical encounters (either via direct collisions or by stimulating mass-transfer to occur by altering the distribution of binary orbital parameters) relative to what is expected by assuming a simple population of binaries evolving in isolation. (4) The agreement between the predictions of our model and the observations can be improved by including blue stragglers that form outside the core but later migrate in due to dynamical friction. (5) Longer blue straggler lifetimes are preferred in models that include blue stragglers formed outside the core since this increases the fraction that will have sufficient time to migrate in via dynamical friction.

We present a homogeneous catalogue for blue straggler, red giant branch, horizontal branch and main-sequence turn-off stars in a sample of 35 clusters taken from the ACS Survey for Globular Clusters. As a result of the superior photometry and relatively large field of view offered by the ACS data, this new catalogue is a significant improvement upon the one presented in Leigh, Sills & knigge (2007). Using our catalogue, we study and compare the radial distributions of the different stellar populations. We have confirmed our previous result (Knigge, Leigh & Sills 2009) that there is a clear, but sub-linear, correlation between the number of blue stragglers found in the cluster core and the total stellar mass contained within it. By considering a larger spatial extent than just the core, our results suggest that mass segregation is not the dominant effect contributing to the observed sub-linearity. We also investigate the radial distributions of the different stellar populations in our sample of clusters. Our results are consistent with a linear relationship between the number of stars in these populations and the total mass enclosed within the same radius. Therefore, we conclude that the cluster dynamics does not significantly affect the relative distributions of these populations in our sample.

We present a technique to identify the most probable dynamical formation scenario for an observed binary or triple system containing one or more merger products or, alternatively, to rule out the possibility of a dynamical origin. Our method relies on an analytic prescription for energy conservation during stellar encounters. With this, observations of the multiple star system containing the merger product(s) can be used to work backwards in order to constrain the initial orbital energies of any single, binary or triple systems that went into the encounter. The initial semi-major axes of the orbits provide an estimate for the collisional cross section and therefore the time-scale for the encounter to occur in its host cluster. We have applied our analytic prescription to observed binary and triple systems containing blue stragglers, in particular the triple system S1082 in M67 and the period distribution of the blue straggler binaries in NGC 188. We have shown that both S1082 and most of the blue straggler binaries in NGC 188 could have a dynamical origin, and that encounters involving triples could be a significant contributor to BS populations in old open clusters. In general, our results suggest that encounters involving triples could make up a significant fraction of those dynamical interactions that result in stellar mergers, in particular encounters that produce multiple star systems containing one or more blue stragglers.

We present a technique to identify the most probable dynamical formation scenario for an observed binary or triple system containing one or more merger products or, alternatively, to rule out the possibility of a dynamical origin. Our method relies on an analytic prescription for energy conservation during stellar encounters. With this, observations of the multiple star system containing the merger product(s) can be used to work backwards in order to constrain the initial orbital energies of any single, binary or triple systems that went into the encounter. The initial semi-major axes of the orbits provide an estimate for the collisional cross section and therefore the time-scale for the encounter to occur in its host cluster. We have applied our analytic prescription to observed binary and triple systems containing blue stragglers, in particular the triple system S1082 in M67 and the period distribution of the blue straggler binaries in NGC 188. We have shown that both S1082 and most of the blue straggler binaries in NGC 188 could have a dynamical origin, and that encounters involving triples are a significant contributor to blue straggler populations in old open clusters. In general, our results suggest that encounters involving triples could make up a significant fraction of those dynamical interactions that result in stellar mergers, in particular encounters that produce multiple star systems containing one or more blue stragglers.

It has become clear in recent years that globular clusters are not simple stellar populations, but may host chemically distinct sub-populations, typically with an enhanced helium abundance. These helium-rich populations can make up a substantial fraction of all cluster stars. One of the proposed formation channels for blue straggler stars is the physical collision and merger of two stars. In the context of multiple populations, collisions between stars with different helium abundances should occur and contribute to the observed blue straggler population. This will affect the predicted blue straggler colour and luminosity function. We quantify this effect by calculating models of mergers resulting from collisions between stars with different helium abundances and using these models to model a merger population. We then compare these results to four observed clusters, NGC 1851, NGC 2808, NGC 5634 and NGC 6093. As in previous studies our models deviate from the observations, particularly in the colour distributions. However, our results are consistent with observations of multiple populations in these clusters. In NGC 2808, our best fitting models include normal and helium enhanced populations, in agreement with helium enhancement inferred in this cluster. The other three clusters show better agreement with models that do not include helium enhancement. We discuss future prospects to improve the modelling of blue straggler populations and the role that the models we present here can play in such a study.

We investigate the relationship between the mass of a globular cluster core and the sizes of its various stellar populations in a sample of 56 globular clusters. The number of core red giant branch stars is found to scale sub-linearly with core mass at the 3-$\sigma$ confidence level, whereas the relation is linear to within one standard deviation for main-sequence turn-off and sub-giant branch stars. We interpret our results as evidence for a surplus of red giant branch stars in the least massive cluster cores which is not seen for main-sequence turn-off and sub-giant branch stars. We explore various possibilities for the source of this discrepancy, discussing our results primarily in terms of the interplay between the cluster dynamics and stellar evolution.

Blue stragglers in globular clusters are abnormally massive stars that should have evolved off the stellar main sequence long ago. There are two known processes that can create these objects: direct stellar collisions and binary evolution. However, the relative importance of these processes has remained unclear. In particular, the total number of blue stragglers found in a given cluster does not seem to correlate with the predicted collision rate, providing indirect support for the binary-evolution model. Yet the radial distributions of blue stragglers in many clusters are bimodal, with a dominant central peak: this has been interpreted as an indication that collisions do dominate blue straggler production, at least in the high-density cluster cores. Here we report that there is a clear, but sublinear, correlation between the number of blue stragglers found in a cluster core and the total stellar mass contained within it. From this we conclude that most blue stragglers, even those found in cluster cores, come from binary systems. The parent binaries, however, may themselves have been affected by dynamical encounters. This may be the key to reconciling all of the seemingly conflicting results found to date.

The formation of blue stragglers is still not completely understood, particularly the relationship between formation environment and mechanism. We use a large, homogeneous sample of blue stragglers in the cores of 57 globular clusters to investigate the relationships between blue straggler populations and their environments. We use a consistent definition of "blue straggler" based on position in the color-magnitude diagram, and normalize the population relative to the number of red giant branch stars in the core. We find that the previously determined anti-correlation between blue straggler frequency and total cluster mass is present in the purely core population. We find some weak correlations with central velocity dispersion and with half-mass relaxation time. The blue straggler frequency does not show any trend with any other cluster parameter. Even though collisions may be expected to be a dominant blue straggler formation process in globular cluster cores, we find no correlation between the frequency of blue stragglers and the collision rate in the core. We also investigated the blue straggler luminosity function shape, and found no relationship between any cluster parameter and the distribution of blue stragglers in the color-magnitude diagram. Our results are inconsistent with some recent models of blue straggler formation that include collisional formation mechanisms, and may suggest that almost all observed blue stragglers are formed in binary systems.