Energy Partitioning at the Termination Shock
Abstract
We show new results of a global 3D magnetohydrodynamic (MHD) simulation of the Boston University outer heliosphere model where we used a newly developed approach that distributes the non-adiabatic shock heating among the cold protons, electrons, and pickup-ions (PUIs), while maintaining total energy conservation of all ions (cold protons and PUIs) and electrons. In our previous simulations ( E.S. Bair et al. 2025; B. van der Holst et al. 2026), all non-adiabatic shock heating was channeled to the cold protons, resulting in a too large temperature jump at the termination shock (TS) for the thermal solar wind. Using a new methodology we improved the simulation results with respect to the temperature jump observed at the TS by Voyager 2 (V2) spacecraft. Our simulations approached the observed jump conditions of a factor 10--20 in the cold solar wind temperature V2 measurements, and we obtain improvements in the simulation results relative to the data. Because we directly estimate in this way the distribution of non-adiabatic heating at the TS, we have the opportunity to study the physical process of heating cold plasma and PUIs in the TS. The results show that having almost 100\% non-adiabatic shock heating going towards PUIs at the TS reproduces the jump conditions observed along the V2 trajectory. This information is key to understanding the physical processes that shape the heliosphere. As shown by M. Opher et al. (2020), PUIs significantly change the shape of the heliosphere, for example, the presence of hot PUIs results in a deflated inner heliosheath. Our work provides the architecture of how energy partitioning at shocks in kinetic simulations (J. Giacalone et al. 2021) can be utilized in global MHD models.
Related
- Laboratory Measurements of Energy Partitioning and Anomalous Electron Heating in Magnetized, Perpendicular Collisionless Shocks
- Energy Partitioning in Dust-catalyzed H2 and HD Formation Revealed by Molecular Simulations Considering Nuclear Quantum Effects
- Little Red Dots as Shock-Powered High-Energy Neutrino Sources
- Water-rich sub-Neptunes and rocky super-Earths around different stars: Radii shaped by volatile partitioning, formation, and evolution
- Cosmological Backreaction and the Separate Universe Conjecture: Dynamical volume partitioning of Newtonian cosmologies
BibTeX
@misc{szente2026energy,
title = {{Energy Partitioning at the Termination Shock}},
author = {Judit Szente and Bart van der Holst and Gabor Toth and Merav Opher},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.05267},
eprint = {2609.05267},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.05267},
note = {CC BY-NC-ND 4.0},
}
RIS
TY - UNPB
TI - Energy Partitioning at the Termination Shock
AU - Judit Szente
AU - Bart van der Holst
AU - Gabor Toth
AU - Merav Opher
PY - 2026
DA - 2026/09/04
PB - arXiv
DO - 10.48550/arxiv.2609.05267
UR - https://arxiv.org/abs/2609.05267
AB - We show new results of a global 3D magnetohydrodynamic (MHD) simulation of the Boston University outer heliosphere model where we used a newly developed approach that distributes the non-adiabatic shock heating among the cold protons, electrons, and pickup-ions (PUIs), while maintaining total energy conservation of all ions (cold protons and PUIs) and electrons. In our previous simulations ( E.S. Bair et al. 2025; B. van der Holst et al. 2026), all non-adiabatic shock heating was channeled to the cold protons, resulting in a too large temperature jump at the termination shock (TS) for the thermal solar wind. Using a new methodology we improved the simulation results with respect to the temperature jump observed at the TS by Voyager 2 (V2) spacecraft. Our simulations approached the observed jump conditions of a factor 10--20 in the cold solar wind temperature V2 measurements, and we obtain improvements in the simulation results relative to the data. Because we directly estimate in this way the distribution of non-adiabatic heating at the TS, we have the opportunity to study the physical process of heating cold plasma and PUIs in the TS. The results show that having almost 100\% non-adiabatic shock heating going towards PUIs at the TS reproduces the jump conditions observed along the V2 trajectory. This information is key to understanding the physical processes that shape the heliosphere. As shown by M. Opher et al. (2020), PUIs significantly change the shape of the heliosphere, for example, the presence of hot PUIs results in a deflated inner heliosheath. Our work provides the architecture of how energy partitioning at shocks in kinetic simulations (J. Giacalone et al. 2021) can be utilized in global MHD models.
ER -