We apply for the first time fission yields determined across the chart of nuclides from the macroscopic-microscopic theory of the Finite Range Liquid Drop Model to simulations of rapid neutron capture ($r$-process) nucleosynthesis. With the fission rates and yields derived within the same theoretical framework utilized for other relevant nuclear data, our results represent an important step toward self-consistent applications of macroscopic-microscopic models in $r$-process calculations. The yields from this model are wide for nuclei with extreme neutron excess. We show that these wide distributions of neutron-rich nuclei, and particularly the asymmetric yields for key species that fission at late times in the $r$ process, can contribute significantly to the abundances of the lighter heavy elements, specifically the light precious metals palladium and silver. Since these asymmetric yields correspondingly also deposit into the lanthanide region, we consider the possible evidence for co-production by comparing our nucleosynthesis results directly with the trends in the elemental ratios of metal-poor stars rich in $r$-process material. We show that for $r$-process enhanced stars palladium over europium and silver over europium display mostly flat trends suggestive of co-production and compare to the lanthanum over europium trend which is often used to justify robustness arguments in the lanthanide region. We find that such robustness arguments may be extendable down to palladium and heavier and demonstrate that fission deposition is a mechanism by which such a universality or robustness can be achieved.

Lanthanide element signatures are key to understanding many astrophysical observables from merger kilonova light curves to stellar and solar abundances. To learn about the lanthanide element synthesis which enriched our solar system, we apply the statistical method of Markov Chain Monte Carlo to examine the nuclear masses capable of forming the $r$-process rare-earth abundance peak. We describe the physical constraints we implement with this statistical approach and demonstrate the use of the parallel chains method to explore the multi-dimensional parameter space. We apply our procedure to three moderately neutron-rich astrophysical outflows with distinct types of $r$-process dynamics. We show that the mass solutions found are dependent on outflow conditions and are related to the $r$-process path. We describe in detail the mechanism behind peak formation in each case. We then compare our mass predictions for neutron-rich neodymium and samarium isotopes to the latest experimental data from the CPT at CARIBU. We find our mass predictions given outflows which undergo an extended (n,$\gamma$)$\rightleftarrows$($\gamma$,n) equilibrium to be those most compatible with both observational solar abundances and neutron-rich mass measurements.

Of the variations in the elemental abundance patterns of stars enhanced with $r$-process elements, the variation in the relative actinide-to-lanthanide ratio is among the most significant. We investigate the source of these actinide differences in order to determine whether these variations are due to natural differences in astrophysical sites, or due to the uncertain nuclear properties that are accessed in $r$-process sites. We find that variations between relative stellar actinide abundances is most likely astrophysical in nature, owing to how neutron-rich the ejecta from an $r$-process event may be. Furthermore, if an $r$-process site is capable of generating variations in the neutron-richness of its ejected material, then only one type of $r$-process site is needed to explain all levels of observed relative actinide enhancements.

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.

The astrophysical production site of the heaviest elements in the universe remains a mystery. Incorporating heavy element signatures of metal-poor, $r$-process enhanced stars into theoretical studies of $r$-process production can offer crucial constraints on the origin of heavy elements. In this study, we introduce and apply the "Actinide-Dilution with Matching" model to a variety of stellar groups ranging from actinide-deficient to actinide-enhanced to empirically characterize $r$-process ejecta mass as a function of electron fraction. We find that actinide-boost stars do not indicate the need for a unique and separate $r$-process progenitor. Rather, small variations of neutron richness within the same type of $r$-process event can account for all observed levels of actinide enhancements. The very low-$Y_e$, fission-cycling ejecta of an $r$-process event need only constitute 10-30% of the total ejecta mass to accommodate most actinide abundances of metal-poor stars. We find that our empirical $Y_e$ distributions of ejecta are similar to those inferred from studies of GW170817 mass ejecta ratios, which is consistent with neutron-star mergers being a source of the heavy elements in metal-poor, $r$-process enhanced stars.

In July 2018 an FRIB Theory Alliance program was held on the implications of GW170817 and its associated kilonova for r-process nucleosynthesis. Topics of discussion included the astrophysical and nuclear physics uncertainties in the interpretation of the GW170817 kilonova, what we can learn about the astrophysical site or sites of the r process from this event, and the advances in nuclear experiment and theory most crucial to pursue in light of the new data. Here we compile a selection of scientific contributions to the workshop, broadly representative of progress in r-process studies since the GW170817 event.

We study the cosmic MeV neutrino background from accretion disks formed during collapsars and the coalescence of compact-object mergers. We provide updated estimates, including detection rates, of relic neutrinos from collapsars, as well as estimates for neutrinos that are produced in mergers. Our results show that diffuse neutrinos detected at HyperK would likely include some that were emitted from binary neutron-star mergers. The collapsar rate is uncertain, but at its upper limit relic neutrinos from these sources would provide a significant contribution to the Cosmic Diffuse Neutrino Background.

We present a covariant ray tracing algorithm for computing high-resolution neutrino distributions in general relativistic numerical spacetimes with hydrodynamical sources. Our formulation treats the very important effect of elastic scattering of neutrinos off of nuclei and nucleons (changing the neutrino's direction but not energy) by incorporating estimates of the background neutrino fields. Background fields provide information about the spectra and intensities of the neutrinos scattered into each ray. These background fields may be taken from a low-order moment simulation or be ignored, in which case the method reduces to a standard state-of-the-art ray tracing formulation. The method handles radiation in regimes spanning optically thick to optically thin. We test the new code, highlight its strengths and weaknesses, and apply it to a simulation of a neutron star merger to compute neutrino fluxes and spectra, and to demonstrate a neutrino flavor oscillation calculation. In that environment, we find qualitatively different fluxes, spectra, and oscillation behaviors when elastic scattering is included.

Neutron star mergers offer unique conditions for the creation of the heavy elements and additionally provide a testbed for our understanding of this synthesis known as the $r$-process. We have performed dynamical nucleosynthesis calculations and identified a single isotope, $^{254}$Cf, which has a particularly high impact on the brightness of electromagnetic transients associated with mergers on the order of 15 to 250 days. This is due to the anomalously long half-life of this isotope and the efficiency of fission thermalization compared to other nuclear channels. We estimate the fission fragment yield of this nucleus and outline the astrophysical conditions under which $^{254}$Cf has the greatest impact to the light curve. Future observations in the middle-IR which are bright during this regime could indicate the production of actinide nucleosynthesis.

This is an exciting time for the study of r-process nucleosynthesis. Recently, a neutron star merger GW170817 was observed in extraordinary detail with gravitational waves and electromagnetic radiation from radio to gamma rays. The very red color of the associated kilonova suggests that neutron star mergers are an important r-process site. Astrophysical simulations of neutron star mergers and core collapse supernovae are making rapid progress. Detection of both, electron neutrinos and antineutrinos from the next galactic supernova will constrain the composition of neutrino-driven winds and provide unique nucleosynthesis information. Finally FRIB and other rare-isotope beam facilities will soon have dramatic new capabilities to synthesize many neutron-rich nuclei that are involved in the r-process. The new capabilities can significantly improve our understanding of the r-process and likely resolve one of the main outstanding problems in classical nuclear astrophysics. However, to make best use of the new experimental capabilities and to fully interpret the results, a great deal of infrastructure is needed in many related areas of astrophysics, astronomy, and nuclear theory. We will place these experiments in context by discussing astrophysical simulations and observations of r-process sites, observations of stellar abundances, galactic chemical evolution, and nuclear theory for the structure and reactions of very neutron-rich nuclei. This review paper was initiated at a three-week International Collaborations in Nuclear Theory program in June 2016 where we explored promising r-process experiments and discussed their likely impact, and their astrophysical, astronomical, and nuclear theory context.

Simulations of neutrino flavor evolution in compact merger environments have shown that neutrino flavor, and hence nucleosynthesis, can be strongly affected by the presence of matter-neutrino resonances (MNRs), where there is a cancelation between the matter and the neutrino potential. Simulations performed thus far follow flavor evolution along a single neutrino trajectory, but self-consistency requires all trajectories to be treated simultaneously, and it has not been known whether MNR phenomena would still occur in multi-angle models. In this paper, we present the first fully multi-angle calculations of MNR. We find that familiar MNR phenomena, where neutrinos transform to a greater extent than anti-neutrinos and a feedback mechanism maintains the cancellation between the matter and neutrino potential, still occurs for a subset of angular bins, although the flavor transformation is not as efficient as in the single-angle case. In addition, we find other types of flavor transformation that are not seen in single-angle simulations. These flavor transformation phenomena appear to be robust and are present for a wide range of model parameters, as long as an MNR is present. Although computational constraints currently limit us to models with spherical symmetry, our results suggest that the presence of an MNR generally leads to large-scale neutrino flavor evolution in multi-angle systems.

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 bulk of the rare earth elements are believed to be synthesized in the rapid neutron capture process or $r$ process of nucleosynthesis. The solar $r$-process residuals show a small peak in the rare earths around $A\sim 160$, which is proposed to be formed dynamically during the end phase of the $r$ process by a pileup of material. This abundance feature is of particular importance as it is sensitive to both the nuclear physics inputs and the astrophysical conditions of the main $r$ process. We explore the formation of the rare earth peak from the perspective of an inverse problem, using Monte Carlo studies of nuclear masses to investigate the unknown nuclear properties required to best match rare earth abundance sector of the solar isotopic residuals. When nuclear masses are changed, we recalculate the relevant $\beta$-decay properties and neutron capture rates in the rare earth region. The feedback provided by this observational constraint allows for the reverse engineering of nuclear properties far from stability where no experimental information exists. We investigate a range of astrophysical conditions with this method and show how these lead to different predictions in the nuclear properties influential to the formation of the rare earth peak. We conclude that targeted experimental campaigns in this region will help to resolve the type of conditions responsible for the production of the rare earth nuclei, and will provide new insights into the longstanding problem of the astrophysical site(s) of the $r$ process.

The Matter-Neutrino Resonance (MNR) phenomenon has the potential to significantly alter the flavor content of neutrinos emitted from compact object mergers. We present the first calculations of MNR transitions using neutrino self interaction potentials and matter potentials generated selfconsistently from a dynamical model of a three-dimensional neutron star merger. In the context of the single angle approximation, we find that Symmetric and Standard MNR transitions occur in both normal and inverted hierarchy scenarios. We examine the spatial regions above the merger remnant where propagating neutrinos will encounter the matter neutrino resonance and find that a significant fraction of the neutrinos are likely to undergo MNR transitions.

This white paper informs the nuclear astrophysics community and funding agencies about the scientific directions and priorities of the field and provides input from this community for the 2015 Nuclear Science Long Range Plan. It summarizes the outcome of the nuclear astrophysics town meeting that was held on August 21-23, 2014 in College Station at the campus of Texas A&M University in preparation of the NSAC Nuclear Science Long Range Plan. It also reflects the outcome of an earlier town meeting of the nuclear astrophysics community organized by the Joint Institute for Nuclear Astrophysics (JINA) on October 9- 10, 2012 Detroit, Michigan, with the purpose of developing a vision for nuclear astrophysics in light of the recent NRC decadal surveys in nuclear physics (NP2010) and astronomy (ASTRO2010). The white paper is furthermore informed by the town meeting of the Association of Research at University Nuclear Accelerators (ARUNA) that took place at the University of Notre Dame on June 12-13, 2014. In summary we find that nuclear astrophysics is a modern and vibrant field addressing fundamental science questions at the intersection of nuclear physics and astrophysics. These questions relate to the origin of the elements, the nuclear engines that drive life and death of stars, and the properties of dense matter. A broad range of nuclear accelerator facilities, astronomical observatories, theory efforts, and computational capabilities are needed. With the developments outlined in this white paper, answers to long standing key questions are well within reach in the coming decade.

The primary astrophysical source of the rare earth elements is the rapid neutron capture process ($r$ process). The rare earth peak that is seen in the solar $r$-process residuals has been proposed to originate as a pile-up of nuclei during the end of the $r$ process. We introduce a new method utilizing Monte Carlo studies of nuclear masses in the rare earth region, that includes self-consistently adjusting $\beta$-decay rates and neutron capture rates, to find the mass surfaces necessary for the formation of the rare earth peak. We demonstrate our method with two types of astrophysical scenarios, one corresponding conditions typical of core-collapse supernova winds and one corresponding to conditions typical of the ejection of the material from the tidal tails of neutron star mergers. In each type of astrophysical conditions, this method successfully locates a region of enhanced stability in the mass surface that is responsible for the rare earth peak. For each scenario, we find that the change in the mass surface has qualitatively different features, thus future measurements can shed light on the type of environment in which the $r$ process occurred.

Neutrinos are copiously emitted from black hole accretion disks playing a fundamental role in their evolution, as well as in the production of gamma ray bursts and r-process nucleosynthesis. The black hole generates a strong gravitational field able to change the properties of the emerging neutrinos. We study the influence of the black hole spin on the structure of the neutrino surfaces, neutrino luminosities, average neutrino energies, and event counts at SuperK. We consider several disk models and provide estimates that cover different black hole efficiency scenarios. We discuss the influence of the detector's inclination with respect to the axis of the torus on neutrino properties. We find that tori around spinning black holes have larger luminosities, energies and rates compared to tori around static black holes, and that the inclination of the observer causes a reduction in the luminosities and detection rates but an increase in the average energies.

Matter Neutrino Resonances (MNRs) can drastically modify neutrino flavor evolution in astrophysical environments and may significantly impact nucleosynthesis. Here we further investigate the underlying physics of MNR type flavor transitions. We provide generalized resonance conditions and make analytical predictions for the behavior of the system. We discuss the adiabatic evolution of these transitions, considering both Symmetric and Standard scenarios. Symmetric MNR transitions differ from Standard MNR transitions in that both neutrinos and antineutrinos can completely transform to other flavors simultaneously. We provide an example of the simplest system in which such transitions can occur with a neutrino and an antineutrino having a single energy and emission angle. We further apply linearized stability analysis to predict the location of self-induced nutation type (or bipolar) oscillations due to neutrino-neutrino interactions in the regions where MNR is ineffective. In all cases, we compare our analytical predictions to numerical calculations.

The astrophysical rapid neutron capture process or `$r$ process' of nucleosynthesis is believed to be responsible for the production of approximately half the heavy element abundances found in nature. This multifaceted problem remains one of the greatest open challenges in all of physics. Knowledge of nuclear physics properties such as masses, $\beta$-decay and neutron capture rates, as well as $\beta$-delayed neutron emission probabilities are critical inputs that go into calculations of $r$-process nucleosynthesis. While properties of nuclei near stability have been established, much still remains unknown regarding neutron-rich nuclei far from stability that may participate in the $r$ process. Sensitivity studies gauge the astrophysical response of a change in nuclear physics input(s) which allows for the isolation of the most important nuclear properties that shape the final abundances observed in nature. This review summarizes the extent of recent sensitivity studies and highlights how these studies play a key role in facilitating new insight into the $r$ process. The development of these tools promotes a focused effort for state-of-the-art measurements, motivates construction of new facilities and will ultimately move the community towards addressing the grand challenge of `How were the elements from iron to uranium made?'.

Matter-neutrino resonances (MNR) can occur in environments where the flux of electron antineutrinos is greater than the flux of electron neutrinos. These resonances may result in dramatic neutrino flavor transformation. Compact object merger disks are an example of an environment where electron antineutrinos outnumber neutrinos. We study MNR resonances in several such disk configurations and find two qualitatively different types of matter-neutrino resonances: a standard MNR and a symmetric MNR. We examine the transformation that occurs in each type of resonance and explore the consequences for nucleosynthesis.

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.

The synthesis of heavy elements in the Universe presents several challenges. From one side the astrophysical site is still undetermined and on other hand the input from nuclear physics requires the knowledge of properties of exotic nuclei, some of them perhaps accessible in ion beam facilities. Black hole accretion disks have been proposed as possible r-process sites. Analogously to Supernovae these objects emit huge amounts of neutrinos. We discuss the neutrino emission from black hole accretion disks. In particular we show the influence that the black hole strong gravitational field has on changing the electron fraction relevant to the synthesis of elements.

Accretion disks arising from neutron star- neutron star mergers or black hole- neutron star mergers produce large numbers of neutrinos and antineutrinos. In contrast to other astrophysical scenarios, like supernovae, in mergers the antineutrinos outnumber the neutrinos. This antineutrino dominance gives neutrinos from merger disks the opportunity to exhibit new oscillation physics, specifically a matter-neutrino resonance. We explore this resonance, finding that consequences can be a large transition of $\nu_e$ to other flavors, while the $\bar{\nu}_e$s return to their initial state. We present numerical calculations of neutrinos from merger disks and compare with a single energy model. We explain both the basic features and the conditions for a transition.

The likely outcome of a compact object merger event is a central black hole surrounded by a rapidly accreting torus of debris. This disk of debris is a rich source of element synthesis, the outcome of which is needed to predict electromagnetic counterparts of individual events and to understand the contribution of mergers to galactic chemical evolution. Here we study disk outflow nucleosynthesis in the context of a two-dimensional, time-dependent black hole-neutron star merger accretion disk model. We use two time snapshots from this model to examine the impact of the evolution of the neutrino fluxes from the disk on the element synthesis. While the neutrino fluxes from the early-time disk snapshot appear to favor neutron-rich outflows, by the late-time snapshot the situation is reversed. As a result we find copious production of Nickel-56 in the outflows.

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.

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 study the sensitivity of the r-process abundance pattern to neutron capture rates along the rare earth region (A~150 to A~180). We introduce the concepts of large nuclear flow and flow saturation which determine the neutron capture rates that are influential in setting the rare earth abundances. We illustrate the value of the two concepts by considering high entropy conditions favorable for rare earth peak production and identifying important neutron capture rates among the rare earth isotopes. We also show how these rates influence nuclear flow and specific sections of the abundance pattern.

The astrophysical site or sites responsible for the r-process of nucleosynthesis still remains an enigma. Since the rare earth region is formed in the latter stages of the r-process it provides a unique probe of the astrophysical conditions during which the r-process takes place. We use features of a successful rare earth region in the context of a high entropy r-process (S>100k_B) and discuss the types of astrophysical conditions that produce abundance patterns that best match meteoritic and observational data. Despite uncertainties in nuclear physics input, this method effectively constrains astrophysical conditions.

We examine the prospects for producing Nickel-56 from black hole accretion disks, by examining a range of steady state disk models. We focus on relatively slowly accreting disks in the range of 0.05 - 1 solar masses per second, as are thought to be appropriate for the central engines of long-duration gamma-ray bursts. We find that significant amounts of Nickel-56 are produced over a wide range of parameter space. We discuss the influence of entropy, outflow timescale and initial disk position on mass fraction of Nickel-56 which is produced. We keep careful track of the weak interactions to ensure reliable calculations of the electron fraction, and discuss the role of the neutrinos.

We study the formation and final structure of the rare earth peak ($A\sim160$) of the $r$-process nucleosynthesis. The rare earth peak forms at late times in the $r$-process after neutron exhaustion (neutron-to-seed ratio unity or R=1) as matter decays back to stability. Since rare earth peak formation does not occur during \nggn \ equilibrium it is sensitive to the strong interplay between late time thermodynamic evolution and nuclear physics input. Depending on the conditions the peak forms either because of the pattern of the neutron capture rates or because of the pattern of the separation energies. We analyze three mass models under different thermodynamic conditions. We find that the subtleties of each mass model, including separation energies and neutron capture rates, influence not only the final shape of the peak but also when it forms. We identify the range of nuclei which are influential in rare earth peak formation.

Black hole accretion disks have been proposed as good candidates for a range of interesting nucleosynthesis, including the r-process. The presence of the black hole influences the neutrino fluxes and affects the nucleosynthesis resulting from the interaction of the emitted neutrinos and hot outflowing material ejected from the disk. We study the impact of general relativistic effects on the neutrinos emitted from black hole accretion disks. We present abundances obtained by considering null geodesics and energy shifts for two different disk models. We find that both the bending of the neutrino trajectories and the energy shifts have important consequences for the nucleosynthetic outcome

Recently, it has been demonstrated that neutrinos in a supernova oscillate collectively. This process occurs much deeper than the conventional matter-induced MSW effect and hence may have an impact on nucleosynthesis. In this paper we explore the effects of collective neutrino oscillations on the r-process, using representative late-time neutrino spectra and outflow models. We find that accurate modeling of the collective oscillations is essential for this analysis. As an illustration, the often-used "single-angle" approximation makes grossly inaccurate predictions for the yields in our setup. With the proper multiangle treatment, the effect of the oscillations is found to be less dramatic, but still significant. Since the oscillation patterns are sensitive to the details of the emitted fluxes and the sign of the neutrino mass hierarchy, so are the r-process yields. The magnitude of the effect also depends sensitively on the astrophysical conditions - in particular on the interplay between the time when nuclei begin to exist in significant numbers and the time when the collective oscillation begins. A more definitive understanding of the astrophysical conditions, and accurate modeling of the collective oscillations for those conditions, is necessary.

We study proton rich nucleosynthesis in windlike outflows from gamma-ray bursts accretion disks with the aim to determine if such outflows are a site of the vp-process. The efficacy of this vp-process depends on thermodynamic and hydrodynamic factors. We discuss the importance of the entropy of the material, the outflow rate, the initial ejection point and accretion rate of the disk. In some cases the vp-process pushes the nucleosynthesis out to A~100 and produces light p-nuclei. However, even when these nuclei are not produced, neutrino induced interactions can significantly alter the abundance pattern and cannot be neglected.

While it is well known that neutrinos are emitted from standard core collapse protoneutron star supernovae, less attention has been focused on neutrinos from accretion disks. These disks occur in some supernovae (i.e. "collapsars") as well as in compact object mergers, and they emit neutrinos with similar properties to those from protoneutron star supernovae. These disks and their neutrinos play an important role in our understanding of gamma ray bursts as well as the nucleosynthesis they produce. We study a disk that forms in the merger of a black hole and a neutron star and examine the neutrino fluxes, luminosities and neutrino surfaces for the disk. We also estimate the number of events that would be registered in current and proposed supernova neutrino detectors if such an event were to occur in the Galaxy.

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.

We consider hot accretion disk outflows from black hole - neutron star mergers in the context of the nucleosynthesis they produce. We begin with a three dimensional numerical model of a black hole - neutron star merger and calculate the neutrino and antineutrino fluxes emitted from the resulting accretion disk. We then follow the element synthesis in material outflowing the disk along parameterized trajectories. We find that at least a weak r-process is produced, and in some cases a main r-process as well. The neutron-rich conditions required for this production of r-process nuclei stem directly from the interactions of the neutrinos emitted by the disk with the free neutrons and protons in the outflow.

We investigate the impact of neutron capture rates near the A=130 peak on the $r$-process abundance pattern. We show that these capture rates can alter the abundances of individual nuclear species, not only in the region of A=130 peak, but also throughout the abundance pattern. We discuss the nonequilibrium processes that produce these abundance changes and determine which capture rates have the most significant impact.

Recently, we suggested that low-energy beta-beam neutrinos can be very useful for the study of supernova neutrino interactions. In this paper, we examine the use of a such experiment for the analysis of a supernova neutrino signal. Since supernova neutrinos are oscillating, it is very likely that the terrestrial spectrum of supernova neutrinos of a given flavor will not be the same as the energy distribution with which these neutrinos were first emitted. We demonstrate the efficacy of the proposed method for untangling multiple neutrino spectra. This is an essential feature of any model aiming at gaining information about the supernova mechanism, probing proto-neutron star physics, and understanding supernova nucleosynthesis, such as the neutrino process and the r-process. We also consider the efficacy of different experimental approaches including measurements at multiple beam energies and detector configurations.

Recent halo star abundance observations exhibit an important feature of consequence to the r-process: the presence of a main r-process between the second and third peaks which is consistent among halo stars. We explore fission cycling and steady-beta flow as the driving mechanisms behind this feature. The presence of fission cycling during the r-process can account for nucleosynthesis yields between the second and third peaks, whereas the presence of steady-beta flow can account for consistent r-process patterns, robust under small variations in astrophysical conditions. We employ the neutrino-driven wind of the core-collapse supernova to examine fission cycling and steady-beta flow in the r-process. As the traditional neutrino-driven wind model does not produce the required very neutron-rich conditions for these mechanisms, we examine changes to the neutrino physics necessary for fission cycling to occur in the neutrino-driven wind environment, and we explore under what conditions steady-beta flow is obtained.

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.

It has long been suggested that fission cycling may play an important role in the r-process. Fission cycling can only occur in a very neutron rich environment. In traditional calculations of the neutrino driven wind of the core-collapse supernova, the environment is not sufficiently neutron rich to produce the r-process elements. However, we show that with a reduction of the electron neutrino flux coming from the supernova, fission cycling does occur and furthermore it produces an abundance pattern which is consistent with observed r-process abundance pattern in halo stars. Such a reduction can be caused by active-sterile neutrino oscillations or other new physics.

Neutrinos from core collapse supernovae can be emitted from a rapidly accreting disk surrounding a black hole, instead of the canonical proto-neutron star. For Galactic events, detector count rates are considerable and in fact can be in the thousands for Super-Kamiokande. The rate of occurrence of these accreting disks in the Galaxy is predicted to be on the order of 10^-5 yr^-1, yet there is little observational evidence to provide an upper limit on their formation rate. It would therefore be useful to discriminate between neutrinos which have been produced in a proto-neutron star and those which have been produced accretion disks. In order to distinguish between the two scenarios, either the time profile of the neutrino luminosity or the relative fluxes of different neutrino flavors may be considered. There are some signals that would clearly point to one scenario or the other.

Only weakly interacting, neutrinos are the principal messengers reaching us from the center of a supernova. Terrestrial neutrino telescopes, such as SNO and SuperKamiokande, can provide precious information about the processes in the core of a collapsing and exploding star. But the information about the supernova that a neutrino detector can supply, is restricted by the fact that little experimental data on the neutrino-nucleus cross sections exists and by the uncertainties in theoretical calculations. In this letter, we propose a novel procedure that determines the response of a target nucleus in a supernova-neutrino detector directly, by using low-energy beta-beams. We show that fitting 'synthetic' spectra, constructed by taking linear combinations of beta-beam spectra, to the original supernova-neutrino spectra reproduces the folded differential cross sections very accurately. Comparing the response in a terrestrial detector to these synthetic responses provides a direct way to determine the main parameters of the supernova-neutrino energy-distribution.

We investigate nucleosynthesis in the supernovae post-core bounce neutrino-driven wind environment in the presence of active-sterile neutrino transformation. We consider active-sterile neutrino oscillations for a range of mixing parameters: vacuum mass-squared differences of 0.1 eV^2 < dm^2 < 100 eV^2, and vacuum mixing angles of sin^2(2 theta_v) > 10^-4. We find a consistent r-process pattern for a large range of mixing parameters that is in rough agreement with the halo star CS 22892-052 abundances and the pattern shape is determined by fission cycling. We find that the allowed region for the formation of the r-process peaks overlaps the LSND and NSBL (3+1) allowed region.

We examine the nucleosynthesis products that are produced in the outflow from rapidly accreting disks. We find that the type of element synthesis varies dramatically with the degree of neutrino trapping in the disk and therefore the accretion rate of the disk. Disks with relatively high accretion rates such as 10 M_solar/s can produce very neutron rich nuclei that are found in the r process. Disks with more moderate accretion rates can produce copious amounts of Nickel as well as the light elements such as Lithium and Boron. Disks with lower accretion rates such as 0.1 M_solar/s produce large amounts of Nickel as well as some unusual nuclei such as Ti-49, Sc-45, Zn-64, and Mo-92. This wide array of potential nucleosynthesis products is due to the varying influence of electron neutrinos and antineutrinos emitted from the disk on the neutron-to-proton ratio in the outflow. We use a parameterization for the outflow and discuss our results in terms of entropy and outflow acceleration.

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.

We discuss the possibility that r-process nucleosynthesis may occur in the winds from gamma ray burst accretion disks. This can happen if the temperature of the disk is sufficiently high that electron antineutrinos are trapped as well as neutrinos. This implies accretion disks with greater than a solar mass per second accretion rate, although lower accretion rates with higher black hole spin parameters may provide viable environments as well. Additionally, the outflow from the disk must either have relatively low entropy, e.g. around s = 10, or the initial acceleration of the wind must be slow enough that it is neutrino and antineutrino capture as opposed to electron and positron capture that sets the electron fraction.

We examine the composition of matter as it flows away from gamma ray burst accretion disks, in order to determine what sort of nucleosynthesis may occur. Since there is a large flux of neutrinos leaving the surface of the disk, the electron fraction of the outflowing material will change due to charged current neutrino interactions. We calculate the electron fraction in the wind using detailed neutrino fluxes from every point on the disk and study a range of trajectories and outflow conditions for several different accretion disk models. We find that low electron fractions, conducive to making r-process elements, only appear in outflows from disks with high accretion rates that have a significant region both of trapped neutrinos and antineutrinos. Disks with lower accretion rates that have only a significant region of trapped neutrinos can have outflows with very high electron fractions, whereas the lowest accretion rate disks with little trapping have outflow electrons fractions of closer to one half.

Two new facilities have recently been proposed to measure low energy neutrino-nucleus cross sections, the nu-SNS (Spallation Neutron Source) and low energy beta beams. The former produces neutrinos by pion decay at rest, while the latter produces neutrinos from the beta decays of accelerated ions. One of the uses of neutrino-nucleus cross section measurements is for supernova studies, where typical neutrino energies are 10s of MeV. In this energy range there are many different components to the nuclear response and this makes the theoretical interpretation of the results of such an experiment complex. Although even one measurement on a heavy nucleus such as lead is much anticipated, more than one data set would be still better. We suggest that this can be done by breaking the electron spectrum down into the parts produced in coincidence with one or two neutrons, running a beta beam at more than one energy, comparing the spectra produced with pions and a beta beam or any combination of these.

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 calculate the nuclear composition of matter in accretion disks surrounding stellar mass black holes as are thought to accompany gamma-ray bursts (GRBs). We follow a mass element in the accretion disk starting at the point of nuclear dissociation and calculate the evolution of the electron fraction due to electron, positron, electron neutrino and electron antineutrino captures. We find that the neutronization of the disk material by electron capture can be reversed by neutrino interactions in the inner regions of disks with accretion rates of 1 M_solar/s and higher. For these cases the inner disk regions are optically thick to neutrinos, and so to estimate the emitted neutrino fluxes we find the surface of last scattering for the neutrinos (the equivalent of the proto-neutron star neutrinosphere) for each optically thick disk model. We also estimate the influence of neutrino interactions on the neutron-to-proton ratio in outflows from GRB accretion disks, and find it can be significant even when the disk is optically thin to neutrinos.

Light curves from a growing number of Gamma Ray Bursts (GRBs) indicate that GRBs copiously produce radioactive Ni moving outward at fractions of the speed of light. We calculate nuclear abundances of elements accompanying the outflowing Ni under the assumption that this Ni originates from a wind blown off of a viscous accretion disk. We also show that GRB's likely contribute appreciably to the galactic inventory of 42Ca, 45Sc, 46Ti, 49Ti, 63Cu, and may be a principal site for the production of 64Zn.

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 examine the prospects for using lead as a supernova-neutrino detector by considering the spectrum of electrons produced, and the number of one- and two-neutron events. We show that the electron energy spectrum from charged-current reactions can be used to extract information about the high-temperature component of the neutrino spectrum. Some degree of electron neutrino oscillation is expected in the supernova envelope. We examine the prospects for untangling the signatures of various oscillation scenarios, including, e.g. normal or inverted hierarchies, and different values for the small mixing angle, theta_13.

Recently observed emission lines in the X-ray afterglow of gamma ray bursts suggest that iron group elements are either produced in the gamma ray burst, or are present nearby. If this material is the product of a thermonuclear burn, then such material would be expected to be rich in Nickel-56. If the nickel remains partially ionized, this prevents the electron capture reaction normally associated with the decay of Nickel-56, dramatically increasing the decay timescale. Here we examine the consequences of rapid ejection of a fraction of a solar mass of iron group material from the center of a collapsar/hypernova. The exact rate of decay then depends on the details of the ionization and therefore the ejection process. Future observations of iron, nickel and cobalt lines can be used to diagnose the origin of these elements and to better understand the astrophysical site of gamma ray bursts. In this model, the X-ray lines of these iron-group elements could be detected in suspected hypernovae that did not produce an observable gamma ray burst due to beaming.

We examine active-sterile neutrino conversion in the late time post-core-bounce supernova environment. By including the effect of feedback on the Mikheyev-Smirnov-Wolfenstein (MSW) conversion potential, we obtain a large range of neutrino mixing parameters which produce a favorable environment for the r-process. We look at the signature of this effect in the current generation of neutrino detectors now coming on line. We also investigate the impact of the neutrino-neutrino forward scattering-induced potential on the MSW conversion.

In the hypernova/collapsar model of gamma-ray bursts, it is natural that radiation is emitted by the inner engine for some time after the burst. This has been discussed as a possible source of the X-ray line emission observed in some afterglows. We show here that the natural geometry of a hypernova --a source of radiation at the bottom of a deep funnel-- has very significant consequences for the shape and central energy of the observed emission lines: the lines acquire a very broad scattering wing on the low-energy side and a characteristic second peak one Compton wavelength away from the initial energy, due to the once- and twice-scattered photons. We suggest that this explains the large width of the observed lines. Furthermore, the downscattering lowers the central line energy (by up to 1 keV in the source rest frame, before the lines become unrecognizable), so that the observed line energies can become consistent with originating from cobalt and nickel, as expected in few-day-old supernova ejecta.

Theories involving extra dimensions, a low (TeV) string scale and bulk singlet neutrinos will produce an effective neutrino magnetic moment which may be large (< 10^{-11} mu_B). The effective magnetic moment increases with neutrino energy, and therefore high energy reactions are most useful for limiting the allowed number of extra dimensions. We examine constraints from both neutrino-electron scattering and also astrophysical environments. We find that supernova energy loss considerations require a number of extra dimensions, n > 1, for an electron neutrino-bulk neutrino Yukawa coupling of order 1.

We study the impact of Standard Model singlet neutrinos on neutrino flavor transformation. We focus on an extension of the Zee model which includes singlet neutrinos, and find that the best limits on the interactions of the singlet neutrinos come from astrophysical phenomena. Singlet neutrino-electron scattering will impact both the mattter enhanced flavor transformation potential as well as detector cross sections. If electron neutrino - singlet neutrino oscillations are responsible for the solar neutrino anomaly, then the limit on the singlet neutrino interaction strength is of order of the weak interaction scale. Zee model modification of nu_tau - e scattering also impacts solar neutrino transformation, although this interaction is more tightly constrained.

We discuss how matter-enhanced active-sterile neutrino transformation in both neutrino and antineutrino channels could enable the production of the rapid neutron capture (r-process) nuclei in neutrino-heated supernova ejecta. In this scheme the lightest sterile neutrino would be heavier than the electron neutrino and split from it by a vacuum mass-squared difference roughly between 3 and 70 eV$^2$ and vacuum mixing angle given by $\sin^2 2\theta_{es} > 10^{-4}$.

We explore neutrino capture during r-process nucleosynthesis in neutrino-driven ejecta from nascent neutron stars. We focus on the interplay between charged-current weak interactions and element synthesis, and we delineate the important role of equilibrium nuclear dynamics. During the period of coexistence of free nucleons and light and/or heavy nuclei, electron neutrino capture inhibits the r-process. At all stages, capture on free neutrons has a larger impact than capture on nuclei. However, neutrino capture on heavey nuclei by itself, if it is very strong, is also detrimental to the r-process until large nuclear equilibrium clusters break down and the classical neutron-capture phase of the r-process begins. The sensitivity of the r-process to neutrino irradiation means that neutrino-capture effects can strongly constrain the r-process site, neutrino physics, or both. These results apply also to r-process scenarios other than neutrino-heated winds.

The neutrinos from a Type II supernova provide perhaps our best opportunity to probe cosmologically interesting muon and/or tauon neutrino masses. This is because matter enhanced neutrino oscillations can lead to an anomalously hot nu_e spectrum, and thus to enhanced charged current cross sections in terrestrial detectors. Two recently proposed supernova neutrino observatories, OMNIS and LAND, will detect neutrons spalled from target nuclei by neutral and charged current neutrino interactions. As this signal is not flavor specific, it is not immediately clear whether a convincing neutrino oscillation signal can be extracted from such experiments. To address this issue we examine the responses of a series of possible light and heavy mass targets, 9Be, 23Na, 35Cl, and 208Pb. We find that strategies for detecting oscillations which use only neutron count rates are problematic at best, even if cross sections are determined by ancillary experiments. Plausible uncertainties in supernova neutrino spectra tend to obscure rate enhancements due to oscillations. However, in the case of 208Pb, a signal emerges that is largely flavor specific and extraordinarily sensitive to the nu_e temperature, the emission of two neutrons. This signal and its flavor specificity are associated with the strength and location of the first-forbidden responses for neutral and charge current reactions, aspects of the 208Pb neutrino cross section that have not been discussed previously. Hadronic spin transfer experiments might be helpful in confirming some of the nuclear structure physics underlying our conclusions.

We study the early evolution of the electron fraction (or, alternatively, the neutron-to-proton ratio) in the region above the hot proto-neutron star formed after a supernova explosion. We study the way in which the electron fraction in this environment is set by a competition between lepton (electron, positron, neutrino, and antineutrino) capture processes on free neutrons and protons and nuclei. Our calculations take explicit account of the effect of nuclear composition changes, such as formation of alpha particles (the alpha effect) and the shifting of nuclear abundances in nuclear statistical equilibrium associated with cooling in near-adiabatic outflow. We take detailed account of the process of weak interaction freeze-out in conjunction with these nuclear composition changes. Our detailed treatment shows that the alpha effect can cause significant increases in the electron fraction, while neutrino and antineutrino capture on heavy nuclei tends to have a buffering effect on this quantity. We also examine the effect on weak rates and the electron fraction of fluctuations in time in the neutrino and antineutrino energy spectra arising from hydrodynamic waves. Our analysis is guided by the Mayle & Wilson supernova code numerical results for the neutrino energy spectra and density and velocity profiles.

We assess the prospects for attaining steady nuclear flow equilibrium in expanding r-process environments where beta decay and/or neutrino capture determine the nuclear charge-changing rates. For very rapid expansions, we find that weak steady flow equilibrium normally cannot be attained. However, even when neutron capture processes freeze out in such nonequilibrium conditions, abundance ratios of nuclear species in the r-process peaks might still mimic those attained in weak steady flow. This result suggests that the r-process yield in a regime of rapid expansion can be calculated reliably only when all neutron capture, photodisintegration, and weak interaction processes are fully coupled in a dynamical calculation. We discuss the implications of these results for models of the r-process sited in rapidly expanding neutrino-heated ejecta.

We investigate the possibility that neutrino capture on heavy nuclei competes with beta decay in the environment where the $r$-Process elements are synthesized. We find that such neutrino capture is not excluded by existing abundance determinations. We show that inclusion of significant neutrino capture on the (neutron number) N=82 waiting point nuclei can allow the inferred abundances of these species to provide a good fit to steady weak (beta decay plus neutrino capture) flow equilibrium. In fact, for particular choices of neutrino flux conditions, this fit is improved over the case where nuclei change their charge by beta decay alone. However, this improved fit can be realized only if neutrino capture plays a negligible role in nuclear decay back toward stability. We discuss the implications of these considerations for current proposed sites and models for $r$-Process nucleosynthesis.