In the presence of an ultralight bosonic field, spinning black holes are unstable to superradiance. The rotational energy of the black hole is converted into a non-axisymmetric, oscillating boson cloud which dissipates through the emission of nearly monochromatic gravitational radiation. Thus, gravitational wave observations by ground- or space-based detectors can be used to probe the existence of dark particles weakly coupled to the Standard Model. In this work, we focus on massive vector bosons, which grow much faster through superradiance, and produce significantly stronger gravitational waves compared to the scalar case. We use techniques from black hole perturbation theory to compute the relativistically-correct gravitational wave signal across the parameter space of different boson masses and black hole masses and spins. This fills in a gap in the literature between flatspace approximations, which underestimate the gravitational wave amplitude in the non-relativistic limit, and overestimate it in the relativistic regime, and time-domain calculations, which have only covered a limited part of the parameter space. We also identify parameter ranges where overtone superradiantly unstable modes will grow faster than the lower frequency fundamental modes. Such cases will produce a distinct gravitational wave signal due to the beating of the simultaneously populated modes, which we compute.

Spin can have significant effects on the electromagnetic transients accompanying binary neutron star mergers. The measurement of spin can provide important information about binary formation channels. In the absence of a strong neutron star spin prior, the degeneracy of spin with other parameters leads to significant uncertainties in their estimation, in particular limiting the power of gravitational waves to place tight constraints on the nuclear equation of state. Thus detailed studies of highly spinning neutron star mergers are essential to understand all aspects of multimessenger observation of such events. We perform a systematic investigation of the impact of neutron star spin---considering dimensionless spin values up to $a_{\rm NS}=0.33$---on the merger of equal mass, quasicircular binary neutron stars using fully general-relativistic simulations. We find that the peak frequency of the post-merger gravitational wave signal is only weakly influenced by the neutron star spin, with cases where the spin is aligned (antialigned) with the orbital angular momentum giving slightly lower (higher) values compared to the irrotational case. We find that the one-arm instability arises in a number of cases, with some dependence on spin. Spin has a pronounced impact on the mass, velocity, and angular distribution of the dynamical ejecta, and the mass of the disk that remains outside the merger remnant. We discuss the implications of these findings on anticipated electromagnetic signals, and on constraints that have been placed on the equation of state based on multimessenger observations of GW170817.

We study the dynamics and observational signatures of a neutron star being consumed by a much less massive black hole residing inside the star. This phenomenon could arise in a variety of scenarios, including after the capture of a primordial black hole, or in some models of asymmetric dark matter where the dark matter particles collect at the center of a neutron star and eventually collapse to form a black hole. However, the details of how the neutron star implodes are not well known, which is crucial to determining the observational implications of such events. We utilize fully general relativistic simulations to follow the evolution of such a black hole as it grows by several orders of magnitude, and ultimately consumes the neutron star. We consider a range of spin values for the neutron star, from non-rotating stars, to those with millisecond periods, as well as different equations of state. We find that as the black hole grows, it obtains a non-negligible spin and induces differential rotation in the core of the neutron star. In contrast to previous studies, we find that the amount of dynamical ejecta is very small, even for rapidly rotating stars, dampening the prospects for producing a kilonova-type electromagnetic signal from such events. We comment on other possible electromagnetic and gravitational signals.

We study the dynamics of small inhomogeneities in an expanding universe collapsing to form bound structures using full solutions of the Einstein-Vlasov (N-body) equations. We compare these to standard Newtonian N-body solutions using quantities defined with respect to fiducial observers in order to bound relativistic effects. We focus on simplified initial conditions containing a limited range of length scales, but vary the inhomogeneities from small magnitude, where the Newtonian and general-relativistic calculations agree quite well, to large magnitude, where the background metric receives an order one correction. For large inhomogeneities, we find that the collapse of overdensities tends to happen faster in Newtonian calculations relative to fully general-relativistic ones. Even in this extreme regime, the differences in the spacetime evolution outside the regions of large gravitational potential and velocity are small. For standard cosmological values, we corroborate the robustness of Newtonian N-body simulations to model large scale perturbations and the related cosmic variance in the local expansion rate.

Future gravitational-wave observations will enable unprecedented and unique science in extreme gravity and fundamental physics answering questions about the nature of dynamical spacetimes, the nature of dark matter and the nature of compact objects.

The physics of neutron stars can be studied with gravitational waves emitted from coalescing binary systems. Tidal effects become significant during the last few orbits and can be visible in the gravitational-wave spectrum above 500 Hz. After the merger, the neutron star remnant oscillates at frequencies above 1 kHz and can collapse into a black hole. Gravitational-wave detectors with a sensitivity of ~10^{-24} strain/sqHz at 2-4 kHz can observe these oscillations from a source which is ~100 Mpc away. The current observatories, such as LIGO and Virgo, are limited by shot noise at high frequencies and have a sensitivity of > 2 * 10^{-23} strain/sqHz at 3 kHz. In this paper, we propose an optical configuration of gravitational-wave detectors which can be set up in present facilities using the current interferometer topology. This scheme has a potential to reach 7 * 10^{-25} strain/sqHz at 2.5 kHz without compromising the detector sensitivity to black hole binaries. We argue that the proposed instruments have a potential to detect similar amount of post-merger neutron star oscillations as the next generation detectors, such as Cosmic Explorer and Einstein Telescope. We also optimise the arm length of the future detectors for neutron star physics and find that the optimal arm length is ~20 km. These instruments have the potential to observe neutron star post-merger oscillations at a rate of ~30 events per year with a signal-to-noise ratio of 5 or more.

The detection of gravitational waves from coalescing binary neutron stars represents another milestone in gravitational-wave astronomy. However, since LIGO is currently not as sensitive to the merger/ringdown part of the waveform, the possibility that such signals are produced by a black hole-neutron star binary can not be easily ruled out without appealing to assumptions about the underlying compact object populations. We review a few astrophysical channels that might produce black holes below 3 $M_{\odot}$ (roughly the upper bound on the maximum mass of a neutron star), as well as existing constraints for these channels. We show that, due to the uncertainty in the neutron star equation of state, it is difficult to distinguish gravitational waves from a binary neutron star system, from those of a black hole-neutron star system with the same component masses, assuming Advanced LIGO sensitivity. This degeneracy can be broken by accumulating statistics from many events to better constrain the equation of state, or by third-generation detectors with higher sensitivity to the late spiral to post-merger signal. We also discuss the possible differences in electromagnetic counterparts between binary neutron star and low mass black hole-neutron star mergers, arguing that it will be challenging to definitively distinguish the two without better understanding of the underlying astrophysical processes.

We study the superradiant instability of a massive boson around a spinning black hole in full general relativity without assuming spatial symmetries. We focus on the case of a rapidly spinning black hole in the presence of a vector boson with a Compton wavelength comparable to the black hole radius, which is the regime where relativistic effects are maximized. We follow the growth of the boson cloud through superradiance into the nonlinear regime as it spins down the black hole, reaches a maximum energy, and begins to dissipate through the emission of gravitational waves. We find that the superradiant instability can efficiently convert a significant fraction of a black hole's rotational energy into gravitational radiation.

We revisit the problem of a spinning black hole immersed in a uniformly magnetized plasma within the context of force-free electrodynamics. Such configurations have been found to relax to stationary jetlike solutions that are powered by the rotational energy of the black hole. We write down an analytic description for the jet solutions in the low black hole spin limit, and demonstrate that it provides a good approximation to the configurations found dynamically. For characterizing the magnetospheres of rapidly spinning black holes, we find that the current sheet which forms in the black hole ergosphere plays an essential role. We study the properties of the current sheet, and its importance in determining the jet solution and the rate at which energy is extracted from the black hole.

This work is the first in a series of studies aimed at understanding the dynamics of highly eccentric binary neutron stars, and constructing an appropriate gravitational-waveform model for detection. Such binaries are possible sources for ground-based gravitational wave detectors, and are expected to form through dynamical scattering and multi-body interactions in globular clusters and galactic nuclei. In contrast to black holes, oscillations of neutron stars are generically excited by tidal effects after close pericenter passage. Depending on the equation of state, this can enhance the loss of orbital energy by up to tens of percent over that radiated away by gravitational waves during an orbit. Under the same interaction mechanism, part of the orbital angular momentum is also transferred to the star. We calculate the impact of the neutron star oscillations on the orbital evolution of such systems, and compare these results to full numerical simulations. Utilizing a Post-Newtonian flux description we propose a preliminary model to predict the timing of different pericenter passages. A refined version of this model (taking into account Post-Newtonian corrections to the tidal coupling and the oscillations of the stars) may serve as a waveform model for such highly eccentric systems.

In the standard approach to studying cosmological structure formation, the overall expansion of the Universe is assumed to be homogeneous, with the gravitational effect of inhomogeneities encoded entirely in a Newtonian potential. A topic of ongoing debate is to what degree this fully captures the dynamics dictated by general relativity, especially in the era of precision cosmology. To quantitatively assess this, we directly compare standard N-body Newtonian calculations to full numerical solutions of the Einstein equations, for cold matter with various magnitude initial inhomogeneities on scales comparable to the Hubble horizon. We analyze the differences in the evolution of density, luminosity distance, and other quantities defined with respect to fiducial observers. This is carried out by reconstructing the effective spacetime and matter fields dictated by the Newtonian quantities, and by taking care to distinguish effects of numerical resolution. We find that the fully general relativistic and Newtonian calculations show excellent agreement, even well into the nonlinear regime. They only notably differ in regions where the weak gravity assumption breaks down, which arise when considering extreme cases with perturbations exceeding standard values.

We study the superradiant instability of massive vector fields, i.e. Proca fields, around spinning black holes in the test field limit. This is motivated by the possibility that observations of astrophysical black holes can probe the existence of ultralight bosons subject to this mechanism. By making use of time-domain simulations, we characterize the growth rate, frequency, spatial distribution, and other properties of the unstable modes, including in the regime where the black hole is rapidly spinning and the Compton wavelength of the Proca field is comparable to the black hole radius. We find that relativistic effects in this regime increase the range of Proca masses that are unstable, as well as the maximum instability rate. We also study the gravitational waves that can be sourced by such an instability, finding that they can be significantly stronger than in the massive scalar field case.

We study the growth and saturation of the superradiant instability of a complex, massive vector (Proca) field as it extracts energy and angular momentum from a spinning black hole, using numerical solutions of the full Einstein-Proca equations. We concentrate on a rapidly spinning black hole ($a=0.99$) and the dominant $m=1$ azimuthal mode of the Proca field, with real and imaginary components of the field chosen to yield an axisymmetric stress-energy tensor and, hence, spacetime. We find that in excess of $9\%$ of the black hole's mass can be transferred into the field. In all cases studied, the superradiant instability smoothly saturates when the black hole's horizon frequency decreases to match the frequency of the Proca cloud that spontaneously forms around the black hole.

A remarkable prediction of the Standard Model is that, in the absence of corrections lifting the energy density, the Higgs potential becomes negative at large field values. If the Higgs field samples this part of the potential during inflation, the negative energy density may locally destabilize the spacetime. We use numerical simulations of the Einstein equations to study the evolution of inflation-induced Higgs fluctuations as they grow towards the true (negative-energy) minimum. These simulations show that forming a single patch of true vacuum in our past light cone during inflation is incompatible with the existence of our Universe; the boundary of the true vacuum region grows outward in a causally disconnected manner from the crunching interior, which forms a black hole. We also find that these black hole horizons may be arbitrarily elongated---even forming black strings---in violation of the hoop conjecture. By extending the numerical solution of the Fokker-Planck equation to the exponentially suppressed tails of the field distribution at large field values, we derive a rigorous correlation between a future measurement of the tensor-to-scalar ratio and the scale at which the Higgs potential must receive stabilizing corrections in order for the Universe to have survived inflation until today.

We continue our investigations of the development and importance of the one-arm spiral instability in long-lived hypermassive neutron stars (HMNSs) formed in dynamical capture binary neutron star mergers. Employing hydrodynamic simulations in full general relativity, we find that the one-arm instability is generic in that it can develop in HMNSs within a few tens of milliseconds after merger for all equations of state in our survey. We find that mergers with stiffer equations of state tend to produce HMNSs with stronger $m=2$ azimuthal mode density deformations, and weaker $m=1$ components, relative to softer equations of state. We also find that for equations of state that can give rise to double-core HMNSs, large $m=1$ density modes can already be present due to asymmetries in the two cores. This results in the generation of $l=2$, $m=1$ gravitational wave modes even before the dominance of a one-arm mode that ultimately arises following merger of the two cores. Our results further suggest that stiffer equations of state give rise to HMNSs generating lower $m=1$ gravitational wave frequencies. Thus, if gravitational waves from the one-arm instability are detected, they could in principle constrain the neutron star equation of state. We estimate that, depending on the equation of state, the one-arm mode could potentially be detectable by second generation gravitational wave detectors at $\sim 10$ Mpc and by third generation ones at $\sim 100$ Mpc. Finally, we provide estimates of the properties of dynamical ejecta, as well as the accompanying kilonovae signatures.

Using numerical solutions of the full Einstein field equations coupled to a scalar inflaton field in 3+1 dimensions, we study the conditions under which a universe that is initially expanding, highly inhomogeneous and dominated by gradient energy can transition to an inflationary period. If the initial scalar field variations are contained within a sufficiently flat region of the inflaton potential, and the universe is spatially flat or open on average, inflation will occur following the dilution of the gradient and kinetic energy due to expansion. This is the case even when the scale of the inhomogeneities is comparable to the initial Hubble length, and overdense regions collapse and form black holes, because underdense regions continue expanding, allowing inflation to eventually begin. This establishes that inflation can arise from highly inhomogeneous initial conditions and solve the horizon and flatness problems, at least as long as the variations in the scalar field do not include values that exceed the inflationary plateau.

Many powerful and variable gamma-ray sources, including pulsar wind nebulae, active galactic nuclei and gamma-ray bursts, seem capable of accelerating particles to gamma-ray emitting energies efficiently over very short time scales. These are likely due to rapid dissipation of electromagnetic energy in a highly magnetized, relativistic plasma. In order to understand the generic features of such processes, we have investigated simple models based on relaxation of unstable force-free magnetostatic equilibria. In this work, we make the connection between the corresponding plasma dynamics and the expected radiation signal, using 2D particle-in-cell simulations that self-consistently include synchrotron radiation reaction. We focus on the lowest order unstable force-free equilibrium in a 2D periodic box. We find that rapid variability, with modest apparent radiation efficiency as perceived by a fixed observer, can be produced during the evolution of the instability. The "flares" are accompanied by an increased polarization degree in the high energy band, with rapid variation in the polarization angle. Furthermore, the separation between the acceleration sites and the synchrotron radiation sites for the highest energy particles facilitates acceleration beyond the synchrotron radiation reaction limit. We also discuss the dynamical consequences of radiation reaction, and some astrophysical applications of this model. Our current simulations with numerically tractable parameters are not yet able to reproduce the most dramatic gamma-ray flares, e.g., from Crab Nebula. Higher magnetization studies are promising and will be carried out in the future.

We present the results of particle-in-cell numerical pair plasma simulations of relativistic 2D magnetostatic equilibria known as the 'ABC' fields. In particular, we focus on the lowest-order unstable configuration consisting of two minima and two maxima of the magnetic vector potential. Breaking of the initial symmetry leads to exponential growth of the electric energy and to the formation of two current layers, which is consistent with the picture of 'X-point collapse' first described by Syrovatskii. Magnetic reconnection within the layers heats a fraction of particles to very high energies. After the saturation of the linear instability, the current layers are disrupted and the system evolves chaotically, diffusing the particle energies in a stochastic second-order Fermi process leading to the formation of power-law energy distributions. The power-law slopes harden with the increasing mean magnetization, but they are significantly softer than those produced in simulations initiated from Harris-type layers. The maximum particle energy is proportional to the mean magnetization, which is attributed partly to the increase of the effective electric field and partly to the increase of the acceleration time scale. We describe in detail the evolving structure of the dynamical current layers, and report on the conservation of magnetic helicity. These results can be applied to highly magnetized astrophysical environments, where ideal plasma instabilities trigger rapid magnetic dissipation with efficient particle acceleration and flares of high-energy radiation.

We perform general-relativistic hydrodynamical simulations of dynamical capture binary neutron star mergers, emphasizing the role played by the neutron star spin. Dynamical capture mergers may take place in globular clusters, as well as other dense stellar systems, where most neutron stars have large spins. We find significant variability in the merger outcome as a function of initial neutron star spin. For cases where the spin is aligned with the orbital angular momentum, the additional centrifugal support in the remnant hypermassive neutron star can prevent the prompt collapse to a black hole, while for antialigned cases the decreased total angular momentum can facilitate the collapse to a black hole. We show that even moderate spins can significantly increase the amount of ejected material, including the amount unbound with velocities greater than half the speed of light, leading to brighter electromagnetic signatures associated with kilonovae and interaction of the ejecta with the interstellar medium. Furthermore, we find that the initial neutron star spin can strongly affect the already rich phenomenology in the postmerger gravitational wave signatures that arise from the oscillation modes of the hypermassive neutron star. In several of our simulations, the resulting hypermassive neutron star develops the one-arm ($m=1$) spiral instability, the most pronounced cases being those with small but non-negligible neutron star spins. For long-lived hypermassive neutron stars, the presence of this instability leads to improved prospects for detecting these events through gravitational waves, and thus may give information about the neutron star equation of state.

Freely decaying relativistic force-free turbulence is studied for the first time. We initiate the magnetic field at a short wavelength and simulate its relaxation toward equilibrium on two and three dimensional periodic domains, in both helical and non-helical settings. Force-free turbulent relaxation is found to exhibit an inverse cascade in all settings, and in 3D to have a magnetic energy spectrum consistent with the Kolmogorov $5/3$ power law. 3D relaxations also obey the Taylor hypothesis; they settle promptly into the lowest energy configuration allowed by conservation of the total magnetic helicity. But in 2D, the relaxed state is a force-free equilibrium whose energy greatly exceeds the Taylor minimum, and which contains persistent force-free current layers and isolated flux tubes. We explain this behavior in terms of additional topological invariants that exist only in two dimensions, namely the helicity enclosed within each level surface of the magnetic potential function. The speed and completeness of turbulent magnetic free energy discharge could help account for rapidly variable gamma-ray emission from the Crab Nebula, gamma-ray bursts, blazars, and radio galaxies.

Jets associated with Active Galactic Nuclei (AGN) have been observed for almost a century, initially at optical and radio wavelengths. They are now widely accepted as "exhausts" produced electromagnetically by the central, spinning, massive black hole and its orbiting, accreting gas. Observations at X-ray and, especially, gamma-ray energies have transformed our understanding of how these jets evolve dynamically, accelerate electrons (and positrons) and radiate throughout the entire electromagnetic spectrum. Some new approaches to modeling the powerful and rapidly variable TeV emission observed from many blazars are sketched. Observations at the highest TeV energies, to which the High Altitude Water Cherenkov Gamma-Ray Observatory (HAWC) will contribute, promise crucial discrimination between rival models of AGN jets.

Using general-relativistic hydrodynamical simulations, we show that merging binary neutron stars can form hypermassive neutrons stars that undergo the one-arm spiral instability. We study the particular case of a dynamical capture merger where the stars have a small spin, as may arise in globular clusters, and focus on an equal-mass scenario where the spins are aligned with the orbital angular momentum. We find that this instability develops when post-merger fluid vortices lead to the generation of a toroidal remnant - a configuration whose maximum density occurs in a ring around the center-of-mass - with high vorticity along its rotation axis. The instability quickly saturates on a timescale of $\sim 10$ ms, with the $m=1$ azimuthal density multipole mode dominating over higher modes. The instability also leaves a characteristic imprint on the post-merger gravitational wave signal that could be detectable if the instability persists in long-lived remnants.

In order to understand the conditions which lead a highly magnetized, relativistic plasma to become unstable, and in such cases how the plasma evolves, we study a prototypical class of magnetostatic equilibria where the magnetic field satisfies $\nabla \times\mathbf B = \alpha \mathbf B$, where $\alpha$ is spatially uniform, on a periodic domain. Using numerical solutions we show that generic examples of such equilibria are unstable to ideal modes (including incompressible ones) which are marked by exponential growth in the linear phase. We characterize the unstable mode, showing how it can be understood in terms of merging magnetic and current structures, and explicitly demonstrate its instability using the energy principle. Following the nonlinear evolution of these solutions, we find that they rapidly develop regions with relativistic velocities and electric fields of comparable magnitude to the magnetic field, liberating magnetic energy on dynamical timescales and eventually settling into a configuration with the largest allowable wavelength. These properties make such solutions a promising setting for exploring the mechanisms behind extreme cosmic sources of gamma rays.

We study dynamical capture binary black hole-neutron star (BH-NS) mergers focusing on the effects of the neutron star spin. These events may arise in dense stellar regions, such as globular clusters, where the majority of neutron stars are expected to be rapidly rotating. We initialize the BH-NS systems with positions and velocities corresponding to marginally unbound Newtonian orbits, and evolve them using general-relativistic hydrodynamical simulations. We find that even moderate spins can significantly increase the amount of mass in unbound material. In some of the more extreme cases, there can be up to a third of a solar mass in unbound matter. Similarly, large amounts of tidally stripped material can remain bound and eventually accrete onto the BH---as much as a tenth of a solar mass in some cases. These simulations demonstrate that it is important to treat neutron star spin in order to make reliable predictions of the gravitational wave and electromagnetic transient signals accompanying these sources.

We use simulations of hydrodynamics coupled with full general relativity to investigate the gravitational waves produced by a star colliding with a massive black hole when the star's tidal disruption radius lies far outside of the black hole horizon. We consider both main-sequence and white-dwarf compaction stars, and nonspinning black holes, as well as those with near-extremal spin. We study the regime in between where the star can be accurately modeled by a point particle, and where tidal effects completely suppress the gravitational wave signal. We find that nonnegligible gravitational waves can be produced even when the star is strongly affected by tidal forces, as well as when it collides with large angular momentum. We discuss the implications that these results have for the potential observation of gravitational waves from these sources with future detectors.

We study the superradiant scattering of gravitational waves by a nearly extremal black hole (dimensionless spin $a=0.99$) by numerically solving the full Einstein field equations, thus including backreaction effects. This allows us to study the dynamics of the black hole as it loses energy and angular momentum during the scattering process. To explore the nonlinear phase of the interaction, we consider gravitational wave packets with initial energies up to $10%$ of the mass of the black hole. We find that as the incident wave energy increases, the amplification of the scattered waves, as well as the energy extraction efficiency from the black hole, is reduced. During the interaction the apparent horizon geometry undergoes sizable nonaxisymmetric oscillations. The largest amplitude excitations occur when the peak frequency of the incident wave packet is above where superradiance occurs, but close to the dominant quasinormal mode frequency of the black hole.

We introduce a new method for numerically evolving the full Einstein field equations in situations where the spacetime is dominated by a known background solution. The technique leverages the knowledge of the background solution to subtract off its contribution to the truncation error, thereby more efficiently achieving a desired level of accuracy. We demonstrate the method by applying it to the radial infall of a solar-type star into supermassive black holes with mass ratios $\geq 10^6$. The self-gravity of the star is thus consistently modeled within the context of general relativity, and the star's interaction with the black hole computed with moderate computational cost, despite the over five orders of magnitude difference in gravitational potential (as defined by the ratio of mass to radius). We compute the tidal deformation of the star during infall, and the gravitational wave emission, finding the latter is close to the prediction of the point-particle limit.

We study the ultrarelativistic head-on collision of equal mass particles, modeled as self-gravitating fluid spheres, by numerically solving the coupled Einstein-hydrodynamic equations. We focus on cases well within the kinetic energy dominated regime, where between 88-92% ($\gamma=8$ to 12) of the initial net energy of the spacetime resides in the translation kinetic energy of the particles. We find that for sufficiently large boosts, black hole formation occurs. Moreover, near yet above the threshold of black hole formation, the collision initially leads to the formation of two distinct apparent horizons that subsequently merge. We argue that this can be understood in terms of a focusing effect, where one boosted particle acts as a gravitational lens on the other and vice versa, and that this is further responsible for the threshold being lower (by a factor of a few) compared to simple hoop conjecture estimates. Cases slightly below threshold result in complete disruption of the model particles. The gravitational radiation emitted when black holes form reaches luminosities of 0.014 $c^5/G$, carrying $16\pm2%$ of the total energy.

We assess the detectability of the gravitational wave signals from highly eccentric compact binaries. We use a simple model for the inspiral, merger, and ringdown of these systems. The model is based on mapping the binary to an effective single black hole system described by a Kerr metric, thereby including certain relativistic effects such as zoom-whirl-type behavior. The resultant geodesics source quadrupolar radiation and, in turn, are evolved under its dissipative effects. At the light ring, we attach a merger model that was previously developed for quasicircular mergers but also performs well for eccentric mergers with little modification. We apply this model to determine the detectability of these sources for initial, Enhanced, and Advanced LIGO across the parameter space of nonspinning close capture compact binaries. We conclude that, should these systems exist in nature, the vast majority will be missed by conventional burst searches or by quasicircular waveform templates in the advanced detector era. Other methods, such as eccentric templates or, more practically, a stacked excess power search, must be developed to avoid losing these sources. These systems would also have been missed frequently in the initial LIGO data analysis. Thus, previous null coincidence results with detected gamma-ray bursts cannot exclude the possibility of coincident gravitational wave signals from eccentric binaries.

We study dynamical capture binary neutron star mergers as may arise in dense stellar regions such as globular clusters. Using general-relativistic hydrodynamics, we find that these mergers can result in the prompt collapse to a black hole or in the formation of a hypermassive neutron star, depending not only on the neutron star equation of state but also on impact parameter. We also find that these mergers can produce accretion disks of up to a tenth of a solar mass and unbound ejected material of up to a few percent of a solar mass. We comment on the gravitational radiation and electromagnetic transients that these sources may produce.

There is great interest in numerical relativity simulations involving matter due to the likelihood that binary compact objects involving neutron stars will be detected by gravitational wave observatories in the coming years, as well as to the possibility that binary compact object mergers could explain short-duration gamma-ray bursts. We present a code designed for simulations of hydrodynamics coupled to the Einstein field equations targeted toward such applications. This code has recently been used to study eccentric mergers of black hole-neutron star binaries. We evolve the fluid conservatively using high-resolution shock-capturing methods, while the field equations are solved in the generalized-harmonic formulation with finite differences. In order to resolve the various scales that may arise, we use adaptive mesh refinement (AMR) with grid hierarchies based on truncation error estimates. A noteworthy feature of this code is the implementation of the flux correction algorithm of Berger and Colella to ensure that the conservative nature of fluid advection is respected across AMR boundaries. We present various tests to compare the performance of different limiters and flux calculation methods, as well as to demonstrate the utility of AMR flux corrections.

There is a high level of interest in black hole-neutron star binaries, not only because their mergers may be detected by gravitational wave observatories in the coming years, but also because of the possibility that they could explain a class of short duration gamma-ray bursts. We study black hole-neutron star mergers that occur with high eccentricity as may arise from dynamical capture in dense stellar regions such as nuclear or globular clusters. We perform general relativistic simulations of binaries with a range of impact parameters, three different initial black hole spins (zero, aligned and anti-aligned with the orbital angular momentum), and neutron stars with three different equations of state. We find a rich diversity across these parameters in the resulting gravitational wave signals and matter dynamics, which should also be reflected in the consequent electromagnetic emission. Before tidal disruption, the gravitational wave emission is significantly larger than perturbative predictions suggest for periapsis distances close to effective innermost stable separations, exhibiting features reflecting the zoom-whirl dynamics of the orbit there. Guided by the simulations, we develop a simple model for the change in orbital parameters of the binary during close encounters. Depending upon the initial parameters of the system, we find that mass transfer during non-merging close encounters can range from essentially zero to a sizable fraction of the initial neutron star mass. The same holds for the amount of material outside the black hole post-merger, and in some cases roughly half of this material is estimated to be unbound. We also see that non-merging close encounters generically excite large oscillations in the neutron star that are qualitatively consistent with f-modes.

Within the next few years gravitational waves (GWs) from merging black holes (BHs) and neutron stars (NSs) may be directly detected, making a thorough theoretical understanding of these systems a high priority. As an additional motivation, these systems may represent a subset of short-duration gamma-ray burst (sGRB) progenitors. BH-NS mergers are expected to result from primordial, quasi-circular inspiral as well as dynamically formed capture binaries. The latter channel allows mergers with high eccentricity, resulting in a richer variety of outcomes. We perform general relativistic simulations of BH-NS interactions with a range of impact parameters, and find significant variation in the properties of these events that have potentially observable consequences, namely the GW signature, remnant accretion disk mass, and amount of unbound material.

In this paper we investigate the feasibility of bremsstrahlung radiation from `nonthermal' electrons as a source of hard X-rays from the intracluster medium of clusters of galaxies. With an exact treatment of the Coulomb collisions in a Fokker-Planck analysis of the electron distribution we find that the severe difficulties with lifetimes of `nonthermal' particles found earlier by Petrosian (2001) using a cold target model remain problematic. We then address possible acceleration of background electrons into a nonthermal tail. We assume a simplified but generic acceleration rate and determine the expected evolution of an initially Maxwellian distribution of electrons. We find that strong nonthermal components arise only for rapid rate of acceleration which also heats up the entire plasma. These results confirm the conclusion that if the observed `nonthermal' excesses are due to some process accelerating the background thermal electrons this process must be short lived.