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.

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.

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

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

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.

The r-mode instability, believed to limit the rotation speed of compact stars, can provide empirical confirmation for the existence of stable deconfined phases of quark matter that are predicted by weak coupling calculations in Quantum Chromodynamics. We construct a model for strange quark stars as heavy as 2 solar masses that are made of superconducting quark matter in the bulk and a thin crystalline quark matter crust. This crystalline quark crust is sufficiently robust to withstand r-mode heating and viscous rubbing for realistic mode amplitudes O(10^{-2}), unlike a crust made of neutron-rich nuclei. The dissipation provided by viscous rubbing at the core-crust boundary is both necessary and sufficient to obtain stable rotation speeds that are consistent with the majority of rapidly spinning pulsars in low mass X-ray binaries. Our analysis implies that while bare strange stars are ruled out by the existence of rapidly spinning pulsars, a strange star with a quark matter crust is a distinct possibility.

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.

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.

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/

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.

The gravitational-wave instability of r-modes in rapidly rotating compact stars is believed to spin them down to angular frequencies of about a tenth of the Kepler frequency soon after their birth in a Supernova. We point out that the r-mode perturbation also impacts the neutrino cooling and viscosity in hot compact stars via processes that restore weak equilibrium. We illustrate this fact with a simple model of spin-down due to gravitational wave emission in compact stars composed entirely of three-flavor degenerate quark matter (a strange quark star). Non-equilibrium neutrino cooling of this oscillating fluid matter is quantified. Our results imply that a consistent treatment of thermal and spin-frequency evolution of a young and hot compact star is a requisite in estimating the persistence of gravitational waves from such a source.

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.

The r-mode instability in rotating compact stars is used to constrain the phase of matter at high density. The color-flavor-locked phase with kaon condensation (CFL-K0) and without (CFL) is considered in the temperature range 10^8K < T <10^{11} K. While the bulk viscosity in either phase is only effective at damping the r-mode at temperatures T > 10^{11} K, the shear viscosity in the CFL-K0 phase is the only effective damping agent all the way down to temperatures T > 10^8 K characteristic of cooling neutron stars. However, it cannot keep the star from becoming unstable to gravitational wave emission for rotation frequencies f ~ 56-11 Hz at T ~ 10^8-10^9 K. Stars composed almost entirely of CFL or CFL-K0 matter are ruled out by observation of rapidly rotating neutron stars, indicating that dissipation at the quark-hadron interface or nuclear crust interface must play a key role in damping the instability.

[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 determine characteristic timescales for the viscous damping of r-mode oscillations in rapidly rotating compact stars that contain quark matter. We present results for the color-flavor-locked (CFL) phase of dense quark matter, in which the up, down and strange quarks are gapped, as well as the normal (ungapped) quark phase. While the ungapped quark phase supports a temperature window between 10^8 K and 5x10^9 K where the r-mode is damped even for rapid rotation, the r-mode in a rapidly rotating pure CFL star is not damped in the temperature range 10^10 K - 10^11 K. Rotating hybrid stars with quark matter cores display an instability window whose width is determined by the amount of quark matter present, and they can have large spin frequencies outside this window. Except at high temperatures T > 10^10 K, the presence of a quark phase allows for larger critical frequencies and smaller spin-periods compared to rotating neutron stars. If low-mass X-ray binaries contain a large amount of ungapped or CFL quark matter, then our estimates of the r-mode instability suggest that there should be a population of rapidly rotating binaries at frequencies greater than 1000 Hz which have not yet been observed.

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.

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.

Under extreme conditions of temperature and/or density, quarks and gluons are expected to undergo a deconfinement phase transition. While this is an ephemeral phenomenon at the ultra-relativistic heavy-ion collider (BNL-RHIC), quark matter may exist naturally in the dense interior of neutron stars. Herein, we present an appraisal of the possible phase structure of dense quark matter inside neutron stars, and the likelihood of its existence given the current status of neutron star observations. We conclude that quark matter inside neutron stars cannot be dismissed as a possibility, although recent observational evidence rules out most soft equations of state.

The existence of deconfined quark matter in the superdense interior of neutron stars is a key question that has drawn considerable attention over the past few decades. Quark matter can comprise an arbitrary fraction of the star, from 0 for a pure neutron star to 1 for a pure quark star, depending on the equation of state of matter at high density. From an astrophysical viewpoint, these two extreme cases are generally expected to manifest different observational signatures. An intermediate fraction implies a hybrid star, where the interior consists of mixed or homogeneous phases of quark and nuclear matter, depending on surface and Coulomb energy costs, as well as other finite size and screening effects. In this brief review article, we discuss what we can deduce about quark matter in neutron stars in light of recent exciting developments in neutron star observations. We state the theoretical ideas underlying the equation of state of dense quark matter, including color superconducting quark matter. We also highlight recent advances stemming from re-examination of an old paradigm for the surface structure of quark stars and discuss possible evolutionary scenarios from neutron stars to quark stars, with emphasis on astrophysical observations.

We investigate the impact of magnetic fields on the electron distribution in the electrosphere of quark stars. For moderately strong magnetic fields $B\sim 10^{13}$G, quantization effects are generally weak due to the large number density of electrons at surface, but can nevertheless affect the spectral features of quark stars. We outline the main observational characteristics of quark stars as determined by their surface emission, and briefly discuss their formation in explosive events termed Quark-Novae, which may be connected to the $r$-process.

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.

If deconfined quark matter exists inside compact stars, the primary cooling mechanism is neutrino radiation via the direct Urca processes d->u+e+antinu_e and u+e->d+nu_e. Below a critical temperature, T_c, quark matter forms a colour superconductor, one possible manifestation of which is a condensate of <ud> quark Cooper pairs in an electric-charge neutralising background of electrons. We compute the neutrino emission rate from such a phase, including charged pair-breaking and recombination effects, and find that on a material temperature domain below T_c the pairing-induced suppression of the neutrino emission rate is not uniformly exponential. If gapless modes are present in the condensed phase, the emissivity at low temperatures is moderately enhanced above that of completely unpaired matter. The importance of charged current pair-breaking processes for neutrino emission both in the fully gapped and partially gapped regimes is emphasised.

A calculation of the current-quark mass dependence of hadron masses can help in using observational data to place constraints on the variation of nature's fundamental parameters. A hadron's sigma-term is a measure of this dependence. The connection between a hadron's sigma-term and the Feynman-Hellmann theorem is illustrated with an explicit calculation for the pion using a rainbow-ladder truncation of the Dyson-Schwinger equations: in the vicinity of the chiral limit sigma_pi = m_pi/2. This truncation also provides a decent estimate of sigma_rho because the two dominant self-energy corrections to the rho-meson's mass largely cancel in their contribution to sigma_rho. The truncation is less accurate for the omega, however, because there is little to compete with an omega->rho+pi self-energy contribution that magnifies the value of sigma_omega by ~25%. A Poincare' covariant Faddeev equation, which describes baryons as composites of confined-quarks and -nonpointlike-diquarks, is solved to obtain the current-quark mass dependence of the masses of the nucleon and Delta, and thereby sigma_N and sigma_Delta. This "quark-core" piece is augmented by the "pion cloud" contribution, which is positive. The analysis yields sigma_N~60MeV and sigma_Delta~50MeV.

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.

We present a calculation of neutrino pair bremsstrahlung due to electron-electron scattering in a relativistic degenerate plasma of electrons. Proper treatment of the in-medium photon propagator, i.e., inclusion of Debye screening of the longitudinal part and Landau damping of the transverse part, leads to a neutrino emissivity which is several orders of magnitude larger than when Debye screening is imposed for the tranverse part. Our results show that this in-medium process can compete with other sources of neutrino radiation and can, in some cases, even be the dominant neutrino emission mechanism. We also discuss the natural extension to quark-quark bremsstrahlung in gapped and ungapped quark matter.

We reexamine the surface composition of strange stars. Strange quark stars are hypothetical compact stars which could exist if strange quark matter was absolutely stable. It is widely accepted that they are characterized by an enormous density gradient ($~10^{26}$ g/cm$^4$) and large electric fields at surface. By investigating the possibility of realizing a heterogeneous crust, comprised of nuggets of strange quark matter embedded in an uniform electron background, we find that the strange star surface has a much reduced density gradient and negligible electric field. We comment on how our findings will impact various proposed observable signatures for strange stars.

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 investigate moduli field dynamics in supergravity/M-theory like set ups where we turn on fluxes along some or all of the extra dimensions. As has been argued in the context of string theory, we observe that the fluxes tend to stabilize the squashing (or shape) modes. Generically we find that at late times the shape is frozen while the radion evolves as a quintessence field. At earlier times we have a phase of radiation domination where both the volume and the shape moduli are slowly evolving. However, depending on the initial conditions and the parameters of the theory, like the value of the fluxes, curvature of the internal manifold and so on, the dynamics of the internal manifold can be richer with interesting cosmological consequences, including inflation.

Observations to date cannot distinguish neutron stars from self-bound bare quark stars on the basis of their gross physical properties such as their masses and radii alone. However, their surface luminosity and spectral characteristics can be significantly different. Unlike a normal neutron star, a bare quark star can emit photons from its surface at super-Eddington luminosities for an extended period of time. We present a calculation of the photon bremsstrahlung rate from the bare quark star's surface, and indicate improvements that are required for a complete characterization of the spectrum. The observation of this distinctive photon spectrum would constitute an unmistakable signature of a strange quark star and shed light on color superconductivity at stellar densities.

The photon emissivity from the bremsstrahlung process ee-> ee\gamma occuring in the electrosphere at the bare surface of a strange quark star is calculated. For surface temperatures T<10^9K, the photon flux exceeds that of e+e- pairs that are produced via the Schwinger mechanism in the presence of a strong electric field that binds electrons to the surface of the quark star. The average energy of photons emitted from the bremsstrahlung process can be 0.5 MeV or more, which is larger than that in e+e- pair annihilation. The observation of this distinctive photon spectrum would constitute an unmistakable signature of a strange quark star and shed light on color superconductivity at stellar densities.

Neutrino emission drives neutron star cooling for the first several hundreds of years after its birth. Given the low energy ($\sim$ keV) nature of this process, one expects very few nonstandard particle physics contributions which could affect this rate. Requiring that any new physics contributions involve light degrees of freedom, one of the likely candidates which can affect the cooling process would be a nonzero magnetic moment for the neutrino. To illustrate, we compute the emission rate for neutrino pair bremsstrahlung in neutron-neutron scattering through photon-neutrino magnetic moment coupling. We also present analogous differential rates for neutrino scattering off nucleons and electrons that determine neutrino opacities in supernovae. Employing current upper bounds from collider experiments on the tau magnetic moment, we find that the neutrino emission rate can exceed the rate through neutral current electroweak interaction by a factor two, signalling the importance of new particle physics input to a standard calculation of relevance to neutron star cooling. However, astrophysical bounds on the neutrino magnetic moment imply smaller effects.

We revisit the emissivity from neutrino pair bremsstrahlung in neutron-neutron scattering, nn -> nn nu nubar, which was calculated from the one-pion exchange potential including correlation effects by Friman and Maxwell. Starting from the free-space low-momentum nucleon-nucleon interaction V_{low k}, we include tensor, spin-orbit and second-order medium-induced non-central contributions to the scattering amplitude in neutron matter. We find that the screening of the nucleon-nucleon interaction reduces the emissivity from neutrino bremsstrahlung for densities below nuclear matter density. We discuss the implications of medium modifications for the cooling of neutron stars via neutrino emission, taking into account recent results for the polarization effects on neutron superfluidity.

We consider the possibility of regarding the maximum temperature $\tm$ of a plasma of inflaton decay products as a Hagedorn temperature connected with the transition from an inflaton-dominated epoch to a radiation-dominated universe. Employing a ``bottom-up'' approach to inflaton decay, we show that the limiting temperature $T_{\rm H}\sim m_{\phi}$, as obtained from a numerical analysis of the evolution equations, can also be understood theoretically in a Hagedorn picture of the phase transition.

We study thermalization of the early Universe when the inflaton can decay into the Standard Model (SM) quarks and gluons, using QCD arguments. We describe the possible formation of the thermal plasma of soft gluons and quarks well before the completion of reheating. Relevant interaction rates of leading order processes and the corresponding thermalization time scale is presented. We discuss hadronization while thermalizing the decay products of the inflaton, when the reheat temperature of the Universe is below the QCD phase transition but above the temperature of the Big Bang nucleosynthesis.

We propose that the origin of ultra-high energy cosmic rays beyond the GZK cutoff and the origin of small cosmological constant can both be explained by vacuum tunneling effects in a theory with degenerate vacua and fermionic doublets. By considering the possibility of tunneling from a particular winding number state, accompanied by violation of some global quantum number of fermions, the small value of the vacuum dark energy and the production of high energy cosmic rays are shown to be related. We predict that the energy of such cosmic rays should be at least $5\times 10^{14}$ GeV.

We calculate neutrino emissivities from the decay and scattering of Goldstone bosons in the color-flavor-locked (CFL) phase of quarks at high baryon density. Interactions in the CFL phase are described by an effective low-energy theory. For temperatures in the tens of keV range, relevant to the long-term cooling of neutron stars, the emissivities involving Goldstone bosons dominate over those involving quarks, because gaps in the CFL phase are $\sim 100$ MeV while the masses of Goldstone modes are on the order of 10 MeV. For the same reason, the specific heat of the CFL phase is also dominated by the Goldstone modes. Notwithstanding this, both the emissivity and the specific heat from the massive modes remain rather small, because of their extremely small number densities. The values of the emissivity and the specific heat imply that the timescale for the cooling of the CFL core in isolation is $\sim 10^{26}$ y, which makes the CFL phase invisible as the exterior layers of normal matter surrounding the core will continue to cool through significantly more rapid processes. If the CFL phase appears during the evolution of a proto-neutron star, neutrino interactions with Goldstone bosons are expected to be significantly more important since temperatures are high enough ($\sim 20-40$ MeV) to admit large number densities of Goldstone modes.

For the low energy Standard Model neutrino-matter interactions, we calculate neutrino pair ($\nu\bar\nu$) emissivites in superfluid quark matter. Just below the critical temperature, Cooper pairs continuously break and recombine, resulting in the emission of $\nu\bar \nu$ pairs with a rate that greatly exceeds the standard quark modified Urca and bremsstrahlung rates. At the same baryon density in baryonic and quark matter, the ratio of baryon to quark $\nu\bar\nu$ emissivites lies in the range 2-5 for the densities of interest in the long-term cooling of solar mass compact stars. We also find that in matter containing hyperons, $\nu\bar\nu$ emission can occur with hyperons of all species.