Research into the processes of photoionised nebulae plays a significant part in our understanding of stellar evolution. It is extremely difficult to visually represent or model ionised nebula, requiring astronomers to employ sophisticated modelling code to derive temperature, density and chemical composition. Existing codes are available that often require steep learning curves and produce models derived from mathematical functions. In this article we will introduce PyCross: PyCloudy Rendering Of Shape Software. This is a pseudo 3D modelling application that generates photoionisation models of optically thin nebulae, created using the Shape software. Currently PyCross has been used for novae and planetary nebulae, and it can be extended to Active Galactic Nuclei or any other type of photoionised axisymmetric nebulae. Functionality, an operational overview, and a scientific pipeline will be described with scenarios where PyCross has been adopted for novae (V5668 Sagittarii (2015) & V4362 Sagittarii (1994)) and a planetary nebula (LoTr1). Unlike the aforementioned photoionised codes this application does not require any coding experience, nor the need to derive complex mathematical models, instead utilising the select features from Cloudy/PyCloudy and Shape. The software was developed using a formal software development lifecycle, written in Python and will work without the need to install any development environments or additional python packages. This application, Shape models and PyCross archive examples are freely available to students, academics and research community on GitHub for download (https://github.com/karolfitzgerald/PyCross_OSX_App).

Classical novae are eruptions on the surface of a white dwarf in a binary system. The material ejected from the white dwarf surface generally forms an axisymmetric shell of gas and dust around the system. The three-dimensional structure of these shells is difficult to untangle when viewed on the plane of the sky. In this work a geometrical model is developed to explain new observations of the 2015 nova V5668 Sagittarii. To understand the ionisation structure in terms of the nova shell morphology and estimate the emission distribution directly following the light-curve's dust-dip. High-cadence optical polarimetry and spectroscopy observations of a nova are presented. The ejecta is modelled in terms of morpho-kinematics and photoionisation structure. Initially observational results are presented, including broadband polarimetry and spectroscopy of V5668 Sgr nova during eruption. Variability over these observations provides clues towards the evolving structure of the nova shell. The position angle of the shell is derived from polarimetry, which is attributed to scattering from small dust grains. Shocks in the nova outflow are suggested in the photometry and the effect of these on the nova shell are illustrated with various physical diagnostics. Changes in density and temperature as the super soft source phase of the nova began are discussed. Gas densities are found to be of the order of 10$^{9}$ cm$^{-3}$ for the nova in its auroral phase. The blackbody temperature of the central stellar system is estimated to be around $2.2\times10^{5}$ K at times coincident with the super soft source turn-on. It was found that the blend around 4640 $\rm{\AA}$ commonly called `nitrogen flaring' is more naturally explained as flaring of the O~{\sc ii} multiplet (V1) from 4638 - 4696 $\rm{\AA}$, i.e. `oxygen flaring'.

To gain a more complete understanding of the dynamics of the GK Per (1901) remnant faint-object high-resolution echelle spectroscopic observations and imaging were undertaken covering the knots which comprise the nova shell and the surrounding nebulosity. New imaging from the Aristarchos telescope in Greece and long-slit spectra from the MES instrument at the San Pedro Martir observatory in Mexico were obtained, supplemented with archival observations from several other optical telescopes. Position-velocity arrays are produced of the shell, and also individual knots, and are then used for morpho-kinematic modelling with the shape code. Evidence is found for the interaction of knots with each other and with a wind component, most likely the periodic fast wind emanating from the central binary system. We find that a cylindrical shell with a lower velocity polar structure gives the best model fit to the spectroscopy and imaging. We show in this work that the previously seen jet-like feature is of low velocity. The individual knots have irregular tail shapes; we propose here that they emanate from episodic winds from ongoing dwarf nova outbursts by the central system. The nova shell is cylindrical and the symmetry axis relates to the inclination of the central binary system. Furthermore, the cylinder axis is aligned with the long axis of the bipolar planetary nebula in which it is embedded. Thus, the central binary system is responsible for the bipolarity of the planetary nebula and the cylindrical nova shell. The gradual planetary nebula ejecta versus sudden nova ejecta is the reason for the different degrees of bipolarity. We propose that the 'jet' feature is an illuminated lobe of the fossil planetary nebula that surrounds the nova shell.

Here we explore the possible preexisting circumstellar debris of supernova type Ia systems. Classical, symbiotic and recurrent novae all accrete onto roughly solar mass white dwarfs from main sequence or Mira type companions and result in thermonuclear runaways and expulsion of the accreted material at high velocity. The expelled material forms a fast moving shell that eventually slows to planetary nebula expansion velocities within several hundred years. All such systems are recurrent and thousands of shells (each of about 0.001 Mo) snow plough into the environment. As these systems involve common envelope binaries the material is distributed in a non-spherical shell. These systems could be progenitors of some SN Ia and thus explode into environments with large amounts of accumulated gas and dust distributed in thin non-spherical shells. Such shells should be observable around 100 years after a SN Ia event in a radio flash as the SN Ia debris meets that of the ejected material of the systems previous incarnation.

Molecules with hyperfine splitting of their rotational line spectra are useful probes of optical depth, via the relative line strengths of their hyperfine components.The hyperfine splitting is particularly advantageous in interpreting the physical conditions of the emitting gas because with a second rotational transition, both gas density and temperature can be derived. For HCN however, the relative strengths of the hyperfine lines are anomalous. They appear in ratios which can vary significantly from source to source, and are inconsistent with local thermodynamic equilibrium. This is the HCN hyperfine anomaly, and it prevents the use of simple LTE models of HCN emission to derive reliable optical depths. In this paper we demonstrate how to model HCN hyperfine line emission, and derive accurate line ratios, spectral line shapes and optical depths. We show that by carrying out radiative transfer calculations over each hyperfine level individually, as opposed to summing them over each rotational level, the anomalous hyperfine emission emerges naturally. To do this requires not only accurate radiative rates between hyperfine states, but also accurate collisional rates. We investigate the effects of different sets of hyperfine collisional rates, derived via the 'proportional method' and through direct recoupling calculations. Through an extensive parameter sweep over typical low mass star forming conditions, we show the HCN line ratios to be highly variable to optical depth. We also reproduce an observed effect whereby the red-blue asymmetry of the hyperfine lines (an infall signature) switches sense within a single rotational transition.

Three bright molecular line sources in G333 have recently been shown to exhibit signatures of infall. We describe a molecular line radiative transfer modelling process which is required to extract the infall signature from Mopra and Nanten2 data. The observed line profiles differ greatly between individual sources but are reproduced well by variations upon a common unified model where the outflow viewing angle is the most significant difference between the sources. The models and data together suggest that the observed properties of the high-mass star-forming regions such as infall, turbulence, and mass are consistent with scaled-up versions of the low-mass case with turbulent velocities that are supersonic and an order of magnitude larger than those found in low-mass star-forming regions. Using detailed radiative transfer modeling, we show that the G333 cores are essentially undergoing a scaled-up version of low mass star formation. This is an extension of earlier work in that the degree of infall and the chemical abundances are constrained by the RT modeling in a way that is not practical with a standard analysis of observational data. We also find high velocity infall and high infall mass rates, possibly suggesting accelerated collapse due to external pressure. Molecular depletion due to freeze-out onto dust grains in central regions of the cores is suggested by low molecular abundances of several species. Strong evidence for a local enhancement of 13C-bearing species towards the outflow cloud cores is discussed, consistent with the presence of shocks caused by the supersonic motions within them.

Abell 14 is a poorly studied object despite being considered a born again planetary nebula. We performed a detailed study of its 3D morphology and ionization structure using the SHAPE and MOCASSIN codes. We found that Abell 14 is a highly evolved, bipolar nebula with a kinematical age of $\sim$19,400 yr for a distance of 4 kpc. The high He abundance, and N/O ratio indicate a progenitor of 5 $M_{\odot}$ that has experienced the third dredge-up and hot bottom burning phases. The stellar parameters of the central source reveal a star at a highly evolved stage near to the white dwarf cooling track, being inconsistent with the born again scenario. The nebula shows unexpectedly strong [N I] $\lambda 5200$ and [O I] $\lambda 6300$ emission lines indicating possible shock interactions. Abell 14 appears to be a member of a small group of highly evolved, extreme Type-I PNe. The members of this group lie at the lower-left corner of the PNe regime on the [N II]/H$\alpha$ vs. [S II]/H$\alpha$ diagnostic diagram, where shock--excited regions/objects are also placed. The low luminosity of their central stars, in conjunction with the large physical size of the nebulae, result in a very low photo-ionization rate, which can make any contribution of shock interaction easily perceptible, even for small velocities.

GK Persei (1901, the 'Firework Nebula') is an old but bright nova remnant that offers a chance to probe the physics and kinematics of nova shells. The kinematics in new and archival longslit optical echelle spectra were analysed using the shape software. New imaging from the Aristarchos telescope continues to track the proper motion, extinction and structural evolution of the knots, which have been observed intermittently over several decades. We present for the first time, kinematical constraints on a large faint 'jet' feature, that was previously detected beyond the shell boundary. These observational constraints allow for the generation of models for individual knots, interactions within knot complexes, and the 'jet' feature. Put together, and taking into account dwarf-nova accelerated winds emanating from the central source, these data and models give a deeper insight into the GK Per nova remnant as a whole.

Context. Complex bipolar shapes can be generated either as a planetary nebula or a symbiotic system. The origin of the material ionised by the white dwarf is very different in these two scenarios, and it complicates the understanding of the morphologies of planetary nebulae. Aims. The physical properties, structure, and dynamics of the bipolar nebulae, M 2-9, Mz 3, and Hen 2-104, are investigated in detail with the aim of understanding their nature, shaping mechanisms, and evolutionary history. Methods. Long-slit optical echelle spectra are used to investigate the morpho-kinematics of M 2-9, Mz 3, and Hen 2-104. Near-infrared (NIR) data, as well as optical, spectra are used to separate Galactic symbiotic-type nebulae from genuine planetary nebulae by means of a 2MASS J-H/H-Ks diagram and a {\lambda}4363/H{\gamma} vs. {\lambda}5007/H\b{eta} diagnostic diagram, respectively. Results. The best-fitted 3-D models for M 2-9, Mz 3, and Hen 2-104 provide invaluable kinematical information on the expansion velocity of its nebular components by means of synthetic spectra. Kinematical ages of the different structures of M 2-9 and Mz 3 have also been determined. Both diagnostic diagrams show M 2-9 and Hen 2-104 to fall well within the category of having a symbiotic source, whereas Mz 3 borders the region of symbiotic and young planetary nebulae in the optical diagram. The optical diagnostic diagram is shown to successfully separate the two types of nebulae. Conclusions. The morphology, kinematics, and evolutionary history of M 2-9, Mz 3, and Hen 2-104 are better understood using the interactive 3-D modelling tool shape. The optical and NIR diagnostic diagrams used are important techniques for separating Galactic symbiotic-type nebulae from genuine planetary nebulae.

We present molecular line imaging observations of three massive molecular outflow sources, G333.6-0.2, G333.1-0.4, and G332.8-0.5, all of which also show evidence for infall, within the G333 giant molecular cloud (GMC). All three are within a beam size (36 arcseconds) of IRAS sources, 1.2-mm dust clumps, various masing species and radio continuum-detected HII regions and hence are associated with high-mass star formation. We present the molecular line data and derive the physical properties of the outflows including the mass, kinematics, and energetics and discuss the inferred characteristics of their driving sources. Outflow masses are of 10 to 40 solar masses in each lobe, with core masses of order 10^3 solar masses. outflow size scales are a few tenth of a parsec, timescales are of several x10^4 years, mass loss rates a few x10^-4 solar masses/year. We also find the cores are turbulent and highly supersonic.

During the late stages of their evolution, Sun-like stars bring the products of nuclear burning to the surface. Most of the carbon in the Universe is believed to originate from stars with masses up to a few solar masses. Although there is a chemical dichotomy between oxygen-rich and carbon-rich evolved stars, the dredge-up itself has never been directly observed. In the last three decades, however, a few stars have been shown to display both carbon- and oxygen-rich material in their circumstellar envelopes. Two models have been proposed to explain this dual chemistry: one postulates that a recent dredge-up of carbon produced by nucleosynthesis inside the star during the Asymptotic Giant Branch changed the surface chemistry of the star. The other model postulates that oxygen-rich material exists in stable keplerian rotation around the central star. The two models make contradictory, testable, predictions on the location of the oxygen-rich material, either located further from the star than the carbon-rich gas, or very close to the star in a stable disk. Using the Faint Object InfraRed CAmera (FORCAST) instrument on board the Stratospheric Observatory for Infrared Astronomy (SOFIA) Telescope, we obtained images of the carbon-rich planetary nebula (PN) BD+30 3639 which trace both carbon-rich polycyclic aromatic hydrocarbons (PAHs) and oxygen-rich silicate dust. With the superior spectral coverage of SOFIA, and using a 3D photoionisation and dust radiative transfer model we prove that the O-rich material is distributed in a shell in the outer parts of the nebula, while the C-rich material is located in the inner parts of the nebula. These observations combined with the model, suggest a recent change in stellar surface composition for the double chemistry in this object. This is evidence for dredge-up occurring ~1000yr ago.

We use the Mopra radio telescope to test for expansion of the molecular gas associated with the bubble HII region RCW120. A ring, or bubble, morphology is common for Galactic HII regions, but the three-dimensional geometry of such objects is still unclear. Detected near- and far-side expansion of the associated molecular material would be consistent with a three-dimensional spherical object. We map the $J = 1\rightarrow 0$ transitions of $^{12}$CO, $^{13}$CO, C$^{18}$O, and C$^{17}$O, and detect emission from all isotopologues. We do not detect the $0_0\rightarrow 1_{-1} E$ masing lines of CH$_3$OH at 108.8939 GHz. The strongest CO emission is from the photodissociation region (PDR), and there is a deficit of emission toward the bubble interior. We find no evidence for expansion of the molecular material associated with RCW120 and therefore can make no claims about its geometry. The lack of detected expansion is roughly in agreement with models for the time-evolution of an HII region like RCW120, and is consistent with an expansion speed of $< 1.5\, {\rm km\, s^{-1}}$. Single-position CO spectra show signatures of expansion, which underscores the importance of mapped spectra for such work. Dust temperature enhancements outside the PDR of RCW120 coincide with a deficit of emission in CO, confirming that these temperature enhancements are due to holes in the RCW120 PDR. H$\alpha$ emission shows that RCW120 is leaking $\sim5\%$ of the ionizing photons into the interstellar medium (ISM) through PDR holes at the locations of the temperature enhancements. H-alpha emission also shows a diffuse "halo" from leaked photons not associated with discrete holes in the PDR. Overall $25\pm10\%$ of all ionizing photons are leaking into the nearby ISM.

Context. The mechanisms that form extreme bipolar planetary nebulae remain unclear. Aims. The physical properties, structure, and dynamics of the bipolar planetary nebula, MyCn 18, are investigated in detail with the aim of understanding the shaping mechanism and evolutionary history of this object. Methods. VLT infrared images, VLT ISAAC infrared spectra, and long-slit optical Echelle spectra are used to investigate MyCn 18. Morpho-kinematic modelling was used to firmly constrain the structure and kinematics of the source. A timescale analysis was used to determine the kinematical age of the nebula and its main components. Results. A spectroscopic study of MyCn 18's central and offset region reveals the detailed make-up of its nebular composition. Molecular hydrogen, atomic helium, and Bracket gamma emission are detected from the central regions of MyCn 18. ISAAC spectra from a slit position along the narrow waist of the nebula demonstrate that the ionised gas resides closer to the centre of the nebula than the molecular emission. A kinematical age of the nebula and its components were obtained by the P-V arrays and timescale analysis. Conclusions. The structure and kinematics of MyCn 18 are better understood using an interactive 3-D modelling tool called shape. A dimensional and timescale analysis of MyCn 18's major components provides a possible mechanism for the nebula's asymmetry. The putative central star is somewhat offset from the geometric centre of the nebula, which is thought to be the result of a binary system. We speculate that the engulfing and destruction of an exoplanet during the AGB phase may have been a key event in shaping MyCn 18 and generating of its hypersonic knotty outflow.

A rich variety of molecular species has now been observed towards hot cores in star forming regions and in the interstellar medium. An increasing body of evidence from millimetre interferometers suggests that many of these form at the interfaces between protostellar outflows and their natal molecular clouds. However, current models have remained unable to explain the origin of the observational bias towards wide-angled "brandy snifter" shaped outflows over narrower "champagne flute" shapes in carbon monoxide imaging. Furthermore, these wide-angled systems exhibit unusually high abundances of the molecular ion HCO$^+$. We present results from a chemo-dynamic model of such regions where a rich chemistry arises naturally as a result of turbulent mixing between cold, dense molecular gas and the hot, ionized outflow material. The injecta drives a rich and rapid ion-neutral chemistry in qualitative and quantitative agreement with the observations. The observational bias towards wide-angled outflows is explained naturally by the geometry-dependent ion injection rate causing rapid dissociation of CO in the younger systems.

Low mass star-forming regions are more complex than the simple spherically symmetric approximation that is often assumed. We apply a more realistic infall/outflow physical model to molecular/continuum observations of three late Class 0 protostellar sources with the aims of (a) proving the applicability of a single physical model for all three sources, and (b) deriving physical parameters for the molecular gas component in each of the sources. We have observed several molecular species in multiple rotational transitions. The observed line profiles were modelled in the context of a dynamical model which incorporates infall and bipolar outflows, using a three dimensional radiative transfer code. This results in constraints on the physical parameters and chemical abundances in each source. Self-consistent fits to each source are obtained. We constrain the characteristics of the molecular gas in the envelopes as well as in the molecular outflows. We find that the molecular gas abundances in the infalling envelope are reduced, presumably due to freeze-out, whilst the abundances in the molecular outflows are enhanced, presumably due to dynamical activity. Despite the fact that the line profiles show significant source-to-source variation, which primarily derives from variations in the outflow viewing angle, the physical parameters of the gas are found to be similar in each core.

HCN is becoming a popular choice of molecule for studying star formation in both low- and high-mass regions and for other astrophysical sources from comets to high-redshift galaxies. However, a major and often overlooked difficulty with HCN is that it can exhibit non-local thermodynamic equilibrium (non-LTE) behaviour in its hyperfine line structure. Individual hyperfine lines can be strongly boosted or suppressed. In low-mass star-forming cloud observations, this could possibly lead to large errors in the calculation of opacity and excitation temperature, while in massive star-forming clouds, where the hyperfine lines are blended due to turbulent broadening, errors will arise in infall measurements that are based on the separation of the peaks in a self-absorbed profile. The underlying line shape cannot be known for certain if hyperfine anomalies are present. We present a first observational investigation of these anomalies across a range of conditions and transitions by carrying out a survey of low-mass starless cores (in Taurus & Ophiuchus) and high-mass protostellar objects (in the G333 giant molecular cloud) using hydrogen cyanide (HCN) J=1-0 and J=3-2 emission lines. We quantify the degree of anomaly in these two rotational levels by considering ratios of individual hyperfine lines compared to LTE values. We find that all the cores observed show some degree of anomaly while many of the lines are severely anomalous. We conclude that HCN hyperfine anomalies are common in both lines in both low-mass and high-mass protostellar objects, and we discuss the differing hypotheses for the generation of the anomalies. In light of the results, we favour a line overlap effect for the origins of the anomalies. We discuss the implications for the use of HCN as a dynamical tracer and suggest in particular that the J=1-0, F=0-1 hyperfine line should be avoided in quantitative calculations.

Although hydrogen cyanide has become quite a common molecular tracing species for a variety of astrophysical sources, it, however, exhibits dramatic non-LTE behaviour in its hyperfine line structure. Individual hyperfine components can be strongly boosted or suppressed. If these so-called hyperfine line anomalies are present in the HCN rotational spectra towards low or high mass cores, this will affect the interpretation of various physical properties such as the line opacity and excitation temperature in the case of low mass objects and infall velocities in the case of their higher mass counterparts. This is as a consequence of the direct effects that anomalies have on the underlying line shape, be it with the line structural width or through the inferred line strength. This work involves the first observational investigation of these anomalies in two HCN rotational transitions, J=1!0 and J=3!2, towards both low mass starless cores and high mass protostellar objects. The degree of anomaly in these two rotational transitions is considered by computing the ratios of neighboring hyperfine lines in individual spectra. Results indicate some degree of anomaly is present in all cores considered in our survey, the most likely cause being line overlap effects among hyperfine components in higher rotational transitions.

Cold massive cores are one of the earliest manifestations of high mass star formation. Following the detection of SiO emission from G333.125-0.562, a cold massive core, further investigations of the physics, chemistry and dynamics of this object has been carried out. Mopra and NANTEN2 molecular line profile observations, Australia Telescope Compact Array (ATCA) line and continuum emission maps, and Spitzer 24 and 70 \mum images were obtained. These new data further constrain the properties of this prime example of the very early stages of high mass star formation. A model for the source was constructed and compared directly with the molecular line data using a 3D molecular line transfer code - MOLLIE. The ATCA data reveal that G333.125-0.562 is composed of two sources. One of the sources is responsible for the previously detected molecular outflow and is detected in the Spitzer 24 and 70 \mum band data. Turbulent velocity widths are lower than other more active regions of G333 which reflects the younger evolutionary stage and/or lower mass of this core. The molecular line modelling requires abundances of the CO isotopes that strongly imply heavy depletion due to freeze-out of this species onto dust grains. The principal cloud is cold, moderately turbulent and possesses an outflow which indicates the presence of a central driving source. The secondary source could be an even less evolved object as no apparent associations with continuum emissions at (far-)infrared wavelengths.

The Z CMa binary is understood to undergo both FU Orionis (FUOR) and EX Orionis (EXOR) type outbursts. While the SE component has been spectro- scopically identified as an FUOR, the NW component, a Herbig Be star, is the source of the EXOR outbursts. The system has been identified as the source of a large outflow, however, previous studies have failed to identify the driver. Here we present adaptive optics (AO) assisted [FeII] spectro-images which reveal for the first time the presence of two jets. Observations made using OSIRIS at the Keck Observatory show the Herbig Be star to be the source of the parsec-scale outflow, which within 2'' of the source shows signs of wiggling and the FUOR to '' be driving a ~ 0.4 jet. The wiggling of the Herbig Be star's jet is evidence for an additional companion which could in fact be generating the EXOR outbursts, the last of which began in 2008 (Grankin & Artemenko 2009). Indeed the dy- namical scale of the wiggling corresponds to a time-scale of 4-8 years which is in agreement with the time-scale of these outbursts. The spectro-images also show a bow-shock shaped feature and possible associated knots. The origin of this structure is as of yet unclear. Finally interesting low velocity structure is also observed. One possibility is that it originates in a wide-angle outflow launched from a circumbinary disk.

The unique Honeycomb nebula, most likely a peculiar supernova remnant, lies in 30 Doradus in the Large Magellanic Cloud. Due to its proximity to SN1987A, it has been serendipitously and intentionally observed at many wavelengths. Here, an optical spectral analysis of forbidden line ratios is performed in order to compare the Honeycomb high-speed gas with supernova remnants in the Galaxy and the LMC, with galactic Wolf-Rayet nebulae and with the optical line emission from the interaction zone of the SS433 microquasar and W50 supernova remnant system. An empirical spatiokinematic model of the images and spectra for the Honeycomb reveals that its striking appearance is most likely due to a fortuitous viewing angle. The Honeycomb nebula is more extended in soft X-ray emission and could in fact be a small part of the edge of a giant LMC shell revealed for the first time in this short wavelength domain. It is also suggested that a previously unnoticed region of optical emission may in fact be an extension of the Honeycomb around the edge of this giant shell. A secondary supernova explosion in the edge of a giant shell is considered for the creation of the Honeycomb nebula. A microquasar origin of the Honeycomb nebula as opposed to a simple supernova origin is also evaluated.

An example of a cold massive core, JCMT 18354-0649S, a possible high mass analogue to a low mass star forming core is studied. Line and continuum observations from JCMT, Mopra Telescope and Spitzer are presented and modelled in detail using a 3D molecular line radiative transfer code. In almost every way JCMT 18354-0649S is a scaled-up version of a typical low mass core with similar temperatures, chemical abundances and densities. The difference is that both the infall velocity and the turbulent width of the line profiles are an order of magnitude larger. While the higher infall velocity is expected due to the large mass of JCMT 18354-0649S, we suggest that the dissipation of this highly supersonic turbulence may lead to the creation of dense clumps of gas that surround the high mass core.

Optical images and high-dispersion spectra have been obtained of the ejected material surrounding the pulsating AGB star Mira A. The two streams of knots on either side of the star, found in far ultra-viollet (FUV) GALEX images, have now been imaged clearly in the light of Halpha. Spatially resolved profiles of the same line reveal that the bulk of these knots form a bi-polar outflow with radial velocity extremes of +- 150 km/s with respect to the central star. The South stream is approaching and the North stream receding from the observer. A displacement away from Mira A between the position of one of the South stream knots in the new Halpha image and its position in the previous Palomar Observatory Sky Survey (POSS I) red plate has been noted. If interpreted as a consequence of expansion proper motions the bipolar outflow is tilted at 69deg +- 2deg to the plane of the sky, has an outflow velocity of 160 +- 10 km/s and is ~1000 y old.

CO isotopes are able to probe the different components in protostellar clouds. These components, core, envelope and outflow have distinct physical conditions and sometimes more than one component contributes to the observed line profile. In this study we determine how CO isotope abundances are altered by the physical conditions in the different components. We use a 3D molecular line transport code to simulate the emission of four CO isotopomers, 12CO J=2-1, 13CO J=2-1, C18O J=2-1 and C17O J=2-1 from the Class 0/1 object L483, which contains a cold quiescent core, an infalling envelope and a clear outflow. Our models replicate JCMT (James Clerk Maxwell Telescope) line observations with the inclusion of freeze-out, a density profile and infall. Our model profiles of 12CO and 13CO have a large linewidth due to a high velocity jet. These profiles replicate the process of more abundant material being susceptible to a jet. C18O and C17O do not display such a large linewidth as they trace denser quiescent material deep in the cloud.

Knots are commonly found in nearby planetary nebulae (PNe) and star forming regions. Within PNe, knots are often found to be associated with the brightest parts of the nebulae and understanding the physics involved in knots may reveal the processes dominating in PNe. As one of the closest PNe, the Helix Nebula (NGC 7293) is an ideal target to study such small-scale (~300 AU) structures. We have obtained infrared integral spectroscopy of a comet-shaped knot in the Helix Nebula using SINFONI on the Very Large Telescope at high spatial resolution (50-125 mas). With spatially resolved 2 micron spectra, we find that the H2 rotational temperature within the cometary knots is uniform. The rotational-vibrational temperature of the cometary knot (situated in the innermost region of the nebula, 2.5 arcmin away from the central star), is 1800 K, higher than the temperature seen in the outer regions (5-6 arcmin from the central star) of the nebula (900 K), showing that the excitation temperature varies across the nebula. The obtained intensities are reasonably well fitted with 27 km s-1 C-type shock model. This ambient gas velocity is slightly higher than the observed [HeII] wind velocity of 13 km s-1. The gas excitation can also be reproduced with a PDR (photo dominant region) model, but this requires an order of magnitude higher UV radiation. Both models have limitations, highlighting the need for models that treats both hydrodynamical physics and the PDR.

We have carried out a systematic search for the molecular ion CO+ in a sample of 8 protoplanetary and planetary nebulae in order to determine the origin of the unexpectedly strong HCO+ emission previously detected in these sources. An understanding of the HCO+ chemistry may provide direct clues to the physical and chemical evolution of planetary nebulae. We find that the integrated intensity of the CO+ line may be correlated with that of HCO+, suggesting that the reaction of CO+ with molecular hydrogen may be an important formation route for HCO+ in these planetary nebulae.

Four mosaics of deep, continuum-subtracted, CCD images have been obtained over the extensive galactic radio continuum shell, W 50, which surrounds the remarkable stellar system SS 433. Two of these mosaics in the Halpha+[N II] and [O III] 5007 A emission lines respectively cover a field of ~2.3 x 2.5 degr^2 which contains all of W 50 but at a low angular resolution of 5 arcsec. The third and fourth mosaics cover the eastern (in [O III] 5007 A) and western (in Halpha+[N II]) filamentary nebulosity respectively but at an angular resolution of 1 arcsec. These observations are supplemented by new low dispersion spectra and longslit, spatially resolved echelle spectra. The [O III] 5007 A images show for the first time the distribution of this emission in both the eastern and western filaments while new Halpha+[N II] emission features are also found in both of these regions. Approaching flows of faintly emitting material from the bright eastern filaments of up 100 km/s in radial velocity are detected. The present observations also suggest that the heliocentric systemic radial velocity of the whole system is 56+-2 km/s. Furthermore, very deep imagery and high resolution spectroscopy of a small part of the northern radio ridge of W 50 has revealed for the first time the very faint optical nebulosity associated with this edge. It is suggested that patchy foreground dust along the ~5 kpc sightline is inhibiting the detection of all of the optical nebulosity associated with W 50. The interaction of the microquasar jets of SS 433 with the W 50 shell is discussed.

A mosaic of six, deep, CCD images in the light of the [O III] 5007 A nebular emission line has been obtained with the 1.3-m Skinakas (Crete) telescope of the filamentary nebulosity surrounding P Cygni. The [O III] 5007 A line discriminates against confusing galactic H II regions along the same sight-lines and the new mosaic did not include the 4.8 mag. central star; a source of artifacts in the previous lower angular resolution observations. New giant `lobes' and `shells' are found to be clustered around P Cygni which must be the relics of historic eruptions between 2400 and up to ~10^5 yr ago.

New molecular line observations of the Bok globule Barnard 68 in HCO+ irrefutably confirm the complex pattern of red and blue asymmetric line profiles seen across the face of the cloud in previous observations of CS. The new observations thus strengthen the previous interpretation that Barnard 68 is undergoing peculiar oscillations. Furthermore, the physical chemistry of B68 indicates that the object is much older than the sound crossing time and is therefore long-lived. A model is presented for the globule in which a modest external pressure perturbation is shown to lead to oscillations about a stable equilibrium configuration. Such oscillations may be present in other stable starless cores as manifested by a similar signature of inward and outward motions.

Optical CCD imaging and spectroscopic observations of three supernova remnants are presented. Optical emission from G 54.4-0.3 and G 59.8+1.2 is detected for the first time, while the first flux calibrated CCD images of the supernova remnant G 126.2+1.6 were performed in the optical emission lines of Halpha+[N II], [O III] and [S II]. A mixture of filamentary and diffuse structures is observed in G 54.4-0.3 and G 59.8+1.2, mainly in Halpha+[N II], while the deep optical images of G 126.2+1.6 reveal several new filamentary and diffuse structures inside the extent of the remnant as defined by its known radio emission. In all cases, the radio emission is found to be well correlated with the optical filaments. [O III] emission was not detected at G 54.4-0.3 and G 59.8+1.2 while in G 126.2+1.6, significant morphological differences between the low and medium ionization images are present suggesting incomplete shock structures. Deep long-slit spectra were taken at different positions of the remnants. Both the flux calibrated images and the long-slit spectra clearly show that the emission originates from shock-heated gas, while some spectra of G 126.2+1.6 are characterized by large [O III]/Hbeta ratios. This remnant's [O III] flux suggests shock velocities into the interstellar "clouds" between 100 and 120 km/s, while the [O III] absence in the other two remnants indicates slower shock velocities. For all remnants, the [S II]6716/6731 ratio indicates electron densities below 600 cm^{-3} with particularly low densities for G 54.4-0.3 (below 50 cm^{-3}). Finally, the Halpha emission has been measured to be between 3.0 to 15.2x10^{-17} erg/s/cm^2/ arcsec^2, 3.2x10^{-17} erg/s/cm^2/ arcsec^2 and between 6.5 to 16.8x10^{-17} erg/s/cm^2/ arcsec^2 for G 54.4-0.3, G 59.8+1.2 and G 126.2+1.6, respectively.

A detailed analysis is presented of the dominant ionised knot in the halo of the planetary nebula NGC 6543. Observations were made at high spectral and spatial resolution of the [OIII] 5007 line using the Manchester echelle spectrometer combined with the 2.1-m San Pedro Martir Telescope. A 20-element multislit was stepped across the field to give almost complete spatial coverage of the knot and surrounding halo. The spectra reveal, for the first time, gas flows around the kinematically inert knot. The gas flows are found to have velocities comparable to the sound speed as gas is photo-evaporated off an ionised surface. No evidence is found of fast wind interaction with the knot and we find it likely that the fast wind is still contained in a pressure-driven bubble in the core of the nebula. This rules out the possibility of the knot having its origin in instabilities at the interface of the fast and AGB winds. We suggest that the knot is embedded in the slowly expanding Red Giant wind and that its surfaces are being continually photoionised by the central star.

A deep, continuum-subtracted, image of NGC 7293 has been obtained in the light of the Halpha+[N II] emission lines. New images of two filamentary halo stuctures have been obtained and the possible detection of a collimated outflow made. Spatially resolved, longslit profiles of the Halpha+[N II] lines have been observed across several of these features with the MES combined with the SPM 2.1m telescope; these are compared with the [N II]6584, [O III]5007, HeII 6560 and Halpha profiles obtained over the nebular core. The central HeII emission is originating in a ~0.34pc diameter spherical volume expanding at <=12km/s which is surrounded, and partially coincident with an [O III] emitting inner shell expanding at 12km/s. The bright helical structure surrounding this inner region is modelled as a bi-polar nebula with lobe expansions of 25km/s whose axis is tilted at 37deg to the sight line but with a toroidal waist itself expanding at 14 km/s. These observations are compared with the expectations of the interacting two winds model for the formation of PNe. Only after the fast wind has switched off could this global velocity structure be generated. Ablated flows must complicate any interpretation. It is suggested that the clumpy nature of much of the material could play a part in creating the radial `spokes' shown here to be apparently present close to the central star. These `spokes' could in fact be the persistant tails of cometary globules whose heads have now photo-evaporated completely. A halo arc projecting from the north-east of the bright core has a conterpart to the south-east. Anomolies in the position-velocity arrays of line profiles could suggest that these are part of an expanding disc not aligned with the central helical structure though expanding bi-polar lobes along a tilted axis are not ruled out.

The search for the onset of star formation in pre-stellar cores has focussed on the identification of an infall signature in the molecular line profiles of tracer species. The classic infall signature is a double peaked line profile with an asymmetry in the strength of the peaks such that the blue peak is stronger. L1689B is a pre-stellar core and infall candidate but new JCMT HCO+ line profile data, presented here, confirms that both blue and red asymmetric line profiles are present in this source. Moreover, a dividing line can be drawn between the locations where each type of profile is found. It is argued that it is unlikely that the line profiles can be interpreted with simple models of infall or outflow and that rotation of the inner regions is the most likely explanation. A rotational model is developed in detail with a new 3D molecular line transport code and it is found that the best type of model is one in which the rotational velocity profile is in between solid body and Keplerian. It is firstly shown that red and blue asymmetric line profiles can be generated with a rotation model entirely in the absence of any infall motion. The model is then quantitively compared with the JCMT data and an iteration over a range of parameters is performed to minmize the difference between the data and model. The results indicate that rotation can dominate the line profile shape even before the onset of infall.

Further deep, narrow-band images in the light of [O III] 5007 A have been added to the previous mosaic of the faint galactic supernova remnant G 65.3+5.7. Additionally longslit spatially resolved [O III] 5007 A line profiles have been obtained at sample positions using the Manchester Echelle Spectrometer at the San Pedro Martir observatory. The remnant is shown to be predominantly bi-lobal with an EW axis for this structure. However, a faint additional northern lobe has now been revealed. Splitting of the profiles along the slit lengths, when extrapolated to the remnant's centre, although uncertain suggests that the expansion velocity of this remnant is between 124 and 187 km/s ie much lower than the 400 km/s previously predicted for the forward shock velocity from the X-ray emission. An expansion proper motion measurement of 2.1+-0.4 arcsec in 48 years for the remnant's filamentary edge in the light of Halpha+[N II] has also been made. When combined with an expansion velocity of ~155 km/s, a distance of ~800 pc to G 65.3+5.7 is derived. Several possibilities are considered for the large difference in the expansion velocity measured here and the 400 km/s shock velocity required to generate the X-ray emission. It is also suggested that the morphology of the remnant may be created by a tilt in the galactic magnetic field in this vicinity.

Deep, wide-field, continuum-subtracted, images in the light of the Halpha+[NII] 6548 & 6584 A and [O III] 5007 A nebular emission lines have been obtained of the environment of the Luminous Blue Variable (LBV) star P Cygni. A previously discovered, receding, nebulous filament along PA 50 deg has now been shown to extend up to 12' from this star. Furthermore, in the light of [O III] 5007 A, a southern counterpart is discovered as well as irregular filaments on the opposite side of P Cygni. Line profiles from this nebulous complex indicate that this extended nebulosity is similar to that associated with middle-aged supernova remnants. However, there are several indications that it has originated in P Cygni and is not just a chance superposition along the same sight-line. This possibility is explored here and comparison is made with a new image of the LBV star R 143 in the LMC from which similar filaments appear to project. The dynamical age of the P Cygni giant lobe of ~5x10^{4} yr is consistent with both the predicted and observed durations of the LBV phases of 50M stars after they have left the main sequence. Its irregular shape may have been determined by the cavity formed in the ambient gas by the energetic wind of the star, and shaped by a dense torus, when on the main sequence. The proper motion and radial velocity of P Cygni, with respect to its local environment, could explain the observed angular and kinematical shifts of the star compared with the giant lobe.

In cold molecular clouds submillimetre emission lines are excited by the ambient radiation field. The pumping is dominated by the cosmic microwave background (CMB). It is usual in molecular line radiative transfer modelling to simply assume that this is the only incident radiation field. In this paper, a molecular line transport code and a dust radiative transfer code are used to explore the effects of the inclusion of a full interstellar radiation field (ISRF) on a simple test molecular cloud. It is found that in many galactic situations, the shape and strength of the line profiles that result are robust to variations in the ISRF and thus that in most cases, it is safe to adopt the CMB radiation field for the molecular line transport calculations. However, we show that in two examples, the inclusion of a plausible radiation field can have a significant effect on the the line profiles. Firstly, in the vicinity of an embedded massive star, there will be an enhanced far infared component to the radiation field. Secondly, for molecular clouds at large redshift, the CMB temperature increases and this of course also alters the radiation field. In both of these cases, the line profiles are weakened significantly compared to a cloud exposed to a standard radiation field. Therefore this effect should be accounted for when investigating prestellar cores in massive star forming regions and when searching for molecular clouds at high redshift.

A plausible model is proposed for the enhancement of the abundance of molecular species in bipolar outflow sources. In this model, levels of HCO+ enhancement are considered based on previous chemical calculations, that are assumed to result from shock-induced desorption and photoprocessing of dust grain ice mantles in the boundary layer between the outflow jet and the surrounding envelope. A radiative transfer simulation that incorporates chemical variations within the flow shows that the proposed abundance enhancements in the boundary layer are capable of reproducing the observed characteristics of the outflow seen in HCO+ emission in the star forming core L1527. The radiative transfer simulation also shows that the emission lines from the enhanced molecular species that trace the boundary layer of the outflow exhibit complex line profiles indicating that detailed spatial maps of the line profiles are essential in any attempt to identify the kinematics of potential infall/outflow sources. This study is one of the first applications of a full three dimensional radiative transfer code which incorporates chemical variations within the source.

Echelle spectroscopy has been carried out towards a sample region of the halo of the giant HII region 30 Doradus in the Large Magellanic Cloud. This new kinematical data is the amongst the most sensitive yet obtained for this nebula and reveals a wealth of faint, complex high speed features. These are interpreted in terms of localised shells due to individual stellar winds and supernova explosions, and collections of discrete knots of emission that still retain the velocity pattern of the giant shells from which they fragmented. The high speed velocity features may trace the base of the superwind that emanates from the 30 Doradus starburst, distributed around the super star cluster R136.

We study the chemistry in the slow wind during the transition from AGB star to pre-planetary nebula (PPN) and planetary nebula (PN). We show that there is a very rich chemistry of degradation products created by photoprocessing, driven by the gradually hardening radiation field of the central star. Most of these products are, however, photodissociated during the PPN phase if the wind is smooth. By contrast, if the wind is clumpy, possibly because of clumpiness in the AGB atmosphere, then many of these degradation products survive into the PN phase. Thus, chemistry may be used to infer the existence of clumpiness in the AGB phase. We identify potential molecular tracers, and we note that, in the case of clumpiness, large molecules may survive the transport from the stellar atmosphere to the interstellar medium. We compare between our model results with observations of three objects at differing evolutionary stages: CRL618, NGC 7027 and the Helix nebula (NGC 7293).

C17O (J=2-1) observations have been carried out towards the pre-stellar core L1689B. By comparing the relative strengths of the hyperfine components of this line, the emission is shown to be optically thin. This allows accurate CO column densities to be determined and, for reference, this calculation is described in detail. The hydrogen column densities that these measurements imply are substantially smaller than those calculated from SCUBA dust emission data. Furthermore, the C17O column densities are approximately constant across L1689B whereas the SCUBA column densities are peaked towards the centre. The most likely explanation is that CO is depleted from the central regions of L1689B. Simple models of pre-stellar cores with an inner depleted region are compared with the results. This enables the magnitude of the CO depletion to be quantified and also allows the spatial extent of the freeze-out to be firmly established. We estimate that within about 5000 AU of the centre of L1689B, over 90% of the CO has frozen onto grains. This level of depletion can only be achieved after a duration that is at least comparable to the free-fall timescale.

The association between the Pencil nebula (RCW 37, NGC 2736), the Vela X-ray fragment D/D' and the recently discovered new X-ray supernova remnant (RX J0852.0-4622) in Vela is investigated. Recently published Chandra images of D/D' are compared with optical images of RCW 37 and confirm the close association of the two objects. New optical line profiles of RCW 37 from an extended slit position passing through this unusual optical nebula are presented. They reveal a partial velocity ellipse with expansion velocities of around 120 km/s. Various scenarios for the origin of the nebula are considered and the evidence of a link with RX J0852.0-4622 is reviewed. A funnel of gas similar to those in the Crab and DEM34a SNRs is not ruled out but a more plausible explanation may be that a `wavy sheet' is responsible. We suggest that RX J0852.0-4622 is located within the older larger Vela SNR and that some of the X-ray gas from RX J0852.0-4622 has collided with the dense HI wall of the older remnant. This gives rise to the morphology and velocity structure of the optical emission and explains the unusual X-ray emission from this portion of the supernova remnant. If our hypothesis is correct, a distance prediction of 250 +/- 30 pc can be made, based on recent measurements of the distance to the old Vela SNR. This is at the lower end of the range of distances quoted in the literature would confirm unusual nature of this young nearby supernova remnant.

A spatially unresolved velocity feature, with an approaching radial velocity of ~100 km/s with respect to the systemic radial velocity, in a position-velocity array of [O III] 5007 line profiles is identified as the kinematical counterpart of a jet from the proplyd LV 5 (158-323) in the core of the Orion Nebula. The only candidate in HST imagery for this jet appears to be a displaced, ionized knot. Also an elongated jet projects from the proplyd GMR 15 (161-307). Its receding radial velocity difference appears at ~80 km/s in the same position-velocity array. A `standard' model for jets from young, low mass stars invokes an accelerating, continuous flow outwards with an opening angle of a few degrees. Here an alternative explanation is suggested which may apply to some, if not all, of the proplyd jets. In this, a `bullet' of dense material is ejected which ploughs through dense circumstellar ambient gas. The decelerating tail of material ablated from the bullet's surface would be indistinguishable from a continuously emitted jet in current observations.

The narrow optical filaments (`strings' or `spikes') emerging from the Homunculus of Eta Carinae are modelled as resulting from the passage of ballistic `bullets' of material through the dense circumstellar environment. In this explanation, the string is the decelerating flow of ablated gas from the bullet. An archive HST image and new forbidden line profiles of the most distinct of the strings are presented and discussed in terms of this simple model.

A deep \ha image of interlocking filamentary arcs of nebulosity has been obtained with a wide-field ($\approx$ 30\degree diameter) narrow-band filter camera combined with a CCD as a detector. The resultant mosaic of images, extending to a galactic latitude of 65$^{o}$, has been corrected for field distortions and had galactic coordinates superimposed on it to permit accurate correlations with the most recent H{\sc i} (21 cm), X-ray (0.75 kev) and FIR (IRAS 100 $\mu$m) maps. Furthermore, an upper limit of 0.13 arcsec/yr to the expansion proper motion of the primary 25\degree long nebulous arc has been obtained by comparing a recent \ha image obtained with the San Pedro Martir telescope of its filamentary edge with that on a POSS E plate obtained in 1951. It is concluded that these filamentary arcs are the superimposed images of separate shells (driven by supernova explosions and/or stellar winds) rather than the edges of a single `superbubble' stretching from Barnard's Arc (and the Orion Nebula) to these high galactic latitudes. The proper motion measurement argues against the primary \ha emitting arc being associated with the giant radio loop (Loop 2) except in extraordinary circumstances.

We have modelled ultracompact HII regions (UCHIIR) in terms of steady subsonic ionized flows in a clumpy medium. Mass loss from neutral clumps allows the regions to be long-lived. We examine the form of global flows for different dependences of the volume mass injection rate, $\dot{q}$, on radius and Mach number, and describe the solutions in detail. We find that three observed UCHIIR morphologies are reproduced with these models. Mach number independent flows that include a radial variation can give centre- brightened core--halo morphologies. Mach number dependent flows reproduce naturally the uniform UCHIIR morphology. In a hybrid model, including subsonic and supersonic flows, we allow a supersonic wind to shock in the ionized region. The ionized subsonic gas has a high density and so dominates the emission. The shell produced has a velocity structure very different from that of fully supersonic models. Several morphologies of spherical UCHIIR can be understood in terms of these various models; however, kinematic data are crucial as a discriminant between them.

Cometary ultracompact HII regions have been modelled as the interaction of the hypersonic wind from a moving star with the molecular cloud which surrounds the star. We here show that a similar morphology can ensue even if the star is stationary with respect to the cloud material. We assume that the HII region is within a stellar wind bubble which is strongly mass loaded: the cometary shape results from a gradient in the distribution of mass loading sources. This model circumvents problems associated with the necessarily high spatial velocities of stars in the moving star models.

We propose that a significant fraction of the ultracompact HII regions found in massive star-forming clouds are the result of the interaction of the wind and ionizing radiation from a young massive star with the clumpy molecular cloud gas in its neighbourhood. Distributed mass loading in the flow allows the compact nebulae to be long-lived. In this paper, we discuss a particularly simple case, in which the flow in the HII region is everywhere supersonic. The line profiles predicted for this model are highly characteristic, for the case of uniform mass loading. We discuss briefly other observational diagnostics of these models.