There are many stars that are rotating spheroids in the Universe, and studying them is of very important significance. Since the times of Newton, many astronomers and physicists have researched gravitational properties of stars by considering the moment equations derived from Eulerian hydrodynamic equations. In this paper we study the scattering of spinors of the Dirac equation, and in particular investigate the scattering issue in the limit case of rotating Maclaurin spheroids. Firstly we give the metric of a rotating ellipsoid star, then write the Dirac equation under this metric, and finally derive the scattering solution to the Dirac equation and establish a relation between differential scattering cross-section, $\sigma$, and stellar matter density, $\mu$. It is found that the sensitivity of $\sigma$ to the change in $\mu$ is proportional to the density $\mu$. Because of weak gravitational field and constant mass density, our results are reasonable. The results can be applied to white dwarfs, main sequence stars, red giants, supergiant stars and so on, as long as their gravitational fields are so weak that they can be treated in the Newtonan approximations, and the fluid is assumed to be incompressible. Notice that we take the star's matter density to be its average density and the star is not taken to be compact. Obviously our results cannot be used to study neutron stars and black holes. In particular, our results are suitable for white dwarfs, which have average densities of about $10^{5}-10^{6}$\,g~cm$^{-3}$, corresponding to a range of mass of about $0.21-0.61 M_{\bigodot}$ and a range of radius of about $6000-10000$\,km.

PSR J1640+2224 is a binary millisecond pulsar (BMSP) with a white dwarf (WD) companion. Recent observations indicate that the WD is very likely to be a $\sim 0.7\,M_{\odot}$ CO WD. Thus the BMSP should have evolved from an intermediate-mass X-ray binary (IMXB). However, previous investigations on IMXB evolution predict that the orbital periods of the resultant BMSPs are generally $<40$ days, in contrast with the 175 day orbital period of PSR J1640+2224. In this paper, we explore the influence of the mass of the neutron star (NS) and the chemical compositions of the companion star on the formation of BMSPs. Our results show that, the final orbital period becomes longer with increasing NS mass, and the WD mass becomes larger with decreasing metallicity. In particular, to reproduce the properties of PSR J1640+2224, the NS was likely born massive ($>2.0\,M_{\odot}$).

According to the recycling model, neutron stars in low-mass X-ray binaries were spun up to millisecond pulsars (MSPs), which indicates that all MSPs in the Galactic plane ought to be harbored in binaries. However, about $20\%$ Galactic field MSPs are found to be solitary. To interpret this problem, we assume that the accreting neutron star in binaries may collapse and become a strange star when it reaches some critical mass limit. Mass loss and a weak kick induced by asymmetric collapse during the phase transition (PT) from neutron star to strange star can result in isolated MSPs. In this work, we use a population-synthesis code to examine the PT model. The simulated results show that a kick velocity of $\sim60~{\rm km~s}^{-1}$ can produce $\sim6\times10^3$ isolated MSPs and birth rate of $\sim6.6\times10^{-7} {\rm ~yr}^{-1}$ in the Galaxy, which is approximately in agreement with predictions from observations. For the purpose of comparisons with future observation, we also give the mass distributions of radio and X-ray binary MSPs, along with the delay time distribution.

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 low braking-index pulsar PSR J1734$-$3333 could be born with superhigh internal magnetic fields $B_{\rm in}\sim10^{15}-10^{16}$ G, and undergo a supercritical accretion soon after its formation in a supernova explosion. The buried multipole magnetic fields will merger into a dipole magnetic field. Since the magnetic flow transfers from the core to the crust of the pulsar, its surface dipole field grows quickly at a power-law form assumed until it saturates at the level of internal dipole field. The increase in surface dipole magnetic field results in the observed low braking index of $n=0.9(2)$. Keeping an average field growth index $\varepsilon=1.34(6)$, this pulsar will become a magnetar with surface dipole magnetic field at the equator $B_{\rm d}\sim 2.6(1)\times 10^{14}$\,G and $\sim 5.3(2)\times 10^{14}$\,G after next 50\,kyrs and 100\,kys, respectively.

In this work, we reinvestigate the electron fraction $Y_{e}$ and electron Fermi energy $E_{F}(e)$ of neutron stars, based on our previous work of Li et al.(2016), in which we firstly deduced a special solution to $E_{F}(e)$, and then obtained several useful analytical formulae for $Y_{\rm e}$ and matter density $\rho$ within classical models and the relativistic mean field(RMF) theory using numerically fitting. The advantages of this work include the following aspects:(1) The linear functions are substituted for the nonlinear exponential functions used in the previous work. This method may be more simple, and closer to realistic equation of state\,(EoS) of a neutron star(NS), because there are linear or quasi-linear relationships between number fractions of leptons and matter density, which can be seen by solving NS EoS; (2)we introduce a dimensionless variable $\varrho$\,($\varrho=\rho/\rho_0$, $\rho_{0}$ is the standard saturated nuclear density), which greatly reduces the scope of the fitting coefficients;(3)we present numerical errors including absolute and relative deviations between the data and fit. By numerically simulating, we have obtained several analytical formulae for $Y_{e}$ and $\rho$ for both APR98 and RMF models. Combining these analytical formulae with the special solution, we can calculate the value of $E_{\rm F}(e)$ for any given matter density. Since $Y_e$ and $E_{ F}(e)$ are important in assessing cooling rate of a NS and the possibility of kaon/pion condensation in the NS interior, this study could be useful in the future study on the thermal evolution of a NS.

In this work, we interpreted the high braking index of PSR J1640$-$4631 with a combination of the magneto-dipole radiation and dipole magnetic field decay models. By introducing a mean rotation energy conversion coefficient $\overline{\zeta}$, the ratio of the total high-energy photon energy to the total rotation energy loss in the whole life of the pulsar, and combining the pulsar's high-energy and timing observations with reliable nuclear equation of state, we estimate the pulsar's initial spin period, $P_{0}\sim (17-44)$ ms, corresponding to the moment of inertia $I\sim (0.8-2.1)\times 10^{45}$ g cm$^{2}$. Assuming that PSR J1640$-$4631 has experienced a long-term exponential decay of the dipole magnetic field, we calculate the true age $t_{\rm age}$, the effective magnetic field decay timescale $\tau_{D}$, and the initial surface dipole magnetic field at the pole $B_{p}(0)$ of the pulsar to be $(2900-3100)$ yrs, $1.07(2)\times10^{5}$ yrs, and $(1.84-4.20)\times10^{13}$ G, respectively. The measured braking index of $n=3.15(3)$ for PSR J1640$-$4631 is attributed to its long-term dipole magnetic field decay and a low magnetic field decay rate, $dB_{\rm p}/dt\sim -(1.66-3.85)\times10^{8}$ G yr$^{-1}$. Our model can be applied to both the high braking index ($n>3$) and low braking index ($n<3$) pulsars, tested by the future polarization, timing, and high-energy observations of PSR J1640$-$4631.

In this paper we discuss in detail the quantization of Landau energy levels of a strongly magnetized and completely degenerate relativistic electron gas in neutron stars. In particular, we focus on the Fermi energy dependence of the magnetic field for a relativistic electron gas in the superstrong magnetic field of magnetars. We would like to point out that some of the results concerning the microscopic number density of states of a strongly magnetized electron gas given by well-known statistical physics text books are incorrect. The repeated use of these results in the last five decades probably seriously affects the physics of neutron stars and magnetars. The quantization of Landau energy levels is accurately delineated in terms of the Dirac (delta) function. Relatively simple calculation shows that the Fermi energy of relativistic electron gas in magnetars with superstrong magnetic fields( B > Bcr , here Bcr is the Landau critical magnetic field) increases with magnetic fieldstrength as B**(1/4). On the basis of this simple but importantnew result we are leading naturally to propose a new mechanism for the production of high X-ray luminosity from magnetars.

eXTP is a science mission designed to study the state of matter under extreme conditions of density, gravity and magnetism. Primary targets include isolated and binary neutron stars, strong magnetic field systems like magnetars, and stellar-mass and supermassive black holes. The mission carries a unique and unprecedented suite of state-of-the-art scientific instruments enabling for the first time ever the simultaneous spectral-timing-polarimetry studies of cosmic sources in the energy range from 0.5-30 keV (and beyond). Key elements of the payload are: the Spectroscopic Focusing Array (SFA) - a set of 11 X-ray optics for a total effective area of about 0.9 m^2 and 0.6 m^2 at 2 keV and 6 keV respectively, equipped with Silicon Drift Detectors offering <180 eV spectral resolution; the Large Area Detector (LAD) - a deployable set of 640 Silicon Drift Detectors, for a total effective area of about 3.4 m^2, between 6 and 10 keV, and spectral resolution <250 eV; the Polarimetry Focusing Array (PFA) - a set of 2 X-ray telescope, for a total effective area of 250 cm^2 at 2 keV, equipped with imaging gas pixel photoelectric polarimeters; the Wide Field Monitor (WFM) - a set of 3 coded mask wide field units, equipped with position-sensitive Silicon Drift Detectors, each covering a 90 degrees x 90 degrees FoV. The eXTP international consortium includes mostly major institutions of the Chinese Academy of Sciences and Universities in China, as well as major institutions in several European countries and the United States. The predecessor of eXTP, the XTP mission concept, has been selected and funded as one of the so-called background missions in the Strategic Priority Space Science Program of the Chinese Academy of Sciences since 2011. The strong European participation has significantly enhanced the scientific capabilities of eXTP. The planned launch date of the mission is earlier than 2025.

In this paper, we investigate the electron Landau-level stability and its influence on the electron Fermi energy, $E_{\rm F}(e)$, in the circumstance of magnetars, which are powered by magnetic field energy. In a magnetar, the Landau levels of degenerate and relativistic electrons are strongly quantized. A new quantity $g_{n}$, the electron Landau-level stability coefficient is introduced. According to the requirement that $g_{n}$ decreases with increasing the magnetic field intensity $B$, the magnetic-field index $\beta$ in the expression of $E_{\rm F}(e)$ must be positive. By introducing the Dirac$-\delta$ function, we deduce a general formulae for the Fermi energy of degenerate and relativistic electrons, and obtain a particular solution to $E_{\rm F}(e)$ in a superhigh magnetic field (SMF). This solution has a low magnetic-field index of $\beta=1/6$, compared with the previous one, and works when $\rho\geq 10^{7}$~g cm$^{-3}$ and $B_{\rm cr}\ll B\leq 10^{17}$~Gauss. By modifying the phase space of relativistic electrons, a SMF can enhance the electron number density $n_e$, and decrease the maximum of electron Landau level number, which results in a redistribution of electrons. According to Pauli exclusion principle, the degenerate electrons will fill quantum states from the lowest Landau level to the highest Landau level. As $B$ increases, more and more electrons will occupy higher Landau levels, though $g_{n}$ decreases with the Landau level number $n$. The enhanced $n_{e}$ in a SMF means an increase in the electron Fermi energy and an increase in the electron degeneracy pressure. The results are expected to facilitate the study of the weak-interaction processes inside neutron stars and the magnetic-thermal evolution mechanism for megnetars.

Based on the basic definition of Fermi energy of degenerate and relativistic electrons, we obtain a special solution to electron Fermi energy, $E_{\rm F}(e)$, and express $E_{\rm F}(e)$ as a function of electron fraction, $Y_{e}$, and matter density, $\rho$. Several useful analytical formulae for $Y_{e}$ and $\rho$ within classical models and the work of Dutra et al. 2014 (Type-2) in relativistic mean field theory are obtained using numerically fitting. When describing the mean-field Lagrangian, density, we adopt the TMA parameter set, which is remarkably consistent with with the updated astrophysical observations of neutron stars. Due to the importance of the density dependence of the symmetry energy, $S$, in nuclear astrophysics, a brief discussion on the symmetry parameters $S_v$ and $L$ (the slope of $S$) is presented. Combining these fit formulae with boundary conditions for different density regions, we can evaluate the value of $E_{\rm F}(e)$ in any given matter density, and obtain a schematic diagram of $E_{\rm F}(e)$ as a continuous function of $\rho$. Compared with previous study on the electron Fermi energy in other models, our methods of calculating $E_{\rm F}(e)$ are more simple and convenient, and can be universally suitable for the relativistic electron regions in the circumstances of common neutron stars. We have deduced a general expression of $E_{\rm F}(e)$ and $n_{e}$, which could be used to indirectly test whether one EoS of a NS is correct in our future studies on neutron star matter properties. Since URCA reactions are expected in the center of a massive star due to high-value electron Fermi energy and electron fraction, this study could be useful in the future studies on the NS thermal evolution.

By introducing Dirac ?-function in superhigh magnetic field, we deduce a general formula for pressure of degenerate and relativistic electrons, Pe, which is suitable for superhigh magnetic fields, discuss the quantization of Landau levels of electrons, and consider the quantum electrodynamic(QED) effects on the equations of states (EOSs) for different matter systems. The main conclusions are as follows: the stronger the magnetic field strength, the higher the electron pressure becomes; compared with a common radio pulsar, a magnetar could be a more compact oblate spheroid-like deformed neutron star due to the anisotropic total pressure; and an increase in the maximum mass of a magnetar is expected because of the positive contribution of the magnetic field energy to the EOS of the star. Since this is an original work in which some uncertainties could exist, to further modify and perfect our theory model should be considered in our future studies.

Due to the lack of long term pulsed emission in quiescence and the strong timing noise, it is impossible to directly measure the braking index $n$ of a magnetar. Based on the estimated ages of their potentially associated supernova remnants (SNRs), we estimate the values of the mean braking indices of eight magnetars with SNRs, and find that they cluster in a range of $1\sim$42. Five magnetars have smaller mean braking indices of $1<n<3$, and we interpret them within a combination of magneto-dipole radiation and wind aided braking, while the larger mean braking indices of $n>3$ for other three magnetars are attributed to the decay of external braking torque, which might be caused by magnetic field decay. We estimate the possible wind luminosities for the magnetars with $1<n<3$, and the dipolar magnetic field decay rates for the magnetars with $n>3$ within the updated magneto-thermal evolution models. Although the constrained range of the magnetars' braking indices is tentative, due to the uncertainties in the SNR ages, which come from distance uncertainties and the unknown conditions of the expanding shells, our method provides an effective way to constrain the magnetars' braking indices if the measurements of the SNRs' ages are reliable, which can be improved by future observations.

In this paper, we consider the effect of Landau levels on the decay of superhigh magnetic fields of magnetars. Applying ${}^3P_2$ anisotropic neutron superfluid theory yield a second-order differential equation for a superhigh magnetic field $B$ and its evolutionary timescale $t$. The superhigh magnetic fields may evolve on timescales $\sim (10^{6}-10^{7})$ yrs for common magnetars. According to our model, the activity of a magnetar may originate from instability caused by the high electron Fermi energy.

In a superhigh magnetic field, direct Urca reactions can proceed for an arbitrary proton concentration. Since only the electrons with high energy $E$ ($E > Q$, $Q$ is the threshold energy of inverse $\beta-$decay) at large Landau levels can be captured, we introduce the Landau level effect coefficient $q$ and the effective electron capture rate $\Gamma_{\rm eff}$. By using $\Gamma_{\rm eff}$, the values of $L_{\rm X}$ and $L_{\rm \nu}$ are calculated, where and $L_{\rm \nu}$, $L_{\rm X}$ are the average neutrino luminosity of Anomalous X-ray Pulsars (AXPs) and the average X-ray luminosity of AXPs, respectively. The complete process of electron capture inside a magnetar is simulated numerically.

As soon as the energy of electrons near the Fermi surface are higher than $Q$, the threshold energy of inverse $\beta-$ decay, the electron capture process will dominate. The resulting high-energy neutrons will destroy anisotropic ${}^3P_2$ neutron superfluid Cooper pairs. By colliding with the neutrons produced in the process $n+ (n\uparrow n\downarrow)\longrightarrow n+ n+ n$, the kinetic energy of the outgoing neutrons will be transformed into thermal energy. The transformed thermal energy would transported from the star interior to the star surface by conduction, then would be transformed into radiation energy as soft X-rays and gamma-rays. After a highly efficient modulation within the pulsar magnetosphere, the surface thermal emission (mainly soft X/$\gamma$-ray emission) has been shaped into a spectrum with the observed characteristics of magnetars. By introducing two important parameters: Landau level-superfluid modified factor and effective X/$\gamma$-ray coefficient, we numerically simulate the process of magnetar cooling and magnetic field decay, and then compute magnetars' soft X/$\gamma$-ray luminosities $L_{X}$. Further, we obtain aschematic diagrams of $L_{X}$ as a function of magnetic field strength $B$. The observations are compared with the calculations.

In this paper, we modify our previous research carefully, and derive a new expression of electron energy density in superhigh magnetic fields. Based on our improved model, we re-compute the electron capture rates and the magnetic fields' evolutionary timescales $t$ of magnetars. According to the calculated results, the superhigh magnetic fields may evolve on timescales $\sim (10^{6}-10^{7})$ yrs for common magnetars, and the maximum timescale of the field decay, $t\approx 2.9507 \times 10^{6}$ yrs, corresponding to an initial internal magnetic field $B_{\rm 0}= 3.0 \times 10^ {15}$ G and an initial inner temperature $T_{\rm 0}= 2.6 \times 10^ {8}$ K. Motivated by the results of the neutron star-supernova remnant(SNR) association of Zhang $\&$ Xie(2011), we calculate the maximum $B_{\rm 0}$ of magnetar progenitors, $B_{\rm max}\sim (2.0\times 10^{14}-2.93 \times 10^{15})$ G when $T_{\rm 0}= 2.6 \times 10^ {8}$ K. When $T_{\rm 0}\sim 2.75 \times 10^ {8}-~1.75 \times 10^ {8}$ K, the maximum $B_{\rm 0}$ will also be in the range of $\sim 10^{14}-10^{15}$ G, not exceeding the upper limit of magnetic field of a magnetar under our magnetar model. We also investigate the relationship between the spin-down ages of magnetars and the ages of their SNRs, and explain why all AXPs associated with SNRs look older than their real ages, whereas all SGRs associated with SNRs appear younger than they are.

Based on our previous work, we deduce a general formula for pressure of degenerate and relativistic electrons,Pe, which is suitable for superhigh magnetic fields, discuss the quantization of Landau levels of electrons, and consider the quantum electrodynam-ic(QED) effects on the equations of states (EOSs) for different matter systems. The main conclusions are as follows:Pe is related to the magnetic field B, matter density ?, and electron fraction Ye ; the stronger the magnetic field, the higher the electron pressure becomes; the high electron pressure could be caused by high Fermi energy of electrons in a superhigh magnetic field; compared with a common radio pulsar, a magnetar could be a more compact oblate spheroid-like deformed neutron star due to the anisotropic total pressure; and an increase in the maximum mass of a magnetar is expected because of the positive contribution of the magnetic field energy to the EOS of the star.

Assuming that the timescale of the magnetic field decay is approximately equal to that of the stellar cooling via neutrino emission, we obtain a one-to-one relationship between the effective surface thermal temperature and the inner temperature. The ratio of the effective neutrino luminosity to the effective X-ray luminosity decreases with decaying magnetic field.

In this paper, an approximate method of calculating the Fermi energy of electrons ($E_{F}(e)$) in a high-intensity magnetic field, based on the analysis of the distribution of a neutron star magnetic field, has been proposed. In the interior of a Neutron star, different forms of intense magnetic field could exist simultaneously and a high electron Fermi energy could be generated by the release of magnetic field energy. The calculation results show that: $E_{F}(e)$ is related to density $\rho$, the mean electron number per baryon $Y_{e}$ and magnetic field strength $B$.

In order to depict the quantization of Landau levels, we introduce Dirac $\delta$ function, and gain a concise expression for the electron Fermi energy, $E_{F}(e) \propto B^{1/4}$. The high soft X-ray luminosities of magnetars may be naturally explained by our theory.

Seven years of pulse time-of-arrival measurements have been collected from observations of the young pulsar PSR B2334+61 using the Nanshan radio telescope of Urumqi Observatory. A phase-connected solution has been obtained over the whole data span, 2002 August to 2009 August. This includes a very large glitch that occurred between 2005 August 26 and September 8 (MJDs 53608 and 53621). The relative increase in rotational frequency for this glitch, $\Delta\nu_{g}/\nu~\sim~20.5\times10^{-6}$, is the largest ever seen. Although accounting for less than 1\% of the glitch, there were two well-defined exponential decay terms with time constants of 21 and 147 days respectively. There was also a large long-term increase in the spindown rate with $\Delta\dot\nu_p/\dot\nu \sim 0.011$ at the time of the glitch. A highly significant oscillation with a period of close to one year is seen in the post-glitch residuals. It is very unlikely that this can be accounted for by a pulsar position error or proper motion -- it appears to result from effects interior to the neutron star. Implications of these results for pulsar glitch models are discussed.

The triaxial-octupole Y$_{32}$ correlation in atomic nuclei has long been expected to exist but experimental evidence has not been clear. We find, in order to explain the very low-lying 2$^-$ bands in the transfermium mass region, that this exotic effect may manifest itself in superheavy elements. Favorable conditions for producing triaxial-octupole correlations are shown to be present in the deformed single-particle spectrum, which is further supported by quantitative Reflection Asymmetric Shell Model calculations. It is predicted that the strong nonaxial-octupole effect may persist up to the element 108. Our result thus represents the first concrete example of spontaneous breaking of both axial and reflection symmetries in the heaviest nuclear systems.

A method for calculation of Gamow-Teller transition rates is developed by using the concept of the Projected Shell Model (PSM). The shell model basis is constructed by superimposing angular-momentum-projected multi-quasiparticle configurations, and nuclear wave functions are obtained by digonalizing the two-body interactions in these projected states. Calculation of transition matrix elements in the PSM framework is discussed in detail, and the effects caused by the Gamow-Teller residual forces and by configuration-mixing are studied. With this method, it may become possible to perform a state-by-state calculation for beta-decay and electron-capture rates in heavy, deformed nuclei at finite temperatures. Our first example indicates that, while experimentally known Gamow-Teller transition rates from the ground state of the parent nucleus are reproduced, stronger transitions from some low-lying excited states are predicted to occur, which may considerably enhance the total decay rates once these nuclei are exposed to hot stellar environments.