We report a new evaluation of the accretion properties of PDS~70b obtained with VLT/MUSE. The main difference from previous studies in Haffert et al. (2019) and Aoyama & Ikoma (2019) is in the mass accretion rate. Simultaneous multiple line observations, such as H$\alpha$ and H$\beta$, can better constrain the physical properties of an accreting planet. While we clearly detected H$\alpha$ emissions from PDS~70b, no H$\beta$ emissions were detected. We estimate the line flux of H$\beta$ with a 3-$\sigma$ upper limit to be 2.3~$\times$~10$^{-16}$~erg~s$^{-1}$~cm$^{-2}$. The flux ratio $F_{\rm H\beta}$/$F_{\rm H\alpha}$ for PDS~70b is $<$~0.28. Numerical investigations by Aoyama et al. (2018) suggest that $F_{\rm H\beta}$/$F_{\rm H\alpha}$ should be close to unity if the extinction is negligible. We attribute the reduction of the flux ratio to the extinction, and estimate the extinction of H$\alpha$ ($A_{\rm H\alpha}$) for PDS~70b to be $>$~2.0~mag using the interstellar extinction value. %The expected $A_{\rm H\alpha}$ value in the gap of the protoplanetary disk at the PDS~70b location is 2.4~mag, which is consistent with the estimated extinction. By combining with the H$\alpha$ linewidth and the dereddening line luminosity of H$\alpha$, %we derive the PDS~70b dynamical mass and mass accretion rate to be \hashimotor{12~$\pm$~3~$M_{\rm Jup}$} and $\gtrsim$~5~$\times$~10$^{-7}$~$M_{\rm Jup}$~yr$^{-1}$, respectively. we derive the PDS~70b mass accretion rate to be $\gtrsim$~5~$\times$~10$^{-7}$~$M_{\rm Jup}$~yr$^{-1}$. The PDS~70b mass accretion rate is an order of magnitude larger than that of PDS~70. We found that the filling factor $f_{\rm f}$ (the fractional area of the planetary surface emitting H$\alpha$) is $\gtrsim$0.01, which is similar to the typical stellar value. The small value of $f_{\rm f}$ indicates that the H$\alpha$ emitting areas are localized at the surface of PDS~70b.

For a better understanding of magnetic field in the solar corona and dynamic activities such as flares and coronal mass ejections, it is crucial to measure the time-evolving coronal field and accurately estimate the magnetic energy. Recently, a new modeling technique called the data-driven coronal field model, in which the time evolution of magnetic field is driven by a sequence of photospheric magnetic and velocity field maps, has been developed and revealed the dynamics of flare-productive active regions. Here we report on the first qualitative and quantitative assessment of different data-driven models using a magnetic flux emergence simulation as a ground-truth (GT) data set. We compare the GT field with those reconstructed from the GT photospheric field by four data-driven algorithms. It is found that, at least, the flux rope structure is reproduced in all coronal field models. Quantitatively, however, the results show a certain degree of model dependence. In most cases, the magnetic energies and relative magnetic helicity are comparable to or at most twice of the GT values. The reproduced flux ropes have a sigmoidal shape (consistent with GT) of various sizes, a vertically-standing magnetic torus, or a packed structure. The observed discrepancies can be attributed to the highly non-force-free input photospheric field, from which the coronal field is reconstructed, and to the modeling constraints such as the treatment of background atmosphere, the bottom boundary setting, and the spatial resolution.

Protostellar flares are rapid magnetic energy release events associated with formation of hot plasma in protostars. In the previous models of protostellar flares, the interaction between a protostellar magnetosphere with the surrounding disk plays crucial roles in building-up and releasing the magnetic energy. However, it remains unclear if protostars indeed have magnetospheres because vigorous disk accretion and strong disk magnetic fields in the protostellar phase may destroy the magnetosphere. Considering this possibility, we investigate the energy accumulation and release processes in the absence of a magnetosphere using a three-dimensional magnetohydrodynamic simulation. Our simulation reveals that protostellar flares are repeatedly produced even in such a case. Unlike in the magnetospheric models, the protostar accumulates magnetic energy by acquiring large-scale magnetic fields from the disk by accretion. Protostellar flares occur when a portion of the large-scale magnetic fields are removed from the protostar as a result of magnetic reconnection. Protostellar flares in the simulation are consistent with observations; the released magnetic energy (up to $\sim 3\times 10^{38}$ erg) is large enough to drive observed flares, and the flares produce hot ejecta. The expelled magnetic fields enhance accretion, and the energy build-up and release processes are repeated as a result. The magnetic flux removal via reconnection leads to redistribution of magnetic fields in the inner disk. We therefore consider that protostellar flares will play an important role in the evolution of the disk magnetic fields in the vicinity of protostars.

A new Riemann solver is presented for the ideal magnetohydrodynamics (MHD) equations with the so-called Boris correction. The Boris correction is applied to reduce wave speeds, avoiding an extremely small timestep in MHD simulations. The proposed Riemann solver, Boris-HLLD, is based on the HLLD solver. As done by the original HLLD solver, (1) the Boris-HLLD solver has four intermediate states in the Riemann fan when left and right states are given, (2) it resolves the contact discontinuity, Alfv\'en waves, and fast waves, and (3) it satisfies all the jump conditions across shock waves and discontinuities except for slow shock waves. The results of a shock tube problem indicate that the scheme with the Boris-HLLD solver captures contact discontinuities sharply and it exhibits shock waves without any overshoot when using the minmod limiter. The stability tests show that the scheme is stable when $|u| \lesssim 0.5c$ for a low Alfv\'en speed ($V_A \lesssim c$), where $u$, $c$, and $V_A$ denote the gas velocity, speed of light, and Alfv\'en speed, respectively. For a high Alfv\'en speed ($V_A \gtrsim c$), where the plasma beta is relatively low in many cases, the stable region is large, $|u| \lesssim (0.6-1) c$. We discuss the effect of the Boris correction on physical quantities using several test problems. The Boris-HLLD scheme can be useful for problems with supersonic flows in which regions with a very low plasma beta appear in the computational domain.

We present the results of a global, three-dimensional magnetohydrodynamics simulation of an accretion disk with a rotating, weakly magnetized central star. The disk is threaded by a weak, large-scale poloidal magnetic field, and the central star has no strong stellar magnetosphere initially. Our simulation investigates the structure of the accretion flows from a turbulent accretion disk onto the star. The simulation reveals that fast accretion onto the star at high latitudes occurs even without a stellar magnetosphere. We find that the failed disk wind becomes the fast, high-latitude accretion as a result of angular momentum exchange mediated by magnetic fields well above the disk, where the Lorentz force that decelerates the rotational motion of gas can be comparable to the centrifugal force. Unlike the classical magnetospheric accretion scenario, fast accretion streams are not guided by magnetic fields of the stellar magnetosphere. Nevertheless, the accretion velocity reaches the free-fall velocity at the stellar surface due to the efficient angular momentum loss at a distant place from the star. This study provides a possible explanation why Herbig Ae/Be stars whose magnetic fields are generally not strong enough to form magnetospheres also show indications of fast accretion. A magnetically driven jet is not formed from the disk in our model. The differential rotation cannot generate sufficiently strong magnetic fields for the jet acceleration because the Parker instability interrupts the field amplification.

A flare kernel associated with a C4 class flare was observed in a spectral window including the He I triplet 10830 A and Si I 10827 A with a spectropolarimeter on the Domeless Solar Telescope at Hida Observatory on August 9th, 2015. Observed Stokes profiles of the He I triplet in the flare kernel in its post-maximum phase are well reproduced through inversions considering the Zeeman and the Paschen-Back effects with a three slab model of the flare kernel, in which two slabs having up and downward velocities produce emissions and one slab produces an absorption. The magnetic field strength inferred from the emission components of the He I line is 1400 G, which is significantly stronger than 690 G that is observed at the same location in the same line 6.5 hours before the flare. In addition, photospheric magnetic field vector derived from the Si I10827 A is similar to that of the flare kernel. To explain this results, we suggest that the emission in the He I triplet during the flare is produced in the deep layer, around which bombardment of non-thermal electrons leads to the formation of a coronal temperature plasma. Assuming a hydrogen column density at the location where the He I emissions are formed, and a power-law index of non-thermal electron energy distribution, we derived the low-energy cutoff of the non-thermal electron as 20 - 30 keV, which is consistent with that inferred from hard X-ray data obtained by RHESSI.

Solar flare emission is detected in all EM bands and variations in flux density of solar energetic particles. Often the EM radiation generated in solar and stellar flares shows a pronounced oscillatory pattern, with characteristic periods ranging from a fraction of a second to several minutes. These oscillations are referred to as quasi-periodic pulsations (QPPs), to emphasise that they often contain apparent amplitude and period modulation. We review the current understanding of quasi-periodic pulsations in solar and stellar flares. In particular, we focus on the possible physical mechanisms, with an emphasis on the underlying physics that generates the resultant range of periodicities. These physical mechanisms include MHD oscillations, self-oscillatory mechanisms, oscillatory reconnection/reconnection reversal, wave-driven reconnection, two loop coalescence, MHD flow over-stability, the equivalent LCR-contour mechanism, and thermal-dynamical cycles. We also provide a histogram of all QPP events published in the literature at this time. The occurrence of QPPs puts additional constraints on the interpretation and understanding of the fundamental processes operating in flares, e.g. magnetic energy liberation and particle acceleration. Therefore, a full understanding of QPPs is essential in order to work towards an integrated model of solar and stellar flares.

Solar active regions (ARs) that produce strong flares and coronal mass ejections (CMEs) are known to have a relatively high non-potentiality and are characterized by delta-sunspots and sheared magnetic structures. In this study, we conduct a series of flux emergence simulations from the convection zone to the corona and model four types of active regions that have been observationally suggested to cause strong flares, namely the Spot-Spot, Spot-Satellite, Quadrupole, and Inter-AR cases. As a result, we confirm that delta-spot formation is due to the complex geometry and interaction of emerging magnetic fields, with finding that the strong-field, high-gradient, highly-sheared polarity inversion line (PIL) is created by the combined effect of the advection, stretching, and compression of magnetic fields. We show that free magnetic energy builds up in the form of a current sheet above the PIL. It is also revealed that photospheric magnetic parameters that predict flare eruptions reflect the stored free energy with high accuracy, while CME-predicting parameters indicate the magnetic relationship between flaring zones and entire ARs.

There is a subclass of the X-ray jets from young stellar objects which are heated very close to the footpoint of the jets, particularly DG Tau jets. Previous models attribute the strong heating to shocks in the jets. However, the mechanism that localizes the heating at the footpoint remains puzzling. We presented a different model of such X-ray jets, in which the disk atmosphere is magnetically heated. Our disk corona model is based on the so-called nanoflare model for the solar corona. We show that the magnetic heating near the disks can result in the formation of a hot corona with a temperature of > 10^6 K even if the average field strength in the disk is moderately weak, > 1 G. We determine the density and the temperature at the jet base by considering the energy balance between the heating and cooling. We derive the scaling relations of the mass loss rate and terminal velocity of jets. Our model is applied to the DG Tau jets. The observed temperature and estimated mass loss rate are consistent with the prediction of our model in the case of the disk magnetic field strength of ~20 G and the heating region of < 0.1 au. The derived scaling relation of the temperature of X-ray jets could be a useful tool to estimate the magnetic field strength. We also found that the jet X-ray can have a significant impact on the ionization degree near the disk surface and the dead-zone size.

Coronal disturbances associated with solar flares, such as H$\alpha$ Moreton waves, X-ray waves, and extreme ultraviolet (EUV) coronal waves are discussed herein in relation to magnetohydrodynamics fast-mode waves or shocks in the corona. To understand the mechanism of coronal disturbances, full-disk solar observations with high spatial and temporal resolution over multiple wavelengths are of crucial importance. We observed a filament eruption, whose shape is like a "dandelion", associated with the M1.6 flare that occurred on 2011 February 16 in the H$\alpha$ images taken by the Flare Monitoring Telescope at Ica University, Peru. We derive the three-dimensional velocity field of the erupting filament. We also identify winking filaments that are located far from the flare site in the H$\alpha$ images, whereas no Moreton wave is observed. By comparing the temporal evolution of the winking filaments with those of the coronal wave seen in the extreme ultraviolet images data taken by the Atmospheric Imaging Assembly on board the {\it Solar Dynamics Observatory} and by the Extreme Ultraviolet Imager on board the {\it Solar Terrestrial Relations Observatory-Ahead}, we confirm that the winking filaments were activated by the EUV coronal wave.

We report a strong association between the particle acceleration and plasma motions found in the 2010 August 18 solar flare. The plasma motions are tracked in the extreme-ultraviolet (EUV) images taken by the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory and the Extreme UltraViolet Imager (EUVI) on the Solar Terrestrial Relation Observatory spacecraft Ahead, and the signature of particle acceleration was investigated by using Nobeyama Radioheliograph data. In our previous paper, we reported that in EUV images many plasma blobs appeared in the current sheet above the flare arcade. They were ejected bidirectionally along the current sheet, and the blobs that were ejected sunward collided with the flare arcade. Some of them collided or merged with each other before they were ejected from the current sheet. We discovered impulsive radio bursts associated with such plasma motions (ejection, coalescence, and collision with the post flare loops). The radio bursts are considered to be the gyrosynchrotron radiation by nonthermal high energy electrons. In addition, the stereoscopic observation by AIA and EUVI suggests that plasma blobs had a three-dimensionally elongated structure. We consider that the plasma blobs were three-dimensional plasmoids (i.e. flux ropes) moving in a current sheet. We believe that our observation provides clear evidence of particle acceleration associated with the plasmoid motions. We discuss possible acceleration mechanisms on the basis of our results.

The role of slow-mode MHD shocks in magnetic reconnection is one of great importance for energy conversion and transport, but in many astrophysical plasmas the plasma is not fully ionised. In this paper, we investigate, using numerical simulations, the role of collisional coupling between a proton-electron charge-neutral fluid and a neutral hydrogen fluid for the 1D Riemann problem initiated in a constant pressure and density background state by a discontinuity in the magnetic field. This system, in the MHD limit, is characterised by two waves: a fast-mode rarefaction wave that drives a flow towards a slow-mode MHD shock. The system evolves through four stage: initiation, weak coupling, intermediate coupling and a quasi steady state. The initial stages are characterised by an over-pressured neutral region that expands with characteristics of a blast wave. In the later stages, the system tends towards a self-similar solution where the main drift velocity is concentrated in the thin region of the shock front. Due to the nature of the system, the neutral fluid is overpressured by the shock when compared to a purely hydrodynamic shock which results in the neutral fluid expanding to form the shock precursor. The thickness of the shockfront once it has formed proportional to the ionisation fraction to the power -1.2, which is a smaller exponent than would be naively expected from simple scaling arguments. One interesting result is that the shock front is a continuous transition of the physical variables for sub-sonic velocity upstream of the shock front (a c-shock) to a sharp jump in the physical variables followed by a relaxation to the downstream values for supersonic upstream velocity (a j-shock). The frictional heating that results from the velocity drift across the shock front can amount to approximately two per cent of the reference magnetic energy.

Recent space based observations of the Sun revealed that magnetic reconnection is ubiquitous in the solar atmosphere, ranging from small scale reconnection (observed as nanoflares) to large scale one (observed as long duration flares or giant arcades). Often the magnetic reconnection events are associated with mass ejections or jets, which seem to be closely related to multiple plasmoid ejections from fractal current sheet. The bursty radio and hard X-ray emissions from flares also suggest the fractal reconnection and associated particle acceleration. We shall discuss recent observations and theories related to the plasmoid-induced-reconnection and the fractal reconnection in solar flares, and their implication to reconnection physics and particle acceleration. Recent findings of many superflares on solar type stars that has extended the applicability of the fractal reconnection model of solar flares to much a wider parameter space suitable for stellar flares are also discussed.

Observations revealed that various kinds of oscillations are excited in solar flare regions. Quasi-periodic pulsations (QPPs) in the flare emissions are commonly observed in a wide range of wavelengths. Recent observations have found that fast-mode magnetohydrodynamic (MHD) waves are quasi-periodically emitted from some flaring sites (quasi-periodic propagating fast-mode magnetoacoustic waves; QPFs). Both of QPPs and QPFs imply a cyclic disturbance originating from the flaring sites. However, the physical mechanisms remain puzzling. By performing a set of two-dimensional MHD simulations of a solar flare, we discovered the local oscillation above the loops filled with evaporated plasma (above-the-loop-top region) and the generation of QPFs from such oscillating regions. Unlike all previous models for QPFs, our model includes essential physics for solar flares, such as magnetic reconnection, heat conduction, and chromospheric evaporation. We revealed that QPFs can be spontaneously excited by the above-the-loop-top oscillation. It was found that this oscillation is controlled by the backflow of the reconnection outflow. The new model revealed that flare loops and the above-the-loop-top region are full of shocks and waves, which is different from the previous expectations based on a standard flare model and previous simulations. In this paper, we will show the QPF generation process based on our new picture of flare loops and will briefly discuss a possible relationship between QPFs and QPPs. Our findings will change the current view of solar flares to a new view in which they are a very dynamic phenomenon with full of shocks and waves.

We carried out a magnetohydrodynamics simulation where a subsurface twisted kink-unstable flux tube emerges from the solar interior to the corona. Unlike the previous expectations based on the bodily emergence of a knotted tube, we found that the kinked tube can spontaneously form a complex quadrupole structure at the photosphere. Due to the development of the kink instability before the emergence, the magnetic twist at the kinked apex of the tube is greatly reduced, although the other parts of the tube is still strongly twisted. This leads to the formation of a complex quadrupole structure: a pair of the coherent, strongly twisted spots and a narrow complex bipolar pair between it. The quadrupole is formed by the submergence of a portion of emerged magnetic fields. This result is relevant for understanding of the origin of the complex multipolar $\delta$-spot regions that have a strong magnetic shear and emerge with polarity orientations not following Hale-Nicholson and Joy Laws.

Magnetic reconnection is an important physical process in various explosive phenomena in the universe. In the previous studies, it was found that fast re- connection takes place when the thickness of a current sheet becomes on the order of a microscopic length such as the ion larmor radius or the ion inertial length. In this study, we investigated the pinching process of a current sheet by the Lorentz force in a low-{\beta} plasma using one-dimensional magnetohydrodynam- ics (MHD) simulations. It is known that there is an exact self-similar solution for this problem that neglects gas pressure. We compared the non-linear MHD dynamics with the analytic self-similar solution. From the MHD simulations, we found that with the gas pressure included the implosion process deviates from the analytic self-similar solution as t {\rightarrow} t 0, where t 0 is the explosion time when the thickness of a current sheet of the analytic solution becomes 0. We also found a pair of MHD fast-mode shocks are generated and propagate after the formation of the pinched current sheet as t {\rightarrow} t 0 . On the basis of the Rankine-Hugoniot relations, we derived the scaling law of the physical quantities with respect to the initial plasma beta in the pinched current sheet. Our study could help us to estimate the physical quantities in the pinched current sheet formed in a low-\b{eta} plasma.

Solar flares are an explosive phenomenon, where super-sonic flows and shocks are expected in and above the post-flare loops. To understand the dynamics of post-flare loops, a two-dimensional magnetohydrodynamic (2D MHD) simulation of a solar flare has been carried out. We found new shock structures in and above the post-flare loops, which were not resolved in the previous work by Yokoyama and Shibata 2001. To study the dynamics of flows along the reconnected magnetic field, kinematics and energetics of the plasma are investigated along selected field lines. It is found that shocks are crucial to determine the thermal and flow structures in the post-flare loops. On the basis of the 2D MHD simulation, we have developed a new post-flare loop model which we call the pseudo-2D MHD model. The model is based on the 1D MHD equations, where all the variables depend on one space dimension and all the three components of the magnetic and velocity fields are considered. Our pseudo-2D model includes many features of the multi-dimensional MHD processes related to magnetic reconnection (particularly MHD shocks), which the previous 1D hydrodynamic models are not able to include. We compare the shock formation and energetics of a specific field line in the 2D calculation with those in our pseudo-2D MHD model, and we found that they give similar results. This model will allow us to study the evolution of the post-flare loops in a wide parameter space without expensive computational cost and without neglecting important physics associated with magnetic reconnection.

We studied the acceleration mechanisms of chromospheric jets associated with emerging flux using a two dimensional magnetohydrodynamic (MHD) simulation. We found that slow mode shock waves generated by magnetic reconnection in the chromosphere and the photosphere play key roles in the acceleration mechanisms of chromospheric jets. An important parameter is the height of magnetic reconnection. When magnetic reconnection takes place near the photosphere, the reconnection out- flow collides with the region where the plasma beta is much larger than unity. Then the plasma moves along a magnetic field. This motion generates a slow mode wave. The slow mode wave develops to a strong slow shock as it propagates upward. When the slow shock crosses the transition region, the transition region is lifted up. As a result, we obtain a chromospheric jet as the lifted transition region. When magnetic reconnection takes place in the upper chromosphere, the chromospheric plasma is accelerated due to the combination of the Lorentz force and the whip-like motion of magnetic field. We found that the chromospheric plasma is further accelerated through the interaction between the transition region (steep density gradient) and a slow shock emanating from the reconnection point. This is an MHD effect which has not been discussed before.

We report the simultaneous extreme ultraviolet observation of magnetic reconnection inflow and outflow in a flare on 2010 August 18 observed by the Atmospheric Imaging Assembly on board the {\it Solar Dynamics Observatory}. We found that during the rise phase of the flare, some plasma blobs appeared in the sheet structure above the hot loops. The plasma blobs were ejected bidirectionally along the sheet structure (outflow), at the same time as the threads visible in extreme ultraviolet images moved toward the sheet structure (inflow). The upward and downward ejection velocities are 220-460 km s$^{-1}$ and 250-280 km s$^{-1}$, respectively. The inflow speed changed from 90 km s$^{-1}$ to 12 km s$^{-1}$ in 5 minutes. By using these velocities, we estimated the nondimensional reconnection rate, which we found to vary during this period from 0.20 to 0.055. We also found that the plasma blobs in the sheet structure collided or merged with each other before they were ejected from the sheet structure. We hypothesize that the sheet structure is the current sheet and that these plasma blobs are plasmoids or magnetic islands, which could be important for understanding the dynamics of the reconnection region.