Abstract In a collisionless plasma, the energy distribution function of plasma particles can be strongly affected by turbulence. In particular, it can develop a nonthermal power-law tail at high energies. We argue that turbulence with initially relativistically strong magnetic perturbations (magnetization parameterσ≫ 1) quickly evolves into a state with ultrarelativistic plasma temperature but mildly relativistic turbulent fluctuations. We present a phenomenological and numerical study suggesting that in this case, the exponentαin the power-law particle-energy distribution function,f(γ)dγ∝γ−αdγ, depends on magnetic compressibility of turbulence. Our analytic prediction for the scaling exponentαis in good agreement with the numerical results.
more »
« less
Universal Nonthermal Power-law Distribution Functions from the Self-consistent Evolution of Collisionless Electrostatic Plasmas
Abstract Collisionless systems often exhibit nonthermal power-law tails in their distribution functions. Interestingly, collisionless plasmas in various physical scenarios (e.g., the ion population of the solar wind) feature av−5tail in their velocity (v) distribution, whose origin has been a long-standing puzzle. We show this power-law tail to be a natural outcome of the collisionless relaxation of driven electrostatic plasmas. Using a quasi-linear analysis of the perturbed Vlasov–Poisson equations, we show that the coarse-grained mean distribution function (DF),f0, follows a quasi-linear diffusion equation with a diffusion coefficientD(v) that depends onvthrough the plasma dielectric constant. If the plasma is isotropically forced on scales larger than the Debye length with a white-noise-like electric field,D(v) ∼v4forσ<v<ωP/k, withσthe thermal velocity,ωPthe plasma frequency, andkthe characteristic wavenumber of the perturbation; the corresponding quasi-steady-statef0develops av−(d+ 2)tail inddimensions (v−5tail in 3D), while the energy (E) distribution develops anE−2tail independent of dimensionality. Any redness of the noise only alters the scaling in the highvend. Nonresonant particles moving slower than the phase velocity of the plasma waves (ωP/k) experience a Debye-screened electric field, and significantly less (power-law suppressed) acceleration than the near-resonant particles. Thus, a Maxwellian DF develops a power-law tail, while its core (v<σ) eventually also heats up but over a much longer timescale. We definitively show that self-consistency (ignored in test-particle treatments) is crucial for the emergence of the universalv−5tail.
more »
« less
- PAR ID:
- 10558761
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 977
- Issue:
- 1
- ISSN:
- 0004-637X
- Format(s):
- Medium: X Size: Article No. 91
- Size(s):
- Article No. 91
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The quasi-steady states of collisionless plasmas in space (e.g., in the solar wind and planetary environments) are governed by the interactions of charged particles with wave fluctuations. These interactions are responsible not only for the dissipation of plasma waves but also for their excitation. The present analysis focuses on two instabilities, mirror and electromagnetic ion cyclotron instabilities, associated with the same proton temperature anisotropyT⊥>T∥(where ⊥, ∥ are directions defined with respect to the local magnetic field vector). Theories relying on standard Maxwellian models fail to link these two instabilities (i.e., predicted thresholds) to the proton quasi-stable anisotropies measured in situ in a completely satisfactory manner. Here we revisit these instabilities by modeling protons with the generalized bi-Kappa (bi-κpower-law) distribution, and by a comparative analysis of a 2D hybrid simulation with the velocity-moment-based quasi-linear (QL) theory. It is shown that the two methods feature qualitative and, even to some extent, quantitative agreement. The reduced QL analysis based upon the assumption of a time-dependent bi-Kappa model thus becomes a valuable theoretical approach that can be incorporated into the present studies of solar wind dynamics.more » « less
-
Abstract Collisionless low-Mach-number shocks are abundant in astrophysical and space plasma environments, exhibiting complex wave activity and wave–particle interactions. In this paper, we present 2D Particle-in-Cell (PIC) simulations of quasi-perpendicular nonrelativistic (vsh≈ (5500–22000) km s−1) low-Mach-number shocks, with a specific focus on studying electrostatic waves in the shock ramp and precursor regions. In these shocks, an ion-scale oblique whistler wave creates a configuration with two hot counterstreaming electron beams, which drive unstable electron acoustic waves (EAWs) that can turn into electrostatic solitary waves (ESWs) at the late stage of their evolution. By conducting simulations with periodic boundaries, we show that the EAW properties agree with linear dispersion analysis. The characteristics of ESWs in shock simulations, including their wavelength and amplitude, depend on the shock velocity. When extrapolated to shocks with realistic velocities (vsh≈ 300 km s−1), the ESW wavelength is reduced to one-tenth of the electron skin depth and the ESW amplitude is anticipated to surpass that of the quasi-static electric field by more than a factor of 100. These theoretical predictions may explain a discrepancy, between PIC and satellite measurements, in the relative amplitude of high- and low-frequency electric field fluctuations.more » « less
-
Abstract The local escape velocity provides valuable inputs to the mass profile of the galaxy, and requires understanding the tail of the stellar speed distribution. Following Leonard & Tremaine, various works have since modeled the tail of the stellar speed distribution as , wherevescis the escape velocity, andkis the slope of the distribution. In such studies, however, these two parameters were found to be largely degenerate and often a narrow prior is imposed onkin order to constrainvesc. Furthermore, the validity of the power-law form can breakdown in the presence of multiple kinematic substructures or other mis-modeled features in the data. In this paper, we introduce a strategy that for the first time takes into account the presence of kinematic substructure. We model the tail of the velocity distribution as a sum of multiple power laws as a way of introducing a more flexible fitting framework. Using mock data and data from FIRE simulations of Milky Way-like galaxies, we show the robustness of this method in the presence of kinematic structure that is similar to the recently discovered Gaia Sausage. In a companion paper, we present the new measurement of the escape velocity and subsequently the mass of the Milky Way using Gaia eDR3 data.more » « less
-
Let (Mn,g) be a complete simply connectedn-dimensional Riemannian manifold with curvature bounds Sectg≤ κ for κ ≤ 0 and Ricg≥ (n− 1)KgforK≤ 0. We prove that for any bounded domain Ω ⊂Mnwith diameterdand Lipschitz boundary, if Ω* is a geodesic ball in the simply connected space form with constant sectional curvature κ enclosing the same volume as Ω, then σ1(Ω) ≤Cσ1(Ω*), where σ1(Ω) and σ1(Ω*) denote the first nonzero Steklov eigenvalues of Ω and Ω* respectively, andC=C(n, κ,K,d) is an explicit constant. When κ =K, we haveC= 1 and recover the Brock–Weinstock inequality, asserting that geodesic balls uniquely maximize the first nonzero Steklov eigenvalue among domains of the same volume, in Euclidean space and the hyperbolic space.more » « less
An official website of the United States government
