We investigate Box/Peanut and bar structures in image data of edge-on and face-on nearby galaxies taken from SDSS (Sloan Digital Sky Survey) to present catalogues containing the surface brightness parameters and the morphology classification. About 1700 edge-on galaxies and 2600 face-on galaxies are selected from SDSS DR7 in g, r and i-band. The images of each galaxy are fitted with the model of 2-dimensional surface brightness of S\'{e}rsic bulge and exponential disk. After removing some irregular data, Box/Peanut, bar and other structures are easily distinguished by eye using residual (observed minus model) images. We find 292 Box/Peanut structures in the 1329 edge-on samples and 630 bar structures in 1890 face-on samples in i-band, after removing some irregular data. Then the fraction of Box/Peanut galaxies is about 22 percent against the edge-on samples, and that of bar is about 33 percent (about 50 percent if 629 elliptical galaxies are removed) against the face-on samples. Furthermore the strengths of the Box/Peanut and bar are evaluated as "Strong", "Standard" and "Weak". We find that the strength increases slightly with increasing $B/T$ (Bulge-to-Total flux ratio), and that the fraction of Box/Peanut is generally about a half of that of bar irrespective of the strength and the $B/T$. Our result supports the idea that the Box/Peanut is bar seen edge-on.

In this work, we investigate the dependence of the covering factor (CF) of active galactic nuclei (AGNs) (i) on the mid-infrared (MIR) luminosity and (ii) on the redshift. We constructed 12- and 22-micron luminosity functions (LFs) at 0.006 < z < 0.3 using the Wide-field Infrared Survey Explorer} (WISE) data. Combining the WISE catalog with the Sloan Digital Sky Survey (SDSS) spectroscopic data, we selected 223,982 galaxies at 12 micron and 25,721 galaxies at 22 micron for spectroscopic classification. We then identified 16,355 AGNs at 12 micron and 4,683 AGNs at 22 micron by their optical emission lines and cataloged classifications in the SDSS. Following that, we estimated the CF as the fraction of type 2 AGN in all AGNs whose MIR emissions are dominated by the active nucleus (not their host galaxies) based on their MIR colors. We found that (i) the CF decreased with increasing MIR luminosity, regardless of the choice of type 2 AGN classification criteria, and (ii) the CF did not change significantly with the redshift for z < 0.2. Furthermore, we carried out various tests to determine the influence of selection bias and confirmed similar dependences exist even when taking these uncertainties into account. The luminosity dependence of the CF can be explained by the receding torus model, but the "modified" receding torus model gives a slightly better fit, as suggested by Simpson.

We present the 9 and 18 micron luminosity functions (LFs) of galaxies at 0.006 < z < 0.8 (with an average redshift of ~ 0.04) using the AKARI mid-infrared all-sky survey catalog. We selected 243 galaxies at 9 micron and 255 galaxies at 18 micron from the Sloan Digital Sky Survey (SDSS) spectroscopy region. These galaxies were then classified by their optical emission lines, such as the line width of H_alpha or by their emission line ratios of [OIII]/H_beta and [NII]/H_alpha into five types: Type 1 active galactic nuclei (AGN) (Type 1); Type 2 AGN (Type 2); low-ionization narrow emission line galaxies (LINER); galaxies with both star formation and narrow-line AGN activity (composite galaxies); and star-forming galaxies (SF). We found that (i) the number density ratio of Type 2 to Type 1 AGNs is 1.73 +/- 0.36, which is larger than a result obtained from the optical LF and (ii) this ratio decreases with increasing 18 micron luminosity.

We have taken a deep, moderate-resolution Keck/Deimos spectra of QSO, CFHQS2329, at z=6.4. At the wavelength of Lya, the spectrum shows a spatially-extended component, which is significantly more extended than a stellar spectrum, and also a continuum part of the spectrum. The restframe line width of the extended component is 21+-7 A, and thus smaller than that of QSO (52+-4 A), where they should be identical if the light is incomplete subtraction of the QSO component. Therefore, these comparisons argue for the detection of a spatially extended Lya nebulae around this QSO. This is the first z>6 QSO that an extended Lya halo has been observed around. Careful subtraction of the central QSO spectrum reveals a lower limit to the Lya luminosity of (1.7+-0.1)x 10^43 erg s^-1. This emission may be from the theoretically predicted infalling gas in the process of forming a primordial galaxy that is ionized by a central QSO. On the other hand, if it is photoionized by the host galaxy, an estimated star-formation rate of >3.0 Msun yr^-1 is required. If we assume the gas is virialized, we obtain dynamical mass estimate of Mdyn=1.2x10^12 Msun. The derived MBH/Mhost is 2.1x10^-4, which is two orders smaller than those from more massive z~6 QSOs, and places this galaxy in accordance with the local M-sigma relation, in contrast to a previous claim on the evolution of M-sigma relation at z~6. We do not claim evolution or non-evolution of the M-sigma relation based on a single object, but our result highlights the importance of investigating fainter QSOs at z~6.

We develop a software kit called "2MASS Catalog Server Kit" to easily construct a high-performance database server for the 2MASS Point Source Catalog (includes 470,992,970 objects) and several all-sky catalogs. Users can perform fast radial search and rectangular search using provided stored functions in SQL similar to SDSS SkyServer. Our software kit utilizes open-source RDBMS, and therefore any astronomers and developers can install our kit on their personal computers for research, observation, etc. Out kit is tuned for optimal coordinate search performance. We implement an effective radial search using an orthogonal coordinate system, which does not need any techniques that depend on HTM or HEALpix. Applying the xyz coordinate system to the database index, we can easily implement a system of fast radial search for relatively small (less than several million rows) catalogs. To enable high-speed search of huge catalogs on RDBMS, we apply three additional techniques: table partitioning, composite expression index, and optimization in stored functions. As a result, we obtain satisfactory performance of radial search for the 2MASS catalog. Our system can also perform fast rectangular search. It is implemented using techniques similar to those applied for radial search. Our way of implementation enables a compact system and will give important hints for a low-cost development of other huge catalog databases.

The AKARI All-Sky Catalogues are an important infrared astronomical database for next-generation astronomy that take over the IRAS catalog. We have developed an online service, AKARI Catalogue Archive Server (AKARI-CAS), for astronomers. The service includes useful and attractive search tools and visual tools. One of the new features of AKARI-CAS is cached SIMBAD/NED entries, which can match AKARI catalogs with other catalogs stored in SIMBAD or NED. To allow advanced queries to the databases, direct input of SQL is also supported. In those queries, fast dynamic cross-identification between registered catalogs is a remarkable feature. In addition, multiwavelength quick-look images are displayed in the visualization tools, which will increase the value of the service. In the construction of our service, we considered a wide variety of astronomers' requirements. As a result of our discussion, we concluded that supporting users' SQL submissions is the best solution for the requirements. Therefore, we implemented an RDBMS layer so that it covered important facilities including the whole processing of tables. We found that PostgreSQL is the best open-source RDBMS products for such purpose, and we wrote codes for both simple and advanced searches into the SQL stored functions. To implement such stored functions for fast radial search and cross-identification with minimum cost, we applied a simple technique that is not based on dividing celestial sphere such as HTM or HEALPix. In contrast, the Web application layer became compact, and was written in simple procedural PHP codes. In total, our system realizes cost-effective maintenance and enhancements.

Understanding the cosmic re-ionization is one of the key goals of the modern observational cosmology. High redshift QSO spectra can be used as background light sources to measure absorption by intervening neutral hydrogen. We investigate neutral hydrogen absorption in a deep, moderate-resolution Keck/Deimos spectrum of QSO CFHQSJ2329-0301 at z=6.4. This QSO is one of the highest redshift QSOs presently known at z=6.4 but is 2.5 mag fainter than a previously well-studied QSO SDSSJ1148+5251 at z=6.4. Therefore, it has a smaller Stromgren sphere, and allows us to probe the highest redshift hydrogen absorption to date. The average transmitted flux at 5.915<z_abs<6.365 (200 comoving Mpc) is consistent with zero, in Ly_alpha, Ly_beta, and Ly_gamma absorption measurements. This corresponds to the lower limit of optical depth, tau_eff>4.9. These results are consistent with strong evolution of the optical depth at z>5.7.

By cross-correlating AKARI infrared (IR) sources with the SDSS galaxies, we identified 2357 infrared galaxies with a spectroscopic redshift. This is not just one of the largest samples of local IR galaxies, but AKARI provides crucial FIR bands (9, 18, 65, 90, 140, and 160um) in accurately measuring galaxy SED across the peak of the dust emission at ~100um. By fitting modern IR SED models to the AKARI photometry, we measured the total infrared luminosity (L_IR) of individual galaxies more accurately. Using this L_IR, we constructed luminosity functions of infrared galaxies at a median redshift of z=0.031, with 4 times larger sample than previous work. The LF agrees well with that at z=0.0082 (RBGS), showing smooth and continuous evolution toward higher redshift LFs measured in the AKARI NEP deep field. The derived local cosmic IR luminosity density is Omega_IR=3.8x10^8 LsunMpc^-3. We separate galaxies into AGN, star-forming, and composite by using the [NII]/Ha vs [OIII]/Hb line ratios. The fraction of AGN shows a continuous increase with increasing L_IR from 25% to 90% at 9<log L_IR<12.5. The SFR_Ha and L_[OIII] show good correlations with L_IR for SFG (star-forming galaxies) and AGN, respectively. The self-absorption corrected Ha/Hb ratio shows a weak increase with L_IR with a substantial scatter. When we separate IR LFs into contributions from AGN and SFG, the AGN contribution becomes dominant at L_IR>10^11Lsun, coinciding the break of the both SFG and AGN IR LFs. At L_IR<10^11Lsun, SFG dominates IR Lfs. Only 1.1% of Omega_IR is produced by LIRG, and only 0.03% is by ULIRG in the local Universe. This work also provides the most accurate infrared luminosity density of the local Universe to date. Compared with high redshift results from the AKARI NEP deep survey, we observed a strong evolution of Omega_IR^SFG ~(1+z)^4.1+-0.4 and Omega_IR^AGN ~(1+z)^4.1+-0.5 (abridged).

Infrared (IR) luminosity is fundamental to understanding the cosmic star formation history and AGN evolution. The AKARI IR space telescope performed all sky survey in 6 IR bands (9, 18, 65, 90, 140, and 160um) with 3-10 times better sensitivity than IRAS, covering the crucial far-IR wavelengths across the peak of the dust emission. Combined with a better spatial resolution, AKARI can much more precisely measure the total infrared luminosity (L_TIR) of individual galaxies, and thus, the total infrared luminosity density in the local Universe. By fitting IR SED models, we have re-measured L_TIR of the IRAS Revised Bright Galaxy Sample. We present mid-IR monochromatic luminosity to L_TIR conversions for Spitzer 8,24um, AKARI 9,18um, IRAS 12um, WISE 12,22um, and ISO 15um filters, with scatter ranging 13-44%. The resulting AKARI IR luminosity function (LF) agrees well with that from the IRAS. We integrate the LF weighted by L_TIR to obtain a cosmic IR luminosity density of Omega_TIR= (8.5^{+1.5}_{-2.3})x 10^7 L Mpc^-3, of which 7+-1% is produced by LIRGs, and only 0.4+-0.1% is from ULIRGs in the local Universe. Once IR contributions from AGN and star-forming galaxies (SFG) are separated, SFG IR LF shows a steep decline at the bright-end. Compared with high-redshift results from the AKARI NEP deep survey, these data show a strong evolution of Omega_TIRSF propto (1+z)^4.0+-0.5, and Omega_TIRAGN propto (1+z)^4.4+-0.4. For Omega_TIRAGN, the ULIRG contribution exceeds that from LIRG already by z~1. A rapid evolution in both Omega_TIRAGN and Omega_TIRSFG suggests the correlation between star formation and black hole accretion rate continues up to higher redshifts. We compare the evolution of Omega_TIRAGN to that of X-ray luminosity density. The Omega_TIRAGN/Omega_X-rayAGN ratio shows a possible increase at z>1, suggesting an increase of obscured AGN at z>1.

Context : AKARI is the first Japanese astronomical satellite dedicated to infrar ed astronomy. One of the main purposes of AKARI is the all-sky survey performed with six infrared bands between 9 and 200um during the period from 2006 May 6 to 2007 August 28. In this paper, we present the mid-infrared part (9um and 18um b ands) of the survey carried out with one of the on-board instruments, the Infrar ed Camera (IRC). Aims : We present unprecedented observational results of the 9 and 18um AKARI al l-sky survey and detail the operation and data processing leading to the point s ource detection and measurements. Methods : The raw data are processed to produce small images for every scan and point sources candidates, above the 5-sigma noise level per single scan, are der ived. The celestial coordinates and fluxes of the events are determined statisti cally and the reliability of their detections is secured through multiple detect ions of the same source within milli-seconds, hours, and months from each other. Results : The sky coverage is more than 90% for both bands. A total of 877,091 s ources (851,189 for 9um, 195,893 for 18um) are confirmed and included in the cur rent release of the point source catalogue. The detection limit for point source s is 50mJy and 90mJy for the 9um and 18um bands, respectively. The position accu racy is estimated to be better than 2". Uncertainties in the in-flight absolute flux calibration are estimated to be 3% for the 9um band and 4% for the 18um ban d. The coordinates and fluxes of detected sources in this survey are also compar ed with those of the IRAS survey and found to be statistically consistent.

We reveal cosmic star-formation history obscured by dust using deep infrared observation with the AKARI. A continuous filter coverage in the mid-IR wavelength (2.4, 3.2, 4.1, 7, 9, 11, 15, 18, and 24um) by the AKARI satellite allows us to estimate restframe 8um and 12um luminosities without using a large extrapolation based on a SED fit, which was the largest uncertainty in previous work. We found that restframe 8um (0.38<z<2.2), 12um (0.15<z<1.16), and total infrared (TIR) luminosity functions (LFs) (0.2<z<1.6) constructed from the AKARI NEP deep data, show a continuous and strong evolution toward higher redshift. In terms of cosmic infrared luminosity density (Omega_IR), which was obtained by integrating analytic fits to the LFs, we found a good agreement with previous work at z<1.2, with Omega_IR propto (1+z)^4.4+-1.0. When we separate contributions to Omega_IR by LIRGs and ULIRGs, we found more IR luminous sources are increasingly more important at higher redshift. We found that the ULIRG (LIRG) contribution increases by a factor of 10 (1.8) from z=0.35 to z=1.4.

We demonstrate the capability of AKARI for mapping diffuse far-infrared emission and achieved reliability of all-sky diffuse map. We have conducted an all-sky survey for more than 94 % of the whole sky during cold phase of AKARI observation in 2006 Feb. -- 2007 Aug. The survey in far-infrared waveband covers 50 um -- 180 um with four bands centered at 65 um, 90 um, 140 um, and 160 um and spatial resolution of 3000 -- 4000 (FWHM).This survey has allowed us to make a revolutionary improvement compared to the IRAS survey that was conducted in 1983 in both spatial resolution and sensitivity after more than a quarter of a century. Additionally, it will provide us the first all-sky survey data with high-spatial resolution beyond 100 um. Considering its extreme importance of the AKARI far-infrared diffuse emission map, we are now investigating carefully the quality of the data for possible release of the archival data. Critical subjects in making image of diffuse emission from detected signal are the transient response and long-term stability of the far-infrared detectors. Quantitative evaluation of these characteristics is the key to achieve sensitivity comparable to or better than that for point sources (< 20 -- 95 [MJy/sr]). We describe current activities and progress that are focused on making high quality all-sky survey images of the diffuse far-infrared emission.

We discuss the capability of AKARI in recovering diffuse far-infrared emission, and examine the achieved reliability. Critical issues in making images of diffuse emission are the transient response and long-term stability of the far-infrared detectors. Quantitative evaluation of these characteristics are the key to achieving sensitivity comparable to or better than that for point sources (< 20 -- 95 MJy sr-1). We describe current activity and progress toward the production of high quality images of the diffuse far-infrared emission using the AKARI all-sky survey data.

Based on our intensive spectroscopic campaign with the GoldCam spectrograph on the Kitt Peak National Observatory (KPNO) 2.1-m telescope, we have constructed the first catalogue of E+A galaxies with spectroscopic companion galaxies, and investigated a probability that an E+A galaxy have close companion galaxies. We selected 660 E+A galaxies with $4.0{\rm \AA} < {\rm H}\delta {\rm EW}$ at a redshift of $<0.167$ from the Data Release 5 of the Sloan Digital Sky Survey (SDSS). We selected their companion candidates from the SDSS imaging data, and classified them into true companions, fore/background galaxies and companion candidates using the SDSS and our KPNO spectra. We observed 26 companion candidates of E+A galaxies at the KPNO to measure their redshifts. Their spectra showed that 17 targets are true companion galaxies. The number of spectroscopically-confirmed E+A's companions are now 34. This becomes the first catalogue of E+A galaxies with spectroscopic companion systems. We found that E+A galaxies have an 54% larger probability of having companion galaxies (7.88%) compared to the comparison sample of normal galaxies (5.12%). A statistical test shows the probabilities are different with 99.7% significance. Our results based on spectroscopy tightens the connection between the dynamical merger/interaction and the origin of E+A galaxies.

We have performed a spatially-resolved medium resolution long-slit spectroscopy of a nearby E+A (post-starburst) galaxy system, SDSSJ161330.18+510335.5. This E+A galaxy has an obvious companion galaxy 14kpc in front with the velocity difference of 61.8 km/s. Both galaxies have obviously disturbed morphology We have found that H$\delta$ equivalent width (EW) of the E+A galaxy is greater than 7\AA galaxy wide (8.5 kpc). The E+A galaxy have a weak [OIII] emission (EW$\sim$1\AA) by $\sim$2.6 kpc offset from the peak of the Balmer absorption lines. We detected a rotational velocity in the companion galaxy of $>$175km/s. The progenitor of the companion may have been a rotationally-supported, but yet passive S0 galaxy. The age of the E+A galaxy after quenching the star formation is estimated to be 100-500 Myr, with its centre having slightly younger stellar population. These findings are inconsistent with a simple picture where the dynamical interaction creates infall of the gas reservoir that causes the central starburst/post-starburst. Instead, our results present an important example where the galaxy-galaxy interaction can trigger a galaxy-wide post-starburst phenomena.

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

We present the spatially resolved observations of IRAS sources from the Japanese infrared astronomy satellite AKARI All-Sky Survey during the performance verification (PV) phase of the mission. We extracted reliable point sources matched with IRAS point source catalogue. By comparing IRAS and AKARI fluxes, we found that the flux measurements of some IRAS sources could have been over or underestimated and affected by the local background rather than the global background. We also found possible candidates for new AKARI sources and confirmed that AKARI observations resolved IRAS sources into multiple sources. All-Sky Survey observations are expected to verify the accuracies of IRAS flux measurements and to find new extragalactic point sources.

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.

We investigated the radial colour gradient and two-dimensional (2D) colour properties of 22 E+A galaxies with 5.5 < Hdelta equivalent width(EW) < 8.5\AA and 49 normal early-type galaxies as a control sample at a redshift of <0.2 in the Second Data Release of the Sloan Digital Sky Survey. We found that a substantial number of E+A galaxies exhibit positive slopes of radial colour gradient (bluer gradients toward the centre) which are seldom seen in normal early-type galaxies. We found irregular 'colour morphologies'-- asymmetrical and clumpy patterns -- at the centre of g-r and r-i 2D colour maps of E+A galaxies with positive slopes of colour gradient. Kolomogorov-Smirnov two-sample tests show that g-r and r-i colour gradient distributions of E+A galaxies differ from those of early-type galaxies with a more than 99.99 per cent significance level. We also found a tight correlation between radial colour gradients and colours, and between radial colour gradients and the 4000-\AA break in the E+A sample; E+A galaxies which exhibit bluer colour or weaker D4000 tend to have positive slopes of radial colour gradient. We compared the GISSEL model and E+A observational quantities, Hdelta EW, D4000 and u-g colour, and found that almost all our E+A galaxies are situated along a single evolutionary track. Therefore, these results are interpreted as E+A galaxies evolving from Hdelta EW ~ 8 \AA to Hdelta EW ~ 5 \AA, with colour gradients changing from positive to negative, and with the irregular 2D colour map becoming smoother, during a time-scale of ~300 Myr. Our results favour the hypothesis that E+A galaxies are post-starburst galaxies caused by merger/interaction, having undergone a centralized violent starburst.

Passive spiral galaxies - unusual galaxies with spiral morphologies but without any sign of on-going star formation - have recently been shown to exist preferentially in cluster infalling regions. This discovery directly connects passive spiral galaxies to cluster galaxy evolution studies, such as the Butcher-Oemler effect or the morphology-density relation. Thus, detailed study of passive spiral galaxies could potentially yield new insight on the underlying physical mechanisms governing cluster galaxy evolution. However, in previous work, passive spiral galaxies were selected from low-resolution optical images with $\sim$1.5 arcsec seeing. Passive spirals could therefore be mis-identified as S0 galaxies, or as dusty-starburst galaxies which are not passive at all. To address this issue we performed deep, high-resolution, near-infrared imaging of 32 passive spiral galaxies with UKIRT.

Passive spiral galaxies -- unusual galaxies with spiral morphology without any sign of on-going star formation -- have recently been discovered to exist preferentially in cluster infalling regions (at about the virial radius, or at a local galaxy density of $\sim 1$ Mpc$^{-2}$). The discovery directly connects the passive spiral galaxies to the cluster galaxy evolution studies such as the Butcher-Oemler effect or the morphology-density relation, i.e., passive spiral galaxies are likely to be transition objects between high-z blue, spiral galaxies and low-z red, cluster early-type galaxies. Thus, detailed study of passive spiral galaxies potentially could bring a new insight on the underlying physical mechanisms governing cluster galaxy evolution. However, in previous work, passive spiral galaxies are selected from the low resolution optical images with $\sim 1.5$ arcsec of seeing. Therefore, passive spirals could be a mis-identification of S0 galaxies; or dusty-starburst galaxies which are not passive at all. To answer these questions, we performed a deep, high-resolution, near-infrared imaging of 32 passive spiral galaxies. Our high resolution $K$ band images show clear spiral arm structures. Thus, passive spirals are not S0s. Optical-infrared colour does not show any signs of dusty-starburst at all. Therefore, it is likely that they are truly ``passive'' and ``spiral'' galaxies in the midst of cluster galaxy evolution.

We have studied the morphology-density relation and morphology--cluster-centric-radius relation using a volume limited sample (0.05$<z<$0.1, $Mr^*<-$20.5) of the Sloan Digital Sky Survey (SDSS) data. We found there are two characteristic changes in both the morphology-density and the morphology-radius relations, suggesting two different mechanisms are responsible for the relations. In the sparsest regions (below 1 Mpc$^{-2}$ or outside of 1 virial radius), both relations become less noticeable, suggesting the responsible physical mechanisms for galaxy morphological change require denser environment. In the intermediate density regions (density between 1 and 6 Mpc$^{-2}$ or virial radius between 0.3 and 1), intermediate-type fractions increase toward denser regions, whereas late-disc fractions decrease. In the densest regions (above 6 Mpc$^{-2}$ or inside of 0.3 virial radius), intermediate-type fractions decrease radically and early-type fractions increase in turn. This is a contrasting result to that in intermediate regions and it suggests that yet another mechanism is more responsible for the morphological change in these regions. We also compared the morphology-density relation from the SDSS (0.01$<z<$0.054) with that of the MORPHS data ($z\sim$0.5). Two relations lie on top of each other, suggesting that the morphology-density relation was already established at $z\sim$0.5 as is in the present universe.

In previous work on galaxy clusters, several authors reported a discovery of an unusual population of galaxies, which have spiral morphologies, but do not show any star formation activity. These galaxies are called ``passive spirals'', and have been interesting since it has been difficult to understand the existence of such galaxies. Using a volume limited sample (0.05<z<0.1 and Mr<-20.5; 25813 galaxies) of the Sloan Digital Sky Survey data, we have found 73 (0.28$\pm$0.03%) passive spiral galaxies and studied their environments. It is found that passive spiral galaxies live in local galaxy density 1-2 Mpc$^{-2}$ and 1-10 cluster-centric virial radius. Thus the origins of passive spiral galaxies are likely to be cluster related. These characteristic environments coincide with the previously reported environment where galaxy star formation rate suddenly declines and the so-called morphology-density relation turns. It is likely that the same physical mechanism is responsible for all of these observational results. The existence of passive spiral galaxies suggests that a physical mechanism that works calmly is preferred to dynamical origins such as major merger/interaction since such a mechanism can destroy spiral arm structures. Compared with observed cluster galaxy evolution such as the Butcher-Oemler effect and the morphological Butcher-Oemler effect, passive spiral galaxies are likely to be a galaxy population in transition between red, elliptical/S0 galaxies in low redshift clusters and blue, spiral galaxies more numerous in higher redshift clusters.