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.

Observations of transient objects, such as short gamma-ray bursts (GRBs) and electromagnetic counterparts of gravitational wave sources, require prompt spectroscopy. To carry out prompt spectroscopy, we developed an optical-fiber integral field unit (IFU) and connected it with an existing optical spectrograph KOOLS. KOOLS with IFU, hereafter KOOLS--IFU, was mounted on the Okayama Astrophysical Observatory 188-cm telescope. The fiber core and cladding diameters of the fiber bundle are 100 $\mu$m and 125 $\mu$m, respectively, and 127 fibers are hexagonally close-packed in the sleeve of the two-dimensional (2D) fiber array. We conducted test observations to measure the KOOLS--IFU performance and obtained the following conclusions: (1) the spatial sampling is 2''.34 $\pm$ 0''.05 per fiber, and the total field of view (FoV) is 30''.4 $\pm$ 0''.65 with 127 fibers; (2) the observable wavelength and the spectral resolving power of the grisms of KOOLS are 4030--7310 A and 400--600; 5020--8830 A and 600--900; 4160--6000 A and 1000--1200; and 6150--7930 A and 1800--2400, respectively; and (3) the estimated limiting magnitude is 18.2--18.7 AB mag during 30-min exposure under the optimal condition.

We report the initial result of an adaptive-optics assisted, optical integral-field-unit observation on IRAS04576+0912, the nearest (z=0.039) active galactic nucleus with a prominent blueshift/tail in [O III] emission from a sample of such objects that we have collected from the literature. We aim at addressing the putative quasar-mode feedback process with Subaru/Kyoto 3D II+AO188. The optical waveband (6400--7500 AA) enables us to measure the gas density via the [S II] doublets, in contrast to earlier Near-IR studies. Since the fast [O III] outflow happens only around rapidly growing central black holes, this object is suitable for investigating the black hole-galaxy coevolution. The obtained data cube exhibits blue tail in the [S II] emission at many lenslets. By fitting the spectrum with the high excess flux at the [S II] blue tail, we find the fast (~ 860 km/s), dense (>3000/cc), wide-angle and offset outflow in central 100-pc scales. Although the large opening angle and the high gas outflow-to-accretion ratio may favour the feedback hypothesis, the inferred kinetic power injection rate of this ionized gas outflow seems insufficient to influence the whole host galaxy. A conventional assumption of a low density must have overestimated the feedback process.

We perform a $z$-band survey for an optical counterpart of a binary neutron star coalescence GW170817 with Subaru/Hyper Suprime-Cam. Our untargeted transient search covers $23.6$ deg$^2$ corresponding to the $56.6\%$ credible region of GW170817 and reaches the $50\%$ completeness magnitude of $20.6$ mag on average. As a result, we find 60 candidates of extragalactic transients, including J-GEM17btc (a.k.a. SSS17a/DLT17ck). While J-GEM17btc is associated with NGC 4993 that is firmly located inside the 3D skymap of GW170817, the other 59 candidates do not have distance information in the GLADE v2 catalog or NASA/IPAC Extragalactic Database (NED). Among 59 candidates, 58 are located at the center of extended objects in the Pan-STARRS1 catalog, while one candidate has an offset. We present location, $z$-band apparent magnitude, and time variability of the candidates and evaluate the probabilities that they are located inside of the 3D skymap of GW170817. The probability for J-GEM17btc is $64\%$ being much higher than those for the other 59 candidates ($9.3\times10^{-3}-2.1\times10^{-1}\%$). Furthermore, the possibility, that at least one of the other 59 candidates is located within the 3D skymap, is only $3.2\%$. Therefore, we conclude that J-GEM17btc is the most-likely and distinguished candidate as the optical counterpart of GW170817.

The first detected gravitational wave from a neutron star merger was GW170817. In this study, we present J-GEM follow-up observations of SSS17a, an electromagnetic counterpart of GW170817. SSS17a shows a 2.5-mag decline in the $z$-band from 1.7 days to 7.7 days after the merger. Such a rapid decline is not comparable with supernovae light curves at any epoch. The color of SSS17a also evolves rapidly and becomes redder for later epochs; the $z-H$ color changed by approximately 2.5 mag in the period of 0.7 days to 7.7 days. The rapid evolution of both the optical brightness and the color are consistent with the expected properties of a kilonova that is powered by the radioactive decay of newly synthesized $r$-process nuclei. Kilonova models with Lanthanide elements can reproduce the aforementioned observed properties well, which suggests that $r$-process nucleosynthesis beyond the second peak takes place in SSS17a. However, the absolute magnitude of SSS17a is brighter than the expected brightness of the kilonova models with the ejecta mass of 0.01 $\Msun$, which suggests a more intense mass ejection ($\sim 0.03 \Msun$) or possibly an additional energy source.

We present an 8 pc x 5 pc resolution view of the central ~ 200 pc region of the nearby starburst galaxy NGC 253, based on ALMA Band 7 ({\lambda} ~ 0.85 mm or {\nu} ~ 350 GHz) observations covering 11 GHz. We resolve the nuclear starburst of NGC 253 into eight dusty star-forming clumps, 10 pc in scale, for the first time. These clumps, each of which contains (4-10) x 10^4 M_sun of dust (assuming that the dust temperature is 25 K) and up to 6 x 10^2 massive (O5V) stars, appear to be aligned in two parallel ridges, while they have been blended in previous studies. Despite the similarities in sizes and dust masses of these clumps, their line spectra vary drastically from clump to clump although they are separated by only ~ 10 pc. Specifically, one of the clumps, Clump 1, exhibits line confusion-limited spectra with at least 36 emission lines from 19 molecules (including CH3OH, HNCO, H2CO, CH3CCH, H2CS, and H3O+) and a hydrogen recombination line (H26{\alpha}), while much fewer kinds of molecular lines are detected in some other clumps where fragile species, such as complex organic molecules and HNCO, completely disappear from their spectra. We demonstrate the existence of hot molecular gas ($T_rot (SO2) = 90 {\pm} 11 K) in the former clump, which suggests that the hot and chemically rich environments are localized within a 10-pc scale star-forming clump.

We present the results of the detailed analysis of an optical imaging survey conducted using the Subaru / Hyper Suprime-Cam (HSC), which aims to identify an optical counterpart to the gravitational wave event GW151226. In half a night, the $i$- and $z$-band imaging survey by HSC covers 63.5deg$^2$ of the error region, which contains about 7\% of the LIGO localization probability, and the same field is observed in three different epochs. The detectable magnitude of the candidates in a differenced image is evaluated as $i \sim 23.2$ mag for the requirement of at least two 5$\sigma$ detections, and 1744 candidates are discovered. Assuming a kilonova as an optical counterpart, we compared the optical properties of the candidates with model predictions. A red and rapidly declining light curve condition enables the discrimination of a kilonova from other transients, and a small number of candidates satisfy this condition. The presence of stellar-like counterparts in the reference frame suggests that the surviving candidates are likely to be flare stars. The fact that most of those candidates are in galactic plane, $|b|<5^{\circ}$, supports this interpretation. We also checked whether the candidates are associated with the nearby GLADE galaxies, which reduces the number of contaminants even with a looser color cut. When a better probability map (with localization accuracy of $\sim50{\rm deg}^2$) is available, kilonova searches of up to approximately $200$ Mpc will become feasible by conducting immediate follow-up observations with an interval of 3--6 days.

We present the H$\alpha$ intensity map of the host galaxy of the repeating fast radio burst FRB 121102 at a redshift of z=0.193 obtained with the AO-assisted Kyoto 3DII optical integral-field unit mounted on the 8.2-m Subaru Telescope. We detected a compact H$\alpha$-emitting (i.e., star-forming) region in the galaxy, which has a much smaller angular size [$< 0".57$ (1.9 kpc) at full width at half maximum (FWHM)] than the extended stellar continuum emission region determined by the Gemini/GMOS z'-band image [$\simeq 1".4$ (4.6 kpc) at FWHM with ellipticity b/a=0.45]. The spatial offset between the centroid of the H$\alpha$ emission region and the position of the radio bursts is $0".08 \pm 0".02$ ($0.26 \pm 0.07$ kpc), indicating that FRB 121102 is located within the star-forming region. This close spatial association of FRB 121102 with the star-forming region is consistent with expectations from young pulsar/magnetar models for FRB 121102, and it also suggests that the observed H$\alpha$ emission region can make a major dispersion measure (DM) contribution to the host galaxy DM component of FRB 121102. Nevertheless, the largest possible value of the DM contribution from the H$\alpha$ emission region inferred from our observations still requires a significant amount of ionized baryons in intergalactic medium (the so-called `missing' baryons) as the DM source of FRB 121102, and we obtain a 90\% confidence level lower limit on the cosmic baryon density in the intergalactic medium in the low-redshift universe as $\Omega_{IGM} > 0.012$.

We report the results of optical--infrared follow-up observations of the gravitational wave (GW) event GW151226 detected by the Advanced LIGO in the framework of J-GEM (Japanese collaboration for Gravitational wave ElectroMagnetic follow-up). We performed wide-field optical imaging surveys with Kiso Wide Field Camera (KWFC), Hyper Suprime-Cam (HSC), and MOA-cam3. The KWFC survey started at 2.26 days after the GW event and covered 778 deg$^2$ centered at the high Galactic region of the skymap of GW151226. We started the HSC follow-up observations from 12 days after the event and covered an area of 63.5 deg$^2$ of the highest probability region of the northern sky with the limiting magnitudes of 24.6 and 23.8 for i band and z band, respectively. MOA-cam3 covered 145 deg$^2$ of the skymap with MOA-red filter 2.5 months after the GW alert. Total area covered by the wide-field surveys was 986.5 deg$^2$. The integrated detection probability of all the observed area was $\sim$29%. We also performed galaxy-targeted observations with six optical and near-infrared telescopes from 1.61 days after the event. Total of 238 nearby (<100 Mpc) galaxies were observed with the typical I band limiting magnitude of $\sim$19.5. We detected 13 supernova candidates with the KWFC survey, and 60 extragalactic transients with the HSC survey. Two third of the HSC transients were likely supernovae and the remaining one third were possible active galactic nuclei. With our observational campaign, we found no transients that are likely to be associated with GW151226.

We present our optical follow-up observations to search for an electromagnetic counterpart of the first gravitational wave source GW150914 in the framework of the Japanese collaboration for Gravitational wave ElectroMagnetic follow-up (J-GEM), which is an observing group utilizing optical and radio telescopes in Japan, as well as those in New Zealand, China, South Africa, Chile, and Hawaii. We carried out a wide-field imaging survey with Kiso Wide Field Camera (KWFC) on the 1.05-m Kiso Schmidt telescope in Japan and a galaxy-targeted survey with Tripole5 on the B&C 61-cm telescope in New Zealand. Approximately 24 deg2 regions in total were surveyed in i-band with KWFC and 18 nearby galaxies were observed with Tripole5 in g-, r-, and i-bands 4-12 days after the gravitational wave detection. Median 5-sigma depths are i~18.9 mag for the KWFC data and g~18.9 mag, r~18.7 mag, and i~18.3 mag for the Tripole5 data. Probability for a counterpart to be in the observed area is 1.2% in the initial skymap and 0.1% in the final skymap. We do not find any transient source associated to an external galaxy with spatial offset from its center, which is consistent with the local supernova rate. We summarize future prospects and ongoing efforts to pin down electromagnetic counterparts of binary black hole mergers as well as neutron star mergers.

We investigated the fraction of [OII] emitters in galaxies at z~0.9 as a function of the local galaxy density in the Hubble Space Telescope (HST) COSMOS 2 square degree field. [OII] emitters are selected by the narrow-band excess technique with the NB711-band imaging data taken with Suprime-Cam on the Subaru telescope. We carefully selected 614 photo-z selected galaxies with M_U3500 < -19.31 at z=0.901-0.920, which includes 195 [OII] emitters, to directly compare results with our previous study at z~1.2. We found that the fraction is almost constant at 0.3 Mpc^-2 < Sigma_10th < 10 Mpc^-2. We also checked the fraction of galaxies with blue rest-frame colors of NUV-R < 2 in our photo-z selected sample, and found that the fraction of blue galaxies does not significantly depend on the local density. On the other hand, the semi-analytic model of galaxy formation predicted that the fraction of star-forming galaxies at z~0.9 decreases with increasing the projected galaxy density even if the effects of the projection and the photo-z error in our analysis were taken into account. The fraction of [OII] emitters decreases from ~60% at z~1.2 to ~30% at z~0.9 independent of the galaxy environment. The decrease of the [OII] emitter fraction could be explained mainly by the rapid decrease of the star formation activity in the universe from z~1.2 to z~0.9.

The M82 `cap' is a gas cloud at a projected radius of 11.6 kpc along the minor axis of this well known superwind source. The cap has been detected in optical line emission and X-ray emission and therefore provides an important probe of the wind energetics. In order to investigate the ionization source of the cap, we observed it with the Kyoto3DII Fabry-Perot instrument mounted on the Subaru Telescope. Deep continuum, Ha, [NII]6583/Ha, and [SII]6716,6731/Ha maps were obtained with sub-arcsecond resolution. The superior spatial resolution compared to earlier studies reveals a number of bright Ha emitting clouds within the cap. The emission line widths (< 100 km s^-1 FWHM) and line ratios in the newly identified knots are most reasonably explained by slow to moderate shocks velocities (v_shock = 40--80 km s^-1) driven by a fast wind into dense clouds. The momentum input from the M82 nuclear starburst region is enough to produce the observed shock. Consequently, earlier claims of photoionization by the central starburst are ruled out because they cannot explain the observed fluxes of the densest knots unless the UV escape fraction is very high (f_esc > 60%), i.e., an order of magnitude higher than observed in dwarf galaxies to date. Using these results, we discuss the evolutionary history of the M82 superwind. Future UV/X-ray surveys are expected to confirm that the temperature of the gas is consistent with our moderate shock model.

We present our new deep optical imaging and long-slit spectroscopy for Arp 220 that is the archetypical ULIRG in the local universe. Our sensitive Ha imaging has newly revealed large-scale, Ha absorption, i.e., post-starburst regions in this merger; one is found in the eastern superbubble and the other is in the two tidal tails that are clearly reveled in our deep optical imaging. The size of Ha absorption region in the eastern bubble is 5 kpc x 7.5 kpc and the observed Ha equivalent widths are ~2 A +- 0.2 A. The sizes of the northern and southern Ha-absorption tidal tails are ~5 kpc x 10 kpc and ~6 kpc x 20 kpc, respectively. The observed Ha equivalent widths range from 4 A to 7 A. In order to explain the presence of the two post-starburst tails, we suggest a possible multiple-merger scenario for Arp 220 in which two post-starburst disk-like structures merged into one, and then caused the two tails. This favors that Arp 220 is a multiple merging system composed of four or more galaxies, arising from a compact group of galaxies. Taking our new results into account, we discuss a star formation history in the last 1 Gyr in Arp 220.

In order to investigate the evolution of E+A galaxies, we observed a galaxy SDSS J160241.00+521426.9, a possible E+A progenitor which shows both emission and strong Balmer absorptions, and its neighbor galaxy. We used the integral field spectroscopic mode of the Kyoto Tridimensional Spectrograph (Kyoto3DII), mounted on the University of Hawaii 88-inch telescope located on Mauna Kea, and the slit-spectroscopic mode of the Faint Object Camera and Spectrograph (FOCAS) on the Subaru Telescope. We found a strong Balmer absorption region in the center of the galaxy and an emission-line region located 2 kpc from the center, in the direction of its neighbor galaxy. The recession velocities of the galaxy and its neighbor galaxy differ only by 100 km s^-1, which suggests that they are a physical pair and would have been interacting. Comparing observed Lick indices of Balmer lines and color indices with those predicted from stellar population synthesis models, we find that a suddenly quenched star-formation scenario is plausible for the star-formation history of the central region. We consider that star formation started in the galaxy due to galaxy interactions and was quenched in the central region, whereas star formation in a region offset from the center still continues or has begun recently. This work is the first study of a possible E+A progenitor using spatially resolved spectroscopy.

We have observed the central region of the nearby starburst galaxy NGC 253 with the Kyoto Tridimensional Spectrograph II (Kyoto3DII) Fabry-Perot mode in order to investigate the properties of its galactic wind. Since this galaxy has a large inclination, it is easy to observe its galactic wind. We produced the Ha, [N II]6583, and [S II]6716,6731 images, as well as those line ratio maps. The [N II]/Ha ratio in the galactic wind region is larger than those in H II regions in the galactic disk. The [N II]/Ha ratio in the southeastern filament, a part of the galactic wind, is the largest and reaches about 1.5. These large [N II]/Ha ratios are explained by shock ionization/excitation. Using the [S II]/Ha ratio map, we spatially separate the galactic wind region from the starburst region. The kinetic energy of the galactic wind can be sufficiently supplied by supernovae in a starburst region in the galactic center. The shape of the galactic wind and the line ratio maps are non-axisymmetric about the galactic minor axis, which is also seen in M82. In the [N II]6583/[S II]6716,6731 map, the positions with large ratios coincide with the positions of star clusters found in the Hubble Space Telescope (HST) observation. This means that intense star formation causes strong nitrogen enrichment in these regions. Our unique data of the line ratio maps including [S II] lines have demonstrated their effectiveness for clearly distinguishing between shocked gas regions and starburst regions, determining the extent of galactic wind and its mass and kinetic energy, and discovering regions with enhanced nitrogen abundance.