skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Search for: All records

Creators/Authors contains: "Zank, G. P."

Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

  1. Context.The heating of the solar corona and solar wind, particularly through suprathermal particles and kinetic Alfvén waves (KAWs) within the 0–10 RSunrange, has been a subject of great interest for many decades. This study investigates and explores the acceleration and heating of charged particles and the role of KAWs in the solar corona. Aims.We investigate how KAWs transport energy and accelerate and heat the charged particles, focusing on the behavior of perturbed electromagnetic (EM) fields, the Poynting flux vectors, net power transfer through the solar flux loop tubes, resonant particles’ speed, group speed, and the damping length of KAWs. The study examines how these elements are influenced by suprathermal particles (κ) and the electron-to-ion temperature ratios (Te/Ti). Methods.We used kinetic plasma theory coupled with the Vlasov-Maxwell model to investigate the dynamics of KAWs and particles. We assumed a collisionless, homogeneous, and low-beta electron-ion plasma in which Alfvén waves travel in the kinetic limits; that is,me/mi ≪ β ≪ 1. Furthermore, the plasma incorporates suprathermal high-energy particles, necessitating an appropriate distribution function to accurately describe the system. We adopted the Kappa distribution function as the most suitable choice for our analysis. Results.The results show that the perturbed EM fields are significantly influenced byκand the effect of Te/Ti. We evaluate both the parallel and perpendicular Poynting fluxes and find that the parallel Poynting flux (Sz) dissipates gradually for lowerκvalues. In contrast, the perpendicular flux (Sx) dissipates quickly over shorter distances. Power deposition in solar flux tubes is significantly influenced byκand Te/Ti. We find that particles can heat the solar corona over long distances (RSun) in the parallel direction and short distances in the perpendicular direction. The group velocity of KAWs increases for lowerκvalues, and the damping length, LG, is enhanced under lowerκ, suggesting longer energy transport distances (RSun). These findings offer a comprehensive understanding of particle-wave interactions in the solar corona and wind, with potential applications for missions such as the Parker Solar Probe, (PSP), and can also apply to other environments where non-Maxwellian particle distributions are frequently observed. 
    more » « less
    Free, publicly-accessible full text available February 1, 2026
  2. Abstract We use in situ measurements from the first 19 encounters of Parker Solar Probe and the most recent five encounters of Solar Orbiter to study the evolution of the turbulent sonic Mach numberMt(the ratio of the amplitude of velocity fluctuations to the sound speed) with radial distance and its relationship to density fluctuations. We focus on the near-Sun region with radial distances ranging from about 11 to 80R. Our results show that (1) the turbulent sonic Mach numberMtgradually moves toward larger values as it approaches the Sun, until at least 11R, whereMtis much larger than the previously observed value of 0.1 at and above 0.3 au; (2) transonic turbulence withMt ∼ 1 is observed in situ for the first time and is found mostly near the Alfvén critical surface; (3) Alfvén Mach number of the bulk flowMAshows a strong correlation with the plasma beta, indicating that most of the observed sub-Alfvénic intervals correspond to a low-beta plasma; (4) the scaling relation between density fluctuations andMtgradually changes from a linear scaling at larger radial distances to a quadratic scaling at smaller radial distances; and (5) transonic turbulence is more compressible than subsonic turbulence, with enhanced density fluctuations and slightly flatter spectra than subsonic turbulence. A systematic understanding of compressible turbulence near the Sun is necessary for future solar wind modeling efforts. 
    more » « less
    Free, publicly-accessible full text available January 13, 2026
  3. Abstract We revisit previous hybrid simulations of the heating and acceleration of interstellar pickup ions (PUIs) at the solar wind termination shock. In previous simulations, a relatively cold initial distribution of PUIs was assumed; and while the resulting shock-heated distribution was consistent with Voyager 2 LECP measurements at about 30 keV, the intensity of the distribution downstream of the shock in the ~1–10 keV energy range was lower than predictions based on analysis of energetic neutral atoms (ENAs) from the Interstellar Boundary Explorer-Hi and Cassini's Ion and Neutral Camera. Here we perform new simulations with more realistic initial PUI distributions. We assume the distribution is a partially filled spherical shell in velocity space with a radius that varies from 320 to 640 km s−1. We then use the distributions downstream of the shock from these new simulations to estimate the ENA flux spectrum and compare with observations. We find that the predicted ENA spectrum from the new simulations much better matches the observations over a broad range of energies. We conclude that the hybrid simulations provide reasonable predictions for the distribution of charged particles in the energy range from ~0.5 to 50 keV. 
    more » « less
    Free, publicly-accessible full text available February 3, 2026
  4. Abstract Magnetic field fluctuations measured in the heliosheath by the Voyager spacecraft are often characterized as compressible, as indicated by a strong fluctuating component parallel to the mean magnetic field. However, the interpretation of the turbulence data faces the caveat that the standard Taylor’s hypothesis is invalid because the solar wind flow velocity in the heliosheath becomes subsonic and slower than the fast magnetosonic speed, given the contributions from hot pickup ions (PUIs) in the heliosheath. We attempt to overcome this caveat by introducing a 4D frequency-wavenumber spectral modeling of turbulence, which is essentially a decomposition of different wave modes following their respective dispersion relations. Isotropic Alfvén and fast mode turbulence are considered to represent the heliosheath fluctuations. We also include two dispersive fast wave modes derived from a three-fluid theory. We find that (1) magnetic fluctuations in the inner heliosheath are less compressible than previously thought, an isotropic turbulence spectral model with about 25% in compressible fluctuation power is consistent with the observed magnetic compressibility in the heliosheath; (2) the hot PUI component and the relatively cold solar wind ions induce two dispersive fast magnetosonic wave branches in the perpendicular propagation limit, PUI fast wave may account for the spectral bump near the proton gyrofrequency in the observable spectrum; (3) it is possible that the turbulence wavenumber spectrum is not Kolmogorov-like although the observed frequency spectrum has a −5/3 power-law index, depending on the partitioning of power among the various wave modes, and this partitioning may change with wavenumber. 
    more » « less
    Free, publicly-accessible full text available September 1, 2025
  5. On electron kinetic scales, ions and electrons decouple, and electron velocity shear on electron inertial length ∼de can trigger electromagnetic (EM) electron Kelvin–Helmholtz instability (EKHI). In this paper, we present an analytic study of EM EKHI in an inviscid collisionless plasma with a step-function electron shear flow. We show that in incompressible collisionless plasma, the ideal electron frozen-in condition E+ve×B/c=0 must be broken for the EM EKHI to occur. In a step-function electron shear flow, the ideal electron frozen-in condition is replaced by magnetic flux conservation, i.e., ∇×(E+ve×B/c)=0, resulting in a dispersion relation similar to that of the standard ideal and incompressible magnetohydrodynamics KHI. The magnetic field parallel to the electron streaming suppresses the EM EKHI due to magnetic tension. The threshold for the EM mode of the EKHI is (k·ΔUe)2>ne1+ne2ne1ne2[ne1(vAe1·k)2+ne2(vAe2·k)2], where vAe=B/(4πmene)1/2, ΔUe, and ne are the electron streaming velocity shear and densities, respectively. The growth rate of the EM mode is γem∼Ωce, which is the electron gyro-frequency. 
    more » « less
  6. Abstract A steady-state, semi-analytical model of energetic particle acceleration in radio-jet shear flows due to cosmic-ray viscosity obtained by Webb et al. is generalized to take into account more general cosmic-ray boundary spectra. This involves solving a mixed Dirichlet–Von Neumann boundary value problem at the edge of the jet. The energetic particle distribution functionf0(r,p) at cylindrical radiusrfrom the jet axis (assumed to lie along thez-axis) is given by convolving the particle momentum spectrum f 0 ( , p ) with the Green’s function G ( r , p ; p ) , which describes the monoenergetic spectrum solution in which f 0 δ ( p p ) asr→ ∞ . Previous work by Webb et al. studied only the Green’s function solution for G ( r , p ; p ) . In this paper, we explore for the first time, solutions for more general and realistic forms for f 0 ( , p ) . The flow velocityu=u(r)ezis along the axis of the jet (thez-axis).uis independent ofz, andu(r) is a monotonic decreasing function ofr. The scattering time τ ( r , p ) = τ 0 ( p / p 0 ) α in the shear flow region 0 <r<r2, and τ ( r , p ) = τ 0 ( p / p 0 ) α ( r / r 2 ) s , wheres> 0 in the regionr>r2is outside the jet. Other original aspects of the analysis are (i) the use of cosmic ray flow lines in (r,p) space to clarify the particle spatial transport and momentum changes and (ii) the determination of the probability distribution ψ p ( r , p ; p ) that particles observed at (r,p) originated fromr→ ∞ with momentum p . The acceleration of ultrahigh-energy cosmic rays in active galactic nuclei jet sources is discussed. Leaky box models for electron acceleration are described. 
    more » « less
  7. Abstract In January 2021, Metis/SolO and PSP formed a quadrature from which the slow solar wind was able to be measured from the extended solar corona (3.5 – 6.3 R ⊙ ) to the very inner heliosphere (23.2 R ⊙ ). Metis/SolO remotely measured the coronal solar wind, finding a speed of 96 – 201 kms −1 , and PSP measured the solar wind in situ, finding a speed of 219.34 kms −1 . Similarly, the normalized cross-helicity and the normalized residual energy measured by PSP are 0.96 and -0.07. In this manuscript, we study the evolution of the proton entropy and the turbulence cascade rate of the outward Elsässer energy during this quadrature. We also study the relationship between solar wind speed, density and temperature, and their relationship with the turbulence energy, the turbulence cascade rate, and the solar wind proton entropy. We compare the theoretical results with the observed results measured by Metis/SolO and PSP. 
    more » « less
  8. Abstract Parker Solar Probe (PSP) observed sub-Alfvénic solar wind intervals during encounters 8–14, and low-frequency magnetohydrodynamic (MHD) turbulence in these regions may differ from that in super-Alfvénic wind. We apply a new mode decomposition analysis to the sub-Alfvénic flow observed by PSP on 2021 April 28, identifying and characterizing entropy, magnetic islands, forward and backward Alfvén waves, including weakly/nonpropagating Alfvén vortices, forward and backward fast and slow magnetosonic (MS) modes. Density fluctuations are primarily and almost equally entropy- and backward-propagating slow MS modes. The mode decomposition provides phase information (frequency and wavenumberk) for each mode. Entropy density fluctuations have a wavenumber anisotropy ofk≫k, whereas slow-mode density fluctuations havek>k. Magnetic field fluctuations are primarily magnetic island modes (δBi) with anO(1) smaller contribution from unidirectionally propagating Alfvén waves (δBA+) giving a variance anisotropy of δ B i 2 / δ B A 2 = 4.1 . Incompressible magnetic fluctuations dominate compressible contributions from fast and slow MS modes. The magnetic island spectrum is Kolmogorov-like k 1.6 in perpendicular wavenumber, and the unidirectional Alfvén wave spectra are k 1.6 and k 1.5 . Fast MS modes propagate at essentially the Alfvén speed with anticorrelated transverse velocity and magnetic field fluctuations and are almost exclusively magnetic due toβp≪ 1. Transverse velocity fluctuations are the dominant velocity component in fast MS modes, and longitudinal fluctuations dominate in slow modes. Mode decomposition is an effective tool in identifying the basic building blocks of MHD turbulence and provides detailed phase information about each of the modes. 
    more » « less
  9. The role of pickup ions (PUIs) in the solar wind interaction with the local interstellar medium is investigated with 3D, multifluid simulations. The flow of the mixture of all charged particles is described by the ideal MHD equations, with the source terms responsible for charge exchange between ions and neutral atoms. The thermodynamically distinct populations of neutrals are governed by individual sets of gas dynamics Euler equations. PUIs are treated as a separate, comoving fluid. Because the anisotropic behavior of PUIs at the heliospheric termination shocks is not described by the standard conservation laws (a.k.a. the Rankine–Hugoniot relations), we derived boundary conditions for them, which are obtained from the dedicated kinetic simulations of collisionless shocks. It is demonstrated that this approach to treating PUIs makes the computation results more consistent with observational data. In particular, the PUI pressure in the inner heliosheath (IHS) becomes higher by ∼40%–50% in the new model, as compared with the solutions where no special boundary conditions are applied. Hotter PUIs eventually lead to charge-exchange-driven cooling of the IHS plasma, which reduces the IHS width by ∼15% (∼8–10 au) in the upwind direction, and even more in the other directions. The density of secondary neutral atoms born in the IHS decreases by ∼30%, while their temperature increases by ∼60%. Simulation results are validated with New Horizons data at distances between 11 and 47 au. 
    more » « less
  10. Abstract We compare hybrid (kinetic proton, fluid electron) and particle-in-cell (kinetic proton, kinetic electron) simulations of the solar wind termination shock with parameters similar to those observed by Voyager 2 during its crossing. The steady-state results show excellent agreement between the downstream variations in the density, plasma velocity, and magnetic field. The quasi-perpendicular shock accelerates interstellar pickup ions to a maximum energy limited by the size of the computational domain, with somewhat higher fluxes and maximal energies observed in the particle-in-cell simulation, likely due to differences in the cross-shock electric field arising from electron kinetic-scale effects. The higher fluxes may help address recent discrepancies noted between observations and large-scale hybrid simulations. 
    more » « less