The age of the universe from nuclear chronometers

An overview is presented of the current situation regarding radioactive dating of the matter of which our Galaxy is comprised. A firm lower bound on the age from nuclear chronometers of ≈9–10 Gyr is entirely consistent with age determinations from globular clusters and white dwarf cooling histories. The reasonable assumption of an approximately uniform nucleosynthesis rate yields an age for the Galaxy of 12.8 ± 3 Gyr, which again is consistent with current determinations from other methods.

X-ray emission from clusters and groups of galaxies

Recent major advances in x-ray imaging and spectroscopy of clusters have allowed the determination of their mass and mass profile out to ≈1/2 the virial radius. In rich clusters, most of the baryonic mass is in the gas phase, and the ratio of mass in gas/stars varies by a factor of 2–4. The baryonic fractions vary by a factor of ≈3 from cluster to cluster and almost always exceed 0.09 h50−[3/2] and thus are in fundamental conflict with the assumption of Ω = 1 and the results of big bang nucleosynthesis. The derived Fe abundances are 0.2–0.45 solar, and the abundances of O and Si for low redshift systems are 0.6–1.0 solar. This distribution is consistent with an origin in pure type II supernova. The amount of light and energy produced by these supernovae is very large, indicating their importance in influencing the formation of clusters and galaxies. The lack of evolution of Fe to a redshift of z ≈ 0.4 argues for very early enrichment of the cluster gas. Groups show a wide range of abundances, 0.1–0.5 solar. The results of an x-ray survey indicate that the contribution of groups to the mass density of the universe is likely to be larger than 0.1 h50−2. Many of the very poor groups have large x-ray halos and are filled with small galaxies whose velocity dispersion is a good match to the x-ray temperatures.

The peculiar Type Ia supernova iPTF14atg: Chandrasekhar-mass explosion or violent merger?

iPTF14atg, a subluminous peculiar Type Ia supernova (SN Ia) similar to SN 2002es, is the first SN Ia for which a strong UV flash was observed in the early-time light curves. This has been interpreted as evidence for a single-degenerate (SD) progenitor system, where such a signal is expected from interactions between the SN ejecta and the non-degenerate companion star. Here, we compare synthetic observables of multidimensional state-of-the-art explosion models for different progenitor scenarios to the light curves and spectra of iPTF14atg. From our models, we have difficulties explaining the spectral evolution of iPTF14atg within the SD progenitor channel. In contrast, we find that a violent merger of two carbon-oxygen white dwarfs with 0.9 and 0.76 M⊙, respectively, provides an excellent match to the spectral evolution of iPTF14atg from 10 d before to several weeks after maximum light. Our merger model does not naturally explain the initial UV flash of iPTF14atg. We discuss several possibilities like interactions of the SN ejecta with the circumstellar medium and surface radioactivity from an He-ignited merger that may be able to account for the early UV emission in violent merger models.

Three-dimensional simulations of gravitationally confined detonations compared to observations of SN 1991T

The gravitationally confined detonation (GCD) model has been proposed as a possible explosion mechanism for Type Ia supernovae in the single-degenerate evolution channel. It starts with ignition of a deflagration in a single off-centre bubble in a near-Chandrasekhar-mass white dwarf. Driven by buoyancy, the deflagration flame rises in a narrow cone towards the surface. For the most part, the main component of the flow of the expanding ashes remains radial, but upon reaching the outer, low-pressure layers of the white dwarf, an additional lateral component develops. This causes the deflagration ashes to converge again at the opposite side, where the compression heats fuel and a detonation may be launched. We first performed five three-dimensional hydrodynamic simulations of the deflagration phase in 1.4 M⊙ carbon/oxygen white dwarfs at intermediate-resolution (2563computational zones). We confirm that the closer the initial deflagration is ignited to the centre, the slower the buoyant rise and the longer the deflagration ashes takes to break out and close in on the opposite pole to collide. To test the GCD explosion model, we then performed a high-resolution (5123 computational zones) simulation for a model with an ignition spot offset near the upper limit of what is still justifiable, 200 km. This high-resolution simulation met our deliberately optimistic detonation criteria, and we initiated a detonation. The detonation burned through the white dwarf and led to its complete disruption. For this model, we determined detailed nucleosynthetic yields by post-processing 106 tracer particles with a 384 nuclide reaction network, and we present multi-band light curves and time-dependent optical spectra. We find that our synthetic observables show a prominent viewing-angle sensitivity in ultraviolet and blue wavelength bands, which contradicts observed SNe Ia. The strong dependence on the viewing angle is caused by the asymmetric distribution of the deflagration ashes in the outer ejecta layers. Finally, we compared our model to SN 1991T. The overall flux level of the model is slightly too low, and the model predicts pre-maximum light spectral features due to Ca, S, and Si that are too strong. Furthermore, the model chemical abundance stratification qualitatively disagrees with recent abundance tomography results in two key areas: our model lacks low-velocity stable Fe and instead has copious amounts of high-velocity 56Ni and stable Fe. We therefore do not find good agreement of the model with SN 1991T.

Photoionization of Co+ and electron-impact excitation of Co2 + using the Dirac R-matrix method

Modelling of massive stars and supernovae (SNe) plays a crucial role in understanding galaxies. From this modelling we can derive fundamental constraints on stellar evolution, mass-loss processes, mixing, and the products of nucleosynthesis. Proper account must be taken of all important processes that populate and depopulate the levels (collisional excitation, de-excitation, ionization, recombination, photoionization, bound–bound processes). For the analysis of Type Ia SNe and core collapse SNe (Types Ib, Ic and II) Fe group elements are particularly important. Unfortunately little data is currently available and most noticeably absent are the photoionization cross-sections for the Fe-peaks which have high abundances in SNe. Important interactions for both photoionization and electron-impact excitation are calculated using the relativistic Dirac atomic R-matrix codes (DARC) for low-ionization stages of Cobalt. All results are calculated up to photon energies of 45 eV and electron energies up to 20 eV. The wavefunction representation of Co III has been generated using GRASP0 by including the dominant 3d⁷, 3d⁶[4s, 4p], 3p⁴3d⁹ and 3p⁶3d⁹ configurations, resulting in 292 fine structure levels. Electron-impact collision strengths and Maxwellian averaged effective collision strengths across a wide range of astrophysically relevant temperatures are computed for Co III. In addition, statistically weighted level-resolved ground and metastable photoionization cross-sections are presented for Co II and compared directly with existing work.

The Status of Multi-Dimensional Core-Collapse Supernova Models

Models of neutrino-driven core-collapse supernova explosions have matured considerably in recent years. Explosions of low-mass progenitors can routinely be simulated in 1D, 2D, and 3D. Nucleosynthesis calculations indicate that these supernovae could be contributors of some lighter neutron-rich elements beyond iron. The explosion mechanism of more massive stars remains under investigation, although first 3D models of neutrino-driven explosions employing multi-group neutrino transport have become available. Together with earlier 2D models and more simplified 3D simulations, these have elucidated the interplay between neutrino heating and hydrodynamic instabilities in the post-shock region that is essential for shock revival. However, some physical ingredients may still need to be added/improved before simulations can robustly explain supernova explosions over a wide range of progenitors. Solutions recently suggested in the literature include uncertainties in the neutrino rates, rotation, and seed perturbations from convective shell burning. We review the implications of 3D simulations of shell burning in supernova progenitors for the ‘perturbations-aided neutrino-driven mechanism,’ whose efficacy is illustrated by the first successful multi-group neutrino hydrodynamics simulation of an 18 solar mass progenitor with 3D initial conditions. We conclude with speculations about the impact of 3D effects on the structure of massive stars through convective boundary mixing.

The role of nuclear physics in supernovae and the evolution of neutron stars, Neutrino Opacities, Equation of State, Transport Coefficients, and Dark Matter Production

Thesis (Ph.D.)--University of Washington, 2016-08

The locations of recent supernovae near the Sun from modelling 60Fe transport

The signature of 60Fe in deep-sea crusts indicates that one or more supernovae exploded in the solar neighbourhood about 2.2 million years ago1–4. Recent isotopic analysis is consistent with a core-collapse or electron-capture supernova that occurred 60 to 130 parsecs from the Sun5. Moreover, peculiarities in the cosmic ray spectrum point to a nearby supernova about two million years ago6. The Local Bubble of hot, diffuse plasma, in which the Solar System is embedded, originated from 14 to 20 supernovae within a moving group, whose surviving members are now in the Scorpius– Centaurus stellar association7,8. Here we report calculations of the most probable trajectories and masses of the supernova progenitors, and hence their explosion times and sites. The 60Fe signal arises from two supernovae at distances between 90 and 100 parsecs. The closest occurred 2.3 million years ago at present-day galactic coordinates l = 327°, b = 11°, and the second-closest exploded about 1.5 million years ago at l = 343°, b = 25°, with masses of 9.2 and 8.8 times the solar mass, respectively. The remaining supernovae, which formed the Local Bubble, contribute to a smaller extent because they happened at larger distances and longer ago (60Fe has a half- life of 2.6 million years9,10). There are uncertainties relating to the nucleosynthesis yields and the loss of 60Fe during transport, but they do not influence the relative distribution of 60Fe in the crust layers, and therefore our model reproduces the measured relative abundances very well.