A HED Laboratory Astrophysics Testbed Comes of Age: JET Deflection via Cross Winds

We present new data from High-Energy Density (HED) laboratory experiments designed to explore the interaction of a heavy hypersonic radiative jet with a cross wind. The jets are generated with the MAGPIE pulsed power machine where converging conical plasma flows are produced from a cylindrically symmetric array of inclined wires. Radiative hypersonic jets emerge from the convergence point. The cross wind is generated by ablation of a plastic foil via soft-X-rays from the plasma convergence region. Our experiments show that the jets are deflected by the action of the cross wind with the angle of deflection dependent on the proximity of the foil. Shocks within the jet beam are apparent in the data. Analysis of the data shows that the interaction of the jet and cross wind is collisional and therefore in the hydrodynamic regime. We consider the astrophysical relevance of these experiments applying published models of jet deflection developed for AGN and YSOs. We also present results of 3-D numerical simulations of jet deflection using a new astrophysical Adaptive Mesh Refinement code. These simulations show highly structured shocks occurring within the beam similar to what was observed in the experiments.     

Formation of ionization-cone structures in active galactic nuclei: II. Nonlinear hydrodynamic modelling

In Part I of this paper (Paper I) we described an equilibrium model of a jet in the gravitational field corresponding to the rigid-rotation region of the galactic disk. We used linear stability analysis to find the waveguide-resonance instability of internal gravity waves due to the superreflection of these waves from the jet boundary. In this part of the paper, we perform nonlinear numerical 2D and 3D simulations of the development of this instability. We show that the shocks produced by this instability in the ambient medium of the jet are localized inside a cone with a large opening angle and are capable of producing features that are morphologically similar to those observed in galaxies with active nuclei (NGC 5252 for example).     

Deflection of jets induced by jet–cloud and jet–galaxy interactions

The model first introduced by Raga & Cantó in which astrophysical jets are deflected on passing through an isothermal high‐density region is generalized by taking into account gravitational effects on the motion of the jet as it crosses the high‐density cloud. The problem is also generalized for relativistic jets in which gravitational effects induced by the cloud are neglected. Two further cases, classical and relativistic, are discussed for the cases in which the jet is deflected on passing through the interstellar gas of a galaxy in which a dark matter halo dominates the gravitational potential. The criteria for the stability of jets due to the formation of internal shocks are also discussed.     

iShocks: X-ray binary jets with an internal shocks model

In the following paper, we present an internal shocks model, iShocks, for simulating a variety of relativistic jet scenarios; these scenarios can range from a single ejection event to an almost continuous jet, and are highly user configurable. Although the primary focus in the following paper is black hole X-ray binary jets, the model is scale and source independent and could be used for supermassive black holes in active galactic nuclei or other flows such as jets from neutron stars. Discrete packets of plasma (or 'shells') are used to simulate the jet volume. A two-shell collision gives rise to an internal shock, which acts as an electron re-energization mechanism. Using a pseudo-random distribution of the shell properties, the results show how for the first time it is possible to reproduce a flat/inverted spectrum (associated with compact radio jets) in a conical jet whilst taking the adiabatic energy losses into account. Previous models have shown that electron re-acceleration is essential in order to obtain a flat spectrum from an adiabatic conical jet: multiple internal shocks prove to be efficient in providing this re-energization. We also show how the high-frequency turnover/break in the spectrum is correlated with the jet power,  ν _b ∝ L ^∼0.6_W  , and the flat-spectrum synchrotron flux is correlated with the total jet power,   F _ν∝ L ^∼1.4_W  . Both the correlations are in agreement with previous analytical predictions.     

Multiple relativistic outbursts of GRS 1915+105: radio emission and internal shocks

We present 5-GHz Multi-Element Radio-Linked Interferometer Network (MERLIN) radio images of the microquasar GRS 1915+105 during two separate outbursts in 2001 March and July, following the evolution of the jet components as they move outwards from the core of the system. Proper motions constrain the intrinsic jet speed to be  >0.57 c   , but the uncertainty in the source distance prevents an accurate determination of the jet speed. No deceleration is observed in the jet components out to an angular separation of ∼300 mas. Linear polarization is observed in the approaching jet component, with a gradual rotation in position angle and a decreasing fractional polarization with time. Our data lend support to the internal shock model whereby the jet velocity increases leading to internal shocks in the pre-existing outflow before the jet switches off. The compact nuclear jet is seen to reestablish itself within 2 d, and is visible as core emission at all epochs. The energetics of the source are calculated for the possible range of distances; a minimum power of 1–10 per cent of the Eddington luminosity ( L _Edd) is required to launch the jet.     

赫比格-哈罗天体高分辨观测研究进展

赫比格－哈罗天体（ＨＨ天体）包含了有关原恒星吸积和抛射过程的许多重要信息,ＨＨ天体高分辨观测研究取得了一系列新进展：分辨出激波峰面、马赫盘和辐射冷却区;分辨出喷流节点的结构,发现它们大多是内工作面,而不是由Ｋｅｌｖｉｎ－Ｈｅｌｍｈｏｌｔｚ不稳定性所产生的斜激波;发现喷流宽度随到激发源距离的减小仅缓慢减小,对喷流的准直和加速模型提供了限制条件;ＨＨ天体在小尺度上尚有复杂的激发结构。对这些进展进行了评     

Experiments and Numerical Simulations on the Mid-Term Evolution of Hypersonic Jets

The experiment described here is focussed to the observation of underexpanded, hypersonic turbulent jets. The experiment is relevant to a few aspects concerning the dynamics of astrophysical phenomena such as the Herbig-Haro jets and to the interaction between the large-scale vortices and the system of shocks that determine the spreading and entrainment properties of highly compressible free-flows. A number of orifice jets with a ratio between the stagnation pressure and the ambient pressure of the order 10³-10⁴ have been studied by changing the stagnation/ambient density ratio. This has been realized using dissimilar gases in the jet and in the ambient medium: by using He, Ar and air the stagnation/ambient density ratio can be changed by one order of magnitude while keeping fixed the pressure ratio. It has been possible to visualize the near and mid-term evolution of the jets and measure the axial and transversal density distributions. A comparison relevant to the shock waves configuration in between the nozzle exit and the first Mach's disk is shown for an air in air laboratory jet and its numerical simulation.     

Modelling jet-driven molecular outflows

We have investigated the basic physical properties of the outflow that is created by a supersonic jet in a dense molecular cloud. We show that the dynamics of the interaction is strongly controlled by the rapid cooling of the post-shock gas at the head of the jet. The velocity of the gas is high in the vicinity of the jet head, but decreases rapidly as more material is swept-up. This type of outflow produces extremely high velocity clumps of post-shock gas which resemble the features seen in outflows. We also show that momentum transfer in bow shocks is more important than entrainment in high Mach number jets, as found in the protostellar environment.     

Recent Experimental Results and Modelling of High-Mach-Number Jets and the Transition to Turbulence

In recent years, we have carried out experiments at the University of Rochester’s Omega laser in which supersonic, dense-plasma jets are formed by the interaction of strong shocks in a complex target assembly (Foster et al., Phys. Plasmas 9 (2002) 2251). We describe recent, significant extensions to this work, in which we consider scaling of the experiment, the transition to turbulence, and astrophysical analogues. In new work at the Omega laser, we are developing an experiment in which a jet is formed by laser ablation of a titanium foil mounted over a titanium washer with a central, cylindrical hole. Some of the resulting shocked titanium expands, cools, and accelerates through the vacuum region (the hole in the washer) and then enters a cylinder of low-density foam as a jet. We discuss the design of this new experiment and present preliminary experimental data and results of simulations using AWE hydrocodes. In each case, the high Reynolds number of the jet suggests that turbulence should develop, although this behaviour cannot be reliably modelled by present, resolution-limited simulations (because of their low-numerical Reynolds number).     

Simulations of the supersonic radiative jet propagation in plasmas

Supersonic plasma jets are ubiquitous in astrophysics. Our study focus on the jets emanated from Herbig-Haro (HH) objects. They have velocities of a few hundred km/s and are extending over the distances more than a parsec. Interaction of the jets with surrounding matter produces two specific structures in the jet head: the bow shock and the Mach disk. The radiative cooling of these shocks affects strongly the jet dynamics. A tool to understand the physics of these jets is the laboratory experiment. A supersonic jet interaction with surrounding plasma was studied on the PALS laser facility. A collimated high-Z plasma jet with a velocity exceeding 400 km/s was generated and propagated over a few millimeters length. Here we report on study the effect of radiative cooling on the head jet structure with a 2D radiative hydrodynamic code. The simulation results demonstrated the scalability of the experimental observations to the HH jets.     

Radio Band Observations of Blazar Variability

The properties of blazar variability in the radio band are studied using the unique combination of temporal resolution from single dish monitoring and spatial resolution from VLBA imaging. Such measurements now available in all four Stokes parameters, together with theoretical simulations, identify the origin of radio band variability and probe the characteristics of the radio jet where the broadband blazar emission originates. Outbursts in total flux density and linear polarization in the optical-to-radio bands are attributed to shocks propagating within the jet spine, in part, based on limited modelling invoking transverse shocks; new radiative transfer simulations allowing for shocks at arbitrary angle to the flow direction confirm this picture by reproducing the observed centimeter-band variations observed more generally, and are of current interest since these shocks may play a role in the ??-ray flaring detected by Fermi. Recent UMRAO multifrequency Stokes V studies of bright blazars identify the spectral variability properties of circular polarization for the first time and demonstrate that polarity flips are relatively common. All-Stokes data are consistent with the production of circular polarization by linear-to-circular mode conversion in a region that is at least partially self-absorbed. Detailed analysis of single-epoch, multifrequency, all-Stokes VLBA observations of 3C?279 support this physical picture and are best explained by emission from an electron-proton plasma.     

VSOP polarization observations of the BL Lacertae object OJ 287

VLBI total intensity and linear polarization images of the BL Lacertae object OJ 287 have been obtained at     using a global ground array and the HALCA orbiting antenna, and at     two weeks earlier using the VLBA. In the ground-based 6-cm images, the source is dominated by a core–jet double structure the components of which are essentially unresolved. The baselines to the orbiting antenna resolve both of these compact components. In the VSOP images, the ground-based 'core' breaks up into several distinct components, demonstrating that this region is dominated by the contribution of bright, optically thin knots of jet emission. A very similar structure is observed in the 1.3-cm image. The magnetic field in the core is transverse, becomes longitudinal in the inner jet, then makes a sharp transition to a region of transverse field further from the core. This suggests that the field in the outer jet has become highly ordered in the transverse direction owing to the action of a shock; the physical nature of the extended region of longitudinal field closer to the core is not clear. The availability of nearly simultaneous observations with comparable resolution at widely spaced frequencies enabled detection of a ≃90° rotation in polarization position angle for the core, owing to the transition from the optically thick (6 cm) to the optically thin (1.3 cm) regime.     

3-D Rpic Simulations of Relativistic Jets: Particle Acceleration, Magnetic Field Generation, and Emission

We have applied numerical simulations and modeling to the particle acceleration, magnetic field generation, and emission from relativistic shocks. We investigate the nonlinear stage of theWeibel instability and compare our simulations with the observed gamma-ray burst emission. In collisionless shocks, plasma waves and their associated instabilities (e.g., the Weibel, Buneman and other two-stream instabilities) are responsible for particle (electron, positron, and ion) acceleration and magnetic field generation. 3-D relativistic electromagnetic particle (REMP) simulations with three different electron-positron jet velocity distributions and also with an electron-ion plasma have been performed and show shock processes including spatial and temporal evolution of shocks in unmagnetized ambient plasmas. The growth time and nonlinear saturation levels depend on the initial jet parallel velocity distributions. Simulations show that the Weibel instability created in the collisionless shocks accelerates jet and ambient particles both perpendicular and parallel to the jet propagation direction. The nonlinear fluctuation amplitude of densities, currents, electric, and magnetic fields in the electron-positron shocks are larger for smaller jet Lorentz factor. This comes from the fact that the growth time of the Weibel instability is proportional to the square of the jet Lorentz factor. We have performed simulations with broad Lorentz factor distribution of jet electrons and positrons, which is assumed to be created by photon annihilation. Simulation results with this broad distribution show that the Weibel instability is excited continuously by the wide-range of jet Lorentz factor from lower to higher values. In all simulations the Weibel instability is responsible for generating and amplifying magnetic fields perpendicular to the jet propagation direction, and contributes to the electron’s (positron’s) transverse deflection behind the jet head. This small scale magnetic field structure contributes to the generation of “jitter” radiation from deflected electrons (positrons), which is different from synchrotron radiation in uniform magnetic fields. The jitter radiation resulting from small scale magnetic field structures may be important for understanding the complex time structure and spectral evolution observed in gamma-ray bursts or other astrophysical sources containing relativistic jets and relativistic collisionless shocks. The detailed studies of shock microscopic process evolution may provide some insights into early and later GRB afterglows.     

The Bastille day Shocks and Merged Interaction Region

Near 1 AU the solar wind structure associated with the solar flare of 14 July 2000 (Bastille Day) consisted of a large high-speed stream of 15 July and five nearby small streams during a 10-day period. At the leading edge of the large high-speed stream, in less than 6 hours, the flow speed increased from 600 km s⁻¹ to 1100 km s⁻¹, the magnetic field intensity increased from 10 nT to 60 nT, and an interaction region was identified. The interaction region was bounded between the pair of a forward shock F and a reverse shock R. Additional forward shocks were also identified at the leading edge of each of the five smaller streams. This paper presents a magnetohydrodynamics (MHD) simulation using ACE plasma and magnetic field data near 1 AU as input to study the radial evolution of the Bastille Day solar wind event. The two shocks, F and R, propagated in opposite directions away from each other in the solar wind frame and interacted with neighboring shocks and streams; the spatial and temporal extent of the interaction region continued to increase with the heliocentric distance. The solar wind was restructured from a series of streams at 1 AU to a huge merged interaction region (MIR) extending over a period of 12 days at 5.5 AU. Throughout the interior of the MIR bounded by the shock pair F and R the magnetic field intensity was a few times stronger than that outside the MIR. The simulation shows how merging of shocks, collision of shocks, and formation of new shocks contributed to the evolution process.     

On the dynamic efficiency of internal shocks in magnetized relativistic outflows

We study the dynamic efficiency of conversion of kinetic-to-thermal/magnetic energy of internal shocks in relativistic magnetized outflows. We model internal shocks as being caused by collisions of shells of plasma with the same energy flux and a non-zero relative velocity. The contact surface, where the interaction between the shells takes place, can break up either into two oppositely moving shocks (in the frame where the contact surface is at rest), or into a reverse shock and a forward rarefaction. We find that for moderately magnetized shocks (magnetization  σ≃ 0.1  ), the dynamic efficiency in a single two-shell interaction can be as large as 40 per cent. Thus, the dynamic efficiency of moderately magnetized shocks is larger than in the corresponding unmagnetized two-shell interaction. If the slower shell propagates with a sufficiently large velocity, the efficiency is only weakly dependent on its Lorentz factor. Consequently, the dynamic efficiency of shell interactions in the magnetized flow of blazars and gamma-ray bursts is effectively the same. These results are quantitatively rather independent on the equation of state of the plasma. The radiative efficiency of the process is expected to be a fraction   f _r < 1  of the estimated dynamic one, the exact value of f _r depending on the particularities of the emission processes which radiate away the thermal or magnetic energy of the shocked states.     

Numerical Simulations of the Interaction of Jets with the Intracluster Medium

We study, by numerical simulations, the propagation of an axisymmetric supersonic jet in an isothermal King atmosphere and we analyse the evolution of the resulting X-ray properties and their dependence on the jet physical parameters. We show the existence of two distinct regimes of interaction, with strong and weak shocks. In the first case shells of enhanced X-ray emission are to be expected, whereas in the second case we expect deficit of X-ray emission coincident with the cocoon. Analysing the results of our simulations we find that the jet kinetic power is the main parameter controlling the transition between the two regimes. We also discuss, in the same scheme, the ICM heating induced by the jet propagation, considering its effects on the observed relations between the cluster X-ray luminosity and temperature and between cluster entropy and temperature.     

Distinguishing Propagation vs. Launch Physics of Astrophysical Jets and the Role of Experiments

The absence of other viable momentum sources for collimated flows leads to the likelihood that magnetic fields play a fundamental role in jet launch and/or collimation in astrophysical jets. To best understand the physics of jets, it is useful to distinguish between the launch region where the jet is accelerated and the larger scales where the jet propagates as a collimated structure. Observations presently resolve jet propagation, but not the launch region. Simulations typically probe the launch and propagation regions separately, but not both together. Here, I IDentify some of the physics of jet launch vs. propagation and what laboratory jet experiments to date have probed. Reproducing an astrophysical jet in the lab is unrealistic, so maximizing the benefit of the experiments requires clarifying the astrophysical connection.     

Stability of Relativistic Hydrodynamical Planar Jets: Linear and Nonlinear Evolution of Kelvin-Helmholtz Modes

Some aspects about the stability of relativistic flows against Kelvin-Helmholtz (KH) perturbations are studied by means of relativistic, hydrodynamical simulations. In particular, we analyze the transition to the fully nonlinear regime and the long-term evolution of two jet models with different specific internal energies.     

Critical Self-Organization of Astrophysical Shocks

There are two distinct regimes of the first-order Fermi acceleration of shocks. The first is a linear (test-particle) regime in which most of the shock energy goes into thermal and bulk motions of the plasma. The second is an efficient regime in which the shock energy goes into accelerated particles. Although the transition region between them is narrow, we identify the factors that drive the system toward a self-organized critical state between those two regimes. Using an analytic solution, we determine this critical state and calculate the spectra and maximum energy of accelerated particles.     

Laboratory Simulations of Supernova Shockwave Propagation

Supernovae launch spherical shocks into the circumstellar medium (CSM). These shocks have high Mach numbers and may be radiative. We have created similar shocks in the laboratory by focusing laser pulses onto the tip of a solid pin surrounded by ambient gas; ablated material from the pin rapidly expands and launches a shock through the surrounding gas. Laser pulses were typically 5 ns in duration with ablative energies ranging from 1–150 J. Shocks in ambient gas pressures of ~1 kPa were observed at spatial scales of up to 5 cm using optical cameras with schlieren. Emission spectroscopy data were obtained to infer electron temperatures (< 10 eV). In this experiment we have observed a new phenomena; at the edge of the radiatively heated gas ahead of the shock, a second shock forms. The two expanding shocks are simultaneously visible for a time, until the original shock stalls from running into the heated gas. The second shock remains visible and continues to expand. A minimum condition for the formation of the second shock is that the original shock is super-critical, i.e., the temperature distribution ahead of the original shock has an inflexion point. In a non-radiative control experiment the second shock does not form. We hypothesize that a second shock could form in the astrophysical case, possibly in radiative supernova remnants such as SN1993J, or in shock-CSM interaction.