Neutron Stars

X-ray Bursting Neutron Stars

Get the Castro code for Compressible ASTROphysics

Our Castro code and our simulation inputs files are open-source and freely available at the Castro GitHub.

Solving nuclear reactions on GPUs enabled accelerated nuclear astrophysics simulations using the Castro compressible hydrodynamics code. For more details about how we accelerated our nuclear reaction solver, see Nuclear Reactions.

My GPU work enabled our nuclear astrophysics collaboration to simulate detailed two and three-dimensional nuclear flames on neutron star surfaces!

Neutron stars in binary star systems can under certain astrophysical conditions accrete material from their companion star. In these cases, the neutron star’s atmosphere can undergo periodic thermonuclear runaway burning visible from Earth as a bright burst of X-rays.

By carefully modeling the X-ray burst physics, we hope to correlate observed X-ray bursts with properties of the unseen, super-dense neutron star core hiding below its atmosphere.

Nuclear energy generation rate within a Helium fusion flame on a neutron star surface rotating with 1ms period.

For details, see Zingale et al. (2019), Eiden et al. (2020), and Harpole et al. (2021).

Particle Acceleration in Neutron Star Environments

Get the WarpX code for Electromagnetic PIC

The WarpX code we used is open-source and freely available at the WarpX GitHub.

At LBL, I also collaborated with scientists who worked on electromagnetic modeling for particle accelerator physics. Together, we started a project to model relativistic particle acceleration from magnetic reconnection in rapidly rotating neutron star magnetospheres.

We worked with a talented intern, Eloise Yang, who set up these simulations in the WarpX electromagnetic particle-in-cell (PIC) code. After Eloise left LBL, postdoc Hannah Klion picked up the project. Our project’s first paper explored new advances in electromagnetic PIC solver methods and their utility for simulating charge acceleration from relativistic magnetic reconnection.

Secondary plasmoid reconnection event shown before (top) and after (bottom) initially evolving from a relativistic transverse current sheet.

For details, see Klion et al. (2023).

This research sets the foundation for applying our numerical methods to simulate high-energy particle emission from not only magnetic reconnection but also cosmic ray acceleration in other astrophysical environments.

Back to top

References

Eiden, K., Zingale, M., Harpole, A., Willcox, D., Cavecchi, Y., & Katz, M. P. 2020, ApJ, 894, 6, https://iopscience.iop.org/article/10.3847/1538-4357/ab80bc
Harpole, A., Ford, N. M., Eiden, K., Zingale, M., Willcox, D. E., Cavecchi, Y., & Katz, M. P. 2021, ApJ, 912, 36, https://iopscience.iop.org/article/10.3847/1538-4357/abee87
Klion, H., Jambunathan, R., Rowan, M. E., Yang, E., Willcox, D., Vay, J.-L., Lehe, R., Myers, A., Huebl, A., & Zhang, W. 2023, ApJ, 952, 8, https://iopscience.iop.org/article/10.3847/1538-4357/acd75b
Zingale, M., Eiden, K., Cavecchi, Y., Harpole, A., Bell, J. B., Chang, M., Hawke, I., Katz, M. P., Malone, C. M., Nonaka, A. J., Willcox, D. E., & Zhang, W. 2019, J Phys: Conf Ser, 1225, 012005, https://iopscience.iop.org/article/10.1088/1742-6596/1225/1/012005