The unique combination of superb spatial resolution, wide-area coverage, and deep depth of the optical imaging from the Hyper Suprime-Cam (HSC) Subaru Strategic Program is utilized to search for dual quasar candidates. Using an automated image analysis routine on 34,476 known SDSS quasars, we identify those with two (or more) distinct optical point sources in HSC images covering 796 deg^2. We find 421 candidates out to a redshift of 4.5 of which one hundred or so are more likely after filtering out contaminating stars. Angular separations of 0.6 - 4.0" correspond to projected separations of 3 - 30 kpc, a range relatively unexplored for population studies of luminous dual quasars. Using Keck-I/LRIS and Gemini-N/NIFS, we spectroscopically confirm three dual quasar systems at z < 1, two of which are previously unknown out of eight observed, based on the presence of characteristic broad emission lines in each component, while highlighting that the continuum of one object in one of the pairs is reddened. In all cases, the [OIII]5007 emission lines have mild velocity offsets, thus the joint [OIII] line profile is not double-peaked. We find a dual quasar fraction of 0.26+/-0.18% and no evidence for evolution. A comparison with the Horizon-AGN simulation seems to support the case of no evolution in the dual quasar fraction when broadly matching the quasar selection. These results may indicate a scenario in which the frequency of the simultaneous triggering of luminous quasars is not as sensitive as expected to the cosmic evolution of the merger rate or gas content of galaxies.

Gravitationally lensed quasars are useful for studying astrophysics and cosmology, and enlarging the sample size of lensed quasars is important for multiple studies. In this work, we develop a lens search algorithm for four-image (quad) lensed quasars based on their time variability. In the development of the lens search algorithm, we constructed a pipeline simulating multi-epoch images of lensed quasars in cadenced surveys, accounting for quasar variabilities, quasar hosts, lens galaxies, and the PSF variation. Applying the simulation pipeline to the Hyper Suprime-Cam (HSC) transient survey, we generated HSC-like difference images of the mock lensed quasars from Oguri & Marshall's lens catalog. We further developed a lens search algorithm that picks out variable objects as lensed quasar candidates based on their spatial extent in the difference images. We tested our lens search algorithm with the mock lensed quasars and variable objects from the HSC transient survey. Using difference images from multiple epochs, our lens search algorithm achieves a high true-positive rate (TPR) of 90.1% and a low false-positive rate (FPR) of 2.3% for the bright quads with wide separation. With a preselection of the number of blobs in the difference image, we obtain a TPR of 97.6% and a FPR of 2.6% for the bright quads with wide separation. Even when difference images are only available in one single epoch, our lens search algorithm can still detect the bright quads with wide separation at high TPR of 97.6% and low FPR of 2.4% in the optimal seeing scenario, and at TPR of $\sim94%$ and FPR of $\sim5%$ in typical scenarios. Therefore, our lens search algorithm is promising and is applicable to ongoing and upcoming cadenced surveys, particularly the HSC transient survey and the Rubin Observatory Legacy Survey of Space and Time, for finding new lensed quasar systems. [abridged]

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

We study variability of active galactic nuclei (AGNs) by using the deep optical multiband photometry data obtained from the Hyper Suprime-Cam Subaru Strategic Program (HSC SSP) survey in the COSMOS field. The images analyzed here were taken with 8, 10, 13, and 15 epochs over three years in the $g$, $r$, $i$, and $z$ bands, respectively. We identified 491 robust variable AGN candidates, down to $i=25$ mag and with redshift up to $4.26$. Ninety percent of the variability-selected AGNs are individually identified with the X-ray sources detected in the Chandra COSMOS Legacy survey. We investigate their properties in variability by using structure function analysis and find that the structure function for low-luminosity AGNs ($L_{\mathrm{bol}}\lesssim10^{45}$ erg s$^{-1}$) shows a positive correlation with luminosity, which is the opposite trend for the luminous quasars. This trend is likely to be caused by larger contribution of the host galaxy light for lower-luminosity AGNs. Using the model templates of galaxy spectra, we evaluate the amount of host galaxy contribution to the structure function analysis and find that dominance of the young stellar population is needed to explain the observed luminosity dependence. This suggests that low-luminosity AGNs at $0.8\lesssim z\lesssim1.8$ are predominantly hosted in star-forming galaxies. The X-ray stacking analysis reveals the significant emission from the individually X-ray undetected AGNs in our variability-selected sample. The stacked samples show very large hardness ratios in their stacked X-ray spectrum, which suggests that these optically variable sources have large soft X-ray absorption by dust-free gas.

With a booming number of Type Ia supernovae (SNe Ia) discovered within a few days of their explosions, a fraction of SNe Ia that show luminosity excess in the early phase (early-excess SNe Ia) have been confirmed. In this article, we report early-phase observations of seven photometrically normal SNe Ia (six early detections and one deep non-detection limit) at the COSMOS field through a half-year transient survey as a part of the Hyper Suprime-Cam Subaru Strategic Program (HSC SSP). In particular, a blue light-curve excess was discovered for HSC17bmhk, a normal SN Ia with rise time longer than 18.8 days, during the first four days after the discovery. The blue early excess in optical wavelength can be explained not only by interactions with a non-degenerate companion or surrounding dense circumstellar matter but also radiation powered by radioactive decays of $^{56}$Ni at the surface of the SN ejecta. Given the growing evidence of the early-excess discoveries in normal SNe Ia that have longer rise times than the average and a similarity in the nature of the blue excess to a luminous SN Ia subclass, we infer that early excess discovered in HSC17bmhk and other normal SNe Ia are most likely attributed to radioactive $^{56}$Ni decay at the surface of the SN ejecta. In order to successfully identify normal SNe Ia with early excess similar to that of HSC17bmhk, early UV photometries or high-cadence blue-band surveys are necessary.

Rapidly evolving transients form a new class of transients which show shorter timescales of the light curves than those of typical core-collapse and thermonuclear supernovae. We performed a systematic search for rapidly evolving transients using the deep data taken with the Hyper Suprime-Cam Subaru Strategic Program Transient Survey. By measuring the timescales of the light curves of 1824 transients, we identified 5 rapidly evolving transients. Our samples are found in a wide range of redshifts (0.3 $\le$ z $\le$ 1.5) and peak absolute magnitudes ($-$17 $\ge$ $M_i$ $\ge$ $-$20). The properties of the light curves are similar to those of the previously discovered rapidly evolving transients. They show a relatively blue spectral energy distribution, with the best-fit blackbody of 8,000 - 18,000 K. We show that some of the transients require power sources other than the radioactive decays of $^{56}$Ni because of their high peak luminosities and short timescales. The host galaxies of all the samples are star-forming galaxies, suggesting a massive star origin for the rapidly evolving transients. The event rate is roughly estimated to be $\sim$4,000 events yr$^{-1}$ Gpc$^{-3}$, which is about 1 $\%$ of core-collapse supernovae.

We perform a high-cadence transient survey with Subaru Hyper Suprime-Cam (HSC), which we call the Subaru HSC survey Optimized for Optical Transients (SHOOT). We conduct HSC imaging observations with time intervals of about one hour on two successive nights, and spectroscopic and photometric follow-up observations. A rapidly declining blue transient SHOOT14di at $z=0.4229$ is found in observations on two successive nights with an image subtraction technique. The rate of brightness change is $+1.28^{+0.40}_{-0.27}~{\rm mag~day^{-1}}$ ($+1.83^{+0.57}_{-0.39}~{\rm mag~day^{-1}}$) in the observer (rest) frame and the rest-frame color between $3400$ and $4400~\unicode[.8,0]{x212B}$ is $M_{\rm 3400\unicode[.8,0]{x212B}}-M_{\rm 4400\unicode[.8,0]{x212B}}=-0.4$. The nature of the object is investigated by comparing its peak luminosity, decline rate, and color with those of transients and variables previously observed, and those of theoretical models. None of the transients or variables share the same properties as SHOOT14di. Comparisons with theoretical models demonstrate that, while the emission from the cooling envelope of a Type IIb supernova shows a slower decline rate than SHOOT14di, and the explosion of a red supergiant star with a dense circumstellar wind shows a redder color than SHOOT14di, the shock breakout at the stellar surface of the explosion of a $25M_{\odot}$ red supergiant star with a small explosion energy of $\leq0.4\times10^{51}$ erg reproduces the multicolor light curve of SHOOT14di. This discovery shows that a high-cadence, multicolor optical transient survey at intervals of about one hour, and continuous and immediate follow-up observations, is important for studies of normal core-collapse supernovae at high redshifts.

We report our observations of HSC16aayt (SN 2016jiu), which was discovered by the Subaru/Hyper Suprime-Cam (HSC) transient survey conducted as part of Subaru Strategic Program (SSP). It shows very slow photometric evolution and its rise time is more than 100 days. The optical magnitude change in 400 days remains within 0.6 mag. Spectra of HSC16aayt show a strong narrow emission line and we classify it as a Type IIn supernova. The redshift of HSC16aayt is 0.6814 +/- 0.0002 from the spectra. Its host galaxy center is at 5 kpc from the supernova location and HSC16aayt might be another example of isolated Type IIn supernovae, although the possible existence of underlying star forming activity of the host galaxy at the supernova location is not excluded.

We study the faint stellar halo of isolated central galaxies, by stacking galaxy images in the HSC survey and accounting for the residual sky background sampled with random points. The surface brightness profiles in HSC $r$-band are measured for a wide range of galaxy stellar masses ($9.2<\log_{10}M_\ast/M_\odot<11.4$) and out to 120 kpc. Failing to account for the stellar halo below the noise level of individual images will lead to underestimates of the total luminosity by $\leq 15\%$. Splitting galaxies according to the concentration parameter of their light distributions, we find that the surface brightness profiles of low concentration galaxies drop faster between 20 and 100 kpc than those of high concentration galaxies. Albeit the large galaxy-to-galaxy scatter, we find a strong self-similarity of the stellar halo profiles. They show unified forms once the projected distance is scaled by the halo virial radius. The colour of galaxies is redder in the centre and bluer outside, with high concentration galaxies having redder and more flattened colour profiles. There are indications of a colour minimum, beyond which the colour of the outer stellar halo turns red again. This colour minimum, however, is very sensitive to the completeness in masking satellite galaxies. We also examine the effect of the extended PSF in the measurement of the stellar halo, which is particularly important for low mass or low concentration galaxies. The PSF-corrected surface brightness profile can be measured down to $\sim$31 $\mathrm{mag}/\mathrm{arcsec}^2$ at 3-$\sigma$ significance. PSF also slightly flattens the measured colour profiles.

We present an overview of a deep transient survey of the COSMOS field with the Subaru Hyper Suprime-Cam (HSC). The survey was performed for the 1.77 deg$^2$ ultra-deep layer and 5.78 deg$^2$ deep layer in the Subaru Strategic Program over 6- and 4-month periods from 2016 to 2017, respectively. The ultra-deep layer shows a median depth per epoch of 26.4, 26.3, 26.0, 25.6, and 24.6 mag in $g$, $r$, $i$, $z$, and $y$ bands, respectively; the deep layer is $\sim0.6$ mag shallower. In total, 1,824 supernova candidates were identified. Based on light curve fitting and derived light curve shape parameter, we classified 433 objects as Type Ia supernovae (SNe); among these candidates, 129 objects have spectroscopic or COSMOS2015 photometric redshifts and 58 objects are located at $z > 1$. Our unique dataset doubles the number of Type Ia SNe at $z > 1$ and enables various time-domain analyses of Type II SNe, high redshift superluminous SNe, variable stars, and active galactic nuclei.

We present Keck spectroscopic observations of three probable high redshift superluminous supernovae (SLSNe) from the Subaru HIgh-Z sUpernova CAmpaign (SHIZUCA), confirming redshifts of 1.851, 1.965 and 2.399. The host galaxies were selected for transient monitoring from multi-band photometric redshifts. The supernovae are detected during their rise, and the classically scheduled spectra are collected near maximum light. The restframe far-ultraviolet (FUV;$\sim$1000A-2500A) spectra include a significant host galaxy flux contribution and we compare our host galaxy subtracted spectra to UV-luminous SNe from the literature. While the signal-to-noise ratios of the spectra presented here are sufficient for redshift confirmation, supernova spectroscopic type confirmation remains inconclusive. The success of the first SHIZUCA Keck spectroscopic follow-up program demonstrates that campaigns such as SHIZUCA are capable of identifying high redshift SLSNe with sufficient accuracy, speed and depth for rapid, well-cadenced and informative follow-up.

We report our first discoveries of high-redshift supernovae from the Subaru HIgh-Z sUpernova CAmpaign (SHIZUCA), the transient survey using Subaru/Hyper Suprime-Cam. We report the discovery of three supernovae at the spectroscopically-confirmed redshifts of 2.399 (HSC16adga), 1.965 (HSC17auzg), and 1.851 (HSC17dbpf), and two supernova candidates with the host-galaxy photometric redshifts of 3.2 (HSC16apuo) and 4.2 (HSC17dsid), respectively. In this paper, we present their photometric properties and the spectroscopic properties of the confirmed high-redshift supernovae are presented in the accompanying paper (Curtin et al. 2019). The supernovae with the confirmed redshifts of z ~ 2 have the rest ultraviolet peak magnitudes close to -21 mag and they are likely superluminous supernovae. The discovery of three supernovae at z ~ 2 roughly corresponds to the approximate event rate of ~ 900 +/- 520 Gpc-3 yr-1 with Poisson error, which is consistent with the total superluminous supernova rate estimated by extrapolating the local rate based on the cosmic star-formation history. Adding unconfirmed superluminous supernova candidates would increase the event rate. Our superluminous supernova candidates at the redshifts of around 3 and 4 indicate the approximate superluminous supernova rates of ~ 400 +/- 400 Gpc-3 yr-1 (z ~ 3) and ~ 500 +/- 500 Gpc-3 yr-1 (z ~ 4) with Poisson errors. Our initial results demonstrate the outstanding capability of Hyper Suprime-Cam to discover high-redshift supernovae.

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

The origin of fast radio bursts (FRBs), bright millisecond radio transients, is still somewhat of a mystery. Several theoretical models expect that the FRB accompanies an optical afterglow (e.g., Totani 2013; Kashiyama et al. 2013). In order to investigate the origin of FRBs, we perform $gri$-band follow-up observations of FRB~151230 (estimated $z \lesssim 0.8$) with Subaru/Hyper Suprime-Cam at $8$, $11$, and $14$~days after discovery. The follow-up observation reaches a $50\%$ completeness magnitude of $26.5$~mag for point sources, which is the deepest optical follow-up of FRBs to date. We find $13$ counterpart candidates with variabilities during the observation. We investigate their properties with multicolor and multi-wavelength observations and archival catalogs. Two candidates are excluded by the non-detection of FRB~151230 in the other radio feed horns that operated simultaneously to the detection, as well as the inconsistency between the photometric redshift and that derived from the dispersion measure of FRB~151230. Eight further candidates are consistent with optical variability seen in AGNs. Two more candidates are well fitted with transient templates (Type IIn supernovae), and the final candidate is poorly fitted with all of our transient templates and is located off-center of an extended source. It can only be reproduced with rapid transients with a faint peak and rapid decline and the probability of chance coincidence is $\sim3.6\%$. We also find that none of our candidates are consistent with Type Ia supernovae, which rules out the association of Type Ia supernovae to FRB~151230 at $z\leq0.6$ and limits the dispersion measure of the host galaxy to $\lesssim300$~pc~cm$^{-3}$ in a Type Ia supernova scenario.

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

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

Supernovae Type-Ia (SNeIa) play a significant role in exploring the history of the expansion of the Universe, since they are the best-known standard candles with which we can accurately measure the distance to the objects. Finding large samples of SNeIa and investigating their detailed characteristics have become an important issue in cosmology and astronomy. Existing methods relied on a photometric approach that first measures the luminance of supernova candidates precisely and then fits the results to a parametric function of temporal changes in luminance. However, it inevitably requires multi-epoch observations and complex luminance measurements. In this work, we present a novel method for classifying SNeIa simply from single-epoch observation images without any complex measurements, by effectively integrating the state-of-the-art computer vision methodology into the standard photometric approach. Our method first builds a convolutional neural network for estimating the luminance of supernovae from telescope images, and then constructs another neural network for the classification, where the estimated luminance and observation dates are used as features for classification. Both of the neural networks are integrated into a single deep neural network to classify SNeIa directly from observation images. Experimental results show the effectiveness of the proposed method and reveal classification performance comparable to existing photometric methods with multi-epoch observations.

Recent detection of gravitational waves from a neutron star (NS) merger event GW170817 and identification of an electromagnetic counterpart provide a unique opportunity to study the physical processes in NS mergers. To derive properties of ejected material from the NS merger, we perform radiative transfer simulations of kilonova, optical and near-infrared emissions powered by radioactive decays of r-process nuclei synthesized in the merger. We find that the observed near-infrared emission lasting for > 10 days is explained by 0.03 Msun of ejecta containing lanthanide elements. However, the blue optical component observed at the initial phases requires an ejecta component with a relatively high electron fraction (Ye). We show that both optical and near-infrared emissions are simultaneously reproduced by the ejecta with a medium Ye of ~ 0.25. We suggest that a dominant component powering the emission is post-merger ejecta, which exhibits that mass ejection after the first dynamical ejection is quite efficient. Our results indicate that NS mergers synthesize a wide range of r-process elements and strengthen the hypothesis that NS mergers are the origin of r-process elements in the Universe.

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

Type Ia supernovae (SNe Ia) arise from the thermonuclear explosion of carbon-oxygen white dwarfs. Though the uniformity of their light curves makes them powerful cosmological distance indicators, long-standing issues remain regarding their progenitors and explosion mechanisms. Recent detection of the early ultraviolet pulse of a peculiar subluminous SN Ia has been claimed as new evidence for the companion-ejecta interaction through the single-degenerate channel. Here, we report the discovery of a prominent but red optical flash at $\sim$ 0.5 days after the explosion of a SN Ia which shows hybrid features of different SN Ia sub-classes: a light curve typical of normal-brightness SNe Ia, but with strong titanium absorptions, commonly seen in the spectra of subluminous ones. We argue that the early flash of such a hybrid SN Ia is different from predictions of previously suggested scenarios such as the companion-ejecta interaction. Instead it can be naturally explained by a SN explosion triggered by a detonation of a thin helium shell either on a near-Chandrasekhar-mass white dwarf ($\gtrsim$ 1.3 M$_{\odot}$) with low-yield $^{56}$Ni or on a sub-Chandrasekhar-mass white dwarf ($\sim$ 1.0 M$_{\odot}$) merging with a less massive white dwarf. This finding provides compelling evidence that one branch of the previously proposed explosion models, the helium-ignition scenario, does exist in nature, and such a scenario may account for explosions of white dwarfs in a wider mass range in contrast to what was previously supposed.

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

Although Type Ia supernova cosmology has now reached a mature state, it is important to develop as many independent methods as possible to understand the true nature of dark energy. Recent studies have shown that Type II supernovae (SNe II) offer such a path and could be used as alternative distance indicators. However, the majority of these studies were unable to extend the Hubble diagram above redshift $z=0.3$ because of observational limitations. Here, we show that we are now ready to move beyond low redshifts and attempt high-redshift ($z \gtrsim 0.3$) SN~II cosmology as a result of new-generation deep surveys such as the Subaru/Hyper Suprime-Cam (HSC) survey. Applying the "standard candle method" to SN$\sim$2016jhj ($z=0.3398 \pm 0.0002$; discovered by HSC) together with a low-redshift sample, we are able to construct the highest-redshift SN II Hubble diagram to date with an observed dispersion of 0.27 mag (i.e., 12-13% in distance). This work demonstrates the bright future of SN~II cosmology in the coming era of large, wide-field surveys like that of the Large Synoptic Survey Telescope.

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

We present the deepest optical images of the COSMOS field based on a joint dataset taken with Hyper Suprime-Cam (HSC) by the HSC Subaru Strategic Program (SSP) team and the University of Hawaii (UH). The COSMOS field is one of the key extragalactic fields with a wealth of deep, multi-wavelength data. However, the current optical data are not sufficiently deep to match with, e.g., the UltraVista data in the near-infrared. The SSP team and UH have joined forces to produce very deep optical images of the COSMOS field by combining data from both teams. The coadd images reach depths of g=27.8, r=27.7, i=27.6, z=26.8, and y=26.2 mag at 5 sigma for point sources based on flux uncertainties quoted by the pipeline and they cover essentially the entire COSMOS 2 square degree field. The seeing is between 0.6 and 0.9 arcsec on the coadds. We perform several quality checks and confirm that the data are of science quality; ~2% photometry and 30 mas astrometry. This accuracy is identical to the Public Data Release 1 from HSC-SSP. We make the joint dataset including fully calibrated catalogs of detected objects available to the community at https://hsc-release.mtk.nao.ac.jp/.

In this paper, we describe the optical imaging data processing pipeline developed for the Subaru Telescope's Hyper Suprime-Cam (HSC) instrument. The HSC Pipeline builds on the prototype pipeline being developed by the Large Synoptic Survey Telescope's Data Management system, adding customizations for HSC, large-scale processing capabilities, and novel algorithms that have since been reincorporated into the LSST codebase. While designed primarily to reduce HSC Subaru Strategic Program (SSP) data, it is also the recommended pipeline for reducing general-observer HSC data. The HSC pipeline includes high level processing steps that generate coadded images and science-ready catalogs as well as low-level detrending and image characterizations.

We report the serendipitous discovery of a quadruply lensed source at $z_{\rm s}=3.76$, HSC~J115252+004733, from the Hyper Suprime-Cam (HSC) Survey. The source is lensed by an early-type galaxy at $z_{\rm l}=0.466$ and a satellite galaxy. Here, we investigate the properties of the source by studying its size and luminosity from the imaging and the luminosity and velocity width of the Ly-$\alpha$ line from the spectrum. Our analyses suggest that the source is most probably a low-luminosity active galactic nucleus (LLAGN) but the possibility of it being a compact bright galaxy (e.g., a Lyman-$\alpha$ emitter or Lyman Break Galaxy) cannot be excluded. The brighter pair of lensed images appears point-like except in the HSC $i$-band (with a seeing $\sim0.5"$). The extended emission in the $i$-band image could be due to the host galaxy underneath the AGN, or alternatively, due to a highly compact lensed galaxy (without AGN) which appears point-like in all bands except in $i$-band. We also find that the flux ratio of the brighter pair of images is different in the Ks-band compared to optical wavelengths. Phenomena such as differential extinction and intrinsic variability cannot explain this chromatic variation. While microlensing from stars in the foreground galaxy is less likely to be the cause, it cannot be ruled out completely. If the galaxy hosts an AGN, then this represents the highest redshift quadruply imaged AGN known to date, enabling study of a distant LLAGN. Discovery of this unusually compact and faint source demonstrates the potential of the HSC survey.

We compare the isophote shape parameter $a_{4}$ of early-type galaxies (ETGs) between $z\sim1$ and 0 as a proxy for dynamics to investigate the epoch at which the dynamical properties of ETGs are established, using cluster ETG samples with stellar masses of $\log(M_{*}/M_{\odot})\geq10.5$ which have spectroscopic redshifts. We have 130 ETGs from the Hubble Space Telescope Cluster Supernova Survey for $z\sim1$ and 355 ETGs from the Sloan Digital Sky Survey for $z\sim0$. We have developed an isophote shape analysis method which can be used for high-redshift galaxies and has been carefully compared with published results. We have applied the same method for both the $z\sim1$ and $0$ samples. We find similar dependence of the $a_{4}$ parameter on the mass and size at $z\sim1$ and 0; the main population of ETGs changes from disky to boxy at a critical stellar mass of $\log(M_{*}/M_{\odot})\sim11.5$ with the massive end dominated by boxy. The disky ETG fraction decreases with increasing stellar mass both at $z\sim1$ and $0$, and is consistent between these redshifts in all stellar mass bins when the Eddington bias is taken into account. Although uncertainties are large, the results suggest that the isophote shapes and probably dynamical properties of ETGs in massive clusters are already in place at $z>1$ and do not significantly evolve in $z<1$, despite significant size evolution in the same galaxy population. The constant disky fraction favors less violent processes than mergers as a main cause of the size and morphological evolution of intermediate mass ETGs in $z<1$.

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

We have analyzed multi-band light curves of 328 intermediate redshift (0.05 <= z < 0.24) type Ia supernovae (SNe Ia) observed by the Sloan Digital Sky Survey-II Supernova Survey (SDSS-II SN Survey). The multi-band light curves were parameterized by using the Multi-band Stretch Method, which can simply parameterize light curve shapes and peak brightness without dust extinction models. We found that most of the SNe Ia which appeared in red host galaxies (u - r > 2.5) don't have a broad light curve width and the SNe Ia which appeared in blue host galaxies (u - r < 2.0) have a variety of light curve widths. The Kolmogorov-Smirnov test shows that the colour distribution of SNe Ia appeared in red / blue host galaxies is different (significance level of 99.9%). We also investigate the extinction law of host galaxy dust. As a result, we find the value of Rv derived from SNe Ia with medium light curve width is consistent with the standard Galactic value. On the other hand, the value of Rv derived from SNe Ia that appeared in red host galaxies becomes significantly smaller. These results indicate that there may be two types of SNe Ia with different intrinsic colours, and they are obscured by host galaxy dust with two different properties.

We present an application of machine-learning (ML) techniques to source selection in the optical transient survey data with Hyper Suprime-Cam (HSC) on the Subaru telescope. Our goal is to select real transient events accurately and in a timely manner out of a large number of false candidates, obtained with the standard difference-imaging method. We have developed the transient selector which is based on majority voting of three ML machines of AUC Boosting, Random Forest, and Deep Neural Network. We applied it to our observing runs of Subaru-HSC in 2015 May and August, and proved it to be efficient in selecting optical transients. The false positive rate was 1.0% at the true positive rate of 90% in the magnitude range of 22.0--25.0 mag for the former data. For the latter run, we successfully detected and reported ten candidates of supernovae within the same day as the observation. From these runs, we learned the following lessons: (1) the training using artificial objects is effective in filtering out false candidates, especially for faint objects, and (2) combination of ML by majority voting is advantageous.

The Prime Focus Spectrograph (PFS) is a wide-field, multi-object spectrograph accommodating 2394 fibers to observe the sky at the prime focus of the Subaru telescope. The software system to operate a spectroscopic survey is structured by the four packages: Instrument control software, exposure targeting software, data reduction pipeline, and survey planning and tracking software. In addition, we operate a database system where various information such as properties of target objects, instrument configurations, and observation conditions is stored and is organized via a standardized data model for future references to update survey plans and to scientific researches. In this article, we present an overview of the software system and describe the workflows that need to be performed in the PFS operation, with some highlights on the database that organizes various information from sub-processes in the survey operation, and on the process of fiber configuration from the software perspectives.

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

A fading radio source, coincident in time and position with the fast radio burst FRB150418, has been associated with the galaxy WISE J071634.59-190039.2. Subsequent observations of this galaxy have revealed that it contains a persistent, but variable, radio source. We present e-MERLIN, VLBA, and ATCA radio observations and Subaru optical observations of WISE J071634.59-190039.2 and find that the persistent radio source is unresolved and must be compact (<0.01 kpc), and that its location is consistent with the optical centre of the galaxy. We conclude that it is likely that WISE J071634.59-190039.2 contains a weak radio AGN.

We study the extent to which very bright (-23.0 < MUV < -21.75) Lyman-break selected galaxies at redshifts z~7 display detectable Lya emission. To explore this issue, we have obtained follow-up optical spectroscopy of 9 z~7 galaxies from a parent sample of 24 z~7 galaxy candidates selected from the 1.65 sq.deg COSMOS-UltraVISTA and SXDS-UDS survey fields using the latest near-infrared public survey data, and new ultra-deep Subaru z'-band imaging (which we also present and describe in this paper). Our spectroscopy has yielded only one possible detection of Lya at z=7.168 with a rest-frame equivalent width EW_0 = 3.7 (+1.7/-1.1) Angstrom. The relative weakness of this line, combined with our failure to detect Lya emission from the other spectroscopic targets allows us to place a new upper limit on the prevalence of strong Lya emission at these redshifts. For conservative calculation and to facilitate comparison with previous studies at lower redshifts, we derive a 1-sigma upper limit on the fraction of UV bright galaxies at z~7 that display EW_0 > 50 Angstrom, which we estimate to be < 0.23. This result may indicate a weak trend where the fraction of strong Lya emitters ceases to rise, and possibly falls between z~6 and z~7. Our results also leave open the possibility that strong Lya may still be more prevalent in the brightest galaxies in the reionization era than their fainter counterparts. A larger spectroscopic sample of galaxies is required to derive a more reliable constraint on the neutral hydrogen fraction at z~7 based on the Lya fraction in the bright galaxies.

We present a new method to effectively select objects which may be low-mass active black holes (BHs) at galaxy centers using high-cadence optical imaging data, and our first spectroscopic identification of an active 2.7x10^6 Msun BH at z=0.164. This active BH was originally selected due to its rapid optical variability, from a few hours to a day, based on Subaru Hyper Suprime-Cam~(HSC) g-band imaging data taken with 1-hour cadence. Broad and narrow H-alpha and many other emission lines are detected in our optical spectra taken with Subaru FOCAS, and the BH mass is measured via the broad H-alpha emission line width (1,880 km s^{-1}) and luminosity (4.2x10^{40} erg s^{-1}) after careful correction for the atmospheric absorption around 7,580-7,720A. We measure the Eddington ratio to be as low as 0.05, considerably smaller than those in a previous SDSS sample with similar BH mass and redshift, which indicates one of the strong potentials of our Subaru survey. The g-r color and morphology of the extended component indicate that the host galaxy is a star-forming galaxy. We also show effectiveness of our variability selection for low-mass active BHs.

In recent years, millisecond duration radio signals originating from distant galaxies appear to have been discovered in the so-called Fast Radio Bursts. These signals are dispersed according to a precise physical law and this dispersion is a key observable quantity which, in tandem with a redshift measurement, can be used for fundamental physical investigations. While every fast radio burst has a dispersion measurement, none before now have had a redshift measurement, due to the difficulty in pinpointing their celestial coordinates. Here we present the discovery of a fast radio burst and the identification of a fading radio transient lasting $\sim 6$ days after the event, which we use to identify the host galaxy; we measure the galaxy's redshift to be $z=0.492\pm0.008$. The dispersion measure and redshift, in combination, provide a direct measurement of the cosmic density of ionised baryons in the intergalactic medium of $\Omega_{\mathrm{IGM}}=4.9 \pm 1.3\%$, in agreement with the expectation from WMAP, and including all of the so-called "missing baryons". The $\sim6$-day transient is largely consistent with a short gamma-ray burst radio afterglow, and its existence and timescale do not support progenitor models such as giant pulses from pulsars, and supernovae. This contrasts with the interpretation of another recently discovered fast radio burst, suggesting there are at least two classes of bursts.

We present rapidly rising transients discovered by a high-cadence transient survey with Subaru telescope and Hyper Suprime-Cam. We discovered five transients at z=0.384-0.821 showing the rising rate faster than 1 mag per 1 day in the restframe near-ultraviolet wavelengths. The fast rising rate and brightness are the most similar to SN 2010aq and PS1-13arp, for which the ultraviolet emission within a few days after the shock breakout was detected. The lower limit of the event rate of rapidly rising transients is ~9 % of core-collapse supernova rates, assuming a duration of rapid rise to be 1 day. We show that the light curves of the three faint objects agree with the cooling envelope emission from the explosion of red supergiants. The other two luminous objects are, however, brighter and faster than the cooling envelope emission. We interpret these two objects to be the shock breakout from dense wind with the mass loss rate of ~10^{-3} Msun yr^{-1}, as also proposed for PS1-13arp. This mass loss rate is higher than that typically observed for red supergiants. The event rate of these luminous objects is >~1 % of core-collapse supernova rate, and thus, our study implies that more than ~1 % of massive stars can experience an intensive mass loss at a few years before the explosion.

We present measurements of the rates of high-redshift Type Ia supernovae derived from the Subaru/XMM-Newton Deep Survey (SXDS). We carried out repeat deep imaging observations with Suprime-Cam on the Subaru Telescope, and detected 1040 variable objects over 0.918 deg$^2$ in the Subaru/XMM-Newton Deep Field. From the imaging observations, light curves in the observed $i'$-band are constructed for all objects, and we fit the observed light curves with template light curves. Out of the 1040 variable objects detected by the SXDS, 39 objects over the redshift range $0.2 < z < 1.4$ are classified as Type Ia supernovae using the light curves. These are among the most distant SN Ia rate measurements to date. We find that the Type Ia supernova rate increase up to $z \sim 0.8$ and may then flatten at higher redshift. The rates can be fitted by a simple power law, $r_V(z)=r_0(1+z)^\alpha$ with $r_0=0.20^{+0.52}_{-0.16}$(stat.)$^{+0.26}_{-0.07}$(syst.)$\times 10^{-4} {\rm yr}^{-1}{\rm Mpc}^{-3}$, and $\alpha=2.04^{+1.84}_{-1.96}$(stat.)$^{+2.11}_{-0.86}$(syst.).

This paper describes the data release of the Sloan Digital Sky Survey-II (SDSS-II) Supernova Survey conducted between 2005 and 2007. Light curves, spectra, classifications, and ancillary data are presented for 10,258 variable and transient sources discovered through repeat ugriz imaging of SDSS Stripe 82, a 300 deg2 area along the celestial equator. This data release is comprised of all transient sources brighter than r~22.5 mag with no history of variability prior to 2004. Dedicated spectroscopic observations were performed on a subset of 889 transients, as well as spectra for thousands of transient host galaxies using the SDSS-III BOSS spectrographs. Photometric classifications are provided for the candidates with good multi-color light curves that were not observed spectroscopically. From these observations, 4607 transients are either spectroscopically confirmed, or likely to be, supernovae, making this the largest sample of supernova candidates ever compiled. We present a new method for SN host-galaxy identification and derive host-galaxy properties including stellar masses, star-formation rates, and the average stellar population ages from our SDSS multi-band photometry. We derive SALT2 distance moduli for a total of 1443 SN Ia with spectroscopic redshifts as well as photometric redshifts for a further 677 purely-photometric SN Ia candidates. Using the spectroscopically confirmed subset of the three-year SDSS-II SN Ia sample and assuming a flat Lambda-CDM cosmology, we determine Omega_M = 0.315 +/- 0.093 (statistical error only) and detect a non-zero cosmological constant at 5.7 sigmas.

Obtaining reliable, matched photometry for galaxies imaged by different observatories represents a key challenge in the era of wide-field surveys spanning more than several hundred square degrees. Methods such as flux fitting, profile fitting, and PSF homogenization followed by matched-aperture photometry are all computationally expensive. We present an alternative solution called "synthetic aperture photometry" that exploits galaxy profile fits in one band to efficiently model the observed, PSF-convolved light profile in other bands and predict the flux in arbitrarily sized apertures. Because aperture magnitudes are the most widely tabulated flux measurements in survey catalogs, producing synthetic aperture magnitudes (SYNMAGs) enables very fast matched photometry at the catalog level, without reprocessing imaging data. We make our code public and apply it to obtain matched photometry between SDSS ugriz and UKIDSS YJHK imaging, recovering red-sequence colors and photometric redshifts with a scatter and accuracy as good as if not better than FWHM-homogenized photometry from the GAMA Survey. Finally, we list some specific measurements that upcoming surveys could make available to facilitate and ease the use of SYNMAGs.

We present observations of SDF-05M05, an unusual optical transient discovered in the Subaru Deep Field (SDF). The duration of the transient is > ~800 d in the observer frame, and the maximum brightness during observation reached approximately 23 mag in the i' and z' bands. The faint host galaxy is clearly identified in all 5 optical bands of the deep SDF images. The photometric redshift of the host yields z~0.6 and the corresponding absolute magnitude at maximum is ~-20. This implies that this event shone with an absolute magnitude brighter than -19 mag for approximately 300 d in the rest frame, which is significantly longer than a typical supernova and ultra-luminous supernova. The total radiated energy during our observation was 1x10^51 erg. The light curves and color evolution are marginally consistent with some of luminous IIn supernova. We suggest that the transient may be a unique and peculiar supernova at intermediate redshift.

We present a Subaru weak lensing measurement of ACT-CL J0022.2-0036, one of the most luminous, high-redshift (z=0.81) Sunyaev-Zel'dovich (SZ) clusters discovered in the 268 deg^2 equatorial region survey of the Atacama Cosmology Telescope. For the weak lensing analysis using i'-band images, we use a model-fitting (Gauss-Laguerre shapelet) method to measure shapes of galaxy images, where we fit galaxy images in different exposures simultaneously to obtain best-fit ellipticities taking into account the different PSFs in each exposure. We also take into account the astrometric distortion effect on galaxy images by performing the model fitting in the world coordinate system. To select background galaxies behind the cluster at z=0.81, we use photometric redshift (photo-z) estimates for every galaxy derived from the co-added images of multi-passband Br'i'z'Y, with PSF matching/homogenization. After a photo-z cut for background galaxy selection, we detect the tangential weak lensing distortion signal with a total signal-to-noise ratio of about 3.7. By fitting a Navarro-Frenk-White model to the measured shear profile, we find the cluster mass to be M_200\bar{\rho}_m = [7.5^+3.2_-2.8(stat.)^+1.3_-0.6(sys.)] x 10^14 M_\odot/h. The weak lensing-derived mass is consistent with previous mass estimates based on the SZ observation, with assumptions of hydrostatic equilibrium and virial theorem, as well as with scaling relations between SZ signal and mass derived from weak lensing, X-ray, and velocity dispersion, within the measurement errors. We also show that the existence of ACT-CL J0022.2-0036 at z=0.81 is consistent with the cluster abundance prediction of the \Lambda-dominated cold dark matter structure formation model. We thus demonstrate the capability of Subaru-type ground-based images for studying weak lensing of high-redshift clusters.

Using the sample of Type Ia supernovae (SNe Ia) discovered by the Hubble Space Telescope (HST) Cluster Supernova Survey and augmented with HST-observed SNe Ia in the GOODS fields, we search for correlations between the properties of SNe and their host galaxies at high redshift. We use galaxy color and quantitative morphology to determine the red sequence in 25 clusters and develop a model to distinguish passively evolving early-type galaxies from star-forming galaxies in both clusters and the field. With this approach, we identify six SN Ia hosts that are early-type cluster members and eleven SN Ia hosts that are early-type field galaxies. We confirm for the first time at z>0.9 that SNe Ia hosted by early-type galaxies brighten and fade more quickly than SNe Ia hosted by late-type galaxies. We also show that the two samples of hosts produce SNe Ia with similar color distributions. The relatively simple spectral energy distributions (SEDs) expected for passive galaxies enable us to measure stellar masses of early-type SN hosts. In combination with stellar mass estimates of late-type GOODS SN hosts from Thomson & Chary (2011), we investigate the correlation of host mass with Hubble residual observed at lower redshifts. Although the sample is small and the uncertainties are large, a hint of this relation is found at z>0.9. By simultaneously fitting the average cluster galaxy formation history and dust content to the red-sequence scatters, we show that the reddening of early-type cluster SN hosts is likely E(B-V) <~ 0.06. The similarity of the field and cluster early-type host samples suggests that field early-type galaxies that lie on the red sequence may also be minimally affected by dust. Hence, the early-type hosted SNe Ia studied here occupy a more favorable environment to use as well-characterized high-redshift standard candles than other SNe Ia.

We report a measurement of the Type Ia supernova (SN Ia) rate in galaxy clusters at 0.9 < z < 1.45 from the Hubble Space Telescope (HST) Cluster Supernova Survey. This is the first cluster SN Ia rate measurement with detected z > 0.9 SNe. Finding 8 +/- 1 cluster SNe Ia, we determine a SN Ia rate of 0.50 +0.23-0.19 (stat) +0.10-0.09 (sys) SNuB (SNuB = 10^-12 SNe L_{sun,B}^-1 yr^-1). In units of stellar mass, this translates to 0.36 +0.16-0.13 (stat) +0.07-0.06 (sys) SNuM (SNuM = 10^-12 SNe M_sun^-1 yr^-1). This represents a factor of approximately 5 +/- 2 increase over measurements of the cluster rate at z < 0.2. We parameterize the late-time SN Ia delay time distribution with a power law (proportional to t^s). Under the assumption of a cluster formation redshift of z_f = 3, our rate measurement in combination with lower-redshift cluster SN Ia rates constrains s = -1.41 +0.47/-0.40, consistent with measurements of the delay time distribution in the field. This measurement is generally consistent with expectations for the "double degenerate" scenario and inconsistent with some models for the "single degenerate" scenario predicting a steeper delay time distribution at large delay times. We check for environmental dependence and the influence of younger stellar populations by calculating the rate specifically in cluster red-sequence galaxies and in morphologically early-type galaxies, finding results similar to the full cluster rate. Finally, the upper limit of one host-less cluster SN Ia detected in the survey implies that the fraction of stars in the intra-cluster medium is less than 0.47 (95% confidence), consistent with measurements at lower redshifts.

We present a measurement of the volumetric Type Ia supernova (SN Ia) rate out to z ~ 1.6 from the Hubble Space Telescope Cluster Supernova Survey. In observations spanning 189 orbits with the Advanced Camera for Surveys we discovered 29 SNe, of which approximately 20 are SNe Ia. Twelve of these SNe Ia are located in the foregrounds and backgrounds of the clusters targeted in the survey. Using these new data, we derive the volumetric SN Ia rate in four broad redshift bins, finding results consistent with previous measurements at z > 1 and strengthening the case for a SN Ia rate that is equal to or greater than ~0.6 x 10^-4/yr/Mpc^3 at z ~ 1 and flattening out at higher redshift. We provide SN candidates and efficiency calculations in a form that makes it easy to rebin and combine these results with other measurements for increased statistics. Finally, we compare the assumptions about host-galaxy dust extinction used in different high-redshift rate measurements, finding that different assumptions may induce significant systematic differences between measurements.

Building on the legacy of the Sloan Digital Sky Survey (SDSS-I and II), SDSS-III is a program of four spectroscopic surveys on three scientific themes: dark energy and cosmological parameters, the history and structure of the Milky Way, and the population of giant planets around other stars. In keeping with SDSS tradition, SDSS-III will provide regular public releases of all its data, beginning with SDSS DR8 (which occurred in Jan 2011). This paper presents an overview of the four SDSS-III surveys. BOSS will measure redshifts of 1.5 million massive galaxies and Lya forest spectra of 150,000 quasars, using the BAO feature of large scale structure to obtain percent-level determinations of the distance scale and Hubble expansion rate at z<0.7 and at z~2.5. SEGUE-2, which is now completed, measured medium-resolution (R=1800) optical spectra of 118,000 stars in a variety of target categories, probing chemical evolution, stellar kinematics and substructure, and the mass profile of the dark matter halo from the solar neighborhood to distances of 100 kpc. APOGEE will obtain high-resolution (R~30,000), high signal-to-noise (S/N>100 per resolution element), H-band (1.51-1.70 micron) spectra of 10^5 evolved, late-type stars, measuring separate abundances for ~15 elements per star and creating the first high-precision spectroscopic survey of all Galactic stellar populations (bulge, bar, disks, halo) with a uniform set of stellar tracers and spectral diagnostics. MARVELS will monitor radial velocities of more than 8000 FGK stars with the sensitivity and cadence (10-40 m/s, ~24 visits per star) needed to detect giant planets with periods up to two years, providing an unprecedented data set for understanding the formation and dynamical evolution of giant planet systems. (Abridged)

The type Ia supernova (SN Ia) rate, when compared to the cosmic star formation history (SFH), can be used to derive the delay-time distribution (DTD) of SNe Ia, which can distinguish among progenitor models. We present the results of a SN survey in the Subaru Deep Field (SDF). Over a period of 3 years, we have observed the SDF on 4 independent epochs with Suprime-Cam on the Subaru 8.2-m telescope, with 2 nights of exposure per epoch, in the R, i', and z' bands. We have discovered 150 SNe out to redshift z~2. Using 11 photometric bands from the observer-frame far-ultraviolet to the near-infrared, we derive photometric redshifts for the SN host galaxies (for 24 we also have spectroscopic redshifts). This information is combined with the SN photometry to determine the type and redshift distribution of the SN sample. Our final sample includes 28 SNe Ia in the range 1.0<z<1.5 and 10 in the range 1.5<z<2.0. Our SN Ia rate measurements are consistent with those derived from the HST/GOODS sample, but the overall uncertainty of our 1.5<z<2.0 measurement is a factor of 2 smaller, of 35-50%. We find that the SN Ia rate evolution levels off at 1.0<z<2.0, but shows no sign of declining. Combining our SN Ia rate measurements and those from the literature, and comparing to a wide range of possible SFHs, the best-fit DTD is a power law of the form Psi(t) ~ t^beta, with beta = -1.1 +/- 0.1 (statistical) +/- 0.17 (systematic). This result is consistent with other recent DTD measurements at various redshifts and environments, and is in agreement with a generic prediction of the double-degenerate progenitor scenario for SNe Ia. By combining the contribution to iron production from core-collapse SNe, based on the wide range of SFHs, with that from SNe Ia, calculated with the best-fit DTD, we predict that the mean present-day cosmic iron abundance is in the range Z(Fe) = (0.09-0.37) Z(Fe,Sun) (abridged).

We present ACS, NICMOS, and Keck AO-assisted photometry of 20 Type Ia supernovae SNe Ia from the HST Cluster Supernova Survey. The SNe Ia were discovered over the redshift interval 0.623 < z < 1.415. Fourteen of these SNe Ia pass our strict selection cuts and are used in combination with the world's sample of SNe Ia to derive the best current constraints on dark energy. Ten of our new SNe Ia are beyond redshift $z=1$, thereby nearly doubling the statistical weight of HST-discovered SNe Ia beyond this redshift. Our detailed analysis corrects for the recently identified correlation between SN Ia luminosity and host galaxy mass and corrects the NICMOS zeropoint at the count rates appropriate for very distant SNe Ia. Adding these supernovae improves the best combined constraint on the dark energy density \rho_{DE}(z) at redshifts 1.0 < z < 1.6 by 18% (including systematic errors). For a LambdaCDM universe, we find \Omega_\Lambda = 0.724 +0.015/-0.016 (68% CL including systematic errors). For a flat wCDM model, we measure a constant dark energy equation-of-state parameter w = -0.985 +0.071/-0.077 (68% CL). Curvature is constrained to ~0.7% in the owCDM model and to ~2% in a model in which dark energy is allowed to vary with parameters w_0 and w_a. Tightening further the constraints on the time evolution of dark energy will require several improvements, including high-quality multi-passband photometry of a sample of several dozen z>1 SNe Ia. We describe how such a sample could be efficiently obtained by targeting cluster fields with WFC3 on HST.

We present the temporal evolution of line profiles ranging from near ultraviolet to optical wavelengths by analyzing 59 Subaru telescope spectra of normal Type Ia Supernovae (SNe Ia) in the intermediate redshift range (0.05 < z < 0.4) discovered by the Sloan Digital Sky Survey-II (SDSS-II) Supernova Survey. We derive line velocities, peak wavelengths and pseudo-equivalent widths (pEWs) of these lines. Additionally, we compare the line profiles around the date of maximum brightness with those from their nearby counterparts. We find that line profiles represented by their velocities and pEWs for intermediate redshift SNe Ia are consistent with their nearby counterparts within 2 $\sigma$. These findings support the picture that SNe Ia are a "standard" candle for the intermediate redshift range as has been shown between SNe Ia at nearby and high redshifts. There is a hint that the "MgII \lambda 4300" pEW distribution for intermediate redshift SNe Ia is larger than for the nearby sample, which could be interpreted as a difference in the progenitor abundance.

The Sloan Digital Sky Survey (SDSS) started a new phase in August 2008, with new instrumentation and new surveys focused on Galactic structure and chemical evolution, measurements of the baryon oscillation feature in the clustering of galaxies and the quasar Ly alpha forest, and a radial velocity search for planets around ~8000 stars. This paper describes the first data release of SDSS-III (and the eighth counting from the beginning of the SDSS). The release includes five-band imaging of roughly 5200 deg^2 in the Southern Galactic Cap, bringing the total footprint of the SDSS imaging to 14,555 deg^2, or over a third of the Celestial Sphere. All the imaging data have been reprocessed with an improved sky-subtraction algorithm and a final, self-consistent photometric recalibration and flat-field determination. This release also includes all data from the second phase of the Sloan Extension for Galactic Understanding and Evolution (SEGUE-2), consisting of spectroscopy of approximately 118,000 stars at both high and low Galactic latitudes. All the more than half a million stellar spectra obtained with the SDSS spectrograph have been reprocessed through an improved stellar parameters pipeline, which has better determination of metallicity for high metallicity stars.

We present the dependences of the properties of type Ia Supernovae (SNe Ia) on their host galaxies by analyzing the multi-band lightcurves of 118 spectroscopically confirmed SNe Ia observed by the Sloan Digital Sky Survey (SDSS) Supernova Survey and the spectra of their host galaxies. We derive the equivalent width of the \rm{H}$\alpha$ emission line, star formation rate, and gas-phase metallicity from the spectra and compare these with the lightcurve widths and colors of SNe Ia. In addition, we compare host properties with the deviation of the observed distance modulus corrected for lightcurve parameters from the distance modulus determined by the best fit cosmological parameters. This allows us to investigate uncorrected systematic effects in the magnitude standardization. We find that SNe Ia in host galaxies with a higher star formation rate have synthesized on average a larger $^{56}$Ni mass and show wider lightcurves. The $^{56}$Ni mass dependence on metallicity is consistent with a prediction of Timmes et al. 2003 based on nucleosynthesis. SNe Ia in metal-rich galaxies ({$\log_{10}(O/H)+12>8.9$) have become 0.13 $\pm$ 0.06 magnitude brighter after corrections for their lightcurve widths and colors, which corresponds to up to 6% uncertainty in the luminosity distance. We investigate whether parameters for standardizing SN Ia maximum magnitude differ among samples with different host characteristics. The coefficient of the color term is larger by 0.67 $\pm$ 0.19 for SNe Ia in metal-poor hosts than those in metal-rich hosts when no color cuts are imposed.

The Sloan Digital Sky Survey II (SDSS-II) Supernova Survey discovered Type Ia supernovae (SNe Ia) in an almost unexplored intermediate redshift range of $0.05 < z < 0.4$ and provided densely sampled multi-color light curves for SN candidates. Followup spectroscopy of this survey was carried out with the Subaru telescope and spectra of 71 SN Ia candidates were obtained. One spectrum was observed per candidate except for a peculiar variable. This paper presents the method for processing these spectra. The observed wavelength ranges of our spectra are 4000 to 9000 {\AA} for Year 2005 and 3600 to 9000 {\AA} for Year 2006. Most SN Ia spectra have signal to noise ratios (S/N) between 4 and 10 per 2 {\AA} averaged over the entire wavelength region. We developed a new code to extract the SN spectral component from spectra contaminated by the host galaxy. Of 71 SN Ia candidates, 59 are identified as normal SNe Ia and 3 are peculiar SNe Ia. The range of spectral phases varies from -7 days to +30 days from maximum brightness. There are also 7 SNe II, 1 possible hypernova and 1 AGN.

We present the result of projected cross correlation analysis of AGNs and galaxies at redshifts from 0.3 to 3.0. The Japanese Virtual Observatory (JVO) was used to obtain the Subaru Suprime-Cam images and UKIDSS catalog data around AGNs. We investigated 1,809 AGNs, which is about ten times larger a sample than that used in previous studies on AGN-galaxy clustering at redshifts larger than 0.6. 90% of the AGN samples are optically-selected AGN from the SDSS and 2dF catalogs. The galaxy samples at low redshift includes many redder objects from UKIDSS survey, while at higher redshift they are mainly blue galaxies from Suprime-Cam. We found significant excess of galaxies around the AGNs at redshifts from 0.3 to 1.8. For the low redshift samples ($z<0.9$), we obtained correlation length of $r_{0} = $5--6 $h^{-1}$Mpc ($\gamma = 1.8$), which indicates that the AGNs at this redshift range reside in a similar environment around typical local galaxies. We also found that AGNs at higher redshift ranges reside in a denser environment than lower redshift AGNs; For $z=1.3 \sim 1.8$ AGNs, the cross correlation length was measured as 11$^{+6}_{-3}$ $h^{-1}$Mpc ($\gamma=1.8$). Considering that our galaxies sample is based on optical observations with Suprime-Cam at the redshift range, it is expected that blue star-forming galaxies comprise the majority of objects that are observed to be clustered around the AGNs. It is successfully demonstrated that the use of the archive through the Virtual Observatory system can provide a powerful tool for investigating the small scale environment of the intermediate redshift AGNs.

We study the photometric properties of stars in the data archive of the Sloan Digital Sky Survey (SDSS), the prime aim being to understand the photometric calibration over the entire data set. It is confirmed that the photometric calibration for point sources has been made overall tightly against the SDSS standard stars. We have also confirmed that photometric synthesis of the SDSS spectrophotometric data gives broad band fluxes that agree with broad band photometry with errors no more than 0.04 mag and little tilt along the wide range of colours, verifying that the response functions of the SDSS 2.5 m telescope system are well characterised. We locate stars in the SDSS photometric system, so that stars can roughly be classified into spectral classes from the colour information. We show how metallicity and surface gravity affect colours, and that stars contained in the SDSS general catalogue, plotted in colour space, show the distribution that matches well with what is anticipated from the variations of metallicity and surface gravity. The colour-colour plots are perfectly consistent among the three samples, stars in the SDSS general catalogue, SDSS standard stars and spectrophotometric stars of Gunn & Stryker, especially when some considerations are taken into account of the differences (primarily metallicity) of the samples. We show that the g-r - inverse temperature relation is tight and can be used as a good estimator of the effective temperature of stars over a fairly wide range of effective temperatures. We also confirm that the colours of G2V stars in the SDSS photometric system match well with the Sun.

We investigate the stellar kinematics and line strength indices in the outer halos of brightest cluster galaxies (BCGs) in the Coma cluster to obtain the outer halo V_rot and sigma profiles and to derive constraints on the formation history of these objects. Methods: We analyzed absorption lines in deep, medium-resolution, long-slit spectra in the wavelength range ~ 4500 - 5900 Angstrom, out to ~50 kpc for NGC 4874 and ~65 kpc for NGC 4889, probing regions with a surface brightness down to mu_R ~24 mag/arcsec^2. Results: These data provide stellar velocity and velocity dispersion profiles along the major axes of both BCGs, and also along the minor axis of NGC 4889. The kinematic properties of NGC 4874 and NGC 4889 halos extend the previous relations of early-type galaxy halos to bright luminosities and indicate that the stars in the outer regions are still bound to these galaxies. For NGC 4889 we also determine Hbeta, Mg and Fe line strength indices, finding strong radial gradients for Mg and Fe. The current dataset for NGC 4889 is one of the most extended in radius, including both stellar kinematics AND line strength index measurements.

We report on work to increase the number of well-measured Type Ia supernovae (SNe Ia) at high redshifts. Light curves, including high signal-to-noise HST data, and spectra of six SNe Ia that were discovered during 2001 are presented. Additionally, for the two SNe with z>1, we present ground-based J-band photometry from Gemini and the VLT. These are among the most distant SNe Ia for which ground based near-IR observations have been obtained. We add these six SNe Ia together with other data sets that have recently become available in the literature to the Union compilation (Kowalski et al. 2008). We have made a number of refinements to the Union analysis chain, the most important ones being the refitting of all light curves with the SALT2 fitter and an improved handling of systematic errors. We call this new compilation, consisting of 557 supernovae, the Union2 compilation. The flat concordance LambdaCDM model remains an excellent fit to the Union2 data with the best fit constant equation of state parameter w=-0.997^{+0.050}_{-0.054} (stat) ^{+0.077}_{-0.082} (stat+sys\ together) for a flat universe, or w=-1.035^{+0.055}_{-0.059} (stat)^{+0.093}_{-0.097} (stat+sys together) with curvature. We also present improved constraints on w(z). While no significant change in w with redshift is detected, there is still considerable room for evolution in w. The strength of the constraints depend strongly on redshift. In particular, at z > 1, the existence and nature of dark energy are only weakly constrained by the data.

The monochromatic illumination system is constructed to carry out in situ measurements of the response function of the mosaicked CCD imager used in the Sloan Digital Sky Survey (SDSS). The system is outlined and the results of the measurements, mostly during the first 6 years of the SDSS, are described. We present the reference response functions for the five colour passbands derived from these measurements, and discuss column to column variations and variations in time, and also their effects on photometry. We also discuss the effect arising from various, slightly different response functions of the associated detector systems that were used to give SDSS photometry. We show that the calibration procedures of SDSS remove these variations reasonably well with the resulting final errors from variant response functions being unlikely to be larger than 0.01 mag for g, r, i, and z bands over the entire duration of the survey. The considerable aging effect is uncovered in the u band, the response function showing a 30% decrease in the throughput in the short wavelength side during the survey years, which potentially causes a systematic error in photometry. The aging effect is consistent with variation of the instrumental sensitivity in u-band, which is calibrated out. The expected colour variation is consistent with measured colour variation in the catalog of repeated photometry. The colour variation is delta (u-g) ~ 0.01 for most stars, and at most delta (u-g) ~ 0.02 mag for those with extreme colours. We verified in the final catalogue that no systematic variations in excess of 0.01 mag are detected in the photometry which can be ascribed to aging and/or seasonal effects except for the secular u-g colour variation for stars with extreme colours.

We search for stars with proper motions in a set of deep Subaru images, covering about 0.48 square degrees to a depth of $i' \simeq 26$, taken over a span of five and a half years. We follow the methods described in \citet{Richmond2009} to reduce and analyze this dataset. We present a sample of 69 stars with motions of high significance, and discuss briefly the populations from which they are likely drawn. Based on photometry and motions alone, we expect that 14 of the candidates may be white dwarfs. Our candidate with the largest proper motion is surprisingly faint and likely to prove interesting: its colors and motions suggest that it might be an M dwarf moving at over 500 km/sec or an L dwarf in the halo.

We apply the Standardized Candle Method (SCM) for Type II Plateau supernovae (SNe II-P), which relates the velocity of the ejecta of a SN to its luminosity during the plateau, to 15 SNe II-P discovered over the three season run of the Sloan Digital Sky Survey - II Supernova Survey. The redshifts of these SNe - 0.027 < z < 0.144 - cover a range hitherto sparsely sampled in the literature; in particular, our SNe II-P sample contains nearly as many SNe in the Hubble flow (z > 0.01) as all of the current literature on the SCM combined. We find that the SDSS SNe have a very small intrinsic I-band dispersion (0.22 mag), which can be attributed to selection effects. When the SCM is applied to the combined SDSS-plus-literature set of SNe II-P, the dispersion increases to 0.29 mag, larger than the scatter for either set of SNe separately. We show that the standardization cannot be further improved by eliminating SNe with positive plateau decline rates, as proposed in Poznanski et al. (2009). We thoroughly examine all potential systematic effects and conclude that for the SCM to be useful for cosmology, the methods currently used to determine the Fe II velocity at day 50 must be improved, and spectral templates able to encompass the intrinsic variations of Type II-P SNe will be needed.

We present spectra of high-redshift supernovae (SNe) that were taken with the Subaru low resolution optical spectrograph, FOCAS. These SNe were found in SN surveys with Suprime-Cam on Subaru, the CFH12k camera on the Canada-France-Hawaii Telescope (CFHT), and the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope (HST). These SN surveys specifically targeted z>1 Type Ia supernovae (SNe Ia). From the spectra of 39 candidates, we obtain redshifts for 32 candidates and spectroscopically identify 7 active candidates as probable SNe Ia, including one at z=1.35, which is the most distant SN Ia to be spectroscopically confirmed with a ground-based telescope. An additional 4 candidates are identified as likely SNe Ia from the spectrophotometric properties of their host galaxies. Seven candidates are not SNe Ia, either being SNe of another type or active galactic nuclei. When SNe Ia are observed within a week of maximum light, we find that we can spectroscopically identify most of them up to z=1.1. Beyond this redshift, very few candidates were spectroscopically identified as SNe Ia. The current generation of super red-sensitive, fringe-free CCDs will push this redshift limit higher.

We report the discovery of a multiply lensed Ly Alpha (Lya) emitter at z = 3.90 behind the massive galaxy cluster WARPS J1415.1+3612 at z = 1.026. Images taken by the Hubble Space Telescope(HST) using ACS reveal a complex lensing system that produces a prominent, highly magnified arc and a triplet of smaller arcs grouped tightly around a spectroscopically confirmed cluster member. Spectroscopic observations using FOCAS on Subaru confirm strong Lya emission in the source galaxy and provide redshifts for more than 21 cluster members, from which we obtain a velocity dispersion of 807+/-185 km/s. Assuming a singular isothermal sphere profile, the mass within the Einstein ring (7.13+/-0.38") corresponds to a central velocity dispersion of 686+15-19 km/s for the cluster, consistent with the value estimated from cluster member redshifts. Our mass profile estimate from combining strong lensing and dynamical analyses is in good agreement with both X-ray and weak lensing results.

We present ugriz light curves for 146 spectroscopically confirmed or spectroscopically probable Type Ia supernovae from the 2005 season of the SDSS-II Supernova survey. The light curves have been constructed using a photometric technique that we call scene modelling, which is described in detail here; the major feature is that supernova brightnesses are extracted from a stack of images without spatial resampling or convolution of the image data. This procedure produces accurate photometry along with accurate estimates of the statistical uncertainty, and can be used to derive photometry taken with multiple telescopes. We discuss various tests of this technique that demonstrate its capabilities. We also describe the methodology used for the calibration of the photometry, and present calibrated magnitudes and fluxes for all of the spectroscopic SNe Ia from the 2005 season.

We present measurements of the Hubble diagram for 103 Type Ia supernovae (SNe) with redshifts 0.04 < z < 0.42, discovered during the first season (Fall 2005) of the Sloan Digital Sky Survey-II (SDSS-II) Supernova Survey. These data fill in the redshift "desert" between low- and high-redshift SN Ia surveys. We combine the SDSS-II measurements with new distance estimates for published SN data from the ESSENCE survey, the Supernova Legacy Survey, the Hubble Space Telescope, and a compilation of nearby SN Ia measurements. Combining the SN Hubble diagram with measurements of Baryon Acoustic Oscillations from the SDSS Luminous Red Galaxy sample and with CMB temperature anisotropy measurements from WMAP, we estimate the cosmological parameters w and Omega_M, assuming a spatially flat cosmological model (FwCDM) with constant dark energy equation of state parameter, w. For the FwCDM model and the combined sample of 288 SNe Ia, we find w = -0.76 +- 0.07(stat) +- 0.11(syst), Omega_M = 0.306 +- 0.019(stat) +- 0.023(syst) using MLCS2k2 and w = -0.96 +- 0.06(stat) +- 0.12(syst), Omega_M = 0.265 +- 0.016(stat) +- 0.025(syst) using the SALT-II fitter. We trace the discrepancy between these results to a difference in the rest-frame UV model combined with a different luminosity correction from color variations; these differences mostly affect the distance estimates for the SNLS and HST supernovae. We present detailed discussions of systematic errors for both light-curve methods and find that they both show data-model discrepancies in rest-frame $U$-band. For the SALT-II approach, we also see strong evidence for redshift-dependence of the color-luminosity parameter (beta). Restricting the analysis to the 136 SNe Ia in the Nearby+SDSS-II samples, we find much better agreement between the two analysis methods but with larger uncertainties.

A bulge-disk decomposition is made for 737 spiral and lenticular galaxies drawn from a SDSS galaxy sample for which morphological types are estimated. We carry out the bulge-disk decomposition using the growth curve fitting method. It is found that bulge properties, effective radius, effective surface brightness, and also absolute magnitude, change systematically with the morphological sequence; from early to late types, the size becomes somewhat larger, and surface brightness and luminosity fainter. In contrast disks are nearly universal, their properties remaining similar among disk galaxies irrespective of detailed morphologies from S0 to Sc. While these tendencies were often discussed in previous studies, the present study confirms them based on a large homogeneous magnitude-limited field galaxy sample with morphological types estimated. The systematic change of bulge-to-total luminosity ratio, $B/T$, along the morphological sequence is therefore not caused by disks but mostly by bulges. It is also shown that elliptical galaxies and bulges of spiral galaxies are unlikely to be in a single sequence. We infer the stellar mass density (in units of the critical mass density) to be $\Omega=$0.0021 for spheroids, i.e., elliptical galaxies plus bulges of spiral galaxies, and $\Omega=$0.00081 for disks.

The intracluster light (ICL) is a faint diffuse stellar component in clusters made of stars not bound to individual galaxies. We have carried out a large scale study of this component in the nearby Virgo cluster. The diffuse light is traced using planetary nebulae (PNe). The PNe are detected in the on-band image due to their strong emission in the [OIII] 5007 line, but disappear in the off-band image. The contribution of Ly-alpha emitters at z=3.14 are corrected statistically using blank field surveys. We have surveyed a total area of 3.3 square degrees in the Virgo cluster with eleven fields located at different radial distances. Those fields located at smaller radii than 80 arcmin from the cluster center contain most of the detected diffuse light. In this central region of the cluster, the ICL has a surface brightness in the range 28.8 - 30 mag per sqarsec in the B band, it is not uniformly distributed, and represents about 7% of the total galaxy light in this area. At distances larger than 80 arcmin the ICL is confined to single fields and individual sub-structures, e.g. in the Virgo sub-clump B, the M60/M59 group. For several fields at 2 and 3 degrees from the Virgo cluster center we set only upper limits. These results indicate that the ICL is not homogeneously distributed in the Virgo core, and it is concentrated in the high density regions of the Virgo cluster, e.g. the cluster core and other sub-structures. Outside these regions, the ICL is confined within areas of 100 kpc in size, where tidal effects may be at work. These observational results link the formation of the ICL with the formation history of the most luminous cluster galaxies.

We present a new survey strategy to discover and study high redshift Type Ia supernovae (SNe Ia) using the Hubble Space Telescope (HST). By targeting massive galaxy clusters at 0.9<z<1.5, we obtain a twofold improvement in the efficiency of finding SNe compared to an HST field survey and a factor of three improvement in the total yield of SN detections in relatively dust-free red-sequence galaxies. In total, sixteen SNe were discovered at z>0.95, nine of which were in galaxy clusters. This strategy provides a SN sample that can be used to decouple the effects of host galaxy extinction and intrinsic color in high redshift SNe, thereby reducing one of the largest systematic uncertainties in SN cosmology.

The sample of 137 low-redshift type Ia supernovae with 0.05 < z < 0.3 obtained from the SDSS-II Supernova Survey for the southern equatorial stripe of 300 square degrees is used to derive the luminosity functions of type Ia supernovae and of their host galaxies in the gri passbands. We show that the luminosity function of type Ia supernova host galaxies matches well with that of galaxies in the general field, suggesting that the occurrence of type Ia supernovae does not favour a particular type of galaxies but is predominantly proportional to the luminosity of galaxies. The only evidence that points to possible correlation between the supernova rate and star formation activity is that the supernova rate in late-type galaxies is higher than that in early-type galaxies by 31+/-35%. The sample contains 8 type Ia supernovae whose host galaxies were not identified, but it is shown that their occurrence is consistent with them occurred in low luminous galaxies beyond the survey. The luminosity function of type Ia supernovae is approximately Gaussian with the full-width half maximum being a factor of 1.4 in luminosity. The Gaussian distribution becomes tighter if the ratio of extinction to reddening, R_V, is lower than the characteristic value for the Milky Way and if luminosity is corrected for the light curve shape. The colour excess is ~0.07 mag which is significantly smaller than reddening expected for field galaxies. This colour excess does not vary with the distance of the supernovae from the centre of the host galaxy to 15 kpc. This suggests that the major part of the colour excess appears to be either intrinsic or reddening that arises in the immediate environment of supernova, rather than interstellar reddening in host galaxies.

The delay time distribution (DTD) of type Ia supernovae (SNe Ia) from star formation is an important clue to reveal the still unknown progenitor system of SNe Ia. Here we report on a measurement of the SN Ia DTD in a delay time range of t_Ia = 0.1-8.0 Gyr by using the faint variable objects detected in the Subaru/XMM-Newton Deep Survey (SXDS) down to i' ~ 25.5. We select 65 SN candidates showing significant spatial offset from nuclei of the host galaxies having old stellar population at z ~ 0.4-1.2, out of more than 1,000 SXDS variable objects. Although spectroscopic type classification is not available for these, we quantitatively demonstrate that more than ~80% of these should be SNe Ia. The DTD is derived using the stellar age estimates of the old galaxies based on 9 band photometries from optical to mid-infrared wavelength. Combined with the observed SN Ia rate in elliptical galaxies at the local universe, the DTD in t_Ia ~ 0.1-10 Gyr is well described by a featureless power-law as f_D(t_Ia) \propto t_Ia^{-1}. The derived DTD is in excellent agreement with the generic prediction of the double-degenerate scenario, giving a strong support to this scenario. In the single-degenerate (SD) scenario, although predictions by simple analytic formulations have broad DTD shapes that are similar to the observation, DTD shapes calculated by more detailed binary population synthesis tend to have strong peaks at characteristic time scales, which do not fit the observation. This result thus indicates either that the SD channel is not the major contributor to SNe Ia in old stellar population, or that improvement of binary population synthesis theory is required. Various sources of systematic uncertainties are examined and tested, but our main conclusions are not affected significantly.

We search for stars with proper motions in a set of twenty deep Subaru images, covering about 0.28 square degrees to a depth of i ~ 25, taken over a span of six years. In this paper, we describe in detail our reduction and techniques to identify moving objects. We present a first sample of 99 stars with motions of high significance, and discuss briefly the populations from which they are likely drawn. Based on photometry and motions alone, we expect that 9 of the candidates may be white dwarfs. We also find a group of stars which may be extremely metal-poor subdwarfs in the halo.

We present observations of SCP 06F6, an unusual optical transient discovered during the Hubble Space Telescope Cluster Supernova Survey. The transient brightened over a period of ~100 days, reached a peak magnitude of ~21.0 in both i_775 and z_850, and then declined over a similar timescale. There is no host galaxy or progenitor star detected at the location of the transient to a 3 sigma upper limit of i_775 = 26.4 and z_850 = 26.1, giving a corresponding lower limit on the flux increase of a factor of ~120. Multiple spectra show five broad absorption bands between 4100 AA and 6500 AA and a mostly featureless continuum longward of 6500 AA. The shape of the lightcurve is inconsistent with microlensing. The transient's spectrum, in addition to being inconsistent with all known supernova types, is not matched to any spectrum in the Sloan Digital Sky Survey (SDSS) database. We suggest that the transient may be one of a new class.

We present a measurement of the rate of type Ia supernovae (SNe Ia) from the first of three seasons of data from the SDSS-II Supernova Survey. For this measurement, we include 17 SNe Ia at redshift $z\le0.12$. Assuming a flat cosmology with $\Omega_m = 0.3=1-\Omega_\Lambda$, we find a volumetric SN Ia rate of $[2.93^{+0.17}_{-0.04}({\rm systematic})^{+0.90}_{-0.71}({\rm statistical})] \times 10^{-5} {\rm SNe} {\rm Mpc}^{-3} h_{70}^3 {\rm year}^{-1}$, at a volume-weighted mean redshift of 0.09. This result is consistent with previous measurements of the SN Ia rate in a similar redshift range. The systematic errors are well controlled, resulting in the most precise measurement of the SN Ia rate in this redshift range. We use a maximum likelihood method to fit SN rate models to the SDSS-II Supernova Survey data in combination with other rate measurements, thereby constraining models for the redshift-evolution of the SN Ia rate. Fitting the combined data to a simple power-law evolution of the volumetric SN Ia rate, $r_V \propto (1+z)^{\beta}$, we obtain a value of $\beta = 1.5 \pm 0.6$, i.e. the SN Ia rate is determined to be an increasing function of redshift at the $\sim 2.5 \sigma$ level. Fitting the results to a model in which the volumetric SN rate, $r_V=A\rho(t)+B\dot \rho(t)$, where $\rho(t)$ is the stellar mass density and $\dot \rho(t)$ is the star formation rate, we find $A = (2.8 \pm 1.2) \times 10^{-14} \mathrm{SNe} \mathrm{M}_{\sun}^{-1} \mathrm{year}^{-1}$, $B = (9.3^{+3.4}_{-3.1})\times 10^{-4} \mathrm{SNe} \mathrm{M}_{\sun}^{-1}$.

We create new U, B, V, R and I-band light curve templates of type Ia Supernovae (SNe Ia) and re-analyze 122 nearby (redshift < 0.11) SNe Ia using a new ``Multi-band Stretch method,'' which is a revised Stretch method (cf. Perlmutter et al. 1997; Goldhaber et al. 2001) extended to five bands. We find (i) our I-band template can fit about 90% of SNe Ia I-band light curves, (ii) relationships between luminosity, colours and stretch factors, (iii) possible sub-groups of SNe Ia, and (iv) the ratio of total to selective extinction R in other galaxies can be consistent with that in the Milky Way under the assumption that SNe Ia have diversity in their intrinsic colour. Based on these results, we discuss how to select subsets of SNe Ia to serve as good distance indicators for cosmology. We find two possibilities: one is to choose ``BV bluest'' (-0.14 < (B-V)_{max} <= -0.10) objects and the other is to use only SNe Ia which occur in E or S0 galaxies. Within these subsets, we find the root mean square (r.m.s.) of peak B-band magnitudes is 0.17 mag (``BV bluest'' sample) and 0.12 mag (E or S0 sample).

This paper presents spectroscopy of supernovae discovered in the first season of the Sloan Digital Sky Survey-II Supernova Survey. This program searches for and measures multi-band light curves of supernovae in the redshift range z = 0.05 - 0.4, complementing existing surveys at lower and higher redshifts. Our goal is to better characterize the supernova population, with a particular focus on SNe Ia, improving their utility as cosmological distance indicators and as probes of dark energy. Our supernova spectroscopy program features rapid-response observations using telescopes of a range of apertures, and provides confirmation of the supernova and host-galaxy types as well as precise redshifts. We describe here the target identification and prioritization, data reduction, redshift measurement, and classification of 129 SNe Ia, 16 spectroscopically probable SNe Ia, 7 SNe Ib/c, and 11 SNe II from the first season. We also describe our efforts to measure and remove the substantial host galaxy contamination existing in the majority of our SN spectra.

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

We report a comprehensive statistical analysis of the observational data of the cosmic evolution of supernova (SN) rate density, to derive constraints on cosmic star formation history and the nature of type Ia supernova (SN Ia) progenitor. We use all available information of magnitude, SN type, and redshift information of both type Ia and core-collapse (CC) SNe in GOODS and SDF, as well as SN Ia rate densities reported in the literature. Furthermore, we also add 157 SN candidates in the past Subaru/Suprime-Cam data that are newly reported here, to increase the statistics. We find that the current data set of SN rate density evolution already gives a meaningful constraint on the evolution of the cosmic star formation rate (SFR) at z <~ 1, though strong constraints cannot be derived for the delay time distribution (DTD) of SNe Ia. We derive a constraint of the evolutionary index of SFR density alpha ~ 3--4 [(1+z)^alpha at z <~ 1] with an evidence for a significant evolution of mean extinction of CC SNe [E(B-V) ~ 0.5 at z ~ 0.5 compared with ~ 0.2 at z = 0], which does not change significantly within a reasonable range of various DTD models. This result is nicely consistent with the systematic trend of alpha estimates based on galactic SFR indicators in different wavelengths (ultraviolet, H_alpha, and infrared), indicating that there is a strong evolution in mean extinction of star forming regions in galaxies at relatively low redshift range of z <~ 0.5. These results are obtained by a method that is completely independent of galaxy surveys, and especially, there is no detection limit about the host galaxy luminosity in our analysis, giving a strong constraint on the star formation activity in high-z dwarf galaxies or intergalactic space.

We present the properties of active galactic nuclei (AGN) selected by optical variability in the Subaru/XMM-Newton Deep Field (SXDF). Based on the locations of variable components and light curves, 211 optically variable AGN were reliably selected. We made three AGN samples; X-ray detected optically non-variable AGN (XA), X-ray detected optically variable AGN (XVA), and X-ray undetected optically variable AGN (VA). In the VA sample, we found a bimodal distribution of the ratio between the variable component flux and the host flux. One of these two components in the distribution, a class of AGN with a faint variable component $i'_{\rm{vari}}\sim25$ mag in bright host galaxies $i'\sim21$ mag, is not seen in the XVA sample. These AGN are expected to have low Eddington ratios if we naively consider a correlation between bulge luminosity and black hole mass. These galaxies have photometric redshifts $z_{\rm{photo}}\sim0.5$ and we infer that they are low-luminosity AGN with radiatively inefficient accretion flows (RIAFs). The properties of the XVA and VA objects and the differences from those of the XA objects can be explained within the unified scheme for AGN. Optical variability selection for AGN is an independent method and could provide a complementary AGN sample which even deep X-ray surveys have not found.

We present our survey for optically faint variable objects using multi-epoch (8-10 epochs over 2-4 years) $i'$-band imaging data obtained with Subaru Suprime-Cam over 0.918 deg$^2$ in the Subaru/XMM-Newton Deep Field (SXDF). We found 1040 optically variable objects by image subtraction for all the combinations of images at different epochs. This is the first statistical sample of variable objects at depths achieved with 8-10m class telescopes or HST. The detection limit for variable components is $i'_{\rm{vari}}\sim25.5$ mag. These variable objects were classified into variable stars, supernovae (SNe), and active galactic nuclei (AGN), based on the optical morphologies, magnitudes, colors, and optical-mid-infrared colors of the host objects, spatial offsets of variable components from the host objects, and light curves. Detection completeness was examined by simulating light curves for periodic and irregular variability. We detected optical variability for $36\pm2%$ ($51\pm3%$ for a bright sample with $i'<24.4$ mag) of X-ray sources in the field. Number densities of variable obejcts as functions of time intervals $\Delta{t}$ and variable component magnitudes $i'_{\rm{vari}}$ are obtained. Number densities of variable stars, SNe, and AGN are 120, 489, and 579 objects deg$^{-2}$, respectively. Bimodal distributions of variable stars in the color-magnitude diagrams indicate that the variable star sample consists of bright ($V\sim22$ mag) blue variable stars of the halo population and faint ($V\sim23.5$ mag) red variable stars of the disk population. There are a few candidates of RR Lyrae providing a possible number density of $\sim10^{-2}$ kpc$^{-3}$ at a distance of $>150$ kpc from the Galactic center.

The Sloan Digital Sky Survey-II Supernova Survey has identified a large number of new transient sources in a 300 sq. deg. region along the celestial equator during its first two seasons of a three-season campaign. Multi-band (ugriz) light curves were measured for most of the sources, which include solar system objects, Galactic variable stars, active galactic nuclei, supernovae (SNe), and other astronomical transients. The imaging survey is augmented by an extensive spectroscopic follow-up program to identify SNe, measure their redshifts, and study the physical conditions of the explosions and their environment through spectroscopic diagnostics. During the survey, light curves are rapidly evaluated to provide an initial photometric type of the SNe, and a selected sample of sources are targeted for spectroscopic observations. In the first two seasons, 476 sources were selected for spectroscopic observations, of which 403 were identified as SNe. For the Type Ia SNe, the main driver for the Survey, our photometric typing and targeting efficiency is 90%. Only 6% of the photometric SN Ia candidates were spectroscopically classified as non-SN Ia instead, and the remaining 4% resulted in low signal-to-noise, unclassified spectra. This paper describes the search algorithm and the software, and the real-time processing of the SDSS imaging data. We also present the details of the supernova candidate selection procedures and strategies for follow-up spectroscopic and imaging observations of the discovered sources.

We present a new measurement of the volumetric rate of Type Ia supernova up to a redshift of 1.7, using the Hubble Space Telescope (HST) GOODS data combined with an additional HST dataset covering the North GOODS field collected in 2004. We employ a novel technique that does not require spectroscopic data for identifying Type Ia supernovae (although spectroscopic measurements of redshifts are used for over half the sample); instead we employ a Bayesian approach using only photometric data to calculate the probability that an object is a Type Ia supernova. This Bayesian technique can easily be modified to incorporate improved priors on supernova properties, and it is well-suited for future high-statistics supernovae searches in which spectroscopic follow up of all candidates will be impractical. Here, the method is validated on both ground- and space-based supernova data having some spectroscopic follow up. We combine our volumetric rate measurements with low redshift supernova data, and fit to a number of possible models for the evolution of the Type Ia supernova rate as a function of redshift. The data do not distinguish between a flat rate at redshift > 0.5 and a previously proposed model, in which the Type Ia rate peaks at redshift >1 due to a significant delay from star-formation to the supernova explosion. Except for the highest redshifts, where the signal to noise ratio is generally too low to apply this technique, this approach yields smaller or comparable uncertainties than previous work.

Large samples of high-redshift supernovae (SNe) are potentially powerful probes of cosmic star formation, metal enrichment, and SN physics. We present initial results from a new deep SN survey, based on re-imaging in the R, i', z' bands, of the 0.25 deg2 Subaru Deep Field (SDF), with the 8.2-m Subaru telescope and Suprime-Cam. In a single new epoch consisting of two nights of observations, we have discovered 33 candidate SNe, down to a z'-band magnitude of 26.3 (AB). We have measured the photometric redshifts of the SN host galaxies, obtained Keck spectroscopic redshifts for 17 of the host galaxies, and classified the SNe using the Bayesian photometric algorithm of Poznanski et al. (2007) that relies on template matching. After correcting for biases in the classification, 55% of our sample consists of Type Ia supernovae and 45% of core-collapse SNe. The redshift distribution of the SNe Ia reaches z ~ 1.6, with a median of z ~ 1.2. The core-collapse SNe reach z ~ 1.0, with a median of z ~ 0.5. Our SN sample is comparable to the Hubble Space Telescope/GOODS sample both in size and redshift range. The redshift distributions of the SNe in the SDF and in GOODS are consistent, but there is a trend (which requires confirmation using a larger sample) for more high-z SNe Ia in the SDF. This trend is also apparent when comparing the SN Ia rates we derive to those based on GOODS data. Our results suggest a fairly constant rate at high redshift that could be tracking the star-formation rate. Additional epochs on this field, already being obtained, will enlarge our SN sample to the hundreds, and determine whether or not there is a decline in the SN Ia rate at z >~ 1.

The Sloan Digital Sky Survey-II (SDSS-II) has embarked on a multi-year project to identify and measure light curves for intermediate-redshift (0.05 < z < 0.35) Type Ia supernovae (SNe Ia) using repeated five-band (ugriz) imaging over an area of 300 sq. deg. The survey region is a stripe 2.5 degrees wide centered on the celestial equator in the Southern Galactic Cap that has been imaged numerous times in earlier years, enabling construction of a deep reference image for discovery of new objects. Supernova imaging observations are being acquired between 1 September and 30 November of 2005-7. During the first two seasons, each region was imaged on average every five nights. Spectroscopic follow-up observations to determine supernova type and redshift are carried out on a large number of telescopes. In its first two three-month seasons, the survey has discovered and measured light curves for 327 spectroscopically confirmed SNe Ia, 30 probable SNe Ia, 14 confirmed SNe Ib/c, 32 confirmed SNe II, plus a large number of photometrically identified SNe Ia, 94 of which have host-galaxy spectra taken so far. This paper provides an overview of the project and briefly describes the observations completed during the first two seasons of operation.

We present extensive ugrizYHJK photometry and optical spectroscopy of SN 2005gj obtained by the SDSS-II and CSP Supernova Projects, which give excellent coverage during the first 150 days after the time of explosion. These data show that SN 2005gj is the second clear case, after SN 2002ic, of a thermonuclear explosion in a dense circumstellar environment. Both the presence of singly and doubly ionized iron-peak elements (FeIII and weak SII, SiII) near maximum light as well as the spectral evolution show that SN 2002ic-like events are Type Ia explosions. Independent evidence comes from the exponential decay in luminosity of SN 2005gj, pointing to an exponential density distribution of the ejecta. The interaction of the supernova ejecta with the dense circumstellar medium is stronger than in SN 2002ic: (1) the supernova lines are weaker; (2) the Balmer emission lines are more luminous; and (3) the bolometric luminosity is higher close to maximum light. The velocity evolution of the Halpha components suggest that the CSM around SN 2005gj is clumpy and it has a flatter density distribution compared with the steady wind solution, in agreement with SN 2002ic. An early X-ray observation with Chandra gives an upper-limit on the mass loss rate from the companion of < 2x10^{-4} Msun/yr.

Galactic extinction is tested using galaxy number counts at low Galactic latitude obtained from five band photometry of the Sloan Digital Sky Survey. The spatial variation of galaxy number counts for low extinction regions of $E(B-V) < 0.15$ is consistent with the all-sky reddening map of Schlegel, Finkbeiner, & Davis (1998) and the standard extinction law. For higher extinction regions of $E(B-V) > 0.15$, however, the map of Schlegel et al.(1998) overestimates the reddening by a factor up to 1.4, which is likely ascribed to the departure from proportionality of reddening to infrared emissivity of dust. This result is consistent with the analysis of Arce & Goodman (1999) for the Taurus dark cloud complex.

We carried out a wide-field V, I imaging survey of the Local Group dwarf spheroidal galaxy Leo II using the Subaru Prime Focus Camera on the 8.2-m Subaru Telescope. The survey covered an area of 26.67 x 26.67 arcmin^2, far beyond the tidal radius of Leo II (8.63 arcmin), down to the limiting magnitude of V ~26, which is roughly 1 mag deeper than the turn-off point of the main sequence stars of Leo II. Radial number density profiles of bright and faint red giant branch (RGB) stars were found to change their slopes at around the tidal radius, and extend beyond the tidal radius with shallower slopes. A smoothed surface brightness map of Leo II suggests the existence of a small substructure of globular cluster luminosity beyond the tidal radius. We investigated the properties of the stellar population by means of the color-magnitude diagram. The horizontal branch (HB) morphology index shows a radial gradient in which red HB stars are more concentrated than blue HB stars, which is common to many Local Group dwarf spheroidal galaxies. The color distribution of RGB stars around the mean RGB sequence shows a larger dispersion at the center than in the outskirts, indicating a mixture of stellar populations at the center and a more homogeneous population in the outskirts. Based on the age estimation using subgiant branch (SGB) stars, we found that although the major star formation took place ~8 Gyr ago, a considerable stellar population younger than 8 Gyr is found at the center; such a younger population is insignificant in the outskirts.

We present a catalogue of morphologically classified bright galaxies in the north equatorial stripe (230 deg$^2$) derived from the Third Data Release of the Sloan Digital Sky Survey (SDSS). Morphological classification is performed by visual inspection of images in the $g$ band. The catalogue contains 2253 galaxies complete to a magnitude limit of $r=16$ after Galactic extinction correction, selected from 2658 objects that are judged as extended in the photometric catalogue in the same magnitude limit. 1866 galaxies in our catalogue have spectroscopic information. A brief statistical analysis is presented for the frequency of morphological types and mean colours in the catalogue. A visual inspection of the images reveals that the rate of interacting galaxies in the local Universe is approximately 1.5% in the $r\le16$ sample. A verification is made for the photometric catalogue generated by the SDSS, especially as to its bright end completeness.

We present Keck diffraction limited H-band photometry of a z~1.3 Type Ia supernova (SN) candidate, first identified in a Hubble Space Telescope (HST) search for SNe in massive high redshift galaxy clusters. The adaptive optics (AO) data were obtained with the Laser Guide Star facility during four observing runs from September to November 2005. In the analysis of data from the observing run nearest to maximum SN brightness, the SN was found to have a magnitude H=23.9 +/- 0.14 (Vega). We present the H-band (approximately rest-frame R) light curve and provide a detailed analysis of the AO photometric uncertainties. By constraining the aperture correction with a nearby (4" separation) star we achieve 0.14 magnitude photometric precision, despite the spatially varying AO PSF.

The Coma cluster is the richest and most compact of the nearby clusters, yet there is growing evidence that its formation is still on-going. With a new multi-slit imaging spectroscopy technique pioneered at the 8.2 m Subaru telescope and FOCAS, we have detected and measured the line-of-sight velocities of 37 intracluster planetary nebulae associated with the diffuse stellar population of stars in the Coma cluster core, at 100 Mpc distance. We detect clear velocity substructures within a 6 arcmin diameter field. A substructure is present at ~5000 km/s, probably from in-fall of a galaxy group, while the main intracluster stellar component is centered around ~6500 km/s, ~700 km/s offset from the nearby cD galaxy NGC 4874. The kinematics and morphology of the intracluster stars show that the cluster core is in a highly dynamically evolving state. In combination with galaxy redshift and X-ray data this argues strongly that the cluster is currently in the midst of a subcluster merger, where the NGC 4874 subcluster core may still be self-bound, while the NGC 4889 subcluster core has probably dissolved. The NGC 4889 subcluster is likely to have fallen into Coma from the eastern A2199 filament, in a direction nearly in the plane of the sky, meeting the NGC 4874 subcluster arriving from the west. The two inner subcluster cores are presently beyond their first and second close passage, during which the elongated distribution of diffuse light has been created. We predict the kinematic signature expected in this scenario, and argue that the extended western X-ray arc recently discovered traces the arc shock generated by the collision between the two subcluster gas halos. Any preexisting cooling core region would have been heated by the subcluster collision.

We present extensive u'g'r'i'BVRIYJHKs photometry and optical spectroscopy of SN 2005hk. These data reveal that SN 2005hk was nearly identical in its observed properties to SN 2002cx, which has been called ``the most peculiar known type Ia supernova.'' Both supernovae exhibited high ionization SN 1991T-like pre-maximum spectra, yet low peak luminosities like SN 1991bg. The spectra reveal that SN 2005hk, like SN 2002cx, exhibited expansion velocities that were roughly half those of typical type Ia supernovae. The R and I light curves of both supernovae were also peculiar in not displaying the secondary maximum observed for normal type Ia supernovae. Our YJH photometry of SN 2005hk reveals the same peculiarity in the near-infrared. By combining our optical and near-infrared photometry of SN 2005hk with published ultraviolet light curves obtained with the Swift satellite, we are able to construct a bolometric light curve from ~10 days before to ~60 days after B maximum. The shape and unusually low peak luminosity of this light curve, plus the low expansion velocities and absence of a secondary maximum at red and near-infrared wavelengths, are all in reasonable agreement with model calculations of a 3D deflagration which produces ~0.25 M_sun of 56Ni.

The Coma cluster is the richest and most compact of the nearby clusters, yet there is growing evidence that its formation is still on-going. A sensitive probe of this evolution is the dynamics of intracluster stars, which are unbound from galaxies while the cluster forms, according to cosmological simulations. With a new multi-slit imaging spectroscopy technique pioneered at the 8.2 m Subaru telescope and FOCAS, we can now detect and measure the line-of-sight velocities of the intracluster planetary nebulae which are associated with the diffuse stellar population of stars, at 100 Mpc distance. We detect significant velocity substructures within a 6 arcmin diameter field, centred on the Coma X-ray cluster emission. One substructure ispresent at ~5000 km/s, probably from infall of a galaxy group, while the main intracluster stellar component moves at ~ 6500 km/s. Hence the ICPNs associated with the diffuse light at the position of the MSIS field are not bound to the nearby cD galaxy NGC 4874, whose radial velocity is ~ 700 km/s greater.

We present our activity on making the Subaru Data Archive accessible through the Japanese Virtual Observatory (JVO) system. There are a lot of demand to use the archived Subaru data from various fields of astronomers. To be used by those who are not familiar with the way to reduce the Subaru data, the data reduction should be made before providing for them, or at least it should be easily done without precise knowledge about instrument's characteristic and environment where data are taken. For those purposes, data quality assessment system NAQATA is developed, which is presented in this meeting by Nakata et al. (2006), and the science-ready image data are provided for some of the SuprimeCam data at SMOKA data service which is presented by Enoki et al (2006). JVO portal will provide the way to access the reduce data, and also provides the way to reduce from raw data with very few efforts through the user-friendly web browser I/F. To provide such a CPU-intensive service, we have developed a GRID computing system. The architecture of this Subaru Data and Analysis system are discussed.

We report the discovery of a cluster-scale lensed quasar, SDSS J1029+2623, selected from the Sloan Digital Sky Survey. The lens system exhibits two lensed images of a quasar at z_s=2.197. The image separation of 22.5" makes it the largest separation lensed quasar discovered to date. The similarity of the optical spectra and the radio loudnesses of the two components support the lensing hypothesis. Images of the field show a cluster of galaxies at z_l~0.55 that is responsible for the large image separation. The lensed images and the cluster light center are not collinear, which implies that the lensing cluster has a complex structure.

We present the design and implementation of the Japanese Virtual Observatory (JVO) system. JVO is a portal site to various kinds of astronomical resources distributed all over the world. We have developed five components for constructing the portal: (1) registry, (2) data service, (3) workflow system, (4) data analysis service (5) portal GUI. Registry services are used for publishing and searching data services in the VO, and they are constructed using an OAI-PMH metadata harvesting protocol and a SOAP web service protocol so that VO standard architecture is applied. Data services are developed based on the Astronomical Data Query Language (ADQL) which is an international VO standard and an extension of the standard SQL. The toolkit for building the ADQL-based service is released to the public on the JVO web site. The toolkit also provides the protocol translation from a Simple Image Access Protocol (SIAP) to ADQL protocol, so that both the VO standard service can be constructed using our toolkit. In order to federate the distributed databases and analysis services, we have designed a workflow language which is described in XML and developed execution system of the workflow. We have succeeded to connect to a hundred of data resources of the world as of April 2006. We have applied this system to the study of QSO environment by federating a QSO database, a Subaru Suprim-Cam database, and some analysis services such a SExtractor and HyperZ web services. These experiences are described is this paper.

We devise improved photometric parameters for the morphological classification of galaxies using a bright sample from the First Data Release of the Sloan Digital Sky Survey. In addition to using an elliptical aperture concentration index for classification, we introduce a new texture parameter, coarseness, which quantifies deviations from smooth galaxy isophotes. The elliptical aperture concentration index produces morphological classifications that are in appreciably better agreement with visual classifications than those based on circular apertures. With the addition of the coarseness parameter,the success rate of classifying galaxies into early and late types increases to \~88% with respect to the reference visual classification. A reasonably high success rate ~68% is also attained in classifying galaxies into three types, early-type galaxies (E+S0), early- (Sa+Sb) and late- (Sc+Sdm+Im) type spiral galaxies.

Radial velocities of 40 intracluster planetary nebulae (ICPNe) in the Virgo cluster were obtained with the new multi-fiber FLAMES spectrograph on UT2 at VLT. For the first time, the lambda 4959 AA line of the [OIII] doublet is seen in a large fraction (50%) of ICPNe spectra, and a large fraction of the photometric candidates with m(5007) <\= 27.2 is spectroscopically confirmed. ICPNe with the velocity dispersion of the Virgo cluster are found in our CORE field 1 deg from M87. These may have originated from tidal mass loss of smaller galaxies in the M87 subcluster halo. In a field 0.25 deg from M87, we see an extended stellar halo of M87 in approximate dynamical equilibrium, but with few ICPNe. Finally, in a field near M84/M86, the ICPNe velocities are highly correlated with the galaxy velocities, showing that any well-mixed intracluster population is yet to form. Overall, the measured velocity distributions confirm the non-uniform dynamical structure and on-going assembly of the Virgo cluster.(Based on data collected with the FLAMES spectrograph at the UT2 of the VLT at Cerro Paranal, Chile, operated by ESO, during observing run 71.B-0147(A))

[OIII] lambda 5007AA emission lines of 16 intracluster planetary nebulae candidates in the Coma cluster were detected with a Multi-Slit Imaging Spectroscopy (MSIS) technique using FOCAS on the Subaru telescope. The identification of these faint emission sources as PNe is supported by (i) their point-like flux distribution in both space and wavelength, with tight limits on the continuum flux; (ii) the identification of the second [OIII] lambda 4959 line in the only object at high enough velocity that this line too falls into the filter bandpass; (iii) emission line fluxes consistent with PNe at 100 Mpc distance, in the range 2.8 x 10^{-19} - 1.7 x 10^{-18} erg/s/cm^2; and (iv) a narrow velocity distribution approximately centered on the systemic velocity of the Coma cluster. Comparing with the velocities of galaxies in our field, we conclude that the great majority of these candidates would be intracluster PNe, free floating in the Coma cluster core. Their velocity dispersion is ~760 km/s, and their mean velocity is lower than that of the galaxies. The velocity distribution suggests that the intracluster stellar population has different dynamics from the galaxies in the Coma cluster core.

We report a serendipitous discovery of six very low-luminosity active galactic nuclei (AGNs) only by optical variability in one-month baseline. The detected flux variability is ~ 1-5% of the total luminosity of host galaxies. Careful subtraction of host galaxy components in nuclear regions indicates that the fractional variability (Delta F / F) of the nuclei is of order unity. At least one of them is showing a compelling flaring activity within just a few days, which appears to be quite different from previously known AGN variability. We obtained spectroscopic data for the one showing the largest flare and confirmed that it is in fact an AGN at z = 0.33 with an estimated black hole mass of ~10^8 M_sun. As a possible interpretation, we suggest that these activities are coming from the region around the black hole event horizon, which is physically similar to the recently discovered near-infrared flares of our Galactic nucleus. It is indicated that our Galaxy is not special, and that surprisingly rapid flaring activity in optical/near-infrared bands may be commonly hidden in nuclei of apparently normal galaxies with low Eddington ratios, in contrast to the variability of well-studied luminous AGNs or quasars.

We report the discovery of a probable large scale structure composed of many galaxy clumps around the known twin clusters at z=1.26 and z=1.27 in the Lynx region. Our analysis is based on deep, panoramic, and multi-colour imaging (26.4'x24.1') in VRi'z'-bands with the Suprime-Cam on the 8.2m Subaru telescope. This unique, deep and wide-field imaging data-set allows us for the first time to map out the galaxy distribution in the highest redshift supercluster known. We apply a photometric redshift technique to extract plausible cluster members at z~1.27 down to i'=26.15 (5sigma) corresponding to \~M*+2.5 at this redshift. From the 2-D distribution of these photometrically selected galaxies, we newly identify seven candidates of galaxy groups or clusters where the surface density of red galaxies is significantly high (>5sigma), in addition to the two known clusters. These candidates show clear red colour-magnitude sequences consistent with a passive evolution model, which suggests the existence of additional high density regions around the Lynx superclusters.

Currently two query languages are defined as standards for the Virtual Observatory (VO). Astronomical Data Query Language (ADQL) is used for catalog data query and Simple Image Access Protocol (SIAP) is for image data query. As a result, when we query each data service, we need to know in advance which language is supported and then construct a query language accordingly. The construct of SIAP is simple, but they have a limited capability. For example, there is no way to specify multiple regions in one query, and it is difficult to specify complex query conditions. In this paper, we propose a unified query language for any kind of astronomical database on the basis of SQL99. SQL is a query language optimized for a table data, so to apply the SQL to the image and spectrum data set, the data structure need to be mapped to a table like structure. We present specification of this query language and an example of the architecture for the database system.

We present VLT FORS1 and FORS2 spectra of 39 candidate high-redshift supernovae that were discovered as part of a cosmological study using Type Ia supernovae (SNe Ia) over a wide range of redshifts. From the spectra alone, 20 candidates are spectrally classified as SNe Ia with redshifts ranging from z=0.212 to z=1.181. Of the remaining 19 candidates, 1 might be a Type II supernova and 11 exhibit broad supernova-like spectral features and/or have supernova-like light curves. The candidates were discovered in 8 separate ground-based searches. In those searches in which SNe Ia at z ~ 0.5 were targeted, over 80% of the observed candidates were spectrally classified as SNe Ia. In those searches in which SNe Ia with z > 1 were targeted, 4 candidates with z > 1 were spectrally classified as SNe Ia and later followed with ground and space based observatories. We present the spectra of all candidates, including those that could not be spectrally classified as supernova.

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

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

We investigate the clustering properties of 2,600 Lyman Break Galaxies (LBGs) at z=3.5-5.2 in two large blank fields, the Subaru Deep Field and the Subaru/XMM Deep Field (600arcmin^2 each). The angular correlation functions of these LBGs show a clear clustering at both z~4 and 5. The correlation lengths are r_0= 4.1^{+0.2}_{-0.2} and 5.9^{+1.3}_{-1.7} h_{100}^{-1} Mpc (r_0= 5.1^{+1.0}_{-1.1} and 5.9^{+1.3}_{-1.7} h_{100}^{-1} Mpc) for all the detected LBGs (for L>L* LBGs) at z~4 and 5, respectively. These correlation lengths correspond to galaxy-dark matter biases of b_g= 2.9^{+0.1}_{-0.1} and 4.6^{+0.9}_{-1.2} (b_g=3.5^{+0.6}_{-0.7} and 4.6^{+0.9}_{-1.2}), for all the detected LBGs (for L>L^* LBGs) at z~4 and 5, respectively. These results, combined with estimates for z~3 LBGs in the literature, show that the correlation length of L>L^* LBGs is almost constant, ~5 h_{100}^{-1} Mpc, over z~3-5, while the bias monotonically increases with redshift at z>3. We also find that for LBGs at z~4 the clustering amplitude increases with UV-continuum luminosity and with the amount of dust extinction. We estimate the mass of dark halos hosting various kinds of high-z galaxies including LBGs with the analytic model given by Sheth & Tormen (1999). We find that the typical mass of dark halos hosting L>L^* LBGs is about 1x10^{12} h_{70}^{-1}Msol over z~3-5, which is comparable to that of the Milky Way Galaxy. A single dark halo with ~10^{12} h_{70}^{-1} Msol is found to host 0.1-0.3 LBG on average but host about four K-band selected galaxies.

(abridged) We investigate photometric properties of Lyman Break Galaxies (LBGs) at z=3.5-5.2 based on large samples of 2,600 LBGs detected in deep (i'~27) and wide-field (1,200 arcmin^2) images taken in the Subaru Deep Field (SDF) and the Subaru/XMM Deep Field (SXDF). The selection criteria for the LBG samples are examined with 85 spectroscopically identified objects and by Monte Carlo simulations. We find in the luminosity functions of LBGs (i) that the number density of bright galaxies (M_{1700}<-22; corresponding to SFR_{corr}>100 Msolar yr^{-1}) decreases significantly from z=4 to 5 and (ii) that the faint-end slope of the luminosity function may become steeper towards higher redshifts. We estimate dust extinction of z=4 LBGs with M<M^* from UV slopes, and obtain E(B-V)=0.15+/-0.03 as the mean value. The dust extinction remains constant with apparent luminosity, but increases with intrinsic luminosity. We find no evolution in dust extinction between LBGs at z=3 and 4. We investigate the evolution of UV-luminosity density at 1700A, rho, and find that rho does not significantly change from z=3 to z=5, i.e., rho(z=4)/rho(z=3)=1.0+/-0.2 and rho(z=5)/rho(z=3)=0.8+/-0.4, thus the cosmic star-formation rate (SFR) density remains constant. We find that the stellar mass density estimated from the cosmic SFR is consistent with those derived directly from the stellar mass function at z=0-1, but exceeds those at z~3 by a factor of 3. We find that the ratio of the UV-luminosity density of Ly-a emitters (LAEs) to that of LBGs is ~60% at z=5, and thus about a half of the star formation at z=5 probably occurs in LAEs. We obtain a constraint on the escape fraction of UV-ionizing photons produced by LBGs, f_{esc}>~0.13.

We study the spatial variation of galaxy number counts using five band photometric images from the Sloan Digital Sky Survey. The spatial variation of this sample of 46 million galaxies collected from 2200 sq. degrees can be understood as the combination of Galactic extinction and large-scale clustering. With the use of the reddening map of Schlegel, Finkbeiner & Davis (1998), the standard extinction law is verified for the colour bands from u to z within 5% in the region of small extinction values, E(B-V)<0.15. The residual spatial variations of the number counts suggests that the error of global calibration for SDSS photometry is smaller than 0.02 mag.

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

We report measurements of $\Omega_M$, $\Omega_\Lambda$, and w from eleven supernovae at z=0.36-0.86 with high-quality lightcurves measured using WFPC-2 on the HST. This is an independent set of high-redshift supernovae that confirms previous supernova evidence for an accelerating Universe. Combined with earlier Supernova Cosmology Project data, the new supernovae yield a flat-universe measurement of the mass density $\Omega_M=0.25^{+0.07}_{-0.06}$ (statistical) $\pm0.04$ (identified systematics), or equivalently, a cosmological constant of $\Omega_\Lambda=0.75^{+0.06}_{-0.07}$ (statistical) $\pm0.04$ (identified systematics). When the supernova results are combined with independent flat-universe measurements of $\Omega_M$ from CMB and galaxy redshift distortion data, they provide a measurement of $w=-1.05^{+0.15}_{-0.20}$ (statistical) $\pm0.09$ (identified systematic), if w is assumed to be constant in time. The new data offer greatly improved color measurements of the high-redshift supernovae, and hence improved host-galaxy extinction estimates. These extinction measurements show no anomalous negative E(B-V) at high redshift. The precision of the measurements is such that it is possible to perform a host-galaxy extinction correction directly for individual supernovae without any assumptions or priors on the parent E(B-V) distribution. Our cosmological fits using full extinction corrections confirm that dark energy is required with $P(\Omega_\Lambda>0)>0.99$, a result consistent with previous and current supernova analyses which rely upon the identification of a low-extinction subset or prior assumptions concerning the intrinsic extinction distribution.

We present measurements of the cosmic shear correlation in the shapes of galaxies in the Suprime-Cam 2.1 deg^2 R_c-band imaging data. As an estimator of the shear correlation originated from the gravitational lensing, we adopt the aperture mass variance. We detect a non-zero E mode variance on scales between 2 and 40arcmin. We also detect a small but non-zero B mode variance on scales larger than 5arcmin. We compare the measured E mode variance to the model predictions in CDM cosmologies using maximum likelihood analysis. A four-dimensional space is explored, which examines sigma_8, Omega_m, Gamma and zs (a mean redshift of galaxies). We include three possible sources of error: statistical noise, the cosmic variance estimated using numerical experiments, and a residual systematic effect estimated from the B mode variance. We derive joint constraints on two parameters by marginalizing over the two remaining parameters. We obtain an upper limit of Gamma<0.5 for zs>0.9 (68% confidence). For a prior Gamma\in[0.1,0.4] and zs\in[0.6,1.4], we find sigma_8=(0.50_{-0.16}^{+0.35})Omega_m^{-0.37} for flat cosmologies and sigma_8=(0.51_{-0.16}^{+0.29})Omega_m^{-0.34}$ for open cosmologies (95% confidence). If we take the currently popular LCDM model, we obtain a one-dimensional confidence interval on sigma_8 for the 95.4% level, 0.62<\sigma_8<1.32 for zs\in[0.6,1.4]. Information on the redshift distribution of galaxies is key to obtaining a correct cosmological constraint. An independent constraint on Gamma from other observations is useful to tighten the constraint.

It has been suggested that cosmic strings produced at a phase transition in the early universe can be the origin of the extremely high energy cosmic rays (EHCR) observed by AGASA above 10^20 eV. Superheavy cosmic strings with linear mass density of 10^22 g/cm can be indirectly observed through the gravitational lensing effect the distant galaxies. The lensing effect by a long straight object can be characterized by a line of double galaxies or quasars with angular separation of about 5 arcsec. We have searched for aligned double objects from the archived data taken by the Subaru Prime Focus Camera (Suprime-Cam). The Suprime-Cam has a great advantage in observing the wide field of view (30x30 arcmin^2) with high sensitivity (R<26 400s exposure), so it is suitable for this research. In this paper, we describe the result of simulation study for developing the method of searching the objects lensed by cosmic strings, and present the observational result obtained by this method.

We report the discovery of a large-scale structure of Lyman alpha emitters (LAEs) at z=4.86 based on wide-field imaging with the prime-focus camera (Suprime-Cam) on the Subaru telescope. We observed a 25' x 45' area of the Subaru Deep Field in a narrow band (NB, lambdaC=7126 A and FWHM=73 A) together with R and i'. We isolate from these data 43 LAE candidates down to NB=25.5 mag using color criteria. Follow-up spectroscopy of five candidates suggests the contamination by low-z objects to be ~ 20%. We find that the LAE candidates are clustered in an elongated region on the sky of 20 Mpc in width and 50 Mpc in length at z=4.86, which is comparable in size to present-day large-scale structures (we adopt H0=70 km/s/Mpc, Omega0=0.3, lambda0=0.7). This elongated region includes a circular region of 12 Mpc radius of higher surface overdensity (delta=2), which may be the progenitor of a cluster of galaxies. Assuming this circular region to be a sphere with a spatial overdensity of 2, we compare our observation with predictions by Cold Dark Matter models. We find that an Omega0=0.3 flat model with sigma8=0.9 predicts the number of such spheres consistent with the observed number (one sphere in our survey volume) if the bias parameter of LAEs is b \simeq 6. This value suggests that the typical mass of dark haloes hosting LAEs at z \simeq 5 is of the order of 10^12 Msolar. Such a large mass poses an interesting question about the nature of LAEs.

The morphological dependence of the luminosity function is studied using a sample containing approximately 1500 bright galaxies classified into Hubble types by visual inspections for a homogeneous sample obtained from the Sloan Digital Sky Survey (SDSS) northern equatorial stripes. Early-type galaxies are shown to have a characteristic magnitude by 0.45 mag brighter than spiral galaxies in the $r^{\ast}$ band, consistent with the `universal characteristic luminosity' in the $B$ band. The shape of the luminosity function differs rather little among different morphological types: we do not see any symptoms of the sharp decline in the faint end for the luminosity function for early-type galaxies at least 2 mag fainter than the characteristic magnitude, although the faint end behaviour shows a slight decline ($\alpha\lsim -1$) compared with the total sample. We also show that a rather flat faint end slope for early-type galaxies is not due to an increasing mixture of the dwarf galaxies which have softer cores. This means that there are numerous faint early-type galaxies with highly concentrated cores.

(Abridged) We present new results on the Hubble diagram of distant type Ia supernovae (SNe Ia) segregated according to the type of host galaxy. This makes it possible to check earlier evidence for a cosmological constant by explicitly comparing SNe residing in galaxies likely to contain negligible dust with the larger sample. The cosmological parameters derived from these SNe Ia hosted by presumed dust-free early-type galaxies supports earlier claims for a cosmological constant, which we demonstrate at 5 sigma significance, and the internal extinction implied is small even for late-type systems (A_B<0.2). Thus, our data demonstrate that host galaxy extinction is unlikely to systematically dim distant SNe Ia in a manner that would produce a spurious cosmological constant. We classify the host galaxies of 39 distant SNe discovered by the Supernova Cosmology Project (SCP) using the combination of HST STIS imaging, Keck ESI spectroscopy and ground-based broad-band photometry. We compare with a low-redshift sample of 25 SNe Ia. The scatter observed in the SNe Ia Hubble diagrams correlates closely with host galaxy morphology. We find the scatter in the SNe Ia Hubble diagram is smallest for SNe occurring in early-type hosts and largest for those occurring in late-type galaxies. Moreover, SNe residing in early-type hosts appear only ~0.14+/-0.09 mag brighter in their light-curve-width-corrected luminosity than those in late-type hosts, implying only a modest amount of dust extinction even in the late-type systems.

We have built an 80-mega pixels (10240 X 8192) mosaic CCD camera, called Suprime-Cam, for the wide-field prime focus of the 8.2 m Subaru telescope. Suprime-Cam covers a field of view 34' X 27', a unique facility among the the 8-10 m class telescopes, with a resolution of 0.''202 per pixel. The focal plane consists of ten high-resistivity 2k X 4k CCDs developed by MIT Lincoln Laboratory, which are cooled by a large stirling-cycle cooler. The CCD readout electronics was designed to be scalable, which allows the multiple read-out of tens of CCDs. It takes 50 seconds to readout entire arrays. We designed a filter-exchange mechanism of the jukebox type that can hold up to ten large filters (205 X 170 X 15 mm^3). The wide-field corrector is basically a three-lens Wynne-type, but has a new type of atmospheric dispersion corrector. The corrector provides a flat focal plane and an un-vignetted field of view of 30' in diameter. The achieved co-planarity of the focal array mosaic is smaller than 30 micron peak-to-peak, which realizes mostly the seeing limited image over the entire field. The median seeing in the Ic-band, measured over one year and a half, is 0.''61. The PSF anisotropy in Suprime-Cam images, estimated by stellar ellipticities, is about 2% under this median seeing condition. At the time of commissioning, Suprime-Cam had the largest survey speed, which is defined as the field of view multiplied by the primary mirror area of the telescope, among those cameras built for sub-arcsecond imaging.

We have identified 38 candidates of intracluster planetary nebulae (ICPNe) in a $34'\times27'$ field in the core of the Virgo cluster based on the Suprime-Cam imaging through two narrow-band filters centered at the redshifted wavelengths of the [OIII]$\lambda=5007\AA$ and the Halpha \lambda=6563\AA lines. Broad-band images in V and R bands are used to check for any emissions in the adjacent continuum. We describe the method briefly and present the list of intracluster planetary nebulae candidates, together with their finding charts. The ICPN candidates show a highly inhomogeneous distribution, which may suggest an association with the M86-M84 subcluster. Fraction of diffuse intracluster light with respect to total light in galaxies is estimated to be about 10%, leading to an estimate of about 20% for the baryon fraction. Spectroscopic follow up and a wider survey are critical to reveal the nature of intracluster stellar population.

We have identified intracluster planetary nebulae in a Virgo cluster core field by imaging with the Subaru Suprime-Cam through two narrow band filters centered at the redshifted wavelengths of the [OIII] $\lambda=5007$ \AA and the H$\alpha$ $\lambda = 6563$ \AA lines; broad-band images in V and R were acquired to check for emission in the adjacent continuum. Emission line objects in Virgo are then selected from the two-color diagram ([OIII] -- Halpha) vs. ([OIII] -- (V+R)), which is calibrated using PNe in M84 (Jacoby et al. 1990). Using both [OIII] and Halpha allows us to distinguish bona-fide planetary nebulae from high redshift emission-line galaxies at the bright end of the [OIII] luminosity function. Spectroscopic observations of a subsample of these objects were made at the TNG and at the VLT, in a region around M84 and in an intracluster field respectively. The observations confirm the efficiency of the combined [OIII]+H$\alpha$ imaging to identify true PNe. We also obtained the first spectrum of an intracluster PN which shows the [OIII] doublet with S/N > 10 and its Halpha emission. From the results based on the spectroscopic follow-up, we derive a lower limit to the fraction of the Virgo cluster light contributed by the intracluster stars at the surveyed position in the cluster core: it amounts to 10%.

We study the evolution of elliptical galaxies at z >~ 1 on the basis of a sample of 247 Extremely Red Objects (EROs) with R-Ks >= 3.35 (AB) and Ks <= 22.1 (AB) constructed from BVRi'z'JHKs multicolor data of a 114 arcmin2 area in the Subaru/XMM Deep Survey Field. By fitting template spectra of old ellipticals (OEs) and young, dusty starbursts (DSs) to the multicolor data, we classify EROs into these two classes and estimate their redshifts. We find that 58% of the EROs in our sample belong to the OE class and that 24% of the OEs are fit by a spectrum having a disk component with the B-band bulge-to-total luminosity ratio of <= 0.9. We compare the observed surface density of OEs with predictions from passive evolution models. We then derive rest-frame B-band luminosity functions (LFs) of OEs in our sample at z=1-1.5 and 1.5-2.5. We find that the LF at z=1-1.5 roughly agrees with the LF of local ellipticals if a dimming of 1.3 mag from z=1.25 to the present epoch is assumed. On the other hand, the amplitude of the LF at z=1.5-2.5 is found to be lower than that of the local LF by a factor of ~ 3 over the whole range of magnitude observed. Taking account of a strong decrease in the number density of morphologically classified early-type galaxies at z >~ 1.5 found by several authors, we conclude that the majority of ellipticals seen at present have not established either a red color or a smooth 1/4-law profile before z ~ 1.5$. (abridged)

We report the first result of weak gravitational lensing survey on a 2.1 sq deg Rc-band image taken with a wide field camera (Suprime-Cam) on the prime focus of 8.2 m Subaru Telescope. The weak lensing mass reconstruction is applied to the data to search for dark matter halos of cluster scale; M >= 10^14 solar mass. The reconstructed convergence field is divided by 1-sigma noise to obtain the signal-to-noise ratio map (S/N-map) of the detection. Local maxima and minima are searched on the S/N-map and the probability distribution function (PDF) of the peaks are created to compare with model predictions. We found excess over noise PDF created from the randomized realization on both positive and negative sides. Negative peaks imply the presence of voids in the dark matter distribution and this is the first report of the detection. Positive peaks, on the other hand, represent the dark matter halos and the number count of the halos on the 2.1 sq deg image is 4.9 +- 2.3 for S/N > 5 where the Gaussian smoothing radius of the convergence map is 1'. The result is consistent with the prediction assuming the Press-Schechter mass function and the NFW halo profile under the cluster normalized CDM cosmology. Although the present statistics is not enough due to the limited field of view, this work demonstrates that dark matter halo count via weak lensing is a promising way to test the paradigm of structure formation and cosmological model.

We present the results of a survey for emission-line objects based on optical intermediate-band ($\lambda_{\rm c}$ = 5736 \AA ~ and $\Delta\lambda$ = 280 \AA) and broad-band ($B$, $V$, $R$, and $i^\prime$) observations of the Subaru/XMM-Newton Deep Field on the 8.2 m Subaru telescope with the Subaru Prime Focus Camera, Suprime-Cam. All the data were obtained during the guaranteed time observations of the Suprime-Cam instrument. The intermediate-band image covered a sky area with 10\minpoint62 $\times 12\minpoint40 \approx 132$ arcmin$^2$ in the Subaru/XMM-Newton Deep Field (Ouchi et al.). Using this image, we have found 23 emission-line sources whose observed emission-line equivalent widths are greater than 250 \AA. Their optical multicolor properties indicate that six emission-line sources are Ly$\alpha$ emitters at $z \approx$ 3.7 ($\Delta z \approx 0.22$). They are either intense starburst galaxies or active galactic nuclei like quasars at $z \approx$ 3.7. Two more emission-line sources may also be Ly$\alpha$ emitters at $z \approx$ 3.7 although their multicolor properties are marginal. Among the remaining 15 emission-line objects, eight objects appear strong emission-line galaxies at lower redshift; e.g., [O {\sc ii}] $\lambda$3727 emitters at $z \approx 0.54$, H$\beta$ at $z \approx 0.18$, or [O {\sc iii}]$\lambda$5007 emitters at $z \approx 0.15$. The remaining seven objects are unclassified because they are too faint to be detected in broad-band images. We discuss observational properties of these strong emission-line sources. In particular, our data allow us to estimate the star formation density at $z \approx 3.7$ for the first time.

We report on early results from deep and wide-field narrow-band imaging on a 543 arcmin^2 area of the Subaru Deep Field. We find 87 Lya emitters (LAEs) at z=4.86(+/-0.03) which are photometrically selected by a combination of two broad bands (R and i') and one narrow band (NB711; 7126A FWHM=73A). We derive the luminosity functions (LFs) of the LAEs at Lya luminosity and at UV-continuum (1700A) luminosity. The LFs show little evolution between z=3.4 and z=4.86 either in Lya or UV-continuum emission. The amplitude of the LAE LF tends to decline at the bright magnitudes more rapidly than that of the LBG LF at similar redshifts. We calculate the angular correlation function of our LAEs up to ~15 arcmin separations. The angular correlation w(theta) is found to increase with decreasing angular separations, showing a clear signal of clustering. It is also found that the distribution of LAEs shows a large density gradient with a scale of >~15 arcmin, which would indicate the existence of a large-scale structure of LAEs on >~20h^(-1) Mpc scales.We fit the observed correlation function by Aw theta^(-0.8) to find Aw=29arcsec^(0.8). The estimated correlation length is r0=3.5(+/-0.3)h^(-1)Mpc in comoving units (Omega_m=0.3 and Omega_Lambda=0.7), which is slightly larger than the value for z~4 LBGs with i'<26. We calculate the angular correlation function for two sub-samples of the 87 LAEs divided by Lya luminosity, UV-continuum luminosity, and Lya equivalent width (EW). The Lya bright sub-sample shows a larger correlation amplitude than the Lya faint sub-sample, while no significant difference is found for the sub-samples divided by UV-continuum luminosity or EW.

We found a very large, 35 kpc, emission-line region around the Seyfert type 2 galaxy NGC 4388, using deep narrow-band imaging with the prime focus camera (Suprime-Cam) of the Subaru telescope. This region consists of many faint gas clouds or filaments, and extends northeastwards from the galaxy. The typical H-alpha luminosity L(H-alpha) of the clouds is 10^37 erg s^-1, and the total L(H-alpha) of the region within 10 kpc from the nucleus is 2 x 10^38 erg s^-1, which corresponds to an ionized gas mass of 10^5 M_solar. The map of the emission-line intensity ratio I([OIII])/I(H-alpha) indicates that the inner (r < 12 kpc) region of the VEELR may be excited by nuclear ionizing radiation. The excitation mechanism of the outer (r > 12 kpc) region is unclear, but it is likely that the nuclear radiation is also a dominant source of its ionization. We discuss the origin of the ionized gas. Two plausible origins of the gas in the VEELR are (i) the tidal debris resulting from a past interaction with a gas-rich dwarf galaxy, i.e., a minor merger, or (ii) the interstellar medium of NGC 4388, stripped by the ram pressure of the hot intracluster medium of the Virgo cluster.

We investigate the radial dependence of the spectroscopic properties, in particular the Mg2, <Fe> and H beta spectroscopic indices, in a sample of galaxies spanning a wide range of absolute luminosity in the Coma cluster. After allowing for the magnitude dependence of these indices, we find a significant gradient in Mg2, in the sense that galaxies in the core of the cluster have stronger Mg2. We find only weak gradients in <Fe> and H beta. Using the model grids presented in an earlier paper in this series, we attribute the Mg2 gradient to changes in metal abundance. One possible mechanism to create this abundance gradient is pressure confinement by the intracluster medium of material from Supernova driven winds early in the history of the galaxies.

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

(abr.) We present the analysis of the spectroscopic catalog of galaxies in the Coma cluster from Mobasher et al. This catalog comprises ~300 spectra of cluster members with magnitudes M_B=-20.5 to -14 in two areas of ~1 X 1.5 Mpc towards the center and the SW region of the cluster. In this paper the analysis is restricted to the 257 galaxies with no emission lines in their spectra. The strength of the age-sensitive indices is found to correlate with galaxy magnitude over the whole magnitude range. Similarly, the metallicity-sensitive indices anticorrelate with magnitude. By comparing the observed indices with model grids based on the Padova isochrones, we derive luminosity-weighted ages and metallicities. We present the distributions of ages and metallicities for galaxies in various magnitude bins. The mean metallicity decreases with galaxy magnitude and, at a given luminosity, appears to be generally lower for galaxies in the SW region of Coma as compared to the center of the cluster. A broad range of ages, from younger than 3 Gyr to older than 9 Gyr, is found in galaxies of any magnitude. However, systematic trends of age with luminosity are present among galaxies in the central field. In the central Mpc of Coma, a large fraction of galaxies at any luminosity (50-60% of the giants, >30 % of the dwarfs) show no evidence in their central regions of star formation occurred at redshift z<2, while the proportion of galaxies with significant star formation occurring at intermediate (0.35<z<2) and low (z<0.35) redshifts is found to depend on galaxy luminosity. Faint galaxies with young luminosity-weighted ages appear to have a bimodal metallicity distribution that, if confirmed, would point to a composite formation scenario. Finally, we present an interpretation of the index-magnitude relations observed.

We study the luminosity function and the correlation function of about 1200 z~4 Lyman break galaxies (LBGs) with i'<26 that are photometrically selected from deep BRi' imaging data of a 618 arcmin^2 area in the Subaru/XMM-Newton Deep Field taken with Subaru Prime Focus Camera. The contamination and completeness of our LBG sample are evaluated, on the basis of the Hubble Deep Field-North (HDF-N) objects, to be 17% and 45%, respectively. We derive the UV (rest 1700A) luminosity functions (LFs) and find a large population of UV-luminous galaxies at z~4. The LFs of the red and blue subsamples imply that the bright LBGs are redder in the UV continuum than the average color of the LBGs. Then we calculate the correlation function over theta = 2''-1000'' and find that it is fitted fairly well by a power law, omega(theta)=A_omega theta^(-0.8), with A_omega=0.71 +/- 0.26. We estimate the correlation length r_0 (in comoving units) of the two-point spatial correlation function xi(r) = (r/r_0)^(-1.8) to be r_0=2.7 +0.5/-0.6 h^(-1) Mpc (Omega_m=0.3 and Omega_Lambda=0.7). The correlation function shows an excess of omega (theta) on small scales (theta < 5''), departing from the power-law fit at > 3 sigma significance level. Interpreting this as being due to galaxy mergers, we evaluate the fraction of galaxies undergoing mergers to be 3.0 +/- 0.9%, which is significantly smaller than those of galaxies at intermediate redshifts.

The ages of stellar populations in 52 elliptical and S0 galaxies in the Coma cluster are investigated using a new spectroscopic sample. More than 40% of the S0s are found to have experienced recent star formation in their central regions during the last ~5 Gyrs, while such activity is absent in the ellipticals. Galaxies in this sample have absolute magnitudes in the range -20.5 < M_B < -17.5, and the fraction of S0 galaxies with recent star formation is higher at fainter luminosities. The observed luminosity range of S0 galaxies with signs of recent star formation activity is consistent with them being the descendants of typical star-forming spirals at intermediate redshift whose star formation has been halted as a consequence of the dense environment.

We study the clustering properties of about 1200 z~4 Lyman Break Galaxy (LBG) candidates with i'<26 which are selected by color from deep BRi' imaging data of a 618 arcmin^2 area in the Subaru/XMM-Newton Deep Field taken with Subaru Prime Focus Camera. The contamination and completeness of our LBG sample are evaluated, on the basis of the Hubble Deep Field North (HDFN) objects, to be 17% and 45%, respectively. We derive the angular correlation function over theta = 2''-1000'', and find that it is fitted fairly well by a power law, omega(theta)=A_omega theta^{-0.8}, with A_omega = 0.71 +/- 0.26. We then calculate the correlation length r0 (in comoving units) of the two-point spatial correlation function xi(r) = (r/r0)^{-1.8} from A_omega using the redshift distribution of LBGs derived from the HDFN, and find r0=2.7 (+0.5/-0.6) h^{-1} Mpc in a Lambda-dominated universe (Omega_m=0.3 and Omega_Lambda=0.7). This is twice larger than the correlation length of the dark matter at z~4 predicted from an analytic model by Peacock & Dodds but about twice smaller than that of bright galaxies predicted by a semi-analytic model of Baugh et al. We find an excess of omega(theta) on small scales (theta < 5'') departing from the power law fit over 3 sigma significance levels. Interpreting this as due to galaxy mergers, we estimate the fraction of galaxies undergoing mergers in our LBG sample to be 3.0 +/- 0.9%, which is significantly smaller than those of galaxies at intermediate redshifts.

This is the second paper in a series studying the photometric and spectroscopic properties of galaxies of different luminosities in the Coma cluster. The sample selection, spectroscopic observations and completeness functions are presented here. To study the spectral properties of galaxies as a function of their local environment, two fields were selected for spectroscopic observations to cover both the core (Coma1) and outskirts (ie. south-west of the core and centered on NGC4839)- (Coma3) of the cluster. Medium resolution spectroscopy (6-9 \AA) was carried out for a total of 490 galaxies in both fields (302 in Coma1 and 188 in Coma3), using the WYFFOS multi-fiber spectrograph on the William Herschel Telescope. The galaxies cover a range of $12 < R < 20$, corresponding to -23 < M_R < -15 (H0=65 km/sec/Mpc). The redshifts are measured with an accuracy of 100 km/sec. The spectral line strengths and equivalent widths are also measured for the same galaxies and analysed in Poggianti et al (2001- paper III). A total of 189 (Coma1) and 90 (Coma3) galaxies are identified as members of the Coma cluster. An analysis of the colors show that only two members of the Coma cluster in our sample have B-R > 2. The completeness functions for the spectroscopic sample is presented.

We present bright galaxy number counts in five broad bands ($u', g', r', i', z'$) from imaging data taken during the commissioning phase of the Sloan Digital Sky Survey (SDSS). The counts are derived from two independent stripes of imaging scans along the Celestial Equator, one each toward the North and the South Galactic cap, covering about 230 and 210 square degrees, respectively. A careful study is made to verify the reliability of the photometric catalog. For galaxies brighter than $r^* = 16$, the catalog produced by automated software is examined against eye inspection of all objects. Statistically meaningful results on the galaxy counts are obtained in the magnitude range $12 \le r^* \le 21$, using a sample of 900,000 galaxies. The counts from the two stripes differ by about 30% at magnitudes brighter than $r^*= 15.5$, consistent with a local $2\sigma$ fluctuation due to large scale structure in the galaxy distribution. The shape of the number counts-magnitude relation brighter than $r^* = 16$ is well characterized by $N \propto 10^{0.6m}$, the relation expected for a homogeneous galaxy distribution in a ``Euclidean'' universe. In the magnitude range $ 16 < r^* < 21$, the galaxy counts from both stripes agree very well, and follow the prediction of the no-evolution model, although the data do not exclude a small amount of evolution. We use empirically determined color transformations to derive the galaxy number counts in the $B$ and $I_{814}$ bands. We compute the luminosity density of the universe at zero redshift in the five SDSS bands and in the $B$ band. We find ${\cal L}_{B} = 2.4 \pm 0.4 \times 10^8L_\odot h $Mpc$^{-3}$, for a reasonably wide range of parameters of the Schechter luminosity function in the $B$ band.

We summarise the properties of the Sloan Digital Sky Survey (SDSS) project, discuss our software infrastructure, and outline the architecture of the SDSS image processing pipelines. We then discuss two of the algorithms used in the SDSS image processing; the KL-transform based modelling of the spatial variation of the PSF, and the use of galaxy models in star/galaxy separation. We conclude with the first author's personal opinions on the challenges that the astronomical community faces with major software projects.

We investigate the photometric properties of 456 bright galaxies using imaging data recorded during the commissioning phase of the Sloan Digital Sky Survey (SDSS). Morphological classification is carried out by correlating results of several human classifiers. Our purpose is to examine the statistical properties of color indices, scale lengths, and concentration indices as functions of morphology for the SDSS photometric system. We find that $u'-g'$, $g'-r'$, and $r'-i'$ colors of SDSS galaxies match well with those expected from the synthetic calculation of spectroscopic energy distribution of template galaxies and with those transformed from $UBVR_CI_C$ color data of nearby galaxies. The agreement is somewhat poor, however, for $i'-z'$ color band with a discrepancy of $0.1-0.2$ mag. With the aid of the relation between surface brightness and radius obtained by Kent (1985), we estimate the averages of the effective radius of early type galaxies and the scale length of exponential disks both to be 2.6 kpc for $L^*$ galaxies. We find that the half light radius of galaxies depends slightly on the color bands, consistent with the expected distribution of star-forming regions for late type galaxies and with the known color gradient for early type galaxies. We also show that the (inverse) concentration index, defined by the ratio of the half light Petrosian radius to the 90% light Petrosian radius, correlates tightly with the morphological type; this index allows us to classify galaxies into early (E/S0) and late (spiral and irregular) types, allowing for a 15-20% contamination from the opposite class compared with eye-classified morphology.

Infrared J and H surface photometry are carried out for nearby 12 galaxies whose distances have been accurately measured via HST Cepheid observations. Using the total, isophotal, and surface-photometric aperture magnitudes we calibrate the infrared luminosity-line width relation (IRTF). It is found that IRTF changes its slope at log W_{20}^c=2.45 in all the examined magnitude systems. The apparent scatter of IRTF is not significantly reduced when surface photometric magnitudes are used instead of the conventionally used synthetic aperture magnitude H_{-0.5}. It is also shown that the color (I-H)_T of the nearby calibrator galaxies is redder by 0.2mag than the Coma cluster galaxies, but such a trend is not clearly visible for the Ursa Major mostly because of poor statistics. The color offset of the Coma is analogous to that previously found in I_T-H_{-0.5}. From the present calibration of H-band IRTF, we obtain the distance to the Coma cluster to be m-M=34.94 +/- 0.13mag, where no account is taken of the I-H color problem. Using the CMB-rest recession velocity of the Coma cluster we obtain H_0=73 +/- 4 km/s/Mpc.

We investigate pure luminosity evolution models for early-type (elliptical and S0) galaxies (i.e., no number density change or morphology transition), and examine whether these models are consistent with observed number counts in the B, I and K bands and redshift distributions of two samples of faint galaxies selected in the I and K bands. The models are characterized by the star formation time scale $\tau_{SF}$ and the time $t_{gw}$ when galactic wind blows in addition to several other conventional parameters. We find the single-burst model ($\tau_{SF}$=0.1 Gyr and $t_{gw}$=0.353 Gyr), which is known to reproduce the photometric properties of early-type galaxies in clusters, is inconsistent with redshift distributions of early-type galaxies in the field environment due to overpredictions of galaxies at $z\gsim1.4$ even with strong extinction which is at work until $t_{gw}$. In order for dust extinction to be more effective, we change $\tau_{SF}$ and $t_{gw}$ as free parameters, and find that models with $\tau_{SF}\gsim0.5$ Gyr and $t_{gw}>1.0$ Gyr can be made consistent with both the observed redshift distributions and number counts, if we introduce strong extinction $(E(B-V)\geq1$ as a peak value). These results suggest that early-type galaxies in the field environment do not have the same evolutionary history as described by the single-burst model.

We show that the growth curves of galaxies can be used to determine their bulge and disk parameters and bulge-to-total luminosity ratios, in addition to their conventional asymptotic magnitudes, provided that the point spread function is accurately known and signal-to-noise ratio is modest (S/N$\gtrsim30$). The growth curve is a fundamental quantity that most future large galaxy imaging surveys will measure. Bulge and disk parameters retrieved from the growth curve will enable us to perform statistical studies of luminosity structure for a large number of galaxies.

We carry out CCD photometry of galaxies in the 5.25 square region centered on Coma cluster down to $M_R=-16.0$, beyond the limit of conventional morphological classification. We use the angular two-point correlation function as well as radial profiles in order to characterize the luminosity segregation. We find strong luminosity segregation for our total sample over the magnitude range of $-20 \leq M_R \leq -16$, which is not entirely accounted for in terms of the morphology-density relation that is known to exist only for bright galaxies. We use a single consistent parameter, the degree of luminosity concentration, to parameterize the morphology of galaxies over the wide magnitude range, where both giant and dwarf galaxies are included. Galaxies with high central concentration (HCC) show strong luminosity segregation, i.e. their clustering strength depends strongly on luminosity while those with low central concentration (LCC) show almost no luminosity segregation. Radial density profile shows that brighter HCC-type galaxies tend to more strongly concentrate near the cluster center while LCC-type galaxies do not show such a dependence on luminosity. We show that these results are tenable against the contamination by field galaxies and uncertainties in our method of classification.