We study properties of magnetohydrodynamic (MHD) eigenmodes by decomposing the data of MHD simulations into linear MHD modes - namely the Alfven, slow magnetosonic, and fast magnetosonic modes. We drive turbulence with a mixture of solenoidal and compressive driving, while varying the Alfven Mach number (MA), plasma beta, and the sonic Mach number from sub-sonic to trans-sonic. We find that the proportion of fast and slow modes in the mode mixture increases with increasing compressive forcing. This proportion of the magnetosonic modes can also become the dominant fraction in the mode mixture. The anisotropy of the modes is analyzed by means of their structure functions. The Alfven mode anisotropy is consistent with the Goldreich-Sridhar theory. We find a transition from weak to strong Alfvenic turbulence as we go from low to high MA. The slow mode properties are similar to the Alfven mode. On the other hand the isotropic nature of fast modes is verified in the cases where the fast mode is a significant fraction of the mode mixture. The fast mode behavior does not show any transition in going from low to high MA. We find indications that there is some interaction between the different modes and the properties of the dominant mode can affect the properties of the weaker modes. This work identifies the conditions under which magnetosonic modes can be a major fraction of turbulent astrophysical plasmas, including the regime of weak turbulence. Important astrophysical implications for cosmic ray transport and magnetic reconnection are discussed.

We review and discuss insights on ideal magnetohydrodynamic (MHD) instabilities that can play a role in destabilizing solar coronal flux rope structures. For single flux ropes, failed or actual eruptions may result from internal or external kink evolutions, or from torus unstable configurations. We highlight recent findings from 3D magnetic field reconstructions and simulations where kink and torus instabilities play a prominent role. For interacting current systems, we critically discuss different routes to coronal dynamics and global eruptions, due to current channel coalescence or to tilt-kink scenarios. These scenarios involve the presence of two nearby current channels and are clearly distinct from the popular kink or torus instability. Since the solar corona is pervaded with myriads of magnetic loops -- creating interacting flux ropes typified by parallel or antiparallel current channels as exemplified in various recent observational studies -- coalescence or tilt-kink evolutions must be very common for destabilizing adjacent flux rope systems. Since they also evolve on ideal MHD timescales, they may well drive many sympathetic eruptions witnessed in the solar corona. Moreover, they necessarily lead to thin current sheets that are liable to reconnection. We review findings from 2D and 3D MHD simulations for tilt and coalescence evolutions, as well as on particle acceleration aspects derived from computed charged particle motions in tilt-kink disruptions and coalescing flux ropes. The latter were recently studied in two-way coupled kinetic-fluid models, where the complete phase-space information of reconnection is incorporated.

The equipartition of magnetic and thermal energy in the interstellar medium (ISM) indicates the magneto-hydrodynamic nature of the interstellar turbulence, which can be decomposed into three wave modes: Alfv\'en, fast and slow magnetosonic modes\cite{Herlofson50,LG01,CL03}. Even for studies performed in the case of subsonic turbulence, the magnetosonic modes have been revealed to be more effective than the Alfv\'enic modes for processes such as cosmic ray (CR) transport and acceleration \cite{YL04,Lynn14}. The multiphase nature of ISM and diversity of driving mechanisms give rise to spatial variation of turbulence properties. Nevertheless, the employed model of magneto-hydrodynamic turbulence is often oversimplified being assumed to be only Alfv\'enic or even hydrodynamic due to a lack of observational evidence. Here we report the employment of our novel method, the signature from polarization analysis (SPA), on unveiling the plasma modes in interstellar turbulence. Its application leads to the first discovery of magnetosonic modes in the Cygnus X region. It is found that the magnetosonic modes overlap to a high degree with Fermi cocoon, consistent with theoretical expectations. Moreover, through comparison with the spectrum at other wavelengths, the plasma modes of turbulence are unveiled in active star formation zones. The SPA casts light on the plasma modes composition of the Galactic turbulence, and marks the onset of a new era in the study of interstellar turbulence and accordingly our understandings of relevant processes including cosmic ray transport and star formation.

Simulations of decaying magnetohydrodynamic (MHD) turbulence are performed with a fluid and a kinetic code. The initial condition is an ensemble of long-wavelength, counter-propagating, shear-Alfv\'{e}n waves, which interact and rapidly generate strong MHD turbulence. The total energy is conserved and the rate of turbulent energy decay is very similar in both codes, although the fluid code has numerical dissipation whereas the kinetic code has kinetic dissipation. The inertial range power spectrum index is similar in both the codes. The fluid code shows a perpendicular wavenumber spectral slope of $k_{\perp}^{-1.3}$. The kinetic code shows a spectral slope of $k_{\perp}^{-1.5}$ for smaller simulation domain, and $k_{\perp}^{-1.3}$ for larger domain. We estimate that collisionless damping mechanisms in the kinetic code can account for the dissipation of the observed nonlinear energy cascade. Current sheets are geometrically characterized. Their lengths and widths are in good agreement between the two codes. The length scales linearly with the driving scale of the turbulence. In the fluid code, their thickness is determined by the grid resolution as there is no explicit diffusivity. In the kinetic code, their thickness is very close to the skin-depth, irrespective of the grid resolution. This work shows that kinetic codes can reproduce the MHD inertial range dynamics at large scales, while at the same time capturing important kinetic physics at small scales.