Search for: All records

Large-amplitude electrostatic fluctuations are routinely observed by spacecraft upon traversal of collisionless shocks in the heliosphere. Kinetic simulations of shocks have struggled to reproduce the amplitude of such fluctuations, complicating efforts to un- derstand their influence on energy dissipation and shock structure. In this paper, 1D particle-in-cell simulations with realistic proton-to-electron mass ratio are used to show that in cases with upstream electron temperature Te exceeding the ion temperature Ti, the magnitude of the fluctuations increases with the electron plasma-to-cyclotron frequency ratio ωpe/Ωce, reaching realistic values at ωpe/Ωce ≳ 30. The large-amplitude fluctuations in the simulations are shown to be associated with electrostatic solitary structures, such as ion phase-space holes. In the cases where upstream temperature ratio is reversed, the magnitude of the fluctuations remains small. 

Abstract We present the HelioCubed, a high-order magnetohydrodynamic (MHD) code designed for modeling the inner heliosphere. The code is designed to achieve 4th order accuracy both in space and in time. In addition, HelioCubed can perform simulations on mapped grids, such as those based on cubed spheres, which makes it possible to overcome stability limitations caused by the geometrical singularity at the polar axis of a spherical grid, thus enabling substantially larger time steps. HelioCubed has been developed using the high-level Proto library, ensures performance portability across CPU and GPU architectures, and supports back-end implementations, e.g., CUDA, HIP, OpenMP, and MPI. The code is compatible with the HDF5 library, which facilitates seamless data handling for simulations and boundary conditions derived from semi-empirical and MHD models of the solar corona. While presenting the results of preliminary simulations, we demonstrate that our simulations are indeed performed with 4th order of accuracy. Our approach ensures that HelioCubed solves the MHD equations preserving the radial flow to machine round-off error even on cubed-sphere grids. Solar wind simulations are performed using the boundary conditions provided by the Wang–Sheeley–Arge coronal model of the ambient solar wind. It also allows us to to simulate coronal mass ejections using observation-driven flux rope models. These capabilities make HelioCubed a versatile and powerful tool to advance heliophysics research and space weather forecasting. 

The interaction of the solar wind with the local interstellar medium (LISM) spans a wide range of interacting particle populations, energies, and scales. Sophisticated models are required to capture the global picture, interpret near-Earth observations, and ultimately understand the properties of the LISM at distances of thousands of AUs, where the medium is presumed to be unperturbed by this interaction. We present a new extension of our MHD-plasma/kinetic-neutral heliospheric model, implemented within the Multi-Scale Fluid- Kinetic Simulation Suite (MS-FLUKSS). The new model treats singly and doubly charged helium ions, pickup protons, and electrons as separate, self-consistently coupled populations, interacting through six charge exchange processes and photoionization with kinetically treated neutral hydrogen and helium atoms. In this paper, we provide detailed information on the implementation, including new fits for the charge-exchange cross sections, and demonstrate the functionality and performance of the new code 

Abstract The stability of a realistic multicomponent pickup ion (PUI) velocity distribution derived from a global model of neutral atoms in the heliosphere, which treats hydrogen and helium atoms self-consistently and includes equations for electrons and helium ions, is investigated using linear instability analysis and hybrid simulations. Linear instability analysis shows that the excited oblique mirror waves and the parallel/quasi-parallel Alfvén-cyclotron (AC) waves have lower growth rates than those obtained previously by A. Mousavi et al. for the PUI velocity distributions given by J. Heerikhuisen et al. The PUI scattering by each of the two modes alone is studied. In contrast to the previous investigations, our current simulations using the updated realistic distributions indicate that mirror waves alone do not effectively scatter PUIs in pitch angle. Instead, they primarily contribute to reducing the thermal spread anisotropy of the PUIs originating from the neutral solar wind. The unstable AC waves exhibit lower growth rates but higher saturation levels than the mirror waves. Two-dimensional (2D) simulation results show that when all unstable waves are present, the predominant contributor to the fluctuating magnetic field energy is the AC mode. The AC waves quickly scatter the PUIs with pitch angles away from 90^∘toward isotropy, while the PUIs near 90^∘pitch angle maintain a degree of anisotropy within our simulation timeframe. Moreover, several 1D and 2D hybrid simulations with different numbers of particles per cell are performed to examine the impact of numerical noise on PUI scattering. Finally, the implications of these results for the Interstellar Boundary Explorer energetic neutral atom ribbon are discussed. 

Abstract Interplanetary coronal mass ejections (ICMEs) are the primary sources of geomagnetic storms at Earth. The negative out-of-ecliptic component (B_z) of magnetic field in the ICME or its associated sheath region is necessary for it to be geoeffective. For this reason, magnetohydrodynamic simulations of CMEs containing data-constrained flux ropes are more suitable for forecasting their geoeffectiveness as compared to hydrodynamic models of the CME. ICMEs observed in situ by radially aligned spacecraft can provide an important setup to validate the physics-based heliospheric modeling of CMEs. In this work, we use the constant-turn flux rope (CTFR) model to study an ICME that was observed in situ by Solar Orbiter (SolO) and at Earth, when they were in a near-radial alignment. This was a stealth CME that erupted on 2020 April 14 and reached Earth on 2020 April 20 with a weak shock and a smoothly rotating magnetic field signature. We found that the CTFR model was able to reproduce the rotating magnetic field signature at both SolO and Earth with very good accuracy. The simulated ICME arrived 5 hr late at SolO and 5 hr ahead at Earth, when compared to the observed ICME. We compare the propagation of the CME front through the inner heliosphere using synthetic J-maps and those observed in the heliospheric imager data and discuss the role of incorrect ambient solar wind background on kinematics of the simulated CME. This study supports the choice of the CTFR model for reproducing the magnetic field of ICMEs. 

Abstract We introduce the first solar-cycle simulations from our 3D, global MHD-plasma/kinetic-neutrals model, where both hydrogen and helium atoms are treated kinetically, while electrons and helium ions are described as individual fluids. Using Voyager/PWS observations of electron density up to 160 au from the Sun for validation of several different global models, we conclude that the current estimates for the proton density in the local interstellar medium (LISM) need a revision. Our findings indicate that the commonly accepted value of 0.054 cm⁻³may need to be increased to values exceeding 0.07 cm⁻³. We also show how different assumptions regarding the proton velocity distribution function in the outer heliosheath may affect the global solution. A new feature revealed by our simulations is that the helium ion flow may be significantly compressed and heated in the heliotail at heliocentric distances exceeding ∼400 au. Additionally, we identify a Kelvin–Helmholtz instability at the boundary of the slow and fast solar wind in the inner heliosheath, which acts as a driver of turbulence in the heliotail. These results are crucial for inferring the properties of the LISM and of the global heliosphere structure. 

Abstract We present recent advancements in our 3D modeling of the interaction between the solar wind and the local interstellar medium (LISM). The latest model results (Fraternale et al., ApJ, 2023) have raised a question about the electron density of the LISM near the heliopause. We have shown that the presence of helium ions leads to a significant underestimation of this parameter compared to the past simulations and Voyager 1 PWS observations. The latter observations, with over 12 years’ worth of LISM data, offers a robust constraint on our models. Here we present additional simulations in support of the idea that the LISM proton density may need to be revised from approximately 0.054 cm^–3to values around 0.07 cm^–3or higher. Additionally, we have developed and successfully tested a new version of the kinetic code suitable for simulating time-dependent solutions. 

Diffusive shock acceleration requires the production of backstreaming superthermal ions (injection) as a first step. Such ions can be generated in the process of scattering of ions in the superthermal tail off the shock front. Knowledge of the scattering of high-energy ions is essential for matching conditions of upstream and downstream distributions at the shock transition. Here we analyze the generation of backstreaming ions as a function of their initial energy in a model stationary shock and in a similar rippled shock. Rippling substantially enhances ion reflection and the generation of backstreaming ions for slightly and moderately superthermal energies, and thus is capable of ensuring ion injection into a further diffusive shock acceleration process. For high-energy ions, there is almost no difference in the fraction of backstreaming ions produced and the ion distributions between the planar stationary shock and the rippled shock. 

Abstract In a collisionless shock the energy of the directed flow is converted to heating and acceleration of charged particles, and to magnetic compression. In low-Mach number shocks the downstream ion distribution is made of directly transmitted ions. In higher-Mach number shocks ion reflection is important. With the increase of the Mach number, rippling develops, which is expected to affect ion dynamics. Using ion tracing in a model shock front, downstream distributions of ions are analyzed and compared for a planar stationary shock with an overshoot and a similar shock with ripples propagating along the shock front. It is shown that rippling results in the distributions, which are substantially broader and more diffuse in the phase space. Gyrotropization is sped up. Rippling is able to generate backstreaming ions, which are absent in the planar stationary case. 

Abstract Using ion tracing in a model shock front we study heating of thermal (Maxwellian) and superthermal (Vasyliunas–Siscoe) populations of protons, singly charged helium, and alpha particles. It is found that heating of thermal and superthermal populations is different, mainly because of substantially higher ion reflection in the superthermal populations. Accordingly, the temperature increase of initially superthermal populations is substantially higher than that of the thermal ions. Heating per mass decreases with the increase of the mass-to-charge ratio because of the reduced effect of the cross-shock potential and, accordingly, weaker ion reflection. The findings are supported by two-dimensional hybrid simulations.