ANGEO Communicates 20 Sep 2019
ANGEO Communicates  20 Sep 2019
Multichannel coupling of decay instability in threedimensional lowbeta plasma
 ^{1}Institut für Theoretische Physik, Technische Universität Braunschweig, Mendelssohnstr. 3, 38106 Braunschweig, Germany
 ^{2}Institute for Space Sciences, Atomiştilor 409, P.O. Box MG23, BucharestMăgurele, 077125, Romania
 ^{3}Space Research Institute, Austrian Academy of Sciences, Schmiedlstr. 6, 8042 Graz, Austria
 ^{4}Institut für Geophysik und extraterrestrische Physik, Technische Universität Braunschweig, Mendelssohnstr. 3, 38106 Braunschweig, Germany
 ^{5}Deutsches Zentrum für Luft und Raumfahrt, Institut für Planetenforschung, Rutherfordstr. 2, 12489 Berlin, Germany
 ^{1}Institut für Theoretische Physik, Technische Universität Braunschweig, Mendelssohnstr. 3, 38106 Braunschweig, Germany
 ^{2}Institute for Space Sciences, Atomiştilor 409, P.O. Box MG23, BucharestMăgurele, 077125, Romania
 ^{3}Space Research Institute, Austrian Academy of Sciences, Schmiedlstr. 6, 8042 Graz, Austria
 ^{4}Institut für Geophysik und extraterrestrische Physik, Technische Universität Braunschweig, Mendelssohnstr. 3, 38106 Braunschweig, Germany
 ^{5}Deutsches Zentrum für Luft und Raumfahrt, Institut für Planetenforschung, Rutherfordstr. 2, 12489 Berlin, Germany
Correspondence: Horia Comişel (h.comisel@tubraunschweig.de)
Hide author detailsCorrespondence: Horia Comişel (h.comisel@tubraunschweig.de)
Threedimensional hybrid simulations have been carried out to verify the hypothesis of simultaneous multichannel decay of a largeamplitude Alfvén wave in a lowbeta plasma, e.g., in the shockupstream region or the solar corona. Obliquely propagating daughter modes are excited along the perpendicular direction to the mean magnetic field at the same parallel wavenumbers and frequencies as the daughter modes driven by the fieldaligned decay. We find that the transversal spectrum of waves is controlled by the multichannel coupling of the decay process in lowbeta plasmas and originates in the dispersion state of the shear Alfvén wave.
Quasiparallelpropagating Alfvén waves with circular polarization are frequently observed in various space plasma domains, in particular in the vicinity of the Earth's bow shock or other planetary bow shocks (see, e.g., Narita et al., 2007), where the shockreflected ions transfer energy from the ion beam into electromagnetic waves through the ion beam instability (often righthand resonance in which the gyration of the beam ions is observed as righthand polarized in the rest frame of the core or thermal plasma). The beamexcited waves attain large amplitudes up to the saturation level of the beam instability and develop into a stage of wave–wave interactions (see, e.g., Gary, 1991; Akimoto et al., 1993; Gomberoff et al., 2000; Wang and Lin, 2003; Li et al., 2013). Parametric instability is thought to play an important role in distributing the fluctuation energy of largeamplitude Alfvén waves into different modes and different wavelengths (see, e.g., Terasawa et al., 1986; Hoshino and Goldstein, 1989; Spangler et al., 1997; Nariyuki and Hada, 2006; Bekhor and Drake, 2003).
The decay instability is a type of parametric instability and is expected in both low and highbeta plasmas (see, e.g., Inhester, 1990; Vasquez, 1995). We point out that the decay instability in threedimensional nature has multiple channels at once and can excite daughter waves in both the Alfvén mode and sound (or ion acoustic) mode without further cascade of the daughter waves or before the cascade occurs and the wave develops into turbulence. The reason for the multichannel coupling lies in the fact that the decay instability in threedimensional lowbeta plasmas occurs simultaneously in various directions imposed by the mean magnetic field and has a degeneration over the directions perpendicular to the mean field. In the analytical dispersion analysis, one typically sets a priori the propagation direction of the daughter wave and studies the growth rate of the decay instability in the wavenumber domain (see, e.g., Derby, 1978; Goldstein, 1978; Longtin and Sonnerup, 1986; Wong and Goldstein, 1986; Hollweg, 1994; Ruderman and Simpson, 2004; Araneda et al., 2007; Brodin and Stenflo, 2015). The situation is different in a higherdimensional setup; that is, the daughter waves can be excited in various directions at once in the spirit of Boltzmann's equal probability over all the possible channels. The degeneration of the decay instability over the perpendicular directions is illustrated in Fig. 1. Wave–wave coupling in the decay instability forms a parallelogram in the frequency–wavenumber domain (parallel to the mean magnetic field) such that both the energy quantum ℏω and the momentum quantum ℏk are conserved during the wave decay. The pump Alfvén wave (propagating parallel to the mean field, denoted by A_{1}, the firstorder magnetic field after the zeroth order mean field), decays into a backward or antiparallelpropagating Alfvén mode (A_{2}, the secondorder field) and a forwardpropagating sound wave S_{2}. Wave decay in the parallel direction has a degeneracy over the perpendicular directions. The degeneracy can be resolved by plotting parallelogramtype wave couplings in the wave–vector domain spanning the parallel component k_{‖} and the perpendicular component k_{⟂}. That is, the momentum quantum is conserved in the vectorial sense, $\mathrm{\hslash}\mathit{k}=\mathrm{const}.$ One may write the individual realization of the wave decay as “0” for the exactly parallel decay, “+” and “−” for a slightly oblique decay, “$++$” and “−−” for an even more oblique decay, and so on.
Earlier studies on the nonlinear interaction between Alfvén waves and obliquely propagating waves (see, e.g., Mjølhus and Hada, 1990; Viñas and Goldstein, 1991a, b; Laveder et al., 2002; Nariyuki et al., 2008) report that the growth rates of the decay instability in the oblique direction to the mean magnetic field are typically smaller than the case of fieldaligned parametric instabilities. Numerical simulations such as multidimensional MHD (magnetohydrodynamics) and hybrid plasma simulations lead to the conclusion with respect to the circular polarization that the parametric decay of largeamplitude Alfvén waves develops particularly along the propagation direction of the pump wave (see, e.g., Del Zanna et al., 2001; Verscharen et al., 2012; Gao et al., 2013).
The nature of Alfvén wave decay into obliquepropagation angles has been studied extensively. Viñas and Goldstein (1991a, b), for example, argue in their analytical dispersion analysis with circularly polarized Alfvén waves that the obliquepropagation decay instability can compete against the fieldaligned decay in the limit of small propagation or decay angles to the mean magnetic field in the lowbeta plasma. Ghosh et al. (1993) confirm in their MHD twodimensional simulations that obliquely propagating daughter waves are excited by the decay of a fieldaligned Alfvén wave for the lowbeta regime. The results by Ghosh et al. (1993) are based on a β value of 0.5, but nevertheless they demonstrate that the fieldaligned decay instability is persistent and dominates the obliquedecay instabilities. Matteini et al. (2010a) discovered by twodimensional hybrid simulations that linearly polarized Alfvén waves with an oblique direction of propagation decay in lowbeta plasmas in a broad spectrum of coupled Alfvén waves and density fluctuations perpendicular to the direction of the mean magnetic field. Matteini et al. (2010a) also noticed that the magnetic daughter modes follow the dispersion relation for a shear Alfvén wave, $\mathit{\omega}={k}_{\parallel}{V}_{\mathrm{A}}$, where ω, k_{∥}, and V_{A} are the frequency, parallel wavenumber, and Alfvén velocity, respectively.
Here we present a scenario that decay instability can occur simultaneously at various angles to the mean magnetic field, generating a number of secondorder fluctuations or waves (after the pump wave as the firstorder fluctuation). We refer to the simultaneous decay as “multichannel couplings” following the notion in scattering theory, for example, coupled channels in nuclear reactions (see, e.g., Tamura, 1969; Tobocman, 1975). Our goal is to study the hypothesis or the scenario of the multichannel coupling by running a threedimensional hybrid plasma simulation in a lowbeta plasma setup.
2.1 Simulation run
A hybrid plasma simulation is carried out to verify the hypothesis of the multichannel coupling in the lowbeta decay instability. We use the AIKEF hybrid code (Müller et al., 2011) in a threedimensional spatial configuration. The size of the simulation box in each direction is L=288d_{i}, and (576)^{3} computational cells are set, containing 1000 superparticles in each cell. The length scale is normalized to the ion inertia length ${d}_{i}={V}_{\mathrm{A}}/{\mathrm{\Omega}}_{\mathrm{p}}$, where V_{A} and Ω_{p} are the Alfvén velocity and the ion gyrofrequency (for protons), respectively. The values of ion and electron betas used in the simulation are ${\mathit{\beta}}_{i}={\mathit{\beta}}_{e}=\mathrm{0.01}$. This low value of beta (the equivalent fluid beta in the simulation is β=0.02) ensures a faster linear growth and an earlier nonlinear saturation of the parametric instabilities. The magneticfield amplitude of the Alfvén pump wave (normalized to the value of the background magnetic field) has a value of 0.2. The simulation is halted at a time of tΩ_{p}=600, shortly after the saturation of the decay instability.
Figure 2 displays the time evolution of the fluctuation energy for the fieldaligned pump wave (Alfvén mode A_{1}) which represents the Fourier mode $({m}_{\Vert},{m}_{\u27c2})=(\mathrm{10},\mathrm{0})$, (${m}_{\parallel (\u27c2)}={k}_{\parallel (\u27c2)}L/\mathrm{2}\mathit{\pi}$), and the daughters ${A}_{\mathrm{2}}^{\left(\mathrm{0}\right)},{A}_{\mathrm{2}}^{()},{A}_{\mathrm{2}}^{()},{A}_{\mathrm{3}}^{()}$ and ${S}_{\mathrm{2}}^{\left(\mathrm{0}\right)}$, ${S}_{\mathrm{2}}^{()}$, ${S}_{\mathrm{2}}^{()}$, ${S}_{\mathrm{3}}^{()}$, respectively, in the parallel and perpendicular wavenumber domains. Counterpropagating Alfvén daughter modes (to the pump wave propagation direction) and sound daughter modes start to develop simultaneously, and the fluctuation energy increases both along the mean magnetic field and in the oblique directions to the mean field. The mode number for the parallelpropagating Alfvén lower sideband mode is ${m}_{\parallel}=\mathrm{9}$, while the corresponding sound mode has a mode number of ${m}_{\parallel}=\mathrm{19}$, satisfying the threewave coupling rule. An exponential growth of the fluctuation energy represents the linear stage of the decay instability. Both parallel (or fieldaligned) propagating waves and obliquely propagating waves grow simultaneously. The slope of the wave growth is nearly the same between the magnetic mode and the sound mode at lower propagation angles to the mean magnetic field (typically up to 40^{∘}). The decay instability then becomes saturated at a time of tΩ_{p}=300. It is interesting to note that the pump Alfvén wave still dominates the fluctuation energy in the system. The low level of the amplitude of the magnetic field and density daughter modes in Fig. 2 could be a consequence of the lower level of the fluctuation background developed in the threedimensional system with respect to that one developed in a twodimensional system with equivalent characteristics. The obliquely propagating modes are strong enough and compete against the parallelpropagating modes in the linear stage. At later times (tΩ_{p}>300) the fluctuation energy decreases at larger propagation angles to the mean field.
Multichannel coupling in the decay instability is identified in the spectral domain over the perpendicular wavenumbers. The energy spectra for the magnetic field fluctuations, δB^{2}(ω,k), and that of the density, δρ^{2}(ω,k), are plotted as slices in the frequency–wavenumber domain parallel to the mean magnetic field (Fig. 3a) and in the wave–vector domain spanning the perpendicular component k_{⟂} and the parallel component k_{∥} (Fig. 3b). The spectral analysis is performed after the saturation of the instability at the latest time tΩ_{p}≈600 of the simulation. The pump wave is in the Alfvén mode and propagates along the mean magnetic field (${k}_{\u27c2}=\mathrm{0}$). Daughter Alfvén waves and daughter sound waves appear in the wavenumber–frequency spectrum at ${k}_{\parallel}{V}_{\mathrm{A}}/{\mathrm{\Omega}}_{\mathrm{p}}\sim \mathrm{0.2}$ and ${k}_{\parallel}{V}_{\mathrm{A}}/{\mathrm{\Omega}}_{\mathrm{p}}\sim \mathrm{0.4}$, respectively. The magnetic field and density fluctuations in Fig. 3 (bottom panel) are analyzed at the frequency of the fieldaligned daughter modes and differ in the perpendicular wavenumbers (or in the propagation angles). The obliquely propagating Alfvén modes share the same frequency as that of the fieldaligned Alfvén daughter wave. These fluctuations along the perpendicular direction to the mean magnetic field are daughter waves in the Alfvén mode resulting from the multichannel coupling of the decay process. The daughter sound modes with oblique propagation have the same frequency as the parallelpropagating daughter mode ($\mathit{\omega}\sim {k}_{\parallel}{v}_{\mathrm{A}}\sqrt{\mathit{\beta}}\approx \mathrm{0.05}{\mathrm{\Omega}}_{\mathrm{p}}$) and thus fulfill the threewave coupling conservation law. One can mention that in the limit of lowbeta plasmas (β≪1), these oblique compressive waves are still normal modes of plasma.
2.2 Growth rate estimate
Growth rates for the individual decay instability can be computed from the hybrid simulation and are compared to that of the analytical dispersion analysis of the twofluid model developed by Viñas and Goldstein (1991a, b). Briefly, the instability analysis of Viñas and Goldstein (1991a, b) consists in the following. By using the bifluid magnetohydrodynamic model, a linear perturbation of linearmode waves is applied to study the stability of the largeamplitude circularly polarized Alfvén wave. Sideband electromagnetic waves are driven at upper and lower sideband frequencies obeying resonant wave–wave interaction rules. A set of linear equations in the perturbed quantities is obtained, satisfying the resonance condition, while higher sidebands are neglected. We consider the oblique modes at nonzero propagation angles to the mean magnetic field. The determinant of the matrix system is on order 78 in terms of the frequency normalized by using the Alfvén speed V_{A} and the ion gyrofrequency Ω_{p} for protons, $\widehat{\mathit{\omega}}=\mathit{\omega}/\left({k}_{\mathrm{0}}{V}_{\mathrm{A}}\right)=\widehat{{\mathit{\omega}}_{r}}+i\widehat{\mathit{\gamma}}$. The elements of sixbysix order matrices can be found in the Appendix in Viñas and Goldstein (1991a). The dispersion relation ω(k) and the growth rate γ(k) are obtained by solving the determinant of the matrix. We choose a lowbeta value of 0.02 in the semianalytical calculation. The solutions are numerically solved by using the symbolic calculator (Mathematica). We denote that the pump wave has the wavenumber k_{0} and the frequency ω_{0}, the longitudinal daughter wave k and ω, and the Alfvén (sideband) waves k^{±} and ω_{±}. We assume that all the participating waves (pump wave, longitudinal daughter wave, and sideband wave) satisfy the coupling rules, that is,
Figure 4 displays the growth rate $\widehat{\mathit{\gamma}}$ as a function of the wavenumber for the parameters used in the hybrid simulation. Seven values of the propagation angle are chosen for the comparison: $\mathit{\theta}=\mathit{\{}\mathrm{0},\mathrm{10},\mathrm{20},\mathrm{30},\mathrm{45},\mathrm{60},\mathrm{65}{}^{\circ}\mathit{\}}$. The maximum growth rate falls at $k/{k}_{\mathrm{0}}\sim \mathrm{1.9}$ for the fieldaligned decay (θ=0^{∘}). The solutions of the dispersion relation indicate that the growth rates slightly decrease at larger propagation angles. The result of the hybrid simulation is overall consistent with the analytical dispersion analysis. The observed growth rates are overplotted in Fig. 4 and fit reasonably well with the analytical prediction. At larger oblique angles, the prediction has larger deviations with respect to the simulation results. Kinetic effects are most probably responsible for such discrepancies between the MHD model and the hybrid simulation. The fieldaligned decay weakly prevails, but the low and moderate oblique decay instabilities have growth rates close to the former one, confirming the increasing trend with beta decreasing (see, e.g., Viñas and Goldstein, 1991b). Therefore, the multichannel coupling decay instability is regarded as a logical extension of the decay instability theory developed by Viñas and Goldstein (1991a, b) to simultaneous, independent, multiple decay channels.
The damping of the oblique daughter waves observed in Fig. 2 can be explained in terms of the waveparticle processes (e.g., pitchangle scattering, cyclotron resonance, or Landau damping) in competition with the plasma turbulence evolution. We have computed the time evolution of the spectrum of the magnetic field and density fluctuations. In the late stages of the fluctuation evolution, the onedimensional reduced spectra (not shown) along the perpendicular wavenumber reveal turbulence cascades (by the fact that a powerlaw spectrum is being formed) starting from maxima at ${k}_{\u27c2}=\mathrm{0}$ specifically for lowbeta plasma turbulence. The spectral slope along the parallel wavenumber is flatter in the low wavenumber range and becomes steeper at larger values beyond the pump wavenumber. This spectral anisotropy at low and moderate wavenumbers favors the parallelpropagating modes which remain persistent at the latest stage of the decay process.
The oblique electromagnetic daughter waves in Fig. 3 are followed by corresponding fluctuations in the wavenumber–frequency density spectrum (not shown). Viñas and Goldstein (1991a, b) in their analytical study found that the coupling between the electrostatic and electromagnetic solutions increases with the increase in the propagation angle of the daughter wave. On the basis of this remark, we naively interpret this result by interchanging the density daughter modes with the sideband electromagnetic waves. At large oblique propagation angles, the density daughter modes are electromagnetic waves with fluctuation in both density and electromagnetic fields, while the sideband waves become electrostatic waves. Thus, the threewave coupling rule is fulfilled. The obliquely propagating modes develop into a turbulent cascade which saturates the instability at a stage when the Alfvén pump wave still dominates the spectrum. A modulationlike process is observed at early times close to the pump wavenumber and saturates before the saturation of the decay instability. The modulation is exclusively oblique and can be reminiscent of the beat instability. This additional parametric instability and the consequences for the decay instability and its early saturation are still under study.
Matteini et al. (2010a) first noticed the occurrence of a perpendicular spectrum of daughter waves developed at the decay of a linearly polarized Alfvén pump wave with oblique direction of propagation to the background magnetic field. The authors concluded that the generation of obliquely propagating daughter waves is a consequence of the obliqueness of the Alfvén pump wave, i.e., its finite perpendicular wavenumber, in the framework of the threewave coupling process. Our threedimensional treatment is simplified with respect to that former study in the sense that the pump wave is circularly polarized and propagates parallel to the mean magnetic field. Thus we can make useful the predictions of the MHD analytical analysis for the dispersion relations and the growth rates of the decay instability at different propagation angles. The dispersion relation for lefthanded circular polarization, $(k{V}_{\mathrm{A}}/{\mathrm{\Omega}}_{\mathrm{p}}{)}^{\mathrm{2}}=(\mathit{\omega}/{\mathrm{\Omega}}_{\mathrm{p}}{)}^{\mathrm{2}}/(\mathrm{1}\mathit{\omega}/{\mathrm{\Omega}}_{\mathrm{p}})$, is satisfied for both the Alfvén pump wave (as an imposed initial condition) and Alfvén daughter modes. The value of the pump wavenumber, ${k}_{\mathrm{0}}{V}_{\mathrm{A}}/{\mathrm{\Omega}}_{\mathrm{p}}=\mathrm{0.21}$, is fairly high with respect to typical MHD scales, but still it can be assumed to be in the dispersionless range on the first order. The analytical dispersion analysis at weak and moderate oblique propagation angles provides the largest growth rates at approximately the same parallel wavenumber as that one predicted for the fieldaligned decay. This result can be easily understood by introducing the threewave coupling equations and the dispersion relations:
where $({k}_{\mathrm{0}},{\mathit{\omega}}_{\mathrm{0}}={k}_{\mathrm{0}}{V}_{\mathrm{A}})$, $({k}_{\parallel}^{},{k}_{\u27c2}^{},{\mathit{\omega}}^{}={k}_{\parallel}^{}{V}_{\mathrm{A}})$, and (${k}_{\parallel},{k}_{\u27c2},\mathit{\omega}$) are the wavenumbers, frequencies, and dispersion relations characterizing the Alfvén pump, Alfvén daughter, and ion acoustic waves, respectively. The MHD slow mode still preserves the fieldaligned dispersion relation at small and moderate propagation angles in the limit of lowbeta plasmas, namely, ${\mathit{\omega}}_{\mathrm{slow}}={k}_{\parallel}{V}_{\mathrm{A}}\sqrt{\mathit{\beta}}$ (here, $\mathit{\beta}={c}_{s}^{\mathrm{2}}/{V}_{\mathrm{A}}^{\mathrm{2}}$). Therefore, the sound daughter modes with oblique propagation are normal modes of plasma like their counterpart oblique Alfvén modes, and they fulfill similar dispersion relations, i.e., $\mathit{\omega}={k}_{\parallel}{c}_{s}$. Equations (3)–(5), the dispersion relations, and some elementary algebra provide the parallel wavenumber of the oblique ion acoustic wave,
where k_{0} is the wavenumber of the sound daughter driven by the fieldaligned decay (see, e.g., Spangler et al., 1997). The overall decay process is not controlled by the fieldaligned decay but by the dispersion relation of the participating waves which drives the oblique decay to share identical parallel wavenumbers to those attained in the parallel decay. The above rough evaluation is validated by both the analytical analysis in the MHD framework and by the hybrid simulation applied in the present study. In the analytical analysis, the evolution of the daughter waves is determined by the interaction with the fieldaligned pump wave by constructing a system of quasilinear equations which exhibits wave–wave couplings generating obliquely propagating waves as the daughter component. The growth rate of the oblique modes generally has a maximum value at the parallel wavenumber prescribed by the maximum growth rate of the fieldaligned decay instability and depends on the propagation angle.
However, we still point out that the fieldaligned decay remains in any case the fastest one, as shown in Fig. 4. Also, this dynamics is different from the case of an oblique pump wave propagating at a given theta angle with respect to B_{0}, where the oblique mode's decay rate gamma is found to scale cos θ_{kB}, so controlled by the ${k}_{\left\right}$ projection of the initial oblique wave vector (Del Zanna, 2001). This suggests then that there are two possible ways of generating oblique modes from the parametric decay: (1) from a purely parallel mother wave, as in this study; and (2) from an oblique pump wave (see, e.g., Matteini et al., 2010a). In both cases the oblique modes grow at a rate that is smaller than the parallel decay. In our study the central role is played by the conservation of frequencies and the dispersion relations of the components involved in the threewave coupling. Matteini et al. (2010a) show that in their analysis, the central role in driving the transversal spectrum is played by the conservation of the momentum and the nonzero perpendicular projection of the pump wave vector. An obliquely propagating (pump) wave is expected to be more compressive, and it can generate a broadband spectrum of compressive fluctuations. Thus, the configuration with an oblique pump wave should be more efficient in driving oblique modes, even though the daughter waves in our study are driven earlier after a short time of nonlinear evolution. On the other hand, one should not ignore the major differences in the involved setups, e.g., the different polarization of the pump wave and different β_{e} values (see discussion below). Besides the plasma beta, the wavenumber of the pump wave may play a significant role in the process of driving the transverse broadband modulation. In both studies, the medium is still weakly dispersive (${k}_{\mathrm{0}}{V}_{\mathrm{A}}/{\mathrm{\Omega}}_{\mathrm{p}}\sim \mathrm{0.2}$). In a more dispersive medium such as that one used by Verscharen et al. (2012), the dispersion relations probably cannot construct a transversal spectrum of waves for the decay of a fieldaligned Alfvén wave.
Unlike MHD simulations, where the instability can saturate only through the steepening of the excited sound waves, in hybrid simulations it saturates via particle trapping and phasespace modulation (see, e.g., Matteini et al., 2010b). A consequence of this dynamics is a significant perturbation of the ion velocity distribution function, leading also to the generation of fieldaligned beams (see, e.g., Araneda et al., 2008). In lowbeta plasmas, the decay instability dominates over the other parametric instabilities. Plasma beta values for both ions and electrons also play an important role in the dynamics of the parametric decay. For instance, β_{e} is responsible for activating different saturation mechanisms. Matteini et al. (2010b) reported that for cold fluid electrons (β_{e}∼0), the MHD saturation mechanism is recovered. At larger values (e.g., β_{e}∼0.1), the trapping and beam formation is in use. The electron plasma beta β_{e} in our simulation has an intermediate value of 0.01. The parallel distribution in the phase space $z{v}_{\parallel}$ (not shown) at a time still close to the linear phase of the instability (tΩ_{p}=300) suggests that particles are confined and accelerated in different regions along the parallel axis to the mean magnetic field. The phasespace modulation can be a signature for the instability saturation. There is no evidence of a fieldaligned velocity beam at the latest time of the simulation. The accelerated particles are thermalized via pitchangle diffusion by the developed oblique modes.
A consequence of the oblique parametric decay in threedimensional plasmas is a more efficient heating of the ions; see, e.g., the pitch angle scattering study by Comişel et al. (2018). Their conclusion based on a hybrid simulation by using similar parameters to the actual run was that the plateau levels observed in the proton velocity distribution functions are driven by the pitchangle diffusion of the ions by obliquely propagating modes. The spectral analysis performed in Fig. 3 in the frequency–wavevector domain over a broad range of oblique propagation angles proves that the multichannel decay of the parallelpropagating Alfvén pump wave is the source of obliquely propagating daughter waves assumed to satisfy the resonance condition with ions. The former and actual results reconfirm in particular the role of oblique compressive daughter waves in heating the protons parallel and perpendicular to the mean magnetic field according to the in situ observation in solar wind plasmas (see, e.g., Marsch and Tu, 2001).
We have studied by threedimensional hybrid simulations the decay of a largeamplitude Alfvén pump wave expected to achieve the largest growth rate for parallel propagation and circular polarization. Our conclusions are drawn below.

The parallelpropagating Alfvén wave decays into a transversal spectrum of daughter waves. The oblique decay is controlled by the growth rates of the decay instability which are significantly larger at small and moderate propagation angles in lowbeta plasmas, in agreement with the analytic analysis.

The transversal spectrum of daughter waves is controlled by the dispersion relation under the constraint of the threewave coupling process. Thus, the perpendicular alignment to the mean magnetic field originates in the dispersion state of the genuine circularly polarized Alfvén wave.
Kinetic effects are assumed to become more important at larger propagation angles. Their role correlated with the value of the plasma β parameter will be the subject of a forthcoming study on the decay instability in threedimensional plasmas.
Data from our hybrid simulations are stored at the Institut fuer Theoretische Physik – Technische Universitaet Braunschweig. Data can be obtained by writing to the following email addresses: h.comisel@tubraunschweig.de or comisel@spacescience.ro.
HC worked on theory development, simulations, and manuscript writing. YN worked on discussion and manuscript writing. UM worked on theory development, discussion, and coordination of the simulation study.
The authors declare that they have no conflict of interest.
Horia Comişel acknowledges Yasuhiro Nariyuki for fruitful discussions and the hospitality at the University of Toyama. We also acknowledge support by the German Research Foundation and the Open Access Publication Funds of the Technical University of Braunschweig. The work conducted by Horia Comişel in Bucharest is supported by ESA project MAGICS, PRODEX contract C4000127660. The authors gratefully acknowledge the Gauss Centre for Supercomputing e.V. for funding this project by providing computing time through the John von Neumann Institute for Computing (NIC) on the GCS Supercomputer JUWELS at the Jülich Supercomputing Centre (JSC). Furthermore, we are grateful to the anonymous referees for their valuable suggestions and evaluation of this paper.
This research has been supported by the DFG (grant no. MO539/201).
This openaccess publication was funded
by Technische Universität Braunschweig.
This paper was edited by Vincent Maget and reviewed by two anonymous referees.
Akimoto, K., Winske, D., Gary, S. P., and Thomsen, M. F.: Nonlinear evolution of electromagnetic ion beam instabilities, J. Geophys. Res., 98, 1419–1433, https://doi.org/10.1029/92JA02345, 1993. a
Araneda, J. A., Marsch, E., and Viñas, A. F.: Collisionless damping of parametrically unstable Alfvén waves, J. Geophys. Res., 112, A04104, https://doi.org/10.1029/2006JA011999, 2007. a
Araneda, J. A., Marsch, E., and Viñas, A. F.: Proton Core Heating and Beam Formation via Parametrically Unstable AlfvénCyclotron Wave, Phys. Rev. Lett., 100, 12, https://doi.org/10.1103/PhysRevLett.100.125003, 2008. a
Bekhor, S. H. and Drake, R. P.: Plasma heating via parametric beating of Alfvén waves, with heliospheric applications, Phys. Plasmas, 10, 4800, https://doi.org/10.1063/1.1619975, 2003. a
Brodin, G. and Stenflo, L.: Threewave coupling coefficients for perpendicular wave propagation in a magnetized plasma, Phys. Plasmas, 22, 104503, https://doi.org/10.1063/1.4934938, 2015. a
Comişel, H., Nariyuki, Y., Narita, Y., and Motschmann, U.: On heating of solar wind protons by the parametric decay of largeamplitude Alfvén waves, Ann. Geophys., 36, 1647–1655, https://doi.org/10.5194/angeo3616472018, 2018. a
Del Zanna, L.: Parametric decay of oblique arcpolarized Alfvén waves, Geophys. Res. Lett., 28, 2585–2588, https://doi.org/10.1029/2001GL012911, 2001. a
Del Zanna, L., Velli, M., and Londrillo, P.: Parametric decay of circularly polarized Alfvén waves: Multidimensional simulations in periodic and open domains, Astron. Astrophys. 367, 705–718, https://doi.org/10.1051/00046361:20000455, 2001. a
Derby, N. F. J.: Modulational instability of finite amplitude circularly polarized Alfvén waves, Astrophys. J., 224, 1013–1016, 1978. a
Gao, X., Lu, Q., Li, X., Shan, L., and Wang, S.: Parametric instability of a monochromatic Alfven wave: Perpendicular decay in low beta plasma, Phys. Plasmas 20, 072902, https://doi.org/10.1063/1.4816703, 2013. a
Gary, S. P.: Electromagnetic ion/ion instabilities and their consequences in space plasmas: a review, Space Sci. Rev., 56, 373–415, https://doi.org/10.1007/BF00196632, 1991. a
Ghosh, S., Viñas, A. F., and Goldstein, M. L.: Parametric instabilities of a largeamplitude circularly polarized Alfvén wave: Linear growth in twodimensional geometries, J. Geophys. Res., 98, 15561–15570, https://doi.org/10.1029/93JA01534, 1993. a, b
Goldstein, M. L.: An instability of finiteamplitude circularly polarized Alfvén waves, Astrophys. J., 219, 700–704, 1978. a
Gomberoff, K., Gomberoff, L., and, Astudillo, H. F.: Ionbeamplasma electromagnetic instabilities, J. Plasma Phys., 64, 75–87, https://doi.org/10.1017/S0022377800008436, 2000. a
Hollweg, J. V.: The beat, modulational and decay instabilities of a circularly polarized Alfvén wave, J. Geophys. Res., 99, 431, https://doi.org/10.1029/94JA02185, 1994. a
Hoshino, H. and Goldstein, M. L.: Time evolution from linear to nonlinear stages in magnetohydrodynamic parametric instabilities, Phys. Fluids B, 1, 1405, https://doi.org/10.1063/1.858971, 1989. a
Inhester, B.: A driftkinetic treatment of the parametric decay of largeamplitude Alfvén waves, J. Geophys. Res., 95, 10525–10539, 1990. a
Laveder, D., Passot, T., and Sulem, P. L.: Transverse dynamics of dispersive Alfvén waves. II. Driving of a reduced magnetohydrodymanic flow, Phys. Plasmas, 9, 305–314, https://doi.org/10.1063/1.1417511, 2002. a
Li, H., Pang, Y., Huang, S., Zhou, M., Deng, X., Yuan, Z., Wang, D., and Li, H. M.: The turbulence evolution in the high β region of the Earth's foreshock, J. Geophys. Res., 118, 7151–7159, https://doi.org/10.1002/2013JA019424, 2013. a
Longtin, M. and Sonnerup, B. U. Ö.: Modulational instability of circularly polarized Alfvén waves, J. Geophys. Res., 91, 6816–6824, https://doi.org/10.1029/JA091iA06p06816, 1986. a
Marsch, E. and Tu, C.Y.: Evidence for pitch angle diffusion of solar wind protons in resonance with cyclotron waves, J. Geophys. Res., 106, 8357, https://doi.org/10.1029/2000JA000414, 2001. a
Matteini, L., Landi, S., Del Zanna, L., Velli, M., and Hellinger, P.: Parametric decay of linearly polarized shear Alfvén waves in oblique propagation: One and twodimensional hybrid simulations, Geophys. Res. Lett., 37, L20101, https://doi.org/10.1029/2010GL044806, 2010. a, b, c, d, e
Matteini, L., Landi, S., Del Zanna, L., Velli, M., and Hellinger, P.: Kinetics of parametric instabilities of Alfvén waves: Evolution of ion distribution functions, J. Geophys. Res., 115, A09106, https://doi.org/10.1029/2009JA014987, 2010. a, b
Mjølhus, E. and Hada, T.: Oblique stability of circularly polarized MHD waves, J. Plasma Phys., 43, 257–268, https://doi.org/10.1017/S002237780001477X, 1990. a
Müller, J., Simon, S., Motschmann, U., Schüle, J., Glassmeier, K.H., and Pringle, G. J.: A.I.K.E.F.: Adaptive hybrid model for space plasma simulations, Comp. Phys. Comm., 182, 946–966, https://doi.org/10.1016/j.cpc.2010.12.033, 2011. a
Narita, Y., Glassmeier, K.H., Fränz, M., Nariyuki, Y., and Hada, T.: Observations of linear and nonlinear processes in the foreshock wave evolution, Nonlin. Processes Geophys., 14, 361–371, https://doi.org/10.5194/npg143612007, 2007. a
Nariyuki, Y. and Hada, T.: Remarks on nonlinear relation among phases and frequencies in modulational instabilities of parallel propagating Alfvén waves, Nonlin. Processes Geophys., 13, 425–441, https://doi.org/10.5194/npg134252006, 2006. a
Nariyuki, Y., Matsukiyo, S., and, Hada, T.: Parametric instabilities of largeamplitude parallel propagating Alfvén waves: 2D PIC simulation, New J. Phys., 10, 083004, https://doi.org/10.1088/13672630/10/8/083004, 2008. a
Ruderman, M. S. and Simpson, D.: The stability of parallelpropagating circularly polarized Alfvén waves revisited, J. Plasma Phys., 70, 143, https://doi.org/10.1017/S0022377803002599, 2004. a
Spangler, S. R., Leckband, J. A., and Cairns, I. H.: Observations of the parametric decay instability of nonlinear magnetohydrodynamic waves, Phys. Plasmas, 4, 846, https://doi.org/10.1063/1.872183, 1997. a, b
Tamura, T.: Coupledchannel approach to nuclear reaction, Ann. Rev. Nucl. Sci., 19, 99–138, https://doi.org/10.1146/annurev.ns.19.120169.000531, 1969. a
Terasawa, T., Hoshino, M., Sakai, J.I., and Hada, T.: Decay instability of finiteamplitude circularly polarized Alfven Waves: A numerical simulation of stimulated Brillouin scattering, J. Geophys. Res., 91, 4171–4187, https://doi.org/10.1029/JA091iA04p04171, 1986. a
Tobocman, W.: New coupledreactionchannels formalism for nuclear reactions, Phys. Rev. C, 11, 43, https://doi.org/10.1103/PhysRevC.11.43, 1975. a
Vasquez, B. J.: Simulation study of the role of ion kinetics in low frequency wave train evolution, J. Geophys. Res., 100, 1779–1792, https://doi.org/10.1029/94JA02724, 1995. a
Verscharen, D., Marsch, E., Motschmann, U., and Müller, J.: Parametric decay of oblique Alfvén waves in twodimensional hybrid simulations, Phys. Rev. E, 86, 027401, https://doi.org/10.1103/PhysRevE.86.027401, 2012. a, b
Viñas, A. F. and Goldstein, M. L.: Parametric instabilities of circularly polarized largeamplitude dispersive Alfvén waves: excitation of parallelpropagating electromagnetic daughter waves, J. Plasma Phys., 46, 107, https://doi.org/10.1017/S0022377800015981, 1991a. a, b, c, d, e, f, g
Viñas, A. F. and Goldstein, M. L.: Parametric instabilities of circularly polarized largeamplitude dispersive Alfvén waves: excitation of obliquelypropagating daughter and sideband waves, J. Plasma Phys., 46, 129, https://doi.org/10.1017/S0022377800015993, 1991b. a, b, c, d, e, f, g
Wang, X. Y. and Lin, Y.: Generation of nonlinear Alfvén and magnetosonic waves by beamplasma interaction, Phys. Plasmas, 10, 9, https://doi.org/10.1063/1.1599359,2003. a
Wong, H. K. and Goldstein, M. L.: Parametric instabilities of circularly polarized Alfvén waves including dispersion, J. Geophys. Res., 91, 5617–5628, https://doi.org/10.1029/JA091iA05p05617, 1986. a
multichannel couplingsfollowing the notion in scattering theory. Our goal is to study the hypothesis of the multichannel coupling by running a threedimensional hybrid plasma simulation.