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.

We investigate the orbital motions of dwarf spheroidal galaxies (dSphs) in the halo of the Milky Way (MW) to understand their possible effects on the diversity of the star formation histories seen in these MW satellites. In this work, we explicitly consider a time-varying gravitational potential due to the growth of the MW's dark halo mass to calculate the long-term orbital evolutions of the dSphs, guided with {\it Gaia} DR2 proper motions, over the past 13.5 billion years. We find that the infall time of a satellite, defined at which the galaxy first crosses within the growing virial radius of the MW's halo, coincides well with the time when the star formation rate (SFR) is peaked for the sample of classical dSphs. On the other hand, ultra-faint dSphs already finished their SF activity prior to their infall times as already suggested in previous works, but there is a signature that their earlier SF histories are affected by interaction with the growing MW's halo to some extent. We also find, for classical dSphs, that the relative fraction of stars formed after the peak of the SFR to the current stellar mass is smaller for the smaller pericentric radius of the galaxy at its first infall. These results suggest that the infalling properties of the dSphs into the MW and the resultant environmental effects such as ram-pressure stripping and/or tidal disturbance in the MW's dark halo containing hot gas play important roles in their star formation histories.

As deformations of the main reflector of a radio telescope directly affect the observations, the evaluation of the deformation is extremely important. Dynamic characteristics of the main reflector of the Nobeyama 45 m radio telescope, Japan, are measured under two conditions: The first is when the pointing observation is in operation, and the second is when the reflector is stationary and is subjected to wind loads when the observation is out of operation. Dynamic characteristics of the main reflector are measured using piezoelectric accelerometers. When the telescope is in operation, a vibration mode with one nodal line horizontally or vertically on the reflector is induced, depending on whether the reflector is moving in the azimuthal or elevational planes, whereas under windy conditions, vibration modes that have two to four nodal lines are simultaneously induced. The predominant mode is dependent on the direction of wind loads.

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.

We investigate the basic properties of voids from high resolution, cosmological N-body simulations of {\Lambda}-dominated cold dark matter ({\Lambda}CDM) models, in order to compare with the analytical model of Sheth and van de Weygaert (SvdW) for void statistics. For the subsample of five dark matter simulations in the {\Lambda}CDM cosmology with box sizes ranging from 1000 Mpc/h to 8 Mpc/h, we find that the standard void-in-cloud effect is too simplified to explain several properties of identified small voids in simulations. (i) The number density of voids is found to be larger than the prediction of the analytical model up to 2 orders of magnitude below 1 Mpc/h scales. The Press-Schechter model with the linear critical threshold of void {\delta}_v = -2.71, or a naive power law, is found to provide an excellent agreement with the void size function, suggesting that the void-in-cloud effect does not suppress as much voids as predicted by the SvdW model. (ii) We then measured the density and velocity profiles of small voids, and find that they are mostly partially collapsing underdensities, instead of being completely crushed in the standard void-in-cloud scenario. (iii) Finally, we measure the void distributions in four different tidal environments, and find that the void-in-void effect alone can explain the correlation between distribution and environments, whereas the void-in-cloud effect is only weakly influencing the abundance of voids, even in filaments and clusters.

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]

We report on the global structure of the Milky Way (MW) stellar halo up to its outer boundary based on the analysis of blue-horizontal branch stars (BHBs). These halo tracers are extracted from the $(g,r,i,z)$ band multi-photometry in the internal data release of the on-going Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) surveyed over $\sim550$~deg$^2$ area. In order to select most likely BHBs by removing blue straggler stars (BSs) and other contamination in a statistically significant manner, we have developed and applied an extensive Bayesian method, instead of the simple color cuts adopted in our previous work, where each of the template BHBs and non-BHBs obtained from the available catalogs is represented as a mixture of multiple Gaussian distributions in the color-color diagrams. We found from the candidate BHBs in the range of 18.5<g<23.5 mag that the radial density distribution over a Galactocentric radius of r=36-360 kpc can be approximated as a single power-law profile with an index of $\alpha=3.74^{+0.21}_{-0.22}$ or a broken power-law profile with an index of $\alpha_{\rm in}=2.92^{+0.33}_{-0.33}$ at $r$ below a broken radius of $r_{\rm b}=160^{+18}_{-19}$~kpc and a very steep slope of $\alpha_{\rm out}=15.0^{+3.7}_{-4.5}$ at $r>r_{\rm b}$. The latter profile with a prolate shape having an axial ratio of $q=1.72^{+0.44}_{-0.28}$ is most likely and this halo may hold a rather sharp boundary at r=160kpc. The slopes of the halo density profiles are compared with those from the suite of hydrodynamical simulations for the formation of stellar halos. This comparison suggests that the MW stellar halo may consist of the two overlapping components: the in situ. inner halo as probed by RR Lyrae stars showing a relatively steep radial density profile and the ex situ. outer halo with a shallow profile probed by BHBs here, which is made by accretion of small stellar systems.

Primordial black holes (PBHs) have long been suggested as a viable candidate for the elusive dark matter (DM). The abundance of such PBHs has been constrained using a number of astrophysical observations, except for a hitherto unexplored mass window of $M_{\rm PBH}=[10^{-14},10^{-9}]M_\odot$. Here we carry out a dense-cadence (2~min sampling rate), 7 hour-long observation of the Andromeda galaxy (M31) with the Subaru Hyper Suprime-Cam to search for microlensing of stars in M31 by PBHs lying in the halo regions of the Milky Way (MW) and M31. Given our simultaneous monitoring of tens of millions of stars in M31, if such light PBHs make up a significant fraction of DM, we expect to find many microlensing events for the PBH DM scenario. However, we identify only a single candidate event, which translates into the most stringent upper bounds on the abundance of PBHs in the mass range $M_{\rm PBH}\simeq [10^{-11}, 10^{-6}]M_\odot$.

We present results from a pilot observation of nearby (~20 Mpc) galaxies with mass similar to that of the Milky Way (MW) to address the missing satellite problem. This is the first paper from an on-going project to address the problem with a statistical sample of galaxies outside of the Local Group (LG) without employing an assumption that the LG is a typical halo in the Universe. Thanks to the close distances of our targets, dwarf galaxies around them can be identified as extended, diffuse galaxies. By applying a surface brightness cut together with a careful visual screening to remove artifacts and background contamination, we construct a sample of dwarf galaxies. The luminosity function (LF) of one of the targets is broadly consistent with that of the MW, but the other has a more abundant dwarf population. Numerical simulations by Okamoto (2013) seem to overpredict the number of dwarfs on average, while more recent predictions from Copernicus Complexio are in a better agreement. In both observations and simulations, there is a large diversity in the LFs, demonstrating the importance of addressing the missing satellite problem with a statistically representative sample. We also characterize the projected spatial distributions of the satellites and do not observe strong evidence for alignments around the central galaxies. Based on this successful pilot observation, we are carrying out further observations to increase the sample of nearby galaxies, which we plan to report in our future paper.

We present the structure of the Milky Way stellar halo beyond Galactocentric distances of $r = 50$ kpc traced by blue horizontal-branch (BHB) stars, which are extracted from the survey data in the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). We select BHB candidates based on $(g,r,i,z)$ photometry, where the $z$-band is on the Paschen series and the colors that involve the $z$-band are sensitive to surface gravity. About 450 BHB candidates are identified between $r = 50$ kpc and 300 kpc, most of which are beyond the reach of previous large surveys including the Sloan Digital Sky Survey. We find that the global structure of the stellar halo in this range has substructures, which are especially remarkable in the GAMA15H and XMM-LSS fields in the HSC-SSP. We find that the stellar halo can be fitted to a single power-law density profile with an index of $\alpha \simeq 3.3$ ($3.5$) with (without) these fields and its global axial ratio is $q \simeq 2.2$ ($1.3$). Thus, the stellar halo may be significantly disturbed and be made in a prolate form by halo substructures, perhaps associated with the Sagittarius stream in its extension beyond $r \sim 100$ kpc. For a broken power-law model allowing different power-law indices inside/outside a break radius, we obtain a steep power-law slope of $\alpha \sim 5$ outside a break radius of $\simeq 100$ kpc ($200$ kpc) for the case with (without) GAMA15H and XMM-LSS. This radius of $200$ kpc might be as close as a halo boundary if there is any, although larger BHB sample is required from further HSC-SSP survey to increase its statistical significance.

Hyper Suprime-Cam (HSC) is a wide-field imaging camera on the prime focus of the 8.2m Subaru telescope on the summit of Maunakea in Hawaii. A team of scientists from Japan, Taiwan and Princeton University is using HSC to carry out a 300-night multi-band imaging survey of the high-latitude sky. The survey includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad bands ($grizy$), with a $5\,\sigma$ point-source depth of $r \approx 26$. The Deep layer covers a total of 26~deg$^2$ in four fields, going roughly a magnitude fainter, while the UltraDeep layer goes almost a magnitude fainter still in two pointings of HSC (a total of 3.5 deg$^2$). Here we describe the instrument, the science goals of the survey, and the survey strategy and data processing. This paper serves as an introduction to a special issue of the Publications of the Astronomical Society of Japan, which includes a large number of technical and scientific papers describing results from the early phases of this survey.

We investigate the formation history of the stellar disk component in the Milky Way (MW) based on our new chemical evolution model. Our model considers several fundamental baryonic processes, including gas infall, re-accretion of outflowing gas, and radial migration of disk stars. Each of these baryonic processes in the disk evolution is characterized by model parameters, which are determined by fitting to various observational data of the stellar disk in the MW, including the radial dependence of the metallicity distribution function (MDF) of the disk stars, which has recently been derived in the APOGEE survey. We succeeded to obtain the best set of model parameters, which well reproduces the observed radial dependences of the mean, standard deviation, skewness, and kurtosis of the MDFs for the disk stars. We analyze the basic properties of our model results in detail to get new insights into the important baryonic processes in the formation history of the MW. One of the remarkable findings is that outflowing gas, containing much heavy elements, preferentially re-accretes onto the outer disk parts, and this recycling process of metal-enriched gas is a key ingredient to reproduce the observed narrower MDFs at larger radii. Moreover, important implications for the radial dependence of gas infall and the influence of radial migration on the MDFs are also inferred from our model calculation. Thus, the MDF of disk stars is a useful clue for studying the formation history of the MW.

We present wide and deep photometry of the northwest part of the halo of the Andromeda galaxy (M31) using Hyper Suprime-Cam on the Subaru Telescope. The survey covers 9.2 deg$^{2}$ field in the $g$, $i$, and $NB515$ bands and shows a clear red giant branch (RGB) of M31's halo stars and a pronounced red clump (RC) feature. The spatial distribution of RC stars shows a prominent stream feature, the North Western (NW) Stream, and a diffuse substructure in the south part of our survey field. We estimate the distances based on the RC method and obtain $(m-M)$ = 24.63$\pm 0.191$(random)$\pm0.057$(systematic) and 24.29$\pm 0.211$(random)$\pm0.057$(systematic) mag for the NW stream and diffuse substructure, respectively, implying that the NW Stream is located behind M31, whereas the diffuse substructure is located in front. We also estimate line-of-sight distances along the NW Stream and find that the south part of the stream is $\sim$20 kpc closer to us relative to the north part. The distance to the NW Stream inferred from the isochrone fitting to the color-magnitude diagram favors the RC-based distance, but the TRGB-based distance estimated for $NB515$-selected RGB stars does not agree with it. The surface number density distribution of RC stars across the NW Stream is found to be approximately Gaussian with a FWHM of $\sim$25 arcmin (5.7 kpc), with a slight skew to the south-west side. That along the NW Stream shows a complicated structure including variations in number density and a significant gap in the stream.

We present the velocity dispersion of red giant branch (RGB) stars in M31's halo, derived by modeling the line of sight velocity distribution of over 5000 stars in 50 fields spread throughout M31's stellar halo. The dataset was obtained as part of the SPLASH (Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo) Survey, and covers projected radii of 9 to 175 kpc from M31's center. All major structural components along the line of sight in both the Milky Way (MW) and M31 are incorporated in a Gaussian Mixture Model, including all previously identified M31 tidal debris features in the observed fields. The probability an individual star is a constituent of M31 or the MW, based on a set of empirical photometric and spectroscopic diagnostics, is included as a prior probability in the mixture model. The velocity dispersion of stars in M31's halo is found to decrease only mildly with projected radius, from 108 km/s in the innermost radial bin (8.2 to 14.1 kpc) to $\sim 80$ to 90 km/s at projected radii of $\sim 40$ to 130 kpc, and can be parameterized with a power-law of slope $-0.12\pm 0.05$. The quoted uncertainty on the power-law slope reflects only the precision of the method, although other sources of uncertainty we consider contribute negligibly to the overall error budget.

The Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) is a three-layered imaging survey aimed at addressing some of the most outstanding questions in astronomy today, including the nature of dark matter and dark energy. The survey has been awarded 300 nights of observing time at the Subaru Telescope and it started in March 2014. This paper presents the first public data release of HSC-SSP. This release includes data taken in the first 1.7 years of observations (61.5 nights) and each of the Wide, Deep, and UltraDeep layers covers about 108, 26, and 4 square degrees down to depths of i~26.4, ~26.5, and ~27.0 mag, respectively (5sigma for point sources). All the layers are observed in five broad bands (grizy), and the Deep and UltraDeep layers are observed in narrow bands as well. We achieve an impressive image quality of 0.6 arcsec in the i-band in the Wide layer. We show that we achieve 1-2 per cent PSF photometry (rms) both internally and externally (against Pan-STARRS1), and ~10 mas and 40 mas internal and external astrometric accuracy, respectively. Both the calibrated images and catalogs are made available to the community through dedicated user interfaces and database servers. In addition to the pipeline products, we also provide value-added products such as photometric redshifts and a collection of public spectroscopic redshifts. Detailed descriptions of all the data can be found online. The data release website is https://hsc-release.mtk.nao.ac.jp/.

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.

We present the first results of the Subaru/Hyper Suprime-Cam (HSC) survey of the interacting galaxy system, NGC4631 and NGC4656. From the maps of resolved stellar populations, we identify 11 dwarf galaxies (including already-known dwarfs) in the outer region of NGC4631 and the two tidal stellar streams around NGC4631, named Stream SE and Stream NW, respectively. This paper describes the fundamental properties of these tidal streams. Based on the tip of red giant branch method and the Bayesian statistics, we find that StreamSE (7.10 Mpc in Expected a posteriori, EAP, with the 90% credible intervals of [6.22, 7.29] Mpc) and StreamNW (7.91 Mpc in EAP with the 90% credible intervals of [6.44, 7.97] Mpc) are located in front of and behind NGC4631, respectively. We also calculate the metallicity distribution of stellar streams by comparing the member stars with theoretical isochrones on the color-magnitude diagram. We find that both streams have the same stellar population based on the Bayesian model selection method, suggesting that they originated from a tidal interaction between NGC4631 and a single dwarf satellite. The expected progenitor has a positively skewed metallicity distribution function with [M/H]_EAP=-0.92 with the 90% credible intervals of [-1.46, -0.51]. The stellar mass of the progenitor is estimated as 3.7 x 10e+8 Msun with the 90% credible intervals of [5.8 x 10e+6, 8.6 x 10e+9] Msun based on the mass-metallicity relation for Local group dwarf galaxies. This is in good agreement with an initial stellar mass of the progenitor presumed in the previous N-body simulation.

Motivated by a recently found interesting property of the dark halo surface density within a radius, $r_{\rm max}$, giving the maximum circular velocity, $V_{\rm max}$, we investigate it for dark halos of the Milky Way's and Andromeda's dwarf satellites based on cosmological simulations. We select and analyze the simulated subhalos associated with Milky Way-sized dark halos and find that the values of their surface densities, $\Sigma_{V_{\rm max}}$, are in good agreement with those for the observed dwarf spheroidal satellites even without employing any fitting procedures. This implies that this surface density would not be largely affected by any baryonic feedbacks and thus universal. Moreover, all subhalos on the small scales of dwarf satellites are expected to obey the relation $\Sigma_{V_{\rm max}}\propto V_{\rm max}$, irrespective of differences in their orbital evolutions, host halo properties, and observed redshifts. Therefore, we find that the universal scaling relation for dark halos on dwarf galaxy mass scales surely exists and provides us important clues to understanding fundamental properties of dark halos. We also investigate orbital and dynamical evolutions of subhalos to understand the origin of this universal dark halo relation and find that most of subhalos evolve generally along the $r_{\rm max}\propto V_{\rm max}$ sequence, even though these subhalos have undergone different histories of mass assembly and tidal stripping. This sequence, therefore, should be the key feature to understand the nature of the universality of $\Sigma_{V_{\rm max}}$.

We report the $4 \, \sigma$ detection of a faint object with a flux of ~ 0.3 mJy, in the vicinity of the quadruply lensed QSO MG0414+0534 using the Atacama Large Millimeter/submillimeter array (ALMA) Band 7. The object is most probably a dusty dark dwarf galaxy, which has not been detected in either the optical, near-infrared (NIR) or radio (cm) bands. An anomaly in the flux ratio of the lensed images observed in Band 7 and the mid-infrared (MIR) band and the reddening of the QSO light color can be simultaneously explained if we consider the object as a lensing substructure with an ellipticity ~ 0.7 at a redshift of $0.5 \lesssim z \lesssim 1$. Using the best-fit lens models with three lenses, we find that the dark matter plus baryon mass associated with the object is $\sim 10^9\, M_{\odot}$, the dust mass is $\sim 10^7\,M_{\odot}$ and the linear size is $\gtrsim 5\,$kpc. Thus our findings suggest that the object is a dusty dark dwarf galaxy. A substantial portion of faint submillimeter galaxies (SMGs) in the universe may be attributed to such dark objects.

We investigate the role of radial migration history of stars in chemical evolution of a disk galaxy, in particular in understanding the origin of their bimodal distribution on the [alpha/Fe]-[Fe/H] plane. For this purpose, we examine the three different models with no, continuous, and discontinuous radial migration, respectively. We find that for the model with radial migration, the [alpha/Fe] ratios of stars in outer disk regions decrease more rapidly with time than the model without radial migration, because the associated net transfer of intermediate and old disk stars from inner to outer disk regions increases the rate of Type Ia relative to that of Type II supernovae in the latter regions. Moreover, in the model assuming rapid and discontinuous radial migration, its effect on the stellar abundances at larger radii is significant enough to provide the large difference in the evolution of stars on the [alpha/Fe]-[Fe/H] plane between inner and outer disk regions. As a result we obtain the bimodal distribution of disk stars on the [alpha/Fe]-[Fe/H] plane as observed in the Galactic stellar disk, thereby implying that the event of discontinuous radial migration may play a key role in reproducing the observed bimodality of stars on this abundance-ratio diagram. We discuss possible mechanisms causing such discontinuous radial migration in the early evolution of the Galactic disk, including the event of minor merging of a relatively massive satellite onto the stellar disk.

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.

PFS (Prime Focus Spectrograph), a next generation facility instrument on the 8.2-meter Subaru Telescope, is a very wide-field, massively multiplexed, optical and near-infrared spectrograph. Exploiting the Subaru prime focus, 2394 reconfigurable fibers will be distributed over the 1.3 deg field of view. The spectrograph has been designed with 3 arms of blue, red, and near-infrared cameras to simultaneously observe spectra from 380nm to 1260nm in one exposure at a resolution of ~1.6-2.7A. An international collaboration is developing this instrument under the initiative of Kavli IPMU. The project is now going into the construction phase aiming at undertaking system integration in 2017-2018 and subsequently carrying out engineering operations in 2018-2019. This article gives an overview of the instrument, current project status and future paths forward.

We investigate the growth histories of dark matter halos associated with dwarf satellites in Local Group galaxies and the resultant evolution of the baryonic component. Our model is based on the recently proposed property that the mean surface density of a dark halo inside a radius at maximum circular velocity V_max is universal over a large range of V_max. Following that this surface density of 20 Msun/pc^2 well explains dwarf satellites in the Milky Way and Andromeda, we find that the evolution of the dark halo in this common surface-density scale is characterized by the rapid increase of the halo mass assembled by the redshift z_TT of the tidal truncation by its host halo, at early epochs of z_TT >~ 6 or V_max <~ 22 km/s. This mass growth of the halo is slow at lower z_TT or larger V_max. Taking into account the baryon content in this dark halo evolution, under the influence of the ionizing background radiation, we find that the dwarf satellites are divided into roughly two families: those with V_max <~ 22 km/s having high star formation efficiency and those with larger V_max having less efficient star formation. This semi-analytical model is in good agreement with the high-resolution numerical simulation for galaxy formation and with the observed star formation histories for Fornax and Leo~II. This suggests that the evolution of a dark halo plays a key role in understanding star formation histories in dwarf satellites.

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 an analysis of the ALMA long baseline science verification data of the gravitational lens system SDP.81. We fit the positions of the brightest clumps at redshift z=3.042 and a possible AGN component of the lensing galaxy at redshift z=0.2999 in the band 7 continuum image using a canonical lens model, a singular isothermal ellipsoid plus an external shear. Then, we measure the ratio of fluxes in some apertures at the source plane where the lensed images are inversely mapped. We find that the aperture flux ratios of band 7 continuum image are perturbed by 10-20 percent with a significance at 2 ~ 3 sigma level. Moreover, we measure the astrometric shifts of multiply lensed images near the caustic using the CO(8-7) line. Using a lens model best-fitted to the band 7 continuum image, we reconstruct the source image of the CO(8-7) line by taking linear combination of inverted quadruply lensed images. At the 50th channel (rest-frame velocity 28.6 km/s) of the CO(8-7) line, we find an imprint of astrometric shifts of the order of 0.01 arcsec in the source image. Based on a semi-analytic calculation, we find that the observed anomalous flux ratios and the astrometric shifts can be explained by intergalactic dark structures in the line of sight. A compensated homogeneous spherical clump with a mean surface mass density of the order of 10^8 solar mass h^-1 arcsec^-2 can explain the observed anomaly and astrometric shifts simultaneously.

We investigate the non-spherical density structure of dark halos of the dwarf spheroidal (dSph) galaxies in the Milky Way and Andromeda galaxies based on revised axisymmetric mass models from our previous work. The models we adopt here fully take into account velocity anisotropy of tracer stars confined within a flattened dark halo. Applying our models to the available kinematic data of the 12 bright dSphs, we find that these galaxies associate with, in general, elongated dark halos, even considering the effect of this velocity anisotropy of stars. We also find that the best-fit parameters, especially for the shapes of dark halos and velocity anisotropy, are susceptible to both the availability of velocity data in the outer regions and the effect of the lack of sample stars in each spatial bin. Thus, to obtain more realistic limits on dark halo structures, we require photometric and kinematic data over much larger areas in the dSphs than previously explored. The results obtained from the currently available data suggest that the shapes of dark halos in the dSphs are more elongated than those of $\Lambda$CDM subhalos. This mismatch needs to be solved by theory including baryon components and the associated feedback to dark halos as well as by further observational limits in larger areas of dSphs. It is also found that more diffuse dark halos may have undergone consecutive star formation history, thereby implying that dark-halo structure plays an important role in star formation activity.

We investigate gas inflow and outflow histories in Milky Way-like disk galaxies, to get new insights into the baryonic processes in galaxy formation and evolution. For this purpose, we solve the equations for the evolutions of the surface mass densities of gas and metals at each radius in a galactic disk, based on the observed structural properties of distant star-forming galaxies, including the redshift evolution of their stellar mass distribution, their scaling relation between the mass of baryonic components, star formation rate (SFR) and chemical abundance, as well as the supposed evolution of their radial metallicity gradients (RMGs). We find that the efficiency of gas inflow for a given SFR decreases with time and that the inflow rate is always nearly proportional to the SFR. For gas outflow, although its efficiency for a given SFR is a decreasing function of time, similarly to gas inflow, the outflow rate is not necessarily proportional to the SFR and the relation between the outflow rate and SFR strongly depends on the evolution of the adopted RMG. We also find that the results on the outflow rate can be reproduced in the framework of momentum-driven (energy-driven) wind mechanism if the RMG is steepening (flattening) with time. Therefore if the well measured RMGs and their evolution for Milky Way-like galaxies are obtained from future observations, then our results will be useful to constrain the main driving mechanism for their galactic outflows.

The TMT Detailed Science Case describes the transformational science that the Thirty Meter Telescope will enable. Planned to begin science operations in 2024, TMT will open up opportunities for revolutionary discoveries in essentially every field of astronomy, astrophysics and cosmology, seeing much fainter objects much more clearly than existing telescopes. Per this capability, TMT's science agenda fills all of space and time, from nearby comets and asteroids, to exoplanets, to the most distant galaxies, and all the way back to the very first sources of light in the Universe. More than 150 astronomers from within the TMT partnership and beyond offered input in compiling the new 2015 Detailed Science Case. The contributing astronomers represent the entire TMT partnership, including the California Institute of Technology (Caltech), the Indian Institute of Astrophysics (IIA), the National Astronomical Observatories of the Chinese Academy of Sciences (NAOC), the National Astronomical Observatory of Japan (NAOJ), the University of California, the Association of Canadian Universities for Research in Astronomy (ACURA) and US associate partner, the Association of Universities for Research in Astronomy (AURA).

The tidal disruption of the Sagittarius dwarf Spheroidal galaxy (Sgr dSph) is producing the most prominent substructure in the Milky Way (MW) halo, the Sagittarius Stream. Aside from field stars, the Sgr dSph is suspected to have lost a number of globular clusters (GC). Many Galactic GC are suspected to have originated in the Sgr dSph. While for some candidates an origin in the Sgr dSph has been confirmed due to chemical similarities, others exist whose chemical composition has never been investigated. NGC 5053 and NGC 5634 are two among these scarcely studied Sgr dSph candidate-member clusters. To characterize their composition we analyzed one giant star in NGC 5053, and two in NGC 5634. We analize high-resolution and signal-to-noise spectra by means of the MyGIsFOS code, determining atmospheric parameters and abundances for up to 21 species between O and Eu. The abundances are compared with those of MW halo field stars, of "unassociated" MW halo globulars, and of the metal poor Sgr dSph main body population. We derive a metallicity of [FeII/H]=-2.26+-0.10 for NGC 5053, and of [FeI/H]=-1.99+-0.075 and -1.97+-0.076 for the two stars in NGC 5634. This makes NGC 5053 one of the most metal poor globular clusters in the MW. Both clusters display an alpha enhancement similar to the one of the halo at comparable metallicity. The two stars in NGC 5634 clearly display the Na-O anticorrelation widespread among MW globulars. Most other abundances are in good agreement with standard MW halo trends. The chemistry of the Sgr dSph main body populations is similar to the one of the halo at low metallicity. It is thus difficult to discriminate between an origin of NGC 5053 and NGC 5634 in the Sgr dSph, and one in the MW. However, the abundances of these clusters do appear closer to that of Sgr dSph than of the halo, favoring an origin in the Sgr dSph system.

In an attempt to place an explicit constraint on dark matter models, we define and estimate a mean surface density of a dark halo within a radius of maximum circular velocity, which is derivable for various galaxies with any dark-matter density profiles. We find that this surface density is generally constant across a wide range of maximum circular velocities of $\sim$ 10 to 400 km s$^{-1}$, irrespective of different density distribution in each of the galaxies. This common surface density at high halo-mass scales is found to be naturally reproduced by both cold and warm dark matter (CDM and WDM) models, even without employing any fitting procedures. However, the common surface density at dwarf-galaxy scales, for which we have derived from the Milky Way and Andromeda dwarf satellites, is reproduced only in a massive range of WDM particle masses, whereas CDM provides a reasonable agreement with the observed constancy. This is due to the striking difference between mass-concentration relations for CDM and WDM halos at low halo-mass scales. In order to explain the universal surface density of dwarf-galaxy scales in WDM models, we suggest that WDM particles need to be heavier than 3 keV.

We present the metallicity distribution of red giant branch (RGB) stars in M31's stellar halo, derived from photometric metallicity estimates for over 1500 spectroscopically confirmed RGB halo stars. The stellar sample comes from 38 halo fields observed with the Keck/DEIMOS spectrograph, ranging from 9 to 175 kpc in projected distance from M31's center, and includes 52 confirmed M31 halo stars beyond 100 kpc. While a wide range of metallicities is seen throughout the halo, the metal-rich peak of the metallicity distribution function becomes significantly less prominent with increasing radius. The metallicity profile of M31's stellar halo shows a continuous gradient from 9 to ~100 kpc, with a magnitude of -0.01 dex/kpc. The stellar velocity distributions in each field are used to identify stars that are likely associated with tidal debris features. The removal of tidal debris features does not significantly alter the metallicity gradient in M31's halo: a gradient is maintained in fields spanning 10 to 90 kpc. We analyze the halo metallicity profile, as well as the relative metallicities of stars associated with tidal debris features and the underlying halo population, in the context of current simulations of stellar halo formation. We argue that the large scale gradient in M31's halo implies M31 accreted at least one relatively massive progenitor in the past, while the field to field variation seen in the metallicity profile indicates that multiple smaller progenitors are likely to have contributed substantially to M31's outer halo.

We present new limits on the global shape of the dark matter halo in the Andromeda galaxy using and generalizing non-spherical mass models developed by Hayashi & Chiba and compare our results with theoretical predictions of cold dark matter (CDM) models. This is motivated by the fact that CDM models predict non-spherical virialized dark halos, which reflect the process of mass assembly in the galactic scale. Applying our models to the latest kinematic data of globular clusters and dwarf spheroidal galaxies in the Andromeda halo, we find that the most plausible cases for Andromeda yield a prolate shape for its dark halo, irrespective of assumed density profiles. We also find that this prolate dark halo in Andromeda is consistent with theoretical predictions in which the satellites are distributed anisotropically and preferentially located along major axes of their host halos. It is a reflection of the intimate connection between galactic dark matter halos and the cosmic web. Therefore, our result is profound in understanding internal dynamics of halo tracers in Andromeda, such as orbital evolutions of tidal stellar streams, which play important roles in extracting the abundance of CDM subhalos through their dynamical effects on stream structures.

We study the detailed properties of the radial metallicity gradient in the stellar disk of our Galaxy to constrain its chemical and structural evolution. For this purpose we select and analyze $\sim$ 18,500 disk stars taken from two datasets, the Sloan Digital Sky Survey (SDSS) and the High-Accuracy Radial velocity Planetary Searcher (HARPS). On these surveys we examine the metallicity gradient, $\Delta$[Fe/H]/$\Delta R_{\rm g}$, along the guiding-center radii, $R_{\rm g}$, of stars and its dependence on the [$\alpha$/Fe] ratios, to infer the original metallicity distribution of the gas disk from which those stars formed and its time evolution. In both sample sources, the thick-disk candidate stars characterized by high [$\alpha$/Fe] ratios ([$\alpha$/Fe] $>$ 0.3 in SDSS, [$\alpha$/Fe] $>$ 0.2 in HARPS) are found to show a positive $\Delta$[Fe/H]/$\Delta R_{\rm g}$, whereas the thin-disk candidate stars characterized by lower [$\alpha$/Fe] ratios show a negative one. Furthermore, we find that the relatively young thin-disk population characterized by much lower [$\alpha$/Fe] ratios ([$\alpha$/Fe] $<$ 0.2 in SDSS, [$\alpha$/Fe] $<$ 0.1 in HARPS) shows notably a flattening $\Delta$[Fe/H]/$\Delta R_{\rm g}$ with decreasing [$\alpha$/Fe], in contrast to the old one with higher [$\alpha$/Fe] ratios ([$\alpha$/Fe] $\sim$ 0.2 in SDSS, [$\alpha$/Fe] $\sim$ 0.1 in HARPS). The possible implication for early disk evolution is discussed, in the context of galaxy formation accompanying the rapid infall of primordial gas on the inner disk region, which can generate a positive metallicity gradient, and the subsequent chemical evolution of the disk, which results in a flattening effect of a metallicity gradient at later epochs.

The Subaru Prime Focus Spectrograph (PFS) is a massively-multiplexed fiber-fed optical and near-infrared 3-arm spectrograph (N_fiber=2400, 380<lambda<1260nm, 1.3 degree diameter FoV), offering unique opportunities in survey astronomy. Here we summarize the science case feasible for a survey of Subaru 300 nights. We describe plans to constrain the nature of dark energy via a survey of emission line galaxies spanning a comoving volume of 9.3 (Gpc/h)^3 in the redshift range 0.8<z<2.4. In each of 6 redshift bins, the cosmological distances will be measured to 3% precision via BAO, and redshift-space distortions will be used to constrain structure growth to 6% precision. In the GA program, radial velocities and chemical abundances of stars in the Milky Way and M31 will be used to infer the past assembly histories of spiral galaxies and the structure of their dark matter halos. Data will be secured for 10^6 stars in the Galactic thick-disk, halo and tidal streams as faint as V~22, including stars with V < 20 to complement the goals of the Gaia mission. A medium-resolution mode with R = 5000 to be implemented in the red arm will allow the measurement of multiple alpha-element abundances and more precise velocities for Galactic stars, elucidating the detailed chemo-dynamical structure and evolution of each of the main stellar components of the Milky Way Galaxy and of its dwarf spheroidal galaxies. For the extragalactic program, our simulations suggest the wide avelength range will be powerful in probing the galaxy population and its clustering over a wide redshift range. We propose to conduct a color-selected survey of 1<z<2 galaxies and AGN over 16 deg^2 to J~23.4, yielding a fair sample of galaxies with stellar masses above ~10^{10}Ms at z~2. A two-tiered survey of higher redshift LBGs and LAEs will quantify the properties of early systems close to the reionization epoch.

We have found a radio-wave-reflection effect in rock salt for the detection of ultra-high energy neutrinos which are expected to be generated in Greisen, Zatsepin, and Kuzmin (GZK) processes in the universe. When an UHE neutrino interacts with rock salt or ice as a detection medium, a shower is generated. That shower is formed by hadronic and electromagnetic avalanche processes. The energy of the UHE neutrino shower converts to thermal energy through ionization processes. Consequently, the temperature rises along the shower produced by the UHE neutrino. The refractive index of the medium rises with temperature. The irregularity of the refractive index in the medium leads to a reflection of radio waves. This reflection effect combined with the long attenuation length of radio waves in rock salt and ice would yield a new method to detect UHE neutrinos. We measured the phase of the reflected radio wave under irradiation with an electron beam on ice and rock salt powder. The measured phase showed excellent consistence with the power reflection fraction which was measured directly. A model taking into account the temperature change explained the phase and the amplitude of the reflected wave. Therefore the reflection mechanism was confirmed. The power reflection fraction was compared with that calculated with the Fresnel equations, the ratio between the measured result and that obtained with the Fresnel equations in ice was larger than that of rock salt.

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 the surface brightness profile of M31's stellar halo out to a projected radius of 175 kpc. The surface brightness estimates are based on confirmed samples of M31 red giant branch stars derived from Keck/DEIMOS spectroscopic observations. A set of empirical spectroscopic and photometric M31 membership diagnostics is used to identify and reject foreground and background contaminants. This enables us to trace the stellar halo of M31 to larger projected distances and fainter surface brightnesses than previous photometric studies. The surface brightness profile of M31's halo follows a power-law with index -2.2 +/- 0.2 and extends to a projected distance of at least ~175 kpc (~ 2/3 of M31's virial radius), with no evidence of a downward break at large radii. The best-fit elliptical isophotes have b/a=0.94 with the major axis of the halo aligned along the minor axis of M31's disk, consistent with a prolate halo, although the data are also consistent with M31's halo having spherical symmetry. The fact that tidal debris features are kinematically cold is used to identify substructure in the spectroscopic fields out to projected radii of 90 kpc, and investigate the effect of this substructure on the surface brightness profile. The scatter in the surface brightness profile is reduced when kinematically identified tidal debris features in M31 are statistically subtracted; the remaining profile indicates a comparatively diffuse stellar component to M31's stellar halo exists to large distances. Beyond 90 kpc, kinematically cold tidal debris features can not be identified due to small number statistics; nevertheless, the significant field-to-field variation in surface brightness beyond 90 kpc suggests that the outermost region of M31's halo is also comprised to a significant degree of stars stripped from accreted objects.

We construct axisymmetric mass models for dwarf spheroidal (dSph) galaxies in the Milky Way to obtain plausible limits on the non-spherical structure of their dark halos. This is motivated by the fact that the observed luminous parts of the dSphs are actually non-spherical and Cold Dark Matter (CDM) models predict non-spherical virialized dark halos. Our models consider velocity anisotropy of stars $\bar{v^2_R} / \bar{v^2_{\phi}}$, which can vary with the adopted cylindrical coordinates under the assumption $\bar{v^2_z}=\bar{v^2_R}$ for simplicity, and also include an inclination of the system as a fitting parameter to explain the observed line-of-sight velocity dispersion profile. Applying these models to six of the bright dSphs in the Milky Way, we find that the best-fitting cases for most of the dSphs yield oblate and flattened dark halos, irrespective of assumed density profiles in their central parts. We also find that the total mass of the dSphs enclosed within a spheroid with major-axis length of 300 pc varies from $10^6M_{\odot}$ to $10^7M_{\odot}$, contrary to the conclusion from spherical models. This suggests the importance of considering shapes of dark halos in mass models of the dSphs. It is also found that dark halos of the Galactic dSphs may be more flattened than N-body predictions, thereby implying our yet incomplete understanding of baryonic and/or non-baryonic dark matter physics in dwarf galaxy scales.

The Andromeda galaxy (M31) shows many tidal features in its halo, including the Giant Southern Stream (GSS) and a sharp ledge in surface density on its western side (the W Shelf). Using DEIMOS on the Keck telescope, we obtain radial velocities of M31's giant stars along its NW minor axis, in a radial range covering the W Shelf and extending beyond its edge. In the space of velocity versus radius, the sample shows the wedge pattern expected from a radial shell, which is detected clearly here for the first time. This confirms predictions from an earlier model of formation of the GSS, which proposed that the W Shelf is a shell from the third orbital wrap of the same tidal debris stream that produces the GSS, with the main body of the progenitor lying in the second wrap. We calculate the distortions in the shelf wedge pattern expected from its outward expansion and angular momentum, and show that these effects are echoed in the data. In addition, a hot, relatively smooth spheroid population is clearly present. We construct a bulge-disk-halo N-body model that agrees with surface brightness and kinematic constraints, and combine it with a simulation of the GSS. From the contrasting kinematic signatures of the hot spheroid and shelf components, we decompose the observed stellar metallicity distribution into contributions from each component using a non-parametric mixture model. The shelf component's metallicity distribution matches previous observations of the GSS superbly, further strengthening the evidence they are connected and bolstering the case for a massive progenitor of this stream.

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.

Carollo et al. have recently resolved the stellar population of the Milky Way halo into at least two distinct components, an inner halo and an outer halo. This result has been criticized by Schoenrich et al., who claim that the retrograde signature associated with the outer halo is due to the adoption of faulty distances. We refute this claim, and demonstrate that the Schoenrich et al. photometric distances are themselves flawed because they adopted an incorrect main-sequence absolute magnitude relationship from the work of Ivezi\'c et al. When compared to the recommended relation from Ivezi\'c et al., which is tied to a Milky Way globular cluster distance scale and accounts for age and metallicity effects, the relation adopted by Schoenrich et al. yields up to 18% shorter distances for stars near the main-sequence turnoff (TO). Use of the correct relationship yields agreement between the distances assigned by Carollo et al. and Ivezi\'{c} et al. for low-metallicity dwarfs to within 6-10%. Schoenrich et al. also point out that intermediate-gravity stars (3.5 <= log g <= 4.0) with colors redder than the TO region are likely misclassified, with which we concur. We implement a new procedure to reassign luminosity classifications for the TO stars that require it. New derivations of the rotational behavior demonstrate that the retrograde signature and high velocity dispersion of the outer-halo population remains. We summarize additional lines of evidence for a dual halo, including a test of the retrograde signature based on proper motions alone, and conclude that the preponderance of evidence strongly rejects the single-halo interpretation.

As part of our survey of galactic stellar halos, we investigate the structure and stellar populations of the northern outer part of the stellar halo in NGC55, a member galaxy of the Sculptor Group, using deep and wide-field V- and I-band images taken with Subaru/Suprime-Cam. Based on the analysis of the color-magnitude diagrams (CMDs) for red-giant-branch (RGB) stars, we derive a tip of RGB (TRGB)-based distance modulus to the galaxy of (m-M)_0 = 26.58 +/- 0.11 (d = 2.1 +/- 0.1 Mpc). From the stellar density maps, we detect the asymmetrically disturbed, thick disk structure and two metal-poor overdense substructures in the north region of NGC55, which may correspond to merger remnants associated with hierarchical formation of NGC55's halo. In addition, we identify a diffuse metal-poor halo extended out to at least z ~ 16 kpc from the galactic plane. The surface-brightness profiles toward the z-direction perpendicular to the galactic plane suggest that the stellar density distribution in the northern outer part of NGC55 is described by a locally isothermal disk at z <~ 6 kpc and a likely diffuse metal-poor halo with V-band surface brightness of mu_V >~ 32 mag arcsec^{-2}, where old RGB stars dominate. We derive the metallicity distributions (MDs) of these structures on the basis of the photometric comparison of RGB stars with the theoretical stellar evolutionary models. The MDs of the thick disk structures show the peak and mean metallicity of [Fe/H]peak ~ -1.4 and [Fe/H]mean ~ -1.7, respectively, while the outer substructures show more metal-poor features than the thick disk structure. Combined with the current results with our previous study for M31's halo, we discuss the possible difference in the formation process of stellar halos among different Hubble types.

The structure and kinematics of the recognized stellar components of the Milky Way are explored, based on well-determined atmospheric parameters and kinematic quantities for 32360 calibration stars from the Sloan Digital Sky Survey (SDSS) and its first extension, (SDSS-II), which included the sub-survey SEGUE: Sloan Extension for Galactic Understanding and Exploration. Full space motions for a sub-sample of 16920 stars, exploring a local volume within 4 kpc of the Sun, are used to derive velocity ellipsoids for the inner- and outer-halo components of the Galaxy, as well as for the canonical thick-disk and proposed metal-weak thick-disk populations. We first examine the question of whether the data require the presence of at least a two-component halo in order to account for the rotational behavior of likely halo stars in the local volume, and whether more than two components are needed. We also address the question of whether the proposed metal-weak thick disk is kinematically and chemically distinct from the canonical thick disk. In addition, we consider the fractions of each component required to understand the nature of the observed kinematic behavior of the stellar populations of the Galaxy as a function of distance from the plane. Scale lengths and scale heights for the thick-disk and metal-weak thick-disk components are determined. Spatial density profiles for the inner- and outer-halo populations are inferred from a Jeans Theorem analysis. The full set of calibration stars (including those outside the local volume) is used to test for the expected changes in the observed stellar metallicity distribution function with distance above the Galactic plane in-situ, due to the changing contributions from the underlying stellar populations. [abridged]

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.

We present a photometric survey of the stellar halo of the Andromeda galaxy, using Suprime-Cam on the Subaru Telescope. A detailed analysis of VI color-magnitude diagrams of the resolved stellar population is used to measure properties such as line-of-sight distance, surface brightness, metallicity, and age, and these are used to isolate and characterize different components of the M31 halo: (1) several substructures, and (2) the smooth halo. First, we study M31's halo substructure along the north-west/south-east minor axis out to R ~ 100 kpc and the south-west major axis region at R ~ 60 kpc. We confirm two substructures in the south-east halo reported by Ibata et al. (2007) and discover two overdense substructures in the north-west halo. We investigate the properties of these four substructures as well as other structures including the western shelf and find that differences in stellar populations among these systems, thereby suggesting each has a different origin. Our statistical analysis implies that the M31 halo as a whole may contain at least 16 substructures, each with a different origin. Second, we investigate the properties of an underlying, smooth and extended halo component out to R > 100 kpc. We find that the surface density of this smooth halo can be fitted to a Hernquist model of scale radius ~ 17 kpc or a power-law profile with ~ R^{-2.17 +/- 0.15}. In contrast to the relative smoothness of the halo density profile, its metallicity distribution appears to be spatially non-uniform with non-monotonic variations with radius, suggesting that the halo population has not had sufficient time to dynamically homogenize the accreted populations. Further implications for the formation of the M31 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.

The presence of the lithium-6 isotope in some metal-poor stars is a matter of surprise because of the high values observed. Non-standard models of Big Bang nucleosynthesis and pre-Galactic cosmic ray fusion and spallation have been proposed to explain these values. However, the observations of this light isotope are challenging which may make some detections disputable. The goal was to determine 6Li/7Li for a sample of metal-poor stars; three of them have been previously studied and the remaining two are new for this type of study. The purpose was to increase, if possible, the number of lithium-6 detections and to confirm previously published results. Spectra of the resonance doublet line of neutral lithium Li I 670.78nm were taken with the High Dispersion Spectrograph at the Subaru 8.2m-telescope for a sample of five metal-poor stars (-3.12 < [Fe/H] < -2.19). The contribution of lithium-6 to the total observed line profile was estimated from the 1D-LTE analysis of the line asymmetry. Observed asymmetries could be reproduced assuming isotopic abundance ratios 6Li/7Li of the order of: 0.004 for BD+26 3578, ~ 0.010 for BD+02 3375 and G 64-37, 0.025 for BD+20 3603 and 0.047 for BD-04 3208. We found that these results were very sensitive to several of the assumptions made in the analysis, in particular, the treatment of the residual structure in the analysed spectra. Our final estimates for the errors are respectively delta_6Li/7Li=+- 0.028, 0.029, 0.039, 0.025 and 0.039. The 6Li/7Li ratios for the sample are comparable to or even lower than these error values, so that detections of lithium-6 can not safely be claimed despite of the high resolving power (R ~ 95 000) and S/N (400-600).

We investigate the dynamical effects of dark matter subhalos on the structure and evolution of a galactic disk, using semi-analytic method that includes approximated and empirical relations as achieved in detailed numerical simulations of the cold dark matter model. We calculate the upper limit for the size of a galactic disk at a specific redshift $z$, based on the orbital properties of subhalos characterized by their pericentric distances from the center of a host halo. We find that this possibly largest size of a disk as determined by the smallest pericentric distances of subhalos shows the characteristic properties, which are basically in agreement with an observed galactic disk at low and high $z$. Namely, it is found that a massive disk can have a larger size than a less massive one, because of its stability against the destruction effect of subhalos. Also, with fixed mass, the size of a galactic disk at low $z$ can be larger than that at high $z$, reflecting the orbital evolution of subhalos with respect to a host halo. These results suggest that the presence and structure of a galactic disk may be dynamically limited by the interaction with dark matter substructures, especially at high $z$.

We investigate structural, kinematical, and chemical properties of stars and gas in the Small Magellanic Cloud (SMC) interacting with the Large Magellanic Cloud (LMC) and the Galaxy based on a series of self-consistent chemodynamical simulations. We adopt a new "dwarf spheroidal model" in which the SMC initially has both old stars with a spherical spatial distribution and an extended HI gas disk. We mainly investigate SMC's evolution for the last 3 Gyr within which the Magellanic stream (MS) and the Magellanic bridge (MB) can be formed as a result of the LMC-SMC-Galaxy interaction. Our principal results, which can be tested against observations, are as follows. The final spatial distribution of the old stars projected onto the sky is spherical even after the strong LMC-SMC-Galaxy interaction, whereas that of the new ones is significantly flattened and appears to form a bar structure. Old stars have the line-of-sight velocity dispersion (sigma) of ~ 30 km/s and slow rotation with the maximum rotational velocity (V) of less than slow rotation with the maximum rotational velocity (V) of less than 20 km/s and show asymmetry in the radial profiles. New stars have a smaller sigma than old ones and a significant amount of rotation (V/sigma >1). HI gas shows velocity dispersions of sigma = 10-40 km/s a high maximum rotational velocity (V ~ 50 km/s), and the spatial distribution similar to that of new stars. The new stars with ages younger than 3 Gyr show a negative metallicity gradient in the sense that more metal-rich stars are located in the inner regions of the SMC.

We present mid-infrared imaging at 11.7 $\mu$m for the quadruple lens systems, MG0414+0534 and Q2237+030, using the cooled mid-infrared camera and spectrometer (COMICS) attached on the Subaru telescope. MG0414+0534 is characterized by a bright pair of lensed images (A1, A2) and their optical flux ratio A2/A1 deviates significantly from the prediction of a smooth lens model. Q2237+030 is comprised of four lensed images, which are significantly affected by microlensing in a foreground lensing galaxy. Our mid-infrared observations of these lensed images have revealed that the mid-infrared flux ratio for A2/A1 of MG0414+0534 is nearly unity (0.90 +- 0.04). We find that this flux ratio is systematically small, at 4 to 5 sigma level, compared with the prediction of a best smooth lens model (1.09) represented by a singular isothermal ellipsoid and external shear, thereby suggesting the presence of substructures to explain our observational result. In contrast, for Q2237+030, our high signal-to-noise observation indicates that the mid-infrared flux ratios between all the four images of Q2237+030 are virtually consistent with the prediction of a smooth lens model. Based on the size estimate of the dust torus surrounding the nuclei of these QSOs, we set limits on the mass of a substructure in these lens systems, which can cause anomalies in the flux ratios. For MG0414+0534, since the required mass of a substructure inside its Einstein radius is > 360 M_sun, millilensing by a CDM substructure is most likely. If it is modeled as a singular isothermal sphere, the mass inside radius of 100 pc is given as > 1.0x10^5 M_sun. For Q2237+030, there is no significant evidence of millilensing, so the reported anomalous flux ratios in shorter wavelengths are entirely caused by microlensing by stars (abridged).

We demonstrate for the first time that gaseous halos of disk galaxies can play a vital role in recycling metal-rich gas ejected from the bulges and thus in promoting chemical evolution of disks. Our numerical simulations show that metal-rich stellar winds from bulges in disk galaxies can be accreted onto the thin disks owing to hydrodynamical interaction between the gaseous ejecta and the gaseous halos, if the mean densities of the halos (rho_ hg) are as high as 10^{-5} cm^{-3}. The total amount of gas that is ejected from a bulge through a stellar wind and then accreted onto the disk depends mainly on rho_ hg and the initial velocity of the stellar wind. About ~ 1% of gaseous ejecta from bulges in disk galaxies of scale length a_d can be accreted onto disks around R ~ 2.5 a_ d for a reasonable set of model parameters. We discuss these results in the context of the origin of the surprisingly high metallicities of the solar neighborhood disk stars in the Galaxy. We also discuss some implications of the present results in terms of chemical evolution of disk galaxies with possibly different rho_ hg in different galaxy environments.

Using the HST ACS, we have obtained deep optical images reaching well below the oldest main sequence turnoff in fields on the southeast minor-axis of the Andromeda Galaxy, 35 kpc from the nucleus. These data probe the star formation history in the extended halo of Andromeda -- that region beyond 30 kpc that appears both chemically and morphologically distinct from the metal-rich, highly-disturbed inner spheroid. The present data, together with our previous data for fields at 11 and 21 kpc, do not show a simple trend toward older ages and lower metallicities, as one might expect for populations further removed from the obvious disturbances of the inner spheroid. Specifically, the mean ages and [Fe/H] values at 11 kpc, 21 kpc, and 35 kpc are 9.7 Gyr and -0.65, 11.0 Gyr and -0.87, and 10.5 Gyr and -0.98, respectively. In the best-fit model of the 35 kpc population, one third of the stars are younger than 10 Gyr, while only ~10% of the stars are truly ancient and metal-poor. The extended halo thus exhibits clear evidence of its hierarchical assembly, and the contribution from any classical halo formed via early monolithic collapse must be small.

Recent observations of the Small Magellanic Cloud (SMC) have revealed that the HI gas shows a significant amount of rotation (V_c 60 km/s), while no or little rotation is evident for the old stellar populations. We suggest that this unique kinematical difference between these components in the SMC can be caused by a major merger event which occurred in the early stage of the SMC formation. Our simulations show that dissipative dwarf-dwarf merging can transform two gas-rich dwarf irregulars into a new dwarf, which consists of a spheroidal stellar component and a rotating extended HI disk. The remnant of this dwarf-dwarf merging shows significantly different kinematics between stars and gas, in the sense that a gas disk rotates rapidly while a stellar component shows little rotation. We thus suggest that the simulated dwarf having a dynamically hot spheroid and an extended gas disk finally evolves into the present SMC after efficient stripping of the outer gas via tidal fields of the Galaxy and the Large Magellanic Cloud. We also suggest that spatial distributions and kinematics of RGB and AGB stars with different ages in the possible spheroidal component of the SMC can provide valuable information on whether and when a past major merger event really occurred in the SMC.

To set useful limits on the abundance of small-scale dark matter halos (subhalos) in a galaxy scale, we have carried out mid-infrared imaging and integral-field spectroscopy for a sample of quadruple lens systems showing anomalous flux ratios. These observations using Subaru have been successful for distinguishing millilensing by subhalos from microlensing by stars. Current status for our lensing analysis of dark matter substructure is reported.

Hierarchical clustering models of cold dark matter (CDM) predict that about 5% - 10% of a galaxy-sized halo with mass ~ 10^12 solar masses (M_sun) resides in substructures (CDM subhalos) with masses <= 10^8 M_sun. To directly identify such substructures, we propose to observe radio continuum emission from multiply imaged QSOs using VSOP-2 with a high angular resolution.

The halo of the Milky Way provides unique elemental abundance and kinematic information on the first objects to form in the Universe, which can be used to tightly constrain models of galaxy formation and evolution. Although the halo was once considered a single component, evidence for its dichotomy has slowly emerged in recent years from inspection of small samples of halo objects. Here we show that the halo is indeed clearly divisible into two broadly overlapping structural components -- an inner and an outer halo -- that exhibit different spatial density profiles, stellar orbits and stellar metallicities (abundances of elements heavier than helium). The inner halo has a modest net prograde rotation, whereas the outer halo exhibits a net retrograde rotation and a peak metallicity one-third that of the inner halo. These properties indicate that the individual halo components probably formed in fundamentally different ways, through successive dissipational (inner) and dissipationless (outer) mergers and tidal disruption of proto-Galactic clumps.

We show that collisions between the outer Galactic HI disk and the leading arms (LAs) of the Magellanic stream (MS) can create giant HI holes and chimney-like structures in the disk. Based on the results of our N-body simulations on the last 2.5 Gyr evolution of the Large and Small Magellanic Clouds (LMC and SMC, respectively) interacting with the Galaxy, we investigate when and where the LAs can pass through the Galactic plane after the MS formation. We then investigate hydrodynamical interaction between LAs and the Galactic HI disk (``the Magellanic impact'') by using our new hydrodynamical simulations with somewhat idealized models of the LAs. We find that about 1-3% of the initial gas mass of the SMC, which consists of the LAs, can pass through the outer part (R=20-35 kpc) of the Galactic HI disk about 0.2 Gyr ago. We also find that the Magellanic impact can push out some fraction (~1%) of the outer Galactic HI disk to form 1-10 kpc-scale HI holes and chimney-like bridges between the LAs and the disk.

Existence of GZK neutrinos (ultra high energy neutrinos) have been justified although the flux is very low. A new method is desired to use a huge mass of a detector medium to detect them. A fundamental study of radar method was carried out to measure microwave reflection from electromagnetic energy deposit by X-ray irradiation in a small rock salt sample. The reflection rate of 1x10^-6 was found at the energy deposit of 1x10^19 eV which was proportional to square of the X-ray intensity suggesting the effect to be coherent scattering. The decay time of the reflection was several seconds. This effect implies a large scale natural rock salt formation could be utilized like a bubble chamber irradiated by radio wave instead of visible light to detect GZK neutrinos.

We first show that a large amount of metal-poor gas is stripped from the Small Magellanic Cloud (SMC) and fallen into the Large Magellanic Cloud (LMC) during the tidal interaction between the SMC, the LMC, and the Galaxy over the last 2 Gyrs. We propose that this metal-poor gas can closely be associated with the origin of LMC's young and intermediate-age stars and star clusters with distinctively low-metallicities with [Fe/H] < -0.6. We numerically investigate whether gas initially in the outer part of the SMC's gas disk can be stripped during the LMC-SMC-Galaxy interaction and consequently can pass through the central region (R<7.5 kpc) of the LMC. We find that about 0.7 % and 18 % of the SMC's gas can pass through the central region of the LMC about 1.3 Gyr ago and 0.2 Gyr ago, respectively. The possible mean metallicity of the replenished gas from the SMC to LMC is about [Fe/H] = -0.9 to -1.0 for the two interacting phases. These results imply that the LMC can temporarily replenish gas supplies through the sporadic accretion and infall of metal-poor gas from the SMC. These furthermore imply that if these gas from the SMC can collide with gas in the LMC to form new stars in the LMC, the metallicities of the stars can be significantly lower than those of stars formed from gas initially within the LMC.

We first present chemodynamical simulations to investigate how stellar winds of massive stars influence early dynamical and chemical evolution of forming globular clusters (GCs). In our numerical models, GCs form in turbulent,high-density giant molecular clouds (GMCs), which are embedded in a massive dark matter halo at high redshifts. We show how high-density, compact stellar systems are formed from GMCs influenced both by physical processes associated with star formation and by tidal fields of their host halos. We also show that chemical pollution of GC-forming GMCs by stellar winds from massive stars can result in star-to-star abundance inhomogeneities among light elements (e.g., C, N, and O) of stars in GCs. The present model with a canonical initial mass function (IMF) also shows a C-N anticorrelation that stars with smaller [C/Fe] have larger [N/Fe] in a GC. Although these results imply that ``self-pollution'' of GC-forming GMCs by stellar winds from massive stars can cause abundance inhomogeneities of GCs, the present models with different parameters and canonical IMFs can not show N-rich stars with [N/Fe] ~ 0.8 observed in some GCs (e.g., NGC 6752). We discuss this apparent failure in the context of massive star formation preceding low-mass one within GC-forming GMCs (``bimodal star formation scenario''). We also show that although almost all stars (~97%) show normal He abundances (Y) of ~0.24 some stars later formed in GMCs can have Y as high as ~0.3 in some models. The number fraction of He-rich stars with Y >0.26 is however found to be small (~10^-3) for most models.

We have used Suprime-Cam on the Subaru Telescope to conduct a V- and I-band imaging survey of fields sampling the spheroid of the Andromeda galaxy along its south-east minor axis. Our photometric data are deep enough to resolve stars down to the red clump. Based on a large and reliable sample of red giant stars available from this deep wide-field imager, we have derived metallicity distributions vs. radius and a surface brightness profile over projected distances of R=23-66 kpc from the galaxy's center. The metallicity distributions across this region shows a clear high mean metallicity and a broad distribution ([Fe/H] ~ -0.6 +/- 0.5), and indicates no metallicity gradient within our observed range. The surface brightness profile at R>40 kpc is found to be flatter than previously thought. It is conceivable that this part of the halo samples as yet unidentified, metal-rich substructure.

We investigate the present distributions of gas and young stars in the Large and Small Magellanic Clouds (LMC and SMC) based on fully self-consistent numerical simulations of the Clouds for the last ~ 0.8 Gyr. Our principal results, which can be tested against observations, are as follows. The last dynamical and hydrodynamical interaction between the Clouds about ~ 0.2 Gyr ago can form the apparently off-center bar and peculiar HI spirals of the LMC. The present spatial distributions of young stars with ages less than 20 Myr in the LMC can be significantly asymmetric and clumpy owing to the interaction. A small but non-negligible fraction of stellar and gaseous components can be transferred from the SMC into the LMC during the interaction to form diffuse halo components around the LMC. The burst of star formation in the SMC can be synchronized with that of the LMC about 0.2 Gyr ago in some models. New stars can form from gas in the SMC's tidal tails, one of which can be observed as the Magellanic Bridge (MB). The metallicity distribution function of new stars in the MB has a peak of [Fe/H] ~ -0.8, which is significantly smaller than the stellar metallicity of the SMC. Based on these results, we discuss the origin of 30 Doradus, the southern molecular ridge of the LMC, the globular cluster ESO 121-SC03, metal-poor inter-Cloud stars within the MB, and giant HI holes of the LMC.

We have observed the quadruply lensed quasar 1RXS J1131-1231 with the integral field spectrograph mode of the Kyoto Tridimensional Spectrograph II mounted on the Subaru telescope. Its field of view has covered simultaneously the three brighter lensed images A, B, and C, which are known to exhibit anomalous flux ratios in their continuum emission. We have found that the [OIII] line flux ratios among these lensed images are consistent with those predicted by smooth-lens models. The absence of both microlensing and millilensing effects on this [OIII] narrow line region sets important limits on the mass of any substructures along the line of sight, which is expressed as M_E < 10^5 M_solar for the mass inside an Einstein radius. In contrast, the H_beta line emission, which originates from the broad line region, shows an anomaly in the flux ratio between images B and C, i.e., a factor two smaller C/B ratio than predicted by smooth-lens models. The ratio of A/B in the H_beta line is well reproduced. We show that the anomalous C/B ratio for the H_beta line is caused most likely by micro/milli-lensing of image C. This is because other effects, such as the differential dust extinction and/or arrival time difference between images B and C, or the simultaneous lensing of another pair of images A and B, are all unlikely. In addition, we have found that the broad H_beta line of image A shows a slight asymmetry in its profile compared with those in the other images, which suggests the presence of a small microlensing effect on this line emitting region of image A.

We investigate the dynamical interaction between a galactic disk and surrounding numerous dark subhalos as expected for a galaxy-sized halo in the cold dark matter (CDM) models. Our particular interest is to what extent accretion events of subhalos into a disk are allowed in light of the observed thinness of a disk. Several models of subhalos are considered in terms of their internal density distribution, mass function, and spatial and velocity distributions. Based on a series of N-body simulations, we find that the disk thickening quantified by the change of its scale height, Delta z_d, depends strongly on the individual mass of an interacting subhalo M_{sub}. This is described by the relation, Delta z_d / R_d = 8 Sum_{j=1}^N (M_{sub,j}/M_d)**2, where R_d is a disk scale length, M_d is a disk mass, and N is the total number of accretion events of subhalos inside a disk region (< 3R_d). Using this relation, we find that an observed thin disk has not ever interacted with subhalos with the total mass of more than 15% disk mass. Also, a less massive disk with smaller circular velocity V_c is more affected by subhalos than a disk with larger V_c, in agreement with the observation. Further implications of our results for the origin of a thick disk component are also discussed.

A galaxy-sized halo may contain a large number of intermediate mass (10^{2-4} solar mass) compact objects (IMCOs), which can be intermediate mass black holes (IMBHs) or the CDM subhalos. We propose to directly detect the IMBHs by observing multiply imaged QSO-galaxy lens systems with a high angular resolution (~ 0.03mas), which would be achieved by the next-VLBI space missions. The silhouette of the IMBHs would appear as an either monopole-like or dipole-like variation at the scale of the Einstein radius against the QSO jets. As a byproduct, we can also directly detect the 10^{4-5} solar mass CDM subhalos. From a measurement of the local distortion in the surface brightness of the QSO jet, we can make a distinction between a point mass (corresponding to an IMBH) and an extended structure (corresponding to a CDM subhalo). It would be a unique probe of the IMCOs whose nature has been under the veil of mistery.

Recent observations have discovered star formation activities in the extreme outer regions of disk galaxies. However it remains unclear what physical mechanisms are responsible for triggering star formation in such low-density gaseous environments of galaxies. In order to understand the origin of these outer star-forming regions, we numerically investigate how the impact of dark matter subhalos orbiting a gas-rich disk galaxy embedded in a massive dark matter halo influences the dynamical evolution of outer HI gas disk of the galaxy. We find that if the masses of the subhalos ($M_{\rm sb}$) in a galaxy with an extended HI gas disk are as large as $10^{-3} \times M_{\rm h}$, where $M_{\rm h}$ is the total mass of the galaxy's dark halo, local fine structures can be formed in the extended HI disk. We also find that the gas densities of some apparently filamentary structures can exceed a threshold gas density for star formation and thus be likely to be converted into new stars in the outer part of the HI disk in some models with larger $M_{\rm sb}$. These results thus imply that the impact of dark matter subhalos (``dark impact'') can be important for better understanding the origin of recent star formation discovered in the extreme outer regions of disk galaxies. We also suggest that characteristic morphologies of local gaseous structures formed by the dark impact can indirectly prove the existence of dark matter subhalos in galaxies. We discuss the origin of giant HI holes observed in some gas-rich galaxies (e.g., NGC 6822) in the context of the dark impact.

We investigate the extended source size effects on gravitational lensing in which a lens consists of a smooth potential and small mass clumps (``substructure lensing''). We first consider a lens model that consists of a clump modeled as a singular isothermal sphere (SIS) and a primary lens modeled as an external background shear and convergence. For this simple model, we derive analytic formulae for (de)magnification of circularly symmetric top-hat sources with three types of parity for their lensed images, namely, positive, negative, and doubly negative parities. Provided that the source size is sufficiently larger than the Einstein radius of the SIS, we find that in the positive (doubly negative) parity case, an extended source is always magnified (demagnified) in comparison with the unperturbed macrolens system, whereas in the negative parity case, the (de)magnification effect, which depends on the sign of convergence minus unity is weaker than those in other parities. It is shown that a measurement of the distortion pattern in a multiply lensed image enables us to break the degeneracy between the lensing effects of clump mass and those of clump distance if lensing parameters of the relevant macrolens model are determined from the position and flux of multiple images. We also show that an actual density profile of a clump can be directly measured by analyzing the ``fine structure'' in a multiply lensed image within the Einstein radius of the clump.

The cold dark matter (CDM) structure formation model predicts that about 5-10 percent of a typical galactic halo of mass $\sim 10^{12} \ms$ is in substructures with masses $\lesssim 10^8 \ms$. To directly detect such substructures, we propose to observe dust continuum emission from a strongly lensed QSO-host galaxy using a large submillimeter interferometer array with a high angular resolution of $\sim 0.01$arcsec such as the planned Atacama Large Submillimeter Array (ALMA). To assess their observational feasibility, we numerically simulate millilensing of an extended circular source by a CDM substructure modeled as a tidally truncated singular isothermal sphere (SIS) embedded in a typical QSO-galaxy lens system, B1422+231, modeled as a singular isothermal ellipsoid (SIE) with an external constant shear and a constant convergence. Assuming an angular resolution of 0.01arcsec, we find that the angular positions of $\sim 10^8 \ms$ substructures at several kpc from the center of the macrolens halo can be directly measured if the size of the dust continuum emission region and the gradient of the surface brightness at the position of the perturber are sufficiently large. From the astrometric shift on a scale of a few times $10~$mas of an image perturbed by a subhalo with respect to an unperturbed macrolensed image, we can break the degeneracy between subhalo mass and distance provided that macrolensing parameters are determined from positions and fluxes of multiple images.

We present Subaru spectroscopy of the early-type lensing galaxies in the lens systems HST14113+5211 (z=0.464) and B2045+265 (z=0.868), being aimed at measuring the velocity dispersion of the lensing galaxies as an important component to their mass distributions and internal dynamics. For HST14113+5211 we have obtained 174+/-20km/s (1-sigma) inside an aperture of 0.9", and for B2045+265 we have obtained 213+/-23km/s inside an aperture of 1.2". To extract the significance of these information on the mass distributions and stellar dynamics of the lensing galaxies, we construct two alternative models for both the mass and velocity distributions. It is found that the mass-to-light ratios as derived from a joint lensing and stellar dynamics analysis are virtually consistent with those obtained from the use of the Fundamental Plane. Also, for HST14113+5211, we find that the total mass distribution is well reproduced by a power-law density profile with an index of 2, thereby suggesting that a singular isothermal model is a good fit to the lensing galaxy. In contrast, the corresponding slope for B2045+265 is flatter than isothermal, suggesting that additional contributions to lensing mass from surrounding structures (possibly the group of galaxies) provide the observed angular separations of multiple images.

Galactic stellar halos have long been considered to contain fossil information on early dynamical and chemical evolution of galaxies. We propose that the surface brightness distributions of old stellar halos contain the influence of reionization on early formation histories of galaxies. By assuming that reionization significantly suppresses star formation in small subgalactic clumps virialized after reionization redshift ($z_{\rm reion}$), we first numerically investigate how structural and kinematical properties of stellar halos formed from merging of these subgalactic clumps depend on $z_{\rm reion}$. We then discuss what observable properties of galactic stellar halos offer us the fossil records of reionization influence on hierarchical formation of halos based on the current results of numerical simulations. We particularly suggest that both the half-light radius of stellar halos and the slope of their surface brightness profile contain useful information on when star formation in subgalactic clumps were significantly influenced by reionization. By using the simulated surface brightness distributions of galactic stellar halos for models with different $z_{\rm reion}$, we also discuss how wide-field imaging studies of extragalactic halos will help us to elucidate the epoch of cosmic reionization.

We present mid-infrared imaging at 11.7 mu m for the quadruple lens systems, PG1115+080 and B1422+231, using the cooled mid-infrared camera and spectrometer (COMICS) attached on the Subaru telescope. These lensed QSOs are characterized by their anomalous optical and radio flux ratios as obtained for (A1, A2) images of PG1115+080 and (A, B, C) images of B1422+231, respectively, i.e., such flux ratios are hardly reproduced by lens models with smooth mass distribution. Our mid-infrared observations for these images have revealed that the mid-infrared flux ratio A2/A1 of PG1115+080 is virtually consistent with smooth lens models (but inconsistent with the optical flux ratio), whereas for B1422+231, the mid-infrared flux ratios among (A, B, C) are in good agreement with the radio flux ratios. We also identify a clear infrared bump in the spectral energy distributions of these QSOs, thereby indicating that the observed mid-infrared fluxes originate from a hot dust torus around a QSO nucleus. Based on the size estimate of the dust torus, we place limits on the mass of a substructure in these lens systems, causing the anomalous optical or radio flux ratios. For PG1115+080, the mass of a substructure inside an Einstein radius, M_E, is < 16 Msun, corresponding to either a star or a low-mass CDM subhalo having the mass of M_{100}^{SIS} < 2.2 * 10^4 Msun inside radius of 100 pc if modeled as a singular isothermal sphere (SIS). For B1422+231, we obtain M_E > 209 Msun, indicating that a CDM subhalo is more likely, having the mass of M_{100}^{SIS} > 7.4 * 10^4 Msun

We investigate the dynamical and chemical evolution of the Large Magellanic Cloud (LMC) interacting with the Galaxy and the Small Magellanic Cloud (SMC) based on a series of self-consistent chemodynamical simulations. Our numerical models are aimed at explaining the entire properties of the LMC, i.e., the observed structure and kinematics of its stellar halo and disk components as well as the populations of the field stars and star clusters. The main results of the present simulations are summarized as follows. (1) Tidal interaction between the Clouds and the Galaxy during the last 9 Gyr transforms the initially thin, non-barred LMC disk into the three different components: the central bar, thick disk, and kinematically hot stellar halo. The central bar is composed both of old field stars and newly formed ones with each fraction being equal in its innermost part. The final thick disk has the central velocity dispersion of $\sim$ 30 km s$^{-1}$ and shows rotationally supported kinematics with $V_{\rm m}/{\sigma}_{0}$ $\sim$ 2.3. (2) The stellar halo is formed during the interaction, consisting mainly of old stars originating from the outer part of the initially thin LMC disk. The outer halo shows velocity dispersion of $\sim$ 40 km s$^{-1}$ at the distance of 7.5 kpc from the LMC center and has somewhat inhomogeneous distribution of stars. The stellar halo contains relatively young, metal-rich stars with the mass fraction of 2 %.

Subaru observations of A1689 (z=0.183) are used to derive an accurate, model-independent mass profile for the entire cluster, r<2 Mpc/h, by combining magnification bias and distortion measurements. The projected mass profile steepens quickly with increasing radius, falling away to zero at r~1.0 Mpc/h, well short of the anticipated virial radius. Our profile accurately matches onto the inner profile, r<200 kpc/h, derived from deep HST/ACS images. The combined ACS and Subaru information is well fitted by an NFW profile with virial mass, (1.93 \pm 0.20)10^15 M_sun, and surprisingly high concentration, c_vir=13.7^{+1.4}_{-1.1}, significantly larger than theoretically expected (c_vir~4), corresponding to a relatively steep overall profile. A slightly better fit is achieved with a steep power-law model that has its 2D logarithmic slope -3 and core radius theta_c~1.7' (r_c~210 kpc/h), whereas an isothermal profile is strongly rejected. These results are based on a reliable sample of background galaxies selected to be redder than the cluster E/S0 sequence. By including the faint blue galaxy population a much smaller distortion signal is found, demonstrating that blue cluster members significantly dilute the true signal for r~400 kpc/h. This contamination is likely to affect most weak lensing results to date.

The Large Magellanic Cloud (LMC) has a unique cluster formation history in that nearly all of its globular clusters were formed either $\sim$ 13 Gyr ago or less than $\sim$ 3 Gyr ago. It is not clear what physical mechanism is responsible for the most recent cluster formation episode and thus the mysterious age gap between the LMC clusters. We first present results of gas dynamical N-body simulations of the evolution of the LMC in the context of its Galactic orbit and interactions with the SMC, paying special attention to the effect of tidal forces. We find that the first close encounter between the LMC and the Small Magellanic Cloud (SMC) about 4 Gyr ago was the beginning of a period of strong tidal interaction which likely induced dramatic gas cloud collisions, leading to an enhancement of the formation of globular clusters which has been sustained by strong tidal interactions to the present day. The tidal interaction results in the formation of a barred, elliptical, thick disk in the LMC. The model also predicts the presence of a large, diffuse stellar stream circling the Galaxy, which originated from the LMC.

The kinematic properties of tidal debris from an orbiting Galactic satellite is presented, on the assumption that its central part once contained the most massive Galactic globular cluster, omega Cen. We simulate dynamical evolution of a satellite galaxy that follows the present-day and likely past orbits of omega Cen and analyze the kinematic nature of their tidal debris and randomly generated Galactic stars comprising spheroidal halo and flat disk components. It is found that the debris stars show a retrograde rotation at about -100 km/s, which may accord with a recently discovered stellar stream at radial velocity of about 300 km/s towards the Galactic longitude of about 270 degree. These stars also contribute, only in part, to a reported retrograde motion of the outer halo in the North Galactic Pole direction, without significantly modifying local halo kinematics near the Sun. The prospects for future debris searches and the implications for the early evolution of the Galaxy are briefly presented.

We first demonstrate that globular cluster G1 could have been formed by tidal interaction between M31 and a nucleated dwarf galaxy (dE,N). Our fully self-consistent numerical simulations show that during tidal interaction between M31 and G1's progenitor dE,N with $M_{\rm B}$ $\sim$ -15 mag and its nucleus mass of $\sim$ $10^7$ $M_{\odot}$, the dark matter and the outer stellar envelope of the dE,N are nearly completely stripped whereas the nucleus can survive the tidal stripping because of its initially compact nature. The naked nucleus (i.e., G1) has orbital properties similar to those of its progenitor dE,N. The stripped stars form a metal-poor ([Fe/H] $\sim$ -1) stellar halo around M31 and its structure and kinematics depend strongly on the initial orbit of G1's progenitor dE,N. We suggest that the observed large projected distance of G1 from M31 ($\sim$ 40 kpc) can give some strong constraints on the central density of the dark matter halo of dE,N. We discuss these results in the context of substructures of M31's stellar halo recently revealed by Ferguson et al. (2002).

A significant fraction of non-baryonic or baryonic dark matter in galactic halos may consist of MASsive Compact Objects (MASCOs) with mass M=10^{1-4}M_{sun}. Possible candidates for such compact objects include primordial black holes or remnants of primordial (Population III) stars. We propose a method for directly detecting MASCOs in extragalactic halos, using the VLBI techniques with extremely high resolution. If a galactic halo comprising a large number of MASCOs produces multiple images of a background radio-loud QSO by gravitational lensing, then a high-resolution radio map of each macro-lensed image should reveal microlensing effects by MASCOs. To assess their observational feasibility, we simulate microlensing of the radio-loud, four-image lensed QSO, B1422+231, assuming angular resolution of ~0.01 mas. MASCOs are represented by point masses. For comparison, we also simulate microlensing of B1422+231 by singular isothermal spheres. We find that the surface brightness of the macro-lensed images shows distinct spatial patterns on the scale of the Einstein radius of the perturbers. In the case of point-mass perturbers, many tiny dark spots also appear in the macro-lensed images associated with a decrease in the surface brightness toward the fringe of the original QSO image, whereas no such spots are available in the SIS models. Based on the size, position and magnified or demagnified patterns of images, we shall be able to determine the mass and density profile of a MASCO as well as its spatial distribution and abundance in a galactic halo.

We present the kinematic properties of a tidally disrupted dwarf galaxy in the Milky Way, based on the hypothesis that its central part once contained the most massive Galactic globular cluster, omega Cen. Dynamical evolution of a self-gravitating progenitor galaxy that follows the present-day and likely past orbits of omega Cen is calculated numerically and the kinematic nature of their tidal debris is analyzed, combined with randomly generated stars comprising spheroidal halo and flat disk components. We show that the retrograde rotation of the debris stars at $\sim -100$ km/s accords with a recently discovered, large radial velocity stream at $\sim 300$ km/s towards the Galactic longitude of $\sim 270^\circ$. These stars also contribute, only in part, to a reported retrograde motion of the outer halo at the North Galactic Pole. The prospects for future debris searches and the implications for the early evolution of the Galaxy are briefly presented.

We present further observations of the extremely red object ERO J094258+4659.2, identified by \citet{iye00} as ERO R1 in their deep images of the cluster A851. We estimate its redshift independently by eight-band photometric redshift determination and cross-correlation of a new H-band spectrum with the optical spectra of local E/S0 galaxies, and conclude that it lies at $z \sim 1.5$. Although its colors are consistent both with an elliptical galaxy and an S0 galaxy at that redshift, its elongated shape and exponential luminosity profile suggest the presence of an evolved stellar disk component. We rule out the possibility that these properties are strongly influenced by gravitational lensing by the foreground cluster, and therefore conclude that this object is more likely to be an S0-like galaxy, rather than a lensed elliptical. The H-band spectrum does not show strong H$\alpha$ emission and the star formation rate therefore appears to be very modest. The presence of such a galaxy with an apparently relaxed disk of stars at this high redshift provides a new and strong constraint on theoretical models which aim to explain the formation and evolution of galaxies.

We set new limits on the mass of the Milky Way, making use of the latest kinematic information for Galactic satellites and halo objects. In order to bind these sample objects to the Galaxy, their rest-frame velocities must be lower than their escape velocities at their estimated distances. This constraint enables us to show that the mass estimate of the Galaxy is largely affected by several high-velocity objects (Leo I, Pal 3, Draco, and a few FHB stars), not by a single object alone (such as Leo I), as has often been the case in past analyses. We also find that a gravitational potential that gives rise to a declining rotation curve is insufficient to bind many of our sample objects to the Galaxy; a possible lower limit on the mass of the Galaxy is about 2.2 x 10^12 Msolar. To be more quantitative, we adopt a Bayesian likelihood approach to reproduce the observed distribution of the current positions and motions of the sample, in a prescribed Galactic potential that yields a flat rotation curve. This method enables a search for the most likely total mass of the Galaxy, without undue influence in the final result arising from the presence or absence of Leo I, provided that both radial velocities and proper motions are used. The most likely total mass derived from this method is 2.5^+0.5_-1.0 x 10^12 Msolar (including Leo I), and 1.8^+0.4_-0.7 x 10^12 Msolar (excluding Leo I).

We report medium-resolution (1-2 A) spectroscopy and broadband (UBV) photometry for a sample of 39 bright stars (the majority of which are likely to be giants) selected as metal-deficient candidates from an objective-prism survey concentrating on Galactic latitudes below |b| = 30 deg, the LSE survey of Drilling & Bergeron. Although the primary purpose of the LSE survey was to select OB stars (hence the concentration on low latitudes), the small number of bright metal-deficient giant candidates noted during this survey provide interesting information on the metal-weak thick disk (MWTD) population. The kinematics of the LSE giants indicate the presence of a rapidly rotating population, even at quite low metallicity. We consider the distribution of orbital eccentricity of the LSE giants as a function of [Fe/H], and conclude that the local fraction (i.e., within 1 kpc from the Sun) of metal-poor stars that might be associated with the MWTD is on the order of 30%-40% at abundances below [Fe/H] = -1.0. Contrary to recent analyses of previous (much larger) samples of non-kinematically selected metal-poor stars, we find that this relatively high fraction of local metal-poor stars associated with the MWTD may extend to metallicities below [Fe/H] = -1.6, much lower than had been considered before. We identify a subsample of 11 LSE stars that are very likely to be members of the MWTD, based on their derived kinematics; the lowest metallicity among these stars is [Fe/H] = -2.35. Implications of these results for the origin of the MWTD and for the formation of the Galaxy are considered. (abridged)

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.

The origin of metal-rich, highly flattened, and rapidly rotating disk globular cluster system in the Galaxy is one of longstanding issues in the context of the Galaxy formation. Our numerical simulations suggest a new ``two-fold'' scenario that the disk globular clusters are firstly formed in the high-pressure, dense central region of a gas-rich dwarf galaxy, as induced during the tidal interaction with the pre-existing, young thin disk of the Galaxy, and then dispersed into the disk region owing to the final tidal destruction of the merging dwarf. We also demonstrate that spatial distribution, total number, and metallicity distribution of the clusters formed in this minor merging depend on the mass ratio of the host to dwarf galaxy and the orbital configuration of merging. Based on these results, we discuss whether a minor merging event about 10 Gyr ago can explain both the Galactic thick disk and the disk globular clusters. Several other implications for the possible relation between the properties of disk galaxies and their disk globular clusters are also discussed.

We investigate the effects of numerous dark matter subhalos in a galaxy-sized halo on the events of strong lensing, to assess their presence as expected from the cold dark matter scenario. Lens galaxies are represented by a smooth ellipsoid in an external shear field and additional cold dark matter subhalos taken from Monte Carlo realizations which accord with recent N-body results. We also consider other possible perturbers, globular clusters and luminous dwarf satellites, for comparison. We then apply the models to the particular lens systems with four images, B1422+231 and PG1115+080, for which smooth lens models are unable to reproduce both the positions of the images and their radio flux ratios or dust-free optical flux ratios simultaneously. We show that the perturbations by both globular clusters and dwarf satellites are too small to change the flux ratios, whereas cold dark matter subhalos are most likely perturbers to reproduce the observed flux ratios in a statistically significant manner. This result suggests us the presence of numerous subhalos in lens galaxies, which is consistent with the results of cosmological N-body simulations.

We perform numerical simulations for the formation of the Galactic stellar halo, based on the currently favored cold dark matter (CDM) theory of galaxy formation. Our numerical models, taking into account both dynamical and chemical evolution processes in a consistent manner, are aimed at explaining observed structure and kinematics of the stellar halo in the context of hierarchical galaxy formation. The main results of the present simulations are summarized as follows. (1) Basic physical processes involved in the formation of the stellar halo, composed of metal-deficient stars with [Fe/H] $\le$ -1.0, are described by both dissipative and dissipationless merging of subgalactic clumps and their resultant tidal disruption in the course of gravitational contraction of the Galaxy at high redshift ($z$ $>$ 1). (2) The simulated halo has the density profile similar to the observed power-law form of $\rho (r)$ $\sim$ $r^{-3.5}$, and has also the similar metallicity distribution to the observations. The halo virtually shows no radial gradient for stellar ages and only small gradient for metallicities. (3) The dual nature of the halo, i.e., its inner flattened and outer spherical density distribution, is reproduced, at least qualitatively, by the present model. The outer spherical halo is formed via essentially dissipationless merging of small subgalactic clumps, whereas the inner flattened one is formed via three different mechanisms, i.e., dissipative merging between larger, more massive clumps, adiabatic contraction due to the growing Galactic disk, and gaseous accretion onto the equatorial plane.

We present a review of recent studies of the formation and evolution of the Milky Way galaxy, in particular based on large samples of non-kinematically selected stars with available proper motions. The Milky Way is argued to be a reasonable template for the formation of large spiral galaxies, the only one in which complete kinematical and abundance information can be readily obtained. Ongoing and future projects to obtain proper motions and spectral information for much larger samples of stars that will sharpen our perspective are discussed.

We consider the distortion in the cosmic microwave background (CMB) due to galactic winds at high redshift. Winds outflowing from galaxies have been hypothesized to be possible sources of metals in the intergalactic medium, which is known to have been enriched to $10^{-2.5}$ Z$_{\odot}$ at $z \sim 3$. We model these winds as functions of mass of the parent galaxy and redshift, assuming they activate at a common initial redshift, $z_{in}$, and calculate the mean $y$ distortion and the angular power spectrum of the distortion in the CMB. We find that the thermal Sunyaev-Zel'dovich (SZ) effect due to the winds to be consistent with previous estimates. The distortion due to the kinetic SZ (kSZ) effect is, however, found to be more important than the thermal SZ (tSZ) effect. We find that the distortion due to galactic winds is an important contribution to the power spectrum of distortion at very small angular scales $(l \sim 10^4)$. We also find that the power spectrum due to clustering dominates the Poisson power spectrum for $l \le 4\hbox{--}5 \times 10^5$. We show explicitly how the combined power spectrum from wind dominates over that of clusters at 217 GHz, relevant for PLANCK. We also show how these constraints change when the efficiency of the winds is varied.

We present a method, based on the kinematic analysis of the Galactic disk stars, to clarify whether the internal motions of the stellar system in spiral arms follow those expected in the density wave theory. The method relies on the comparison with the linear relation between the phases of spatial positions and epicyclic motions of stars, as drawn from the theory. The application of the method to the 78 Galactic Cepheids near the Sun, for which accurate proper motions are available from the {\it Hipparcos} Catalogue, has revealed that these Cepheids hold no correlation between both phases, thereby implying that their motions are in contradiction with the theoretical predictions. Possible reasons for this discrepancy are discussed and future prospects are outlined.

We develop a method for recovering the global density distribution of the ancient Galactic stellar halo prior to disk formation, based on the present orbits of metal-poor stars observed in the solar neighborhood. The method relies on the adiabatic invariance of the action integrals of motion for the halo population during the slow accumulation of a disk component, subsequent to earlier halo formation. The method is then applied to a sample of local stars with [Fe/H]<=-1.5, likely to be dominated by the halo component, taken from Beers et al.'s recently revised and supplemented catalog of metal-poor stars selected without kinematic bias. We find that even if the Galactic potential is made spherical by removing the disk component in an adiabatic manner, the halo density distribution in the inner halo region (R <= 15 kpc) remains moderately flattened, with axial ratio of about 0.8 for stars in the abundance range [Fe/H]<=-1.8 and about 0.7 for the more metal-rich interval -1.8<[Fe/H]<=-1.5. The outer halo remains spherical for both abundance intervals. We also find that this initial flattening of the inner halo is caused by the anisotropic velocity dispersions of the halo stars. These results suggest that the two-component nature of the present-day stellar halo, characterized by a highly flattened inner halo and nearly spherical outer halo, is a consequence of both an initially two-component density distribution of the halo (perhaps a signature of dissipative halo formation) and of the adiabatic flattening of the inner part by later disk formation. Further implications of our results for the formation of the Galaxy are also discussed, in particular in the context of the hierarchical clustering scenario of galaxy formation.

Motivated by the recently improved knowledge on the kinematic and chemical properties of the Galactic metal-poor stars, we present the numerical simulation for the formation of the Galactic stellar halo to interpret the observational results. As a model for the Galaxy contraction, we adopt the currently standard theory of galaxy formation based on the hierarchical assembly of the cold dark matter fluctuations. We find, for the simulated stars with [Fe/H]$\le-1.0$, that there is no strong correlation between metal abundances and orbital eccentricities, in good agreement with the observations. Moreover, the observed fraction of the low eccentricity stars is reproduced correctly for [Fe/H]$\le-1.6$ and approximately for the intermediate abundance range of $-1.6<$[Fe/H]$\le-1.0$. We show that this successful reproduction of the kinematics of the Galactic halo is a natural consequence of the hierarchical evolution of the subgalactic clumps seeded from the cold dark matter density fluctuations.

We present a revised catalog of 2106 Galactic stars, selected without kinematic bias, and with available radial velocities, distance estimates, and metal abundances in the range 0.0 <= [Fe/H] <= -4.0. This update of the Beers and Sommer-Larsen (1995) catalog includes newly-derived homogeneous photometric distance estimates, revised radial velocities for a number of stars with recently obtained high-resolution spectra, and refined metallicities for stars originally identified in the HK objective-prism survey (which account for nearly half of the catalog) based on a recent re-calibration. A subset of 1258 stars in this catalog have available proper motions, based on measurements obtained with the Hipparcos astrometry satellite, or taken from the updated Astrographic Catalogue (AC 2000; second epoch positions from either the Hubble Space Telescope Guide Star Catalog or the Tycho Catalogue), the Yale/San Juan Southern Proper Motion (SPM) Catalog 2.0, and the Lick Northern Proper Motion (NPM1) Catalog. Our present catalog includes 388 RR Lyrae variables (182 of which are newly added), 38 variables of other types, and 1680 non-variables, with distances in the range 0.1 to 40 kpc.

(Abbreviated) We present a detailed analysis of the space motions of 1203 solar-neighborhood stars with metal abundances [Fe/H] <= -0.6, on the basis of a recently revised and supplemented catalog of metal-poor stars selected without kinematic bias (Beers et al. 2000). This sample, having available proper motions, radial velocities, and distance estimates for stars with a wide range of metal abundances, is by far the largest such catalog to be assembled to date. Unlike essentially all previous kinematically selected catalogs, the metal-poor stars in our sample exhibit a diverse distribution of orbital eccentricities, e, with no apparent correlation between [Fe/H] and e. This demonstrates, clearly and convincingly, that the evidence offered by Eggen, Lynden-Bell, and Sandage (1962) for a rapid collapse of the Galaxy, an apparent correlation between the orbital eccentricity of halo stars with metallicity, is basically the result of their proper-motion selection bias. However, even in our non-kinematically selected sample, we have identified a small concentration of high-e stars at [Fe/H] = -1.7, which may originate, in part, from infalling gas during the early formation of the Galaxy. The implications of our results for the formation of the Galaxy are also discussed, in particular in the context of the currently favored CDM theory of hierarchical galaxy formation.

We present new limits on cosmological parameters from the statistics of gravitational lensing, based on the recently revised knowledge of the luminosity function and internal dynamics of E/S0 galaxies that are essential in lensing high-redshift QSOs. We find that the lens models using updated Schechter parameters for such galaxies, derived from the recent redshift surveys combined with morphological classification, are found to give smaller lensing probabilities than earlier calculated. Inconsistent adoption of these parameters from a mixture of various galaxy surveys gives rise to systematic biases in the results. We also show that less compact dwarf-type galaxies which largely dominate the faint part of the Schechter-form luminosity function contribute little to lensing probabilities, so that earlier lens models overestimate incidents of small separation lenses. Applications of the lens models to the existing lens surveys indicate that reproduction of both the lensing probability of optical sources and the image separations of optical and radio lenses is significantly improved in the revised lens models. The likelihood analyses allow us to conclude that a flat universe with Omega=0.3(+0.2-0.1) and Omega+Lambda=1 is most preferable, and a matter-dominated flat universe with Lambda=0 is ruled out at 98 % confidence level. These new limits are unaffected by inclusion of uncertainties in the lens properties.

The kinematics of 122 red giants and 124 RR Lyrae variables in the solar neighborhood is studied using accurate measurements of their proper motions by the Hipparcos astrometry satellite, combined with the published photometric distances, metal abundances and radial velocities. A majority of these sample stars have metal abundances with [Fe/H]<-1 and thus represent the old stellar populations in the Galaxy. The halo component with [Fe/H]<-1.6 is characterized by no systemic rotation and a radially elongated velocity ellipsoid. About 16% of such metal-poor stars have low orbital eccentricities e<0.4, and we see no evidence for the correlation between [Fe/H] and e. Based on the model for the e distribution of orbits, we show that this fraction of low e stars for [Fe/H]<-1.6 is explained from the halo component alone, without introducing the extra disk component claimed by recent workers. This is also supported by no significant change of the e distribution with the height from the Galactic plane. This metal-weak thick disk component appears to comprise only about 10% for -1.6< [Fe/H]<-1 and 20% for -1.4<[Fe/H]<-1. It is also verified that the disk has the mean rotation of about 195 km/s and the vertical extent of 1 kpc, which is consistent with the thick disk dominating at [Fe/H]=-0.6 to -1. The implications of these results for the early evolution of the Galaxy are also presented.

We present the three-dimensional orbital motions of metal-poor stars in conjunction with their metal abundances, for the purpose of getting insight into the formation process of the Galaxy. Our sample stars, which include metal-deficient red giants and RR Lyrae variables observed by the Hipparcos satellite, are least affected by known systematics, stemmed from kinematic bias, metallicity calibration, and secondary metal contamination of stellar surface. We find, for the stars in the metallicity range of [Fe/H]<-1, that there is no evidence for the correlation between [Fe/H] and their orbital eccentricities e. Even for [Fe/H]<-1.6, about 16% of the stars have e less than 0.4. We show that the e distribution of orbits for [Fe/H]<-1.6 is independent of the height |z| away from the Galactic plane, whereas for [Fe/H]>-1.6 the stars at |z|>1 kpc are systematically devoid of low-e orbits with e<0.6. This indicates that low-e stars with [Fe/H]<-1.6 belong to the halo component, whereas the rapidly-rotating thick disk with a scale height about 1 kpc has a metal-weak tail in the range of -1.6<[Fe/H]<-1. The fraction of this metal-weak thick disk appears to be only less than 20%. The significance of these results for the early evolution of the Galaxy is briefly discussed.

We investigate the hypothesis that $\lyal$ absorption lines arise in two populations of halos --- minihalos of small circular velocity ($V_c \la 55$ km s$^{-1}$) in which star formation and metal production are inhibited by photoionization of the UV background radiation, and large galactic halos ($55 \la V_c \la 250$ km s$^{-1}$) which possess stars and metals. Based on the model of $\lyal$-absorbing gas confined in both populations of halos, we attempt to explain the recent observations of (1) associations of visible galaxies with $\lyal$ lines at low redshifts $z \la 1$, and (2) metal lines associated with a non-negligible fraction of low H I column density $\lyal$ lines at $z \sim 3$. For galactic halos, we find that photoionized gas clouds confined in the pressure of ambient hot gas can produce $\lyal$ absorptions with H I column density as low as $10^{14}$ cm$^{-2}$, and that the impact parameter of a sightline for such absorptions matches well with the observed radius of gaseous envelope in a typical luminous galaxy. Using the Press- Schechter prescription for the mass function of halos, we also show that the fraction of $\lyal$ lines with associated metal lines can be understood in terms of the fraction of $\lyal$ absorbers that are associated with galactic halos. In particular, the reported fraction of $\sim 0.5--0.75$ at $z \sim 3$ is reproduced when the boundary value of $V_c$ to separate mini or galactic halos is $40 \sim 60$ km s$^{-1}$, which is consistent with the theoretical prediction of galaxy formation under photoionization. .....

Following the recent suggestion that it is dwarf galaxies in clusters -- as opposed to large ellipticals -- that provide the intracluster gas, we estimate the metallicity of the intracluster medium (ICM) in such a case. We derive analytical expressions for the fraction of mass of dwarf galaxies that is ejected, and estimate the metallicity of the resulting intracluster gas. We find that the metallicity resulting from this hypothesis is adequate only for clusters with low-metallicity gas. Since galactic winds from dwarf galaxies occur much earlier than those from ellipticals, due to the smaller binding energy of the gas, we show that the gas ejected by dwarf galaxies is enriched mostly by Type II supernovae, as opposed to Type I supernovae in the case of large galaxies. We further point out that the gas in small scale structures, which never cools and forms stars due to low temperatures and, consequently, large cooling time scale, gets dispersed in the process of hierarchical clustering, and is incorporated as the diffuse intracluster gas. We show that this process can provide enough hot gas to be compatible with X-ray observations in rich clusters.

The evolution of large-scale magnetic fields in disk galaxies is investigated numerically. The gasdynamical simulations in a disk perturbed by spiral or bar potential are incorporated into the kinematic calculations of induction equations to elucidate the effects of non-axisymmetric disk structure on magnetic fields. The effects of interstellar turbulence are given as the turbulent diffusion of magnetic fields. The usually adopted dynamo mechanism of alpha-effect is not considered in our computations, because it is not obvious about the actual existence of the effect in a galaxy. Our principal concern is to clear how observationally and theoretically well-established gas flow affects the magnetic-field structure and evolution, without putting a lot of artificial parameters in the model. We have found that the density-wave streaming motion of gas has a significant influence on the distribution of magnetic fields: the lines of force are well aligned with spiral arms due to the compressional and additional shearing flow of gas in these regions.

We show that the relatively strong magnetic fields ($\ge 1 \mu$G) in high redshift objects can be explained by the combined action of an evolving protogalactic fluctuation and electrodynamic processes providing the magnetic seed fields. Three different seed field mechanisms are reviewed and incorporated into a spherical "top-hat" model and tidal torque theory for the fate of a forming galaxy in an expanding universe. Very weak fields $10^{-19} \sim 10^{-23}$G created in an expanding over-dense region are strongly enhanced due to the dissipative disk formation by a factor $\sim 10^4$, and subsequently amplified by strong non-axisymmetric flow by a factor $\sim 10^{6-10}$, depending on the cosmological parameters and the epoch of galaxy formation. The resulting field strength at $z \sim 0.395$ can be of the order of a few $\mu$G and be close to this value at $z \sim 2$.

We propose a method of distance determination based on the internal structure and dynamics of disk galaxies. The method relies on the universal luminosity profile of a stellar disk represented by an exponential law. Calibrating nearby galaxies with known distances, it is found that the scale length of the disk is tightly correlated with the specific combination of central surface brightness {\it and} rotational velocity at a characteristic radius of 2.15 scale lengths from the center. This suggests that the scale length of the disk may be used as an indicator for extragalactic distance scale. The application of this relation to M51 and M100 allows us to arrive at the distances of about 6 Mpc and 14 Mpc, respectively, implying a Hubble constant of $H_0 = 92 \sim 94$ km s$^{-1}$ Mpc$^{-1}$.

We investigate the effects of the observed UV background radiation on galaxy formation. Photoionization by UV radiation decreases the cooling rate of the gas in halos, so that objects with only large density contrasts can self-shield against the background radiation, thereby allowing the shielded neutral cores to cool and form stars. In the context of the CDM model, we use the criterion that self-shielding is essential for star formation to calculate the mass function of galaxies based on both Press-Schechter and the peaks formalisms. The ionizing UV radiation results in the inhibition of galaxy formation --- we show the decrease in the number density of the galaxies quantitatively. We also find that the merging in general is made inefficient by the UV photons through photoionization of the gas in the bigger system into which small objects are incorporated. The latter means that, in a merging dominated region, where the number density at the low-mass end ($M_b \la 10^{10}M_{\odot}$) is usually expected to decrease with time, the trend is {\it reversed} (the number of low-mass galaxies is more at lower redshifts $z$) due to the decreasing UV flux with time after $z \sim 2$. We further discuss the implication of our results for the number counts of galaxies and possible evolution of the luminosity function of the galaxies.