Clues to the formation and evolution of Nuclear Star Clusters (NSCs) lie in their stellar populations. However, these structures are often very faint compared to their host galaxy, and spectroscopic analysis of NSCs is hampered by contamination of light from the rest of the system. With the introduction of wide-field IFU spectrographs, new techniques have been developed to model the light from different components within galaxies, making it possible to cleanly extract the spectra of the NSCs and study their properties with minimal contamination from the light of the rest of the galaxy. This work presents the analysis of the NSCs in a sample of 12 dwarf galaxies in the Fornax Cluster observed with MUSE. Analysis of the stellar populations and star-formation histories reveal that all the NSCs show evidence of multiple episodes of star formation, indicating that they have built up their mass further since their initial formation. The NSCs were found to have systematically lower metallicities than their host galaxies, which is consistent with a scenario for mass-assembly through mergers with infalling globular clusters, while the presence of younger stellar populations and gas emission in the core of two galaxies is indicative of in-situ star formation. We conclude that the NSCs in these dwarf galaxies likely originated as globular clusters that migrated to the core of the galaxy which have built up their mass mainly through mergers with other infalling clusters, with gas-inflow leading to in-situ star formation playing a secondary role.

Virtual Reality (VR) technology has been subject to a rapid democratization in recent years, driven in large by the entertainment industry, and epitomized by the emergence of consumer-grade, plug-and-play, room-scale VR devices. To explore the scientific potential of this technology for the field of observational astrophysics, we have created an experimental VR application: E0102-VR. The specific scientific goal of this application is to facilitate the characterization of the 3D structure of the oxygen-rich ejecta in the young supernova remnant 1E 0102.2-7219 in the Small Magellanic Cloud. Using E0102-VR, we measure the physical size of two large cavities in the system, including a (7.0$\pm$0.5) pc-long funnel structure on the far-side of the remnant. The E0102-VR application, albeit experimental, demonstrates the benefits of using human depth perception for a rapid and accurate characterization of complex 3D structures. Given the implementation costs (time-wise) of a dedicated VR application like E0102-VR, we conclude that the future of VR for scientific purposes in astrophysics most likely resides in the development of a robust, generic application dedicated to the exploration and visualization of 3D observational datasets, akin to a ``ds9-VR''.

Hickson Compact Group (HCG) 16 is a prototypical compact group of galaxies, apparently in an intermediate stage of evolution, where its galaxies are losing gas to the intra-group medium (IGrM). The wide variety of ongoing process in HCG 16 make it an ideal case study for exploring which processes are likely to dominate the late stages of evolution in compact groups. In order to build a coherent picture of the evolution of this group we make use of the multi-wavelength data available, but focus particularly on HI as a tracer of interactions and evolutionary phase. We reprocess archival VLA L-band data of HCG 16 using the multi-scale CLEAN algorithm to accurately recover diffuse features. Despite the clear disruption of the HI component of HCG 16 we find that it is not globally HI deficient, even though HCG 16a and b have lost the majority of their HI and almost 50% of the group's HI is in the IGrM, including a 160 kpc long tail extending towards the SE. This study indicates that in the recent past (~1 Gyr) galaxies HCG 16a and b likely underwent major interactions that unbound gas without triggering significant star formation. This gas was then swept away by a high speed, close encounter with NGC 848. The starburst events HCG 16c and d have triggered galactic winds which, in the case of HCG 16d, appears to have disrupted its HI reservoir. The tidal features still connected to these galaxies indicate that more HI will soon be lost to the IGrM, while that which remains in the discs will likely be consumed by star formation episodes. This is expected to result in a collection of gas-poor galaxies embedded in a diffuse HI structure, which will gradually be evaporated by the UV background, resembling the final stage of the evolutionary model of compact groups. This work is accompanied by a complete workflow to reproduce the final data products, intended to improve reproducibility. (Abridged)

Raman scattering enables unforeseen uses for the laser guide-star system of the Very Large Telescope. Here, we present the observation of one up-link sodium laser beam acquired with the ESPRESSO spectrograph at a resolution $\lambda/\Delta\lambda \sim 140'000$. In 900s on-source, we detect the pure rotational Raman lines of $^{16}$O$_2$, $^{14}$N$_2$, and $^{14}$N$^{15}$N (tentatively) up to rotational quantum numbers $J$ of 27, 24, and 9, respectively. We detect the $^{16}$O$_2$ fine-structure lines induced by the interaction of the electronic spin \textbf{S} and end-over-end rotational angular momentum \textbf{N} in the electronic ground state of this molecule up to $N=9$. The same spectrum also reveals the $\nu_{1\leftarrow0}$ rotational-vibrational Q-branch for $^{16}$O$_2$ and $^{14}$N$_2$. These observations demonstrate the potential of using laser guide-star systems as accurate calibration sources for characterizing new astronomical spectrographs.

We report the discovery of optical emission from the non-radiative shocked ejecta of three young Type Ia supernova remnants (SNRs): SNR 0519-69.0, SNR 0509-67.5, and N103B. Deep integral field spectroscopic observations reveal broad and spatially resolved [Fe XIV] 5303{\AA} emission. The width of the broad line reveals, for the first time, the reverse shock speeds. For two of the remnants we can constrain the underlying supernova explosions with evolutionary models. SNR 0519-69.0 is well explained by a standard near-Chandrasekhar mass explosion, whereas for SNR 0509-67.5 our analysis suggests an energetic sub-Chandrasekhar mass explosion. With [S XII], [Fe IX], and [Fe XV] also detected, we can uniquely visualize different layers of the explosion. We refer to this new analysis technique as "supernova remnant tomography".

We present a low-redshift (z=0.029) Damped Lyman-alpha (DLA) system in the spectrum of a background Quasi-Stellar Object (QSO). The DLA is associated with an interacting galaxy pair within a galaxy group. We detected weak Lyman-alpha emission centered at the absorption trough of the DLA. The emission was likely tracing the neutral HI reservoir around the galaxies in the interacting pair, which scattered the Lyman-alpha generated by star formation within those galaxies. We also found that the interacting pair is enveloped by a large HI cloud with $M(HI)=2\times 10^{10}M_{\odot}$. We discovered blueshifted 21cm HI emission, corresponding to M(HI)=$ 2\times10^{9}~M_{\odot}$, associated with J151225.15+012950.4 - one of the galaxies in the interacting pair. The blueshifted HI was tracing gas flowing into the galaxy from behind and towards us. Gas at similar blueshifted velocities was seen in the QSO sightline thus suggesting the presence of a filamentary structure of the order of 100kpc feeding the galaxy. We estimated a mass inflow rate of $2 M_{\odot}~yr^{-1}$ into the galaxy, which matches the star formation rate estimated from H-alpha emission. It is likely that the inflow of enormous amounts of gas has triggered star formation in this galaxy. The sudden acquisition of cold gas may lead to a starburst in this galaxy like those commonly seen in simulations.

(Abr.) Laser guide stars employed at astronomical observatories provide artificial wavefront reference sources to help correct (in part) the impact of atmospheric turbulence on astrophysical observations. Following the recent commissioning of the 4 Laser Guide Star Facility (4LGSF) on UT4 at the VLT, we characterize the spectral signature of the uplink beams from the 22W lasers to assess the impact of laser scattering from the 4LGSF on science observations. We use the MUSE optical integral field spectrograph to acquire spectra at a resolution of R~3000 of the uplink laser beams over the wavelength range of 4750\AA\ to 9350\AA. We report the first detection of laser-induced Raman scattering by N2, O2, CO2, H2O and (tentatively) CH4 molecules in the atmosphere above the astronomical observatory of Cerro Paranal. In particular, our observations reveal the characteristic spectral signature of laser photons -- but 480\AA\ to 2210\AA\ redder than the original laser wavelength of 5889.959\AA\ -- landing on the 8.2m primary mirror of UT4 after being Raman-scattered on their way up to the sodium layer. Laser-induced Raman scattering is not unique to the observatory of Cerro Paranal, but common to any astronomical telescope employing a laser-guide-star (LGS) system. It is thus essential for any optical spectrograph coupled to a LGS system to handle thoroughly the possibility of a Raman spectral contamination via a proper baffling of the instrument and suitable calibrations procedures. These considerations are particularly applicable for the HARMONI optical spectrograph on the upcoming Extremely Large Telescope. At sites hosting multiple telescopes, laser collision prediction tools also ought to account for the presence of Raman emission from the uplink laser beam(s) to avoid the unintentional contamination of observations acquired with telescopes in the vicinity of a LGS system.

Oxygen-rich young supernova remnants are valuable objects for probing the outcome of nucleosynthetic processes in massive stars, as well as the physics of supernova explosions. Observed within a few thousand years after the supernova explosion, these systems contain fast-moving oxygen-rich and hydrogen-poor filaments visible at optical wavelengths: fragments of the progenitor's interior expelled at a few 1000 km/s during the supernova explosion. Here we report the first identification of the compact object in 1E0102.2-7219 in reprocessed Chandra X-ray Observatory data, enabled via the discovery of a ring-shaped structure visible primarily in optical recombination lines of Ne I and O I. The optical ring, discovered in integral field spectroscopy observations from the Multi Unit Spectroscopic Explorer (MUSE) at the Very Large Telescope, has a radius of $(2.10\pm0.35)$ arcsec = $(0.63\pm0.11)$ pc, and is expanding at a velocity of $90.5_{-30}^{+40}$ km/s. It surrounds an X-ray point source with an intrinsic X-ray luminosity $L_{i}$ (1.2--2.0 keV)= $(1.4\pm0.2)\times10^{33}$ erg/s. The energy distribution of the source indicates that this object is an isolated neutron star: a Central Compact Object akin to those present in the Cas A and Puppis A supernova remnants, and the first of its kind to be identified outside of our Galaxy.

Type Ia supernovae are vital to our understanding of the Universe due to their use in measuring cosmological distances and their significance in enriching the interstellar medium with heavy elements. They are understood to be the thermonuclear explosions of white dwarfs, but the exact mechanism(s) leading to these explosions remains unclear. The two competing models are the single degenerate scenario, wherein a white dwarf accretes material from a companion star and explodes when it reaches the Chandrasekhar limit, and the double degenerate scenario, wherein the explosion results from a merger of two white dwarfs. Here we report results which rule out hot, luminous progenitors consistent with the single degenerate scenario for four young Type Ia supernova remnants in the Large Magellanic Cloud. Using the integral field spectrograph WiFeS, we have searched these remnants for relic nebulae ionized by the progenitor, which would persist for up to $\sim 10^5$ years after the explosion. We detected no such nebula around any of the remnants. By comparing our upper limits with photoionization simulations performed using Cloudy, we have placed stringent upper limits on the luminosities of the progenitors of these supernova remnants. Our results add to the growing evidence disfavouring the single degenerate scenario.

Using integral field data we extract the optical spectra of shocked interstellar clouds in Kepler's supernova remnant located in the inner regions of our Galaxy, as well as in the Large Magellanic Cloud (LMC), the Small Magellanic Cloud (SMC), NGC6822 and IC 1613. Using self-consistent shock modelling, we make a new determination of the chemical composition of the interstellar medium (ISM) in N, O, Ne, S, Cl and Ar in these galaxies and obtain accurate estimates of the fraction of refractory grains destroyed in the shock. By comparing our derived abundances with those obtained in recent works using observations of B stars, F supergiant stars and HII regions, we provide a new calibration for abundance scaling in the range $7.9 \lesssim 12+\log {\mathrm {O/H}} \lesssim 9.1$.

The first observations of laser guide star photons Raman-scattered by air molecules above the Very Large Telescope (VLT) were reported in June 2017. The initial detection came from the Multi-Unit Spectroscopic Explorer (MUSE) optical integral field spectrograph, following the installation of the 4 Laser Guide Star Facility (4LGSF) on the Unit Telescope 4 (UT4) of the VLT. In this Letter, we delve further into the symbiotic relationship between the 4LGSF laser guide star system, the UT4 telescope, and MUSE by monitoring the spectral contamination of MUSE observations by Raman photons over a 27 month period. This dataset reveals that dust particles deposited on the primary and tertiary mirrors of UT4 -- responsible for a reflectivity loss of ~8% at 6000{\AA} -- contribute (60$\pm5)% to the laser line fluxes detected by MUSE. The flux of Raman lines, contaminating scientific observations acquired with optical spectrographs, thus provides a new, non-invasive means to monitor the evolving scatter properties of the mirrors of astronomical telescopes equipped with laser guide star systems.

fcmaker is a python module that creates astronomical finding charts for Observing Blocks (OBs) on the p2 web server from the European Southern Observatory (ESO). It provides users with the ability to automate the creation of ESO-compliant finding charts for Service Mode and/or Visitor Mode OBs at the Very Large Telescope (VLT). The design of the fcmaker finding charts, based on an intimate knowledge of VLT observing procedures, is fine-tuned to best support night time operations. As an automated tool, fcmaker also provides observers with the means to independently check visually the observing sequence coded inside an OB. This includes, for example, the signs of telescope and position angle offsets. VLT instruments currently supported by fcmaker include MUSE (WFM-AO, WFM-NOAO, NFM), HAWK-I (AO, NOAO), and X-shooter (full support). The fcmaker code is published on a dedicated Github repository under the GNU General Public License, and is also available via pypi.

From integral field data we extract the optical spectra of 20 shocked clouds in the supernova remnant N132D in the Large Magellanic Cloud (LMC). Using self-consistent shock modelling, we derive the shock velocity, pre-shock cloud density and shock ram pressure in these clouds. We show that the [Fe X] and [Fe XIV] emission arises in faster, partially radiative shocks moving through the lower density gas near the periphery of the clouds. In these shocks dust has been effectively destroyed, while in the slower cloud shocks the dust destruction is incomplete until the recombination zone of the shock has been reached. These dense interstellar clouds provide a sampling of the general interstellar medium (ISM) of the LMC. Our shock analysis allows us to make a new determination of the ISM chemical composition in N, O, Ne, S, Cl and Ar, and to obtain accurate estimates of the fraction of refractory grains destroyed. From the derived cloud shock parameters, we estimate cloud masses and show that the clouds previously existed as typical self-gravitating Bonnor-Ebert spheres into which converging cloud shocks are now being driven.

We report the discovery of an unusual stellar system, SMSS J130522.47-293113.0. The optical spectrum is dominated by a blue continuum together with emission lines of hydrogen, neutral and ionized helium, and the N III, C III blend at around 4640-4650 Angstrom. The emission line profiles vary in strength and position on timescales as short as 1 day, while optical photometry reveals fluctuations of as much as ~0.2 mag in g on timescales as short as 10-15 min. The system is a weak X-ray source (f_{0.3-10} = (1.2 +/- 0.1) x 10^{-13} ergs cm^{-2} s^{-1} in the 0.3-10 keV band) but is not detected at radio wavelengths (3-sigma upper limit of 50 microJy at 5.5 GHz). The most intriguing property of the system, however, is the existence of two "blobs", a few arcsec in size, that are symmetrically located 3.8 arcmin (2.2 pc for our preferred system distance of ~2 kpc) each side of the central object. The blobs are detected in optical and near-IR broadband images but do not show any excess emission in Halpha images. We discuss the interpretation of the system, suggesting that the central object is most likely a nova-like CV, and that the blobs are relics of a pc-scale accretion-powered collimated outflow.

We study the optical emission from heavy element ejecta in the oxygen-rich young supernova remnant (SNR) 1E 0102.2-7219 (1E 0102) in the Small Magellanic Cloud. We have used the Multi-Unit Spectroscopic Explorer (MUSE) optical integral field spectrograph at the Very Large Telescope (VLT) on Cerro Paranal and the wide field spectrograph (WiFeS) at the ANU 2.3 m telescope at Siding Spring Observatory to obtain deep observations of 1E 0102. Our observations cover the entire extent of the remnant from below 3500{\AA} to 9350{\AA}. Our observations unambiguously reveal the presence of fast-moving ejecta emitting in [S II], [S III], [Ar III], and [Cl II]. The sulfur-rich ejecta appear more asymmetrically distributed compared to oxygen or neon, a product of carbon-burning. In addition to the forbidden line emission from products of oxygen burning (S, Ar, Cl), we have also discovered H{\alpha} and H{\beta} emission from several knots of low surface brightness, fast-moving ejecta. The presence of fast-moving hydrogen points towards a progenitor that had not entirely shed its hydrogen envelope prior to the supernova. The explosion that gave rise to 1E 0102 is therefore commensurate with a Type IIb supernova.

We present results of integral field spectroscopy of Balmer-dominated shocks in the LMC supernova remnant (SNR) N103B, carried out using the Wide Field Integral Spectrograph (WiFeS) on the 2.3 m telescope at Siding Spring Observatory in Australia. Existing X-ray studies of N103B have indicated a Type Ia supernova origin. Radiative shock emission from clumpy material surrounding the SNR may result from interaction of the forward shock with relic stellar wind material, possibly implicating a thermonuclear explosion in a single-degenerate binary system. The recently discovered Balmer-dominated shocks mark the impact of the forward shock with low density, partially neutral CSM gas, and form a partial shell encircling clumps of material exhibiting radiative shocks. The WiFeS spectra of N103B cover the range 3500 - 7000 Angstroms and reveal broad Halpha profiles up to 2360 km/s in width along the northern rim, and both Halpha and Hbeta broad profiles approximately 1300 km/s in width along the southern rim. Our fits to the Halpha line profiles indicate that in addition to from the usual broad and narrow emission components, a third component of intermediate width exists in the Balmer-dominated shocks, ranging from around 125 km/s up to 225 km/s in width. Recent Balmer-dominated shock models have predicted such an intermediate-width component will be generated in a fast neutral precursor. We derive a Sedov age of approximately 685 years for N103B from the Balmer-dominated spectra, consistent with the young age of 380 - 860 years estimated from light echo studies.

We have observed the oxygen-rich SNR 1E 0102.2-7219 with the integral field spectrograph WiFeS at Siding Spring Observatory and discovered sulfur-rich ejecta for the first time. Follow-up deep DDT observations with MUSE on the VLT (8100 s on source) have led to the additional discovery of fast- moving hydrogen as well as argon-rich and chlorine-rich material. The detection of fast-moving hydrogen knots challenges the interpretation that the progenitor of 1E 0102 was a compact core of a Wolf-Rayet star that had shed its entire envelope. In addition to the detection of hydrogen and the products of oxygen-burning, this unprecedented sharp (0.2" spaxel size at ~0.7" seeing) and deep MUSE view of an oxygen-rich SNR in the Magellanic Clouds reveals further exciting discoveries, including [Fe xiv]{\lambda}5303 and [Fe xi]{\lambda}7892 emission, which we associate with the forward shock. We present this exciting data set and discuss some of its implications for the explosion mechanism and nucleosynthesis of the associated supernova.

Aims. We study the forward shock in the oxygen-rich young supernova remnant (SNR) 1E0102.2-7219 (1E0102 in short) via optical coronal emission from [Fe XIV] and [Fe XI]: emission lines which offer an alternative method to X-rays to do so. Methods. We have used the Multi-Unit Spectroscopic Explorer (MUSE) optical integral field spectrograph at the Very Large Telescope (VLT) on Cerro Paranal to obtain deep observations of SNR 1E0102 in the Small Magellanic Cloud. Our observations cover the entire extent of the remnant with a seeing limited spatial resolution of 0.7" = 0.2 pc at the distance of 1E 0102. Results. Our MUSE observations unambiguously reveal the presence of [Fe XIV] and [Fe XI] emission in 1E0102. The emission largely arises from a thin, partial ring of filaments surrounding the fast moving O-rich ejecta in the system. The brightest [Fe XIV] and [Fe XI] emission is found along the eastern and north-western sides of 1E0102, where shocks are driven into denser ISM material, while fainter emission along the northern edge reveals the location of the forward shock in lower density gas, possibly the relic stellar wind cavity. Modeling of the eastern shocks and the photoionization precursor surrounding 1E0102, we derive a pre-shock density $n_H$ = (7.4 +-1.5) cm$^{-3}$, and a shock velocity 330 km/s < $v_s$ < 350 km/s.

Context. The distribution of elements in galaxies forms an important diagnostic tool to characterize the system's formation and evolution. This tool is however complex to use in practice, as galaxies are subject to a range of simultaneous physical processes active from pc to kpc scales. This renders observations of the full optical extent of galaxies down to sub-kpc scales essential. Aims. Using the WiFeS integral field spectrograph, we previously detected abrupt and localized variations in the gas-phase oxygen abundance of the spiral galaxy HCG91c. Here, we follow-up on these observations to map HCG91c's disk out to ~2Re at a resolution of 600pc, and characterize the non-radial variations of the gas-phase oxygen abundance in the system. Methods. We obtained deep MUSE observations of the target under ~0.6 arcsec seeing conditions. We perform both a spaxel-based and aperture-based analysis of the data to map the spatial variations of 12+log(O/H) across the disk of the galaxy. Results. We confirm the presence of rapid variations of the oxygen abundance across the entire extent of the galaxy previously detected with WiFeS, for all azimuths and radii. The variations can be separated in two categories: a) localized and associated with individual HII regions, and b) extended over kpc scales, and occurring at the boundaries of the spiral structures in the galaxy. Conclusions. Our MUSE observations suggest that the enrichment of the interstellar medium in HGC91c has proceeded preferentially along spiral structures, and less efficiently across them. Our dataset highlights the importance of distinguishing individual star-forming regions down to scales of a few 100pc when using integral field spectrographs to spatially resolve the distribution of oxygen abundances in a given system, and accurately characterize azimuthal variations and intrinsic scatter.

The concept of the X3D pathway was introduced by Vogt et al. (2016) as a new approach to sharing and publishing 3-D structures interactively in online scientific journals. The core characteristics of the X3D pathway are that: 1) it does not rely on specific software, but rather a file format (X3D), 2) it can be implemented using fully open-source tools, and 3) article readers can access the interactive models using most main stream web browsers without the need for any additional plugins. In this article, we further demonstrate the potential of the X3D pathway to visualize datasets from optical integral field spectrographs. We use recent observations of the oxygen-rich young supernova remnant 1E0102.2-7219 in the Small Magellanic Cloud to implement additional X3DOM tools & techniques and expand the range of interactions that can be offered to article readers. In particular, we present a set of javascript functions allowing the creation and interactive handling of clip planes, effectively allowing users to take measurements of distances and angles directly from the interactive model itself.

Aims. The structural inhomogeneities and kinematics of massive star nebulae are tracers of their mass-loss history. We conduct a three-dimensional morpho-kinematic analysis of the ejecta of eta Car outside its famous Homunculus nebula. Methods. We carried out the first large-scale integral field unit observations of eta Car in the optical, covering a field of view of 1'x1' centered on the star. Observations with the Multi Unit Spectroscopic Explorer (MUSE) at the Very Large Telescope (VLT) reveal the detailed three-dimensional structure of eta Car's outer ejecta. Morpho-kinematic modeling of these ejecta is conducted with the code SHAPE. Results. The largest coherent structure in eta Car's outer ejecta can be described as a bent cylinder with roughly the same symmetry axis as the Homunculus nebula. This large outer shell is interacting with the surrounding medium, creating soft X-ray emission. We establish the shape and extent of the ghost shell in front of the southern Homunculus lobe and confirm that the NN condensation can best be modeled as a bowshock in the orbital/equatorial plane. Conclusions. The SHAPE modeling of the MUSE observations indicates that the kinematics of the outer ejecta measured with MUSE can be described by a spatially coherent structure, and this structure also correlates with the extended soft X-ray emission associated with the outer debris field. The ghost shell just outside the southern Homunculus lobe hints at a sequence of eruptions within the time frame of the Great Eruption from 1837-1858 or possibly a later shock/reverse shock velocity separation. Our 3D morpho-kinematic modeling and the MUSE observations constitute an invaluable dataset to be confronted with future radiation-hydrodynamics simulations. Such a comparison may shed light on the yet elusive physical mechanism responsible for eta Car-like eruptions.

We present results of a complete integral field survey of the bright SNR N49 in the LMC, obtained with the WiFeS instrument mounted on the ANU 2.3m telescope at Siding Spring Observatory. From theoretical shock modelling with the new MAPPINGS 5.1 code we have, for the first time, subjected the optical Fe emission line spectrum of a supernova remnant to a detailed abundance and dynamical analysis covering 8 separate stages of ionisation. This allows us to derive the dust depletion factors as a function of ionisation stage. We have shown that there is substantial (30% - 90%) destruction of Fe-bearing dust grains in these fast shocks ($v_s \sim 250$ km/s), and we have confirmed that the dominant dust destruction is through the non-thermal sputtering and grain-grain collision mechanisms developed in a number of theoretical works.

Most modern astrophysical datasets are multi-dimensional; a characteristic that can nowadays generally be conserved and exploited scientifically during the data reduction/simulation and analysis cascades. Yet, the same multi-dimensional datasets are systematically cropped, sliced and/or projected to printable two-dimensional (2-D) diagrams at the publication stage. In this article, we introduce the concept of the "X3D pathway" as a mean of simplifying and easing the access to data visualization and publication via three-dimensional (3-D) diagrams. The X3D pathway exploits the facts that 1) the X3D 3-D file format lies at the center of a product tree that includes interactive HTML documents, 3-D printing, and high-end animations, and 2) all high-impact-factor & peer-reviewed journals in Astrophysics are now published (some exclusively) online. We argue that the X3D standard is an ideal vector for sharing multi-dimensional datasets, as it provides direct access to a range of different data visualization techniques, is fully-open source, and is a well defined ISO standard. Unlike other earlier propositions to publish multi-dimensional datasets via 3-D diagrams, the X3D pathway is not tied to specific software (prone to rapid and unexpected evolution), but instead compatible with a range of open-source software already in use by our community. The interactive HTML branch of the X3D pathway is also actively supported by leading peer-reviewed journals in the field of Astrophysics. Finally, this article provides interested readers with a detailed set of practical astrophysical examples designed to act as a stepping stone towards the implementation of the X3D pathway for any other dataset.

Here we aim to study the physical and kinematical characteristics of the unstudied old planetary nebula (PN) PN G342.0-01.7, which shows evidence of interaction with its surrounding interstellar medium. We used Integral Field Spectra from the Wide Field Spectrograph on the ANU 2.3 m telescope to provide spectroscopy across the whole object covering the spectral range 3400-7000 {\AA}. We formed narrow-band images to investigate the excitation structure. The spectral analysis shows that the object is a distant Peimbert Type I PN of low excitation, formally of excitation class of 0.5. The low electron density, high dynamical age, and low surface brightness of the object confirm that it is observed fairly late in its evolution. It shows clear evidence for dredge-up of CN-processed material characteristic of its class. In addition, the low peculiar velocity of 7 km s$^{-1}$ shows it to be a member of the young disk component of our Galaxy. We built a self-consistent photoionisation model for the PNe matching the observed spectrum, the H$\beta$ luminosity, and the diameter. On the basis of this we derive an effective temperature $\log T_{\rm eff} \sim 5.05$ and luminosity $1.85 < \log L < 2.25$. The temperature is much higher than might have been expected using the excitation class, proving that this can be misleading in classifying evolved PNe. PN G342.0-01.7 is in interaction with its surrounding interstellar medium through which the object is moving in the south-west direction. This interaction drives a slow shock into the outer PN ejecta. A shock model suggests that it only accounts for about 10\% of the total luminosity, but has an important effect on the global spectrum of the PN.

Galaxies in Hickson Compact Group 91 (HCG 91) were observed with the WiFeS integral field spectrograph as part of our ongoing campaign targeting the ionized gas physics and kinematics inside star forming members of compact groups. Here, we report the discovery of HII regions with abundance and kinematic offsets in the otherwise unremarkable star forming spiral HCG 91c. The optical emission line analysis of this galaxy reveals that at least three HII regions harbor an oxygen abundance ~0.15 dex lower than expected from their immediate surroundings and from the abundance gradient present in the inner regions of HCG 91c. The same star forming regions are also associated with a small kinematic offset in the form of a lag of 5-10 km/s with respect to the local circular rotation of the gas. HI observations of HCG 91 from the Very Large Array and broadband optical images from Pan-STARRS suggest that HCG 91c is caught early in its interaction with the other members of HCG 91. We discuss different scenarios to explain the origin of the peculiar star forming regions detected with WiFeS, and show that evidence point towards infalling and collapsing extra-planar gas clouds at the disk-halo interface, possibly as a consequence of long-range gravitational perturbations of HCG 91c from the other group members. As such, HCG 91c provides evidence that some of the perturbations possibly associated with the early phase of galaxy evolution in compact groups impact the star forming disk locally, and on sub-kpc scales.

The details of bulge formation via collapse, mergers, secular processes or their interplay remain unresolved. To start answering this question and quantify the importance of distinct mechanisms, we mapped a sample of three galactic bulges using data from the integral field spectrograph WiFeS on the ANU 2.3m telescope in Siding Spring Observatory. Its high resolution gratings (R=7000) allow us to present a detailed kinematic and stellar population analysis of their inner structures with classical and novel techniques. The comparison of those techniques calls for the necessity of inversion algorithms in order to understand complex substructures and separate populations. We use line-strength indices to derive SSP-equivalent ages and metallicities. Additionally, we use full spectral fitting methods, here the code STECKMAP, to extract their star formation histories. The high quality of our data allows us to study the 2D distribution of different stellar populations (i.e. young, intermediate, and old). We can identify their dominant populations based on these age-discriminated 2D light and mass contribution. In all galactic bulges studied, at least 50% of the stellar mass already existed 10 Gyrs ago, more than currently predicted by simulations. A younger component (age between 1 to 8 Gyrs) is also prominent and its present day distribution seems to be affected much more strongly by morphological structures, especially bars, than the older one. This in-depth analysis of the three bulges supports the notion of increasing complexity in their evolution, likely to be found in numerous bulge structures if studied at this level of detail, which cannot be achieved by mergers alone and require a non-negligible contribution of secular evolution.

We present a new method for inferring the metallicity (Z) and ionization parameter (q) of HII regions and star-forming galaxies using strong nebular emission lines (SEL). We use Bayesian inference to derive the joint and marginalized posterior probability density functions for Z and q given a set of observed line fluxes and an input photo-ionization model. Our approach allows the use of arbitrary sets of SELs and the inclusion of flux upper limits. The method provides a self-consistent way of determining the physical conditions of ionized nebulae that is not tied to the arbitrary choice of a particular SEL diagnostic and uses all the available information. Unlike theoretically calibrated SEL diagnostics the method is flexible and not tied to a particular photo-ionization model. We describe our algorithm, validate it against other methods, and present a tool that implements it called IZI. Using a sample of nearby extra-galactic HII regions we assess the performance of commonly used SEL abundance diagnostics. We also use a sample of 22 local HII regions having both direct and recombination line (RL) oxygen abundance measurements in the literature to study discrepancies in the abundance scale between different methods. We find that oxygen abundances derived through Bayesian inference using currently available photo-ionization models in the literature can be in good (~30%) agreement with RL abundances, although some models perform significantly better than others. We also confirm that abundances measured using the direct method are typically 0.2 dex lower than both RL and photo-ionization model based abundances.

Two-dimensional (2D) line ratio diagnostic diagrams have become a key tool in understanding the excitation mechanisms of galaxies. The curves used to separate the different regions - HII-like or else excited by an active galactic nucleus (AGN) - have been refined over time but the core technique has not evolved significantly. However, the classification of galaxies based on their emission line ratios really is a multi-dimensional problem. Here we exploit recent software developments to explore the potential of three-dimensional (3D) line ratio diagnostic diagrams. We introduce a specific set of 3D diagrams, the ZQE diagrams, which separate the oxygen abundance and the ionisation parameter of HII region-like spectra, and which also enable us to probe the excitation mechanism of the gas. By examining these new 3D spaces interactively, we define a new set of 2D diagnostics, the ZE diagnostics, which can provide the metallicity of objects excited by hot young stars, and which cleanly separate HII region-like objects from the different classes of AGNs. We show that these ZE diagnostics are consistent with the key log[NII]/H$\alpha$ vs. log[OIII]/H$\beta$ diagnostic currently used by the community. They also have the advantage of attaching a probability that a given object belongs to one class or to the other. Finally, we discuss briefly why ZQE diagrams can provide a new way to differentiate and study the different classes of AGNs in anticipation of a dedicated follow-up study.

This paper presents a detailed spectral pixel (spaxel) analysis of the ten Luminous Infrared Galaxies (LIRGs) previously observed with the Wide Field Spectrograph (WiFeS), an integral field spectrograph mounted on the ANU 2.3m telescope, and for which an abundance gradient analysis has already been presented by Rich et al. (2012). Here we use the strong emission line analysis techniques developed by Dopita et al. (2013) to measure the ionisation parameter and the oxygen abundance in each spaxel. In addition, we use the observed H$\alpha$ flux to determine the surface rate of star formation ($M_{\odot}$ yr$^{-1}$ kpc$^{-2}$) and use the [\ion{S}{2}] $\lambda\lambda6717/6731$ ratio to estimate the local pressure in the ionised plasma. We discuss the correlations discovered between these physical quantities, and use them to infer aspects of the physics of star formation in these extreme star forming environments. In particular, we find a correlation between the star formation rate and the inferred ionisation parameter. We examine the possible reasons for this correlation, and determine that the most likely explanation is that the more active star forming regions have a different distribution of molecular gas which favour higher ionisation parameters in the ionised plasma.

We present PyWiFeS, a new Python-based data reduction pipeline for the Wide Field Spectrograph (WiFeS). PyWiFeS consists of a series of core data processing routines built on standard scientific Python packages commonly used in astronomical applications. Included in PyWiFeS is an implementation of a new global optical model of the spectrograph which provides wavelengths solutions accurate to $\sim$0.05 \AA\ (RMS) across the entire detector. The core PyWiFeS package is designed to be scriptable to enable batch processing of large quantities of data, and we present a default format for handling of observation metadata and scripting of data reduction.

The optical spectroscopic followup of 27 sources belonging to a sample of 30 high-energy objects selected by positionally cross correlating the first Fermi/LAT Catalog and the ROSAT All-Sky Survey Bright Source Catalog is presented here. It has been found or confirmed that 25 of them are BL Lacertae objects (BL Lacs), while the remaining two are Galactic cataclysmic variables (CVs). This strongly suggests that the sources in the first group are responsible for the GeV emission detected with Fermi, while the two CVs most likely represent spurious associations. We thus find an 80% a posteriori probability that the sources selected by matching GeV and X-ray catalogs belong to the BL Lac class. We also show suggestions that the BL Lacs selected with this approach are probably high-synchrotron-peaked sources and in turn good candidates for the detection of ultra-high-energy (TeV) photons from them.

Recently, \citet{Nicholls12}, inspired by \emph{in situ} observations of solar system astrophysical plasmas, suggested that the electrons in \HII regions are characterised by a $\kappa$-distribution of electron energies rather than by a simple Maxwell-Boltzmann distribution. Here we have collected together the new atomic data within a modified photoionisation code to explore the effects of both the new atomic data and the $\kappa$-distribution on the strong-line techniques used to determine chemical abundances in \HII regions. By comparing the recombination temperatures ($T_{\rm rec}$) with the forbidden line temperatures ($T_{\rm FL}$) we conclude that $ \kappa \sim 20$. While representing only a mild deviation from equilibrium, this is sufficient to strongly influence abundances determined using methods which depend on measurements of the electron temperature from forbidden lines. We present a number of new emission line ratio diagnostics which cleanly separate the two parameters determining the optical spectrum of \HII regions - the ionisation parameter $q$ or $\cal{U}$ and the chemical abundance; 12+log(O/H). An automated code to extract these parameters is presented. Using the homogeneous dataset from \citet{vanZee98}, we find self-consistent results between all these different diagnostics. The systematic errors between different line ratio diagnostics are much smaller than was found in the earlier strong line work. Overall the effect of the $\kappa$-distribution on the strong line abundances derived solely on the basis of theoretical models is rather small.

Augmented Reality consists of merging live images with virtual layers of information. The rapid growth in the popularity of smartphones and tablets over recent years has provided a large base of potential users of Augmented Reality technology, and virtual layers of information can now be attached to a wide variety of physical objects. In this article, we explore the potential of Augmented Reality for astrophysical research with two distinct experiments: (1) Augmented Posters and (2) Augmented Articles. We demonstrate that the emerging technology of Augmented Reality can already be used and implemented without expert knowledge using currently available apps. Our experiments highlight the potential of Augmented Reality to improve the communication of scientific results in the field of astrophysics. We also present feedback gathered from the Australian astrophysics community that reveals evidence of some interest in this technology by astronomers who experimented with Augmented Posters. In addition, we discuss possible future trends for Augmented Reality applications in astrophysics, and explore the current limitations associated with the technology. This Augmented Article, the first of its kind, is designed to allow the reader to directly experiment with this technology.

We present 65 optical spectra of the Type Ia supernova SN 2012fr, of which 33 were obtained before maximum light. At early times SN 2012fr shows clear evidence of a high-velocity feature (HVF) in the Si II 6355 line which can be cleanly decoupled from the lower velocity "photospheric" component. This Si II 6355 HVF fades by phase -5; subsequently, the photospheric component exhibits a very narrow velocity width and remains at a nearly constant velocity of v~12,000 km/s until at least 5 weeks after maximum brightness. The Ca II infrared (IR) triplet exhibits similar evidence for both a photospheric component at v~12,000 km/s with narrow line width and long velocity plateau, as well as a high-velocity component beginning at v~31,000 km/s two weeks before maximum. SN 2012fr resides on the border between the "shallow silicon" and "core-normal" subclasses in the Branch et al. (2009) classification scheme, and on the border between normal and "high-velocity" SNe Ia in the Wang et al. (2009a) system. Though it is a clear member of the "low velocity gradient" (LVG; Benetii et al., 2005) group of SNe Ia and exhibits a very slow light-curve decline, it shows key dissimilarities with the overluminous SN 1991T or SN 1999aa subclasses of SNe Ia. SN 2012fr represents a well-observed SN Ia at the luminous end of the normal SN Ia distribution, and a key transitional event between nominal spectroscopic subclasses of SNe Ia.

Using the WiFeS integral field spectrograph, we have undertaken a series of observations of star-forming galaxies in Compact Groups. In this first paper dedicated to the project, we present the analysis of the spiral galaxy NGC838, a member of the Hickson Compact Group 16, and of its galactic wind. Our observations reveal that the wind forms an asymmetric, bipolar, rotating structure, powered by a nuclear starburst. Emission line ratio diagnostics indicate that photoionization is the dominant excitation mechanism at the base of the wind. Mixing from slow shocks (up to 20%) increases further out along the outflow axis. The asymmetry of the wind is most likely caused by one of the two lobes of the wind bubble bursting out of its HI envelope, as indicated by line ratios and radial velocity maps. The characteristics of this galactic wind suggest that it is caught early (a few Myr) in the wind evolution sequence. The wind is also quite different to the galactic wind in the partner galaxy NGC839 which contains a symmetric, shock-excited wind. Assuming that both galaxies have similar interaction histories, the two different winds must be a consequence of the intrinsic properties of NGC838 and NGC839 and their starbursts.

We demonstrate a newly developed mid-infrared planetary nebula (PN) selection technique. It is designed to enable efficient searches for obscured, previously unknown, PN candidates present in the photometric source catalogues of Galactic plane MIR sky surveys. Such selection is now possible via new, sensitive, high-to-medium resolution, MIR satellite surveys such as those from the Spitzer Space Telescope and the all-sky Wide-Field Infrared Survey Explorer (WISE) satellite missions. MIR selection is based on how different colour-colour planes isolate zones (sometimes overlapping) that are predominately occupied by different astrophysical object types. These techniques depend on the reliability of the available MIR source photometry. In this pilot study we concentrate on MIR point source detections and show that it is dangerous to take the MIR GLIMPSE (Galactic Legacy Infrared Mid-Plane Survey Extraordinaire) photometry from Spitzer for each candidate at face value without examining the actual MIR image data. About half of our selected sources are spurious detections due to the applied source detection algorithms being affected by complex MIR backgrounds and the de-blending of diffraction spikes around bright MIR point sources into point sources themselves. Nevertheless, once this additional visual diagnostic checking is performed, valuable MIR selected PN candidates are uncovered. Four turned out to have faint, compact, optical counterparts in our H-alpha survey data missed in previous optical searches. We confirm all of these as true PNe via our follow-up optical spectroscopy. This lends weight to the veracity of our MIR technique. It demonstrates sufficient robustness that high-confidence samples of new Galactic PN candidates can be extracted from these MIR surveys without confirmatory optical spectroscopy and imaging. This is problematic or impossible when the extinction is large.

We present the analysis of an IRS 5-38 {\mu}m spectrum and MIPS photometric measurements of an infrared echo near the Cassiopeia A supernova remnant observed with the Spitzer Space Telescope. We have modeled the recorded echo accounting for PAHs, quantum-heated carbon and silicate grains, as well as thermal carbon and silicate particles. Using the fact that optical light echo spectroscopy has established that Cas A originated from a type IIb supernova explosion showing an optical spectrum remarkably similar to the prototypical type IIb SN 1993J, we use the latter to construct template data input for our simulations. We are then able to reproduce the recorded infrared echo spectrum by combining the emission of dust heated by the UV burst produced at the shock breakout after the core-collapse and dust heated by optical light emitted near the visual maximum of the supernova light curve, where the UV burst and optical light curve characteristics are based on SN 1993J. We find a mean density of \sim680 H cm^{-3} for the echo region, with a size of a few light years across. We also find evidence of dust processing in the form of a lack of small PAHs with less than \sim300 carbon atoms, consistent with a scenario of PAHs destruction by the UV burst via photodissociation at the estimated distance of the echo region from Cas A. Furthermore, our simulations suggest that the weak 11 {\mu}m features of our recorded infrared echo spectrum are consistent with a strong dehydrogenated state of the PAHs. This exploratory study highlights the potential of investigating dust processing in the interstellar medium through infrared echoes.

Detection and characterization of exoplanets by direct imaging requires a coronagraph designed to deliver high contrast at small angular separation. To achieve this, an accurate control of low order aberrations, such as pointing and focus errors, is essential to optimize coronagraphic rejection and avoid the possible confusion between exoplanet light and coronagraphic leaks in the science image. Simulations and laboratory prototyping have shown that a Coronagraphic Low Order Wave-Front Sensor (CLOWFS), using a single defocused image of a reflective focal plane ring, can be used to control tip-tilt to an accuracy of 10^{-3} lambda/D. This paper demonstrates that the data acquired by CLOWFS can also be used in post-processing to calibrate residual coronagraphic leaks from the science image. Using both the CLOWFS camera and the science camera in the system, we quantify the accuracy of the method and its ability to successfully remove light due to low order errors from the science image. We also report the implementation and performance of the CLOWFS on the Subaru Coronagraphic Extreme AO (SCExAO) system and its expected on-sky performance. In the laboratory, with a level of disturbance similar to what is encountered in a post Adaptive Optics beam, CLOWFS post-processing has achieved speckle calibration to 1/300 of the raw speckle level. This is about 40 times better than could be done with an idealized PSF subtraction that does not rely on CLOWFS.

Stereoscopic visualization is seldom used in Astrophysical publications and presentations compared to other scientific fields, e.g., Biochemistry, where it has been recognized as a valuable tool for decades. We put forth the view that stereo pairs can be a useful tool for the Astrophysics community in communicating a truer representation of astrophysical data. Here, we review the main theoretical aspects of stereoscopy, and present a tutorial to easily create stereo pairs using Python. We then describe how stereo pairs provide a way to incorporate 3D data in 2D publications of standard journals. We illustrate the use of stereo pairs with one conceptual and two Astrophysical science examples: an integral field spectroscopy study of a supernova remnant, and numerical simulations of a relativistic AGN jet. We also use these examples to make the case that stereo pairs are not merely an ostentatious way to present data, but an enhancement in the communication of scientific results in publications because they provide the reader with a realistic view of multi-dimensional data, be it of observational or theoretical nature. In recognition of the ongoing 3D expansion in the commercial sector, we advocate an increased use of stereo pairs in Astrophysics publications and presentations as a first step towards new interactive and multi-dimensional publication methods.

We have used the Wide Field Spectrograph (WiFeS) on the 2.3m telescope at Siding Spring Observatory to map the [O III] 5007{\AA} dynamics of the young oxygen-rich supernova remnant N132D in the Large Magellanic Cloud. From the resultant data cube, we have been able to reconstruct the full 3D structure of the system of [O III] filaments. The majority of the ejecta form a ring of ~12pc in diameter inclined at an angle of 25 degrees to the line of sight. We conclude that SNR N132D is approaching the end of the reverse shock phase before entering the fully thermalized Sedov phase of evolution. We speculate that the ring of oxygen-rich material comes from ejecta in the equatorial plane of a bipolar explosion, and that the overall shape of the SNR is strongly influenced by the pre-supernova mass loss from the progenitor star. We find tantalizing evidence of a polar jet associated with a very fast oxygen-rich knot, and clear evidence that the central star has interacted with one or more dense clouds in the surrounding ISM.