Source author record

Huiyuan Wang

Huiyuan Wang 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

48works
7topics
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

48 published item(s)

preprint2026arXiv

Foundation Models to Unlock Real-World Evidence from Nationwide Medical Claims

Evidence derived from large-scale real-world data (RWD) is increasingly informing regulatory evaluation and healthcare decision-making. Administrative claims provide population-scale, longitudinal records of healthcare utilization, expenditure, and detailed coding of diagnoses, procedures, and medications, yet their potential as a substrate for healthcare foundation models remains largely unexplored. Here we present ReClaim, a generative transformer trained from scratch on 43.8 billion medical events from more than 200 million enrollees in the MarketScan claims data spanning 2008-2022. ReClaim models longitudinal trajectories across diagnoses, procedures, medications, and expenditure, and was scaled to 140 million, 700 million, and 1.7 billion parameters. Across over 1,000 disease-onset prediction tasks, ReClaim achieved a mean AUC of 75.6%, substantially outperforming disease-specific LightGBM (66.3%) and the transformer-based Delphi model (69.4%), with the largest gains for rare diseases. These advantages held across retrospective and prospective evaluations and in external validation on two independent datasets. Performance improved monotonically with scale, and post-training added 13.8 percentage points over pre-training alone. Beyond disease prediction, ReClaim captured financial outcomes and improved real-world evidence (RWE) analyses: for healthcare expenditure forecasting it increased explained variance from 0.28 to 0.37 relative to LightGBM, and in a target trial emulation it reduced systematic bias by 72% on average relative to Delphi. Together, these results establish administrative claims as a scalable substrate for healthcare foundation models and show that learned representations generalize across time periods and data sources, supporting disease surveillance, expenditure forecasting, and RWE generation.

preprint2026arXiv

Radio AGN feedback sustains quiescence only in a minority of massive galaxies

Radio active galactic nuclei (AGNs) eject a huge amount of energy into the surrounding medium and are thought to potentially prevent gas cooling and maintain the quiescence of massive galaxies. The short-lived, sporadic, and anisotropic nature of radio activities, coupled with the detection of abundant cold gas around some massive quiescent galaxies, raise questions about the efficiency of radio feedback in massive galaxies. Here we present an innovative method rooted in artificial intelligence to separate galaxies in which radio feedback is effective (RFE), regardless of current radio emission, from those in which radio feedback is ineffective (RFI), according to their optical images. Galaxies categorized as RFE are all dynamically hot, whereas quiescent RFI (RFI-Q) galaxies usually have extended cold-disk components. At given stellar mass, dark matter halos hosting RFE galaxies are between four to ten times more massive than those of RFI-Q galaxies. We find, for the first time, that almost all RFE galaxies have scant cold gas, irrespective of AGN activity. In contrast, many RFI-Q galaxies are surrounded by substantial amounts of condensed atomic gas, indicating a different evolutionary path from RFE galaxies. Our finding provides direct and compelling evidence that a radio AGN has gone through about 300 on-off cycles and that radio feedback can prevent gas cooling over a timescale much longer than that of radio activity. Contrary to general belief, our analysis shows that only a small fraction of massive galaxies are influenced by strong radio AGNs, suggesting that current galaxy formation models need serious revision.

preprint2022arXiv

An Extended Halo-based Group/Cluster finder: application to the DESI legacy imaging surveys DR8

We extend the halo-based group finder developed by \citet[][]{Yang2005a} to use data {\it simultaneously} with either photometric or spectroscopic redshifts. A mock galaxy redshift survey constructed from a high-resolution N-body simulation is used to evaluate the performance of this extended group finder. For galaxies with magnitude ${\rm z\le 21}$ and redshift $0<z\le 1.0$ in the DESI legacy imaging surveys (the Legacy Surveys), our group finder successfully identifies more than 60\% of the members in about $90\%$ of halos with mass $\ga 10^{12.5}\msunh$. Detected groups with mass $\ga 10^{12.0}\msunh$ have a purity (the fraction of true groups) greater than 90\%. The halo mass assigned to each group has an uncertainty of about 0.2 dex at the high mass end $\ga 10^{13.5}\msunh$ and 0.40 dex at the low mass end. Groups with more than 10 members have a redshift accuracy of $\sim 0.008$. We apply this group finder to the Legacy Surveys DR8 and find 5.2 Million groups with at least 3 members. About 387,000 of these groups have at least 10 members. The resulting catalog containing 3D coordinates, richness, halo masses, and total group luminosities, is made publicly available.

preprint2022arXiv

ELUCID VII: Using Constrained Hydro Simulations to Explore the Gas Component of the Cosmic Web

Using reconstructed initial conditions in the SDSS survey volume, we carry out constrained hydrodynamic simulations in three regions representing different types of the cosmic web: the Coma cluster of galaxies; the SDSS great wall; and a large low-density region at $z\sim 0.05$. These simulations, which include star formation and stellar feedback but no AGN formation and feedback, are used to investigate the properties and evolution of intergalactic and intra-cluster media. About half of the warm-hot intergalactic gas is associated with filaments in the local cosmic web. Gas in the outskirts of massive filaments and halos can be heated significantly by accretion shocks generated by mergers of filaments and halos, respectively, and there is a tight correlation between gas temperature and the strength of the local tidal field. The simulations also predict some discontinuities associated with shock fronts and contact edges, which can be tested using observations of the thermal SZ effect and X-rays. A large fraction of the sky is covered by Ly$α$ and OVI absorption systems, and most of the OVI systems and low-column density HI systems are associated with filaments in the cosmic web. The constrained simulations, which follow the formation and heating history of the observed cosmic web, provide an important avenue to interpret observational data. With full information about the origin and location of the cosmic gas to be observed, such simulations can also be used to develop observational strategies.

preprint2022arXiv

Elucidating Galaxy Assembly Bias in SDSS

We investigate the level of galaxy assembly bias in the Sloan Digital Sky Survey (SDSS) main galaxy sample using ELUCID, a state-of-the-art constrained simulation that accurately reconstructed the initial density perturbations within the SDSS volume. On top of the ELUCID haloes, we develop an extended HOD model that includes the assembly bias of central and satellite galaxies, parameterized as $\mathcal{Q}_\mathrm{cen}$ and $\mathcal{Q}_\mathrm{sat}$, respectively, to predict a suite of one- and two-point observables. In particular, our fiducial constraint employs the probability distribution of the galaxy number counts measured on $8\,\mathrm{Mpc}\,h^{-1}$ scales $N_8^g$ and the projected cross-correlation functions of quintiles of galaxies selected by $N_8^g$ with our entire galaxy sample. We perform extensive tests of the efficacy of our method by fitting the same observables to mock data using both constrained and non-constrained simulations. We discover that in many cases the level of cosmic variance between the two simulations can produce biased constraints that lead to an erroneous detection of galaxy assembly bias if the non-constrained simulation is used. When applying our method to the SDSS data, the ELUCID reconstruction effectively removes an otherwise strong degeneracy between cosmic variance and galaxy assembly bias in SDSS, enabling us to derive an accurate and stringent constraint on the latter. Our fiducial ELUCID constraint, for galaxies above a stellar mass threshold $M_*{=}10^{10.2}\,h^{-2}\,M_\odot$, is $\mathcal{Q}_\mathrm{cen}{=}{-}0.09\pm{0.05}$ and $\mathcal{Q}_\mathrm{sat}{=}0.09\pm{0.10}$, indicating no evidence for a significant~($>2σ$) galaxy assembly bias in the local Universe probed by SDSS. Finally, our method provides a promising path to the robust modelling of the galaxy-halo connection within future surveys like DESI and PFS.

preprint2022arXiv

Evidence for quasar fast outflows being accelerated at the scale of tens of parsecs

Quasar outflows may play a crucial role in regulating the host galaxy, although the spatial scale of quasar outflows remain a major enigma, with their acceleration mechanism poorly understood. The kinematic information of outflow is the key to understanding its origin and acceleration mechanism. Here, we report the galactocentric distances of different outflow components for both a sample and an individual quasar. We find that the outflow distance increases with velocity, with a typical value from several parsecs to more than one hundred parsecs, providing direct evidence for an acceleration happening at a scale of the order of 10 parsecs. These outflows carry ~1% of the total quasar energy, while their kinematics are consistent with a dust driven model with a launching radius comparable to the scale of a dusty torus, indicating that the coupling between dust and quasar radiation may produce powerful feedback that is crucial to galaxy evolution.

preprint2022arXiv

Massive Star-Forming Galaxies Have Converted Most of Their Halo Gas into Stars

In the local Universe, the efficiency for converting baryonic gas into stars is very low. In dark matter halos where galaxies form and evolve, the average efficiency varies with galaxy stellar mass and has a maximum of about twenty percent for Milky-Way-like galaxies. The low efficiency at higher mass is believed to be produced by some quenching processes, such as the feedback from active galactic nuclei. We perform an analysis of weak lensing and satellite kinematics for SDSS central galaxies. Our results reveal that the efficiency is much higher, more than sixty percent, for a large population of massive star-forming galaxies around $10^{11}M_{\odot}$. This suggests that these galaxies acquired most of the gas in their halos and converted it into stars without being affected significantly by quenching processes. This population of galaxies is not reproduced in current galaxy formation models, indicating that our understanding of galaxy formation is incomplete. The implications of our results on circumgalactic media, star formation quenching and disc galaxy rotation curves are discussed. We also examine systematic uncertainties in halo-mass and stellar-mass measurements that might influence our results.

preprint2020arXiv

Detection of missing baryons in galaxy groups with kinetic Sunyaev-Zel'dovich effect

We present the detection of the kinetic Sunyaev-Zel'dovich effect (kSZE) signals from groups of galaxies as a function of halo mass down to $\log (M_{500}/{\rm M_\odot}) \sim 12.3$, using the {\it Planck} CMB maps and stacking about $40,000$ galaxy systems with known positions, halo masses, and peculiar velocities. The signals from groups of different mass are constrained simultaneously to take care of projection effects of nearby halos. The total kSZE flux within halos estimated implies that the gas fraction in halos is about the universal baryon fraction, even in low-mass halos, indicating that the `missing baryons' are found. Various tests performed show that our results are robust against systematic effects, such as contamination by infrared/radio sources and background variations, beam-size effects and contributions from halo exteriors. Combined with the thermal Sunyaev-Zel'dovich effect, our results indicate that the `missing baryons' associated with galaxy groups are contained in warm-hot media with temperatures between $10^5$ and $10^6\,{\rm K}$.

preprint2020arXiv

Detection of missing baryons in galaxy groups with kinetic Sunyaev-Zel'dovich effect

We present the detection of the kinetic Sunyaev-Zel'dovich effect (kSZE) signals from groups of galaxies as a function of halo mass down to $\log (M_{500}/{\rm M_\odot}) \sim 12.3$, using the {\it Planck} CMB maps and stacking about $40,000$ galaxy systems with known positions, halo masses, and peculiar velocities. The signals from groups of different mass are constrained simultaneously to take care of projection effects of nearby halos. The total kSZE flux within halos estimated implies that the gas fraction in halos is about the universal baryon fraction, even in low-mass halos, indicating that the `missing baryons' are found. Various tests performed show that our results are robust against systematic effects, such as contamination by infrared/radio sources and background variations, beam-size effects and contributions from halo exteriors. Combined with the thermal Sunyaev-Zel'dovich effect, our results indicate that the `missing baryons' associated with galaxy groups are contained in warm-hot media with temperatures between $10^5$ and $10^6\,{\rm K}$.

preprint2020arXiv

Probing Primordial Chirality with Galaxy Spins

Chiral symmetry is maximally violated in weak interactions, and such microscopic asymmetries in the early Universe might leave observable imprints on astrophysical scales without violating the cosmological principle. In this Letter, we propose a helicity measurement to detect primordial chiral violation. We point out that observations of halo-galaxy angular momentum directions (spins), which are frozen in during the galaxy formation process, provide a fossil chiral observable. From the clustering mode of large scale structure of the Universe, we construct a spin mode in Lagrangian space and show in simulations that it is a good probe of halo-galaxy spins. In standard model, a strong symmetric correlation between the left and right helical components of this spin mode and galaxy spins is expected. Measurements of these correlations will be sensitive to chiral breaking, providing a direct test of chiral symmetry breaking in the early Universe.

preprint2020arXiv

Relating the structure of dark matter halos to their assembly and environment

We use a large $N$-body simulation to study the relation of the structural properties of dark matter halos to their assembly history and environment. The complexity of individual halo assembly histories can be well described by a small number of principal components (PCs), which, compared to formation times, provide a more complete description of halo assembly histories and have a stronger correlation with halo structural properties. Using decision trees built with the random ensemble method, we find that about $60\%$, $10\%$, and $20\%$ of the variances in halo concentration, axis ratio, and spin, respectively, can be explained by combining four dominating predictors: the first PC of the assembly history, halo mass, and two environment parameters. Halo concentration is dominated by halo assembly. The local environment is found to be important for the axis ratio and spin but is degenerate with halo assembly. The small percentages of the variance in the axis ratio and spin that are explained by known assembly and environmental factors suggest that the variance is produced by many nuanced factors and should be modeled as such. The relations between halo intrinsic properties and environment are weak compared to their variances, with the anisotropy of the local tidal field having the strongest correlation with halo properties. Our method of dimension reduction and regression can help simplify the characterization of the halo population and clarify the degeneracy among halo properties.

preprint2020arXiv

The Breakdown Scale of HI Bias Linearity

The 21 cm intensity mapping experiments promise to obtain the large-scale distribution of HI gas at the post-reionization epoch. In order to reveal the underlying matter density fluctuations from the HI mapping, it is important to understand how HI gas traces the matter density distribution. Both nonlinear halo clustering and nonlinear effects modulating HI gas in halos may determine the scale below which the HI bias deviates from linearity. We employ three approaches to generate the mock HI density from a large-scale N-body simulation at low redshifts, and demonstrate that the assumption of HI linearity is valid at the scale corresponding to the first peak of baryon acoustic oscillations, but breaks down at $k \gtrsim 0.1\,h\, {\rm Mpc}^{-1}$. The nonlinear effects of halo clustering and HI content modulation counteract each other at small scales, and their competition results in a model-dependent "sweet-spot" redshift near $z$=1 where the HI bias is scale-independent down to small scales. We also find that the linear HI bias scales approximately linearly with redshift for $z\le 3$.

preprint2020arXiv

The Formation History of Subhalos and the Evolution of Satellite Galaxies

Satellites constitute an important fraction of the overall galaxy population and are believed to form in dark matter subhalos. Here we use the cosmological hydrodynamic simulation TNG100 to investigate how the formation histories of subhalos affect the properties and evolution of their host galaxies. We use a scaled formation time ($a_{\rm nf}$) to characterize the mass assembly histories of the subhalos before they are accreted by massive host halos. We find that satellite galaxies in young subhalos (low $a_{\rm nf}$) are less massive and more gas rich, and have stronger star formation and a higher fraction of ex situ stellar mass than satellites in old subhalos (high $a_{\rm nf}$). Furthermore, these low $a_{\rm nf}$ satellites require longer timescales to be quenched as a population than the high $a_{\rm nf}$ counterparts. We find very different merger histories between satellites in fast accretion (FA, $a_{\rm nf}<1.3$) and slow accretion (SA, $a_{\rm nf}>1.3$) subhalos. For FA satellites, the galaxy merger frequency dramatically increases just after accretion, which enhances the star formation at accretion. While, for SA satellites, the mergers occur smoothly and continuously across the accretion time. Moreover, mergers with FA satellites happen mainly after accretion, while a contrary trend is found for SA satellites. Our results provide insight into the evolution and star formation quenching of the satellite population.

preprint2019arXiv

The Dearth of Difference between Central and Satellite Galaxies III. Environmental Dependence of Mass-Size and Mass-Structure Relations

As demonstrated in Paper I, the quenching properties of central and satellite galaxies are quite similar as long as both stellar mass and halo mass are controlled. Here we extend the analysis to the size and bulge-to-total light ratio (B/T) of galaxies. In general central galaxies have size-stellar mass and B/T-stellar mass relations different from satellites. However, the differences are eliminated when halo mass is controlled. We also study the dependence of size and B/T on halo-centric distance and find a transitional stellar mass (M$_{*,t}$) at given halo mass (M$_h$), which is about one fifth of the mass of the central galaxies in halos of mass M$_h$. The transitional stellar masses for size, B/T and quenched fraction are similar over the whole halo mass range, suggesting a connection between the quenching of star formation and the structural evolution of galaxies. Our analysis further suggests that the classification based on the transitional stellar mass is more fundamental than the central-satellite dichotomy, and provide a more reliable way to understand the environmental effects on galaxy properties. We compare the observational results with the hydro-dynamical simulation, EAGLE and the semi-analytic model, L-GALAXIES. The EAGLE simulation successfully reproduces the similarities of size for centrals and satellites and even M$_{*,t}$, while L-GALAXIES fails to recover the observational results.

preprint2016arXiv

Alignments of dark matter halos with large-scale tidal fields: mass and redshift dependence

Large scale tidal field estimated directly from the distribution of dark matter halos is used to investigate how halo shapes and spin vectors are aligned with the cosmic web. The major, intermediate and minor axes of halos are aligned with the corresponding tidal axes, and halo spin axes tend to be parallel with the intermediate axes and perpendicular to the major axes of tidal field. The strengths of these alignments generally increase with halo mass and redshift, but the dependencies are only through the peak height, ν. The scaling relations of the alignment strengths with the value of ν indicate that the alignment strengths remain roughly constant when the structures within which the halos reside are still in quasi-linear regime, but decreases as nonlinear evolution becomes more important. We also calculate the alignments in projection so that our results can be compared directly with observations. Finally, we investigate the alignments of tidal tensors on large scales, and use the results to understand alignments of halo pairs separated at various distances. Our results suggest coherent structure of the tidal field is the underlying reason for the alignments of halos and galaxies seen in numerical simulations and in observations.

preprint2016arXiv

An empirical model to form and evolve galaxies in dark matter halos

Based on the star formation histories (SFH) of galaxies in halos of different masses, we develop an empirical model to grow galaxies in dark mattet halos. This model has very few ingredients, any of which can be associated to observational data and thus be efficiently assessed. By applying this model to a very high resolution cosmological $N$-body simulation, we predict a number of galaxy properties that are a very good match to relevant observational data. Namely, for both centrals and satellites, the galaxy stellar mass function (SMF) up to redshift $z\simeq4$ and the conditional stellar mass functions (CSMF) in the local universe are in good agreement with observations. In addition, the 2-point correlation is well predicted in the different stellar mass ranges explored by our model. Furthermore, after applying stellar population synthesis models to our stellar composition as a function of redshift, we find that the luminosity functions in $^{0.1}u$, $^{0.1}g$, $^{0.1}r$, $^{0.1}i$ and $^{0.1}z$ bands agree quite well with the SDSS observational results down to an absolute magnitude at about -17.0. The SDSS conditional luminosity functions (CLF) itself is predicted well. Finally, the cold gas is derived from the star formation rate (SFR) to predict the HI gas mass within each mock galaxy. We find a remarkably good match to observed HI-to-stellar mass ratios. These features ensure that such galaxy/gas catalogs can be used to generate reliable mock redshift surveys.

preprint2016arXiv

Differences in Halo-Scale Environments between Type 1 and Type 2 AGNs at Low Redshift

Using low-redshift (z<0.09) samples of AGNs, normal galaxies and groups of galaxies selected from the Sloan Digital Sky Survey (SDSS), we study the environments of type 1 and type 2 AGNs both on small and large scales. Comparisons are made for galaxy samples matched in redshift, $r$-band luminosity, [OIII] luminosity, and also the position in groups (central or satellite). We find that type 2 AGNs and normal galaxies reside in similar environments. Type 1 and type 2 AGNs have similar clustering properties on large scales ($\gtrsim1$Mpc), but at scales smaller than 100 kpc, type 2s have significant more neighbors than type 1s ($3.09\pm0.69$ times more for central AGNs at $\lesssim30$kpc). These results suggest that type 1 and type 2 AGNs are hosted by halos of similar masses, as is also seen directly from the mass distributions of their host groups ($\sim10^{12}h^{-1} M_{\odot}$ for centrals and $\sim10^{13}h^{-1} M_{\odot}$ for satellites). Type~2s have significantly more satellites around them, and the distribution of their satellites is also more centrally concentrated. The host galaxies of both types of AGNs have similar optical properties, but their infrared colors are significantly different. Our results suggest that the simple unified model based solely on torus orientation is not sufficient, but that galaxy interactions in dark matter halos must have played an important role in the formation of the dust structure that obscures AGNs.

preprint2016arXiv

ELUCID - Exploring the Local Universe with reConstructed Initial Density field III: Constrained Simulation in the SDSS Volume

A method we developed recently for the reconstruction of the initial density field in the nearby Universe is applied to the Sloan Digital Sky Survey Data Release 7. A high-resolution N-body constrained simulation (CS) of the reconstructed initial condition, with $3072^3$ particles evolved in a 500 Mpc/h box, is carried out and analyzed in terms of the statistical properties of the final density field and its relation with the distribution of SDSS galaxies. We find that the statistical properties of the cosmic web and the halo populations are accurately reproduced in the CS. The galaxy density field is strongly correlated with the CS density field, with a bias that depend on both galaxy luminosity and color. Our further investigations show that the CS provides robust quantities describing the environments within which the observed galaxies and galaxy systems reside. Cosmic variance is greatly reduced in the CS so that the statistical uncertainties can be controlled effectively even for samples of small volumes.

preprint2016arXiv

Evidence for Fluorescent Fe II Emission from Extended Low Ionization Outflows in Obscured Quasars

Recent studies have shown that outflows in at least some broad absorption line (BAL) quasars are extended well beyond the putative dusty torus. Such outflows should be detectable in obscured quasars. We present four WISE selected infrared red quasars with very strong and peculiar ultraviolet Fe ii emission lines: strong UV Fe II UV arising from transitions to ground/low excitation levels, and very weak Fe II at wavelengths longer than 2800 Å. The spectra of these quasars display strong resonant emission lines, such as C IV, Al III and Mg II but sometimes, a lack of non-resonant lines such as C III], S III and He II. We interpret the Fe II lines as resonantly scattered light from the extended outflows that are viewed nearly edge-on, so that the accretion disk and broad line region are obscured by the dusty torus, while the extended outflows are not. We show that dust free gas exposed to strong radiation longward of 912 Å produces Fe II emission very similar to that observed. The gas is too cool to collisionally excite Fe II lines, accounting for the lack of optical emission. The spectral energy distribution from the UV to the mid-infrared can be modeled as emission from a clumpy dusty torus, with UV emission being reflected/scattered light either by the dusty torus or the outflow. Within this scenario, we estimate a minimum covering factor of the outflows from a few to 20% for the Fe II scattering region, suggesting that Fe II BAL quasars are at a special stage of quasar evolution.

preprint2016arXiv

Galaxy groups in the 2MASS Redshift Survey

A galaxy group catalog is constructed from the 2MASS Redshift Survey (2MRS) with the use of a halo-based group finder. The halo mass associated with a group is estimated using a `GAP' method based on the luminosity of the central galaxy and its gap with other member galaxies. Tests using mock samples shows that this method is reliable, particularly for poor systems containing only a few members. On average 80% of all the groups have completeness >0.8, and about 65% of the groups have zero contamination. Halo masses are estimated with a typical uncertainty $\sim 0.35\,{\rm dex}$. The application of the group finder to the 2MRS gives 29,904 groups from a total of 43,246 galaxies at $z \leq 0.08$, with 5,286 groups having two or more members. Some basic properties of this group catalog is presented, and comparisons are made with other groups catalogs in overlap regions. With a depth to $z\sim 0.08$ and uniformly covering about 91% of the whole sky, this group catalog provides a useful data base to study galaxies in the local cosmic web, and to reconstruct the mass distribution in the local Universe.

preprint2016arXiv

Mapping the real space distributions of galaxies in SDSS DR7: I. Two Point Correlation Functions

Using a method to correct redshift space distortion (RSD) for individual galaxies, we mapped the real space distributions of galaxies in the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7). We use an ensemble of mock catalogs to demonstrate the reliability of our method. Here as the first paper in a series, we mainly focus on the two point correlation function (2PCF) of galaxies. Overall the 2PCF measured in the reconstructed real space for galaxies brighter than $^{0.1}{\rm M}_r-5\log h=-19.0$ agrees with the direct measurement to an accuracy better than the measurement error due to cosmic variance, if the reconstruction uses the correct cosmology. Applying the method to the SDSS DR7, we construct a real space version of the main galaxy catalog, which contains 396,068 galaxies in the North Galactic Cap with redshifts in the range $0.01 \leq z \leq 0.12$. The Sloan Great Wall, the largest known structure in the nearby Universe, is not as dominant an over-dense structure as appears to be in redshift space. We measure the 2PCFs in reconstructed real space for galaxies of different luminosities and colors. All of them show clear deviations from single power-law forms, and reveal clear transitions from 1-halo to 2-halo terms. A comparison with the corresponding 2PCFs in redshift space nicely demonstrates how RSDs boost the clustering power on large scales (by about $40-50\%$ at scales $\sim 10 h^{-1}{\rm {Mpc}}$) and suppress it on small scales (by about $70-80\%$ at a scale of $0.3 h^{-1}{\rm {Mpc}}$).

preprint2016arXiv

SDSS-IV MaNGA: environmental dependence of stellar age and metallicity gradients in nearby galaxies

We present a study on the stellar age and metallicity distributions for 1105 galaxies using the STARLIGHT software on MaNGA integral field spectra. We derive age and metallicity gradients by fitting straight lines to the radial profiles, and explore their correlations with total stellar mass M*, NUV-r colour and environments, as identified by both the large scale structure (LSS) type and the local density. We find that the mean age and metallicity gradients are close to zero but slightly negative, which is consistent with the inside-out formation scenario. Within our sample, we find that both the age and metallicity gradients show weak or no correlation with either the LSS type or local density environment. In addition, we also study the environmental dependence of age and metallicity values at the effective radii. The age and metallicity values are highly correlated with M* and NUV-r and are also dependent on LSS type as well as local density. Low-mass galaxies tend to be younger and have lower metallicity in low-density environments while high-mass galaxies are less affected by environment.

preprint2016arXiv

Strong Ly alpha Emission in the Proximate Damped Ly alpha Absorption Trough toward the Quasar SDSS J095253.83+011422.0

SDSS J095253.83+011422.0 (SDSS J0952+0114) was reported by Hall et al. (2004) as an exotic quasar at $z_{\rm em}=3.020$. In contrast to prominent broad metal--line emissions with FWHM~9000 km/s, only a narrow Ly αemission line is present with FWHM~1000 km/s. The absence of broad Ly alpha emission line has been a mystery for more than a decade. In this paper, we demonstrate that this is due to dark Proximate Damped Ly alpha Absorption (PDLA) at $z_{\rm abs}=3.010$ by identifying associated Lyman absorption line series from the damped Ly beta up to Ly9, as well as the Lyman limit absorption edge. The PDLA cloud has a column density of $\log N_{\rm H\,I}({\rm cm}^{-2})=21.8\pm0.2$, a metallicity of [Zn/H]$>-1.0$, and a spatial extent exceeding the Narrow Emission Line Region (NELR) of the quasar. With a luminosity of $L_{{\rm Ly}α}\sim10^{45}$ erg s$^{-1}$, the residual Ly alpha emission superposed on the PDLA trough is of two orders of magnitude stronger than previous reports. This is best explained as re-radiated photons arising from the quasar outflowing gas at a scale larger than the NELR. The PDLA here, acting like a natural coronagraph, provides us with a good insight into the illuminated gas in the vicinity of the quasar, which are usually hard to resolve due to their small size and "seeing fuzz" of bright quasars. Notably, SDSS J0952+0114 analogs might be easily omitted in the spectroscopic surveys of DLAs and PDLAs, as their damped Ly alpha troughs can be fully filled by additional strong Ly alpha emissions. Our preliminary survey shows that such systems are not very rare. They are potentially a unique sample for probing strong quasar feedback phenomena in the early universe.

preprint2015arXiv

A Comprehensive Study of Broad Absorption Line Quasars: I. Prevalence of HeI* Absorption Line Multiplets in Low-Ionization Objects

Neutral Helium multiplets, HeI*3189,3889,10830 are very useful diagnostics to the geometry and physical conditions of the absorbing gas in quasars. So far only a handful of HeI* detections have been reported. Using a newly developed method, we detected HeI*3889 absorption line in 101 sources of a well-defined sample of 285 MgII BAL quasars selected from the SDSS DR5. This has increased the number of HeI* BAL quasars by more than one order of magnitude. We further detected HeI*3189 in 50% (52/101) quasars in the sample. The detection fraction of HeI* BALs in MgII BAL quasars is about 35% as a whole, and increases dramatically with increasing spectral signal-to-noise ratios, from 18% at S/N <= 10 to 93% at S/N >= 35. This suggests that HeI* BALs could be detected in most MgII LoBAL quasars, provided spectra S/N is high enough. Such a surprisingly high HeI* BAL fraction is actually predicted from photo-ionization calculations based on a simple BAL model. The result indicates that HeI* absorption lines can be used to search for BAL quasars at low-z, which cannot be identified by ground-based optical spectroscopic survey with commonly seen UV absorption lines. Using HeI*3889, we discovered 19 BAL quasars at z<0.3 from available SDSS spectral database. The fraction of HeI* BAL quasars is similar to that of LoBAL objects.

preprint2015arXiv

Detection of the Intermediate-width Emission Line Region in Quasar OI 287 with the Broad Emission Line Region Obscured by the Dusty Torus

The existence of intermediate-width emission line regions (IELRs) in active galactic nuclei has been discussed for over two decades. A consensus, however, is yet to be arrived at due to the lack of convincing evidence for their detection. We present a detailed analysis of the broadband spectrophotometry of the partially obscured quasar OI 287. The ultraviolet intermediate-width emission lines (IELs) are very prominent, in high contrast to the corresponding broad emission lines (BELs) which are heavily suppressed by dust reddening. Assuming that the IELR is virialized, we estimated its distance to the central black hole of $\sim 2.9$ pc, similar to the dust sublimation radius of $\sim 1.3$ pc. Photo-ionization calculations suggest that the IELR has a hydrogen density of $\sim 10^{8.8}-10^{9.4} ~ \rm cm^{-3}$, within the range of values quoted for the dusty torus near the sublimation radius. Both its inferred location and physical conditions suggest that the IELR originates from the inner surface of the dusty torus. In the spectrum of this quasar, we identified only one narrow absorption-line system associated with the dusty material. With the aid of photo-ionization model calculations, we found that the obscuring material might originate from an outer region of the dusty torus. We speculate that the dusty torus, which is exposed to the central ionizing source, may produce IELs through photo-ionization processes, while also obscure BELs as a natural "coronagraph". Such a "coronagraph" could be found in a large number of partially obscured quasars and be a useful tool to study IELRs.

preprint2015arXiv

Evidence for Photoionization Driven Broad Absorption Line Variability

We present a qualitative analysis of the variability of quasar broad absorption lines using the large multi-epoch spectroscopic dataset of the Sloan Digital Sky Survey Data Release 10. We confirm that variations of absorption lines are highly coordinated among different components of the same ion or the same absorption component of different ions for C IV, Si IV and N V. Furthermore, we show that the equivalent widths of the lines decrease or increase statistically when the continuum brightens or dims. This is further supported by the synchronized variations of emission and absorption line equivalent width, when the well established intrinsic Baldwin effect for emission lines is taken into account. We find that the emergence of an absorption component is usually accompanying with dimming of the continuum while the disappearance of an absorption line component with brightening of the continuum. This suggests that the emergence or disappearance of a C IV absorption component is only the extreme case, when the ionic column density is very sensitive to continuum variations or the continuum variability amplitude is larger. These results support the idea that absorption line variability is driven mainly by changes in the gas ionization in response to continuum variations, that the line-absorbing gas is highly ionized, and in some extreme cases, too highly ionized to be detected in UV absorption lines. Due to uncertainties in the spectroscopic flux calibration, we cannot quantify the fraction of quasars with asynchronized continuum and absorption line variations.

preprint2015arXiv

Flow Patterns around Dark Matter Halos: the Link between Halo Dynamical Properties and Large Scale Tidal Field

We study how halo intrinsic dynamical properties are linked to their formation processes for halos in two mass ranges, $10^{12}-10^{12.5}h^{-1}{\rm M_\odot}$ and $\ge 10^{13}h^{-1}{\rm M_\odot}$, and how both are correlated with the large scale tidal field within which the halos reside at present. Halo merger trees obtained from cosmological $N$-body simulations are used to identify infall halos that are about to merge with their hosts. We find that the tangential component of the infall velocity increases significantly with the strength of the local tidal field, but no strong correlation is found for the radial component. These results can be used to explain how the internal velocity anisotropy and spin of halos depend on environment. The position vectors and velocities of infall halos are aligned with the principal axes of the local tidal field, and the alignment depends on the strength of the tidal field. Opposite accretion patterns are found in weak and strong tidal fields, in the sense that in a weak field the accretion flow is dominated by radial motion within the local structure, while a large tangential component is present in a strong field. These findings can be used to understand the strong alignments we find between the principal axes of the internal velocity ellipsoids of halos and the local tidal field, and their dependence on the strength of tidal field. They also explain why halo spin increases with the strength of local tidal field, but only in weak tidal fields does the spin-tidal field alignment follow the prediction of the tidal torque theory. We discuss how our results may be used to understand the spins of disk galaxies and velocity structures of elliptical galaxies and their correlations with large-scale structure.

preprint2015arXiv

Strong Variability of Overlapping Iron Broad Absorption Lines in Five Radio-selected Quasars

We present the variability study of broad absorption lines (BALs) in a uniformly radio-selected sample of 28 BAL quasars using the archival data from the FIRST Bright Quasar Survey (FBQS) and the Sloan Digital Sky Survey (SDSS), as well as those obtained by ourselves, covering time scales $\sim 1-10$ years in the quasar's rest-frame. To our surprise, 5 quasars showing strong variations are all belong to a special subclass of overlapping iron low ionization BAL (OFeLoBAL) quasars, however, other 4 non-overlapping FeLoBALs (non-OFeLoBALs) are invariable except one case with weak optical depth change. Meanwhile, we also find 6 typical variations of high-ionization and low-ionization BALs in this BAL quasar sample. Photoionization models suggest that OFeLoBALs are formed in a relative dense ($n_e>10^6$ cm$^{-3}$) outflows at a distance from the subparsec to the dozens of parsecs from the continuum source. They differ from those of non-OFeLoBALs, which are likely produced by low-density gas, locating at a distance of hundreds to thousands parsecs. Thus, OFeLoBALs and non-OFeLoBALs, i.e., FeLoBALs with/without strong BAL variations, perhaps represent the bimodality of Fe II absorption, the former is located in the active galactic nucleus environment rather than the host galaxy. We suggest that high density and small distance are the necessary conditions that cause OFeLoBALs. As suggested in the literature, strong BAL variability is possibly due to variability of the covering factor of BAL regions caused by clouds transiting across the line of sight rather than ionization variations.

preprint2014arXiv

ELUCID - Exploring the Local Universe with reConstructed Initial Density field I: Hamiltonian Markov Chain Monte Carlo Method with Particle Mesh Dynamics

Simulating the evolution of the local universe is important for studying galaxies and the intergalactic medium in a way free of cosmic variance. Here we present a method to reconstruct the initial linear density field from an input non-linear density field, employing the Hamiltonian Markov Chain Monte Carlo (HMC) algorithm combined with Particle Mesh (PM) dynamics. The HMC+PM method is applied to cosmological simulations, and the reconstructed linear density fields are then evolved to the present day with N-body simulations. The constrained simulations so obtained accurately reproduce both the amplitudes and phases of the input simulations at various $z$. Using a PM model with a grid cell size of 0.75 Mpc/h and 40 time-steps in the HMC can recover more than half of the phase information down to a scale k~0.85 h/Mpc at high z and to k~3.4 h/Mpc at z=0, which represents a significant improvement over similar reconstruction models in the literature, and indicates that our model can reconstruct the formation histories of cosmic structures over a large dynamical range. Adopting PM models with higher spatial and temporal resolutions yields even better reconstructions, suggesting that our method is limited more by the availability of computer resource than by principle. Dynamic models of structure evolution adopted in many earlier investigations can induce non-Gaussianity in the reconstructed linear density field, which in turn can cause large systematic deviations in the predicted halo mass function. Such deviations are greatly reduced or absent in our reconstruction.

preprint2014arXiv

Outflow and hot dust emission in broad absorption line quasars

We have investigated a sample of 2099 broad absorption line (BAL) quasars with z=1.7-2.2 built from the Sloan Digital Sky Survey Data Release Seven and the Wide-field Infrared Survey. This sample is collected from two BAL quasar samples in the literature, and refined by our new algorithm. Correlations of outflow velocity and strength with hot dust indicator (beta_NIR) and other quasar physical parameters, such as Eddington ratio, luminosity and UV continuum slope, are explored in order to figure out which parameters drive outflows. Here beta_NIR is the near-infrared continuum slope, a good indicator of the amount of hot dust emission relative to accretion disk emission. We confirm previous findings that outflow properties moderately or weakly depends on Eddington ratio, UV slope and luminosity. For the first time, we report moderate and significant correlations of outflow strength and velocity with beta_NIR in BAL quasars. It is consistent with the behavior of blueshifted broad emission lines in non-BAL quasars. The statistical analysis and composite spectra study both reveal that outflow strength and velocity are more strongly correlated with beta_NIR than Eddington ratio, luminosity and UV slope. In particular, the composites show that the entire C IV absorption profile shifts blueward and broadens as beta_NIR increases, while Eddington ratio and UV slope only affect the high and low velocity part of outflows, respectively. We discuss several potential processes and suggest that dusty outflow scenario, i.e. dust is intrinsic to outflows and may contribute to the outflow acceleration, is most likely. The BAL quasar catalog is available from the authors upon request.

preprint2014arXiv

Spin alignments of spiral galaxies within the large-scale structure from SDSS DR7

Using a sample of spiral galaxies selected from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) and Galaxy Zoo 2 (GZ2), we investigate the alignment of spin axes of spiral galaxies with their surrounding large scale structure, which is characterized by the large-scale tidal field reconstructed from the data using galaxy groups above a certain mass threshold. We find that the spin axes of only have weak tendency to be aligned with (or perpendicular to) the intermediate (or minor) axis of the local tidal tensor. The signal is the strongest in a \cluster environment where all the three eigenvalues of the local tidal tensor are positive. Compared to the alignments between halo spins and local tidal field obtained in N-body simulations, the above observational results are in best agreement with those for the spins of inner regions of halos, suggesting that the disk material traces the angular momentum of dark matter halos in the inner regions.

preprint2014arXiv

Unshifted Metastable He I* Mini-Broad Absorption Line System in the Narrow Line Type 1 Quasar SDSS J080248.18$+$551328.9

We report the identification of an unusual absorption line system in the quasar SDSS J080248.18$+$551328.9 and present a detailed study of the system, incorporating follow-up optical and NIR spectroscopy. A few tens of absorption lines are detected, including He I*, Fe II* and Ni II* that arise from metastable or excited levels, as well as resonant lines in Mg I, Mg II, Fe II, Mn II, and Ca II. All of the isolated absorption lines show the same profile of width $Δv\sim 1,500$km s$^{-1}$ centered at a common redshift as that of the quasar emission lines, such as [O II], [S II], and hydrogen Paschen and Balmer series. With narrow Balmer lines, strong optical Fe II multiplets, and weak [O III] doublets, its emission line spectrum is typical for that of a narrow-line Seyfert 1 galaxy (NLS1). We have derived reliable measurements of the gas-phase column densities of the absorbing ions/levels. Photoionization modeling indicates that the absorber has a density of $n_{\rm H} \sim (1.0-2.5)\times 10^5~ {\rm cm}^{-3}$ and a column density of $N_{\rm H} \sim (1.0-3.2)\times 10^{21} \sim {\rm cm}^{-2}$, and is located at $R\sim100-250$ pc from the central super-massive black hole. The location of the absorber, the symmetric profile of the absorption lines, and the coincidence of the absorption and emission line centroid jointly suggest that the absorption gas is originated from the host galaxy and is plausibly accelerated by stellar processes, such as stellar winds \zhy{and/or} supernova explosions. The implications for the detection of such a peculiar absorption line system in an NLS1 are discussed in the context of co-evolution between super-massive black hole growth and host galaxy build-up.

preprint2013arXiv

Alignments of galaxies within cosmic filaments from SDSS DR7

Using a sample of galaxy groups selected from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7), we examine the alignment between the orientation of galaxies and their surrounding large scale structure in the context of the cosmic web. The latter is quantified using the large-scale tidal field, reconstructed from the data using galaxy groups above a certain mass threshold. We find that the major axes of galaxies in filaments tend to be preferentially aligned with the directions of the filaments, while galaxies in sheets have their major axes preferentially aligned parallel to the plane of the sheets. The strength of this alignment signal is strongest for red, central galaxies, and in good agreement with that of dark matter halos in N-body simulations. This suggests that red, central galaxies are well aligned with their host halos, in quantitative agreement with previous studies based on the spatial distribution of satellite galaxies. There is a luminosity and mass dependence that brighter and more massive galaxies in filaments and sheets have stronger alignment signals. We also find that the orientation of galaxies is aligned with the eigenvector associated with the smallest eigenvalue of the tidal tensor. These observational results indicate that galaxy formation is affected by large-scale environments, and strongly suggests that galaxies are aligned with each other over scales comparable to those of sheets and filaments in the cosmic web.

preprint2013arXiv

Anomalously Steep Reddening Law in Quasars: An Exceptional Example Observed in IRAS14026+4341

A fraction of the heavily reddened quasars require a reddening curve which is even steeper than that of the Small Magellanic Cloud. In this paper, we thoroughly characterize the anomalously steep reddening law in quasars, via an exceptional example observed in IRAS 14026+4341. By comparing the observed spectrum to the quasar composite spectrum, we derive a reddening curve in the rest-frame wavelength range of 1200 Å--10000 Å. It is featured with a steep rise at wavelengths shorter than 3000 Å, but no significant reddening at longer wavelengths. The absence of dust reddening in optical continuum is confirmed by the normal broad-line Balmer decrement (the H$α$/H$β$ ratio) in IRAS 14026+4341. The anomalous reddening curve can be satisfactorily reproduced by a dust model containing silicate grains in a power-law size distribution, $dn(a)/da \propto a^{-1.4}$, truncated at a maximum size $a_{max}=70 {\rm nm}$. The unusual size distribution may be caused by the destruction of large "stardust" grains by quasar activities or a different dust formation mechanism (i.e., the in situ formation of dust grains in quasar outflows). It is also possible that the analogies of the dust grains observed toward the Galactic center is responsible for the steep reddening curve. In addition, we find that IRAS 14026+4341 is a weak emission-line quasar (i.e., PHL 1811 analogies) with heavy dust reddening and blueshifted broad absorption lines.

preprint2013arXiv

Outflow and hot dust emission in high redshift quasars

Correlations of hot dust emission with outflow properties are investigated, based on a large z~2 non-broad absorption lines quasar sample built from the Wide-field Infrared Survey and the Sloan Digital Sky Survey data releases. We use the near infrared slope and the infrared to UV luminosity ratio to indicate the hot dust emission relative to the emission from the accretion disk. In our luminous quasars, these hot dust emission indicators are almost independent of the fundamental parameters, such as luminosity, Eddington ratio and black hole mass, but moderately dependent on the blueshift and asymmetry index (BAI) and full width at half-maximum (FWHM) of CIV lines. Interestingly, the latter two correlations dramatically strengthen with increasing Eddington ratio. We suggest that, in high Eddington ratio quasars, CIV regions are dominated by outflows so the BAI and FWHM(CIV) can reliably reflect the general properties and velocity of outflows, respectively. While in low Eddington ratio quasars, CIV lines are primarily emitted by virialized gas so the BAI and FWHM(CIV) become less sensitive to outflows. Therefore, the correlations for the highest Eddington ratio quasars are more likely to represent the true dependence of hot dust emission on outflows and the correlations for the entire sample are significantly diluted by the low Eddington ratio quasars. Our results show that an outflow with a large BAI or velocity can double the hot dust emission on average. We suggest that outflows either contain hot dust in themselves or interact with the dusty interstellar medium or torus.

preprint2013arXiv

Reconstructing the Initial Density Field of the Local Universe: Method and Test with Mock Catalogs

Our research objective in this paper is to reconstruct an initial linear density field, which follows the multivariate Gaussian distribution with variances given by the linear power spectrum of the current CDM model and evolves through gravitational instability to the present-day density field in the local Universe. For this purpose, we develop a Hamiltonian Markov Chain Monte Carlo method to obtain the linear density field from a posterior probability function that consists of two components: a prior of a Gaussian density field with a given linear spectrum, and a likelihood term that is given by the current density field. The present-day density field can be reconstructed from galaxy groups using the method developed in Wang et al. (2009a). Using a realistic mock SDSS DR7, obtained by populating dark matter haloes in the Millennium simulation with galaxies, we show that our method can effectively and accurately recover both the amplitudes and phases of the initial, linear density field. To examine the accuracy of our method, we use $N$-body simulations to evolve these reconstructed initial conditions to the present day. The resimulated density field thus obtained accurately matches the original density field of the Millennium simulation in the density range 0.3 <= rho/rho_mean <= 20 without any significant bias. Especially, the Fourier phases of the resimulated density fields are tightly correlated with those of the original simulation down to a scale corresponding to a wavenumber of ~ 1 h/Mpc, much smaller than the translinear scale, which corresponds to a wavenumber of ~ 0.15 h\Mpc.

preprint2012arXiv

Discovery of Balmer Broad Absorption Lines in the Quasar LBQS 1206+1052

We report the discovery of Balmer broad absorption lines (BALs) in the quasar LBQS 1206+1052 and present a detailed analysis of the peculiar absorption line spectrum. Besides Mg II $λλ2796, 2803$ doublet, BALs are also detected in He I* multiplet at $λλ2946, 3189, 3889$ Åarising from metastable helium $2^3S$ level, and in H$α$ and H$β$ from excited hydrogen H I* $n=2$ level, which are rarely seen in quasar spectra. We identify two components in the BAL troughs of $Δv\sim$2000 km s$^{-1}$ width: One component shows an identical profile in H I*, \hei* and \mgii with its centroid blueshifted by $-v_{\rm c}\approx 726$ km\ s$^{-1}$. The other component is detected in \hei* and \mgii with $-v_{\rm c}\approx 1412$ km s$^{-1}$. We estimate the column densities of H I*, He I*, and Mg II, and compare them with possible level population mechanisms. Our results favor the scenario that the Balmer BALs originate in a partially ionized region with a column density of $N_{\rm H}\sim 10^{21-22}$ cm$^{-2}$ for an electron density of $n_e\sim 10^{6-8} $cm$^{-3}$ via Ly$α$ resonant scattering pumping. The harsh conditions needed may help to explain the rarity of Balmer absorption line systems in quasar spectra. With an $i$-band PSF magnitude of 16.50, LBQS 1206+1052 is the brightest Balmer-BAL quasar ever reported. Its high brightness and unique spectral properties make LBQS 1206+1052 a promising candidate for follow-up high-resolution spectroscopy, multi-band observations, and long-term monitoring.

preprint2012arXiv

Extreme Coronal Line Emitters: Tidal Disruption of Stars by Massive Black Holes in Galactic Nuclei?

Tidal disruption of stars by supermassive black holes at the centers of galaxies is expected to produce unique emission line signatures, which have not yet been explored adequately. Here we report the discovery of extremely strong coronal lines from [Fe X] up to [Fe XIV] in a sample of seven galaxies (including two recently reported cases), that we interpret as such signatures. This is the first systematic search for objects of this kind, by making use of the immense database of the Sloan Digital Sky Survey. The galaxies, which are non-active as evidenced by the narrow-line ratios, show broad emission lines of complex profiles in more than half of the sample. Both the high ionization coronal lines and the broad lines turn out to be fading on time scales of years in objects observed with spectroscopic follow-ups, suggesting their transient nature. Variations of inferred non-stellar continua, which have absolute magnitudes of at least -16 to -18 mag in the g band, are also detected in more than half of the sample. These extreme coronal line emitters reside in sub-L_* disk galaxies (-21.3 < M_i < -18.5) with small stellar velocity dispersions. The sample seems to form two distinct types based on the presence or absence of the [Fe VII] lines, with the latter having relatively low luminosities of [O III], [Fe XI], and the host galaxies. These characteristics can most naturally be understood in the context of transient accretion onto intermediate mass black holes at galactic centers following tidal disruption of stars in a gas-rich environment. We estimate the incidence of such events to be around 10^-5 per year for a galaxy with -21.5 < M_i < -18.5.

preprint2012arXiv

Metallicity and Quasar Outflows

Correlations are investigated of the outflow strength of quasars, as measured by the blueshift and asymmetry index (BAI) of the CIV line (Wang et al. 2011), with intensities and ratios of broad emission lines, based on composite quasar spectra built from the Sloan Digital Sky Survey. We find that most of the line ratios of other ions to CIV prominently increases with BAI. These behaviors can be well understood in the context of increasing metallicity with BAI. The strength of dominant coolant, CIV line, decreases and weak collisionally excited lines increase with gas metallicity as a result of the competition between different line coolants. Using SiIV+OIV]/CIV as an indicator of gas metallicity, we present, for the first time, a strong correlation between the metallicitiy and the outflow strength of quasars over a wide range of 1.7 to 6.9 times solar abundance. Our result implies that the metallicity plays an important role in the formation of quasar outflows, likely via affecting outflow acceleration. This effect may have a profound impact on galaxy evolution via momentum feedback and chemical enrichment.

preprint2012arXiv

Resolving Subhaloes' Lives with the Hierarchical Bound-Tracing Algorithm

We develop a new code, the Hierarchical Bound-Tracing (HBT for short) code, to find and trace dark matter subhaloes in simulations based on the merger hierarchy of dark matter haloes. Application of this code to a recent benchmark test of finding subhaloes demonstrates that HBT stands as one of the best codes to trace the evolutionary history of subhaloes. The success of the code lies in its careful treatment of the complex physical processes associated with the evolution of subhaloes and in its robust unbinding algorithm with an adaptive source subhalo management. We keep a full record of the merger hierarchy of haloes and subhaloes, and allow growth of satellite subhaloes through accretion from its "satellite-of-satellites", hence allowing mergers among satellites. Local accretion of background mass is omitted, while rebinding of stripped mass is allowed. The justification of these treatments is provided by case studies of the lives of individual subhaloes and by the success in finding the complete subhalo catalogue. We compare our result to other popular subhalo finders and show that HBT is able to well resolve subhaloes in high density environment and keep strict physical track of subhaloes' merger history. This code is fully parallelized and freely available upon request to the authors.

preprint2011arXiv

Coexistence of Gravitationally Bound and Radiation Driven CIV Emission Line Regions in Active Galactic Nuclei

There are mutually contradictory views in the literature of the kinematics and structure of high-ionization line (e.g. CIV) emitting regions in active galactic nuclei (AGNs). Two kinds of broad emission line region (BELR) models have been proposed, outflow and gravitationally bound BELR, which are supported respectively by blueshift of the CIV line and reverberation mapping observations. To reconcile these two apparently different models, we present a detailed comparison study between the CIV and MgII lines using a sample of AGNs selected from the Sloan Digital Sky Survey. We find that the kinematics of the CIV region is different from that of MgII, which is thought to be controlled by gravity. A strong correlation is found between the blueshift and asymmetry of the CIV profile and the Eddington ratio. This provides strong observational support for the postulation that the outflow is driven by radiation pressure. In particular, we find robust evidence that the CIV line region is largely dominated by outflow at high Eddington ratios, while it is primarily gravitationally bounded at low Eddington ratios. Our results indicate that these two emitting regions coexist in most of AGNs. The emission strength from these two gases varies smoothly with Eddington ratio in opposite ways. This explanation naturally reconciles the apparently contradictory views proposed in previous studies. Finally, candidate models are discussed which can account for both, the enhancement of outflow emission and suppression of normal BEL, in AGN with high Eddington ratios.

preprint2011arXiv

Dust reddening in star-forming galaxies

We present empirical relations between the global dust reddening and other physical galaxy properties including the Halpha luminosity, Halpha surface brightness, metallicity and axial ratio for star-forming disc galaxies. The study is based on a large sample of ~22 000 well-defined star-forming galaxies selected from the Sloan Digital Sky Survey (SDSS). The reddening parameterized by color excess E(B-V) is derived from the Balmer decrement. Besides the dependency of reddening on Halpha luminosity / surface brightness and gas phase metallicity, it is also correlated with the galaxy inclination, in the sense that edge-on galaxies are more attenuated than face-on galaxies at a give intrinsic luminosity. In light of these correlations, we present the empirical formulae of E(B-V) as a function of these galaxy properties, with a scatter of only 0.07 mag. The empirical relation can be reproduced if most dust attenuation to the HII region is due to diffuse interstellar dust distributing in a disc thicker than that of HII regions. The empirical formulae can be incorporated into semi-analytical models of galaxy formation and evolution to estimate the dust reddening and enable comparison with observations more practically.

preprint2011arXiv

Reconstructing the Cosmic Velocity and Tidal Fields with Galaxy Groups Selected from the Sloan Digital Sky Survey

[abridge]Cosmic velocity and tidal fields are important for the understanding of the cosmic web and the environments of galaxies, and can also be used to constrain cosmology. In this paper, we reconstruct these two fields in SDSS volume from dark matter halos represented by galaxy groups. Detailed mock catalogues are used to test the reliability of our method against uncertainties arising from redshift distortions, survey boundaries, and false identifications of groups by our group finder. We find that both the velocity and tidal fields, smoothed on a scale of ~2Mpc/h, can be reliably reconstructed in the inner region (~66%) of the survey volume. The reconstructed tidal field is used to split the cosmic web into clusters, filaments, sheets, and voids, depending on the sign of the eigenvalues of tidal tensor. The reconstructed velocity field nicely shows how the flows are diverging from the centers of voids, and converging onto clusters, while sheets and filaments have flows that are convergent along one and two directions, respectively. We use the reconstructed velocity field and the Zel'dovich approximation to predict the mass density field in the SDSS volume as function of redshift, and find that the mass distribution closely follows the galaxy distribution even on small scales. We find a large-scale bulk flow of about 117km/s in a very large volume, equivalent to a sphere with a radius of ~170Mpc/h, which seems to be produced by the massive structures associated with the SDSS Great Wall. Finally, we discuss potential applications of our reconstruction to study the environmental effects of galaxy formation, to generate initial conditions for simulations of the local Universe, and to constrain cosmological models. The velocity, tidal and density fields in the SDSS volume, specified on a Cartesian grid with a spatial resolution of ~700kpc/h, are available from the authors upon request.

preprint2011arXiv

What controls the Fe II strength in active galactic nuclei?

We used a large, homogeneous sample of 4178 z <= 0.8 Seyfert 1 galaxies and QSOs selected from the Sloan Digital Sky Survey to investigate the strength of Fe II emission and its correlation with other emission lines and physical parameters of active galactic nuclei. We find that the strongest correlations of almost all the emission-line intensity ratios and equivalent widths (EWs) are with the Eddington ratio (L/L_{Edd}), rather than with the continuum luminosity at 5100Å (L_{5100}) or black hole mass (M_{BH}); the only exception is the EW of ultraviolet Fe II emission, which does not correlate at all with broad-line width, L_{5100}, M_{BH}, or L/L_{Edd}. By contrast, the intensity ratios of both the ultraviolet and optical Fe II emission to Mg II λ2800 correlate quite strongly with L/L_{Edd}. Interestingly, among all the emission lines in the near-UV and optical, the EW of narrow optical Fe II emission has the strongest correlation with L/L_{Edd}. We suggest that the variation of the emission-line strength in active galaxies is regulated by L/L_{Edd} because it governs the global distribution of the column density of the clouds gravitationally bound in the line-emitting region, as well as its overall gas supply. The systematic dependence on L/L_{Edd} must be corrected when using the FeII/MgII intensity ratio as a measure of the Fe/Mg abundance ratio to study the history of chemical evolution in QSO environments.

preprint2010arXiv

Low-z MgII Broad Absorption-Line Quasars from the Sloan Digital Sky Survey

We present a sample of 68 low-z MgII low-ionization broad absorption-line (loBAL) quasars. The sample is uniformly selected from the SDSS5 according to the following criteria: (1) 0.4<z<0.8, (2) median S/N>7, and (3) MgII absorption-line width > 1600 \kms. The last criterion is a trade-off between the completeness and consistency with respect to the canonical definition of BAL quasars that have the `balnicity index' BI>0 in CIV BAL. We adopted such a criterion to ensure that ~90% of our sample are classical BAL quasars and the completeness is ~80%, based on extensive tests using high-z quasar samples with measurements of both CIV and MgII BALs. We found (1) MgII BAL is more frequently detected in quasars with narrower Hbeta emission-line, weaker [OIII] emission-line, stronger optical FeII multiplets and higher luminosity. In term of fundamental physical parameters of a black hole accretion system, loBAL fraction is significantly higher in quasars with a higher Eddington ratio than those with a lower Eddington ratio. The fraction is not dependent on the black hole mass in the range concerned. The overall fraction distribution is broad, suggesting a large range of covering factor of the absorption material. (2) [OIII]-weak loBAL quasars averagely show undetected [NeV] emission line and a very small line ratio of [NeV] to [OIII]. However, the line ratio in non-BAL quasars, which is much larger than that in [OIII]-weak loBAL quasars, is independent of the strength of the [OIII] line. (3) loBAL and non-loBAL quasars have similar colors in near-infrared to optical band but different colors in ultraviolet. (4) Quasars with MgII absorption lines of intermediate width are indistinguishable from the non-loBAL quasars in optical emission line properties but their colors are similar to loBAL quasars, redder than non-BAL quasars. We also discuss the implication of these results.

preprint2010arXiv

The Prevalence of Narrow Optical Fe II Emission Lines in Type 1 Active Galactic Nuclei

From detailed spectral analysis of a large sample of low-redshift active galactic nuclei (AGNs) selected from the Sloan Digital Sky Survey, we demonstrate---statistically for the first time---that narrow optical Fe II emission lines, both permitted and forbidden, are prevalent in type 1 AGNs. Remarkably, these optical lines are completely absent in type 2 AGNs, across a wide luminosity range, from Seyfert 2 galaxies to type 2 quasars. We suggest that the narrow FeII-emitting gas is confined to a disk-like geometry in the innermost regions of the narrow-line region on physical scales smaller than the obscuring torus.

preprint2009arXiv

Estimating Black Hole Masses in Active Galactic Nuclei Using the MgII 2800 Emission Line

We investigate the relationship between the linewidths of broad Mg II \lambda2800 and Hbeta in active galactic nuclei (AGNs) to refine them as tools to estimate black hole (BH) masses. We perform a detailed spectral analysis of a large sample of AGNs at intermediate redshifts selected from the Sloan Digital Sky Survey, along with a smaller sample of archival ultraviolet spectra for nearby sources monitored with reverberation mapping. Careful attention is devoted to accurate spectral decomposition, especially in the treatment of narrow-line blending and Fe II contamination. We show that, contrary to popular belief, the velocity width of Mg II tends to be smaller than that of Hbeta, suggesting that the two species are not cospatial in the broad-line region. Using these findings and recently updated BH mass measurements from reverberation mapping, we present a new calibration of the empirical prescriptions for estimating virial BH masses for AGNs using the broad Mg II and Hbeta lines. We show that the BH masses derived from our new formalisms show subtle but important differences compared to some of the mass estimators currently used in the literature.

preprint2009arXiv

The distribution of ejected subhalos and its implication for halo assembly bias

Using a high-resolution cosmological $N$-body simulation, we identify the ejected population of subhalos, which are halos at redshift $z=0$ but were once contained in more massive `host' halos at high redshifts. The fraction of the ejected subhalos in the total halo population of the same mass ranges from 9% to 4% for halo masses from $\sim 10^{11}$ to $\sim 10^{12}\msun$. Most of the ejected subhalos are distributed within 4 times the virial radius of their hosts. These ejected subhalos have distinct velocity distribution around their hosts in comparison to normal halos. The number of subhalos ejected from a host of given mass increases with the assembly redshift of the host. Ejected subhalos in general reside in high-density regions, and have a much higher bias parameter than normal halos of the same mass. They also have earlier assembly times, so that they contribute to the assembly bias of dark matter halos seen in cosmological simulations. However, the assembly bias is {\it not} dominated by the ejected population, indicating that large-scale environmental effects on normal halos are the main source for the assembly bias.