The interaction of the heliosphere with the local interstellar medium (LISM) results in a complicated series of boundary regions. The
Voyager 1
and
2
spacecraft are exploring these distant boundaries ...in situ, as is the Interstellar Boundary Explorer from 1 AU, which measures energetic neutral atoms created in the distant reaches of the heliosphere and LISM. Lyman-α absorption and backscatter measurements also probe the structure and physics of the interface of the heliosphere and LISM. We survey the suite of observations, the underlying theory, and the resulting models that describe the boundary regions of the solar wind and LISM.
The high temperatures and strong magnetic fields of the solar corona form streams of solar wind that expand through the Solar System into interstellar space. At 09:33 UT on 28 April 2021 Parker Solar ...Probe entered the magnetized atmosphere of the Sun 13 million km above the photosphere, crossing below the Alfvén critical surface for five hours into plasma in casual contact with the Sun with an Alfvén Mach number of 0.79 and magnetic pressure dominating both ion and electron pressure. The spectrum of turbulence below the Alfvén critical surface is reported. Magnetic mapping suggests the region was a steady flow emerging on rapidly expanding coronal magnetic field lines lying above a pseudostreamer. The sub-Alfvénic nature of the flow may be due to suppressed magnetic reconnection at the base of the pseudostreamer, as evidenced by unusually low densities in this region and the magnetic mapping.
The theory of nearly incompressible magnetohydrodynamics (NI MHD) was developed largely in the early 1990s, together with an important extension to inhomogeneous flows in 2010. Much of the focus in ...the earlier work was to understand the apparent incompressibility of the solar wind and other plasma environments, and the relationship of density fluctuations to apparently incompressible manifestations of turbulence in the solar wind and interstellar medium. Further important predictions about the "dimensionality" of solar wind turbulence and its relationship to the plasma beta were made and subsequently confirmed observationally. However, despite the initial success of NI MHD in describing fluctuations in the solar wind, a detailed application to solar wind turbulence has not been undertaken. Here, we use the equations of NI MHD to describe solar wind turbulence, rewriting the NI MHD system in terms of Elsässer variables. Distinct descriptions of 2D and slab turbulence emerge naturally from the Elsässer formulation, as do the nonlinear couplings between 2D and slab components. For plasma beta order 1 or less regions, predictions for 2D and slab spectra result from the NI MHD description, and predictions for the spectral characteristics of density fluctuations can be made. We conclude by presenting a NI MHD formulation describing the transport of majority 2D and minority slab turbulence throughout the solar wind. A preliminary comparison of theory and observations is presented.
Numerous problems in space physics and astrophysics require a detailed understanding of the transport and dissipation of low-frequency turbulence in an expanding magnetized flow. We employ a ...scale-separated decomposition of the incompressible MHD equations (based on an Elssasser description) and develop a moment hierarchy to describe the transport of the total energy density in fluctuations, the cross-helicity, the energy difference, and correlation lengths corresponding to forward- and backward-propagating modes and to the energy difference. The dissipation terms for the various transport equations are derived. One-point closure schemes are utilized. The technical elements of this work that distinguish it from previous studies are (1) the inclusion of the large-scale background inhomogeneous Alfvenic velocity VA at a level of detail greater than before, (2) the introduction of a tractable slow timescale closure to eliminate high-frequency interference terms that is likely to prove a useful approximation for practical problems related to the transport of turbulence in an inhomogeneous flow such as the solar wind or solar corona, and finally, (3) we develop a simplified phenomenology for the energy difference or equivalently residual energy that may be useful for practical applications. This yields a coupled system of six equations that describes the transport of turbulence in inhomogeneous sub-Alfvenic and super-Alfvenic flows. The turbulence transport equations are quasi-linear in their spatial evolution operators and nonlinear in the dissipation terms, making the model equations relatively tractable to analysis.
Abstract
The origin, structure, and propagation characteristics of a switchback are compelling questions posed by Parker Solar Probe (PSP) observations of velocity spikes and magnetic field ...reversals. By assuming interchange reconnection between coronal loop and open magnetic field, we show that this results in the generation of upward (into the heliosphere) and downward complex structures propagating at the fast magnetosonic speed (i.e., the Alfvén speed in the low plasma beta corona) that can have an arbitrary radial magnetic field deflection, including “S-shaped.” We derive the evolution equation for the switchback radial magnetic field as it propagates through the inhomogeneous supersonic solar corona. An analytic solution for arbitrary initial conditions is used to investigate the properties of a switchback propagating from launch ∼6 to ∼35
R
⊙
where PSP observed switchbacks during its first encounter. We provide a detailed comparison to an example event, showing that the magnetic field and plasma solutions are in accord with PSP observations. For a simple single switchback, the model predicts either a single or a double-humped structure; the former corresponding to PSP observing either the main body or the flanks of the switchback. The clustering of switchbacks and their sometimes complicated structure may be due to the formation of multiple closely spaced switchbacks created by interchange reconnection with numerous open and loop magnetic field lines over a short period. We show that their evolution yields a complex, aggregated group of switchbacks that includes “sheaths” with large-amplitude radial magnetic field and velocity fluctuations.
Magnetic reconnection is believed to be responsible for the acceleration of energetic electrons with a power-law spectrum in the solar wind and solar flares. However, recent studies of the leading ...mechanism of electron acceleration in reconnection, namely the acceleration by tearing instability induced multi-islands, demonstrates that this mechanism suffers from an "injection problem" for mildly relativistic reconnection acceleration. In this paper, we investigate a new type of reconnection acceleration in which an electron Kelvin-Helmholtz instability (EKHI) is driven as the current sheet reaches electron inertial length scales during magnetic reconnection with a strong guide field. Electrons are accelerated by stochastic electric fields, induced by the EKHI generated vortices that expand rapidly, and a power-law electron energy spectrum with index ∼ 3.5 is produced (W is the electron kinetic energy and f(W) is the energy distribution function). We show that the mechanism is a second-order Fermi acceleration process, and the index where a = Bg/B0, which is determined by the ratio of the spatial scale of the inductive electric field D to that of vortices R and the ratio of guide field Bg to asymptotic magnetic field B0.
The nearly incompressible theory of magnetohydrodynamics (MHD) is formulated in the presence of a static large-scale inhomogeneous background. The theory is an inhomogeneous generalization of the ...homogeneous nearly incompressible MHD description of Zank and Matthaeus and a polytropic equation of state is assumed. The theory is primarily developed to describe solar wind turbulence where the assumption of a composition of two-dimensional (2D) and slab turbulence with the dominance of the 2D component has been used for some time. It was however unclear, if in the presence of a large-scale inhomogeneous background, the dominant component will also be mainly 2D and we consider three distinct MHD regimes for the plasma beta {beta} << 1, {beta} {approx} 1, and{beta} >> 1. For regimes appropriate to the solar wind ({beta} << 1, {beta} {approx} 1), compared to the homogeneous description of Zank and Matthaeus, the reduction of dimensionality for the leading-order description from three dimensional (3D) to 2D is only weak, with the parallel component of the velocity field proportional to the large-scale gradients in density and the magnetic field. Close to the Sun, however, where the large-scale magnetic field can be considered as purely radial, the collapse of dimensionality to 2D is complete. Leading-order density fluctuations are shown to be of the order of the sonic Mach number O(M) and evolve as a passive scalar mixed by the turbulent velocity field. It is emphasized that the usual 'pseudosound' relation used to relate density and pressure fluctuations through the sound speed as {delta}{rho} = c {sup 2} {sub s{delta}}p is not valid for the leading-order O(M) density fluctuations, and therefore in observational studies, the density fluctuations should not be analyzed through the pressure fluctuations. The pseudosound relation is valid only for higher order O(M{sup 2}) density fluctuations, and then only for short-length scales and fast timescales. The spectrum of the leading-order density fluctuations should be modeled as k {sup -5/3} in the inertial range, followed by a Bessel function solution K {sub {nu}}(k), where for stationary turbulence {nu} = 1, in the viscous-convective and diffusion range. Other implications for solar wind turbulence with an emphasis on the evolution of density fluctuations are also discussed.
The Voyager 1 magnetometer has detected several shock waves in the very local interstellar medium (VLISM). Interplanetary shock waves can be transmitted across the heliopause (HP) into the VLISM. The ...first in situ shock observed by Voyager 1 inside the VLISM was remarkably broad and had properties different than those of shocks inside the heliosphere. We present a model of the 2012 VLISM shock, which was observed to be a weak, quasi-perpendicular, low magnetosonic Mach number, low beta, and subcritical shock. Although the heliosphere is a collisionless environment, we show that the VLISM is collisional with respect to the thermal plasma, and that the thermal collisions introduce dissipative terms such as heat conduction and viscosity. The structure of the VLISM shock is determined by thermal proton-proton collisions. VLISM pickup ions (PUIs) do not introduce a significant pressure or dissipation through the shock transition, meaning that the VLISM shock is not mediated by PUIs but only by the thermal gas and magnetic field. Therefore, VLISM shocks are controlled by particle collisions and not by wave-particle interactions. We find that the weak VLISM shock is very broad with a thickness of about 0.12 au, corresponding to the characteristic thermal heat conduction scale length.
Previous application of our focused transport equation for energetic ion test particle acceleration by numerous active small-scale flux ropes to solar wind conditions near 1 au yielded the formation ...of hard power-law spectra with high particle pressure. We present an extended theory where the focused transport equation is coupled to a new MHD turbulence transport equation for coherent, quasi-2D magnetic island structures, based on nearly incompressible (N i) MHD turbulence theory. The latter equation includes new expressions for the magnetic island damping/growth rates that enable a self-consistent description of energy exchange between energetic particles and flux ropes during flux-rope acceleration for four flux-rope acceleration scenarios identified in focused transport theory. Revised, more detailed expressions for coherent acceleration in response to mean dynamic flux-rope properties and for stochastic acceleration due to fluctuations in dynamic flux-rope properties are presented. A comparison is made between the efficiencies of the different flux-rope acceleration scenarios for suprathermal protons in the solar wind near 1 au. Dynamic flux-rope-induced pitch-angle scattering and stochastic acceleration rates are compared with the corresponding rates generated by interaction with parallel-propagating Alfvén waves. The results stress the importance of parallel guiding center motion acceleration by the parallel reconnection electric field formed in merging flux ropes, combined curvature drift and generalized betatron acceleration in contracting/merging flux ropes in the compressible limit (flux-rope compression acceleration), and the fluctuating magnetic mirroring force in flux ropes for pitch-angle scattering.
Abstract
We propose a turbulence-driven solar wind model for a fast solar wind flow in an open coronal hole where the solar wind flow and the magnetic field are highly aligned. We compare the ...numerical results of our model with Parker Solar Probe measurements of the fast solar wind flow and find good agreement between them. We find that (1) the majority quasi-2D turbulence is mainly responsible for coronal heating, raising the temperature to about
K within a few solar radii, which leads in turn to the acceleration of the solar wind; (2) the heating rate due to quasi-2D turbulence near the coronal base is larger than that due to nearly incompressible/slab turbulence; (3) the quasi-2D energy in forward-propagating modes decreases with increasing distance, while the nearly incompressible/slab energy in forward-propagating modes increases, reaching a peak value at ∼11.7
before decreasing with increasing heliocentric distance; (4) the correlation length increases with increasing distance from the coronal base; and (5) the variance of the density fluctuations decreases as a function of heliocentric distance.