Both remote sensing and in situ measurements show that the fast solar wind plasma significantly deviates from thermal equilibrium and is strongly permeated by turbulent electromagnetic waves, which ...regulate the ion temperature anisotropies and relative drifts. Thus, the ion kinetics is governed by heating and cooling related to absorption and emission of ion‐acoustic and ion‐cyclotron waves, as well as nonresonant pitch angle scattering and diffusion in phase space. Additionally, the solar wind properties are affected by its nonadiabatic expansion as the wind travels away from the Sun. In this study we present results from 1.5‐D hybrid simulations to investigate the effects of a nonlinear turbulent spectrum of Alfvén‐cyclotron waves and the solar wind expansion on the anisotropic heating and differential acceleration of protons and He++ ions. We compare the different heating and acceleration by turbulent Alfvén‐cyclotron wave spectra and by pure monochromatic waves. For the waves and the wave spectra used in our model, we find that the He++ ions are preferentially heated and by the end of the simulations acquire much more than mass‐proportional temperature ratios, Tα/Tp>mα/mp. The differential acceleration between the two species strongly depends on the initial wave amplitude and the related spectral index and is often suppressed by the solar wind expansion. We also find that the expansion leads to perpendicular cooling for both species, and depending on the initial wave spectra, it can either heat or cool the ions in parallel direction. Despite the cooling effect of the expansion in perpendicular direction, the wave‐particle interactions provide an additional heating source, and the perpendicular temperature components remain higher than the adiabatic predictions.
Key Points
He++ ions heating and acceleration by Alfven-cyclotron wave-spectra
Solar wind expansion ‐ perpendicular cooling, parallel heating with wave-spectra
Ion trapping results in prominent proton and He++ ion beam formations
ABSTRACT We perform 2.5D hybrid simulations with massless fluid electrons and kinetic particle-in-cell ions to study the temporal evolution of ion temperatures, temperature anisotropies, and velocity ...distribution functions in relation to the dissipation and turbulent evolution of a broadband spectrum of parallel and obliquely propagating Alfvén-cyclotron waves. The purpose of this paper is to study the relative role of parallel versus oblique Alfvén-cyclotron waves in the observed heating and acceleration of alpha particles in the fast solar wind. We consider collisionless homogeneous multi-species plasma, consisting of isothermal electrons, isotropic protons, and a minor component of drifting particles in a finite-β fast stream near the Earth. The kinetic ions are modeled by initially isotropic Maxwellian velocity distribution functions, which develop nonthermal features and temperature anisotropies when a broadband spectrum of low-frequency nonresonant, ≤ 0.34 p, Alfvén-cyclotron waves is imposed at the beginning of the simulations. The initial plasma parameter values, such as ion density, temperatures, and relative drift speeds, are supplied by fast solar wind observations made by the Wind spacecraft at 1 AU. The imposed broadband wave spectra are left-hand polarized and resemble Wind measurements of Alfvénic turbulence in the solar wind. The imposed magnetic field fluctuations for all cases are within the inertial range of the solar wind turbulence and have a Kraichnan-type spectral slope = −3/2. We vary the propagation angle from θ = 0° to θ = 30° and θ = 60°, and find that the heating of alpha particles is most efficient for the highly oblique waves propagating at 60°, whereas the protons exhibit perpendicular cooling at all propagation angles.
In order to study chromospheric magnetosonic wave propagation including, for the first time, the effects of ion-neutral interactions in the partially ionized solar chromosphere, we have developed a ...new multi-fluid computational model accounting for ionization and recombination reactions in gravitationally stratified magnetized collisional media. The two-fluid model used in our 2D numerical simulations treats neutrals as a separate fluid and considers charged species (electrons and ions) within the resistive MHD approach with Coulomb collisions and anisotropic heat flux determined by Braginskiis transport coefficients. The electromagnetic fields are evolved according to the full Maxwell equations and the solenoidality of the magnetic field is enforced with a hyperbolic divergence-cleaning scheme. The initial density and temperature profiles are similar to VAL III chromospheric model in which dynamical, thermal, and chemical equilibrium are considered to ensure comparison to existing MHD models and avoid artificial numerical heating. In this initial setup we include simple homogeneous flux tube magnetic field configuration and an external photospheric velocity driver to simulate the propagation of MHD waves in the partially ionized reactive chromosphere. In particular, we investigate the loss of chemical equilibrium and the plasma heating related to the steepening of fast magnetosonic wave fronts in the gravitationally stratified medium.
The power spectra of magnetic field fluctuations in the solar wind typically follow a power-law dependence with respect to the observed frequencies and wave-numbers. The background magnetic field ...often influences the plasma properties, setting a preferential direction for plasma heating and acceleration. At the same time the evolution of the solar-wind turbulence at the ion and electron scales is influenced by the plasma properties through local micro-instabilities and wave-particle interactions. The solar-wind-plasma temperature and the solar-wind turbulence at sub- and sup-ion scales simultaneously show anisotropic features, with different components and fluctuation power in parallel with and perpendicular to the orientation of the background magnetic field. The ratio between the power of the magnetic field fluctuations in parallel and perpendicular direction at the ion scales may vary with the heliospheric distance and depends on various parameters, including the local wave properties and nonthermal plasma features, such as temperature anisotropies and relative drift speeds. In this work we have performed two-and-a-half-dimensional hybrid simulations to study the generation and evolution of anisotropic turbulence in a drifting multi-ion species plasma. We investigate the evolution of the turbulent spectral slopes along and across the background magnetic field for the cases of initially isotropic and anisotropic turbulence. Finally, we show the effect of the various turbulent spectra for the local ion heating in the solar wind.
We study the preferential heating and differential acceleration of minor ions by dissipation of ion-acoustic waves (IAWs) generated by parametric instabilities of a finite-amplitude monochromatic ...Alfven-cyclotron pump wave. We consider the associated kinetic effects of Landau damping and nonlinear pitch-angle scattering of protons and a particles in the tenuous plasma of coronal holes and the fast solar wind. Various data collected by Wind spacecraft show signatures for a local transverse heating of the minor ions, presumably by Alfven-cyclotron wave dissipation, and an unexpected parallel heating by a so far unknown mechanism. Here, we present the results from a set of 1.5 dimensional hybrid simulations in search for a plausible explanation for the observed field-aligned kinetic features in the fast solar wind minor ions. We investigate the origin and regulation of ion relative drifts and temperature anisotropies in low plasma beta , fast solar wind conditions. Depending on their initial drifts, both ion species can heat up not only transversely through cyclotron resonance and non-resonant wave-particle interactions, but also strongly in the parallel direction by Landau damping of the daughter IAWs. We discuss the dependence of the relative ion drifts and temperature anisotropies on the plasma beta of the individual species and we describe the effect of the pump wave amplitude on the ion heating and acceleration.
Alfvenic fluctuations in the solar wind exhibit a high degree of velocities and magnetic field correlations consistent with Alfven waves propagating away and toward the Sun. Two remarkable properties ...of these fluctuations are the tendencies to have either positive or negative magnetic helicity (-1 < or =, slant sigma sub(m) < or =, slant +1) associated with either left- or right-topological handedness of the fluctuations and to have a constant magnetic field magnitude. This paper provides, for the first time, a theoretical framework for reconstructing both the magnetic and velocity field fluctuations with a divergence-free magnetic field, with any specified power spectral index and normalized magnetic- and cross-helicity spectrum field fluctuations for any plasma species. The spectrum is constructed in the Fourier domain by imposing two conditions-a divergence-free magnetic field and the preservation of the sense of magnetic helicity in both spaces-as well as using Parseval's theorem for the conservation of energy between configuration and Fourier spaces. Applications to die one-dimensional spatial Alfvenic propagation are presented. The theoretical construction is in agreement with typical time series and power spectra properties observed in the solar wind. The theoretical ideas presented in this spectral reconstruction provide a foundation for more realistic simulations of plasma waves, solar wind turbulence, and the propagation of energetic particles in such fluctuating fields.
Absorption of Alfvén waves is considered to be the main mechanism of heating in the solar corona. It is concluded that the sharp increase of the plasma temperature by two orders of magnitude is ...related to a self-induced opacity with respect to Alfvén waves. The maximal frequency for propagation of Alfvén waves is determined by the strongly temperature dependent kinematic viscosity. In such a way the temperature jump is due to absorption of high frequency Alfvén waves in a narrow layer above the solar surface. It is calculated the power per unit area dissipated in this layer due to damping of Alfvén waves that blows up the plasma and gives birth to the solar wind. A model short wave-length (WKB) evaluation takes into account the 1/f2 frequency dependence of the transversal magnetic field and velocity spectral densities. Such spectral densities agree with old magnetometric data taken by Voyager 1 and recent theoretical calculations in the framework of Langevin-Burgers MHD. The presented theory predicts existence of intensive high frequency MHD Alfvén waves in the cold layer beneath the corona. It is briefly discussed how this statement can be checked experimentally. It is demonstrated that the magnitude of the Alfvén waves generating random noise and the solar wind velocity can be expressed only in terms of satellite experimental data. It is advocated that investigation of properties of the solar surface as a random driver by optical methods is an important task for future solar physics.
Aims.
In the presence of a sufficient amount of target material,
γ
-rays can be used as a tracer in the search for sources of Galactic cosmic rays (CRs). Here we present deep observations of the ...Galactic center (GC) region with the MAGIC telescopes and use them to infer the underlying CR distribution and to study the alleged PeV proton accelerator at the center of our Galaxy.
Methods.
We used data from ≈100 h observations of the GC region conducted with the MAGIC telescopes over five years (from 2012 to 2017). Those were collected at high zenith angles (58−70 deg), leading to a larger energy threshold, but also an increased effective collection area compared to low zenith observations. Using recently developed software tools, we derived the instrument response and background models required for extracting the diffuse emission in the region. We used existing measurements of the gas distribution in the GC region to derive the underlying distribution of CRs. We present a discussion of the associated biases and limitations of such an approach.
Results.
We obtain a significant detection for all four model components used to fit our data (Sgr A*, “Arc”, G0.9+0.1, and an extended component for the Galactic Ridge). We observe no significant difference between the
γ
-ray spectra of the immediate GC surroundings, which we model as a point source (Sgr A*) and the Galactic Ridge. The latter can be described as a power-law with index 2 and an exponential cut-off at around 20 TeV with the significance of the cut-off being only 2
σ
. The derived cosmic-ray profile hints to a peak at the GC position and with a measured profile index of 1.2 ± 0.3 is consistent with the 1/
r
radial distance scaling law, which supports the hypothesis of a CR accelerator at the GC. We argue that the measurements of this profile are presently limited by our knowledge of the gas distribution in the GC vicinity.
Starburst galaxies and star-forming active galactic nuclei are among the candidate sources thought to contribute appreciably to the extragalactic gamma-ray and neutrino backgrounds. NGC 1068 is the ...brightest of the star-forming galaxies found to emit gamma-rays from 0.1 to 50 GeV. Precise measurements of the high-energy spectrum are crucial to study the particle accelerators and probe the dominant emission mechanisms. We have carried out 125 hr of observations of NGC 1068 with the MAGIC telescopes in order to search for gamma-ray emission in the very-high-energy band. We did not detect significant gamma-ray emission, and set upper limits at the 95% confidence level to the gamma-ray flux above 200 GeV f < 5.1 × 10−13 cm−2 s−1. This limit improves previous constraints by about an order of magnitude and allows us to put tight constraints on the theoretical models for the gamma-ray emission. By combining the MAGIC observations with the Fermi-LAT spectrum we limit the parameter space (spectral slope, maximum energy) of the cosmic ray protons predicted by hadronuclear models for the gamma-ray emission, while we find that a model postulating leptonic emission from a semi-relativistic jet is fully consistent with the limits. We provide predictions for IceCube detection of the neutrino signal foreseen in the hadronic scenario. We predict a maximal IceCube neutrino event rate of 0.07 yr−1.
PDF
