Source author record

Takaaki Yokoyama

Takaaki Yokoyama appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

25works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

25 published item(s)

preprint2022arXiv

Mean-field analysis on large-scale magnetic fields at high Reynolds numbers

Solar magnetic fields comprise an 11-year activity cycle, represented by the number of sunspots. The maintenance of such a solar magnetic field can be attributed to fluid motion in the convection zone, i.e. a dynamo. This study conducts the mean-field analyses of the global solar dynamo simulation presented by Hotta et al. (2016). Although the study succeeds in producing coherent large-scale magnetic fields at high Reynolds numbers, the detailed physics of the maintenance of this field have not been fully understood. This study extracts the alpha-tensor and the turbulent magnetic diffusivity tensor through mean-field analyses. The turbulent magnetic diffusivity exhibits a significant decrease towards high Reynolds numbers. The decrease in the turbulent magnetic diffusivity suppresses the energy conversion of large-scale field to small-scale field. This implies that the decrease in the turbulent magnetic diffusivity contributes to the maintenance of a large-scale magnetic field at high Reynolds numbers. A significant downward turbulent pumping is observed; it is enhanced in the weak phase of the large-scale field. This study proposes a cyclic reversal process of a large-scale field which is dominantly driven by the alpha-effect and is possibly triggered by downward pumping.

preprint2020arXiv

Inverse FIP Effects in Giant Solar Flares Found from Earth X-Ray Albedo with Suzaku/XIS

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$α$, Si Ly$α$, S He$α$, S Ly$α$, Ar He$α$, and Ca He$α$ originating from solar flares as well as fluorescent Ar K$α$ and Ar K$β$ 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.

preprint2020arXiv

Simulation of Alfven wave propagation in magnetic chromosphere with radiative loss: effects of non-linear mode coupling on chromospheric heating

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.

preprint2019arXiv

Estimating the temperature and density of a spicule from 100 GHz data obtained with ALMA

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Å 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.

preprint2016arXiv

Modelling observed decay-less oscillations as resonantly enhanced Kelvin-Helmholtz vortices from transverse MHD waves and their seismological application

In the highly structured solar corona, resonant absorption is an unavoidable mechanism of energy transfer from global transverse MHD waves to local azimuthal Alfvén 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én 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.

preprint2016arXiv

Nonlinear reflection process of linearly-polarized, broadband Alfvén waves in the fast solar wind

Using one-dimensional numerical simulations, we study the elementary process of Alfvén wave reflection in a uniform medium, including nonlinear effects. In the linear regime, Alfvé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én wave. In this study we consider a linearly polarized Alfvén 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én wave to the backscattered one. Such nonlinear reflection explains the observed increasing energy ratio of the sunward to the anti-sunward Alfvénic fluctuations in the solar wind with distance against the dynamical alignment effect.

preprint2015arXiv

Electron density of active region outflows measured by the EUV Imaging Spectrometer onboard Hinode

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.

preprint2015arXiv

Numerical Study on In-Situ Prominence Formation by Radiative Condensation in the Solar Corona

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.

preprint2015arXiv

Resonant Absorption of Transverse Oscillations and Associated Heating in a Solar Prominence. I- Observational aspects

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.

preprint2015arXiv

Resonant Absorption of Transverse Oscillations and Associated Heating in a Solar Prominence. II- Numerical aspects

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.

preprint2015arXiv

Three-Dimensional MHD Magnetic Reconnection Simulations with Finite Guide Field: Proposal of the Shock-Evoking Positive-Feedback Model

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.

preprint2014arXiv

Statistical Analysis of the Horizontal Divergent Flow in Emerging Solar Active Regions

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.

preprint2013arXiv

Probing the Shallow Convection Zone: Rising Motion of Subsurface Magnetic Fields in the Solar Active Region

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.

preprint2013arXiv

Three-dimensional magnetohydrodynamic simulation of the solar magnetic flux emergence: Parametric study on the horizontal divergent flow

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.

preprint2011arXiv

Dependence of the Magnetic Energy of Solar Active Regions on the Twist Intensity of the Initial Flux Tubes

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.

preprint2011arXiv

Modeling of differential rotation in rapidly rotating solar-type stars

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 $ΔΩ$. (4) The results of our mean field model can explain observations of stellar differential rotation.

preprint2011arXiv

Numerical Experiments on the Two-step Emergence of Twisted Magnetic Flux Tubes in the Sun

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.

preprint2010arXiv

Generation of Alfven Waves by Magnetic Reconnection

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.

preprint2010arXiv

Mode identification of MHD waves in an active region observed with Hinode/EIS

In order to better understand the possibility of coronal heating by MHD waves, we analyze Fe xii 195.12Å 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.

preprint2010arXiv

Two-step Emergence of the Magnetic Flux Sheet from the Solar Convection Zone

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.

preprint2009arXiv

Imaging Spectroscopy on Preflare Coronal Nonthermal Sources Associated with the 2002 July 23 Flare

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.

preprint2009arXiv

Prominence Formation Associated with an Emerging Helical Flux Rope

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.

preprint2007arXiv

Hinode Observations of Vector Magnetic Field Change Associated with a Flare on 2006 December 13

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.

preprint2005arXiv

GEOTAIL observation of the SGR1806-20 Giant Flare: The first 600 ms

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.

preprint2002arXiv

A HR-like Diagram for Solar/Stellar Flares and Corona -- Emission Measure vs Temperature Diagram

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.