Search

Among the growing list of confirmed exoplanets, the number of planets identified in dense star clusters remains sparse. Previous analyses have suggested that this may be due in part to dynamical interactions that can unbind planets from their host stars, limiting the survival of planetary systems in clusters. Thus, alternative detection strategies may be necessary to study planets in clusters that may no longer be bound to a host. Here, we use N-body models to explore the evolution of planetary systems in dense star clusters. Depending on various initial conditions, we show that 10%–50% of primordial planetary systems are broken through dynamical encounters over a cluster’s full lifetime, populating clusters with “free-floating” planets. Furthermore, a large number (30%–80%) of planets are ejected from their host cluster through strong dynamical encounters and/or tidal loss. Additionally, we show that planets naturally mix with stellar-mass black holes (BHs) in the cores of their host cluster. As a consequence, up to a few hundred planets will be tidally disrupted through close passages of BHs. We show that these BH–planet tidal disruption events (TDEs) occur in clusters at a rate of up to 10−5 yr−1 in a Milky-Way-type galaxy. In principle, these BH–planet TDEs may be detected by upcoming transient surveys such as the Large Synoptic Survey Telescope at a rate of a few events per year, although identification of these events may prove challenging. The observed rate of BH–planet TDEs could place new constraints upon the formation and survival of planetary systems and BHs in dense star clusters. [ABSTRACT FROM AUTHOR]

In a dense stellar environment, such as the core of a globular cluster (GC), dynamical interactions with black holes (BHs) are expected to lead to a variety of astrophysical transients. Here we explore tidal disruption events (TDEs) of stars by stellar-mass BHs through collisions and close encounters. Using state-of-the-art cluster simulations, we show that these TDEs occur at significant rates throughout the evolution of typical GCs and we study how their relative rates relate to cluster parameters such as mass and size. By incorporating a realistic cosmological model of GC formation, we predict a BH–main-sequence-star TDE rate of approximately 3 Gpc−3 yr−1 in the local universe (z < 0.1) and a cosmological rate that peaks at roughly 25 Gpc−3 yr−1 for redshift 3. Furthermore, we show that the ejected mass associated with these TDEs could produce optical transients of luminosity ∼1041−1044 erg s−1 with timescales of about a day to a month. These should be readily detectable by optical transient surveys such as the Zwicky Transient Facility. Finally, we comment briefly on BH–giant encounters and discuss how these events may contribute to the formation of BH–white-dwarf binaries. [ABSTRACT FROM AUTHOR]

Over 100 millisecond radio pulsars (MSPs) have been observed in globular clusters (GCs), motivating theoretical studies of the formation and evolution of these sources through stellar evolution coupled to stellar dynamics. Here we study MSPs in GCs using realistic N-body simulations with our Cluster Monte Carlo code. We show that neutron stars (NSs) formed in electron-capture supernovae (including both accretion-induced and merger-induced collapse of white dwarfs) can be spun up through mass transfer to form MSPs. Both NS formation and spin-up through accretion are greatly enhanced through dynamical interaction processes. We find that our models for average GCs at the present day with masses ≈2 × 105M⊙ can produce up to 10–20 MSPs, while a very massive GC model with mass ≈106M⊙ can produce close to 100. We show that the number of MSPs is anti-correlated with the total number of stellar-mass black holes (BHs) retained in the host cluster. The radial distributions are also affected: MSPs are more concentrated toward the center in a host cluster with a smaller number of retained BHs. As a result, the number of MSPs in a GC could be used to place constraints on its BH population. Some intrinsic properties of MSP systems in our models (such as the magnetic fields and spin periods) are in good overall agreement with observations, while others (such as the distribution of binary companion types) are less so, and we discuss the possible reasons for such discrepancies. Interestingly, our models also demonstrate the possibility of dynamically forming NS–NS and NS–BH binaries in GCs, although the predicted numbers are very small. [ABSTRACT FROM AUTHOR]

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 . 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 M⊙, 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 it 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. [ABSTRACT FROM AUTHOR]

Globular clusters (GCs) in the Milky Way exhibit a well-observed bimodal distribution in core radii separating the so-called core-collapsed and non-core-collapsed clusters. Here, we use our Hénon-type Monte Carlo code, CMC, to explore initial cluster parameters that map into this bimodality. Remarkably, we find that by varying the initial size of clusters (specified in our initial conditions in terms of the initial virial radius, rv) within a relatively narrow range consistent with the measured radii of young star clusters in the local universe (rv ≈ 0.5–5 pc), our models reproduce the variety of present-day cluster properties. Furthermore, we show that stellar-mass black holes (BHs) play an intimate role in this mapping from initial conditions to the present-day structural features of GCs. We identify “best-fit” models for three GCs with known observed BH candidates, NGC 3201, M22, and M10, and show that these clusters harbor populations of ∼50–100 stellar-mass BHs at present. As an alternative case, we also compare our models to the core-collapsed cluster NGC 6752 and show that this cluster likely contains few BHs at present. Additionally, we explore the formation of BH binaries in GCs and demonstrate that these systems form naturally in our models in both detached and mass-transferring configurations with a variety of companion stellar types, including low-mass main-sequence stars, white dwarfs, and sub-subgiants. [ABSTRACT FROM AUTHOR]

We explore the formation of mass-transferring binary systems containing black holes (BHs) within globular clusters (GC). We show that it is possible to form mass-transferring BH binaries with main sequence, giant, and white dwarf companions with a variety of orbital parameters in GCs spanning a large range in present-day properties. All mass-transferring BH binaries found in our models at late times are dynamically created. The BHs in these systems experienced a median of ∼30 dynamical encounters within the cluster before and after acquiring the donor. Furthermore, we show that the presence of mass-transferring BH systems has little correlation with the total number of BHs within the cluster at any time. This is because the net rate of formation of BH–non-BH binaries in a cluster is largely independent of the total number of retained BHs. Our results suggest that the detection of a mass-transferring BH binary in a GC does not necessarily indicate that the host cluster contains a large BH population. [ABSTRACT FROM AUTHOR]

We explore the long-term evolution of mass-transferring white dwarf (WD) binaries undergoing both direct-impact and disk accretion and explore implications of such systems to gravitational-wave (GW) astronomy. We cover a broad range of initial component masses and show that these systems, the majority of which lie within the Laser Interferometer Space Antenna (LISA) sensitivity range, exhibit prominent negative orbital frequency evolution (chirp) for a significant fraction of their lifetimes. Using a galactic population synthesis, we predict ∼2700 of these systems will be observable with a negative chirp of 0.1 yr−2 by a space-based GW detector like LISA. We also show that detections of mass-transferring double WD systems by LISA may provide astronomers with unique ways of probing the physics governing close compact object binaries. [ABSTRACT FROM AUTHOR]

We calculate the long-term evolution of angular momentum in double white dwarf binaries undergoing direct impact accretion over a broad range of parameter space. We allow the rotation rate of both components to vary and account for the exchange of angular momentum between the spins of the white dwarfs and the orbit, while conserving the total angular momentum. We include gravitational, tidal, and mass transfer effects in the orbital evolution, and allow the Roche radius of the donor star to vary with both the stellar mass and the rotation rate. We examine the long-term stability of these systems, focusing in particular on those systems that may be progenitors of AM CVn or SNe Ia. We find that our analysis yields an increase in the predicted number of stable systems compared to that in previous studies. Additionally, we find that by properly accounting for the effects of asynchronism between the donor and the orbit on the Roche-lobe size, we eliminate oscillations in the orbital parameters, which were found in previous studies. Removing these oscillations can reduce the peak mass transfer rate in some systems, keeping them from entering an unstable mass transfer phase. [ABSTRACT FROM AUTHOR]