Haumea is the only known trans-Neptunian object with a ring. The ring is in a position that causes a 3:1 spin-orbit resonance with the rotational period of Haumea, which has a triaxial shape. The non-axisymmetric gravitational field around Haumea is thought to affect the dynamics of the ring; however, the process of ring formation has not been elucidated. In this study, we analyze in some detail a potential ring formation scenario for Haumea. We first calculated the gravitational field around the triaxial ellipsoid and estimated the distance at which an object revolving around Haumea can exist stably using simulation that incorporated the time-varying gravitational field. The results of this simulation showed that the trajectory of the object became unstable just inside its current ring position. Next, we analytically derived the Roche radius for a rigid body revolving around a triaxial ellipsoid and showed that the Roche radius could be near the current ring position. Furthermore, as a parameter study, we performed N-body simulations using the coefficient of rubble pile restitution as a variable. Results demonstrated that, according to numerous parameters, the position of the Roche radius was near the current position of Haumea's ring. Based on these findings, we can assume that there is a high possibility that the ring formed in the region between the boundary of the unstable region of the orbit and the Roche radius. The scenario presented in this study could help explain the process by which Haumea's ring formed.

The ice-giant planet Uranus likely underwent a giant impact, given that its spin axis is tilted by 98 degrees. That its satellite system is equally inclined and prograde suggests that it was formed as a consequence of the impact. However, the disks predicted by the impact simulations generally have sizes one order smaller and masses two orders larger than those of the observed system at present. Here we show, by means of a theoretical model, that the Uranian satellite formation is regulated by the evolution of the impact-generated disk. Because the vaporization temperature of water ice is low and both Uranus and the impactor are assumed to be ice-dominated, we can conclude that the impact-generated disk has mostly vaporized. We predict that the disk lost a significant amount of water vapour mass and spread to the levels of the current system until the disk cooled down enough for ice condensation and accretion of icy particles to begin. From the predicted distribution of condensed ices, our N-body simulation is able to reproduce the observed mass-orbit configuration of Uranian satellites. This scenario contrasts with the giant-impact model for the Earth's Moon, in which about half of the compact, impact-generated, solid or liquid disk is immediately incorporated into the Moon on impact.

It is generally accepted that the four major (Galilean) satellites formed out of the gas disk that accompanied Jupiter's formation. However, understanding the specifics of the formation process is challenging as both small particles (pebbles) as well as the satellites are subject to fast migration processes. Here, we hypothesize a new scenario for the origin of the Galilean system, based on the capture of several planetesimal seeds and subsequent slow accretion of pebbles. To halt migration, we invoke an inner disk truncation radius, and other parameters are tuned for the model to match physical, dynamical, compositional, and structural constraints. In our scenario it is natural that Ganymede's mass is determined by pebble isolation. Our slow-pebble-accretion scenario then reproduces the following characteristics: (1) the mass of all the Galilean satellites; (2) the orbits of Io, Europa, and Ganymede captured in mutual 2:1 mean motion resonances; (3) the ice mass fractions of all the Galilean satellites; (4) the unique ice-rock partially differentiated Callisto and the complete differentiation of the other satellites. Our scenario is unique to simultaneously reproduce these disparate properties.

Hydrodynamic simulations of a giant impact to proto-Uranus indicated that such an impact could tilt its rotational axis and produce a circumplanetary debris disk beyond the corotation radius of Uranus. However, whether Uranian satellites can actually be formed from such a wide disk remains unclear. Herein, we modeled a wide debris disk of solids with several initial conditions inferred from the hydrodynamic simulations, and performed $N$-body simulations to investigate in-situ satellite formation from the debris disk. We also took account of orbital evolutions of satellites due to the planetary tides after the growth of satellites. We found that, in any case, the orbital distribution of the five major satellites could not be reproduced from the disk as long as the power index of its surface density is similar to that of the disk generated just after the giant impact: Satellites in the middle region obtained much larger masses than Ariel or Umbriel, while the outermost satellite did not grow to the mass of Oberon. Our results indicate that we should consider the thermal and viscous evolution of the evaporated disk after the giant impact to form the five major satellites through the in-situ formation scenario. On the other hand, the small inner satellites would be formed from the rings produced by the disrupted satellites which migrated from around the corotation radius of Uranus due to the planetary tides.

Sub-Jupiter classed circumbinary planets discovered in close-in binary systems have orbits just beyond the dynamically unstable region, which is determined by the eccentricity and mass ratio of the host binary stars. These planets are assumed to have formed beyond the snow line and migrated to the current orbits rather than forming in situ. We propose a scenario in which a planet formed beyond the snow line and migrated to the inner edge of the circumbinary disk, which was within the unstable area, and then moved to the current orbit through outward transportation. This outward transportation is driven by the balance of orbital excitation of the central stars inside the gravitationally unstable region and damping by the gas-drag force. We carried out N-body simulations with a dissipating circumbinary protoplanetary disk for binary systems with different eccentricities and mass ratios. Planets are more likely to achieve a stable orbit just beyond the unstable region in less eccentric binary systems. This result is not as sensitive to mass ratio as it is to eccentricity. These dependencies are consistent with the data from observed binary systems hosting circumbinary planets. We find CBPs' orbits close to the instability boundaries are explained by our orbital evolution scenario.

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 present the results from the Hitomi Soft Gamma-ray Detector (SGD) observation of the Crab nebula. The main part of SGD is a Compton camera, which in addition to being a spectrometer, is capable of measuring polarization of gamma-ray photons. The Crab nebula is one of the brightest X-ray / gamma-ray sources on the sky, and, the only source from which polarized X-ray photons have been detected. SGD observed the Crab nebula during the initial test observation phase of Hitomi. We performed the data analysis of the SGD observation, the SGD background estimation and the SGD Monte Carlo simulations, and, successfully detected polarized gamma-ray emission from the Crab nebula with only about 5 ks exposure time. The obtained polarization fraction of the phase-integrated Crab emission (sum of pulsar and nebula emissions) is (22.1 $\pm$ 10.6)% and, the polarization angle is 110.7$^o$ + 13.2 / $-$13.0$^o$ in the energy range of 60--160 keV (The errors correspond to the 1 sigma deviation). The confidence level of the polarization detection was 99.3%. The polarization angle measured by SGD is about one sigma deviation with the projected spin axis of the pulsar, 124.0$^o$ $\pm$0.1$^o$.

The future detection projects of gravitational waves are expected to have the sensitivity of detecting the primordial gravitational waves, which may be useful to get new insight into the very early universe. It is essential to analyze the evolution equation of the gravitational waves to estimate the present field strength of the primordial gravitational waves. In this paper, we obtain analytic solutions of the gravitational wave equation in the presence of non-relativistic matters as well as the cosmological constant. Although it is difficult to obtain the solution directly, we find that the equation for the square of the amplitudes has a simple polynomial solution. This quantity, which is directly related to the energy density of the gravitational waves, turns out to be useful to construct analytic solutions for the amplitudes by using Weierstrass's elliptic functions.

The formation of the Moon from the circumterrestrial disk has been investigated by using $N$-body simulations with the number $N$ of particles limited from $10^4$ to $10^5$. We develop an $N$-body simulation code on multiple Pezy-SC processors and deploy FDPS (Framework for Developing Particle Simulators) to deal with large number of particles. We execute several high- and extra-high-resolution $N$-body simulations of lunar accretion from a circumterrestrial disk of debris generated by a giant impact on Earth. The number of particles is up to $10^7$, in which 1 particle corresponds to a 10 km-size satellitesimal. We find that the spiral structures inside the Roche limit radius differ between low-resolution simulations ($N \leq10^5$) and high-resolution simulations ($N \geq10^6$). According to this difference, angular momentum fluxes, which determine the accretion timescale of the Moon also depend on the numerical resolution.

We present results from the Hitomi X-ray observation of a young composite-type supernova remnant (SNR) G21.5$-$0.9, whose emission is dominated by the pulsar wind nebula (PWN) contribution. The X-ray spectra in the 0.8-80 keV range obtained with the Soft X-ray Spectrometer (SXS), Soft X-ray Imager (SXI) and Hard X-ray Imager (HXI) show a significant break in the continuum as previously found with the NuSTAR observation. After taking into account all known emissions from the SNR other than the PWN itself, we find that the Hitomi spectra can be fitted with a broken power law with photon indices of $\Gamma_1=1.74\pm0.02$ and $\Gamma_2=2.14\pm0.01$ below and above the break at $7.1\pm0.3$ keV, which is significantly lower than the NuSTAR result ($\sim9.0$ keV). The spectral break cannot be reproduced by time-dependent particle injection one-zone spectral energy distribution models, which strongly indicates that a more complex emission model is needed, as suggested by recent theoretical models. We also search for narrow emission or absorption lines with the SXS, and perform a timing analysis of PSR J1833$-$1034 with the HXI and SGD. No significant pulsation is found from the pulsar. However, unexpectedly, narrow absorption line features are detected in the SXS data at 4.2345 keV and 9.296 keV with a significance of 3.65 $\sigma$. While the origin of these features is not understood, their mere detection opens up a new field of research and was only possible with the high resolution, sensitivity and ability to measure extended sources provided by an X-ray microcalorimeter.

The present paper investigates the temperature structure of the X-ray emitting plasma in the core of the Perseus cluster using the 1.8--20.0 keV data obtained with the Soft X-ray Spectrometer (SXS) onboard the Hitomi Observatory. A series of four observations were carried out, with a total effective exposure time of 338 ks and covering a central region $\sim7'$ in diameter. The SXS was operated with an energy resolution of $\sim$5 eV (full width at half maximum) at 5.9 keV. Not only fine structures of K-shell lines in He-like ions but also transitions from higher principal quantum numbers are clearly resolved from Si through Fe. This enables us to perform temperature diagnostics using the line ratios of Si, S, Ar, Ca, and Fe, and to provide the first direct measurement of the excitation temperature and ionization temperature in the Perseus cluster. The observed spectrum is roughly reproduced by a single temperature thermal plasma model in collisional ionization equilibrium, but detailed line ratio diagnostics reveal slight deviations from this approximation. In particular, the data exhibit an apparent trend of increasing ionization temperature with increasing atomic mass, as well as small differences between the ionization and excitation temperatures for Fe, the only element for which both temperatures can be measured. The best-fit two-temperature models suggest a combination of 3 and 5 keV gas, which is consistent with the idea that the observed small deviations from a single temperature approximation are due to the effects of projection of the known radial temperature gradient in the cluster core along the line of sight. Comparison with the Chandra/ACIS and the XMM-Newton/RGS results on the other hand suggests that additional lower-temperature components are present in the ICM but not detectable by Hitomi SXS given its 1.8--20 keV energy band.

The Hitomi SXS spectrum of the Perseus cluster, with $\sim$5 eV resolution in the 2-9 keV band, offers an unprecedented benchmark of the atomic modeling and database for hot collisional plasmas. It reveals both successes and challenges of the current atomic codes. The latest versions of AtomDB/APEC (3.0.8), SPEX (3.03.00), and CHIANTI (8.0) all provide reasonable fits to the broad-band spectrum, and are in close agreement on best-fit temperature, emission measure, and abundances of a few elements such as Ni. For the Fe abundance, the APEC and SPEX measurements differ by 16%, which is 17 times higher than the statistical uncertainty. This is mostly attributed to the differences in adopted collisional excitation and dielectronic recombination rates of the strongest emission lines. We further investigate and compare the sensitivity of the derived physical parameters to the astrophysical source modeling and instrumental effects. The Hitomi results show that an accurate atomic code is as important as the astrophysical modeling and instrumental calibration aspects. Substantial updates of atomic databases and targeted laboratory measurements are needed to get the current codes ready for the data from the next Hitomi-level mission.

We present Hitomi observations of N132D, a young, X-ray bright, O-rich core-collapse supernova remnant in the Large Magellanic Cloud (LMC). Despite a very short observation of only 3.7 ks, the Soft X-ray Spectrometer (SXS) easily detects the line complexes of highly ionized S K and Fe K with 16-17 counts in each. The Fe feature is measured for the first time at high spectral resolution. Based on the plausible assumption that the Fe K emission is dominated by He-like ions, we find that the material responsible for this Fe emission is highly redshifted at ~800 km/s compared to the local LMC interstellar medium (ISM), with a 90% credible interval of 50-1500 km/s if a weakly informative prior is placed on possible line broadening. This indicates (1) that the Fe emission arises from the supernova ejecta, and (2) that these ejecta are highly asymmetric, since no blue-shifted component is found. The S K velocity is consistent with the local LMC ISM, and is likely from swept-up ISM material. These results are consistent with spatial mapping that shows the He-like Fe concentrated in the interior of the remnant and the S tracing the outer shell. The results also show that even with a very small number of counts, direct velocity measurements from Doppler-shifted lines detected in extended objects like supernova remnants are now possible. Thanks to the very low SXS background of ~1 event per spectral resolution element per 100 ks, such results are obtainable during short pointed or slew observations with similar instruments. This highlights the power of high-spectral-resolution imaging observations, and demonstrates the new window that has been opened with Hitomi and will be greatly widened with future missions such as the X-ray Astronomy Recovery Mission (XARM) and Athena.

We report a Hitomi observation of IGR J16318-4848, a high-mass X-ray binary system with an extremely strong absorption of N_H~10^{24} cm^{-2}. Previous X-ray studies revealed that its spectrum is dominated by strong fluorescence lines of Fe as well as continuum emission. For physical and geometrical insight into the nature of the reprocessing material, we utilize the high spectroscopic resolving power of the X-ray microcalorimeter (the soft X-ray spectrometer; SXS) and the wide-band sensitivity by the soft and hard X-ray imager (SXI and HXI) aboard Hitomi. Even though photon counts are limited due to unintended off-axis pointing, the SXS spectrum resolves Fe K{\alpha_1} and K{\alpha_2} lines and puts strong constraints on the line centroid and width. The line width corresponds to the velocity of 160^{+300}_{-70} km s^{-1}. This represents the most accurate, and smallest, width measurement of this line made so far from any X-ray binary, much less than the Doppler broadening and shift expected from speeds which are characteristic of similar systems. Combined with the K-shell edge energy measured by the SXI and HXI spectra, the ionization state of Fe is estimated to be in the range of Fe I--IV. Considering the estimated ionization parameter and the distance between the X-ray source and the absorber, the density and thickness of the materials are estimated. The extraordinarily strong absorption and the absence of a Compton shoulder component is confirmed. These characteristics suggest reprocessing materials which are distributed in a narrow solid angle or scattering primarily with warm free electrons or neutral hydrogen.

The origin of the narrow Fe-K{\alpha} fluorescence line at 6.4 keV from active galactic nuclei has long been under debate; some of the possible sites are the outer accretion disk, the broad line region, a molecular torus, or interstellar/intracluster media. In February-March 2016, we performed the first X-ray microcalorimeter spectroscopy with the Soft X-ray Spectrometer (SXS) onboard the Hitomi satellite of the Fanaroff-Riley type I radio galaxy NGC 1275 at the center of the Perseus cluster of galaxies. With the high energy resolution of ~5 eV at 6 keV achieved by Hitomi/SXS, we detected the Fe-K{\alpha} line with ~5.4 {\sigma} significance. The velocity width is constrained to be 500-1600 km s$^{-1}$ (FWHM for Gaussian models) at 90% confidence. The SXS also constrains the continuum level from the NGC 1275 nucleus up to ~20 keV, giving an equivalent width ~20 eV of the 6.4 keV line. Because the velocity width is narrower than that of broad H{\alpha} line of ~2750 km s$^{-1}$, we can exclude a large contribution to the line flux from the accretion disk and the broad line region. Furthermore, we performed pixel map analyses on the Hitomi/SXS data and image analyses on the Chandra archival data, and revealed that the Fe-K{\alpha} line comes from a region within ~1.6 kpc from the NGC 1275 core, where an active galactic nucleus emission dominates, rather than that from intracluster media. Therefore, we suggest that the source of the Fe-K{\alpha} line from NGC 1275 is likely a low-covering fraction molecular torus or a rotating molecular disk which probably extends from a pc to hundreds pc scale in the active galactic nucleus system.

Extending the earlier measurements reported in Hitomi collaboration (2016, Nature, 535, 117), we examine the atmospheric gas motions within the central 100~kpc of the Perseus cluster using observations obtained with the Hitomi satellite. After correcting for the point spread function of the telescope and using optically thin emission lines, we find that the line-of-sight velocity dispersion of the hot gas is remarkably low and mostly uniform. The velocity dispersion reaches maxima of approximately 200~km~s$^{-1}$ toward the central active galactic nucleus (AGN) and toward the AGN inflated north-western `ghost' bubble. Elsewhere within the observed region, the velocity dispersion appears constant around 100~km~s$^{-1}$. We also detect a velocity gradient with a 100~km~s$^{-1}$ amplitude across the cluster core, consistent with large-scale sloshing of the core gas. If the observed gas motions are isotropic, the kinetic pressure support is less than 10\% of the thermal pressure support in the cluster core. The well-resolved optically thin emission lines have Gaussian shapes, indicating that the turbulent driving scale is likely below 100~kpc, which is consistent with the size of the AGN jet inflated bubbles. We also report the first measurement of the ion temperature in the intracluster medium, which we find to be consistent with the electron temperature. In addition, we present a new measurement of the redshift to the brightest cluster galaxy NGC~1275.

We present a joint X-ray, optical and weak-lensing analysis for X-ray luminous galaxy clusters selected from the MCXC (Meta-Catalog of X-Ray Detected Clusters of Galaxies) cluster catalog in the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) survey field with S16A data, As a pilot study of our planned series papers, we measure hydrostatic equilibrium (H.E.) masses using XMM-Newton data for four clusters in the current coverage area out of a sample of 22 MCXC clusters. We additionally analyze a non-MCXC cluster associated with one MCXC cluster. We show that H.E. masses for the MCXC clusters are correlated with cluster richness from the CAMIRA catalog (Oguri et al. 2017), while that for the non-MCXC cluster deviates from the scaling relation. The mass normalization of the relationship between the cluster richness and H.E. mass is compatible with one inferred by matching CAMIRA cluster abundance with a theoretical halo mass function. The mean gas mass fraction based on H.E. masses for the MCXC clusters is $\langle f_{\rm gas} \rangle = 0.125\pm0.012$ at spherical overdensity $\Delta=500$, which is $\sim80-90$ percent of the cosmic mean baryon fraction, $\Omega_b/\Omega_m$, measured by cosmic microwave background experiments. We find that the mean baryon fraction estimated from X-ray and HSC-SSP optical data is comparable to $\Omega_b/\Omega_m$. A weak-lensing shear catalog of background galaxies, combined with photometric redshifts, is currently available only for three clusters in our sample. Hydrostatic equilibrium masses roughly agree with weak-lensing masses, albeit with large uncertainty. This study demonstrates that further multiwavelength study for a large sample of clusters using X-ray, HSC-SSP optical and weak lensing data will enable us to understand cluster physics and utilize cluster-based cosmology.

Thanks to its high spectral resolution (~5 eV at 6 keV), the Soft X-ray Spectrometer (SXS) on board Hitomi enables us to measure the detailed structure of spatially resolved emission lines from highly ionized ions in galaxy clusters for the first time. In this series of papers, using the SXS we have measured the velocities of gas motions, metallicities and the multi-temperature structure of the gas in the core of the Perseus cluster. Here, we show that when inferring physical properties from line emissivities in systems like Perseus, the resonant scattering (RS) effect should be taken into account. In the Hitomi waveband, RS mostly affects the FeXXV He$\alpha$ line ($w$) - the strongest line in the spectrum. The flux measured by Hitomi in this line is suppressed by a factor ~1.3 in the inner ~30 kpc, compared to predictions for an optically thin plasma; the suppression decreases with the distance from the center. The $w$ line also appears slightly broader than other lines from the same ion. The observed distortions of the $w$ line flux, shape and distance dependence are all consistent with the expected effect of the resonant scattering in the Perseus core. By measuring the ratio of fluxes in optically thick ($w$) and thin (FeXXV forbidden, He$\beta$, Ly$\alpha$) lines, and comparing these ratios with predictions from Monte Carlo radiative transfer simulations, the velocities of gas motions have been obtained. The results are consistent with the direct measurements of gas velocities from line broadening described elsewhere in this series, although the systematic and statistical uncertainties remain significant. Further improvements in the predictions of line emissivities in plasma models, and deeper observations with future X-ray missions will enable RS measurements to provide powerful constraints on the amplitude and anisotropy of clusters gas motions.

To search for giant X-ray pulses correlated with the giant radio pulses (GRPs) from the Crab pulsar, we performed a simultaneous observation of the Crab pulsar with the X-ray satellite Hitomi in the 2 -- 300 keV band and the Kashima NICT radio observatory in the 1.4 -- 1.7 GHz band with a net exposure of about 2 ks on 25 March 2016, just before the loss of the Hitomi mission.The timing performance of the Hitomi instruments was confirmed to meet the timing requirement and about 1,000 and 100 GRPs were simultaneously observed at the main and inter-pulse phases, respectively, and we found no apparent correlation between the giant radio pulses and the X-ray emission in either the main or inter-pulse phases.All variations are within the 2 sigma fluctuations of the X-ray fluxes at the pulse peaks, and the 3 sigma upper limits of variations of main- or inter- pulse GRPs are 22\% or 80\% of the peak flux in a 0.20 phase width, respectively, in the 2 -- 300 keV band.The values become 25\% or 110\% for main or inter-pulse GRPs, respectively, when the phase width is restricted into the 0.03 phase.Among the upper limits from the Hitomi satellite, those in the 4.5-10 keV and the 70-300 keV are obtained for the first time, and those in other bands are consistent with previous reports.Numerically, the upper limits of main- and inter-pulse GRPs in the 0.20 phase width are about (2.4 and 9.3) $\times 10^{-11}$ erg cm$^{-2}$, respectively. No significant variability in pulse profiles implies that the GRPs originated from a local place within the magnetosphere and the number of photon-emitting particles temporally increases.However, the results do not statistically rule out variations correlated with the GRPs, because the possible X-ray enhancement may appear due to a $>0.02$\% brightening of the pulse-peak flux under such conditions.

The icy satellites around Jupiter are considered to have formed in a circumplanetary disk. While previous models focused on the formation of satellites starting from satellitesimals, the question of how satellitesimals form from smaller dust particles has not been addressed so far. In this work, we study the possibility that satellitesimals form in situ in a circumplanetary disk. We calculate the radial distribution of the surface density and representative size of icy dust particles that grow by colliding with each other and drift toward the central planet in a steady circumplanetary disk with a continuous supply of gas and dust from the parent protoplanetary disk. The radial drift barrier is overcome if the ratio of the dust to gas accretion rates onto the circumplanetary disk, $\dot{M}_{\mathrm{d}}/\dot{M}_{\mathrm{g}}$, is high and the strength of turbulence, $\alpha$, is not too low. The collision velocity is lower than the critical velocity of fragmentation when $\alpha$ is low. Taken together, we find that the conditions for satellitesimal formation via dust coagulation are given by $\dot{M}_{\mathrm{d}}/\dot{M}_{\mathrm{g}}\ge1$ and $10^{-4}\le\alpha<10^{-2}$. The former condition is generally difficult to achieve, suggesting that the in-situ satellitesimal formation via particle sticking is viable only under an extreme condition. We also show that neither satellitesimal formation via the collisional growth of porous aggregates nor via streaming instability is viable as long as $\dot{M}_{\mathrm{d}}/\dot{M}_{\mathrm{g}}$ is low.

The Crab nebula originated from a core-collapse supernova (SN) explosion observed in 1054 A.D. When viewed as a supernova remnant (SNR), it has an anomalously low observed ejecta mass and kinetic energy for an Fe-core collapse SN. Intensive searches were made for a massive shell that solves this discrepancy, but none has been detected. An alternative idea is that the SN1054 is an electron-capture (EC) explosion with a lower explosion energy by an order of magnitude than Fe-core collapse SNe. In the X-rays, imaging searches were performed for the plasma emission from the shell in the Crab outskirts to set a stringent upper limit to the X-ray emitting mass. However, the extreme brightness of the source hampers access to its vicinity. We thus employed spectroscopic technique using the X-ray micro-calorimeter onboard the Hitomi satellite. By exploiting its superb energy resolution, we set an upper limit for emission or absorption features from yet undetected thermal plasma in the 2-12 keV range. We also re-evaluated the existing Chandra and XMM-Newton data. By assembling these results, a new upper limit was obtained for the X-ray plasma mass of <~ 1Mo for a wide range of assumed shell radius, size, and plasma temperature both in and out of the collisional equilibrium. To compare with the observation, we further performed hydrodynamic simulations of the Crab SNR for two SN models (Fe-core versus EC) under two SN environments (uniform ISM versus progenitor wind). We found that the observed mass limit can be compatible with both SN models if the SN environment has a low density of <~ 0.03 cm-3 (Fe core) or <~ 0.1 cm-3 (EC) for the uniform density, or a progenitor wind density somewhat less than that provided by a mass loss rate of 10-5 Mo yr-1 at 20 km s-1 for the wind environment.

The Earth was born in violence. Many giant collisions of protoplanets are thought to have occurred during the terrestrial planet formation. Here we investigated the giant impact stage by using a hybrid code that consistently deals with the orbital evolution of protoplanets around the Sun and the details of processes during giant impacts between two protoplanets. A significant amount of materials (up to several tens of percent of the total mass of the protoplanets) is ejected by giant impacts. We call these ejected fragments the giant-impact fragments (GIFs). In some of the erosive hit-and-run and high-velocity collisions, metallic iron is also ejected, which comes from the colliding protoplanets' cores. From ten numerical simulations for the giant impact stage, we found that the mass fraction of metallic iron in GIFs ranges from ~ 1wt% to ~ 25wt%. We also discussed the effects of the GIFs on the dynamical and geochemical characteristics of formed terrestrial planets. We found that the GIFs have the potential to solve the following dynamical and geochemical conflicts: (1) The Earth, currently in a near circular orbit, is likely to have had a highly eccentric orbit during the giant impact stage. The GIFs are large enough in total mass to lower the eccentricity of the Earth to its current value via their dynamical friction. (2) The concentrations of highly siderophile elements (HSEs) in the Earth's mantle are greater than what was predicted experimentally. Re-accretion of the iron-bearing GIFs onto the Earth can contribute to the excess of HSEs. In addition, the estimated amount of iron-bearing GIFs provides significant reducing agent that could transform primitive CO2-H2O atmosphere and ocean into more reducing H2-bearing atmosphere. Thus, GIFs are important for the origin of Earth's life and its early evolution.

High-resolution X-ray spectroscopy with Hitomi was expected to resolve the origin of the faint unidentified E=3.5 keV emission line reported in several low-resolution studies of various massive systems, such as galaxies and clusters, including the Perseus cluster. We have analyzed the Hitomi first-light observation of the Perseus cluster. The emission line expected for Perseus based on the XMM-Newton signal from the large cluster sample under the dark matter decay scenario is too faint to be detectable in the Hitomi data. However, the previously reported 3.5 keV flux from Perseus was anomalously high compared to the sample-based prediction. We find no unidentified line at the reported high flux level. Taking into account the XMM measurement uncertainties for this region, the inconsistency with Hitomi is at a 99% significance for a broad dark-matter line and at 99.7% for a narrow line from the gas. We do not find anomalously high fluxes of the nearby faint K line or the Ar satellite line that were proposed as explanations for the earlier 3.5 keV detections. We do find a hint of a broad excess near the energies of high-n transitions of Sxvi (E=3.44 keV rest-frame) -- a possible signature of charge exchange in the molecular nebula and another proposed explanation for the unidentified line. While its energy is consistent with XMM pn detections, it is unlikely to explain the MOS signal. A confirmation of this interesting feature has to wait for a more sensitive observation with a future calorimeter experiment.

We report $Suzaku$ X-ray observations of the dark subhalo associated with the merging group of NGC 4839 in the Coma cluster. The X-ray image exhibits an elongated tail toward the southwest. The X-ray peak shifts approximately $1'$ away from the weak-lensing mass center toward the opposite direction of the Coma cluster center. We investigated the temperature, normalization, pressure, and entropy distributions around the subhalo. Excluding the X-ray tail, the temperature beyond the truncation radius is 8-10$~\rm keV$, which is two times higher than that of the subhalo and the X-ray tail. The pressure is nearly uniform excluding southern part of the subhalo at two times of the truncation radius. We computed the gas mass within the truncation radius and the X-ray tail. While the gas fraction within the truncation radius is about 5 times smaller than that of regular groups, the gas mass in the subhalo and the X-ray tail to weak-lensing mass ratio is consistent with that of regular groups. Assuming the infall velocity, $2000~\rm km~s^{-1}$, the ram pressure is 1.4 times greater than gravitational force per unit area. Assuming the Kelvin-Helmholtz instabilities, the total lost mass is approximately $3\times10^{11}~M_{\odot}$. If this gas had originally been within the truncation radius, the gas mass fraction of the subhalo would have been comparable with those of regular groups before infalling to the Coma cluster.

Population synthesis studies into planet formation have suggested that distributions consistent with observations can only be reproduced if the actual Type I migration timescale is at least an order of magnitude longer than that deduced from linear theories. Although past studies considered the effect of the Type I migration of protoplanetary embryos, in most cases they used a conventional formula based on static torques in isothermal disks, and employed a reduction factor to account for uncertainty in the mechanism details. However, in addition to static torques, a migrating planet experiences dynamic torques that are proportional to the migration rate. These dynamic torques can impact on planet migration and predicted planetary populations. In this study, we derived a new torque formula for Type I migration by taking into account dynamic corrections. This formula was used to perform population synthesis simulations with and without the effect of dynamic torques. In many cases, inward migration was slowed significantly by the dynamic effects. For the static torque case, gas giant formation was effectively suppressed by Type I migration; however, when dynamic effects were considered, a substantial fraction of cores survived and grew into gas giants.

We perform a weak-lensing study of the nearby cool-core galaxy clusters, Hydra A ($z=0.0538$) and A478 ($z=0.0881$), of which brightest cluster galaxies (BCGs) host powerful activities of active galactic nuclei (AGNs). For each cluster, the observed tangential shear profile is well described either by a single Navarro--Frenk--White model or a two-component model including the BCG as an unresolved point mass. For A478, we determine the BCG and its host-halo masses from a joint fit to weak-lensing and stellar photometry measurements. We find that the choice of initial mass functions (IMFs) can introduce a factor of two uncertainty in the BCG mass, whereas the BCG host halo mass is well constrained by data. We perform a joint analysis of weak-lensing and stellar kinematics data available for the Hydra A cluster, which allows us to constrain the central mass profile without assuming specific IMFs. We find that the central mass profile ($r<300$kpc) determined from the joint analysis is in excellent agreement with those from independent measurements, including dynamical masses estimated from the cold gas disk component, X-ray hydrostatic total mass estimates, and the central stellar mass estimated based on the Salpeter IMF. The observed dark-matter fraction around the BCG for Hydra A is found to be smaller than those predicted by adiabatic contraction models, suggesting the importance of other physical processes, such as the the AGN feedback and/or dissipationless mergers.

We obtain a constraint on the parameters of a generalized cubic Galileon gravity model exhibiting the Vainshtein mechanism by using multi-wavelength observations of the Coma Cluster. The generalized cubic Galileon model is characterized by three parameters of the turning scale associated with the Vainshtein mechanism, and the amplitude of modifying a gravitational potential and a lensing potential. X-ray and Sunyaev-Zeldovich (SZ) observations of the intra-cluster medium are sensitive to the gravitational potential, while the weak-lensing (WL) measurement is specified by the lensing potential. A joint fit of a complementary multi-wavelength dataset of X-ray, SZ and WL measurements enables us to simultaneously constrain these three parameters of the generalized cubic Galileon model for the first time. We also find a degeneracy between the cluster mass parameters and the gravitational modification parameters, which is influential in the limit of the weak screening of the fifth force.

No Hadean rocks have ever been found on Earth's surface except for zircons---evidence of continental crust, suggesting that Hadean continental crust existed but later disappeared. One hypothesis for the disappearance of the continental crust is excavation/melting by the Late Heavy Bombardment (LHB), a concentration of impacts in the last phase of the Hadean eon. In this paper, we calculate the effects of LHB on Hadean continental crust in order to investigate this hypothesis. Approximating the size-frequency distribution of the impacts by a power-law scaling with an exponent {\alpha} as a parameter, we have derived semi-analytical expressions for the effects of LHB impacts. We calculated the total excavation/melting volume and area affected by the LHB from two constraints of LHB on the moon, the size of the largest basin during LHB, and the density of craters larger than 20 km. We also investigated the effects of the value of {\alpha}. Our results show that LHB does not excavate/melt all of Hadean continental crust directly, but over 70% of the Earth's surface area can be covered by subsequent melts in a broad range of {\alpha}. If there have been no overturns of the continental crust until today, LHB could be responsible for the absence of Hadean rocks because most of Hadean continental crust is not be exposed on the Earth's surface in this case.

We observed three massive subhalos in the Coma cluster with {\it Suzaku}. These subhalos, labeled "ID 1", "ID 2", and "ID 32", were detected with a weak-lensing survey using the Subaru/Suprime-Cam (Okabe et al. 2014a), and are located at the projected distances of 1.4 $r_{500}$, 1.2 $r_{500}$, and 1.6 $r_{500}$ from the center of the Coma cluster, respectively. The subhalo "ID 1" has a compact X-ray excess emission close to the center of the weak-lensing mass contour, and the gas mass to weak-lensing mass ratio is about 0.001. The temperature of the emission is about 3 keV, which is slightly lower than that of the surrounding intracluster medium (ICM) and that expected for the temperature vs. mass relation of clusters of galaxies. The subhalo "ID 32" shows an excess emission whose peak is shifted toward the opposite direction from the center of the Coma cluster. The gas mass to weak-lensing mass ratio is also about 0.001, which is significantly smaller than regular galaxy groups. The temperature of the excess is about 0.5 keV and significantly lower than that of the surrounding ICM and far from the temperature vs. mass relation of clusters. However, there is no significant excess X-ray emission in the "ID 2" subhalo. Assuming an infall velocity of about 2000 $\rm km~s^{-1}$, at the border of the excess X-ray emission, the ram pressures for "ID 1" and "ID 32" are comparable to the gravitational restoring force per area. We also studied the effect of the Kelvin-Helmholtz instability to strip the gas. Although we found X-ray clumps associated with the weak-lensing subhalos, their X-ray luminosities are much lower than the total ICM luminosity in the cluster outskirts.

The diversity and quantity of moons in the Solar System suggest a manifold population of natural satellites exist around extrasolar planets. Of peculiar interest from an astrobiological perspective, the number of sizable moons in the stellar habitable zones may outnumber planets in these circumstellar regions. With technological and theoretical methods now allowing for the detection of sub-Earth-sized extrasolar planets, the first detection of an extrasolar moon appears feasible. In this review, we summarize formation channels of massive exomoons that are potentially detectable with current or near-future instruments. We discuss the orbital effects that govern exomoon evolution, we present a framework to characterize an exomoon's stellar plus planetary illumination as well as its tidal heating, and we address the techniques that have been proposed to search for exomoons. Most notably, we show that natural satellites in the range of 0.1 - 0.5 Earth mass (i) are potentially habitable, (ii) can form within the circumplanetary debris and gas disk or via capture from a binary, and (iii) are detectable with current technology.

We conduct a joint X-ray and weak-lensing study of four relaxed galaxy clusters (Hydra A, A478, A1689 and A1835) observed by both Suzaku and Subaru out to virial radii, with an aim to understand recently-discovered unexpected feature of the ICM in cluster outskirts. We show that the average hydrostatic-to-lensing total mass ratio for the four clusters decreases from \sim 70% to \sim 40% as the overdensity contrast decreases from 500 to the virial value.The average gas mass fraction from lensing total mass estimates increases with cluster radius and agrees with the cosmic mean baryon fraction within the virial radius, whereas the X-ray-based gas fraction considerably exceeds the cosmic values due to underestimation of the hydrostatic mass. We also develop a new advanced method for determining normalized cluster radial profiles for multiple X-ray observables by simultaneously taking into account both their radial dependence and multivariate scaling relations with weak-lensing masses. Although the four clusters span a range of halo mass, concentration, X-ray luminosity and redshift, we find that the gas entropy, pressure, temperature and density profiles are all remarkably self-similar when scaled with the lensing M_200 mass and r_200 radius.The entropy monotonically increases out to \sim 0.5r_200 following the accretion shock heating model K(r)\propto r^1.1, and flattens at \simgt 0.5r_200.The universality of the scaled entropy profiles indicates that the thermalization mechanism over the entire cluster region (>0.1r_200) is controlled by gravitation in a common to all clusters, although the heating efficiency in the outskirts needs to be modified from the standard law.The bivariate scaling functions of the gas density and temperature reveal that the flattening of the outskirts entropy profile is caused by the steepening of the temperature, rather than the flattening of the gas density.

We studied the distributions of metal abundances and metal-mass-to-light ratios in the intracluster medium (ICM) of four galaxy groups, MKW 4, HCG 62, the NGC 1550 group, and the NGC 5044 group, out to 0.5 r180 observed with Suzaku. The Fe abundance decreases with radius, and about 0.2-0.4 solar beyond 0.1 r180. At a given radius in units of r180, the Fe abundance in the ICM of the four galaxy groups were consistent or smaller than those of clusters of galaxies. The Mg/Fe and Si/Fe ratios in the ICM are nearly constant at the solar ratio out to 0.5 r180. We also studied systematic uncertainties in the derived metal abundances comparing the results from two versions of atomic data for astrophysicists (ATOMDB) and single- and two temperature model fits. Since the metals have been synthesized in galaxies, we collected K-band luminosities of galaxies from Two Micron All Sky Survey catalogue (2MASS) and calculated the integrated iron-mass-to-light-ratios (IMLR), or the ratios of the iron mass in the ICM to light from stars in galaxies. The groups with smaller gas mass to light ratios have smaller IMLR values and the IMLR inversely correlated with the entropy excess. Based on these abundance features, we discussed the past history of metal enrichment process in groups of galaxies.

Many extrasolar (bound) terrestrial planets and free-floating (unbound) planets have been discovered. The existence of bound and unbound terrestrial planets with liquid water is an important question, and of particular importance is the question of their habitability. Even for a globally ice-covered planet, geothermal heat from the planetary interior may melt the interior ice, creating an internal ocean covered by an ice shell. In this paper, we discuss the conditions that terrestrial planets must satisfy for such an internal ocean to exist on the timescale of planetary evolution. The question is addressed in terms of planetary mass, distance from a central star, water abundance, and abundance of radiogenic heat sources. In addition, we investigate the structures of the surface-H2O layers of ice-covered planets by considering the effects of ice under high pressure (high-pressure ice). As a fiducial case, 1M$\oplus$ planet at 1 AU from its central star and with 0.6 to 25 times the H2O mass of Earth could have an internal ocean. We find that high-pressure ice layers may appear between the internal ocean and the rock portion on a planet with an H2O mass over 25 times that of Earth. The planetary mass and abundance of surface water strongly restrict the conditions under which an extrasolar terrestrial planet may have an internal ocean with no high-pressure ice under the ocean. Such high-pressure-ice layers underlying the internal ocean are likely to affect the habitability of the planet.

We determine the quark-hadron transition line in the whole region of temperature (T) and baryon-number chemical potential (\mu_B) from lattice QCD results and neutron-star mass measurements, making the quark-hadron hybrid model that is consistent with the two solid constraints. The quark part of the hybrid model is the Polyakov-loop extended Nambu-Jona-Lasinio (PNJL) model with entanglement vertex that reproduces lattice QCD results at mu_B/T=0, while the hadron part is the hadron resonance gas model with volume-exclusion effect that reproduces neutron-star mass measurements and the neutron-matter equation of state calculated from two- and three-nucleon forces based on the chiral effective field theory. The lower bound of the critical mu_B of the quark-hadron transition at zero T is mu_B = 1.6 GeV. The interplay between the heavy-ion collision physics around mu_B/T =6 and the neutron-star physics where mu_B/T is infinity is discussed.

We present results from a large mosaic of Suzaku observations of the Coma Cluster, the nearest and X-ray brightest hot, dynamically active, non-cool core system, focusing on the thermodynamic properties of the ICM on large scales. For azimuths not aligned with an infalling subcluster towards the southwest, our measured temperature and X-ray brightness profiles exhibit broadly consistent radial trends, with the temperature decreasing from about 8.5 keV at the cluster center to about 2 keV at a radius of 2 Mpc, which is the edge of our detection limit. The SW merger significantly boosts the surface brightness, allowing us to detect X-ray emission out to ~2.2 Mpc along this direction. Apart from the southwestern infalling subcluster, the surface brightness profiles show multiple edges around radii of 30-40 arcmin. The azimuthally averaged temperature profile, as well as the deprojected density and pressure profiles, all show a sharp drop consistent with an outward propagating shock front located at 40 arcmin, corresponding to the outermost edge of the giant radio halo observed at 352 MHz with the WSRT. The shock front may be powering this radio emission. A clear entropy excess inside of r_500 reflects the violent merging events linked with these morphological features. Beyond r_500, the entropy profiles of the Coma Cluster along the relatively relaxed directions are consistent with the power-law behavior expected from simple models of gravitational large-scale structure formation. The pressure is also in agreement at these radii with the expected values measured from SZ data from the Planck satellite. However, due to the large uncertainties associated with the Coma Cluster measurements, we cannot yet exclude an entropy flattening in this system consistent with that seen in more relaxed cool core clusters.

We analyzed XMM-Newton data of the Perseus cluster out to $\sim$1 Mpc, or approximately half the virial radius. Using the flux ratios of Lyalpha lines of H-like Si and S to Kalpha line of He-like Fe, the abundance ratios of Si/Fe and S/Fe of the intracluster medium (ICM) were derived using the APEC plasma code v2.0.1. The temperature dependence of the line ratio limits the systematic uncertainty in the derived abundance ratio. The Si/Fe and S/Fe in the ICM of the Perseus cluster show no radial gradient. The emission-weighted averages of the Si/Fe and S/Fe ratios outside the cool core are 0.91 +- 0.08 and 0.93 +- 0.10, respectively, in solar units according to the solar abundance table of Lodders (2003). These ratios indicate that most Fe was synthesized by supernovae Ia. We collected K-band luminosities of galaxies and calculated the ratio of Fe and Si mass in the ICM to K-band luminosity, iron-mass-to-light ratio (IMLR) and silicon-mass-to-light ratio (SMLR). Within $\sim$1 Mpc, the cumulative IMLR and SMLR increase with radius. Using Suzaku data for the northwest and east directions, we also calculated the IMLR out to $\sim$ 1.8 Mpc, or about the virial radius. We constrained the SMLR out to this radius and discussed the slope of the initial mass function of stars in the cluster. Using the cumulative IMLR profile, we discuss the past supernova Ia rate.

We report Suzaku observations of the northern half of the Hydra A cluster out to ~1.4 Mpc, reaching the virial radius. This is the first Suzaku observations of a medium-size (kT ~3 keV) cluster out to the virial radius. Two observations were conducted, north-west and north-east offsets, which continue in a filament direction and a void direction of the large-scale structure of the Universe, respectively. The X-ray emission and distribution of galaxies elongate in the filament direction. The temperature profiles in the two directions are mostly consistent with each other within the error bars and drop to 1.5 keV at 1.5 r_500. As observed by Suzaku in hot clusters, the entropy profile becomes flatter beyond r_500, in disagreement with the r^1.1 relationship that is expected from accretion shock heating models. When scaled with the average intracluster medium (ICM) temperature, the entropy profiles of clusters observed with Suzaku are universal and do not depend on system mass. The hydrostatic mass values in the void and filament directions are in good agreement, and the Navarro, Frenk, and White universal mass profile represents the hydrostatic mass distribution up to ~ 2 r_500. Beyond r_500, the ratio of gas mass to hydrostatic mass exceeds the result of the Wilkinson microwave anisotropy probe, and at r_100, these ratios in the filament and void directions reach 0.4 and 0.3, respectively. We discuss possible deviations from hydrostatic equilibrium at cluster outskirts. We derived radial profiles of the gasmass- to-light ratio and iron-mass-to-light ratio out to the virial radius. Within r_500, the iron-mass-to-light ratio of the Hydra A cluster was compared with those in other clusters observed with Suzaku.

We formulated tidal decay lifetimes for hypothetical moons orbiting extrasolar planets with both lunar and stellar tides. Previous work neglected the effect of lunar tides on planet rotation, and are therefore applicable only to systems in which the moon's mass is much less than that of the planet. This work, in contrast, can be applied to the relatively large moons that might be detected around newly-discovered Neptune-mass and super-Earth planets. We conclude that moons are more stable when the planet/moon systems are further from the parent star, the planets are heavier, or the parent stars are lighter. Inclusion of lunar tides allows for significantly longer lifetimes for a massive moon relative to prior formulations. We expect that the semi-major axis of the planet hosting the first detected exomoon around a G-type star is 0.4-0.6 AU and is 0.2-0.4 AU for an M-type star.

We study the multi-dimensional geometry of supernova (SN) explosions by means of spectropolarimetric observations of stripped-envelope SNe, i.e., SNe without a H-rich layer. We perform spectropolarimetric observations of 2 stripped-envelope SNe, the Type Ib SN 2009jf and the Type Ic SN 2009mi. Both objects show non-zero polarization at the wavelength of the strong lines. They also show a loop in the Stokes Q-U diagram, which indicates a non-axisymmetric, three-dimensional ion distribution in the ejecta. We show that five out of six stripped-envelope SNe which have been observed spectropolarimetrically so far show such a loop. This implies that a three-dimensional geometry is common in stripped-envelope SNe. We find that stronger lines tend to show higher polarization. This effect is not related to the geometry, and must be corrected to compare the polarization of different lines or different objects. Even after the correction, however, there remains a dispersion of polarization degree among different objects. Such a dispersion might be caused by three-dimensional clumpy ion distributions viewed from different directions.

The Jovian regular satellite system mainly consists of four Galilean satellites that have similar masses and are trapped in mutual mean motion resonances except for the outer satellite, Callisto. On the other hand, the Saturnian regular satellite system has only one big icy body, Titan, and a population of much smaller icy moons. We have investigated the origin of these major differences between the Jovian and Saturnian satellite systems by semi-analytically simulating the growth and orbital migration of proto-satellites in an accreting proto-satellite disk. We set up two different disk evolution/structure models that correspond to Jovian and Saturnian systems, by building upon previously developed models of an actively-supplied proto-satellite disk, the formation of gas giants, and observations of young stars. Our simulations extend previous models by including the (1) different termination timescales of gas infall onto the proto-satellite disk and (2) different evolution of a cavity in the disk, between the Jovian and Saturnian systems. We have performed Monte Carlo simulations and show that in the case of the Jovian systems, four to five similar-mass satellites are likely to remain trapped in mean motion resonances. This orbital configuration is formed by type I migration, temporal stopping of the migration near the disk inner edge, and quick truncation of gas infall caused by Jupiter opening a gap in the Solar nebula. The Saturnian systems tend to end up with one dominant body in the outer regions caused by the slower decay of gas infall associated with global depletion of the Solar nebula. The total mass and compositional zoning of the predicted Jovian and Saturnian satellite systems are consistent with the observed satellite systems.

We present the first spectropolarimetric observations of a candidate of super-Chandrasekhar mass Type Ia supernova (SN): SN 2009dc. The observations were performed at 5.6 and 89.5 days after the B-band maximum. The data taken at the later epoch are used to determine the interstellar polarization. Continuum polarization is found to be small (<0.3 %), indicating that the explosion is nearly spherically symmetric. This fact suggests that a very aspherical explosion is not a likely scenario for SN 2009dc. Polarization at the Si II and Ca II lines clearly shows a loop in the Q-U plane, indicating a non-axisymmetric, clumpy distribution of intermediate-mass elements. The degree of line polarization at the Si and Ca line is moderate (0.5% +- 0.1% and 0.7% +- 0.1%, respectively), but it is higher than expected from the trend of other Type Ia SNe. This may suggest that there are thick enough, clumpy Si-rich layers above the thick 56Ni-rich layers (>~ 1.2 Msun). The observed spectropolarimetric properties, combined with the photometric and spectroscopic properties, suggest that the progenitor of SN 2009dc has a super-Chandrasekhar mass, and that the explosion geometry is globally spherically symmetric, with clumpy distribution of intermediate-mass elements.

Core-collapse supernovae (CC-SNe) are the explosions that announce the death of massive stars. Some CC-SNe are linked to long-duration gamma-ray bursts (GRBs) and are highly aspherical. One important question is to what extent asphericity is common to all CC-SNe. Here we present late-time spectra for a number of CC-SNe from stripped-envelope stars, and use them to explore any asphericity generated in the inner part of the exploding star, near the site of collapse. A range of oxygen emission-line profiles is observed, including a high incidence of double-peaked profiles, a distinct signature of an aspherical explosion. Our results suggest that all CC-SNe from stripped-envelope stars are aspherical explosions and that SNe accompanied by GRBs exhibit the highest degree of asphericity.

We present multi-waveband optical imaging data obtained from observations of the Subaru/XMM-Newton Deep Survey (SXDS). The survey field, centered at R.A.=02:18:00, decl.=-05:00:00, has been the focus of a wide range of multi-wavelength observing programs spanning from X-ray to radio wavelengths. A large part of the optical imaging observations are carried out with Suprime-Cam on Subaru Telescope at Mauna Kea in the course of Subaru Telescope Observatory Projects. This paper describes our optical observations, data reduction and analysis procedures employed, and the characteristics of the data products. A total area of 1.22 sqdeg is covered in five contiguous sub-fields, each of which corresponds to a single Suprime-Cam field of view (34'x27'), in five broad-band filters B, V, Rc, i', z' to the depths of B=28.4, V=27.8, Rc=27.7, i'=27.7 and z'=26.6 (AB, 3-sigma, 2-arcsec aperture). The data are reduced and compiled into five multi-waveband photometric catalogs, separately for each Suprime-Cam pointing. The i'-band catalogs contain about 900,000 objects, making the SXDS catalogs one of the largest multi-waveband catalogs in corresponding depth and area coverage. The SXDS catalogs can be used for an extensive range of astronomical applications such as the number density of the Galactic halo stars to the large scale structures at the distant universe. The number counts of galaxies are derived and compared with those of existing deep extragalactic surveys. The optical data, the source catalogs, and configuration files used to create the catalogs are publicly available via the SXDS web page (http://www.naoj.org/Science/SubaruProject/SXDS/index.html)

We report an extensive search for Lyman-alpha emitters (LAEs) at z=6.5 in the Subaru Deep Field. Subsequent spectroscopy with Subaru and Keck identified eight more LAEs, giving a total of 17 spectroscopically confirmed LAEs at z=6.5. Based on this spectroscopic sample of 17, complemented by a photometric sample of 58 LAEs, we have derived a more accurate Lyman-alpha luminosity function of LAEs at z=6.5, which reveals an apparent deficit at the bright end of ~0.75 mag fainter L*, compared with that observed at z=5.7. The difference in the LAE luminosity functions between z=5.7 and 6.5 is significant at the 3-sigma level, which is reduced to 2-sigma when cosmic variance is taken into account. This result may imply that the reionization of the universe has not been completed at z=6.5. We found that the spatial distribution of LAEs at z=6.5 was homogeneous over the field. We discuss the implications of these results for the reionization of the universe.

We present the properties of Ly Alpha emitters (LAEs) at z = 5.7 in the Subaru Deep Field. A photometric sample of 89 LAE candidates is constructed from narrow-band (NB816) data down to NB816 = 26.0 (AB) in a continuous 725 arcmin^2 area. Spectra of 39 objects satisfying the photometric selection criteria for LAEs were obtained with Subaru and Keck II Telescopes, among which 28 were confirmed LAEs, one was a nearby galaxy, and eight were unclassified. We also obtained spectra of another 24 NB816-excess objects in the field, identifying six additional LAEs. We find that the Ly Alpha luminosity function derived from the photometric sample is reproduced well by a Schechter function with L* = (7.9+3.0-2.2) x 10^42 erg/s and phi* = (6.3+3.0-2.0) x 10^-4 Mpc^-3 for alpha = -1.5 (fixed) over the whole luminosity range of L ~= 3x10^42 - 3x10^43 erg/s. We then measure rest-frame Ly Alpha equivalent widths (EWs) for the confirmed LAEs, to find that the median among the 28 objects satisfying the photometric selection criteria is W_0^i = 233 A. We infer that 30% - 40% of LAEs at z=5.7 exceed W_0^i = 240 A. These large-EW objects probably cannot be accounted for by ordinary star-forming populations with a Salpeter IMF. We also find that LAEs with fainter far-UV luminosities have larger EWs. Finally, we derive the far-UV luminosity function of LAEs down to M_UV ~= -19.6 using the photometric sample, and compare it with that of Lyman-break galaxies (LBGs). We find that as high as about 80% of LBGs at z ~ 6 have W_0^i >= 100 A, in sharp contrast to lower-z counterparts.

We report angular correlation function (ACF) of Lyman Break Galaxies (LBGs) with unprecedented statistical quality on the basis of 16,920 LBGs at z=4 detected in the 1 deg^2 sky of the Subaru/XMM-Newton Deep Field. The ACF significantly departs from a power law, and shows an excess on small scale. Particularly, the ACF of LBGs with i'<27.5 have a clear break between the small and large-scale regimes at the angular separation of ~7'' whose projected length corresponds to the virial radius of dark halos with a mass of 10^11-12 Mo, indicating multiple LBGs residing in a single dark halo. Both on small (2''<theta<3'') and large (40''<theta<400'') scales, clustering amplitudes monotonically increase with luminosity for the magnitude range of i'=24.5-27.5, and the small-scale clustering shows a stronger luminosity dependence than the large-scale clustering. The small-scale bias reaches b~10-50, and the outskirts of small-scale excess extend to a larger angular separation for brighter LBGs. The ACF and number density of LBGs can be explained by the cold dark matter model.

We report the discovery of primeval large-scale structures (LSSs) including two proto-clusters in a forming phase at z=5.7. We carried out extensive deep narrow-band imaging in the 1 deg^2 sky of the Subaru/XMM-Newton Deep Field, and obtained a cosmic map of 515 Ly alpha emitters (LAEs) in a volume with a transverse dimension of 180 Mpc x 180 Mpc and a depth of ~40 Mpc in comoving units. This cosmic map shows filamentary LSSs, including clusters and surrounding 10-40 Mpc scale voids, similar to the present-day LSSs. Our spectroscopic follow-up observations identify overdense regions in which two dense clumps of LAEs with a sphere of 1-Mpc diameter in physical units are included. These clumps show about 130 times higher star-formation rate density, mainly due to a large overdensity, ~80, of LAEs. These clumps would be clusters in a formation phase involving a burst of galaxy formation.

The Subaru Deep Field (SDF) project is a program of Subaru Observatory to carry out a deep galaxy survey over a blank field as large as 34'x27'. The program consists of very deep multi-band optical imaging, near infrared imaging for smaller portions of the field and follow-up optical spectroscopy. Major scientific goals of the project are to construct large samples of Lyman-break galaxies at z~4-5 and Lyman alpha emitters at z~5.7 and 6.6, and to make detailed studies these very high-redshift galaxy populations. In this paper, we describe the optical imaging observations and data reduction, presenting mosaicked images and object catalogs in seven bandpasses.The optical imaging was made through five broad-band filters, B, V, R, i', z', and two narrow-band filters, NB816 (lambda_c=8150A) and NB921 (lambda_c=9196A) with almost 10 hours long integrations for each band. The limiting magnitudes measured at 3-sigma on a 2" aperture are B=28.45, V=27.74, R=27.80, i'=27.43, z'=26.62, NB816=26.63, and NB921=26.54 in the AB system. The object catalog constructed for each of the seven bands contains more than 10^5 objects. The galaxy number counts corrected for detection incompleteness and star count contribution are found to be consistent with previous results in the literature. The mosaicked images and catalogs of all the bands have been made open to the public on Oct. 1, 2004 on the SDF project website at http://soaps.naoj.org/sdf/.

We report on the discovery of a star-forming galaxy at z=6.33 in the Subaru Deep Field. This object is selected as a candidate of an i'-dropout, high-redshift galaxy around z=6 because of its red i'-z' color in our deep optical imaging survey in the Subaru Deep Field. Our follow up optical spectroscopy reveals that this object is a strong Ly-alpha emitter with only very faint ultraviolet continuum. The rest-frame equivalent width of the detected Ly-alpha emission is as much as 130 A. Thus the light detected in our z' image is largely attributed to the Ly-alpha emission, i.e., ~40% of the z'-band flux is the strong Ly-alpha emission, giving a very red i'-z' color. This is consistent with the photometric property of this object because the narrow-band data obtained with the NB921 filter shows a significant depression, z'-NB921 = -0.54 mag. By using the photometric data, we show that some other objects among the 48 i'-dropout high-redshift galaxy candidates found in the Subaru Deep Field also show a significant NB921 depression. We briefly discuss the nature of these NB921-depressed objects.

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.

We report here the results of deep optical spectroscopy of the very extended emission-line region (VEELR) found serendipitously around the Seyfert 2 galaxy NGC 4388 in the Virgo cluster using the Subaru Telescope. The H-alpha recession velocities of most of the filaments of the region observed are highly blue-shifted with respect to the systemic velocity of the galaxy. The velocity field is complicated, and from the kinematic and morphological points of view, there seem to be several streams of filaments: low velocity (v = -100 km/s) filaments, high velocity (v = -300 km/s) filaments, and a very high velocity (v > -500 km/s) cloud complex. The emission-line ratios of the VEELR filaments are well explained by power-law photoionization models with solar abundances, suggesting that the Seyfert nucleus of NGC 4388 is the dominant ionization source of the VEELR and that the VEELR gas has moderate metallicity. In addition to photoionization, shock heating probably contributes to the ionization of the gas. In particular, the filaments outside the ionization cone of the Seyfert nucleus are mainly excited by shocks. We conclude that the VEELR was formerly the disk gas of NGC 4388, which has been stripped by ram pressure due to the interaction between the hot intra-cluster medium (ICM) and the galaxy. The velocity field and the morphology of the VEELR closely resemble snapshots from some numerical simulations of this process. In the case of NGC 4388, the ram pressure-stripped gas, which is normally seen as extended HI filaments, happens to be exposed and ionized by the radiation from the AGN, and so can be seen as optical emission-line gas.

We performed deep K'-band imaging observations of 2 massive clusters, MS 0451.6-0305 at z = 0.55 and MS 0440.5+0204 at z = 0.19, for searching counterparts of the faint sub-mm sources behind these clusters, which would provide one of the deepest extremely red object(ERO) samples. Comparing our near-infrared images with optical images taken by the Hubble Space Telescope and by the Subaru Telescope, we identified 13 EROs in these fields. The sky distributions of EROs are consistent with the previous results, that there is a sign of strong clustering among detected EROs. Also, the surface density with corrected lensing amplification factors in both clusters are in good agreement with that derived from previous surveys. We found 7 EROs and 3 additional very red objects in a small area (\sim 0.6 arcmin^2) of the MS 0451.6-0305 field around an extended SCUBA source. Many of their optical and near-infrared colors are consistent with dusty star-forming galaxies at high redshifts(z \sim 1.0-4.0), and they may be constituting a cluster of dusty starburst galaxies and/or lensed star-forming galaxies at high redshift. Their red J-K' colors and faint optical magnitudes suggest they are relatively old massive stellar systems with ages(>300 Mega years) suffering from dust obscuration. We also found a surface-density enhancement of EROs around the SCUBA source in the MS 0440.5+0204 field.

Spectroscopic and spectropolarimetric observations of SN 2003dh/GRB 030329 obtained in 2003 May using the Subaru 8.2 m telescope are presented. The properties of the SN are investigated through a comparison with spectra of the Type Ic hypernovae SNe 1997ef and 1998bw. (Hypernovae being a tentatively defined class of SNe with very broad absorption features: these features suggest a large velocity of the ejected material and possibly a large explosion kinetic energy.) Comparison with spectra of other hypernovae shows that the spectrum of SN 2003dh obtained on 2003 May 8 and 9, i.e., 34-35 rest-frame days after the GRB (for z=0.1685), are similar to those of SN 1997ef obtained ~34-42 days after the fiducial time of explosion of that SN. The match with SN 1998bw spectra is not as good (at rest 7300-8000 A, but again spectra obtained ~33-43 days after GRB 980425 are preferred. This indicates that the SN may have intermediate properties between SNe 1997ef and 1998bw. Based on the analogy with the other hypernovae, the time of explosion of SN 2003dh is then constrained to be between -8 and +2 days of the GRB. The Si and O P-Cygni lines of SN 2003dh seem comparable to those of SN 1997ef, which suggests that the ejected mass in SN 2003dh may match that in SN 1997ef. Polarization was marginally detected at optical wavelengths. This is consistent with measurements of the late afterglow, implying that it mostly originated in the interstellar medium of the host galaxy.

We have performed a detailed analysis of 3 optically normal galaxies extracted from the XMM Bright Serendipitous Source Sample. Thanks to the good statistics of the XMM-Newton data, we have unveiled the presence of an AGN in all of them. In particular, we detect both X-ray obscured (N_H>10^{22} cm^{-2}) and unobscured (N_H<10^{22} cm^{-2}) AGN with intrinsic 2--10 keV luminosities in the range between 10^{42} -- 10^{43} erg s^{-1}. We find that the X-ray and optical properties of the sources discussed here could be explained assuming a standard AGN hosted by galaxies with magnitudes M_R<M^*, taking properly into account the absorption associated with the AGN, the optical faintness of the nuclear emission with respect to the host galaxy, and the inadequate set--up and atmospheric conditions during the optical spectroscopic observations. Our new spectroscopic observations have revealed the expected AGN features also in the optical band. These results clearly show that optical spectroscopy sometimes can be inefficient in revealing the presence of an AGN, which instead is clearly found from an X-ray spectroscopic investigation. This remarks the importance of being careful in proposing the identification of X-ray sources (especially at faint fluxes) when only low quality optical spectra are in hand. This is particularly important for faint surveys (such as those with XMM-Newton and Chandra), in which optically dull but X-ray active objects are being found in sizeable numbers.

We present spectropolarimetry of the Type Ic supernova SN 2002ap and give a preliminary analysis: the data were taken at two epochs, close to and one month later than the visual maximum (2002 February 8). In addition we present June 9 spectropolarimetry without analysis. The data show the development of linear polarization. Distinct polarization profiles were seen only in the O I \lambda 7773 multiplet/Ca II IR triplet absorption trough at maximum light and in the Ca II IR triplet absorption trough a month later, with the latter showing a peak polarization as high as ~2 %. The intrinsic polarization shows three clear position angles: 80 degs for the February continuum, 120 degs for the February line feature, and 150 degs for the March data. We conclude that there are multiple asymmetric components in the ejecta. We suggest that the supernova has a bulk asymmetry with an axial ratio projected on the sky that is different from 1 by of order 10 %. Furthermore, we suggest very speculatively that a high velocity ejecta component moving faster than ~0.115c (e.g., a jet) contributes to polarization in the February epoch.

We have constructed a very deep K'-selected multicolor BVRIz'JK' sample of 439 field galaxies. Photometric redshifts for sample galaxies were estimated. The overall redshift distribution N(z) for the K'<21.0 sample is consistent with previous observations, and for the first time we derive the N(z) down to K'=24.0. After taking account of the dust extinction and selection effects of the sample, the observed N(z) distribution is well described with the PLE model, while the hierarchical galaxy formation model shows an apparent deficiency of galaxies especially at z>2. The photometric redshift and the best-fit SED model evaluations allow the derivation of the rest-frame K', B, and UV(2000A) luminosity functions (LFs) and their evolutions. The rest-K' LF shows almost no evolution up to z=3, while the rest B LF shows mild luminosity evolution, and the rest UV LF shows strong luminosity evolution. This trend corresponds to the evolution of the rest UV-K' color, which gets bluer with increasing redshift. We also found that more massive galaxies are redder in this rest-frame color in every epoch. The rest-UV LF of our K'-selected galaxies shows a much shallower faint end slope at z=3 than that of previous estimations for rest-UV selected Lyman break galaxies. As a consequence, the contribution to the global star formation rate of our K'-selected galaxies is about 42% of that derived from the integration of LF of Lyman break galaxies at z=3. This result suggests that a large fraction of the star formation rate density at z>1.5 comes from the contribution from the faint (M_2000A>-20) blue galaxy population at high redshift universe that have not yet obviously been identified. (abridged)

We present new optical images ($B$, $V$, and H$\alpha$) of the archetypical starburst/superwind galaxy M82 obtained with the 8.2 m Subaru Telescope to reveal new detailed structures of the superwind-driven nebula and the high-latitude dark lanes. The emission-line nebula is decomposed into (1) a ridge-dominated component comprising numerous filament/loop sub-structures whose overall morphology appears as a pair of narrow cylinders, and (2) a diffuse component extended over much wider opening angle from the nucleus. We suggest that these two components have different origins. The ridge-dominated component appears as a pair of cylinders rather than a pair of cones. Since this morphological property is similar to that of hot plasma probed by soft X-ray, this component seems to surround the hot plasma. On the other hand, the diffuse component may arise from dust grains which scatter stellar light from the galaxy. Since inner region of this component is seen over the prominent ^^ ^^ X"-shaped dark lanes streaming out from the nuclear region and they can be reproduced as a conical distribution of dust grains, there seems to be a dusty cold outflow as well as the hot one probed by soft X-ray and shock-excited optical emission lines. If this is the case, the presence of such high-latitude dust grains implies that neutral gaseous matter is also blown out during the course of the superwind activity.

During the course of our deep optical imaging survey for Ly alpha emitters at z approximately 5.7 in the field around the z=5.74 quasar SDSSp J104433.04-012502.2, we have found a candidate strong emission-line source. Follow-up optical spectroscopy shows that the emission line profile of this object is asymmetric, showing excess red-wing emission. These properties are consistent with an identification of Ly alpha emission at a redshift of z=5.687 +/- 0.002. The observed broad line width, Delta V_{FWHM} ~= 340 km s^{-1} and excess red-wing emission also suggest that this object hosts a galactic superwind.

We present new optical spectroscopy of the lens elliptical galax in the ``Einstein Cross'' lens system HST 14176+5226, using the Faint Object Camera and Spectrograph (FOCAS) of the Subaru t Our spectroscopic observations are aimed at measuring the stella dispersion of the lens galaxy, located at high redshift of z_L= as an important component to lens models. We have measured this 230 +- 14 km s^{-1} (1 sigma) inside 0.35 effective radi based on the comparison between the observed galaxy spectrum and templates of three G-K giants by means of the Fourier cross-corr To extract the significance of this information on the geometry universe which also affects the lensing of the background image, to fit three different lens models to the available data of the Provided that the lens galaxy has the structural and dynamical p (i.e., its radial density profile, core radius, and velocity ani similar to those of local elliptical galaxies, we calculate the function for the simultaneous reproduction of both the observed and newly measured velocity dispersion of the lens. Although the interval depends rather sensitively on the adopted lens models o parameters, our experiments suggest the larger likelihood for a cosmological constant, Omega_Lambda: formal 1 sigma lower Omega_Lambda in the flat universe ranges 0.73 to 0.97, where lower limit is basically unavailable. This method for determinin model is thus dependent on lens models but is insensitive to oth ambiguities, such as the dust absorption or the evolutionary eff galaxies. Exploring spectroscopic observations of more lens gala redshift may minimize the model uncertainties and thus place a m constraint on Omega_Lambda.

We discuss the properties of galaxies around the radio galaxy 3C 324 at z=1.2 based on BVRIK' multi-band imaging data. We have applied a photometric-redshift technique to objects in the 3C 324 field, and identified 35 objects as plausible cluster members. We have found that red and luminous members are concentrated in a small region enclosed by a circle of 40" radius (0.33 Mpc at z=1.2 for \Omega_0=0.3, \lambda_0=0.7, H_0=70km/s/Mpc cosmology) from the 3C 324 galaxy. The 3C 324 cluster is probably much more compact in size compared with the local clusters. We constructed a K'-band luminosity function of the cluster members and fit a Schechter function, and found the characteristic magnitude to be K*_{AB}=20.2+-0.6. This value is consistent with the extrapolation of the pure passive evolution seen for z<1 clusters. We have identified eight bright galaxies which form a red color-magnitude sequence. The slope of the sequence is consistent with the passive evolution model down to K'_{AB}<22; we also found that there is no clear age variation in these bright red galaxies. However, seven out of these eight galaxies exhibit a significant excess in the rest UV light with respect to the passive evolution model. This may suggest that the massive early-type galaxies in this high-redshift cluster are still forming stars to some extent. We have confirmed a truncation of the color-magnitude sequence at K'_{AB}\sim22; faint passively-evolving galaxies may not yet be present in this cluster at z\sim1.2. The overall color distribution of the cluster members, selected by the photometric redshift technique, is found to be very broad. We derived the fraction of blue galaxies in this cluster following a

We present new optical spectroscopy of the high-redshift powerful radio galaxy MRC 0406$-$244 at redshift of 2.429. We find that the two extensions toward NW and SE probed in the rest-frame ultraviolet image are heated mainly by the nonthermal continuum of the active galactic nucleus. However, each extension shows a shell-like morphology, suggesting that they are a pair of superbubbles induced by the superwind activity rather than by the interaction between the radio jet and the ambient gas clouds. If this is the case, the intense starburst responsible for the formation of superbubbles could occur $\sim 1 \times 10^9$ yr ago. On the other hand, the age of the radio jets may be of the order of $\sim 10^6$ yr, being much shorter than the starburst age. Therefore, the two events, i.e., the starburst and the radio-jet activities, are independent phenomena. However, their directions of the expanding motions could be governed by the rotational motion of the gaseous component in the host galaxy. This idea appears to explain the alignment effect of MRC 0406$-$244.

Deep near-infrared images of a blank 2'x2' section of sky near the Galactic north pole taken by Subaru Telescope are presented. The total integration times of the J and K' bands are 12.1 hours and 9.7 hours, resulting in 5-sigma limiting magnitudes of 25.1 and 23.5 mag, respectively. The numbers of sources within these limiting magnitudes found with an automated detection procedure are 385 in the J band and 350 in K'. Based on photometric measurements of these sources, we present number count vs. magnitude relations, color vs. magnitude diagrams, size vs. color relationships, etc. The slope of the galaxy number count plotted against the AB magnitude scale is about 0.23 in the 22 to 26 AB magnitude range of both bands. The spatial number density of galaxies as well as the slopes in the faint-end region given by the Subaru Deep Field (SDF) survey is consistent with those given by HST-NICMOS surveys as expressed on the AB magnitude diagram. Several sources having very large J-K' color are found including a few K' objects without detection at J. In addition, a number of faint Galactic stars are also detected, most of which are assigned to M-subdwarfs, together with a few brown dwarf candidates.

We have obtained high-resolution K'-band images of the powerful z=1.206 radio galaxy 3C 324 with the Subaru telescope under seeing conditions of 0.3--0.4 arcsec. We clearly resolved the galaxy and directly compared it to the optical images obtained with the Hubble Space Telescope. The host galaxy of 3C 324 is revealed to be a moderately luminous elliptical galaxy with a smooth light profile. The effective radius of the galaxy, as determined by profile fitting, is 1.3+-0.1 arcsec (1.2 kpc), which is significantly smaller than the value of 2.2 arcsec, published in Best et al. (1998, MNRAS, 292, 758). The peak of the K'-band light coincides with the position of the radio core, which implies that the powerful AGN lies at the nucleus of the host galaxy. The peak also coincides with the gap in the optical knotty structures which may be a dust lane hiding the UV-optical emission of the AGN from our line of sight; it is very likely that we are seeing the obscuring structure almost edge-on. We clearly detected the `aligned component' in the K'-band image by subtracting a model elliptical galaxy from the observed image. The red R_F702W-K color of the outer region of the galaxy avoiding the aligned component indicates that the near infrared light of the host galaxy is dominated by an old stellar population.

We have investigated the optical and near-infrared colors of K'-selected galaxies in clusters at z ~ 1.2 near to the radio galaxy 3C 324 using images obtained with the Subaru telescope and archival HST data. The distribution of colors of the galaxies in the cluster region is found to be fairly broad, and it may imply significant scatter in their star-formation histories, although the effect of contamination of field galaxies is uncertain. The red sequence of galaxies whose R-K colors are consistent with passive evolution models for old galaxies is found to be truncated at K' ~ 20 mag, and there are few fainter galaxies with similar red colors in the cluster region. We find that the bulge-dominated galaxies selected by quantitative morphological classification form a broad sequence in the color-magnitude diagram, whose slope is much steeper than that expected from metallicity variations within a passively evolving coeval galaxy population. We argue that the observed color-magnitude sequence can be explained by metallicity and age variations, and the fainter galaxies with K' > 20 mag may be 1-2 Gyr younger than the brighter galaxies. Some spatial segregation of the color and K'-band luminosity is seen in the sky distribution; the redder and the brighter objects tend to be located near 3C 324.

We investigate the $K$-band luminosity distribution of galaxies in the region of clusters at $z\sim1.2$ near to the radio galaxy 3C 324. The imaging data were obtained during the commissioning period of the Subaru telescope. There is a significant excess of the surface number density of the galaxies with $K =$ 17--20 mag in the region within $\sim$ 40'' from 3C 324. At this bright end, the measured luminosity distribution shows a drop, which can be represented by the exponential cut off of the Schechter-function formula; the best-fitted value of the characteristic magnitude, $K^{*}$, is $\sim 18.4\pm0.8$. This measurement follows the evolutionary trend of the $K^*$ of the rich clusters observed at an intermediate redshift, which is consistent with passive evolution models with a formation redshift $z_f \gtsim 2$. At $K \gtsim 20$ mag, however, the excess of the galaxy surface density in the region of the clusters decreases abruptly, which may imply that the luminosity function of the cluster galaxies has a negative slope at the faint end. This may imply strong luminosity segregation between the inner and outer parts of the clusters, or some deficit of faint galaxies in the cluster central region of the cluster.

We have performed a long-slit K band spectroscopic observation of the luminous infrared galaxy NGC 6240. The peak position of the H_2 v=1-0 S(1) emission in the slit is located ~0.3" - 0.4" north of the southern nucleus. It is almost the midpoint between the southern nucleus and the peak position of the ^12CO J=1-0 emission. Based on the line-ratio analyses, we suggest the excitation mechanism of H_2 is pure thermal at most positions. In the southern region we find the following three velocity components in the H_2 emission: the blueshifted shell component (~-250 km s^-1 with respect to V_sys) which is recognized as a distinct C-shape distortion in the velocity field around the southern nucleus, the high-velocity blueshifted ``wing'' component (~-1000 km s^-1 with respect to V_sys), and the component indicating possible line splitting of ~500 km s^-1. The latter two components are extended to the south from the southern nucleus. We show that the kinematic properties of these three components can be reproduced by expanding motion of a shell-like structure around the southern nucleus. The offset peak position of the H_2 emission can be understood if we assume that the shell expanding to the north interacts with the extragalactic molecular gas. At the interface between the shell and the molecular gas concentration the cloud-crushing mechanism proposed by Cowie et al. (1981) may work efficiently, and the intense H_2 emission is thus expected there. All these findings lead us to propose a model that the most H_2 emission is attributed to the shock excitation driven by the superwind activity of the southern nucleus.

We report on infrared imaging observations of the z=2.429 radio galaxy B3 0731+438 with the Subaru telescope. The images were taken with the K'-band filter and the 2.25 um narrow-band filter to examine the structure and properties of the Ha+[N II] 6548,6583 emission-line components. The Ha+[N II] emission-line image shows biconical lobes with an extent of 40 kpc, which are aligned with the radio axis. The rest-frame equivalent widths of the emission lines at these cones are as large as 1100 AA, and can be well explained by a gas-cloud model photoionized by power-law continuum radiation. The isotropic ionizing photon luminosity necessary to ionize the hydrogen gas in these cones amounts to 1e57(photons/s), which is larger than that in the majority of radio-loud QSOs. From these results, we propose that the Ha alignment effect in this object is produced by biconical gas clouds, which are swept up by the passage of radio jets, and are ionized by strong UV radiation from a hidden AGN. The continuum image consists of two components, a stellar-like point source and an extended diffuse galaxy. These are supposed to be a type-2 AGN and its host galaxy. The SED is fitted by a combination of spectra of a reddened dust-scattered AGN and an instantaneous starburst population of 500 Myr old. The stellar mass of the galaxy is estimated to be 3e11 M_solar, which is as large as that of typical 3C radio galaxies at z=1.