Source author record

Yu Lu

Yu Lu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

55works
17topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

55 published item(s)

preprint2026arXiv

Power Reinforcement Post-Training of Text-to-Image Models with Super-Linear Advantage Shaping

Recently, post-training methods based on reinforcement learning, with a particular focus on Group Relative Policy Optimization (GRPO), have emerged as the robust paradigm for further advancement of text-to-image (T2I) models. However, these methods are often prone to reward hacking, wherein models exploit biases in imperfect reward functions rather than yielding genuine performance gains. In this work, we identify that normalization could lead to miscalibration and directly removing the prompt-level standard deviation term yields an optimal policy ascent direction that is linear in the advantage but still limits the separation of genuine signals from noise. To mitigate the above issues, we propose Super-Linear Advantage Shaping (SLAS) by revisiting the functional update from an information geometry perspective. By extending the Fisher-Rao information metric with advantage-dependent weighting, SLAS introduces a non-linear geometric structure that reshapes the local policy space. This design relaxes constraints along high-advantage directions to amplify informative updates, while tightening those in low-advantage regions to suppress illusory gradients. In addition, batch-level normalization is applied to stabilize training under varying reward scales. Extensive evaluations demonstrate that SLAS consistently surpasses the DanceGRPO baseline across multiple backbones and benchmarks. In particular, it yields faster training dynamics, improved out-of-domain performance on GenEval and UniGenBench++, and enhanced robustness to model scaling, while mitigating reward hacking and preserving semantic and compositional fidelity in generations.

preprint2026arXiv

PsyCLIENT: Client Simulation via Conversational Trajectory Modeling for Trainee Practice and Model Evaluation in Mental Health Counseling

LLM-based client simulation has emerged as a promising tool for training novice counselors and evaluating automated counseling systems. However, existing client simulation approaches face three key challenges: (1) limited diversity and realism in client profiles, (2) the lack of a principled framework for modeling realistic client behaviors, and (3) a scarcity in Chinese-language settings. To address these limitations, we propose PsyCLIENT, a novel simulation framework grounded in conversational trajectory modeling. By conditioning LLM generation on predefined real-world trajectories that incorporate explicit behavior labels and content constraints, our approach ensures diverse and realistic interactions. We further introduce PsyCLIENT-CP, the first open-source Chinese client profile dataset, covering 60 distinct counseling topics. Comprehensive evaluations involving licensed professional counselors demonstrate that PsyCLIENT significantly outperforms baselines in terms of authenticity and training effectiveness. Notably, the simulated clients are nearly indistinguishable from human clients, achieving an about 95\% expert confusion rate in discrimination tasks. These findings indicate that conversational trajectory modeling effectively bridges the gap between theoretical client profiles and dynamic, realistic simulations, offering a robust solution for mental health education and research. Code and data will be released to facilitate future research in mental health counseling.

preprint2022arXiv

Coupled Channel Effects for the Charmed-Strange Mesons

We make a systematic calculation of the spectra and hadronic decays of the $D_s$ system in a coupled channel framework, where the unquenched effects are induced by the $^3P_0$ model. In the calculation, the wave functions are obtained by using a nonrelativistic potential model and are handled precisely with Gaussian Expansion Method. Even though the fitting mainly focuses on the spectrum, our model agrees well with the experiments on both the spectra and the hadronic decays, suggesting that the coupled channel effect could result in a reasonable and coherent description of the $D_s$ mesons. Based on the calculation, we give a detailed analysis on various aspects of the excited states, especially $D_s(2317), D_s(2460), D_s(2536), D_s(2860), D_s(3040)$. We also predict that $D_s(2{}^3P_0)$ should be a $D^*K^*$ dominant molecule with mass $2894$~MeV, which is only $5$~MeV below the $D^*K^*$ threshold.

preprint2022arXiv

CRIS: CLIP-Driven Referring Image Segmentation

Referring image segmentation aims to segment a referent via a natural linguistic expression.Due to the distinct data properties between text and image, it is challenging for a network to well align text and pixel-level features. Existing approaches use pretrained models to facilitate learning, yet separately transfer the language/vision knowledge from pretrained models, ignoring the multi-modal corresponding information. Inspired by the recent advance in Contrastive Language-Image Pretraining (CLIP), in this paper, we propose an end-to-end CLIP-Driven Referring Image Segmentation framework (CRIS). To transfer the multi-modal knowledge effectively, CRIS resorts to vision-language decoding and contrastive learning for achieving the text-to-pixel alignment. More specifically, we design a vision-language decoder to propagate fine-grained semantic information from textual representations to each pixel-level activation, which promotes consistency between the two modalities. In addition, we present text-to-pixel contrastive learning to explicitly enforce the text feature similar to the related pixel-level features and dissimilar to the irrelevances. The experimental results on three benchmark datasets demonstrate that our proposed framework significantly outperforms the state-of-the-art performance without any post-processing. The code will be released.

preprint2022arXiv

Interpretable Fault Diagnosis of Rolling Element Bearings with Temporal Logic Neural Network

Machine learning-based methods have achieved successful applications in machinery fault diagnosis. However, the main limitation that exists for these methods is that they operate as a black box and are generally not interpretable. This paper proposes a novel neural network structure, called temporal logic neural network (TLNN), in which the neurons of the network are logic propositions. More importantly, the network can be described and interpreted as a weighted signal temporal logic. TLNN not only keeps the nice properties of traditional neuron networks but also provides a formal interpretation of itself with formal language. Experiments with real datasets show the proposed neural network can obtain highly accurate fault diagnosis results with good computation efficiency. Additionally, the embedded formal language of the neuron network can provide explanations about the decision process, thus achieve interpretable fault diagnosis.

preprint2022arXiv

Proposal for the search for exotic spin-spin interactions at the micrometer scale using functionalized cantilever force sensors

Spin-dependent exotic interactions can be generated by exchanging hypothetical bosons, which were introduced to solve some puzzles in physics. Many precision experiments have been performed to search for such interactions, but no confirmed observation has been made. Here, we propose new experiments to search for the exotic spin-spin interactions that can be mediated by axions or Z$^\prime$ bosons. A sensitive functionalized cantilever is utilized as a force sensor to measure the interactions between the spin-polarized electrons in a periodic magnetic source structure and a closed-loop magnetic structure integrated on the cantilever. The source is set to oscillate during data acquisition to modulate the exotic force signal to high harmonics of the oscillating frequency. This helps to suppress the spurious signals at the signal frequency. Different magnetic source structures are designed for different interaction detections. A magnetic stripe structure is designed for Z$^\prime$-mediated interaction, which is insensitive to the detection of axion-mediated interaction. This allows us to measure the coupling constant of both if we assume both exist. With the force sensitivity achievable at low temperature, the proposed experiments are expected to search for the parameter spaces with much smaller coupling constant than the current stringent constraints from micrometer to millimeter range. Specifically, the lower bound of the parameter space will be seven orders of magnitude lower than the stringent constraints for Z$^\prime$-mediated interaction, and an order of magnitude lower for axion-mediated interaction, at the interaction range of $10\, μ$m.

preprint2021arXiv

Development of a GPU-accelerated Monte Carlo dose calculation module for nuclear medicine, ARCHER-NM: Demonstration for a PET/CT imaging procedure

This paper describes the development and validation of a Monte Carlo (MC) dose computing module dedicated to organ dose calculations of patients undergoing nuclear medicine (NM) internal radiation exposures involving 18F-FDG PET/CT examination. This new module extends the more-than-10-years-long ARCHER project that developed a GPU-accelerated MC dose engine by adding dedicated NM source-definition features. To validate the code, we compared dose distributions from the 0.511-MeV point photon source calculated for a water phantom as well as a patient PET/CT phantom against a well-tested MC code, GATE. The water-phantom results show excellent agreement, suggesting that the radiation physics module in the new NM code is adequate. To demonstrate the clinical utility and advantage of ARCHER-NM, one set of PET/CT data for an adult male NM patient is calculated using the new code. Radiosensitive organs in the CT dataset are segmented using a CNN-based tool called DeepViewer. The PET image intensity maps are converted to radioactivity distributions to allow for MC radiation transport dose calculations at the voxel level. The dose rate maps and corresponding statistical uncertainties were calculated for the duration of PET image acquisition. The dose rate results of the 18F-FDG PET imaging patient show that ARCHER-NM's results agree very well with those of the GATE within 0.58% to 4.11%. Most impressively, ARCHER-NM obtains such results in less than 0.5 minutes while it takes GATE as much as 376 minutes. This is the first study presenting GPU-accelerated patient-specific MC internal radiation dose rate calculations for clinically realistic 18F-FDG PET/CT imaging cases involving auto-segmentation of whole-body PET/CT images. This study suggests that modern computing tools -- ARCHER-NM and DeepViewer -- are accurate and fast enough for routine internal dosimetry in NM clinics.

preprint2020arXiv

A random-walk model for dark matter halo spins

We extend the random-walk model of Vitvitska et al. for predicting the spins of dark matter halos from their merger histories. Using updated merger rates, orbital parameter distributions, and N-body constraints we show that this model can accurately reproduce the distribution of spin parameters measured in N-body simulations when we include a weak correlation between the spins of halos and the angular momenta of infalling subhalos. We further show that this model is in approximate agreement with the correlation of the spin magnitude over time as determined from N-body simulations, while it slightly underpredicts the correlation in the direction of the spin vector measured from the same simulations. This model is useful for predicting spins from merger histories derived from non-N-body sources, thereby circumventing the need for very high resolution simulations to permit accurate measurements of spins. It may be particularly relevant to modeling systems which accumulate angular momentum from halos over time (such as galactic discs) - we show that this model makes small but significant changes in the distribution of galactic disc sizes computed using the Galacticus semi-analytic galaxy formation model.

preprint2020arXiv

C-DLinkNet: considering multi-level semantic features for human parsing

Human parsing is an essential branch of semantic segmentation, which is a fine-grained semantic segmentation task to identify the constituent parts of human. The challenge of human parsing is to extract effective semantic features to resolve deformation and multi-scale variations. In this work, we proposed an end-to-end model called C-DLinkNet based on LinkNet, which contains a new module named Smooth Module to combine the multi-level features in Decoder part. C-DLinkNet is capable of producing competitive parsing performance compared with the state-of-the-art methods with smaller input sizes and no additional information, i.e., achiving mIoU=53.05 on the validation set of LIP dataset.

preprint2020arXiv

GINet: Graph Interaction Network for Scene Parsing

Recently, context reasoning using image regions beyond local convolution has shown great potential for scene parsing. In this work, we explore how to incorporate the linguistic knowledge to promote context reasoning over image regions by proposing a Graph Interaction unit (GI unit) and a Semantic Context Loss (SC-loss). The GI unit is capable of enhancing feature representations of convolution networks over high-level semantics and learning the semantic coherency adaptively to each sample. Specifically, the dataset-based linguistic knowledge is first incorporated in the GI unit to promote context reasoning over the visual graph, then the evolved representations of the visual graph are mapped to each local representation to enhance the discriminated capability for scene parsing. GI unit is further improved by the SC-loss to enhance the semantic representations over the exemplar-based semantic graph. We perform full ablation studies to demonstrate the effectiveness of each component in our approach. Particularly, the proposed GINet outperforms the state-of-the-art approaches on the popular benchmarks, including Pascal-Context and COCO Stuff.

preprint2020arXiv

Object Instance Mining for Weakly Supervised Object Detection

Weakly supervised object detection (WSOD) using only image-level annotations has attracted growing attention over the past few years. Existing approaches using multiple instance learning easily fall into local optima, because such mechanism tends to learn from the most discriminative object in an image for each category. Therefore, these methods suffer from missing object instances which degrade the performance of WSOD. To address this problem, this paper introduces an end-to-end object instance mining (OIM) framework for weakly supervised object detection. OIM attempts to detect all possible object instances existing in each image by introducing information propagation on the spatial and appearance graphs, without any additional annotations. During the iterative learning process, the less discriminative object instances from the same class can be gradually detected and utilized for training. In addition, we design an object instance reweighted loss to learn larger portion of each object instance to further improve the performance. The experimental results on two publicly available databases, VOC 2007 and 2012, demonstrate the efficacy of proposed approach.

preprint2020arXiv

Quark structure of the $χ_{\rm c}(3P)$ and $X(4274)$ resonances and their strong and radiative decays

We calculate the masses of $χ_{\rm c}(3P)$ states with threshold corrections in a coupled-channel model. The model was recently applied to the description of the properties of $χ_{\rm c}(2P)$ and $χ_{\rm b}(3P)$ multiplets [Phys.\ Lett.\ B {\bf 789}, 550 (2019)]. We also compute the open-charm strong decay widths of the $χ_{\rm c}(3P)$ states and their radiative transitions. According to our predictions, the $χ_{\rm c}(3P)$ states should be dominated by the charmonium core, but they may also show small meson-meson components. The $X(4274)$ is interpreted as a $c \bar c$ $χ_{\rm c1}(3P)$ state. More informations on the other members of the $χ_{\rm c}(3P)$ multiplet, as well as a more rigorous analysis of the $X(4274)$'s decay modes, are needed to provide further indications on the quark structure of the previous resonance.

preprint2020arXiv

Towards Interpretable Deep Learning Models for Knowledge Tracing

As an important technique for modeling the knowledge states of learners, the traditional knowledge tracing (KT) models have been widely used to support intelligent tutoring systems and MOOC platforms. Driven by the fast advancements of deep learning techniques, deep neural network has been recently adopted to design new KT models for achieving better prediction performance. However, the lack of interpretability of these models has painfully impeded their practical applications, as their outputs and working mechanisms suffer from the intransparent decision process and complex inner structures. We thus propose to adopt the post-hoc method to tackle the interpretability issue for deep learning based knowledge tracing (DLKT) models. Specifically, we focus on applying the layer-wise relevance propagation (LRP) method to interpret RNN-based DLKT model by backpropagating the relevance from the model's output layer to its input layer. The experiment results show the feasibility using the LRP method for interpreting the DLKT model's predictions, and partially validate the computed relevance scores from both question level and concept level. We believe it can be a solid step towards fully interpreting the DLKT models and promote their practical applications in the education domain.

preprint2020arXiv

Unveiling the secrets of the mid-infrared Moon

The Moon's optical characteristics in visible and long-wavelength infrared (LWIR) have long been observed with our eyes or with instruments. What the mid-infrared (MIR) Moon looks like is still a mystery. For the first time we present detailed appearance of the MIR Moon observed by a high-resolution geostationary satellite and reveal the essence behind its appearance. The appearance of the MIR Moon is opposite to its normal visible appearance. In addition the MIR Moon shows limb darkening. Both the absolute and the relative brightness distribution of the MIR lunar disk changes with the solar incidence angle. The signatures of the MIR Moon are controlled by both the reflection and emission of the lunar surface. We also show first-ever brightness temperature maps of the lunar disk without needing a mosaic, which better show the temperature variation across the lunar disk. They reveal that the relationship between brightness temperature and solar incidence angle i is cos1/bi, and the power parameter is smaller than the Lambertian temperature model of cos1/4i observed for lunar orbit-based measurements. The slower decrease of the brightness temperature when moving away from the sub-solar point than the Lambertian model is due to topographic effects. The brightness temperature is dominated by albedo and the solar incidence angle and influenced by the topography. Our results indicate that the Moon in the MIR exhibits many interesting phenomena which were previously unknown, and contains abundant information about lunar reflection and thermal emission for future study.

preprint2016arXiv

Coupled-Channel Effects for the Bottomonium with Realistic Wave Functions

With Gaussian expansion method (GEM), realistic wave functions are used to calculate coupled-channel effects for the bottomonium under the framework of ${}^3P_0$ model. The simplicity and accuracy of GEM are explained. We calculate the mass shifts, probabilities of the $B$ meson continuum, $S-D$ mixing angles, strong and dielectric decay widths. Our calculation shows that both $S-D$ mixing and the $B$ meson continuum can contribute to the suppression of the vector meson's dielectric decay width. We suggest more precise measurements on the radiative decays of $Υ(10580)$ and $Υ(11020)$ to distinguish these two effects. The above quantities are also calculated with simple harmonic oscillator (SHO) wave function approximation for comparison. The deviation between GEM and SHO indicates that it is essential to treat the wave functions accurately for near threshold states.

preprint2016arXiv

Exact Exponent in Optimal Rates for Crowdsourcing

In many machine learning applications, crowdsourcing has become the primary means for label collection. In this paper, we study the optimal error rate for aggregating labels provided by a set of non-expert workers. Under the classic Dawid-Skene model, we establish matching upper and lower bounds with an exact exponent $mI(π)$ in which $m$ is the number of workers and $I(π)$ the average Chernoff information that characterizes the workers' collective ability. Such an exact characterization of the error exponent allows us to state a precise sample size requirement $m>\frac{1}{I(π)}\log\frac{1}ε$ in order to achieve an $ε$ misclassification error. In addition, our results imply the optimality of various EM algorithms for crowdsourcing initialized by consistent estimators.

preprint2016arXiv

Stellar Mass--Gas-phase Metallicity Relation at $0.5\leq z\leq0.7$: A Power Law with Increasing Scatter toward the Low-mass Regime

We present the stellar mass ($M_{*}$)--gas-phase metallicity relation (MZR) and its scatter at intermediate redshifts ($0.5\leq z\leq0.7$) for 1381 field galaxies collected from deep spectroscopic surveys. The star formation rate (SFR) and color at a given $M_{*}$ of this magnitude-limited ($R\lesssim24$ AB) sample are representative of normal star-forming galaxies. For masses below $10^9 M_\odot$, our sample of 237 galaxies is $\sim$10 times larger than those in previous studies beyond the local universe. This huge gain in sample size enables superior constraints on the MZR and its scatter in the low-mass regime. We find a power-law MZR at $10^{8} M_\odot < M_{*} < 10^{11} M_\odot$: ${12+log(O/H) = (5.83\pm0.19) + (0.30\pm0.02)log(M_{*}/M_\odot)}$. Our MZR shows good agreement with others measured at similar redshifts in the literature in the intermediate and massive regimes, but is shallower than the extrapolation of the MZRs of others to masses below $10^{9} M_\odot$. The SFR dependence of the MZR in our sample is weaker than that found for local galaxies (known as the Fundamental Metallicity Relation). Compared to a variety of theoretical models, the slope of our MZR for low-mass galaxies agrees well with predictions incorporating supernova energy-driven winds. Being robust against currently uncertain metallicity calibrations, the scatter of the MZR serves as a powerful diagnostic of the stochastic history of gas accretion, gas recycling, and star formation of low-mass galaxies. Our major result is that the scatter of our MZR increases as $M_{*}$ decreases. Our result implies that either the scatter of the baryonic accretion rate or the scatter of the $M_{*}$--$M_{halo}$ relation increases as $M_{*}$ decreases. Moreover, our measures of scatter at $z=0.7$ appears consistent with that found for local galaxies.

preprint2016arXiv

The connection between the host halo and the satellite galaxies of the Milky Way

Many properties of the Milky Way's dark matter halo, including its mass assembly history, concentration, and subhalo population, remain poorly constrained. We explore the connection between these properties of the Milky Way and its satellite galaxy population, especially the implication of the presence of the Magellanic Clouds for the properties of the Milky Way halo. Using a suite of high-resolution $N$-body simulations of Milky Way-mass halos with a fixed final Mvir ~ 10^{12.1}Msun, we find that the presence of Magellanic Cloud-like satellites strongly correlates with the assembly history, concentration, and subhalo population of the host halo, such that Milky Way-mass systems with Magellanic Clouds have lower concentration, more rapid recent accretion, and more massive subhalos than typical halos of the same mass. Using a flexible semi-analytic galaxy formation model that is tuned to reproduce the stellar mass function of the classical dwarf galaxies of the Milky Way with Markov-Chain Monte-Carlo, we show that adopting host halos with different mass-assembly histories and concentrations can lead to different best-fit models for galaxy-formation physics, especially for the strength of feedback. These biases arise because the presence of the Magellanic Clouds boosts the overall population of high-mass subhalos, thus requiring a different stellar-mass-to-halo-mass ratio to match the data. These biases also lead to significant differences in the mass--metallicity relation, the kinematics of low-mass satellites, the number counts of small satellites associated with the Magellanic Clouds, and the stellar mass of Milky Way itself. Observations of these galaxy properties can thus provide useful constraints on the properties of the Milky Way halo.

preprint2016arXiv

The Evolution of the Galaxy Stellar Mass Function at z= 4-8: A Steepening Low-mass-end Slope with Increasing Redshift

We present galaxy stellar mass functions (GSMFs) at $z=$ 4-8 from a rest-frame ultraviolet (UV) selected sample of $\sim$4500 galaxies, found via photometric redshifts over an area of $\sim$280 arcmin$^2$ in the CANDELS/GOODS fields and the Hubble Ultra Deep Field. The deepest Spitzer/IRAC data yet-to-date and the relatively large volume allow us to place a better constraint at both the low- and high-mass ends of the GSMFs compared to previous space-based studies from pre-CANDELS observations. Supplemented by a stacking analysis, we find a linear correlation between the rest-frame UV absolute magnitude at 1500 Å ($M_{\rm UV}$) and logarithmic stellar mass ($\log M_*$) that holds for galaxies with $\log(M_*/M_{\odot}) \lesssim 10$. We use simulations to validate our method of measuring the slope of the $\log M_*$-$M_{\rm UV}$ relation, finding that the bias is minimized with a hybrid technique combining photometry of individual bright galaxies with stacked photometry for faint galaxies. The resultant measured slopes do not significantly evolve over $z=$ 4-8, while the normalization of the trend exhibits a weak evolution toward lower masses at higher redshift. We combine the $\log M_*$-$M_{\rm UV}$ distribution with observed rest-frame UV luminosity functions at each redshift to derive the GSMFs, finding that the low-mass-end slope becomes steeper with increasing redshift from $α=-1.55^{+0.08}_{-0.07}$ at $z=4$ to $α=-2.25^{+0.72}_{-0.35}$ at $z=8$. The inferred stellar mass density, when integrated over $M_*=10^8$-$10^{13} M_{\odot}$, increases by a factor of $10^{+30}_{-2}$ between $z=7$ and $z=4$ and is in good agreement with the time integral of the cosmic star formation rate density.

preprint2015arXiv

A Critical Assessment of Stellar Mass Measurement Methods

In this paper we perform a comprehensive study of the main sources of random and systematic errors in stellar mass measurement for galaxies using their Spectral Energy Distributions (SEDs). We use mock galaxy catalogs with simulated multi-waveband photometry (from U-band to mid-infrared) and known redshift, stellar mass, age and extinction for individual galaxies. Given different parameters affecting stellar mass measurement (photometric S/N ratios, SED fitting errors, systematic effects, the inherent degeneracies and correlated errors), we formulated different simulated galaxy catalogs to quantify these effects individually. We studied the sensitivity of stellar mass estimates to the codes/methods used, population synthesis models, star formation histories, nebular emission line contributions, photometric uncertainties, extinction and age. For each simulated galaxy, the difference between the input stellar masses and those estimated using different simulation catalogs, $Δ\log(M)$, was calculated and used to identify the most fundamental parameters affecting stellar masses. We measured different components of the error budget, with the results listed as follows: (1). no significant bias was found among different codes/methods, with all having comparable scatter; (2). A source of error is found to be due to photometric uncertainties and low resolution in age and extinction grids; (3). The median of stellar masses among different methods provides a stable measure of the mass associated with any given galaxy; (4). The deviations in stellar mass strongly correlate with those in age, with a weaker correlation with extinction; (5). the scatter in the stellar masses due to free parameters are quantified, with the sensitivity of the stellar mass to both the population synthesis codes and inclusion of nebular emission lines studied.

preprint2015arXiv

An analytical model for galaxy metallicity: What do metallicity relations tell us about star formation and outflow?

We develop a simple analytical model that tracks galactic metallicities governed by star formation and feedback to gain insight from the observed galaxy stellar mass-metallicity relations over a large range of stellar masses and redshifts. The model reveals the following implications of star formation and feedback processes in galaxy formation. First, the observed metallicity relations provide a stringent upper limit for the averaged outflow mass-loading factors of local galaxies, which is ~20 for M_*~10^9Msun galaxies and monotonically decreases to ~1 for M_*~10^{11}Msun galaxies. Second, the inferred upper-limit for the outflow mass-loading factor sensitively depends on whether the outflow is metal-enriched with respect to the ISM metallicity. If half of the metals ejected from SNe leave the galaxy in metal-enriched winds, the outflow mass-loading factor for galaxies at any mass can barely be higher than ~10, which puts strong constraints on galaxy formation models. Third, the relatively lower stellar-phase to gas-phase metallicity ratio for lower-mass galaxies indicate that low-mass galaxies are still rapidly enriching their metallicities in recent times, while high-mass galaxies are more settled, which seems to show a downsizing effect in the metallicity evolution of galaxies. The analysis presented in the paper demonstrates the importance of accurate measurements of galaxy metallicities and the cold gas fraction of galaxies at different redshifts for constraining star formation and feedback processes, and demonstrates the power of these relations in constraining the physics of galaxy formation.

preprint2015arXiv

Clumpy Galaxies in CANDELS. I. The Definition of UV Clumps and the Fraction of Clumpy Galaxies at 0.5<z<3

Although giant clumps of stars are crucial to galaxy formation and evolution, the most basic demographics of clumps are still uncertain, mainly because the definition of clumps has not been thoroughly discussed. In this paper, we study the basic demographics of clumps in star-forming galaxies (SFGs) at 0.5<z<3, using our proposed physical definition that UV-bright clumps are discrete star-forming regions that individually contribute more than 8% of the rest-frame UV light of their galaxies. Clumps defined this way are significantly brighter than the HII regions of nearby large spiral galaxies, either individually or blended, when physical spatial resolution and cosmological dimming are considered. Under this definition, we measure the fraction of SFGs that contain at least one off-center clump (Fclumpy) and the contributions of clumps to the rest-frame UV light and star formation rate of SFGs in the CANDELS/GOODS-S and UDS fields, where our mass-complete sample consists of 3239 galaxies with axial ratio q>0.5. The redshift evolution of Fclumpy changes with the stellar mass (M*) of the galaxies. Low-mass (log(M*/Msun)<9.8) galaxies keep an almost constant Fclumpy of about 60% from z~3.0 to z~0.5. Intermediate-mass and massive galaxies drop their Fclumpy from 55% at z~3.0 to 40% and 15%, respectively, at z~0.5. We find that (1) the trend of disk stabilization predicted by violent disk instability matches the Fclumpy trend of massive galaxies; (2) minor mergers are a viable explanation of the Fclumpy trend of intermediate-mass galaxies at z<1.5, given a realistic observability timescale; and (3) major mergers are unlikely responsible for the Fclumpy trend in all masses at z<1.5. The clump contribution to the rest-frame UV light of SFGs shows a broad peak around galaxies with log(M*/Msun)~10.5 at all redshifts, possibly linked to the molecular gas fraction of the galaxies. (Abridged)

preprint2015arXiv

Estimating stellar atmospheric parameters based on LASSO and support-vector regression

A scheme for estimating atmospheric parameters T$_{eff}$, log$~g$, and [Fe/H] is proposed on the basis of Least Absolute Shrinkage and Selection Operator (LASSO) algorithm and Haar wavelet. The proposed scheme consists of three processes. A spectrum is decomposed using the Haar wavelet transform and low-frequency components at the fourth level are considered as candidate features. Then, spectral features from the candidate features are detected using the LASSO algorithm to estimate the atmospheric parameters. Finally, atmospheric parameters are estimated from the extracted spectral features using the support-vector regression (SVR) method. The proposed scheme was evaluated using three sets of stellar spectra respectively from Sloan Digital Sky Survey (SDSS), Large Sky Area Multi-object Fiber Spectroscopic Telescope (LAMOST), and Kurucz's model, respectively. The mean absolute errors are as follows: for 40~000 SDSS spectra, 0.0062 dex for log~T$_{eff}$ (85.83 K for T$_{eff}$), 0.2035 dex for log$~g$ and 0.1512 dex for [Fe/H]; for 23963 LAMOST spectra, 0.0074 dex for log~T$_{eff}$ (95.37 K for T$_{eff}$), 0.1528 dex for log~$g$, and 0.1146 dex for [Fe/H]; and for 10469 synthetic spectra, 0.0010 dex for log T$_{eff}$(14.42K for T$_{eff}$), 0.0123 dex for log~$g$, and 0.0125 dex for [Fe/H].

preprint2015arXiv

Galaxy Ecosystems: gas contents, inflows and outflows

We use a set of observational data for galaxy cold gas mass fraction and gas phase metallicity to constrain the content, inflow and outflow of gas in central galaxies hosted by halos with masses between $10^{11} M_{\odot}$ to $10^{12} M_{\odot}$. The gas contents in high redshift galaxies are obtained by combining the empirical star formation histories of Lu et al. (2014) and star formation models that relate star formation rate with the cold gas mass in galaxies. We find that the total baryon mass in low-mass galaxies is always much less than the universal baryon mass fraction since $z = 2$, regardless of star formation model adopted. The data for the evolution of the gas phase metallicity require net metal outflow at $z\lesssim 2$, and the metal loading factor is constrained to be about $0.01$, or about $60\%$ of the metal yield. Based on the assumption that galactic outflow is more enriched in metal than both the interstellar medium and the material ejected at earlier epochs, we are able to put stringent constraints on the upper limits for both the net accretion rate and the net mass outflow rate. The upper limits strongly suggest that the evolution of the gas phase metallicity and gas mass fraction for low-mass galaxies at $z < 2$ is not compatible with strong outflow. We speculate that the low star formation efficiency of low-mass galaxies is owing to some preventative processes that prevent gas from accreting into galaxies in the first place.

preprint2015arXiv

Gas and Metal Contents of Galaxies and Gaseous Halos: Preventive versus Ejective Feedback

Using a semi-analytical approach we investigate the characteristics of predictions for the masses and metallicities of the baryonic matter in and around galaxies made by three galaxy formation models. These models represent three different feedback scenarios: one model with purely ejective feedback, one model with ejective feedback with reincorporation of ejected gas, and one preventative model. We find that, when the model parameters are adjusted to predict the correct stellar masses for a range of halo masses between 10^{10} to 10^{12}Msun, these three scenarios have very different predictions for the masses and metallicities of the interstellar and circum-galactic media. Compared with current observational data, the model implementing preventative feedback has a large freedom to match a broad range of observational data, while the ejective models have difficulties to match a number of observational constraints simultaneously, independent of how the ejection and reincorporation are implemented. Our results suggest that the feedback process which regulates the amounts of stars and cold gas in low-mass galaxies is preventative in nature.

preprint2015arXiv

Rate-optimal graphon estimation

Network analysis is becoming one of the most active research areas in statistics. Significant advances have been made recently on developing theories, methodologies and algorithms for analyzing networks. However, there has been little fundamental study on optimal estimation. In this paper, we establish optimal rate of convergence for graphon estimation. For the stochastic block model with $k$ clusters, we show that the optimal rate under the mean squared error is $n^{-1}\log k+k^2/n^2$. The minimax upper bound improves the existing results in literature through a technique of solving a quadratic equation. When $k\leq\sqrt{n\log n}$, as the number of the cluster $k$ grows, the minimax rate grows slowly with only a logarithmic order $n^{-1}\log k$. A key step to establish the lower bound is to construct a novel subset of the parameter space and then apply Fano's lemma, from which we see a clear distinction of the nonparametric graphon estimation problem from classical nonparametric regression, due to the lack of identifiability of the order of nodes in exchangeable random graph models. As an immediate application, we consider nonparametric graphon estimation in a Hölder class with smoothness $α$. When the smoothness $α\geq1$, the optimal rate of convergence is $n^{-1}\log n$, independent of $α$, while for $α\in(0,1)$, the rate is $n^{-2α/(α+1)}$, which is, to our surprise, identical to the classical nonparametric rate.

preprint2015arXiv

Satellite Quenching and Galactic Conformity at 0.3 < z < 2.5

We measure the evolution of the quiescent fraction and quenching efficiency of satellites around star-forming and quiescent central galaxies with stellar mass $\log(M_{\mathrm{cen}}/M_{\odot})>10.5$ at $0.3<z<2.5$. We combine imaging from three deep near-infrared-selected surveys (ZFOURGE/CANDELS, UDS, and UltraVISTA), which allows us to select a stellar-mass complete sample of satellites with $\log(M_{\mathrm{sat}}/M_{\odot})>9.3$. Satellites for both star-forming and quiescent central galaxies have higher quiescent fractions compared to field galaxies matched in stellar mass at all redshifts. We also observe "galactic conformity": satellites around quiescent centrals are more likely to be quenched compared to the satellites around star-forming centrals. In our sample, this conformity signal is significant at $\gtrsim3σ$ for $0.6<z<1.6$, whereas it is only weakly significant at $0.3<z<0.6$ and $1.6<z<2.5$. Therefore, conformity (and therefore satellite quenching) has been present for a significant fraction of the age of the universe. The satellite quenching efficiency increases with increasing stellar mass of the central, but does not appear to depend on the stellar mass of the satellite to the mass limit of our sample. When we compare the satellite quenching efficiency of star-forming centrals with stellar masses 0.2 dex higher than quiescent centrals (which should account for any difference in halo mass), the conformity signal decreases, but remains statistically significant at $0.6<z<0.9$. This is evidence that satellite quenching is connected to the star-formation properties of the central as well as to the mass of the halo. We discuss physical effects that may contribute to galactic conformity, and emphasize that they must allow for continued star-formation in the central galaxy even as the satellites are quenched.

preprint2015arXiv

SDSS/SEGUE Spectral Feature Analysis For Stellar Atmospheric Parameter Estimation

Large-scale and deep sky survey missions are rapidly collecting a large amount of stellar spectra, which necessitate the estimation of atmospheric parameters directly from spectra and makes it feasible to statistically investigate latent principles in a large dataset. We present a technique for estimating parameters $T_{eff}$, log$~g$ and [Fe/H] from stellar spectra. With this technique, we first extract features from stellar spectra using the LASSO algorithm; then, the parameters are estimated from the extracted features using the SVR. On a subsample of 20~000 stellar spectra from SDSS with reference parameters provided by SDSS/SEGUE Pipeline SSPP, estimation consistency are 0.007458 dex for log$~T_{eff}$ (101.609921 K for $T_{eff}$), 0.189557 dex for log$~g$ and 0.182060 for [Fe/H], where the consistency is evaluated by mean absolute error. Prominent characteristics of the proposed scheme are sparseness, locality, and physical interpretability. In this work, every spectrum consists of 3821 fluxes, and 10, 19, and 14 typical wavelength positions are detected respectively for estimating $T_{eff}$, log$~g$ and [Fe/H]. It is shown that the positions are related to typical lines of stellar spectra. This characteristic is important in investigating physical indications from analysis results. Then, stellar spectra can be described by the individual fluxes on the detected positions (PD) or local integration of fluxes near them (LI). The abovementioned consistency is the result based on features described by LI. If features are described by PD, consistency are 0.009092 dex for log$~T_{eff}$ (124.545075 K for $T_{eff}$), 0.198928 dex for log$~g$, and 0.206814 dex for [Fe/H].

preprint2015arXiv

Star Formation and Stellar Mass Assembly in Dark Matter Halos: From Giants to Dwarfs

The empirical model of Lu et al. 2014 is updated with recent data and used to study galaxy star formation and assembly histories. At $z > 2$, the predicted galaxy stellar mass functions are steep, and a significant amount of star formation is hosted by low-mass haloes that may be missed in current observations. Most of the stars in cluster centrals formed earlier than $z\approx 2$ but have been assembled much later. Milky Way mass galaxies have had on-going star formation without significant mergers since $z\approx 2$, and are thus free of significant (classic) bulges produced by major mergers. In massive clusters, stars bound in galaxies and scattered in the halo form a homogeneous population that is old and with solar metallicity. In contrast, in Milky Way mass systems the two components form two distinct populations, with halo stars being older and poorer in metals by a factor of $\approx 3$. Dwarf galaxies in haloes with $M_{\rm h} < 10^{11}h^{-1}M_{\odot}$ have experienced a star formation burst accompanied by major mergers at $z > 2$, followed by a nearly constant star formation rate after $z = 1$. The early burst leaves a significant old stellar population that is distributed in spheroids.

preprint2015arXiv

The Relation Between SFR and Stellar Mass for Galaxies at 3.5 $\le z\le$ 6.5 in CANDELS

Distant star-forming galaxies show a correlation between their star formation rates (SFR) and stellar masses, and this has deep implications for galaxy formation. Here, we present a study on the evolution of the slope and scatter of the SFR-stellar mass relation for galaxies at $3.5\leq z\leq 6.5$ using multi-wavelength photometry in GOODS-S from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) and Spitzer Extended Deep Survey. We describe an updated, Bayesian spectral-energy distribution fitting method that incorporates effects of nebular line emission, star formation histories that are constant or rising with time, and different dust attenuation prescriptions (starburst and Small Magellanic Cloud). From $z$=6.5 to $z$=3.5 star-forming galaxies in CANDELS follow a nearly unevolving correlation between stellar mass and SFR that follows SFR $\sim$ $M_\star^a$ with $a = 0.54 \pm 0.16$ at $z\sim 6$ and $0.70 \pm 0.21$ at $z\sim 4$. This evolution requires a star formation history that increases with decreasing redshift (on average, the SFRs of individual galaxies rise with time). The observed scatter in the SFR-stellar mass relation is tight, $σ(\log \mathrm{SFR}/\mathrm{M}_\odot$ yr$^{-1})< 0.3\ - $ 0.4 dex, for galaxies with $\log M_\star/\mathrm{M}_\odot > 9$ dex. Assuming that the SFR is tied to the net gas inflow rate (SFR $\sim$ $\dot{M}_\mathrm{gas}$), then the scatter in the gas inflow rate is also smaller than 0.3$-$0.4 dex for star-forming galaxies in these stellar mass and redshift ranges, at least when averaged over the timescale of star formation. We further show that the implied star formation history of objects selected on the basis of their co-moving number densities is consistent with the evolution in the SFR-stellar mass relation.

preprint2015arXiv

Using Galaxy Pairs to Probe Star Formation During Major Halo Mergers

Currently-proposed galaxy quenching mechanisms predict very different behaviours during major halo mergers, ranging from significant quenching enhancement (e.g., clump-induced gravitational heating models) to significant star formation enhancement (e.g., gas starvation models). To test real galaxies' behaviour, we present an observational galaxy pair method for selecting galaxies whose host haloes are preferentially undergoing major mergers. Applying the method to central L* (10^10 Msun < M_* < 10^10.5 Msun) galaxies in the Sloan Digital Sky Survey (SDSS) at z<0.06, we find that major halo mergers can at most modestly reduce the star-forming fraction, from 59% to 47%. Consistent with past research, however, mergers accompany enhanced specific star formation rates for star-forming L* centrals: ~10% when a paired galaxy is within 200 kpc (approximately the host halo's virial radius), climbing to ~70% when a paired galaxy is within 30 kpc. No evidence is seen for even extremely close pairs (<30 kpc separation) rejuvenating star formation in quenched galaxies. For galaxy formation models, our results suggest: (1) quenching in L* galaxies likely begins due to decoupling of the galaxy from existing hot and cold gas reservoirs, rather than a lack of available gas or gravitational heating from infalling clumps, (2) state-of-the-art semi-analytic models currently over-predict the effect of major halo mergers on quenching, and (3) major halo mergers can trigger enhanced star formation in non-quenched central galaxies.

preprint2014arXiv

An Empirical Model for the Star Formation History in Dark Matter Halos

We develop an empirical approach to infer the star formation rate in dark matter halos from the galaxy stellar mass function (SMF) at different redshifts and the local cluster galaxy luminosity function (CGLF), which has a steeper faint end relative to the SMF of local galaxies. As satellites are typically old galaxies which have been accreted earlier, this feature can cast important constraint on the formation of low-mass galaxies at high-redshift. The evolution of the SMFs suggests the star formation in high mass halos ($>10^{12}M_{\odot}h^{-1}$) has to be boosted at high redshift beyond what is expected from a simple scaling of the dynamical time. The faint end of the CGLF implies a characteristic redshift $z_c\approx2$ above which the star formation rate in low mass halos with masses $< 10^{11}M_{\odot}h^{-1}$ must be enhanced relative to that at lower z. This is not directly expected from the standard stellar feedback models. Also, this enhancement leads to some interesting predictions, for instance, a significant old stellar population in present-day dwarf galaxies with $M_* < 10^8 M_{\odot}h^{-2}$ and steep slopes of high redshift stellar mass and star formation rate functions.

preprint2014arXiv

Bayesian inferences of galaxy formation from the K-band luminosity and HI mass functions of galaxies: constraining star formation and feedback

We infer mechanisms of galaxy formation for a broad family of semi-analytic models (SAMs) constrained by the K-band luminosity function and HI mass function of local galaxies using tools of Bayesian analysis. Even with a broad search in parameter space the whole model family fails to match to constraining data. In the best fitting models, the star formation and feedback parameters in low-mass haloes are tightly constrained by the two data sets, and the analysis reveals several generic failures of models that similarly apply to other existing SAMs. First, based on the assumption that baryon accretion follows the dark matter accretion, large mass-loading factors are required for haloes with circular velocities lower than 200 km/s, and most of the wind mass must be expelled from the haloes. Second, assuming that the feedback is powered by Type-II supernovae with a Chabrier IMF, the outflow requires more than 25% of the available SN kinetic energy. Finally, the posterior predictive distributions for the star formation history are dramatically inconsistent with observations for masses similar to or smaller than the Milky-Way mass. The inferences suggest that the current model family is still missing some key physical processes that regulate the gas accretion and star formation in galaxies with masses below that of the Milky Way.

preprint2014arXiv

De-Confusing blended field images using graphs and bayesian priors

We present a new technique for overcoming confusion noise in deep far-infrared \Herschel space telescope images making use of prior information from shorter $λ<2$\micron wavelengths. For the deepest images obtained by \Herschels, the flux limit due to source confusion is about a factor of three brighter than the flux limit due to instrumental noise and (smooth) sky background. We have investigated the possibility of de-confusing simulated \Herschel PACS-160\micron images by using strong Bayesian priors on the positions and weak priors on the flux of sources. We find the blended sources and group them together and simultaneously fit their fluxes. We derive the posterior probability distribution function of fluxes subject to these priors through Monte Carlo Markov Chain (MCMC) sampling by fitting the image. Assuming we can predict FIR flux of sources based on ultraviolet-optical part of their SEDs to within an order of magnitude, the simulations show that we can obtain reliable fluxes and uncertainties at least a factor of three fainter than the confusion noise limit of $3σ_{c} $=2.7 mJy in our simulated PACS-160 image. This technique could in principle be used to mitigate the effects of source confusion in any situation where one has prior information of positions and plausible fluxes of blended sources. For \Herschel, application of this technique will improve our ability to constrain the dust content in normal galaxies at high redshift.

preprint2014arXiv

Formal Specification and Quantitative Analysis of a Constellation of Navigation Satellites

Navigation satellites are a core component of navigation satellite based systems such as GPS, GLONASS and Galileo which provide location and timing information for a variety of uses. Such satellites are designed for operating on orbit to perform tasks and have lifetimes of 10 years or more. Reliability, availability and maintainability (RAM) analysis of systems has been indispensable in the design phase of satellites in order to achieve minimum failures or to increase mean time between failures (MTBF) and thus to plan maintenance strategies, optimise reliability and maximise availability. In this paper, we present formal models of both a single satellite and a navigation satellite constellation and logical specification of their reliability, availability and maintainability properties respectively. The probabilistic model checker PRISM has been used to perform automated analysis of these quantitative properties.

preprint2014arXiv

Formation of disk galaxies in preheated media: a preventative feedback model

We introduce a semi-analytic galaxy formation model implementing a self-consistent treatment for the hot halo gas configuration and the assembly of central disks. Using the model, we explore a preventative feedback model, in which the circum-halo medium is assumed to be preheated up to a certain entropy level by early starbursts or other processes, and compare it with an ejective feedback model, in which baryons are first accreted into dark matter halos and subsequently ejected out by feedback. The model demonstrates that when the medium is preheated to an entropy comparable to the halo virial entropy the baryon accretion can be largely reduced and delayed. In addition, the preheated medium can establish an extended low density gaseous halo when it accretes into the dark matter halos, and result in a specific angular momentum of the cooling gas large enough to form central disks as extended as those observed. Combined with simulated halo assembly histories, the preventative feedback model can reproduce remarkably well a number of observational scaling relations. These include the cold baryon (stellar plus cold gas) mass fraction-halo mass relations, star formation histories, disk size-stellar mass relation and its evolution, and the number density of low-mass galaxies as a function of redshift. In contrast, the conventional ejective feedback model fails to reproduce these observational trends. Using the model, we demonstrate that the properties of disk galaxies are closely tied to the thermal state of hot halo gas and even possibly the circum-halo medium, which suggests that observational data for the disk properties and circum-galactic hot/warm medium may jointly provide interesting constraints for galaxy formation models.

preprint2014arXiv

Individualized Rank Aggregation using Nuclear Norm Regularization

In recent years rank aggregation has received significant attention from the machine learning community. The goal of such a problem is to combine the (partially revealed) preferences over objects of a large population into a single, relatively consistent ordering of those objects. However, in many cases, we might not want a single ranking and instead opt for individual rankings. We study a version of the problem known as collaborative ranking. In this problem we assume that individual users provide us with pairwise preferences (for example purchasing one item over another). From those preferences we wish to obtain rankings on items that the users have not had an opportunity to explore. The results here have a very interesting connection to the standard matrix completion problem. We provide a theoretical justification for a nuclear norm regularized optimization procedure, and provide high-dimensional scaling results that show how the error in estimating user preferences behaves as the number of observations increase.

preprint2014arXiv

Properties of Submillimeter Galaxies in the CANDELS GOODS-S Field

We derive physical properties of 10 submillimeter galaxies located in the CANDELS coverage of the GOODS-S field. The galaxies were first identified as submillimeter sources with the LABOCA bolometer and subsequently targeted for 870um continuum observation with ALMA. The high angular resolution of the ALMA imaging allows secure counterparts to be identified in the CANDELS multiband dataset. The CANDELS data provide deep photometric data from UV through near-infrared wavelengths. Using synthetic spectral energy distributions, we derive photometric redshifts, stellar masses, extinction, ages, and the star formation history. The redshift range is z=1.65-4.76, with two of the galaxies located at z>4. Two SMG counterparts have stellar masses 2-3 orders of magnitude lower than the rest. The remaining SMG counterparts have stellar masses around 1x10^11 Msun. The stellar population in the SMGs is typically older than the expected duration of the submillimeter phase, suggesting that the star formation history of submillimeter galaxies is more complex than a single burst. Non-parametric morphology indices suggest that the SMG counterparts are among the most asymmetric systems compared with galaxies of the same stellar mass and redshift. The HST images shows that 3 of the SMGs are associated with on-going mergers. The remaining counterparts are isolated. Estimating the dust and molecular gas mass from the submm fluxes, and comparing with our stellar masses shows that the molecular gas mass fraction of SMGs is ~28% and that the final stellar mass is likely to be (1-2)x10^11 Msun.

preprint2014arXiv

The distribution of satellites around massive galaxies at 1<z<3 in ZFOURGE/CANDELS: dependence on star formation activity

We study the statistical distribution of satellites around star-forming and quiescent central galaxies at 1<z<3 using imaging from the FourStar Galaxy Evolution Survey (ZFOURGE) and the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). The deep near-IR data select satellites down to $\log(M/M_\odot)>9$ at z<3. The radial satellite distribution around centrals is consistent with a projected NFW profile. Massive quiescent centrals, $\log(M/M_\odot)>10.78$, have $\sim$2 times the number of satellites compared to star-forming centrals with a significance of 2.7$σ$ even after accounting for differences in the centrals' stellar-mass distributions. We find no statistical difference in the satellite distributions of intermediate-mass quiescent and star-forming centrals, $10.48<\log(M/M_\odot)<10.78$. Comparing to the Guo2011 semi-analytic model, the excess number of satellites indicates that quiescent centrals have halo masses 0.3 dex larger than star-forming centrals, even when the stellar-mass distributions are fixed. We use a simple toy model that relates halo mass and quenching, which roughly reproduces the observed quenched fractions and the differences in halo mass between star-forming and quenched galaxies only if galaxies have a quenching probability that increases with halo mass from $\sim$0 for $\log(M_h/M_\odot)\sim$11 to $\sim$1 for $\log(M_h/M_\odot)\sim$13.5. A single halo-mass quenching threshold is unable to reproduce the quiescent fraction and satellite distribution of centrals. Therefore, while halo quenching may be an important mechanism, it is unlikely to be the only factor driving quenching. It remains unclear why a high fraction of centrals remain star-forming even in relatively massive halos.

preprint2014arXiv

The mass evolution of the first galaxies: stellar mass functions and star formation rates at $4 < z < 7$ in the CANDELS GOODS-South field

We measure new estimates for the galaxy stellar mass function and star formation rates for samples of galaxies at $z \sim 4,~5,~6~\&~7$ using data in the CANDELS GOODS South field. The deep near-infrared observations allow us to construct the stellar mass function at $z \geq 6$ directly for the first time. We estimate stellar masses for our sample by fitting the observed spectral energy distributions with synthetic stellar populations, including nebular line and continuum emission. The observed UV luminosity functions for the samples are consistent with previous observations, however we find that the observed $M_{UV}$ - M$_{*}$ relation has a shallow slope more consistent with a constant mass to light ratio and a normalisation which evolves with redshift. Our stellar mass functions have steep low-mass slopes ($α\approx -1.9$), steeper than previously observed at these redshifts and closer to that of the UV luminosity function. Integrating our new mass functions, we find the observed stellar mass density evolves from $\log_{10} ρ_{*} = 6.64^{+0.58}_{-0.89}$ at $z \sim 7$ to $7.36\pm0.06$ $\text{M}_{\odot} \text{Mpc}^{-3}$ at $z \sim 4$. Finally, combining the measured UV continuum slopes ($β$) with their rest-frame UV luminosities, we calculate dust corrected star-formation rates (SFR) for our sample. We find the specific star-formation rate for a fixed stellar mass increases with redshift whilst the global SFR density falls rapidly over this period. Our new SFR density estimates are higher than previously observed at this redshift.

preprint2013arXiv

A Comparison between Semi-Analytic Model Predictions for the CANDELS Survey

We compare the predictions of three independently developed semi-analytic galaxy formation models that are being used to aid in the interpretation of results from the CANDELS survey. These models are each applied to the same set of halo merger trees extracted from the "Bolshoi" simulation and are carefully tuned to match the local galaxy stellar mass function using the powerful method of Bayesian Inference coupled with MCMC or by hand. The comparisons reveal that in spite of the significantly different parameterizations for star formation and feedback processes, the three models yield qualitatively similar predictions for the assembly histories of galaxy stellar mass and star formation over cosmic time. We show that the SAMs generally require strong outflows to suppress star formation in low-mass halos to match the present day stellar mass function. However, all of the models considered produce predictions for the star formation rates and metallicities of low-mass galaxies that are inconsistent with existing data and diverge between the models. We suggest that large differences in the metallicity relations and small differences in the stellar mass assembly histories of model galaxies stem from different assumptions for the outflow mass-loading factor. Importantly, while more accurate observational measurements for stellar mass, SFR and metallicity of galaxies at 1<z<5 will discriminate between models, the discrepancies between the models and existing data of these observables have already revealed challenging problems in understanding star formation and its feedback in galaxy formation. The three sets of models are being used to construct catalogs of mock galaxies on light cones that have the same geometry as the CANDELS survey, which should be particularly useful for quantifying the biases and uncertainties on measurements and inferences from the real observations. -ABRIDGED

preprint2013arXiv

Constraining The Assembly Of Normal And Compact Passively Evolving Galaxies From Redshift z=3 To The Present With CANDELS

We study the evolution of the number density, as a function of the size, of passive early-type galaxies with a wide range of stellar masses 10^10<M*/Msun<10^11.5) from z~3 to z~1, exploiting the unique dataset available in the GOODS-South field, including the recently obtained WFC3 images as a part of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). In particular, we select a sample of 107 massive (M*>10^10 M_sun), passive (SSFR<10^-2 Gyr^-1) and morphologically spheroidal galaxies at 1.2<z<3, taking advantage of the panchromatic dataset available for GOODS, including VLT, CFHT, Spitzer, Chandra and HST ACS+WFC3 data. We find that at 1<z<3 the passively evolving early-type galaxies are the reddest and most massive objects in the Universe, and we prove that a correlation between mass, morphology, color and star-formation activity is already in place at that epoch. We measure a significant evolution in the mass-size relation of passive early-type galaxies (ETGs) from z~3 to z~1, with galaxies growing on average by a factor of 2 in size in a 3 Gyr timescale only. We witness also an increase in the number density of passive ETGs of 50 times over the same time interval. We find that the first ETGs to form at z>2 are all compact or ultra-compact, while normal sized ETGs (meaning ETGs with sizes comparable to those of local counterparts of the same mass) are the most common ETGs only at z<1. The increase of the average size of ETGs at 0<z<1 is primarily driven by the appearance of new large ETGs rather than by the size increase of individual galaxies.

preprint2013arXiv

Constraining the Star Formation Histories in Dark Matter Halos: I. Central Galaxies

Using the self-consistent modeling of the conditional stellar mass functions across cosmic time by Yang et al. (2012), we make model predictions for the star formation histories (SFHs) of {\it central} galaxies in halos of different masses. The model requires the following two key ingredients: (i) mass assembly histories of central and satellite galaxies, and (ii) local observational constraints of the star formation rates of central galaxies as function of halo mass. We obtain a universal fitting formula that describes the (median) SFH of central galaxies as function of halo mass, galaxy stellar mass and redshift. We use this model to make predictions for various aspects of the star formation rates of central galaxies across cosmic time. Our main findings are the following. (1) The specific star formation rate (SSFR) at high $z$ increases rapidly with increasing redshift [$\propto (1+z)^{2.5}$] for halos of a given mass and only slowly with halo mass ($\propto M_h^{0.12}$) at a given $z$, in almost perfect agreement with the specific mass accretion rate of dark matter halos. (2) The ratio between the star formation rate (SFR) in the main-branch progenitor and the final stellar mass of a galaxy peaks roughly at a constant value, $\sim 10^{-9.3} h^2 {\rm yr}^{-1}$, independent of halo mass or the final stellar mass of the galaxy. However, the redshift at which the SFR peaks increases rapidly with halo mass. (3) More than half of the stars in the present-day Universe were formed in halos with $10^{11.1}\msunh < M_h < 10^{12.3}\msunh$ in the redshift range $0.4 < z < 1.9$. (4) ... [abridged]

preprint2013arXiv

Cosmological Constraints from Galaxy Clustering and the Mass-to-Number Ratio of Galaxy Clusters: Marginalizing over the Physics of Galaxy Formation

Many approaches to obtaining cosmological constraints rely on the connection between galaxies and dark matter. However, the distribution of galaxies is dependent on their formation and evolution as well as the cosmological model, and galaxy formation is still not a well-constrained process. Thus, methods that probe cosmology using galaxies as a tracer for dark matter must be able to accurately estimate the cosmological parameters without knowing the details of galaxy formation a priori. We apply this reasoning to the method of obtaining $Ω_m$ and $σ_8$ from galaxy clustering combined with the mass-to-number ratio of galaxy clusters. To test the sensitivity of this method to variations due to galaxy formation, we consider several different models applied to the same cosmological dark matter simulation. The cosmological parameters are then estimated using the observables in each model, marginalizing over the parameters of the Halo Occupation Distribution (HOD). We find that for models where the galaxies can be well represented by a parameterized HOD, this method can successfully extract the desired cosmological parameters for a wide range of galaxy formation prescriptions.

preprint2013arXiv

Mergers and Mass Accretion for Infalling Halos Both End Well Outside Cluster Virial Radii

We find that infalling dark matter halos (i.e., the progenitors of satellite halos) begin losing mass well outside the virial radius of their eventual host halos. The peak mass occurs at a range of clustercentric distances, with median and 68th percentile range of 1.8 +2.3/-1.0 R_(vir,host) for progenitors of z=0 satellites. The peak circular velocity for infalling halos occurs at significantly larger distances (3.7 +3.3/-2.2 R_(vir,host) at z=0). This difference arises because different physical processes set peak circular velocity (typically, ~1:5 and larger mergers which cause transient circular velocity spikes) and peak mass (typically, smooth accretion) for infalling halos. We find that infalling halos also stop having significant mergers well before they enter the virial radius of their eventual hosts. Mergers larger than a 1:40 ratio in halo mass end for infalling halos at similar clustercentric distances (~ 1.9 R_(vir,host)) as the end of overall mass accretion. However, mergers larger than 1:3 typically end for infalling halos at more than 4 virial radial away from their eventual hosts. This limits the ability of mergers to affect quenching and morphology changes in clusters. We also note that the transient spikes which set peak circular velocity may lead to issues with abundance matching on that parameter, including unphysical galaxy stellar mass growth profiles near clusters; we propose a simple observational test to check if a better halo proxy for galaxy stellar mass exists.

preprint2013arXiv

Preliminary Application of Reinsch Splines to Cosmology: Transition Redshift Determination with Simulated OHD

Many schemes have been proposed to perform a model-independent constraint on cosmological dynamics, such as nonparametric dark energy equation of state (EoS) ω(z) or the deceleration parameter q(z). These methods usually contain derivative processes with respect to observational data with noise. However, it still remains remarkably uncertain when one estimates the numerical differentiation, especially the corresponding truncation errors. In this work, we introduce a global numerical differentiation method, first formulated by Reinsch(1967), which is smoothed by cubic spline functions. The optimal solution is obtained by minimizing the functional Φ(f). To investigate the potential of the algorithm further, we apply it to the estimation of the transition redshift z_{t} with simulated expansion rate E(z) based on observational Hubble parameter data(OHD). An effective method to determine the free parameter S appearing in Reinsch Splines is provided.

preprint2013arXiv

Submm/mm Galaxy Counterpart Identification Using a Characteristic Density Distribution

We present a new submm/mm galaxy counterpart identification technique which builds on the use of Spitzer IRAC colors as discriminators between likely counterparts and the general IRAC galaxy population. Using 102 radio- and SMA-confirmed counterparts to AzTEC sources across three fields (GOODS-N, GOODS-S, and COSMOS), we develop a non-parametric IRAC color-color characteristic density distribution (CDD), which, when combined with positional uncertainty information via likelihood ratios, allows us to rank all potential IRAC counterparts around SMGs and calculate the significance of each ranking via the reliability factor. We report all robust and tentative radio counterparts to SMGs, the first such list available for AzTEC/COSMOS, as well as the highest ranked IRAC counterparts for all AzTEC SMGs in these fields as determined by our technique. We demonstrate that the technique is free of radio bias and thus applicable regardless of radio detections. For observations made with a moderate beamsize (~18"), this technique identifies ~85 per cent of SMG counterparts. For much larger beamsizes (>30"), we report identification rates of 33-49 per cent. Using simulations, we demonstrate that this technique is an improvement over using positional information alone for observations with facilities such as AzTEC on the LMT and SCUBA-2 on JCMT.

preprint2012arXiv

Bayesian inference of galaxy formation from the K-band luminosity function of galaxies: tensions between theory and observation

We conduct Bayesian model inferences from the observed K-band luminosity function of galaxies in the local Universe, using the semi-analytic model (SAM) of galaxy formation introduced in Lu et al (2011). The prior distributions for the 14 free parameters include a large range of possible models. We find that some of the free parameters, e.g. the characteristic scales for quenching star formation in both high-mass and low-mass halos, are already tightly constrained by the single data set. The posterior distribution includes the model parameters adopted in other SAMs. By marginalising over the posterior distribution, we make predictions that include the full inferential uncertainties for the colour-magnitude relation, the Tully-Fisher relation, the conditional stellar mass function of galaxies in halos of different masses, the HI mass function, the redshift evolution of the stellar mass function of galaxies, and the global star formation history. Using posterior predictive checking with the available observational results, we find that the model family (i) predicts a Tully-Fisher relation that is curved; (ii) significantly over predicts the satellite fraction; (iii) vastly over predicts the HI mass function; (iv) predicts high-z stellar mass functions that have too many low mass galaxies and too few high mass ones. and (v) predicts a redshift evolution of the stellar mass density and the star formation history that are in moderate disagreement. These results suggest that some important processes are still missing in the current model family and we discuss a number of possible solutions to solve the discrepancies, such as interactions between galaxies and dark matter halos, tidal stripping, the bimodal accretion of gas, preheating, and a redshift-dependent initial mass function.

preprint2011arXiv

A Bayesian approach to the semi-analytic model of galaxy formation: methodology

We believe that a wide range of physical processes conspire to shape the observed galaxy population but we remain unsure of their detailed interactions. The semi-analytic model (SAM) of galaxy formation uses multi-dimensional parameterisations of the physical processes of galaxy formation and provides a tool to constrain these underlying physical interactions. Because of the high dimensionality, the parametric problem of galaxy formation may be profitably tackled with a Bayesian-inference based approach, which allows one to constrain theory with data in a statistically rigorous way. In this paper we develop a SAM in the framework of Bayesian inference. We show that, with a parallel implementation of an advanced Markov-Chain Monte-Carlo algorithm, it is now possible to rigorously sample the posterior distribution of the high-dimensional parameter space of typical SAMs. As an example, we characterise galaxy formation in the current $Λ$CDM cosmology using the stellar mass function of galaxies as an observational constraint. We find that the posterior probability distribution is both topologically complex and degenerate in some important model parameters, suggesting that thorough explorations of the parameter space are needed to understand the models. We also demonstrate that because of the model degeneracy, adopting a narrow prior strongly restricts the model. Therefore, the inferences based on SAMs are conditional to the model adopted. Using synthetic data to mimic systematic errors in the stellar mass function, we demonstrate that an accurate observational error model is essential to meaningful inference.

preprint2011arXiv

On the algorithms of radiative cooling in semi-analytic models

We study the behaviour of multiple radiative cooling algorithms implemented in seven Semi-Analytic Models (SAMs) of galaxy formation, including a new model we propose in this paper. We use versions of the models without feedback and apply them to dark matter haloes growing in a cosmological context, which have final masses that range from 10^{11}Msun to 10^{14}Msun. First, using simplified smoothly-growing halo models, we demonstrate that the different algorithms predict cooling rates and final cold gas masses that differ by a factor of ~5 for massive haloes (>10^{12}Msun). The algorithms are in better agreement for less massive haloes because they cool efficiently and, therefore, their cooling rates are largely limited by the halo accretion rate. However, for less massive haloes, all the SAMs predict less cooling than corresponding 1D hydrodynamic models. Second, we study the gas accretion history of the central galaxies of dark matter haloes using merger trees. The inclusion of mergers alters the cooling history of haloes by locking up gas in galaxies within small haloes at early times. For realistic halo models, the dispersion in the cold gas mass predicted by the algorithms is 0.5 dex for high mass haloes and 0.1 dex for low mass haloes, while the dispersion in the accretion rate is about two times larger. Comparing to cosmological SPH simulations, we find that most SAMs systematically under-predict the gas accretion rates for low-mass haloes but over-predict the gas accretion rates for massive haloes. Although the models all include both "rapid" and "slow" mode accretion, the transition between the two accretion modes varies between models and also differs from the simulations. Finally, we construct a new model that explicitly incorporates cold halo gas to illustrate that such a class of models can better match the results from cosmological hydrodynamic simulations.

preprint2011arXiv

The 2009 outburst from the new X-ray transient and black-hole candidate XTE J1652$-$453

The RXTE and Swift observations on the 2009 outburst from a newly discovered transient and black-hole candidate XTE J1652$-$453 are analyzed. The source was observed by RXTE to behave a sequence of spectral states that are typical to the black hole XRBs. During the first 7 observations, the source is diagnosed to stay in a high/soft state, from the spectrum dominated by soft thermal component which contributes an average of $\sim$ 85% to the X-ray flux at 2-20keV, and from the hardness $\sim$ 0.1 showing up in the hardness-intensity diagram (HID). For the last 20 observations the spectral state is classified as {\sl low/hard} according to an average hardness of $\sim$ 0.8 and a balance between the thermal and the non-thermal components: a power-law component takes ${\geq}80%$ of total 2-20keV flux. Located in between is an {\sl intermediate} state that the source might have experienced. The usual relationship between rms and hardness presents in HRD as well. Throughout the outburst no QPOs are found in XTE J1652$-$453.

preprint2010arXiv

On the selection effect of radio quasars in the Sloan Digital Sky Survey

We identified a large sample of radio quasars, including those with complex radio morphology, from the Sloan Digital Sky Survey (SDSS) and the Faint Images of Radio Sky at Twenty-cm (FIRST). Using this sample, we inspect previous radio quasar samples for selection effects resulting from complex radio morphologies and adopting positional coincidence between radio and optical sources alone. We find that 13.0% and 8.1% radio quasars do not show a radio core within 1.2 and 2 arcsecs of their optical position, and thus are missed in such samples. Radio flux is under-estimated by a factor of more than 2 for an additional 8.7% radio quasars. These missing radio extended quasars are more radio loud with a typical radio-to-optical flux ratio namely radio loudness RL >100, and radio power P >10^{25} W/Hz. They account for more than one third of all quasars with RL>100. The color of radio extended quasars tends to be bluer than the radio compact quasars. This suggests that radio extended quasars are more radio powerful sources, e.g., Fanaroff-Riley type 2 (FR-II) sources, rather than the compact ones viewed at larger inclination angles. By comparison with the radio data from the NRAO VLA Sky Survey (NVSS), we find that for sources with total radio flux less than 3 mJy, low surface brightness components tend to be underestimated by FIRST, indicating that lobes in these faint radio sources are still missed.

preprint2010arXiv

Type 1 Active Galactic Nucleus Fraction in SDSS/FIRST Survey

In the unification scheme, narrow-lined (type 2) active galactic nuclei (AGN) are intrinsically similar to broad-lined (type 1) AGN with the exception that the line of sight to the broad emission line region and accretion disk is blocked by a dusty torus. The fraction of type 1 AGN measures the average covering factor of the torus. In this paper, we explore the dependence of this fraction on nuclear properties for a sample of low redshift (z <0.35) radio strong (P_{1.4GHz} >10^{23}W/Hz) AGN selected by matching the spectroscopic catalog of Sloan Digital Sky Survey and the radio source catalog of Faint Image of Radio Sky at Twenty cm. After correcting for several selection effects, we find that : (1) type 1 fraction $f_1$ keeps at a constant of ~20 per cent in the [O III] 5007 luminosity range of 40.7< log(L_{[O III]}/ erg/s) <43.5 . This result is significantly different from previous studies, and the difference can be explained by extinction correction and different treatment of selection effects. (2) $f_1$ rises with black hole mass from ~20 per cent (M_bh below 10^8 Msun) to ~30 per cent (M_bh above that). This coincides with the decrease of the fraction of highly-inclined disk galaxies with black hole mass, implying a population of Seyfert galaxies seen as type 2 due to galaxy-scale obscuration in disk when the host galaxy type transfer from bulge-dominant to disk-dominant. (3) $f_1$ is independent of the Eddington ratio for its value between 0.01 and 1; (4) $f_1$ ascends from 15 per cent to 30 per cent in the radio power range of 23< log(P_{1.4GHz}/ W/Hz) <24, then remain a constant at ~30 per cent up to 10^{26} W/Hz.

preprint2009arXiv

Phase-resolved spectral analysis of 4U 1901+03 during its outburst

The high mass X-ray binary 4U 1901+03 was reported to have the pulse profile evolving with the X-ray luminosity and energy during its outburst in February-July 2003: the pulse peak changed from double to single along with the decreasing luminosity. We have carried out a detailed analysis on the contemporary phase-resolved energy spectrum of 4U 1901+03 as observed by Rossi X-ray Timing Explorer (RXTE). We find that, both the continuum and the pulse spectra are phase dependent. The optical depth derived from the pulse spectrum is in general larger than that from the continuum. Fe Ka emission line is only detected in the spectrum of the continuum and is missing in the pulse spectrum. This suggests an origin of Fe emission from the accretion disk but not the surface of the neutron star.

preprint2008arXiv

Feedback from galactic stellar bulges and hot gaseous haloes of galaxies

We demonstrate that the feedback from stellar bulges can play an essential role in shaping the halo gas of galaxies with substantial bulge components by conducting 1-D hydrodynamical simulations. The feedback model we consider consists of two distinct phases: 1) an early starburst during the bulge formation and 2) a subsequent long-lasting mass and energy injection from stellar winds of low-mass stars and Type Ia SNe. An energetic outward blastwave is initiated by the starburst and is maintained and enhanced by the long-lasting stellar feedback. For a MW-like galactic bulge, this blastwave sweeps up the halo gas in the proto-galaxy and heats up the surrounding medium to a scale much beyond the virial radius of the halo, thus the accretion of the halo hot gas can be completely stopped. In addition, the long-lasting feedback in the later phase powers a galactic bulge wind that is reverse-shocked at a large radius in the presence of surrounding intergalactic medium and hence maintains a hot gaseous halo. As the mass and energy injection decreases with time, the feedback evolves to a subsonic and quasi-stable outflow, which is enough to prevent halo gas from cooling. The two phases of the feedback thus re-enforce each-other's impact on the gas dynamics. The simulation results demonstrate that the stellar bulge feedback may provide a plausible solution to the long-standing problems in understanding the MW type galaxies, such as the "missing stellar feedback" problem and the "over-cooling" problem. The simulations also show that the properties of the hot gas in the subsonic outflow state depend sensitively on the environment and the formation history of the bulge. This dependence and variance may explain the large dispersion in the X-ray to B-band luminosity ratio of the low $L_X/L_B$ Es.