Theory and observations suggest that single-star evolution is not able to produce black holes (BHs) with masses in the range $3-5M_{\odot}$ and above $\sim 45M_{\odot}$, referred to as the lower mass gap (LMG) and the upper mas gap (UMG), respectively. However, it is possible to form BHs in these gaps through merger of compact objects in dense clusters, e.g. the LMG and the UMG can be populated through binary neutron star- and BBH mergers, respectively. This implies that if binary mergers are observed in gravitational waves (GWs) with at least one mass gap object, then either clusters are effective in assembling binary mergers, or our single-star models have to be revised. Understanding how effective clusters are at populating both mass gaps have therefore major implications for both stellar- and GW astrophysics. In this paper we present a systematic study on how efficient stellar clusters are at populating both mass gaps through in-cluster GW mergers. For this, we derive a set of closed form relations for describing the evolution of compact object binaries undergoing dynamical interactions and GW merger inside their cluster. By considering both static and time evolving populations, we find in particular that globular clusters are clearly inefficient at populating the LMG in contrast to the UMG. We further describe how these results relate to the characteristic mass, time, and length scales associated with the problem.

We study the formation rate of binary black hole mergers formed through gravitational-wave emission between unbound, single black holes in globular clusters. While the formation of these binaries in very dense systems such as galactic nuclei has been well studied, we show here that this process can operate in lower-density stellar systems as well, forming binaries at a rate similar to other proposed pathways for creating eccentric mergers. Recent advances in post-Newtonian cluster dynamics indicate that a large fraction of dynamically-assembled binary black holes merge inside their host clusters during weak and strong binary-single and binary-binary interactions, and that these systems may retain measurable eccentricities as they travel through the LIGO and LISA sensitivity bands. Using an analytic approach to modeling binary black holes from globular clusters, we show that the formation of merging binaries from previously unbound black holes can operate at a similar rate to mergers forming during strong binary encounters, and that these binaries inhabit a unique region of the gravitational-wave frequency space which can be identified by proposed deci-Hertz space-based detectors.

We estimate the rate of gravitational microlensing events of cluster stars due to black holes (BHs) in the globular cluster NGC 5139 ($\omega Cen$). Theory and observations both indicate that $\omega Cen$ may contain thousands of BHs, but their mass spectrum and exact distribution are not well constrained. In this Letter we show that one may observe microlensing events on a timescale of years in $\omega Cen$, and such an event sample can be used to infer the BH distribution. Direct detection of BHs will, in the near future, play a major role in distinguishing binary BH merger channels. Here we explore how gravitational microlensing can be used to put constraints on BH populations in globular clusters.

Binary-binary interactions are important in a number of astrophysical contexts including dense stellar systems such as globular clusters. Although less frequent than binary-single encounters, binary-binary interactions lead to a much richer range of possibilities such as the formation of stable triple systems. Here, we focus on the regime of distant binary-binary encounters, i.e., two binaries approaching each other on an unbound orbit with a periapsis distance Q much larger than the internal binary separations. This `secular' regime gives rise to changes in the orbital eccentricities and orientations, which we study using analytic considerations and numerical integrations. We show that `direct' interactions between the three orbits only occur starting at a high expansion order of the Hamiltonian (hexadecupole order), and that the backreaction of the outer orbit on the inner two orbits at lower expansion orders is weak. Therefore, to good approximation, one can obtain the changes of each orbit by using previously-known analytic results for binary-single interactions, and replacing the mass of the third body with the total mass of the companion binary. Nevertheless, we find some dependence of the `binarity' of the companion binary, and derive explicit analytic expressions for the secular changes that are consistent with numerical integrations. In particular, the eccentricity and inclination changes of orbit 1 due to orbit 2 scale as \eps_{SA1} (a_2/Q)^2 [m_3 m_4/(m_3+m_4)^2], where \eps_{SA1} is the approximate quadrupole-order change, and a_2 and (m_3,m_4) are the companion binary orbital semimajor axis and component masses, respectively. Our results are implemented in several Python scripts that are freely available.

One of the proposed channels of binary black hole mergers involves dynamical interactions of three black holes. In such scenarios, it is possible that all three black holes merge in a so-called hierarchical merger chain, where two of the black holes merge first and then their remnant subsequently merges with the remaining single black hole. Depending on the dynamical environment, it is possible that both mergers will appear within the observable time window. Here we perform a search for such merger pairs in the public available LIGO and Virgo data from the O1/O2 runs. Using a frequentist p-value assignment statistics we do not find any significant merger pair candidates. Assuming no observed candidates in O3/O4, we derive upper limits on merger pairs to be $\sim11-110\ {\rm year^{-1}Gpc^{-3}}$, corresponding to a rate that relative to the total merger rate is $\sim 0.1-1.0$. From this we argue that both a detection and a non-detection within the next few years can be used to put useful constraints on some dynamical progenitor models.

Among the growing list of confirmed exoplanets, the number of planets identified in dense star clusters remains sparse. Previous analyses have suggested this may be due in part to dynamical interactions that unbind large fractions of 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 star. Here, we use the cluster Monte Carlo code CMC to explore the evolution of planetary systems in dense star clusters. Depending on a number of 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 central regions of their host cluster. As a consequence, up to a few hundreds of planets will be tidally disrupted through close passages of BHs. We show these BH-planet tidal disruption events (TDEs) occur in clusters at a rate of up to $10^{-5}\,\rm{yr}^{-1}$ in a Milky Way-type galaxy. We predict BH-planet TDEs should be detected by upcoming transient surveys such as LSST at a rate of roughly a few events per year. The observed rate of BH-planet TDEs would place new constraints upon the formation and survival of both planetary systems and BHs in dense star clusters. Additionally, depending on various assumptions including the initial number of planets and their orbital properties, we predict that typical globular clusters may contain a few dynamically-formed NS-planet systems at present as well as up to roughly 100 dynamically-formed WD-planet systems.

Dense stellar systems such as globular clusters are believed to harbor merging binary black holes (BHs). The evolution of such binaries is driven by interactions with other stars, most notably, binary-single interactions. Traditionally, so-called `strong' interactions are believed to be the driving force in this evolution. However, we recently showed that more distant, i.e., "weak" or "secular" encounters, can have important implications for the properties of merging BH binaries in globular clusters. This motivates more detailed understanding of the effects of secular encounters on a binary. In another previous paper, we analytically calculated expressions for the changes of the eccentricity and angular-momentum vectors taking into account second-order (SO) perturbation theory, and showed that, for highly eccentric binaries, the new expressions give rise to behavior that is not captured by first-order (FO) theory. Here, we extend our previous work to third order (TO) perturbation theory. We also include terms up to and including octupole order. The latter are nonzero for binaries with unequal component masses. In addition, we consider the effects of post-Newtonian terms, and we determine the steady-state distribution due to the cumulative effect of secular encounters by computing the associated angular-momentum diffusion coefficients, and applying the Fokker-Planck equation. Together with our previous work, the results in this paper provide a framework for incorporating the effects of distant encounters on binaries in models of cluster evolution, such as Monte Carlo codes.

Binary-single interactions play a crucial role in the evolution of dense stellar systems such as globular clusters. In addition, they are believed to drive black hole (BH) binary mergers in these systems. A subset of binary-single interactions are secular encounters, for which the third body approaches the binary on a relatively wide orbit, and such that it is justified to average the equations of motion over the binary's orbital phase. Previous works used first-order perturbation theory to compute the effects of such secular encounters on the binary. However, this approach can break down for highly eccentric binaries, which are important for BH binary mergers and gravitational wave sources. Here, we present an analytic computation using second-order perturbation techniques, valid to the quadrupole-order approximation. In our calculation, we take into account the instantaneous back-reaction of the binary to the third body, and compute corrections to previous first-order results. Using singly-averaged and direct 3-body integrations, we demonstrate the validity of our expressions. In particular, we show that the eccentricity change for highly eccentric binaries can reach a plateau, associated with a large inclination change, and can even reverse sign. These effects are not captured by previous first-order results. We provide a simple script to conveniently evaluate our analytic expressions, including routines for numerical integration and verification.

We present the first systematic study on how distant weak interactions impact the dynamical evolution of merging binary black holes (BBHs) in dense stellar clusters. Recent studies indicate that dense clusters are likely to significantly contribute to the rate of merging BBHs observable through gravitational waves (GWs), and that many of these mergers will appear with notable eccentricities measurable in the LISA and LIGO sensitivity bands. This is highly interesting, as eccentricity can be used to distinguish between different astrophysical merger channels. However, all of these recent studies are based on various Monte Carlo (MC) techniques that only include strong interactions for the dynamical evolution of BBHs, whereas any binary generally undergoes orders-of-magnitude more weak interactions than strong. It is well known that weak interactions primarily lead to a change in the binary's eccentricity, which for BBHs implies that weak interactions can change their GW inspiral time and thereby their merger probability. With this motivation, we perform MC simulations of BBHs evolving in dense clusters under the influence of both weak and strong interactions. We find that including weak interactions leads to a notable increase in the number of BBHs that merge inside their cluster, which correspondingly leads to a higher number of eccentric LISA sources. These preliminary results illustrate the importance of including weak interactions for accurately modeling how BBHs merge in clusters, and how to link their emitted GW signals to their astrophysical environment.

We derive the observable gravitational wave (GW) peak frequency ($f$) distribution of binary black holes (BBHs) that currently reside inside their globular clusters (GCs), with and without 2.5 Post-Newtonian (2.5PN) effects included in the dynamical evolution of the BBHs. Recent Newtonian studies have reported that a notable number of nearby non-merging BBHs, i.e. those BBHs that are expected to undergo further dynamical interactions before merger, in GCs are likely to be observable by LISA. However, our 2.5PN calculations show that the distribution of $\log f$ for the non-merging BBH population above $\sim 10^{-3.5}$ Hz scales as $f^{-34/9}$ instead of the $f^{-2/3}$ scaling found in the Newtonian case. This leads to an approximately two-orders-of-magnitude reduction in the expected number of GW sources at $\sim 10^{-3}$ Hz, which lead us to conclude that observing nearby BBHs with LISA is not as likely as has been claimed in the recent literature. In fact, our results suggest that it might be more likely that LISA detects the population of BBHs that will merge before undergoing further interactions. This interestingly suggests that the BBH merger rate derived from LIGO can be used to forecast the number of nearby LISA sources, as well as providing insight into the fraction of BBH mergers forming in GCs.

The dynamical processing of black holes in the dense cores of globular clusters (GCs), makes them efficient factories for producing binary black holes (BBHs). Here we explore the population of BBHs that form dynamically in GCs and may be observable at mHz frequencies or higher with LISA. We use our Monte Carlo stellar dynamics code, which includes gravitational radiation reaction effects for all BH encounters. By creating a representative local universe of GCs, we show that up to dozens of these systems may be resolvable by LISA with signal-to-noise ratios of at least 5. Approximately one third of these binaries will have measurable eccentricities ($e > 10^{-3}$) in the LISA band and a small number ($\lesssim 5$) may evolve from the LISA band to the LIGO band during the LISA mission.

We present the first systematic study of strong binary-single and binary-binary black hole interactions with the inclusion of general relativity. When including general relativistic effects in strong encounters, dissipation of orbital energy from gravitational waves (GWs) can lead to captures and subsequent inspirals with appreciable eccentricities when entering the sensitive frequency ranges of the LIGO and Virgo GW detectors. In this study, we perform binary-binary and binary-single scattering experiments with general relativistic dynamics up through the 2.5 post-Newtonian order included, both in a controlled setting to gauge the importance of non-dissipative post-Newtonian terms and derive scaling relations for the cross-section of GW captures, as well as experiments tuned to the strong interactions from state-of-the art globular cluster models to assess the relative importance of the binary-binary channel at facilitating GW captures and the resultant eccentricity distributions of inspiral from channel. Although binary-binary interactions are 10-100 times less frequent in globular clusters than binary-single interactions, their longer lifetime and more complex dynamics leads to a higher probability for GW captures to occur during the encounter. We find that binary-binary interactions contribute 25-45% of the eccentric mergers which occur during strong black hole encounters in globular clusters, regardless of the properties of the cluster environment. The inclusion of higher multiplicity encounters in dense star clusters therefore have major implications on the predicted rates of highly eccentric binaries potentially detectable by the LIGO/Virgo network. As gravitational waveforms of eccentric inspirals are distinct from those generated by merging binaries which have circularized, measurements of eccentricity in such systems would highly constrain their formation scenario.

The dynamical assembly of binary black holes (BBHs) in dense star clusters (SCs) is one of the most promising pathways for producing observable gravitational wave (GW) sources, however several other formation scenarios likely operate as well. One of the current outstanding questions is how these different pathways may be distinguished apart. In this paper we suggest a new multi-messenger observable that can be used to constrain the formation of BBH mergers originating from SCs: the electromagnetic signal from tidal disruptions (TDs) of stars by BBHs. Such TDs will show variability in their light curve from the orbital motion of the disruptive BBHs, and can therefore be used to map the BBH orbital period distribution, and thereby also the dynamical mechanisms that eventually drive the BBHs to merger. Using an analytical approach including General Relativistic effects, we find that the orbital period distribution of BBHs within globular clusters peaks on timescales of days, which we argue is unique to this assembly pathway. We propose that the search for variable TDs in current and future EM transient surveys might be used to constrain the merger history of BBHs in SCs.

Binary black holes (BBHs) appear to be widespread and are able to merge through the emission of gravitational waves, as recently illustrated by LIGO. The spin of the BBHs is one of the parameters that LIGO can infer from the gravitational wave signal and can be used to constrain their production site. If BBHs are assembled in stellar clusters they are likely to interact with stars, which could occasionally lead to a tidal disruption event (TDE). When a BBH tidally disrupts a star it can accrete a significant fraction of the debris, effectively altering the spins of the BHs. Therefore, although dynamically formed BBHs are expected to have random spin orientations, tidal stellar interactions can significantly alter their birth spins both in direction and magnitude. Here we investigate how TDEs by BBHs can affect the properties of the BH members as well as exploring the characteristics of the resulting electromagnetic signatures. We conduct hydrodynamic simulations with a Lagrangian Smoothed Particle Hydrodynamics code of a wide range of representative tidal interactions. We find that both spin magnitude and orientation can be altered and temporarily aligned or anti-aligned through accretion of stellar debris, with a significant dependence on the mass ratio of the disrupted star and the BBH members. These tidal interactions feed material to the BBH at very high accretion rates, with the potential to launch a relativistic jet. The corresponding beamed emission is a beacon to an otherwise quiescent BBH.

Using state-of-the-art dynamical simulations of globular clusters, including radiation reaction during black hole encounters and a cosmological model of star cluster formation, we create a realistic population of dynamically-formed binary black hole mergers across cosmic space and time. We show that in the local universe, 10% of these binaries form as the result of gravitational-wave emission between unbound black holes during chaotic resonant encounters, with roughly half of those events having eccentricities detectable by current ground-based gravitational-wave detectors. The mergers that occur inside clusters typically have lower masses than binaries that were ejected from the cluster many Gyrs ago. Gravitational-wave captures from globular clusters contribute 1-2 Gpc^-3 yr^-1 to the binary merger rate in the local universe, increasing to ~10 Gpc^-3 yr^-1 at z~3. Finally, we discuss some of the technical difficulties associated with post-Newtonian scattering encounters, and how care must be taken when measuring the binary parameters during a dynamical capture.

In paper I of this series we showed that a large percentage of the binary black hole (BBH) mergers that form through dynamical interactions in globular clusters will have significant eccentricity in the ~10^{-3}-10^{-1} Hz LISA band. In this work we quantify the evolution of these highly eccentric binaries through the LISA and LIGO bands, and compute the stochastic gravitational wave background from the merging, eccentric population. We find that the population of BBHs that merge in-between three-body encounters inside their cluster (~50% of all cluster-formed BBH mergers) will have measurable eccentricity for their entire lifetime in the LISA band. The population of BBHs that merge during three-body encounters (~5% of all cluster-formed BBH mergers), will be detectable by LIGO with eccentricities of e~0.1. The gravitational wave background from dynamically assembled BBHs encodes a characteristic bump due to the high initial eccentricities of these systems. The location and amplitude of this bump depends on globular cluster properties.

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.

In this paper we study the dynamical outcome in which black hole (BH) binary-single interactions lead to two successive gravitational wave (GW) mergers; a scenario we refer to as a `double GW merger'. The first GW merger happens during the three-body interaction through a two-body GW capture, where the second GW merger is between the BH formed in the first GW merger and the remaining bound single BH. We estimate the probability for observing both GW mergers, and for observing only the second merger that we refer to as a `prompt second-generation (2G) merger'. We find that the probability for observing both GW mergers is only notable for co-planar interactions with low GW kicks ($\lesssim 10^{1}-10^{2}$ kms$^{-1}$), which suggests that double GW mergers can be used to probe environments facilitating such interactions. For isotropic encounters, such as the one found in globular clusters, the probability for prompt 2G mergers to form is only at the percent level, suggesting that second-generation mergers are most likely to be between BHs which have swapped partners at least once.

We show that nearly half of all binary black hole (BBH) mergers dynamically assembled in globular clusters have measurable eccentricities ($e>0.01$) in the LISA band ($10^{-2}$ Hz), when General Relativistic corrections are properly included in the $N$-body evolution. If only Newtonian gravity is included, the derived fraction of eccentric LISA sources is significantly lower, which explains why recent studies all have greatly underestimated this fraction. Our findings have major implications for how to observationally distinguish between BBH formation channels using eccentricity with LISA, which is one of the key science goals of the mission. We illustrate that the relatively large population of eccentric LISA sources reported here originates from BBHs that merge between hardening binary-single interactions inside their globular cluster. These results indicate a bright future for using LISA to probe the origin of BBH mergers.

We derive the probability for a newly formed binary black hole (BBH) to undergo an eccentric gravitational wave (GW) merger during binary-single interactions inside a stellar cluster. By integrating over the hardening interactions such a BBH must undergo before ejection, we find that the observable rate of BBH mergers with eccentricity $>0.1$ at $10\ \rm{Hz}$ relative to the rate of circular mergers can be as high as $\sim 5\%$ for a typical globular cluster (GC). This further suggests that BBH mergers forming through GW captures in binary-single interactions, eccentric or not, are likely to constitute $\sim 10\%$ of the total BBH merger rate from GCs. Such GW capture mergers can only be probed with an $N$-body code that includes General Relativistic corrections, which explains why recent Newtonian cluster studies not have been able to resolve this population. Finally, we show that the relative rate of eccentric BBH mergers depends on the compactness of their host cluster, suggesting that an observed eccentricity distribution can be used to probe the origin of BBH mergers.

We study the gravitational wave (GW) frequency and chirp mass distribution of binary black hole (BBH) mergers assembled through three-body interactions in globular clusters (GCs), when GW emission at the 2.5 post-Newtonian (PN) level is included in the $N$-body equation-of-motion (EOM). From performing $\sim 2.5\times10^{6}$ PN binary-single interactions based on GC data from the `MOCCA-Survey Database I' project, and by the use of analytical methods, we find that $5-10\%$ of all the three-body assembled GC BBH mergers have a GW frequency at formation that is $\gtrsim 10^{-1}$ Hz, implying they enter the LIGO band without having drifted through the LISA band first. If PN terms are not included in the EOM, one finds instead that all BBH mergers drifts through both LISA and LIGO. As the fraction of BBH mergers that only show up in LIGO is expected to be $\sim 0\%$ for standard field binary BBH mergers, future joint measurements with LISA and LIGO can be used to gain insight into the formation of BBH mergers.

We estimate the population of eccentric gravitational wave (GW) binary black hole (BBH) mergers forming during binary-single interactions in globular clusters (GCs), using ~ 800 GC models that were evolved using the MOCCA code for star cluster simulations as part of the MOCCA-Survey Database I project. By re-simulating binary-single interactions (only involving 3 BHs) extracted from this set of GC models using an N-body code that includes GW emission at the 2.5 post-Newtonian level, we find that ~ 10% of all the BBHs assembled in our GC models that merge at present time form during chaotic binary-single interactions, and that about half of this sample have an eccentricity > 0.1 at 10 Hz. We explicitly show that this derived rate of eccentric mergers is ~ 100 times higher than one would find with a purely Newtonian N-body code. Furthermore, we demonstrate that the eccentric fraction can be accurately estimated using a simple analytical formalism when the interacting BHs are of similar mass; a result that serves as the first successful analytical description of eccentric GW mergers forming during three-body interactions in realistic GCs.

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 present a study on how the outcomes of co-planar binary-single black hole (BH) interactions distribute as a function of the orbital initial conditions. We refer to this distribution as the topology. Using an $N$-body code that includes BH finite sizes and gravitational wave (GW) emission in the equation-of-motion (EOM), we perform more than a million three-body scatterings to explore the topology of both the classical and the relativistic limit. We further describe how the inclusion of GW emission in the EOM naturally leads to scenarios where the binary-single system undergoes two successive GW mergers. The time span from the first to the second GW merger is shortest in co-planar interactions, and are for some configurations even short enough for LIGO to observe both events. We note that co-planar interactions could be frequent in environments such as active galactic nuclei discs.

We present a study on binary-single interactions with energy loss terms such as tidal dissipation and gravitational wave emission added to the equation-of-motion. The inclusion of such terms leads to the formation of compact binaries that form during the three-body interaction through two-body captures. These binaries predominantly merge relative promptly at high eccentricity, with several observable and dynamical consequences to follow. Despite their possibility for being observed in both present and upcoming transient surveys, their outcomes are not firmly constrained. In this paper we present an analytical framework that allows to estimate the cross section of such two-body captures, which permits us to study how the corresponding rates depends on the initial orbital parameters, the mass hierarchy, the type of interacting objects, and the energy dissipation mechanism. This formalism is applied here to study the formation of two-body gravitational wave captures, for which we estimate absolute and relative rates relevant to Advanced LIGO detections. It is shown that two-body gravitational wave captures should have compelling observational implications if a sizable fraction of detected compact binaries are formed via dynamical interactions.

In this {\it Letter} we calculate the fraction of highly eccentric binary black hole (BBH) mergers resulting from binary-single interactions. Using an $N$-body code that includes post-Newtonian correction terms, we show that $\gtrsim 1\%$ of all BBH mergers resulting from this channel will have an eccentricity $e>0.1$ when coming into the LIGO frequency band. As the majority of BBH mergers forming in globular clusters are assembled through three-body encounters, we suggest that such interactions are likely to dominate the population of high eccentricity BBH mergers detectable by LIGO. The relative frequency of highly eccentric events could eventually help to constrain the astrophysical origin of BBH mergers.

We perform the first systematic study on how dynamical stellar tides and general relativistic (GR) effects affect the dynamics and outcomes of binary-single interactions. For this, we have constructed an N-body code that includes tides in the affine approximation, where stars are modeled as self-similar ellipsoidal polytropes, and GR corrections using the commonly-used post-Newtonian formalism. Using this numerical formalism, we are able resolve the leading effect from tides and GR across several orders of magnitude in both stellar radius and initial target binary separation. We find that the main effect from tides is the formation of two-body tidal captures that form during the chaotic and resonant evolution of the triple system. The two stars undergoing the capture spiral in and merge. The inclusion of tides can thus lead to an increase on the stellar coalescence rate. We also develop an analytical framework for calculating the cross section of tidal inspirals between any pair of objects with similar mass. From our analytical and numerical estimates we find that the rate of tidal inspirals relative to collisions increases as the initial semi-major axis of the target binary increases and the radius of the interacting tidal objects decreases. The largest effect is therefore found for triple systems hosting white dwarfs and neutron stars. In this case, we find the rate of highly eccentric white dwarf - neutron star mergers to likely be dominated by tidal inspirals. While tidal inspirals occur rarely, we note that they can give rise to a plethora of thermonuclear transients such as Ca-rich transients.

We present a simple and powerful method for extracting transit signals associated with a known transiting planet from noisy light curves. Assuming the orbital period of the planet is known and the signal is periodic, we illustrate that systematic noise can be removed in Fourier space at all frequencies, by only using data within a fixed time frame with a width equal to an integer number of orbital periods. This results in a reconstruction of the full transit signal which on average is unbiased, despite that no prior knowledge of either the noise or the transit signal itself is used in the analysis. The method has therefore clear advantages over standard phase folding, which normally requires external input such as nearby stars or noise models for removing systematic components. In addition, we can extract the full orbital transit signal ($360$ degrees) simultaneously, and \emph{Kepler} like data can be analyzed in just a few seconds. We illustrate the performance of our method by applying it to a dataset composed of light curves from \emph{Kepler} with a fake injected signal emulating a planet with rings. For extracting periodic transit signals, our presented method is in general the optimal and least biased estimator and could therefore lead the way toward the first detections of, e.g., planet rings and exo-trojan asteroids.

We present a dynamical model describing how halo particles can receive a significant energy kick from the merger between their own host halo and a target halo. This is highly relevant for understanding the growth of cosmological halos, and could especially provide an explanation for some high velocity objects. The model we present includes a \emph{double scattering mechanism}, where a halo particle is given a significant energy kick by undergoing two subsequent gravitational deflections during the merger. The first deflection is by the potential of the target halo, whereas the second is by the potential of the particle's original host halo. The resultant energy kick arises because the two halos move relative to each other during the two deflections. To our knowledge, this mechanism has never been characterized in this context before. We derive analytically a halo particle's total kick energy, which is composed of energy from the double scattering mechanism and energy release from tidal fields, as a function of its position in its original host halo. In the case of a $1:10$ merger between two Hernquist halos, we estimate that the presented mechanisms can generate particles with a velocity $\sim 2$ times the virial velocity of the target halo measured at its virial sphere. This motivates us to suggest that the high velocity of the recently discovered globular cluster HVCG-1 \citep{2014ApJ...787L..11C} can be explained by a head-on halo merger. Finally, we illustrate the orbital evolution of particles outside the virial sphere of the target halo, by solving the equation of motion in an expanding universe. We find a 'sweet spot' around a scale factor of $0.3-0.5$ for ejecting particles into large orbits, which easily can reach beyond $\sim 5$ virial radii.

White dwarfs (WDs) can be tidally disrupted only by massive black holes (MBHs) with masses less than $\sim10^5 M_\odot$. These tidal interactions feed material to the MBH well above its Eddington limit, with the potential to launch a relativistic jet. The corresponding beamed emission is a promising signpost to an otherwise quiescent MBH of relatively low mass. We show that the mass transfer history, and thus the lightcurve, are quite different when the disruptive orbit is parabolic, eccentric, or circular. The mass lost each orbit exponentiates in the eccentric-orbit case leading to the destruction of the WD after several tens of orbits. We examine the stellar dynamics of clusters surrounding MBHs to show that single-passage WD disruptions are substantially more common than repeating encounters. The $10^{49}$ erg s$^{-1}$ peak luminosity of these events makes them visible to cosmological distances. They may be detectible at rates of as many as tens per year by instruments like Swift. In fact, WD-disruption transients significantly outshine their main-sequence star counterparts, and are the most likely tidal interaction to be detected arising from MBHs with masses less than $10^5 M_\odot$. The detection or non-detection of such WD-disruption transients by Swift is, therefore, a powerful tool to constrain lower end of the MBH mass function. The emerging class of ultra-long gamma ray bursts all have peak luminosities and durations reminiscent of WD disruptions, offering a hint that WD-disruption transients may already be present in existing datasets.

The inspiral and merger of eccentric binaries leads to gravitational waveforms distinct from those generated by circularly merging binaries. Dynamical environments can assemble binaries with high eccentricity and peak frequencies within the {\it LIGO} band. In this paper, we study binary-single stellar scatterings occurring in dense stellar systems as a source of eccentrically-inspiraling binaries. Many interactions between compact binaries and single objects are characterized by chaotic resonances in which the binary-single system undergoes many exchanges before reaching a final state. During these chaotic resonances, a pair of objects has a non-negligible probability of experiencing a very close passage. Significant orbital energy and angular momentum are carried away from the system by gravitational wave (GW) radiation in these close passages and in some cases this implies an inspiral time shorter than the orbital period of the bound third body. We derive the cross section for such dynamical inspiral outcomes through analytical arguments and through numerical scattering experiments including GW losses. We show that the cross section for dynamical inspirals grows with increasing target binary semi-major axis, $a$, and that for equal-mass binaries it scales as $a^{2/7}$. Thus, we expect wide target binaries to predominantly contribute to the production of these relativistic outcomes. We estimate that eccentric inspirals account for approximately one percent of dynamically assembled non-eccentric merging binaries. While these events are rare, we show that binary-single scatterings are a more effective formation channel than single-single captures for the production of eccentrically-inspiraling binaries, even given modest binary fractions.

Cosmic microwave background measurements have pushed to higher resolution, lower noise, and more sky coverage. These data enable a unique test of the early universe's expansion rate and constituents such as effective number of relativistic degrees of freedom and dark energy. Using the most recent data from Planck and WMAP9, we constrain the expansion history in a model independent manner from today back to redshift z=10^5. The Hubble parameter is mapped to a few percent precision, limiting early dark energy and extra relativistic degrees of freedom within a model independent approach to 2-16% and 0.71 equivalent neutrino species respectively (95% CL). Within dark radiation, barotropic aether, and Doran-Robbers models, the early dark energy constraints are 3.3%, 1.9%, 1.2% respectively.

Redshift space distortions in galaxy clustering offer a promising technique for probing the growth rate of structure and testing dark energy properties and gravity. We consider the issue of to what accuracy they need to be modeled in order not to unduly bias cosmological conclusions. Fitting for nonlinear and redshift space corrections to the linear theory real space density power spectrum in bins in wavemode, we analyze both the effect of marginalizing over these corrections and of the bias due to not correcting them fully. While naively subpercent accuracy is required to avoid bias in the fixed case, in the fitting approach the Kwan-Lewis-Linder reconstruction function for redshift space distortions is found to be accurately selfcalibrated with little degradation in dark energy and gravity parameter estimation for a next generation galaxy redshift survey such as BigBOSS.

We examine model independent constraints on the high redshift and prerecombination expansion history from cosmic microwave background observations, using a combination of principal component analysis and other techniques. This can be translated to model independent limits on early dark energy and the number of relativistic species $N_{\rm eff}$. Models such as scaling (Doran-Robbers), dark radiation ($\Delta N_{\rm eff}$), and barotropic aether fall into distinct regions of eigenspace and can be easily distinguished from each other. Incoming CMB data will map the expansion history from $z=0$--$10^5$, achieving subpercent precision around recombination, and enable determination of the amount of early dark energy and valuable guidance to its nature.

Observations and numerical simulations of galaxy clusters strongly indicate that the hot intracluster x-ray emitting gas is not spherically symmetric. In many earlier studies spherical symmetry has been assumed partly because of limited data quality, however new deep observations and instrumental designs will make it possible to go beyond that assumption. Measuring the temperature and density profiles are of interest when observing the x-ray gas, however the spatial shape of the gas itself also carries very useful information. For example, it is believed that the x-ray gas shape in the inner parts of galaxy clusters is greatly affected by feedback mechanisms, cooling and rotation, and measuring this shape can therefore indirectly provide information on these mechanisms. In this paper we present a novel method to measure the three-dimensional shape of the intracluster x-ray emitting gas. We can measure the shape from the x-ray observations only, i.e. the method does not require combination with independent measurements of e.g. the cluster mass or density profile. This is possible when one uses the full spectral information contained in the observed spectra. We demonstrate the method by measuring radial dependent shapes along the line of sight for CHANDRA mock data. We find that at least 10^6 photons are required to get a 5-{\sigma} detection of shape for an x-ray gas having realistic features such as a cool core and a double powerlaw for the density profile. We illustrate how Bayes' theorem is used to find the best fitting model of the x-ray gas, an analysis that is very important in a real observational scenario where the true spatial shape is unknown. Not including a shape in the fit may propagate to a mass bias if the x-ray is used to estimate the total cluster mass. We discuss this mass bias for a class of spacial shapes.

The accuracy in the photometry of a point source depends on the point-spread function (PSF), detector pixelization, and observing strategy. The PSF and pixel response describe the spatial blurring of the source, the pixel scale describes the spatial sampling of a single exposure, and the observing strategy determines the set of dithered exposures with pointing offsets from which the source flux is inferred. In a wide-field imaging survey, sources of interest are randomly distributed within the field of view and hence are centered randomly within a pixel. A given hardware configuration and observing strategy therefore have a distribution of photometric uncertainty for sources of fixed flux that fall in the field. In this article we explore the ensemble behavior of photometric and position accuracies for different PSFs, pixel scales, and dithering patterns. We find that the average uncertainty in the flux determination depends slightly on dither strategy, whereas the position determination can be strongly dependent on the dithering. For cases with pixels much larger than the PSF, the uncertainty distributions can be non-Gaussian, with rms values that are particularly sensitive to the dither strategy. We also find that for these configurations with large pixels, pointings dithered by a fractional pixel amount do not always give minimal average uncertainties; this is in contrast to image reconstruction for which fractional dithers are optimal. When fractional pixel dithering is favored, a pointing accuracy of better than $\sim 0.15$ pixel width is required to maintain half the advantage over random dithers.

Some functions entering cosmological analysis, such as the dark energy equation of state or systematic uncertainties, are unknown functions of redshift. To include them without assuming a particular form we derive an efficient method for generating realizations of all possible functions subject to certain bounds or physical conditions, e.g. w\in[-1,+1] as for quintessence. The method is optimal in the sense that it is both pure and complete in filling the allowed space of principal components. The technique is applied to propagation of systematic uncertainties in supernova population drift and dust corrections and calibration through to cosmology parameter estimation and bias in the magnitude-redshift Hubble diagram. We identify specific ranges of redshift and wavelength bands where the greatest improvements in supernova systematics due to population evolution and dust correction can be achieved.

We present the first automated spectroscopic search for disk-galaxy lenses, using the Sloan Digital Sky Survey database. We follow up eight gravitational lens candidates, selected among a sample of ~40000 candidate massive disk galaxies, using a combination of ground-based imaging and long-slit spectroscopy. We confirm two gravitational lens systems: one probable disk galaxy, and one probable S0 galaxy. The remaining systems are four promising disk-galaxy lens candidates, as well as two probable gravitational lenses whose lens galaxy might be an S0 galaxy. The redshifts of the lenses are z_lens ~ 0.1. The redshift range of the background sources is z_source ~ 0.3 - 0.7. The systems presented here are (confirmed or candidate) galaxy-galaxy lensing systems, that is, systems where the multiple images are faint and extended, allowing an accurate determination of the lens galaxy mass and light distributions without contamination from the background galaxy. Moreover, the low redshift of the (confirmed or candidates) lens galaxies is favorable for measuring rotation points to complement the lensing study. We estimate the rest-frame total mass-to-light ratio within the Einstein radius for the two confirmed lenses: we find M_tot/L_I = 5.4 +- 1.5 within 3.9 +- 0.9 kpc for SDSS J081230.30+543650.9, and M_tot/L_I = 1.5 +- 0.9 within 1.4 +- 0.8 kpc for SDSS J145543.55+530441.2 (all in solar units). Hubble Space Telescope or Adaptive Optics imaging is needed to further study the systems.