We explore the induction of circular polarization (CP) of near-infrared light in star-forming regions using three-dimensional radiative transfer calculations. The simulations trace the change of Stokes parameters at each scattering/absorption process in a dusty gas slab composed of aligned grains. We find that the CP degree enlarges significantly according as the size of dust grains increases and exceeds $\sim 20$ percent for micron-size grains. Therefore, if micron-size grains are dominant in a dusty gas slab, the high CP observed around luminous young stellar objects can be accounted for. The distributions of CP show the asymmetric quadrupole patters regardless of the grain sizes. Also, we find that the CP degree depends on the relative position of a dusty gas slab. If a dusty gas slab is located behind a star-forming region, the CP reaches $\sim 60$ percent in the case of $1.0~{\rm \mu m}$ dust grains. Hence, we suggest that the observed variety of CP maps can be explained by different size distributions of dust grains and the configuration of aligned grains.

We have initiated a new survey for local extremely metal-poor galaxies (EMPGs) with Subaru/Hyper Suprime-Cam (HSC) large-area (~500 deg^2) optical images reaching a 5 sigma limit of ~26 magnitude, about 100 times deeper than the Sloan Digital Sky Survey (SDSS). To select Z/Z_sun<0.1 EMPGs from ~40 million sources detected in the Subaru images, we first develop a machine-learning (ML) classifier based on a deep neural network algorithm with a training data set consisting of optical photometry of galaxy, star, and QSO models. We test our ML classifier with SDSS objects having spectroscopic metallicity measurements, and confirm that our ML classifier accomplishes 86%-completeness and 46%-purity EMPG classifications with photometric data. Applying our ML classifier to the photometric data of the Subaru sources as well as faint SDSS objects with no spectroscopic data, we obtain 27 and 86 EMPG candidates from the Subaru and SDSS photometric data, respectively. We conduct optical follow-up spectroscopy for 10 out of our EMPG candidates with Magellan/LDSS-3+MagE, Keck/DEIMOS, and Subaru/FOCAS, and find that the 10 EMPG candidates are star-forming galaxies at z=0.007-0.03 with large H_beta equivalent widths of 104-265 A, stellar masses of log(M*/M_sun)=5.0-7.1, and high specific star-formation rates of ~300 Gyr^{-1}, which are similar to those of early galaxies at z>6 reported recently. We spectroscopically confirm that 3 out of 10 candidates are truly EMPGs with Z/Z_sun<0.1, one of which is HSC J1631+4426, the most metal-poor galaxy with Z/Z_sun=0.016 reported ever.

We present element abundance ratios and ionizing radiation of local young low-mass (~$10^{6}$ M_sun) extremely metal poor galaxies (EMPGs) with a 2% solar oxygen abundance (O/H)_sun and a high specific star-formation rate (sSFR~300 Gyr$^{-1}$), and other (extremely) metal poor galaxies, which are compiled from Extremely Metal-Poor Representatives Explored by the Subaru Survey (EMPRESS) and the literature. Weak emission lines such as [FeIII]4658 and HeII4686 are detected in very deep optical spectra of the EMPGs taken with 8m-class telescopes including Keck and Subaru (Kojima et al. 2019, Izotov et al. 2018), enabling us to derive element abundance ratios with photoionization models. We find that neon- and argon-to-oxygen ratios are comparable to those of known local dwarf galaxies, and that the nitrogen-to-oxygen abundance ratios (N/O) are lower than 20% (N/O)_sun consistent with the low oxygen abundance. However, the iron-to-oxygen abundance ratios (Fe/O) of the EMPGs are generally high; the EMPGs with the 2%-solar oxygen abundance show high Fe/O ratios of ~90-140% (Fe/O)_sun, which are unlikely explained by suggested scenarios of Type Ia supernova iron productions, iron's dust depletion, and metal-poor gas inflow onto previously metal-riched galaxies with solar abundances. Moreover, these EMPGs have very high HeII4686/H$\beta$ ratios of ~1/40, which are not reproduced by existing models of high-mass X-ray binaries whose progenitor stellar masses are less than 120 M_sun. Comparing stellar-nucleosynthesis and photoionization models with a comprehensive sample of EMPGs identified by this and previous EMPG studies, we propose that both the high Fe/O ratios and the high HeII4686/H$\beta$ ratios are explained by the past existence of super massive ($>$300 M_sun) stars, which may evolve into intermediate-mass black holes ($\gtrsim$100 M_sun).

We investigate the distribution of metals in the cosmological volume at $z\sim3$, in particular, provided by massive population III (Pop III) stars using a cosmological $N$-body simulation in which a model of Pop III star formation is implemented. Owing to the simulation, we can choose minihaloes where Pop III star formation occurs at $z>10$ and obtain the spatial distribution of the metals at lower-redshifts. To evaluate the amount of heavy elements provided by Pop III stars, we consider metal yield of pair-instability or core-collapse supernovae (SNe) explosions of massive stars. By comparing our results to the Illustris-1 simulation, we find that heavy elements provided by Pop III stars often dominate those from galaxies in low density regions. The median value of the volume averaged metallicity is $Z\sim 10^{-4.5 - -2} Z_{\odot}$ at the regions. Spectroscopic observations with the next generation telescopes are expected to detect the metals imprinted on quasar spectra.

The survey of Lyman $\alpha$ emitters (LAEs) with Subaru Hyper Suprime-Cam, called SILVERRUSH (Ouchi et al.), is producing massive data of LAEs at $z\gtrsim6$. Here we present LAE simulations to compare the SILVERRUSH data. In 162$^3$ comoving Mpc$^3$ boxes, where numerical radiative transfer calculations of reionization were performed, LAEs have been modeled with physically motivated analytic recipes as a function of halo mass. We have examined $2^3$ models depending on the presence or absence of dispersion of halo Ly$\alpha$ emissivity, dispersion of the halo Ly$\alpha$ optical depth, $\tau_\alpha$, and halo mass dependence of $\tau_\alpha$. The unique free parameter in our model, a pivot value of $\tau_\alpha$, is calibrated so as to reproduce the $z=5.7$ Ly$\alpha$ luminosity function (LF). We compare our model predictions with Ly$\alpha$ LFs at $z=6.6$ and $7.3$, LAE angular auto-correlation functions (ACFs) at $z=5.7$ and $6.6$, and LAE fractions in Lyman break galaxies at $5<z<7$. The Ly$\alpha$ LFs and ACFs are reproduced by multiple models, but the LAE fraction turns out to be the most critical test. The dispersion of $\tau_\alpha$ and the halo mass dependence of $\tau_\alpha$ are essential to explain all observations reasonably. Therefore, a simple model of one-to-one correspondence between halo mass and Ly$\alpha$ luminosity with a constant Ly$\alpha$ escape fraction has been ruled out. Based on our best model, we present a formula to estimate the intergalactic neutral hydrogen fraction, $x_{\rm HI}$, from the observed Ly$\alpha$ luminosity density at $z\gtrsim6$. We finally obtain $x_{\rm HI}=0.5_{-0.3}^{+0.1}$ as a volume-average at $z=7.3$.

Recently, advanced Laser Interferometer Gravitational-Wave Observatory (aLIGO) has detected black hole (BH) merger events, most of which are sourced by BHs more massive than $30~M_\odot$. Especially, the observation of GW170104 suggests dynamically assembled binaries favoring a distribution of misaligned spins. It has been argued that mergers of unassociated BHs can be engendered through a "chance meeting" in a multiple BH system under gas-rich environments. In this paper, we consider the merger of unassociated BHs, concordant with the massive BH merger events. To that end, we simulate a multiple BH system with a post-Newtonian $N$-body code incorporating gas accretion and general relativistic effects. As a result, we find that gas dynamical friction effectively promotes three-body interaction of BHs in dense gas of $n_\mathrm{gas}\gtrsim 10^6\,\mathrm{cm}^{-3}$, so that BH mergers can take place within $30$ Myr. This scenario predicts an isotropic distribution of spin tilts. In the concordant models with GW150914, the masses of seed BHs are required to be $\gtrsim 25~M_\odot$. The potential sites of such "chance meeting" BH mergers are active galactic nucleus (AGN) disks and dense interstellar clouds. Assuming the LIGO O1, we roughly estimate the event rates for PopI BHs and PopIII BHs in AGN disks to be respectively $\simeq 1-2\,\mathrm{yr}^{-1}$ and $\simeq 1\,\mathrm{yr}^{-1}$. Multiple episodes of AGNs may enhance the rates by roughly an order of magnitude. For massive PopI BHs in dense interstellar clouds, the rate is $\simeq 0.02\,\mathrm{yr}^{-1}$. Hence, high-density AGN disks are a more plausible site for mergers of chance meeting BHs.

We present a novel Lyman alpha (Ly$\alpha$) radiative transfer code, SEURAT, where line scatterings are solved adaptively with the resolution of the smoothed particle hydrodynamics (SPH). The radiative transfer method implemented in SEURAT is based on a Monte Carlo algorithm in which the scattering and absorption by dust are also incorporated. We perform standard test calculations to verify the validity of the code; (i) emergent spectra from a static uniform sphere, (ii) emergent spectra from an expanding uniform sphere, and (iii) escape fraction from a dusty slab. Thereby we demonstrate that our code solves the Ly$\alpha$ radiative transfer with sufficient accuracy. We emphasise that SEURAT can treat the transfer of Ly$\alpha$ photons even in highly complex systems that have significantly inhomogeneous density fields. The high adaptivity of SEURAT is desirable to solve the propagation of Ly$\alpha$ photons in the interstellar medium of young star-forming galaxies like Ly$\alpha$ emitters (LAEs). Thus, SEURAT provides a powerful tool to model the emergent spectra of Ly$\alpha$ emission, which can be compared to the observations of LAEs.

The coagulation of dust particles is a key process in planetesimal formation. However, the radial drift and bouncing barriers are not completely resolved, especially for silicate dust. Since the collision velocities of dust particles are regulated by turbulence in a protoplanetary disk, the turbulent clustering should be properly treated. To that end, direct numerical simulations (DNSs) of the Navier Stokes equations are requisite. In a series of papers, Pan & Padoan used a DNS with the Reynolds number Re~1000. Here, we perform DNSs with up to Re=16100, which allow us to track the motion of particles with Stokes numbers of 0.01<~St<~0.2 in the inertial range. By the DNSs, we confirm that the rms relative velocity of particle pairs is smaller by more than a factor of two, compared to those by Ormel & Cuzzi (2007). The distributions of the radial relative velocities are highly non-Gaussian. The results are almost consistent with those by Pan & Padoan or Pan et al. at low-Re. Also, we find that the sticking rates for equal-sized particles are much higher than those for different-sized particles. Even in the strong-turbulence case with alpha-viscosity of 10^{-2}, the sticking rates are as high as >~50% and the bouncing probabilities are as low as ~10% for equal-sized particles of St<~0.01. Thus, the turbulent clustering plays a significant role for the growth of cm-sized compact aggregates (pebbles) and also enhances the solid abundance, which may lead to the streaming instability in a disk.

Providing initial conditions is an essential procedure for numerical simulations of galaxies. The initial conditions for idealised individual galaxies in $N$-body simulations should resemble observed galaxies and be dynamically stable for time scales much longer than their characteristic dynamical times. However, generating a galaxy model ab initio as a system in dynamical equilibrium is a difficult task, since a galaxy contains several components, including a bulge, disc, and halo. Moreover, it is desirable that the initial-condition generator be fast and easy to use. We have now developed an initial-condition generator for galactic $N$-body simulations that satisfies these requirements. The developed generator adopts a distribution-function-based method, and it supports various kinds of density models, including custom-tabulated inputs and the presence of more than one disc. We tested the dynamical stability of systems generated by our code, representing early- and late-type galaxies, with $N=$~2,097,152 and 8,388,608 particles, respectively, and we found that the model galaxies maintain their initial distributions for at least 1~Gyr. The execution times required to generate the two models were $8.5$ and $221.7$ seconds, respectively, which is negligible compared to typical execution times for $N$-body simulations. The code is provided as open-source software and is publicly and freely available at \url{https://bitbucket.org/ymiki/magi}.

We present the SILVERRUSH program strategy and clustering properties investigated with $\sim 2,000$ Ly$\alpha$ emitters at $z=5.7$ and $6.6$ found in the early data of the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey exploiting the carefully designed narrowband filters. We derive angular correlation functions with the unprecedentedly large samples of LAEs at $z=6-7$ over the large total area of $14-21$ deg$^2$ corresponding to $0.3-0.5$ comoving Gpc$^2$. We obtain the average large-scale bias values of $b_{\rm avg}=4.1\pm 0.2$ ($4.5\pm 0.6$) at $z=5.7$ ($z=6.6$) for $\gtrsim L^*$ LAEs, indicating the weak evolution of LAE clustering from $z=5.7$ to $6.6$. We compare the LAE clustering results with two independent theoretical models that suggest an increase of an LAE clustering signal by the patchy ionized bubbles at the epoch of reionization (EoR), and estimate the neutral hydrogen fraction to be $x_{\rm HI}=0.15^{+0.15}_{-0.15}$ at $z=6.6$. Based on the halo occupation distribution models, we find that the $\gtrsim L^*$ LAEs are hosted by the dark-matter halos with the average mass of $\log (\left < M_{\rm h} \right >/M_\odot) =11.1^{+0.2}_{-0.4}$ ($10.8^{+0.3}_{-0.5}$) at $z=5.7$ ($6.6$) with a Ly$\alpha$ duty cycle of 1 % or less, where the results of $z=6.6$ LAEs may be slightly biased, due to the increase of the clustering signal at the EoR. Our clustering analysis reveals the low-mass nature of $\gtrsim L^*$ LAEs at $z=6-7$, and that these LAEs probably evolve into massive super-$L^*$ galaxies in the present-day universe.

The tree method is a widely implemented algorithm for collisionless $N$-body simulations in astrophysics well suited for GPU(s). Adopting hierarchical time stepping can accelerate $N$-body simulations; however, it is infrequently implemented and its potential remains untested in GPU implementations. We have developed a Gravitational Oct-Tree code accelerated by HIerarchical time step Controlling named \texttt{GOTHIC}, which adopts both the tree method and the hierarchical time step. The code adopts some adaptive optimizations by monitoring the execution time of each function on-the-fly and minimizes the time-to-solution by balancing the measured time of multiple functions. Results of performance measurements with realistic particle distribution performed on NVIDIA Tesla M2090, K20X, and GeForce GTX TITAN X, which are representative GPUs of the Fermi, Kepler, and Maxwell generation of GPUs, show that the hierarchical time step achieves a speedup by a factor of around 3--5 times compared to the shared time step. The measured elapsed time per step of \texttt{GOTHIC} is 0.30~s or 0.44~s on GTX TITAN X when the particle distribution represents the Andromeda galaxy or the NFW sphere, respectively, with $2^{24} =$~16,777,216 particles. The averaged performance of the code corresponds to 10--30\% of the theoretical single precision peak performance of the GPU.

We explore the possibility of the formation of globular clusters under ultraviolet (UV) background radiation. One-dimensional spherical symmetric radiation hydrodynamics (RHD) simulations by Hasegawa et al. have demonstrated that the collapse of low-mass (10^6-10^7 solar masses) gas clouds exposed to intense UV radiation can lead to the formation of compact star clusters like globular clusters (GCs) if gas clouds contract with supersonic infall velocities. However, three-dimensional effects, such as the anisotropy of background radiation and the inhomogeneity in gas clouds, have not been studied so far. In this paper, we perform three-dimensional RHD simulations in a semi-cosmological context, and reconsider the formation of compact star clusters in strong UV radiation fields. As a result, we find that although anisotropic radiation fields bring an elongated shadow of neutral gas, almost spherical compact star clusters can be procreated from a "supersonic infall" cloud, since photo-dissociating radiation suppresses the formation of hydrogen molecules in the shadowed regions and the regions are compressed by UV heated ambient gas. The properties of resultant star clusters match those of GCs. On the other hand, in weak UV radiation fields, dark matter-dominated star clusters with low stellar density form due to the self-shielding effect as well as the positive feedback by ionizing photons. Thus, we conclude that the "supersonic infall" under a strong UV background is a potential mechanism to form GCs.

We present post-Newtonian $N$-body simulations on mergers of accreting stellar-mass black holes (BHs), where such general relativistic effects as the pericenter shift and gravitational wave (GW) emission are taken into consideration. The attention is concentrated on the effects of the dynamical friction and the Hoyle-Lyttleton mass accretion by ambient gas. We consider a system composed of ten BHs with initial mass of $30~M_\odot$. As a result, we show that mergers of accreting stellar-mass BHs are classified into four types: a gas drag-driven, an interplay-driven, a three body-driven, or an accretion-driven merger. We find that BH mergers proceed before significant mass accretion, even if the accretion rate is $\sim10$ Eddington accretion rate, and then all BHs can merge into one heavy BH. Using the simulation results for a wide range of parameters, we derive a critical accretion rate ($\dot{m}_{\rm c}$), below which the BH growth is promoted faster by mergers. Also, it is found that the effect of the recoil by the GW emission can reduce $\dot{m}_{\rm c}$ especially in gas number density higher than $10^8~{\rm cm}^{-3}$, and enhance the escape probability of merged BHs. Very recently, a gravitational wave event, GW150914, as a result of the merger of a $\sim 30~M_\odot$ BH binary has been detected (Abbott et al. 2016). Based on the present simulations, the BH merger in GW150914 is likely to be driven by three-body encounters accompanied by a few $M_\odot$ of gas accretion, in high-density environments like dense interstellar clouds or galactic nuclei.

We explore the gas dynamics near the dust sublimation radius of active galactic nucleus (AGN). For the purpose, we perform axisymmetric radiation hydrodynamic simulations of a dusty gas disk of radius $\approx 1\,\mathrm{pc}$ around a supermassive black hole of mass $10^{7}\,\mathrm{M_{\odot}}$ taking into account (1) anisotropic radiation of accretion disk, (2) X-ray heating by corona, (3) radiative transfer of infrared (IR) photons reemitted by dust, (4) frequency dependency of direct and IR radiations, and (5) separate temperatures for gas and dust. As a result, we find that for Eddington ratio $\approx 0.77$, a nearly neutral, dense ($\approx 10^{6\operatorname{-}8}\;\mathrm{cm^{-3}}$), geometrically-thin ($h/r<0.06$) disk forms with a high velocity ($\approx 200 \sim 3000\;\mathrm{km/s}$) dusty outflow launched from the disk surface. The disk temperature is determined by the balance between X-ray heating and various cooling, and the disk is almost supported by thermal pressure. Contrary to \citet{krolik07:_agn}, the radiation pressure by IR photons is not effective to thicken the disk, but rather compresses it. Thus, it seems difficult for a radiation-supported, geometrically-thick, obscuring torus to form near the dust sublimation radius as far as the Eddington ratio is high ($\sim 1$). The mass outflow rate is $0.05$-$0.1\;\mathrm{M_{\odot}}/\mathrm{yr}$ and the column density of the outflow is $N_{\mathrm{H}}\lesssim 10^{21}\;\mathrm{cm^{-2}}$. To explain observed type-II AGN fraction, it is required that outflow gas is extended to larger radii ($r\gtrsim 10\;\mathrm{pc}$) or that a denser dusty wind is launched from smaller radii ($r\sim 10^{4}\;R_{g}$).

Powerful relativistic jets in radio galaxies are capable of driving strong outflows but also inducing star-formation by pressure-triggering collapse of dense clouds. We review theoretical work on negative and positive active galactic nuclei feedback, discussing insights gained from recent hydrodynamical simulations of jet-driven feedback on galaxy scales that are applicable to compact radio sources. The simulations show that the efficiency of feedback and the relative importance of negative and positive feedback depends strongly on interstellar medium properties, especially the column depth and spatial distribution of clouds. Negative feedback is most effective if clouds are distributed spherically and individual clouds have small column depths, while positive feedback is most effective if clouds are predominantly in a disc-like configuration.

We present statistical properties of diffuse Lyman-alpha halos (LAHs) around high-$z$ star-forming galaxies with large Subaru samples of Lyman-alpha emitters (LAEs) at $z=2.2$. We make subsamples defined by the physical quantities of LAEs' central Lyman-alpha luminosities, UV magnitudes, Lyman-alpha equivalent widths, and UV slopes, and investigate LAHs' radial surface brightness (SB) profiles and scale lengths $r_n$ as a function of these physical quantities. We find that there exist prominent LAHs around LAEs with faint Lyman-alpha luminosities, bright UV luminosities, and small Lyman-alpha equivalent widths in cumulative radial Lyman-alpha SB profiles. We confirm this trend with the anti-correlation between $r_n$ and Lyman-alpha luminosities (equivalent widths) based on the Spearman's rank correlation coefficient that is $\rho=-0.9$ ($-0.7$) corresponding to the $96\%$ ($93\%$) confidence level, although the correlation between $r_n$ and UV magnitudes is not clearly found in the rank correlation coefficient. Our results suggest that LAEs with properties similar to typical Lyman-break galaxies (with faint Lyman-alpha luminosities and small equivalent widths) possess more prominent LAHs. We investigate scenarios for the major physical origins of LAHs with our results, and find that the cold stream scenario is not preferred, due to the relatively small equivalent widths up to $77$\AA\ in LAHs that include LAEs' central components. There remain two possible scenarios of Lyman-alpha scattering in circum-galactic medium and satellite galaxies that cannot be tested with our observational data.

We perform numerical simulations on the merger of multiple black holes (BHs) in primordial gas at early cosmic epochs. We consider two cases of BH mass: $M_{BH} = 30 M_{\odot}$ and $M_{BH} = 10^4 M_{\odot}$. Attention is concentrated on the effect of the dynamical friction by gas in a host object. The simulations incorporate such general relativistic effects as the pericentre shift and gravitational wave emission. As a result, we find that multiple BHs are able to merge into one BH within 100 Myr in a wide range of BH density. The merger mechanism is revealed to be categorized into three types: gas-drag-driven merger (type A), interplay-driven merger (type B), and three-body-driven merger (type C). We find the relation between the merger mechanism and the ratio of the gas mass within the initial BH orbit ($M_{gas}$) to the total BH mass (${\Sigma}M_{BH}$). Type A merger occurs if $M_{gas} \gtrsim 10^5 {\Sigma}M_{BH}$, type B if $M_{gas} \lesssim 10^5 {\Sigma}M_{BH}$, and type C if $M_{gas} \ll 10^5 {\Sigma}M_{BH}$. Supposing the gas and BH density based on the recent numerical simulations on first stars, all the BH remnants from first stars are likely to merge into one BH through the type B or C mechanism. Also, we find that multiple massive BHs ($M_{BH} = 10^4 M_{\odot}$) distributed over several parsec can merge into one BH through the type B mechanism, if the gas density is higher than $5\times 10^6$ cm$^{-3}$. The present results imply that the BH merger may contribute significantly to the formation of supermassive BHs at high redshift epochs.

We perform three-dimensional radiation hydrodynamic simulations of uniform dusty gas clouds irradiated by an active galactic nucleus (AGN) to investigate the dependence of evolution of clouds on the ionization parameter $\mathcal{U}$ and the Str{\"o}mgren number $\mathcal{N}_{S}$. We find that the evolution can be classified into two cases depending on $\mathcal{U}$. In low $\mathcal{U}$ cases ($\mathcal{U}\approx 10^{-2}$), the evolution is mainly driven by photo-evaporation. A approximately spherically-symmetric evaporation flow with velocity of $100\operatorname{-}150\;\mathrm{km\;s^{-1}}$ is launched from the irradiated face. The cloud is compressed by a D-type shock with losing its mass due to photo-evaporation and is finally turned into a dense filament by $t\lesssim 1.5t_{\mathrm{sc}}$. In high $\mathcal{U}$ cases ($\mathcal{U}\approx 5\times 10^{-2}$), radiation pressure suppresses photo-evaporation from the central part of the irradiated face, reducing photo-evaporation rate. A evaporation flow from the outskirts of the irradiated face is turned into a high velocity ($\lesssim 500\;\mathrm{km\;s^{-1}}$) gas wind because of radiation pressure on dust. The cloud is swept by a radiation pressure-driven shock and becomes a dense gas disk by $t\approx t_{\mathrm{sweep}}$. Star formation is expected in these dense regions for both cases of $\mathcal{U}$. We discuss the influences of the AGN radiation on the clumpy torus. A simple estimate suggests that the clumps are destroyed in timescales shorter than their orbital periods. For the clumpy structure to be maintained over long period, the incident radiation field needs to be sufficiently weaken for most of the clumps, or, some mechanism that creates the clumps continuously is needed.

We present a statistical study of velocities of Lya, interstellar (IS) absorption, and nebular lines and gas covering fraction for Lya emitters (LAEs) at z~2. We make a sample of 22 LAEs with a large Lya equivalent width (EW) of > 50A based on our deep Keck/LRIS observations, in conjunction with spectroscopic data from the Subaru/FMOS program and the literature. We estimate the average velocity offset of Lya from a systemic redshift determined with nebular lines to be dv_Lya=234+-9 km s-1. Using a Kolmogorv-Smirnov test, we confirm the previous claim of Hashimoto et al. (2013) that the average dv_Lya of LAEs is smaller than that of LBGs. Our LRIS data successfully identify blue-shifted multiple IS absorption lines in the UV continua of four LAEs on an individual basis. The average velocity offset of IS absorption lines from a systemic redshift is dv_IS=204+-27 km s-1, indicating LAE's gas outflow with a velocity comparable to typical LBGs. Thus, the ratio, R^Lya_ IS = dv_Lya/dv_IS of LAEs, is around unity, suggestive of low impacts on Lya transmission by resonant scattering of neutral hydrogen in the IS medium. We find an anti-correlation between Lya EW and the covering fraction, f_c, estimated from the depth of absorption lines, where f_c is an indicator of average neutral hydrogen column density, N_HI. The results of our study support the idea that N_HI is a key quantity determining Lya emissivity.

We present diffuse Lyman-alpha halos (LAHs) identified in the composite Subaru narrowband images of 100-3600 Lyman-alpha emitters (LAEs) at z=2.2, 3.1, 3.7, 5.7, and 6.6. First, we carefully examine potential artifacts mimicking LAHs that include a large-scale point-spread function (PSF) made by instrumental and atmospheric effects. Based on our critical test with composite images of non-LAE samples whose narrowband-magnitude and source-size distributions are the same as our LAE samples, we confirm that no artifacts can produce a diffuse extended feature similar to our LAHs. After this test, we measure the scale lengths of exponential profile for the LAHs estimated from our z=2.2-6.6 LAE samples of L(Lyman-alpha) > 2 x 10^42 erg s^-1. We obtain the scale lengths of ~ 5-10 kpc at z=2.2-5.7, and find no evolution of scale lengths in this redshift range beyond our measurement uncertainties. Combining this result and the previously-known UV-continuum size evolution, we infer that the ratio of LAH to UV-continuum sizes is nearly constant at z=2.2-5.7. The scale length of our z=6.6 LAH is larger than 5-10 kpc just beyond the error bar, which is a hint that the scale lengths of LAHs would increase from z=5.7 to 6.6. If this increase is confirmed by future large surveys with significant improvements of statistical and systematical errors, this scale length change at z > 6 would be a signature of increasing fraction of neutral hydrogen scattering Lyman-alpha photons, due to cosmic reionization.

We perform N-body simulations on a multiple massive black hole (MBH) system in a host galaxy to derive the criteria for successive MBH merger. The calculations incorporate the dynamical friction by stars and general relativistic effects as pericentre shift and gravitational wave recoil. The orbits of MBHs are pursed down to ten Schwarzschild radii (~ 1 AU). As a result, it is shown that about a half of MBHs merge during 1 Gyr in a galaxy with mass $10^{11}M_{\odot}$ and stellar velocity dispersion 240 km/s, even if the recoil velocity is two times as high as the stellar velocity dispersion. The dynamical friction allows a binary MBH to interact frequently with other MBHs, and then the decay of the binary orbits leads to the merger through gravitational wave radiation, as shown by Tanikawa & Umemura (2011). We derive the MBH merger criteria for the masses, sizes, and luminosities of host galaxies. It is found that the successive MBH mergers are expected in bright galaxies, depending on redshifts. Furthermore, we find that the central stellar density is reduced by the sling-shot mechanism and that high-velocity stars with ~ 1000 km/s are generated intermittently in extremely radial orbits.

We present the results of structure analyses for a large sample of 426 Lya emitters (LAEs) at z~2.2 that are observed with HST/ACS and WFC3-IR by deep extra-galactic legacy surveys. We confirm that the merger fraction and the average ellipticity of LAE's stellar component are 10-30 % and 0.4-0.6, respectively, that are comparable with previous study results. We successfully identify that some LAEs have a spatial offset between Lya and stellar-continuum emission peaks, d_Lya, by ~2.5-4 kpc beyond our statistical errors. To uncover the physical origin of strong Lya emission found in LAEs, we investigate Lya equivalent width (EW) dependences of these three structural parameters, merger fraction, d_Lya, and ellipticity of stellar distribution in the range of EW(Lya)=20-250A. Contrary to expectations, we find that merger fraction does not significantly increase with Lya EW. We reveal an anti-correlation between d_Lya and EW(Lya) by Kolmogorov-Smirnov (KS) test. There is a trend that the LAEs with a large Lya EW have a small ellipticity. This is consistent with the recent theoretical claims that Lya photons can more easily escape from face-on disks having a small ellipticity, due to less inter-stellar gas along the line of sight, although our KS test indicates that this trend is not statistically significant. Our results of Lya-EW dependence generally support the idea that an HI column density is a key quantity determining Lya emissivity.

We present the first systematic search for extended metal-line [OII]{\lambda}{\lambda}3726,3729 nebulae, or [OII] blobs (OIIBs), at z=1.2 using deep narrowband imaging with a survey volume of 1.9x10^5 Mpc^3 on the 0.62 deg^2 sky of Subaru-XMM Deep Survey (SXDS) field. We discover a giant OIIB, dubbed 'OIIB 1', with a spatial extent over ~75 kpc at a spectroscopic redshift of z=1.18, and also identify a total of twelve OIIBs with a size of >30 kpc. Our optical spectrum of OIIB 1 presents [NeV]{\lambda}3426 line at the 6{\sigma} level, indicating that this object harbors an obscured type-2 AGN. The presence of gas outflows in this object is suggested by two marginal detections of FeII{\lambda}2587 absorption and FeII*{\lambda}2613 emission lines both of which are blueshifted at as large as 500-600 km/s, indicating that the heating source of OIIB 1 is AGN or associated shock excitation rather than supernovae produced by starbursts. The number density of OIIB 1-type giant blobs is estimated to be ~5x10^{-6} Mpc^{-3} at z~1.2, which is comparable with that of AGNs driving outflow at a similar redshift, suggesting that giant OIIBs are produced only by AGN activity. On the other hand, the number density of small OIIBs, 6x10^{-5} Mpc^{-3}, compared to that of z~1 galaxies in the blue cloud in the same M_B range, may imply that 3% of star-forming galaxies at z~1 are quenching star formation through outflows involving extended [OII] emission.

We show, using global 3D grid-based hydrodynamical simulations, that Ultra Fast Outflows (UFOs) from Active Galactic Nuclei (AGN) result in considerable feedback of energy and momentum into the interstellar medium (ISM) of the host galaxy. The AGN wind interacts strongly with the inhomogeneous, two-phase ISM consisting of dense clouds embedded in a tenuous hot hydrostatic medium. The outflow floods through the inter-cloud channels, sweeps up the hot ISM, and ablates and disperses the dense clouds. The momentum of the UFO is primarily transferred to the dense clouds via the ram pressure in the channel flow, and the wind-blown bubble evolves in the energy-driven regime. Any dependence on UFO opening angle disappears after the first interaction with obstructing clouds. On kpc scales, therefore, feedback by UFOs operates similarly to feedback by relativistic AGN jets. Negative feedback is significantly stronger if clouds are distributed spherically, rather than in a disc. In the latter case the turbulent backflow of the wind drives mass inflow toward the central black hole. Considering the common occurrence of UFOs in AGN, they are likely to be important in the cosmological feedback cycles of galaxy formation.

First, the formation of first objects driven by dark matter is revisited by high-resolution hydrodynamic simulations. It is revealed that dark matter haloes of ~10^4M_sun can produce first luminous objects with the aid of dark matter cusps. Therefore, the mass of first objects is smaller by roughly two orders of magnitude than in the previous prediction. This implies that the number of Pop III stars formed in the early universe could be significantly larger than hitherto thought. Secondly, the feedback by photo-ionization and photo-dissociation photons in the first objects is explored with radiation hydrodynamic simulations, and it is demonstrated that multiple stars can form in a 10^5M_sun halo. Thirdly, the fragmentation of an accretion disk around a primordial protostar is explored with photo-dissociation feedback. As a result, it is found that the photo-dissociation can reduce the mass accretion rate onto protostars. Also, protostars as small as 0.8M_sun may be ejected and evolve with keeping their mass, which might be detected as "real first stars" in the Galactic halo. Finally, state-of-the-art radiation hydrodynamic simulations are performed to investigate the internal ionization of first galaxies and the escape of ionizing photons. We find that UV feedback by forming massive stars enhances the escape fraction even in a halo as massive as > 6* 10^9M_sun, while it reduces the star formation rate significantly. This may have a momentous impact on the cosmic reionization.

We examine the detailed physics of the feedback mechanism by relativistic AGN jets interacting with a two-phase fractal interstellar medium in the kpc-scale core of galaxies using 29 3D grid-based hydrodynamical simulations. The feedback efficiency, as measured by the amount of cloud-dispersal generated by the jet-ISM interactions, is sensitive to the maximum size of clouds in the fractal cloud distribution but not to their volume filling factor. Feedback ceases to be efficient for Eddington ratios P_jet/L_edd<10^-4, although systems with large cloud complexes ~50 pc require jets of Eddington ratio in excess of 10^-2 to disperse the clouds appreciably. Based on measurements of the bubble expansion rates in our simulations we argue that sub-grid AGN prescriptions resulting in negative feedback in cosmological simulations without a multi-phase treatment of the ISM are good approximations if the volume filling factor of warm phase material is less than 0.1 and the cloud complexes are smaller than ~25 pc. We find that the acceleration of the dense embedded clouds is provided by the ram pressure of the high velocity flow through the porous channels of the warm phase, flow that has fully entrained the shocked hot-phase gas it has swept up, and is additionally mass-loaded by ablated cloud material. This mechanism transfers 10% to 40% of the jet energy to the cold and warm gas, accelerating it within 10 to 100 Myr to velocities that match those observed in a range of high and low redshift radio galaxies hosting powerful radio jets.

We present a scheme for numerical simulations of collisionless self-gravitating systems which directly integrates the Vlasov--Poisson equations in six-dimensional phase space. By the results from a suite of large-scale numerical simulations, we demonstrate that the present scheme can simulate collisionless self-gravitating systems properly. The integration scheme is based on the positive flux conservation method recently developed in plasma physics. We test the accuracy of our code by performing several test calculations including the stability of King spheres, the gravitational instability and the Landau damping. We show that the mass and the energy are accurately conserved for all the test cases we study. The results are in good agreement with linear theory predictions and/or analytic solutions. The distribution function keeps the property of positivity and remains non-oscillatory. The largest simulations are run on 64^6 grids. The computation speed scales well with the number of processors, and thus our code performs efficiently on massively parallel supercomputers.

Using stacks of Ly-a images of 2128 Ly-a emitters (LAEs) and 24 protocluster UV-selected galaxies (LBGs) at z=3.1, we examine the surface brightness profiles of Ly-a haloes around high-z galaxies as a function of environment and UV luminosity. We find that the slopes of the Ly-a radial profiles become flatter as the Mpc-scale LAE surface densities increase, but they are almost independent of the central UV luminosities. The characteristic exponential scale lengths of the Ly-a haloes appear to be proportional to the square of the LAE surface densities (r(Lya) \propto Sigma(LAE)^2). Including the diffuse, extended Ly-a haloes, the rest-frame Ly-a equivalent width of the LAEs in the densest regions approaches EW_0(Lya) ~ 200 A, the maximum value expected for young (< 10^7 yr) galaxies. This suggests that Ly-a photons formed via shock compression by gas outflows or cooling radiation by gravitational gas inflows may partly contribute to illuminate the Ly-a haloes; however, most of their Ly-a luminosity can be explained by photo-ionisation by ionising photons or scattering of Ly-a photons produced in HII regions in and around the central galaxies. Regardless of the source of Ly-a photons, if the Ly-a haloes trace the overall gaseous structure following the dark matter distributions, it is not surprising that the Ly-a spatial extents depend more strongly on the surrounding Mpc-scale environment than on the activities of the central galaxies.

We present the results of the observations of the Ly\alpha\ line profiles of 91 emission-line galaxies at z=3.1 with the spectral resolution of \lambda/\delta\lambda (FWHM) = 1700, or 180 km/s. A significant fraction, ~50% of the observed objects show the characteristic double peaks in their Ly\alpha profile. The red peak is much stronger than the blue one for most of the cases. The red peaks themselves also show weak but significant asymmetry and their widths are correlated with the velocity separation of the red and the blue peaks, which implies that the peaks are not isolated multiple components with different velocities but the parts of the single line which is modified by the absorption and/or scattering by the associated neutral hydrogen gas. The characteristic profile can be naturally explained by the scattering in the expanding shell of neutral hydrogen surrounding the Ly\alpha\ emitting region while the attenuation by the inter-galactic medium should also be considered. Our results suggest that the star-formation in these Ly\alpha\ emitters are dominated by the young burst-like events which produce the intrinsic Ly\alpha\ emission as well as the gas outflow.

We present the results of the extensive narrow-band survey of Ly\alpha\ emission-line objects at z=3.1 in the 1.38 deg^2 area surrounding the high density region of star-forming galaxies at z=3.09 in the SSA22 field, as well as in the 1.04 deg^2 area of the three separated general blank fields. In total of 2161 Ly\alpha\ emitters, 1394 in the SSA22 fields and 767 in the general fields, respectively, are detected to the narrow-band AB magnitude limit of 25.73, which corresponds to the line flux of 1.8 x 10^{-17} erg s^{-1} cm^{-2} or luminosity of 1.5 x 10^{42} erg s^{-1} at z=3.1, above the observed equivalent width threshold, 190\AA . The average surface number density of the emitters at z=3.1 in the whole general fields above the thresholds is 0.20+-0.01 arcmin^{-2}. The SSA22 high-density region at z=3.09 whose peak local density is 6 times the average is found to be the most prominent outstanding structure in the whole surveyed area and is firmly identified as a robust `protocluster' with the enough large sample. We also compared the overdensity of the 100 arcmin^2 and 700 arcmin^2 areas which contain the protocluster with the expected fluctuation of the dark matter as well as those of the model galaxies in cosmological simulations. We found that the peak height values of the overdensity correspond to be 8-10 times and 3-4 times of the expected standard deviations of the counts of Ly\alpha emitters at z=3.1 in the corresponding volume, respectively. We conclude that the structure at z=3.09 in the SSA22 field is a very significant and rare density peak up to the scale of 60 Mpc.

Based on a three-dimensional model of an early star-forming galaxy, we explore the evolution of the sub-millimeter brightness. The model galaxy is employed from an ultra-high-resolution chemodynamic simulation of a primordial galaxy by Mori & Umemura, where the SFR is ~ 10 Msun/yr at t<0.3 Gyr and several Msun/yr at t>0.3 Gyr. The former phase well reproduces the observed properties of LAEs and the latter does LBGs. We solve the three-dimensional radiative transfer in the clumpy interstellar media in this model galaxy, taking the size distributions of dust grains into account, and calculate the dust temperature as a function of galactic evolutionary time. We find that the clumpiness of interstellar media plays an important role for the sub-millimeter brightness. In the LAE phase, dust grains are concentrated on clumpy star-forming regions that are distributed all over the galaxy, and the grains can effectively absorb UV radiation from stars. As a result, the dust is heated up to T>35 K. In the LBG phase, the continuous supernovae drive dust grains far away from star-forming regions. Then, the grains cannot absorb much radiation from stars, and becomes into a cold state close to the CMB temperature. Consequently, the dust temperature decreases with the evolutionary time, where the mass-weighted mean temperature is T=26 K at t=0.1 Gyr and T=21 K at t=1.0 Gyr. By this analysis, it turns out that the sub-millimeter brightness is higher in the LAE phase than that in the LBG phase, although the dust-to-gas ratio increases monotonically as a function of time. We derive the spectral energy distributions by placing the model galaxy at a given redshift. The peak flux at 850 micron is found to be S_850 ~ 0.2 - 0.9 mJy if the model galaxy is placed at 6>z>2. This means that ALMA can detect an early star-forming galaxy with SFR of ~ 10 Msun/yr by less than one hour integration with 16 antennas.

We present two types of numerical prescriptions that accelerate the radiative transfer calculation around point sources within a three-dimensional Cartesian grid by using the oct-tree structure for the distribution of radiation sources. In one prescription, distant radiation sources are grouped as a bright extended source when the group's angular size, $\theta_{\rm s}$, is smaller than a critical value, $\theta_{\rm crit}$, and radiative transfer is solved on supermeshes whose angular sizes are similar to that of the group of sources. The supermesh structure is constructed by coarse-graining the mesh structure. With this method, the computational time scales with $N_{\rm m} \log(N_{\rm m}) \log(N_{\rm s})$ where $N_{\rm m}$ and $N_{\rm s}$ are the number of meshes and that of radiation sources, respectively. While this method is very efficient, it inevitably overestimates the optical depth when a group of sources acts as an extended powerful radiation source and affects distant meshes. In the other prescription, a distant group of sources is treated as a bright point source ignoring the spatial extent of the group and the radiative transfer is solved on the meshes rather than the supermeshes. This prescription is simply a grid-based version of {\scriptsize START} by Hasegawa & Umemura and yields better results in general with slightly more computational cost ($\propto N_{\rm m}^{4/3} \log(N_{\rm s})$) than the supermesh prescription. Our methods can easily be implemented to any grid-based hydrodynamic codes and are well-suited to the adaptive mesh refinement methods.

Using highly-accurate $N$-body simulations, we explore the evolution of multiple massive black holes (hereafter, MBHs) in a primordial galaxy that is composed of stars and MBHs. The evolution is pursed with a fourth-order Hermite scheme, where not only three-body interaction of MBHs but also dynamical friction by stars are incorporated. Initially, ten MBHs with equal mass of $10^7M_\odot$ are set in a host galaxy with $10^{11}M_\odot$. It is found that 4 - 6 MBHs merge successively within 1 Gyr, emitting gravitational wave radiation. The key process for the successive merger of MBHs is the dynamical friction by field stars, which enhances three-body interactions of MBHs when they enter the central regions of the galaxy. The heaviest MBH always composes a close binary at the galactic center, which shrinks owing to the angular momentum transfer by the third MBH and eventually merges. The angular momentum transfer by the third MBH is due to the sling-shot mechanism. We find that the secular Kozai mechanism does not work for a binary to merge if we include the relativistic pericenter shift. The simulations show that a multiple MBH system can produce a heavier MBH at the galactic center purely through $N$-body process. This merger path can be of great significance for the growth of MBHs in a primordial galaxy. The merger of multiple MBHs may be a potential source of gravitational waves for the Laser Interferometer Space Antenna (LISA) and pulsar timing.

We present results of a survey for giant Ly-alpha nebulae (LABs) at z=3 with Subaru/Suprime-Cam. We obtained Ly-alpha imaging at z=3.09+-0.03 around the SSA22 protocluster and in several blank fields. The total survey area is 2.1 square degrees, corresponding to a comoving volume of 1.6 x 10^6 Mpc^3. Using a uniform detection threshold of 1.4 x 10^{-18} erg s^{-1} cm^{-2} arcsec^{-2} for the Ly-alpha images, we construct a sample of 14 LAB candidates with major-axis diameters larger than 100 kpc, including five previously known blobs and two known quasars. This survey triples the number of known LABs over 100 kpc. The giant LAB sample shows a possible "morphology-density relation": filamentary LABs reside in average density environments as derived from compact Ly-alpha emitters, while circular LABs reside in both average density and overdense environments. Although it is hard to examine the formation mechanisms of LABs only from the Ly-alpha morphologies, more filamentary LABs may relate to cold gas accretion from the surrounding inter-galactic medium (IGM) and more circular LABs may relate to large-scale gas outflows, which are driven by intense starbursts and/or by AGN activities. Our survey highlights the potential usefulness of giant LABs to investigate the interactions between galaxies and the surrounding IGM from the field to overdense environments at high-redshift.

We present a novel radiation hydrodynamics code, START, which is a smoothed particle hydrodynamics (SPH) scheme coupled with accelerated radiative transfer. The basic idea for the acceleration of radiative transfer is parallel to the tree algorithm that is hitherto used to speed up the gravitational force calculation in an N-body system. It is demonstrated that the radiative transfer calculations can be dramatically accelerated, where the computational time is scaled as Np log Ns for Np SPH particles and Ns radiation sources. Such acceleration allows us to readily include not only numerous sources but also scattering photons, even if the total number of radiation sources is comparable to that of SPH particles. Here, a test simulation is presented for a multiple source problem, where the results with START are compared to those with a radiation SPH code without tree-based acceleration. We find that the results agree well with each other if we set the tolerance parameter as < 1.0, and then it demonstrates that START can solve radiative transfer faster without reducing the accuracy. One of important applications with START is to solve the transfer of diffuse ionizing photons, where each SPH particle is regarded as an emitter. To illustrate the competence of START, we simulate the shadowing effect by dense clumps around an ionizing source. As a result, it is found that the erosion of shadows by diffuse recombination photons can be solved. Such an effect is of great significance to reveal the cosmic reionization process.

From multi-wavelength observations of LAEs,we know that while many LAEs appear to be young and less massive,a noticeable fraction of LAEs possess much older populations of stars and larger stellar mass.How these two classes of LAEs are concordant with the hierarchical galaxy formation scenario has not been understood clearly so far.In this paper,we model LAEs by three-dimensional cosmological simulations of dark halo merger in a CDM universe.As a result,it is shown that the age of simulated LAEs can spread over a wide range from 2*10^6yr to 9*10^8yr.Also,we find that there are two types of LAEs, in one of which the young half-mass age is comparable to the mean age of stellar component,and in the other of which the young half-mass age is appreciably shorter than the mean age.We define the former as Type 1 LAEs and the latter as Type 2 LAEs.A Type 1 corresponds to early starburst in a young galaxy,whereas a Type 2 does to delayed starburst in an evolved galaxy,as a consequence of delayed accretion of a subhalo onto a larger parent halo.Thus,the same halo can experience a Type 2 LAE-phase as well as a Type 1 LAE-phase in the merger history.Type 1s are expected to be younger than 1.5*10^8yr,less dusty,and less massive with stellar mass M*<10^8 Msun,while Type 2s are older than 1.5*10^8yr,even dustier,and as massive as M*~10^8-10^10Msun.The fraction of Type 2s in all LAEs is a function of redshift.Type 2s discriminated clearly from Type 1s in two color diagram of z'-H vs J-K.We find that the brightness distribution of LyA in Type 2s is more extended than the main stellar component,in contrast to Type 1s.This is not only because delayed starbursts tend to occur in the outskirts of a parent galaxy,but also because LyA photons are effectively absorbed by dust in an evolved galaxy.Hence,the extent of LyA emission may be an additional measure to distinguish Type 2s from Type 1s

We explore the formation of secondary Population III (Pop III) stars under radiation hydrodynamic (RHD) feedback by a preformed massive star. To properly treat RHD feedback, we perform three-dimensional RHD simulations incorporating the radiative transfer of ionizing photons as well as H_2 dissociating photons from a preformed star. A collapsing gas cloud is settled at a given distance from a 120Msun Pop III star, and the evolution of the cloud is pursued including RHD feedback. We derive the threshold density depending on the distance, above which the cloud can keep collapsing owing to the shielding of H_2 dissociating radiation. We find that an H_2 shell formed ahead of an ionizing front works effectively to shield the H_2 dissociating radiation, leading to the positive feedback for the secondary Pop III star formation. Also, near the threshold density, the envelope of gas cloud is stripped significantly by a shock associated with an ionizing front. By comparing the mass accretion timescale with the Kelvin-Helmholtz timescale, we estimate the mass of secondary Pop III stars. It turns out that the stripping by a shock can reduce the mass of secondary Pop III stars down to \approx 20Msum.

We present a novel scenario for globular cluster (GC) formation, where the ultraviolet (UV) background radiation effectively works so as to produce compact star clusters. Here, we explore the formation of GCs in UV radiation fields. For this purpose, we calculate baryon and dark matter (DM) dynamics in spherical symmetry, incorporating the self-shielding effects by solving the radiative transfer of UV radiation. In addition, we prescribe the star formation in cooled gas components and pursue the dynamics of formed stars. As a result, we find that the evolution of subgalactic objects in UV background radiation are separated into three types, that is, (1) prompt star formation, where less massive clouds ~10^{5-8} M_sun are promptly self-shielded and undergo star formation, (2) delayed star formation, where photoionized massive clouds >10^8 M_sun collapse despite high thermal pressure and are eventually self-shielded to form stars in a delayed fashion, and (3) supersonic infall, where photoionized less massive clouds ~10^{5-8} M_sun contract with supersonic infall velocity and are self-shielded when a compact core forms. In particular, the type (3) is a novel type found in the present simulations, and eventually produces a very compact star cluster. The resultant mass-to-light ratios, half-mass radii, and velocity dispersions for the three types are compared to the observations of GCs, dwarf spheroidals (dSphs), and ultra-compact dwarfs (UCDs). It turns out that the properties of star clusters resulting from supersonic infall match well with those of observed GCs, whereas the other two types are distinct from GCs. Hence, we conclude that supersonic infall in a UV background is a promising mechanism to form GCs.

By solving radiative transfer equations, we examine three-dimensional radiative properties of a magnetohydrodynamic accretion flow model confronting with the observed spectrum of Sgr A*, in the vicinity of supermassive black hole at the Galactic centre. As a result, we find that the core of radio emission is larger than the size of the event horizon shadow and its peak location is shifted from the gravitational centre. We also find that the self-absorbed synchrotron emissions by the superposition of thermal electrons within a few tens of the Schwartzschild radius can account for low-frequency spectra below the critical frequency $\nu_{c}\approx 10^{12}$ Hz. Above the critical frequency, the synchrotron self-Compton emission by thermal electrons can account for variable emissions in recent near-infrared observations. In contrast to the previous study by Ohsuga et al. (2005), we found that the X-ray spectra by Bremsstrahlung emission of thermal electrons for the different mass accretion rates can be consistent with both the flaring state and the quiescent state of Sgr A* observed by {\it Chandra}.

In order to assess the contribution of Lyman break galaxies (LBGs) and Lyman alpha emitters (LAEs) at redshifts 3<z<7 to the ionization of intergalactic medium (IGM), we investigate the escape fractions of ionizing photons from supernova-dominated primordial galaxies by solving the three-dimensional radiative transfer. The model galaxy is employed from an ultra-high-resolution chemodynamic simulation of a primordial galaxy by Mori & Umemura (2006),which well reproduces the observed properties of LAEs and LBGs. The total mass of model galaxy is 10^{11}M_sun. We solve not only photo-ionization but also collisional ionization by shocks. In addition, according to the chemical enrichment, we incorporate the effect of dust extinction, taking the size distributions of dust into account. As a result, we find that dust extinction reduces the escape fractions by a factor 1.5-8.5 in the LAE phase and by a factor 2.5-11 in the LBG phase, while the collisional ionization by shocks increases the escape fractions by a factor ~ 2. The resultant escape fractions are 0.07-0.47 in the LAE phase and 0.06-0.17 in the LBG phase. These results are well concordant with the recent estimations derived from the flux ratio at 1500 angstrom to 900 angstrom of LAEs and LBGs. Combining the resultant escape fractions with the luminosity functions of LAEs and LBGs, we find that high-z LAEs and LBGs can ionize the IGM at z=3-5. However, ionizing radiation from LAEs as well as LBGs falls short to ionize the IGM at z>6. That implies that additional ionization sources may required at z>6.

The development of radiation hydrodynamical methods that are able to follow gas dynamics and radiative transfer self-consistently is key to the solution of many problems in numerical astrophysics. Such fluid flows are highly complex, rarely allowing even for approximate analytical solutions against which numerical codes can be tested. An alternative validation procedure is to compare different methods against each other on common problems, in order to assess the robustness of the results and establish a range of validity for the methods. Previously, we presented such a comparison for a set of pure radiative transfer tests (i.e. for fixed, non-evolving density fields). This is the second paper of the Cosmological Radiative Transfer (RT) Comparison Project, in which we compare 9 independent RT codes directly coupled to gasdynamics on 3 relatively simple astrophysical hydrodynamics problems: (5) the expansion of an H II region in a uniform medium; (6) an ionization front (I-front) in a 1/r^2 density profile with a flat core, and (7), the photoevaporation of a uniform dense clump. Results show a broad agreement between the different methods and no big failures, indicating that the participating codes have reached a certain level of maturity and reliability. However, many details still do differ, and virtually every code has showed some shortcomings and has disagreed, in one respect or another, with the majority of the results. This underscores the fact that no method is universal and all require careful testing of the particular features which are most relevant to the specific problem at hand.

We explore the impact of ultraviolet (UV) radiation from massive Population III (Pop III) stars of 25, 40, 80, and 120 M_sun on the subsequent Pop III star formation. In this paper, particular attention is paid to the dependence of radiative feedback on the mass of source Pop III star. UV radiation from the source star can work to impede the secondary star formation through the photoheating and photodissociation processes. Recently, Susa & Umemura (2006) have shown that the ionizing radiation alleviates the negative effect by H_2-dissociating radiation from 120$M_sun PopIII star, since an H_2 shell formed ahead of an ionizing front can effectively shield H_2-dissociating radiation. On the other hand, it is expected that the negative feedback by H_2-dissociating radiation can be predominant if a source star is less massive, since a ratio of the H_2-dissociating photon number to the ionizing photon number becomes higher. In order to investigate the radiative feedback effects from such less massive stars, we perform three-dimensional radiation hydrodynamic simulations, incorporating the radiative transfer effect of ionizing and H_2-dissociating radiation. As a result, we find that if a source star is less massive than ~25M_sun, the ionizing radiation cannot suppress the negative feedback of H_2-dissociating radiation. Therefore, the fate of the neighboring clouds around such less massive stars is determined solely by the flux of H_2-dissociating radiation from source stars. With making analytic estimates of H_2 shell formation and its shielding effect, we derive the criteria for radiation hydrodynamic feedback depending on the source star mass.

We consider the effect of a supernova (SN) explosion in a very massive binary that is expected to form in a portion of Population III stars with the mass higher than 100$M_\odot$. In a Population III binary system, a more massive star can result in the formation of a BH and a surrounding accretion disc. Such BH accretion could be a significant source of the cosmic reionization in the early universe. However, a less massive companion star evolves belatedly and eventually undergoes a SN explosion, so that the accretion disc around a BH might be blown off in a lifetime of companion star. In this paper, we explore the dynamical impact of a SN explosion on an accretion disc around a massive BH, and elucidate whether the BH accretion disc is totally demolished or not. For the purpose, we perform three-dimensional hydrodynamic simulations of a very massive binary system, where we assume a BH of $10^3 M_{\odot}$ that results from a direct collapse of a very massive star and a companion star of $100 M_{\odot}$ that undergoes a SN explosion. We calculate the remaining mass of a BH accretion disc as a function of time. As a result, it is found that a significant portion of gas disc can survive through three-dimensional geometrical effects even after the SN explosion of a companion star. Even if the SN explosion energy is higher by two orders of magnitude than the binding energy of gas disc, about a half of disc can be left over. The results imply that the Population III BH accretion disc can be a long-lived luminous source, and therefore could be an important ionizing source in the early universe.

The current understanding of galaxy formation is that it proceeds in a 'bottom up' way, with the formation of small clumps of gas and stars that merge hierarchically until giant galaxies are built up. The baryonic gas loses the thermal energy by radiative cooling and falls towards the centres of the new galaxies, while supernovae (SNe) blow gas out. Any realistic model therefore requires a proper treatment of these processes, but hitherto this has been far from satisfactory. Here we report an ultra-high-resolution simulation that follows evolution from the earliest stages of galaxy formation through the period of dynamical relaxation. The bubble structures of gas revealed in our simulation ($< 3\times10^8$ years) resemble closely the high-redshift Lyman $\alpha$ emitters (LAEs). After $10^9$ years these bodies are dominated by stellar continuum radiation and look like the Lyman break galaxies (LBGs) known as the high-redshift star-forming galaxies at which point the abundance of elements heavier than helium ("metallicity") appears to be solar. After $1.3\times10^{10}$ years, these galaxies resemble present-day ellipticals.

Recently discovered large-scale structure of Ly\alpha Emitters (LAEs) raises a novel challenge to the cold dark matter (CDM) cosmology. The structure is extended over more than 50 Mpc at redshift z=3.1, and exhibits a considerably weak angular correlation. Such properties of LAE distributions appear to be incompatible with the standard biased galaxy formation scenario in the CDM cosmology. In this paper, by considering the possibility that LAEs are short-lived events, we attempt to build up the picture of LAEs concordant with the CDM cosmology. We find that if the lifetime of LAEs is as short as (6.7 \pm 0.6) \times 10^7 yr, the distributions of simulated galaxies successfully match the extension and morphology of large-scale structure of LAEs at z=3.1, and also the weak angular correlation function. This result implies that LAEs at z=3.1 do not necessarily reside in high density peaks, but tends to be located in less dense regions, in a different way from the expectation by the standard biased galaxy formation scenario. In addition, we make a prediction for the angular correlation function of LAEs at redshifts higher than 3. It is found that the prediction deviates from that by the standard biased galaxy formation scenario even at redshifts 4 < z < 6.

We report an ultra-high-resolution simulation that follows evolution from the earliest stages of galaxy formation through the period of dynamical relaxation. The bubble structures of gas revealed in our simulation ($< 3\times 10^8$ years) resemble closely the high-redshift Lyman $\alpha$ emitters (LAEs). After $10^9$ years these bodies are dominated by stellar continuum radiation and look like the Lyman break galaxies (LBGs) known as the high-redshift star-forming galaxies at which point the abundance of elements appears to be solar. After $1.3\times10^{10}$ years, these galaxies resemble present-day ellipticals. The comparisons of simulation results with the observations of elliptical galaxies allow us to conclude that LAEs and LBGs are infants of elliptical galaxies or bulge systems in the nearby universe.

We propose a novel method to extract the imprint of gravitational lensing by Pop III stars in the light curves of Gamma Ray Bursts (GRBs). Significant portions of GRBs can originate in hypernovae of Pop III stars and be gravitationally lensed by foreground Pop III stars or their remnants. If the lens mass is on the order of $10^2-10^3M_\odot$ and the lens redshift is greater than 10, the time delay between two lensed images of a GRB is $\approx 1$s and the image separation is $\approx 10 \mu$as. Although it is difficult to resolve the two lensed images spatially with current facilities, the light curves of two images are superimposed with a delay of $\approx 1$ s. GRB light curves usually exhibit noticeable variability, where each spike is less than 1s. If a GRB is lensed, all spikes are superimposed with the same time delay. Hence, if the autocorrelation of light curve with changing time interval is calculated, it should show the resonance at the time delay of lensed images. Applying this autocorrelation method to GRB light curves which are archived as the {\it BATSE} catalogue, we demonstrate that more than half light curves can show the recognizable resonance, if they are lensed. Furthermore, in 1821 GRBs we actually find one candidate of GRB lensed by a Pop III star, which may be located at redshift 20-200. The present method is quite straightforward and therefore provides an effective tool to search for Pop III stars at redshift greater than 10. Using this method, we may find more candidates of GRBs lensed by Pop III stars in the data by the {\it Swift} satellite.

We report an extensive search for Lyman-alpha emitters (LAEs) at z=6.5 in the Subaru Deep Field. Subsequent spectroscopy with Subaru and Keck identified eight more LAEs, giving a total of 17 spectroscopically confirmed LAEs at z=6.5. Based on this spectroscopic sample of 17, complemented by a photometric sample of 58 LAEs, we have derived a more accurate Lyman-alpha luminosity function of LAEs at z=6.5, which reveals an apparent deficit at the bright end of ~0.75 mag fainter L*, compared with that observed at z=5.7. The difference in the LAE luminosity functions between z=5.7 and 6.5 is significant at the 3-sigma level, which is reduced to 2-sigma when cosmic variance is taken into account. This result may imply that the reionization of the universe has not been completed at z=6.5. We found that the spatial distribution of LAEs at z=6.5 was homogeneous over the field. We discuss the implications of these results for the reionization of the universe.

We explore the possibility of subsequent star formation after a first star forms in a Pop III object, by focusing on the radiation hydrodynamic (RHD) feedback brought by ionizing photons as well as H2 dissociating photons. For the purpose, we perform three-dimensional RHD simulations, where the radiative transfer of ionizing photons and H2 dissociating photons from a first star is self-consistently coupled with hydrodynamics based on a smoothed particle hydrodynamics method. As a result, it is shown that density peaks above a threshold density can keep collapsing owing to the shielding of H2 dissociating radiation by an H2 shell formed ahead of a D-type ionization front. But, below the threshold density, an M-type ionization front accompanied by a shock propagates, and density peaks are radiation hydrodynamically evaporated by the shock. The threshold density is dependent on the distance from a source star, which is $\approx 10^2 cm^{-3}$ for the source distance of 30pc. Taking into consideration that the extent of a Pop III object is $\approx 100$pc and density peaks within it have the density of $10^{2-4}$cm$^{-3}$, it is concluded that the secondary star formation is allowed in the broad regions in a Pop III object.

Radiative transfer simulations are now at the forefront of numerical astrophysics. They are becoming crucial for an increasing number of astrophysical and cosmological problems; at the same time their computational cost has come to the reach of currently available computational power. Further progress is retarded by the considerable number of different algorithms (including various flavours of ray-tracing and moment schemes) developed, which makes the selection of the most suitable technique for a given problem a non-trivial task. Assessing the validity ranges, accuracy and performances of these schemes is the main aim of this paper, for which we have compared 11 independent RT codes on 5 test problems: (0) basic physics, (1) isothermal H II region expansion and (2) H II region expansion with evolving temperature, (3) I-front trapping and shadowing by a dense clump, (4) multiple sources in a cosmological density field. The outputs of these tests have been compared and differences analyzed. The agreement between the various codes is satisfactory although not perfect. The main source of discrepancy appears to reside in the multi-frequency treatment approach, resulting in different thicknesses of the ionized-neutral transition regions and different temperature structure. The present results and tests represent the most complete benchmark available for the development of new codes and improvement of existing ones. To this aim all test inputs and outputs are made publicly available in digital form.

We examine whether the ultraluminous infrared galaxies that contain a type I Seyfert nucleus (a type I ULIRG) are in the transition stage from ULIRGs to quasi-stellar objects (QSOs). To inspect this issue, we compare the black hole (BH) mass, the bulge luminosity and the far infrared luminosity among type I ULIRGs, QSOs and elliptical galaxies. As a result, we find the following results; (1) The type I ULIRGs have systematically smaller BH masses in spite of the comparable bulge luminosity relative to QSOs and elliptical galaxies. (2) The far-infrared luminosity of most type I ULIRGs is larger than the Eddington luminosity. We show that above results do not change significantly for 3 type I ULIRGs that we can estimate the visual extinction from the column density. Also, for all 8 type I ULIRGs, we investigate the effect of uncertainties of BH mass measurments and our sample bias, so that it turns out that our results do not alter even if we consider above two effects. In addition, Anabuki (2004) revealed that their X-ray properties are similar to those of the narrow line Seyfert 1 galaxies. These would indicate that active galactic nuclei (AGNs) with a high mass accretion rate exist in the type I ULIRGs. Based on all of these findings, we conclude that it would be a natural interpretation that type I ULIRGs are the early phase of BH growth, namely the missing link between ULIRGs and QSOs. Moreover, by comparing our results with a theoretical model of a coevolution scenario of a QSO BH and a galactic bulge, we show clearly that this explanation would be valid.

Based on the gamma-ray burst (GRB) event rate at redshifts of $4 \leq z \leq 12$, which is assessed by the spectral peak energy-to-luminosity relation recently found by Yonetoku et al., we observationally derive the star formation rate (SFR) for Pop III stars in a high redshift universe. As a result, we find that Pop III stars could form continuously at $4 \leq z \leq 12$. Using the derived Pop III SFR, we attempt to estimate the ultraviolet (UV) photon emission rate at $7 \leq z \leq 12$ in which redshift range no observational information has been hitherto obtained on ionizing radiation intensity. We find that the UV emissivity at $7 \leq z \leq 12$ can make a noticeable contribution to the early reionization. The maximal emissivity is higher than the level required to keep ionizing the intergalactic matter at $7 \leq z \leq 12$. However, if the escape fraction of ionizing photons from Pop III objects is smaller than 10%, then the IGM can be neutralized at some redshift, which may lead to the double reionization. As for the enrichment, the ejection of all metals synthesized in Pop III objects is marginally consistent with the IGM metallicity, although the confinement of metals in Pop III objects can reduce the enrichment significantly.

he formation and growth of supermassive black holes (SMBHs) physically linked with bulges are considered. We focus on the radiation hydrodynamic process for the growth of SMBH in the optically thick starburst phase, where radiation from bulge stars drives the mass accretion on to a galactic center through radiation drag effect. In the present scenario, the AGN luminosity-dominant phase (QSO phase) is preceded by the host luminosity-dominat phase, which is called ``proto-QSO phase''. In this phase, there exists the massive dusty disks within younger bulges. Also, the proto-QSO phase is anticipated by an optically-thick ultraluminous infrared galaxy (ULIRG) phase. Furthermore, such radiation hydrodynamic model has been also applied to disk galaxies. It turns out that the mass of a SMBH primarily correlates with a bulge component even in a disk galaxy. Thus, by analogy to proto-QSOs, the BH growing phase in disk galaxies may have massive dusty disks within younger bulges.

The hydrodynamic formation of massive black holes (BHs) in globular clusters is considered. In particular, we examine the possibility of the BH formation induced by the radiation drag that is exerted on the interstellar matter by stellar radiation in globular clusters. The radiation drag extracts angular momentum from interstellar gas and thus allows the gas to accrete on to the cluster center.By incorporating the realistic chemical evolution of globular clusters, we scrutinize the efficiency of the radiation drag to assess the total accreted mass. As a result, we find that if a globular cluster is more massive than $\approx 6\times 10^{6}M_{\odot}$ in the present-day stellar component, a massive BH with $>260M_{\odot}$ can form within it. But, the BH-to-cluster mass ratio s considerably smaller than the BH-to-bulge mass ratio ($\approx 10^{-3}$) found in galactic bulges. The results are not sensitive to the assumed stellar initial mass function and star formation rate in the cluster, as long as the resultant color-magnitude relation, metallicity, and mass-to-luminosity ratio satisfy those observed in globular clusters. Hence, the putative linear relation between BH mass and bulge mass($M_{\rm BH}-M_{\rm bulge}$ relation) cannot be extrapolated to globular cluster systems. In the present regime, we discuss the BH formation in M15, G1, $\omega$ Cen, the M33 nucleus, and the compact X-ray sources in M82. Finally, we argue observational indications of the formation process of massive BHs in globular clusters. We find that the final phase of BH growth due to the radiation drag can be observed as ultraluminous X-ray sources (ULXs) with $\sim 10^{41}{\rm erg} {\rm s}^{-1}$.

We explore the possibility for the formation of Population III binaries. The collapse of a rotating cylinder is simulated with a three-dimensional, high-resolution nested grid, assuming the thermal history of primordial gas. The simulations are done with dimensionless units, and the results are applicable to low-mass as well as massive systems by scaling with the initial density. We find that if the initial angular momentum is as small as $\beta \approx 0.1$, where $\beta$ is the ratio of centrifugal force to pressure force, then the runaway collapse of the cloud stops to form a rotationally-supported disk. After the accretion of the envelope, the disk undergoes a ring instability, eventually fragmenting into a binary. If the initial angular momentum is relatively large, a bar-type instability arises, resulting in the collapse into a single star through rapid angular momentum transfer. The present results show that a significant fraction of Pop III stars are expected to form in binary systems, even if they are quite massive or less massive. The cosmological implications of Population III binaries are briefly discussed.

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 examine the possibility of the active galactic nucleus (AGN) obscuration by dusty gas clouds that spurt out from circumnuclear starburst regions. For the purpose, the dynamical evolution of gas clouds is pursued, including the effects of radiation forces by an AGN as well as a starburst. Here, we solve the radiative transfer equations for clouds, taking into consideration the growth of clouds by inelastic cloud-cloud collisions and the resultant change in optical depth. As a result, it is shown that if the starburst is more luminous than the AGN, gas clouds are distributed extensively above a galactic disk with the assistance of radiation pressure from the starburst. The total covering factor of gas clouds reaches a maximum of around 20%. After several $10^{7}$yr, gas clouds with larger optical depth form by cloud-cloud collisions and thereafter the clouds fall back due to weakened radiation pressure. The larger clouds undergo runaway growth and are eventually distributed around the equatorial plane on the inner sides of circumnuclear starburst regions. These clouds have an optical depth of several tens. The result is qualitatively consistent with the putative tendency that Seyfert 2 galaxies appear more frequently associated with starbursts than Seyfert 1s. On the other hand, if the AGN luminosity overwhelms that of the starburst, almost all clouds are ejected from the galaxy due to the radiation pressure from the AGN, resulting in the formation of a quasar-like object. The origin of obscuration of AGNs is discussed with relevant observations.

We consider a forming galaxy undergoing multitudinous supernova (SN) explosions, as a possible model of \lya blobs (LABs). For this purpose, an ultra--high resolution hydrodynamic simulation is performed using $1024^3$ grid points, where SN remnants are resolved with sufficient accuracy. It is found that multiple SN explosions produce kpc--size expanding hot bubbles, which drive cool, dense shells by strong shock. The colliding high--density cooling shells radiate intensive \lya emission, resulting in a high \lya luminosity of $\sim10^{43}$ erg s$^{-1}$, comparable to the observed level in LABs. Also, recently discovered bubbly features in some LABs are quite similar to the structure predicted in the present simulation. Furthermore, the result demonstrates that LABs are representative of evolving primordial galaxies; they could hold direct information on the early chemical enrichment of galaxies, contrary to present--day galaxies which have undergone intense recycling of interstellar matter, thus erasing most of the early chemical history. It turns out that the metal mixing proceeds in a very inhomogeneous fashion, so that there appears a large spread of metallicity, that is,$[{\rm Fe/H}] \approx 0 {\rm ~to} -5$ or $[{\rm O/H}] \approx 1 {\rm ~to} -4$. Hence, the early galactic chemical evolution may have proceeded in a different manner from that hitherto considered in one--zone models.

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 perform two-dimensional hydrodynamic simulation on the gas stripping by radiation drag from an interstellar cloud moving in uniform radiation fields. To properly include relativistic radiation drag, the radiation hydrodynamic equation is solved with taking into account the dilution of radiation fields by optical depth of the cloud. As a result, it is found that the optically-thin surface layers are effectively stripped by radiation drag from an optically-thick gas cloud, and simultaneously stripped gas loses momentum. The momentum loss time-scale is found to be on the order of $10^{8}$ years under intensive radiation fields which are expected in the early phase of galaxy evolution. The present results show that the radiation drag is an effective mechanism to extract angular momentum from interstellar medium and allows it to accrete onto the galactic center. The mass accretion driven by radiation drag may lead to the formation of a central supermassive black hole.

Recent observations on the cosmic microwave background by Wilkinson Microwave Anisotropy Probe(WMAP) strongly suggest that the reionization of the universe took place quite early (z \sim 17). On the other hand, it has been pointed out that the cold dark matter cosmology suffers from the substructure problem that subgalactic halos are overproduced than the observed dwarf galaxies in the Local Group. In this paper, as a potential mechanism to solve this problem, we consider the feedback effects of the early reionization on the formation of small-scale structures. For this purpose, we perform 3D radiation hydrodynamic simulations with incorporating the radiative transfer for ionizing photons. As a result, it is found that the early reionization is so devastating for low-mass systems with M_{vir}\la 10^8 M_\odot or v_{circ}\la 20km/s, and almost all gas is photo-evaporated in more than 95% of low-mass systems. Such a strong negative feedback on the formation of low-mass galaxies may solve the substructure problem and support the picture that Local Group dwarf galaxies are descendants of the more massive halos that experienced and survived tidal stripping.

We study the formation and evolution of HII regions around the first stars formed at z=10-30. We use a one-dimensional Lagrangian hydrodynamics code which self-consistently incorporates radiative transfer and non-equilibrium primordial gas chemistry. The star-forming region is defined as a spherical gas cloud with a Population III star embedded at the center. We explore a large parameter space of host halo mass, gas density profile, and stellar luminosity. The formation of the HII region is characterized by initial slow expansion of a weak D-type ionization front near the center, followed by rapid propagation of an R-type front throughout the outer gas envelope. We find that the transition between the two front types is indeed a critical condition for the complete ionization of halos of cosmological interest. In small mass (< 10^6 M_sun) halos, the transition takes place within a few 10^5 yr, yielding high escape fractions (>80%) of both ionizing and photodissociating photons. In larger mass (> 10^7 M_sun) halos, the ionization front remains to be of D-type over the lifetime of the massive star, the HII region is confined well inside the virial radius, and the escape fractions are essentially zero. We derive an analytic formula that reproduces well the results of our simulations. We discuss immediate implications of the present results for the star formation history and early reionization of the Universe.

A potential mechanism is proposed to account for the fact that supermassive black holes (SMBHs) in disk galaxies appear to be smaller than those in elliptical galaxies in the same luminosity range. We consider the formation of SMBHs by radiation drag (Poynting-Robertson effect), which extracts angular momentum from interstellar matter and thereby drives the mass accretion onto a galactic center. In particular, we quantitatively scrutinize the efficiency of radiation drag in a galaxy composed of bulge and disk, and elucidate the relation between the final mass of SMBH and the bulge-to-disk ratio of the galaxy. As a result, it is found that the BH-to-galaxy mass ratio, $M_{\rm BH}/M_{\rm galaxy}$, decreases with a smaller bulge-to-disk ratio, due to the effects of geometrical dilution and opacity, and is reduced maximally by two orders of magnitude, resulting in $M_{\rm BH}/M_{\rm galaxy}\approx 10^{-5}$. In contrast, if only the bulge components in galaxies are focused, the BH-to-bulge mass ratio becomes $M_{\rm BH}/M_{\rm bulge}\approx 10^{-3}$, which is similar to that found in elliptical galaxies. Thus, it turns out that the mass of SMBH primarily correlates with a bulge, not with a disk, consistently with observational data.

Recently, ``CDM crisis'' is under discussion. The main point of this crisis is that number of substructures presented by cosmological N-body simulations based on CDM scenario for structure formation is much larger than observed substructures. Therefore, it is crucial for this crisis to discriminate whether expected number of CDM substructures really exist but non-luminous or do not exist. In this paper, we present a new idea to detect such invisible substructures by utilizing a gravitational lensing. Here, we consider quasars that are gravitationally lensed by a foreground galaxy. A substructure around the lensing galaxy may superposed on one of the lensed images of such quasars. In this situation, additional image splitting should occur in the image behind the substructure, and further multiple images are created. This is ``quasar mesolensing''. We estimate separation and time delay between further multiple images due to quasar mesolensing. The expected value is $1 \sim 30$ milli-arcsecond for the separation and future fine resolution imaging enable us to find invisible substructures, and is $1 \sim 10^3$ second for the time delay and high-speed monitoring of such quasar will be able to find ``echo''-like variation due to quasar mesolensing in intrinsic variability of the quasar. Furthermore, we evaluate that the optical depth for the quasar mesolensing is $\sim 0.1$. Consequently, if we monitor a few multiple quasars, we can find ``echo''-like variation in one of the images after intrinsic flux variations of quasars.

We reanalyze the photoevaporation problem of subgalactic objects irradiated by UVB radiation in a reionized universe. For the purpose, we perform three-dimensional radiation smoothed-particle-hydrodynamics calculations, where the radiative transfer is solved by a direct method and also the non-equilibrium chemistry of primordial gas including H$_2$ molecules is incorporated. Attention is concentrated on radiative transfer effects against UVB for the formation of subgalactic objects with $T_{\rm vir}\la 10^4$K. We consider the reionization model with $z_{\rm reion} \approx 7$ and also the earlier reionization model ($z_{\rm reion} \approx 17$) inferred by the WMAP. We find that the star formation is suppressed appreciably by UVB, but baryons at high-density peaks are self-shielded even during the reionization, forming some amount of stars eventually. In that sense, the photoevaporation for subgalactic systems is not so perfect as argued by one-dimensional spherical calculations. The final stellar fraction depends on the collapse epoch and the mass of system, but almost regardless of the reionization epoch. For instance, a few tenths of formed stars are born after the cosmic reionization in $z_{\rm reion} \approx 7$ model, while more than 90% stars are born after the reionization in the WMAP reionization model. Thus, effects of UVB feedback on the substructure problem with CDM scenario should be evaluated with careful treatment of the radiative transfer. The star clusters formed at high-density peaks coalesce with each other in a dissipationless fashion in a dark matter potential, resultantly forming a spheroidal system. As a result, these low-mass galaxies have large mass-to-light ratios such as observed in dwarf spheroidals in the Local Group.

A novel mechanism to build up a supermassive black hole (SMBH) is proposed. Here, relativistic radiation-hydrodynamic processes are considered from a galactic scale to a horizon scale. In the present scenario, the mass of a SMBH is predicted to be in proportion to the bulge mass, and it turns out that the SMBH-to-bulge mass ratio is basically determined by the nuclear energy conversion efficiency from hydrogen to helium, $\epsilon=0.007$. The BH hierarchy is also addressed in relation to the BH formation in globular clusters. Futhermore, a new picture for QSO formation is proposed based on the present SMBH formation scenario. It is predicted that a host luminosity-dominant ``proto-QSO phase'' exists before the QSO phase, and the proto-QSO phase is preceded by an optically-thick ultraluminous infrared galaxy (ULIRG) phase.

A new method which enables us to estimate the inclination of Class I young stellar objects is proposed. Since Class I objects are not spherically symmetric, it is likely that the observed feature is sensitive to the inclination of the system. Thus, we construct a protostar model by carefully treating two-dimensional (2D) radiative transfer and radiative equilibrium. We show from the present 2D numerical simulations that the emergent luminosity L_SED,which is the frequency integration of spectral energy distribution (SED), depends strongly on the inclination of the system i, whereas the peak flux is insensitive to i. Based on this result, we introduce a novel indicator f_L, which is the ratio of L_SED to the peak flux, as a good measure for the inclination. By using f_L, we can determine the inclination regardless of the other physical parameters. The inclination would be determined by f_L within the accuracy of +- 5 degree, if the opening angle of bipolar outflows is specified by any other procedure. Since this spectrophotometric method is easier than a geometrical method or a full SED fitting method, this method could be a powerful tool to investigate the feature of protostars statistically with observational data which will be provided by future missions, such as SIRTF, ASTRO-F, and ALMA.

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.

Based on the radiation hydrodynamical model for the black hole(BH) growth,incorporated with the chemical evolution of the early-type host galaxy, we construct the coevolution model of a QSO BH and the host galaxy. As a result, it is found that after a galactic wind epoch,the luminosity is shifted from the host-dominant phase to the AGN-dominant phase (QSO phase) in the timescale of a few $10^{8}$ years.The former phase corresponds to the early stage of growing BH, and can be regarded as a ``proto-QSO'' phase. It has observable characteristic properties (detail inthis paper).By comparing these predictions with recent observations, radio galaxies are a possible candidate for proto-QSOs.Also, it is anticipated that the proto-QSO phase is preceded by an optically thick phase, which may correspond to ULIRGs.In this phase, $M_{\rm BH}/M_{\rm bulge}$ is predicted to be much less than $10^{-3}$ and grow with metallicity.Moreover, as precursors of ULIRGs, optically-thin star-forming galaxies are predicted. These may be in the assembly phase of Lyman break galaxies (LBGs) or Ly$\alpha$ emitters.

We investigate the photon-trapping effects in the super-critical black hole accretion flows by solving radiation transfer as well as the energy equations of radiation and gas. It is found that the slim-disk model generally overestimates the luminosity of the disk at around the Eddington luminosity (L_E) and is not accurate in describing the effective temperature profile, since it neglects time delay between energy generation at deeper inside the disk and energy release at the surface. Especially, the photon-trapping effects are appreciable even below L ~ L_E, while they appear above ~ 3L_E according to the slim disk. Through the photon-trapping effects, the luminosity is reduced and the effective temperature profile becomes flatter than r^{-3/4} as in the standard disk. In the case that the viscous heating is effective only around the equatorial plane, the luminosity is kept around the Eddington luminosity even at very large mass accretion rate, Mdot>>L_E/c^2. The effective temperature profile is almost flat, and the maximum temperature decreases in accordance with rise in the mass accretion rate. Thus, the most luminous radius shifts to the outer region when Mdot/(L_E/c^2) >> 10^2. In the case that the energy is dissipated equally at any heights, the resultant luminosity is somewhat larger than in the former case, but the energy-conversion efficiency still decreases with increase of the mass accretion rate, as well. The most luminous radius stays around the inner edge of the disk in the latter case. Hence, the effective temperature profile is sensitive to the vertical distribution of energy production rates, so is the spectral shape. Future observations of high L/L_E objects will be able to test our model.

In the context of star formation through fragmentation of an extremely metal-deficient protogalactic cloud, the gravitational collapse of filamentary gas clouds is explored with 1D numerical hydrodynamics coupled with non-equilibrium chemistry of H$_2$ and HD. It is found that the cloud evolution is governed mainly by the initial central density ($n_{\rm c,0}$) and H$_2$ abundance ($x_{\rm H_2,0}$). In particular, the evolution of low-density filaments bifurcates at a threshold H$_2$ abundance of $x_{\rm H_2,cr}\simeq 3 \times 10^{-3}$, beyond which HD cooling overwhelms H$_2$ cooling. The contraction of a low density filament with $x_{\rm H_2, 0}\gtrsim x_{\rm H_2,cr}$ is strongly decelerated when $n_{\rm c}$ reaches a critical density of HD, and the filament is expected to fragment at $\sim 10^{7}$ cm$^{-3}$. The fragment mass is lowered to be $\approx 10M_\odot$. In contrast, the contraction of a low density filament with $x_{\rm H_2, 0}\lesssim x_{\rm H_2,cr}$ is regulated by H$_2$ cooling. In this case, the fragment mass is as high as $\approx 10^2M_\odot$. For a high-density filament, the cloud evolution is governed by H$_2$ cooling. The fragmentation is not expected to take place until the cloud becomes opaque in H$_2$ lines at $n_{\rm c,0}\sim 10^{12-13}$ cm$^{-3}$, so that the fragment mass is reduced to 1-2 M$_\odot$. As a result, the stellar IMF could be bimodal and deficient in sub-solar mass stars, where the high mass peak is around $10M_\odot$ or $10^2M_\odot$, dependently on $n_{\rm c,0}$ and $x_{\rm H_2,0}$.

We quantitatively scrutinize the effects of the radiation drag arising from the radiation fields in a galactic bulge in order to examine the possibility that the radiation drag could be an effective mechanism to extract angular momentum in a spheroidal system like a bulge and allow plenty of gas to accrete onto the galactic center. For this purpose, we numerically solve the relativistic radiation hydrodynamical equation coupled with the accurate radiative transfer and quantitatively assess the radiation drag efficiency. As a result, we find that in an optically thick regime the radiation drag efficiency is sensitively dependent on the density distributions of interstellar medium (ISM). The efficiency drops according to $\tau_{\rm T}^{-2}$ in an optically thick {\it uniform} ISM, where $\tau_{\rm T}$ is the total optical depth of the dusty ISM, whereas the efficiency remains almost constant at a high level if the ISM is {\it clumpy}. Hence, if the bulge formation begins with a star formation event in a clumpy ISM, the radiation drag will effectively work to remove the angular momentum and the accreted gas may form a supermassive black hole. As a natural consequence, this mechanism reproduces a putative linear relation between the mass of a supermassive black hole and the mass of a galactic bulge, although further detailed modeling for stellar evolution is required for the more precise prediction.

As a potential mechanism to build up supermassive black holes (BHs) in a spheroidal system, we consider the radiation drag effect by bulge stars, which extracts angular momentum from interstellar gas and thus allows the gas to accrete onto the galactic center. With incorporating radiation hydrodynamical equation with simple stellar evolution, it is shown that the BH-to-bulge mass ratio, $f_{BH}$, is basically determined by a fundamental constant, that is, the energy conversion efficiency for nuclear fusion of hydrogen to helium, $\epsilon=0.007$. More specifically, $f_{BH}$ is predicted to be $0.3\epsilon -0.5\epsilon$. Based on the present model for BH growth, a scenario for quasar formation is addressed in relation to ultraluminous infrared galaxies.

Using 3D radiative transfer calculations on the reionization of an inhomogeneous universe, QSO absorption line systems are simulated and they are compared with observations of Ly$\alpha$ continuum depression at high redshifts. By this comparison, it is found that the metagalactic UV intensity decreases rapidly with $z$ at $z>4$ as $I_{21}=0.5\exp [3(4-z)]$, and the reionization must have taken place between $z=6$ and 10. Based on this time-dependence of UV background intensity, we explore the collapse of pregalactic clouds in the UV background, and find that the self-shielding is prominent above a mass scale as $ M_{BIF}=3.0\times 10^{11} M_\odot [(1+z_c)/5]^{-4.2}(I_{21}/0.5)^{0.6}. $ This mass scale predicts the bifurcation of galactic morphology, and by confrontation with observations it turns out that the bifurcation mass successfully discriminates between elliptical and spiral galaxies.

Full radiative transfer in 3D space has been solved to pursue the reionization history in an inhomogeneous universe. It has been shown that the reionization of an inhomogeneous universe is not a prompt event, but a fairly slow process. Also, QSO absorption line systems are simulated with using the results of radiative transfer calculations. Ly$\alpha$ continuum depression implies that the metagalactic UV intensity decreases rapidly with $z$ at $z>5$, and the reionization must have taken place between $z=6$ and 10. Finally, it is stressed that H$\alpha$ forest is a more powerful tool to probe the reionization history and the density fluctuations in the universe at $z>5$.

Cosmological evolution of the QSO luminosity functions (LFs) at NIR/optical/X-ray bands for 1.3 < z < 3.5 is investigated based on the realistic QSO spectra. The accretion-disk theory predicts that although QSO luminosities only depend on mass-accretion rate, \Mdot, QSO spectra have a dependence on black-hole mass, M_{BH}, as well. The smaller M_{BH} is and/or the larger \Mdot is, the harder becomes the QSO NIR/optical/UV spectrum. We model disk spectra which can reproduce these features and calculated LFs for redshift z ~ 3 with the assumption of new-born QSOs being shining at the Eddington luminosity. The main results are: (i) the observed LFs at optical and X-rays can be simultaneously reproduced. (ii) LFs at optical and X-ray bands are not sensitive to M_{BH}, while LFs at NIR bands are; about one order of magnitude difference is expected in volume number densities at L_{I, J} ~ 10^{46} erg s^{-1} between the case that all QSOs would have the same spectral shape as that of M_{BH} = 10^{9} M_sun and the case with M_{BH} = 10^{11} M_sun. (iii) The resultant LFs at NIR are dominated by 10^{7} M_sun black-holes at L_{I, J} ~ 10^{44} erg s^{-1}, and by 10^{11} M_sun black-holes at L_{I, J} \~ 10^{46} erg s^{-1}. Future infrared observations from space(e.g.NGST) will probe cosmological evolution of black hole masses. For redshift z < 3, on the other hand, the observed optical/X-ray LFs can be fitted, if the initial QSO luminosity L_0 is below the Eddington luminosity. Interestingly, the best fitting values of l = L_0/L_{Edd} are different in B- and X-ray bands; l_B ~ 2.5 l_X. The reason for this discrepancy is briefly discussed.

In the context of the star formation through the fragmentation of an extremely metal-deficient protogalactic cloud, the gravitational collapse of filamentary gas clouds is explored with H$_2$ and HD chemistry. It is found by 1D hydrodynamical simulations that the cloud evolution is prescribed mainly by the initial density ($n_0$) and H$_2$ abundance ($x_{\rm H_2,0}$). In particular, it turns out that the evolution of low-density filaments ($n_0 \lesssim 10^5$ cm$^{-3}$) bifurcates at a critical H$_2$ abundance of $x_{\rm H_2,cr}\simeq 3\times 10^{-3}$, beyond which HD cooling overwhelms H$_2$ cooling. The numerical results indicate that the stellar IMF is likely to be double-peaked and deficient in sub-solar mass stars, where the high mass peak of the IMF is around $10M_\odot$ or $10^2M_\odot$, dependently on the initial density and H$_2$ abundance. If the gas in protogalactic clouds is photoionized by UV radiation or shock-heated, the H$_2$ abundance could exceed $x_{\rm H_2,cr}\simeq 3\times 10^{-3}$ by H$^-$ reactions. Then, the high mass peak would be $O(10) M_\odot$.

We explore the formation of dusty gas walls induced by a circumnuclear starburst around an AGN. We concentrate our attention on the role of the radiation force by a starburst as well as an AGN, where the effects of optical depth of dusty gas are taken into consideration. In two-dimensional axisymmetric space, we analyze the configuration and the stability of geometrically thin walls which are in balance between radiation pressure and gravity. As a result, it is shown that the radiation force by the circumnuclear starburst works to stabilize optically thick walls surrounding the nucleus. In the case of a brighter starburst with a fainter AGN, there form double walls, an inner one of which is located between the nucleus and the circumnuclear starburst, and an outer one of which enshrouds both the starburst regions and the nucleus. The total extinction (Av) of both walls turns out to be ~10 for a brighter starburst. As a consequence, double walls could heavily obscure the nucleus to make the AGN type 2. The outer wall may provide an explanation for the recent indications for large-scale obscuring materials in Seyfert 2's. Also, it is predicted that the AGN type is shifts from type 2 to type 1 in several times 10^7 yr according to the stellar evolution in the starburst. In contrast, if the AGN itself is much brighter than the starburst as a quasar is, then neither wall forms regardless of the starburst activity and the nucleus is likely to be identified as type 1. To conclude, the radiatively-supported gas walls could be responsible for the putative correlation between AGN type and the starbursts, whereby Seyfert 2 galaxies are more frequently associated with circumnuclear starbursts than type 1, whereas quasars are mostly observed as type 1 regardless of star-forming activity in the host galaxies.

We explore possibilities of collapse and star formation in Population III objects exposed to the external ultraviolet background (UVB) radiation. Assuming spherical symmetry, we solve self-consistently radiative transfer of photons, non-equilibrium H2 chemistry, and gas hydrodynamics. Although the UVB does suppress the formation of low mass objects, the negative feedback turns out to be weaker than previously suggested. In particular, the cut-off scale of collapse drops significantly below the virial temperature 10^4 K at weak UV intensities, due to both self-shielding of the gas and H2 cooling. Clouds above this cut-off tend to contract highly dynamically, further promoting self-shielding and H2 formation. For plausible radiation intensities and spectra, the collapsing gas can cool efficiently to temperatures well below 10^4 K before rotationally supported and the final H2 fraction reaches 10^{-3}. Our results imply that star formation can take place in low mass objects collapsing in the UVB. The threshold baryon mass for star formation is \sim 10^9 solar mass for clouds collapsing at redshifts z \simlt 3, but drops significantly at higher redshifts. In a conventional cold dark matter universe, the latter coincides roughly with that of the 1 \sigma density fluctuations. Objects near and above this threshold can thus constitute `building blocks' of luminous structures, and we discuss their links to dwarf spheroidal/elliptical galaxies and faint blue objects. These results suggest that the UVB can play a key role in regulating the star formation history of the Universe.

We consider the radiation-hydrodynamical formation of a dusty wall in the circumnuclear regions of the Circinus Seyfert 2 galaxy. We focus on the radiative flux force due to the circumnuclear starburst, the nuclear starburst and the active galactic nucleus (AGN), and analyze the equilibrium configuration and stability of the dusty gas in the circumnuclear regions. The mass distributions used to model the gravitational fields are determined by the observed rotation velocities of the stars and of the gas. Moreover, by using the simple stellar evolution in the circumnuclear starburst, the bolometric luminosity of the starburst is estimated as a function of time. As a result, it is shown that the radiation force by the circumnuclear starburst does play an important role to build up a radiatively-supported dusty wall of $\la$ 100 pc, which obscures the circumnuclear regions with $A_V$ of a few magnitudes. It is also found that the age of a circumnuclear starburst when the circumnuclear regions are enshrouded by the dusty wall is constrained to be $\la 10^8$ yr, which is consistent with the age observationally inferred for the starburst, $4 10^7$ yr $-1.5 10^8$ yr.

The collapse and fragmentation of filamentary primordial gas clouds are explored using 1D and 2D hydrodynamical simulations coupled with the nonequilibrium processes of H2 formation. The simulations show that depending upon the initial density,there are two occasions for the fragmentation of primordial filaments. If a filament has relatively low initial density, the radial contraction is slow due to less effective H2 cooling. This filament tends to fragment into dense clumps before the central density reaches $10^{8-9}$ cm$^{-3}$, where H2 cooling by three-body reactions is effective and the fragment mass is more massive than some tens $M_\odot$. In contrast, if a filament is initially dense, the more effective H2 cooling with the help of three-body reactions allows the filament to contract up to $n\sim 10^{12}$ cm$^{-3}$. After the density reaches $n\sim 10^{12}$ cm$^{-3}$, the filament becomes optically thick to H2 lines and the radial contraction subsequently almost stops. At this final hydrostatic stage, the fragment mass is lowered down to $\approx 1M_\odot$ because of the high density of the filament. The dependence of the fragment mass upon the initial density could be translated into the dependence on the local amplitude of random Gaussian density fields or the epoch of the collapse of a parent cloud. Hence, it is predicted that the initial mass function of Population III stars is likely to be bimodal with peaks of $\approx 10^2 M_\odot$ and $\approx 1M_\odot$, where the relative heights could be a function of the collapse epoch.

Based upon a novel paradigm of the galaxy formation under UV background, the evolutionary bifurcation of pregalactic clouds is confronted with observations on elliptical and spiral galaxies. The theory predicts that the dichotomy between the dissipational and dissipationless galaxy formation stems from the degree of self-shielding from the UV background radiation. This is demonstrated on a bifurcation diagram of collapse epochs versus masses of pregalactic clouds. Using the observed properties, the collapse epochs are assessed for each type of galaxies with attentive mass estimation. By the direct comparison of the theory with the observations, it turns out that the theoretical bifurcation branch successfully discriminates between elliptical and spiral galaxies. This suggests that the UV background radiation could play a profound role for the differentiation of the galactic morphology into the Hubble sequence.

To elucidate the effects of the UV background radiation on the collapse of pregalactic clouds, we implement a radiation-hydrodynamical calculation, combining one-dimensional spherical hydrodynamics with an accurate treatment of the radiative transfer of ionizing photons. Both absorption and scattering of UV photons are explicitly taken into account. It turns out that a gas cloud contracting within the dark matter potential does not settle into hydrostatic equilibrium, but undergoes run-away collapse even under the presence of the external UV field. The cloud center is shown to become self-shielded against ionizing photons by radiative transfer effects before shrinking to the rotation barrier. Based on our simulation results, we further discuss the possibility of H2 cooling and subsequent star formation in a run-away collapsing core. The present results are closely relevant to the survival of subgalactic Population III objects as well as to metal injection into intergalactic space.

The pancake collapse of pregalactic clouds under UV background radiation is explored with a one-dimensional sheet model. Here, attention is concentrated on elucidating the basic physics on the thermal evolution of pregalactic clouds exposed to diffuse UV radiation. So, we treat accurately the radiation transfer for the ionizing photons, with solving chemical reactions regarding hydrogen molecules as well as atoms. The self-shielding against UV radiation by H$_2$ Lyman-Werner bands, which regulates the photo-dissociation of hydrogen molecules, is also taken into account. As a result, it is found that when the UV background radiation is at a level of $10^{-22} (\nu/\nu_L)^{-1} {\rm erg~ s^{-1} cm^{-2} Hz^{-1} str^{-1}}$, the cloud evolution bifurcates with a critical mass as $M_{\rm SB} = 2.2\times 10^{11} M_\odot [(1+z_c)/5]^{-4.2}$, where $z_c$ is the final collapse epoch. A cloud more massive than $M_{\rm SB}$ cools below $5\times 10^3$K due to H$_2$ line emission at the pancake collapse and would undergo the initial starburst. The pancake possibly evolves into a virialized system in a dissipationless fashion. Consequently, this leads to the dissipationless galaxy formation at $3\lsim z_c \lsim 10$. A cloud less massive than $M_{\rm SB}$ cannot cool by H$_2$ emission shortly after the pancake collapse, but could cool in the course of shrinking to the rotation barrier. This is likely to lead to the dissipational galaxy formation at relatively low redshifts as $0\lsim z_c \lsim 4$. The present results provide a solid physical mechanism which controls the star formation efficiency in the pregalactic clouds. In the context of a standard CDM cosmology, $M_{\rm SB}$ lies between 1$\sigma$ and 2$\sigma$ density fluctuations.

By considering the radiative force by a circumnuclear starburst as well as an AGN, we analyze the equilibrium configuration and the stability of dusty gas in the circumnuclear regions. It is found that the radiative force by an intensive starburst can support a stable gaseous wall with a scale-height of several hundred parsecs. Moreover, by taking the simple stellar evolution in the starburst into account, we find that the covering factor of the wall decreases on a time-scale of several $10^7$ yr. The large-scale wall, if formed, works to obscure the nucleus due to the dust opacity. Hence, it is anticipated that the index of AGN type tends to shift from higher to lower in several $10^7$ yr according as the circumnuclear starburst becomes dimmer. On the other hand, if the AGN itself is brighter than the circumnuclear starburst (e.g. quasar case), no stable large-scale wall forms. In that case, the AGN is highly probably identified as type 1. The present mechanism may provide a physical explanation for the putative correlation between AGN type and host properties that Sy2's are more frequently associated with circumnuclear starbursts than Sy1's, whereas quasars are mostly observed as type 1 regardless of star-forming activity in the host galaxies.

Performing 1D hydrodynamical calculations coupled with non-equilibrium processes for H2 formation, we pursue the thermal and dynamical evolution of filamentary primordial clouds and attempt to make an estimate on the mass of population III stars. It is found that, almost independent of initial conditions, a filamentary cloud continues to collapse nearly isothermally due to H_2 cooling until the cloud becomes optically thick against the H_2 lines. During the collapse the cloud structure separates into two parts, i.e., a denser spindle and a diffuse envelope. The spindle contracts quasi-statically, and thus the line mass of the spindle keeps a characteristic value determined solely by the temperature ($\sim 800$ K). Applying a linear theory, we find that the spindle is unstable against fragmentation during the collapse. The wavelength of the fastest growing perturbation lessens as the collapse proceeds. Consequently, successive fragmentation could occur. When the central density exceeds $n_c \sim 10^{10-11} cm^{-3}$, the successive fragmentation may cease since the cloud becomes opaque against the H_2 lines and the collapse decelerates appreciably. The mass of the first star is then expected to be typically $\sim 3 M_\odot$, which may grow up to $\sim 16 M_\odot$ by accreting the diffuse envelope. Thus, the first-generation stars are anticipated to be massive but not supermassive.

The optical-ultraviolet continuum from active galactic nuclei (AGN) seem to originate from optically thick and/or thin disks, and occasionally from associated circumnuclear starburst regions. These different possible origins can, in principle, be discriminated by observations of gravitational microlensing event. We performed numerical simulations of microlensing of an AGN disk by a single lensing star passing in front of the AGN. Calculated spectral variations and light curves show distinct behavior, depending on the nature of emitting region; the time variation of a few months with strong wavelength dependence is expected in the case of an optically thick disk ({\it standard disk}), while an optically thin disk ({\it advection-dominated disk}) will produce shorter, nearly wavelength-independent variation. In the case of a starburst region much slower variations (over a year) will be superposed on the shorter variations.

To elucidate the permeation of cosmic ultraviolet (UV) background radiation into a pregalactic cloud and subsequent ionization, the frequency-dependent radiative transfer equation is solved, coupled with the ionization process. A spherical top-hat cloud which is composed of pure hydrogen is assumed. As a result, it is shown that the self-shielding, although it is often disregarded in cosmological hydrodynamic simulations, could start to emerge shortly after the maximum expansion stages of density fluctuations. The present analysis reveals that the Stro"mgren approximation leads to overestimation of the photoionization effects. Also, the self-shielded neutral core is no longer sharply separated from surrounding ionized regions; a low but noticeable degree of ionization is caused by high energy photons even in the self-shielded core. The present results may be substantial on considering the biasing by photoionization against low-mass galaxy formation.

We examine the radiatively driven mass accretion onto galactic nuclei due to the intensive radiation from circumnuclear starbursts. The radiation from a starburst not only contracts an inner gas disk by the radiation flux force, but also extracts angular momenta due to the relativistic radiation drag, thereby inducing an avalanche of the surface layer of disk. The present radiatively driven mass accretion may provide a physical mechanism which enables mass accretion from 100 pc scales down to pc, and it can eventually link to an advection-dominated viscous accretion onto a massive black hole. The radiation-hydrodynamical and self-gravitational instabilities of the disk are briefly discussed. In particular, the radiative acceleration possibly builds up a dusty wall, which provides another version for the formation of an obscuring torus.

The innermost regions of quasars can be resolved by a gravitational-lens {\lq}telescope{\rq} on scales down to a few AU. For the purpose, X-ray observations are most preferable, because X-rays originating from the innermost regions, can be selectively amplified by microlensing due to the so-called `caustic crossing'. If detected, X-ray variations will constrain the size of the X-ray emitting region down to a few AU. The maximum attainable resolution depends mainly on the monitoring intervals of lens events, which should be much shorter than the crossing time. On the basis of this idea, we performe numerical simulations of microlensing of an optically-thick, standard-type disk as well as an optically-thin, advection-dominated accretion flow (ADAF). Calculated spectral variations and light curves show distinct behaviors, depending on the photon energy. X-ray radiation which is produced in optically thin region, exhibits intensity variation over a few tens of days. In contrast, optical-UV fluxes, which are likely to come from optically thick region, exhibit more gradual light changes, which is consistent with the microlensing events so far observed in Q2237+0305. Currently, Q2237+0305 is being monitored in the optical range at Apache Point Observatory. Simultaneous multi-wavelength observations by X-ray sattelites (e.g., ASCA, AXAF, XMM) as well as HST at the moment of a microlens event enable us to reveal an AU scale structure of the central accretion disk around black hole.

Aperiodic optical variability is a common property of Active Galactic Nuclei (AGNs), though its physical origin is still open to question. To study the origin of the optical -- ultraviolet variability in AGN, we compare light curves of two models to observations of quasar 0957+561 in terms of a structure function analysis. In the starburst (SB) model, random superposition of supernovae in the nuclear starburst region produce aperiodic luminosity variations, while in the disk-instability (DI) model, variability is caused by instabilities in the accretion disk around a supermassive black hole. We calculate fluctuating light curves and structure functions, $V(\tau)$, by simple Monte-Carlo simulations on the basis of the two models. Each resultant $V(\tau)$ possesses a power-law portion, $[V(\tau)]^{1/2} \propto \tau^{\beta}$, at short time lags ($\tau$). The two models can be distinguished by the logarithmic slope, $\beta$; $\beta \sim$ 0.74--0.90 in the SB model and $\beta \sim$ 0.41--0.49 in the DI model, while the observed light curves exhibit $\beta \sim$ 0.35. Therefore, we conclude that the DI model is favored over the SB model to explain the slopes of the observational structure function, in the case of 0957+561, though this object is a radio-loud object and thus not really a fair test for the SB model. In addition, we examine the time-asymmetry of the light curves by calculating $V(\tau)$ separately for brightening and decaying phases. The two models exhibit opposite trends of time-asymmetry to some extent, although the present observation is not long enough to test this prediction.

A self-similar solution for time evolution of isothermal, self-gravitating viscous disks is found under the condition that $\alpha' \equiv \alpha (H/r)$ is constant in space (where $\alpha$ is the viscosity parameter and $H/r$ is the ratio of a half-thickness to radius of the disk). This solution describes a homologous collapse of a disk via self-gravity and viscosity. The disk structure and evolution is distinct in the inner and outer parts. There is a constant mass inflow in the outer portions so that the disk has flat rotation velocity, constant accretion velocity, and surface density decreasing outward as $\Sigma \propto r^{-1}$. In the inner portions, in contrast, mass is accumulated near the center owing to the boundary condition of no radial velocity at the origin, thereby a strong central concentration being produced; surface density varies as $\Sigma \propto r^{-5/3}$. Moreover, the transition radius separating the inner and outer portions increases linearly with time. The consequence of such a high condensation is briefly discussed in the context of formation of a quasar black hole.

The evolution of nonlinear density fluctuations around the Jeans mass shortly after cosmological recombination is analyzed using a 3D hydrodynamics/dark--matter code. The Cosmic Background Radiation (CBR) exerts Compton friction on free electrons due to peculiar velocities. The dynamics therefore depends strongly on the gas ionization history. Under a variety of ionization conditions and in systems with or without non-baryonic components, the baryons lose angular momentum efficiently and collapse to form a compact optically--thick object which would probably quickly evolve into a massive black hole. Attention is concentrated on elucidating some of the novel physical effects in early cosmological collapses, but ways in which more realistic calculations might be made and in which the scenario could be incorporated into a more complete cosmogonic model are discussed.