We report multi-wavelength monitoring observations of an M-dwarf flare star AD Leonis with Seimei Telescope (6150--7930 {\AA}), SCAT (Spectroscopic Chuo-university Astronomical Telescope; 3700--7500 {\AA}), 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$\alpha$, 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$\alpha$ emission line full width at 1/8 maximum dramatically increases to 14 {\AA} from 8 {\AA} in the low-resolution spectra (R$\sim$ 2000) accompanied with the large white-light flares, (2) some weak H$\alpha$/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$\alpha$ 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$\alpha$ 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$\alpha$ 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.

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.

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.

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{\deg} 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.

Quasi-periodic pulsations (QPPs) appear to be a common feature observed in the light curves of both solar and stellar flares. However, their quasi-periodic nature, along with the fact that they can be small in amplitude and short-lived, makes QPPs difficult to unequivocally detect. In this paper, we test the strengths and limitations of state-of-the-art methods for detecting QPPs using a series of hare-and-hounds exercises. The hare simulated a set of flares, both with and without QPPs of a variety of forms, while the hounds attempted to detect QPPs in blind tests. We use the results of these exercises to create a blueprint for anyone who wishes to detect QPPs in real solar and stellar data. We present eight clear recommendations to be kept in mind for future QPP detections, with the plethora of solar and stellar flare data from new and future satellites. These recommendations address the key pitfalls in QPP detection, including detrending, trimming data, accounting for colored noise, detecting stationary-period QPPs, detecting QPPs with nonstationary periods, and ensuring that detections are robust and false detections are minimized. We find that QPPs can be detected reliably and robustly by a variety of methods, which are clearly identified and described, if the appropriate care and due diligence are taken.

High-energy radiation caused by exoplanetary space weather events from planet-hosting stars can play a crucial role in conditions promoting or destroying habitability in addition to the conventional factors. In this paper, we present the first quantitative impact evaluation system of stellar flares on the habitability factors with an emphasis on the impact of Stellar Proton Events. We derive the maximum flare energy from stellar starspot sizes and examine the impacts of flare associated ionizing radiation on CO$_2$, H$_2$, N$_2$+O$_2$ --rich atmospheres of a number of well-characterized terrestrial type exoplanets. Our simulations based on the Particle and Heavy Ion Transport code System [PHITS] suggest that the estimated ground level dose for each planet in the case of terrestrial-level atmospheric pressure (1 bar) for each exoplanet does not exceed the critical dose for complex (multi-cellular) life to persist, even for the planetary surface of Proxima Centauri b, Ross-128 b and TRAPPIST-1 e. However, when we take into account the effects of the possible maximum flares from those host stars, the estimated dose reaches fatal levels at the terrestrial lowest atmospheric depth on TRAPPIST-1 e and Ross-128 b. Large fluxes of coronal XUV radiation from active stars induces high atmospheric escape rates from close-in exoplanets suggesting that the atmospheric depth can be substantially smaller than that on the Earth. In a scenario with the atmospheric thickness of 1/10 of Earth's, the radiation dose from close-in planets including Proxima Centauri b and TRAPPIST-1 e reach near fatal dose levels with annual frequency of flare occurrence from their hoststars.

The current progress in the detection of terrestrial type exoplanets has opened a new avenue in the characterization of exoplanetary atmospheres and in the search for biosignatures of life with the upcoming ground-based and space missions. To specify the conditions favorable for the origin, development and sustainment of life as we know it in other worlds, we need to understand the nature of astrospheric, atmospheric and surface environments of exoplanets in habitable zones around G-K-M dwarfs including our young Sun. Global environment is formed by propagated disturbances from the planet-hosting stars in the form of stellar flares, coronal mass ejections, energetic particles, and winds collectively known as astrospheric space weather. Its characterization will help in understanding how an exoplanetary ecosystem interacts with its host star, as well as in the specification of the physical, chemical and biochemical conditions that can create favorable and/or detrimental conditions for planetary climate and habitability along with evolution of planetary internal dynamics over geological timescales. A key linkage of (astro) physical, chemical, and geological processes can only be understood in the framework of interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary and Earth sciences. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets will significantly expand the current definition of the habitable zone to the biogenic zone and provide new observational strategies for searching for signatures of life. The major goal of this paper is to describe and discuss the current status and recent progress in this interdisciplinary field and to provide a new roadmap for the future development of the emerging field of exoplanetary science and astrobiology.

We report the latest view of Kepler solar-type (G-type main-sequence) superflare stars, including recent updates with Apache Point Observatory (APO) 3.5m telescope spectroscopic observations and Gaia-DR2 data. First, we newly conducted APO3.5m spectroscopic observations of 18 superflare stars found from Kepler 1-min time cadence data. More than half (43 stars) are confirmed to be "single" stars, among 64 superflare stars in total that have been spectroscopically investigated so far in this APO3.5m and our previous Subaru/HDS observations. The measurements of $v\sin i$ (projected rotational velocity) and chromospheric lines (Ca II H\&K and Ca II 8542\AA) support the brightness variation of superflare stars is caused by the rotation of a star with large starspots. We then investigated the statistical properties of Kepler solar-type superflare stars by incorporating Gaia-DR2 stellar radius estimates. As a result, the maximum superflare energy continuously decreases as the rotation period $P_{\mathrm{rot}}$ increases. Superflares with energies $\lesssim 5\times10^{34}$ erg occur on old, slowly-rotating Sun-like stars ($P_{\mathrm{rot}}\sim$25 days) approximately once every 2000--3000 years, while young rapidly-rotating stars with $P_{\mathrm{rot}}\sim$ a few days have superflares up to $10^{36}$ erg. The maximum starspot area does not depend on the rotation period when the star is young, but as the rotation slows down, it starts to steeply decrease at $P_{\mathrm{rot}}\gtrsim$12 days for Sun-like stars. These two decreasing trends are consistent since the magnetic energy stored around starspots explains the flare energy, but other factors like spot magnetic structure should also be considered.

The field of exoplanetary science is making rapid progress both in statistical studies of exoplanet properties as well as in individual characterization. As space missions provide an emerging picture of formation and evolution of exoplanetary systems, the search for habitable worlds becomes one of the fundamental issues to address. To tackle such a complex challenge, we need to specify the conditions favorable for the origin, development and sustainment of life as we know it. This requires the understanding of global (astrospheric) and local (atmospheric, surface and internal) environments of exoplanets in the framework of the physical processes of the interaction between evolving planet-hosting stars along with exoplanetary evolution over geological timescales, and the resulting impact on climate and habitability of exoplanets. Feedbacks between astrophysical, physico-chemical atmospheric and geological processes can only be understood through interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary, Earth sciences, astrobiology, and the origin of life communities. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets and potential exomoons around them may significantly modify the extent and the location of the habitable zone and provide new directions for searching for signatures of life. Thus, characterization of stellar ionizing outputs becomes an important task for further understanding the extent of habitability in the universe. The goal of this white paper is to identify and describe promising key research goals to aid the theoretical characterization and observational detection of ionizing radiation from quiescent and flaring upper atmospheres of planet hosts as well as properties of stellar coronal mass ejections and stellar energetic particle events.

Active solar-type stars show large quasi-periodic brightness variations caused by stellar rotations with star spots, and the amplitude changes as the spots emerge and decay. The Kepler data are suitable for investigations on the emergence and decay processes of star spots, which are important to understand underlying stellar dynamo and stellar flares. In this study, we measured temporal evolutions of star spot area with Kepler data by tracing local minima of the light curves. In this analysis, we extracted temporal evolutions of star spots showing clear emergence and decay without being disturbed by stellar differential rotations. We applied this method to 5356 active solar-type stars observed by Kepler and obtained temporal evolutions of 56 individual star spots. We calculated lifetimes, emergence and decay rates of the star spots from the obtained temporal evolutions of spot area. As a result, we found that lifetimes ($T$) of star spots are ranging from 10 to 350 days when spot areas ($A$) are 0.1-2.3 percent of the solar hemisphere. We also compared them with sunspot lifetimes, and found that the lifetimes of star spots are much shorter than those extrapolated from an empirical relation of sunspots ($T\propto A$), while being consistent with other researches on star spot lifetimes. The emerging and decay rates of star spots are typically $5 \times 10^{20}$ $\rm Mx\cdot h^{-1}$ ($8$ $\rm MSH\cdot h^{-1}$) with the area of 0.1-2.3 percent of the solar hemisphere and are mostly consistent with those expected from sunspots, which may indicate the same underlying processes.

We have performed 5 night spectroscopic observation of the Halpha line of EV Lac with a medium wavelength resolution (R~ 10,000) using the 2m Nayuta telescope at the Nishi-Harima Astronomical Observatory. EV Lac always possesses the Halpha emission line; however, its intensity was stronger on August 15, 2015 than during other four-night periods. On this night, we observed a rapid rise (~ 20min) and a subsequent slow decrease (~ 1.5h) of the emission-line intensity of Halpha, which was probably caused by a flare. We also found an asymmetrical change in the Halpha line on the same night. The enhancement has been observed in the blue wing of the Ha line during each phase of this flare (from the flare start to the flare end), and absorption components were present in its red wing during the early and later phases of the flare. Such blue enhancement (blue asymmetry) of the Halpha line is sometimes seen during solar flares, but only during the early phases. Even for solar flares, little is known about the origin of the blue asymmetry. Compared with solar-flare models, the presented results can lead to the understanding of the dynamics of stellar flares.

Recently, many superflares on solar-type stars were discovered as white-light flares (WLFs). A correlation between the energies (E) and durations (t) of superflares is derived as $t\propto E^{0.39}$, and this can be theoretically explained by magnetic reconnection ($t\propto E^{1/3}$). In this study, we carried out a statistical research on 50 solar WLFs with SDO/HMI to examine the t-E relation. As a result, the t-E relation on solar WLFs ($t\propto E^{0.38}$) is quite similar stellar superflares, but the durations of stellar superflares are much shorter than those extrapolated from solar WLFs. We present the following two interpretations; (1) in solar flares, the cooling timescale of WL emission may be longer than the reconnection one, and the decay time can be determined by the cooling timescale; (2) the distribution can be understood by applying a scaling law $t\propto E^{1/3}B^{-5/3}$ derived from the magnetic reconnection theory.

It is currently unknown how common life is on exoplanets, or how long planets can remain viable for life. To date, we have a superficial notion of habitability, a necessary first step, but so far lacking an understanding of the detailed interaction between stars and planets over geological timescales, dynamical evolution of planetary systems, and atmospheric evolution on planets in other systems. A planet mass, net insolation, and atmospheric composition alone are insufficient to determine the probability that life on a planet could arise or be detected. The latter set of planetary considerations, among others, underpin the concept of the habitable zone (HZ), defined as the circumstellar region where standing bodies of liquid water could be supported on the surface of a rocky planet. However, stars within the same spectral class are often treated in the same way in HZ studies, without any regard for variations in activity among individual stars. Such formulations ignore differences in how nonthermal emission and magnetic energy of transient events in different stars affect the ability of an exoplanet to retain its atmosphere.In the last few years there has been a growing appreciation that the atmospheric chemistry, and even retention of an atmosphere in many cases, depends critically on the high-energy radiation and particle environments around these stars. Indeed, recent studies have shown stellar activity and the extreme space weather, such as that created by the frequent flares and coronal mass ejections (CMEs) from the active stars and young Sun, may have profoundly affected the chemistry and climate and thus habitability of the early Earth and terrestrial type exoplanets. The goal of this white paper is to identify and describe promising key research goals to aid the field of the exoplanetary habitability for the next 20 years.

The search of life in the Universe is a fundamental problem of astrobiology and a major priority for NASA. A key area of major progress since the NASA Astrobiology Strategy 2015 (NAS15) has been a shift from the exoplanet discovery phase to a phase of characterization and modeling of the physics and chemistry of exoplanetary atmospheres, and the development of observational strategies for the search for life in the Universe by combining expertise from four NASA science disciplines including heliophysics, astrophysics, planetary science and Earth science. The NASA Nexus for Exoplanetary System Science (NExSS) has provided an efficient environment for such interdisciplinary studies. Solar flares, coronal mass ejections and solar energetic particles produce disturbances in interplanetary space collectively referred to as space weather, which interacts with the Earth upper atmosphere and causes dramatic impact on space and ground-based technological systems. Exoplanets within close in habitable zones around M dwarfs and other active stars are exposed to extreme ionizing radiation fluxes, thus making exoplanetary space weather (ESW) effects a crucial factor of habitability. In this paper, we describe the recent developments and provide recommendations in this interdisciplinary effort with the focus on the impacts of ESW on habitability, and the prospects for future progress in searching for signs of life in the Universe as the outcome of the NExSS workshop held in Nov 29 - Dec 2, 2016, New Orleans, LA. This is one of five Life Beyond the Solar System white papers submitted by NExSS to the National Academy of Sciences in support of the Astrobiology Science Strategy for the Search for Life in the Universe.

Recently, many superflares on solar-type stars have been discovered as white-light flares (WLFs). The statistical study found a correlation between their energies ($E$) and durations ($\tau$): $\tau \propto E^{0.39}$ (Maehara et al. 2017 $EP\& S$, 67, 59), similar to those of solar hard/soft X-ray flares: $\tau \propto E^{0.2-0.33}$. This indicates a universal mechanism of energy release on solar and stellar flares, i.e., magnetic reconnection. We here carried out a statistical research on 50 solar WLFs observed with \textit{SDO}/HMI and examined the correlation between the energies and durations. As a result, the $E$--$\tau$ relation on solar WLFs ($\tau \propto E^{0.38}$) is quite similar to that on stellar superflares ($\tau \propto E^{0.39}$). However, the durations of stellar superflares are one order of magnitude shorter than those expected from solar WLFs. We present the following two interpretations for the discrepancy. (1) In solar flares, the cooling timescale of WLFs may be longer than the reconnection one, and the decay time of solar WLFs can be elongated by the cooling effect. (2) The distribution can be understood by applying a scaling law ($\tau \propto E^{1/3}B^{-5/3}$) derived from the magnetic reconnection theory. In this case, the observed superflares are expected to have 2-4 times stronger magnetic field strength than solar flares.

We analyze the correlation between starspots and superflares on solar-type stars using observations from the Kepler mission. The analysis shows that the observed fraction of stars with superflares decreases as the rotation period increases and as the amplitude of photometric variability associated with rotation decreases. We found that the fraction of stars with superflares among the stars showing large-amplitude rotational variations, which are thought to be the signature of the large starspots, also decreases as the rotation period increases. The small fraction of superflare stars among the stars with large starspots in the longer-period regime suggests that some of the stars with large starspots show a much lower flare activity than the superflare stars with the same spot area. Assuming simple relations between spot area and lifetime and between spot temperature and photospheric temperature, we compared the size distribution of large starspot groups on slowly-rotating solar-type stars with that of sunspot groups. The size distribution of starspots shows the power-law distribution and the size distribution of larger sunspots lies on this power-law line. We also found that frequency-energy distributions for flares originating from spots with different sizes are the same for solar-type stars with superflares and the Sun. These results suggest that the magnetic activity we observe on solar-type stars with superflares and that on the Sun is caused by the same physical processes.

Recent studies of solar-type superflare stars have suggested that even old slowly rotating stars similar to the Sun can have large starspots and superflares. We conducted high dispersion spectroscopy of 49 nearby solar-analog stars (G-type main sequence stars with $T_{\rm{eff}}\approx5,600\sim6,000$ K) identified as ROSAT soft X-ray sources, which are not binary stars from the previous studies. We expected that these stars can be used as a proxy of bright solar-analog superflare stars, since superflare stars are expected to show strong X-ray luminosity. More than half (37) of the 49 target stars show no evidence of binarity, and their atmospheric parameters ($T_{\rm{eff}}$, $\log g$, and [Fe/H]) are within the range of ordinary solar-analog stars. We measured Ca II 8542 and H$\alpha$ lines, which are good indicators of the chromospheric activity. The intensity of these lines indicates that all the target stars have large starspots. We also measured $v\sin i$ (projected rotational velocity) and Lithium abundance for the target stars. Li abundance is a key to understanding the evolution of the stellar convection zone, which reflects the stellar age, mass and rotational history. We confirmed that many of the target stars rapidly rotate and have high Li abundance, compared with the Sun, as suggested by many previous studies. There are, however, also some target stars that rotate slowly ($v\sin i$=2$\sim$3 km s$^{-1}$) and have low Li abundance like the Sun. These results support that old and slowly rotating stars similar to the Sun could have high activity level and large starspots. This is consistent with the results of our previous studies of solar-type superflare stars. In the future, it is important to conduct long-term monitoring observations of these active solar-analog stars in order to investigate detailed properties of large starspots from the viewpoint of stellar dynamo theory.

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.

We searched for superflares on solar-type stars using the Kepler short-cadence (1-min sampling) data 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. 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 ^{-\alpha}$) with $\alpha=1.5$ for $10^{33}<E<10^{36}$ erg. 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 amplitude of brightness variations. We also found that the duration of superflares ($\tau$) increases with the flare energy ($E$) as $\tau \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.

We carried out spectroscopic observations with Subaru/HDS of 50 solar-type superflare stars found from Kepler data. More than half (34 stars) of the target stars show no evidence of the binary system, and we confirmed atmospheric parameters of these stars are roughly in the range of solar-type stars. We then conducted the detailed analyses for these 34 stars. First, the value of the "$v\sin i$" (projected rotational velocity) measured from spectroscopic results is consistent with the rotational velocity estimated from the brightness variation. Second, there is a correlation between the amplitude of the brightness variation and the intensity of Ca II IR triplet line. All the targets expected to have large starspots because of their large amplitude of the brightness variation show high chromospheric activities compared with the Sun. These results support that the brightness variation of superflare stars is explained by the rotation of a star with large starspots.

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.

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^{-\alpha}$) with $\alpha \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 ($\tau$) increases with the flare energy ($E$) as $\tau \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.

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 ($(\Delta T_{\rm{eff}})_{\rm{rms}} \sim 219$K, $(\Delta \log g)_{\rm{rms}} \sim 0.37$ dex, and $(\Delta\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.

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$\alpha$ 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.

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.

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^{-\alpha}, where \alpha ~ 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.

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.

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.

Recent observations of solar type stars with the Kepler satellite by Maehara et al. have revealed the existence of superflares (with energy of 10^33 - 10^35 erg) on Sun-like stars, which are similar to our Sun in their surface temperature (5600 K - 6000 K) and slow rotation (rotational period > 10 days). From the statistical analysis of these superflares, it was found that superflares with energy 10^34 erg occur once in 800 years and superflares with 10^35 erg occur once in 5000 years on Sun-like stars. In this paper, we examine whether superflares with energy of 10^33 - 10^35 erg could occur on the present Sun through the use of simple order-of-magnitude estimates based on current ideas relating to the mechanisms of the solar dynamo.