ST1.4 | Theory and Simulation of Solar System Plasmas
Theory and Simulation of Solar System Plasmas
Convener: Patricio A. Munoz | Co-conveners: Philippa Browning, Maria Elena Innocenti, Giovanni Lapenta, Shangbin Yang
Orals
| Mon, 15 Apr, 10:45–12:30 (CEST)
 
Room 0.16
Posters on site
| Attendance Mon, 15 Apr, 16:15–18:00 (CEST) | Display Mon, 15 Apr, 14:00–18:00
 
Hall X3
Posters virtual
| Attendance Mon, 15 Apr, 14:00–15:45 (CEST) | Display Mon, 15 Apr, 08:30–18:00
 
vHall X3
Orals |
Mon, 10:45
Mon, 16:15
Mon, 14:00
The "Theory and Simulation of Solar System Plasmas" session aims at presenting recent results related to theoretical investigations and/or numerical modelling of plasmas processes and dynamics of heliospheric plasmas: the Sun and its corona, the solar wind or planetary magnetospheres. Those plasma processes include magnetic reconnection, turbulence, shock waves, plasma instabilities, plasma heating, particle acceleration, radiation, etc. Results using all kind of plasma models are welcome, including fluid or kinetic models, from global to local modelling, but also multi-scale approaches. Of particular interest are applications relevant for the interpretation of in-situ measurements and remote observations from current and future space missions such as MMS, Parker Solar Probe, Solar Orbiter and ASO-S.
The focus of this year's session is space plasma turbulence: their origin, effects and links with other plasma processes, from large/fluid to small/kinetic scales.

Orals: Mon, 15 Apr | Room 0.16

Chairpersons: Shangbin Yang, Patricio A. Munoz, Maria Elena Innocenti
10:45–10:50
Kinetics and small-scales
10:50–11:00
|
EGU24-6557
|
ECS
|
On-site presentation
Mihailo Martinović and Kristopher Klein

Linear theory is a well-established theoretical framework proven to accurately characterize instabilities in the solar wind weakly collisional plasma. We aim to describe the statistical properties of linear ion-driven instabilities between 0.3 and 1 au. We analyzed ∼1.5M proton and alpha particle Velocity Distribution Functions (VDFs) observed by Helios I & II, and ~5M VDFs observed by Wind, using Plasma in a Linear Uniform Magnetized Environment (PLUME) dispersion solver to calculate growth rate, frequency, wavevector, and the power emitted or absorbed by each VDF component. The descriptive statistical analysis shows that the stability of the solar wind is primarily determined by the collisional processing, rather than the distance from the Sun. We use this data set to train Stability Analysis Vitalizing Instability Classification (SAVIC) Machine Learning algorithm capable to classify the predicted unstable modes into physically meaningful “textbook” types. This method enables us to map the instability properties in multi-dimensional phase space. The proton-core-induced Ion Cyclotron (IC) mode dominates the collisionally young solar wind, while the alpha population plays more important role in the wave energy dynamics of the older wind. We also demonstrate that the proton beam population is not affected by the collisions and has the core-beam drift as the main source of free energy that determines its overall behavior in the solar wind. SAVIC code used here is publicly available for the community.

How to cite: Martinović, M. and Klein, K.: Overview of Ion-Driven Instabilities in the Inner Heliosphere, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-6557, https://doi.org/10.5194/egusphere-egu24-6557, 2024.

11:00–11:10
|
EGU24-17196
|
On-site presentation
Jinsong Zhao

How the solar wind heat flux is constrained and how strahl electrons are scattered are two fundamental problems in the study of the electron dynamics in the solar wind. Recently, much attention has been paid to the role of the electron heat flux instability on these two problems. We have performed the instability analyses for the electron heat flux instability in the near-Sun solar wind where the plasma beta is much lower than one. We found that lower-hybrid waves and oblique Alfvén waves can be primarily triggered by the electron heat flux in the low-beta plasma environment, and we also explored the wave-particle interactions in each type instability through the energy transfer rate method. Moreover, using PSP observations, we will show the possible observation signature for the lower-hybrid waves in the near-Sun solar wind.

How to cite: Zhao, J.: Electron heat flux instability in the near-Sun solar wind, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-17196, https://doi.org/10.5194/egusphere-egu24-17196, 2024.

11:10–11:20
|
EGU24-2175
|
On-site presentation
Seiji Zenitani and Tsunehiko N Kato

We propose a family of numerical solvers for the nonrelativistic Newton-Lorentz equation in particle-in-cell (PIC) simulation. The new solvers extend a popular 4-step procedure, which has second-order accuracy in time, in several ways. First, we repeat the 4-step procedure n cycles, using an n-times smaller timestep (Delta_t/n). To speed up the calculation, we derive a polynomial formula for an arbitrary cycling number n, based on our earlier work (Zenitani & Kato 2020, Comput. Phys. Commun.). Second, prior to the 4-step procedure, we apply Boris-type gyrophase corrections to the electromagnetic field. In addition to the magnetic field, we amplify the electric field in an anisotropic manner to achieve higher-order (N=2,4,6... th order) accuracy. Finally, we construct a hybrid solver of the n-cycle solver and the Nth-order solver. We call it the hyper Boris solver. The (n,N) hyper Boris solver gives a numerical error of ~ (Delta_t/n)N at affordable computational cost.

How to cite: Zenitani, S. and Kato, T. N.: Hyper Boris integrators for particle-in-cell simulation, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-2175, https://doi.org/10.5194/egusphere-egu24-2175, 2024.

11:20–11:30
|
EGU24-7222
|
ECS
|
Highlight
|
On-site presentation
Hayato Higuchi, Juan Pedersen, and Akimasa Yoshikawa

The space plasma environment, extending from the Sun to the Earth, includes regions of frozen conditions, zones of anomalous resistance caused by electromagnetic turbulence, interconnected regions characterized by weakly ionized gas systems in strong magnetic fields, coupled neutral-atmosphere chemical processes, and pure neutral-atmosphere collision systems. Owing to their complex interactions, an inclusive understanding and forecasting of the space environment remains an elusive goal, even with the advancements in high-performance instrumentation and in-situ observation of satellites. Therefore, it is imperative to develop space plasma simulations capable of providing comprehensive insights, ranging from local spatial domains to the global schematic.

Historically, the development of space plasma simulations has been constrained by computational time, memory capacity, and data storage limitations, resolving complex phenomena with restricted physics at local space scales.

Recently, as quantum computing hardware has advanced, quantum algorithms have proven to benefit exponential speedups. There has been a focus on the practical applications of quantum computing in finance, chemistry, fluids, and a variety of fields (e.g., Bouland et al.,[2020], Cao et al.,[2019], Egger et al.,[2020] and Budinski, [2022]). Then a quantum algorithm for the collisionless Boltzmann equation using the discrete velocity method was developed by (Todorova and Steijl, [2020]).

We have developed a quantum algorithm for the six-dimensional Boltzmann-Maxwell equations for collisionless plasmas with reference to (Higuchi, et al.,[2023]). We applied this methodology to propose a quantum approach that predicts kinetic and multiscale plasma dynamics.

In this presentation, we will introduce a quantum algorithm for performing super large- scale plasma simulations in a quantum computer and estimate the performance required by this task for future quantum error tolerant large-scale quantum computers. We will also discuss the prospects for applications that can be expected in our field, based on trends in the development of quantum computer hardware and software.

How to cite: Higuchi, H., Pedersen, J., and Yoshikawa, A.: Quantum Computing for Future Super Large-Scale Plasma Simulations: A Novel Approach for Simulating the Vlasov-Maxwell System, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-7222, https://doi.org/10.5194/egusphere-egu24-7222, 2024.

11:30–11:40
|
EGU24-6152
|
ECS
|
Virtual presentation
Sebastian Echeverria Veas, Pablo Moya, Marian Lazar, Stefaan Poedts, and Felipe Asenjo

Multiscale modeling of expanding plasmas is crucial for understanding the dynamics and evolution of various astrophysical plasma systems. In this context, the Expanding Box Model (EBM) was used to add the expansion into the kinetic equations, allowing us to describe plasma physics in a new system of reference non-expanding and co-moving with the plasma. This system allows us to maintain a constant volume through non-inertial forces, and its interpretation is fundamental to describing plasma physics.

We have employed the EBM formalism to incorporate the expanding properties of the system into the plasma dynamics, which mainly affects transverse coordinates (i.e., y y/o z). Coordinate transformations were introduced within the co-moving frame system to obtain the modified Vlasov equation. Our main goal is to develop a plasma physics theory through a novel first principle description in the expanding frame, which is fundamentally based on the (collisionless) Vlasov equation for the evolution of the velocity distribution functions. Based on this, the expanding moments, such as the continuity, momentum, and energy equations, can then be derived, and an MHD model of the plasma expansion can be developed. Finally, coupling the obtained moments and Maxwell equations, a CGL-like plasma description is developed in the EB frame to study the evolution of macroscopic quantities (temperature, magnetic field, parallel beta, and anisotropy).

Our results show the expansion affecting the kinetic and fluid equations through non-inertial and fictitious forces in the transverse directions, which contain all the information related to the expansion. These are thus reflected by the equations derived for the expanding moments of the distribution function, including density, bulk (drift) velocity, and pressure (or temperature). Furthermore, we developed an ideal expanding-MHD model based on these modified moments, providing a new interpretation and comparison with the existing results when expansion is considered. The EBM modifies the conservative form of the two adiabatic invariants in the CGL approximation. Equations are solved for radially decreasing magnetic fields and density profiles to study the relations between plasma parallel beta and anisotropy within the expansion. 

How to cite: Echeverria Veas, S., Moya, P., Lazar, M., Poedts, S., and Asenjo, F.: First Principle Description of Plasma Expansion using the Expanding Box Model. Implications for CGL Theory in the Solar Wind., EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-6152, https://doi.org/10.5194/egusphere-egu24-6152, 2024.

MHD and large-scales
11:40–11:50
|
EGU24-6886
|
solicited
|
On-site presentation
Liping Yang, Jiansen He, Daniel Verscharen, Hui Li, Trevor A. Bowen, Stuart D. Bale, Honghong Wu, Wenya Li, Ying Wang, Lei Zhang, Xueshang Feng, and Ziqi Wu

Imbalanced Alfvénic turbulence is a universal process playing a crucial role in energy transfer in space, astrophysical, and laboratory plasmas. A funda-
mental and long-lasting question about the imbalanced Alfvénic turbulence is how and through which mechanism the energy transfers between scales. Here, we show that the energy transfer of imbalanced Alfvénic turbulence is completed by coherent interactions between Alfvén waves and co-propagating anomalous fluctuations. These anomalous fluctuations are generated by nonlinear couplings instead of linear reflection. We also reveal that the energy transfer of the waves and the anomalous fluctuations is carried out mainly through local-scale and large-scale nonlinear interactions, respectively, responsible for their bifurcated power-law spectra. This work unveils the energy transfer physics of imbalanced Alfvénic turbulence, and advances the understanding of imbalanced Alfvénic turbulence observed by Parker Solar Probe in the inner heliosphere.

How to cite: Yang, L., He, J., Verscharen, D., Li, H., Bowen, T. A., Bale, S. D., Wu, H., Li, W., Wang, Y., Zhang, L., Feng, X., and Wu, Z.: Energy transfer of imbalanced Alfvénic turbulence in the heliosphere, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-6886, https://doi.org/10.5194/egusphere-egu24-6886, 2024.

11:50–12:00
|
EGU24-11324
|
On-site presentation
Bart van der Holst, Enrico Landi, and Tamas Gombosi

We present a multi-ion generalization of the AWSoM model, a 3D global magnetohydrodynamic (MHD) solar corona and inner heliosphere model with incompressible turbulence. The AWSoM model with charge states is further extended to include temperatures for the various heavy ion species. The coronal heating is addressed via outward propagating low-frequency Alfven waves that are partially reflected by Alfven speed gradients. The nonlinear interaction of these counter-propagating waves results in turbulent energy cascade. To apportion the cascaded turbulent energies to the electron and ion temperatures, we employ the results of the theories of linear wave damping and nonlinear stochastic heating. This heat partitioning results in a mass proportional heating among ions.

How to cite: van der Holst, B., Landi, E., and Gombosi, T.: AWSoM solar wind model with charge states and ion temperatures, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-11324, https://doi.org/10.5194/egusphere-egu24-11324, 2024.

12:00–12:10
|
EGU24-3315
|
On-site presentation
Gwangson Choe, Minseon Lee, and Sibaek Yi

It is generally accepted that solar prominences, also known as filaments, are formed through thermal condensation instability from hot coronal plasmas. Within solar prominences (filaments), counter-streaming subscale flows are frequently observed. It is quite surprising that thermal instability in magnetized plasmas with such shear flows has never been earnestly studied despite extensive research on thermal instability in various astrophysical contexts. In this paper, we have investigated this unexplored territory of thermal instability and have unexpectedly gained hints on why solar filaments are filamentary.

We have performed linear stability analysis of magnetized plasmas with shear flows within the framework of magnetohydrodynamics (MHD), incorporating radiative cooling, phenomenological plasma heating, and anisotropic thermal conduction. Our approach formulates an eigenvalue problem, which we solve numerically to derive eigenfrequencies and eigenfunctions.

Our findings reveal that, for shear speeds less than the Alfven speed of the background plasma, the dominant mode corresponds to an isobaric thermal condensation mode. Most remarkably, the eigenfunctions associated with this mode display a distinctive, discrete structure resembling delta functions, especially when the shear velocity in the k-direction exceeds 10−5 of the Alfven speed. We identify that these delta function-like spikes coincide with the zeroes of the coefficients of the second-order derivative terms in the differential equation of the eigenvalue problem.

In contrast, for shear speeds exceeding the Alfven speed (a rare occurrence in reality), we observe an isentropic Kelvin-Helmholtz instability, incompatible with thermal condensation.

Our investigation underscores that any non-uniform velocity field with a magnitude surpassing 10−5 of the Alfven speed triggers the discrete eigenfunction characteristic of the condensation mode. Consequently, filamentary condensation at discrete layers or threads, emerges as a natural and universal process whenever thermal condensation instability arises in magnetized plasmas with shear flows.

How to cite: Choe, G., Lee, M., and Yi, S.: A Clue to the Filamentary Nature of Solar Filaments, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-3315, https://doi.org/10.5194/egusphere-egu24-3315, 2024.

12:10–12:20
|
EGU24-2303
|
ECS
|
Virtual presentation
Mayank Kumar, Kris Murawski, Blazej Kuźma, Luis Kadowaki, Emilia Kilpua, Stephan Poedts, and Robertus Erdelyi

Context. The heating of the solar chromosphere, the associated plasma outflows, and the origin of the solar wind are key issues in heliophysics. In this paper, we provide a new perspective on their connection to the propagation and dissipation of waves generated by solar granulation.


Aims. The primary objective of this paper is to conduct 2.5-D numerical simulations of the partially ionized lower solar atmosphere, investigating the propagation and dissipation of granulation-generated waves in the context of plasma outflows and the related heating of the chromosphere, which is due to ion-neutral collisions.


Methods. We use the JOint ANalytical and Numerical Approach (JOANNA) code to solve the two-fluid model equations. We take into account partially ionized hydrogen plasma composed of ions (protons) and neutrals (H atoms), which are coupled via ion-neutral collisions. We focus on a quiet region of the solar chromosphere which is gravitationally stratified and magnetically constrained by an initially set magnetic Plasma flows and solar chromosphere heating arcade. Solar convection situated beneath the photosphere is the main source of these waves that are propagating through the simulated solar atmosphere.


Results. The numerical results obtained in our study reveal an important process in the lower solar atmosphere. The naturally evolving convection generates waves and a portion of the wave energy is dissipated due to ion-neutral collisions in the solar chromosphere. This dissipation of waves, in turn, leads to the release of thermal energy, resulting in the heating of the solar atmosphere. This phenomenon is also associated with upward-directed plasma flows, which may play a role in the formation of the solar wind. Furthermore, our analysis of the dominant wave periods in the various layers
of the solar atmosphere closely aligns with observational data from Wiśniewska et al. (2016) and Kayshap et al. (2018). This alignment underscores the crucial role ion-neutral collisions play in facilitating the energy release process, shedding light on the intricate dynamics of the solar atmosphere.


Conclusions. Based on our numerical simulations, we can draw the following conclusions: The dissipation of waves in the two-fluid plasma caused by ion-neutral collisions in the two-fluid plasma model leads to plasma outflows and increased heating in the chromosphere. These plasma outflows may play a role in the generation of the solar wind and the accompanying heating.

How to cite: Kumar, M., Murawski, K., Kuźma, B., Kadowaki, L., Kilpua, E., Poedts, S., and Erdelyi, R.: Solar granulation-generated chromospheric heating and plasma outflows in two-fluid magnetic arcade, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-2303, https://doi.org/10.5194/egusphere-egu24-2303, 2024.

12:20–12:30
|
EGU24-2367
|
On-site presentation
Xiaoshuai Zhu

Modeling the three-dimensional (3D) magnetic fields of the solar active region in multiple layers is very important. The main approach is to extrapolate the magnetic field from magnetograms measured in the photosphere. A basic assumption of the modeling in the past several decades was to completely neglect all plasma effects and to perform the so-called force-free field (FFF) extrapolations. While the force-free assumption is well justified in the solar corona, it is not the case in the photosphere and chromosphere. A magneto-hydro-static (MHS) equilibrium which takes into account plasma forces, such as pressure gradient and gravitational force, is considered to be more appropriate to describe the lower atmosphere and has been developing rapidly during the past several years. In this talk, I am going to review various MHS modeling methods, including tests of these methods with known reference models and applications to real data. Recent developments of the MHS modeling will also be presented.

How to cite: Zhu, X.: On the current state of the development of the magneto-hydro-static extrapolation method, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-2367, https://doi.org/10.5194/egusphere-egu24-2367, 2024.

Posters on site: Mon, 15 Apr, 16:15–18:00 | Hall X3

Display time: Mon, 15 Apr 14:00–Mon, 15 Apr 18:00
X3.1
|
EGU24-4321
Liping Yang, Jiansen He, Xueshang Feng, Hui Li, Fan Guo, Hui Tian, and Fang Shen

Alfven waves contribute significantly to the solar coronal heating, the solar wind acceleration, as well as Alfv\'enic turbulence formation. As a universal process, magnetic reconnection has long been credited as a potentially crucial source of Alfven waves, but how magnetic reconnection trigger Alfven waves remains elusive. Here, with simulations of three-dimensional bursty interchange magnetic reconnection in the solar corona, for the first time, we find that Alfven waves are spontaneously excited in the reconnection sheet and propagate bi-directionally even along the unreconnected magnetic fields. The enhanced total pressure inherently carried by flux ropes gives kicks to the magnetic fields, and the nearly same propagation speed of the flux ropes and the kicks makes the kicks growing into the observed Alfven waves, which have large amplitudes and high frequencies, carrying substantial energy for heating the quiet corona and accelerating the solar wind. Our findings demonstrate that Alfven waves are natural products of three-dimensional intermittent magnetic reconnection, bringing its fundamental significance for energy release, transport, and conversion occurring in the plasma system.

How to cite: Yang, L., He, J., Feng, X., Li, H., Guo, F., Tian, H., and Shen, F.: Spontaneous Generation of Alfven Waves by Bursty Interchange Magnetic Reconnection in the Solar Corona, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-4321, https://doi.org/10.5194/egusphere-egu24-4321, 2024.

X3.2
|
EGU24-11114
|
ECS
Rong Lin, Giovanni Lapenta, and Jiansen He

Recent observations of the interplanetary magnetic field have enlarged the magnetic field topology family in the inner heliosphere, and raises interest in the relationship between those topologies and the evolution of solar wind. For example, switchbacks that are kinks of magnetic field lines accompanied by velocity spikes and enhancement of the proton parallel temperature, have free energy both of the field and of the particles, and thus have potential for a series of plasma instabilities. It is favorable to extract elementary topologies from complicated structures like switchbacks before we look further in. Based on Particle-in-Cell simulations, we discuss from elementary topologies on the effect they have on kinetic processes of the solar wind plasma.

How to cite: Lin, R., Lapenta, G., and He, J.: How will interplanetary magnetic field topology modify solar wind kinetics?, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-11114, https://doi.org/10.5194/egusphere-egu24-11114, 2024.

X3.3
|
EGU24-13696
|
Highlight
An Automatic Approach for Grouping Sunspots and Calculating Relative Sunspot Number on SDO/HMI Continuum Images
(withdrawn)
Cui Zhao and Shangbin Yang
X3.4
|
EGU24-18356
|
ECS
Bing Ma, Ling Chen, and De-Jin Wu

    Solar and interplanetary radio bursts are important phenomena to reflect the electron acceleration and kinetic process in the corona and space plasmas. The near-Sun radio observation by Parker Solar Probe (PSP) provides an good chance to make the approach and in situ measurements for the emission source of radio bursts. According to the radio dynamic spectrum, we found that a large number of weak radio bursts with higher cutoff frequency flo can only be detected by PSP. These bursts have several obvious characteristics: (1) short duration (1-5 min); (2) narrow frequency band (0.5-15 MHz); (3) weak peak intensity (~ 10-15 V2/Hz). Their relative frequency drift rate decreases from > 0.01 s-1 to < 0.01 s-1 implies that they are not the typical type III and II radio bursts. Based on the plasma empirical models and the data from in situ detection by PSP, the fitted models indicate that the radiation of these bursts might be generated by the electron cyclotron maser emission in the acceleration region of solar wind (1.1-10 RS). The spectral characteristics of these bursts manifest that these bursts come from small-scale emission source, which have experienced strong kinetic evolution. We propose that these weak bursts are possibly the solitary wave radiation (SWR) generated from electron cyclotron maser instability with the energetic electrons trapped and accelerated by solitary kinetic Alfvén wave (SKAW) in the close magnetic structure. The decay of SKAW can well explain the deceleration of the emitting sources and the short duration.

How to cite: Ma, B., Chen, L., and Wu, D.-J.: Emission Mechanism of radio bursts from the Solar Wind Acceleration Region observed by Parker Solar Probe in near-Sun plasmas, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-18356, https://doi.org/10.5194/egusphere-egu24-18356, 2024.

X3.5
|
EGU24-18387
|
ECS
Minseon Lee, Gwangson Choe, and Sibaek Yi

A nonlinear force-free field is solely determined by the normal components of magnetic field and current density on the entire boundary of the domain. Methods using three components of magnetic field suffer from either overspecification of boundary conditions and/or a nonzero divergence B problem. A vector potential formulation eliminates the latter issue, yet poses difficulties in imposing normal components of both magnetic field and current density on the boundary. This challenge arises due to the inability to fix all three components of the vector potential at the boundary while the vector potential within the interior domain undergoes iterative changes.

This paper explores four distinct boundary treatment approaches within the vector potential formulation. In two methods, the normal component of the vector potential is adjusted to satisfy the prescribed boundary-normal component of current density at each iteration, while the tangential components remain fixed throughout. In the other two methods, the tangential components of the vector potential are modified at each iteration, while the normal component remains fixed. Each group comprises both first-order and second-order boundary treatments. We conduct a comparative analysis of these methods against our established poloidal-toroidal formulation code and the optimization code in Solar Software.

While none of the four methods outperform our poloidal-toroidal formulation code, they demonstrate comparable or superior performance compared to the latter. This research is regarded to provide insights into optimizing boundary conditions for data-driven simulations using vector potentials. 

How to cite: Lee, M., Choe, G., and Yi, S.: Implementation of Boundary Conditions for Nonlinear Force-Free Field Computations Using Vector Potentials, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-18387, https://doi.org/10.5194/egusphere-egu24-18387, 2024.

X3.6
|
EGU24-7217
Shaaban Mohammed Shaaban Hamd, Rodrigo A. López, Marian Lazar, and Stefaan Poedts

It is believed that the suprathermal populations of electrons and protons in the solar wind and terrestrial magnetosphere (with energies exceeding those of the bulk population up to several keV) can offer key answers to many major problems, such as the origin of energetic solar particles from interplanetary shocks, particle acceleration by the energy dissipation of the small-scale wave fluctuations, but also a certain level of kinetic turbulence which can explain the non-equilibrium quasi-stable states of space plasmas. Due to their low density and relatively high energy, these populations are collisionless, which means that their dynamics is governed by wave-particle interactions, especially resonant Landau or cyclotron interactions. We present here a series of recent results that demonstrate that the suprathermal populations, whether in the form of a less-drifting halo or the field-aligned beams/strahls, are mainly responsible for the resonant excitation of kinetic wave instabilities. The in-situ observations confirm both the electromagnetic fluctuations triggered by temperature anisotropy and the predominantly electrostatic excitations induced for instance by electron beams.

How to cite: Hamd, S. M. S., López, R. A., Lazar, M., and Poedts, S.: Resonant Wave Excitations of Suprathermal Populations in the Solar Wind and Terrestrial Magnetosphere, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-7217, https://doi.org/10.5194/egusphere-egu24-7217, 2024.

X3.7
|
EGU24-19524
|
ECS
Studying a Non-Collisional Plasma with Superstatistics considerations.
(withdrawn)
Abiam Tamburrini, Pablo Moya, and Sergio Davis
X3.8
|
EGU24-8120
Daniel Verscharen, Alfredo Micera, Maria Elena Innocenti, Jesse Coburn, Elisabetta Boella, Viviane Pierrard, Jingting Liu, Christopher J. Owen, Georgios Nicolaou, and Kristopher G. Klein

The electrons in the solar wind often exhibit non-equilibrium velocity distribution functions. Observed non-equilibrium electron features in the inner heliosphere include a field-aligned beam (called "strahl"), a suprathermal halo population, a sunward deficit in the distribution, and temperature anisotropy. These features are the result of a complex interplay between global expansion effects, collisions, and local interactions between the particles and the electromagnetic fields. Global effects create, for example, the strahl via the mirror force in the decreasing magnetic field and the sunward deficit via reflection effects in the interplanetary electrostatic potential. Local wave-particle interactions such as instabilities and wave damping change the shape of these signatures and thus the overall properties and moments of the electron distribution.

We discuss the formation of the relevant features in the electron distribution and analyse their impact on the linear stability of whistler waves in the inner heliosphere. We then present results from our numerical ALPS code that is capable of evaluating the linear stability of plasma with arbitrary background distributions. With results from our ALPS code, we show that the strahl-core-deficit configuration near the Sun drives oblique whistler waves unstable. However, it leads to enhanced damping of parallel whistler waves compared to a Maxwellian configuration. As the distribution evolves, the sunward deficit fills with electrons, at which point the plasma becomes unstable and drives parallel whistler waves. Our results highlight the need to treat electrons statistically as a globally inhomogeneous plasma component and to account for the detailed shape of their distribution in the evaluation of the plasma's linear stability.

How to cite: Verscharen, D., Micera, A., Innocenti, M. E., Coburn, J., Boella, E., Pierrard, V., Liu, J., Owen, C. J., Nicolaou, G., and Klein, K. G.: Statistical mechanics of the electrons in the solar wind: stability and instability of whistler waves in the inner heliosphere, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-8120, https://doi.org/10.5194/egusphere-egu24-8120, 2024.

X3.9
|
EGU24-8536
|
ECS
Nicolás Francisco Sepúlveda, Pablo S. Moya, Rodrigo A. López, and Daniel Verscharen

The kinetic Alfvén wave (KAW) is an extension of the Alfvén mode to the kinetic scales that propagates at quasi-perpendicular angles with respect to the background magnetic field in a magnetized plasma. The KAW is characterized by a right-handed polarization in the plane orthogonal to the background magnetic field. This allows the KAW to resonate with electrons, as contrary to the electromagnetic ion-cyclotron (EMIC) wave, which resonates with positive ions due to it’s its left-handed polarization. The EMIC mode is also a kinetic extension of the classical Alfvén wave, and propagates at quasi-parallel wavenormal angles. Due to the relevance and overall presence of both of these modes in space plasmas, understanding the nature of the transition from the EMIC mode to the KAW is a matter of great interest for the study of the micro-scale physics of these systems.

The transition from left-hand to right-hand polarized Alfvén waves depends on the wavenumber, plasma beta, temperature anisotropy, and ion composition of the plasma. Along with the temperature anisotropy, the electron-to-proton temperature ratio Te/Tp is of great relevance for the characterization of the thermal properties of a plasma. This ratio varies significantly between different space plasma environments. Thus, studying how variations on this ratio affect the polarization properties of electromagnetic waves becomes of high relevance highly relevant
for our understanding of the dynamics of space plasmas.

In this work, we present an extensive study on the effect of the thermal properties of electrons on the behaviour and characteristics of Alfvénic waves in fully kinetic linear theory, as well as on the transition from EMIC to KAW. We show that the temperature ratio Te/Tp has strong and non-trivial effects on the polarization of the Alfvénic modes, especially at kinetic scales (kρL>1, where k=k sinθ and ρL= cs/Ωp, with cs the plasma sound speed and Ωp the proton’s gyrofrequency) and βep>0.5, where βs=8πnkTs/B2 is the ratio between thermal and magnetic pressure. We conclude that electron inertia plays an important role in the kinetic scale physics of the KAW in the warm plasma regime, and thus cannot be excluded in hybrid models for computer simulations.

How to cite: Sepúlveda, N. F., Moya, P. S., López, R. A., and Verscharen, D.: The effect of the thermal properties of electrons on the dispersion of kinetic Alfvén waves, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-8536, https://doi.org/10.5194/egusphere-egu24-8536, 2024.

X3.10
|
EGU24-8669
Relative Alpha in the MHD with open magnetic field boundary and its application to the solar activity
(withdrawn)
Shangbin Yang, Joerg Buchner, Jean Santos, Jan Skala, and Hongqi Zhang
X3.11
|
EGU24-13542
|
ECS
Quan Wang, Shangbin Yang, Mei Zhang, Xiao Yang, and Xiaoshuai Zhu

Magnetic helicity is an important concept in solar physics, with a number of theoretical statements pointing out the important role of magnetic helicity in solar flares and coronal mass ejections (CMEs). Here we construct a sample of 47 solar flares, which contains 18 no-CME-associated confined flares and 29 CME-associated eruptive flares. We calculate the change ratios of magnetic helicity and magnetic free energy before and after these 47 flares. Our calculations show that the change ratios of magnetic helicity and magnetic free energy show distinct different distributions in confined flares and eruptive flares. The median value of the change ratios of magnetic helicity in confined flares is −0.8%, while this number is −14.5% for eruptive flares. For the magnetic free energy, the median value of the change ratios is −4.3% for confined flares, whereas this number is −14.6% for eruptive flares. This statistical result, using observational data, is well consistent with the theoretical understandings that magnetic helicity is approximately conserved in the magnetic reconnection, as shown by confined flares, and the CMEs take away magnetic helicity from the corona, as shown by eruptive flares.

How to cite: Wang, Q., Yang, S., Zhang, M., Yang, X., and Zhu, X.:  Change Ratios of Magnetic Helicity and Magnetic Free Energy During MajorSolar Flares , EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-13542, https://doi.org/10.5194/egusphere-egu24-13542, 2024.

X3.12
|
EGU24-15727
|
ECS
Isabella Kraus and Philippe Bourdin

The exact coronal heating mechanism remains a riddle, but magnetically active regions are known to host coronal loops with extreme-UV emission. But also at much smaller sizes, up to 10 Mm, there are bipolar regions that can be associated with UV emission in coronal bright points (CBPs). We study the statistical properties of CBPs with continuous data from the SDO spacecraft to track the lifetime of CBPs. We use their tracking data to verify that the lower corona co-rotates with the photosphere. In a next step, we aim to reproduce an isolated CBP in a 3D magneto-hydrodynamic (MHD) simulation. To this end, we need observational data to drive the simulation from both, SDO/HMI and Hinode/NFI. As these instruments feature different resolutions and field-of-views, they also detect different levels of small-scale and large-scale magnetic structures in the photosphere. As we know, these magnetic patches are advected from photospheric horizontal motions and create the necessary Poynting flux at the base of the corona. Combining these data from two very different instruments is a task that needs careful overlaying, so that not artificial effects would appear in the MHD simulation. We use a multi-scale overlaying method to enlarge the field-of-view of Hinode with SDO data, to drive the simulation with consistent photospheric magnetic fields. The bottom and top boundaries are fully closed for any mass and heat flows. The output of the simulation will allow us to compute synthetic UV emission maps that we may compare directly to SDO and Hinode observations. With our model we test if the field-line braiding mechanism is sufficient to heat a CPB to the required temperature. We find that loop-like CBPs usually originate from bipolar regions. Weaker magnetic polarities produce fainter and hence cooler CBPs. And the same time, we find that lifetimes of typical CPBs are easily more than 6 hours. This supports the theory that the heating of CBP is mainly based on magnetic energy dissipation through a relatively steady and slow magnetic reconnection process.

How to cite: Kraus, I. and Bourdin, P.: Coronal Bright Points observed in the corona and photosphere with SDO and Hinode, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-15727, https://doi.org/10.5194/egusphere-egu24-15727, 2024.

X3.13
|
EGU24-17085
|
ECS
Ahmed Houeibib, Filippo Pantellini, and Léa Griton

We simulate the propagation of relativistic test particles within the field of a 3D MHD simulation of the solar wind. The adiabatic ideal MHD equations are integrated numerically using the MPI-AMRVAC code. Test particles are initialized within the MHD simulation grid and advanced in time according to the guiding center equations, we employ a third-order accurate time prediction-correction method from Mignone et al. 2023 for integration. We also include possibility of diffusion in velocity space based on a particle-turbulence mean free path λ∥ along the magnetic field line. One of the first results where we consider 81 keV electrons injected at 0.139 AU heliocentric distance and mean free path λ∥ = 0.5 AU is in good qualitative agreement with measurements at 1 AU.

How to cite: Houeibib, A., Pantellini, F., and Griton, L.: Exploring Solar Energetic Particles Transport in the Inner Heliosphere, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-17085, https://doi.org/10.5194/egusphere-egu24-17085, 2024.

X3.14
|
EGU24-18011
|
ECS
Jingting Liu, Daniel Verscharen, Jesse Coburn, Jeffersson Agudelo, Kai Germaschewski, Hamish Reid, Georgios Nicolaou, and Christopher Owen

Langmuir waves are frequently detected in the solar wind and affect the energetics of the plasma electrons. Previous observational studies have found Langmuir waves associated with magnetic holes in the solar wind. In our work, we aim to understand the connection between magnetic holes and these waves.

The Langmuir instability is a well-known consequence of electron beams in plasmas. We use particle-in-cell (PIC) simulations of a collisionless electron-ion plasma in a magnetic hole structure. We study the triggering of Langmuir waves by inhomogeneous beam instabilities in this configuration.

Our simulations reveal patterns that suggest that injecting a beam into the magnetic hole has the potential to create spatially inhomogeneous Langmuir wave packets by instability. We support our PIC simulations with linear stability calculations and discuss the conditions required for magnetic holes to emit Langmuir waves through our proposed mechanism. Our simulation results will be used in comparative studies with observational data of Langmuir-wave emitting magnetic holes.

 

 

How to cite: Liu, J., Verscharen, D., Coburn, J., Agudelo, J., Germaschewski, K., Reid, H., Nicolaou, G., and Owen, C.: Simulations of Langmuir-wave emission by magnetic holes in the solar wind, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-18011, https://doi.org/10.5194/egusphere-egu24-18011, 2024.

X3.15
|
EGU24-8372
|
ECS
Hannah Theresa Rüdisser, Andreas J. Weiss, Christian Möstl, Ute V. Amerstorfer, Emma E. Davies, and Eva Weiler

A special observational signature in ICMEs that has puzzled researchers for a long time are so-called ”back regions” or ”magnetic-cloud-like (MCL)” parts of ICMEs which follow after the main flux rope rotation has passed the observer. These regions, occurring after the main flux rope rotation, display a peculiar behaviour where the magnetic field remains elevated with minimal rotation. Two proposed explanations involve magnetic reconnection during CME propagation or a purely geometric effect as the observer traverses the flux rope . 

We investigate in detail whether MCLs can be explained by an effect of the trajectory of an observer through a 3D expanding magnetic flux rope, without invoking magnetic reconnection as an explanation for those signatures. To this end, we employ the 3D coronal rope ejection (3DCORE) model, which has proven its ability to fit in situ magnetic fields of CME flux ropes. 

The model assumes an empirically motivated torus-like flux rope structure that expands self-similarly within the heliosphere, is influenced by a simplified interaction with the solar wind environment, and carries along an embedded analytical magnetic field. Developed for fitting, generating synthetic signatures, comparing models to observations, and analysing results, 3DCOREweb accelerates the determination of physical parameters, fostering research on the global magnetic structure of CMEs. Additionally, the interface aids in advancing our understanding of magnetic flux rope signatures arising from diverse 3D trajectories through CMEs.

How to cite: Rüdisser, H. T., Weiss, A. J., Möstl, C., Amerstorfer, U. V., Davies, E. E., and Weiler, E.: Understanding the effects of spacecraft trajectories through solar coronal mass ejection flux ropes using 3DCOREweb, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-8372, https://doi.org/10.5194/egusphere-egu24-8372, 2024.

X3.16
|
EGU24-20270
|
ECS
|
Luca Barbieri, Lapo Casetti, Andrea Verdini, Simone Landi, Emanuele Papini, and Pierfrancesco Di Cintio

The temperature of the solar atmosphere increases from thousands to millions of degrees moving from the lower layer (chromosphere) to the outermost one (corona), while the density drops accordingly. The mechanism behind this phenomenon, known as a temperature inversion, is still unknown. In this work, we model a coronal loop as a collisionless plasma confined in a semicircular tube that is subject to the Sun's gravity and in thermal contact with a fully collisional chromosphere behaving as a thermostat at the loop's feet. By using kinetic N-particle simulations and analytical calculations, we show that rapid, intermittent, and short-lived heating events in the chromosphere drive the coronal plasma towards a non-equilibrium stationary state. The latter is characterized by suprathermal tails in the particles velocity distribution functions, exhibiting temperature and density profiles strikingly similar to those observed in the atmosphere of the Sun. These results suggest that a million-Kelvin solar corona can be produced without the local deposition of heat in the upper layer of the atmosphere that is typically assumed by standard approaches. We find that suprathermal distribution functions in the corona are self-consistently produced instead of postulated a priori, in contrast to classical kinetic models based on a velocity filtration mechanism.

How to cite: Barbieri, L., Casetti, L., Verdini, A., Landi, S., Papini, E., and Di Cintio, P.: Temperature inversion in a gravitationally bound plasma: Case of the solar corona, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-20270, https://doi.org/10.5194/egusphere-egu24-20270, 2024.

X3.17
|
EGU24-18200
|
ECS
Bahaeddine Gannouni, Victor Réville, and Alexis Rouillard

Our research delves into solar wind's mesoscale features known as microstreams—periods of heightened speed and temperature lasting hours. Initially observed in Helios and Ulysses data, they're now prevalent in the 'young' solar wind by Parker Solar Probe and Solar Orbiter. Recent findings unveil microstreams' carriage of abundant Alfvénic perturbations—velocity spikes and magnetic switchbacks.

Employing a high-resolution 2.5 MHD model, we scrutinize the genesis of microstreams from emerging bipoles interacting with the ambient corona. Our simulations reveal the tearing-mode instability, forming plasmoids released into the solar wind. Our domain spans the lower corona to 20 Rs, enabling observation of plasmoid formation and their evolution into Alfvénic spikes (Gannouni et al. 2023). The magnetic emergence rate modulates microstream characteristics.

Additionally, 3D MHD simulations explore intermittent interchange reconnection in a 24x24x30Mm domain with a flux emergence rate of 1.38G/s. Reconnection spawns plasma jets and twisted magnetic field bundles, releasing hot plasma into the solar wind, inducing propagating waves and twists along magnetic field lines.

How to cite: Gannouni, B., Réville, V., and Rouillard, A.: Modelling the formation and  evolution of solar wind microstreams:from coronal plumes to propagating Alfvénic velocity spikes., EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-18200, https://doi.org/10.5194/egusphere-egu24-18200, 2024.

Posters virtual: Mon, 15 Apr, 14:00–15:45 | vHall X3

Display time: Mon, 15 Apr 08:30–Mon, 15 Apr 18:00
vX3.1
|
EGU24-6171
TerraVirtuale: a ERC AdG funded project to model planetary magnetospheres with full electron physics. 
(withdrawn)
Giovanni Lapenta
vX3.2
|
EGU24-16790
|
ECS
|
Soumyaranjan Khuntia and Wageesh Mishra

Coronal Mass Ejections (CMEs), huge magnetized plasma erupting from the Sun, pose potential risks to space weather and space-based infrastructure. While extensive research has focused on examining the kinematics of CMEs, there has been limited study of their thermodynamic evolution, particularly at specific heliocentric distances closer to the Sun. Acknowledging that variations in internal plasma properties can impact the overall evolution of CMEs and vice versa is crucial. This study investigates diverse kinematic profiles and associated thermodynamic changes in nine fast CMEs at coronal heights where measuring thermodynamics is challenging. We estimated the distance-dependent evolution of various internal parameters, including polytropic index, temperature, heating rate, pressure, and internal forces driving CME expansion by leveraging the improved Flux Rope Internal State (FRIS) model. The FRIS model utilizes the 3D kinematics derived from the Graduated Cylindrical Shell (GCS) model as input. Our findings reveal that CMEs can maintain their temperature above the adiabatic cooling threshold despite expansion, progressing towards an isothermal state during later propagation phases. The fast CMEs maintaining higher expansion speeds exhibit less pronounced temperature decreases. We found that CME's expansion speed and acceleration correlate well with its maximum temperature drop to reach the isothermal state. Multi-wavelength observations of flux ropes at the source region support the FRIS model-derived results at lower coronal heights. Our analysis elucidates that the primary forces influencing CME radial expansion are the centrifugal and thermal pressure forces, while the Lorentz force acts as a constraining factor. Notably, the thermal pressure force governs expansion at higher heights and is solely responsible for radial expansion. This study enhances our comprehension of the thermodynamic properties of fast CMEs, offering valuable insights for refining assumptions of the polytropic index value in different magnetohydrodynamics (MHD) models to improve the prediction of CME properties at 1 AU.

How to cite: Khuntia, S. and Mishra, W.: Understanding the Thermal Properties of Fast CMEs by Integrating White-light Observations and Analytical Modeling, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-16790, https://doi.org/10.5194/egusphere-egu24-16790, 2024.

vX3.3
|
EGU24-7061
|
ECS
|
Neetasha Arya and Amar Kakad

One of the most important diamagnetic current driven instability transverse to the magnetic field is lower hybrid drift instability (LHDI) which excites lower hybrid drift waves (LHDW). LHDI gives rise to anomalous resistivity which further leads to onset of magnetic reconnection. Because of its high anomalous collision frequency, LHDI enhances the rate of transverse diffusion. The lower hybrid frequency (LHF) ranges between electron and ion cyclotron frequency which is a natural resonance. LHDW is generally observed in transition layer regions and magnetic reconnection sites, where the gradient in density occurs. We are presenting electromagnetic kinetic model including gradients in density and magnetic field, finite parallel wavenumber and non-thermal particle distribution function or kappa velocity distribution function.  The effect of the aforementioned factors on the growth rate of LHDI in different plasma beta circumstances has been thoroughly investigated and will be discussed. Space observation of drift driven wave using MMS spacecraft will also be discussed.

How to cite: Arya, N. and Kakad, A.: Observational and Theoretical Study of Lower Hybrid Drift Waves, EGU General Assembly 2024, Vienna, Austria, 14–19 Apr 2024, EGU24-7061, https://doi.org/10.5194/egusphere-egu24-7061, 2024.