Fusion of very neutron rich nuclei may be important to determine the composition and heating of the crust of accreting neutron stars. Fusion cross sections are calculated using time-dependent Hartree-Fock theory coupled with density-constrained Hartree-Fock calculations to deduce an effective potential. Systems studied include 16O+16O, 16O+24O, 24O+24O, 12C+16O, and 12C+24O. We find remarkable agreement with experimental cross sections for the fusion of stable nuclei. Our simulations use the SLy4 Skyrme force that has been previously fit to the properties of stable nuclei, and no parameters have been fit to fusion data. We compare our results to the simple S\~{a}o Paulo static barrier penetration model. For the asymmetric systems 12C+24O or 16O+24O we predict an order of magnitude larger cross section than those predicted by the S\~{a}o Paulo model. This is likely due to the transfer of neutrons from the very neutron rich nucleus to the stable nucleus and dynamical rearrangements of the nuclear densities during the collision process. These effects are not included in potential models. This enhancement of fusion cross sections, for very neutron rich nuclei, can be tested in the laboratory with radioactive beams.

We couple two-dimensional hydrodynamics to realistic one-dimensional multigroup flux-limited diffusion neutrino transport to investigate protoneutron star convection in core collapse supernovae, and more specifically, the interplay between its development and neutrino transport, for both 15 and 25 solar mass models. In the presence of neutrino transport, protoneutron star convection velocities are too small relative to bulk inflow velocities to result in any significant convective transport of entropy and leptons. A simple analytical model supports our numerical results, indicating that the inclusion of neutrino transport reduces the entropy-driven (lepton-driven) convection growth rates and asymptotic velocities by a factor of 3 (50) at the neutrinosphere and a factor 250 (1000) at a density of 10^{12} g/cm^{3}, for both our 15 and 25 solar mass models. Moreover, when transport is included, the initial postbounce entropy gradient is smoothed out by neutrino diffusion, whereas the initial lepton gradient is maintained by electron capture and neutrino escape near the neutrinosphere. Despite the maintenance of the lepton gradient, protoneutron star convection does not develop over the 100 ms duration typical of all our simulations, except in the instance where ``low-test'' initial conditions are used, which are generated by core collapse and bounce simulations that neglect neutrino-electron scattering and ion-ion screening corrections to neutrino- nucleus elastic scattering.

We investigate neutrino-driven convection in core collapse supernovae and its ramifications for the explosion mechanism, for a 15 solar mass model. Our two-dimensional simulation begins at 12 ms after bounce and proceeds for 500 ms. We couple two-dimensional PPM hydrodynamics to precalculated one-dimensional MGFLD neutrino transport. (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 (a) is multigroup and (b) simulates with greater 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 212 ms after bounce, we see large-scale, neutrino-driven convection beneath the shock, characterized by higher-entropy, expanding upflows and denser, lower-entropy, finger-like downflows. The radial convection velocities at this time become supersonic just below the shock, reaching magnitudes in excess of 10^{9} cm/sec. Eventually, however, the shock recedes to smaller radii, and at 500 ms after bounce there is no evidence in our simulation of an explosion. 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 investigate neutrino-driven convection in core collapse supernovae and its ramifications for the explosion mechanism. We begin with an ``optimistic'' 15 solar mass precollapse model, which is representative of the class of stars with compact iron cores. This model is evolved through core collapse and bounce in one dimension using multigroup (neutrino-energy--dependent) flux-limited diffusion (MGFLD) neutrino transport and Lagrangian hydrodynamics, providing realistic initial conditions for the postbounce convection and evolution. Our two-dimensional simulation begins at 106 ms after bounce at a time when there is a well-developed gain region, and proceeds for 400 ms. We couple two-dimensional (PPM) hydrodynamics to one-dimensional MGFLD neutrino transport. At 225 ms after bounce we see large-scale convection behind the shock, characterized by high-entropy, mushroom-like, expanding upflows and dense, low-entropy, finger-like downflows. The upflows reach the shock and distort it from sphericity. The radial convection velocities become supersonic just below the shock, reaching magnitudes in excess of 10^9 cm/sec. Eventually, however, the shock recedes to smaller radii, and at about 500 ms after bounce there is no evidence in our simulation of an explosion or of a developing explosion. Failure in our ``optimistic'' 15 solar mass 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.

We couple two-dimensional hydrodynamics to detailed one-dimensional multigroup flux-limited diffusion neutrino transport to investigate prompt convection in core collapse supernovae. Our initial conditions, time-dependent boundary conditions, and neutrino distributions for computing neutrino heating, cooling, and deleptonization rates are obtained from one-dimensional simulations that implement multigroup flux-limited diffusion neutrino transport and one-dimensional hydrodynamics. The development and evolution of prompt convection and its ramifications for the shock dynamics are investigated for both 15 and 25 solar mass models, representative of the two classes of stars with compact and extended iron cores, respectively. In the absence of neutrino transport, prompt convection develops and dissipates on a time scale $\sim$15 ms for both models. Prompt convection seeds convection behind the shock, which causes distortions in the shock's sphericity, but on the average, the shock radius is not boosted significantly. In the presence of neutrino transport, prompt convection velocities are too small relative to bulk inflow velocities to result in any significant convective transport of entropy and leptons. A simple analytical model supports our numerical results, indicating that the inclusion of transport reduces the convection growth rates and asymptotic velocities by factors of 4--250.