Observations of several protostellar jets show systematic differences in radial velocity transverse to the jet propagation direction, which have been interpreted as evidence of rotation in the jets. In this paper we discuss the origin of these velocity shifts, and show that they could be originated by rotation in the flow, or by side to side asymmetries in the shock velocity, which could be due to asymmetries in the jet ejection velocity/density or in the ambient medium. For typical poloidal jet velocities (~ 100-200 km/s), an asymmetry >~ 10% can produce velocity shifts comparable to those observed. We also present three dimensional numerical simulations of rotating, precessing and asymmetric jets, and show that, even though for a given jet there is a clear degeneracy between these effects, a statistical analysis of jets with different inclination angles can help to distinguish between the alternative origins of transverse velocity shifts. Our analysis indicate that side to side velocities asymmetries could represent an important contribution to transverse velocity shifts, being the most important contributor for large jet inclination angles (with respect the the plane of the sky), and can not be neglected when interpreting the observations.

A very important issue concerning protostellar jets is the mechanism behind their formation. Obtaining information on the region at the base of a jet can shed light into the subject and some years ago this has been done through a search for a rotational signature at the jet line spectrum. The existence of such signatures, however, remains controversial. In order to contribute to the clarification of this issue, in this paper we show that the Principal Component Analysis (PCA) can potentially help to distinguish between rotation and precession effects in protostellar jet images. We apply the PCA to synthetic spectro-imaging datacubes generated as an output of numerical simulations of protostellar jets. In this way we generate a benchmark to which a PCA diagnostics of real observations can be confronted. Using the computed emission line profiles for [O I]6300A and [S II]6716A, we recover and analyze the effects of rotation and precession in tomograms generated by PCA. We show that different combinations of the eigenvectors can be used to enhance and to identify the rotation features present in the data. Our results indicate that the PCA can be useful for disentangling rotation from precession in jets with an inclination of the jet with respect to the plane of the sky as high as 45 degrees. We have been able to recover the initially imposed rotation jet profile for models at moderate inclination angle (< 15 degrees) and without precession (abridged).

We present new optical observations of the HH 111 Herbig-Haro jet using the Gemini Multi Object Spectrograph in its Integral Field Unit mode. Eight fields of 5" x 3.5" have been positioned along and across the HH 111 jet, covering the spatial region from knot E to L in HH 111 (namely, knots E, F, G, H, J, K and L). We present images and velocity channel maps for the [O I] 6300+6360, Halpha, [N II] 6548+6583 and [S II] 6716+6730 lines, as well as for the [S II]6716/6730 line ratio. We find that the HH 111 jet has an inner region with lower excitation and higher radial velocity, surrounded by a broader region of higher excitation and lower radial velocity. Also, we find higher electron densities at lower radial velocities. These results imply that the HH 111 jet has a fast, axial region with lower velocity shocks surrounded by a lower velocity sheath with higher velocity shocks.

The EUV photoionizing radiation and FUV dissociating radiation from newly born stars photoevaporate their parental neutral cloud, leading to the formation of dense clumps that could eventually form additional stars. We study the effects of including a photodissociating FUV flux in models of the fragmentation of a photoevaporating, self-gravitating molecular cloud. We compute 3D simulations of the interaction of an inhomogeneous, neutral, self-gravitating cloud with external EUV and FUV radiation fields, and calculate the number of collapsing clumps and their mass. We find that the presence of an outer photodissociation region has an important effect on the formation of dense structures due to the expansion of an HII region. In particular, including a FUV field leads to the earlier formation of a larger number of dense clumps, which might lead to the formation of more stars.

This article combines new CCD polarimetric data with previous information about protostellar objects in a search for correlations involving the interstellar magnetic field. Specifically, we carried out an optical polarimetric study of a sample of 28 fields of 10 X 10 arcmin^2 located in the neighborhood of protostellar jets and randomly spread over the Galaxy. The polarimetry of a large number of field stars is used to estimate both the average and dispersion of the interstellar magnetic field (ISMF) direction in each region. The results of the applied statistical tests are as follows. Concerning the alignment between the jet direction and the interstellar magnetic field, the whole sample does not show alignment. There is, however, a statistically significant alignment for objects of Classes 0 and I. Regarding the interstellar magnetic field dispersion, our sample presents values slightly larger for regions containing T Tauri objects than for those harboring younger protostars. Moreover the ISMF dispersion in regions containing high-mass objects tends to be larger than in those including only low-mass protostars. In our sample, the mean interstellar polarization as a function of the average interstellar extinction in a region reaches a maximum value around 3% for A(V) = 5, after which it decreases. Our data also show a clear correlation of the mean value of the interstellar polarization with the dispersion of the interstellar magnetic field: the larger the dispersion, the smaller the polarization. Based on a comparison of our and previous results, we suggest that the dispersion in regions forming stars is larger than in quiescent regions.

In this work we study the influence of the UV radiation field of a massive star on the evolution of a disklike mass of gas and dust around a nearby star. This system has similarities with the Orion proplyds. We study disks with different inclinations and distances from the source, performing 3D numerical simulations. We use the YGUAZ\'U-A adaptative grid code modified to account for EUV/FUV fluxes and non-spherical mass distributions. We treat H and C photoionization to reproduce the ionization fronts and photodissociation regions observed in proplyds. We also incorporate a wind from the ionizing source, to investigate the formation of the bow shock observed in several proplyds. Our results show that a photoevaporated wind propagates from the disk surface and becomes ionized after an ionization front (IF) seen as a bright peak in Ha maps. We follow the development of an HI region inside the photoevaporated wind which corresponds to a photodissociated region (PDR) for most of our models, except those without a FUV flux. For disks that are at a distance from the source d \geq 0.1 pc, the PDR is thick and the IF is detached from the disk surface. In contrast, for disks that are closer to the source, the PDR is thin and not resolved in our simulations. The IF then coincides with the first grid points of the disk that are facing the ionizing photon source. In both cases, the photoevaporated wind shocks (after the IF) with the wind that comes from the ionizing source, and this interaction region is bright in Ha. Our 3D models produce two emission features: a hemispherically shaped structure (associated with the IF) and a detached bow shock where both winds collide. A photodissociated region develops in all of the models exposed to the FUV flux. More importantly, disks with different inclinations with respect to the ionizating source have relatively similar photodissociation regions. (abridged)

We present 3D hydrodynamical simulations of a precessing jet propagating inside a supernova remnant (SNR) shell, particularly applied to the W50-SS433 system in a search for the origin of its peculiar elongated morphology. Several runs were carried out with different values for the mass loss rate of the jet, the initial radius of the SNR, and the opening angle of the precession cone. We found that our models successfully reproduce the scale and morphology of W50 when the opening angle of the jets is set to 10$\degr$ or if this angle linearly varies with time. For these models, more realistic runs were made considering that the remnant is expanding into an interstellar medium (ISM) with an exponential density profile (as HI observations suggest). Taking into account all these ingredients, the large scale morphology of the W50-SS 433 system, including the asymmetry between the lobes (formed by the jet-SNR interaction), is well reproduced.

We present optical linear polarimetry in the line of sight to HH135/HH136. The polarimetry of the field stars reveals two populations: one corresponds to a foreground interstellar component; the other originates in the interstellar medium in the vicinity of the Herbig-Haro pair and, therefore, can be used to study the magnetic field in the star forming region. Its direction is aligned with the jet of HH135/HH136, which could be an indication that the interstellar magnetic field is important in the outflow collimation. The interstellar magnetic field magnitude was estimated to be of order 90 uG. According to recent numerical simulations, an interstellar magnetic field of such strength can be important in the definition of the outflow direction. There is also evidence that the associated dark cloud has an elongation parallel to the magnetic field. Our image polarimetry of the extended emission associated with HH135/HH136 shows a centro-symmetric pattern pointing to the knot E of HH136. Previous near infrared polarimetry traces a different illumination center, namely IRAS 11101-5829 - the probable exciting source of the system. This discrepancy can be explained if the YSO emission is completely blocked in optical wavelengths and the dominant optical source in the region is the knot E, whose nature is uncertain. A discussion of the spectral energy distributions of HH136-E and IRAS 11101-5829 is presented.

Neutral clumps immersed in HII regions are frequently found in star formation regions. We investigate here the formation of tail of neutral gas, which are not reached by the direct ionizing flux coming from two massive stars, using both an analytical approximation, that allows us to estimate the shadow geometry behind the clumps for different initial geometric configurations, and three-dimensional numerical simulations. We found a good agreement between both approaches to this theoretically interesting problem. A particularly important application could be the proplyds that are found in the Trapezium cluster in Orion, which are being photoevaporated primordially by the O stars theta-1 Ori C and theta-2 Ori A.

We study in the present paper how Hall effect modifies the quenching process of the electromotive force (e.m.f.) in Mean Field Dynamo (MFD) theories. We write down the evolution equations for the e.m.f. and for the large and small scale magnetic helicity, treat Hall effect as a perturbation and integrate the resulting equations assuming boundary conditions such that the total divergencies vanish. For force-free large scale magnetic fields, Hall effect acts by coupling the small scale velocity and magnetic fields. For the range of parameters considered, the overall effect is a stronger quenching of the e.m.f. than in standard MHD and a damping of the inverse cascade of magnetic helicity. In astrophysical environments characterized by the parameters considered here, Hall effect would produce an earlier quenching of the e.m.f. and consequently a weaker large scale magnetic field.

Using the Yguazu-a three-dimensional hydrodynamic code, we have computed a set of numerical simulations of heavy, supersonic, radiatively cooling jets including variabilities in both the ejection direction (precession) and the jet velocity (intermittence). In order to investigate the effects of jet rotation on the shape of the line profiles, we also introduce an initial toroidal rotation velocity profile, in agreement with some recent observational evidence found in jets from T Tauri stars which seems to support the presence of a rotation velocity pattern inside the jet beam, near the jet production region. Since the Yguazu-a code includes an atomic/ionic network, we are able to compute the emission coefficients for several emission lines, and we generate line profiles for the H, [O I]6300, [S II]6716 and [N II]6548 lines. Using initial parameters that are suitable for the DG Tau microjet, we show that the computed radial velocity shift for the medium-velocity component of the line profile as a function of distance from the jet axis is strikingly similar for rotating and non-rotating jet models. These findings lead us to put forward some caveats on the interpretation of the observed radial velocity distribution from a few outflows from young stellar objects, and we claim that these data should not be directly used as a doubtless confirmation of the magnetocentrifugal wind acceleration models.

We present high spatial resolution spectroscopic observations of the proplyd 167-317 (LV2) near the Trapezium cluster in the Orion nebula, obtained during the System Verification run of the Gemini Multi Object Spectrograph (GMOS) Integral Field Unit (IFU) at the Gemini South Observatory. We have detected 38 forbidden and permitted emission lines associated with the proplyd and its redshifted jet. We have been able to detect three velocity components in the profiles of some of these lines: a peak with a 28-33 km/s systemic velocity that is associated with the photoevaporated proplyd flow, a highly redshifted component associated with a previously reported jet (which has receding velocities of about 80-120 km/s with respect to the systemic velocity and is spatially distributed to the southeast of the proplyd) and a less obvious, approaching structure, which may possibly be associated with a faint counter-jet with systemic velocity of (-75 +/- 15) km/s. We find evidences that the redshifted jet has a variable velocity, with slow fluctuations as a function of the distance from the proplyd. We present several background subtracted, spatially distributed emission line maps and we use this information to obtain the dynamical characteristics over the observed field. Using a simple model and with the extinction corrected Halpha fluxes, we estimate the mass loss rate for both the proplyd photoevaporated flow and the redshifted microjet, obtaining (6.2 +/- 0.6) x 10^{-7} M_sun/year and (2.0 +/- 0.7) x 10^{-8} M_sun/year, respectively.

We here present results of three-dimensional Smoothed Particle hydro and magnetohydrodynamics simulations of rotating jets, also including the effects of radiative cooling, precession and velocity variability. Using initial conditions and parameters which are particularly suitable for the DG Tau microjet, we have been able to approximately reproduce its complex knotty morphology and kinematics. We have also obtained radial velocity maps which are in good agreement with the data obtained by Bacciotti et al., thus indicating that their interpretation that the DG Tau microjet is rotating is correct. Finally, we have found that a magnetic field of the order of ~ 0.5 mG is sufficient to collimate the jet against the lateral expansion that is caused by the centrifugal forces.

(Abridged) We here investigate, by means of fully 3-D Smoothed Particle Magnetohydrodynamic numerical simulations, the effects of magnetic fields on overdense, radiatively cooling, pulsed jets, using different initial magnetic field topologies and strengths ($B \simeq 260 \mu$G-0). The relative differences that have been previously detected in 2-D simulations involving distinct magnetic field configurations are diminished in the 3-D flows. While the presence of toroidal magnetic components can modify the morphology close to the jet head inhibiting its fragmentation in the early jet evolution, as previously reported in the literature, the impact of the pulsed-induced internal knots causes the appearance of a complex morphology at the jet head (as required by the observations of H-H jets) even in the MHD jet models with toroidal components. The detailed structure and emission properties of the internal working surfaces can be also significantly altered by the presence of magnetic fields. The increase of the magnetic field strength improves the jet collimation, and amplifies the density (by factors up to 1.4, and 4) and the H\alpha$ intensity (by factors up to 4, and 5) behind the knots of jets with helical field and $\beta \simeq 1-0.1$ (respectively), relative to a non magnetic jet. As a consequence, the corresponding $I_{[SII]}}/I_{H}\alpha}$ ratio (which is frequently used to determine the excitation level of HH objects) can be decreased in the MHD models with toroidal components relative to non-magnetic calculations. We also find that the helical mode of the K-H instability can be triggered in MHD models with helical magnetic fields, causing some jet wiggling. No evidence for the formation of the nose cones is found in the 3-D flows, nor even in the $\beta \simeq 0.1$ case.

We here present the first results of fully three-dimensional (3-D) MHD simulations of radiative cooling pulsed (time-variable) jets for a set of parameters which are suitable for protostellar outflows. Considering different initial magnetic field topologies in approximate $equipartition$ with the thermal gas, i.e., (i) a longitudinal, and (ii) a helical field, both of which permeating the jet and the ambient medium; and (iii) a purely toroidal field permeating only the jet, we find that the overall morphology of the pulsed jet is not very much affected by the presence of the different magnetic field geometries in comparison to a nonmagnetic calculation. Instead, the magnetic fields tend to affect essentially the detailed structure and emission properties behind the shocks at the head and at the pulse-induced internal knots, particularly for the helical and toroidal geometries. In these cases, we find, for example, that the $H_\alpha$ emissivity behind the internal knots can be about three to four times larger than that of the purely hydrodynamical jet. We also find that some features, like the nose cones that often develop at the jet head in 2-D calculations involving toroidal magnetic fields, are smoothed out or absent in the 3-D calculations.

We here review the results of recent studies of magnetic field effects on the structure and evolution of overdense, radiative cooling, steady and pulsed jets which were carried out with the help of fully three-dimensional MHD simulations. Three initial magnetic field topologies are considered: (i) a helical and (ii) a longitudinal field, both of which permeate the jet and the ambient medium, and (iii) a pure toroidal field in the jet. Using a set of parameters which is particularly suitable for protostellar jets, we examine the emission structure of the internal knots and the leading working surface and compare with a nonmagnetic baseline and also with 2-D MHD calculations.

We investigate the effect of magnetic fields on the propagation dynamics and morphology of overdense, radiatively cooling, supermagnetosonic jets, with the help of fully three-dimensional SPMHD simulations. Evaluated for a set of parameters which are mainly suitable for protostellar jets (with density ratios between the jet and the ambient medium 3-10, and ambient Mach number ~ 24), these simulations are also compared with baseline non-magnetic and adiabatic calculations. We find that, after amplification by compression and re-orientation in nonparallel shocks at the working surface, the magnetic field that is carried backward with the shocked gas into the cocoon improves the jet collimation relative to the purely hydrodynamic (HD) systems. Low-amplitude, approximately equally spaced internal shocks (which are absent in the HD systems) are produced by MHD K-H reflection pinch modes. The longitudinal field geometry also excites non-axisymmetric helical modes which cause some beam wiggling. The strength and amount of these modes are, however, reduced (by ~ twice) in the presence of radiative cooling relative to the adiabatic cases. Besides, a large density ratio between the jet and the ambient medium also reduces, in general, the number of the internal shocks. As a consequence, the weakness of the induced internal shocks makes it doubtful that the magnetic pinches could produce by themselves the bright knots observed in the overdense, radiatively cooling protostellar jets.

In the present work, with the help of fully 3-D smooth particle magnetohydodynamic (SPMHD) simulations, we explore the nonlinear effects of magnetic fields (in close equipartition with the gas) in the structure of radiatively cooling, overdense jets and test their importance on the dynamics of protostellar jets.

In this work, we sumarize the results of three-dimensional smooth particle magnetohydrodynamic numerical simulations of supermagnetosonic, radiatively cooling jets.

We present the results of 3-D SPMHD numerical simulations of supermagnetosonic, overdense, radiatively cooling jets. Two initial magnetic configurations are considered: (i) a helical and (ii) a longitudinal field. We find that magnetic fields have important effects on the dynamics and structure of radiative cooling jets, especially at the head. The presence of a helical field suppresses the formation of the clumpy structure which is found to develop at the head of purely hydrodynamical jets. On the other hand, a cooling jet embedded in a longitudinal magnetic field retains clumpy morphology at its head. This fragmented structure resembles the knotty pattern commonly observed in HH objects behind the bow shocks of HH jets. This suggests that a strong (equipartition) helical magnetic field configuration is ruled out at the jet head. Therefore, if strong magnetic fields are present, they are probably predominantly longitudinal in those regions. In both magnetic configurations, we find that the confining pressure of the cocoon is able to excite short-wavelength MHD K-H pinch modes that drive low-amplitude internal shocks along the beam. These shocks are not strong however, and it likely that they could only play a secondary role in the formation of the bright knots observed in HH jets.

We present fully three-dimensional hydrodynamical and magnetohydrodynamical simulations of supersonic, radiative cooling jets propagating into homogeneous and stratified environments based on Smoothed Particle Hydrodynamics technique.