A comprehensive introduction to the theory of the solar neutrino problem is given that is aimed at instructors who are not experts in quantum field theory but would like to incorporate these ideas into instruction of advanced undergraduate or beginning graduate students in courses like astrophysics or quantum mechanics, it is also aimed at the inquisitive student who would like to learn this topic on their own. The presentation assumes as theoretical preparation only that the reader is familiar with the basics of quantum mechanics in Dirac notation and elementary differential equations and matrices.

The recent discovery of merging black holes suggests that a stochastic gravitational-wave background is within reach of the advanced detector network operating at design sensitivity. However, correlated magnetic noise from Schumann resonances threatens to contaminate observation of a stochastic background. In this paper, we report on the first effort to eliminate intercontinental correlated noise from Schumann resonances using Wiener filtering. Using magnetometers as proxies for gravitational-wave detectors, we demonstrate as much as a factor of two reduction in the coherence between magnetometers on different continents. While much work remains to be done, our results constitute a proof-of-principle and motivate follow-up studies with a dedicated array of magnetometers.

We demonstrate the first implementation of recently-developed fast explicit kinetic integration algorithms on modern graphics processing unit (GPU) accelerators. Taking as a generic test case a Type Ia supernova explosion with an extremely stiff thermonuclear network having 150 isotopic species and 1604 reactions coupled to hydrodynamics using operator splitting, we demonstrate the capability to solve of order 100 realistic kinetic networks in parallel in the same time that standard implicit methods can solve a single such network on a CPU. This orders-of-magnitude decrease in compute time for solving systems of realistic kinetic networks implies that important coupled, multiphysics problems in various scientific and technical fields that were intractible, or could be simulated only with highly schematic kinetic networks, are now computationally feasible.

We show that, even for extremely stiff systems, explicit integration may compete in both accuracy and speed with implicit methods if algebraic methods are used to stabilize the numerical integration. The required stabilizing algebra depends on whether the system is well-removed from equilibrium or near equilibrium. This paper introduces a quantitative distinction between these two regimes and addresses the former case in depth, presenting explicit asymptotic methods appropriate when the system is extremely stiff but only weakly equilibrated. A second paper examines quasi-steady-state methods as an alternative to asymptotic methods in systems well away from equilibrium and a third paper extends these methods to equilibrium conditions in extremely stiff systems using partial equilibrium methods. All three papers present systematic evidence for timesteps competitive with implicit methods. Because an explicit method can execute a timestep faster than an implicit method, algebraically-stabilized explicit algorithms might permit integration of larger networks than have been feasible before in various disciplines.

In contrast to the prevailing view in the literature, it is shown that even extremely stiff sets of ordinary differential equations may be solved efficiently by explicit methods if limiting algebraic solutions are used to stabilize the numerical integration. The stabilizing algebra differs essentially for systems well-removed from equilibrium and those near equilibrium. Explicit asymptotic and quasi-steady-state methods that are appropriate when the system is only weakly equilibrated are examined first. These methods are then extended to the case of close approach to equilibrium through a new implementation of partial equilibrium approximations. Using stringent tests with astrophysical thermonuclear networks, evidence is provided that these methods can deal with the stiffest networks, even in the approach to equilibrium, with accuracy and integration timestepping comparable to that of implicit methods. Because explicit methods can execute a timestep faster and scale more favorably with network size than implicit algorithms, our results suggest that algebraically-stabilized explicit methods might enable integration of larger reaction networks coupled to fluid dynamics than has been feasible previously for a variety of disciplines.

A preceding paper demonstrated that explicit asymptotic methods generally work much better for extremely stiff reaction networks than has previously been shown in the literature. There we showed that for systems well removed from equilibrium explicit asymptotic methods can rival standard implicit codes in speed and accuracy for solving extremely stiff differential equations. In this paper we continue the investigation of systems well removed from equilibrium by examining quasi-steady-state (QSS) methods as an alternative to asymptotic methods. We show that for systems well removed from equilibrium, QSS methods also can compete with, or even exceed, standard implicit methods in speed, even for extremely stiff networks, and in many cases give somewhat better integration speed than for asymptotic methods. As for asymptotic methods, we will find that QSS methods give correct results, but with non-competitive integration speed as equilibrium is approached. Thus, we shall find that both asymptotic and QSS methods must be supplemented with partial equilibrium methods as equilibrium is approached to remain competitive with implicit methods.

In two preceding papers we have shown that, when reaction networks are well-removed from equilibrium, explicit asymptotic and quasi-steady-state approximations can give algebraically-stabilized integration schemes that rival standard implicit methods in accuracy and speed for extremely stiff systems. However, we also showed that these explicit methods remain accurate but are no longer competitive in speed as the network approaches equilibrium. In this paper we analyze this failure and show that it is associated with the presence of fast equilibration timescales that neither asymptotic nor quasi-steady-state approximations are able to remove efficiently from the numerical integration. Based on this understanding, we develop a partial equilibrium method to deal effectively with the approach to equilibrium and show that explicit asymptotic methods, combined with the new partial equilibrium methods, give an integration scheme that plausibly can deal with the stiffest networks, even in the approach to equilibrium, with accuracy and speed competitive with that of implicit methods. Thus we demonstrate that such explicit methods may offer alternatives to implicit integration of even extremely stiff systems, and that these methods may permit integration of much larger networks than have been possible before in a number of fields.

Proton capture by 17F plays an important role in the synthesis of nuclei in nova explosions. A revised rate for this reaction, based on a measurement of the 1H(17F,p)17F excitation function using a radioactive 17F beam at ORNL's Holifield Radioactive Ion Beam Facility, is used to calculate the nucleosynthesis in nova outbursts on the surfaces of 1.25 and 1.35 solar mass ONeMg white dwarfs and a 1.00 solar mass CO white dwarf. We find that the new 17F(p,gamma)18Ne reaction rate changes the abundances of some nuclides (e.g., 17O) synthesized in the hottest zones of an explosion on a 1.35 solar mass white dwarf by more than a factor of 10,000 compared to calculations using some previous estimates for this reaction rate, and by more than a factor of 3 when the entire exploding envelope is considered. In a 1.25 solar mass white dwarf nova explosion, this new rate changes the abundances of some nuclides synthesized in the hottest zones by more than a factor of 600, and by more than a factor of 2 when the entire exploding envelope is considered. Calculations for the 1.00 solar mass white dwarf nova show that this new rate changes the abundance of 18Ne by 21%, but has negligible effect on all other nuclides. Comparison of model predictions with observations is also discussed.

Based on available experimental data, a new set of Nilsson parameters is proposed for proton-rich nuclei with proton or neutron numbers $28\leq N\leq 40$. The resulting single-particle spectra are compared with those from relativistic and non-relativistic mean field theories. Collective excitations in some even--even proton-rich nuclei in the upper $pf$ shell are investigated using the Projected Shell Model with the new Nilsson basis. It is found that the regular bands are sharply disturbed by band crossings involving $1g_{9/2}$ neutrons and protons. Physical quantities for exploring the nature of the band disturbance and the role of the $1g_{9/2}$ single-particle are predicted, which may be tested by new experiments with radioactive beams.

We compare Newtonian three-flavor multigroup Boltzmann (MGBT) and (Bruenn's) multigroup flux-limited diffusion (MGFLD) neutrino transport in postbounce core collapse supernova environments. We focus our study on quantities central to the postbounce neutrino heating mechanism for reviving the stalled shock. Stationary-state three-flavor neutrino distributions are developed in thermally and hydrodynamically frozen time slices obtained from core collapse and bounce simulations that implement Lagrangian hydrodynamics and MGFLD neutrino transport. Most important, we find, for a region above the gain radius, net heating rates for MGBT that are as much as ~2 times the corresponding MGFLD rates, and net cooling rates below the gain radius that are typically ~0.8 times the MGFLD rates. These differences stem from differences in the neutrino luminosities and mean inverse flux factors, which can be as much as 11% and 24%, respectively. They are greatest at earlier postbounce times for a given progenitor mass and, for a given postbounce time, greater for greater progenitor mass. We discuss the ramifications these new results have for the supernova mechanism.

Core-collapse supernovae are, despite their spectacular visual display, neutrino events. Virtually all of the 10^53 ergs of gravitational binding energy released in the formation of the nascent neutron star is carried away in the form of neutrinos and antineutrinos of all three flavors, and these neutrinos are primarily responsible for powering the explosion. This mechanism depends sensitively on the neutrino transport between the neutrinospheres and the shock. In light of this, we have performed a comparison of multigroup Boltzmann neutrino transport (MGBT) and multigroup flux-limited diffusion (MGFLD) in post-core bounce environments. Differences in the mean inverse flux factors, luminosities, and RMS energies translate to heating rates that are up to 2 times larger for Boltzmann transport, with net cooling rates below the gain radius that are typically 0.8 times the MGFLD rates. These differences are greatest at earlier postbounce times for a given progenitor mass, and for a given postbounce time, greater for greater progenitor mass. The increased differences with increased progenitor mass suggest that the net heating enhancement from MGBT is potentially robust and self-regulated.

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 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. 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.