GenASiS Basics provides Fortran 2003 classes furnishing extensible object-oriented utilitarian functionality for large-scale physics simulations on distributed memory supercomputers. This functionality includes physical units and constants; display to the screen or standard output device; message passing; I/O to disk; and runtime parameter management and usage statistics. This revision---Version 3 of Basics---includes a significant name change, some minor additions to functionality, and a major addition to functionality: infrastructure facilitating the offloading of computational kernels to devices such as GPUs.

We use OpenMP to target hardware accelerators (GPUs) on Summit, a newly deployed supercomputer at the Oak Ridge Leadership Computing Facility (OLCF), demonstrating simplified access to GPU devices for users of our astrophysics code GenASiS and useful speedup on a sample fluid dynamics problem. We modify our workhorse class for data storage to include members and methods that significantly streamline the persistent allocation of and association to GPU memory. Users offload computational kernels with OpenMP target directives that are rather similar to constructs already familiar from multi-core parallelization. In this initial example we ask, "With a given number of Summit nodes, how fast can we compute with and without GPUs?", and find total wall time speedups of $\sim 12\mathrm{X}$. We also find reasonable weak scaling up to 8000 GPUs (1334 Summit nodes). We make available the source code from this work at https://github.com/GenASiS/GenASiS_Basics.

A kinetic theory of classical particles serves as a unified basis for developing a geometric $3+1$ spacetime perspective on fluid dynamics capable of embracing both Minkowski and Galilei/Newton spacetimes. Parallel treatment of these cases on as common a footing as possible reveals that the particle four-momentum is better regarded as comprising momentum and \textit{inertia} rather than momentum and energy; and consequently, that the object now known as the stress-energy or energy-momentum tensor is more properly understood as a stress-\textit{inertia} or \textit{inertia}-momentum tensor. In dealing with both fiducial and comoving frames as fluid dynamics requires, tensor decompositions in terms of the four-velocities of observers associated with these frames render use of coordinate-free geometric notation not only fully viable, but conceptually simplifying. A particle number four-vector, three-momentum $(1,1)$ tensor, and kinetic energy four-vector characterize a simple fluid and satisfy balance equations involving spacetime divergences on both Minkowski and Galilei/Newton spacetimes. Reduced to a fully $3+1$ form, these equations yield the familiar conservative formulations of special relativistic and non-relativistic hydrodynamics as partial differential equations in inertial coordinates, and in geometric form will provide a useful conceptual bridge to arbitrary-Lagrange-Euler and general relativistic formulations.

In Newtonian fluid dynamics simulations in which composition has been tracked by a nuclear reaction network, energy generation due to composition changes has generally been handled as a separate source term in the energy equation. A relativistic equation in conservative form for total fluid energy, obtained from the spacetime divergence of the stress-energy tensor, in principle encompasses such energy generation; but it is not explicitly manifest. An alternative relativistic energy equation in conservative form---in which the nuclear energy generation appears explicitly, and that reduces directly to the Newtonian internal + kinetic energy in the appropriate limit---emerges naturally and self-consistently from the difference of the equation for total fluid energy and the equation for baryon number conservation multiplied by the average baryon mass $m$, when $m$ is expressed in terms of contributions from the nuclear species in the fluid, and allowed to be mutable.

The large-scale computer simulation of a system of physical fields governed by partial differential equations requires some means of approximating the mathematical limit of continuity. For example, conservation laws are often treated with a `finite-volume' approach in which space is partitioned into a large number of small `cells,' with fluxes through cell faces providing an intuitive discretization modeled on the mathematical definition of the divergence operator. Here we describe and make available Fortran 2003 classes furnishing extensible object-oriented implementations of simple meshes and the evolution of generic conserved currents thereon, along with individual `unit test' programs and larger example problems demonstrating their use. These classes inaugurate the Mathematics division of our developing astrophysics simulation code GenASiS (General Astrophysical Simulation System), which will be expanded over time to include additional meshing options, mathematical operations, solver types, and solver variations appropriate for many multiphysics applications.

We present an initial report on 160 simulations of a highly simplified model of the post-bounce core-collapse supernova environment in three spatial dimensions (3D). We set different values of a parameter characterizing the impact of nuclear dissociation at the stalled shock in order to regulate the post-shock fluid velocity, thereby determining the relative importance of convection and the stationary accretion shock instability (SASI). While our convection-dominated runs comport with the paradigmatic notion of a `critical neutrino luminosity' for explosion at a given mass accretion rate (albeit with a nontrivial spread in explosion times just above threshold), the outcomes of our SASI-dominated runs are much more stochastic: a sharp threshold critical luminosity is `smeared out' into a rising probability of explosion over a $\sim 20\%$ range of luminosity. We also find that the SASI-dominated models are able to explode with 3 to 4 times less efficient neutrino heating, indicating that progenitor properties, and fluid and neutrino microphysics, conducive to the SASI would make the neutrino-driven explosion mechanism more robust.

Core-collapse supernovae are among the most powerful explosions in the Universe, releasing about $10^{53}~\mbox{erg}$ of energy on timescales of a few tens of seconds. These explosion events are also responsible for the production and dissemination of most of the heavy elements, making life as we know it possible. Yet exactly how they work is still unresolved. One reason for this is the sheer complexity and cost of a self-consistent, multi-physics, and multi-dimensional core-collapse supernova simulation, which is impractical, and often impossible, even on the largest supercomputers we have available today. To advance our understanding we instead must often use simplified models, teasing out the most important ingredients for successful explosions, while helping us to interpret results from higher fidelity multi-physics models. In this paper we investigate the role of instabilities in the core-collapse supernova environment. We present here simulation and visualization results produced by our code GenASiS.

In this paper, the second in a series, we document the algorithms and solvers for compressible nonrelativistic hydrodynamics implemented in GenASiS (General Astrophysical Simulation System)---a new code being developed initially and primarily, though by no means exclusively, for the simulation of core-collapse supernovae. In the Mathematics division of GenASiS we introduce Solvers, which includes finite-volume updates for generic hyperbolic BalanceEquations and ordinary differential equation integration Steps. We also introduce the Physics division of GenASiS; this extends the Manifolds division of Mathematics into physical Spaces, defines StressEnergies, and combines these into Universes. We benchmark the hydrodynamics capabilities of GenASiS against many standard test problems; the results illustrate the basic competence of our implementation, demonstrate the manifest superiority of the HLLC over the HLL Riemann solver in a number of interesting cases, and provide preliminary indications of the code's ability to scale and to function with cell-by-cell fixed-mesh refinement.

GenASiS (General Astrophysical Simulation System) is a new code being developed initially and primarily, though by no means exclusively, for the simulation of core-collapse supernovae on the world's leading capability supercomputers. This paper---the first in a series---demonstrates a centrally refined coordinate patch suitable for gravitational collapse and documents methods for compressible nonrelativistic hydrodynamics. We benchmark the hydrodynamics capabilities of GenASiS against many standard test problems; the results illustrate the basic competence of our implementation, demonstrate the strengths and limitations of the HLLC relative to the HLL Riemann solver in a number of interesting cases, and provide preliminary indications of the code's ability to scale and to function with cell-by-cell fixed-mesh refinement.

We present conservative 3+1 general relativistic variable Eddington tensor radiation transport equations, including greater elaboration of the momentum space divergence (that is, the energy derivative term) than in previous work. These equations are intended for use in simulations involving numerical relativity, particularly in the absence of spherical symmetry. The independent variables are the lab frame coordinate basis spacetime position coordinates and the particle energy measured in the comoving frame. With an eye towards astrophysical applications---such as core-collapse supernovae and compact object mergers---in which the fluid includes nuclei and/or nuclear matter at finite temperature, and in which the transported particles are neutrinos, we pay special attention to the consistency of four-momentum and lepton number exchange between neutrinos and the fluid, showing the term-by-term cancellations that must occur for this consistency to be achieved.

We present a new derivation of the conservative form of the general relativistic Boltzmann equation and specialize it to the 3+1 metric. The resulting transport equation is intended for use in simulations involving numerical relativity, particularly in the absence of spherical symmetry. The independent variables are lab frame coordinate basis spacetime position components and comoving frame curvilinear momentum space coordinates. With an eye towards astrophysical applications---such as core-collapse supernovae and compact object mergers---in which the fluid includes nuclei and/or nuclear matter at finite temperature, and in which the transported particles are neutrinos, we examine the relationship between lepton number and four-momentum exchange between neutrinos and the fluid.

We derive conservative, multidimensional, energy-dependent moment equations for neutrino transport in core-collapse supernovae and related astrophysical systems, with particular attention to the consistency of conservative four-momentum and lepton number transport equations. After taking angular moments of conservative formulations of the general relativistic Boltzmann equation, we specialize to a conformally flat spacetime, which also serves as the basis for four further limits. Two of these---the multidimensional special relativistic case, and a conformally flat formulation of the spherically symmetric general relativistic case---are given in appendices for the sake of comparison with extant literature. The third limit is a weak-field, `pseudo-Newtonian' approach \citep{kim_etal_2009,kim_etal_2012} in which the source of the gravitational potential includes the trace of the stress-energy tensor (rather than just the mass density), and all orders in fluid velocity $v$ are retained. Our primary interest here is in the fourth limit: `$\mathcal{O}(v)$' moment equations for use in conjunction with Newtonian self-gravitating hydrodynamics. We show that the concept of `$\mathcal{O}(v)$' transport requires care when dealing with both conservative four-momentum and conservative lepton number transport, and present two self-consistent options: `$\mathcal{O}(v)$-plus' transport, in which an $\mathcal{O}(v^2)$ energy equation combines with an $\mathcal{O}(v)$ momentum equation to give an $\mathcal{O}(v^2)$ number equation; and `$\mathcal{O}(v)$-minus' transport, in which an $\mathcal{O}(v)$ energy equation combines with an $\mathcal{O}(1)$ momentum equation to give an $\mathcal{O}(v)$ number equation.

The stationary accretion shock instability (SASI) plays a central role in modern simulations of the explosion phase of core-collapse supernovae (CCSNe). It may be key to realizing neutrino powered explosions, and possibly links birth properties of pulsars (e.g., kick, spin, and magnetic field) to supernova dynamics. Using high-resolution magnetohydrodynamic simulations, we study the development of turbulence, and subsequent amplification of magnetic fields in a simplified model of the post-bounce core-collapse supernova environment. Turbulence develops from secondary instabilities induced by the SASI. Our simulations suggest that the development of turbulence plays an important role for the subsequent evolution of the SASI. The turbulence also acts to amplify weak magnetic fields via a small-scale dynamo.

We describe the initial implementation of magnetohydrodynamics (MHD) in our astrophysical simulation code \genasis. Then, we present MHD simulations exploring the capacity of the stationary accretion shock instability (SASI) to generate magnetic fields by adding a weak magnetic field to an initially spherically symmetric fluid configuration that models a stalled shock in the post-bounce supernova environment. Upon perturbation and nonlinear SASI development, shear flows associated with the spiral SASI mode contributes to a widespread and turbulent field amplification mechanism. While the SASI may contribute to neutron star magnetization, these simulations do not show qualitatively new features in the global evolution of the shock as a result of SASI-induced magnetic field amplification.

We extend our investigation of magnetic field evolution in three-dimensional flows driven by the stationary accretion shock instability (SASI) with a suite of higher-resolution idealized models of the post-bounce core-collapse supernova environment. Our magnetohydrodynamic simulations vary in initial magnetic field strength, rotation rate, and grid resolution. Vigorous SASI-driven turbulence inside the shock amplifies magnetic fields exponentially; but while the amplified fields reduce the kinetic energy of small-scale flows, they do not seem to affect the global shock dynamics. The growth rate and final magnitude of the magnetic energy are very sensitive to grid resolution, and both are underestimated by the simulations. Nevertheless our simulations suggest that neutron star magnetic fields exceeding $10^{14}$ G can result from dynamics driven by the SASI, \emph{even for non-rotating progenitors}.

We begin an exploration of the capacity of the stationary accretion shock instability (SASI) to generate magnetic fields by adding a weak, stationary, and radial (but bipolar) magnetic field, and in some cases rotation, to an initially spherically symmetric fluid configuration that models a stalled shock in the post-bounce supernova environment. In axisymmetric simulations we find that cycles of latitudinal flows into and radial flows out of the polar regions amplify the field parallel to the symmetry axis, typically increasing the total magnetic energy by about two orders of magnitude. Nonaxisymmetric calculations result in fundamentally different flows and a larger magnetic energy increase: shearing associated with the SASI spiral mode contributes to a widespread and turbulent field amplification mechanism, boosting the magnetic energy by almost four orders of magnitude (a result which remains very sensitive to the spatial resolution of the numerical simulations). While the SASI may contribute to neutron star magnetization, these simulations do not show qualitatively new features in the global evolution of the shock as a result of SASI-induced magnetic field amplification.

The calculation of neutrino decoupling from nuclear matter requires a transport formalism capable of handling both collisions and flavor mixing. The first steps towards such a formalism are the construction of neutrino and antineutrino "distribution matrices," and a determination of the Liouville equations they satisfy in the noninteracting case. These steps are accomplished through study of a Wigner-transformed "density function," the mean value of paired neutrino quantum field operators.

Formulae for the energies of degenerate non-relativistic and ultra-relativistic Fermi gases play multiple roles in simple arguments related to the collapse of a stellar core to a neutron star. These formulae, deployed in conjunction with the virial theorem and a few other basic physical principles, provide surprisingly good estimates of the temperature, mass, and radius (and therefore also density and entropy) of the core at the onset of collapse; the final radius and composition of the cold compact remnant; and the total energy lost to neutrino emission during collapse.

The classical notion of a single-particle scalar distribution function or phase space density can be generalized to a matrix in order to accommodate superpositions of states of discrete quantum numbers, such as neutrino mass/flavor. Such a `neutrino distribution matrix' is thus an appropriate construct to describe a neutrino gas that may vary in space as well as time and in which flavor mixing competes with collisions. The Liouville equations obeyed by relativistic neutrino distribution matrices, including the spatial derivative and vacuum flavor mixing terms, can be explicitly but elegantly derived in two new ways: from a covariant version of the familiar simple model of flavor mixing, and from the Klein-Gordon equations satisfied by a quantum `density function' (mean value of paired quantum field operators). Associated with the latter derivation is a case study in how the joint position/momentum dependence of a classical gas (albeit with Fermi statistics) emerges from a formalism built on quantum fields.

By adding a weak magnetic field to a spherically symmetric fluid configuration that caricatures a stalled shock in the post-bounce supernova environment, we explore the capacity of the stationary accretion shock instability (SASI) to generate magnetic fields. The SASI develops upon perturbation of the initial condition, and the ensuing flow generates--{\em in the absence of rotation}--dynamically significant magnetic fields ($\sim 10^{15}$ G) on a time scale that is relevant for the explosion mechanism of core-collapse supernovae. We describe our model, present some recent results, and discuss their potential relevance for supernova models.

A nascent neutron star resulting from stellar collapse is a prodigious source of neutrinos of all flavors. While the most basic features of this neutrino emission can be estimated from simple considerations, the detailed simulation of the neutrinos' decoupling from the hot neutron star is not yet computationally tractable in its full glory, being a time-dependent six-dimensional transport problem. Nevertheless, supernova neutrino fluxes are of great interest in connection with the core-collapse supernova explosion mechanism and supernova nucleosynthesis, and as a potential probe of the supernova environment and of some of the neutrino mixing parameters that remain unknown; hence a variety of approximate transport schemes have been used to obtain results with reduced dimensionality. However, none of these approximate schemes have addressed a recent challenge to the conventional wisdom that neutrino flavor mixing cannot impact the explosion mechanism or r-process nucleosynthesis.

While an understanding of supernova explosions will require sophisticated large-scale simulations, it is nevertheless possible to outline the most basic features of the neutrino emission resulting from stellar core collapse with a pedestrian account that, through reliance upon broadly accessible physical ideas, remains simple and largely self-contained.

The computational difficulty of six-dimensional neutrino radiation hydrodynamics has spawned a variety of approximations, provoking a long history of uncertainty in the core-collapse supernova explosion mechanism. Under the auspices of the Terascale Supernova Initiative, we are honoring the physical complexity of supernovae by meeting the computational challenge head-on, undertaking the development of a new adaptive mesh refinement code for self-gravitating, six-dimensional neutrino radiation magnetohydrodynamics. This code--called {\em GenASiS,} for {\em Gen}eral {\em A}strophysical {\em Si}mulation {\em S}ystem--is designed for modularity and extensibility of the physics. Presently in use or under development are capabilities for Newtonian self-gravity, Newtonian and special relativistic magnetohydrodynamics (with `realistic' equation of state), and special relativistic energy- and angle-dependent neutrino transport--including full treatment of the energy and angle dependence of scattering and pair interactions.

Many astrophysical phenomena exhibit relativistic radiative flows. While velocities in excess of $v \sim 0.1c$ can occur in these systems, it has been common practice to approximate radiative transfer to $\cO(v/c)$. In the case of neutrino transport in core-collapse supernovae, this approximation gives rise to an inconsistency between the lepton number transfer and lab frame energy transfer, which have different $\cO(v/c)$ limits. A solution used in spherically symmetric $\cO(v/c)$ simulations has been to retain, for energy accounting purposes, the $\cO(v^2/c^2)$ terms in the lab frame energy transfer equation that arise from the $\cO(v/c)$ neutrino number transport equation. Avoiding the proliferation of such ``extra'' $\cO(v^2/c^2)$ terms in the absence of spherical symmetry motivates a special relativistic formalism, which we exhibit in coordinates sufficiently general to encompass Cartesian, spherical, and cylindrical coordinate systems.

The computational difficulty of six-dimensional neutrino radiation hydrodynamics has spawned a variety of approximations, provoking a long history of uncertainty in the core-collapse supernova explosion mechanism. Under the auspices of the Terascale Supernova Initiative, we are honoring the physical complexity of supernovae by meeting the computational challenge head-on, undertaking the development of a new adaptive mesh refinement code for self-gravitating, six-dimensional neutrino radiation magnetohydrodynamics. This code--called {\em GenASiS}, for {\em Gen}eral {\em A}strophysical {\em Si}mulation {\em S}ystem--is designed for modularity and extensibility of the physics. Presently in use or under development are capabilities for Newtonian self-gravity, Newtonian and special relativistic magnetohydrodynamics (with `realistic' equation of state), and special relativistic energy- and angle-dependent neutrino transport--including full treatment of the energy and angle dependence of scattering and pair interactions.

A new, very fast method for 3D radiative transfer on fully threaded grids with arbitrarily high angular resolution is presented. The method uses completely cell-based discretization, and is ideally suited for problems with diffuse background radiation, often encountered in cosmological and star formation models. We find that for accurate statistical study of intergalactic Ly-alpha absorption lines one needs of order of few hundred angular discretization elements even for models without radiative feedback from star forming galaxies.

A principal `supernova neutrino challenge' is the computational difficulty of six-dimensional neutrino radiation hydrodynamics. The variety of resulting approximations has provoked a long history of uncertainty in the core-collapse supernova explosion mechanism, but recent work highlighting low-mode convection and a newly-recognized instability in spherical accretion shocks may signal (yet another) resolution. As part of its goal of elucidating the explosion mechanism, the Terascale Supernova Initiative is committed to meeting the full complexity of the computational challenge. The understanding of supernova neutrino emission gained in detailed simulations provides a potential basis for learning about two major remaining unknowns in neutrino flavor mixing: the value of the mixing angle $\theta_{13}$, and distinguishing between ``normal'' and ``inverted'' mass hierarchies.

To date, detailed core-collapse supernova simulations have not succeeded in elucidating the explosion mechanism. But whatever the mechanism turns out to be, there will be a neutrino-heated outflow between the hot, newly-born neutron star and the outgoing supernova shock wave whose neutron richness, low density, and high temperature could provide promising conditions for the synthesis of heavy nuclei via rapid neutron capture (the r-process). General relativistic effects improve the prospects for this outflow as an r-process site. Among relativistic effects, the enhanced "gravitational potential" is more important than the gravitational redshift or trajectory bending of neutrinos.

Neutrinos dominate the energetics of core-collapse supernovae, and are believed to play an important role in driving the explosion. The development of a spatially multidimensional neutrino radiative transfer code is a key part of the Terascale Supernova Initiative, a multi-year effort aimed at understanding the supernova explosion mechanism through supercomputer simulations. I will describe the conservative formulations of general relativistic radiative transfer we have developed for this application, and the algorithm we are implementing in our parallel code.

Accurate accounting of particle number and 4-momentum in radiative transfer may be facilitated by the use of transport equations that allow transparent conversion between volume and surface integrals in both spacetime and momentum space. Such conservative formulations of general relativistic radiative transfer in multiple spatial dimensions are presented, and their relevance to core-collapse supernova simulations described.

Core-collapse supernovae are dramatic events with a rich phenomenology, including gravitational radiation. Simulations of these events in multiple spatial dimensions with energy- and angle-dependent neutrino transport are still in their infancy. Core collapse and bounce in a supernova in our galaxy may well be visible by first-generation LIGO, but detailed understanding waits on improvements in modeling both stellar progenitors and the supernova process.

We present an implicit finite difference representation for general relativistic radiation hydrodynamics in spherical symmetry. Our code, Agile-Boltztran, solves the Boltzmann transport equation for the angular and spectral neutrino distribution functions in self-consistent simulations of stellar core collapse and postbounce evolution. It implements a dynamically adaptive grid in comoving coordinates. Most macroscopically interesting physical quantities are defined by expectation values of the distribution function. We optimize the finite differencing of the microscopic transport equation for a consistent evolution of important expectation values. We test our code in simulations launched from progenitor stars with 13 solar masses and 40 solar masses. ~0.5 s after core collapse and bounce, the protoneutron star in the latter case reaches its maximum mass and collapses further to form a black hole. When the hydrostatic gravitational contraction sets in, we find a transient increase in electron flavor neutrino luminosities due to a change in the accretion rate. The muon- and tauon-neutrino luminosities and rms energies, however, continue to rise because previously shock-heated material with a non-degenerate electron gas starts to replace the cool degenerate material at their production site. We demonstrate this by supplementing the concept of neutrinospheres with a more detailed statistical description of the origin of escaping neutrinos. We compare the evolution of the 13 solar mass progenitor star to simulations with the MGFLD approximation, based on a recently developed flux limiter. We find similar results in the postbounce phase and validate this MGFLD approach for the spherically symmetric case with standard input physics.

Experience with core-collapse supernova simulations shows that accurate accounting of total particle number and 4-momentum can be a challenge for computational radiative transfer. This accurate accounting would be facilitated by the use of particle number and 4-momentum transport equations that allow transparent conversion between volume and surface integrals in both configuration and momentum space. Such conservative formulations of general relativistic kinetic theory in multiple spatial dimensions are presented in this paper, and their relevance to core-collapse supernova simulations is described.

Core-collapse supernovae can be used to place limits on dark matter candidate particles, but the strength of these limits depends on the depth of our theoretical understanding of these astrophysical events. To date, limitations on computing power have prevented inclusion of all the physics that would constitute a realistic simulation. The TeraScale Supernova Initiative (TSI) will overcome these obstacles in the next few years, elucidating the explosion mechanism and other phenomena closely associated with the core collapse of massive stars.

We propose an approximation to full relativity that captures the main gravitational effects of dynamical importance in supernovae. The conceptual link between this formalism and the Newtonian limit is such that it could likely be implemented relatively easily in existing multidimensional Newtonian gravitational hydrodynamics codes employing a Poisson solver. As a test of the formalism's utility, we display results for rapidly rotating (and therefore highly deformed) neutron stars.

We propose an approximation to general relativity that captures the main gravitational effects of dynamical importance in supernovae. The conceptual link between this formalism and the Newtonian limit is such that it could likely be implemented in existing multidimensional Newtonian gravitational hydrodynamics codes employing a Poisson solver. As a test of the formalism's utility, we display results for rapidly rotating (and therefore highly deformed) neutron stars.

We study steady state winds from compact objects in the regime where the wind velocity at infinity is ultra-relativistic. This may have relevance to some models of Gamma-Ray-Bursts (GRB's). Particular attention is paid to the case where neutrinos provide the heating. Unless the neutrino luminosity is very large, L>10^{54} erg/s, the only allowed steady state solutions are those where energy deposition is dominated by neutrino-antineutrino annihilation at the sonic point. In this case, the matter temperature near the neutron star surface is low, less than 1MeV for typical neutrino luminosities. This is in contrast to the case for sub-relativistic winds discussed in the context of supernovae where the matter temperature near the neutron star approximates the temperature characterizing the neutrinos. We also investigate the setting of the neutron to proton ratio (n/p) in these winds and find that only for large (>10 MeV) electron neutrino or electron anti-neutrino temperatures is n/p entirely determined by neutrino capture. Otherwise, n/p retains an imprint of conditions in the neutron star.

We study static neutron stars with poloidal magnetic fields and a simple class of electric current distributions consistent with the requirement of stationarity. For this class of electric current distributions, we find that magnetic fields are too large for static configurations to exist when the magnetic force pushes a sufficient amount of mass off-center that the gravitational force points outward near the origin in the equatorial plane. (In our coordinates an outward gravitational force corresponds to $\partial\ln g_{tt}/\partial r>0$, where $t$ and $r$ are respectively time and radial coordinates and $g_{tt}$ is coefficient of $dt^2$ in the line element.) For the equations of state (EOSs) employed in previous work, we obtain configurations of higher mass than had been reported; we also present results with more recent EOSs. For all EOSs studied, we find that the maximum mass among these static configurations with magnetic fields is noticeably larger than the maximum mass attainable by uniform rotation, and that for fixed values of baryon number the maximum mass configurations are all characterized by an off-center density maximum.

We show in detail the general relationship between the Schr\"{o}dinger equation approach to calculating the MSW effect and the quantum field theoretical S-matrix approach. We show the precise form a generic neutrino propagator must have to allow a physically meaningful ``oscillation probability'' to be decoupled from neutrino production fluxes and detection cross-sections, and explicitly list the conditions---not realized in cases of current experimental interest---in which the field theory approach would be useful.

It has recently been suggested that cosmologically significant numbers of black holes could form during a first-order QCD phase transition. Further, it has been asserted that these black holes would have masses corresponding naturally to the inferred mass ($\sim 1 M_{\odot}$) of the MACHOs responsible for the observed gravitational microlensing events. In this model, the underlying spectrum of primordial density perturbations provides the fluctuations that give rise to black holes at the epoch of the QCD transition. We employ a simplified model to estimate the reduction in the critical overdensity of a horizon-sized primordial perturbation required for collapse to a black hole. We find that a first-order QCD transition does indeed produce a sharp peak in the black hole mass spectrum, but that this peak corresponds to the horizon mass at an epoch somewhat earlier than the cosmological transition itself. Assuming a COBE normalized primordial density perturbation spectrum with constant spectral index, for the black holes so produced to be cosmologically significant would require an extremely finely tuned ``blue'' primordial density perturbation spectrum. Specifically, in the context of our simplified model, a spectral index in the range $n=1.37-1.42$ corresponds to the range $\Omega \sim 10^{-5}-10^3$ of the black hole contribution to the present-day density parameter.

Limits can be placed on nonstandard neutrino physics when big bang nucleosynthesis (BBN) calculations employing standard neutrino physics agree with the observationally inferred primordial abundances of deuterium (D), $^3$He, $^4$He, and $^7$Li. These constraints depend most sensitively on the abundances of D and $^4$He. New observational determinations of the primordial D and/or $^4$He abundances could force revisions in BBN constraints on nonstandard neutrino physics.

We discuss general relativistic effects in the steady-state neutrino-driven ``wind'' which may arise from nascent neutron stars. In particular, we generalize previous analytic estimates of the entropy per baryon $S$, the mass outflow rate $\dot M$, and the dynamical expansion time scale $\tau_{dyn}$. We show that $S$ increases and $\tau_{dyn}$ decreases with increasing values of the mass-to-radius ratio describing the supernova core. Both of these trends indicate that a more compact core will lead to a higher number of neutrons per iron peak seed nucleus. Such an enhancement in the neutron-to-seed ratio may be required for successful r-process nucleosynthesis in neutrino-heated supernova ejecta.

We discuss neutrino oscillations in curved spacetime. Our heuristic approach can accomodate matter effects and gravitational contributions to neutrino spin precession in the presence of a magnetic field. By way of illustration, we perform explicit calculations in the Schwarzschild geometry. In this case, gravitational effects on neutrino oscillations are intimately related to the redshift. We discuss how spacetime curvature could affect the resonance position and adiabaticity of matter-enhanced neutrino flavor conversion.

Recent determinations of the deuterium abundance, $^2$H/H, in high redshift Lyman limit hydrogen clouds challenge the usual picture of primordial nucleosynthesis based on \lq\lq concordance\rq\rq\ of the calculated light element ($^2$H, $^3$He, $^4$He, $^7$Li) nucleosynthesis yields with the observationally-inferred abundances of these species. Concordance implies that all light element yields can be made to agree with the observationally-inferred abundances (within errors) for single global specifications of the baryon-to-photon ratio, $\eta$; lepton number; neutron lifetime; and expansion rate (or equivalently, effective number of light neutrino degrees of freedom $N_{\nu} $). Though one group studying Lyman limit systems obtains a high value of $^2$H/H ($\sim 2\times {10}^{-4}$), another group finds consistently low values ($\sim 2\times {10}^{-5}$). In the former case, concordance for $N_{\nu} =3$ is readily attained for the current observationally-inferred abundances of $^4$He and $^7$Li. But if the latter case represents the primordial deuterium abundance, then concordance for {\it any} $N_{\nu}$ is impossible unless the primordial value of $^7$Li/H is considerably larger than the abundance of lithium as measured in old, hot Pop II halo stars. Furthermore, concordance with $N_{\nu}=3$ is possible for low $^2$H/H only if either (1) the primordial $^4$He abundance has been significantly underestimated, or (2) new neutrino sector physics is invoked. We argue that systematic underestimation of both the $^7$Li and $^4$He primordial abundances is the likely resolution of this problem, a conclusion which is strengthened by new results on $^4$He.

The reported signal at the LSND experiment, when interpreted as neutrino mixing with $\delta m^2 = 6 \;\rm{eV}^2$, provides evidence for neutrinos with a cosmologically significant mass. However, attempts to reconcile this interpretation of the experiment with other hints about neutrino properties require a (sterile) fourth neutrino and/or an ``inverted'' neutrino mass hierarchy. An interpretation of the LSND experiment employing $\delta m^2 = 0.3\; \rm{eV}^2$, with three-generation mixing and a ``normal'' neutrino mass hierarchy, can just barely be reconciled with the negative results of other laboratory neutrino oscillation experiments and the positive hints of neutrino oscillation from the solar and atmospheric neutrino problems. Though subject to test by by future experiments, such a solution allows (but does not demand) neutrino masses relevant for dark matter.

Studies of limits on active-sterile neutrino mixing derived from big bang nucleosynthesis considerations are extended to consider the dependance of these constraints on the primordial deuterium abundance. This study is motivated by recent measurements of D/H in quasar absorption systems, which at present yield discordant results. Limits on active-sterile mixing are somewhat relaxed for high D/H. For low D/H ($\approx 2 \times 10^{-5}$), no active-sterile neutrino mixing is allowed by currently popular upper limits on the primordial $^4$He abundance $Y$. For such low primordial D/H values, the observational inference of active-sterile neutrino mixing by upcoming solar neutrino experiments would imply that $Y$ has been systematically underestimated, unless there is new physics not included in standard BBN.

Two dimensional concordance plots involving the baryon-to-photon ratio, $\eta$, and an effective number of light neutrinos, $N_{\nu}$, are used to discuss the overall consistency of standard big-bang nucleosynthesis in light of recent determinations of the primordial deuterium ($^2$H) abundance. Observations of high-redshift Ly-$\alpha$ clouds have provided discordant $^2$H/H determinations: one cloud with $^2$H/H high compared with the previously accepted upper limit on ($^2$H + $^3$He)/H; and one system with a significantly lower upper bound on $^2$H/H than those previously obtained. The high value of $^2$H/H agrees well with the current observationally-inferred primordial abundances of $^4$He and $^7$Li for $N_{\nu}=3$. The low value of $^2$H/H does not fit well with the current observationally-inferred primordial abundance of $^4$He for $N_{\nu}=3$. In addition, if the low value of $^2$H/H is indicative of the primordial deuterium abundance, then significant depletion of $^7$Li in old, hot Pop II halo stars is probably required to obtain a concordant range of $\eta$, for {\em any} effective number of neutrino flavors. The use of conservative ranges for the primordial abundances of $^2$H, $^4$He, and $^7$Li allow success of the standard picture for $N_{\nu}=3$.

We examine the potential for a ``natural'' three-neutrino mixing scheme to satisfy available data and astrophysical arguments. By ``natural'' we mean no sterile neutrinos, and a neutrino mass hierarchy similar to that of the charged leptons. We seek to satisfy (or solve): 1. Accelerator and reactor neutrino oscillation constraints, including LSND; 2. The atmospheric muon neutrino deficit problem; 3. The solar neutrino problem; 4. Supernova $r$-process nucleosynthesis in neutrino-heated supernova ejecta; 5. Cold+hot dark matter models. We argue that putative supernova $r$-process nucleosynthesis bounds on two-neutrino flavor mixing can be applied directly to three-neutrino mixing in the case where one vacuum neutrino mass eigenvalue difference dominates the others. We show that in this ``one mass scale dominance'' limit, a natural three-neutrino oscillation solution meeting all the above constraints exists only if the atmospheric neutrino data {\em and} the LSND data can be explained with one neutrino mass difference. In this model, an explanation for the solar neutrino data can be effected by employing the {\em other} independent neutrino mass difference. Such a solution is only marginally allowed by the current data, and proposed long-baseline neutrino oscillation experiments can definitively rule it out. If it were ruled out, the simultaneous solution of the above constraints by neutrino oscillations would then require sterile neutrinos and/or a neutrino mass hierarchy of a different nature than that of the charged leptons.