We present results of two-and three-dimensional Particle-In-Cell simulations of magnetic-turbulence production by isotropic cosmic-ray ions drifting upstream of supernova remnant shocks. The studies aim at testing recent predictions of a strong amplification of short-wavelength magnetic field and at studying the subsequent evolution of the magnetic turbulence and its backreaction on cosmic ray trajectories. We observe that an oblique filamentary mode grows more rapidly than the non-resonant parallel modes found in analytical theory, and the growth rate of the field perturbations is much slower than is estimated for the parallel plane-wave mode, possibly because in our simulations we cannot maintain ω ≪ Ω i , the ion gyrofrequency, to the degree required for the plane-wave mode to emerge. The evolved oblique filamentary mode was also observed in MHD simulations to dominate in the non-linear phase, when the structures are already isotropic. We thus confirm the generation of the turbulent magnetic field due to the drift of cosmic-ray ions in the upstream plasma, but as our main result find that the amplitude of the turbulence saturates at about δB/B ∼ 1. The backreaction of the magnetic turbulence on the particles leads to an alignment of the bulk-flow velocities of the cosmic rays and the background medium. This is an essential characteristic of cosmic-ray modified shocks: the upstream flow speed is continuously changed by the cosmic rays. The deceleration of the cosmic-ray drift

We have investigated the development of current-driven (CD) kink instability through three-dimensional relativistic MHD simulations. A static force-free equilibrium helical magnetic configuration is considered in order to study the influence of the initial configuration on the linear and nonlinear evolution of the instability. We found that the initial configuration is strongly distorted but not disrupted by the kink instability. The instability develops as predicted by linear theory. In the non-linear regime the kink amplitude continues to increase up to the terminal simulation time, albeit at different rates, for all but one simulation. The growth rate and nonlinear evolution of the CD kink instability depends moderately on the density profile and strongly on the magnetic pitch profile. The growth rate of the kink mode is reduced in the linear regime by an increase in the magnetic pitch with radius and the non-linear regime is reached at a later time than for constant helical pitch. On the other hand, the growth rate of the kink mode is increased in the linear regime by a decrease in the magnetic pitch with radius and reaches the non-linear regime sooner than the case with constant magnetic pitch. Kink amplitude growth in the non-linear regime for decreasing magnetic pitch leads to a slender helically twisted column wrapped by magnetic field. On the other hand, kink amplitude growth in the non-linear regime nearly ceases for increasing magnetic pitch.

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.

Shock acceleration is a ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., Buneman, Weibel, and other two-stream instabilities) created in collisionless shocks are responsible for particle (electron, positron, and ion) acceleration. Using a three-dimensional relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic electron-positron jet front propagating into an ambient electron-positron plasma with and without initial magnetic fields. We find small differences in the results for no ambient and modest ambient magnetic fields. New simulations show that the Weibel instability created in the collisionless shock front accelerates jet and ambient particles both perpendicular and parallel to the jet propagation direction. Furthermore, the nonlinear fluctuation amplitudes of densities, currents, and electric and magnetic fields in the electron-positron shock are larger than those found in the electron-ion shock studied in a previous paper at a comparable simulation time. This comes from the fact that both electrons and positrons contribute to generation of the Weibel instability. In addition, we have performed simulations with different electron skin depths. We find that growth times scale inversely with the plasma frequency, and the sizes of structures created by the Weibel instability scale proportionally to the electron skin depth. This is the expected result and indicates that the simulations have sufficient grid resolution. While some Fermi acceleration may occur at the jet front, the majority of electron and positron acceleration takes place behind the jet front and cannot be characterized as Fermi acceleration. The simulation results show that the Weibel instability is responsible for generating and amplifying nonuniform, small-scale magnetic fields, which contribute to the electron's (positron's) transverse deflection behind the jet head. This smallscale magnetic field structure is appropriate to the generation of ''jitter'' radiation from deflected electrons (positrons) as opposed to synchrotron radiation. The jitter radiation has different properties than synchrotron radiation calculated assuming a uniform magnetic field. 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.

Numerical simulations of weakly magnetized and strongly magnetized relativistic jets embedded in a weakly magnetized and strongly magnetized stationary or weakly relativistic (v = c/2) sheath have been performed. A magnetic field parallel to the flow is used in these simulations performed by the new GRMHD numerical code RAISHIN used in its RMHD configuration. In the numerical simulations the Lorentz factor γ = 2.5 jet is precessed to break the initial equilibrium configuration. In the simulations sound speeds are c/ √ 3 in the weakly magnetized simulations and 0.3 c in the strongly magnetized simulations. The Alfvén wave speed is 0.07 c in the weakly magnetized simulations and 0.56 c in the strongly magnetized simulations. The results of the numerical simulations are compared to theoretical predictions from a normal mode analysis of the linearized relativistic magnetohydrodynamic (RMHD) equations capable of describing a uniform axially magnetized cylindrical relativistic jet embedded in a uniform axially magnetized relativistically moving sheath. The theoretical dispersion relation allows investigation of effects associated with maximum possible sound speeds, Alfvén wave speeds near light speed and relativistic sheath speeds. The prediction of increased stability of the weakly magnetized system resulting from c/2 sheath speeds and the stabilization of the strongly magnetized system resulting from c/2 sheath speeds is verified by the numerical simulation results.

We have performed two-dimensional special-relativistic magnetohydrodynamic simulations of nonequilibrium 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.

The enhancement effect of nucleation in immiscible blend systems has recently attracted interest. Although several authors have reported that the effect occurs at the phase interface, little is known about the mechanism involved. We focused on poly(L-lactide) (PLLA)/poly(ɛ-caprolactone) (PCL) immiscible blend systems in which the presence of PCL enhanced the nucleation of PLLA at low temperature. We investigated the nucleation behavior of PLLA during aging at temperatures below T g . Generally, neat polymers, including PLLA, seldom generate nuclei below T g due to restrictions in chain mobility. However, through DSC analysis of the crystallization behavior following an aging process, we revealed that the nucleation of PLLA occurs during aging even at temperatures below T g in the PLLA/PCL blend. Since the nuclei density became saturated with increasing aging time, the nucleation behavior was regarded as heterogeneous nucleation. The asymptotic density of nuclei depended on the PCL content, indicating that dispersed PCL acted as active sites for nucleation. The nucleation rate R was almost independent of the aging temperature, suggesting that the marked decrease in chain mobility due to the glass transition is locally evaded at the active sites. Nucleation was observed even at temperatures as much as 40 °C lower than T g following the addition of only 1 wt % PCL, while the T g obtained by a DSC heating scan showed a subtle decrease. This suggests that the limited miscibility of PLLA/PCL leads to the aggregation of PCL and induces local and deep depression of T g at the interface of the PCL domains, resulting in marked enhancement of PLLA nucleation.

No abstract