At the present time even the most sophisticated, multi-dimensional simulations of core-collapse supernovae do not (self-consistently) include neutrino flavor transformation. This physics is missing despite the importance of neutrinos in the core-collapse explosion paradigm. Because of this dependence, any flavor transformation that occurs in the region between the proto-neutron star and the shock could result in major effects upon the dynamics of the explosion. We present the first hydrodynamic core-collapse supernova simulation which simultaneously includes flavor transformation of the free-streaming neutrinos in the neutrino transport. These oscillation calculations are dynamically updated and evolve self-consistently alongside the hydrodynamics. Using a $M=20\;{\rm M_{\odot}}$ progenitor, we find that while the oscillations have an effect on the neutrino emission and the heating rates, flavor transformation alone does not lead to a successful explosion of this progenitor in spherical symmetry.

Neutrinos play a critical role of transporting energy and changing the lepton density within core-collapse supernovae and neutron star mergers. The quantum kinetic equations (QKEs) combine the effects of neutrino-matter interactions treated in classical Boltzmann transport with the neutrino flavor-changing effects treated in neutrino oscillation calculations. We present a method for extending existing neutrino interaction rates to full QKE source terms for use in numerical calculations. We demonstrate the effects of absorption and emission by nucleons and nuclei, electron scattering, electron-positron pair annihilation, nucleon-nucleon bremsstrahlung, neutrino-neutrino scattering. For the first time, we include all these collision terms self-consistently in a simulation of the full isotropic QKEs in conditions relevant to core-collapse supernovae and neutron star mergers. For our choice of parameters, the long-term evolution of the neutrino distribution function proceeds similarly with and without the oscillation term, though with measurable differences. We demonstrate that electron scattering, nucleon-nucleon bremsstrahlung processes, and four-neutrino processes dominate flavor decoherence in the protoneutron star (PNS), absorption dominates near the shock, and all of the considered processes except elastic nucleon scattering are relevant in the decoupling region. Finally, we propose an effective decoherence opacity that at most energies predicts decoherence rates to within a factor of 10 in our model PNS and within 20% outside of the PNS.

It has recently been demonstrated analytically that the two-point correlation function for pairs of neutrinos may contain information about the size of the proto-neutron star formed in a Galactic core-collapse supernova. The information about the size of the source emerges via the neutrino equivalent of intensity interferometry originally used by Hanbury-Brown and Twiss with photons to measure the radii of stars. However the analytic demonstration of neutrino intensity interferometry with supernova neutrinos made a number of approximations: that the two neutrinos had equal energies, the neutrinos were emitted at simultaneous times from two points and were detected simultaneously at two detection points that formed a plane with the emission points. These approximations need to be relaxed in order to better determine the feasibility of neutrino intensity interferometry for supernovae neutrinos in a more realistic scenario. In this paper we further investigate the feasibility of intensity interferometry for supernova neutrinos by relaxing all the approximations made in the earlier study. We find that, while relaxing any one assumption reduces the correlation signal, the relaxation of the assumption of equal times of detection is by far the largest detrimental factor. For neutrino energies of order $\sim$15 MeV and a supernova distance of L = 10 kpc, we show that in order to observe the interference pattern in the two-point correlation function of the neutrino pairs, the timing resolution of a detector needs to be on the order of $\lesssim 10^{-21}\;{\rm s}$ if the initial neutrino wave packet has a size of $\sigma_x \sim 10^{-11}\;{\rm cm}$.

We explore the possibility of probing the nonstandard interactions between the neutrino and a hypothetical massive scalar or pseudoscalar via neutrino flavor transformation in supernovae. We find that in the ultrarelativistic limit, the effective interaction between the neutrinos vanishes if neutrinos are Dirac fermions but not if they are Majorana fermions. The impact of the new neutrino interaction upon the flavor transformation above the neutrinosphere is calculated in the context of the multi-angle "neutrino bulb model". We find that the addition of the nonstandard neutrino self-interaction (NSSI) to the ordinary V-A self-interaction between neutrinos is capable of dramatically altering the collective oscillations when its strength is comparable to the standard, V-A, interaction. The effect of flavor-preserving (FP) NSSI is generally to suppress flavor transformation, while the flavor-violating (FV) interactions are found to promote flavor transformations. If the neutrino signal from a Galactic supernova can be sufficiently well understood, supernova neutrinos can provide complimentary constraints on scalar/pseudoscalar interactions of neutrinos as well as distinguishing whether the neutrino is a Majorana or Dirac fermion.

This year marks the thirtieth anniversary of the only supernova from which we have detected neutrinos - SN 1987A. The twenty or so neutrinos that were detected were mined to great depth in order to determine the events that occurred in the explosion and to place limits upon all manner of neutrino properties. Since 1987 the scale and sensitivity of the detectors capable of identifying neutrinos from a Galactic supernova have grown considerably so that current generation detectors are capable of detecting of order ten thousand neutrinos for a supernova at the Galactic Center. Next generation detectors will increase that yield by another order of magnitude. Simultaneous with the growth of neutrino detection capability, our understanding of how massive stars explode and how the neutrino interacts with hot and dense matter has also increased by a tremendous degree. The neutrino signal will contain much information on all manner of physics of interest to a wide community. In this review we describe the expected features of the neutrino signal, the detectors which will detect it, and the signatures one might try to look for in order to get at these physics.

Intensity interferometry is a technique that has been used to measure the size of sources ranging from the quark-gluon plasma formed in heavy ion collisions to the radii of stars. We investigate using the same technique to measure proto-neutron star (PNS) radii with the neutrino signal received from a core-collapse supernovae. Using a full wave-packet analysis, including the neutrino mass for the first time, we derive criteria where the effect can be expected to provide the desired signal, and find that neutrinos from the next Galactic supernova should contain extractable PNS radius information.

A very massive star with a carbon-oxygen core in the range of $64$ M$_{\odot}<M_{\mathrm{CO}}<133$ M$_{\odot}$ is expected to undergo a very different kind of explosion known as a pair instability supernova. Pair instability supernovae are candidates for superluminous supernovae due to the prodigious amounts of radioactive elements they create. While the basic mechanism for the explosion is understood, how a star reaches a state is not, thus observations of a nearby pair-instability supernova would allow us to test current models of stellar evolution at the extreme of stellar masses. Much will be sought within the electromagnetic radiation we detect from such a supernova but we should not forget that the neutrinos from a pair-instability supernova contain unique signatures of the event that unambiguously identify this type of explosion. We calculate the expected neutrino flux at Earth from two, one-dimensional pair-instability supernova simulations which bracket the mass range of stars which explode by this mechanism taking into account the full time and energy dependence of the neutrino emission and the flavor evolution through the outer layers of the star. We calculate the neutrino signals in five different detectors chosen to represent present or near future designs. We find the more massive progenitors explode as pair-instability supernova which can easily be detected in multiple different neutrino detectors at the `standard' supernova distance of $10\;{\rm kpc}$ producing several events in DUNE, JUNE and SuperKamiokande, while the lightest progenitors only produce a handful of events (if any) in the same detectors. The proposed HyperKamiokande detector would detect neutrinos from a large pair-instability supernova as far as $\sim 50\;{\rm kpc}$ allowing it to reach the Megallanic Clouds and the several very high mass stars known to exist there.

The strong gravitational field around a proto-neutron star can modify the neutrino flavor transformations that occur above the neutrinosphere via three General Relativistic (GR) effects: time dilation, energy redshift, and trajectory bending. Depending on the compactness of the central object, the neutrino self-interaction potential is up to three times as large as that without GR principally due to trajectory bending which increases the intersection angles between different neutrino trajectories, and time dilation which changes the fluxes. We determine whether GR effects are important for flavor transformation during the different epochs of a supernova by using multi-angle flavor transformation calculations and consider a density profile and neutrino spectra representative of both the accretion and cooling phases. We find the GR effects are smaller during the accretion phase due to low compactness of the proto-neutron star and merely delay the decoherence; the neutrino bipolar oscillations during the cooling phase are also delayed due to the GR effects but the delay may be more important because the delay occurs at radii where it might alter the nucleosynthesis in the neutrino driven wind.

A neutrino subject to an external, time-dependent perturbing potential can be forced to make transitions between its flavor states. A neutrino with three (or more) flavors can exhibit phenomena that cannot occur if the neutrino had just two. We present an approximate analytic solution for the temporal evolution of a multi-flavor neutrino in response to an arbitrary perturbing Hamiltonian that has been decomposed into its Fourier modes. We impose no restriction upon the number of flavors nor upon the structure of the perturbing Hamiltonian, the number of Fourier modes, their amplitude or their frequencies. We apply the theory to study three-flavor neutrino transformation due to perturbations built from two and three Fourier modes. For the case of two Fourier modes we observe the equivalent of "induced transparency" from quantum optics whereby transitions between a given pair of states are suppressed due to the presence of a resonant mode between another pair. When we add a third Fourier mode we find a new effect whereby the third mode can manipulate the transition probabilities of the two mode case so as to force complete transparency or, alternatively, restore "opacity" meaning the perturbative potential regains its ability to induce neutrino transitions. In both applications we demonstrate how the analytic solutions are able to match the amplitude and wavenumber of the numerical results to within a few percent.

During the accretion phase of a core-collapse supernovae, large amplitude turbulence is generated by the combination of the standing accretion shock instability and convection driven by neutrino heating. The turbulence directly affects the dynamics of the explosion, but there is also the possibility of an additional, indirect, feedback mechanism due to the effect turbulence can have upon neutrino flavor evolution and thus the neutrino heating. In this paper we consider the effect of turbulence during the accretion phase upon neutrino evolution, both numerically and analytically. Adopting representative supernova profiles taken from the accretion phase of a supernova simulation, we find the numerical calculations exhibit no effect from turbulence. We explain this absence using two analytic descriptions: the Stimulated Transition model and the Distorted Phase Effect model. In the Stimulated Transition model turbulence effects depend upon six different lengthscales, and three criteria must be satisfied between them if one is to observe a change in the flavor evolution due to Stimulated Transition. We further demonstrate that the Distorted Phase Effect depends upon the presence of multiple semi-adiabatic MSW resonances or discontinuities that also can be expressed as a relationship between three of the same lengthscales. When we examine the supernova profiles used in the numerical calculations we find the three Stimulated Transition criteria cannot be satisfied, independent of the form of the turbulence power spectrum, and that the same supernova profiles lack the multiple semi-adiabatic MSW resonances or discontinuities necessary to produce a Distorted Phase Effect. Thus we conclude that even though large amplitude turbulence is present in supernova during the accretion phase, it has no effect upon neutrino flavor evolution.

Despite their use as cosmological distance indicators and their importance in the chemical evolution of galaxies, the unequivocal identification of the progenitor systems and explosion mechanism of normal type Ia supernova (SN Ia) remains elusive. The leading hypothesis is that such a supernova is a thermonuclear explosion of a carbon-oxygen white dwarf, but the exact explosion mechanism is still a matter of debate. Observation of a galactic SN Ia would be of immense value in answering the many open questions related to these events. One potentially useful source of information about the explosion mechanism and progenitor is the neutrino signal. In this paper we compute the expected neutrino signal from a gravitationally confined detonation (GCD) explosion scenario for a SN~Ia and show how the flux at Earth contains features in time and energy unique to this scenario. We then calculate the expected event rates in the Super-K, Hyper-K, JUNO, DUNE, and IceCube detectors and find both Hyper-K and IceCube would see a few events for a GCD supernova at 1 kpc or closer, while Super-K, JUNO, and DUNE would see a events if the supernova were closer than ${\sim}0.3$ kpc. The distance and detector criteria needed to resolve the time and spectral features arising from the explosion mechanism, neutrino production, and neutrino oscillation processes are also discussed. The neutrino signal from the GCD is then compared with the signal from a deflagration-to-detonation transition (DDT) explosion model computed previously. We find the overall event rate is the most discriminating feature between the two scenarios followed by the event rate time structure. Using the event rate in the Hyper-K detector alone, the DDT can be distinguished from the GCD at 2$\sigma$ if the distance to the supernova is less than $2.3\;{\rm kpc}$ for a normal mass ordering and $3.6\;{\rm kpc}$ for an inverted ordering.

Nonstandard interactions (NSI) of neutrinos with matter can significantly alter neutrino flavor evolution in supernovae with the potential to impact explosion dynamics, nucleosynthesis, and the neutrinos signal. In this paper, we explore, both numerically and analytically, the landscape of neutrino flavor transformation effects in supernovae due to NSI and find a new, heretofore unseen transformation processes can occur. These new transformations can take place with NSI strengths well below current experimental limits. Within a broad swath of NSI parameter space, we observe symmetric and standard matter-neutrino resonances for supernovae neutrinos, a transformation effect previously only seen in compact object merger scenarios; in another region of the parameter space we find the NSI can induce neutrino collective effects in scenarios where none would appear with only the standard case of neutrino oscillation physics; and in a third region the NSI can lead to the disappearance of the high density Mikheyev-Smirnov-Wolfenstein resonance. Using a variety of analytical tools, we are able to describe quantitatively the numerical results allowing us to partition the NSI parameter according to the transformation processes observed. Our results indicate nonstandard interactions of supernova neutrinos provide a sensitive probe of beyond the Standard Model physics complementary to present and future terrestrial experiments.

The DUNE/LBNF program aims to address key questions in neutrino physics and astroparticle physics. Realizing DUNE's potential to reconstruct low-energy particles in the 10-100 MeV energy range will bring significant benefits for all DUNE's science goals. In neutrino physics, low-energy sensitivity will improve neutrino energy reconstruction in the GeV range relevant for the kinematics of DUNE's long-baseline oscillation program. In astroparticle physics, low-energy capabilities will make DUNE's far detectors the world's best apparatus for studying the electron-neutrino flux from a supernova. This will open a new window to unrivaled studies of the dynamics and neutronization of a star's central core in real time, the potential discovery of the neutrino mass hierarchy, provide new sensitivity to physics beyond the Standard Model, and evidence of neutrino quantum-coherence effects. The same capabilities will also provide new sensitivity to `boosted dark matter' models that are not observable in traditional direct dark matter detectors.

It has long been recognized that the neutrinos detected from the next core-collapse supernova in the Galaxy have the potential to reveal important information about the dynamics of the explosion and the nucleosynthesis conditions as well as allowing us to probe the properties of the neutrino itself. The neutrinos emitted from thermonuclear - type Ia - supernovae also possess the same potential, although these supernovae are dimmer neutrino sources. For the first time, we calculate the time, energy, line of sight, and neutrino-flavor-dependent features of the neutrino signal expected from a three-dimensional delayed-detonation explosion simulation, where a deflagration-to-detonation transition triggers the complete disruption of a near-Chandrasekhar mass carbon-oxygen white dwarf. We also calculate the neutrino flavor evolution along eight lines of sight through the simulation as a function of time and energy using an exact three-flavor transformation code. We identify a characteristic spectral peak at $\sim 10$ MeV as a signature of electron captures on copper. This peak is a potentially distinguishing feature of explosion models since it reflects the nucleosynthesis conditions early in the explosion. We simulate the event rates in the Super-K, Hyper-K, JUNO, and DUNE neutrino detectors with the SNOwGLoBES event rate calculation software and also compute the IceCube signal. Hyper-K will be able to detect neutrinos from our model out to a distance of $\sim 10$ kpc. At 1 kpc, JUNO, Super-K, and DUNE would register a few events while IceCube and Hyper-K would register several tens of events.

The neutrinos emitted from the proto-neutron star created in a core-collapse supernova must run through a significant amount of turbulence before exiting the star. Turbulence can modify the flavor evolution of the neutrinos imprinting itself upon the signal detected here at Earth. The turbulence effect upon individual neutrinos, and the correlation between pairs of neutrinos, might exhibit sensitivity to the power spectrum of the turbulence and recent analysis of the turbulence in a two-dimensional hydrodynamical simulation of a core-collapse supernova indicates the power spectrum may not be the Kolmogorov 5/3 inverse power law as has been previously assumed. In this paper we study the effect of non-Kolmogorov turbulence power spectra upon neutrinos from a point source as a function of neutrino energy and turbulence amplitude at a fixed post-bounce epoch. We find the two effects of turbulence upon the neutrinos - the distorted phase effect and the stimulated transitions - both possess strong and weak limits in which dependence upon the power spectrum is absent or evident respectively. Since neutrinos of a given energy will exhibit these two effects at different epochs of the supernova each with evolving strength, we find there is sensitivity to the power spectrum present in the neutrino burst signal from a Galactic supernova.

On February 23, 1987 we collected 24 neutrinos from the explosion of a blue super-giant star in the Large Magellanic Cloud confirming the basic paradigm of core-collapse supernova. During the many years we have been waiting for a repeat of that momentous day, the number and size of neutrino detectors around the world has grown considerably. If the neutrinos from the next supernova in our Galaxy arrive tomorrow we shall collect upwards of tens of thousands of events and next generation detectors will increase the amount of data we collect by more than an order of magnitude. But it is also now apparent that the message is much more complex than previously thought because many time, energy and neutrino flavor dependent features are imprinted upon the signal either at emission or by the passage through the outer layers of the star. These features arise due to the explosion dynamics, the physics of nuclei at high temperatures and densities, and the properties of neutrinos. In this proceedings I will present some aspects of the physics of supernova neutrino oscillations and what we should expect to observe when the neutrinos from the next Galactic supernova (eventually) arrive.

We apply the model of stimulated neutrino transitions to neutrinos traveling through turbulence on a non-constant density profile. We describe a method to predict the location of large amplitude transitions and demonstrate the effectiveness of this method by comparing to numerical calculations using a model supernova (SN) profile. The important wavelength scales of turbulence, both those that stimulate neutrino transformations and those that suppress them, are presented and discussed. We then examine the effects of changing the parameters of the turbulent spectrum, specifically the root-mean-square amplitude and cutoff wavelength, and show how the stimulated transitions model offers an explanation for the increase in both the amplitude and number of transitions with large amplitude turbulence, as well as a suppression or absence of transitions for long cutoff wavelengths. The method can also be used to predict the location of transitions between antineutrino states which, in the normal hierarchy we are using, will not undergo Mikheev-Smirnov-Wolfenstein (MSW) transitions. Finally, the stimulated neutrino transitions method is applied to the turbulence derived found in a 2D supernova simulation and explains the minimal observed effect on neutrino oscillations in the simulation as as being due to excessive long wavelength modes suppressing transitions and the absence of modes that fulfill the parametric resonance condition.

(Abridged) In order to decode the neutrino burst signal from a Galactic core-collapse supernova and reveal the complicated inner workings of the explosion we need a thorough understanding of the neutrino flavor evolution from the proto-neutron star outwards. The flavor content of the signal evolves due to both neutrino collective effects and matter effects which can lead to a highly interesting interplay and distinctive spectral features. In this paper we investigate the supernova neutrino flavor evolution in three different progenitors and include collective flavor effects, the evolution of the Mikheyev, Smirnov & Wolfenstein conversion due to the shock wave passage through the star, and the impact of turbulence. In the Oxygen-Neon-Magnesium supernova we find that the impact of turbulence is both brief and slight during a window of 1-2 seconds post bounce. Thus the spectral features of collective and shock effects in the neutrino signals from ONeMg supernovae may be almost turbulence free making them the easiest to interpret. For the more massive progenitors we again find that small amplitude turbulence, up to 10%, leads to a minimal modification of the signal, and the emerging neutrino spectra retain both collective and MSW features. However, when larger amounts of turbulence is added, 30% and 50%, the features of collective and shock wave effects in the high density resonance channel are almost completely obscured at late times. Yet at the same time we find the other mixing channels - the low density resonance channel and the non-resonant channels - begin to develop turbulence signatures. Large amplitude turbulent motions in the outer layers of more massive, iron core-collapse supernovae may obscure the most obvious fingerprints of collective and shock wave effects in the neutrino signal but cannot remove them completely, and additionally bring about new features in the signal.

Large amplitude oscillations between the states of a quantum system can be stimulated by sinusoidal external potentials with frequencies that are similar to the energy level splitting of the states or a fraction thereof. Situations when the applied frequency is equal to an integer fraction of the energy level splittings are known as parametric resonances. We investigate this effect for neutrinos both analytically and numerically for the case of arbitrary numbers of neutrino flavors. We look for environments where the effect may be observed and find that supernova are the one realistic possibility due to the necessity of both large densities and large amplitude fluctuations. The comparison of numerical and analytic results of neutrino propagation through a model supernova reveals it is possible to predict the locations and strengths of the stimulated transitions that occur.

The set of transition probabilities for a single neutrino emitted from a point source after passage through a turbulent supernova density profile have been found to be random variates drawn from parent distributions whose properties depend upon the stage of the explosion, the neutrino energy and mixing parameters, the observed channel, and the properties of the turbulence such as the amplitude C*. In this paper we examine the consequences of the recently measured mixing angle \theta_13 upon the neutrino flavor transformation in supernova when passing through turbulence. We find the measurements of a relatively large value of \theta_13 means the neutrinos are relatively immune to small, C* < 1%, amplitude turbulence but as C* increases the turbulence effects grow rapidly and spread to all mixing channels. For C* > 10% the turbulence effects in the high (H) density resonance mixing channels are independent of \theta_13 but non-resonant mixing channels are more sensitive to turbulence when \theta_13 is large.

Black hole accretion disks can form through the collapse of rotating massive stars. These disks produce large numbers of neutrinos and antineutrinos of electron flavor that can influence energetics and nucleosynthesis. Neutrinos are produced in sufficient numbers that, after they are emitted, they can undergo flavor transformation facilitated by the neutrino self interaction. We show that some of the neutrino flavor transformation phenomenology for accretion disks is similar to that of the supernova case, but also, we find the disk geometry lends itself to different transformation behaviors. These transformations strongly influence the nucleosynthetic outcome of disk winds.

We consider the Hamiltonian for neutrino oscillations in matter in the case of arbitrary potentials including off-diagonal complex terms. We derive the 'matter basis' Hamiltonian in terms of quantities one can derive from the flavor basis Hamiltonian and its derivative, for an arbitrary number of neutrino flavors. We make our expressions explicit for the two-neutrino flavor case and apply our results to the situation where the Hamiltonian includes both coupling to matter and to neutrinos, which describes neutrino propagation in core-collapse supernovae. We show that the neutrino flavour evolution depends on the mixing matrix derivatives involving not only the derivative of the matter mixing angles but also of the phases. In particular, we point out the important role of the phase derivatives, that appear due to the neutrino-neutrino interaction, and show how it can cause an oscillating degeneracy between the diagonal elements of the matter basis Hamiltonian. Our results also reveal, that the end of the synchronization regime is due to a rapid increase of the phase derivative, and identify the condition to be fulfilled for the onset of bipolar oscillations involving both the off-diagonal neutrino-neutrino interaction contributions and the vacuum terms.

It is becoming ever clearer that the neutrino signal from the next supernova in our Galaxy can reveal missing information about the neutrino as well as allowing us to probe the explosion of the star by decoding the temporal and spectral evolution of the flavor composition of the signal. But this information may be lost if turbulence in the supernova `depolarizes' the neutrinos so that the observed flux for each flavor is an equal mixture of the initial - unencoded - spectra. Determining if depolarization occurs is one of the most pressing issues of this field. The most difficult aspect of studying the effect of turbulence upon the neutrinos is the lack of any theoretical models that allow us to understand the results of numerical studies. This paper makes the suggestion that Random Matrix Theory (RMT) may shine some light in this direction and presents support for this the possibility by comparing the distribution of crossing and survival probabilities obtained numerically for some `test case' calculations with the distributions one expects from RMT in the calculable limit of depolarization of N neutrino flavors.

Multi-dimensional core-collapse supernova simulations exhibit turbulence of large amplitude and over large scales. As neutrinos pass through the supernova mantle the turbulence is expected to modify their evolution compared to the case where the explosion is free of turbulence. In this paper we study this turbulence effect upon the neutrinos modelling the turbulence expected from multi-dimensional simulations by adding matter density fluctuations to density profiles taken from one-dimensional hydrodynamical simulations. We investigate the impact upon the supernova neutrino transition probabilities as a function of the neutrino mixing angle theta_13 and turbulence amplitude. In the high (H) resonant channel and with large theta_13 values we find that turbulence is effectively two flavor for fluctuation amplitudes <~ 1% and have identified a new effect due to the combination of turbulence and multiple H resonances that leads to a sensitivity to fluctuations amplitudes as small as ~ 0.001%. At small values of theta_13, beyond the range achievable in Earth based experiments, we find that turbulence leads to new flavor transient effects in the channel where the MSW H resonance occurs. Finally, we investigate large amplitude fluctuations which lead to three flavor effects due to broken HL factorization and significant non-resonant transitions and identify two non-resonant turbulence effects, one depending on the theta_13, and the other independent of this angle and due to the low (L) MSW resonance.

We present the first calculations with three flavors of collective and shock wave effects for neutrino propagation in core-collapse supernovae using hydroynamical density profiles and the S matrix formalism. We explore the interplay between the neutrino-neutrino interaction and the effects of multiple resonances upon the time signal of positrons in supernova observatories. A specific signature is found for the inverted hierarchy and a large third neutrino mixing angle and we predict, in this case, a dearth of lower energy positrons in Cherenkov detectors midway through the neutrino signal and the simultaneous revelation of valuable information about the original fluxes. We show that this feature is also observable with current generation neutrino detectors at the level of several sigmas.

Rapid progress has been made during recent years in the understanding of the flavour oscillations that occur as neutrinos traverse through supernova. The previous paradigm has given way and it is now clear that the neutrino signals we shall receive from future Galactic supernovae will allow us both to peer inside these extraordinary cosmic events and to probe some of the fundamental properties of these elusive particles. In this review we aim to distill the progress that has been made focusing upon the effects of the dynamic density profile and the emergence of collective flavour oscillations due to neutrino self-interactions.

The neutrinos detected from the next Galactic core-collapse supernova will contain valuable information on the internal dynamics of the explosion. One mechanism leading to a temporal evolution of the neutrino signal is the variation of the induced neutrino flavor mixing driven by changes in the density profile. With one and two dimensional hydrodynamical simulations we identify the behavior and properties of prominent features of the explosion. Using these results we demonstrate the time variation of the neutrino crossing probabilities due to changes in the MSW neutrino transformations as the star explodes by using the S-matrix - Monte Carlo - approach to neutrino propagation. After adopting spectra for the neutrinos emitted from the proto-neutron star we calculate for a Galactic supernova the evolution of the positron spectra within a water Cerenkov detector and the ratio of charged current to neutral current event rates for a heavy water - SNO like - detector and find that these detector signals are feasible probes of a number of explosion features.

The central engine that drives gamma ray burst (GRB) explosions may derive from the ability of electrons/positrons and nucleons to tap into the momentum and energy from the large neutrino luminosity emitted by an accretion disk surrounding a black hole. This transfer of momentum and energy occurs due to neutrino absorption, scattering, and annihilation and the non-spherical geometry of the source both increases the annihilation efficiency and, close to the black hole, directs the momentum transfer towards the disk axis. We present annihilation efficiencies and the momentum/energy transfers for a number of accretion disk models and compute the critical densities of infalling material below which the transfer of neutrino momentum/energy will lead to an explosion. Models in which the neutrinos and antineutrinos become trapped within the disk have noticeably different momentum and energy deposition structure compared to thin disk models that may lead to significant differences in the explosion dynamics.

The simplest, `standard' model of Big Bang Nucleosynthesis (SBBN) assumes three light neutrinos (N_nu = 3) and no significant electron neutrino asymmetry, leaving only one adjustable parameter: the baryon to photon ratio eta. The primordial abundance of any one nuclide can, therefore, be used to measure the baryon abundance and the value derived from the observationally inferred primordial abundance of deuterium closely matches that from current, non-BBN data, primarily from the WMAP survey. However, using this same estimate there is a tension between the SBBN-predicted 4He and 7Li abundances and their current, observationally inferred primordial abundances, suggesting that N_nu may differ from the standard model value of three and/or that there may be a non-zero neutral lepton asymmetry (or, that systematic errors in the abundance determinations have been underestimated or overlooked). The differences are not large and the allowed ranges of the BBN parameters permitted by the data are quite small. Within these ranges, the BBN-predicted abundances of D, 3He, 4He, and 7Li are very smooth, monotonic functions of eta, N_nu, and the lepton asymmetry. It is possible to describe the dependencies of these abundances (or powers of them) upon the three parameters by simple, linear fits which, over their ranges of applicability, are accurate to a few percent or better. The fits presented here have not been maximized for their accuracy but, for their simplicity. To identify the ranges of applicability and relative accuracies, they are compared to detailed BBN calculations; their utility is illustrated with several examples. Given the tension within BBN, these fits should prove useful in facilitating studies of the viability of proposals for non-standard physics and cosmology, prior to undertaking detailed BBN calculations.

We have examined the effects of a bound dineutron, n2, upon big bang nucleosynthesis (BBN) as a function of its binding energy B_n2. We find a weakly bound dineutron has little impact but as B_n2 increases its presence begins to alter the flow of free nucleons to helium-4. Due to this disruption, and in the absence of changes to other binding energies or fundamental constants, BBN sets a reliable upper limit of B_n2 <~ 2.5 MeV in order to maintain the agreement with the observations of the primordial helium-4 mass fraction and D/H abundance.

We quantify the extent to which extra relativistic energy density can be concealed by a neutrino asymmetry without conflicting with the baryon asymmetry measured by the Wilkinson Microwave Anisotropy Probe (WMAP). In the presence of a large electron neutrino asymmetry, slightly more than seven effective neutrinos are allowed by Big Bang Nucleosynthesis (BBN) and WMAP at 2\sigma. The same electron neutrino degeneracy that reconciles the BBN prediction for the primordial helium abundance with the observationally inferred value also reconciles the LSND neutrino with BBN by suppressing its thermalization prior to BBN.

Big Bang Nucleosynthesis (BBN) has increasingly become the tool of choice for investigating the permitted variation of fundamental constants during the earliest epochs of the Universe. Here we present a BBN calculation that has been modified to permit changes in the QCD scale, Lambda_QCD. The primary effects of changing the QCD scale upon BBN are through the deuteron binding energy, B_D, and the neutron-proton mass difference, delta-m_np, which both play crucial roles in determining the primordial abundances. In this paper, we show how a simplified BBN calculation allows us to restrict the nuclear data we need to just B_D and delta-m_np yet still gives useful results so that any variation in Lambda_QCD may be constrained via the corresponding shifts in B_D and delta-m_np by using the current estimates of the primordial deuterium abundance and helium mass fraction. The simplification predicts the helium-4 and deuterium abundances to within 1% and 50% respectively when compared with the results of a standard BBN code. But Lambda_QCD also affects much of remaining required nuclear input so this method introduces a systematic error into the calculation and we find a degeneracy between B_D and delta-m_np. We show how increased understanding of the relationship of the pion mass and/or B_D to other nuclear parameters, such as the binding energy of tritium and the cross section of T + D -> He4 + n, would yield constraints upon any change in B_D and delta-m_np at the 10% level.

We place constraints on the number of relativistic degrees of freedom and on the baryon asymmetry at the epoch of Big Bang Nucleosynthesis (BBN) and at recombination, using cosmic background radiation (CBR) data from the Wilkinson Microwave Anisotropy Probe (WMAP), complemented by the Hubble Space Telescope (HST) Key Project measurement of the Hubble constant, along with the latest compilation of deuterium abundances and measurements of the primordial helium abundance. The agreement between the derived values of these key cosmological and particle physics parameters at these widely separated (in time or redshift) epochs is remarkable. From the combination of CBR and BBN data, we find the 2\sigma ranges for the effective number of neutrinos and for the baryon asymmetry (baryon to photon number ratio \eta) to be 1.7-3.0 and 5.53-6.76 \times 10^{-10}, respectively.

A common property of popular models of quintessence dark energy is the convergence to a common solution from a large range of the initial conditions. We re-examine the popular inverse power-law model of quintessence (where the common solution is dubbed as the 'tracker') with particular attention to the initial conditions for the field and their influence on the evolution. We find that previously derived limits on the parameters of the potential in this model are valid only in a range of initial conditions. A reasonably sharp boundary lies where the initial energy density of the scalar field is equal to that of the background radiation component. An initial quintessence energy density above this equipartition value lead to a solution that will not have joined the tracker solution by the present epoch. These non-tracker solutions possess the property that their present equation of state is very compatible with the observed bounds and independent of the exponent of the potential.

The spallative production rates of Lithium, Beryllium and Boron (LiBeB) are a necessary component in any calculation of the evolution of these nuclei in the Galaxy. Previous calculations of these rates relied on two assumptions relating to the nuclear physics aspects: the straight-ahead approximation that describes the distribution of fragment energies and the assumption that the major contributor to the production rate arises from single-step reactions between primary cosmic ray projectiles and interstellar medium targets. We examine both assumptions by using a semi-empirical description for the spall's energy distribution and by including the reactions that proceed via intermediary fragments. After relaxing the straight-ahead approximation we find the changes in the production rates and emerging fluxes are small and do not warrant rejection of this approximation. In contrast we discover that two-step reactions can alter the production rate considerably leading to noticeable increases in the efficiency of producing the LiBeB nuclei. Motivated by this result we introduce a cascade technique to compute the production rates exactly and find that the results differ only slightly from those of our two-step calculations. We thus conclude that terminating the reaction network at the two-step order is sufficiently accurate for current studies of spallation.

Current observational data favor cosmological models which differ from the standard model due to the presence of some form of dark energy and, perhaps, by additional contributions to the more familiar dark matter. Primordial nucleosynthesis provides a window on the very early evolution of the universe and constraints from Big Bang Nucleosynthesis (BBN) can bound the parameters of models for dark matter/energy at redshifts of order ten billion. The spectrum of temperature fluctuations imprinted on the Cosmic Microwave Background (CMB) radiation opens a completely different window on the universe at epochs from redshifts of order ten thousand to nearly the present. The CMB anisotropy spectrum provides constraints on new physics which are independent of, and complementary to those from BBN. Here we consider three classes of models for the dark matter/energy: extra particles which were relativistic during the early evolution of the universe (X); Quintessence models involving a minimally-coupled scalar field (Q); models with a non-minimally coupled scalar field which modify the strength of gravity during the early evolution of the universe (G). We constrain the parameters of these models using data from BBN and the CMB and identify the allowed regions in their parameter spaces consistent with the more demanding joint BBN and CMB constraints. For X and Q such consistency is relatively easy to find; it is more difficult for the G models with an inverse power law potential for the scalar field.

Recent cosmic microwave background (CMB) results from BOOMERANG, MAXIMA, and DASI provide cosmological constraints on new physics that can be competitive with those derived from Big Bang Nucleosynthesis (BBN). In particular, both CMB and BBN can be used to place limits on models involving neutrino degeneracy and additional relativistic degrees of freedom. However, for the case of the CMB, these constraints are, in general, sensitive to the assumed priors. We examine the CMB and BBN constraints on such models and study the sensitivity of ``new physics" to the assumed priors. If we add a constraint on the age of the universe ($t_0 \ga 11$ Gyr), then for models with a cosmological constant, the range of baryon densities and neutrino degeneracy parameters allowed by the CMB and BBN is fairly robust: $\eta_{10} = 6.0 \pm 0.6$, \deln $\la 6$, $\xi_e \la 0.3$. In the absence of new physics, models without a cosmological constant are only marginally compatible with recent CMB observations (excluded at the 93% confidence level).

Primordial nucleosynthesis provides a probe of the universal abundance of baryons when the universe was only a few minutes old. Recent observations of anisotropy in the cosmic microwave background (CMB) probe the baryon abundance when the universe was several hundred thousand years old. Observations of type Ia supernovae and clusters of galaxies in the very recent past, when the universe is several billion years old and older, provide a complementary measure of the baryon density in excellent agreement with the early universe values. The general agreement among the three measurements represents an impressive confirmation of the standard model of cosmology. However, there is a hint that the CMB observations may not be in perfect agreement with those from big bang nucleosynthesis (BBN). If this "tension" between BBN and the CMB persists, the standard model of cosmology may need to be modified. Here, in a contribution dedicated to Silvia Torres-Peimbert and Manuel Peimbert, we describe how an asymmetry between neutrinos and antineutrinos ("neutrino degeneracy") has the potential for resolving this possible conflict between BBN and the CMB.

We calculate the production rates of Li6, Li7, Be9, B10 and B11 via spallation of Carbon, Nitrogen and Oxygen nuclei by protons and alpha-particles and by alpha-alpha fusion reactions. We include recent measurements of the cross sections of alpha-alpha fusion reactions and find that the computations yield rates of Li6 and Li7 production that are nearly a factor of two smaller than previously calculated. We begin by using the `straight ahead' approximation for the fragment energy and the `leaky-box' model for product capture in the Galaxy. In addition we test the straight ahead approximation by recalculating the production rates using an empirical description of the fragment energy distribution and find that the results closely match. We have also calculated the rates for various cosmic ray spectra and find that the hardest spectra tested decrease the rates with CR CNO by approximately an order of magnitude relative to our chosen standard. Finally we have computed the Population I elemental ratios and the Population II scaling relations for our standard and find that our computations predict an abundance of Lithium for a given abundance of Beryllium that is 1/4 smaller than previously derived.