No Arabic abstract
The double pulsar system PSR J0737-3039A/B is a double neutron star binary, with a 2.4-hour orbital period, which has allowed measurement of relativistic orbital perturbations to high precision. The low mass of the second-formed neutron star, as well as the low system eccentricity and proper motion, point to a different evolutionary scenario compared to most other known double neutron star systems. We describe analysis of the pulse profile shape over 6 years of observations, and present the resulting constraints on the system geometry. We find the recycled pulsar in this system, PSR J0737-3039A, to be a near-orthogonal rotator, with an average separation between its spin and magnetic axes of 90 +/- 11 +/- 5 deg. Furthermore, we find a mean 95% upper limit on the misalignment between its spin and orbital angular momentum axes of 3.2 deg, assuming that the observed emission comes from both magnetic poles. This tight constraint lends credence to the idea that the supernova that formed the second pulsar was relatively symmetric, possibly involving electron-capture onto an O-Ne-Mg core.
Most supernova explosions accompany the death of a massive star. These explosions give birth to neutron stars and black holes and eject solar masses of heavy elements. However, determining the mechanism of explosion has been a half-century journey of great complexity. In this paper, we present our perspective of the status of this theoretical quest and the physics and astrophysics upon which its resolution seems to depend. The delayed neutrino-heating mechanism is emerging as a robust solution, but there remain many issues to address, not the least of which involves the chaos of the dynamics, before victory can unambiguously be declared. It is impossible to review in detail all aspects of this multi-faceted, more-than-half-century-long theoretical quest. Rather, we here map out the major ingredients of explosion and the emerging systematics of the observables with progenitor mass, as we currently see them. Our discussion will of necessity be speculative in parts, and many of the ideas may not survive future scrutiny. Some statements may be viewed as informed predictions concerning the numerous observables that rightly exercise astronomers witnessing and diagnosing the supernova Universe. Importantly, the same explosion in the inside, by the same mechanism, can look very different in photons, depending upon the mass and radius of the star upon explosion. A 10$^{51}$-erg (one Bethe) explosion of a red supergiant with a massive hydrogen-rich envelope, a diminished hydrogen envelope, no hydrogen envelope, and, perhaps, no hydrogen envelope or helium shell all look very different, yet might have the same core and explosion evolution.
We present results from an ab initio three-dimensional, multi-physics core collapse supernova simulation for the case of a 15 M progenitor. Our simulation includes multi-frequency neutrino transport with state-of-the-art neutrino interactions in the ray-by-ray approximation, and approximate general relativity. Our model exhibits a neutrino-driven explosion. The shock radius begins an outward trajectory at approximately 275 ms after bounce, giving the first indication of a developing explosion in the model. The onset of this shock expansion is delayed relative to our two-dimensional counterpart model, which begins at approximately 200 ms after core bounce. At a time of 441 ms after bounce, the angle-averaged shock radius in our three-dimensional model has reached 751 km. Further quantitative analysis of the outcomes in this model must await further development of the post-bounce dynamics and a simulation that will extend well beyond 1 s after stellar core bounce, based on the results for the same progenitor in the context of our two-dimensional, counterpart model. This more complete analysis will determine whether or not the explosion is robust and whether or not observables such as the explosion energy, 56Ni mass, etc. are in agreement with observations. Nonetheless, the onset of explosion in our ab initio three-dimensional multi-physics model with multi-frequency neutrino transport and general relativity is encouraging.
We explore with self-consistent 2D F{sc{ornax}} simulations the dependence of the outcome of collapse on many-body corrections to neutrino-nucleon cross sections, the nucleon-nucleon bremsstrahlung rate, electron capture on heavy nuclei, pre-collapse seed perturbations, and inelastic neutrino-electron and neutrino-nucleon scattering. Importantly, proximity to criticality amplifies the role of even small changes in the neutrino-matter couplings, and such changes can together add to produce outsized effects. When close to the critical condition the cumulative result of a few small effects (including seeds) that individually have only modest consequence can convert an anemic into a robust explosion, or even a dud into a blast. Such sensitivity is not seen in one dimension and may explain the apparent heterogeneity in the outcomes of detailed simulations performed internationally. A natural conclusion is that the different groups collectively are closer to a realistic understanding of the mechanism of core-collapse supernovae than might have seemed apparent.
Close double neutron stars have been observed as Galactic radio pulsars, while their mergers have been detected as gamma-ray bursts and gravitational-wave sources. They are believed to have experienced at least one common-envelope episode during their evolution prior to double neutron star formation. In the last decades there have been numerous efforts to understand the details of the common-envelope phase, but its computational modelling remains challenging. We present and discuss the properties of the donor and the binary at the onset of the Roche-lobe overflow leading to these common-envelope episodes as predicted by rapid binary population synthesis models. These properties can be used as initial conditions for detailed simulations of the common-envelope phase. There are three distinctive populations, classified by the evolutionary stage of the donor at the moment of the onset of the Roche-lobe overflow: giant donors with fully-convective envelopes, cool donors with partially-convective envelopes, and hot donors with radiative envelopes. We also estimate that, for standard assumptions, tides would not circularise a large fraction of these systems by the onset of Roche-lobe overflow. This makes the study and understanding of eccentric mass-transferring systems relevant for double neutron star populations.
We have been working within the fundamental paradigm that core collapse supernovae (CCSNe) may be neutrino driven, since the first suggestion of this by Colgate and White nearly five decades ago. Computational models have become increasingly sophisticated, first in one spatial dimension assuming spherical symmetry, then in two spatial dimensions assuming axisymmetry, and now in three spatial dimensions with no imposed symmetries. The increase in the number of spatial dimensions has been accompanied by an increase in the physics included in the models, and an increase in the sophistication with which this physics has been modeled. Computation has played an essential role in the development of CCSN theory, not simply for the obvious reason that such multidimensional, multi-physics, nonlinear events cannot possibly be fully captured analytically, but for its role in discovery. In particular, the discovery of the standing accretion shock instability (SASI) through computation about a decade ago has impacted all simulations performed since then. Today, we appear to be at a threshold, where neutrinos, neutrino-driven convection, and the SASI, working together over time scales significantly longer than had been anticipated in the past, are able to generate explosions, and in some cases, robust explosions, in a number of axisymmetric models. But how will this play out in three dimensions? Early results from the first three-dimensional (3D), multi-physics simulation of the Oak Ridge group are promising. I will discuss the essential components of todays models and the requirements of realistic CCSN modeling, present results from our one-, two-, and three-dimensional models, place our models in context with respect to other efforts around the world, and discuss short- and long-term next steps.