Showing 1–50 of 50 results for author: Lü, G

Novae: An Important Source of Lithium in the Galaxy

Authors: Jun Gao, Chunhua Zhu, Guoliang Lü, Jinlong Yu, Lin Li, Helei Liu, Sufen Guo

Abstract: The source of the Galactic Lithium (Li) has long been a puzzle. With the discovery of Li in novae, extensive research has been conducted. However, there still exists a significant disparity between the observed abundance of lithium in novae and the existing theoretical predictions. Using the Modules for Experiments in Stellar Astrophysics (MESA), we simulate the evolution of nova with element diff… ▽ More The source of the Galactic Lithium (Li) has long been a puzzle. With the discovery of Li in novae, extensive research has been conducted. However, there still exists a significant disparity between the observed abundance of lithium in novae and the existing theoretical predictions. Using the Modules for Experiments in Stellar Astrophysics (MESA), we simulate the evolution of nova with element diffusion and appropriately increased the amount of 3^He in the mixtures. Element diffusion enhances the transport efficiency between the nuclear reaction zone and the convective region on the surface of the white dwarf during nova eruptions, which results in more 7^Be to be transmitted to the white dwarf surface and ultimately ejected. Compared to the previous predictions, the abundance of 7^Be in novae simulated in our model significantly increases. And the result is able to explain almost all observed novae. Using the method of population synthesis, we calculate Li yield in the Galaxy. We find that the Galactic occurrence rate of nova is about 130 yr^{-1}, and about 110M Li produced by nova eruption is ejected into the interstellar medium (ISM). About 73\% of Li in the Galactic ISM originates from novae, and approximately 15\%-20\% of the entire Galaxy. It means that novae are the important source of Li in the Galactic. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Comments: 12 pages, 4 figures. Accepted for publication in Astrophysical Journal

The population synthesis of Wolf-Rayet stars involving binary merger channels

Authors: Zhuowen Li, Chunhua Zhu, Guoliang Lü, Lin Li, Helei Liu, Sufen Guo, Jinlong Yu, Xizhen Lu

Abstract: Wolf-Rayet stars (WRs) are very important massive stars. However, their origin and the observed binary fraction within the entire WR population are still debated. We investigate some possible merger channels for the formation of WRs, including main sequence (MS)/ Hertzsprung Gap (HG) + MS, He + HG/ Giant Branch (GB). We find that many products produced via binary merger can evolve into WRs, the MS… ▽ More Wolf-Rayet stars (WRs) are very important massive stars. However, their origin and the observed binary fraction within the entire WR population are still debated. We investigate some possible merger channels for the formation of WRs, including main sequence (MS)/ Hertzsprung Gap (HG) + MS, He + HG/ Giant Branch (GB). We find that many products produced via binary merger can evolve into WRs, the MS/ HG + MS merger channel can explain WRs with luminosities higher than $\sim 10^{5.4}$\,L$_{\odot}$, while the He + HG/ GB merger channel can explain low-luminosity WRs in the range of $10^{4.7}$\,L$_{\odot}$\,$\sim$\,$10^{5.5}$\,L$_{\odot}$. In the population synthesis analysis of WRs, we assume an initial binary fraction ($f_{\rm ini,bin}$) of 50\% and 100\% for massive stars. We also assume that MS/ HG + MS merger products are non-rotating or rapidly rotating ($ω/ω_{\rm crit}=0.8$). In different cases, the calculated single fractions of WRs range from $22.2\%$ to $60.6\%$ in the Milky Way (MW) and from $8.3\%$ to $70.9\%$ in the Large Magellanic Cloud (LMC). The current observations fall within the range of our calculations. When the merger product of MS/HG+MS rotates rapidly, we estimate that there are approximately 1015 to 1396 WRs in the MW and 128 to 204 WRs in the LMC. Our model also roughly reproduces the observed single-peak luminosity distribution of WRs in the MW. However, the weak bimodal luminosity distribution observed in the LMC is not reproduced in our model. We assess that this may be due to the model underestimating the mass-loss rate in the LMC. In conclusion, we consider that the binary merger is significant formation channel for WR formation, and can explain the observed high fraction of the single WRs in the total population. △ Less

Submitted 19 May, 2024; originally announced May 2024.

Comments: 14 pages, 8 figures , Accepted to APJ

Reconstructing the Baryonic Acoustic Oscillations in the presence of photo-$z$ uncertainties

Authors: Kwan Chuen Chan, Guoyuan Lu, Xin Wang

Abstract: The reconstruction method has been widely employed to improve the Baryon Acoustic Oscillations (BAO) measurement in spectroscopic survey data analysis. In this study, we explore the reconstruction of the BAO signals in the realm of photometric data. By adapting the Zel'dovich reconstruction technique, we develop a formalism to reconstruct the transverse BAO in the presence of photo-$z$ uncertainti… ▽ More The reconstruction method has been widely employed to improve the Baryon Acoustic Oscillations (BAO) measurement in spectroscopic survey data analysis. In this study, we explore the reconstruction of the BAO signals in the realm of photometric data. By adapting the Zel'dovich reconstruction technique, we develop a formalism to reconstruct the transverse BAO in the presence of photo-$z$ uncertainties \change{under the plane-parallel approximation}. We access the performance of the BAO reconstruction through comoving $N$-body simulations. The transverse reconstruction potential can be derived by solving a 2D potential equation, with the surface density and the radial potential contribution acting as the source terms. The solution is predominantly determined by the surface density. As is evident in dense samples, such as the matter field, the transverse BAO reconstruction can enhance both the strength of the BAO signals and their cross correlation with the initial conditions. At $z=0$, the cross correlation is increased by a factor of 1.2 at $ k_\perp = 0.2 \, \mathrm{Mpc}^{-1}h $ and 1.4 at $ k_\perp = 0.3 \, \mathrm{Mpc}^{-1}h $, respectively. We contrast the 2D potential results with the 3D Poisson equation solution, wherein we directly solve the potential equation using the position in photo-$z$ space, and find good agreement. Additionally, we examine the impact of various conditions, such as the smoothing scales and the level of photo-$z$ uncertainties, on the reconstruction results. We envision the straightforward application of this method to survey data. △ Less

Submitted 13 March, 2024; v1 submitted 21 November, 2023; originally announced November 2023.

Comments: 15 pages, 13 figures, matched to the published version

Evolutionary tracks of massive stars with different rotation and metallicity in neutrino H-R diagram

Authors: Hao Wang, Chunhua Zhu, Helei Liu, Sufen Guo, Guoliang Lü

Abstract: Neutrino losses play a crucial role in the evolution of massive stars. We study the neutrino luminosity of stars ranging from 20 to 90 M_{\odot} from Zero Age Main Sequence (ZAMS) to Fe Core Collapse (FeCC) with different rotation and metallicity in a neutrino Hertzsprung-Russell diagram. In our simulations, we consider ω/ωcrit = 0 and 0.7 to represent non-rotation and high rotation, respectively,… ▽ More Neutrino losses play a crucial role in the evolution of massive stars. We study the neutrino luminosity of stars ranging from 20 to 90 M_{\odot} from Zero Age Main Sequence (ZAMS) to Fe Core Collapse (FeCC) with different rotation and metallicity in a neutrino Hertzsprung-Russell diagram. In our simulations, we consider ω/ωcrit = 0 and 0.7 to represent non-rotation and high rotation, respectively, and set the metallicities to 0.014, 0.001, and 0.0001. During hydrogen burning stages, neutrino luminosity primarily originates from CNO cycle, and increases with higher stellar mass while decreasing with increasing metallicity. For the high metallicity models (Z = 0.014) during the helium burning stage, the reduction of the hydrogen envelope caused by a larger mass loss rate leads to a gradual decrease in neutrino luminosity. The rapid rotation results in extra mixing inside massive stars, which increases the neutrino luminosity during main sequence (MS), while decreases the neutrino luminosity during helium burning phase. Simultaneously, the rapid rotation also increases CO core mass, which enhances the neutrino luminosity during C and O burning phase. We also investigate the effect of neutrino magnetic moment (NMM) on the massive stars. We find that the energy loss caused by the NMM does not have effects on the evolutionary destiny of massive stars, and it does not significant change the compactness at the time of Fe core collapse. △ Less

Submitted 16 October, 2023; v1 submitted 4 October, 2023; originally announced October 2023.

Comments: This paper has been accepted by Monthly Notices of the Royal Astronomical Society on October 3, 2023

An Alternative Formation Scenario for Uranium-rich Giants: Engulfing a Earth-like Planet

Authors: Dian Xie, Chunhua Zhu, Sufen Guo, Helei Liu, Guoliang Lü

Abstract: The actinides, such as the uranium (U) element, are typically synthesized through the rapid neutron-capture process (r-process), which can occur in core-collapse supernovae or double neutron star mergers. There exist nine r-process giant stars exhibiting conspicuousUabundances, commonly referred to as U-rich giants. However, the origins of these U-rich giants remain ambiguous. We propose an altern… ▽ More The actinides, such as the uranium (U) element, are typically synthesized through the rapid neutron-capture process (r-process), which can occur in core-collapse supernovae or double neutron star mergers. There exist nine r-process giant stars exhibiting conspicuousUabundances, commonly referred to as U-rich giants. However, the origins of these U-rich giants remain ambiguous. We propose an alternative formation scenario for these U-rich giants whereby a red giant (RG) engulfs an Earth-like planet. To approximate the process of a RG engulfing an Earth-like planet, we employ an accretion model wherein the RG assimilates materials from said planet. Our findings demonstrate that this engulfment event can considerably enhance the presence of heavy elements originating from Earth-like planets on the surfaces of very metal-poor stars (Z = 0.00001), while its impact on solar-metallicity stars is comparatively modest. Importantly, the structural and evolutionary properties of both very metalpoor and solar-metallicity stars remain largely unaffected. Notably, our engulfment model effectively accounts for the observed U abundances in known U-rich giants. Furthermore, the evolutionary trajectories of U abundances on the surfaces of RGs subsequent to the engulfment of Earth-like planets encompass all known U-rich giants. Therefore, it is plausible that U-rich giants are formed when a RG engulfs an Earth-like planet. △ Less

Submitted 12 July, 2023; originally announced July 2023.

Comments: 9 pages, 8 figures, accepted 2023 July 10 by MNRAS

The effects of rotation, metallicity and magnetic field on the islands of failed supernovae

Authors: Lei Li, Chunhua Zhu, Sufen Guo, Helei Liu, Guoliang Lu

Abstract: Failed supernovae (FSN) are a possible channel for the formation of heavy stellar-mass black holes ($M_{ BH}>\sim 30$ M$_\odot$). However, the effects of metallicity, rotation and magnetic field on the islands of explodabilty of massive stars are not clear. Here, we simulate the stellar structure and evolution in the mass range between 6 and 55 $M_{\odot}$ with different initial rotational velocit… ▽ More Failed supernovae (FSN) are a possible channel for the formation of heavy stellar-mass black holes ($M_{ BH}>\sim 30$ M$_\odot$). However, the effects of metallicity, rotation and magnetic field on the islands of explodabilty of massive stars are not clear. Here, we simulate the stellar structure and evolution in the mass range between 6 and 55 $M_{\odot}$ with different initial rotational velocities, metallicities, and magnetic fields from zero-age main sequence (ZAMS) to pre-collapse. We find that the rapid rotating stars can remain lower $\rm ^{12}C$ mass fraction at the time of C ignition, which allows the transition, from convective carbon burning to radiative burning, to occur at lower $M_{\rm ZAMS}$ than those from stars without rotation. However, the rapid rotation is unfavorable for FSN occurring but is conducive to long gamma-ray bursts (lGRBs) because it results in the specific angular momentum in the CO core is greater than the last stable orbit at core collapse. The increasing metallicity does not affect FSN islands, but high metallicity inhibits rotational mixing and is unfavorable for producing lGRBs. A magnetic field can constrain the mass-loss rate even for rapid rotating stars, resulting in higher mass at pre-collapse. The magnetic braking triggered by the magnetic field can reduce the rotation velocity for high-metallicity models, which decreases the specific angular momentum in the CO core and is favorable for FSN occurring. We suggest that the heavy-mass black holes detected by LIGO may originate from rapidly rotating massive stars with strong magnetic fields, rather than those with very low metallicity. △ Less

Submitted 27 June, 2023; originally announced June 2023.

Comments: 16 pages, 11 figures

The First GECAM Observation Results on Terrestrial Gamma-ray Flashes and Terrestrial Electron Beams

Authors: Y. Zhao, J. C. Liu, S. L. Xiong, W. C. Xue, Q. B. Yi, G. P. Lu, W. Xu, F. C. Lyu, J. C. Sun, W. X. Peng, C. Zheng, Y. Q. Zhang, C. Cai, S. Xiao, S. L. Xie, C. W. Wang, W. J. Tan, Z. H. An, G. Chen, Y. Q. Du, Y. Huang, M. Gao, K. Gong, D. Y. Guo, J. J. He , et al. (37 additional authors not shown)

Abstract: Gravitational-wave high-energy Electromagnetic Counterpart All-sky Monitor (GECAM) is a space-borne instrument dedicated to monitoring high-energy transients, including Terrestrial Gamma-ray Flashes (TGFs) and Terrestrial Electron Beams (TEBs). We implemented a TGF/TEB search algorithm for GECAM, with which 147 bright TGFs, 2 typical TEBs and 2 special TEB-like events are identified during an effe… ▽ More Gravitational-wave high-energy Electromagnetic Counterpart All-sky Monitor (GECAM) is a space-borne instrument dedicated to monitoring high-energy transients, including Terrestrial Gamma-ray Flashes (TGFs) and Terrestrial Electron Beams (TEBs). We implemented a TGF/TEB search algorithm for GECAM, with which 147 bright TGFs, 2 typical TEBs and 2 special TEB-like events are identified during an effective observation time of $\sim$9 months. We show that, with gamma-ray and charged particle detectors, GECAM can effectively identify and distinguish TGFs and TEBs, and measure their temporal and spectral properties in detail. A very high TGF-lightning association rate of $\sim$80\% is obtained between GECAM and GLD360 in east Asia region. △ Less

Submitted 17 June, 2023; originally announced June 2023.

Comments: The paper was accepted by Geophysical Research Letters on June 16th, 2023

RRAT J1913+1330: an extremely variable and puzzling pulsar

Authors: S. B. Zhang, J. J. Geng, J. S. Wang, X. Yang, J. Kaczmarek, Z. F. Tang, S. Johnston, G. Hobbs, R. Manchester, X. F. Wu, P. Jiang, Y. F. Huang, Y. C. Zou, Z. G. Dai, B. Zhang, D. Li, Y. P. Yang, S. Dai, C. M. Chang, Z. C. Pan, J. G. Lu, J. J. Wei, Y. Li, Q. W. Wu, L. Qian , et al. (4 additional authors not shown)

Abstract: Rotating Radio Transients (RRATs) are neutron stars that emit sporadic radio bursts. We detected 1955 single pulses from RRAT J1913+1330 using the 19-beam receiver of the Five-hundred-meter Aperture Spherical Radio Telescope (FAST). These pulses were detected in 19 distinct clusters, with 49.4% of them occurring sequentially with a waiting time of one rotation period. The energy distribution of th… ▽ More Rotating Radio Transients (RRATs) are neutron stars that emit sporadic radio bursts. We detected 1955 single pulses from RRAT J1913+1330 using the 19-beam receiver of the Five-hundred-meter Aperture Spherical Radio Telescope (FAST). These pulses were detected in 19 distinct clusters, with 49.4% of them occurring sequentially with a waiting time of one rotation period. The energy distribution of these individual pulses exhibited a wide range, spanning three orders of magnitude, reminiscent of repeating fast radio bursts (FRBs). Furthermore, we observed abrupt variations in pulse profile, width, peak flux, and fluence between adjacent sequential pulses. These findings suggest that this RRAT could be interpreted as pulsars with extreme pulse-to-pulse modulation. The presence of sequential pulse trains during active phases, along with significant pulse variations in profile, fluence, flux, and width, should be intrinsic to a subset of RRATs. Our results indicate that J1913+1330 represents a peculiar source that shares certain properties with populations of nulling pulsars, giant pulses, and FRBs from different perspectives. The dramatic pulse-to-pulse variation observed in J1913+1330 could be attributed to unstable pair creation above the polar cap region and the variation of the site where streaming pairs emit coherently. Exploring a larger sample of RRATs exhibiting similar properties to J1913+1330 has the potential to significantly advance our understanding of pulsars, RRATs, and FRBs. △ Less

Submitted 5 June, 2023; originally announced June 2023.

Comments: 12 pages, 10 figures, 1 table

Hydrogen-free Wolf-Rayet stars: Helium stars with envelope-inflation structure and rotation

Authors: Xizhen Lu, Chunhua Zhu, Helei Liu, Sufen Guo, Jinlong Yu, Guoliang Lü

Abstract: Observations have shown that the effective temperature of hydrogen-free Wolf-Rayet (WR) stars is considerably lower than that of the standard model, which means that the radius of the observed H-free WR stars is several times larger than that estimated by the standard model. The envelope inflation structure (EIS) caused by the radiation luminosity being close to the Eddington luminosity in the iro… ▽ More Observations have shown that the effective temperature of hydrogen-free Wolf-Rayet (WR) stars is considerably lower than that of the standard model, which means that the radius of the observed H-free WR stars is several times larger than that estimated by the standard model. The envelope inflation structure (EIS) caused by the radiation luminosity being close to the Eddington luminosity in the iron opacity peak region of H-free WR stars may be the key to resolve the radius problem of H-free WR stars. We try to explain the H-free WR stars observed in the Milk Way (MW) and the Large Magellanic Cloud (LMC) by the He stars. Using the Modules for Experiments in Stellar Astrophysics code, we compute the evolution of He stars with and without MLT++ prescriptions and discuss their effects on the EIS. We have calculated the evolution of He stars using a new mass-loss rate formula and three different relative rotational velocity and compared our results with observations on Hertzsprung-Russell diagrams. The low luminosity (log$(L/L_{\odot})\leq5.2$) H-free WR stars in the MW and the LMC can be explained by the helium giant phase in low-mass He stars, the high $X_{C}$ and $X_{O}$ in WC stars can only evolve through low-mass He stars with a rapid rotation. High-mass He stars with the EIS can explain H-free WR stars with a luminosity exceeding $10^{5.7} L_{\odot}$ and an effective temperature above $10^{4.7}$ K in the MW. They can also explain H-free WR stars on the right-hand side of the He zero-age main sequence in the LMC. High-mass stars with the EIS evolve into WO stars at the final evolution stage, and the shorter lifetime fraction is consistent with the small number of observed WO stars. △ Less

Submitted 12 April, 2023; originally announced April 2023.

Comments: 9 pages, 7 figures 1 tables, Accepted to A&A

Journal ref: A&A 674, A216 (2023)

Formation and Destiny of White Dwarf and Be Star Binaries

Authors: ChunHua Zhu, GuoLiang Lü, Xizhen Lu, Jie He

Abstract: The binary systems consisting of a Be star and a white dwarf (BeWDs) are very interesting.They can originate from the binaries composed of a Be star and a subdwarf O or B star (BesdOBs), and they can merge into red giants via luminous red nova or can evolve into double WD potentially detected by $LISA$ mission. Using the method of population synthesis, we investigate the formation and the destiny… ▽ More The binary systems consisting of a Be star and a white dwarf (BeWDs) are very interesting.They can originate from the binaries composed of a Be star and a subdwarf O or B star (BesdOBs), and they can merge into red giants via luminous red nova or can evolve into double WD potentially detected by $LISA$ mission. Using the method of population synthesis, we investigate the formation and the destiny of BeWDs,and discuss the effects of the metallicity ($Z$) and the common envelope evolution parameters. We find that BesdOBs are significant progenitors of BeWDs. About 30\% ($Z=0.0001$)-50\% ($Z=0.02$) of BeWDs come from BesdOBs. About 60\% ($Z=0.0001$) -70\% ($Z=0.02$) of BeWDs turn into red giants via a merger between a WD and a non-degenerated star. About 30\% ($Z=0.0001$) -40\% ($Z=0.02$) of BeWDs evolve into double WDs which are potential gravitational waves of $LISA$ mission at a frequency band between about $3\times10^{-3}$ and $3\times10^{-2}$ Hz. The common envelope evolution parameter introduces an uncertainty with a factor of about 1.3 on BeWD populations in our simulations. △ Less

Submitted 5 April, 2023; originally announced April 2023.

Comments: 17 pages, 12 figures, 2 table, accepted for publication in RAA

Journal ref: Research in Astronomy and Astrophysics, Volume 23, Issue 2, id.025021, 13 pp, February 2023

First Detailed Analysis of a Relatively Deep, Low Mass-ratio Contact Binary: ATO J108.6991+27.8306

Authors: Shuo Ma, Jinzhong Liu, Yu Zhang, Guoliang Lü, Ting Wu, Chenyang He

Abstract: We present the first detailed photometric analysis of ATO J108.6991+27.8306 (hereinafter as J108). The short-period close binary J108 was observed by the Nanshan 1 m Wide Field Telescope of the Xinjiang Astronomical Observatory. The obtained BVRI-band light curves were used to determine the photometric solution by using the 2003 version of the Wilson-Devinney code. J108 is a typical deep ( f > 50%… ▽ More We present the first detailed photometric analysis of ATO J108.6991+27.8306 (hereinafter as J108). The short-period close binary J108 was observed by the Nanshan 1 m Wide Field Telescope of the Xinjiang Astronomical Observatory. The obtained BVRI-band light curves were used to determine the photometric solution by using the 2003 version of the Wilson-Devinney code. J108 is a typical deep ( f > 50%), low mass ratio (q < 0.25) overcontact binary system with a mass ratio of q = 0.1501 and a fill-out factor of f = 50.1 %, suggesting that it is in the late evolutionary stage of contact binary systems. We found the target to be a W-type W UMa binary and provided evidence for the presence of starspots on both components. From the temperature-luminosity diagram, the main component is the evolved main sequence star with an evolutionary age of about 7.94 Gyr. △ Less

Submitted 14 March, 2023; originally announced March 2023.

Comments: 7 pages, 6 figures

GECAM Localization of High Energy Transients and the Systematic Error

Authors: Yi Zhao, Wang-Chen Xue, Shao-Lin Xiong, Yuan-Hao Wang, Jia-Cong Liu, Qi Liuo, Yan-Qiu Zhang, Jian-Chao Sun, Xiao-Yun Zhao, Ce Cai, Shuo Xiao, Yue Huang, Xiao-Bo Li, Zhen Zhang, Jin-Yuan Liao, Sheng Yang, Rui Qiao, Dong-Ya Guo, Chao Zheng, Qi-Bin Yi, Sheng-Lun Xie, Zhi-Wei Guo, Chao-Yang Li, Chen-Wei Wang, Wen-Jun Tan , et al. (41 additional authors not shown)

Abstract: Gravitational wave high-energy Electromagnetic Counterpart All-sky Monitor (GECAM) is a pair of microsatellites (i.e. GECAM-A and GECAM-B) dedicated to monitoring gamma-ray transients including gravitational waves high-energy electromagnetic counterparts, Gamma-ray Bursts, Soft Gamma-ray Repeaters, Solar Flares and Terrestrial Gamma-ray Flashes. Since launch in December 2020, GECAM-B has detected… ▽ More Gravitational wave high-energy Electromagnetic Counterpart All-sky Monitor (GECAM) is a pair of microsatellites (i.e. GECAM-A and GECAM-B) dedicated to monitoring gamma-ray transients including gravitational waves high-energy electromagnetic counterparts, Gamma-ray Bursts, Soft Gamma-ray Repeaters, Solar Flares and Terrestrial Gamma-ray Flashes. Since launch in December 2020, GECAM-B has detected hundreds of astronomical and terrestrial events. For these bursts, localization is the key for burst identification and classification as well as follow-up observations in multi-wavelength. Here, we propose a Bayesian localization method with Poisson data with Gaussian background profile likelihood to localize GECAM bursts based on the burst counts distribution in detectors with different orientations. We demonstrate that this method can work well for all kinds of bursts, especially for extremely short ones. In addition, we propose a new method to estimate the systematic error of localization based on a confidence level test, which can overcome some problems of the existing method in literature. We validate this method by Monte Carlo simulations, and then apply it to a burst sample with accurate location and find that the mean value of the systematic error of GECAM-B localization is $\sim 2.5^{\circ}$. By considering this systematic error, we can obtain a reliable localization probability map for GECAM bursts. Our methods can be applied to other gamma-ray monitors. △ Less

Submitted 23 December, 2022; v1 submitted 28 November, 2022; originally announced November 2022.

Comments: The paper has been accepted by Astrophysical Journal Supplement Series

Li-rich and super Li-rich giants produced by element diffusion

Authors: Jun Gao, Chunhua Zhu, Jinlong Yu, Helei Liu, Xizhen Lu, Jianrong Shi, Guoliang Lü

Abstract: Context. About 0.2-2% of giant stars are Li-rich, whose lithium abundance (A(Li)) is higher than 1.5 dex. Among them, near 6% are super Li-rich with A(Li) exceeding 3.2 dex. Meanwhile, the formation mechanism of these Li-rich and super Li-rich giants is still under debate. Aims. Considering the compact He core of red giants, attention is paid to the effect of element diffusion on A(Li). In particu… ▽ More Context. About 0.2-2% of giant stars are Li-rich, whose lithium abundance (A(Li)) is higher than 1.5 dex. Among them, near 6% are super Li-rich with A(Li) exceeding 3.2 dex. Meanwhile, the formation mechanism of these Li-rich and super Li-rich giants is still under debate. Aims. Considering the compact He core of red giants, attention is paid to the effect of element diffusion on A(Li). In particular, when the He core flash occurs, the element diffusion makes the thermohaline mixing zone extend inward and connect to the inner convection region of stars. Then, a large amount of 7Be produced by the He flash can be transferred to stellar surface, finally turning into 7Li. Thus, the goal of this work is to propose the mechanism of A(Li) enrichment and achieve the consistency between the theoretical and observation data. Methods. Using the Modules for Experiments in Stellar Astrophysics (MESA), we simulate the evolution of low-mass stars, with considering the effects of element diffusion on the Li abundances. The timescale ratio of Li-rich giants to normal giants is estimated by population synthesis method. Then we get the theoretical value of A(Li) and make a comparison with observations. Results. Considering the influence of element diffusion in the model results in the increase of lithium abundance up to about 1.8dex, which can reveal Li-rich giants. Simultaneously, introducing high constant diffusive mixing coefficients (Dmix) with the values from 10e11 to 10e15in the model allows A(Li) to increase from 2.4 to 4.5dex, which can explain the most of Li-rich and super Li-rich giant stars. The population synthesis method reveals that the amount of Li-rich giants among giants is about 0.2-2%, which is consistent with observation estimated levels. △ Less

Submitted 24 October, 2022; originally announced October 2022.

Report number: AA/2022/43871

Journal ref: A&A 668, A126 (2022)

A Localization Method of High Energy Transients for All-Sky Gamma-Ray Monitor

Authors: Yi Zhao, Wangchen Xue, Shaolin Xiong, Qi Luo, Yuanhao Wang, Jiacong Liu, Heng Yu, Xiaoyun Zhao, Yue Huang, Jinyuan Liao, Jianchao Sun, Xiaobo Li, Qibin Yi, Ce Cai, Shuo Xiao, Shenglun Xie, Chao Zheng, Yanqiu Zhang, Chenwei Wang, Wenjun Tan, Zhiwei Guo, Chaoyang Li, Zhenghua An, Gang Chen, Yanqi Du , et al. (40 additional authors not shown)

Abstract: Fast and reliable localization of high-energy transients is crucial for characterizing the burst properties and guiding the follow-up observations. Localization based on the relative counts of different detectors has been widely used for all-sky gamma-ray monitors. There are two major methods for this counts distribution localization: $χ^{2}$ minimization method and the Bayesian method. Here we pr… ▽ More Fast and reliable localization of high-energy transients is crucial for characterizing the burst properties and guiding the follow-up observations. Localization based on the relative counts of different detectors has been widely used for all-sky gamma-ray monitors. There are two major methods for this counts distribution localization: $χ^{2}$ minimization method and the Bayesian method. Here we propose a modified Bayesian method that could take advantage of both the accuracy of the Bayesian method and the simplicity of the $χ^{2}$ method. With comprehensive simulations, we find that our Bayesian method with Poisson likelihood is generally more applicable for various bursts than $χ^{2}$ method, especially for weak bursts. We further proposed a location-spectrum iteration approach based on the Bayesian inference, which could alleviate the problems caused by the spectral difference between the burst and location templates. Our method is very suitable for scenarios with limited computation resources or time-sensitive applications, such as in-flight localization software, and low-latency localization for rapid follow-up observations. △ Less

Submitted 22 February, 2023; v1 submitted 26 September, 2022; originally announced September 2022.

Searching for Variable Stars in the Open Cluster NGC 2355 and Its Surrounding Region

Authors: Hong Wang, Yu Zhang, Xiangyun Zeng, Qingshun Hu, Jinzhong Liu, Mingfeng Qin, Guoliang Lü

Abstract: We have investigated the variable stars in the field surrounding NGC 2355 based on the time-series photometric observation data. More than 3000 CCD frames were obtained in the V band spread over 13 nights with the Nanshan One-meter Wide-field Telescope. We have detected 88 variable stars, containing 72 new variable stars and 16 known variable stars. By analyzing these light curves, we classified t… ▽ More We have investigated the variable stars in the field surrounding NGC 2355 based on the time-series photometric observation data. More than 3000 CCD frames were obtained in the V band spread over 13 nights with the Nanshan One-meter Wide-field Telescope. We have detected 88 variable stars, containing 72 new variable stars and 16 known variable stars. By analyzing these light curves, we classified the variable stars as follows: 26 eclipsing binaries, 52 pulsating stars, 4 rotating variables, and 6 unclear type variable stars for which their periods are much longer than the time baseline chosen. Employing Gaia DR2 parallax, kinematics, and photometry, the cluster membership of these variable stars were also analyzed for NGC 2355. In addition to the 11 variable members reported by Cantat-Gaudin et al. (2018), we identify 4 more variable member candidates located at the outer region of NGC 2355 and showed homogeneity in space positions and kinematic properties with the cluster members. The main physical parameters of NGC 2355 estimated from the two-color and color-magnitude diagrams are log(age/yr) = 8.9, E(B - V) = 0.24 mag, and [Fe/H] = - 0.07 dex. △ Less

Submitted 13 June, 2022; originally announced June 2022.

Comments: 15 pages, 11 figures. 6 tables,Accepted for publication in AJ

Optical Properties of C$-$rich ($^{12}$C, SiC and FeC) Dust Layered Structure of Massive Stars

Authors: Ruiqing Wu, Mengqiu Long, Xiaojiao Zhang, Yunpeng Wang, Mengli Yao, Mingming Li, Chunhua Zhu, Guoliang Lü, Zhaojun Wang, Jujia Zhang, Zhao Wang, Wujin Chen

Abstract: The composition and structure of interstellar dust are important and complex for the study of the evolution of stars and the \textbf{interstellar medium} (ISM). However, there is a lack of corresponding experimental data and model theories. By theoretical calculations based on ab-initio method, we have predicted and geometry optimized the structures of Carbon-rich (C-rich) dusts, carbon ($^{12}$C)… ▽ More The composition and structure of interstellar dust are important and complex for the study of the evolution of stars and the \textbf{interstellar medium} (ISM). However, there is a lack of corresponding experimental data and model theories. By theoretical calculations based on ab-initio method, we have predicted and geometry optimized the structures of Carbon-rich (C-rich) dusts, carbon ($^{12}$C), iron carbide (FeC), silicon carbide (SiC), even silicon ($^{28}$Si), iron ($^{56}$Fe), and investigated the optical absorption coefficients and emission coefficients of these materials in 0D (zero$-$dimensional), 1D, and 2D nanostructures. Comparing the \textbf{nebular spectra} of the supernovae (SN) with the coefficient of dust, we find that the optical absorption coefficient of the 2D $^{12}$C, $^{28}$Si, $^{56}$Fe, SiC and FeC structure corresponds to the absorption peak displayed in the infrared band (5$-$8) $μ$$m$ of the spectrum at 7554 days after the SN1987A explosion. And it also corresponds to the spectrum of 535 days after the explosion of SN2018bsz, when the wavelength in the range of (0.2$-$0.8) and (3$-$10) $μ$$m$. Nevertheless, 2D SiC and FeC corresponds to the spectrum of 844 days after the explosion of SN2010jl, when the wavelength is within (0.08$-$10) $μ$$m$. Therefore, FeC and SiC may be the second type of dust in SN1987A corresponding to infrared band (5$-$8) $μ$$m$ of dust and may be in the ejecta of SN2010jl and SN2018bsz. △ Less

Submitted 1 February, 2022; originally announced February 2022.

Report number: RAA-2021-0337.R1

A single pulse study of a millisecond pulsar PSR J0621+1002

Authors: S. Q. Wang, J. B. Wang, N. Wang, Y. Feng, S. B. Zhang, K. J. Lee, D. Li, J. G. Lu, J. T. Xie, D. J. Zhou, L. Zhang

Abstract: We present radio observation of a millisecond pulsar PSR J0621+1002 using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The pulsar shows periodic pulse intensity modulations for both the first and the third pulse components. The fluctuation spectrum of the first pulse component has one peak of 3.0$\pm$0.1 pulse periods, while that of the third pulse component has two diffused p… ▽ More We present radio observation of a millisecond pulsar PSR J0621+1002 using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The pulsar shows periodic pulse intensity modulations for both the first and the third pulse components. The fluctuation spectrum of the first pulse component has one peak of 3.0$\pm$0.1 pulse periods, while that of the third pulse component has two diffused peaks of 3.0$\pm$0.1 and 200$\pm$1 pulse periods. The single pulse timing analysis is carried out for this pulsar and the single pulses can be divided into two classes based on the post-fit timing residuals. We examined the achievable timing precision using only the pulses in one class or bright pulses. However, the timing precision improvement is not achievable. △ Less

Submitted 17 April, 2021; originally announced April 2021.

Comments: accepted for publication in ApJ. 7 pages, 8 figures

Quiescent luminosities of accreting neutron stars with different equation of states

Authors: Helei Liu, Akira Dohi, Masa-aki Hashimoto, Yasuhide Matsuo, Guo-Liang Lü, Tsuneo Noda

Abstract: We model the quiescent luminosity of accreting neutron stars with several equation of states (EOSs), including the effect of pion condensation and superfluidity. As a consequence of comparison with the observations, we show that the results with Togashi EoS (the strong direct Urca process is forbidden) and TM1e EoS (mass at direct Urca process is $2.06 M_\odot$) can explain the observations well b… ▽ More We model the quiescent luminosity of accreting neutron stars with several equation of states (EOSs), including the effect of pion condensation and superfluidity. As a consequence of comparison with the observations, we show that the results with Togashi EoS (the strong direct Urca process is forbidden) and TM1e EoS (mass at direct Urca process is $2.06 M_\odot$) can explain the observations well by considering pion condensation and the effect of superfluidity, while LS220 EoS and TM1 EoS can explain the observations well by considering the baryon direct Urca process and the effect of superfluidity. Besides, we compare the results with the observations of a neutron star RX J0812.4-3114 which has the low average mass accretion rate ($\langle\dot{M}\rangle\sim(4-15)\times 10^{-12}~M_\odot ~\rm yr^{-1}$) but high thermal luminosity ($L_q^\infty\sim(0.6-3)\times10^{33}~\rm erg ~ s^{-1}$), and we suggest that a low-mass neutron star ($<1M_\odot$) with minimum cooling can explain the lower limit of the observation of thermal luminosity of RX J0812.4-3114, which is qualitatively consistent with the previous work~\cite{Zhao2019}. However, to explain its upper limit, some other heating mechanisms besides standard deep crustal heating may be needed. △ Less

Submitted 8 February, 2021; originally announced February 2021.

Comments: 12 pages, 8 figures, accepted by Phys. Rev. D

Report number: RIKEN-iTHEMS-Report-21

Heavy element contributions of rotating massive stars to Interstellar Medium

Authors: R. Q. Wu, C. H. Zhu, G. L. Lü, Z. J. Wang, H. L. Liu

Abstract: Employing the the stellar evolution code (Modules for Experiments in Stellar Astrophysics), we calculate yields of heavy elements from massive stars via stellar wind and core-collapse supernovae (CCSN) ejecta to interstellar medium (ISM). In our models, the initial masses ($M_{\rm ini}$) of massive stars are taken from 13 to 80 $M_\odot$, their initial rotational velocities (V) are 0, 300 and 500… ▽ More Employing the the stellar evolution code (Modules for Experiments in Stellar Astrophysics), we calculate yields of heavy elements from massive stars via stellar wind and core-collapse supernovae (CCSN) ejecta to interstellar medium (ISM). In our models, the initial masses ($M_{\rm ini}$) of massive stars are taken from 13 to 80 $M_\odot$, their initial rotational velocities (V) are 0, 300 and 500 km s$^{-1}$, and their metallicities are [Fe/H] = -3, -2, -1, and 0. The yields of heavy elements coming from stellar winds are mainly affected by the stellar rotation which changes the chemical abundances of stellar surfaces via chemically homogeneous evolution, and enhances mass-loss rate. We estimate that the stellar wind can produce heavy element yields of about $10^{-2}$ (for low metallicity models) to several $M_\odot$ (for low metallicity and rapid rotation models) mass. The yields of heavy element produced by CCSN ejecta also depend on the remnant mass of massive mass which is mainly determined by the mass of CO-core. Our models calculate that the yields of heavy elements produced by CCSN ejecta can get up to several $M_\odot$. Compared with stellar wind, CCSN ejecta has a greater contribution to the heavy elements in ISM. We also compare the $^{56}$Ni yields by calculated in this work with observational estimate. Our models only explain the $^{56}$Ni masses produced by faint SNe or normal SNe with progenitor mass lower than about 25 $M_\odot$, and greatly underestimate the $^{56}$Ni masses produced by stars with masses higher than about 30 $M_\odot$. △ Less

Submitted 2 February, 2021; originally announced February 2021.

Comments: 11 pages,9 figures

Report number: RAA-2020-0450.R1

Periodic and Phase-locked Modulation in PSR B1929+10 Observed with FAST

Authors: F. F. Kou, W. M. Yan, B. Peng, J. G. Lu, K. Liu, C. M. Zhang, R. G. Strom, L. Wang, J. P. Yuan, Rai Yuen, Y. Z. Yu, J. M. Yao, B. Liu, J. Yan, P. Jiang, C. J. Jin, D. Li, L. Qian, Y. L. Yue, Y. Zhu

Abstract: We present a detailed single-pulse analysis for PSR B1929+10 based on observations with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The main pulse and interpulse are found to be modulated with a periodicity of $\sim12$ times the pulsar's rotational period ($P$). The $\sim12P$ modulation is confirmed as a periodic amplitude modulation instead of systematic drifting. The period… ▽ More We present a detailed single-pulse analysis for PSR B1929+10 based on observations with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The main pulse and interpulse are found to be modulated with a periodicity of $\sim12$ times the pulsar's rotational period ($P$). The $\sim12P$ modulation is confirmed as a periodic amplitude modulation instead of systematic drifting. The periodic amplitude modulation in the IP is found to be anti-correlated with that in the weak preceding component of the MP (MP_I), but correlated with that in the first two components of the MP (MP_II), which implies that the modulation patterns in the IP and the MP are phase-locked. What is more interesting is that the modulation in MP_II is delayed that in the IP by about 1P. Furthermore, high sensitivity observations by FAST reveal that weak emission exists between the MP and the IP. In addition, we confirm that the separation between the IP and the MP is independent of radio frequency. The above results are a conundrum for pulsar theories and cannot be satisfactorily explained by the current pulsar models. Therefore, our results observed with FAST provide an opportunity to probe the structure of pulsar emission and the neutron star's magnetosphere. △ Less

Submitted 18 December, 2020; originally announced December 2020.

Comments: 13 pages 12 figures

Formation and Evolution of Ultraluminous X-ray Pulsar Binaries to Pulsar-Neutron Star and Pulsar-White Dwarf Systems

Authors: Kunduz Abdusalam, Iminhaji Ablimit, P. Hashim, G. -L Lü, M. K. Mardini, Z. -J Wang

Abstract: Recent observational and theoretical results have suggested that some of ultraluminous X-ray (ULX) sources may contain neutron star (NS) accretors. However, the formation channel and properties of donor stars of NS ULXs remain uncertain. By adopting the non-conservative and rotation-dependent mass transfer model in the primordial binary evolution, we investigate the way to form pulsar ULXs like ob… ▽ More Recent observational and theoretical results have suggested that some of ultraluminous X-ray (ULX) sources may contain neutron star (NS) accretors. However, the formation channel and properties of donor stars of NS ULXs remain uncertain. By adopting the non-conservative and rotation-dependent mass transfer model in the primordial binary evolution, we investigate the way to form pulsar ULXs like observed pulsar ULXs in a systematic way. Our simulation results indicate that pulsar ULXs with Be stars and intermediate or/and high mass donors match observed apparent luminosities, orbital periods and observationally indicated donor masses of known pulsar ULXs. ULXs with Be and intermediate donors are main contributors. The route of accretion-induced collapse of WDs has 4.5% contribution to the NS ULXs, 4.0% of NSs in ULXs are formed through electron-capture supernovae (SNe), and 91.5% of NSs in ULXs are born with core-collapse SNe. We also studied the evolution of pulsar ULXs to double compact star systems. We do not find NS-black hole systems (merging in a Hubble time) that evolved from pulsar ULXs. Pulsar-white dwarf (WD) cases that evolve through pulsar ULXs have significant contributions to the whole NS-WD gravitational wave sources. Contributions of pulsar-WD and pulsar-NS cases that experienced pulsar ULXs are $\sim$40% and 11% among all LISA NS-WD and NS-NS sources, respectively. Monte Carlo simulation noise with different models give a non-negligible uncertainty. △ Less

Submitted 23 October, 2020; v1 submitted 28 September, 2020; originally announced September 2020.

Comments: 6 Figures, 2 tables; Published in ApJ (new version)

Journal ref: The Astrophysical Journal, 902:125 (10pp), 2020 October 20

Possible Formation Scenarios of ZTF J153932.16+502738.8-A Gravitational Source Close to the Peak of LISA's Sensitivity

Authors: Guoliang Lü, Chunhua Zhu, Zhaojun Wang, Helei Liu, Lin Li, Dian Xie, Jinzhong Liu

Abstract: ZTF J153932.16+502738.8 (ZTFJ1539) is an eclipsing double-white-dwarf system with an orbital period of 6.91 minutes, and is a significant source of LISA detecting gravitational wave. However, the massive white dwarf (WD) with mass of about 0.61 M$_\odot$ has a high effective temperature (48900 K), and the lower mass WD with mass of about 0.21 M$_{\odot}$ has a low effective temperature($<$10000 K)… ▽ More ZTF J153932.16+502738.8 (ZTFJ1539) is an eclipsing double-white-dwarf system with an orbital period of 6.91 minutes, and is a significant source of LISA detecting gravitational wave. However, the massive white dwarf (WD) with mass of about 0.61 M$_\odot$ has a high effective temperature (48900 K), and the lower mass WD with mass of about 0.21 M$_{\odot}$ has a low effective temperature($<$10000 K). It is challenging the popular theory of binary evolution. We investigate the formation of ZTFJ1539 via nova and Algol scenarios. Assuming that the massive WD in ZTFJ1539 just experiences a thermalnuclear runaway, nova scenario can explain the effective temperatures of two WDs in ZTFJ1539. However, in order to enlarging a semi-detached orbit of about 4---5 minutes to a detached orbit of about 7 minutes, nova scenario needs a much high kick velocity of about 200 km s$^{-1}$ during nova eruption. The high kick velocity can result in high eccentricity of about 0.2---0.6. Algol scenario can also produce ZTFJ1539 if we take a high efficient parameter for ejecting common envelope and enhance the mass-loss rate via stellar wind trigger by tidal effect. △ Less

Submitted 6 January, 2020; originally announced January 2020.

Comments: 10 pages, 8 figures, Accepted for publication in ApJ

The formation of subdwarf A-type stars

Authors: Jinlong Yu, Zhenwei Li, Chunhua Zhu, Zhaojun Wang, Helei Liu, Sufen Guo, Zhanwen Han, Xuefei Chen, Guoliang Lü

Abstract: Subdwarf A-type stars (sdAs) are objects that have hydrogen-rich spectra with surface gravity similar to that of hot subdwarf stars but effective temperature below the zero-age horizontal branch (ZAHB). They are considered to be metal-poor main sequence (MS) stars or extremely low-mass white dwarfs (ELM WDs). In this work, using the stellar evolution code, Modules for Experiments in Stellar Astrop… ▽ More Subdwarf A-type stars (sdAs) are objects that have hydrogen-rich spectra with surface gravity similar to that of hot subdwarf stars but effective temperature below the zero-age horizontal branch (ZAHB). They are considered to be metal-poor main sequence (MS) stars or extremely low-mass white dwarfs (ELM WDs). In this work, using the stellar evolution code, Modules for Experiments in Stellar Astrophysics (MESA), we investigate the sdAs formed both by the evolution of (pre-)ELM WDs in double degenerate systems (DDs) and metal-poor main sequence with the single evolution models. We find that both of the evolutionary tracks of ELM WDs and metal-poor MSs can explain the observation properties of sdAs. However, the proportion between these two populations are uncertain. In this work, we adopt the method of binary population synthesis of both ELM WDs in the disk and metal-poor MS in the halo to obtain the populations of ELM WDs and metal-poor MSs at different stellar population ages and calculate their proportions. We find that the proportion of metal-poor MSs to sdA population for a stellar population of 10 Gyr is $\sim$98.5 percent, which is consistent with the conclusion that most sdAs ($>$95 percent) are metal-poor MSs. And the proportion of ELM WDs (metal-poor MSs) to sdA population increases (decreases) from 0.1\% (99.9\%) to 20\% (80\%) with the stellar population ages from 5 to 13.7 Gyr. △ Less

Submitted 17 September, 2019; originally announced September 2019.

Comments: 9 pages, 6 figures, 1 table,Accepter for publication in ApJ

Quark Novae: An Alternative Channel For the Formation of Isolated Millisecond Pulsars

Authors: Nurimangul Nurmamat, Chunhua Zhu, Guoliang Lü, Zhaojun Wang, Lin LI, Helei Liu

Abstract: Isolated millisecond pulsars (IMSPs) are a topic of academic contention. There are various models to explain their formation. We explore the formation of IMSP via quark novae (QN). During this formation process, low-mass X-ray binaries (LMXBs) are disrupted when the mass of the neutron star (NS) reaches 1.8 $\rm M_\odot$. Using population synthesis, this work estimates that the Galactic birthrate… ▽ More Isolated millisecond pulsars (IMSPs) are a topic of academic contention. There are various models to explain their formation. We explore the formation of IMSP via quark novae (QN). During this formation process, low-mass X-ray binaries (LMXBs) are disrupted when the mass of the neutron star (NS) reaches 1.8 $\rm M_\odot$. Using population synthesis, this work estimates that the Galactic birthrate of QN-produced IMSPs lies between $\sim 9.5\times10^{-6}$ and $\sim1.7\times10^{-4}$ $\rm yr^{-1}$. The uncertainties shown in our experiment model is due to the QN's kick velocity. Furthermore, our findings not only show that QN-produced IMSPs are statistically more significant than those produced by mergers, but also that millisecond pulsar binaries with a high eccentricity may originate from LMXBs that have been involved in, yet not disrupted by, a QN. △ Less

Submitted 19 June, 2019; originally announced June 2019.

Comments: 7 pages, 2 figures, accepted for publication in the Journal of Astrophysics and Astronomy

The Formation of Bimodal Dust Species in Nova Ejecta I: Chemical Conditions

Authors: Chunhua Zhu, Helei Liu, Guoliang Lü, Zhaojun Wang, Lin Li

Abstract: The formation of bimodal dust species (that is, both of silicate and amorphous carbon dust grains are observed in a nova eruption) in nova ejecta is still debated. Using the Modules for Experiments in Stellar Astrophysics code and considering the effects of WD's mass, mass-accretion rate and the chemical profiles of WD which are described by new parameter --- mixing depth on the chemical abundance… ▽ More The formation of bimodal dust species (that is, both of silicate and amorphous carbon dust grains are observed in a nova eruption) in nova ejecta is still debated. Using the Modules for Experiments in Stellar Astrophysics code and considering the effects of WD's mass, mass-accretion rate and the chemical profiles of WD which are described by new parameter --- mixing depth on the chemical abundances of nova ejecta, we investigate the possibility that bimodal dust species are produced in a nova eruption. We find that $C/O$ (the ratio of carbon number density to oxygen number density) of nova ejecta is affected by the mixing depth. For the model with a small mixing depth, the $C/O$ of nova ejecta can evolute from lager than 1.0 to smaller than 1.0 in a whole eruption, which provides the chemical condition for the formation of bimodal dust species. △ Less

Submitted 10 June, 2019; originally announced June 2019.

Comments: 13 pages, 7 figures, 4 Tables, accepted for publication in MNRAS

On Isolated Millisecond Pulsars Formed by the Coalescence of Neutron Stars and Massive White Dwarfs

Authors: Shengnan Sun, Lin Li, Helei Liu, Guoliang Lü, Zhaojun Wang, Chunhua Zhu

Abstract: This paper uses population synthesis to investigate the possible origin of isolated millisecond pulsars (IMSPs) as born from the coalescence of a neutron star (NS) and a white dwarf (WD). Results show that the galactic birth - rate of IMSPs is likely to lie between $5.8\times10^{-5}$ $\rm yr^{-1}$ and $2.0\times10^{-4}$ $\rm yr^{-1}$, depending on critical variables, such as the stability of mass… ▽ More This paper uses population synthesis to investigate the possible origin of isolated millisecond pulsars (IMSPs) as born from the coalescence of a neutron star (NS) and a white dwarf (WD). Results show that the galactic birth - rate of IMSPs is likely to lie between $5.8\times10^{-5}$ $\rm yr^{-1}$ and $2.0\times10^{-4}$ $\rm yr^{-1}$, depending on critical variables, such as the stability of mass transfer via the Roche lobe and the value of kick velocity. In addition to this, this paper estimates that the solar mass of IMSPs can range from $1.5$ and 2.0 $M_{\odot}$, making them more massive than other 'normal' pulsars. Finally, the majority of IMSPs in our simulations have spin periods ranging from several to 20 milliseconds, which is consistent with previous observations. △ Less

Submitted 20 December, 2018; v1 submitted 3 December, 2018; originally announced December 2018.

Comments: 9 pages, 3 figures, 1 table, Accepted for publication in PASA

Magnetic-oscillation mechanism for understanding periodic or quasi-periodic modulation of the timing residuals from pulsars

Authors: Zhaojun Wang, Guoliang Lü, Chunhua Zhu, Lin Li, Anzhong Wang

Abstract: Highly precise pulsar timing is very important for understanding the nature of neutron stars and the implied physics, even being used to detect directly gravitational waves. Unfortunately, the accuracy of the pulsar timing is seriously affected by the spin-down irregularities, such as spin fluctuations with a manifestation of low frequency structures (the so-called red noise processes), and variou… ▽ More Highly precise pulsar timing is very important for understanding the nature of neutron stars and the implied physics, even being used to detect directly gravitational waves. Unfortunately, the accuracy of the pulsar timing is seriously affected by the spin-down irregularities, such as spin fluctuations with a manifestation of low frequency structures (the so-called red noise processes), and various activities of magnetospheres. Except from the random timing noise, significant periodic or quasi-periodic components could also be found in some timing residuals of pulsars, indicating the presence of unmodelled deterministic effects. The physical origins of these effects still remain unexplained. In this paper, we suggest a new mechanism involving the de Haas-van Alphen magnetic oscillation, which could trigger the observed low frequency structures. We find that the proposed magnetic oscillation period is about 1-$10^{2}$ yr, which is about $10^{-4}$ times as long as the classical characteristic time scale of interior magnetic field evolution for a normal neutron star. Due to the magnetic oscillation, we estimate both the braking index range between $10^{-5}$ and $10^{5}$ and the residuals range between 4 to 906 ms for some specific samples of quasi-periodic pulsars. Those are consistent with the pulsar timing observations. Nonlinear phenomena of de Haas-van Alphen effect and the Condon domain structure are considered. Avalanche-like magnetized process of domain wall motion (rapidly magnetic energy release) could associated with some special emitting manifestations of neuron star radiation. △ Less

Submitted 8 November, 2018; originally announced November 2018.

Comments: 10 pages, 1 figure, submitted. arXiv admin note: text overlap with arXiv:1801.01832

The Giant Radio Array for Neutrino Detection (GRAND): Science and Design

Authors: GRAND Collaboration, Jaime Alvarez-Muniz, Rafael Alves Batista, Aswathi Balagopal V., Julien Bolmont, Mauricio Bustamante, Washington R. Carvalho, Didier Charrier, Ismael Cognard, Valentin Decoene, Peter B. Denton, Sijbrand J. De Jong, Krijn D. De Vries, Ralph Engel, Ke Fang, Chad Finley, Stefano Gabici, QuanBu Gou, Junhua Gu, Claire Guépin, Hongbo Hu, Yan Huang, Kumiko Kotera, Sandra Le Coz, Jean-Philippe Lenain , et al. (26 additional authors not shown)

Abstract: 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 th… ▽ More 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. △ Less

Submitted 18 July, 2019; v1 submitted 23 October, 2018; originally announced October 2018.

Comments: 45 pages, 34 figures. Matches published version

Journal ref: Sci. China-Phys. Mech. Astron. 63, 219501 (2020)

An Alternative Channel For High-Mass Binary Black Holes-----Dark Matter Accretion Onto Black Hole

Authors: Tongzheng Wang, Lin Li, Chunhua Zhu, Zhaojun Wang, Anzhong Wang, Qiang Wu, Hai-Shan Liu, Guoliang Lü

Abstract: By a method of population synthesis we construct a model of dark matter (DM) accretion onto binary black holes (BHs), and investigate the merger rate of the binary BHs. We find that the merger rate can weakly increase (less than 10\%). However, the DM accretion can efficiently enhance the masses of binary BHs. In our model, the result for $Z = 0.01$ without the DM accretion can not explain the GW1… ▽ More By a method of population synthesis we construct a model of dark matter (DM) accretion onto binary black holes (BHs), and investigate the merger rate of the binary BHs. We find that the merger rate can weakly increase (less than 10\%). However, the DM accretion can efficiently enhance the masses of binary BHs. In our model, the result for $Z = 0.01$ without the DM accretion can not explain the GW170104, GW170814, and GW150914, while with the DM accretion, it can cover all observations well. For the higher metallicity ($Z = 0.02$), our model can not explain the mergers of high-mass binary BHs like GW170104, GW170814, and GW150914. We estimate that the merger rate of binary BHs lies between $\rm55Gpc^{-3}yr^{-1}$ to $\rm197Gpc^{-3}yr^{-1}$. △ Less

Submitted 26 June, 2018; originally announced June 2018.

Comments: 8 pages, 5 figures, accepted for the ApJ

A Study of Wolf-Rayet stars Formed via Chemically Homogeneous Evolution

Authors: Zhe Cui, Zhaojun Wang, Chunhua Zhu, Guoliang Lu, Hailiang Chen, Zhanwen Han

Abstract: 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 o… ▽ More 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). △ Less

Submitted 22 May, 2018; originally announced May 2018.

Comments: 13 pages, 4 figures, Accepted for publication in PASP

Magnetic-oscillation Mechanism for Understanding Timing Noise of Pulsar Observation

Authors: Zhaojun Wang, Guoliang Lü, Chunhua Zhu, Lin Li, Anzhong Wang

Abstract: Highly precise pulsar timing is very important for understanding the nature of a neutron star, and it can even be used to detect gravitational waves. Unfortunately, the accuracy of the pulsar timing is seriously affected by the spin-down irregularities of pulsars. Highly precise pulsar timing is very important for understanding the nature of a neutron star, and it can even be used to detect gravitational waves. Unfortunately, the accuracy of the pulsar timing is seriously affected by the spin-down irregularities of pulsars. △ Less

Submitted 12 November, 2018; v1 submitted 4 January, 2018; originally announced January 2018.

Comments: The manuscript is not accepted by any journal, and there are some mistakes in it

Effects of Irradiation on the Evolution of Ultra-compact X-ray Binaries

Authors: Guoliang Lü, Chunhua Zhu, Zhaojun Wang, Hoernisa Iminniyaz

Abstract: Using the Modules for Experiments in Stellar Astrophysics code, we investigate the influences of irradiation on ultra-compact X-ray binary (UCXB) evolution. Although the persistent UCXBs have short orbital periods which result in high irradiation flux, the irradiation hardly affects the evolution of persistent sources because the WDs in these binaries have large masses which lead to very low irrad… ▽ More Using the Modules for Experiments in Stellar Astrophysics code, we investigate the influences of irradiation on ultra-compact X-ray binary (UCXB) evolution. Although the persistent UCXBs have short orbital periods which result in high irradiation flux, the irradiation hardly affects the evolution of persistent sources because the WDs in these binaries have large masses which lead to very low irradiation depth. The irradiation has a significant effect on the transient sources during outburst phase. At the beginning of the outburst, high X-ray luminosity produces high radiation flux, which results in the significant expansion of WD. Then, the irradiation triggers high mass-transfer rates, which can last several days for the transient sources with WDs whose masses are larger than $\sim0.015 M_\odot$ or several hundred years for these sources with WDs whose masses are less than $\sim0.012 M_\odot$. The observed three persistent UCXBs, XTE J0929-314, 4U 1916-05 and SWIFT J1756.9-2508, may belong to the latter. △ Less

Submitted 4 September, 2017; originally announced September 2017.

Comments: 11 pages, 7 figures, 1 table, accepted for publication in ApJ

Pulsar Glitches in a Strangeon Star Model

Authors: X. Y. Lai, C. A. Yun, J. G. Lu, G. L. Lv, Z. J. Wang, R. X. Xu

Abstract: Pulsar-like compact stars provide us a unique laboratory to explore properties of dense matter at supra-nuclear densities. One of the models for pulsar-like stars is that they are totally composed of "strangeons", and in this paper we studied the pulsar glitches in a strangeon star model. Strangeon stars would be solidified during cooling, and the solid stars would be natural to have glitches as t… ▽ More Pulsar-like compact stars provide us a unique laboratory to explore properties of dense matter at supra-nuclear densities. One of the models for pulsar-like stars is that they are totally composed of "strangeons", and in this paper we studied the pulsar glitches in a strangeon star model. Strangeon stars would be solidified during cooling, and the solid stars would be natural to have glitches as the result of starquakes. Based on the starquake model established before, we proposed that when the starquake occurs, the inner motion of the star which changes the moment of inertia and has impact on the glitch sizes, is divided into plastic flow and elastic motion. The plastic flow which is induced in the fractured part of the outer layer, would move tangentially to redistribute the matter of the star and would be hard to recover. The elastic motion, on the other hand, changes its shape and would recover significantly. Under this scenario, we could understand the behaviors of glitches without significant energy releasing, including the Crab and the Vela pulsars, in an uniform model. We derive the recovery coefficient as a function of glitch size, as well as the time interval between two successive glitches as the function of the released stress. Our results show consistency with observational data under reasonable ranges of parameters. The implications on the oblateness of the Crab and the Vela pulsars are discussed. △ Less

Submitted 12 March, 2018; v1 submitted 24 July, 2017; originally announced July 2017.

Comments: MNRAS, accepted

Investigating multi-frequency pulse profiles of PSRs B0329+54 and B1642-03 in an inverse Compton scattering (ICS) model

Authors: L. H. Shang, J. G. Lu, Y. J. Du, L. F. Hao, D. Li, K. J. Lee, Bin Li, L. X. Li, G. J. Qiao, Z. Q. Shen, D. H. Wang, M. Wang, X. J. Wu, Y. J. Wu, R. X. Xu, Y. L. Yue, Z. Yan, Q. J. Zhi, R. B. Zhao, R. S. Zhao

Abstract: The emission geometries, e.g. the emission region height, the beam shape, and radius-to-frequency mapping, are important predictions of pulsar radiation model. The multi-band radio observations carry such valuable information. In this paper, we study two bright pulsars, (PSRs B0329+54 and B1642-03) and observe them in high frequency (2.5 GHz, 5 GHz, and 8 GHz). The newly acquired data together wit… ▽ More The emission geometries, e.g. the emission region height, the beam shape, and radius-to-frequency mapping, are important predictions of pulsar radiation model. The multi-band radio observations carry such valuable information. In this paper, we study two bright pulsars, (PSRs B0329+54 and B1642-03) and observe them in high frequency (2.5 GHz, 5 GHz, and 8 GHz). The newly acquired data together with historical archive provide an atlas of multi-frequency profiles spanning from 100 MHz to 10 GHz. We study the frequency evolution of pulse profiles and the radiation regions with the these data. We firstly fit the pulse profiles with Gaussian functions to determine the phase of each component, and then calculate the radiation altitudes of different emission components and the radiation regions. We find that the inverse Compton scattering (ICS) model can reproduce the radiation geometry of these two pulsars. But for PSR B0329+54 the radiation can be generated in either annular gap (AG) or core gap (CG), while the radiation of PSR B1642-03 can only be generated in the CG. This difference is caused by the inclination angle and the impact angle of these two pulsars. The relation of beaming angle (the angle between the radiation direction and the magnetic axis) and the radiation altitudes versus frequency is also presented by modelling the beam-frequency evolution in the ICS model. The multi-band pulse profiles of these two pulsars can be described well by the ICS model combined with the CG and AG. △ Less

Submitted 10 March, 2017; originally announced March 2017.

Effects of the core-collapse supernova ejecta impact on a rapidly rotating massive companion star

Authors: Chunhua Zhu, Guoliang Lu, Zhaojun Wang

Abstract: 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… ▽ More 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. △ Less

Submitted 3 January, 2017; originally announced January 2017.

Comments: 11 pages, 10 figures, accepted for publication in ApJ

The Contribution of Chemical Abundances in Nova Ejecta to the Interstellar Medium

Authors: Fanger Li, Chunhua Zhu, Guoliang Lu, Zhaojun Wang

Abstract: 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,… ▽ More 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. △ Less

Submitted 30 May, 2016; originally announced May 2016.

Comments: 15 pages, 5 fiures, 4 tables, Accepted for publication in PASJ

The diamagnetic phase transition in Magnetars

Authors: Zhaojun Wang, Guoliang Lu, Chunhua Zhu, Baoshan Wu

Abstract: 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… ▽ More 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. △ Less

Submitted 30 May, 2016; originally announced May 2016.

Comments: 15pages, 4 fiures, accepted for publication in PASP

Population Synthesis of Millisecond X-ray Pulsars

Authors: Chunhua Zhu, Guoliang Lu, Zhaojun Wang

Abstract: 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)… ▽ More 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. △ Less

Submitted 15 September, 2015; originally announced September 2015.

Comments: 12 pages, 7 figures, 3 tables, Accepted for publication in MNRAS

Monte Carlo Simulation of HERD Calorimeter

Authors: M. Xu, G. M. Chen, Y. W. Dong, J. G. Lu, Z. Quan, L. Wang, Z. G. Wang, B. B. Wu, S. N. Zhang

Abstract: The High Energy cosmic-Radiation Detection (HERD) facility onboard China's Space Station is planned for operation starting around 2020 for about 10 years. It is designed as a next generation space facility focused on indirect dark matter search, precise cosmic ray spectrum and composition measurements up to the knee energy, and high energy gamma-ray monitoring and survey. The calorimeter plays an… ▽ More The High Energy cosmic-Radiation Detection (HERD) facility onboard China's Space Station is planned for operation starting around 2020 for about 10 years. It is designed as a next generation space facility focused on indirect dark matter search, precise cosmic ray spectrum and composition measurements up to the knee energy, and high energy gamma-ray monitoring and survey. The calorimeter plays an essential role in the main scientific objectives of HERD. A 3-D cubic calorimeter filled with high granularity crystals as active material is a very promising choice for the calorimeter. HERD is mainly composed of a 3-D calorimeter (CALO) surrounded by silicon trackers (TK) from all five sides except the bottom. CALO is made of 9261 cubes of LYSO crystals, corresponding to about 55 radiation lengths and 3 nuclear interaction lengths, respectively. Here the simulation results of the performance of CALO with GEANT4 and FLUKA are presented: 1) the total absorption CALO and its absorption depth for precise energy measurements (energy resolution: 1\% for electrons and gamma-rays beyond 100 GeV, 20\% for protons from 100 GeV to 1 PeV); 2) its granularity for particle identification (electron/proton separation power better than $10^{-5}$); 3) the homogenous geometry for detecting particles arriving from every unblocked direction for large effective geometrical factor ($>$3 ${\rm m}^{2}{\rm sr}$ for electron and diffuse gamma-rays, $>$2 $ {\rm m}^{2}{\rm sr}$ for cosmic ray nuclei); 4) expected observational results such as gamma-ray line spectrum from dark matter annihilation and spectrum measurement of various cosmic ray chemical components. △ Less

Submitted 16 July, 2014; originally announced July 2014.

Journal ref: Proceedings of the SPIE, 9144-145, 2014

New Vacuum Solar Telescope and Observations with High Resolution

Authors: Zhong Liu, Jun Xu, Bo-Zhong Gu, Sen Wang, Jian-Qi You, Long-Xiang Shen, Ru-Wei Lu, Zhen-Yu Jin, Lin-Fei Chen, Ke Lou, Zhi Li, Guang-Qian Lu, Zhi Xu, Chang-Hui Rao, Qi-Qian Hu, Ru-Feng Li, Hao-Wen Fu, Men-Xian Bao, Ming-Chan Wu, Bo-Rong Zhang

Abstract: 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,… ▽ More 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. △ Less

Submitted 26 March, 2014; originally announced March 2014.

Comments: 16 pages, 14 figures, accepted by RAA (Research in Astronomy and Astrophysics)

The nature of [S III]λλ9096, 9532 emitters at z = 1.34 and 1.23

Authors: F. X. An, X. Z. Zheng, Y. Z. Meng, Y. Chen, Z. Z. Wen, G. L. Lü

Abstract: A study of [S III]$λ\lambda9096,9532$ emitters at $z$ = 1.34 and 1.23 is presented using our deep narrow-band $H_2S1$ (centered at 2.13 $μ$m) imaging survey of the Extended Chandra Deep Field South (ECDFS). We combine our data with multi-wavelength data of ECDFS to build up spectral energy distributions (SEDs) from the $U$ to the $K_{s}$-band for emitter candidates selected with strong excess in… ▽ More A study of [S III]$λ\lambda9096,9532$ emitters at $z$ = 1.34 and 1.23 is presented using our deep narrow-band $H_2S1$ (centered at 2.13 $μ$m) imaging survey of the Extended Chandra Deep Field South (ECDFS). We combine our data with multi-wavelength data of ECDFS to build up spectral energy distributions (SEDs) from the $U$ to the $K_{s}$-band for emitter candidates selected with strong excess in $H_2S1 - K_{s}$ and derive photometric redshifts, line luminosities, stellar masses and extinction. A sample of 14 [S III] emitters are identified with $H_2S1<22.8$ and $K_{\rm s}<24.8$ (AB) over 381 arcmin$^{2}$ area, having [S III] line luminosity $L_{[SIII]}= \sim 10^{41.5-42.6}$erg s$^{-1}$. None of the [S III] emitters is found to have X-ray counterpart in the deepest Chandra 4 Ms observation, suggesting that they are unlikely powered by AGN. HST/ACS F606W and HST/WFC3 F160W images show their rest-frame UV and optical morphologies. About half of the [S III] emitters are mergers and at least one third are disk-type galaxies. Nearly all [S III] emitters exhibit a prominent Balmer break in their SEDs, indicating the presence of a significant post-starburst component. Taken together, our results imply that both shock heating in post-starburst and photoionization caused by young massive stars are likely to excite strong [S III] emission lines. We conclude that the emitters in our sample are dominated by star-forming galaxies with stellar mass $8.7<\log (M/M_{sun})<9.9$. △ Less

Submitted 14 February, 2014; originally announced February 2014.

Comments: 9 pages, 3 figures, 4 tables, published in SCIENCE CHINA Physics, Mechanics & Astronomy

Journal ref: SCIENCE CHINA Physics, Mechanics & Astronomy, 2013, 56(11): 2226-2235

Origin of Dust around V1309 Sco

Authors: Chunhua Zhu, Guoliang Lu, Zhaojun Wang

Abstract: 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… ▽ More 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. △ Less

Submitted 26 August, 2013; originally announced August 2013.

Comments: 7 pages, 4 figures, Accepted for publication in the Astrophysical Journal

Dust Formation in the Ejecta of Common Envelope Systems

Authors: Guoliang Lu, Chunhua Zhu, Philipp Podsiadlowski

Abstract: 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 t… ▽ More 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. △ Less

Submitted 27 March, 2013; originally announced March 2013.

Comments: 11pages, 7 figures, accepted for publication by ApJ

Common Envelope Evolution: Where we stand and how we can move forward

Authors: N. Ivanova, S. Justham, X. Chen, O. De Marco, C. L. Fryer, E. Gaburov, H. Ge, E. Glebbeek, Z. Han, X. -D. Li, G. Lu, T. Marsh, Ph. Podsiadlowski, A. Potter, N. Soker, R. Taam, T. M. Tauris, E. P. J. van den Heuvel, R. F. Webbink

Abstract: 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 observation… ▽ More 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. △ Less

Submitted 6 December, 2012; v1 submitted 19 September, 2012; originally announced September 2012.

Comments: 83 pages (11 chapters, 2 appendixes), 10 figures, Astronomy & Astrophysics Review format. This version contains revisions as a response to community comments

Population Synthesis for Symbiotic X-ray Binaries

Authors: G. -L. Lu, C. -H. Zhu, K. A. Postnov, L. R. Yungelson, A. G. Kuranov, N. Wang

Abstract: 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 obs… ▽ More 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. △ Less

Submitted 25 May, 2012; originally announced May 2012.

Comments: 12 pages, 6 figures, 3 tables (accepted for publication in MNRAS)

Wind accretion in symbiotic X-ray binaries

Authors: K. Postnov, N. Shakura, A. González-Galán, E. Kuulkers, P. Kretschmar, S. Larsson, M. H. Finger, A. Kochetkova, G. Lü, L. Yungelson

Abstract: The properties of wind accretion in symbiotic X-ray binaries (SyXBs) consisting of red-giant and magnetized neutron star (NS) are discussed. The spin-up/spin-down torques applied to NS are derived based on a hydrodynamic theory of quasi-spherical accretion onto magnetized NSs. In this model, a settling subsonic accretion proceeds through a hot shell formed around the NS magnetosphere. The accretio… ▽ More The properties of wind accretion in symbiotic X-ray binaries (SyXBs) consisting of red-giant and magnetized neutron star (NS) are discussed. The spin-up/spin-down torques applied to NS are derived based on a hydrodynamic theory of quasi-spherical accretion onto magnetized NSs. In this model, a settling subsonic accretion proceeds through a hot shell formed around the NS magnetosphere. The accretion rate onto the NS is determined by the ability of the plasma to enter the magnetosphere.Due to large Reynolds numbers in the shell, the interaction of the rotating magnetosphere with plasma initiates a subsonic turbulence. The convective motions are capable of carrying the angular momentum through the shell. We carry out a population synthesis of SyXBs in the Galaxy with account for the spin evolution of magnetized NS. The Galactic number of SyXBs with bright (M_v<1) low-mass red-giant companion is found to be from \sim 40 to 120, and their birthrate is \sim 5\times 10^{-5}-10^{-4} per year. According to our model, among known SyXBs, Sct X-1 and IRXS J180431.1-273932 are wind-fed accretors. GX 1+4 lies in the transition from the wind-fed SyXBs to SyXBs in which the giants overflow their Roche lobe. The model successfully reproduces very long NS spins (such as in IGR J16358-4724 and 4U 1954+31) without the need to invoke very strong magnetic fields. △ Less

Submitted 6 October, 2011; originally announced October 2011.

Comments: 7 pages, 2 figures, Proceedings of the 8th INTEGRAL Workshop, Dublin; http://pos.sissa.it/archive/conferences/115/015/INTEGRAL%202010_015.pdf

Journal ref: PoS (INTEGRAL 2010)015

Gamma-ray sources like V407 Cygni in Symbiotic Stars

Authors: Guoliang Lu, Chunhua Zhu, Zhaojun Wang, Wensheng Huo, Yuangui Yang

Abstract: Using a simple accelerating model and an assumption that $γ$-rays originate from $p-p$ collisions for a $π^0$ model, we investigate $γ$-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 $γ$-rays. The maximum energies of the accelerated protons main… ▽ More Using a simple accelerating model and an assumption that $γ$-rays originate from $p-p$ collisions for a $π^0$ model, we investigate $γ$-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 $γ$-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 $γ$-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 $γ$-ray emission like that in V407 Cygni. △ Less

Submitted 31 January, 2011; originally announced January 2011.

Comments: 4pages, 3 figures, Accepted for publication in MNRAS Letter

Soft and Supersoft X-ray Sources in Symbiotic Stars

Authors: Chunhua Zhu, Guoliang Lu, Zhaojun Wang, Jun Zhang

Abstract: 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… ▽ More 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. △ Less

Submitted 23 July, 2009; originally announced July 2009.

Comments: 14 pages, 4 figures, accepted for publication in New Astronomy

Chemical Abundances in Symbiotic Stars

Authors: Guoliang Lu, Chunhua Zhu, Zhanwen Han, Zhaojun Wang

Abstract: 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 $λ$ and the… ▽ More 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 $λ$ 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 $λ$, 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. △ Less

Submitted 6 May, 2008; originally announced May 2008.

Comments: 28 pages, 6 figure, accepted for publication in ApJ

Population synthesis for symbiotic stars with white dwarf accretors

Authors: Guoliang Lu, Lev Yungelson, Zhanwen Han

Abstract: 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 a… ▽ More 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. △ Less

Submitted 22 August, 2006; originally announced August 2006.

Comments: 18 pages, 14 figures, accepted for publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.372:1389-1406,2006