Source author record

Yuta Notsu

Yuta Notsu 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

18works
7topics
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

18 published item(s)

preprint2023arXiv

Multiwavelength observation of an active M-dwarf star EV Lac and its stellar flare accompanied by a delayed prominence eruption

We conducted 4-night multiwavelength observations of an active M-dwarf star EV Lac on 2022 October 24$-$27 with simultaneous coverage of soft X-rays (NICER; 0.2$-$12 $\mathrm{keV}$, Swift XRT; 0.2$-$10 $\mathrm{keV}$), near-ultraviolet (Swift UVOT/UVW2; 1600$-$3500 Å), optical photometry (TESS; 6000$-$10000 Å), and optical spectroscopy (Nayuta/MALLS; 6350$-$6800 Å). During the campaign, we detected a flare starting at 12:28 UTC on October 25 with its white-light bolometric energy of $3.4 \times 10^{32}$ erg. At about 1 hour after this flare peak, our $\mathrm{Hα}$ spectrum showed a blue-shifted excess component at its corresponding velocity of $\sim 100 \: \mathrm{km \: s^{-1}}$. This may indicate that the prominence erupted with a 1-hour delay of the flare peak. Furthermore, the simultaneous 20-second cadence near-ultraviolet and white-light curves show gradual and rapid brightening behaviors during the rising phase at this flare. The ratio of flux in NUV to white light at the gradual brightening was $\sim 0.49$, which may suggest that the temperature of the blackbody is low ($< 9000 \: \mathrm{K}$) or the maximum energy flux of a nonthermal electron beam is less than $5\times10^{11} \: \mathrm{erg \: cm^{-2} \: s^{-1}}$. Our simultaneous observations of NUV and white-light flare raise the issue of a simple estimation of UV flux from optical continuum data by using a blackbody model.

preprint2022arXiv

Discovery of a Long-Duration Superflare on a Young Solar-Type Star EK Draconis with Nearly Similar Time Evolution for H$α$ and White-Light Emissions

Young solar-type stars are known to show frequent "superflares", which may severely influence the habitable worlds on young planets via intense radiations and coronal mass ejections. Here we report an optical spectroscopic and photometric observation of a long-duration superflare on the young solar-type star EK Draconis (50-120 Myr age) with the Seimei telescope and $Transiting$ $Exoplanet$ $Survey$ $Satellite$ ($TESS$). The flare energy 2.6$\times$10$^{34}$ erg and white-light flare duration 2.2 hr are much larger than those of the largest solar flares, and this is the largest superflare on a solar-type star ever detected by optical spectroscopy. The H$α$ emission profile shows no significant line asymmetry, meaning no signature of a filament eruption, unlike the only previous detection of a superflare on this star (Namekata et al. 2021, $Nat.Astron$). Also, it did not show significant line broadening, indicating that the non-thermal heating at the flare footpoints are not essential or that the footpoints are behind the limb. The time evolution and duration of the H$α$ flare are surprisingly almost the same as those of the white-light flare, which is different from general M-dwarf (super-)flares and solar flares. This unexpected time evolution may suggest that different radiation mechanisms than general solar flares are predominant, as follows: (1) radiation from (off-limb) flare loops, and (2) re-radiation via radiative backwarming, in both of which the cooling timescales of flare loops could determine the timescales of H$α$ and white light.

preprint2020arXiv

Cost estimation for alternative aviation plans against potential radiation exposure associated with solar proton events for the airline industry

We present a systematic approach to effectively evaluate potential risk cost caused by exposure to solar proton events (SPEs) from solar flares for the airline industry. We also evaluate associated health risks from radiation, to provide relevant alternative ways to minimize economic loss and opportunity. The estimated radiation dose induced by each SPE for the passengers of each flight is calculated using ExoKyoto and PHITS. We determine a few scenarios for the estimated dose limit at 1 and 20mSv, corresponding to the effective dose limit for the general public and occupational exposure, respectively, as well as a higher dose induced an extreme superflare. We set a hypothetical airline shutdown scenario at 1mSv for a single flight per passenger, due to legal restrictions under the potential radiation dose. In such a scenario, we calculate the potential loss in direct and opportunity cost under the cancelation of the flight. At the same time, we considered that, even under such a scenario, if the airplane flies at a slightly lower altitude (from 12 to 9.5km: atmospheric depth from 234 to 365g/cm$^{2}$), the total loss becomes much smaller than flight cancelation, and the estimated total dose goes down from 1.2 to 0.45mSv, which is below the effective dose limit for the general public. In case of flying at an even lower altitude (7km: atmospheric depth 484g/cm$^{2}$), the estimated total dose becomes much smaller, 0.12 mSv. If we assume the increase of fuel cost is proportional to the increase in atmospheric depth, the increase in cost becomes 1.56 and 2.07 for the case of flying at 9.5 km and at 7 km, respectively. Lower altitude flights provide more safety for the potential risk of radiation doses induced by severe SPEs. At the same time, since there is total loss caused by flight cancelation, we propose that considering lower flight altitude is the best protection against solar flares.

preprint2020arXiv

Optical and X-ray observations of stellar flares on an active M dwarf AD Leonis with Seimei Telescope, SCAT, NICER and OISTER

We report multi-wavelength monitoring observations of an M-dwarf flare star AD Leonis with Seimei Telescope (6150--7930 Å), SCAT (Spectroscopic Chuo-university Astronomical Telescope; 3700--7500 Å), NICER (Neutron Star Interior Composition Explorer; 0.2--12.0 keV), and collaborations of OISTER (Optical and Infrared Synergetic Telescopes for Education and Research) program. Twelve flares are detected in total which include ten H$α$, four X-ray, and four optical-continuum flares; one of them is a superflare with the total energy of $\sim$ 2.0$\times$10$^{33}$ erg. We found that (1) during the superflare, the H$α$ emission line full width at 1/8 maximum dramatically increases to 14 Å from 8 Å in the low-resolution spectra (R$\sim$ 2000) accompanied with the large white-light flares, (2) some weak H$α$/X-ray flares are not accompanied with white-light emissions, and (3) the non-flaring emissions show clear rotational modulations in X-ray and H$α$ intensity in the same phase. To understand these observational features, one-dimensional hydrodynamic flare simulations are performed by using the RADYN code. As a result of simulations, we found the simulated H$α$ line profiles with hard and high-energy non-thermal electron beams are consistent with that of the initial phase line profiles of the superflares, while those with more soft- and/or weak-energy beam are consistent with those in decay phases, indicating the changes in the energy fluxes injected to the lower atmosphere. Also, we found that the relation between optical continuum and H$α$ intensity is nonlinear, which can be one cause of the non-white-light flares. The flare energy budget exhibits diversity in the observations and models, and more observations of stellar flares are necessary for constraining the occurrence of various emission line phenomena in stellar flares.

preprint2020arXiv

Solar Analogs as a Tool to Understand the Sun

Solar analogs, broadly defined as stars similar to the Sun in mass or spectral type, provide a useful laboratory for exploring the range of Sun-like behaviors and exploring the physical mechanisms underlying some of the Sun's most elusive processes like coronal heating and the dynamo. We describe a series of heliophysics-motivated, but astrophysics-like studies of solar analogs. We argue for a range of stellar observations, including (a) the identification and fundamental parameter determination of new solar analogs, and (b) characterizing emergent properties like activity, magnetism, and granulation. These parameters should be considered in the framework of statistical studies of the dependences of these observables on fundamental stellar parameters like mass, metallicity, and rotation.

preprint2020arXiv

Statistical Properties of Superflares on Solar-type Stars: Results Using All of the Kepler Primary Mission Data

We report the latest statistical analyses of superflares on solar-type (G-type main-sequence; effective temperature is 5100 - 6000 K) stars using all of the $Kepler$ primary mission data, and $Gaia$-DR2 (Data Release 2) catalog. We updated the flare detection method from our previous studies by using high-pass filter to remove rotational variations caused by starspots. We also examined the sample biases on the frequency of superflares, taking into account gyrochronology and flare detection completeness. The sample size of solar-type stars and Sun-like stars (effective temperature is 5600 - 6000 K and rotation period is over 20 days in solar-type stars) are $\sim$4 and $\sim$12 times, respectively, compared with Notsu et al. (2019, ApJ, 876, 58). As a result, we found 2341 superflares on 265 solar-type stars, and 26 superflares on 15 Sun-like stars: the former increased from 527 to 2341 and the latter from 3 to 26 events compared with our previous study. This enabled us to have a more well-established view on the statistical properties of superflares. The observed upper limit of the flare energy decreases as the rotation period increases in solar-type stars. The frequency of superflares decreases as the stellar rotation period increases. The maximum energy we found on Sun-like stars is $4 \times 10^{34}$ erg. Our analysis of Sun-like stars suggest that the Sun can cause superflares with energies of $\sim 7 \times 10^{33}$ erg ($\sim$X700-class flares) and $\sim 1 \times 10^{34}$ erg ($\sim$X1000-class flares) once every $\sim$3,000 years and $\sim$6,000 years, respectively.

preprint2020arXiv

Temporal Evolution of Spatially-Resolved Individual Star Spots on a Planet-Hosting Solar-type Star: Kepler 17

Star spot evolution is visible evidence of the emergence/decay of the magnetic field on stellar surface, and it is therefore important for the understanding of the underlying stellar dynamo and consequential stellar flares. In this paper, we report the temporal evolution of individual star spot area on the hot-Jupiter-hosting active solar-type star Kepler 17 whose transits occur every 1.5 days. The spot longitude and area evolution are estimated (1) from the stellar rotational modulations of Kepler data and (2) from the brightness enhancements during the exoplanet transits caused by existence of large star spots. As a result of the comparison, number of spots, spot locations, and the temporal evolution derived from the rotational modulations is largely different from those of in-transit spots. We confirm that although only two light curve minima appear per rotation, there are clearly many spots present on the star. We find that the observed differential intensity changes are sometimes consistent with the spot pattern detected by transits, but they sometimes do not match with each other. Although the temporal evolution derived from the rotational modulation differs from those of in-transit spots to a certain degree, the emergence/decay rates of in-transit spots are within an order of magnitude of those derived for sunspots as well as our previous research based only on rotational modulations. This supports a hypothesis that the emergence/decay of sunspots and extremely-large star spots on solar-type stars occur through the same underlying processes.

preprint2020arXiv

The Extreme Space Weather Event in 1903 October/November: An Outburst from the Quiet Sun

While the Sun is generally more eruptive during its maximum and declining phases, observational evidence shows certain cases of powerful solar eruptions during the quiet phase of the solar activity. Occurring in the weak Solar Cycle 14 just after its minimum, the extreme space weather event in 1903 October -- November was one of these cases. Here, we reconstruct the time series of geomagnetic activity based on contemporary observational records. With the mid-latitude magnetograms, the 1903 magnetic storm is thought to be caused by a fast coronal mass ejection (~1500 km/s) and is regarded as an intense event with an estimated minimum Dst' of ~-513 nT The reconstructed time series has been compared with the equatorward extension of auroral oval (~44.1° in invariant latitude) and the time series of telegraphic disturbances. This case study shows that potential threats posed by extreme space weather events exist even during weak solar cycles or near their minima.

preprint2020arXiv

Time-resolved spectroscopy and photometry of an M dwarf flare star YZ Canis Minoris with OISTER and TESS: Blue asymmetry in H$α$ line during the non-white light flare

In this paper, we present the results from spectroscopic and photometric observations of the M-type flare star YZ CMi in the framework of the Optical and Infrared Synergetic Telescopes for Education and Research (OISTER) collaborations during the Transiting Exoplanet Survey Satellite (TESS) observation period. We detected 145 white-light flares from the TESS light curve and 4 H$α$ flares from the OISTER observations performed between 2019-01-16 and 2019-01-18. Among them, 3 H$α$ flares were associated with white-light flares. However, one of them did not show clear brightening in continuum; during this flare, the H$α$ line exhibited blue-asymmetry which has lasted for $\sim 60$ min. The line of sight velocity of the blue-shifted component is $-80$ - $-100$ km s$^{-1}$. This suggests that there can be upward flows of chromospheric cool plasma even without detectable red/NIR continuum brightening. By assuming that the blue-asymmetry in H$α$ line was caused by a prominence eruption on YZ CMi, we estimated the mass and kinetic energy of the upward-moving material to be $10^{16}$ - $10^{18}$ g and $10^{29.5}$ - $10^{31.5}$ erg, respectively. The estimated mass is comparable to expectations from the empirical relation between the flare X-ray energy and mass of upward-moving material for stellar flares and solar CMEs. In contrast, the estimated kinetic energy for the non-white-light flare on YZ CMi is roughly $2$ orders of magnitude smaller than that expected from the relation between flare X-ray energy and kinetic energy for solar CMEs. This could be understood by the difference in the velocity between CMEs and prominence eruptions.

preprint2016arXiv

Statistical properties of superflares on solar-type stars with Kepler data

Superflares are flares that release total energy 10$\sim$10$^{4}$ times greater than that of the biggest solar flares with energy of $\sim$10$^{32}$ erg. We searched superflares on solar-type stars (G-type main sequence stars) using the Kepler 30-min (long) and 1-min (short) cadence data. We found more than 1500 superflares on 279 stars from 30-min cadence data (Q0-6) and 187 superflares on 23 stars from 1-min cadence data (Q0-17). The bolometric energy of detected superflares ranges from the order of 10$^{32}$ erg to 10$^{36}$ erg. Using these data, we found that the occurrence frequency ($dN/dE$) of superflares is expressed as a power-law function of flare energy ($E$) with the index of -1.5 for $10^{33}<E<10^{36}$erg. Most of the superflare stars show quasi-periodic light variations with the amplitude of a few percent, which can be explained by the rotation of the star with large starspots. The bolometric energy released by flares is consistent with the magnetic energy stored around such large starspots. Furthermore, our analyses indicate that the occurrence frequency of superflares depends on the rotation period, and that the flare frequency increases as the rotation period decreases. However, the energy of the largest flares observed in a given period bin does not show any clear correlation with the rotation period. We also found that the duration of superflares increases with the flare energy as $E^{0.39+/-0.03}$. This can be explained if we assume the time-scale of flares is determined by the Alfv$\acute{\rm{e}}$n time.

preprint2015arXiv

High Dispersion Spectroscopy of Solar-type Superflare Stars. I. Temperature, Surface Gravity, Metallicity, and $v \sin i$

We conducted high dispersion spectroscopic observations of 50 superflare stars with Subaru/HDS, and measured the stellar parameters of them. These 50 targets were selected from the solar-type (G-type main sequence) superflare stars that we had discovered from the Kepler photometric data. As a result of these spectroscopic observations, we found that more than half (34 stars) of our 50 targets have no evidence of binary system. We then estimated effective temperature ($T_{\rm{eff}}$), surface gravity ($\log g$), metallicity ([Fe/H]), and projected rotational velocity ($v\sin i$) of these 34 superflare stars on the basis of our spectroscopic data. The accuracy of our estimations is higher than that of Kepler Input Catalog (KIC) values, and the differences between our values and KIC values ($(ΔT_{\rm{eff}})_{\rm{rms}} \sim 219$K, $(Δ\log g)_{\rm{rms}} \sim 0.37$ dex, and $(Δ\rm{[Fe/H]})_{\rm{rms}} \sim 0.46$ dex) are comparable to the large uncertainties and systematic differences of KIC values reported by the previous researches. We confirmed that the estimated $T_{\rm{eff}}$ and $\log g$ values of the 34 superflare stars are roughly in the range of solar-type stars. In particular, these parameters and the brightness variation period ($P_{0}$) of 9 stars are in the range of "Sun-like" stars ($5600\leq T_{\rm{eff}}\leq 6000$K, $\log g\geq$4.0, and $P_{0}>$10 days). Five of the 34 target stars are fast rotators ($v \sin i \geq 10$km s$^{-1}$), while 22 stars have relatively low $v \sin i$ values ($v \sin i<5$km s$^{-1}$). These results suggest that stars whose spectroscopic properties similar to the Sun can have superflares, and this supports the hypothesis that the Sun might cause a superflare.

preprint2015arXiv

High Dispersion Spectroscopy of Solar-type Superflare Stars. III. Lithium Abundances

We report on the abundance analysis of Li in solar-type (G-type main sequence) superflare stars which were found by the analysis of Kepler photometric data. Li is a key element to understand the evolution of the stellar convection zone which reflects the age of solar-type stars. We performed the high dispersion spectroscopy of solar-type superflare stars with Subaru/HDS, and confirmed that 34 stars show no evidence of binarity in our previous study. In this study, we derived the Li abundances of these 34 objects. We investigate correlations of Li abundance with stellar atmospheric parameters, rotational velocity, and superflare activities to understand the nature of superflare stars and the possibility of the nucleosynthesis of Li by superflares. We confirm the large dispersion in the Li abundance, and the correlation with stellar parameters is not seen. As compared with the Li abundance in Hyades cluster which is younger than the Sun, it is suggested that half of the observed stars are younger than Hyades cluster. The measured value of $v \sin i$ (projected rotational velocity) supports those objects are younger than the Sun. However, there are some objects which show the low Li abundance and slowly rotate on the basis of the estimated $v \sin i$ and $P$ (period of brightness variation). This result indicates that the superflare stars are not only young stars but also old stars like our Sun. In our observations, we could not find any evidence of Li productions by superflares. Further research on Li isotope abundances of superflare stars would clarify the Li production by stellar flares.

preprint2015arXiv

Statistical properties of superflares on solar-type stars based on 1-min cadence data

We searched for superflares on solar-type stars using Kepler data with 1 min sampling in order to detect superflares with short duration. We found 187 superflares on 23 solar-type stars whose bolometric energy ranges from the order of $10^{32}$ erg to $10^{36}$ erg. Some superflares show multiple peaks with the peak separation of the order of $100$-$1000$ seconds which is comparable to the periods of quasi-periodic pulsations in solar and stellar flares. Using these new data combined with the results from the data with 30 min sampling, we found the occurrence frequency ($dN/dE$) of superflares as a function of flare energy ($E$) shows the power-law distribution ($dN/dE \propto E^{-α}$) with $α\sim -1.5$ for $10^{33}<E<10^{36}$ erg which is consistent with the previous results. The average occurrence rate of superflares with the energy of $10^{33}$ erg which is equivalent to X100 solar flares is about once in 500-600 years. The upper limit of energy released by superflares is basically comparable to a fraction of the magnetic energy stored near starspots which is estimated from the photometry. We also found that the duration of superflares ($τ$) increases with the flare energy ($E$) as $τ\propto E^{0.39\pm 0.03}$. This can be explained if we assume the time-scale of flares is determined by the Alfv$\acute{\rm e}$n time.

preprint2014arXiv

High Dispersion Spectroscopy of Solar-type Superflare Stars. II. Stellar Rotation, Starspots, and Chromospheric Activities

We conducted high dispersion spectroscopic observations of 50 superflare stars with Subaru/HDS. These 50 stars were selected from the solar-type superflare stars that we had discovered from the Kepler data. More than half (34 stars) of these 50 target superflare stars show no evidence of binarity, and we estimated stellar parameters of these 34 stars in our previous study (Notsu et al. 2015, hereafter called Paper I). According to our previous studies using Kepler data, superflare stars show quasi-periodic brightness variations whose amplitude (0.1-10\%) is much larger than that of the solar brightness variations (0.01-0.1\%) caused by the existence of sunspots on the rotating solar surface. In this study, we investigated whether these quasi-periodic brightness variations of superflare stars are explained by the rotation of a star with fairly large starspots, by using stellar parameters derived in Paper I. First, we confirmed that the value of the projected rotational velocity $v \sin i$ is consistent with the rotational velocity estimated from the period of the brightness variation. Next, we measured the intensity of Ca II infrared triplet lines and H$α$ line, good indicators of the stellar chromospheric activity, and compared them with other stellar properties. The intensity of Ca II infrared triplet lines indicates that the mean magnetic field strength ($\langle fB\rangle$) of the target superflare stars can be higher than that of the Sun. A correlation between the amplitude of the brightness variation and the intensity of Ca II triplet line was found. All the targets expected to have large starspots because of their large amplitude of the brightness variation show high chromospheric activities compared to the Sun. These results support that the brightness variation of superflare stars is due to the rotation with large starspots.

preprint2014arXiv

Two Sun-like Superflare Stars Rotating as Slow as the Sun

We report on the results of high dispersion spectroscopy of two `superflare stars', KIC 9766237, and KIC 9944137 with Subaru/HDS. Superflare stars are G-type main sequence stars, but show gigantic flares compared to the Sun, which have been recently discovered in the data obtained with the Kepler spacecraft. Though most of these stars are thought to have a rotation period shorter than 10 days on the basis of photometric variabilities, the two targets of the present paper are estimated to have a rotation period of 21.8 d, and 25.3 d. Our spectroscopic results clarified that these stars have stellar parameters similar to those of the Sun in terms of the effective temperature, surface gravity, and metallicity. The projected rotational velocities derived by us are consistent with the photometric rotation period, indicating a fairy high inclination angle. The average strength of the magnetic field on the surface of these stars are estimated to be 1-20 G, by using the absorption line of Ca II 8542. We could not detect any hint of binarity in our spectra, although more data are needed to firmly rule out the presence of an unseen low-mass companion. These results claim that the spectroscopic properties of these superflare stars are very close to those of the Sun, and support the hypothesis that the Sun might cause a superflare.

preprint2013arXiv

High Dispersion Spectroscopy of the Superflare Star KIC6934317

We conducted the high-resolution spectroscopic observation with Subaru/HDS for a G-type star (KIC6934317). We selected this star from the data of the Kepler spacecraft. This star produces a lot of superflares, and the total energy of the largest superflare on this star is ~10^3 times larger (~2.2*10^35 erg) than that of the most energetic flare on the Sun (~10^32 erg). The core depth and emission flux of Ca II infrared triplet lines and Ha line show high chromospheric activity in this star, in spite of its low lithium abundance and the small amplitude of the rotational modulation. Using the empirical relations between emission flux of chromospheric lines and X-ray flux, this star is considered to show much higher coronal activity than that of the Sun. It probably has large starspots which can store a large amount of magnetic energy enough to give rise to superflares. We also estimated the stellar parameters, such as effective temperature, surface gravity, metallicity, projected rotational velocity (v sin_i), and radial velocity. KIC6934317 is then confirmed to be an early G-type main sequence star. The value of v sin_i is estimated to be ~1.91 km s^-1. In contrast, the rotational velocity is calculated to be ~20 km s^-1 by using the period of the brightness variation as the rotation period. This difference can be explained by its small inclination angle (nearly pole-on). The small inclination angle is also supported by the contrast between the large superflare amplitude and the small stellar brightness variation amplitude. The lithium abundance and isochrones implies that the age of this star is more than about a few Gyr, though a problem why this star with such an age has a strong activity remains unsolved.

preprint2013arXiv

Superflares on Solar Type Stars Observed with Kepler I. Statistical Properties of Superflares

By extending our previous study by Maehara et al. (2012), we searched for superflares on G-type dwarfs (solar type stars) using Kepler data for a longer period (500 days) than that (120 days) in our previous study. As a result, we found 1547 superflares on 279 G-type dwarfs, which are much more than previous 365 superflares on 148 stars. Using these new data, we studied the statistical properties of occurrence frequency of superflares, and basically confirmed the previous results, i.e., the occurrence frequency (dN/dE) of superflares vs flare energy (E) shows power-law distribution with dN/dE \propto E^{-α}, where α~ 2. It is interesting that this distribution is roughly on the same line as that for solar flares. In the case of the Sun-like stars (with surface temperature 5600-6000K and slowly rotating with a period longer than 10 days), the occurrence frequency of superflares with energy of 10^34 -10^35 erg is once in 800-5000 years. We also studied long term (500 days) stellar brightness variation of these superflare stars, and found that in some G-type dwarfs the occurrence frequency of superflares was extremely high, ~ 57 superflares in 500 days (i.e., once in 10 days). In the case of Sun-like stars, the most active stars show the frequency of one superflares (with 10^34 erg) in 100 days. There is evidence that these superflares have extremely large starspots with a size about 10 times larger than that of the largest sunspot. We argue that the physical origin of extremely high occurrence frequency of superflares in these stars may be attributed to the existence of extremely large starspots.

preprint2013arXiv

Superflares on Solar-Type Stars Observed with Kepler II. Photometric Variability of Superflare-Generating Stars : A Signature of Stellar Rotation and Starspots

We performed simple spot-model calculations for quasi-periodic brightness variations of solar-type stars showing superflares, by using Kepler photometric data. Most of superflare stars show quasi-periodic brightness modulations with the typical period of one to a few tens of days. Our results indicate that these brightness variations of superflare stars can be explained by the rotation of the star with fairly large starspots. Using the result of the period analysis, we investigated the relation between the energy and frequency of superflares and the rotation period. Stars with relatively slower rotation rates can still produce flares that are as energetic as those of more rapidly rotating stars, although the average flare frequency is lower for more slowly rotating stars. We found that the energy of superflares are related to the total coverage of starspots. The correlation between the spot coverage and the flare energy in superflares is similar to that in solar flares. These results suggest that the energy of superflares can be explained by the magnetic energy stored around starspots.