We investigate the dynamics of a closed corona cartesian reduced magnetohydrodynamic (MHD) model where photospheric vortices twist the coronal magnetic field lines. We consider two corotating or counter-rotating vortices localized at the center of the photospheric plate, and additionally more corotating vortices that fill the plate entirely. Our investigation is specifically devoted to study the fully nonlinear stage, after the linear stage during which the vortices create laminar and smoothly twisting flux tubes. Our main goal is to understand the dynamics of photospheric vortices twisting the field lines of a coronal magnetic field permeated by finite amplitude broadband fluctuations. We find that depending on the arrangement and handedness of the photospheric vortices an inverse cascade storing a significant amount of magnetic energy may occur or not. In the first case a reservoir of magnetic energy available to large events such as destabilization of a pre-CME configuration develops, while in the second case the outcome is a turbulent heated corona. Although our geometry is simplified our simulations are shown to have relevant implications for coronal dynamics and CME initiation.

The heating of coronal loops is investigated to understand the observational consequences in terms of the thermodynamics and radiative losses from the Sun as well as the magnetized coronae of stars with an outer convective envelope. The dynamics of the Parker coronal heating model are studied for different ratios of the photospheric forcing velocity timescale $t_p$ to the Alfv\'en crossing time along a loop $t_A$. It is shown that for $t_p/t_A \gtrsim$ 10--24 the heating rate and maximum temperature are largest and approximately independent of $t_p/t_A$, leading to a strong emission in X-rays and EUV. On the opposite decreasing $t_p/t_A$ to smaller values leads to lower heating rates and plasma temperatures, and consequently fading high-energy radiative emission once $t_p/t_A \lesssim$ 1--3. The average volumetric loop heating rate is shown to scale as $\ell_p u_p B_0^2/4\pi L^2$, where $\ell_p$ and $u_p$ are respectively the convective granule length-scale and velocity, $B_0$ is the intensity of the strong magnetic field threading the loop, and $L$ the loop length. These findings support a recent hypothesis explaining ultracool dwarf observations of stars with similar magnetic field strength but radically different topologies displaying different radiative emission.

The magnetic topology and field line random walk properties of a nanoflare-heated and magnetically confined corona are investigated in the reduced magnetohydrodynamic regime. Field lines originating from current sheets form coherent structures, called Current Sheet Connected (CSC) regions, extended around them. CSC field line random walk is strongly anisotropic, with preferential diffusion along the current sheets' in-plane length. CSC field line random walk properties remain similar to those of the entire ensemble but exhibit enhanced mean square displacements and separations due to the stronger magnetic field intensities in CSC regions. The implications for particle acceleration and heat transport in the solar corona and wind, and for solar moss formation are discussed.

We investigate the scaling properties of the long-range temporal evolution and intermittency of SDO/AIA intensity observations in four solar environments: an active region core, a weak emission region, and two core loops. We use two approaches: the probability distribution function (PDF) of time series increments, and multifractal detrended fluctuation analysis (MF-DFA). Noise taints the results, so we focus on the 171 Angstrom waveband , which has the highest signal-to-noise ratio. The lags between pairs of wavebands distinguish between coronal versus transition region (TR) emission. In all physical regions studied, scaling in the range 15-45 min is multifractal, and the time series are anti-persistent on the average. The degree of anti-correlation in the TR time series is greater than for coronal emission. The multifractality stems from long term correlations in the data rather than the wide distribution of intensities. Observations in the 335 Angstrom waveband can be described in terms of a multifractal with added noise. The multiscaling of the EUV data agrees qualitatively with the radiance from a phenomenological model of impulsive bursts plus noise, and also from ohmic dissipation in a reduced magnetohydrodynamics (RMHD) model for coronal loop heating. The parameter space must be further explored to seek quantitative agreement. Thus, the observational signatures obtained by the combined tests of the PDF of increments and the MF-DFA offer strong constraints which can systematically discriminate among models for coronal heating.

Magnetic reconnection is thought to be the dynamical mechanism underlying many explosive phenomena observed both in space and in the laboratory, though the question of how fast magnetic reconnection is triggered in such high Lundquist ($S$) number plasmas has remained elusive. It has been well established that reconnection can develop over timescales faster than those predicted traditionally once kinetic scales are reached. It has also been shown that, within the framework of resistive Magnetohydrodynamics (MHD), fast reconnection is achieved for thin enough sheets via the onset of the so-called plasmoid instability. The latter was discovered in studies specifically devoted to the Sweet-Parker current sheet, either as an initial condition or an apparent transient state developing in nonlinear studies. On the other hand, a fast tearing instability can grow on an ideal, i.e., $S$-independent, timescale (dubbed "ideal" tearing) within current sheets whose aspect ratio scales with the macroscopic Lundquist number as $L/a\sim S^{1/3}$ -- much smaller than the Sweet-Parker one -- suggesting a new way to approach to the initiation of fast reconnection in collapsing current configurations. Here we present an overview of what we have called "ideal" tearing in resistive MHD, and discuss how the same reasoning can be extended to other plasma models commonly used that include electron inertia and kinetic effects. We then discuss a scenario for the onset of "ideal" fast reconnection via collapsing current sheets and describe a quantitative model for the interpretation of the nonlinear evolution of "ideally" unstable sheets in two dimensions.

The evolution of a coronal loop is studied by means of numerical simulations of the fully compressible three-dimensional magnetohydrodynamic equations using the HYPERION code. The footpoints of the loop magnetic field are advected by random motions. As a consequence the magnetic field in the loop is energized and develops turbulent nonlinear dynamics characterized by the continuous formation and dissipation of field-aligned current sheets: energy is deposited at small scales where heating occurs. Dissipation is non-uniformly distributed so that only a fraction of the coronal mass and volume gets heated at any time. Temperature and density are highly structured at scales which, in the solar corona, remain observationally unresolved: the plasma of our simulated loop is multi-thermal, where highly dynamical hotter and cooler plasma strands are scattered throughout the loop at sub-observational scales. Numerical simulations of coronal loops of 50000 km length and axial magnetic field intensities ranging from 0.01 to 0.04 Tesla are presented. To connect these simulations to observations we use the computed number densities and temperatures to synthesize the intensities expected in emission lines typically observed with the Extreme ultraviolet Imaging Spectrometer (EIS) on Hinode. These intensities are used to compute differential emission measure distributions using the Monte Carlo Markov Chain code, which are very similar to those derived from observations of solar active regions. We conclude that coronal heating is found to be strongly intermittent in space and time, with only small portions of the coronal loop being heated: in fact, at any given time, most of the corona is cooling down.

The dynamics of magnetic fields in closed regions of solar and stellar coronae are investigated with a reduced magnetohydrodynamic (MHD) model in the framework of Parker scenario for coronal heating. A novel analysis of reduced MHD equilibria shows that their magnetic fields have an asymmetric structure in the axial direction with variation length-scale $z_\ell \sim \ell B_0/b$, where $B_0$ is the intensity of the strong axial guide field, $b$ that of the orthogonal magnetic field component, and $\ell$ the scale of $\mathbf{b}$. Equilibria are then quasi-invariant along the axial direction for variation scales larger than approximatively the loop length $z_\ell \gtrsim L_z$, and increasingly more asymmetric for smaller variation scales $z_\ell \lesssim L_z$. The $critical$ $length$ $z_\ell \sim L_z$ corresponds to the magnetic field intensity threshold $b \sim \ell B_0/L_z$. Magnetic fields stressed by photospheric motions cannot develop strong axial asymmetries. Therefore fields with intensities below such threshold evolve quasi-statically, readjusting to a nearby equilibrium, without developing nonlinear dynamics nor dissipating energy. But stronger fields cannot access their corresponding asymmetric equilibria, hence they are out-of-equilibrium and develop nonlinear dynamics. The subsequent formation of current sheets and energy dissipation is $necessary$ for the magnetic field to relax to equilibrium, since dynamically accessible equilibria have variation scales larger than the loop length $z_\ell \gtrsim L_z$, with intensities smaller than the threshold $b \lesssim \ell B_0/L_z$. The dynamical implications for magnetic fields of interest to solar and stellar coronae are investigated numerically and the impact on coronal physics discussed.

The acceleration of charged particles is relevant to the solar corona over a broad range of scales and energies. High-energy particles are usually detected in concomitance with large energy release events like solar eruptions and flares, nevertheless acceleration can occur at smaller scales, characterized by dynamical activity near current sheets. To gain insight into the complex scenario of coronal charged particle acceleration, we investigate the properties of acceleration with a test-particle approach using three-dimensional magnetohydrodynamic (MHD) models. These are obtained from direct solutions of the reduced MHD equations, well suited for a plasma embedded in a strong axial magnetic field, relevant to the inner heliosphere. A multi-box, multi-scale technique is used to solve the equations of motion for protons. This method allows us to resolve an extended range of scales present in the system, namely from the ion inertial scale of the order of a meter up to macroscopic scales of the order of $10\,$km ($1/100$th of the outer scale of the system). This new technique is useful to identify the mechanisms that, acting at different scales, are responsible for acceleration to high energies of a small fraction of the particles in the coronal plasma. We report results that describe acceleration at different stages over a broad range of time, length and energy scales.

We investigate the dynamical evolution of magnetic fields in closed regions of solar and stellar coronae. To understand under which conditions current sheets form, we examine dissipative and ideal reduced magnetohydrodynamic models in cartesian geometry, where two magnetic field components are present: the strong guide field $B_0$, extended along the axial direction, and the dynamical orthogonal field $\mathbf{b}$. Magnetic field lines thread the system along the axial direction, that spans the length $L$, and are line-tied at the top and bottom plates. The magnetic field $b$ initially has only large scales, with its gradient (current) length-scale of order $\ell_b$. We identify the magnetic intensity threshold $b/B_0 \sim \ell_b/L$. For values of $b$ below this threshold, field-line tension inhibits the formation of current sheets, while above the threshold they form quickly on fast ideal timescales. In the ideal case, above the magnetic threshold, we show that current sheets thickness decreases in time until it becomes smaller than the grid resolution, with the analyticity strip width $\delta$ decreasing at least exponentially, after which the simulations become under-resolved.

We present simulations modeling closed regions of the solar corona threaded by a strong magnetic field where localized photospheric vortical motions twist the coronal field lines. The linear and nonlinear dynamics are investigated in the reduced magnetohydrodynamic regime in Cartesian geometry. Initially the magnetic field lines get twisted and the system becomes unstable to the internal kink mode, confirming and extending previous results. As typical in this kind of investigations, where initial conditions implement smooth fields and flux-tubes, we have neglected fluctuations and the fields are laminar until the instability sets in. But previous investigations indicate that fluctuations, excited by photospheric motions and coronal dynamics, are naturally present at all scales in the coronal fields. Thus, in order to understand the effect of a photospheric vortex on a more realistic corona, we continue the simulations after kink instability sets in, when turbulent fluctuations have already developed in the corona. In the nonlinear stage the system never returns to the simple initial state with ordered twisted field lines, and kink instability does not occur again. Nevertheless field lines get twisted, but in a disordered way, and energy accumulates at large scales through an inverse cascade. This energy can subsequently be released in micro-flares or larger flares, when interaction with neighboring structures occurs or via other mechanisms. The impact on coronal dynamics and CMEs initiation is discussed.

Magnetic reconnection at the interface between coronal holes and loops, so-called interchange reconnection, can release the hotter, denser plasma from magnetically confined regions into the heliosphere, contributing to the formation of the highly variable slow solar wind. The interchange process is often thought to develop at the apex of streamers or pseudo-streamers, near Y and X-type neutral points, but slow streams with loop composition have been recently observed along fanlike open field lines adjacent to closed regions, far from the apex. However, coronal heating models, with magnetic field lines shuffled by convective motions, show that reconnection can occur continuously in unipolar magnetic field regions with no neutral points: photospheric motions induce a magnetohydrodynamic turbulent cascade in the coronal field that creates the necessary small scales, where a sheared magnetic field component orthogonal to the strong axial field is created locally and can reconnect. We propose that a similar mechanism operates near and around boundaries between open and closed regions inducing a continual stochastic rearrangement of connectivity. We examine a reduced magnetohydrodynamic model of a simplified interface region between open and closed corona threaded by a strong unipolar magnetic field. This boundary is not stationary, becomes fractal, and field lines change connectivity continuously, becoming alternatively open and closed. This model suggests that slow wind may originate everywhere along loop-coronal hole boundary regions, and can account naturally and simply for outflows at and adjacent to such boundaries and for the observed diffusion of slow wind around the heliospheric current sheet.

Context. Photospheric motions shuffle the footpoints of the strong axial magnetic field that threads coronal loops giving rise to turbulent nonlinear dynamics characterized by the continuous formation and dissipation of field-aligned current sheets where energy is deposited at small-scales and the heating occurs. Previous studies show that current sheets thickness is orders of magnitude smaller than current state of the art observational resolution (~700 km). Aim. In order to understand coronal heating and interpret correctly observations it is crucial to study the thermodynamics of such a system where energy is deposited at unresolved small-scales. Methods. Fully compressible three-dimensional magnetohydrodynamic simulations are carried out to understand the thermodynamics of coronal heating in the magnetically confined solar corona. Results. We show that temperature is highly structured at scales below observational resolution and nonhomogeneously distributed so that only a fraction of the coronal mass and volume gets heated at each time. Conclusions. This is a multi-thermal system where hotter and cooler plasma strands are found one next to the other also at sub-resolution scales and exhibit a temporal dynamics.

Evidence for inhomogeneous heating in the interplanetary plasma near current sheets dynamically generated by magnetohydrodynamic (MHD) turbulence is obtained using measurements from the ACE spacecraft. These coherent structures only constitute 19% of the data, but contribute 50% of the total plasma internal energy. Intermittent heating manifests as elevations in proton temperature near current sheets, resulting in regional heating and temperature enhancements extending over several hours. The number density of non-Gaussian structures is found to be proportional to the mean proton temperature and solar wind speed. These results suggest magnetofluid turbulence drives intermittent dissipation through a hierarchy of coherent structures, which collectively could be a significant source of coronal and solar wind heating.

The Parker model for coronal heating is investigated through a high resolution simulation. An inertial range is resolved where fluctuating magnetic energy E_M (k_perp) \propto k_\perp^{-2.7} exceeds kinetic energy E_K (k_\perp) \propto k_\perp^{-0.6}. Increments scale as \delta b_\ell \simeq \ell^{-0.85} and \delta u_\ell \simeq \ell^{+0.2} with velocity increasing at small scales, indicating that magnetic reconnection plays a prime role in this turbulent system. We show that spectral energy transport is akin to standard magnetohydrodynamic (MHD) turbulence even for a system of reconnecting current sheets sustained by the boundary. In this new MHD turbulent cascade, kinetic energy flows are negligible while cross-field flows are enhanced, and through a series of "reflections" between the two fields, cascade more than half of the total spectral energy flow.

We present a series of numerical simulations aimed at understanding the nature and origin of turbulence in coronal loops in the framework of the Parker model for coronal heating. A coronal loop is studied via reduced magnetohydrodynamics simulations in Cartesian geometry. A uniform and strong magnetic field threads the volume between the two photospheric planes, where a velocity field in the form of a 1D shear flow pattern is present. Initially the magnetic field which developes in the coronal loop is a simple map of the photospheric velocity field. This initial configuration is unstable to a multiple tearing instability which develops islands with X and O points in the plane orthogonal to the axial field. Once the nonlinear stage sets in the system evolution is characterized by a regime of MHD turbulence dominated by magnetic energy. A well developed power law in energy spectra is observed and the magnetic field never returns to the simple initial state mapping the photospheric flow. The formation of X and O points in the planes orthogonal to the axial field allows the continued and repeated formation and dissipation of small scale current sheets where the plasma is heated. We conclude that the observed turbulent dynamics are not induced by the complexity of the pattern that the magnetic field-lines footpoints follow but they rather stem from the inherent nonlinear nature of the system.

The plasma density enhancements recently observed by the Large-Angle Spectrometric Coronagraph (LASCO) instrument onboard the Solar and Heliospheric Observatory (SOHO) spacecraft have sparked considerable interest. In our previous theoretical study of the formation and initial motion of these density enhancements it is found that beyond the helmet cusp of a coronal streamer the magnetized wake configuration is resistively unstable, that a traveling magnetic island develops at the center of the streamer, and that density enhancements occur within the magnetic islands. As the massive magnetic island travels outward, both its speed and width increase. The island passively traces the acceleration of the inner part of the wake. In the present paper a few spherical geometry effects are included, taking into account either the radial divergence of the magnetic field lines and the average expansion suffered by a parcel of plasma propagating outward, using the Expanding Box Model (EBM), and the diamagnetic force due to the overall magnetic field radial gradients, the so-called melon-seed force. It is found that the values of the acceleration and density contrasts can be in good agreement with LASCO observations, provided the spherical divergence of the magnetic lines starts beyond a critical distance from the Sun and the initial stage of the formation and acceleration of the plasmoid is due to the cartesian evolution of MHD instabilities. This result provides a constraint on the topology of the magnetic field in the coronal streamer.

To understand the nonlinear dynamics of the Parker scenario for coronal heating, long-time high-resolution simulations of the dynamics of a coronal loop in cartesian geometry are carried out. A loop is modeled as a box extended along the direction of the strong magnetic field $B_0$ in which the system is embedded. At the top and bottom plates, which represent the photosphere, velocity fields mimicking photospheric motions are imposed. We show that the nonlinear dynamics is described by different regimes of MHD anisotropic turbulence, with spectra characterized by intertial range power laws whose indexes range from Kolmogorov-like values ($\sim 5/3$) up to $\sim 3$. We briefly describe the bearing for coronal heating rates.

The Parker or field line tangling model of coronal heating is investigated through long-time high-resolution simulations of the dynamics of a coronal loop in cartesian geometry within the framework of reduced magnetohydrodynamics (RMHD). Slow photospheric motions induce a Poynting flux which saturates by driving an anisotropic turbulent cascade dominated by magnetic energy and characterized by current sheets elongated along the axial magnetic field. Increasing the value of the axial magnetic field different regimes of MHD turbulence develop with a bearing on coronal heating rates. In physical space magnetic field lines at the scale of convection cells appear only slightly bended in agreement with observations of large loops of current (E)UV and X-ray imagers.

This work consists of two parts: the first devoted to the study of the heating of the magnetically confined Solar Corona, and the second to the acceleration of the Slow Solar Wind. Direct 3D reduced MHD simulations are presented. They model the heating of coronal loops in the solar atmosphere via the tangling of coronal field lines by photospheric footpoints motions within the framework of the "Parker scenario". We have derived scalings of physical quantities with loop length, and the ratio of photospheric to coronal Alfven velocities. The development of a turbulent dynamics makes the dissipation rate independent of the Reynolds number. The dynamics in physical space are desribed by weak turbulence, which develops when an MHD system is embedded in a strong axial magnetic field. The slow wind originates in and around the coronal streamer belt. The LASCO instrument onboard the SOHO spacecraft has observed plasma density enhancements forming beyond the cusp of a helmet streamer. Previous theoretical models for the formation and initial motion of these density enhancements are improved. The average expansion suffered by a parcel of plasma propagating outward, and the diamagnetic force due to the overall magnetic field radial gradients are now included. It is found that the magnetized wake configuration is resistively unstable, that an outward accelerating magnetic island develops at the center of the streamer, and that density enhancements occur within the magnetic islands. The values of the acceleration and density contrasts can be in good agreement with LASCO observations, provided the spherical divergence of the magnetic lines starts beyond a critical distance from the Sun. This result provides a constraint on the topology of the magnetic field.

MHD turbulence has long been proposed as a mechanism for the heating of coronal loops in the framework of the Parker scenario for coronal heating. So far most of the studies have focused on its dynamical properties without considering its thermodynamical and radiative features, because of the very demanding computational requirements. In this paper we extend this previous research to the compressible regime, including an energy equation, by using HYPERION, a new parallelized, viscoresistive, three-dimensional compressible MHD code. HYPERION employs a Fourier collocation -- finite difference spatial discretization, and uses a third-order Runge-Kutta temporal discretization. We show that the implementation of a thermal conduction parallel to the DC magnetic field induces a radiative emission concentrated at the boundaries, with properties similar to the chromosphere--transition region--corona system.

The Parker or field line tangling model of coronal heating is studied comprehensively via long-time high-resolution simulations of the dynamics of a coronal loop in cartesian geometry within the framework of reduced magnetohydrodynamics (RMHD). Slow photospheric motions induce a Poynting flux which saturates by driving an anisotropic turbulent cascade dominated by magnetic energy. In physical space this corresponds to a magnetic topology where magnetic field lines are barely entangled, nevertheless current sheets (corresponding to the original tangential discontinuities hypothesized by Parker) are continuously formed and dissipated. Current sheets are the result of the nonlinear cascade that transfers energy from the scale of convective motions ($\sim 1,000 km$) down to the dissipative scales, where it is finally converted to heat and/or particle acceleration. Current sheets constitute the dissipative structure of the system, and the associated magnetic reconnection gives rise to impulsive ``bursty'' heating events at the small scales. This picture is consistent with the slender loops observed by state-of-the-art (E)UV and X-ray imagers which, although apparently quiescent, shine bright in these wavelengths with little evidence of entangled features. The different regimes of weak and strong MHD turbulence that develop, and their influence on coronal heating scalings, are shown to depend on the loop parameters, and this dependence is quantitatively characterized: weak turbulence regimes and steeper spectra occur in {\it stronger loop fields} and lead to {\it larger heating rates} than in weak field regions.

Long-time high-resolution simulations of the dynamics of a coronal loop in cartesian geometry are carried out, within the framework of reduced magnetohydrodynamics (RMHD), to understand coronal heating driven by motion of field lines anchored in the photosphere. We unambiguously identify MHD anisotropic turbulence as the physical mechanism responsible for the transport of energy from the large scales, where energy is injected by photospheric motions, to the small scales, where it is dissipated. As the loop parameters vary different regimes of turbulence develop: strong turbulence is found for weak axial magnetic fields and long loops, leading to Kolmogorov-like spectra in the perpendicular direction, while weaker and weaker regimes (steeper spectral slopes of total energy) are found for strong axial magnetic fields and short loops. As a consequence we predict that the scaling of the heating rate with axial magnetic field intensity $B_0$, which depends on the spectral index of total energy for given loop parameters, must vary from $B_0^{3/2}$ for weak fields to $B_0^{2}$ for strong fields at a given aspect ratio. The predicted heating rate is within the lower range of observed active region and quiet Sun coronal energy losses.