This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.

We perform magnetohydrodynamic (MHD) simulations to investigate the propagation of Alfven wave in magnetic chromosphere. We use the 1.5-dimensional expanding flux tube geometry setting and transverse perturbation at the bottom to generate the Alfven wave. Compared with previous studies, our expansion is that we include the radiative loss term introduced by Carlsson & Leenaarts (2012). We find that when an observation based transverse wave generator is applied, the spatial distribution of the time-averaged radiative loss profile in our simulation is consistent with that in the classic atmospheric model. In addition, the energy flux in the corona is larger than the required value for coronal heating in the quiet region. Our study shows that the Alfven wave driven model has the potential to explain chromospheric heating and transport enough energy to the corona simultaneously.

We report X-ray spectroscopic results for four giant solar flares occurred on 2005 September 7 (X17.0), 2005 September 8 (X5.4), 2005 September 9 (X6.2), and 2006 December 5 (X9.0), obtained from Earth albedo data with the X-ray imaging spectrometer (XIS) onboard Suzaku. The good energy resolution of the XIS (FWHM$\sim$100 eV) enables us to separate a number of line-like features and detect the underlying continuum emission. These features include Si He$\alpha$, Si Ly$\alpha$, S He$\alpha$, S Ly$\alpha$, Ar He$\alpha$, and Ca He$\alpha$ originating from solar flares as well as fluorescent Ar K$\alpha$ and Ar K$\beta$ from the Earth atmosphere. Absolute elemental abundances (X/H) averaged over the four flares are obtained to be $\sim$2.0 (Ca), $\sim$0.7 (Si), $\sim$0.3 (S), and $\sim$0.9 (Ar) at around flare peaks. This abundance pattern is similar to those of active stars' coronae showing inverse first ionization potential (i-FIP) effects, i.e., elemental abundances decrease with decreasing FIP with a turnover at the low end of FIP. The abundances are almost constant during the flares, with an exception of Si which increases by a factor of $\sim$2 in the decay phase. The evolution of the Si abundance is consistent with a picture that the i-FIP plasma originates from the chromosphere evaporation and then mixes with the surrounding low-FIP biased materials. Flare-to-flare abundance varied by a factor of 2, agreeing with past observations of solar flares. Finally, we emphasize that Earth albedo data acquired by X-ray astronomy satellites like Suzaku and XRISM can significantly contribute to studies of solar physics.

We succeeded in observing two large spicules simultaneously with the Atacama Large Millimeter/submillimeter Array (ALMA), the Interface Region Imaging Spectrograph (IRIS), and the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory. One is a spicule seen in the IRIS Mg II slit-jaw images and AIA 304\AA\ images (MgII/304A spicule). The other one is a spicule seen in the 100GHz images obtained with ALMA (100GHz spicule). Although the 100GHz spicule overlapped with the MgII/304A spicule in the early phase, it did not show any corresponding structures in the IRIS Mg II and AIA 304A images after the early phase. It suggests that the spicules are individual events and do not have a physical relationship. To obtain the physical parameters of the 100GHz spicule, we estimate the optical depths as a function of temperature and density using two different methods. One is using the observed brightness temperature by assuming a filling factor, and the other is using an emission model for the optical depth. As a result of comparing them, the kinetic temperature of the plasma and the number density of ionized hydrogens in the 100GHz spicule are ~6800 K and 2.2 x 10^10 cm^-3. The estimated values can explain the absorbing structure in the 193A image, which appear as a counterpart of the 100GHz spicule. These results suggest that the 100GHz spicule presented in this paper is classified to a macrospicule without a hot sheath in former terminology.

Using a three-dimensional compressible magnetohydrodynamic (MHD) simulation, we have reproduced the fast solar wind in a direct and self-consistent manner, based on the wave/turbulence driven scenario. As a natural consequence of Alfvenic perturbations at its base, highly compressional and turbulent fluctuations are generated, leading to heating and acceleration of the solar wind. The analysis of power spectra and structure functions reveals that the turbulence is characterized by its imbalanced (in the sense of outward Alfvenic fluctuations) and anisotropic nature. The density fluctuation originates from the parametric decay instability of outwardly propagating Alfven waves and plays a significant role in the Alfven wave reflection that triggers turbulence. Our conclusion is that the fast solar wind is heated and accelerated by compressible MHD turbulence driven by parametric decay instability and resultant Alfven wave reflection.

The understanding of magnetic reconnection in three-dimension (3D) is far shallower than its counterpart in two-dimension (2D) due to its potential complication, not to mention the evolving of the spontaneously growing turbulence. We investigate the reason for reconnection acceleration on the characters and development of diffusion regions and sheared 3D energy modes (energy modes that are not parallel to the anti-parallel magnetic fields) during the turbulence building stage. We found that multiple reconnection layers emerge due to the growth of 3D sheared tearing instability. Diffusion regions on adjacent reconnection layers form an inflow-outflow coupling that enhances the local reconnection. Further coupling of the existing energy modes breeds new energy modes near the current sheet edge. As reconnection layers span and interact with each other across the whole current sheet, global magnetic energy consumption accelerates. The significant contribution of 3D energy modes and their interaction to the reconnection rate enhancement seems to be independent on magnetic diffusivity. On the other hand, the global guide field changes the layout of 3D reconnection layer thus determines whether the system is fast-reconnection-preferable.

The structure of the photospheric vector magnetic field below a dark filament on the Sun is studied using the observations of the Spectro-Polarimeter attached to the Solar Optical Telescope onboard Hinode. Special attention is paid to discriminate the two suggested models, a flux rope or a bent arcade. "Inverse-polarity" orientation is possible below the filament in a flux rope, whereas "normal-polarity" can appear in both models. We study a filament in active region NOAA 10930, which appeared on the solar disk during 2006 December. The transverse field perpendicular to the line of sight has a direction almost parallel to the filament spine with a shear angle of 30 deg, whose orientation includes the 180-degree ambiguity. To know whether it is in the normal orientation or in the inverse one, the center-to-limb variation is used for the solution under the assumption that the filament does not drastically change its magnetic structure during the passage. When the filament is near the east limb, we found that the line-of-site magnetic component below is positive, while it is negative near the west limb. This change of sign indicates that the horizontal photospheric field perpendicular to the polarity inversion line beneath the filament has an "inverse-polarity", which indicates a flux-rope structure of the filament supporting field.

Chondrites are undifferentiated sediments of material left over from the earliest solar system and are widely considered as representatives of the unprocessed building blocks of the terrestrial planets. The chondrites, along with processed igneous meteorites, have been divided into two broad categories based upon their isotopic signatures; these have been termed the CC and NC groups and have been interpreted as reflecting their distinctive birth places within the solar system. The isotopic distinctiveness of NC and CC meteorites document limited radial-mixing in the accretionary disk. The enstatite and ordinary chondrites are NC-type and likely represent samples from inner solar system (likely $<$4 AU). Measurement and modeling of ratios of refractory lithophile elements (RLE) in enstatite chondrites establish these meteorites as the closest starting materials for the bulk of the silicate Earth and the core. Comparing chondritic and terrestrial RLE ratios demonstrate that the Bulk Silicate Earth, not the core, host the Earth's inventory of Ti, Zr, Nb, and Ta, but not the full complement of V.

The interiors of quiescent prominences are filled with turbulent flows. The evolution of upflow plumes, descending pillars, and vortex motions has been clearly detected in high-resolution observations. The Rayleigh-Taylor instability is thought to be a driver of such internal flows. Descending pillars are related to the mass budgets of prominences. There is a hypothesis of dynamic equilibrium where the mass drainage via descending pillars and the mass supply via radiative condensation are balanced to maintain the prominence mass; however, the background physics connecting the two different processes is poorly understood. In this study, we reproduced the dynamic interior of a prominence via radiative condensation and the mechanism similar to the Rayleigh-Taylor instability using a three-dimensional magnetohydrodynamic simulation including optically thin radiative cooling and nonlinear anisotropic thermal conduction. The process to prominence formation in the simulation follows the reconnection-condensation model, where topological change in the magnetic field caused by reconnection leads to radiative condensation. Reconnection is driven by converging motion at the footpoints of the coronal arcade fields. In contrast to the previous model, by randomly changing the speed of the footpoint motion along a polarity inversion line, the dynamic interior of prominence is successfully reproduced. We find that the mass condensation rate of the prominence is enhanced in the case with dynamic state. Our results support the observational hypothesis that the condensation rate is balanced with the mass drainage rate and suggest that a self-induced mass maintenance mechanism exists.

We report the results of the Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of the solar chromosphere on the southern polar limb. Coordinated observations with the Interface Region Imaging Spectrograph (IRIS) are also conducted. ALMA provided unprecedented high spatial resolution in the millimeter band ($\approx$ 2.0 arcsec) at 100 GHz frequency with a moderate cadence (20 s). The results are as follows: (1) The ALMA 100 GHz images show saw-tooth patterns on the limb, and a comparison with SDO/AIA 171\AA\ images shows a good correspondence of the limbs with each other. (2) The ALMA 100 GHz movie shows a dynamic thorn-like structure elongating from the saw-tooth patterns on the limb, with lengths reaching at least 8 arcsec, thus suggesting jet-like activity in the ALMA microwave range. These ALMA jets are in good correspondence with IRIS jet clusters. (3) A blob ejection event is observed. By comparing with the IRIS Mg II slit-jaw images, the trajectory of the blob is located along the spicular patterns.

We conduct a three-dimensional magnetohydrodynamic (MHD) simulation of the parametric decay instability of Alfv\'en waves and resultant compressible MHD turbulence, which is likely to develop in the solar wind acceleration region. Because of the presence of the mean magnetic field, the nonlinear stage is characterized by filament-like structuring and anisotropic cascading. By calculating the timescales of phase mixing and the evolution of Alfv\'en wave turbulence, we have found that the early nonlinear stage is dominated by phase mixing, while the later phase is dominated by imbalanced Alfv\'en wave turbulence. Our results indicate that the regions in the solar atmosphere with large density fluctuation, such as the coronal bottom and wind acceleration region, are heated by phase-mixed Alfv\'en waves, while the other regions are heated by Alfv\'en wave turbulence.

Using numerical simulations, we investigate the onset and suppression of parametric decay instability (PDI) in the solar wind, focusing on the suppression effect by the wind acceleration and expansion. We also discuss the role of PDI in generating the density fluctuation in the solar wind and the cross-helicity evolution up to $1 {\rm \ au}$. Numerical simulations from the coronal base to 1 au are conducted in a self-consistent manner; we take into account the feedback of wave energy and pressure in the background. Monochromatic waves with various injection frequencies $f_0$ are injected to discuss the suppression of PDI, while broadband waves are applied to compare the numerical results with observation. We find, via the simulations with monochromatic waves, that high-frequency ($f_0 \gtrsim10^{-3} {\rm \ Hz}$) Alfv\'en waves are subject to PDI. Meanwhile, the maximum growth rate of the PDI of low-frequency ($f_0 \lesssim10^{-4} {\rm \ Hz}$) Alfv\'en waves becomes negative due to acceleration and expansion effects. Medium-frequency ($f_0 \approx 10^{-3.5} {\rm \ Hz}$) Alfv\'en waves have a positive growth rate but do not show the signature of PDI up to 1 au because the growth rate is too small. The medium-frequency waves experience neither PDI nor reflection so they propagate through the solar wind most efficiently. This is called the frequency-filtering mechanism of the solar wind. The simulations with broadband waves indicate that the observed trend of the density fluctuation is well explained by the evolution of PDI while the observed cross-helicity evolution is in agreement with low-frequency wave propagation.

We propose a novel one-dimensional model that includes both shock and turbulence heating and qualify how these processes contribute to heating the corona and driving the solar wind. Compressible MHD simulations allow us to automatically consider shock formation and dissipation, while turbulent dissipation is modeled via a one-point closure based on Alfv\'en wave turbulence. Numerical simulations were conducted with different photospheric perpendicular correlation lengths $\lambda_0$, which is a critical parameter of Alfv\'en wave turbulence, and different root-mean-square photospheric transverse-wave amplitudes $\delta v_0$. For the various $\lambda_0$, we obtain a low-temperature chromosphere, high-temperature corona, and supersonic solar wind. Our analysis shows that turbulence heating is always dominant when $\lambda_0 \lesssim 1{\rm \ Mm}$. This result does not mean that we can ignore the compressibility because the analysis indicates that the compressible waves and their associated density fluctuations enhance the Alfv\'en wave reflection and therefore the turbulence heating. The density fluctuation and the cross helicity are strongly affected by $\lambda_0$, while the coronal temperature and mass loss rate depend weakly on $\lambda_0$.

Spicule oscillations involve high-frequency components with a typical period approximately corresponding to $40-50$ s. The typical time scale of the photospheric oscillation is a few minutes, and thus, the origin of this high-frequency component is not trivial. In this study, a one-dimensional numerical simulation is performed to demonstrate that the observed spicule oscillations originate from longitudinal-to-transverse mode conversion that occurs around the equipartition layer in the chromosphere. Calculations are conducted in a self-consistent manner with the exception of additional heating to maintain coronal temperature. The analyses indicate the following features: (1) mode conversion efficiently excites high-frequency transverse waves; (2) the typical period of the high-frequency waves corresponds to the sound-crossing time of the mode conversion region; and (3) simulated root-mean-square velocity of the high-frequency component is consistent with the observed value. These results indicate that the observation of spicule oscillation provides direct evidence of mode conversion in the chromosphere.

It has recently been recognized that the convective velocities achieved in the current solar convection simulations might be over-estimated. The newly-revealed effects of the prevailing small-scale magnetic field within the convection zone may offer possible solutions to this problem. The small-scale magnetic fields can reduce the convective amplitude of small-scale motions through the Lorentz-force feedback, which concurrently inhibits the turbulent mixing of entropy between upflows and downflows. As a result, the effective Prandtl number may exceed unity inside the solar convection zone. In this paper, we propose and numerically confirm a possible suppression mechanism of convective velocity in the effectively high-Prandtl number regime. If the effective horizontal thermal diffusivity decreases (the Prandtl number accordingly increases), the subadiabatic layer which is formed near the base of the convection zone by continuous depositions of low entropy transported by adiabatically downflowing plumes is enhanced and extended. The global convective amplitude in the high-Prandtl thermal convection is thus reduced especially in the lower part of the convection zone via the change in the mean entropy profile which becomes more subadiabatic near the base and less superadiabatic in the bulk.

This paper presents a three-dimensional simulation of chromospheric jets with twisted magnetic field lines. Detailed treatments of the photospheric radiative transfer and the equation of states allow us to model realistic thermal convection near the solar surface, which excites various MHD waves and produces chromospheric jets in the simulation. A tall chromospheric jet with a maximum height of 10--11 Mm and lifetime of 8--10 min is formed above a strong magnetic field concentration. The magnetic field lines are strongly entangled in the chromosphere, which helps the chromospheric jet to be driven by the Lorentz force. The jet exhibits oscillatory motion as a natural consequence of its generation mechanism. We also find that the produced chromospheric jet forms a cluster with a diameter of several Mm with finer strands. These results imply a close relationship between the simulated jet and solar spicules.

We propose a reconnection-condensation model in which topological change in a coronal magnetic field via reconnection triggers radiative condensation, thereby resulting in prominence formation. Previous observational studies have suggested that reconnection at a polarity inversion line of a coronal arcade field creates a flux rope that can sustain a prominence; however, they did not explain the origin of cool dense plasmas of prominences. Using three-dimensional magnetohydrodynamic simulations including anisotropic nonlinear thermal conduction and optically thin radiative cooling, we demonstrate that reconnection can lead not only to flux rope formation but also to radiative condensation under a certain condition. In our model, this condition is described by the Field length, which is defined as the scale length for thermal balance between radiative cooling and thermal conduction. This critical condition depends weakly on the artificial background heating. The extreme ultraviolet emissions synthesized with our simulation results have good agreement with observational signatures reported in previous studies.

MHD waves permeate the solar atmosphere and constitute potential coronal heating agents. Yet, the waves detected so far may be but a small subset of the true existing wave power. Detection is limited by instrumental constraints, but also by wave processes that localise the wave power in undetectable spatial scales. In this study we conduct 3D MHD simulations and forward modelling of standing transverse MHD waves in coronal loops with uniform and non-uniform temperature variation in the perpendicular cross-section. The observed signatures are largely dominated by the combination of the Kelvin-Helmholtz instability (KHI), resonant absorption and phase mixing. In the presence of a cross-loop temperature gradient we find that emission lines sensitive to the loop core catch different signatures than those more sensitive to the loop boundary and the surrounding corona, leading to an out-of-phase intensity modulation produced by the KHI mixing. Common signatures to all considered models include an intensity and loop width modulation at half the kink period, fine strand-like structure, a characteristic arrow-shaped structure in the Doppler maps, overall line broadening in time but particularly at the loop edges. For our model, most of these features can be captured with a spatial resolution of $0.33\arcsec$ and spectral resolution of 25~km~s$^{-1}$, although severe over-estimation of the line width is obtained. Resonant absorption leads to a significant decrease of the observed kinetic energy from Doppler motions over time, which is not recovered by a corresponding increase in the line width from phase mixing and the KHI motions. We estimate this hidden wave energy to be a factor of $5-10$ of the observed value.

Meteoritic abundances of r-process elements are analyzed to deduce the history of chemical enrichment by r-process from the beginning of disk formation to the present time in the solar vicinity, by combining the abundance information from short-lived radioactive nuclei such as 244Pu with that from stable r-process nuclei such as Eu. These two types of nuclei can be associated with one r-process event and cumulation of events till formation of the solar system, respectively. With help of the observed local star formation history, we deduce the chemical evolution of 244Pu and obtain three main results: (i) the last r-process event occurred 130-140 Myr before formation of the solar system, (ii) the present-day low 244Pu abundance as measured in deep sea reservoirs results from the low recent star formation rate compared to ~4.5 - 5 Gyr ago, and (iii) there were ~15 r-process events in the solar vicinity from formation of the Galaxy to the time of solar system formation and ~30 r-process events to the present time. Then, adopting a reasonable hypothesis that a neutron star merger is the r-process production site, we find that the ejected r-process elements are extensively spread out and mixed with interstellar matter with a mass of ~3.5 million solar masses, which is about 100 times larger than that for supernova ejecta. In addition, the event frequency of r-process production is estimated to be one per about 1400 core-collapse supernovae, which is identical to the frequency of neutron star mergers estimated from the analysis of stellar abundances.

The main object of the paper is to present the condition of the non-diffusive part of the Reynolds stress for driving the double-cell structure of the solar meridional circulation, which has been revealed by recent helioseismic observations. By conducting a set of mean-field hydrodynamic simulations, we confirm for the first time that the double-cell meridional circulation can be achieved along with the solar-like differential rotation when the Reynolds stress transports the angular momentum upward in the lower part and downward in the upper part of the convection zone. It is concluded that, in a stationary state, the accumulated angular momentum via the Reynolds stress in the middle layer is advected to both the upper and lower parts of the convection zone by each of the two meridional circulation cells, respectively.

In the highly structured solar corona, resonant absorption is an unavoidable mechanism of energy transfer from global transverse MHD waves to local azimuthal Alfv\'en waves. Due to its localised nature, a direct detection of this mechanism is extremely difficult. Yet, it is the leading theory explaining the observed fast damping of the global transverse waves. However, at odds with this theoretical prediction, recent observations indicate that in the low amplitude regime such transverse MHD waves can also appear decay-less, a yet unsolved phenomenon. Recent numerical work has shown that Kelvin-Helmholtz instabilities (KHI) often accompany transverse MHD waves. In this work, we combine 3D MHD simulations and forward modelling to show that for currently achieved spatial resolution and observed small amplitudes, an apparent decay-less oscillation is obtained. This effect results from the combination of periodic brightenings produced by the KHI and the coherent motion of the KHI vortices amplified by resonant absorption. Such effect is especially clear in emission lines forming at temperatures that capture the boundary dynamics rather than the core, and reflects the low damping character of the local azimuthal Alfv\'en waves resonantly coupled to the kink mode. Due to phase mixing, the detected period can vary depending on the emission line, with those sensitive to the boundary having shorter periods than those sensitive to the loop core. This allows to estimate the density contrast at the boundary.

Using one-dimensional numerical simulations, we study the elementary process of Alfv\'{e}n wave reflection in a uniform medium, including nonlinear effects. In the linear regime, Alfv\'{e}n wave reflection is triggered only by the inhomogeneity of the medium, whereas in the nonlinear regime, it can occur via nonlinear wave-wave interactions. Such nonlinear reflection (backscattering) is typified by decay instability. In most studies of decay instabilities, the initial condition has been a circularly polarized Alfv\'{e}n wave. In this study we consider a linearly polarized Alfv\'en wave, which drives density fluctuations by its magnetic pressure force. For generality, we also assume a broadband wave with a red-noise spectrum. In the data analysis, we decompose the fluctuations into characteristic variables using local eigenvectors, thus revealing the behaviors of the individual modes. Different from circular-polarization case, we find that the wave steepening produces a new energy channel from the parent Alfv\'en wave to the backscattered one. Such nonlinear reflection explains the observed increasing energy ratio of the sunward to the anti-sunward Alfv\'{e}nic fluctuations in the solar wind with distance against the dynamical alignment effect.

Spatially averaged (> 50'') EUV spectral lines in the transition region of solar quiet regions are known to be redshifted. Because the mechanism underlying this phenomenon is unclear, we require additional physical information on the lower corona for limiting the theoretical models. To acquire this information, we measured the Doppler shifts over a wide coronal temperature range (log T[K]=5.7--6.3) using the spectroscopic data taken by the Hinode EUV Imaging Spectrometer. By analyzing the data over the center-to-limb variations covering the meridian from the south to the north pole, we successfully measured the velocity to an accuracy of 3 km/s. Below log T[K] = 6.0, the Doppler shifts of the emission lines were almost zero with an error of 1--3 km/s; above this temperature, they were blueshifted with a gradually increasing magnitude, reaching - 6.3 +/- 2.1 km/s at log T[K]=6.25.

In this second paper in the series, we investigate occurrence frequencies of apparent unipolar processes, cancellation, and emergence of patch structures in quiet regions. Apparent unipolar events are considerably more frequent than cancellation and emergence as per our definition, which is consistent with Lamb et al. (2013). Furthermore, we investigate the frequency distributions of changes in flux during apparent unipolar processes are and found that they concentrate around the detection limit of the analysis. Combining these findings with the results of our previous paper, Iida et al. (2012), that merging and splitting are more dominant than emergence and cancellation, these results support the understanding that apparent unipolar processes are actually interactions with and among patches below the detection limit and that there still are numerous flux interactions between the flux range in this analysis and below the detection limit. We also investigate occurrence frequency distributions of flux decrease during cancellation. We found a relatively strong dependence, 2.48$\pm$0:26 as a power-law index. This strong dependence on flux is consistent with the model, which is suggested in the previous paper.

We investigate the effect of coronal temperature on the formation process of solar chromospheric jets using two-dimensional magnetohydrodynamic simulations of the region from the upper convection zone to the lower corona. We develop a new radiative magnetohydrodynamic code for the dynamic modeling of the solar atmosphere, employing a LTE equation of state, optically thick radiative loss in the photosphere, optically thin radiative loss in the chromosphere and the corona, and thermal conduction along the magnetic field lines. Many chromospheric jets are produced in the simulations by shock waves passing through the transition region. We find that these jets are projected farther outward when the coronal temperature is lower (similar to that in coronal holes) and shorter when the coronal temperature is higher (similar to that in active regions). When the coronal temperature is high, the deceleration of the chromospheric jets is consistent with the model in which deceleration is determined by the periodic chromospheric shock waves. However, when the coronal temperature is low, the gravitational deceleration becomes more important and the chromospheric jets approach ballistic motion.

In this paper we show that the phase mixing of continuum Alfv\'{e}n waves and/or continuum slow waves in magnetic structures of the solar atmosphere as, e.g., coronal arcades, can create the illusion of wave propagation across the magnetic field. This phenomenon could be erroneously interpreted as fast magnetosonic waves. The cross-field propagation due to phase mixing of continuum waves is apparent because there is no real propagation of energy across the magnetic surfaces. We investigate the continuous Alfv\'{e}n and slow spectra in 2D Cartesian equilibrium models with a purely poloidal magnetic field. We show that apparent superslow propagation across the magnetic surfaces in solar coronal structures is a consequence of the existence of continuum Alfv\'{e}n waves and continuum slow waves that naturally live on those structures and phase mix as time evolves. The apparent cross-field phase velocity is related to the spatial variation of the local Alfv\'{e}n/slow frequency across the magnetic surfaces and is slower than the Alfv\'{e}n/sound velocities for typical coronal conditions. Understanding the nature of the apparent cross-field propagation is important for the correct analysis of numerical simulations and the correct interpretation of observations.

Using a three-dimensional magnetohydrodynamic model, we simulate the magnetic reconnection in a single current sheet. We assume a finite guide field, a random perturbation on the velocity field and uniform resistivity. Our model enhances the reconnection rate relative to the classical Sweet-Parker model in the same configuration. The efficiency of magnetic energy conversion is increased by interactions between the multiple tearing layers coexisting in the global current sheet. This interaction, which forms a positive-feedback system, arises from coupling of the inflow and outflow regions in different layers across the current sheet. The coupling accelerates the elementary reconnection events, thereby enhancing the global reconnection rate. The reconnection establishes flux tubes along each tearing layer. Slow-mode shocks gradually form along the outer boundaries of these tubes, further accelerating the magnetic energy conversion.Such positive-feedback system is absent in two-dimensional simulation, three-dimensional reconnection without a guide field and a reconnection under a single perturbation mode. We refer to our model as the "shock-evoking positive-feedback" model.

Transverse magnetohydrodynamic (MHD) waves have been shown to be ubiquitous in the solar atmosphere and can in principle carry sufficient energy to generate and maintain the Sun's million-degree outer atmosphere or corona. However, direct evidence of the dissipation process of these waves and subsequent heating has not yet been directly observed. Here we report on high spatial, temporal, and spectral resolution observations of a solar prominence that show a compelling signature of so-called resonant absorption, a long hypothesized mechanism to efficiently convert and dissipate transverse wave energy into heat. Aside from coherence in the transverse direction, our observations show telltale phase differences around 180 degrees between transverse motions in the plane-of-sky and line-of-sight velocities of the oscillating fine structures or threads, and also suggest significant heating from chromospheric to higher temperatures. Comparison with advanced numerical simulations support a scenario in which transverse oscillations trigger a Kelvin-Helmholtz instability (KHI) at the boundaries of oscillating threads via resonant absorption. This instability leads to numerous thin current sheets in which wave energy is dissipated and plasma is heated. Our results provide direct evidence for wave-related heating in action, one of the candidate coronal heating mechanisms.

Transverse magnetohydrodynamic (MHD) waves are ubiquitous in the solar atmosphere and may be responsible for generating the Sun's million-degree outer atmosphere. However, direct evidence of the dissipation process and heating from these waves remains elusive. Through advanced numerical simulations combined with appropriate forward modeling of a prominence flux tube, we provide the observational signatures of transverse MHD waves in prominence plasmas. We show that these signatures are characterized by thread-like substructure, strong transverse dynamical coherence, an out-of-phase difference between plane-of-the-sky motions and LOS velocities, and enhanced line broadening and heating around most of the flux tube. A complex combination between resonant absorption and Kelvin-Helmholtz instabilities (KHI) takes place in which the KHI extracts the energy from the resonant layer and dissipates it through vortices and current sheets, which rapidly degenerate into turbulence. An inward enlargement of the boundary is produced in which the turbulent flows conserve the characteristic dynamics from the resonance, therefore guaranteeing detectability of the resonance imprints. We show that the features described in the accompanying paper (Okamoto et al. 2015) through coordinated Hinode and IRIS observations match well the numerical results.

We propose an in-situ formation model for inverse-polarity solar prominence and demonstrate it using self-consistent 2.5-dimensional magnetohydrodynamics simulations, including thermal conduction along magnetic fields and optically thin radiative cooling. The model enables us to form cool dense plasma clouds inside a flux rope by radiative condensation, which is regarded as an inverse-polarity prominence. Radiative condensation is triggered by changes in the magnetic topology, i.e., formation of the flux rope from the sheared arcade field, and by thermal imbalance due to the dense plasma trapped inside the flux rope. The flux rope is created by imposing converging and shearing motion on the arcade field. Either when the footpoint motion is in the anti-shearing direction or when heating is proportional to local density, the thermal state inside the flux rope becomes cooling-dominant, leading to radiative condensation. By controlling the temperature of condensation, we investigate the relationship between the temperature and density of prominence and derive a scaling formula for this relationship. This formula suggests that the proposed model reproduce the observed density of prominence, which is 10-100 times larger than the coronal density. Moreover, the time evolution of the extreme ultraviolet emission synthesized by combining our simulation results with the response function of the Solar Dynamics Observatory Atmospheric Imaging Assembly filters agrees with the observed temporal and spatial intensity shift among multi-wavelength during in-situ condensation.

In order to better understand the nature of active region outflows, the electron density was measured by using a density-sensitive line pair Fe xiv 264.78A/274.20A.Since coronal line profiles of the outflow region are composed of a major component with a Doppler shift of < 10 km s^-1 and a minor component (enhanced blue wing: EBW) blueshifted by up to 100 km s^-1, we extracted EBW from the line profiles through double-Gaussian fitting. We tried applying the simultaneous fitting to those two Fe xiv lines with several physical restrictions. Electron density for both components (n_Major and n_EBW, respectively) was calculated by referring to the theoretical intensity ratio as a function of electron density as per the CHIANTI database. We studied six locations in the outflow regions around NOAA AR10978. The average electron density was n_Major = 10^(9.16 +- 0.16) cm^-3 and n_EBW = 10^(8.74 +- 0.29) cm^-3. The magnitude relationship between n_Major and n_EBW was opposite in the eastern and western outflow regions. The column depth was also calculated for each component, which leads to the result that the outflows possess only a small fraction (~ 0.1) in the eastern region, while they dominate over the major component in the line profiles by a factor of five in the western region. When taking into account the extending coronal structures, the western region can be thought to represent the mass leakage. In contrast, we suggest a possibility that the eastern region actually contributes to the mass supply to coronal loops.

We investigate small-scale dynamo action in the solar convection zone through a series of high resolution MHD simulations in a local Cartesian domain with 1$R_\odot$ (solar radius) of horizontal extent and a radial extent from 0.715 to 0.96$R_\odot$. The dependence of the solution on resolution and diffusivity is studied. For a grid spacing of less than 350 km, the root mean square magnetic field strength near the base of the convection zone reaches 95% of the equipartition field strength (i.e. magnetic and kinetic energy are comparable). For these solutions the Lorentz force feedback on the convection velocity is found to be significant. The velocity near the base of the convection zone is reduced to 50% of the hydrodynamic one. In spite of a significant decrease of the convection velocity, the reduction in the enthalpy flux is relatively small, since the magnetic field also suppresses the horizontal mixing of the entropy between up- and downflow regions. This effect increases the amplitude of the entropy perturbation and makes convective energy transport more efficient. We discuss potential implications of these results for solar global convection and dynamo simulations.

We present a high-resolution, highly stratified numerical simulation of rotating thermal convection in a spherical shell. Our aim is to study in detail the processes that can maintain a near surface shear layer (NSSL) as inferred from helioseismology. Using the reduced speed of sound technique we can extend our global convection simulation to $0.99\,R_{\odot}$ and include near the top of our domain small-scale convection with short time scales that is only weakly influenced by rotation. We find the formation of a NSSL preferentially in high latitudes in the depth range $r=0.95-0.975R_\odot$. The maintenance mechanisms are summarized as follows. Convection under weak influence of rotation leads to Reynolds stresses that transport angular momentum radially inward in all latitudes. This leads to the formation of a strong poleward directed meridional flow and a NSSL, which is balanced in the meridional plane by forces resulting from the $\langle v'_r v'_\theta\rangle$ correlation of turbulent velocities. The origin of the required correlations depends to some degree on latitude. In high latitudes a positive correlation $\langle v'_rv'_\theta\rangle$ is induced in the NSSL by the poleward meridional flow whose amplitude increases with the radius, while a negative correlation is generated by the Coriolis force in bulk of the convection zone. In low latitudes a positive correlation $\langle v'_rv'_\theta\rangle$ results from rotationally aligned convection cells ("banana cells"). The force caused by these Reynolds stresses is in balance with the Coriolis force in the NSSL.

We investigate the triggering mechanisms of plasma eruptions in the solar atmosphere due to interactions between emerging flux and coronal arcade fields by using two-dimensional MHD simulations. We perform parameter surveys with respect to arcade field height, magnetic field strength, and emerging flux location. Our results show that two possible mechanisms exist, and which mechanism is dominant depends mostly on emerging flux location. One mechanism appears when the location of emerging flux is close to the polarity inversion line (PIL) of an arcade field. This mechanism requires reconnection between the emerging flux and the arcade field, as pointed out by previous studies. The other mechanism appears when the location of emerging flux is around the edge of an arcade field. This mechanism does not require reconnection between the emerging flux and the arcade field but does demand reconnection in the arcade field above the PIL. Furthermore, we found that the eruptive condition for this mechanism can be represented by a simple formula.

Solar active regions (ARs) are thought to be formed by magnetic fields from the convection zone. Our flux emergence simulations revealed that a strong horizontal divergent flow (HDF) of unmagnetized plasma appears at the photosphere before the flux begins to emerge. In our earlier study, we analyzed HMI data for a single AR and confirmed presence of this precursor plasma flow in the actual Sun. In this paper, as an extension of our earlier study, we conducted a statistical analysis of the HDFs to further investigate their characteristics and better determine the properties. From SDO/HMI data, we picked up 23 flux emergence events over a period of 14 months, the total flux of which ranges from 10^{20} to 10^{22} Mx. Out of 23 selected events, 6 clear HDFs were detected by the method we developed in our earlier study, and 7 HDFs detected by visual inspection were added to this statistic analysis. We found that the duration of the HDF is on average 61 minutes and the maximum HDF speed is on average 3.1 km s^{-1}. We also estimated the rising speed of the subsurface magnetic flux to be 0.6-1.4 km s^{-1}. These values are highly consistent with our previous one-event analysis as well as our simulation results. The observation results lead us to the conclusion that the HDF is rather a common feature in the earliest phase of AR emergence. Moreover, our HDF analysis has capability of determining the subsurface properties of emerging fields that cannot be directly measured.

The Nice model of the dynamical instability and migration of the giant planets can explain many properties of the present Solar System, and can be used to constrain its early architecture. In the jumping-Jupiter version of the Nice model, required from the terrestrial planet constraint and dynamical structure of the asteroid belt, Jupiter has encounters with an ice giant. Here we study the survival of the Galilean satellites in the jumping-Jupiter model. This is an important concern because the ice-giant encounters, if deep enough, could dynamically perturb the orbits of the Galilean satellites, and lead to implausible results. We performed numerical integrations where we tracked the effect of planetary encounters on the Galilean moons. We considered three instability cases from Nesvorny & Morbidelli (2012) that differed in the number and distribution of encounters. We found that in one case, where the number of close encounters was relatively small, the Galilean satellite orbits were not significantly affected. In the other two, the orbital eccentricities of all moons were excited by encounters, Callisto's semimajor axis changed, and, in a large fraction of trials, the Laplace resonance of the inner three moons was disrupted. The subsequent evolution by tides damps eccentricities and can recapture the moons in the Laplace resonance. A more important constraint is represented by the orbital inclinations of the moons, which can be excited during the encounters and not appreciably damped by tides. We find that one instability case taken from Nesvorny & Morbidelli (2012) clearly fails this constraint. This shows how the regular satellites of Jupiter can be used to set limits on the properties of encounters in the jumping-Jupiter model, and help us to better understand how the early Solar System evolved.

Current analytical and numerical modelling suggest the existence of ubiquitous thin current sheets in the corona that could explain the observed heating requirements. On the other hand, new high resolution observations of the corona indicate that its magnetic field may tend to organise itself in fine strand-like structures of few hundred kilometres widths. The link between small structure in models and the observed widths of strand-like structure several orders of magnitude larger is still not clear. A popular theoretical scenario is the nanoflare model, in which each strand is the product of an ensemble of heating events. Here, we suggest an alternative mechanism for strand generation. Through forward modelling of 3D MHD simulations we show that small amplitude transverse MHD waves can lead in a few periods time to strand-like structure in loops in EUV intensity images. Our model is based on previous numerical work showing that transverse MHD oscillations can lead to Kelvin-Helmholtz instabilities that deform the cross-sectional area of loops. While previous work has focused on large amplitude oscillations, here we show that the instability can occur even for low wave amplitudes for long and thin loops, matching those presently observed in the corona. We show that the vortices generated from the instability are velocity sheared regions with enhanced emissivity hosting current sheets. Strands result as a complex combination of the vortices and the line-of-sight angle, last for timescales of a period and can be observed for spatial resolutions of a tenth of loop radius.

We carry out non-rotating high-resolution calculations of the solar global convection, which resolve convective scales of less than 10 Mm. To cope with the low Mach number conditions in the lower convection zone, we use the reduced speed of sound technique (RSST), which is simple to implement and requires only local communication in the parallel computation. In addition, the RSST allows us to expand the computational domain upward to about $0.99 R_{\odot}$ as it can also handle compressible flows. Using this approach, we study the solar convection zone on the global scale, including small-scale near-surface convection. In particular, we investigate the influence of the top boundary condition on the convective structure throughout the convection zone as well as on small-scale dynamo action. Our main conclusions are: 1. The small-scale downflows generated in the near-surface layer penetrate into deeper layers to some extent and excite small-scale turbulence in the region $>0.9R_\odot$, where $R_\odot$ is the solar radius. 2. In the deeper convection zone ($<0.9R_\odot$), the convection is not influenced by the location of the upper boundary. 3. Using an LES approach, we can achieve small-scale dynamo action and maintain a field of about $0.15-0.25 B_\mathrm{eq}$ throughout the convection zone, where $B_\mathrm{eq}$ is the equipartition magnetic field to the kinetic energy. 4. The overall dynamo efficiency varies significantly in the convection zone as a consequence of the downward directed Poynting flux and the depth variation of the intrinsic convective scales.

In this Letter we present a seismological detection of a rising motion of magnetic flux in the shallow convection zone of the Sun, and show estimates of the emerging speed and its decelerating nature. In order to evaluate the speed of subsurface flux that creates an active region, we apply six Fourier filters to the Doppler data of NOAA AR 10488, observed with SOHO/MDI, to detect the reduction of acoustic power at six different depths from -15 to -2 Mm. All the filtered acoustic powers show reductions, up to 2 hours before the magnetic flux first appears at the visible surface. The start times of these reductions show a rising trend with a gradual deceleration. The obtained velocity is first several km s^-1 in a depth range of 15--10 Mm, then ~1.5 km s^-1 at 10-5 Mm, finally ~0.5 km s^-1 at 5-2 Mm. If we assume that the power reduction is actually caused by the magnetic field, the velocity of order of 1 km s^-1 is well in accordance with previous observations and numerical studies. Moreover, the gradual deceleration strongly supports the theoretical model that the emerging flux slows down in the uppermost convection zone before it expands into the atmosphere to build an active region.

Solar active regions are formed through the emergence of magnetic flux from the deeper convection zone. Recent satellite observations have shown that a horizontal divergent flow (HDF) stretches out over the solar surface just before the magnetic flux appearance. The aims of this study are to investigate the driver of the HDF and to see the dependency of the HDF on the parameters of the magnetic flux in the convection zone. We conduct three-dimensional magnetohydrodynamic (3D MHD) numerical simulations of the magnetic flux emergence and vary the parameters in the initial conditions. An analytical approach is also taken to explain the dependency. The horizontal gas pressure gradient is found to be the main driver of the HDF. The maximum HDF speed shows positive correlations with the field strength and twist intensity. The HDF duration has a weak relation with the twist, while it shows negative dependency on the field strength only in the case of the stronger field regime. Parametric dependencies analyzed in this study may allow us to probe the structure of the subsurface magnetic flux by observing properties of the HDF.

We studied electron spectral indices of nonthermal emissions seen in hard X-rays (HXRs) and in microwaves. We analyzed 12 flares observed by the Hard X-ray Telescope aboard {\it Yohkoh}, Nobeyama Radio Polarimeters (NoRP), and the Nobeyama Radioheliograph (NoRH), and compared the spectral indices derived from total fluxes of hard X-rays and microwaves. Except for four events, which have very soft HXR spectra suffering from the thermal component, these flares show a gap $\Delta\delta$ between the electron spectral indices derived from hard X-rays $\delta_{X}$ and those from microwaves $\delta_{\mu}$ ($\Delta\delta = \delta_{X} - \delta_{\mu}$) of about 1.6. Furthermore, from the start to the peak times of the HXR bursts, the time profiles of the HXR spectral index $\delta_{X}$ evolve synchronously with those of the microwave spectral index $\delta_{\mu}$, keeping the constant gap. We also examined the spatially resolved distribution of the microwave spectral index by using NoRH data. The microwave spectral index $\delta_{\mu}$ tends to be larger, which means a softer spectrum, at HXR footpoint sources with stronger magnetic field than that at the loop tops. These results suggest that the electron spectra are bent at around several hundreds of keV, and become harder at the higher energy range that contributes the microwave gyrosynchrotron emission.

Using two-dimensional magnetohydrodynamics calculations, we investigate a twist gen- eration mechanism on a magnetic flux tube at the base of the solar convection zone based on the idea of Choudhuri, 2003, Sol. Phys., 215, 31 in which a toroidal mag- netic field is wrapped by a surrounding mean poloidal field. During generation of the twist, the flux tube follows four phases. (1) It quickly splits into two parts with vortex motions rolling up the poloidal magnetic field. (2) Owing to the physical mechanism similar to that of the magneto-rotational instability, the rolled-up poloidal field is bent and amplified. (3) The magnetic tension of the disturbed poloidal magnetic field re- duces the vorticity, and the lifting force caused by vortical motion decreases. (4) The flux tube gets twisted and begins to rise again without splitting. Investigation of these processes is significant because it shows that a flux tube without any initial twist can rise to the surface in relatively weak poloidal fields.

We investigate an intensification mechanism for the magnetic field near the base of the solar convection zone that does not rely on differential rotation. Such mechanism in addition to differential rotation has been suggested by studies of flux emergence, which typically require field strength in excess of those provided by differential rotation alone. We study here a process in which potential energy of the superadiabatically stratified convection zone is converted into magnetic energy. This mechanism, know as explosion of magnetic flux tubes, has been previously studied in the thin flux tube approximation as well as two-dimensional MHD simulations, we expand the investigation to three-dimensional MHD simulations. Our main result is that enough intensification can be achieved in a three-dimensional magnetic flux sheet as long as the spatial scale of the imposed perturbation normal to the magnetic field is sufficiently large. When this spatial scale is small, the flux sheet tends to rise toward the surface, resulting in a significant decrease of the magnetic field amplification.

We have studied the relationship between the velocity and temperature of a solar EUV jet. The highly accelerated jet occurred in the active region NOAA 10960 on 2007 June 5. Multi-wavelength spectral observations with EIS/Hinode allow us to investigate Doppler velocities at the wide temperature range. We analyzed the three-dimensional angle of the jet from the stereoscopic analysis with STEREO. Using this angle and Doppler velocity, we derived the true velocity of the jet. As a result, we found that the cool jet observed with \ion{He}{2} 256 \AA $\log_{10}T_e[\rm{K}] = 4.9$ is accelerated to around $220 \rm{km/s}$ which is over the upper limit of the chromospheric evaporation. The velocities observed with the other lines are under the upper limit of the chromospheric evaporation while most of the velocities of hot lines are higher than that of cool lines. We interpret that the chromospheric evaporation and magnetic acceleration occur simultaneously. A morphological interpretation of this event based on the reconnection model is given by utilizing the multi-instrumental observations.

Frequencies of magnetic patch processes on supergranule boundary, namely flux emergence, splitting, merging, and cancellation, are investigated through an automatic detection. We use a set of line of sight magnetograms taken by the Solar Optical Telescope (SOT) on board Hinode satellite. We found 1636 positive patches and 1637 negative patches in the data set, whose time duration is 3.5 hours and field of view is 112" \times 112". Total numbers of magnetic processes are followed: 493 positive and 482 negative splittings, 536 positive and 535 negative mergings, 86 cancellations, and 3 emergences. Total numbers of emergence and cancellation are significantly smaller than those of splitting and merging. Further, frequency dependences of merging and splitting processes on flux content are investigated. Merging has a weak dependence on flux content only with a power- law index of 0.28. Timescale for splitting is found to be independent of parent flux content before splitting, which corresponds to \sim 33 minutes. It is also found that patches split into any flux contents with a same probability. This splitting has a power-law distribution of flux content with an index of -2 as a time independent solution. These results support that the frequency distribution of flux content in the analyzed flux range is rapidly maintained by merging and splitting, namely surface processes. We suggest a model for frequency distributions of cancellation and emergence based on this idea.

We investigate the value of horizontal turbulent diffusivity {\eta} by numerical calculation of thermal convection. In this study, we introduce a new method whereby the turbulent diffusivity is estimated by monitoring the time devel- opment of the passive scalar, which is initially distributed in a given Gaussian function with a spatial scale d0. Our conclusions are as follows: (1) Assuming the relation {\eta} = Lcvrms/3 where vrms is the RMS velocity, the characteristic length Lc is restricted by the shortest one among the pressure (density) scale height and the region depth. (2) The value of turbulent diffusivity becomes greater with the larger initial distribution scale d0. (3) The approximation of turbulent diffusion holds better when the ratio of the initial distribution scale d0 to the characteristic length Lc is larger.

It is widely accepted that solar active regions including sunspots are formed by the emerging magnetic flux from the deep convection zone. In previous numerical simulations, we found that the horizontal divergent flow (HDF) occurs before the flux emergence at the photospheric height. This Paper reports the HDF detection prior to the flux emergence of NOAA AR 11081, which is located away from the disk center. We use SDO/HMI data to study the temporal changes of the Doppler and magnetic patterns from those of the reference quiet Sun. As a result, the HDF appearance is found to come before the flux emergence by about 100 minutes. Also, the horizontal speed of the HDF during this time gap is estimated to be 0.6 to 1.5 km s^-1, up to 2.3 km s^-1. The HDF is caused by the plasma escaping horizontally from the rising magnetic flux. And the interval between the HDF and the flux emergence may reflect the latency during which the magnetic flux beneath the solar surface is waiting for the instability onset to the further emergence. Moreover, SMART Halpha images show that the chromospheric plages appear about 14 min later, located co-spatial with the photospheric pores. This indicates that the plages are caused by plasma flowing down along the magnetic fields that connect the pores at their footpoints. One importance of observing the HDF may be the possibility to predict the sunspot appearances that occur in several hours.

We investigate the frequency of magnetic activities, namely flux emergence, splitting, merging, and cancellation, through an automatic detection in order to understand the generation of the power-law distribution of magnetic flux reported by Parnell et al. (2009). The quiet Sun magnetograms observed in the Na I line by Hinode/ SOT is used in this study. Investigated patches range from 10^{17}Mx to 10^{19}Mx. Emergence and cancellation are much less frequent than merging and splitting. The time scale for splitting is found to be ~33 minutes and is independent of the flux contained in the splitting patch. Moreover magnetic patches split into any flux contents with an equal probability. It is shown that such a fragmentation process leads to a distribution with a power-law index -2. Merging has a very weak dependence on flux content only with a power-law index -0.33. These results suggest that 1) magnetic patches are fragmented by splitting, merging, and tiny cancellation; 2) flux is removed from the photosphere through tiny cancellations after these fragmentations.

We have performed a three-dimensional magnetohydrodynamic simulation to study the emergence of a twisted magnetic flux tube from -20,000 km of the solar convection zone to the corona through the photosphere and the chromosphere. The middle part of the initial tube is endowed with a density deficit to instigate a buoyant emergence. As the tube approaches the surface, it extends horizontally and makes a flat magnetic structure due to the photosphere ahead of the tube. Further emergence to the corona breaks out via the interchange-mode instability of the photospheric fields, and eventually several magnetic domes build up above the surface. What is new in this three-dimensional experiment is, multiple separation events of the vertical magnetic elements are observed in the photospheric magnetogram, and they reflect the interchange instability. Separated elements are found to gather at the edges of the active region. These gathered elements then show shearing motions. These characteristics are highly reminiscent of active region observations. On the basis of the simulation results above, we propose a theoretical picture of the flux emergence and the formation of an active region that explains the observational features, such as multiple separations of faculae and the shearing motion.

Context. The anelastic approximation is often adopted in numerical calculation with low Mach number, such as stellar internal convection. This approximation requires frequent global communication, because of an elliptic partial differential equation. Frequent global communication is negative factor for the parallel computing with a large number of CPUs. Aims. The main purpose of this paper is to test the validity of a method that artificially reduces the speed of sound for the compressible fluid equations in the context of stellar internal convection. The reduction of speed of sound allows for larger time steps in spite of low Mach number, while the numerical scheme remains fully explicit and the mathematical system is hyperbolic and thus does not require frequent global communication. Methods. Two and three dimensional compressible hydrodynamic equations are solved numerically. Some statistical quantities of solutions computed with different effective Mach numbers (due to reduction of speed of sound) are compared to test the validity of our approach. Results. Numerical simulations with artificially reduced speed of sound are a valid approach as long as the effective Mach number (based on the reduced speed of sound) remains less than 0.7.

We investigate differential rotation in rapidly rotating solar-type stars by means of an axisymmetric mean field model that was previously applied to the sun. This allows us to calculate the latitudinal entropy gradient with a rea- sonable physical basis. Our conclusions are as follows: (1) Differential rotation approaches the Taylor-Proudman state when stellar rotation is faster than so- lar rotation. (2) Entropy gradient generated by the attached subadiabatic layer beneath the convection zone becomes relatively small with a large stellar angu- lar velocity. (3) Turbulent viscosity and turbulent angular momentum transport determine the spatial difference of angular velocity $\Delta \Omega$. (4) The results of our mean field model can explain observations of stellar differential rotation.

We present the new results of the two-dimensional numerical experiments on the cross-sectional evolution of a twisted magnetic flux tube rising from the deeper solar convection zone (-20,000 km) to the corona through the surface. The initial depth is ten times deeper than most of previous calculations focusing on the flux emergence from the uppermost convection zone. We find that the evolution is illustrated by the two-step process described below: the initial tube rises due to its buoyancy, subject to aerodynamic drag due to the external flow. Because of the azimuthal component of the magnetic field, the tube maintains its coherency and does not deform to become a vortex roll pair. When the flux tube approaches the photosphere and expands sufficiently, the plasma on the rising tube accumulates to suppress the tube's emergence. Therefore, the flux decelerates and extends horizontally beneath the surface. This new finding owes to our large scale simulation calculating simultaneously the dynamics within the interior as well as above the surface. As the magnetic pressure gradient increases around the surface, magnetic buoyancy instability is triggered locally and, as a result, the flux rises further into the solar corona. We also find that the deceleration occurs at a higher altitude than in our previous experiment using magnetic flux sheets (Toriumi and Yokoyama). By conducting parametric studies, we investigate the conditions for the two-step emergence of the rising flux tube: field strength > 1.5x10^4 G and the twist > 5.0x10^-4 km^-1 at -20,000 km depth.

We present a series of numerical experiments that model the evolution of magnetic flux tubes with a different amount of initial twist. As a result of calculations, tightly twisted tubes reveal a rapid two-step emergence to the atmosphere with a slight slowdown at the surface, while weakly twisted tubes show a slow two-step emergence waiting longer the secondary instability to be triggered. This picture of the two-step emergence is highly consistent with recent observations. These tubes show multiple magnetic domes above the surface, indicating that the secondary emergence is caused by interchange mode of magnetic buoyancy instability. As for the weakest twist case, the tube exhibits an elongated photospheric structure and never rises into the corona. The formation of the photospheric structure is due to inward magnetic tension force of the azimuthal field component of the rising flux tube (i.e., tube's twist). When the twist is weak, azimuthal field cannot hold the tube's coherency, and the tube extends laterally at the subadiabatic surface. In addition, we newly find that the total magnetic energy measured above the surface depends on the initial twist. Strong twist tubes follow the initial relation between the twist and the magnetic energy, while weak twist tubes deviates from this relation, because these tubes store their magnetic energy in the photospheric structures.

In order to better understand the possibility of coronal heating by MHD waves, we analyze Fe xii 195.12{\AA} data observed with EUV Imaging Spectrometer (EIS) onboard Hinode. We performed a Fourier analysis of EUV intensity and Doppler velocity time series data in the active region corona. Notable intensity and Doppler velocity oscillations were found for two moss regions out of the five studied, while only small oscillations were found for five apexes of loops. The amplitudes of the oscillations were 0.4 - 5.7% for intensity and 0.2 - 1.2 kms-1 for Doppler velocity. In addition, oscillations of only Doppler velocity were seen relatively less often in the data. We compared the amplitudes of intensity and those of Doppler velocity in order to identify MHD wave modes, and calculated the phase delays between Fourier components of intensity and those of Doppler velocity. The results are interpreted in terms of MHD waves as follows: (1) few kink modes or torsional Alfv'en mode waves were seen in both moss regions and the apexes of loops; (2) upwardly propagating and standing slowmode waves were found inmoss regions; and (3) consistent with previous studies, estimated values of energy flux of the waves were several orders of magnitude lower than that required for heating active regions.

In this paper, results of 2.5-dimensional magnetohydrodynamical simulations are reported for the magnetic reconnection of non-perfectly antiparallel magnetic fields. The magnetic field has a component perpendicular to the computational plane, that is, guide field. The angle theta between magnetic field lines in two half regions is a key parameter in our simulations whereas the initial distribution of the plasma is assumed to be simple; density and pressure are uniform except for the current sheet region. Alfven waves are generated at the reconnection point and propagate along the reconnected field line. The energy fluxes of the Alfven waves and magneto-acoustic waves (slow mode and fast mode) generated by the magnetic reconnection are measured. Each flux shows the similar time evolution independent of theta. The percentage of the energies (time integral of energy fluxes) carried by the Alfven waves and magneto-acoustic waves to the released magnetic energy are calculated. The Alfven waves carry 38.9%, 36.0%, and 29.5% of the released magnetic energy at the maximum (theta=80^\circ) in the case of beta=0.1, 1, and 20 respectively, where beta is the plasma beta (the ratio of gas pressure to magnetic pressure). The magneto-acoustic waves carry 16.2% (theta=70^\circ), 25.9% (theta=60^\circ), and 75.0% (theta=180^\circ) of the energy at the maximum. Implications of these results for solar coronal heating and acceleration of high-speed solar wind are discussed.

We investigated the dependence of the solar magnetic parity between the hemispheres on two important parameters, the turbulent diffusivity and the meridional flow, by means of axisymmetric kinematic dynamo simulations based on the flux-transport dynamo model. It is known that the coupling of the magnetic field between hemispheres due to turbulent diffusivity is an important factor for the solar parity issue, but the detailed criterion for the generation of the dipole field has not been investigated. Our conclusions are as follows. (1) The stronger diffusivity near the surface is more likely to cause the magnetic field to be a dipole. (2) The thinner layer of the strong diffusivity near the surface is also more apt to generate a dipolar magnetic field. (3) The faster meridional flow is more prone to cause the magnetic field to be a quadrupole, i.e., symmetric about the equator. These results show that turbulent diffusivity and meridional flow are crucial for the configuration of the solar global magnetic field.

We perform two-dimensional MHD simulations on the solar flux emergence. We set the initial magnetic flux sheet at z=-20,000 km in the convection zone. The flux sheet rises through the convective layer due to the Parker instability, however, decelerates beneath the photosphere because the plasma on the flux sheet piles up owing to the convectively stable photosphere above. Meanwhile, the flux sheet becomes locally unstable to the Parker instability within the photosphere, and the further evolution to the corona occurs (two-step emergence model). We carry out a parameter survey to investigate the condition for this two-step model. We find that magnetic fluxes which form active regions are likely to have undergone the two-step emergence. The condition for the two-step emergence is 10^21 - 10^22 Mx with 10^4 G at z=-20,000 km in the convection zone.

A cancellation is thought to be a basic process of the photospheric magnetic field and plays an important role in magnetic flux budget and in various solar activities. There are two major theoretical scenarios for this phenomena, i.e. the "U-loop emergence" and the Omega-loop submergence models. It is important to clarify which is the dominant process during the cancellation for the estimation of the solar magnetic flux transport through the surface. We study the vector magnetic field and velocity structures around a quiet Sun cancellation by using the Solar Optical Telescope on board Hinode satellite. Transverse magnetic field connecting the canceling magnetic features and strong long-lasting Doppler red-shift signal are found. The transverse field is observed in the first spectropolarimetric observation after the occurrence of the cancellation while the red-shift clearly delayed to the cancellation by 20 minutes. These results indicate that the observed cancellation is an Omega-loop submergence.

We numerically investigate the long-term dynamics of the Saturn's small satellites Methone (S/2004 S1), Anthe (S/2007 S4) and Pallene (S/2004 S2). In our numerical integrations, these satellites are disturbed by non-spherical shape of Saturn and the six nearest regular satellites. The stability of the small bodies is studied here by analyzing long-term evolution of their orbital elements. We show that long-term evolution of Pallene is dictated by a quasi secular resonance involving the ascending nodes ($\Omega$) and longitudes of pericentric distances ($\varpi$) of Mimas (subscript 1) and Pallene (subscript 2), which critical argument is $\varpi_2-\varpi_1-\Omega_1+\Omega_2$. Long-term orbital evolution of Methone and Anthe are probably chaotic since: i) their orbits randomly cross the orbit of Mimas in time scales of thousands years); ii) numerical simulations involving both small satellites are strongly affected by small changes in the initial conditions.

The formation and evolution process and magnetic configuration of solar prominences remain unclear. In order to study the formation process of prominences, we examine continuous observations of a prominence in NOAA AR 10953 with the Solar Optical Telescope on the \emph{Hinode} satellite. As reported in our previous Letter, we find a signature suggesting that a helical flux rope emerges from below the photosphere under a pre-existing prominence. Here we investigate more detailed properties and photospheric indications of the emerging helical flux rope, and discuss their relationship to the formation of the prominence. Our main conclusions are: (1) A dark region with absence of strong vertical magnetic fields broadens and then narrows in Ca \textsc{ii} H-line filtergrams. This phenomenon is consistent with the emergence of the helical flux rope as photospheric counterparts. The size of the flux rope is roughly 30,000 km long and 10,000 km wide. The width is larger than that of the prominence. (2) No shear motion or converging flows are detected, but we find diverging flows such as mesogranules along the polarity inversion line. The presence of mesogranules may be related to the emergence of the helical flux rope. (3) The emerging helical flux rope reconnects with magnetic fields of the pre-existing prominence to stabilize the prominence for the next several days. We thus conjecture that prominence coronal magnetic fields emerge in the form of helical flux ropes that contribute to the formation and maintenance of the prominence.

We present a multi-wavelength observation of a solar flare occurring on 2006 December 13 with Hinode, RHESSI, and the Nobeyama Radio Observatory, to study the electron acceleration site and mechanism. The Solar Optical Telescope (SOT) on board Hinode observed elongated flare ribbons, and RHESSI observed double-footpoint hard X-ray (HXR) sources appearing in part of the ribbons. A photospheric vector magnetogram obtained from SOT reveals that the HXR sources are located at the region where horizontal magnetic fields change the direction. The region is interpreted as the footpoint of magnetic separatrix. Microwave images taken with the Nobeyama Radioheliograph show a loop structure connecting the HXR sources. The brighter parts of the microwave intensity are located between the top and footpoints of the loop. We consider these observations as an evidence of the electron acceleration near the magnetic separatrix and injection parallel to the field line.

We present a detailed examination on the coronal nonthermal emissions during the preflare phase of the X4.8 flare that occurred on 2002 July 23. The microwave (17 GHz and 34 GHz) data obtained with Nobeyama Radioheliograph, at Nobeyama Solar Radio Observatory and the hard X-ray (HXR) data taken with {\it Reuven Ramaty High Energy Solar Spectroscopic Imager} obviously showed nonthermal sources that are located above the flare loops during the preflare phase. We performed imaging spectroscopic analyses on the nonthermal emission sources both in microwaves and in HXRs, and confirmed that electrons are accelerated from several tens of keV to more than 1 MeV even in this phase. If we assume the thin-target model for the HXR emission source, the derived electron spectral indices ($\sim 4.7$) is the same value as that from microwaves ($\sim 4.7$) within the observational uncertainties, which implies that the distribution of the accelerated electrons follows a single power-law. The number density of the microwave-emitting electrons is, however, larger than that of the HXR-emitting electrons, unless we assume low ambient plasma density of about $1.0 \times 10^9$ cm$^{-3}$ for the HXR-emitting region. If we adopt the thick-target model for the HXR emission source, on the other hand, the electron spectral index ($\sim 6.7$) is much different, while the gap of the number density of the accelerated electrons is somewhat reduced.

In a previous work (Callegari and Yokoyama 2007, Celest. Mech. Dyn. Astr. vol. 98), the main features of the motion of the pair Enceladus-Dione were analyzed in the frozen regime, i.e., without considering the tidal evolution. Here, the results of a great deal of numerical simulations of a pair of satellites similar to Enceladus and Dione crossing the 2:1 mean-motion resonance are shown. The resonance crossing is modeled with a linear tidal theory, considering a two-degrees-of-freedom model written in the framework of the general three-body planar problem. The main regimes of motion of the system during the passage through resonance are studied in detail. We discuss our results comparing them with classical scenarios of tidal evolution of the system. We show new scenarios of evolution of the Enceladus-Dione system through resonance not shown in previous approaches of the problem.

We analytically and numerically study the motion of electrons along a magnetic loop, to compare with the observation of the propagating feature of the non-thermal microwave source in the 1999 August 28 solar flare reported by Yokoyama et al. (2002). We model the electron motion with the Fokker-Planck equation and calculate the spatial distribution of the gyrosynchrotron radiation. We find that the microwave propagating feature does not correspond to the motion of electrons with a specific initial pitch angle. This apparent propagating feature is a consequence of the motion of an ensemble of electrons with different initial pitch angles, which have different time and position to produce strong radiation in the loop. We conclude that the non-thermal electrons in the 1999 August 28 flare were isotropically accelerated and then are injected into the loop.

Continuous observations were obtained of active region 10953 with the Solar Optical Telescope (SOT) on board the \emph{Hinode} satellite during 2007 April 28 to May 9. A prominence was located over the polarity inversion line (PIL) in the south-east of the main sunspot. These observations provided us with a time series of vector magnetic fields on the photosphere under the prominence. We found four features: (1) The abutting opposite-polarity regions on the two sides along the PIL first grew laterally in size and then narrowed. (2) These abutting regions contained vertically-weak, but horizontally-strong magnetic fields. (3) The orientations of the horizontal magnetic fields along the PIL on the photosphere gradually changed with time from a normal-polarity configuration to a inverse-polarity one. (4) The horizontal-magnetic field region was blueshifted. These indicate that helical flux rope was emerging from below the photosphere into the corona along the PIL under the pre-existing prominence. We suggest that this supply of a helical magnetic flux into the corona is associated with evolution and maintenance of active-region prominences.

We study the non-thermal emissions in a solar flare occurring on 2003 May 29 by using RHESSI hard X-ray (HXR) and Nobeyama microwave observations. This flare shows several typical behaviors of the HXR and microwave emissions: time delay of microwave peaks relative to HXR peaks, loop-top microwave and footpoint HXR sources, and a harder electron energy distribution inferred from the microwave spectrum than from the HXR spectrum. In addition, we found that the time profile of the spectral index of the higher-energy ($\gsim 100$ keV) HXRs is similar to that of the microwaves, and is delayed from that of the lower-energy ($\lsim 100$ keV) HXRs. We interpret these observations in terms of an electron transport model called {\TPP}. We numerically solved the spatially-homogeneous {\FP} equation to determine electron evolution in energy and pitch-angle space. By comparing the behaviors of the HXR and microwave emissions predicted by the model with the observations, we discuss the pitch-angle distribution of the electrons injected into the flare site. We found that the observed spectral variations can qualitatively be explained if the injected electrons have a pitch-angle distribution concentrated perpendicular to the magnetic field lines rather than isotropic distribution.

Continuous observations of a flare productive active region 10930 were successfully carried out with the Solar Optical Telescope onboard the Hinode spacecraft during 2007 December 6 to 19. We focus on the evolution of photospheric magnetic fields in this active region, and magnetic field properties at the site of the X3.4 class flare, using a time series of vector field maps with high spatial resolution. The X3.4 class flare occurred on 2006 December 13 at the apparent collision site between the large, opposite polarity umbrae. Elongated magnetic structures with alternatingly positive and negative polarities resulting from flux emergence appeared one day before the flare in the collision site penumbra. Subsequently, the polarity inversion line at the collision site became very complicated. The number of bright loops in Ca II H increased during the formation of these elongated magnetic structures. The flare ribbons and bright loops evolved along the polarity inversion line and one footpoint of the bright loop was located in a region having a large departure of field azimuth angle with respect to its surroundings. The SOT observations with high spatial resolution and high polarization precision reveal temporal change in fine structure of magnetic fields at the flare site: some parts of the complicated polarity inversion line then disappeared, and in those regions the azimuth angle of photospheric magnetic field changed by about 90 degrees, becoming more spatially uniform within the collision site.

The Solar Optical Telescope (SOT) aboard HINODE successfully and continuously observed a formation process of a light bridge in a matured sunspot of the NOAA active region 10923 for several days with high spatial resolution. During its formation, many umbral dots were observed emerging from the leading edges of penumbral filaments, and intruding into the umbra rapidly. The precursor of the light bridge formation was also identified as the relatively slow inward motion of the umbral dots which emerged not near the penumbra, but inside the umbra. The spectro-polarimeter on SOT provided physical conditions in the photosphere around the umbral dots and the light bridges. We found the light bridges and the umbral dots had significantly weaker magnetic fields associated with upflows relative to the core of the umbra, which implies that there was hot gas with weak field strength penetrating from subphotosphere to near the visible surface inside those structures. There needs to be a mechanism to drive the inward motion of the hot gas along the light bridges. We suggest that the emergence and the inward motion are triggered by a buoyant penumbral flux tube as well as the subphotospheric flow crossing the sunspot.

On 2005 September 13 a filament eruption accompanied by a halo CME occurred in the most flare-productive active region NOAA 10808 in Solar Cycle 23. Using multi-wavelength observations before the filament eruption on Sep. 13th, we investigate the processes leading to the catastrophic eruption. We find that the filament slowly ascended at a speed of 0.1km/s over two days before the eruption. During slow ascending, many small flares were observed close to the footpoints of the filament, where new magnetic elements were emerging. On the basis of the observational facts we discuss the triggering mechanism leading to the filament eruption. We suggest the process toward the eruption as follows: First, a series of small flares played a role in changing the topology of the loops overlying the filament. Second, the small flares gradually changed the equilibrium state of the filament and caused the filament to ascend slowly over two days. Finally, a C2.9 flare that occurred when the filament was close to the critical point for loss of equilibrium directly led to the catastrophic filament eruption right after itself.

We consider MHD waves as a heating source of cool cores of galaxy clusters. In particular, we focus on transverse waves (Alfven waves), because they can propagate a longer distance than longitudinal waves (sound waves). Using MHD simulations, we found that the transverse waves can stably heat a cool core if the wave period is large enough (>~ 10^8 yr). Moreover, the longitudinal waves that are created as a by-product of the nonlinear evolution of the transverse waves could be observed as the 'ripples' found in cool cores.

It has been recognized that magnetic reconnection process is of great importance in high-energy astrophysics. We develop a new two-dimensional relativistic resistive magnetohydrodynamic (R$^2$MHD) code, and carry out numerical simulations of magnetic reconnection. We find that outflow velocity reaches Alfv\'{e}n velocity in the inflow region, and that higher Alfv\'{e}n velocity provides higher reconnection rate. We also find Lorentz contraction plays an important role in enhancement of reconnection rate.

We report a statistical study of flares observed with the Soft X-ray Telescope (SXT) onboard Yohkoh in the year of 2000. We measure physical parameters of 77 flares, such as the temporal scale, the size, and the magnetic flux density and find that the sizes of flares tend to be distributed more broadly as the GOES class becomes weaker and that there is a lower limit of magnetic flux density that depends on the GOES class. We also examine the relationship between these parameters and find weak correlation between temporal and spatial scales of flares. We estimate reconnection inflow velocity, coronal Alfven velocity, and reconnection rate using above observed values. The inflow velocities are distributed from a few km/s to several tens km/s and the Alfven velocities in the corona are in the range from 10^3 to 10^4 km/s. Hence the reconnection rate is 10^-3 - 10^-2. We find that the reconnection rate in a flare tends to decrease as the GOES class of the flare increases. This value is within one order of magnitude from the theoretical maximum value predicted by the Petschek model, although the dependence of the reconnection rate on the magnetic Reynolds number tends to be stronger than that in the Petschek model.

On December 27, 2004, plasma particle detectors on the GEOTAIL spacecraft detected an extremely strong signal of hard X-ray photons from the giant flare of SGR1806-20, a magnetar candidate. While practically all gamma-ray detectors on any satellites were saturated during the first ~500 ms interval after the onset, one of the particle detectors on GEOTAIL was not saturated and provided unique measurements of the hard X-ray intensity and the profile for the first 600 ms interval with 5.48 ms time resolution. After ~50 ms from the initial rapid onset, the peak photon flux (integrated above ~50 keV) reached the order of 10^7 photons sec^{-1} cm^{-2}. Assuming a blackbody spectrum with kT=175 keV, we estimate the peak energy flux to be 21 erg sec^{-1} cm^{-2} and the fluence (for 0-600 ms) to be 2.4 erg cm^{-2}. The implied energy release comparable to the magnetic energy stored in a magnetar (~10^{47} erg) suggests an extremely efficient energy release mechanism.

We present the relation of the spatial distribution of H-alpha kernels with the distribution of hard X-ray (HXR) sources seen during the 2001 April 10 solar flare. This flare was observed in H-alpha with the (Sartorius) telescope at Kwasan Observatory, Kyoto University, and in hard X-rays (HXRs) with the Hard X-ray Telescope (HXT) onboard Yohkoh. We compared the spatial distribution of the HXR sources with that of the H-alpha kernels. While many H-alpha kernels are found to brighten successively during the evolution of the flare ribbons, only a few radiation sources are seen in the HXR images. We measured the photospheric magnetic field strengths at each radiation source in the H-alpha images, and found that the H-alpha kernels accompanied by HXR radiation have magnetic strengths about 3 times larger than those without HXR radiation. We also estimated the energy release rates based on the magnetic reconnection model. The release rates at the H-alpha kernels with accompanying HXR sources are 16-27 times larger than those without HXR sources. These values are sufficiently larger than the dynamic range of HXT, which is about 10, so that the difference between the spatial distributions of the H$\alpha$ kernels and the HXR sources can be explained.

We propose the Galactic flare model for the origin of the X-ray gas in the Galactic halo. For this purpose, we examine the magnetic reconnection triggered by Parker instability (magnetic buoyancy instability), by performing the two-dimensional resistive numerical magnetohydrodynamic simulations. As a result of numerical simulations, the system evolves as following phases: Parker instability occurs in the Galactic disk. In the nonlinear phase of Parker instability, the magnetic loop inflates from the Galactic disk into the Galactic halo, and collides with the anti-parallel magnetic field, so that the current sheets are created in the Galactic halo. The tearing instability occurs, and creates the plasmoids (magnetic islands). Just after the plasmoid ejection, further current-sheet thinning occurs in the sheet, and the anomalous resistivity sets in. Petschek reconnection starts, and heats the gas quickly in the Galactic halo. It also creates the slow and fast shock regions in the Galactic halo. The magnetic field ($B\sim 3 \mu$G), for example, can heat the gas ($n\sim 10^{-3}$ cm$^{-3}$) to temperature of $\sim 10^6$ K via the reconnection in the Galactic halo. The gas is accelerated to Alfv\'en velocity ($\sim 300$ km s$^{-1}$). Such high velocity jets are the evidence of the Galactic flare model we present in this paper, if the Doppler shift of the bipolar jet is detected in the Galactic halo. Full size figures are available at http://www.kwasan.kyoto-u.ac.jp/~tanuma/study/ApJ2002/ApJ2002.html

In our previous paper, we have presented a theory to explain the observed universal correlation between the emission measure ($EM=n^2 V$) and temperature (T) for solar/stellar flares on the basis of the magnetic reconnection model with heat conduction and chromospheric evaporation. Here n is the electron density and V is the volume. By extending our theory to general situations, we examined the EM-T diagram in detail, and found the following properties: 1) The universal correlation sequence (``main sequence flares'') with $EM \propto T^{17/2}$ corresponds to the case of constant heating flux or equivalently the case of constant magnetic field strength in the reconnection model. 2) The EM-T diagram has a forbidden region, where gas pressure of flares exceeds magnetic pressure. 3) There is a coronal branch with $EM \propto T^{15/2}$ for $T< 10^7$ K, and $EM \propto T^{13/2}$ for $T> 10^7$ K. This branch is situated left side of the main sequence flares in the EM-T diagram. 4) There is another forbidden region determined by the length of flare loop; a lower limit of flare loop is $10^7$ cm. Small flares near this limit correspond to nanoflares observed by SOHO/EIT. 5) We can plot flare evolution track on the EM-T diagram. A flare evolves from the coronal branch to main sequence flares, then returns to the coronal branch eventually. These properties of the EM-T diagram are similar to those of the HR diagram for stars, and thus we propose that the EM-T diagram is quite useful to estimate the physical quantities (loop length, heating flux, magnetic field strength, total energy and so on) of flares and corona when there is no spatially resolved imaging observations.

We present an examination of the multi-wavelength observation of a C7.9 flare which occurred on 1998 November 10. This is the first time of imaging observation of the quasi-periodic pulsations (QPPs). Four bursts were observed with the hard X-ray telescope aboard Yohkoh and the Nobeyama Radioheliograph during the impulsive phase of the flare. In the second burst, the hard X-ray and microwave time profiles clearly showed a QPP. We estimated the Alfven transit time along the flare loop using the images of the soft X-ray telescope aboard it Yohkoh and the photospheric magnetgrams, and found that the transit time was almost equal to the period of the QPP. We therefore suggest, based on a shock acceleration model, that variations of macroscopic magnetic structures, such as oscillations of coronal loops, affect the efficiency of particle injection/acceleration.

We examine the magnetic reconnection triggered by a supernova (or a point explosion) in interstellar medium, by performing two-dimensional resistive magnetohydrodynamic (MHD) numerical simulations with high spatial resolution. We found that the magnetic reconnection starts long after a supernova shock (fast-mode MHD shock) passes a current sheet. The current sheet evolves as follows: (i) Tearing-mode instability is excited by the supernova shock, and the current sheet becomes thin in its nonlinear stage. (ii) The current-sheet thinning is saturated when the current-sheet thickness becomes comparable to that of Sweet-Parker current sheet. After that, Sweet-Parker type reconnection starts, and the current-sheet length increases. (iii) ``Secondary tearing-mode instability'' occurs in the thin Sweet-Parker current sheet. (iv) As a result, further current-sheet thinning occurs and anomalous resistivity sets in, because gas density decreases in the current sheet. Petschek type reconnection starts and heats interstellar gas. Magnetic energy is released quickly while magnetic islands are moving in the current sheet during Petschek type reconnection. The released magnetic energy is determined by the interstellar magnetic field strength, not energy of initial explosion nor distance to explosion. We suggest that magnetic reconnection is a possible mechanism to generate X-ray gas in Galaxy.