We present multi-wavelength images observed with SOFIA-FORCAST from ~10 to 40 $\mu$m of 14 protostars, selected as intermediate-mass protostar candidates, as part of the SOFIA Massive (SOMA) Star Formation Survey. We build protostellar spectral energy distributions (SEDs) with the SOFIA observations, together with archival data from Spitzer, Herschel and IRAS. We then fit the SEDs with radiative transfer (RT) models of Zhang & Tan (2018), based on Turbulent Core Accretion theory, to estimate key properties of the protostars. The SEDs generally indicate the validity of these RT models down to intermediate-mass and/or early-stage protostars. The protostars analyzed so far in the SOMA survey span a range of luminosities from ~$10^{2}$ to ~$10^{6} L_{\odot}$, current protostellar masses from ~0.5 to ~30 $M_{\odot}$, and ambient clump mass surface densities, $\Sigma_{\rm cl}$ of 0.1 - 3 g cm$^{-2}$. A wide range of evolutionary states of the individual protostars and of the protocluster environments are also probed. The 19 to 37 $\mu$m spectral index of the sources correlates with outflow cavity opening angle, ratio of this angle to viewing angle, and evolutionary stage. We have also added a sample of ~50 protostellar sources identified from within Infrared Dark Clouds and expected to be at the earliest stages of their evolution. With this global sample, most of the evolutionary stages of high- and intermediate-mass protostars are probed. From the best fitting models of the protostars, there is no evidence of a threshold value of protocluster clump mass surface density being needed to form protostars up to ~25 $M_\odot$. However, to form more massive protostars, there is tentative evidence that $\Sigma_{\rm cl}$ needs to be at least 1 g cm$^{-2}$. We discuss how this is consistent with expectations from core accretion models that include internal feedback from the forming massive star.

We perform a sequence of 3D magnetohydrodynamic (MHD) simulations of the outflow-core interaction for a massive protostar forming via collapse of an initial cloud core of $60~{M_\odot}$. This allows us to characterize the properties of disk wind driven outflows from massive protostars, which can allow testing of different massive star formation theories. It also enables us to assess quantitatively the impact of outflow feedback on protostellar core morphology and overall star formation efficiency. We find that the opening angle of the flow increases with increasing protostellar mass, in agreement with a simple semi-analytic model. Once the protostar reaches $\sim24~{M_\odot}$ the outflow's opening angle is so wide that it has blown away most of the envelope, thereby nearly ending its own accretion. We thus find an overall star formation efficiency of $\sim50\%$, similar to that expected from low-mass protostellar cores. Our simulation results therefore indicate that the MHD disk wind outflow is the dominant feedback mechanism for helping to shape the stellar initial mass function from a given prestellar core mass function.

We present multi-wavelength images observed with SOFIA-FORCAST from $\sim$10 to 40 $\mu$m of seven high luminosity massive protostars, as part of the SOFIA Massive (SOMA) Star Formation Survey. Source morphologies at these wavelengths appear to be influenced by outflow cavities and extinction from dense gas surrounding the protostars. Using these images, we build spectral energy distributions (SEDs) of the protostars, also including archival data from Spitzer, Herschel and other facilities. Radiative transfer (RT) models of Zhang & Tan (2018), based on Turbulent Core Accretion theory, are then fit to the SEDs to estimate key properties of the protostars. Considering the best five models fit to each source, the protostars have masses $m_{*} \sim 12-64 \: M_{\odot}$ accreting at rates of $\dot{m}_{*} \sim 10^{-4}-10^{-3} \: M_{\odot} \: \rm yr^{-1}$ inside cores of initial masses $M_{c} \sim 100-500 \: M_{\odot}$ embedded in clumps with mass surface densities $\Sigma_{\rm cl} \sim 0.1-3 \: \rm g \: cm^{-2}$ and span a luminosity range of $10^{4} -10^{6} \: L_{\odot}$. Compared with the first eight protostars in Paper I, the sources analyzed here are more luminous, and thus likely to be more massive protostars. They are often in a clustered environment or have a companion protostar relatively nearby. From the range of parameter space of the models, we do not see any evidence that $\Sigma_{\rm cl}$ needs to be high to form these massive stars. For most sources the RT models provide reasonable fits to the SEDs, though the cold clump material often influences the long wavelength fitting. However, for sources in very clustered environments, the model SEDs may not be such a good description of the data, indicating potential limitations of the models for these regions.

We present the results of an investigation of the dredge-up and mixing during the merger of two white dwarfs with different chemical compositions by conducting hydrodynamic simulations of binary mergers for three representative mass ratios. In all the simulations, the total mass of the two white dwarfs is $\lesssim1.0~{\rm M_\odot}$. Mergers involving a CO and a He white dwarf have been suggested as a possible formation channel for R Coronae Borealis type stars, and we are interested in testing if such mergers lead to conditions and outcomes in agreement with observations. Even if the conditions during the merger and subsequent nucleosynthesis favor the production of $^{18}{\mathrm O}$, the merger must avoid dredging up large amounts of $^{16}{\mathrm O}$, or else it will be difficult to produce sufficient $^{18}{\mathrm O}$ to explain the oxygen ratio observed to be of order unity. We performed a total of 9 simulations using two different grid-based hydrodynamics codes using fixed and adaptive meshes, and one smooth particle hydrodynamics (SPH) code. We find that in most of the simulations, $>10^{-2}~{\rm M_\odot}$ of $^{16}{\mathrm O}$ is indeed dredged up during the merger. However, in SPH simulations where the accretor is a hybrid He/CO white dwarf with a $\sim 0.1~{\rm M_\odot}$ layer of helium on top, we find that no $^{16}{\mathrm O}$ is being dredged up, while in the $q=0.8$ simulation $<10^{-4}~{\rm M_\odot}$ of $^{16}{\mathrm O}$ has been brought up, making a WD binary consisting of a hybrid CO/He WD and a companion He WD an excellent candidate for the progenitor of RCB stars.

The common envelope binary interaction remains one of the least understood phases in the evolution of compact binaries, including those that result in Type Ia supernovae and in mergers that emit detectable gravitational waves. In this work we continue the detailed and systematic analysis of 3D hydrodynamic simulations of the common envelope interaction aimed at understanding the reliability of the results. Our first set of simulations replicate the 5 simulations of Passy et al. (a 0.88Msun, 90Rsun RGB primary with companions in the range 0.1 to 0.9Msun) using a new AMR gravity solver implemented on our modified version of the hydrodynamic code Enzo. Despite smaller final separations obtained, these more resolved simulations do not alter the nature of the conclusions that are drawn. We also carry out 5 identical simulations but with a 2.0Msun primary RGB star with the same core mass as the Passy et al. simulations, isolating the effect of the envelope binding energy. With a more bound envelope all the companions in-spiral faster and deeper though relatively less gas is unbound. Even at the highest resolution, the final separation attained by simulations with a heavier primary is similar to the size of the smoothed potential even if we account for the loss of some angular momentum by the simulation. As a result we suggest that a ~2.0$Msun RGB primary may possibly end in a merger with companions as massive as 0.6Msun, something that would not be deduced using analytical arguments based on energy conservation.

The discovery of the first electromagnetic counterpart to a gravitational wave signal has generated follow-up observations by over 50 facilities world-wide, ushering in the new era of multi-messenger astronomy. In this paper, we present follow-up observations of the gravitational wave event GW170817 and its electromagnetic counterpart SSS17a/DLT17ck (IAU label AT2017gfo) by 14 Australian telescopes and partner observatories as part of Australian-based and Australian-led research programs. We report early- to late-time multi-wavelength observations, including optical imaging and spectroscopy, mid-infrared imaging, radio imaging, and searches for fast radio bursts. Our optical spectra reveal that the transient source afterglow cooled from approximately 6400K to 2100K over a 7-day period and produced no significant optical emission lines. The spectral profiles, cooling rate, and photometric light curves are consistent with the expected outburst and subsequent processes of a binary neutron star merger. Star formation in the host galaxy probably ceased at least a Gyr ago, although there is evidence for a galaxy merger. Binary pulsars with short (100 Myr) decay times are therefore unlikely progenitors, but pulsars like PSR B1534+12 with its 2.7 Gyr coalescence time could produce such a merger. The displacement (about 2.2 kpc) of the binary star system from the centre of the main galaxy is not unusual for stars in the host galaxy or stars originating in the merging galaxy, and therefore any constraints on the kick velocity imparted to the progenitor are poor.

In this series of papers, we model the formation and evolution of the photoionized region and its observational signatures during massive star formation. Here we focus on the early break out of the photoionized region into the outflow cavity. Using results of 3-D magnetohydrodynamic-outflow simulations and protostellar evolution calculations, we perform post-processing radiative-transfer. The photoionized region first appears at a protostellar mass of 10Msun in our fiducial model, and is confined to within 10-100AU by the dense inner outflow, similar to some observed very small hypercompact HII regions. Since the ionizing luminosity of the massive protostar increases dramatically as Kelvin-Helmholz (KH) contraction proceeds, the photoionized region breaks out to the entire outflow region in <10,000yr. Accordingly, the radio free-free emission brightens significantly in this stage. In our fiducial model, the radio luminosity at 10 GHz changes from 0.1 mJy kpc2 at m=11Msun to 100 mJy kpc2 at 16Msun, while the infrared luminosity increases by less than a factor of two. The radio spectral index also changes in the break-out phase from the optically thick value of 2 to the partially optically thin value of 0.6. Additionally, we demonstrate that short-timescale variation in free-free flux would be induced by an accretion burst. The outflow density is enhanced in the accretion burst phase, which leads to a smaller ionized region and weaker free-free emission. The radio luminosity may decrease by one order of magnitude during such bursts, while the infrared luminosity is much less affected, since internal protostellar luminosity dominates over accretion luminosity after KH contraction starts. Such variability may be observable on timescales as short 10-100 yr, if accretion bursts are driven by disk instabilities.

We present an overview and first results of the Stratospheric Observatory For Infrared Astronomy Massive (SOMA) Star Formation Survey, which is using the FORCAST instrument to image massive protostars from $\sim10$--$40\:\rm{\mu}\rm{m}$. These wavelengths trace thermal emission from warm dust, which in Core Accretion models mainly emerges from the inner regions of protostellar outflow cavities. Dust in dense core envelopes also imprints characteristic extinction patterns at these wavelengths, causing intensity peaks to shift along the outflow axis and profiles to become more symmetric at longer wavelengths. We present observational results for the first eight protostars in the survey, i.e., multiwavelength images, including some ancillary ground-based MIR observations and archival {\it{Spitzer}} and {\it{Herschel}} data. These images generally show extended MIR/FIR emission along directions consistent with those of known outflows and with shorter wavelength peak flux positions displaced from the protostar along the blueshifted, near-facing sides, thus confirming qualitative predictions of Core Accretion models. We then compile spectral energy distributions and use these to derive protostellar properties by fitting theoretical radiative transfer models. Zhang and Tan models, based on the Turbulent Core Model of McKee and Tan, imply the sources have protostellar masses $m_*\sim10$--50$\:M_\odot$ accreting at $\sim10^{-4}$--$10^{-3}\:M_\odot\:{\rm{yr}}^{-1}$ inside cores of initial masses $M_c\sim30$--500$\:M_\odot$ embedded in clumps with mass surface densities $\Sigma_{\rm{cl}}\sim0.1$--3$\:{\rm{g\:cm}^{-2}}$. Fitting Robitaille et al. models typically leads to slightly higher protostellar masses, but with disk accretion rates $\sim100\times$ smaller. We discuss reasons for these differences and overall implications of these first survey results for massive star formation theories.

The common envelope binary interaction occurs when a star transfers mass onto a companion that cannot fully accrete it. The interaction can lead to a merger of the two objects or to a close binary. The common envelope interaction is the gateway of all evolved compact binaries, all stellar mergers and likely many of the stellar transients witnessed to date. Common envelope simulations are needed to understand this interaction and to interpret stars and binaries thought to be the byproduct of this stage. At this time, simulations are unable to reproduce the few observational data available and several ideas have been put forward to address their shortcomings. The need for more definitive simulation validation is pressing, and is already being fulfilled by observations from time-domain surveys. In this article, we present an initial method and its implementation for post-processing grid-based common envelope simulations to produce the light-curve so as to compare simulations with upcoming observations. Here we implemented a zeroth order method to calculate the light emitted from common envelope hydrodynamic simulations carried out with the 3D hydrodynamic code Enzo used in uni-grid mode. The code implements an approach for the computation of luminosity in both optically thick and optically thin regimes and is tested using the first 135 days of the common envelope simulation of Passy et al. (2012), where a 0.8 solar masses red giant branch star interacts with a 0.6 solar masses companion. This code is used to highlight two large obstacles that need to be overcome before realistic light curves can be calculated. We explain the nature of these problems and the attempted solutions and approximations in full detail to enable the next step to be identified and implemented. We also discuss our simulation in relation to recent data of transients identified as common envelope interactions.

Currently there is great interest in the outcomes and astrophysical implications of mergers of double degenerate binaries. In a commonly adopted approximation, the components of such binaries are represented by polytropes with an index n=3/2. We present detailed comparisons of stellar mass-transfer and merger simulations of polytropic binaries that have been carried out using two very different numerical algorithms --- a finite-volume "grid" code and a smoothed-particle hydrodynamics (SPH) code. We find that there is agreement in both the ultimate outcomes of the evolutions and the intermediate stages if the initial conditions for each code are chosen to match as closely as possible. We find that even with closely matching initial setups, the time it takes to reach a concordant evolution differs between the two codes because the initial depth of contact cannot be matched exactly. There is a general tendency for SPH to yield higher mass transfer rates and faster evolution to the final outcome. We also present comparisons of simulations calculated from two different energy equations: in one series we assume a polytropic equation of state and in the other series an ideal gas equation of state. In the latter series of simulations an atmosphere forms around the accretor which can exchange angular momentum and cause a more rapid loss of orbital angular momentum. In the simulations presented here, the effect of the ideal equation of state is to de-stabilize the binary in both SPH and grid simulations, but the effect is more pronounced in the grid code.

We present hydrodynamic simulations of the common envelope binary interaction between a giant star and a compact companion carried out with the adaptive mesh refinement code ENZO and the smooth particle hydrodynamics code PHANTOM. These simulations mimic the parameters of one of the simulations by Passy et al., but assess the impact of a larger, more realistic initial orbital separation on the simulation outcome. We conclude that for both codes the post-common envelope separation is somewhat larger and the amount of unbound mass slightly greater when the initial separation is wide enough that the giant does not yet overflow or just overflows its Roche lobe. PHANTOM has been adapted to the common envelope problem here for the first time and a full comparison with ENZO is presented, including an investigation of convergence as well as energy and angular momentum conservation. We also set our simulations in the context of past simulations. This comparison reveals that it is the expansion of the giant before rapid in-spiral and not spinning up of the star that causes a larger final separation. We also suggest that the large range in unbound mass for different simulations is difficult to explain and may have something to do with simulations that are not fully converged.

The common envelope interaction is thought to be the gateway to all evolved compact binaries and mergers. Hydrodynamic simulations of the common envelope interaction between giant stars and their companions are restricted to the dynamical, fast, in-spiral phase. They find that the giant envelope is lifted during this phase, but remains mostly bound to the system. At the same time, the orbital separation is greatly reduced, but in most simulations it levels off? at values larger than measured from observations. We conjectured that during the post-in-spiral phase the bound envelope gas will return to the system. Using hydrodynamic simulations, we generate initial conditions for our simulation that result in a fall-back disk with total mass and angular momentum in line with quantities from the simulations of Passy et al. We find that the simulated fall-back event reduces the orbital separation efficiently, but fails to unbind the gas before the separation levels off once again. We also find that more massive fall-back disks reduce the orbital separation more efficiently, but the efficiency of unbinding remains invariably very low. From these results we deduce that unless a further energy source contributes to unbinding the envelope (such as was recently tested by Nandez et al.), all common envelope interactions would result in mergers. On the other hand, additional energy sources are unlikely to help, on their own, to reduce the orbital separation. We conclude by discussing our dynamical fall-back event in the context of a thermally-regulated post-common envelope phase.

We present the results of hydrodynamic simulations of the interaction between a 10 Jupiter mass planet and a red or asymptotic giant branch stars, both with a zero-age main sequence mass of 3.5 $M_\odot$. Dynamic in-spiral timescales are of the order of few years and a few decades for the red and asymptotic giant branch stars, respectively. The planets will eventually be destroyed at a separation from the core of the giants smaller than the resolution of our simulations, either through evaporation or tidal disruption. As the planets in-spiral, the giant stars' envelopes are somewhat puffed up. Based on relatively long timescales and even considering the fact that further in-spiral should take place before the planets are destroyed, we predict that the merger would be difficult to observe, with only a relatively small, slow brightening. Very little mass is unbound in the process. These conclusions may change if the planet's orbit enhances the star's main pulsation modes. Based on the angular momentum transfer, we also suspect that this star-planet interaction may be unable to lead to large scale outflows via the rotation-mediated dynamo effect of Nordhaus and Blackman. Detectable pollution from the destroyed planets would only result for the lightest, lowest metallicity stars. We furthermore find that in both simulations the planets move through the outer stellar envelopes at Mach-3 to Mach-5, reaching Mach-1 towards the end of the simulations. The gravitational drag force decreases and the in-spiral slows down at the sonic transition, as predicted analytically.

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.

The Rotten Egg Nebula has at its core a binary composed of a Mira star and an A-type companion at a separation >10 au. It has been hypothesized to have formed by strong binary interactions between the Mira and a companion in an eccentric orbit during periastron passage ~800 years ago. We have performed hydrodynamic simulations of an asymptotic giant branch star interacting with companions with a range of masses in orbits with a range of initial eccentricities and periastron separations. For reasonable values of the eccentricity, we find that Roche lobe overflow can take place only if the periods are <<100 years. Moreover, mass transfer causes the system to enter a common envelope phase within several orbits. Since the central star of the Rotten Egg nebula is an AGB star, we conclude that such a common envelope phase must have lead to a merger, so the observed companion must have been a tertiary companion of a binary that merged at the time of nebula ejection. Based on the mass and timescale of the simulated disc formed around the companion before the common envelope phase, we analytically estimate the properties of jets that could be launched. Allowing for super-Eddington accretion rates, we find that jets similar to those observed are plausible, provided that the putative lost companion was relatively massive.

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.

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.

The R Coronae Borealis (RCB) stars are hydrogen-deficient, variable stars that are most likely the result of He-CO WD mergers. They display extremely low oxygen isotopic ratios, 16O/18O ~ 1 - 10, 12C/13C>=100, and enhancements up to 2.6dex in F and in s-process elements from Zn to La, compared to solar. These abundances provide stringent constraints on the physical processes during and after the double-degenerate merger. As shown before O-isotopic ratios observed in RCB stars cannot result from the dynamic double-degenerate merger phase, and we investigate now the role of the long-term 1D spherical post-merger evolution and nucleosynthesis based on realistic hydrodynamic merger progenitor models. We adopt a model for extra envelope mixing to represent processes driven by rotation originating in the dynamical merger. Comprehensive nucleosynthesis post-processing simulations for these stellar evolution models reproduce, for the first time, the full range of the observed abundances for almost all the elements measured in RCB stars: 16O/18O ratios between 9 and 15, C-isotopic ratios above 100, and ~1.4 - 2.35dex F enhancements, along with enrichments in s-process elements. The nucleosynthesis processes in our models constrain the length and temperature in the dynamic merger shell-of-fire feature as well as the envelope mixing in the post-merger phase. s-process elements originate either in the shell-of-fire merger feature or during the post-merger evolution, but the contribution from the AGB progenitors is negligible. The post-merger envelope mixing must eventually cease ~ 10^6yr after the dynamic merger phase, before the star enters the RCB phase.

We present 30 and 40 micron imaging of the massive protostar G35.20-0.74 with SOFIA-FORCAST. The high surface density of the natal core around the protostar leads to high extinction, even at these relatively long wavelengths, causing the observed flux to be dominated by that emerging from the near-facing outflow cavity. However, emission from the far-facing cavity is still clearly detected. We combine these results with fluxes from the near-infrared to mm to construct a spectral energy distribution (SED). For isotropic emission the bolometric luminosity would be 3.3x10^4 Lsun. We perform radiative transfer modeling of a protostar forming by ordered, symmetric collapse from a massive core bounded by a clump with high mass surface density, Sigma_cl. To fit the SED requires protostellar masses ~20-34 Msun depending on the outflow cavity opening angle (35 - 50 degrees), and Sigma_cl ~ 0.4-1 g cm-2. After accounting for the foreground extinction and the flashlight effect, the true bolometric luminosity is ~ (0.7-2.2)x10^5 Lsun. One of these models also has excellent agreement with the observed intensity profiles along the outflow axis at 10, 18, 31 and 37 microns. Overall our results support a model of massive star formation involving the relatively ordered, symmetric collapse of a massive, dense core and the launching bipolar outflows that clear low density cavities. Thus a unified model may apply for the formation of both low and high mass stars.

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.

A leading formation scenario for R Coronae Borealis (RCB) stars invokes the merger of degenerate He and CO white dwarfs (WD) in a binary. The observed ratio of 16O/18O for RCB stars is in the range of 0.3-20 much smaller than the solar value of ~500. In this paper, we investigate whether such a low ratio can be obtained in simulations of the merger of a CO and a He white dwarf. We present the results of five 3-dimensional hydrodynamic simulations of the merger of a double white dwarf system where the total mass is 0.9 Mdot and the initial mass ratio (q) varies between 0.5 and 0.99. We identify in simulations with $q\lesssim0.7$ a feature around the merged stars where the temperatures and densities are suitable for forming 18O. However, more 16O is being dredged-up from the C- and O-rich accretor during the merger than the amount of 18O that is produced. Therefore, on a dynamical time scale over which our hydrodynamics simulation runs, a 16O/18O ratio of ~2000 in the "best" case is found. If the conditions found in the hydrodynamic simulations persist for 10^6 seconds the oxygen ratio drops to 16 in one case studied, while in a hundred years it drops to ~4 in another case studied, consistent with the observed values in RCB stars. Therefore, the merger of two white dwarfs remains a strong candidate for the formation of these enigmatic stars.

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.

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 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 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.

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.

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.

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 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.

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.

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 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.

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.

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.