Source author record

Nathan W. C. Leigh

Nathan W. C. Leigh appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

28works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

28 published item(s)

preprint2022arXiv

Tidal Disruption Events by Compact Supermassive Black Hole Binaries

Stars can be tidally destroyed or swallowed by supermassive black hole binaries. Using a large number of accurate few-body simulations, we investigate the enhancement and suppression of full and partial disruption and direct capture events by hard supermassive black hole binaries with a wide ranges of key parameters, i.e., the primary black hole mass ($10^{5}-10^{8}M_{\odot}$), the binary mass ratio ($10^{-3}-1$), the ratio of the binary semimajor axis to the hardening radius ($10^{-4}-1$), the binary eccentricity ($0.0-0.9$) and the stellar mass ($0.3-3M_{\odot}$). This is a significant extension of the parameter space compared to previous work. We show that the encounter probabilities of all three events are well-described by the encounter cross section, which is proportional to the pericenter distance. The probability of full disruptions by supermassive black hole binaries can be enhanced by up to a factor of $40-50$ or suppressed by up to a factor of $10$, relative to that by single black holes, depending on the binary parameters. Relativistic effects are not important for the primary black hole mass $\leq 10^{7}M_{\odot}$, but can provide an additional enhancement of the full disruption probability by less than a factor of $2-3$ for higher primary black hole masses. We provide a fitting formula for the full disruption probability by the hard supermassive black hole binaries that works for a wide range of parameters. We also find that partial disruption events can occur multiple times before full disruptions or direct captures, and their probabilities can be greater than that of full disruption events by a factor of three. Because partial disruption events can induce stellar spins and mass loss and change the orbits, it can significantly affect the overall full disruption event rate and the shape of the light curves.

preprint2021arXiv

Evolution of Retrograde Orbiters in an AGN Disk

AGN disks have been proposed as promising locations for the mergers of stellar mass black hole binaries (BBHs). Much recent work has been done on this merger channel, but the majority focuses on stellar mass black holes (BHs) orbiting in the prograde direction. Little work has been done to examine the impact of retrograde orbiters (ROs) on the formation and mergers of BBHs in AGN disks. Quantifying the retrograde contribution is important, since roughly half of all orbiters should initially be on retrograde orbits when the disk forms. We perform an analytic calculation of the evolution of ROs in an AGN disk. Because this evolution could cause the orbits of ROs to cross those of prograde BBHs, we derive the collision rate between a given RO and a given BBH orbiting in the prograde direction. In the examples given here, ROs in the inner region of the disk experience a rapid decrease in the semimajor axis of their orbits while also becoming highly eccentric in less than a million years. This rapid orbital evolution could lead to extreme mass ratio inspirals detectable by the Laser Interferometer Space Antenna. The collision rates of our example ROs with prograde BBHs in the migration trap depend strongly on the volume of the inner radiation-pressure-dominated region which depends on the mass of the supermassive black hole (SMBH). Rates are lowest for larger mass SMBHs, which dominate the AGN merger channel, suggesting that merger rates for this channel may not be significantly altered by ROs.

preprint2021arXiv

Symmetry Breaking in Dynamical Encounters in the Disks of Active Galactic Nuclei

Active galactic nucleus (AGN) disks may be important sites of binary black hole (BBH) mergers. Here we show via numerical experiments with the high-accuracy, high precision code {\tt SpaceHub} that broken symmetry in dynamical encounters in AGN disks can lead to an asymmetry between prograde and retrograde BBH mergers. The direction of the hardening asymmetry depends on the initial binary semi-major axis. An asymmetric distribution of mass-weighted projected spin $χ_{\rm eff}$ should therefore be expected in LIGO-Virgo detections of BBH mergers from AGN disks. This channel further predicts that negative $χ_{\rm eff}$ BBH mergers are most likely for massive binaries.

preprint2020arXiv

Chaos and Lévy Flights in the Three-Body Problem

We study chaos and Lévy 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évy 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évy flight interactions follow a power-law distribution. Lévy 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évy 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évy 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.

preprint2020arXiv

Electromagnetic transients and gravitational waves from white dwarf disruptions by stellar black holes in triple systems

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.

preprint2020arXiv

Filament rotation in the California L1482 cloud

We analyze the gas mass distribution, the gas kinematics, and the young stellar object (YSO) content of the California Molecular Cloud (CMC) L1482 filament. We derive a Gaia DR2 YSO distance of 511$^{+17}_{-16}$ pc. We derive scale-free power-laws for the mean gas line-mass (M/L) profiles; we calculate the gravitational potential and field profiles consistent with these. We present IRAM 30 m C$^{18}$O (1-0) (and other tracers) position-velocity (PV) diagrams that exhibit complex velocity twisting and turning structures. We find a rotational profile in C$^{18}$O perpendicular to the southern filament ridgeline. The profile is regular, confined ($r\lesssim0.4$ pc), anti-symmetric, and to first order linear with a break at $r\sim0.25$ pc. The timescales of the inner (outer) gradients are $\sim$0.7 (6.0) Myr. We show that the centripetal force, compared to gravity, increases toward the break; when the ratio of forces approaches unity, the profile turns over, just before filament breakup is achieved. The timescales and relative roles of gravity to rotation indicate that the structure is stable, long lived ($\sim$ a few times 6 Myr), and undergoing outside-in evolution. Moreover, this filament has practically no star formation, a perpendicular Planck plane-of-the-sky (POS) magnetic field morphology, and POS "zig-zag" morphology, which together with the rotation profile lead to the suggestion that the 3D shape is a corkscrew filament with a helical magnetic field. These results, combined with results in Orion and G035.39-00.33, suggest evolution toward higher densities as rotating filaments shed angular momentum. Thus, magnetic fields may be an essential feature of high-mass (M $\sim10^5$ M$_{\odot}$) cloud filament evolution toward cluster formation.

preprint2020arXiv

Giant planet swaps during close stellar encounters

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 $Γ\sim 4.4(M_{\rm c}/100M_\odot)^{-2}(a_{\rm p}/{\rm AU}) (σ/1\,{\rm km}\,{\rm s}^{-1})^{5}$Gyr$^{-1}$, where $a_{\rm p}$ is the planet semi-major axis around the birth star, $σ$ 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.

preprint2020arXiv

Hot Jupiter and ultra-cold Saturn formation in dense star clusters

The discovery of high incidence of hot Jupiters in dense clusters challenges the field-based hot Jupiter formation theory. In dense clusters, interactions between planetary systems and flyby stars are relatively common. This has a significant impact on planetary systems, dominating hot Jupiter formation. In this paper, we perform high precision, few-body simulations of stellar flybys and subsequent planet migration in clusters. A large parameter space exploration demonstrates that close flybys that change the architecture of the planetary system can activate high eccentricity migration mechanisms: Lidov-Kozai and planet-planet scattering, leading to high hot Jupiter formation rate in dense clusters. Our simulations predict that many of the hot Jupiters are accompanied by "ultra-cold Saturns", expelled to apastra of thousands of AU. This increase is particularly remarkable for planetary systems originally hosting two giant planets with semi-major axis ratios $\sim$ 4 and the flyby star approaching nearly perpendicular to the planetary orbital plane. The estimated lower limit to the hot Jupiter formation rate of a virialized cluster is $\sim 1.6\times10^{-4}(σ/{\rm 1kms^{-1}})^5({a_{\rm p}}/{\rm 20 AU})({M_{\rm c}}/{\rm 1000M_\odot})^{-2}$Gyr$^{-1}$ per star, where $σ$ is the cluster velocity dispersion, $a_{\rm p}$ is the size of the planetary system and $M_{\rm c}$ is the mass of the cluster. Our simulations yield a hot Jupiter abundance which is $\sim$ 50 times smaller than that observed in the old open cluster M67. We expect that interactions involving binary stars, as well as a third or more giant planets, will close the discrepancy.

preprint2020arXiv

Kinematic study of the association Cyg OB3 with Gaia DR2

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.

preprint2020arXiv

Massive Stellar Triples Leading to Sequential Binary Black-Hole Mergers in the Field

Stellar triples with massive stellar components are common, and can lead to sequential binary black-hole mergers. Here, we outline the evolution towards these sequential mergers, and explore these events in the context of gravitational-wave astronomy and the pair-instability mass gap. We find that binary black-hole mergers in the pair-instability mass gap can be of triple origin and therefore are not exclusively formed in dense dynamical environments. We discuss the sequential merger scenario in the context of the most massive gravitational-wave sources detected to date: GW170729 and GW190521. We propose that the progenitor of GW170729 is a low-metallicity field triple. We support the premise that GW190521 could not have been formed in the field. We conclude that triple stellar evolution is fundamental in the understanding of gravitational-wave sources, and likely, other energetic transientsas well.

preprint2020arXiv

On the eccentricity evolution of massive black hole binaries in stellar backgrounds

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$^*$).

preprint2020arXiv

Planetary Architectures in Interacting Stellar Environments

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.

preprint2019arXiv

A Statistical Solution to the Chaotic, Non-Hierarchical Three-Body Problem

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.

preprint2019arXiv

Constraining Mass Transfer Histories of Blue Straggler Stars with COS Spectroscopy of White Dwarf Companions

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.

preprint2019arXiv

The Keplerian three-body encounter II. Comparisons with isolated encounters and impact on gravitational wave merger timescales

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.

preprint2016arXiv

Interrupted Binary Mass Transfer in Star Clusters

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

preprint2016arXiv

On the origins of enigmatic stellar populations in Local Group galactic nuclei

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.

preprint2016arXiv

The chaotic four-body problem in Newtonian gravity I: Identical point-particles

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.

preprint2016arXiv

When does a star cluster become a multiple star system? I. Lifetimes of equal-mass small-N systems

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.

preprint2015arXiv

Back to the Future: Estimating Initial Globular Cluster Masses from their Present Day Stellar Mass Functions

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.

preprint2015arXiv

Interrupted Stellar Encounters in Star Clusters

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.

preprint2015arXiv

Small-N collisional dynamics II: Roaming the realm of not-so-small-N

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.

preprint2014arXiv

Modifying two-body relaxation in N-body systems by gas accretion

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.

preprint2014arXiv

On the coexistence of stellar-mass and intermediate-mass black holes in globular clusters

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.

preprint2014arXiv

Stellar dynamics in gas: The role of gas damping

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.

preprint2014arXiv

The state of globular clusters at birth II: primordial binaries

(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}^α$ with $α<$ 1/3.

preprint2014arXiv

The state of globular clusters at birth: emergence from the gas-embedded phase

In this paper, we discuss the origin of the observed correlation between cluster concentration c and present-day mass function (PDMF) slope α 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-α 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.

preprint2013arXiv

The Dynamical Significance of Triple Star Systems in Star Clusters

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.