The Late Devonian was a protracted period of low speciation resulting in biodiversity decline, culminating in extinction events near the Devonian-Carboniferous boundary. Recent evidence indicates that the final extinction event may have coincided with a dramatic drop in stratospheric ozone, possibly due to a global temperature rise. Here we study an alternative possible cause for the postulated ozone drop: a nearby supernova explosion that could inflict damage by accelerating cosmic rays that can deliver ionizing radiation for up to $\sim 100$ kyr. We therefore propose that end-Devonian extinction was triggered by one or more supernova explosions at $\sim 20 \ \rm pc$, somewhat beyond the ``kill distance'' that would have precipitated a full mass extinction. Nearby supernovae are likely due to core-collapses of massive stars in clusters in the thin Galactic disk in which the Sun resides. Detecting any of the long-lived radioisotopes \sm146, \u235 or \pu244 in one or more end-Devonian extinction strata would confirm a supernova origin, point to the core-collapse explosion of a massive star, and probe supernova nucleosythesis. Other possible tests of the supernova hypothesis are discussed.

Motivated by the occurrence of a moderately nearby supernova near the beginning of the Pleistocene, possibly as part of a long-term series beginning in the Miocene, we investigate whether nitrate rainout resulting from the atmospheric ionization of enhanced cosmic ray flux could have, through its fertilizer effect, initiated carbon dioxide drawdown. Such a drawdown could possibly reduce the greenhouse effect and induce the climate change that led to the Pleistocene glaciations. We estimate that the nitrogen flux enhancement onto the surface from an event at 50 pc would be of order 10%, probably too small for dramatic changes. We estimate deposition of iron (another potential fertilizer) and find it is also too small to be significant. There are also competing effects of opposite sign, including muon irradiation and reduction in photosynthetic yield caused by UV increase from stratospheric ozone layer depletion, leading to an ambiguous result. However, if the atmospheric ionization induces a large increase in the frequency of lightning, as argued elsewhere, the amount of nitrate synthesis should be much larger, dominate over the other effects, and induce the climate change. More work needs to be done to clarify effects on lightning frequency.

Motivated by recent measurements of deposits of $^{60}$Fe on the ocean floor and the lunar surface, we model the transport of dust grains containing $^{60}$Fe from a near-Earth (i.e., within 100 pc) supernova (SN). We inject dust grains into the environment of an SN remnant (SNR) and trace their trajectories by applying a 1D hydrodynamic description assuming spherical symmetry to describe the plasma dynamics, and include a rudimentary, 3D magnetic field description to examine its influence on charged dust grains. We assume the interstellar medium (ISM) magnetic fields are turbulent, and are amplified by the SNR shock, while the SN wind and ejecta fields are negligible. We examine the various influences on the dust grains within the SNR to determine when/if the dust decouples from the plasma, how much it is sputtered, and where within the SNR the dust grains are located. We find that Rayleigh-Taylor instabilities are important for dust survival, as they influence the location of the SN's reverse shock. We find that the presence of a magnetic field within the shocked ISM material limits the passage of SN dust grains, with the field either reflecting or trapping the grains within the heart of the SNR. These results have important implications for {\it in situ} $^{60}$Fe measurements, and for dust evolution in SNRs generally.

Neutron star mergers (NSMs) are rapid neutron capture ($\textit{r}$-process) nucleosynthesis sites that expel matter at high velocities, from $0.1c$ to as high as $0.6c$. Nuclei ejected at these speeds are sufficiently energetic to initiate spallation nuclear reactions with interstellar medium particles. We adopt a thick-target model for the propagation of high-speed heavy nuclei in the interstellar medium, similar to the transport of cosmic rays. We find that spallation may create observable perturbations to NSM isotopic abundances, particularly around the low-mass edges of the $\textit{r}$-process peaks where neighboring nuclei have very different abundances. The extent to which spallation modifies the final NSM isotopic yields depends on: (1) the ejected abundances, which are determined by the NSM astrophysical conditions and the properties of nuclei far from stability, (2) the ejecta velocity distribution and propagation in interstellar matter, and (3) the spallation cross-sections. Observed solar and stellar $\textit{r}$-process yields could thus constrain the velocity distribution of ejected neutron star matter, assuming NSMs are the dominant $\textit{r}$-process source. We suggest avenues for future work, including measurement of relevant cross sections.

Neutron star mergers (NSMs) are rapid neutron capture (r-process) nucleosynthesis sites, which eject materials at high velocities, from 0.1c to as high as 0.6c. Thus the r-process nuclei ejected from a NSM event are sufficiently energetic to initiate spallation reactions with the interstellar medium (ISM) particles. With a thick-target model for the propagation of high-speed heavy nuclei in the ISM, we find that spallation reactions may shift the r-process abundance patterns towards solar data, particularly around the low-mass edges of the r-process peaks where neighboring nuclei have very different abundances. The spallation effects depend both on the astrophysical conditions of the r-process nuclei and nuclear physics inputs for the nucleosynthesis calculations and the propagation process. This work extends that of [Wang et al.(2019)] by focusing on the influence of nuclear physics variations on spallation effects.

We assess the status of big-bang nucleosynthesis (BBN) in light of the final Planck data release and other recent developments, and in anticipation of future measurements. Planck data fix the cosmic baryon density to 0.9% precision, and determine the helium abundance and effective number of neutrinos with precision approaching that of astronomical and BBN determinations respectively. In addition, new high-redshift measurements give D/H to better precision than theoretical predictions, and new Li/H data reconfirm the lithium problem. We present new ${}^{7}{\rm Be}(n,p){}^{7}{\rm Li}$ rates using new neutron capture measurements; we have also examined the effect of proposed changes in the $d(p,\gamma){}^{3}{\rm He}$ rates. Using these results we perform a series of likelihood analyses. We assess BBN/CMB consistency, with attention to how our results depend on the choice of Planck data, as well as how the results depend on the choice of non-BBN, non-Planck data sets. Most importantly the lithium problem remains, and indeed is more acute given the very tight D/H observational constraints; new neutron capture data reveals systematics that somewhat increases uncertainty and thus slightly reduces but does not essentially change the problem. We confirm that $d(p,\gamma){}^{3}{\rm He}$ theoretical rates brings D/H out of agreement and slightly increases 7Li; new experimental data are needed at BBN energies. Setting the lithium problem aside, we find the effective number of neutrino species at BBN is $N_\nu = 2.86 \pm 0.15$. Future CMB Stage-4 measurements promise substantial improvements in BBN parameters: helium abundance determinations will be competitive with the best astronomical determinations, and $N_{\rm eff}$ will approach sensitivities capable of detecting the effects of Standard Model neutrino heating of the primordial plasma. (Abridged)

The All-sky Medium Energy Gamma-ray Observatory (AMEGO) is a probe class mission concept that will provide essential contributions to multimessenger astrophysics in the late 2020s and beyond. AMEGO combines high sensitivity in the 200 keV to 10 GeV energy range with a wide field of view, good spectral resolution, and polarization sensitivity. Therefore, AMEGO is key in the study of multimessenger astrophysical objects that have unique signatures in the gamma-ray regime, such as neutron star mergers, supernovae, and flaring active galactic nuclei. The order-of-magnitude improvement compared to previous MeV missions also enables discoveries of a wide range of phenomena whose energy output peaks in the relatively unexplored medium-energy gamma-ray band.

A Milky-Way Type Ia Supernova (SNIa) could be unidentified or even initially unnoticed, being dim in radio, X-rays, and neutrinos, and suffering large optical/IR extinction in the Galactic plane. But SNIa emit nuclear gamma-ray lines from $^{56}{\rm Ni}\to ^{56}{\rm Co}\to ^{56}{\rm Fe}$ radioactive decays. These lines fall within the Fermi/GBM energy range, and the $^{56}{\rm Ni}$ 158 keV line is detectable by Swift/BAT. Both instruments frequently monitor the Galactic plane, which is transparent to gamma rays. Thus GBM and BAT are ideal Galactic SNIa early warning systems. We simulate SNIa MeV light curves and spectra to show that GBM and BAT could confirm a Galactic SNIa explosion, followed by Swift localization and observation in X-rays and UVOIR band. The time of detection depends sensitively on the $^{56}{\rm Ni}$ distribution, and can be as early as a few days if $\gtrsim$ 10% of the $^{56}{\rm Ni}$ is present in the surface as suggested by SN2014J gamma data.

There is now solid experimental evidence of at least one supernova explosion within 100 pc of Earth within the last few million years, from measurements of the short-lived isotope 60Fe in widespread deep-ocean samples, as well as in the lunar regolith and cosmic rays. This is the first established example of a specific dated astrophysical event outside the Solar System having a measurable impact on the Earth, offering new probes of stellar evolution, nuclear astrophysics, the astrophysics of the solar neighborhood, cosmic-ray sources and acceleration, multi-messenger astronomy, and astrobiology. Interdisciplinary connections reach broadly to include heliophysics, geology, and evolutionary biology. Objectives for the future include pinning down the nature and location of the established near-Earth supernova explosions, seeking evidence for others, and searching for other short-lived isotopes such as 26Al and 244Pu. The unique information provided by geological and lunar detections of radioactive 60Fe to assess nearby supernova explosions make now a compelling time for the astronomy community to advocate for supporting multi-disciplinary, cross-cutting research programs.

The high-energy universe has revealed that energetic particles are ubiquitous in the cosmos and play a vital role in the cultivation of cosmic environments on all scales. Though they play a key role in cultivating the cosmological environment and/or enabling our studies of it, there is still much we do not know about AGNs and GRBs, particularly the avenue in which and through which they supply radiation and energetic particles, namely their jets. This White Paper is the second of a two-part series highlighting the most well-known high-energy cosmic accelerators and contributions that MeV gamma-ray astronomy will bring to understanding their energetic particle phenomena. The focus of this white paper is active galactic nuclei and gamma-ray bursts.

The high-energy universe has revealed that energetic particles are ubiquitous in the cosmos and play a vital role in the cultivation of cosmic environments on all scales. Energetic particles in our own galaxy, galactic cosmic rays (GCRs), engage in a complex interplay with the interstellar medium and magnetic fields in the galaxy, giving rise to many of its key characteristics. This White Paper is the first of a two-part series highlighting the most well-known high-energy cosmic accelerators and contributions that MeV gamma-ray astronomy will bring to understanding their energetic particle phenomena. The focus of this white paper is galactic cosmic rays, supernova remnants, protostellar jets and superbubbles, and colliding wind binaries.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

We argue that the exclusion of the Galactic Plane and Bulge from the uniform wide-fast-deep (WFD) LSST survey cadence is fundamentally inconsistent with two of the main science drivers of LSST: Mapping the Milky Way and Exploring the Transient Optical Sky. We outline the philosophical basis for this claim and then describe a number of important science goals that can only be addressed by WFD-like coverage of the Plane and Bulge.

Several starburst galaxies have been observed in the GeV and TeV bands. In these dense environments, gamma-ray emission should be dominated by cosmic-ray interactions with the interstellar medium ($p_{\rm cr}p_{\rm ism} \to \pi^{0} \to \gamma\gamma$). Indeed, starbursts may act as proton "calorimeters" where a substantial fraction of cosmic-ray energy input is emitted in gamma rays. Here we build a one-zone, "thick-target" model implementing calorimetry and placing a firm upper bound on gamma-ray emission from cosmic-ray interactions. The model assumes that cosmic rays are accelerated by supernovae (SNe), and all suffer nuclear interactions rather than escape. Our model has only two free parameters: the cosmic-ray proton acceleration energy per supernova $\epsilon_{\rm cr}$, and the proton injection spectral index $s$. We calculate the pionic gamma-ray emission from 10 MeV to 10 TeV, and derive thick-target parameters for six galaxies with Fermi, H.E.S.S., and/or VERITAS data. Our model provides good fits for the M82 and NGC 253, and yields $\epsilon_{\rm cr}$ and $s$ values suggesting that supernova cosmic-ray acceleration is similar in starbursts and in our Galaxy. We find that these starbursts are indeed nearly if not fully proton calorimeters. For NGC 4945 and NGC 1068, the models are consistent with calorimetry but are less well-constrained due to the lack of TeV data. However, the Circinus galaxy and the ultraluminous infrared galaxy Arp 220 exceed our pionic upper-limit, possible explanations are discussed.

Several searches have found evidence of $^{60}$Fe deposition, presumably from a near-Earth supernova (SN), with concentrations that vary in different locations on Earth. This paper examines various influences on the path of interstellar dust carrying $^{60}$Fe from a SN through the heliosphere, with the aim of estimating the final global distribution on the ocean floor. We study the influences of magnetic fields, angle of arrival, wind and ocean cycling of SN material on the concentrations at different locations. We find that the passage of SN material through the mesosphere/lower thermosphere (MLT) is the greatest influence on the final global distribution, with ocean cycling causing lesser alteration as the SN material sinks to the ocean floor. SN distance estimates in previous works that assumed a uniform distribution are a good approximation. Including the effects on surface distributions, we estimate a distance of $46^{+10}_{-6}$ pc for a $8-10 \ M_{\odot}$ SN progenitor. This is consistent with a SN occurring within the Tuc-Hor stellar group $\sim$2.8 Myr ago with SN material arriving on Earth $\sim$2.2 Myr ago. We note that the SN dust retains directional information to within $1^{\circ}$ through its arrival in the inner Solar System, so that SN debris deposition on inert bodies such as the Moon will be anisotropic, and thus could in principle be used to infer directional information. In particular, we predict that existing lunar samples should show measurable $^{60}$Fe differences.

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.

We perform spatial and spectral analyses of the Large Magellanic Cloud (LMC) gamma-ray emission collected over 66 months by the Fermi Gamma-ray Space Telescope. In our spatial analysis, we model the LMC cosmic-ray distribution and gamma-ray production using observed maps of the LMC interstellar medium, star formation history, interstellar radiation field, and synchrotron emission. We use bootstrapping of the data to quantify the robustness of spatial model performance. We model the LMC gamma-ray spectrum using fitting functions derived from the physics of $\pi^0$ decay, bremsstrahlung, and inverse Compton scattering. We find the integrated gamma-ray flux of the LMC from 200 MeV to 20 GeV to be $1.38 \pm 0.02 \times 10^{-7}$ ph cm$^{-2}$ s$^{-1}$, of which we attribute about 10% to inverse Compton scattering and 40% to bremsstrahlung. From our work, we conclude that the spectral index of the LMC cosmic-ray proton population is 2.4$\pm$0.2, and we find that cosmic-ray energy loss through gamma-ray production is concentrated within a few hundred pc of acceleration sites. Assuming cosmic-ray energy equipartition with magnetic fields, we estimate LMC cosmic rays encounter an average magnetic field strength ~3 $\mu$G.

We present a new method to measure or constrain p-wave-suppressed cross sections for dark matter (DM) annihilations inside the steep density spikes induced by supermassive black holes. We demonstrate that the high DM densities, together with the increased velocity dispersion, within such spikes combine to make thermal p-wave annihilation cross-sections potentially visible in gamma-ray observations of the Galactic center (GC). The resulting DM signal is a bright central point source with emission originating from DM annihilations in the absence of a detectable spatially-extended signal from the halo. We define two simple reference theories of DM with a thermal p-wave annihilation cross-section and establish new limits on the combined particle and astrophysical parameter space of these models, demonstrating that Fermi is currently sensitive to thermal p-wave DM over a wide range of possible scenarios for the DM distribution in the GC.

Big-bang nucleosynthesis (BBN) describes the production of the lightest nuclides via a dynamic interplay among the four fundamental forces during the first seconds of cosmic time. We briefly overview the essentials of this physics, and present new calculations of light element abundances through li6 and li7, with updated nuclear reactions and uncertainties including those in the neutron lifetime. We provide fits to these results as a function of baryon density and of the number of neutrino flavors, N_nu. We review recent developments in BBN, particularly new, precision Planck cosmic microwave background (CMB) measurements that now probe the baryon density, helium content, and the effective number of degrees of freedom, n_eff. These measurements allow for a tight test of BBN and of cosmology using CMB data alone. Our likelihood analysis convolves the 2015 Planck data chains with our BBN output and observational data. Adding astronomical measurements of light elements strengthens the power of BBN. We include a new determination of the primordial helium abundance in our likelihood analysis. New D/H observations are now more precise than the corresponding theoretical predictions, and are consistent with the Standard Model and the Planck baryon density. Moreover, D/H now provides a tight measurement of N_nu when combined with the CMB baryon density, and provides a 2sigma upper limit N_nu < 3.2. The new precision of the CMB and of D/H observations together leave D/H predictions as the largest source of uncertainties. Future improvement in BBN calculations will therefore rely on improved nuclear cross section data. In contrast with D/H and he4, li7 predictions continue to disagree with observations, perhaps pointing to new physics.

We consider the production and deposition on Earth of isotopes with half-lives in the range 10$^{5}$ to 10$^{8}$ years that might provide signatures of nearby stellar explosions, extending previous analyses of Core-Collapse Supernovae (CCSNe) to include Electron-Capture Supernovae (ECSNe), Super-Asymptotic Giant Branch (SAGBs) stars, Thermonuclear/Type Ia Supernovae (TNSNe), and Kilonovae/Neutron Star Mergers (KNe). We revisit previous estimates of the $^{60}$Fe and $^{26}$Al signatures, and extend these estimates to include $^{244}$Pu and $^{53}$Mn. We discuss interpretations of the $^{60}$Fe signals in terrestrial and lunar reservoirs in terms of a nearby stellar ejection ~2.2 Myr ago, showing that (i) the $^{60}$Fe yield rules out the TNSN and KN interpretations, (ii) the $^{60}$Fe signals highly constrain a SAGB interpretation but do not completely them rule out, (iii) are consistent with a CCSN origin, and (iv) are highly compatible with an ECSN interpretation. Future measurements could resolve the radioisotope deposition over time, and we use the Sedov blast wave solution to illustrate possible time-resolved profiles. Measuring such profiles would independently probe the blast properties including distance, and would provide additional constraints the nature of the explosion.

Emission from blazar jets in the ultraviolet, optical, and infrared is polarized. If these low-energy photons were inverse-Compton scattered, the upscattered high-energy photons retain a fraction of the polarization. Current and future X-ray and gamma-ray polarimeters such as INTEGRAL-SPI, PoGOLITE, X-Calibur, Gamma-Ray Burst Polarimeter, GEMS-like missions, ASTRO-H, and POLARIX have the potential to discover polarized X-rays and gamma-rays from blazar jets for the first time. Detection of such polarization will open a qualitatively new window into high-energy blazar emission; actual measurements of polarization degree and angle will quantitatively test theories of jet emission mechanisms. We examine the detection prospects of blazars by these polarimetry missions using examples of 3C 279, PKS 1510-089, and 3C 454.3, bright sources with relatively high degrees of low-energy polarization. We conclude that while balloon polarimeters will be challenged to detect blazars within reasonable observational times (with X-Calibur offering the most promising prospects), space-based missions should detect the brightest blazars for polarization fractions down to a few percent. Typical flaring activity of blazars could boost the overall number of polarimetric detections by nearly a factor of five to six purely accounting for flux increase of the brightest of the comprehensive, all-sky, Fermi-LAT blazar distribution. The instantaneous increase in the number of detections is approximately a factor of two, assuming a duty cycle of 20% for every source. The detectability of particular blazars may be reduced if variations in the flux and polarization fraction are anticorrelated. Simultaneous use of variability and polarization trends could guide the selection of blazars for high-energy polarimetric observations.

We give updated constraints on hypothetical light bosons with a two-photon coupling such as axions or axion-like particles (ALPs). We focus on masses and lifetimes where decays happen near big bang nucleosynthesis (BBN), thus altering the baryon-to-photon ratio and number of relativistic degrees of freedom between the BBN epoch and the cosmic microwave background (CMB) last scattering epoch, in particular such that $N_{\rm eff}^{\rm CMB} < N_{\rm eff}^{\rm BBN}$ and $\eta^{\rm CMB} < \eta^{\rm BBN}$. New constraints presented here come from Planck measurements of the CMB power spectrum combined with the latest inferences of primordial $^4$He and D/H abundances. We find that a previously allowed region in parameter space near $m=1\,\rm MeV$ and $\tau=100\,\rm ms$, consistent with a QCD axion arising from a symmetry breaking near the electroweak scale, is now ruled out at $>3\sigma$ by the combination of CMB+D/H measurements if only ALPs and three thermalized neutrino species contribute to $N_{\rm eff}$. The bound relaxes if there are additional light degrees of freedom present which, in this scenario, have their contribution limited to $\Delta N_{\rm eff}=1.1\pm0.3$. We give forecasts showing that a number of experiments are expected to reach the sensitivity needed to further test this region, such as Stage-IV CMB and SUPER-KEKB, the latter a direct test insensitive to any extra degrees of freedom.

A critical review is given of the current status of cosmological nucleosynthesis. In the framework of the Standard Model with 3 types of relativistic neutrinos, the baryon-to-photon ratio, $\eta$, corresponding to the inferred primordial abundances of deuterium and helium-4 is consistent with the independent determination of $\eta$ from observations of anisotropies in the cosmic microwave background. However the primordial abundance of lithium-7 inferred from observations is significantly below its expected value. Taking systematic uncertainties in the abundance estimates into account, there is overall concordance in the range $\eta = (5.7-6.7)\times 10^{-10}$ at 95% CL (corresponding to a cosmological baryon density $\Omega_B h^2 = 0.021 - 0.025$). The D and He-4 abundances, when combined with the CMB determination of $\eta$, provide the bound $N_\nu=3.28 \pm 0.28$ on the effective number of neutrino species. Other constraints on new physics are discussed briefly.

If the supermassive black hole Sgr A* at the center of the Milky Way grew adiabatically from an initial seed embedded in an NFW dark matter (DM) halo, then the DM profile near the hole has steepened into a spike. We calculate the dramatic enhancement to the gamma ray flux from the Galactic center (GC) from such a spike if the 1-3 GeV excess observed in Fermi data is due to DM annihilations. We find that for the parameter values favored in recent fits, the point source-like flux from the spike is 35 times greater than the flux from the inner $1^\circ$ of the halo, far exceeding all Fermi point source detections near the GC. We consider the dependence of the spike signal on astrophysical and particle parameters and conclude that if the GC excess is due to DM, then a canonical adiabatic spike is disfavored by the data. We discuss alternative Galactic histories that predict different spike signals, including: (i) the nonadiabatic growth of the black hole, possibly associated with halo and/or black hole mergers, (ii) gravitational interaction of DM with baryons in the dense core, such as heating by stars, or (iii) DM self-interactions. We emphasize that the spike signal is sensitive to a different combination of particle parameters than the halo signal, and that the inclusion of a spike component to any DM signal in future analyses would provide novel information about both the history of the GC and the particle physics of DM annihilations.

We studied previously the impact on light-element abundances of gravitinos decaying during or after Big-Bang nucleosynthesis (BBN). We found regions of the gravitino mass m_{3/2} and abundance zeta_{3/2} plane where its decays could reconcile the calculated abundance of Li7 with observation without perturbing the other light-element abundances unacceptably. Here we revisit this issue in light of LHC measurements of the Higgs mass and constraints on supersymmetric model parameters, as well as updates in the astrophysical measurements of light-element abundances. In addition to the constrained minimal supersymmetric extension of the Standard Model with universal soft supersymmetry-breaking masses at the GUT scale (the CMSSM) studied previously, we also study models with universality imposed below the GUT scale and models with non-universal Higgs masses (NUHM1). We calculate the total likelihood function for the light-element abundances, taking into account the observational uncertainties. We find that gravitino decays provide a robust solution to the cosmological Li7 problem along strips in the (m_{3/2}, zeta_{3/2}) plane along which the abundances of deuterium, He4 and Li7 may be fit with chi^2_min < 3, compared with chi^2 ~ 34 if the effects of gravitino decays are unimportant. The minimum of the likelihood function is reduced to chi^2 < 2 when the uncertainty on D/H is relaxed and < 1 when the lithium abundance is taken from globular cluster data.

The primordial abundances of light elements produced in the standard theory of Big Bang nucleosynthesis (BBN) depend only on the cosmic ratio of baryons to photons, a quantity inferred from observations of the microwave background. The predicted primordial 7Li abundance is four times that measured in the atmospheres of Galactic halo stars. This discrepancy could be caused by modification of surface lithium abundances during the stars' lifetimes or by physics beyond the Standard Model that affects early nucleosynthesis. The lithium abundance of low-metallicity gas provides an alternative constraint on the primordial abundance and cosmic evolution of lithium that is not susceptible to the in situ modifications that may affect stellar atmospheres. Here we report observations of interstellar 7Li in the low-metallicity gas of the Small Magellanic Cloud, a nearby galaxy with a quarter the Sun's metallicity. The present-day 7Li abundance of the Small Magellanic Cloud is nearly equal to the BBN predictions, severely constraining the amount of possible subsequent enrichment of the gas by stellar and cosmic-ray nucleosynthesis. Our measurements can be reconciled with standard BBN with an extremely fine-tuned depletion of stellar Li with metallicity. They are also consistent with non-standard BBN.

We consider the effects of metastable charged sparticles on Big-Bang Nucleosynthesis (BBN), including bound-state reaction rates and chemical effects. We make a new analysis of the bound states of negatively-charged massive particles with the light nuclei most prominent in BBN, and present a new code to track their abundances, paying particular attention to that of Li7. Assuming, as an example, that the gravitino is the lightest supersymmetric particle (LSP), and that the lighter stau slepton, stau_1, is the metastable next-to-lightest sparticle within the constrained minimal supersymmetric extension of the Standard Model (CMSSM), we analyze the possible effects on the standard BBN abundances of stau_1 bound states and decays for representative values of the gravitino mass. Taking into account the constraint on the CMSSM parameter space imposed by the discovery of the Higgs boson at the LHC, we delineate regions in which the fit to the measured light-element abundances is as good as in standard BBN. We also identify regions of the CMSSM parameter space in which the bound state properties, chemistry and decays of metastable charged sparticles can solve the cosmological Li7 problem.

Fermi has resolved several star-forming galaxies, but the vast majority of the star-forming universe is unresolved and thus contributes to the extragalactic gamma ray background (EGB). Here, we calculate the contribution from star-forming galaxies to the EGB in the Fermi range from 100 MeV to 100 GeV, due to inverse-Compton (IC) scattering of the interstellar photon field by cosmic-ray electrons. We first construct a one-zone model for a single star-forming galaxy, assuming supernovae power the acceleration of cosmic rays. The same IC interactions leading to gamma rays also substantially contribute to the energy loss of the high-energy cosmic-ray electrons. Consequently, a galaxy's IC emission is determined by the relative importance of IC losses in the cosmic-ray electron energy budget ("partial calorimetry"). We use our template for galactic IC luminosity to find the cosmological contribution of star-forming galaxies to the EGB. For all of our models, we find the IC EGB contribution is almost an order of magnitude less than the peak of the emission due to cosmic-ray ion interactions (mostly pionic p_cr p_ism \rightarrow \pi_0 \rightarrow \gamma \gamma); even at the highest Fermi energies, IC is subdominant. Moreover, the flatter IC spectrum increases the high-energy signal of the pionic+IC sum, bringing it into better agreement with the EGB spectral index observed by Fermi . Partial calorimetry ensures that the overall IC signal is well constrained, with only modest uncertainties in the amplitude and spectral shape for plausible model choices. Partial calorimetry of cosmic-ray electrons should hold true in both normal and starburst galaxies, and thus we include starbursts in our calculation. We conclude with a brief discussion on how the pionic spectral feature and other methods can be used to measure the star-forming component of the EGB.

Big-bang nucleosynthesis (BBN) theory, together with the precise WMAP cosmic baryon density, makes tight predictions for the abundances of the lightest elements. Deuterium and 4He measurements agree well with expectations, but 7Li observations lie a factor 3-4 below the BBN+WMAP prediction. This 4-5\sigma\ mismatch constitutes the cosmic "lithium problem," with disparate solutions possible. (1) Astrophysical systematics in the observations could exist but are increasingly constrained. (2) Nuclear physics experiments provide a wealth of well-measured cross-section data, but 7Be destruction could be enhanced by unknown or poorly-measured resonances, such as 7Be + 3He -> 10C^* -> p + 9B. (3) Physics beyond the Standard Model can alter the 7Li abundance, though D and 4He must remain unperturbed; we discuss such scenarios, highlighting decaying Supersymmetric particles and time-varying fundamental constants. Present and planned experiments could reveal which (if any) of these is the solution to the problem.

The origin of the diffuse extragalactic gamma-ray background (EGB) has been intensively studied but remains unsettled. Current popular source candidates include unresolved star-forming galaxies, starburst galaxies, and blazars. In this paper we calculate the EGB contribution from the interactions of cosmic rays accelerated by Type Ia supernovae (SNe), extending earlier work which only included core-collapse SNe. We consider Type Ia events in star-forming galaxies, but also in quiescent galaxies that lack star formation. For star-forming galaxies, consistently including Type Ia events makes little change to the star-forming EGB prediction, so long as both SN types have the same cosmic-ray acceleration efficiencies in star-forming galaxies. Thus, our updated EGB estimate continues to show that star-forming galaxies can represent a substantial portion of the signal measured by Fermi. For quiescent galaxies, conversely, we find a wide range of possibilities for the EGB contribution. The dominant uncertainty we investigated comes from the mass in hot gas, which provides targets for cosmic rays; total gas masses are as yet poorly known, particularly at larger radii. Additionally, the EGB estimation is very sensitive to the cosmic-ray acceleration efficiency and confinement, especially in quiescent galaxies. In the most optimistic allowed scenarios, quiescent galaxies can be an important source of the EGB. In this case, star-forming galaxies and quiescent galaxies together will dominate the EGB and leave little room for other contributions. If other sources, such as blazars, are found to have important contributions to the EGB, then either the gas mass or cosmic-ray content of quiescent galaxies must be significantly lower than in their star-forming counterparts. In any case, improved Fermi EGB measurements will provide important constraints on hot gas and cosmic rays in quiescent galaxies.

Residual late-time dark matter particle annihilations during and after Big Bang Nucleosynthesis (BBN) may alter the predicted cosmological abundances of the light elements. Within the constrained minimal supersymmetric extension of the Standard Model (the CMSSM) with a neutralino LSP, we find negligible effects on the abundances of Deuterium, He3, He4 and Li7 predicted by homogeneous BBN, but potentially a large enhancement in the predicted abundance of Li6. This enhancement may be as much as two orders of magnitude in the focus-point WMAP strip and in the coannihilation and funnel regions for large tan beta for small m_{1/2}, and the effect is still significant at large m_{1/2}. However, the potential Li6 enhancement is negligible in the part of the coannihilation strip for tan beta = 10 that survives the latest LHC constraints. A similar enhancement of the \li6 abundance may also be found in a model with common, non-universal Higgs masses (the NUHM1).

Radio properties of supernova outbursts remain poorly understood despite longstanding campaigns following events discovered at other wavelengths. After ~ 30 years of observations, only ~ 50 supernovae have been detected at radio wavelengths, none of which are Type Ia. Even the most radio-loud events are ~ 10^4 fainter in the radio than in the optical; to date, such intrinsically dim objects have only been visible in the very local universe. The detection and study of radio supernovae (RSNe) will be fundamentally altered and dramatically improved as the next generation of radio telescopes comes online, including EVLA, ASKAP, and MeerKAT, and culminating in the Square Kilometer Array (SKA); the latter should be > 50 times more sensitive than present facilities. SKA can repeatedly scan large (> 1 deg^2) areas of the sky, and thus will discover RSNe and other transient sources in a new, automatic, untargeted, and unbiased way. We estimate SKA will be able to detect core-collapse RSNe out to redshift z ~ 5, with an all-redshift rate ~ 620 events yr^-1 deg^-2, assuming a survey sensitivity of 50 nJy and radio lightcurves like those of SN 1993J. Hence SKA should provide a complete core-collapse RSN sample that is sufficient for statistical studies of radio properties of core-collapse supernovae. EVLA should find ~ 160 events yr^-1 deg^-2 out to redshift z ~ 3, and other SKA precursors should have similar detection rates. We also provided recommendations of the survey strategy to maximize the RSN detections of SKA. This new radio core-collapse supernovae sample will complement the detections from the optical searches, such as the LSST, and together provide crucial information on massive star evolution, supernova physics, and the circumstellar medium, out to high redshift. Additionally, SKA may yield the first radio Type Ia detection via follow-up of nearby events discovered at other wavelengths.

We have conducted Combined Array for Research in Millimeter-wave Astronomy (CARMA) observations of LiH, in absorption, toward three quasars. These quasars, B0218+357, PKS1830-211, and PKS0201+113, have redshifts of z = 0.685 - 3.387$, and shift the LiH J=1-0 transition to the 1 mm and 3mm wavelength bands, where atmospheric absorption is sharply reduced from that predominating near the rest frequency of 443 GHz. We report a 3$\sigma$ detection of LiH toward B0218+357 with a column density of 1.4x10^{12} cm^{-2} and place an upper limit on the ^6Li/^7Li ratio of <0.16. LiH was not detected toward any other source.

We explore a nuclear physics resolution to the discrepancy between the predicted standard big-bang nucleosynthesis (BBN) abundance of 7Li and its observational determination in metal-poor stars. The theoretical 7Li abundance is 3-4 times greater than the observational values, assuming the baryon-to-photon ratio, eta_wmap, determined by WMAP. The 7Li problem could be resolved within the standard BBN picture if additional destruction of A=7 isotopes occurs due to new nuclear reaction channels or upward corrections to existing channels. This could be achieved via missed resonant nuclear reactions, which is the possibility we consider here. We find some potential candidate resonances which can solve the lithium problem and specify their required resonant energies and widths. For example, a 1^- or 2^- excited state of 10C sitting at approximately 15.0 MeV above its ground state with an effective width of order 10 keV could resolve the 7Li problem; the existence of this excited state needs experimental verification. Other examples using known states include 7Be+t \rightarrow 10B(18.80 MeV), and 7Be+d \rightarrow 9B(16.71 MeV). For all of these states, a large channel radius (a > 10 fm) is needed to give sufficiently large widths. Experimental determination of these reaction strengths is needed to rule out or confirm these nuclear physics solutions to the lithium problem.

The origin of the extragalactic gamma-ray background is a pressing cosmological mystery. The Fermi Gamma-Ray Space Telescope has recently measured the intensity and spectrum of this background; both are substantially different from previous measurements. We present a novel calculation of the gamma-ray background from normal star-forming galaxies. Contrary to long-standing expectations, we find that numerous but individually faint normal galaxies may comprise the bulk of the Fermi signal, rather than rare but intrinsically bright active galaxies. This result has wide-ranging implications, including: the possibility to probe the cosmic star formation history with gamma rays; the ability to infer the cosmological evolution of cosmic rays and galactic magnetic fields; and an increased likelihood of identifying subdominant components from rare sources (e.g., dark matter clumps) through their large anisotropy.

We explore the consequences of an early population of intermediate mass stars in the 2 - 8 M\odot range on cosmic chemical evolution. We discuss the implications of this population as it pertains to several cosmological and astrophysical observables. For example, some very metal-poor galactic stars show large enhancements of carbon, typical of the C-rich ejecta of low-mass stars but not of supernovae; moreover, halo star carbon and oxygen abundances show wide scatter, which imply a wide range of star-formation and nucleosynthetic histories contributed to the first generations of stars. Also, recent analyses of the 4He abundance in metal-poor extragalactic H II regions suggest an elevated abundance Yp \simeq 0.256 by mass, higher than the predicted result from big bang nucleosynthesis assuming the baryon density determined by WMAP, Yp = 0.249. Although there are large uncertainties in the observational determination of 4He, this offset may suggest a prompt initial enrichment of 4He in early metal-poor structures. We also discuss the effect of intermediate mass stars on global cosmic evolution, the reionization of the Universe, the density of white dwarfs, as well as SNII and SNIa rates at high redshift. We also comment on the early astration of D and 7Li. We conclude that if intermediate mass stars are to be associated with Population III stars, their relevance is limited (primarily from observed abundance patterns) to low mass structures involving a limited fraction of the total baryon content of the Universe.

Observational and theoretical evidence suggests that high-energy Galactic cosmic rays are primarily accelerated by supernova remnants. If also true for low-energy cosmic rays, the ionization rate near a supernova remnant should be higher than in the general Galactic interstellar medium (ISM). We have searched for H3+ absorption features in 6 sight lines which pass through molecular material near IC 443---a well-studied case of a supernova remnant interacting with its surrounding molecular material---for the purpose of inferring the cosmic-ray ionization rate in the region. In 2 of the sight lines (toward ALS 8828 and HD 254577) we find large H3+ column densities, N(H3+)~3*10^14 cm^-2, and deduce ionization rates of zeta_2~2*10^-15 s^-1, about 5 times larger than inferred toward average diffuse molecular cloud sight lines. However, the 3 sigma upper limits found for the other 4 sight lines are consistent with typical Galactic values. This wide range of ionization rates is likely the result of particle acceleration and propagation effects, which predict that the cosmic-ray spectrum and thus ionization rate should vary in and around the remnant. While we cannot determine if the H3+ absorption arises in post-shock (interior) or pre-shock (exterior) gas, the large inferred ionization rates suggest that IC 443 is in fact accelerating a large population of low-energy cosmic rays. Still, it is unclear whether this population can propagate far enough into the ISM to account for the ionization rate inferred in diffuse Galactic sight lines.

We consider the effects of uncertainties in nuclear reaction rates on the cosmological constraints on the decays of unstable particles during or after Big-Bang nucleosynthesis (BBN). We identify the nuclear reactions due to non-thermal hadrons that are the most important in perturbing standard BBN, then quantify the uncertainties in these reactions and in the resulting light-element abundances. These results also indicate the key nuclear processes for which improved cross section data would allow different light-element abundances to be determined more accurately, thereby making possible more precise probes of BBN and evaluations of the cosmological constraints on unstable particles. Applying this analysis to models with unstable gravitinos decaying into neutralinos, we calculate the likelihood function for the light-element abundances measured currently, taking into account the current experimental errors in the determinations of the relevant nuclear reaction rates. We find a region of the gravitino mass and abundance in which the abundances of deuterium, He4 and Li7 may be fit with chi^2 = 5.5, compared with chi^2 = 31.7 if the effects of gravitino decays are unimportant. The best-fit solution is improved to chi^2 ~ 2.0 when the lithium abundance is taken from globular cluster data. Some such re-evaluation of the observed light-element abundances and/or nuclear reaction rates would be needed if this region of gravitino parameters is to provide a complete solution to the cosmological Li7 problem.

As the Galaxy evolves, the abundance of deuterium in the interstellar medium (ISM) decreases from its primordial value: deuterium is "astrated". The deuterium astration factor, f_D, the ratio of the primordial D abundance (the D to H ratio by number) to the ISM D abundance, is determined by the competition between stellar destruction and infall, providing a constraint on models of the chemical evolution of the Galaxy. Although conventional wisdom suggests that the local ISM (i.e., within ~1-2 kpc of the Sun) should be well mixed and homogenized on timescales short compared to the chemical evolution timescale, the data reveal gas phase variations in the deuterium, iron, and other metal abundances as large as factors of ~4-5 or more, complicating the estimate of the "true" ISM D abundance and of the deuterium astration factor. Here, assuming that the variations in the observationally inferred ISM D abundances result entirely from the depletion of D onto dust, rather than from unmixed accretion of nearly primordial material, a model-independent, Bayesian approach is used to determine the undepleted abundance of deuterium in the ISM (or, a lower limit to it). We find the best estimate for the undepleted, ISM deuterium abundance to be D/H_ISM >= (2.0+-0.1)x10^(-5). This result is used to provide an estimate of (or, an upper bound to) the deuterium astration factor, f_D = (D/H)_P/D/H_ISM <= 1.4+-0.1.

The cumulative (anti)neutrino production from all core-collapse supernovae within our cosmic horizon gives rise to the diffuse supernova neutrino background (DSNB), which is on the verge of detectability. The observed flux depends on supernova physics, but also on the cosmic history of supernova explosions; currently, the cosmic supernova rate introduces a substantial (+/-40%) uncertainty, largely through its absolute normalization. However, a new class of wide-field, repeated-scan (synoptic) optical sky surveys is coming online, and will map the sky in the time domain with unprecedented depth, completeness, and dynamic range. We show that these surveys will obtain the cosmic supernova rate by direct counting, in an unbiased way and with high statistics, and thus will allow for precise predictions of the DSNB. Upcoming sky surveys will substantially reduce the uncertainties in the DSNB source history to an anticipated +/-5% that is dominated by systematics, so that the observed high-energy flux thus will test supernova neutrino physics. The portion of the universe (z < 1) accessible to upcoming sky surveys includes the progenitors of a large fraction (~ 87%) of the expected 10-26 MeV DSNB event rate. We show that precision determination of the (optically detected) cosmic supernova history will also make the DSNB into a strong probe of an extra flux of neutrinos from optically invisible supernovae, which may be unseen either due to unexpected large dust obscuration in host galaxies, or because some core-collapse events proceed directly to black hole formation and fail to give an optical outburst.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

The decays of massive gravitinos into neutralino dark matter particles and Standard Model secondaries during or after Big-Bang nucleosynthesis (BBN) may alter the primordial light-element abundances. We present here details of a new suite of codes for evaluating such effects, including a new treatment based on PYTHIA of the evolution of showers induced by hadronic decays of massive, unstable particles such as a gravitino. We also develop an analytical treatment of non-thermal hadron propagation in the early universe, and use this to derive analytical estimates for light-element production and in turn on decaying particle lifetimes and abundances. We then consider specifically the case of an unstable massive gravitino within the constrained minimal supersymmetric extension of the Standard Model (CMSSM). We present upper limits on its possible primordial abundance before decay for different possible gravitino masses, with CMSSM parameters along strips where the lightest neutralino provides all the astrophysical cold dark matter density. We do not find any CMSSM solution to the cosmological Li7 problem for small m_{3/2}. Discounting this, for m_{1/2} ~ 500 GeV and tan beta = 10 the other light-element abundances impose an upper limit m_{3/2} n_{3/2}/n_\gamma < 3 \times 10^{-12} GeV to < 2 \times 10^{-13} GeV for m_{3/2} = 250 GeV to 1 TeV, which is similar in both the coannihilation and focus-point strips and somewhat weaker for tan beta = 50, particularly for larger m_{1/2}. The constraints also weaken in general for larger m_{3/2}, and for m_{3/2} > 3 TeV we find a narrow range of m_{3/2} n_{3/2}/n_\gamma, at values which increase with m_{3/2}, where the Li7 abundance is marginally compatible with the other light-element abundances.

The identification of point sources poses a great challenge for the high energy community. We present a new approach to evaluate the likelihood of a set of sources being a Galactic population based on the simple assumption that galaxies similar to the Milky Way host comparable populations of gamma-ray emitters. We propose a luminosity constraint on Galactic source populations which complements existing approaches by constraining the abundance and spatial distribution of any objects of Galactic origin, rather than focusing on the properties of a specific candidate emitter. We use M31 as a proxy for the Milky Way, and demonstrate this technique by applying it to the unidentified EGRET sources. We find that it is highly improbable that the majority of the unidentified EGRET sources are members of a Galactic halo population (e.g., dark matter subhalos), but that current observations do not provide any constraints on all of these sources being Galactic objects if they reside entirely in the disk and bulge. Applying this method to upcoming observations by the Fermi Gamma-ray Space Telescope has the potential to exclude association of an even larger number of unidentified sources with any Galactic source class.

Large synoptic (repeated scan) imaging sky surveys are poised to observe enormous numbers of core-collapse supernovae. We quantify the discovery potential of upcoming projects, including DES, Pan-STARRS, and LSST. These surveys will map out the cosmic core-collapse supernova redshift distribution via direct counting, with very small statistical uncertainties out to a redshift depth which is a strong function of the survey limiting magnitude. This supernova redshift history encodes rich information about cosmology, star formation, and supernova astrophysics and phenomenology; the large statistics of the supernova sample will be crucial to disentangle possible degeneracies among these issues. For example, the cosmic supernova rate can be measured to high precision out to z ~ 0.5 for all core-collapse types, and out to redshift z ~ 1 for Type IIn events. Because of the tight link between supernovae and star formation, synoptic sky surveys will provide precision measurements of the normalization and z < 1 history of cosmic star-formation rate in a manner independent of and complementary to current methods. Furthermore, Type II supernovae can serve as distance indicators and would independently cross-check Type Ia distances measured in the same surveys. Arguably the largest and least-controlled uncertainty in these efforts comes from dust obscuration of supernovae in their host galaxies; we outline a strategy to determine empirically the obscuration properties by leveraging the large supernova samples over a broad range of redshift. [abstract abridged]

Diffuse interstellar clouds show large abundances of H_3^+ which can be maintained only by a high ionization rate of H_2. Cosmic rays are the dominant ionization mechanism in this environment, so the large ionization rate implies a high cosmic-ray flux, and a large amount of energy residing in cosmic rays. In this paper we find that the standard propagated cosmic-ray spectrum predicts an ionization rate much lower than that inferred from H_3^+. Low-energy (~10 MeV) cosmic rays are the most efficient at ionizing hydrogen, but cannot be directly detected; consequently, an otherwise unobservable enhancement of the low-energy cosmic-ray flux offers a plausible explanation for the H_3^+ results. Beyond ionization, cosmic rays also interact with the interstellar medium by spalling atomic nuclei and exciting atomic nuclear states. These processes produce the light elements Li, Be, and B, as well as gamma-ray lines. To test the consequences of an enhanced low-energy cosmic-ray flux, we adopt two physically-motivated cosmic-ray spectra which by construction reproduce the ionization rate inferred in diffuse clouds, and investigate the implications of these spectra on dense cloud ionization rates, light element abundances, gamma-ray fluxes, and energetics. One spectrum proposed here provides an explanation for the high ionization rate seen in diffuse clouds while still appearing to be broadly consistent with other observables, but the shape of this spectrum suggests that supernovae remnants may not be the predominant accelerators of low-energy cosmic rays.

Because of the dynamic nature of the interstellar medium, the Sun should have encountered a variety of different interstellar environments in its lifetime. As the solar wind interacts with the surrounding interstellar medium to form a heliosphere, different heliosphere shapes, sizes, and particle contents result from the different environments. Some of the large possible interstellar parameter space (density, velocity, temperature) is explored here with the help of global heliosphere models, and the features in the resulting heliospheres are compared and discussed. The heliospheric size, expressed as distance of the nose of the heliopause to the Sun, is set by the solar wind - interstellar pressure balance, even for extreme cases. Other heliospheric boundary locations and neutral particle results correlate with the interstellar parameters as well. If the H0 clouds identified in the Millennium Arecibo survey are typical of clouds encountered by the Sun, then the Sun spends ~99.4% of the time in warm low density ISM, where the typical upwind heliosphere radii are up to two orders of magnitude larger than at present.

The lithium problem arises from the significant discrepancy between the primordial 7Li abundance as predicted by BBN theory and the WMAP baryon density, and the pre-Galactic lithium abundance inferred from observations of metal-poor (Population II) stars. This problem has loomed for the past decade, with a persistent discrepancy of a factor of 2--3 in 7Li/H. Recent developments have sharpened all aspects of the Li problem. Namely: (1) BBN theory predictions have sharpened due to new nuclear data, particularly the uncertainty on 3He(alpha,gamma)7Be, has reduced to 7.4%, and with a central value shift of ~ +0.04 keV barn. (2) The WMAP 5-year data now yields a cosmic baryon density with an uncertainty reduced to 2.7%. (3) Observations of metal-poor stars have tested for systematic effects, and have reaped new lithium isotopic data. With these, we now find that the BBN+WMAP predicts 7Li/H = (5.24+0.71-0.67) 10^{-10}. The Li problem remains and indeed is exacerbated; the discrepancy is now a factor 2.4--4.3 or 4.2sigma (from globular cluster stars) to 5.3sigma (from halo field stars). Possible resolutions to the lithium problem are briefly reviewed, and key nuclear, particle, and astronomical measurements highlighted.

Measurements of deuterium in the local interstellar medium have revealed large variations in D/H along different lines of sight. Moreover, recent {\it Far Ultraviolet Spectroscopic Explorer} (FUSE) measurements find D/H to be anticorrelated with several indicators of dust formation and survival, suggesting that interstellar deuterium suffers significant depletion onto dust grains. This in turn implies that the total deuterium abundance in the local Galactic disk could be as high as $\sim 84 %$ of the primordial D abundance. It was proposed that the infall/accretion of pristine gas is needed to explain such a high deuterium abundance. However, we point out that the infall needed to maintain a high present-day D/H is, within the preferred models, in tension with observations that gas represents only some $\sim 20%$ of Galactic baryons, with the balance in stars. We study this tension in the context of a wide class of Galactic evolution models for baryonic processing through stars, which show that deuterium destruction is strongly and cleanly correlated with the drop in the gas fraction. We find that FUSE deuterium observations and Galactic gas fraction estimates can be reconciled in some models; these demand a significant infall rate of pristine material that almost completely balances the rate of star formation. These successful models also require that the average fraction of gas that is returned by dying stars is less than 40% of the initial stellar mass. Cosmological implications of dust depletion of D in high-redshift systems are discussed.

The large majority of EGRET point sources remain without an identified low-energy counterpart, and a large fraction of these sources are most likely extragalactic. Whatever the nature of the extragalactic EGRET unidentified sources, faint unresolved objects of the same class must have a contribution to the diffuse extragalactic gamma-ray background (EGRB). Understanding this component of the EGRB, along with other guaranteed contributions from known sources, is essential if we are to use this emission to constrain exotic high-energy physics. Here, we follow an empirical approach to estimate whether a potential contribution of unidentified sources to the EGRB is likely to be important, and we find that it is. Additionally, we show how upcoming GLAST observations of EGRET unidentified sources, as well as of their fainter counterparts, can be combined with GLAST observations of the Galactic and extragalactic diffuse backgrounds to shed light on the nature of the EGRET unidentified sources even without any positional association of such sources with low-energy counterparts.

Recently Supernova 2006gy was noted as the most luminous ever recorded, with a total radiated energy of ~10^44 Joules. It was proposed that the progenitor may have been a massive evolved star similar to eta Carinae, which resides in our own galaxy at a distance of about 2.3 kpc. eta Carinae appears ready to detonate. Although it is too distant to pose a serious threat as a normal supernova, and given its rotation axis is unlikely to produce a Gamma-Ray Burst oriented toward the Earth, eta Carinae is about 30,000 times nearer than 2006gy, and we re-evaluate it as a potential superluminous supernova. We find that given the large ratio of emission in the optical to the X-ray, atmospheric effects are negligible. Ionization of the atmosphere and concomitant ozone depletion are unlikely to be important. Any cosmic ray effects should be spread out over ~10^4 y, and similarly unlikely to produce any serious perturbation to the biosphere. We also discuss a new possible effect of supernovae, endocrine disruption induced by blue light near the peak of the optical spectrum. This is a possibility for nearby supernovae at distances too large to be considered "dangerous" for other reasons. However, due to reddening and extinction by the interstellar medium, eta Carinae is unlikely to trigger such effects to any significant degree.

Cosmic structure formation leads to large-scale shocked baryonic flows which are expected to produce a cosmological population of structure-formation cosmic rays (SFCRs). Interactions between SFCRs and ambient baryons will produce lithium isotopes via \alpha+\alpha \to ^{6,7}Li. This pre-Galactic (but non-primordial) lithium should contribute to the primordial 7Li measured in halo stars and must be subtracted in order to arrive to the true observed primordial lithium abundance. In this paper we point out that the recent halo star 6Li measurements can be used to place a strong constraint to the level of such contamination, because the exclusive astrophysical production of 6Li is from cosmic-ray interactions. We find that the putative 6Li plateau, if due to pre-Galactic cosmic-ray interactions, implies that SFCR-produced lithium represents Li_{SFCR}/Li_{plateau}\approx 15% of the observed elemental Li plateau. Taking the remaining plateau Li to be cosmological 7Li, we find a revised (and slightly worsened) discrepancy between the Li observations and Big Bang Nucleosynthesis predictions by a factor of ^7Li_{BBN}/^7Li_{plateau} \approx 3.7. Moreover, SFCRs would also contribute to the extragalactic gamma-ray background (EGRB) through neutral pion production. This gamma-ray production is tightly related to the amount of lithium produced by the same cosmic rays; the 6Li plateau limits the pre-Galactic (high-redshift) SFCR contribution to be at the level of I_{\pi_{\gamma}SFCR}/I_{EGRB} < 5% of the currently observed EGRB.

We present an analytical description of the energetics of the population of cosmic accretion shocks, for a concordance cosmology. We calculate how the shock-processed accretion power and mass current are distributed among different shock Mach numbers, and how they evolve with cosmic time. We calculate the cumulative energy input of cosmic accretion shocks of any Mach number to the intergalactic medium as a function of redshift, and we compare it with the energy output of supernova explosions as well as with the energy input required to reionize the universe. In addition, we investigate and quantify the effect of environmental factors, such as local clustering properties and filament preheating on the statistical properties of these shocks. We find that the energy processed by accretion shocks is higher than the supernova energy output for z<3 and that it becomes more than an order of magnitude higher in the local universe. The energy processed by accretion shocks alone becomes comparable to the energy required to reionize the universe by z~3.5. Finally, we establish both qualitative and quantitatively that both local clustering as well as filament compression and preheating are important factors in determining the statistical properties of the cosmic accretion shock population.

The large majority of EGRET point sources remain to this day without an identified low-energy counterpart. Whatever the nature of the EGRET unidentified sources, faint unresolved objects of the same class must have a contribution to the diffuse gamma-ray background: if most unidentified objects are extragalactic, faint unresolved sources of the same class contribute to the background, as a distinct extragalactic population; on the other hand, if most unidentified sources are Galactic, their counterparts in external galaxies will contribute to the unresolved emission from these systems. Understanding this component of the gamma-ray background, along with other guaranteed contributions from known sources, is essential in any attempt to use gamma-ray observations to constrain exotic high-energy physics. Here, we follow an empirical approach to estimate whether a potential contribution of unidentified sources to the extragalactic gamma-ray background is likely to be important, and we find that it is. Additionally, we comment on how the anticipated GLAST measurement of the diffuse gamma-ray background will change, depending on the nature of the majority of these sources.

The third EGRET catalog contains a large number of unidentified sources. This subset of objects is expected to include known gamma-ray emitters of Galactic origin such as pulsars and supernova remnants, in addition to an extragalactic population of blazars. However, current data allows the intriguing possibility that some of these objects may represent a new class of yet undiscovered gamma-ray sources. Many theoretically motivated candidate emitters (e.g. clumps of annihilating dark matter particles) have been suggested to account for these detections. We take a new approach to determine to what extent this population is Galactic and to investigate the nature of the possible Galactic component. By assuming that galaxies similar to the Milky Way should host comparable populations of objects, we constrain the allowed Galactic abundance and distribution of various classes of gamma-ray sources using the EGRET data set. We find it is highly improbable that a large number of the unidentified sources are members of a Galactic halo population, but that a distribution of the sources entirely in the disk and bulge is plausible. Finally, we discuss the additional constraints and new insights that GLAST will provide.

A critical review is given of the current status of cosmological nucleosynthesis. In the framework of the Standard Model with 3 types of relativistic neutrinos, the baryon-to-photon ratio, $\eta$, corresponding to the inferred primordial abundances of deuterium and helium-4 is consistent with the independent determination of $\eta$ from WMAP observations of anisotropies in the cosmic microwave background. However the primordial abundance of lithium-7 inferred from observations is significantly below its expected value. Taking systematic uncertainties in the abundance estimates into account, there is overall concordance in the range $\eta = (4.7 - 6.5) x 10^{-10}$ @ 95% c.l. (corresponding to a cosmological baryon density $\Omega_B h^2$ = 0.017 - 0.024). The D and He-4 abundances, together with the CMB determination of $\eta$, provide the bound $N_\nu = 3.24 \pm 1.2$ @ 95% c.l. on the effective number of neutrino species. Other constraints on new physics are discussed briefly.

If the gravitino is the lightest supersymmetric particle and the long-lived next-to-lightest sparticle (NSP) is the stau, the charged partner of the tau lepton, it may be metastable and form bound states with several nuclei. These bound states may affect the cosmological abundances of Li6 and Li7 by enhancing nuclear rates that would otherwise be strongly suppressed. We consider the effects of these enhanced rates on the final abundances produced in Big-Bang nucleosynthesis (BBN), including injections of both electromagnetic and hadronic energy during and after BBN. We calculate the dominant two- and three-body decays of both neutralino and stau NSPs, and model the electromagnetic and hadronic decay products using the PYTHIA event generator and a cascade equation. Generically, the introduction of bound states drives light element abundances further from their observed values; however, for small regions of parameter space bound state effects can bring lithium abundances in particular in better accord with observations. We show that in regions where the stau is the NSP with a lifetime longer than 10^3-10^4 s, the abundances of Li6 and Li7 are far in excess of those allowed by observations. For shorter lifetimes of order 1000 s, we comment on the possibility in minimal supersymmetric and supergravity models that stau decays could reduce the Li7 abundance from standard BBN values while at the same time enhancing the Li6 abundance.

We construct a simple model for radioisotopic enrichment of the protosolar nebula by injection from a nearby supernova, based on the inverse square law for ejecta dispersion. We find that the presolar radioisotopes abundances (i.e., in solar masses) demand a nearby supernova: its distance can be no larger than 66 times the size of the protosolar nebula, at a 90% confidence level, assuming 1 solar mass of protosolar material. The relevant size of the nebula depends on its state of evolution at the time of radioactivity injection. In one scenario, a collection of low-mass stars, including our sun, formed in a group or cluster with an intermediate- to high-mass star that ended its life as a supernova while our sun was still a protostar, a starless core, or perhaps a diffuse cloud. Using recent observations of protostars to estimate the size of the protosolar nebula constrains the distance of the supernova at 0.02 to 1.6 pc. The supernova distance limit is consistent with the scales of low-mass stars formation around one or more massive stars, but it is closer than expected were the sun formed in an isolated, solitary state. Consequently, if any presolar radioactivities originated via supernova injection, we must conclude that our sun was a member of such a group or cluster that has since dispersed, and thus that solar system formation should be understood in this context. In addition, we show that the timescale from explosion to the creation of small bodies was on the order of 1.8 Myr (formal 90% confidence range of 0 to 2.2 Myr), and thus the temporal choreography from supernova ejecta to meteorites is important. Finally, we can not distinguish between progenitor masses from 15 to 25 solar masses in the nucleosynthesis models; however, the 20 solar mass model is somewhat preferred.

Cosmic-ray interactions with interstellar gas produces both 6Li, which accumulates in the interstellar medium (ISM), and $\pi^0$ mesons, which decay to gamma-rays which propagate throughout the cosmos. Local 6Li abundances and extragalactic gamma-rays thus have a common origin which tightly links them. We exploit this connection to use gamma-ray observations to infer the contribution to 6Li nucleosynthesis by standard Galactic cosmic-ray (GCR) interactions with the ISM. Our calculation uses a carefully propagated cosmic-ray spectrum and accounts for 6Li production from both fusion reactions ($\alpha \alpha \to ^6Li$) as well as from spallation channels (${p,\alpha+CNO \to ^6Li$). We find that although extreme assumptions yield a consistent picture, more realistic ones indicate that solar 6Li cannot be produced by standard GCRs alone without overproducing the hadronic gamma rays. Implications for the primordial 6Li production by decaying dark matter and cosmic rays from cosmological structure formation are discussed. Upcoming gamma-ray observations by GLAST will be crucial for determining the resolution of this problem.

Pion decay gamma rays have long been recognized as a unique signature of hadronic cosmic rays and their interactions with the interstellar medium. We present a model-independent way of constraining this signal with observations of the Galactic Plane in diffuse gamma rays. We combine detections by the EGRET instrument at GeV energies and the Milagro Cherenkov detector at TeV energies with upper limits from KASCADE and CASA-MIA ground arrays at PeV energies. Such a long "lever arm", spanning at least six orders of magnitude in energy, reveals a "TeV excess" in the diffuse Galactic Plane gamma-ray spectrum. While the origin of this excess is unknown, it likely implies also enhanced TeV neutrino fluxes, significantly improving the prospects for their detection. We show that unresolved point sources are a possible source of the TeV excess. In fact, the spectra of the unidentified EGRET sources in the Milagro region must break between ~10 GeV and ~1 TeV to avoid strongly overshooting the Milagro measurement; this may have important implications for cosmic-ray acceleration. Finally, we use our approach to examine the recent suggestion that dark-matter annihilation may account for the observed excess in diffuse Galactic gamma-rays detected by EGRET at energies above 1 GeV. Within our model-independent approach, current data cannot rule this possibility in or out; however we point out how a long "lever arm" can be used to constrain the pionic gamma-ray component and in turn limit the "GeV excess" and its possible sources. Experiments such as HESS and MAGIC, and the upcoming VERITAS and GLAST, should be able to finally disentangle the main sources of the Galactic gamma rays.

The pre-Galactic abundance of lithium offers a unique window into non-thermal cosmological processes. The primordial Li abundance is guaranteed to be present and probes big bang nucleosynthesis (BBN), while an additional Li component is likely to have been produced by cosmic rays accelerated in large scale structure formation. Pre-Galactic Li currently can only be observed in low metallicity Galactic halo stars, but abundance measurements are plagued with systematic uncertainties due to modeling of stellar atmospheres and convection. We propose a new site for measuring pre-Galactic Li: low-metallicity, high-velocity clouds (HVCs) which are likely to be extragalactic gas accreted onto the Milky Way, and which already have been found to have deuterium abundances consistent with primordial. A Li observation in such an HVC would provide the first extragalactic Li measurement, and could shed new light on the apparent discrepancy between BBN predictions and halo star Li abundance determinations. Furthermore, HVC Li could at the same time test for the presence of non-primordial Li due to cosmic rays. The observability of elemental and isotopic Li abundances is discussed, and candidate sites identified.

Big bang nucleosynthesis (BBN) and the cosmic baryon density from cosmic microwave background anisotropies together predict a primordial Li7 abundance a factor of 2--3 higher than that observed in galactic halo dwarf stars. A recent analysis of Li7 observations in halo stars, using significantly higher surface temperature for these stars, found a higher Li plateau abundance. These results go a long way towards resolving the discrepancy with BBN. Here, we examine the implications of the higher surface temperatures on the abundances of Be and B which are thought to have been produced in galactic cosmic-ray nucleosynthesis by spallation of CNO together with Li (produced in alpha + alpha collisions). While the Be abundance is not overly sensitive to the surface temperature, the derived B abundances and more importantly the derived oxygen abundances are very temperature dependent. If the new temperature scale is correct, the implied increased abundances of these elements poses a serious challenge to models of galactic cosmic ray nucleosynthesis and galactic chemical evolution.

Live 60Fe has recently been detected in a deep-ocean ferromanganese crust, isolated in layers dating from about 3 Myr ago. Since 60Fe has a mean life of 2.2 Myr, a near-Earth supernova is the only likely source for such a signal, and we explore here the consequences of a supernova origin. We combine the 60Fe data with several supernova nucleosynthesis models to calculate the supernova distance as a function of progenitor mass, finding an allowed range of 15-120 pc. We also predict the signals expected for several other radioisotopes, which are independent of the supernova distance. Species likely to be present near or above background levels are 10Be, 26Al, 53Mn, 182Hf and 244Pu. Of these, 182Hf and 244Pu are nearly background-free, presenting the best opportunities to provide strong confirmation of the supernova origin of the 60Fe signal, and to demonstrate that at least some supernovae are the source for the r-process. The accuracies of our predictions are hampered by large uncertainties in the predicted 60Fe yields for supernovae of different masses, so the new crust data motivate a redoubled theoretical attack on this problem.

We present a double distribution function of dark matter halos, with respect to both object mass and local over- (or under-) density. This analytical tool provides a statistical treatment of the properties of matter surrounding collapsed objects, and can be used to study environmental effects on hierarchical structure formation. The size of the "local environment" of a collapsed object is defined to depend on the mass of the object. The Press-Schechter mass function is recovered by integration of our double distribution over the density contrast. We also present a detailed treatment of the evolution of overdensities and underdensities in Einstein-deSitter and flat LCDM universes, according to the spherical evolution model. We explicitly distinguish between true and linearly extrapolated overdensities and provide conversion relations between the two quantities.

A recent analysis of the He4 abundance determined from observations of extragalactic HII regions indicates a significantly greater uncertainty for the He4 mass fraction. The derived value is now in line with predictions from big bang nucleosynthesis when the baryon density determined by WMAP is assumed. Based on this new analysis of He4, we derive constraints on a host of particle properties which include: limits on the number of relativistic species at the time of BBN (commonly taken to be the limit on neutrino flavors), limits on the variations of fundamental couplings such as alpha_{em} and G_N, and limits on decaying particles.

The rare isotope 6Li is made only by cosmic rays, dominantly in alpha+alpha fusion reactions with ISM helium. Consequently, this nuclide provides a unique diagnostic of the history of cosmic rays in our Galaxy. The same hadronic cosmic-ray interactions also produce high-energy gamma rays (mostly via neutral pion production). Thus, hadronic gamma-rays and 6Li are intimately linked. Specifically, 6Li directly encodes the local cosmic-ray fluence over cosmic time, while extragalactic hadronic gamma rays encode an average cosmic-ray fluence over lines of sight out to the horizon. We examine this link and show how 6Li and gamma-rays can be used together to place important model-independent limits on the cosmic-ray history of our Galaxy and the universe. We first constrain gamma-ray production from ordinary Galactic cosmic rays, using the local 6Li abundance. We find that the solar 6Li abundance demands an accompanying extragalactic pionic gamma-ray intensity which exceeds that of the entire observed EGRB by a factor of 2-6. Possible explanations for this discrepancy are discussed. We then constrain Li production using recent determinations of extragalactic gamma-ray background (EGRB). We note that cosmic rays created during cosmic structure formation would lead to pre-Galactic Li production, which would act as a "contaminant" to the primordial 7Li content of metal-poor halo stars. We find the uncertainties in the observed EGRB are so large that we cannot exclude a pre-Galactic Li which is comparable to primordial 7Li. Our limits and their more model-dependent extensions will improve significantly with additional observations of 6Li in halo stars, and with improved measurements of the EGRB spectrum by GLAST. (Abriged abstract)

A critical review is given of the current status of cosmological nucleosynthesis. In the framework of the standard model with 3 types of relativistic neutrinos, the baryon-to-photon ratio, \eta, corresponding to the inferred primordial abundances of helium-4 and lithium-7 is presently ~2 \sigma below the value implied by the abundance of deuterium. The latter value is also coincident with the independent determination of \eta from WMAP observations of CMB anisotropy. However taking systematic uncertainties in the abundance estimates into account, there is overall concordance in the range \eta = (3.4 - 6.9) x 10^{-10} @ 95% c.l. corresponding to a cosmological baryon density \Omega_B h^2 = 0.012 - 0.025. If the above discrepancy is due to a neutrino chemical potential, then upto 7.1 effective neutrino species are allowed by nucleosynthesis. Other constraints on new physics are briefly discussed.

Decays of radionuclides throughout the Earth's interior produce geothermal heat, but also are a source of antineutrinos. The (angle-integrated) geoneutrino flux places an integral constraint on the terrestrial radionuclide distribution. In this paper, we calculate the angular distribution of geoneutrinos, which opens a window on the differential radionuclide distribution. We develop the general formalism for the neutrino angular distribution, and we present the inverse transformation which recovers the terrestrial radioisotope distribution given a measurement of the neutrino angular distribution. Thus, geoneutrinos not only allow a means to image the Earth's interior, but offering a direct measure of the radioactive Earth, both (1) revealing the Earth's inner structure as probed by radionuclides, and (2) allowing for a complete determination of the radioactive heat generation as a function of radius. We present the geoneutrino angular distribution for the favored Earth model which has been used to calculate geoneutrino flux. In this model the neutrino generation is dominated by decays in the Earth's mantle and crust; this leads to a very ``peripheral'' angular distribution, in which 2/3 of the neutrinos come from angles > 60 degrees away from the downward vertical. We note the possibility of that the Earth's core contains potassium; different geophysical predictions lead to strongly varying, and hence distinguishable, central intensities (< 30 degrees from the downward vertical). Other uncertainties in the models, and prospects for observation of the geoneutrino angular distribution, are briefly discussed. We conclude by urging the development and construction of antineutrino experiments with angular sensitivity. (Abstract abridged.)

Diffuse gamma rays probe the highest-energy processes at the largest scales. Here we derive model-independent constraints on the hadronic contribution to the Galactic and extragalactic gamma-ray spectra in the energy range 50 MeV < E_gamma < 10 GeV. The hadronic component is dominated by emission from neutral pions, with a characteristic spectrum symmetric about m_{pi^0}/2. We exploit the well-defined properties of the pion decay spectrum to quantify the maximum pionic fraction of the observed gamma-ray intensity. We find that the Galactic spectrum above 30 MeV can be at most about 50% pionic. The maximum pionic contribution to the extragalactic spectrum is energy dependent; it also depends on the redshift range over which the sources are distributed, ranging from as low as about 20% for pions generated very recently, to as much as 90% if the pions are generated around redshift 10. The implications of these constraints for models of gamma-ray and neutrino emission are briefly discussed.

Using the recent WMAP determination of the baryon-to-photon ratio, 10^{10} \eta = 6.14 to within a few percent, big bang nucleosynthesis (BBN) calculations can make relatively accurate predictions of the abundances of the light element isotopes which can be tested against observational abundance determinations. At this value of \eta, the Li7 abundance is predicted to be significantly higher than that observed in low metallicity halo dwarf stars. Among the possible resolutions to this discrepancy are 1) Li7 depletion in the atmosphere of stars; 2) systematic errors originating from the choice of stellar parameters - most notably the surface temperature; and 3) systematic errors in the nuclear cross sections used in the nucleosynthesis calculations. Here, we explore the last possibility, and focus on possible systematic errors in the He3(\alpha,\gamma)Be7 reaction, which is the only important Li7 production channel in BBN. The absolute value of the cross section for this key reaction is known relatively poorly both experimentally and theoretically. The agreement between the standard solar model and solar neutrino data thus provides additional constraints on variations in the cross section (S_{34}). Using the standard solar model of Bahcall, and recent solar neutrino data, we can exclude systematic S_{34} variations of the magnitude needed to resolve the BBN Li7 problem at > 95% CL. Additional laboratory data on He3(\alpha,\gamma)Be7 will sharpen our understanding of both BBN and solar neutrinos, particularly if care is taken in determining the absolute cross section and its uncertainties. Nevertheless, it already seems that this ``nuclear fix'' to the Li7 BBN problem is unlikely; other possible solutions are briefly discussed.

Big bang nucleosynthesis has long provided the primary determination of the cosmic baryon density $\omb h^2$, or equivalently the baryon-to-photon ratio, \eta. Recently, data on CMB anisotropies have become increasingly sensitive to \eta. The comparison of these two independent measures provides a key test for big bang cosmology. The first release of results from the Wilkinson Microwave Anisotropy Probe (WMAP) marks a milestone in this test. With the precision of WMAP, the CMB now offers a significantly stronger constraint on \eta. We discuss the current state of BBN theory and light element observations (including their possible lingering systematic errors). The resulting BBN baryon density prediction is in overall agreement with the WMAP prediction, an important and non-trivial confirmation of hot big bang cosmology. Going beyond this, the powerful CMB baryometer can be used as an input to BBN and one can accurately predict the primordial light element abundances. By comparing these with observations one can obtain new insight into post-BBN nucleosynthesis processes and associated astrophysics. Finally, one can test the possibility of nonstandard physics at the time of BBN, now with all light elements available as probes. Indeed, with the WMAP precision \eta, deuterium is already beginning to rival \he4's sensitivity to nonstandard physics, and additional D/H measurements can improve this further.

Deuterium plays a crucial role in cosmology because the primordial D/H abundance, in the context of big bang nucleosynthesis (BBN) theory, yields a precise measure of the cosmic baryon content. Observations of D/H can limit or measure the true primordial abundance because D is thought to be destroyed by stars and thus D/H monotonically decreases after BBN. Recently, however, Mullan & Linsky have pointed out that D arises as a secondary product of neutrons in stellar flares which then capture on protons via n+p \to d + gamma, and that this could dominate over direct D production in flares. Mullan & Linsky note that if this process is sufficiently vigorous in flaring dwarf stars, it could lead to significant non-BBN D production. We have considered the production of D in stellar flares, both directly and by n capture. We find that for reasonable flare spectra, n/d < 10 and (n+d)/6Li < 400, both of which indicate that the n-capture channel does not allow for Galactic D production at a level which will reverse the monotonic decline of D. We also calculate the 2.22 MeV gamma-ray line production associated with n capture, and find that existing COMPTEL limits also rule out significant D production in the Galaxy today. Thus, we find flares in particular, and neutron captures in general, are not an important Galactic source of D. On the other hand, we cannot exclude that flare production might contribute to recent FUSE observations of large variations in the local interstellar D/H abundance; we do, however, offer important constraints on this possibility. Finally, since flare stars should inevitably produce some n-capture events, a search for diffuse 2.22 MeV gamma-rays by INTEGRAL can further constrain (or measure!) Galactic deuterium production via n-capture.

We revisit the upper limits on the abundance of unstable massive relic particles provided by the success of Big-Bang Nucleosynthesis calculations. We use the cosmic microwave background data to constrain the baryon-to-photon ratio, and incorporate an extensively updated compilation of cross sections into a new calculation of the network of reactions induced by electromagnetic showers that create and destroy the light elements deuterium, he3, he4, li6 and li7. We derive analytic approximations that complement and check the full numerical calculations. Considerations of the abundances of he4 and li6 exclude exceptional regions of parameter space that would otherwise have been permitted by deuterium alone. We illustrate our results by applying them to massive gravitinos. If they weigh ~100 GeV, their primordial abundance should have been below about 10^{-13} of the total entropy. This would imply an upper limit on the reheating temperature of a few times 10^7 GeV, which could be a potential difficulty for some models of inflation. We discuss possible ways of evading this problem.

Recent observations by Bania et al. (2002) measure He3 versus oxygen in Galactic HII regions, finding that He3/H is within a factor of 2 of the solar abundance for [O/H] > -0.6. These results are consistent with a flat behavior in this metallicity range, tempting one to deduce from these observations a primordial value for the He3 abundance, which could join D and Li7 as an indicator of the cosmic baryon density. However, using the same data, we show that it is not possible to obtain a strong constraint on the baryon density range. This is due to (i) the intrinsically weak sensitivity of the primordial He3 abundance to the baryon density; (ii) the limited range in metallicity of the sample; (iii) the intrinsic scatter in the data; and (iv) our limited understanding of the chemical and stellar evolution of this isotope. Consequently, the He3 observations correspond to an extended range of baryon-to-photon ratio, eta = (2.2 - 6.5) x 10^{-10}, which diminishes the role of He3 as a precision baryometer. On the other hand, once the baryon-to-photon ratio is determined by the CMB, D/H, or Li7/H, the primordial value of He3/H can be inferred. Henceforth new observations of Galactic He3, can in principle greatly improve our understanding of stellar and/or chemical evolution and reconcile the observations of the HII regions and those of the planetary nebulae.

The diffuse extragalactic gamma-ray background (EGRB) above 100 MeV encodes unique information about high-energy processes in the universe. Numerous sources for the EGRB have been proposed, but the two systems which are certain to make some contribution are active galaxies (blazars) as well as normal galaxies. In this paper, we evaluate the contribution to the background from both sources. The active galaxy contribution arises from unresolved blazars. We compute this contribution using the Stecker-Salamon model. For normal galaxies, the emission is due to cosmic-ray interactions with diffuse gas. Our key assumption here is that the cosmic-ray flux in a galaxy is proportional to the supernova rate and thus the massive star formation rate, quantified observationally by the cosmic star formation rate (CSFR). In addition, the existence of stars today requires a considerably higher ISM mass in the past. Using the CSFR to compute both these effects, we find that normal galaxies are responsible for a significant portion (\sim 1/3) of the EGRB near 1 GeV, but make a smaller contribution at other energies. Finally, we present a "minimal" two-component model which includes contributions from both normal galaxies and blazars. We show that the spectrum of the diffuse radiation is a key constraint on this model: while neither the blazar spectra, nor the galactic spectra, are separately optimal fits to the observed spectrum, the combined emission provides an excellent fit. We close by noting key observational tests of this two-component model, which can be probed by future gamma-ray observatories such as GLAST.

We present images of the supernova remnant N206 in the LMC, taken with ATCA at wavelengths of 3 and 6 cm. Based on our data and previously published flux densities, the spectral index of N206 is -0.20 +- 0.07. The 6-cm radio morphology shows a filled center. Most interesting is the discovery of a peculiar linear feature previously undetected at any wavelength. The feature lies to the east of the center of the remnant, stretching from about one-fourth to three-fourths of the remnant's radius. It is wedge-shaped, with a steady opening angle from an apex on the eastern side. The feature resembles the disturbance expected from an object moving through the material supersonically at about 800 km/s. We present arguments suggesting that the linear feature might have been produced by a low-mass star or compact object ejected from a binary system which may have led to a Type Ic supernova.

Self-interacting dark matter (SIDM) was introduced by Spergel & Steinhardt to address possible discrepancies between collisionless dark matter simulations and observations on scales of less than 1 Mpc. We examine the case in which dark matter particles not only have strong self-interactions but also have strong interactions with baryons. The presence of such interactions will have direct implications for nuclear and particle astrophysics. Among these are a change in the predicted abundances from big bang nucleosynthesis (BBN) and the flux of gamma-rays produced by the decay of neutral pions which originate in collisions between dark matter and Galactic cosmic rays (CR). From these effects we constrain the strength of the baryon--dark matter interactions through the ratio of baryon - dark matter interaction cross section to dark matter mass, $s$. We find that BBN places a weak upper limit to this ratio $< 10^8 cm^2/g$. CR-SIDM interactions, however, limit the possible DM-baryon cross section to $< 5 \times 10^{-3} cm^2/g$; this rules out an energy-independent interaction, but not one which falls with center-of-mass velocity as $s \propto 1/v$ or steeper.

We investigate a possible nucleosynthetic signature of highly energetic explosions of C-O cores ("hypernovae," HNe) which might be associated with gamma-ray bursts (GRBs). We note that the direct impact of C- and O-enriched hypernova ejecta on the ambient hydrogen and helium leads to spallation reactions which can produce large amounts of the light nuclides lithium, beryllium, and boron (LiBeB). Using analytic velocity spectra of the hypernova ejecta, we calculate the LiBeB yields of different exploding C-O cores associated with observed hypernovae. The deduced yields are $\sim 10^3$ times higher than those produced by similar (direct) means in normal Type II supernovae, and are higher than the commonly used ones arising from shock wave acceleration induced by Type II supernova (SN) explosions. To avoid overproduction of these elements in our Galaxy, hypernovae should be rare events, with $\la 10^{-3}$ hypernovae per supernova, assuming a constant HN/SN ratio over time. This rate is in good agreement with that of long duration GRBs if we assume that the gamma-ray emission is focussed with a beaming factor $\Omega/4\pi \la 10^{-2}$. This encouraging result supports the possible HN-GRB association. Thus, Galactic LiBeB abundance measurements offer a promising way to probe the HN rate history and the possible HN-GRB correlation. On the other hand, if hypernovae are associated to very massive pregalactic stars (Population III) they would produce a LiBeB pre-enrichment in proto-galactic gas, which could show up as a plateau in the lowest metallicities of the Be-Fe relation in halo stars.

Recent observations suggest the existence of a white dwarf population in the Galactic halo, while others suggest that deuterium has been astrated in systems at high redshift and low metallicity. We propose that these observations could be signatures of an early population of intermediate-mass stars. Such a population requires a Population III initial mass function different from that of the solar neighborhood, as perhaps also suggested by the observed cosmic infrared background. Also, to avoid overproduction of C and N, it is required that the Z=0 yields of these stars have low ($\sim 10^{-3}$ solar) abundances as suggested by some recent calculations. Under these assumptions, we present a model which reproduces the observed D vs Si trend, and predicts a high cosmic Type Ia supernova rate, while producing a white dwarf population that accounts for only $\sim$ 1.5% of the dark halo. This scenario can be tested by observations of the cosmic supernova rate, and by confirmation and further studies of the putative white dwarf halo population.

Cosmic microwave background (CMB) determinations of the baryon-to-photon ratio $\eta \propto \Omega_{\rm baryon} h^2$ will remove the last free parameter from (standard) big bang nucleosynthesis (BBN) calculations. This will make BBN a much sharper probe of early universe physics, for example, greatly refining the BBN measurement of the effective number of light neutrino species, $N_{\nu,eff}$. We show how the CMB can improve this limit, given current light element data. Moreover, it will become possible to constrain $N_{\nu,eff}$ independent of \he4, by using other elements, notably deuterium; this will allow for sharper limits and tests of systematics. For example, a 3% measurement of $\eta$, together with a 10% (3%) measurement of primordial D/H, can measure $N_{\nu,eff}$ to a 95% confidence level of $\sigma_{95%}(N_\nu) = 1.8$ (1.0) if $\eta \sim 6.0\times 10^{-10}$. If instead, one adopts the standard model value $N_{\nu,eff}=3$, then one can use $\eta$ (and its uncertainty) from the CMB to make accurate predictions for the primordial abundances. These determinations can in turn become key inputs in the nucleosynthesis history (chemical evolution) of galaxies thereby placing constraints on such models.

The theoretical predictions of big bang nucleosynthesis (BBN) are dominated by uncertainties in the input nuclear reaction cross sections. In this paper, we examine the impact on BBN of the recent compilation of nuclear data and thermonuclear reactions rates by the NACRE collaboration. We confirm that the adopted rates do not make large overall changes in central values of predictions, but do affect the magnitude of the uncertainties in these predictions. Therefore, we then examine in detail the uncertainties in the individual reaction rates considered by NACRE. When the error estimates by NACRE are treated as 1\sigma limits, the resulting BBN error budget is similar to those of previous tabulations. We propose two new procedures for deriving reaction rate uncertainties from the nuclear data: one which sets lower limits to the error, and one which we believe is a reasonable description of the present error budget. We propagate these uncertainty estimates through the BBN code, and find that when the nuclear data errors are described most accurately, the resulting light element uncertainties are notably smaller than in some previous tabulations, but larger than others. Using these results, we derive limits on the cosmic baryon-to-photon ratio $\eta$, and compare this to independent limits on $\eta$ from recent balloon-borne measurements of the cosmic microwave background radiation (CMB). We discuss means to improve the BBN results via key nuclear reaction measurements and light element observations.

Diffuse gamma-ray radiation in galaxies is produced by cosmic ray interactions with the interstellar medium. With the completion of EGRET observations, the only extragalactic object from which there has been a positive detection of diffuse gamma-ray emission is the Large Magellanic Cloud. We systematically estimate the expected diffuse gamma-ray flux from Local Group galaxies, and determine their detectability by new generation gamma-ray observatories such as GLAST. For each galaxy, the expected gamma-ray flux depends only on its total gas content and its cosmic ray flux. We present a method for calculating cosmic ray flux in these galaxies in terms of the observed rate of supernova explosions, where cosmic ray acceleration is believed to take place. The difficulty in deriving accurate supernova rates from observational data is a dominant uncertainty in our calculations. We estimate the gamma-ray flux for Local Group galaxies and find that our predictions are consistent with the observations for the LMC and with the observational upper limits for the Small Magellanic Cloud and M31. Both the Andromeda galaxy, with a flux of $\sim 1.0 \times 10^{-8}$ photons sec$^{-1}$ cm$^{-2}$ above 100 MeV, and the SMC, with a flux of $\sim 1.7 \times 10^{-8}$ photons sec$^{-1}$ cm$^{-2}$ above 100 MeV, are expected to be observable by GLAST. M33 is at the limit of detectability with a flux of $\sim 0.11 \times 10^{-8}$ sec$^{-1}$ cm$^{-2}$. Other Local Group galaxies are at least two orders of magnitude below GLAST sensitivity.

Recent observations of heavy elements produced by rapid neutron capture (r-process) in the halo have shown a striking and unexpected behavior: within a single star, the relative abundances of r-process elements heavier than Eu are the same as the same as those of solar system matter, while across stars with similar metallicity Fe/H, the r/Fe ratio varies over two orders of magnitude. In this paper we present a simple analytic model which describes a star's abundances in terms of its ``ancestry,'' i.e., the number of nucleosynthesis events (e.g., supernova explosions) which contributed to the star's composition. This model leads to a very simple analytic expression for the abundance scatter versus Fe/H, which is in good agreement with the data and with more sophisticated numerical models. We investigate two classes of scenarios for r-process nucleosynthesis, one in which r-process synthesis events occur in only \sim 4% of supernovae but iron synthesis is ubiquitous, and one in which iron nucleosynthesis occurs in only about 9% of supernovae. (the Wasserburg- Qian model). We find that the predictions in these scenarios are similar for [Fe/H] \ga -2.5, but that these models can be readily distinguished observationally by measuring the dispersion in r/Fe at [Fe/H] \la -3.

We review recent observational studies of heavy element abundances in low metallicity stars and explore some implications of these results for nucleosynthesis and early Galactic chemical evolution.

The recent observations of Be and B in metal poor stars has led to a reassessment of the origin of the light elements in the early Galaxy. At low it is metallicity ([O/H] < -1.75), it is necessary to introduce a production mechanism which is independent of the interstellar metallicity (primary). At higher metallicities, existing data might indicate that secondary production is dominant. In this paper, we focus on the secondary process, related to the standard Galactic cosmic rays, and we examine the cosmic ray energy requirements for both present and past epochs. We find the power input to maintain the present-day Galactic cosmic ray flux is about 1.5e41 erg/s = 5e50 erg/century. This implies that, if supernovae are the sites of cosmic ray acceleration, the fraction of explosion energy going to accelerated particles is about 30%, a value which we obtain consistently both from considering the present cosmic ray flux and confinement and from the present 9Be and 6Li abundances. Using the abundances of 9Be (and 6Li) in metal-poor halo stars, we extend the analysis to show the effect of the interstellar gas mass on the standard galactic cosmic ray energetic constraints on models of Li, Be, and B evolution. The efficiency of the beryllium production per erg may be enhanced in the past by a factor of about 10; thus the energetic requirement by itself cannot be used to rule out a secondary origin of light elements. Only a clear and undisputable observational determination of the O-Fe relation in the halo will discriminate between the two processes. (abridged)

The nature of the dark matter in the Universe is one of the outstanding questions in astrophysics. In this talk, I address possible stellar baryonic contributions to the 50-90% of our Galaxy that is made of unknown dark matter. First I show that faint stars and brown dwarfs constitute only a few percent of the mass of the Galaxy. Next, I show that stellar remnants, including white dwarfs and neutron stars, are also insufficient in abundance to explain all the dark matter of the Galaxy. High energy gamma-rays observed in HEGRA data place the most robust constraints, $\Omega_{WD} < 3 \times 10^{-3} h^{-1}$, where $h$ is the Hubble constant in units of 100 km s$^{-1}$ Mpc$^{-1}$. Overproduction of chemical abundances (carbon, nitrogen, and helium) provide the most stringent constraints, $\Omega_{WD} < 2 \times 10^{-4} h^{-1}$. Comparison with recent updates of microlensing data are also made. According to the gamma-ray limit, all Massive Compact Halo Objects seen by the experiments (Machos) can be white dwarfs if one takes the extreme numbers; however, from chemical overproduction limits, NOT all Machos can be white dwarfs. Comments on recent observations of the infrared background and of white dwarfs are also made. In conclusion, a nonbaryonic component in the Halo seems to be required.

Project GRAND measures the asymmetries of secondary cosmic muons by calculating the differences between the muons from the east and the muons from the west for complete solar days. A similar north-south difference is also presented. These muons are identified in single track data taken with the 64 stations of Project GRAND; the data cover a two year span with a total muon sample of 30 billion muons.

The nucleosynthesis of Be and B by spallation processes provides unique insight into the origin of cosmic rays. Namely, different spallation schemes predict sharply different trends for the growth of LiBeB abundances with respect to oxygen. ``Primary'' mechanisms predict BeB $\propto$ O, and are well motivated by the data if O/Fe is constant at low metallicity. In contrast, ``secondary'' mechanisms predict BeB $\propto$ O$^2$ and are consistent with the data if O/Fe increases towards low metallicity as some recent data suggest. Clearly, any primary mechanism, if operative, will dominate early in the history of the Galaxy. In this paper, we fit the BeB data to a two-component scheme which includes both primary and secondary trends. In this way, the data can be used to probe the period in which primary mechanisms are effective. We analyze the data using consistent stellar atmospheric parameters based on Balmer line data and the continuum infrared flux. Results depend sensitively on Pop II O abundances and, unfortunately, on the choice of stellar parameters. When using recent results which show O/Fe increasing toward lower metallicity, a two-component Be-O fits indicates that primary and secondary components contribute equally at [O/H]$_{eq}$ = -1.8 for Balmer line data; and [O/H]$_{eq}$ = -1.4 to -1.8 for IRFM. We apply these constraints to recent models for LiBeB origin. The Balmer line data does not show any evidence for primary production. On the other hand, the IRFM data does indicate a preference for a two-component model, such as a combination of standard GCR and metal-enriched particles accelerated in superbubbles. These conclusions rely on a detailed understanding of the abundance data including systematic effects which may alter the derived O-Fe and BeB-Fe relations.

In the standard hot big bang nucleosynthesis (BBN) model the primordial abundances of H, H2, He3, He4, and Li7, fix the baryon density of the universe, $\Omega_b$, via the baryon-to-photon ratio, $\eta$, for a given Hubble parameter. Recent observations of Li show that its intrinsic dispersion in metal-poor stars is essentially zero, and the random error in the mean Li abundance is negligible. However, a decreasing trend in the Li abundance towards lower metallicity, plus Li6 detections, indicate that its primordial abundance can be inferred only after allowing for nucleosynthesis processes in the Galaxy (Galactic chemical evolution (GCE)). We show that the Li vs Fe trend provides a tough discriminant between alternative models for GCE of light-elements. We critically assess current systematic uncertainties, and determine the primordial Li abundance within new, much tighter limits. We show that the Li constraint on $\Omega_b$ is now limited principally by uncertainties in the nuclear cross-sections used in BBN {\it calculations}, not by the abundance itself. A clearer understanding of systematics allows for a much more accurate inference of the primordial Li abundance, sharpening the comparison with He4 and deuterium and the resulting test of BBN. We find that the Li data are in good agreement with He4 and with ``high'' deuterium values, but that low deuterium abundances are at best marginally within the Li range.

The nature of the dark matter in the Halo of our Galaxy remains a mystery. Arguments are presented that the dark matter does not consist of ordinary stellar or substellar objects, i.e., the dark matter is not made of faint stars, brown dwarfs, white dwarfs, or neutron stars. In fact, faint stars and brown dwarfs constitute no more than a few percent of the mass of our Galaxy, and stellar remnants must satisfy $\Omega_{WD} \leq 3 \times 10^{-3} h^{-1}$, where $h$ is the Hubble constant in units of 100 km/s Mpc^{-1}. On theoretical grounds one is then pushed to more exotic explanations. Indeed a nonbaryonic component in the Halo seems to be required.

We examine the chemical abundance constraints on a population of white dwarfs in the Halo of our Galaxy. We are motivated by microlensing evidence for massive compact halo objects (Machos) in the Galactic Halo, but our work constrains white dwarfs in the Halo regardless of what the Machos are. We focus on the composition of the material that would be ejected as the white dwarfs are formed; abundance patterns in the ejecta strongly constrain white dwarf production scenarios. Using both analytical and numerical chemical evolution models, we confirm that very strong constraints come from Galactic Pop II and extragalactic carbon abundances. We also point out that depending on the stellar model, significant nitrogen is produced rather than carbon. The combined constraints from C and N give $\Omega_{WD} h < 2 \times 10^{-4}$ from comparison with the low C and N abundances in the Ly$\alpha$ forest. We note, however, that these results are subject to uncertainties regarding the nucleosynthesis of low-metallicity stars. We thus investigate additional constraints from D and $^4$He, finding that these light elements can be kept within observational limits only for $\Omega_{WD} \la 0.003$ and for a white dwarf progenitor initial mass function sharply peaked at low mass (2$M_\odot$). Finally, we consider a Galactic wind, which is required to remove the ejecta accompanying white dwarf production from the galaxy. We show that such a wind can be driven by Type Ia supernovae arising from the white dwarfs themselves, but these supernovae also lead to unacceptably large abundances of iron. We conclude that abundance constraints exclude white dwarfs as Machos. (abridged)

The galactic chemical evolution of Be and B provides unique information about the origin and history of cosmic rays. The available Pop II data demonstrate that Be and B have a Galactic source, probably in one or more kinds of spallation processes. However, the data are not unequivocal about the nature of Be and B origin, as encoded in the primary or secondary (linear or quadratic) scaling with metallicity. We summarize a careful analysis of the trends among Be, B, Fe, and O observations. We show that if O/Fe is constant, some other cosmic ray origin or component is needed. On the other hand, if O/Fe is not constant, as recent data suggest, then the data could indicate a standard cosmic ray origin, wherein the abundances of cosmic rays scale with those of the ISM. We suggest future observational tests which will distinguish several proposed scenarios of LiBeB and cosmic ray origin.

The dual origin of population II Li7, in both big bang nucleosynthesis and galactic cosmic-ray nucleosynthesis is discussed. It is argued that with additional Li6 data, stringent limits on the degree of Li7 depletion can be obtained. Li7 depletion is also constrained by the concordance of big bang predictions with observational determinations of light element abundances. Stringent limits can also be obtained for a fixed primordial D/H abundance.

The nature of the Massive Compact Halo objects seen in microlensing experiments and interpreted as dark matter in the Halo of our Galaxy remains a mystery. Arguments are presented that these events are probably not ordinary stellar or substellar objects, i.e., they are probably not faint stars, brown dwarfs, white dwarfs, or neutron stars. On theoretical grounds one is then pushed to either exotic explanations or a "no-Macho" Halo (in which the Machos reside elsewhere than in the Halo). Indeed a nonbaryonic component in the Halo seems to be required.

Live $^{60}$Fe has recently been reported in a deep-ocean ferromanganese crust. Analysis of the isotopic ratios in the sample suggests that the measured $^{60}$Fe abundance exceeds the levels generated by terrestrial and cosmogenic sources, and it has been proposed that the excess of $^{60}$Fe is a signature of a supernova that exploded near the earth several Myr ago. In this paper, we consider the possible background sources, and confirm that the measured $^{60}$Fe is significantly higher than all known backgrounds, in contrast with the reported abundance of live $^{53}$Mn. We discuss scenarios in which the data are consistent with a supernova event at a distance $D \sim 30$ pc and an epoch $t_{\rm SN} \sim 5$ Myr ago. We propose tests that could confirm or refute the interpretation of the $^{60}$Fe discovery, including searches for $^{10}$Be, $^{129}$I and $^{146}$Sm. Such a nearby supernova event might have had some impact on the earth's biosphere, principally by enhancing the cosmic-ray flux. This might have damaged the earth's ozone layer, enhancing the penetration of solar ultraviolet radiation. In this connection, we comment on the Middle Miocene and Pliocene mini-extinction events. We also speculate on the possibility of a supernova-induced "cosmic-ray winter," if cosmic rays play a significant role in seeding cloud formation.

We review the Galactic chemical evolution of Li6 and compare these results with recent observational determinations of the lithium isotopic ratio. In particular, we concentrate on so-called standard Galactic cosmic-ray nucleosynthesis in which Li, Be, and B are produced (predominantly) by the inelastic scattering of accelerated protons and \alpha's off of CNO nuclei in the ambient interstellar medium. If O/Fe is constant at low metallicities, then the Li6 vs Fe/H evolution-as well as Be and B vs Fe/H-has difficulty in matching the observations. However, recent determinations of Population II oxygen abundances, as measured via OH lines, indicate that O/Fe increases at lower metallicity; if this trend is confirmed, then the Li6 evolution in a standard model of cosmic-ray nucleosynthesis is consistent with the data. We also show that another key indicator of Li6BeB origin is the Li6/Be ratio which also fits the available data if O/Fe is not constant at low metallicity. Finally we note that Li6 evolution in this scenario can strongly constrain the degree to which Li6 and Li7 are depleted in halo stars.

Because of the roughly linear correlation between Be/H and Fe/H in low metallicity halo stars, it has been argued that a ``primary'' component in the nucleosynthesis of Be must be present in addition to the ``secondary'' component from standard Galactic cosmic ray nucleosynthesis. In this paper we critically re-evaluate the evidence for the primary versus secondary character of Li, Be, and B evolution, analyzing both in the observations and in Galactic chemical evolution models. While it appears that [Be/H] versus [Fe/H] has a logarithmic slope near 1, it is rather the Be-O trend that directly arises from the physics of spallation production. Using new abundances for oxygen in halo stars based on UV OH lines, we find that the Be-O slope has a large uncertainty due to systematic effects, rendering it difficult to distinguish from the data between the secondary slope of 2 and the primary slope of 1. The possible difference between the Be-Fe and Be-O slopes is a consequence of the variation in O/Fe versus Fe: recent data suggests a negative slope rather than zero (i.e., Fe $\propto$ O) as is often assumed. In addition to a phenomenological analysis of Be and B evolution, we have also examined the predicted LiBeB, O, and Fe trends in Galactic chemical evolution models which include outflow. Based on our results, it is possible that a good fit to the LiBeB evolution requires only traditional the Galactic cosmic ray spallation, and the (primary) neutrino-process contribution to B11. We thus suggest that these two processes might be sufficient to explain Li6, Be, and B evolution in the Galaxy, without the need for an additional primary source of Be and B.

Data in the 2090 A B region of HD 76932 have been obtained at high S/N using the HST GHRS echelle at a resolution of 90,000. This wavelength region has been previously identified as a likely candidate for observing the B11/B10 isotopic splitting. The observations do not match a calculated line profile extremely well at any abundance for any isotopic ratio. If the B abundance previously determined from observations at 2500 A is assumed, the calculated line profile is too weak, indicating a possible blending line. Assuming that the absorption at 2090 A is entirely due to boron, the best-fit total B abundance is higher than but consistent with that obtained at 2500 A, and the best-fit isotopic ratio (B11/B10) is in the range ~10:1 to ~4:1. If the absorption is not entirely due to B and there is an unknown blend, the best-fit isotopic ratio may be closer to 1:1. Future observations of a similar metal-poor star known to have unusually low B should allow us to distinguish between these two possibilities. The constraints that can be placed on the isotopic ratio based on comparisons with similar observations of HD 102870 and HD 61421 (Procyon) are also discussed.

Recently, there has been great progress toward observationally determining the mean star formation history of the universe. When accurately known, the cosmic star formation rate could provide much information about Galactic evolution, if the Milky Way's star formation rate is representative of the average cosmic star formation history. A simple hypothesis is that our local star formation rate is proportional to the cosmic mean. In addition, to specify a star formation history, one must also adopt an initial mass function (IMF); typically it is assumed that the IMF is a smooth function which is constant in time. We show how to test directly the compatibility of all these assumptions, by making use of the local (solar neighborhood) star formation record encoded in the present-day stellar mass function. Present data suggests that at least one of the following is false: (1) the local IMF is constant in time; (2) the local IMF is a smooth (unimodal) function; and/or (3) star formation in the Galactic disk was representative of the cosmic mean. We briefly discuss how to determine which of these assumptions fail, and improvements in observations which will sharpen this test.

[Abridged] We estimate the contribution of Massive Compact Halo Objects (Machos) and their stellar progenitors to the mass density of the Universe. If the Machos that have been detected reside in the Halo of our Galaxy, then a simple extrapolation of the Galactic population (out to 50 kpc) of Machos to cosmic scales gives a cosmic density $\rho_{Macho} = (1-5) \times 10^9 h \msun \Mpc^{-3}$, which in terms of the critical density corresponds to $\Omega_{Macho}=(0.0036-0.017) h^{-1}$. Such a mass density is comparable to the baryon density implied by Big Bang Nucleosynthesis. If we take the central values of the estimates, then Machos dominate the baryonic content of the Universe today, with $\Omega_{Macho}/\Omega_{Baryon} \sim 0.7 h$. However, the cumulative uncertainties in the density determinations only require that $\Omega_{Macho}/\Omega_{Baryon} \geq 1/6 h f_{gal}$, where the fraction of galaxies that contain Machos $f_{gal} > 0.17$, and $h$ is the Hubble constant in units of 100 km s$^{-1}$ Mpc$^{-1}$. Our best estimate for $\Omega_{Macho}$ is hard to reconcile with the current best estimates of the baryonic content of the intergalactic medium indicated by measurements of the Lyman-$\alpha$ forest. We explore the addition constraints that arise if the Machos are white dwarfs as suggested by the present microlensing data. We discuss the challenges this scenario presents at both the local and cosmic scales, emphasizing in particular the constraints on the required mass budget and nucleosynthesis products (particularly carbon).

We discuss the chemical evolution of dwarf irregular and blue compact galaxies in light of recent data, new stellar yields and chemical evolution models. We examine the abundance data for evidence of HII region self-enrichment effects, which would lead to correlations in the scatter of helium, nitrogen, and oxygen abundances around their mean trends. The observed helium abundance trends show no such correlations, though the nitrogen--oxygen trend does show strong evidence for real scatter beyond observational error. We construct simple models for the chemical evolution of these galaxies, using the most recent yields of \he4, C, N and O in intermediate- and high-mass stars. The effects of galactic outflows, which can arise both from bulk heating and evaporation of the ISM, and from the partial escape of enriched supernova ejecta are included. In agreement with other studies, we find that supernova-enriched outflows can roughly reproduce the observed He, C, N, and O trends; however, in models that fit N versus O, the slopes $\Delta Y/\Delta$O and $\Delta Y/\Delta$N consistently fall more than $2\sigma$ below the fit to observations. We discuss the role of the models and their uncertainties in the extrapolation of primordial helium from the data. We also explore the model dependence arising nucleosynthesis uncertainties associated with nitrogen yields in intermediate mass stars, the fate of $8-11 \msol$ stars, and massive star winds.

We present a schematic model for the formation of baryonic galactic halos and hot gas in the Local Group and the intergalactic medium. We follow the dynamics, chemical evolution, heat flow and gas flows of a hierarchy of scales, including: protogalactic clouds, galactic halos, and the Local Group itself. Within this hierarchy, the Galaxy is built via mergers of protogalactic fragments. We find that early bursts of star formation lead to a large population of remnants (mostly white dwarfs), which would reside presently in the halo and contribute to the dark component observed in the microlensing experiments. The hot, metal-rich gas from early starbursts and merging evaporates from the clouds and is eventually incorporated into the intergalactic medium. The model thus suggests that most microlensing objects could be white dwarfs ($m \sim 0.5 \msol$), which comprise a significant fraction of the halo mass. Furthermore, the Local Group could have a component of metal-rich hot gas similar to, although less than, that observed in larger clusters. We discuss the known constraints on such a scenario and show that all local observations can be satisfied with present data in this model. The best-fit model has a halo that is 40% baryonic, with an upper limit of 77%.

Nearby supernova explosions may cause geological isotope anomalies via the direct deposition of debris or by cosmic-ray spallation in the earth's atmosphere. We estimate the mass of material deposited terrestrially by these two mechanisms, showing the dependence on the supernova distance. A number of radioactive isotopes are identified as possible diagnostic tools, such as Be-10, Al-26, Cl-36, Mn-53, Fe-60, and Ni-59, as well as the longer-lived I-129, Sm-146, and Pu-244. We discuss whether the 35 and 60 kyr-old Be-10 anomalies observed in the Vostok antarctic ice cores could be due to supernova explosions. Combining our estimates for matter deposition with results of recent nucleosynthesis yields, we calculate the expected signal from nearby supernovae using ice cores back to $\sim 300$ kyr ago, and we discuss using deep ocean sediments back to several hundred Myr. In particular, we examine the prospects for identifying isotope anomalies due to the Geminga supernova explosion, and signatures of the possibility that supernovae might have caused one or more biological mass extinctions.

We examine in detail how BBN theory is constrained, and what predictions it can make, when using only the most model-independent observational constraints. We avoid the uncertainties and model-dependencies that necessarily arise when solar neighborhood D and \he3 abundances are used to infer primordial D and \he3 via chemical and stellar evolution models. Instead, we use \he4 and \li7, thoroughly examining the effects of possible systematic errors in each. Via a likelihood analysis, we find near perfect agreement between BBN theory and the most model-independent data. Given this agreement, we then {\it assume} the correctness of BBN to set limits on the single parameter of standard BBN, the baryon-to-photon ratio, and to predict the primordial D and \he3 abundances. We also repeat our analysis including recent measurements of D/H from quasar absorption systems and find that the near perfect agreement between theory and observation of the three isotopes, D, \he4 and \li7 is maintained. These results have strong implications for the chemical and stellar evolution of the light elements, in particular for \he3. In addition, our results (especially if the D/H measurements are confirmed) have implications for the stellar depletion of \li7. Finally, we set limits on the number \nnu\ of neutrino flavors, using an analysis which carefully and systematically includes all available experimental constraints. The value \nnu = 3.0 fits best with BBN and a 95\% CL upper limit of \nnu $\la 4$ is established.

The fastest ejecta of supernova explosions propagate as a precursor to the main supernova shock wave, and can be quite energetic. The spectrum of such fast ejecta is estimated based on recent analytic and numerical supernova models, and found to be a power law having a cutoff at an energy of order 10 MeV/nucl, the precise value of which depends on the supernova mass and energetics. With cutoffs in this range there can be significant flux with energies above the thresholds for $\gamma$-ray and Li, Be, and B production. These nuclear interactions should inevitably accompany the passage of prompt supernova ejecta through the surrounding medium. The fast particle composition is that of the outermost layers of the progenitor; if the progenitor experienced significant loss of its envelope to a companion or by mass loss, then the composition is nonsolar, and in particular, metal-rich. Such a composition and spectrum of fast particles is required to explain the recent COMPTEL observations of $\gamma$-ray emission in Orion. Supernova ejecta from one progenitor having lost a large amount of mass are shown to quantitatively account for the Orion $\gamma$ rays, and to imply the presence of other, weaker lines and of a supernova explosion in Orion in the last $\la 10^4$ years. Implications for this mechanism as a mode of nucleosynthesis, and as a potential supernova diagnostic, are discussed. Model dependences and uncertainties are noted, and the need is shown for accurate measurements of $\gamma$-ray and spallation production cross sections at energies near threshold.

We compute the nucleosynthesis bounds on the masses of stable Dirac and Majorana neutrinos by solving an evolution equation network comprising of all neutrino species which in the Dirac case includes different helicity states as separate species. We will not commit ourselves to any particular value of the nucleosynthesis bound on effective number of light neutrino degrees of freedom $N_\nu$, but present all our mass bounds as {\em functions} of $\dNu$. For example, we find that the excluded region in the mass of a Majorana $\mu$- or $\tau$- neutrino, $0.31$ MeV $< m_{\nu}^M < 52$ MeV corresponding to a bound $\dNu < 0.3$ gets relaxed to $0.93$ MeV $< m_{\nu}^M < 31$ MeV if $\dNu < 1.0$ is used instead. For the Dirac neutrinos this latter constraint gives the upper limits (for $T_{\rm QCD} = 100$ MeV): $\mnm < 0.31$ MeV and $\mnt < 0.37$ MeV. Also, the constraint $\dNu <1$ {\it allows} a stable Dirac neutrino with $m^D_{\nu_\tau} > 22$ MeV.

The light elements D, \he3, \he4, and \li7 are produced in big bang nucleosynthesis and undergo changes in their abundances due to galactic processing. Since one may observe most of these elements only in contemporary environments, knowledge of the intervening evolution is necessary for determining the observational constraints on primordial nucleosynthesis. Chemical and stellar evolution model dependences in light element evolution are systematically investigated via a comparison of 1460 possible chemical evolution scenarios and of stellar nucleosynthesis yields, all of which have been selected to fit solar neighborhood C, N, O, and Fe abundances as well as the observed local gas density and gas mass fraction. The light element evolution and solar system yields in these models are found to span a wide range, explicitly demonstrating the model dependence. The range of model dependence for D, \he3, and \he4 solar abundances is calculated, and its sensitivity to the heavy element constraints is noted. The chemical evolution contribution to the uncertainty in the observed primordial light element abundances is estimated, and the effects of this uncertainty on big bang nucleosynthesis results are discussed. The predictions for the light elements are found to be correlated; the extent and physical origin of these correlations is discussed. D and \he3 evolution is found to have significant model dependence, however, the dominant factor determining their solar and interstellar abundances is their primordial abundance. In addition, \he3 is found to be very sensitive to the details of processing in low mass stars. \he4 yields are shown to be very model dependent; in particular, both the introduction of mass loss and the possibly very high \he4 stellar yields in the

The production of LiBeB isotopes by nuclear and neutrino spallation are compared in the framework of galactic evolutionary models. As motivated by $\gamma$-ray observations of Orion, different possible sources of low-energy C and O nuclei are considered, such as supernovae of various masses and WC stars. We confirm that the low energy nuclei (LEN), injected in molecular clouds by stellar winds and type II supernovae originating from the most massive progenitors, can very naturally reproduce the observed Be and B evolution in the early galaxy (halo phase). Assuming the global importance of the LEN component, we compute upper and lower bounds to the neutrino process contribution corresponding to limiting cases of LEN particle spectra. A consistent solution is found with a spectrum of the kind proposed by Ramaty \etal (1995a,b), e.g. flat up to $E_c=30$ MeV/n and decreasing abruptly above. This solution fulfills the challenge of explaining at the same time the general Be and B evolution, and their solar system abundances without overproducing \li7 at very low metallicities, and the meteoritic \b11/\b10 ratio. In this case, neutrino spallation is constrained to play a limited role in the genesis of the solar system \b11. Galactic cosmic rays (GCR) become operative late in the evolution of the disk ([Fe/H]$>$-1), but their contribution to the solar abundances of \be9, \b10 and \b11 is not dominant (35\%, 30\% and 20\% respectively). Thus, with this LEN spectrum, GCR are {\it not}\ the main source of \be9 and B in the Galaxy. The most favorable case for neutrinos, (adopting the same kind of spectrum) has $E_c=20$ MeV/n. Even in this case, the neutrino yields of Woosley and Weaver (1995) must to be reduced by a factor of 5 to avoid \b11 overproduction. Furthermore, this solution leads to a high B/Be

We derive constraints from standard (with $N_\nu = 3$) BBN arising solely from two cosmologically produced nuclides, \he4 and \li7, from which the extrapolation to primordial abundances is straightforward. The abundances of D and \he3 are at present only inferred from their solar and local interstellar medium values using models of galactic chemical evolution. However, our knowledge of chemical evolution suffers from large uncertainties, and so it is of use to take an approach which minimizes the reliance on chemical evolution in determining the consistency of BBN. Using only the data on \he4 and \li7, and reasonable estimates of possible systematic errors in their abundance determinations, concordance is found if the baryon-to-photon ratio lies in the range $1.4 < 10^{10} \eta < 3.8$ (95\% CL). The \he4 and \li7 abundances are also used to predict the primordial abundances D and \he3, which provides an initial condition for their chemical evolution.

The observed B/Be ratio in extreme Pop II stars has been interpreted as evidence of Be and B synthesis by early galactic cosmic rays. However, a recent reanalysis of the boron abundance in the Pop II halo star HD140283 suggests that B/H may be larger than previously reported, by as much as a factor of 4. This would yield a B/Be ratio lying in the range $14 \la {\rm B/Be} \la 50$. The possibility of a high Pop II B/Be ratio stresses the importance of the upper limit to the B/Be ratio arising from cosmic ray production. It is found that the limit to cosmic ray-produced B/Be depends upon the assumed cosmic ray spectrum. For any Pop II comic ray spectrum that is a single power law in either total energy per nucleon or in momentum the B/Be ratio is constrained to lie in the range $7.6 \la$ B/Be $\la 14$. Thus, if the new B/Be ratio is correct, it requires either a bimodal cosmic ray flux with a large low energy component, or for another B source, possibly the proposed $\nu$-process in supernovae, either of which may be helpful in explaining the observed \bor11/\bor10 ratio. Finally, it is noted that the boron reanalysis highlights the uncertainty in our knowledge of the B/Be ratio, and the need for additional data on Be and B abundances.

To understand better the early galactic production of Li, Be, and B by cosmic ray spallation and fusion reactions, the dependence of these production rates on cosmic ray models and model parameters is examined. The sensitivity of elemental and isotopic production to the cosmic ray pathlength magnitude and energy dependence, source spectrum, spallation kinematics, and cross section uncertainties is studied. Changes in these model features, particularly those features related to confinement, are shown to alter the Be- and B-versus-Fe slopes {}from a na\"{\i}ve quadratic relation. The implications of our results for the diffuse $\gamma$-ray background are examined, and the role of chemical evolution and its relation to our results is noted. It is also noted that the unmeasured high energy behavior of $\alpha+\alpha$ fusion can lead to effects as large as a factor of 2 in the resultant yields. Future data should enable Population II Li, Be, and B abundances to constrain cosmic ray models for the early Galaxy.

We investigate the possibility that inhomogeneous nucleosynthesis may eventually be used to explain the abundances of \li6, \be9 and B in population II stars. The present work differs from previous studies in that we have used a more extensive reaction network. It is demonstrated that in the simplest scenario the abundances of the light elements with $A\le7$ constrain the separation of inhomogeneities to sufficiently small scales that the model is indistinguishable from homogeneous nucleosynthesis and that the abundances of \li6, \be9\ and B are then below observations by several orders of magnitude. This conclusion does not depend on the \li7 constraint. We also examine alternative scenarios which involve a post-nucleosynthesis reprocessing of the light elements to reproduce the observed abundances of Li and B, while allowing for a somewhat higher baryon density (still well below the cosmological critical density). Future B/H measurements may be able to exclude even this exotic scenario and further restrict primordial nucleosynthesis to approach the homogeneous model conclusions.

We discuss in detail the effect of neutrino oscillations in Big Bang nucleosynthesis, between active and sterile neutrinos, as well as between active and active neutrinos. We calculate the constraints on mixings between active and sterile neutrinos from the present observation of the primordial helium abundance and discuss the potential implications on various astrophysical and cosmological problems of such oscillations. In particular, we show that large angle sterile neutrino mixing seems to be excluded as a MSW solution to the solar neutrino situation or a solution to the atmospheric neutrino mixing hinted at in some underground experiments. We show how with this constraint, the next generation of solar neutrino experiments should be able to determine the resolution of the solar neutrino problem. It is also shown how sterile neutrinos remain a viable dark matter candidate.

We have modified the standard code for primordial nucleosynthesis to include the effect of the slight heating of neutrinos by $e^\pm$ annihilations. There is a small, systematic change in the $^4$He yield, $\Delta Y \simeq +1.5\times 10^{-4}$, which is insensitive to the value of the baryon-to-photon ratio $\eta$ for $10^{-10}\la \eta \la 10^{-9}$. We also find that the baryon-to-photon ratio decreases by about 0.5\% less than the canonical factor of 4/11 because some of the entropy in $e^\pm$ pairs is transferred to neutrinos. These results are in accord with recent analytical estimates.

The ability to now make measurements of Be and B as well as put constraints on \lisix\ abundances in metal-poor stars has led to a detailed reexamination of Big Bang Nucleosynthesis in the $A\groughly6$ regime. The nuclear reaction network has been significantly expanded with many new rates added. It is demonstrated that although a number of $A>7$ reaction rates are poorly determined, even with extreme values chosen, the standard homogeneous model is unable to produce significant yields (Be/H and B/H $<10^{-17}$ when $A\le7$ abundances fit) above $A=7$ and the \liseven/\lisix\ ratio always exceeds 500. We also preliminarily explore inhomogeneous models, such as those inspired by a first order quark-hadron phase transition, where regions with high neutron/proton ratios can allow some leakage up to $A>7$. However models that fit the $A\le7$ abundances still seem to have difficulty in obtaining significant $A>7$ yields.