Search

Stellar mergers are one important path to highly magnetised stars. Mergers of two low-mass white dwarfs may create up to every third hot subdwarf star. The merging process is usually assumed to dramatically amplify magnetic fields. However, so far only four highly magnetised hot subdwarf stars have been found, suggesting a fraction of less than $1\%$. We present two high-resolution magnetohydrodynamical (MHD) simulations of the merger of two helium white dwarfs in a binary system with the same total mass of $0.6\,M_\odot$. We analysed an equal-mass merger with two $0.3\,M_\odot$ white dwarfs, and an unequal-mass merger with white dwarfs of $0.25\,M_\odot$ and $0.35\,M_\odot$. We simulated the inspiral, merger, and further evolution of the merger remnant for about $50$ rotations. We found efficient magnetic field amplification in both mergers via a small-scale dynamo, reproducing previous results of stellar merger simulations. The magnetic field saturates at a similar strength for both simulations. We then identified a second phase of magnetic field amplification in both merger remnants that happens on a timescale of several tens of rotational periods of the merger remnant. This phase generates a large-scale ordered azimuthal field via a large-scale dynamo driven by the magneto-rotational instability. Finally, we speculate that in the unequal-mass merger remnant, helium burning will initially start in a shell around a cold core, rather than in the centre. This forms a convection zone that coincides with the region that contains most of the magnetic energy, and likely destroys the strong, ordered field. Ohmic resistivity might then quickly erase the remaining small-scale field. Therefore, the mass ratio of the initial merger could be the selecting factor that decides if a merger remnant will stay highly magnetised long after the merger., Comment: 13 pages, 10 figures, accepted by A&A, comments welcome

We performed 3D hydrodynamic simulations of the inner $\approx 50\%$ radial extent of a $25\ \mathrm{M_\odot}$ star in the early phase of the main sequence and investigate core convection and internal gravity waves in the core-envelope boundary region. Simulations for different grid resolutions and driving luminosities establish scaling relations to constrain models of mixing for 1D applications. As in previous works, the turbulent mass entrainment rate extrapolated to nominal heating is unrealistically high ($1.58\times 10^{-4}\ \mathrm{M_\odot/yr}$), which is discussed in terms of the non-equilibrium response of the simulations to the initial stratification. We measure quantitatively the effect of mixing due to internal gravity waves excited by core convection interacting with the boundary in our simulations. The wave power spectral density as a function of frequency and wavelength agrees well with the GYRE eigenmode predictions based on the 1D spherically averaged radial profile. A diffusion coefficient profile that reproduces the spherically averaged abundance distribution evolution is determined for each simulation. Through a combination of eigenmode analysis and scaling relations it is shown that in the $N^2$-peak region, mixing is due to internal gravity waves and follows the scaling relation $D_\mathrm{IGW-hydro} \propto L^{4/3}$ over a $\gtrapprox 2\ \mathrm{dex}$ range of heating factors. Different extrapolations of the mixing efficiency down to nominal heating are discussed. If internal gravity wave mixing is due to thermally-enhanced shear mixing, an upper limit is $D_\mathrm{IGW} \lessapprox 2$ to $3\times 10^4\ \mathrm{cm^2/s}$ at nominal heating in the $N^2$-peak region above the convective core., Comment: MNRAS. Clarifications and improvements following referee recommendations. Movies and data access https://www.ppmstar.org. 38 Figs. In v2 the bottom panel of Fig 28 is showing the wrong case (M116 instead of M114). Corrected in v3

We present two mixing models for post-processing of 3D hydrodynamic simulations applied to convective-reactive i-process nucleosynthesis in a rapidly accreting white dwarf (RAWD) with $\mathrm{[Fe/H]} = -2.6$, in which H is ingested into a convective He shell. A 1D advective two-stream model adopts physically motivated radial and horizontal mixing coefficients constrained by 3D hydrodynamic simulations. A simpler approach uses diffusion coefficients calculated from the same simulations. All 3D simulations include the energy feedback of the $^{12}$C(p,$\gamma$)$^{13}$N reaction from the H entrainment. Global oscillations of shell H ingestion in two of the RAWD simulations cause bursts of entrainment of H and non-radial hydrodynamic feedback. With the same nuclear network as in the 3D simulations, the 1D advective two-stream model reproduces the rate and location of the H burning within the He shell closely matching the 3D simulation predictions, as well as qualitatively displaying the asymmetry of the $X_{\mathrm{H}}$ profiles between the up- and downstream. With a full i-process network the advective mixing model captures the difference in the n-capture nucleosynthesis in the up- and downstream. For example, $^{89}$Kr and $^{90}$Kr with half-lives of 3.18 min and 32.3 s differ by a factor 2-10 in the two streams. In this particular application the diffusion approach provides globally the same abundance distribution as the advective two-stream mixing model. The resulting i-process yields are in excellent agreement with observations of the exemplary CEMP-r/s star CS31062-050., Comment: 18 pages, 17 figures, 2 tables, accepted for publication in MNRAS

The special computational challenges of simulating 3-D hydrodynamics in deep stellar interiors are discussed, and numerical algorithmic responses described. Results of recent simulations carried out at scale on the NSF's Blue Waters machine at the University of Illinois are presented, with a special focus on the computational challenges they address. Prospects for future work using GPU-accelerated nodes such as those on the DoE's new Summit machine at Oak Ridge National Laboratory are described, with a focus on numerical algorithmic accommodations that we believe will be necessary., Comment: 14 pages, 7 figures

We report on our ongoing investigation into the nucleosynthetic and hydrodynamic nature of mixing at the interface between the H- and He-convection zones in massive Pop III stars. Studying a grid of 26 1D stellar evolution simulations with different mixing assumptions, we find that H-He interactions occur in 23/26 cases. We demonstrate the nucleosynthesis expected in a H-He interaction in an 80M$_\odot$. Finally, we describe our progress in simulating a Pop III double convection zone in the PPMStar hydrodynamics code., Comment: NIC 2018 Proceeding

We have modelled the multicycle evolution of rapidly-accreting CO white dwarfs (RAWDs) with stable H burning intermittent with strong He-shell flashes on their surfaces for $0.7\leq M_\mathrm{RAWD}/M_\odot\leq 0.75$ and [Fe/H] ranging from $0$ to $-2.6$. We have also computed the i-process nucleosynthesis yields for these models. The i process occurs when convection driven by the He-shell flash ingests protons from the accreted H-rich surface layer, which results in maximum neutron densities $N_\mathrm{n,max}\approx 10^{13}$-$10^{15}\ \mathrm{cm}^{-3}$. The H-ingestion rate and the convective boundary mixing (CBM) parameter $f_\mathrm{top}$ adopted in the one-dimensional nucleosynthesis and stellar evolution models are constrained through 3D hydrodynamic simulations. The mass ingestion rate and, for the first time, the scaling laws for the CBM parameter $f_\mathrm{top}$ have been determined from 3D hydrodynamic simulations. We confirm our previous result that the high-metallicity RAWDs have a low mass retention efficiency ($\eta < 10\%$). A new result is that RAWDs with [Fe/H]$< -2$ have $\eta > 20\%$, therefore their masses may reach the Chandrasekhar limit and they may eventually explode as SNeIa. This result and the good fits of the i-process yields from the metal-poor RAWDs to the observed chemical composition of the CEMP-r/s stars suggest that some of the present-day CEMP-r/s stars could be former distant members of triple systems, orbiting close binary systems with RAWDs that may have later exploded as SNeIa., Comment: 15 pages, 12 figures, 2 tables, accepted for publication in MNRAS

Interactions between convective shells in evolved massive stars have been linked to supernova impostors, to the production of the odd-Z elements Cl, K, and Sc, and they might also help generate the large-scale asphericities that are known to facilitate shock revival in supernova explosion models. We investigate the process of ingestion of C-shell material into a convective O-burning shell, including the hydrodynamic feedback from the nuclear burning of the ingested material. Our 3D hydrodynamic simulations span almost 3 dex in the total luminosity $L_\mathrm{tot}$. All but one of the simulations reach a quasi-stationary state with the entrainment rate and convective velocity proportional to $L_\mathrm{tot}$ and $L_\mathrm{tot}^{1/3}$, respectively. Carbon burning provides $14\,$--$\,33\%$ of the total luminosity, depending on the set of reactions considered. Equivalent simulations done on $768^3$ and $1152^3$ grids are in excellent quantitative agreement. The flow is dominated by a few large-scale convective cells. An instability leading to large-scale oscillations with Mach numbers in excess of $0.2$ develops in an experimental run with the energy yield from C burning increased by a factor of 10. This run represents most closely the conditions expected in a violent O-C shell merger, which is a potential production site for odd-Z elements such as K and Sc and which may seed asymmetries in the supernova progenitor. 1D simulations may underestimate the energy generation from the burning of ingested material by as much as a factor two owing to their missing the effect of clumpiness of entrained material on the nuclear reaction rate., Comment: 24 pages, 27 figures, accepted for publication in MNRAS

Collaborations in astronomy and astrophysics are faced with numerous cyber infrastructure challenges, such as large data sets, the need to combine heterogeneous data sets, and the challenge to effectively collaborate on those large, heterogeneous data sets with significant processing requirements and complex science software tools. The cyberhubs system is an easy-to-deploy package for small to medium-sized collaborations based on the Jupyter and Docker technology, that allows web-browser enabled, remote, interactive analytic access to shared data. It offers an initial step to address these challenges. The features and deployment steps of the system are described, as well as the requirements collection through an account of the different approaches to data structuring, handling and available analytic tools for the NuGrid and PPMstar collaborations. NuGrid is an international collaboration that creates stellar evolution and explosion physics and nucleosynthesis simulation data. The PPMstar collaboration performs large-scale 3D stellar hydrodynamics simulation of interior convection in the late phases of stellar evolution. Examples of science that is presently performed on cyberhubs, in the areas 3D stellar hydrodynamic simulations, stellar evolution and nucleosynthesis and Galactic chemical evolution, are presented., Comment: 19 pages, 9 figures, ApJ Supplement Special Issue on Data, accepted

This work investigates the properties of convection in stars with particular emphasis on entrainment across the upper convective boundary (CB). Idealised simulations of turbulent convection in the O-burning shell of a massive star are performed in $4\pi$ geometry on $768^3$ and $1536^3$ grids, driven by a representative heating rate. A heating series is also performed on the $768^3$ grid. The $1536^3$ simulation exhibits an entrainment rate at the upper CB of $1.33\times10^{-6}~M_\odot~\mathrm{s}^{-1}$. The $768^3$ simulation with the same heating rate agrees within 17 per cent. The entrainment rate at the upper convective boundary is found to scale linearly with the driving luminosity and with the cube of the shear velocity at the upper boundary, while the radial RMS fluid velocity scales with the cube root of the driving luminosity, as expected. The mixing is analysed in a 1D diffusion framework, resulting in a simple model for CB mixing. The analysis confirms previous findings that limiting the MLT mixing length to the distance to the CB in 1D simulations better represents the spherically-averaged radial velocity profiles from the 3D simulations and provides an improved determination of the reference diffusion coefficient $D_0$ for the exponential diffusion CB mixing model in 1D. From the 3D simulation data we adopt as the convective boundary the location of the maximum gradient in the horizontal velocity component which has $2\sigma$ spatial fluctuations of $\approx0.17 H_P$ . The exponentially decaying diffusion CB mixing model with $f = 0.03$ reproduces the spherically-averaged 3D abundance profiles., Comment: 21 pages, 24 figures, 1 table. Accepted for publication in MNRAS