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.

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}}$.

After the death of a runaway massive star, its supernova shock wave interacts with the bow shocks produced by its defunct progenitor, and may lose energy, momentum, and its spherical symmetry before expanding into the local interstellar medium (ISM). We investigate whether the initial mass and space velocity of these progenitors can be associated with asymmetric supernova remnants. We run hydrodynamical models of supernovae exploding in the pre-shaped medium of moving Galactic core-collapse progenitors. We find that bow shocks that accumulate more than about 1.5 Mo generate asymmetric remnants. The shock wave first collides with these bow shocks 160-750 yr after the supernova, and the collision lasts until 830-4900 yr. The shock wave is then located 1.35-5 pc from the center of the explosion, and it expands freely into the ISM, whereas in the opposite direction it is channelled into the region of undisturbed wind material. This applies to an initially 20 Mo progenitor moving with velocity 20 km/s and to our initially 40 Mo progenitor. These remnants generate mixing of ISM gas, stellar wind and supernova ejecta that is particularly important upstream from the center of the explosion. Their lightcurves are dominated by emission from optically-thin cooling and by X-ray emission of the shocked ISM gas. We find that these remnants are likely to be observed in the [OIII] lambda 5007 spectral line emission or in the soft energy-band of X-rays. Finally, we discuss our results in the context of observed Galactic supernova remnants such as 3C391 and the Cygnus Loop.

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.

Massive young star clusters contain dozens or hundreds of massive stars that inject mechanical energy in the form of winds and supernova explosions, producing an outflow which expands into their surrounding medium, shocking it and forming structures called superbubbles. The regions of shocked material can have temperatures in excess of 10$^6$ K, and emit mainly in thermal X-rays (soft and hard). This X-ray emission is strongly affected by the action of thermal conduction, as well as by the metallicity of the material injected by the massive stars. We present three-dimensional numerical simulations exploring these two effects, metallicity of the stellar winds and supernova explosions, as well as thermal conduction.

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

We present 3D hydrodynamical simulations of a jet launched from the secondary star of a binary system inside a proto-planetary nebula. The secondary star moves around the primary in a close eccentric orbit. From the gasdynamic simulations we compute synthetic [NII] 6583 emission maps. Different jet axis inclinations with respect to the orbital plane, as well as different orientations of the flow with respect to the observer are considered. For some parameter combinations, we obtain structures that show point- or mirror-symmetric morphologies depending on the orientation of the flow with respect to the observer. Furthermore, our models can explain some of the emission distribution asymmetries that are summarized in the classification given by Soker & hadar (2002).

We present 3D hydrodynamical simulations of the superbubble M17, also known as the Omega nebula, carried out with the adaptive grid code yguazu'-a, which includes radiative cooling. The superbubble is modelled considering the winds of 11 individual stars from the open cluster inside the nebula (NGC 6618), for which there are estimates of the mass loss rates and terminal velocities based on their spectral types. These stars are located inside a dense interstellar medium, and they are bounded by two dense molecular clouds. We carried out three numerical models of this scenario, considering different line of sight positions of the stars (the position in the plane of the sky is known, thus fixed). Synthetic thermal X-ray emission maps are calculated from the numerical models and compared with ROSAT observations of this astrophysical object. Our models reproduce successfully both the observed X-ray morphology and the total X-ray luminosity, without taking into account thermal conduction effects.

We present new XMM-Newton observations of the supernova remnant N120 in the LMC, and numerical simulations on the evolution of this supernova remnant which we compare with the X-ray observations. The supernova remnant N120, together with several HII regions, forms a large nebular complex5D (also called N120) whose shape resembles a semicircular ring. From the XMM-Newton data we generate images and spectra of this remnant in the energy band between 0.2 to 2.0 keV. The images show that the X-ray emission is brighter towards the east (i.e., towards the rim of the large nebular complex). The EPIC/MOS1 and MOS2 data reveal a thermal spectrum in soft X-rays. 2D axisymmetric numerical simulations with the Yguaz\'u-a code were carried out assuming that the remnant is expanding into an inhomogeneous ISM with an exponential density gradient and showing that thermal conduction effects are negligible. Simulated X-ray emission maps were obtained from the numerical simulations in order to compare them with the observations. We find good agreement between the XMM-Newton data, previous optical kinematic data, and the numerical simulations; the simulations reproduce the observed X-ray luminosity and surface brightness distribution. We have also detected more extended diffuse X-ray emission probably due to the N120 large HII complex or superbubble.

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.

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.

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.

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.

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.

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 report new high resolution and high sensitivity radio observations of the extended supernova remnant (SNR) CTB 80 (G69.0+2.7) at 240 MHz, 324 MHz, 618 MHz, and 1380 MHz. The imaging of CTB 80 at 240 MHz and 618 MHz was performed using the Giant Metrewave Radio Telescope (GMRT) in India. The observations at 324 MHz and 1380 MHz were obtained using the Very Large Array (VLA, NRAO) in its C and D configurations. The new radio images reveal faint extensions for the asymmetric arms of CTB 80. The arms are irregular with filaments and clumps of size 1' (or 0.6 pc at a distance of 2 kpc). The radio image at 1380 MHz is compared with IR and optical emission. The correspondence IR/radio is excellent along the N arm of CTB 80. Ionized gas observed in the [SII] line perfectly matches the W and N edges of CTB 80. The central nebula associated with the pulsar PSR B1951+32 was investigated with an angular resolution of 10" x 6". The new radio image obtained at 618 MHz shows with superb detail structures in the 8' x 4' E-W ``plateau'' nebula that hosts the pulsar on its western extreme. A twisted filament, about 6' in extent (~3.5 pc), trails behind the pulsar in an approximate W-E direction. In the bright ``core'' nebula (size \~45"), located to the W of the plateau, the images show a distortion in the morphology towards the W; this feature corresponds to the direction in which the pulsar escapes from the SNR with a velocity of ~240 km/s. Based on the new observations, the energetics of the SNR and of the PWN are investigated.

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

The distribution and kinematics of neutral hydrogen have been studied in a wide area around the supernova remnant W28. A 2.5 x 2.5 arcdeg field centered at l = 6.5 arcdeg, b = 0 arcdeg was surveyed using the Parkes 64-m radio telescope (HPBW 14.7 arcmin at lambda 21 cm). Even though W28 is located in a complex zone of the Galactic plane, we have found different HI features which are evidence of the interaction between W28 and its surrounding gas.