Search

3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures and understand the physical conditions in the compact region of 3C84. We used EHT 228GHz observations and, given the limited (u,v)-coverage, applied geometric model fitting to the data. We also employed quasi-simultaneously observed, multi-frequency VLBI data for the source in order to carry out a comprehensive analysis of the core structure. We report the detection of a highly ordered, strong magnetic field around the central, SMBH of 3C84. The brightness temperature analysis suggests that the system is in equipartition. We determined a turnover frequency of $\nu_m=(113\pm4)$GHz, a corresponding synchrotron self-absorbed magnetic field of $B_{SSA}=(2.9\pm1.6)$G, and an equipartition magnetic field of $B_{eq}=(5.2\pm0.6)$G. Three components are resolved with the highest fractional polarisation detected for this object ($m_\textrm{net}=(17.0\pm3.9)$%). The positions of the components are compatible with those seen in low-frequency VLBI observations since 2017-2018. We report a steeply negative slope of the spectrum at 228GHz. We used these findings to test models of jet formation, propagation, and Faraday rotation in 3C84. The findings of our investigation into different flow geometries and black hole spins support an advection-dominated accretion flow in a magnetically arrested state around a rapidly rotating supermassive black hole as a model of the jet-launching system in the core of 3C84. However, systematic uncertainties due to the limited (u,v)-coverage, however, cannot be ignored., Comment: 15 pages, 6 figures, published in A&A

The Next Generation Very Large Array (ngVLA) is a planned radio interferometer providing unprecedented sensitivity at wavelengths between 21 cm and 3 mm. Its 263 antenna element array will be spatially distributed across North America to enable both superb low surface brightness recovery and sub-milliarcsecond angular resolution imaging. The project was developed by the international astronomy community under the lead of the National Radio Astronomy Observatory (NRAO), and is anticipated to be built between 2027 and 2037. Two workshops have been held in 2022 and 2023 with the goal to discuss and consolidate the scientific interests in the ngVLA within the German astronomical community. This community paper constitutes a collection of 48 science ideas which the German community aims to pursue with the ngVLA in the 2030s. This is not a complete list and the ideas are not developed at the level of a "Science Book", such that the present document is mainly meant provide a basis for further discussion within the community. As such, additional contributions are welcome, and will be considered for inclusion in future revisions., Comment: Version 2.0 (status June 18, 2024): 169 pages, comments and future contributions welcome [v2.0: 7 new science cases added, some minor revisions to other chapters]

The high stellar density in the central parsecs around the Galactic Centre makes it a seemingly favourable environment for finding relativistic binary pulsars. These include pulsars orbiting other neutron stars, stellar-mass black holes or the central supermassive black hole, Sagittarius A*. Here we present multi-epoch pulsar searches of the Galactic Centre at four observing frequencies, (4.85, 8.35, 14.6 18.95) GHz, using the Effelsberg 100-m radio telescope. Observations were conducted one year prior to the discovery of, and during monitoring observations of, the Galactic Centre magnetar PSR J1745-2900. Our data analysis features acceleration searches on progressively shorter time series to maintain sensitivity to relativistic binary pulsars. The multi-epoch observations increase the likelihood of discovering transient or nulling pulsars, or ensure orbital phases are observed at which acceleration search methods work optimally. In ~147 h of separate observations, no previously undiscovered pulsars have been detected. Through calibration observations, we conclude this might be due to insufficient instantaneous sensitivity; caused by the intense continuum emission from the Galactic Centre, its large distance and, at higher frequencies, the aggregate effect of steep pulsar spectral indices and atmospheric contributions to the system temperature. Additionally we find that for millisecond pulsars in wide circular orbits ~<800 d around Sagittarius A*, linear acceleration effects cannot be corrected in deep observations (9 h) with existing software tools. Pulsar searches of the Galactic Centre with the next generation of radio telescopes - such as MeerKat, ngVLA and SKA1-mid - will have improved chances of uncovering this elusive population., Comment: 17 pages, 12 figures, Accepted for publication in Monthly Notices of the Royal Astronomical Society

THESEUS, one of the two space mission concepts being studied by ESA as candidates for next M5 mission within its Comsic Vision programme, aims at fully exploiting Gamma-Ray Bursts (GRB) to solve key questions about the early Universe, as well as becoming a cornerstone of multi-messenger and time-domain astrophysics. By investigating the first billion years of the Universe through high-redshift GRBs, THESEUS will shed light on the main open issues in modern cosmology, such as the population of primordial low mass and luminosity galaxies, sources and evolution of cosmic re-ionization, SFR and metallicity evolution up to the "cosmic dawn" and across Pop-III stars. At the same time, the mission will provide a substantial advancement of multi-messenger and time-domain astrophysics by enabling the identification, accurate localisation and study of electromagnetic counterparts to sources of gravitational waves and neutrinos, which will be routinely detected in the late '20s and early '30s by the second and third generation Gravitational Wave (GW) interferometers and future neutrino detectors, as well as of all kinds of GRBs and most classes of other X/gamma-ray transient sources. In all these cases, THESEUS will provide great synergies with future large observing facilities in the multi-messenger domain. A Guest Observer programme, comprising Target of Opportunity (ToO) observations, will expand the science return of the mission, to include, e.g., solar system minor bodies, exoplanets, and AGN., Comment: Submitted to Experimental Astronomy. arXiv admin note: text overlap with arXiv:2102.08702

The Event Horizon Telescope (EHT) will soon provide the first high-resolution images of the Galactic Centre supermassive black hole (SMBH) candidate Sagittarius A* (Sgr A*), enabling us to probe gravity in the strong-field regime. Besides studying the accretion process in extreme environments, the obtained data and reconstructed images could be used to investigate the underlying spacetime structure. In its current configuration, the EHT is able to distinguish between a rotating Kerr black hole and a horizon-less object like a boson star. Future developments can increase the ability of the EHT to tell different spacetimes apart. We investigate the capability of an advanced EHT concept, including an orbiting space antenna, to image and distinguish different spacetimes around Sgr A*. We use GRMHD simulations of accreting compact objects (Kerr and dilaton black holes, as well as boson stars) and compute their radiative signatures via general relativistic radiative transfer calculations. To facilitate comparison with upcoming and future EHT observations we produce realistic synthetic data including the source variability, diffractive and refractive scattering while incorporating the observing array, including a space antenna. From the generated synthetic observations we dynamically reconstructed black hole shadow images using regularised Maximum Entropy methods. We employ a genetic algorithm to optimise the orbit of the space antenna with respect to improved imaging capabilities and u-v-plane coverage of the combined array (ground array and space antenna and developed a new method to probe the source variability in Fourier space. The inclusion of an orbiting space antenna improves the capability of the EHT to distinguish the spin of Kerr black holes and dilaton black holes based on reconstructed radio images and complex visibilities., Comment: 15 pages, 14 figures, A&A accepted

We present a new method for the numerical solution of the radiative-transfer equation (RTE) in multidimensional scenarios commonly encountered in computational astrophysics. The method is based on the direct solution of the Boltzmann equation via an extension of the Lattice Boltzmann (LB) equation and allows to model the evolution of the radiation field as it interacts with a background fluid, via absorption, emission, and scattering. As a first application of this method, we restrict our attention to a frequency independent ("grey") formulation within a special-relativistic framework, which can be employed also for classical computational astrophysics. For a number of standard tests that consider the performance of the method in optically thin, optically thick and intermediate regimes with a static fluid, we show the ability of the LB method to produce accurate and convergent results matching the analytic solutions. We also contrast the LB method with commonly employed moment-based schemes for the solution of the RTE, such as the M1 scheme. In this way, we are able to highlight that the LB method provides the correct solution for both non-trivial free-streaming scenarios and the intermediate optical-depth regime, for which the M1 method either fails or provides inaccurate solutions. When coupling to a dynamical fluid, on the other hand, we present the first self-consistent solution of the RTE with LB methods within a relativistic-hydrodynamic scenario. Finally, we show that besides providing more accurate results in all regimes, the LB method features smaller or comparable computational costs compared to the M1 scheme., Comment: 22 pages, 16 figures, matches version accepted in MNRAS

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images., Comment: 20 pages, 15 figures, accepted for publication in A&A

The galaxy M 87 harbors a kiloparsec-scale relativistic jet, whose origin coincides with a supermassive black hole. Observational mm-VLBI campaigns are capable of resolving the jet-launching region at the scale of the event horizon. In order to provide a context for interpreting these observations, realistic general-relativistic magnetohydrodynamical (GRMHD) models of the accretion flow are constructed. The characteristics of the observed spectral-energy distribution (SED) depend on the shape of the electrons' energy-distribution function (eDF). The dependency on the eDF is omitted in the modeling of the first Event Horizon Telescope results. In this work, we aim to model the M 87 SED from radio up to NIR/optical frequencies using a thermal-relativistic Maxwell- J\"uttner distribution, as well as a relativistic $\kappa$-distribution function. The electrons are injected based on sub-grid, particle-in-cell parametrizations for sub-relativistic reconnection. A GRMHD simulation in Cartesian-Kerr-Schild coordinates, using eight levels of adaptive mesh refinement (AMR), forms the basis of our model. To obtain spectra and images, the GRMHD data is post-processed with the ray-tracing code RAPTOR, which is capable of ray tracing through AMR GRMHD simulation data. We obtain radio spectra in both the thermal-jet and $\kappa$-jet models consistent with radio observations. Additionally, the $\kappa$-jet models also recover the NIR/optical emission. The models recover the observed source sizes and core shifts and obtain a jet power of $\approx 10^{43}$ ergs/s. In the $\kappa$-jet models, both the accretion rates and jet powers are approximately two times lower than the thermal-jet model. The frequency cut-off observed at $\nu \approx 10^{15}$ Hz is recovered when the accelerator size is $10^6$ - $10^8$ cm, this could potentially point to an upper limit for plasmoid sizes in the jet of M 87., Comment: 15 pages, 10 figures

High-resolution Very-Long-Baseline Interferometry observations of NGC 1052 show a two sided jet with several regions of enhanced emission and a clear emission gap between the two jets.This gap shrinks with increasing frequency and vanishes around $\nu\sim43$ GHz. The observed structures are due to both the macroscopic fluid dynamics interacting with the surrounding ambient medium including an obscuring torus and the radiation microphysics. In this paper we investigate the possible physical conditions in relativistic jets of NGC 1052 by directly modelling the observed emission and spectra via state-of-the-art special-relativistic hydrodynamic (SRHD) simulations and radiative transfer calculations. To investigate the physical conditions in the relativistic jet we coupled our radiative transfer code to evolutionary algorithms and performed simultaneous modelling of the observed jet structure and the broadband radio spectrum. During the calculation of the radiation we consider both thermal and non-thermal emission. In order to compare our model to VLBI observations we take into account the sparse sampling of the u-v plane, the array properties and the imaging algorithm. We present for the first time an end-to-end pipeline for fitting numerical simulations to VLBI observations of relativistic jets taking into account the macrophysics including fluid dynamics and ambient medium configurations together with thermal/non-thermal emission and the properties of the observing array. The detailed analysis of our simulations shows that the structure and properties of the observed relativistic jets in NGC 1052 can be reconstructed by a slightly over-pressured jet ($d_k\sim1.5$) embedded in a decreasing pressure ambient medium, Comment: submitted to A&A, 24 pages, 25 figures

The recent discovery of the electromagnetic counterpart of the gravitational wave source GW170817, has demonstrated the huge informative power of multi-messenger observations. During the next decade the nascent field of multi-messenger astronomy will mature significantly. Around 2030, third generation gravitational wave detectors will be roughly ten times more sensitive than the current ones. At the same time, neutrino detectors currently upgrading to multi km^3 telescopes, will include a 10 km^3 facility in the Southern hemisphere that is expected to be operational around 2030. In this review, we describe the most promising high frequency gravitational wave and neutrino sources that will be detected in the next two decades. In this context, we show the important role of the Transient High Energy Sky and Early Universe Surveyor (THESEUS), a mission concept proposed to ESA by a large international collaboration in response to the call for the Cosmic Vision Programme M5 missions. THESEUS aims at providing a substantial advancement in early Universe science as well as playing a fundamental role in multi-messenger and time-domain astrophysics, operating in strong synergy with future gravitational wave and neutrino detectors as well as major ground- and space-based telescopes. This review is an extension of the THESEUS white paper (Amati et al. 2017), also in light of the discovery of GW170817/GRB170817A that was announced on October 16th, 2017., Comment: 25 pages, 13 figures, accepted to Advances in Space Research with minor revisions. Details on the THESEUS instrumentation, science case and expected performances can also be found in Amati et al. 2017 (arXiv:1710.04638) and in the presentations of the THESEUS Workshop 2017 (http://www.isdc.unige.ch/theseus/workshop2017-programme.html). v2 few typos corrected

THESEUS is a space mission concept aimed at exploiting Gamma-Ray Bursts for investigating the early Universe and at providing a substantial advancement of multi-messenger and time-domain astrophysics. These goals will be achieved through a unique combination of instruments allowing GRB and X-ray transient detection over a broad field of view (more than 1sr) with 0.5-1 arcmin localization, an energy band extending from several MeV down to 0.3 keV and high sensitivity to transient sources in the soft X-ray domain, as well as on-board prompt (few minutes) follow-up with a 0.7 m class IR telescope with both imaging and spectroscopic capabilities. THESEUS will be perfectly suited for addressing the main open issues in cosmology such as, e.g., star formation rate and metallicity evolution of the inter-stellar and intra-galactic medium up to redshift $\sim$10, signatures of Pop III stars, sources and physics of re-ionization, and the faint end of the galaxy luminosity function. In addition, it will provide unprecedented capability to monitor the X-ray variable sky, thus detecting, localizing, and identifying the electromagnetic counterparts to sources of gravitational radiation, which may be routinely detected in the late '20s / early '30s by next generation facilities like aLIGO/ aVirgo, eLISA, KAGRA, and Einstein Telescope. THESEUS will also provide powerful synergies with the next generation of multi-wavelength observatories (e.g., LSST, ELT, SKA, CTA, ATHENA)., Comment: Accepted for publication in Advances in Space Research. Partly based on the proposal submitted on October 2016 in response to the ESA Call for next M5 mission, with expanded and updated science sections

High-resolution Very-Long-Baseline Interferometry observations of active galactic nuclei have revealed asymmetric structures in the jets of radio galaxies. These asymmetric structures may be due to internal asymmetries in the jet, could be induced by the different conditions in the surrounding ambient medium including the obscuring torus, or a combination of the two. In this paper we investigate the influence of the ambient medium (including the obscuring torus) on the observed properties of jets from radio galaxies. We performed special-relativistic hydrodynamic (RHD) simulations of over-pressured and pressure-matched jets using the special-relativistic hydrodynamics code \texttt{Ratpenat}, which is based on a second-order accurate finite-volume method and an approximate Riemann solver. Using a newly developed emission code to compute the electromagnetic emission, we have investigated the influence of different ambient medium and torus configurations on the jet structure and subsequently computed the non-thermal emission produced by the jet and the thermal absorption due to the torus. To better compare the emission simulations with observations we produced synthetic radio maps, taking into account the properties of the observatory. The detailed analysis of our simulations shows that the observed asymmetries can be produced by the interaction of the jet with the ambient medium and by the absorption properties of the obscuring torus., Comment: 14 pages, 17 figures, submitted to A&A

Einstein's General Theory of Relativity (GR) successfully describes gravity. The most fundamental predictions of GR are black holes (BHs), but in spite of many convincing BH candidates in the Universe, there is no conclusive experimental proof of their existence using astronomical observations in the electromagnetic spectrum. Are BHs real astrophysical objects? Does GR hold in its most extreme limit or are alternatives needed? The prime target to address these fundamental questions is in the center of our own Galaxy, which hosts the closest and best-constrained supermassive BH candidate in the Universe, Sagittarius A* (Sgr A*). Three different types of experiments hold the promise to test GR in a strong-field regime using observations of Sgr A* with new-generation instruments. The first experiment aims to image the relativistic plasma emission which surrounds the event horizon and forms a "shadow" cast against the background, whose predicted size (~50 microarcseconds) can now be resolved by upcoming VLBI experiments at mm-waves such as the Event Horizon Telescope (EHT). The second experiment aims to monitor stars orbiting Sgr A* with the upcoming near-infrared interferometer GRAVITY at the Very Large Telescope (VLT). The third experiment aims to time a radio pulsar in tight orbit about Sgr A* using radio telescopes (including the Atacama Large Millimeter Array or ALMA). The BlackHoleCam project exploits the synergy between these three different techniques and aims to measure the main BH parameters with sufficient precision to provide fundamental tests of GR and probe the spacetime around a BH in any metric theory of gravity. Here, we review our current knowledge of the physical properties of Sgr A* as well as the current status of such experimental efforts towards imaging the event horizon, measuring stellar orbits, and timing pulsars around Sgr A*., Comment: review paper, 36 pages, 12 figures (v2 is the published version)

A large international effort is under way to assess the presence of a shadow in the radio emission from the compact source at the centre of our Galaxy, Sagittarius A$^*$ (Sgr A$^*$). If detected, this shadow would provide the first direct evidence of the existence of black holes and that Sgr A$^*$ is a supermassive black hole. In addition, the shape of the shadow could be used to learn about extreme gravity near the event horizon and to determine which theory of gravity better describes the observations. The mathematical description of the shadow has so far used a number of simplifying assumptions that are unlikely to be met by the real observational data. We here provide a general formalism to describe the shadow as an arbitrary polar curve expressed in terms of a Legendre expansion. Our formalism does not presume any knowledge of the properties of the shadow, e.g. the location of its centre, and offers a number of routes to characterize the distortions of the curve with respect to reference circles. These distortions can be implemented in a coordinate independent manner by different teams analysing the same data. We show that the new formalism provides an accurate and robust description of noisy observational data, with smaller error variances when compared to previous approaches for the measurement of the distortion., Comment: 12 pages, 8 figures, accepted for publication on MNRAS

We have studied the dynamics of an equal-mass magnetized neutron-star binary within a resistive magnetohydrodynamic (RMHD) approach in which the highly conducting stellar interior is matched to an electrovacuum exterior. Because our analysis is aimed at assessing the modifications introduced by resistive effects on the dynamics of the binary after the merger and through to collapse, we have carried out a close comparison with an equivalent simulation performed within the traditional ideal magnetohydrodynamic approximation. We have found that there are many similarities between the two evolutions but also one important difference: the survival time of the hyper massive neutron star increases in a RMHD simulation. This difference is due to a less efficient magnetic-braking mechanism in the resistive regime, in which matter can move across magnetic-field lines, thus reducing the outward transport of angular momentum. Both the RMHD and the ideal magnetohydrodynamic simulations carried here have been performed at higher resolutions and with a different grid structure than those in previous work of ours [L. Rezzolla, B. Giacomazzo, L. Baiotti, J. Granot, C. Kouveliotou, and M. A. Aloy, Astrophys. J. Letters 732, L6 (2011)], but confirm the formation of a low-density funnel with an ordered magnetic field produced by the black hole--torus system. In both regimes the magnetic field is predominantly toroidal in the highly conducting torus and predominantly poloidal in the nearly evacuated funnel. Reconnection processes or neutrino annihilation occurring in the funnel, none of which we model, could potentially increase the internal energy in the funnel and launch a relativistic outflow, which, however, is not produced in these simulations., Comment: 26 pages, 17 figures; animations available at http://www.southampton.ac.uk/~kd10g13/movies/index.shtml

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study., Comment: Proc. SPIE 9144, Space Telescopes and Instrumentation 2014: Ultraviolet to Gamma Ray, 91442T

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants, Comment: Replaced figures 2 and 3: corrected position SRT. Corrected minor typo in 5.1

The qualitative and quantitative understanding of near-horizon gravitational dynamics in the strong-field regime represents a challenge both at a fundamental level and in astrophysical applications. Recent advances in numerical relativity and in the geometric characterization of black hole horizons open new conceptual and technical avenues into the problem. We discuss here a research methodology in which spacetime dynamics is probed through the cross-correlation of geometric quantities constructed on the black hole horizon and on null infinity. These two hypersurfaces respond to evolving gravitational fields in the bulk, providing canonical "test screens" in a "scattering"-like perspective onto spacetime dynamics. More specifically, we adopt a 3+1 Initial Value Problem approach to the construction of generic spacetimes and discuss the role and properties of dynamical trapping horizons as canonical inner "screens" in this context. We apply these ideas and techniques to the study of the recoil dynamics in post-merger binary black holes, an important issue in supermassive galactic black hole mergers., Comment: 16 pages, 5 figures, contribution to the proceedings volume of the Spanish Relativity Meeting ERE2011: "Towards new paradigms", Madrid, Spain, 29 Aug-2 Sep 2011

Einstein gravitational-wave Telescope (ET) is a design study funded by the European Commission to explore the technological challenges of and scientific benefits from building a third generation gravitational wave detector. The three-year study, which concluded earlier this year, has formulated the conceptual design of an observatory that can support the implementation of new technology for the next two to three decades. The goal of this talk is to introduce the audience to the overall aims and objectives of the project and to enumerate ET's potential to influence our understanding of fundamental physics, astrophysics and cosmology., Comment: Conforms to conference proceedings, several author names corrected

We discuss different ways that neutron stars can generate gravitational waves, describe recent improvements in modelling the relevant scenarios in the context of improving detector sensitivity, and show how observations are beginning to test our understanding of fundamental physics. The main purpose of the discussion is to establish promising science goals for third-generation ground-based detectors, like the Einstein Telescope, and identify the various challenges that need to be met if we want to use gravitational-wave data to probe neutron star physics., Comment: review article for Einstein Telescope special issue in Gen Rel Grav, 39 pages, 7 figures

(Abridged.) The accretion-induced collapse (AIC) of a white dwarf (WD) may lead to the formation of a protoneutron star and a collapse-driven supernova explosion. This process represents a path alternative to thermonuclear disruption of accreting white dwarfs in Type Ia supernovae. Neutrino and gravitational-wave (GW) observations may provide crucial information necessary to reveal a potential AIC. Motivated by the need for systematic predictions of the GW signature of AIC, we present results from an extensive set of general-relativistic AIC simulations using a microphysical finite-temperature equation of state and an approximate treatment of deleptonization during collapse. Investigating a set of 114 progenitor models in rotational equilibrium, with a wide range of rotational configurations, temperatures and central densities, we extend previous Newtonian studies and find that the GW signal has a generic shape akin to what is known as a "Type III" signal in the literature. We discuss the detectability of the emitted GWs, showing that the signal-to-noise ratio for current or next-generation interferometer detectors could be high enough to detect such events in our Galaxy. Some of our AIC models form massive quasi-Keplerian accretion disks after bounce. In rapidly differentially rotating models, the disk mass can be as large as ~0.8-Msun. Slowly and/or uniformly rotating models produce much smaller disks. Finally, we find that the postbounce cores of rapidly spinning white dwarfs can reach sufficiently rapid rotation to develop a nonaxisymmetric rotational instability., Comment: Submitted to PRD, 31 pages, 23 figures

An increase in the central density of a neutron star may trigger a phase transition from hadronic matter to deconfined quark matter in the core, causing it to collapse to a more compact hybrid-star configuration. We present a study of this, building on previous work by Lin et al. (2006). We follow them in considering a supersonic phase transition and using a simplified equation of state, but our calculations are general relativistic (using 2D simulations in the conformally flat approximation) as compared with their 3D Newtonian treatment. We also improved the treatment of the initial phase transformation, avoiding the introduction of artificial convection. As before, we find that the emitted gravitational-wave spectrum is dominated by the fundamental quasi-radial and quadrupolar pulsation modes but the strain amplitudes are much smaller than suggested previously, which is disappointing for the detection prospects. However, we see significantly smaller damping and observe a nonlinear mode resonance which substantially enhances the emission in some cases. We explain the damping mechanisms operating, giving a different view from the previous work. Finally, we discuss the detectability of the gravitational waves, showing that the signal-to-noise ratio for current or second generation interferometers could be high enough to detect such events in our Galaxy, although third generation detectors would be needed to observe them out to the Virgo cluster, which would be necessary for having a reasonable event rate., Comment: 28 pages, 27 figures. Minor changes to be consistent with published version

We have studied extreme mass-ratio inspirals (EMRIs) in spacetimes containing a rotating black hole and a non self-gravitating torus with constant specific angular momentum. We have found that the effect of the hydrodynamic drag exerted by the torus on the satellite is much smaller than the corresponding one due to radiation reaction, for systems such as those generically expected in AGNs and at distances from the SMBH which can be probed with LISA. However, given the uncertainty on the parameters of these systems, there exist configurations in which the effect of the hydrodynamic drag can be comparable to the radiation-reaction one in phases of the inspiral which are detectable by LISA. This is the case, for instance, for a 10^6 M_sun SMBH surrounded by a corotating torus of comparable mass and with radius of 10^3-10^4 gravitational radii, or for a 10^5 M_sun SMBH surrounded by a corotating 10^4 M_sun torus with radius of 10^5 gravitational radii. Should these conditions be met in astrophysical systems, EMRI-gravitational waves could provide a characteristic signature of the presence of the torus. In fact, while radiation reaction always increases the inclination of the orbit with respect to the equatorial plane, the hydrodynamic drag from a torus corotating with the SMBH always decreases it. However, even when initially dominating over radiation reaction, the influence of the hydrodynamic drag decays very rapidly as the satellite moves into the very strong-field region of the SMBH (i.e., p <~ 5M). Although our results have been obtained for a specific class of tori, we argue that they will be qualitatively valid also for more generic distributions of the specific angular momentum., Comment: 21 pages, 8 figures. Phys. Rev. D, in press. Typos fixed and minor changes; added discussion on the gravitational pull exerted by the torus

During the first 40 s after their birth, proto-neutron stars are expected to be subject to at least two types of instability: the convective instability and the neutron-finger one. Both instabilities involve convective motions and hence can trigger dynamo actions which may be responsible for the large magnetic fields in neutron stars and magnetars. We have solved the mean-field induction equation in a simplified one-dimensional model of both the convective and the neutron-finger instability zones. Although very idealized, the model includes the nonlinearities introduced by the feedback processes which tend to saturate the growth of the magnetic field (alpha-quenching) and suppress its turbulent diffusion (eta-quenching). The possibility of a dynamo action is studied within a dynamical model of turbulent diffusivity where the boundary of the unstable zone is allowed to move. We show that the dynamo action can be operative and that the amplification of the magnetic field can still be very effective. Furthermore, we confirm the existence of a critical spin-period, below which the dynamo is always excited independently of the degree of differential rotation, and whose value is related to the size of the neutron-finger instability zone. Finally we provide a relation for the intensity of the final field as a function of the spin of the star and of its differential rotation. Although they were obtained by using a toy model, we expect that our results are able to capture the qualitative and asymptotic behaviour of a mean-field dynamo action developing in the neutron-finger instability zone. Overall, we find that such a dynamo is very efficient in producing magnetic fields well above equipartition and thus that it could represent a possible explanation for the large surface magnetic fields observed in neutron stars., Comment: Accepted for publication in A&A - 10 pages; corrections after language editing included

Gravitational waveforms from the inspiral and ring-down stages of the binary black hole coalescences can be modelled accurately by approximation/perturbation techniques in general relativity. Recent progress in numerical relativity has enabled us to model also the non-perturbative merger phase of the binary black-hole coalescence problem. This enables us to \emph{coherently} search for all three stages of the coalescence of non-spinning binary black holes using a single template bank. Taking our motivation from these results, we propose a family of template waveforms which can model the inspiral, merger, and ring-down stages of the coalescence of non-spinning binary black holes that follow quasi-circular inspiral. This two-dimensional template family is explicitly parametrized by the physical parameters of the binary. We show that the template family is not only \emph{effectual} in detecting the signals from black hole coalescences, but also \emph{faithful} in estimating the parameters of the binary. We compare the sensitivity of a search (in the context of different ground-based interferometers) using all three stages of the black hole coalescence with other template-based searches which look for individual stages separately. We find that the proposed search is significantly more sensitive than other template-based searches for a substantial mass-range, potentially bringing about remarkable improvement in the event-rate of ground-based interferometers. As part of this work, we also prescribe a general procedure to construct interpolated template banks using non-spinning black hole waveforms produced by numerical relativity., Comment: A typo fixed in Eq.(B11)

We present new results on dynamical instabilities in rapidly rotating neutron-stars. In particular, using numerical simulations in full General Relativity, we analyse the effects that the stellar compactness has on the threshold for the onset of the dynamical bar-mode instability, as well as on the appearance of other dynamical instabilities. By using an extrapolation technique developed and tested in our previous study [1], we explicitly determine the threshold for a wide range of compactnesses using four sequences of models of constant baryonic mass comprising a total of 59 stellar models. Our calculation of the threshold is in good agreement with the Newtonian prediction and improves the previous post-Newtonian estimates. In addition, we find that for stars with sufficiently large mass and compactness, the m=3 deformation is the fastest growing one. For all of the models considered, the non-axisymmetric instability is suppressed on a dynamical timescale with an m=1 deformation dominating the final stages of the instability. These results, together with those presented in [1], suggest that an m=1 deformation represents a general and late-time feature of non-axisymmetric dynamical instabilities both in full General Relativity and in Newtonian gravity., Comment: To appear on CQG, NFNR special issue. 16 pages, 5 color figures, movies from http://www.fis.unipr.it/numrel/BarMode/ResearchBarMode.html

Recent progress in numerical relativity has enabled us to model the non-perturbative merger phase of the binary black-hole coalescence problem. Based on these results, we propose a phenomenological family of waveforms which can model the inspiral, merger, and ring-down stages of black hole coalescence. We also construct a template bank using this family of waveforms and discuss its implementation in the search for signatures of gravitational waves produced by black-hole coalescences in the data of ground-based interferometers. This template bank might enable us to extend the present inspiral searches to higher-mass binary black-hole systems, i.e., systems with total mass greater than about 80 solar masses, thereby increasing the reach of the current generation of ground-based detectors., Comment: Minor changes, Submitted to Class. Quantum Grav. (Proc. GWDAW11)

The accurate modelling of astrophysical scenarios involving compact objects and magnetic fields, such as the collapse of rotating magnetized stars to black holes or the phenomenology of gamma-ray bursts, requires the solution of the Einstein equations together with those of general-relativistic magnetohydrodynamics. We present a new numerical code developed to solve the full set of general-relativistic magnetohydrodynamics equations in a dynamical and arbitrary spacetime with high-resolution shock-capturing techniques on domains with adaptive mesh refinements. After a discussion of the equations solved and of the techniques employed, we present a series of testbeds carried out to validate the code and assess its accuracy. Such tests range from the solution of relativistic Riemann problems in flat spacetime, over to the stationary accretion onto a Schwarzschild black hole and up to the evolution of oscillating magnetized stars in equilibrium and constructed as consistent solutions of the coupled Einstein-Maxwell equations., Comment: minor changes to match the published version

We discuss the procedure for the exact solution of the Riemann problem in special relativistic magnetohydrodynamics (MHD). We consider both initial states leading to a set of only three waves analogous to the ones in relativistic hydrodynamics, as well as generic initial states leading to the full set of seven MHD waves. Because of its generality, the solution presented here could serve as an important test for those numerical codes solving the MHD equations in relativistic regimes., Comment: 36 pages, 13 figures. Minor changes to match published version

We present the first calculation of gravitational wave emission produced in the gravitational collapse of uniformly rotating neutron stars to black holes in fully three-dimensional simulations. The initial stellar models are relativistic polytropes which are dynamically unstable and with angular velocities ranging from slow rotation to the mass-shedding limit. An essential aspect of these simulations is the use of progressive mesh-refinement techniques which allow to move the outer boundaries of the computational domain to regions where gravitational radiation attains its asymptotic form. The waveforms have been extracted using a gauge-invariant approach in which the numerical spacetime is matched with the non-spherical perturbations of a Schwarzschild spacetime. Overall, the results indicate that the waveforms have features related to the properties of the initial stellar models (in terms of their w-mode oscillations) and of the newly produced rotating black holes (in terms of their quasi-normal modes). While our waveforms are in good qualitative agreement with those computed by Stark and Piran in two-dimensional simulations, our amplitudes are about one order of magnitude smaller and this difference is mostly likely due to our less severe pressure reduction. For a neutron star rotating uniformly near mass-shedding and collapsing at 10 kpc, the signal-to-noise ratio computed uniquely from the burst is S/N ~ 0.25, but this grows to be S/N <~ 4 in the case of LIGO II., Comment: 4 pages, Phys. Rev. Lett., in press

We present a comprehensive numerical study of the dynamics of relativistic axisymmetric accretion tori with a power-law distribution of specific angular momentum orbiting in the background spacetime of a Kerr black hole. By combining general relativistic hydrodynamics simulations with a linear perturbative approach we investigate the main dynamical properties of these objects over a large parameter space. The astrophysical implications of our results extend and improve two interesting results that have been recently reported in the literature. Firstly, the induced quasi-periodic variation of the mass quadrupole moment makes relativistic tori of nuclear matter densities, as those formed during the last stages of binary neutron star mergers, promising sources of gravitational radiation, potentially detectable by interferometric instruments. Secondly, $p$-mode oscillations in relativistic tori of low rest-mass densities could be used to explain high frequency quasi-periodic oscillations observed in X-ray binaries containing a black hole candidate under conditions more generic than those considered so far., Comment: 12 pages, 7 figures, accepted for publication by MNRAS

We present a new three-dimensional fully general-relativistic hydrodynamics code using high-resolution shock-capturing techniques and a conformal traceless formulation of the Einstein equations. Besides presenting a thorough set of tests which the code passes with very high accuracy, we discuss its application to the study of the gravitational collapse of uniformly rotating neutron stars to Kerr black holes. The initial stellar models are modelled as relativistic polytropes which are either secularly or dynamically unstable and with angular velocities which range from slow rotation to the mass-shedding limit. We investigate the gravitational collapse by carefully studying not only the dynamics of the matter, but also that of the trapped surfaces, i.e. of both the apparent and event horizons formed during the collapse. The use of these surfaces, together with the dynamical horizon framework, allows for a precise measurement of the black-hole mass and spin. The ability to successfully perform these simulations for sufficiently long times relies on excising a region of the computational domain which includes the singularity and is within the apparent horizon. The dynamics of the collapsing matter is strongly influenced by the initial amount of angular momentum in the progenitor star and, for initial models with sufficiently high angular velocities, the collapse can lead to the formation of an unstable disc in differential rotation. All the simulations performed with uniformly rotating initial data and a polytropic or ideal-fluid equation of state show no evidence of shocks or of the presence of matter on stable orbits outside the black hole., Comment: 28 pages. Animations can be found at http://www.sissa.it/~rezzolla/Whisky/WhiskyI/

Observations of X-ray emissions from binary systems have always been considered important tools to test the validity of General Relativity in strong-field regimes. The pairs and triplets of high frequency quasi-periodic oscillations observed in binaries containing a black hole candidate, in particular, have been proposed as a means to measure more directly the black hole properties such as its mass and spin. Numerous models have been suggested over the years to explain the QPOs and the rich phenomenology accompanying them. Many of these models rest on a number of assumptions and are at times in conflict with the most recent observations. We here propose a new, simple model in which the QPOs result from basic p-mode oscillations of a non-Keplerian disc of finite size. We show that within this new model all of the key properties of the QPOs: a) harmonic ratios of frequencies even as the frequencies change; b) variations in the relative strength of the frequencies with spectral energy distribution and with photon energy; c) small and systematic changes in the frequencies, can all be explained simply given a single reasonable assumption., Comment: To appear in Proceedings of "X-ray Timing 2003: Rossi and Beyond", ed. P. Kaaret, F.K. Lamb, & J.H. Swank

We present general relativistic hydrodynamics simulations of constant specific angular momentum tori orbiting a Schwarzschild black hole. These tori are expected to form as a result of stellar gravitational collapse, binary neutron star merger or disruption, can reach very high rest-mass densities and behave effectively as neutron stars but with a toroidal topology (i.e. ``toroidal neutron stars''). Our attention is here focussed on the dynamical response of these objects to axisymmetric perturbations. We show that, upon the introduction of perturbations, these systems either become unstable to the runaway instability or exhibit a regular oscillatory behaviour resulting in a quasi-periodic variation of the accretion rate as well as of the mass quadrupole. The latter, in particular, is responsible for the emission of intense gravitational radiation whose signal-to-noise ratio at the detector is comparable or larger than the typical one expected in stellar-core collapse, making these new sources of gravitational waves potentially detectable. We discuss a systematic investigation of the parameter space both in the linear and nonlinear regimes, providing estimates of how the gravitational radiation emitted depends on the mass of the torus and on the strength of the perturbation., Comment: 17 pages, 12 figures. Animations can be found at http://www.sissa.it/~rezzolla/movies.html

We extend our approach for the exact solution of the Riemann problem in relativistic hydrodynamics to the case in which the fluid velocity has components tangential to the initial discontinuity. As in one-dimensional flows, we here show that the wave-pattern produced in a multidimensional relativistic Riemann problem can be predicted entirely by examining the initial conditions. Our method is logically very simple and allows for a numerical implementation of an exact Riemann solver which is both straightforward and computationally efficient. The simplicity of the approach is also important for revealing special relativistic effects responsible for a smooth transition from one wave-pattern to another when the tangential velocities in the initial states are suitably varied. While the content of this paper is focussed on a flat spacetime, the local Lorentz invariance allows its use also in fully general relativistic calculations., Comment: 20 pages, 8 figures

In Newtonian and relativistic hydrodynamics the Riemann problem consists of calculating the evolution of a fluid which is initially characterized by two states having different values of uniform rest-mass density, pressure and velocity. When the fluid is allowed to relax, one of three possible wave-patterns is produced, corresponding to the propagation in opposite directions of two nonlinear hydrodynamical waves. New effects emerge in a special relativistic Riemann problem when velocities tangential to the initial discontinuity surface are present. We show that a smooth transition from one wave-pattern to another can be produced by varying the initial tangential velocities while otherwise maintaining the initial states unmodified. These special relativistic effects are produced by the coupling through the relativistic Lorentz factors and do not have a Newtonian counterpart., Comment: 4 pages, 5 figures

We have solved numerically the general relativistic induction equations in the interior background spacetime of a slowly rotating magnetized neutron star. The analytic form of these equations was discussed in a recent paper (Rezzolla et al 2001a), where corrections due both to the spacetime curvature and to the dragging of reference frames were shown to be present. Through a number of calculations we have investigated the evolution of the magnetic field with different rates of stellar rotation, different inclination angles between the magnetic moment and the rotation axis, as well as different values of the electrical conductivity. All of these calculations have been performed for a constant temperature relativistic polytropic star and make use of a consistent solution of the initial value problem which avoids the use of artificial analytic functions. Our results show that there exist general relativistic effects introduced by the rotation of the spacetime which tend to decrease the decay rate of the magnetic field. The rotation-induced corrections are however generally hidden by the high electrical conductivity of the neutron star matter and when realistic values for the electrical conductivity are considered, these corrections become negligible even for the fastest known pulsar., Comment: 13 pages, 5 figures. Accepted for publication by MNRAS. Replaces previous version without unnecessary mn.sty

This is the second in a series of papers on the construction and validation of a three-dimensional code for the solution of the coupled system of the Einstein equations and of the general relativistic hydrodynamic equations, and on the application of this code to problems in general relativistic astrophysics. In particular, we report on the accuracy of our code in the long-term dynamical evolution of relativistic stars and on some new physics results obtained in the process of code testing. The tests involve single non-rotating stars in stable equilibrium, non-rotating stars undergoing radial and quadrupolar oscillations, non-rotating stars on the unstable branch of the equilibrium configurations migrating to the stable branch, non-rotating stars undergoing gravitational collapse to a black hole, and rapidly rotating stars in stable equilibrium and undergoing quasi-radial oscillations. The numerical evolutions have been carried out in full general relativity using different types of polytropic equations of state using either the rest-mass density only, or the rest-mass density and the internal energy as independent variables. New variants of the spacetime evolution and new high resolution shock capturing (HRSC) treatments based on Riemann solvers and slope limiters have been implemented and the results compared with those obtained from previous methods. Finally, we have obtained the first eigenfrequencies of rotating stars in full general relativity and rapid rotation. A long standing problem, such frequencies have not been obtained by other methods. Overall, and to the best of our knowledge, the results presented in this paper represent the most accurate long-term three-dimensional evolutions of relativistic stars available to date., Comment: 19 pages, 17 figures

Following the general formalism presented by Rezzolla, Ahmedov and Miller (MNRAS, 322, 723 2001), we here derive analytic solutions of the electromagnetic fields equations in the internal and external background spacetime of a slowly rotating highly conducting magnetized neutron star. The star is assumed to be isolated and in vacuum, with a dipolar magnetic field not aligned with the axis of rotation. Our results indicate that the electromagnetic fields of a slowly rotating neutron star are modified by general relativistic effects arising from both the monopolar and the dipolar parts of the gravitational field. The results presented here differ from the ones discussed by Rezzolla, Ahmedov and Miller (MNRAS, 322, 723 2001) mainly in that we here consider the interior magnetic field to be dipolar with the same radial dependence as the external one. While this assumption might not be a realistic one, it should be seen as the application of our formalism to a case often discussed in the literature., Comment: 12 pages, no figures

A Riemann problem with prescribed initial conditions will produce one of three possible wave patterns corresponding to the propagation of the different discontinuities that will be produced once the system is allowed to relax. In general, when solving the Riemann problem numerically, the determination of the specific wave pattern produced is obtained through some initial guess which can be successively discarded or improved. We here discuss a new procedure, suitable for implementation in an exact Riemann solver in one dimension, which removes the initial ambiguity in the wave pattern. In particular we focus our attention on the relativistic velocity jump between the two initial states and use this to determine, through some analytic conditions, the wave pattern produced by the decay of the initial discontinuity. The exact Riemann problem is then solved by means of calculating the root of a nonlinear equation. Interestingly, in the case of two rarefaction waves, this root can even be found analytically. Our procedure is straightforward to implement numerically and improves the efficiency of numerical codes based on exact Riemann solvers., Comment: Extended the discussion on the computational efficiency with improved estimates; matches the version in press

We present analytic solutions of Maxwell equations in the internal and external background spacetime of a slowly rotating magnetized neutron star. The star is considered isolated and in vacuum, with a dipolar magnetic field not aligned with the axis of rotation. With respect to a flat spacetime solution, general relativity introduces corrections related both to the monopolar and the dipolar parts of the gravitational field. In particular, we show that in the case of infinite electrical conductivity general relativistic corrections due to the dragging of reference frames are present, but only in the expression for the electric field. In the case of finite electrical conductivity, however, corrections due both to the spacetime curvature and to the dragging of reference frames are shown to be present in the induction equation. These corrections could be relevant for the evolution of the magnetic fields of pulsars and magnetars. The solutions found, while obtained through some simplifying assumption, reflect a rather general physical configuration and could therefore be used in a variety of astrophysical situations., Comment: A few typos corrected; matches the versions in MNRAS

We report on our numerical implementation of fully relativistic hydrodynamics coupled to Einstein's field equations in three spatial dimensions. We briefly review several steps in our code development, including our recasting of Einstein's equations and several tests which demonstrate its advantages for numerical integrations. We outline our implementation of relativistic hydrodynamics, and present numerical results for the evolution of both stable and unstable Oppenheimer-Volkov equilibrium stars, which represent a very promising first test of our code., Comment: 5 Pages, 4 Figures, submitted to Proceedings of the 8th Canadian Conference on General Relativity and Relativistic Astrophysics

We present a new set of 3.5 Post-Newtonian equations in which Newtonian hydrodynamics is coupled to the nonconservative effects of gravitational radiation emission. Our formalism differs in two significant ways from a similar 3.5 Post-Newtonian approach proposed by Blanchet (1993, 1997). Firstly we concentrate only on the radiation-reaction effects produced by a time-varying mass-current quadrupole $S_{ij}$. Secondly, we adopt a gauge in which the radiation-reaction force densities depend on the fourth time derivative of $S_{ij}$, rather than on the fifth, as in Blanchet's approach. This difference makes our formalism particularly well-suited to numerical implementation and could prove useful in performing fully numerical simulations of the recently discovered $r$-mode instability for rotating neutron stars subject to axial perturbations., Comment: 17 pages, Latex2e, emulateapj.sty. To appear in ApJ

We present a new method of extracting gravitational radiation from three-dimensional numerical relativity codes and providing outer boundary conditions. Our approach matches the solution of a Cauchy evolution of Einstein's equations to a set of one-dimensional linear wave equations on a curved background. We illustrate the mathematical properties of our approach and discuss a numerical module we have constructed for this purpose. This module implements the perturbative matching approach in connection with a generic three-dimensional numerical relativity simulation. Tests of its accuracy and second-order convergence are presented with analytic linear wave data., Comment: 13 pages, 6 figures, RevTeX

We report new results which establish that the accurate 3-dimensional numerical simulation of generic single-black-hole spacetimes has been achieved by characteristic evolution with unlimited long term stability. Our results cover a selection of distorted, moving and spinning single black holes, with evolution times up to 60,000M., Comment: 4 pages, 3 figures

We compute the propagation and scattering of linear gravitational waves off a Schwarzschild black hole using a numerical code which solves a generalization of the Zerilli equation to a three dimensional cartesian coordinate system. Since the solution to this problem is well understood it represents a very good testbed for evaluating our ability to perform three dimensional computations of gravitational waves in spacetimes in which a black hole event horizon is present., Comment: 13 pages, RevTeX, to appear in Phys. Rev. D

Binary black hole interactions provide potentially the strongest source of gravitational radiation for detectors currently under development. We present some results from the Binary Black Hole Grand Challenge Alliance three- dimensional Cauchy evolution module. These constitute essential steps towards modeling such interactions and predicting gravitational radiation waveforms. We report on single black hole evolutions and the first successful demonstration of a black hole moving freely through a three-dimensional computational grid via a Cauchy evolution: a hole moving ~6M at 0.1c during a total evolution of duration ~60M.