skip to main content

Title: Particle-in-Cell Simulations of Sunward and Anti-sunward Whistler Waves in the Solar Wind

We present particle-in-cell simulations of a combined whistler heat flux and temperature anisotropy instability that is potentially operating in the solar wind. The simulations are performed in a uniform plasma and initialized with core and halo electron populations typical of the solar wind beyond about 0.3 au. We demonstrate that the instability produces whistler-mode waves propagating both along (anti-sunward) and opposite (sunward) to the electron heat flux. The saturated amplitudes of both sunward and anti-sunward whistler waves are strongly correlated with their initial linear growth rates,Bw/B0(γ/ωce)ν, where for typical electron betas we have 0.6 ≲ν≲ 0.9. We show that because of the relatively large spectral width of the whistler waves, the instability saturates through the formation of quasi-linear plateaus around the resonant velocities. The revealed correlations of whistler wave amplitudes and spectral widths with electron beta and temperature anisotropy are consistent with solar wind observations. We show that anti-sunward whistler waves result in an electron heat flux decrease, while sunward whistler waves actually lead to an electron heat flux increase. The net effect is the electron heat flux suppression, whose efficiency is larger for larger electron betas and temperature anisotropies. The electron heat flux suppression can be up to 10%–60% provided that the saturated whistler wave amplitudes exceed about 1% of the background magnetic field. The experimental applications of the presented results are discussed.

more » « less
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Medium: X Size: Article No. 65
["Article No. 65"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We measure the thermal electron energization in 1D and 2D particle-in-cell simulations of quasi-perpendicular, low-beta (βp= 0.25) collisionless ion–electron shocks with mass ratiomi/me= 200, fast Mach numberMms=1–4, and upstream magnetic field angleθBn= 55°–85° from the shock normalnˆ. It is known that shock electron heating is described by an ambipolar,B-parallel electric potential jump, Δϕ, that scales roughly linearly with the electron temperature jump. Our simulations haveΔϕ/(0.5miush2)0.1–0.2 in units of the pre-shock ions’ bulk kinetic energy, in agreement with prior measurements and simulations. Different ways to measureϕ, including the use of de Hoffmann–Teller frame fields, agree to tens-of-percent accuracy. Neglecting off-diagonal electron pressure tensor terms can lead to a systematic underestimate ofϕin our low-βpshocks. We further focus on twoθBn= 65° shocks: aMs=4(MA=1.8) case with a long, 30diprecursor of whistler waves alongnˆ, and aMs=7(MA=3.2) case with a shorter, 5diprecursor of whistlers oblique to bothnˆandB;diis the ion skin depth. Within the precursors,ϕhas a secular rise toward the shock along multiple whistler wavelengths and also has localized spikes within magnetic troughs. In a 1D simulation of theMs=4,θBn= 65° case,ϕshows a weak dependence on the electron plasma-to-cyclotron frequency ratioωpece, andϕdecreases by a factor of 2 asmi/meis raised to the true proton–electron value of 1836.

    more » « less
  2. 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 ofkk, 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δBi2/δBA2=4.1. Incompressible magnetic fluctuations dominate compressible contributions from fast and slow MS modes. The magnetic island spectrum is Kolmogorov-likek1.6in perpendicular wavenumber, and the unidirectional Alfvén wave spectra arek1.6andk1.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
  3. Abstract

    The Parker Solar Probe (PSP) entered a region of sub-Alfvénic solar wind during encounter 8, and we present the first detailed analysis of low-frequency turbulence properties in this novel region. The magnetic field and flow velocity vectors were highly aligned during this interval. By constructing spectrograms of the normalized magnetic helicity, cross-helicity, and residual energy, we find that PSP observed primarily Alfvénic fluctuations, a consequence of the highly field-aligned flow that renders quasi-2D fluctuations unobservable to PSP. We extend Taylor’s hypothesis to sub- and super-Alfvénic flows. Spectra for the fluctuating forward and backward Elsässer variables (z±, respectively) are presented, showing thatz+modes dominatezby an order of magnitude or more, and thez+spectrum is a power law in frequency (parallel wavenumber)f−3/2(k3/2) compared to the convexzspectrum withf−3/2(k3/2) at low frequencies, flattening around a transition frequency (at which the nonlinear and Alfvén timescales are balanced) tof−1.25at higher frequencies. The observed spectra are well fitted using a spectral theory for nearly incompressible magnetohydrodynamics assuming a wavenumber anisotropykk3/4, that thez+fluctuations experience primarily nonlinear interactions, and that the minorityzfluctuations experience both nonlinear and Alfvénic interactions withz+fluctuations. The density spectrum is a power law that resembles neither thez±spectra nor the compressible magnetic field spectrum, suggesting that these are advected entropic rather than magnetosonic modes and not due to the parametric decay instability. Spectra in the neighboring modestly super-Alfvénic intervals are similar.

    more » « less
  4. Abstract

    The 3D radial escape-velocity profile of galaxy clusters has been suggested to be a promising and competitive tool for constraining mass profiles and cosmological parameters in an accelerating universe. However, the observed line-of-sight escape profile is known to be suppressed compared to the underlying 3D radial (or tangential) escape profile. Past work has suggested that velocity anisotropy in the phase-space data is the root cause. Instead, we find that the observed suppression is from the statistical undersampling of the phase spaces and that the 3D radial escape edge can be accurately inferred from projected data. We build an analytical model for this suppression that only requires the number of observed galaxiesNin the phase-space data within the sky-projected range 0.3 ≤r/R200,critical≤ 1. The radially averaged suppression function is an inverse power lawZv=1+(N0/N)λwithN0= 17.818 andλ= 0.362. We test our model withN-body simulations, using dark matter particles, subhalos, and semianalytic galaxies as the phase-space tracers, and find excellent agreement. We also assess the model for systematic biases from cosmology (ΩΛ,H0), cluster mass (M200,critical), and velocity anisotropy (β). We find that varying these parameters over large ranges can impart a maximal additional fractional change in 〈Zv〉 of 2.7%. These systematics are highly subdominant (by at least a factor of 13.7) to the suppression fromN.

    more » « less
  5. Abstract

    We conduct a systematic tidal disruption event (TDE) demographics analysis using the largest sample of optically selected TDEs. A flux-limited, spectroscopically complete sample of 33 TDEs is constructed using the Zwicky Transient Facility over 3 yr (from 2018 October to 2021 September). We infer the black hole (BH) mass (MBH) with host galaxy scaling relations, showing that the sampleMBHranges from 105.1Mto 108.2M. We developed a survey efficiency corrected maximum volume method to infer the rates. The rest-frameg-band luminosity function can be well described by a broken power law ofϕ(Lg)Lg/Lbk0.3+Lg/Lbk2.61, withLbk= 1043.1erg s−1. In the BH mass regime of 105.3≲ (MBH/M) ≲ 107.3, the TDE mass function followsϕ(MBH)MBH0.25, which favors a flat local BH mass function (dnBH/dlogMBHconstant). We confirm the significant rate suppression at the high-mass end (MBH≳ 107.5M), which is consistent with theoretical predictions considering direct capture of hydrogen-burning stars by the event horizon. At a host galaxy mass ofMgal∼ 1010M, the average optical TDE rate is ≈3.2 × 10−5galaxy−1yr−1. We constrain the optical TDE rate to be [3.7, 7.4, and 1.6] × 10−5galaxy−1yr−1in galaxies with red, green, and blue colors.

    more » « less