The Giant Radio Array for Neutrino Detection (GRAND) is a planned large-scale observatory of ultra-high-energy (UHE) cosmic particles, with energies exceeding 10^8 GeV. Its goal is to solve the long-standing mystery of the origin of UHE cosmic rays. To do this, GRAND will detect an unprecedented number of UHE cosmic rays and search for the undiscovered UHE neutrinos and gamma rays associated to them with unmatched sensitivity. GRAND will use large arrays of antennas to detect the radio emission coming from extensive air showers initiated by UHE particles in the atmosphere. Its design is modular: 20 separate, independent sub-arrays, each of 10 000 radio antennas deployed over 10 000 km^2. A staged construction plan will validate key detection techniques while achieving important science goals early. Here we present the science goals, detection strategy, preliminary design, performance goals, and construction plans for GRAND.

Using the stellar evolution code---Modules for Experiments in Stellar Astrophysics (MESA), we investigate the evolution of massive stars with different rotational velocities and metallicities towards Wolf-Rayet stars. In our simulations, the initial rotating velocities are taken as 0, 250, 500 and 650 km s$^{-1}$, and the metallicities equal to 0.02, 0.014, 0.008, 0.006, 0.004 and 0.002. We show our rapid rotation models in the HR diagram compared with the observations. We find that the rotational mixing is less efficient at high metallicity, and these stars become Wolf-Rayet (WR) stars when the helium in their center is ignited. However, rapid rotating massive stars at low metallicity can easily evolve into WR stars due to the rotation resulted in chemically homogeneous evolution. This can explain the origin of single WR stars in galaxy at low metallicity. In our models, the observed SMC WR stars are consistent with the single-star evolution models. However at higher metallicities our single-star evolution models can only explain the luminous, hydrogen-rich WN stars and O stars (which are classified as WR stars previously).

We investigate the effects of the core-collapse supernova ejecta on a rapidly rotating and massive companion star. We show that the stripped mass raises by twice when compare with a massive but non-rotating companion star. In close binaries with orbital periods of about 1 day, the stripped masses reach up to $\sim 1 M_\odot$. By simulating the evolutions of the rotational velocities of the massive companion stars based on different stripped masses, we find that the rotational velocity decreases greatly for stripped mass that is higher than about $1 M_\odot$. Of all the known high mass X-ray binaries (HMXBs), Cygnus X-3 and 1WGA J0648.024418 have the shortest orbital periods of 0.2 and 1.55 days, respectively. The optical counterpart of the former is a Wolf-Rayet star, whereas it is a hot subdwarf for the latter. Applying our model to the two HMXBs, we suggest that the hydrogen-rich envelopes of their optical counterparts may have been stripped by CCSN ejecta.

According to the nova model from \citet{Yaron2005} and \citet{Jose1998} and using Monte Carlo simulation method, we investigate the contribution of chemical abundances in nova ejecta to the interstellar medium (ISM) of the Galaxy. We find that the ejected mass by classical novae (CNe) is about $2.7\times10^{-3}$ $ \rm M_\odot\ {\rm yr^{-1}}$. In the nova ejecta, the isotopic ratios of C, N and O, that is, $^{13}$C/$^{12}$C, $^{15}$N/$^{14}$N and $^{17}$O/$^{16}$O, are higher about one order of magnitude than those in red giants. We estimate that about 10$\%$, 5$\%$ and 20$\%$ of $^{13}$C, $^{15}$N and $^{17}$O in the ISM of the Galaxy come from nova ejecta, respectively. However, the chemical abundances of C, N and O calculated by our model can not cover all of observational values. This means that there is still a long way to go for understanding novae.

Neutron stars are ideal astrophysical laboratories for testing theories of the de Haas-van Alphen (dHvA) effect and diamagnetic phase transition which is associated with magnetic domain formation. The "magnetic interaction" between delocalized magnetic moments of electrons (the Shoenberg effect), can result in an effect of the diamagnetic phase transition into domains of alternating magnetization (Condon's domains). Associated with the domain formation are prominent magnetic field oscillation and anisotropic magnetic stress which may be large enough to fracture the crust of magnetar with a super-strong field. Even if the fracture is impossible as in "low-field" magnetar, the depinning phase transition of domain wall motion driven by low field rate (mainly due to the Hall effect) in the randomly perturbed crust can result in a catastrophically variation of magnetic field. This intermittent motion, similar to the avalanche process, makes the Hall effect be dissipative. These qualitative consequences about magnetized electron gas are consistent with observations of magnetar emission, and especially the threshold critical dynamics of driven domain wall can partially overcome the difficulties of "low-field" magnetar bursts and the heating mechanism of transient, or "outbursting" magnetar.

As the evolutionary link between the radio millisecond pulsars (MSPs) and the low mass X-ray binaries or intermediate mass X-ray binaries, the millisecond X-ray pulsars (MSXPs) are important objects in testing theories of pulsar formation and evolution. In general, neutron stars in MSXPs can form via core collapse supernova (CC channel) of massive stars or accretion induced collapse (AIC channel) of an accreting ONeMg WD whose mass reaches the Chandrasekhar limit. Here, in addition to CC and AIC channels we also consider another channel, i.e., evolution induced collapse (EIC channel) of a helium star with mass between $1.4$ and $2.5 M_{\odot}$. Using a population synthesis code, we have studied MSXPs arising from three different evolutionary channels. We find that the Galactic birthrates of transient MSXPs and persistent MSXPs are about 0.7---$1.4\times 10^{-4}$ yr$^{-1}$. Our population synthesis calculations have shown that about 50\%---90\% of the MSXPs have undergone CC channel, about 10\%---40\% of them have undergone EIC channel, and the MSXPs via AIC channel are the least.

The New Vacuum Solar Telescope (NVST) is a 1 meter vacuum solar telescope that aims to observe the fine structures on the Sun. The main tasks of NVST are high resolution imaging and spectral observations, including the measurements of solar magnetic field. NVST is the primary ground-based facility of Chinese solar community in this solar cycle. It is located by the Fuxian Lake of southwest China, where the seeing is good enough to perform high resolution observations. In this paper, we first introduce the general conditions of Fuxian Solar Observatory and the primary science cases of NVST. Then, the basic structures of this telescope and instruments are described in detail. Finally, some typical high resolution data of solar photosphere and chromosphere are also shown.

The origin of dust grains in the interstellar medium is still open problem. \cite{Nicholls2013} found the presence of a significant amount of dust around V1309 Sco which maybe originate from the merger of a contact binary. We investigate the origin of dust around V1309 Sco, and suggest that these dust grains are efficiently produced in the binary-merger ejecta. By means of \emph{AGBDUST} code, we estimate that $\sim 5.2\times10^{-4} M_\odot$ of dust grains are produced, and their radii are $\sim 10^{-5}$ cm. These dust grains mainly are composed of silicate and iron grains. Because the mass of the binary-merger ejecta is very small, the contribution of dust produced by binary-merger ejecta to the overall dust production in the interstellar medium is negligible. However, it is the most important that the discovery of a significant amount of dust around V1309 Sco offers a direct support for the idea---common-envelope ejecta provides an ideal environment for dust formation and growth. Therefore, we confirm that common-envelope ejecta can be important source of cosmic dust.

The material that is ejected in a common-envelope (CE) phase in a close binary system provides an ideal environment for dust formation. By constructing a simple toy model to describe the evolution of the density and the temperature of CE ejecta and using the \emph{AGBDUST} code to model dust formation, we show that dust can form efficiently in this environment. The actual dust masses produced in the CE ejecta depend strongly on their temperature and density evolution. We estimate the total dust masses produced by CE evolution by means of a population synthesis code and show that, compared to dust production in AGB stars, the dust produced in CE ejecta may be quite significant and could even dominate under certain circumstances.

This work aims to present our current best physical understanding of common-envelope evolution (CEE). We highlight areas of consensus and disagreement, and stress ideas which should point the way forward for progress in this important but long-standing and largely unconquered problem. Unusually for CEE-related work, we mostly try to avoid relying on results from population synthesis or observations, in order to avoid potentially being misled by previous misunderstandings. As far as possible we debate all the relevant issues starting from physics alone, all the way from the evolution of the binary system immediately before CEE begins to the processes which might occur just after the ejection of the envelope. In particular, we include extensive discussion about the energy sources and sinks operating in CEE, and hence examine the foundations of the standard energy formalism. Special attention is also given to comparing the results of hydrodynamic simulations from different groups and to discussing the potential effect of initial conditions on the differences in the outcomes. We compare current numerical techniques for the problem of CEE and also whether more appropriate tools could and should be produced (including new formulations of computational hydrodynamics, and attempts to include 3D processes within 1D codes). Finally we explore new ways to link CEE with observations. We compare previous simulations of CEE to the recent outburst from V1309 Sco, and discuss to what extent post-common-envelope binaries and nebulae can provide information, e.g. from binary eccentricities, which is not currently being fully exploited.

Symbiotic X-ray binaries (SyXBs) comprise a rare class of low-mass X-ray binaries. We study the Galactic SyXBs, which we consider as detached binaries composed of low-mass giants and wind-fed neutron star companions, by simulation of the interaction of a magnetized neutron star (NS) with its environment and utilizing a population synthesis code. We focus mainly on the parameters that influence observational appearance of the SyXB: the donor wind velocity (vw) and the angular momentum distribution in the shell of matter settling onto NS. We estimate the birthrate of SyXB as $\sim 4.1\times 10^{-5}$ yr$^{-1}$ to $ \sim 6.6\times 10^{-6}$ yr$^{-1}$ and their number in the Galaxy as $\sim$(100 -- 1000). Assumed stellar wind velocity from cool giants is the input parameter that influences the model SyXBs population most. Among known SyXBs or candidate systems, 4U 1954+31 and IGR J16358-4724 in which NS have very long spin periods may host quasi-spherically accreting NSs. GX 1+4 has a peculiar long-term spin behaviour and it may also be a quasi-spherical wind-accreting source. We cannot identify whether there are wind-fed accretion disks in 4U 1700+24, Sct X-1, IRXS J180431.1-273932 and 2XMM J174016.0-290337.

Using a simple accelerating model and an assumption that $\gamma$-rays originate from $p-p$ collisions for a $\pi^0$ model, we investigate $\gamma$-ray sources like V407 Cygni in symbiotic stars. The upper limit of their occurrence rate in the Galaxy is between 0.5 and 5 yr$^{-1}$, indicating that they may be an important source of the high-energy $\gamma$-rays. The maximum energies of the accelerated protons mainly distribute around $10^{11}$ eV, and barely reach $10^{15}$ eV. The novae occurring in D-type SSs with ONe WDs and long orbital periods are good candidates for $\gamma$-ray sources. Due to a short orbital period which results in a short acceleration duration, the nova occurring in symbiotic star RS Oph can not produce the $\gamma$-ray emission like that in V407 Cygni.

Assuming that soft X-ray sources in symbiotic stars result from strong thermonuclear runaways, and supersoft X-ray sources from weak thermonuclear runaways or steady hydrogen burning symbiotic stars, we investigate the Galactic soft and supersoft X-ray sources in symbiotic stars by means of population synthesis. The Galactic occurrence rates of soft X-ray sources and supersoft X-ray sources are from $\sim$ 2 to 20 $\rm yr^{-1}$, and $\sim$ 2 to 17 $\rm yr^{-1}$, respectively. The numbers of X-ray sources in symbiotic stars range from 2390 to 6120. We simulate the distribution of X-ray sources over orbital periods, masses and mass-accretion rates of white dwarfs. The agreement with observations is reasonable.

We have carried out a study of the chemical abundances of $^1$H, $^4$He, $^{12}$C, $^{13}$C, $^{14}$N, $^{15}$N, $^{16}$O, $^{17}$O, $^{20}$Ne and $^{22}$Ne in symbiotic stars (SSs) by means of a population synthesis code. We find that the ratios of the number of O-rich SSs to that of C-rich SSs in our simulations are between 3.4 and 24.1, depending on the third dredge-up efficiency $\lambda$ and the terminal velocity of the stellar wind $v(\infty)$. The fraction of SSs with $extrinsic$ C-rich cool giants in C-rich cool giants ranges from 2.1% to 22.7%, depending on $\lambda$, the common envelope algorithm and the mass-loss rate. Compared with the observations, the distributions of the relative abundances of $^{12}$C/$^{13}$C vs. [C/H] of the cool giants in SSs suggest that the thermohaline mixing in low-mass stars may exist. The distributions of the relative abundances of C/N vs. O/N, Ne/O vs. N/O and He/H vs. N/O in the symbiotic nebulae indicate that it is quite common that the nebular chemical abundances in SSs are modified by the ejected materials from the hot components. Helium overabundance in some symbiotic nebulae may be relevant to a helium layer on the surfaces of white dwarf accretors.

We have carried out a detailed study of symbiotic stars with white dwarf accretors by means of a population synthesis code. We estimate the total number of symbiotic stars with white dwarf accretors in the Galaxy as 1,200 - 15,000. This range is compatible with observational estimates. Two crucial physical parameters that define the birthrate and number of symbiotic stars are the efficiency of accretion by white dwarfs (which greatly depends on the separation of components after common envelope stage and stellar wind velocity) and the mass of the hydrogen layer which the white dwarf can accumulate prior to the hydrogen ignition. The theoretical estimate of the Galactic occurrence rate of symbiotic novae ranges from about 1.3 to about 13.5 per year out of which weak symbiotic novae comprise about 0.5 to 6.0 per year, depending on the model assumptions. We simulate the distributions of symbiotic stars over orbital periods, masses of components, mass-loss rates of cool components, mass-accretion rates of hot components and luminosity of components. Agreement with observations is reasonable.