Publications
2025
- A&ALine detections in photospheric radius expansion bursts from 4U 1820-303F. Barra, D. Barret, C. Pinto, and 3 more authorsarXiv e-prints, Jan 2025
Context: NICER (Neutron star Interior Composition ExploreR) is the instrument of choice for the spectral analysis of type I X-ray bursts, as it provides high throughput at X-ray CCD resolution, down to 0.3 keV. Aims: This study investigates whether the energies of absorption lines detected in photospheric radius expansion (PRE) bursts correlate with the inferred blackbody radius. Previous reports suggested such a correlation, attributed to a combination of weaker gravitational redshift and higher blueshifts in bursts with larger radii. Methods: The analysis reexamines four previously studied PRE bursts and examines eight additional bursts from 4U 1820-303, evidencing PRE. Spectral evolution is tracked on the shortest possible timescales (tenth of a second) adopting two parallel continuum descriptions to characterise the photospheric expansion and line evolution. Applying the accretion-enhanced model, maximum blackbody radii of up to ∼900 km are inferred, with peak bolometric luminosities exceeding the Eddington limit of an Helium accretor. Absorption lines are assessed for significance using Monte Carlo simulations, and spectral lines are characterised using the state-of-art plasma codes available within \scspex with a phenomenological continuum. A thorough parameter search explores Doppler shifts to avoid local minima. Results: Several significant (> 99.9%) absorption lines, including the previously reported 2.97 keV line, are detected. While no consistent correlation between line energies and blackbody radii is confirmed, bursts with larger radii exhibit up to four lines and the line strength is higher. The modelling suggests that the observed lines mostly originate from slightly redshifted (almost rest-frame) photo-/collisionally ionised gas in emission. For the burst with the largest PRE, a combination of photo-ionised plasma in both emission and absorption is preferred.
2024
- AstroNumThe Challenges of Modeling Astrophysical Reacting FlowsMichael Zingale, Khanak Bhargava, Ryan Brady, and 5 more authorsarXiv e-prints, Nov 2024
Stellar evolution is driven by the changing composition of a star from nuclear reactions. At the late stages of evolution and during explosive events, the timescale can be short and drive strong hydrodynamic flows, making simulations of astrophysical reacting flows challenging. Over the past decades, the standard approach to modeling reactions in simulation codes has been operator splitting, using implicit integrators for reactions. Here we explore some of the assumptions in this standard approach and describe some techniques for improving the efficiency and accuracy of astrophysical reacting flows.
- ApJHydrodynamical Simulations of Proton Ingestion Flashes in Type I X-Ray BurstsSimon Guichandut, Michael Zingale, and Andrew CummingThe Astrophysical Journal, Nov 2024
We perform the first multidimensional fluid simulations of thermonuclear helium ignition underneath a hydrogenrich shell. This situation is relevant to Type I X-ray bursts on neutron stars that accrete from a hydrogen-rich companion. Using the low-Mach number fluid code MAESTROeX, we investigate the growth of the convection zone due to nuclear burning, and the evolution of the chemical abundances in the atmosphere of the star. We also examine the convective boundary mixing processes that cause the evolution to differ significantly from previous one-dimensional simulations that rely on mixing-length theory. We find that the convection zone grows outward as penetrating fluid elements cool the overlying radiative layer, rather than directly from the increasing entropy of the convection zone itself. Simultaneously, these flows efficiently mix composition, carrying carbon out of and protons into the convection zone even before contact with the hydrogen shell. We discuss the implications of these effects for future modeling of these events and observations.
2023
- ApJThe Imprint of Convection on Type I X-Ray Bursts: Pauses in Photospheric Radius Expansion LightcurvesSimon Guichandut, and Andrew CummingThe Astrophysical Journal, Aug 2023
Motivated by the recent observation by NICER of a type I X-ray burst from SAX J1808.4–3658 with a distinct “pause” feature during its rise, we show that bursts which ignite in a helium layer underneath a hydrogen-rich shell naturally give rise to such pauses, as long as enough energy is produced to eject the outer layers of the envelope by super-Eddington winds. The length of the pause is determined by the extent of the convection generated after ignition, while the rate of change of luminosity following the pause is set by the hydrogen gradient left behind by convection. Using the MESA stellar evolution code, we simulate the accumulation, nuclear burning, and convective mixing prior to and throughout the ignition of the burst, followed by the hydrodynamic wind. We show that the results are sensitive to the treatment of convection adopted within the code. In particular, the efficiency of mixing at the H/He interface plays a key role in determining the shape of the lightcurve. The data from SAX J1808.4–3658 favor strong mixing scenarios. Multidimensional simulations will be needed to properly model the interaction between convection and nuclear burning during these bursts, which will then enable a new way to use X-ray burst lightcurves to study neutron star surfaces.
2021
- ApJExpanded Atmospheres and Winds in Type I X-Ray Bursts from Accreting Neutron StarsSimon Guichandut, Andrew Cumming, Maurizio Falanga, and 2 more authorsThe Astrophysical Journal, Jun 2021
We calculate steady-state models of radiation-driven super-Eddington winds and static expanded envelopes of neutron stars caused by high luminosities in type I X-ray bursts. We use flux-limited diffusion to model the transition from optically thick to optically thin, and include effects of general relativity, allowing us to study the photospheric radius close to the star as the hydrostatic atmosphere evolves into a wind. We find that the photospheric radius evolves monotonically from static envelopes (rph ≲ 50-70 km) to winds (rph ≈ 100-1000 km). Photospheric radii of less than 100 km, as observed in most photospheric radius expansion bursts, can be explained by static envelopes, but only in a narrow range of luminosity. In most bursts, we would expect the luminosity to increase further, leading to a wind with photospheric radius ≳100 km. In the contraction phase, the expanded envelope solutions show that the photosphere is still ≈1 km above the surface when the effective temperature is only 3% away from its maximum value. This is a possible systematic uncertainty when interpreting the measured Eddington fluxes from bursts at touchdown. We also discuss the applicability of steady-state models to describe the dynamics of bursts. In particular, we show that the sub- to super-Eddington transition during the burst rise is rapid enough that static models are not appropriate. Finally, we analyze the strength of spectral shifts in our models. Expected shifts at the photosphere are dominated by gravitational redshift, and are therefore predicted to be less than a few percent.