Galaxies change their properties as they assemble into clusters. In order to understand the physics behind that, we need to go back in time and observe directly what is occurring in galaxies as they fall into a cluster. We have conducted a narrow-band and $J$-band imaging survey on a cluster CL1604-D at $z=0.923$ using a new infrared instrument SWIMS installed at the Subaru Telescope. The narrow-band filter, NB1261, matches to H$\alpha$ emission from the cluster at $z=0.923$. Combined with a wide range of existing data from various surveys, we have investigated galaxy properties in and around this cluster in great detail. We have identified 27 H$\alpha$ emitters associated with the cluster. They have significant overlap with MIPS 24$\mu$m sources and are located exclusively in the star forming regime on the rest-frame $UVJ$ diagram. We have identified two groups of galaxies near the cluster in the 2D spatial distribution and the phase-space diagram, which are likely to be in-falling to the cluster main body. We have compared various physical properties of star forming galaxies, such as specific star formation rates (burstiness) and morphologies (merger) as a function of environment; cluster center, older group, younger group, and the field. As a result, a global picture has emerged on how the galaxy properties are altered as they assemble into a denser region. This includes the occurrence of mergers, enhancement of star formation activity, excursion to the dusty starburst phase, and eventual quenching to a passive phase.

We present the results of a dust-reverberation survey of quasars at redshifts z<0.6. We found a delayed response of the K-band flux variation after the optical flux variation in 25 out of 31 targets, and obtained the lag time between them for 22 targets. Combined with the results for nearby Seyfert galaxies, we provide the largest homogeneous collection of K-band dust-reverberation data for 36 type 1 active galactic nuclei (AGNs). This doubles the sample and includes the most distant AGN and the largest lag so far measured. We estimated the optical luminosity of the AGN component of each target using three different methods: spectral decomposition, the flux-variation-gradient method, and image decomposition. We found a strong correlation between the reverberation radius for the innermost dust torus and the optical luminosity over a range of approximately four orders of magnitude in luminosity, as is already known for Seyfert galaxies. We estimated the luminosity distances of the AGNs based on their dust-reverberation lags, and found that the data in the redshift-distance diagram are consistent with the current standard estimates of the cosmological parameters. We also present the radius-luminosity relations for isotropic luminosity indicators such as the hard X-ray (14--195 keV), [OIV] 25.89 um, and mid-infrared (12 um) continuum luminosities, which are applicable to obscured AGNs.

We present properties of two types of bulges (classical- and pseudo- bulges) in 20 luminous infrared galaxies (LIRGs) observed in the near infrared of the $H$, $K_s$ and 1.91$\mu$m narrow-band targeting at the hydrogen Pa$\alpha$ emission line by the University of Tokyo Atacama Observatory (TAO) 1.0 m telescope. To classify the two types of bulges, we first perform a two-dimensional bulge-disk decomposition analysis in the $K_\mathrm{s}$-band images. The result shows a tentative bimodal distribution of S\'ersic indices with a separation at $\log(n_b)\sim0.5$, which is consistent with that of classical and normal galaxies. We next measure extents of the distribution of star forming regions in Pa$\alpha$ emission line images, normalized with the size of the bulges, and find that they decrease with increasing S\'ersic indices. These results suggest that star-forming galaxies with classical bulges have compact star forming regions concentrated within the bulges, while those with pseudobulges have extended star forming regions beyond the bulges, suggesting that there are different formation scenarios at work in classical and pseudobulges.

The formal division of the distance indicators into primary and secondary leads to difficulties in description of methods which can actually be used in two ways: with, and without the support of the other methods for scaling. Thus instead of concentrating on the scaling requirement we concentrate on all methods of distance determination to extragalactic sources which are designated, at least formally, to use for individual sources. Among those, the Supernovae Ia is clearly the leader due to its enormous success in determination of the expansion rate of the Universe. However, new methods are rapidly developing, and there is also a progress in more traditional methods. We give a general overview of the methods but we mostly concentrate on the most recent developments in each field, and future expectations.

In Yoshii et al. (2014), we described a new method for measuring extragalactic distances based on dust reverberation in active galactic nuclei (AGNs), and we validated our new method with Cepheid variable stars. In this paper, we validate our new method with Type Ia supernovae (SNe Ia) which occurred in two of the AGN host galaxies during our AGN monitoring program: SN 2004bd in NGC 3786 and SN 2008ec in NGC 7469. Their multicolor light curves were observed and analyzed using two widely accepted methods for measuring SN distances, and the distance moduli derived are $\mu=33.47\pm 0.15$ for SN 2004bd and $33.83\pm 0.07$ for SN 2008ec. These results are used to obtain independently the distance measurement calibration factor, $g$. The $g$ value obtained from the SN Ia discussed in this paper is $g_{\rm SN} = 10.61\pm 0.50$ which matches, within the range of 1$\sigma$ uncertainty, $g_{\rm DUST} = 10.60$, previously calculated ab initio in Yoshii et al. (2014). Having validated our new method for measuring extragalactic distances, we use our new method to calibrate reverberation distances derived from variations of H$\beta$ emission in the AGN broad line region (BLR), extending the Hubble diagram to $z\approx 0.3$ where distinguishing between cosmologies is becoming possible.

We present an analysis of Mg II $\lambda2798$ and Fe II UV emission lines for archival Sloan Digital Sky Survey (SDSS) quasars to explore diagnostics of the magnesium-to-iron abundance ratio in a broad-line region cloud. Our sample consists of 17,432 quasars selected from the SDSS Data Release 7 with a redshift range of $0.72 < z < 1.63$. A strong anticorrelation between Mg II equivalent width (EW) and the Eddington ratio is found, while only a weak positive correlation is found between Fe II EW and the Eddington ratio. To investigate the origin of these differing behaviors of Mg II and Fe II emission lines, we have performed photoionization calculations using the Cloudy code, where constraints from recent reverberation mapping studies are considered. We find from calculations that (i) Mg II and Fe II emission lines are created at different regions in a photoionized cloud, and (ii) their EW correlations with the Eddington ratio can be explained by just changing the cloud gas density. These results indicate that the Mg II/Fe II flux ratio, which has been used as a first-order proxy for the Mg/Fe abundance ratio in chemical evolution studies with quasar emission lines, depends largely on the cloud gas density. By correcting this density dependence, we propose new diagnostics of the Mg/Fe abundance ratio for a broad line region cloud. Comparing the derived Mg/Fe abundance ratios with chemical evolution models, we suggest that $\alpha$-enrichment by mass loss from metal-poor intermediate-mass stars occurred at $z\sim2$ or earlier.

We have developed a near-infrared camera called ANIR (Atacama Near-InfraRed camera) for the University of Tokyo Atacama Observatory 1.0m telescope (miniTAO) installed at the summit of Cerro Chajnantor (5640 m above sea level) in northern Chile. The camera provides a field of view of 5'.1 $\times$ 5'.1 with a spatial resolution of 0".298 /pixel in the wavelength range of 0.95 to 2.4 $\mu$m. Taking advantage of the dry site, the camera is capable of hydrogen Paschen-$\alpha$ (Pa$\alpha$, $\lambda=$1.8751 $\mu$m in air) narrow-band imaging observations, at which wavelength ground-based observations have been quite difficult due to deep atmospheric absorption mainly from water vapor. We have been successfully obtaining Pa$\alpha$ images of Galactic objects and nearby galaxies since the first-light observation in 2009 with ANIR. The throughputs at the narrow-band filters ($N1875$, $N191$) including the atmospheric absorption show larger dispersion (~10%) than those at broad-band filters (a few %), indicating that they are affected by temporal fluctuations in Precipitable Water Vapor (PWV) above the site. We evaluate the PWV content via the atmospheric transmittance at the narrow-band filters, and derive that the median and the dispersion of the distribution of the PWV are 0.40+/-0.30 mm for $N1875$ and 0.37+/-0.21 mm for $N191$, which are remarkably smaller (49+/-38% for $N1875$ and 59+/-26% for $N191$) than radiometry measurements at the base of Cerro Chajnantor (5100 m alt.). The decrease in PWV can be explained by the altitude of the site when we assume that the vertical distribution of the water vapor is approximated at an exponential profile with scale heights within 0.3-1.9 km (previously observed values at night). We thus conclude that miniTAO/ANIR at the summit of Cerro Chajnantor indeed provides us an excellent capability for a "ground-based" Pa$\alpha$ observation.

Luminous infrared galaxies (LIRGs) are enshrouded by a large amount of dust, produced by their active star formation, and it is difficult to measure their activity in the optical wavelength. We have carried out Pa$\alpha$ narrow-band imaging observations of 38 nearby star-forming galaxies including 33 LIRGs listed in $IRAS$ RBGS catalog with the Atacama Near InfraRed camera (ANIR) on the University of Tokyo Atacama Observatory (TAO) 1.0 m telescope (miniTAO). Star formation rates (SFRs) estimated from the Pa$\alpha$ fluxes, corrected for dust extinction using the Balmer Decrement Method (typically $A_V$ $\sim$ 4.3 mag), show a good correlation with those from the bolometric infrared luminosity of $IRAS$ data within a scatter of 0.27 dex. This suggests that the correction of dust extinction for Pa$\alpha$ flux is sufficient in our sample. We measure the physical sizes and the surface density of infrared luminosities ($\Sigma_{L(\mathrm{IR})}$) and $SFR$ ($\Sigma_{SFR}$) of star-forming region for individual galaxies, and find that most of the galaxies follow a sequence of local ultra luminous or luminous infrared galaxies (U/LIRGs) on the $L(\mathrm{IR})$-$\Sigma_{L(\mathrm{IR})}$ and $SFR$-$\Sigma_{SFR}$ plane. We confirm that a transition of the sequence from normal galaxies to U/LIRGs is seen at $L(\mathrm{IR})=8\times10^{10}$ $L_{\odot}$. Also, we find that there is a large scatter in physical size, different from those of normal galaxies or ULIRGs. Considering the fact that most of U/LIRGs are merging or interacting galaxies, this scatter may be caused by strong external factors or differences of their merging stage.

A large reverberation mapping study of the Seyfert 1 galaxy NGC 7469 has yielded emission-line lags for Hbeta 4861 and He II 4686 and a central black hole mass measurement of about 10 million solar masses, consistent with previous measurements. A very low level of variability during the monitoring campaign precluded meeting our original goal of recovering velocity-delay maps from the data, but with the new Hbeta measurement, NGC 7469 is no longer an outlier in the relationship between the size of the Hbeta-emitting broad-line region and the AGN luminosity. It was necessary to detrend the continuum and Hbeta and He II 4686 line light curves and those from archival UV data for different time-series analysis methods to yield consistent results.

We present the results of a dust reverberation survey for 17 nearby Seyfert 1 galaxies, which provides the largest homogeneous data collection for the radius of the innermost dust torus. A delayed response of the K-band light curve after the V-band light curve was found for all targets, and 49 measurements of lag times between the flux variation of the dust emission in the K band and that of the optical continuum emission in the V band were obtained. The lag times strongly correlated with the optical luminosity in the luminosity range of M_V=-16 to -22 mag, and the regression analysis was performed to obtain the correlation log $\Delta t$ (days) = -2.11 -0.2 M_V assuming $\Delta t \propto L^{0.5}$, which was theoretically expected. We discuss the possible origins of the intrinsic scatter of the dust lag-luminosity correlation, which was estimated to be about 0.13 dex, and we find that the difference of internal extinction and delayed response of changes in lag times to the flux variations could have partly contributed to intrinsic scatter. However, we could not detect any systematic change of the correlation with the subclass of the Seyfert type or the Eddington ratio. Finally, we compare the dust reverberation radius with the near-infrared interferometric radius of the dust torus and the reverberation radius of broad Balmer emission lines. The interferometric radius in the K band was found to be systematically larger than the dust reverberation radius in the same band by about a factor of two, which could be interpreted by the difference between the flux-weighted radius and the response-weighted radius of the innermost dust torus. The reverberation radius of the broad Balmer emission lines was found to be systematically smaller than the dust reverberation radius by about a factor of 4-5, which strongly supports the unified scheme of the Seyfert type of active galactic nuclei. (Abridged)

We have pioneered a new method for the measurement of extragalactic distances. This method uses the time-lag between variations in the short wavelength and long wavelength light from an active galactic nucleus (AGN), based on a quantitative physical model of dust reverberation that relates the time-lag to the absolute luminosity of the AGN. We use the large homogeneous data set from intensive monitoring observations in optical and near-infrared wavelength bands with the dedicated 2-m MAGNUM telescope to obtain the distances to 17 AGNs in the redshift range z=0.0024 to z=0.0353. These distance measurements are compared with distances measured using Cepheid variable stars, and are used to infer that H_0= 73 +- 3 (random) km/s/Mpc. The systematic error in H_0 is examined, and the uncertainty in the size distribution of dust grains is the largest source of the systematic error, which is much reduced for a sample of AGNs for which their parameter values in the model of dust reverberation are individually measured. This AGN time-lag method can be used beyond 30 Mpc, the farthest distance reached by extragalactic Cepheids, and can be extended to high-redshift quasi-stellar objects.

The entire evolution of the Milky Way, including its mass-assembly and star-formation history, is imprinted onto the chemo-dynamical distribution function of its member stars, f(x, v, [X/H]), in the multi-dimensional phase space spanned by position, velocity, and elemental abundance ratios. In particular, the chemo-dynamical distribution functions for low-mass stars (e.g., G- or K-type dwarfs) are precious tracers of the earliest stages of the Milky Way's formation, since their main-sequence lifetimes approach or exceed the age of the universe. A basic tenet of essentially all previous analyses is that the stellar metallicity, usually parametrized as [Fe/H], is conserved over time for main-sequence stars (at least those that have not been polluted due to mass transfer from binary companions). If this holds true, any correlations between metallicity and kinematics for long-lived main-sequence stars of different masses, effective temperatures, or spectral types must strictly be the same, since they reflect the same mass-assembly and star-formation histories. By analyzing a sample of nearby metal-poor halo and thick-disk stars on the main sequence, taken from Data Release 8 of the Sloan Digital Sky Survey, we find that the median metallicity of G-type dwarfs is systematically higher (by about 0.2 dex) than that of K-type dwarfs having the same median rotational velocity about the Galactic center. If it can be confirmed, this finding may invalidate the long-accepted assumption that the atmospheric metallicities of long-lived stars are conserved over time.

We present new data for five under-luminous type II-plateau supernovae (SNe IIP), namely SN 1999gn, SN 2002gd, SN 2003Z, SN 2004eg and SN 2006ov. This new sample of low-luminosity SNe IIP (LL SNe IIP) is analyzed together with similar objects studied in the past. All of them show a flat light curve plateau lasting about 100 days, an under luminous late-time exponential tail, intrinsic colours that are unusually red, and spectra showing prominent and narrow P-Cygni lines. A velocity of the ejected material below 10^3 km/s is inferred from measurements at the end of the plateau. The 56Ni masses ejected in the explosion are very small (less than 10^-2 solar masses). We investigate the correlations among 56Ni mass, expansion velocity of the ejecta and absolute magnitude in the middle of the plateau, confirming the main findings of Hamuy (2003), according to which events showing brighter plateau and larger expansion velocities are expected to produce more 56Ni. We propose that these faint objects represent the low luminosity tail of a continuous distribution in parameters space of SNe IIP. The physical properties of the progenitors at the explosion are estimated through the hydrodynamical modeling of the observables for two representative events of this class, namely SN 2005cs and SN 2008in. We find that the majority of LL SNe IIP, and quite possibly all, originate in the core-collapse of intermediate mass stars, in the mass range 10-15 solar masses.

The near-Earth asteroid (308635) 2005 YU55 is a potentially hazardous asteroid which was discovered in 2005 and passed Earth on November 8th 2011 at 0.85 lunar distances. This was the closest known approach by an asteroid of several hundred metre diameter since 1976 when a similar size object passed at 0.5 lunar distances. We observed 2005 YU55 from ground with a recently developed mid-IR camera (miniTAO/MAX38) in N- and Q-band and with the Submillimeter Array (SMA) at 1.3 mm. In addition, we obtained space observations with Herschel/PACS at 70, 100, and 160 micron. Our thermal measurements cover a wide range of wavelengths from 8.9 micron to 1.3 mm and were taken after opposition at phase angles between -97 deg and -18 deg. We performed a radiometric analysis via a thermophysical model and combined our derived properties with results from radar, adaptive optics, lightcurve observations, speckle and auxiliary thermal data. We find that (308635) 2005 YU55 has an almost spherical shape with an effective diameter of 300 to 312 m and a geometric albedo pV of 0.055 to 0.075. Its spin-axis is oriented towards celestial directions (lam_ecl, beta_ecl) = (60 deg +/- 30deg, -60 deg +/- 15 deg), which means it has a retrograde sense of rotation. The analysis of all available data combined revealed a discrepancy with the radar-derived size. Our radiometric analysis of the thermal data together with the problem to find a unique rotation period might be connected to a non-principal axis rotation. A low to intermediate level of surface roughness (r.m.s. of surface slopes in the range 0.1 - 0.3) is required to explain the available thermal measurements. We found a thermal inertia in the range 350-800 Jm^-2s^-0.5K^-1, very similar to the rubble-pile asteroid (25143) Itokawa and indicating a mixture of low conductivity fine regolith with larger rocks and boulders of high thermal inertia on the surface.

Flat space models with multiply connected topology, which have compact dimensions, are tested against the distribution of high-redshift ($z \geq 4$) quasars of the Sloan Digital Sky Survey (SDSS). When the compact dimensions are smaller in size than the observed universe, topological lensing occurs, in which multiple images of single objects (ghost images) are observed. We improve on the recently introduced method to identify ghost images by means of four-point statistics. Our method is valid for any of the 17 multiply connected flat models, including nontoroial ones that are compactified by screw motions or glide reflection. Applying the method to the data revealed one possible case of topological lensing caused by sixth-turn screw motion, however, it is consistent with the simply connected model by this test alone. Moreover, simulations suggest that we cannot exclude the other space models despite the absence of their signatures. This uncertainty mainly originates from the patchy coverage of SDSS in the South Galactic cap, and the situation will be improved by future wide-field spectroscopic surveys.

High resolution (0.4 arcsec) Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 0 observations of HCO+(4-3) and HCN(4-3) toward a mid-stage infrared bright merger VV114 have revealed compact nuclear (<200 pc) and extended (3 - 4 kpc) dense gas distribution across the eastern part of the galaxy pair. We find a significant enhancement of HCN(4-3) emission in an unresolved compact and broad (290km/s) component found in the eastern nucleus of VV114, and we suggest dense gas associated with the surrounding material around an Active Galactic Nucleus (AGN), with a mass upper limit of < 4 x 10^8 Msun. The extended dense gas is distributed along a filamentary structure with resolved dense gas concentrations (230pc; 10^6 Msun) separated by a mean projected distance of 600 pc, many of which are generally consistent with the location of star formation traced in Pa alpha emission. Radiative transfer calculations suggest moderately dense (10^5 - 10^6 cm^-3) gas averaged over the entire emission region. These new ALMA observations demonstrate the strength of the dense gas tracers in identifying both the AGN and star formation activity in a galaxy merger, even in the most dust enshrouded environments in the local universe.

The orbital motions of halo stars in the Milky Way reflect the orbital motions of the progenitor systems in which they formed, making it possible to trace the mass-assembly history of the Galaxy. Direct measurement of three-dimensional velocities, based on accurate proper motions and line-of-sight velocities, has revealed that the majority of halo stars in the inner-halo region move on eccentric orbits. However, our understanding of the motions of distant, in-situ halo-star samples is still limited, due to the lack of accurate proper motions for these stars. Here we explore a model-independent analysis of the line-of-sight velocities and spatial distribution of a recent sample of 1865 carefully selected halo blue horizontal-branch (BHB) stars within 30 kpc of the Galactic center. We find that the mean rotational velocity of the very metal-poor ([Fe/H] < -2.0) BHB stars significantly lags behind that of the relatively more metal-rich ([Fe/H] > -2.0) BHB stars. We also find that the relatively more metal-rich BHB stars are dominated by stars with eccentric orbits, as previously observed for other stellar samples in the inner-halo region. By contrast, the very metal-poor BHB stars are dominated by stars on rounder, lower-eccentricity orbits. Our results indicate that the motion of the progenitor systems of the Milky Way that contributed to the stellar populations found within 30 kpc correlates directly with their metal abundance, which may be related to their physical properties such as gas fractions. These results are consistent with the existence of an inner/outer halo structure for the halo system, as advocated by Carollo et al. (2010).

In order to test a recent hypothesis that the dispersion in the Schmidt-Kennicutt law arises from variations in the evolutionary stage of star forming molecular clouds, we compared molecular gas and recent star formation in an early-phase merger galaxy pair, Taffy I (UGC\ 12915/UGC\ 12914, VV\ 254) which went through a direct collision 20 Myr ago and whose star forming regions are expected to have similar ages. Narrow-band Pa-alpha image is obtained using the ANIR near-infrared camera on the mini-TAO 1m telescope. The image enables us to derive accurate star formation rates within the galaxy directly. The total star formation rate, 22.2 M_sun/yr, was found to be much higher than previous estimates. Ages of individual star forming blobs estimated from equivalent widths indicate that most star forming regions are ~7 Myr old, except for a giant HII region at the bridge which is much younger. Comparison between star formation rates and molecular gas masses for the regions with the same age exhibits a surprisingly tight correlation, a slope of unity, and star formation efficiencies comparable to those of starburst galaxies. These results suggest that Taffy I has just evolved into a starburst system after the collision, and the star forming sites are at a similar stage in their evolution from natal molecular clouds except for the bridge region. The tight Schmidt-Kennicutt law supports the scenario that dispersion in the star formation law is in large part due to differences in evolutionary stage of star forming regions.

We present theoretical calculations for the differential distribution of stellar orbital eccentricity for a sample of solar-neighbour halo stars. Two types of static, spherical gravitational potentials are adopted to define the eccentricity e for given energy E and angular momentum L, such as an isochrone potential and a Navarro-Frenk-White potential that can serve as two extreme ends covering in-between any realistic potential of the Milky Way halo. The solar-neighbour eccentricity distribution \Delta N(e) is then formulated, based on a static distribution function of the form f(E,L) in which the velocity anisotropy parameter \beta monotonically increases in the radial direction away from the galaxy center, such that beta is below unity (near isotropic velocity dispersion) in the central region and asymptotically approaches \sim 1 (radially anisotropic velocity dispersion) in the far distant region of the halo. We find that \Delta N(e) sensitively depends upon the radial profile of \beta, and this sensitivity is used to constrain such profile in comparison with some observational properties of \Delta N_{obs}(e) recently reported by Carollo et al. (2010). Especially, the linear e-distribution and the fraction of higher-e stars for their sample of solar-neighbour inner-halo stars rule out a constant profile of \beta, contrary to the opposite claim by Bond et al. (2010). Our constraint of \beta \lesssim 0.5 at the galaxy center indicates that the violent relaxation that has acted on the inner halo is effective within a scale radius of \sim 10 kpc from the galaxy center. We discuss that our result would help understand the formation and evolution of the Milky Way halo.

We present a method to constrain cosmic topology from the distribution of astronomical objects projected on the celestial sphere. This is an extension of the 3D method introduced in Fujii & Yoshii (2011) that is to search for a pair of pairs of observed objects (quadruplet) linked by a holonomy, i.e., the method we present here is to search for a pair of celestial sphere $n$-tuplets for $n \geq 3$. We find, however, that this method is impractical to apply in realistic situations due to the small signal to noise ratio. We conclude therefore that it is unrealistic to constrain the topology of the Universe from the celestial distribution, and the 3D catalogs are necessary for the purpose.

A new, improved version of a cosmic crystallography method for constraining cosmic topology is introduced. Like the circles-in-the-sky method using CMB data, we work in a thin, shell-like region containing plenty of objects. Two pairs of objects (quadruplet) linked by a holonomy show a specific distribution pattern, and three filters of \emph{separation, vectorial condition}, and \emph{lifetime of objects} extract these quadruplets. Each object $P_i$ is assigned an integer $s_i$, which is the number of candidate quadruplets including $P_i$ as their members. Then an additional device of $s_i$-histogram is used to extract topological ghosts, which tend to have high values of $s_i$. In this paper we consider flat spaces with Euclidean geometry, and the filters are designed to constrain their holonomies. As the second filter, we prepared five types that are specialized for constraining specific holonomies: one for translation, one for half-turn corkscrew motion and glide reflection, and three for $n$-th turn corkscrew motion for $n=4, 3,$ and 6. {Every multiconnected space has holonomies that are detected by at least one of these five filters.} Our method is applied to the catalogs of toy quasars in flat $\Lambda$-CDM universes whose typical sizes correspond to $z\sim 5$. With these simulations our method is found to work quite well. {These are the situations in which type-II pair crystallography methods are insensitive because of the tiny number of ghosts. Moreover, in the flat cases, our method should be more sensitive than the type-I pair (or, in general, $n$-tuplet) methods because of its multifilter construction and its independence from $n$.}

We examine whether the spectral energy distribution of UV continuum emission of active galactic nuclei changes during flux variation. We used multi-epoch photometric data of QSOs in the Stripe 82 observed by the SDSS Legacy Survey and selected 10 bright QSOs observed with high photometric accuracies, in the redshift range of z = 1.0-2.4 where strong broad emission lines such as Ly\alpha and CIV do not contaminate SDSS filters, to examine spectral variation of the UV continuum emission with broadband photometries. All target QSOs showed clear flux variations during the monitoring period 1998-2007, and the multi-epoch flux data in two different bands obtained on the same night showed a linear flux-to-flux relationship for all target QSOs. Assigning the flux in the longer wavelength to the x-axis in the flux-to-flux diagram, the x-intercept of the best-fit linear regression line was positive for most targets, which means that their colors in the observing bands become bluer as they become brighter. Then, the host-galaxy flux was estimated on the basis of the correlation between the stellar mass of the bulge of the host galaxy and the central black hole mass. We found that the longer-wavelength flux of the host galaxy was systematically smaller than that of the fainter extension of the best-fit regression line at the same shorter-wavelength flux for most targets. This result strongly indicates that the spectral shape of the continuum emission of QSOs in the UV region usually becomes bluer as it becomes brighter. We found that the multi-epoch flux-to-flux plots could be fitted well with the standard accretion disk model changing the mass accretion rate with a constant black hole mass for most targets. This finding strongly supports the standard accretion disk model for UV continuum emission of QSOs.

We present the final results from a high sampling rate, multi-month, spectrophotometric reverberation mapping campaign undertaken to obtain either new or improved Hbeta reverberation lag measurements for several relatively low-luminosity AGNs. We have reliably measured thetime delay between variations in the continuum and Hbeta emission line in six local Seyfert 1 galaxies. These measurements are used to calculate the mass of the supermassive black hole at the center of each of these AGNs. We place our results in context to the most current calibration of the broad-line region (BLR) R-L relationship, where our results remove outliers and reduce the scatter at the low-luminosity end of this relationship. We also present velocity-resolved Hbeta time delay measurements for our complete sample, though the clearest velocity-resolved kinematic signatures have already been published.

We present theoretical calculations for the differential distribution of stellar orbital eccentricity in a galaxy halo, assuming that the stars constitute a spherical, collisionless system in dynamical equilibrium with a dark matter halo. In order to define the eccentricity e of a halo star for given energy E and angular momentum L, we adopt two types of gravitational potential, such as an isochrone potential and a Navarro-Frenk-White potential, that could form two ends covering in-between any realistic potential of dark matter halo. Based on a distribution function of the form f(E,L) that allows constant anisotropy in velocity dispersions characterized by a parameter \beta, we find that the eccentricity distribution is a monotonically increasing function of e for the case of highly radially anisotropic velocity dispersions (\beta > 0.6), while showing a hump-like shape for the cases from radial through tangential velocity anisotropy (\beta < 0.6). We also find that when the velocity anisotropy agrees with that observed for the Milky Way halo stars (\beta = 0.5-0.7), a nearly linear eccentricity distribution of N(e) \alpha e results at e < 0.7, largely independent of the potential adopted. Our theoretical eccentricity distribution would be a vital tool of examining how far out in the halo the dynamical equilibrium has been achieved, through comparison with kinematics of halo stars sampled at greater distances. Given that large surveys of the SEGUE and Gaia projects would be in progress, we discuss how our results would serve as a new guide in exploring the formation and evolution of the Milky Way halo.

We have recently completed a 64-night spectroscopic monitoring campaign at the Lick Observatory 3-m Shane telescope with the aim of measuring the masses of the black holes in 12 nearby (z < 0.05) Seyfert 1 galaxies with expected masses in the range ~10^6-10^7M_sun and also the well-studied nearby active galactic nucleus (AGN) NGC 5548. Nine of the objects in the sample (including NGC 5548) showed optical variability of sufficient strength during the monitoring campaign to allow for a time lag to be measured between the continuum fluctuations and the response to these fluctuations in the broad Hbeta emission, which we have previously reported. We present here the light curves for the Halpha, Hgamma, HeII 4686, and HeI 5876 emission lines and the time lags for the emission-line responses relative to changes in the continuum flux. Combining each emission-line time lag with the measured width of the line in the variable part of the spectrum, we determine a virial mass of the central supermassive black hole from several independent emission lines. We find that the masses are generally consistent within the uncertainties. The time-lag response as a function of velocity across the Balmer line profiles is examined for six of the AGNs. Finally we compare several trends seen in the dataset against the predictions from photoionization calculations as presented by Korista & Goad. We confirm several of their predictions, including an increase in responsivity and a decrease in the mean time lag as the excitation and ionization level for the species increases. Further confirmation of photoionization predictions for broad-line gas behavior will require additional monitoring programs for these AGNs while they are in different luminosity states. [abridged]

We examine whether the spectral energy distribution of optical continuum emission of active galactic nuclei (AGNs) changes during flux variation, based on accurate and frequent monitoring observations of 11 nearby Seyfert galaxies and QSOs carried out in the B, V, and I bands for seven years by the MAGNUM telescope. The multi-epoch flux data in any two different bands obtained on the same night show a very tight linear flux to flux relationship for all target AGNs. The flux of the host galaxy within the photometric aperture is carefully estimated by surface brightness fitting to available high-resolution HST images and MAGNUM images. The flux of narrow emission lines in the photometric bands is also estimated from available spectroscopic data. We find that the non-variable component of the host galaxy plus narrow emission lines for all target AGNs is located on the fainter extension of the linear regression line of multi-epoch flux data in the flux to flux diagram. This result strongly indicates that the spectral shape of AGN continuum emission in the optical region does not systematically change during flux variation. The trend of spectral hardening that optical continuum emission becomes bluer as it becomes brighter, which has been reported by many studies, is therefore interpreted as the domination of the variable component of the nearly constant spectral shape of an AGN as it brightens over the non-variable component of the host galaxy plus narrow lines, which is usually redder than AGN continuum emission.

We present observations and analysis of the broadband afterglow of Swift GRB 071025. Using optical and infrared (RIYJHK) photometry, we derive a photometric redshift of 4.4 < z < 5.2; at this redshift our simultaneous multicolour observations begin at ~30 s after the GRB trigger in the host frame and during the initial rising phase of the afterglow. We associate the light curve peak at 580 s in the observer frame with the formation of the forward shock, giving an estimate of the initial Lorentz factor Gamma_0 ~ 200. The red spectral energy distribution (even in regions not affected by the Lyman-alpha break) provides secure evidence of a large dust column. However, the inferred extinction curve shows a prominent flat component between 2000-3000 Angstroms in the rest-frame, inconsistent with any locally observed template but well-fit by models of dust formed by supernovae. Time-dependent fits to the extinction profile reveal no evidence of dust destruction and limit the decrease in the extinction column to Delta A_3000 < 0.54 mag after t = 50 s in the rest frame. Our observations provide evidence of a transition in dust properties at z~5, in agreement with studies of high-z quasars, and suggest that SN-formed dust continues to dominate the opacity of typical galaxies at this redshift.

(Abridged) We present a two month Suzaku X-ray monitoring of the Seyfert 1 galaxy NGC 5548. The campaign consists of 7 observations. We analyze the response in the opacity of the gas that forms the ionized absorber to ionizing flux variations. Despite variations by a factor of 4 in the impinging continuum, the soft X-ray spectra of the source show little spectral variations, suggesting no response from the ionized absorber. A detailed time modeling confirms the lack of opacity variations for an absorbing component with high ionization. Instead, the models tentatively suggest that the ionization parameter of a low ionization absorbing component might be changing with the ionizing flux, as expected for gas in photoionization equilibrium. Using the lack of variations, we set an upper limit of n_e <2.0E7 cm-3 for the electron density of the gas forming the high ionization, high velocity component. This implies a large distance from the continuum source (R > 0.033 pc). If the variations in the low ionization component are real, they imply n_e >9.8E4 cm-3 and R < 3 pc. We discuss our results in terms of two different scenarios: a large scale outflow originating in the inner parts of the accretion disk, or a thermally driven wind originating much farther out. Given the large distance of the wind, the implied mass outflow rate is also large (Mw > 0.08 Maccr). The associated total kinetic energy deployed by the wind in the host galaxy (>1.2E56 erg) can be enough to disrupt the interstellar medium, possibly regulating large scale star formation. The total mass and energy ejected by the wind is still lower than the one required for cosmic feedback, even when extrapolated to quasar luminosities. Such feedback would require that we are observing the wind before it is fully accelerated.

We present 7 sequential weekly observations of NGC 5548 conducted in 2007 with the \textit{Suzaku} X-ray Imaging Spectrometer (XIS) in the 0.2-12 keV band and Hard X-ray Detector (HXD) in 10-600 keV band. The iron K$\alpha$ line is well detected in all seven observations and K$\beta$ line is also detected in four observations. In this paper, we investigate the origin of the Fe K lines using both the width of the line and the reverberation mapping method. With the co-added XIS and HXD spectra, we identify Fe K$\alpha$ and K$\beta$ line at 6.396$_{-0.007}^{+0.009}$ keV and 7.08$_{-0.05}^{+0.05}$ keV, respectively. The width of line obtained from the co-added spectra is 38$_{-18}^{+16}$ eV ($\textrm{FWHM}=4200_{-2000}^{+1800}$ km/s) which corresponds to a radius of 20$_{-10}^{+50}$ light days, for the virial production of $1.220\times10^7$ M$_{\odot}$ in NCG 5548. To quantitatively investigate the origin of the narrow Fe line by the reverberation mapping method, we compare the observed light curves of Fe K$\alpha$ line with the predicted ones, which are obtained by convolving the continuum light curve with the transfer functions in a thin shell and an inclined disk. The best-fit result is given by the disk case with $i=30^\circ$ which is better than a fit to a constant flux of the Fe K line at the 92.7% level (F-test). We find that the emitting radius obtained from the light curve is 25-37 light days, which is consistent with the radius derived from the Fe K line width. Combining the results of the line width and variation, the most likely site for the origin of the narrow iron lines is 20-40 light days away from the central engine, though other possibilities are not completely ruled out. (abridged)

We discuss the origin of the optical variations in the Narrow line Seyfert 1 galaxy NGC 4051 and present the results of a cross-correlation study using X-ray and optical light curves spanning more than 12 years. The emission is highly variable in all wavebands, and the amplitude of the optical variations is found to be smaller than that of the X-rays, even after correcting for the contaminating host galaxy flux falling inside the photometric aperture. The optical power spectrum is best described by an unbroken power law model with slope $\alpha=1.4^{+0.6}_{-0.2}$ and displays lower variability power than the 2-10 keV X-rays on all time-scales probed. We find the light curves to be significantly correlated at an optical delay of $1.2^{+1.0}_{-0.3}$ days behind the X-rays. This time-scale is consistent with the light travel time to the optical emitting region of the accretion disc, suggesting that the optical variations are driven by X-ray reprocessing. We show, however, that a model whereby the optical variations arise from reprocessing by a flat accretion disc cannot account for all the optical variability. There is also a second significant peak in the cross-correlation function, at an optical delay of $39^{+2.7}_{-8.4}$ days. The lag is consistent with the dust sublimation radius in this source, suggesting that there is a measurable amount of optical flux coming from the dust torus. We discuss the origin of the additional optical flux in terms of reprocessing of X-rays and reflection of optical light by the dust.

We present the detailed optical to far-infrared observations of SST J1604+4304, an ULIRG at z = 1.135. Analyzing the stellar absorption lines, namely, the CaII H & K and Balmer H lines in the optical spectrum, we derive the upper limits of an age for the stellar population. Given this constraint, the minimum {chi}^2 method is used to fit the stellar population models to the observed SED from 0.44 to 5.8um. We find the following properties. The stellar population has an age 40 - 200 Myr with a metallicity 2.5 Z_{sun}. The starlight is reddened by E(B-V) = 0.8. The reddening is caused by the foreground dust screen, indicating that dust is depleted in the starburst site and the starburst site is surrounded by a dust shell. The infrared (8-1000um) luminosity is L_{ir} = 1.78 +/- 0.63 * 10^{12} L_{sun}. This is two times greater than that expected from the observed starlight, suggesting either that 1/2 of the starburst site is completely obscured at UV-optical wavelengths, or that 1/2 of L_{ir} comes from AGN emission. The inferred dust mass is 2.0 +/- 1.0 * 10^8 M_{sun}. This is sufficient to form a shell surrounding the galaxy with an optical depth E(B-V) = 0.8. From our best stellar population model - an instantaneous starburst with an age 40 Myr, we infer the rate of 19 supernovae(SNe) per year. Simply analytical models imply that 2.5 Z_{sun} in stars was reached when the gas mass reduced to 30% of the galaxy mass. The gas metallcity is 4.8 Z_{sun} at this point. The gas-to-dust mass ratio is then 120 +/- 73. The inferred dust production rate is 0.24 +/- 0.12 M_{sun} per SN. If 1/2 of L_{ir} comes from AGN emission, the rate is 0.48 +/- 0.24 M_{sun} per SN. We discuss the evolutionary link of SST J1604+4304 to other galaxy populations in terms of the stellar masses and the galactic winds.

The Lick AGN Monitoring Project targeted 13 nearby Seyfert 1 galaxies with the intent of measuring the masses of their central black holes using reverberation mapping. The sample includes 12 galaxies selected to have black holes with masses roughly in the range 10^6-10^7 solar masses, as well as the well-studied AGN NGC 5548. In conjunction with a spectroscopic monitoring campaign, we obtained broad-band B and V images on most nights from 2008 February through 2008 May. The imaging observations were carried out by four telescopes: the 0.76-m Katzman Automatic Imaging Telescope (KAIT), the 2-m Multicolor Active Galactic Nuclei Monitoring (MAGNUM) telescope, the Palomar 60-in (1.5-m) telescope, and the 0.80-m Tenagra II telescope. Having well-sampled light curves over the course of a few months is useful for obtaining the broad-line reverberation lag and black hole mass, and also allows us to examine the characteristics of the continuum variability. In this paper, we discuss the observational methods and the photometric measurements, and present the AGN continuum light curves. We measure various variability characteristics of each of the light curves. We do not detect any evidence for a time lag between the B- and V-band variations, and we do not find significant color variations for the AGNs in our sample.

A detailed analysis of the data from a high sampling rate, multi-month reverberation mapping campaign, undertaken primarily at MDM Observatory with supporting observations from telescopes around the world, reveals that the Hbeta emission region within the broad line regions (BLRs) of several nearby AGNs exhibit a variety of kinematic behaviors. While the primary goal of this campaign was to obtain either new or improved Hbeta reverberation lag measurements for several relatively low luminosity AGNs (presented in a separate work), we were also able to unambiguously reconstruct velocity-resolved reverberation signals from a subset of our targets. Through high cadence spectroscopic monitoring of the optical continuum and broad Hbeta emission line variations observed in the nuclear regions of NGC 3227, NGC 3516, and NGC 5548, we clearly see evidence for outflowing, infalling, and virialized BLR gas motions, respectively.

We have recently completed a 64-night spectroscopic monitoring campaign at the Lick Observatory 3-m Shane telescope with the aim of measuring the masses of the black holes in 12 nearby (z < 0.05) Seyfert 1 galaxies with expected masses in the range ~10^6-10^7 M_sun and also the well-studied nearby active galactic nucleus (AGN) NGC 5548. Nine of the objects in the sample (including NGC 5548) showed optical variability of sufficient strength during the monitoring campaign to allow for a time lag to be measured between the continuum fluctuations and the response to these fluctuations in the broad Hbeta emission. We present here the light curves for the objects in this sample and the subsequent Hbeta time lags for the nine objects where these measurements were possible. The Hbeta lag time is directly related to the size of the broad-line region, and by combining the lag time with the measured width of the Hbeta emission line in the variable part of the spectrum, we determine the virial mass of the central supermassive black hole in these nine AGNs. The absolute calibration of the black hole masses is based on the normalization derived by Onken et al. We also examine the time lag response as a function of velocity across the Hbeta line profile for six of the AGNs. The analysis of four leads to ambiguous results with relatively flat time lags as a function of velocity. However, SBS 1116+583A exhibits a symmetric time lag response around the line center reminiscent of simple models for circularly orbiting broad-line region (BLR) clouds, and Arp 151 shows an asymmetric profile that is most easily explained by a simple gravitational infall model. Further investigation will be necessary to fully understand the constraints placed on physical models of the BLR by the velocity-resolved response in these objects.

We present the first results from a high sampling rate, multi-month reverberation mapping campaign undertaken primarily at MDM Observatory with supporting observations from telescopes around the world. The primary goal of this campaign was to obtain either new or improved Hbeta reverberation lag measurements for several relatively low luminosity AGNs. We feature results for NGC 4051 here because, until now, this object has been a significant outlier from AGN scaling relationships, e.g., it was previously a ~2-3sigma outlier on the relationship between the broad-line region (BLR) radius and the optical continuum luminosity - the R_BLR-L relationship. Our new measurements of the lag time between variations in the continuum and Hbeta emission line made from spectroscopic monitoring of NGC 4051 lead to a measured BLR radius of R_BLR = 1.87 (+0.54 -0.50) light days and black hole mass of M_BH = 1.73 (+0.55 -0.52) x 10^6 M_sun. This radius is consistent with that expected from the R_BLR-L relationship, based on the present luminosity of NGC 4051 and the most current calibration of the relation by Bentz et al. (2009a). We also present a preliminary look at velocity-resolved Hbeta light curves and time delay measurements, although we are unable to reconstruct an unambiguous velocity-resolved reverberation signal.

The long-term optical and near infrared monitoring observations for a type 1 act ive galactic nucleus NGC 4151 were carried out for six years from 2001 to 2006 b y using the MAGNUM telescope, and delayed response of flux variations in the $K(2.2\mu m)$ band to those in the $V(0.55\mu m)$ band was clearly detected. Based on cross correlation analysis, we precisely measured a lag time $\Delta t$ for eight separate periods, and we found that $\Delta t$ is not constant changing be tween 30 and 70 days during the monitoring period. Since $\Delta t$ is the ligh t travel time from the central energy source out to the surrounding dust torus, this is the first convincing evidence that the inner radius of dust torus did ch ange in an individual AGN. In order to relate such a change of $\Delta t$ with a change of AGN luminosity $L$, we presented a method of taking an average of th e observed $V$-band fluxes that corresponds to the measured value of $\Delta t$, and we found that the time-changing track of NGC 4151 in the $\Delta t$ versus $L$ diagram during the monitoring period deviates from the relation of $\Delta t \propto L^{0.5}$ expected from dust reverberation. This result, combined with t he elapsed time from period to period for which $\Delta t$ was measured, indicat es that the timescale of dust formation is about one year, which should be taken into account as a new constraint in future studies of dust evolution in AGNs.

We present a theoretical model for supernova (SN) 2008D associated with the luminous X-ray transient 080109. The bolometric light curve and optical spectra of the SN are modelled based on the progenitor models and the explosion models obtained from hydrodynamic/nucleosynthetic calculations. We find that SN 2008D is a more energetic explosion than normal core-collapse supernovae, with an ejecta mass of Mej = 5.3 +- 1.0 Msun and a kinetic energy of E = 6.0 +- 2.5 x 10^{51} erg. The progenitor star of the SN has a 6-8 Msun He core with essentially no H envelope (< 5 x 10^{-4} Msun) prior to the explosion. The main-sequence mass of the progenitor is estimated to be Mms =20-25 Msun, with additional systematic uncertainties due to convection, mass loss, rotation, and binary effects. These properties are intermediate between those of normal SNe and hypernovae associated with gamma-ray bursts. The mass of the central remnant is estimated as 1.6 - 1.8 Msun, which is near the boundary between neutron star and black hole formation.

We present spectra of six luminous quasars at z ~ 2, covering rest wavelengths 1600-3200 A. The fluxes of the UV Fe II emission lines and Mg II 2798 doublet, the line widths of Mg II, and the 3000 A luminosity were obtained from the spectra. These quantities were compared with those of low-redshift quasars at z = 0.06 - 0.55 studied by Tsuzuki et al. In a plot of the Fe II(UV)/Mg II flux ratio as a function of the cental black hole mass, Fe II(UV)/Mg II in our z ~ 2 quasars is systematically greater than in the low-redshift quasars. We confermed that luminosity is not responsible for this excess. It is unclear whether this excess is caused by rich Fe abundance at z ~ 2 over low-redshift or by non-abundance effects such as high gas density, strong radiation field, and high microturbulent velocity.

We present the results of concurrent X-ray and optical monitoring of the Seyfert 1 galaxy Mrk 79 over a period of more than five years. We find that on short to medium time-scales (days to a few tens of days) the 2-10 keV X-ray and optical u and V band fluxes are significantly correlated, with a delay between the bands consistent with zero days. We show that most of these variations may be well reproduced by a model where the short-term optical variations originate from reprocessing of X-rays by an optically thick accretion disc. The optical light curves, however, also display long time-scale variations over thousands of days, which are not present in the X-ray light curve. These optical variations must originate from an independent variability mechanism and we show that they can be produced by variations in the (geometrically) thin disc accretion rate as well as by varying reprocessed fractions through changes in the location of the X-ray corona.

We have recently completed a 64-night spectroscopic monitoring campaign at the Lick Observatory 3-m Shane telescope with the aim of measuring the masses of the black holes in 13 nearby (z < 0.05) Seyfert 1 galaxies with expected masses in the range ~10^6-10^7 M_sun. We present here the first results from this project -- the mass of the central black hole in Arp 151. Strong variability throughout the campaign led to an exceptionally clean Hbeta lag measurement in this object of 4.25(+0.68/-0.66) days in the observed frame. Coupled with the width of the Hbeta emission line in the variable spectrum, we determine a black hole mass of (7.1 +/- 1.2)x10^6 M_sun, assuming the Onken et al. normalization for reverberation-based virial masses. We also find velocity-resolved lag information within the Hbeta emission line which clearly shows infalling gas in the Hbeta-emitting region. Further detailed analysis may lead to a full model of the geometry and kinematics of broad line region gas around the central black hole in Arp 151.

We investigate the dynamical response, in terms of disc size and rotation velocity, to mass loss by supernovae in the evolution of spiral galaxies. A thin baryonic disc having the Kuzmin density profile embedded in a spherical dark matter halo having a density profile proposed by Navarro, Frenk & White is considered. For a purpose of comparison, we also consider the homogeneous and $r^{-1}$ profiles for dark matter in a truncated spherical halo. Assuming for simplicity that the dark matter distribution is not affected by mass loss from discs and the change of baryonic disc matter distribution is homologous, we evaluate the effects of dynamical response in the resulting discs. We found that the dynamical response only for an adiabatic approximation of mass loss can simultaneously account for the rotation velocity and disc size as observed particularly in dwarf spiral galaxies, thus reproducing the Tully-Fisher relation and the size versus magnitude relation over the full range of magnitude. Furthermore, we found that the mean specific angular momentum in discs after the mass loss becomes larger than that before the mass loss, suggesting that the mass loss would occur most likely from the central disc region where the specific angular momentum is low.

We present a theoretical model for Type Ib supernova (SN) 2006jc. We calculate the evolution of the progenitor star, hydrodynamics and nucleosynthesis of the SN explosion, and the SN bolometric light curve (LC). The synthetic bolometric LC is compared with the observed bolometric LC constructed by integrating the UV, optical, near-infrared (NIR), and mid-infrared (MIR) fluxes. The progenitor is assumed to be as massive as $40M_\odot$ on the zero-age main-sequence. The star undergoes extensive mass loss to reduce its mass down to as small as $6.9M_\odot$, thus becoming a WCO Wolf-Rayet star. The WCO star model has a thick carbon-rich layer, in which amorphous carbon grains can be formed. This could explain the NIR brightening and the dust feature seen in the MIR spectrum. We suggest that the progenitor of SN 2006jc is a WCO Wolf-Rayet star having undergone strong mass loss and such massive stars are the important sites of dust formation. We derive the parameters of the explosion model in order to reproduce the bolometric LC of SN 2006jc by the radioactive decays: the ejecta mass $4.9M_\odot$, hypernova-like explosion energy $10^{52}$ ergs, and ejected $^{56}$Ni mass $0.22M_\odot$. We also calculate the circumstellar interaction and find that a CSM with a flat density structure is required to reproduce the X-ray LC of SN 2006jc. This suggests a drastic change of the mass-loss rate and/or the wind velocity that is consistent with the past luminous blue variable (LBV)-like event.

We present new spectroscopic and photometric observations of the HAT-P-1 planetary system. Spectra obtained during three transits exhibit the Rossiter-McLaughlin effect, allowing us to measure the angle between the sky projections of the stellar spin axis and orbit normal, \lambda = 3.7 +/- 2.1 degrees. The small value of \lambda for this and other systems suggests that the dominant planet migration mechanism preserves spin-orbit alignment. Using two new transit light curves, we refine the transit ephemeris and reduce the uncertainty in the orbital period by an order of magnitude. We find a upper limit on the orbital eccentricity of 0.067, with 99% confidence, by combining our new radial-velocity measurements with those obtained previously.

We present our latest results on near- to mid- infrared observation of SN2006jc at 200 days after the discovery using the Infrared Camera (IRC) on board $AKARI$. The near-infrared (2--5$\mu$m) spectrum of SN2006jc is obtained for the first time and is found to be well interpreted in terms of the thermal emission from amorphous carbon of 800$\pm 10$K with the mass of $6.9\pm 0.5 \times 10^{-5}M_{\odot}$ that was formed in the supernova ejecta. This dust mass newly formed in the ejecta of SN 2006jc is in a range similar to those obtained for other several dust forming core collapse supernovae based on recent observations (i.e., $10^{-3}$--$10^{-5}$$M_{\odot}$). Mid-infrared photometric data with {\it{AKARI}}/IRC MIR-S/S7, S9W, and S11 bands have shown excess emission over the thermal emission by hot amorphous carbon of 800K. This mid-infrared excess emission is likely to be accounted for by the emission from warm amorphous carbon dust of 320$\pm 10$K with the mass of 2.7$^{+0.7}_{-0.5} \times 10^{-3}M_{\odot}$ rather than by the band emission of astronomical silicate and/or silica grains. This warm amorphous carbon dust is expected to have been formed in the mass loss wind associated with the Wolf-Rayet stellar activity before the SN explosion. Our result suggests that a significant amount of dust is condensed in the mass loss wind prior to the SN explosion. A possible contribution of emission bands by precursory SiO molecules in 7.5--9.5$\mu$m is also suggested.

We report a measurement of the Rossiter--McLaughlin effect in the transiting extrasolar planetary system TrES-1, via simultaneous spectroscopic and photometric observations with the Subaru and MAGNUM telescopes. By modeling the radial velocity anomaly that was observed during a transit, we determine the sky-projected angle between the stellar spin axis and the planetary orbital axis to be $\lambda = 30 \pm 21$ [deg]. This is the third case for which $\lambda$ has been measured in a transiting exoplanetary system, and the first demonstration that such measurements are possible for relatively faint host stars ($V \sim 12$, as compared to $V \sim 8$ for the other systems). We also derive a time of mid-transit, constraints on the eccentricity of the TrES-1b orbit ($e = 0.048 \pm 0.025$), and upper limits on the mass of the Trojan companions ($\lesssim$14 $M_{\oplus}$) at the 3$\sigma$ level.

We present the most successful infrared cloud monitor for a robotic telescope. This system was originally developed for the MAGNUM 2-m telescope, which has been achieving unmanned and automated monitoring observation of active galactic nuclei at Haleakala on the Hawaiian island of Maui since 2001. Using a thermal imager and two aspherical mirrors, it at once sees almost the whole sky at a wavelength of $\lambda\sim 10\mu{\rm m}$. Its outdoor part is weather-proof and is totally maintenance-free. The images obtained every one or two minutes are analysed immediately into several ranks of weather condition, from which our automated observing system not only decides to open or close the dome, but also selects what types of observations should be done. The whole-sky data accumulated over four years show that 50$-$60 % of all nights are photometric, and about 75 % are observable with respect to cloud condition at Haleakala. Many copies of this system are now used all over the world such as Mauna Kea in Hawaii, Atacama in Chile, and Okayama and Kiso in Japan.

We describe results from a new ground-based monitoring campaign on NGC 5548, the best studied reverberation-mapped AGN. We find that it was in the lowest luminosity state yet recorded during a monitoring program, namely L(5100) = 4.7 x 10^42 ergs s^-1. We determine a rest-frame time lag between flux variations in the continuum and the Hbeta line of 6.3 (+2.6/-2.3) days. Combining our measurements with those of previous campaigns, we determine a weighted black hole mass of M_BH = 6.54 (+0.26/-0.25) x 10^7 M_sun based on all broad emission lines with suitable variability data. We confirm the previously-discovered virial relationship between the time lag of emission lines relative to the continuum and the width of the emission lines in NGC 5548, which is the expected signature of a gravity-dominated broad-line region. Using this lowest luminosity state, we extend the range of the relationship between the luminosity and the time lag in NGC 5548 and measure a slope that is consistent with alpha = 0.5, the naive expectation for the broad line region for an assumed form of r ~ L^alpha. This value is also consistent with the slope recently determined by Bentz et al. for the population of reverberation-mapped AGNs as a whole.

We present photometry of HD 189733 during eight transits of its close-in giant planet, and out-of-transit photometry spanning two years. Using the transit photometry, we determine the stellar and planetary radii and the photometric ephemeris. Outside of transits, there are quasiperiodic flux variations with a 13.4 day period that we attribute to stellar rotation. In combination with previous results, we derive upper limits on the orbital eccentricity, and on the true angle between the stellar rotation axis and planetary orbit (as opposed to the angle between the projections of those axes on the sky).

We propose a new method of analysing a variable component for type 1 active galactic nuclei (AGNs) in the near-infrared wavelength region. This analysis uses a multiple regression technique and divides the variable component into two components originating in the accretion disk at the center of AGNs and from the dust torus that far surrounds the disk. Applying this analysis to the long-term $VHK$ monitoring data of MCG+08-11-011 that were obtained by the MAGNUM project, we found that the $(H-K)$-color temperature of the dust component is $T = 1635$K $\pm20$K, which agrees with the sublimation temperature of dust grains, and that the time delay of $K$ to $H$ variations is $\Delta t\approx 6$ days, which indicates the existence of a radial temperature gradient in the dust torus. As for the disk component, we found that the power-law spectrum of $f_\nu \propto \nu^\alpha$ in the $V$ to near-infrared $HK$ bands varies with a fixed index of $\alpha\approx -0.1$ -- +0.4, which is broadly consistent with the irradiated standard disk model. The outer part of the disk therefore extends out to a radial distance where the temperature decreases to radiate the light in the near-infrared.

We present new observations leading to an improved black hole mass estimate for the Seyfert 1 galaxy NGC 4593 as part of a reverberation-mapping campaign conducted at the MDM Observatory. Cross-correlation analysis of the H_beta emission-line light curve with the optical continuum light curve reveals an emission-line time delay of 3.73 (+-0.75) days. By combining this time delay with the H_beta line width, we derive a central black hole mass of M_BH = 9.8(+-2.1)x10^6 M_sun, an improvement in precision of a factor of several over past results.

We have undertaken a new ground-based monitoring campaign to improve the estimates of the mass of the central black hole in NGC 4151. We measure the lag time of the broad H beta line response compared to the optical continuum at 5100 A and find a lag of 6.6 (+1.1/-0.8) days. We combine our data with the recent reanalysis of UV emission lines by Metzroth et al. to calculate a weighted mean of the black hole mass, M_BH = 4.57 (+0.57/-0.47) x 10^7 M_sun. The absolute calibration of the black hole mass is based on normalization of the AGN black hole mass - stellar velocity dispersion (M_BH - sigma_*) relationship to that of quiescent galaxies by Onken et al. The scatter in the M_BH - sigma_* relationship suggests that reverberation-mapping based mass measurements are typically uncertain by a factor of 3-4.

We present the spectra of 14 quasars with a wide coverage of rest wavelengths from 1000 to 7300 A. The redshift ranges from z = 0.061 to 0.555 and the luminosity from M_{B} = -22.69 to -26.32. We describe the procedure of generating the template spectrum of Fe II line emission from the spectrum of a narrow-line Seyfert 1 galaxy I Zw 1 that covers two wavelength regions of 2200-3500 A and 4200-5600 A. Our template Fe II spectrum is semi-empirical in the sense that the synthetic spectrum calculated with the CLOUDY photoionization code is used to separate the Fe II emission from the Mg II line. The procedure of measuring the strengths of Fe II emission lines is twofold; (1) subtracting the continuum components by fitting models of the power-law and Balmer continua in the continuum windows which are relatively free from line emissions, and (2) fitting models of the Fe II emission based on the Fe II template to the continuum-subtracted spectra. From 14 quasars, we obtained the Fe II fluxes in five wavelength bands, the total flux of Balmer continuum, and the fluxes of Mg II, Halpha, and other emission lines, together with the full width at half maxima (FWHMs) of these lines. Regression analysis was performed by assuming a linear relation between any two of these quantities. Eight correlations were found with a confidence level higher than 99%. The fact that six of these eight are related to FWHM or M_{BH} may imply that M_{BH} is a fundamental quantity that controls Gamma or the spectral energy distribution (SED) of the incident continuum, which in turn controls the Fe II emission. Furthermore, it is worthy of noting that Fe II(O1)/Fe II(U1) is found to tightly correlate with Fe II(O1)/Mg II, but not with Fe II(U1)/Mg II.

We present results from Chandra observations of SDSS J1004+4112, a strongly lensed quasar system with a maximum image separation of 15". All four bright images of the quasar, as well as resolved X-ray emission originating from the lensing cluster, are clearly detected. The emission from the lensing cluster extends out to approximately 1.5 arcmin. We measure the bolometric X-ray luminosity and temperature of the lensing cluster to be 4.7e44 erg s^-1 and 6.4 keV, consistent with the luminosity-temperature relation for distant clusters. The mass estimated from the X-ray observation shows excellent agreement with the mass derived from gravitational lensing. The X-ray flux ratios of the quasar images differ markedly from the optical flux ratios, and the combined X-ray spectrum of the images possesses an unusually strong Fe Kalpha emission line, both of which are indicative of microlensing.

We carried out a one-night optical V and near-infrared JHK monitoring observation of the least luminous Seyfert 1 galaxy, NGC4395, on 2004 May 1, and detected for the first time the intraday flux variations in the J and H bands, while such variation was not clearly seen for the K band. The detected J and H variations are synchronized with the flux variation in the V band, which indicates that the intraday-variable component of near-infrared continuum emission of the NGC4395 nucleus is an extension of power-law continuum emission to the near-infrared and originates in an outer region of the central accretion disk. On the other hand, from our regular program of long-term optical BVI and near-infrared JHK monitoring observation of NGC4395 from 2004 February 12 until 2005 January 22, we found large flux variations in all the bands on time scales of days to months. The optical BVI variations are almost synchronized with each other, but not completely with the near-infrared JHK variations. The color temperature of the near-infrared variable component is estimated to be T=1320-1710 K, in agreement with thermal emission from hot dust tori in active galactic nuclei (AGNs). We therefore conclude that the near-infrared variation consists of two components having different time scales, so that a small K-flux variation on a time scale of a few hours would possibly be veiled by large variation of thermal dust emission on a time scale of days.

We report near simultaneous multi-color (RIYJHK) observations made with the MAGNUM 2m telescope of the gamma ray burst GRB 050904 detected by the SWIFT satellite. The spectral energy distribution shows a very large break between the I and J bands. Using intergalactic transmissions measured from high redshift quasars we show that the observations place a 95% confidence lower limit of z=6.18 on the object, consistent with a later measured spectroscopic redshift of 6.29 obtained by Kawai et al. (2005) with the Subaru telescope. We show that the break strength in the R and I bands is consistent with that measured in the quasars. Finally we consider the implications for the star formation history at high redshift.

Late-time BVRIJHK photometry of the peculiar Type Ic SN 2002ap, taken between 2002 June 12 and 2003 August 29 with the MAGNUM telescope, is presented. The light curve decline rate is derived in each band and the color evolution is studied through comparison with nebular spectra and with SN 1998bw. Using the photometry, the OIR bolometric light curve is built, extending from before light maximum to day 580 after explosion. The light curve has a late-time shape strikingly similar to that of the hypernova SN 1998bw. The decline rate changes from 0.018 mag/day between day 130 and 230 to 0.014 mag/day between day 270 and 580. To reproduce the late-time light curve, a dense core must be added to the 1-D hypernova model that best fits the early-time observations, bringing the ejecta mass from 2.5 Msun to 3 Msun without much change in the kinetic energy, which is 4 times 10^51 ergs. This is similar to the case of other hypernovae and suggests asymmetry. A large H-band bump developed in the spectral energy distribution after about day 300, probably caused by strong [Si I] 1.646 micron and 1.608 micron emissions. The near-infrared flux contribution increased simultaneously from <30% to >50% at day 580. The near-infrared light curves were compared with those of other Type Ib/c supernovae, among which SN 1983I seems similar to SN 2002ap both in the near-infrared and in the optical.

The most intense monitoring observations yet made in the optical (UBV) and near-infrared (JHK) wave bands were carried out for nearby Seyfert1 galaxies of NGC 5548, NGC 4051, NGC 3227, and NGC 7469. Over three years of observations with MAGNUM telescope since early 2001, clear time-delayed response of the K-band flux variations to the V-band flux variations was detected for all of these galaxies. Their H-K color temperature was estimated to be 1500-1800 K from the observed flux variation gradients, which supports a view that the bulk of the K flux should originate in the thermal radiation of hot dust that surrounds the central engine. Cross-correlation analysis was performed to quantify the lag time corresponding to the light-travel distance of the hot dust region from the central engine. The measured lag time is 47-53 days for NGC 5548, 11-18 days for NGC 4051, about 20 days for NGC 3227, and 65-87 days for NGC 7469. We found that the lag time is tightly correlated with the optical luminosity as expected from dust reverberation ($\Delta t \propto L^{0.5}$), while only weakly with the central virial mass, which suggests that an inner radius of the dust torus around the active nucleus has a one-to-one correspondence to central luminosity. In the lag time versus central luminosity diagram, the K-band lag times place an upper boundary on the similar lag times of broad-emission lines in the literature. This not only supports the unified scheme of AGNs, but also implies a physical transition from the BLR out to the dust torus that encircles the BLR. Furthermore, our V-band flux variations of NGC 5548 on timescales of up to 10 days are found to correlate with X-ray variations and delay behind them by one or two days, indicating the thermal reprocessing of X-ray emission by the central accretion flow.

Observations and modeling for the light curve (LC) and spectra of supernova (SN) 2005bf are reported. This SN showed unique features: the LC had two maxima, and declined rapidly after the second maximum, while the spectra showed strengthening He lines whose velocity increased with time. The double-peaked LC can be reproduced by a double-peaked $^{56}$Ni distribution, with most $^{56}$Ni at low velocity and a small amount at high velocity. The rapid post-maximum decline requires a large fraction of the $\gamma$-rays to escape from the $^{56}$Ni-dominated region, possibly because of low-density ``holes''. The presence of Balmer lines in the spectrum suggests that the He layer of the progenitor was substantially intact. Increasing $\gamma$-ray deposition in the He layer due to enhanced $\gamma$-ray escape from the $^{56}$Ni-dominated region may explain both the delayed strengthening and the increasing velocity of the He lines. The SN has massive ejecta ($\sim6-7\Msun$), normal kinetic energy ($\sim 1.0-1.5\times 10^{51}$ ergs), high peak bolometric luminosity ($\sim 5\times 10^{42}$ erg s$^{-1}$) for an epoch as late as $\sim$ 40 days, and a large $^{56}$Ni mass ($\sim0.32\Msun$). These properties, and the presence of a small amount of H suggest that the progenitor was initially massive (M$\sim 25-30 \Msun$) and had lost most of its H envelope, and was possibly a WN star. The double-peaked $^{56}$Ni distribution suggests that the explosion may have formed jets that did not reach the He layer. The properties of SN 2005bf resemble those of the explosion of Cassiopeia A.

We describe the discovery of HE 1327-2326, a dwarf or subgiant with [Fe/H]}=-5.4. The star was found in a sample of bright metal-poor stars selected from the Hamburg/ESO survey. Its abundance pattern is characterized by very high C and N abundances. The detection of Sr which is overabundant by a factor of 10 as compared to iron and the Sun, suggests that neutron-capture elements had already been produced in the very early Galaxy. A puzzling Li depletion is observed in this unevolved star which contradicts the value of the primordial Li derived from WMAP and other Li studies. Possible scenarios for the origin of the abundance pattern (Pop. II or Pop. III) are presented as well as an outlook on future observations.

We present the elemental abundances of HE1327-2326, the most iron-deficient star known, determined from a comprehensive analysis of spectra obtained with the Subaru Telescope High Dispersion Spectrograph.

We construct the Numerical Galaxy Catalog ($\nu$GC), based on a semi-analytic model of galaxy formation combined with high-resolution N-body simulations in a $\Lambda$-dominated flat cold dark matter ($\Lambda$CDM) cosmological model. The model includes several essential ingredients for galaxy formation, such as merging histories of dark halos directly taken from N-body simulations, radiative gas cooling, star formation, heating by supernova explosions (supernova feedback), mergers of galaxies, population synthesis, and extinction by internal dust and intervening HI clouds. As the first paper in a series using this model, we focus on basic photometric, structural and kinematical properties of galaxies at present and high redshifts. Two sets of model parameters are examined, strong and weak supernova feedback models, which are in good agreement with observational luminosity functions of local galaxies in a range of observational uncertainty. Both models agree well with many observations such as cold gas mass-to-stellar luminosity ratios of spiral galaxies, HI mass functions, galaxy sizes, faint galaxy number counts and photometric redshift distributions in optical pass-bands, isophotal angular sizes, and cosmic star formation rates. In particular, the strong supernova feedback model is in much better agreement with near-infrared (K'-band) faint galaxy number counts and redshift distribution than the weak feedback model and our previous semi-analytic models based on the extended Press-Schechter formalism. (Abridged)

By observing the near-infrared spectrum of the quasar PG 1116+215 at z = 0.176 and combining with the HST/FOS spectrum, we obtained the relative strengths of three permitted OI lines ({lambda}1304, {lambda}8446, and {lambda}11287) in a quasar for the first time. The photon flux ratios of the OI lines of the quasar were compared with those previously measured in a Seyfert 1 and six narrow-line Seyfert 1s. No significant differences were found in the OI line flux ratios between the quasar and the other Seyferts, suggesting that the gas density in the OI and FeII line-emitting regions in the quasar is of the same order as those in low-luminosity AGNs. It was also found that the line width of OI {lambda}11287 is significantly narrower than that of Ly{alpha}, which is consistent with OI and FeII emission occurring in the partly ionized regions at the outermost portion of the broad-line region where velocities are small.

We present a deep $K^{\prime}$-band (2.12$\mu$m) imaging of 1\arcmin\ $\times$ 1\arcmin Subaru Super Deep Field (SSDF) taken with the Subaru adaptive optics (AO) system. Total integration time of 26.8 hours results in the limiting magnitude of $K^{\prime} \sim 24.7$ (5$\sigma$, 0\farcs2 aperture) for point sources and $K^{\prime} \sim 23.5$ (5$\sigma$, 0\farcs6 aperture) for galaxies, which is the deepest limit ever achieved in the $K^{\prime}$ band. The average stellar FWHM of the co-added image is 0\farcs18. Based on the photometric measurements of detected galaxies, we obtained the differential galaxy number counts, for the first time, down to $K^{\prime} \sim 25$, which is more than 0.5 mag deeper than the previous data. We found that the number count slope $d\log N/dm$ is about 0.15 at $22 < K^{\prime} < 25$, which is flatter than the previous data. Therefore, detected galaxies in the SSDF have only negligible contribution to the near-infrared extragalactic background light (EBL), and the discrepancy claimed so far between the diffuse EBL measurements and the estimated EBL from galaxy count integration has become more serious . The size distribution of detected galaxies was obtained down to the area size of less than 0.1 arcsec$^2$, which is less than a half of the previous data in the $K^{\prime}$ band. We compared the observed size-magnitude relation with a simple pure luminosity evolution model allowing for intrinsic size evolution, and found that a model with no size evolution gives the best fit to the data. It implies that the surface brightness of galaxies at high redshift is not much different from that expected from the size-luminosity relation of present-day galaxies.

We present optical, near-IR, and radio follow up of sixteen Swift bursts, including our discovery of nine afterglows and a redshift determination for three. These observations, supplemented by data from the literature, provide an afterglow recovery rate of 60% in the optical/near-IR, much higher than in previous missions (BeppoSAX, HETE-2, INTEGRAL, and IPN). The optical/near-IR afterglows of Swift events are on average 1.7 mag fainter at t=12 hr than those of previous missions. The X-ray afterglows are similarly fainter compared to those of pre-Swift bursts. In the radio the limiting factor is the VLA threshold and the detection rate for Swift bursts is similar to that for past missions. The redshift distribution of pre-Swift bursts peaked at z~1, whereas the five Swift bursts with measured redshifts are distributed evenly between 1.3 and 3.2. From these results we conclude that (i) the pre-Swift distributions were biased in favor of bright events and low redshift events, (ii) the higher sensitivity and accurate positions of Swift result in a better representation of the true burst redshift and brightness distributions (which are higher and dimmer, respectively), and (iii) as many as 1/3 of the bursts can be optically dark, as a result of a high redshift and/or dust extinction. We remark that the apparent lack of low redshift, low luminosity Swift bursts, and the lower event rate compared to pre-launch estimates (90 vs. 150 per year), are the result of a threshold that is similar to that of BATSE. In view of these inferences, afterglow observers may find it advisable to make significant changes in follow up strategies of Swift events. [abridged]

The chemically most primitive stars provide constraints on the nature of the first stellar objects that formed in the Universe; elements other than hydrogen, helium and traces of lithium within these objects were generated by nucleosynthesis in the very first stars. The relative abundances of elements in the surviving primitive stars reflect the masses of the first stars, because the pathways of nucleosynthesis are quite sensitive to stellar masses. Several models have been suggested to explain the origin of the abundance pattern of the giant star HE 0107-5240, which hitherto exhibited the highest deficiency of heavy elements known. Here we report the discovery of HE 1327-2326, a subgiant or main-sequence star with an iron abundance about a factor of two lower than that of HE 0107-5240. Both stars show extreme overabundances of carbon and nitrogen with respect to iron, suggesting a similar origin of the abundance patterns. The unexpectedly low Li and high Sr abundances of HE 1327-2326, however, challenge existing theoretical understanding: none predicts the high Sr abundance or provides a Li depletion mechanism consistent with data available for the most metal-poor stars.

Near-infrared spectra of four QSOs located at $z>6$ are obtained with the OH-airglow suppressor mounted on the Subaru telescope. The FeII/MgII emission-line ratios of these QSOs are examined by the same fitting algorithm as in our previous study of $z<5.3$ QSOs. The fitting results show that two out of the four $z>6$ QSOs have significant FeII emission in their rest-UV spectra, while the other two have almost no FeII features. We also applied our fitting algorithm to more than 10,000 SDSS QSOs and found two trends in the distribution of FeII/MgII against redshift: (1) the upper envelope of the FeII/MgII distribution at $z>3$ shows a probable declination toward high redshift, and (2) the median distribution settles into lower ratios at $z\sim 1.5$ with small scatter compared to the other redshift. We discuss an Fe/Mg abundance evolution of QSOs with a substantial contribution from the diverse nature of the broad-line regions in high-redshift QSOs.

Optical and near-infrared observations of the gamma-ray burst GRB 031203, at z = 0.1055, are reported. A very faint afterglow is detected superimposed to the host galaxy in our first infrared JHK observations, carried out ~9 hours after the burst. Subsequently, a rebrightening is detected in all bands, peaking in the R band about 18 rest-frame days after the burst. The rebrightening closely resembles the light curve of a supernova like SN 1998bw, assuming that the GRB and the SN went off almost simultaneously, but with a somewhat slower evolution. Spectra taken close to the maximum of the rebrightening show extremely broad features as in SN 1998bw. The determination of the absolute magnitude of this SN (SN 2003lw) is difficult owing to the large and uncertain extinction, but likely this event was brighter than SN 1998bw by 0.5 mag in the VRI bands, reaching an absolute magnitude M_V = -19.75+-0.15.

We reanalyze the formation and evolution of galaxies in the hierarchical clustering scenario. Using a semi-analytic model (SAM) of galaxy formation described in this paper, which we hereafter call the Mitaka model, we extensively investigate the observed scaling relations of galaxies among photometric, kinematic, structural and chemical characteristics. In such a scenario, spheroidal galaxies are assumed to be formed by major merger and subsequent starburst, in contrast to the traditional scenario of monolithic cloud collapse. As a new ingredient of SAMs, we introduce the effects of dynamical response to supernova-induced gas removal on size and velocity dispersion, which play an important role on dwarf galaxy formation. In previous theoretical studies of dwarf galaxies based on the monolithic cloud collapse given by Yoshii & Arimoto and Dekel & Silk, the dynamical response was treated in the extremes of a purely baryonic cloud and a baryonic cloud fully supported by surrounding dark matter. To improve this simple treatment, in our previous paper, we formulated the dynamical response in more realistic, intermediate situations between the above extremes. While the effects of dynamical response depend on the mass fraction of removed gas from a galaxy, how much amount of the gas remains just after major merger depends on the star formation history. A variety of star formation histories are generated through the Monte Carlo realization of merging histories of dark halos, and it is found that our SAM naturally makes a wide variety of dwarf galaxies and their dispersed characteristics as observed. (Abridged)

Optical spectroscopy of the Type Ia supernova SN 2002ic obtained on 2003 June 27.6 UT, i.e., about 222 rest-frame days after explosion, is presented. Strong H emission indicates an interaction between the expanding SN ejecta and an H-rich circumstellar medium (CSM). The spectrum of SN 2002ic resembles those of SNe 1997cy and 1999E. The three SNe also have similar luminosities, suggesting that they are the same phenomenon and that the CSM is also similar. We propose a new classification, Type IIa SNe, for these events. The observed line profiles and line ratios are measured and discussed within the ejecta-CSM interaction scenario. The emission in H Balmer, [O III], and He I lines, and in permitted Fe II blends, resembles the spectra of the Type IIn SN 1987F and of Seyfert 1 galaxies. A high-density, clumpy CSM is inferred. Strong, very broad [Ca II]/Ca II and [O I]/O I missions imply that not all the outer SN ejecta were decelerated in the interaction, suggesting that the CSM is aspherical.

The mass function of dark halos in a Lambda-dominated cold dark matter (LambdaCDM) universe is investigated. 529 output files from five runs of N-body simulations are analyzed using the friends-of-friends cluster finding algorithm. All the runs use 512^3 particles in the box size of 35 h^{-1}Mpc to 140 h^{-1}Mpc. Mass of particles for 35 h^{-1} Mpc runs is 2.67 times 10^7 h^{-1} M_{solar}. Because of the high mass resolution of our simulations, the multiplicity function in the low-mass range, where the mass is well below the characteristic mass and $\nu = \delta_c / \sigma \lesssim 1.0$, is evaluated in the present work, and is well fitted by the functional form proposed by Sheth & Tormen (ST). However, the maximum value of the multiplicity function from our simulations at $\nu \sim 1$ is smaller, and its low mass tail is shallower when compared with the ST multiplicity function.

The most intense monitoring observations yet made were carried out on the Seyfert 1 galaxy NGC 4151 in the optical and near-infrared wave-bands. A lag from the optical light curve to the near-infrared light curve was measured. The lag-time between the V and K light curves at the flux minimum in 2001 was precisely 48+2-3 days, as determined by a cross-correlation analysis. The correlation between the optical luminosity of an active galactic nucleus (AGN) and the lag-time between the UV/optical and the near-infrared light curves is presented for NGC 4151 in combination with previous lag-time measurements of NGC 4151 and other AGNs in the literature. This correlation is interpreted as thermal dust reverberation in an AGN, where the near-infrared emission from an AGN is expected to be the thermal re-radiation from hot dust surrounding the central engine at a radius where the temperature equals to that of the dust sublimation temperature. We find that the inner radius of the dust torus in NGC 4151 is $\sim $ 0.04 pc corresponding to the measured lag-time, well outside the broad line region (BLR) determined by other reverberation studies of the emission lines.

The initial mass function (IMF) of metal-free stars that form in the initial starburst of massive (virial temperatures >10^4K) metal-free protogalaxies is studied. In particular, we focus on the effect of H2 photodissociation by pre-existing stars on the fragmentation mass scale, presumedly determined by the Jeans mass at the end of the initial free-fall phase, i.e., at the so-called ``loitering phase,'' characterized by the local temperature minimum. Photodissociation diminishes the Jeans mass at the loitering phase, thereby reducing the fragmentation mass scale of primordial clouds. Thus, in a given cloud, far ultraviolet (FUV) radiation from the first star, which is supposedly very massive (about 10^3Msun), reduces the mass scale for subsequent fragmentation. Through a series of similar processes the IMF for metal-free stars is established. If FUV radiation exceeds a threshold level, the star-forming clumps collapse solely through atomic cooling. Correspondingly, the fragmentation scale drops discontinuously from a few time 10Msun to sub-solar scales. In compact clouds (>1.6kpc for clouds of gas mass 10^8Msun), this level of radiation field is attained, and sub-solar mass stars are formed even in a metal-free environment. Consequently, the IMF becomes bi-modal, with peaks at a few tenths of Msun and a few times 10Msun. The high-mass portion of the IMF is found to be a very steep function of the stellar mass, xi_high(m) being proportinal to m^{-5}. Therefore, the typical mass scale of metal-free stars is significantly smaller than that of the very first stars. Also we study the thermal instability in collapsing primordial prestellar clumps, and discuss why the thermal instability occuring during the three-body H2 formation does not appear to manifest itself in causing further fragmentation of such clumps.

The light curves of 'hypernovae', i.e. very energetic supernovae with $E_{51} \equiv E/10^{51}$ergs $\gsim 5-10$ are characterized at epochs of a few months by a phase of linear decline. Classical, one-dimensional explosion models fail to simultaneously reproduce the light curve near peak and at the linear decline phase. The evolution of these light curves may however be explained by a simple model consisting of two concentric components. The outer component is responsible for the early part of the light curve and for the broad absorption features observed in the early spectra of hypernovae, similar to the one-dimensional models. In addition, a very dense inner component is added, which reproduces the linear decline phase in the observed magnitude-versus-time relation for SNe 1998bw, 1997ef, and 2002ap. This simple approach does contain one of the main features of jet-driven, asymmetric explosion models, namely the presence of a dense core. Although the total masses and energies derived with the two-component model are similar to those obtained in previous studies which also adopted spherical symmetry, this study suggests that the ejecta are aspherical, and thus the real energies and masses may deviate from those derived assuming spherical symmetry. The supernovae which were modeled are divided into two groups, according to the prominence of the inner component: the inner component of SN 1997ef is denser and more $^{56}$Ni-rich, relative to the outer component, than the corresponding inner components of SNe 1998bw and 2002ap. These latter objects have a similar inner-to-outer component ratio, although they have very different global values of mass and energy.

Supernova (SN) 2002ap in M74 was observed in the $UBVRIJHK$ bands for the first 40 days following its discovery (2002 January 29) until it disappeared because of solar conjunction, and then in June after it reappeared. The magnitudes and dates of peak brightness in each band were determined. While the rate of increase of the brightness before the peak is almost independent of wavelength, the subsequent rate of decrease becomes smaller with wavelength from the $U$ to the $R$ band, and is constant at wavelengths beyond $I$. The photometric evolution is faster than in the well-known ``hypernovae'' SNe~1998bw and 1997ef, indicating that SN 2002ap ejected less mass. The bolometric light curve of SN 2002ap for the full period of observations was constructed. The absolute magnitude is found to be much fainter than that of SN 1998bw, but is similar to that of SN 1997ef, which lies at the faint end of the hypernova population. The bolometric light curve at the early epochs was best reproduced with the explosion of a C+O star that ejects $2.5~M_\sun$ with kinetic energy $E_{\rm K}=4\times 10^{51}~{\rm ergs}$. A comparison of the predicted brightness of SN 2002ap with that observed after solar conjunction may imply that $\gamma$-ray deposition at the later epochs was more efficient than in the model. This may be due to an asymmetric explosion.

The external pollution of the first stars in the Galaxy is investigated. The first stars were born in clouds composed of the pristine gas without heavy elements. These stars accreted gas polluted with heavy elements while they still remained in the cloud. As a result, it is found that they exhibit a distribution with respect to the surface metallicity. We have derived the actual form of this distribution function. This metallicity distribution function strongly suggests that the recently discovered most metal-deficient star HE0107-5240 with [Fe/H]=-5.3 was born as a metal-free star and accreted gas polluted with heavy elements. Thus the heavy elements such as Fe in HE0107-5240 must have been supplied from supernovae of later generations exploding inside the cloud in which the star had been formed. The elemental abundance pattern on the surface of stars suffering from such an external pollution should not be diverse but exhibit the average pattern of numerous supernovae. Future observations for a number of metal-deficient stars with [Fe/H]<-5 will be able to prove or disprove this external pollution scenario. Other possibilities to produce a star with this metallicity are also discussed.

We investigate dynamical response on size and velocity dispersion to mass loss by supernovae in formation of two-component spherical galaxies composed of baryon and dark matter. Three-dimensional deprojected de Vaucouleurs-like and exponential-like profiles for baryon, embedded in truncated singular isothermal and homogeneous profiles for dark matter, are considered. As a more realistic case, we also consider a dark matter profile proposed by Navarro, Frenk & White. For simplicity we assume that dark matter distribution is not affected by mass loss and that the change of baryonic matter distribution is homologous. We found that the degree of the response depends on the fraction of dark matter in the region where baryon is distributed, so that dwarf spheroidal galaxies would be affected even in a dark halo if they are formed by galaxy mergers in the envelope of the dark halo. Our results suggest that this scenario, combined with dynamical response, would make not only the observed trends but the dispersed characteristics of dwarf spheroidals.

Number counts of galaxies are re-analyzed using a semi-analytic model (SAM) of galaxy formation based on the hierarchical clustering scenario. Faint galaxies in the Subaru Deep Field (SDF) and the Hubble Deep Field (HDF) are compared with our model galaxies. We have determined the astrophysical parameters in the SAM that reproduce observations of nearby galaxies, and used them to predict the number counts and redshifts of faint galaxies for three cosmological models, the standard cold dark matter (CDM) universe, a flat lambda-CDM, and an open CDM. The novelty of our SAM analysis is the inclusion of selection effects arising from the cosmological dimming of surface brightness of high-z galaxies, and from the absorption of visible light by internal dust and intergalactic HI clouds. As was found in our previous work, in which the UV/optical HDF galaxies were compared with our model galaxies, we find that our SAM reproduces counts of near-IR SDF galaxies in low-density models, and that the standard CDM universe is not preferred, as suggested by other recent studies. Moreover, we find that simple prescriptions for (1) the timescale of star formation being proportional to the dynamical time scale of the formation of galactic disks, (2) the size of galactic disks being rotationally supported with the same specific angular momentum as that of surrounding dark halo, and (3) the dust optical depth being proportional to the metallicity of cold gas, cannot completely explain all of observed data. Improved prescriptions incorporating mild z-dependence for those are suggested from our SAM analysis.

Photometric and spectroscopic data of the energetic Type Ic supernova (SN) 2002ap are presented, and the properties of the SN are investigated through models of its spectral evolution and its light curve. The SN is spectroscopically similar to the "hypernova" SN 1997ef. However, its kinetic energy [$\sim (4-10) \times 10^{51}$ erg] and the mass ejected (2.5-5 $M_{\odot}$) are smaller, resulting in a faster-evolving light curve. The SN synthesized $\sim 0.07 M_{\odot}$ of $^{56}$Ni, and its peak luminosity was similar to that of normal SNe. Brightness alone should not be used to define a hypernova, whose defining character, namely very broad spectral features, is the result of a high kinetic energy. The likely main-sequence mass of the progenitor star was 20-25 $M_{\odot}$, which is also lower than that of both hypernovae SNe 1997ef and 1998bw. SN 2002ap appears to lie at the low-energy and low-mass end of the hypernova sequence as it is known so far. Observations of the nebular spectrum, which is expected to dominate by summer 2002, are necessary to confirm these values.

Variability of 226 AGNs in the near-infrared $J$, $H$, and $K'$ bands is analyzed and discussed. An ensemble average for measured variabilities was obtained for various samples of the AGNs divided by absolute $B$-magnitude $M_B$, redshift $z$, and radio strength. All the samples in the $J$, $H$, and $K'$ bands are found to give significant ensemble variability, but no significant wavelength dependence is found. The ensemble variability in the entire sample combining the $J$, $H$, and $K'$ samples is $\Delta m\approx 0.22$ mag, while $\Delta m\approx 0.18$ mag for the radio-quiet AGNs and $\Delta m\approx 0.26$ mag for radio-loud AGNs. The ensemble variability for the radio-quiet AGNs shows no significant $M_B$-dependence, while showing positive $M_B$-dependence for the radio-loud AGNs. In any samples the measured variability shows positive correlation among different passbands, with the correlation coefficients of $r_{JH}$, $r_{HK'}$, and $r_{JK'}$ ranging from 0.6 to 0.9. For radio-quiet AGNs, the coefficient $r_{HK'}$ in a redshift range of $0.1<z<0.3$ is significantly higher than $r_{JH}$ or $r_{JK'}$. The coefficient for the radio-loud AGNs with $0.6<z<1.0$ is as high as 0.95, irrespective of the passband. However, for the radio-quiet AGNs with $z>0.3$ and radio-loud AGNs with $z<0.3$, we cannot confirm such strong correlation among different passbands. All the features of near-infrared variability for the radio-quiet AGNs are consistent with a simple dust reverberation model of the central regions of AGNs. However, the features for the radio-loud AGNs are not fully explained by such a model, and a non-thermal variable component is suggested as a viable candidate for causing their large and fast variability in the near-infrared region.

Near-infrared $JHK'$ imaging photometry was obtained of 331 AGNs consisting mainly of Seyfert 1 AGNs and quasars (QSOs). This sample was selected to cover a range of radio emission strength, redshift from $z=0$ to 1, and absolute $B$-magnitude from $M_B=-29$ mag to -18 mag. Among low-$z$ AGNs with $z<0.3$, Seyfert $1-1.5$ AGNs are distributed over a region from a location typical of ``galaxies'' to a location typical of ``QSOs'' in the two-color $J-H$ to $H-K'$ diagram, but Seyfert $1.8-2$ AGNs are distributed around the location of ``galaxies''. Moreover, bright AGNs with respect to absolute $B$-magnitude are distributed near the location of ``QSOs'', while faint AGNs are near the location of ``galaxies''. The distribution of such low-$z$ AGNs in this diagram was found to have little dependence on their 6 cm radio flux. The near-infrared colors of the AGNs observed with an aperture of 7 pixels ($7.49''$) are more QSO-like than those observed with larger apertures up to 15 pixels ($16.1''$). This aperture effect may be explained by contamination from the light of host galaxies within larger apertures. This effect is more prominent for less luminous AGNs.

Observations of 226 AGNs in the near-infrared $J$, $H$, and $K'$ bands are presented along with the analysis of the observations for variability. Our sample consists mainly of Seyfert 1 AGNs and QSOs. About a quarter of the objects in each category are radio loud. The AGNs in the entire sample have the redshifts spanning the range from $z=0$ to 1, and the absolute magnitudes from $M_B=-29$ to -18. All the objects were observed twice and their variability was measured by differential photometry. A reduction method of differential photometry, optimized to the analysis of extended images, has been developed. The systematic error in variability arising from AGNs of highly extended images is estimated to be less than 0.01 mag in each of the $J$, $H$, and $K'$ bands. The systematic error arising from the flat fielding is negligible for most AGNs, although it is more than 0.1 mag for some particular cases. The overall average flat fielding error is 0.03 mag for the image pairs. We find that these systematic errors are superseded by statistical errors, and the overall average total systematic and statistical errors amounts to 0.05 mag in the measured variability in each band. We find that 58% of all the AGNs in the entire sample show variability of more than $2\sigma$, and 44% of more than $3\sigma$. This result holds independent of the $J$, $H$, and $K'$ bands. The detection rate of variability is higher for a subsample of higher photometric accuracy, and there appears no limit to this tendency. In particular, when we consider a subsample with small photometric errors of $\sigma<0.03$ mag, the rate of $2\sigma$ detection is 80%, and 64% for $3\sigma$ detection. This suggests that most AGNs are variable in the near-infrared.

We investigate the chemical evolution of odd-numbered elements such as sodium (Na) and aluminum (Al) during the early epochs of the Galactic halo with the use of a model that reproduces the observed box-shaped distribution of extremely metal-poor stars in the [Na, Al/Mg] versus [Mg/H] plane. Our model is constructed under the assumptions that those stars retain the elemental abundance patterns produced by individual Type II supernovae (SNe), and that the yields of the odd elements depend on the initial metallicity, z, of their SN progenitors. As a result, recent abundance determinations that clarify how the [Na, Al/Mg] ratios of field stars have evolved to the solar values enable us to deduce how the yields of these odd elements depend on z. The observed trends in these abundances, in particular the very large scatter (over 1 dex in [Al/Mg]) requires that the Al yield scales as m_Al proportional to z^0.6 for [Mg/H]<-1.8, while the observed [Na/Mg] trend requires that the Na yield scales as m_Na proportional to z^-0.4 for [Mg/H]<-1.8 and m_Na proportional to z^0.4 for [Mg/H]>-1.8. It is found that the predicted frequency distribution of stars in the [Na/Mg] versus [Mg/H] diagram is very sensitive to the assumed form of the primordial IMF, and that its slope is steeper than the Salpeter IMF. The necessity to match the observed abundance patterns of odd elements and the frequency distribution of extremely metal-poor stars should provide useful constraints on nucleosynthesis calculations of metal-free massive stars as well as on theories of their formation.

Number counts of galaxies are re-analyzed using a semi-analytic model (SAM) of galaxy formation based on the hierarchical clustering scenario. We have determined the astrophysical parameters in the SAM that reproduce observations of nearby galaxies, and used them to predict the number counts and redshifts of faint galaxies for three cosmological models for (1) the standard cold dark matter (CDM) universe, (2) a low-density flat universe with nonzero cosmological constant, and (3) a low-density open universe with zero cosmological constant. The novelty of our SAM analysis is the inclusion of selection effects arising from the cosmological dimming of surface brightness of high-redshift galaxies, and also from the absorption of visible light by internal dust and intergalactic \ion{H}{1} clouds. Contrary to previous SAM analyses which do not take into account such selection effects, we find, from comparison with observed counts and redshifts of faint galaxies in the Hubble Deep Field (HDF), that the standard CDM universe is {\it not} preferred, and a low-density universe either with or without cosmological constant is favorable, as suggested by other recent studies. Moreover, we find that a simple prescription for the time scale of star formation (SF), being proportional to the dynamical time scale of the formation of the galactic disk, is unable to reproduce the observed number- redshift relation for HDF galaxies, and that the SF time scale should be nearly independent of redshift, as suggested by other SAM analyses for the formation of quasars and the evolution of damped Ly-$\alpha$ systems.

We report observational analyses and theoretical interpretations of unusually red galaxies in the Subaru Deep Field (SDF). A careful analysis of the SDF data revealed a population with unusually red near-infrared (NIR) colors of J - K >~ 3-4, with higher confidence than the previous SDF result. Their surface number density drastically increases at K >~ 22 and becomes roughly the same with that of dusty starburst galaxies detected by submillimeter observations in recent years. These colors are even redder than the known population of the extremely red objects (EROs), and too red to explain by passively evolving elliptical galaxies which are the largest population of EROs. Hence these hyper extremely red objects (HEROs) should be considered as a distinct population from EROs. We discuss several possible interpretations of these enigmatic objects, and we show that these red NIR colors, K-band and sub-mm flux, and surface number density are quantitatively best explained by primordial elliptical galaxies reddened by dust, still in the starburst phase of their formation at z ~ 3.

Galaxy counts in the K band, (J-K)-colors, and apparent size distributions of faint galaxies in the Subaru Deep Field (SDF) down to K~24.5 were studied in detail. Special attention has been paid to take into account various selection effects including the cosmological dimming of surface brightness, to avoid any systematic bias which may be the origin of controversy in previously published results. We also tried to be very careful about systematic model uncertainties; we present a comprehensive surveys of these systematic uncertainties and dependence on various parameters. We found that the pure luminosity evolution (PLE) model is well consistent with all the SDF data down to K~22.5, without any evidence for number or size evolution in a low-density, Lambda-dominated flat universe which is now favored by various cosmological observations. If the popular Lambda-dominated universe is taken for granted, our result then gives a strong constraint on the number evolution of giant elliptical or early-type galaxies to z~1-2 which must be met by any models in the hierarchically clustering universe, since such galaxies are the dominant population in this magnitude range (K<~22.5). In the fainter magnitude range of K>~22.5, we found a slight excess of observed counts over the prediction of the PLE model when elliptical galaxies are treated as a single population. We suggest that this discrepancy reflects some number evolution of dwarf galaxies and/or the distinct populations of giant and dwarf elliptical galaxies which have been known for local elliptical galaxies.

We have developed a simulation code with the techniques which enhance both spatial and time resolution of the PM method for which the spatial resolution is restricted by the spacing of structured mesh. The adaptive mesh refinement (AMR) technique subdivides the cells which satisfy the refinement criterion recursively. The hierarchical meshes are maintained by the special data structure and are modified in accordance with the change of particle distribution. In general, as the resolution of the simulation increases, its time step must be shortened and more computational time is required to complete the simulation. Since the AMR enhances the spatial resolution locally, we reduce the time step locally also, instead of shortening it globally. For this purpose we used a technique of hierarchical time steps (HTS) which changes the time step, from particle to particle, depending on the size of the cell in which particles reside. Some test calculations show that our implementation of AMR and HTS is successful. We have performed cosmological simulation runs based on our code and found that many of halo objects have density profiles which are well fitted to the universal profile proposed by Navarro, Frenk, & White (1996) over the entire range of their radius.

Deep optical and near-infrared galaxy counts are utilized to estimate the extragalactic background light (EBL) coming from normal galactic light in the universe. Although the slope of number-magnitude relation of the faintest counts is flat enough for the count integration to converge, considerable fraction of EBL from galaxies could still have been missed in deep galaxy surveys because of various selection effects including the cosmological dimming of surface brightness of galaxies. Here we give an estimate of EBL from galaxy counts, in which these selection effects are quantitatively taken into account for the first time, based on reasonable models of galaxy evolution which are consistent with all available data of galaxy counts, size, and redshift distributions. We show that the EBL from galaxies is best resolved into discrete galaxies in the near-infrared bands (J, K) by using the latest data of the Subaru Deep Field; more than 80-90% of EBL from galaxies has been resolved in these bands. Our result indicates that the contribution by missing galaxies cannot account for the discrepancy between the count integration and recent tentative detections of diffuse EBL in the K-band (2.2 micron), and there may be a very diffuse component of EBL which has left no imprints in known galaxy populations.

We present predictions of the evolution of the light elements, Li, Be, and B, in the early epochs of the Galactic halo, using a model of supernova-induced chemical evolution based on contributions from supernovae (SNe) and cosmic rays (CRs). The most important prediction from our model is that the abundances of light elements in extremely metal-poor stars might be used as age indicators in the very early stages of an evolving halo population, at times when metallicity in most stars are dominated by local metal enrichment due to nearby SN events, and is poorly correlated with age. Our results show that the best ``cosmic clock'' is the $^6$Li abundance. We have derived relationships among various cosmic-ray parameters. Although our model is successful for certain sets of cosmic-ray parameters, larger energy should be absorbed by energetic particles from each SN than required to the current situation of Galactic disk. We discuss an alternative hypothesis of AGN activity in the early Galaxy as another possible accelerator of CRs.

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

(Abridged) We present a detailed analysis of the number count and photometric redshift distribution of faint galaxies in the Hubble Deep Field (HDF), paying a special attention to the selection effects including the cosmological dimming of surface brightness of galaxies. We find a considerably different result from previous studies ignoring the selection effects, and these effects should therefore be taken into account in the analysis. We find that the model of pure luminosity evolution (PLE) of galaxies in the Einstein-de Sitter (EdS) universe predicts much smaller counts than those observed at faint magnitude limits by a factor of more than 10, so that a very strong number evolution of galaxies with \eta > 3-4 must be invoked to reproduce the I_{814} counts, when parametrized as \phi^* \propto (1+z)^\eta. However we show that such a strong number evolution under realistic merging processes of galaxies can not explain the steep slope of the B_{450} and V_{606} counts, and it is seriously inconsistent with their photometric redshift distribution. We find that these difficulties still persist in an open universe with \Omega_0 > 0.2, but are resolved only when we invoke a $\Lambda$-dominated flat universe, after examining various systematic uncertainties in modeling the formation and evolution of galaxies. The present analysis revitalizes the practice of using faint number counts as an important cosmological test, giving one of the arguments against the EdS universe and suggests acceleration of the cosmic expansion by vacuum energy density. While a modest number evolution of galaxies with \eta ~ 1 is still necessary even in a Lambda-dominated universe, a stronger number evolution with \eta > 1 is rejected from the HDF data, giving a strong constraint on the merger history of galaxies.

We have refined the estimate of the primordial level of $^7$Li abundance to an accuracy better than 10%, based on high-precision Li abundances for metal-poor halo stars, and a recent model of post-BBN (Big Bang Nucleosynthesis) chemical evolution that provides a quantitative explanation of the detected gentle ascent of the Spite Plateau for stars with metallicities [Fe/H]$ > -3$. Our maximum likelihood analysis obtains an estimate for the primordial Li abundance of $A({\rm Li})_p=2.07^{+0.16}_{-0.04}$, after taking into account possible systematic errors in the estimation of Li abundances, with the exception of a still-controversial issue regarding stellar depletion. The inferred value of $\eta$ (the baryon-to-photon number-density ratio in the universe) based on this estimate is more consistent with that derived from the set of reported ``low He'' + ``high D'' from extragalactic sites than that derived from reported ``high He'' + ``low D'' measurements. Since, within current models of stellar depletion processes, it is difficult to account for the observed very small scatter of Li abundance in metal-poor stars, our estimate of $A({\rm Li})_p$ should be taken as an independent constraint on the baryonic mass density parameter in the universe, giving $\Omega_b h^2=(0.64-1.4)\times 10^{-2}$ with $h=H_0/100$ km s$^{-1}$Mpc$^{-1}$.

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.

Recent detailed simulations of Galactic Chemical Evolution have shown that the heavy elements, in particular [Fe/H], are expected to exhibit a weak, or absent, correlation with stellar ages in the early Galaxy due to the lack of efficient mixing of interstellar material enriched by individual Type II supernovae. A promising alternative "chronometer" of stellar ages is suggested, based on the expectation that the light elements Be and B are formed primarily as spallation products of Galactic Cosmic Rays.

We suggest that if the astrophysical site for r-process nucleosynthesis in the early Galaxy is confined to a narrow mass range of Type II supernova (SN II) progenitors, with a lower mass limit of Mms = 20 Msun, a unique feature in the observed distribution of [Ba/Mg] vs.[Mg/H] for extremely metal-poor stars can be adequately reproduced. We associate this feature, a bifurcation of the observed elemental ratios into two branches in the Mg abundance interval -2.7 < [Mg/H] < -2.3, with two distinct processes. The first branch, which we call the ``y''-branch, is associated with the production of Ba and Mg from individual massive supernovae. We conclude that SNe II with Mms = 20 Msun are the dominant source of r-process nucleosynthesis in the early Galaxy. An SN-induced chemical evolution model with this Mms-dependent Ba yield creates the y-branch, reflecting the different nucleosynthesis yields of [Ba/Mg] for each SN II with Mms > 20 Msun. The second branch, which we call the ``i''-branch, is associated with the elemental abundance ratios of stars which were formed in the dense shells of the interstellar medium swept up by SNe II with Mms < 20 Msun that do not synthesize r-process elements, and applies to stars with observed Mg abundances in the range [Mg/H] < -2.7. The Ba abundances in these stars reflect those of the interstellar gas at the (later) time of their formation. The existence of a [Ba/Mg] i-branch strongly suggests that SNe II which are associated with stars of progenitor mass Mms < 20 Msun are infertile sources for the production of r-process elements. We predict the existence of this i-branch for other r-process elements, such as europium (Eu), to the extent that their production site is in common with Ba.

A model of supernova-driven chemical evolution of the Galactic halo, recently proposed by Tsujimoto, Shigeyama, & Yoshii (1999, ApJL, 519, 64), is extended in order to investigate the evolution of light elements such as Be and B (BeB), which are produced mainly through spallative reactions with Galactic cosmic rays. In this model each supernova sweeps up the surrounding interstellar gas into a dense shell and directly enriches it with ejecta which consist of heavy elements produced in each Type II supernova with different progenitor masses. We propose a two-component source for GCRs such that both interstellar gas and fresh SN ejecta engulfed in the shell are accelerated by the shock wave. Our model results include: (1) a prediction of the intrinsic scatter in BeB and [Fe/H] abundances within the model, (2) a successful prediction of the observed linear trend between BeB and [Fe/H], (3) a proposal for using BeB as a cosmic clock, as an alternative to [Fe/H], and (4) a method for possibly constraining the BBN model from future observations of metal-poor stars.

A model for Galactic chemical evolution, driven by supernova-induced star formation, is formulated and used to examine the nature of the Galactic halo at early epochs. In this model, new stars are formed following each supernova event, thus their abundance pattern is determined by the combination of heavy elements ejected from the supernova itself and those elements which are already present in the interstellar gas swept up by the supernova remnant. The end result is a prediction of large scatter in the abundance ratios among low-metallicity stars, reflecting a different nucleosynthesis yield for each Type II supernova with a different progenitor mass. Formation of new stars is terminated when supernova remnants sweep up too little gas to form shells. We show from calculations based on the above scenario that (i) the observed [Fe/H] distribution for the Galactic halo field stars can be reproduced without effectively decreasing the heavy-element yields from Type II supernovae by some manipulation required by previous models (e.g., via mass loss from the early Galaxy, or later mixing with ``pristine'' hydrogen clouds), (ii) the large observed scatter in the abundance ratio [Eu/Fe] for the most metal-poor stars can also be reproduced, and (iii) the frequency distribution of stars in the [Eu/Fe]-[Fe/H] plane can be predicted. Our model suggests that the probability of identifying essentially metal-free stars (Population III) in the local halo is around one in 10^{3-4}, provided that star formation in the halo is confined to individual gas clouds with mass of 10^{6-7} Msun and that the initial mass function of metal-free stars is not significantly different from the Salpeter mass function.

We have developed a new parallel tree method which will be called the forest method hereafter. This new method uses the sectional Voronoi tessellation (SVT) for the domain decomposition. The SVT decomposes a whole space into polyhedra and allows their flat borders to move by assigning different weights. The forest method determines these weights based on the load balancing among processors by means of the over-load diffusion (OLD). Moreover, since all the borders are flat, before receiving the data from other processors, each processor can collect enough data to calculate the gravity force with precision. Both the SVT and the OLD are coded in a highly vectorizable manner to accommodate on vector parallel processors. The parallel code based on the forest method with the Message Passing Interface is run on various platforms so that a wide portability is guaranteed. Extensive calculations with 15 processors of Fujitsu VPP300/16R indicate that the code can calculate the gravity force exerted on 10^5 particles in each second for some ideal dark halo. This code is found to enable an N-body simulation with 10^7 or more particles for a wide dynamic range and is therefore a very powerful tool for the study of galaxy formation and large-scale structure in the universe.

We simulate the dynamical, chemical, and spectro-photometric evolution of dwarf and normal elliptical galaxies embedded in a dark matter halo, using a three-dimensional N-body/SPH simulation code. For the forming dwarf elliptical galaxies, supernova-driven winds propagating outwards from inside the system collide with the infalling gas and produce the super-shell in which stars are formed. The resulting stellar system forms a loosely bound virialized system due to the significant mass loss and has a large velocity dispersion and a large core. Consequently the surface brightness distribution shows an exponential profile and the color distribution shows a positive gradient such that the colors become redder away from the galaxy center in agreement with observations. On the other hand, for the normal elliptical galaxies with deep gravitational potential, the mass loss out of the system does not have a significant dynamical effect on the stellar system. The resulting surface brightness distribution has a large central concentration like de Vaucouleurs' $r^{1/4}$-profile and the color distribution shows a negative gradient as observed. Our simulation shows that different features between dwarf and normal elliptical galaxies stem from different cooling efficiencies for their respective protogalaxies in a standard CDM universe.

The observed Fe II(UV+optical)/Mg II lambda lambda 2796,2804 flux ratio from a gravitationally lensed quasar B1422+231 at z=3.62 is interpreted in terms of detailed modeling of photoionization and chemical enrichment in the broad-line region (BLR) of the host galaxy. The delayed iron enrichment by Type Ia supernovae is used as a cosmic clock. Our standard model, which matches the Fe II/Mg II ratio, requires the age of 1.5 Gyr for B1422+231 with a lower bound of 1.3 Gyr, which exceeds the expansion age of the Einstein-de Sitter Omega_0=1 universe at a redshift of 3.62 for any value of the Hubble constant in the currently accepted range, H_0=60-80 km,s^{-1},Mpc^{-1}. This problem of an age discrepancy at z=3.62 can be unraveled in a low-density Omega_0<0.2 universe, either with or without a cosmological constant, depending on the allowable redshift range of galaxy formation. However, whether the cosmological constant is a required option in modern cosmology awaits a thorough understanding of line transfer processes in the BLRs.

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.

We have performed a detailed V/Vmax test for a sample of the Canada-France Redshift Survey (CFRS) for the purpose of examining whether the comoving number density of field galaxies changes significantly at redshifts of z<1. Taking into account the luminosity evolution of galaxies which depends on their morphological type through different history of star formation, we obtain <V/Vmax> \sim 0.5 in the range of 0.3<z<0.8, where reliable redshifts were secured by spectroscopy of either absorption or emission lines for the CFRS sample. This indicates that a picture of mild evolution of field galaxies without significant mergers is consistent with the CFRS data. Early-type galaxies, selected by their (V-I)_{AB} color, become unnaturally deficient in number at z>0.8 due to the selection bias, thereby causing a fictitious decrease of <V/Vmax>. We therefore conclude that a reasonable choice of upper bound of redshift z \sim 0.8 in the V/Vmax test saves the picture of passive evolution for field ellipticals in the CFRS sample, which was rejected by Kauffman, Charlot, & White (1996) without confining the redshift range. However, about 10% of the CFRS sample consists of galaxies having colors much bluer than predicted for irregular galaxies, and their \avmax is significantly larger than 0.5. We discuss this population of extremely blue galaxies in terms of starburst that has just turned on at their observed redshifts.

The present determination of the absolute magnitude $M_V(RR)$ of RR Lyrae stars is twofold, relying upon Hipparcos proper motions and trigonometric parallaxes separately. First, applying the statistical parallax method to the proper motions, we find $<M_V(RR)>=0.69\pm0.10$ for 99 halo RR Lyraes with $<$[Fe/H]$>$ =--1.58. Second, applying the Lutz-Kelker correction to the RR Lyrae HIP95497 with the most accurately measured parallax, we obtain $M_V(RR)$=(0.58--0.68)$^{+0.28}_{-0.31}$ at [Fe/H]=--1.6. Furthermore, allowing full use of low accuracy and negative parallaxes as well for 125 RR Lyraes with -- 2.49$\leq$[Fe/H]$\leq$0.07, the maximum likelihood estimation yields the relation, $M_V(RR)$=(0.59$\pm$0.37)+(0.20$\pm$0.63)([Fe/H]+1.60), which formally agrees with the recent preferred relation. The same estimation yields again $<M_V(RR)>$ = $0.65\pm0.33$ for the 99 halo RR Lyraes. Although the formal errors in the latter three parallax estimates are rather large, all of the four results suggest the fainter absolute magnitude, $M_V(RR)$$\approx$0.6--0.7 at [Fe/H]=--1.6. The present results still provide the lower limit on the age of the universe which is inconsistent with a flat, matter-dominated universe and current estimates of the Hubble constant.

We present new deep $K'$ counts of stars at the South Galactic Pole (SGP) taken with the NAOJ PICNIC camera to $K'=17.25$. Star-galaxy separation to $K'=17.5$ was accomplished effectively using image profiles because the pixel size we used is 0.509 arcsec. We interpret these counts using the SKY (Cohen 1994) model of the Galactic point source sky and determine the relative normalization of halo-to-disk populations, and the location of the Sun relative to the Galactic plane, within the context of this model. The observed star counts constrain these parameters to be: halo/disk $\sim$ 1/900 and z$_\odot$=16.5$\pm$2.5 pc. These values have been used to correct our SGP galaxy counts for contamination by the point source Galactic foreground.

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 present new K-band galaxy number counts from K=13 to 20.5 obtained from $K'$-band surveys in the south galactic pole region, which cover 180.8 arcmin$^2$ to a limiting magnitude of K=19, and 2.21 arcmin$^2$ to K=21. These are currently the most precise K-band galaxy counts at $17.5<K<19.0$ because the area of coverage is largest among the existing surveys for this magnitude range. The completeness and photometry corrections are estimated from the recovery of simulated galaxy and stellar profiles added to the obtained field image. Many simulations were carried out to construct a probability matrix which corrects the galaxy counts at the faint-end magnitudes of the surveys so the corrected counts can be compared with other observations. The K-band star counts in the south galactic pole region to $K=17.25$ are also presented for use to constrain the vertical structure of the Galaxy.

New X-ray observations of the sparse NGC 3258 group of galaxies made by the ASCA satellite with good spectral and spatial resolution has revealed that this group has a gravitational potential deep enough to prohibit significant mass removal from the system. The baryonic fraction within 240 kpc is found to be 0.065 +0.051/-0.020 for h_{50}=1, where h_{50}=H_0/50 km/s/Mpc, in good agreement with the universal value of 0.05 +/-0.01 predicted by standard Big Bang nucleosynthesis for a Universe with Omega_0=1 and h_{50}=1. Since the deep potential of the NGC 3258 group ensures that all pristine intragroup gas has been retained, the baryonic fraction of the NGC 3258 group is indicative of the universal value. Consequently it seems premature to rule out a critical Universe.

We show that evolution of the luminosity density of galaxies in the universe provides a powerful test for the geometry of the universe. Using reasonable galaxy evolution models of population synthesis which reproduce the colors of local galaxies of various morphological types, we have calculated the luminosity density of galaxies as a function of redshift $z$. Comparison of the result with recent measurements by the Canada-France Redshift Survey in three wavebands of 2800{\AA}, 4400{\AA}, and 1 micron at z<1 indicates that the \Lambda-dominated flat universe with \lambda_0 \sim 0.8 is favored, and the lower limit on \lambda_0 yields 0.37 (99% C.L.) or 0.53 (95% C.L.) if \Omega_0+\lambda_0=1. The Einstein-de Sitter universe with (\Omega_0, \lambda_0)=(1, 0) and the low-density open universe with (0.2, 0) are however ruled out with 99.86% C.L. and 98.6% C.L., respectively. The confidence levels quoted apply unless the standard assumptions on galaxy evolution are drastically violated. We have also calculated a global star formation rate in the universe to be compared with the observed rate beyond z \sim 2. We find from this comparison that spiral galaxies are formed from material accretion over an extended period of a few Gyrs, while elliptical galaxies are formed from initial star burst at z >~ 5 supplying enough amount of metals and ionizing photons in the intergalactic medium.

Oxygen to iron abundance ratios of metal-poor stars provide information on nucleosynthesis yields from massive stars which end in Type II supernova explosions. Using a standard model of chemical evolution of the Galaxy we have reproduced the solar neighborhood abundance data and estimated the oxygen and iron yields of genuine SN II origin. The estimated yields are compared with the theoretical yields to derive the relation between the lower and upper mass limits in each generation of stars and the IMF slope. Independently of this relation, we furthermore derive the relation between the lower mass limit and the IMF slope from the stellar mass to light ratio in the solar neighborhood. These independent relations unambiguously determine the upper mass limit of $m_u=50 \pm 10 M_sun$ and the IMF slope index of 1.3 - 1.6 above 1 M_sun. This upper mass limit corresponds to the mass beyond which stars end as black holes without ejecting processed matter into the interstellar medium. We also find that the IMF slope index below 0.5 M_sun cannot be much shallower than 0.8.

We simulate the dynamical and chemical evolution of a dwarf galaxy embedded in a dark matter halo, using a three-dimensional N-body/SPH simulation code combined with stellar population synthesis. The initial condition is adopted in accord with a $10^{10}M_\odot$ virialized sphere in a 1-sigma CDM perturbation which contains 10% baryonic mass. A supersonic spherical outflow is driven by the first star burst near the center of the galaxy and produces an expanding super shell in which stars are subsequently formed. Consecutive formation of stars in the expanding shell makes the stellar system settled with the exponential brightness profile, the positive metallicity gradient, and the inverse color gradient in agreement with observed features of dwarf galaxies. We therefore propose that the energy feedback via stellar winds and supernovae is a decisive mechanism for formation of less compact, small systems like dwarf galaxies.

We combine deep K-band (Keck) with V- and I-band (NTT) observations of two high-Galactic latitude fields, surveying a total of ~2 sq. arcmin. The K-band galaxy counts continue to rise above K=22, reaching surface densities of few x 10^5 per sq. degree. The slope for the counts is (d log(N) per mag per sq. degree) = 0.23 +/- 0.02 between K=18-23, consistent with other deep K surveys. The numbers of galaxies in each mag bin is about two times greater than the galaxy counts of Djorgovski et al. (1995). The optical and near infrared magnitudes of all objects detected in the V+I+K image are discussed in the context of grids of isochrone synthesis galaxy evolutionary models (Bruzual & Charlot 1993, 1995). The colors of most of the observed galaxies are consistent with a population drawn from a broad redshift distribution. A few galaxies at K=19-20 are red in both colors (V-I>3; I-K>2), consistent with being early-type galaxies having undergone a burst of star formation at z>5 and viewed at z~1. At K>20, we find ~8 ``red outlier'' galaxies with I-K>4 and V-I<2.5, whose colors are difficult to mimic by a single evolving or non-evolving stellar population at any redshift. They are likely either low-metallicity, dusty dwarf galaxies, or old galaxies at high redshift (z>1.2). Their surface density is several per square arcminute, which is so high that they are probably common objects of low luminosity $L<L_*$.

A new solution is presented for the puzzling, observed universality of the exponential luminosity profiles, perpendicular to the disk plane of spiral and lenticular galaxies. It is shown that such exponential $z$-profiles result naturally from gaseous protodisks which settle into isothermal equilibrium prior to star formation. Subsequent cooling leads to a gravitational contraction of the gas towards the equatorial plane and to a stellar exponential $z$-profile if the star formation rate is assumed to be comparable to the cooling rate of the gas. The final stellar scale height depends only on the initial gas temperature and local surface density. This model therefore provides a new method to investigate the early energetic state of galactic protodisks with measured scale heights and surface densities along the disk plane.

We report the discovery of a group of galaxies at redshift 2.38. We imaged about 10% of a claimed supercluster of QSO absorption-lines at z=2.38 (Francis & Hewett 1993). In this small field (2 arcmin radius) we detect two Ly-alpha emitting galaxies. The discovery of two such galaxies in our tiny field supports Francis & Hewett's interpretation of the absorption-line supercluster as a high redshift "Great Wall". One of the Ly-alpha galaxies lies 22 arcsec from a background QSO, and may be associated with a multi-component Ly-alpha absorption complex seen in the QSO spectrum. This galaxy has an extended (50kpc) lumpy Ly-alpha morphology, surrounding a compact IR-bright nucleus. The nucleus shows a pronounced break in its optical-UV colors at about 4000 A (rest-frame), consistent with a stellar population of mass about 7E11 solar masses, an age of more than 500 Myr, and little on-going star-formation. C IV emission is detected, suggesting that a concealed AGN is present. Extended H-alpha emission is also detected; the ratio of Ly-alpha flux to H-alpha is abnormally low (about 0.7), probable evidence for extended dust. This galaxy is surrounded by a number of very red (B-K>5) objects, some of which have colors suggesting that they too are at z=2.38. We hypothesize that this galaxy, its neighbors and a surrounding lumpy gas cloud may be a giant elliptical galaxy in the act of bottom-up formation.

The spectrum of the supernova relic neutrino background (SRN) from collapse-driven supernovae ever occurred in the universe is calculated by using a realistic, time-dependent supernova rate derived from a standard model of galaxy evolution based on the population synthesis method. The SRN spectrum we show here is the most realistic at present, because the largest uncertainty in previous theoretical predictions has come from unrealistic assumptions of the supernova rate so far made. The SRN is one of the targets of the Superkamiokande (SK) detector which will be constructed in a year and the SRN, if at all detected, would provide a new tool to probe the history of supernova explosions in the universe. The expected event rate at the SK is therefore calculated in this paper. Our major results include: (1) the supernova rate is much higher in the early phase of evolution of galaxies and there appears a hump in the SRN spectrum in the low-energy region of $\ltilde 5$ MeV, (2) the SRN flux depends on the Hubble constant ($H_0$) in a way approximately proportional to $H_0^2$ and only weakly on the density parameter of the universe ($\Omega_0$) and a cosmological constant ($\lambda_0$), and (3) the plausible event rate at the SK is 1.2 yr$^{-1}$ in the observable energy range. Such a low event rate is due mainly to a quite low supernova rate at present which is averaged over the morphological types of galaxies. The most optimistic rate in our model is found to be 4.7 yr$^{-1}$, and if more events are detected, we will have to reconsider our current understanding of collapse-driven supernovae and evolution of galaxies.

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}$.