We use the multi-spacecraft capabilities of the Cluster and THEMIS missions to show that two types of foreshock may be detected in spacecraft data. One is the global foreshock that appears upstream of the Earth's quasi-parallel bow-shock under steady or variable interplanetary magnetic field. Another type is a traveling foreshock that is bounded by two rotational discontinuities in the interplanetary magnetic field and propagates along the bow-shock. Foreshock compressional boundaries are found at the edges of both types of foreshock. We show that isolated foreshock cavities are a subset of the traveling foreshock that form when two bounding rotational discontinuities are so close that the ultra-low frequency waves do not develop in the region between them. We also report observations of a spontaneous hot flow anomaly inside a traveling foreshock. This means that other phenomena, such as foreshock cavitons, may also exist inside this type of foreshock. In the second part of this work we present statistical properties of phenomena related to the foreshock, namely foreshock cavities, cavitons, spontaneous hot flow anomalies and foreshock compressional boundaries. We show that spontaneous hot flow anomalies are the most depleted transient structures in terms of the B-field and plasma density inside them and that the foreshock compressional boundaries and foreshock cavities are closely related structures.

We show for the first time that different types of suprathermal ion distributions may exist upstream of a single interplanetary shock. ACE and the two ARTEMIS satellites observed a shock on 8 October 2013. The ARTEMIS P1 and P2 spacecraft first observed field-aligned ions (P1) and gyrating ions (P2) arriving from the shock. These were followed by intermediate ions and later by a diffuse population. At the location of the P2 the shock exhibited an Alfv\'enic Mach number of M$_A$=5.7 and was marginally quasi-perpendicular, ($\theta_{Bn}$=47$^\circ$). At P1 spacecraft the shock was weaker (M$_A$=4.9) and more perpendicular ($\theta_{Bn}$=61$^\circ$). Consequently the observed suprathermal ion and ultra low frequency wave properties were somewhat different. At P2 the ULF waves are more intense and extend farther upstream from the shock. The energies of field aligned and gyrating ions in the shock rest frame were $\sim$20~keV, which is much more than in the case of the stronger (M$_A$=6-7) Earth's bow-shock, where they are less than 10~keV

We perform the first statistical study of the effects of the interaction of suprathermal electrons with narrow-band whistler mode waves in the solar wind. We show that this interaction does occur and that it is associated with enhanced widths of the so called strahl component. The latter is directed along the inter- planetary magnetic field away from the Sun. We do the study by comparing the strahl pitch angle widths in the solar wind at 1AU in the absence of large scale discontinuities and transient structures, such as interplanetary shocks, interplanetary coronal mass ejections, stream interaction regions, etc. during times when the whistler mode waves were present and when they were absent. This is done by using the data from two Cluster instruments: STAFF data in frequency range between ~0.1 Hz and ~200 Hz were used for determining the wave properties and PEACE datasets at twelve central energies between ~57 eV (equivalent to ~10 typical electron thermal energies in the solar wind, E_T ) and ~676 eV (~113 E_T ) for pitch angle measurements. Statistical analysis shows that during the inter- vals with the whistler waves the strahl component on average exhibits pitch angle widths between 2 and 12 degrees larger than during the intervals when these waves are not present. The largest difference is obtained for the electron central energy of ~344 eV (~57 E_T ).

In this work we study Herbig-Haro objects located in the region around the head of the cometary globule CG 30. Two sets of optical images are presented. The first set was obtained with the 3.5 m New Technology Telescope in 1995 in three emission lines: Halpha, [SII]6731,6716 A and [OII]3729 A. The second set is an Halpha image of the CG 30/31/38 complex obtained in 2006 with the 8 m Subaru telescope. A proper motion study of the HH objects in the region was performed using the Halpha images from both epochs. Due to the high resolution of our images we were able to, for the first time, resolve the HH 120 object into ten knots and measure proper motions for some of them. We discover several new HH objects and a large bipolar jet, HH 950, emerging from the head of CG 30. We suggest that two previously known submillimeter sources are the driving sources for the HH 120 and HH 950 flows.

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).

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.

We present 3D gasdynamic simulations of the Herbig Haro object HH 555. HH 555 is a bipolar jet emerging from the tip of an elephant trunk entering the Pelican Nebula from the adjacent molecular cloud. Both beams of HH 555 are curved away from the center of the H II region. This indicates that they are being deflected by a side-wind probably coming from a star located inside the nebula or by the expansion of the nebula itself. HH 555 is most likely an irradiated jet emerging from a highly embedded protostar, which has not yet been detected. In our simulations we vary the incident photon flux, which in one of our models is equal to the flux coming from a star 1 pc away emitting 5x10^48 ionizing (i. e., with energies above the H Lyman limit) photons per second. An external, plane-parallel flow (a ``side-wind'') is coming from the same direction as the photoionizing flux. We have made four simulations, decreasing the photon flux by a factor of 10 in each simulation. We discuss the properties of the flow and we compute Halpha emission maps (integrated along lines of sight). We show that the level of the incident photon flux has an important influence on the shape and visibility of the jet. If the flux is very high, it causes a strong evaporation of the neutral clump, producing a photoevaporated wind traveling in the direction opposite to the incident flow. The interaction of the two flows creates a double shock ``working surface'' around the clump protecting it and the jet from the external flow. The jet only starts to curve when it penetrates through the working surface.