To date no metal-free stars have been identified by direct observations. The most common method of constraining their properties is searching the spectra of the most metal-poor stars for the chemical elements created in the first stars and their supernova. In this approach, modelled supernova yields are compared to the observed abundance patterns in extremely metal-poor stars. The method typically only uses the abundance ratios, i.e., the yields are diluted to the observed level. Following the usual assumption of spherical symmetry we compute a simple lower limit of the mass a supernova can mix with and find that it is consistent with all published simulations of early chemical enrichment in the interstellar medium. For three different cases, we demonstrate that this dilution limit can change the conclusions from the abundance fitting. There is a large discrepancy between the dilution found in simulations of SN explosions in minihaloes and the dilution assumed in many abundance fits. Limiting the dilution can significantly alter the likelihood of which supernovae are possible progenitors of observed CEMP-no stars. In particular, some of the faint, very low-yield SNe, which have been suggested as models for the abundance pattern of SMSS0313-6708, cannot explain the measured metal abundances, as their predicted metal yields are too small by two orders of magnitude. Altogether, the new dilution model presented here emphasizes the need to better understand the mixing and dilution behaviour of aspherical SNe.

We combine asteroseismology, optical high-resolution spectroscopy, and kinematic analysis for 26 halo red giant branch stars in the \textit{Kepler} field in the range of $-2.5<[\mathrm{{Fe}/{H}}]<-0.6$. After applying theoretically motivated corrections to the seismic scaling relations, we obtain an average mass of $0.96\pm 0.02\,\mathrm{M_{\odot}}$ for our sample of halo stars. This maps into an age of $\sim 8\,\mathrm{Gyr}$, significantly younger than independent age estimates of the Milky Way stellar halo. Although our results suggest that further improvements are needed to accurately determine ages of metal-poor stars, asteroseismology still provides robust relative ages. The lack of a mass dispersion among our sample constrains a relative age dispersion to $<18\%$, corresponding to $<2\,\mathrm{Gyr}$. The precise chemical abundances allow us to separate most of the stars with [{Fe}/{H}]$>-1.7$ into two [{Mg}/{Fe}] groups. While [$\alpha$/{Fe}] and [{Eu}/{Mg}] ratios are different between the two subsamples, [$s$/Eu], where $s$ stands for Ba, La, Ce, and Nd, does not show a significant difference. These abundance ratios suggest that the chemical evolution of the low-Mg population is contributed by type~Ia supernovae, but not by low-to-intermediate mass asymptotic giant branch stars, providing a constraint on its star formation timescale as $100\,\mathrm{Myr}<\tau<300\,\mathrm{Myr}$. We also do not detect any significant mass difference between the two [{Mg}/{Fe}] groups, thus suggesting that their formation epochs are not separated by more than 1.5 Gyr.

We perform an extensive review of the numerous studies and methods used to determine the total mass of the Milky Way. We group the various methods into seven broad classes, including: i) estimating Galactic escape velocity using high velocity objects; ii) measuring the rotation curve through terminal and circular velocities; iii) modeling halo stars, globular clusters and satellite galaxies with the Spherical Jeans equation and iv) with phase-space distribution functions; v) simulating and modeling the dynamics of stellar streams and their progenitors; vi) modeling the motion of the Milky Way, M31 and other distant satellites under the framework of Local Group timing argument; and vii) measurements made by linking the brightest Galactic satellites to their counterparts in simulations. For each class of methods, we introduce their theoretical and observational background, the method itself, the sample of available tracer objects, model assumptions, uncertainties, limits and the corresponding measurements that have been achieved in the past. Both the measured total masses within the radial range probed by tracer objects and the extrapolated virial masses are discussed and quoted. We also discuss the role of modern numerical simulations in terms of helping to validate model assumptions, understanding systematic uncertainties and calibrating the measurements. While measurements in the last two decades show a factor of two scatters, recent measurements using \textit{Gaia} DR2 data are approaching a higher precision. We end with a detailed discussion of future developments, especially as the size and quality of the observational data will increase tremendously with current and future surveys. In such cases, the systematic uncertainties will be dominant and thus will necessitate a much more rigorous testing and characterization of the various mass determination methods.

We present proper motion measurements for more than $0.55$ million main-sequence stars, by comparing astrometric positions of matched stars between the multi-band imaging datasets from the Hyper Suprime-Cam (HSC) Survey and the SDSS Stripe 82. In doing this we use $3$ million galaxies to recalibrate the astrometry and set up a common reference frame between the two catalogues. The exquisite depth and the nearly $12$ years of time baseline between HSC and SDSS enable high-precision measurements of statistical proper motions for stars down to $i\simeq 24$. A validation of our method is demonstrated by the agreement with the $Gaia$ proper motions, to the precision better than $0.1$ mas yr$^{-1}$. To retain the precision, we make a correction of the subtle effects due to the differential chromatic refraction in the SDSS images based on the comparison with the $Gaia$ proper motions against colour of stars, which is validated using the SDSS spectroscopic quasars. Combining with the photometric distance estimates for individual stars based on the precise HSC photometry, we show a significant detection of the net proper motions for stars in each bin of distance out to $100$ kpc. The two-component tangential velocities after subtracting the apparent motions due to our own motion display rich phase-space structures including a clear signature of the Sagittarius stream in the halo region of distance range $[10,\ 35]$ kpc. We also measure the tangential velocity dispersion in the distance range $5-20$ kpc and find that the data are consistent with a constant isotropic dispersion of $80\pm 10\ {\rm km/s}$. More distant stars appear to have random motions with respect to the Galactic centre on average.

Dwarf spheroidal galaxies (dSphs) are promising targets of indirect detection experiments searching for dark matter (DM) at present universe. Toward robust prediction for the amount of signal flux originating in DM annihilation inside dSphs, a precise determination of DM distributions as well as $J$-factors of the dSphs is particularly important. In this work, we estimate those of Draco, Sculptor, and Ursa Minor dSphs by an improved statistical method in which both foreground stars and dSph member stars are simultaneously taken into account. We define the likelihood function of the method as the so-called conditional one to remove sampling bias of observed stellar data. This improved method enables us to estimate DM distributions and $J$-factors of the dSphs directly from observed stellar data contaminated by foreground stars without imposing stringent membership criteria on the measured quantities.

The nature of the first massive stars may be inferred by investigating the origin of the extremely metal-poor (EMP) stars, likely formed from the ejecta of one or a few previous massive stars. We investigate the rotational properties of early massive stars by comparing the abundance patterns of EMP stars with rotating massive stellar models. Low metallicity 20 $M_{\odot}$ stellar models with initial rotation rates between 0 and $70~\%$ of the critical velocity are computed. Explosions with strong fallback are assumed. The ejected material is considered to fit individually the abundance patterns of 272 EMP stars with $-4<$ [Fe/H] $<-3$. With increasing initial rotation, the [C/H], [N/H], [O/H], [Na/H], [Mg/H] and [Al/H] ratios in the massive star ejecta are gradually increased. Among the 272 EMP stars considered, $\sim 40-50~\%$ are consistent with our models. About $60 - 70~\%$ of the CEMP star sample is reproduced against $\sim 20 - 30~\%$ for the C-normal EMP star sample. The CEMP stars are preferentially reproduced with a material coming from mid to fast rotating massive stars. The velocity distribution derived from the best massive star models increases from no rotation to fast rotation. The maximum is reached for massive stars having initial equatorial velocities of $\sim 550 - 640$ km~s$^{-1}$. Although subject to significant uncertainties, these results suggest that the rotational mixing operating in between the H-burning shell and the He-burning core of early massive stars played an important role in the early chemical enrichment of the Universe. The comparison of the velocity distribution derived from the best massive star models with velocity distributions of nearby OB stars suggests a greater amount of massive fast rotators in the early Universe. This may have important consequences for reionization or integrated light from high redshift galaxies.

The Milky Way is a unique laboratory, where stellar properties can be measured and analyzed in detail. In particular, stars in the older populations encode information on the mechanisms that led to the formation of our Galaxy. In this article, we analyze the kinematics, spatial distribution, and chemistry of a large number of stars in the Solar Neighborhood, where all of the main Galactic components are well-represented. We find that the thick disk comprises two distinct and overlapping stellar populations, with different kinematic properties and chemical compositions. The metal-weak thick disk (MWTD) contains two times less metal content than the canonical thick disk, and exhibits enrichment of light elements typical of the oldest stellar populations of the Galaxy. The rotational velocity of the MWTD around the Galactic center is ~ 150 km s^(-1), corresponding to a rotational lag of 30 km s^(-1) relative to the canonical thick disk (~ 180 km s^(-1)), with a velocity dispersion of 60 km s^(-1). This stellar population likely originated from the merger of a dwarf galaxy during the early phases of our Galaxy's assembly, or it is a precursor disk, formed in the inner Galaxy and brought into the Solar Neighborhood by bar instability or spiral-arm formation mechanisms.

I would like to review recent efforts of detailed chemical abundance measurements for field Milky Way halo stars. Thanks to the advent of wide-field spectroscopic surveys up to a several kpc from the Sun, large samples of field halo stars with detailed chemical measurements are continuously expanding. Combination of the chemical information and full six dimensional phase-space information is now recognized as a powerful tool to identify cosmological accretion events that have built a sizable fraction of the present-day stellar halo. Future observational prospects with wide-field spectroscopic surveys and theoretical prospects with supernova nucleosynthetic yields are also discussed.

Geodetic Very Long Baseline Interferometry (VLBI) data are capable of measuring the light deflection caused by the gravitational field of the Sun and large planets with high accuracy. The parameter $\gamma$ of the parametrized Post-Newtonian (PPN) formalism estimated using observations of reference radio sources near the Sun should be equal to unity in the general relativity. We have run several VLBI experiments tracking reference radio sources from 1 to 3 degrees from the Sun. The best formal accuracy of the parameter $\gamma$ achieved in the single-session mode is less than 0.01 percent, or better than the formal accuracy obtained with a global solution included all available observations at arbitrary elongation from the Sun. We are planning more experiments starting from 2020 using better observing conditions near the minimum of the Solar activity cycle.

We report on the discovery of a new Milky Way (MW) satellite in Bo\"otes based on data from the on-going Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP). This satellite, named Bo\"otes IV, is the third ultra-faint dwarf that we have discovered in the HSC-SSP. We have identified a statistically significant (32.3$\sigma$) overdensity of stars having characteristics of a metal-poor, old stellar population. The distance to this stellar system is $D_{\odot}=209^{+20}_{-18}$ kpc with a $V$-band absolute magnitude of $M_V=-4.53^{+0.23}_{-0.21}$ mag. Bo\"otes IV has a half-light radius of $r_h=462^{+98}_{-84}$ pc and an ellipticity of $0.64^{+0.05}_{-0.05}$, which clearly suggests that this is a dwarf satellite galaxy. We also found another overdensity that appears to be a faint globular cluster with $M_V=-0.20^{+0.59}_{-0.83}$ mag and $r_h=5.9^{+1.5}_{-1.3}$ pc located at $D_{\odot}=46^{+4}_{-4}$ kpc. Adopting the recent prediction for the total population of satellites in a MW-sized halo by Newton et al. (2018), which combined the characteristics of the observed satellites by SDSS and DES with the subhalos obtained in $\Lambda$CDM models, we estimate that there should be about two MW satellites at $M_V\le0$ in the $\sim676$ deg$^2$ covered by HSC-SSP, whereas that area includes six satellites. Thus, the observed number of satellites is larger than the theoretical prediction. On the face of it, we have a problem of too many satellites, instead of the well-known missing satellites problem whereby the $\Lambda$CDM theory overpredicts the number of satellites in a MW-sized halo. This may imply that the models need more refinements for the assignment of subhalos to satellites such as considering those found by the current deeper survey. [abridged]

Small stellar systems like dwarf galaxies are suggested to be the main building blocks of our Galaxy by numerical simulations in Lambda CDM models. The existence of star streams like Sagittarius tidal stream indicates that dwarf galaxies play a role in the formation of the Milky Way. However, it is unclear how many and what kind of stars in our Galaxy are originated from satellite dwarf galaxies, which could be constrained by chemical abundances of metal-poor stars. Here we report on the discovery of a metal-poor star with an extreme r-process enhancement and alpha-element deficiency. In this star, the abundance ratio of the r-process element Eu with respect to Fe is more than one order of magnitude higher than the Sun and the metallicity is 1/20 of the solar one. Such kind of stars have been found in present-day dwarf galaxies, providing the clearest chemical signature of past accretion events. The long timescale of chemical evolution of the host dwarf galaxy expected from the abundance of alpha element with respect to Fe suggests that the accretion occurred in a relatively late phase compared to most of the accretions that formed the bulk of the Milky Way halo.

We present observational evidence that an aspherical supernova explosion could have occurred in the First stars in the early universe. Our results are based on the First determination of a Zn abundance in a Hubble Space Telescope/Cosmic Origins Spectrograph high-resolution UV spectrum of a hyper metal-poor (HMP) star, HE1327-2326, with [Fe/H](NLTE) = -5.2. We determine [Zn/Fe] = 0.80$\pm$0.25 from a UV Zn I line at 2138 detected at $3.4\sigma$. Yields of a 25M$_{\odot}$ aspherical supernova model with artificially modified densities exploding with E = 5x10$^{51}$ ergs best match the entire abundance pattern of HE1327-2326. Such high-entropy hypernova explosions are expected to produce bipolar outfows which could facilitate the external enrichment of small neighboring galaxies. This has already been predicted by theoretical studies of the earliest star forming minihalos. Such a scenario would have significant implications for the chemical enrichment across the early universe as HMP Carbon Enhanced Metal-Poor (CEMP) stars such as HE1327-2326 might have formed in such externally enriched environments.

High-resolution spectroscopy of the core of the Perseus Cluster of galaxies, using the $Hitomi$ satellite above 2 keV and the $XMM$-$Newton$ Reflection Grating Spectrometer at lower energies, provides reliable constraints on the abundances of O, Ne, Mg, Si, S, Ar, Ca, Cr, Mn, Fe, and Ni. Accounting for all known systematic uncertainties, the Ar/Fe, Ca/Fe, and Ni/Fe ratios are determined with a remarkable precision of less than 10%, while the constraints on Si/Fe, S/Fe, and Cr/Fe are at the 15% level, and Mn/Fe is measured with a 20% uncertainty. The average biases in determining the chemical composition using archival CCD spectra from $XMM$-$Newton$ and $Suzaku$ range typically from 15-40%. A simple model in which the enrichment pattern in the Perseus Cluster core and the proto-solar nebula are identical gives a surprisingly good description of the high-resolution X-ray spectroscopy results, with $\chi^2=10.7$ for 10 d.o.f. However, this pattern is challenging to reproduce with linear combinations of existing supernova nucleosynthesis calculations, particularly given the precise measurements of intermediate $\alpha$-elements enabled by $Hitomi$. We discuss in detail the degeneracies between various supernova progenitor models and explosion mechanisms, and the remaining uncertainties in these theoretical models. We suggest that including neutrino physics in the core-collapse supernova yield calculations may improve the agreement with the observed pattern of $\alpha$-elements in the Perseus Cluster core. Our results provide a complementary benchmark for testing future nucleosynthesis calculations required to understand the origin of chemical elements.

Extremely metal poor (EMP) stars in the Milky Way inherited the chemical composition of the gas out of which they formed. They therefore carry the chemical fingerprint of the first stars in their spectral lines. It is commonly assumed that EMP stars form from gas that was enriched by only one progenitor supernova ('mono-enriched'). However, recent numerical simulations show that the first stars form in small clusters. Consequently, we expect several supernovae to contribute to the abundances of an EMP star ('multi-enriched'). We analyse seven recently observed EMP stars from the TOPoS survey by applying the divergence of the chemical displacement and find that J1035+0641 is mono-enriched ($p_{mono}=53\%$) and J1507+0051 is multi-enriched ($p_{mono}=4\%$). For the remaining five stars we can not make a distinct prediction ($p_{mono} \lesssim 50\%$) due to theoretical and observational uncertainties. Further observations in the near-UV will help to improve our diagnostic and therefore contribute to constrain the nature of the first stars.

We present a study of the high-resolution optical spectrum for the hot post-asymptotic giant branch (post-AGB) star, Hen 3-1013, identified as the optical counterpart of the infrared sourceIRAS 14331-6435. For the first time the detailed identifications of the observed absorption and emission features in the wavelength range 3700-9000 \AA\ is carried out. Absorption lines of HI, HeI, CI, NI, OI, NeI, CII, NII, OII, SiII, SII, ArII, FeII, MnII, CrII, TiII, CoII, NiII, SIII, FeIII and SIV were detected. From the absorption lines, we derived heliocentric radial velocities of $V_r=-29.6\pm0.4$ km/s. We have identified emission permitted lines of OI, NI, FeII, MgII, SiII and AlII. The forbidden lines of [NI], [FeII], [CrII] and [NiII] have been identified also. Analysis of [NiII] lines in the gaseous shell gives an estimate for the electron density $N_e\sim10^7$ cm$^{-3}$ and the expansion velocity of the nebula $V_{exp}=12$ km/s. The mean radial velocity as measured from emission features of the envelope is $V_r=-36.0\pm0.4$ km/s. The Balmer lines H$\alpha$, H$\beta$ and H$\gamma$ show P Cyg behaviour which indicate ongoing post-AGB mass-loss. Based on ASAS and ASAS-SN data, we have detected rapid photometric variability in Hen 3-1013 with an amplitude up to 0.2 mag in the V band. The star's low-resolution spectrum underwent no significant changes from 1994 to 2012. Based on archival data, we have traced the photometric history of the star over more than 100 years. No significant changes in the star brightness have been found.

Dwarf spheroidal galaxies (dSphs) are promising targets for the gamma-ray dark matter (DM) search. In particular, DM annihilation signal is expected to be strong in some of the recently discovered nearby ultra-faint dSphs, which potentially give stringent constraints on the $\mathcal{O}(1)$ TeV WIMP DM. However, various non-negligible systematic uncertainties complicate the estimation of the astrophysical factors relevant for the DM search in these objects. Among them, the effects of foreground stars particularly attract attention because the contamination is unavoidable even for the future kinematical survey. In this article, we assess the effects of the foreground contamination on the astrophysical $J$-factor estimation by generating mock samples of stars in the four ultra-faint dSphs and using a model of future spectrographs. We investigate various data cuts to optimize the quality of the data and apply a likelihood analysis which takes member and foreground stellar distributions into account. We show that the foreground star contaminations in the signal region (the region of interest) and their statistical uncertainty can be estimated by interpolating the foreground star distribution in the control region where the foreground stars dominate the member stars. Such regions can be secured at future spectroscopic observations utilizing a multiple object spectrograph with a large field of view; e.g. the Prime Focus Spectrograph mounted on Subaru Telescope. The above estimation has several advantages: The data-driven estimation of the contamination makes the analysis of the astrophysical factor stable against the complicated foreground distribution. Besides, foreground contamination effect is considered in the likelihood analysis.

We report on the elemental abundances of the carbon-enhanced metal-poor (CEMP) star J2217+2104 discovered by our metal-poor star survey with LAMOST and Subaru. This object is a red giant having extremely low Fe abundance ([Fe/H]=-4.0) and very large enhancement of C, N, and O with excesses of Na, Mg, Al, and Si. This star is a new example of a small group of such CEMP stars identified by previous studies. We find a very similar abundance pattern for O-Zn in this class of objects that shows enhancement of elements up to Si and normal abundance of Ca and Fe-group elements. Whereas the C/N ratio is different among these stars, the (C+N)/O ratio is similar. This suggests that C was also yielded with similar abundance ratios relative to O-Zn in progenitors, and was later affected by the CN-cycle. By contrast, the heavy neutron-capture elements Sr and Ba are deficient in J2217+2104, compared to the four objects in this class previously studied. This indicates that the neutron-capture process in the early Galaxy, presumably the r-process, has no direct connection to the phenomenon that has formed such CEMP stars. Comparisons of the abundance pattern well determined for such CEMP stars with those of supernova nucleosynthesis models constrain the progenitor mass to be about 25Msun, which is not particularly different from typical mass of progenitors expected for extremely metal-poor stars in general.

Extremely metal-poor (EMP) stars in the Milky Way (MW) allow us to infer the properties of their progenitors by comparing their chemical composition to the metal yields of the first supernovae. This method is most powerful when applied to mono-enriched stars, i.e. stars that formed from gas that was enriched by only one previous supernova. We present a novel diagnostic to identify this subclass of EMP stars. We model the first generations of star formation semi-analytically, based on dark matter halo merger trees that yield MW-like halos at the present day. Radiative and chemical feedback are included self-consistently and we trace all elements up to zinc. Mono-enriched stars account for only $\sim 1\%$ of second generation stars in our fiducial model and we provide an analytical formula for this probability. We also present a novel analytical diagnostic to identify mono-enriched stars, based on the metal yields of the first supernovae. This new diagnostic allows us to derive our main results independently from the specific assumptions made regarding Pop III star formation, and we apply it to a set of observed EMP stars to demonstrate its strengths and limitations. Our results may provide selection criteria for current and future surveys and therefore contribute to a deeper understanding of EMP stars and their progenitors.

We report chemical abundances of 14 young $\alpha$-rich stars including neutron-capture elements based on high-quality optical spectra from HIRES/Keck I and differential line-by-line analysis. From the comparison of the abundance patterns of young $\alpha$-rich stars to those of nearby bright red giants with a similar metallicity range ($-0.7<[\mathrm{Fe/H}]<-0.2$), we confirm their high $\alpha$-element abundances reported by previous studies based on near-infrared spectroscopy. We reveal for the first time low abundances of $s$-process elements and high abundances of $r$-process elements. All the abundances are consistent with those seen in typical $\alpha$-rich population of the Galactic disk, and no abundance anomalies are found except for Li-enhancement in one object previously reported and mild enhancement of Na in two stars. In particular, the lack of $s$-process enhancement excludes the hypothesis that mass transfer from asymptotic giant branch stars plays an important role in the formation of young $\alpha$-rich stars. The high frequency of radial velocity variation (more than $50\%$) is also confirmed. We argue that mass transfer from low-mass red giants is the likely dominant formation mechanism for young $\alpha$-rich stars.

We compare elemental abundance patterns of $\sim 200$ extremely metal-poor (EMP; [Fe/H]$<-3$) stars with supernova yields of metal-free stars in order to obtain insights into the characteristic masses of the first (Population III or Pop III) stars in the Universe. Supernova yields are prepared with nucleosynthesis calculations of metal-free stars with various initial masses ($M=$13, 15, 25, 40 and 100 $M_{\odot}$) and explosion energies ($E_{51}=E/10^{51}$[erg]$=0.5-60$) to include low-energy, normal-energy, and high-energy explosions. We adopt the mixing-fallback model to take into account possible asymmetry in the supernova explosions and the yields that best-fit the observed abundance patterns of the EMP stars are searched by varying the model parameters. We find that the abundance patterns of the EMP stars are predominantly best-fitted with the supernova yields with initial masses $M<40 M_{\odot}$, and that more than than half of the stars are best fitted with the $M=25 M_\odot$ hypernova ($E_{51}=10$) models. The results also indicate that the majority of the primordial supernovae have ejected $10^{-2}-10^{-1} M_\odot$ of $^{56}$Ni leaving behind a compact remnant, either a neutron star or a black hole, with mass in a range of $\sim 1.5-5 M_{\odot}$. The results suggest that the masses of the first stars responsible for the first metal-enrichment are predominantly $< 40 M_{\odot}$. This implies that the higher mass first stars were either less abundant or directly collapsing into a blackhole without ejecting heavy elements or that a supernova explosion of a higher-mass first star inhibits the formation of the next generation of low-mass stars at [Fe/H]$<-3$.

We construct $z\sim6-7$, 8, and 9 faint Lyman break galaxy samples (334, 61, and 37 galaxies, respectively) with accurate size measurements with the software $\texttt{glafic}$ from the complete Hubble Frontier Fields (FF) cluster and parallel fields data. These are the largest samples hitherto and reach down to the faint ends of recently obtained deep luminosity functions. At faint magnitudes, however, these samples are highly incomplete for galaxies with large sizes, implying that derivation of the luminosity function (LF) sensitively depends on the intrinsic size--luminosity (RL) relation. We thus conduct simultaneous maximum-likelihood estimation of LF and RL relation parameters from the observed distribution of galaxies on the RL plane with help of a completeness map as a function of size and luminosity. At $z\sim6-7$, we find that the intrinsic RL relation expressed as $r_\textrm{e} \propto L^\beta$ has a notably steeper slope of $\beta=0.46^{+0.08}_{-0.09}$ than those at lower redshifts, which in turn implies that the LF has a relatively shallow faint-end slope of $\alpha=-1.86^{+0.17}_{-0.18}$. This steep $\beta$ can be reproduced by a simple analytical model in which smaller galaxies have lower specific angular momenta. The $\beta$ and $\alpha$ values for the $z\sim8$ and 9 samples are consistent with those for $z\sim6-7$ but with larger errors. For all three samples there is a large, positive covariance between $\beta$ and $\alpha$, implying that the simultaneous determination of these two parameters is important. We also provide new strong lens mass models of Abell S1063 and Abell 370 as well as updated mass models of Abell 2744 and MACS J0416.1$-$2403.

We present UV luminosity functions of dropout galaxies at $z\sim6-10$ with the complete Hubble Frontier Fields data. We obtain a catalog of $\sim450$ dropout-galaxy candidates (350, 66, and 40 at $z\sim6-7$, 8, and 9, respectively), whose UV absolute magnitudes reach $\sim-14$ mag, $\sim2$ mag deeper than the Hubble Ultra Deep Field detection limits. We carefully evaluate number densities of the dropout galaxies by Monte-Carlo simulations, including all lensing effects such as magnification, distortion, and multiplication of images as well as detection completeness and contamination effects in a self-consistent manner. We find that UV luminosity functions at $z\sim6-8$ have steep faint-end slopes, $\alpha\sim-2$, and likely steeper slopes, $\alpha\lesssim-2$ at $z\sim9-10$. We also find that the evolution of UV luminosity densities shows a non-accelerated decline beyond $z\sim8$ in the case of $M_\mathrm{trunc}=-15$, while an accelerated in the case of $M_\mathrm{trunc}=-17$. We examine whether our results are consistent with the Thomson scattering optical depth from the Planck satellite and the ionized hydrogen fraction $Q_\mathrm{HII}$ at $z\lesssim7$ based on the standard analytic reionization model. We find that there exist reionizaiton scenarios that consistently explain all the observational measurements with the allowed parameters of $f_{\rm esc}=0.17^{+0.07}_{-0.03}$ and $M_\mathrm{trunc}>-14.0$ for $\log\xi_\mathrm{ion}/[\mathrm{erg}^{-1}\ \mathrm{Hz}]=25.34$, where $f_{\rm esc}$ is the escape fraction, $M_\mathrm{trunc}$ is the faint limit of the UV luminosity function, and $\xi_\mathrm{ion}$ is the conversion factor of the UV luminosity to the ionizing photon emission rate. The length of the reionization period is estimated to be $\Delta z=3.9^{+2.0}_{-1.6}$ (for $0.1<Q_\mathrm{HII}<0.99$), consistent with the recent estimate from Planck.

The stellar record of elemental abundances in satellite galaxies is important to identify the origin of r-process because such a small stellar system could have hosted a single r-process event, which would distinguish member stars that are formed before and after the event through the evidence of a considerable difference in the abundances of r-process elements, as found in the ultra-faint dwarf galaxy Reticulum II (Ret II). However, the limited mass of these systems prevents us from collecting information from a sufficient number of stars in individual satellites. Hence, it remains unclear whether the discovery of a remarkable r-process enrichment event in Ret II explains the nature of r-process abundances or is an exception. We perform high-resolution spectroscopic measurements of r-process abundances for twelve metal-poor stars in the Draco dwarf galaxy in the metallicity range of -2.5<[Fe/H]<-2. We found that these stars are separated into two groups with r-process abundances differing by one order of magnitude. A group of stars with high abundances of r-process elements was formed by a single r-process event that corresponds to the event evidenced in Ret II. On the other hand, the low r-process abundance group was formed by another sporadic enrichment channel producing a far fewer r-process elements, which is clearly identified for the first time. Accordingly, we identified two populations of stars with different r-process abundances, which are built by two r-process events that enriched gases at levels that differ by more than one order of magnitude.

We present the results from a search for new Milky Way (MW) satellites from the first two years of data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) $\sim 300$~deg$^2$ and report the discovery of a highly compelling ultra-faint dwarf galaxy candidate in Cetus. This is the second ultra-faint dwarf we have discovered after Virgo~I reported in our previous paper. This satellite, Cetus~III, has been identified as a statistically significant (10.7$\sigma$) spatial overdensity of star-like objects, which are selected from a relevant isochrone filter designed for a metal-poor and old stellar population. This stellar system is located at a heliocentric distance of 251$^{+24}_{-11}$~kpc with a most likely absolute magnitude of $M_V = -2.4 \pm 0.6$~mag estimated from a Monte Carlo analysis. Cetus~III is extended with a half-light radius of $r_h = 90^{+42}_{-17}$~pc, suggesting that this is a faint dwarf satellite in the MW located beyond the detection limit of the Sloan Digital Sky Survey. Further spectroscopic studies are needed to assess the nature of this stellar system. We also revisit and update the parameters for Virgo~I finding $M_V = -0.33^{+0.75}_{-0.87}$~mag and $r_h = 47^{+19}_{-13}$~pc. Using simulations of $\Lambda$-dominated cold dark matter models, we predict that we should find one or two new MW satellites from $\sim 300$~deg$^2$ HSC-SSP data, in rough agreement with the discovery rate so far. The further survey and completion of HSC-SSP over $\sim 1,400$~deg$^2$ will provide robust insights into the missing satellites problem.

The gamma-ray observation of the dwarf spheroidal galaxies (dSphs) is a promising approach to search for the dark matter annihilation (or decay) signal. The dSphs are the nearby satellite galaxies with a clean environment and dense dark matter halo so that they give stringent constraints on the 1 TeV dark matter. However, recent studies have revealed that current estimation of astrophysical factors relevant for the dark matter searches are not conservative, where the various non-negligible systematic uncertainties are not taken into account. Among them, the effect of foreground stars on the astrophysical factors has not been paid much attention, which becomes more important for deeper and wider stellar surveys in the future. In this article, we assess the effects of the foreground contamination by generating the mock samples of stars and using a model of future spectrographs. We investigate various data cuts to optimize the quality of the data and find that the cuts on the velocity and surface gravity can efficiently eliminate the contamination. We also propose a new likelihood function which includes the foreground distribution function. We apply this likelihood function to the fit of the three types of the mock data (Ursa Minor, Draco with large dark matter halo, and Draco with small halo) and three cases of the observation. The likelihood successfully reproduces the input $J$-factor value while the fit without considering the foreground distribution gives large deviation from the input value by a factor of three.

We report the discovery of a new ultra-faint dwarf satellite companion of the Milky Way based on the early survey data from the Hyper Suprime-Cam Subaru Strategic Program. This new satellite, Virgo I, which is located in the constellation of Virgo, has been identified as a statistically significant (5.5 sigma) spatial overdensity of star-like objects with a well-defined main sequence and red giant branch in their color-magnitude diagram. The significance of this overdensity increases to 10.8 sigma when the relevant isochrone filter is adopted for the search. Based on the distribution of the stars around the likely main sequence turn-off at r ~ 24 mag, the distance to Virgo I is estimated as 87 kpc, and its most likely absolute magnitude calculated from a Monte Carlo analysis is M_V = -0.8 +/- 0.9 mag. This stellar system has an extended spatial distribution with a half-light radius of 38 +12/-11 pc, which clearly distinguishes it from a globular cluster with comparable luminosity. Thus, Virgo I is one of the faintest dwarf satellites known and is located beyond the reach of the Sloan Digital Sky Survey. This demonstrates the power of this survey program to identify very faint dwarf satellites. This discovery of VirgoI is based only on about 100 square degrees of data, thus a large number of faint dwarf satellites are likely to exist in the outer halo of the Milky Way.

The dwarf spheroidal galaxies (dSphs) in the Milky Way are the primary targets in the indirect searches for particle dark matter. To set robust constraints on candidate of dark matter particles, understanding the dark halo structure of these systems is of substantial importance. In this paper, we first evaluate the astrophysical factors for dark matter annihilation and decay for 24 dSphs, with taking into account a non-spherical dark halo, using generalized axisymmetric mass models based on axisymmetric Jeans equations. First, from a fitting analysis of the most recent kinematic data available, our axisymmetric mass models are a much better fit than previous spherical ones, thus, our work should be the most realistic and reliable estimator for astrophysical factors. Secondly, we find that among analysed dSphs, the ultra-faint dwarf galaxies Triangulum II and Ursa Major II are the most promising but large uncertain targets for dark matter annihilation while the classical dSph Draco is the most robust and detectable target for dark matter decay. It is also found that the non-sphericity of luminous and dark components influences the estimate of astrophysical factors, even though these factors largely depend on the sample size, the prior range of parameters and the spatial extent of the dark halo. Moreover, owing to these effects, the constraints on the dark matter annihilation cross-section are more conservative than those of previous spherical works. These results are important for optimizing and designing dark matter searches in current and future multi-messenger observations by space and ground-based telescopes.

Gravitational lensing by clusters of galaxies offers a powerful probe of their structure and mass distribution. Deriving a lens magnification map for a galaxy cluster is a classic inversion problem and many methods have been developed over the past two decades to solve it. Several research groups have developed techniques independently to map the predominantly dark matter distribution in cluster lenses. While these methods have all provided remarkably high precision mass maps, particularly with exquisite imaging data from the Hubble Space Telescope (HST), the reconstructions themselves have never been directly compared. In this paper, we report the results of comparing various independent lens modeling techniques employed by individual research groups in the community. Here we present for the first time a detailed and robust comparison of methodologies for fidelity, accuracy and precision. For this collaborative exercise, the lens modeling community was provided simulated cluster images -- of two clusters Ares and Hera -- that mimic the depth and resolution of the ongoing HST Frontier Fields. The results of the submitted reconstructions with the un-blinded true mass profile of these two clusters are presented here. Parametric, free-form and hybrid techniques have been deployed by the participating groups and we detail the strengths and trade-offs in accuracy and systematics that arise for each methodology. We note in conclusion that lensing reconstruction methods produce reliable mass distributions that enable the use of clusters as extremely valuable astrophysical laboratories and cosmological probes.

We present Subaru/FOCAS and Keck/DEIMOS medium-resolution spectroscopy of a tidally disrupting Milky Way (MW) globular cluster Palomar 5 and its tidal stream. The observed fields are located to cover an angular extent of $\sim 17$\arcdeg along the stream, providing an opportunity to investigate a trend in line-of-sight velocities ($V_{\rm los}$) along the stream, which is essential to constrain its orbit and underlying gravitational potential of the Milky Way's dark matter halo. A spectral fitting technique is applied to the observed spectra to obtain stellar parameters and metallicities ([Fe/H]) of the target stars. The 19 stars most likely belonging to the central Pal 5 cluster have a mean $V_{\rm los}$ of $-58.1\pm 0.7$ km s$^{-1}$ and metallicity [Fe/H]$=-1.35\pm 0.06$ dex, both of which are in good agreement with those derived in previous high-resolution spectroscopic studies. Assuming that the stream stars have the same [Fe/H] as the progenitor cluster, the derived [Fe/H] and $V_{\rm los}$ values are used to estimate the possible $V_{\rm los}$ range of the member stars at each location along the stream. Because of the heavy contamination of the field MW stars, estimated $V_{\rm los}$ range depends on prior assumptions about the stream's $V_{\rm los}$, which highlights the importance of more definitely identifying the member stars using proper motion and chemical abundances to obtain an unbiased information of $V_{\rm los}$ in the outer part of the Pal 5 stream. The models for the gravitational potential of the MW's dark matter halo that are compatible with the estimated $V_{\rm los}$ range are discussed.

We present clustering analysis results from 10,381 Lyman break galaxies (LBGs) at z~ 4-7, identified in the Hubble legacy deep imaging and new complimentary large-area Subaru/Hyper Suprime-Cam data. We measure the angular correlation functions (ACFs) of these LBGs at z~4, 5, 6, and 7, and fit these measurements using halo occupation distribution (HOD) models that provide an estimate of halo masses, M_h~(1-20)x10^11 Msun. Our M_h estimates agree with those obtained by previous clustering studies in a UV-magnitude vs. M_h plane, and allow us to calculate stellar-to-halo mass ratios (SHMRs) of LBGs. By comparison with the z~0 SHMR, we identify evolution of the SHMR from z~0 to z~4, and z~4 to z~7 at the >98% confidence levels. The SHMR decreases by a factor of ~2 from z~0 to 4, and increases by a factor of ~4 from z~4 to 7. We compare our SHMRs with results of a hydrodynamic simulation and a semi-analytic model, and find that these theoretical studies do not predict the SHMR increase from z~4 to 7. We obtain the baryon conversion efficiency (BCE) of LBGs at z~4, and find that the BCE increases with increasing dark matter halo mass. Finally, we compare our clustering+HOD estimates with results from abundance matching techniques, and conclude that the M_h estimates of the clustering+HOD analyses agree with those of the simple abundance matching within a factor of 3, and that the agreement improves when using more sophisticated abundance matching techniques that include subhalos, incompleteness, and/or evolution in the star formation and stellar mass functions.

We report the first identification of a compact dense galaxy overdensity at $z\sim8$ dubbed A2744z8OD. A2744z8OD consists of eight $Y$-dropout galaxies behind Abell 2744 that is originally pinpointed by Hubble Frontier Fields studies. However, no studies have, so far, derived basic physical quantities of structure formation or made comparisons with theoretical models. We obtain a homogeneous sample of dropout galaxies at $z\sim8$ from eight field data of Hubble legacy images that are as deep as the A2744z8OD data. Using the sample, we find that a galaxy surface overdensity value of A2744z8OD is very high $\delta\simeq130$, where $\delta$ is defined by an overdensity in a small circle of $6"$ ($\simeq30$ physical kpc) radius. Because there is no such a large $\delta$ value reported for high-$z$ overdensities to date, A2744z8OD is a system clearly different from those found in previous high-$z$ overdensity studies. In the galaxy+structure formation models of Henriques et al. (2015), there exist a very similar overdensity, Modelz8OD, that is made of eight model dropout galaxies at $z\sim8$ in a $6"$-radius circle. Modelz8OD is a progenitor of a today's $10^{14}M_\odot$ cluster, and more than a half of the seven Modelz8OD galaxies are merged into the brightest cluster galaxy (BCG) of the cluster. If Modelz8OD is a counterpart of A2744z8OD, the models suggest that A2744z8OD would be a part of a forming cluster core of a today's $10^{14}M_\odot$ cluster that started star formation at $z>12$.

The properties of the first generation of stars and their supernova (SN) explosions remains unknown due to the lack of their actual observations. Recently many transient surveys are conducted and the feasibility of the detection of supernovae (SNe) of Pop III stars is growing. In this paper we study the multicolor light curves for a number of metal-free core-collapse SN models (25-100 M$_{\odot}$) to provide the indicators for finding and identification of first generation SNe. We use mixing-fallback supernova explosion models which explain the observed abundance patterns of metal poor stars. Numerical calculations of the multicolor light curves are performed using multigroup radiation hydrodynamic code STELLA. The calculated light curves of metal-free SNe are compared with non-zero metallicity models and several observed SNe. We have found that the shock breakout characteristics, the evolution of the photosphere's velocity, the luminosity, the duration and color evolution of the plateau - all the SN phases are helpful to estimate the parameters of SN progenitor: the mass, the radius, the explosion energy and the metallicity. We conclude that the multicolor light curves can be potentially used to identify first generation SNe in the current (Subaru/HSC) and future transient surveys (LSST, JWST). They are also suitable for identification of the low-metallicity SNe in the nearby Universe (PTF, Pan-STARRS, Gaia).

Na I D absorbing systems toward Type Ia supernovae (SNe Ia) have been intensively studied over the last decade with the aim of finding circumstellar material (CSM), which is an indirect probe of the progenitor system. However, it is difficult to deconvolve CSM components from non-variable, and often dominant, components created by interstellar material (ISM). We present a series of high-resolution spectra of SN Ia 2014J from before maximum brightness to ~250 days after maximum brightness. The late-time spectrum provides unique information for determining the origin of the Na I D absorption systems. The deep late-time observation allows us to probe the environment around the SN at a large scale, extending to ~40 pc. We find that a spectrum of diffuse light in the vicinity, but not directly in the line-of-sight, of the SN has absorbing systems nearly identical to those obtained for the "pure" SN line-of-sight. Therefore, basically all Na I D systems seen toward SN 2014J must originate from foreground material that extends to at least ~40 pc in projection and none at the CSM scale. A fluctuation in the column densities at a scale of ~20 pc is also identified. After subtracting the diffuse, "background" spectrum, the late-time SN Na I D profile along the SN line-of-sight is consistent with the profile in the near-maximum brightness spectra. The lack of variability on a ~1 year timescale is consistent with the ISM interpretation for the gas.

We conduct precise strong lensing mass modeling of four ${\it Hubble}$ Frontier Fields (HFF) clusters, Abell$~$2744, MACS$~$J0416.1$-$2403, MACS$~$J0717.5$+$3745, and MACS$~$J1149.6$+$2223, for which HFF imaging observations are completed. We construct a refined sample of more than 100 multiple images for each cluster by taking advantage of the full depth HFF images, and conduct mass modeling using the ${\small \rm GLAFIC}$ software, which assumes simply parametrized mass distributions. Our mass modeling also exploits a magnification constraint from the lensed Type Ia supernova HFF14Tom for Abell$~$2744 and positional constraints from the multiple images S1$-$S4 of the lensed supernova SN Refsdal for MACS$~$J1149.6$+$2223. We find that our best-fitting mass models reproduce the observed image positions with RMS errors of $\sim 0.4$ arcsec, which are smaller than RMS errors in previous mass modeling that adopted similar numbers of multiple images. Our model predicts a new image of SN Refsdal with a relative time delay and magnification that are fully consistent with a recent detection of reappearance. We then construct catalogs of $z\sim 6-9$ dropout galaxies behind the four clusters and estimate magnification factors for these dropout galaxies with our best-fitting mass models. The dropout sample from the four cluster fields contains $\sim 120$ galaxies at $z\gtrsim 6$, about 20 of which are predicted to be magnified by a factor of more than 10. Some of the high-redshift galaxies detected in the HFF have lensing-corrected magnitudes of $M_{\rm UV}\sim -15$ to $-14$. Our analysis demonstrates that the HFF data indeed offer an ideal opportunity to study faint high-redshift galaxies. All lensing maps produced from our mass modeling will be made available on the STScI website.

Supernova "Refsdal," multiply imaged by cluster MACSJ1149.5+2223, represents a rare opportunity to make a true blind test of model predictions in extragalactic astronomy, on a time scale that is short compared to a human lifetime. In order to take advantage of this event, we produced seven gravitational lens models with five independent methods, based on Hubble Space Telescope (HST) Hubble Frontier Field images, along with extensive spectroscopic follow-up observations by HST, the Very Large and the Keck Telescopes. We compare the model predictions and show that they agree reasonably well with the measured time delays and magnification ratios between the known images, even though these quantities were not used as input. This agreement is encouraging, considering that the models only provide statistical uncertainties, and do not include additional sources of uncertainties such as structure along the line of sight, cosmology, and the mass sheet degeneracy. We then present the model predictions for the other appearances of SN "Refsdal." A future image will reach its peak in the first half of 2016, while another image appeared between 1994 and 2004. The past image would have been too faint to be detected in existing archival images. The future image should be approximately one third as bright as the brightest known image (i.e., H_AB~25.7 mag at peak and H_AB~26.7 mag six months before peak), and thus detectable in single-orbit HST images. We will find out soon whether our predictions are correct.

We present statistics of 133 faint 1.2-mm continuum sources detected in about 120 deep ALMA pointing data that include all the archival deep data available by 2015 June. We derive number counts of 1.2 mm continuum sources down to 0.02 mJy partly with the assistance of gravitational lensing, and find that the total integrated 1.2 mm flux of the securely identified sources is 22.9^(+6.7)_(-5.6) Jy deg^(-2) that corresponds to 104^(+31)_(-25) % of the extragalactic background light (EBL) measured by COBE observations. These results suggest that the major 1.2 mm EBL contributors are sources with 0.02 mJy, and that very faint 1.2 mm sources with ~< 0.02 mJy contribute negligibly to the EBL with the possible flattening and/or truncation of number counts in this very faint flux regime. To understand the physical origin of our faint ALMA sources, we measure the galaxy bias bg by the counts-in-cells technique, and place a stringent upper limit of bg < 3.5 that is not similar to bg values of massive DRGs and SMGs but comparable to those of UV-bright sBzKs and LBGs. Moreover, in optical and near-infrared (NIR) deep fields, we identify optical-NIR counterparts for 59% of our faint ALMA sources, majority of which have luminosities, colors, and the IRX-beta relation same as sBzKs and LBGs. We thus conclude that about a half of our faint ALMA sources are dust-poor high-z galaxies as known as sBzKs and LBGs in optical studies, and that these faint ALMA sources are not miniature (U)LIRGs simply scaled down with the infrared brightness.

We report the late-time evolution of Type IIb Supernova (SN IIb) 2013df. SN 2013df showed a dramatic change in its spectral features at ~1 year after the explosion. Early on it showed typical characteristics shared by SNe IIb/Ib/Ic dominated by metal emission lines, while later on it was dominated by broad and flat-topped Halpha and He I emissions. The late-time spectra are strikingly similar to SN IIb 1993J, which is the only previous example clearly showing the same transition. This late-time evolution is fully explained by a change in the energy input from the $^{56}$Co decay to the interaction between the SN ejecta and dense circumstellar matter (CSM). The mass loss rate is derived to be ~(5.4 +- 3.2) x 10^{-5} Msun/yr (for the wind velocity of ~20 km/s), similar to SN 1993J but larger than SN IIb 2011dh by an order of magnitude. The striking similarity between SNe 2013df and 1993J in the (candidate) progenitors and the CSM environments, and the contrast in these natures to SN 2011dh, infer that there is a link between the natures of the progenitor and the mass loss: SNe IIb with a more extended progenitor have experienced a much stronger mass loss in the final centuries toward the explosion. It might indicate that SNe IIb from a more extended progenitor are the explosions during a strong binary interaction phase, while those from a less extended progenitor have a delay between the strong binary interaction and the explosion.

Chemical abundance of r-process elements in nearby dwarf spheroidal (dSph) galaxies is a powerful tool to probe the site of r-process since their small-mass scale can sort out individual events producing r-process elements. A merger of binary neutron stars is a promising candidate of this site. In faint, or less massive dSph galaxies such as the Draco, a few binary neutron star mergers are expected to have occurred at most over the whole past. We have measured chemical abundances including Eu and Ba of three red giants in the Draco dSph by Subaru/HDS observation. The Eu detection for one star with [Fe/H]=-1.45 confirms a broadly constant [Eu/H] of ~-1.3 for stars with [Fe/H]>-2. This feature is shared by other dSphs with similar masses, i.e., the Sculptor and the Carina, and suggests that neutron star merger is the origin of r-process elements in terms of its rarity. In addition, two very metal-poor stars with [Fe/H]=-2.12 and -2.51 are found to exhibit very low Eu abundances such as [Eu/H]<-2 with an implication of a sudden increase of Eu abundance by more than 0.7 dex at [Fe/H] ~ -2.2 in the Draco dSph. In addition, the detection of Ba abundances for these stars suggests that the r-process enrichment initiated no later than the time when only a few % of stars in the present-day Draco dSph was formed. Though an identification of the origin of an early Eu production inside the Draco dSph should be awaited until more abundance data of stars with [Fe/H]<-2 in the Draco as well as other faint dSphs become available, the implied early emergence of Eu production event might be reconciled with the presence of extremely metal-poor stars enriched by r-process elements in the Galactic halo.

We investigate sizes of $z\sim6-8$ dropout galaxies using the complete data of the Abell 2744 cluster and parallel fields in the Hubble Frontier Fields program. By directly fitting light profiles of observed galaxies with lensing-distorted S\'ersic profiles on the image plane with the \texttt{glafic} software, we accurately measure intrinsic sizes of 31 $z\sim6-7$ and eight $z\sim8$ galaxies, including those as faint as $M_{\mathrm{UV}}\simeq-16.6$. We find that half-light radii $r_\mathrm{e}$ positively correlates with UV luminosity at each redshift, although the correlation is not very tight. Largest ($r_\mathrm{e}>0.8$ kpc) galaxies are mostly red in UV color while smallest ($r_\mathrm{e} < 0.08$ kpc) ones tend to be blue. We also find that galaxies with multiple cores tend to be brighter. Combined with previous results at $2.5\lesssim z\lesssim12$, our result confirms that the average $r_{\mathrm{e}}$ of bright ($(0.3-1)L^\star_{z=3}$) galaxies scales as $r_{\mathrm{e}}\propto(1+z)^{-m}$ with $m=1.24\pm0.1$. We find that the ratio of $r_\mathrm{e}$ to virial radius is virtually constant at $3.3\pm0.1\%$ over a wide redshift range, where the virial radii of hosting dark matter halos are derived based on the abundance matching. This constant ratio is consistent with the disk formation model by Mo et al. (1998) with $j_\mathrm{d}\sim m_\mathrm{d}$, where $j_\mathrm{d}$ and $m_\mathrm{d}$ are the fractions of the angular momentum and mass within halos confined in the disks. A comparison with various types of local galaxies indicates that our galaxies are most similar to circumnuclear star-forming regions of barred galaxies in the sense that a sizable amount of stars are forming in a very small area.

We present the comprehensive analyses of faint dropout galaxies up to $z\sim10$ with the first full-depth data set of Abell 2744 lensing cluster and parallel fields observed by the Hubble Frontier Fields (HFF) program. We identify $54$ dropouts at $z\sim5-10$ in the HFF fields, and enlarge the size of $z\sim9$ galaxy sample obtained to date. Although the number of highly magnified ($\mu\sim10$) galaxies is small due to the tiny survey volume of strong lensing, our study reaches the galaxies' intrinsic luminosities comparable to the deepest-field HUDF studies. We derive UV luminosity functions with these faint dropouts, carefully evaluating the combination of observational incompleteness and lensing effects in the image plane by intensive simulations including magnification, distortion, and multiplication of images, with the evaluations of mass model dependences. Our results confirm that the faint-end slope, $\alpha$, is as steep as $-2$ at $z\sim6-8$, and strengthen the evidence of the rapid decrease of UV luminosity densities, $\rho_\mathrm{UV}$, at $z>8$ from the large $z\sim9$ sample. We examine whether the rapid $\rho_\mathrm{UV}$ decrease trend can reconcile with the large Thomson scattering optical depth, $\tau_\mathrm{e}$, measured by CMB experiments allowing a large space of free parameters such as average ionizing photon escape fraction and stellar-population dependent conversion factor. No parameter set can reproduce both the rapid $\rho_\mathrm{UV}$ decrease and the large $\tau_\mathrm{e}$. It is possible that the $\rho_\mathrm{UV}$ decrease moderates at $z\gtrsim11$, that the free parameters significantly evolve towards high-$z$, or that there exist additional sources of reionization such as X-ray binaries and faint AGNs.

The most iron-deficient stars in the Milky Way provide important observational constraints on the nature of astrophysical objects that have enriched the primordial gas with heavy elements from which these stars were formed. Among them, the recently discovered iron-deficient star SMSS J031300.36-670839.3 shows a remarkable chemical composition with non-detection of iron ([Fe/H]$<-7.1$) and large enhancement of carbon and magnesium relative to calcium. We investigate the supernova yields of metal-free (Population III) stars to interpret the observed abundance pattern for this star. We report that the high [C/Ca] and [C/Mg] ratios and upper limits determined for other elemental abundances are well reproduced with the yields of core-collapse supernovae (that have normal kinetic energies of explosion $E$ of $E_{51}=E/10^{51}$erg$=1$) or hypernovae ($E_{51}\geq 10$) of the Population III 25$M_{\odot}$ or 40$M_{\odot}$ stars. The best-fit model assumes that the explosion of the Population III progenitor undergoes extensive matter mixing and fallback, leaving behind a black hole remnant. In these models, Ca is produced by static/explosive O burning and incomplete Si burning in the Population III supernova/hypernova, in contrast to the suggestion that Ca is originated from the hot CNO cycle during the presupernova evolution. Among the five most iron-poor stars with [Fe/H]$<-4.5$, four carbon-enhanced stars are consistent with the faint supernova models with the ejected mass of $^{56}$Ni less than 10$^{-3}M_{\odot}$.

We show that the low ratios of $\alpha$ elements (Mg, Si, and Ca) to Fe recently found for a small fraction of extremely metal-poor stars can be naturally explained with the nucleosynthesis yields of core-collapse supernovae, i.e., $13-25M_\odot$ supernovae, or hypernovae. For the case without carbon enhancement, the ejected iron mass is normal, consistent with observed light curves and spectra of nearby supernovae. On the other hand, the carbon enhancement requires much smaller iron production, and the low [$\alpha$/Fe] of carbon enhanced metal-poor stars can also be reproduced with $13-25M_\odot$ faint supernovae or faint hypernovae. Iron-peak element abundances, in particular Zn abundances, are important to put further constraints on the enrichment sources from galactic archaeology surveys.

Ultra-faint dwarf galaxies recently discovered around the Milky Way (MW) contain extremely metal-poor stars, and might represent the building blocks of low-metallicity components of the MW. Among them, the Bo\"otes I dwarf spheroidal galaxy is of particular interest because of its exclusively old stellar population. We determine chemical compositions of six red giant stars in Bo\"otes I, based on the high-resolution spectra obtained with the High Dispersion Spectrograph mounted on the Subaru Telescope. Abundances of 12 elements, including C, Na, alpha, Fe-peak, and neutron capture elements, were determined for the sample stars. The abundance results were compared to those in field MW halo stars previously obtained using an abundance analysis technique similar to the present study. We confirm the low metallicity of Boo-094 ([Fe/H]=-3.4). Except for this star, the abundance ratios ([X/Fe]) of elements lighter than zinc are generally homogeneous with small scatter around the mean values in the metallicities spanned by the other five stars (-2.7<[Fe/H]<-1.8). Specifically, all of the sample stars with [Fe/H]>-2.7 show no significant enhancement of carbon. The [Mg/Fe] and [Ca/Fe] ratios are almost constant with a modest decreasing trend with increasing [Fe/H] and are slightly lower than the field halo stars. The [Sr/Fe] and [Sr/Ba] ratios also tend to be lower in the Bo\"otes I stars than in the halo stars. Our results of small scatter in the [X/Fe] ratios for elements lighter than zinc suggest that these abundances were homogeneous among the ejecta of prior generation(s) of stars in this galaxy.

We present chemical abundance analyses of sodium, iron-peak and neutron-capture elements for 97 kinematically selected thick disk, inner halo and outer halo stars with metallicities -3.3<[Fe/H]<-0.5. The main aim of this study is to examine chemical similarities and differences among metal-poor stars belonging to these old Galactic components as a clue to determine their early chemodynamical evolution. In our previous paper, we obtained abundances of alpha elements by performing a one-dimensional LTE abundance analysis based on the high-resolution (R~50000) spectra obtained with the Subaru/HDS. In this paper, a similar analysis is performed to determine abundances of an additional 17 elements. We show that, in metallicities below [Fe/H]~-2, the abundance ratios of many elements in the thick disk, inner halo, and outer halo subsamples are largely similar. In contrast, in higher metallicities ([Fe/H]>-1.5), differences in some of the abundance ratios among the three subsamples are identified. Specifically, the [Na/Fe], [Ni/Fe], [Cu/Fe], and [Zn/Fe] ratios in the inner and outer halo subsamples are found to be lower than those in the thick disk subsample. In contrast to what was observed for [Mg/Fe] in our previous paper, [Eu/Fe] ratios are more enhanced in the two halo subsamples rather than in the thick disk subsample. The observed distinct chemical abundances of some elements between the thick disk and inner/outer halo subsamples with [Fe/H]>-1.5 support the hypothesis that these components formed through different mechanisms. In particular, our results favor the scenario that the inner and outer halo components formed through an assembly of multiple progenitor systems that experienced various degrees of chemical enrichments, while the thick disk formed through rapid star formation with an efficient mixing of chemical elements.

We present an analysis of a high resolution (R~50,000) optical spectrum of the central region of the proto-planetary nebula CRL 2688. This object is thought to have recently moved off the AGB, and display abundance patterns of CNO and heavy elements that can provide us with important clues to understand the nucleosynthesis, dredge-up and mixing experienced by the envelope of the central star during its AGB stage of evolution. The analysis of the molecular features, presumably originated from the circumstellar matter provides further constraints on the chemistry and velocity of the expanding shell, expelled as a consequence of the strong mass loss experienced by the central star. We confirm that the central star shows a spectrum typical of an F-type supergiant with Teff=7250 K, log g=0.5 and [Fe/H]=-0.3 dex. We find that the abundance pattern of this object is characterized by enhancements of Carbon ([C/Fe]=0.6), Nitrogen ([N/Fe]=1.0) and Na ([Na/Fe]=0.7), similar to other previously known carbon-rich post-AGB stars. Yttrium is also enhanced while the [Ba/Y] ratio is very low (-1.0), indicating that only the light s-process elements are enhanced. The Zinc abundance is found to be normal, [Zn/Fe]=0.0, suggesting that there is no depletion of refractory elements. The Halpha, Na I and K I resonance lines show prominent emission components, whose helio-centric radial velocities are offsetted by -41 km/s relative to the photospheric metal-absorption lines. The molecular features of C_2 and CN also show emission components, whose velocities are consistent with the emission components of the Halpha, Na I, and K I lines. On the other hand, their absorption components are more highly blue shifted than the corresponding emission components, which suggests that the regions where the emission and absorption components arise are expanding at different velocities.

We present the abundance analysis of 97 nearby metal-poor (-3.3<[Fe/H]<-0.5) stars having kinematics characteristics of the Milky Way (MW) thick disk, inner, and outer stellar halos. The high-resolution, high-signal-to-noise optical spectra for the sample stars have been obtained with the High Dispersion Spectrograph mounted on the Subaru Telescope. Abundances of Fe, Mg, Si, Ca and Ti have been derived using a one-dimensional LTE abundance analysis code with Kurucz NEWODF model atmospheres. By assigning membership of the sample stars to the thick disk, inner or outer halo components based on their orbital parameters, we examine abundance ratios as a function of [Fe/H] and kinematics for the three subsamples in wide metallicity and orbital parameter ranges. We show that, in the metallicity range of -1.5<[Fe/H]<= -0.5, the thick disk stars show constantly high mean [Mg/Fe] and [Si/Fe] ratios with small scatter. In contrast, the inner, and the outer halo stars show lower mean values of these abundance ratios with larger scatter. The [Mg/Fe], [Si/Fe] and [Ca/Fe] for the inner and the outer halo stars also show weak decreasing trends with [Fe/H] in the range [Fe/H]$>-2$. These results favor the scenarios that the MW thick disk formed through rapid chemical enrichment primarily through Type II supernovae of massive stars, while the stellar halo has formed at least in part via accretion of progenitor stellar systems having been chemically enriched with different timescales.

We present chemical abundances of 57 metal-poor stars that are likely constituents of the outer stellar halo in the Milky Way. Almost all of the sample stars have an orbit reaching a maximum vertical distance (Z_max) of >5 kpc above and below the Galactic plane. High-resolution, high signal-to-noise spectra for the sample stars obtained with Subaru/HDS are used to derive chemical abundances of Na, Mg, Ca, Ti, Cr, Mn, Fe, Ni, Zn, Y and Ba with an LTE abundance analysis code. The resulting abundance data are combined with those presented in literature that mostly targeted at smaller Z_max stars, and both data are used to investigate any systematic trends in detailed abundance patterns depending on their kinematics. It is shown that, in the metallicity range of -2<[Fe/H]<-1, the [Mg/Fe] ratios for the stars with Z_max>5 kpc are systematically lower (~0.1 dex) than those with smaller Z_max. This result of the lower [alpha/Fe] for the assumed outer halo stars is consistent with previous studies that found a signature of lower [alpha/Fe] ratios for stars with extreme kinematics. A distribution of the [Mg/Fe] ratios for the outer halo stars partly overlaps with that for stars belonging to the Milky Way dwarf satellites in the metallicity interval of -2<[Fe/H]<-1 and spans a range intermediate between the distributions for the inner halo stars and the stars belonging to the satellites. Our results confirm inhomogeneous nature of chemical abundances within the Milky Way stellar halo depending on kinematic properties of constituent stars as suggested by earlier studies. Possible implications for the formation of the Milky Way halo and its relevance to the suggested dual nature of the halo are discussed.

Chemical abundances of 26 metal-poor dwarfs and giants are determined from high-resolution and high signal-to-noise ratio spectra obtained with Subaru/HDS. The sample is selected so that most of the objects have outer-halo kinematics. Self-consistent atmospheric parameters were determined by an iterative procedure based on spectroscopic analysis. Abundances of 13 elements, including $\alpha$-elements (Mg, Si, Ca, Ti), odd-Z light elements (Na, Sc), iron-peak elements (Cr, Mn, Fe, Ni, Zn) and neutron-capture elements (Y, Ba), are determined by two independent data reduction and LTE analysis procedures, confirming the consistency of the stellar parameters and abundances results. We find a decreasing trend of [$\alpha$/Fe] with increasing [Fe/H] for the range of $-3.5 <$ [Fe/H]$ < -1$, as found by Stephens and Boesgaard (2002). [Zn/Fe] values of most objects in our sample are slightly lower than the bulk of halo stars previously studied. These results are discussed as possible chemical properties of the outer halo in the Galaxy.

Passive spiral galaxies, despite their spiral morphological appearance, do not have any emission lines indicative of ongoing star formation in their optical spectra. Previous studies have suggested that passive spiral galaxies preferentially exist in infall regions of galaxy clusters, suggesting that the cluster environment is likely to be responsible for creating these galaxies. By carrying out spatially resolved long-slit spectroscopy on four nearby passive spiral galaxies with the Apache Point Observatory 3.5-m telescope, we investigated the stellar populations of passive spiral galaxies separately for their inner and outer regions. In the two unambiguously passive spiral galaxies among the four observed galaxies, H$\delta$ absorption lines are more prominent in the outer regions of the galaxies, whereas the 4000-{\AA} breaks (D$_{4000}$) are strongest in the inner regions of the galaxies. A comparison with a simple stellar population model for the two passive spiral galaxies indicates that the outer regions of the galaxies tend to harbour younger populations of stars. The strong H$\delta$ absorption observed in the outer regions of the sample galaxies is consistent with that of galaxies whose star formation ceased a few Gyrs ago. Because of the large uncertainty in the absorption indices in our samples, further observations are needed in order to place constraints on the mechanisms that quench star formation in passive spiral galaxies.

The Infrared Camera (IRC) is one of the two instruments on board the AKARI satellite. In addition to deep imaging from 1.8-26.5um for the pointed observation mode of the AKARI, it has a spectroscopic capability in its spectral range. By replacing the imaging filters by transmission-type dispersers on the filter wheels, it provides low-resolution (lambda/d_lambda ~ 20-120) spectroscopy with slits or in a wide imaging field-of-view (approximately 10'X10'). The IRC spectroscopic mode is unique in space infrared missions in that it has the capability to perform sensitive wide-field spectroscopic surveys in the near- and mid-infrared wavelength ranges. This paper describes specifications of the IRC spectrograph and its in-orbit performance.

AKARI, the first Japanese satellite dedicated to infrared astronomy, was launched on 2006 February 21, and started observations in May of the same year. AKARI has a 68.5 cm cooled telescope, together with two focal-plane instruments, which survey the sky in six wavelength bands from the mid- to far-infrared. The instruments also have the capability for imaging and spectroscopy in the wavelength range 2 - 180 micron in the pointed observation mode, occasionally inserted into the continuous survey operation. The in-orbit cryogen lifetime is expected to be one and a half years. The All-Sky Survey will cover more than 90 percent of the whole sky with higher spatial resolution and wider wavelength coverage than that of the previous IRAS all-sky survey. Point source catalogues of the All-Sky Survey will be released to the astronomical community. The pointed observations will be used for deep surveys of selected sky areas and systematic observations of important astronomical targets. These will become an additional future heritage of this mission.

The Infrared Camera (IRC) is one of two focal-plane instruments on the AKARI satellite. It is designed for wide-field deep imaging and low-resolution spectroscopy in the near- to mid-infrared (1.8--26.5um) in the pointed observation mode of AKARI. IRC is also operated in the survey mode to make an all-sky survey at 9 and 18um. It comprises three channels. The NIR channel (1.8--5.5um) employs a 512 x 412 InSb array, whereas both the MIR-S (4.6--13.4um) and MIR-L (12.6--26.5um) channels use 256 x 256 Si:As impurity band conduction arrays. Each of the three channels has a field-of-view of about 10' x 10' and are operated simultaneously. The NIR and MIR-S share the same field-of-view by virtue of a beam splitter. The MIR-L observes the sky about $25' away from the NIR/MIR-S field-of-view. IRC gives us deep insights into the formation and evolution of galaxies, the evolution of planetary disks, the process of star-formation, the properties of interstellar matter under various physical conditions, and the nature and evolution of solar system objects. The in-flight performance of IRC has been confirmed to be in agreement with the pre-flight expectation. This paper summarizes the design and the in-flight operation and imaging performance of IRC.