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.

Hadron-quark combustion in dense matter is a central topic in the study of phases in compact stars and their high-energy astrophysics. We critically review the literature on hadron-quark combustion, dividing them into a "first wave" that treats the problem as a steady-state burning with or without constraints of mechanical equilibrium, and a "second wave" which uses numerical techniques to resolve the burning front and solves the underlying Partial Differential Equations for the chemistry of the burning front under less restrictive conditions. We detail the inaccuracies that the second wave amends over the first wave, and highlight crucial differences between various approaches in the second wave. We also include results from time-dependent simulations of the reaction zone that include a hadronic EOS, neutrinos, and self-consistent thermodynamics without using parameterized shortcuts.

The hadron-quark combustion front is a system that couples various processes, such as chemical reactions, hydrodynamics, diffusion, and neutrino transport. Previous numerical work has shown that this system is very nonlinear, and can be very sensitive to some of these processes. In these proceedings, we contextualize the hadron-quark combustion as a nonlinear system, subject to dramatic feedback triggered by leptonic weak decays and neutrino transport.

We simulate and study the microphysics of combustion (flame burning) of two flavored quark matter (u, d) to three flavoured quark matter (u,d,s) in a trapped neutrino regime applicable to conditions prevailing in a hot proto-neutron star. The reaction-diffusion-advection equations for (u,d) to (u,d,s) combustion are coupled with neutrino transport, which is modelled through a flux-limited diffusion scheme. The flame speed is proportional to initial lepton fraction because of the release of electron chemical potential as heat, and reaches a steady-state burning speed of (0.001-0.008)c. We find that the burning speed is ultimately driven by the neutrino pressure gradient, given that the pressure gradient induced by quarks is opposed by the pressure gradients induced by electrons. This suggests, somewhat counter-intuitively, that the pressure gradients that drive the interface are controlled primarily by leptonic weak decays rather than by the quark Equation of State (EOS). In other words, the effects of the leptonic weak interaction, including the corresponding weak decay rates and the EOS of electrons and neutrinos, are at least as important as the uncertainties related to the EOS of high density matter. We find that for baryon number densities nB <= 0.35 fm-3, strong pressure gradients induced by leptonic weak decays drastically slow down the burning speed, which is thereafter controlled by the much slower burning process driven by backflowing downstream matter. We discuss the implications of our findings to proto-neutron stars.

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.

Ouyed et al. (1998) proposed Deuterium (DD) fusion at the core-mantle interface of giant planets as a mechanism to explain their observed heat excess. But rather high interior temperatures (~10^5 K) and a stratified D layer are needed, making such a scenario unlikely. In this paper, we re-examine DD fusion, with the addition of screening effects pertinent to a deuterated core containing ice and some heavy elements. This alleviates the extreme temperature constraint and removes the requirement of a stratified D layer. As an application, we propose that, if their core temperatures are a few times 10^4 K and core composition is chemically inhomogeneous, the observed inflated size of some giant exoplanets ("hot Jupiters") may be linked to screened DD fusion occurring deep in the interior. Application of an analytic evolution model suggests that the amount of inflation from this effect can be important if there is sufficient rock-ice in the core, making DD fusion an effective extra internal energy source for radius inflation. The mechanism of screened DD fusion, operating in the above temperature range, is generally consistent with the trend in radius anomaly with planetary equilibrium temperature $T_{\rm eq}$, and also depends on planetary mass. Although we do not consider the effect of incident stellar flux, we expect that a minimum level of irradiation is necessary to trigger core erosion and subsequent DD fusion inside the planet. Since DD fusion is quite sensitive to the screening potential inferred from laboratory experiments, observations of inflated hot Jupiters may help constrain screening effects in the cores of giant planets.

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.

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.

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.

Current r-process models under-produce A<130 nuclei. P-process models either underproduce A<100 and 150<A<165 p-nuclei, or have other limitations. We argue that the puzzles in these two mechanisms can be solved by heavy-ion spallation. Spallation happens in the explosive phase transition from a neutron star to a quark star - the quark-nova. The quark-nova triggers the r-process which creates abundant A>130 nuclei, and spallation fragments these isotopes into A<130 nuclei and all 35 p-nuclei. Our model is universal in relation to a star's age, metallicity, and chemistry.

In 1982, Monique and Francois Spite discovered that the 7Li abundance in the atmosphere of old metal-poor dwarf stars in the galactic halo was independent of metallicity and temperature. Since then, 7Li abundance in the Universe has become a subject of intrigue, because there is less of it in Population II dwarf stars (by a factor of 3) than standard big bang nucleosynthesis predicts. Here we show how quark-novae (QNe) occurring in the wake of Pop. III stars, can elegantly produce an A(Li) ~ 2.2 Lithium plateau in Pop. II (low-mass) stars formed in the pristine cloud swept up by the mixed SN+QN ejecta. We also find an increase in the scatter as well as an eventual drop in A(Li) below the Spite plateau values for very low metallicity ([Fe/H] < -3) in excellent agreement with observations. We propose a solution to the discrepancy between the Big Bang Nucleosynthesis 7Li abundance and the Spite plateau and list some implications and predictions of our model.

[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/

[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/

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 review aspects of the hydrodynamical combustion of nuclear matter to strange quark matter in a neutron star. Numerical studies on non-premixed combustion that consistently include hydrodynamical flows in a reactive-diffusive setup show that in 1D, the conversion (burning) front moves at sub-sonic speeds. Although the front eventually comes to a halt (in the star's reference frame), it still has ud matter falling onto it, and so is capable of converting the rest of the star. However, in the process, we also find that neutrino cooling of the interface causes it to wrinkle, laying a platform for a deflagrative-to-detonative transition (DDT). We outline progress on improvements in the burning code (Burn-UD: http://quarknova.ucalgary.ca/software/Burn-UD/) that will ultimately reveal the mechanism that can explode the outermost layers of even a dense compact object like a neutron star.

Super-luminous supernova (SLSN) are supernovae showing extreme properties in their light-curves: high peak luminosities (more than 10 times brighter than bright SN Ia), and long durations. Several mechanisms have been proposed for SLSN, such as pair instability SN of a massive progenitor, interaction of the ejecta with a massive circumstellar shell, and the dual-shock quark nova (dsQN) model. The dual-shock quark nova model is unique in that it predicts a normal SN event will be seen about 10 days prior to the main SLSN event. The dsQN model is described here and shown that it is consistent with the light curve of the one currently known double-humped SLSN, 2006oz.

We show that appealing to a Quark-Nova (QN) in a tight binary system containing a massive neutron star and a CO white dwarf (WD), a Type Ia explosion could occur. The QN ejecta collides with the WD driving a shock that triggers Carbon burning under degenerate conditions (the QN-Ia). The conditions in the compressed low-mass WD (M_WD < 0.9M_sun) in our model mimics those of a Chandrasekhar mass WD. The spin-down luminosity from the QN compact remnant (the quark star) provides additional power that makes the QN-Ia light-curve brighter and broader than a standard SN-Ia with similar 56Ni yield. In QNe-Ia, photometry and spectroscopy are not necessarily linked since the kinetic energy of the ejecta has a contribution from spin-down power and nuclear decay. Although QNe-Ia may not obey the Phillips relationship, their brightness and their relatively "normal looking" light-curves means they could be included in the cosmological sample. Light-curve fitters would be confused by the discrepancy between spectroscopy at peak and photometry and would correct for it by effectively brightening or dimming the QNe-Ia apparent magnitudes. Thus over- or under-estimating the true magnitude of these spin-down powered SNe-Ia. Contamination of QNe-Ia in samples of SNe-Ia used for cosmological analyses could systematically bias measurements of cosmological parameters if QNe-Ia are numerous enough at high-redshift. The strong mixing induced by spin-down wind combined with the low 56Ni yields in QNe-Ia means that these would lack a secondary maximum in the i-band despite their luminous nature. We discuss possible QNe-Ia progenitors.

The genesis and chemical patterns of the metal poor stars in the galactic halo remains an open question. Current models do not seem to give a satisfactory explanation for the observed abundances of Lithium in the galactic metal-poor stars and the existence of carbon-enhanced metal-poor (CEMP) and Nitrogen-enhanced metal-poor (NEMP) stars. In order to deal with some of these theoretical issues, we suggest an alternative explanation, where some of the Pop. III SNe are followed by the detonation of their neutron stars (Quark-Novae; QNe). In QNe occurring a few days to a few weeks following the preceding SN explosion, the neutron-rich relativistic QN ejecta leads to spallation of 56Ni processed in the ejecta of the preceding SN explosion and thus to "iron/metal impoverishment" of the primordial gas swept by the combined SN+QN ejecta. We show that the generation of stars formed from fragmentation of pristine clouds swept-up by the combined SN+QN ejecta acquire a metallicity with -7.5 < Fe/H] < -1.5 for dual explosions with 2 < t_delay (days) < 30. Spallation leads to the depletion of 56Ni and formation of sub-Ni elements such as Ti, V, Cr, and Mn providing a reasonable account of the trends observed in galactic halo metal-poor stars. CEMP stars form in dual explosions with short delays (t_delay < 5 days). These lead to important destruction of 56Ni (and thus to a drastic reduction of the amount of Fe in the swept up cloud) while preserving the carbon processed in the outer layers of the SN ejecta. Lithium is produced from the interaction of the neutron-rich QN ejecta with the outer (oxygen-rich) layers of the SN ejecta. A Lithium plateau with 2 < A(Li) < 2.4 can be produced in our model as well as a corresponding 6Li plateau with 6Li/7Li < 0.3.

Current cc-SN models predict overproduction of 44Ti compared to observations. We present a model for an alternative channel where a cc-SN explosion is followed by a neutron star detonation (Quark Nova or QN), resulting in a spallation reaction of SN ejecta that produces 44Ti. We can achieve a 44Ti production of ~ 10^-4 Msun with our model under the right time delay between the QN and the SN. Our model also produces unique signals not found in standard, cc-SN nucleosynthesis models. Some of these unique signals include a significantly large production of 7Be and 22Na. We discuss some of these signals by analyzing the late time light curve and gamma spectroscopy of our model.

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.

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 study the spindown of isolated neutron stars from initially rapid rotation rates, driven by two factors: (i) gravitational wave emission due to r-modes and (ii) magnetic braking. In the context of isolated neutron stars, we present the first study including self-consistently the magnetic damping of r-modes in the spin evolution. We track the spin evolution employing the RNS code, which accounts for the rotating structure of neutron stars for various equations of state. We find that, despite the strong damping due to the magnetic field, r-modes alter the braking rate from pure magnetic braking for B<10^{13}G. For realistic values of the saturation amplitude, the r-mode can also decrease the time to reach the threshold central density for quark deconfinement. Within a phenomenological model, we assess the gravitational waveform that would result from r-mode driven spindown of a magnetized neutron star. To contrast with the persistent signal during the spindown phase, we also present a preliminary estimate of the transient gravitational wave signal from an explosive quark-hadron phase transition, which can be a signal for the deconfinement of quarks inside neutron stars.

SN2006oz is a super-luminous supernova with a mysterious bright precursor that has resisted explanation in standard models. However, such a precursor has been predicted in the dual-shock quark nova (dsQN) model of super-luminous supernovae -- the precursor is the SN event while the main light curve of the SLSN is powered by the Quark-Nova (QN; explosive transition of the neutron star to a quark star). As the SN is fading, the QN re-energizes the SN ejecta, producing a "double-humped" light curve. In this paper, we show the dsQN model successfully reproduces the observed light curve of SN2006oz.

Titanium-rich subluminous supernovae are rare and challenge current SN nucleosynthesis models. We present a model in which ejecta from a standard Supernova is impacted by a second explosion of the neutron star (a Quark-nova), resulting in spallation reactions that lead to 56Ni destruction and 44Ti creation under the right conditions. Basic calculations of the spallation products shows that a delay between the two explosions of ~ 5 days reproduces the observed abundance of 44Ti in Cas A and explains its low luminosity as a result of the destruction of 56Ni. Our results could have important implications for lightcurves of subluminous as well as superluminous supernovae.

We propose a simple model explaining two outstanding astrophysical problems related to compact objects: (1) that of stars such as G87-7 (alias EG 50) that constitute a class of relatively low-mass white dwarfs which nevertheless fall away from the C/O composition and (2) that of GRB 110328A/Swift J164449.3+57345 which showed spectacularly long-lived strong X-ray flaring, posing a challenge to standard GRB models. We argue that both these observations may have an explanation within the unified framework of a Quark-Nova occurring in a low-mass X-ray binary (neutron star- white dwarf). For LMXBs where the binary separation is sufficiently tight, ejecta from the exploding Neutron Star triggers nuclear burning in the white dwarf on impact, possibly leading to Fe-rich composition compact white dwarfs with mass 0.43M_sun < M_WD < 0.72M_sun, reminiscent of G87-7. Our results rely on the assumption, which ultimately needs to be tested by hydrodynamic and nucleosynthesis simulations, that under certain circumstances the WD can avoid the thermonuclear runaway. For heavier white dwarfs (i.e. M_WD > 0.72M_sun) experiencing the QN shock, degeneracy will not be lifted when Carbon burning begins, and a sub-Chandrasekhar Type Ia Supernovae may result in our model. Under slightly different conditions, and for pure He white dwarfs (i.e. M_WD < 0.43M_sun), the white dwarf is ablated and its ashes raining down on the Quark star leads to accretion-driven X-ray luminosity with energetics and duration reminiscent of GRB 110328A. We predict additional flaring activity towards the end of the accretion phase if the Quark star turns into a Black Hole.

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/

Soft gamma repeaters and anomalous X-ray pulsars are believed to be magnetars, i.e. neutron stars powered by extreme magnetic fields, B~10^(14)-10^(15) Gauss. The recent discovery of a soft gamma repeater with low magnetic field (< 7.5x10^(12) Gauss), SGR 0418+5729, which shows bursts similar to those of SGRs, implies that a high surface dipolar magnetic field might not be necessary for magnetar-like activity. We show that the quiescent and bursting properties of SGR 0418+5729 find natural explanations in the context of low-magnetic field Quark-Nova (detonative transition from a neutron star to a quark star) remnants, i.e. an old quark star surrounded by degenerate (iron-rich) Keplerian ring/debris ejected during the Quark-Nova explosion. We find that a 16 Myr old quark star surrounded by a ~ 10^(-10)xM_sun ring, extending in radius from ~ 30 km to 60 km, reproduces many observed properties of SGR 0418+5729. The SGR-like burst is caused by magnetic penetration of the inner part of the ring and subsequent accretion. Radiation feedback results in months-long accretion from the ring's non-degenerate atmosphere which matches well the observed decay phase. We make specific predictions (such as an accretion glitch of Delta P/P ~ - 2x10^(-11) during burst and a sub-keV proton cyclotron line from the ring) that can be tested by sensitive observations.

We show that several features reminiscent of short-hard Gamma-ray Bursts (GRBs) arise naturally when Quark-Novae occur in low-mass X-ray binaries born with massive neutron stars (> 1.6M_sun) and harboring a circumbinary disk. Near the end of the first accretion phase, conditions are just right for the explosive conversion of the neutron star to a quark star (Quark-Nova). In our model, the subsequent interaction of material from the neutron star's ejected crust with the circumbinary disk explains the duration, variability and near-universal nature of the prompt emission in short-hard GRBs. We also describe a statistical approach to ejecta break-up and collision to obtain the photon spectrum in our model, which turns out remarkably similar to the empirical Band function (Band 1993). We apply the model to the fluence and spectrum of GRB 000727, GRB 000218, and GRB980706A obtaining excellent fits. Extended emission (spectrum and duration) is explained by shock-heating and ablation of the white dwarf by the highly energetic ejecta. Depending on the orbital separation when the Quark-Nova occurs, we isolate interesting regimes within our model when both prompt and extended emission can occur. We find that the spectrum can carry signatures typical of Type Ib/c SNe although these should appear less luminous than normal type Ib/c SNe. Late X-ray activity is due to accretion onto the quark star as well as its spin-down luminosity. Afterglow activity arise from the expanding shell of material from the shock-heated expanding circumbinary disk. We find a correlation between the duration and spectrum of short-hard GRBs as well as modest hard-to-soft time evolution of the peak energy.

We present results from a numerical solution to the burning of neutron matter inside a cold neutron star into stable (u,d,s) quark matter. Our method solves hydrodynamical flow equations in 1D with neutrino emission from weak equilibrating reactions, and strange quark diffusion across the burning front. We also include entropy change due to heat released in forming the stable quark phase. Our numerical results suggest burning front laminar speeds of 0.002-0.04 times the speed of light, much faster than previous estimates derived using only a reactive-diffusive description. Analytic solutions to hydrodynamical jump conditions with a temperature dependent equation of state agree very well with our numerical findings for fluid velocities. The most important effect of neutrino cooling is that the conversion front stalls at lower density (below approximately 2 times saturation density). In a 2-dimensional setting, such rapid speeds and neutrino cooling may allow for a flame wrinkle instability to develop, possibly leading to detonation.

We perform state-of-the-art, 3D, time-dependent simulations of magnetized disk winds, carried out to simulation scales of 60 Astronomical Units, in order to confront optical HST observations of protostellar jets. We ``observe'' the optical forbidden line emission produced by shocks within our simulated jets and compare these with actual observations. Our simulations reproduce the rich structure of time varying jets, including jet rotation far from the source, an inner (up to 400 km/s) and outer (less than 100 km/s) component of the jet, and jet widths of up to 20 Astronomical Units in agreement with observed jets. These simulations when compared with the data are able to constrain disk wind models. In particular, models featuring a disk magnetic field with a modest radial spatial variation across the disk are favored.

Compact stars made of quark matter rather than confined hadronic matter, are expected to form a color superconductor. This superconductor ought to be threaded with rotational vortex lines within which the star's interior magnetic field is confined. The vortices (and thus magnetic flux) would be expelled from the star during stellar spin-down, leading to magnetic reconnection at the surface of the star and the prolific production of thermal energy. In this Letter, we show that this energy release can re-heat quark stars to exceptionally high temperatures, such as observed for Soft Gamma Repeaters (SGRs), Anomalous X-Ray pulsars (AXPs), and X-ray dim isolated neutron stars (XDINs). Moreover, our numerical investigations of the temperature evolution, spin-down rate, and magnetic field behavior of such superconducting quark stars suggest that SGRs, AXPs, and XDINs may be linked ancestrally. Finally, we discuss the possibility of a time delay before the star enters the color superconducting phase, which can be used to estimate the density at which quarks deconfine. We find this density to be five times that of nuclear saturation.

In this paper we consider a quark nova occurring inside an exploding star. The quark nova ejecta will shock when interacting with the stellar envelope. When this shock reaches the surface of the star, the energy is radiated away. We suggest that this energy may be seen in X-rays, and show here that this may explain some flares seen in the X-ray afterglow of long gamma ray bursts (GRBs). A quark nova inside an exploding star need not be followed by a GRB, or the GRB may not be beamed towards us. However, the shock breakout is likely not beamed and could be seen even in the absence of a GRB. We suggest that XRO 080109 is such an event in which a quark nova occurs inside an exploding star. No GRB is formed, but the break out of the shock leads to the XRO.

(Abridged) XTE J1810-197 and 1E 1547.0-5408 are two transient AXPs exhibiting radio emission with unusual properties. In addition, their spin down rates during outburst show opposite trends, which so far has no explanation. Here, we extend our quark-nova model for AXPs to include transient AXPs, in which the outbursts are caused by transient accretion events from a Keplerian (iron-rich) degenerate ring. For a ring with inner and outer radii of 23.5 km and 26.5 km, respectively, our model gives a good fit to the observed X-ray outburst from XTE J1810-197 and the behavior of temperature, luminosity, and area of the two X-ray blackbodies with time. The two blackbodies in our model are related to a heat front (i.e. Bohm diffusion front) propagating along the ring's surface and an accretion hot spot on the quark star surface. Radio pulsations in our model are caused by dissipation at the light cylinder of magnetic bubbles, produced near the ring during the X-ray outburst. The delay between X-ray peak emission and radio emission in our model is related to the propagation time of these bubbles to the light cylinder. We predict a ~1 year and ~1 month delay for XTE J1810-197 and 1E 1547.0-5408, respectively. The observed flat spectrum, erratic pulse profile, and the pulse duration are all explained in our model as a result of X-point reconnection events induced by the dissipation of the bubbles at the light cylinder. The spin down rate of the central quark star can either increase or decrease depending on how the radial drift velocity of the magnetic islands changes with distance from the central star. We suggest an evolutionary connection between transient AXPs and typical AXPs in our model.

Compact stars made of quark matter rather than confined hadronic matter, are expected to form a color superconductor. This superconductor ought to be threaded with rotational vortex lines, within which the star's interior magnetic field is at least partially confined. The vortices (and thus magnetic flux) would be expelled from the star during stellar spin-down, leading to magnetic reconnection at the surface of the star and the prolific production of thermal energy. In this paper, we show that this energy release can re-heat quark stars to exceptionally high temperatures, such as observed for Soft Gamma Repeaters (SGRs), Anomalous X-Ray pulsars (AXPs), and X-ray dim isolated neutron stars (XDINs). Moreover, our numerical investigations of the temperature evolution, spin-down rate, and magnetic field behavior of such superconducting quark stars suggest that SGRs, AXPs, and XDINs may be linked ancestrally. Finally, we discuss the possibility of a time delay before the star enters the color superconducting phase, which can be used to estimate the density at which quarks deconfine. From observations, we find this density to be of the order of five times that of nuclear saturation.

[Abridged] Superluminous Supernovae (SN2006gy, SN2005gj, SN2005ap, SN2008fz, SN2003ma) have been a challenge to explain by standard models. We present an alternative scenario involving a quark-nova (QN), an explosive transition of the newly born neutron star to a quark star in which a second explosion (delayed) occurs inside the already expanding ejecta of a normal SN. The reheated SN ejecta can radiate at higher levels for longer periods of time primarily due to reduced adiabatic expansion losses, unlike the standard SN case. Our model is successfully applied to SN2006gy, SN2005gj, SN2005ap, SN2008fz, SN2003ma with encouraging fits to the lightcurves. There are four predictions in our model: (i) superluminous SNe optical lightcurves should show a double-hump with the SN hump at weaker magnitudes occurring days to weeks before the QN; (ii) Two shock breakouts should be observed vis-a-vis one for a normal SN. Depending on the time delay, this would manifest as two distinct spikes in the X-ray region or a broadening of the first spike for extremely short delays; (iii) The QN deposits heavy elements of mass number A> 130 at the base of the preceeding SN ejecta. These QN r-processed elements should be visible in the late spectrum (few days-weeks in case of strong ejecta mixing) of the superluminous SN; (iv) The QN yield will also contain lighter elements (Hydrogen and Helium). We expect the late spectra to include H_alpha emission lines that should be distinct in their velocity signature from standard H_alpha emission.

We examine the case for Quark-Novae (QNe) as possible sources for the reionization and early metal enrichment of the universe. Quark-Novae are predicted to arise from the explosive collapse (and conversion) of sufficiently massive neutron stars into quark stars. A Quark-Nova (QN) can occur over a range of time scales following the supernova event. For QNe that arise days to weeks after the supernovae, we show that dual-shock that arises as the QN ejecta encounter the supernova ejecta can produce enough photons to reionize hydrogen in most of the Inter-Galactic medium (IGM) by z ~ 6. Such events can explain the large optical depth tau_e ~ 0.1 as measured by WMAP, if the clumping factor, C, of the material being ionized is smaller than 10. We suggest a way in which a normal initial mass function (IMF) for the oldest stars can be reconciled with a large optical depth as well as the mean metallicity of the early IGM post reionization. We find that QN also make a contribution to r-process element abundances for atomic numbers A > 130. We predict that the main cosmological signatures of Quark-Novae are the gamma-ray bursts that announce their birth. These will be clustered at redshifts in the range z ~ 7-8 in our model.

We present three-dimensional magnetohydrodynamic (MHD) simulations of superbubbles, to study the importance of MHD effects in the interpretation of images from recent surveys of the Galactic plane. These simulations focus mainly on atmospheres defined by an exponential density distribution and the Dickey & Lockman (1990) density distribution. In each case, the magnetic field is parallel to the Galactic plane and we investigate cases with either infinite scale height (constant magnetic field) or a constant ratio of gas pressure to magnetic pressure. The three-dimensional structure of superbubbles in these simulations is discussed with emphasis on the axial ratio of the cavity as a function of magnetic field strength and the age of the bubble. We investigate systematic errors in the age of the bubble and scale height of the surrounding medium that may be introduced by modeling the data with purely hydrodynamic models. Age estimates derived with symmetric hydrodynamic models fitted to an asymmetric magnetized superbubble can differ by up to a factor of four, depending on the direction of the line of sight. The scale height of the surrounding medium based on the Kompaneets model may be up to 50% lower than the actual scale height. We also present the first ever predictions of Faraday rotation by a magnetized superbubble based on three-dimensional MHD simulations. We emphasize the importance of MHD effects in the interpretation of observations of superbubbles.

If a quark-nova occurs inside a collapsar, the interaction between the quark-nova ejecta (relativistic iron-rich chunks) and the collapsar envelope, leads to features indicative of those observed in Gamma Ray Bursts. The quark-nova ejecta collides with the stellar envelope creating an outward moving cap (Gamma ~ 1-10) above the polar funnel. Prompt gamma-ray burst emission from internal shocks in relativistic jets (following accretion onto the quark star) become visible after the cap becomes optically thin. Model features include: (i) precursor activity (optical, X-ray, gamma-ray), (ii) prompt gamma-ray emission, and (iii) afterglow emission. We discuss SN-less long duration GRBs, short hard GRBs (including association and non-association with star forming regions), dark GRBs, the energetic X-ray flares detected in Swift GRBs, and the near-simultaneous optical and gamma-ray prompt emission observed in GRBs in the context of our model.

We explore the dynamical and thermal evolution of the ejected neutron star crust in a Quark-Nova explosion. Typical explosion energies and ejected crust masses result in relativistic ejection with Lorentz factors of a few to a few hundred. The ejecta undergoes a rapid cooling and stretching resulting in break up into many small pieces (clumps) when the ejecta is only ~ 100 km from the explosion site. The number and size of the clumps depends on whether the breakup occurs in the liquid or solid phase. For these two cases, the clump number is ~ 10^3 (liquid phase) or ~ 10^7 (solid phase) and, at break up, are spherical (size ~ 10^4 cm; liquid phase) or needle shaped (~ 10^4x10^2 cm; solid phase).

Astrophysical jets are associated with the formation of young stars of all masses, stellar and massive black holes, and perhaps even with the formation of massive planets. Their role in the formation of planets, stars, and galaxies is increasingly appreciated and probably reflects a deep connection between the accretion flows - by which stars and black holes may be formed - and the efficiency by which magnetic torques can remove angular momentum from such flows. We compare the properties and physics of jets in both non-relativistic and relativistic systems and trace, by means of theoretical argument and numerical simulations, the physical connections between these different phenomena. We discuss the properties of jets from young stars and black holes, give some basic theoretical results that underpin the origin of jets in these systems, and then show results of recent simulations on jet production in collapsing star-forming cores as well as from jets around rotating Kerr black holes.

Magnetic reconnection plays a critical role in many astrophysical processes where high energy emission is observed, e.g. particle acceleration, relativistic accretion powered outflows, pulsar winds and probably in dissipation of Poynting flux in GRBs. The magnetic field acts as a reservoir of energy and can dissipate its energy to thermal and kinetic energy via the tearing mode instability. We have performed 3d nonlinear MHD simulations of the tearing mode instability in a current sheet. Results from a temporal stability analysis in both the linear regime and weakly nonlinear (Rutherford) regime are compared to the numerical simulations. We observe magnetic island formation, island merging and oscillation once the instability has saturated. The growth in the linear regime is exponential in agreement with linear theory. In the second, Rutherford regime the island width grows linearly with time. We find that thermal energy produced in the current sheet strongly dominates the kinetic energy. Finally preliminary analysis indicates a P(k) 4.8 power law for the power spectral density which suggests that the tearing mode vortices play a role in setting up an energy cascade.

The most luminous Supernova SN2006gy (more than a 100 times brighter than a typical supernova) has been a challenge to explain by standard models. For example, pair instability supernovae which are luminous enough seem to have too slow a rise, and core collapse supernovae do not seem to be luminous enough. We present an alternative scenario involving the quark-nova phenomenon (an explosive transition of the newly born neutron star to a quark star) in which a second explosion (delayed) occurs inside the ejecta of a normal supernova. The reheated supernova ejecta can radiate at higher levels for longer periods of time primarily due to reduced adiabatic expansion losses, unlike the standard supernova case. We find an encouraging match between the resulting lightcurve and that observed in the case of SN2006gy suggesting that we might have at hand the first ever signature of a quark-nova. Successful application of our model to SN2005gj and SN2005ap is also presented.

This is the third paper of a series of papers where we explore the evolution of iron-rich ejecta from quark-novae. In the first paper, we explored the case where a quark-nova ejecta forms a degenerate shell, supported by the star's magnetic field, with applications to SGRs. In the second paper we considered the case where the ejecta would have sufficient angular momentum to form a degenerate Keplerian torus and applied such a system to two AXPs, namely 1E2259 + 586 and 4U0142 + 615. Here we explore the late evolution of the degenerate torus and find that it can remain unchanged for $\sim 10^{6}$ years before it becomes non-degenerate. This transition from a degenerate torus (accretion dominated) to a non-degenerate disk (no accretion), occurs about $10^{6}$ years following the quark-nova, and exhibits features that are reminiscent of observed properties of RRATs. Using this model, we can account for the duration of both the radio bursts and the quiet phase, as well as the observed radio flux from RRATs. The unique on and off activity of the radio pulsar PSR B1931 + 24 is similar to that of "old RRATs" in our model. For old RRATs, in our model, the spin-down rate during the radio loud phase is about a factor 1.6 larger than when it is quiet, remarkably similar to what has been measured for PSR B1931 + 24. We discuss a connection between XDINs and RRATs and argue that some XDINs may be ``dead RRATs'' that have already consumed their non-degenerate disk.

We examine neutron star properties based on a model of dense matter composed of B=1 skyrmions immersed in a mesonic mean field background. The model realizes spontaneous chiral symmetry breaking non-linearly and incorporates scale-breaking of QCD through a dilaton VEV that also affects the mean fields. Quartic self-interactions among the vector mesons are introduced on grounds of naturalness in the corresponding effective field theory. Within a plausible range of the quartic couplings, the model generates neutron star masses and radii that are consistent with a preponderance of observational constraints, including recent ones that point to the existence of relatively massive neutron stars with mass M 1.7 Msun and radius R (12-14) km. If the existence of neutron stars with such dimensions is confirmed, matter at supra-nuclear density is stiffer than extrapolations of most microscopic models suggest.

The effects of temperature on strange stars are studied and it is found that the maximum mass of the star decreases with the increase of temperature since at high temperatures the equations of state become softer. Moreover, if the temperature of a strange star increases, keeping its baryon number fixed, its gravitational mass increases and radius decreases. This leads to a limiting temperature where it turns into a black hole. These features are due to a combined effect of the change of gluon mass and the quark distribution with temperature. We report a new kind of radial oscillations of strange stars driven by what we call chromo-thermal instability. We also discuss the relevance of our findings in the astrophysics of core collapse supernovae and gamma ray bursts.

We describe a model within the ``Quark-nova'' scenario to interpret the recent observations of early X-ray afterglows of long Gamma-Ray Bursts (GRB) with the Swift satellite. This is a three-stage model within the context of a core-collapse supernova. STAGE 1 is an accreting (proto-) neutron star leading to a possible delay between the core collapse and the GRB. STAGE 2 is accretion onto a quark-star, launching an ultrarelativistic jet generating the prompt GRB. This jet also creates the afterglow as the jet interacts with the surrounding medium creating an external shock. Slower shells ejected from the quark star (during accretion), can re-energize the external shock leading to a flatter segment in the X-ray afterglow. STAGE 3, which occurs only if the quark-star collapses to form a black-hole, consists of an accreting black-hole. The jet launched in this accretion process interacts with the preceding quark star jet, and could generate the flaring activity frequently seen in early X-ray afterglows. Alternatively, a STAGE 2b can occur in our model if the quark star does not collapse to a black hole. The quark star in this case can then spin down due to magnetic braking, and the spin down energy may lead to flattening in the X-ray afterglow as well. This model seems to account for both the energies and the timescales of GRBs, in addition to the newly discovered early X-ray afterglow features.

We present new determination of the birth rate of AXPs and SGRS and their associated SNRs. We find a high birth rate of 1/(500 yr) for AXPs/SGRs and 1/(1700 yr) for associated SNRs. These high rates suggest that all massive stars (greater than ~ 25 M_sun) give rise to remnants with magnetar-like fields. Recent observations indicate that fossil fields cannot explain such high fields in the progenitor stars. Dynamo mechanisms during the birth of the neutron stars require spin rates much faster than either observations or theory indicate. Here, we propose the neutron stars form with normal (~ 10^{12} G) magnetic fields, which are then amplified to 10^{14}-10^{15} G after a delay of a few hundred years. The amplification is speculated to be a consequence of color ferromagnetism and occurs after the neutron star core reaches quark-deconfinement density. This delayed amplification alleviates many difficulties in interpreting simultaneously the high birth rate and high magnetic fields of AXPs/SGRs and their link to massive stars.

We present a three-stage model for a long GRB inner engine to explain the prompt gamma ray emission, and interpret recent Swift satellite observations of early X-ray afterglow plateaus followed by a sharp drop off or a shallow power law decay. The three stages involves a neutron star phase, a quark star (QS) and a black hole phase as described in Staff et al. (2007). We find that the QS stage allows for more energy to be extracted from neutron star to QS conversion as well as from ensuing accretion onto the QS. The QS accretion phase naturally extends the engine activity and can account for both the prompt emission and irregular early X-ray afterglow activity. Following the accretion phase, the QS can spin-down by emission of a baryon-free outflow. The magnetar-like magnetic field strengths resulting from the NS to QS transition provide enough spin-down energy, for the correct amount of time, to account for the plateau in the X-ray afterglow. In our model, a sharp drop-off following the plateau occurs when the QS collapses to a BH during the spin-down, thus shutting-off the secondary outflow. We applied our model to GRB 070110 and GRB 060607A and found that we can consistently account for the energetics and duration during the prompt and plateau phases.

In the first paper of this series we explored the case where a quark-nova ejecta forms a degenerate shell, supported by the star's magnetic field. Herein, we consider the case where the ejecta has sufficient angular momentum to form a Keplerian torus, and we show that the density and temperature of the torus are such that it will remain degenerate throughout it's lifetime. We go on to discuss the evolution of such a torus and apply it to two AXPs, namely 1E2259$+$586 and 4U0142$+$615. As it turns out, using our model we can account for many of the observations of these objects including the quiescent phase luminosity, and blackbody temperatures during quiescence and bursting phases. Furthermore, for 1E2259$+$586 our model explains the steep and slow decay components seen in the burst lightcurve, as well as the rotation period glitches and enhanced spin-down rate. We also estimate the mass of the degenerate torus to be of the order of $10^{-6}M_{\odot}$, and speculate that the observed optical/infrared emission from 4U0142$+$615 might be a signature of the thin degenerate torus we describe here.

We explore the formation and evolution of debris ejected around quark stars in the Quark Nova scenario, and the application to Soft Gamma-ray Repeaters (SGRs) and Anomolous X-ray Pulsars (AXPs). If an isolated neutron star explodes as a Quark Nova, an Iron-rich shell of degenerate matter forms out of the fall-back (crust) material. Our model can account for many of the observed features of SGRs and AXPs such as: (i) the two types of bursts (giant and regular); (ii) the spin-up and spin-down episodes during and following the bursts with associated persistant increases in $\dot{P}$; (iii) the energetics of the boxing day burst, SGR1806$+$20; (iv) the presence of an Iron line as observed in SGR1900$+$14; (v) the correlation between the far-Infrared and the X-ray fluxes during the bursting episode and the quiescent phase; (vi) the hard X-ray component observed in SGRs during the giant bursts, and (vii) the discrepancy between the ages of SGRs/AXPs and their supernova remnants. We also find a natural evolutionary relationship between SGRs and AXPs in our model which predicts that only the youngest SGRs/AXPs are most likely to exhibit strong bursting. Many features of X-ray Dim Isolated Neutron stars (XDINs) are also accounted for in our model such as, (i) the two-component blackbody spectra; (ii) the absorption lines around 300 eV; and (iii) the excess optical emission.

Relativistic outflows carrying large scale magnetic fields have large inductive potential and may accelerate protons to ultra high energies. We discuss a novel scheme of Ultra-High Energy Cosmic Ray (UHECR) acceleration due to drifts in magnetized, cylindrically collimated, sheared jets of powerful active galaxies (with jet luminosity $\geq 10^{46}$ erg s$^{-1}$). A positively charged particle carried by such a plasma is in an unstable equilibrium if ${\bf B} \cdot \nabla \times {\bf v}< 0$, so that kinetic drift along the velocity shear would lead to fast, regular energy gain. The highest rigidity particles are accelerated most efficiently implying the dominance of light nuclei for energies above the ankle in our model: from a mixed population of pre-accelerated particle the drift mechanism picks up and boosts protons preferably.

The recently suggested correlation between the surface temperature and the magnetic field in isolated neutron stars does not seem to work well for SGRs, AXPs and X-ray dim isolated neutron stars (XDINs; specifically the Magnificent Seven or M7). Instead by appealing to a Color-Flavor Locked Quark Star (CFLQS) we find a more natural explanation. In this picture, the heating is provided by magnetic flux expulsion from a crust-less superconducting quark star. Combined with our previous studies concerning the possibility of SGRs, AXPs, and XDINs as CFLQSs, this provides another piece of evidence that these objects are all related. Specifically, we propose that XDINs are the descendants of SGRs and AXPs.

We explore heavy-element nucleosynthesis by rapid neutron capture (r-process) in the decompressing ejecta from the surface of a neutron star. The decompression is triggered by a violent phase transition to strange quark matter (quark-nova scenario). The presence of neutron-rich large Z nuclei (40,95) < (Z,A) < (70,177), the large neutron-to-seed ratio, and the low electron fraction Ye ~ 0.03 in the decompressing ejecta present favorable conditions for the r-process. We perform network calculations that are adapted to the quark-nova conditions, and which mimic usual (n-gamma) equilibrium r-process calculations during the initially cold decompression phase. They match to dynamical r-process calculations at densities below neutron drip (4.10^11 g cm-3). We present results for the final element abundance distribution with and without heating from nuclear reactions, and compare to the solar abundance pattern of r-process elements. We highlight the distinguishing features of quark-novae by contrasting it with conventional nucleosynthetic sites such as type II supernovae and neutron star mergers, especially in the context of heavy-element compositions of extremely metal-deficient stars.

The highest energy cosmic rays could be produced by drifts in magnetized, cylindrically collimated, sheared jets of powerful active galaxies (i.e. FR II radiogalaxies; radio loud quasars and high power BL Lacs). We show that in such scenarios proton synchrotron radiation can give rise to detectable photon fluxes at energies ranging from hundreds of keV to tens of MeV.

We propose a new model within the ``Quark-nova'' scenario to interpret the recent observations of early afterglows of long Gamma-Ray Bursts (GRB) with the Swift satellite. This is a three-stage model within the context of a core-collapse supernova. Stage 1 is an accreting (proto-) neutron star leading to a possible delay between the core collapse and the GRB. Stage 2 is an accreting quark-star, generating the prompt GRB. Stage 3, which occurs only if the quark-star collapses to form a black-hole, consists of an accreting black-hole. The jet launched in this accretion process interacts with the ejecta from stage 2, and could generate the flaring activity frequently seen in X-ray afterglows. This model may be able to account for both the energies and the timescales of GRBs, in addition to the newly discovered early X-ray afterglow features.

(Abridged)Relativistic outflows carrying large scale magnetic fields have large inductive potential and may accelerate protons to ultra high energies. We discuss a novel scheme of Ultra-High Energy Cosmic Ray (UHECR) acceleration due to drifts in magnetized, cylindrically collimated, sheared jets of powerful active galaxies (with jet luminosity $\geq 10^{46}$ erg s$^{-1}$). We point out that a positively charged particle carried by such a plasma is in an unstable equilibrium if ${\bf B} \cdot \nabla \times {\bf v}< 0$, so that kinetic drift along the velocity shear would lead to fast, {\it regular} energy gain. We show that if a seed of pre-accelerated particles with energy below the ankle $\leq 10^{18}$ eV is present, these particles can be boosted to energies above $ 10^{19}$ eV. A key feature of the mechanism is that the highest rigidity particles are accelerated most efficiently implying the dominance of light nuclei for energies above the ankle in our model: from a mixed population of pre-accelerated particle the drift mechanism picks up and boosts protons preferably. In addition, after a particle traversed large fraction of the available potential, its Larmor radius becomes of the order of the jet thickness. In this case, {\it the maximum possible acceleration rate of inverse relativistic gyro-frequency is achieved} and a particle finally become unconfined and leave the jet.

We discuss a model wherein soft gamma-ray repeaters (SGRs), anomalous X-ray pulsars (AXPs), and radio quiet isolated neutron stars (RQINSs) are all compact objects exhibiting superconductivity, namely color-flavor locked quark stars. In particular we calculate the magnetic field decay due to the expulsion of spin-induced vortices from the star's superfluid-superconducting interior, and the resultant spin-down rate. We find that, for initial parameters characteristic of AXPs/SGRs ($10^{13}<B<10^{14}~\rm{G}$; $3<P<12~\rm{s}$), the magnetic field strengths and periods remain unchanged within a factor of two for timescales of the order of $5\times 10^{5} - 5\times 10^{7}~\rm{yrs}$ given a quark star of radius $10~\rm{km}$. Within these timescales, we show that the observed period clustering in RQINSs can be explained by compactness, as well as calculate how the magnetic field and period evolve in a manner concurrent with RQINS observations.

We study the role of spin-down in driving quark deconfinement in the high density core of isolated neutron stars. Assuming spin-down to be solely due to magnetic braking, we obtain typical timescales to quark deconfinement for neutron stars that are born with Keplerian frequencies. Employing different equations of state (EOS), we determine the minimum and maximum neutron star masses that will allow for deconfinement via spin-down only. We find that the time to reach deconfinement is strongly dependent on the magnetic field and that this time is least for EOS that support the largest minimum mass at zero spin, unless rotational effects on stellar structure are large. For a fiducial critical density of $5\rho_0$ for the transition to the quark phase ($\rho_0=2.5\times10^{14}$g/cm$^3$ is the saturation density of nuclear matter), we find that neutron stars lighter than $1.5M_{\odot}$ cannot reach a deconfined phase. Depending on the EOS, neutron stars of more than $1.5M_{\odot}$ can enter a quark phase only if they are spinning faster than about 3 milliseconds as observed now, whereas larger spin periods imply that they are either already quark stars or will never become one.

We review advances in the theoretical and computational studies of disk winds, jets and outflows including: the connection between accretion and jets, the launch of jets from magnetized disks, the coupled evolution of jets and disks, the interaction of magnetized young stellar objects with their surrounding disks and the relevance to outflows, and finally, the link between jet formation and gravitational collapse. We also address the predictions that the theory makes about jet kinematics, collimation, and rotation, that have recently been confirmed by high spatial and spectral resolution observations. Disk winds have a universal character that may account for jets and outflows during the formation of massive stars as well as brown dwarfs.

We have obtained maps of the large scale outflow associated with the UCHII regionG5.89-0.39 in CO and 13CO (J=3-2), SiO (J=8-7,J=5-4), SO2 (13(2,12)-13(1,13)) and H13CO+(J=4-3). From these maps we have been able to determine the mass (3.3 M_sun), momentum (96 M_sun km/s), energy (3.5x10^{46} erg), mechanical luminosity (141 L_sun), and mass loss rate (~1x10^{-3}M_sun/yr) in the large scale outflow. The observationally derived parameters were used toguide 3D magnetohydrodynamic models of the jet entrained outflow. Through the combination of observations and simulations, we suggest that the large scale outflowmay be inclined by approximately 45 deg to the line of sight, and that the jet entraining the observed molecular outflow may have been active for as little as 1000 years, half the kinematic age of the outflow.

We investigate possible formation sites of the cannonballs (in the gamma ray bursts context) by calculating their physical parameters, such as density, magnetic field and temperature close to the origin. Our results favor scenarios where the cannonballs form as instabilities (knots) within magnetized jets from hyperaccreting disks. These instabilities would most likely set in beyond the light cylinder where flow velocity with Lorentz factors as high as 2000 can be achieved. Our findings challenge the cannonball model of gamma ray bursts if these indeed form inside core-collapse supernovae (SNe) as suggested in the literature; unless hyperaccreting disks and the corresponding jets are common occurrences in core-collapse SNe.

The cooling history of a quark star in the color superconductive phase is investigated. Here we specifically focus on the 2-flavour color (2SC) phase where novel process of photon generation via glueball (GLB) decay have been already investigated (Ouyed & Sannino 2001). The picture we present here can in principle be generalized to quark stars entering a superconductive phase where similar photon generation mechanisms are at play. As much as 10^{45}-10^{47} erg of energy is provided by the GLB decay in the 2SC phase. The generated photons slowly diffuse out of the quark star keeping it hot and radiating as a black-body (with possibly a Wien spectrum in gamma-rays) for millions of years. We discuss hot radio-quiet isolated neutron stars in our picture (such as RX J185635-3754 and RX J0720.4-3125) and argue that their nearly blackbody spectra (with a few broad features) and their remarkably tiny hydrogen atmosphere are indications that these might be quark stars in the color superconductive phase where some sort of photon generation mechanism (reminiscent of the GLB decay) has taken place. Fits to observed data of cooling compact stars favor models with superconductive gaps of Delta_2SC = 15-35 MeV and densities rho_2SC=(2.5-3.0)rho_N (rho_N being the nuclear matter saturation density) for quark matter in the 2SC phase. If correct, our model combined with more observations of isolated compact stars could provide vital information to studies of quark matter and its exotic phases.

We study mass-radius relations for compact stars employing an equation of state (EOS) of dense matter based on a Skyrme fluid. The zero-temperature mean-field model is based on mesonic excitations, incorporates the scale breaking of QCD, and accommodates baryons (nucleons) which arise as a solitonic configuration of mesonic fields. Stable configurations are obtained for central densities $\rho_c/\rho_n\leq 5.0$ where $\rho_n=2.575\times 10^{14}$ g/cc is the nuclear saturation density. These {\it Skyrmion Stars} are mostly fluid, with a crust which we describe by the EOS of Baym, Pethick and Sutherland. Their masses and radii are in the range $0.4 \le M/M_{\odot} \le 3.6$ and $13 {\rm km} \le R \le 23 {\rm km}$, respectively. The minimum spin period is computed to be between $0.7 {\rm ms}\le P \le 2.1 {\rm ms}$. They appear to have a mass-radius curve quite different from either neutron or quark stars, and provide a suitable description of the heavier mass neutron stars discovered recently due to the inherently stiff EOS. Within the same model, we compute the dominant neutrino emissivity in neutron-rich $\beta$-equilibrated matter, and determine the cooling behaviour of Skyrmion stars.

In a previous paper (Ouyed et al. 2004) we presented a new model for soft gamma-ray repeaters (SGR), based on the onset of colour superconductivity in quark stars. In this model, the bursts result from the reorganization of the exterior magnetic field following the formation of vortices that confine the internal magnetic field (the \mz effect). Here we extend the model by presenting full 3-dimensional simulations of the evolution of the inclined exterior magnetic field immediately following vortex formation. The simulations capture the violent reconnection events in the entangled surface magnetic field as it evolves into a smooth, more stable, configuration which consists of a dipole field aligned with the star's rotation axis. The total magnetic energy dissipated in this process is found to be of the order of $10^{44} \mathrm{erg}$ and, if it is emitted as synchrotron radiation, peaks typically at $280 \mathrm{keV}$. The intensity decays temporally in a way resembling SGRs and AXPs (anomalous X-ray pulsars), with a tail lasting from a few to a few hundred times the rotation period of the star, depending on the initial inclination between the rotation and dipole axis. One of the obvious consequences of our model's final state (aligned rotator) is the suppression of radio-emission in SGRs and AXPs following their bursting era. We suggest that magnetar-like magnetic field strength alone cannot be responsible for the properties of SGRs and AXPs, while a quark star entering the ``Meissner phase'' is compatible with the observational facts. We compare our model to observations and highlight our predictions.

(Abriged) We present a comprehensive set of axisymmetric, time-dependent simulations of jets from Keplerian disks whose mass loading as a function of disk radius is systematically changed. For a reasonable model for the density structure and injection speed of the underlying accretion disk, mass loading is determined by the radial structure of the disk's magnetic field structure. We vary this structure by using four different magnetic field configurations, ranging from the "potential" configuration (Ouyed&Pudritz 1997), to the increasingly more steeply falling Blandford&Payne (1982) and Pelletier&Pudritz (1992) models, and ending with a quite steeply raked configuration that bears similarities to the Shu X-wind model. We find that the radial distribution of the mass load has a profound effect on both the rotational profile of the underlying jet as well as the degree of collimation of its outflow velocity and magnetic field lines. We show analytically, and confirm by our simulations, that the collimation of a jet depends on its radial current distribution, which in turn is prescribed by the mass load. Models with steeply descending mass loads have strong toroidal fields, and these collimate to cylinders (this includes the Ouyed-Pudritz and Blandford-Payne outflows). On the other hand, the more gradually descending mass load profiles (the PP92 and monopolar distributions) have weaker toroidal fields, and these result in wide-angle outflows with parabolic collimation. We also present detailed structural information about jets such as their radial profiles of jet density, toroidal magnetic field, and poloidal jet speed, as well as an analysis of the bulk energetics of our different simulations.

Recent studies of photon-generation mechanisms in the color-superconducting Color-Flavor Locked (CFL) phase of dense quark matter have found gamma-ray emissivities in excess of ~ 10^{50} erg cm^{-3} s^{-1} for temperatures in the 10-30 MeV range. We suggest that this property can trigger gamma-ray bursts (GRBs) and associated fireballs at the surface of hypothetical hot (newly born) quark stars with an energy release of up to 10^{48}-10^{50} erg within a fraction of a millisecond. If surrounded by an accretion disk following its formation, the star's bursting activity can last from tens of milliseconds to hundreds of seconds releasing up to 10^{52} erg in total energy. We discuss typical features of observed GRBs within our model and explain how quark stars in the CFL phase might constitute natural candidates for corresponding inner engines.

We have carried out 2D and 3D general relativistic magnetohydrodynamic simulations of jets launched self-consistently from accretion disks orbiting Kerr black holes and applied the results to the inner engine of the collapsar model of gamma-ray bursts. The accretion flow launches energetic jets in the axial funnel region of the disk/jet system, as well as a substantial coronal wind. The jets feature knot-like structures of extremely hot, ultra-relativistic gas; the gas in these knots begins at moderate velocities near the inner engine, and is accelerated to ultra-relativistic velocities (Lorentz factors of 50, and higher) by the Lorentz force in the axial funnel. The increase in jet velocity takes place in an acceleration zone extending to at least a few hundred gravitational radii from the inner engine. The overall energetics of the jets are strongly spin-dependent, with high-spin black holes producing the highest energy and mass fluxes. In addition, with high-spin black holes, the ultra-relativistic outflow is cylindrically collimated within a few hundred gravitational radii of the black hole, whereas in the zero- spin case the jet retains a constant opening angle of approximately 16 degrees. The simulations also show that the coronal wind, though considerably slower and colder than the jets, also carries a significant amount of mass and energy. When simulation data is scaled to the physical dimensions of a collapsar the jets operate for a period ranging from 0.1 to 1.4 seconds, until the accretion disk is depleted, delivering $10^{48}$ to $10^{49}$ erg.

We explore acceleration of ions in the Quark Nova (QN) scenario, where a neutron star experiences an explosive phase transition into a quark star (born in the propeller regime). In this picture, two cosmic ray components are isolated: one related to the randomized pulsar wind and the other to the propelled wind, both boosted by the ultra-relativistic Quark Nova shock. The latter component acquires energies $10^{15} {\rm eV}<E<10^{18} {\rm eV}$ while the former, boosted pulsar wind, achieves ultra-high energies $E> 10^{18.6}$ eV. The composition is dominated by ions present in the pulsar wind in the energy range above $10^{18.6}$ eV, while at energies below $10^{18}$ eV the propelled ejecta, consisting of the fall-back neutron star crust material from the explosion, is the dominant one. Added to these two components, the propeller injects relativistic particles with Lorentz factors $\Gamma_{\rm prop.} \sim 1-1000$, later to be accelerated by galactic supernova shocks. The QN model appears to be able to account for the extragalactic cosmic rays above the ankle and to contribute a few percent of the galactic cosmic rays below the ankle. We predict few hundred ultra-high energy cosmic ray events above $10^{19}$ eV for the Pierre Auger detector per distant QN, while some thousands are predicted for the proposed EUSO and OWL detectors.

We explore the role of neutrinos in a Quark Nova explosion. We study production of neutrinos during this event, their propagation and their interactions with the surrounding quark matter and neutron-rich envelope. We address relevant physical issues such as the timescale for the initial core collapse, the total energy emitted in neutrinos and their effect on the low density matter surrounding the core. We find that it is feasible that the neutrino burst can lead to significant mass ejection of the nuclear envelope.

We present a new model for soft gamma-ray repeaters based on a quark star born with temperatures above the critical value (T_c) for the onset of the colour-flavor locked superconductivity. The quark star then quickly cools below T_c, expelling a fraction of the surface magnetic field via the Meissner effect. We show that if a small fraction (\leq 10%) of the surface magnetic field (10^{14} - 10^{15} {\rm G}) is expelled, it quickly decays via magnetic reconnection and heats up the quark star surface to temperatures > 10^9 {\rm K}. Created (e^{+},e^{-}) pairs annihilate into gamma rays emitted in a giant burst (the first burst in our model), with a luminosity of \sim 10^{45} {\rm ergs} {\rm s}^{-1}. Subsequent bursts result from the restructuring of the surface magnetic field following the formation and relaxation of a vortex lattice which confines the internal magnetic field. During this phase, energy is sporadically released as a consequence of magnetic reconnection events in the entangled surface magnetic field as it evolves into a smooth, more stable, configuration. The star eventually enters a quiescent phase in which energy is continuously supplied by vortex annihilation at the surface. As the star spins down, the outermost vortex lines will be pushed to the surface where they annihilate and release their confined magnetic field. We show that the corresponding luminosity is L_v \sim 10^{36} {\rm ergs} {\rm s}^{-1} for a typical soft gamma-ray repeater spinning with a period of 8 {\rm s} and a surface magnetic field not exceeding 10^{15} {\rm G}. Our model can be applied to any situation where a T>T_{\rm c} quark star is generated. We discuss the connection between anomalous X-ray pulsars and soft gamma-ray repeaters in the context of our model.

Thermal emission rates and mean free paths of photons in a color-flavor locked (CFL) phase of quark matter at high densities and moderate temperatures are evaluated. Our calculations are based on a low-energy effective theory for CFL matter describing Goldstone boson excitations and their electromagnetic as well as strong interactions. In-medium coupling strengths of vector mesons are found to be smaller than in vacuum. As a consequence of in-medium modified pion dispersion relations, novel processes such as pi+ pi- -> gamma and gamma -> pi+ pi- become possible. The total photon emissivity is found to be very large, exceeding contributions from thermal e+ e- annihilation above temperatures of about 5 MeV. At the same time, the corresponding mean free paths become very small. Our results imply that the photon flux from the surface of a (hypothetical) CFL star in its early hot stages saturates the black-body limit. Estimates for the early thermal evolution of the star are also presented.

In simulations of (magnetized-)fluid dynamics in physics and astrophysics, the visualization techniques are so frequently applied to analyse data that they have become a fundamental part of the research. Data produced is often a multi-dimensional set with several physical quantities, that are usually complex to manage and analyse. JETGET is a visualization and analysis tool we developed for accessing data stored in Hierarchical Data Format (HDF) and ASCII files. Although JETGET has been optimized to handle data output from jet simulations using the Zeus code from NCSA, it is also capable of analysing other data output from simulations using other codes. JETGET can select variables from the data files, render both two- and three-dimensional graphics and analyse and plot important physical quantities. Graphics can be saved in encapsulated Postscript, JPEG, VRML or saved into an MPEG for later visualization and/or presentations. An example of use of JETGET in analysing a 3-dimensional simulation of jets emanating from accretion disks surrounding a protostar is shown. The strength of JETGET in extracting the physics underlying such phenomena is demonstrated as well as its capabilities in visualizing the 3-dimensional features of the simulated magneto-hydrodynamic jets. The JETGET tool is written in Interactive Data Language (IDL) and uses a graphical user interface to manipulate the data. The tool was developed on a LINUX platform and can be run on any platform that supports IDL. JETGET can be downloaded (including more information about its utilities) from http://www.capca.ucalgary.ca/software.

We discuss the Skyrme model for strong interactions and the concept of a Skyrmion fluid. The corresponding compact objects, namely Skyrmion stars, constructed from the equation of state describing such a fluid are presented and compared to models of neutron stars based on modern equation of states. We suggest plausible Skyrmion star candidates in the 4U 1636-53 and 4U 1820-30 low mass X-ray binary systems where the suggested masses of the accreting compact companion (\sim 2.0M_{\odot}) remain a challenge for neutron star models.

We present a detailed numerical study of the dynamics and evolution of ultrarelativistic magnetohydrodynamic jets in the black hole-disk system under extreme magnetization conditions. We find that Lorentz factors of up to 3000 are achieved and derived a modifiedMichel scaling (Gamma ~ sigma) which allows for a wide variation in the flow Lorentz factor. Pending contamination induced by mass-entrainment, the linear Michel scaling links modulations in the ultrarelativistic wind to variations in mass accretion in the disk for a given magnetization. The jet is asymptotically dominated by the toroidal magnetic field allowing for efficient collimation. We discuss our solutions (jets) in the context of Gamma ray bursts and describe the relevant features such as the high variability in the Lorentz factor and how high collimation angles (~ 0-5 degrees), or cylindrical jets, can be achieved. We isolate a jet instability mechanism we refer to as the "bottle-neck" instability which essentially relies on a high magnetization and a recollimation of the magnetic flux surfaces. The instability occurs at large radii where any dissipation of the magnetic energy into radiation would in principle result in an optically thin emission.

We explore planet formation in the Quark Nova scenario. If a millisecond pulsar explodes as a Quark Nova, a protoplanetary disk can be formed out of the metal rich fall-back material. The propeller mechanism transfers angular momentum from the born quark star to the disk that will go through viscous evolution with later plausible grain condensation and planet formation. As a result, earth-size planets on circular orbits may form within short radii from the central quark star. The planets in the PSR1257+12 system can be explained by our model if the Quark Nova compact remnant is born with a period of $\sim 0.5$ ms following the explosion. We suggest that a good portion of the Quark Nova remnants may harbour planetary systems.

We explore the scenario where the core of a neutron star (having experienced a transition to an up and down quark phase) shrinks into the equilibrated quark object after reaching strange quark matter saturation density (where a composition of up, down and strange quarks is the favored state of matter). The overlaying (envelope) material free-falls following the core contraction releasing upto 10^{53} {\rm ergs} in energy as radiation, partly as a result of the conversion of envelope material to quarks. This phenomena, we named Quark-Nova, leads to a wide variety of ejectae ranging form the Newtonian, "dirty" to the ultra-relativistic fireball. The mass range of the corresponding compact remnant (the quark star) ranges from less than 0.3M_{\odot} up to a solar mass. We discuss the connection between Quark-Novae and Gamma ray bursts and suggest the recently studied GRB011211 event as a plausible Quark-Nova candidate.

We present the extension of previous two-dimensional simulations of the time-dependent evolution of non-relativistic outflows from the surface of Keplerian accretion disks, to three dimensions. The accretion disk itself is taken to provide a set of fixed boundary conditions for the problem. The 3-D results are consistent with the theory of steady, axisymmetric, centrifugally driven disk winds up to the Alfv\'en surface of the outflow. Beyond the Alfv\'en surface however, the jet in 3-D becomes unstable to non-axisymmetric, Kelvin-Helmholtz instabilities. We show that jets maintain their long-term stability through a self-limiting process wherein the average Alfv\'enic Mach number within the jet is maintained to order unity. This is accomplished in at least two ways. First, poloidal magnetic field is concentrated along the central axis of the jet forming a ``backbone'' in which the Alfv\'en speed is sufficiently high to reduce the average jet Alfv\'enic Mach number to unity. Second, the onset of higher order Kelvin-Helmholtz ``flute'' modes (m \ge 2) reduce the efficiency with which the jet material is accelerated, and transfer kinetic energy of the outflow into the stretched, poloidal field lines of the distorted jet. This too has the effect of increasing the Alfv\'en speed, and thus reducing the Alfv\'enic Mach number. The jet is able to survive the onset of the more destructive m=1 mode in this way. Our simulations also show that jets can acquire corkscrew, or wobbling types of geometries in this relatively stable end-state, depending on the nature of the perturbations upon them. Finally, we suggest that jets go into alternating periods of low and high activity as the disappearance of unstable modes in the sub-Alfv\'enic regime enables another cycle of acceleration to super-Alfv\'enic speeds.

A model for Gamma ray bursts inner engine based on quark stars (speculated to exist in nature) is presented. We describe how and why these objects might constitute new candidates for GRB inner engines. At the heart of the model is the onset of exotic phases of quark matter at the surface of such stars, in particular the 2-flavor color superconductivity. A novel feature of such a phase is the generation of particles which are unstable to photon decay providing a natural mechanism for a fireball generation; an approach which is fundamentally different from models where the fireball is generated during collapse or conversion of neutron star to quark star processes. The model is capable of reproducing crucial features of Gamma ray bursts, such as the episodic activity of the engine (multiple and random shell emission) and the two distinct categories of the bursts (two regimes are isolated in the model with \sim 2 s and \sim 81 s burst total duration).

At this conference many interesting talks were presented on the plausible existence of Quark Stars. Other talks dealt with the exotic new phases of quark matter at very high density. Here, I show how combining these two elements might offer a new way of tackling the Gamma Ray Burst puzzle.

In a previous paper, using an equation of state of dense matter representing a fluid of Skyrmions we constructed the corresponding non-rotating compact-star models in hydrostatic equilibrium; these are mostly fluid stars (the Skyrmion fluid) thus naming them {\it Skyrmion Stars}. Here we generalize our previous calculations by constructing equilibrium sequences of rotating Skyrmion stars in general relativity using the computer code {\it RNS} developed by Stergioulas. We calculated their masses and radii to be 0.4 \le M/M_{\odot} \le 3.45, and 13.0 {\rm km}\le R\le 23.0 {\rm km}, respectively (R being the circumferential radius of the star). The period of the maximally rotating Skyrmion stars is calculated to be 0.8 {\rm ms}\le P \le 2.0 {\rm ms}. We find that a gap (the height between the star surface and the inner stable circular orbit) starts to appear for M\sim 2.0M_{\odot}. Specifically, the Skyrmion star mass range with an existing gap is calculated to be 1.8 < M/ M_{\odot} < 3.0 with the corresponding orbital frequency 0.8 {\rm kHz} < \nu_{\rm ISCO} < 1.3 {\rm kHz}. We apply our model to the 4U 1820-30 low mass X-ray binary and suggest a plausible Skyrmion star candidate in the 4U 1636-53 system. We discuss the difficulties encountered by our model in the 4U 0614+09 case with the highest known Quasi-Periodic Oscillation frequency of 1329 Hz. A comparative study of Skyrmion stars and models of neutron stars based on recent/modern equations of state is also presented.

We solve, and provide analytical expressions, for current-free magnetic configurations in the context of initial setups of 3-dimensional simulations of astrophysical jets involving an accretion disk corona in hydrostatic balance around a central object. These configurations which thread through the accretion disk and its corona preserve the initial hydrostatic state. This work sets stage for future 3-dimensional jet simulations (including disk rotation and mass-load) where launching, acceleration and collimation mechanisms can be investigated.

Properties of inertial modes of Jupiter are investigated for an n=1 polytropic description of the planet interior. We use the anelastic approximation to overcome the usual handicap of a severe spherical harmonics truncation. A powerful iterative solver then allows us to compute the frequencies of the most promising low-order modes using as many spherical harmonics as necessary. The induced O(1%) errors of our model are now within observational limits. A plausible seismological model of Jupiter might thus be at hand provided the use of a more realistic description of the planet interior.

Observational advances over the last decade reveal that star formation is associated with the simultaneous presence of gravitationally collapsing gas, bipolar outflow, and an accretion disk. Two theoretical views of star formation suppose that either stellar mass is determined from the outset by gravitational instability, or by the outflow which sweeps away the collapsing envelope of initially singular density distributions. Neither picture appears to explain all of the facts. This contribution examines some of the key issues facing star formation theory.