General Relativistic Effects in the Core Collapse Supernova Mechanism
We apply our recently developed code for spherically symmetric, fully general relativistic (GR) Lagrangian hydrodynamics and multigroup flux-limited diffusion (MGFLD) neutrino transport to examine the effects of GR on the hydrodynamics and transport during collapse, bounce, and the critical shock reheating phase of core collapse supernovae. GR effects were examined by performing core collapse simulations from several precollapse models in the Newtonian limit, in a hybrid limit consisting of GR hydrodynamics and Newtonian transport, and in the fully GR limit. Comparisons of models computed with GR versus Newtonian hydrodynamics show that collapse to bounce takes slightly less time in the GR limit and that the shock propagates slightly farther out in radius before receding. After a secondary quasi-static rise in the shock radius, the shock radius declines considerably more rapidly in the GR simulations than in the corresponding Newtonian simulations. During the shock reheating phase, core collapse computed with GR hydrodynamics results in a substantially more compact structure from the center out to the stagnated shock, the shock radius being reduced by a factor of 2 after 300 ms for a 25 M model and 600 ms for a 15 M model, times being measured from bounce. The inflow speedmore » of material behind the shock is also increased by about a factor of 2 throughout most of the evolution as a consequence of GR hydrodynamics. Regarding neutrino transport, comparisons show that the luminosity and rms energy of any neutrino flavor during the shock reheating phase increases when switching from Newtonian to GR hydrodynamics. This arises from the close coupling of the hydrodynamics and transport and the effect of GR hydrodynamics to produce more compact core structures, hotter neutrinospheres at smaller radii. On additionally switching from Newtonian to GR transport, gravitational time dilation and redshift effects decrease the luminosities and rms energies of all neutrino flavors. This decrease is less in magnitude than the increase in neutrino luminosities and rms energies that arise when switching from Newtonian to GR hydrodynamics, with the result that a fully GR simulation gives higher neutrino luminosities and harder neutrino spectra than a fully Newtonian simulation of the same precollapse model.« less

title: Stability of Standing Accretion Shocks, with an Eye toward Core-Collapse Supernovae
abstract: We examine the stability of standing, spherical accretion shocks. Accretion shocks arise in core-collapse supernovae (the focus of this paper), star formation, and accreting white dwarfs and neutron stars. We present a simple analytic model and use time-dependent hydrodynamics simulations to show that this solution is stable to radial perturbations. In two dimensions we show that small perturbations to a spherical shock front can lead to rapid growth of turbulence behind the shock, driven by the injection of vorticity from the now nonspherical shock. We discuss the ramifications this instability may have for the supernova mechanism.

title: Cosmic Closure: Relating the Ultimate Fate of Massive Stars and the Ultimate Fate of the Universe
abstract: We give a brief overview of the status of core collapse supernova modeling, particularly as it pertains to predictions of neutrino signatures for the next galactic or near extragalactic supernova. We also consider the implications of neutrino mass for both the supernova mechanism and neutrino signature predictions.

title: An Investigation of Neutrino-driven Convection and the Core Collapse Supernova Mechanism Using Multigroup Neutrino Transport
abstract: We investigate neutrino-driven convection in core collapse supernovae and its ramifications for the explosion mechanism. We begin with a postbounce model that is optimistic in two important respects: (1) we begin with a 15 M☉ precollapse model, which is representative of the class of stars with compact iron cores; (2) we implement Newtonian gravity. Our precollapse model is evolved through core collapse and bounce in one dimension using multigroup (neutrino energy-dependent) flux-limited diffusion (MGFLD) neutrino transport and Newtonian Lagrangian hydrodynamics, providing realistic initial conditions for the postbounce convection and evolution. Our two-dimensional simulation begins at 12 ms after bounce and proceeds for 500 ms. We couple two-dimensional piecewise parabolic method (PPM) hydrodynamics to precalculated one-dimensional MGFLD neutrino transport. (The neutrino distributions used for matter heating and deleptonization in our two-dimensional run are obtained from an accompanying one-dimensional simulation. The accuracy of this approximation is assessed.) For the moment, we sacrifice dimensionality for realism in other aspects of our neutrino transport. MGFLD is an implementation of neutrino transport that simultaneously (1) is multigroup and (2) simulates with sufficient realism the transport of neutrinos in opaque, semitransparent, and transparent regions. Both are crucial to the accurate determination of postshock neutrino heating, which sensitively depends on the luminosities, spectra, and flux factors of the electron neutrinos and antineutrinos emerging from their respective neutrinospheres. By 137 ms after bounce, we see neutrino-driven convection rapidly developing beneath the shock. By 212 ms after bounce, this convection becomes large scale, characterized by higher entropy, expanding upflows and lower entropy, denser, finger-like downflows. The upflows reach the shock and distort it from sphericity. The radial convection velocities at this time become supersonic just below the shock, reaching magnitudes in excess of 109 cm s-1. Eventually, however, the shock recedes to smaller radii, and at ~500 ms after bounce there is no evidence in our simulation of an explosion or of a developing explosion. Our angle-averaged density, entropy, electron fraction, and radial velocity profiles in our two-dimensional model agree well with their counterparts in our accompanying one-dimensional MGFLD run above and below the neutrino-driven convection region. In the convection region, the one-dimensional and angle-averaged profiles differ somewhat because (1) convection tends to flatten the density, entropy, and electron fraction profiles, and (2) the shock radius is boosted somewhat by convection. However, the differences are not significant, indicating that, while vigorous, neutrino-driven convection in our model does not have a significant impact on the overall shock dynamics. The differences between our results and those of other groups are considered. These most likely result from differences in (1) numerical hydrodynamics methods; (2) initial postbounce models, and, most important; (3) neutrino transport approximations. We have compared our neutrino luminosities, rms energies, and inverse flux factors with those from the exploding models of other groups. Above all, we find that the neutrino rms energies computed by our multigroup (MGFLD) transport are significantly lower than the values obtained by Burrows and coworkers, who specified their neutrino spectra by tying the neutrino temperature to the matter temperature at the neutrinosphere and by choosing the neutrino degeneracy parameter arbitrarily, and by Herant and coworkers in their transport scheme, which (1) is gray and (2) patches together optically thick and thin regions. The most dramatic difference between our results and those of Janka and Muller is exhibited by the difference in the net cooling rate below the gain radii: Our rate is 2-3 times greater during the critical 50-100 ms after bounce. We have computed the mass and internal energy in the gain region as a function of time. Up to ~150 ms after bounce, we find that both increase as a result of the increasing gain region volume, as the gain and shock radii diverge. However, at all subsequent times, we find that the mass and internal energy in the gain region decrease with time in accordance with the density falloff in the preshock region and with the flow of matter into the gain region at the shock and out of the gain region at the gain radius. Therefore, we see no evidence in the simulations presented here that neutrino-driven convection leads to mass and energy accumulation in the gain region. We have compared our one- and two-dimensional densities, temperatures, and electron fractions in the region below the electron neutrino and antineutrino gain radii, above which the neutrino luminosities are essentially constant (i.e., the neutrino sources are entirely enclosed), in an effort to assess how spherically symmetric our neutrino sources remain during our two-dimensional evolution, and therefore, in an effort to assess our use of precalculated one-dimensional MGFLD neutrino distributions in calculating the matter heating and deleptonization. We find no difference below the neutrinosphere radii. Between the neutrinosphere and gain radii we find no differences with obvious ramifications for the supernova outcome. We note that the interplay between neutrino transport and convection below the neutrinospheres is a delicate matter and is discussed at greater length in another paper (Mezzacappa and coworkers). However, the results presented therein do support our use of precalculated one-dimensional MGFLD in the present context. Failure in our "optimistic" 15 M☉ Newtonian model leads us to conclude that it is unlikely, at least in our approximation, that neutrino-driven convection will lead to explosions for more massive stars with fatter iron cores or in cases in which general relativity is included.

title: Nuclear cross sections, nuclear structure and stellar nucleosynthesis☆
abstract: The role of nuclear reactions (strong, weak and electromagnetic) and nuclear structure effects are discussed in a number of stellar applications. We address fusion cross sections in stellar evolution, neutrino-induced reactions in type II supernovae, electron captures in type Ia supernovae and fission in the r-process. All of this is discussed in the context of nucleosynthesis products and their role in galactic chemical evolution.

title: ASCERTAINING THE CORE COLLAPSE SUPERNOVA MECHANISM: The State of the Art and the Road Ahead
abstract: ▪ Abstract More than four decades have elapsed since modeling of the core collapse supernova mechanism began in earnest. To date, the mechanism remains elusive, at least in detail, although significant progress has been made in understanding these multiscale, multiphysics events. One-, two-, and three-dimensional simulations of or relevant to core collapse supernovae have shown that (a) neutrino transport, (b) fluid instabilities, (c) rotation, and (d) magnetic fields, together with proper treatments of (e) the sub- and super- nuclear density stellar core equation of state, (f) the neutrino interactions, and (g) gravity are all important. The importance of these ingredients applies to both the explosion mechanism and to phenomena directly associated with the mechanism, such as neutron star kicks, supernova neutrino and gravitational wave emission, and supernova spectropolarimetry. Not surprisingly, current two- and three-dimensional models have yet to include (a)–(d) with sufficient realism. One-dimension...

title: Metastability of the false vacuum in a Higgs-seesaw model of dark energy
abstract: In a recently proposed Higgs-Seesaw model the observed scale of dark energy results from a metastable false vacuum energy associated with mixing of the standard model Higgs particle and a scalar associated with new physics at the GUT or Planck scale. Here we address the issue of how to ensure metastability of this state over cosmological time. We consider new tree-level operators, the presence of a thermal bath of hidden sector particles, and quantum corrections to the effective potential. We find that in the thermal scenario many additional light degrees of freedom are typically required unless coupling constants are somewhat fine-tuned. However quantum corrections arising from as few as one additional light scalar field can provide the requisite support. We also briefly consider implications of late-time vacuum decay for the perdurance of observed structures in the universe in this model.

title: Reionization of Hydrogen and Helium by Early Stars and Quasars
abstract: We compute the reionization histories of hydrogen and helium caused by the ionizing radiation fields produced by stars and quasars. For the quasars we use a model based on halo-merger rates that reproduces all known properties of the quasar luminosity function at high redshifts. The less constrained properties of the ionizing radiation produced by stars are modeled with two free parameters: (i) a transition redshift, ztran, above which the stellar population is dominated by massive, zero-metallicity stars and below which it is dominated by a Scalo mass function; and (ii) the product of the escape fraction of stellar ionizing photons from their host galaxies and the star formation efficiency, fescf*. We constrain the allowed range of these free parameters at high redshifts on the basis of the lack of the H I Gunn-Peterson trough at z 6 and the upper limit on the total intergalactic optical depth for electron scattering, τes < 0.18, from recent cosmic microwave background (CMB) experiments. We find that quasars ionize helium by a redshift z ~ 4, but cannot reionize hydrogen by themselves before z ~ 6. A major fraction of the allowed combinations of fescf* and ztran leads to an early peak in the ionized fraction because of the presence of metal-free stars at high redshifts. This sometimes results in two reionization epochs, namely, an early H II or He III overlap phase followed by recombination and a second overlap phase. Even if early overlap is not achieved, the peak in the visibility function for scattering of the CMB often coincides with the early ionization phase rather than with the actual reionization epoch. Consequently, τes does not correspond directly to the reionization redshift. We generically find values of τes 7%, which should be detectable by the MAP satellite.

title: Final state interactions and new physics in B ---> pi K decays
abstract: Within the standard model, and if one assumes that soft rescattering effects are negligible, the $\mathrm{CP}$ asymmetry ${\mathcal{A}}_{\mathrm{CP}}^{\mathrm{dir}}{(B}^{\ifmmode\pm\else\textpm\fi{}}\ensuremath{\rightarrow}{\ensuremath{\pi}}^{\ifmmode\pm\else\textpm\fi{}}K)$ is predicted to be very small and the ratio ${R=B(B}_{d}\ensuremath{\rightarrow}{\ensuremath{\pi}}^{\ensuremath{\mp}}{K}^{\ifmmode\pm\else\textpm\fi{}}{)/B(B}^{\ifmmode\pm\else\textpm\fi{}}\ensuremath{\rightarrow}{\ensuremath{\pi}}^{\ifmmode\pm\else\textpm\fi{}}K)$ provides a bound on the angle $\ensuremath{\gamma}$ of the unitarity triangle, ${\mathrm{sin}}^{2}\ensuremath{\gamma}l~R.$ We estimate the corrections from soft rescattering effects using an approach based on Regge phenomenology, and find effects of order $10%$ with large uncertainties. In particular, we conclude that ${\mathcal{A}}_{\mathrm{CP}}^{\mathrm{dir}}\ensuremath{\sim}0.2$ and ${\mathrm{sin}}^{2}\ensuremath{\gamma}\ensuremath{\sim}1.2R$ could not be taken unambiguously to signal new physics. Using SU(3) relations, we suggest experimental tests that could constrain the size of the soft rescattering effects thus reducing the related uncertainty. Finally, we study the effect of various models of new physics on ${\mathcal{A}}_{\mathrm{CP}}^{\mathrm{dir}}$ and on $R.$

title: X-ray orbital modulations in intermediate polars
abstract: We present an analysis of 30 archival ASCA and RXTE X-ray observations of 16 intermediate polars to investigate the nature of their orbital modulation. We show that X-ray orbital modulation is widespread amongst these systems, but not ubiquitous as indicated by previous studies that included fewer objects. Only seven of the sixteen systems show a clearly statistically significant modulation depth whose amplitude decreases with increasing X-ray energy. Interpreting this as due to photoelectric absorption in material at the edge of an accretion disc would imply that such modulations are visible for all system inclination angles in excess of 60 degrees. However, it is also apparent that the presence of an X-ray orbital modulation can appear and disappear on a timescale of ∼years or months in an individual system. This may be evidence for the presence of a precessing, tilted accretion disc, as inferred in some low mass X-ray binaries.

title: Broad-band X-ray measurements of the black hole candidate XTE J1908+094
abstract: XTE J1908+094 is an X-ray transient that went into outburst in February 2002. After two months it reached a 2- 250 keV peak fl ux of 1t o 2 10 8 erg cm 2 s 1 . Circumstantial evidence points to an accreting galactic black hole as the origin of the X-radiation: pulsations nor thermonuclear flashes were detected that would identify a neutron star and the spectrum was unusually hard for a neutron star at the outburst onset. We report on BeppoSAX and RXTE All Sky Monitor observations of the broad-band spectrum of XTE J1908+094. The spectrum is consistent with a model consisting of a Comptonization component by a40 keV plasma, a multicolor accretion disk blackbody component with a temperature just below 1 keV and a broad emission line at about 6 keV. The spectrum is heavily absorbed by cold interstellar matter with an equivalent hydrogen column density of 2:5 10 22 cm 2 , which makes it dicult to study the black body component in detail. The black body component exhibits strong evolution about 6 weeks into the outburst. Two weeks later this is followed by a swift decay of the power law component. The broadness of the 6 keV feature may be due to relativistic broadening or Compton scattering of a narrow Fe-K line.