We model the time dependent radio emission from a disk accretion event in a T-Tauri star using 3D, ideal magnetohydrodynamic simulations combined with a gyrosynchrotron emission and radiative transfer model. We predict for the first time, the multi-frequency (1$-$1000 GHz) intensity and circular polarisation from a flaring T-Tauri star. A flux tube, connecting the star with its circumstellar disk, is populated with a distribution of non-thermal electrons which is allowed to decay exponentially after a heating event in the disk and the system is allowed to evolve. The energy distribution of the electrons, as well as the non-thermal power law index and loss rate, are varied to see their effect on the overall flux. Spectra are generated from different lines of sight, giving different views of the flux tube and disk. The peak flux typically occurs around 20$-$30 GHz and the radio luminosity is consistent with that observed from T-Tauri stars. For all simulations, the peak flux is found to decrease and move to lower frequencies with elapsing time. The frequency-dependent circular polarisation can reach 10$-$30$\%$ but has a complex structure which evolves as the flare evolves. Our models show that observations of the evolution of the spectrum and its polarisation can provide important constraints on physical properties of the flaring environment and associated accretion event.

The goal of this study is to explore a novel method for the solar photospheric magnetic field diagnostics using Stokes V widths of different magnetosensitive Fe~I spectral lines. We calculate Stokes I and V profiles of several Fe I lines based on a one-dimensional photospheric model VAL C using the NICOLE radiative transfer code. These profiles are used to produce calibration curves linking the intrinsic magnetic field values with the widths of blue peaks of Stokes V profiles. The obtained calibration curves are then tested using the Stokes profiles calculated for more realistic photospheric models based on MHD models of magneto-convection. It is shown that the developed Stokes V widths (SVW) method can be used with various optical and near-infrared lines. Out of six lines considered in this study, FeI 6301 line appears to be the most effective: it is sensitive to fields over ~200G and does not show any saturation up to ~2kG. Other lines considered can also be used for the photospheric field diagnostics with this method, however, only in narrower field value ranges, typically from about 100G to 700-1000G. The developed method can be a useful alternative to the classical magnetic line ratio method, particularly when the choice of lines is limited.

The observed properties (i.e., source size, source position, time duration, decay time) of solar radio emission produced through plasma processes near the local plasma frequency, and hence the interpretation of solar radio bursts, are strongly influenced by propagation effects in the inhomogeneous turbulent solar corona. In this work, a 3D stochastic description of the propagation process is presented, based on the Fokker-Planck and Langevin equations of radio-wave transport in a medium containing anisotropic electron density fluctuations. Using a numerical treatment based on this model, we investigate the characteristic source sizes and burst decay times for Type III solar radio bursts. Comparison of the simulations with the observations of solar radio bursts shows that predominantly perpendicular density fluctuations in the solar corona are required, with an anisotropy factor $\sim 0.3$ for sources observed at around 30~MHz. The simulations also demonstrate that the photons are isotropized near the region of primary emission, but the waves are then focused by large-scale refraction, leading to plasma radio emission directivity that is characterized by a half-width-half-maximum of about 40~degrees near 30~MHz. The results are applicable to various solar radio bursts produced via plasma emission.

Forced magnetic reconnection, a reconnection event triggered by external perturbation, should be ubiquitous in the solar corona. Energy released during such cases can be much greater than that which was introduced by the perturbation. It is unclear how the properties of the external perturbation and the initial current sheet affect the reconnection region properties, and thereby the reconnection dynamics and energy release profile. We investigate the effect of the form of the external perturbation and initial current sheet on the evolution of the reconnection region and the energy release process. Chiefly we explore the non-linear interactions between multiple, simultaneous perturbations, which represent more realistic scenarios. Future work will use these results in test particle simulations to investigate particle acceleration over multiple reconnection events. Simulations are performed using Lare2d, a 2.5D Lagrangian-remap solver for the visco-resistive MHD equations. The model of forced reconnection is extended to include superpositions of sinusoidal driving disturbances, including localised Gaussian perturbations. A transient perturbation is applied to the boundaries of a region containing a force-free current sheet. The simulation domain is sufficiently wide to allow multiple magnetic islands to form and coalesce. Island coalescence contributes significantly to energy release and involves rapid reconnection. Long wavelength modes in perturbations dominate the evolution, without the presence of which reconnection is either slow, as in the case of short wavelength modes, or the initial current sheet remains stable, as in the case of noise perturbations. Multiple perturbations combine in a highly non-linear manner: reconnection is typically faster than when either disturbance is applied individually, with multiple low-energy events contributing to the same total energy release.

Solar flares are explosive events in the solar corona, representing fast conversion of magnetic energy into thermal and kinetic energy, and hence radiation, due to magnetic reconnection. Modelling is essential for understanding and predicting these events. However, self-consistent modelling is extremely difficult due to the vast spatial and temporal scale separation between processes involving thermal plasma (normally considered using magnetohydrodynamic (MHD) approach) and non-thermal plasma (requiring a kinetic approach). In this mini-review we consider different approaches aimed at bridging the gap between fluid and kinetic modelling of solar flares. Two types of approaches are discussed: combined MHD/test-particle (MHDTP) models, which can be used for modelling the flaring corona with relatively small numbers of energetic particles, and hybrid fluid-kinetic methods, which can be used for modelling stronger events with higher numbers of energetic particles. Two specific examples are discussed in more detail: MHDTP models of magnetic reconnection and particle acceleration in kink-unstable twisted coronal loops, and a novel reduced-kinetic model of particle transport.

Magnetic reconnection and particle acceleration due to the kink instability in twisted coronal loops can be a viable scenario for confined solar flares. Detailed investigation of this phenomenon requires reliable methods for observational detection of magnetic twist in solar flares, which may not be possible solely through extreme UV and soft X-ray thermal emission. Polarisation of microwave emission in flaring loops can be used as one of the detection criteria. The aim of this study is to investigate the effect of magnetic twist in flaring coronal loops on the polarisation of gyro-synchrotron microwave (GSMW) emission, and determine whether it could provide a means for magnetic twist detection. We use time-dependent magnetohydrodynamic and test-particle models developed using LARE3D and GCA codes to investigate twisted coronal loops relaxing following the kink-instability. Synthetic GSMW emission maps (I and V Stokes components) are calculated using GX simulator. It is found that flaring twisted coronal loops produce GSMW radiation with a gradient of circular polarisation across the loop. However, these patterns may be visible only for a relatively short period of time due to fast magnetic reconfiguration after the instability. Their visibility also depends on the orientation and position of the loop on solar disk. Typically, it would be difficult to see these characteristic polarisation pattern in a twisted loop seen from the top (close to the centre of the solar disk), but easier in a twisted loop seen from the side (i.e. observed very close to the limb).

Observation of coronal EUV spectral lines sensitive to different temperatures offers an opportunity to evaluate the thermal structure and flows in flaring atmospheres. This can be used to estimate the partitioning between the thermal and kinetic energies. Our aim is to model large-scale (50-10000km) velocity distributions in order to interpret non-thermal broadening of different spectral EUV lines. The developed models allow us to understand the origin of the observed spectral line shifts and broadening, and link these features to particular physical effects in flaring atmospheres. We use ideal MHD to derive twisted magnetic fluxtube configurations in a stratified atmosphere. The evolution of these fluxtubes is followed using resistive MHD, with anomalous resistivity depending on the local density and temperature. The model also takes into account the thermal conduction and radiative losses. The model allows us to evaluate average velocities and velocity dispersions, which would be interpreted as `non-thermal' velocities in observations, at different temperatures for different parts of the models. Our models show qualitative and quantitative agreement with observations. The line-of-sight (LOS) velocity dispersions demonstrate substantial correlation with the temperature, increasing from about 20-30 km/s around 1 MK to about 200-400 km/s near 10-20 MK. The average LOS velocities also correlate with velocity dispersions, although they demonstrate a very strong scattering, compared to observations. We also note that near foot-points the velocity dispersions across the magnetic field are systematically lower than those along the field. We conclude, that the correlation between the flow velocities, velocity dispersions and temperatures are likely to indicate that the same heating mechanism is responsible for heating the plasma, its turbulisation and expansion/evaporation.

The goal of this study is to develop an approximate self-consistent description of particle motion in strongly magnetised solar corona. We derive a set of reduced drift-kinetic equations based on the assumption that the gyro-velocity distribution is Maxwellian. The equations are tested using simple 1D models.

Twisted magnetic fields should be ubiquitous in flare-producing active regions where the magnetic fields are strongly non-potential. It has been shown that reconnection in helical magnetic coronal loops results in plasma heating and particle acceleration distributed within a large volume, including the lower coronal and chromospheric sections of the loops. This scenario can be an alternative to the standard flare model, where particles are accelerated only in a small volume located in the upper corona. We use a combination of MHD simulations and test-particle methods, which describe the development of kink instability and magnetic reconnection in twisted coronal loops using resistive compressible MHD, and incorporate atmospheric stratification and large-scale loop curvature. The resulting distributions of hot plasma let us estimate thermal X-ray emission intensities. The electric and magnetic fields obtained are used to calculate electron trajectories using the guiding-centre approximation. These trajectories combined with the MHD plasma density distributions let us deduce synthetic HXR bremsstrahlung intensities. Our simulations emphasise that the geometry of the emission patterns produced by hot plasma in flaring twisted coronal loops can differ from the actual geometry of the underlying magnetic fields. The twist angles revealed by the emission threads (SXR) are consistently lower than the field-line twist present at the onset of the kink-instability. HXR emission due to the interaction of energetic electrons with the stratified background are concentrated at the loop foot-points in these simulations, even though the electrons are accelerated everywhere within the coronal volume of the loop. The maximum of HXR emission consistently precedes that of SXR emission, with the HXR light-curve being approximately proportional to the temporal derivative of the SXR light-curve.

In the present study we investigate magnetic reconnection in twisted magnetic fluxtubes with different initial configurations. In all considered cases, energy release is triggered by the ideal kink instability, which is itself the result of applying footpoint rotation to an initially potential field. The main goal of this work is to establish the influence of the field topology and various thermodynamic effects on the energy release process. Specifically, we investigate convergence of the magnetic field at the loop footpoints, atmospheric stratification, as well as thermal conduction. In all cases, the application of vortical driving at the footpoints of an initally potential field leads to an internal kink instability. With the exception of the curved loop with high footpoint convergence, the global geometry of the loop change little during the simulation. Footpoint convergence, curvature and atmospheric structure clearly influences the rapidity with which a loop achieves instability as well as the size of the subsequent energy release. Footpoint convergence has a stabilising influence and thus the loop requires more energy for instability, which means that the subsequent relaxation has a larger heating effect. Large-scale curvature has the opposite result: less energy is needed for instability and so the amount of energy released from the field is reduced. Introducing a stratified atmosphere gives rise to decaying wave phenomena during the driving phase, and also results in a loop that is less stable.

Twisted coronal loops should be ubiquitous in the solar corona. Twisted magnetic fields contain excess magnetic energy, which can be released during magnetic reconnection, causing solar flares. The aim of this work is to investigate magnetic reconnection, and particle acceleration and transport in kink-unstable twisted coronal loops, with a focus on the effects of resistivity, loop geometry and atmospheric stratification. Another aim is to perform forward-modelling of bremsstrahlung emission and determine the structure of hard X-ray sources. We use a combination of magnetohydrodynamic (MHD) and test-particle methods. First, the evolution of the kinking coronal loop is considered using resistive MHD model, incorporating atmospheric stratification and loop curvature. Then, the obtained electric and magnetic fields and density distributions are used to calculate electron and proton trajectories using a guiding-centre approximation, taking into account Coulomb collisions. It is shown that electric fields in twisted coronal loops can effectively accelerate protons and electrons to energies up to 10 MeV. High-energy particles have hard, nearly power-law energy spectra. The volume occupied by high-energy particles demonstrates radial expansion, which results in the expansion of the visible hard X-ray loop and a gradual increase in hard X-ray footpoint area. Synthesised hard X-ray emission reveals strong footpoint sources and the extended coronal source, whose intensity strongly depends on the coronal loop density.

We study a model of particle acceleration coupled with an MHD model of magnetic reconnection in unstable twisted coronal loops. The kink instability leads to the formation of helical currents with strong parallel electric fields resulting in electron acceleration. The motion of electrons in the electric and magnetic fields of the reconnecting loop is investigated using a test-particle approach taking into account collisional scattering. We discuss the effects of Coulomb collisions and magnetic convergence near loop footpoints on the spatial distribution and energy spectra of high-energy electron populations and possible implications on the hard X-ray emission in solar flares.

The structure of the photospheric magnetic field during solar flares is examined using echelle spectropolarimetric observations. The study is based on several Fe I and Cr I lines observed at locations corresponding to brightest H$\alpha$ emission during thermal phase of flares. The analysis is performed by comparing magnetic field values deduced from lines with different magnetic sensitivities, as well as by examining the fine structure of $I\pm V$ Stokes profiles splitting. It is shown that the field has at least two components, with stronger unresolved flux tubes embedded in weaker ambient field. Based on a two-component magnetic field model, we compare observed and synthetic line profiles and show that the field strength in small-scale flux tubes is about $2-3$ kG. Furthermore, we find that the small-scale flux tubes are associated with flare emission, which may have implications for flare phenomenology.

In the merging-compression method of plasma start-up, two flux-ropes with parallel toroidal current are formed around in-vessel poloidal field coils, before merging to form a spherical tokamak plasma. This start-up method, used in the Mega-Ampere Spherical Tokamak (MAST), is studied as a high Lundquist number and low plasma-beta magnetic reconnection experiment. In this paper, 2D fluid simulations are presented of this merging process in order to understand the underlying physics, and better interpret the experimental data. These simulations examine the individual and combined effects of tight-aspect ratio geometry and two-fluid physics on the merging. The ideal self-driven flux-rope dynamics are coupled to the diffusion layer physics, resulting in a large range of phenomena. For resistive MHD simulations, the flux-ropes enter the sloshing regime for normalised resistivity eta < 1E-5. In Hall-MHD three regimes are found for the qualitative behaviour of the current sheet, depending on the ratio of the current sheet width to the ion-sound radius. These are a stable collisional regime, an open X-point regime, and an intermediate regime that is highly unstable to tearing-type instabilities. In toroidal axisymmetric geometry, the final state after merging is a MAST-like spherical tokamak with nested flux-surfaces. It is also shown that the evolution of simulated 1D radial density profiles closely resembles the Thomson scattering electron density measurements in MAST. An intuitive explanation for the origin of the measured density structures is proposed, based upon the results of the toroidal Hall-MHD simulations.