-
Pressure anisotropy and viscous heating in weakly collisional plasma turbulence
Authors:
Jonathan Squire,
Matthew W Kunz,
Lev Arzamasskiy,
Zade Johnston,
Eliot Quataert,
Alexander A Schekochihin
Abstract:
Pressure anisotropy can strongly influence the dynamics of weakly collisional, high-beta plasmas, but its effects are missed by standard magnetohydrodynamics (MHD). Small changes to the magnetic-field strength generate large pressure-anisotropy forces, heating the plasma, driving instabilities, and rearranging flows, even on scales far above the particles' gyroscales where kinetic effects are trad…
▽ More
Pressure anisotropy can strongly influence the dynamics of weakly collisional, high-beta plasmas, but its effects are missed by standard magnetohydrodynamics (MHD). Small changes to the magnetic-field strength generate large pressure-anisotropy forces, heating the plasma, driving instabilities, and rearranging flows, even on scales far above the particles' gyroscales where kinetic effects are traditionally considered important. Here, we study the influence of pressure anisotropy on turbulent plasmas threaded by a mean magnetic field (Alfvénic turbulence). Extending previous results that were concerned with Braginskii MHD, we consider a wide range of regimes and parameters using a simplified fluid model based on drift kinetics with heat fluxes calculated using a Landau-fluid closure. We show that viscous (pressure-anisotropy) heating dissipates between a quarter and half of the turbulent cascade power injected at large scales; this does not depend strongly on either plasma beta or the ion-to-electron temperature ratio. This will in turn influence the plasma's thermodynamics by regulating energy partition between different dissipation channels (e.g., electron and ion heat). Due to the pressure anisotropy's rapid dynamical feedback onto the flows that create it -- an effect we term `magneto-immutability' -- the viscous heating is confined to a narrow range of scales near the forcing scale, supporting a nearly conservative, MHD-like inertial-range cascade, via which the rest of the energy is transferred to small scales. Despite the simplified model, our results -- including the viscous heating rate, distributions, and turbulent spectra -- compare favourably to recent hybrid-kinetic simulations. This is promising for the more general use of extended-fluid (or even MHD) approaches to model weakly collisional plasmas such as the intracluster medium, hot accretion flows, and the solar wind.
△ Less
Submitted 6 July, 2023; v1 submitted 20 February, 2023;
originally announced March 2023.
-
Robust inference of neutron-star parameters from thermonuclear burst observations
Authors:
D. K. Galloway,
Z. Johnston,
A. J. Goodwin,
C. -C. He
Abstract:
Thermonuclear (type-I) bursts arise from unstable ignition of accumulated fuel on the surface of neutron stars in low-mass X-ray binaries. Measurements of burst properties in principle enable observers to infer the properties of the host neutron star and mass donors, but a number of confounding astrophysical effects contribute to systematic uncertainties. Here we describe some commonly-used approa…
▽ More
Thermonuclear (type-I) bursts arise from unstable ignition of accumulated fuel on the surface of neutron stars in low-mass X-ray binaries. Measurements of burst properties in principle enable observers to infer the properties of the host neutron star and mass donors, but a number of confounding astrophysical effects contribute to systematic uncertainties. Here we describe some commonly-used approaches for determining system parameters, including composition of the burst fuel, and introduce a new suite of software tools, concord, intended to fully account for astrophysical uncertainties. Comparison of observed burst properties with the predictions of numerical models is a complementary method of constraining host properties, and the tools presented here are intended to make comprehensive model-observation comparisons straightforward. When combined with the extensive samples of burst observations accumulated by X-ray observatories, these software tools will provide a step-change in the amount of information that can be inferred about typical burst sources.
△ Less
Submitted 7 October, 2022;
originally announced October 2022.
-
On the properties of Alfvénic switchbacks in the expanding solar wind: the influence of the Parker spiral
Authors:
Jonathan Squire,
Zade Johnston,
Alfred Mallet,
Romain Meyrand
Abstract:
Switchbacks -- rapid, large deflections of the solar wind's magnetic field -- have generated significant interest as possible signatures of the key mechanisms that heat the corona and accelerate the solar wind. In this context, an important task for theories of switchback formation and evolution is to understand their observable distinguishing features, allowing them to be assessed in detail using…
▽ More
Switchbacks -- rapid, large deflections of the solar wind's magnetic field -- have generated significant interest as possible signatures of the key mechanisms that heat the corona and accelerate the solar wind. In this context, an important task for theories of switchback formation and evolution is to understand their observable distinguishing features, allowing them to be assessed in detail using spacecraft data. Here, we work towards this goal by studying the influence of the Parker spiral on the evolution of Alfvénic switchbacks in an expanding plasma. Using simple analytic arguments based on the physics of one-dimensional spherically polarized (constant-field-magnitude) Alfvén waves, we find that, by controlling the wave's obliquity, a Parker spiral strongly impacts switchback properties. Surprisingly, the Parker spiral can significantly enhance switchback formation, despite normalized wave amplitudes growing more slowly in its presence. In addition, switchbacks become strongly asymmetric: large switchbacks preferentially involve magnetic-field rotation in the plane of the Parker spiral (tangential deflections) rather than perpendicular (normal) rotations, and such deflections are strongly "tangentially skewed," meaning switchbacks always involve field rotations in the same direction (towards the positive-radial direction for an outwards mean field). In a companion paper, we show that these properties also occur in turbulent 3-D fields with switchbacks, with various caveats. These results demonstrate that substantial care is needed in assuming that specific features of switchbacks can be used to infer properties of the low corona; asymmetries and nontrivial correlations can develop as switchbacks propagate due to the interplay between expansion and spherically polarized, divergence-free magnetic fields.
△ Less
Submitted 23 May, 2022; v1 submitted 19 May, 2022;
originally announced May 2022.
-
On the properties of Alfvénic switchbacks in the expanding solar wind: three-dimensional numerical simulations
Authors:
Zade Johnston,
Jonathan Squire,
Alfred Mallet,
Romain Meyrand
Abstract:
Switchbacks -- abrupt reversals of the magnetic field within the solar wind -- have been ubiquitously observed by Parker Solar Probe (PSP). Their origin, whether from processes near the solar surface or within the solar wind itself, remains under debate, and likely has key implications for solar wind heating and acceleration. Here, using three-dimensional expanding box simulations, we examine the…
▽ More
Switchbacks -- abrupt reversals of the magnetic field within the solar wind -- have been ubiquitously observed by Parker Solar Probe (PSP). Their origin, whether from processes near the solar surface or within the solar wind itself, remains under debate, and likely has key implications for solar wind heating and acceleration. Here, using three-dimensional expanding box simulations, we examine the properties of switchbacks arising from the evolution of outwards-propagating Alfvén waves in the expanding solar wind in detail. Our goal is to provide testable predictions that can be used to differentiate between properties arising from solar surface processes and those from the in-situ evolution of Alfvén waves in switchback observations by PSP. We show how the inclusion of the Parker spiral causes magnetic field deflections within switchbacks to become asymmetric, preferentially deflecting in the plane of the Parker spiral and rotating in one direction towards the radial component of the mean field. The direction of the peak of the magnetic field distribution is also shown to be different from the mean field direction due to its highly skewed nature. Compressible properties of switchbacks are also explored, with magnetic-field-strength and density fluctuations being either correlated or anticorrelated depending on the value of $β$, agreeing with predictions from theory. We also measure dropouts in magnetic-field strength and density spikes at the boundaries of these synthetic switchbacks, both of which have been observed by PSP. The agreement of these properties with observations provide further support for the Alfvén wave model of switchbacks.
△ Less
Submitted 2 December, 2022; v1 submitted 19 May, 2022;
originally announced May 2022.
-
Horizons: Nuclear Astrophysics in the 2020s and Beyond
Authors:
H. Schatz,
A. D. Becerril Reyes,
A. Best,
E. F. Brown,
K. Chatziioannou,
K. A. Chipps,
C. M. Deibel,
R. Ezzeddine,
D. K. Galloway,
C. J. Hansen,
F. Herwig,
A. P. Ji,
M. Lugaro,
Z. Meisel,
D. Norman,
J. S. Read,
L. F. Roberts,
A. Spyrou,
I. Tews,
F. X. Timmes,
C. Travaglio,
N. Vassh,
C. Abia,
P. Adsley,
S. Agarwal
, et al. (140 additional authors not shown)
Abstract:
Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilit…
▽ More
Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilities across an ever growing number of disciplines and subfields that need to be integrated. We take a holistic view of the field discussing the unique challenges and opportunities in nuclear astrophysics in regards to science, diversity, education, and the interdisciplinarity and breadth of the field. Clearly nuclear astrophysics is a dynamic field with a bright future that is entering a new era of discovery opportunities.
△ Less
Submitted 16 May, 2022;
originally announced May 2022.
-
Comparison of Electron Capture Rates in the N=50 Region using 1D Simulations of Core-collapse Supernovae
Authors:
Zac Johnston,
Sheldon Wasik,
Rachel Titus,
MacKenzie L. Warren,
Evan P. O'Connor,
Remco Zegers,
Sean M. Couch
Abstract:
Recent studies have highlighted the sensitivity of core-collapse supernovae (CCSNe) models to electron-capture (EC) rates on neutron-rich nuclei near the N=50 closed-shell region. In this work, we perform a large suite of one-dimensional CCSN simulations for 200 stellar progenitors using recently updated EC rates in this region. For comparison, we repeat the simulations using two previous implemen…
▽ More
Recent studies have highlighted the sensitivity of core-collapse supernovae (CCSNe) models to electron-capture (EC) rates on neutron-rich nuclei near the N=50 closed-shell region. In this work, we perform a large suite of one-dimensional CCSN simulations for 200 stellar progenitors using recently updated EC rates in this region. For comparison, we repeat the simulations using two previous implementations of EC rates: a microphysical library with parameterized N=50 rates (LMP), and an older independent-particle approximation (IPA). We follow the simulations through shock revival up to several seconds post-bounce, and show that the EC rates produce a consistent imprint on CCSN properties, often surpassing the role of the progenitor itself. Notable impacts include the timescale of core collapse, the electron fraction and mass of the inner core at bounce, the accretion rate through the shock, the success or failure of revival, and the properties of the central compact remnant. We also compare the observable neutrino signal of the neutronization burst in a DUNE-like detector, and find consistent impacts on the counts and mean energies. Overall, the updated rates result in properties that are intermediate between LMP and IPA, and yet slightly more favorable to explosion than both.
△ Less
Submitted 11 November, 2022; v1 submitted 18 February, 2022;
originally announced February 2022.
-
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
Authors:
Yi Hua Lam,
Zi Xin Liu,
Alexander Heger,
Ning Lu,
Adam Michael Jacobs,
Zac Johnston
Abstract:
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 mu…
▽ More
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.
△ Less
Submitted 11 January, 2022; v1 submitted 26 October, 2021;
originally announced October 2021.
-
Advancement of Photospheric Radius Expansion and Clocked Type-I X-Ray Burst Models with the New $^{22}$Mg$(α,p)^{25}$Al Reaction Rate Determined at Gamow Energy
Authors:
J. Hu,
H. Yamaguchi,
Y. H. Lam,
A. Heger,
D. Kahl,
A. M. Jacobs,
Z. Johnston,
S. W. Xu,
N. T. Zhang,
S. B. Ma,
L. H. Ru,
E. Q. Liu,
T. Liu,
S. Hayakawa,
L. Yang,
H. Shimizu,
C. B. Hamill,
A. St J. Murphy,
J. Su,
X. Fang,
K. Y. Chae,
M. S. Kwag,
S. M. Cha,
N. N. Duy,
N. K. Uyen
, et al. (12 additional authors not shown)
Abstract:
We report the first (in)elastic scattering measurement of $^{25}\mathrm{Al}+p$ with the capability to select and measure in a broad energy range the proton resonances in $^{26}$Si contributing to the $^{22}$Mg$(α,p)$ reaction at type I x-ray burst energies. We measured spin-parities of four resonances above the $α$ threshold of $^{26}$Si that are found to strongly impact the $^{22}$Mg$(α,p)$ rate.…
▽ More
We report the first (in)elastic scattering measurement of $^{25}\mathrm{Al}+p$ with the capability to select and measure in a broad energy range the proton resonances in $^{26}$Si contributing to the $^{22}$Mg$(α,p)$ reaction at type I x-ray burst energies. We measured spin-parities of four resonances above the $α$ threshold of $^{26}$Si that are found to strongly impact the $^{22}$Mg$(α,p)$ rate. The new rate advances a state-of-the-art model to remarkably reproduce light curves of the GS 1826$-$24 clocked burster with mean deviation $<9$ % and permits us to discover a strong correlation between the He abundance in the accreting envelope of photospheric radius expansion burster and the dominance of $^{22}$Mg$(α,p)$ branch.
△ Less
Submitted 20 October, 2021; v1 submitted 10 August, 2021;
originally announced August 2021.
-
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
Authors:
Yi Hua Lam,
Ning Lu,
Alexander Heger,
Adam Michael Jacobs,
Nadezda A. Smirnova,
Teresa Kurtukian Nieto,
Zac Johnston,
Shigeru Kubono
Abstract:
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 importa…
▽ More
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.
△ Less
Submitted 15 January, 2022; v1 submitted 24 July, 2021;
originally announced July 2021.
-
The efficiency of nuclear burning during thermonuclear (Type I) bursts as a function of accretion rate
Authors:
Y. Cavecchi,
D. K. Galloway,
A. J. Goodwin,
Z. Johnston,
A. Heger
Abstract:
We measured the thermonuclear burning efficiency as a function of accretion rate for the Type I X-ray bursts of five low-mass X-ray binary systems. We chose sources with measured neutron star spins and a substantial population of bursts from a large observational sample. The general trend for the burst rate is qualitatively the same for all sources; the burst rate first increases with the accretio…
▽ More
We measured the thermonuclear burning efficiency as a function of accretion rate for the Type I X-ray bursts of five low-mass X-ray binary systems. We chose sources with measured neutron star spins and a substantial population of bursts from a large observational sample. The general trend for the burst rate is qualitatively the same for all sources; the burst rate first increases with the accretion rate up to a maximum, above which the burst rate declines, despite the increasing accretion rate. At higher accretion rates, when the burst rate decreases, the α-value (the ratio of accretion energy and burst energy) increases by up to a factor of 10 above that in the rising burst rate regime. These observations are contrary to the predictions of 1D numerical models, but can be explained as the consequence of a zone of stable burning on the neutron star surface, which expands with increasing accretion rate. The stable burning also "pollutes" the unstable burning layer with ashes, contributing to the change in burst properties measured in the falling burst rate regime. We find that the mass accretion rate at which the burst rate begins to decrease is anti-correlated with the spin of the neutron star. We conclude that the neutron star spin is a key factor, moderating the nuclear burning stability, via the local accretion rate and fuel composition over the star.
△ Less
Submitted 15 July, 2020;
originally announced July 2020.
-
Modelling Thermonuclear X-ray Bursts on Accreting Neutron Stars
Authors:
Zac Johnston
Abstract:
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 simulati…
▽ More
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.
△ Less
Submitted 31 March, 2020;
originally announced April 2020.
-
Multi-epoch X-ray burst modelling: MCMC with large grids of 1D simulations
Authors:
Zac Johnston,
Alexander Heger,
Duncan K. Galloway
Abstract:
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…
▽ More
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}$.
△ Less
Submitted 16 April, 2020; v1 submitted 17 September, 2019;
originally announced September 2019.
-
A Bayesian Approach to Matching Thermonuclear X-ray Burst Observations with Models
Authors:
A. J. Goodwin,
D. K. Galloway,
A. Heger,
A. Cumming,
Z. Johnston
Abstract:
We present a new method of matching observations of Type I (thermonuclear) X-ray bursts with models, comparing the predictions of a semi-analytic ignition model with X-ray observations of the accretion-powered millisecond pulsar SAX J1808.4--3658 in outburst. We used a Bayesian analysis approach to marginalise over the parameters of interest and determine parameters such as fuel composition, dista…
▽ More
We present a new method of matching observations of Type I (thermonuclear) X-ray bursts with models, comparing the predictions of a semi-analytic ignition model with X-ray observations of the accretion-powered millisecond pulsar SAX J1808.4--3658 in outburst. We used a Bayesian analysis approach to marginalise over the parameters of interest and determine parameters such as fuel composition, distance/anisotropy factors, neutron star mass and neutron star radius. Our study includes a treatment of the system inclination effects, inferring that the rotation axis of the system is inclined $\left(69^{+4}_{-2}\right)^\circ$ from the observers line of sight, assuming a flat disc model. This method can be applied to any accreting source that exhibits Type I X-ray bursts. We find a hydrogen mass fraction of $0.57^{+0.13}_{-0.14}$ and CNO metallicity of $0.013^{+0.006}_{-0.004}$ for the accreted fuel is required by the model to match the observed burst energies, for a distance to the source of $3.3^{+0.3}_{-0.2}\,\mathrm{kpc}$. We infer a neutron star mass of $1.5^{+0.6}_{-0.3}\,\mathrm{M}_{\odot}$ and radius of $11.8^{+1.3}_{-0.9}\,\mathrm{km}$ for a surface gravity of $1.9^{+0.7}_{-0.4}\times10^{14}\,\mathrm{cm}\,\mathrm{s}^{-2}$ for SAX J1808.4--3658.
△ Less
Submitted 17 September, 2019; v1 submitted 1 July, 2019;
originally announced July 2019.
-
High-energy transients: thermonuclear (type-I) X-ray bursts
Authors:
Duncan K. Galloway,
Zac Johnston,
Adelle Goodwin,
Alexander Heger
Abstract:
Many distinct classes of high-energy variability have been observed in astrophysical sources, on a range of timescales. The widest range (spanning microseconds-decades) is found in accreting, stellar-mass compact objects, including neutron stars and black holes. Neutron stars are of particular observational interest, as they exhibit surface effects giving rise to phenomena (thermonuclear bursts an…
▽ More
Many distinct classes of high-energy variability have been observed in astrophysical sources, on a range of timescales. The widest range (spanning microseconds-decades) is found in accreting, stellar-mass compact objects, including neutron stars and black holes. Neutron stars are of particular observational interest, as they exhibit surface effects giving rise to phenomena (thermonuclear bursts and pulsations) not seen in black holes. Here we briefly review the present understanding of thermonuclear (type-I) X-ray bursts. These events are powered by an extensive chain of nuclear reactions, which are in many cases unique to these environments. Thermonuclear bursts have been exploited over the last few years as an avenue to measure the neutron star mass and radius, although the contribution of systematic errors to these measurements remains contentious. We describe recent efforts to better match burst models to observations, with a view to resolving some of the astrophysical uncertainties related to these events. These efforts have good prospects for providing complementary information to nuclear experiments.
△ Less
Submitted 5 September, 2019; v1 submitted 1 March, 2018;
originally announced March 2018.
-
Simulating X-ray bursts during a transient accretion event
Authors:
Zac Johnston,
Alexander Heger,
Duncan K. Galloway
Abstract:
Modelling of thermonuclear X-ray bursts on accreting neutron stars has to date focused on stable accretion rates. However, bursts are also observed during episodes of transient accretion. During such events, the accretion rate can evolve significantly between bursts, and this regime provides a unique test for burst models. The accretion-powered millisecond pulsar SAX J1808.4-3658 exhibits accretio…
▽ More
Modelling of thermonuclear X-ray bursts on accreting neutron stars has to date focused on stable accretion rates. However, bursts are also observed during episodes of transient accretion. During such events, the accretion rate can evolve significantly between bursts, and this regime provides a unique test for burst models. The accretion-powered millisecond pulsar SAX J1808.4-3658 exhibits accretion outbursts every 2-3 years. During the well-sampled month-long outburst of 2002 October, four helium-rich X-ray bursts were observed. Using this event as a test case, we present the first multi-zone simulations of X-ray bursts under a time-dependent accretion rate. We investigate the effect of using a time-dependent accretion rate in comparison to constant, averaged rates. Initial results suggest that using a constant, average accretion rate between bursts may underestimate the recurrence time when the accretion rate is decreasing, and overestimate it when the accretion rate is increasing. Our model, with an accreted hydrogen fraction of $X=0.44$ and a CNO metallicity of $Z_\mathrm{CNO}=0.02$, reproduces the observed burst arrival times and fluences with root mean square (RMS) errors of $2.8\,\mathrm{h}$, and $0.11\times 10^{-6}\,\mathrm{erg\, cm^{-2}}$, respectively. Our results support previous modelling that predicted two unobserved bursts, and indicate that additional bursts were also missed by observations.
△ Less
Submitted 20 March, 2018; v1 submitted 10 November, 2017;
originally announced November 2017.