Collisionally excited filaments in HST H$alpha$ and H$beta$ images of HH~1/2

We present new H$alpha$ and H$beta$ images of the HH~1/2 system, and we find that the H$alpha$/H$beta$ ratio has high values in ridges along the leading edges of the HH~1 bow shock and of the brighter condensations of HH~2. These ridges have H$alpha$/H$beta=4to 6$, which is consistent with collisional excitation from the $n=1$ to the $n=3$ and 4 levels of hydrogen in a gas of temperatures $T=1.5to 10times 10^4$~K. This is therefore the first direct proof that the collisional excitation/ionization region of hydrogen right behind Herbig-Haro shock fronts is detected.

An asymmetric jet launching model for the protoplanetary nebula CRL 618

We propose an asymmetrical jet ejection mechanism in order to model the mirror symmetry observed in the lobe distribution of some protoplanetary nebulae (pPNe), such as the pPN CRL 618. 3D hydrodynamical simulations of a precessing jet launched from an orbiting source were carried out including an alternation in the ejections of the two outflow lobes, depending on which side of the precessing accretion disk is hit by the accretion column from a Roche lobe-filling binary companion. Both synthetic optical emission maps and position-velocity (PV) diagrams were obtained from the numerical results with the purpose of carrying out a direct comparison with observations. Depending on the observers point of view, multipolar morphologies are obtained which exhibit a mirror symmetry at large distances from the central source. The obtained lobe sizes and their spatial distribution are in good agreement with the observed morphology of the pPN CRL 618. We also obtain that the kinematic ages of the fingers are similar to those obtained in the observations.

Gravitational drag on a point mass in hypersonic motion within a Gaussian disk

We develop an analytical model for the accretion and gravitational drag on a point mass that moves hypersonically in the midplane of a gaseous disk with a Gaussian vertical density stratification. Such a model is of interest for studying the interaction between a planet and a protoplanetary disk, as well as the dynamical decay of massive black holes in galactic nuclei. The model considers that the flow is ballistic, and gives fully analytical expressions for both the accretion rate onto the point mass, and the gravitational drag it suffers. The expressions are further simplified by taking the limits of a thick, and of a thin disk. The results for the thick disk reduce correctly to those for a uniform density environment (Canto et al. 2011). We find that for a thin disk (small vertical scaleheight compared to the gravitational radius) the accretion rate is proportional to the mass of the moving object and to the surface density of the disk, while the drag force is independent of the velocity of the object. The gravitational deceleration of the hypersonic perturber in a thin disk was found to be independent of its parameters (i.e. mass or velocity) and depends only on the surface mass density of the disk. The predictions of the model are compared to the results of three-dimensional hydrodynamical simulations, with a reasonable agreement.

Proper Motions of Young Stellar Outflows in the Mid-Infrared with Spitzer (IRAC). I. The NGC 1333 region

We use two 4.5micron Spitzer (IRAC) maps of the NGC 1333 region taken over approx. 7 yr interval to determine proper motions of its associated outflows. This is a first, successful attempt at obtaining proper motions of stellars outflow from Spitzer observations. For the outflow formed by the Herbig-Haro objects HH7, 8 and 10, we find proper motions of approx. 9-13 km/s, which are consistent with previously determined optical proper motions of these objects. We determine proper motions for a total of 8 outflows, ranging from approx. 10 to 100 km/s. The derived proper motions show that out of these 8 outflows, 3 have tangential velocities less or equal to 20 km/s. This result shows that a large fraction of the observed outflows have low intrinsic velocities, and that the low proper motions are not merely a projection effect.

Spitzer observations of the HH 1/2 system. The discovery of the counterjet

We present unpublished Spitzer IRAC observations of the HH 1/2 young stellar outflow processed with a high angular resolution deconvolution algorithm that produces sub-arcsecond (approx. 0.6 - 0.8) images. In the resulting mid-infrared images the optically invisible counterjet is detected for the first time. The counterjet is approximately half as bright as the jet at 4.5 micron (the IRAC band that best traces young stellar outflows) and has a length of approx. 10. The NW optical jet itself can be followed back in the mid-IR to the position of the exciting VLA 1 source. An analysis of the IRAC colors indicates that the jet/counterjet emission is dominated by collisionally excited H2 pure rotational lines arising from a medium with a neutral Hydrogen gas density of 1000-2000 per cubic cm and a temperature of 1500 K. The observed jet/counterjet brightness asymmetry is consistent with an intrinsically symmetric outflow with extinction from a dense, circumstellar structure of 6 size (along the outflow axis), and with a mean visual extinction of Av=11 mag.

The precession of the HH 111 flow in the infrared

We present Spitzer IRAC images of the HH 111 outflow, that show a wealth of condensations/knots in both jet and counterjet. Studying the positional distribution of these knots, we find very suggestive evidence of a mirror symmetric pattern in the jet/counterjet flow. We model this pattern as the result of an orbital motion of the jet source around a binary companion. From a fit of an analytic, ballistic model to the observed path of the HH 111 system, we find that the motion in a binary with two approx. 1 Msolar stars (one of them being the HH 111 source), in a circular orbit with a separation of approx. 186 AU would produce the mirror symmetric pattern seen in the outflow.

A model of Miras cometary head/tail entering the Local Bubble

We model the cometary structure around Mira as the interaction of an AGB wind from Mira A, and a streaming environment. Our simulations introduce the following new element: we assume that after 200 kyr of evolution in a dense environment Mira entered the Local Bubble (low density coronal gas). As Mira enters the bubble, the head of the comet expands quite rapidly, while the tail remains well collimated for a 100 kyr timescale. The result is a broad-head/narrow-tail structure that resembles the observed morphology of Miras comet. The simulations were carried out with our new adaptive grid code WALICXE, which is described in detail.

Kinematics of the H2O masers at the centre of the PN K3-35

We have studied the kinematics traced by the water masers located at the centre of the planetary nebula (PN) K3-35, using data from previous Very Large Array (VLA) observations. An analysis of the spatial distribution and line-of-sight velocities of the maser spots allows us to identify typical patterns of a rotating and expanding ring in the position-velocity diagrams, according to our kinematical model. We find that the distribution of the masers is compatible with tracing a circular ring with a ~0.021 arcsec (~100 AU) radius, observed with an inclination angle with respect to the line of sight of 55 degrees. We derive expansion and rotation velocities of 1.4 and 3.1 km/s, respectively. The orientation of the ring projected on the plane of the sky, at PA 158 degrees, is almost orthogonal to the direction of the innermost region of the jet observed in K3-35, suggesting the presence of a disc or torus that may be related to the collimation of the outflow.

A latitude-dependent wind model for Miras 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 Miras comet. We therefore conclude that a time-dependence in the ejected wind is not required for reproducing the observed double bow shock.