Periodic Variation Studies of the Two Short Period W UMa-type Eclipsing Binaries: LX Lyn and V0853 Aur

In this paper,new light curves(LCs) of contact eclipsing binary(CEB) systems LX Lyn and V0853 Aur are presented and analyzed by using the 2015 version of the Wilson-Devinney(W-D) code.In order to explain their asymmetric LCs,cool starspots on the components were employed.It is suggested that their fill-out degrees are f=12.0%(LX Lyn) and f=26.3%(V0853 Aur).At the same time,we found that LX Lyn is a W-type eclipsing binary(EB) with an orbital inclination of i=84°.88 and a mass ratio of q=2.31.V0853 Aur is also a W-type CEB with a mass ratio of q=2.77 and an orbital inclination of i= 79°.26.Based on all available times of light minimum,their orbital period changes are studied by using the O-C method.The O-C diagram of LX Lyn reveals a cyclic oscillation with a period of about 14.84 yr and an amplitude of 0.0019 days,which can be explained by the light-travel time effect(LTTE) due to the presence of a third body with a minimum mass of0.06M_⊙.For V0853 Aur,it is discovered that the O-C diagram of the system also shows a cyclic oscillation with a period of 9.64 yr and an amplitude of 0.03365 days.The cyclic oscillation of V0853 Aur can be attributed to the LTTE by means of a third body with a mass no less than 3.77M_⊙.The third body may play an important role in the formation and evolution of these systems.

Time-resolved Spectral Properties of Fermi-GBM Bright Long Gamma-Ray Bursts

The prompt emission mechanism of gamma-ray bursts(GRBs)is still unclear,and the time-resolved spectral analysis of GRBs is a powerful tool for studying their underlying physical processes.We performed a detailed time-resolved spectral analysis of 78 bright long GRB samples detected by Fermi/Gamma-ray Burst Monitor.A total of 1490 spectra were obtained and their properties were studied using a typical Band-shape model.First,the parameter distributions of the time-resolved spectrum are given as follows:the low-energy spectral indexα~-0.72,high-energy spectral indexβ~2.42,the peak energy E_(p)~221.69 keV,and the energy flux F~7.49×10^(-6)erg cm^(-2)s^(-1).More than 80%of the bursts exhibit the hardest low-energy spectral indexα_(max),exceeding the synchrotron limit(-2/3).Second,the evolution patterns of a and E_(p)were statistically analyzed.The results show that for multi-pulse GRBs the intensity-tracking pattern is more common than the hard-to-soft pattern in the evolution of both E_(p)andα.The hard-to-soft pattern is generally shown in single-pulse GRBs or in the initial pulse of multi-pulse GRBs.Finally,we found a significant positive correlation between F and E_(p),with half of the samples exhibiting a positive correlation between F andα.We discussed the spectral evolution of different radiation models.The diversity of spectral evolution patterns indicates that there may be more than one radiation mechanism occurring in the GRB radiation process,including photo spheric radiation and synchrotron radiation.However,it may also involve only one radiation mechanism,but more complicated physical details need to be considered.

Luminosity of a radio pulsar and its new emission death line

We investigated the pulsar radio luminosity(L),emission efficiency(ratio of radio luminosity to its spin-down powerE),and death line in the magnetic field(B)versus spin period(P)diagram.We found that the dependence of pulsar radio luminosity on its spin-down power(L-E)is very weak,shown as L^E0.06,which deduces an equivalent inverse correlation between emission efficiency and spin-down power asξ~E-0.94.Furthermore,we examined the distributions of radio luminosity of millisecond and normal pulsars and found that for the similar spin-down powers,the radio luminosity of millisecond pulsars is about one order of magnitude lower than that of the normal pulsars.The analysis of pulsar radio flux suggests that these correlations are not due to a selective effect but are intrinsic to the pulsar radio emission physics.Their radio radiations may be dominated by the different radiation mechanisms.The cutoff phenomenon of currently observed radio pulsars in B-P diagram is usually referred as the"pulsar death line",which corresponds toE≈1030 erg s-1 and is obtained by the cut-off voltage of electron acceleration gap in the polar cap model of pulsar proposed by Ruderman and Sutherland.Observationally,this death line can be inferred by the actual observed pulsar flux S≥1 mJy and 1 kpc distance,together with the maximum radio emission efficiency of 1%.However,the observation data show that the 37 pulsars pass over the death line,including the recently observed two pulsars with long periods of 23.5 s and 12.1 s,which violate the prediction of the polar cap model.At present,the actual observed pulsar flux can reach0.01 mJy by FAST telescope.This will arise the observational limit of spin-down power of pulsars as low as E˙≈1028 erg s-1.This means that the new death line is downward shifted two orders of magnitude,which might be favorably referred as the"observational limit–line".Accordingly,the pulsar theoretical model for the cut-off voltage of gap should be heavily modified.

The “Bi-drifting” Subpulses of PSR J0815+0939 Observed with the Fivehundred-meter Aperture Spherical Radio Telescope

We report the "Bi-drifting" subpulses observed in PSR J0815+0939 using the Five-hundred-meter Aperture Spherical radio Telescope(FAST). The observation at band from 1050 to 1450 MHz is evenly divided into two bands, i.e., the bands at center frequencies of 1150 and 1350 MHz. The mean pulse profiles and the "Bi-drifting"subpulses at these two bands are investigated. It is found that the pulse profiles at these two frequencies show four emission components, and the peak separations between four emission components decrease with the increase of frequency. In addition, the ratio of peak intensity of each component to the intensity of component Ⅳ at 1150 MHz is larger than that at 1350 MHz. We carry out an analysis of the longitude-resolved fluctuation spectrum and twodimensional fluctuation spectrum for each emission component, and find that the P3 of components Ⅰ,Ⅱ and Ⅲ is about 10.56, 10.57 and 10.59 s at 1150 and 1350 MHz. However, the reliable measurements of P3 of component IV and P2 for these four components were not obtained due to the low signal-to-noise ratio of observation data.The pulse energy distributions at frequencies 1150 and 1350 MHz are presented, and it is found that no nulling phenomenon has been found in this pulsar. With our observation from the FAST, the "Bi-drifting" subpulse phenomenon of PSR J0815+0939 is expanded from 400 to 1350 MHz, which is helpful for the relevant researchers to test and constrain the pulsar emission model, especially the model of "Bi-drifting" subpulse.

Investigation of Mode Changing,Pulse Nulling and Subpulse Drifting Properties in the Asymmetric Conal Triple Radio Pulsar B2319+60

We report on a detailed analysis of the mode changing,nulling and subpulse drifting behavior of the conal triple pulsar B2319+60 at 1.5 GHz observed with the Nanshan 25 m radio telescope.The pulsar’s profile can be interpreted as resulting from a sightline traverse which cuts across an outer cone and tangentially grazes an inner cone.About 30 per cent of nulls are found to create alternating bunches of nulls and emission in a quasi-periodic manner with an averaged fluctuation rate of about four rotation periods(P_(1)).The presence of two distinct drift modes(A and ABN)plus a phase-stationary non-drift emission mode(B)is confirmed,and each corresponds to a different pulse profile.The mode A is dominated by a phase modulation in the trailing component with a circulation time of around 7P_(1).In addition to a 3P_(1)phase modulation in the leading component,mode ABN presents an amplitude modulation in the leading and trailing components with a period of around 40P_(1).The emission region and geometry remain constant during mode changing.The diversity seen in the individual pulse behavior of the pulsar B2319+60 provides a unique window into the emission physics.

Results of 23 yr of Pulsar Timing of PSR J1453-6413

In this paper,we presented the 23.3 yr of pulsar timing results of PSR J1456-6413 based on the observations of Parkes 64 m radio telescope.We detected two new glitches at MJD 57093(3)and 59060(12)and confirmed its first glitch at MJD 54554(10).The relative sizes(Δν/ν)of these two new glitches are 0.9×10^(-9)and 1.16×10^(-9),respectively.Using the“Cholesky”timing analysis method,we have determined its position,proper motion,and two-dimensional transverse velocities from the data segments before and after the second glitch,respectively.Furthermore,we detected exponential recovery behavior after the first glitch,with a recovery timescale of approximately 200 days and a corresponding exponential recovery factor Q of approximately 0.15(2),while no exponential recovery was detected for the other two glitches.More interestingly,we found that the leading component of the integral pulse profile after the second glitch became stronger,while the main component became weaker.Our results will expand the sample of pulsars with magnetosphere fluctuation triggered by the glitch event.

Applying Hybrid Clustering in Pulsar Candidate Sifting with Multi-modality for FAST Survey

Pulsar search is always the basis of pulsar navigation,gravitational wave detection and other research topics.Currently,the volume of pulsar candidates collected by the Five-hundred-meter Aperture Spherical radio Telescope(FAST)shows an explosive growth rate that has brought challenges for its pulsar candidate filtering system.Particularly,the multi-view heterogeneous data and class imbalance between true pulsars and non-pulsar candidates have negative effects on traditional single-modal supervised classification methods.In this study,a multi-modal and semi-supervised learning based on a pulsar candidate sifting algorithm is presented,which adopts a hybrid ensemble clustering scheme of density-based and partition-based methods combined with a feature-level fusion strategy for input data and a data partition strategy for parallelization.Experiments on both High Time Resolution Universe SurveyⅡ(HTRU2)and actual FAST observation data demonstrate that the proposed algorithm could excellently identify pulsars:On HTRU2,the precision and recall rates of its parallel mode reach0.981 and 0.988 respectively.On FAST data,those of its parallel mode reach 0.891 and 0.961,meanwhile,the running time also significantly decreases with the increment of parallel nodes within limits.Thus,we can conclude that our algorithm could be a feasible idea for large scale pulsar candidate sifting for FAST drift scan observation.

Ancient subsurface structure beneath crater Clavius:constraint by recent high-precision gravity and topography data

With the increasing precision of the GRAIL gravity field models and topography from LOLA, it is possible to investigate the substructure beneath crater Clavius. An admittance between gravity and topography data is commonly used to estimate selenophysical parameters, including load ratio, crustal thickness and density, and elastic thickness. Not only a surface load, but also a subsurface load is considered in estimation. The algorithm of particle swarm optimization(PSO) with a swarm size of 400 is employed as well.Results indicate that the observed admittance is best-fitted by the modeled admittance based on a spherical shell model, which was proved to be unsatisfactory in the previous study. The best-fitted load ratio f is around-0.194. Such a small load ratio conforms to the direct proportion between the nearly uncompensated topography and its corresponding negative gravity anomaly. It also indicates that a surface load dominates all the loads. Constrained within 2σSTD, a small crustal thickness(～30 km) and a crustal density of ～2587 kg m-3are found, quite close to the results from previous GRAIL research. Considering the well constrained crustal thickness and density, the best-fitted elastic thickness(～7 km) is rational. This result is slightly smaller than the previous study(～12 km). Such difference can be attributed to the difference in crustal density used and the precision of gravity and topography data. Considering that the small difference between the modeled gravity anomaly and observations is quite small, a parameter inversed here could be an indicator of the subsurface structure beneath Clavius.

Effects of the momentum dependence of nuclear symmetry potential on pion observables in Sn+Sn collisions at 270 MeV/nucleon

Within a transport model,we investigated the effects of the momentum dependence of the nuclear symmetry potential on the pion observables in central Sn+Sn collisions at 270 MeV/nucleon.To this end,the quantity U^(∞)_(sym)(ρ_(0))(i.e.,the value of the nuclear symmetry potential at the saturation densityρ_(0) and infinitely large nucleon momentum)was used to characterize the momentum dependence of the nuclear symmetry potential.With a certain L(i.e.,the slope of the nuclear symmetry energy at(ρ_(0))),the characteristic parameter U^(∞)_(sym)(ρ_(0))of the symmetry potential significantly affects the production of−and+and their pion ratios.Moreover,by comparing the charged pion yields,pion ratios,and spectral pion ratios of the theoretical simulations for the reactions ^(108) Sn+^(112) Sn and ^(132)Sn+^(124)Sn with the corresponding data in the SRIT experiments,we found that our results favor a constraint on U^(∞)_(sym)(ρ_(0))(i.e.,−160_(−9)^(+18) MeV),and L is also suggested within a range of 62.7 MeV<L<93.1 MeV.In addition,the pion observable for^(197)Au+^(197)Au collisions at 400 MeV/nucleon also supports the extracted value for U^(∞)_(sym)(ρ_(0)).

Examination of an isospin-dependent single-nucleon momentum distribution for isospin-asymmetric nuclear matter in heavy-ion collisions

Within a transport model using nucleon momentum profiles as the input from a parameterized isospin-dependent single-nucleon momentum distribution,with a high momentum tail induced by short-range correlations,we employ197 Au+197 Au collisions at 400 MeV/nucleon to examine the effects of the short-range correlations on the pion and flow observables in probing the nuclear symmetry energy.We investigate how reliable this isospin-dependent single-nucleon momentum distribution is and determine the corresponding parameter settings.Apart from the significant effects of the short-range correlations on the pion and flow observables that are observed,we also find that the theoretical simulations of the197 Au+197 Au collisions with this momentum distribution using two sets of parameters,extracted from the experimental analysis and the self-consistent Green’s function prediction,can reproduce the neutron ellipticflows of the FOPI-LAND experiment and the pà=ptratios of the FOPI experiment,respectively,under the symmetry energy setting in a particular range.Therefore,we conclude that this parameterized isospin-dependent single-nucleon momentum distribution is reliable for isospin-asymmetric nuclear matter.Correspondingly,two sets of parameters extracted from both the experimental analysis and the selfconsistent Green’s function prediction cannot be excluded according to the available experimental information at present.

On the intrinsic shape of the gamma-ray spectrum for Fermi blazars

The curvature of the γ-ray spectrum in blazars may reflect the intrinsic distribution of emitting electrons, which will further give some information on the possible acceleration and cooling processes in the emitting region. The γ-ray spectra of Fermi blazars are normally fitted either by a single power-law（PL） or a log-normal（call Logarithmic Parabola, LP） form. The possible reason for this difference is not clear. We statistically explore this issue based on the different observational properties of 1419 Fermi blazars in the 3 LAC Clean Sample. We find that the γ-ray flux（100 Me V–100 Ge V） and variability index follow bimodal distributions for PL and LP blazars, where the γ-ray flux and variability index show a positive correlation. However, the distributions of γ-ray luminosity and redshift follow a unimodal distribution. Our results suggest that the bimodal distribution of γ-ray fluxes for LP and PL blazars may not be intrinsic and all blazars may have an intrinsically curved γ-ray spectrum, and the PL spectrum is just caused by the fitting effect due to less photons.

Nucleonic ~1S_0 Superfluidity Induced by a Soft Pion in Neutron Star Matter with Antikaon Condensations

The nucleonic ~1S_0 superfluidity is investigated by solving the gap equation for the Reid soft-core potential as the nucleon-nucleon interaction in neutron star(NS) matter which is considered to be made up of p, e, μ and condensed antikaon matter. We mainly study the influence of the soft pion-induced potential on the nucleonic^1S_0 pairing gaps in the above NS matter. It is found that the intensities of the nucleonic ~1S_0 pairing gaps including the soft pion-induced potential are smaller than those calculated in the case of not including the soft pion-induced potential. Furthermore, the nucleonic ~1S_0 pairing gaps with the soft pion-induced potential fall into decline with the deepening of the optical potential of antikaons in the above NS matter, whereas they increase with the parameter η for the fixed optical potential of antikaons. Due to the appearance of the soft pion-induced potential, the maximum values of nucleonic ~1S_0 pairing gaps at parameter η = 0.20,0.55 are suppressed by1.7%-6.8% with respect to the case without soft pion-induced potential in the above NS matter.

Subpulse Drifting of PSR J1110–5637

We report a detailed study of polarization characteristics and subpulse drifting in PSR J1110-5637 with the observations of the Parkes 64 m radio telescope at 1369 MHz.The observations revealed that the trailing component of the pulse profile has obvious subpulse drifting,while the leading component has no subpulse drifting.Using the two-dimensional fluctuation spectrum(2DFS),we detected three distinct emission modes in the trailing component(modes A,B and C).The emission in mode A is chaotic and indistinguishable,while modes B and C have obvious subpulse drifting.The vertical modulation periods P3of modes B and C are around the mean values of 12 P and 8 P,respectively.The subpulse drifting of PSR J1110-5637 will expand the pulsar sample with multiple subpulse drifting rates,and this will help future systematic studies on the physical origin of the subpulse drifting phenomenon.

The Distribution of UV Radiation Field in the Molecular Clouds of Gould Belt

The distribution of ultraviolet(UV)radiation field provides critical constraints on the physical environments of molecular clouds.Within 1 kpc of our solar system and fostering protostars of different masses,the giant molecular clouds in the Gould Belt present an excellent opportunity to resolve the UV field structure in star-forming regions.We performed spectral energy distribution(SED)fitting of the archival data from the Herschel Gould Belt Survey(HGBS).Dust radiative transfer analysis with the DUSTY code was applied to 23 regions in 14 molecular complexes of the Gould Belt,resulting in the spatial distribution of the radiation field in these regions.For 10 of 15 regions with independent measurements of star formation rate,their star formation rate and UV radiation intensity largely conform to a linear correlation found in previous studies.

A discovery of two slow pulsars with FAST:“Ronin”from the globular cluster M15

Globular clusters harbor numerous millisecond pulsars,but long-period pulsars(P 100 ms)are rarely found.In this study,we employed a fast folding algorithm to analyze observational data from multiple globular clusters obtained by the Five-hundredmeter Aperture Spherical radio Telescope(FAST),aiming to detect the existence of long-period pulsars.We estimated the impact of the median filtering algorithm in eliminating red noise on the minimum detectable flux density(S_(min))of pulsars.Subsequently,we successfully discovered two isolated long-period pulsars in M15 with periods approximately equal to 1.928451 and3.960716 s,respectively.On the P-˙P diagram,both pulsars are positioned below the spin-up line,suggesting a possible history of partial recycling in X-ray binary systems disrupted by dynamical encounters later on.According to timing results,these two pulsars exhibit remarkably strong magnetic fields.If the magnetic fields were weakened during the accretion process,then a short duration of accretion might explain the strong magnetic fields of these pulsars.

Long-term variations of X-ray pulse profiles for the Crab pulsar:data analysis and modeling

The pulse profiles of the Crab pulsar(as well as some other pulsars)vary with time.They can lead to a major source of intrinsic timing noise,which lacks a detailed physical model.The phase separation?between the first left peak(P1)and the second right peak(P2)is a key parameter that shows the variations of pulse profiles for the Crab pulsar.It was found that the evolution of?has a tendency with increasing rates of 0.82?±0.25?,0.80?±0.54?,and 0.77?±0.28?per century for the 2-6,6-15,and15-60 ke V bands,respectively.Furthermore,the flux ratios(P2/P1)of X-ray pulse profiles in the three bands were calculated,and the derived flux ratios were consistent with the radio and X-ray measurements of the Insight-HXMT.In addition to discovering the physical origin of the pulse changes,the high-SNR X-ray pulse profiles were simulated in the annular gap model,and two model parameters(e.g.,the maximum emission heights of the two peaks)were observed to slightly affect the variations of peak separation.We fitted the long-term variations of emission heights of the two peaks and discovered that the emission heights showed increasing tendencies with time.Variations of these emission heights induced a characteristic period derivative,and a complete formula for both the magnetic dipole radiation and wind-particle-induced variations of the moment of inertia was used for the pulsar’s spin-down to obtain the variation rate˙αof the magnetic inclination angle,which was-1.60?per century.Intrinsic timing noise is observed to be mainly induced by the variations of pulse profiles,which might correlate with a characteristic spin period derivative arising from the fluctuations of the emission regions.This work will lay a foundation for understanding the origin of intrinsic timing noise and making high-precision timing models in the future.

CSS J075415.6+191052 and NW Leo: Two Contact Binaries at Different Evolutionary Stages

Multi-color light curves of CSS J075415.6+191052 and NW Leo are presented and the photometric solutions suggest that CSS J075415.6+191052 is an A-subtype contact binary with low mass ratio(q = 0.178) while NW Leo has a high mass ratio(q = 0.707). For CSS J075415.6+191052, the RI light curves show weakening around the left shoulder of the secondary minimum, which indicates that there may be a dark spot on the secondary component. However, the light curves of BV bands are totally symmetric. It is unreasonable if the dark spot is caused by magnetic activity or mass transfer between the two components. Therefore, weakening of the light curves in this contact binary is caused by something else. A possible explanation is mass transferring from the primary component to the common convective envelope through the inner Lagrangian point, and this part of the mass, for some reason, weakens the RI bands of light from the secondary component. O-C analysis of NW Leo reveals a cyclic period change with a modulation period of 4.7 yr, which may be caused by the light travel time effect of a third body. The positions of CSS J075415.6+191052 and NW Leo in the P–J′_(orb)diagram indicate that CSS J075415.6+191052 mainly abides by the angular momentum loss theory while NW Leo is dominated by the thermal relaxation oscillation theory.

1SWASP J034439.97+030425.5:a short-period eclipsing binary system with a close-in stellar companion

First multi-wavelength photometric light curves(LCs)of the short-period eclipsing binary(EB)1 SWASP J034439.97+030425.5(hereafter J0344)are presented and analyzed by using the 2013 version of the Wilson-Devinney(W-D)code.To explain the asymmetric LCs of J0344,a cool star-spot on the less massive component was employed.The photometric solutions suggest that J0344 is a W-subtype shallow contact EB with a contact degree of f=4.9%±3.0%and a mass ratio of q=2.456±0.013.Moreover,an obvious third light was detected in our analysis.We calculated the average luminosity contribution of the third light to the total light,and that value reaches up to 49.78%.Based on the O-C method,the variations of the orbital period were studied for the first time.Our O-C diagram reveals a secular decrease superimposed on a cyclic oscillation.The orbital period decreases at a rate of d P/dt=-6.07×10-7 d yr-1,which can be explained by the mass transfer from the more massive component to the less massive one.Besides,its O-C diagram also shows a cyclic oscillation with an amplitude of 0.0030 d and a period about 7.08 yr,which can be explained by the presence of a third body with a minimum mass of M3 min=0.15±0.02 M⊙.The third component may play an important role in the formation and evolution of J0344 by drawing angular momentum from the central system.

Spin-down and emission variations for PSR J0742-2822

PSR J0742-2822 is known for its quasi-periodic changes in the observed pulse profile and spindown rate.In this paper,we analyzed 13 years of timing data obtained with the Nanshan 25-m radio telescope and the Parkes 64-m radio telescope.We found that the average values of the spin-down rate((v))of this pulsar changed in four different states.We investigated the correlation between v and W50,and ascertained that the correlation changed in different states.Moreover,not all the changes in states and correlation can be associated with glitch activities.We examined the long term evolution ofγ-ray flux(0.1-300 GeV)and the pulse profiles corresponding to the four different states using Fermi-LAT Pass8(P8 R3)data from 2008 August 5 to 2019 October 1.We did not detect a significant change inγ-ray flux or the pulse profile.Our results suggest that the connection between pulsar rotation and emission is more complex than previously reported for this pulsar.

The contribution of the first forbidden transitions to the nuclear β^--decay half-life

β-decay half-life is a key quantity for nuclear structure and nucleosynthesis studies.There exist large uncertainties in the contributions of allowed and forbidden transitions to the totalβ-decay life,which limits the resolution of the predictedβ-decay half-life.We systematically study the contribution of the first forbidden(FF)transitions to theβ--decay half-life,and quantify it with a formula based on simple physics considerations.We also propose a new formula for calculation of theβ--decay half-life that includes the FF contribution.It is shown that the inclusion of the contribution of FF transitions significantly improves the precision of calculations of theβ--decay half-life.By fitting of the RQRPA results for neutron-rich Z=47,57 isotopes and N=80,94 isotones,the formula for the contribution of the FF transitions gives similar results as the RQRPA calculations.However,because of limited experimental data for the branching ratios of unstable nuclei,the fit parameters are not fully constrained.Therefore,the proposed formula for theβ--decay half-life is more suitable for calculations of half-lives than of the FF contributions.The formula could be used to predict theβ?-decay half-life in nuclear structure studies as well as nucleosynthesis calculations in stars.