Researcher profile

Zac Johnston

Zac Johnston contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 17 - UnverifiedVerification L1Unclaimed author
4works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

4 published item(s)

preprint2022arXiv

The Impact of the New $^{65\!}$As(p,$γ$)$^{66\!}$Se Reaction Rate on the Two-Proton Sequential Capture of $^{64}\!$Ge, Weak GeAs Cycles, and Type-I X-Ray Bursts such as the Clocked Burster GS 1826$-$24

We re-assess $^{65}$As(p,$γ$)$^{66}$Se reaction rates based on a set of proton thresholds of $^{66}$Se, $S_\mathrm{p}$($^{66}$Se), estimated from the experimental mirror nuclear masses, theoretical mirror displacement energies, and full $pf$-model space shell-model calculation. The self-consistent relativistic Hartree-Bogoliubov theory is employed to obtain the mirror displacement energies with much reduced uncertainty, and thus reducing the proton-threshold uncertainty up to 161 keV compared to the AME2020 evaluation. Using the simulation instantiated by the one-dimensional multi-zone hydrodynamic code, KEPLER, that closely reproduces the observed GS 1826$-$24 clocked bursts, the present forward and reverse $^{65}$As(p,$γ$)$^{66}$Se reaction rates based on a selected $S_\mathrm{p}$($^{66}$Se) = 2.469$\pm$0.054 MeV, and the latest $^{22}$Mg($α$,p)$^{25}$Al, $^{56}$Ni(p,$γ$)$^{57}$Cu(p,$γ$)$^{58}$Zn, $^{55}$Ni(p,$γ$)$^{56}$Cu, and $^{64}$Ge(p,$γ$)$^{65}$As reaction rates, we find that though the GeAs cycles is weakly established in the rapid-proton capture process path, the $^{65}$As(p,$γ$)$^{66}$Se reaction still strongly characterizes the burst tail end due to the two-proton sequential capture on $^{64}$Ge, not found by Cyburt et al. (2016) sensitivity study. The $^{65}$As(p,$γ$)$^{66}$Se reaction influences the abundances of nuclei $A$ = 64, 68, 72, 76, and 80 up to a factor of 1.4. The new $S_\mathrm{p}$($^{66}$Se) and the inclusion of the updated $^{22}$Mg($α$,p)$^{25}$Al reaction rate increases the production of $^{12}$C up to a factor of $4.5$ that is not observable and could be the main fuel for superburst. The waiting point status of and two-proton sequential capture on $^{64}$Ge, weak-cycle feature of GeAs at region heavier than $^{64}$Ge, and impact of other possible $S_\mathrm{p}$($^{66}$Se) are also discussed.

preprint2022arXiv

The Regulated NiCu Cycles with the new $^{57}$Cu(p,$γ$)$^{58}$Zn reaction rate and the Influence on Type-I X-Ray Bursts: GS 1826$-$24 Clocked Burster

During the X-ray bursts of GS 1826$-$24, "clocked burster", the nuclear reaction flow that surges through the rapid-proton capture process path has to pass through the NiCu cycles before reaching the ZnGa cycles that moderate the further extent of hydrogen burning in the region above germanium and selenium isotopes. The $^{57}$Cu(p,$γ$)$^{58}$Zn reaction located in the NiCu cycles plays an important role in influencing the burst light curves as found by Cyburt et al. (2016). We deduce the $^{57}$Cu(p,$γ$)$^{58}$Zn reaction rate based on the experimentally determined important nuclear structure information, isobaric-multiplet-mass equation, and large-scale shell model calculations. Based on the isobaric-multiplet-mass equation, we propose a possible order of $1^+_1$ and $2^+_3$ dominant resonance states and constrain the resonance energy of the $1^+_2$ state. The latter reduces the contribution of the $1^+_2$ dominant resonance state. The new reaction rate is up to a factor of four lower than the Forstner et al. (2001) rate recommended by JINA REACLIB v2.2 at the temperature regime sensitive to clocked bursts of GS 1826$-$24. Using the simulation from the one-dimensional implicit hydrodynamic code, KEPLER, to model the thermonuclear X-ray bursts of GS 1826$-$24 clocked burster, we find that the new $^{57}$Cu(p,$γ$)$^{58}$Zn coupled with the latest $^{56}$Ni(p,$γ$)$^{57}$Cu and $^{55}$Ni(p,$γ$)$^{56}$Cu reaction rates redistributes the reaction flow in the NiCu cycles and strongly influences the burst ash composition, whereas the $^{59}$Cu(p,$α$)$^{56}$Ni and $^{59}$Cu(p,$γ$)$^{60}$Zn reactions suppress the influence of the $^{57}$Cu(p,$γ$)$^{58}$Zn reaction and diminish the impact of nuclear reaction flow that by-passes the important $^{56}$Ni waiting point induced by the $^{55}$Ni(p,$γ$)$^{56}$Cu reaction on burst light curve.

preprint2020arXiv

Modelling Thermonuclear X-ray Bursts on Accreting Neutron Stars

In low-mass X-ray binaries, the accretion of stellar material onto a neutron star can fuel unstable thermonuclear flashes known as Type I X-ray bursts. Simulating these events using computational models can provide valuable information about the nature of the accreting system. One-dimensional (1D) astrophysics codes with large nuclear reaction networks are the current state-of-the-art for simulating X-ray bursts. These codes can track the evolution of isotopes through thousands of nuclear reaction pathways, to predict the released nuclear energy and final composition of the ashes. In this thesis, I make extensive use of KEPLER, a 1D code at the forefront of these efforts. I first present improvements to the setup and analysis of KEPLER burst models. By accounting for nuclear heating in the initial conditions, I shorten the thermal burn-in time, thereby reducing computational expense and producing more consistent burst trains. To model bursts fueled by transient accretion events, I perform the first such simulations with fully time-dependent accretion rates. Building upon previous efforts to model the "Clocked Burster", GS 1826$-$238, I precompute a grid of 3840 simulations and sample the interpolated results using Markov Chain Monte Carlo (MCMC) methods. By comparing the predictions to multi-epoch observations, I obtain posterior probability distributions for the system parameters. I then extend these MCMC methods to the pure-helium burster, 4U 1820$-$30, using a grid of 168 simulations. Finally, I discuss potential improvements for future studies, to further develop the computational modelling of accreting neutron stars.

preprint2020arXiv

Multi-epoch X-ray burst modelling: MCMC with large grids of 1D simulations

Type-I X-ray bursts are recurring thermonuclear explosions on the surface of accreting neutron stars. Matching observed bursts to computational models can help to constrain system properties, such as the neutron star mass and radius, crustal heating rates, and the accreted fuel composition, but systematic parameter studies to date have been limited. We apply Markov chain Monte Carlo methods to 1D burst models for the first time, and obtain system parameter estimations for the `Clocked Burster', GS 1826$-$238, by fitting multiple observed epochs simultaneously. We explore multiple parameters which are often held constant, including the neutron star mass, crustal heating rate, and hydrogen composition. To improve the computational efficiency, we precompute a grid of 3840 KEPLER models - the largest set of 1D burst simulations to date - and by interpolating over the model grid, we can rapidly sample burst predictions. We obtain estimates for a CNO metallicity of $Z_\mathrm{CNO} = 0.010^{+0.005}_{-0.004}$, a hydrogen fraction of $X_0 = 0.74^{+0.02}_{-0.03}$, a distance of $d \sqrt{ξ_\mathrm{b}} = 6.5^{+0.4}_{-0.6}\, \mathrm{kpc}$, and a system inclination of $i = {69^{+2}_{-3}}^{\circ}$.