We present Spitzer (IRAC) images observations and a VLT 2.1micron image of the HH 212 outflow. We find that this outflow has a strong symmetry, with jet/counterjet knot pairs with Delta x less than 1 arcsec position offsets. We deduce that the jet/counterjet knots are ejected with time differences Delta tau_0 approx. 6 yr and velocity differences Delta v_0~ 2 km/s. We also analyze the deviations of the knot positions perpendicular to the outflow axis, and interpret them in terms of a binary orbital motion of the outflow source. Through this model, we deduce a ~0.7M_solar mass for the outflow source, and a separation of ~80 AU between the components of the binary (assuming equal masses for the two components). Finally, using the IRAC data and the VLT 2.1micron image we have measured the proper motion velocities, obtaining values from 50 to 170km/s.

The bipolar jet HH 212, among the finest collimated jets known, has so far been detected only in near-infrared H$_2$ emission. Here we present deep optical images that show two of the major bow shocks weakly detected in optical [SII] emission, as expected for a bona fide Herbig-Haro jet. We present widefield H$_2$ images which reveal two more bow shocks located symmetrically around the source and along the main jet axis. Additionally, examination of Spitzer 4.5 $\mu$m images reveals yet another bright bow shock further to the north along the jet axis; no corresponding bow shock is seen to the south. In total, the HH 212 flow has an extent of 1050 arcsec, corresponding to a projected dimension of 2.0 pc. HH 212 thus joins the growing group of parsec-scale Herbig-Haro jets. Proper motion measurements indicate a velocity of about 170 km/sec, highly symmetric around the source, with an uncertainty of $\sim$30 km/sec, suggesting a probable age of the giant HH 212 flow of about 7000 yr. The jet is driven by a deeply embedded source, known as IRAS 05413-0104. We draw attention to a Spitzer near- and mid-infrared source, which we call IRS-B, located only 7 arcsec from the driving source, towards the outskirts of the dense cloud core. Infrared photometry and spectroscopy suggests that IRS-B is a K-type star with a substantial infrared excess, except that for an extinction of A$_V$ = 44 the star would have only a weakinfrared excess, and so in principle it could be a K-giant at a distance of about 2 kpc.

During star formation, the accretion disk drives fast MHD winds which usually contain two components, a collimated jet and a radially distributed wide-angle wind. These winds entrain the surrounding ambient gas producing molecular outflows. We report recent observation of 12CO (2-1) emission of the HH 46/47 molecular outflow by the Atacama Large Millimeter/sub-millimeter Array, in which we identify multiple wide-angle outflowing shell structures in both the blue and red-shifted outflow lobes. These shells are highly coherent in position-position-velocity space, extending to >40-50 km/s in velocity and 10^4 au in space with well defined morphology and kinematics. We suggest these outflowing shells are the result of the entrainment of ambient gas by a series of outbursts from an intermittent wide-angle wind. Episodic outbursts in collimated jets are commonly observed, yet detection of a similar behavior in wide-angle winds has been elusive. Here we show clear evidence that the wide-angle component of the HH 46/47 protostellar outflows experiences similar variability seen in the collimated component.

The molecular richness of fast protostellar jets within 20-100 au of their source, despite strong ultraviolet irradiation, remains a challenge for the models investigated so far. We aim to investigate the effect of interaction between a time-variable jet and a surrounding steady disk wind, to assess the possibility of jet chemical enrichement by the wind, and the characteristic signatures of such a configuration. We have constructed an analytic model of a jet bow shock driven into a surrounding slower disk wind in the thin shell approximation. The refilling of the post bow shock cavity from below by the disk wind is also studied. An extension of the model to the case of two or more successive internal working surfaces (IWS) is made. We then compared this analytic model with numerical simulations with and without a surrounding disk wind. We find that at early times (of order the variability period), jet bow shocks travel in refilled pristine disk-wind material, before interacting with the cocoon of older bow shocks. This opens the possibility of bow shock chemical enrichment (if the disk wind is molecular and dusty) and of probing the unperturbed disk wind structure near the jet base. Several distinctive signatures of the presence of a surrounding disk wind are identified, in the bow shock morphology and kinematics. Numerical simulations validate our analytical approach and further show that at large scale, the passage of many jet IWS inside a disk wind produces a stationary V-shaped cavity, closing down onto the axis at a finite distance from the source.

We present a comparison between the time-evolution over the past $\sim 20$ years of the radio continuum and H$\alpha$ emission of HH~1 and 2. We find that the radio continuum and the H$\alpha$ emission of both objects show very similar trends, with HH~1 becoming fainter and HH~2 brightening quite considerably (about a factor of 2). We also find that the $F_{\rm H\alpha}/F_{ff}$ (H$\alpha$ to free-free continuum) ratio of HH~1 and 2 has higher values than the ones typically found in planetary nebulae (PNe) which we interpret as an indication that the H$\alpha$ and free-free emission of HH~1/2 is produced in emitting regions with lower temperatures ($\sim 2000$~K) than the emission of PNe (with $\sim 10^4$~K).

We describe a new method for determining proper motions of extended objects, and a pipeline developed for the application of this method. We then apply this method to an analysis of four epochs of [S~II] HST images of the HH~1 jet (covering a period of $\sim 20$~yr). We determine the proper motions of the knots along the jet, and make a reconstruction of the past ejection velocity time-variability (assuming ballistic knot motions). This reconstruction shows an "acceleration" of the ejection velocities of the jet knots, with higher velocities at more recent times. This acceleration will result in an eventual merging of the knots in $\sim 450$~yr and at a distance of $\sim 80"$ from the outflow source, close to the present-day position of HH~1.

We have analyzed four epochs of H$\alpha$ and [S~II] HST images of the HH~1/2 outflow (covering a time interval from 1994 to 2014) to determine proper motions and emission line fluxes of the knots of HH~2. We find that our new proper motions agree surprisingly well with the motions measured by Herbig \& Jones (1981), although there is partial evidence for a slight deceleration of the motion of the HH~2 knots from 1945 to 2014. We also measure the time-variability of the H$\alpha$ intensities and the [S~II]/H$\alpha$ line ratios, and find that knots H and A have the largest intensity variabilities (in $1994\to 2014$). Knot H (which now dominates the HH~2 emission) has strengthened substantially, while keeping an approximately constant [S~II]/H$\alpha$ ratio. Knot A has dramatically faded, and at the same time has had a substantial increase in its [S~II]/H$\alpha$ ratio. Possible interpretations of these results are discussed.

A significant fraction of OB-type, main-sequence massive stars are classified as runaway and move supersonically through the interstellar medium (ISM). Their strong stellar winds interact with their surroundings where the typical strength of the local ISM magnetic field is about 3.5-7 micro-G, which can result in the formation of bow shock nebulae. We investigate the effects of such magnetic fields, aligned with the motion of the flow, on the formation and emission properties of these circumstellar structures. Our axisymmetric, magneto-hydrodynamical simulations with optically-thin radiative cooling, heating and anisotropic thermal conduction show that the presence of the background ISM magnetic field affects the projected optical emission our bow shocks at Ha and [OIII] lambda 5007 which become fainter by about 1-2 orders of magnitude, respectively. Radiative transfer calculations against dust opacity indicate that the magnetic field slightly diminishes their projected infrared emission and that our bow shocks emit brightly at 60 micron. This may explain why the bow shocks generated by ionizing runaway massive stars are often difficult to identify. Finally, we discuss our results in the context of the bow shock of Zeta Ophiuchi and we support the interpretation of its imperfect morphology as an evidence of the presence of an ISM magnetic field not aligned with the motion of its driving star.

We present high-resolution, long-slit spectroscopic observations of two planetary nebulae with [WC] central stars located near the galactic bulge, M 1-32 and M 3-15. The observations were obtained with the 2.1-m telescope at the Observatorio Astron\'omico Nacional, San Pedro M\'artir. M 1-32 shows wide wings on the base of its emission lines and M 3-15 has two very faint high-velocity knots. In order to model both PNe we built a three-dimensional model consisting of a jet interacting with an equatorially concentrated slow wind, emulating the presence of a dense torus, using the Yguaz\'u hydrodynamical code. From our hydrodynamical models, we obtained position-velocity (PV) diagrams in the [N II]$\lambda$6583 line for comparison with the observations. We find that the spectral characteristics of M 1-32 and M 3-15 can be explained with the same physical model -a jet moving inside an AGB wind- using different parameters (physical conditions and position angles of the jet). In agreement with our model and observations, these objects contain a dense torus seeing pole-on and a bipolar jet escaping thorough the poles. Then we propose to classify this kind of objects as spectroscopic bipolar nebulae, although they have been classified morphologically as compact, round, or elliptical nebulae or with "close collimated lobes".

We present ALMA Cycle 1 observations of the HH46/47 molecular outflow using combined 12m array and ACA observations. The improved angular resolution and sensitivity of our multi-line maps reveal structures that help us study the entrainment process in much more detail and allow us to obtain more precise estimates of outflow properties than previous observations. We use 13CO(1-0) and C18O(1-0) emission to correct for the 12CO(1-0) optical depth to accurately estimate the outflow mass, momentum and kinetic energy. This correction increases the estimates of the mass, momentum and kinetic energy by factors of about 9, 5 and 2, respectively, with respect to estimates assuming optically thin emission. The new 13CO and C18O data also allow us to trace denser and slower outflow material than that traced by the 12CO maps, and they reveal an outflow cavity wall at very low velocities (as low as 0.2km/s with respect to the cores central velocity). Adding with the slower material traced only by 13CO and C18O, there is another factor of 3 increase in the mass estimate and 50% increase in the momentum estimate. The estimated outflow properties indicate that the outflow is capable of dispersing the parent core within the typical lifetime of the embedded phase of a low-mass protostar, and that it is responsible for a core-to-star efficiency of 1/4 to 1/3. We find that the outflow cavity wall is composed of multiple shells associated with a series of jet bow-shock events. Within about 3000AU of the protostar the 13CO and C18O emission trace a circumstellar envelope with both rotation and infall motions, which we compare with a simple analytic model. The CS(2-1) emission reveals tentative evidence of a slowly-moving rotating outflow, which we suggest is entrained not only poloidally but also toroidally by a disk wind that is launched from relatively large radii from the source.

We present an analysis of four epochs of H$\alpha$ and [S II] $\lambda\lambda$ 6716/6731 HST images of HH 1. For determining proper motions we explore a new method based on analysis of spatially degraded images obtained convolving the images with wavelet functions of chosen widths. With this procedure we are able to generate maps of proper motion velocities along and across the outflow axis, as well as (angularly integrated) proper motion velocity distributions. From the four available epochs, we find the time evolution of the velocities, intensities and spatial distribution of the line emission. We find that over the last two decades HH 1 shows a clear acceleration. Also, the H$\alpha$ and [S II] intensities have first dropped, and then recovered in the more recent (2014) images. Finally, we show a comparison between the two available HST epochs of [O III] $\lambda$ 5007 (1994 and 2014), in which we see a clear drop in the value of the [O III]/H$\alpha$ ratio.

Close-in hot Jupiters are exposed to a tremendous photon flux that ionizes the neutral escaping material from the planet leaving an observable imprint that makes them an interesting laboratory for testing theoretical models. In this work we present 3D hydrodynamic simulations with radiation transfer calculations of a close-in exoplanet in a blow-off state. We calculate the Ly-$\alpha$ absorption and compare it with observations of HD 209458b an previous simplified model results.Our results show that the hydrodynamic interaction together with a proper calculation of the photoionization proccess are able to reproduce the main features of the observed Ly-$\alpha$ absorption, in particular at the blue-shifted wings of the line. We found that the ionizing stellar flux produce an almost linear effect on the amount of absorption in the wake. Varying the planetary mass loss rate and the radiation flux, we were able to reproduce the $10\%$ absorption observed at $-100~\mathrm{km~s^{-1}}$.

StarBench is a project focused on benchmarking and validating different star-formation and stellar feedback codes. In this first StarBench paper we perform a comparison study of the D-type expansion of an HII region. The aim of this work is to understand the differences observed between the twelve participating numerical codes against the various analytical expressions examining the D-type phase of HII region expansion. To do this, we propose two well-defined tests which are tackled by 1D and 3D grid- and SPH- based codes. The first test examines the `early phase' D-type scenario during which the mechanical pressure driving the expansion is significantly larger than the thermal pressure of the neutral medium. The second test examines the `late phase' D-type scenario during which the system relaxes to pressure equilibrium with the external medium. Although they are mutually in excellent agreement, all twelve participating codes follow a modified expansion law that deviates significantly from the classical Spitzer solution in both scenarios. We present a semi-empirical formula combining the two different solutions appropriate to both early and late phases that agrees with high-resolution simulations to $\lesssim2\%$. This formula provides a much better benchmark solution for code validation than the Spitzer solution. The present comparison has validated the participating codes and through this project we provide a dataset for calibrating the treatment of ionizing radiation hydrodynamics codes.

We present narrow-band optical and near-IR imaging and optical long-slit spectroscopic observations of Hu1-2, a Galactic planetary nebula (PN) with a pair of [N II]-bright, fast-moving (> 340 km/s) bipolar knots. Intermediate-dispersion spectra are used to derive physical conditions and abundances across the nebula, and high-dispersion spectra to study the spatio-kinematical structure. Generally Hu1-2 has high He/H (~0.14) and N/O ratios (~0.9), typical of Type I PNe. On the other hand, its abundances of O, Ne, S, and Ar are low as compared with the average abundances of Galactic bulge and disc PNe. The position-velocity maps can be generally described as an hour-glass shaped nebula with bipolar expansion, although the morphology and kinematics of the innermost regions cannot be satisfactorily explained with a simple, tilted equatorial torus. The spatio-kinematical study confines the inclination angle of its major axis to be within 10 degrees of the plane of sky. As in the irradiated bow-shocks of IC4634 and NGC7009, there is a clear stratification in the emission peaks of [O III], H_alpha, and [N II] in the northwest (NW) knot of Hu1-2. Fast collimated outflows in PNe exhibit higher excitation than other low-ionization structures. This is particularly the case for the bipolar knots of Hu1-2, with He II emission levels above those of collimated outflows in other Galactic PNe. The excitation of the knots in Hu1-2 is consistent with the combined effects of shocks and UV radiation from the central star. The mechanical energy and luminosity of the knots are similar to those observed in the PNe known to harbor a post-common envelope (post-CE) close binary central star.

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.

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=4\to 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.5\to 10\times 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.

Dwarf Galaxies are the most common objects in the Universe and are believed to contain large amounts of dark matter. There are mainly three morphologic types of dwarf galaxies: dwarf ellipticals, dwarf spheroidals and dwarf irregulars. Dwarf irregular galaxies are particularly interesting in dwarf galaxy evolution, since dwarf spheroidal predecessors could have been very similar to them. Therefore, a mechanism linked to gas-loss in dwarf irregulars should be observed, i.e. ram pressure stripping. In this paper, we study the interaction between the ISM of a dwarf galaxy, and a flowing IGM. We derive the weak-shock, plasmon solution corresponding to the balance between the post-bow shock pressure and the pressure of the stratified ISM (which we assume follows the fixed stratification of a gravitationally dominant dark matter halo). We compare our model with previously published numerical simulations and with the observed shape of the HI cloud around the Ho II and Pegasus dwarf irregular galaxies. We show that such a comparison provides a straightforward way for estimating the Mach number of the impinging flow.

We propose a new Planetary Nebula Luminosity Function (PNLF) that includes two populations in the distribution. Our PNLF is a direct extension of the canonical function proposed by Jacoby et al. (1987), in order to avoid problems related with the histogram construction, it is cast in terms of cumulative functions. We are interested in recovering the shape of the faint part of the PNLF in a consistent manner, for galaxies with and without a dip in their PN luminosity functions. The parameters for the two mode PNLF are obtained with a genetic algorithm, which obtains a best fit to the PNLF varying all of the parameters simultaneously in a broad parameter space. We explore a sample of 9 galaxies with various Hubble types and construct their PNLF. All of the irregular galaxies, except one, are found to be consistent with a two-mode population, while the situation is less clear for ellipticals and spirals.For the case of NGC\, 6822, we show that the two-mode PNLF is consistent with previous studies of the star formation history within that galaxy. Our results support two episodes of star formation, in which the latter is significantly stronger.

We report deep long-slit emission-line spectra, the line flux ratios, and Doppler profile shapes of various bright optical lines. The low-ionization lines (primarily [N I], [O I], [S II], and [N II]) originate in shocked knots, as reported by many previous observers. Dust-scattered lines of higher ionization are seen throughout the lobes but do not peak in the knots. Our analysis of these line profiles and the readily discernible stellar continuum shows that (1) the central star is an active symbiotic (whose spectrum resembles the central stars of highly bipolar and young PNe such as M2-9 and Hen2-437) whose compact companion shows a WC8-type spectrum, (2) extended nebular lines of [O III] and He I originate in the heavily obscured nuclear H II region, and (3) the Balmer lines observed throughout the lobes are dominated by reflected H{\alpha} emission from the symbiotic star. Comparing our line ratios with those observed historically shows that (1) the [O III]/H\b{eta} and He I/H\b{eta} ratios have been steadily rising by large amounts throughout the nebula, (2) the H{\alpha}/H\b{eta} ratio is steadily decreasing while H{\gamma}/H\b{eta} remains nearly constant, and (3) the low-ionization line ratios formed in the shocked knots have been in decline in different ways at various locations. We show that the first two of these results might be expected if the symbiotic central star has been active and if its bright H{\alpha} line has faded significantly in the past 20 years.

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 observer's 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.

We have used multiple mid-infrared observations at 4.5 micron obtained with the Infrared Array Camera, of the compact (~1.4 arcmin) young stellar bipolar outflow Cep E to measure the proper motion of its brightest condensations. The images span a period of ~6 yr and have been reprocessed to achieve a higher angular resolution (~0.8 arcsec) than their normal beam (2 arcsec). We found that for a distance of 730 pc, the tangential velocities of the North and South outflow lobes are 62+/-29 and 94+/-6 km/s respectively, and moving away from the central source roughly along the major axis of the flow. A simple 3D hydrodynamical simulation of the H2 gas in a precessing outflow supports this idea. Observations and model confirm that the molecular Hydrogen gas, traced by the pure rotational transitions, moves at highly supersonic velocities without being dissociated. This suggests either a very efficient mechanism to reform H2 molecules along these shocks or the presence of some other mechanism (e.g. strong magnetic field) that shields the H2 gas.

This work aims at improving the current understanding of the interaction between H ii regions and turbulent molecular clouds. We propose a new method to determine the age of a large sample of OB associations by investigating the development of their associated H ii regions in the surrounding turbulent medium. Using analytical solutions, one-dimensional (1D), and three-dimensional (3D) simulations, we constrained the expansion of the ionized bubble depending on the turbulent level of the parent molecular cloud. A grid of 1D simulations was then computed in order to build isochrone curves for H ii regions in a pressure-size diagram. This grid of models allowed to date large sample of OB associations and was used on the H ii Region Discovery Survey (HRDS). Analytical solutions and numerical simulations showed that the expansion of H ii regions is slowed down by the turbulence up to the point where the pressure of the ionized gas is in a quasi-equilibrium with the turbulent ram pressure. Based on this result, we built a grid of 1D models of the expansion of H ii regions in a profile based on Larson laws. The 3D turbulence is taken into account by an effective 1D temperature profile. The ages estimated by the isochrones of this grid agree well with literature values of well-known regions such as Rosette, RCW 36, RCW 79, and M16. We thus propose that this method can be used to give ages of young OB associations through the Galaxy such as the HRDS survey and also in nearby extra-galactic sources.

We present 3D, adaptive mesh refinement simulations of G2, a cloud of gas moving in a highly eccentric orbit towards the galactic center. We assume that G2 originates from a stellar wind interacting with the environment of the Sgr A* black hole. The stellar wind forms a cometary bubble which becomes increasingly elongated as the star approaches periastron. A few months after periastron passage, streams of material begin to accrete on the central black hole with accretion rates $\dot{M} \sim 10^{-8}$ M$_\odot$ yr$^{-1}$. Predicted Br$\gamma$ emission maps and position-velocity diagrams show an elongated emission resembling recent observations of G2. A large increase in luminosity is predicted by the emission coming from the shocked wind region during periastron passage. The observations, showing a constant Br$\gamma$ luminosity, remain puzzling, and are explained here assuming that the emission is dominated by the free-wind region. The observed Br$\gamma$ luminosity ($\sim 8 \times 10^{30}$ erg s$^{-1}$) is reproduced by a model with a $v_w=50$ km s$^{-1}$ wind velocity and a $10^{-7}$ M$_\odot$ yr$^{-1}$ mass loss rate if the emission comes from the shocked wind. A faster and less dense wind reproduces the Br$\gamma$ luminosity if the emission comes from the inner, free wind region. The extended cometary wind bubble, largely destroyed by the tidal interaction with the black hole, reforms a few years after periastron passage. As a result, the Br$\gamma$ emission is more compact after periastron passage.

We present 3D hydrodynamical simulations of an isotropic fast wind interacting with a previously ejected toroidally-shaped slow wind in order to model both the observed morphology and the kinematics of the planetary nebula (PN) NGC 6302. This source, also known as the Butterfly nebula, presents one of the most complex morphologies ever observed in PNe. From our numerical simulations, we have obtained an intensity map for the H$\alpha$ emission to make a comparison with the Hubble Space Telescope (HST) observations of this object. We have also carried out a proper motion (PM) study from our numerical results, in order to compare with previous observational studies. We have found that the two interacting stellar wind model reproduces well the morphology of NGC 6302, and while the PM in the models are similar to the observations, our results suggest that an acceleration mechanism is needed to explain the Hubble-type expansion found in HST observations.

In this review we focus on the role jets and outflows play in the star and planet formation process. Our essential question can be posed as follows: are jets/outflows merely an epiphenomenon associated with star formation or do they play an important role in mediating the physics of assembling stars both individually and globally? We address this question by reviewing the current state of observations and their key points of contact with theory. Our review of jet/outflow phenomena is organized into three length-scale domains: Source and Disk Scales ($0.1-10^2$ au) where the connection with protostellar and disk evolution theories is paramount; Envelope Scales ($10^2-10^5$ au) where the chemistry and propagation shed further light on the jet launching process, its variability and its impact on the infalling envelope; Parent Cloud Scales ($10^5-10^6$ au) where global momentum injection into cluster/cloud environments become relevant. Issues of feedback are of particular importance on the smallest scales where planet formation regions in a disk may be impacted by the presence of disk winds, irradiation by jet shocks or shielding by the winds. Feedback on envelope scales may determine the final stellar mass (core-to-star efficiency) and envelope dissipation. Feedback also plays an important role on the larger scales with outflows contributing to turbulent support within clusters including alteration of cluster star formation efficiencies (feedback on larger scales currently appears unlikely). A particularly novel dimension of our review is that we consider results on jet dynamics from the emerging field of High Energy Density Laboratory Astrophysics (HEDLA). HEDLA is now providing direct insights into the 3-D dynamics of fully magnetized, hypersonic, radiative outflows.

We present a 3D hydrodynamic study of the effects that different stellar wind conditions and planetary wind structures have on the calculated Ly-$\alpha$ absorptions produced during the transit of HD 209458b. Considering a range of stellar wind speeds $\sim$[350-800] km s$^{-1}$, coronal temperature $\sim$[3-7] $\times10^{6}$ K and two values of the polytropic index $\Gamma$ $\sim$[1.01-1.13], while keeping fixed the stellar mass loss rate, we found a that a $\dot M_p$ range between $\sim$[3-5] $\times 10^{10}$g s$^{-1}$ give account for the observational absorption in Ly-$\alpha$ measured for the planetary system. Also, several models with anisotropic evaporation profiles for the planetary escaping atmosphere were carried out, showing that both, the escape through polar regions and through the night side yields larger absorptions than an isotropic planetary wind.

We present multi-epoch Very Long Baseline Array (VLBA) H$_2$O maser observations toward the massive young stellar objects (YSOs) VLA 2 and VLA 3 in the star-forming region AFGL 2591. Through these observations, we have extended the study of the evolution of the masers towards these objects up to a time span of $\sim$ 10 yrs, measuring their radial velocities and proper motions. The H$_2$O masers in VLA 3, the most massive YSO in AFGL 2591 ($\sim$ 30--40~M$_{\odot}$), are grouped within projected distances of $\lesssim$ 40 mas ($\lesssim$ 130 AU) from VLA 3. In contrast to other H$_2$O masers in AFGL 2591, the masers associated with VLA 3 are significantly blueshifted (up to $\sim$ 30 km s$^{-1}$) with respect to the velocity of the ambient molecular cloud. We find that the H$_2$O maser cluster as a whole, has moved westwards of VLA~3 between the 2001 and 2009 observations, with a proper motion of $\sim$ 1.2 mas yr$^{-1}$ ($\sim$ 20 km s$^{-1}$). We conclude that these masers are tracing blueshifted outflowing material, shock excited at the inner parts of a cavity seen previously in ammonia molecular lines and infrared images, and proposed to be evacuated by the outflow associated with the massive VLA 3 source. The masers in the region of VLA 2 are located at projected distances of $\sim$ 0.7$''$ ($\sim$ 2300 AU) north from this source, with their kinematics suggesting that they are excited by a YSO other than VLA 2. This driving source has not yet been identified.

We present a model of a "thermal" (i.e., a hot bubble) rising within an exponentially stratified region of the ISM. This model includes terms representing the ram pressure braking and the entrainment of environmental gas into the thermal. We then calibrate the free parameters associated with these two terms through a comparison with 3D numerical simulations of a rising bubble. Finally, we apply our "thermal" model to the case of a hot bubble produced by a SN within the stratified ISM of the Galactic disk.

The morphology, kinematics and entrainment mechanism of the HH 46/47 molecular outflow were studied using new ALMA Cycle 0 observations. Results show that the blue and red lobes are strikingly different. We argue that these differences are partly due to contrasting ambient densities that result in different wind components having a distinct effect on the entrained gas in each lobe. A 29-point mosaic, covering the two lobes at an angular resolution of about 3", detected outflow emission at much higher velocities than previous observations, resulting in significantly higher estimates of the outflow momentum and kinetic energy than previous studies of this source, using the CO(1-0) line. The morphology and the kinematics of the gas in the blue lobe are consistent with models of outflow entrainment by a wide-angle wind, and a simple model describes the observed structures in the position-velocity diagram and the velocity-integrated intensity maps. The red lobe exhibits a more complex structure, and there is evidence that this lobe is entrained by a wide-angle wind and a collimated episodic wind. Three major clumps along the outflow axis show velocity distribution consistent with prompt entrainment by different bow shocks formed by periodic mass ejection episodes which take place every few hundred years. Position-velocity cuts perpendicular to the outflow cavity show gradients where the velocity increases towards the outflow axis, inconsistent with outflow rotation. Additionally, we find evidence for the existence of a small outflow driven by a binary companion.

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.

In this paper we analyze multi-epoch Very Long Baseline Interferometry (VLBI) water maser observations carried out with the Very Long Baseline Array (VLBA) toward the high-mass star-forming region AFGL 2591. We detected maser emission associated with the radio continuum sources VLA 2 and VLA 3. In addition, a water maser cluster, VLA 3-N, was detected ~ 0.5" north of VLA 3. We concentrate the discussion of this paper on the spatio-kinematical distribution of the water masers towards VLA 3-N. The water maser emission toward the region VLA 3-N shows two bow shock-like structures, Northern and Southern, separated from each other by ~ 100 mas (~ 330 AU). The spatial distribution and kinematics of the water masers in this cluster have persisted over a time span of seven years. The Northern bow shock has a somewhat irregular morphology, while the Southern one has a remarkably smooth morphology. We measured the proper motions of 33 water maser features, which have an average proper motion velocity of ~ 1.3 mas/yr (~ 20 km/s). The morphology and the proper motions of this cluster of water masers show systematic expanding motions that could imply one or two different centers of star formation activity. We made a detailed model for the Southern structure, proposing two different kinematic models to explain the 3-dimensional spatio-kinematical distribution of the water masers: (1) a static central source driving the two bow-shock structures; (2) two independent driving sources, one of them exciting the Northern bow-shock structure, and the other one, a young runaway star moving in the local molecular medium exciting and molding the remarkably smoother Southern bow-shock structure. Future observations will be necessary to discriminate between the two scenarios, in particular by identifying the still unseen driving source(s).

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 (Cant\'o 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.

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.

The classic HII region M17 is one of the best studied across the electromagnetic spectrum. We present sensitive, high angular resolution observations made with the Jansky Very Large Array (JVLA) at 4.96 and 8.46 GHz that reveal the presence of 38 compact radio sources, in addition to the well known hypercompact cometary HII region M17 UC1. For this last source we find that its spectral index of value $\sim$1 is due to a gradient in opacity across its face. Of the 38 compact radio sources detected, 19 have stellar counterparts detected in the infrared, optical, or X-rays. Finally, we discuss the nature of the radio emission from the massive binary system CEN 1a and 1b, concluding that both are most probably non-thermal emitters, although the first is strongly time variable and the second is steady.

DG Tau B is a Class I young stellar source that drives the asymmetric HH 159 bipolar jet. At optical wavelengths it is obscured by circumstellar optically-thick material. Using VLA and JVLA observations, we determine for the first time the proper motions of this source and find them to be consistent, within error, with those of the nearby young star DG Tau. We also discuss an ejection event that is evident in the 1994 VLA data. As the optical and molecular outflows, this ejection traced in the radio continuum is markedly asymmetric and was detected only to the NW of the star. We propose that this knot, no longer detectable in the radio, could be observed in future optical images of DG Tau B. The positions of the VLA source and of a nearby infrared object are not coincident and we suggest that the VLA source traces the exciting object, while the infrared source could be a reflection lobe.

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.

Context: HH 158, the jet from the young star DG Tau, is one of the few sources of its type where jet knots have been detected at optical and X-ray wavelengths. Aims: To search, using Very Large Array observations of this source, radio knots and if detected, compare them with the optical and X-ray knots. To model the emission from the radio knots. Methods: We analyzed archive data and also obtained new Very Large Array observations of this source, as well as an optical image, to measure the present position of the knots. We also modeled the radio emission from the knots in terms of shocks in a jet with intrinsically time-dependent ejection velocities. Results: We detected radio knots in the 1996.98 and 2009.62 VLA data. These radio knots are,within error, coincident with optical knots. We also modeled satisfactorily the observed radio flux densities as shock features from a jet with intrinsic variability. All the observed radio, optical, and X-ray knot positions can be intepreted as four successive knots, ejected with a period of 4.80 years and traveling away from the source with a velocity of 198 km s$^{-1}$ in the plane of the sky. Conclusions: The radio and optical knots are spatially correlated and our model can explain the observed radio flux densities. However, the X-ray knots do not appear to have optical or radio counterparts and their nature remains poorly understood.

We explore a ballistic orbit model to infer the gravitational drag force on an accreting point mass M, such as a black hole, moving at a hypersonic velocity v_{0} through a gaseous environment of density \rho_{0}. The streamlines blend in the flow past the body and transfer momentum to it. The total drag force acting on the body, including the nonlinear contribution of those streamlines with small impact parameter that bend significantly and pass through a shock, can be calculated by imposing conservation of momentum. In this fully analytic approach, the ambiguity in the definition of the lower cut-off distance $r_{\rm min}$ in calculations of the effect of dynamical friction is removed. It turns out that $r_{\rm min}=\sqrt{e}GM/2v_{0}^{2}$. Using spherical surfaces of control of different sizes, we carry out a successful comparison between the predicted drag force and the one obtained from a high resolution, axisymmetric, isothermal flow simulation. We demonstrate that ballistic models are reasonably successful in accounting for both the accretion rate and the gravitational drag.

We present a "thin shell" model of the interaction of a biconical outflow and a streaming environment (aligned with the direction of the flow), as well as numerical (axisymmetric) simulations of such an interaction. A similar situation, although in a more complex setup, takes place at the head of the cometary structure of Mira. Thus, for most of the numerical simulations we explore parameters consistent with the observed bipolar outflow from Mira B. For these parameters, the interaction is non-radiative, so that a rather broad jet/streaming environment interaction region is formed. In spite of this, a reasonable agreement between the thin-shell analytic model and the numerical simulations is obtained.

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.

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 results from models of galactic winds driven by energy injected by nuclear starbursts. The total energy of the starburst is provided by young central stellar clusters and parts of the galactic interstellar medium are pushed out as part of the galactic wind (in some cases the galactic wind contains an important part of the metals produced in the new generation of stars). We have performed adiabatic and radiative 3D N-Body/Smooth Particle Hydrodynamics simulations of galactic winds using the GADGET-2 code. The numerical models cover a wide range of starburst (from $\sim10^2$ to $\sim10^7$ M$_\odot$) and galactic gas masses (from $\sim6\times10^6$ to $\sim10^{11}$ M$_\odot$). The concentrated central starburst regions are an efficient engine for producing of the mass and metal loss in galaxies, and also for driving the metal redistribution in the galaxies.

We present a new Spitzer IRAC images of the HH 34 outflow. These are the first images that detect both the knots along the southern jet and the northern counterjet (the counterjet knots were only detected previously in a long slit spectrum). This result removes the problem of the apparent coexistence of a large scale symmetry (at distances of up to approx.1 pc) and a complete lack of symmetry close to the source (at distances of approx. 1.e+17 cm) for this outflow. We present a quantitative evaluation of the newly found symmetry between the HH 34 jet and counterjet, and show that the observed degree of symmetry implies that the jet production region has a characteristic size <2.8 AU. This is the strongest constraint yet derived for the size of the region in which HH jets are produced.

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.

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.

We present five epochs of VLBI water maser observations around the massive protostar Cepheus A HW2 with 0.4 mas (0.3 AU) resolution. The main goal of these observations was to follow the evolution of the remarkable water maser linear/arcuate structures found in earlier VLBI observations. Comparing the data of our new epochs of observation with those observed five years before, we find that at "large" scales of > 1" (700 AU) the main regions of maser emission persist, implying that both the surrounding medium and the exciting sources of the masers have been relatively stable during that time span. However, at smaller scales of < 0.1" (70 AU) we see large changes in the maser structures, particularly in the expanding arcuate structures R4 and R5. R4 traces a nearly elliptical patchy ring of ~ 70 mas size (50 AU) with expanding motions of ~ 5 mas/yr (15 km/s). This structure is probably driven by the wind of a still unidentified YSO located at the centre of the ring (~ 0.18" south of HW2). On the other hand, the R5 expanding bubble structure (driven by the wind of a previously identified YSO located ~ 0.6" south of HW2) is currently dissipating in the circumstellar medium and losing its previous degree of symmetry, indicating a very short-lived event. In addition, our results reveal, at scales of ~ 1" (700 AU), the simultaneous presence of a relatively slow (~ 10-70 km/s) wide-angle outflow (opening angle of ~ 102 deg, traced by the masers, and the fast (~ 500~km/s) highly collimated radio jet associated with HW2 (opening angle of ~ 18 deg, previously observed with the VLA. This simultaneous presence of a wide-angle outflow and a highly collimated jet associated with a massive protostar is similar to what is found in some low-mass YSOs. The implications of these results in the study of the formation of high-mass stars are discussed.

We present two Very Large Array observations of the pre-main-sequence star PV Cephei, taken with a separation of 10.5 years. These data show that the proper motions of this star are $\mu_\alpha \cos \delta = +10.9 \pm 3.0$ mas yr$^{-1}$; $\mu_\delta = +0.2 \pm 1.8$ mas yr$^{-1}$, very similar to those -previously known- of HD 200775, the B2Ve star that dominates the illumination of the nearby reflection nebula NGC 7023. This result suggests that PV Cephei is not a rapidly moving run-away star as suggested by previous studies. The large velocity of PV Cephei had been inferred from the systematic eastward displacement of the bisectors of successive pairs of Herbig Haro knots along its flow. These systematic shifts might instead result from an intrinsic dissymmetry in the ejection mechanisms, or from an asymmetric distribution of the circumstellar material.

The physical parameters of Herbig-Haro jets are usually determined from emission line ratios, obtained from spectroscopy or narrow band imaging, assuming that the emitting region is homogeneous along the line of sight. Under the more general hypothesis of axisymmetry, we apply tomographic reconstruction techniques to the analysis of Herbig-Haro jets. We use data of the HH30 jet taken by Hartigan & Morse (2007) with the Hubble space telescope using the slitless spectroscopy technique. Using a non-parametric Tikhonov regularization technique, we determine the volumetric emission line intensities of the [SII]6716,6731, [OI]6300 and [NII]6583 forbidden emission lines. From our tomographic analysis of the corresponding line ratios, we produce "three-dimensional" images of the physical parameters. The reconstructed density, temperature and ionization fraction present much steeper profiles than those inferred using the assumption of homogeneity. Our technique reveals that the reconstructed jet is much more collimated than the observed one close to the source (a width ~ 5 AU vs. ~ 20 AU at a distance of 10 AU from the star), while they have similar widths at larger distances. In addition, our results show a much more fragmented and irregular jet structure than the classical analysis, suggesting that the the ejection history of the jet from the star-disk system has a shorter timescale component (~ some months) superimposed on a longer, previously observed timescale (of a few years). Finally, we discuss the possible application of the same technique to other stellar jets and planetary nebulae.

We present an analytic model for the interaction between planetary atmospheres and stellar winds from main sequence M stars, with the purpose of obtaining a quick test-model that estimates the timescale for total atmospheric mass loss due to this interaction. Planets in the habitable zone of M dwarfs may be tidally locked and may have weak magnetic fields, because of this we consider the extreme case of planets with no magnetic field. The model gives the planetary atmosphere mass loss rate as a function of the stellar wind and planetary properties (mass, atmospheric pressure and orbital distance) and an entrainment efficiency coefficient $\alpha$. We use a mixing layer model to explore two different cases: a time-independent stellar mass loss and a stellar mass loss rate that decreases with time. For both cases we consider planetary masses within the range of $1\to10$ M$_{\oplus}$ and atmospheric pressures with values of 1, 5 and 10 atm. For the time dependent case, planets without magnetic field in the habitable zone of M dwarfs with initial stellar mass losses of $\leq \dot{M}_{w} < 10^{-11}$ M$_{\odot}$ yr$^{-1}$, may retain their atmospheres for at least 1 Gyr. This case may be applied to early spectral type M dwarfs (earlier than M5). Studies have shown that late type M dwarfs (later than M5) may be active for long periods of time ($\geq 4$Gyr), and because of that our model with constant stellar mass loss rate may be more accurate. For these stars most planets may have lost their atmospheres in 1 Gyr or less because most of the late type M dwarfs are expected to be active. We emphasize that our model only considers planets without magnetic fields. Clearly we must expect a higher resistance to atmospheric erosion if we include the presence of a magnetic field.

We have done a statistical analysis of Very Long Baseline Array (VLBA) data of water masers in the star-forming regions (SFRs) Cepheus A and W75 N, using correlation functions to study the spatial clustering and Doppler-velocity distribution of these masers. Two-point spatial correlation functions show a characteristic scale size for clusters of water maser spots < or ~1 AU, similar to the values found in other SFRs. This suggests that the scale for water maser excitation tends to be < or ~1 AU. Velocity correlation functions show power-law dependences with indices that can be explained by regular velocity fields, such as expansion and/or rotation. These velocity fields are similar to those indicated by the water maser proper-motion measurements; therefore, the velocity correlation functions appear to reveal the organized motion of water maser spots on scales larger than 1 AU.

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 here new results of two-dimensional hydrodynamical simulations of the eruptive events of the 1840s (the great) and the 1890s (the minor) eruptions suffered by the massive star $\eta$ Car. The two bipolar nebulae commonly known as the Homunculus and the little Homunculus were formed from the interaction of these eruptive events with the underlying stellar wind. As in previous work (Gonzalez et al. 2004a, 2004b), we assume here an interacting, nonspherical multiple-phase wind scenario to explain the shape and the kinematics of both Homunculi, but adopt a more realistic parametrization of the phases of the wind. During the 1890s eruptive event, the outflow speed {\it decreased} for a short period of time. This fact suggests that the little Homunculus is formed when the eruption ends, from the impact of the post-outburst $\eta$ Car wind (that follows the 1890s event) with the eruptive flow (rather than by the collision of the eruptive flow with the pre-outburst wind, as claimed in previous models; Gonzalez et al. 2004a, 2004b). Our simulations reproduce quite well the shape and the observed expansion speed of the large Homunculus. The little Homunculus (which is embedded within the large Homunculus) becomes Rayleigh-Taylor unstable and develop filamentary structures that resembles the spatial features observed in the polar caps. In addition, we find that the interior cavity between the two Homunculi is partially filled by material that is expelled during the decades following the great eruption. This result may be connected with the observed double-shell structure in the polar lobes of the $\eta$ Car nebula. Finally, as in previous work, we find the formation of tenuous, equatorial, high-speed features that seem to be related to the observed equatorial skirt of $\eta$ Car.

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 are proposing to conduct a multicolor, synoptic infrared (IR) imaging survey of the Northern sky with a new, dedicated 6.5-meter telescope at San Pedro M\'artir (SPM) Observatory. This initiative is being developed in partnership with astronomy institutions in Mexico and the University of California. The 4-year, dedicated survey, planned to begin in 2017, will reach more than 100 times deeper than 2MASS. The Synoptic All-Sky Infrared (SASIR) Survey will reveal the missing sample of faint red dwarf stars in the local solar neighborhood, and the unprecedented sensitivity over such a wide field will result in the discovery of thousands of z ~ 7 quasars (and reaching to z > 10), allowing detailed study (in concert with JWST and Giant Segmented Mirror Telescopes) of the timing and the origin(s) of reionization. As a time-domain survey, SASIR will reveal the dynamic infrared universe, opening new phase space for discovery. Synoptic observations of over 10^6 supernovae and variable stars will provide better distance measures than optical studies alone. SASIR also provides significant synergy with other major Astro2010 facilities, improving the overall scientific return of community investments. Compared to optical-only measurements, IR colors vastly improve photometric redshifts to z ~ 4, enhancing dark energy and dark matter surveys based on weak lensing and baryon oscillations. The wide field and ToO capabilities will enable a connection of the gravitational wave and neutrino universe - with events otherwise poorly localized on the sky - to transient electromagnetic phenomena.

The radiation from newly born stars photoevaporates their parental neutral cloud, leading to the formation of dense clumps that will eventually form stars. We present 3D simulations of the interaction of a neutral cloud with an external ionising radiation field, and compute the angular momenta of these collapsing clumps. The angular momenta of these collapsing clumps show that they have preferential orient mostly perpendicular to the direction of the incident ionising photon field. Therefore, the axes of the jet systems that will be eventually ejected (from the star + accretion disk systems that will form) will be oriented approximately perpendicular to the direction to the photoionising source.

The well-established correlations between the mass of massive black holes (BHs) in the nuclei of most studied galaxies and various global properties of their hosting galaxy lend support to the idea that dwarf galaxies and globular clusters could also host a BH in their centers. Direct kinematic detection of BHs in dwarf spheroidal (dSph) galaxies are seriously hindered by the small number of stars inside the gravitational influence region of the BH. The aim of this Letter is to establish an upper dynamical limit on the mass of the putative BH in the Ursa Minor (UMi) dSph galaxy. We present direct N-body simulations of the tidal disruption of the dynamical fossil observed in UMi, with and without a massive BH. We find that the observed substructure is incompatible with the presence of a massive BH of (2-3)x10^4 Msun within the core of UMi. These limits are consistent with the extrapolation of the M_{BH}-sigma relation to the M_{BH}<10^6 Msun regime. We also show that the BH may be off-center with respect to the center of symmetry of the whole galaxy.

Intermediate and high velocity HI clouds rain onto the plane of our Galaxy. They are observed at heights of between 500 and 1500 pc, falling onto the Galactic plane at velocities from 50 to 140 km s$^{-1}$. To explain the origin of these clouds, we present a galactic fountain model, driven by the wind from a super stellar cluster (SSC). We solve the equations for a steady, radiative de Laval nozzle flow. We consider two effects not considered previously in astrophysical nozzle flow models: cooling functions for different metallicities, and the direct action of the galactic gravitational field on the gas flowing along the nozzle. For an adiabatic nozzle flow, the gravity acting directly on the gas within the nozzle "stalls" the nozzle flow for initial wind velocities lower than the escape velocity from the Galaxy. For the same wind velocity, a radiative nozzle flow stalls at lower altitudes above the galactic plane. We find that SSC winds with velocities of $v_w=500 - 800$ km s$^{-1}$ produce nozzles stall at heights of $x_m=1 - 15$ kpc. The stalled nozzle flow then rains back onto the galactic plane at velocities in the range observed in intermediate and high velocity HI clouds. We study a nozzle flow driven by a wind from a SSC close to the Galactic centre. We find that for velocities within the range expected for a SSC wind, we can produce nozzle flows that stall above the galactic plane. These stalled flows produce cool, infalling clouds with velocities similar to those of intermediate and high velocity HI clouds.

Context. SuperclusterA in the extragalactic HII region NGC2363 is remarkable for the hypersonic gas seen as faint extended broad emission lines with a full width zero intensity of 7000km/s. Aims. We explore the possibility that the observed broad profiles are the result of the interaction of a high velocity cluster wind with dense photoionized clumps. Methods. The geometry considered is that of near static photoionized condensations at the surface of which turbulent mixing layers arise as a result of the interaction with the hot wind. The approximative treatment of turbulence is carried out using the mixing length approach of Canto & Raga. The code mappings Ic is used to derive the mean quantities describing the flow and to compute the line emissivities within the turbulent layers. The velocity projection in three dimensions of the line sources is carried out analytically. Results. A fast entraining wind of up to ~4300km/s appears to be required to reproduce the faint wings of the broad H-alpha and [O III] profiles. A slower wind of 3500km/s, however, can still reproduce the bulk of the broad component and does provide a better fit than an ad hoc Gaussian profile. Conclusions. Radial acceleration in 3D (away from supercluster A) of the emission gas provides a reasonable first order fit to the broad line component. No broad component is predicted for the [N II] and [S II] lines, as observed. The wind velocity required is uncomfortably high and alternative processes that would provide comparable constant acceleration of the emission gas up to 4000km/s might have to be considered.

It is generally believed that O stars, confined near the galactic midplane, are somehow able to photoionize a significant fraction of what is termed the "diffuse ionized gas" (DIG) of spiral galaxies, which can extend up to 1-2 kpc above the galactic midplane. The heating of the DIG remains poorly understood, however, as simple photoionization models do not reproduce the observed line ratio correlations well or the DIG temperature. We present turbulent mixing layer models in which warm photoionized condensations are immersed in a hot supersonic wind. Turbulent dissipation and mixing generate an intermediate region where the gas is accelerated, heated and mixed. The emission spectrum of such layers are compared with observations of Rand (ApJ 462, 712) of the DIG in the edge-on spiral NGC2363. We generate two sequence of models that fit the line ratio correlations between [SII]/H-alpha, [OI]/H-alpha, [NII]/[SII] and [OIII]/H-beta reasonably well. In one sequence of models the hot wind velocity increases while in the other the ionization parameter and layer opacity increases. Despite the success of the mixing layer models, the overall efficiency in reprocessing the stellar UV is much too low, much less than 1%, which compels us to reject the TML model in its present form.

We review here the effects of supernovae (SNe) explosions on the environment of star-forming galaxies. Randomly distributed, clustered SNe explosions cause the formation of hot superbubbles that drive either galactic fountains or supersonic winds out of the galactic disk. In a galactic fountain, the ejected gas is re-captured by the gravitational potential and falls back onto the disk. From 3D non-equilibrium radiative cooling hydrodynamical simulations of these fountains, we find that they may reach altitudes smaller than 5 kpc in the halo and hence explain the formation of the so-called intermediate-velocity-clouds (IVCs). On the other hand, the high-velocity-clouds (HVCs) that are observed at higher altitudes (of up to 12 kpc) require another mechanism to explain their production. We argue that they could be formed either by the capture of gas from the intergalactic medium and/or by the action of magnetic fields that are carried out to the halo with the gas in the fountains. Due to angular momentum losses (of 10-15%) to the halo, we find that the fountain material falls back to smaller radii and is not largely spread over the galactic disk, as previously expected. This result is consistent with the metal distribution derived from recent chemical models of the galaxy. We also find that after about 150 Myr, the gas circulation between the halo and the disk in the fountains reaches a steady state regime (abridged).

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.

The dynamical effects of magnetic fields in models of radiative, Herbig-Haro (HH) jets have been studied in a number of papers. For example, magnetized, radiative jets from variable sources have been studied with axisymmetric and 3D numerical simulations. In this paper, we present an analytic model describing the effect of a toroidal magnetic field on the internal working surfaces that result from a variability in the ejection velocity. We find that for parameters appropriate for HH jets the forces associated with the magnetic field dominate over the gas pressure force within the working surfaces. Depending on the ram pressure radial cross section of the jet, the magnetic field can produce a strong axial pinch, or, alternatively, a broadening of the internal working surfaces. We check the validity of the analytic model with axisymmetric numerical simulations of variable, magnetized jets.

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.

We explore the dynamics of a ``cluster wind'' flow in the regime in which the shocks resulting from the interaction of winds from nearby stars are radiative. We first show that for a cluster with T Tauri stars and/or Herbig Ae/Be stars, the wind interactions are indeed likely to be radiative. We then compute a set of four, three dimensional, radiative simulations of a cluster of 75 young stars, exploring the effects of varying the wind parameters and the density of the initial ISM that permeates the volume of the cluster. These simulations show that the ISM is compressed by the action of the winds into a structure of dense knots and filaments. The structures that are produced resemble in a qualitative way the observations of the IRAS 18511+0146 of Vig et al. (2007).

CONTEXT: In the interpretation of stellar jet observations, the physical parameters are usually determined from emission line ratios, obtained from spectroscopic observations or using the information contained in narrow band images. The basic hypothesis in the interpretation of the observations is that the emitting region is homogeneous along the line of sight. Actually, stellar jets are in general not homogeneous, and therefore line of sight convolution effects may lead to the main uncertainty in the determination of the physical parameters. AIMS: This paper is aimed at showing the systematic errors introduced when assuming an homogeneous medium, and studying the effect of an inhomogeneous medium on plasma diagnostics for the case of a stellar jet. In addition, we explore how to reconstruct the volumetric physical parameters of the jet (i. e., with dependence both across and along the line of sight). METHODS: We use standard techniques to determine the physical parameters, i. e., the electron density, temperature and hydrogen ionisation fraction across the jet, and a multi-Gaussian method to invert the Abel transform and determine the reconstructed physical structure. RESULTS: When assuming an homogeneous medium the physical parameters, integrated along the line of sight, do not represent the average of the true values, and do not have a clear physical interpretation. We show that when some information is available on the emissivity profile across the jet, it is then possible to obtain appropriate derivations of the electron density, temperature and ionisation fraction.

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.

With its remarkable double-S shape, IC 4634 is an archetype of point-symmetric planetary nebulae (PN). In this paper, we present a detailed study of this PN using archival HST WFPC2 and ground-based narrow-band images to investigate its morphology, and long-slit spectroscopic observations to determine its kinematics and to derive its physical conditions and excitation. The data reveal new structural components, including a distant string of knots distributed along an arc-like feature 40"-60" from the center of the nebula, a skin of enhanced [O III]/H-alpha ratio enveloping the inner shell and the double-S feature, and a triple-shell structure. The spatio-kinematical study also finds an equatorial component of the main nebula that is kinematically independent from the bright inner S-shaped arc. We have investigated in detail the bow shock-like features in IC 4634 and found that their morphological, kinematical and emission properties are consistent with the interaction of a collimated outflow with surrounding material. Indeed, the morphology and kinematics of some of these features can be interpreted using a 3D numerical simulation of a collimated outflow precessing at a moderate, time-dependent velocity. Apparently, IC 4634 has experienced several episodes of point-symmetric ejections oriented at different directions with the outer S-shaped feature being related to an earlier point-symmetric ejection and the outermost arc-like string of knots being the relic of an even much earlier point-symmetric ejection. There is tantalizing evidence that the action of these collimated outflows has also taken part in the shaping of the innermost shell and inner S-shaped arc of IC 4634.

We here explore the effects of the SN explosions into the environment of star-forming galaxies like the Milky Way. Successive randomly distributed and clustered SNe explosions cause the formation of hot superbubbles that drive either fountains or galactic winds above the galactic disk, depending on the amount and concentration of energy that is injected by the SNe. In a galactic fountain, the ejected gas is re-captured by the gravitational potential and falls back onto the disk. From 3D nonequilibrium radiative cooling hydrodynamical simulations of these fountains, we find that they may reach altitudes up to about 5 kpc in the halo and thus allow for the formation of the so called intermediate-velocity-clouds (IVCs) which are often observed in the halos of disk galaxies. The high-velocity-clouds that are also observed but at higher altitudes (of up to 12 kpc) require another mechanism to explain their production. We argue that they could be formed either by the capture of gas from the intergalactic medium and/or by the action of magnetic fields that are carried to the halo with the gas in the fountains. Due to angular momentum losses to the halo, we find that the fountain material falls back to smaller radii and is not largely spread over the galactic disk. Instead, the SNe ejecta fall nearby the region where the fountain was produced, a result which is consistent with recent chemical models of the galaxy. The fall back material leads to the formation of new generations of molecular clouds and to supersonic turbulence feedback in the disk.

In the present study, considering the physical conditions that are relevant interactions between supernova remnants (SNRs) and dense molecular clouds for triggering star formation we have built a diagram of SNR radius versus cloud density in which the constraints above delineate a shaded zone where star formation is allowed. We have also performed fully 3-D radiatively cooling numerical simulations of the impact between SNRs and clouds under different initial conditions in order to follow the initial steps of these interactions. We determine the conditions that may lead either to cloud collapse and star formation or to complete cloud destruction and find that the numerical results are consistent with those of the SNR-cloud density diagram. Finally, we have applied the results above to the $\beta-$ Pictoris stellar association which is composed of low mass Post-T Tauri stars with an age of 11 Myr. It has been recently suggested that its formation could have been triggered by the shock wave produced by a SN explosion localized at a distance of about 62 pc that may have occurred either in the Lower Centaurus Crux (LCC) or in the Upper Centaurus Lupus (UCL) which are both nearby older subgroups of that association (Ortega and co-workers).Using the results of the analysis above we have shown that the suggested origin for the young association at the proposed distance is plausible only for a very restricted range of initial conditions for the parent molecular cloud, i.e., a cloud with a radius of the order of 10 pc, a density of the order of 10$-$20 cm$^{-3}$, and a temperature of the order of 10$-$100 K }.

The bipolar morphology of the planetary nebula (PN) K 3-35 observed in radio-continuum images was modelled with 3D hydrodynamic simulations with the adaptive grid code yguazu-a. We find that the observed morphology of this PN can be reproduced considering a precessing jet evolving in a dense AGB circumstellar medium, given by a mass loss rate \dot{M}_{csm}=5x10^{-5}M_{\odot}/yr and a terminal velocity v_{w}=10 km/s. Synthetic thermal radio-continuum maps were generated from numerical results for several frequencies. Comparing the maps and the total fluxes obtained from the simulations with the observational results, we find that a model of precessing dense jets, where each jet injects material into the surrounding CSM at a rate \dot{M}_j=2.8x10^{-4} {M_{\odot}/yr (equivalent to a density of 8x10^{4} {cm}^{-3}, a velocity of 1500 km/s, a precession period of 100 yr, and a semi-aperture precession angle of 20 degrees agrees well with the observations.

We present analytic and numerical models of the `cluster wind' resulting from the multiple interactions of the winds ejected by the stars of a dense cluster of massive stars. We consider the case in which the distribution of stars (i.e., the number of stars per unit volume) within the cluster is spherically symmetric, has a power-law radial dependence, and drops discontinuously to zero at the outer radius of the cluster. We carry out comparisons between an analytic model (in which the stars are considered in terms of a spatially continuous injection of mass and energy) and 3D gasdynamic simulations (in which we include 100 stars with identical winds, located in 3D space by statistically sampling the stellar distribution function). From the analytic model, we find that for stellar distributions with steep enough radial dependencies the cluster wind flow develops a very high central density and a non-zero central velocity, and for steeper dependencies it becomes fully supersonic throughout the volume of the cluster (these properties are partially reproduced by the 3D numerical simulations). Therefore, the wind solutions obtained for stratified clusters can differ dramatically from the case of a homogeneous stellar distribution (which produces a cluster wind with zero central velocity, and a fully subsonic flow within the cluster radius). Finally, from our numerical simulations we compute predictions of X-ray emission maps and luminosities, which can be directly compared with observations of cluster wind flows.

We have studied the maser emission from a thin, planar, gaseous ring in Keplerian rotation around a central mass observed edge-on. The absorption coefficient within the ring is assumed to follow a power law dependence with the distance from the central mass as, k=k0r^{-q}. We have calculated position-velocity diagrams for the most intense maser features, for different values of the exponent q. We have found that, depending on the value of q, these diagrams can be qualitatively different. The most intense maser emission at a given velocity can either come mainly from regions close to the inner or outer edges of the amplifying ring or from the line perpendicular to the line of sight and passing through the central mass (as is commonly assumed). Particularly, when q>1 the position-velocity diagram is qualitatively similar to the one observed for the water maser emission in the nucleus of the galaxy NGC 4258. In the context of this simple model, we conclude that in this object the absorption coefficient depends on the radius of the amplifying ring as a decreasing function, in order to have significant emission coming from the inner edge of the ring.

We present [SII] images of the HH 30 and HL/XZ Tau region obtained at two epochs, as well as long-slit optical spectroscopy of the HH 30 jet. We measured proper motions of about 100-300 km/s for the HH 30 jet and counterjet, and of about 120 km/s for the HL Tau jet. Inclination angles with respect to the plane of the sky are 0-40 deg for the HH 30 jet and 60 deg for the HL Tau jet. Comparison with previous observations suggests that most of the jet knots consist of persisting structures. Also, we corroborate that the HH 30-N knots correspond to the head of the HH 30 jet. The overall HH 30 jet structure can be well described by a wiggling ballistic jet, arising either by the orbital motion of the jet source around a primary or by precession of the jet axis because of the tidal effects of a companion. In the first scenario, the orbital period would be 53 yr and the total mass 0.25-2 solar masses. In the precession scenario, the mass of the jet source would be 0.1-1 solar masses, the orbital period <1 yr, and the mass of the companion less than a few times 0.01 solar masses, thus being a substellar object or a giant exoplanet. In both scenarios a binary system with a separation <18 AU (<0.13 arcsec) is required. Since the radius of the flared disk observed with the HST is about 250 AU, we conclude that this disk appears to be circumbinary rather than circumstellar, suggesting that the search for the collimating agent of the HH 30 jet should be carried out at much smaller scales.

We present high spatial resolution optical integral field spectroscopy of a collimated Herbig-Haro jet viewed nearly edge-on. Maps of the line emission, velocity centroid, and velocity dispersion were generated for the H$\alpha$ and [S II] emission features from the inner collimated jet and exciting source region of the HH 34 outflow. The kinematic structure of the jet shows several maxima and minima in both velocity centroid value and velocity dispersion along the jet axis. Perpendicular to the flow direction the velocity decreases outward from the axis to the limb of the jet, but the velocity dispersion increases. Maps of the electron density structure were derived from the line ratio of [S II] 6731/6716 emission. We have found that the jet exhibits a pronounced ``striped'' pattern in electron density; the high $n_e$ regions are at the leading side of each of the emission knots in the collimated jet, and low $n_e$ regions in the down-flow direction. On average, the measured electron density decreases outward from the inner regions of the jet, but the highest $n_e$ found in the outflow is spatially offset from the nominal position of the exciting star. The results of our high spatial resolution optical integral field spectroscopy show very good agreement with the kinematics and electron density structure predicted by the existing internal working surface models of the HH~34 outflow.

In the present study, considering the physical conditions that are relevant in interactions between supernova remnants (SNRs) and dense molecular clouds for triggering star formation we have built a diagram of SNR radius versus cloud density in which the constraints above delineate a shaded zone where star formation is allowed. We have also performed fully 3-D radiatively cooling numerical simulations of the impact between SNRs and clouds under different initial conditions in order to follow the initial steps of these interactions. We determine the conditions that may lead either to cloud collapse and star formation or to complete cloud destruction and find that the numerical results are consistent with those of the SNR-cloud density diagram. Finally, we have applied the results above to the $\beta-$Pictoris stellar association which is composed of low mass Post-T Tauri stars with an age of 11 Myr. It has been recently suggested that its formation could have been triggered by the shock wave produced by a SN explosion localized at a distance of about 62 pc that may have occurred either in the Lower Centaurus Crux (LCC) or in the Upper Centaurus Lupus (UCL) which are both nearby older subgroups of that association (Ortega and co-workers). Using the results of the analysis above we have shown that the suggested origin for the young association at the proposed distance is plausible only for a very restricted range of initial conditions for the parent molecular cloud, i.e., a cloud with a radius of the order of 10 pc and density of the order of 20 cm$^{-3}$ and a temperature of the order of 50$-$100 K.

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 present new optical STIS HST spectroscopic observations of the jets of the proto-planetary nebula Hen 3-1475. The excitation conditions of the knots of Hen 3-1475 are derived from the observed optical spectra, confirming that the knots are shock excited. The shocked spectra are qualitatively reproduced by simple ``3/2''D bow shock models. We present a set of bow shock models devoted to planetary nebulae, and discuss the effects of the pre-ionization conditions, the bow shock velocity, the bow shock shape and the chemical abundances on the predicted spectra. To explore the reliability of the ``3/2''D bow shock models, we also compare the observed spectra of other three proto-planetary nebulae (M 1-92, M 2-56 and CRL 618) to the predicted spectra.

We study the H$\alpha$ emission from jets using two-dimensional axisymmetrical simulations. We compare the emission obtained from hydrodynamic (HD) simulations with that obtained from magnetohydrodynamics (MHD) simulations. The magnetic field is supposed to be present in the jet only, and with a toroidal configuration. The simulations have time-dependent ejection velocities and different intensities for the initial magnetic field. The results show an increase in the H$\alpha$ emission along the jet for the magnetized cases with respect to the HD case. The increase in the emission is due to a better collimation of the jet in the MHD case, and to a small increase in the shock velocity. These results could have important implications for the interpretation of the observations of jets from young stellar objects.

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 hydrodynamical calculations carried out with the 3D yguazu-a code of a precessing jet model, which interacts with a plane parallel wind. This scenario describes an extragalactic jet, in which the jet source is in motion with respect to the surrounding intra-cluster medium. From the numerical results, synthetic emission maps and spectra in X-ray band were obtained. We compare these predictions with observations of the radio jets emanating from the radio-galaxy 4C 26.42 (in the Abell 1795 galaxy cluster). We find that the general morphology of the radio jets can be described by a point-symmetric precessing jet system interacting with a plane parallel wind (i.e., the intra-cluster medium flowing past the galaxy). We also find that our synthetic X-ray emission maps reproduce the observed large scale structures (with sizes of the order of tens of kpc).

In this work we are interested to study the by-products of SNR-clouds in a starburst (SB) system. These interactions can have an important role in the recycling of matter from the clouds to the ISM and vice-versa. In the present work, we have focused our attention on the global effects of the interactions between clouds and SN shock waves in the ISM of SB environments, and performed 3-D radiative cooling hydrodynamical simulations with the adaptive YGUAZU grid code. We have also considered the effects of the photo-evaporation due to the presence of a high number of UV photons from hot stars and supernovae (SNe). The results have shown that, in the presence of radiative cooling, instead of an efficient gas mixing with the diffuse ISM, the interactions cause the fragmentation of the clouds into smaller ones. The results have also revealed that the SNR-clouds interactions are less efficient at producing substantial mass loss from the clouds to the diffuse ISM than mechanisms such as the photo-evaporation caused by the UV flux from the hot stars.

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.

The propagation of cosmological ionization fronts (I-fronts) during reionization is strongly influenced by small-scale structure. Here we summarize our recent attemps to understand the effect of this small-scale structure. We present high resolution cosmological N-body simulations at high-z (z>6) which resolve a wide range of halo mass, from minihalos to clusters of large, rare halos. We also study how minihalos affect I-fronts, through simulations of minihalo photoevaporation. Furthermore, we modify the I-front propagation equations to account for evolving small-scale structure, and incorporate these results into a semi-analytical reionization model. When intergalactic medium clumping and minihalo clustering around sources are included, small-scale structure affects reionization by slowing it down and extending it in time. This helps to explain observations by the Wilkinson Microwave Anisotropy Probe, which imply an early and extended reionization epoch. We also study how source clustering affects the evolution and size of H II regions, finding, in agreement with simulations, that H II regions usually expand, rarely shrinking. Hence, "relic H II regions'' are an exception, rather than the rule. When the suppression of small-mass sources in already-ionized regions by Jeans-mass filtering is accounted for, H II regions are smaller, delaying overlap. We also present a new numerical method for radiative transfer which is fast, efficient, and easily coupled to hydrodynamics and N-body codes, along with sample tests and applications.

We present numerical simulations of the interaction between a collimated, bipolar ``pulse'' ejected from a star and a continuous wind ejected from a stellar companion. We explore the characteristics of the predicted H$\alpha$ intensity maps by varying selected input parameters. We find that the asymmetry (in size and strength) between the two lobes of the proto-planetary nebula OH 231.8+4.2 is reproduced in this scenario if the wind ejected by the companion star has a strong latitude dependence.

We present new results on the kinematics of the jet HH 110. New proper motion measurements have been calculated from [SII] CCD images obtained with a time baseline of nearly fifteen years. HH 110 proper motions show a strong asymmetry with respect to the outflow axis, with a general trend of pointing towards the west of the axis direction. Spatial velocities have been obtained by combining the proper motions and radial velocities from Fabry-Perot data. Velocities decrease by a factor ~3 over a distance of ~10$^{18}$ cm, much shorter than the distances expected for the braking caused by the jet/environment interaction. Our results show evidence of an anomalously strong interaction between the outflow and the surrounding environment, and are compatible with the scenario in which HH 110 emerges from a deflection in a jet/cloud collision.

First evidences of IR excess and disk mass accretion (strong H$\alpha$ emission) around brown dwarfs seem to indicate the existence of circumstellar disks around these sub-stellar objects. Nothing is known at the present time about outflows which potentially might be launched from brown dwarfs, although jets are typically associated with the accretion in standard T Tauri star disks. In this paper we calculate the H$\alpha$ emission of internal working surfaces produced by a radiative jet in a neutral and in a photoionized environment as a function of the jet parameters (the ejection velocity $v_j$, shock velocity $v_s$, mass loss rate ${\dot M}$ and radius $r_j$ of the jet) and we provide estimates of the H$\alpha$ luminosity for the parameters of ``standard'' Herbig-Haro (HH) jets from T Tauri stars and for the parameters expected for jets from BDs. Interestingly, we find that while the mass loss rates associated with jets from BDs are found to be two orders of magnitude lower than the mass loss rates associated with ``standard'' HH jets (from T Tauri stars), their velocities are likely to be similar. Based on our calculations, we discuss the conditions in which jets from BDs can be detected and we conclude that the H$\alpha$ luminosities of internal working surfaces of jets from BDs in a photoionized environment should have only one order of magnitude lower than the H$\alpha$ luminosities of T Tauri jets in a neutral environment.

In this paper, we present two-dimensional gas dynamic simulations of the formation and evolution of the eta-Car bipolar outflows. Adopting the interacting nonspherical winds model, we have carried out high-resolution numerical simulations, which include explicitly computed time-dependent radiative cooling, for different possible scenarios of the colliding winds. In our simulations, we consider different degrees of non-spherical symmetry for the pre-outburst wind and the great eruption of the 1840s presented by the eta-Car wind. From these models, we obtain important differences in the shape and kinematical properties of the Homunculus structure. In particular, we find an appropriate combination of the wind parameters (that control the degree of non-spherical symmetry) and obtain numerical experiments that best match both the observed morphology and the expansion velocity of the eta-Car bipolar shell. In addition, our numerical simulations show the formation of a bipolar nebula embedded within the Homunculus (the little Homunculus) developed from a secondary eruptive event suffered by the star in the 1890s, and also the development of tenuous, high velocity ejections in the equatorial region that result from the impact of the eruptive wind of the 1840s with the pre-outburst wind and that could explain some of the high speed features observed in the equatorial ejecta. The models were, however, unable to produce equatorial ejections associated to the second eruptive event.

The weak, R-type ionization fronts (I-fronts) which swept across the intergalactic medium (IGM) during the reionization of the universe often found their paths blocked by cosmological minihaloes. When this happened, the neutral gas which filled each minihalo was photoevaporated; as the I-front burned its way through the halo, decelerating from R-type to D-type, all the gas was blown back into the IGM as an ionized, supersonic wind. In a previous paper (Shapiro, Iliev and Raga 2004), we described this process and presented our results of the first simulations of it by numerical gas dynamics with radiation transport in detail. For illustration we focused on the particular case of a 10^7 solar masses minihalo overrun at z=9 by an intergalactic I-front caused by a distant source of ionizing radiation, for different source spectra (either stellar from massive Pop. II or III stars, or QSO-like) and a flux level typical of that expected during reionization. In a LambdaCDM universe, minihaloes formed in abundance before and during reionization and, thus, their photoevaporation is an important, possibly dominant, feature of reionization, which slowed it down and wasted ionizing photons. We have now performed a larger set of high-resolution simulations to determine and quantify the dependence of minihalo photoevaporation times and photon consumption rates on halo mass, redshift, ionizing flux level and spectrum. We find that the average number of ionizing photons each minihalo atom absorbs during its photoevaporation is typically in the range 2-10. For the collapsed fraction in minihaloes expected during reionization, this can add about 1 photon per total atom to the requirements for completing reionization, potentially doubling the minimum number of photons required to reionize the universe.

Molecular clouds in the interstellar medium suffer gravitational instabilities that lead to the formation of one or multiple stars. A recently formed star inside a cold cloud communicates its gravitational force to the surrounding environment and soon an accretion flow falling into the star develops. After their formation, all stars soon eject a wind of gas that interacts with the external accretion flow. This interaction produces a shock wave that evolves with time. The work presented in this article formulates a simple prescription for the evolution of this interaction. With the aid of this model we construct a few radio continuum maps of the source.

We present 2D axisymmetric and 3D numerical simulations of the proto-planetary nebula Hen 3-1475, which is characterized by a remarkably highly collimated optical jet, formed by a string of shock-excited knots along the axis of the nebula. It has recently been suggested that the kinematical and morphological properties of the Hen 3-1475 jet could be the result of an ejection variability of the central source (Riera et al. 2003). The observations suggest a periodic variability of the ejection velocity superimposed on a smoothly increasing ejection velocity ramp. From our numerical simulations, we have obtained intensity maps (for different optical emission lines) and position-velocity diagrams, in order to make a direct comparison with the HST observations of this object. Our numerical study allows us to conclude that a model of a precessing jet with a time-dependent ejection velocity, which is propagating into an ISM previously perturbed by an AGB wind, can succesfully explain both the morphological and the kinematical characteristics of this proto-planetary nebula.

We here present two-dimensional, time-dependent radiatively cooling hydrodynamical simulations of the large and little Homunculus nebulae around $\eta$ Carinae. We employ an alternative scenario to previous interacting stellar wind models which is supported by both theoretical and observational evidence, where a non-spherical outburst wind (with a latitudinal velocity dependence that matches the observations of the large Homunculus), which is expelled for 20 years, interacts with a pre-eruptive slow wind also with a toroidal density distribution, but with a much smaller equator-to-polar density contrast than that assumed in previous models. A second eruptive wind with spherical shape is ejected about 50 years after the first outburst, and causes the development of the little internal nebula. We find that, as a result of an appropriate combination of the parameters that control the degree of asymmetry of the interacting winds, we are able to produce not only the structure and kinematics of both Homunculus, but also the high-velocity equatorial ejecta. These arise from the impact between the non-spherical outburst and the pre-outburst winds in the equatorial plane.

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.

We present the first gas dynamical simulations of the photoevaporation of cosmological minihalos overtaken by the ionization fronts which swept through the IGM during reionization in a LCDM universe, including the effects of radiative transfer. We demonstrate the phenomenon of I-front trapping inside minihalos, in which the weak, R-type fronts which traveled supersonically across the IGM decelerated when they encountered the dense, neutral gas inside minihalos, becoming D-type I-fronts, preceded by shock waves. For a minihalo with virial temperature T_vir < 10^4 K, the I-front gradually burned its way through the minihalo which trapped it, removing all of its baryonic gas by causing a supersonic, evaporative wind to blow backwards into the IGM, away from the exposed layers of minihalo gas just behind the advancing I-front. Such hitherto neglected feedback effects were widespread during reionization. N-body simulations and analytical estimates of halo formation suggest that sub-kpc minihalos such as these, with T_vir < 10^4 K, were so common as to cover the sky around larger-mass source halos and possibly dominate the absorption of ionizing photons. This means that previous estimates of the number of ionizing photons per H atom required to complete reionization which neglected this effect may be too low. Regardless of their effect on the progress of reionization, however, the minihalos were so abundant that random lines of sight thru the high-z universe should encounter many of them, which suggests that it may be possible to observe the processes described here in the absorption spectra of distant sources.

We present an analytical model for the cross section of a turbulent, radiative jet or wake. This model is appropriate for modeling HH jets, or "wakes" left behind by "astrophysical bullets". Even though the model is very simple, it has the benign property of only having four free parameters (the outer radius of the beam, the axial velocity, the velocity at the edge of the beam, and the turbulent velocity width), which can be derived by fitting the radial velocity and line width cross sections of an observed outflow. We illustrate how to do such fits using previously published spectroscopic data of the HH 110 jet.

We have obtained a Halpha position-velocity cube from Fabry-Perot interferometric observations of the HH 110 flow. We analyze the results in terms of anisotropic wavelet transforms, from which we derive the spatial distribution of the knots as well as their characteristic sizes (along and across the outflow axis). We then study the spatial behaviour of the line width and the central radial velocity. The results are interpreted in terms of a simple ``mean flow+turbulent eddy'' jet/wake model. We find that most of the observed kinematics appear to be a direct result of the mean flow, on which are superposed low amplitude (35 km/s) turbulent velocities.

We present the first gas dynamical simulations of the photoevaporation of cosmological minihalos overtaken by the ionization fronts which swept through the IGM during reionization in a LCDM universe, including the effects of radiative transfer. We demonstrate the phenomenon of I-front trapping inside minihalos, in which the weak, R-type fronts which traveled supersonically across the IGM decelerated when they encountered the dense, neutral gas inside minihalos, becoming D-type I-fronts, preceded by shock waves. For a minihalo with virial temperature T_vir < 10,000 K, the I-front gradually burned its way through the minihalo which trapped it, removing all of its baryonic gas by causing a supersonic, evaporative wind to blow backwards into the IGM, away from the exposed layers of minihalo gas just behind the advancing I-front. Such hitherto neglected feedback effects were widespread during reionization. N-body simulations and analytical estimates of halo formation suggest that sub-kpc minihalos such as these, with T_vir < 10,000 K, were so common as to dominate the absorption of ionizing photons. This means that previous estimates of the number of ionizing photons per H atom required to complete reionization which neglected this effect may be too low. Regardless of their effect on the progress of reionization, however, the minihalos were so abundant that random lines of sight thru the high-z universe should encounter many of them, which suggests that it may be possible to observe the processes described here in the absorption spectra of distant sources.

We present long-slit spectroscopic observations of the HH 110 jet obtained with the 4.2~m William Herschel Telescope. We have obtained for the first time, spectra for slit positions along and across the jet axis (at the position of knots B, C, I, J and P) to search for the observational signatures of entrainment and turbulence by studying the kinematics and the excitation structure. We find that the HH 110 flow accelerates from a velocity of 35 km/s in knot A up to 110 km/s in knot P. We find some systematic trends for the variation of the emission line ratios along the jet. No clear trends for the variation of the radial velocity are seen across the width of the jet beam. The cross sections of the jet show complex radial velocity and line emission structures which differ quite strongly from each other.

The magnitude of the temperature fluctuations (t^2) required to explain the observed inconsistencies between metallicities inferred from recombination lines and from forbidden lines cannot be attained by steady-state equilibrium photoionization models. If on the other hand the nebular ionizing source was variable, the temperature fluctuations t^2 would be significantly larger. We investigate the time-dependent response of the nebular ionization and temperature structure when photoionized by a periodically varying source. We study how the asymptotic mean value, <t^2>, behaves as a function of the period or amplitude of the source variability. We find that the temperature fluctuations occur only in the outer section of the nebula, close to the ionization front, within a zone corresponding to 8-20% of the ionized layer's thickness. We conclude that the amplitude of the exciting star variations required to achieve a <t^2> = 0.025 (as in the Orion nebula) is unacceptably large. Source variability is therefore not a viable mechanism to explain the observed values of t^2. We reach a similar conclusion from studies of the temporal variability resulting from intermittent shadows behind opaque condensations. We find that photoionized nebulae are on average less massive but somewhat hotter in the case of cyclicly variable ionizing sources.

Numerical simulations of the interaction between supernova ejecta and a stellar wind are presented. We follow the temporal evolution of the shock fronts that are formed through such an interaction and determine the velocities, temperatures and densities. We model the X-ray emission from the SNR-stellar wind collision region and we compare it with recent results from X-ray observations carried out with the Chandra satellite of the SMC supernova remnant SNR 0057-7226 which could be interacting with the wind of the Wolf-Rayet system HD 5980. The simulations predict the presence of shell-like regions of enhanced X-ray emission which are consistent with the presence of X-ray emitting arcs in the Chandra image. Also the observed X-ray luminosity is comparable to the X-ray luminosities we obtain from the simulations for a supernova with an initial energy in the (1-5)E50 erg range.

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

The recent detection of X-ray emission from HH 2 and HH 154 with the Chandra and XMM-Newton satellites (respectively) have opened up an interesting, new observational possibility in the field of Herbig-Haro objects. In order to be able to plan further X-ray observations of other HH objects, it is now of interest to be able to estimate their X-ray luminosities in order to choose which objects to observe. This paper describes a simple, analytic model for predicting the X-ray luminosity of a bow shock from the parameters of the flow (i.e., the size of the bow shock, its velocity, and the pre-shock density). The accuracy of the analytic model is analyzed through a comparison with the predictions obtained from axisymmetric, gasdynamic simulations of the leading working surface of an HH jet. We find that our analytic model reproduces the observed X-ray luminosities of HH 2 and HH 154, and we propose that HH~80/81 is a good candidate for future observations with Chandra.

We present 3D, gasdynamic simulations of jet/cloud collisions, with the purpose of modelling the HH 270/110 system. From the models, we obtain predictions of H$\alpha$ and H$_2$ 1-0 s(1) emission line maps, which qualitatively reproduce some of the main features of the corresponding observations of HH 110. We find that the model that better reproduces the observed structures corresponds to a jet that was deflected at the surface of the cloud $\sim 1000$ yr ago, but is now boring a tunnel directly into the cloud. This model removes the apparent contradiction between the jet/cloud collision model and the lack of detection of molecular emission in the crossing region of the HH 270 and HH 110 axes.

The kinematics and physical conditions of the bipolar planetary nebula M 2-48 are analysed from high and low dispersion long-slit spectra. Previous CCD narrow-band optical observations have suggested that this nebula is mainly formed by a pair of symmetric bow-shocks, an off-center semi-circular shell, and an internal bipolar structure. The bipolar outflow has a complex structure, characterised by a series of shocked regions located between the bright core and the polar tips. There is an apparent kinematic discontinuity between the bright bipolar core and the outer regions. The fragmented ring around the bright bipolar region presents a low expansion velocity and could be associated to ejection in the AGB-PN transition phase, although its nature remains unclear. The chemical abundances of the central region are derived, showing that M 2-48 is a Type I planetary nebula (PN).

We present a long-slit spectrum of the Herbig-Haro object HH 111, obtained with SITS on the Hubble Space Telescope. This spectrum has a spectral resolution of approx. 50 km/s, so it gives a good picture of the kinematical properties of the observed object at very high angular resolution. We find that some of the knots along the jet are associated with sudden drops in the radial velocity (modulus), confirming that the emission from the knots is formed in shocks. We interprt the observed position-velocity diagrams in terms of a model of a jet from a variable source, and we attempt to carry out a reconstruction of the ejection variability necessary for reproducing the observed kinematical structure.

We have surveyed the inner 5' of the Orion Nebula by means of Halpha and [NII] Fabry-Perot imaging spectroscopy to present a kinematical study of the Herbig-Haro objects in the nebula. The objects studied in this work are HH 202, 203, 204, 529, 269 and other associated features. For HH 202 we find new features that, because of their high velocities (up to 100 km/seg) indicate the presence of an outflow that probably is a HH flow not catalogued previously. HH 202 could be only a part of this larger outflow. Large internal motions are found in the fainter regions of HH 203-204, as well as evidence of transverse density gradients that could account for the asymmetry in the brightness distribution of HH 204. We report for the first time a high blueshifted velocity (-118 km/seg) associated with HH 204, and show that the apex of HH 204 is indeed the zone of maximum velocity, in agreement with bow shock models. We also studied the radial velocity field of HH 269 finding features associated with the HH object. From our studies, we find kinematic evidence that suggests that HH 203-204 and HH 202 are part of a big (aprox. 0.55 pc) bipolar HH outflow.

We explain the Ha emission of the cometary knots in the Helix Nebula (NGC 7293) with an analytical model that describes the emission of the head of the globules as a photoevaporated flow produced by the incident ionizing radiation of the central star.We compare these models with the Ha emission obtained from the HST archival images of the Helix Nebula. From a comparison of the Ha emission with the predictions of the analytical model we obtain a rate of ionizing photons from the central star of about 5e45 s^-1, which is consistent with estimates based on the total Hb flux of the nebula. We also model the tails of the cometary knots as a photoevaporated wind from a neutral shadow region produced by the diffuse ionizing photon field of the nebula. A comparison with the HST images allows us to obtain a direct determination of the value of the diffuse ionizing flux. We compare the ratio of diffuse to direct stellar flux as a function of radius inside an HII region with those obtained from the observational data through the analytical tail and head wind model. The agreement of this model with the values determined from the observations of the knots is excellent.

The first sources of ionizing radiation to condense out of the dark and neutral IGM sent ionization fronts sweeping outward through their surroundings, overtaking other condensed objects and photoevaporating them. This feedback of universal reionization on cosmic structure formation is demonstrated here by gas dynamical simulations, including radiative transfer, for a cosmological minihalo of dark matter and baryons exposed to an external source of ionizing radiation, either starlight or quasar light, just after the passage of the global ionization front created by the source.

The first sources of ionizing radiation to condense out of the dark and neutral IGM sent ionization fronts sweeping outward through their surroundings, overtaking other condensed objects and photoevaporating them. This feedback of universal reionization on cosmic structure formation is demonstrated here by gas dynamical simulations, including radiative transfer, for the case of a cosmological minihalo of dark matter and baryons exposed to an external source of ionizing starlight, just after the passage of the global ionization front created by the source.

We present images and spectra of the Cepheus E (Cep E) region at both optical and infrared wavelengths. Only the brightest region of the southern lobe of the Cep E outflow reveals optical emission, suggesting that the extinction close to the outflow source plays an important r\^ole in the observed difference between the optical and IR morphologies. Cep E is a unique object since it provides a link between the spectroscopic properties of the optical Herbig-Haro (HH) objects and those of deeply embedded outflows.

The first sources of ionizing radiation to condense out of the dark and neutral IGM sent ionization fronts sweeping outward through their surroundings, overtaking other condensed objects and photoevaporating them. This feedback effect of universal reionization on cosmic structure formation is demonstrated here for the case of a cosmological minihalo of dark matter and baryons exposed to an external source of ionizing radiation with a quasar-like spectrum, just after the passage of the global ionization front created by the source. We model the pre-ionization minihalo as a truncated, nonsingular isothermal sphere in hydrostatic equilibrium following its collapse out of the expanding background universe and virialization. Results are presented of the first gas dynamical simulations of this process, including radiative transfer. A sample of observational diagnostics is also presented, including the spatially-varying ionization levels of C, N, and O in the flow if a trace of heavy elements is present and the integrated column densities of H I, He I and II, and C IV thru the photoevaporating gas at different velocities which would be measured in absorption against a background source like that responsible for the ionization.

The postcollapse structure of objects which form by gravitational condensation out of the expanding cosmological background universe is a key element in the theory of galaxy formation. Towards this end, we have reconsidered the outcome of the nonlinear growth of a uniform, spherical density perturbation in an unperturbed background universe - the cosmological ``top-hat'' problem. We adopt the usual assumption that the collapse to infinite density at a finite time predicted by the top-hat solution is interrupted by a rapid virialization caused by the growth of small-scale inhomogeneities in the initial perturbation. We replace the standard description of the postcollapse object as a uniform sphere in virial equilibrium by a more self-consistent one as a truncated, nonsingular, isothermal sphere in virial and hydrostatic equilibrium, including for the first time a proper treatment of the finite-pressure boundary condition on the sphere. The results differ significantly from both the uniform sphere and the singular isothermal sphere approximations for the postcollapse objects. These results will have a significant effect on a wide range of applications of the Press-Schechter and other semi-analytical models to cosmology. The truncated isothermal sphere solution presented here predicts the virial temperature and integrated mass distribution of the X-ray clusters formed in the CDM model as found by detailed, 3D, numerical gas and N-body dynamical simulations remarkably well. This solution allows us to derive analytically the numerically-calibrated mass-temperature and radius-temperature scaling laws for X-ray clusters which were derived empirically by Evrard, Metzler and Navarro from simulation results for the CDM model. (Shortened)

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.

The first sources of ionizing radiation to condense out of the dark and neutral IGM sent ionization fronts sweeping outward through their surroundings, overtaking other primordial gas-clouds and photoevaporating them. Results are presented of the first gas dynamical simulations of this process, including radiative transfer, along with some observational diagnostics.

Energy released by a small fraction of the baryons in the universe, which condensed out of while the IGM was cold, dark, and neutral,reheated and reionized it, exposing gas clouds within it to the glare of ionizing radiation. The first gas dynamical simulations of the photoevaporation of an intergalactic cloud by a quasar, including radiative transfer, are presented, along with a few observational diagnostics.

We study the evolution of dense neutral clumps located in the outer parts of planetary nebulae. These clumps will be photo-ionized by the ionizing radiation from the central star and change their structure in the process. The main effect of the ionization process is the setting up of a photo-evaporation flow and a shock running through the clump. Once this shock has moved through the entire clump it starts to accelerate because of the `rocket effect'. This continues until the entire clump has been photo-ionized. We present an analytic model for the shock and accelerating phases and also the results of numerical simulations which include detailed microphysics. We find a good match between the analytic description and the numerical results and use the numerical results to produce some of the clump's observational characteristics at different phases of its evolution. We compare the results with the properties of the fast moving low ionization knots (ansae or FLIERs) seen in a number of planetary nebulae. We find that the models match many of the kinematic and emission properties of FLIERs.

We study the kinematical and brightness evolution of emission line clouds in the narrow line region (NLR) of Seyfert galaxies during the passage of a jet. We derive a critical density above which a cloud remains radiative after compression by the jet cocoon. The critical density depends mainly on the cocoon pressure. Super-critical clouds increase in emission line brightness, while sub-critical clouds generally are highly overheated reducing their luminosity below that of the inter-cloud medium. Due to the pressure stratification in the bow-shock of the jet, a cylindrical structure of nested shells develops around the jet. The most compact and brightest compressed clouds surround the cloud-free channel of the radio jet. To support our analytical model we present a numerical simulation of a supersonic jet propagating into a clumpy NLR. The position-velocity diagram of the simulated H_alpha emission shows total line widths of the order of 500 km/s with large-scale variations in the radial velocities of the clouds due to the stratified pressure in the bow-shock region of the jet. Most of the luminosity is concentrated in a few dense clouds surrounding the jet. These morphological and kinematic signatures are all found in the well observed NLR of NGC1068 and other Seyfert galaxies.

We present new narrow-band near-infrared images together with K band spectra of highly collimated bipolar jets close to the IRAS 05487+0255 source. The jets are located at 50" West of the Herbig-Haro 110 outflow. The jets are not visible at optical wavelengths, and therefore, do not fall into the `standard' Herbig-Haro object classification scheme. Nevertheless, they belong to an ever growing group of molecular hydrogen jets associated with YSOs which are optically undetected. The jets are very well collimated, with a length-to-width ratio of 10-20. The spectra of the jet and counter-jet in the K-band show a limited number of molecular hydrogen emission lines which makes it difficult to obtain an accurate excitation temperature. We estimate Tex = 1104+/-67 K and Tex = 920+/- 156 K for the red and blue jet components respectively. The radial velocities of the jet and counter-jet, based on the shift of the (1,0) S(1) 2.121 micron line, are -275+/- 50 km/s and 180+/- 50 km/s respectively, suggesting an angle of 30 to 45 degrees between the jet and the line of sight. The molecular hydrogen emission of the entire jet extends for at least 40" or 0.1 pc at the distance of Orion. If the flow velocity is comparable to that of the radial velocities, then the dynamical age of the system is quite short (about 500 yrs), consistent with a young jet arising from an embedded source. Entrainment in a turbulent mixing layer may explain this morphology and spectral character.

We present non-adiabatic hydrodynamic simulations of a supersonic light jet propagating into a fully ionized medium of uniform density on a scale representative of the narrow line region (NLR) in Seyfert galaxies with associated radio jets. In this regime the cooling distance of the swept up gas in the bowshock of the jet is of the same order as the transverse extent of the jet bowshock, as opposed to the more extreme regimes found for more powerful adiabatic large scale jets or the slow galactic jets which have been simulated previously. We calculate the emissivity for the H_alpha line and radio synchrotron emission. We find that the structure of the line emitting cold envelope of the jet cocoon is strongly dependent on the non-stationary dynamics of the jet head as it propagates through the ambient medium. We observe the formation of cloud-like high density regions which we associate with NLR clouds and filaments. We find that some of these clouds might be partially neutral and represent sites of jet induced star formation. The calculated H_alpha flux and the spectral line width are consistent with NLR observations. The simulation of the radio-optical emission with radiative cooling confirms the basic result of the geometric bowshock model developed by Taylor et al. (1989) that the start of noticeable optical line emission can be significantly offset from the hotspot of the radio emission. However, the time-dependent nature of the jet dynamics implies significant differences from their geometric bowshock model.

We have chosen a reduced set of 18 ionization rate equations (for ions of H, C, N, O, S and Ne), which allow us to obtain a moderately accurate estimate of the non-equilibrium radiative cooling function. We evaluate the accuracy of this approach by comparing our cooling function with previous calculations, computed with a more extended set of ions, for the case of gas that cools from a high temperature at constant density. We also compute steady, plane shock models, which we find to compare well with models calculated with much more detailed microphysics. Using our reduced set of rate equations, we present a simulation of a radiative bow shock formed by a supersonic, plane stream impinging on a rigid sphere. The parameters for the calculation are chosen as to approximately represent typical values found for Herbig-Haro objects, and to give a cooling distance to bow shock radius ratio of 1/10. This simulation is done with an adaptive grid code, which allows a reasonably high resolution (with greater than or equal to 25 points) of the standoff distance between the bow shock and the rigid obstacle. Contrary to the standard expectation, we find that the bow shock standoff distance in the stagnation region is considerably shorter than the cooling distance behind a plane-parallel shock. Also, the centrifugal pressure is found to play an important role in determining the structure of the recombin- ation region. This appears to partially invalidate previous bow shock models based on a ``quasi-1D'' approach, at least for the particular parameters chosen for the present simulation. Finally, we present tabulations of the cooling rates that have been used (for the different species), in order to facilitate the inclusion of this treatment in other gasdynamic codes.

Detailed comparisons are made between observations of scattered light upstream of the head of the HH~1 jet and predictions of simple scattering models. It is shown that, in order to unambiguously determine the velocity of the head of the jet (bow shock) with respect to the upstream dust, existing spectroscopic observations are insufficient and that spectropolarimetric observations of the scattered light are necessary. Such an independent measure of the bow shock velocity is important in order to test ``multiple outflow'' theories of Herbig-Haro jets. It is also shown 2that the scattering dust must have a very forward-throwing scattering phase function ($\langle\cos\theta\rangle\msim 0.7$) and slight evidence is found for a dust-gas ratio that is higher than average.