Existing models of galaxy formation have not yet explained striking correlations between structure and star-formation activity in galaxies, notably the sloped and moving boundaries that divide star-forming from quenched galaxies in key structural diagrams. This paper uses these and other relations to ``reverse-engineer'' the quenching process for central galaxies. The basic idea is that star-forming galaxies with larger radii (at a given stellar mass) have lower black-hole masses due to lower central densities. Galaxies cross into the green valley when the cumulative effective energy radiated by their black hole equals $\sim4\times$ their halo-gas binding energy. Since larger-radii galaxies have smaller black holes, one finds they must evolve to higher stellar masses in order to meet this halo-energy criterion, which explains the sloping boundaries. A possible cause of radii differences among star-forming galaxies is halo concentration. The evolutionary tracks of star-forming galaxies are nearly parallel to the green-valley boundaries, and it is mainly the sideways motions of these boundaries with cosmic time that cause galaxies to quench. BH-scaling laws for star-forming, quenched, and green-valley galaxies are different, and most BH mass growth takes place in the green valley. Implications include: the radii of star-forming galaxies are an important second parameter in shaping their black holes; black holes are connected to their halos but in different ways for star-forming, quenched, and green-valley galaxies; and the same BH-halo quenching mechanism has been in place since $z \sim 3$. We conclude with a discussion of black hole-galaxy co-evolution, the origin and interpretation of BH scaling laws.

We propose and analyze a class of robust, uniformly high-order accurate discontinuous Galerkin (DG) schemes for multidimensional relativistic magnetohydrodynamics (RMHD) on general meshes. A distinct feature of the schemes is their physical-constraint-preserving (PCP) property, i.e., they are proven to preserve the subluminal constraint on the fluid velocity and the positivity of density, pressure, and specific internal energy. Developing PCP high-order schemes for RMHD is highly desirable but remains a challenging task, especially in the multidimensional cases, due to the inherent strong nonlinearity in the constraints and the effect of the magnetic divergence-free condition. Inspired by some crucial observations at the PDE level, we construct the provably PCP schemes by using the locally divergence-free DG schemes of the recently proposed symmetrizable RMHD equations as the base schemes, a limiting technique to enforce the PCP property of the DG solutions, and the strong-stability-preserving methods for time discretization. We rigorously prove the PCP property by using a novel "quasi-linearization" approach to handle the highly nonlinear physical constraints, technical splitting to offset the influence of divergence error, and sophisticated estimates to analyze the beneficial effect of the additional source term in the symmetrizable RMHD system. Several two-dimensional numerical examples are provided to confirm the PCP property and to demonstrate the accuracy, effectiveness and robustness of the proposed PCP schemes.

This paper presents entropy symmetrization and high-order accurate entropy stable schemes for the relativistic magnetohydrodynamic (RMHD) equations. It is shown that the conservative RMHD equations are not symmetrizable and do not possess an entropy pair. To address this issue, a symmetrizable RMHD system, which admits a convex entropy pair, is proposed by adding a source term into the equations. Arbitrarily high-order accurate entropy stable finite difference schemes are then developed on Cartesian meshes based on the symmetrizable RMHD system. The crucial ingredients of these schemes include (i) affordable explicit entropy conservative fluxes which are technically derived through carefully selected parameter variables, (ii) a special high-order discretization of the source term in the symmetrizable RMHD system, and (iii) suitable high-order dissipative operators based on essentially non-oscillatory reconstruction to ensure the entropy stability. Several benchmark numerical tests demonstrate the accuracy and robustness of the proposed entropy stable schemes of the symmetrizable RMHD equations.

This paper proposes and analyzes arbitrarily high-order discontinuous Galerkin (DG) and finite volume methods which provably preserve the positivity of density and pressure for the ideal MHD on general meshes. Unified auxiliary theories are built for rigorously analyzing the positivity-preserving (PP) property of MHD schemes with a HLL type flux on polytopal meshes in any space dimension. The main challenges overcome here include establishing relation between the PP property and discrete divergence of magnetic field on general meshes, and estimating proper wave speeds in the HLL flux to ensure the PP property. In 1D case, we prove that the standard DG and finite volume methods with the proposed HLL flux are PP, under condition accessible by a PP limiter. For multidimensional conservative MHD system, standard DG methods with a PP limiter are not PP in general, due to the effect of unavoidable divergence-error. We construct provably PP high-order DG and finite volume schemes by proper discretization of symmetrizable MHD system, with two divergence-controlling techniques: locally divergence-free elements and a penalty term. The former leads to zero divergence within each cell, while the latter controls the divergence error across cell interfaces. Our analysis reveals that a coupling of them is important for positivity preservation, as they exactly contribute the discrete divergence-terms absent in standard DG schemes but crucial for ensuring the PP property. Numerical tests confirm the PP property and the effectiveness of proposed PP schemes. Unlike conservative MHD system, the exact smooth solutions of symmetrizable MHD system are proved to retain the positivity even if the divergence-free condition is not satisfied. Our analysis and findings further the understanding, at both discrete and continuous levels, of the relation between the PP property and the divergence-free constraint.

The density and pressure are positive physical quantities in magnetohydrodynamics (MHD). Design of provably positivity-preserving (PP) numerical schemes for ideal compressible MHD is highly desirable, but remains a challenge especially in the multidimensional cases. In this paper, we first develop uniformly high-order discontinuous Galerkin (DG) schemes which provably preserve the positivity of density and pressure for multidimensional ideal MHD. The schemes are constructed by using the locally divergence-free DG schemes for the symmetrizable ideal MHD equations as the base schemes, a PP limiter to enforce the positivity of the DG solutions, and the strong stability preserving methods for time discretization. The significant innovation is that we discover and rigorously prove the PP property of the proposed DG schemes by using a novel equivalent form of the admissible state set and very technical estimates. Several two-dimensional numerical examples further confirm the PP property, and demonstrate the accuracy, effectiveness and robustness of the proposed PP methods.

As the second paper of a series on studying galaxy-galaxy lensing signals using the Sloan Digital Sky Survey Data Release 7 (SDSS DR7), we present our measurement and modelling of the lensing signals around groups of galaxies. We divide the groups into four halo mass bins, and measure the signals around four different halo-center tracers: brightest central galaxy (BCG), luminosity-weighted center, number-weighted center and X-ray peak position. For X-ray and SDSS DR7 cross identified groups, we further split the groups into low and high X-ray emission subsamples, both of which are assigned with two halo-center tracers, BCGs and X-ray peak positions. The galaxy-galaxy lensing signals show that BCGs, among the four candidates, are the best halo-center tracers. We model the lensing signals using a combination of four contributions: off-centered NFW host halo profile, sub-halo contribution, stellar contribution, and projected 2-halo term. We sample the posterior of 5 parameters i.e., halo mass, concentration, off-centering distance, sub halo mass, and fraction of subhalos via a MCMC package using the galaxy-galaxy lensing signals. After taking into account the sampling effects (e.g. Eddington bias), we found the best fit halo masses obtained from lensing signals are quite consistent with those obtained in the group catalog based on an abundance matching method, except in the lowest mass bin. Subject headings: (cosmology:) gravitational lensing, galaxies: clusters: general

After thirty years of the discovery of the fundamental plane, explanations to the tilt of the fundamental plane with respect to the virial plane still suffer from the need of fine-tuning. In this paper, we try to explore the origin of this tilt from the perspective of modified Newtonian dynamics (MOND) by applying the 16 Coma galaxies available in Thomas et al.[1]. Based on the mass models that can reproduce de Vaucouleurs' law closely, we find that the tilt of the traditional fundamental plane is naturally explained by the simple form of the MONDian interpolating function, if we assume a well motivated choice of anisotropic velocity distribution, and adopt the Kroupa or Salpeter stellar mass-to-light ratio. Our analysis does not necessarily rule out a varying stellar mass-to-light ratio.

We give the three-dimensional dynamical autonomous systems for most of the popular scalar field dark energy models including (phantom) quintessence, (phantom) tachyon, k-essence and general non-canonical scalar field models, change the dynamical variables from variables $(x, y, \lambda)$ to observable related variables $(w_{\phi}, \Omega_{\phi}, \lambda)$, and show the intimate relationships between those scalar fields that the three-dimensional system of k-essence can reduce to (phantom) tachyon, general non-canonical scalar field can reduce to (phantom) quintessence and k-essence can also reduce to (phantom) quintessence for some special cases. For the applications of the three-dimensional dynamical systems, we investigate several special cases and give the exactly dynamical solutions in detail. In the end of this paper, we argue that, it is more convenient and also has more physical meaning to express the differential equations of dynamical systems in $(w_{\phi}, \Omega_{\phi}, \lambda)$ instead of variables $(x, y, \lambda)$ and to investigate the dynamical system in 3-Dimension instead of 2-Dimension. We also raise a question about the possibility of the chaotic behavior in the spatially flat single scalar field FRW cosmological models in the presence of ordinary matter.

As the first paper in a series on the study of the galaxy-galaxy lensing from Sloan Digital Sky Survey Data Release 7 (SDSS DR7), we present our image processing pipeline that corrects the systematics primarily introduced by the Point Spread Function (PSF). Using this pipeline, we processed SDSS DR7 imaging data in $r$ band and generated a background galaxy catalog containing the shape information of each galaxy. Based on our own shape measurements of the galaxy images from SDSS DR7, we extract the galaxy-galaxy (GG) lensing signals around foreground spectroscopic galaxies binned in different luminosity and stellar mass. The overall signals are in good agreement with those obtained by \citet{Mandelbaum2005, Mandelbaum2006} from the SDSS DR4. The results in this paper with higher signal to noise ratio is due to the larger survey area than SDSS DR4, confirm that more luminous/massive galaxies bear stronger GG lensing signal. We also divide the foreground galaxies into red/blue and star forming/quenched subsamples and measured their GG lensing signals, respectively. We find that, at a specific stellar mass/luminosity, the red/quenched galaxies have relatively stronger GG lensing signals than their counterparts especially at large radii. These GG lensing signals can be used to probe the galaxy-halo mass relations and their environmental dependences in the halo occupation or conditional luminosity function framework.

$f(R)$ gravity is one of the simplest theories of modified gravity to explain the accelerated cosmic expansion. Although it is usually assumed that the quasi-Newtonian approach (a combination of the quasi-static approximation and sub-Hubble limit) for cosmic perturbations is good enough to describe the evolution of large scale structure in $f(R)$ models, some studies have suggested that this method is not valid for all $f(R)$ models. Here, we show that in the matter-dominated era, the pressure and shear equations alone, which can be recast into four first-order equations to solve for cosmological perturbations exactly, are sufficient to solve for the Newtonian potential, $\Psi$, and the curvature potential, $\Phi$. Based on these two equations, we are able to clarify how the exact linear perturbations fit into different limits. We find that the Compton length controls the quasi-static behaviours in $f(R)$ gravity. In addition, regardless the validity of quasi-static approximation, a strong version of the sub-Hubble limit alone is sufficient to reduce the exact linear perturbations in any viable $f(R)$ gravity to second order. Our findings disagree with some previous studies where we find little difference between our exact and quasi-Newtonian solutions even up to $k=10 c^{-1} \mathcal{H}_0$.

We give a full investigation on the dynamics of power-law kinetic quintessence $L(X, \phi)=V(\phi)(-X+X^2)$ by considering the potential related parameter $\Gamma$($=\frac{V V''}{V'^2}$) as a function of another potential parameter $\lambda$($=\frac{V'}{\kappa V^{3/2}}$), which correspondingly extends the analysis of the dynamical system of our universe from two-dimension to three-dimension. Beside the critical points found in previous papers, we find a new de-Sitter-like dominant attractor(\textit{cp$6$}) and give its stable condition using the center manifold theorem. For the dark energy dominant solution(\textit{cp$6$} and \textit{cp$7$}), it could be distinguished from canonical quintessence and tachyon models since the sound speed $c_s^2=0$ or $c_s^2\ll 1$. For the scaling solution (\textit{cp$8$}), it is very interesting that the sound speed $c_s^2=1/5$ while it behaves as ordinary matter. We therefore point out that the power-law kinetic quintessence should have different signatures on cold dark matter power spectrum and cosmic microwave background both at early time when this scalar field is an early dark energy with $\Omega_\phi$ being non-negligible at high redshift and at late time when it drives the accelerating expansion. We even do not know whether there are any degeneracies of the impacts between these two epoches. They are expected to be investigated in future.

We use a suite of cosmological hydrodynamic simulations, run by two fixed grid codes, to investigate the properties of solenoidal and dilatational motions of the intergalactic medium (IGM), and the impact of numerical viscosity on turbulence in a LCDM universe. The codes differ only in the spatial difference discretization. We find that (1) The vortical motion grows rapidly since $z=2$, and reaches $\sim 10 km/s -90 km/s$ at $z=0$. Meanwhile, the small-scale compressive ratio $r_{CS}$ drops from 0.84 to 0.47, indicating comparable vortical and compressive motions at present. (2) Power spectra of the solenoidal velocity possess two regimes, $\propto k^{-0.89}$ and $\propto k^{-2.02}$, while the total and dilatational velocity follow the scaling $k^{-1.88}$ and $k^{-2.20}$ respectively in the turbulent range. The IGM turbulence may contain two distinct phases, the supersonic and post-supersonic phases. (3) The non-thermal pressure support, measured by the vortical kinetic energy, is comparable with the thermal pressure for $\rho_b \simeq 10-100$, or $T <10^{5.5} K$ at $z=0.0$. The deviation of the baryon fraction from the cosmic mean shows a preliminary positive correlation with the turbulence pressure support. (4) A relatively higher numerical viscosity would dissipate both the compressive and vortical motions of the IGM into thermal energy more effectively, resulting in less developed vorticity, remarkably shortened inertial range, and leading to non-negligible uncertainty in the thermal history of gas accretion. Shocks in regions outside of clusters are significantly suppressed by numerical viscosity since $z=2$, which may directly cause the different levels of turbulence between two codes.

After carefully cross-identifying a previously discovered GALEX-selected Lyman Break Galaxy (LBG) candidates one-to-one with their optical counterparts in the field of the CDF-S, we re-estimate their photometric redshifts using multi-wavelength data from UV, optical to NIR. We refine a new updated sample of 383 LBGs at $0.7\la z \la 1.4$. Most LBGs are classified as starburst and irregular types. Ages spread from several Myr to 1.5Gyr. Their dust-corrected star formation rates (SFRs) and stellar masses ($M_*$) are from $4\my$ to $220\my$ and from $2.3\times 10^8 \msun$ to $4 \times 10^{11} \msun$. The rest-frame FUV luminosity function of LBGs are presented. LBGs of irregular types mainly distribute along the "main sequence" of star forming galaxies while most LBGs of starburst types locate in the starburst region. A "downsizing" effect is clearly found and LBGs distribute in the "blue" cloud. HST images in F606W ($V$ band) and F850LP ($z$ band) are taken from the GEMS and GOODS-S surveys. SExtractor and GALFIT are applied to get their morphological parameters. A morphological sample of 142 LBGs with reliable results of $\sersic$ and sizes in both bands is defined. We find that LBGs at $z\sim 1 $ are dominated by disk-like galaxies. Correlations between photometric and morphological properties of LBGs are investigated. Strong correlations between their half-light radii and $M_*$, i.e., size-stellar mass relations, are found in both bands. Physical connections between correlations and the "downsizing" effect are discussed.

We present a multi-wavelength study of a 3.6 $\mu$m-selected galaxy sample in the Extended Groth strip. The sample is complete for galaxies with stellar mass $>10^{9.5}$ \Msun and redshift $0.4<z<1.2$. In this redshift range, the IRAC 3.6 $\mu$m band measures the rest-frame near-infrared band, permitting nearly unbiased selection with respect to both quiescent and star-forming galaxies. The numerous spectroscopic redshifts available in the EGS are used to train an Artificial Neural Network to estimate photometric redshifts. The distribution of photometric redshift errors is Gaussian with standard deviation ${\sim}0.025(1+z)$, and the fraction of redshift failures (${>}3\sigma$ errors) is about 3.5%. A new method of validation based on pair statistics confirms the estimate of standard deviation even for galaxies lacking spectroscopic redshifts. Basic galaxy properties measured include rest-frame $U-B$ colors, $B$- and $K$-band absolute magnitudes, and stellar masses. We divide the sample into quiescent and star-forming galaxies according to their rest-frame $U-B$ colors and 24 to 3.6 \micron\ flux density ratios and derive rest $K$-band luminosity functions and stellar mass functions for quiescent, star forming, and all galaxies. The results show that massive, quiescent galaxies were in place by $z\approx1$, but lower mass galaxies generally ceased their star formation at later epochs.

We investigate the effects of dust on Ly{\alpha} photons emergent from an optically thick medium by solving the integro-differential equation of the radiative transfer of resonant photons. To solve the differential equations numerically we use the Weighted Essentially Non-oscillatory method (WENO). Although the effects of dust on radia-tive transfer is well known, the resonant scattering of Ly{\alpha} photons makes the problem non-trivial. For instance, if the medium has the optical depth of dust absorption and scattering to be {\tau}a>>1, {\tau}>>1, and {\tau}>>{\tau}a, the effective absorption optical depth in a random walk scenario would be equal to \surd {\tau}a({\tau}a+{\tau}). We show, however, that for a resonant scattering at frequency {\nu}0, the effective absorption optical depth would be even larger than {\tau}({\nu}0). If the cross section of dust scattering and absorption is frequency-independent, the double-peaked structure of the frequency profile given by the resonant scattering is basically dust-independent. That is, dust causes neither narrowing nor widening of the width of the double peaked profile. One more result is that the time scales of the Ly{\alpha} photon transfer in the optically thick halo are also basically independent of the dust scattering, even when the scattering is anisotropic. This is because those time scales are mainly determined by the transfer in the frequency space, while dust scattering, either isotropic or anisotropic, does not affect the behavior of the transfer in the frequency space when the cross section of scattering is wavelength-independent. This result does not support the speculation that dust will lead to the smoothing of the brightness distribution of Ly{\alpha} photon source with optical thick halo.

We introduce a simple model to explore the star formation histories of disk galaxies. We assume that the disk origins and grows by continuous gas infall. The gas infall rate is parametrized by the Gaussian formula with one free parameter: infall-peak time $t_p$. The Kennicutt star formation law is adopted to describe how much cold gas turns into stars. The gas outflow process is also considered in our model. We find that, at given galactic stellar mass $M_*$, model adopting late infall-peak time $t_p$ results in blue colors, low metallicity, high specific star formation rate and high gas fraction, while gas outflow rate mainly influences the gas-phase metallicity and star formation efficiency mainly influences the gas fraction. Motivated by the local observed scaling relations, we construct a mass-dependent model by assuming low mass galaxy has later infall-peak time $t_p$ and larger gas outflow rate than massive systems. It is shown that this model can be in agreement with not only the local observations, but also the observed correlations between specific star formation rate and galactic stellar mass $SFR/M_* \sim M_*$ at intermediate redshift $z<1$. Comparison between the Gaussian-infall model and exponential-infall model is also presented. It shows that the exponential-infall model predicts higher star formation rate at early stage and lower star formation rate later than that of Gaussian-infall. Our results suggest that the Gaussian infall rate may be more reasonable to describe the gas cooling process than the exponential infall rate, especially for low-mass systems.

With a state-of-the-art numerical method for solving the integral-differential equation of radiative transfer, we investigate the flux of the Ly$\alpha$ photon $\nu_0$ emergent from an optically thick halo containing a central light source. Our focus is on the time-dependent effects of the resonant scattering. We first show that the frequency distribution of photons in the halo are quickly approaching to a locally thermalized state around the resonant frequency, even when the mean intensity of the radiation is highly time-dependent. Since initial conditions are forgotten during the thermalization, some features of the flux, such as the two peak structure of its profile, actually are independent of the intrinsic width and time behavior of the central source, if the emergent photons are mainly from photons in the thermalized state. In this case, the difference $|\nu_{\pm}-\nu_0|$, where $\nu_{\pm}$ are the frequencies of the two peaks of the flux, cannot be less than $2$ times of Doppler broadening. We then study the radiative transfer in the case where the light emitted from the central source is a flash. We calculate the light curves of the flux from the halo. It shows that the flux is still a flash. The time duration of the flash for the flux, however, is independent of the original time duration of the light source but depends on the optical depth of the halo. Therefore, the spatial transfer of resonant photons is a diffusion process, even though it is not a purely Brownian diffusion. This property enables an optically thick halo to trap and store thermalized photons around $\nu_0$ for a long time after the cease of the central source emission. The photons trapped in the halo can yield delayed emission, of which the profile also shows typical two peak structure as that from locally thermalized photons. Possible applications of these results are addressed.

The Tully-Fisher relation of spiral galaxies shows notable dependence on morphological types, with earlier type spirals having systematically lower luminosity at fixed maximum rotation velocity $V_{max}$. This decrement of luminosity is more significant in shorter wavelengths. By modeling the rotation curve and stellar population of different morphological type spiral galaxies in combination, we find the $V_{max}$ of spiral galaxies is weakly dependent on the morphological type, whereas the difference of the stellar population originating from the bulge disk composition effect mainly account for the morphological type dependence of the Tully-Fisher relation.

The 21 cm emission and absorption from gaseous halos around the first generation of star depend on the Wouthuysen-Field (W-F) coupling, which relates the spin temperature with the kinetic temperature of hydrogen gas via the resonant scattering between Lyman alpha photons and neutral hydrogen. Although the center object generally is a strong source of these photons, the transfer of these photons in the 21 cm region is inefficient, as the optical depth of the photons is large. Consequently, these photons from the source may not be able to transfer to the entire 21 cm region timely to provide the W-F coupling. This problem is important because the lifetime of first stars generally is short. The problem is investigated with numerical solution of the integro-differential equation, which describes the kinetics of these resonant photons in both physical and frequency spaces. We show that the photon transfer process in the physical space is actually coupled to that in the frequency space. Firstly diffusion in the frequency space provides a shortcut for the diffusion in the physical space. It makes the mean time for the escape of the resonant photon in optical depth \tau media roughly proportional to the optical depth \tau, not \tau^2. Secondly the resonant scattering is effective in bouncing photons with a frequency which is not equal to initial frequency back to the initial frequency. This process can restore initial frequency photons and establish the local Boltzmann distribution of the photon spectrum around the initial frequency. Therefore, the mechanism of 'escape via shortcut' plus 'bounce back' enables W-F coupling to be properly realized in the 21 cm region around first stars. This mechanism also works for photons injected into the 21 cm region by redshift.

We develop a numerical solver for the integral-differential equations, which describes the radiative transfer of photon distribution in the frequency space with resonant scattering of Lyalpha photons by hydrogen gas in the early universe. The time-dependent solutions of this equation is crucial to the estimation of the effect of the Wouthuysen-Field (WF) coupling in relation to the 21 cm emission and absorption at the epoch of reionization. The resonant scattering leads to the photon distribution in the frequency space to be piecewise smooth containing sharp changes. The weighted essentially nonoscillatory (WENO) scheme is suitable to handle this problem, as this algorithm has been found to be highly stable and robust for solving Boltzmann equation. We test this numerical solver against analytic solutions of the evolution of the photon distribution in rest background, analytic solution in expanding background without resonant scattering and formation of local Boltzmann distribution around the resonant frequency with the temperature same as that of atom for recoil. We find that evolution of photon distribution undergoes three phases; profile is similar to the initial one, a flat plateau (without recoil) or local Boltzmann distribution (with recoil) forms around the resonant frequency, and finally the distribution around the resonant frequency is saturated when the photons from the source is balanced by the redshift of the expansion. This result indicates that the onset of the W-F coupling should not be determined by the third phase, but by the time scale of the second phase. We found that the time scale of the W-F coupling is equal to about a few hundreds of the mean free flight time of photons with resonant frequency, and is independent of the Sobolev parameter if this parameter is much less than 1.

We study the Wouthuysen-Field coupling at early universe with numerical solutions of the integrodifferential equation describing the kinetics of photons undergoing resonant scattering. The numerical solver is developed based on the weighted essentially non-oscillatory (WENO) scheme for the Boltzmann-like integrodifferential equation. We focus on the time evolution of the Wouthuysen-Field (W-F) coupling in relation to the 21 cm emission and absorption at the epoch of reionization. We show that a local Boltzmann distribution will be formed if photons with frequency \sim \nu_0 have undergone a ten thousand or more times of scattering, which corresponds to the order of 10^3 yrs for neutral hydrogen density of the concordance \Lambda CDM model. The time evolution of the shape and width of the local Boltzmann distribution actually doesn't dependent on the details of atomic recoil, photon sources, or initial conditions very much. However, the intensity of photon flux at the local Boltzmann distribution is substantially time-dependent. The time scale of approaching the saturated intensity can be as long as 10^5-10^6 yrs for typical parameters of the \Lambda CDM model. The intensity of the local Boltzmann distribution at time less than 10^5 yrs is significantly lower than that of the saturation state. Therefore, it may not be always reasonable to assume that the deviation of the spin temperature of 21 cm energy states from cosmic background temperature is mainly due to the W-F coupling if first stars or their emission/absorption regions evolved with a time scale equal to or less than Myrs.

While both isothermal and NFW-based mass models for galaxy clusters are widely adopted in strong-lensing studies, they cannot easily be distinguished based solely on observed positions of arcs and arclets. We compare the magnifications predicted for giant arcs obtained from isothermal and NFW profiles, taking axially-symmetric and asymmetric mass distributions into account. We find that arc magnifications can differ strongly between the two types of density profiles even if the image morphology is well reproduced. Magnifications by lenses with NFW density profiles are usually larger than those for lenses with singular or nearly singular isothermal density profiles, unless the latter have large cores. Asymmetries play an important role. We illustrate our results with the two well-studied clusters MS 2137 and A~370. We confirm earlier results showing that both isothermal and NFW mass models can very well reproduce the observed arcs, radial arcs and other arclets. While the mass model for MS 2137 is not very well constrained, the two types of mass models produce strongly differing critical curves and caustics for A~370. We find that the NFW mass profile is preferred for A~370. We identify new arclet candidates in the field of A~370. Redshift estimates allowed by the lens model are consistent with results in the literature, where available. Three newly found counter-images are suggested to arise from an active, dust-enshrouded star-forming galaxy at $z\approx1.1$.

The strong lensing modelling of gravitational ``rings'' formed around massive galaxies is sensitive to the amplitude of the external shear and convergence produced by nearby mass condensations. In current wide field surveys, it is now possible to find out a large number of rings, typically 10 gravitational rings per square degree. We propose here, to systematically study gravitational rings around galaxy clusters to probe the cluster mass profile beyond the cluster strong lensing regions. For cluster of galaxies with multiple arc systems, we show that rings found at various distances from the cluster centre can improve the modelling by constraining the slope of the cluster mass profile. We outline the principle of the method with simple numerical simulations and we apply it to 3 rings discovered recently in Abell~1689. In particular, the lens modelling of the 3 rings confirms that the cluster is bimodal, and favours a slope of the mass profile steeper than isothermal at a cluster radius $\sim 300 \kpc$. These results are compared with previous lens modelling of Abell~1689 including weak lensing analysis. Because of the difficulty arising from the complex mass distribution in Abell~1689, we argue that the ring method will be better implemented on simpler and relaxed clusters.

We investigate the volume growth of ionized regions around UV photon sources with the WENO algorithm, which is an effective solver of photon kinetics in the phase space described by the radiative transfer equation. We show that the volume growth rate, either of isolated ionized regions or of clustered regions in merging, generally consists of three phases: fast or relativistic growth phase at the early stage, slow growth phase at the later stage, and a transition phase between the fast and slow phases. We also show that the volume growth of ionized regions around clustered sources with intensity $\dot{E}_i$ ($i=1, 2, ...$) would have the same behavior as a single source with intensity $\dot{E}=\sum_i\dot{E}_i$, if all the distances between nearest neighbor sources $i$ and $j$ are smaller than $c(t^i_c+t^j_c)$, $t^i_c$ being the time scale $t_c$ of source $i$. Therefore, a tightly clustered UV photon sources would lead to a slow growth of ionized volume. This effect would be important for studying the redshift-dependence of 21cm signals from the reionization epoch.

Based on the catalog of merging galaxies in the Early Data Release (EDR) of the Sloan Digital Sky Survey (SDSS), the differences of star formation history between merging galaxies and field galaxies are studied statistically by means of three spectroscopic indicators the 4000-\r{A} break strength, the Balmer absorption-line index, and the specific star formation rate. It is found that for early-type merging galaxies the interactions will not induce significant enhancement of the star-formation activity because of its stability and lack of cool gas. On the other hand, late-type merging galaxies always in general display more active star formation than field galaxies on different timescales within about 1Gyr. We also conclude that the mean stellar ages of late-type merging galaxies are younger than those of late-type field galaxies.

We investigate the 21 cm signals from the UV ionizing sources in the reionization epoch. The formation and evolution of 21 cm emission and absorption regions depend essentially on the kinetics of photons in the physical and frequency spaces. To solve the radiative transfer equation, we use the WENO algorithm, which is effective to capture the sharp ionization profile and the cut-off at the front of light $(r=ct)$ and to handle the small fraction of neutral hydrogen and helium in the ionized sphere. We show that a spherical shell of 21 cm emission and absorption will develop around a point source once the speed of the ionization front (I-front) is significantly lower than the speed of light. The 21 cm shell extends from the I-front to the front of light; its inner part is the emission region and its outer part is the absorption region. The 21 cm emission region depends strongly on the intensity, frequency-spectrum and life-time of the UV ionizing source. For a source of short life-time, no 21 cm emission region can be formed if the source dies out before the I-front speed is significantly lower than the speed of light. Yet, a 21 cm absorption region can form and develop even after the emission of the source ceases.

We develop a numerical solver for radiative transfer problems based on the weighted essentially nonoscillatory (WENO) scheme modified with anti-diffusive flux corrections, in order to solve the temperature and ionization profiles around a point source of photons in the reionization epoch. Algorithms for such simulation must be able to handle the following two features: 1. the sharp profiles of ionization and temperature at the ionizing front (I-front) and the heating front (T-front), and 2. the fraction of neutral hydrogen within the ionized sphere is extremely small due to the stiffness of the rate equations of atom processes. The WENO scheme can properly handle these two features, as it has been shown to have high order of accuracy and good convergence in capturing discontinuities and complicated structures in fluid as well as to be significantly superior over piecewise smooth solutions containing discontinuities. With this algorithm, we show the time-dependence of the preheated shell around a UV photon source. In the first stage the I-front and T-front are coincident, and propagate with almost the speed of light. In later stage, when the frequency spectrum of UV photons is hardened, the speeds of propagation of the ionizing and heating fronts are both significantly less than the speed of light, and the heating front is always beyond the ionizing front. In the spherical shell between the I- and T-fronts, the IGM is heated, while atoms keep almost neutral. The time scale of the preheated shell evolution is dependent on the intensity of the photon source. We also find that the details of the pre-heated shell and the distribution of neutral hydrogen remained in the ionized sphere are actually sensitive to the parameters used. The WENO algorithm can provide stable and robust solutions to study these details.

We report the detection of rest--frame 6.2 and 7.7 \micron emission features arising from Polycyclic Aromatic Hydrocarbons (PAH) in the Spitzer/IRS spectrum of an infrared-luminous Lyman break galaxy at z=3.01. This is currently the highest redshift galaxy where these PAH emission features have been detected. The total infrared luminosity inferred from the MIPS 24 \micron and radio flux density is 2$\times10^{13}$ L$_{\odot}$, which qualifies this object as a so--called hyperluminous infrared galaxy (HyLIRG). However, unlike local HyLIRGs which are generally associated with QSO/AGNs and have weak or absent PAH emission features, this HyLIRG has very strong 6.2 and 7.7 \micron PAH emission. We argue that intense star formation dominates the infrared emission of this source, although we cannot rule out the presence of a deeply obscured AGN. This LBG appears to be a distorted system in the HST ACS F606W and F814W images, possibly indicating that a significant merger or interaction is driving the large IR luminosity.

We use a large sample of galaxies drawn from the Sloan Digital Sky Survey (SDSS) and Two Micro All Sky Survey (2MASS) to present Color-Magnitude Relations (CMRs) for late-type galaxies in both optical and optical-infrared bands. A sample from SDSS Data Release 4 (DR4) is selected to investigate the optical properties. Optical-infrared colors are estimated from a position matched sample of DR4 and 2MASS, in which the photometric aperture mismatch between these two surveys is carefully corrected. It is shown that, after correcting the dust attenuation, the optical colors for faint galaxies (i.e. $M_r > -21$) have very weak correlation with the luminosity, while the optical colors for bright galaxies (i.e. $M_r < -21$) are redder for more luminous galaxies. All (optical, optical-infrared and infrared) colors show similar but stronger correlations with stellar mass than with absolute magnitude. The optical colors correlate more strongly with stellar mass surface density than with stellar mass, while optical-infrared and infrared colors show stronger correlations with stellar mass. By comparing the observed colors of our sample galaxies with the colors predicted by stellar population synthesis model, we find that massive late-type galaxies have older and higher metallicity stellar population than less massive galaxies. This suggests that CMRs for late-type galaxies are trends defined by the combination of stellar mean age and metallicity. Moreover, our results suggest that the stellar mean metallicity of late-type galaxy is mainly determined by its stellar mass, while the star formation history is mainly regulated by the stellar mass surface density.

We describe the spectral energy distributions (SEDs) of Lyman Break Galaxies (LBGs) at z~3 using deep mid-infrared and optical observations of the Extended Groth Strip, obtained with IRAC and MIPS on board Spitzer and from the ground, respectively. We focus on LBGs with detections at all four IRAC bands, in particular the 26 galaxies with IRAC 8 micron band (rest--frame K-band) detections. We use stellar population synthesis models and probe the stellar content of these galaxies. Based on best--fit continuous star-formation models we derive estimates of the stellar mass for these LBGs. As in previous studies, we find that a fraction of LBGs have very red colors and large estimated stellar masses (M_stellar > 5 x 10 ^{10} M_solar): the present Spitzer data allow us, for the first time, to study these massive LBGs in detail. We discuss the link between these LBGs and submm-luminous galaxies.We find that the number density of these massive LBGs at high redshift is higher than predicted by current semi-analytic models of galaxy evolution.

We show that the algorithm based on the weighted essentially nonoscillatory (WENO) scheme with anti-diffusive flux corrections can be used as a solver of the radiative transfer equations. This algorithm is highly stable and robust for solving problems with both discontinuities and smooth solution structures. We test this code with the ionized sphere around point sources. It shows that the WENO scheme can reveal the discontinuity of the radiative or ionizing fronts as well as the evolution of photon frequency spectrum with high accuracy on coarse meshes and for a very wide parameter space. This method would be useful to study the details of the ionized patch given by individual source in the epoch of reionization. We demonstrate this method by calculating the evolution of the ionized sphere around point sources in physical and frequency spaces. It shows that the profile of the fraction of neutral hydrogen and the ionized radius are sensitively dependent on the intensity of the source.

We investigate the statistical properties of cosmic baryon fluid in the nonlinear regime, which is crucial for understanding the large-scale structure formation of the universe. With the hydrodynamic simulation sample of the Universe in the cold dark matter model with a cosmological constant, we show that the intermittency of the velocity field of cosmic baryon fluid at redshift z=0 in the scale range from the Jeans length to about 16 Mpc/h can be extremely well described by She-Leveque's universal scaling formula. The baryon fluid also possesses the following features: (1) for volume weight statistics, the dissipative structures are dominated by sheets, and (2) the relation between the intensities of fluctuations is hierarchical. These results imply that the evolution of highly evolved cosmic baryon fluid is similar to a fully developed turbulence.

A deep mid- and far-infrared survey in the Extended Groth Strip (EGS) area gives 3.6 to 8\micron flux densities or upper limits for 253 Lyman Break Galaxies (LBGs). The LBGs are a diverse population but with properties correlated with luminosity. The LBGs show a factor of 30 range in indicated stellar mass and a factor of 10 in apparent dust content relative to stellar mass. About 5% of LBGs are luminous at all wavelengths with powerful emission at rest 6\micron. In the rest 0.9 to 2\micron spectral range these galaxies have stellar spectral slopes with no sign of an AGN power law component, suggesting that their emission is mainly powered by intensive star formation. Galaxies in this luminous population share the infrared properties of cold SCUBA sources: both are massive and dusty starburst galaxies at $2<z<3$; their stellar mass is larger than $10^{11} M_{\odot}$. We suggest that these galaxies are the progenitors of present-day giant elliptical galaxies, with a substantial fraction of their stars already formed at $z \approx 3$.

We have examined some basic properties of damped Ly$\alpha$ systems(DLAs) by semi-analytic model. We assume that DLA hosts are disk galaxies whose mass function is generated by Press-Schechter formulism at redshift 3. Star formation and chemical evolution undergo in the disc. We select modelled DLAs according to their observational criterion by Monte Carlo simulation using random line of sights and disk inclinations. The DLA ages are set to be 1 to 3 Gyr. By best-fitting the predicted metallicity distribution to the observed ones, we get the effective yield for DLAs about $0.25Z_{\odot}$. On the basis of this constrain, we further compared our model predictions with observations at redshift 3 in the following items: number density; gas content; HI frequency distribution; star formation rate density; relationship between metallicity and HI column density. We found that the predicted number density at redshift 3 agree well with the observed value, but the gas content $\Omega_{DLA}$ is about 3 times larger than observed since our model predicts more DLA systems with higher column density. The frequency distribution at higher HI column density is quite consistent with observation while some difference exists at lower HI end. The predicted star formation rate density contributed by DLAs is consistent with the most recent observations. Meanwhile, the connection between DLAs and Lyman Break galaxies(LBGs) is discussed by comparing their UV luminosity functions which shows that the DLAs host galaxies are much fainter than LBGs. However, there is a discrepancy between model prediction and observation in the correlation between metallicity and HI column density for DLAs. Further investigations are needed for the star formation mode at high redshift environments.

Based on the disk galaxy formation theory within the framework of standard LCDM hierarchical picture, we selected modelled DLAs, according to their observational criterion, by Monte Carlo simulation with the random inclinations being considered, to examine their observed properties. By best-fitting the predicted metallicity distribution to the observed ones, we get the effective yield for DLAs about 0.25Z_sun, which is comparable to those for SMC and LMC. And the predicted distribution is the same as that of observation at the significant level higher than 60%. The predicted column density distribution of modelled DLAs is compared with the observed ones with the corresponding number density, gas content being discussed. We found that the predicted number density n(z) at redshift 3 agree well with the observed value, but the gas content \Omega_DLA is about 3 times larger than observed since our model predicts more DLA systems with higher column density. It should be noted that the predicted star formation rate density contributed by DLAs is consistent with the most recent observations if the star formation timescale in DLAs is assumed to be 1 to 3 Gyr. Meanwhile, the connection between DLAs and LBGs is discussed by comparing their UV luminosity functions which shows that the DLAs host galaxies are much fainter than LBGs. We also predict that only few percent of DLAs can host LBGs which is also consistent with current observations. However, there is a discrepancy between model prediction and observation in the correlation between metallicity and HI column density for DLAs. We suggest that this could result from either the inadequacy of Schmidt-type star formation law at high redshift, the diversities of DLA populations, or the model limitations.

We analyzed the shapes of Galactic open clusters by the star counting technique with the 2MASS star catalog database. Morphological parameters such as the ellipticity and size have been derived via stellar density distribution, weighed by clustering probability. We find that most star clusters are elongated, even for the youngest star clusters of a few million years old, which are located near to the Galactic disk. The shapes of young star clusters must reflect the conditions in the parental molecular clouds and during the cluster formation process. As an open cluster ages, stellar dynamics cause the inner part of the cluster to circularize, but the overall radius gets larger and the stellar density becomes sparser. We discuss how internal relaxation process competes with Galactic external perturbation during cluster evolution.

We present a newly developed cosmological hydrodynamics code based on weighted essentially non-oscillatory (WENO) schemes for hyperbolic conservation laws. WENO is a higher order accurate finite difference scheme designed for problems with piecewise smooth solutions containing discontinuities, and has been successfully applied for problems involving both shocks and complicated smooth solution structures. We couple hydrodynamics based on the WENO scheme with standard Poisson solver - particle-mesh (PM) algorithm for evolving the self-gravitating system. The third order low storage total variation diminishing (TVD) Runge-Kutta scheme has been used for the time integration of the system. To test accuracy and convergence rate of the code, we subject it to a number of typical tests including the Sod shock tube in multidimensions, the Sedov blast wave and formation of the Zeldovich pancake. These tests validate the WENO hydrodynamics with fast convergence rate and high accuracy. We also evolve a low density flat cosmological model ($\Lambda$CDM) to explore the validity of the code in practical simulations.

We investigate some basic properties of Damped Lyman alpha systems based on the Semi-Analytical model of disk galaxy formation theory. We derive the DLA metallicity, column density, number density, gas content and cosmic star formation rate by assuming that disks form at the center of dark halos, and the modelled DLAs are selected by Monte Carlo simulation according to the distributions of halo properties. We find that DLA hosts are dominated by small galaxies and biased to extended galaxies. In terms of model results, DLAs could naturally arise in a $\Lambda$CDM universe from radiatively cooled gas in dark matter halos. However, model predicts a reverse correlation between metallicity and the column density when compared with observations, regardless of the proposed observational bias. We argue that this could be resulted from the model limitations, or the inadequacy of Schmidt-type star formation mode at high redshift, or/and the diversities of DLA populations.

Galactic winds and mass outflows are observed both in nearby starburst galaxies and in high-redshift star-forming galaxies. In this paper we develop a simple analytic model to understand the observed superwind phenomenon. Our model is built upon the model of McKee & Ostriker (1977) for the interstellar medium. It allows one to predict how properties of a superwind, such as wind velocity and mass outflow rate, are related to properties of its star-forming host galaxy, such as size, gas density and star formation rate. The model predicts a threshold of star formation rate density for the generation of observable galactic winds. Galaxies with more concentrated star formation produce superwinds with higher velocities. The predicted mass outflow rates are comparable to (or slightly larger than) the corresponding star formation rates. We apply our model to both local starburst galaxies and high-redshift Lyman break galaxies, and find its predictions to be in good agreement with current observations. Our model is simple, and so can be easily incorporated into numerical simulations and semi-analytical models of galaxy formation.

We describe a newly developed cosmological hydrodynamics code based on the weighted essentially non-oscillatory (WENO) schemes for hyperbolic conservation laws. High order finite difference WENO schemes are designed for problems with piecewise smooth solutions containing discontinuities, and have been successful in applications for problems involving both shocks and complicated smooth solution structures. We couple hydrodynamics based on the WENO scheme with standard Poisson solver - particle-mesh (PM) algorithm for evolving the self-gravitating system. A third order total variation diminishing (TVD) Runge-Kutta scheme has been used for time-integration of the system. We brief the implementation of numerical technique. The cosmological applications in simulating intergalactic medium and Ly$\alpha$ forest in the CDM scenario are also presented.

Current theory of disk galaxy formation is used to study fundamental-plane (FP) type of relations for disk galaxies. We examine how the changes in model parameters affect these relations and explore the possibility of using such relations to constrain theoretical models. The distribution of galaxy disks in the space of their fundamental properties are predicted to be concentrated in a plane, with the Tully-Fisher (TF) relation (a relation between luminosity $L$ and maximum rotation velocity $V_m$) being an almost edge-on view. Using rotation velocities at larger radii generally leads to larger TF scatter. In searching for a third parameter, we find that both the disk scale-length $R_d$ (or surface brightness) and the rotation-curve shape are correlated with the TF scatter. The FP relation in the $(\Log L, \Log V_m, \Log R_d)$-space obtained from the theory is $L\propto R_d^{\alpha'} V_m^{\beta'}$, with ${\alpha'}\sim 0.50$ and ${\beta'}\sim 2.60$, consistent with the preliminary result we obtain from observational data. Among the model parameters we probe, variation in any of them can generate significant scatter in the TF relation, but the effects of the spin parameter and halo concentration can be reduced significantly by introducing $R_d$ while the scatter caused by varying $m_d$ (the ratio between disk mass and halo mass) is most effectively reduced by introducing the parameters which describes the rotation-curve shape. The TF and FP relations combined should therefore provide useful constraints on models of galaxy formation.

In this paper, we investigate the star formation and chemical evolution of damped Lyman-$\alpha$ systems (DLAs) based on the disc galaxy formation model which is developed by Mo, Mao & White. We propose that the DLAs are the central galaxies of less massive dark haloes present at redshifts $z\sim 3$, and they should inhabit haloes of moderately low circular velocity. The empirical Schmidt law of star formation rates, and closed box model of chemical evolution that an approximation known as instantaneous recycling is assumed, are adopted. In our models, when the predicted distribution of metallicity for DLAs is calculated, two cases are considered. One is that, using the closed box model, empirical Schmidt law and star formation time, the distribution of metallicity can be directly calculated. The other is that, when the simple gravitational instability of a thin isothermal gas disc as first discussed by Toomre is considered, the star formation occurs only in the region where the surface density of gas satisfies the critical value, not everywhere of a gas disc. We assume that star formation in each DLA lasts for a period of 1 Gyr from redshifts $z=3$. There is only one output parameter in our models, i.e. the stellar yield, which relates to the time of star formation history and is obtained by normalizing the predicted distribution of metallicity to the mean value of 1/13 $Z_{\odot}$ as presented by Pettini et al. The predicted metallicity distribution is consistent with the current (rather limited) observational data. A random distribution of galactic discs is taken into account.

We use the observed sizes and star formation rates of Lyman-break galaxies (LBGs) to understand their properties expected in the hierarchical scenario of galaxy formation. The observed size distribution constrains the masses of the host haloes of LBGs from below, while the observed star-formation distribution constrains the masses from above. Assuming a flat CDM model with a cosmological constant ($\Omega_\Lambda=0.7$), we find that consistency with observations requires the circular velocities of LBG host haloes to be in the range 100--$300 {\rm km s^{-1}}$. The predicted comoving correlation length of these objects is $\sim 3 h^{-1}{\rm Mpc}$, and the predicted velocity dispersion of their stellar contents is typically $70 {\rm km s^{-1}}$. The same prescription applied to larger haloes in the CDM cosmogony predicts the existence of galaxies with star formation rates $\sim 1000 {\rm M_\odot yr^{-1}}$ at redshift $z\sim 3$. We explore the possibility of identifying these galaxies to be the bright sub-millimeter (sub-mm) sources detected by SCUBA. The model predicts that the host haloes of these sub-mm sources are massive. These galaxies are predicted to be strongly correlated, with a comoving correlation length of $\sim 7 h^{-1}{\rm Mpc}$. The descendants of the bright sub-mm sources should reside in clusters of galaxies at the present time, and it is likely that these objects are the progenitors of giant ellipticals. The model predicts the existence of a relatively bright and red population of galaxies at $z\sim 3$, which may be observed in the K-band.

The normalized inclination distributions are presented for the spiral galaxies in RC3. The results show that, except for the bin of $81^{\circ}$-$90^{\circ}$, in which the apparent minor isophotal diameters that are used to obtain the inclinations, are affected by the central bulges, the distributions for Sa, Sab, Scd and Sd are well consistent with the Monte-Carlo simulation of random inclinations within 3-$\sigma$, and Sb and Sbc almost, but Sc is different. One reason for the difference between the real distribution and the Monte-Carlo simulation of Sc may be that some quite inclined spirals, the arms of which are inherently loosely wound on the galactic plane and should be classified to Sc galaxies, have been incorrectly classified to the earlier ones, because the tightness of spiral arms which is one of the criteria of the Hubble classification in RC3 is different between on the galactic plane and on the tangent plane of the celestial sphere. Our result also implies that there might exist biases in the luminosity functions of individual Hubble types if spiral galaxies are only classified visually.

This paper presents some statistical correlations of 72 northern spiral galaxies. The results show that early-type spirals that are brighter, and thicker, and the axis ratios ($H_{z}/H_r$) of the disk tend to be smaller along the Hubble sequence. We also find that $H_{z}/H_r$ correlates strongly with the galaxy's color, and early-type spirals have larger values of $H_{z}/H_r$. The inclinations obtained by fitting the pattern of a spiral structure with a logarithmic spiral form are nearly the same as those obtained by using the formulas of Aaronson et al. (1980). Finally, the mean measured pitch angles for the different Hubble sequences in the Third Catalogue of Bright Galaxies by de Vaucouleurs et al. (1991) are derived.

Photographic photometries of 789 stars in the region of the open clusters NGC1750 and NGC1758 are derived from a set of 8 photometric plates. According to the astronomical data including proper motions, positions and membership probabilities of individual stars, several astrophysical parameters for these two clusters are determined, such as HR diagrams, ages, distances, luminosity functions, masses and kinematics, etc. It is found that the distances are ($525\pm 48$pc) and ($794\pm 73$pc) with ages of $1.5\times10^8$yr and $6.3\times10^8$yr for NGC1750 and NGC1758 respectively. Both clusters show no significant mass segregation effects because of their relatively young ages. The analysis of the proper motions imply that their velocity distributions are isotropy.

The number density and clustering properties of Lyman-break galaxies (LBGs) observed at redshift $z\sim 3$ are best explained by assuming that they are associated with the most massive haloes at $z\sim 3$ predicted in hierarchical models of structure formation. In this paper we study, under the same assumption, how star formation and chemical enrichment may have proceeded in the LBG population. A consistent model, in which the amount of cold gas available for star formation must be regulated, is suggested. It is found that gas cooling in dark haloes provides a natural regulation process. In this model, the star formation rate in an LBG host halo is roughly constant over about 1 Gyr. The predicted star formation rates and effective radii are consistent with observations. The metallicity of the gas associated with an LBG is roughly equal to the chemical yield, or about the order of $1 Z_{\odot}$ for a Salpeter IMF. The contribution to the total metals of LBGs is roughly consistent with that obtained from the observed cosmic star formation history. The model predicts a marked radial metallicity gradient in a galaxy, with the gas in the outer region having much lower metallicity. As a result, the metallicities for the damped Lyman-alpha absorption systems expected from the LBG population are low. Since LBG halos are filled with hot gas in this model, their contributions to the soft X-ray background and to the UV ionization background are calculated and discussed.

In the present paper, we introduce a two-component model of the Galactic disk to investigate its chemical evolution. The formation of the thick and thin disks occur in two main accretion episodes with both infall rates to be Gaussian. Both the pre-thin and post-thin scenarios for the formation of the Galactic disk are considered. The best-fitting is obtained through $\chi^2$-test between the models and the new observed metallicity distribution function of G dwarfs in the solar neighbourhood (Hou et al 1998). Our results show that post-thin disk scenario for the formation of the Galactic disk should be preferred. Still, other comparison between model predictions and observations are given.