The spectral energy distribution (SED) sequence for type Ia supernovae (SN Ia) is modeled by an artificial neural network. The SN Ia luminosity is characterized as a function of phase, wavelength, a color parameter and a decline rate parameter. After training and testing the neural network, the SED sequence could give both the spectrum with wavelength range from 3000\AA~to 8000\AA~ and the light curve with phase from 20 days before to 50 days after the maximum luminosity for the supernovae with different colors and decline rates. Therefore, we call this the Artificial Neural Network Spectral Light Curve Template (ANNSLCT) model. We retrain the Joint Light-curve Analysis (JLA) supernova sample by using the ANNSLCT model and obtain the parameters for each supernova to make a constraint on the cosmological $\Lambda$CDM model. We find that the best fitting values of these parameters are almost the same as those from the JLA sample trained with the Spectral Adaptive Lightcurve Template 2 (SALT2) model. So we believe that the ANNSLCT model could be used to analyze a large number of SN Ia multi-color light curves measured in the current and future observational projects.

Glitches correspond to sudden jumps of rotation frequency ($\nu$) and its derivative ($\dot{\nu}$) of pulsars, the origin of which remains not well understood yet, partly because the jump processes of most glitches are not well time-resolved. There are three large glitches of the Crab pulsar, detected in 1989, 1996 and 2017, which were found to have delayed spin-up processes before the normal recovery processes. Here we report two additional glitches of the Crab pulsar occurred in 2004 and 2011 for which we discovered delayed spin up processes, and present refined parameters of the largest glitch occurred in 2017. The initial rising time of the glitch is determined as $<0.48$ hour. We also carried out a statistical study of these five glitches with observed spin-up processes. The two glitches occurred in 2004 and 2011 have delayed spin-up time scales ($\tau_{1}$) of $1.7\pm0.8$\,days and $1.6\pm0.4$\,days, respectively. We find that the $\Delta{\nu}$ vs. $|\Delta{\dot\nu}|$ relation of these five glitches is similar to those with no detected delayed spin-up process, indicating that they are similar to the others in nature except that they have larger amplitudes. For these five glitches, the amplitudes of the delayed spin-up process ($|\Delta{\nu}_{\rm d1}|$) and recovery process ($\Delta{\nu}_{\rm d2}$), their time scales ($\tau_{1}$, $\tau_{2}$), and permanent changes in spin frequency ($\Delta{\nu}_{\rm p}$) and total frequency step ($\Delta{\nu}_{\rm g}$) have positive correlations. From these correlations, we suggest that the delayed spin-up processes are common for all glitches, but are too short and thus difficult to be detected for most glitches.

Gravitational wave emissions (GWEs) of pulsars could not only make them promising targets for continuous gravitational wave searches but also leave imprints in their timing data. We interpret the measured braking index of PSR J1640-4631 with a model involving both the GWE and dipole magnetic field decay. Combining the timing data of PSR J1640-4631 and the theory of magnetic field decay, we propose a new approach of constraining the number of precession cycles, $\xi$, which is highly uncertain currently but can be tightly related to the interior physics of a neutron star and its GWE. We suggest that future observation of the tilt angle $\chi$ of PSR J1640-4631 would not merely help to constrain $\xi$ but also possibly provide information about the internal magnetic field configuration of this pulsar. We find that $\xi$ would be larger than previous estimates unless a tiny angle $\chi\lesssim5^\circ$ is observed. Furthermore, a measured angle $\chi\gtrsim12^\circ$ would indicate $\xi\gtrsim 10^6$, which is at least ten times larger than that suggested previously.

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.

We construct and train an artificial neural network called the back-propagation neural network to describe the evolution of the type Ia supernova spectrum by using the data from the CfA Supernova Program. This network method has many attractive features, and one of them is that the constructed model is differentiable. Benefitting from this, we calculate the absorption velocity and its variation. The model we constructed can well describe not only the spectrum of SNe Ia with wavelength range from $3500\AA$ to $8000\AA$, but also the light-curve evolution with phase time from $-15$ to $50$ with different colors. Moreover, the number of parameters needed during the training process is much less than the usual methods.

Since a large population of massive O/B stars and putative neutron stars (NSs) located in the vicinity of the Galactic center (GC), intermediate-mass X-ray binaries (IMXBs) constituted by a NS and a B-type star probably exist here. We investigate the evolutions of accreting NSs in IMXBs (similar to M82 X-2) with a $\sim5.2M_\odot$ companion, and orbit period $\simeq2.53$ day. By adopting an mildly super-Eddington rate $\dot{M}=6\times10^{-8}M_\odot~{\rm yr}^{-1}$ for the early Case B Roche-lobe overflow (RLOF) accretion, we find only in accreting NSs with quite elastic crusts (slippage factor $s=0.05$), the toroidal magnetic fields can be amplified within 1 Myrs, which is assumed to be the longest duration of the RLOF. These IMXBs will evolve into NS+white dwarf (WD) binaries if they are dynamical stable. However, before the formation of NS+WD binaries, the high stellar density in the GC will probably lead to frequent encounters between the NS+evolved star binaries (in post-early Case B mass transfer phase) and NSs or exchange encounters with other stars, which may produce single NSs. These NSs will evolve into magnetars when the amplified poloidal magnetic fields diffuse out to the NS surfaces. Consequently, our results provide a possible explanation for the origin of the GC magnetar SGR 1745-2900. Moreover, the accreting NSs with $s>0.05$ will evolve into millisecond pulsars (MSPs). Therefore, our model reveals the GC magnetars and MSPs could both originate from the special kind of IMXBs.

Newly born massive magnetars are generally considered to be produced by binary neutron star (NS) mergers, which could give rise to short gamma-ray bursts (SGRBs). The strong magnetic fields and fast rotation of these magnetars make them promising sources for gravitational wave (GW) detection using ground based GW interferometers. Based on the observed masses of Galactic NS-NS binaries, by assuming different equations of state (EOSs) of dense matter, we investigate the stochastic gravitational wave background (SGWB) produced by an ensemble of newly born massive magnetars. The massive magnetar formation rate is estimated through: (i) the SGRB formation rate (hereafter entitled as MFR1); (ii) the NS-NS merger rate (hereafter entitled as MFR2). We find that for massive magnetars with masses $M_{\rm mg}=2.4743 M_\odot$, if EOS CDDM2 is assumed, the resultant SGWBs may be detected by the future Einstein Telescope (ET) even for MFR1 with minimal local formation rate, and for MFR2 with a local merger rate $\dot{\rho}_{\rm c}^{\rm o}(0)\lesssim 10$ ${\rm Mpc}^{-3}{\rm Myr}^{-1}$. However, if EOS BSk21 is assumed, the SGWB may be detectable by the ET for MFR1 with the maximal local formation rate. Moreover, the background spectra show cutoffs at about 350 Hz in the case of EOS BSk21, and at 124 Hz for CDDM2, respectively. We suggest that if the cutoff at $\sim100$ Hz in the background spectrum from massive magnetars could be detected, then the quark star EOS CDDM2 seems to be favorable. Moreover, the EOSs, which present relatively small TOV maximum masses, would be excluded.

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.

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.

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.

We have performed a linear pulsational stability survey of 6 series of long period variable models with M=1.0 Msun, L=3000 - 8000Lsun, and (X,Z)= (0.700,0.020),(0.735,0.005). The dynamic and thermodynamic couplings between convection and oscillations are treated by using a statistical theory of nonlocal and time-dependent convection. The results show that the fundamental and all the low overtones are always pulsationally unstable for the low temperature models when the coupling between convection and oscillations is ignored. When the coupling is considered, there is indeed a "Mira" pulsationally instability region outside of the Cepheid instability strip on the H-R diagram. The coolest models near the Hayashi track are pulsationally stable. Towards high temperature the fundamental mode first becomes unstable, and then the first overtone. Some one of the 2nd -4th overtone may become unstable for the hotter models. All the modes higher than 4th (n > 4) are pulsationally stable. The position and the width of such an instability region on the H-R diagram critically depends on the mass, luminosity and metal abundance of the star.

Using a nonlocal time-dependent theory of convection, we have calculated the linear non-adiabatic oscillations of the Horizontal Branch (HB) stars, with both the dynamic and thermodynamic coupling between convection and oscillations carefully treated. Turbulent pressure and turbulent viscosity have been included consistently in our equations of non-adiabatic pulsation. When the coupling between convection and oscillations is ignored, for all models with Teff <= 7350, the fundamental through the second overtone are pulsationally unstable; while for Teff <= 6200 all the models are unstable up to (at least) the 9th overtone. When the coupling between convection and oscillations is included, the RR Lyrae instability strip is very well predicted. Within the strip the most models are pulsationally unstable only for the fundamental and the first few overtones. The turbulent viscosity is an important damping mechanism. Being exclusively distinct from the luminous red variables (long period variables), the HB stars to the right of the RR strip are pulsationally stable for the fundamental and low-order overtones, but become unstable for some oen of the high-order overtones. This may provide a valuable clue for the short period, low amplitude red variables found outside the red edge of the RR strip on the H-R diagram of globular clusters. Moreover, we present a new radiation modulated excitation mechanism functioning in a zone of radiation flux gradient. The effects of nonlocal convection and the dynamic coupling between convection and oscillations are discussed. The spatial oscillations of the thermal variables in the pulsational calculations have been effectively suppressed.

We relate the MSW effect to the efective absorption of the electronic collective motion energy by retaining the imaginary part of the index of refraction associated with the charged-current scattering and show that the small angle MSW solution to the solar neutrino anomaly can be used as a probe of the physical conditions of the solar interior if it is correct. We find that the constraint on the absorption imposed by the small angle MSW solution and the theoretical estimate of the absorption by the Boltzmann kinetic theory are consistent, which shows that a consistent theoretical picture can be developed when plasma absorption processes are taken into account.