Because of the development of many-body theories of nuclear matter, the long-standing, open problem of the equation of state (EOS) of dense matter may be understood in the near future through the confrontation of theoretical calculations with laboratory measurements of nuclear properties \& reactions and increasingly accurate observations in astronomy. In this review, we focus on the following six aspects: 1) providing a survey of the quark mean-field (QMF) model, which consistently describes a nucleon and many-body nucleonic system from a quark potential; 2) applying QMF to both nuclear matter and neutron stars; 3) extending QMF formalism to the description of hypernuclei and hyperon matter, as well as hyperon stars; 4) exploring the hadron-quark phase transition and hybrid stars by combining the QMF model with the quark matter model characterized by the sound speed; 5) constraining interquark interactions through both the gravitational wave signals and electromagnetic signals of binary merger event GW170817; and 6) discussing further opportunities to study dense matter EOS from compact objects, such as neutron star cooling and pulsar glitches.

We investigate the stability and $e^+e^-$ pair creation of supercritically charged superheavy nuclei, $ud$QM nuggets, strangelets, and strangeon nuggets based on Thomas-Fermi approximation. The model parameters are fixed by reproducing their masses and charge properties reported in earlier publications. It is found that $ud$QM nuggets, strangelets, and strangeon nuggets may be more stable than ${}^{56}$Fe at $A\gtrsim 315$, $5\times10^4$, and $1.2\times10^8$, respectively. For those stable against neutron emission, the most massive superheavy element has a baryon number $\sim$965, while $ud$QM nuggets, strangelets, and strangeon nuggets need to have baryon numbers larger than $39$, 433, and $2.7\times10^5$. The $e^+e^-$ pair creation will inevitably start for superheavy nuclei with charge numbers $Z\geq177$, $ud$QM nuggets with $Z\geq163$, strangelets with $Z\geq 192$, and strangeon nuggets with $Z\geq 212$. A universal relation $Q/R_e = \left(m_e - \bar{\mu}_e\right)/\alpha$ is obtained at a given electron chemical potential $\bar{\mu}_e$, where $Q$ is the total charge and $R_e$ the radius of electron cloud. This predicts the maximum charge number by taking $\bar{\mu}_e=-m_e$. For supercritically charged objects with $\bar{\mu}_e<-m_e$, the decay rate for $e^+e^-$ pair production is estimated based on the JWKB approximation. It is found that most positrons are emitted at $t\lesssim 10^{-15}$ s, while a long lasting positron emission is observed for large objects with $R\gtrsim 1000$ fm. The emission and annihilation of positrons from supercritically charged objects may be partially responsible for the short $\gamma$-ray burst during the merger of binary compact stars, the 511 keV continuum emission, as well as the narrow faint emission lines in X-ray spectra from galaxies and galaxy clusters.

The 'standard solar' flare model collects all physical ingredients identified by multi-wavelength observations of our Sun: magnetic reconnection, fast particle acceleration and the resulting emission at various wavelengths, especially in soft to hard X-ray channels. Its cartoon representation is found throughout textbooks on solar and plasma astrophysics, and guides interpretations of unresolved energetic flaring events on other stars, accretion disks and jets. To date, a fully self-consistent model that reproduces the standard scenario in all its facets is lacking, since this requires the combination of a large scale, multi-dimensional magnetohydrodynamic (MHD) plasma description with a realistic fast electron treatment. Here, we demonstrate such a novel combination, where MHD combines with an analytic fast electron model, adjusted to handle time-evolving, reconnecting magnetic fields and particle trapping. This allows to study (1) the role of fast electron deposition in the triggering of chromospheric evaporation flows; (2) the physical mechanisms that generate various hard X-ray sources at chromospheric footpoints or looptops; and (3) the relationship between soft X-ray and hard X-ray fluxes throughout the entire flare loop evolution. For the first time, this self-consistent solar flare model demonstrates the observationally suggested relationship between flux swept out by the hard X-ray footpoint regions, and the actual reconnection rate at the X-point, which is a major unknown in flaring scenarios. We also demonstrate that a looptop hard X-ray source can result from fast electron trapping.

We report on the latest additions to our open-source, block-grid adaptive framework MPI-AMRVAC, which is a general toolkit for especially hyperbolic/parabolic partial differential equations (PDEs). Applications traditionally focused on shock-dominated, magnetized plasma dynamics described by either Newtonian or special relativistic (magneto)hydrodynamics, but its versatile design easily extends to different PDE systems. Here, we demonstrate applications covering any-dimensional scalar to system PDEs, with e.g. Korteweg-de Vries solutions generalizing early findings on soliton behaviour, shallow water applications in round or square pools, hydrodynamic convergence tests as well as challenging computational fluid and plasma dynamics applications. The recent addition of a parallel multigrid solver opens up new avenues where also elliptic constraints or stiff source terms play a central role. This is illustrated here by solving several multi-dimensional reaction-diffusion-type equations. We document the minimal requirements for adding a new physics module governed by any nonlinear PDE system, such that it can directly benefit from the code flexibility in combining various temporal and spatial discretisation schemes. Distributed through GitHub, MPI-AMRVAC can be used to perform 1D, 1.5D, 2D, 2.5D or 3D simulations in Cartesian, cylindrical or spherical coordinate systems, using parallel domain-decomposition, or exploiting fully dynamic block quadtree-octree grids.

Prominences show a surprising amount of fine structure and it is widely believed that their threads, as seen in H$\alpha$ observations, provide indirect information concerning magnetic field topology. We investigate the spontaneous emergence and evolution of fine structure in high-density condensations formed through the process of thermal instability under typical solar coronal conditions. Our study reveals intricate multidimensional processes that occur through in situ condensations in a representative coronal volume in a low plasma beta regime. We performed 2D and 3D numerical simulations of interacting slow magnetohydrodynamic (MHD) wave modes when all relevant non-adiabatic effects are included. We show that the interaction of multiple slow MHD wave modes in a regime unstable to the thermal mode leads to thermal instability. This initially forms pancake-like structures almost orthogonal to the local magnetic field, while low-pressure induced inflows of matter generate rebound shocks. This is succeeded by the rapid disruption of these pancake sheets through thin-shell instabilities evolving naturally from minute ram pressure imbalances. This eventually creates high-density blobs accompanied by thread-like features from shear flow effects. The further evolution of the blobs follows the magnetic field lines, such that a dynamical realignment with the background magnetic field appears. However, the emerging thread-like features are not at all field-aligned, implying only a very weak link between fine structure orientation and magnetic field topology which has far-reaching implications for field topology interpretations based on H$\alpha$ observations.

The symmetry energy obtained with the effective Skyrme energy density functional is related to the values of isoscalar effective mass and isovector effective mass, which is also indirectly related to the incompressibility of symmetric nuclear matter. In this work, we analyze the values of symmetry energy and its related nuclear matter parameters in five-dimensional parameter space by describing the heavy ion collision data, such as isospin diffusion data at 35 MeV/u and 50 MeV/u, neutron skin of $^{208}$Pb, and tidal deformability and maximum mass of neutron star. We obtain the parameter sets which can describe the isospin diffusion, neutron skin, tidal deformability and maximum mass of neutron star, and give the incompressibility $K_0$=250.23$\pm$20.16 MeV, symmetry energy coefficient $S_0$=31.35$\pm$2.08 MeV, the slope of symmetry energy $L$=59.57$\pm$10.06 MeV, isoscalar effective mass $m_s^*/m$=0.75$\pm$0.05 and quantity related to effective mass splitting $f_I$=0.005$\pm$0.170. At two times normal density, the symmetry energy we obtained is in 35-55 MeV. To reduce the large uncertainties of $f_I$, more critical works in heavy ion collisions at different beam energies are needed.

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 phase state in the intermediate density range of dense matter ($\sim1-10$ times of nuclear saturation density) is both intriguing and unclear, and could have important observable effects in the present gravitational wave era of neutron stars. Since the matter density in neutron stars is in the nonperturbative interaction region, the sound velocity is expected to approach the conformal limit ($c_s/c=1/\sqrt{3}$) at high densities, and it should also fulfill the causality limit ($c_s/c<1$). However, its detailed behavior remains a hot topic of debate. We explore the general properties of the sound velocity and the adiabatic index in hybrid stars, as well as in neutron stars and quark stars. For this purpose, the bag model, the perturbation model, the equivparticle model, and the quasiparticle model are employed for the quark phase. One representative effective field theory model is used for the hardon phase. Various conditions are employed for hadron-quark phase transition and we employ an interface tension in a preferred range of $1-50~\rm MeV/fm^2$. The results are compared with various ab-initio calculations. We find a characteristic behavior of dynamical rescaling of the bag constant on the sound velocity in quark matter, which resembles that of quark deconfinement phase transition at $\sim3-7$ times of nuclear saturation density. And it leads to a much compact star with similar mass. We also propose a new class of quark star equation of states attributing the feature. The quark star equation of state model, as well as the in-medium scaling of the bag parameter in our calculation can be tested by future high-precision radius measurements of pulsar-like objects.

We study the interface effects of quark-hadron mixed phase in compact stars. The properties of nuclear matter are obtained based on the relativistic-mean-field model. For the quark phase, we adopt perturbation model with running quark masses and coupling constant. At certain choices of parameter sets, it is found that varying the quark-hadron interface tension will have sizable effects on the radii ($\Delta R \approx 600$ m) and tidal deformabilities ($\Delta \Lambda/\Lambda \approx 50\%$) of hybrid stars. These provide possibilities for us to constrain the quark-hadron interface tension with future gravitational wave observations as well as the ongoing NICER mission.

We conduct forward-modeling analysis based on our 2.5 dimensional magnetohydrodynamics (MHD) simulation of magnetic flux rope (MFR) formation and eruption driven by photospheric converging motion. The current sheet (CS) evolution during the MFR formation and eruption process in our MHD simulation can be divided into four stages. The current sheet (CS) evolution during the MFR formation and eruption process in our MHD simulation can be divided into four stages. The first stage shows the CS forming and gradually lengthening. Resistive instabilities that disrupt the CS mark the beginning of the second stage. Magnetic islands disappear in the third stage and reappear in the fourth stage. Synthetic images and light curves of the seven Solar Dynamics Observatory/Atmospheric Imaging Assembly (AIA) channels, i.e., 9.4 nm, 13.1 nm, 17.1 nm, 19.3 nm, 21.1 nm, 30.4 nm, and 33.5 nm, and the 3-25 keV thermal X-ray are obtained with forward-modeling analysis. The loop-top source and the coronal sources of the soft X-ray are reproduced in forward modeling. The light curves of the seven SDO/AIA channels start to rise once resistive instabilities develop. The light curve of the 3-25 keV thermal X-ray starts to go up when the reconnection rate reaches one of its peaks. Quasiperiodic pulsations (QPPs) appear twice in the SDO/AIA 17.1 nm, 21.1 nm, and 30.4 nm channels, corresponding to the period of chaotic (re)appearance and CS-guided displacements of the magnetic islands. QPPs appear once in the SDO/AIA 9.4 nm and 33.5 nm channels after the disruption of the CS by resistive instabilities and in the 19.3 nm channel when the chaotic motion of the magnetic islands reappears.

The interface effects play important roles for the properties of strange quark matter (SQM) and the related physical processes. We show several examples on the implications of interface effects for both stable and unstable SQM. Based on an equivparticle model and adopting mean-field approximation (MFA), the surface tension and curvature term of SQM can be obtained, which are increasing monotonically with the density of SQM at zero external pressure. For a parameter set constrained according to the 2$M_\odot$ strange star, we find the surface tension is $\sim$2.4 MeV/fm${}^2$, while it is larger for other cases.

We analyze particle acceleration in explosive reconnection events in magnetically dominated proton-electron plasmas. Reconnection is driven by large-scale magnetic stresses in interacting current-carrying flux tubes. Our model relies on development of current-driven instabilities on macroscopic scales. These tilt-kink instabilities develop in an initially force-free equilibrium of repelling current channels. Using MHD methods we study a 3D model of repelling and interacting flux tubes in which we simultaneously evolve test particles, guided by electromagnetic fields obtained from MHD. We identify two stages of particle acceleration; Initially particles accelerate in the current channels, after which the flux ropes start tilting and kinking and particles accelerate due to reconnection processes in the plasma. The explosive stage of reconnection produces non-thermal energy distributions with slopes that depend on plasma resistivity and the initial particle velocity. We also discuss the influence of the length of the flux ropes on particle acceleration and energy distributions. This study extends previous 2.5D results to 3D setups, providing all ingredients needed to model realistic scenarios like solar flares, black hole flares and particle acceleration in pulsar wind nebulae: formation of strong resistive electric fields, explosive reconnection and non-thermal particle distributions. By assuming initial energy equipartition between electrons and protons, applying low resistivity in accordance with solar corona conditions and limiting the flux rope length to a fraction of a solar radius we obtain realistic energy distributions for solar flares with non-thermal power law tails and maximum electron energies up to 11 MeV and maximum proton energies up to 1 GeV.

We compare relativistic particle integrators commonly used in plasma physics showing several test cases relevant for astrophysics. Three explicit particle pushers are considered, namely the Boris, Vay, and Higuera-Cary schemes. We also present a new relativistic fully implicit particle integrator that is energy conserving. Furthermore, a method based on the relativistic guiding center approximation is included. The algorithms are described such that they can be readily implemented in magnetohydrodynamics codes or Particle-in-Cell codes. Our comparison focuses on the strengths and key features of the particle integrators. We test the conservation of invariants of motion, and the accuracy of particle drift dynamics in highly relativistic, mildly relativistic, and non-relativistic settings. The methods are compared in idealized test cases, i.e., without considering feedback on the electrodynamic fields, collisions, pair creation, or radiation. The test cases include uniform electric and magnetic fields, $\mathbf{E}\times\mathbf{B}$-fields, force-free fields, and setups relevant for high-energy astrophysics, e.g., a magnetic mirror, a magnetic dipole, and a magnetic null. These tests have direct relevance for particle acceleration in shocks and in magnetic reconnection.

Employing Solar Dynamics Observatory/Atmospheric Imaging Assembly (AIA) multi-wavelength images, we have presented coronal condensations caused by magnetic reconnection between a system of open and closed solar coronal loops. In this Letter, we report the quasi-periodic fast magnetoacoustic waves propagating away from the reconnection region upward across the higher-lying open loops during the reconnection process. On 2012 January 19, reconnection between the higher-lying open loops and lower-lying closed loops took place, and two sets of newly reconnected loops formed. Thereafter, cooling and condensations of coronal plasma occurred in the magnetic dip region of higher-lying open loops. During the reconnection process, disturbances originating from the reconnection region propagate upward across the magnetic dip region of higher-lying loops with the mean speed and mean speed amplitude of 200 and 30 km s$^{-1}$, respectively. The mean speed of the propagating disturbances decreases from $\sim$230 km s$^{-1}$ to $\sim$150 km s$^{-1}$ during the coronal condensation process, and then increases to $\sim$220 km s$^{-1}$. This temporal evolution of the mean speed anti-correlates with the light curves of the AIA 131 and 304 \AA~channels that show the cooling and condensation process of coronal plasma. Furthermore, the propagating disturbances appear quasi-periodically with a peak period of 4 minutes. Our results suggest that the disturbances represent the quasi-periodic fast propagating magnetoacoustic (QFPM) waves originating from the magnetic reconnection between coronal loops.

We explore the effects of strangeness and $\Delta$ resonance in baryonic matter and compact stars within the relativistic-mean-field (RMF) models. The covariant density functional PKDD is adopted for $N$-$N$ interaction, parameters fixed based on finite hypernuclei and neutron stars are taken for the hyperon-meson couplings, and the universal baryon-meson coupling scheme is adopted for the $\Delta$-meson couplings. In light of the recent observations of GW170817 with the dimensionless combined tidal deformability $197 \leq \bar{\Lambda}\leq 720$, we find it is essential to include the $\Delta$ resonances in compact stars, and small $\Delta$-$\rho$ coupling $g_{\rho \Delta}$ is favored if the mass $2.27{}_{-0.15}^{+0.17}\ M_\odot$ of PSR J2215+5135 is confirmed.

The combination of magnetohydrodynamic (MHD) simulation and multi-wavelength observations is an effective way to study mechanisms of magnetic flux rope eruption. We develop a data-driven MHD model using the zero-$\beta$ approximation. The initial condition is provided by nonlinear force-free field derived by the magneto-frictional method based on vector magnetic field observed by the Helioseismic and Magnetic Imager (HMI) aboard the Solar Dynamics Observatory (SDO). The bottom boundary uses observed time series of the vector magnetic field and the vector velocity derived by the Differential Affine Velocity Estimator for Vector Magnetograms (DAVE4VM). We apply the data-driven model to active region 11123 observed from 06:00 UT on 2011 November 11 to about 2 hours later. The evolution of the magnetic field topology coincides with the flare ribbons observed in the 304 and 1600 \AA $ $ wavebands by the Atmospheric Imaging Assembly. The morphology, propagation path, and propagation range of the flux rope are comparable with the observations in 304 \AA . We also find that a data-constrained boundary condition, where the bottom boundary is fixed to the initial values, reproduces a similar simulation result. This model can reproduce the evolution of a magnetic flux rope in its dynamic eruptive phase.

Hard X-ray (HXR) sources are frequently observed near the top of solar flare loops, and the emission is widely ascribed to bremsstrahlung. We here revisit an alternative scenario which stresses the importance of inverse Compton processes and the Kelvin- Helmholtz instability (KHI) proposed by Fang et al. (2016). This scenario adds a novel ingredient to the standard flare model, where evaporation flows from flare-impacted chromospheric foot-points interact with each other near the loop top and produce turbulence via KHI. The turbulence can act as a trapping region and as an efficient accelerator to provide energetic electrons, which scatter soft X-ray (SXR) photons to HXR photons via the inverse Compton mechanism. This paper focuses on the trigger of the KHI and the resulting turbulence in this new scenario. We perform a parameter survey to investigate the necessary ingredients to obtain KHI through interaction of chromospheric evaporation flows. When turbulence is produced in the loop apex, an index of -5/3 can be found in the spectra of velocity and magnetic field fluctuations. The KHI development and the generation of turbulence are controlled by the amount of energy deposited in the chromospheric foot-points and the time scale of its energy deposition, but typical values for M class flares show the KHI development routinely. Asymmetry of energy deposition determines the location where the turbulence is produced, and the synthesized SXR light curve shows a clear periodic signal related to the sloshing motion of the vortex pattern created by the KHI.

Employing Solar Dynamics Observatory (SDO)/Atmospheric Imaging Assembly (AIA) multi-wavelength images, we report the coronal condensation during the magnetic reconnection (MR) between a system of open and closed coronal loops. Higher-lying magnetically open structures, observed in AIA 171 A images above the solar limb, move downward and interact with the lower-lying closed loops, resulting in the formation of dips in the former. An X-type structure forms at the interface. The interacting loops reconnect and disappear. Two sets of newly-reconnected loops then form and recede from the MR region. During the MR process, bright emission appears sequentially in the AIA 131 A and 304 A channels repeatedly in the dips of higher-lying open structures. This indicates the cooling and condensation process of hotter plasma from ~0.9 MK down to ~0.6 MK, and then to ~0.05 MK, also supported by the light curves of the AIA 171 A, 131 A, and 304 A channels. The part of higher-lying open structures supporting the condensations participate in the successive MR. The condensations without support by underlying loops then rain back to the solar surface along the newly-reconnected loops. Our results suggest that the MR between coronal loops leads to the condensation of hotter coronal plasma and its downflows. MR thus plays an active role in the mass cycle of coronal plasma because it can initiate the catastrophic cooling and condensation. This underlines that the magnetic and thermal evolution has to be treated together and cannot be separated, even in the case of catastrophic cooling.

Solar prominences are subject to all kinds of perturbations during their lifetime, and frequently demonstrate oscillations. The study of prominence oscillations provides an alternative way to investigate their internal magnetic and thermal structures as the oscillation characteristics depend on their interplay with the solar corona. Prominence oscillations can be classified into longitudinal and transverse types. We perform three-dimensional ideal magnetohydrodynamic simulations of prominence oscillations along a magnetic flux rope, with the aim to compare the oscillation periods with those predicted by various simplified models and to examine the restoring force. We find that the longitudinal oscillation has a period of about 49 minutes, which is in accordance with the pendulum model where the field-ligned component of gravity serves as the restoring force. In contrast, the horizontal transverse oscillation has a period of about 10 minutes and the vertical transverse oscillation has a period of about 14 minutes, and both of them can be nicely fitted with a two-dimensional slab model. We also find that the magnetic tension force dominates most of the time in transverse oscillations, except for the first minute when magnetic pressure overwhelms.

We report on the development of MPI-AMRVAC version 2.0, which is an open-source framework for parallel, grid-adaptive simulations of hydrodynamic and magnetohydrodynamic (MHD) astrophysical applications. The framework now supports radial grid stretching in combination with adaptive mesh refinement (AMR). The advantages of this combined approach are demonstrated with one-dimensional, two-dimensional and three-dimensional examples of spherically symmetric Bondi accretion, steady planar Bondi-Hoyle-Lyttleton flows, and wind accretion in Supergiant X-ray binaries. Another improvement is support for the generic splitting of any background magnetic field. We present several tests relevant for solar physics applications to demonstrate the advantages of field splitting on accuracy and robustness in extremely low plasma $\beta$ environments: a static magnetic flux rope, a magnetic null-point, and magnetic reconnection in a current sheet with either uniform or anomalous resistivity. Our implementation for treating anisotropic thermal conduction in multi-dimensional MHD applications is also described, which generalizes the original slope limited symmetric scheme from 2D to 3D. We perform ring diffusion tests that demonstrate its accuracy and robustness, and show that it prevents the unphysical thermal flux present in traditional schemes. The improved parallel scaling of the code is demonstrated with 3D AMR simulations of solar coronal rain, which show satisfactory strong scaling up to 2000 cores. Other framework improvements are also reported: the modernization and reorganization into a library, the handling of automatic regression tests, the use of inline/online Doxygen documentation, and a new future-proof data format for input/output

We intend to investigate the underlying physics for the coronal rain phenomenon in a representative bipolar magnetic field, including the formation and the dynamics of coronal rain blobs. With the MPI-AMRVAC code, we performed three dimensional radiative magnetohydrodynamic (MHD) simulation with strong heating localized on footpoints of magnetic loops after a relaxation to quiet solar atmosphere. Progressive cooling and in-situ condensation starts at the loop top due to radiative thermal instability. The first large-scale condensation on the loop top suffers Rayleigh-Taylor instability and becomes fragmented into smaller blobs. The blobs fall vertically dragging magnetic loops until they reach low beta regions and start to fall along the loops from loop top to loop footpoints. A statistic study of the coronal rain blobs finds that small blobs with masses of less than 10^10 g dominate the population. When blobs fall to lower regions along the magnetic loops, they are stretched and develop a non-uniform velocity pattern with an anti-parallel shearing pattern seen to develop along the central axis of the blobs. Synthetic images of simulated coronal rain with Solar Dynamics Observatory Atmospheric Imaging Assembly well resemble real observations presenting dark falling clumps in hot channels and bright rain blobs in a cool channel. We also find density inhomogeneities during a coronal rain "shower", which reflects the observed multi-stranded nature of coronal rain.

Magnetic reconnection and non-thermal particle distributions associated with current-driven instabilities are investigated by means of resistive magnetohydrodynamics (MHD) simulations combined with relativistic test particle methods. We propose a system with two parallel, repelling current channels in an initially force-free equilibrium, as a simplified representation of flux ropes in a stellar magnetosphere. The current channels undergo a rotation and separation on Alfv\'enic timescales, forming secondary islands and (up to tearing unstable) current sheets in which non-thermal energy distributions are expected to develop. Using the recently developed particle module of our open-source grid-adaptive MPI-AMRVAC software, we simulate MHD evolution combined with test particle treatments in MHD snapshots. We explore under which plasma-$\beta$ conditions the fastest reconnection occurs in two-and-a-half dimensional (2.5D) scenarios and in these settings test particles are evolved. We quantify energy distributions, acceleration mechanisms, relativistic corrections to the particle equations of motion and effects of resistivity in magnetically dominated proton-electron plasmas. Due to large resistive electric fields and indefinite acceleration of particles in the infinitely long current channels, hard energy spectra are found in 2.5D configurations. Solutions to these numerical artifacts are proposed for both 2.5D setups and future 3D work. We discuss the magnetohydrodynamics of an additional kink instability in 3D setups and the expected effects on energy distributions. The obtained results hold as a proof-of-principle for test particle approaches in MHD simulations, relevant to explore less idealised scenarios like solar flares and more exotic astrophysical phenomena, like black hole flares, magnetar magnetospheres and pulsar wind nebulae.

Solar prominences are long-lived cool and dense plasma curtains in the hot and rarefied outer solar atmosphere or corona. The physical mechanism responsible for their formation and especially for their internal plasma circulation has been uncertain for decades. The observed ubiquitous down flows in quiescent prominences are difficult to interpret as plasma with high conductivity seems to move across horizontal magnetic field lines. Here we present three-dimensional numerical simulations of prominence formation and evolution in an elongated magnetic flux rope as a result of in-situ plasma condensations fueled by continuous plasma evaporation from the solar chromosphere. The prominence is born and maintained in a fragmented, highly dynamic state with continuous reappearance of multiple blobs and thread structures that move mainly downward dragging along mass-loaded field lines. The prominence plasma circulation is characterized by the dynamic balance between the drainage of prominence plasma back to the chromosphere and the formation of prominence plasma via continuous condensation. Plasma evaporates from the chromosphere, condenses into the prominence in the corona, and drains back to the chromosphere, establishing a stable chromosphere-corona plasma cycle. Synthetic images of the modeled prominence with the Solar Dynamics Observatory Atmospheric Imaging Assembly closely resemble actual observations, with many dynamical threads underlying an elliptical coronal cavity.

In this paper, we produce synthetic radio views on simulated flux ropes in the solar corona, where finite-beta magnetohydrodynamic (MHD) simulations serve to mimic the flux rope formation stages, as well as their stable endstates. These endstates represent twisted flux ropes where balancing Lorentz forces, gravity and pressure gradients determine the full thermodynamic variation throughout the flux rope. The obtained models are needed to quantify radiative transfer in radio bands, and allow us to contrast weak to strong magnetic field conditions. Field strengths of up to 100 G in the flux rope yield the radio views dominated by optically thin free-free emission. The forming flux rope shows clear morphological changes in its emission structure as it deforms from an arcade to a flux rope, both on disk and at the limb. For an active region filament channel with a field strength of up to 680 G in the flux rope, gyroresonance emission (from the third-fourth gyrolayers) can be detected and even dominates over free-free emission at the frequencies of up to 7 GHz. Finally, we also show synthetic views on a simulated filament embedded within a (weak-field) flux rope, resulting from an energetically consistent MHD simulation. For this filament, synthetic views at the limb show clear similarities with actual observations, and the transition from optically thick (below 10 GHz) to optically thin emission can be reproduced. On the disk, its dimension and temperature conditions are as yet not realistic enough to yield the observed radio brightness depressions.

Quasi-periodic propagating intensity disturbances have been observed in large coronal loops in EUV images over a decade, and are widely accepted to be slow magnetosonic waves. However, spectroscopic observations from Hinode/EIS revealed their association with persistent coronal upflows, making this interpretation debatable. We perform a 2.5D magnetohydrodynamic simulation to imitate the chromospheric evaporation and the following reflected patterns in a flare loop. Our model encompasses the corona, transition region, and chromosphere. We demonstrate that the quasi periodic propagating intensity variations captured by the synthesized \textit{Solar Dynamics Observatory}/Atmospheric Imaging Assembly (AIA) 131, 94~\AA~emission images match the previous observations well. With particle tracers in the simulation, we confirm that these quasi periodic propagating intensity variations consist of reflected slow mode waves and mass flows with an average speed of 310 km/s in an 80 Mm length loop with an average temperature of 9 MK. With the synthesized Doppler shift velocity and intensity maps of the \textit{Solar and Heliospheric Observatory}/Solar Ultraviolet Measurement of Emitted Radiation (SUMER) Fe XIX line emission, we confirm that these reflected slow mode waves are propagating waves.

We extend our earlier multidimensional, magnetohydrodynamic simulations of coronal rain occurring in magnetic arcades with higher resolution, grid-adaptive computations covering a much longer ($>6$ hour) timespan. We quantify how in-situ forming blob-like condensations grow along and across field lines and show that rain showers can occur in limit cycles, here demonstrated for the first time in 2.5D setups. We discuss dynamical, multi-dimensional aspects of the rebound shocks generated by the siphon inflows and quantify the thermodynamics of a prominence-corona-transition-region like structure surrounding the blobs. We point out the correlation between condensation rates and the cross-sectional size of loop systems where catastrophic cooling takes place. We also study the variations of the typical number density, kinetic energy and temperature while blobs descend, impact and sink into the transition region. In addition, we explain the mechanisms leading to concurrent upflows while the blobs descend. As a result, there are plenty of shear flows generated with relative velocity difference around 80 km s$^{-1}$ in our simulations. These shear flows are siphon flows set up by multiple blob dynamics and they in turn affect the deformation of the falling blobs. In particular, we show how shear flows can break apart blobs into smaller fragments, within minutes.

Observations revealed rich dynamics within prominences, the cool 10,000 K, macroscopic (sizes of order 100 Mm) "clouds" in the million degree solar corona. Even quiescent prominences are continuously perturbed by hot, rising bubbles. Since prominence matter is hundredfold denser than coronal plasma, this bubbling is related to Rayleigh-Taylor instabilities. Here we report on true macroscopic simulations well into this bubbling phase, adopting a magnetohydrodynamic description from chromospheric layers up to 30 Mm height. Our virtual prominences rapidly establish fully non-linear (magneto)convective motions where hot bubbles interplay with falling pillars, with dynamical details including upwelling pillars forming within bubbles. Our simulations show impacting Rayleigh-Taylor fingers reflecting on transition region plasma, ensuring that cool, dense chromospheric material gets mixed with prominence matter up to very large heights. This offers an explanation for the return mass cycle mystery for prominence material. Synthetic views at extreme ultraviolet wavelengths show remarkable agreement with observations, with clear indications of shear-flow induced fragmentations.

We present numerical simulations in 3D settings where coronal rain phenomena take place in a magnetic configuration of a quadrupolar arcade system. Our simulation is a magnetohydrodynamic simulation including anisotropic thermal conduction, optically thin radiative losses, and parametrised heating as main thermodynamical features to construct a realistic arcade configuration from chromospheric to coronal heights. The plasma evaporation from chromospheric and transition region heights eventually causes localised runaway condensation events and we witness the formation of plasma blobs due to thermal instability, that evolve dynamically in the heated arcade part and move gradually downwards due to interchange type dynamics. Unlike earlier 2.5D simulations, in this case there is no large scale prominence formation observed, but a continuous coronal rain develops which shows clear indications of Rayleigh-Taylor or interchange instability, that causes the denser plasma located above the transition region to fall down, as the system moves towards a more stable state. Linear stability analysis is used in the non-linear regime for gaining insight and giving a prediction of the system's evolution. After the plasma blobs descend through interchange, they follow the magnetic field topology more closely in the lower coronal regions, where they are guided by the magnetic dips.

We present a numerical study in resistive magnetohydrodynamics where the initial equilibrium configuration contains adjacent, oppositely directed, parallel current channels. Since oppositely directed current channels repel, the equilibrium is liable to an ideal magnetohydrodynamic tilt instability. This tilt evolution, previously studied in planar settings, involves two magnetic islands or fluxropes, which on Alfvenic timescales undergo a combined rotation and separation. This in turn leads to the creation of (near) singular current layers, posing severe challenges to numerical approaches. Using our open-source grid-adaptive MPI-AMRVAC software, we revisit the planar evolution case in compressible MHD, as well as its extension to 2.5D and full 3D scenarios. As long as the third dimension remains ignorable, pure tilt evolutions result which are hardly affected by out of plane magnetic field components. In all 2.5D runs, our simulations do show secondary tearing type disruptions throughout the near singular current sheets in the far nonlinear saturation regime. In full 3D, both current channels can be liable to additional ideal kink deformations. We discuss the effects of having both tilt and kink instabilities acting simultaneously in the violent, reconnection dominated evolution. In 3D, both the tilt and the kink instabilities can be stabilized by tension forces. As a concrete space plasma application, we argue that interacting tilt-kink instabilities in repelling current channels provide a novel route to initiate solar coronal mass ejections, distinctly different from currently favored pure kink or torus instability routes.

Prominences in the solar corona are hundredfold cooler and denser than their surroundings, with a total mass of 1.e13 up to 1.e15 g. Here we report on the first comprehensive simulations of three-dimensional, thermally and gravitationally stratified magnetic flux ropes, where in situ condensation to a prominence happens due to radiative losses. After a gradual thermodynamic adjustment, we witness a phase where runaway cooling happens while counter-streaming shearing flows drain off mass along helical field lines. After this drainage, a prominence-like condensation resides in concave upward field regions, and this prominence retains its overall characteristics for more than two hours. While condensing, the prominence establishes a prominence-corona transition region, where magnetic field-aligned thermal conduction is operative during the runaway cooling. The prominence structure represents a force-balanced state in a helical flux rope. The simulated condensation demonstrates a right-bearing barb, as a remnant of the drainage. Synthetic images at extreme ultraviolet wavelengths follow the onset of the condensation, and confirm the appearance of horns and a three-part structure for the stable prominence state, as often seen in erupting prominences. This naturally explains recent Solar Dynamics Observatory views with the Atmospheric Imaging Assembly on prominences in coronal cavities demonstrating horns.

In this paper we present an update on the open source MPI-AMRVAC simulation toolkit where we focus on solar- and non-relativistic astrophysical magneto-fluid dynamics. We highlight recent developments in terms of physics modules such as hydrodynamics with dust coupling and the conservative implementation of Hall magnetohydrodynamics. A simple conservative high-order finite difference scheme that works in combination with all available physics modules is introduced and demonstrated at the example of monotonicity preserving fifth order reconstruction. Strong stability preserving high order Runge-Kutta time steppers are used to obtain stable evolutions in multidimensional applications realizing up to fourth order accuracy in space and time. With the new distinction between active and passive grid cells, MPI-AMRVAC is ideally suited to simulate evolutions where parts of the solution are controlled analytically, or have a tendency to progress into or out of a stationary state. Typical test problems and representative applications are discussed, with an outlook to follow-up research. Finally, we discuss the parallel scaling of the code and demonstrate excellent weak scaling up to 30 000 processors allowing to exploit modern petascale infrastructure.

We present numerical simulations in 2.5D settings where large scale prominences form in situ out of coronal condensation in magnetic dips, in close agreement with early as well as recent reporting of `funnel prominences'. Our simulation uses full thermodynamic MHD with anisotropic thermal conduction, optically thin radiative losses, and parametrized heating as main ingredients to establish a realistic arcade configuration from chromosphere to corona. The chromospheric evaporation from especially transition region heights ultimately causes thermal instability and we witness the growth of a prominence suspended well above the transition region, continuously gaining mass and cross-sectional area. Several hours later, the condensation has grown into a structure connecting the prominence-corona transition region with the underlying transition region, and a continuous downward motion from the accumulated mass represents a drainage that matches observational findings. A more dynamic phase is found as well, with coronal rain, induced wave trains, and even a reconnection event when the core prominence plasma weighs down the fieldlines until a fluxrope gets formed. The upper part of the prominence is then trapped in a fluxrope structure, and we argue for its violent kink-unstable eruption as soon as the (ignored) length dimension would allow for ideal kink deformations.

The magnetic configuration hosting prominences and their surrounding coro- nal structure is a key research topic in solar physics. Recent theoretical and observational studies strongly suggest that a helical magnetic flux rope is an es- sential ingredient to fulfill most of the theoretical and observational requirements for hosting prominences. To understand flux rope formation details and obtain magnetic configurations suitable for future prominence formation studies, we here report on three-dimensional isothermal magnetohydrodynamic simulations including finite gas pressure and gravity. Starting from a magnetohydrostatic corona with a linear force-free bipolar magnetic field, we follow its evolution when introducing vortex flows around the main polarities and converging flows towards the polarity inversion line near the bottom of the corona. The con- verging flows bring feet of different loops together at the polarity inversion line and magnetic reconnection and flux cancellation happens. Inflow and outflow signatures of the magnetic reconnection process are identified, and the thereby newly formed helical loops wind around pre-existing ones so that a complete flux rope grows and ascends. When a macroscopic flux rope is formed, we switch off the driving flows and find that the system relaxes to a stable state containing a helical magnetic flux rope embedded in an overlying arcade structure. A major part of the formed flux rope is threaded by dipped field lines which can stably support prominence matter, while the total mass of the flux rope is in the order of 4-5.e14 g.

We present the first multidimensional, magnetohydrodynamic simulations which capture the initial formation and the long-term sustainment of the enigmatic coronal rain phenomenon. We demonstrate how thermal instability can induce a spectacular display of in-situ forming blob-like condensations which then start their intimate ballet on top of initially linear force-free arcades. Our magnetic arcades host chromospheric, transition region, and coronal plasma. Following coronal rain dynamics for over 80 minutes physical time, we collect enough statistics to quantify blob widths, lengths, velocity distributions, and other characteristics which directly match with modern observational knowledge. Our virtual coronal rain displays the deformation of blobs into $V$-shaped like features, interactions of blobs due to mostly pressure-mediated levitations, and gives the first views on blobs which evaporate in situ, or get siphoned over the apex of the background arcade. Our simulations pave the way for systematic surveys of coronal rain showers in true multidimensional settings, to connect parametrized heating prescriptions with rain statistics, ultimately allowing to quantify the coronal heating input.

It is found that both microflare-sized impulsive heating at one leg of the loop and a suddenly imposed velocity perturbation can propel the prominence to oscillate along the magnetic dip. An extensive parameter survey results in a scaling law, showing that the period of the oscillation, which weakly depends on the length and height of the prominence, and the amplitude of the perturbations, scales with $\sqrt{R/g_\odot}$, where $R$ represents the curvature radius of the dip, and $g_\odot$ is the gravitational acceleration of the Sun. This is consistent with the linear theory of a pendulum, which implies that the field-aligned component of gravity is the main restoring force for the prominence longitudinal oscillations, as confirmed by the force analysis. However, the gas pressure gradient becomes non-negligible for short prominences. The oscillation damps with time in the presence of non-adiabatic processes. Compared to heat conduction, the radiative cooling is the dominant factor leading to the damping. A scaling law for the damping timescale is derived, i.e., $\tau\sim l^{1.63} D^{0.66}w^{-1.21}v_{0}^{-0.30}$, showing strong dependence on the prominence length $l$, the geometry of the magnetic dip (characterized by the depth $D$ and the width $w$), and the velocity perturbation amplitude $v_0$. The larger the amplitude, the faster the oscillation damps. It is also found that mass drainage significantly reduces the damping timescale when the perturbation is too strong.

Filament longitudinal oscillations have been observed on the solar disk in H$\alpha$. We intend to find an example of the longitudinal oscillations of a prominence, where the magnetic dip can be seen directly, and examine what is the restoring force of such kind of oscillations. We carry out a multiwavelength data analysis of the active region prominence oscillations above the western limb on 2007 February 8. Besides, we perform a one-dimensional hydrodynamic simulation of the longitudinal oscillations. The high-resolution observations by Hinode/SOT indicate that the prominence, seen as a concave-inward shape in lower-resolution Extreme Ultraviolet (EUV) images, actually consists of many concave-outward threads, which is indicative of the existence of magnetic dips. After being injected into the dip region, a bulk of prominence material started to oscillate for more than 3.5 hours, with the period being 52 min. The oscillation decayed with time, with the decay timescale being 133 min. Our hydrodynamic simulation can well reproduce the oscillation period, but the damping timescale in the simulation is 1.5 times as long as the observations. The results clearly show the prominence longitudinal oscillations around the dip of the prominence and our study suggests that the restoring force of the longitudinal oscillations might be the gravity. Radiation and heat conduction are insufficient to explain the decay of the oscillations. Other mechanisms, such as wave leakage and mass accretion, have to be considered. The possible relation between the longitudinal oscillations and the later eruption of a prominence thread, as well as a coronal mass ejection (CME), is also discussed.

Starting from a realistically sheared magnetic arcade connecting chromospheric, transition region to coronal plasma, we simulate the in-situ formation and sustained growth of a quiescent prominence in the solar corona. Contrary to previous works, our model captures all phases of the prominence formation, including the loss of thermal equilibrium, its successive growth in height and width to macroscopic dimensions, and the gradual bending of the arched loops into dipped loops, as a result of the mass accumulation. Our 2.5-dimensional, fully thermodynamically and magnetohydrodynamically consistent model mimics the magnetic topology of normal-polarity prominences above a photospheric neutral line, and results in a curtain-like prominence above the neutral line through which the ultimately dipped magnetic field lines protrude at a finite angle. The formation results from concentrated heating in the chromosphere, followed by plasma evaporation and later rapid condensation in the corona due to thermal instability, as verified by linear instability criteria. Concentrated heating in the lower atmosphere evaporates plasma from below to accumulate at the top of coronal loops and supply mass to the later prominence constantly. This is the first evaporation-condensation model study where we can demonstrate how the formed prominence stays in a force balanced state, which can be compared to the Kippenhahn-Schluter type magnetohydrostatic model, all in a finite low-beta corona.

It has been established that cold plasma condensations can form in a magnetic loop subject to localized heating of the footpoints. In this paper, we use grid-adaptive numerical simulations of the radiative hydrodynamic equations to parametrically investigate the filament formation process in a pre-shaped loop with both steady and finite-time chromospheric heating. Compared to previous works, we consider low-lying loops with shallow dips, and use a more realistic description for the radiative losses. We demonstrate for the first time that the onset of thermal instability satisfies the linear instability criterion. The onset time of the condensation is roughly \sim 2 hr or more after the localized heating at the footpoint is effective, and the growth rate of the thread length varies from 800 km hr-1 to 4000 km hr-1, depending on the amplitude and the decay length scale characterizing this localized chromospheric heating. We show how single or multiple condensation segments may form in the coronal portion. In the asymmetric heating case, when two segments form, they approach and coalesce, and the coalesced condensation later drains down into the chromosphere. With a steady heating, this process repeats with a periodicity of several hours. While our parametric survey confirms and augments earlier findings, we also point out that steady heating is not necessary to sustain the condensation. Once the condensation is formed, it can keep growing also when the localized heating ceases. Finally, we show that the condensation can survive continuous buffeting by perturbations resulting from the photospheric p-mode waves.