Researcher profile

Philip E. Hardee

Philip E. Hardee contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
13works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2016arXiv

Microscopic Processes in Global Relativistic Jets Containing Helical Magnetic Fields

In the study of relativistic jets one of the key open questions is their interaction with the environment on the microscopic level. Here, we study the initial evolution of both electron$-$proton ($e^{-}-p^{+}$) and electron$-$positron ($e^{\pm}$) relativistic jets containing helical magnetic fields, focusing on their interaction with an ambient plasma. We have performed simulations of "global" jets containing helical magnetic fields in order to examine how helical magnetic fields affect kinetic instabilities such as the Weibel instability, the kinetic Kelvin-Helmholtz instability (kKHI) and the Mushroom instability (MI). In our initial simulation study these kinetic instabilities are suppressed and new types of instabilities can grow. In the $e^{-}-p^{+}$ jet simulation a recollimation-like instability occurs and jet electrons are strongly perturbed. In the $e^{\pm}$ jet simulation a recollimation-like instability occurs at early times followed by a kinetic instability and the general structure is similar to a simulation without helical magnetic field. Simulations using much larger systems are required in order to thoroughly follow the evolution of global jets containing helical magnetic fields.

preprint2016arXiv

Observations of the Structure and Dynamics of the Inner M87 Jet

M87 is the best source in which to study a jet at high resolution in gravitational units because it has a very high mass black hole and is nearby. The angular size of the black hole is second only to Sgr A*, which does not have a strong jet. The jet structure is edge brightened with a wide opening angle base and a weak counterjet. We have roughly annual observations for 17 years plus intensive monitoring at three week intervals for a year and five day intervals for 2.5 months made with the Very Long Baseline Array (VLBA) at 43 GHz. The inner jet shows very complex dynamics, with apparent motions both along and across the jet. Speeds from zero to over 2c are seen, with acceleration observed over the first 3 milli-arcseconds. The counterjet decreases in brightness much more rapidly than the main jet, as is expected from relativistic beaming in an accelerating jet oriented near the line-of-sight. Details of the structure and dynamics are discussed. The roughly annual observations show side-to-side motion of the whole jet with a characteristic time scale of about 9 years.

preprint2015arXiv

Recollimation Shocks in Magnetized Relativistic Jets

We have performed two-dimensional special-relativistic magnetohydrodynamic simulations of non-equilibrium over-pressured relativistic jets in cylindrical geometry. Multiple stationary recollimation shock and rarefaction structures are produced along the jet by the nonlinear interaction of shocks and rarefaction waves excited at the interface between the jet and the surrounding ambient medium. Although initially the jet is kinematically dominated, we have considered axial, toroidal and helical magnetic fields to investigate the effects of different magnetic-field topologies and strengths on the recollimation structures. We find that an axial field introduces a larger effective gas-pressure and leads to stronger recollimation shocks and rarefactions, resulting in larger flow variations. The jet boost grows quadratically with the initial magnetic field. On the other hand, a toroidal field leads to weaker recollimation shocks and rarefactions, modifying significantly the jet structure after the first recollimation rarefaction and shock. The jet boost decreases systematically. For a helical field, instead, the behaviour depends on the magnetic pitch, with a phenomenology that ranges between the one seen for axial and toroidal magnetic fields, respectively. In general, however, a helical magnetic field yields a more complex shock and rarefaction substructure close to the inlet that significantly modifies the jet structure. The differences in shock structure resulting from different field configurations and strengths may have observable consequences for disturbances propagating through a stationary recollimation shock.

preprint2014arXiv

Magnetic Field Amplification and Saturation in Turbulence Behind a Relativistic Shock

We have investigated via two-dimensional relativistic MHD simulations the long-term evolution of turbulence created by a relativistic shock propagating through an inhomogeneous medium. In the postshock region, magnetic field is strongly amplified by turbulent motions triggered by preshock density inhomogeneities. Using a long-simulation box we have followed the magnetic-field amplification until it is fully developed and saturated. The turbulent velocity is sub-relativistic even for a strong shock. Magnetic-field amplification is controled by the turbulent motion and saturation occurs when the magnetic energy is comparable to the turbulent kinetic energy. Magnetic-field amplification and saturation depend on the initial strength and direction of the magnetic field in the preshock medium, and on the shock strength. If the initial magnetic field is perpendicular to the shock normal, the magnetic field is first compressed at the shock and then can be amplified by turbulent motion in the postshock region. Saturation occurs when the magnetic energy becomes comparable to the turbulent kinetic energy in the postshock region. If the initial magnetic field in the preshock medium is strong, the postshock region becomes turbulent but significant field amplification does not occur. If the magnetic energy after shock compression is larger than the turbulent kinetic energy in the postshock region, significant field amplification does not occur. We discuss possible applications of our results to gamma-ray bursts and active galactic nuclei.

preprint2014arXiv

On the deceleration of FRI jets: mass loading by stellar-winds

Jets in low-luminosity radio galaxies are known to decelerate from relativistic speeds on parsec scales to mildly or sub-relativistic speeds on kiloparsec scales. Several mechanisms have been proposed to explain this effect, including strong reconfinement shocks and the growth of instabilities (both leading to boundary-layer entrainment) and mass loading from stellar winds or molecular clouds. We have performed a series of axisymmetric simulations of the early evolution of jets in a realistic ambient medium to probe the effects of mass loading from stellar winds using the code Ratpenat. We study the evolution of Fanaroff-Riley Class I (FRI) jets, with kinetic powers L_j \sim 1.e41-1.e44 erg/s, within the first 2 kpc of their evolution, where deceleration by stellar mass loading should be most effective. Mass entrainment rates consistent with present models of stellar mass loss in elliptical galaxies produce deceleration and effective decollimation of weak FRI jets within the first kiloparsec. However, powerful FRI jets are not decelerated significantly. In those cases where the mass loading is important, the jets show larger opening angles and decollimate at smaller distances, but the overall structure and dynamics of the bow-shock are similar to those of unloaded jets with the same power and thrust. According to our results, the flaring observed on kpc scales is initiated by mass loading in the weaker FRI jets and by reconfinement shocks or the growth of instabilities in the more powerful jets. The final mechanism of decollimation and deceleration is always the development of disruptive pinching modes.

preprint2014arXiv

Spatial Growth of the Current-Driven Instability in Relativistic Jets

We have investigated the influence of velocity shear and a radial density profile on the spatial development of the current driven kink instability along helically magnetized relativistic jets via three-dimensional relativistic magnetohydrodynamic simulations. In this study, we use a non-periodic computational box, the jet flow is initially established across the computational grid, and a precessional perturbation at the inlet triggers growth of the kink instability. If the velocity shear radius is located inside the characteristic radius of the helical magnetic field, a static non-propagating current driven kink is excited as the perturbation propagates down the jet. Temporal growth disrupts the initial flow across the computational grid not too far from the inlet. On the other hand, if the velocity shear radius is outside the characteristic radius of the helical magnetic field, the kink is advected with the flow and grows spatially down the jet. In this case flow is maintained to much larger distances from the inlet. The effect of different radial density profiles is more subtle. When the density increases with radius, the kink appears to saturate by the end of the simulation without apparent disruption of the helical twist. This behavior suggests that relativistic jets consisting of a tenuous spine surrounded by a denser medium with a velocity shear radius outside the radius of maximum toroidal magnetic field have a relatively stable configuration.

preprint2012arXiv

Anatomy of helical relativistic jets: The case of S5 0836+710

Helical structures are common in extragalactic jets. They are usually attributed in the literature to periodical phenomena in the source (e.g., precession). In this work, we use VLBI data of the radio-jet in the quasar S5 0836+710 and hypothesize that the ridge-line of helical jets like this corresponds to a pressure maximum in the jet and assume that the helically twisted pressure maximum is the result of a helical wave pattern. For our study, we use observations of the jet in S5 0836+710 at different frequencies and epochs. The results show that the structures observed are physical and not generated artificially by the observing arrays. Our hypothesis that the observed intensity ridge-line can correspond to a helically twisted pressure maximum is confirmed by our observational tests. This interpretation allows us to explain jet misalignment between parsec and kiloparsec scales when the viewing angle is small, and also brings us to the conclusion that high-frequency observations may show only a small region of the jet flow concentrated around the maximum pressure ridge-line observed at low frequencies. Our work provides a potential explanation for the apparent transversal superluminal speeds observed in several extragalactic jets by means of transversal shift of an apparent core position with time.

preprint2012arXiv

S5 0836+710: An FRII jet disrupted by the growth of a helical instability?

The remarkable stability of extragalactic jets is surprising, given the reasonable possibility of the growth of instabilities. In addition, much work in the literature has invoked this possibility in order to explain observed jet structures and obtain information from these structures. For example, it was recently shown that the observed helical structures in the jet in S5 0836+710 could be associated with helical pressure waves generated by Kelvin-Helmholtz instability. Our aim is to resolve the arc-second structure of the jet in the quasar S5 0836+710 and confirm the lack of a hot-spot (reverse jet-shock) found by present observing arrays, as this lack implies a loss of jet collimation before interaction with the intergalactic medium. In this work, we use an observation performed in 2008 using EVN and MERLIN. The combined data reduction has provided a complete image of the object at arc-second scales. The lack of a hot-spot in the arc-second radio structure is taken as evidence that the jet losses its collimation between the VLBI region and the region of interaction with the ambient medium. This result, together with the previous identification of the helical structures in the jet with helical pressure waves that grow in amplitude with distance, allow us to conclude that the jet is probably disrupted by the growth of Kelvin-Helmholtz instability. This observational evidence confirms that the physical parameters of jets can be extracted using the assumption that instability is present in jets and can be the reason for many observed structures. Interestingly, the observed jet is classified as a FRII object in terms of its luminosity, but its large-scale morphology does not correspond to this classification. The implications of this fact are discussed.

preprint2012arXiv

Three-Dimensional Relativistic Magnetohydrodynamic Simulations of Current-Driven Instability. III. Rotating Relativistic Jets

We have investigated the influence of jet rotation and differential motion on the linear and nonlinear development of the current-driven (CD) kink instability of force-free helical magnetic equilibria via three-dimensional relativistic magnetohydrodynamic simulations. In this study, we follow the temporal development within a periodic computational box. Displacement of the initial helical magnetic field leads to the growth of the CD kink instability. We find that, in accord with linear stability theory, the development of the instability depends on the lateral distribution of the poloidal magnetic field. If the poloidal field significantly decreases outwards from the axis, the initial small perturbations grow strongly, and if multiple wavelengths are excited non-linear interaction eventually disrupts the initial cylindrical configuration. When the profile of the poloidal field is shallow, the instability develops slowly and eventually saturates. We briefly discuss implications of our findings for Poynting dominated jets.

preprint2011arXiv

Three-Dimensional Relativistic Magnetohydrodynamic Simulations of Current-Driven Instability with A Sub-Alfvenic Jet: Temporal Properties

We have investigated the influence of a velocity shear surface on the linear and non-linear development of the CD kink instability of force-free helical magnetic equilibria in 3D. In this study we follow the temporal development within a periodic computational box and concentrate on flows that are sub-Alfvenic on the cylindrical jet's axis. Displacement of the initial force-free helical magnetic field leads to the growth of CD kink instability. We find that helically distorted density structure propagates along the jet with speed and flow structure dependent on the radius of the velocity shear surface relative to the characteristic radius of the helically twisted force-free magnetic field. At small velocity shear surface radius the plasma flows through the kink with minimal kink propagation speed. The kink propagation speed increases as the velocity shear radius increases and the kink becomes more embedded in the plasma flow. A decreasing magnetic pitch profile and faster flow enhance the influence of velocity shear. Simulations show continuous transverse growth in the nonlinear phase of the instability. The growth rate of the CD kink instability and the nonlinear behavior also depend on the velocity shear surface radius and flow speed, and the magnetic pitch radial profile. Larger velocity shear radius leads to slower linear growth, makes a later transition to the nonlinear stage, and with larger maximum amplitude than occur for a static plasma column. However, when the velocity shear radius is much greater than the characteristic radius of the helical magnetic field, linear and non-linear development can be similar to the development of a static plasma column.

preprint2011arXiv

Using Twisted Filaments to Model the Inner Jet in M87

Radio and optical images of the M87 jet show bright filaments, twisted into an apparent double helix, extending from HST-1 to knot A. Proper motions within the jet suggest a decelerating jet flow passing through a slower, accelerating wave pattern. We use these observations to develop a mass and energy flux conserving model describing the jet flow and conditions along the jet. We determine the cocoon conditions required if the twisted filaments are the result of the Kelvin-Helmholtz (KH) unstable elliptical mode. We find that the cocoon must be cooler than the jet at HST-1 but must be about as hot as the jet at knot A. Under these conditions we find that the observed filament wavelength is near the elliptical mode maximum growth rate and growth is rapid enough for the filaments to develop and saturate well before HST-1. We generate a pseudo-synchrotron image of a model jet carrying a combination of normal modes of the KH instability. The pseudo-synchrotron image of the jet reveals: (1) that a slow decline in the model jet's surface brightness is still about five times faster than the real jet; (2) that KH produced dual helically twisted filaments can appear qualitatively similar to those on the real jet if any helical perturbation to the jet is very small or nonexistent inside knot A; (3) that the knots in the real jet cannot be associated with the twisted filamentary features and are unlikely to be the result of a KH instability. The existence of the knots in the real jet, the limb brightening of the real jet in the radio, and the slower decline of the surface brightness of the real jet indicate that additional processes --- such as unsteady jet flow and internal particle acceleration --- are occurring within the jet. Disruption of the real jet beyond knot A by KH instability is consistent with the jet and cocoon conditions we find at knot A.

preprint2010arXiv

Magnetic Field Amplification by Turbulence in A Relativistic Shock Propagating through An Inhomogeneous Medium

We perform two-dimensional relativistic magnetohydrodynamic simulations of a mildly relativistic shock propagating through an inhomogeneous medium. We show that the postshock region becomes turbulent owing to preshock density inhomogeneity, and the magnetic field is strongly amplified due to the stretching and folding of field lines in the turbulent velocity field. The amplified magnetic field evolves into a filamentary structure in two-dimensional simulations. The magnetic energy spectrum is flatter than the Kolmogorov spectrum and indicates that a so-called small-scale dynamo is occurring in the postshock region. We also find that the amount of magnetic-field amplification depends on the direction of the mean preshock magnetic field, and the time scale of magnetic-field growth depends on the shock strength.

preprint2010arXiv

Three-Dimensional Relativistic Magnetohydrodynamic Simulations of Current-Driven Instability. II. Relaxation of Pulsar Wind Nebula

We have investigated the relaxation of a hydrostatic hot plasma column containing toroidal magnetic field by the Current-Driven (CD) kink instability as a model of pulsar wind nebulae. In our simulations the CD kink instability is excited by a small initial velocity perturbation and develops turbulent structure inside the hot plasma column. We demonstrate that, as envisioned by Begelman, the hoop stress declines and the initial gas pressure excess near the axis decreases. The magnetization parameter σ, the ratio of the Poynting to the kinetic energy flux, declines from an initial value of 0.3 to about 0.01 when the CD kink instability saturates. Our simulations demonstrate that axisymmetric models strongly overestimate the elongation of the pulsar wind nebulae. Therefore, the previous requirement for an extremely low pulsar wind magnetization can be abandoned. The observed structure of the pulsar wind nebulae do not contradict the natural assumption that the magnetic energy flux still remains a good fraction of the total energy flux after dissipation of alternating fields.