The vertical structure and dynamics of stars in our local Galactic neighbourhood contains much information about the local distribution of visible and dark matter and of perturbations to the Milky Way disc. We use data on the positions and velocities of stars in the solar neighbourhood from \gaia\ DR2 and large spectroscopic surveys to investigate the vertical number counts and mean-velocity trend as a function of distance from the local Galactic mid-plane. We perform a detailed measurement of the wave-like North-South asymmetry in the vertical number counts, which reveals a number of deficits at heights $\approx 0.4\,\mathrm{kpc}$, $\approx 0.9\,\mathrm{kpc}$, and $\approx 1.5\,\mathrm{kpc}$, and peaks at $\approx 0.2\,\mathrm{kpc}$, $\approx 0.7\,\mathrm{kpc}$, and $\approx 1.1\,\mathrm{kpc}$. We find that the asymmetry pattern is independent of colour. The mean vertical velocity is almost constant to $<1\,\mathrm{km\,s}^{-1}$ within a few 100 pc from the mid-plane and then displays a North-South symmetric dip at $\approx0.5\,\mathrm{kpc}$ with an amplitude of $\approx 2\,\mathrm{km\,s}^{-1}$ that is a plausible velocity counterpart to the main number-count dip at a similar height. Thus, with \gaia\ DR2 we confirm at high fidelity that the local Galactic disc is undergoing a wave-like oscillation and a dynamically-consistent observational picture of the perturbed local vertical structure emerges for the first time. We also present the most precise and accurate determination of the Sun's height above the local Galactic mid-plane, correcting for any asymmetry in the vertical density: $z_\odot = 20.8 \pm 0.3\,\mathrm{pc}$.

We provide a set of stellar evolution and nucleosynthesis calculations that applies established physics assumptions simultaneously to low- and intermediate-mass and massive star models. Our goal is to provide an internally consistent and comprehensive nuclear production and yield data base for applications in areas such as pre-solar grain studies. Our non-rotating models assume convective boundary mixing where it has been adopted before. We include 8 (12) initial masses for $Z = 0.01$ ($0.02$). Models are followed either until the end of the asymptotic giant branch phase or the end of Si burning, complemented by a simple analytic core-collapse supernova models with two options for fallback and shock velocities. The explosions show which pre-supernova yields will most strongly be effected by the explosive nucleosynthesis. We discuss how these two explosion parameters impacts the light elements and the $s$ and $p$ process. For low- and intermediate-mass models our stellar yields from H to Bi include the effect of convective boundary mixing at the He-intershell boundaries and the stellar evolution feedback of the mixing process that produces the $^{13}$C pocket. All post-processing nucleosynthesis calculations use the same nuclear reaction rate network and nuclear physics input. We provide a discussion of the nuclear production across the entire mass range organized by element group. All our stellar nucleosynthesis profile and time evolution output is available electronically, and tools to explore the data on the NuGrid VOspace hosted by the Canadian Astronomical Data Centre are introduced.

High-Reynolds-number turbulence driven by stellar convection in main-sequence stars generates stochastic gravitational radiation. We calculate the wave-strain power spectral density as a function of the zero-age main-sequence mass for an individual star and for an isotropic, universal stellar population described by the Salpeter initial mass function and redshift-dependent Hopkins-Beacom star formation rate. The spectrum is a broken power law, which peaks near the turnover frequency of the largest turbulent eddies. The signal from the Sun dominates the universal background. For the Sun, the far-zone power spectral density peaks at $S(f_\mathrm{peak}) = 5.2 \times 10^{-52}$ Hz$^{-1}$ at frequency $f_\mathrm{peak} = 2.3 \times 10^{-7}$ Hz. However, at low observing frequencies $f < 3 \times 10^{-4}$ Hz, the Earth lies inside the Sun's near zone and the signal is amplified to $S_\mathrm{near}(f_\mathrm{peak}) = 4.1 \times 10^{-27}$ Hz$^{-1}$ because the wave strain scales more steeply with distance ($\propto d^{-5}$) in the near zone than in the far zone ($\propto d^{-1}$). Hence the Solar signal may prove relevant for pulsar timing arrays. Other individual sources and the universal background fall well below the projected sensitivities of the Laser Interferometer Space Antenna and next-generation pulsar timing arrays. Stellar convection sets a fundamental noise floor for more sensitive stochastic gravitational-wave experiments in the more distant future.

Classical novae are expected to contribute to the 1809-keV Galactic $\gamma$-ray emission by producing its precursor $^{26}$Al, but the yield depends on the thermonuclear rate of the unmeasured $^{25}$Al($p,\gamma$)$^{26}$Si reaction. Using the $\beta$ decay of $^{26}$P to populate the key $J^{\pi}=3^+$ resonance in this reaction, we report the first evidence for the observation of its exit channel via a $1741.6 \pm 0.6 (\textrm{stat}) \pm 0.3 (\textrm{syst})$ keV primary $\gamma$ ray, where the uncertainties are statistical and systematic, respectively. By combining the measured $\gamma$-ray energy and intensity with other experimental data on $^{26}$Si, we find the center-of-mass energy and strength of the resonance to be $E_r = 414.9 \pm 0.6(\textrm{stat}) \pm 0.3 (\textrm{syst}) \pm 0.6(\textrm{lit.})$ keV and $\omega\gamma = 23 \pm 6 (\textrm{stat})^{+11}_{-10}(\textrm{lit.})$ meV, respectively, where the last uncertainties are from adopted literature data. We use hydrodynamic nova simulations to model $^{26}$Al production showing that these measurements effectively eliminate the dominant experimental nuclear-physics uncertainty and we estimate that novae may contribute up to 30% of the Galactic $^{26}$Al.

Hydrodynamic turbulence driven by crust-core differential rotation imposes a fundamental noise floor on gravitational wave observations of neutron stars. The gravitational wave emission peaks at the Kolmogorov decoherence frequency which, for reasonable values of the crust-core shear, \Delta\Omega, occurs near the most sensitive part of the frequency band for ground-based, long-baseline interferometers. We calculate the energy density spectrum of the stochastic gravitational wave background from a cosmological population of turbulent neutron stars generalising previous calculations for individual sources. The spectrum resembles a piecewise power law, \Omega_{gw}(\nu)=\Omega_{\alpha}\nu^{\alpha}, with \alpha=-1 and 7 above and below the decoherence frequency respectively, and its normalisation scales as \Omega_{\alpha}\propto(\Delta\Omega)^{7}. Non-detection of a stochastic signal by Initial LIGO implies an upper limit on \Delta\Omega and hence by implication on the internal relaxation time-scale for the crust and core to come into co-rotation, \tau_{d}=\Delta\Omega/\dot{\Omega}, where \dot{\Omega} is the observed electromagnetic spin-down rate, with \tau_{d}\lesssim 10^{7} yr for accreting millisecond pulsars and \tau_{d}\lesssim 10^{5} yr for radio-loud pulsars. Target limits on \tau_{d} are also estimated for future detectors, namely Advanced LIGO and the Einstein Telescope, and are found to be astrophysically interesting.

We propose a new method of gravitational wave detection using a modified form of higher criticism, a statistical technique introduced by Donoho & Jin (2004). Higher criticism is designed to detect a group of sparse, weak sources, none of which are strong enough to be reliably estimated or detected individually. We apply higher criticism as a second-pass method to synthetic F-statistic and C-statistic data for a monochromatic periodic source in a binary system and quantify the improvement relative to the first-pass methods. We find that higher criticism on C-statistic data is more sensitive by ~6% than the C-statistic alone under optimal conditions (i.e. binary orbit known exactly) and the relative advantage increases as the error in the orbital parameters increases. Higher criticism is robust even when the source is not monochromatic (e.g. phase wandering in an accreting system). Applying higher criticism to a phase-wandering source over multiple time intervals gives a >30% increase in detectability with few assumptions about the frequency evolution. By contrast, in all-sky searches for unknown periodic sources, which are dominated by the brightest source, second-pass higher criticism does not provide any benefits over a first pass search

Despite much effort in the past decades, the C-burning reaction rate is uncertain by several orders of magnitude, and the relative strength between the different channels 12C(12C,alpha)20Ne, 12C(12C,p)23Na and 12C(12C,n)23Mg is poorly determined. Additionally, in C-burning conditions a high 12C+12C rate may lead to lower central C-burning temperatures and to 13C(alpha,n)16O emerging as a more dominant neutron source than 22Ne(alpha,n)25Mg, increasing significantly the s-process production. This is due to the rapid decrease of the 13N(gamma,p)12C with decreasing temperature, causing the 13C production via 13N(beta+)13C. Presented here is the impact of the 12C+12C reaction uncertainties on the s-process and on explosive p-process nucleosynthesis in massive stars, including also fast rotating massive stars at low metallicity. Using various 12C+12C rates, in particular an upper and lower rate limit of ~ 50000 higher and ~ 20 lower than the standard rate at 5*10^8 K, five 25 Msun stellar models are calculated. The enhanced s-process signature due to 13C(alpha,n)16O activation is considered, taking into account the impact of the uncertainty of all three C-burning reaction branches. Consequently, we show that the p-process abundances have an average production factor increased up to about a factor of 8 compared to the standard case, efficiently producing the elusive Mo and Ru proton-rich isotopes. We also show that an s-process being driven by 13C(alpha,n)16O is a secondary process, even though the abundance of 13C does not depend on the initial metal content. Finally, implications for the Sr-peak elements inventory in the Solar System and at low metallicity are discussed.

[Shortened] The 12C + 12C fusion reaction has been the subject of considerable experimental efforts to constrain uncertainties at temperatures relevant for stellar nucleosynthesis. In order to investigate the effect of an enhanced carbon burning rate on massive star structure and nucleosynthesis, new stellar evolution models and their yields are presented exploring the impact of three different 12C + 12C reaction rates. Non-rotating stellar models were generated using the Geneva Stellar Evolution Code and were later post-processed with the NuGrid Multi-zone Post-Processing Network tool. The enhanced rate causes core carbon burning to be ignited more promptly and at lower temperature. This reduces the neutrino losses, which increases the core carbon burning lifetime. An increased carbon burning rate also increases the upper initial mass limit for which a star exhibits a convective carbon core. Carbon shell burning is also affected, with fewer convective-shell episodes and convection zones that tend to be larger in mass. Consequently, the chance of an overlap between the ashes of carbon core burning and the following carbon shell convection zones is increased, which can cause a portion of the ashes of carbon core burning to be included in the carbon shell. Therefore, during the supernova explosion, the ejecta will be enriched by s-process nuclides synthesized from the carbon core s process. The yields were used to estimate the weak s-process component in order to compare with the solar system abundance distribution. The enhanced rate models were found to produce a significant proportion of Kr, Sr, Y, Zr, Mo, Ru, Pd and Cd in the weak component, which is primarily the signature of the carbon-core s process. Consequently, it is shown that the production of isotopes in the Kr-Sr region can be used to constrain the 12C + 12C rate using the current branching ratio for a- and p-exit channels.

Depending on mass and metallicity as well as evolutionary phase, stars occasionally experience convective-reactive nucleosynthesis episodes. We specifically investigate the situation when nucleosynthetically unprocessed, H-rich material is convectively mixed with a He-burning zone, for example in convectively unstable shell on top of electron-degenerate cores in AGB stars, young white dwarfs or X-ray bursting neutron stars. Such episodes are frequently encountered in stellar evolution models of stars of extremely low or zero metal content [...] We focus on the convective-reactive episode in the very-late thermal pulse star Sakurai's object (V4334 Sagittarii). Asplund etal. (1999) determined the abundances of 28 elements, many of which are highly non-solar, ranging from H, He and Li all the way to Ba and La, plus the C isotopic ratio. Our simulations show that the mixing evolution according to standard, one-dimensional stellar evolution models implies neutron densities in the He that are too low to obtain a significant neutron capture nucleosynthesis on the heavy elements. We have carried out 3D hydrodynamic He-shell flash convection [...] we assume that the ingestion process of H into the He-shell convection zone leads only after some delay time to a sufficient entropy barrier that splits the convection zone [...] we obtain significantly higher neutron densities (~few 10^15 1/cm^3) and reproduce the key observed abundance trends found in Sakurai's object. These include an overproduction of Rb, Sr and Y by about 2 orders of magnitude higher than the overproduction of Ba and La. Such a peculiar nucleosynthesis signature is impossible to obtain with the mixing predictions in our one-dimensional stellar evolution models. [...] We determine how our results depend on uncertainties of nuclear reaction rates, for example for the C13(\alpha, n)O16 reaction.

The contribution by massive stars (M > 9 solar masses) to the weak s-process component of the solar system abundances is primarily due to the 22Ne neutron source, which is activated near the end of helium-core burning. The residual 22Ne left over from helium-core burning is then reignited during carbon burning, initiating further s-processing that modifies the isotopic distribution. This modification is sensitive to the stellar structure and the carbon burning reaction rate. Recent work on the 12C + 12C reaction suggests that resonances located within the Gamow peak may exist, causing a strong increase in the astrophysical S-factor and consequently the reaction rate. To investigate the effect of an increased rate, 25 solar mass stellar models with three different carbon burning rates, at solar metallicity, were generated using the Geneva Stellar Evolution Code (GENEC) with nucleosynthesis post-processing calculated using the NuGrid Multi-zone Post-Processing Network code (MPPNP). The strongest rate caused carbon burning to occur in a large convective core rather than a radiative one. The presence of this large convective core leads to an overlap with the subsequent convective carbon-shell, significantly altering the initial composition of the carbon-shell. In addition, an enhanced rate causes carbon-shell burning episodes to ignite earlier in the evolution of the star, igniting the 22Ne source at lower temperatures and reducing the neutron density.

The physical mechanisms responsible for pulsar timing glitches are thought to excite quasi-normal mode oscillations in their parent neutron star that couple to gravitational wave emission. In August 2006, a timing glitch was observed in the radio emission of PSR B0833-45, the Vela pulsar. At the time of the glitch, the two co-located Hanford gravitational wave detectors of the Laser Interferometer Gravitational-wave observatory (LIGO) were operational and taking data as part of the fifth LIGO science run (S5). We present the first direct search for the gravitational wave emission associated with oscillations of the fundamental quadrupole mode excited by a pulsar timing glitch. No gravitational wave detection candidate was found. We place Bayesian 90% confidence upper limits of 6.3e-21 to 1.4e-20 on the peak intrinsic strain amplitude of gravitational wave ring-down signals, depending on which spherical harmonic mode is excited. The corresponding range of energy upper limits is 5.0e44 to 1.3e45 erg.

Nonaxisymmetric, meridional circulation inside a neutron star, excited by a glitch and persisting throughout the post-glitch relaxation phase, emits gravitational radiation. Here, it is shown that the current quadrupole contributes more strongly to the gravitational wave signal than the mass quadrupole evaluated in previous work. We calculate the signal-to-noise ratio for a coherent search and conclude that a large glitch may be detectable by second-generation interferometers like the Laser Interferometer Gravitational-Wave Observatory. It is shown that the viscosity and compressibility of bulk nuclear matter, as well as the stratification length-scale and inclination angle of the star, can be inferred from a gravitational wave detection in principle.

The slow neutron capture process in massive stars (the weak s-process) produces most of the s-only isotopes in the mass region 60 < A < 90. The nuclear reaction rates used in simulations of this process have a profound effect on the final s-process yields. We generated 1D stellar models of a 25 solar mass star varying the 12C + 12C rate by a factor of 10 and calculated full nucleosynthesis using the post-processing code PPN. Increasing or decreasing the rate by a factor of 10 affects the convective history and nucleosynthesis, and consequently the final yields.

We review some of the uncertainties in calculating nucleosynthetic yields, focusing on the explosion mechanism. Current yield calculations tend to either use a piston, energy injection, or enhancement of neutrino opacities to drive an explosion. We show that the energy injection, or more accurately, an entropy injection mechanism is best-suited to mimic our current understanding of the convection-enhanced supernova engine. The enhanced neutrino-opacity technique is in qualitative disagreement with simulations of core-collapse supernovae and will likely produce errors in the yields. But piston-driven explosions are the most discrepant. Piston-driven explosion severely underestimate the amount of fallback, leading to order-of-magnitude errors in the yields of heavy elements. To obtain yields accurate to the factor of a few level, we must use entropy or energy injection and this has become the NuGrid collaboration approach.

The "collapsar" engine for gamma-ray bursts invokes as its energy source the failure of a normal supernova and the formation of a black hole. Here we present the results of the first three-dimensional simulation of the collapse of a massive star down to a black hole, including the subsequent accretion and explosion. The explosion differs significantly from the axisymmetric scenario obtained in two-dimensional simulations; this has important consequences for the nucleosynthetic yields. We compare the nucleosynthetic yields to those of hypernovae. Calculating yields from three-dimensional explosions requires new strategies in post-process nucleosynthesis; we discuss NuGrid's plan for three-dimensional yields.

We investigate the physical conditions where 44Ti and 56Ni are created in core-collapse supernovae. In this preliminary work we use a series of post-processing network calculations with parametrized expansion profiles that are representative of the wide range of temperatures, densities and electron-to-baryon ratios found in 3D supernova simulations. Critical flows that affect the final yields of 44Ti and 56Ni are assessed.

Simulations of nucleosynthesis in astrophysical environments are at the intersection of nuclear physics reaction rate research and astrophysical applications, for example in the area of galactic chemical evolution or near-field cosmology. Unfortunately, at present the available yields for such applications are based on heterogeneous assumptions between the various contributing nuclear production sites, both in terms of modeling the thermodynamic environment itself as well as the choice of specifc nuclear reaction rates and compilations. On the other side, new nuclear reaction rate determinations are often taking a long time to be included in astrophysical applications. The NuGrid project addresses these issues by providing a set of codes and a framework in which these codes interact. In this contribution we describe the motivation, goals and first results of the NuGrid project. At the core is a new and evolving post-processing nuclesoynthesis code (PPN) that can follow quiescent and explosive nucleosynthesis following multi-zone 1D-stellar evolution as well as multi-zone hydrodynamic input, including explosions. First results are available in the areas of AGB and massive stars.

The s-process production in massive stars at very low metallicities is expected to be negligible due to the low abundance of the neutron source 22Ne, to primary neutron poisons and decreasing iron seed abundances. However, recent models of massive stars including the effects of rotation show that a strong production of 22Ne is possible in the helium core, as a consequence of the primary nitrogen production (observed in halo metal poor stars). Using the PPN post-processing code, we studied the impact of this primary 22Ne on the s process. We find a large production of s elements between strontium and barium, starting with the amount of primary 22Ne predicted by stellar models. There are several key reaction rate uncertainties influencing the s-process efficiency. Among them, 17O(alpha,gamma) may play a crucial role strongly influencing the s process efficiency, or it may play a negligible role, according to the rate used in the calculations. We also report on the development of a new parallel (MPI) post-processing code (MPPNP) designed to follow the complete nucleosynthesis in stars on highly resolved grids. We present here the first post-processing run from the ZAMS up to the end of helium burning for a 15 solar mass model.

We examine observed heavy element abundances in the Cassiopeia A supernova remnant as a constraint on the nature of the Cas A supernova. We compare bulk abundances from 1D and 3D explosion models and spatial distribution of elements in 3D models with those derived from X-ray observations. We also examine the cospatial production of 26Al with other species. We find that the most reliable indicator of the presence of 26Al in unmixed ejecta is a very low S/Si ratio (~0.05). Production of N in O/S/Si-rich regions is also indicative. The biologically important element P is produced at its highest abundance in the same regions. Proxies should be detectable in supernova ejecta with high spatial resolution multiwavelength observations.

Many nucleosynthesis and mixing processes of low-mass stars as they evolve from the Main Sequence to the thermal-pulse Asymptotic Giant Branch phase (TP-AGB) are well understood (although of course important physics components, e.g. rotation, magnetic fields, gravity wave mixing, remain poorly known). Nevertheless, in the last years presolar grain measurements with high resolution have presented new puzzling problems and strong constraints on nucleosynthesis processes in stars. The goal of the NuGrid collaboration is to present uniform yields for a large range of masses and metallicities, including low$-$mass stars and massive stars and their explosions. Here we present the first calculations of stellar evolution and high-resolution, post-processing simulations of an AGB star with an initial mass of 2 M_sun and solar-like metallicity (Z=0.01), based on the post-processing code PPN. In particular, we analyze the formation and evolution of the radiative 13C-pocket between the 17th TP and the 18th TP. The s-process nucleosynthesis profile of a sample of heavy isotopes is also discussed, before the next convective TP occurrence.

We present preliminary results from recent high-resolution double-degenerate merger simulations with the Smooth Particle Hydrodynamics (SPH) technique. We put particular emphasis on verification and validation in our effort and show the importance of details in the initial condition setup for the final outcome of the simulation. We also stress the dynamical importance of including shocks in the simulations. These results represent a first step toward a suite of simulations that will shed light on the question whether double-degenerate mergers are a viable path toward type 1a supernovae. In future simulations, we will make use of the capabilities of the NuGrid collaboration in post-processing SPH particle trajectories with a complete nuclear network to follow the detailed nuclear reactions during the dynamic merger phase.

The nucleosynthetic yield from a supernova explosion depends upon a variety of effects: progenitor evolution, explosion process, details of the nuclear network, and nuclear rates. Especially in studies of integrated stellar yields, simplifications reduce these uncertainties. But nature is much more complex, and to actually study nuclear rates, we will have to understand the full, complex set of processes involved in nucleosynthesis. Here we discuss a few of these complexities and detail how the NuGrid collaboration will address them.