Transport properties of microstructured ultrathin films of La0.67Ca0.33MnO3 on SrTiO3

We have investigated the electrical transport properties of 8 nm thick La0.67Ca0.33MnO3 films, sputter-deposited on SrTiO3 (STO), and etched into 5 micrometer-wide bridges by Ar-ion etching. We find that even slight overetching of the film leads to conductance of the STO substrate, and asymmetric and non-linear current-voltage (I-V) characteristics. However, a brief oxygen plasma etch allows full recovery of the insulating character of the substrate. The I-V characteristics of the bridges are then fully linear over a large range of current densities. We find colossal magnetoresistance properties typical for strained LCMO on STO but no signature of non-linear effects (so-called electroresistance) connected to electronic inhomogeneites. In the metallic state below 150 K, the highest current densities lead to heating effects and non-linear I-V characteristics.Comment: 3 pages, 5 figure

Tearing instability in relativistic magnetically dominated plasmas

Many astrophysical sources of high energy emission, such as black hole magnetospheres, superstrongly magnetized neutron stars (magnetars), and probably relativistic jets in Active Galactic Nuclei and Gamma Ray Bursts involve relativistically magnetically dominated plasma. In such plasma the energy density of magnetic field greatly exceeds the thermal and the rest mass energy density of particles. Therefore the magnetic field is the main reservoir of energy and its dissipation may power the bursting emission from these sources, in close analogy to Solar flares. One of the principal dissipative instabilities that may lead to release of magnetic energy is the tearing instability. In this paper we study, both analytically and numerically, the development of tearing instability in relativistically magnetically-dominated plasma using the framework of resistive magnetodynamics. We confirm and elucidate the previously obtained result on the growth rate of the tearing mode: the shortest growth time is the same as in the case of classical non-relativistic MHD, namely $\tau =\sqrt{\tau_a \tau_d}$ where $\tau_a$ is the \Alfven crossing time and $\tau_d$ is the resistive time of a current layer.Comment: Submitted to MNRAS, few typos correcte

On the origin of variable gamma-ray emission from the Crab Nebula

The oblique geometry of pulsar wind termination shock ensures that the Doppler beaming has a strong impact on the shock emission. We illustrate this using recent relativistic MHD simulations of the Crab Nebula and also show that the observed size, shape, and distance from the pulsar of the Crab Nebula inner knot are consistent with its interpretation as a Doppler-boosted emission from the termination shock. If the electrons responsible for the synchrotron gamma-rays are accelerated only at the termination shock then their short life-time ensures that these gamma-rays originate close to the shock and are also strongly effected by the Doppler beaming. As the result, bulk of the observed synchrotron gamma-rays of the Crab Nebula around 100 MeV may come from its inner knot. This hypothesis is consistent with the observed optical flux of the inner knot provided its optical-gamma spectral index is the same as the injection spectral index found in the Kennel & Coroniti model of the nebula spectrum. The observed variability of synchrotron gamma-ray emission can be caused by the instability of the termination shock discovered in recent numerical simulations. Given the small size of the knot, it is possible that the September 2010 gamma-ray flare of the Crab Nebula also came from the knot, though the actual mechanism remains unclear. The model predicts correlation of the temporal variability of the synchrotron gamma-ray flux in the Fermi and AGILE windows with the variability of the unpulsed optical flux from within 1 arcsec of the Crab pulsar.Comment: submitted to MNRAS, typos corrected, new references added, additional issues discusse

Do the mildly superluminal VLBI knots exclude ultrarelativistic blazar jets?

We compute the effective values of apparent transverse velocity and flux boosting factors for the VLBI radio knots of blazar jets, by integrating over the angular distributions of these quantities across the widths of jets with finite opening angles but constant velocities. For high bulk Lorentz factors (Gamma > 10) variations across the jet can be quite large if the opening angle, omega, is even a few degrees on sub-parsec scales. The resulting apparent speeds are often much lower than those obtained from the usual analyses that ignore the finite jet opening angles. We can thus reconcile the usually observed subluminal or mildly superluminal speeds with the very high (>~ 20) Gamma factors, required by the inverse Compton origin and rapid variability of TeV fluxes, as well as by intraday radio variability. Thus it is possible to associate the VLBI radio knots directly with shocks in the ultra-relativistic main jet flow, without invoking very rapid jet deceleration on parsec scales, or extremely unlikely viewing angles.Comment: 10 pages, 1 figure, to appear in ApJ Letters, Nov. 10 2004 issu

Relativistic MHD Simulations of Jets with Toroidal Magnetic Fields

This paper presents an application of the recent relativistic HLLC approximate Riemann solver by Mignone & Bodo to magnetized flows with vanishing normal component of the magnetic field. The numerical scheme is validated in two dimensions by investigating the propagation of axisymmetric jets with toroidal magnetic fields. The selected jet models show that the HLLC solver yields sharper resolution of contact and shear waves and better convergence properties over the traditional HLL approach.Comment: 12 pages, 5 figure

A multi-dimensional numerical scheme for two-fluid Relativistic MHD

The paper describes an explicit multi-dimensional numerical scheme for Special Relativistic Two-Fluid Magnetohydrodynamics of electron-positron plasma and a suit of test problems. The scheme utilizes Cartesian grid and the third order WENO interpolation. The time integration is carried out using the third order TVD method of Runge-Kutta type, thus ensuring overall third order accuracy on smooth solutions. The magnetic field is kept near divergence-free by means of the method of generalized Lagrange multiplier. The test simulations, which include linear and non-linear continuous plasma waves, shock waves, strong explosions and the tearing instability, show that the scheme is sufficiently robust and confirm its accuracy

3C454.3 reveals the structure and physics of its 'blazar zone'

Recent multi-wavelength observations of 3C454.3, in particular during its giant outburst in 2005, put severe constraints on the location of the 'blazar zone', its dissipative nature, and high energy radiation mechanisms. As the optical, X-ray, and millimeter light-curves indicate, significant fraction of the jet energy must be released in the vicinity of the millimeter-photosphere, i.e. at distances where, due to the lateral expansion, the jet becomes transparent at millimeter wavelengths. We conclude that this region is located at ~10 parsecs, the distance coinciding with the location of the hot dust region. This location is consistent with the high amplitude variations observed on ~10 day time scale, provided the Lorentz factor of a jet is ~20. We argue that dissipation is driven by reconfinement shock and demonstrate that X-rays and gamma-rays are likely to be produced via inverse Compton scattering of near/mid IR photons emitted by the hot dust. We also infer that the largest gamma-to-synchrotron luminosity ratio ever recorded in this object - having taken place during its lowest luminosity states - can be simply due to weaker magnetic fields carried by a less powerful jet.Comment: 19 pages, 3 figures, accepted for publication in Ap

Synchrotron Self-Compton Model for Rapid Nonthermal Flares in Blazars with Frequency-Dependent Time Lags

We model rapid variability of multifrequency emission from blazars occurring across the electromagnetic spectrum (from radio to gamma-rays). Lower energy emission is produced by the synchrotron mechanism, whereas higher energy emission is due to inverse Compton scattering of the synchrotron emission. We take into account energy stratification established by particle acceleration at shock fronts and energy losses due to synchrotron emission. We also consider the effect of light travel delays for the synchrotron emission that supplies the seed photons for inverse Compton scattering. The production of a flare is caused by the collision between a relativistic shock wave and a stationary feature in the jet (e.g., a Mach disk). The collision leads to the formation of forward and reverse shocks, which confine two contiguous emission regions resulting in complex profiles of simulated flares. Simulations of multifrequency flares indicate that relative delays between the inverse Compton flares and their synchrotron counterparts are dominated by energy stratification and geometry of the emitting regions, resulting in both negative and positive time delays depending on the frequency of observation. Light travel effects of the seed photons may lead to a noticeable delay of the inverse Compton emission with respect to synchrotron variability if the line of sight is almost perfectly aligned with the jet. We apply the model to a flare in 3C 273 and derive the properties of shocked plasma responsible for the flare. We show that the pronounced negative time delay between the X-ray and IR light curves (X-rays peak after the maximum in the synchrotron emission) can be accounted for if both forward and reverse shocks are considered.Comment: 48 pages, 18 figures, accepted for publication in The Astrophysical Journa

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