There are expected to be physical relationships between the globular clusters (GCs) and stellar substructures in the Milky Way, not all of which have yet been found. We search for such substructures from a combined halo sample of SDSS blue horizontal-branch and SDSS+LAMOST RR Lyrae stars, cross-matched with astrometric information from $Gaia$ DR2. This is a sample of old stars which are also excellent tracers of structures, ideal for searching for ancient relics in the outer stellar halo. By applying the neural-network-based method StarGO to the full 4D dynamical space of our sample, we rediscover the Sagittarius Stream, and find the debris of the $Gaia$-Enceladus-Sausage (GES) and the Sequoia events in the outer halo, as well as their linkages with several GCs. Most importantly, we find a new, low-mass, debris stream associated with a pair of GCs (NGC 5024 and NGC 5053), which we dub LMS-1. This stream has a very polar orbit, and occupies a region between 10 to 20 kpc from the Galactic center. NGC 5024 (M53), the more-massive of the associated GC pair, is very likely the nuclear star cluster of a now-disrupted dwarf galaxy progenitor, based on the results from N-body simulations.

Binary neutron star mergers (BNSMs) is currently the most promising source of \textsl{r}-process thanks to the detection of GW170817. The estimated occurring frequency and the amount of mass ejected per merger indicate that BNSMs by itself can account for all the \textsl{r}-process enrichment in the Galaxy. However, the decreasing trend of [Eu/Fe] vs [Fe/H] of disk stars for [Fe/H]$\gtrsim -1$ in the solar neighborhood is inconsistent with the flat trend expected from BNSMs with a standard delay time distribution (DTD) $\propto t^{-1}$. This has led to the suggestion that either additional sources or modification to the DTD of BNSMs is required to match the observations. We investigate the effects of natal kicks received during the birth of neutron star binaries on the chemical evolution of \textsl{r}-process element Eu in the Milky Way by combining the results from the galactic dynamics code \textsc{galpy} with a one-zone Galactic chemical evolution model \textsc{omega}. We show that when key inputs from simulations of the inside-out disk evolution are combined with natal kicks, BNSMs can naturally reproduce the observed decreasing trend of [Eu/Fe] with [Fe/H] in the solar neighborhood without the need for modification to the DTD or additional \textsl{r}-process sources.

We report on the spectroscopic analysis of RAVE J183013.5-455510, an extremely metal-poor star, highly enhanced in CNO, and with discernible contributions from the rapid neutron-capture process. There is no evidence of binarity for this object. At [Fe/H]=-3.57, this is one of the lowest metallicity stars currently observed, with 18 measured abundances of neutron-capture elements. The presence of Ba, La, and Ce abundances above the Solar System r-process predictions suggest that there must have been a non-standard source of r-process elements operating at such low metallicities. One plausible explanation is that this enhancement originates from material ejected at unusually fast velocities in a neutron star merger event. We also explore the possibility that the neutron-capture elements were produced during the evolution and explosion of a rotating massive star. In addition, based on comparisons with yields from zero-metallicity faint supernova, we speculate that RAVE J1830-4555 was formed from a gas cloud pre-enriched by both progenitor types. From analysis based on Gaia DR2 measurements, we show that this star has orbital properties similar to the Galactic metal-weak thick-disk stellar population.

We propose a flexible method for estimating luminosity functions (LFs) based on kernel density estimation (KDE), the most popular nonparametric density estimation approach developed in modern statistics, to overcome issues surrounding binning of LFs. One challenge in applying KDE to LFs is how to treat the boundary bias problem, since astronomical surveys usually obtain truncated samples predominantly due to the flux-density limits of surveys. We use two solutions, the transformation KDE method ($\hat{\phi}_{\mathrm{t}}$), and the transformation-reflection KDE method ($\hat{\phi}_{\mathrm{tr}}$) to reduce the boundary bias. We develop a new likelihood cross-validation criterion for selecting optimal bandwidths, based on which, the posterior probability distribution of bandwidth and transformation parameters for $\hat{\phi}_{\mathrm{t}}$ and $\hat{\phi}_{\mathrm{tr}}$ are derived within a Markov chain Monte Carlo (MCMC) sampling procedure. The simulation result shows that $\hat{\phi}_{\mathrm{t}}$ and $\hat{\phi}_{\mathrm{tr}}$ perform better than the traditional binned method, especially in the sparse data regime around the flux-limit of a survey or at the bright-end of the LF. To further improve the performance of our KDE methods, we develop the transformation-reflection adaptive KDE approach ($\hat{\phi}_{\mathrm{tra}}$). Monte Carlo simulations suggest that it has a good stability and reliability in performance, and is around an order of magnitude more accurate than using the binned method. By applying our adaptive KDE method to a quasar sample, we find that it achieves estimates comparable to the rigorous determination by a previous work, while making far fewer assumptions about the LF. The KDE method we develop has the advantages of both parametric and non-parametric methods.

A large number of new members ($\sim$150) of the Cetus Stream (CS) were identified from their clustering features in dynamical space using 6D kinematic data by combining LAMOST DR5 and Gaia DR2 surveys. They map a diffuse structure that extends over at least 100 degrees in the northern and southern Galactic hemispheres, at heliocentric distances between 20 to 50 kpc. Taking advantage of this expanded dataset, we model the stream with a suite of tailored N-body simulations. Our findings exclude the possibility that the NGC 5824 globular cluster is the core of the progenitor of the stream, as postulated by previous studies. Our best models, which successfully reproduce the features of the CS indicate that the progenitor is likely a dwarf galaxy of $\sim$ 2$\times$10$^9$M$_{\odot}$, with a diffuse disc morphology. The merger occured $\sim$ 5 Gyr ago and since then it has experienced approximately eight apo-center passages. Our results suggest that NGC 5824 was either a globular cluster situated off-centre in the dwarf progenitor or, alternatively, it was the nuclear star cluster of another dwarf galaxy that has very similar orbit as the progenitor of the CS. In both scenarios, the progenitor systems would leave streams around NGC 5824, but with distinct distance distributions. To discriminate between these scenarios, the detection and accurate distance measurements of the predicted stream around the GC are crucial, which will be possible in the upcoming LSST era. Our simulations also predict that part of the Southern Cetus stream is very likely the newly discovered Palca stream, and possibly related to another, more diffuse Southern substructure, the Eridanus-Pheonix overdensity.

We search for dynamical substructures in the LAMOST DR3 very metal-poor (VMP) star catalog. After cross-matching with Gaia DR2, there are 3300 VMP stars with available high-quality astrometric information that have halo-like kinematics. We apply a method based on self-organizing maps to find groups clustered in the 4D space of orbital energy and angular momentum. We identify 57 dynamically tagged groups, which we label DTG-1 to DTG-57. Most of them belong to existing substructures in the nearby halo, such as the $Gaia$ Sausage or Sequoia. The stream identified by Helmi et al. is recovered, but the two disjoint portions of the substructure have distinct dynamical properties. The very retrograde substructure Rg5 found previously by Myeong et al. is also retrieved. We report 6 new DTGs with highly retrograde orbits, 2 with very prograde orbits, and 12 with polar orbits. By mapping other datasets (APOGEE halo stars, and catalogs of r-process-enhanced and CEMP stars) onto the trained neuron map, we can associate stars with detailed chemical abundances to the DTGs, and look for associations with chemically peculiar stars. The highly eccentric $Gaia$ Sausage groups contain representatives both of debris from the satellite itself (which is $\alpha$-poor) and the Splashed Disk, sent up into eccentric halo orbits from the encounter (and is $\alpha$-rich). The new prograde substructures also appear to be associated with the Splashed Disk. The DTGs belonging to the $Gaia$ Sausage host two relatively metal-rich $r$-II stars and six CEMP stars in different sub-classes, consistent with the idea that the $Gaia$ Sausage progenitor is a massive dwarf galaxy. Rg5 is dynamically associated with two highly $r$-process-enhanced stars with [Fe/H] $\sim -$3. This finding indicates that its progenitor might be an ultra-faint dwarf galaxy that has experienced $r$-process enrichment from neutron star mergers.

We use a novel cluster identification tool StarGO to explore the metal poor ([Fe/H] $<$ -1.5) outer stellar halo (d $>$ 15 kpc) of the Milky Way using data from Gaia, LAMOST and SDSS. Our method is built using an unsupervised learning algorithm, a self-organizing map, which trains a 2-D neural network to learn the topological structures of a data set from an n-D input space. Using a 4-D space of angular momentum and orbital energy, we identify three distinct groups corresponding to the Sagittarius, Orphan, and Cetus Streams. For the first time we are able to discover a northern counterpart to the Cetus stream. We test the robustness of this new detection using mock data and find that the significance is more than 5-sigma. We also find that the existing southern counterpart bifurcates into two clumps with different radial velocities. By exploiting the visualization power of StarGO, we attach MW globular clusters to the same trained neural network. The Sagittarius stream is found to have five related clusters, confirming recent literature studies, and the Cetus stream has one associated cluster, NGC 5824. This latter association has previously been postulated, but can only now be truly confirmed thanks to the high-precision Gaia proper motions and large numbers of stellar spectra from LAMOST. The large metallicity dispersion of the stream indicates that the progenitor cannot be a globular cluster. Given the mean metallicity of the stream, we propose that the stream is the result of a merger of a low-mass dwarf galaxy that hosted a large nuclear star cluster (NGC 5824).

The dense clusters within the Serpens Molecular Cloud are among the most active regions of nearby star formation. In this paper, we use Gaia DR2 parallaxes and proper motions to statistically measure 1167 kinematic members of Serpens, few of which were previously identified, to evaluate the star formation history of the complex. The optical members of Serpens are concentrated in three distinct groups located at 380 to 480 pc; the densest clusters are still highly obscured by optically-thick dust and have few optical members. The total population of young stars and protostars in Serpens is at least 2000 stars, including past surveys that were most sensitive to protostars and disks, and may be far higher. Distances to dark clouds measured from deficits in star counts are consistent with the distances to the optical star clusters. The Serpens Molecular Cloud is seen in the foreground of the Aquila Rift, dark clouds located at 600 to 700 pc, and behind patchy extinction, here called the Serpens Cirrus, located at ~250 pc. Based on the lack of a distributed population of older stars, the star formation rate throughout the Serpens Molecular Cloud increased by at least a factor of 20 within the past ~5 Myr. The optically bright stars in Serpens Northeast are visible because their natal molecular cloud has been eroded and not because they were flung outwards from a central factory of star formation. The separation between subclusters of 20 to 100 pc and the absence of an older population leads to speculation that an external forcing was needed to trigger the active star formation.

We report the discovery of tidal structures around the intermediate-aged ($\sim$ 700--800~Myr), nearby ($\sim85$~pc) star cluster Coma Berenices. The spatial and kinematic grouping of stars is determined with the {\it Gaia} DR2 parallax and proper motion data, by a clustering analysis tool, \textsc{StarGO}, to map 5D parameters ($X, Y, Z$, $\mu_\alpha \cos\delta, \mu_\delta$) onto a 2D neural network. A leading and a trailing tails, each with an extension of $\sim50$~pc are revealed for the first time around this disrupting star cluster. The cluster members, totaling to $\sim115^{+5}_{-3}\,\rm {M_\odot}$, are clearly mass segregated, and exhibit a flat mass function with $\alpha \sim 0.79\pm0.16$, in the sense of $dN/dm \propto m^{-\alpha}$, where $N$ is the number of member stars and $m$ is stellar mass, in the mass range of $m=0.25$--$2.51~{\rm M_\odot}$. Within the tidal radius of $\sim$6.9~pc, there are 77 member candidates with an average position, i.e., as the cluster center, of R.A.= 186.8110~deg, and decl.= 25.8112~deg, and an average distance of 85.8~pc. Additional 120 member candidates reside in the tidal structures, i.e., outnumbering those in the cluster core. The expansion of escaping members lead to an anisotropy in the velocity field of the tidal tails. Our analysis also serendipitously uncovers an adjacent stellar group, part of which has been cataloged in the literature. We identify 218 member candidates, 10 times more than previously known. This star group is some 65~pc away from, and $\sim400$~Myr younger than, Coma Ber, but is already at the final stage of disruption.

We use a sample of ~3000 galaxies from the MaNGA MPL-7 internal data release to study the alpha abundance distribution within low-redshift galaxies. We use the Lick index ratio Mgb/<Fe> as an alpha abundance indicator to study relationships between the alpha abundance distribution and galaxy properties such as effective stellar velocity dispersion within 0.3 effective radii (sigma_*), galaxy environment, and dark matter halo formation time (z_f). We find that (1) all galaxies show a tight correlation between Mgb/<Fe> and sigma_*; (2) `old' (H_beta < 3) low-sigma_* galaxies in high local density environment and inner regions within galaxy groups are enhanced in Mgb/<Fe>, while `young' (H_beta>3) galaxies and high-mass galaxies show no or less environmental dependence; (3) `old' galaxies with high-z_f show enhanced Mgb/<Fe> over low- and medium-z_f; (4) Mgb/<Fe> gradients are close to zero and show dependence on sigma_* but no obvious dependence on the environment or z_f. Our study indicates that stellar velocity dispersion or galaxy mass is the main parameter driving the Mgb/<Fe> enhancement, although environments appear to have modest effects, particularly for low- and medium-mass galaxies.

The radio luminosity functions (RLFs) of active galactic nuclei (AGNs) are traditionally measured based on total emission, which doesn't reflect the current activity of the central black hole. The increasing interest in compact radio cores of AGNs requires determination of the RLF based on core emission (i.e., core RLF). In this work we have established a large sample (totaling 1207) of radio-loud AGNs, mainly consisting of radio galaxies (RGs) and steep-spectrum radio quasars (SSRQs). Based on the sample, we explore the relationship between core luminosity ($L_c$) and total luminosity ($L_t$) via a powerful statistical tool called "Copula". The conditional probability distribution $p(\log L_{c} \mid \log L_{t})$ is obtained. We derive the core RLF as a convolution of $p(\log L_{c} \mid \log L_{t})$ with the total RLF which was determined by previous work. We relate the separate RG and SSRQ core RLFs via a relativistic beaming model and find that SSRQs have an average Lorentz factor of $\gamma=9.84_{-2.50}^{+3.61}$, and that most are seen within $8^{\circ} \lesssim \theta \lesssim 45^{\circ}$ of the jet axis. Compared with the total RLF which is mainly contributed by extended emission, the core RLF shows a very weak luminosity-dependent evolution, with the number density peaking around $z\thicksim 0.8$ for all luminosities. Differences between core and total RLFs can be explained in a framework involving a combination of density and luminosity evolutions where the cores have significantly weaker luminosity evolution than the extended emission.

Following previous work, we further confirm that the cosmic evolution of steep-spectrum radio-loud AGNs (active galactic nuclei) can be reproduced by a simple combination of density evolution (DE) and luminosity evolution (LE). This mixture evolution scenario can naturally explain the luminosity-dependent evolution of radio-loud AGNs. Our models successfully fitted a large amount of data on radio luminosity functions (RLFs) of steep-spectrum sources and multi-frequency source counts. The modeling indicates that the DE slowly increase as $(1+z)^{0.3 \thicksim 1.3}$ out to $z \thicksim 0.8$, and then rapidly decreases as $(1+z)^{-6.8 \thicksim -5.7}$, while the LE rapidly increase as $(1+z)^{4.8}$ out to a higher redshift (at least $z>3.5$). We find a high-redshift decline (i.e. redshift cutoff) in the number density of steep-spectrum radio sources, but we cannot conclude whether such decline is sharp or shallow. We believe that whether a redshift cutoff occurs or not depends mainly on DE, while its steepness is decided by LE, which, however, cannot be well constrained due to the lack of high-redshift samples. Most intriguingly, according to our mixture evolution scenario, there appears to be no need for different evolution for the low- and high-power radio-loud AGNs. Both types of sources experience the same combined evolution of DE and LE.

We develop a new method StarGO (Stars' Galactic Origin) to identify the galactic origins of halo stars using their kinematics. Our method is based on self-organizing map (SOM), which is one of the most popular unsupervised learning algorithms. StarGO combines SOM with a novel adaptive group identification algorithm with essentially no free parameters. In order to evaluate our model, we build a synthetic stellar halo from mergers of nine satellites in the Milky Way. We construct the mock catalogue by extracting a heliocentric volume of 10 kpc from our simulations and assigning expected observational uncertainties corresponding to bright stars from Gaia DR2 and LAMOST DR5. We compare the results from StarGO against that from a Friends-of-Friends (FoF) based method in the space of orbital energy and angular momentum. We show that StarGO is able to systematically identify more satellites and achieve higher number fraction of identified stars for most of the satellites within the extracted volume. When applied to data from Gaia DR2, StarGO will enable us to reveal the origins of the inner stellar halo in unprecedented detail.

We fit the spectral energy distributions (SEDs) of 46 GeV - TeV BL Lac objects in the frame of leptonic one-zone synchrotron self-Compton (SSC) model and investigate the physical properties of these objects. We use the Markov Chain Monte Carlo (MCMC) method to obtain the basic parameters, such as magnetic field (B), the break energy of the relativistic electron distribution ($\gamma'_{\rm{b}}$) and the electron energy spectral index. Based on the modeling results, we support the following scenarios on GeV-TeV BL Lac objects: (1) Some sources have large Doppler factors, implying other radiation mechanism should be considered. (2) Comparing with FSRQs, GeV-TeV BL Lac objects have weaker magnetic field and larger Doppler factor, which cause the ineffective cooling and shift the SEDs to higher bands. Their jet powers are around $4.0\times 10^{45}~\rm{ erg\cdot s}^{-1}$, comparing with radiation power, $5.0\times 10^{42}~\rm{ erg\cdot s}^{-1}$, indicating that only a small fraction of jet power is transformed into the emission power. (3) For some BL Lacs with large Doppler factors, their jet components could have two substructures, e.g., the fast core and the slow sheath. For most GeV-TeV BL Lacs, Kelvin-Helmholtz instabilities are suppressed by their higher magnetic fields, leading few micro-variability or intro-day variability in the optical bands. (4) Combined with a sample of FSRQs, an anti-correlation between the peak luminosity $L_{\rm {pk}}$ and the peak frequency $\nu_{\rm {pk}}$ is obtained, favoring the blazar sequence scenario. In addition, an anti-correlation between the jet power $P_{\rm {jet}}$ and the break Lorentz factor $\gamma_{\rm {b}}$ also supports the blazar sequence.

We show that the mass of a dark matter halo can be inferred from the dynamical status of its satellite galaxies. Using 9 dark-matter simulations of halos like the Milky Way (MW), we find that the present-day substructures in each halo follow a characteristic distribution in the phase space of orbital binding energy and angular momentum, and that this distribution is similar from halo to halo but has an intrinsic dependence on the halo formation history. We construct this distribution directly from the simulations for a specific halo and extend the result to halos of similar formation history but different masses by scaling. The mass of an observed halo can then be estimated by maximizing the likelihood in comparing the measured kinematic parameters of its satellite galaxies with these distributions. We test the validity and accuracy of this method with mock samples taken from the simulations. Using the positions, radial velocities, and proper motions of 9 tracers and assuming observational uncertainties comparable to those of MW satellite galaxies, we find that the halo mass can be recovered to within $\sim$40%. The accuracy can be improved to within $\sim$25% if 30 tracers are used. However, the dependence of the phase-space distribution on the halo formation history sets a minimum uncertainty of $\sim$20% that cannot be reduced by using more tracers. We believe that this minimum uncertainty also applies to any mass determination for a halo when the phase space information of other kinematic tracers is used.

We present a model-dependent method to estimate the redshifts of three TeV BL Lac objects (BL Lacs) through fltting their (quasi-) simultaneous multi-waveband spectral energy distributions (SEDs) by one-zone leptonic synchrotron self-Compton (SSC) model. Considering the impact of electron energy distributions (EEDs) on the results, we use three types of EEDs, such as the power-law EED with exponential cut-ofi (PLC), the log-parabola (PLLP) and the broken power-law (BPL) EEDs, to flt the SEDs. We also use a parameter fi to describe the uncertainties of the extragalactic background light (EBL) models, as in (Abdo et al. 2010b). We then use Markov Chain Monte Carlo (MCMC) method to explore multi-dimensional parameter space and obtain the uncertainties of the model parameters based on the observational data. We apply our method to obtain the redshifts of three TeV BL Lac objects in the marginalized 68% confldence, and flnd that the PLC EED does not flt the SEDs. For 3C66A, the redshift is 0.14 - 0.31 and 0.16 - 0.32 in the BPL and PLLP EEDs; for PKS1424+240, the redshift is 0.55 - 0.68 and 0.55 - 0.67 in the BPL and PLLP EEDs; for PG1553+113, the redshift is 0.22 - 0.48 and 0.22 - 0.39 in the BPL and PLLP EEDs. We also estimate the redshift of PKS1424+240 in the high stage to be 0.46 - 0.67 in the PLLP EED, roughly consistent with that in the low stage.

We report the observations of an electron vortex magnetic hole corresponding to a new type of coherent structures in the magnetosheath turbulent plasma using the Magnetospheric Multiscale (MMS) mission data. The magnetic hole is characterized by a magnetic depression, a density peak, a total electron temperature increase (with a parallel temperature decrease but a perpendicular temperature increase), and strong currents carried by the electrons. The current has a dip in the center of the magnetic hole and a peak in the outer region of the magnetic hole. The estimated size of the magnetic hole is about 0.23 \r{ho}i (~ 30 \r{ho}e) in the circular cross-section perpendicular to its axis, where \r{ho}i and \r{ho}e are respectively the proton and electron gyroradius. There are no clear enhancement seen in high energy electron fluxes, but an enhancement in the perpendicular electron fluxes at ~ 90{\deg} pitch angles inside the magnetic hole is seen, implying that the electron are trapped within it. The variations of the electron velocity components Vem and Ven suggest that an electron vortex is formed by trapping electrons inside the magnetic hole in the circular cross-section (in the M-N plane). These observations demonstrate the existence of a new type of coherent structures behaving as an electron vortex magnetic hole in turbulent space plasmas as predicted by recent kinetic simulations.

In the solar wind, power spectral density (PSD) of the magnetic field fluctuations generally follow the so-called Kolmogorov spectrum f^-5/3 in the inertial range, where the dynamics is thought to be dominated by nonlinear interactions between counter-propagating incompressible Alfv\'en wave parquets. These features are thought to be ubiquitous in space plasmas. The present study gives a new and more complex picture of magnetohydrodynamics (MHD) turbulence as observed in the terrestrial magnetosheath. The study uses three years of in-situ data from the Cluster mission to explore the nature of the magnetic fluctuations at MHD scales in different locations within the magnetosheath, including flanks and subsolar regions. It is found that the magnetic field fluctuations at MHD scales generally have a PSD close to f^-1 (shallower than the Kolmogorov one f^-5/3) down to the ion characteristic scale, which recalls the energy containing scales of solar wind turbulence. The Kolmogorov spectrum is observed only away from the bow shock toward the flank and the magnetopause regions in 17% of the analyzed time intervals. Measuring the magnetic compressibility, it is shown that only a fraction (35%) of the observed Kolmogorov spectra were populated by shear Alfv\'enic fluctuations, whereas the majority of the events (65%) was found to be dominated by compressible magnetosonic-like fluctuations, which contrasts with well-known turbulence properties in the solar wind. This study gives a first comprehensive view of the origin of the f^-1 and the transition to the Kolmogorov inertial range; both questions remain controversial in solar wind turbulence.

In this paper, we scrutinize the effect of spectral index distribution on estimating the AGN (active galactic nucleus) radio luminosity function (RLF) by a Monte Carlo method. We find that the traditional bivariate RLF estimators can cause bias in varying degree. The bias is especially pronounced for the flat-spectrum radio sources whose spectral index distribution is more scattered. We believe that the bias is caused because the $K$-corrections complicate the truncation boundary on the $L-z$ plane of the sample, but the traditional bivariate RLF estimators have difficulty in dealing with this boundary condition properly. We suggest that the spectral index distribution should be incorporated into the RLF analysis process to obtain a robust estimation. This drives the need for a trivariate function of the form $\Phi(\alpha,z,L)$ which we show provides an accurate basis for measuring the RLF.

By cross-matching the currently largest optical catalog of galaxy clusters and the NVSS radio survey database, we obtain the largest complete sample of brightest cluster galaxies (BCGs) in the redshift range of 0.05<z<0.45, which have radio emission and redshift information. We confirm that more powerful radio BCGs tend to be these optically very bight galaxies located in more relaxed clusters. We derived the radio luminosity functions of BCGs from the largest complete sample of BCGs, and find that the functions depend on the optical luminosity of BCGs and the dynamical state of galaxy clusters. However, the radio luminosity function does not show significant evolution with redshift.

We propose a mixture evolution scenario to model the evolution of the radio luminosity function (RLF) of steep spectrum AGNs (active galactic nuclei) based on a Bayesian method. In this scenario, the shape of the RLF is determined by both the density and luminosity evolution. Our models indicate that the density evolution is positive until a redshift of $\thicksim 0.9$, at which point it becomes negative, while the luminosity evolution is positive to a higher redshift ($z \thicksim 5$ for model B and $z \thicksim 3.5$ for model C), where it becomes negative. Our mixture evolution model works well, and the modeled RLFs are in good agreement with previous determinations. The mixture evolution scenario can naturally explain the luminosity-dependent evolution of the RLFs.

Fornax is the brightest Milky Way (MW) dwarf spheroidal galaxy and its star formation history (SFH) has been derived from observations. We estimate the time evolution of its gas mass and net inflow and outflow rates from the SFH using a simple star formation law that relates the star formation rate to the gas mass. We present a chemical evolution model on a 2D mass grid with supernovae (SNe) as sources of metal enrichment. We find that a key parameter controlling the enrichment is the mass M_x of the gas to mix with the ejecta from each SN. The choice of M_x depends on the evolution of SN remnants and on the global gas dynamics. It differs between the two types of SNe involved and between the periods before and after Fornax became an MW satellite at time t = t_sat . Our results indicate that due to the global gas outflow at t > t_sat , part of the ejecta from each SN may directly escape from Fornax. Sample results from our model are presented and compared with data.

We propose that detailed data on the star formation history of a dwarf spheroidal galaxy (dSph) may be used to estimate the evolution of the total mass M_g(t) for cold gas in its star-forming disk. Using Fornax as an illustrative example, we estimate its M_g(t) and the corresponding net gas flow rate Delta F(t) assuming a global star formation rate psi(t) = lambda_*(t)[M_g(t)/M_sun]^alpha consistent with observations of nearby galaxies. We identify the onset of the transition in Delta F(t) from a net inflow to a net outflow as the time t_sat at which the Fornax halo became a Milky Way satellite and estimate the evolution of its total mass M_h(t) at t < t_sat using the median halo growth history in the LambdaCDM cosmology and its present mass within the half-light radius derived from observations. We examine three different cases of alpha = 1, 1.5, and 2, and justify the corresponding lambda_*(t) by comparing the gas mass fraction f_g(t) = M_g(t)/M_h(t) at t < t_sat with results from simulations of gas accretion by halos in a reionized universe. We find that the Fornax halo grew to M_h(t_sat) ~ 2x10^9 M_sun at t_sat ~ 5 or 8 Gyr, in broad agreement with previous studies using data on its stellar kinematics and its orbital motion. We describe qualitatively the evolution of Fornax as a satellite and discuss potential extension of our approach to other dSphs.

Diffuse radio emission in galaxy clusters is known to be related to cluster mass and cluster dynamical state. We collect the observed fluxes of radio halos, relics, and mini-halos for a sample of galaxy clusters from the literature, and calculate their radio powers. We then obtain the values of cluster mass or mass proxies from previous observations, and also obtain the various dynamical parameters of these galaxy clusters from optical and X-ray data. The radio powers of relics, halos, and mini-halos are correlated with the cluster masses or mass proxies, as found by previous authors, with the correlations concerning giant radio halos being, in general, the strongest ones. We found that the inclusion of dynamical parameters as the third dimension can significantly reduce the data scatter for the scaling relations, especially for radio halos. We therefore conclude that the substructures in X-ray images of galaxy clusters and the irregular distributions of optical brightness of member galaxies can be used to quantitatively characterize the shock waves and turbulence in the intracluster medium responsible for re-accelerating particles to generate the observed diffuse radio emission. The power of radio halos and relics is correlated with cluster mass proxies and dynamical parameters in the form of a fundamental plane.

BL Lacertae (BL Lac) objects at high redshifts ($z\geq 2$) are rarely detected. Through careful analysis of the SDSS spectrum, BZB J1450+5201 is confirmed to be a high-$z$ BL Lac object with $z\geq$ 2.471 by identifying the Ly$\alpha$ 1216 and CIV 1548/1550 absorption lines. This indicates that BZB J1450+5201 is the most distant BL Lac object discovered to date. Careful analysis of the five-year \fermi data of 2FGL J1451.0+5159 shows that its $\gamma$-ray emission is robust with confidence level of 6.2$\sigma$ at 1-3 GeV and 6.7$\sigma$ at 3-10 GeV, and that the confusion of bright neighbor is negligible, which can not be fixed in the analysis of the two-year data. Meanwhile, 2FGL J1451.0+5159 is confirmed to be associated with BZB J1450+5201 using the five-year data. The analysis of multiwavelength data, from radio to $\gamma$-ray energies, indicates BZB J1450+5201 is an intermediate synchrotron peaked (ISP) source and consistent with distributions of other ISP sources at lower redshifts in the second LAT AGN catalog. The pure SSC model seems to be disfavoured, while scattering of weak external emission plus SSC process can provide a satisfactory description of the broadband emission.

We present a first statistical study of subproton and electron scales turbulence in the terrestrial magnetosheath using the Cluster Search Coil Magnetometer (SCM) waveforms of the STAFF instrument measured in the frequency range [1,180] Hz. It is found that clear spectral breaks exist near the electron scale, which separate two power-law like frequency bands referred to as the dispersive and the electron dissipation ranges. The frequencies of the breaks f_b are shown to be well correlated with the electron gyroscale \rho_e rather than with the electron inertial length de. The distribution of the slopes below fb was found to be narrow and peaks near -2.9, while that of the slopes above fb was found broader, peaks near -5.2 and has values as low as -7.5. This is the first time that such steep power-law spectra are reported in space plasma turbulence. These observations provide strong constraints on theoretical modeling of kinetic turbulence and dissipation in collisionless magnetized plasmas.

The classical 1/Va and PC methods of constructing binned luminosity functions (LFs) are revisited and compared by graphical analysis. Using both theoretical analysis and illustration with an example, we show why the two methods give different results for the bins which are crossed by the flux limit curves $L=L_{lim}(z)$. Based on a combined sample simulated by a Monte Carlo method, the estimate $\phi$ of two methods are compared with the input model LFs. The two methods give identical and ideal estimate for the high luminosity points of each redshift interval. However, for the low luminosity bins of all the redshift intervals both methods give smaller estimate than the input model. We conclude that once the LF is evolving with redshift, the classical binned methods will unlikely give an ideal estimate over the total luminosity range. Page & Carrera (2000) noticed that for objects close to the flux limit $\phi_{1/V_{a}}$ nearly always to be too small. We believe this is due to the arbitrary choosing of redshift and luminosity intervals. Because $\phi_{1/V_{a}}$ is more sensitive to how the binning are chosen than $\phi_{PC}$. We suggest a new binning method, which can improve the LFs produced by the 1/Va method significantly, and also improve the LFs produced by the PC methods. Our simulations show that after adopting this new binning, both the 1/Va and PC methods have comparable results.

The work in this paper aims at determining the evolution and possible co-evolution of radio-loud active galactic nuclei (AGNs) and their cores via their radio luminosity functions (i.e., total and core RLF respectively). Using a large combined sample of 1063 radio-loud AGNs selected at low radio frequency, we investigate the radio luminosity function (RLF) at 408 MHz of steep-spectrum radio sources. Our results support a luminosity-dependent evolution. Using core flux density data of the complete sample 3CRR, we investigate the core RLF at 5.0 GHz. Based on the combined sample with incomplete core flux data, we also estimate the core RLF using a modified factor of completeness. Both results are consistent and show that the comoving number density of radio cores displays a persistent decline with redshift, implying a negative density evolution. We find that the core RLF is obviously different from the total RLF at 408 MHz band which is mainly contributed by extended lobes, implying that the cores and extended lobes could not be co-evolving at radio emission.