The complexity of physico-chemical models of star formation is increasing, with models that take into account new processes and more realistic setups. These models allow astrochemists to compute the evolution of chemical species throughout star formation. Hence, comparing the outputs of such models to observations allows to bring new constraints on star formation. The work presented in this paper is based on the recent public release of a database of radiation hydrodynamical low-mass star formation models. We used this database as physical parameters to compute the time dependent chemical composition of collapsing cores with a 3-phase gas-grain model. The results are analyzed to find chemical tracers of the initial physical parameters of collapse such as the mass, radius, temperature, density, and free-fall time. They are also compared to observed molecular abundances of Class 0 protostars. We find numerous tracers of the initial parameters of collapse, except for the initial mass. More particularly, we find that gas phase CH3CN, NS and OCS trace the initial temperature, while H2CS trace the initial density and free-fall time of the parent cloud. The comparison of our results with a sample of 12 Class 0 low mass protostars allows us to constrain the initial parameters of collapse of low-mass prestellar cores. We find that low-mass protostars are preferentially formed within large cores with radii greater than 20000 au, masses between 2 and 4 Msol, temperatures lower or equal to 15 K, and densities between 6e4 and 2.5e5 part.cm-3, corresponding to free-fall times between 100 and 200 kyrs.

From molecular clouds to protoplanetary disks, non-ideal magnetic effects are important in many astrophysical environments. Indeed, in star and disk formation processes, it has become clear that these effects are critical to the evolution of the system. The efficacy of non-ideal effects are, however, determined by the complex interplay between magnetic fields, ionising radiation, cosmic rays, microphysics, and chemistry. In order to understand these key microphysical parameters, we present a one-dimensional non-ideal magnetohydrodynamics code and apply it to a model of a time-dependent, oblique, magnetic shock wave. By varying the microphysical ingredients of the model, we find that cosmic rays and dust play a major role, and that, despite the uncertainties, the inclusion of microphysics is essential to obtain a realistic outcome in magnetic astrophysical simulations.

The current study was developed to provide a database of relatively simple numerical simulations of protostellar collapse, as a template library for observations of cores and very young protostars, and for researchers who wish to test their chemical modeling under dynamic astrophysical conditions. It was also designed to identify statistical trends that may appear when running many models of the formation of low-mass stars by varying the initial conditions. A large set of 143 calculations of the gravitational collapse of an isolated sphere of gas with uniform temperature and a Bonnor-Ebert like density profile was undertaken using a 1D fully implicit Lagrangian radiation hydrodynamics code. The parameter space covered initial masses from 0.2 to 8 Msun, temperatures of 5-30 K and radii between 3000 and 30,000 AU. A spread in the thermal evolutionary tracks of the runs was found, due to differing initial conditions and optical depths. Within less than an order of magnitude, all first and second Larson cores had masses and radii independent of the initial conditions. The time elapsed between the formation of the first and second cores was found to strongly depend on the first core mass accretion rate, and no first core in our grid of models lived for longer than 2000 years, before the onset of the second collapse. The end product of a protostellar cloud collapse, the second Larson core, is, at birth, a canonical object with a mass and radius of about 3 Mjup and 8 Rjup, independent of its initial conditions. The evolution sequence which brings the gas to stellar densities can however proceed in a variety of scenarios, on different timescales, along different isentropes, but each story line can largely be predicted by the initial conditions. All the data from the simulations are publicly available at this address: http://starformation.hpc.ku.dk/grid-of-protostars.

The transport of angular momentum is capital during the formation of low-mass stars; too little removal and rotation ensures stellar densities are never reached, too much and the absence of rotation means no protoplanetary disks can form. Magnetic diffusion is seen as a pathway to resolving this long-standing problem. We investigate the impact of including resistive MHD in simulations of the gravitational collapse of a 1 solar mass gas sphere, from molecular cloud densities to the formation of the protostellar seed; the second Larson core. We used the AMR code RAMSES to perform two 3D simulations of collapsing magnetised gas spheres, including self-gravity, radiative transfer, and a non-ideal gas equation of state to describe H2 dissociation which leads to the second collapse. The first run was carried out under the ideal MHD approximation, while ambipolar and ohmic diffusion was incorporated in the second calculation. In the ideal MHD simulation, the magnetic field dominates the energy budget everywhere inside and around the first core, fueling interchange instabilities and driving a low-velocity outflow. High magnetic braking removes essentially all angular momentum from the second core. On the other hand, ambipolar and ohmic diffusion create a barrier which prevents amplification of the magnetic field beyond 0.1 G in the first Larson core which is now fully thermally supported. A significant amount of rotation is preserved and a small Keplerian-like disk forms around the second core. When studying the radiative efficiency of the first and second core accretion shocks, we found that it can vary by several orders of magnitude over the 3D surface of the cores. Magnetic diffusion is a pre-requisite to star-formation; it enables the formation of protoplanetary disks in which planets will eventually form, and also plays a determinant role in the formation of the protostar itself.

Here we report that the most rapidly recurring nova, M31N 2008-12a, which erupts annually, is surrounded by a "nova super-remnant" which demonstrates that M31N 2008-12a has erupted with high frequency for millions of years.

We develop a detailed chemical network relevant to the conditions characteristic of prestellar core collapse. We solve the system of time-dependent differential equations to calculate the equilibrium abundances of molecules and dust grains, with a size distribution given by size-bins for these latter. These abundances are used to compute the different non-ideal magneto-hydrodynamics resistivities (ambipolar, Ohmic and Hall), needed to carry out simulations of protostellar collapse. For the first time in this context, we take into account the evaporation of the grains, the thermal ionisation of Potassium, Sodium and Hydrogen at high temperature, and the thermionic emission of grains in the chemical network, and we explore the impact of various cosmic ray ionisation rates. All these processes significantly affect the non-ideal magneto-hydrodynamics resistivities, which will modify the dynamics of the collapse. Ambipolar diffusion and Hall effect dominate at low densities, up to n_H = 10^12 cm^-3, after which Ohmic diffusion takes over. We find that the time-scale needed to reach chemical equilibrium is always shorter than the typical dynamical (free fall) one. This allows us to build a large, multi-dimensional multi-species equilibrium abundance table over a large temperature, density and ionisation rate ranges. This table, which we make accessible to the community, is used during first and second prestellar core collapse calculations to compute the non-ideal magneto-hydrodynamics resistivities, yielding a consistent dynamical-chemical description of this process.

In this paper, we provide a more accurate description of the evolution of the magnetic flux redistribution during prestellar core collapse by including resistive terms in the magnetohydrodynamics (MHD) equations. We focus more particularly on the impact of ambipolar diffusion. We use the adaptive mesh refinement code RAMSES to carry out such calculations. The resistivities required to calculate the ambipolar diffusion terms were computed using a reduced chemical network of charged, neutral and grain species. The inclusion of ambipolar diffusion leads to the formation of a magnetic diffusion barrier in the vicinity of the core, preventing accumulation of magnetic flux in and around the core and amplification of the field above 0.1G. The mass and radius of the first Larson core remain similar between ideal and non-ideal MHD models. This diffusion plateau has crucial consequences on magnetic braking processes, allowing the formation of disk structures. Magnetically supported outflows launched in ideal MHD models are weakened when using non-ideal MHD. Contrary to ideal MHD misalignment between the initial rotation axis and the magnetic field direction does not significantly affect the results for a given mu, showing that the physical dissipation truly dominate over numerical diffusion. We demonstrate severe limits of the ideal MHD formalism, which yield unphysical behaviours in the long-term evolution of the system. This includes counter rotation inside the outflow, interchange instabilities, and flux redistribution triggered by numerical diffusion, none observed in non-ideal MHD. Disks with Keplerian velocity profiles form in all our non-ideal MHD simulations, with final mass and size which depend on the initial magnetisation. This ranges from a few 0.01 solar masses and 20-30 au for the most magnetised case (mu=2) to 0.2 solar masses and 40-80 au for a lower magnetisation (mu=5).

Radiative transfer plays a key role in the star formation process. Due to a high computational cost, radiation-hydrodynamics simulations performed up to now have mainly been carried out in the grey approximation. In recent years, multi-frequency radiation-hydrodynamics models have started to emerge, in an attempt to better account for the large variations of opacities as a function of frequency. We wish to develop an efficient multigroup algorithm for the adaptive mesh refinement code RAMSES which is suited to heavy proto-stellar collapse calculations. Due to prohibitive timestep constraints of an explicit radiative transfer method, we constructed a time-implicit solver based on a stabilised bi-conjugate gradient algorithm, and implemented it in RAMSES under the flux-limited diffusion approximation. We present a series of tests which demonstrate the high performance of our scheme in dealing with frequency-dependent radiation-hydrodynamic flows. We also present a preliminary simulation of a three-dimensional proto-stellar collapse using 20 frequency groups. Differences between grey and multigroup results are briefly discussed, and the large amount of information this new method brings us is also illustrated. We have implemented a multigroup flux-limited diffusion algorithm in the RAMSES code. The method performed well against standard radiation-hydrodynamics tests, and was also shown to be ripe for exploitation in the computational star formation context.

Hydrogen molecules (H2) come in two forms in the interstellar medium, ortho- and para-hydrogen, corresponding to the two different spin configurations of the two hydrogen atoms. The relative abundances of the two flavours in the interstellar medium are still very uncertain, and this abundance ratio has a significant impact on the thermal properties of the gas. In the context of star formation, theoretical studies have recently adopted two different strategies when considering the ortho:para ratio (OPR) of H2 molecules; the first considers the OPR to be frozen at 3:1 while the second assumes that the species are in thermal equilibrium. As the OPR potentially affects the protostellar cores which form as a result of the gravitational collapse of a dense molecular cloud, the aim of this paper is to quantify precisely what role the choice of OPR plays in the properties and evolution of the cores. We used two different ideal gas equations of state for a hydrogen and helium mix in a radiation hydrodynamics code to simulate the collapse of a dense cloud and the formation of the first and second Larson cores; the first equation of state uses a fixed OPR of 3:1 while the second assumes thermal equilibrium. Simulations using an equilibrium ratio collapse faster at early times and show noticeable oscillations around hydrostatic equilibrium, to the point where the core expands for a short time right after its formation before resuming its contraction. In the case of a fixed 3:1 OPR, the core's evolution is a lot smoother. The OPR was however found to have little impact on the size, mass and radius of the two Larson cores. We conclude that if one is solely interested in the final properties of the cores when they are formed, it does not matter which OPR is used. On the other hand, if one's focus lies primarily in the evolution of the first core, the choice of OPR becomes important.

Star formation begins with the gravitational collapse of a dense core inside a molecular cloud. As the collapse progresses, the centre of the core begins to heat up as it becomes optically thick. The temperature and density in the centre eventually reach high enough values where fusion reactions can ignite; the protostar is born. This sequence of events entail many physical processes, of which radiative transfer is of paramount importance. Many simulations of protostellar collapse make use of a grey treatment of radiative transfer coupled to the hydrodynamics. However, interstellar gas and dust opacities present large variations as a function of frequency. In this paper, we follow-up on a previous paper on the collapse and formation of Larson's first core using multigroup radiation hydrodynamics (Paper I) by extending the calculations to the second phase of the collapse and the formation of Larson's second core. We have made the use of a non-ideal gas equation of state as well as an extensive set of spectral opacities in a spherically symmetric fully implicit Godunov code to model all the phases of the collapse of a 0.1, 1 and 10 solar mass cloud cores. We find that, for a same central density, there are only small differences between the grey and multigroup simulations. The first core accretion shock remains supercritical while the shock at the second core border is found to be strongly subcritical with all the accreted energy being transfered to the core. The size of the first core was found to vary somewhat in the different simulations (more unstable clouds form smaller first cores) while the size, mass and temperature of the second cores are independent of initial cloud mass, size and temperature. Our simulations support the idea of a standard (universal) initial second core size of 0.003 AU and mass 0.0014 solar masses.

We present the first detailed spatio-kinematical analysis and modelling of the planetary nebula Abell 65, which is known to host a post-common envelope, binary, central star system. As such, this object is of great interest in studying the link between nebular morphology and central star binarity. [OIII]5007A and H-alpha+[NII]6584A longslit spectra and imagery of Abell 65 were obtained with the Manchester Echelle Spectrometer on the 2.1-m telescope at the San Pedro Martir Observatory (MES-SPM). Further [OIII]5007A longslit spectra were obtained with the Ultraviolet and Visual Echelle Spectrograph on the Very Large Telescope (VLT-UVES). These data were used to develop a spatio-kinematical model for the [OIII]5007A emission from Abell 65. A "best-fit" model was found by comparing synthetic spectra and images rendered from the model to the data. The model comprises an outer shell and an inner shell, with kinematical ages of 15000 +/- 5000 yr kpc^-1 and 8000 +/- 3000 yr kpc^-1, respectively. Both shells have peanut-shaped bipolar structures with symmetry axes at inclinations of 55 +/- 10 deg (to the line-of-sight) for the outer shell and 68 +/- 10 deg for the inner shell. The near-alignment between the nebular shells and the binary orbital inclination (of 68 +/- 2 deg) is strongly indicative that the binary is responsible for shaping the nebula. Abell 65 is one of a growing number of planetary nebulae (seven to date, including Abell 65 itself) for which observations and modelling support the shaping influence of a central binary.

Radiative shocks are shocks in a gas where the radiative energy and flux coming from the very hot post-shock material are non-negligible in the shock's total energy budget, and are often large enough to heat the material ahead of the shock. Many simulations of radiative shocks, both in the contexts of astrophysics and laboratory experiments, use a grey treatment of radiative transfer coupled to the hydrodynamics. However, the opacities of the gas show large variations as a function of frequency and this needs to be taken into account if one wishes to reproduce the relevant physics. We have performed radiation hydrodynamics simulations of radiative shocks in Ar using multigroup (frequency dependent) radiative transfer with the HERACLES code. The opacities were taken from the ODALISC database. We show the influence of the number of frequency groups used on the dynamics and morphologies of subcritical and supercritical radiative shocks in Ar gas, and in particular on the extent of the radiative precursor. We find that simulations with even a low number of groups show significant differences compared to single-group (grey) simulations, and that in order to correctly model such shocks, a minimum number of groups is required. Results appear to eventually converge as the number of groups increases above 50. We were also able to resolve in our simulations of supercritical shocks the adaptation zones which connect the cooling layer to the final post-shock state and the precursor. Inside these adaptation zones, we find that the radiative flux just ahead of the shock in one or several high-opacity groups can heat the gas to a temperature higher than the post-shock temperature. Through the use of Hugoniot curves, we have checked the consistency of our radiation hydrodynamics scheme by showing that conservation of mass, momentum and energy (including radiative flux) holds. ABRIDGED.

Radiative transfer plays a major role in the process of star formation. Many simulations of gravitational collapse of a cold gas cloud followed by the formation of a protostellar core use a grey treatment of radiative transfer coupled to the hydrodynamics. However, dust opacities which dominate extinction show large variations as a function of frequency. In this paper, we used frequency-dependent radiative transfer to investigate the influence of the opacity variations on the properties of Larson's first core. We used a multigroup M1 moment model in a 1D radiation hydrodynamics code to simulate the spherically symmetric collapse of a 1 solar mass cloud core. Monochromatic dust opacities for five different temperature ranges were used to compute Planck and Rosseland means inside each frequency group. The results are very consistent with previous studies and only small differences were observed between the grey and multigroup simulations. For a same central density, the multigroup simulations tend to produce first cores with a slightly higher radius and central temperature. We also performed simulations of the collapse of a 10 and 0.1 solar mass cloud, which showed the properties of the first core to be independent of the initial cloud mass, with again no major differences between grey and multigroup models. For Larson's first collapse, where temperatures remain below 2000 K, the vast majority of the radiation energy lies in the IR regime and the system is optically thick. In this regime, the grey approximation does a good job reproducing the correct opacities, as long as there are no large opacity variations on scales much smaller than the width of the Planck function. The multigroup method is however expected to yield more important differences in the later stages of the collapse when high energy (UV and X-ray) radiation is present and matter and radiation are strongly decoupled.

We present the first detailed spatio-kinematical analysis and modelling of the planetary nebula Shapley 1 (Sp 1), which is known to contain a close-binary central star system. Close-binary central stars have been identified as a likely source of shaping in planetary nebulae, but with little observational support to date. Deep narrowband imaging in the light of [O III] {\lambda}5007A suggests the presence of a large bow-shock to the west of the nebula, indicating that it is undergoing the first stages of an interaction with the interstellar medium. Further narrowband imaging in the light of H{\alpha}+[NII] {\lambda}6584A combined with longslit observations of the H{\alpha} emission have been used to develop a spatio-kinematical model of Sp 1. The model clearly reveals Sp 1 to be a bipolar, axisymmetric structure viewed almost pole-on. The symmetry axis of the model nebula is within a few degrees of perpendicular to the orbital plane of the central binary system - strong evidence that the central close-binary system has played an important role in shaping the nebula. Sp 1 is one of very few nebulae to have this link, between nebular symmetry axis and binary plane, shown observationally.

Following the Swift X-ray observations of the 2006 outburst of the recurrent nova RS Ophiuchi, we developed hydrodynamical models of mass ejection from which the forward shock velocities were used to estimate the ejecta mass and velocity. In order to further constrain our model parameters, here we present synthetic X-ray spectra from our hydrodynamical calculations which we compare to the Swift data. An extensive set of simulations was carried out to find a model which best fits the spectra up to 100 days after outburst. We find a good fit at high energies but require additional absorption to match the low energy emission. We estimate the ejecta mass to be in the range (2-5) x 10^{-7} solar masses and the ejection velocity to be greater than 6000 km/s (and probably closer to 10,000 km/s). We also find that estimates of shock velocity derived from gas temperatures via standard model fits to the X-ray spectra are much lower than the true shock velocities.

We present in this paper a multigroup model for radiation hydrodynamics to account for variations of the gas opacity as a function of frequency. The entropy closure model (M1) is applied to multigroup radiation transfer in a radiation hydrodynamics code. In difference from the previous grey model, we are able to reproduce the crucial effects of frequency-variable gas opacities, a situation omnipresent in physics and astrophysics. We also account for the energy exchange between neighbouring groups which is important in flows with strong velocity divergence. These terms were computed using a finite volume method in the frequency domain. The radiative transfer aspect of the method was first tested separately for global consistency (reversion to grey model) and against a well established kinetic model through Marshak wave tests with frequency dependent opacities. Very good agreement between the multigroup M1 and kinetic models was observed in all tests. The successful coupling of the multigroup radiative transfer to the hydrodynamics was then confirmed through a second series of tests. Finally, the model was linked to a database of opacities for a Xe gas in order to simulate realistic multigroup radiative shocks in Xe. The differences with the previous grey models are discussed.

We present the first detailed spatio-kinematical analysis and modelling of the planetary nebula Abell~41, which is known to contain the well-studied close-binary system MT Ser. This object represents an important test case in the study of the evolution of planetary nebulae with binary central stars as current evolutionary theories predict that the binary plane should be aligned perpendicular to the symmetry axis of the nebula. Longslit observations of the \NII\ emission from Abell~41 were obtained using the Manchester Echelle Spectrometer on the 2.1-m San Pedro M\'artir Telescope. These spectra, combined with deep, narrowband imagery acquired using ACAM on the William Herschel Telescope, were used to develop a spatio-kinematical model of \NII\ emission from Abell~41. The best fitting model reveals Abell~41 to have a waisted, bipolar structure with an expansion velocity of $\sim$40\kms{} at the waist. The symmetry axis of the model nebula is within 5$^\circ$ of perpendicular to the orbital plane of the central binary system. This provides strong evidence that the close-binary system, MT Ser, has directly affected the shaping of its host nebula, Abell~41.

We present the first detailed spatio-kinematical analysis and modelling of the planetary nebula Abell 41, which is known to contain the well-studied close-binary system MT Ser. This object represents an important test case in the study of the evolution of planetary nebulae with binary central stars as current evolutionary theories predict that the binary plane should be aligned perpendicular to the symmetry axis of the nebula. Deep narrowband imaging in the light of [NII], [OIII] and [SII], obtained using ACAM on the William Herschel Telescope, has been used to investigate the ionisation structure of Abell 41. Longslit observations of the H-alpha and [NII] emission were obtained using the Manchester Echelle Spectrometer on the 2.1-m San Pedro M\'artir Telescope. These spectra, combined with the narrowband imagery, were used to develop a spatio-kinematical model of [NII] emission from Abell 41. The best fitting model reveals Abell 41 to have a waisted, bipolar structure with an expansion velocity of ~40km\s at the waist. The symmetry axis of the model nebula is within 5$\degr$ of perpendicular to the orbital plane of the central binary system. This provides strong evidence that the close-binary system, MT Ser, has directly affected the shaping of its nebula, Abell 41. Although the theoretical link between bipolar planetary nebulae and binary central stars is long established, this nebula is only the second to have this link, between nebular symmetry axis and binary plane, proved observationally.

We present the first detailed spatio-kinematical analysis and modelling of the Southern planetary nebula SuWt 2. This object presents a problem for current theories of planetary nebula formation and evolution, as it is not known to contain a central post-main sequence star. Deep narrowband [NII]6584 images reveal the presence of faint bipolar lobes emanating fromthe edges of the nebular ring. Longslit observations of the H-alpha and [NII]6584 emission lines were obtained using EMMI on the 3.6-m ESO-NTT. The spectra reveal the nebular morphology as a bright torus encircling the waist of an extended bipolar structure. By deprojection, the inclination of the ring is found to be 68$\degr$ $\pm$ 2$\degr$ (c.f. ~90$\degr$ for the double A-type binary believed to lie at the centre of the nebula), and the ring expansion velocity is found to be 28 km/s. Our findings are discussed with relation to possible formation scenarios for SuWt 2. Through comparison of the nebular heliocentric systemic velocity, found here to be -25 $\pm$ 5 km/s, and the heliocentric systemic velocity of the double A-type binary, we conclude that neither component of the binary could have been the nebular progenitor. However, we are unable to rule out the presence of a third component to the system, which would have been the nebula progenitor.

We present a detailed kinematical analysis of the young compact hourglass-shaped planetary nebula Hb 12. We performed optical imaging and longslit spectroscopy of Hb 12 using the Manchester echelle spectrometer with the 2.1m San Pedro Martir telescope. We reveal, for the first time, the presence of end caps (or knots) aligned with the bipolar lobes of the planetary nebula shell in a deep [NII]6584 image of Hb 12. We measured from our spectroscopy radial velocities of 120 km/s for these knots. We have derived the inclination angle of the hourglass shaped nebular shell to be 65 degrees to the line of sight. It has been suggested that Hb 12's central star system is an eclipsing binary (Hsia et al. 2006) which would imply a binary inclination of at least 80 degrees. However, if the central binary has been the major shaping influence on the nebula then both nebula and binary would be expected to share a common inclination angle. Finally, we report the discovery of high-velocity knots with Hubble-type velocities, close to the core of Hb 12, observed in Halpha and oriented in the same direction as the end caps. Very different velocities and kinematical ages were calculated for the outer and inner knots showing that they may originate from different outburst events.

The ouflowing proper motions of fifteen knots in the dominant northwestern lobe of the high-excitation poly-polar planetary nebula NGC 6302 have been determined by comparing their positions relative to those of faint stars in an image taken at the San Pedro Martir Observatory in 2007 to those in a South African Astronomical Observatory archival plate obtained by Evans in 1956. The Hubble-type expansion of this lobe is now directly confirmed in a model independent way from these measurements. Furthermore, an unambiguous distance to NGC 6302 of 1.17 +/- 0.14 kpc is now determined. Also all the velocity vectors of the fifteen knots (and two others) point back to the central source. An eruptive event from within the central torus, approximately 2200 years previously must have created the high speed lobes of NGC 6302.

High-resolution longslit Halpha spectra of the shell of the old nova DQ Her have been obtained with the William Herschel Telescope using the ISIS spectrograph. An equatorial expansion velocity of 370+/-14 km/s is derived from the spectra which, in conjunction with a narrowband Halpha image of the remnant, allows a distance estimate of 525+/-28 pc. An equatorial ring which exhibits enhanced [NII] emission has also been detected and the inclination angle of the shell is found to be 86.8+/-0.2 degrees with respect to the line of sight. The spectra also reveal tails extending from the clumps in the shell, which have a radial velocity increasing along their length. This suggests the presence of a stellar wind, collimated in the polar direction, which ablates fragments of material from the clumps and accelerates them into its stream up to a terminal velocity of order 800-900 km/s.

Following the early Swift X-ray observations of the latest outburst of the recurrent nova RS Ophiuchi in February 2006 (Paper I), we present new 1D hydrodynamical models of the system which take into account all three phases of the remnant evolution. The models suggest a novel way of modelling the system by treating the outburst as a sudden increase then decrease in wind mass-loss rate and velocity. The differences between this wind model and previous Primakoff-type simulations are described. A more complex structure, even in 1D, is revealed through the presence of both forward and reverse shocks, with a separating contact discontinuity. The effects of radiative cooling are investigated and key outburst parameters such as mass-loss rate, ejecta velocity and mass are varied. The shock velocities as a function of time are compared to the ones derived in Paper I. We show how the manner in which the matter is ejected controls the evolution of the shock and that for a well-cooled remnant, the shock deceleration rate depends on the amount of energy that is radiated away.

We present the first detailed kinematical analysis of the planetary nebula Abell 63, which is known to contain the eclipsing close-binary nucleus UU Sge. Abell 63 provides an important test case in investigating the role of close-binary central stars on the evolution of planetary nebulae. Longslit observations were obtained using the Manchester echelle spectrometer combined with the 2.1-m San Pedro Martir Telescope. The spectra reveal that the central bright rim of Abell 63 has a tube-like structure. A deep image shows collimated lobes extending from the nebula, which are shown to be high-velocity outflows. The kinematic ages of the nebular rim and the extended lobes are calculated to be 8400+/-500 years and 12900+/-2800 years, respectively, which suggests that the lobes were formed at an earlier stage than the nebular rim. This is consistent with expectations that disk-generated jets form immediately after the common envelope phase. A morphological-kinematical model of the central nebula is presented and the best-fit model is found to have the same inclination as the orbital plane of the central binary system; this is the first proof that a close-binary system directly affects the shaping of its nebula. A Hubble-type flow is well-established in the morphological-kinematical modelling of the observed line profiles and imagery. Two possible formation models for the elongated lobes of Abell 63 are considered (1) a low-density, pressure-driven jet excavates a cavity in the remnant AGB envelope; (2) high-density bullets form the lobes in a single ballistic ejection event.