All ten LIGO/Virgo binary black hole (BH-BH) coalescences reported from the O1/O2 runs have near zero effective spins. There are only three potential explanations of this fact. If the BH spin magnitudes are large then (i) either both BH spin vectors must be nearly in the orbital plane or (ii) the spin angular momenta of the BHs must be oppositely directed and similar in magnitude. Or, (iii) the BH spin magnitudes are small. We test the third hypothesis within the framework of the classical isolated binary evolution scenario of the BH-BH merger formation. We test three models of angular momentum transport in massive stars: a mildly efficient transport by meridional currents (as employed in the Geneva code), an efficient transport by the Tayler-Spruit magnetic dynamo (as implemented in the MESA code), and a very-efficient transport (as proposed by Fuller et al.) to calculate natal BH spins. We allow for binary evolution to increase the BH spins through accretion and account for the potential spin-up of stars through tidal interactions. Additionally, we update the calculations of the stellar-origin BH masses, include revisions to the history of star formation and to the chemical evolution across cosmic time. We find that we can match simultaneously the observed BH-BH merger rate density, BH masses, and effective spins. Models with efficient angular momentum transport are favored. The updated stellar-mass weighted gas-phase metallicity evolution now used in our models appears to be a key in better reproducing the LIGO/Virgo merger rate estimate. Mass losses during the pair-instability pulsation supernova phase are likely overestimated if the merger GW170729 hosts a BH more massive than 50 Msun. We also estimate rate of BH-NS mergers from recent LIGO/Virgo observations. Our updated models of BH-BH, BH-NS and NS-NS mergers are now publicly available at www.syntheticuniverse.org.

Our understanding of the engines and progenitors of gamma-ray bursts has expanded through the ages as a broader set of diagnostics has allowed us to test our understanding of these objects. Here we review the history of the growth in our understanding, focusing on 3 leading engines and 9 potential progenitors. The first gravitational wave detection of a short burst is the latest in a series of breakthrough observations shaping this understanding and we study the importance of multi-diagnostic, multi-messenger observations on these engines and their progenitors. Our understanding based on a detailed study of nearby bursts can be applied to make predictions for the trends expected as we begin to observe high redshift bursts and we discuss these trends.

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.

KIC 08626021 is a pulsating DB white dwarf of considerable recent interest, and first of its class to be extensively monitored by Kepler for its pulsation properties. Fitting the observed oscillation frequencies of KIC 08626021 to a model can yield insights into its otherwise-hidden internal structure. Template-based white dwarf models choose a luminosity profile where the luminosity is proportional to the enclosed mass, $L_r \propto M_r$, independent of the effective temperature $T_{\rm eff}$. Evolutionary models of young white dwarfs with $T_{\rm eff} \gtrsim$ 25,000 K suggest neutrino emission gives rise to luminosity profiles with $L_r$ $\not\propto$ $M_r$. We explore this contrast by comparing the oscillation frequencies between two nearly identical white dwarf models: one with an enforced $L_r \propto M_r$ luminosity profile and the other with a luminosity profile determined by the star's previous evolution history. We find the low order g-mode frequencies differ by up to $\simeq$ 70 $\mu$Hz over the range of Kepler observations for KIC 08626021. This suggests that by neglecting the proper thermal structure of the star (e.g., accounting for the effect of plasmon neutrino losses), the model frequencies calculated by using an $L_r \propto M_r$ profile may have uncorrected, effectively-random errors at the level of tens of $\mu$Hz. A mean frequency difference of 30 $\mu$Hz, based on linearly extrapolating published results, suggests a template model uncertainty in the fit precision of $\simeq$ 12% in white dwarf mass, $\simeq$ 9% in the radius, and $\simeq$ 3% in the central oxygen mass fraction.

We explore properties of core-collapse supernova progenitors with respect to the composite uncertainties in the thermonuclear reaction rates by coupling the reaction rate probability density functions provided by the STARLIB reaction rate library with $\texttt{MESA}$ stellar models. We evolve 1000 15 $M_{\odot}$ models from the pre main-sequence to core O-depletion at solar and subsolar metallicities for a total of 2000 Monte Carlo stellar models. For each stellar model, we independently and simultaneously sample 665 thermonuclear reaction rates and use them in a $\texttt{MESA}$ in situ reaction network that follows 127 isotopes from $^{1}$H to $^{64}$Zn. With this framework we survey the core mass, burning lifetime, composition, and structural properties at five different evolutionary epochs. At each epoch we measure the probability distribution function of the variations of each property and calculate Spearman Rank-Order Correlation coefficients for each sampled reaction rate to identify which reaction rate has the largest impact on the variations on each property. We find that uncertainties in $^{14}$N$(p,\gamma)^{15}$O, triple-$\alpha$, $^{12}$C$(\alpha,\gamma)^{16}$O, $^{12}$C($^{12}$C,$p$)$^{23}$Na, $^{12}$C($^{16}$O,$p$)$^{27}$Al, $^{16}$O($^{16}$O,$n$)$^{31}$S, $^{16}$O($^{16}$O,$p$)$^{31}$P, and $^{16}$O($^{16}$O,$\alpha$)$^{28}$Si reaction rates dominate the variations of the properties surveyed. We find that variations induced by uncertainties in nuclear reaction rates grow with each passing phase of evolution, and at core H-, He-depletion are of comparable magnitude to the variations induced by choices of mass resolution and network resolution. However, at core C-, Ne-, and O-depletion, the reaction rate uncertainties can dominate the variation causing uncertainty in various properties of the stellar model in the evolution towards iron core-collapse.

We explore the variation in single star 15-30 $\rm{M}_{\odot}$, non-rotating, solar metallicity, pre-supernova MESA models due to changes in the number of isotopes in a fully-coupled nuclear reaction network and adjustments in the mass resolution. Within this two-dimensional plane we quantitatively detail the range of core masses at various stages of evolution, mass locations of the main nuclear burning shells, electron fraction profiles, mass fraction profiles, burning lifetimes, stellar lifetimes, and compactness parameter at core-collapse for models with and without mass loss. Up to carbon burning we generally find mass resolution has a larger impact on the variations than the number of isotopes, while the number of isotopes plays a more significant role in determining the span of the variations for neon, oxygen and silicon burning. Choice of mass resolution dominates the variations in the structure of the intermediate convection zone and secondary convection zone during core and shell hydrogen burning respectively, where we find a minimum mass resolution of $\approx$ 0.01 $\rm{M}_{\odot}$ is necessary to achieve convergence in the helium core mass at the $\approx$5% level. On the other hand, at the onset of core-collapse we find $\approx$30% variations in the central electron fraction and mass locations of the main nuclear burning shells, a minimum of $\approx$127 isotopes is needed to attain convergence of these values at the $\approx$10% level.

We investigate properties of carbon-oxygen white dwarfs with respect to the composite uncertainties in the reaction rates using the stellar evolution toolkit, Modules for Experiments in Stellar Astrophysics (MESA) and the probability density functions in the reaction rate library STARLIB. These are the first Monte Carlo stellar evolution studies that use complete stellar models. Focusing on 3 M$_{\odot}$ models evolved from the pre main-sequence to the first thermal pulse, we survey the remnant core mass, composition, and structure properties as a function of 26 STARLIB reaction rates covering hydrogen and helium burning using a Principal Component Analysis and Spearman Rank-Order Correlation. Relative to the arithmetic mean value, we find the width of the 95\% confidence interval to be $\Delta M_{{\rm 1TP}}$ $\approx$ 0.019 M$_{\odot}$ for the core mass at the first thermal pulse, $\Delta$$t_{\rm{1TP}}$ $\approx$ 12.50 Myr for the age, $\Delta \log(T_{{\rm c}}/{\rm K}) \approx$ 0.013 for the central temperature, $\Delta \log(\rho_{{\rm c}}/{\rm g \ cm}^{-3}) \approx$ 0.060 for the central density, $\Delta Y_{\rm{e,c}} \approx$ 2.6$\times$10$^{-5}$ for the central electron fraction, $\Delta X_{\rm c}(^{22}\rm{Ne}) \approx$ 5.8$\times$10$^{-4}$, $\Delta X_{\rm c}(^{12}\rm{C}) \approx$ 0.392, and $\Delta X_{\rm c}(^{16}\rm{O}) \approx$ 0.392. Uncertainties in the experimental $^{12}$C($\alpha,\gamma)^{16}\rm{O}$, triple-$\alpha$, and $^{14}$N($p,\gamma)^{15}\rm{O}$ reaction rates dominate these variations. We also consider a grid of 1 to 6 M$_{\odot}$ models evolved from the pre main-sequence to the final white dwarf to probe the sensitivity of the initial-final mass relation to experimental uncertainties in the hydrogen and helium reaction rates.

We explore the detailed and broad properties of carbon burning in Super Asymptotic Giant Branch (SAGB) stars with 2755 MESA stellar evolution models. The location of first carbon ignition, quenching location of the carbon burning flames and flashes, angular frequency of the carbon core, and carbon core mass are studied as a function of the ZAMS mass, initial rotation rate, and mixing parameters such as convective overshoot, semiconvection, thermohaline and angular momentum transport. In general terms, we find these properties of carbon burning in SAGB models are not a strong function of the initial rotation profile, but are a sensitive function of the overshoot parameter. We quasi-analytically derive an approximate ignition density, $\rho_{ign} \approx 2.1 \times 10^6$ g cm$^{-3}$, to predict the location of first carbon ignition in models that ignite carbon off-center. We also find that overshoot moves the ZAMS mass boundaries where off-center carbon ignition occurs at a nearly uniform rate of $\Delta M_{\rm ZAMS}$/$\Delta f_{\rm{ov}}\approx$ 1.6 $M_{\odot}$. For zero overshoot, $f_{\rm{ov}}$=0.0, our models in the ZAMS mass range $\approx$ 8.9 to 11 $M_{\odot}$ show off-center carbon ignition. For canonical amounts of overshooting, $f_{\rm{ov}}$=0.016, the off-center carbon ignition range shifts to $\approx$ 7.2 to 8.8 $M_{\odot}$. Only systems with $f_{\rm{ov}}$ $\geq 0.01$ and ZAMS mass $\approx$ 7.2-8.0 $M_{\odot}$ show carbon burning is quenched a significant distance from the center. These results suggest a careful assessment of overshoot modeling approximations on claims that carbon burning quenches an appreciable distance from the center of the carbon core.

A re-analysis of data from electron-pair production following 160 A GeV 207Pb bombardment of nuclear emulsions indicates the production and decay of neutral particles of significantly lower invariant mass and shorter lifetimes than previously claimed (J. Phys. G: Nucl. Part. Phys. 34 (2007) 129-138).