Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Neutron Source”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,099 records · Page 61

Dynamics of Laterally Propagating Flames in X-Ray Bursts. I. Burning Front Structure

We investigate the structure of laterally propagating flames through the highly stratified burning layer in an X-ray burst. Two-dimensional hydrodynamics simulations of flame propagation are performed through a rotating plane-parallel atmosphere, exploring the structure of the flame. We discuss the approximations needed to capture the length and timescales at play in an X-ray burst and describe the flame acceleration observed. Finally, our studies complement other multidimensional studies of burning in X-ray bursts.

79 ASTRONOMY AND ASTROPHYSICS↗

Could a Kilonova Kill: A Threat Assessment

Abstract Binary neutron star mergers produce high-energy emissions from several physically different sources, including a gamma-ray burst (GRB) and its afterglow, a kilonova (KN), and, at late times, a remnant many parsecs in size. Ionizing radiation from these sources can be dangerous for life on Earth-like planets when located too close. Work to date has explored the substantial danger posed by the GRB to on-axis observers; here we focus instead on the potential threats posed to nearby off-axis observers. Our analysis is based largely on observations of the GW170817/GRB 170817A multi-messenger event, as well as theoretical predictions. For baseline KN parameters, we find that the X-ray emission from the afterglow may be lethal out to ∼1 pc and the off-axis gamma-ray emission may threaten a range out to ∼4 pc, whereas the greatest threat comes years after the explosion, from the cosmic rays accelerated by the KN blast, which can be lethal out to distances up to ∼11 pc. The distances quoted here are typical, but the values have significant uncertainties and depend on the viewing angle, ejected mass, and explosion energy in ways we quantify. Assessing the overall threat to Earth-like planets, KNe have a similar kill distance to supernovae, but are far less common. However, our results rely on the scant available KN data, and multi-messenger observations will clarify the danger posed by such events.

Astronomy & Astrophysics↗

Connecting Nuclear Structure to Stellar Astrophysics: Neutron Skin in Tin Isotopes

The first observation of a neutron star merger by the LIGO-Virgo collaboration in 2017 highlights the need to improve our fundamental understanding of the equation of state of dense, neutron rich matter. The origin of heavy elements in the r-process and the structure of neutron stars are governed by the properties of neutron rich matter, for which experimental data is limited. Further analysis of this historic event and all future neutron star mergers relies on constraining the nuclear equation of state with experimental observables. We propose a novel method for systematically studying the evolution of the neutron skin in stable tin isotopes, by measuring the low-energy nuclear dipole strength over the broadest possible range of neutron-to-proton ratios in a single element. Nuclear resonance fluorescence with 100% linearly polarized photons from the High Intensity Gamma-ray Source (HIGS) facility was used to selectively measure the E1 photoabsorption strength of 112 Sn and 124 Sn at excitation energies from 3.5 MeV up to neutron separation, where the Pygmy Dipole Resonance dominates. The dipole polarizabilities of 112 Sn and 124 Sn were measured to be 9.03 ± 0.23 fm 3 and 9.11 ± 0.24 fm 3 , respectively. These uniquely systematically consistent measurements provide highly accurate experimental data for improving microscopic nuclear models as well as calculations of astrophysical nucleosynthesis and neutron star structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The HTR-Proteus Benchmark: Analysis and Use as a Verification and Validation Case

This presentation outlines the evaluation and application of the HTR-Proteus benchmark as a verification and validation (V&V) case for advanced reactor modeling tools. The work supports the U.S. Department of Energy’s HALEU Availability Program (HAP) and the joint DOE/NRC DNCSH project, which aims to reduce criticality safety uncertainties in commercial-scale HALEU fuel cycle and transportation systems. The HTR-Proteus experiments, conducted at the Paul Scherrer Institute, provide high-fidelity data for TRISO-fueled, graphite-moderated pebble bed reactors with high neutron leakage—conditions relevant to HALEU transport scenarios. This study focuses on Cores 4.2 and 4.3 of the HTR-Proteus benchmark, analyzing key sources of uncertainty including pebble packing, TRISO particle positioning, and core height. Using Project Chrono for realistic pebble geometries and Serpent, SHIFT, and MCNP for neutronics simulations, the study quantifies the impact of these uncertainties on the effective multiplication factor (keff). Results show that a sample size of 110 pebble configurations is sufficient to converge keff, with ±30 pcm uncertainty due to packing randomness. TRISO positioning and core height variations also significantly influence keff, highlighting the importance of detailed modeling in V&V efforts. The benchmark serves as a valuable test case for validating the Griffin reactor physics code and improving confidence in HALEU system simulations.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Comprehensive Analysis of Streaming and Shutdown Dose Rate Experiments at JET with ORNL Fusion Neutronics Workflows

Current experimental fusion systems and conceptual designs of fusion pilot plants (FPPs) are growing in complexity and size. Several radiation metrics are crucial to the safe operation of fusion machines, including neutron flux streaming through openings and the shutdown dose rate (SDDR). Most current designs of advanced experimental fusion systems—and the most probable candidates for FPPs—are based on the tokamak concept, which is prone to neutron streaming through the myriad openings needed for diagnostic and support systems. SDDR is caused by decay gamma rays from radionuclides that become activated by neutrons during the operation of a fusion system that use deuterium-deuterium (DD), tritium-tritium, or deuterium-tritium plasma. Because computational tools have become essential for determining these radiation metrics, they must be validated against reliable and applicable experimental data. Experiments at the Joint European Torus (JET) provide a unique source of experimental data for validating computational tools and nuclear data used to determine SDDR and neutron fluxes in streaming-dominated geometries. Here, this paper presents the comprehensive analysis of the high-performance DD JET SDDR, and streaming experiments performed using Oak Ridge National Laboratory (ORNL) fusion workflows. The computational results were compared with experimental results that consist of online SDDR measurements with ionization chambers and neutron fluence streaming measurements using thermoluminescent detectors. The ratio of calculated-to-experimental SDDR values ranges from 0.6 to 2.5, and the streaming results range from 0.5 to 8.0. Future work will include analyzing the JET 2021 DTE2 campaign alongside the integration of the Shift Monte Carlo transport code into all ORNL fusion neutronics workflows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

UFCC Calibration Measurements for Bulk U 3 O 8 Powder Samples

The uranium UF 6 coincidence counter (UFCC) was designed to measure 5A/B cylinders containing UF 6 in the passive neutron self-interrogation mode. The purpose of the present report is to extend the UFCC capability to other sample types such as oxides and liquids that might be found during verification activities. A 252 Cf source was introduced into the detector to make the active mode measurements possible, because only the UF 6 samples have adequate neutron emission for self-interrogation. The sample enrichments varied from 0.72% to 91% to establish the shape of the calibration curve and the measurement uncertainties. The measurements included the singles, doubles, and triples rates. The time correlations between the 252 Cf source and the induced fission reactions reduced the statistical error in the doubles and triples rates to be small enough for useful measurements in 5-10 minutes.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

HPC-driven computational reproducibility in numerical relativity codes: a use case study with IllinoisGRMHD

Abstract Reproducibility of results is a cornerstone of the scientific method. Scientific computing encounters two challenges when aiming for this goal. Firstly, reproducibility should not depend on details of the runtime environment, such as the compiler version or computing environment, so results are verifiable by third-parties. Secondly, different versions of software code executed in the same runtime environment should produceconsistent numerical results for physical quantities. In this manuscript, we test the feasibility of reproducing scientific results obtained using theIllinoisGRMHDcode that is part of an open-source community software for simulation in relativistic astrophysics, theEinstein Toolkit. We verify that numerical results of simulating a single isolated neutron star withIllinoisGRMHDcan be reproduced, and compare them to results reported by the code authors in 2015. We use two different supercomputers: Expanse at SDSC, and Stampede2 at TACC. By compiling the source code archived along with the paper on both Expanse and Stampede2, we find thatIllinoisGRMHDreproduces results published in its announcement paper up to errors comparable to round-off level changes in initial data parameters. We also verify that a current version ofIllinoisGRMHDreproduces these results once we account for bug fixes which have occurred since the original publication.

Astronomy & Astrophysics↗

The Impact of Classical and General Relativistic Obliquity Precessions on the Habitability of Circumstellar Neutron Stars’ Planets

Recently, it has been shown that rocky planets orbiting neutron stars can be habitable under plausible circumstances. If a distant, point-like source of visible light, such as a Sun-like main-sequence star or the gravitationally-lensed accretion disk of a supermassive black hole is present, possible temporal variations Δε{sub p}(t) of the planet’s axial tilt ε {sub p} to the ecliptic plane should be included in the overall habitability budget since the obliquity determines the insolation at a given latitude on a body’ s surface. I point out that, for rather generic initial spin–orbit initial configurations, general relativistic and classical spin variations induced by the post-Newtonian de Sitter and Lense–Thirring components of the field of the host neutron star and by its pull to the planetary oblateness J{sub 2}{sup p} may induce huge and very fast variations of ε {sub p} that would likely have an impact on the habitability of such worlds. In particular, for a planet’s distance of, say, 0.005 au from a 1.4 M {sub ⊙} neutron star corresponding to an orbital period P {sub b} = 0.109 day, obliquity shifts Δε {sub p} as large as ε{sub p}{sup max}−ε{sub p}{sup min}≃50{sup ∘}--100{sup ∘} over characteristic timescales as short as 10 days (J{sub 2}{sup p}) to 3 Myr (Lense–Thirring) may occur for arbitrary orientations of the orbital and spin angular momenta L, S {sub ns}, S {sub p} of the planet-neutron star system. In view of this feature of their spins, I dub such hypothetical planets as “nethotrons.”.

79 ASTRONOMY AND ASTROPHYSICS↗

FENIX: An Open-Source Multiphysics Integrated Framework Enabling Collaborative Development of Plasma Facing Component Modeling Capabilities

Advanced modeling and simulation tools have a crucial role to play in accelerating fusion energy deployment as a sustainable power source. Multiphysics, high-fidelity computational tools can help understand, model, and quantify the complex interactions between materials performance, plasma and neutron exposure, and engineering processes. As a result, they accelerate the design, safety analysis, and performance evaluation of fusion systems. This webinar introduces the Fusion ENergy Integrated multiphys-X (FENIX) framework, an open-source multiphysics tool for plasma facing component modeling. FENIX leverages the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, which has been developed by the United States Department of Energy Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. FENIX couples various MOOSE capabilities such as heat transfer, thermomechanics, thermal hydraulics, electromagnetics, and plasma kinetics with the MOOSE-based applications Cardinal (neutronics) and TMAP8 (tritium transport). During the webinar, we will present FENIX and discuss how its modularity, openness, software quality assurance processes, and licensing approach supports effective collaborations, including public-private partnerships.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Complete β-decay patterns of 142 Cs, 142 Ba, and 142 La determined using total absorption spectroscopy

Background: The β decays of fission products produced in nuclear fuel are important for nuclear energy applications and fundamental science of reactor antineutrinos. In particular, nuclear reactor safety is related to the decay modes of radioactive neutron-rich nuclei, primarily via the emission of γ rays, neutrons, and electrons. Additionally nuclear reactors are the most powerful man-made source of antineutrinos emitted during the β decay of fission products. These antineutrinos are used to inspect fundamental properties of leptons as well as informing reactor operation. However, the majority of data on complex decays of fission products collected in the evaluated nuclear data repositories like Evaluated Nuclear Structure Data File (ENSDF) and Evaluated Nuclear Data Files (ENDF) are based on low-efficiency and often incomplete measurements resulting in questionable reference reactor antineutrino flux predictions, see the analysis by [Nichols, J. Nucl. Sci. Technol. 52, 17 (2015)]. Various assessments like the one done under the auspices of the [Yoshida et al., Assessment of Fission Product Decay Data for Decay Heat Calculations: A report by the Working Party on International Evaluation Co-operation of the Nuclear Energy Agency Nuclear Science Committee (Nuclear Energy Agency, Organization for Economic Co-operation and Development, Paris, France, 2007), Vol. 25], as well as by [Sonzogni, Johnson, and McCutchan, Phys. Rev. C 91, 011301(R) (2015)] and [Dwyer and Langford, Phys. Rev. Lett. 114, 012502 (2015)], list the A = 142 isobars with high cumulative fission yield among the important nuclei where data for reactor decay heat and/or antineutrino production should be verified and/or improved. Purpose: Here, our goal is to improve the quality of β -decay measurements and evaluate the impact of modified decay schemes on reactor decay heat and antineutrino energy spectra, for fission products along the A = 142 isobaric chain. This work is an in depth follow-up on [Rasco et al., Phys. Rev. Lett. 117, 092501 (2016)]. which presented briefly the impact of the corrected decay scheme of 142 Cs . Here, we extend the data to full isobaric decay chain including the daughter nuclei, 142 Ba and 142 La, and present more details on the 142 Cs results. Method: The decays of neutron-rich isobars of mass A = 142 produced by means of proton-induced fission of 238 U were measured using the Modular Total Absorption Spectrometer (MTAS) array on-line at the mass separator and Tandem accelerator at Oak Ridge National Laboratory. Results: The β -decay schemes for 142 Cs and 142 La were modified with respect to the nuclear data repositories. A small β-delayed neutron branching ratio for 142 Cs emitter was remeasured as $0.10^{+5}_{–3}% %. Improved precision on the measured half-lives is reported. Small corrections to the low-energy decay of 142 Ba are made. The β-decay patterns for 142 La and 142 Cs are presented. The decay heat release and cross section for the detection of reactor antineutrinos are deduced and compared to earlier results. Conclusions: The β-feeding pattern for 142 Cs having decay energy value $Q_β$ of over 7 MeV was substantially modified with respect to the current ENSDF entry. Smaller changes were encountered for 142 La, but since this A = 142 isobar also has a large cumulative yield in fission, the changes influence both decay heat and the antineutrino spectra. The previously known β intensities for 142 Ba decay ($Q_β$ value of 2.2 MeV) were verified and slightly modified. Overall, increased decay heat values and lower flux of antineutrinos interacting with matter are presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Gamma Source Verification for GAMSRC and GAMSOR

Nuclear reactors that rely upon the fission reaction have two modes of thermal energy deposition in the reactor system: neutron absorption and gamma absorption. The gamma rays are typically generated by neutron capture reactions or during the fission process which means the primary driver of energy production is of course the neutron interactions. The GAMSOR program was first built in the mid 1980s to properly account for the gamma heating in an operating reactor core on core internals. The GAMSOR code is sequence of DIF3D calculations to compute the neutron and gamma flux and combine them to define both the neutron and gamma heating throughout the modeled domain. The goal of this manuscript is to present the software verification of GAMSOR. The first step of the GAMSOR sequence of calculations involves of a modified version of DIF3D (called DIF3D-GAMSOR) which generates a gamma source distribution for the follow-on DIF3D gamma transport calculation (step 2). This modified version of DIF3D increases the burden of maintenance and verification work on GAMSOR as one must reverify the DIF3D capabilities which is undesirable. Because the calculation of the gamma source is the only unique aspect of GAMSOR beyond the regular DIF3D capabilities, that part was put in a standalone code called GAMSRC such that one can use the verified DIF3D code in step 1 followed by GAMSRC to carry out the same GAMSOR calculation step. As a consequence, this manuscript is focused on verification of the gamma source files generated by GAMSRC. The verification of the modified version of DIF3D (DIF3D-GAMSOR) will be done less rigorously in that it will be verified that it produces the same output that GAMSRC does and thus GAMSRC is equivalent to GAMSOR on the problems studied here. Hand calculations and independent numerical calculations of the gamma source generation are used for the verification work. This work follows the same methodology of GAMSRC. As expected, the results agree well with those calculated by GAMSRC as will be shown.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Gamma Source Verification for GAMSRC and GAMSOR

Nuclear reactors that rely upon the fission reaction have two modes of thermal energy deposition in the reactor system: neutron absorption and gamma absorption. The gamma rays are typically generated by neutron capture reactions or during the fission process which means the primary driver of energy production is of course the neutron interactions. The GAMSOR program was first built in the mid 1980s to properly account for the gamma heating in an operating reactor core on core internals. The GAMSOR code is sequence of DIF3D calculations to compute the neutron and gamma flux and combine them to define both the neutron and gamma heating throughout the modeled domain. The goal of this manuscript is to present the software verification of GAMSOR. The first step of the GAMSOR sequence of calculations involves of a modified version of DIF3D (called DIF3D-GAMSOR) which generates a gamma source distribution for the follow-on DIF3D gamma transport calculation (step 2). This modified version of DIF3D increases the burden of maintenance and verification work on GAMSOR as one must reverify the DIF3D capabilities which is undesirable. Because the calculation of the gamma source is the only unique aspect of GAMSOR beyond the regular DIF3D capabilities, that part was put in a standalone code called GAMSRC such that one can use the verified DIF3D code in step 1 followed by GAMSRC to carry out the same GAMSOR calculation step. As a consequence, this manuscript is focused on verification of the gamma source files generated by GAMSRC. The verification of the modified version of DIF3D (DIF3D-GAMSOR) will be done less rigorously in that it will be verified that it produces the same output that GAMSRC does and thus GAMSRC is equivalent to GAMSOR on the problems studied here. Hand calculations and independent numerical calculations of the gamma source generation are used for the verification work. This work follows the same methodology of GAMSRC. As expected, the results agree well with those calculated by GAMSRC as will be shown.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Towards the Development of Molten Salt Reactor Multiphysics Modeling Capabilities: Depletion, Thermochemistry, Reactor Dynamics, and Species Tracking

This report documents the simulation and modeling work performed under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, utilizing the tools and codes developed under the program. The focus in this Fiscal Year (FY) is on the depletion, thermochemistry, reactor dynamics, and species transport in MSRs. The fission products travelling in the entire primary loop is unique in MSRs. There are many physical/chemical mechanisms that can significantly affect the fission products concentration in the fuel salt, e.g., diffusion into gas bubbles, diffusion in the core moderators, deposition on the surface of structures/pipes, etc. The fission product concentration in the fuel salt can potentially affect the neutronics, system dynamics, chemical potentials, etc., so accurately predicting the fission products distribution in the primary loop is important. Under the main objective, different activities are pursued to assess the NEAMS tools applicability for a variety of areas on species transport, including the depletion analyses using Griffin (fission product source), the redox potential assessment at different depletion points, the assessment of xenon transport effect on the neutronics and reactor dynamics, a survey of correlations for calculating the mass transfer coefficient of different fission products, etc.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

SAM Code Development for Source Term Modeling in Fluoride-salt-cooled High-temperature Reactors

The SAM code is under development and supported by DOE-NE’s Nuclear Energy Advanced Modeling and Simulation program as a modern system-level modeling and simulation tool for advanced non-light water reactor safety analyses. These advanced reactor concepts incorporate novel and improved approaches to achieve safety and economic feasibility. Many concepts leverage higher operating temperatures for improved efficiency with a variety of coolants and structural materials to support those needs. Such design choices may also pose unique radiological source term risks. There are continuous efforts to incorporate new physics and phenomena relevant to advanced reactor concepts, including capabilities to enable system-level source term modeling. The Pebble-Bed Fluoride-salt-cooled High-temperature Reactor (PB-FHR) is a promising candidate among advanced nuclear reactor concepts with its improved passive safety characteristics and high thermal efficiency. In addition to past efforts to support the development and utilization of SAM for PB-FHR safety analysis, the species transport modeling capabilities in SAM have been extended to simulate source term phenomena in concepts like the PB-FHR. The PB-FHR concept utilizes pebble-form TRISO fuel and FLiBe salt coolant to provide robust barriers to the release of almost all fission products and radiological source terms. However, tritium poses a unique radiological risk due to its significant production from neutron interactions with 6 Li and 9 Be in the FLiBe as well as its high mobility at elevated temperatures where it can permeate through structural metals. Secondly, the use of graphite at elevated temperatures poses the risk of oxidation damage to any structural components or fuel pebbles that are exposed to accidental air ingress. This report summarizes progress made in modeling these phenomena in SAM, which leverages the system-level multi-physics thermal hydraulic simulation to support an effective engineering-scale source term modeling capability. A tritium transport model is developed to simulate the various transport pathways in FHR and MSR concepts and is presented through verification and validation examples as well as demonstrations simulating experimental test loops and reference plant FHR models. A graphite oxidation model is also introduced to calculate local oxidation rates and is presented with initial validation results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The Legendre Polynomial Axial Expansion Method

This work presents a new formulation of the axial expansion transport method explicitly using Legendre polynomials for arbitrarily high-order expansions. This new formulation also features an alternative method of axial leakage calculation to allow for nonextruded flat source region meshes. This alternative axial leakage is introduced alongside a balance equation requirement to ensure that neutron balance is preserved in the coarse mesh for a given axial leakage formulation, which allows for effective coarse mesh finite difference acceleration. A matrix exponential table method is derived to allow for fast computations of arbitrarily high-order matrix exponentials for this work and precludes the need for further research into matrix exponential calculations for this method. Numerical results are presented that demonstrate the stability of the axial expansion method in systems with voidlike regions, showcase the speedup from matrix exponential tables, and investigate the axial convergence of the method in terms of both expansion order and mesh size.

Herring, Nicholas↗

Enhanced Neutron and γ-Ray Detection via 6 Li Substitution in Undoped and Tl-Doped Zero-Dimensional Perovskite Cs 3 Cu 2 I 5 Scintillators

Radiation detectors are crucial in a wide variety of research and commercial applications, such as oil and gas exploration, medical imaging, nuclear nonproliferation, and homeland security. Neutron and gamma-ray detectors are fundamental components in portal monitors at ports and border crossings, bolstering national security against radiological threats. This study presents a dual-mode scintillator, undoped and Tl-doped 6 Li-Cs 3 Cu 2 I 5 , and demonstrates its potential as a promising material for simultaneous thermal neutron and gamma-ray detection. We explore the Bridgman growth of both undoped and thallium doped Li → Cu and Li → Cs substitutional systems with various Li doping levels and assess their impact on scintillation properties. Under 662 keV gamma-ray excitation, the undoped crystals had light yields up to 35,900 ph/MeV, with energy resolutions down to 4.5%. The Tl-doped crystals performed better than the undoped crystals with light yields peaking at 65,900 ph/MeV and energy resolutions as low as 3.5%. When exposed to a moderated 252 Cf excitation source, our crystals had light yields between 102,900 and 167,200 photons per thermal neutron capture, with a full energy thermal neutron peak reaching 3 MeV in gamma equivalent energy. Pulse shape discrimination studies reveal well-separated gamma and neutron events, resulting in Figure-Of-Merit (FOM) as high as 3.7. Furthermore, these findings highlight the potential of Li-doped Cs 3 Cu 2 I 5 as a viable candidate for next-generation dual-mode scintillators.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spinteract: a program to refine magnetic interactions to diffuse scattering data

Magnetic diffuse scattering—the broad magnetic scattering features observed in neutron-diffraction data above a material's magnetic ordering temperature—provides a rich source of information about the material's magnetic Hamiltonian. However, this information has often remained under-utilised due to a lack of available computer software that can fit values of magnetic interaction parameters to such data. Here, an open-source computer program, Spinteract, is presented, which enables straightforward refinement of magnetic interaction parameters to powder and single-crystal magnetic diffuse scattering data. Here, the theory and implementation of this approach are summarised. Examples are presented of refinements to published experimental diffuse-scattering data sets for the canonical antiferromagnet MnO and the highly-frustrated classical spin liquid Gd 3 Ga 5 O 12 . Guidelines for data collection and refinement are outlined, and possible developments of the approach are discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Pulse profile modelling of thermonuclear burst oscillations – II. Handling variability

ABSTRACT Pulse profile modelling is a relativistic ray-tracing technique that can be used to infer masses, radii, and geometric parameters of neutron stars. In a previous study, we looked at the performance of this technique when applied to thermonuclear burst oscillations from accreting neutron stars. That study showed that ignoring the variability associated with burst oscillation sources resulted in significant biases in the inferred mass and radius, particularly for the high count rates that are nominally required to obtain meaningful constraints. In this follow-on study, we show that the bias can be mitigated by slicing the bursts into shorter segments where variability can be neglected, and jointly fitting the segments. Using this approach, the systematic uncertainties on the mass and radius are brought within the range of the statistical uncertainty. With about 106 source counts, this yields uncertainties of approximately 10 per cent for both the mass and radius. However, this modelling strategy requires substantial computational resources. We also confirm that the posterior distributions of the mass and radius obtained from multiple bursts of the same source can be merged to produce outcomes comparable to that of a single burst with an equivalent total number of counts.

Astronomy & Astrophysics↗