A latitude-dependent wind model for Mira's cometary head

We present a 3D numerical simulation of the recently discovered cometary structure produced as Mira travels through the galactic ISM. In our simulation, we consider that Mira ejects a steady, latitude-dependent wind, which interacts with a homogeneous, streaming environment. The axisymmetry of the problem is broken by the lack of alignment between the direction of the relative motion of the environment and the polar axis of the latitude-dependent wind. With this model, we are able to produce a cometary head with a ``double bow shock'' which agrees well with the structure of the head of Mira's comet. We therefore conclude that a time-dependence in the ejected wind is not required for reproducing the observed double bow shock.Comment: 4 pages, 4 figures, accepted for publication in ApJ

Influence of He++ and Shock Geometry on Interplanetary Shocks in the Solar Wind: 2D Hybrid Simulations

After protons, alpha particles (He$^{++}$) are the most important ion species in the solar wind, constituting typically about 5\% of the total ion number density. Due to their different charge-to-mass ratio protons and He$^{++}$ particles are accelerated differently when they cross the electrostatic potential in a collisionless shock. This behavior can produce changes in the velocity distribution function (VDF) for both species generating anisotropy in the temperature which is considered to be the energy source for various phenomena such as ion cyclotron and mirror mode waves. How these changes in temperature anisotropy and shock structure depend on the percentage of He$^{++}$ particles and the geometry of the shock is not completely understood. In this paper we have performed various 2D local hybrid simulations (particle ions, massless fluid electrons) with similar characteristics (e.g., Mach number) to interplanetary shocks for both quasi-parallel and quasi-perpendicular geometries self-consistently including different percentages of He$^{++}$ particles. We have found changes in the shock transition behavior as well as in the temperature anisotropy as functions of both the shock geometry and He$^{++}$ particle abundance: The change of the initial $\theta_{Bn}$ leads to variations of the efficiency with which particles can escape to the upstream region facilitating or not the formation of compressive structures in the magnetic field that will produce increments in perpendicular temperature. The regions where both temperature anisotropy and compressive fluctuations appear tend to be more extended and reach higher values as the He$^{++}$ content in the simulations increases.Data set in h5 format corresponding to each panel of the figures of the publicatio

Numerical simulations of radiative magnetized Herbig-Haro jets: the influence of pre-ionization from X-rays on emission lines

We investigate supersonic, axisymmetric magnetohydrodynamic (MHD) jets with a time-dependent injection velocity by numerical simulations with the PLUTO code. Using a comprehensive set of parameters, we explore different jet configurations in the attempt to construct models that can be directly compared to observational data of microjets. In particular, we focus our attention on the emitting properties of traveling knots and construct, at the same time, accurate line intensity ratios and surface brightness maps. Direct comparison of the resulting brightness and line intensity ratios distributions with observational data of microjets shows that a closer match can be obtained only when the jet material is pre-ionized to some degree. A very likely source for a pre-ionized medium is photoionization by X-ray flux coming from the central object.Comment: Accepted for publication in Ap

Variability of the Magnetic Field Power Spectrum in the Solar Wind at Electron Scales

At electron scales, the power spectrum of solar-wind magnetic fluctuations can be highly variable and the dissipation mechanisms of the magnetic energy into the various particle species is under debate. In this paper, we investigate data from the Cluster mission's STAFF Search Coil magnetometer when the level of turbulence is sufficiently high that the morphology of the power spectrum at electron scales can be investigated. The Cluster spacecraft sample a disturbed interval of plasma where two streams of solar wind interact. Meanwhile, several discontinuities (coherent structures) are seen in the large-scale magnetic field, while at small scales several intermittent bursts of wave activity (whistler waves) are present. Several different morphologies of the power spectrum can be identified: (1) two power laws separated by a break, (2) an exponential cutoff near the Taylor shifted electron scales, and (3) strong spectral knees at the Taylor shifted electron scales. These different morphologies are investigated by using wavelet coherence, showing that, in this interval, a clear break and strong spectral knees are features that are associated with sporadic quasi parallel propagating whistler waves, even for short times. On the other hand, when no signatures of whistler waves at similar to 0.1-0.2f(ce) are present, a clear break is difficult to find and the spectrum is often more characteristic of a power law with an exponential cutoff.Peer reviewe

Shaping point- and mirror-symmetric proto-planetary nebulae by the orbital motion of the central binary system

We present 3D hydrodynamical simulations of a jet launched from the secondary star of a binary system inside a proto-planetary nebula. The secondary star moves around the primary in a close eccentric orbit. From the gasdynamic simulations we compute synthetic [NII] 6583 emission maps. Different jet axis inclinations with respect to the orbital plane, as well as different orientations of the flow with respect to the observer are considered. For some parameter combinations, we obtain structures that show point- or mirror-symmetric morphologies depending on the orientation of the flow with respect to the observer. Furthermore, our models can explain some of the emission distribution asymmetries that are summarized in the classification given by Soker & hadar (2002).Comment: 15 pages, 3 figures, 2 tables, Accepted in Apj Letter

Jets Downstream of Collisionless Shocks

The magnetosheath flow may take the form of large amplitude, yet spatially localized, transient increases in dynamic pressure, known as "magnetosheath jets" or "plasmoids" among other denominations. Here, we describe the present state of knowledge with respect to such jets, which are a very common phenomenon downstream of the quasi-parallel bow shock. We discuss their properties as determined by satellite observations (based on both case and statistical studies), their occurrence, their relation to solar wind and foreshock conditions, and their interaction with and impact on the magnetosphere. As carriers of plasma and corresponding momentum, energy, and magnetic flux, jets bear some similarities to bursty bulk flows, which they are compared to. Based on our knowledge of jets in the near Earth environment, we discuss the expectations for jets occurring in other planetary and astrophysical environments. We conclude with an outlook, in which a number of open questions are posed and future challenges in jet research are discussed.Peer reviewe

ULF Wave Transmission Across Collisionless Shocks : 2.5D Local Hybrid Simulations

We study the interaction of upstream ultralow frequency (ULF) waves with collisionless shocks by analyzing the outputs of 11 2D local hybrid simulation runs. Our simulated shocks have Alfvenic Mach numbers between 4.29 and 7.42 and their theta BN angles are 15 degrees, 30 degrees, 45 degrees, and 50 degrees. The ULF wave foreshocks develop upstream of all of them. The wavelength and the amplitude of the upstream waves exhibit a complex dependence on the shock's MA and theta BN. The wavelength positively correlates with both parameters, with the dependence on theta BN being much stronger. The amplitude of the ULF waves is proportional to the product of the reflected beam velocity and density, which also depend on MA and theta BN. The interaction of the ULF waves with the shock causes large-scale (several tens of upstream ion inertial lengths) shock rippling. The properties of the shock ripples are related to the ULF wave properties, namely their wavelength and amplitude. In turn, the ripples have a large impact on the ULF wave transmission across the shock because they change local shock properties (theta BN, strength), so that different sections of the same ULF wavefront encounter shock with different characteristics. Downstream fluctuations do not resemble the upstream waves in terms the wavefront extension, orientation or their wavelength. However, some features are conserved in the Fourier spectra of downstream compressive waves that present a bump or flattening at wavelengths approximately corresponding to those of the upstream ULF waves. In the transverse downstream spectra, these features are weaker.Peer reviewe

Transient Foreshock Structures Upstream of Mars: Implications of the Small Martian Bow Shock

We characterize the nature of magnetic structures in the foreshock region of Mars associated with discontinuities in the solar wind. The structures form at the upstream edge of moving foreshocks caused by slow rotations in the interplanetary magnetic field (IMF). The solar wind plasma density and the IMF strength noticeably decrease inside the structures' core, and a compressional shock layer is present at their sunward side, making them consistent with foreshock bubbles (FBs). Ion populations responsible for these structures include backstreaming ions that only appear within the moving foreshock, and accelerated reflected ions from the quasi-perpendicular bow shock. Both ion populations accumulate near the upstream edge of the moving foreshock which facilitates FB formation. Reflected ions with hybrid trajectories that straddle between the quasi-perpendicular and quasi-parallel bow shocks during slow IMF rotations contribute to formation of foreshock transients.Comment: Submitted to Geophysical Research Letter