A pure two-body problem has seven integrals including the Kepler energy, the Laplace vector, and the angular momentum vector. However, only five of them are independent. When the five independent integrals are preserved, the two other dependent integrals are naturally preserved from a theoretical viewpoint; but they may not be either from a numerical computational viewpoint. Because of this, we use seven scale factors to adjust the integrated positions and velocities so that the adjusted solutions strictly satisfy the seven constraints. Noticing the existence of the two dependent integrals, we adopt the Newton iterative method combined with the singular value decomposition to calculate these factors. This correction scheme can be applied to perturbed two-body and N-body problems in the solar system. In this case, the seven quantities of each planet slowly vary with time. More accurate values can be given to the seven slowly-varying quantities by integrating the integral invariant relations of these quantities and the equations of motion. They should be satisfied with the adjusted solutions. Numerical tests show that the new method can significantly reduce the rapid growth of numerical errors of all orbital elements.

While a growing body of research indicates that relativistic magnetic reconnection is a prodigious source of particle acceleration in high-energy astrophysical systems, the dominant acceleration mechanism remains controversial. Using a combination of fully kinetic simulations and theoretical analysis, we demonstrate that Fermi-type acceleration within the large-scale motional electric fields dominates over direct acceleration from non-ideal electric fields within small-scale diffusion regions. This result has profound implications for modeling particle acceleration in large-scale astrophysical problems, since it opens up the possiblity of modeling the energetic spectra without resolving microscopic diffusion regions.

Recently, Su et al. (2016) propose a telescope design, referred as the "SYZ" design, for Chinese new project of a 12m optical-infrared telescope. The SYZ telescope design consists of three aspheric mirrors with non-zero power, including a relay mirror below the primary mirror. SYZ design yields a good imaging quality and has a relatively flat field curvature at Nasmyth focus. To evaluate the science-compatibility of this three-mirror telescope, in this paper, we thoroughly compare the performance of SYZ design with that of Ritchey-Chr\'etien (RC) design, a conventional two-mirror telescope design. Further, we propose the Observing Information Throughput ($OIT$) as a metric for quantitatively evaluating the telescopes' science performance. We find that although a SYZ telescope yields a superb imaging quality over a large field of view, a two-mirror (RC) telescope design holds a higher overall throughput, a better diffraction-limited imaging quality in the central field of view (FOV$<5'$) which is better for the performance of extreme Adaptive Optics (AO), and a generally better scientific performance with a higher $OIT$ value.

N-body simulations of dark matter halos show that the density is cusped near the center of the halo. The density profile behaves as $r^{-\gamma}$ in the inner parts, where $\gamma \simeq 1$ for the NFW model and $\gamma \simeq 1.5$ for the Moore's model, but in the outer parts, both models agree with each other in the asymptotic behavior of the density profile. The simulations also show the information about anisotropy parameter $\beta(r)$ of velocity distribution. $\beta\approx 0$ in the inner part and $\beta\approx 0.5$ (radially anisotropic) in the outer part of the halo. We provide some distribution functions $F(E,L)$ with the constant anisotropy parameter $\beta$ for the two spherical models of dark matter halos: a new generalized NFW model and a generalized Moore model. There are two parameters $\alpha$ and $\epsilon$ for those two generalized models to determine the asymptotic behavior of the density profile. In this paper, we concentrate on the situation of $\beta(r)=1/2$ from the viewpoint of the simulation.

N-body simulations of dark matter halos show that the density profiles of halos behave as $\rho(r)\propto r^{-\alpha(r)}$, where the density logarithmic slope $\alpha \simeq 1\sim1.5$ in the center and $\alpha \simeq 3\sim 4$ in the outer parts of halos. However, some observations are not in agreement with simulations in the very central region of halos. The simulations also show that velocity dispersion anisotropy parameter $\beta\approx 0$ in the inner part of the halo and the so called "pseudo phase-space density" $\rho/\sigma^3$ behaves as a power-law in radius $r$. With these results in mind, we study the distribution function and the pseudo phase-space density $\rho/\sigma^3$ of the center of dark matter halos and find that they are closely-related.

In this paper we implement the new agegraphic dark energy model with quintessence field. We demonstrate that the new agegraphic evolution of the universe can be described completely by a single quintessence field. Its potential as a function of the quintessence field is reconstructed numerically. In particular, the analytical solution of the new agegraphic quintessence dark energy model (NAQDE) is approximately obtained in the matter-dominated epoch. Furthermore, we investigate the evolution of the NAQDE model in the $\omega-\omega'$ phase plane. It turns out that by quantum corrections, the trajectory of this model lies outside the thawing and freezing regions at early times. But at late times, it enters the freezing regions and gradually approaches to a static cosmological constant state in the future. Therefore the NAQDE should belong to the freezing model at late times. For comparison, we further extend this model by including the interaction between the NADE and DM and discuss its evolution in the $\omega-\omega'$ phase plane.