Numerical Simulations of HH 555

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.Comment: 14 pages, 4 figures, accepted by Ap

Side-entrainment in a jet embedded in a sidewind

Numerical simulations of HH jets never show side-entrainment of environmental material into the jet beam. This is because the bow shock associated with the jet head pushes the surrounding environment into a dense shell, which is never in direct contact with the sides of the jet beam. We present 3D simulations in which a side-streaming motion (representing the motion of the outflow source through the surrounding medium) pushes the post-bow shock shell into direct contact with the jet beam. This is a possible mechanism for modelling well collimated "molecular jets" as an atomic/ionic flow which entrains molecules initially present only in the surrounding environment.Comment: 8 pages, 12 figures, 1 table, accepted for publication in Ap

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

A 3-mode, Variable Velocity Jet Model for HH 34

Variable ejection velocity jet models can qualitatively explain the appearance of successive working surfaces in Herbig-Haro (HH) jets. This paper presents an attempt to explore which features of the HH-34 jet can indeed be reproduced by such a model. From previously published data on this object, we find evidence for the existence of a 3-mode ejection velocity variability, and then explore the implications of such a variability. From simple, analytic considerations it is possible to show that the longer period modes produce a modulation on the shorter period modes, resulting in the formation of ``trains'' of multiple knots. The knots observed close to the source of HH-34 could correspond to such a structure. Finally, a numerical simulation with the ejection velocity variability deduced from the HH-34 data is computed. This numerical simulation shows a quite remarkable resemblance with the observed properties of the HH-34 jet.Comment: 28 pages LaTex, 10 postscript figure

A model of Mira's 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 Mira's comet. The simulations were carried out with our new adaptive grid code WALICXE, which is described in detail.Comment: 12 pages, 8 figures (4 in color). Accepted for publication in The Astrophysical Journa

The 3-Dimensional Structure of HH 32 from GMOS IFU Spetroscopy

We present new high resolution spectroscopic observations of the Herbig-Haro object HH 32 from System Verification observations made with the GMOS IFU at Gemini North Observatory. The 3D spectral data covers a 8''.7 x 5''.85 spatial field and 4820 - 7040 Angstrom spectral region centered on the HH~32 A knot complex. We show the position-dependent line profiles and radial velocity channel maps of the Halpha line, as well as line ratio velocity channel maps of [OIII]5007/Halpha, [OI]6300/Halpha, [NII]6583/Halpha, [SII](6716+6730)/Halpha and [SII]6716/6730. We find that the line emission and the line ratios vary significantly on spatial scales of ~1'' and over velocities of ~50 km/s. A ``3/2-D'' bow shock model is qualitatively successful at reproducing the general features of the radial velocity channel maps, but it does not show the same complexity as the data and it fails to reproduce the line ratios in our high spatial resolution maps. The observations of HH 32 A show two or three superimposed bow shocks with separations of ~3'', which we interpret as evidence of a line of sight superposition of two or three working surfaces located along the redshifted body of the HH 32 outflow.Comment: Accepted for Publication in the Astronomical Journal (January 2004

Filaments in Galactic Winds Driven by Young Stellar Clusters

The starburst galaxy M82 shows a system of H$\alpha$-emitting filaments which extend to each side of the galactic disk. We model these filaments as the result of the interaction between the winds from a distribution of Super Stellar Clusters (SSCs). We first derive the condition necessary for producing a radiative interaction between the cluster winds (a condition which is met by the SSC distribution of M82). We then compute 3D simulations for SSC wind distributions which satisfy the condition for a radiative interaction, and also for distributions which do not satisfy this condition. We find that the highly radiative models, that result from the interaction of high metallicity cluster winds, produce a structure of H$\alpha$ emitting filaments, which qualitatively agrees with the observations of the M82, while the non-radiative SSC wind interaction models do not produce filamentary structures. Therefore, our criterion for radiative interactions (which depends on the mass loss rate and the terminal velocity of the SSC winds, and the mean separation between SSCs) can be used to predict whether or not an observed galaxy should have associated H$\alpha$ emitting filaments.Comment: 10 pages, 6 Figures. ApJ Accepted, August 7, 200

Emission lines from rotating proto-stellar jets with variable velocity profiles. I. Three-dimensional numerical simulation of the non-magnetic case

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.Comment: 15 pages, 8 figures. Accepted to publication in Astronomy and Astrophysic

The precession of the giant HH34 outflow: a possible jet deceleration mechanism

The giant jets represent a fundamental trace of the historical evolution of the outflow activity over timescales which are comparable to the accretion time of the outflow sources in their main protostellar phase. The study of such huge jets provides the possibility of retrieving important elements related to the life of the outflow sources. In this paper, we study the role of precession (combined with jet velocity-variability and the resulting enhanced interaction with the surrounding environment) as a deceleration mechanism for giant jets using a numerical approach. We obtain predictions of H alpha intensity maps and position-velocity diagrams from 3D simulations of the giant HH 34 jet (including an appropriate ejection velocity time-variability and a precession of the outflow axis), and we compare them with previously published observations of this object. Our simulations represent a step forward from previous numerical studies of HH objects, in that the use of a 7-level, binary adaptive grid has allowed us to compute models which appropiately cover all relevant scales of a giant jet, from the ~ 100 AU jet radius close to the source to the ~ 1 pc length of the outflow. A good qualitative and quantitative agreement is found between the model predictions and the observations. Moreover, we show that a critical parameter for obtaining a better or worse agreement with the observations is the ratio rho_j/rho_a between the jet and the environmental densities. The implications of this result in the context of the current star formation models are discussed (ABRIDGED).Comment: 19 pages, 8 eps figs.,uses aaspp4; accepted by the Ap