Theoretical modeling of massive stars predicts a gap in the black hole (BH) mass function above $\sim 40-50\,M_{\odot}$ for BHs formed through single star evolution, arising from (pulsational) pair-instability supernovae. However, in dense star clusters, dynamical channels may exist that allow construction of BHs with masses in excess of those allowed from single star evolution. The detection of BHs in this so-called "upper-mass gap" would provide strong evidence for the dynamical processing of BHs prior to their eventual merger. Here, we explore in detail the formation of BHs with masses within or above the pair-instability gap through pre-collapse stellar collisions in dense star clusters. We run a suite of 68 independent cluster simulations, exploring a variety of physical assumptions pertaining to growth through stellar collisions, including primordial cluster mass segregation and the efficiency of envelope stripping during collisions. We find that as many as $\sim20\%$ of all BH progenitors undergo one or more collisions prior to stellar collapse and up to $\sim1\%$ of all BHs reside within or above the pair-instability gap through the effects of these collisions. We show that these BHs readily go on to merge with other BHs in the cluster, creating a population of massive BH mergers at a rate that may compete with the "multiple-generation" merger channel described in other analyses. Finally, we describe how stellar collisions in clusters may provide a unique pathway to pair-instability supernovae and briefly discuss the expected rate of these events and other electromagnetic transients.

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 explore the possibility that GW190412, a binary black hole merger with a non-equal-mass ratio and significantly spinning primary, was formed through repeated black hole mergers in a dense super star cluster. Using a combination of semi-analytic prescriptions for the remnant spin and recoil kick of black hole mergers, we show that the mass ratio and spin of GW190412 are consistent with a binary black hole whose primary component has undergone two successive mergers from a population of $\sim 10M_{\odot}$ black holes in a high-metallicity environment. We then explore the production of GW190412-like analogs in the CMC Cluster Catalog, a grid of 148 $N$-body star cluster models, as well as a new model, behemoth, with nearly $10^7$ particles and initial conditions taken from a cosmological MHD simulation of galaxy formation. We show that the production of binaries with GW190412-like masses and spins is dominated by massive super star clusters with high metallicities and large central escape speeds. While many are observed in the local universe, our results suggest that a careful treatment of these massive clusters, many of which may have been disrupted before the present day, is necessary to characterize the production of unique gravitational-wave events produced through dynamics.

We present a set of 148 independent $N$-body simulations of globular clusters (GCs) computed using the code $\texttt{CMC}$ ($\texttt{Cluster Monte Carlo}$). At an age of $\sim10-13\,$Gyr, the resulting models cover nearly the full range of cluster properties exhibited by the Milky Way GCs, including total mass, core and half-light radii, metallicity, and galactocentric distance. We use our models to investigate the role that stellar-mass black holes play in the process of core collapse. Furthermore, we study how dynamical interactions affect the formation and evolution of several important types of sources in GCs, including low-mass X-ray binaries, millisecond pulsars, blue stragglers, cataclysmic variables, Type Ia supernovae, calcium-rich transients, and merging compact binaries. While our focus here is on old, low-metallicity GCs, our $\texttt{CMC}$ simulations follow the evolution of clusters over a Hubble time, and they include a wide range of metallicities (up to solar), so that our results can also be used to study younger and higher-metallicity star clusters.

The hierarchical three-body problem has many applications in relativistic astrophysics, and can potentially explain all of the binary black hole mergers detected by LIGO/Virgo. However, the majority of studies have only included relativistic corrections to the two-body equations of motion, and have ignored the relativistic effects that arise in the presence of a third body. We revisit this problem and develop a fully consistent derivation of the secular three-body problem to first post-Newtonian order. We start with the Einstein-Infeld-Hoffman equations for a three-body system and expand the accelerations as a power series in the ratio of the semi-major axes of the inner ($a_1$) and outer ($a_2$) binary. We then perform a post-Keplerian, two-parameter expansion of the single orbit-averaged Lagrange planetary equations in $\delta = v^2/c^2$ and $\epsilon = a_1/a_2$ using the method of multiple scales. It is well established that eccentricity growth through the Lidov-Kozai (LK) mechanism can be suppressed or amplified by two-body 1pN (2BpN) precession. In this study, we derive and investigate three-body 1pN (3BpN) effects. We find that these effects can lead to eccentricity growth in triples, occurring over hundreds of LK cycles. The octupole terms can enhance these effects and lead to even greater eccentricities. In such cases, inclusion of these effects can substantially alter the evolution of three-body systems as compared to an analysis in which they are neglected. Careful analysis of post-Newtonian three-body effects will be important to understand the formation and properties of coalescing binaries that form via three-body dynamical processes.

The first detection of gravitational waves from a neutron star - neutron star (NS-NS) merger, GW170817, and the increasing number of observations of short gamma-ray bursts (SGRBs) have greatly motivated studies of the origins of NS-NS and neutron star - black hole (NS-BH) binaries. We calculate the merger rates of NS-NS and NS-BH binaries from globular clusters (GCs) using realistic GC simulations with the \texttt{CMC} cluster catalog. We use a large sample of models with a range of initial numbers of stars, metallicities, virial radii and galactocentric distances, representative of the present-day Milky Way GCs, to quantify the inspiral times and volumetric merger rates as a function of redshift, both inside and ejected from clusters. We find that over the complete lifetime of most GCs, stellar BHs dominate the cluster cores and prevent the mass segregation of NSs, thereby reducing the dynamical interaction rates of NSs so that at most a few NS binary mergers are ever produced. We estimate the merger rate in the local universe to be $\sim\rm{0.02\,Gpc^{-3}\,yr^{-1}}$ for both NS-NS and NS-BH binaries, or a total of $\sim 0.04$~Gpc$^{-3}$~yr$^{-1}$ for both populations. These rates are about 5 orders of magnitude below the current empirical merger rate from LIGO/Virgo. We conclude that dynamical interactions in GCs do not play a significant role in enhancing the NS-NS and NS-BH merger rates.

The formation and evolution of binary stars is a critical component of several fields in astronomy. The most numerous sources for gravitational wave observatories are inspiraling and/or merging compact binaries, while binary stars are present in nearly every electromagnetic survey regardless of the target population. Simulations of large binary populations serve to both predict and inform observations of electromagnetic and gravitational wave sources. Binary population synthesis is a tool that balances physical modeling with simulation speed to produce large binary populations on timescales of days. We present a community-developed binary population synthesis suite: COSMIC which is designed to simulate compact-object binary populations and their progenitors. As a proof of concept, we simulate the Galactic population of compact binaries and their gravitational wave signal observable by the Laser Interferometer Space Antenna (LISA). We find that $\sim10^8$ compact binaries reside in the Milky Way today, while $\sim10^4$ of them may be resolvable by LISA.

A measurement of the history of cosmic star formation is central to understand the origin and evolution of galaxies. The measurement is extremely challenging using electromagnetic radiation: significant modeling is required to convert luminosity to mass, and to properly account for dust attenuation, for example. Here we show how detections of gravitational waves from inspiraling binary black holes made by proposed third-generation detectors can be used to measure the star formation rate of massive stars with high precision up to redshifts of ~10. Depending on the time-delay model, the predicted detection rates ranges from ~1400 to ~16000 per month with the current measurement of local merger rate density. With three months of observations, parameters describing the volumetric star formation rate can be constrained at the few percent level, and the volumetric merger rate can be directly measured to 3% at z~2. Given a parameterized star formation rate, the characteristic delay time between binary formation and merger can be measured to ~60%.

When two black holes merge in a dense star cluster, they form a new black hole with a well-defined mass and spin. If that "second-generation" black hole remains in the cluster, it will continue to participate in dynamical encounters, form binaries, and potentially merge again. Using a grid of 96 dynamical models of dense star clusters and a cosmological model of cluster formation, we explore the production of binary black hole mergers where at least one component of the binary was forged in a previous merger. We create four hypothetical universes where every black hole born in the collapse of a massive star has a dimensionless Kerr spin parameter, $\chi_{\rm birth}$, of 0.0, 0.1, 0.2, or 0.5. We show that if all stellar-born black holes are non-spinning ($\chi_{\rm birth} = 0.0$), then more than 10% of merging binary black holes from clusters have components formed from previous mergers, accounting for more than 20% of the mergers from globular clusters detectable by LIGO/Virgo. Furthermore, nearly 7% of detectable mergers would have a component with a mass $\gtrsim 55M_{\odot}$, placing them clearly in the mass "gap" region where black holes cannot form from isolated collapsing stars due to the pulsational-pair instability mechanism. On the other hand, if black holes are born spinning, then the contribution from these second-generation mergers decreases, making up as little as 1% of all detections from globular clusters when $\chi_{\rm birth} = 0.5$. We make quantitative predictions for the detected masses, mass ratios, and spin properties of first- and second-generation mergers from dense star clusters, and show how these distributions are highly sensitive to the birth spins of black holes.

As self-gravitating systems, dense star clusters exhibit a natural diffusion of energy from their innermost to outermost regions, which leads to a slow and steady contraction of the core until it ultimately collapses under gravity. However, in spite of the natural tendency toward so-called "core collapse," the globular clusters (GCs) in the Milky Way exhibit a well-observed bimodal distribution in core radii separating the core-collapsed and non-core-collapsed clusters. This suggests an internal energy source is at work, delaying the onset of core collapse in many clusters. Primordial binary stars have been thought for a long time to provide this energy source, but recent analyses have cast doubt upon the corresponding "binary-burning" mechanism as a viable explanation. Over the past decade, a large amount of both observational and theoretical work has suggested that many stellar-mass black holes (BHs) are retained in typical clusters today and that they play a dynamically-significant role in these clusters throughout their entire lifetimes. Here we review our latest understanding of the formation and evolution of BH populations in GCs and demonstrate that, through their dynamical interaction with their host cluster, BHs can naturally explain the distinction between core-collapsed and non-core-collapsed clusters through a process we call "black hole burning."

The Wide Aperture Exoplanet Telescope (WAET) is a new ground-based optical telescope layout with an extremely asymmetric aperture, which results in new exoplanet imaging reach at very low cost. We suggest that hWAET, a 100x2 telescope, can be built for $150M in the 2020s, and >300m versions merit further R&D.

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

The Galactic Center (GC) is dominated by the gravity of a super-massive black hole (SMBH), Sagittarius A$^*$, and is suspected to contain a sizable population of binary stars. Such binaries form hierarchical triples with the SMBH, undergoing Eccentric Kozai-Lidov (EKL) evolution, which can lead to high eccentricity excitations for the binary companions' mutual orbit. This effect can lead to stellar collisions or Roche-lobe crossings, as well as orbital shrinking due to tidal dissipation. In this work we investigate the dynamical and stellar evolution of such binary systems, especially with regards to the binaries' post-main-sequence evolution. We find that the majority of binaries (~75%) is eventually separated into single stars, while the remaining binaries (~25%) undergo phases of common-envelope evolution and/or stellar mergers. These objects can produce a number of different exotic outcomes, including rejuvenated stars, G2-like infrared-excess objects, stripped giant stars, Type Ia supernovae (SNe), cataclysmic variables (CVs), symbiotic binaries (SBs), or compact object binaries. We estimate that, within a sphere of 250 Mpc radius, about 7.5 to 15 Type Ia SNe per year should occur in galactic nuclei due to this mechanism, potentially detectable by ZTF and ASAS-SN. Likewise we estimate that, within a sphere of 1 Gpc$^3$ volume, about 10 to 20 compact object binaries form per year that could become gravitational wave sources. Based on results of EKL-driven compact object binary mergers in galactic nuclei by Hoang at al. (2018), this compact object binary formation rate translates to about 15 to 30 events per year detectable by Advanced LIGO.

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

We present new measurements of the flux power-spectrum P(k) of the $z<0.5$ HI Lyman-$\alpha$ forest spanning scales k ~ 0.001-0.1 s/km. These results were derived from 65 far ultraviolet quasar spectra (resolution R~18000) observed with the Cosmic Origin Spectrograph (COS) on board the Hubble Space Telescope. The analysis required careful masking of all contaminating, coincident absorption from HI and metal-line transitions of the Galactic interstellar medium and intervening absorbers as well as proper treatment of the complex COS line-spread function. From the P(k) measurements, we estimate the HI photoionization rate ($\Gamma_{\rm HI}$) in the z<0.5 intergalactic medium. Our results confirm most of the previous $\Gamma_{\rm HI}$ estimates. We conclude that previous concerns of a photon underproduction crisis are now resolved by demonstrating that the measured $\Gamma_{\rm HI}$ can be accounted for by ultraviolet emission from quasars alone. In a companion paper, we will present constraints on the thermal state of the $z<0.5$ intergalactic medium from the P(k) measurements presented here.

As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations.

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.

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

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.

The "gravitational million-body problem," to model the dynamical evolution of a self-gravitating, collisional N-body system with ~10^6 particles over many relaxation times, remains a major challenge in computational astrophysics. Unfortunately, current techniques to model such systems suffer from severe limitations. A direct N-body simulation with more than 10^5 particles can require months or even years to complete, while an orbit-sampling Monte Carlo approach cannot adequately model the dynamics in a dense cluster core, particularly in the presence of many black holes. We have developed a new technique combining the precision of a direct N-body integration with the speed of a Monte Carlo approach. Our Rapid And Precisely Integrated Dynamics code, the RAPID code, statistically models interactions between neighboring stars and stellar binaries while integrating directly the orbits of stars or black holes in the cluster core. This allows us to accurately simulate the dynamics of the black holes in a realistic globular cluster environment without the burdensome N^2 scaling of a full N-body integration. We compare RAPID models of idealized globular clusters to identical models from the direct N-body and Monte Carlo methods. Our tests show that RAPID can reproduce the half-mass radii, core radii, black hole ejection rates, and binary properties of the direct N-body models far more accurately than a standard Monte Carlo integration while remaining significantly faster than a full N-body integration. With this technique, it will be possible to create more realistic models of Milky Way globular clusters with sufficient rapidity to explore the full parameter space of dense stellar clusters.

As the sensitivity of current and future gravitational-wave detectors improves, it will become possible to measure the evolution of the binary black hole merger rate with redshift. Here, we combine detailed fits to state-of-the-art dynamical models of binary black hole formation in dense star clusters with a cosmological model of cluster formation across cosmic time. We find a typical merger rate of 14 $\rm{Gpc}^{-3} \rm{yr}^{-1}$ in the local universe, with a reasonable range of 4-18 $\rm{Gpc}^{-3} \rm{yr}^{-1}$, depending on the rate of cluster disruption and the cluster initial mass function. This rate increases by a factor of 6 to redshift $z=2.7$ before declining at higher redshifts. We compare the merger rate from binaries produced in clusters to similar estimates from isolated binaries and triples in galactic fields, and discuss various ways that these different formation channels could add up to the current merger rate observed by LIGO/Virgo.

We explore the masses, merger rates, eccentricities, and spins for field binary black holes driven to merger by a third companion through the Lidov-Kozai mechanism. Using a population synthesis approach, we model the creation of stellar-mass black hole triples across a range of different initial conditions and stellar metallicities. We find that the production of triple-mediated mergers is enhanced at low metallicities by a factor of ~100 due to the lower black hole natal kicks and reduced stellar mass loss. These triples naturally yield heavy binary black holes with near-zero effective spins, consistent with most of the mergers observed to date. This process produces a merger rate of between 2 and 25 Gpc^-3 yr^-1 in the local universe, suggesting that the Lidov-Kozai mechanism can potentially explain all of the low-spin, heavy black hole mergers observed by Advanced LIGO/Virgo. Finally, we show that triples admit a unique eccentricity and spin distribution that will allow this model to be tested in the near future.

Recent discoveries of black hole (BH) candidates in Galactic and extragalactic globular clusters (GCs) have ignited interest in understanding how BHs dynamically evolve in a GC and the number of BHs ($N_{\rm{BH}}$) that may still be retained by today's GCs. Numerical models show that even if stellar-mass BHs are retained in today's GCs, they are typically in configurations that are not directly detectable. We show that a suitably defined measure of mass segregation ($\Delta$) between, e.g., giants and low-mass main-sequence stars, can be an effective probe to indirectly estimate $N_{\rm{BH}}$ in a GC aided by calibrations from numerical models. Using numerical models including all relevant physics we first show that $N_{\rm{BH}}$ is strongly anticorrelated with $\Delta$ between giant stars and low-mass main-sequence stars. We apply the distributions of $\Delta$ vs $N_{\rm{BH}}$ obtained from models to three Milky Way GCs to predict the $N_{\rm{BH}}$ retained by them at present. We calculate $\Delta$ using the publicly available ACS survey data for 47 Tuc, M 10, and M 22, all with identified stellar-mass BH candidates. Using these measured $\Delta$ and distributions of $\Delta$ vs $N_{\rm{BH}}$ from models as calibration we predict distributions for $N_{\rm{BH}}$ expected to be retained in these GCs. For 47 Tuc, M 10, and M 22 our predicted distributions peak at $N_{\rm{BH}}\approx20$, $24$, and $50$, whereas, within the $2\sigma$ confidence level, $N_{\rm{BH}}$ can be up to $\sim150$, $50$, and $200$, respectively.

The binary black hole mergers detected by Advanced LIGO/Virgo have shown no evidence of large black hole spins. However, because LIGO/Virgo best measures the effective combination of the two spins along the orbital angular momentum ($\chi_{\rm eff}$), it is difficult to distinguish between binaries with slowly-spinning black holes and binaries with spins lying in the orbital plane. Here, we study the spin dynamics for binaries with a distant black hole companion. For spins initially aligned with the orbital angular momentum of the binary, we find that $\chi_{\rm eff}$ "freezes" near zero as the orbit decays through the emission of gravitational waves. Through a population study, we show that this process predominantly leads to merging black hole binaries with near-zero $\chi_{\rm eff}$. We conclude that if the detected black hole binaries were formed in triples, then this would explain their low $\chi_{\rm eff}$ without the need to invoke near-zero spins or initially large spin-orbit angles.

We explore the formation of double-compact-object binaries in Milky Way (MW) globular clusters (GCs) that may be detectable by the Laser Interferometer Space Antenna (LISA). We use a set of 137 fully evolved GC models that, overall, effectively match the properties of the observed GCs in the MW. We estimate that, in total, the MW GCs contain $\sim21$ sources that will be detectable by LISA. These detectable sources contain all combinations of black hole (BH), neutron star, and white dwarf components. We predict $\sim7$ of these sources will be BH-BH binaries. Furthermore, we show that some of these BH-BH binaries can have signal-to-noise ratios large enough to be detectable at the distance of the Andromeda galaxy or even the Virgo cluster.

We present models of realistic globular clusters with post-Newtonian dynamics for black holes. By modeling the relativistic accelerations and gravitational-wave emission in isolated binaries and during three- and four-body encounters, we find that nearly half of all binary black hole mergers occur inside the cluster, with about 10% of those mergers entering the LIGO/Virgo band with eccentricities greater than 0.1. In-cluster mergers lead to the birth of a second generation of black holes with larger masses and high spins, which, depending on the black hole natal spins, can sometimes be retained in the cluster and merge again. As a result, globular clusters can produce merging binaries with detectable spins regardless of the birth spins of black holes formed from massive stars. These second-generation black holes would also populate any upper mass gap created by pair-instability supernovae.

Numerical simulations have shown that black holes (BHs) can strongly influence the evolution and present-day observational properties of globular clusters (GCs). Using a Monte Carlo code, we construct GC models that match the Milky Way (MW) cluster NGC 3201, the first cluster in which a stellar-mass BH was identified through radial-velocity measurements. We predict that NGC 3201 contains $\gtrsim 200$ stellar-mass BHs. Furthermore, we explore the dynamical formation of main sequence-BH binaries and demonstrate that systems similar to the observed BH binary in NGC 3201 are produced naturally. Additionally, our models predict the existence of bright blue-straggler-BH binaries unique to core-collapsed clusters, which otherwise retain few BHs.

We explore the population of mass-transferring binaries ejected from globular clusters (GCs) with both black hole (BH) and neutron star (NS) accretors. We use a set of 137 fully evolved globular cluster models which span a large range in cluster properties and, overall, match very well the properties of old GCs observed in the Milky Way. We identify all binaries ejected from our set of models that eventually undergo mass-transfer. These binaries are ejected from their host clusters over a wide range of ejection times and include white dwarf, giant, and main sequence donors. We calculate the orbits of these ejected systems in the Galactic potential to determine their present-day positions in the Galaxy and compare to the distribution of observed low-mass X-ray binaries (XRBs) in the Milky Way. We estimate $\sim 300$ mass-transferring NS binaries and $\sim 180$ mass-transferring BH binaries may currently be present in the Milky Way that originated from within GCs. Of these, we estimate, based on mass-transfer rates and duty cycles at the present time, at most a few would be observable as BH--XRBs and NS--XRBs at the present day. Based on our results, XRBs that originated from GCs are unlikely to contribute significantly to the total population of low-mass XRBs in the Galactic field.

The Wide Aperture Exoplanet Telescope (WAET) is a ground-based optical telescope layout in which one dimension of a filled aperture can be made very, very large (beyond 100 m) at low cost and complexity. With an unusual beam path but an otherwise-conventional optics, we obtain a fully-steerable telescope on a low-rise mount with a fixed gravity vector on key components. Numerous design considerations and scaling laws suggest that WAET can be far less expensive than other giant segmented mirror telescopes.

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

Gravitational waves (GWs) from binary black hole (BBH) mergers provide a new probe of massive-star evolution and the formation channels of binary compact objects. By coupling the growing sample of BBH systems with population synthesis models, we can begin to constrain the parameters of such models and glean unprecedented knowledge about the inherent physical processes that underpin binary stellar evolution. In this study, we apply a hierarchical Bayesian model to mass measurements from a synthetic GW sample to constrain the physical prescriptions in population models and the relative fraction of systems generated from various channels. We employ population models of two canonical formation scenarios in our analysis --- isolated binary evolution involving a common-envelope phase and dynamical formation within globular clusters --- with model variations for different black hole natal kick prescriptions. We show that solely with chirp mass measurements, it is possible to constrain natal kick prescriptions and the relative fraction of systems originating from each formation channel with $\mathcal{O}(100)$ of confident detections. This framework can be extended to include additional formation scenarios, model parameters, and measured properties of the compact binary.

Using numerical models for star clusters spanning a wide range in ages and metallicities (Z) we study the masses of binary black holes (BBHs) produced dynamically and merging in the local universe ($z\lesssim0.2$). After taking into account cosmological constraints on star-formation rate and metallicity evolution, which realistically relate merger delay times obtained from models with merger redshifts, we show here for the first time that while old, metal-poor globular clusters can naturally produce merging BBHs with heavier components, as observed in GW150914, lower-mass BBHs like GW151226 are easily formed dynamically in younger, higher-metallicity clusters.More specifically, we show that the mass of GW151226 is well within $1\sigma$ of the mass distribution obtained from our models for clusters with Z/Zsun$\gtrsim0.5$. Indeed dynamical formation of a system like GW151226 likely requires a cluster that is younger and has a higher metallicity than typical Galactic globular clusters. The LVT151012 system, if real, could have been created in any cluster with Z/Zsun$\lesssim0.25$. On the other hand, GW150914 is more massive (beyond $1\sigma$) than typical BBHs from even the lowest-metallicity (Z/Zsun=0.005) clusters we consider, but is within $2\sigma$ of the intrinsic mass distribution from our cluster models with Z/Zsun$\lesssim0.05$; of course detection biases also push the observed distributions towards higher masses.

Recent N-body simulations predict that large numbers of stellar black holes (BHs) could remain bound to globular clusters (GCs) at present, and merging BH--BH binaries are produced dynamically in significant numbers. We systematically vary "standard" assumptions made by numerical simulations related to, e.g., BH formation, stellar winds, binary properties of high-mass stars, and IMF within existing uncertainties, and study the effects on the evolution of the structural properties of GCs, and the BHs in GCs. We find that variations in initial assumptions can set otherwise identical initial clusters on completely different evolutionary paths, significantly affecting their present observable properties, or even affecting the cluster's very survival to the present. However, these changes usually do not affect the numbers or properties of local BH--BH mergers. The only exception is that variations in the assumed winds and IMF can change the masses and numbers of local BH--BH mergers, respectively. All other variations (e.g., in initial binary properties and binary fraction) leave the masses and numbers of locally merging BH--BH binaries largely unchanged. This is in contrast to binary population synthesis models for the field, where results are very sensitive to many uncertain parameters in the initial binary properties and binary stellar-evolution physics. We find that weak winds are required for producing GW150914-like mergers from GCs at low redshifts. LVT151012 can be produced in GCs modeled both with strong and weak winds. GW151226 is lower-mass than typical mergers from GCs modeled with weak winds, but is similar to mergers from GCs modeled with strong winds.

The recent detections of the binary black hole mergers GW150914 and GW151226 have inaugurated the field of gravitational-wave astronomy. For the two main formation channels that have been proposed for these sources, isolated binary evolution in galactic fields and dynamical formation in dense star clusters, the predicted masses and merger rates overlap significantly, complicating any astrophysical claims that rely on measured masses alone. Here, we examine the distribution of spin- orbit misalignments expected for binaries from the field and from dense star clusters. Under standard assumptions for black-hole natal kicks, we find that black-hole binaries similar to GW150914 could be formed with significant spin-orbit misalignment only through dynamical processes. In particular, these heavy-black-hole binaries can only form with a significant spin-orbit anti-alignment in the dynamical channel. Our results suggest that future detections of merging black hole binaries with measurable spins will allow us to identify the main formation channel for these systems.

The recent detections of GW150914 and GW151226 imply an abundance of stellar-mass binary-black-hole mergers in the local universe. While ground-based gravitational-wave detectors are limited to observing the final moments before a binary merges, space-based detectors, such as the Laser Interferometer Space Antenna (LISA), can observe binaries at lower orbital frequencies where such systems may still encode information about their formation histories. In particular, the orbital eccentricity and mass of binary black holes in the LISA frequency band can be used together to discriminate between binaries formed in isolation in galactic fields and those formed in dense stellar environments such as globular clusters. In this letter, we explore the orbital eccentricity and mass of binary-black-hole populations as they evolve through the LISA frequency band. Overall we find that there are two distinct populations discernible by LISA. We show that up to ~90% of binaries formed either dynamically or in isolation have eccentricities measurable by LISA. Finally, we note how measured eccentricities of low-mass binary black holes evolved in isolation could provide detailed constraints on the physics of black-hole natal kicks and common-envelope evolution.

We explore the possibility that GW150914, the binary black hole merger recently detected by Advanced LIGO, was formed by gravitational interactions in the core of a dense star cluster. Using models of globular clusters with detailed $N$-body dynamics and stellar evolution, we show that a typical cluster with a mass of $3\times10^5M_{\odot}$ to $6\times10^5M_{\odot}$ is optimal for forming GW150914-like binary black holes that will merge in the local universe. We identify the most likely dynamical processes for forming GW150914 in such a cluster, and we show that the detection of GW150914 is consistent with the masses and merger rates expected for binary black hole mergers from globular clusters. Our results show that dynamical processes provide a significant and well-understood pathway for forming binary black hole mergers in the local universe. Understanding the contribution of dynamics to the binary black hole merger problem is a critical step in unlocking the full potential of gravitational-wave astronomy.

The recent discovery of GW150914, the binary black hole merger detected by Advanced LIGO, has the potential to revolutionize observational astrophysics. But to fully utilize this new window into the universe, we must compare these new observations to detailed models of binary black hole formation throughout cosmic time. Expanding upon our previous work (Rodriguez et al., 2015), we study merging binary black holes formed in globular clusters using our Monte Carlo approach to stellar dynamics. We have created a new set of 52 cluster models with different masses, metallicities, and radii to fully characterize the binary black hole merger rate. These models include all the relevant dynamical processes (such as two-body relaxation, strong encounters, and three-body binary formation) and agree well with detailed direct N-body simulations. In addition, we have enhanced our stellar evolution algorithms with updated metallicity-dependent stellar wind and supernova prescriptions, allowing us to compare our results directly to the most recent population synthesis predictions for merger rates from isolated binary evolution. We explore the relationship between a cluster's global properties and the population of binary black holes that it produces. In particular, we derive a numerically calibrated relationship between the merger times of ejected black hole binaries and a cluster's mass and radius. We explore the masses and mass ratios of these binaries as a function of redshift, and find a merger rate of ~5 Gpc$^{-3}$ yr$^{-1}$ in the local universe, with 80% of sources having total masses from $32M_{\odot}$ to $64M_{\odot}$. Under standard assumptions, approximately 1 out of every 7 binary black hole mergers in the local universe will have originated in a globular cluster, but we also explore the sensitivity of this result to different assumptions for binary stellar evolution. [Abridged]

We present the first detailed comparison between million-body globular cluster simulations computed with a H\'enon-type Monte Carlo code, CMC, and a direct $N$-body code, NBODY6++GPU. Both simulations start from an identical cluster model with $10^6$ particles, and include all of the relevant physics needed to treat the system in a highly realistic way. With the two codes "frozen" (no fine-tuning of any free parameters or internal algorithms of the codes) we find excellent agreement in the overall evolution of the two models. Furthermore, we find that in both models, large numbers of stellar-mass black holes (> 1000) are retained for 12 Gyr. Thus, the very accurate direct $N$-body approach confirms recent predictions that black holes can be retained in present-day, old globular clusters. We find only minor disagreements between the two models and attribute these to the small-$N$ dynamics driving the evolution of the cluster core for which the Monte Carlo assumptions are less ideal. Based on the overwhelming general agreement between the two models computed using these vastly different techniques, we conclude that our Monte Carlo approach, which is more approximate, but dramatically faster compared to the direct $N$-body, is capable of producing a very accurate description of the long-term evolution of massive globular clusters even when the clusters contain large populations of stellar-mass black holes.

The predicted rate of binary black hole mergers from galactic fields can vary over several orders of magnitude and is extremely sensitive to the assumptions of stellar evolution. But in dense stellar environments such as globular clusters, binary black holes form by well-understood gravitational interactions. In this letter, we study the formation of black hole binaries in an extensive collection of realistic globular cluster models. By comparing these models to observed Milky Way and extragalactic globular clusters, we find that the mergers of dynamically-formed binaries could be detected at a rate of ~100 per year, potentially dominating the binary black hole merger rate. We also find that a majority of cluster-formed binaries are more massive than their field-formed counterparts, suggesting that Advanced LIGO could identify certain binaries as originating from dense stellar environments.

Hierarchical triple-star systems are expected to form frequently via close binary-binary encounters in the dense cores of globular clusters. In a sufficiently inclined triple, gravitational interactions between the inner and outer binary can cause large-amplitude oscillations in the eccentricity of the inner orbit ("Lidov-Kozai cycles"), which can lead to a collision and merger of the two inner components. In this paper we use Monte Carlo models of dense star clusters to identify all triple systems formed dynamically and we compute their evolution using a highly accurate three-body integrator which incorporates relativistic and tidal effects. We find that a large fraction of these triples evolve through a non-secular dynamical phase which can drive the inner binary to higher eccentricities than predicted by the standard secular perturbation theory (even including octupole-order terms). We place constraints on the importance of Lidov-Kozai-induced mergers for producing: (i) gravitational wave sources detectable by Advanced LIGO (aLIGO), for triples with an inner pair of stellar black holes; and (ii) blue straggler stars, for triples with main-sequence-star components. We find a realistic aLIGO detection rate of black hole mergers due to the Lidov-Kozai mechanism of 1yr^-1, with about 20% of these having a finite eccentricity when they first chirp into the aLIGO frequency band. While rare, these events are likely to dominate among eccentric compact object inspirals that are potentially detectable by aLIGO. For blue stragglers, we find that the Lidov-Kozai mechanism can contribute only up to ~10% of their total numbers in globular clusters.

The Advanced LIGO and Advanced Virgo gravitational wave (GW) detectors will begin operation in the coming years, with compact binary coalescence events a likely source for the first detections. The gravitational waveforms emitted directly encode information about the sources, including the masses and spins of the compact objects. Recovering the physical parameters of the sources from the GW observations is a key analysis task. This work describes the LALInference software library for Bayesian parameter estimation of compact binary signals, which builds on several previous methods to provide a well-tested toolkit which has already been used for several studies. We show that our implementation is able to correctly recover the parameters of compact binary signals from simulated data from the advanced GW detectors. We demonstrate this with a detailed comparison on three compact binary systems: a binary neutron star, a neutron star black hole binary and a binary black hole, where we show a cross-comparison of results obtained using three independent sampling algorithms. These systems were analysed with non-spinning, aligned spin and generic spin configurations respectively, showing that consistent results can be obtained even with the full 15-dimensional parameter space of the generic spin configurations. We also demonstrate statistically that the Bayesian credible intervals we recover correspond to frequentist confidence intervals under correct prior assumptions by analysing a set of 100 signals drawn from the prior. We discuss the computational cost of these algorithms, and describe the general and problem-specific sampling techniques we have used to improve the efficiency of sampling the compact binary coalescence parameter space.

Our current understanding of the stellar initial mass function and massive star evolution suggests that young globular clusters may have formed hundreds to thousands of stellar-mass black holes, the remnants of stars with initial masses from $\sim 20 - 100\, M_\odot$. Birth kicks from supernova explosions may eject some black holes from their birth clusters, but most should be retained. Using a Monte Carlo method we investigate the long-term dynamical evolution of globular clusters containing large numbers of stellar black holes. We describe numerical results for 42 models, covering a range of realistic initial conditions, including up to $1.6\times10^6$ stars. In almost all models we find that significant numbers of black holes (up to $\sim10^3$) are retained all the way to the present. This is in contrast to previous theoretical expectations that most black holes should be ejected dynamically within a few Gyr. The main reason for this difference is that core collapse driven by black holes (through the Spitzer "mass segregation instability") is easily reverted through three-body processes, and involves only a small number of the most massive black holes, while lower-mass black holes remain well-mixed with ordinary stars far from the central cusp. Thus the rapid segregation of stellar black holes does not lead to a long-term physical separation of most black holes into a dynamically decoupled inner core, as often assumed previously. Combined with the recent detections of several black hole X-ray binary candidates in Galactic globular clusters, our results suggest that stellar black holes could still be present in large numbers in many globular clusters today, and that they may play a significant role in shaping the long-term dynamical evolution and the present-day dynamical structure of many clusters.

Searches for a stochastic gravitational-wave background (SGWB) using terrestrial detectors typically involve cross-correlating data from pairs of detectors. The sensitivity of such cross-correlation analyses depends, among other things, on the separation between the two detectors: the smaller the separation, the better the sensitivity. Hence, a co-located detector pair is more sensitive to a gravitational-wave background than a non-co-located detector pair. However, co-located detectors are also expected to suffer from correlated noise from instrumental and environmental effects that could contaminate the measurement of the background. Hence, methods to identify and mitigate the effects of correlated noise are necessary to achieve the potential increase in sensitivity of co-located detectors. Here we report on the first SGWB analysis using the two LIGO Hanford detectors and address the complications arising from correlated environmental noise. We apply correlated noise identification and mitigation techniques to data taken by the two LIGO Hanford detectors, H1 and H2, during LIGO's fifth science run. At low frequencies, 40 - 460 Hz, we are unable to sufficiently mitigate the correlated noise to a level where we may confidently measure or bound the stochastic gravitational-wave signal. However, at high frequencies, 460-1000 Hz, these techniques are sufficient to set a $95%$ confidence level (C.L.) upper limit on the gravitational-wave energy density of \Omega(f)<7.7 x 10^{-4} (f/ 900 Hz)^3, which improves on the previous upper limit by a factor of $\sim 180$. In doing so, we demonstrate techniques that will be useful for future searches using advanced detectors, where correlated noise (e.g., from global magnetic fields) may affect even widely separated detectors.

In 2009-2010, the Laser Interferometer Gravitational-wave Observa- tory (LIGO) operated together with international partners Virgo and GEO600 as a network to search for gravitational waves of astrophysical origin. The sensitiv- ity of these detectors was limited by a combination of noise sources inherent to the instrumental design and its environment, often localized in time or frequency, that couple into the gravitational-wave readout. Here we review the performance of the LIGO instruments during this epoch, the work done to characterize the de- tectors and their data, and the effect that transient and continuous noise artefacts have on the sensitivity of LIGO to a variety of astrophysical sources.

We present the results of searches for gravitational waves from a large selection of pulsars using data from the most recent science runs (S6, VSR2 and VSR4) of the initial generation of interferometric gravitational wave detectors LIGO (Laser Interferometric Gravitational-wave Observatory) and Virgo. We do not see evidence for gravitational wave emission from any of the targeted sources but produce upper limits on the emission amplitude. We highlight the results from seven young pulsars with large spin-down luminosities. We reach within a factor of five of the canonical spin-down limit for all seven of these, whilst for the Crab and Vela pulsars we further surpass their spin-down limits. We present new or updated limits for 172 other pulsars (including both young and millisecond pulsars). Now that the detectors are undergoing major upgrades, and, for completeness, we bring together all of the most up-to-date results from all pulsars searched for during the operations of the first-generation LIGO, Virgo and GEO600 detectors. This gives a total of 195 pulsars including the most recent results described in this paper.

Cosmic strings can give rise to a large variety of interesting astrophysical phenomena. Among them, powerful bursts of gravitational waves (GWs) produced by cusps are a promising observational signature. In this Letter we present a search for GWs from cosmic string cusps in data collected by the LIGO and Virgo gravitational wave detectors between 2005 and 2010, with over 625 days of live time. We find no evidence of GW signals from cosmic strings. From this result, we derive new constraints on cosmic string parameters, which complement and improve existing limits from previous searches for a stochastic background of GWs from cosmic microwave background measurements and pulsar timing data. In particular, if the size of loops is given by the gravitational backreaction scale, we place upper limits on the string tension $G\mu$ below $10^{-8}$ in some regions of the cosmic string parameter space.

Within the next five years, it is expected that the Advanced LIGO/Virgo network will have reached a sensitivity sufficient to enable the routine detection of gravitational waves. Beyond the initial detection, the scientific promise of these instruments relies on the effectiveness of our physical parameter estimation capabilities. The majority of this effort has been towards the detection and characterization of gravitational waves from compact binary coalescence, e.g. the coalescence of binary neutron stars. While several previous studies have investigated the accuracy of parameter estimation with advanced detectors, the majority have relied on approximation techniques such as the Fisher Matrix. Here we report the statistical uncertainties that will be achievable for optimal detection candidates (SNR = 20) using the full parameter estimation machinery developed by the LIGO/Virgo Collaboration via Markov-Chain Monte Carlo methods. We find the recovery of the individual masses to be fractionally within 9% (15%) at the 68% (95%) credible intervals for equal-mass systems, and within 1.9% (3.7%) for unequal-mass systems. We also find that the Advanced LIGO/Virgo network will constrain the locations of binary neutron star mergers to a median uncertainty of 5.1 deg^2 (13.5 deg^2) on the sky. This region is improved to 2.3 deg^2 (6 deg^2) with the addition of the proposed LIGO India detector to the network. We also report the average uncertainties on the luminosity distances and orbital inclinations of ideal detection candidates that can be achieved by different network configurations.

Long gamma-ray bursts (GRBs) have been linked to extreme core-collapse supernovae from massive stars. Gravitational waves (GW) offer a probe of the physics behind long GRBs. We investigate models of long-lived (~10-1000s) GW emission associated with the accretion disk of a collapsed star or with its protoneutron star remnant. Using data from LIGO's fifth science run, and GRB triggers from the swift experiment, we perform a search for unmodeled long-lived GW transients. Finding no evidence of GW emission, we place 90% confidence level upper limits on the GW fluence at Earth from long GRBs for three waveforms inspired by a model of GWs from accretion disk instabilities. These limits range from F<3.5 ergs cm^-2 to $F<1200 ergs cm^-2, depending on the GRB and on the model, allowing us to probe optimistic scenarios of GW production out to distances as far as ~33 Mpc. Advanced detectors are expected to achieve strain sensitivities 10x better than initial LIGO, potentially allowing us to probe the engines of the nearest long GRBs.

Gravitational waves emitted during compact binary coalescences are a promising source for gravitational-wave detector networks. The accuracy with which the location of the source on the sky can be inferred from gravitational wave data is a limiting factor for several potential scientific goals of gravitational-wave astronomy, including multi-messenger observations. Various methods have been used to estimate the ability of a proposed network to localize sources. Here we compare two techniques for predicting the uncertainty of sky localization -- timing triangulation and the Fisher information matrix approximations -- with Bayesian inference on the full, coherent data set. We find that timing triangulation alone tends to over-estimate the uncertainty in sky localization by a median factor of $4$ for a set of signals from non-spinning compact object binaries ranging up to a total mass of $20 M_\odot$, and the over-estimation increases with the mass of the system. We find that average predictions can be brought to better agreement by the inclusion of phase consistency information in timing-triangulation techniques. However, even after corrections, these techniques can yield significantly different results to the full analysis on specific mock signals. Thus, while the approximate techniques may be useful in providing rapid, large scale estimates of network localization capability, the fully coherent Bayesian analysis gives more robust results for individual signals, particularly in the presence of detector noise.

Compact binary systems with neutron stars or black holes are one of the most promising sources for ground-based gravitational wave detectors. Gravitational radiation encodes rich information about source physics; thus parameter estimation and model selection are crucial analysis steps for any detection candidate events. Detailed models of the anticipated waveforms enable inference on several parameters, such as component masses, spins, sky location and distance that are essential for new astrophysical studies of these sources. However, accurate measurements of these parameters and discrimination of models describing the underlying physics are complicated by artifacts in the data, uncertainties in the waveform models and in the calibration of the detectors. Here we report such measurements on a selection of simulated signals added either in hardware or software to the data collected by the two LIGO instruments and the Virgo detector during their most recent joint science run, including a "blind injection" where the signal was not initially revealed to the collaboration. We exemplify the ability to extract information about the source physics on signals that cover the neutron star and black hole parameter space over the individual mass range 1 Msun - 25 Msun and the full range of spin parameters. The cases reported in this study provide a snap-shot of the status of parameter estimation in preparation for the operation of advanced detectors.

During the LIGO and Virgo joint science runs in 2009-2010, gravitational wave (GW) data from three interferometer detectors were analyzed within minutes to select GW candidate events and infer their apparent sky positions. Target coordinates were transmitted to several telescopes for follow-up observations aimed at the detection of an associated optical transient. Images were obtained for eight such GW candidates. We present the methods used to analyze the image data as well as the transient search results. No optical transient was identified with a convincing association with any of these candidates, and none of the GW triggers showed strong evidence for being astrophysical in nature. We compare the sensitivities of these observations to several model light curves from possible sources of interest, and discuss prospects for future joint GW-optical observations of this type.