We show that old isolated neutron stars in groups and clusters of galaxies experiencing a Quark-Nova phase (QN: an explosive transition to a quark star) may be the sources of FRBs. Each fragment ("chunk") of the ultra-relativistic QN ejecta provides a collisionless plasma for which the ambient medium (galactic/halo, the intra-group/intra-cluster medium) acts as a relativistic plasma beam. Plasma instabilities (the Buneman and the Buneman-induced thermal Weibel instabilities, successively) are induced by the beam in the chunk. These generate particle bunching and observed coherent emission at GHz frequency with a corresponding fluence in the Jy ms range. The duration (from micro-seconds to hundreds of milli-seconds), repeats (on timescales of minutes to months), frequency drift and the high occurrence rate of FRBs (a few per thousand years per galaxy) in our model are in good agreement with observed properties of FRBs. All FRBs intrinsically repeat in our model and non-repetition (i.e. the non detection of the fainter QN chunks) is detector-dependent and an artifact of the bandwidth and of the fluence sensitivity threshold. Key properties of FRB 121102 (its years of activity) and of FRB 180916.J0158+65 (its 16 day period) are recovered in our model. We give specific predictions, notably: (i) because of the viewing angle (Doppler) effect, sub-GHz detectors (CHIME) will be associated with dimmer and longer duration FRBs than GHz detectors (e.g. Parkes and ASKAP); (ii) CHIME should detect on average 5 times more FRBs from a given QN than ASKAP and Parkes; (iii) super FRBs (i.e. tens of thousands of Jy ms fluence) should be associated with intra-cluster medium QNe; (iii) monster FRBs (i.e. millions of Jy ms fluence) associated with inter-galactic medium QNe might plausibly occur with frequencies at the lower limit of the LOFAR's low-band antenna.

[Abridged] By appealing to a Quark-Nova (QN; the explosive transition of a neutron star to a quark star) in the wake of a core-collapse Supernova explosion of a massive star, we develop a unified model for long duration Gamma-ray Bursts (LGRBs) and Fast Radio Bursts (FRBs). The time delay (years to decades) between the SN and the QN and, the fragmented nature (i.e. millions of chunks) of the relativistic QN ejecta are key to yielding a robust LGRB engine. In our model, a LGRB light-curve exhibits the interaction of the fragmented QN ejecta with a turbulent (i.e. filamentary and magnetically saturated) SN ejecta which is shaped by its interaction with an underlying pulsar wind nebula (PWN). The afterglow is due to the interaction of the QN chunks, exiting the SN ejecta, with the surrounding medium. Our model can fit BAT/XRT prompt and afterglow light-curves, simultaneously with their spectra, thus yielding the observed properties of LGRBs (e.g. the Band function and the X-ray flares). We find that the Yonetoku law and the Amati law are not fundamental but phenomenological. FRBs result from coherent synchrotron emission when the QN chunks interact with non-turbulent weakly magnetized PWN-SN ejecta, where conditions are prone to the Weibel instability. Magnetic field amplification induced by the Weibel instability sets the bunching length for electrons and pairs to radiate coherently. The resulting emission frequency, luminosity, duration and dispersion measure in our model are consistent with FRB data. We find a natural unification of high-energy burst phenomena from FRBs to LGRBs including X-ray Flashes and X-ray rich GRBs as well as Super-Luminous SNe. We find a possible connection between Ultra-High Energy Cosmic Rays and FRBs and propose that a QN following a binary neutron star merger can yield a short GRB (SGRB) with fits to BAT/XRT light-curves.

We demonstrate the potential for research and outreach of mixed polygon and hydrodynamic modeling and multi-waveband rendering in the interactive 3-D astrophysical virtual laboratory Shape. In 3-D special effects and animation software for the mass media, computer graphics techniques that mix polygon and numerical hydrodynamics have become common place. In astrophysics, however, interactive modeling with polygon structures has only become available with the software Shape. Numerical hydrodynamic simulations and their visualization are usually separate, while in Shape it is integrated with the polygon modeling approach that requires no programming by the user. With two generic examples, we demonstrate that research and outreach modeling can be achieved with techniques similar to those used in the media industry with the added capability for physical rendering at any wavelength band, yielding more realistic radiation modeling. Furthermore, we show how the hydrodynamics and the polygon mesh modeling can be mixed to achieve results that are superior to those obtained using either one of these modeling techniques alone.

Any simulation of the r-process is affected by uncertainties in our present knowledge of nuclear physics quantities and astrophysical conditions. It is common to quantify the impact of these uncertainties through a global sensitivity metric, which is then used to identify specific nuclides that would be most worthwhile to measure experimentally. Using descriptive statistics, we assess a set of metrics used in previous sensitivity studies, as well as a new logarithmic measure. For certain neutron-rich nuclides lying near the r-process path for the typical hot-wind scenario, we find opposing conclusions on their relative sensitivity implied by different metrics, although they all generally agree which ones are the most sensitive nuclei. The underlying reason is that sensitivity metrics which simply sum over variations in the r-process distribution depend on the scaling used in the baseline, which often varies between simulations. We show that normalization of the abundances causes changes in the reported sensitivity factors and recommend reporting a minimized F statistic in addition to a scale estimation for rough calibration to be used when comparing tables of sensitivity factors from different studies.

We show that a Quark-Nova (QN; the explosive transition of a neutron star to a quark star) occurring a few days following the supernova explosion of an Oxygen-type Wolf-Rayet (WO) star can account for the intriguing features of ASASSN-15lh, including its extreme energetics, its double-peaked light-curve and the evolution of its photospheric radius and temperature. A two-component configuration of the homologously expanding WO remnant (an extended envelope and a compact core) is used to harness the kinetic energy (>10^52 ergs) of the QN ejecta. The delay between the WO SN and the QN yields a large (~ 10^4 Rsun) envelope which when energized by the QN ejecta/shock gives the first peak in our model. As the envelope's photosphere recedes into the slowly expanding, hot and insulated, denser core (initially heated by the QN shock) a second hump emerges. The spectrum in our model should reflect the composition of an WO SN remnant re-heated by a QN going off in its wake.

Quark-novae leave behind quark stars with a surrounding metal-rich fall-back (ring-like) material. These compact remnants have high magnetic fields and are misconstrued as magnetars; however, several observational features allow us to distinguish a quark star (left behind by a quark-nova) from a neutron star with high magnetic field. In our model, bursting activity is expected from intermittent accretion events from the surrounding fall-back debris leading to X-ray bursts (in the case of a Keplerian ring) or gamma ray bursts (in the case of a co-rotating shell). The details of the spectra are described by a constant background X-ray luminosity from the expulsion of magnetic flux tubes which will be temporarily buried by bursting events caused by accretion of material onto the quark star surface. These accretion events emit high energy photons and heat up the quark star and surrounding debris leading to hot spots which may be observable as distinct blackbodies. Additionally, we explain observed spectral line features as atomic lines from r-process material and explain an observed anti-glitch in an AXP as the transfer of angular momentum from a surrounding Keplerian disk to the quark star.

The explosive transition of a massive neutron star to a quark star (the Quark-Nova, QN) releases in excess of ~ 10^52 erg in kinetic energy which can drastically impact the surrounding environment of the QN. A QN is triggered when a neutron star gains enough mass to reach the critical value for quark deconfinement to happen in the core. In binaries, a neutron star has access to mass reservoirs (e.g. accretion from a companion or from a Common Envelope, CE). We explain observed light-curves of hydrogen-poor superluminous Supernovae (SLSNe Ia) in the context of a QN occurring in the second CE phase of a massive binary. In particular this model gives good fits to light-curves of SLSNe with double-humped light-curves. Our model suggests the QN as a mechanism for CE ejection and that they be taken into account during binary evolution. In a short period binary with a white dwarf companion, the neutron star can quickly grow in mass and experience a QN event. Part of the QN ejecta collides with the white dwarf, shocking, compressing, and heating it to driving a thermonuclear runaway producing a SN Ia impostor (a QN-Ia). Unlike "normal" Type Ia supernovae where no compact remnant is formed, a QN-Ia produces a quark star undergoing rapid spin-down providing additional power along with the 56Ni decay energy. Type Ia SNe are used as standard candles and contamination of this data by QNe-Ia can infer an incorrect cosmology.

A Quark-Nova (QN, the sudden transition from a neutron star into a quark star) which occurs in the second common envelope (CE) phase of a massive binary (Ouyed et al., 2015a&b), gives excellent fits to super-luminous, hydrogen-poor, Supernovae (SLSNe) with double-peaked light curves including DES13S2cmm, SN 2006oz and LSQ14bdq (http://www.quarknova.ca/LCGallery.html). In our model, the H envelope of the less massive companion is ejected during the first CE phase while the QN occurs deep inside the second, He-rich, CE phase after the CE has expanded in size to a radius of a few tens to a few thousands solar radii, this yields the first peak in our model. The ensuing merging of the quark star with the CO core leads to black hole formation and accretion explaining the second long-lasting peak. We study a sample of 8 SLSNe Ic with double-humped light-curves. Our model provides good fits to all of these with a universal explosive energy of 2x10^52 erg (which is the kinetic energy of the QN ejecta) for the first hump. The late-time emissions seen in iPTF13ehe and LSQ14bdq are fit with a shock interaction between the outgoing He-rich (i.e second) CE and the previously ejected H-rich (i.e first) CE.

FRBs are puzzling, millisecond, energetic radio transients with no discernible source; observations show no counterparts in other frequency bands. The birth of a quark star from a parent neutron star experiencing a quark nova - previously thought undetectable when born in isolation - provides a natural explanation for the emission characteristics of FRBs. The generation of unstable r-process elements in the quark nova ejecta provides millisecond exponential injection of electrons into the surrounding strong magnetic field at the parent neutron star's light cylinder via $\beta$-decay. This radio synchrotron emission has a total duration of hundreds of milliseconds and matches the observed spectrum while reducing the inferred dispersion measure by approximately 200 cm$^{-3}$ pc. The model allows indirect measurement of neutron star magnetic fields and periods in addition to providing astronomical measurements of $\beta$-decay chains of unstable neutron rich nuclei. Using this model, we can calculate expected FRB average energies ($\sim$ 10$^{41}$ ergs), spectra shapes and provide a theoretical framework for determining distances.

We show that by appealing to a Quark-Nova (the explosive transition of a neutron star to a quark star) occurring in an He-HMXB system we can account for the lightcurve of the first superluminous SN, DES13S2cmm, discovered by the Dark-Energy Survey. The neutron star's explosive conversion is triggered as a result of accretion during the He-HMXB second Common Envelope phase. The dense, relativistic, Quark-Nova ejecta in turn energizes the extended He-rich Common Envelope in an inside-out shock heating process. We find an excellent fit (reduced chi^2 of 1.09) to the bolometric light-curve of SN DES13S2cmm including the late time emission, which we attribute to Black Hole accretion following the conversion of the Quark Star to a Black hole.

LSQ14bdq and SN 2006oz are super-luminous, hydrogen-poor, SNe with double-humped light curves. We show that a Quark-Nova (QN; explosive transition of the neutron star to a quark star) occurring in a massive binary, experiencing two Common Envelope (CE) phases, can quantitatively explain the light curves of LSQ14bdq and SN 2006oz. The more massive component (A) explodes first as a normal SN, yielding a Neutron Star which ejects the hydrogen envelope of the companion when the system enters its first CE phase. During the second CE phase, the NS spirals into and inflates the second He-rich CE. In the process it gains mass and triggers a Quark-Nova, outside of the CO core, leaving behind a Quark Star. The first hump in our model is the QN shock re-energizing the expanded He-rich CE. The QN occurs when the He-rich envelope is near maximum size (~ 1000R_sun) and imparts enough energy to unbind and eject the envelope. Subsequent merging of the Quark Star with the CO core of component B, driven by gravitational radiation, turns the Quark star to a Black Hole. The ensuing Black Hole accretion provides sufficient power for the second brighter and long lasting hump. Our model suggests a possible connection between SLSNe-I and type Ic-BL SNe which occur when the Quark Nova is triggered inside the CO core. We estimate the rate of QNe in massive binaries during the second CE phase to be ~ 5x10^(-5) of that of core-collapse SNe.

Burn-UD is a hydrodynamic combustion code used to model the phase transition of hadronic to quark matter with particular application to the interior of neutron stars. Burn-UD models the flame micro-physics for different equations of state (EoS) on both sides of the interface, i.e. for both the ash (up-down-strange quark phase) and the fuel (up-down quark phase). It also allows the user to explore strange quark seeding produced by different processes including DM annihilation inside neutron stars. The simulations provide a physical window to diagnose whether the combustion process will simmer quietly and slowly, lead to a transition from deflagration to detonation or a (quark) core-collapse explosion. Such an energetic phase transition (a Quark-Nova) would have consequences in high-energy astrophysics and could aid in our understanding of many still enigmatic astrophysical transients. Furthermore, having a precise understanding of the phase transition dynamics for different EoSs could aid further in constraining the nature of the non-perturbative regimes of QCD in general. We hope that Burn-UD will evolve into a platform/software to be used and shared by the QCD community exploring the phases of Quark Matter and astrophysicists working on Compact Stars.

We show that the low-velocity 56Ni decay lines detected earlier than expected in the type Ia SN 2014J find an explanation in the Quark-Nova Ia model which involves the thermonuclear explosion of a tidally disrupted sub-Chandrasekhar White Dwarf in a tight Neutron-Star-White-Dwarf binary system. The explosion is triggered by impact from the Quark-Nova ejecta on the WD material; the Quark-Nova is the explosive transition of the Neutron star to a Quark star triggered by accretion from a CO torus (the circularized WD material). The presence of a compact remnant (the Quark Star) provides: (i) an additional energy source (spin-down power) which allows us to fit the observed light-curve including the steep early rise; (ii) a central gravitational potential which slows down some of the 56Ni produced to velocities of a few 1000 km/s. In our model, the 56Ni decay lines become optically visible at ~20 days from explosion time in agreement with observations. We list predictions that can provide important tests for our model.

We present the results of large scale, three-dimensional magneto-hydrodynamics simulations of disc-winds for different initial magnetic field configurations. The jets are followed from the source to 90 AU scale, which covers several pixels of HST images of nearby protostellar jets. Our simulations show that jets are heated along their length by many shocks. We compute the emission lines that are produced, and find excellent agreement with observations. The jet width is found to be between 20 and 30 AU while the maximum velocities perpendicular to the jet is found to be up to above 100 km/s. The initially less open magnetic field configuration simulations results in a wider, two-component jet; a cylindrically shaped outer jet surrounding a narrow and much faster, inner jet. These simulations preserve the underlying Keplerian rotation profile of the inner jet to large distances from the source. However, for the initially most open magnetic field configuration the kink mode creates a narrow corkscrew-like jet without a clear Keplerian rotation profile and even regions where we observe rotation opposite to the disc (counter-rotating). The RW Aur jet is narrow, indicating that the disc field in that case is very open meaning the jet can contain a counter-rotating component that we suggests explains why observations of rotation in this jet has given confusing results. Thus magnetized disc winds from underlying Keplerian discs can develop rotation profiles far down the jet that are not Keplerian.

Modern instrumentation in radioastronomy constitutes a valuable tool for studying the Universe: ALMA has reached unprecedented sensitivities and spatial resolution, while Herschel/HIFI has opened a new window for probing molecular warm gas (~50-1000 K). On the other hand, the software SHAPE has emerged in the past few years as a standard tool for determining the morphology and velocity field of different kinds of gaseous emission nebulae via spatio-kinematical modelling. SHAPE implements radiative transfer solving, but it is only available for atomic species and not for molecules. Being aware of the growing importance of the development of tools for simplifying the analyses of molecular data, we introduce shapemol, a complement to SHAPE, with which we intend to fill the so-far under-developed molecular niche. shapemol enables user-friendly, spatio-kinematic modelling with accurate non-LTE calculations of excitation and radiative transfer in CO lines. It allows radiative transfer solving in the 12CO and 13CO J=1-0 to J=17-16 lines, but its implementation permits easily extending the code to different molecular species. shapemol allows easily generating synthetic maps and line profiles to match against interferometric or single-dish observations. We fully describe shapemol and discuss its limitations and the sources of uncertainty to be expected in the final synthetic profiles or maps. As an example of the power and versatility of shapemol, we build a model of the molecular envelope of the planetary nebula NGC 6302 and compare it with 12CO and 13CO J=2-1 interferometric maps from SMA and high-J transitions from Herschel/HIFI. We find the molecular envelope to have a complex, broken ring-like structure with an inner, hotter region and several 'fingers' and high-velocity blobs, emerging outwards from the plane of the ring. We derive a mass of 0.11 Msun for the molecular envelope.

Radio observations of novae in outburst are of particular interest due to the physical parameters that may be retrieved from fitting the radio light curves. Most models that have fitted previous data assumed spherical symmetry however, it is becoming more and more clear that this is not the case. We explore morpho-kinematical techniques to retrieve the free-free radio light curves of non-spherical models and explore the effects of a non-spherical outburst on the physical parameters. In particular, we find that we may have been over estimating the ejected masses in the outburst of non-spherical novae.

We show that the explosive transition of the neutron star (NS) to a quark star (QS) (a Quark Nova) in Cassiopeia A (Cas A) a few days following the SN proper can account for several of the puzzling kinematic and nucleosynthetic features observed. The observed decoupling between Fe and 44Ti and the lack of Fe emission within 44Ti regions is expected in the QN model owing to the spallation of the inner SN ejecta by the relativistic QN neutrons. Our model predicts the 44Ti to be more prominent to the NW of the central compact object (CCO) than in the SE and little of it along the NE-SW jets, in agreement with NuStar observations. Other intriguing features of Cas A such as the lack of a pulsar wind nebula (PWN) and the reported a few percent drop of the CCO temperature over a period of 10 years are also addressed.

Observations of novae at radio frequencies provide us with a measure of the total ejected mass, density profile and kinetic energy of a nova eruption. The radio emission is typically well characterized by the free-free emission process. Most models to date have assumed spherical symmetry for the eruption, although it has been known for as long as there have been radio observations of these systems, that spherical eruptions are to simplistic a geometry. In this paper, we build bipolar models of the nova eruption, assuming the free-free process, and show the effects of varying different parameters on the radio light curves. The parameters considered include the ratio of the minor- to major-axis, the inclination angle and shell thickness (further parameters are provided in the appendix). We also show the uncertainty introduced when fitting spherical model synthetic light curves to bipolar model synthetic light curves. We find that the optically thick phase rises with the same power law ($S_{\nu} \propto t^2$) for both the spherical and bipolar models. In the bipolar case there is a "plateau" phase -- depending on the thickness of the shell as well as the ratio of the minor- to major-axis -- before the final decline, that follows the same power law ($S_{\nu} \propto t^{-3}$) as in the spherical case. Finally, fitting spherical models to the bipolar model synthetic light curves requires, in the worst case scenario, doubling the ejected mass, more than halving the electron temperature and reducing the shell thickness by nearly a factor of 10. This implies that in some systems we have been over predicting the ejected masses and under predicting the electron temperature of the ejecta.

In the Quark-Nova model, AXPs are quark stars surrounded by a degenerate iron-rich Keplerian ring (a few stellar radii away). AXP bursts are caused by accretion of chunks from the inner edge of the ring following magnetic field penetration. For bright bursts, the inner disk is prone to radiation induced warping which can tilt it into counter-rotation (i.e. retrograde). For AXP 1E2259+586, the 2002 burst satisfies the condition for the formation of a retrograde inner ring. We hypothesize the 2002 burst reversed the inner ring setting the scene for the 2012 outburst and "anti-glitch" when the retrograde inner ring was suddenly accreted leading to the basic observed properties of the 2012 event.

[Aims:] We present r-Java 2.0, a nucleosynthesis code for open use that performs r-process calculations as well as a suite of other analysis tools. [Methods:] Equipped with a straightforward graphical user interface, r-Java 2.0 is capable of; simulating nuclear statistical equilibrium (NSE), calculating r-process abundances for a wide range of input parameters and astrophysical environments, computing the mass fragmentation from neutron-induced fission as well as the study of individual nucleosynthesis processes. [Results:] In this paper we discuss enhancements made to this version of r-Java, paramount of which is the ability to solve the full reaction network. The sophisticated fission methodology incorporated into r-Java 2.0 which includes three fission channels (beta-delayed, neutron-induced and spontaneous fission) as well as computation of the mass fragmentation is compared to the upper limit on mass fission approximation. The effects of including beta-delayed neutron emission on r-process yield is studied. The role of coulomb interactions in NSE abundances is shown to be significant, supporting previous findings. A comparative analysis was undertaken during the development of r-Java 2.0 whereby we reproduced the results found in literature from three other r-process codes. This code is capable of simulating the physical environment of; the high-entropy wind around a proto-neutron star, the ejecta from a neutron star merger or the relativistic ejecta from a quark nova. As well the users of r-Java 2.0 are given the freedom to define a custom environment. This software provides an even platform for comparison of different proposed r-process sites and is available for download from the website of the Quark-Nova Project: http://quarknova.ucalgary.ca/

[Context:] This article is the second in a two part series introducing r-Java 2.0, a nucleosynthesis code for open use that performs r-process calculations and provides a suite of other analysis tools. [Aims:] The first paper discussed the nuclear physics inherent to r-Java 2.0 and in this article the astrophysics incorporated into the software will be detailed. [Methods:] R-Java 2.0 allows the user to specify the density and temperature evolution for an r-process simulation. Defining how the physical parameters (temperature and density) evolve can effectively simulate the astrophysical conditions for the r-process. Within r-Java 2.0 the user has the option to select astrophysical environments which have unique sets of input parameters available for the user to adjust. In this work we study three proposed r-process sites; neutrino-driven winds around a proto-neutron star, ejecta from a neutron star merger and ejecta from a quark nova. The underlying physics that define the temperature and density evolution for each site is described in this work. [Results:] In this paper a survey of the available parameters for each astrophysical site is undertaken and the effect on final r-process abundance is compared. The resulting abundances for each site are also compared to solar observations both independently and in concert. R-Java 2.0 is available for download from the website of the Quark-Nova Project: http://quarknova.ucalgary.ca/

The exact mechanism behind Type Ia supernovae (SNe-Ia) and the nature of the progenitors is poorly understood, although several theories vie for supremacy. Their secured importance to the study of cosmology necessitates a resolution to this problem. Observations of nearby SNe-Ia are therefore critical to determine which theory, if any, is the correct one. SN 2014J discovered at a relatively close 3.5 Mpc in the galaxy M82 provides such an opportunity. In this paper we give specific predictions for SN 2014J in the context of the Quark Nova Ia (QN-Ia) model. Predictions include X-ray luminosities just prior and hundreds of days after the explosion, light curve "glitches", neutrino emission, heavy element creation, and gravitational signatures.

The accelerated expansion of the Universe was proposed through the use of Type-Ia SNe as standard candles. The standardization depends on an empirical correlation between the stretch/color and peak luminosity of the light curves. The use of Type Ia SN as standard candles rests on the assumption that their properties (and this correlation) do not vary with red-shift. We consider the possibility that the majority of Type-Ia SNe are in fact caused by a Quark-Nova detonation in a tight neutron-star-CO-white-dwarf binary system; a Quark-Nova Ia. The spin-down energy injected by the Quark Nova remnant (the quark star) contributes to the post-peak light curve and neatly explains the observed correlation between peak luminosity and light curve shape. We demonstrate that the parameters describing Quark-Novae Ia are NOT constant in red-shift. Simulated Quark-Nova Ia light curves provide a test of the stretch/color correlation by comparing the true distance modulus with that determined using SN light curve fitters. We determine a correction between the true and fitted distance moduli which when applied to Type-Ia SNe in the Hubble diagram recovers the Omega_M = 1 cosmology. We conclude that Type-Ia SNe observations do not necessitate the need for an accelerating expansion of the Universe (if the observed SNe-Ia are dominated by QNe-Ia) and by association the need for Dark Energy.

In recent years a number of double-humped supernovae have been discovered. This is a feature predicted by the dual-shock Quark-Nova model where a SN explosion is followed (a few days to a few weeks later) by a Quark-Nova explosion. SN 2009ip and SN 2010mc are the best observed examples of double-humped SNe. Here, we show that the dual-shock Quark-Nova model naturally explains their lightcurves including the late time emission, which we attribute to the interaction between the mixed SN and QN ejecta and the surrounding CSM. Our model applies to any star (O-stars, LBVs, WRs etc.) provided that the SN explosion mass is ~ 20M_sun which point to the conditions for forming a Quark-Nova.

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

We report the results of radio interferometric observations of the 21-micron source IRAS 22272+5435 in the CO J=2-1 line. 21-micron sources are carbon-rich objects in the post-AGB phase of evolution which show an unidentified emission feature at 21 micron. Since 21-micron sources usually also have circumstellar molecular envelopes, the mapping of CO emission from the envelope will be useful in tracing the nebular structure. From observations made with the Combined Array for Research in Millimeter-wave Astronomy (CARMA), we find that a torus and spherical wind model can explain only part of the CO structure. An additional axisymmetric region created by the interaction between an invisible jet and ambient material is suggested.

Recent observational surveys have uncovered the existence of super-luminous supernovae (SLSNe). While several possible explanations have been put forth, a consensus description for SLSNe has yet to be found. In this work we study the light curves of eight SLSNe in the context of dual-shock quark novae. We find that progenitor stars in the range of 25-35 $M_{\sun}$ provide ample energy to power each light curve. An examination into the effects of varying the physical properties of a dual-shock quark nova on light curve composition is undertaken. We conclude that the wide variety of SLSN light curve morphologies can be explained predominantly by variations in the length of time between supernova and quark nova. Our analysis shows that a singular H$\alpha$ spectral profile found in three SLSNe can be naturally described in the dual-shock quark nova scenario. Predictions of spectral signatures unique to the dual-shock quark nova are presented.

During 2003 and 2004 the Anomalous X-Ray Pulsar XTE J1810-197 went through a series of four bursts. The spectrum in the tail of one of these bursts shows a strong, significant emission feature ~13 keV, thereby encoding a wealth of information about the environment surrounding this object. In this paper we analyse this emission feature considering both cyclotron and atomic emission processes and weigh our findings against three leading AXP models: the Magnetar model, Fall-back disk model and the Quark nova model. We find that atomic emission from Rubidium within a Keplerian ring ($\sim$15 km from a compact object of $\sim 2M_\odot$) is the most consistent scenario with the observations, supporting the Quark nova model. Cyclotron emission from an atmosphere a few hundred meters thick also fits the feature well, but is ruled out on account of its positional coincidence in three separate AXP sources.

NGC 2392 is a young double-shell planetary nebula (PN). Its intrinsic structure and shaping mechanism are still not fully understood. In this paper we present new spectroscopic observations of NGC 2392. The slits were placed at two different locations to obtain the spectra of the inner and outer regions. Several [Fe III] lines are clearly detected in the inner region. We infer that NGC 2392 might have an intrinsic structure similar to the bipolar nebula Mz 3, which also exhibits a number of [Fe III}] lines arising from the central regions. In this scenario, the inner and outer regions of NGC 2392 correspond to the inner lobes and the outer outflows of Mz 3, respectively. We construct a three-dimensional morpho-kinematic model to examine our hypothesis. We also compare the physical conditions and chemical composition of the inner and outer regions, and discuss the implications on the formation of this type of PN.

Interactive visualization and simulation of astrophysical phenomena help astronomers and enable digital planetariums and television documentaries to take their spectators on a journey into deep space to explore the astronomical wonders of our universe in 3D.

We observed four "water fountain" sources in the CO J=3-2 line emission with the Atacama Submillimeter Telescope Experiment (ASTE) 10 m telescope in 2010-2011. The water fountain sources are evolved stars that form high-velocity collimated jets traced by water maser emission. The CO line was detected only from IRAS 16342-3814. The present work confirmed that the ^{12}CO to ^{13}CO line intensity ratio is ~1.5 at the systemic velocity. We discuss the origins of the very low ^{12}CO to ^{13}CO intensity ratio, as possible evidence for the "hot-bottom burning" in an oxygen-rich star, and the CO intensity variation in IRAS 16342-3814.

The extremely luminous supernova 2006gy (SN 2006gy) is among the most energetic ever observed. The peak brightness was 100 times that of a typical supernova and it spent an unheard of 250 days at magnitude -19 or brighter. Efforts to describe SN 2006gy have pushed the boundaries of current supernova theory. In this work we aspire to simultaneously reproduce the photometric and spectroscopic observations of SN 2006gy using a quark nova model. This analysis considers the supernova explosion of a massive star followed days later by the quark nova detonation of a neutron star. We lay out a detailed model of the interaction between the supernova envelope and the quark nova ejecta paying special attention to a mixing region which forms at the inner edge of the supernova envelope. This model is then fit to photometric and spectroscopic observations of SN 2006gy. This QN model naturally describes several features of SN 2006gy including the late stage light curve plateau, the broad H{\alpha} line and the peculiar blue H{\alpha} absorption. We find that a progenitor mass between 20Msun and 40Msun provides ample energy to power SN 2006gy in the context of a QN.

We present r-Java, an r-process code for open use, that performs r-process nucleosynthesis calculations. Equipped with a simple graphical user interface, r-Java is capable of carrying out nuclear statistical equilibrium (NSE) as well as static and dynamic r-process calculations for a wide range of input parameters. In this introductory paper, we present the motivation and details behind r-Java, and results from our static and dynamic simulations. Static simulations are explored for a range of neutron irradiation and temperatures. Dynamic simulations are studied with a parameterized expansion formula. Our code generates the resulting abundance pattern based on a general entropy expression that can be applied to degenerate as well as non-degenerate matter, allowing us to track the rapid density and temperature evolution of the ejecta during the initial stages of ejecta expansion. At present, our calculations are limited to the waiting-point approximation. We encourage the nuclear astrophysics community to provide feedback on the code and related documentation, which is available for download from the website of the Quark-Nova Project: http://quarknova.ucalgary.ca/

High-resolution HST imaging of the compact planetary nebula NGC 6644 has revealed two pairs of bipolar lobes and a central ring lying close to the plane of the sky. From mid-infrared imaging obtained with the Gemini Telescope, we have found a dust torus which is oriented nearly perpendicular to one pair of the lobes. We suggest that NGC 6644 is a multipolar nebula and have constructed a 3-D model which allows the visualization of the object from different lines of sight. These results suggest that NGC 6644 may have similar intrinsic structures as other multipolar nebulae and the phenomenon of multipolar nebulosity may be more common than previously believed.

We present a new analysis of kinematic data of the young planetary nebula BD+30 3639. The data include spectroscopic long-slit and internal proper motion measurements. In this paper we also introduce a new type of mapping of kinematic proper motion data that we name "criss-cross" mapping. It basically consists of finding all points where extended proper motion vectors cross converge. From the crossing points a map is generated which helps to interpret the kinematic data. From the criss-cross mapping of BD+30 3639, we conclude that the kinematic center is approximately 0.5 arcsec offset to the South-East from the central star. The mapping does also show evidence for a non-homologous expansion of the nebula that is consistent with a disturbance aligned with the bipolar molecular bullets.

We report on the results of a radio interferometric observation of NGC7027 in the CO J=2-1 and 13CO J=2-1 lines. The results are analyzed with morpho-kinematic models developed from the software tool Shape. Our goal is to reveal the morpho-kinematic properties of the central region of the nebula, and to explore the nature of unseen high-velocity jets that may have created the characteristic structure of the central region consisting of molecular and ionized components. A simple ellipsoidal shell model explains the intensity distribution around the systemic velocity, but the high velocity features deviate from the ellipsoidal model. Through the Shape automatic reconstruction model, we found a possible trail of a jet only in one direction, but no other possible holes were created by the passage of a jet.

We present a flexible interactive 3D morpho-kinematical modeling application for astrophysics. Compared to other systems, our application reduces the restrictions on the physical assumptions, data type and amount that is required for a reconstruction of an object's morphology. It is one of the first publicly available tools to apply interactive graphics to astrophysical modeling. The tool allows astrophysicists to provide a-priori knowledge about the object by interactively defining 3D structural elements. By direct comparison of model prediction with observational data, model parameters can then be automatically optimized to fit the observation. The tool has already been successfully used in a number of astrophysical research projects.

We observed CO J=3-2 emission from the "water fountain" sources, which exhibit high-velocity collimated stellar jets traced by water maser emission, with the Atacama Submillimeter Telescope Experiment (ASTE) 10 m telescope. We detected the CO emission from two sources, IRAS 16342-3814 and IRAS 18286-0959. The IRAS 16342-3814 CO emission exhibits a spectrum that is well fit to a Gaussian profile, rather than to a parabolic profile, with a velocity width (FWHM) of 158+/-6 km/s and an intensity peak at VLSR = 50+/-2 km/s. The mass loss rate of the star is estimated to be ~2.9x10^-5 M_sun/yr. Our morpho-kinematic models suggest that the CO emission is optically thin and associated with a bipolar outflow rather than with a (cold and relatively small) torus. The IRAS 18286-0959 CO emission has a velocity width (FWHM) of 3.0+/-0.2 km/s, smaller than typically seen in AGB envelopes. The narrow velocity width of the CO emission suggests that it originates from either an interstellar molecular cloud or a slowly-rotating circumstellar envelope that harbors the water maser source.

Spectral line surveys are useful since they allow identification of new molecules and new lines in uniformly calibrated data sets. Nonetheless, large portions of the sub-millimetre spectral regime remain unexplored due to severe absorptions by H2O and O2 in the terrestrial atmosphere. The purpose of the measurements presented here is to cover wavelength regions at and around 0.55 mm -- regions largely unobservable from the ground. Using the Odin astronomy/aeronomy satellite, we performed the first spectral survey of the Orion KL molecular cloud core in the bands 486--492 and 541--576 GHz with rather uniform sensitivity (22--25 mK baseline noise). Odin's 1.1 m size telescope, equipped with four cryo-cooled tuneable mixers connected to broad band spectrometers, was used in a satellite position-switching mode. Two mixers simultaneously observed different 1.1 GHz bands using frequency steps of 0.5 GHz (25 hours each). An on-source integration time of 20 hours was achieved for most bands. The entire campaign consumed ~1100 orbits, each containing one hour of serviceable astro-observation. We identified 280 spectral lines from 38 known interstellar molecules (including isotopologues) having intensities in the range 80 to 0.05 K. An additional 64 weak lines remain unidentified. Apart from the ground state rotational 1(1,0)--1(0,1) transitions of ortho-H2O, H218O and H217O, the high energy 6(2,4)--7(1,7) line of para-H2O and the HDO(2,0,2--1,1,1) line have been observed, as well as the 1,0--0,1 lines from NH3 and its rare isotopologue 15NH3. We suggest assignments for some unidentified features, notably the new interstellar molecules ND and SH-. Severe blends have been detected in the line wings of the H218O, H217O and 13CO lines changing the true linewidths of the outflow emission.

We investigate the physical and chemical conditions in a typical star forming region, including an unbiased search for new molecules in a spectral region previously unobserved. Due to its proximity, the Orion KL region offers a unique laboratory of molecular astrophysics in a chemically rich, massive star forming region. Several ground-based spectral line surveys have been made, but due to the absorption by water and oxygen, the terrestrial atmosphere is completely opaque at frequencies around 487 and 557 GHz. To cover these frequencies we used the Odin satellite to perform a spectral line survey in the frequency ranges 486-492 GHz and 541-577 GHz, filling the gaps between previous spectral scans. Odin's high main beam efficiency and observations performed outside the atmosphere make our intensity scale very well determined. We observed 280 spectral lines from 38 molecules including isotopologues, and, in addition, 64 unidentified lines. The beam-averaged emission is dominated by CO, H2O, SO2, SO, 13CO and CH3OH. Species with the largest number of lines are CH3OH, (CH33)2O, SO2, 13CH3OH, CH3CN and NO. Six water lines are detected including the ground state rotational transition o-H2O, its isotopologues o-H218O and o-H217O, the Hot Core tracing p-H2O transition 6(2,4)-7(1,7), and the 2(0, 2)-1(1,1) transition of HDO. Other lines of special interest are the 1_0-0_0 transition of NH3 and its isotopologue 15NH3. Isotopologue abundance ratios of D/H, 12C/13C, 32S/34S, 34S/33S, and 18O/17O are estimated. The temperatures, column densities and abundances in the various subregions are estimated, and we find very high gas-phase abundances of H2O, NH3, SO2, SO, NO, and CH3OH. A comparison with the ice inventory of ISO sheds new light on the origin of the abundant gas-phase molecules.

In search for deviations from homologous expansion in planetary nebulae we present a 3D morphokinematical model of NGC 7009. The model has been constructed with {\em Shape} based on PV diagrams from the literature and HST images. We find that the data are consistent with a radial velocity field with increased gradient at high latitudes compared to the equatorial region (Model 1). In a second model we assume a linearly increasing radial velocity component with an added poloidal component of order 10 km/s at latitudes around $70^\circ$. The true velocity field is likely to be in between these two limiting cases. We also find that the expansion of the ansae is non-radial with reference to the central star. Their velocity field is focused near the apparent exit points from the main shell. We predict the proper motion pattern for the model with a non-zero poloidal velocity component.

We report the results of a Submillimeter Array (SMA) interferometric observation of 21-micron source IRAS 07134+1005 in the CO J=3-2 line. In order to determine the morpho-kinematic properties of the molecular envelope of the object, we constructed a model using the Shape software to model the observed CO map. We find that the molecular gas component of the envelopes can be interpreted as a geometrically thick expanding torus with an expanding velocity of 8 km/s. The inner and outer radii of the torus determined by fitting Shape models are 1.2" and 3.0", respectively. The inner radius is consistent with the previous values determined by radiative transfer modeling of the spectral energy distribution and mid-infrared imaging of the dust component. The radii and expansion velocity of the torus suggest that the central star has left the asymptotic giant branch about 1140-1710 years ago, and that the duration of the equatorial enhanced mass loss is about 2560-3130 years. From the absence of an observed jet, we suggest that the formation of a bipolar outflow may lack behind in time from the creation of the equatorial torus.

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

We report on the results of a Submillimeter Array interferometric observation of the proto-planetary nebula CRL 618 in the 12CO J=6-5 line. With the new capability of SMA enabling us to use two receivers at a time, we also observed simultaneously in the 12CO J=2-1 and 13CO J=2-1 lines. The 12CO J=6-5 and 13CO J=2-1 lines were first interferometrically observed toward CRL 618. The flux of the high velocity component of the 12CO J=6-5 line is almost fully recovered, while roughly 80% of the flux of the low velocity component is resolved out. The low recovery rate suggests that the emission region of the low velocity component of the 12CO J=6-5 line is largely extended. Continuum emission is detected both at 230 and 690 GHz. The flux of the 690 GHz continuum emission seems to be partially resolved out, suggesting dust emission partly contaminates the 690 GHz continuum flux. The cavity structure, which has been confirmed in a previous observation in the 12CO J=2-1 line, is not clearly detected in the 12CO J=6-5 line, and only the south wall of the cavity is detected. This result suggests that the physical condition of the molecular envelope of CRL 618 is not exactly axial symmetric.

Using the 3.6/4.5/5.8/8.0 micron images with 1.2 arcsec pixel resolution from the Spitzer/GLIMPSE survey, we investigated 23 masing and 18 very red objects that were not identified in the 2MASS survey. Counterparts for all selected objects were found in the GLIMPSE images. Color indices in these IR bands suggest the presence of a high-extinction layer of more than a few tenths of a solar mass in front of the central star. Furthermore, radio observations in the SiO and H2O maser lines found characteristic maser-line spectra of the embedded objects, e.g., the SiO J=1-0 line intensity in the v=2 state stronger than that of the v=1 state, or very widespread H2O maser emission spectra. This indicates that these objects are actually enshrouded by very thick circumstellar matter, some of which cannot be ascribed to the AGB wind of the central star. Individually interesting objects are discussed, including two newly found water fountains and an SiO source with nebulosity.