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 thermonuclear explosions that occur on the surfaces of white dwarf stars in interacting binary systems (Bode & Evans 2008). It has long been thought that the luminosity of classical novae is powered by continued nuclear burning on the surface of the white dwarf after the initial runaway (Gallaher & Starrfield 1978). However, recent observations of GeV $\gamma$-rays from classical novae have hinted that shocks internal to the nova ejecta may dominate the nova emission. Shocks have also been suggested to power the luminosity of events as diverse as stellar mergers (Metzger & Pejcha 2017), supernovae (Moriya et al. 2018), and tidal disruption events (Roth et al. 2016), but observational confirmation has been lacking. Here we report simultaneous space-based optical and $\gamma$-ray observations of the 2018 nova V906 Carinae (ASASSN-18fv), revealing a remarkable series of distinct correlated flares in both bands. The optical and $\gamma$-ray flares occur simultaneously, implying a common origin in shocks. During the flares, the nova luminosity doubles, implying that the bulk of the luminosity is shock-powered. Furthermore, we detect concurrent but weak X-ray emission from deeply embedded shocks, confirming that the shock power does not appear in the X-ray band and supporting its emergence at longer wavelengths. Our data, spanning the spectrum from radio to $\gamma$-ray, provide direct evidence that shocks can power substantial luminosity in classical novae and other optical transients.

Ultra-cool white dwarfs are among the oldest stellar remnants in the Universe. Their efficient gravitational settling and low effective temperatures are responsible for the smooth spectra they exhibit. For that reason, it is not possible to derive their radial velocities or to find the chemistry of the progenitors. The best that can be done is to infer such properties from associated sources, which are coeval. The simplest form of such a system is a common proper motion pair where one star is an evolved degenerate and the other a main-sequence star. In this work, we present the discovery of the first of such a system, the M dwarf LHS 6328 and the ultra-cool white dwarf PSO J1801+625, from the Pan-STARRS 1 3$\pi$ survey and the $Gaia$ Data Release 2. Follow-up spectra were collected covering a usable wavelength range of 3750-24500 {\AA}. Their spectra show that the white dwarf has an effective temperature of 3550K and surface gravity of $\log g = 7.45 \pm 0.13$ or $\log g = 7.49 \pm 0.13$ for a CO or He core, respectively, when compared against synthetic spectra of ultra-cool white dwarf atmosphere models. The system has slightly subsolar metallicity with -0.25 < [Fe/H] < 0.0, and a spatial velocity of (U, V, W) = ($-114.26 \pm 0.24$, $222.94 \pm 0.60$, $10.25 \pm 0.34$) km s $^{-1}$ the first radial velocity and metallicity measurements of an ultra-cool white dwarf. This makes it the first and only benchmark of its kind to date.

We present a detailed study of the 2017 eruption of the classical nova ASASSN-17pf (LMCN 2017-11a), which is located in the Large Magellanic Cloud, including data from AAVSO, ASAS-SN, SALT, SMARTS, SOAR, and the Neil Gehrels \textit{Swift} Observatory. The optical light-curve is characterized by multiple maxima (flares) on top of a slowly evolving light-curve (with a decline time, $t_2>$ 100 d). The maxima correlate with the appearance of new absorption line systems in the optical spectra characterized by increasing radial velocities. We suggest that this is evidence of multiple episodes of mass-ejection with increasing expansion velocities. The line profiles in the optical spectra indicate very low expansion velocities (FWHM $\sim$ 190 km s$^{-1}$), making this nova one of the slowest expanding ever observed, consistent with the slowly evolving light-curve. The evolution of the colors and spectral energy distribution show evidence of decreasing temperatures and increasing effective radii for the pseudo-photosphere during each maximum. The optical and infrared light-curves are consistent with dust formation 125 days post-discovery. We speculate that novae showing several optical maxima have multiple mass-ejection episodes leading to shocks that may drive $\gamma$-ray emission and dust formation.

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.

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.

V1280 Sco is one of the slowest dust-forming nova ever historically observed. We performed multi-epoch high-spatial resolution observations of the circumstellar dusty environment of V1280 Sco to investigate the level of asymmetry of the ejecta We observed V1280 Sco in 2009, 2010 and 2011 using unprecedented high angular resolution techniques. We used the NACO/VLT adaptive optics system in the J, H and K bands, together with contemporaneous VISIR/VLT mid-IR imaging that resolved the dust envelope of V1280 Sco, and SINFONI/VLT observations secured in 2011. We report the discovery of a dusty hourglass-shaped bipolar nebula. The apparent size of the nebula increased from 0.30" x 0.17" in July 2009 to 0.64" x 0.42" in July 2011. The aspect ratio suggests that the source is seen at high inclination. The central source shines efficiently in the K band and represents more than 56+/-5% of the total flux in 2009, and 87+/-6% in 2011. A mean expansion rate of 0.39+/-0.03 mas per day is inferred from the VISIR observations in the direction of the major axis, which represents a projected upper limit. Assuming that the dust shell expands in that direction as fast as the low-excitation slow ejecta detected in spectroscopy, this yields a lower limit distance to V1280 Sco of 1kpc; however, the systematic errors remain large due to the complex shape and velocity field of the dusty ejecta. The dust seems to reside essentially in the polar caps and no infrared flux is detected in the equatorial regions in the latest dataset. This may imply that the mass-loss was dominantly polar.

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.