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.

Measurements of the gas mass are necessary to determine the planet formation potential of protoplanetary disks. Observations of rare CO isotopologues are typically used to determine disk gas masses; however, if the line emission is optically thick this will result in an underestimated disk mass. With ALMA we have detected the rarest stable CO isotopologue, 13C17O, in a protoplanetary disk for the first time. We compare our observations with the existing detections of 12CO, 13CO, C18O and C17O in the HD163296 disk. Radiative transfer modelling using a previously benchmarked model, and assuming interstellar isotopic abundances, significantly underestimates the integrated intensity of the 13C17O J=3-2 line. Reconciliation between the observations and the model requires a global increase in CO gas mass by a factor of 3.5. This is a factor of 2-6 larger than previous gas mass estimates using C18O. We find that C18O emission is optically thick within the CO snow line, while the 13C17O emission is optically thin and is thus a robust tracer of the bulk disk CO gas mass.

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.

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.

In this paper, we analyse the upper limit fluxes of sub-millimeter ortho-H$_{2}$$^{16}$O 321 GHz, para-H$_{2}$$^{18}$O 322 GHz, and HDO 335 GHz lines from the protoplanetary disk around the Herbig Ae star HD 163296, using the Atacama Large Millimeter/Submillimeter Array (ALMA). These water lines are considered to be the best candidate sub-millimeter lines to locate the position of the H$_{2}$O snowline, on the basis of our previous model calculations. We compare the upper limit fluxes with the values calculated by our models with dust emission included, and we constrain the line emitting region and the dust opacity from the observations. We conclude that, if the outer edge of the region with high water vapor abundance and if the position of the water snowline are beyond 8 au also, the mm dust opacity will have a value larger than 2.0 cm$^{2}$ g$^{-1}$. In addition, the position of the water snowline will be inside 20 au, if the mm dust opacity is 2.0 cm$^{2}$ g$^{-1}$. Future observations of the dust continuum emission at higher angular resolution and sub-millimeter water lines with longer observation time are required to clarify the detailed structures and the position of the H$_{2}$O snowline in the disk midplane.

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.

In this paper, we extend the results presented in our former papers (Notsu et al. 2016, 2017) on using ortho-$\mathrm{H_2}$$^{16}\mathrm{O}$ line profiles to constrain the location of the $\mathrm{H_2O}$ snowline in T Tauri and Herbig Ae disks, to include sub-millimeter para-$\mathrm{H_2}$$^{16}\mathrm{O}$ and ortho- and para-$\mathrm{H_2}$$^{18}\mathrm{O}$ lines. Since the number densities of the ortho- and para-H$_{2}$$^{18}$O molecules are about 560 times smaller than their $^{16}$O analogues, they trace deeper into the disk than the ortho-H$_{2}$$^{16}$O lines (down to $z=0$, i.e., the midplane). Thus these H$_{2}$$^{18}$O lines are potentially better probes of the position of the H$_{2}$O snowline at the disk midplane, depending on the dust optical depth. The values of the Einstein $A$ coefficients of sub-millimeter candidate water lines tend to be lower (typically $<$$10^{-4}$ s$^{-1}$) than infrared candidate water lines (Notsu et al. 2017). Thus in the sub-millimeter candidate water line cases, the local intensity from the outer optically thin region in the disk is around $10^{4}$ times smaller than that in the infrared candidate water line cases. Therefore, in the sub-millimeter lines, especially H$_{2}$$^{18}$O and para-H$_{2}$$^{16}$O lines with relatively lower upper state energies ($\sim$ a few 100K) can also locate the position of the $\mathrm{H_2O}$ snowline. We also investigate the possibility of future observations with ALMA to identify the position of the water snowline. There are several candidate water lines that trace the hot water vapor inside the $\mathrm{H_2O}$ snowline in ALMA Bands $5-10$.

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.

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.

Observationally measuring the location of the H$_{2}$O snowline is crucial for understanding the planetesimal and planet formation processes, and the origin of water on Earth. In disks around Herbig Ae stars ($T_{\mathrm{*}}\sim$ 10,000K, $M_{\mathrm{*}}\gtrsim$ 2.5$M_{\bigodot}$), the position of the H$_{2}$O snowline is further from the central star compared with that around cooler, and less massive T Tauri stars. Thus, the H$_{2}$O emission line fluxes from the region within the H$_{2}$O snowline are expected to be stronger. In this paper, we calculate the chemical composition of a Herbig Ae disk using chemical kinetics. Next, we calculate the H$_{2}$O emission line profiles, and investigate the properties of candidate water lines across a wide range of wavelengths (from mid-infrared to sub-millimeter) that can locate the position of the H$_{2}$O snowline. Those line identified have small Einstein $A$ coefficients ($\sim 10^{-6} -10^{-3}$ s$^{-1}$) and relatively high upper state energies ($\sim$ 1000K). The total fluxes tend to increase with decreasing wavelengths. We investigate the possibility of future observations (e.g., ALMA, SPICA/SMI-HRS) to locate the position of the H$_{2}$O snowline. Since the fluxes of those identified lines from Herbig Ae disks are stronger than those from T Tauri disks, the possibility of a successful detection is expected to increase for a Herbig Ae disk.

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.

Inside the H$_{2}$O snowline of protoplanetary disks, water evaporates from the dust-grain surface into the gas phase, whereas it is frozen out on to the dust in the cold region beyond the snowline. H$_{2}$O ice enhances the solid material in the cold outer part of a disk, which promotes the formation of gas-giant planet cores. We can regard the H$_{2}$O snowline as the surface that divides the regions between rocky and gaseous giant planet formation. Thus observationally measuring the location of the H$_{2}$O snowline is crucial for understanding the planetesimal and planet formation processes, and the origin of water on Earth. In this paper, we find candidate water lines to locate the H$_{2}$O snowline through future high-dispersion spectroscopic observations. First, we calculate the chemical composition of the disk and investigate the abundance distributions of H$_{2}$O gas and ice, and the position of the H$_{2}$O snowline. We confirm that the abundance of H$_{2}$O gas is high not only in the hot midplane region inside the H$_{2}$O snowline but also in the hot surface layer of the outer disk. Second, we calculate the H$_{2}$O line profiles and identify those H$_{2}$O lines which are promising for locating the H$_{2}$O snowline: the identified lines are those which have small Einstein $A$ coefficients and high upper state energies. The wavelengths of the candidate H$_{2}$O lines range from mid-infrared to sub-millimeter, and they overlap with the regions accessible to ALMA and future mid-infrared high dispersion spectrographs (e.g., TMT/MICHI, SPICA).

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.