MHD interaction of pulsar wind nebulae with SNRs and with the ISM

In the late 1960s the discovery of the Crab pulsar in its associated supernova remnant, launched a new field in supernova remnant research: the study of pulsar-driven or plerionic supernova remnants. In these type of remnants, the relativistic wind emitted by the pulsar, blows a pulsar wind nebula into the interior of its supernova remnant. Now, more then forty years after the discovery of the Crab pulsar, there are more then fifty plerionic supernova remnants known, due to the ever-increasing capacity of observational facilities. These observational studies reveal a Zoo of complex morphologies over a wide range of frequencies, indicating the significance of the interaction between a pulsar wind nebula with its surrounding supernova remnant. A pulsar which gained a kick velocity at birth, will ultimately break outside of its remnant, after which the pulsar wind nebula interacts directly with the interstellar medium. In general these pulsar wind nebulae are bounded by a bow shock, due to the supersonic motion of the pulsar. There are a few examples known of these pulsar-powered bow shocks, a number which is slowly increasing.I will review our current understanding of the different evolutionary stages of a pulsar wind nebula as it is interacting with its associated supernova remnant. Therefore, I will discuss both analytical and more recent numerical (M)HD models. The four main stages of a pulsar wind nebula are: the supersonic expansion stage, the reverse shock interaction stage, the subsonic expansion stage and ultimately the stage when the head of the bubble is bounded by a bow shock, due to the supersonic motion of the pulsar. Ultimately this pulsar wind nebula bow shock will break through its associated remnant, after which the pulsar-powered bow shock will interact directly with the interstellar medium. I will discuss recent numerical models from these type of pulsar wind nebulae and their morphology.     

TeV gamma-ray observations of pulsar wind nebulae with the HESS detector

High energy stereoscopic system (HESS) is a recent operational detector dedicated to the observation of γ-rays in the very high energy range. Situated in Namibia, it is composed of four Imaging Atmospheric Cherenkov Telescopes and gives a significant improvement in sensitivity and in accuracy of the reconstructed γ-ray parameters. First results on observations of pulsar wind nebulae are reported here. The binary system PSR B1259-63/SS 2883 has been detected in 2004 around the periastron, showing clear flux variations. The pulsar wind nebula around PSR B1706-44 has been observed and upper limits on its flux have been derived.     

The hard X-ray emission from the Vela pulsar wind nebula

We present the results of a preliminary spectral analysis performed on the BeppoSAX and XMM observations of the Vela plerion. The broad energy range covered by the instruments on board the two observatories allows an evaluation of the spectral parameters of the high energy emission model and provides an indication on the morphology of the source emission above 10 keV. We confirm the softening of the PWN spectrum (3–10 keV band) at distances greater than 4′ from the pulsar and estimate the diameter of the high energy (>10 keV) emission region to be on the order of 25′–30′.     

The termination shock in a striped pulsar wind

I present observational and theoretical evidence that most of the pulsar spin-down energy is transferred away as a striped pulsar wind and that this energy is released by annihilation of the alternating magnetic field at the pulsar wind termination shock. One-dimensional particle-in-cells (PIC) simulations show that the alternating fields do annihilate at the termination shock in a striped wind. The particle acceleration should be studied in multidimensional simulations. As a first step, I simulated driven annihilation of alternating fields undergoing compression by an external force. It is shown that that in the course of this process, a particle distribution function is formed, which resembles that observed in plerions.     

A VLBI search for pulsar wind nebulae in supernovae

We examine recent supernovae which have been observed with very-long-baseline interferometry in order to detect or limit the emission from a possible compact remnant of the explosion. Such a remnant could be a neutron star, generating a pulsar wind nebula, or a black hole with an accretion disk and jets. Four supernovae, and also more than a dozen supernovae or their young remnants in M82, have structure sufficiently resolved to allow useful conclusions as to the strength of the emission from such young neutron stars or black holes. We recently discovered a compact component in the center of SN 1986J’s shell with a spectral luminosity at 15 GHz 200 times that of the Crab Nebula. This is most likely the compact remnant of the explosion, the first and only one found in any modern supernova. For other modern supernovae, the upper limits on the radio spectral luminosities of such young compact remnants range from 180 times that of the Crab Nebula for SN 1979C in M100 in the Virgo cluster to 0.001 times that of the Crab Nebula for SN 1987A in the Large Magellanic Cloud.     

Solar wind interaction with comet Halley

In March 6 and 9, 1986 the spacecrafts ‘Vega-1’ and ‘Vega-2’ have flown through the coma of comet Halley and have carried measurements of plasma, energetic particles, magnetic field and plasma waves along its trajectory. A short review of these measurements and its comparison with theoretical models of solar wind interaction with comets are given.

The spacecrafts ‘Vega-1’ and ‘Vega-2’ have studied the solar wind loading by cometary ions, the structure of cometary bow shock and the processes in the inner coma of comet Halley. Exactly in this sequence we discuss the results of measurements and compare them with the theory.     

The devil is in the details: Compact structures in pulsar wind nebulae

The large-scale structure of pulsar wind nebulae (PWNe) tells us a considerable amount about their average magnetic fields, the total particle input from the pulsar winds, and the confining pressure at their outer boundaries. However, the details of the pulsar outflow, the sites of shocks and particle acceleration, the effects of instabilities in the magnetic field, and the interaction between the relativistic wind and the surrounding ejecta are contained in small-scale structures, where we observe jets and toroidal structures, time-varying emission from compact clumps, and filaments in both the inner and outer regions of the nebulae. Here, I review recent observational studies of compact structures in PWNe and present current scenarios (and questions) regarding their origin.     

Simulation of cometary jets in interaction with the solar wind

The influence of cometary jets on the solar wind interaction is studied with a 3D hybrid simulation. Anisotropic outgassing patterns were until recently not considered in cometary simulations, despite strong anisotropies found at observations. Comet 67P Churyumov–Gerasimenko, the target of the ROSETTA mission, was chosen as a case study for a simulation series. The cometary outgassing at 2.7 AU is modeled to originate from a single sun-facing jet with different levels of collimation, from isotropy to extremely thin jets. As no bow shock is present at this distance, solar wind patterns resulting from the anisotropic outgassing become more apparent. We find narrower jets to increase the standoff distance of the plasma interaction structures. Also, the Mach cone is wider and stronger for certain jet profiles. The magnetic field remains unable to propagate through the coma, resulting in strong draping patterns for narrow jets due to the increased standoff distance.     

Solar wind interaction with lunar crustal magnetic anomalies

Using Lunar Prospector data, we review the magnetic field and electron signatures of solar wind interaction with lunar crustal magnetic sources. Magnetic field amplifications, too large to represent direct measurements of crustal fields, appear in the solar wind over strong crustal sources, with the chance of observing these amplifications depending on upstream solar wind parameters. We often observe increases in low-energy (?100 eV) electron energy fluxes simultaneously with large magnetic field amplifications, consistent with an increase in plasma density across a shock surface. We also often observe low frequency wave activity in the magnetic field data (both broadband turbulence and monochromatic waves), often associated with electron energization, sometimes up to keV energies. Electron energization appears to be correlated more closely with wave activity than with magnetic amplifications. Detailed studies of the interaction region will be necessary in order to understand the physics of the Moon–solar wind interaction. At present, the Moon represents the only natural laboratory available to us to study solar wind interaction with small-scale crustal magnetic fields, though simulation results and theoretical work can also help us understand the physical processes at work.     

Real-time global MHD simulation of the solar wind interaction with the earth’s magnetosphere

We have developed a real-time global MHD (magnetohydrodynamics) simulation of the solar wind interaction with the earth’s magnetosphere. By adopting the real-time solar wind parameters and interplanetary magnetic field (IMF) observed routinely by the ACE (Advanced Composition Explorer) spacecraft, responses of the magnetosphere are calculated with MHD code. The simulation is carried out routinely on the super computer system at National Institute of Information and Communications Technology (NICT), Japan. The visualized images of the magnetic field lines around the earth, pressure distribution on the meridian plane, and the conductivity of the polar ionosphere, can be referred to on the web site (http://www2.nict.go.jp/y/y223/simulation/realtime/).The results show that various magnetospheric activities are almost reproduced qualitatively. They also give us information how geomagnetic disturbances develop in the magnetosphere in relation with the ionosphere. From the viewpoint of space weather, the real-time simulation helps us to understand the whole image in the current condition of the magnetosphere. To evaluate the simulation results, we compare the AE indices derived from the simulation and observations. The simulation and observation agree well for quiet days and isolated substorm cases in general.     

Geminga pulsar observations with gamma-telescope Gamma-1

The Geminga light curve obtained with the “Gamma-1” telescope features two peaks separated by 0.5 ± 0.03 period. The light curve is pronounced for γ-quanta energies higher than 400 MeV. The pulsed flux upper limit (1σ) in the energy interval 50 – 300 MeV is 6·10⁻⁷ cm⁻²sec⁻¹. For energies >300 MeV the pulsed component power law spectrum has an exponent 1.1 _−0.3^+1.1 and an integral flux (1.1±0.3)·10⁻⁶ cm⁻²sec⁻¹.     

Chandra observation of the interaction between the plasma nebula RCW89 and the pulsar jet of PSR B1509–58

We present a Chandra observation of the H II region RCW89. The nebula lies 10′ north from the central pulsar PSR B1509–58, and it has been suggested that the nebula is irradiated by the pulsar jet. We performed a spectral analysis of the seven brightest emitting regions aligned in a horse-shoe like shape, and found that the temperature of the knots increases along the horse-shoe in the clockwise direction, while, in contrast, the ionization parameter n_et decreases. This result implies that RCW89 was heated in sequence. We examined the energy budget assuming that RCW89 is powered by the pulsar jet. The rate of energy injection into RCW89 by the jet was estimated from the synchrotron radiation flux. We obtained a heating time-scale of 1400 yr, which is consistent with the pulsar characteristic age of 1700 yr. To explain the temperature gradient, we discuss the cooling process for plasma clouds in RCW89. We argue that the plasma clumps can be cooled down by the adiabatic expansion within 250 yr, and form the temperature gradient reflecting the sequential heating by the precessing pulsar jet.     

Challenges to cosmic ray modeling: From beyond the solar wind termination shock

Voyager 1 crossed the solar wind termination shock on December 16, 2004 at a distance of 94 AU from the Sun, to become the first spacecraft to explore the termination shock region and to enter the heliosheath, the final heliospheric frontier. By the end of 2006, Voyager 1 will be at ∼101 AU, with Voyager 2 at ∼81 AU and still approaching the termination shock. Both spacecraft have been observing the modulation of galactic and anomalous cosmic rays since their launch in 1977. The recent observations close to or inside the heliosheath have provided several interesting ‘surprises’ with subsequent theoretical and modeling challenges. Examples are: what does the modulation of galactic cosmic rays amount to in this region?; how do the anomalous cosmic rays get accelerated and modulated?; why are there ‘breaks’ in the power-law slopes of the spectra of accelerated particles? Several numerical models have been applied to most of these topics over the years and comprehensive global predictions have been made the past decade, thought to be based on reasonable assumptions about the termination shock and the heliosheath. Examples of these predictions and assumptions are concisely discussed within the context of the main observed features of cosmic rays in the vicinity of the termination shock, ending with a discussion of some of the issues and challenges to cosmic ray modeling in particular.     

Giant radio pulses from the Crab pulsar

Individual giant radio pulses (GRPs) from the Crab pulsar last only a few microseconds. However, during that time they rank among the brightest objects in the radio sky reaching peak flux densities of up to 1500 Jy even at high radio frequencies. Our observations show that GRPs can be found in all phases of ordinary radio emission including the two high frequency components (HFCs) visible only between 5 and 9 GHz [Moffett, D.A., Hankins, T.H. Multifrequency radio observations of the Crab pulsar. Astrophys. J. 468, 779–783, 1996]. This leads us to believe that there is no difference in the emission mechanism of the main pulse (MP), inter pulse (IP) and HFCs. High resolution dynamic spectra from our recent observations of giant pulses with the Effelsberg telescope at a center frequency of 8.35 GHz show distinct spectral maxima within our observational bandwidth of 500 MHz for individual pulses. Their narrow band components appear to be brighter at higher frequencies (8.6 GHz) than at lower ones (8.1 GHz). Moreover, there is an evidence for spectral evolution within and between those structures. High frequency features occur earlier than low frequency ones. Strong plasma turbulence might be a feasible mechanism for the creation of the high energy densities of ∼6.7 × 10⁴ erg cm⁻³ and brightness temperatures of ∼10³¹ K.     

Numerical MHD modeling of the solar driver for the destabilization of a coronal helmet streamer

We investigate the forms of the solar driver which cause the destabilization of helmet streamers. Two forms of solar drivers are considered; (i) emergence of a flux-rope from sub-photospheric levels and (ii) application of a photospheric shear motion to a streamer-flux rope system. Numerical results showed that both cases exhibit the characteristics of commonly observed coronal mass ejections (CMEs), but the propagation speed of the CME is higher than the background solar wind speed when the solar driver is the emerging magnetic flux and is the same as the solar wind speed when the photospheric shear is used as the solar driver. The energy constraint allowing the magnetic field transition from a closed to an open configuration is also addressed.     

Propagation and interaction of interplanetary transient disturbances. Numerical simulations

We study the heliocentric evolution of ICME-like disturbances and their associated transient forward shocks (TFSs) propagating in the interplanetary (IP) medium comparing the solutions of a hydrodynamic (HD) and magnetohydrodynamic (MHD) models using the ZEUS-3D code [Stone, J.M., Norman, M.L., 1992. Zeus-2d: a radiation magnetohydrodynamics code for astrophysical flows in two space dimensions. i – the hydrodynamic algorithms and tests. Astrophysical Journal Supplement Series 80, 753–790]. The simulations show that when a fast ICME and its associated IP shock propagate in the inner heliosphere they have an initial phase of about quasi-constant propagation speed (small deceleration) followed, after a critical distance (deflection point), by an exponential deceleration. By combining white light coronograph and interplanetary scintillation (IPS) measurements of ICMEs propagating within 1 AU [Manoharan, P.K., 2005. Evolution of coronal mass ejections in the inner heliosphere: a study using white-light and scintillation images. Solar Physics 235 (1–2), 345–368], such a critical distance and deceleration has already been inferred observationally. In addition, we also address the interaction between two ICME-like disturbances: a fast ICME 2 overtaking a previously launched slower ICME 1. After interaction, the leading ICME 1 accelerates and the tracking ICME 2 decelerates and both ICMEs tend to arrive at 1 AU having similar speeds. The 2-D HD and MHD models show similar qualitative results for the evolution and interaction of these disturbances in the IP medium.     

Numerical modeling of inhomogeneous orographic wave influence on planetary waves in the middle atmosphere

Parameterization of dynamical and thermal effects of stationary orographic gravity waves (OGWs) generated by the Earth’s surface topography is incorporated into a numerical model of general circulation of the middle and upper atmosphere. Responses of atmospheric general circulation and characteristics of planetary waves at altitudes from the troposphere up to the thermosphere to the effects of OGWs propagating from the earth surface are studied. Changes in atmospheric circulation and amplitudes of planetary waves due to variations of OGW generation and propagation in different seasons are considered. It is shown that during solstices the main OGW dynamical and heat effects occur in the middle atmosphere of winter hemispheres, where changes in planetary wave amplitudes due to OGWs may reach up to 50%. During equinoxes OGW effects are distributed more homogeneously between northern and southern hemispheres.     

The neutral gas in the environs of the Geminga gamma-ray pulsar

We present a high-resolution (24″) study of the HI interstellar gas distribution around the radio-quiet neutron star Geminga. Based on very large array and MPIfR Effelsberg telescope data, we analyzed a 40′ × 40′ field around Geminga. These observations have revealed the presence of a neutral gas shell, 0.4 pc in radius, with an associated HI mass of 0.8M_⊙, which surrounds Geminga at a radial velocity compatible with the kinematical distance of the neutron star. In addition, morphological agreement is observed between the internal face of the HI shell and the brightest structure of Geminga’s tail observed in X-rays. We explore the possibility that this morphological agreement is the result of a physical association.     

Numerical modeling of propagation of breaking nonlinear acoustic-gravity waves from the lower to the upper atmosphere

Acoustic-gravity waves (AGWs) observed in the upper atmosphere may be generated near the Earth’s surface due to a variety of meteorological sources. Two-dimensional simulations of vertical propagation and breaking of nonlinear AGWs in the atmosphere are performed. Forcing near the Earth’s surface is used as the AGW source in the model. We use a numerical method based on finite-difference analogues of fundamental conservation laws for solving atmospheric hydrodynamic equations. This approach selects physically correct generalized solutions of the wave hydrodynamic equations. Numerical simulations are performed in a representative region of the Earth’s atmosphere up to altitude 500 km. Vertical profiles of temperature, density, molecular viscosity and heat conductivity were taken from the standard atmosphere model MSIS-90 for January. Calculations were made for different amplitudes and frequencies of lower boundary wave forcing. It is shown that after activating the tropospheric wave forcing, the initial pulse of AGWs may very quickly propagate to altitudes of 100 km and above and relatively slowly dissipate due to molecular viscosity and heat conduction. This may increase the role of transient nonstationary waves in effective energy transport and variations of atmospheric parameters and gas admixtures in a broad altitude range.