Source author record

X. Kang

X. Kang 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

25works
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

25 published item(s)

preprint2025arXiv

All-sky neutrino point-source search with IceCube combined track and cascade data

Despite extensive efforts, discovery of high-energy astrophysical neutrino sources remains elusive. We present an event-level simultaneous maximum likelihood analysis of tracks and cascades using IceCube data collected from 04/06/2008 to 05/23/2022 to search the whole sky for neutrino sources and, using a source catalog, for coincidence of neutrino emission with gamma-ray emission. This is the first time a simultaneous fit of different detection channels is used to conduct a time-integrated all-sky scan with IceCube. Combining all-sky tracks, with superior pointing-power and sensitivity in the northern sky, with all-sky cascades, with good energy-resolution and sensitivity in the southern sky, we have developed the most sensitive point-source search to date by IceCube which targets the entire sky. The most significant point in the northern sky aligns with NGC 1068, a Seyfert II galaxy, which, from the catalog search, shows a 3.5$σ$ excess over background after accounting for trials. The most significant point in the southern sky does not align with any source in the catalog and is not significant after accounting for trials. A search for the single most significant Gaussian flare at the locations of NGC 1068, PKS 1424+240, and the southern highest significance point shows results consistent with expectations for steady emission. Notably, this is the first time that a flare shorter than four years has been excluded as being responsible for NGC 1068's emergence as a neutrino source. Our results show that combining tracks and cascades when conducting neutrino source searches improves sensitivity and can lead to new discoveries.

preprint2022arXiv

Density of GeV muons in air showers measured with IceTop

We present a measurement of the density of GeV muons in near-vertical air showers using three years of data recorded by the IceTop array at the South Pole. Depending on the shower size, the muon densities have been measured at lateral distances between 200 m and 1000 m. From these lateral distributions, we derive the muon densities as functions of energy at reference distances of 600 m and 800 m for primary energies between 2.5 PeV and 40 PeV and between 9 PeV and 120 PeV, respectively. The muon densities are determined using, as a baseline, the hadronic interaction model Sibyll 2.1 together with various composition models. The measurements are consistent with the predicted muon densities within these baseline interaction and composition models. The measured muon densities have also been compared to simulations using the post-LHC models EPOS-LHC and QGSJet-II.04. The result of this comparison is that the post-LHC models together with any given composition model yield higher muon densities than observed. This is in contrast to the observations above 1 EeV where all model simulations yield for any mass composition lower muon densities than the measured ones. The post-LHC models in general feature higher muon densities so that the agreement with experimental data at the highest energies is improved but the muon densities are not correct in the energy range between 2.5 PeV and about 100 PeV.

preprint2022arXiv

First all-flavor search for transient neutrino emission using 3-years of IceCube DeepCore data

Since the discovery of a flux of high-energy astrophysical neutrinos, searches for their origins have focused primarily at TeV-PeV energies. Compared to sub-TeV searches, high-energy searches benefit from an increase in the neutrino cross section, improved angular resolution on the neutrino direction, and a reduced background from atmospheric neutrinos and muons. However, the focus on high energy does not preclude the existence of sub-TeV neutrino emission where IceCube retains sensitivity. Here we present the first all-flavor search from IceCube for transient emission of low-energy neutrinos, focusing on the energy region of 5.6-100 GeV using three years of data obtained with the IceCube-DeepCore detector. We find no evidence of transient neutrino emission in the data, thus leading to a constraint on the volumetric rate of astrophysical transient sources in the range of $\sim 705-2301\, \text{Gpc}^{-3}\, \text{yr}^{-1}$ for sources following a subphotospheric energy spectrum with a mean energy of 100 GeV and a bolometric energy of $10^{52}$ erg.

preprint2022arXiv

Precision tests of Quantum Mechanics and CPT symmetry with entangled neutral kaons at KLOE

The quantum interference between the decays of entangled neutral kaons is studied in the process $ϕ\rightarrow K_S K_L \rightarrowπ^+π^-π^+π^-$, which exhibits the characteristic Einstein--Podolsky--Rosen correlations that prevent both kaons to decay into $π^+π^-$ at the same time. This constitutes a very powerful tool for testing at the utmost precision the quantum coherence of the entangled kaon pair state, and to search for tiny decoherence and CPT violation effects, which may be justified in a quantum gravity framework. The analysed data sample was collected with the KLOE detector at DA$Φ$NE, the Frascati $ϕ$-factory, and corresponds to an integrated luminosity of about 1.7 fb$^{-1}$, i.e. to about $1.7 \times 10^9$ $ϕ\rightarrow K_S K_L$ decays produced. From the fit of the observed $Δt$ distribution, being $Δt$ the difference of the kaon decay times, the decoherence and CPT violation parameters of various phenomenological models are measured with a largely improved accuracy with respect to previous analyses. The results are consistent with no deviation from quantum mechanics and CPT symmetry, while for some parameters the precision reaches the interesting level at which -- in the most optimistic scenarios -- quantum gravity effects might show up. They provide the most stringent limits up to date on the considered models.

preprint2022arXiv

Search for Relativistic Magnetic Monopoles with Eight Years of IceCube Data

We present an all-sky 90\% confidence level upper limit on the cosmic flux of relativistic magnetic monopoles using 2886 days of IceCube data. The analysis was optimized for monopole speeds between 0.750$c$ and 0.995$c$, without any explicit restriction on the monopole mass. We constrain the flux of relativistic cosmic magnetic monopoles to a level below $2.0\times 10^{-19} {\textrm{cm}}^{-2} {\textrm{s}}^{-1} {\textrm{sr}}^{-1}$ over the majority of the targeted speed range. This result constitutes the most strict upper limit to date for magnetic monopoles above the Cherenkov threshold and up to $β\sim 0.995$ and fills the gap between existing limits on the cosmic flux of non-relativistic and ultrarelativistic magnetic monopoles

preprint2022arXiv

Search for Spatial Correlations of Neutrinos with Ultra-High-Energy Cosmic Rays

For several decades, the origin of ultra-high-energy cosmic rays (UHECRs) has been an unsolved question of high-energy astrophysics. One approach for solving this puzzle is to correlate UHECRs with high-energy neutrinos, since neutrinos are a direct probe of hadronic interactions of cosmic rays and are not deflected by magnetic fields. In this paper, we present three different approaches for correlating the arrival directions of neutrinos with the arrival directions of UHECRs. The neutrino data is provided by the IceCube Neutrino Observatory and ANTARES, while the UHECR data with energies above $\sim$50 EeV is provided by the Pierre Auger Observatory and the Telescope Array. All experiments provide increased statistics and improved reconstructions with respect to our previous results reported in 2015. The first analysis uses a high-statistics neutrino sample optimized for point-source searches to search for excesses of neutrinos clustering in the vicinity of UHECR directions. The second analysis searches for an excess of UHECRs in the direction of the highest-energy neutrinos. The third analysis searches for an excess of pairs of UHECRs and highest-energy neutrinos on different angular scales. None of the analyses has found a significant excess, and previously reported over-fluctuations are reduced in significance. Based on these results, we further constrain the neutrino flux spatially correlated with UHECRs.

preprint2021arXiv

A search for neutrino emission from cores of Active Galactic Nuclei

The sources of the majority of the high-energy astrophysical neutrinos observed with the IceCube neutrino telescope at the South Pole are unknown. So far, only a gamma-ray blazar was compellingly associated with the emission of high-energy neutrinos. In addition, several studies suggest that the neutrino emission from the gamma-ray blazar population only accounts for a small fraction of the total astrophysical neutrino flux. In this work we probe the production of high-energy neutrinos in the cores of Active Galactic Nuclei (AGN), induced by accelerated cosmic rays in the accretion disk region. We present a likelihood analysis based on eight years of IceCube data, searching for a cumulative neutrino signal from three AGN samples created for this work. The neutrino emission is assumed to be proportional to the accretion disk luminosity estimated from the soft X-ray flux. Next to the observed soft X-ray flux, the objects for the three samples have been selected based on their radio emission and infrared color properties. For the largest sample in this search, an excess of high-energy neutrino events with respect to an isotropic background of atmospheric and astrophysical neutrinos is found, corresponding to a post-trial significance of 2.60$σ$. Assuming a power-law spectrum, the best-fit spectral index is $2.03^{+0.14}_{-0.11}$, consistent with expectations from particle acceleration in astrophysical sources. If interpreted as a genuine signal with the assumptions of a proportionality of X-ray and neutrino fluxes and a model for the sub-threshold flux distribution, this observation implies that at 100 TeV, 27$\%$ - 100$\%$ of the observed neutrinos arise from particle acceleration in the core of AGN.

preprint2021arXiv

Detection of a particle shower at the Glashow resonance with IceCube

The Glashow resonance describes the resonant formation of a $W^-$ boson during the interaction of a high-energy electron antineutrino with an electron, peaking at an antineutrino energy of 6.3 petaelectronvolts (PeV) in the rest frame of the electron. Whereas this energy scale is out of reach for currently operating and future planned particle accelerators, natural astrophysical phenomena are expected to produce antineutrinos with energies beyond the PeV scale. Here we report the detection by the IceCube neutrino observatory of a cascade of high-energy particles (a particle shower) consistent with being created at the Glashow resonance. A shower with an energy of $6.05 \pm 0.72$ PeV (determined from Cherenkov radiation in the Antarctic Ice Sheet) was measured. Features consistent with the production of secondary muons in the particle shower indicate the hadronic decay of a resonant $W^-$ boson, confirm that the source is astrophysical and provide improved directional localization. The evidence of the Glashow resonance suggests the presence of electron antineutrinos in the astrophysical flux, while also providing further validation of the standard model of particle physics. Its unique signature indicates a method of distinguishing neutrinos from antineutrinos, thus providing a way to identify astronomical accelerators that produce neutrinos via hadronuclear or photohadronic interactions, with or without strong magnetic fields. As such, knowledge of both the flavour (that is, electron, muon or tau neutrinos) and charge (neutrino or antineutrino) will facilitate the advancement of neutrino astronomy.

preprint2021arXiv

IceCube Data for Neutrino Point-Source Searches Years 2008-2018

IceCube has performed several all-sky searches for point-like neutrino sources using track-like events, including a recent time-integrated analysis using 10 years of IceCube data. This paper accompanies the public data release of these neutrino candidates detected by IceCube between April 6, 2008 and July 8, 2018. The selection includes through-going tracks, primarily due to muon neutrino candidates, that reach the detector from all directions, as well as neutrino track events that start within the instrumented volume. An updated selection and reconstruction for data taken after April 2012 slightly improves the sensitivity of the sample. While more than 80% of the sample overlaps between the old and new versions, differing events can lead to changes relative to the previous 7 year event selection. An a posteriori estimate of the significance of the 2014-2015 TXS flare is reported with an explanation of observed discrepancies with previous results. This public data release, which includes 10 years of data and binned detector response functions for muon neutrino signal events, shows improved sensitivity in generic time-integrated point source analyses and should be preferred over previous releases.

preprint2021arXiv

Improved Characterization of the Astrophysical Muon-Neutrino Flux with 9.5 Years of IceCube Data

We present a measurement of the high-energy astrophysical muon-neutrino flux with the IceCube Neutrino Observatory. The measurement uses a high-purity selection of ~650k neutrino-induced muon tracks from the Northern celestial hemisphere, corresponding to 9.5 years of experimental data. With respect to previous publications, the measurement is improved by the increased size of the event sample and the extended model testing beyond simple power-law hypotheses. An updated treatment of systematic uncertainties and atmospheric background fluxes has been implemented based on recent models. The best-fit single power-law parameterization for the astrophysical energy spectrum results in a normalization of $ϕ_{\mathrm{@100TeV}}^{ν_μ+\barν_μ} = 1.44_{-0.26}^{+0.25} \times 10^{-18}\,\mathrm{GeV}^{-1}\mathrm{cm}^{-2}\mathrm{s}^{-1}\mathrm{sr}^{-1}$ and a spectral index $γ_{\mathrm{SPL}} = 2.37_{-0.09}^{+0.09}$, constrained in the energy range from $15\,\mathrm{TeV}$ to $5\,\mathrm{PeV}$. The model tests include a single power law with a spectral cutoff at high energies, a log-parabola model, several source-class specific flux predictions from the literature and a model-independent spectral unfolding. The data is well consistent with a single power law hypothesis, however, spectra with softening above one PeV are statistically more favorable at a two sigma level.

preprint2021arXiv

Search for High-Energy Neutrinos from Ultra-Luminous Infrared Galaxies with IceCube

Ultra-luminous infrared galaxies (ULIRGs) have infrared luminosities $L_{\mathrm{IR}} \geq 10^{12} L_{\odot}$, making them the most luminous objects in the infrared sky. These dusty objects are generally powered by starbursts with star-formation rates that exceed $100~ M_{\odot}~ \mathrm{yr}^{-1}$, possibly combined with a contribution from an active galactic nucleus. Such environments make ULIRGs plausible sources of astrophysical high-energy neutrinos, which can be observed by the IceCube Neutrino Observatory at the South Pole. We present a stacking search for high-energy neutrinos from a representative sample of 75 ULIRGs with redshift $z \leq 0.13$ using 7.5 years of IceCube data. The results are consistent with a background-only observation, yielding upper limits on the neutrino flux from these 75 ULIRGs. For an unbroken $E^{-2.5}$ power-law spectrum, we report an upper limit on the stacked flux $Φ_{ν_μ+ \barν_μ}^{90\%} = 3.24 \times 10^{-14}~ \mathrm{TeV^{-1}~ cm^{-2}~ s^{-1}}~ (E/10~ \mathrm{TeV})^{-2.5}$ at 90% confidence level. In addition, we constrain the contribution of the ULIRG source population to the observed diffuse astrophysical neutrino flux as well as model predictions.

preprint2021arXiv

Upper limit on the $η\toπ^{+}π^{-}$ branching fraction with the KLOE experiment

Based on an integrated luminosity of 1.61 fb$^{-1}$ $e^+e^-$ collision data collected with the KLOE detector at DA$Φ$NE, the Frascati $ϕ$-factory, a search for the $P$- and $CP$-violating decay $η\toπ^{+}π^{-}$ has been performed. Radiative $ϕ\toηγ$ decay is exploited to access the $η$ mesons. No signal is observed in the $π^{+}π^{-}$ invariant mass spectrum, and the upper limit on the branching fraction at 90\% confidence level is determined to be ${\mathcal B}(η\toπ^{+}π^{-})<4.9\times10^{-6}$, which is approximately three times smaller than the previous KLOE result. From the combination of these two measurements we get ${\mathcal B}(η\toπ^{+}π^{-}) < 4.4\times10^{-6}$ at 90\% confidence level.

preprint2020arXiv

Measurement of the branching fraction for the decay $K_S \to πμν$ with the KLOE detector

Based on a sample of 300 million $K_S$ mesons produced in $ϕ\to K_L K_S$ decays recorded by the KLOE experiment at the DA$Φ$NE $e^+e^-$ collider we have measured the branching fraction for the decay $K_S \to πμν$. The $K_S$ mesons are identified by the interaction of $K_L$ mesons in the detector. The $K_S \to πμν$ decays are selected by a boosted decision tree built with kinematic variables and by a time-of-flight measurement. Signal efficiencies are evaluated with data control samples of $K_L \to πμν$ decays. A fit to the reconstructed muon mass distribution finds $7223 \pm 180$ signal events. Normalising to the $K_S \to π^+ π^-$ decay events the result for the branching fraction is $\mathcal{B}(K_S \to πμν) = (4.56 \pm 0.11_{\rm stat} \pm 0.17_{\rm syst})\times10^{-4}$.

preprint2020arXiv

The Roles of Mass and Environment in the Quenching of Galaxies. II

We take advantage of an analytic model of galaxy formation coupled to the merger tree of an N-body simulation to study the roles of environment and stellar mass in the quenching of galaxies. The model has been originally set in order to provide the observed evolution of the stellar mass function as well as reasonable predictions of the star formation rate-stellar mass relation, from high redshift to the present time. We analyse the stellar mass and environmental quenching efficiencies and their dependence on stellar mass, halo mass (taken as a proxy for the environment) and redshift. Our analysis shows that the two quenching efficiencies are redshift, stellar and halo mass dependent, and that the halo mass is also a good proxy for the environment. The environmental quenching increases with decreasing redshift and is inefficient below $\log M_* \sim 9.5$, reaches the maximum value at $\log M_* \sim 10.5$, and decreases again, becoming poorly efficient at very high stellar mass ($\log M_* \gtrsim 11.5$). Central and satellites galaxies are mass quenched differently: for the former, the quenching efficiency depends very weakly on redshift, but strongly on stellar mass; for the latter, it strongly depends on both stellar mass and redshift in the range $10\lesssim \log M_* \lesssim 11$. According to the most recent observational results, we find that the two quenching efficiencies are not separable: intermediate mass galaxies are environmental quenched faster, as well as intermediate/massive galaxies in more massive haloes. At stellar masses lower than $\log M_* \lesssim 9.5$ both quenching mechanisms become inefficient, independently of the redshift.

preprint2016arXiv

NIHAO VI. The hidden discs of simulated galaxies

Detailed studies of galaxy formation require clear definitions of the structural components of galaxies. Precisely defined components also enable better comparisons between observations and simulations. We use a subsample of eighteen cosmological zoom-in simulations from the NIHAO project to derive a robust method for defining stellar kinematic discs in galaxies. Our method uses Gaussian Mixture Models in a 3D space of dynamical variables. The NIHAO galaxies have the right stellar mass for their halo mass, and their angular momenta and Sérsic indices match observations. While the photometric disc-to-total ratios are close to 1 for all the simulated galaxies, the kinematic ratios are around ~0.5. Thus, exponential structure does not imply a cold kinematic disc. Above log(M*)~9.5, the decomposition leads to thin discs and spheroids that have clearly different properties, in terms of angular momentum, rotational support, ellipticity, [Fe/H] and [O/Fe]. At log(M*)<9.5, the decomposition selects discs and spheroids with less distinct properties. At these low masses, both the discs and spheroids have exponential profiles with high minor-to-major axes ratios, i.e. thickened discs.

preprint2015arXiv

Semi-Analytic Model Predictions of the Galaxy Population in Proto-clusters

We investigate the galaxy population in simulated proto-cluster regions using a semi-analytic model of galaxy formation, coupled to merger trees extracted from N-body simulations. We select the most massive clusters at redshift $z=0$ from our set of simulations, and follow their main progenitors back in time. The analysis shows that proto-cluster regions are dominated by central galaxies and their number decreases with time as many become satellites, clustering around the central object. In agreement with observations, we find an increasing velocity dispersion with cosmic time, the increase being faster for satellites. The analysis shows that proto-clusters are very extended regions, $\gtrsim 20 \, Mpc$ at $z \gtrsim 1$. The fraction of galaxies in proto-cluster regions that are not progenitor of cluster galaxies varies with redshift, stellar mass and area considered. It is about 20-30 per cent for galaxies with stellar mass $\sim 10^9\,{\rm M}_{\sun}$, while negligible for the most massive galaxies considered. Nevertheless, these objects have properties similar to those of progenitors. We investigate the building-up of the passive-sequence in clusters, and find that their progenitors are on average always active at any redshift of interest of proto-clusters. The main mechanism which quenches their star formation is the removal of the hot gas reservoir at the time of accretion. The later galaxies are accreted (become satellite), and the more the cold gas available, the longer the time spent as active. Central galaxies are active over all redshift range considered, although a non-negligible fraction of them become passive at redshift $z<1$, due to strong feedback from Active Galactic Nuclei.

preprint2014arXiv

Galaxy alignment on large and small scales

Galaxies are not randomly distributed across the universe but showing different kinds of alignment on different scales. On small scales satellite galaxies have a tendency to distribute along the major axis of the central galaxy, with dependence on galaxy properties that both red satellites and centrals have stronger alignment than their blue counterparts. On large scales, it is found that the major axes of Luminous Red Galaxies (LRGs) have correlation up to 30Mpc/h. Using hydro-dynamical simulation with star formation, we investigate the origin of galaxy alignment on different scales. It is found that most red satellite galaxies stay in the inner region of dark matter halo inside which the shape of central galaxy is well aligned with the dark matter distribution. Red centrals have stronger alignment than blue ones as they live in massive haloes and the central galaxy-halo alignment increases with halo mass. On large scales, the alignment of LRGs is also from the galaxy-halo shape correlation, but with some extent of mis-alignment. The massive haloes have stronger alignment than haloes in filament which connect massive haloes. This is contrary to the naive expectation that cosmic filament is the cause of halo alignment.

preprint2014arXiv

Radiation fields by intermediate-age stellar populations with binaries as ionizing sources of H II regions

Radiation fields emitted by O, B-type stars or young stellar populations (SPs) are generally considered as significant central ionizing sources (CISs) of classic H II regions. In our previous studies, we show that the inclusion of binary interactions in stellar population synthesis models can significantly increase the ultraviolet spectrum hardness and the number of ionizing photons of intermediate-age (IA, 7<~log(t/yr)<~8) SPs. In this work, we present photoionization models of H II regions ionized by radiation fields emitted by IA SPs, and show that radiation fields of IA SPs, including binary systems, are in theory possible candidates of significant CISs of classic H II regions. When radiation fields of IA SPs comprising binary systems are used as the CISs of classic H II regions, the theoretical strengths of a number of lines (such as [O III]4959', [S II]6716', etc.), weaker than observations, are raised; the border /selection-criterion lines, between star-forming galaxies and AGNs in the diagnostic diagrams (for example, [N II]6583/Ha~[O III]5007/Hb), move into the region occupied originally by AGNs; and He II 1640 line, observed in Lyman break and high-redshift gravitationally lensed galaxies, also can be produced.

preprint2013arXiv

Multipole Gravitational Lensing and High-order Perturbations on the Quadrupole Lens

An arbitrary surface mass density of gravitational lens can be decomposed into multipole components. We simulate the ray-tracing for the multipolar mass distribution of generalized SIS (Singular Isothermal Sphere) model, based on the deflection angles which are analytically calculated. The magnification patterns in the source plane are then derived from inverse shooting technique. As have been found, the caustics of odd mode lenses are composed of two overlapping layers for some lens models. When a point source traverses such kind of overlapping caustics, the image numbers change by \pm 4, rather than \pm 2. There are two kinds of images for the caustics. One is the critical curve and the other is the transition locus. It is found that the image number of the fold is exactly the average value of image numbers on two sides of the fold, while the image number of the cusp is equal to the smaller one. We also focus on the magnification patterns of the quadrupole (m = 2) lenses under the perturbations of m = 3, 4 and 5 mode components, and found that one, two, and three butterfly or swallowtail singularities can be produced respectively. With the increasing intensity of the high-order perturbations, the singularities grow up to bring sixfold image regions. If these perturbations are large enough to let two or three of the butterflies or swallowtails contact, eightfold or tenfold image regions can be produced as well. The possible astronomical applications are discussed.

preprint2012arXiv

Exploring Galaxy Formation Models and Cosmologies with Galaxy Clustering

Using N-body simulations and galaxy formation models, we study the galaxy stellar mass correlation and the two-point auto-correlation. The simulations are run with cosmological parameters from the WMAP first, third and seven year results, which mainly differ in the perturbation amplitude of σ_{8}. The stellar mass of galaxies are determined using either a semi-analytical galaxy formation model or a simple empirical abundance matching method. Compared to the SDSS DR7 data at z=0 and the DEEP2 results at z=1, we find that the predicted galaxy clusterings from the semi-analytical model are higher than the data at small scales, regardless of the adopted cosmology. Conversely, the abundance matching method predicts good agreement with the data at both z=0 and z=1 for high σ_8 cosmologies (WMAP1 & WMAP7), but the predictions from a low σ_8 cosmology (WMAP3) are significantly lower than the data at z=0. We find that the excess clustering at small-scales in the semi-analytical model mainly arises from satellites in massive haloes, indicating that either the star formation is too efficient in low-mass haloes or tidal stripping is too inefficient at high redshift. Our results show that galaxy clustering is strongly affected by the models for galaxy formation, thus can be used to constrain the baryonic physics. The weak dependence of galaxy clustering on cosmological parameters makes it difficult to constrain the WMAP1 and WMAP7 cosmologies.

preprint2012arXiv

Yunnan-III models for Evolutionary population synthesis

We build the Yunnan-III evolutionary population synthesis (EPS) models by using the MESA stellar evolution code, BaSeL stellar spectra library and the initial mass functions (IMFs) of Kroupa and Salpeter, and present colours and integrated spectral energy distributions (ISEDs) of solar-metallicity stellar populations (SPs) in the range of 1Myr-15 Gyr. The main characteristic of the Yunnan-III EPS models is the usage of a set of self-consistent solar-metallicity stellar evolutionary tracks (the masses of stars are from 0.1 to 100Msun). This set of tracks is obtained by using the state-of-the-art MESA code. MESA code can evolve stellar models through thermally pulsing asymptotic giant branch (TP-AGB) phase for low- and intermediate-mass stars. By comparisons, we confirm that the inclusion of TP-AGB stars make the V-K, V-J and V-R colours of SPs redder and the infrared flux larger at ages log(t/yr)>7.6 (the differences reach the maximum at log(t/yr)~8.6, ~0.5-0.2mag for colours, ~2 times for K-band flux). The stellar evolutionary tracks, isochrones, colours and ISEDs can be obtained on request from the first author or from our website (http://www1.ynao.ac.cn/~zhangfh/). Using the isochrones, you can build your EPS models. Now the format of stellar evolutionary tracks is the same as that in the STARBURST99 code, you can put them into the STARBURST99 code and get the SP's results. Moreover, the colours involving other passbands or on other systems (for example, HST $F439W-F555W$ colour on AB system) can also be obtained on request.

preprint2011arXiv

Binary interactions on the calibrations of star formation rate

Using the evolutionary population synthesis (EPS) models with and without binary interactions (BIs), we present L_Ha, L_[OII], L_{i,UV} and L_FIR for Burst, S0, Sa-Sd and Irr galaxies, and present the calibrations of star formation rate (SFR) in terms of these diagnostics. By comparison, we find that BIs lower the SFR.vs.L_Ha and SFR.vs.L_[OII] conversion factors by ~0.2 dex. Moreover, BIs do not significantly vary the calibrations of SFR in terms of L_{i,UV}. In addition, BIs have little effect on the flux at 2800A. At last, the calibration of SFR from L_FIR is almost unaffected by BIs. This is caused by the fact that BIs almost do not affect the bolometric magnitudes of SPs. We also discuss the effects of initial mass function (IMF), gas-recycle assumption (GRA) and EPS models on these SFR calibrations. Comparing the results by using Salpeter (S55) IMF with those by using Miller & Scalo (MS79) IMF, we find that the SFR.vs.L_Ha and SFR.vs.L_[OII] conversion factors by using S55 IMF are greater by 0.4 and 0.2 dex than those by using MS79 IMF for the models with and without BIs, respectively. The SFR.vs.L_{i,UV} and SFR.vs.L_FIR conversion factors by using S55 IMF are larger by an amount of 0.2 dex than the corresponding ones by using MS79 IMF. The inclusion of GRA only lowers these SFR calibrations at faint SFR. Moreover, comparing the results when using different EPS models, we find that the differences in the SFR.vs.L_Ha and SFR.vs.L_[OII] conversion factors reach ~ 0.7 and 0.9 dex, the difference in the SFR.vs.L_FIR conversion factor reaches 0.4 and 0.8 dex, and the differences in the SFR.vs.L_{i,UV} conversion factors reach 0.3 and 0.2 dex when using S55 and NON-S55 IMF (including Cha03, K01, K93' and MS79 IMFs, partly caused by the difference in the IMF), respectively. At last, we give the conversion coefficients between SFR and these diagnostics for all models.

preprint2009arXiv

Properties of Long Gamma-Ray Burst Host Galaxies in Cosmological Simulation

We use galaxy catalogues constructed by combining high-resolution N-body simulations with semi-analytic models of galaxy formation to study the properties of Long Gamma-Ray Burst (LGRB) host galaxies. We assume that LGRBs originate from the death of massive young stars and analyse how results are affected by different metallicity constraints on the progenitor stars. As expected, the host sample with no metallicity restriction on the progenitor stars provides a perfect tracer of the cosmic star formation history. When LGRBs are required to be generated by low-metallicity stars, they trace a decreasing fraction of the cosmic star formation rate at lower redshift, as a consequence of the global increase in metallicity. We study the properties of host galaxies up to high redshift (~9), finding that they typically have low-metallicity (Z<0.5 Z_sun) and that they are small (M<10^9 M_sun), bluer and younger than the average galaxy population, in agreement with observational data. They are also less clustered than typical L_* galaxies in the Universe, and their descendents are massive, red and reside in groups of galaxies with halo mass between 10^{13} M_sun to 10^{14} M_sun.

preprint2008arXiv

New Constraints on the Efficiencies of Ram-Pressure Stripping and the Tidal Disruption of Satellite Galaxies

Using data from the Sloan Digital Sky Survey (SDSS) it has recently been shown that the red fraction of satellite galaxies increases with stellar mass. Semi-analytical models, however, predict red satellite fractions that are independent of stellar mass, and much higher than observed. It has been argued that this discrepancy owes to the fact that the models assume that satellite galaxies are instantaneously stripped of their hot gas reservoirs at the moment they are accreted into a bigger halo. In this letter we show that the fraction of red satellites can be brought in better agreement with the data by simply decreasing this stripping efficiency. However, this also results in a red fraction of massive centrals that is much too low. This owes to the fact that the massive centrals now accrete satellite galaxies that are bluer and more gas-rich. However, if a significant fraction of low mass satellite galaxies is tidally disrupted before being accreted by their central host galaxy, as suggested by recent studies, the red fractions of both centrals and satellites can be reproduced reasonably well. A problem remains with the red fraction of centrals of intermediate mass, which is likely to reflect an oversimplified treatment of AGN feedback.

preprint2006arXiv

Smoothing Algorithms and High-order Singularities in Gravitational Lensing

We propose a new smoothing method for obtaining surface densities from discrete particle positions from numerical simulations. This is an essential step for many applications in gravitational lensing. This method is based on the ``scatter'' interpretation of the discrete density field in the Smoothed Particle Hydrodynamics. We use Monte Carlo simulations of uniform density fields and one isothermal ellipsoid to empirically derive the noise properties, and best smoothing parameters (such as the number of nearest neighbors used). A cluster from high-resolution simulations is then used to assess the reality of high-order singularities such as swallowtails and butterflies in caustics, which are important for the interpretation of substructures in gravitational lenses. We also compare our method with the Delaunay tesselation field estimator using the galaxy studied by Bradac et al. (2004), and find good agreements. We show that higher order singularities are not only connected with bound subhaloes but also with the satellite streams. However, the presence of high-order singularities are sensitive to not only the fluctuation amplitude of the surface density, but also the detailed form of the underlying smooth lensing potential (such as ellipticity and external shear).