Showing 1–50 of 113 results for author: Zhu, Q

Ultralight scalar dark matter versus non-adiabatic perfect fluid dark matter in pulsar timing

Authors: Qing-Hua Zhu

Abstract: Recent evidence of direct detection of stochastic gravitational waves reported by pulsar timing array collaborations might open a new window for studying cosmology and astrophysical phenomena. In addition to signals from gravitational waves, there is motivation to explore residual signals from oscillating dark matter, which might partially comprise the galactic halo. We investigate fluctuations in… ▽ More Recent evidence of direct detection of stochastic gravitational waves reported by pulsar timing array collaborations might open a new window for studying cosmology and astrophysical phenomena. In addition to signals from gravitational waves, there is motivation to explore residual signals from oscillating dark matter, which might partially comprise the galactic halo. We investigate fluctuations in pulsar timing originating from the coherent oscillation of scalar dark matter up to the subleading order of $\mathcal{O}(k/m)$, as well as from acoustic oscillations of non-adiabatic perfect fluid dark matter. Both types of dark matter can generate oscillating Newtonian potential perturbations and curvature perturbations, thereby affecting pulsar timing. Our results show distinctive signatures in pulsar timing residuals and angular correlations for these dark matters. Specifically, pulsar timing residuals from non-adiabatic perfect fluid dark matter exhibit different directional dependence and are shown to be more sensitive to the distance to a pulsar. We also study the angular correlation patterns from these dark matters in the NANOGrav 15-year data set. The best fit might suggest that the composition of non-adiabatic perfect fluid dark matter in our galaxy is much greater than that of ultralight scalar dark matter. △ Less

Submitted 3 April, 2024; originally announced April 2024.

Comments: 20 pages, 10 figures

AT2023lli: A Tidal Disruption Event with Prominent Optical Early Bump and Delayed Episodic X-ray Emission

Authors: Shifeng Huang, Ning Jiang, Jiazheng Zhu, Yibo Wang, Tinggui Wang, Shan-Qin Wang, Wen-Pei Gan, En-Wei Liang, Yu-Jing Qin, Zheyu Lin, Lin-Na Xu, Min-Xuan Cai, Ji-An Jiang, Xu Kong, Jiaxun Li, Long Li, Jian-Guo Wang, Ze-Lin Xu, Yongquan Xue, Ye-Fei Yuan, Jingquan Cheng, Lulu Fan, Jie Gao, Lei Hu, Weida Hu , et al. (20 additional authors not shown)

Abstract: High-cadence, multiwavelength observations have continuously revealed the diversity of tidal disruption events (TDEs), thus greatly advancing our knowledge and understanding of TDEs. In this work, we conducted an intensive optical-UV and X-ray follow-up campaign of TDE AT2023lli, and found a remarkable month-long bump in its UV/optical light curve nearly two months prior to maximum brightness. The… ▽ More High-cadence, multiwavelength observations have continuously revealed the diversity of tidal disruption events (TDEs), thus greatly advancing our knowledge and understanding of TDEs. In this work, we conducted an intensive optical-UV and X-ray follow-up campaign of TDE AT2023lli, and found a remarkable month-long bump in its UV/optical light curve nearly two months prior to maximum brightness. The bump represents the longest separation time from the main peak among known TDEs to date. The main UV/optical outburst declines as $t^{-4.10}$, making it one of the fastest decaying optically selected TDEs. Furthermore, we detected sporadic X-ray emission 30 days after the UV/optical peak, accompanied by a reduction in the period of inactivity. It is proposed that the UV/optical bump could be caused by the self-intersection of the stream debris, whereas the primary peak is generated by the reprocessed emission of the accretion process. In addition, our results suggest that episodic X-ray radiation during the initial phase of decline may be due to the patched obscurer surrounding the accretion disk, a phenomenon associated with the inhomogeneous reprocessing process. The double TDE scenario, in which two stars are disrupted in sequence, is also a possible explanation for producing the observed early bump and main peak. We anticipate that the multicolor light curves of TDEs, especially in the very early stages, and the underlying physics can be better understood in the near future with the assistance of dedicated surveys such as the deep high-cadence survey of the 2.5-meter Wide Field Survey Telescope (WFST). △ Less

Submitted 26 March, 2024; v1 submitted 3 March, 2024; originally announced March 2024.

Comments: 14 pages, 8 figures,accepted for publication by ApJL

The Jiao Tong University Spectroscopic Telescope Project

Authors: JUST Team, Chengze Liu, Ying Zu, Fabo Feng, Zhaoyu Li, Yu Yu, Hua Bai, Xiangqun Cui, Bozhong Gu, Yizhou Gu, Jiaxin Han, Yonghui Hou, Zhongwen Hu, Hangxin Ji, Yipeng Jing, Wei Li, Zhaoxiang Qi, Xianyu Tan, Cairang Tian, Dehua Yang, Xiangyan Yuan, Chao Zhai, Congcong Zhang, Jun Zhang, Haotong Zhang , et al. (6 additional authors not shown)

Abstract: The Jiao Tong University Spectroscopic Telescope (JUST) is a 4.4-meter f/6.0 segmentedmirror telescope dedicated to spectroscopic observations. The JUST primary mirror is composed of 18 hexagonal segments, each with a diameter of 1.1 m. JUST provides two Nasmyth platforms for placing science instruments. One Nasmyth focus fits a field of view of 10 arcmin and the other has an extended field of vie… ▽ More The Jiao Tong University Spectroscopic Telescope (JUST) is a 4.4-meter f/6.0 segmentedmirror telescope dedicated to spectroscopic observations. The JUST primary mirror is composed of 18 hexagonal segments, each with a diameter of 1.1 m. JUST provides two Nasmyth platforms for placing science instruments. One Nasmyth focus fits a field of view of 10 arcmin and the other has an extended field of view of 1.2 deg with correction optics. A tertiary mirror is used to switch between the two Nasmyth foci. JUST will be installed at a site at Lenghu in Qinghai Province, China, and will conduct spectroscopic observations with three types of instruments to explore the dark universe, trace the dynamic universe, and search for exoplanets: (1) a multi-fiber (2000 fibers) medium-resolution spectrometer (R=4000-5000) to spectroscopically map galaxies and large-scale structure; (2) an integral field unit (IFU) array of 500 optical fibers and/or a long-slit spectrograph dedicated to fast follow-ups of transient sources for multimessenger astronomy; (3) a high-resolution spectrometer (R~100000) designed to identify Jupiter analogs and Earth-like planets, with the capability to characterize the atmospheres of hot exoplanets. △ Less

Submitted 29 February, 2024; v1 submitted 22 February, 2024; originally announced February 2024.

Comments: 28 pages, 6 figures

Intrinsic bispectrum of the scalar-induced gravitational waves

Authors: Qing-Hua Zhu

Abstract: Recent pulsar timing array collaborations have reported evidence of the stochastic gravitational wave background. The gravitational waves induced by primordial curvature perturbations, referred to as scalar-induced gravitational waves (SIGWs), could potentially be the physical origins of this gravitational wave background. In addition to the statistical properties of SIGWs derived from the primord… ▽ More Recent pulsar timing array collaborations have reported evidence of the stochastic gravitational wave background. The gravitational waves induced by primordial curvature perturbations, referred to as scalar-induced gravitational waves (SIGWs), could potentially be the physical origins of this gravitational wave background. In addition to the statistical properties of SIGWs derived from the primordial fluctuations, SIGWs also have intrinsic non-Gaussianity originating from the non-linear interactions of Einstein's gravity, which, however, is rarely explored. In this paper, we study the intrinsic non-Gaussianity of SIGWs with the bispectrum and the skewness under the assumption of the Gaussian primordial curvature perturbations, phenomenologically modeled as a lognormal spectrum. The bispectrum is shown to be vanishing in the collinear limit, which is independent of the initial conditions and the dynamics of SIGWs. For the SIGWs generated during the radiation-dominated era, the bispectrum is of flatten-type non-Gaussianity, with only four polarization components left to be non-vanishing. Additionally, in order to obtain the correct bispectrum, we also propose a time-oscillation average scheme and a regularization scheme. Utilizing the skewness for quantifying the degree of non-Gaussianity, it is found that the curvature power spectrum with a narrow width can result in an enhancement of the third-order non-Gaussianity. The conclusion holds for both the SIGWs generated in the radiation-dominated era and the matter-dominated era. △ Less

Submitted 4 February, 2024; originally announced February 2024.

Comments: 31 pages, 11 figures

Physical Parameters of 11,100 Short-Period ASAS-SN Eclipsing Contact Binaries

Authors: Xu-Zhi Li, Qing-Feng Zhu, Xu Ding, Xiao-Hui Xu, Hang Zheng, Jin-Sheng Qiu, Ming-Chao Liu

Abstract: Starting from more than 11,200 short-period (less than 0.5 days) EW-type eclipsing binary candidates with the All-Sky Automated Survey for Supernovae (ASAS-SN) V-band light curves, we use MCMC and neural networks (NNs) to obtain the mass ratio ($q$), orbital inclination ($incl$), fill-out factor ($f$) and temperature ratio ($T_s/T_p$). After cross-matching with the Gaia DR3 database, the final sam… ▽ More Starting from more than 11,200 short-period (less than 0.5 days) EW-type eclipsing binary candidates with the All-Sky Automated Survey for Supernovae (ASAS-SN) V-band light curves, we use MCMC and neural networks (NNs) to obtain the mass ratio ($q$), orbital inclination ($incl$), fill-out factor ($f$) and temperature ratio ($T_s/T_p$). After cross-matching with the Gaia DR3 database, the final sample contains parameters of 2,399 A-type and 8,712 W-type contact binaries (CBs). We present the distributions of parameters of these 11,111 short-period CBs. The mass ratio ($q$) and fill-out factor ($f$) are found to obey log-normal distributions, and the remaining parameters obey normal distributions. There is a significant period-temperature correlation of these CBs. Additionally, the temperature ratio (${T_s}$/${T_p}$) tends to increase as the orbital period decreases for W-type CBs. There is no significant correlation between them for A-type CBs. The mass ratio and fill-out factor ($q-f$) diagram suggest there is no significant correlation between these two parameters. A clear correlation exists between the mass ratio and radius ratio. The radius ratio increases with the mass ratio. Moreover, the deep fill-out CBs tend to fall on the upper boundary of the $q$$-$${R_s}$/${R_p}$ distribution, while the shallow fill-out CBs fall on the lower boundary. △ Less

Submitted 29 January, 2024; originally announced January 2024.

Comments: 15 pages, 8 figures, 1 table.Accepted for publication in ApJS. The full dataset will be available with the ApJS article

New constraints on primordial non-Gaussianity from missing two-loop contributions of scalar induced gravitational waves

Authors: Zhe Chang, Yu-Ting Kuang, Di Wu, Jing-Zhi Zhou, Qing-Hua Zhu

Abstract: We analyze the energy density spectrum of \acp{SIGW} using the NANOGrav 15-year data set, thereby constraining the primordial non-Gaussian parameter $f_{\mathrm{NL}}$. For the first time, we calculate the seventeen missing two-loop diagrams proportional to $f_{\mathrm{NL}}A_ζ^3$ that correspond to the two-point correlation function $\langle h^{λ,(3)}_{\mathbf{k}} h^{λ',(2)}_{\mathbf{k}'} \rangle$… ▽ More We analyze the energy density spectrum of \acp{SIGW} using the NANOGrav 15-year data set, thereby constraining the primordial non-Gaussian parameter $f_{\mathrm{NL}}$. For the first time, we calculate the seventeen missing two-loop diagrams proportional to $f_{\mathrm{NL}}A_ζ^3$ that correspond to the two-point correlation function $\langle h^{λ,(3)}_{\mathbf{k}} h^{λ',(2)}_{\mathbf{k}'} \rangle$ for local-type primordial non-Gaussianity. The total energy density spectrum of \acp{SIGW} can be significantly suppressed by these two-loop diagrams. If \acp{SIGW} dominate the \acp{SGWB} observed in \ac{PTA} experiments, the parameter interval $f_{\mathrm{NL}}\in [-5,-1]$ is notably excluded based on NANOGrav 15-year data set. After taking into account abundance of \acp{PBH} and the convergence of the cosmological perturbation expansion, we find that the only possible parameter range for $f_{\mathrm{NL}}$ might be $-1\le f_{\mathrm{NL}}< 0$. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: 8 pages, 4 figures

A Gigantic Mid-Infrared Outburst in an Embedded Class-I Young Stellar Object J064722.95+031644.6

Authors: Tinggui Wang, Jiaxun Li, Gregory M. Mace, Tuo Ji, Ning Jiang, Qingfeng Zhu, Min Fang

Abstract: We report the serendipitous discovery of a giant mid-infrared (MIR) outburst from a previously unknown source near a star-forming region in the constellation Monoceros. The source gradually brightened by a factor of 5 from 2014 to 2016 before an abrupt rise by a factor of more than 100 in 2017. A total amplitude increase of >500 at 4.5 microns has since faded by a factor of about 10. Prior to the… ▽ More We report the serendipitous discovery of a giant mid-infrared (MIR) outburst from a previously unknown source near a star-forming region in the constellation Monoceros. The source gradually brightened by a factor of 5 from 2014 to 2016 before an abrupt rise by a factor of more than 100 in 2017. A total amplitude increase of >500 at 4.5 microns has since faded by a factor of about 10. Prior to the outburst, it was only detected at wavelengths longer than 1.8 microns in UKIDSS, Spitzer, and Herschel with a spectral energy distribution of a Class I Young Stellar Object (YSO). It has not been detected in recent optical surveys, suggesting that it is deeply embedded. With a minimum distance of 3.5 kpc, the source has a bolometric luminosity of at least 9 $L_\odot$ in the quiescent state and 400 $L_\odot$ at the peak of the eruption. The maximum accretion rate is estimated to be at least a few $10^{-5}$ $M_\odot$ year$^{-1}$. It shares several common properties with another eruptive event, WISE~J142238.82-611553.7: exceptionally large amplitude, featureless near-infrared spectrum with the exception of H_2 lines, intermediate eruption duration, an embedded Class I YSO, and a low radiative temperature (<600-700 K) in outburst. We interpret that the radiation from the inner accretion disk and young star is obscured and reprocessed by either an inflated outer disk or thick dusty outflow on scales > 6.5 AU during the outburst. △ Less

Submitted 19 September, 2023; originally announced September 2023.

Comments: 15 pages, 7 figures, accepted to ApJ

Unraveling the early universe's equation of state and primordial black hole production with PTA, BBN, and CMB observations

Authors: Qing-Hua Zhu, Zhi-Chao Zhao, Sai Wang, Xin Zhang

Abstract: Pulsar timing array (PTA) data releases showed strong evidence for a stochastic gravitational-wave background in the nanohertz band. When the signal is interpreted by a scenario of scalar-induced gravitational waves (SIGWs), we encounter overproduction of primordial black holes (PBHs). We wonder if varying the equation of state (EoS) of the early Universe can resolve this issue and thereby lead to… ▽ More Pulsar timing array (PTA) data releases showed strong evidence for a stochastic gravitational-wave background in the nanohertz band. When the signal is interpreted by a scenario of scalar-induced gravitational waves (SIGWs), we encounter overproduction of primordial black holes (PBHs). We wonder if varying the equation of state (EoS) of the early Universe can resolve this issue and thereby lead to a consistent interpretation of the PTA data. Analyzing a data combination of PTA, big-bang nucleosynthesis, and cosmic microwave background, we find that an epoch with EoS $w\sim\mathcal{O}(10^{-2})$ between the end of inflation and the onset of radiation domination can significantly suppress the production of PBHs, leading to alleviation of the PBH-overproduction issue. With the inferred interval $w=0.44_{-0.40}^{+0.52}$ at 95\% confidence level, our scenario can interpret the PTA data just as well as the conventional scenario of SIGWs produced during the radiation domination. △ Less

Submitted 5 June, 2024; v1 submitted 25 July, 2023; originally announced July 2023.

Comments: 15 pages, 3 figures, 1 column, major revision

Constraints On Scalar-Induced Gravitational Waves Up To Third Order From Joint Analysis of BBN, CMB, And PTA Data

Authors: Sai Wang, Zhi-Chao Zhao, Qing-Hua Zhu

Abstract: Recently, strong evidence for a gravitational wave background has been reported by collaborations of pulsar timing arrays (PTA). In the framework of scalar-induced gravitational waves (SIGWs), we concurrently investigate the second and third order gravitational waves by jointly analyzing PTA data, alongside big-bang nucleosynthesis (BBN), and cosmic microwave background (CMB) datasets. We determin… ▽ More Recently, strong evidence for a gravitational wave background has been reported by collaborations of pulsar timing arrays (PTA). In the framework of scalar-induced gravitational waves (SIGWs), we concurrently investigate the second and third order gravitational waves by jointly analyzing PTA data, alongside big-bang nucleosynthesis (BBN), and cosmic microwave background (CMB) datasets. We determine the primordial curvature spectral amplitude as $0.021<A_ζ<0.085$ and the spectral peak frequency as $10^{-7.3}\ \mathrm{Hz}<f_\ast<10^{-6.3}\ \mathrm{Hz}$ at a 95\% confidence interval, pointing towards a mass range for primordial black holes of $10^{-4.5}M_\odot<m_{\mathrm{PBH}}<10^{-2.5}M_\odot$. Our findings suggest that third order gravitational waves contribute more significantly to the integrated energy density than the second order ones when $A_ζ\gtrsim0.06$. Furthermore, we expect future PTA projects to validate these findings and provide robust means to investigate the genesis and evolution of the universe, especially inflation. △ Less

Submitted 24 January, 2024; v1 submitted 6 July, 2023; originally announced July 2023.

Comments: 10 pages, 2 columns, 3 figures. Version accepted by Physical Review Research

Journal ref: Phys. Rev. Research 6, 013207 (2024)

Primordial black holes from second order density perturbations as probes of the small-scale primordial power spectrum

Authors: Yu-Ting Kuang, Jing-Zhi Zhou, Zhe Chang, Xukun Zhang, Qing-Hua Zhu

Abstract: We investigate the second order energy density perturbation $δ^{(2)}$ induced by small-scale Gaussian and local-type non-Gaussian primordial curvature perturbations. The relative abundance of primordial black hole is calculated in terms of the probability density function of total energy density perturbation $δ_r=δ^{(1)}+\frac{1}{2}δ^{(2)}$. The effects of second order density perturbation greatly… ▽ More We investigate the second order energy density perturbation $δ^{(2)}$ induced by small-scale Gaussian and local-type non-Gaussian primordial curvature perturbations. The relative abundance of primordial black hole is calculated in terms of the probability density function of total energy density perturbation $δ_r=δ^{(1)}+\frac{1}{2}δ^{(2)}$. The effects of second order density perturbation greatly reduce the upper bounds of small-scale power spectra of primordial curvature perturbations by one to two orders of magnitude. For log-normal primordial power spectrum, its amplitude $A_ζ$ is constrained to be about $A_ζ\sim 3\times10^{-3}$. And for local-type non-Gaussianity with $f_{\mathrm{NL}}=10$, the upper bound of $A_ζ$ is about $2.5\times10^{-4}$. △ Less

Submitted 7 July, 2023; v1 submitted 5 July, 2023; originally announced July 2023.

Comments: 5 pages, 4 figures

Implications of Pulsar Timing Array Data for Scalar-Induced Gravitational Waves and Primordial Black Holes: Primordial Non-Gaussianity $f_{\mathrm{NL}}$ Considered

Authors: Sai Wang, Zhi-Chao Zhao, Jun-Peng Li, Qing-Hua Zhu

Abstract: Multiple pulsar-timing-array collaborations have reported strong evidence for the existence of a gravitational-wave background. We study physical implications of this signal for cosmology, assuming that it is attributed to scalar-induced gravitational waves. By incorporating primordial non-Gaussianity $f_{\mathrm{NL}}$, we specifically examine the nature of primordial curvature perturbations and p… ▽ More Multiple pulsar-timing-array collaborations have reported strong evidence for the existence of a gravitational-wave background. We study physical implications of this signal for cosmology, assuming that it is attributed to scalar-induced gravitational waves. By incorporating primordial non-Gaussianity $f_{\mathrm{NL}}$, we specifically examine the nature of primordial curvature perturbations and primordial black holes. We find that the signal allows for a primordial non-Gaussianity $f_{\mathrm{NL}}$ in the range of $-4.1\lesssim f_{\mathrm{NL}} \lesssim 4.1$ (68\% confidence intervals) and a mass range for primordial black holes $m_{\mathrm{pbh}}$ spanning from $\sim10^{-5}M_{\odot}$ to $\sim10^{-2}M_{\odot}$. Furthermore, we find that the signal favors a negative non-Gaussianity, which can suppress the abundance of primordial black holes. We also demonstrate that the anisotropies of scalar-induced gravitational waves serve as a powerful tool to probe the non-Gaussianity $f_{\mathrm{NL}}$. We conduct a comprehensive analysis of the angular power spectrum within the nano-Hertz band. Looking ahead, we anticipate that future projects, such as the Square Kilometre Array, will have the potential to measure these anisotropies and provide further insights into the primordial universe. △ Less

Submitted 27 September, 2023; v1 submitted 2 July, 2023; originally announced July 2023.

Comments: Major revisions, 2 columns, 8 pages, 4 figures, 1 appendix

Journal ref: Phys. Rev. Research 6, L012060 (2024)

Origins of the shocks in high-mass starless clump candidates

Authors: Feng-Yao Zhu, Junzhi Wang, Yaoting Yan, Qing-Feng Zhu, Juan Li

Abstract: Shocks are abundant in star-forming regions, and are often related with star formation. In our previous observations toward 100 starless clump candidates (SCCs) in the Galaxy, a sample of 34 SCCs associated with shocks is identified. In this work, we perform mapping observations of the SiO 2-1, 3-2, HC$_3$N 10-9, HCO$^+$ 1-0, H$^{13}$CO$^+$ 1-0, and H41$α$ lines toward 9 out of the detected source… ▽ More Shocks are abundant in star-forming regions, and are often related with star formation. In our previous observations toward 100 starless clump candidates (SCCs) in the Galaxy, a sample of 34 SCCs associated with shocks is identified. In this work, we perform mapping observations of the SiO 2-1, 3-2, HC$_3$N 10-9, HCO$^+$ 1-0, H$^{13}$CO$^+$ 1-0, and H41$α$ lines toward 9 out of the detected sources by using IRAM 30-m radio telescope to study the origins of the shocks in the SCCs. We find shocks in three sources (BGPS 3110, 3114, and 3118) are produced by collisions between the expanding ionized gas and ambient molecular gas, instead of by the star formation activity inside SCCs. On the other hand, shocks in the other six sources are related to star formation activity of SCCs. The signatures of protostellar outflows are clearly shown in the molecular lines toward BGPS 4029, 4472, 5064. Comparing our results with the previous ALMA observations performed in the same region, the shocks in BGPS 3686 and 5114 are also likely to be due to protostellar activity. The origin of shock in BGPS 5243 is still unclear although some features in the SiO spectra imply the presence of protostellar activity. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Comments: 21 pages, 21 figues

Journal ref: MNRAS 523, 2770 (2023)

Spatial distributions and kinematics of shocked and ionized gas in M17

Authors: Feng-Yao Zhu, Junzhi Wang, Yaoting Yan, Qing-Feng Zhu, Juan Li

Abstract: Massive stars are formed in molecular clouds, and produce H II regions when they evolve onto the main sequence. The expansion of H II region can both suppress and promote star formation in the vicinity. M17 H II region is a giant cometary H II region near many massive clumps containing starless and protostellar sources. It is an appropriate target to study the effect of feedback from previously fo… ▽ More Massive stars are formed in molecular clouds, and produce H II regions when they evolve onto the main sequence. The expansion of H II region can both suppress and promote star formation in the vicinity. M17 H II region is a giant cometary H II region near many massive clumps containing starless and protostellar sources. It is an appropriate target to study the effect of feedback from previously formed massive stars on the nearby star-forming environments. Observations of SiO 2-1, HCO$^+$ 1-0, H$^{13}$CO$^+$ 1-0, HC$_3$N 10-9, and H41$α$ lines are performed toward M17 H II region with ambient candidates of massive clumps. In the observations, the widespread shocked gas surrounding M17 H II region is detected: it probably originates from the collision between the expanding ionized gas and the ambient neutral medium. Some massive clumps are found in the overlap region of the shock and dense-gas tracing lines while the central velocities of shocked and high-density gases are similar. This suggests that part of massive clumps are located in the shell of H II region, and may be formed from the accumulated neutral materials in the shell. In addition, by comparing the observations toward M17 H II region with the simulation of cometary H II region, we infer the presence of one or more massive stars travelling at supersonic velocity with respect to the natal molecular cloud in the H II region. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Comments: 13 pages, 11 figures

Journal ref: MNRAS 522, 503 (2023)

Sciences with the 2.5-meter Wide Field Survey Telescope (WFST)

Authors: WFST Collaboration, Tinggui Wang, Guilin Liu, Zhenyi Cai, Jinjun Geng, Min Fang, Haoning He, Ji-an Jiang, Ning Jiang, Xu Kong, Bin Li, Ye Li, Wentao Luo, Zhizheng Pan, Xuefeng Wu, Ji Yang, Jiming Yu, Xianzhong Zheng, Qingfeng Zhu, Yi-Fu Cai, Yuanyuan Chen, Zhiwei Chen, Zigao Dai, Lulu Fan, Yizhong Fan , et al. (38 additional authors not shown)

Abstract: The Wide Field Survey Telescope (WFST) is a dedicated photometric surveying facility being built jointly by the University of Science and Technology of China and the Purple Mountain Observatory. It is equipped with a 2.5-meter diameter primary mirror, an active optics system, and a mosaic CCD camera with 0.73 gigapixels on the primary focal plane for high-quality image capture over an FOV of 6.5-s… ▽ More The Wide Field Survey Telescope (WFST) is a dedicated photometric surveying facility being built jointly by the University of Science and Technology of China and the Purple Mountain Observatory. It is equipped with a 2.5-meter diameter primary mirror, an active optics system, and a mosaic CCD camera with 0.73 gigapixels on the primary focal plane for high-quality image capture over an FOV of 6.5-square-degree. It is anticipated that WFST will be set up at the Lenghu site in the summer of 2023 and begin to observe the northern sky in four optical bands (u, g, r, and i) with a range of cadences, from hourly/daily in the Deep High-Cadence Survey (DHS) program to semiweekly in the Wide-Field Survey (WFS) program, three months later. During a photometric night, a nominal 30 s exposure in the WFS program will reach a depth of 22.27, 23.32, 22.84, and 22.31 (AB magnitudes) in these four bands, respectively, allowing for the detection of a tremendous amount of transients in the low-z universe and a systematic investigation of the variability of Galactic and extragalactic objects. In the DHS program, intranight 90 s exposures as deep as 23 (u) and 24 mag (g), in combination with target of opportunity follow-ups, will provide a unique opportunity to explore energetic transients in demand for high sensitivities, including the electromagnetic counterparts of gravitational wave events, supernovae within a few hours of their explosions, tidal disruption events and fast, luminous optical transients even beyond a redshift of unity. In addition, the final 6-year co-added images, anticipated to reach g=25.8 mag in WFS or 1.5 mags deeper in DHS, will be of fundamental importance to general Galactic and extragalactic science. The highly uniform legacy surveys of WFST will serve as an indispensable complement to those of LSST that monitor the southern sky. △ Less

Submitted 14 September, 2023; v1 submitted 13 June, 2023; originally announced June 2023.

Comments: 48 pages

Journal ref: SCPMA-Vol. 66 No. 10: 109512 (2023)

Limiting Magnitudes of the Wide Field Survey Telescope (WFST)

Authors: Lei Lei, Qing-Feng Zhu, Xu Kong, Ting-Gui Wang, Xian-Zhong Zheng, Dong-Dong Shi, Lu-Lu Fan, Wei Liu

Abstract: Expected to be of the highest survey power telescope in the northern hemisphere, the Wide Field Survey Telescope (WFST) will begin its routine observations of the northern sky since 2023. WFST will produce a lot of scientific data to support the researches of time-domain astronomy, asteroids and the solar system, galaxy formation and cosmology and so on. We estimated that the 5 $σ$ limiting magnit… ▽ More Expected to be of the highest survey power telescope in the northern hemisphere, the Wide Field Survey Telescope (WFST) will begin its routine observations of the northern sky since 2023. WFST will produce a lot of scientific data to support the researches of time-domain astronomy, asteroids and the solar system, galaxy formation and cosmology and so on. We estimated that the 5 $σ$ limiting magnitudes of WFST with 30 second exposure are $u=22.31$ mag, $g=23.42$ mag, $r=22.95$ mag, $i=22.43$ mag, $z=21.50$ mag, $w=23.61$ mag. The above values are calculated for the conditions of $airmass=1.2$, seeing = 0.75 arcsec, precipitable water vapour (PWV) = 2.5 mm and Moon-object separation = $45^{\circ}$ at the darkest New Moon night of the Lenghu site (V=22.30 mag, Moon phase $θ=0^{\circ}$). The limiting magnitudes in different Moon phase conditions are also calculated. The calculations are based on the empirical transmittance data of WFST optics, the vendor provided CCD quantum efficiency, the atmospherical model transmittance and spectrum of the site. In the absence of measurement data such as sky transmittance and spectrum, we use model data. △ Less

Submitted 24 September, 2023; v1 submitted 8 January, 2023; originally announced January 2023.

Comments: Figure 5 (a,b,c) has been updated in this latest version. 12 pages, 5 figures, published in RAA (Research in Astronomy and Astrophysics)

Journal ref: Research in Astronomy and Astrophysics, 23:035013 (8pp), 2023 March

Photon ring autocorrelations from gravitational fluctuations around a black hole

Authors: Qing-Hua Zhu

Abstract: The images of supermassive black holes in M87 and our galaxy captured by the Event Horizon Telescope (EHT) might open up a new way for exploring black hole physics at the horizon scale. Theoretically, this could provide insights into studying the emission sources around a black hole or the geometries of the black hole itself. This paper investigates the two-point correlations of intensity fluctuat… ▽ More The images of supermassive black holes in M87 and our galaxy captured by the Event Horizon Telescope (EHT) might open up a new way for exploring black hole physics at the horizon scale. Theoretically, this could provide insights into studying the emission sources around a black hole or the geometries of the black hole itself. This paper investigates the two-point correlations of intensity fluctuations on the photon ring, resulting from the existence of shock waves or gravitational fluctuations around a black hole. Following approaches used in the field of gravitational wave detectors, we introduce response functions of very long baseline interferometry (VLBI) for detecting gravitational fluctuations and study the shape of the overlap reduction functions. It is found that the shape of the correlations here differs from that resulting from stochastic emission sources. By providing an order-of-magnitude estimate for the signal-to-noise ratio (SNR), we obtain sensitivity curves for the detection of gravitational fluctuations. This might reveal the potential for detecting gravitational fluctuations with future VLBI observations in the LIGO frequency band. △ Less

Submitted 14 March, 2024; v1 submitted 2 January, 2023; originally announced January 2023.

Comments: v2: accepted version, 24 pages, 11 figures, 11pt

Journal ref: Phys. Rev. D 109, 064031 (2024)

Non-linear corrections of overlap reduction functions for pulsar timing arrays

Authors: Qing-Hua Zhu

Abstract: The signals from international pulsar timing arrays have presented a hint of gravitational stochastic background in nHz band frequency. Further confirmation will be based on whether the signals follow the angular correlation curves formulated by the overlap reduction functions, known as Hellings-Downs curves. This paper investigates the non-linear corrections of overlap reduction functions in the… ▽ More The signals from international pulsar timing arrays have presented a hint of gravitational stochastic background in nHz band frequency. Further confirmation will be based on whether the signals follow the angular correlation curves formulated by the overlap reduction functions, known as Hellings-Downs curves. This paper investigates the non-linear corrections of overlap reduction functions in the present of non-Gaussianity, in which the self-interaction of gravity is first taken into considerations. Based on perturbed Einstein field equations for the second order metric perturbations, and perturbed geodesic equations to the second order, we obtain non-linear corrections for the timing residuals of pulsar timing, and theoretically study corresponding overlap reduction functions for pulsar timing arrays. There is order-one correction for the overlap reduction functions from the three-point correlations of gravitational waves, and thus the shapes of the overlap reduction functions with non-linear corrections can be distinguished from the Hellings-Downs curves. △ Less

Submitted 4 February, 2024; v1 submitted 31 December, 2022; originally announced January 2023.

Comments: v1: 20 pages, 8 figures. v2: accepted version, 23 pages, 11 figures, added the discussion about scalar-induced GWs

A New Period Determination Method for Periodic Variable Stars

Authors: Xiao-Hui Xu, Qing-Feng Zhu, Xu-Zhi Li, Bin Li, Hang Zheng, Jin-Sheng Qiu, Hai-Bin Zhao

Abstract: Variable stars play a key role in understanding the Milky Way and the universe. The era of astronomical big data presents new challenges for quick identification of interesting and important variable stars. Accurately estimating the periods is the most important step to distinguish different types of variable stars. Here, we propose a new method of determining the variability periods. By combining… ▽ More Variable stars play a key role in understanding the Milky Way and the universe. The era of astronomical big data presents new challenges for quick identification of interesting and important variable stars. Accurately estimating the periods is the most important step to distinguish different types of variable stars. Here, we propose a new method of determining the variability periods. By combining the statistical parameters of the light curves, the colors of the variables, the window function and the GLS algorithm, the aperiodic variables are excluded and the periodic variables are divided into eclipsing binaries and NEB variables (other types of periodic variable stars other than eclipsing binaries), the periods of the two main types of variables are derived. We construct a random forest classifier based on 241,154 periodic variables from the ASAS-SN and OGLE datasets of variables. The random forest classifier is trained on 17 features, among which 11 are extracted from the light curves and 6 are from the Gaia Early DR3, ALLWISE and 2MASS catalogs. The variables are classified into 7 superclasses and 17 subclasses. In comparison with the ASAS-SN and OGLE catalogs, the classification accuracy is generally above approximately 82% and the period accuracy is 70%-99%. To further test the reliability of the new method and classifier, we compare our results with the results of Chen et al. (2020) for ZTF DR2. The classification accuracy is generally above 70%. The period accuracy of the EW and SR variables is 50% and 53%, respectively. And the period accuracy of other types of variables is 65%-98%. △ Less

Submitted 20 December, 2022; originally announced December 2022.

Comments: 23 pages, 10 figures

Journal ref: PASP2022

Overlap reduction function for gravitational wave detectors in an expanding Universe

Authors: Qing-Hua Zhu

Abstract: Since it was confirmed two decades ago that the expansion of the Universe is accelerating, it would be of theoretical interests to figure out what is the influence from cosmological constant on detection of stochastic gravitational wave background. This paper studies the overlap reduction functions in de-Sitter space-time for a pair of one-way tracking gravitational wave detectors. It is shown to… ▽ More Since it was confirmed two decades ago that the expansion of the Universe is accelerating, it would be of theoretical interests to figure out what is the influence from cosmological constant on detection of stochastic gravitational wave background. This paper studies the overlap reduction functions in de-Sitter space-time for a pair of one-way tracking gravitational wave detectors. It is shown to be non-trivial in an expanding Universe, because the propagation of light along line of sight also has effect on the response of GW detectors. It is found that the expansion of the Universe can enhance the value of magnitude of the overlap reduction functions, when the detector pairs are close to each other. For nanohertz gravitational waves, this effect can dominate the values of overlap reduction functions when the galactic pulsar pairs are separated by milliarcsecond. △ Less

Submitted 5 July, 2022; originally announced July 2022.

Comments: 10 pages, 3 figures

Journal ref: Eur. Phys. J. C 84, 2 (2024)

Inverted level populations of hydrogen atoms in ionized gas

Authors: Feng-Yao Zhu, Junzhi Wang, Qing-Feng Zhu, Jiang-Shui Zhang

Abstract: Context. Level population inversion of hydrogen atoms in ionized gas may lead to stimulated emission of hydrogen recombination lines, and the level populations can in turn be affected by powerful stimulated emissions. Aims. In this work the interaction of the radiation fields and the level population inversion of hydrogen atoms is studied. The effect of the stimulated emissions on the line profile… ▽ More Context. Level population inversion of hydrogen atoms in ionized gas may lead to stimulated emission of hydrogen recombination lines, and the level populations can in turn be affected by powerful stimulated emissions. Aims. In this work the interaction of the radiation fields and the level population inversion of hydrogen atoms is studied. The effect of the stimulated emissions on the line profiles is also investigated. Methods. Our previous nl-model for calculating level populations of hydrogen atoms and hydrogen recombination lines is improved. The effects of line and continuum radiation fields on the level populations are considered in the improved model. By using this method the properties of simulated hydrogen recombination lines and level populations are used in analyses. Results. The simulations show that hydrogen radio recombination lines are often emitted from the energy level with an inverted population. The widths of Hn$α$ lines can be significantly narrowed by strong stimulated emissions to be even less than 10 km s$^{-1}$. The amplification of hydrogen recombination lines is more affected by the line optical depth than by the total optical depth. The influence of stimulated emission on the estimates of electron temperature and density of ionized gas is evaluated. We find that comparing multiple line-to-continuum ratios is a reliable method for estimating the electron temperature, while the effectiveness of the estimation of electron density is determined by the relative significance of stimulated emission. △ Less

Submitted 19 June, 2022; originally announced June 2022.

Comments: Accepted for published in A&A. 25 pages, 13 figures

Journal ref: A&A 665, A94 (2022)

The formation of low surface brightness galaxies in the IllustrisTNG simulation

Authors: Luis Enrique Pérez-Montaño, Vicente Rodriguez-Gomez, Bernardo Cervantes Sodi, Qirong Zhu, Annalisa Pillepich, Mark Vogelsberger, Lars Hernquist

Abstract: We explore the nature of low surface brightness galaxies (LSBGs) in the hydrodynamic cosmological simulation TNG100 of the IllustrisTNG project, selecting a sample of LSBGs ($r$-band effective surface brightness $μ_r > 22.0$ mag arcsec$^{-2}$) at $z=0$ over a wide range of stellar masses ($M_{\ast} = 10^{9}$-$10^{12}$ M$_\odot$). We find LSBGs of all stellar masses, although they are particularly… ▽ More We explore the nature of low surface brightness galaxies (LSBGs) in the hydrodynamic cosmological simulation TNG100 of the IllustrisTNG project, selecting a sample of LSBGs ($r$-band effective surface brightness $μ_r > 22.0$ mag arcsec$^{-2}$) at $z=0$ over a wide range of stellar masses ($M_{\ast} = 10^{9}$-$10^{12}$ M$_\odot$). We find LSBGs of all stellar masses, although they are particularly prevalent at $M_{\ast} < 10^{10}$ M$_\odot$. We show that the specific star formation rates of LSBGs are not significantly different from those of high surface brightness galaxies (HSBGs) but, as a population, LSBGs are systematically less massive and more extended than HSBGs, and tend to display late-type morphologies according to a kinematic criterion. At fixed stellar mass, we find that haloes hosting LSBGs are systematically more massive and have a higher baryonic fraction than those hosting HSBGs. We find that LSBGs have higher stellar specific angular momentum and halo spin parameter values compared to HSBGs, as suggested by previous works. We track the evolution of these quantities back in time, finding that the spin parameters of the haloes hosting LSBGs and HSBGs exhibit a clear bifurcation at $z \sim 2$, which causes a similar separation in the evolutionary tracks of other properties such as galactic angular momentum and effective radius, ultimately resulting in the values observed at $z =$ 0. The higher values of specific stellar angular momentum and halo spin in LSBGs seem to be responsible for their extended nature, preventing material from collapsing into the central regions of the galaxies, also causing LSBGs to host less massive black holes at their centres. △ Less

Submitted 17 June, 2022; originally announced June 2022.

Comments: 14 pages, 10 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

Evolution of Massive Black Hole Binaries in Collisionally Relaxed Nuclear Star Clusters -- Impact of Mass Segregation

Authors: Diptajyoti Mukherjee, Qirong Zhu, Go Ogiya, Carl L. Rodriguez, Hy Trac

Abstract: Massive Black Hole (MBH) binaries are considered to be one of the most important sources of Gravitational Waves (GW) that can be detected by GW detectors like LISA. However, there are a lot of uncertainties in the dynamics of MBH binaries in the stages leading up to the GW-emission phase. It has been recently suggested that Nuclear Star Clusters (NSCs) could provide a viable route to overcome the… ▽ More Massive Black Hole (MBH) binaries are considered to be one of the most important sources of Gravitational Waves (GW) that can be detected by GW detectors like LISA. However, there are a lot of uncertainties in the dynamics of MBH binaries in the stages leading up to the GW-emission phase. It has been recently suggested that Nuclear Star Clusters (NSCs) could provide a viable route to overcome the final parsec problem for MBH binaries at the center of galaxies. NSCs are collisional systems where the dynamics would be altered by the presence of a mass spectrum. In this study, we use a suite of high-resolution $N$-body simulations with over 1 million particles to understand how collisional relaxation under the presence of a mass spectrum of NSC particles affects the dynamics of the MBH binary under the merger of two NSCs. We consider MBH binaries with different mass ratios and additional non-relaxed models. We find that mass-segregation driven by collisional relaxation can lead to accelerated hardening in lower mass ratio binaries but has the opposite effect in higher mass ratio binaries. Crucially, the relaxed models also demonstrate much lower eccentricities at binary formation and negligible growth during hardening stages leading to longer merger timescales. The results are robust and highlight the importance of collisional relaxation on changing the dynamics of the binary. Our models are state-of-the-art, use zero softening, and high enough particle numbers to model NSCs realistically. △ Less

Submitted 16 November, 2022; v1 submitted 24 May, 2022; originally announced May 2022.

Comments: 18 pages, 14 figures. New version with higher resolution simulations. Accepted in MNRAS

A deep investigation of two poorly studied open clusters Haffner 22 and Melotte 71 in Gaia era

Authors: D. Bisht, Quingfeng Zhu, R. K. S. Yadav, Geeta Rangwal, Devesh P. Sariya, Alok Durgapal, Ing-Guey Jiang

Abstract: This paper presents a deep investigation of two open clusters, Haffner 22 and Melotte 71, using astrometric and photometric data from Gaia EDR3. We identified 382 and 597 most probable cluster members with membership probability higher than 50 percent. Mean proper motion in RA and DEC are estimated as -1.63 and 2.889 respectively for Haffner 22 and -2.398 and 4.210 for Melotte 71. A comparison of… ▽ More This paper presents a deep investigation of two open clusters, Haffner 22 and Melotte 71, using astrometric and photometric data from Gaia EDR3. We identified 382 and 597 most probable cluster members with membership probability higher than 50 percent. Mean proper motion in RA and DEC are estimated as -1.63 and 2.889 respectively for Haffner 22 and -2.398 and 4.210 for Melotte 71. A comparison of observed CMDs with theoretical isochrones leads to an age of 2.25 and 1.27 Gyr for these clusters. The distances 2.88 and 2.28 kpc based on parallax are comparable with the values derived by isochrone fitting method. Five and four blue straggler stars are identified as cluster members in Haffner 22 and Melotte 71 respectively. Based on the relative number of high velocity (binary) and single stars, we inferred binary fractions for both clusters. We found binary content is larger in core region. Mass function slope is in good agreement with the Salpeter's value while it is flat for Haffner 22. Evidence for the existence of mass segregation effect is observed in both clusters. Using the Galactic potential models, Galactic orbits are derived, indicating that both clusters follow a circular path around the Galactic center, evolving slowly. △ Less

Submitted 30 March, 2022; originally announced March 2022.

Comments: This article contains 15 figures and 9 tables. This paper has been accepted for the publication in Publications of the Astronomical Society of the Pacific. arXiv admin note: text overlap with arXiv:2103.04596

Magnetic Field Reversal around an Active Fast Radio Burst

Authors: S. Dai, Y. Feng, Y. P. Yang, Y. K. Zhang, D. Li, C. H. Niu, P. Wang, M. Y. Xue, B. Zhang, S. Burke-Spolaor, C. J. Law, R. S. Lynch, L. Connor, R. Anna-Thomas, L. Zhang, R. Duan, J. M. Yao, C. W. Tsai, W. W. Zhu, M. Cruces, G. Hobbs, C. C. Miao, J. R. Niu, M. D. Filipovic, S. Q. Zhu

Abstract: The environment of actively repeating fast radio bursts (FRBs) has been shown to be complex and varying. The recently localized FRB 20190520B is extremely active, has the largest confirmed host dispersion measure, and is only the second FRB source associated with a compact, persistent radio source (PRS). The main tracer of the magneto-ionic environments is the rotation measure (RM), a path-integra… ▽ More The environment of actively repeating fast radio bursts (FRBs) has been shown to be complex and varying. The recently localized FRB 20190520B is extremely active, has the largest confirmed host dispersion measure, and is only the second FRB source associated with a compact, persistent radio source (PRS). The main tracer of the magneto-ionic environments is the rotation measure (RM), a path-integral of the line-of-sight component of magnetic field strength (B) and electron density, which does not allow a direct probe of the B-field configuration. Here we report direct evidence for a B-field reversal based on the observed sign change and extreme variation of FRB 20190520B's RM, which changed from $\sim10000$ rad m$^{-2}$ to $\sim-16000$ rad m$^{-2}$ between June 2021 and January 2022. Such extreme RM reversal has never been observed before in any FRB nor in any astronomical object. The implied short-term change of the B-field configuration in or around the FRB could be due to the vicinity of massive black holes, or a magnetized companion star in binary systems, or a young supernova remnant along the line of sight. △ Less

Submitted 11 May, 2023; v1 submitted 15 March, 2022; originally announced March 2022.

Comments: It was merged into arXiv:2202.11112. Science380,599-603(2023). DOI:10.1126/science.abo6526

New infrared spectral indices of luminous cold stars: from early K to M-types

Authors: Maria Messineo, Donald F. Figer, Rolf-Peter Kudritzki, Qingfeng Zhu, Karl M. Menten, Valentin D. Ivanov, C. -H. Rosie Chen

Abstract: We present infrared spectral indices (1.0-2.3 um) of Galactic late-type giants and red supergiants (RSGs). We used existing and new spectra obtained at resolution power R=2000 with SpeX on the IRTF telescope. While a large CO equivalent width (EW), at 2.29 um ([CO, 2.29]>45 AA) is a typical signature of RSGs later than spectral type M0, [CO] of K-type RSGs and giants are similar. In the [CO, 2.29]… ▽ More We present infrared spectral indices (1.0-2.3 um) of Galactic late-type giants and red supergiants (RSGs). We used existing and new spectra obtained at resolution power R=2000 with SpeX on the IRTF telescope. While a large CO equivalent width (EW), at 2.29 um ([CO, 2.29]>45 AA) is a typical signature of RSGs later than spectral type M0, [CO] of K-type RSGs and giants are similar. In the [CO, 2.29] versus [Mg I, 1.71] diagram, RSGs of all spectral types can be distinguished from red giants, because the Mg I line weakens with increasing temperature and decreasing gravity. We find several lines that vary with luminosity, but not temperature: Si I (1.59 um), Sr (1.033 um), Fe+Cr+Si+CN (1.16 um), Fe+Ti (1.185 um), Fe+Ti (1.196 um), Ti+Ca (1.28 um), and Mn (1.29 um). Good markers of CN enhancement are the Fe+Si+CN line at 1.087 um and CN line at 1.093 um. Using these lines, at the resolution of SpeX, it is possible to separate RSGs and giants. Contaminant O-rich Mira and S-type AGBs are recognized by strong molecular features due to water vapor features, TiO band heads, and/or ZrO absorption. Among the 42 candidate RSGs that we observed, all but one were found to be late-types. 21 have EWs consistent with those of RSGs, 16 with those of O-rich Mira AGBs, and one with an S-type AGB. These infrared results open new, unexplored, potential for searches at low-resolution of RSGs in the highly obscured innermost regions of the Milky Way. △ Less

Submitted 8 July, 2021; originally announced July 2021.

Comments: 48 pages, 25 figures, accepted for publication in The Astronomical Journal

The Third Order Scalar Induced Gravitational Waves

Authors: Jing-Zhi Zhou, Xukun Zhang, Qing-Hua Zhu, Zhe Chang

Abstract: Since the gravitational waves were detected by LIGO and Virgo, it has been promising that lots of information about the primordial Universe could be learned by further observations on stochastic gravitational waves background. The studies on gravitational waves induced by primordial curvature perturbations are of great interest. The aim of this paper is to investigate the third order induced gravi… ▽ More Since the gravitational waves were detected by LIGO and Virgo, it has been promising that lots of information about the primordial Universe could be learned by further observations on stochastic gravitational waves background. The studies on gravitational waves induced by primordial curvature perturbations are of great interest. The aim of this paper is to investigate the third order induced gravitational waves. Based on the theory of cosmological perturbations, the first order scalar induces the second order scalar, vector and tensor perturbations. At the next iteration, the first order scalar, the second order scalar, vector and tensor perturbations all induce the third order tensor perturbations. We present the energy density spectrum of the third order gravitational waves for a monochromatic primordial power spectrum. The shape of the energy density spectrum of the third order gravitational waves is different from that of the second order scalar induced gravitational waves. And it is found that the third order gravitational waves sourced by the second order scalar perturbations dominate the energy density spectrum. △ Less

Submitted 9 May, 2022; v1 submitted 3 June, 2021; originally announced June 2021.

Comments: 33 pages, 4 figures

Detailed analysis of the poorly studied northern open cluster NGC 1348 using multi-color photometry and GAIA EDR3 astrometry

Authors: D. Bisht, Qingfeng Zhu, W. H. Elsanhoury, Devesh P. Sariya, Geeta Rangwal, R. K. S. Yadav, Alok Durgapal, Ing-Guey Jiang

Abstract: The membership determination for open clusters in noisy environments of the Milky Way is still an open problem. In this paper, our main aim is provide the membership probability of stars using proper motions and parallax values of stars using Gaia EDR3 astrometry. Apart from the Gaia astrometry, we have also used other photometric data sets like UKIDSS, WISE, APASS and Pan-STARRS1 in order to unde… ▽ More The membership determination for open clusters in noisy environments of the Milky Way is still an open problem. In this paper, our main aim is provide the membership probability of stars using proper motions and parallax values of stars using Gaia EDR3 astrometry. Apart from the Gaia astrometry, we have also used other photometric data sets like UKIDSS, WISE, APASS and Pan-STARRS1 in order to understand cluster properties from optical to mid-infrared regions. We selected 438 likely members with membership probability higher than $50\%$ and G$\le$20 mag. We obtained the mean value of proper motion as $μ_{x}=1.27\pm0.001$ and $μ_{y}=-0.73\pm0.002$ mas yr$^{-1}$. The cluster's radius is determined as 7.5 arcmin (5.67 pc) using radial density profile. Our analysis suggests that NGC 1348 is located at a distance of $2.6\pm0.05$ kpc. The mass function slope is found to be $1.30\pm0.18$ in the mass range 1.0$-$4.1 $M_\odot$, which is in fair agreement with Salpeter's value within the 1$σ$ uncertainty. The present study validates that NGC 1348 is a dynamically relaxed cluster. We computed the apex coordinates $(A, D)$ for NGC 1348 as $(A_\circ, D_\circ)$ = $(-23^{\textrm{o}}.815\pm 0^{\textrm{o}}.135$, $-22^{\textrm{o}}.228\pm 0^{\textrm{o}}.105)$. In addition, calculations of the velocity ellipsoid parameters (VEPs), matrix elements $μ_{ij}$, direction cosines ($l_j$, $m_j$, $n_j$) and the Galactic longitude of the vertex have been also conducted in this analysis. △ Less

Submitted 2 April, 2021; originally announced April 2021.

Comments: 14 Figures, 4 Tables, Accepted for publication by the Publication of the Astronomical Society of Japan (PASJ)

Multi-colour photometry and Gaia EDR3 astrometry of two couples of binary clusters (NGC 5617 and Trumpler 22) and (NGC 3293 and NGC 3324)

Authors: D. Bisht, Qingfeng Zhu, R. K. S. Yadav, Shashikiran Ganesh, Geeta Rangwal, Alok Durgapal, Devesh P. Sariya, Ing-Guey Jiang

Abstract: This paper presents a comprehensive analysis of two pairs of binary clusters (NGC 5617 and Trumpler 22) and (NGC 3293 and NGC 3324) located in the fourth quadrant of our Galaxy. For this purpose we use different data taken from VVV survey, WISE, VPHAS, APASS, GLIMPSE along with Gaia~EDR3 astrometric data. We identified 584, 429, 692 and 273 most probable cluster members with membership probability… ▽ More This paper presents a comprehensive analysis of two pairs of binary clusters (NGC 5617 and Trumpler 22) and (NGC 3293 and NGC 3324) located in the fourth quadrant of our Galaxy. For this purpose we use different data taken from VVV survey, WISE, VPHAS, APASS, GLIMPSE along with Gaia~EDR3 astrometric data. We identified 584, 429, 692 and 273 most probable cluster members with membership probability higher than 80 % towards the region of clusters NGC 5617, Trumpler 22, NGC 3293 and NGC 3324. We estimated the value of R as ~ 3.1 for clusters NGC 5617 and Trumpler 22, which indicates normal extinction law. The value of R ~ 3.8 and 1.9 represent the abnormal extinction law towards the clusters NGC 3293 and NGC 3324. Our Kinematical analysis show that all these clusters have circular orbits. Ages are found to be 90\pm10 and 12\pm3 Myr for the cluster pairs (NGC 5617 and Trumpler 22) and (NGC 3293 and NGC 3324), respectively. The distances of 2.43\pm0.08, 2.64\pm0.07, 2.59\pm0.1 and 2.80\pm0.2 kpc estimated using parallax are alike to the values calculated by using the distance modulus. We have also identified 18 and 44 young stellar object candidates present in NGC 5617 and Trumpler 22, respectively. Mass function slopes are found to be in fair agreement with the Salpeter's value. The dynamical study of these objects shows a lack of faint stars in their inner regions, which leads to the mass-segregation effect. Our study indicates that NGC 5617 and Trumpler 22 are dynamically relaxed but the other pair of clusters are not. Orbital alongwith the physical parameters show that the clusters in both pairs are physically connected. △ Less

Submitted 8 March, 2021; originally announced March 2021.

Comments: 21 pages, 20 figures, 7 tables, the article has been accepted for the publication MNRAS

A deep study of an intermediate age open cluster SAI 35 (Juchert 20) using ground based imaging and Gaia EDR3 astrometry

Authors: D. Bisht, Qingfeng Zhu, R. K. S. Yadav, Geeta Rangwal, Alok Durgapal, Devesh P. Sariya, Ing-Guey Jiang

Abstract: We present CCD $UBVI$ photometric study of poorly studied intermediate age open cluster SAI 35 (Juchert 20) for the first time. To accomplish this study, we also used LAMOST DR5, 2MASS, and Gaia EDR3 databases. We identified 214 most probable cluster members with membership probability higher than 50%. The mean proper motion of the cluster is found as μ_αcosδ=1.10 \pm 0.01 and μ_δ=-1.66 \pm 0.01 m… ▽ More We present CCD $UBVI$ photometric study of poorly studied intermediate age open cluster SAI 35 (Juchert 20) for the first time. To accomplish this study, we also used LAMOST DR5, 2MASS, and Gaia EDR3 databases. We identified 214 most probable cluster members with membership probability higher than 50%. The mean proper motion of the cluster is found as μ_αcosδ=1.10 \pm 0.01 and μ_δ=-1.66 \pm 0.01 mas/yr. We find the normal interstellar extinction law using the various two-color diagrams. The age, distance, reddening, and radial velocity of the cluster are estimated to be 360 \pm 40 Myr, 2.9 \pm 0.15 kpc, 0.72 \pm 0.05 mag and -91.62 \pm 6.39 km/sec. The overall mass function slope for main-sequence stars is found to be 1.49\pm0.16 within the mass range 1.1-3.1 M_\odot, which is in agreement with Salpeter's value within uncertainty. The present study demonstrates that SAI 35 is a dynamically relaxed. Galactic orbital parameters are determined using Galactic potential models. We found that this object follows a circular path around the Galactic center. △ Less

Submitted 7 February, 2021; originally announced February 2021.

Comments: 23 figures, 9 tables, This manuscript is accepted for publication in Astronomical journal

Fast Multipole Methods for $N$-body Simulations of Collisional Star Systems

Authors: Diptajyoti Mukherjee, Qirong Zhu, Hy Trac, Carl L. Rodriguez

Abstract: Direct $N$-body simulations of star clusters are accurate but expensive, largely due to the numerous $\mathcal{O} (N^2)$ pairwise force calculations. To solve the post-million-body problem, it will be necessary to use approximate force solvers, such as tree codes. In this work, we adapt a tree-based, optimized Fast Multipole Method (FMM) to the collisional $N$-body problem. The use of a rotation-a… ▽ More Direct $N$-body simulations of star clusters are accurate but expensive, largely due to the numerous $\mathcal{O} (N^2)$ pairwise force calculations. To solve the post-million-body problem, it will be necessary to use approximate force solvers, such as tree codes. In this work, we adapt a tree-based, optimized Fast Multipole Method (FMM) to the collisional $N$-body problem. The use of a rotation-accelerated translation operator and an error-controlled cell opening criterion leads to a code that can be tuned to arbitrary accuracy. We demonstrate that our code, Taichi, can be as accurate as direct summation when $N> 10^4$. This opens up the possibility of performing large-$N$, star-by-star simulations of massive stellar clusters, and would permit large parameter space studies that would require years with the current generation of direct summation codes. Using a series of tests and idealized models, we show that Taichi can accurately model collisional effects, such as dynamical friction and the core-collapse time of idealized clusters, producing results in strong agreement with benchmarks from other collisional codes such as NBODY6++GPU or PeTar. Parallelized using OpenMP and AVX, Taichi is demonstrated to be more efficient than other CPU-based direct $N$-body codes for simulating large systems. With future improvements to the handling of close encounters and binary evolution, we clearly demonstrate the potential of an optimized FMM for the modeling of collisional stellar systems, opening the door to accurate simulations of massive globular clusters, super star clusters, and even galactic nuclei. △ Less

Submitted 22 August, 2021; v1 submitted 3 December, 2020; originally announced December 2020.

Comments: Published in ApJ. This is the updated version of the paper after revisions. Updated title

The Formation of the First Quasars. I. The Black Hole Seeds, Accretion and Feedback Models

Authors: Qirong Zhu, Yuexing Li, Yiting Li, Moupiya Maji, Hidenobu Yajima, Raffaella Schneider, Lars Hernquist

Abstract: Supermassive black holes (SMBHs) of $\sim 10^9\, M_\odot$ are generally believed to be the central engines of the luminous quasars observed at $z\gtrsim6$, but their astrophysical origin remains elusive. The $z\gtrsim$ quasars reside in rare density peaks, which poses several challenges to uniform hydrodynamic simulations. To investigate the formation of these distant quasars, we perform a suite o… ▽ More Supermassive black holes (SMBHs) of $\sim 10^9\, M_\odot$ are generally believed to be the central engines of the luminous quasars observed at $z\gtrsim6$, but their astrophysical origin remains elusive. The $z\gtrsim$ quasars reside in rare density peaks, which poses several challenges to uniform hydrodynamic simulations. To investigate the formation of these distant quasars, we perform a suite of zoom-in simulations on a favorable halo, with a mass of $\sim 10^{13}\, M_\odot$ at $z = 6$ and a history of multiple major mergers, ideal for BH growth. We test BH seeds of $10 - 10^6\, M_\odot$, and various accretion and feedback models, including thin-disk and slim-disk accretion. We find, contrary to previous studies, that light seeds of $\lesssim 10^3\, M_\odot$ fail to grow to $10^8\, M_\odot$ by $z\sim 6$ even with super-critical accretion; that the hyper-Eddington mode leads to lower accretion rates than the Eddington-limited case due to stronger feedback, resulting in significantly smaller BHs by two orders of magnitude; and that while the super-critical model boosts the growth of low-spin BHs, for high-spin BHs the mass may be reduced due to increased radiative feedback. Our simulations show that the first $10^8 - 10^9\, M_\odot$ SMBHs may grow from heavy seeds of $\gtrsim 10^4\, M_\odot$ via Eddington-limited or mild super-critical accretion facilitated by gas-rich mergers and self-regulated by feedback, and they co-evolve with their host galaxies, producing bright quasars such as those at $z\sim$6 and ULAS J1342+0928, currently the most distant quasar at z = 7.54. △ Less

Submitted 2 December, 2020; originally announced December 2020.

Comments: 21 pages, 15 figures, comments welcome

Spatial variation of the fine-structure constant and Hubble's law in anisotropic coordinate of Friedmann-Lemaitre-Robertson-Walker space-time

Authors: Zhe Chang, Qing-Hua Zhu

Abstract: Recent updated results of quasar spectra suggested a 3.9$σ$ significance of spatial variation of the fine-structure constant. Theoretically, it is important to examine whether the fine-structure constant, as a fundamental constant in quantum theory, is possible varying with space and time. In this paper, we explore the possibility that spatial variation of the fine-structure constant could be comp… ▽ More Recent updated results of quasar spectra suggested a 3.9$σ$ significance of spatial variation of the fine-structure constant. Theoretically, it is important to examine whether the fine-structure constant, as a fundamental constant in quantum theory, is possible varying with space and time. In this paper, we explore the possibility that spatial variation of the fine-structure constant could be compatible with Einstein's general relativity. Namely, the spatially dependent fine-structure constant in the Universe could be originated in different values of the speed of light in separate local frames that are far away from us, since we have known that light rays must be bending in the present of gravity or non-inertial motions. In addition, to learn more about the anisotropic coordinate of FLRW space-time, we also study luminosity distance-redshift relation. It is found that there is a dipole structure in high redshift regime, while in low redshift regime, there is not such dipole. △ Less

Submitted 16 November, 2020; originally announced November 2020.

Comments: 22 pages, 9 figures

Searching for shocks in high-mass starless clump candidates

Authors: Feng-Yao Zhu, Jun-Zhi Wang, Tie Liu, Kee-Tae Kim, Qing-Feng Zhu, Fei Li

Abstract: In order to search for shocks in the very early stage of star formation, we performed single-point surveys of SiO J=1-0, 2-1 and 3-2 lines and the H$_2$CO $2_{12}-1_{11}$ line toward a sample of 100 high-mass starless clump candidates (SCCs) by using the Korean VLBI Network (KVN) 21-m radio telescopes. The detection rates of the SiO J=1-0, 2-1, 3-2 lines and the H$_2$CO line are $31.0\%$,… ▽ More In order to search for shocks in the very early stage of star formation, we performed single-point surveys of SiO J=1-0, 2-1 and 3-2 lines and the H$_2$CO $2_{12}-1_{11}$ line toward a sample of 100 high-mass starless clump candidates (SCCs) by using the Korean VLBI Network (KVN) 21-m radio telescopes. The detection rates of the SiO J=1-0, 2-1, 3-2 lines and the H$_2$CO line are $31.0\%$, $31.0\%$, $19.5\%$ and $93.0\%$, respectively. Shocks seem to be common in this stage of massive star formation. The widths of the observed SiO lines (full width at zero power (FWZP)) range from 3.4 to 55.1 km s$^{-1}$. A significant fraction ($\sim29\%$) of the detected SiO spectra have broad line widths (FWZP $>20~km~s^{-1}$), which are very likely associated with fast shocks driven by protostellar outflows. This result suggests that about one third of the SiO-detected SCCs are not really starless but protostellar. On the other hand, about 40$\%$ of the detected SiO spectra show narrow line widths (FWZP<10 $km~s^{-1}$) probably associated with low-velocity shocks which are not necessarily protostellar in origin. The estimated SiO column densities are mostly $0.31-4.32\times10^{12}~cm^{-2}$. Comparing the SiO column densities derived from SiO J=1-0 and 2-1 lines, we suggest that the SiO molecules in the SCCs may be in the non-LTE condition. The SiO abundances to H$_2$ are usually $0.20-10.92\times10^{-10}$. △ Less

Submitted 5 October, 2020; originally announced October 2020.

Comments: 19 pages, 13 figures

On the Gauge Invariance of Scalar Induced Gravitational Waves: Gauge Fixings Considered

Authors: Zhe Chang, Sai Wang, Qing-Hua Zhu

Abstract: The energy density spectrum is an observable of gravitational waves. Divergence has appeared in the energy density spectra of the scalar induced gravitational waves for different gauge fixings. To resolve the discrepancy, we investigate the gauge invariance of the scalar induced gravitational waves. It is shown that the gauge invariant induced gravitational waves can be obtained by subtracting the… ▽ More The energy density spectrum is an observable of gravitational waves. Divergence has appeared in the energy density spectra of the scalar induced gravitational waves for different gauge fixings. To resolve the discrepancy, we investigate the gauge invariance of the scalar induced gravitational waves. It is shown that the gauge invariant induced gravitational waves can be obtained by subtracting the fictitious tensor perturbations via introducing the counter term composed of the first order scalar perturbations. The kernel function uniquely determines the energy density spectrum of the scalar induced gravitational waves. We explicitly calculate the gauge invariant kernel functions in the Newtonian gauge and the uniform density gauge, respectively. The discrepancy between the energy density spectra upon the Newtonian gauge and the uniform density gauge is shown to be eliminated in the gauge invariant framework. In fact, the gauge invariant approach is also available to other kinds of gauge fixings. △ Less

Submitted 4 October, 2020; originally announced October 2020.

Comments: 9 pages, 5 figures, sequal to arXiv: 2009.11994 [gr-qc]

Gauge Invariant Second Order Gravitational Waves

Authors: Zhe Chang, Sai Wang, Qing-Hua Zhu

Abstract: We investigate the gauge invariance of the second order gravitational waves induced by the first order scalar perturbations by following the Lie derivative method. It is shown explicitly that the second order gravitational waves are gauge invariant in the synchronous frame. In the gauge invariant framework, we derive the equation of motion of the second order gravitational waves and show that the… ▽ More We investigate the gauge invariance of the second order gravitational waves induced by the first order scalar perturbations by following the Lie derivative method. It is shown explicitly that the second order gravitational waves are gauge invariant in the synchronous frame. In the gauge invariant framework, we derive the equation of motion of the second order gravitational waves and show that the second order gravitational waves are sourced from the first order scalar perturbations described well in the gauge invariant Newtonian frame. Since the observables of gravitational waves are measured in the synchronous frame, we define the energy density spectrum of the second order gravitational waves in terms of the gauge invariant synchronous variables. This way guarantees no fictitious tensor perturbations. It is shown that the gauge invariant energy density spectrum of the second order gravitational waves coincides with the one in the Newtonian gauge. △ Less

Submitted 29 September, 2020; v1 submitted 24 September, 2020; originally announced September 2020.

Comments: Major revisions. We further analyze the relation between the gauge invariant synchronous gravitational waves and the gravitational waves fixed in the Newtonian gauge. In this way, we could explain why one obtained correctly the energy density spectrum from the Newtonian gauge in previous studies (our method is available to other gauge fixings). Additional citations are added. Typos are corrected

Note on gauge invariance of second order cosmological perturbations

Authors: Zhe Chang, Sai Wang, Qing-Hua Zhu

Abstract: We study the gauge invariant cosmological perturbations up to second order. We show that there are infinite families of gauge invariant variables at both of the first and second orders. The conversion formulae among different families are shown to be described by a finite number of bases that are gauge invariant. For the second order cosmological perturbations induced by the first order scalar per… ▽ More We study the gauge invariant cosmological perturbations up to second order. We show that there are infinite families of gauge invariant variables at both of the first and second orders. The conversion formulae among different families are shown to be described by a finite number of bases that are gauge invariant. For the second order cosmological perturbations induced by the first order scalar perturbations, we explicitly represent the equations of motion of them in terms of the gauge invariant Newtonian, synchronous and hybrid variables, respectively. △ Less

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

Comments: v1: 31 pages. v2: Publised in Chinese Physics C

Journal ref: Chinese Physics C, Vol. 45. No. 9 (2021) 095101

An investigation of poorly studied open cluster NGC 4337 using multi-color photometric and Gaia DR2 astrometric data

Authors: D. Bisht, W. H. Elsanhoury, Qingfeng Zhu, Devesh P. Sariya, R. K. S. Yadav, Geeta Rangwal, Alok Durgapal, Ing-Guey Jiang

Abstract: We present a comprehensive analysis (photometric and kinematical) of poorly studied open cluster NGC 4337 using 2MASS, WISE, APASS, and Gaia~DR2 database. By determining the membership probabilities of stars, we identified 624 most probable members with membership probability higher than $50\%$ by using proper motion and parallax data taken from Gaia~DR2. The mean proper motion of the cluster is o… ▽ More We present a comprehensive analysis (photometric and kinematical) of poorly studied open cluster NGC 4337 using 2MASS, WISE, APASS, and Gaia~DR2 database. By determining the membership probabilities of stars, we identified 624 most probable members with membership probability higher than $50\%$ by using proper motion and parallax data taken from Gaia~DR2. The mean proper motion of the cluster is obtained as $μ_{x}=-8.83\pm0.01$ and $μ_{y}=1.49\pm0.006$ mas yr$^{-1}$. We find the normal interstellar extinction towards the cluster region. The radial distribution of members provides a cluster radius of 7.75 arcmin (5.63 pc). The estimated age of $1600\pm180$ Myr indicates that NGC 4337 is an old open cluster with a bunch of red giant stars. The overall mass function slope for main-sequence stars is found as $1.46\pm0.18$ within the mass range 0.75$-$2.0 $M_\odot$, which is in fair agreement with Salpeter's value (x=1.35) within uncertainty. The present study demonstrates that NGC 4337 is a dynamically relaxed open cluster. Using the Galactic potential model, Galactic orbits are obtained for NGC 4337. We found that this object follows a circular path around the Galactic center. Under the kinematical analysis, we compute the apex coordinates $(A, D)$ by using two methods: (i) the classical convergent point method and (ii) the AD-diagram method. The obtained coordinates are: $(A_{conv}, D_{conv})$ = (96$^{\textrm{o}}$.27 $\pm$ 0$^{\textrm{o}}$.10, 13$^{\textrm{o}}$.14 $\pm$ 0$^{\textrm{o}}$.27) $\&$ $(A_\circ, D_\circ)$ = (100$^{\textrm{o}}$.282 $\pm$ 0$^{\textrm{o}}$.10, 9$^{\textrm{o}}$.577 $\pm$ 0$^{\textrm{o}}$.323) respectively. We also computed the Velocity Ellipsoid Parameters (VEPs), matrix elements ($μ_{ij}$), direction cosines ($l_j$, $m_j$, $n_j$) and the Galactic longitude of the vertex ($l_2$). △ Less

Submitted 24 June, 2020; originally announced June 2020.

Comments: 24 pages, 16 figures and 4 table. The article has been accepted for publication in AJ

A Comprehensive study of open clusters Czernik 14, Haffner 14, Haffner 17 and King 10 using multicolour photometry and Gaia DR2 astrometry

Authors: D. Bisht, Qingfeng Zhu, R. K. S. Yadav, Alok Durgapal, Geeta Rangwal

Abstract: This paper presents an investigation on the four open clusters Czernik 14, Haffner 14, Haffner 17 and King 10 located near the Perseus arm of Milky Way Galaxy using Gaia DR2, 2MASS, WISE, APASS and Pan-STARRS1 data sets. We find normal interstellar extinction in twelve photometric bands for these clusters. Likely cluster members are identified as 225, 353, 350 and 395 for Czernik 14, Haffner 14, H… ▽ More This paper presents an investigation on the four open clusters Czernik 14, Haffner 14, Haffner 17 and King 10 located near the Perseus arm of Milky Way Galaxy using Gaia DR2, 2MASS, WISE, APASS and Pan-STARRS1 data sets. We find normal interstellar extinction in twelve photometric bands for these clusters. Likely cluster members are identified as 225, 353, 350 and 395 for Czernik 14, Haffner 14, Haffner 17 and King 10, respectively by using Gaia DR2 proper motion data. Radii are determined as 3.5, 3.7, 6.2 and 5.7 arcmin for Czernik 14, Haffner 14, Haffner 17 and King 10 respectively. Mean proper motions in RA and DEC are estimated as (-0.42 \pm 0.02, -0.38 \pm 0.01), (-1.82 \pm 0.009, 1.73 \pm 0.008), (-1.17 \pm 0.007, 1.88 \pm 0.006) and (-2.75 \pm 0.008, -2.04 \pm 0.006) mas/yr for Czernik 14, Haffner 14, Haffner 17 and King 10, respectively. The comparison of observed CMDs with solar metallicity isochrones leads to an age of 570 \pm 60, 320 \pm 35, 90 \pm 10 and 45 \pm 5 Myr for these clusters. The distances 2.9 \pm 0.1, 4.8 \pm 0.4, 3.6 \pm 0.1 and 3.8 \pm 0.1 kpc determined using parallax are comparable with the values derived by the isochrone fitting method. Mass function slopes are found to be in good agreement with the Salpeter value. The total masses are derived as 348, 595, 763 and 1088 solar mass for the clusters Czernik 14, Haffner 14, Haffner 17 and King 10, respectively. Evidence for the existence of mass-segregation effect is observed in each cluster. Using the Galactic potential model, Galactic orbits are derived for the clusters. The present study indicates that all clusters under study fallow a circular path around the Galactic center. △ Less

Submitted 5 March, 2020; originally announced March 2020.

Comments: 17 figures, 7 tables. This article is accepted for publication in MNRAS

ART$^2$: A 3D Parallel Multi-wavelength Radiative Transfer Code for Continuum and Atomic and Molecular Lines

Authors: Yuexing Li, Ming F. Gu, Hidenobu Yajima, Qirong Zhu, Moupiya Maji

Abstract: ART$^2$ is a 3D multi-wavelength Monte Carlo radiative transfer (RT) code that couples continuum and emission lines to track the propagation of photons and their interactions with the interstellar medium (ISM). The previous version of ART$^2$, which included continuum and Ly$α$ line, has been extensively applied to hydrodynamics simulations to study multi-band properties of galaxies and ISM. Here,… ▽ More ART$^2$ is a 3D multi-wavelength Monte Carlo radiative transfer (RT) code that couples continuum and emission lines to track the propagation of photons and their interactions with the interstellar medium (ISM). The previous version of ART$^2$, which included continuum and Ly$α$ line, has been extensively applied to hydrodynamics simulations to study multi-band properties of galaxies and ISM. Here, we describe new implementations of non-local thermodynamic equilibrium RT of molecular and atomic fine structure emission lines, and the parallelization of the code using a number of novel methods. The new ART$^2$ can efficiently and self-consistently produce a full spectrum that includes both continuum and lines such as [CII], [NII], [OIII], Ly$α$, and CO. These essential features, together with the multi-phase ISM model and the adaptive grid, make ART$^2$ a multi-purpose code to study multi-wavelength properties of a wide range of astrophysical systems from planetary disks to large-scale structures. To demonstrate the capability of the new ART$^2$, we applied it to two hydrodynamics simulations: the zoom-in Milky Way Simulation to obtain panchromatic properties of individual galaxies, and the large-scale IllustrisTNG100 Simulation to obtain global properties such as the line intensity mappings. These products are essential for a broad array of studies such as the correlations between physical and panchromatic properties and their evolution. By enabling direct comparison between numerical simulations and multi-band observations, ART$^2$ provides a crucial theoretical framework for the interpretations of existing observations, the plan for future surveys, and the synergy between multi-band galaxy surveys and line intensity mappings. Therefore, ART$^2$ is a powerful and versatile tool to bridge the gap between theories and observations of cosmic structures. △ Less

Submitted 12 March, 2020; v1 submitted 29 January, 2020; originally announced January 2020.

Comments: Accepted for publication in MNRAS; 20 pages, 26 figures

Determining electron temperature and density in a H II region by using the relative strengths of hydrogen radio recombination lines

Authors: F. -Y. Zhu, Q. -F. Zhu, J. -Z. Wang, J. -S Zhang

Abstract: We have introduced a new method of estimating the electron temperature and density of H II regions by using single dish observations. In this method, multiple hydrogen radio recombination lines of different bands are computed under the assumption of low optical depth. We use evolutionary hydrodynamical models of H II regions to model hydrogen recombination line emission from a variety of H II regi… ▽ More We have introduced a new method of estimating the electron temperature and density of H II regions by using single dish observations. In this method, multiple hydrogen radio recombination lines of different bands are computed under the assumption of low optical depth. We use evolutionary hydrodynamical models of H II regions to model hydrogen recombination line emission from a variety of H II regions and assess the reliability of the method. According to the simulated results, the error of the estimated temperature is commonly < 13%, and that of the estimated density is < 25% for a < 1% uncertainty of the observed line fluxes. A reasonable estimated value of electron density can be achieved if the uncertainty of the line fluxes are lower than 3%. In addition, the estimated values are more representative of the properties in the relatively high-density region if the gas density gradient is present in the H II region. Our method can be independent of the radio continuum observations. But the accuracy will be improved if a line-to-continuum ratio at millimeter wavelengths is added to the estimation. Our method provides a way to measure the temperature and density in ionized regions without interferometers. △ Less

Submitted 13 November, 2019; originally announced November 2019.

Comments: 17 pages, 12 figures

Journal ref: The Astrophysical Journal, Volume 881, Issue 1, article id. 14, 14 pp. (2019)

Laboratory formation and photo-chemistry of fullerene/anthracene cluster cations

Authors: Junfeng Zhen, Weiwei Zhang, YuanYuan Yang, Qingfeng Zhu, Alexander G. G. M. Tielens

Abstract: Besides buckminsterfullerene (C60), other fullerenes and their derivatives may also reside in space. In this work, we study the formation and photo-dissociation processes of astronomically relevant fullerene/anthracene (C14H10) cluster cations in the gas phase. Experiments are carried out using a quadrupole ion trap (QIT) in combination with time-of-flight (TOF) mass spectrometry. The results show… ▽ More Besides buckminsterfullerene (C60), other fullerenes and their derivatives may also reside in space. In this work, we study the formation and photo-dissociation processes of astronomically relevant fullerene/anthracene (C14H10) cluster cations in the gas phase. Experiments are carried out using a quadrupole ion trap (QIT) in combination with time-of-flight (TOF) mass spectrometry. The results show that fullerene (C60, and C70)/anthracene (i.e., [(C14H10)nC60]+ and [(C14H10)nC70]+), fullerene (C56 and C58)/anthracene (i.e., [(C14H10)nC56]+ and [(C14H10)nC58]+) and fullerene (C66 and C68)/anthracene (i.e., [(C14H10)nC66]+ and [(C14H10)nC68]+) cluster cations, are formed in the gas phase through an ion-molecule reaction pathway. With irradiation, all the fullerene/anthracene cluster cations dissociate into mono$-$anthracene and fullerene species without dehydrogenation. The structure of newly formed fullerene/anthracene cluster cations and the bonding energy for these reaction pathways are investigated with quantum chemistry calculations. Our results provide a growth route towards large fullerene derivatives in a bottom-up process and insight in their photo-evolution behavior in the ISM, and clearly, when conditions are favorable, fullerene/PAH clusters can form efficiently. In addition, these clusters (from 80 to 154 atoms or ~ 2 nm in size) offer a good model for understanding the physical-chemical processes involved in the formation and evolution of carbon dust grains in space, and provide candidates of interest for the DIBs that could motivate spectroscopic studies. △ Less

Submitted 24 October, 2019; originally announced October 2019.

Comments: 10 pages, 8 figures, accepted

Journal ref: The Astrophysical Journal, 2019

Reference frame dependence of the local measurement of the Hubble constant

Authors: Zhe Chang, Qing-Hua Zhu

Abstract: Observations recently suggest $4.4σ$ discrepancy between Hubble constant measured locally and that inferred from cosmic microwave background in standard cosmology. Either standard cosmology or local measurement might expect something new. We investigate the possibility that the value of Hubble constants might be affected by observers' motional status in the local measurement. Using adapted coordin… ▽ More Observations recently suggest $4.4σ$ discrepancy between Hubble constant measured locally and that inferred from cosmic microwave background in standard cosmology. Either standard cosmology or local measurement might expect something new. We investigate the possibility that the value of Hubble constants might be affected by observers' motional status in the local measurement. Using adapted coordinate for geodesic observers in FLRW space-time and constraints inferred from observation of cosmic shear, we find that the motional status of reference frame could contribute to about $1.1 \pm 0.3\%$ discrepancy of these Hubble constants. Namely, the difference is $ΔH_0= 0.81\pm0.22$ km s$^{-1}$ Mpc$^{-1}$, if the locally measured Hubble constant is $ 74$ km s$^{-1}$ Mpc$^{-1}$. This effect seems not negligible as an uncertainty for local measurement of Hubble constant △ Less

Submitted 28 May, 2020; v1 submitted 29 September, 2019; originally announced September 2019.

Comments: v1: 12 pages; v2: 19 pages, 2 figures. A revised version for publishing in PRD

Journal ref: Phys. Rev. D 101, 103530 (2020)

Statistical analysis on X-ray flares from the nucleus and HST-1 knot in the M87 jet

Authors: Shenbang Yang, Dahai Yan, Benzhong Dai, Pengfei Zhang, Qianqian Zhu, Jiancheng Wang, Li Zhang

Abstract: The statistical properties of X-ray flares from two separate locations (nucleus and HST-1) in the M87 jet are investigated to reveal the physical origin of the flares. We analyse the archival \textit{Chandra} data for M87, and identify 14 flares in the nucleus and 9 flares in HST-1. The peak intensity ($I_{\rm{P}}$) and the flaring duration time ($T_{\rm{fl}}$) for each flare are obtained. It is f… ▽ More The statistical properties of X-ray flares from two separate locations (nucleus and HST-1) in the M87 jet are investigated to reveal the physical origin of the flares. We analyse the archival \textit{Chandra} data for M87, and identify 14 flares in the nucleus and 9 flares in HST-1. The peak intensity ($I_{\rm{P}}$) and the flaring duration time ($T_{\rm{fl}}$) for each flare are obtained. It is found that the distributions of both $I_{\rm{P}}$ and $T_{\rm{fl}}$ for the nucleus obey a power-law form with a similar index. A similar result is also obtained for HST-1, and no significant inconsistency between the nucleus and HST-1 is found for the indices. Similar to solar X-ray flares, the power-law distributions of the flare event parameters can be well explained by a self-organized criticality (SOC) system, which are triggered by magnetic reconnection. Our results suggest that the flares from nucleus and HST-1 are possibly triggered by magnetic reconnection process. The consistent indices for the distributions of $I_{\rm{P}}$ and $T_{\rm{fl}}$ in the CORE and HST-1 indicate that the dimensions of the energy dissipation of the magnetic reconnection is identical in the two regions. A strong correlation between the flares in the two regions also suggests a similar physical origin for the flares. △ Less

Submitted 16 August, 2019; originally announced August 2019.

Comments: 10 pages, 6 figures, accepted for publication in MNRAS

A SiO J = 5 - 4 Survey Toward Massive Star Formation Regions

Authors: Shanghuo Li, Junzhi Wang, Min Fang, Qizhou Zhang, Fei Li, Zhi-Yu Zhang, Juan Li, Qingfeng Zhu, Shaoshan Zeng

Abstract: We performed a survey in the SiO $J=5\rightarrow4$ line toward a sample of 199 Galactic massive star-forming regions at different evolutionary stages with the SMT 10 m and CSO 10.4 m telescopes. The sample consists of 44 infrared dark clouds (IRDCs), 86 protostellar candidates, and 69 young \HII\ regions. We detected SiO $J=5\rightarrow4$ line emission in 102 sources, with a detection rate of 57\%… ▽ More We performed a survey in the SiO $J=5\rightarrow4$ line toward a sample of 199 Galactic massive star-forming regions at different evolutionary stages with the SMT 10 m and CSO 10.4 m telescopes. The sample consists of 44 infrared dark clouds (IRDCs), 86 protostellar candidates, and 69 young \HII\ regions. We detected SiO $J=5\rightarrow4$ line emission in 102 sources, with a detection rate of 57\%, 37\%, and 65\% for IRDCs, protostellar candidates, and young \HII\ regions, respectively. We find both broad line with Full Widths at Zero Power (FWZP) $>$ 20 \kms and narrow line emissons of SiO in objects at various evolutionary stages, likely associated with high-velocity shocks and low-velocity shocks, respectively. The SiO luminosities do not show apparent differences among various evolutionary stages in our sample. We find no correlation between the SiO abundance and the luminosity-to-mass ratio, indicating that the SiO abundance does not vary significantly in regions at different evolutionary stages of star formation. △ Less

Submitted 20 May, 2019; originally announced May 2019.

Comments: 25 pages, 9 figures, 5 tables, accepted for publication in ApJ

The Potential of FAST in Detecting Celestial Hydroxyl Masers and Related Science Topics

Authors: J. S. Zhang, D. Li, J. Z. Wang, Q. F. Zhu, J. Li

Abstract: The Five Hundred Meter Aperture Spherical Radio Telescope (FAST) will make contributions to the study of Galactic and extragalactic masers. The telescope, now finished construction and commissioning in China, has an innovative design that leads to the highest sensitivity of any single dish telescope in the world. FAST's potential for OH megamaser research is discussed, including the sky density of… ▽ More The Five Hundred Meter Aperture Spherical Radio Telescope (FAST) will make contributions to the study of Galactic and extragalactic masers. The telescope, now finished construction and commissioning in China, has an innovative design that leads to the highest sensitivity of any single dish telescope in the world. FAST's potential for OH megamaser research is discussed, including the sky density of masers detectable in surveys. The scientific impact expected from FAST maser studies is also discussed. △ Less

Submitted 14 April, 2019; originally announced April 2019.

Comments: 9 pages, 2 figures. Published in Research in Astronomy and Astrophysics

Journal ref: Res. Astron. Astrophys. 19 022 (2019)

Search for gamma-ray emission from the Sun during solar minimum with the ARGO-YBJ experiment

Authors: B. Bartoli, P. Bernardini, X. J. Bi, Z. Cao, S. Catalanotti, S. Z. Chen, T. L. Chen, S. W. Cui, B. Z. Dai, A. D'Amone, Danzengluobu, I. De Mitri, B. D'Ettorre Piazzoli, T. Di Girolamo, G. Di Sciascio, C. F. Feng, Zhaoyang Feng, Zhenyong Feng, W. Gao, Q. B. Gou, Y. Q. Guo, H. H. He, Haibing Hu, Hongbo Hu, M. Iacovacci , et al. (48 additional authors not shown)

Abstract: The hadronic interaction of cosmic rays with solar atmosphere can produce high energy gamma rays. The gamma-ray luminosity is correlated both with the flux of primary cosmic rays and the intensity of the solar magnetic field. The gamma rays below 200 GeV have been observed by $Fermi$ without any evident energy cutoff. The bright gamma-ray flux above 100 GeV has been detected only during solar mini… ▽ More The hadronic interaction of cosmic rays with solar atmosphere can produce high energy gamma rays. The gamma-ray luminosity is correlated both with the flux of primary cosmic rays and the intensity of the solar magnetic field. The gamma rays below 200 GeV have been observed by $Fermi$ without any evident energy cutoff. The bright gamma-ray flux above 100 GeV has been detected only during solar minimum. The only available data in TeV range come from the HAWC observations, however outside the solar minimum. The ARGO-YBJ dataset has been used to search for sub-TeV/TeV gamma rays from the Sun during the solar minimum from 2008 to 2010, the same time period covered by the Fermi data. A suitable model containing the Sun shadow, solar disk emission and inverse-Compton emission has been developed, and the chi-square minimization method was used to quantitatively estimate the disk gamma-ray signal. The result shows that no significant gamma-ray signal is detected and upper limits to the gamma-ray flux at 0.3$-$7 TeV are set at 95\% confidence level. In the low energy range these limits are consistent with the extrapolation of the Fermi-LAT measurements taken during solar minimum and are compatible with a softening of the gamma-ray spectrum below 1 TeV. They provide also an experimental upper bound to any solar disk emission at TeV energies. Models of dark matter annihilation via long-lived mediators predicting gamma-ray fluxes > $10^{-7}$ GeV $cm^{-2}$ $s^{-1}$ below 1 TeV are ruled out by the ARGO-YBJ limits. △ Less

Submitted 14 January, 2019; originally announced January 2019.

Impact of Plasma Instability on Constraint of the Intergalactic Magnetic Field

Authors: Dahai Yan, Jianeng Zhou, Pengfei Zhang, Qianqian Zhu, Jiancheng Wang

Abstract: A relativistic electron-positron pair beam can be produced in the interaction of TeV photons from a blazar with the extragalactic background light (EBL). The relativistic $e^{\pm}$ pairs would loss energy through inverse-Compton scattering (ICS) photons of cosmic microwave background (CMB) or plasma instabilities. The dominant energy-loss process is under debate. Based on the assumption that the d… ▽ More A relativistic electron-positron pair beam can be produced in the interaction of TeV photons from a blazar with the extragalactic background light (EBL). The relativistic $e^{\pm}$ pairs would loss energy through inverse-Compton scattering (ICS) photons of cosmic microwave background (CMB) or plasma instabilities. The dominant energy-loss process is under debate. Based on the assumption that the dominant energy-loss process is ICS, the resulted cascade GeV radiation is usually used to constrain the intergalactic magnetic field (IGMF). Here, we include the energy-loss due to plasma oblique instability in the calculation of cascade gamma-ray flux, and investigate the impact of the plasma instability on the constraint of IGMF. The up-to-date GeV data and archival TeV data of the blazar 1ES 0229+200 are used. The results indicate that even if the oblique instability cooling is dominating over ICS cooling, the cascade flux could be still used to constrain the IGMF. It is found that with the ratio between the cooling rates of the oblique instability and the ICS varying from 0.1, 1 to 10, the lower limit of the IGMF putted by the cascade flux and the gamma-ray data changes from $8\times10^{-18}\ $G, $5\times10^{-18}\ $G to $10^{-18}\ $G. If the ratio between the two cooling rates is 30, the estimate of IGMF based on the cascade flux is invalid. △ Less

Submitted 6 November, 2018; v1 submitted 16 October, 2018; originally announced October 2018.

Comments: Accepted by ApJ

Design of Near Infrared Sky Brightness Monitor and Test Running at Ngari Observatory in Tibet

Authors: Qi-Jie Tang, Jian Wang, Shu-cheng Dong, Jin-ting Chen, Yi-hao Zhang, Feng-xin Jiang, Zhi-yue Wang, Ya-qi Chen, Ming-hao Jia, Jie Chen, Hong-fei Zhang, Qing-feng Zhu, Peng Jiang, Tuo Ji, Shao-hua Zhang, Yong-qiang Yao, Yun-he Zhou, Hong-shuai Wang, Peng Tang

Abstract: Tibet is known as the third pole of the earth, as high as the South Pole and North Pole. The Ngari (Ali) observatory in Tibet has the advantage of plenty of photometric night, low precipitable water vapor, high transmittance, good seeing. It is a good site, and promising to be one of the best place for infrared and submillimeter observations in the world. However, there is no data available for sk… ▽ More Tibet is known as the third pole of the earth, as high as the South Pole and North Pole. The Ngari (Ali) observatory in Tibet has the advantage of plenty of photometric night, low precipitable water vapor, high transmittance, good seeing. It is a good site, and promising to be one of the best place for infrared and submillimeter observations in the world. However, there is no data available for sky background brightness in such place, which restrict the astronomical development of the sites. In the near infrared band of J, H, Ks, a NIR sky brightness monitor (NISBM) is designed based on InGaAs photoelectric diode. By using the method of chopper modulation and digital lock-in amplifier processing, the SNR (Signal Noise Ratio), detectivity and the data acquisition speed of the device is greatly improved. For each band of J, H, Ks, an independent instrument is designed and calibrated in laboratory. The NISBM has been installed in Ngari observatory in July of 2017 and obtained the first data of NIR sky brightness at Ngari observatory. △ Less

Submitted 20 June, 2018; originally announced June 2018.

Comments: 19 pages,18 figures

Formation of a Malin 1 analogue in IllustrisTNG by stimulated accretion

Authors: Qirong Zhu, Dandan Xu, Massimo Gaspari, Vicente Rodriguez-Gomez, Dylan Nelson, Mark Vogelsberger, Paul Torrey, Annalisa Pillepich, Jolanta Zjupa, Rainer Weinberger, Federico Marinacci, Rüdiger Pakmor, Shy Genel, Yuexing Li, Volker Springel, Lars Hernquist

Abstract: The galaxy Malin 1 contains the largest stellar disk known but the formation mechanism of this structure has been elusive. In this paper, we report a Malin 1 analogue in the 100 Mpc IllustrisTNG simulation and describe its formation history. At redshift zero, this massive galaxy, having a maximum circular velocity $V_{\rm max}$ of 430 ${\rm km\ s^{-1}}$, contains a 100 kpc gas/stellar disk with mo… ▽ More The galaxy Malin 1 contains the largest stellar disk known but the formation mechanism of this structure has been elusive. In this paper, we report a Malin 1 analogue in the 100 Mpc IllustrisTNG simulation and describe its formation history. At redshift zero, this massive galaxy, having a maximum circular velocity $V_{\rm max}$ of 430 ${\rm km\ s^{-1}}$, contains a 100 kpc gas/stellar disk with morphology similar to Malin 1. The simulated galaxy reproduces well many observed features of Malin 1's vast disk, including its stellar ages, metallicities, and gas rotation curve. We trace the extended disk back in time and find that a large fraction of the cold gas at redshift zero originated from the cooling of hot halo gas, triggered by the merger of a pair of intruding galaxies. Our finding provides a novel way to form large galaxy disks as extreme as Malin 1 within the current galaxy formation framework. △ Less

Submitted 23 May, 2018; originally announced May 2018.

Comments: 5 page, 5 figures; resubmitted to MNRAS Letters after first referee report; Comments welcome. Animations of gas/stars can be found here at https://youtu.be/Ibl4Cyybw2Q and https://youtu.be/3aIpCu6e9dQ, or downloaded at http://www.tng-project.org/movies/tng/TNG_Malin_1_gas.mpeg and http://www.tng-project.org/movies/tng/TNG_Malin_1_star.mpeg

Galactic Cosmic-Ray Anisotropy in the Northern hemisphere from the ARGO-YBJ Experiment during 2008-2012

Authors: B. Bartoli, P. Bernardini, X. J. Bi, Z. Cao, S. Catalanotti, S. Z. Chen, T. L. Chen, S. W. Cui, B. Z. Dai, A. D'Amone, Danzengluobu, I. De Mitri, B. D'Ettorre Piazzoli, T. Di Girolamo, G. Di Sciascio, C. F. Feng, Z. Y. Feng, W. Gao, Q. B. Gou, Y. Q. Guo, H. H. He, Haibing Hu, Hongbo Hu, M. Iacovacci, R. Iuppa , et al. (47 additional authors not shown)

Abstract: This paper reports on the observation of the sidereal large-scale anisotropy of cosmic rays using data collected by the ARGO-YBJ experiment over 5 years (2008$-$2012). This analysis extends previous work limited to the period from 2008 January to 2009 December,near the minimum of solar activity between cycles 23 and 24.With the new data sample the period of solar cycle 24 from near minimum to maxi… ▽ More This paper reports on the observation of the sidereal large-scale anisotropy of cosmic rays using data collected by the ARGO-YBJ experiment over 5 years (2008$-$2012). This analysis extends previous work limited to the period from 2008 January to 2009 December,near the minimum of solar activity between cycles 23 and 24.With the new data sample the period of solar cycle 24 from near minimum to maximum is investigated. A new method is used to improve the energy reconstruction, allowing us to cover a much wider energy range, from 4 to 520 TeV. Below 100 TeV, the anisotropy is dominated by two wide regions, the so-called "tail-in" and "loss-cone" features. At higher energies, a dramatic change of the morphology is confirmed. The yearly time dependence of the anisotropy is investigated. Finally, no noticeable variation of cosmic-ray anisotropy with solar activity is observed for a median energy of 7 TeV. △ Less

Submitted 23 May, 2018; originally announced May 2018.