Afterglows of gamma-ray bursts often show flares, plateaus, and sudden intensity drops: these temporal features are difficult to explain as coming from the forward shock. We calculate radiative properties of early GRB afterglows with the dominant contribution from the reverse shock (RS) propagating in an ultra-relativistic (pulsar-like) wind produced by the long-lasting central engine. RS emission occurs in the fast cooling regime - this ensures high radiative efficiency and allows fast intensity variations. We demonstrate that: (i) mild wind power, of the order of $\sim 10^{46}$ erg s$^{-1}$, can reproduce the afterglows' plateau phase; (ii) termination of the wind can produce sudden steep decays; (iii) mild variations in the wind luminosity can produce short-duration afterglow flares.

Coulomb screening and weak interactions in a hot, magnetized plasma are investigated. Coulomb screening is evaluated in a relativistic thermal plasma in which electrons and positrons are in equilibrium. In addition to temperature effects, effects on weak screening from a strong external magnetic field are evaluated. In high fields, the electron transverse momentum components are quantized into Landau levels. The characteristic plasma screening length at high temperatures and at high magnetic fields is explored. In addition to changes to the screening length, changes in weak interaction rates are estimated. It is found that high fields can result in increased $\beta$-decay rates as the electron and positron spectra are dominated by Landau levels. Finally, the effects studied here are evaluated in a simple r-process model. It is found that relativistic Coulomb screening has a small effect on the final abundance distribution. While changes in weak interaction rates in strong magnetic fields can have an effect on the r-process evolution and abundance distribution, the field strength required to have a significant effect may be larger than what is currently thought to be typical of the r-process environment in collapsar jets or neutron star mergers. If r-process sites exist in fields $\gtrsim 10^{14}$ G effects from fields on weak decays could be significant.

It is well known from the Butcher-Oemler effect that galaxies in dense environment are mostly red with little star formation and the fraction of blue galaxies in galaxy groups/clusters also declines rapidly with redshifts. A recent work by Hashimoto et al. reported a local 'blue cluster' with high fraction of blue galaxies ($\sim 0.57$), higher than the model predictions. They ascribed this blue cluster to the feeding of gas along a filamentary structure around the cluster. In this work we use group catalog from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) and the state-of-art of semi-analytic model (SAM) to investigate the formation of blue clusters in local Universe. In total, we find four blue clusters with halo mass $\sim 10^{14}M_{\odot}$ at $0.02 < z < 0.082$, while only the one found by Hashimoto et al. is in a filamentary structure. The SAM predicts that blue clusters have later formation time and most blue satellite galaxies are recently accreted. We conclude that the formation of blue clusters is mainly governed by newly accreted blue satellites, rather than the effect of large-scale environment.

Combining the LAMOST radial velocities with Gaia parallaxes and proper motions, we presented 3D Galactic space motions and the orbits of 182 single-lined hot subdwarf stars. These stars have been identified by Lei et al. (2020) in Gaia DR2 with LAMOST DR6 and DR7 spectra. He-rich hot subdwarf stars with log(y)>0 show the largest standard deviations of the Galactic velocity components and orbital parameters, while those with -1<log(y)<0 exhibit the second largest standard deviations. The two groups of He-deficient stars with log(y)<-1 show similar standard deviations, which is systematically lower compared to He-rich stars. We also presented a kinematic population classification of the four hot subdwarf helium groups based on their positions in the U-V velocity diagram, J_z-eccentricity diagram and their Galactic orbits. The overall tendency of the fractional distributions of the four hot subdwarf helium groups in the halo, thin disk and thick disk is largely consistent with the findings reported by Luo et al.(2019) based on LAMOST DR5, which appears to support the predictions of binary population synthesis (Han et al. 2003; 2008). He-deficient stars with -2.2<log(y)<-1 likely origin from stable the Roche lobe overflow channel, He-deficient stars with log(y)<-2.2 from the common envelope ejection channel, and He-rich stars with log(y)>0 from the merger channel of double He white dwarf stars. The fraction of He-rich hot subdwarf stars with -1<log(y)<0 in the thin disk and the halo are far higher than in the thick disk, which implies that these stars have different formation channels in the thin disk and in the halo.

Existing models of galaxy formation have not yet explained striking correlations between structure and star-formation activity in galaxies, notably the sloped and moving boundaries that divide star-forming from quenched galaxies in key structural diagrams. This paper uses these and other relations to ``reverse-engineer'' the quenching process for central galaxies. The basic idea is that star-forming galaxies with larger radii (at a given stellar mass) have lower black-hole masses due to lower central densities. Galaxies cross into the green valley when the cumulative effective energy radiated by their black hole equals $\sim4\times$ their halo-gas binding energy. Since larger-radii galaxies have smaller black holes, one finds they must evolve to higher stellar masses in order to meet this halo-energy criterion, which explains the sloping boundaries. A possible cause of radii differences among star-forming galaxies is halo concentration. The evolutionary tracks of star-forming galaxies are nearly parallel to the green-valley boundaries, and it is mainly the sideways motions of these boundaries with cosmic time that cause galaxies to quench. BH-scaling laws for star-forming, quenched, and green-valley galaxies are different, and most BH mass growth takes place in the green valley. Implications include: the radii of star-forming galaxies are an important second parameter in shaping their black holes; black holes are connected to their halos but in different ways for star-forming, quenched, and green-valley galaxies; and the same BH-halo quenching mechanism has been in place since $z \sim 3$. We conclude with a discussion of black hole-galaxy co-evolution, the origin and interpretation of BH scaling laws.

Using the galaxy catalog built from ELUCID N-body simulation and the semi-analytical galaxy formation model, we have built a mock HI intensity mapping map. We have implemented the Finger-of-God (FoG) effect in the map by considering the galaxy HI gas velocity dispersion. By comparing the HI power spectrum in the redshift space with the measurement from IllustrisTNG simulation, we have found that such FoG effect can explain the discrepancy between current mock map built from N-body simulation and Illustris TNG simulation. Then we built a parameter-free FoG model and a shot-noise model to calculate the HI power spectrum. We found that our model can accurately fit both the monopole and quadrupole moments of the HI matter power spectrum. Our method of building the mock HI intensity map and the parameter-free FoG model will be widely useful for the up-coming 21cm intensity mapping experiments, such as CHIME, Tianlai, BINGO, FAST and SKA. It is also crucial for us to study the non-linear effects in 21cm intensity mapping.

Quenched massive spiral galaxies have attracted great attention recently, as more data is available to constrain their environment and cold gas content. However, the quenching mechanism is still uncertain, as it depends on the mass range and baryon budget of the galaxy. In this letter, we report the identification of a rare population of very massive, quenched spiral galaxies with stellar mass $\gtrsim10^{11}{\rm~M_\odot}$ and halo mass $\gtrsim10^{13}{\rm~M_\odot}$ from the Sloan Digital Sky Survey at redshift $z\sim0.1$. Our CO observations using the IRAM-30m telescope show that these galaxies contain only a small amount of molecular gas. Similar galaxies are also seen in the state-of-the-art semi-analytical models and hydro-dynamical simulations. It is found from these theoretical models that these quenched spiral galaxies harbor massive black holes, suggesting that feedback from the central black holes has quenched these spiral galaxies. This quenching mechanism seems to challenge the popular scenario of the co-evolution between massive black holes and massive bulges.

We construct a turbulent model of the Crab Nebula's non-thermal emission. The present model resolves a number of long-standing problems of the Kennel-Coroniti (1984) model: (i) the sigma problem; (ii) the hard spectrum of radio electrons; (iii) the high peak energy of gamma-ray flares; (iv) and the spacial evolution of the infrared (IR) emission. The Nebula contains two populations of injected particles: Component-I accelerated at the wind termination shock via Fermi-I mechanism, and Component-II accelerated in reconnecting turbulence in highly magnetized ($\sigma$ $\gg 1$) plasma in the central part of the Crab Nebula. The reconnecting turbulence Component-II extends from radio to gamma rays: it accelerate radio electrons with a hard spectrum, destroy the large scale magnetic flux (and thus resolves the sigma-problem), and occasionally produces gamma-ray flares (from the largest scale reconnection events). The model reproduces the broad-band spectrum of the Crab Nebula, from low-frequency synchrotron emission in radio to inverse-Compton emission at TeV energies, as well as spatially resolved evolution of the spectral indices in IR and optical bands.

We study the optical light curves - primarily probing the variable emission from the accretion disk - of ~ 900 extreme variability quasars (EVQs, with maximum flux variations more than 1 mag) over an observed-frame baseline of ~ 16 years using public data from the SDSS Stripe 82, PanSTARRS-1 and the Dark Energy Survey. We classify the multi-year long-term light curves of EVQs into three categories roughly in the order of decreasing smoothness: monotonic decreasing or increasing (3.7%), single broad peak and dip (56.8%), and more complex patterns (39.5%). The rareness of monotonic cases suggests that the major mechanisms driving the extreme optical variability do not operate over timescales much longer than a few years. Simulated light curves with a damped random walk model generally under-predict the first two categories with smoother long-term trends. Despite the different long-term behaviors of these EVQs, there is little dependence of the long-term trend on the physical properties of quasars, such as their luminosity, BH mass, and Eddington ratio. The large dynamic range of optical flux variability over multi-year timescales of these EVQs allows us to explore the ensemble correlation between the short-term (< 6 months) variability and the seasonal-average flux across the decade-long baseline (the rms-mean flux relation). We find that unlike the results for X-ray variability studies, the linear short-term flux variations do not scale with the seasonal-average flux, indicating different mechanisms that drive the short-term flickering and long-term extreme variability of accretion disk emission. Finally, we present a sample of 16 EVQs, where the approximately bell-shaped large amplitude variation in the light curve can be reasonably well fit by a simple microlensing model.

We explore screening effects arising from a relativistic magnetized plasma with applications to Big Bang nucleosynthesis (BBN). %Specifically, due to their small magnetic moments, energies of electrons and positrons can be easily quantized via Landau quantization. The screening potential which depends on the thermodynamics of charged particles in the plasma is altered by the magnetic field. We focus on the impact of screening on the electron capture interaction. Taking into account the correction in BBN arising from a homogeneous primordial magnetic field (PMF), we constrain the epoch at which the PMF was generated and its strength during BBN. Considering such screening corrections to the electron capture rates and using up-to-date observations of primordial elemental abundances, we also discuss the possibility of solving the problem of under-estimation of the deuterium abundance. We find for certain values of the PMF strength predicted D and $^4$He abundances are both consistent with the observational constraints.

This paper studies pseudo-bulges (P-bulges) and classical bulges (C-bulges) in Sloan Digital Sky Survey central galaxies using the new bulge indicator $\Delta\Sigma_1$, which measures relative central stellar-mass surface density within 1 kpc. We compare $\Delta\Sigma_1$ to the established bulge-type indicator $\Delta\langle\mu_e\rangle$ from Gadotti (2009) and show that classifying by $\Delta\Sigma_1$ agrees well with $\Delta\langle\mu_e\rangle$. $\Delta\Sigma_1$ requires no bulge-disk decomposition and can be measured on SDSS images out to $z = 0.07$. Bulge types using it are mapped onto twenty different structural and stellar-population properties for 12,000 SDSS central galaxies with masses 10.0 < log $M_*$/$M_{\odot}$ < 10.4. New trends emerge from this large sample. Structural parameters show fairly linear log-log relations vs. $\Delta\Sigma_1$ and $\Delta\langle\mu_e\rangle$ with only moderate scatter, while stellar-population parameters show a highly non-linear "elbow" in which specific star-formation rate remains roughly flat with increasing central density and then falls rapidly at the elbow, where galaxies begin to quench. P-bulges occupy the low-density end of the horizontal arm of the elbow and are universally star-forming, while C-bulges occupy the elbow and the vertical branch and exhibit a wide range of star-formation rates at fixed density. The non-linear relation between central density and star-formation rate has been seen before, but this mapping onto bulge class is new. The wide range of star-formation rates in C-bulges helps to explain why bulge classifications using different parameters have sometimes disagreed in the past. The elbow-shaped relation between density and stellar indices suggests that central structure and stellar-populations evolve at different rates as galaxies begin to quench.

Observations of high-redshift quasars imply the presence of supermassive black holes already at z~ 7.5. An appealing and promising pathway to their formation is the direct collapse scenario of a primordial gas in atomic-cooling haloes at z ~ 10 - 20, when the $\rm H_2$ formation is inhibited by a strong background radiation field, whose intensity exceeds a critical value, $J_{\rm crit}$. To estimate $J_{\rm crit}$, typically, studies have assumed idealized spectra, with a fixed ratio of $\rm H_{2}$ photo-dissociation rate $k_{\rm H_2}$ to the $\rm H^-$ photo-detachment rate $k_{\rm H^-}$. This assumption, however, could be too narrow in scope as the nature of the background radiation field is not known precisely. In this work we argue that the critical condition for suppressing the $\rm H_2$ cooling in the collapsing gas could be described in a more general way by a combination of $k_{\rm H_2}$ and $k_{\rm H^-}$ parameters. By performing a series of cosmological zoom-in simulations with an encompassing set of $k_{\rm H_2}$ and $k_{\rm H^-}$, we examine the gas flow by following evolution of basic parameters of the accretion flow. We test under what conditions the gas evolution is dominated by $\rm H_{2}$ and/or atomic cooling. We confirm the existence of a critical curve in the $k_{\rm H_2}-k_{\rm H^-}$ plane, and provide an analytical fit to it. This curve depends on the conditions in the direct collapse, and reveals domains where the atomic cooling dominates over the molecular cooling. Furthermore, we have considered the effect of $\rm H_{2}$ self-shielding on the critical curve, by adopting three methods for the effective column density approximation in $\rm H_{2}$. We find that the estimate of the characteristic length-scale for shielding can be improved by using $\lambda_{\rm Jeans25}$, which is 0.25 times that of the local Jeans length.

Aims. Due to the limitations of current computational technology, the fragmentation and isomerization products of vibrationally-excited polycyclic aromatic hydrocarbon (PAH) molecules and their derivatives are poorly studied. In this work, we investigate the intermediate products of PAHs and their derivatives as well as the gas-phase reactions relevant to the interstellar medium, with coronene as a case study. Methods. Based on the semi-empirical method of PM3 as implemented in the CP2K program, molecular dynamics simulations are performed to model the major processes (e.g., vibrations, fragmentations, and isomerizations) of coronene and its derivatives (e.g., methylated coronene, hydrogenated coronene, dehydrogenated coronene, nitrogen-substituted coronene, and oxygen-substituted coronene) at temperatures of 3000 K and 4000 K. Results. We find that the anharmonic effects are crucial for the simulation of vibrational excitation. For the molecules studied here, H2, CO, HCN, and CH2 are the major fragments. Following the dissociation of these small units, most of the molecules could maintain their ring structures, but a few molecules would break completely into carbon chains. The transformation from hexagon to pentagon or heptagon may occur and the heteroatomic substitutions (e.g., N- or O-substitutions) facilitate the transformation.

We use bulge-type classifications of 809 representative SDSS galaxies by Gadotti (2009) to classify a large sample of galaxies into real bulges (classical or elliptical) and pseudobulges using Random Forest. We use structural and stellar population predictors that can easily be measured without image decomposition. Multiple parameters such as the central mass density with 1 kpc, concentration index, S\'{e}rsic index and velocity dispersion do result in accurate bulge classifications when combined together. We classify $\sim 44,500$ face-on galaxies above stellar mass of 10$^{10}$ M$_\odot$ and redshift $ 0.02 < z < 0.07$ into real bulges or pseudobulges with $93 \pm 2$\% accuracy. We show that $\sim 75 - 90\%$ of AGNs identified by the optical line ratio diagnostic are hosted by real bulges. The pseudobulge fraction significantly decreases with AGN signature as the line ratios change from indicating pure star formation ($\sim 54 \pm 4$ \%), to composite of star formation and AGN ($\sim 18 \pm 3$\%), and to AGN-dominated galaxies ($\sim 5 \pm 3$\%). Using the dust-corrected [\ion{O}{3}] luminosity as an AGN accretion indicator, and the stellar mass and radius as proxies for a black hole mass, we find that AGNs in real bulges have lower Eddington ratios than AGNs in pseudobulges. Real bulges have a wide range of AGN and star formation activities, although most of them are weak AGNs. For both bulge-types, their Eddington ratios are correlated with specific star formation rates (SSFR). Real bulges have lower specific accretion rate but higher AGN fraction than pseudobulges do at similar SSFRs.

The splashback radius ($R_{\rm sp}$) of dark matter halos has recently been detected using weak gravitational lensing and cross-correlations with galaxies. However, different methods have been used to measure $R_{\rm sp}$ and to assess the significance of its detection. In this paper, we use simulations to study the precision and accuracy to which we can detect the splashback radius with 3D density, 3D subhalo, and weak lensing profiles. We study how well various methods and tracers recover $R_{\rm sp}$ by comparing it with the value measured directly from particle dynamics. We show that estimates of $R_{\rm sp}$ from density and subhalo profiles correspond to different percentiles of the underlying $R_{\rm sp}$ distribution of particle orbits. At low accretion rates, a second caustic appears and can bias results. Finally, we show that upcoming lensing surveys may be able to constrain the splashback-accretion rate relation directly.

This paper is one in a series that explores the importance of radius as a second parameter in galaxy evolution. The topic investigated here is the relationship between star formation rate (SFR) and galaxy radius ($R_{\rm e}$) for main-sequence star-forming galaxies. The key observational result is that, over a wide range of stellar mass and redshift in both CANDELS and SDSS, there is little trend between SFR and $R_{\rm e}$ at fixed stellar mass. The Kennicutt-Schmidt law, or any similar density-related star formation law, then implies that smaller galaxies must have lower gas fractions than larger galaxies (at fixed $M_{\ast}$), and this is supported by observations of local star-forming galaxies. We investigate the implication by adopting the equilibrium "bathtub" model: the ISM gas mass is assumed to be constant over time and the net star formation rate is the difference between the accretion rate of gas onto the galaxy from the halo and the outflow rate due to winds. To match the observed null correlation between SFR and radius, the bathtub model requires that smaller galaxies at fixed mass have weaker galactic winds. Our hypothesis is that galaxies are a 2-dimensional family whose properties are set mainly by halo mass and concentration. Galaxy radius and accretion rate plausibly both depend on halo concentration, which predicts how wind strength should vary with $R_{\rm e}$ and SFR.

Context. The 3.3 $\mu$m aromatic C-H stretching band of polycyclic aromatic hydrocarbon (PAH) molecules seen in a wide variety of astrophysical regions is often accompanied by a series of weak satellite bands at ~3.4-3.6 $\mu$m. One of these sources, IRAS 21282+5050, a planetary nebula, also exhibits a weak band at ~1.68 $\mu$m. While the satellite features at ~3.4-3.6 $\mu$m are often attributed to the anharmonicities of PAHs, it is not clear whether overtones or combination bands dominate the 1.68 $\mu$m feature. Aims. In this work, we examine the anharmonic spectra of eight PAH molecules, including anthracene, tetracene, pentacene, phenanthrene, chrysene, benz[a]anthracene, pyrene, and perylene, to explore the origin of the infrared bands in the 1.6-1.7 $\mu$m waveelngth region. Methods. Density Functional Theory (DFT) in combination with the vibrational second-order perturbation theory (VPT2) is utilized for computing the anharmonic spectra of PAHs. To simulate the vibrational excitation process of PAHs, the Wang-Landau random walk technique is employed. Results. All the dominant bands in the 1.6-1.7 $\mu$m wavelength range and in the 3.1-3.5 $\mu$m C-H stretching region are calculated and tabulated. It is demonstrated that combination bands dominate the 1.6-1.7 $\mu$m region, while overtones are rare and weak in this region. We also calculate the intensity ratios of the 3.1-3.5 $\mu$m C-H stretching features to the bands in the 1.6-1.7 $\mu$m region, $I_{3.1-3.5}/I_{1.6-1.7}$, for both ground and vibrationally excited states. On average, we obtain $\langle I_{3.1-3.5}/I_{1.6-1.7} \rangle$ $\approx$ 12.6 and $\langle I_{3.1-3.5}/I_{1.6-1.7} \rangle$ $\approx$ 17.6 for PAHs at ground states and at vibrationally excited states, respectively.

The abundance of primordial lithium is derived from the observed spectroscopy of metal-poor stars in the galactic halo. However, the observationally inferred abundance remains at about a factor of three below the abundance predicted by standard big bang nucleosynthesis (BBN). The resolution of this dilemma can be either astrophysical (stars destroy lithium after BBN), nuclear (reactions destroy lithium during BBN), or cosmological, i.e. new physics beyond the standard BBN is responsible for destroying lithium. Here, we overview a variety of possible cosmological solutions, and their shortcomings. On the one hand, we examine the possibility of physical processes that modify the velocity distribution of particles from the usually assumed Maxwell-Boltzmann statistics. A physical justification for this is an inhomogeneous spatial distribution of domains of primordial magnetic field strength as a means to reduce the primordial lithium abundance. Another possibility is that scattering with the mildly relativistic electrons in the background plasma alters the baryon distribution to one resembling a Fermi-Dirac distribution. We show that neither of these possibilities can adequately resolve the lithium problem. A number of alternate hybrid models are discussed including a mix of neutrino degeneracy, unified dark matter, axion cooling, and the presence of decaying and/or charged supersymmetric particles.

Combing Gaia DR2 with LAMOST DR5, we spectroscopically identified 924 hot subdwarf stars, among which 32 stars exhibit strong double-lined composite spectra. We measured the effective temperature $T_{\rm eff}$, surface gravity $\log\,g$, helium abundance $y=n{\rm He}/n{\rm H}$, and radial velocities of 892 non-composite spectra hot subdwarf stars by fitting LAMOST observations with Tlusty/Synspec non-LTE synthetic spectra. We outlined four different groups in the $T_{\rm eff}-\log\,g$ diagram with our helium abundance classification scheme and two nearly parallel sequences in the $T_{\rm eff}-\log(y)$ diagram. 3D Galactic space motions and orbits of 747 hot subdwarf stars with $(G_{BP}-G_{RP})_{0}<-0.36$ mag were computed using LAMOST radial velocities and Gaia parallaxes and proper motions. Based on the $U-V$ velocity diagram, $J_{z}-$eccentricity diagram, and Galactic orbits, we derived Galactic population classifications and the fractional distributions of the four hot subdwarf helium groups in the halo, thin disk and thick disk. Comparisons with the predictions of binary population synthesis calculations (Han 2008) suggest that He-rich hot subdwarf stars with $\log(y)\ge0$ are from the double helium white dwarfs merger, He-deficient hot subdwarf stars with $-2.2\le\log(y)<-1$ from the common envelope ejection, and He-deficient hot subdwarf stars with $\log(y)<-2.2$ from the stable Roche lobe overflow channels. The relative number of He-rich hot subdwarf stars with $-1\le\log(y)<0$ and $\log(y)\ge0$ in the halo is more than twice the prediction of Zhang et al.(2017), even more than six times in the thin disk, which implies that the mergers of helium white dwarfs with low mass main sequence stars may not be the main formation channel of He-rich hot subdwarf stars with $-1\le\log(y)<0$, specially in younger environments.

Observational work have shown that the two brightest satellite galaxies of the Milky Way (MW), the Large Magellanic Cloud (LMC) and the Small Magellanic Cloud (SMC), are rare amongst MW analogues. It is then interesting to know whether the presence of massive satellite has any effect on the whole satellite population in MW analogues. In this article, we investigate this problem using a semi-analytical model combined with the Millennium-II Simulation. MW-analogous galaxies are defined to have similar stellar mass or dark matter halo mass to the MW. We find that, in the first case, the halo mass is larger and there are, on average, twice as many satellites in Milky Way analogs if there is a massive satellite galaxy in the system. This is mainly from the halo formation bias. The difference is smaller if MW analogues are selected using halo mass. We also find that the satellites distribution is slightly asymmetric, being more concentrated on the line connecting the central galaxy and the massive satellite and that, on average, LMC have brought in 14.7 satellite galaxies with $M_{r}<0$ at its accretion, among which 4.5 satellites are still within a distance of 50kpc from the LMC. Considering other satellites, we predict that thereare 7.8 satellites with 50kpc of the LMC. By comparing our model with the early data of Satellites Around Galactic Analogs (SAGA), a survey to observe satellite galaxies around 100 Milky Way analogues, we find that SAGA has more bright satellites and less faint satellites than our model predictions. A future comparison with the final SAGA data is needed.

The advent of new deep+wide photometric lensing surveys will open up the possibility of direct measurements of the dark matter halos of dwarf galaxies. The HSC wide survey will be the first with the statistical capability of measuring the lensing signal with high signal-to-noise at log(M*)=8. At this same mass scale, LSST will have the most overall constraining power with a predicted signal-to-noise for the galaxy-galaxy lensing signal around dwarfs of SN=200. WFIRST and LSST will have the greatest potential to push below the log(M*) = 7 mass scale thanks to the depth of their imaging data. Studies of the dark matter halos of dwarf galaxies at z=0.1 with gravitational lensing are soon within reach. However, further work will be required to develop optimized strategies for extracting dwarfs samples from these surveys, determining redshifts, and accurately measuring lensing on small radial scales. Dwarf lensing will be a new and powerful tool to constrain the halo masses and inner density slopes of dwarf galaxies and to distinguish between baryonic feedback and modified dark matter scenarios.

The design of low-thrust-based multitarget interplanetary missions requires a method to quickly and accurately evaluate the low-thrust transfer between any two visiting targets. Complete evaluation of the low-thrust transfer includes not only the estimation of the optimal fuel consumption but also the judgment of transfer feasibility. In this paper, a deep neural network (DNN)-based method is proposed for quickly evaluating low-thrust transfer. An efficient database generation method is developed for obtaining both the infeasible and optimal transfers. A classification DNN and a regression DNN are trained based on the infeasible and optimal transfers to judge the transfer feasibility and estimate the optimal fuel consumption, respectively. The simulation results show that the well-trained DNNs are capable of quickly determining the transfer feasibility with a correct rate of greater than 98% and approximating the optimal transfer fuel consumption with a relative estimation error of less than 0.4%. The tests on two asteroid chains further show the superiority of the DNN-based method for application to the design of low-thrust-based multitarget interplanetary missions

We investigate the effect on the Big Bang Nucleosynthesis (BBN) from the presence of a stochastic primordial magnetic field (PMF) whose strength is spatially inhomogeneous. We assume a uniform total energy density and a gaussian distribution of field strength. In this case, domains of different temperatures exist in the BBN epoch due to variations in the local PMF. We show that in such case, the effective distribution function of particle velocities averaged over domains of different temperatures deviates from the Maxwell-Boltzmann distribution. This deviation is related to the scale invariant strength of the PMF energy density $\rho_{\rm Bc}$ and the fluctuation parameter $\sigma_{\rm B}$. We perform BBN network calculations taking into account the PMF strength distribution, and deduce the element abundances as functions of the baryon-to-photon ratio $\eta$, $\rho_{\rm Bc}$, and $\sigma_{\rm B}$. We find that the fluctuations of the PMF reduces the $^7$Be production and enhances D production. We analyze the averaged thermonuclear reaction rates compared with those of a single temperature, and find that the averaged charged-particle reaction rates are very different. Finally, we constrain the parameters $\rho_{\rm Bc}$ and $\sigma_{\rm B}$ from observed abundances of $^4$He and D, and find that the $^7$Li abundance is significantly reduced. We also find that if the $\eta$ value during BBN was larger than the present-day value due to a dissipation of the PMF or a radiative decay of exotic particles after BBN or if the stellar depletion of $^7$Li occurred, abundances of all light elements can be consistent with observational constraints.

We present long-term, multi-color photometric monitoring and spectroscopic observations of the short-period M-type eclipsing binary BX Tri. Six flare events were recorded in 4 years from 2014 to 2017. Three of them were detected on one night within an orbital cycle. The strongest one was identified on December 23, 2014. With the amplitudes of $\Delta$B=0.48 mag, $\Delta$V=0.28 mag, $\Delta$R=0.10 mag, $\Delta$I=0.02 mag, the total energy due to the flare event was measured to be 4.08 ($\pm$0.24) $\times$$10^{34}$ erg, exceeding the superflare level ($10^{34}$). Based on the observations, the evolutionary status of the binary system as well as the long-term orbital period changes were analyzed. It reveals that BX Tri is probably a semi-detached system with the primary component filling its Roche lobe. The extremely high occurrence of flare events of the binary could be related to the rapid mass transfer between components.

The distribution function of the relative velocity in a two-body reaction of nonrelativistic uncorrelated particles is derived for general cases of given distribution functions of single particle velocities. The distribution function is then used in calculations of thermonuclear reaction rates. As an example, we take the Tsallis non-Maxwellian distribution, and show that the distribution function of the relative velocity is different from the Tsallis distribution. We identify an inconsistency in previous studies of nuclear reaction rates within Tsallis statistics, and derive revised nuclear reaction rates. Utilizing the revised rates, accurate results of big bang nucleosynthesis are obtained for the Tsallis statistics. For this application it is more difficult to reduce the primordial 7Li abundance while keeping other nuclear abundances within the observational constraints. A small deviation from a Maxwell-Boltzmann distribution can increase the D abundance and slightly reduce 7Li abundance. Although it is impossible to realize a 7Li abundance at the level observed in metal-poor stars, a significant decrease is possible while maintaining a consistency with the observed D abundance.

We have modeled direct collapse of a primordial gas within dark matter halos in the presence of radiative transfer, in high-resolution zoom-in simulations in a cosmological framework, down to the formation of the photosphere and the central object. Radiative transfer has been implemented in the flux-limited diffusion (FLD) approximation. Adiabatic models were run for comparison. We find that (a) the FLD flow forms an irregular central structure and does not exhibit fragmentation, contrary to adiabatic flow which forms a thick disk, driving a pair of spiral shocks, subject to Kelvin-Helmholtz shear instability forming fragments; (b) the growing central core in the FLD flow quickly reaches ~10 Mo and a highly variable luminosity of 10^{38}-10^{39} erg/s, comparable to the Eddington luminosity. It experiences massive recurrent outflows driven by radiation force and thermal pressure gradients, which mix with the accretion flow and transfer the angular momentum outwards; and (c) the interplay between these processes and a massive accretion, results in photosphere at ~10 AU. We conclude that in the FLD model (1) the central object exhibits dynamically insignificant rotation and slower than adiabatic temperature rise with density; (2) does not experience fragmentation leading to star formation, thus promoting the fast track formation of a supermassive black hole (SMBH) seed; (3) inclusion of radiation force leads to outflows, resulting in the mass accumulation within the central 10^{-3} pc, which is ~100 times larger than characteristic scale of star formation. The inclusion of radiative transfer reveals complex early stages of formation and growth of the central structure in the direct collapse scenario of SMBH seed formation.

We performed calculations of the late radio and X-ray afterglow of GRB/GW170817 in the cocoon-jet paradigm, predicting the appearance of a second peak in the afterglow light curve ~ one-three years after the explosion. The model assumes that the prompt emission and early afterglows originate from a cocoon generated during break-out of the delayed magnetically powered jet. As the jet breaks out from the torus-generated wind, a nearly isotropic mildly relativistic outflow is generated; at the same time the primary jet accelerates to high Lorentz factors and avoids detection. As the fast jet slows down, it should become visible to the off-axis observer. Thus, the model has a clear prediction: the X-ray and radio afterglows should first experience a decay, as the cocoon slows down, followed by a rebrightening when the primary jet starts emitting toward an observer.

The intrinsic alignment of galaxies is an important systematic effect in weak-lensing surveys, which can affect the derived cosmological parameters. One direct way to distinguish different alignment models and quantify their effects on the measurement is to produce mocked weak-lensing surveys. In this work, we use full-sky ray-tracing technique to produce mock images of galaxies from the ELUCID $N$-body simulation run with the WMAP9 cosmology. In our model we assume that the shape of central elliptical galaxy follows that of the dark matter halo, and spiral galaxy follows the halo spin. Using the mocked galaxy images, a combination of galaxy intrinsic shape and the gravitational shear, we compare the predicted tomographic shear correlations to the results of KiDS and DLS. It is found that our predictions stay between the KiDS and DLS results. We rule out a model in which the satellite galaxies are radially aligned with the center galaxy, otherwise the shear-correlations on small scales are too high. Most important, we find that although the intrinsic alignment of spiral galaxies is very weak, they induce a positive correlation between the gravitational shear signal and the intrinsic galaxy orientation (GI). This is because the spiral galaxy is tangentially aligned with the nearby large-scale overdensity, contrary to the radial alignment of elliptical galaxy. Our results explain the origin of detected positive GI term from the weak-lensing surveys. We conclude that in future analysis, the GI model must include the dependence on galaxy types in more detail.

The alignment between satellites and central galaxies has been studied in detail both in observational and theoretical works. The widely accepted fact is that the satellites preferentially reside along the major axis of their central galaxy. However, the origin and large-scale environment dependence of this alignment are still unknown. In an attempt to figure out those, we use data constructed from SDSS DR7 to investigate the large-scale environmental dependence of this alignment with emphasis on examining the alignments' dependence on the colour of the central galaxy. We find a very strong large-scale environmental dependence of the satellite-central alignment in groups with blue centrals. Satellites of blue centrals in knots are preferentially located perpendicular to the major axis of the centrals, and the alignment angle decreases with environment namely when going from knots to voids. The alignment angle strongly depend on the ${}^{0.1}(g-r)$ colour of centrals. We suggest that the satellite-central alignment is the result of a competition between satellite accretion within large scale-structure and galaxy evolution inside host haloes. For groups containing red central galaxies, the satellite-central alignment is mainly determined by the evolution effect, while for blue central dominated groups, the effect of large-scale structure plays a more important role, especially in knots. Our results provide an explanation for how the satellite-central alignment forms within different large-scale environments. The perpendicular case in groups and knots with blue centrals may also provide insight into understanding similar polar arrangements such the formation of the Milky Way and Centaurus A's satellite system.

Direct collapse within dark matter (DM) halos is a promising path to form supermassive black hole (SMBH) seeds at high redshifts. The outer part of this collapse remains optically thin, and has been studied intensively using numerical simulations. However, the innermost region of the collapse is expected to become optically thick and requires us to follow the radiation field in order to understand its subsequent evolution. So far, the adiabatic approximation has been used exclusively for this purpose. We apply radiative transfer in the flux-limited diffusion (FLD) approximation to solve the evolution of coupled gas and radiation, for isolated halos. For direct collapse within isolated DM halos, we find that (1) the photosphere forms at ~10^{-6} pc and rapidly expands outward. (2) A central core forms, with a mass of ~1 Mo, supported by thermal gas pressure gradients and rotation. (3) Growing thermal gas and radiation pressure gradients dissolve it. (4) This process is associated with a strong anisotropic outflow, and another core forms nearby and grows rapidly. (5) Typical radiation luminosity emerging from the photosphere encompassing these cores is ~5 x 10^{37}-5 x 10^{38} erg/s, of order the Eddington luminosity. (6) Two variability timescales are associated with this process: a long one, which is related to the accretion flow within the central ~10^{-4}-10^{-3} pc, and ~0.1 yr, which is related to radiation diffusion. (7) Adiabatic models have been run for comparison and their evolution differs profoundly from that of the FLD models, by forming a central geometrically-thick disk. Overall, an adiabatic equation of state is not a good approximation to the advanced stage of direct collapse, mainly because the radiation in the FLD is capable of escaping due to anisotropy in the optical depth and associated gradients.

Stellar fundamental parameters are important in the asteroseismic study of Kepler light curves. However, the most used estimates in the Kepler Input Catalog (KIC) are not accurate enough for hot stars. Using a sample of B stars from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) spectral survey, we confirmed the systematic underestimation in KIC effective temperature and overestimation in KIC surface gravity. The good agreement between LAMOST and other follow-up observations proved the accuracy of effective temperature and surface gravity of B stars derived from LAMOST low-resolution spectra. By searching through LAMOST data, we found four misclassified main-sequence B stars in the Kepler field, which had been previously classified as A-type variables. We present spectroscopic and detailed frequency analysis of these four stars based on LAMOST spectra and Kepler photometry.

The detection of highly ionized metal absorption lines in the X-ray spectra of the Galactic X-ray binaries (XRBs) implies the distribution of hot gas along the sightline toward the background sources. However, the origin of this hot gas is still unclear: it can arise in the hot interstellar medium (ISM), or is intrinsic to the XRBs. In this paper, we present an XMM-Newton survey of the OVII absorption lines in the spectra of Galactic XRBs. A total of 33 XRBs were selected, with 29 low mass XRBs and 4 high mass XRBs. At more than 3$\sigma$ threshold, OVII absorption line was detected in 16 targets, among which 4 were newly discovered in this work. The average line equivalent width (EW) is centered around $\sim$ 20 m\AA. Additionally we do not find strong correlations between the OVII EWs and the Galactic neutral absorption N$_{\rm H}$, the Galactic coordinates or the distance of background targets. Such non-correlation may suggest the contamination of the circumstellar material, or the lack of constrains on the line Doppler-b parameter. We also find regardless the direction of the XRBs, the OVII absorption lines were always detected when the flux of the background XRBs reach a certain level, suggesting a uniform distribution of this hot gas. We estimate a ratio of 0.004 --- 0.4 between hot and neutral phase of the interstellar medium. This is the second paper in the series following Fang et al. 2015, in which we focused on the local OVII absorption lines detected in the background AGN spectra. Detailed modeling of the hot ISM distribution will be investigated in a future paper.

In this paper we use high-resolution cosmological simulations to study halo intrinsic alignment and its dependence on mass, formation time and large-scale environment. In agreement with previous studies using N-body simulations, it is found that massive halos have stronger alignment. For given mass, older halos have stronger alignment than younger ones. By identifying the cosmic environment of halo using Hessian matrix, we find that for given mass, halos in cluster regions also have stronger alignment than those in filament. The existing theory has not addressed these dependencies explicitly. In this work we extend the linear alignment model with inclusion of halo bias and find that the halo alignment with its mass and formation time dependence can be explained by halo bias. However, the model can not account for the environment dependence, as it is found that halo bias is lower in cluster and higher in filament. Our results suggest that halo bias and environment are independent factors in determining halo alignment. We also study the halo alignment correlation function and find that halos are strongly clustered along their major axes and less clustered along the minor axes. The correlated halo alignment can extend to scale as large as $100h^{-1}$Mpc where its feature is mainly driven by the baryon acoustic oscillation effect.

Time-series, multi-color photometry and high-resolution spectra of the short period eclipsing binary V Tri were obtained by observations. The completely covered light and radial velocity curves of the binary system are presented. All times of light minima derived from both photoelectric and CCD photometry were used to calculate the orbital period and new ephemerides of the eclipsing system. The analysis of $O-C$ diagram reveals that the orbital period is $0.58520481\ days$, decreasing at a rate of $dP/dt=-7.80\times10^{-8} d\ yr^{-1} $. The mass transfer between the two components and the light time-travel effect due to a third body could be used to explain the period decrease. However, a semidetached configuration with the less-mass component filling and the primary nearly filling each of their Roche lobes was derived from the synthesis of the light and radial velocity curves by using the 2015 version of the Wilson-Devinney code. We consider the period decrease to be the nonconservative mass transfer from the secondary component to the primary and the mass loss of the system, which was thought to be an EB type while it should be an EA type (semi-detached Algol-type) from our study. The masses, radii and luminosities of the primary and secondary are $1.60\pm0.07 M_\odot$, $1.64\pm0.02 R_\odot$, $14.14\pm0.73 L_\odot$ and $0.74\pm0.02 M_\odot$, $1.23\pm0.02 R_\odot$, $1.65\pm0.05 L_\odot$, respectively.

Semi-analytic models of galaxy formation are powerful tools to study the evolution of galaxy population in a cosmological context. However, most models over-predict the number of low-mass galaxies at high redshifts and the color of model galaxies are not right in the sense that low-mass satellite galaxies are too red and centrals are too blue. The recent version of the L-Galaxies model by Henriques et al.(H15) is a step forward to solve these problems by reproducing the evolution of stellar mass function and the overall fraction of red galaxies. In this paper we compare the two model predictions of L-Galaxies (the other is Guo et al. , G13) to the SDSS data in detail. We find that in the H15 model the red fraction of central galaxies now agrees with the data due to their implementation of strong AGN feedback, but the stellar mass of centrals in massive haloes is now slightly lower than the data. For satellite galaxies, the red fraction of low-mass galaxies ($\log M_{*}/M_{\odot} < 10$) also agrees with the data, but the color of massive satellites ($10 < \log M_{*}/M_{\odot} < 11$) is slightly bluer. The correct color of centrals and bluer color of massive satellites indicate that the quenching in massive satellites are not strong enough. We also find that there are too much red spirals and less bulge-dominated galaxies in both H15 and G13 models. Our results suggest that additional mechanisms, such as more minor merger or disk instability, are needed to slightly increase the stellar mass of central galaxy in massive galaxies, mainly in the bulge component, and the bulge dominated galaxy will be quenched or then be quenched by any other mechanisms.

The Milky Way presents the too-big-to-fail (TBTF) problem that there are two observed satellite galaxies with maximum circular velocity larger than 55km/s, and others have velocity less than 25km/s, but the cold dark matter model predicts there should be more than 10 subhaloes with velocity larger than 25km/s. Those massive subhaloes with $25km/s<V_{max}<55km/s$ should not have failed to form stars. The TBTF problem severely challenges the CDM model. Most efforts are seeking the effects of baryonic feedback, decreasing the mass of the Milky Way, changing the properties of dark matter, so as to assign the observed low-velocity satellites into the massive subhaloes found in simulations. However, the TBTF problem can be avoided if the MW have not accreted subhaloes with velocity between $25km/s<V_{max}<55km/s$ although the probability of such a gap is lower as $\sim 1\%$ and can not be tested against observations. In this work we study the gap in stellar mass of satellite galaxies using the SDSS group catalogue and a semi-analytical model. We find that there are 1-2\% of galaxy groups with a large gap in the stellar mass of their satellites. These 'big gap' groups have accreted less massive subhaloes in their formation history and naturally display a gap between their satellite galaxies. If extrapolating our results to the Milky Way is appropriate, we conclude that it is very likely that our Milky Way has not accreted enough massive subhaloes to host those low-velocity satellites, and the TBTF problem is naturally avoided.

We provide detailed comparison between the AMR code Enzo-2.4 and the SPH/N- body code GADGET-3 in the context of isolated or cosmological direct baryonic collapse within dark matter (DM) halos to form supermassive black holes. Gas flow is examined by following evolution of basic parameters of accretion flows. Both codes show an overall agreement in the general features of the collapse, however, many subtle differences exist. For isolated models, the codes increase their spatial and mass resolutions at different pace, which leads to substantially earlier collapse in SPH than in AMR cases due to higher gravitational resolution in GADGET-3. In cosmological runs, the AMR develops a slightly higher baryonic resolution than SPH during halo growth via cold accretion permeated by mergers. Still, both codes agree in the buildup of DM and baryonic structures. However, with the onset of collapse, this difference in mass and spatial resolution is amplified, so evolution of SPH models begins to lag behind. Such a delay can have effect on formation/destruction rate of H2 due to UV background, and on basic properties of host halos. Finally, isolated non-cosmological models in spinning halos, with spin parameter {\lambda} ~ 0.01 - 0.07, show delayed collapse for greater {\lambda}, but pace of this increase is faster for AMR. Within our simulation setup, GADGET-3 requires significantly larger computational resources than Enzo- 2.4 during collapse, and needs similar resources, during the pre-collapse, cosmological structure formation phase. Yet it benefits from substantially higher gravitational force and hydrodynamic resolutions, except at the end of collapse.

The quenching of star formation in satellite galaxies is observed over a wide range of dark matter halo masses and galaxy environments. In the recent Guo et al (2011) and Fu et al (2013) semi-analytic + N-body models, the gaseous environment of the satellite galaxy is governed by the properties of the dark matter subhalo in which it resides. This quantity depends of the resolution of the N-body simulation, leading to a divergent fraction of quenched satellites in high- and low-resolution simulations. Here, we incorporate an analytic model to trace the subhaloes below the resolution limit. We demonstrate that we then obtain better converged results between the Millennium I and II simulations, especially for the satellites in the massive haloes ($\rm log M_{halo}=[14,15]$). We also include a new physical model for the ram-pressure stripping of cold gas in satellite galaxies. However, we find very clear discrepancies with observed trends in quenched satellite galaxy fractions as a function of stellar mass at fixed halo mass. At fixed halo mass, the quenched fraction of satellites does not depend on stellar mass in the models, but increases strongly with mass in the data. In addition to the over-prediction of low-mass passive satellites, the models also predict too few quenched central galaxies with low stellar masses, so the problems in reproducing quenched fractions are not purely of environmental origin. Further improvements to the treatment of the gas-physical processes regulating the star formation histories of galaxies are clearly necessary to resolve these problems.

We present a catalog of 166 spectroscopically identified hot subdwarf stars from LAMOST DR1, 44 of which show the characteristics of cool companions in their optical spectra. Atmospheric parameters of 122 non-composite spectra subdwarf stars were measured by fitting the profiles of hydrogen (H) and helium (He) lines with synthetic spectra from non-LTE model atmospheres. Most of the sdB stars scatter near the Extreme Horizontal Branch in the $T_{\rm eff}-\log{g}$ diagram and two well defined groups can be outlined. A clustering of He-enriched sdO stars appears near $T_{\rm eff}=45\,000$ K and $\log(g) = 5.8$. The sdB population separates into several nearly parallel sequences in the $T_{\rm eff}-{\rm He}$ abundance diagram with clumps corresponding to those in the $T_{\rm eff}-\log{g}$ diagram. Over $38\,000$ K (sdO) stars show abundance extremes, they are either He-rich or He-deficient and we observe only a few stars in the $ -1 < \log(y) < 0$ abundance range. With increasing temperature these extremes become less prominent and the He abundance approaches to $\log(y)\sim-0.5$. A unique property of our sample is that it covers a large range in apparent magnitudes and galactic latitudes, therefore it contains a mix of stars from different populations and galactic environments. Our results are consistent with the findings of Hirsch (2009) and we conclude that He-rich and He-deficient sdB stars ($\log(y) < 1$) probably origin from different populations. We also find that most sdO and sdB stars lie in a narrow strip in the luminosity and helium abundance plane, which suggests that these atmospheric parameters are correlated.

We present multi-color CCD photometry of the neglected contact binary XZ Leo. Completely covered VRI band light curves and four times of minimum light were obtained. Combining the photometric and previously published radial velocity data, a revised photometric analysis was carried out for the binary system by applying the Wilson-Devinney code. With a hot spot placed on the massive primary component near the neck region of the common envelope, the light curves were satisfactorily modeled. The photometric solution combined with the radial velocity solution reveals that XZ Leo is an A-type contact binary with a degree of contact of 24($\pm1)\%$. The absolute parameters of the components were determined as M_1 = 1.74($\pm$0.06)M_\odot, M_2 = 0.61($\pm$0.02)M_\odot, R_1 = 1.69($\pm$0.01)R_\odot, R_2 = 1.07($\pm0.01$)R_\odot, L_1 = 6.73($\pm0.08$) L_\odot, L_2 = 2.40($\pm$0.04)L_\odot. Based on all the available data, the long-term orbital period behavior of the system was investigated. It indicates that the binary system was undergoing continuous orbital period increase in the past three decades with a rate of dP/dt = + 6.12 \times {10^{-8}} days yr^{-1}, which suggests a probable mass transfer from the secondary to the primary component at a rate of dM/dt= 3.92\times 10^{-8} M_\odot yr^{-1}. The binary system is expected to evolve into the broken-contact stage in $1.56 \times 10^6$ years. This could be evidence supporting the Thermal Relaxation Oscillation theory.

The Large sky Area Multi-Object Spectroscopic Telescope (LAMOST) General Survey is a spectroscopic survey that will eventually cover approximately half of the celestial sphere and collect 10 million spectra of stars, galaxies and QSOs. Objects both in the pilot survey and the first year general survey are included in the LAMOST First Data Release (DR1). The pilot survey started in October 2011 and ended in June 2012, and the data have been released to the public as the LAMOST Pilot Data Release in August 2012. The general survey started in September 2012, and completed its first year of operation in June 2013. The LAMOST DR1 includes a total of 1202 plates containing 2,955,336 spectra, of which 1,790,879 spectra have observed signal-to-noise S/N >10. All data with S/N>2 are formally released as LAMOST DR1 under the LAMOST data policy. This data release contains a total of 2,204,696 spectra, of which 1,944,329 are stellar spectra, 12,082 are galaxy spectra and 5,017 are quasars. The DR1 includes not only spectra, but also three stellar catalogues with measured parameters: AFGK-type stars with high quality spectra (1,061,918 entries), A-type stars (100,073 entries), and M stars (121,522 entries). This paper introduces the survey design, the observational and instrumental limitations, data reduction and analysis, and some caveats. Description of the FITS structure of spectral files and parameter catalogues is also provided.

We report the results of a search for variable stars in the open cluster NGC 2141. Ten variable stars are detected, among which nine are new variable stars and they are classified as three short period W UMa type eclipsing binaries, two EA type eclipsing binaries, one EB type eclipsing binary, one very short period RS CVn type eclipsing binary, one d type RR Lyrae variable star, and one unknown type variable star. The membership and physical properties are discussed, based on their light curves, positions in the CMDs, spatial locations and periods. A known EB type eclipsing binary is also identified as a blue struggler candidate of the cluster. Furthermore, we find that all eclipsing contact binaries have prominently asymmetric eclipses and O Connell effect (O Connell 1951) which increases with the decrease of the orbital periods. This suggests that the O Connell effect is probably related to the evolution of the orbital period in short period eclipsing binary systems.

X-ray observations of the Galactic X-ray binaries (XRB) revealed numerous highly ionized metal absorption lines. However, it is unclear whether such lines are produced by the hot interstellar medium (ISM) or the circumstellar medium (CSM) intrinsic to the binaries. Here we present a Chandra X-ray absorption line study of 28 observations on 12 X-ray binaries, with a focus on the NeIX and FeXVII lines. We report the first detections of these lines in a significant amount of observations. We do not find significant dependence of the line equivalent width on the distance of the XRBs, but we do see weak dependence on the source X-ray luminosity. We also find two out of twelve selected targets show strong temporal variation of the NeIX absorbers. While the line ratio between the two ion species suggests a temperature consistent with the previous predictions of the ISM, comparing with two theoretical models of the ISM shows the observed column densities are significantly higher than predictions. On the other hand, photoionzation by the XRBs provides reasonably good fit to the data. Our findings suggest that a significant fraction of these X-ray absorbers may originate in the hot gas intrinsic to the X-ray binaries, and that the ISM makes small, if not negligible, contribution. We briefly discuss the implications to the study of the Milky Way hot gas content.

Some experiments have found multi-core events in cosmic rays extensive air showers which should be interpreted by hadronic interaction theory. In this paper, the multi-core events are reproduced by Monte Carlo simulation with CORSIKA. The origin of each sub-cores is tracked back from the observation level. The interaction mechanism and original particles of sub-core are studied in this paper.

The neutrino-cooled accretion disk, which was proposed to work as the central engine of gamma-ray bursts, encounters difficulty in interpreting the X-ray flares after the prompt gamma-ray emission. In this paper, the magnetic coupling between the inner disk and the central black hole is taken into consideration. For mass accretion rates around $0.001 \sim 0.1 M_{\sun} s^{-1}$, our results show that the luminosity of neutrino annihilation can be significantly enhanced due to the coupling effects. As a consequence, after the gamma-ray emission, a remnant disk with mass $M_{disk} \la 0.5 M_{\sun}$ may power most of the observed X-ray flares with the rest frame duration less than 100 seconds. In addition, a comparison between the magnetic coupling process and the Blandford-Znajek mechanism is shown on the extraction of black hole rotational energy.

We present new quasars discovered in the vicinity of the Andromeda and Triangulum galaxies with the LAMOST during the 2010 and 2011 observational seasons. Quasar candidates are selected based on the available SDSS, KPNO 4 m telescope, XSTPS optical, and WISE near infrared photometric data. We present 509 new quasars discovered in a stripe of ~135 sq. deg from M31 to M33 along the Giant Stellar Stream in the 2011 pilot survey datasets, and also 17 new quasars discovered in an area of ~100 sq. deg that covers the central region and the southeastern halo of M31 in the 2010 commissioning datasets. These 526 new quasars have i magnitudes ranging from 15.5 to 20.0, redshifts from 0.1 to 3.2. They represent a significant increase of the number of identified quasars in the vicinity of M31 and M33. There are now 26, 62 and 139 known quasars in this region of the sky with i magnitudes brighter than 17.0, 17.5 and 18.0 respectively, of which 5, 20 and 75 are newly-discovered. These bright quasars provide an invaluable collection with which to probe the kinematics and chemistry of the ISM/IGM in the Local Group of galaxies. A total of 93 quasars are now known with locations within 2.5 deg of M31, of which 73 are newly discovered. Tens of quasars are now known to be located behind the Giant Stellar Stream, and hundreds behind the extended halo and its associated substructures of M31. The much enlarged sample of known quasars in the vicinity of M31 and M33 can potentially be utilized to construct a perfect astrometric reference frame to measure the minute PMs of M31 and M33, along with the PMs of substructures associated with the Local Group of galaxies. Those PMs are some of the most fundamental properties of the Local Group.

The propagation of emissions in the relativistic streaming pair plasmas of pulsars and magnetars is studied within the framework of Gordon effective metric theory, and the effect of vacuum polarization on the radiation spectra is examined particularly. It is found that the electromagnetic and kinetic effects of streaming magnetized plasmas and the vacuum polarization effect can shift the radiation spectra of stars. The modification of redshifts of spectra can reach the same magnitude as the gravitational redshift in particular cases. Moreover, the redshifts induced by the media are anisotropic and associated with wave modes. For the O-mode, the modification of redshifts is dependent on the frequency of the radiation. For lower-frequency radiation, the modification is dominated by the plasma effect, while for higher-frequency radiations, by the vacuum polarization effect. For the X-mode, the modification is non-dispersive and dominated by the vacuum polarization effect. The vacuum polarization also has a significant effect on the redshifts of X-ray emissions in magnetars.

Galaxies observed today are likely to have evolved from density perturbations in the early universe. Perturbations that exceeded some critical threshold are conjectured to have undergone gravitational collapse to form primordial black holes (PBHs) at a range of masses. Such PBHs serve as candidates for cold dark matter and their detection would shed light on conditions in the early universe. Here we propose a mechanism to search for transits of PBHs through/nearby Earth by studying the associated seismic waves. Using a spectral-element method, we simulate and visualize this seismic wave field in Earth's interior. We predict the emergence of two unique signatures, namely, a wave that would arrive almost simultaneously everywhere on Earth's free surface and the excitation of unusual spheroidal modes with a characteristic frequency-spacing in free oscillation spectra. These qualitative characteristics are unaffected by the speed or proximity of the PBH trajectory. The seismic energy deposited by a proximal ${M^{PBH} = 10^{15}}$ g PBH is comparable to a magnitude $M_w=4$ earthquake. The non-seismic collateral damage due to the actual impact of such small PBHs with Earth would be negligible. Unfortunately, the expected collision rate is very low even if PBHs constituted all of dark matter, at ${\sim 10^{-7} {yr}^{-1}}$, and since the rate scales as ${1/M^{PBH}}$, fortunately encounters with larger, Earth-threatening PBHs are exceedingly unlikely. However, the rate at which non-colliding close encounters of PBHs could be detected by seismic activity alone is roughly two orders of magnitude larger --- that is once every hundred thousand years --- than the direct collision rate.

Type Ia supernovae (SNe Ia) play an important role in diverse areas of astrophysics, from the chemical evolution of galaxies to observational cosmology. However, the nature of the progenitors of SNe Ia is still unclear. In this paper, according to a detailed binary population synthesis study, we obtained SN Ia birthrates and delay times from different progenitor models, and compared them with observations. We find that the Galactic SN Ia birthrate from the double-degenerate (DD) model is close to those inferred from observations, while the birthrate from the single-degenerate (SD) model accounts for only about 1/2-2/3 of the observations. If a single starburst is assumed, the distribution of the delay times of SNe Ia from the SD model is a weak bimodality, where the WD + He channel contributes to the SNe Ia with delay times shorter than 100Myr, and the WD + MS and WD + RG channels to those with age longer than 1Gyr.

With the development of network and the World Wide Web (WWW), the Internet has been growing and changing dramatically. More and more on-line database systems and different kinds of services are available for astronomy research. How to help users find their way through the jungle of information services becomes an important challenge. Although astronomers have been aware of the importance of interoperability and introduced the concept of Virtual Observatory as a uniform environment for future astronomical on-line resources and services, transparent access to heterogeneous on-line information is still difficult. SkyMouse is a lightweight interface for distributed astronomical on-line resources and services, which is designed and developed by us, i.e., Chinese Virtual Observatory Project. Taking advantage of screen word-capturing technology, different kinds of information systems can be queried through simple mouse actions, and results are returned in a uniform web page. SkyMouse is an easy to use application, aiming to show basic information or to create a comprehensive overview of a specific astronomical object. In this paper current status of on-line resources and services access is reviewed; system architecture, features and functions of SkyMouse are described; challenges for intelligent interface for on-line astronomical resources and services are discussed.