High--Resolution 3D Simulations of Relativistic Jets

We have performed high-resolution 3D simulations of relativistic jets with beam flow Lorentz factors up to 7, a spatial resolution of 8 cells per beam radius, and for up to 75 normalized time units to study the morphology and dynamics of 3D relativistic jets. Our simulations show that the coherent fast backflows found in axisymmetric models are not present in 3D models. We further find that when the jet is exposed to non-axisymmetric perturbations, (i) it does not display the strong perturbations found for 3D classical hydrodynamic and MHD jets (at least during the period of time covered by our simulations), and (ii) it does propagate according to the 1D estimate. Small 3D effects in the relativistic beam give rise to a lumpy distribution of apparent speeds like that observed in M87. The beam is surrounded by a boundary layer of high specific internal energy. The properties of this layer are briefly discussed.Comment: 15 pages, 4 figures. Accepted to be publish in the ApJ Letters. Tar+gzip documen

Introducing PHAEDRA: a new spectral code for simulations of relativistic magnetospheres

We describe a new scheme for evolving the equations of force-free electrodynamics, the vanishing-inertia limit of magnetohydrodynamics. This pseudospectral code uses global orthogonal basis function expansions to take accurate spatial derivatives, allowing the use of an unstaggered mesh and the complete force-free current density. The method has low numerical dissipation and diffusion outside of singular current sheets. We present a range of one- and two-dimensional tests, and demonstrate convergence to both smooth and discontinuous analytic solutions. As a first application, we revisit the aligned rotator problem, obtaining a steady solution with resistivity localised in the equatorial current sheet outside the light cylinder.Comment: 23 pages, 18 figures, accepted for publication in MNRA

The pulsar force-free magnetosphere linked to its striped wind: time-dependent pseudo-spectral simulations

(abridged) Pulsar activity and its related radiation mechanism are usually explained by invoking some plasma processes occurring inside the magnetosphere. Despite many detailed local investigations, the global electrodynamics around those neutron stars remains poorly described. Better understanding of these compact objects requires a deep and accurate knowledge of their immediate electromagnetic surrounding within the magnetosphere and its link to the relativistic pulsar wind. The aim of this work is to present accurate solutions to the nearly stationary force-free pulsar magnetosphere and its link to the striped wind, for various spin periods and arbitrary inclination. To this end, the time-dependent Maxwell equations are solved in spherical geometry in the force-free approximation using a vector spherical harmonic expansion of the electromagnetic field. An exact analytical enforcement of the divergenceless of the magnetic part is obtained by a projection method. Special care has been given to design an algorithm able to look deeply into the magnetosphere with physically realistic ratios of stellar $R_*$ to light-cylinder \rlight radius. We checked our code against several analytical solutions, like the Deutsch vacuum rotator solution and the Michel monopole field. We also retrieve energy losses comparable to the magneto-dipole radiation formula and consistent with previous similar works. Finally, for arbitrary obliquity, we give an expression for the total electric charge of the system. It does not vanish except for the perpendicular rotator. This is due to the often ignored point charge located at the centre of the neutron star. It is questionable if such solutions with huge electric charges could exist in reality except for configurations close to an orthogonal rotator. The charge spread over the stellar crust is not a tunable parameter as is often hypothesized.Comment: 16 pages, 13 figures, accepted by MNRA

Foundations of Black Hole Accretion Disk Theory

This review covers the main aspects of black hole accretion disk theory. We begin with the view that one of the main goals of the theory is to better understand the nature of black holes themselves. In this light we discuss how accretion disks might reveal some of the unique signatures of strong gravity: the event horizon, the innermost stable circular orbit, and the ergosphere. We then review, from a first-principles perspective, the physical processes at play in accretion disks. This leads us to the four primary accretion disk models that we review: Polish doughnuts (thick disks), Shakura-Sunyaev (thin) disks, slim disks, and advection-dominated accretion flows (ADAFs). After presenting the models we discuss issues of stability, oscillations, and jets. Following our review of the analytic work, we take a parallel approach in reviewing numerical studies of black hole accretion disks. We finish with a few select applications that highlight particular astrophysical applications: measurements of black hole mass and spin, black hole vs. neutron star accretion disks, black hole accretion disk spectral states, and quasi-periodic oscillations (QPOs).Comment: 91 pages, 23 figures, final published version available at http://www.livingreviews.org/lrr-2013-

Lepton Acceleration in Pulsar Wind Nebulae

Pulsar Wind Nebulae (PWNe) act as calorimeters for the relativistic pair winds emanating from within the pulsar light cylinder. Their radiative dissipation in various wavebands is significantly different from that of their pulsar central engines: the broadband spectra of PWNe possess characteristics distinct from those of pulsars, thereby demanding a site of lepton acceleration remote from the pulsar magnetosphere. A principal candidate for this locale is the pulsar wind termination shock, a putatively highly-oblique, ultra-relativistic MHD discontinuity. This paper summarizes key characteristics of relativistic shock acceleration germane to PWNe, using predominantly Monte Carlo simulation techniques that compare well with semi-analytic solutions of the diffusion-convection equation. The array of potential spectral indices for the pair distribution function is explored, defining how these depend critically on the parameters of the turbulent plasma in the shock environs. Injection efficiencies into the acceleration process are also addressed. Informative constraints on the frequency of particle scattering and the level of field turbulence are identified using the multiwavelength observations of selected PWNe. These suggest that the termination shock can be comfortably invoked as a principal injector of energetic leptons into PWNe without resorting to unrealistic properties for the shock layer turbulence or MHD structure.Comment: 19 pages, 5 figures, invited review to appear in Proc. of the inaugural ICREA Workshop on "The High-Energy Emission from Pulsars and their Systems" (2010), eds. N. Rea and D. Torres, (Springer Astrophysics and Space Science series

A search for steep spectrum radio relics and halos with the GMRT

Context: Diffuse radio emission, in the form of radio halos and relics, traces regions in clusters with shocks or turbulence, probably produced by cluster mergers. Some models of diffuse radio emission in clusters indicate that virtually all clusters should contain diffuse radio sources with a steep spectrum. External accretion shocks associated with filamentary structures of galaxies could also accelerate electrons to relativistic energies and hence produce diffuse synchrotron emitting regions. Here we report on Giant Metrewave Radio Telescope (GMRT) observations of a sample of steep spectrum sources from the 74 MHz VLSS survey. These sources are diffuse and not associated with nearby galaxies. Aims: The main aim of the observations is to search for diffuse radio emission associated with galaxy clusters or the cosmic web. Methods: We carried out GMRT 610 MHz continuum observations of unidentified diffuse steep spectrum sources. Results: We have constructed a sample of diffuse steep spectrum sources, selected from the 74 MHz VLSS survey. We identified eight diffuse radio sources probably all located in clusters. We found five radio relics, one cluster with a giant radio halo and a radio relic, and one radio mini-halo. By complementing our observations with measurements from the literature we find correlations between the physical size of relics and the spectral index, in the sense that smaller relics have steeper spectra. Furthermore, larger relics are mostly located in the outskirts of clusters while smaller relics are located closer to the cluster center.Comment: 20 pages, 26 figures, accepted for publication in A&A on October 7, 200

Three little pieces for computer and relativity

Numerical relativity has made big strides over the last decade. A number of problems that have plagued the field for years have now been mostly solved. This progress has transformed numerical relativity into a powerful tool to explore fundamental problems in physics and astrophysics, and I present here three representative examples. These "three little pieces" reflect a personal choice and describe work that I am particularly familiar with. However, many more examples could be made.Comment: 42 pages, 11 figures. Plenary talk at "Relativity and Gravitation: 100 Years after Einstein in Prague", June 25 - 29, 2012, Prague, Czech Republic. To appear in the Proceedings (Edition Open Access). Collects results appeared in journal articles [72,73, 122-124

Magnetic Reconnection in Extreme Astrophysical Environments

Magnetic reconnection is a basic plasma process of dramatic rearrangement of magnetic topology, often leading to a violent release of magnetic energy. It is important in magnetic fusion and in space and solar physics --- areas that have so far provided the context for most of reconnection research. Importantly, these environments consist just of electrons and ions and the dissipated energy always stays with the plasma. In contrast, in this paper I introduce a new direction of research, motivated by several important problems in high-energy astrophysics --- reconnection in high energy density (HED) radiative plasmas, where radiation pressure and radiative cooling become dominant factors in the pressure and energy balance. I identify the key processes distinguishing HED reconnection: special-relativistic effects; radiative effects (radiative cooling, radiation pressure, and Compton resistivity); and, at the most extreme end, QED effects, including pair creation. I then discuss the main astrophysical applications --- situations with magnetar-strength fields (exceeding the quantum critical field of about 4 x 10^13 G): giant SGR flares and magnetically-powered central engines and jets of GRBs. Here, magnetic energy density is so high that its dissipation heats the plasma to MeV temperatures. Electron-positron pairs are then copiously produced, making the reconnection layer highly collisional and dressing it in a thick pair coat that traps radiation. The pressure is dominated by radiation and pairs. Yet, radiation diffusion across the layer may be faster than the global Alfv\'en transit time; then, radiative cooling governs the thermodynamics and reconnection becomes a radiative transfer problem, greatly affected by the ultra-strong magnetic field. This overall picture is very different from our traditional picture of reconnection and thus represents a new frontier in reconnection research.Comment: Accepted to Space Science Reviews (special issue on magnetic reconnection). Article is based on an invited review talk at the Yosemite-2010 Workshop on Magnetic Reconnection (Yosemite NP, CA, USA; February 8-12, 2010). 30 pages, no figure

Relics as probes of galaxy cluster mergers

Galaxy clusters grow by mergers with other clusters and galaxy groups. These mergers create shocks within the intracluster medium (ICM). It is proposed that within the shocks particles can be accelerated to extreme energies. In the presence of a magnetic field these particles should then form large regions emitting synchrotron radiation, creating so-called radio relics. An example of a cluster with relics is CIZA J2242.8+5301. Here we present hydrodynamical simulations of idealized binary cluster collisions with the aim of constraining the merger scenario for this cluster. We conclude that by using the location, size and width of double radio relics we can set constraints on the mass ratios, impact parameters, timescales, and viewing geometries of binary cluster merger events.Comment: Accepted for publication in special issue of Journal of Astrophysics and Astronomy: conference proceedings of "Diffuse Relativistic Plasmas" conference, Bangalore, 1-4 March 2011, 4 pages, 2 figure

Paracrine cyclooxygenase-2 activity by macrophages drives colorectal adenoma progression in the Apc Min/+ mouse model of intestinal tumorigenesis

Genetic deletion or pharmacological inhibition of cyclooxygenase (COX)-2 abrogates intestinal adenoma development at early stages of colorectal carcinogenesis. COX-2 is localised to stromal cells (predominantly macrophages) in human and mouse intestinal adenomas. Therefore, we tested the hypothesis that paracrine Cox-2-mediated signalling from macrophages drives adenoma growth and progression in vivo in the ApcMin/+ mouse model of intestinal tumorigenesis. Using a transgenic C57Bl/6 mouse model of Cox-2 over-expression driven by the chicken lysozyme locus (cLys-Cox-2), which directs integration site-independent, copy number-dependent transgene expression restricted to macrophages, we demonstrated that stromal macrophage Cox-2 in colorectal (but not small intestinal) adenomas from cLys-Cox-2 x ApcMin/+ mice was associated with significantly increased tumour size (P = 0.025) and multiplicity (P = 0.025), compared with control ApcMin/+ mice. Transgenic macrophage Cox-2 expression was associated with increased dysplasia, epithelial cell Cox-2 expression and submucosal tumour invasion, as well as increased nuclear β-catenin translocation in dysplastic epithelial cells. In vitro studies confirmed that paracrine macrophage Cox-2 signalling drives catenin-related transcription in intestinal epithelial cells. Paracrine macrophage Cox-2 activity drives growth and progression of ApcMin/+ mouse colonic adenomas, linked to increased epithelial cell β-catenin dysregulation. Stromal cell (macrophage) gene regulation and signalling represent valid targets for chemoprevention of colorectal cancer