We examine the conditions for the revival of the stalled accretion shock in core-collapse supernovae, in the context of the neutrino heating mechanism. We combine one dimensional simulations of the shock revival process with a derivation of a quasi-stationary approximation, which is both accurate and efficient in predicting the flow. In particular, this approach is used to explore how the evolution of the system depends on the shock radius, $R_S$, and velocity, $V_S$ (in addition to other global properties of the system). We do so through a phase space analysis of the shock acceleration, $a_S$, in the $R_S-V_S$ plane, shown to provide quantitative insights into the initiation of runaway expansion and its nature. In the particular case of an initially stationary ($V_S=0,\;a_S=0$) profile, the prospects for an explosion can be reasonably assessed by the initial signs of the partial derivatives of the shock acceleration, in analogy to a linear damped/anti-damped oscillator. If $\partial a_S/\partial R_S<0$ and $\partial a_S/\partial V_S>0$, runaway expansion will likely occur after several oscillations, while if $\partial a_S/\partial R_S>0$, runaway expansion will commence in a non-oscillatory fashion. These two modes of runaway correspond to low and high mass accretion rates, respectively. We also use the quasi-stationary approximation to assess the advection-to-heating timescale ratio in the gain region, often used as an explosion proxy. Indeed, this ratio does tend to $\sim1$ in conjunction with runaway conditions, but neither this unit value nor the specific choice of the gain region as a point of reference appear to be distinct conditions in this regard.

Galaxy evolution depends strongly on the physics of the interstellar medium (ISM). Motivated by the need to incorporate the properties of the ISM in cosmological simulations we construct a simple method to include the contribution of non-thermal components in the calculation of pressure of interstellar gas. In our method we treat three non-thermal components - turbulence, magnetic fields and cosmic rays - and effectively parametrize their amplitude. We assume that the three components settle into a quasi-steady-state that is governed by the star formation rate, and calibrate their magnitude and density dependence by the observed Radio-FIR correlation, relating synchrotron radiation to star formation rates of galaxies. We implement our model in single cell numerical simulation of a parcel of gas with constant pressure boundary conditions and demonstrate its effect and potential. Then, the non-thermal pressure model is incorporated into RAMSES and hydrodynamic simulations of isolated galaxies with and without the non-thermal pressure model are presented and studied. Specifically, we demonstrate that the inclusion of realistic non-thermal pressure reduces the star formation rate by an order of magnitude and increases the gas depletion time by as much. We conclude that the non-thermal pressure can prolong the star formation epoch and achieve consistency with observations without invoking artificially strong stellar feedback.

Recent observations suggest that some supernovae may be the result of an explosion into an optically thick circumstellar material, the product of pre-explosion mass-loss (wind) by the progenitor star. This scenario has been studied previously both analytically and numerically. However, many previous studies base their analysis on the diffusion approximation for radiation transfer, which is inappropriate in the optically thin outer layers of the wind. Here we study the deviations from diffusion, and calculate light curves more accurately using a Monte Carlo approach to photon transfer. We distinguish between "compact" winds, for which the diffusion approximation is appropriate, and "extended" winds, which require a more delicate treatment of the radiation. We show that this effect is more significant than that of the light travel time difference to a distant observer, which has a secondary influence on the light curves of extended-wind systems. We also comment on the applicability of the widely used flux-limited diffusion approximation in this context: we find that it generally does not reproduce the Monte Carlo results. The flux-limited diffusion approximation leads to results which are not only quantitatively, but also qualitatively wrong, in the extended-wind regime.

When a binary star system is tidally disrupted by a supermassive black hole at a galactic nucleus, one star is ejected at a high speed while the other remains in a tightly bound orbit around the black hole. The cluster of tightly bound stars builds over time, eventually creating a steady state in which the rate of collisions between these stars is similar to the rate of capturing new stars. A large fraction of the collisions occur near the periapsis of the orbits around the black hole, where the kinetic energies are sufficient to generate an explosive disruption of the two stars involved. The typical flare brightens for several days, with a peak luminosity that is comparable to the lower-luminosity end of known supernovae. The explosion lightcurve is followed by a longer flare due to accretion of ejected matter onto the black hole. Dedicated searches in the near universe could observe several such "collisional-supernovae" per year.

The explosion of a core-collapse supernova can be approximated by the breakdown of steady-state solutions for accretion onto a proto-neutron star (PNS). We analytically show that as the neutrino luminosity exceeds a critical value L_c, the neutrinosphere pressure exceeds the hydrostatic limit even for an optimal shock radius R. This yields L_c \propto M^2 T^2 (with logarithmic corrections) and R \propto M/T, in agreement with numerical results, where M, T are the PNS mass, neutrino temperature. The near-critical flow can be approximated as a ballistic shell on top of an isothermal layer.

Observations from the last decade have indicated the existence of a general class of superluminous supernovae (SLSNe), in which the peak luminosity exceeds 10^{44} erg/s. Here we focus on a subclass of these events, where the light curve is also tens of days wide, so the total radiated energy is order 10^{51} erg. If the origin of these SLSNe is a core-collapse-driven explosion of a massive star, then the mechanism which converts the explosion energy into radiation must be very efficient (much more than in typical core collapse SNe, where this efficiency is of order one percent). We examine the scenario where the radiated luminosity is due to efficient conversion of kinetic energy of the ejected stellar envelope into radiation by interaction with an optically thick, pre-existing circumstellar material (CSM), presumably the product of a steady wind from the progenitor. We base the analysis on a simple, numerically solved, hydrodynamic diffusion model, which allows us to identify the qualitative behavior of the observable light curves, and to relate them to the parameters of the wind. We specifically show that a wide and superluminous supernova requires the mass of the relevant wind material to be comparable to that of ejected material from the exploding progenitor. We find the wind parameters which explain the peak luminosity and width of the bolometric light curves of three particular SLSNe, namely, SN 2005ap, SN 2006gy, and SN 2010gx, and show that they are best fitted with a wind that extends to a radius of order 10^{15} cm. These results serve as an additional indication that at least some SLSNe are powered by interaction of the ejected material with a steady wind of similar mass.

The breakout of a supernova shock wave through the progenitor star's outer envelope is expected to appear as an X-ray flash. However, if the supernova explodes inside an optically-thick wind, the breakout flash is delayed. We present a simple model for estimating the conditions at shock breakout in a wind based on the general observable quantities in the X-ray flash lightcurve: the total energy E_X, and the diffusion time after the peak, t_diff. We base the derivation on the self-similar solution for the forward-reverse shock structure expected for an ejecta plowing through a pre-existing wind at large distances from the progenitor's surface. We find simple quantitative relations for the shock radius and velocity at breakout. By relating the ejecta density profile to the pre-explosion structure of the progenitor, the model can also be extended to constrain the combination of explosion energy and ejecta mass. For the observed case of XRO08109/SN2008D, our model provides reasonable constraints on the breakout radius, explosion energy, and ejecta mass, and predicts a high shock velocity which naturally accounts for the observed non-thermal spectrum.

Madsen and Larsen proposed (astro-ph/0211597) that strangelets - stable lumps of strange quark matter (SQM) - are candidates for the highest energy cosmic rays. They point out that the expected properties of strangelets make them easier to accelerate and less prone to energy losses than protons or ions. The purpose of this Comment is to point out that if ultra-high-energy cosmic rays (UHECRs) are indeed strangelets, their flux would guarantee that all neutron stars must be strange quark (SQ) stars, while no ``conventional'' (i.e., composed of nucleons and other baryons) neutron stars can coexist. This seems highly unlikely in view of various observed phenomena in neutron stars.

A central supermassive black hole (SMBH) with a mass $10^6-10^9 M_\odot$ appears to be a common feature in nearby galaxies and the likely power source in quasars and active galactic nuclei. We demonstrate that the formation of a central black hole is a natural and inevitable consequence of the gravothermal catastrophe in a self-interacting dark matter (SIDM) halo. Through gravothermal evolution driven by collisional relaxation, an SIDM halo will form a massive inner core whose density and velocity dispersion will increase secularly in time. Eventually, the inner core arrives at a relativistic radial instability and undergoes dynamical collapse to a black hole. The initial mass of the black hole will be $10^{-8}-10^{-6}$ of the total mass of the halo. We show that if at formation the overdensity in the SIDM halo is not too large, SMBHs in the observed mass range can form directly in very massive halos following core collapse. Alternatively, with large overdensities, moderate mass halos undergo core collapse to form central seed black holes of intermediate mass, and these holes can then merge and/or accrete to reach the SMBH range. Forming SMBHs by core collapse in SIDM halos requires no baryons, no prior epoch of star formation and no other mechanism of forming black holes seeds.

We study the evolution of an isolated, spherical halo of self-interacting dark matter (SIDM) in the gravothermal fluid formalism. We show that the thermal relaxation time, $t_r$, of a SIDM halo with a central density and velocity dispersion of a typical dwarf galaxy is significantly shorter than its age. We find a self-similar solution for the evolution of a SIDM halo in the limit where the mean free path between collisions, $\lambda$, is everywhere longer than the gravitational scale height, $H$. Typical halos formed in this long mean free path regime relax to a quasistationary gravothermal density profile characterized by a nearly homogeneous core and a power-law halo where $\rho \propto r^{-2.19}$. We solve the more general time-dependent problem and show that the contracting core evolves to sufficiently high density that $\lambda$ inevitably becomes smaller than $H$ in the innermost region. The core undergoes secular collapse to a singular state (the ``gravothermal catastrophe'') in a time $t_{coll} \approx 290 t_r$, which is longer than the Hubble time for a typical dark matter-dominated galaxy core at the present epoch. Our model calculations are consistent with previous, more detailed, N-body simulations for SIDM, providing a simple physical interpretation of their results and extending them to higher spatial resolution and longer evolution times. At late times, mass loss from the contracting, dense inner core to the ambient halo is significantly moderated, so that the final mass of the inner core may be appreciable when it becomes relativistic and radially unstable to dynamical collapse to a black hole.

A newly formed black hole may be directly identified if late-time accretion of material from the base of the ejected envelope generates a luminosity that is observable in the tail of the supernova light curve. In this work we estimate the rate at which events where the black hole ``emerges'' in the supernova light curve can be detected with present capabilities. Our investigation is based on an analytical model of the accretion luminosity at emergence as a function of progenitor mass, coupled to the inferred rate of observed Type II supernovae in nearby galaxies. We find through a parameter survey that under optimistic assumptions the potential rate of observable events can be as high as several per year. However, supernovae which produce black holes are also likely to be low energy explosions and therefore subluminous, as was the case for the best candidate to date, SN1997D. If black hole-forming supernovae are underdetected owing to lower luminosities, the rate of observing black hole emergence is probably not larger than once every few years. We therefore emphasize the importance of dedicated searches for nearby supernovae as well as faint supernovae projects for improving the prospects of observationally certifying the supernova--black hole connection.

If a light gluino exists, the lightest gluino-containing baryon, the \OSO, is a possible candidate for self-interacting dark matter. In this scenario, the simplest explanation for the observed ratio $\Omega_{dm}/\Omega_b \approx 6-10$ is that $m_{S^0} \sim 900$\MeVcs; this is not at present excluded by particle physics. Such an \OSO could be present in neutron stars, with hyperon formation serving as an intermediate stage. We calculate equilibrium compositions and equation of state for high density matter with the \OSO, and find that for a wide range of parameters the properties of neutron stars with the \OSO are consistent with observations. In particular, the maximum mass of a nonrotating star is $1.7-1.8 M_\odot$, and the presence of the \OSO is helpful in reconciling observed cooling rates with hyperon formation.

We present a new set of spectroscopic and photometric data extending the observations of SN 1997D to over 400 days after the explosion. These observations confirm the peculiar properties of SN 1997D, such as the very low abundance of 56Co (0.002 Msolar) and the low expansion velocity of the ejecta (~1000 km/s). We discuss the implications of these observations for the character of the progenitor and the nature of the remnant showing that a Crab-like pulsar or an accreting neutron star formed in the explosion of a low mass progenitor should already have produced a detectable luminosity at this epoch, in contrast with photometric data. On the other hand, the explosion of a high mass progenitor with the formation of a black hole is consistent with the available observations. The consequences of this conclusion regarding the nature of the explosion and the prospects of directly identifying the black hole are also addressed.

White dwarfs and neutron stars are stellar objects with masses comparable to that of our sun. However, as the endpoint stages of stellar evolution, these objects do not sustain any thermonuclear burning and therefore can no longer support the gravitational load of their own mass by generating thermal pressure. Rather, matter in their interiors is compressed to much higher densities than commonly found in normal stars, and pressure is created by degenerate fermion kinetic energy and particle interactions. As a result, white dwarfs and neutron stars offer unique cosmic laboratories for studying matter at very high densities. In this review we discuss the basic properties of condensed matter at extreme densities and summarize the extent to which these properties can be examined by observations of compact objects.

If a black hole formed in a core-collapse supernova is accreting material from the base of the envelope, the accretion luminosity could be observable in the supernova light curve. Here we continue the study of matter fall back onto a black hole in the wake of a supernova and examine realistic supernovae models which allow for an early emergence of the accretion luminosity. Such cases may provide a direct observational identification of the black hole formed in the aftermath of the explosion. Our approach combines analytic estimates and fully relativistic, radiation-hydrodynamic numerical computations. We employ a numerical hydrodynamical scaling technique to accommodate the diverse range of dynamical time scales in a single simulation. We find that while in typical Type II supernovae heating by radioactive decays dominates the late-time light curve, low-energy explosions of more massive stars should provide an important exception where the accretion luminosity will emerge while it is still relatively large. Our main focus is on the only current candidate for such an observation, the very unusual SN1997D. Due to the low energy of the explosion and the very small ($2\times10^{-3} M_\sun$) inferred mass of Co56 in the ejected envelope, we find that accretion should become the dominant source of its luminosity during the year 2000. The total luminosity at emergence is expected to lie in the range $0.5-3\times10^{36} $ ergs/s, potentially detectable with HST. We also discuss the more favorable case of explosions which eject negligible amounts of radioactive isotopes and find that the black hole is likely to emerge a few tens of days after the explosion, with a luminosity of $\sim 10^{37} $\ergss.

If a black hole formed in a core-collapse supernova is accreting material from the base of the envelope, the accretion luminosity could be observable in the supernova light curve. We present results of a fully relativistic numerical investigation of the fallback of matter onto a black hole in a supernova and examine conditions which would be favorable for detection of the black hole. In general, heating by radioactive decays is likely to prevent practical detection of the black hole, but we show that low energy explosions of more massive stars may provide an important exception. We emphasize the particular case of SN1997D in NGC1536, for which we predict that the presence of a black hole could be inferred observationally within the next year.

We examine the roles the presence of hyperons in the cores of neutron stars may play in determining global properties of these stars. The study is based on estimates that hyperons appear in neutron star matter at about twice the nuclear saturation density, and emphasis is placed on effects that can be attributed to the general multi-species composition of the matter, hence being only weakly dependent on the specific modeling of strong interactions. Our analysis indicates that hyperon formation not only softens the equation of state but also severely constrains its values at high densities. Correspondingly, the valid range for the maximum neutron star mass is limited to about 1.5-1.8 $M_\odot$, which is a much narrower range than available when hyperon formation is ignored. Effects concerning neutron star radii and rotational evolution are suggested, and we demonstrate that the effect of hyperons on the equation of state allows a reconciliation of observed pulsar glitches with a low neutron star maximum mass. We discuss the effects hyperons may have on neutron star cooling rates, including recent results which indicate that hyperons may also couple to a superfluid state in high density matter. We compare nuclear matter to matter with hyperons and show that once hyperons accumulate in neutron star matter they reduce the likelihood of a meson condensate, but increase the susceptibility to baryon deconfinement, which could result in a mixed baryon-quark matter phase.

The implications of hyperon pairing for the thermal evolution of neutron stars containing hyperons are investigated. The outcome of cooling simulations are compared for neutron star models composed only of nucleons and leptons, models including hyperons, and models including pairing of hyperons. We show that lambda and neutron pairing suppresses all possible fast neutrino emission processes in not too massive neutron stars. The inclusion of lambda pairing yields better agreement with X-ray observations of pulsars. Particularly, the surface temperatures deduced from X-ray observations within the hydrogen atmosphere model are more consistent with the thermal history of neutron stars containing hyperons, if the critical temperature for the onset of lambda and nucleon pairing is not too small.

In this work we calculate the $^1S_0$ gap energies of $\Lambda$ hyperons in neutron star matter. The calculation is based on a solution of the BCS gap equation for an effective G-matrix parameterization of the $\Lambda-\Lambda$ interaction with a nuclear matter background, presented recently by Lanskoy and Yamamoto. We find that a gap energy of a few tenths of MeV is expected for $\Lambda$ Fermi momenta up to about 1.3 fm$^{-1}$. Implications for neutron star matter are examined, and suggest the existence of a $\Lambda$ $^1S_0$ superfluid between the threshold baryon density for $\Lambda$ formation and the baryon density where the $\Lambda$ fraction reaches $15-20%$.

Creation of strange quark stars through strong interaction deconfinement is studied based on modern estimates of hyperon formation in neutron stars. The hyperon abundance is shown to be large enough so that if strange quark matter (SQM) is the true ground state of matter, the deconfinement density should be at most 2.5-3 times the nuclear saturation density. If so, deconfinement occurs in neutron stars at birth, and all neutron stars must be strange quark stars. Alternatively, sould observation indicate that some neutron stars have a baryonic interior, SQM is unlikely to be absolutely stable.

An effective equation of state (EoS) which generalizes the Lattimer-Swesty equation for nuclear matter is presented for matter at supernuclear densities including strange baryons. It contains an adjustable baryon potential energy density, based on models of local potentials for the baryon-baryon interactions. The features of the equation rely on the properties of nuclei for the nucleon-nucleon interactions, and mainly on experimental data from hypernuclei for the hyperon-nucleon and hyperon-hyperon interactions. The EoS is used to calculate equilibrium compositions and thermodynamic properties of high density matter with strangeness in two astrophysical contexts: neutron star matter (transparent to neutrinos) and proto-neutron star matter (opaque to neutrinos). The effective EoS reproduces typical properties of high density matter found in theoretical microscopic models. Of these, the main result is that hyperons appear in both types of matter at about twice the nuclear saturation density, and that their appearance significantly softens the EoS. The range of maximal masses of neutron stars found in a comprehensive parameter survey is 1.4-1.7 $M_\odot$. Another typical result is that the maximal mass of a proto-neutron star with strange baryons is higher than that of an evolved neutron star (opposite to the case of nuclear matter), setting the stage for a "delayed collapse" scenario.