Search

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., Comment: 37 pages, 12 figures, 1 table. Accepted for publication in The Astrophysical Journal

Observationally locating the position of the H\(_{2}\)O snowline in protoplanetary disks is important to understand the processes of planetesimal and planet formation, and the origin of water on terrestrial planets including the Earth. In this thesis, we proposed the method to locate the H\(_{2}\)O snowline position more directly by analyzing the Keplerian water line profiles which can be obtained by high dispersion spectroscopic observations across a wide range of wavelengths (from mid-infrared to sub-millimeter, e.g., ALMA, SPICA) and selected based on specific criteria. In addition, we also reported the results of water line observations for the protoplanetary disk around the Herbig Ae star HD 163296, and constrained the dust opacity and the line emitting region from the observations. Most contents of this thesis is based on our refereed papers which have been published (Notsu et al. 2016, ApJ, 827, 113; 2017, ApJ, 836, 118; 2018, ApJ, 855, 62; 2019, ApJ, 875, 96).

Observationally locating the position of the H\(_{2}\)O snowline in protoplanetary disks is important to understand the processes of planetesimal and planet formation, and the origin of water on terrestrial planets including the Earth. In our studies, first we calculated chemical structures of the disk using the self-consistent physical models of a typical Herbig Ae disk. Next, on the basis of our calculations of disk chemical structures and water line profiles, we proposed how to identify the H\(_{2}\)O snowline positions directly by analyzing the Keplerian water line profiles which can be obtained by high-dispersion spectroscopic observations across a wide range of wavelengths (from mid-infrared to sub-millimeter wavelengths). We selected candidate water lines to locate the H\(_{2}\)O snowline based on specific criteria. We concluded that water lines with small Einstein A coefficients (A\(_{\mathrm {ul}} = 10^{-6}{\sim } 10^{-3}\) s\(^{-1}\)) and relatively high upper state energies (E\(_{\mathrm {up}}\) \(\sim \) 1000 K) trace the hot water vapor within the \({\mathrm {H_2O}}\) snowline, and can locate the H\(_{2}\)O snowline positions. In these candidate water lines, the contribution of the optically thick hot midplane within the H\(_{2}\)O snowline is large compared with that of the outer optically thin surface layer. This is because the line intensities from the optically thin region are proportional to the Einstein A coefficient. In addition, the contribution of the cold water reservoir outside the H\(_{2}\)O snowline is also small, since lines with high excitation energies are not emitted from the regions at a low temperature. The H\(_{2}\)O snowline positions of a Herbig Ae disk exists at a larger radius compared with that around cooler and less massive T Tauri stars. Moreover, the H\(_{2}\)O snowline position migrates closer to the star as the disk becomes older and mass accretion rate to the central star becomes smaller. Thus, observing the candidate water lines and locating the H\(_{2}\)O snowline positions, in Herbig Ae disks and younger T Tauri stars, is expected to be easier. We investigate the possibility of future observations (e.g., ALMA, SPICA/SMI-HRS) to locate the \({\mathrm {H_2O}}\) snowline position. Most contents of this chapter is based on our refereed paper that has been published (Notsu et al. 2017, ApJ, 836, 118).

Starspots are thought to be regions of locally strong magnetic fields, similar to sunspots, and they can generate photometric brightness modulations. To deduce stellar and spot properties, such as spot emergence and decay rates, we implement computational code for starspot modeling. It is implemented with an adaptive parallel tempering algorithm and an importance sampling algorithm for parameter estimation and model selection in the Bayesian framework. For evaluating the performance of the code, we apply it to synthetic light curves produced with 3 spots. The light curves are specified in the spot parameters, such as the radii, intensities, latitudes, longitudes, and emergence/decay durations. The spots are circular with specified radii and intensities relative to the photosphere, and the stellar differential rotation coefficient is also included in the light curves. As a result, stellar and spot parameters are uniquely deduced. The number of spots is correctly determined: the 3-spot model is preferable because the model evidence is much greater than that of 2-spot models by orders of magnitude and more than that of 4-spot model by a more modest factor, whereas the light curves are produced to have 2 or 1 local minimum during one equatorial rotation period by adjusting the values of longitude. The spot emergence and decay rates can be estimated with error less than an order of magnitude, considering the difference of the number of spots., Comment: 27 pages, 14 figures, 2 tables, accepted for publication in ApJ

We found candidate water lines to locate the \(\mathrm {H_2O}\) snowline position through future high-dispersion spectroscopic observations. As a first step, we calculated chemical structures of the disk using the self-consistent physical models of a typical T Tauri disk. We confirmed that the water gas abundance is high not only in the hot disk midplane within the \(\mathrm {H_2O}\) snowline, but also in the outer hot surface and photodesorption region. Next, we calculated the profiles of water lines, and find the lines which are best to locate the position of the \(\mathrm {H_2O}\) snowline. The lines we identified are those with high upper state energies and small Einstein A coefficients. The wavelengths of the candidate water lines range from mid-infrared to sub-millimeter, and they overlap with the wavelength coverages of ALMA and future mid-infrared high dispersion spectrographs (e.g., TMT/MICHI, SPICA). Most contents of this chapter is based on our refereed paper that has been published (Notsu et al. 2016, ApJ, 827, 113).

In this chapter, we extend the results presented in our former papers [41, 42] 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 1/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\)), and lines with relatively smaller upper state energies (\(\sim \)a few 100 K) can also locate the \(\mathrm {H_2O}\) snowline positions. Thus these H\(_{2}\) \(^{18}\)O lines are potentially better probes of the H\(_{2}\)O snowline positions 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. 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 100 K) can also locate the position of the \(\mathrm {H_2O}\) snowline. We also investigate the possibility of future observations with ALMA to identify the water snowline position. There are several candidate water lines that trace the hot water vapor inside the \(\mathrm {H_2O}\) snowline in ALMA Bands 5–10. Most contents of this chapter is based on our refereed paper that has been published (Notsu et al. 2018, ApJ, 855, 62).

We have conducted the water line observations for the protoplanetary disk around the Herbig Ae star HD 163296, and partial data were delivered. We analyzed 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, which are considered to be the best candidate water lines available at sub-millimeter wavelengths to locate the \(\mathrm {H_2O}\) snowline positions, according to our previous model studies. We compared the upper limit fluxes with the values obtained using our models with dust emission, and constrained the line emitting region and the dust opacity from observations. We concluded that, if the outer edge of the water vapor abundant region and also the water snowline position is beyond 8 au, the 1 mm dust opacity \(\kappa _{\mathrm {mm}}\) will have a value larger than 2.0 cm\(^{2}\) g\(^{-1}\). Moreover, if \(\kappa _{\mathrm {mm}}\) is 2.0 cm\(^{2}\) g\(^{-1}\), the water snowline position will be inside 20 au. The multiple gap and ring structures in 0.9 mm (Band 7) dust continuum emission with the spatial resolution of 15 au are also reported. The gap and ring positions are consistent with those estimated in the previous observations. Most contents of this chapter is based on our refereed paper that has been published (Notsu et al. 2019, ApJ, 875, 96).

The radial-dependent positions of snowlines of abundant oxygen- and carbon-bearing molecules in protoplanetary discs will result in systematic radial variations in the C/O ratios in the gas and ice. This variation is proposed as a tracer of the formation location of gas-giant planets. However, disc chemistry can affect the C/O ratios in the gas and ice, thus potentially erasing the chemical fingerprint of snowlines in gas-giant atmospheres. We calculate the molecular composition of hot Jupiter atmospheres using elemental abundances extracted from a chemical kinetics model of a disc midplane where we have varied the initial abundances and ionization rates. The models predict a wider diversity of possible atmospheres than those predicted using elemental ratios from snowlines only. As found in previous work, as the C/O ratio exceeds the solar value, the mixing ratio of CH$_{4}$ increases in the lower atmosphere, and those of C$_{2}$H$_{2}$ and HCN increase mainly in the upper atmosphere. The mixing ratio of H$_{2}$O correspondingly decreases. We find that hot Jupiters with C/O$>1$ can only form between the CO$_{2}$ and CH$_{4}$ snowlines. Moreover, they can only form in a disc which has fully inherited interstellar abundances, and where negligible chemistry has occurred. Hence, carbon-rich planets are likely rare, unless efficient transport of hydrocarbon-rich ices via pebble drift to within the CH$_{4}$ snowline is a common phenomenon. We predict combinations of C/O ratios and elemental abundances that can constrain gas-giant planet formation locations relative to snowline positions, and that can provide insight into the disc chemical history., 18 pages, 7 figures, and 2 tables are contained in this paper. It was accepted for publication in MNRAS on September 18th, 2020

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., 4 pages, 1 figures, paper of the Proceedings of the IAUS340, accepted

Observationally measuring the location of the H2O snowline is crucial for understanding the planetesimal and planet formation processes, and the origin of water on Earth. The velocity profiles of emission lines from protoplanetary disks are usually affected by Doppler shift due to Keplerian rotation and thermal broadening. Therefore, the velocity profiles are sensitive to the radial distribution of the line-emitting regions. In our work (Notsu et al. 2016, 2017), we found candidate water lines to locate the position of the H2O snowline through future high-dispersion spectroscopic observations. First, we calculated the chemical composition of the disks around a T Tauri star and a Herbig Ae star using chemical kinetics. We confirmed that the abundance of H2O gas is high not only in the hot midplane region inside the H2O snowline but also in the hot surface layer and the photodesorption region of the outer disk. The position of the H2O snowline in the Herbig Ae disk exists at a larger radius from the central star than that in the T Tauri disk. Second, we calculated the H2O line profiles and identified that H2O emission lines with small Einstein A coefficients (∼10−6 − 10−3 s−1) and relatively high upper state energies (∼ 1000K) are dominated by emission from the hot midplane region inside the H2O snowline, and therefore their profiles potentially contain information which can be used to locate the position of the H2O snowline. The wavelengths of the H2O lines which are the best candidates to locate the position of the H2O snowline range from mid-infrared to sub-millimeter, and the total line fluxes tend to increase with decreasing wavelengths. We investigated the possibility of future observations using the ALMA and mid-infrared high-dispersion spectrographs (e.g., SPICA/SMI-HRS). Since the fluxes of those identified lines from a Herbig Ae disk are stronger than those of a T Tauri disk, the possibility of a successful detection is expected to increase for a Herbig Ae disk.

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., Comment: Accepted ApJL

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., 37 pages, 19 figures, 4 tables. Accepted for publication in The Astrophysical Journal (on June 16, 2019), Version 2 (fixed typo)

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., 17 pages, 8 figures, and 3 tables are contained in this paper. It was accepted for publication in The Astrophysical Journal (ApJ) on February 26th, 2019

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., 17 Pages, 9 figures, Accepted for Publication in PASJ

In this paper, we extend the results presented in our former papers on using ortho-H216O line profiles to constrain the location of the H2O snowline in T Tauri and Herbig Ae disks, to include submillimeter para-H216O and ortho- and para-H218O lines. Since the number densities of the ortho- and para-H218O molecules are about 560 times smaller than their 16O analogs, they trace deeper into the disk than the ortho-H216O lines (down to z = 0, i.e., the midplane). Thus these H218O lines are potentially better probes of the position of the H2O snowline at the disk midplane, depending on the dust optical depth. The values of the Einstein A coefficients of submillimeter candidate water lines tend to be lower (typically

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., Comment: 19 pages, 11 figures, accepted for publication in ApJ

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., Accepted for publication in PASJ

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}, Comment: 6 pages, 4 figures, submitted to the proceedings of the IAU Symposium 320 "Solar and Stellar Flares and Their Effects on Planets" (eds. A. G. Kosovichev, S. L. Hawley, P. Heinzel)

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., Comment: 52 pages, 18 figures, 6 tables, accepted for publication in PASJ. arXiv admin note: text overlap with arXiv:1412.8243

Observations of superflares on solar-type stars indicate that they are associated with much larger starspots than appear on the Sun, occur more frequently on rapidly rotating stars and, contrary to a previous proposal, are not frequently associated with hot Jupiters. Solar flares are the most energetic explosions in the solar atmosphere, and similar flares occur on many stars. 'Superflares' many thousands of times more energetic than the average solar flare have been observed from a variety of stars, but the relatively small number observed on solar-type stars has hitherto precluded a detailed study of them. Now, on the basis of an analysis of data from the Kepler satellite, Maehara et al. report observations of 365 superflares, including 101 from slowly rotating solar-type stars, from a sample of around 83,000 stars observed over 120 days. The data suggest that superflares occur more frequently on rapidly rotating stars and on stars with 'starspots' much larger than the sunspots with which we are familiar. There is no historical record of superflares on the Sun in the past 2,000 years, and it is probable that none has occurred in the past one billion years. Bradley Schaefer discusses these findings in an accompanying News and Views, and concludes that it is extremely unlikely that the Sun will host a superflare. Solar flares are caused by the sudden release of magnetic energy stored near sunspots. They release 1029 to 1032 ergs of energy on a timescale of hours1. Similar flares have been observed on many stars, with larger ‘superflares’ seen on a variety of stars2,3, some of which are rapidly rotating4,5 and some of which are of ordinary solar type3,6. The small number of superflares observed on solar-type stars has hitherto precluded a detailed study of them. Here we report observations of 365 superflares, including some from slowly rotating solar-type stars, from about 83,000 stars observed over 120 days. Quasi-periodic brightness modulations observed in the solar-type stars suggest that they have much larger starspots than does the Sun. The maximum energy of the flare is not correlated with the stellar rotation period, but the data suggest that superflares occur more frequently on rapidly rotating stars. It has been proposed that hot Jupiters may be important in the generation of superflares on solar-type stars7, but none have been discovered around the stars that we have studied, indicating that hot Jupiters associated with superflares are rare.

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., Comment: 43 pages, 15 figures, accepted for publication in ApJ

We performed the high dispersion spectroscopy of solar-type superflare stars by Subaru/HDS, and estimate the stellar parameters and lithium abundance of the stars to compare with the Sun. Our spectroscopic analysis of superflare stars show more than half of targets have no evidence of binary system and the stellar parameters are in the range of solar-type stars (Notsu et al. 2015a&b). We also investigate the correlations of Lithium abundance with stellar atmospheric parameters, rotational velocity, and superflare activities to understand the nature of superflare stars and the possibility of the nucleosynthesis of lithium by superflares. The derived lithium abundance in superflare stars do not show the correlation with stellar parameters. As compared with the lithium abundance in Hyades cluster which is younger than the sun, it is suggested that half of observed stars are young. However, there are some objects which show the low lithium and slowly rotate from the estimated v sin(i) and period of brightness variation. These results indicate that the superflare stars are not only young stars but also old stars like our sun. In our observations, we could not find the any evidence of lithium productions by superflare.

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., Comment: 6 pages, 3 figures, to appear in the Proceedings of IAU Symposium No. 320 "Solar and Stellar Flares and Their Effects on Planets". The main results of this manuscript are based on Notsu et al. (2015a, PASJ, 67, 32) and (2015b, PASJ, 67, 33)

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., Comment: 19 pages, 6 figures, 2 tables, accepted for publication in PASJ

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, Comment: 45 pages, 15 figures, 5 tables, PASJ in press (Already published online on February 22, 2015. doi:10.1093/pasj/psv001 )

We present the results of an extensive survey of superflares on late-type stars (G, K, and M-type main sequence stars) using the Kepler satellite data. Wefound about 6,800 superflares on late-type stars from the data of about 120,000 stars observed over 500 days. The total bolometric energy of superflares in oursample ranges from 1032 erg to 1036 erg. Our data suggest that the occurrencefrequency of superflares depends on the surface temperature and the rotationperiod of stars. Superflares on M-type stars occur about 10-100 times morefrequently than those on G-type stars. Our results suggest that the average frequency ofsuperflares releasing 1034–1035 erg of energy (100-1,000 times larger than the largestsolar flares) on M-type stars and Sun-like stars is once in 10 years and once in a few thousand years respectively.

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., Comment: 11 pages, 6 figures, accepted for publication in PASJ

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., Comment: Accepted to the ApJ Supplement in August 2013

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., Comment: 33 pages, 7 figures, accepted by Publications of the Astronomical Society of Japan (PASJ) on July 18, 2013 (to be published on PASJ vol. 65, No. 5 (October 25, 2013))

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)., 33 pages, 9 figures, and 2 tables are contained in this paper. It was received by The Astrophysical Journal (ApJ) on December 30th, 2015, and was accepted on June 16th, 2016

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., Accepted by Publ. Astron. Soc. Japan on Dec. 6, 2012 (to be published on PASJ vol. 65, No. 3, (2013) June 25)

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., Comment: 25 pages, 9 figures, 1 tables, accepted for publication in PASJ

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., Comment: 34 pages, 10 figures, Published. (2013, ApJ, 771, 127)