Using the 2D multi-group, flux-limited diffusion version of the code VULCAN/2D, that also incorporates rotation, we have calculated the collapse, bounce, shock formation, and early post-bounce evolutionary phases of a core-collapse supernova for a variety of initial rotation rates. This is the first series of such multi-group calculations undertaken in supernova theory with fully multi-D tools. We find that though rotation generates pole-to-equator angular anisotropies in the neutrino radiation fields, the magnitude of the asymmetries is not as large as previously estimated. Moreover, we find that the radiation field is always more spherically symmetric than the matter distribution, with its plumes and convective eddies. We present the dependence of the angular anisotropy of the neutrino fields on neutrino species, neutrino energy, and initial rotation rate. Only for our most rapidly rotating model do we start to see qualitatively different hydrodynamics, but for the lower rates consistent with the pre-collapse rotational profiles derived in the literature the anisotropies, though interesting, are modest. This does not mean that rotation does not play a key role in supernova dynamics. The decrease in the effective gravity due to the centripetal effect can be quite important. Rather, it means that when a realistic mapping between initial and final rotational profiles and 2D multi-group radiation-hydrodynamics are incorporated into collapse simulations the anisotropy of the radiation fields may be only a secondary, not a pivotal factor, in the supernova mechanism.

We have developed a time-dependent, multi-energy-group, and multi-angle (S$_n$) Boltzmann transport scheme for radiation hydrodynamics simulations, in one and two spatial dimensions. The implicit transport is coupled to both 1D (spherically-symmetric) and 2D (axially-symmetric) versions of the explicit Newtonian hydrodynamics code VULCAN. The 2D variant, VULCAN/2D, can be operated in general structured or unstructured grids and though the code can address many problems in astrophysics it was constructed specifically to study the core-collapse supernova problem. Furthermore, VULCAN/2D can simulate the radiation/hydrodynamic evolution of differentially rotating bodies. We summarize the equations solved and methods incorporated into the algorithm and present results of a time-dependent 2D test calculation. A more complete description of the algorithm is postponed to another paper. We highlight a 2D test run that follows for 22 milliseconds the immediate post-bounce evolution of a collapsed core. We present the relationship between the anisotropies of the overturning matter field and the distribution of the corresponding flux vectors, as a function of energy group. This is the first 2D multi-group, multi-angle, time-dependent radiation/hydro calculation ever performed in core collapse studies. Though the transport module of the code is not gray and does not use flux limiters (however, there is a flux-limited variant of VULCAN/2D), it still does not include energy redistribution and most velocity-dependent terms.

We investigate the action of the magnetorotational instability (MRI) in the context of iron-core collapse. Exponential growth of the field on the rotation time scale by the MRI will dominate the linear growth process of field line "wrapping" with the same characteristic time. We examine a variety of initial rotation states, with solid body rotation or a gradient in rotational velocity, that correspond to models in the literature. A relatively modest value of the initial rotation, a period of ~ 10 s, will give a very rapidly rotating PNS and hence strong differential rotation with respect to the infalling matter. We assume conservation of angular momentum on spherical shells. Results are discussed for two examples of saturation fields, a fiducial field that corresponds to Alfven velocity = rotational velocity and a field that corresponds to the maximum growing mode of the MRI. Modest initial rotation velocities of the iron core result in sub-Keplerian rotation and a sub-equipartition magnetic field that nevertheless produce substantial MHD luminosity and hoop stresses: saturation fields of order 10^{15} - 10^{16} G develop within 300 msec after bounce with an associated MHD luminosity of about 10^{52} erg/s. Bi-polar flows driven by this MHD power can affect or even cause the explosions associated with core-collapse supernovae.

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.