Axions are some of the best motivated particles beyond the Standard Model. We show how the attractive self-interactions of dark matter (DM) axions over a broad range of masses, from $10^{-22}$ eV to $10^7$ GeV, can lead to nongravitational growth of density fluctuations and the formation of bound objects. This structure formation enhancement is driven by parametric resonance when the initial field misalignment is large, and it affects axion density perturbations on length scales of order the Hubble horizon when the axion field starts oscillating, deep inside the radiation-dominated era. This effect can turn an otherwise nearly scale-invariant spectrum of adiabatic perturbations into one that has a spike at the aforementioned scales, producing objects ranging from dense DM halos to scalar-field configurations such as solitons and oscillons. We call this class of cosmological scenarios for axion DM production "the large-misalignment mechanism." We explore observational consequences of this mechanism for axions with masses up to $10$ eV. For axions heavier than $10^{-5}$ eV, the compact axion halos are numerous enough to significantly impact Earth-bound direct detection experiments, yielding intermittent but coherent signals with repetition rates exceeding one per decade and crossing times less than a day. These episodic increases in the axion density and kinematic coherence suggest new approaches for axion DM searches, including for the QCD axion. Dense structures made up of axions from $10^{-22}$ eV to $10^{-5}$ eV are detectable through gravitational lensing searches, and their gravitational interactions can also perturb baryonic structures and alter star formation. At very high misalignment amplitudes, the axion field can undergo self-interaction-induced implosions long before matter-radiation equality, producing potentially-detectable low-frequency stochastic gravitational waves.

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.

The origin of many elements of the periodic table remains an unsolved problem. While many nucleosynthetic channels are broadly understood, significant uncertainties remain regarding certain groups of elements such as the intermediate and rapid neutron-capture processes, the p-process, or the origin of odd-Z elements in the most metal-poor stars. Canada has a long tradition of leadership in nuclear astrophysics, dating back to the work of Alastair Cameron in the 1950s. Recent faculty hires have further boosted activity in the field, including transient observation and theory, survey science on galactic nucleosynthesis, and nuclear experiments. This white paper contains a brief overview of recent activity in the community, highlighting strengths in each sub-field, and provides recommendations to improve interdisciplinary collaboration. Sustaining Canadian leadership in the next decade will require, on the observational side, access to transient and non-transient surveys like LSST, SKA, or MSE, support for target-of-opportunity observing in current and future Canadian telescopes, and participation in next-generation X-ray telescopes such as ATHENA. State-of-the-art theoretical predictions will require an ambitious succession plan for the Niagara supercomputer to support large parallel jobs. We propose a funding instrument for postdoctoral training that reflects the interdisciplinary nature of nuclear astrophysics research, and the creation of a national collaborative funding program that allows for joint projects and workshop organization.

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.

We explore the question of obtaining global solutions in Horndeski's theories of gravity. Towards this end, we study a relevant set of the theory and, by employing the Einstein frame we simplify the analysis by exploiting known results on global solutions of wave equations. We identify conditions for achieving global solutions as well as obstacles that can arise to spoil such goal. We illustrate such problems via numerical simulations.

The identification of the electromagnetic counterpart of the gravitational wave event, GW170817, and discovery of neutrinos and gamma-rays from TXS 0506+056 heralded the new era of multi-messenger astrophysics. As the number of multi-messenger events rapidly grow over the next decade, the cyberinfrastructure requirements to handle the increase in data rates, data volume, need for event follow up, and analysis across the different messengers will also explosively grow. The cyberinfrastructure requirements to enhance multi-messenger astrophysics will both be a major challenge and opportunity for astronomers, physicists, computer scientists and cyberinfrastructure specialists. Here we outline some of these requirements and argue for a distributed cyberinfrastructure institute for multi-messenger astrophysics to meet these challenges.

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.

Testing General Relativity and exploring possible departures has received further input with the possibility to do so through gravitational waves emitted in strongly gravitating/highly dynamical scenarios and also through the availability of exquisitely sensitive cosmological observations. However, most extensions suffer from severe pathologies at the mathematical level which have stymied a thorough exploration of putative theories. With the aid of a model problem which captures typical pathologies, we explore suggested methods to control them. We find that the approach that modifies the equations to control higher-order gradients is both robust and efficient.

The past year has witnessed discovery of the first identified counterparts to a gravitational wave transient (GW 170817A) and a very high-energy neutrino (IceCube-170922A). These source identifications, and ensuing detailed studies, have realized longstanding dreams of astronomers and physicists to routinely carry out observations of cosmic sources by other than electromagnetic means, and inaugurated the era of "multi-messenger" astronomy. While this new era promises extraordinary physical insights into the universe, it brings with it new challenges, including: highly heterogeneous, high-volume, high-velocity datasets; globe-spanning cross-disciplinary teams of researchers, regularly brought together into transient collaborations; an extraordinary breadth and depth of domain-specific knowledge and computing resources required to anticipate, model, and interpret observations; and the routine need for adaptive, distributed, rapid-response observing campaigns to fully exploit the scientific potential of each source. We argue, therefore, that the time is ripe for the community to conceive and propose an Institute for Multi-Messenger Astrophysics that would coordinate its resources in a sustained and strategic fashion to efficiently address these challenges, while simultaneously serving as a center for education and key supporting activities. In this fashion, we can prepare now to realize the bright future that we see, beyond, through these newly opened windows onto the universe.

Observations of gravitational radiation from compact binary systems provide an unprecedented opportunity to test General Relativity in the strong field dynamical regime. In this paper, we investigate how future observations of gravitational radiation from binary neutron star mergers might provide constraints on finite-range forces from a universally coupled massive scalar field. Such scalar degrees of freedom are a characteristic feature of many extensions of General Relativity. For concreteness, we work in the context of metric $f(R)$ gravity, which is equivalent to General Relativity and a universally coupled scalar field with a non-linear potential whose form is fixed by the choice of $f(R)$. In theories where neutron stars (or other compact objects) obtain a significant scalar charge, the resulting attractive finite-range scalar force has implications for both the inspiral and merger phases of binary systems. We first present an analysis of the inspiral dynamics in Newtonian limit, and forecast the constraints on the mass of the scalar and charge of the compact objects for the Advanced LIGO gravitational wave observatory. We then perform a comparative study of binary neutron star mergers in General Relativity with those of a one-parameter model of $f(R)$ gravity using fully relativistic hydrodynamical simulations. These simulations elucidate the effects of the scalar on the merger and post-merger dynamics. We comment on the utility of the full waveform (inspiral, merger, post-merger) to probe different regions of parameter space for both the particular model of $f(R)$ gravity studied here and for finite-range scalar forces more generally.

The ability of a plasma surrounding spinning black holes to extract rotational energy and power energetic emissions has been recognized as a key astrophysical phenomenon. Important insights into the nature of this process are obtained through the analysis of the interplay between a force-free magnetosphere and the black hole. This task involves solving a complicated system of equations, often requiring complex numerical simulations. Recent analytical attempts at tackling this problem have exploited the fact that the near horizon region of extreme Kerr (NHEK) is endowed with an enhanced symmetry group. We continue in this direction and show that for some conformally self-similar solutions, the NHEK force-free equations reduce to a single non-linear ordinary differential equation which is difficult to solve with straightforward integration. We here introduce a new approach specifically tailored to this problem and describe how one can obtain physically meaningful solutions.

Binary black hole mergers are among the most violent events in the Universe, leading to extreme warping of spacetime and copious emission of gravitational radiation. Even though black holes are the most compact objects they are not necessarily the most efficient emitters of gravitational radiation in binary systems. The final black hole resulting from a binary black hole merger retains a significant fraction of the pre-merger orbital energy and angular momentum. A non-vacuum system can in principle shed more of this energy than a black hole merger of equivalent mass. We study these super-emitters through a toy model that accounts for the possibility that the merger creates a compact object that retains a long-lived time-varying quadrupole moment. This toy model can capture the merger of neutron stars, but it can also be used to consider more exotic compact binaries. We hope that this toy model can serve as a guide to more rigorous numerical investigations into these systems.

We examine the development and detectability of the $m=1$ instability in the remnant of binary neutron star mergers. The detection of the gravitational mode associated with the $m=1$ degree of freedom could potentially reveal details of the equation of state. We analyze the post-merger epoch of simulations of both equal and non-equal mass neutron star mergers using three realistic, microphysical equations of state and neutrino cooling. Our studies show such an instability develops generically and within a short dynamical time to strengths that are comparable or stronger than the $m=2$ mode which is the strongest during the early post-merger stage. We estimate the signal to noise ratio that might be obtained for the $m=1$ mode and discuss the prospects for observing this signal with available Earth-based detectors. Because the $m=1$ occurs at roughly half the frequency of the more powerful $m=2$ signal and because it can potentially be long-lived, targeted searches could be devised to observe it. We estimate that with constant amplitude direct detection of the mode could occur up to a distance of roughly $14\,\mathrm{Mpc}$ whereas a search triggered by the inspiral signal could extend this distance to roughly $100\,\mathrm{Mpc}$.

We describe the full nonlinear development of the superradiant instability for a charged massless scalar field, coupled to general relativity and electromagnetism, in the vicinity of a Reissner-Nordstrom-AdS black hole. The presence of the negative cosmological constant provides a natural context for considering perfectly reflecting boundary conditions and studying the dynamics as the scalar field interacts repeatedly with the black hole. At early times, small superradiant perturbations grow as expected from linearized studies. Backreaction then causes the black hole to lose charge and mass until the perturbation becomes nonsuperradiant, with the final state described by a stable hairy black hole. For large gauge coupling, the instability extracts a large amount of charge per unit mass, resulting in greater entropy increase. We discuss the implications of the observed behavior for the general problem of superradiance in black hole spacetimes.

We study the merger of binary neutron stars with different mass ratios adopting three different realistic, microphysical nuclear equations of state, as well as incorporating neutrino cooling effects. In particular, we concentrate on the influence of the equation of state on the gravitational wave signature and also on its role, in combination with neutrino cooling, in determining the properties of the resulting hypermassive neutron star, of the neutrinos produced, and of the ejected material. The ejecta we find are consistent with other recent studies that find that small mass ratios produce more ejecta than equal mass cases (up to some limit) and this ejecta is more neutron rich. This trend indicates the importance with future kilonovae observations of measuring the individual masses of an associated binary neutron star system, presumably from concurrent gravitational wave observations, in order to be able to extract information about the nuclear equation of state

We derive exact scaling relations for two-dimensional relativistic hydrodynamic turbulence in the inertial range of scales. We consider both the energy cascade towards large scales and the enstrophy cascade towards small scales. We illustrate these relations by numerical simulations of turbulent weakly compressible flows. Intriguingly, the fluid-gravity correspondence implies that the gravitational field in spacetimes with anti-de Sitter asymptotics should exhibit similar scaling relations.

This work studies jet-like electromagnetic configurations surrounding a slowly-spinning black-hole immersed in a uniformly magnetized force-free plasma. In the first part of this work, we present a family of stationary solutions that are jet-capable. While these solutions all satisfy the force-free equations and the appropriate boundary conditions, our numerical experiments show a unique relaxed state starting from different initial data, and so one member of the family is likely preferred over the others. In the second part of this work, we analyze the perturbations of this family of jet-like solutions, and show that the perturbative modes exhibit a similar split into the trapped and traveling categories previously found for perturbed Blandford-Znajek solutions. In the eikonal limit, the trapped modes can be identified with the fast magnetosonic waves in the force-free plasma and the traveling waves are essentially the Alfven waves. Moreover, within the scope of our analysis, we have not seen signs of unstable modes at the light-crossing timescale of the system, within which the numerical relaxation process occurs. This observation disfavors mode instability as the selection mechanism for picking out a preferred solution. Consequently, our analytical study is unable to definitively select a particular solution out of the family to serve as the aforementioned preferred final state. This remains an interesting open problem.

We study the merger of binary neutron stars using different realistic, microphysical nuclear equations of state, as well as incorporating magnetic field and neutrino cooling effects. In particular, we concentrate on the influence of the equation of state on the gravitational wave signature and also on its role, in combination with cooling and electromagnetic effects, in determining the properties of the hypermassive neutron star resulting from the merger, the production of neutrinos, and the characteristics of ejecta from the system. The ejecta we find are consistent with other recent studies that find soft equations of state produce more ejecta than stiffer equations of state. Moreover, the degree of neutron richness increases for softer equations of state. In light of reported kilonova observations (associated to GRB~130603B and GRB~060614) and the discovery of relatively low abundances of heavy, radioactive elements in deep sea deposits (with respect to possible production via supernovae), we speculate that a soft EoS might be preferred---because of its significant production of sufficiently neutron rich ejecta---if such events are driven by binary neutron star mergers. We also find that realistic magnetic field strengths, obtained with a sub-grid model tuned to capture magnetic amplification via the Kelvin-Helmholtz instability at merger, are generally too weak to affect the gravitational wave signature post-merger within a time scale of $\approx 10$~ms but can have subtle effects on the post-merger dynamics.

As we showed in previous work, the dynamics and gravitational emission of binary neutron star systems in certain scalar-tensor theories can differ significantly from that expected from General Relativity in the coalescing stage. In this work we examine whether the characteristics of the electromagnetic counterparts to these binaries -- driven by magnetosphere interactions prior to the merger event -- can provide an independent way to test gravity in the most strongly dynamical stages of binary mergers. We find that the electromagnetic flux emitted by binaries in these scalar-tensor theories can show deviations from the GR prediction in particular cases. These differences are quite subtle, thus requiring delicate measurements to differentiate between GR and the type of scalar-tensor theories considered in this work using electromagnetic observations alone. However, if coupled with a gravitational-wave detection, electromagnetic measurements might provide a way to increase the confidence with which GR will be confirmed (or ruled out) by gravitational observations.

This article is an overview of the contributions numerical relativity has made to our understanding of strong field gravity, to be published in the book "General Relativity and Gravitation: A Centennial Perspective", commemorating the 100th anniversary of general relativity.

We study the dependence of the electromagnetic luminosity --produced by interactions of force-free magnetospheres-- on dipole inclinations in binary neutron star systems. We show that this interaction extracts kinetic energy from the system and powers a Poynting flux with a strong dependence on the dipole orientations. This dependence can be linked to the reconnection and redistribution of magnetic field as the stars interact. Although the details of the Poynting luminosity are very much dependent on the orientation, all the cases considered here nevertheless radiate a large Poynting flux. This robust emission suggests that the pre-merger stage of binary neutron star systems can yield interesting electromagnetic counterparts to gravitational wave events.

The theory of eternal inflation in an inflaton potential with multiple vacua predicts that our universe is one of many bubble universes nucleating and growing inside an ever-expanding false vacuum. The collision of our bubble with another could provide an important observational signature to test this scenario. We develop and implement an algorithm for accurately computing the cosmological observables arising from bubble collisions directly from the Lagrangian of a single scalar field. We first simulate the collision spacetime by solving Einstein's equations, starting from nucleation and ending at reheating. Taking advantage of the collision's hyperbolic symmetry, simulations are performed with a 1+1-dimensional fully relativistic code that uses adaptive mesh refinement. We then calculate the comoving curvature perturbation in an open Friedmann-Robertson-Walker universe, which is used to determine the temperature anisotropies of the cosmic microwave background radiation. For a fiducial Lagrangian, the anisotropies are well described by a power law in the cosine of the angular distance from the center of the collision signature. For a given form of the Lagrangian, the resulting observational predictions are inherently statistical due to stochastic elements of the bubble nucleation process. Further uncertainties arise due to our imperfect knowledge about inflationary and pre-recombination physics. We characterize observational predictions by computing the probability distributions over four phenomenological parameters which capture these intrinsic and model uncertainties. This represents the first fully-relativistic set of predictions from an ensemble of scalar field models giving rise to eternal inflation, yielding significant differences from previous non-relativistic approximations. Thus, our results provide a basis for a rigorous confrontation of these theories with cosmological data.

Throughout the Universe many powerful events are driven by strong gravitational effects that require general relativity to fully describe them. These include compact binary mergers, black hole accretion and stellar collapse, where velocities can approach the speed of light, and extreme gravitational fields --$\Phi_{\rm Newt}/c^2 \simeq 1$-- mediate the interactions. Many of these processes trigger emission across a broad range of the electromagnetic spectrum. Compact binaries further source strong gravitational wave emission that could directly be detected in the near future. This feat will open up a gravitational wave window into our Universe and revolutionize its understanding. Describing these phenomena requires general relativity, and --where dynamical effects strongly modify gravitational fields-- the full Einstein equations coupled to matter sources. Numerical relativity is a field within general relativity concerned with studying such scenarios that cannot be accurately modeled via perturbative or analytical calculations. In this review, we examine results obtained within this discipline, with a focus on its impact in astrophysics.

We incorporate realistic, tabulated equations of state into fully relativistic simulations of magnetized neutron stars along with a neutrino leakage scheme which accounts for cooling via neutrino emission. Both these improvements utilize open-source code (GR1D) and tables from http://stellarcollapse.org. Our implementation makes use of a novel method for the calculation of the optical depth which simplifies its use with distributed adaptive mesh refinement. We present various tests with and without magnetization and preliminary results both from single stars and from the merger of a binary system.

We present a framework to study generic neutron-star binaries in scalar-tensor theories of gravity. Our formalism achieves this goal by suitably interfacing a post-Newtonian orbital evolution (described by a set of ordinary differential equations) with a set of non-linear algebraic equations, which provide a description of the scalar charge of each binary's component along the evolution in terms of isolated-star data. We validate this semi-analytical procedure by comparing its results to those of fully general-relativistic simulations, and use it to investigate the behavior of binary systems in large portions of the parameter space of scalar-tensor theories. This allows us to shed further light on the phenomenon of "dynamical scalarization", which we uncovered in [Barausse, Palenzuela, Ponce and Lehner, Phys. Rev. D 87, 081506(R) (2013)] and which takes place in tight binaries, even for stars that have exactly zero scalar charge in isolation. We also employ our formalism to study representative binary systems, obtain their gravitational-wave signals and discuss the extent to which deviations from General Relativity can be detected. The insights gained by this framework allow us to additionally show that eccentric binaries can undergo scalarization/de-scalarization phenomena.

We show that rapidly-spinning black holes can display turbulent gravitational behavior which is mediated by a new type of parametric instability. This instability transfers energy from higher temporal and azimuthal spatial frequencies to lower frequencies--- a phenomenon reminiscent of the inverse energy cascade displayed by 2+1-dimensional turbulent fluids. Our finding reveals a path towards gravitational turbulence for perturbations of rapidly-spinning black holes, and provides the first evidence for gravitational turbulence in an asymptotically flat spacetime. Interestingly, this finding predicts observable gravitational wave signatures from such phenomena in black hole binaries with high spins and gives a gravitational description of turbulence relevant to the fluid-gravity duality.

We study the dynamics of a 2+1 dimensional relativistic viscous conformal fluid in Minkowski spacetime. Such fluid solutions arise as duals, under the "gravity/fluid correspondence", to 3+1 dimensional asymptotically anti-de Sitter (AAdS) black brane solutions to the Einstein equation. We examine stability properties of shear flows, which correspond to hydrodynamic quasinormal modes of the black brane. We find that, for sufficiently high Reynolds number, the solution undergoes an inverse turbulent cascade to long wavelength modes. We then map this fluid solution, via the gravity/fluid duality, into a bulk metric. This suggests a new and interesting feature of the behavior of perturbed AAdS black holes and black branes, which is not readily captured by a standard quasinormal mode analysis. Namely, for sufficiently large perturbed black objects (with long-lived quasinormal modes), nonlinear effects transfer energy from short to long wavelength modes via a turbulent cascade within the metric perturbation. As long wavelength modes have slower decay, this lengthens the overall lifetime of the perturbation. We also discuss various implications of this behavior, including expectations for higher dimensions, and the possibility of predicting turbulence in more general gravitational scenarios.

We expand on our study of the gravitational and electromagnetic emissions from the late stage of an inspiraling neutron star binary as presented in Ref. \cite{Palenzuela:2013hu}. Interactions between the stellar magnetospheres, driven by the extreme dynamics of the merger, can yield considerable outflows. We study the gravitational and electromagnetic waves produced during the inspiral and merger of a binary neutron star system using a full relativistic, resistive MHD evolution code. We show that the interaction between the stellar magnetospheres extracts kinetic energy from the system and powers radiative Poynting flux and heat dissipation. These features depend strongly on the configuration of the initial stellar magnetic moments. Our results indicate that this power can strongly outshine pulsars in binaries and have a distinctive angular and time-dependent pattern. Our discussion provides more detail than Ref. \cite{Palenzuela:2013hu}, showing clear evidence of the different effects taking place during the inspiral. Our simulations include a few milliseconds after the actual merger and study the dynamics of the magnetic fields during the formation of the hypermassive neutron star. We also briefly discuss the possibility of observing such emissions.

The late stage of an inspiraling neutron star binary gives rise to strong gravitational wave emission due to its highly dynamic, strong gravity. Moreover, interactions between the stellar magnetospheres can produce considerable electromagnetic radiation. We study this scenario using fully general relativistic, resistive magneto-hydrodynamics simulations. We show that these interactions extract kinetic energy from the system, dissipate heat, and power radiative Poynting flux, as well as develop current sheets. Our results indicate that this power can: (i) outshine pulsars in binaries, (ii) display a distinctive angular- and time-dependent pattern, and (iii) radiate within large opening angles. These properties suggest that some binary neutron star mergers are ideal candidates for multimessenger astronomy.

This work devises a formalism to obtain the equations of motion for a black hole-fluid configuration. Our approach is based on a Post-Newtonian expansion and adapted to scenarios where obtaining the relevant dynamics requires long time-scale evolutions. These systems are typically studied with Newtonian approaches, which have the advantage that larger time-steps can be employed than in full general-relativistic simulations, but have the downside that important physical effects are not accounted for. The formalism presented here provides a relatively straightforward way to incorporate those effects in existing implementations, up to 2.5PN order, with lower computational costs than fully relativistic simulations.

Interferometric detectors will very soon give us an unprecedented view of the gravitational-wave sky, and in particular of the explosive and transient Universe. Now is the time to challenge our theoretical understanding of short-duration gravitational-wave signatures from cataclysmic events, their connection to more traditional electromagnetic and particle astrophysics, and the data analysis techniques that will make the observations a reality. This paper summarizes the state of the art, future science opportunities, and current challenges in understanding gravitational-wave transients.

Scalar-tensor theories of gravity are natural phenomenological alternatives to General Relativity, where the gravitational interaction is mediated by a scalar degree of freedom, besides the usual tensor gravitons. In regions of the parameter space of these theories where constraints from both solar system experiments and binary-pulsar observations are satisfied, we show that binaries of neutron stars present marked differences from General Relativity in both the late-inspiral and merger phases. In particular, phenomena related to the spontaneous scalarization of isolated neutron stars take place in the late stages of the evolution of binary systems, with important effects in the ensuing dynamics. We comment on the relevance of our results for the upcoming Advanced LIGO/Virgo detectors.

We study the gravitational collapse of a magnetized neutron star using a novel numerical approach able to capture both the dynamics of the star and the behavior of the surrounding plasma. In this approach, a fully general relativistic magnetohydrodynamics implementation models the collapse of the star and provides appropriate boundary conditions to a force-free model which describes the stellar exterior. We validate this strategy by comparing with known results for the rotating monopole and aligned rotator solutions and then apply it to study both rotating and non-rotating stellar collapse scenarios, and contrast the behavior with what is obtained when employing the electrovacuum approximation outside the star. The non-rotating electrovacuum collapse is shown to agree qualitatively with a Newtonian model of the electromagnetic field outside a collapsing star. We illustrate and discuss a fundamental difference between the force-free and electrovacuum solutions, involving the appearance of large zones of electric-dominated field in the vacuum case. This provides a clear demonstration of how dissipative singularities appear generically in the non-linear time-evolution of force-free fluids. In both the rotating and non-rotating cases, our simulations indicate that the collapse induces a strong electromagnetic transient. In the case of sub-millisecond rotation, the magnetic field experiences strong winding and the transient carries much more energy. This result has important implications for models of gamma-ray bursts.

We perform fully non-linear numerical simulations of charged-black-hole collisions, described by the Einstein-Maxwell equations, and contrast the results against analytic expectations. We focus on head-on collisions of non-spinning black holes, starting from rest and with the same charge to mass ratio, Q/M. The addition of charge to black holes introduces a new interesting channel of radiation and dynamics, most of which seem to be captured by Newtonian dynamics and flat-space intuition. The waveforms can be qualitatively described in terms of three stages; (i) an infall phase prior to the formation of a common apparent horizon; (ii) a nonlinear merger phase which corresponds to a peak in gravitational and electromagnetic energy; (iii) the ringdown marked by an oscillatory pattern with exponentially decaying amplitude and characteristic frequencies that are in good agreement with perturbative predictions. We observe that the amount of gravitational-wave energy generated throughout the collision decreases by about three orders of magnitude as the charge-to-mass ratio Q/M is increased from 0 to 0.98. We interpret this decrease as a consequence of the smaller accelerations present for larger values of the charge. In contrast, the ratio of energy carried by electromagnetic to gravitational radiation increases, reaching about 22% for the maximum Q/M ratio explored, which is in good agreement with analytic predictions.

Cosmic bubble collisions provide an important possible observational window on the dynamics of eternal inflation. In eternal inflation, our observable universe is contained in one of many bubbles formed from an inflating metastable vacuum. The collision between bubbles can leave a detectable imprint on the cosmic microwave background radiation. Although phenomenological models of the observational signature have been proposed, to make the theory fully predictive one must determine the bubble collision spacetime, and thus the cosmological observables, from a scalar field theory giving rise to eternal inflation. Because of the intrinsically non-linear nature of the bubbles and their collision, this requires a numerical treatment incorporating General Relativity. In this paper, we present results from numerical simulations of bubble collisions in full General Relativity. These simulations allow us to accurately determine the outcome of bubble collisions, and examine their effect on the cosmology inside a bubble universe. We confirm the validity of a number of approximations used in previous analytic work, and identify qualitatively new features of bubble collision spacetimes. Both vacuum bubbles and bubbles containing a realistic inflationary cosmology are studied. We identify the constraints on the scalar field potential that must be satisfied in order to obtain collisions that are consistent with our observed cosmology, yet leave detectable signatures.

Physics in curved spacetime describes a multitude of phenomena, ranging from astrophysics to high energy physics. The last few years have witnessed further progress on several fronts, including the accurate numerical evolution of the gravitational field equations, which now allows highly nonlinear phenomena to be tamed. Numerical relativity simulations, originally developed to understand strong field astrophysical processes, could prove extremely useful to understand high-energy physics processes like trans-Planckian scattering and gauge-gravity dualities. We present a concise and comprehensive overview of the state-of-the-art and important open problems in the field(s), along with guidelines for the next years. This writeup is a summary of the "NR/HEP Workshop" held in Madeira, Portugal from August 31st to September 3rd 2011.

The Blandford-Znajek mechanism has long been regarded as a key ingredient in models attempting to explain powerful jets in AGNs, quasars, blazzars etc. In such mechanism, energy is extracted from a rotating black hole and dissipated at a load at far distances. In the current work we examine the behaviour of the BZ mechanism with respect to different boundary conditions, revealing the mechanism robustness upon variation of these conditions. Consequently, this work closes a gap in our understanding of this important scenario.

The extraction of rotational energy from a spinning black hole via the Blandford-Znajek mechanism has long been understood as an important component in models to explain energetic jets from compact astrophysical sources. Here we show more generally that the kinetic energy of the black hole, both rotational and translational, can be tapped, thereby producing even more luminous jets powered by the interaction of the black hole with its surrounding plasma. We study the resulting Poynting jet that arises from single boosted black holes and binary black hole systems. In the latter case, we find that increasing the orbital angular momenta of the system and/or the spins of the individual black holes results in an enhanced Poynting flux.

We investigate the merger of a neutron star (of compaction ratio $0.1$) in orbit about a spinning black hole in full general relativity with a mass ratio of $5:1$, allowing for the star to have an initial magnetization of $10^{12} {\rm Gauss}$. We present the resulting gravitational waveform and analyze the fallback accretion as the star is disrupted. The evolutions suggest no significant effects from the initial magnetization. We find that only a negligible amount of matter becomes unbound; $99%$ of the neutron star material has a fallback time of 10 seconds or shorter to reach the region of the central engine and that $99.99%$ of the star will interact with the central disk and black hole within 3 hours.

The interaction of black holes with ambient magnetic fields is important for a variety of highly energetic astrophysical phenomena. We study this interaction within the force-free approximation in which a tenuous plasma is assumed to have zero inertia. Blandford and Znajek (BZ) used this approach to demonstrate the conversion of some of the black hole's energy into electromagnetic Poynting flux in stationary and axisymmetric single black hole systems. We adopt this approach and extend it to examine asymmetric and, most importantly, dynamical systems by implementing the fully nonlinear field equations of general relativity coupled to Maxwell's equations. For single black holes, we study in particular the dependence of the Poynting flux and show that, even for misalignments between the black hole spin and the direction of the asymptotic magnetic field, a Poynting flux is generated with a luminosity dependent on such misalignment. For binary black hole systems, we show both in the head-on and orbiting cases that the moving black holes generate a Poynting flux.

We study the evolution of magnetized and rigidly rotating neutron stars within a fully general relativistic implementation of ideal magnetohydrodynamics with no assumed symmetries in three spatial dimensions. The stars are modeled as rotating, magnetized polytropic stars and we examine diverse scenarios to study their dynamics and stability properties. In particular we concentrate on the stability of the stars and possible critical behavior. In addition to their intrinsic physical significance, we use these evolutions as further tests of our implementation which incorporates new developments to handle magnetized systems.

Supermassive black holes are found at the centers of most galaxies and their inspiral is a natural outcome when galaxies merge. The inspiral of these systems is of utmost astrophysical importance as prodigious producers of gravitational waves and in their possible role in energetic electromagnetic events. We study such binary black hole coalescence under the influence of an external magnetic field produced by the expected circumbinary disk surrounding them. Solving the Einstein equations to describe the spacetime and using the force-free approach for the electromagnetic fields and the tenuous plasma, we present numerical evidence for possible jets driven by these systems. Extending the process described by Blandford and Znajek for a single spinning black hole, the picture that emerges suggests the electromagnetic field extracts energy from the orbiting black holes, which ultimately merge and settle into the standard Blandford-Znajek scenario. Emissions along dual and single jets would be expected that could be observable to large distances.

As one step towards a systematic modeling of the electromagnetic (EM) emission from an inspiralling black hole binary we consider a simple scenario in which the binary moves in a uniform magnetic field anchored to a distant circumbinary disc. We study this system by solving the Einstein-Maxwell equations in which the EM fields are chosen with astrophysically consistent strengths. We consider binaries with spins aligned or anti-aligned with the orbital angular momentum and study the dependence of gravitational and EM signals with these spin configurations. Overall we find that the EM radiation in the lowest l=2, m=2 multipole accurately reflects the gravitational one, with identical phase evolutions and amplitudes that differ only by a scaling factor. We also compute the efficiency of the energy emission in EM waves and find that it is given by E^rad_EM/M ~ 10^-15 (M/10^8 M_Sun)^2 (B/10^4 G)^2, hence 13 orders of magnitude smaller than the gravitational energy for realistic magnetic fields. The corresponding luminosity is much smaller than the accretion luminosity if the system is accreting at near the Eddington rate. Most importantly, this EM emission is at frequencies of 10^-4 (10^8 M_Sun/M) Hz, well outside those accessible to astronomical radio observations. As a result, it is unlikely that the EM emission discussed here can be detected directly and simultaneously with the gravitational-wave one. However, indirect processes, driven by changes in the EM fields behavior could yield observable events. In particular if the accretion rate of the circumbinary disc is small and sufficiently stable over the timescale of the final inspiral, then the EM emission may be observable indirectly as it will alter the accretion rate through the magnetic torques exerted by the distorted magnetic field lines.

In addition to producing loud gravitational waves (GW), the dynamics of a binary black hole system could induce emission of electromagnetic (EM) radiation by affecting the behavior of plasmas and electromagnetic fields in their vicinity. We here study how the electromagnetic fields are affected by a pair of orbiting black holes through the merger. In particular, we show how the binary's dynamics induce a variability in possible electromagnetically induced emissions as well as an enhancement of electromagnetic fields during the late-merge and merger epochs. These time dependent features will likely leave their imprint in processes generating detectable emissions and can be exploited in the detection of electromagnetic counterparts of gravitational waves.

We simulate the possible emission from a disk perturbed by a recoiling super-massive black hole. To this end, we study radiation transfer from the system incorporating bremsstrahlung emission from a Maxwellian plasma and absorption given by Kramer's opacity law modified to incorporate blackbody effects. We employ this model in the radiation transfer integration to compute the luminosity at several frequencies, and compare with previous bremsstrahlung luminosity estimations from a transparent limit (in which the emissivity is integrated over the computational domain and over all frequencies) and with a simple thermal emission model. We find close agreement between the radiation transfer results and the estimated bremsstrahlung luminosity from previous work for electromagnetic signals above $10^{14}$ Hz. For lower frequencies, we find a self-eclipsing behavior in the disk, resulting in a strong intensity variability connected to the orbital period of the disk.

In addition to producing gravitational waves (GW), the dynamics of a binary black hole system could induce emission of electromagnetic (EM) radiation by affecting the behavior of plasmas and electromagnetic fields in their vicinity. We here study how the electromagnetic fields are affected by a pair of orbiting black holes through the merger. In particular, we show how the binary's dynamics induce a variability in possible electromagnetically induced emissions as well as a possible enhancement of electromagnetic fields during the late-merge and merger epochs. These time dependent features will likely leave their imprint in processes generating detectable emissions and can be exploited in the detection of electromagnetic counterparts of gravitational waves.

The merger process of a binary black hole system can have a strong impact on a circumbinary disk. In the present work we study the effect of both central mass reduction (due to the energy loss through gravitational waves) and a possible black hole recoil (due to asymmetric emission of gravitational radiation). For the mass reduction case and recoil directed along the disk's angular momentum, oscillations are induced in the disk which then modulate the internal energy and bremsstrahlung luminosities. On the other hand, when the recoil direction has a component orthogonal to the disk's angular momentum, the disk's dynamics are strongly impacted, giving rise to relativistic shocks. The shock heating leaves its signature in our proxies for radiation, the total internal energy and bremsstrahlung luminosity. Interestingly, for cases where the kick velocity is below the smallest orbital velocity in the disk (a likely scenario in real AGN), we observe a common, characteristic pattern in the internal energy of the disk. Variations in kick velocity simply provide a phase offset in the characteristic pattern implying that observations of such a signature could yield a measure of the kick velocity through electromagnetic signals alone.

Many astrophysical processes involving magnetic fields and quasi-stationary processes are well described when assuming the fluid as a perfect conductor. For these systems, the ideal-magnetohydrodynamics (MHD) description captures the dynamics effectively and a number of well-tested techniques exist for its numerical solution. Yet, there are several astrophysical processes involving magnetic fields which are highly dynamical and for which resistive effects can play an important role. The numerical modeling of such non-ideal MHD flows is significantly more challenging as the resistivity is expected to change of several orders of magnitude across the flow and the equations are then either of hyperbolic-parabolic nature or hyperbolic with stiff terms. We here present a novel approach for the solution of these relativistic resistive MHD equations exploiting the properties of implicit-explicit (IMEX) Runge Kutta methods. By examining a number of tests we illustrate the accuracy of our approach under a variety of conditions and highlight its robustness when compared with alternative methods, such as the Strang-splitting. Most importantly, we show that our approach allows one to treat, within a unified framework, both those regions of the flow which are fluid-pressure dominated (such as in the interior of compact objects) and those which are instead magnetic-pressure dominated (such as in their magnetospheres)

We investigate the influence of magnetic fields upon the dynamics of and resulting gravitational waves from a binary neutron star merger in full general relativity coupled to ideal magnetohydrodynamics (MHD). We consider two merger scenarios, one where the stars begin with initially aligned poloidal magnetic fields and one with no magnetic field. Both mergers result in a strongly differentially rotating object. In comparison to the non-magnetized scenario, the aligned magnetic fields delay the final merger of the two stars. During and after merger we observe phenomena driven by the magnetic field, including Kelvin-Helmholtz instabilities in shear layers, winding of the field lines, and transition from poloidal to toroidal fields. These effects not only produce electromagnetic radiation, but also can have a strong influence on the gravitational waves. Thus, there are promising prospects for studying such systems with both types of waves.

We study and develop constraint preserving boundary conditions for the Newtonian magnetohydrodynamic equations and analyze the behavior of the numerical solution upon considering different possible options.

We model two mergers of orbiting binary neutron stars, the first forming a black hole and the second a differentially rotating neutron star. We extract gravitational waveforms in the wave zone. Comparisons to a post-Newtonian analysis allow us to compute the orbital kinematics, including trajectories and orbital eccentricities. We verify our code by evolving single stars and extracting radial perturbative modes, which compare very well to results from perturbation theory. The Einstein equations are solved in a first order reduction of the generalized harmonic formulation, and the fluid equations are solved using a modified convex essentially non-oscillatory method. All calculations are done in three spatial dimensions without symmetry assumptions. We use the \had computational infrastructure for distributed adaptive mesh refinement.

We present a straightforward approach for estimating the final black hole spin of a binary black hole coalescence with arbitrary initial masses and spins. Making some simple assumptions, we estimate the final angular momentum to be the sum of the individual spins plus the orbital angular momentum of a test particle orbiting at the last stable orbit around a Kerr black hole with a spin parameter of the final black hole. The formula we obtain is able to reproduce with reasonable accuracy the results from available numerical simulations, but, more importantly, it can be used to investigate what configurations might give rise to interesting dynamics. In particular, we discuss scenarios which might give rise to a ``flip'' in the direction of the total angular momentum of the system. By studying the dependence of the final spin upon the mass ratio and initial spins we find that our simple approach suggests that it is not possible to spin-up a black hole to extremal values through merger scenarios irrespective of the mass ratio of the objects involved.

We extend our previous studies of head-on collisions of boson stars by considering orbiting binary boson stars. We concentrate on equal mass binaries and study the dynamical behavior of boson/boson and boson/antiboson pairs. We examine the gravitational wave output of these binaries and compare with other compact binaries. Such a comparison lets us probe the apparent simplicity observed in gravitational waves produced by black hole binary systems. In our system of interest however, there is an additional internal freedom which plays a significant role in the system's dynamics, namely the phase of each star. Our evolutions show rather simple behavior at early times, but large differences occur at late times for the various initial configurations.

We study head-on collisions of boson stars in three dimensions. We consider evolutions of two boson stars which may differ in their phase or have opposite frequencies but are otherwise identical. Our studies show that these phase differences result in different late time behavior and gravitational wave output.

Computer simulations are enabling researchers to investigate systems which are extremely difficult to handle analytically. In the particular case of General Relativity, numerical models have proved extremely valuable for investigations of strong field scenarios and been crucial to reveal unexpected phenomena. Considerable efforts are being spent to simulate astrophysically relevant simulations, understand different aspects of the theory and even provide insights in the search for a quantum theory of gravity. In the present article I review the present status of the field of Numerical Relativity, describe the techniques most commonly used and discuss open problems and (some) future prospects.

We present the first simulations of non-headon (grazing) collisions of binary black holes in which the black hole singularities have been excised from the computational domain. Initially two equal mass black holes $m$ are separated a distance $\approx10m$ and with impact parameter $\approx2m$. Initial data are based on superposed, boosted (velocity $\approx0.5c$) solutions of single black holes in Kerr-Schild coordinates. Both rotating and non-rotating black holes are considered. The excised regions containing the singularities are specified by following the dynamics of apparent horizons. Evolutions of up to $t \approx 35m$ are obtained in which two initially separate apparent horizons are present for $t\approx3.8m$. At that time a single enveloping apparent horizon forms, indicating that the holes have merged. Apparent horizon area estimates suggest gravitational radiation of about 2.6% of the total mass. The evolutions end after a moderate amount of time because of instabilities.