The morphology of planetary nebulae emerging from the common envelope phase of binary star evolution is investigated. Using initial conditions based on the numerical results of hydrodynamical simulations of the common envelope phase it is found that the shapes and sizes of the resulting nebula are very sensitive to the effective temperature of the remnant core, the mass-loss rate at the onset of the common envelope phase, and the mass ratio of the binary system. These parameters are related to the efficiency of the mass ejection after the spiral-in phase, the stellar evolutionary phase (i.e., RG, AGB or TP-AGB), and the degree of departure from spherical symmetry in the stellar wind mass loss process itself respectively. It is found that the shapes are mostly bipolar in the early phase of evolution, but can quickly transition to elliptical and barrel-type shapes. Solutions for nested lobes are found where the outer lobes are usually bipolar and the inner lobes are elliptical, bipolar or barrel-type, a result due to the flow of the photo-evaporated gas from the equatorial region. It is found that the lobes can be produced without the need for two distinct mass ejection events. In all the computations, the bulk of the mass is concentrated in the orbital or equatorial plane, in the form of a large toroid, which can be either neutral (early phases) or photoionized (late phases), depending of the evolutionary state of the system.

High resolution observations of the extended atmospheres of asymptotic giant branch (AGB) stars can now directly confront the theories that describe stellar mass loss. Using Atacama Large Millimeter/submillimeter Array (ALMA) high angular resolution ($30\times42$~mas) observations we have, for the first time, resolved stellar rotation of an AGB star, R~Dor. We measure an angular rotation velocity of $\omega_R\sin{i}=(3.5\pm0.3)\times10^{-9}$~rad~s$^{-1}$ which indicates a rotational velocity of $|\upsilon_{\rm rot}\sin{i}|=1.0\pm0.1$~km~s$^{-1}$ at the stellar surface ($R_*=31.2$~mas at $214$~GHz). The rotation axis projected on the plane of the sky has a position angle $\Phi=7\pm6^\circ$. We find that the rotation of R Dor is two orders of magnitude faster than expected for a solitary AGB star that will have lost most of its angular momentum. Its rotational velocity is consistent with angular momentum transfer from a close companion. As a companion has not been directly detected we thus suggest R~Dor has a low-mass, close-in, companion. The rotational velocity approaches the critical velocity, set by the local sound speed in the extended envelope, and is thus expected to affect the mass loss characteristics of R~Dor.

We present new binary stellar evolution models that include the effects of tidal forces, rotation, and magnetic torques with the goal of testing Planetary Nebulae (PNe) shaping via binary interaction. We explore whether tidal interaction with a companion can spin up the AGB envelope. To do so we have selected binary systems with main sequence masses of 2.5 \Mo and of 0.8 \Mo and evolve them allowing initial separations of 5, 6, 7, and 8 AU. The binary stellar evolution models have been computed all the way to the PNe formation phase or until Roche lobe overflow (RLOF) is reached, whatever happens first. We show that with initial separations of 7 and 8 AU, the binary avoids entering into RLOF, and the AGB star reaches moderate rotational velocities at the surface ($\sim 3.5 $ and $\sim 2 $ \kms respectively) during the inter-pulse phases, but after the thermal pulses it drops to a final rotational velocity of only $\sim 0.03 $ \kms. For the closest binary separations explored, 5 and 6 AU, the AGB star reaches rotational velocities of $\sim 6 $ and $\sim 4 $ \kms respectively when the RLOF is initiated. We conclude that the detached binary models that avoid entering the RLOF phase during the AGB will not shape bipolar PNe, since the acquired angular momentum is lost via the wind during the last two thermal pulses. This study rules out tidal spin-up in non-contact binaries as a sufficient condition to form bipolar PNe.

We present high spatial resolution ($\approx$ 60--90 milliarcseconds) images of the molecular hydrogen emission in the Planetary Nebula (PN) NGC 2346. The data were acquired during the System Verification of the Gemini Multi-Conjugate Adaptive Optics System + Gemini South Adaptive Optics Imager. At the distance of NGC 2346, 700 pc, the physical resolution corresponds to $\approx$ 56 AU, which is slightly higher than that an [N II] image of NGC 2346 obtained with HST/WFPC2. With this unprecedented resolution we were able to study in detail the structure of the H$_2$ gas within the nebula for the first time. We found it to be composed of knots and filaments, which at lower resolution had appeared to be a uniform torus of material. We explain how the formation of the clumps and filaments in this PN is consistent with a mechanism in which a central hot bubble of nebular gas surrounding the central star has been depressurized, and the thermal pressure of the photoionized region drives the fragmentation of the swept-up shell.

We have computed new stellar evolution models that include the effects of rotation and magnetic torques under different hypothesis. The goal is to test if a single star can sustain in the envelope the rotational velocities needed for the magneto hydrodynamical (MHD) simulations to shape bipolar Planetary Nebulae (PNe) when the high mass-loss rates take place. Stellar evolution models with main sequence masses of 2.5 and 5 Mo, and initial rotational velocities of 250 km/s have been followed all the way to the PNe formation phase. We find that stellar cores have to be spun down using magnetic torques in order to reproduce the rotation rates observed for white dwarfs. During the asymptotic giant branch phase and beyond, the magnetic braking of the core has a practically null effect in increasing the rotational velocity of the envelope since the stellar angular momentum is removed efficiently by the wind. We have, as well, tested best possible case scenarios in rather non-physical contexts to give enough angular momentum to the envelope. We find that we cannot get the envelope of a single star rotating at the speeds needed by the MHD simulations to form bipolar PNe. We conclude that single stellar rotators are unlikely to be the progenitors of bipolar PNe under the current MHD model paradigm.

We explore the formation of multipolar structures in planetary and pre-planetary nebulae from the interaction of a fast post-AGB wind with a highly inhomogeneous and filamentary shell structure assumed to form during the final phase of the high density wind. The simulations were performed with a new hydrodynamics code integrated in the interactive framework of the astrophysical modeling package SHAPE. In contrast to conventional astrophysical hydrodynamics software, the new code does not require any programming intervention by the user for setting up or controlling the code. Visualization and analysis of the simulation data has been done in SHAPE without external software. The key conclusion from the simulations is that secondary lobes in planetary nebulae, such as Hubble 5 and K3-17, can be formed through the interaction of a fast low-density wind with a complex high density environment, such as a filamentary circumstellar shell. The more complicated alternative explanation of intermittent collimated outflows that change direction, in many cases may therefore not be necessary. We consider that the wind-shell interaction scenario is more likely since the bow-shock shape expected from a strongly cooling bow-shock from jets is different from that of the observed bubbles. Furthermore, the timescales of the wind-wind interaction suggest that the progenitor star was rather massive.

We study the hydrodynamical behavior of the gas expelled by moving Asymptotic Giant Branch Stars interacting with the ISM. Our models follow the wind modulations prescribed by stellar evolution calculations, and we cover a range of expected relative velocities (10 to 100 km/s), ISM densities (between 0.01 and 1 cm-3), and stellar progenitor masses (1 and 3.5 Msun). We show how and when bow-shocks, and cometary-like structures form, and in which regime the shells are subject to instabilities. Finally, we analyze the results of the simulations in terms of the different kinematical stellar populations expected in the Galaxy.

We explore the effects of photodissociation at the stages of post-asymptotic giant branch stars to find a mechanism able to produce multi-polar shapes. We perform two-dimensional gasdynamical simulations to model the effects of photodissociation in proto-planetary nebulae. We find that post-asymptotic giant branch stars with 7,000 K or hotter are able to photodissociate a large amount of the circumstellar gas. We compute several solutions for nebulae with low-velocity multi-lobes. We find that the early expansion of a dissociation front is crucial to understand the number of lobes in proto-planetary nebulae. A dynamical instability appears when cooling is included in the swept-up molecular shell. This instability is similar to the one found in photoionization fronts, and it is associated with the thin-shell Vishniac instability. The dissociation front exacerbates the growth of the thin-shell instability, creating a fast fragmentation in shells expanding into media with power-law density distributions such as r^-2.

Multi wavelength observations of Cassiopeia A (Cas A) have provided us with a strong evidence for the presence of circumstellar material surrounding the progenitor star. It has been suggested that its progenitor was a massive star with a strong mass loss. But, despite the large amount of observational data from optical, IR, radio and x-ray observations, the identity of Cas A progenitor is still elusive. In this work, we compute stellar and circumstellar numerical models to look for the progenitor of Cas A. The models will be compared with the observational constraints. We have computed stellar evolution models to get time-dependent wind parameters and surface abundances. We have chosen a set of probable progenitor stars, with initial masses of 23, 28, 29, 30 and 33 \Mo, with mass loss. The derived mass loss rates and wind terminal velocities are used to simulate the evolution of the circumstellar medium. Our stellar set gives distinct SN progenitors: RSG, luminous blue super giants (LBSGs) and Wolf-Rayet (WR) stars. The 23 \Mo star explodes as a RSG in a $\rm \rho \sim r^{-2}$ dense, free streaming wind surrounded by a thin, compressed, RSG shell. The 28 $\Mo$ star explodes as a LBSG, and the SN blast wave interacts with a low density, free streaming wind surrounded by an unstable and massive ''RSG+LBSG'' shell. Finally, the 30 and 33 \Mo stars explode as WR stars surrounded by fast, WR winds that terminate in highly fragmented ''WR+RSG shell''. We have compared the surface chemical abundances of our stellar models with the observational abundances in Cas A. The abundance analysis shows that the progenitor was a star with an initial mass of the order of 30 \Mo, while the hydrodynamical analysis favors progenitors of initial masses around 23.

Based on axi-symmetric hydrodynamical simulations and 3D reconstructions with Shape, we investigate the kinematic signatures of deviations from homologous ("Hubble-type") outflows in some typical shapes of planetary nebulae. We find that, in most situations considered in our simulations, the deviations from a Hubble-type flow are significant and observable. The deviations are systematic and a simple parameterization of them considerably improves morpho-kinematical models of the simulations. We describe such extensions to a homologous expansion law that capture the global velocity structure of hydrodynamical axi-symmetric nebulae during their wind-blown phase. It is the size of the poloidal velocity component that strongly influences the shape of the position velocity diagrams that are obtained, not so much the variation of the radial component. The deviations increase with the degree of collimation of the nebula and they are stronger at intermediate latitudes. We describe potential deformations which these deviations might produce in 3D reconstructions that assume "Hubble-type" outflows. The general conclusion is that detailed morpho-kinematic observations and modeling of planetary nebulae can reveal whether a nebula is still in a hydrodynamically active stage (windy phase) or whether it has reached ballistic expansion.

In the course of the life of a massive star, wind-wind interaction can give rise to the formation of circumstellar nebulae which are both predicted and observed in the nature. We present generic model calculations to predict the properties of such nebulae for blue supergiants. From stellar evolution calculations including rotation, we obtain the time dependence of the stellar wind properties and of the stellar radiation field. These are used as input for hydro-calculations of the circumstellar medium throughout the star's life. Here, we present the results for a rapidly rotating 12 solar masses single star. This star undergoes a blue loop during its post main sequence evolution, at the onset of which its contraction spins it up close to critical rotation. Due to the consequent anisotropic mass loss, the blue supergiant wind sweeps up the preceding slow wind into an hour glass structure. Its collision with the previously formed spherical red supergiant wind shell forms a short-lived luminous nebula consisting of two polar caps and a central inner ring. With time, the polar caps evolve into mid-latitude rings which gradually move toward the equatorial plane while the central ring is fading. These structures are reminiscent to the observed nebulae around the blue supergiant Sher 25 and the progenitor of SN 1987A. The simple model of an hour glass colliding with a spherical shell retrieves most of the intriguing nebula geometries discovered around blue supergiants, and suggests them to form an evolutionary sequence. Our results indicate that binarity is not required to obtain them.

The central star of this nebula has an observed intense magnetic field and the fast wind is no longer present, indicating that a back flow process has probably developed. Long-slit, spatially resolved echelle spectra have been obtained across the main body of NGC 1360 and over its system of bipolar jets. Deep images of the knotty structures of the jets have also been obtained. The data allow a detailed study of the structure and kinematics of this object and the results are modeled considering the effects of a magnetic collimation process in the development of the nebula and then switching off the fast stellar wind to follow its evolution to its current state. The model is able to successfully reproduce many of the key features of NGC 1360 under these premises.

Rapidly rotating, chemically homogeneously evolving massive stars are considered to be progenitors of long gamma-ray bursts. We present numerical simulations of the evolution of the circumstellar medium around a rapidly rotating 20 Msol star at a metallicity of Z=0.001. Its rotation is fast enough to produce quasi-chemically homogeneous evolution. While conventionally, a star of 20 Msol would not evolve into a Wolf-Rayet stage, the considered model evolves from the main sequence directly to the helium main sequence. We use the time-dependent wind parameters, such as mass loss rate, wind velocity and rotation-induced wind anisotropy from the evolution model as input for a 2D hydrodynamical simulation. While the outer edge of the pressure-driven circumstellar bubble is spherical, the circumstellar medium close to the star shows strong non-spherical features during and after the periods of near-critical rotation. We conclude that the circumstellar medium around rapidly rotating massive stars differs considerably from the surrounding material of non-rotating stars of similar mass. Multiple blue-shifted high velocity absorption components in gamma-ray burst afterglow spectra are predicted. As a consequence of near critical rotation and short stellar evolution time scales during the last few thousand years of the star's life, we find a strong deviation of the circumstellar density profile in the polar direction from the 1/R^2 density profile normally associated with stellar winds close to the star

Van Marle et al. (2005) showed that circumstellar absorption lines in early Type Ib/c supernova and gamma-ray burst afterglow spectra may reveal the progenitor evolution of the exploding Wolf-Rayet star. While the quoted paper deals with Wolf-Rayet stars which evolved through a red supergiant stage, we investigate here the initially more massive Wolf-Rayet stars which are thought to evolve through a Luminous Blue Variable (LBV) stage. We perform hydrodynamic simulations of the evolution of the circumstellar medium around a 60 Msol star, from the main sequence through the LBV and Wolf-Rayet stages, up to core collapse. We then compute the column density of the circumstellar matter as a function of radial velocity, time and angle. This allows a comparison with the number and blue-shifts, of absorption components in the spectra of LBVs, Wolf-Rayet stars, Type Ib/c supernovae and gamma-ray burst afterglows. Our simulation for the post-LBV stage shows the formation of various absorption components, which are, however, rather short lived; they dissipate on time scales shorter than 50,000yr. As the LBV stage is thought to occur at the beginning of core helium burning, the remaining Wolf-Rayet life time is expected to be one order of magnitude larger. When interpreting the absorption components in the afterglow spectrum of GRB-021004 as circumstellar, it can be concluded that the progenitor of this source did most likely not evolve through an LBV stage. However, a close binary with late common-envelope phase (Case C) may produce a circumstellar medium that closely resembles the LBV to Wolf-Rayet evolution, but with a much shorter Wolf-Rayet period.

We present hydrodynamical models of circumstellar medium (CSM) of long gamma-ray burst (GRB) progenitor candidates. These are massive stars that have lost a large amount of mass in the form of stellar wind during their evolution. There are two possible ways to probe the CSM of long GRB progenitors. Firstly, the GRB afterglow consists of synchrotron radiation, emitted when the GRB jet sweeps up the surrounding medium. Therefore, the lightcurve is directly related to the density profile of the CSM. The density can either decrease with the radius squared (as is the case for a freely expanding stellar wind) or be constant (as we would expect for shocked wind or the interstellar medium). Secondly, material between the GRB and the observer will absorb part of the afterglow radiation, causing absorption lines in the afterglow spectrum. In some cases, such absorption lines are blue-shifted relative to the source indicating that the material is moving away from the progenitor star. This can be explained in terms of wind interactions in the CSM. We can use the CSM of these stars to investigate their prior evolutionary stage.

The progenitor stars of long Gamma-Ray Bursts (GRBs) are thought to be Wolf-Rayet stars, which generate a massive and energetic wind. Nevertheless, about 25 percent of all GRB afterglows light curves indicate a constant density medium close to the exploding star. We explore various ways to produce this, by creating situations where the wind termination shock arrives very close to the star, as the shocked wind material has a nearly constant density. Typically, the distance between a Wolf-Rayet star and the wind termination shock is too large to allow afterglow formation in the shocked wind material. Here, we investigate possible causes allowing for a smaller distance: A high density or a high pressure in the surrounding interstellar medium (ISM), a weak Wolf-Rayet star wind, the presence of a binary companion, and fast motion of the Wolf-Rayet star relative to the ISM. We find that all four scenarios are possible in a limited parameter space, but that none of them is by itself likely to explain the large fraction of constant density afterglows. A low GRB progenitor metallicity, and a high GRB energy make the occurrence of a GRB afterglow in a constant density medium more likely. This may be consistent with constant densities beingpreferentially found for energetic, high redshift GRBs.

In this paper we explore the dynamics of ionization bounded planetary nebulae after the termination of the fast stellar wind. When the stellar wind becomes negligible, the hot, shocked bubble depressurizes and the thermal pressure of the photoionized region, at the inner edge of the swept-up shell, becomes dominant. At this stage the shell tends to fragment creating clumps with comet-like tails and long, photoionized trails in between, while the photoionized material expands back towards the central stars as a rarefaction wave. Once that the photoionized gas fills the inner cavity, it develops a kinematical pattern of increasing velocity from the center outwards with a typical range of velociti es starting from the systemic velocity to 50 Km/s at the edges. The Helix nebula is a clear example of a planetary nebula at this late evolutionary stage.

The dynamical evolution of HII regions with and without stellar motion in dense, structured molecular clouds is studied. Clouds are modeled in hydrostatic equilibrium, with gaussian central cores and external halos that obey r**-2 and r**-3 density power laws. The cloud gravity is included as a time-independent, external force. Stellar velocities of 0, 2, 8, and 12 km/s are considered. When stellar motion is included, stars move from the central core to the edge of the cloud, producing transitions from ultracompact to extended HII regions as the stars move into lower density regions. The opposite behavior occurs when stars move toward the cloud cores. The main conclusion of our study is that ultracompact HII regions are pressure-confined entities while they remain embedded within dense cores. The confinement comes from ram and/or ambient pressures. The survival of ultracompact regions depends on the position of the star with respect to the core, the stellar life-time, and the core crossing time. Stars with velocities less than the cloud dispersion velocity can produce cometary shapes smaller than 0.1 pc at times of 20,000 yr or more. The sequence Ultracompact to Compact to Extended HII region shows a variety of unpredictable structures due to ionization-shock front instability. Some ultracompact HII regions with a core-halo morphology might be explained by self-blocking effects, when stars overtake and ionize leading, piled-up clumps of neutral gas. We use thermal energy to support the cloud against gravity; the results remain the same if other types of isotropic cloud support are used.

Long gamma-ray bursts are thought to be caused by a subset of exploding Wolf-Rayet stars. We argue that the circumstellar absorption lines in early supernova and in gamma-ray burst afterglow spectra may allow us to determine the main properties of the Wolf-Rayet star progenitors which can produce those two events. To demonstrate this, we first simulate the hydrodynamic evolution of the circumstellar medium around a 40 Msun star up to the time of the supernova explosion. Knowledge of density, temperature and radial velocity of the circumstellar matter as function of space and time allows us to compute the column density in the line of sight to the centre of the nebula, as a function of radial velocity, angle, and time. Our column density profiles indicate the possible number, strengths, widths and velocities of absorption line components in supernova and gamma-ray burst afterglow spectra. Our example calculation shows four distinct line features during the Wolf-Rayet stage, at about 0, 50, 150-700 and 2200 km/s, with only those of the lowest and highest velocity present at all times. The 150-700 km/s feature decays rapidly as function of time after the onset of the Wolf-Rayet stage. It consists of a variable number of components, and, especially in its evolved stage, is depending strongly on the particular line of sight. A comparison with absorption lines detected in the afterglow of GRB 021004 suggests that the high velocity absorption component in GRB 021004 may be attributed to the free streaming Wolf-Rayet wind, which is consistent with the steep density drop indicated by the afterglow light curve. The presence of the intermediate velocity components implies that the duration of the Wolf-Rayet phase of the progenitor of GRB 021004 was much smaller than the average Wolf-Rayet life time.

We present a numerical investigation of the contribution of the presupernova ejecta of Wolf-Rayet stars to the environment surrounding gamma-ray bursts (GRBs), and describe how this external matter can affect the observable afterglow characteristics. An implicit hydrodynamic calculation for massive stellar evolution is used here to provide the inner boundary conditions for an explicit hydrodynamical code to model the circumstellar gas dynamics. The resulting properties of the circumstellar medium are then used to calculate the deceleration of a relativistic, gas-dynamic jet and the corresponding afterglow light curve produced as the shock wave propagates through the shocked-wind medium. We find that variations in the stellar wind drive instabilities that may produce radial filaments in the shocked-wind region. These comet-like tails of clumps could give rise to strong temporal variations in the early afterglow lightcurve. Afterglows may be expected to differ widely among themselves, depending on the angular anisotropy of the jet and the properties of the stellar progenitor; a wide diversity of behaviors may be the rule, rather than the exception.

This paper explores the effects of post-AGB winds driven solely by magnetic pressure from the stellar surface. It is found that winds can reach high speeds under this assumption, and lead to the formation of highly collimated proto-planetary nebulae. Bipolar knotty jets with periodic features and constant velocity are well reproduced by the models. Several wind models with terminal velocities from a few tens of $\kms$ up to $10^3$ $\kms$ are calculated, yielding outflows with linear momenta in the range $10^{36}-10^{40} \gcms$, and kinetic energies in the range $10^{42}-10^{47} $ erg. These results are in accord with recent observations of proto-planetary nebulae that have pointed out serious energy and momentum deficits if radiation pressure is considered as the only driver for these outflows. Our models strengthen the notion that the large mass-loss rates of post-AGB stars, together with the short transition times from the late AGB to the planetary nebula stage, could be directly linked with the generation of strong magnetic fields during this transition stage.

High ambient interstellar pressure is suggested as a possible factor to explain the ubiquitous observed growth-rate discrepancy for supernova-driven superbubbles and stellar wind bubbles. Pressures of P/k ~ 1e5 cm-3 K are plausible for regions with high star formation rates, and these values are intermediate between the estimated Galactic mid-plane pressure and those observed in starburst galaxies. High-pressure components also are commonly seen in Galactic ISM localizations. We demonstrate the sensitivity of shell growth to the ambient pressure, and suggest that superbubbles ultimately might serve as ISM barometers.

S308 is a circumstellar bubble blown by the WN4 star HD50896. It is one of the only two single-star bubbles that show detectable diffuse X-ray emission. We have obtained XMM-Newton EPIC observations of the northwest quadrant of S308. The diffuse X-ray emission shows a limb-brightened morphology, with a clear gap extending from the outer edge of the diffuse X-ray emission to the outer rim of the nebular shell. The X-ray spectrum of the diffuse emission is very soft, and is well fitted by an optically thin plasma model for a N-enriched plasma at temperatures of ~1.1x10^6 K. A hotter gas component may exist but its temperature is not well constrained as it contributes less than 6% of the observed X-ray flux. The total X-ray luminosity of S308, extrapolated from the bright northwest quadrant, is <=(1.2+-0.5)x10^{34} ergs/s. We have used the observed bubble dynamics and the physical parameters of the hot interior gas of S308 in conjunction with the circumstellar bubble model of Garcia-Segura & Mac Low (1995) to demonstrate that the X-ray-emitting gas must be dominated by mixed-in nebular material.

We present two-dimensional numerical simulations of the evolution of a low-mass star moving supersonically through its surrounding interstellar medium (ISM). We show that the ejecta of a moving star with a systemic velocity of 20 km/s will interact with the ISM and will form bow-shock structures qualitatively similar to what is observed. We find that, due to ram-pressure stripping, most of the mass ejected during the AGB phase is left downstream of the moving star. As a consequence, the formation of the PN is highly influenced, even at the low relative velocity of the star. The models are based on the predictions of stellar evolution calculations. Therefore, the density and velocity of the AGB and post-AGB winds are time dependent and give rise to the formation of shock regions inside the cavity formed by the previous winds. As a result, the stand-off distance is also time dependent and cannot be determined by simple analytical arguments.

We have studied the effect of the mass of the central star (CS) on the gas evolution during the planetary nebula (PN) phase. We have performed numerical simulations of PN formation using CS tracks for six stellar core masses corresponding to initial masses from 1 to 5 Msun. The gas structure resulting from the previous asymptotic giant branch (AGB) evolution is used as the starting configuration. The formation of multiple shells is discussed in the light of our models, and the density, velocity and halpha emission brightness profiles are shown for each stellar mass considered. We have computed the evolution of the different shells in terms of radius, expansion velocity, and halpha peak emissivity. We find that the evolution of the main shell is controlled by the ionization front rather than by the thermal pressure provided by the hot bubble during the early PN stages. This effect explains why the kinematical ages overestimate the age in young CSs. At later stages in the evolution and for low mass progenitors the kinematical ages severely underestimate the CS age. Large (up to 2.3 pc), low surface brightness shells (less than 2000 times the brightness of the main shell) are formed in all of our models (with the exception of the 5 Msun model). These PN halos contain most of the ionized mass in PNe, which we find is greatly underestimated by the observations because of the low surface brightness of the halos.

We have investigated the dynamical interaction of low- and-intermediate mass stars (from 1 to 5 Msun) with their interstellar medium (ISM). In this first paper, we examine the structures generated by the stellar winds during the Asymptotic Giant Branch (AGB) phase, using a numerical code and the wind history predicted by stellar evolution. The influence of the external ISM is also taken into account. We find that the wind variations associated with the thermal pulses lead to the formation of transient shells with an average lifetime of 20,000 yr, and consequently do not remain recorded in the density or velocity structure of the gas. The formation of shells that survive at the end of the AGB occurs via two main processes: shocks between the shells formed by two consecutive enhancements of the mass-loss or via continuous accumulation of the material ejected by the star in the interaction region with the ISM. Our models show that the mass of the circumstellar envelope increases appreciably due to the ISM material swept up by the wind (up to 70 % for the 1 Msun stellar model). We also point out the importance of the ISM on the deceleration and compression of the external shells. According to our simulations, large regions (up to 2.5 pc) of neutral gas surrounding the molecular envelopes of AGB stars are expected. These large regions of gas are formed from the mass-loss experienced by the star during the AGB evolution.

The Hubble Space Telescope has revealed the existence of multiple, regularly spaced, and faint concentric shells around some planetary nebulae. Here we present 2(1/2)D magnetohydrodinamic numerical simulations of the effects of a solar-like magnetic cycle, with periodic polarity inversions, in the slow wind of an AGB star. The stellar wind is modeled with a steady mass-loss at constant velocity. This simple version of a solar-like cycle, without mass-loss variations, is able to reproduce many properties of the observed concentric rings. The shells are formed by pressure oscillations, which drive compressions in the magnetized wind. These pressure oscillations are due to periodic variations in the field intensity. The periodicity of the shells, then, is simply a half of the magnetic cycle since each shell is formed when the magnetic pressure goes to zero during the polarity inversion. As a consequence of the steady mass-loss rate, the density of the shells scales as $r^{-2}$, and their surface brightness has a steeper drop-off, as observed in the shells of NGC 6543, the best documented case of these HST rings. Deviations from sphericity can be generated by changing the strength of the magnetic field. For sufficiently strong fields, a series of symmetric and equisdistant blobs are formed at the polar axis, resembling the ones observed in He 2-90. These blobs are originated by magnetic collimation within the expanding AGB wind.

Multiple shell planetary nebulae (MSPNe) are most likely a consequence of the modulation of the mass-loss rate during the thermal pulses on the AGB. By using numerical simulations of their formation, the models' predictions and the history of the winds can be investigated. As their halos have low densities and expansion velocities, MSPNe are expected to be affected by interaction with the ISM . This would then give us evidence of the local ISM conditions. In order to study the formation of MSPNe, we have performed numerical simulations following the evolution of the stellar winds for a 1 solar mass star from the AGB to the post-AGB stages. Without invoking any asymmetry for the stellar wind and taking into account the effects of a moving central star, an asymmetric halo is formed as a consequence of the interaction with the ISM. We found that the asymmetries caused by the interaction take place from the beginning of the evolution and have an enormous influence on the formation of the halo.

The observed giant outbursts of Luminous Blue Variables (LBVs) may occur when these massive stars approach their Eddington limits. When this happens, they must reach a point where the centrifugal force and the radiative acceleration cancel out gravity at the equator. We call this the Omega-limit. When stars are close to the Omega-limit, strong non-spherical mass loss should occur. This suggests a scenario where a slow and very dense wind, strongly confined to the equatorial plane, is followed by a fast and almost spherical wind. We compute two-dimensional hydrodynamic models of the evolution of the nebula formed from such interacting winds, using parameters consistent with the outburst of eta Carinae in the last century. This outburst gave birth to the Homunculus, the hourglass-shaped inner part of a highly structured circumstellar nebula. Assuming the star was very close to the Omega-limit during outburst, our models produce gas distributions that strongly resemble the Homunculus on large and small scale. This supports the general conjecture that giant outbursts in LBVs occur when they approach the Eddington limit. Our models constrains the average mass loss rate since the outburst to values smaller than the present-day mass loss rate and suggest that eta Car is approaching another outburst. Our models imply that the occurrence of giant LBV outbursts depends on the initial stellar rotation rate, and that the initial angular momentum is as important to the evolution of very massive stars as their initial mass or metallicity.

Here we discuss the interaction of massive stars with their parental molecular clouds. A summary of the dynamical evolution of HII regions and wind-driven bubbles in high-pressure cloud cores is given. Both ultracompact HII regions and ultracompact wind-driven bubbles can reach pressure equilibrium with their surrounding medium. The structures stall their expansion and become static and, as long as the ionization sources and the ambient densities remain about constant, the resulting regions are stable and long lived. For cases with negative density gradients, and depending on the density distribution, some regions never reach the static equilibrium condition. For power-law density stratifications, $\rho \propto r^{-w}$, the properties of the evolution depend on a critical exponent, $w_{crit}$, above which the ionization front cannot be slowed down by recombinations or new ionizations, and the cloud becomes fully ionized. This critical exponent is $w_{crit}=3/2$ during the expansion phase. For $w>3/2 $ the gas expands supersonically into the surrounding ionized medium, and there are two regimes separated by $w=3$. For $3/2<w\leq 3$, the slow regime, the inner region drives a weak shock moving with almost constant velocity through the cloud. For $w>3$, the fast regime, the shock becomes strong and accelerates with time. Finally, the evolution of slow winds in highly pressurized region is described briefly.

We present the main features in the evolution of the gas pressure in star forming regions, from the formation of the parental cloud to the moment when the region is pressurized by interacting stellar winds. The main processes for cloud destruction and the self-limiting properties of star formation are described. The high star forming efficiency in nuclear starbursts is a consequence of the high gas pressures. The evolution of slow winds in highly pressurized region is also sketched.