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 optical observations of type Ia supernova (SN) 2019ein, starting at 2 days after the estimated explosion date. The spectra and the light curves show that SN 2019ein belongs to the High-Velocity (HV) and Bload Line groups with relatively rapid decline in the light curves (Delta m15(B) = 1.36 +- 0.02 mag) and the short rise time (15.37 +- 0.55 days). The Si II 6355 velocity, associated with a photospheric component but not with a detached high-velocity feature, reached ~ 20,000 km s-1 at 12 days before the B-band maximum. The line velocity however decreased very rapidly and smoothly toward the maximum light, where it was ~ 13,000 km s-1 as relatively low among HV SNe. This indicates that the speed of the spectral evolution of HV SNe Ia is correlated not only to the velocity at the maximum light, but also to the light curve decline rate like the case for Normal-Velocity (NV) SNe Ia. Spectral synthesis modeling shows that the outermost layer at > 17,000 km s-1 is well described by the O-Ne-C burning layer extending to at least 25,000 km s-1, and there is no unburnt carbon below 30,000 km s-1; these properties are largely consistent with the delayed detonation scenario, and are shared with the prototypical HV SN 2002bo despite the large difference in Delta m15(B). This structure is strikingly different from that derived for the well-studied NV SN 2011fe. We suggest that the relation between the mass of 56Ni (or Delta m15) and the extent of the O-Ne-C burning layer provides an important constraint on the explosion mechanism(s) of HV and NV SNe.

The coagulation of dust particles is a key process in planetesimal formation. However, the radial drift and bouncing barriers are not completely resolved, especially for silicate dust. Since the collision velocities of dust particles are regulated by turbulence in a protoplanetary disk, the turbulent clustering should be properly treated. To that end, direct numerical simulations (DNSs) of the Navier Stokes equations are requisite. In a series of papers, Pan & Padoan used a DNS with the Reynolds number Re~1000. Here, we perform DNSs with up to Re=16100, which allow us to track the motion of particles with Stokes numbers of 0.01<~St<~0.2 in the inertial range. By the DNSs, we confirm that the rms relative velocity of particle pairs is smaller by more than a factor of two, compared to those by Ormel & Cuzzi (2007). The distributions of the radial relative velocities are highly non-Gaussian. The results are almost consistent with those by Pan & Padoan or Pan et al. at low-Re. Also, we find that the sticking rates for equal-sized particles are much higher than those for different-sized particles. Even in the strong-turbulence case with alpha-viscosity of 10^{-2}, the sticking rates are as high as >~50% and the bouncing probabilities are as low as ~10% for equal-sized particles of St<~0.01. Thus, the turbulent clustering plays a significant role for the growth of cm-sized compact aggregates (pebbles) and also enhances the solid abundance, which may lead to the streaming instability in a disk.

We investigate cosmological consequences of an extended gravity model which belongs to the same class studied by Accetta and Steinhardt in an extended inflationary scenario. But we do not worry about inflation in our model; instead, we focus on a topological object formed during cosmological phase transitions. Although domain walls appear during first-order phase transitions such as QCD transition, they decay at the end of the phase transition. Therefore the "domain wall problem" does not exist in the suitable range of pameters and, on the contrary, the "fragments" of walls may become seeds of dark matter. A possible connection to "oscillating universe" model offered by Morikawa et al. is also discussed.

We construct the higher order terms of curvatures in Lagrangians of the scale factor for the Friedmann-Lemaitre-Robertson-Walker universe, which are linear in the second derivative of the scale factor with respect to cosmic time. It is shown that they are composed from the Lovelock tensors at the first step; iterative construction yields arbitrarily high order terms. The relation to the former work on higher order gravity is discussed. Despite the absence of scalar degrees of freedom in cosmological models which come from our Lagrangian, it is shown that an inflationary behavior of the scale factor can be found. The application to the thick brane solutions is also studied.

We consider the finite density effects of the fermion with U(1) gauge symmetry in Hosotani mechanism. We construct a vacuum gauge ball, a new kind of non-topological soliton, and investigate their properties numerically. We find the relations between the physical quantities.

We present the results of indirect searches for Weakly Interacting Massive Particles (WIMPs) with 1679.6 live days of data from the Super-Kamiokande detector using neutrino-induced upward through-going muons. The search is performed by looking for an excess of high energy muon neutrinos from WIMP annihilations in the Sun, the core of the Earth, and the Galactic Center, as compared to the number expected from the atmospheric neutrino background. No statistically significant excess was seen. We calculate flux limits in various angular cones around each of the above celestial objects. We obtain conservative model-independent upper limits on WIMP-nucleon cross-section as a function of WIMP mass and compare these results with the corresponding results from direct dark matter detection experiments.

Beginning with the Pauli-Fierz theory, we construct a model for multi-graviton theory. Couplings between gravitons belonging to nearest-neighbor ``theory spaces'' lead to a discrete mass spectrum. Our model coincides with the Kaluza-Klein theory whose fifth dimension is latticized. We evaluate one-loop vacuum energy in models with a circular latticized extra dimension as well as with compact continuous dimensions. We find that the vacuum energy can take a positive value, if the dimension of the continuous space time is $6, 10,...$. Moreover, since the amount of the vacuum energy can be an arbitrary small value according to the choice of parameters in the model, our models is useful to explain the small positive dark energy in the present universe.

An exact solution for a nonrotating boson star in (2+1) dimensional gravity with a negative cosmological constant is found. The relations among mass, particle number, and radius of the (2+1) dimensional boson star are shown.