Showing 1–50 of 57 results for author: Shu, C

A comparative study of ultraluminous infrared galaxies in the IRAS and SDSS Surveys

Authors: Shaohua Zhang, Zhijian Luo, Xiheng Shi, Chenggan Shu, Hubing Xiao, Hongyan Zhou

Abstract: We present a comprehensive study of Ultraluminous Infrared Galaxies (ULIRGs), leveraging data from the IRAS Faint Source Catalogue (FSC) and the spectroscopic catalog in the Sloan Digital Sky Survey (SDSS) DR16. Our meticulous cross-matching technique significantly enhances the reliability of ULIRG identification, resulting in the identification of 283 reliable ULIRGs, including 102 new detections… ▽ More We present a comprehensive study of Ultraluminous Infrared Galaxies (ULIRGs), leveraging data from the IRAS Faint Source Catalogue (FSC) and the spectroscopic catalog in the Sloan Digital Sky Survey (SDSS) DR16. Our meticulous cross-matching technique significantly enhances the reliability of ULIRG identification, resulting in the identification of 283 reliable ULIRGs, including 102 new detections, while discarding 120 previously reported false sources. Covering a redshift range of $z = 0.018 - 0.996$, with a median redshift of $\bar{z} = 0.259$, our uniform sample reveals apparent interaction features in approximately 40\% of ULIRGs, increasing to 92\% for those with $z < 0.1$. Through optical spectra analysis, it is indicated that over 58\% of ULIRGs host an AGN, which is twice as high as the detections based solely on infrared colors. Moreover, a pronounced excess of radio emissions associated with AGN activity results in a steeper radio-far-infrared correlation. Notably, Type I ULIRGs exhibit properties similar to those of narrow-line Seyfert 1 galaxies (NLS1s), with an elevated incidence rate of \ion{Mg}{2} BALs (16.7\%), surpassing that of typical optically selected quasars by over tenfold, consistent with current evolutionary models. We anticipate that forthcoming telescopes such as the China Space Station Telescope (CSST) and Leighton Chajnantor Telescope (LCT) will provide deeper insights into ULIRG morphology, dust distribution, molecular gas, and AGN activity. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: Accepted for publication in ApJS, 33 pages, 13 figures, and 1 table

Acceleration of the Universe without the Hubble tension with Kaniadakis holographic dark energy using the Hubble horizon as the IR cut-off

Authors: Wei Fang, Guo Chen, Chao-Jun Feng, Wei Du, Chenggang Shu

Abstract: We introduce a holographic dark energy model that incorporates the first-order approximate Kaniadaski entropy, utilizing the Hubble horizon, $1/H$, as the infrared cutoff. We investigate the cosmological evolution within this framework. The model introduces an extra parameter relative to the $Λ$CDM model. It posits a Universe that is initially dominated by dark matter, which then evolves to a phas… ▽ More We introduce a holographic dark energy model that incorporates the first-order approximate Kaniadaski entropy, utilizing the Hubble horizon, $1/H$, as the infrared cutoff. We investigate the cosmological evolution within this framework. The model introduces an extra parameter relative to the $Λ$CDM model. It posits a Universe that is initially dominated by dark matter, which then evolves to a phase where dark energy becomes the predominant component, with this transition occurring at a redshift of approximately $z \sim 0.419$. The energy density of dark energy is ultimately expected to become constant, thereby circumventing the potential issue of a "big rip". Employing the most recent Type Ia supernova and Hubble parameter data, we constrain the model's parameters and find a Hubble constant of $H_0=72.8$ km/s/Mpc, thereby resolving the Hubble tension issue. The estimated age of the Universe, based on the best-fit parameter values, is $14.2$ Gyr. Furthermore, we predict the number of strong gravitational lenses and conduct statefinder and $Om$ diagnostic analyses to validate and characterize the model. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: 9 pages, 5 figures

Imputation of Missing Photometric Data and Photometric Redshift Estimation for CSST

Authors: Zhijian Luo, Zhirui Tang, Zhu Chen, Liping Fu, Wei Du, Shaohua Zhang, Yan Gong, Chenggang Shu, Junhao Lu, Yicheng Li, Xian-Min Meng, Xingchen Zhou, Zuhui Fan

Abstract: Accurate photometric redshift (photo-$z$) estimation requires support from multi-band observational data. However, in the actual process of astronomical observations and data processing, some sources may have missing observational data in certain bands for various reasons. This could greatly affect the accuracy and reliability of photo-$z$ estimation for these sources, and even render some estimat… ▽ More Accurate photometric redshift (photo-$z$) estimation requires support from multi-band observational data. However, in the actual process of astronomical observations and data processing, some sources may have missing observational data in certain bands for various reasons. This could greatly affect the accuracy and reliability of photo-$z$ estimation for these sources, and even render some estimation methods unusable. The same situation may exist for the upcoming Chinese Space Station Telescope (CSST). In this study, we employ a deep learning method called Generative Adversarial Imputation Networks (GAIN) to impute the missing photometric data in CSST, aiming to reduce the impact of data missing on photo-$z$ estimation and improve estimation accuracy. Our results demonstrate that using the GAIN technique can effectively fill in the missing photometric data in CSST. Particularly, when the data missing rate is below 30\%, the imputation of photometric data exhibits high accuracy, with higher accuracy in the $g$, $r$, $i$, $z$, and $y$ bands compared to the $NUV$ and $u$ bands. After filling in the missing values, the quality of photo-$z$ estimation obtained by the widely used Easy and Accurate Zphot from Yale (EAZY) software is notably enhanced. Evaluation metrics for assessing the quality of photo-$z$ estimation, including the catastrophic outlier fraction ($f_{out}$), the normalized median absolute deviation ($\rm {σ_{NMAD}}$), and the bias of photometric redshift ($bias$), all show some degree of improvement. Our research will help maximize the utilization of observational data and provide a new method for handling sample missing values for applications that require complete photometry data to produce results. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Estimating Photometric Redshift from Mock Flux for CSST Survey by using Weighted Random Forest

Authors: Junhao Lu, Zhijian Luo, Zhu Chen, Liping Fu, Wei Du, Yan Gong, Yicheng Li, Xian-Min Meng, Zhirui Tang, Shaohua Zhang, Chenggang Shu, Xingchen Zhou, Zuhui Fan

Abstract: Accurate estimation of photometric redshifts (photo-$z$) is crucial in studies of both galaxy evolution and cosmology using current and future large sky surveys. In this study, we employ Random Forest (RF), a machine learning algorithm, to estimate photo-$z$ and investigate the systematic uncertainties affecting the results. Using galaxy flux and color as input features, we construct a mapping bet… ▽ More Accurate estimation of photometric redshifts (photo-$z$) is crucial in studies of both galaxy evolution and cosmology using current and future large sky surveys. In this study, we employ Random Forest (RF), a machine learning algorithm, to estimate photo-$z$ and investigate the systematic uncertainties affecting the results. Using galaxy flux and color as input features, we construct a mapping between input features and redshift by using a training set of simulated data, generated from the Hubble Space Telescope Advanced Camera for Surveys (HST-ACS) and COSMOS catalogue, with the expected instrumental effects of the planned China Space Station Telescope (CSST). To improve the accuracy and confidence of predictions, we incorporate inverse variance weighting and perturb the catalog using input feature errors. Our results show that weighted RF can achieve a photo-$z$ accuracy of $\rm σ_{NMAD}=0.025$ and an outlier fraction of $\rm η=2.045\%$, significantly better than the values of $\rm σ_{NMAD}=0.043$ and $\rm η=6.45\%$ obtained by the widely used Easy and Accurate Zphot from Yale (EAZY) software which uses template-fitting method. Furthermore, we have calculated the importance of each input feature for different redshift ranges and found that the most important input features reflect the approximate position of the break features in galaxy spectra, demonstrating the algorithm's ability to extract physical information from data. Additionally, we have established confidence indices and error bars for each prediction value based on the shape of the redshift probability distribution function, suggesting that screening sources with high confidence can further reduce the outlier fraction. △ Less

Submitted 25 December, 2023; v1 submitted 20 December, 2023; originally announced December 2023.

Black holes regulate cold gas accretion in massive galaxies

Authors: Tao Wang, Ke Xu, Yuxuan Wu, Yong Shi, David Elbaz, Luis C. Ho, Zhi-Yu Zhang, Qiusheng Gu, Yijun Wang, Chenggang Shu, Feng Yuan, Xiaoyang Xia, Kai Wang

Abstract: Nearly every massive galaxy contains a supermassive black hole (BH) at its center. For decades, both theory and numerical simulations have indicated that BHs play a central role in regulating the growth and quenching of galaxies. Specifically, BH feedback by heating or blowing out the interstellar medium (ISM) serves as the groundwork for current models of massive galaxy formation. However, direct… ▽ More Nearly every massive galaxy contains a supermassive black hole (BH) at its center. For decades, both theory and numerical simulations have indicated that BHs play a central role in regulating the growth and quenching of galaxies. Specifically, BH feedback by heating or blowing out the interstellar medium (ISM) serves as the groundwork for current models of massive galaxy formation. However, direct evidence for such an impact on the galaxy-wide ISM from BHs has only been found in some extreme objects. For general galaxy populations, it remains unclear whether and how BHs impact the ISM. Here based on a large sample of nearby galaxies with measurements of masses of both black holes and atomic hydrogen, the major component of cold ISM, we reveal that the atomic hydrogen content ($f_{\rm HI} = M_{\rm HI}/M_{\star}$) is tightly and anti-correlated with black hole mass ($M_{\rm BH}$) with $f_{\rm HI} \propto M^{-α}_{\rm BH}$ ($α\sim 0.5-0.6$). This correlation is valid across five orders of magnitude in $M_{\rm BH}$. Once this correlation is taken into account, $f_{\rm HI}$ loses dependence on other galactic parameters, demonstrating that $M_{\rm BH}$ serves as the primary driver of $f_{\rm HI}$. These findings provide critical evidence for how the accumulated energy from BH accretion impacts galaxy-wide ISM, representing a crucial step forward in our understanding on the role of BHs in regulating the growth and quenching of massive galaxies. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Comments: 24 pages, 7 figures. Submitted to Nature

Mass Reconstruction of Galaxy-scale Strong Gravitational Lenses Using a Broken Power-law Model

Authors: Wei Du, Liping Fu, Yiping Shu, Ran Li, Zuhui Fan, Chenggang Shu

Abstract: With mock strong gravitational lensing images, we investigate the performance of the broken power-law (BPL) model proposed by \citet{2020ApJ...892...62D} on the mass reconstruction of galaxy-scale lenses. An end-to-end test is carried out, including the creation of mock strong lensing images, the subtraction of lens light, and the reconstruction of lensed images, where the lenses are selected from… ▽ More With mock strong gravitational lensing images, we investigate the performance of the broken power-law (BPL) model proposed by \citet{2020ApJ...892...62D} on the mass reconstruction of galaxy-scale lenses. An end-to-end test is carried out, including the creation of mock strong lensing images, the subtraction of lens light, and the reconstruction of lensed images, where the lenses are selected from the galaxies in the Illustris-1 simulation. We notice that, regardless of the adopted mass models (the BPL model or its special cases), the Einstein radius can be robustly determined from imaging data alone, and the median bias is typically less than $1\%$. Away from the Einstein radius, the lens mass distribution tends to be harder to measure, especially at radii where there are no lensed images detected. We find that, with rigid priors, the BPL model can clearly outperform the single power-law models by achieving $<5\%$ median bias on the radial convergence profile within the Einstein radius. As for the source light reconstructions, they are found to be sensitive to both lens light contamination and lens mass models, where the BPL model with rigid priors still performs best when there is no lens light contamination. We show that, by correcting for the projection effect, the BPL model can estimate the aperture and luminosity weighted line-of-sight velocity dispersions to an accuracy of $\sim6\%$ scatter. These results highlight the great potential of the BPL model in strong lensing related studies. △ Less

Submitted 19 August, 2023; v1 submitted 9 February, 2023; originally announced February 2023.

Comments: Accepted for publication in ApJ, 24 pages, 13 figures, 2 tables

Provably Positive Central DG Schemes via Geometric Quasilinearization for Ideal MHD Equations

Authors: Kailiang Wu, Haili Jiang, Chi-Wang Shu

Abstract: In the numerical simulation of ideal MHD, keeping the pressure and density positive is essential for both physical considerations and numerical stability. This is a challenge, due to the underlying relation between such positivity-preserving (PP) property and the magnetic divergence-free (DF) constraint as well as the strong nonlinearity of the MHD equations. This paper presents the first rigorous… ▽ More In the numerical simulation of ideal MHD, keeping the pressure and density positive is essential for both physical considerations and numerical stability. This is a challenge, due to the underlying relation between such positivity-preserving (PP) property and the magnetic divergence-free (DF) constraint as well as the strong nonlinearity of the MHD equations. This paper presents the first rigorous PP analysis of the central discontinuous Galerkin (CDG) methods and constructs arbitrarily high-order PP CDG schemes for ideal MHD. By the recently developed geometric quasilinearization (GQL) approach, our analysis reveals that the PP property of standard CDG methods is closely related to a discrete DF condition, whose form was unknown and differs from the non-central DG and finite volume cases in [K. Wu, SIAM J. Numer. Anal. 2018]. This result lays the foundation for the design of our PP CDG schemes. In 1D case, the discrete DF condition is naturally satisfied, and we prove the standard CDG method is PP under a condition that can be enforced with a PP limiter. However, in the multidimensional cases, the discrete DF condition is highly nontrivial yet critical, and we prove the the standard CDG method, even with the PP limiter, is not PP in general, as it fails to meet the discrete DF condition. We address this issue by carefully analyzing the structure of the discrete divergence and then constructing new locally DF CDG schemes for Godunov's modified MHD equations with an additional source. The key point is to find out the suitable discretization of the source term such that it exactly offsets all the terms in the discrete DF condition. Based on the GQL approach, we prove the PP property of the new multidimensional CDG schemes. The robustness and accuracy of PP CDG schemes are validated by several demanding examples, including the high-speed jets and blast problems with very low plasma beta. △ Less

Submitted 28 March, 2022; originally announced March 2022.

The 2175 Å~ Bump Features in FeLoBAL Quasars: One Indicator of MW-like Dust in the Nuclear Region of Quasar

Authors: Shaohua Zhang, Jian Ge, Tuo Ji, Peng Jiang, Zhijian Luo, Xiang Pan, Xiheng Shi, Chenggang Shu, Jianguo Wang, Hubing Xiao, Hongyan Zhou

Abstract: To investigate the properties of dust in the nuclear region of quasars, we explored the extinction curves of the iron low-ionization broad absorption line (FeLoBAL) quasar SDSS J163004.29+311957.6 and its two analogues. The parameterized extinction curves indicated the Milky Way-like 2175 Å~ bump features in underlying extinction, which are similar to those seen in the Local Group and a subset of… ▽ More To investigate the properties of dust in the nuclear region of quasars, we explored the extinction curves of the iron low-ionization broad absorption line (FeLoBAL) quasar SDSS J163004.29+311957.6 and its two analogues. The parameterized extinction curves indicated the Milky Way-like 2175 Å~ bump features in underlying extinction, which are similar to those seen in the Local Group and a subset of high-redshift star-forming galaxies. Compared to the bump features in the Large Magellanic Clouds (LMC), the detections in this work are much closer to those in the Milky Way (MW). These bump features, as well as those in the high- and low-ionization broad absorption line (BAL) quasars of Zhang et al., are probably the counterpart of the 2175 Å~ bump features in the quasar environment. This type of dust grain is generally small, easily disrupted by high-energy photons and has difficulty surviving in the radiation field of the active galactic nucleus (AGN). However, due to the presence of absorption-line outflows, the 2175 Å~ bump feature in quasars, which should be rare, is seen many times in BAL quasars. The shielding effect of outflow clouds allows the MW-like dust grains to be assembled or extends the survival period in the quasar nuclear region. The process, and physical and chemical conditions deserve further observational study and investigation. △ Less

Submitted 28 March, 2022; v1 submitted 7 March, 2022; originally announced March 2022.

Comments: 13 pages, 10 figures. Accepted by Astronomy & Astrophysics

Journal ref: A&A 663, A63 (2022)

Provably Physical-Constraint-Preserving Discontinuous Galerkin Methods for Multidimensional Relativistic MHD Equations

Authors: Kailiang Wu, Chi-Wang Shu

Abstract: 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… ▽ More 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. △ Less

Submitted 9 February, 2020; originally announced February 2020.

Quenching as a Contest between Galaxy Halos and their Central Black Holes

Authors: Zhu Chen, S. M. Faber, David. C. Koo, Rachel S. Somerville, Joel R. Primack, Avishai Dekel, Aldo Rodríguez-Puebla, Yicheng Guo, Guillermo Barro, Dale D. Kocevski, A. van der Wel, Joanna Woo, Eric F. Bell, Jerome J. Fang, Henry C. Ferguson, Mauro Giavalisco, Marc Huertas-Company, Fangzhou Jiang, Susan Kassin, Lin Lin, Fengshan Liu, Yifei Luo, Zhijian Luo, Camilla Pacifici, Viraj Pandya , et al. (5 additional authors not shown)

Abstract: 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-formi… ▽ More 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. △ Less

Submitted 26 May, 2020; v1 submitted 24 September, 2019; originally announced September 2019.

Comments: 45 pages, 20 figures, 2 tables, accepted for publication in ApJ

Entropy Symmetrization and High-Order Accurate Entropy Stable Numerical Schemes for Relativistic MHD Equations

Authors: Kailiang Wu, Chi-Wang Shu

Abstract: 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 admit a thermodynamic entropy pair. To address this issue, a symmetrizable RMHD system, equipped with a convex thermodynamic entropy pair, is proposed by adding a source… ▽ More 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 admit a thermodynamic entropy pair. To address this issue, a symmetrizable RMHD system, equipped with a convex thermodynamic entropy pair, is proposed by adding a source term into the equations, providing an analogue to the non-relativistic Godunov--Powell system. Arbitrarily high-order accurate entropy stable finite difference schemes are 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 numerical tests demonstrate the accuracy and robustness of the proposed entropy stable schemes. △ Less

Submitted 9 July, 2021; v1 submitted 17 July, 2019; originally announced July 2019.

Comments: 32 pages; 11 figures; add some references on entropy symmetrization

Journal ref: SIAM Journal on Scientific Computing, 42(4): A2230-A2261, 2020

Provably Positive High-Order Schemes for Ideal Magnetohydrodynamics: Analysis on General Meshes

Authors: Kailiang Wu, Chi-Wang Shu

Abstract: 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… ▽ More 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. △ Less

Submitted 30 July, 2018; originally announced July 2018.

Comments: 49 pages, 11 figures

Journal ref: Numerische Mathematik 142, 995-1047, 2019

Provably Positive Discontinuous Galerkin Methods for Multidimensional Ideal Magnetohydrodynamics

Authors: Kailiang Wu, Chi-Wang Shu

Abstract: 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 de… ▽ More 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. △ Less

Submitted 30 June, 2018; originally announced July 2018.

Comments: Submitted for publication on Jan. 31, 2018

Journal ref: SIAM Journal on Scientific Computing, 40(5), B1302-B1329, 2018

Galaxy-galaxy weak-lensing measurement from SDSS: II. host halo properties of galaxy groups

Authors: Wentao Luo, Xiaohu Yang, Tianhuan Lu, Feng Shi, Jun Zhang, H. J. Mo, Chenggang Shu, Liping Fu, Mario Radovich, Jiajun Zhang, Nan Li, Tomomi Sunayama, Lei Wang

Abstract: 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,… ▽ More 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 △ Less

Submitted 25 December, 2017; originally announced December 2017.

Comments: 12 pages, 7 figures, submitted to ApJ

Origin of the fundamental plane of elliptical galaxies in the Coma Cluster without fine-tuning

Authors: Mu-Chen Chiu, Chung-Ming Ko, Chenggang Shu

Abstract: 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 ca… ▽ More 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. △ Less

Submitted 17 March, 2017; v1 submitted 15 March, 2017; originally announced March 2017.

Comments: The paper is accepted for Physical Review D

Galaxy-Galaxy Weak Lensing Measurements from SDSS: I. Image Processing and Lensing signals

Authors: Wentao Luo, Xiaohu Yang, Jun Zhang, Dylan Tweed, Liping Fu, H. J. Mo, Frank C. van den Bosch, Chenggang Shu, Ran Li, Nan Li, Xiangkun Liu, Chuzhong Pan, Yiran Wang, Mario Radovich

Abstract: 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… ▽ More 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. △ Less

Submitted 19 July, 2016; originally announced July 2016.

Comments: 19pages, 12 figures ApJ submitted

Cosmological Perturbations and Quasi-Static Assumption in $f(R)$ Theories

Authors: Mu-Chen Chiu, Andy Taylor, Chenggang Shu, Hong Tu

Abstract: $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)… ▽ More $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, $Ψ$, and the curvature potential, $Φ$. 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$. △ Less

Submitted 7 August, 2015; v1 submitted 12 May, 2015; originally announced May 2015.

Full Investigation on the Dynamics of Power-Law Kinetic Quintessence

Authors: Wei Fang, Hong Tu, Ying Li, Jiasheng Huang, Chenggang Shu

Abstract: We give a full investigation on the dynamics of power-law kinetic quintessence $L(X, φ)=V(φ)(-X+X^2)$ by considering the potential related parameter $Γ$($=\frac{V V''}{V'^2}$) as a function of another potential parameter $λ$($=\frac{V'}{κ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 f… ▽ More We give a full investigation on the dynamics of power-law kinetic quintessence $L(X, φ)=V(φ)(-X+X^2)$ by considering the potential related parameter $Γ$($=\frac{V V''}{V'^2}$) as a function of another potential parameter $λ$($=\frac{V'}{κ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 $Ω_φ$ 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. △ Less

Submitted 1 June, 2014; originally announced June 2014.

Comments: 17 pages, 4 figures, accepted by Physical Review D, in press

Journal ref: Phys. Rev. D 89, 123514(2014)

Dynamical System of Scalar Field from 2-Dimension to 3-D and its Cosmological Implication

Authors: Wei Fang, Hong Tu, Jiasheng Huang, Chenggang Shu

Abstract: 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, λ)$ to observable related variables $(w_φ, Ω_φ, λ)$, and show the intimate relationships between those scalar fields that… ▽ More 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, λ)$ to observable related variables $(w_φ, Ω_φ, λ)$, 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_φ, Ω_φ, λ)$ instead of variables $(x, y, λ)$ 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. △ Less

Submitted 12 September, 2016; v1 submitted 17 February, 2014; originally announced February 2014.

Comments: 20 pages, 8 figures,some references added. Minor changes according to the suggestions from referees

Journal ref: Eur. Phys. J. C (2016) 76:492

Turbulence in the Inter-galactic Medium: Solenoidal and Dilatational Motions, and the Impact of Numerical Viscosity

Authors: Weishan Zhu, Long-Long Feng, Yinhua Xia, Chi-Wang Shu, Qiu-Sheng Gu, Li-Zhi Fang

Abstract: 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 reache… ▽ More 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 $ρ_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. △ Less

Submitted 21 August, 2013; originally announced August 2013.

Comments: 36 pages, 14 figures, Accepted for publication in The Astrophysical Journal

Journal ref: ApJ, 777, 48 (2013)

SED-inferred properties and morphology of Lyman-break galaxies at $z\sim 1$ in the CDF-S

Authors: Z. Chen, C. G. Shu, D. Burgarella, V. Buat, J. -S. Huang, Z. J. Luo

Abstract: 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… ▽ More 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. △ Less

Submitted 16 February, 2013; originally announced February 2013.

Comments: 26 pages, 30 figures, 6 tables, accepted by MNRAS

Multi-Wavelength Study of a Complete IRAC 3.6micron-Selected Galaxy Sample: a Fair Census of Red and Blue Populations at Redshifts 0.4-1

Authors: J. -S. Huang, S. M. Faber, C. N. A. Willmer, D. Rigopoulou, D. Koo, J. Newman, C. Shu, M. L. N. Ashby, P. Barmby, A. Coil, Z. Luo, G. Magdis, T. Wang, B. Weiner, S. P. Willner, X. Z. Zheng, G. G. Fazio

Abstract: We present a multi-wavelength study of a 3.6 $μ$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 $μ$m band measures the rest-frame near-infrared band, permitting nearly unbiased selection with respect to both quiescent and star-forming galaxies. The num… ▽ More We present a multi-wavelength study of a 3.6 $μ$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 $μ$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σ$ 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. △ Less

Submitted 5 December, 2012; originally announced December 2012.

Comments: 64 pages, 24 figures

New determination of the 13C(a, n)16O reaction rate and its influence on the s-process nucleosynthesis in AGB stars

Authors: B. Guo, Z. H. Li, M. Lugaro, J. Buntain, D. Y. Pang, Y. J. Li, J. Su, S. Q. Yan, X. X. Bai, Y. S. Chen, Q. W. Fan, S. J. Jin, A. I. Karakas, E. T. Li, Z. C. Li, G. Lian, J. C. Liu, X. Liu, J. R. Shi, N. C. Shu, B. X. Wang, Y. B. Wang, S. Zeng, W. P. Liu

Abstract: We present a new measurement of the $α$-spectroscopic factor ($S_α$) and the asymptotic normalization coefficient (ANC) for the 6.356 MeV 1/2$^+$ subthreshold state of $^{17}$O through the $^{13}$C($^{11}$B, $^{7}$Li)$^{17}$O transfer reaction and we determine the $α$-width of this state. This is believed to have a strong effect on the rate of the $^{13}$C($α$, $n$)$^{16}$O reaction, the main neut… ▽ More We present a new measurement of the $α$-spectroscopic factor ($S_α$) and the asymptotic normalization coefficient (ANC) for the 6.356 MeV 1/2$^+$ subthreshold state of $^{17}$O through the $^{13}$C($^{11}$B, $^{7}$Li)$^{17}$O transfer reaction and we determine the $α$-width of this state. This is believed to have a strong effect on the rate of the $^{13}$C($α$, $n$)$^{16}$O reaction, the main neutron source for {\it slow} neutron captures (the $s$-process) in asymptotic giant branch (AGB) stars. Based on the new width we derive the astrophysical S-factor and the stellar rate of the $^{13}$C($α$, $n$)$^{16}$O reaction. At a temperature of 100 MK our rate is roughly two times larger than that by \citet{cau88} and two times smaller than that recommended by the NACRE compilation. We use the new rate and different rates available in the literature as input in simulations of AGB stars to study their influence on the abundances of selected $s$-process elements and isotopic ratios. There are no changes in the final results using the different rates for the $^{13}$C($α$, $n$)$^{16}$O reaction when the $^{13}$C burns completely in radiative conditions. When the $^{13}$C burns in convective conditions, as in stars of initial mass lower than $\sim$2 $M_\sun$ and in post-AGB stars, some changes are to be expected, e.g., of up to 25% for Pb in our models. These variations will have to be carefully analyzed when more accurate stellar mixing models and more precise observational constraints are available. △ Less

Submitted 3 August, 2012; originally announced August 2012.

Journal ref: Astrophys. J. 756, 193, 2012

Effect of Dust on Lyman-alpha Photon Transfer in Optically Thick Halo

Authors: Yang Yang, Ishani Roy, Chi-Wang Shu, Li-Zhi Fang

Abstract: We investigate the effects of dust on Lyα 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α photons m… ▽ More We investigate the effects of dust on Lyα 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α photons makes the problem non-trivial. For instance, if the medium has the optical depth of dust absorption and scattering to be τa>>1, τ>>1, and τ>>τa, the effective absorption optical depth in a random walk scenario would be equal to \surd τa(τa+τ). We show, however, that for a resonant scattering at frequency ν0, the effective absorption optical depth would be even larger than τ(ν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α 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α photon source with optical thick halo. △ Less

Submitted 9 July, 2011; originally announced July 2011.

Comments: 24 pages, 8 figures

The mass-dependent star formation histories of disk galaxies: infall model versus observations

Authors: R. X. Chang, J. L. Hou, S. Y. Shen, C. G. Shu

Abstract: 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… ▽ More 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. △ Less

Submitted 17 November, 2010; v1 submitted 18 August, 2010; originally announced August 2010.

Comments: 12 pages, 6 figures, ApJ, 2010, 722, 380

Resonant Scattering and Ly-alpha Radiation Emergent from Neutral Hydrogen Halos

Authors: Ishani Roy, Chi-Wang Shu, Li-Zhi Fang

Abstract: With a state-of-the-art numerical method for solving the integral-differential equation of radiative transfer, we investigate the flux of the Ly$α$ photon $ν_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… ▽ More With a state-of-the-art numerical method for solving the integral-differential equation of radiative transfer, we investigate the flux of the Ly$α$ photon $ν_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 $|ν_{\pm}-ν_0|$, where $ν_{\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 $ν_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. △ Less

Submitted 3 May, 2010; originally announced May 2010.

Comments: 25 pages, 10 figures, accepted for publication in ApJ

The morphological dependent Tully-Fisher relation of spiral galaxies

Authors: Shiyin Shen, Caihong Wang, Ruixiang Chang, Zhengyi Shao, Jinliang Hou, Chenggang Shu

Abstract: 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… ▽ More 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. △ Less

Submitted 26 September, 2009; originally announced September 2009.

Comments: 8 pages, 3 figures, ApJ accepted

Journal ref: Astrophys.J.705:1496-1502,2009

Wouthuysen-Field Coupling in the 21 cm Region Around High Redshift Sources

Authors: I. Roy, W. Xu, J. M. Qiu, C. W. Shu, L. Z. Fang

Abstract: 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… ▽ More 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 τmedia roughly proportional to the optical depth τ, not τ^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. △ Less

Submitted 12 August, 2009; originally announced August 2009.

Comments: 26 pages, 11 figures. Accepted by the Astrophysical journal

A WENO algorithm for radiative transfer with resonant scattering: the time scale of the Wouthuysen-Field Coupling

Authors: Ishani Roy, Jing-Mei Qiu, Chi-Wang Shu, Li-Zhi Fang

Abstract: 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 emissio… ▽ More 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. △ Less

Submitted 14 January, 2009; originally announced January 2009.

Comments: 20 pages, 13 figures

Time evolution of Wouthuysen-Field coupling

Authors: I. Roy, W. Xu, J. -M. Qiu, C. -W. Shu, L. -Z. Fang

Abstract: 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) c… ▽ More 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 ν_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 Λ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 Λ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. △ Less

Submitted 22 December, 2008; originally announced December 2008.

Comments: 22 pages, 25 figures, accepted to the Astrophysical Journal

Journal ref: Astrophys.J.694:1121-1130,2009

Comparisons between isothermal and NFW mass profiles for strong-lensing galaxy clusters

Authors: Chenggang Shu, Binglu Zhou, Matthias Bartelmann, Julia M. Comerford, J. -S. Huang, Yannick Mellier

Abstract: 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 magn… ▽ More 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$. △ Less

Submitted 8 May, 2008; originally announced May 2008.

Comments: 11 pages, 15 figures; accepted for publication in ApJ

Probing the Slope of Cluster Mass Profile with Gravitational Einstein Rings: Application to Abell 1689

Authors: H. Tu, M. Limousin, B. Fort, C. G. Shu, J. F. Sygnet, E. Jullo, J. -P. Kneib, J. Richard

Abstract: 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 arou… ▽ More 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. △ Less

Submitted 10 January, 2008; v1 submitted 11 October, 2007; originally announced October 2007.

Comments: Accepted for publication in MNRAS. Substantial modification after referee's report

A WENO Algorithm for the Growth of Ionized Regions at the Reionization Epoch

Authors: Jing-Mei Qiu, Chi-Wang Shu, Ji-Ren Liu, Li-Zhi Fang

Abstract: 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 th… ▽ More 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. △ Less

Submitted 3 June, 2007; originally announced June 2007.

Comments: Elsart Latex file, 24 pages, 11 figures included, accepted for publication in New Astronomy

The Differences of Star Formation History Between Merging Galaxies and Field Galaxies in the EDR of the SDSS

Authors: Z. J. Luo, C. G. Shu, J. S. Huang

Abstract: 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-Å break strength, the Balmer absorption-line index, and the specific star formation rate. It is found that for early-type… ▽ More 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-Å 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. △ Less

Submitted 14 March, 2007; originally announced March 2007.

Comments: 9 pages, 4 figures, accepted for publication in PASJ

21 cm Signals from Early Ionizing Sources

Authors: Jiren Liu, Jing-Mei Qiu, Long-Long Feng, Chi-Wang Shu, Li-Zhi Fang

Abstract: 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… ▽ More 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. △ Less

Submitted 13 March, 2007; originally announced March 2007.

Comments: 25 pages, 9 figures, accepted by ApJ

Journal ref: Astrophys.J.663:1-9,2007

A WENO Algorithm of the Temperature and Ionization Profiles around a Point Source

Authors: Jing-Mei Qiu, Long-Long Feng, Chi-Wang Shu, Li-Zhi Fang

Abstract: 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 i… ▽ More 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. △ Less

Submitted 28 December, 2006; originally announced December 2006.

Comments: 24 pages, 7 figures, accepted in New Astronomy

Report number: NEWAST-D-06-00121R1

Journal ref: NewAstron.12:398-409,2007

AEGIS: Infrared Spectroscopy of An Infrared Luminous Lyman Break Galaxy at z=3.01

Authors: J. -S. Huang, D. Rigopoulou, C. Papovich, M. L. N. Ashby, S. P. Willner, R. Ivison, E. S. Laird, T. Webb, G. Wilson, P. Barmby, S. Chapman, C. Conselice, B. Mcleod, C. G. Shu, H. A. Smith, E. Le Floc'h, E. Egami, C. A. N. Willmer, G. Fazio

Abstract: 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… ▽ More 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. △ Less

Submitted 21 August, 2006; originally announced August 2006.

Color-magnitude relations of late-type galaxies

Authors: Ruixiang Chang, Shiyin Shen, Jinliang Hou, Chenggang Shu, Zhengyi Shao

Abstract: 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 an… ▽ More 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. △ Less

Submitted 20 July, 2006; originally announced July 2006.

Comments: 16 pages, 8 figures, accepted by MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.372:199-208,2006

SPITZER Observations of z~3 Lyman Break Galaxies: stellar masses and mid-infrared properties

Authors: D. Rigopoulou, J. -S. Huang, C. Papovich, M. L. N. Ashby, P. Barmby, C. Shu, K. Bundy, E. Egami, G. Magdis, H. Smith, S. P. Willner, G. Wilson, G. G. Fazio

Abstract: 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 us… ▽ More 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. △ Less

Submitted 15 May, 2006; originally announced May 2006.

Comments: accepted for publication by Astrophysical Journal

Journal ref: Astrophys.J.648:81-94,2006

A WENO Algorithm for the Radiative Transfer and Ionized Sphere at Reionization

Authors: Jing-Mei Qiu, Chi-Wang Shu, Long-Long Feng, Li-Zhi Fang

Abstract: 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 s… ▽ More 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. △ Less

Submitted 12 April, 2006; originally announced April 2006.

Comments: Elsart Latex file, 20 pages, 8 figures included, accepted for publication in New Astronomy

Low-Redshift Cosmic Baryon Fluid on Large Scales and She-Leveque Universal Scaling

Authors: Ping He, Jiren Liu, Long-Long Feng, Chi-Wang Shu, Li-Zhi Fang

Abstract: 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… ▽ More 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. △ Less

Submitted 9 February, 2006; v1 submitted 10 January, 2006; originally announced January 2006.

Comments: 4 pages, 3 figures

Journal ref: Phys.Rev.Lett. 96 (2006) 051302

Infrared Luminous Lyman Break Galaxies: A Population that Bridges LBGs and SCUBA Galaxies

Authors: J. -S. Huang, D. Rigopoulou, S. P. Willner, C. Papovich, C. Shu, M. L. N. Ashby, P. Barmby, K. Bundy, C. Conselice, E. Egami, P. G. Perez-Gonzalez, J. L. Rosenberg, H. A. Smith, G. Wilson, G. G. Fazio

Abstract: 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… ▽ More 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$. △ Less

Submitted 28 July, 2005; originally announced July 2005.

Comments: Accepted by ApJ

Journal ref: Astrophys.J.634:137-141,2005

Star Formation and Chemical Evolution of DLAs with Semi-Analytical Model

Authors: J. L. Hou, C. G. Shu, S. Y. Shen, R. X. Chang

Abstract: We have examined some basic properties of damped Ly$α$ 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 d… ▽ More We have examined some basic properties of damped Ly$α$ 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 $Ω_{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. △ Less

Submitted 7 May, 2005; originally announced May 2005.

Comments: 6 pages, 3 figures, To appear in "Probing Galaxies through Quasar Absorption Lines", Proceedings IAU Colloquium 199, 2005, Eds. P. R. Williams, C. Shu, and B. Menard

Journal ref: IAUSymp.199:225-230,2005

Metallicity and HI Column Density Properties of Damped Lyman-$α$ Systems

Authors: J. L. Hou, C. G. Shu, S. Y. Shen, R. X. Chang, W. P. Chen, C. Q. Fu

Abstract: 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 a… ▽ More 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 Ω_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. △ Less

Submitted 27 January, 2005; originally announced January 2005.

Comments: 31 pages including 7 figures, will appear in ApJ

Journal ref: Astrophys.J. 624 (2005) 561-570

Morphology of Galactic Open Clusters

Authors: W. P. Chen, C. W. Chen, C. G. Shu

Abstract: 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 t… ▽ More 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. △ Less

Submitted 2 August, 2004; originally announced August 2004.

Comments: 13 pages; 10 pages; accepted Astronomical Journal

A Hybrid Cosmological Hydrodynamic/N-body Code Based on a Weighted Essentially Non-Oscillatory Scheme

Authors: Long-Long Feng, Chi-Wang Shu, Meng-Ping Zhang

Abstract: 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 st… ▽ More 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 ($Λ$CDM) to explore the validity of the code in practical simulations. △ Less

Submitted 14 May, 2004; v1 submitted 10 March, 2004; originally announced March 2004.

Comments: AAS Latex file, 34 pages, 12 figures included, revised version accepted for publication in ApJ

Journal ref: Astrophys.J. 612 (2004) 1-13

Physical properties of DLA metallicity and neutral hydrogen column density

Authors: J. L. Hou, C. G. Shu, W. P. Chen, R. X. Chang, C. Q. Fu

Abstract: 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… ▽ More 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 $Λ$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. △ Less

Submitted 15 October, 2003; originally announced October 2003.

Comments: 6 pages, 6 figures, to appear in `Recycling Intergalactic and Interstellar Matter', IAU Symp. Ser., vol. 217, eds. P.-A. Duc, J. Braine and E. Brinks, 2004

An analytic model for the galactic winds and mass outflows

Authors: Chenggang Shu, H. J. Mo, Shude Mao

Abstract: 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… ▽ More 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. △ Less

Submitted 3 January, 2003; originally announced January 2003.

Comments: 26 pages, 10 figures, submitted to MNRAS

Journal ref: Chin.J.Astron.Astrophys.5:327-346,2005

A Hybrid Cosmological Hydrodynamic/N-body Code Based on the Weighted Essentially Non-Oscillatory Scheme

Authors: Long-Long Feng, Chi-Wang Shu, Meng-Ping Zhang

Abstract: 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 s… ▽ More 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$α$ forest in the CDM scenario are also presented. △ Less

Submitted 31 December, 2002; originally announced December 2002.

Comments: 6 pages, 2 figures included, in Proceeding of 5th Sino-Germany Workshop on Cosmology and Galaxy Formation

The Fundamental Plane of Spiral Galaxies: Theoretical Expectations

Authors: Shiyin Shen, H. J. Mo, Chenggang Shu

Abstract: 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, wit… ▽ More 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^{α'} V_m^{β'}$, with ${α'}\sim 0.50$ and ${β'}\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. △ Less

Submitted 6 May, 2001; originally announced May 2001.

Comments: 25 pages, 11 figures, 4 tables; submitted to MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc. 331 (2002) 259