Using the NOvA neutrino detectors, a broad search has been performed for any signal coincident with 28 gravitational wave events detected by the LIGO/Virgo Collaboration between September 2015 and July 2019. For all of these events, NOvA is sensitive to possible arrival of neutrinos and cosmic rays of GeV and higher energies. For five (seven) events in the NOvA Far (Near) Detector, timely public alerts from the LIGO/Virgo Collaboration allowed recording of MeV-scale events. No signal candidates were found.

The integrated Sachs-Wolfe (ISW) effect is caused by the decay of cosmological gravitational potential, and is therefore a unique probe of dark energy. However, its robust detection is still problematic. Various tensions between different data sets, different large scale structure (LSS) tracers, and between data and the $\Lambda$CDM theory prediction, exist. We propose a novel method of ISW measurement by cross correlating CMB and the LSS traced by "low-density-position" (LDP, \citet{2019ApJ...874....7D}). It isolates the ISW effect generated by low-density regions of the universe, but insensitive to selection effects associated with voids. We apply it to the DR8 galaxy catalogue of the DESI Legacy imaging surveys, and obtain the LDPs at $z\leq 0.6$ over $\sim$ 20000 $deg^2$ sky coverage. We then cross correlate with the Planck temperature map, and detect the ISW effect at $3.2\sigma$. We further compare the measurement with numerical simulations of the concordance $\Lambda$CDM cosmology, and find the ISW amplitude parameter $A_{ISW}=1.14\pm0.38$ when we adopt a LDP definition radius $R_s=3^{'}$, fully consistent with the prediction of the standard $\Lambda$CDM cosmology ($A_{ISW}=1$). This agreement with $\Lambda$CDM cosmology holds for all the galaxy samples and $R_s$ that we have investigated. Furthermore, the S/N is comparable to that of galaxy ISW measurement. These results demonstrate the LDP method as a competitive alternative to existing ISW measurement methods, and provide independent checks to existing tensions.

Non-Gaussianity of the cosmological matter density field can be largely reduced by a local Gaussianization transformation (and its approximations such as the logrithmic transformation). Such behavior can be recasted as the Gaussian copula hypothesis, and has been verified to very high accuracy at two-point level. On the other hand, statistically significant non-Gaussianities in the Gaussianized field have been detected in simulations. We point out that, this apparent inconsistency is caused by the very limited degrees of freedom in the copula function, which make it misleading as a diagnosis of residual non-Gaussianity in the Gaussianized field. Using the copula density, we highlight the departure from Gaussianity. We further quantify its impact in the predicted n-point correlation functions. We explore a remedy of the Gaussian copula hypothesis, which alleviates but not completely solves the above problems.

The joint observation of GW170817 and GRB 170817A indicated that short gamma-ray bursts (SGRBs) can originate from binary neutron star mergers. Moreover, some SGRBs could be detected off-axis, while the SGRB jets are highly structured. Then, by assuming an universal angular distribution of the jet emission for all SGRBs, we re-produce the flux and redshift distributions of the cosmological SGRBs detected by {\it Swift} and {\it Fermi}. For self-consistency, this angular distribution is simultaneously constrained by the luminosity and event rate of GRB 170817A. As a result, it is found that the universal jet structure of SGRBs could approximately have a two-Gaussian profile. Meanwhile, the intrinsic luminosity function (LF) of the on-axis emission of the jets can be simply described by a single power law with a low-luminosity exponential cutoff. The usually discovered broken-power-law apparent LF for relatively high luminosities can naturally result from the coupling of the intrinsic LF with the angular distribution of the jet emission, as the viewing angles to the SGRBs are arbitrarily distributed.

The recent discovery of a fast radio burst (FRB 200428) from the Galactic magnetar SGR 1935+2154 robustly indicated that FRB phenomena can sometimes be produced by magnetars, although it is uncertain whether the cosmological FRBs can share the same origin with this Galactic event. The association of FRB 200428 with an X-ray burst (XRB) further offers important implications for the physical processes responsible for the FRB phenomena. By assuming that the XRB emission is produced in the magnetosphere of the magnetar, we investigate the possibility of that the FRB emission is produced by the synchrotron maser (SM) mechanism, which is powered by a shock due to the collision of an $e^{\pm}$ ejecta with a baryonic cloud. It is found that this shock-powered SM model can in principle account for the FRB 200428 observations, if the collision just occurred on the line of sight and the ejecta lunched by magnetar bursts can have appropriate ingredients and structures. To be specific, a burst ejecta should consist of an ultra-relativistic and extremely highly collimated $e^{\pm}$ component and a sub-relativistic and wide-spreading baryonic component. The cloud blocking the $e^{\pm}$ ejecta is just a remnant of a previous baryonic ejecta. Meanwhile, as a result of the synchrotron emission of the shocked material, an intense millisecond X-ray pulse is predicted to overlap the magnetosphere XRB emission, which in principle provides a way to test the model. Additionally, the peak frequency of the SM radiation is constrained to be about a few hundred MHz and the radiation efficiency is around $10^{-4}$.

In this paper, we present a numerical method to study the predicted lightcurves as a function of viewing angle. We extrapolate the fitting formulae for the mass and velocity of tidal dynamical ejecta across a wide range of mass ratio validated with 66 simulations and use them in the calculations of the kilonova lightcurves. The calculated peak luminosity of a BH-NS merger kilonova is typically about a few times $10^{41}\ {\rm erg\ s^{-1}}$, which is always $\lesssim4.5\times10^{41}\ {\rm erg\ s^{-1}}$. This corresponds to the AB absolute magnitudes fainter than $\sim -15\ {\rm mag}$ in optical and $\sim -16\ {\rm mag}$ in infrared. Since the projected photosphere area of the dynamical ejecta is much larger than that of the disk wind outflows, the dynamical ejecta usually contribute to the majority of the kilonova emission from BH-NS mergers. The fitted blackbody temperature and the shape of the observed multi-band lightcurves are insensitive to the line of sight. The peak time of the observed multi-band lightcurves, affected by the light propagation effect, is related to the relative motion direction between the dynamical ejecta and the observer. The observed luminosity varies with the projected photosphere area determined by the viewing angles. However, the predicted peak luminosity only varies by a factor of $\sim (2 - 3)$ (or by $\sim1\ {\rm mag}$) for different viewing angles. When the short-duration gamma-ray burst afterglow is taken into account, for an on-axis geometry, the kilonova emission is usually outshone by the afterglow emission and can be only observed in the redder bands, especially in the $K$-band at late times. Compared with GW170817/AT2017gfo, the BH-NS merger kilonovae are optically dim but possibly infrared bright. At the same epoch after the merger, the blackbody fitting temperature of the BH-NS merger kilonovae is lower than that of GW170817/AT2017gfo.

It is believed that the relativistic jets of gamma-ray bursts (GRBs) should initially propagate through a heavy envelope of the massive progenitor stars or through a merger ejecta formed from the compact binary mergers. The interaction of a jet with a stellar envelope or a merger ejecta can lead to the deceleration of the head material of the jet and, simultaneously, the formation of a hot cocoon. However, this jet-envelope/ejecta interaction is actually undetectable with electromagnetic radiation and can only be inferred indirectly by the structure of the breakout jet. Therefore, as a way out of this difficulty, it is suggested that the jet-envelope/ejecta interaction can produce a gravitational wave memory of an amplitude of $h\sim10^{-26}-10^{-23}$, which could be detected with some future gravitational wave detectors sensitive in the frequency range from sub-Hertz to several tens of Hertz. This provides a potential direct way to probe the jet propagation and then the interior of the GRB progenitors.

In this work, we propose a powerful probe of neutrino effects on the large-scale structure (LSS) of the Universe, i.e., Minkowski functionals (MFs). The morphology of LSS can be fully described by four MFs. This tool, with strong statistical power, is robust to various systematics and can comprehensively probe all orders of N-point statistics. By using a pair of high-resolution N-body simulations, for the first time, we comprehensively studied the subtle neutrino effects on the morphology of LSS. For an ideal LSS survey of volume $\sim1.73$ Gpc$^3$/$h^3$, neutrino signals are mainly detected from void regions with a significant level up to $\thicksim 10\sigma$ and $\thicksim 300\sigma$ for CDM and total matter density fields, respectively. This demonstrates its enormous potential for much improving the neutrino mass constraint in the data analysis of up-coming ambitious LSS surveys.

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

We report the discovery of a binary millisecond pulsar (namely PSR J1641+3627F or M13F) in the globular cluster M13 (NGC 6205) and timing solutions of M13A to F using observations made with the Five-hundred-metre Aperture Spherical radio Telescope (FAST). PSR J1641+3627F has a spin period of 3.00 ms and an orbital period of 1.4 days. The most likely companion mass is 0.16 M$_{\odot}$. M13A to E all have short spin periods and small period derivatives. We also confirm that the binary millisecond pulsar PSR J1641$+$3627E (also M13E) is a black widow with a companion mass around 0.02 M$_{\odot}$. We find that all the binary systems have low eccentricities compared to those typical for globular cluster pulsars and that they decrease with distance from the cluster core. This is consistent with what is expected as this cluster has a very low encounter rate per binary.

Young neutron stars (NSs) born in core-collapse explosions are promising candidates for the central engines of fast radio bursts (FRBs), since the first localized repeating burst FRB 121102 happens in a star forming dwarf galaxy, which is similar to the host galaxies of superluminous supernovae (SLSNe) and long gamma-ray bursts (LGRBs). However, FRB 180924 and FRB 190523 are localized to massive galaxies with low rates of star formation, compared with the host of FRB 121102. Meanwhile, the offsets between the bursts and host centers are about 4 kpc and 29 kpc for FRB 180924 and FRB 190523, respectively. These properties of hosts are similar to short gamma-ray bursts \textbf{(SGRBs)}, which are produced by mergers of binary neutron star (BNS) or neutron star-black hole (NS-BH). Therefore, the NSs powering FRBs may be formed in BNS mergers. In this paper, we study the BNS merger rates, merger times, and predict their most likely merger locations for different types of host galaxies using population synthesis method. We find that the BNS merger channel is consistent with the recently reported offsets of FRB 180924 and FRB 190523. The offset distribution of short GRBs is well reproduced by population synthesis using galaxy model which is similar to GRB hosts. The event rate of FRBs (including non-repeating and repeating), is larger than those of BNS merger and short GRBs, which requires a large fraction of observed FRBs emitting several bursts. Using curvature radiation by bunches in NS magnetospheres, we also predict the observational properties of FRBs from BNS mergers, including the dispersion measure, and rotation measure. At late times ($t\geq1$yr), the contribution to dispersion measure and rotation measure from BNS merger ejecta could be neglected.

We apply the inverse-Gaussianization method proposed in \citealt{arXiv:1607.05007} to fast produce weak lensing convergence maps and investigate the peak statistics, including the peak height counts and peak steepness counts, in these mocks. We find that the distribution of peak height and steepness is in good agreement with the simulation. The difference is $\lesssim 20\%$ for these peak statistics in the maps at source redshift $z_s=1$. Besides, the loss of off-diagonal elements in peak covariance motivates us to consider the super sample variance in weak lensing peak statistics. We propose correction methods to effectively recover the (anti-)correlation among different bins by adding scatters in the mean value of these mocks. Finally, as an example of the application, we adopt the improved inverse-Gaussianization method with the above improvement to fast generate 40,000 mocks to calculate precision matrices between the power spectrum and peak statistics.

We propose a new velocity reconstruction method based on the displacement estimation by recently developed methods. The velocity is first reconstructed by transfer functions in Lagrangian space and then mapped into Eulerian space. High resolution simulations are used to test the performance. We find that the new reconstruction method outperforms the standard velocity reconstruction in the sense of better cross-correlation coefficient, less velocity misalignment and smaller amplitude difference. We conclude that this new method has the potential to improve the large-scale structure sciences involving a velocity reconstruction, such as kinetic Sunyaev-Zel'dovich measurement and supernova cosmology.

Using the integral field unit (IFU) data from Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, we study the kinematics of gas and stellar components in an edge-on Seyfert 2 galaxy, SDSS J171359.00+333625.5, with X-shaped bi-conical outflows. The gas and stars therein are found to be counter-rotating, indicating that the collision between the inner and external gas might be an effective way to dissipate the angular momentum, which leads to remarkable gas accretion into the galaxy center. Large [OIII]$\lambda$5007 equivalent width and AGN-like line ratio in the large bi-conical region suggest that the gas is ionized by the central AGN. The gas velocity in the bi-cone region shows that ionized gas is receding relative to the galaxy center, which could be the joint effect of inflows, outflows and disk rotation. We are probably witnessing the case where a great amount of gas in the disk is being efficiently accreted into the central black hole, and the AGN-driven galactic winds are blown out along the bi-cone. The kinematics of oxygen, including rotation velocity and velocity dispersion, is different from other elements, like hydrogen, nitrogen and sulfur. The rotation velocity estimated from oxygen is slower than from other elements. The velocity dispersion of other elements follows galactic gravitational potential, while the velocity dispersion of oxygen stays roughly constant along the galactic major-axis. The further advanced observations, e.g. of cold gas or with an IFU of higher spatial resolution, are required to better understand this object.

Recently, GRB 170817A was confirmed to be associated with GW 170817, which was produced by a neutron star - neutron star merger. It indicates that at least some short gamma-ray bursts come from binary neutron star mergers. Theoretically, it is widely accepted that short gamma-ray bursts can be produced by two distinctly different mechanisms, binary neutron star mergers and neutron star - black hole mergers. These two kinds of bursts should be different observationally due to their different trigger mechanisms. Motivated by this idea, we collect a universal data set constituted of 51 short gamma-ray bursts observed by $Swift$/BAT, among which 14 events have extended emission component. We study the observational features of these 51 events statistically. It is found that our samples are consisted of two distinct groups. They clearly show a bimodal distribution when their peak photon fluxes at 15-150 keV band are plotted against the corresponding fluences. Most interestingly, all the 14 short bursts with extended emission lie in a particular region. When the fluences are plotted against the burst durations, short bursts with extended emission again tend to concentrate in the long duration segment. These features strongly indicate that short gamma-ray bursts really may come from two distinct types of progenitors. We argue that those short gamma-ray bursts with extended emission come from the coalescence of neutron stars, while the short gamma-ray bursts without extended emission come from neutron star - black hole mergers.

The precise measurement of the spectrum of protons, the most abundant component of the cosmic radiation, is necessary to understand the source and acceleration of cosmic rays in the Milky Way. This work reports the measurement of the cosmic ray proton fluxes with kinetic energies from 40 GeV to 100 TeV, with two and a half years of data recorded by the DArk Matter Particle Explorer (DAMPE). This is the first time an experiment directly measures the cosmic ray protons up to ~100 TeV with a high statistics. The measured spectrum confirms the spectral hardening found by previous experiments and reveals a softening at ~13.6 TeV, with the spectral index changing from ~2.60 to ~2.85. Our result suggests the existence of a new spectral feature of cosmic rays at energies lower than the so-called knee, and sheds new light on the origin of Galactic cosmic rays.

The nature of fast radio bursts (FRBs), which occurs on millisecond time scales in the radio band, has not been well-understood. Among their unknown observational properties are their broadband spectra and persistent and transient multi-wavelength counterparts. Well-localized FRBs provide us the opportunity to address these issues in archival observations. We have performed searches for 15-150 keV hard X-ray bursts on time scales as short as millisecond in the direction of the repeating FRB 121102 (with a spacial resolution of a few arcminutes) in the archival Swift/BAT data during the period between October 2016 and September 2017. We have found no significant (5 $\sigma$) hard X-ray bursts in the direction of the repeating FRB. We have derived an upper limit of the hard X-ray (15--150 keV) flux of any X-ray bursts on 1 ms time scale of around $1.01 \times 10^{-7}$erg cm$^{-2}$s$^{-1}$, if assuming a photo-index of 2 for potential X-ray flares in X-ray band. A plausible scenario for the repeating FRB as being associated with \emph{magnetar giant flare} is still far below the upper limit.

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

PSR B1259-63/LS 2883 is an elliptical pulsar/Be star binary and emits broadband emissions from radio to TeV $\gamma$-rays. The massive star possesses an equatorial disc, which is inclined with the orbital plane of the pulsar. The non-thermal emission from the system is believed to be produced by the pulsar wind shock and the double-peak profiles in the X-ray and TeV $\gamma$-ray light curves are related to the phases of the pulsar passing through the disc region of the star. In this paper, we investigate the interactions between the pulsar wind and stellar outflows, especially with the presence of the disc, and present a multi-wavelength modelling of the emission from this system. We show that the double-peak profiles of X-ray and TeV $\gamma$-ray light curves are caused by the enhancements of the magnetic field and the soft photons at the shock during the disc passages. As the pulsar is passing through the equatorial disc, the additional pressure of the disc pushes the shock surface closer to the pulsar, which causes the enhancement of magnetic field in the shock, and thus increases the synchrotron luminosity. The TeV $\gamma$-rays due to the inverse-Compton (IC) scattering of shocked electrons with seed photons from the star is expected to peak around periastron which is inconsistent with observations. However, the shock heating of the stellar disc could provide additional seed photons for IC scattering during the disc passages, and thus produces the double-peak profiles as observed in the TeV $\gamma$-ray light curve. Our model can possibly be examined and applied to other similar gamma-ray binaries, such as PSR J2032+4127/MT91 213, HESS J0632+057, and LS I+61$^{\circ}$303.

The DArk Matter Particle Explorer (DAMPE), a satellite-based cosmic ray and gamma-ray detector, was launched on December 17, 2015, and began its on-orbit operation on December 24, 2015. In this work we document the on-orbit calibration procedures used by DAMPE and report the calibration results of the Plastic Scintillator strip Detector (PSD), the Silicon-Tungsten tracKer-converter (STK), the BGO imaging calorimeter (BGO), and the Neutron Detector (NUD). The results are obtained using Galactic cosmic rays, bright known GeV gamma-ray sources, and charge injection into the front-end electronics of each sub-detector. The determination of the boundary of the South Atlantic Anomaly (SAA), the measurement of the live time, and the alignments of the detectors are also introduced. The calibration results demonstrate the stability of the detectors in almost two years of the on-orbit operation.

The accretion-induced collapse (AIC) of a white dwarf in a binary with a nondegenerate companion can sometimes lead to the formation of a rapidly rotating and highly magnetized neutron star (NS). The spin-down of this NS can drive a powerful pulsar wind (PW) and bring out some detectable multi-wavelength emissions. On the one hand, the PW can evaporate the companion in a few days to form a torus surrounding the NS. Then, due to the blockage of the PW by the torus, a reverse shock can be formed in the wind to generate intense hard X-rays. This emission component disappears in a few weeks' time, after the torus is broken down at its inner boundary and scoured into a very thin disk. On the other hand, the interaction between the PW with an AIC ejecta can lead to a termination shock of the wind, which can produce a long-lasting soft X-ray emission component. In any case, the high-energy emissions from deep inside the system can be detected only after the AIC ejecta becomes transparent for X-rays. Meanwhile, by absorbing the X-rays, the AIC ejecta can be heated effectively and generate a fast-evolving and luminous ultraviolet (UV)/optical transient. Therefore, the predicted hard and soft X-ray emissions, associated by an UV/optical transient, provide a clear observational signature for identifying AIC events in current and future observations (e.g., AT 2018cow).

Super-Chandrasekhar remnants of double white dwarf mergers could sometimes collapse into a rapidly rotating neutron star (NS), accompanying with a mass ejection of a few times $0.01M_{\odot}$. Bright optical transient emission can be produced by the ejecta due to heating by radioactivities and particularly by energy injection from the NS. Since the merger remnants before collapse resemble a star evolving from the asymptotic giant branch phase to the planetary nebula phase, an intense dusty wind is considered to be driven about several thousand years ago before the collapse and surround the remnant at large radii. Therefore, the optical transient emission can be somewhat absorbed and scattered by the dusty wind, which can suppress the peak emission and cause a scattering plateau in optical light curves. Several years later, as the ejecta finally catches up with the wind material, the shock interaction between them can further give rise to a detectable radio transient emission on a timescale of tens of days. Discovery of and observations to such dust-affected optical transients and shock-driven radio transients can help to explore the nature of super-Chandrasekhar merger remnants and as well as the density and type ratios of double white dwarf systems, which is beneficial in assessing their gravitational wave contributions.

During the in-spiral stage of a compact binary, a wind bubble could be blown into interstellar medium, if the electromagnetic radiation due to the binary orbital motion is strong enough. Therefore, short gamma-ray bursts (SGRBs) due to double neutron star mergers would in principle happen in a wind bubble environment, which can influence the propagation of the SGRB jet and the consequent afterglow emission. By calculating the dynamics and synchrotron radiation of the jet-driven external shock, we reveal that an abrupt jump could appear in the afterglow light curves of SGRBs and the observational time of the jump is dependent on the viewing angle. This light curve jump provides an observational signature to constrain the radius of the wind bubble and then the power of the binary electromagnetic radiation, by combining with gravitational wave detection.

The evaporation of primordial black holes with a mass in the $1\ {\rm gram}\lesssim M_{\rm PBH}\lesssim$1000 kg range can lead to the production of dark matter particles of almost any mass in the range $0.1\ {\rm MeV}\lesssim m_{\rm DM}\lesssim 10^{18}$ GeV with the right relic density at very early times, $\tau\lesssim 10^{-10}$ s. We calculate, as a function of the primordial black holes mass and initial abundance, the combination of dark matter particle masses and number of effective dark degrees of freedom leading to the right abundance of dark matter today, whether or not evaporation stops around the Planck scale. In addition, since black hole evaporation can also lead to the production of a baryon asymmetry, we calculate where dark matter production and baryogenesis can concurrently happen, under a variety of assumptions: baryogenesis via grand unification boson decay, via leptogenesis, or via asymmetric co-genesis of dark matter and ordinary matter. Finally, we comment on possible ways to test this scenario.

Over the past decade, many rapidly evolving optical transients (REOTs), whose rise and decline timescales are significantly shorter than those of canonical supernovae (SNe), have been discovered and studied. Some REOTs have high peak luminosity ($\gtrsim 10^{43}$ erg s$^{-1}$), disfavoring the radioactivity-powered-SN model that has been widely adopted to explain normal SNe. In this paper, we study three luminous REOTs (PS1-10bjp, PS1-11bbq, and PS1-13ess) and use a model involving magnetar energy input to fit their bolometric light curves and temperature evolution. We find that core-collapse SNe (CCSNe) powered by magnetars with $P_0 \approx 18$--34 ms and $B_p \approx$ (2.5--5.8) $\times10^{15}$ G can reproduce their bolometric light curves as well as the temperature evolution. The inferred values of ejecta mass are $\sim 0.40$--0.46 M$_\odot$, slightly smaller than that of the Type Ic SN 1994I, indicating that they can be explained by the magnetar-powered stripped CCSN model.

Low density regions are less affected by the nonlinear structure formation and baryonic physics. They are ideal places for probing the nature of dark energy, a possible explanation for the cosmic acceleration. Unlike void lensing, which requires identifications of individual voids, we study the stacked lensing signals around the low-density-positions (LDP), defined as places that are devoid of foreground bright galaxies in projection. The method allows a direct comparison with numerical results by drawing correspondence between the bright galaxies with halos. It leads to lensing signals that are significant enough for differentiating several dark energy models. In this work, we use the CFHTLenS catalogue to define LDPs, as well as measuring their background lensing signals. We consider several different definitions of the foreground bright galaxies (redshift range \& magnitude cut). Regarding the cosmological model, we run six simulations: the first set of simulations have the same initial conditions, with $\rm{w_{de}=-1,-0.5,-0.8,-1.2}$; the second set of simulations include a slightly different $\Lambda$CDM model and a w(z) model from \cite{2017NatAs...1..627Z}. The lensing results indicate that the models with $\rm{w_{de}=-0.5,-0.8}$ are not favored, and the other four models all achieve comparable agreement with the data.

The formation of neutron stars (NSs), both from collapses of massive stars and mergers of compact objects, can be usually indicated by bright transients emitted from explosively-ejected material. In particular, if the newborn NSs can rotate at a millisecond period and have a sufficiently high magnetic field, then the spin-down of the NSs would provide a remarkable amount of energy to the emitting material. As a result, super-luminous supernovae could be produced in the massive stellar collapse cases, while some unusual fast evolving and luminous optical transients could arise from the cases of NS mergers and accretion-induced collapses of white dwarfs. In all cases, if the dipolar magnetic fields of the newborn NSs can be amplified to be as high as $10^{15}$ G, a relativistic jet could be launched and then a gamma-ray burst can be produced as the jet successfully breaks out from the surrounding nearly-isotropic ejected material.

PSR B0919+06 is known for its abnormal emission phenomenon, where the pulse emission window occasionally shifts progressively in longitude and returns afterwards. The physical mechanism behind this phenomenon is still under investigation. In this paper, we present our ultra-wideband observation of this pulsar using the Five-hundred-meter Aperture Spherical radio Telescope (FAST), with simultaneous measurements in the frequency ranges 280-780 MHz and 1250-1550 MHz. We have identified three abnormal events, each of which becomes less apparent as the frequency decreases. At 1400 MHz, the averaged profile slightly shifted after the first and third abnormal events, implying a relationship between abnormal event and profile variation. We also found a linear trend in the left-edge position of the averaged profiles between the first and third event as well as after the third event, suggesting the existence of a slow-drifting mode between the two major events. The second event has a comparatively small shift in phase and is thus categorized as a `small flare state'. During the third event, a sequence of approximately nine pulses was seen to significantly weaken in all frequency bands, likely associated with the pseudo-nulling observed at 150 MHz. A three-component de-composition analysis of the normal averaged profiles shows that the trailing component is dominant at our observing frequencies, while the centre component has a comparatively steeper spectrum. We found the overall flux density in an abnormal event to slightly differ from that in an ordinary state, and the difference shows a frequency dependence. A comparison of the normal, abnormal and dimmed averaged profile indicates that the leading component is likely to be stable in all states.

With the largest dish Five-hundred-meter Aperture Spherical radio Telescope (FAST), both the mean and single pulses of PSR B2016$+$28, especially including the single-pulse structure, are investigated in detail in this study. The mean pulse profiles at different frequencies can be well fitted in a conal model, and the peak separation of intensity-dependent pulse profiles increases with intensity. The integrated pulses are obviously frequency dependent (pulse width decreases by $\sim\,20\%$ as frequency increases from 300 MHz to 750 MHz), but the structure of single pulses changes slightly (the corresponding correlation scale decreases by only $\sim\,1\%$). This disparity between mean and single pulses provides independent evidence for the existence of the RS-type vacuum inner gap, indicating a strong bond between particles on the pulsar surface. Diffused drifting sub-pulses are analyzed. The results show that the modulation period along pulse series ($P_3$) is positively correlated to the separation between two adjacent sub-pulses ($P_2$). This correlation may hint a rough surface on the pulsar, eventually resulting in the irregular drift of sparks. All the observational results may have significant implications in the dynamics of pulsar magnetosphere and are discussed extensively in this paper.

Rotating radio transients (RRATs) are peculiar astronomical objects whose emission mechanism remains under investigation. In this paper, we present observations of three RRATs, J1538+2345, J1854+0306 and J1913+1330, observed with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Specifically, we analyze the mean pulse profiles and temporal flux density evolutions of the RRATs. Owing to the high sensitivity of FAST, the derived burst rates of the three RRATs are higher than those in previous reports. RRAT J1854+0306 exhibited a time-dynamic mean pulse profile, whereas RRAT J1913+1330 showed distinct radiation and nulling segments on its pulse intensity trains. The mean pulse profile variation with frequency is also studied for RRAT J1538+2345 and RRAT J1913+1330, and the profiles at different frequencies could be well fitted with a cone-core model and a conal-beam model, respectively.

21 cm intensity mapping has emerged as a promising technique to map the large-scale structure of the Universe, at redshifts $z$ from 1 to 10. Unfortunately, many of the key cross-correlations with the cosmic microwave background and photo-$z$ galaxies have been thought to be impossible due to the foreground contamination for radial modes with small wave numbers. In this paper, we apply tidal reconstruction to the simulated 21 cm fields and recover the lost large-scale radial modes successfully. We estimate the detectability of the cross-correlation signals and find they can be detected at high significance with current 21 cm experiments. The tidal field reconstruction method opens up a new set of possibilities to probe the Universe and is extremely valuable not only for 21 cm surveys but also for cosmic microwave background and photometric-redshift observations.

One of the main purposes of the DArk Matter Particle Explorer (DAMPE) is to measure the cosmic ray nuclei up to several tens of TeV or beyond, whose origin and propagation remains a hot topic in astrophysics. The Plastic Scintillator Detector (PSD) on top of DAMPE is designed to measure the charges of cosmic ray nuclei from H to Fe and serves as a veto detector for discriminating gamma-rays from charged particles. We propose in this paper a charge reconstruction procedure to optimize the PSD performance in charge measurement. Essentials of our approach, including track finding, alignment of PSD, light attenuation correction, quenching and equalization correction are described detailedly in this paper after a brief description of the structure and operational principle of the PSD. Our results show that the PSD works very well and almost all the elements in cosmic rays from H to Fe are clearly identified in the charge spectrum.

DArk Matter Particle Explorer (DAMPE) is a space-borne apparatus for detecting the high-energy cosmic-rays like electrons, $\gamma$-rays, protons and heavy-ions. Plastic Scintillator Detector (PSD) is the top-most sub-detector of the DAMPE. The PSD is designed to measure the charge of incident high-energy particles and it also serves as a veto detector for discriminating $\gamma$-rays from charged particles. In this paper, PSD on-orbit calibration procedure is described, which includes five steps of pedestal, dynode correlation, response to minimum-ionizing particles (MIPs), light attenuation function and energy reconstruction. A method for reconstructing the charge of incident high energy cosmic-ray particles is introduced. The detection efficiency of each PSD strip is verified to be above 99.5%, the total efficiency of the PSD for charged particles is above 99.99%.

The galaxy intrinsic alignment (IA) is a major challenge of weak lensing cosmology. To alleviate this problem, Zhang (2010, MNRAS, 406, L95) proposed a self-calibration method, independent of IA modeling. This proposal relies on several scaling relations between two-point clustering of IA and matter/galaxy fields, which were previously only tested with analytical IA models. In this paper, these relations are tested comprehensively with an $N$-body simulation of $3072^3$ simulation particles and boxsize 600 $h^{-1} \, \mathrm{Mpc}$. They are verified at the accuracy level of $\mathcal{O}(1)\%$ over angular scales and source redshifts of interest. We further confirm that these scaling relations are generic, insensitive to halo mass, weighting in defining halo ellipticities, photo-$z$ error, and misalignment between galaxy ellipticities and halo ellipticities. We also present and verify three new scaling relations on the B-mode IA. These results consolidate and complete the theory side of the proposed self-calibration technique.

The Plastic Scintillator Detector (PSD) of the DArk Matter Particle Explorer (DAMPE) is designed to measure cosmic ray charge (Z) and to act as a veto detector for gamma-ray identification. In order to fully exploit the charge identification potential of the PSD and to enhance its capability to identify the gamma ray events, we develop a PSD detector alignment method. The path length of a given track in the volume of a PSD bar is derived taking into account the shift and rotation alignment corrections. By examining energy spectra of corner-passing events and fully contained events, position shifts and rotations of all PSD bars are obtained, and found to be on average about 1mm and 0.0015 radian respectively. To validate the alignment method, we introduce the artificial shifts and rotations of PSD bars in the detector simulation. These shift and rotation parameters can be recovered successfully by the alignment procedure. As a result of the PSD alignment procedure, the charge resolution of the PSD is improved from $4\%$ to $8\%$ depending on the nuclei.

The groundbreaking discovery of the optical transient AT2017gfo associated with GW170817 opens a unique opportunity to study the physics of double neutron star (NS) mergers. We argue that the standard interpretation of AT2017gfo as being powered by radioactive decays of r-process elements faces the challenge of simultaneously accounting for the peak luminosity and peak time of the event, as it is not easy to achieve the required high mass, and especially the low opacity of the ejecta required to fit the data. A plausible solution would be to invoke an additional energy source, which is probably provided by the merger product. We consider energy injection from two types of the merger products: (1) a post-merger black hole powered by fallback accretion; and (2) a long-lived NS remnant. The former case can only account for the early emission of AT2017gfo, with the late emission still powered by radioactive decay. In the latter case, both early- and late-emission components can be well interpreted as due to energy injection from a spinning-down NS, with the required mass and opacity of the ejecta components well consistent with known numerical simulation results. We suggest that there is a strong indication that the merger product of GW170817 is a long-lived (supramassive or even permanently stable), low magnetic field NS. The result provides a stringent constraint on the equations of state of NSs.

A long-standing issue in peculiar velocity cosmology is whether the halo/galaxy velocity bias $b_v=1$ at large scale. The resolution of this important issue must resort to high precision cosmological simulations. However, this is hampered by another long-standing `sampling artifact' problem in volume weighted velocity measurement. We circumvent this problem with a hybrid approach. We first measure statistics free of sampling artifact, then link them to volume weighted statistics in theory, finally solve for the velocity bias. $b_v$ determined by our method is not only free of sampling artifact, but also free of cosmic variance. We apply this method to a $\Lambda$CDM N-body simulation of $3072^3$ particles and $1200 Mpc/{\rm h}$ box size. For the first time, we determine the halo velocity bias to $0.1\%$-$1\%$ accuracy. Our major findings are as follows: (1) $b_v\neq 1$ at $k>0.1 h/{\rm Mpc}$. The deviation from unity ($|b_v-1|$) increases with $k$. Depending on halo mass and redshift, it may reach $\mathcal{O}(0.01)$ at $k=0.2 h/{\rm Mpc}$ and $\mathcal{O}(0.05)$ at $k\sim 0.3 h/{\rm Mpc}$. The discovered $b_v\neq 1$ has statistically significant impact on structure growth rate measurement by spectroscopic redshift surveys, including DESI, Euclid and SKA. (2) Both the sign and the amplitude of $b_v-1$ depend on mass and redshift. These results disagree with the peak model prediction in that $b_v$ has much weaker deviation from unity, varies with redshift, and can be bigger than unity. (3) Most of the mass and redshift dependences can be compressed into a single dependence on the halo density bias. Based on this finding, we provide an approximate two-parameter fitting formula.

In paper I of this series (Yang et al. 2017, ApJ), we proposed an analytical method of blind separation ({\bf ABS}) to extract the cosmic magnification signal in galaxy number distribution and reconstruct the weak lensing power spectrum. Here we report a new version of the ABS method, with significantly improved performance. This version is characterized by a shift parameter $\mathcal{S}$, with the special case of $\mathcal{S}=0$ corresponding to the original ABS method. We have tested this new version, compared to the previous one, and confirmed its supreme performance in all investigated situations. Therefore it supersedes the previous version. The proof of concept studies presented in this paper demonstrate that it may enable surveys such as LSST and SKA to reconstruct the lensing power spectrum at $z\simeq 1$ with $1\%$ accuracy. We will test with more realistic simulations to verify its applicability in real data.

The successful joint observation of the gravitational wave event GW170817 and its multi-wavelength electromagnetic counterparts first enables human to witness a definite merger event of two neutron stars (NSs). This historical event confirms the origin of short-duration gamma-ray bursts (GRBs), and in particular, identifies the theoretically-predicted kilonova phenomenon that is powered by radioactive decays of $r$-process heavy elements. However, whether a long-lived remnant NS could be formed during this merger event remains unknown, although such a central engine has been suggested by afterglow observations of some short-duration GRBs. By invoking this long-lived remnant NS, we here propose a model of hybrid energy sources for the kilonova AT2017gfo associated with GW 170817. While the early emission of AT2017gfo is still powered radioactively as usually suggested, its late emission is primarily caused by delayed energy injection from the remnant NS. In our model, only one single opacity is required and an intermediate value of $\kappa\simeq0.97\,\rm cm^2g^{-1}$ is revealed, which could be naturally provided by lanthanide-rich ejecta that is deeply ionized by the emission from a wind of the NS. These self-consistent results indicate that a long-lived remnant NS, which must own a very stiff equation of state, had been formed during the merger event of GW170817. This provides a very stringent constraint on the strong interaction in nuclear-quark matter. It is further implied that such GW events could provide a probe of the early spin and magnetic evolutions of NSs, e.g., the burying of surface magnetic fields.

The central engines of some superluminous supernovae (SLSNe) are generally suggested to be newly born fast rotating magnetars, which spin down mainly through magnetic dipole radiation and gravitational wave emission. We calculate the magnetar-powered SLSNe light curves (LCs) with the tilt angle evolution of newly born magnetars involved. We show that, depending on the internal toroidal magnetic fields ${\bar B}_{\rm t}$, the initial spin periods $P_{\rm i}$, and the radii $R_{\rm DU}$ of direct Urca (DU) cores of newly born magnetars, as well as the critical temperature $T_{\rm c}$ for $^3P_2$ neutron superfluidity, bumps could appear in the SLSNe LCs after the maximum lights when the tilt angles grow to $\pi/2$. The value of $T_{\rm c}$ determines the arising time and the relative amplitude of a bump. The quantity $R_{\rm DU}$ can affect the arising time and the luminosity of a bump, as well as the peak luminosity of a LC. Moreover, it is interesting that a stronger ${\bar B}_{\rm t}$ will lead to both a brighter peak and a brighter bump in a LC. While keeping other quantities unchanged, the bump in the LC disappears for the magnetar with smaller $P_{\rm i}$. We suggest that, once the SLSNe LCs with such kinds of bumps are observed, by fitting these LCs with our model, not only $B_{\rm d}$ and $P_{\rm i}$ of newly born magnetars but also the crucial physical quantities ${\bar B}_{\rm t}$, $R_{\rm DU}$, and $T_{\rm c}$ could be determined. Nonobservation of SLSNe LCs with such kinds of bumps hitherto may already put some (\textit{though very rough}) constraints on ${\bar B}_{\rm t}$, $P_{\rm i}$, $R_{\rm DU}$, and $T_{\rm c}$. Therefore, observation of SLSNe LCs may provide a new approach to probe the physics of newly born magnetars.

The brightness of the multi-wavelength afterglow of GRB 170817A is increasing unexpectedly even $\sim$ 160 days after the associated gravitational burst. Here we suggest that the brightening can be caused by a late-time energy injection process. We use an empirical expression to mimic the evolution of the injection luminosity, which is consisted of a power-law rising phase and a power-law decreasing phase. It is found that the power-law indices of the two phases are $0.92$ and $-2.8$, respectively, with the peak time of the injection being $\sim 110$ days. The energy injection could be due to some kind of accretion, with the total accreted mass being $\sim 0.006 M_\odot$. However, normal fall-back accretion, which usually lasts for a much shorter period, cannot provide a natural explanation. Our best-fit decaying index of $-2.8$ is also at odds with the expected value of $-5/3$ for normal fall-back accretion. Noting that the expansion velocities of the kilonova components associated with GW170817 is $0.1-0.3\,c$, we argue that there should also be some ejecta with correspondingly lower velocities during the coalescence of the double neutron star system. They are bound by the gravitational well of the remnant central compact object and might be accreted at a timescale of about 100 days, providing a reasonable explanation for the energy injection. Detailed studies on the long-lasting brightening of GRB 170817A thus may provide useful information on the matter ejection during the merger process of binary neutron stars.

The intra-cluster light (ICL) in observations is usually identified through the surface brightness limit method. In this paper, for the first time we produce the mock images of galaxy groups and clusters using a cosmological hydro- dynamical simulation, to investigate the ICL fraction and focus on its dependence on observational parameters, e.g., the surface brightness limit (SBL), the effects of cosmological redshift dimming, point spread function and CCD pixel size. Detailed analyses suggest that the width of point spread function has a significant effect on the measured ICL fraction, while the relatively small pixel size shows almost no influence. It is found that the measured ICL fraction depends strongly on the SBL. At a fixed SBL and redshift, the measured ICL fraction decreases with increasing halo mass, while with a much faint SBL, it does not depend on halo mass at low redshifts. In our work, the measured ICL fraction shows clear dependence on the cosmological redshift dimming effect. It is found that there are more mass locked in ICL component than light, suggesting that the use of a constant mass-to-light ratio at high surface brightness levels will lead to an underestimate of ICL mass. Furthermore, it is found that the radial profile of ICL shows a characteristic radius which is almost independent of halo mass. The current measurement of ICL from observations has a large dispersion due to different methods, and we emphasize the importance of using the same definition when observational results are compared with the theoretical predictions.

Fast radio bursts (FRBs) are usually suggested to be associated with mergers of compact binaries consisting of white dwarfs (WDs), neutron stars (NSs), or black holes (BHs). We test these models by fitting the observational distributions in both redshift and isotropic energy of 22 Parkes FRBs, where, as usual, the rates of compact binary mergers (CBMs) are connected with cosmic star formation rates by a power-law distributed time delay. It is found that the observational distributions can well be produced by the CBM model with a characteristic delay time from several ten to several hundred Myr and an energy function index $1.2\lesssim\gamma\lesssim1.7$, where a tentative fixed spectral index $\beta=0.8$ is adopted for all FRBs. Correspondingly, the local event rate of FRBs is constrained to $(3-6)\times10^4f_{\rm b}^{-1}(\mathcal T/270\rm s)^{-1}(\mathcal A/2\pi)^{-1}\rm ~Gpc^{-3}yr^{-1}$ for an adopted minimum FRB energy of $E_{\min}=3\times10^{39}$ erg, where $f_{\rm b}$ is the beaming factor of the radiation, $\mathcal T$ is the duration of each pointing observation, and $\mathcal A $ is the sky area of the survey. This event rate, about an order of magnitude higher than the rates of NS-NS/NS-BH mergers, indicates that the most promising origin of FRBs in the CBM scenario could be mergers of WD-WD binaries. Here a massive WD could be produced since no FRB was found to be associated with a type Ia supernova. Alternatively, if actually all FRBs can repeat on a timescale much longer than the period of current observations, then they could also originate from a young active NS that forms from relatively rare NS-NS mergers and accretion-induced collapses of WD-WD binaries.

The state of supranuclear matter in compact star remains puzzling, and it is argued that pulsars could be strangeon stars. What if binary strangeon stars merge? This kind of merger could result in the formation of a hyper-massive strangeon star, accompanied by bursts of gravitational waves and electromagnetic radiation (and even strangeon kilonova explained in the paper). The tidal polarizability of binary strangeon stars is different from that of binary neutron stars, because a strangeon star is self-bound on surface by fundamental strong force while a neutron star by the gravity, and their equations of state are different. Our calculation shows that the tidal polarizability of merging binary strangeon stars is favored by GW170817. Three kinds of kilonovae (i.e., of neutron, quark and strangeon) are discussed, and the light curve of the kilonova AT 2017gfo following GW170817 could be explained by considering the decaying strangeon nuggets and remnant star spin-down. Additionally, the energy ejected to the fireball around the nascent remnant strangeon star, being manifested as a Gamma-ray burst (GRB), is calculated. It is found that, after a promote burst, an X-ray plateau could follow in a timescale of $10^{2-3}$ s. Certainly, the results could be tested also by further observational synergies between gravitational wave detectors (e.g., aLIGO) and X-ray telescopes (e.g., Chinese HXMT and eXTP), and especially if the detected gravitational wave form is checked by peculiar equation of state provided by the numerical relativistical simulation.

We apply nonlinear reconstruction to the dark matter density field in redshift space and solve for the nonlinear mapping from the initial Lagrangian position to the final redshift space position. The reconstructed anisotropic field inferred from the nonlinear displacement correlates with the linear initial conditions to much smaller scales than the redshift space density field. The number of linear modes in the density field is improved by a factor of 30-40 after reconstruction. We thus expect this reconstruction approach to substantially expand the cosmological information including baryon acoustic oscillations and redshift space distortions for dense low-redshift large scale structure surveys including for example SDSS main sample, DESI BGS, and 21 cm intensity mapping surveys.

The gravitational coupling of a long-wavelength tidal field with small-scale density fluctuations leads to anisotropic distortions of the locally measured small-scale matter correlation function. Since the local correlation function is known to be statistically isotropic in the absence of such tidal interactions, the tidal distortions can be used to reconstruct the long-wavelength tidal field and large-scale density field in analogy with the cosmic microwave background lensing reconstruction. In this paper we present the theoretical framework of cosmic tidal reconstruction and test the reconstruction in numerical simulations. We find that the density field on large scales can be reconstructed with good accuracy and the cross-correlation coefficient between the reconstructed density field and the original density field is greater than 0.9 on large scales $(k\lesssim0.1\ h/\mathrm{Mpc})$, with the filter scale $\sim1.25\ \mathrm{Mpc}/h$. This is useful in the 21cm intensity mapping survey, where the long-wavelength radial modes are lost due to a foreground subtraction process.

We present a direct approach to nonparametrically reconstruct the linear density field from an observed nonlinear map. We solve for the unique displacement potential consistent with the nonlinear density and positive definite coordinate transformation using a multigrid algorithm. We show that we recover the linear initial conditions up to the nonlinear scale ($r_{\delta_r\delta_L}>0.5$ for $k\lesssim1\ h/\mathrm{Mpc}$) with minimal computational cost. This reconstruction approach generalizes the linear displacement theory to fully nonlinear fields, potentially substantially expanding the baryon acoustic oscillations and redshift space distortions information content of dense large scale structure surveys, including for example SDSS main sample and 21cm intensity mapping initiatives.

High energy cosmic ray electrons plus positrons (CREs), which lose energy quickly during their propagation, provide an ideal probe of Galactic high-energy processes and may enable the observation of phenomena such as dark-matter particle annihilation or decay. The CRE spectrum has been directly measured up to $\sim 2$ TeV in previous balloon- or space-borne experiments, and indirectly up to $\sim 5$ TeV by ground-based Cherenkov $\gamma$-ray telescope arrays. Evidence for a spectral break in the TeV energy range has been provided by indirect measurements of H.E.S.S., although the results were qualified by sizeable systematic uncertainties. Here we report a direct measurement of CREs in the energy range $25~{\rm GeV}-4.6~{\rm TeV}$ by the DArk Matter Particle Explorer (DAMPE) with unprecedentedly high energy resolution and low background. The majority of the spectrum can be properly fitted by a smoothly broken power-law model rather than a single power-law model. The direct detection of a spectral break at $E \sim0.9$ TeV confirms the evidence found by H.E.S.S., clarifies the behavior of the CRE spectrum at energies above 1 TeV and sheds light on the physical origin of the sub-TeV CREs.

Photo-z error is one of the major sources of systematics degrading the accuracy of weak lensing cosmological inferences. Zhang et al. (2010) proposed a self-calibration method combining galaxy-galaxy correlations and galaxy-shear correlations between different photo-z bins. Fisher matrix analysis shows that it can determine the rate of photo-z outliers at a level of 0.01-1% merely using photometric data and do not rely on any prior knowledge. In this paper, we develop a new algorithm to implement this method by solving a constrained nonlinear optimization problem arising in the self-calibration process. Based on the techniques of fixed-point iteration and non-negative matrix factorization, the proposed algorithm can efficiently and robustly reconstruct the scattering probabilities between the true-z and photo-z bins. The algorithm has been tested extensively by applying it to mock data from simulated stage IV weak lensing projects. We find that the algorithm provides a successful recovery of the scatter rates at the level of 0.01-1%, and the true mean redshifts of photo-z bins at the level of 0.001, which may satisfy the requirements in future lensing surveys.

We apply the nonlinear reconstruction method to simulated halo fields. For halo number density $2.77\times 10^{-2}$ $(h^{-1} {\rm Mpc})^{-3}$ at $z=0$, corresponding to the SDSS main sample density, we find the scale where the noise saturates the linear signal is improved to $k\gtrsim0.36\ h {\rm Mpc}^{-1}$, a factor of $2.29$ improvement in scale, or $12$ in number of linear modes. The improvement is less for higher redshift or lower halo density. We expect this to substantially improve the BAO accuracy of dense, low redshift surveys, including the SDSS main sample, 6dFGS and 21cm intensity mapping initiatives.

The DArk Matter Particle Explorer (DAMPE), one of the four scientific space science missions within the framework of the Strategic Pioneer Program on Space Science of the Chinese Academy of Sciences, is a general purpose high energy cosmic-ray and gamma-ray observatory, which was successfully launched on December 17th, 2015 from the Jiuquan Satellite Launch Center. The DAMPE scientific objectives include the study of galactic cosmic rays up to $\sim 10$ TeV and hundreds of TeV for electrons/gammas and nuclei respectively, and the search for dark matter signatures in their spectra. In this paper we illustrate the layout of the DAMPE instrument, and discuss the results of beam tests and calibrations performed on ground. Finally we present the expected performance in space and give an overview of the mission key scientific goals.

The true ground state of hadronic matter may be strange quark matter (SQM). Consequently, the observed pulsars may actually be strange quark stars, but not neutron stars. However, proving or disproving the SQM hypothesis still remains to be a difficult problem, due to the similarity between the macroscopical characteristics of strange quark stars and neutron stars. Here we propose a hopeful method to probe the existence of strange quark matter. In the frame work of the SQM hypothesis, strange quark dwarfs and even strange quark planets can also stably exist. Noting that SQM planets will not be tidally disrupted even when they get very close to their host stars due to their extreme compactness, we argue that we could identify SQM planets by searching for very close-in planets among extrasolar planetary systems. Especially, we should keep our eyes on possible pulsar planets with orbital radius less than $\sim 5.6 \times 10^{10}$~cm and period less than $\sim 6100$~s. A thorough search in the currently detected $\sim 2950$ exoplanets around normal main sequence stars has failed to identify any stable close-in objects that meet the SQM criteria, i.e. lying in the tidal disruption region for normal matter planets. However, the pulsar planet PSR J1719-1438B, with an orbital radius of $\sim 6 \times 10^{10}$~cm and orbital period of $7837$~s, is encouragingly found to be a good candidate.

We propose an Analytical method of Blind Separation (ABS) of cosmic magnification from the intrinsic fluctuations of galaxy number density in the observed galaxy number density distribution. The ABS method utilizes the different dependences of the signal (cosmic magnification) and contamination (galaxy intrinsic clustering) on galaxy flux, to separate the two. It works directly on the measured cross galaxy angular power spectra between different flux bins. It determines/reconstructs the lensing power spectrum analytically, without assumptions of galaxy intrinsic clustering and cosmology. It is unbiased in the limit of infinite number of galaxies. In reality the lensing reconstruction accuracy depends on survey configurations, galaxy biases, and other complexities, due to finite number of galaxies and the resulting shot noise fluctuations in the cross galaxy power spectra. We estimate its performance (systematic and statistical errors) in various cases. We find that, stage IV dark energy surveys such as SKA and LSST are capable of reconstructing the lensing power spectrum at $z\simeq 1$ and $\ell\la 5000$ accurately. This lensing reconstruction only requires counting galaxies, and is therefore highly complementary to the cosmic shear measurement by the same surveys.

Binary near-Earth asteroid (65803) Didymos is the target of the proposed NASA Double Asteroid Redirection Test (DART), part of the Asteroid Impact & Deflection Assessment (AIDA) mission concept. In this mission, the DART spacecraft is planned to impact the secondary body of Didymos, perturbing mutual dynamics of the system. The primary body is currently rotating at a spin period close to the spin barrier of asteroids, and materials ejected from the secondary due to the DART impact are likely to reach the primary. These conditions may cause the primary to reshape, due to landslides, or internal deformation, changing the permanent gravity field. Here, we propose that if shape deformation of the primary occurs, the mutual orbit of the system would be perturbed due to a change in the gravity field. We use a numerical simulation technique based on the full two-body problem to investigate the shape effect on the mutual dynamics in Didymos after the DART impact. The results show that under constant volume, shape deformation induces strong perturbation in the mutual motion. We find that the deformation process always causes the orbital period of the system to become shorter. If surface layers with a thickness greater than ~0.4 m on the poles of the primary move down to the equatorial region due to the DART impact, a change in the orbital period of the system and in the spin period of the primary will be detected by ground-based measurement.

Solar system small bodies come in a wide variety of shapes and sizes, which are achieved following very individual evolutional paths through billions of years. This paper focuses on the reshaping process of rubble-pile asteroids driven by meteorite impacts. In our study, numerous possible equilibrium configurations are obtained via Monte Carlo simulation, and the structural stability of these configurations is determined via eigen analysis of the geometric constructions. The eigen decomposition reveals a connection between the cluster's reactions and the types of external disturbance. Numerical simulations are performed to verify the analytical results. The gravitational N-body code pkdgrav is used to mimic the responses of the cluster under intermittent non-dispersive impacts. We statistically confirm that the stability index, the total gravitational potential and the volume of inertia ellipsoid show consistent tendency of variation. A common regime is found in which the clusters tend towards crystallization under intermittent impacts, i.e., only the configurations with high structural stability survive under the external disturbances. The results suggest the trivial non-disruptive impacts might play an important role in the rearrangement of the constituent blocks, which may strengthen these rubble piles and help to build a robust structure under impacts of similar magnitude. The final part of this study consists of systematic simulations over two parameters, the projectile momentum and the rotational speed of the cluster. The results show a critical value exists for the projectile momentum, as predicted by theory, below which all clusters become responseless to external disturbances; and the rotation proves to be significant for it exhibits an "enhancing" effect on loose-packed clusters, which coincides with the observation that several fast-spinning asteroids have low bulk densities.

This paper presents a quantitative study of the evolution of the ejecta cloud released from a hypervelocity impact on a binary asteroid. We performed numerical simulations of the post-impact dynamics of the ejecta cloud in the framework of the current mission scenario of AIDA mission project. A grid search of launching sites of ejecta was defined over the globe of Didymoon, and considering a wide range of possible ejection speeds, we determined the dependency of ejecta fate on launching sites (projectile impact sites) and speeds. This range allows us to track all the complex cases that include different types of dynamical fates. Two major mechanisms are found to be working broadly during the post-ejection evolution of the ejecta cloud: 1) ejecta on mean motion resonance orbits with Didymoon produce long-term quasi-periodic showers onto Didymoon over at least a couple of weeks after the projectile impact, 2) ejecta on non-resonant orbits produce a rapid and high re-accretion flux. This rapid and high flux occurs just once because ejecta on such orbits leave the system unless they experience a collision during their first encounter. For the second part of this study, we performed full-scale simulations of the ejecta cloud released from 6 hypothetical impact sites. We considered two kinds of material composing Didymoon's subsurface and then combined a power-law size distribution of the ejecta with an ejection speed distribution. We find that the ejecta cloud evolution can be divided in two periods. It starts with a first violent period (<10 hr) with fast re-accretion or ejection from the system. A second period is found to be more sensitive to the launching site than the first one. During this second period, ejecta will either re-accrete or being ejected from the system, depending both on their sizes and on their average survival time in close proximity of the binary components.

We present analysis of near-infrared, high-resolution spectroscopy towards the Flat-spectrum YSO (Young Stellar Object) ESO H$\alpha$ 279a (1.5 solar mass) in the Serpens star forming region, at the distance of 429 pc. Using the Immersion GRating INfrared Spectrometer (IGRINS, R=45,000), we detect emission lines originating from the accretion channel flow, jet, and inner disk. Specifically, we identify hydrogen Brackett series recombination, [Fe II], [Fe III], [Fe IV], Ca I, Na I, H2, H2O and CO overtone emission lines. By modeling five bands of CO overtone emission lines, and the symmetric double-peaked line profile for Na I emission lines, we find that ESO H$\alpha$ 279a has an actively accreting Keplerian disk. From our Keplerian disk model, we find that Na I emission lines originate between 0.04 AU and 1.00 AU, while CO overtone emission lines are from the outer part of disk, in the range between 0.22 AU and 3.00 AU. It reveals that the neutral atomic Na gas is a good tracer of the innermost region of the actively accreting disk. We derive a mass accretion rate of 2-10x10^{-7} M_solar/yr from the measured Br_gamma emission luminosity of 1.78(+-0.31)x10^{31} erg/s.

In this paper, we report a significant recovery of the linear baryonic acoustic oscillation (BAO) signature by applying the isobaric reconstruction algorithm to the non-linear matter density field. Assuming only the longitudinal component of the displacement being cosmologically relevant, this algorithm iteratively solves the coordinate transform between the Lagrangian and Eulerian frames without requiring any specific knowledge of the dynamics. For dark matter field, it produces the non-linear displacement potential with very high fidelity. The reconstruction error at the pixel level is within a few percent, and is caused only by the emergence of the transverse component after the shell-crossing. As it circumvents the strongest non-linearity of the density evolution, the reconstructed field is well-described by linear theory and immune from the bulk-flow smearing of the BAO signature. Therefore this algorithm could significantly improve the measurement accuracy of the sound horizon scale. For a perfect large-scale structure survey at redshift zero without Poisson or instrumental noise, the fractional error is reduced by a factor of 2.7, very close to the ideal limit with linear power spectrum and Gaussian covariance matrix.

Significant undulations appear in the light curve of a recently discovered super-luminous supernova (SLSN) SN 2015bn after the first peak, while the underlying profile of the light curve can be well explained by a continuous energy supply from a central engine, possibly the spin-down of a millisecond magnetar. We propose that these undulations are caused by an intermittent pulsed energy supply, indicating an energetic flare activity of the central engine of the SLSN. Many post-burst flares were discovered during X-ray afterglow observations of Gamma-Ray Bursts (GRBs). We find that the SLSN flares described here approximately obey the empirical correlation between the luminosity and time scale of GRB flares, extrapolated to the relevant longer time scales of SLSN flares. This confirms the possible connection between these two different phenomena as previously suggested.

By fitting the bolometric light curves of 31 super-luminous supernovae (SLSNe) with the magnetar engine model, we derive the ejecta masses and magnetar parameters for these SLSNe. The lower boundary of magnetic field strengths of SLSN magnetars can be set just around the critical field strength $B_{\rm c}$ of electron Landau quantization. In more details, SLSN magnetars can further be divided into two subclasses of magnetic fields of $\sim(1-5)B_{\rm c}$ and $\sim(5-10)B_{\rm c}$, respectively. It is revealed that these two subclasses of magnetars are just associated with the slow-evolving and fast-evolving bolometric light curves of SLSNe. In comparison, the magnetars harbored in gamma-ray bursts (GRBs) and associated hypernovae are usually inferred to have much higher magnetic fields with a lower boundary about $\sim10B_{\rm c}$. This robustly suggests that it is the magnetic fields that play the crucial role in distinguishing SLSNe from GRBs/hypernovae. The rotational energy of SLSN magnetars are found to be correlated with the masses of supernova ejecta, which provides a clue to explore the nature of their progenitors. Moreover, the distribution of ejecta masses of SLSNe is basically intermediate between those of normal core-collapse supernovae and hypernovae. This could indicate an intrinsic connection among these different stellar explosions.

The similarity of the host galaxy of FRB 121102 with those of long gamma-ray bursts and Type I super-luminous supernovae suggests that this FRB could be associated with a young magnetar. By assuming the FRB emission to be produced within the magnetosphere, we derive a lower limit on the age of the magnetar, after which GHz emission is able to escape freely from the dense relativistic wind of the magnetar. Another lower limit is obtained by requiring the dispersion measure contributed by the electron/positron pair wind to be consistent with the observations of the host galaxy. Furthermore, we also derive some upper limits on the magnetar age with discussions on possible energy sources of the FRB emission and the recently-discovered persistent radio counterpart. As a result, some constraints on model parameters are addressed by reconciling the lower limits with the possible upper limits that are derived with an assumption of rotational energy source.

We have analyzed the temperature, velocity and density of H2 gas in NGC 7023 with a high-resolution near-infrared spectrum of the northwestern filament of the reflection nebula. By observing NGC 7023 in the H and K bands at R ~ 45,000 with the Immersion GRating INfrared Spectrograph (IGRINS), we detected 68 H2 emission lines within the 1" x 15" slit. The diagnostic ratios of 2-1 S(1)/1-0 S(1) is 0.41-0.56. In addition, the estimated ortho-to-para ratios (OPR) is 1.63-1.82, indicating that the H2 emission transitions in the observed region arises mostly from gas excited by UV fluorescence. Gradients in the temperature, velocity, and OPR within the observed area imply motion of the photodissociation region (PDR) relative to the molecular cloud. In addition, we derive the column density of H2 from the observed emission lines and compare these results with PDR models in the literature covering a range of densities and incident UV field intensities. The notable difference between PDR model predictions and the observed data, in high rotational J levels of v = 1, is that the predicted formation temperature for newly-formed H2 should be lower than that of the model predictions. To investigate the density distribution, we combine pixels in 1" x 1" areas and derive the density distribution at the 0.002 pc scale. The derived gradient of density suggests that NGC 7023 has a clumpy structure, including a high clump density of ~10^5 cm^-3 with a size smaller than ~5 x 10^-3 pc embedded in lower density regions of 10^3-10^4 cm^-3.

The gamma-ray binary 1FGL J1018.6-5856 is known to be composed of a massive star with a compact object with an orbital period of 16.54 days. We study mechanisms of high-energy emissions (X-ray, GeV, and TeV $\gamma$-rays) from this binary system. The paper is composed of two parts. In the first part, we re-analyze the GeV emission from 1FGL J1018.6-5856 using 7 yrs of the observational data of \textit{Fermi}/LAT. The data analysis of phase-resolved spectrum suggests that the 0.1-10 GeV emissions may contain two different components: (1) a simple power-law with exponential cutoff component in the LOW state; (2) an additional component extended to $\sim $10 GeV in the HIGH state. In the second part, we develop a model to explain the multi-wavelength emissions of this binary. Assuming that 1FGL J1018.6-5856 includes a pulsar, the GeV spectrum in the LOW state can be well explained by the two-layer outer gap model of pulsar magnetosphere. As the pulsar is orbiting around the star, the pulsar is considered to be wrapped by two kinds of termination shock: a bow shock due to the interaction between pulsar wind and stellar wind, and a back shock due to the effect of orbital motion. We proposed that X-rays are mainly produced by the synchrotron radiation at the apex of the bow shock, and TeV $\gamma$-rays are mainly produced by the inverse-Compton scattering at the back shock. The sharp peak in X-ray and TeV light curves are produced by in the shock tail due to the Doppler boosting effect at inferior-conjunction while another broad sinusoidal modulations are possibly produced by the apex of the bow shock and the back shock due to variations of shock conditions.

In this paper, we propose a new scenario in which a rapidly-rotating strongly-magnetized pulsar without any surrounding supernova ejecta produces fast radio bursts (FRBs) repeatedly via some mechanisms, and meanwhile, an ultra-relativistic electron/positron pair wind from the pulsar sweeps up its ambient dense interstellar medium, giving rise to a non-relativistic pulsar wind nebula (PWN). We show that the synchrotron radio emission from such a PWN is bright enough to account for the recently-discovered persistent radio source associated with the repeating FRB 121102 in reasonable ranges of the model parameters. In addition, our PWN scenario is consistent with the non-evolution of the dispersion measure inferred from all the repeating bursts observed in four years.

Measuring the volume-weighted peculiar velocity statistics from inhomogeneously and sparsely distributed galaxies/halos, by existing velocity assignment methods, suffers from a significant sampling artifact. As an alternative, the Kriging interpolation based on Gaussian processes was introduced and evaluated [Y. Yu, J. Zhang, Y. Jing, and P. Zhang, Phys. Rev. D 92, 083527 (2015)]. Unfortunately, the most straightforward application of Kriging does not perform better than the existing methods in the literature. In this work, we investigate two physically motivated extensions. The first takes into account of the anisotropic velocity correlations. The second introduces the nugget effect, on account of multi-streaming of the velocity field. We find that the performance is indeed improved. For sparsely sampled data [$n_P\lesssim 6\times 10^{-3}(h^{-1} {\rm Mpc})^{-3}$] where the sampling artifact is the most severe, the improvement is significant and is two-fold: 1) The scale of reliable measurement of the velocity power spectrum is extended by a factor $\sim 1.6$, and 2) the dependence on the velocity correlation prior is weakened by a factor of $\sim 2$. We conclude that such extensions are desirable for accurate velocity assignment by Kriging.

he DArk Matter Particle Explorer (DAMPE) is a general purposed satellite-borne high energy $\gamma-$ray and cosmic ray detector, and among the scientific objectives of DAMPE are the searches for the origin of cosmic rays and an understanding of Dark Matter particles. As one of the four detectors in DAMPE, the Plastic Scintillator Detector (PSD) plays an important role in the particle charge measurement and the photons/electrons separation. The PSD has 82 modules, each consists of a long organic plastic scintillator bar and two PMTs at both ends for readout, in two layers and covers an overall active area larger than 82 cm $\times$ 82 cm. It can identify the charge states for relativistic ions from H to Fe, and the detector efficiency for Z=1 particles can reach 0.9999. The PSD has been successfully launched with DAMPE on Dec. 17, 2015. In this paper, the design, the assembly, the qualification tests of the PSD and some of the performance measured on the ground have been described in detail.

We report a study of the three-dimensional (3D) outflow structure of a 15$\arcsec$ $\times$ 13$\arcsec$ area around H$_{2}$ peak 1 in Orion KL with slit-scan observations (13 slits) using the Immersion Grating Infrared Spectrograph. The datacubes, with high velocity-resolution ($\sim$ 7.5 {\kms}) provide high contrast imaging within ultra-narrow bands, and enable the detection of the main stream of the previously reported H$_{2}$ outflow fingers. We identified 31 distinct fingers in H$_{2}$ 1$-$0 S(1) $\lambda$2.122 $\micron$ emission. The line profile at each finger shows multiple-velocity peaks with a strong low-velocity component around the systemic velocity at ${\VLSR}$ = $+$8 {\kms} and high velocity emission ($|$${\VLSR}$$|$ = 45$-$135 {\kms}) indicating a typical bow-shock. The observed radial velocity gradients of $\sim$ 4 {\kms} arcsec$^{-1}$ agree well with the velocities inferred from large-scale proper motions, where the projected motion is proportional to distance from a common origin. We construct a conceptual 3D map of the fingers with the estimated inclination angles of 57$\degree$$-$74$\degree$. The extinction difference ($\Delta$$A_{\rm v}$ $>$ 10 mag) between blueshifted and redshifted fingers indicates high internal extinction. The extinction, the overall angular spread and scale of the flow argue for an ambient medium with very high density (10$^{5}$$-$10$^{6}$ cm$^{-3}$), consistent with molecular line observations of the OMC core. The radial velocity gradients and the 3D distributions of the fingers together support the hypothesis of simultaneous, radial explosion of the Orion KL outflow.

To take full advantage of the unprecedented power of upcoming weak lensing surveys, understanding the noise, such as cosmic variance and geometry/mask effects, is as important as understanding the signal itself. Accurately quantifying the noise requires a large number of statistically independent mocks for a variety of cosmologies. This is impractical for weak lensing simulations, which are costly for simultaneous requirements of large box size (to cover a significant fraction of the past light cone) and high resolution (to robustly probe the small scale where most lensing signal resides). Therefore fast mock generation methods are desired and are under intensive investigation. We propose a new fast weak lensing map generation method, named the inverse-Gaussianization method, based on the finding that a lensing convergence field can be Gaussianized to excellent accuracy by a local transformation [Yu et al, Phys. Rev. D 84, 023523 (2011)]. Given a simulation, it enables us to produce as many as infinite statistically independent lensing maps as fast as producing the simulation initial conditions. The proposed method is tested against simulations for each tomography bin centered at lens redshift $z \sim 0.5$, 1, and 2, with various statistics. We find that the lensing maps generated by our method have reasonably accurate power spectra, bispectra, and power spectrum covariance matrix. Therefore, it will be useful for weak lensing surveys to generate realistic mocks. As an example of application, we measure the probability distribution function of the lensing power spectrum, from 16384 lensing maps produced by the inverse-Gaussianization method.

To take advantage of the astrophysical potential of Gamma-Ray Bursts (GRBs), Chinese and French astrophysicists have engaged the SVOM mission (Space-based multi-band astronomical Variable Objects Monitor). Major advances in GRB studies resulting from the synergy between space and ground observations, the SVOM mission implements space and ground instrumentation. The scientific objectives of the mission put a special emphasis on two categories of GRBs: very distant GRBs at z$>$5 which constitute exceptional cosmological probes, and faint/soft nearby GRBs which allow probing the nature of the progenitors and the physics at work in the explosion. These goals have a major impact on the design of the mission: the on-board hard X-ray imager is sensitive down to 4 keV and computes on line image and rate triggers, and the follow-up telescopes on the ground are sensitive in the NIR. At the beginning of the next decade, SVOM will be the main provider of GRB positions and spectral parameters on very short time scale. The SVOM instruments will operate simultaneously with a wide range of powerful astronomical devices. This rare instrumental conjunction, combined with the relevance of the scientific topics connected with GRB studies, warrants a remarkable scientific return for SVOM. In addition, the SVOM instrumentation, primarily designed for GRB studies, composes a unique multi-wavelength observatory with rapid slew capability that will find multiple applications for the whole astronomy community beyond the specific objectives linked to GRBs. This report lists the scientific themes that will benefit from observations made with SVOM, whether they are specific GRB topics, or more generally all the issues that can take advantage of the multi-wavelength capabilities of SVOM.

We present a new technique to measure neutrino masses using their flow field relative to dark matter. Present day streaming motions of neutrinos relative to dark matter and baryons are several hundred km/s, comparable with their thermal velocity dispersion. This results in a unique dipole anisotropic distortion of the matter-neutrino cross power spectrum, which is observable through the dipole distortion in the cross correlation of different galaxy populations. Such a dipole vanishes if not for this relative velocity and so it is a clean signature for neutrino mass. We estimate the size of this effect and find that current and future galaxy surveys may be sensitive to these signature distortions.

The shallow decay phase and flares in the afterglows of gamma-ray bursts (GRBs) is widely believed to be associated with the later activation of central engine. Some models of energy injection involve with a continuous energy flow since the GRB trigger time, such as the magnetic dipole radiation from a magnetar. However, in the scenario involving with a black hole accretion system, the energy flow from the fall-back accretion may be delayed for a fall-back time $\sim t_{\rm fb}$. Thus we propose a delayed energy injection model, the delayed energy would cause a notable rise to the Lorentz factor of the external shock, which will "generate" a bump in the multiple band afterglows. If the delayed time is very short, our model degenerates to the previous models. Our model can well explain the significant re-brightening in the optical and infrared light curves of GRB 081029 and GRB 100621A. A considerable fall-back mass is needed to provide the later energy, this indicates GRBs accompanied with fall-back material may be associated with a low energy supernova so that fraction of the envelope can be survived during eruption. The fall-back time can give meaningful information of the properties of GRB progenitor stars.

An understanding of the post-impact dynamics of ejecta clouds are crucial to the planning of a kinetic impact mission to an asteroid, and also has great implications for the history of planetary formation. The purpose of this article to track the evolution of ejecta produced by AIDA mission, which targets for kinetic impact the secondary of near-Earth binary asteroid 65803 Didymos on 2022, and to feedback essential informations to AIDA's ongoing phase-A study. We present a detailed dynamic model for the simulation of an ejecta cloud from a binary asteroid that synthesizes all relevant forces based on a previous analysis of the mechanical environment. We apply our method to gain insight into the expected response of Didymos to the AIDA impact, including the subsequent evolution of debris and dust. The crater scaling relations from laboratory experiments are employed to approximate the distributions of ejecta mass and launching speed. The size composition of fragments is modeled with a power law fitted from observations of real asteroid surface. A full-scale demonstration is simulated using parameters specified by the mission. We report the results of the simulation, which include the computed spread of the ejecta cloud and the recorded history of ejecta accretion and escape. The violent period of the ejecta evolution is found to be short, and is followed by a stage where the remaining ejecta is gradually cleared. Solar radiation pressure proves to be efficient in cleaning dust-size ejecta, and the simulation results after two weeks shows that large debris on polar orbits (perpendicular to the binary orbital plane) has a survival advantage over smaller ejecta and ejecta that keep to low latitudes.

We define a set of procedures to numerically study the fate of ejecta produced by the impact of an artificial projectile with the aim of deflecting an asteroid. Here we develop a simplified, idealized model of impact conditions that can be adapted to fit the details of specific deflection-test scenarios, such as what is being proposed for the AIDA project. Ongoing studies based upon the methodology described here can be used to inform observational strategies and safety conditions for an observing spacecraft. To account for ejecta evolution, the numerical strategies we are employing are varied and include a large N-Body component, a smoothed-particle hydrodynamics (SPH) component, and an application of impactor scaling laws. Simulations that use SPH-derived initial conditions show high-speed ejecta escaping at low angles of inclination, and very slowly moving ejecta lofting off the surface at higher inclination angles, some of which re-impacts the small-body surface. We are currently investigating the realism of this and other models' behaviors. Next steps will include the addition of solar perturbations to the model and applying the protocol developed here directly to specific potential mission concepts such as the proposed AIDA scenario.

A supra-massive neutron star (NS) spinning extremely rapidly could survive from a merger of NS-NS binary. The spin-down of this remnant NS that is highly magnetized would power the isotropic merger ejecta to produce a bright mergernova emission in ultraviolet/optical bands. Before the mergernova, the early interaction between the NS wind and the ejecta can drive a forward shock propagating outwards into the ejecta. As a result, a remarkable amount of heat can be accumulated behind the shock front and the final escaping of this heat can produce a shock breakout emission. We describe the dynamics and thermal emission of this shock with a semi-analytical model. It is found that sharp and luminous breakout emission appears mainly in soft X-rays with a luminosity of $\sim10^{45}~\rm erg~s^{-1}$ at a few hours after the merger, by leading the mergernova emission as a precursor. Therefore, detection of such an X-ray precursor would provide a smoking-gun evidence for identifying NS-powered mergernovae and distinguishing them from the radioactive-powered ones (i.e., kilonovae or macronovae). The discovery of NS-powered mergernovae would finally help to confirm the gravitational wave signals due to the mergers and the existence of supra-massive NSs.

PSR B0919+06 generally radiates radio pulses in a normal phase range. It has been known for its occasional perplexing abnormal emission events wherein individual pulses come to an earlier phase range for a few tens of periods and then returns to its usual phase. Heretofore, only a few such events have been available for study. We observed PSR B0919+06 for about 30 hours using the Jiamusi 66-m telescope at Jiamusi Deep Space Station at S-band, and detected 92 abnormal emission events. We identify four types of events based on the abrupted or gradual phase-shifting of individual pulses. The abnormal emission events are seen to occur randomly some every 1000 to 3000 periods, and they affect the leading edge of the mean profile by up to 2\% in amplitude. The abnormal emission events are probably related to gradual changes of emission processing in the pulsar magnetosphere.

We present the results of high-resolution near-IR spectroscopy toward the multiple outflows around the Herbig Be star Lk{\Ha} 234 using the Immersion Grating Infrared Spectrograph (IGRINS). Previous studies indicate that the region around Lk{\Ha} 234 is complex, with several embedded YSOs and the outflows associated with them. In simultaneous H$-$ and K$-$band spectra from HH 167, we detected 5 {\FeII} and 14 H$_{2}$ emission lines. We revealed a new {\FeII} jet driven by radio continuum source VLA 3B. Position-velocity diagrams of H$_{2}$ 1$-$0 S(1) $\lambda$2.122 $\micron$ line show multiple velocity peaks. The kinematics may be explained by a geometrical bow shock model. We detected a component of H$_{2}$ emission at the systemic velocity (V$_{LSR}$ $=$ $-$10.2 {\kms}) along the whole slit in all slit positions, which may arise from the ambient photodissociation region. Low-velocity gas dominates the molecular hydrogen emission from knots A and B in HH 167, which is close to the systemic velocity, {\FeII} emission lines are detected at farther from the systemic velocity, at V$_{LSR}$ $=$ $-$100 $-$ $-$130 {\kms}. We infer that the H$_{2}$ emission arises from shocked gas entrained by a high-velocity outflow. Population diagrams of H$_{2}$ lines imply that the gas is thermalized at a temperature of 2,500 $-$ 3,000 K and the emission results from shock excitation.

In this paper we show that the most luminous supernova discovered very recently, ASASSN-15lh, could have been powered by a newborn ultra-strongly-magnetized pulsar, which initially rotates near the Kepler limit. We find that if this pulsar is a neutron star, its rotational energy could be quickly lost as a result of gravitational-radiation-driven r-mode instability; if it is a strange quark star, however, this instability is highly suppressed due to a large bulk viscosity associated with the nonleptonic weak interaction among quarks and thus most of its rotational energy could be extracted to drive ASASSN-15lh. Therefore, we conclude that such an ultra-energetic supernova provides a possible signature for the birth of a strange quark star.

[abridged] Volume-weighted statistics of large scale peculiar velocity is preferred by peculiar velocity cosmology, since it is free of uncertainties of galaxy density bias entangled in mass-weighted statistics. However, measuring the volume-weighted velocity statistics from galaxy (halo/simulation particle) velocity data is challenging. For the first time, we apply the Kriging interpolation to obtain the volume-weighted velocity field. Kriging is a minimum variance estimator. It predicts the most likely velocity for each place based on the velocity at other places. We test the performance of Kriging quantified by the E-mode velocity power spectrum from simulations. Dependences on the variogram prior used in Kriging, the number $n_k$ of the nearby particles to interpolate and the density $n_P$ of the observed sample are investigated. First, we find that Kriging induces $1\%$ and $3\%$ systematics at $k\sim 0.1h{\rm Mpc}^{-1}$ when $n_P\sim 6\times 10^{-2} ({\rm Mpc}/h)^{-3}$ and $n_P\sim 6\times 10^{-3} ({\rm Mpc}/h)^{-3}$, respectively. The deviation increases for decreasing $n_P$ and increasing $k$. When $n_P\lesssim 6\times 10^{-4} ({\rm Mpc}/h)^{-3}$, a smoothing effect dominates small scales, causing significant underestimation of the velocity power spectrum. Second, increasing $n_k$ helps to recover small scale power. However, for $n_P\lesssim 6\times 10^{-4} ({\rm Mpc}/h)^{-3}$ cases, the recovery is limited. Finally, Kriging is more sensitive to the variogram prior for lower sample density. The most straightforward application of Kriging on the cosmic velocity field does not show obvious advantages over the nearest-particle method (Zheng et al. 2013) and could not be directly applied to cosmology so far. However, whether potential improvements may be achieved by more delicate versions of Kriging is worth further investigation.

Newly born magnetars are promising sources for gravitational wave (GW) detection due to their ultra-strong magnetic fields and high spin frequencies. Within the scenario of a growing tilt angle between the star's spin and magnetic axis, due to the effect of internal viscosity, we obtain improved estimates of the stochastic gravitational wave backgrounds (SGWBs) from magnetic deformation of newly born magnetars. We find that the GW background spectra contributed by the magnetars with ultra-strong toroidal magnetic fields of 10^{17} G could roughly be divided into four segments. Most notably, in contrast to the background spectra calculated by assuming constant tilt angles \chi=\pi/2, the background radiation above 1000 Hz are seriously suppressed. However, the background radiation at the frequency band \sim100-1000 Hz are moderately enhanced, depending on the strengths of the dipole magnetic fields. We suggest that if all newly born magnetars indeed have toroidal magnetic fields of 10^{17} G, the produced SGWBs should show sharp variations with the observed frequency at several tens to about 100 hertz. If these features could be observed through sophisticated detection of the SGWB using the proposed Einstein Telescope, it will provide us a direct evidence of the tilt angle evolutions and further some deep understandings about the properties of newly born magnetars.

The order and chaos of the motion near equilibrium points in the potential of a rotating highly irregular-shaped celestial body are investigated from point of view of the dynamical system theory. The positions of the non-degenerate equilibrium points vary continuously when the parameter changes. The topological structures in the vicinity of equilibrium points are classified into several different cases. Bifurcations at equilibrium points and the topological transfers between different cases for equilibrium points are also discussed. The conclusions can be applied to all kinds of rotating celestial bodies, simple-shaped or highly irregular-shaped, including asteroids, comets, planets and satellites of planets to help one to understand the dynamical behaviors around them. Applications to asteroids 216 Kleopatra, 2063 Bacchus, and 25143 Itokawa are significant and interesting: eigenvalues affiliated to the equilibrium points for the asteroid 216 Kleopatra move and always belong to the same topological cases; while eigenvalues affiliated to two different equilibrium points for the asteroid 2063 Bacchus and 25143 Itokawa move through the resonant cases of equilibrium points, and the collision of eigenvalues in the complex plane occurs. Poincar\'e sections in the potential of the asteroid 216 Kleopatra show the chaos behaviors of the orbits in large scale.

We provide a general analysis on the properties of emitting material of some rapidly evolving and luminous transients discovered recently with the Pan-STARRS1 Medium Deep Survey. It is found that these transients are probably produced by a low-mass non-relativistic outflow that is continuously powered by a newly born, rapidly spinning, and highly magnetized neutron star. Such a system could originate from an accretion-induced collapse of a white dwarf or a merger of a neutron star-neutron star binary. Therefore, observations to these transients would be helpful for constraining white dwarf and neutron star physics and/or for searching and identifying gravitational wave signals from the mergers.

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.

For the sake of complete theoretical research of atmospheric refraction, the atmospheric refraction under the condition of lower angles of elevation is still worthy to be analyzed and explored. In some engineering applications, the objects with larger zenith distance must be observed sometimes. Carrying out observational research of the atmospheric refraction at lower angles of elevation has an important significance. It has been considered difficult to measure the atmospheric refraction at lower angles of elevation. A new idea for determining atmospheric refraction by utilizing differential measurement with double fields of view is proposed. Taking the observational principle of HIPPARCOS satellite as a reference, a schematic prototype with double fields of view was developed. In August of 2013, experimental observations were carried out and the atmospheric refractions at lower angles of elevation can be obtained by the schematic prototype. The measured value of the atmospheric refraction at the zenith distance of 78.8 degree is $240.23"\pm0.27"$, and the feasibility of differential measurement of atmospheric refraction with double fields of view was justified. The limitations of the schematic prototype such as inadequate ability of gathering light, lack of accurate meteorological data recording and lower automatic level of observation and data processing were also pointed out, which need to be improved in subsequent work.

Cosmic shear can only be measured where there are galaxies. This source-lens clustering (SLC) effect has two sources, intrinsic source clustering and cosmic magnification (magnification/size bias). Lensing tomography can suppress the former. However, this reduction is limited by the existence of photo-z error and nonzero redshift bin width. Furthermore, SLC induced by cosmic magnification cannot be reduced by lensing tomography. Through N-body simulations, we quantify the impact of SLC on the lensing power spectrum in the context of lensing tomography. We consider both the standard estimator and the pixel-based estimator. We find that none of them can satisfactorily handle both sources of SLC. (1) For the standard estimator, SLC induced by both sources can bias the lensing power spectrum by O(1)-O(10)%. Intrinsic source clustering also increases statistical uncertainties in the measured lensing power spectrum. However, the standard estimator suppresses intrinsic source clustering in the cross-spectrum. (2) In contrast, the pixel-based estimator suppresses SLC through cosmic magnification. However, it fails to suppress SLC through intrinsic source clustering and the measured lensing power spectrum can be biased low by O(1)-O(10)%. In short, for typical photo-z errors (sigma_z/(1+z)=0.05) and photo-z bin sizes (Delta_z^P=0.2), SLC alters the lensing E-mode power spectrum by 1-10%, with ell~10^3$ and z_s~1 being of particular interest to weak lensing cosmology. Therefore the SLC is a severe systematic for cosmology in Stage-IV lensing surveys. We present useful scaling relations to self-calibrate the SLC effect.

In recent years, more and more gamma-ray bursts with late rebrightenings in multi-band afterglows unveil the late-time activities of the central engines. GRB 100814A is a special one among the well-sampled events, with complex temporal and spectral evolution. The single power-law shallow decay index of the optical light curve observed by GROND between 640 s and 10 ks is $\alpha_{\rm opt} = 0.57 \pm 0.02$, which apparently conflicts with the simple external shock model expectation. Especially, there is a remarkable rebrightening in the optical to near infrared bands at late time, challenging the external shock model with synchrotron emission coming from the interaction of the blast wave with the surrounding interstellar medium. In this paper, we invoke a magnetar with spin evolution to explain the complex multi-band afterglow emission of GRB 100814A. The initial shallow decay phase in optical bands and the plateau in X-ray can be explained as due to energy injection from a spin-down magnetar. At late time, with the falling of materials from the fall-back disk onto the central object of the burster, angular momentum of the accreted materials is transferred to the magnetar, which leads to a spin-up process. As a result, the magnetic dipole radiation luminosity will increase, resulting in the significant rebrightening of the optical afterglow. It is shown that the observed multi-band afterglow emission can be well reproduced by the model.

Long term observations by Brook et al. reveal that the derivative of rotational frequency of PSR J0738-4042 changed abruptly in 2005. Originally, the spin-down rate was relatively stable, with the rotational frequency derivative of $-1.14 \times 10^{-14}~\rm s^{-2}$. After September 2005, the derivative began to rise up. About 1000 days later, it arrived at another relatively stable value of about $-0.98 \times 10^{-14}~\rm s^{-2}$, indicating that the pulsar is spinning-down relatively slowly. To explain the observed spin-down rate change, we resort to an asteroid disrupted by PSR J0738-4042. In our model, the orbital angular momentum of the asteroid is assumed to be parallel to that of the rotating pulsar, so that the pronounced reduction in the spin-down rate can be naturally explained as due to the transfer of the angular momentum from the disrupted material to the central pulsar. The derived magnetospheric radius is about $4.0 \times 10^{9}$ cm, which is smaller than the tidal disruption radius ($4.9 \times 10^{10}$ cm). Our model is self-consistent. It is shown that the variability of the spin-down rate of PSR J0738-4042 can be quantitatively accounted for by the accretion from the asteroid disrupted by the central pulsar.

This paper studies the distribution of characteristic multipliers, the structure of submanifolds, the phase diagram, bifurcations and chaotic motions in the potential field of rotating highly irregular-shaped celestial bodies (hereafter called irregular bodies). The topological structure of the submanifolds for the orbits in the potential field of an irregular body is shown to be classified into 34 different cases, including 6 ordinary cases, 3 collisional cases, 3 degenerate real saddle cases, 7 periodic cases, 7 period-doubling cases, 1 periodic and collisional case, 1 periodic and degenerate real saddle case, 1 period-doubling and collisional case, 1 period-doubling and degenerate real saddle case, and 4 periodic and period-doubling cases. The different distribution of the characteristic multipliers has been shown to fix the structure of the submanifolds, the type of orbits, the dynamical behaviour and the phase diagram of the motion. Classifications and properties for each case are presented. Moreover, tangent bifurcations, period-doubling bifurcations, Neimark-Sacker bifurcations and the real saddle bifurcations of periodic orbits in the potential field of an irregular body are discovered. Submanifolds appear to be Mobius strips and Klein bottles when the period-doubling bifurcation occurs.

Absolute proper motions for $\sim$ 7.7 million objects were derived based on data from the South Galactic Cap u-band Sky Survey (SCUSS) and astrometric data derived from uncompressed Digitized Sky Surveys that the Space Telescope Science Institute (STScI) created from the Palomar and UK Schmidt survey plates. We put a great deal of effort into correcting the position-, magnitude-, and color-dependent systematic errors in the derived absolute proper motions. The spectroscopically confirmed quasars were used to test the internal systematic and random error of the proper motions. The systematic errors of the overall proper motions in the SCUSS catalog are estimated as -0.08 and -0.06 mas/yr for {\mu}{\alpha} cos {\delta} and {\mu}{\delta}, respectively. The random errors of the proper motions in the SCUSS catalog are estimated independently as 4.2 and 4.4 mas/yr for {\mu}{\alpha} cos {\delta} and {\mu}{\delta}. There are no obvious position-, magnitude-, and color-dependent systematic errors of the SCUSS proper motions. The random error of the proper motions goes up with the magnitude from about 3 mas/yr at u < 18.0 mag to about 7 mas/yr at u = 22.0 mag. The proper motions of stars in SCUSS catalog are compared with those in the SDSS catalog, and they are highly consistent.

Intense flares that occur at late times relative to the prompt phase have been observed by the $Swift$ satellite in the X-ray afterglows of gamma-ray bursts (GRBs). Here, we present a detailed analysis on the fall back accretion process to explain the intense flare phase in the very early X-ray afterglow light curves. To reproduce the afterglow at late times, we resort to the external shock by engaging energy injections. By applying our model to GRBs 080810, 081028 and 091029, we show that their X-ray afterglow light curves can be reproduced well. We then apply our model to the ultra-long $Swift$ GRB 111209A, which is the longest burst ever observed. The very early X-ray afterglow of GRB 111209A showed many interesting features, such as a significant bump observed at around 2000 s after the $Swift$/BAT trigger. We assume two constant energy injection processes in our model. These can explain the observed plateau at X-ray wavelength in the relatively early stage ($8.0\times10^{3}$ s) and a second X-ray plateau and optical rebrightening at about $10^{5}$ s. Our analysis supports the scenario that a significant amount of material may fall back toward the central engine after the prompt phase, causing an enhanced and long lived mass accretion rate powering a Poynting-flux-dominated outflow.

By employing a simple semi-analytical star formation model where the formation rates of Population (Pop) I/II and III stars can be calculated, respectively, we account for the number distribution of gamma-ray bursts (GRBs) with high pseudo-redshifts that was derived from an empirical luminosity-indictor relationship. It is suggested that a considerable number of Pop III GRBs could exist in the present sample of Swift GRBs. By further combining the implication for the star formation history from the optical depth of the CMB photons, it is also suggested that only a very small fraction 0.6% of Pop III GRBs could have triggered the Swift BAT. These results could provide an useful basis for estimating future detectability of Pop III stars and their produced transient phenomena.

There is accumulating evidence that at least a fraction of binary neutron star mergers result in rapidly spinning magnetars, with subrelativistic neutron-rich ejecta as massive as a small fraction of solar mass. The ejecta could be heated continuously by the Poynting flux emanated from the central magnetars. Such Poynting flux could become lepton-dominated so that a reverse shock develops. It was demonstrated that such a picture is capable of accounting for the optical transient PTF11agg (Wang & Dai 2013b). In this paper we investigate the X-ray and ultraviolet (UV) radiation as well as the optical and radio radiation studied by Wang & Dai (2013b). UV emission is particularly important because it has the right energy to ionize the hot ejecta at times $t\lesssim 600$ s. It is thought that the ejecta of binary neutron star mergers are a remarkably pure sample of r-process material, about which our understanding is still incomplete. In this paper we evaluate the possibility of observationally determining the bound-bound and bound-free opacities of the r-process material by timing the X-ray, UV, and optical radiation. It is found that these timings depend on the opacities weakly and therefore only loose constraints on the opacities can be obtained.

Some fast radio bursts (FRBs) are expected to be associated with the afterglow emission of gamma-ray bursts (GRBs), while a short-lived, supermassive neutron star (NS) forms during the GRBs. I investigate the possible contributions to the dispersion measure (DM) of the FRBs from the GRB ejecta and the wind blown from the precollapsing NS. On the one hand, sometimes an internal X-ray plateau afterglow could be produced by the NS wind, which indicates that a great number of electron-positron pairs are carried by the wind. If the pair-generation radius satisfies a somewhat rigorous condition, the relativistic and dense wind would contribute a high DM to the associated FRB, which can be comparable to and even exceed the DM contributed by the intergalactic medium. On the other hand, if the wind only carries a Goldreich-Julian particle flux, its DM contribution would become negligible; meanwhile, the internal plateau afterglow would not appear. Alternatively, the FRB should be associated with a GRB afterglow produced by the GRB external shock, i.e., an energy-injection-caused shallow-decay afterglow or a normal single-power-law afterglow if the impulsive energy release of the GRB is high enough. In the latter case, the DM contributed by the high-mass GRB ejecta could be substantially important, in particular, for an environment of main-sequence stellar wind. In summary, a careful assessment on the various DM contributors could be required for the cosmological application of the expected FRB-GRB association. The future DM measurements of GRB-associated FRBs could provide a constraint on the physics of NS winds.

Highly beamed, short-duration electromagnetic bursts could be produced by superconducting cosmic string (SCS) loops oscillating in cosmic magnetic fields. We demonstrated that the basic characteristics of SCS bursts such as the electromagnetic frequency and the energy release could be consistently exhibited in the recently discovered fast radio bursts (FRBs). Moreover, it is first showed that the redshift distribution of the FRBs can also be well accounted for by the SCS burst model. Such agreements between the FRBs and SCS bursts suggest that the FRBs could originate from SCS bursts and thus they could provide an effective probe to study SCSs. The obtained values of model parameters indicate that the loops generating the FRBs have a small length scale and they are mostly formed in the radiation-dominated cosmological epoch.

Our knowledge of the strengths of small bodies in the Solar System is limited by our poor understanding of their internal structures, and this, in turn, clouds our understanding of the formation and evolution of these bodies. Observations of the rotational states of asteroids whose diameters are larger than a few hundreds of meters have revealed that they are dominated by gravity and that most are unlikely to be monoliths; however, there is a wide range of plausible internal structures. Numerical and analytical studies of shape and spin limits of gravitational aggregates and their collisional evolution show a strong dependence on shear strength. In order to study this effect, we carry out a systematic exploration of the dependence of collision outcomes on dissipation and friction parameters of the material components making up the bodies. We simulate the catastrophic disruption (leading to the largest remnant retaining 50% of the original mass) of km-size asteroids modeled as gravitational aggregates using pkdgrav, a cosmology N-body code adapted to collisional problems and recently enhanced with a new soft-sphere collision algorithm that includes more realistic contact forces. We find that for a range of three different materials, higher friction and dissipation values increase the catastrophic disruption threshold by about half a magnitude. Furthermore, we find that pre-impact rotation systematically increases mass loss on average, regardless of the target's internal configuration. Our results have important implications for the efficiency of planet formation via planetesimal growth, and also more generally to estimate the impact energy threshold for catastrophic disruption, as this generally has only been evaluated for non-spinning bodies without detailed consideration of material properties.

Asteroid (99942) Apophis' close approach in 2029 will be one of the most significant small-body encounter events in the near future and offers a good opportunity for in situ exploration to determine the asteroid's surface properties and measure any tidal effects that might alter its regolith configuration. Resurfacing mechanics has become a new focus for asteroid researchers due to its important implications for interpreting surface observations, including space weathering effects. This paper provides a prediction for the tidal effects during the 2029 encounter, with an emphasis on whether surface refreshing due to regolith movement will occur. The potential shape modification of the object due to the tidal encounter is first confirmed to be negligibly small with systematic simulations, thus only the external perturbations are taken into account for this work (despite this, seismic shaking induced by shifting blocks might still play a weak role and we will look into this mechanism in future work). A two-stage approach is developed to model the responses of asteroid surface particles (the regolith) based on the soft-sphere implementation of the parallel N-body gravity tree code pkdgrav. A full-body model of Apophis is sent past the Earth on the predicted trajectory to generate the data of all forces acting at a target point on the surface. A sandpile constructed in the local frame is then used to approximate the regolith materials; all the forces the sandpile feels during the encounter are imposed as external perturbations to mimic the regolith's behavior in the full scenario. The local mechanical environment on the asteroid surface is represented in detail, leading to an estimation of the change in global surface environment due to the encounter. Typical patterns of perturbation are presented that depend on the asteroid orientation and sense of rotation at perigee.

In a newly born (high-temperature and Keplerian rotating) neutron star, r-mode instability can lead to stellar differential rotation, which winds the seed poloidal magnetic field ($\sim 10^{11}$ G) to generate an ultra-high ($\sim 10^{17}$ G) toroidal field component. Subsequently, by succumbing to the Tayler instability, the toroidal field could be partially transformed into a new poloidal field. Through such dynamo processes, the newly born neutron star with sufficiently rapid rotation could become a magnetar on a timescale of $\sim 10^{2-3}$ s, with a surface dipolar magnetic field of $\sim10^{15}$ G. Accompanying the field amplification, the star could spin down to a period of $\sim5$ ms through gravitational wave radiation due to the r-mode instability and, in particular, the non-axisymmetric stellar deformation caused by the toroidal field. This scenario provides a possible explanation for why the remnant neutron stars formed in gamma-ray bursts and superluminous supernovae could be millisecond magnetars.

Observations have shown that the spatial distribution of satellite galaxies is not random, but aligned with the major axes of central galaxies. This alignment is dependent on galaxy properties, such that red satellites are more strongly aligned than blue satellites. Theoretical work done to interpret this phenomena has found that it is due to the non-spherical nature of dark matter halos. However, most studies over-predict the alignment signal under the assumption that the central galaxy shape follows the shape of the host halo. It is also not clear whether the color dependence of alignment is due to an assembly bias or an evolution effect. In this paper we study these problems using a cosmological N-body simulation. Subhalos are used to trace the positions of satellite galaxies. It is found that the shape of dark matter halos are mis-aligned at different radii. If the central galaxy shares the same shape as the inner host halo, then the alignment effect is weaker and agrees with observational data. However, it predicts almost no dependence of alignment on the color of satellite galaxies, though the late accreted subhalos show stronger alignment with the outer layer of the host halo than their early accreted counterparts. We find that this is due to the limitation of pure N-body simulations that satellites galaxies without associated subhalos ('orphan galaxies') are not resolved. These orphan (mostly red) satellites often reside in the inner region of host halos and should follow the shape of the host halo in the inner region.

Laser frequency stabilization is a critical part of the interferometry measurement system of space-based gravitational wave observatories such as the Laser Interferometer Space Antenna (LISA). Arm locking as a proposed frequency stabilization technique, transfers the stability of the long arm lengths to the laser frequency. The arm locking sensor synthesizes an adequately filtered linear combination of the inter-spacecraft phase measurements to estimate the laser frequency noise, which can be used to control the laser frequency. At the University of Florida we developed the hardware-based University of Florida LISA Interferometer Simulator (UFLIS) to study and verify laser frequency noise reduction and suppression techniques under realistic LISA-like conditions. These conditions include the variable Doppler shifts between the spacecraft, LISA-like signal travel times, optical transponders, realistic laser frequency and timing noise. We review the different types of arm locking sensors and discuss their expected performance in LISA. The presented results are supported by results obtained during experimental studies of arm locking under relevant LISA-like conditions. We measured the noise suppression as well as initial transients and frequency pulling in the presence of Doppler frequency errors. This work has demonstrated the validity and \ feasibility of arm locking in LISA.

We show that the peculiar early optical and in particular X-ray afterglow emission of the short duration burst GRB 130603B can be explained by continuous energy injection into the blastwave from a supra-massive magnetar central engine. The observed energetics and temporal/spectral properties of the late infrared bump (i.e., the "kilonova") are also found consistent with emission from the ejecta launched during an NS-NS merger and powered by a magnetar central engine. The isotropic-equivalent kinetic energies of both the GRB blastwave and the kilonova are about $E_{\rm k}\sim 10^{51}$ erg, consistent with being powered by a near-isotropic magnetar wind. However, this relatively small value demands that most of the initial rotational energy of the magnetar $(\sim {\rm a~ few \times 10^{52}~ erg})$ is carried away by gravitational wave radiation. Our results suggest that (i) the progenitor of GRB 130603B would be a NS-NS binary system, whose merger product would be a supra-massive neutron star that lasted for about $\sim 1000$ seconds; (ii) the equation-of-state of nuclear matter would be stiff enough to allow survival of a long-lived supra-massive neutron star, so that it is promising to detect bright electromagnetic counterparts of gravitational wave triggers without short GRB associations in the upcoming Advanced LIGO/Virgo era.

A massive millisecond magnetar may survive a merger of a neutron star (NS) binary, which would continuously power the merger ejecta. We develop a generic dynamic model for the merger ejecta with energy injection from the central magnetar. The ejecta emission (the "merger-nova") powered by the magnetar peaks in the UV band and the peak of lightcurve progressively shifts to an earlier epoch with increasing frequency. A magnetar-powered mergernova could have an optical peak brightness comparable to a supernova, which is a few tens or hundreds times brighter than the radioactive-powered merger-novae (the so-called macro-nova or kilo-nova). On the other hand, such a merger-nova would peak earlier and have a significantly shorter duration than that of a supernova. An early collapse of the magnetar could suppress the brightness of the optical emission and shorten its duration. Such millisecond-magnetar-powered merger-novae may be detected from NS-NS merger events without an observed short gamma-ray burst, and could be a bright electromagnetic counterpart for gravitational wave bursts due to NS-NS mergers. If detected, it suggests that the merger leaves behind a massive NS, which has important implications for the equation-of-state of nuclear matter.

The determination of luminosity function (LF) of gamma-ray bursts (GRBs) is of an important role for the cosmological applications of the GRBs, which is however hindered seriously by some selection effects due to redshift measurements. In order to avoid these selection effects, we suggest to calculate pseudo-redshifts for Swift GRBs according to the empirical L-E_p relationship. Here, such a $L-E_p$ relationship is determined by reconciling the distributions of pseudo- and real redshifts of redshift-known GRBs. The values of E_p taken from Butler's GRB catalog are estimated with Bayesian statistics rather than observed. Using the GRB sample with pseudo-redshifts of a relatively large number, we fit the redshift-resolved luminosity distributions of the GRBs with a broken-power-law LF. The fitting results suggest that the LF could evolve with redshift by a redshift-dependent break luminosity, e.g., L_b=1.2\times10^{51}(1+z)^2\rm erg s^{-1}. The low- and high-luminosity indices are constrained to 0.8 and 2.0, respectively. It is found that the proportional coefficient between GRB event rate and star formation rate should correspondingly decrease with increasing redshifts.

The existence of galaxy intrinsic clustering severely hampers the weak lensing reconstruction from cosmic magnification. In paper I \citep{Yang2011}, we proposed a minimal variance estimator to overcome this problem. By utilizing the different dependences of cosmic magnification and galaxy intrinsic clustering on galaxy flux, we demonstrated that the otherwise overwhelming galaxy intrinsic clustering can be significantly suppressed such that lensing maps can be reconstructed with promising accuracy. This procedure relies heavily on the accuracy of determining the galaxy bias from the same data. Paper I adopts an iterative approach, which degrades toward high redshift. The current paper presents an alternative method, improving over paper I. We prove that the measured galaxy clustering between flux bins allows for simultaneous determination of the lensing power spectrum and the flux dependence of galaxy bias, at this redshift bin. Comparing to paper I, the new approach is not only more straightforward, but also more robust. It identifies an ambiguity in determining the galaxy bias and further discovers a mathematically robust way to suppress this ambiguity to non-negligible level ($\sim 0.1%$). The accurately determined galaxy bias can then be applied to the minimal variance estimator proposed in paper I to improve the lensing map-making. The gain at high redshift is significant. These maps can be used to measure other statistics, such as cluster finding and peak statistics. Furthermore, by including galaxy clustering measurement between different redshift bins, we can also determine the lensing cross power spectrum between these bins, up to a small and correctable multiplicative factor.

The last century has seen enormous progress in our understanding of the Universe. We know the life cycles of stars, the structure of galaxies, the remnants of the big bang, and have a general understanding of how the Universe evolved. We have come remarkably far using electromagnetic radiation as our tool for observing the Universe. However, gravity is the engine behind many of the processes in the Universe, and much of its action is dark. Opening a gravitational window on the Universe will let us go further than any alternative. Gravity has its own messenger: Gravitational waves, ripples in the fabric of spacetime. They travel essentially undisturbed and let us peer deep into the formation of the first seed black holes, exploring redshifts as large as z ~ 20, prior to the epoch of cosmic re-ionisation. Exquisite and unprecedented measurements of black hole masses and spins will make it possible to trace the history of black holes across all stages of galaxy evolution, and at the same time constrain any deviation from the Kerr metric of General Relativity. eLISA will be the first ever mission to study the entire Universe with gravitational waves. eLISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using gravitational waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the early processes at TeV energies, has guaranteed sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales around singularities and black holes, all the way to cosmological dimensions.

By phenomenologically describing the high-redshift star formation history, i.e., $\dot{\rho}_{*}(z)\propto[(1+z)/4.5]^{-\alpha}$, and semi-analytically calculating the fractions of high-redshift Pop I/II and Pop III stars, we investigate the contributions from both high-redshfit Pop I/II and Pop III stars to the observed near-infrared ($3 \mu\rm m<\lambda<5 \mu m$) excess in the cosmic infrared background emission. In order to account for the observational level of the near-infrared excess, the power-law index $\alpha$ of the assumed star formation history is constrained to within the range of $0\lesssim\alpha\lesssim1$. Such a constraint is obtained under the condition that the virial temperature of dark matter halos belongs to the range of $500 {\rm K}\leq T_{\rm vir}\leq10^4$ K.

With assistance of the identified/associated sources in the second Fermi Large Area Telescope (LAT) catalog, we analyze and resolve the spatial distribution and the distributions of the gamma-ray spectral and variability indices of the remaining 575 unassociated Fermi LAT sources. Consequently, it is suggested that the unassociated sources could statistically consist of Galactic supernova remnants/pulsar wind nebulae, BL Lacertae objects, flat spectrum radio quasars, and other types of active galaxies with fractions of 25%, 29%, 41%, and 5%, respectively.

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 thermal evolution of strange stars in both normal and color-flavor-locked (CFL) phases are investigated together with the evolutions of the stellar rotation and the r-mode instability. The heating effects due to the deconfinement transition of the stellar crust and the dissipation of the r-modes are considered. As a result, the cooling of the stars in the normal phase is found to be not very different from the standard one. In contrast, for the stars in the CFL phase, a big bump during the first hundred years and a steep decay ($\sim$7% in ten years) at the ages of $\sim10^{4-6}$ yrs are predicted in their thermal evolution curves. These unique features provide an effective observational test for determining whether or not the CFL phase is reached in strange stars. This thermal test method is independent of and complementary to the rotational test method, which is a direct consequence of the r-mode instability.

The afterglow of GRB 081029 showed unusual behavior, with a significant rebrightening being observed at optical wavelength at about 3000 s after the burst. One possible explanation is that the rebrightening is resulted from energy injection. Here, we present a detailed numerical study of the energy injection process and interpret the X-ray and optical afterglow light curves of GRB 081029. In our model, we have assumed two periods of energy injection, each with a constant injection power. One injection starts at $2.8\times10^{3}$ s and lasts for about 2500 s, with a power of $7.0\times10^{47}$ ${\rm erg}$ ${\rm s^{-1}}$. This energy injection is mainly engaged to account for the rapid rebrightening at about 3000 s. The other injection starts at $8.0\times10^{3}$ s and lasts for about 5000 s. The injection power is $3.5\times10^{47}$ ${\rm erg}$ ${\rm s^{-1}}$. This energy injection can help to explain the slight rebrightening at about 10000 s. It is shown that the observed optical afterglow, especially the marked rebrightening at about 3000 s, can be well reproduced. In X-ray band, the predicted amplitude of the rebrightening is much shallower, which is also consistent with the observed X-ray afterglow light curve. It is argued that the two periods of energy injection can be produced by the falling of clumpy materials onto the central compact object of the burster, which leads to an enhancement of accretion and gives birth to a strong outflow temporarily.

Sunspot numbers form a comprehensive, long-duration proxy of solar activity and have been used numerous times to empirically investigate the properties of the solar cycle. A number of correlations have been discovered over the 24 cycles for which observational records are available. Here we carry out a sophisticated statistical analysis of the sunspot record that reaffirms these correlations, and sets up an empirical predictive framework for future cycles. An advantage of our approach is that it allows for rigorous assessment of both the statistical significance of various cycle features and the uncertainty associated with predictions. We summarize the data into three sequential relations that estimate the amplitude, duration, and time of rise to maximum for any cycle, given the values from the previous cycle. We find that there is no indication of a persistence in predictive power beyond one cycle, and conclude that the dynamo does not retain memory beyond one cycle. Based on sunspot records up to October 2011, we obtain, for Cycle 24, an estimated maximum smoothed monthly sunspot number of 97 +- 15, to occur in January--February 2014 +- 6 months.

With a simple power-law approximation of high-redshift ($\gtrsim3.5$) star formation history, i.e., $\dot{\rho}_*(z)\propto [(1+z)/4.5]^{-\alpha}$, we investigate the reionization of intergalactic medium (IGM) and the consequent Thomson scattering optical depth for cosmic microwave background (CMB) photons. A constraint on the evolution index $\alpha$ is derived from the CMB optical depth measured by the {\it Wilkinson Microwave Anisotropy Probe} (WMAP) experiment, which reads $\alpha\approx2.18\lg{\mathscr{N}_{\gamma}}-3.89$, where the free parameter $\mathscr{N}_\gamma$ is the number of the escaped ionizing ultraviolet photons per baryon. Moreover, the redshift for full reionization, $z_f$, can also be expressed as a function of $\alpha$ as well as $\mathscr{N}_{\gamma}$. By further taking into account the implication of the Gunn-Peterson trough observations to quasars for the full reionization redshift, i.e., $6\lesssim z_f \lesssim7$, we obtain $0.3\lesssim\alpha\lesssim1.3$ and $80\lesssim\mathscr{N}_{\gamma}\lesssim230$. For a typical number of $\sim4000$ of ionizing photons released per baryon of normal stars, the fraction of these photons escaping from the stars, $f_{\rm esc}$, can be constrained to within the range of $(2.0-5.8)%$.

In paper I (Yu et al. [1]), we show through N-body simulation that a local monotonic Gaussian transformation can significantly reduce non-Gaussianity in a noise-free lensing convergence field. This makes the Gaussianization a promising theoretical tool to understand high-order lensing statistics. Here we present a study of its applicability in lensing data analysis, in particular when shape measurement noise is presented in lensing convergence maps. (i) We find that shape measurement noise significantly degrades the Gaussianization performance and the degradation increases for shallower surveys. (ii) The Wiener filter is efficient in reducing the impact of shape measurement noise. The Gaussianization of the Wiener-filtered lensing maps is able to suppress skewness, kurtosis, and the 5th- and 6th-order cumulants by a factor of 10 or more. It also works efficiently to reduce the bispectrum to zero.

Motivated by recent works of Neyrinck et al. 2009 and Scherrer et al. 2010, we proposed a Gaussianization transform to Gaussianize the non-Gaussian lensing convergence field $\kappa$. It performs a local monotonic transformation $\kappa\rightarrow y$ pixel by pixel to make the unsmoothed one-point probability distribution function of the new variable $y$ Gaussian. We tested whether the whole $y$ field is Gaussian against N-body simulations. (1) We found that the proposed Gaussianization suppresses the non-Gaussianity by orders of magnitude, in measures of the skewness, the kurtosis, the 5th- and 6th-order cumulants of the $y$ field smoothed over various angular scales relative to that of the corresponding smoothed $\kappa$ field. The residual non-Gaussianities are often consistent with zero within the statistical errors. (2) The Gaussianization significantly suppresses the bispectrum. Furthermore, the residual scatters around zero, depending on the configuration in the Fourier space. (3) The Gaussianization works with even better performance for the 2D fields of the matter density projected over $\sim 300 \mpch$ distance interval centered at $z\in(0,2)$, which can be reconstructed from the weak lensing tomography. (4) We identified imperfectness and complexities of the proposed Gaussianization. We noticed weak residual non-Gaussianity in the $y$ field. We verified the widely used logarithmic transformation as a good approximation to the Gaussianization transformation. However, we also found noticeable deviations.

The GRB outflow driven by superconducting cosmic strings is likely to be an arc rather than a usually-considered spherical cap. In such a case, the afterglows of the cosmic string GRBs could be basically consistent with the observation of the high-redshift GRBs.

The accumulation of {\it Swift} observed gamma-ray bursts (GRBs) gradually makes it possible to directly derive a GRB luminosity function (LF) from observational luminosity distribution, where however two complexities must be involved as (i) the evolving connection between GRB rate and cosmic star formation rate and (ii) observational selection effects due to telescope thresholds and redshift measurements. With a phenomenological investigation on these two complexities, we constrain and discriminate two popular competitive LF models (i.e., broke-power-law LF and single-power-law LF with an exponential cutoff at low luminosities). As a result, we find that the broken-power-law LF could be more favored by the observation, with a break luminosity $L_b=2.5\times10^{52}\rm erg s^{-1}$ and prior- and post-break indices $\nu_1=1.72$ and $\nu_2=1.98$. For an extra evolution effect expressed by a factor $(1+z)^{\delta}$, if the matallicity of GRB progenitors is lower than $\sim0.1Z_{\odot}$ as expected by some collapsar models, then there may be no extra evolution effect other than the metallicity evolution (i.e., $\delta$ approaches to be zero). Alternatively, if we remove the theoretical metallicity requirement, then a relationship between the degenerate parameters $\delta$ and $Z_{\max}$ can be found, very roughly, $\delta\sim2.4(Z_{\max}/Z_{\odot}-0.06)$. This indicates that an extra evolution could become necessary for relatively high metallicities.

PSR B1259-63/SS 2883 is a binary system in which a 48-ms pulsar orbits around a Be star in a high eccentric orbit with a long orbital period of about 3.4 yr. Extensive broadband observational data are available for this system from radio band to very high energy (VHE) range. The multi-frequency emission is unpulsed and nonthermal, and is generally thought to be related to the relativistic electrons accelerated from the interaction between the pulsar wind and the stellar wind, where X-ray emission is from the synchrotron process and the VHE emission is from the inverse Compton (IC) scattering process. Here a shocked wind model with variation of the magnetic parameter $\sigma$ is developed for explaining the observations. By choosing proper param- eters, our model could reproduce two-peak profile in X-ray and TeV light curves. The effect of the disk exhibits an emission and an absorption components in the X-ray and TeV bands respectively. We suggest that some GeV flares will be produced by Doppler boosting the synchrotron spectrum. This model can possibly be used and be checked in other similar systems such as LS I+61o303 and LS 5039.

The fluorescent X-ray emission from neutral iron in the molecular clouds (Sgr B) indicates that the clouds are being irradiated by an external X-ray source. The source is probably associated with the Galactic central black hole (Sgr A*), which triggered a bright outburst one hundred years ago. We suggest that such an outburst could be due to a partial capture of a star by Sgr A*, during which a jet was generated. By constraining the observed flux and the time variability ($\sim$ 10 years) of the Sgr B's fluorescent emission, we find that the shock produced by the interaction of the jet with the dense interstellar medium represents a plausible candidate for the X-ray source emission.

Properties of proton resonances in $^{18}$Ne have been investigated efficiently by utilizing a technique of proton resonant elastic scattering with a $^{17}$F radioactive ion (RI) beam and a thick proton target. A 4.22~MeV/nucleon $^{17}$F RI beam was produced via a projectile-fragmentation reaction, and subsequently separated by a Radioactive Ion Beam Line in Lanzhou ({\tt RIBLL}). Energy spectra of the recoiled protons were measured by two sets of $\Delta$E-E silicon telescope at center-of-mass scattering angles of $\theta_{c.m.}$$\approx$175${^\circ}$$\pm$5${^\circ}$, $\theta_{c.m.}$$\approx$152${^\circ}$$\pm$8${^\circ}$, respectively. Several proton resonances in $^{18}$Ne were observed, and their resonant parameters have been determined by an $R$-matrix analysis of the differential cross sections in combination with the previous results. The resonant parameters are related to the reaction-rate calculation of the stellar $^{14}$O($\alpha$,$p$)$^{17}$F reaction, which was thought to be the breakout reaction from the hot CNO cycles into the $rp$-process in x-ray bursters. Here, $J^\pi$=(3$^-$, 2$^-$) are tentatively assigned to the 6.15-MeV state which was thought the key 1$^-$ state previously. In addition, a doublet structure at 7.05 MeV are tentatively identified, and its contribution to the resonant reaction rate of $^{14}$O($\alpha$,$p$)$^{17}$F could be enhanced by at least factors of about 4$\sim$6 in comparison with the previous estimation involving only a singlet. The present calculated resonant rates are much larger than those previous values, and it may imply that this breakout reaction could play a crucial role under x-ray bursters conditions.

Context: Swift observations suggest that the X-ray afterglow emission of some gamma-ray bursts (GRB) may have internal origins, and the conventional external shock (ES) cannot be the exclusive source of the afterglow emission. Aims: If the central compact objects of some GRBs are millisecond magentars, the magnetar winds could play an important role in the (internal) X-ray afterglow emission, which is our focus here. Methods: The dynamics and the synchrotron radiation of the termination shock (TS) of the magmnetar winds, as well as the simultaneous GRB ES, are investigated by considering the magnetization of the winds. Results: As a result of the competition between the emission of the wind TS and the GRB ES, two basic types of X-ray afterglows are predicted, i.e., the TS-dominated and the ES-dominated types. Moreover, our results also show that both of the two types of afterglows have a shallow-decay phase and a normal-decay one, as observed by the \textit{Swift} satellite. This indicates that some observed X-ray afterglows could be (internally) produced by the magnetar winds, but not necessarily GRB ESs.

The high-redshift gamma-ray bursts (GRBs), GRBs 080913 and 090423, challenge the conventional GRB progenitor models by their short durations, typical for short GRBs, and their high energy releases, typical for long GRBs. Meanwhile, the GRB rate inferred from high-redshift GRBs also remarkably exceeds the prediction of the collapsar model, with an ordinary star formation history. We show that all these contradictions could be eliminated naturally, if we ascribe some high-redshift GRBs to electromagnetic bursts of superconducting cosmic strings. High-redshift GRBs could become a reasonable way to test the superconducting cosmic string model, because the event rate of cosmic string bursts increases rapidly with increasing redshifts, whereas the collapsar rate decreases.

Swift observations suggest that the central compact objects of some gamma-ray bursts (GRBs) could be newly born millisecond magnetars. Therefore, by considering the spin evolution of the magnetars against r-mode instability, we investigate the role of the magnetars in GRB X-ray afterglow emission. Besides modifying the conventional energy injection model, we pay particular attention to the internal X-ray afterglow emission, whose luminosity is assumed to track the magnetic dipole luminosity of the magentars with a certain fraction. Following a comparison between the model and some selected observational samples, we suggest that some so-called "canonical" X-ray afterglows including the shallow decay, normal decay, and steeper-than-normal decay phases could be internally produced by the magnetars (possibly through some internal dissipations of the magnetar winds), while the (energized) external shocks are associated with another type of X-ray afterglows. If this is true, from those internal X-ray afterglows, we can further determine the magnetic field strengths and the initial spin periods of the corresponding magnetars.

The thermal evolution of neutron stars (NSs) is investigated by coupling with the evolution of $\textit{r}$-mode instability that is described by a second order model.The heating effect due to shear viscous damping of the $\textit{r}$-modes enables us to understand the high temperature of two young pulsars (i.e., PSR B0531+21 and RX J0822-4300) in the framework of the simple $npe$ NS model, without superfluidity or exotic particles.Moreover, the light curves predicted by the model within an acceptable parameter regime may probably cover all of the young and middle-aged pulsars in the $\lg T_s^{\infty}-\lg t$ panel, and an artificially strong $p$ superfluidity invoked in some early works is not needed here. Additionally, by considering the radiative viscous damping of the $\textit{r}$-modes, a surprising extra cooling effect is found, which can even exceed the heating effect sometimes although plays an ignorable role in the thermal history.

The observations of gamma-ray bursts (GRBs) such as 980425, 031203 and 060218, with luminosities much lower than those of other classic bursts, lead to the definition of a new class of GRBs -- low-luminosity GRBs. The nature of the outflow responsible for them is not clear yet. Two scenarios have been suggested: one is the conventional relativistic outflow with initial Lorentz factor of order of $\Gamma_0\ga 10$ and the other is a trans-relativistic outflow with $\Gamma_0\simeq 1-2$. Here we compare the high energy gamma-ray afterglow emission from these two different models, taking into account both synchrotron self inverse-Compton scattering (SSC) and the external inverse-Compton scattering due to photons from the cooling supernova or hypernova envelope (SNIC). We find that the conventional relativistic outflow model predicts a relatively high gamma-ray flux from SSC at early times ($<10^4 {\rm s}$ for typical parameters) with a rapidly decaying light curve, while in the trans-relativistic outflow model, one would expect a much flatter light curve of high-energy gamma-ray emission at early times, which could be dominated by both the SSC emission and SNIC emission, depending on the properties of the underlying supernova and the shock parameter $\epsilon_e$ and $\epsilon_B$. The Fermi Gamma-ray Space Telescope should be able to distinguish between the two models in the future.

The internal-plateau X-ray emission of gamma-ray bursts (GRBs) indicates that a newly born magnetar could be the central object of some GRBs. The observed luminosity and duration of the plateaus suggest that, for such a magnetar, a rapid spin with a sub- or millisecond period is sometimes able to last thousands of seconds. In this case, the conventional neutron star (NS) model for the magnetar may be challenged, since the rapid spin of nascent NSs would be remarkably decelerated within hundreds of seconds due to r-mode instability. In contrast, the r-modes can be effectively suppressed in nascent strange stars (SSs). In other words, to a certain extent, only SSs can keep nearly-constant extremely-rapid spin for a long period of time during the early ages of the stars. We thus propose that the sample of the GRB rapidly-spinning magnetars can be used to test the SS hypothesis based on the distinct spin limits of NSs and SSs.

In a second-order r-mode theory, S'a & Tom'e found that the r-mode oscillation in neutron stars (NSs) could induce stellar differential rotation, which leads to a saturation state of the oscillation spontaneously. Based on a consideration of the coupling of the r-modes and the stellar spin and thermal evolutions, we carefully investigate the influences of the r-mode-induced differential rotation on the long-term evolutions of isolated NSs and NSs in low-mass X-ray binaries, where the viscous damping of the r-modes and its resultant effects are taken into account. The numerical results show that, for both kinds of NSs, the differential rotation can prolong the duration of the r-mode saturation state significantly. As a result, the stars can keep nearly constant temperature and angular velocity over a thousand years. Moreover, due to the long-term steady rotation of the stars, persistent quasi-monochromatic gravitational wave radiation could be expected, which increases the detectibility of gravitational waves from both nascent and accreting old NSs.

X-rays from planetary nebulae (PNs) are believed to originate from a shock driven into the fast stellar wind (v ~ 1000 km/s) as it collides with an earlier circumstellar slow wind (v ~ 10 km/s). In theory, the shocked fast wind (hot bubble) and the ambient cold nebula can remain separated by magnetic fields along a surface referred to as the contact discontinuity (CD) that inhibits diffusion and heat conduction. The CD region is extremely difficult to probe directly owing to its small size and faint emission. This has largely left the study of CDs, stellar-shocks, and the associated micro-physics in the realm of theory. This paper presents spectroscopic evidence for ions from the hot bubble (kT ~ 100 eV) crossing the CD and penetrating the cold nebular gas (kT ~ 1 eV). Specifically, a narrow radiative recombination continuum (RRC) emission feature is identified in the high resolution X-ray spectrum of the PN BD+30 3639 indicating bare C VII ions are recombining with cool electrons at kT_e=1.7+-1.3 eV. An upper limit to the flux of the narrow RRC of H-like C VI is obtained as well. The RRCs are interpreted as due to C ions from the hot bubble of BD+30 3639 crossing the CD into the cold nebula, where they ultimately recombine with its cool electrons. The RRC flux ratio of C VII to C VI constrains the temperature jump across the CD to Delta kT > 80 eV, providing for the first time direct evidence for the stark temperature disparity between the two sides of an astrophysical CD, and constraining the role of magnetic fields and heat conduction accordingly. Two colliding-wind binaries are noted to have similar RRCs suggesting a temperature jump and CD crossing by ions may be a common feature of stellar wind shocks.

In the popular internal shock model for the prompt emission of gamma-ray bursts (GRBs), collisions between a series of relativistic shells generate lots of paired forward and reverse shocks. We show that the synchrotron emission produced by the forward and reverse shocks respectively could peak at two quite different energy bands if the Lorentz factors of these two types of shocks are significantly different with each other (e.g., one shock is relativistic and the other is Newtonian). We then investigate whether this scenario is applicable to the case of GRB 080319B and find that a bimodal distribution of the shell Lorentz factors, peaking at $\sim400$ and $\sim10^5$, is required. In addition, this scenario predicts an accompanying inverse-Compton (IC) GeV emission with a luminosity comparable to (not much higher than) that of the synchrotron MeV emission, which can be tested with future \textit{Fermi} observations.

We study afterglow flares of gamma-ray bursts (GRBs) in the framework of the late internal shock (LIS) model based on a careful description for the dynamics of a pair of shocks generated by a collision between two homogeneous shells,. First, by confronting the model with some fundamental observational features of X-ray flares, we find some constraints on the properties of the pre-collision shells that are directly determined by the central engine of GRBs. Second, high energy emission associated with X-ray flares, which arises from synchrotron self-Compton (SSC) emission of LISs, is investigated in a wide parameter space. The predicted flux of high energy flares may reach as high as $\sim 10^{-8}\rm erg cm^{-2}s^{-1}$, which is likely to be detectable with the Large Area Telescope (LAT) aboard \textit{the Fermi Space Telescope}

We present the results of the first X-ray gratings spectroscopy observations of a planetary nebula (PN), the X-ray-bright, young BD+30 3639. We observed BD+30 3639 for a total of 300 ks with the Chandra X-ray Observatory's Low Energy Transmission Gratings in combination with its Advanced CCD Imaging Spectrometer(LETG/ACIS-S). The LETG/ACIS-S spectrum of BD+30 3639 is dominated by H-like resonance lines of O viii and C sc vi and the He-like triplet line complexes of Ne ix and O vii. Other H-like resonance lines, such as N vii, as well as lines of highly ionized Fe, are weak or absent. Continuum emission is evident over the range 6-18 A. Spectral modeling indicates the presence of a range of plasma temperatures from T~1.7x10^6 K to 2.9x10^6 K and an intervening absorbing column N_H~2.4x10^21 cm-2. The same modeling conclusively demonstrates that C and Ne are highly enhanced, with abundance ratios of C/O~15-45 and Ne/O~3.3-5.0 (90% confidence ranges, relative to the solar ratios), while N and Fe are depleted, N/O~0.0-1.0 and Fe/O~0.1-0.4. The intrinsic luminosity of the X-ray source determined from the modeling and the measured flux (F_X = 4.1x10^-13 ergs cm-2 s-1) is L_X~8.6x10^32 erg s-1(assuming D = 1.2kpc). These gratings spectroscopy results are generally consistent with earlier results obtained from X-ray CCD imaging spectroscopy of BD+30 3639, but are far more precise. The tight constraints placed on the (nonsolar) abundances directly implicate the present-day central star -- hence, ultimately, the intershell region of the progenitor asymptotic giant branch star -- as the origin of the shocked plasma now emitting in X-rays.

The observations of a nearby low-luminosity gamma-ray burst (GRB) 060218 associated with supernova SN 2006aj may imply an interesting astronomical picture where a supernova shock breakout locates behind a relativistic GRB jet. Based on this picture, we study neutrino emission for early afterglows of GRB 060218-like GRBs, where neutrinos are expected to be produced from photopion interactions in a GRB blast wave that propagates into a dense wind. Relativistic protons for the interactions are accelerated by an external shock, while target photons are basically provided by the incoming thermal emission from the shock breakout and its inverse-Compton scattered component. Because of a high estimated event rate of low-luminosity GRBs, we would have more opportunities to detect afterglow neutrinos from a single nearby GRB event of this type by IceCube. Such a possible detection could provide evidence for the picture described above.

The widely existing shallow decay phase of the X-ray afterglows of gamma-ray bursts (GRBs) is generally accepted to be due to long-lasting energy injection. The outflows carrying the injecting energy, based on the component that is dominative in energy, fall into two possible types: baryon-dominated and lepton-dominated ones. The former type of outflow could be ejecta that is ejected during the prompt phase of a GRB and consists of a series of baryonic shells with a distribution of Lorentz factors, and the latter type could be an electron-positron-pair wind that is driven by the post-burst central engine. We here provide a unified description for the dynamics of fireballs based on these two types of energy injection, and calculate the corresponding high-energy photon emission by considering synchrotron radiation and inverse Compton scattering (including synchrotron self-Compton and combined inverse-Compton) of electrons. We find that, in the two energy-injection models, there is a plateau (even a hump) in high-energy light curves during the X-ray shallow decay phase. In particular, a considerable fraction of the injecting energy in the lepton-dominated model can be shared by the long-lasting reverse shock since it is relativistic. Furthermore, almost all of the energy of the reverse shock is carried by leptons, and thus the inverse-Compton emission is enhanced dramatically. Therefore, this model predicts more significant high-energy afterglow emission than the baryon-dominated model. We argue that these observational signatures would be used to discriminate between different energy-injection models in the upcoming {\em Gamma-ray Large Area Space Telescope} (GLAST) era.

The postburst object of a GRB is likely to be a highly magnetized, rapidly rotating compact object (e.g., a millisecond magnetar), which could produce an ultrarelativistic electron-positron-pair wind. The interaction of such a wind with an outwardly expanding fireball ejected during the burst leads to a relativistic wind bubble (RWB). We numerically calculate the dynamics and radiative properties of RWBs and use this model to explain the shallow decay phase of the early X-ray afterglows observed by Swift. We find that RWBs can fall into two types: forward-shock-dominated and reverse-shock-dominated bubbles. Their radiation during a period of $\sim 10^{2}-10^{5}$ seconds is dominated by the shocked medium and the shocked wind, respectively, based on different magnetic energy fractions of the shocked materials. For both types, the resulting light curves always have a shallow decay phase. In addition, we provide an example fit to the X-ray afterglows of GRB 060813 and GRB 060814 and show that they could be produced by forward-shock-dominated and reverse-shock-dominated bubbles, respectively. This implies that, for some early afterglows (e.g., GRB 060814), the long-lasting reverse shock emission is strong enough to explain their shallow decay phase.

We present preliminary results from the first X-ray gratings spectrometer observations of a planetary nebula (PN). We have used the Chandra X-ray Observatory Low Energy Transmission Gratings Spectrometer (LETGS) to observe the bright, diffuse X-ray source within the well-studied BD +30 3639. The LETGS spectrum of BD +30 3639 displays prominent and well-resolved emission lines of H-like C, O, and Ne and He-like O and Ne. Initial modeling indicates a plasma temperature T_x ~ 2.5x10^6 K and abundance ratios of C/O ~ 20, N/O <~ 1, Ne/O ~ 4, and Fe/O <~ 0.1. These results suggest that the X-ray-emitting plasma is dominated by the shocked fast wind from the emerging PN core, where this wind gas likely originated from the intershell region of the progenitor asymptotic giant branch star.

Under the combination effect of the recommencement heating due to spin-down of strange stars and the heat perseveration due to weak conduct heat of the crust, the Cooper pair breaking and formation(PBF) in color superconduction quark matter arises. We investigated the cooling of the strange stars with a crust in color superconductivity phase including both decomfinement heating and PBF process. We find that deconfinement heating can delay the thermal evolution of strange stars and the PBF process suppresses the early temperature rise of the stars. The cooling strange stars behave within the brightness constraint of young compact objects when the color superconductivity gap is small enough.

We studied the evolution of isolated strange stars synthetically, considering the influence of {\it r-}mode instability. Our results show that the cooling of strange stars with non-ultra strong magnetic fields is delayed by the heating due to the {\it r-}modes damping during million years, while the spin-down of the stars is dominated by gravitational radiation. Especially for the strange stars in a possible existing color-flavor locked phase, the effect of the {\it r-}mode instability on the evolution of the stars becomes extremely important since the viscosity, the neutrino emissivity, and the specific heat involving paring quarks are blocked. It leads to the cooling of these color superconducting stars is very slow, and the stars can remain high temperature within million years differing completely from previous understanding. In this case, a strange star in color-flavor locked phase can be located at the bottom of its {\it r-}mode instability window for a long time, but does not spin down to a very low frequency within hours.

We study in a fully selfconsistent approach the structure of a vortex in low density superfluid neutron matter. We determine that the matter density profile of a vortex shows a significant depletion in the region of the core, a feature never reported for a vortex state in a Fermi superfluid.