Showing 1–50 of 71 results for author: Yu, Q

CO Observations of Early-mid Stage Major-mergers in MaNGA Survey

Authors: Qingzheng Yu, Taotao Fang, Cong Kevin Xu, Shuai Feng, Siyi Feng, Yu Gao, Xue-Jian Jiang, Ute Lisenfeld

Abstract: We present a study of the molecular gas in early-mid stage major-mergers, with a sample of 43 major-merger galaxy pairs selected from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey and a control sample of 195 isolated galaxies selected from the xCOLD GASS survey. Adopting kinematic asymmetry as a new effective indicator to describe the merger stage, we aim to study the role… ▽ More We present a study of the molecular gas in early-mid stage major-mergers, with a sample of 43 major-merger galaxy pairs selected from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey and a control sample of 195 isolated galaxies selected from the xCOLD GASS survey. Adopting kinematic asymmetry as a new effective indicator to describe the merger stage, we aim to study the role of molecular gas in the merger-induced star formation enhancement along the merger sequence of galaxy pairs. We obtain the molecular gas properties from CO observations with the James Clerk Maxwell Telescope (JCMT), Institut de Radioastronomie Milimetrique (IRAM) 30-m telescope, and the MASCOT survey. Using these data, we investigate the differences in molecular gas fraction ($f_{\rm H_{2}}$), star formation rate (SFR), star formation efficiency (SFE), molecular-to-atomic gas ratio ($M_{\rm H_{2}}/M_{\rm HI}$), total gas fraction ($f_{\rm gas}$), and the star formation efficiency of total gas (${\rm SFE_{gas}}$) between the pair and control samples. In the full pair sample, our results suggest the $f_{\rm H_{2}}$ of paired galaxies is significantly enhanced, while the SFE is comparable to that of isolated galaxies. We detect significantly increased $f_{\rm H_{2}}$ and $M_{\rm H_{2}}/M_{\rm HI}$ in paired galaxies at the pericenter stage, indicating an accelerated transition from atomic gas to molecular gas due to interactions. Our results indicate that the elevation of $f_{\rm H_{2}}$ plays a major role in the enhancement of global SFR in paired galaxies at the pericenter stage, while the contribution of enhanced SFE in specific regions requires further explorations through spatially resolved observations of a larger sample spanning a wide range of merger stages. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Comments: 25 pages, 12 figures, 5 tables, accepted for publication in ApJS

Close Major-merger Pairs at $z=0$: Star-forming Galaxies with Pseudobulges

Authors: Chuan He, Cong Kevin Xu, Ute Lisenfeld, Y Sophia Dai, Taotao Fang, Jia-Sheng Huang, Wei Wang, Qingzheng Yu

Abstract: We present a study of star-forming galaxies (SFGs) with pseudobulges (bulges with Sérsic index $\rm n < 2$) in a local close major-merger galaxy pair sample (H-KPAIR). With data from new aperture photometries in the optical and near-infrared bands (aperture size of 7\;kpc) and from the literature, we find that the mean Age of central stellar populations in Spirals with pseudobulges is consistent w… ▽ More We present a study of star-forming galaxies (SFGs) with pseudobulges (bulges with Sérsic index $\rm n < 2$) in a local close major-merger galaxy pair sample (H-KPAIR). With data from new aperture photometries in the optical and near-infrared bands (aperture size of 7\;kpc) and from the literature, we find that the mean Age of central stellar populations in Spirals with pseudobulges is consistent with that of disky galaxies and is nearly constant against the bulge-to-total ratio (B/T). Paired Spirals have a slightly lower fraction of pure disk galaxies ($\rm B/T \leq 0.1$) than their counterparts in the control sample. Compared to SFGs with classical bulges, those with pseudobulges have a higher ($>2\;σ$) mean of specific star formation rate (sSFR) enhancement ($\rm sSFR_{enh} = 0.33\pm0.07$ vs $\rm sSFR_{enh} = 0.12\pm0.06$) and broader scatter (by $\sim 1$\;dex). The eight SFGs that have the highest $\rm sSFR_{enh}$ in the sample all have pseudobulges. A majority (69\%) of paired SFGs with strong enhancement (having sSFR more than 5 times the median of the control galaxies) have pseudobulges. The Spitzer data show that the pseudobulges in these galaxies are tightly linked to nuclear/circum-nuclear starbursts. Pseudobulge SFGs in S+S and in S+E pairs have significantly ($>3\;σ$) different sSFR enhancement, with the means of $\rm sSFR_{enh} = 0.45\pm0.08$ and $-0.04\pm0.11$, respectively. We find a decrease in the sSFR enhancements with the density of the environment for SFGs with pseudobulges. Since a high fraction (5/11) of pseudobulge SFGs in S+E pairs are in rich groups/clusters (local density $\rm N_{1Mpc} \geq 7$), the dense environment might be the cause for their low $\rm sSFR_{enh}$. △ Less

Submitted 28 March, 2024; originally announced March 2024.

Comments: Accepted for publication in RAA, ?? pages, 10 figures, 4 tables

Low- and High-velocity \ion{O}{6} in Milky Way-like Galaxies: the Role of Stellar Feedback

Authors: Zhijie Zhang, Xiaoxia Zhang, Hui Li, Taotao Fang, Qingzheng Yu, Yang Luo, Federico Marinacci, Laura V. Sales, Paul Torrey, Mark Vogelsberger

Abstract: Milky Way-type galaxies are surrounded by a warm-hot gaseous halo containing a considerable amount of baryons and metals. The kinematics and spatial distribution of highly-ionized ion species such as \ion{O}{6} can be significantly affected by supernova (SN) explosions and early (pre-SN) stellar feedback (e.g., stellar winds, radiation pressure). Here, we investigate effects of stellar feedback on… ▽ More Milky Way-type galaxies are surrounded by a warm-hot gaseous halo containing a considerable amount of baryons and metals. The kinematics and spatial distribution of highly-ionized ion species such as \ion{O}{6} can be significantly affected by supernova (SN) explosions and early (pre-SN) stellar feedback (e.g., stellar winds, radiation pressure). Here, we investigate effects of stellar feedback on \ion{O}{6} absorptions in Milky Way-like galaxies by analyzing the suites of high-resolution hydrodynamical simulations under the framework of {\it SMUGGLE}, a physically motivated subgrid interstellar medium and stellar feedback model for the moving-mesh code {\sc Arepo}. We find that the fiducial run with the full suite of stellar feedback and moderate star formation activities can reasonably reproduce Galactic \ion{O}{6} absorptions observed by space telescopes such as {\it FUSE}, including the scale height of low-velocity ($|v_{\rm LSR}|< 100\, \rm km~s^{-1}$) \ion{O}{6}, the column density $-$ line width relation for high-velocity ($100 \leq |v_{\rm LSR}|< 400\, \rm km~s^{-1}$) \ion{O}{6}, and the cumulative \ion{O}{6} column densities. In contrast, model variations with more intense star formation activities deviate from observations further. Additionally, we find that the run considering only SN feedback is in broad agreement with the observations, whereas in runs without SN feedback this agreement is absent, which indicates a dominant role of SN feedback in heating and accelerating interstellar \ion{O}{6}. This is consistent with the current picture that interstellar \ion{O}{6} is predominantly produced by collisional ionization where mechanical feedback can play a central role. In contrast, photoionization is negligible for \ion{O}{6} production due to the lack of high-energy ($\gtrsim114\ {\rm eV}$) photons required. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: 20 pages, 14 figures, ApJ accepted

Pulsar timing array detections of supermassive binary black holes: implications from the detected common process signal and beyond

Authors: Yunfeng Chen, Qingjuan Yu, Youjun Lu

Abstract: Pulsar timing arrays (PTAs) are anticipated to detect the stochastic gravitational wave background (GWB) from supermassive binary black holes (BBHs) as well as the gravitational waves from individual BBHs. Recently, a common process signal was reported by several PTAs. In this paper, we investigate the constraints on the BBH population model(s) by current PTA observations and further study the det… ▽ More Pulsar timing arrays (PTAs) are anticipated to detect the stochastic gravitational wave background (GWB) from supermassive binary black holes (BBHs) as well as the gravitational waves from individual BBHs. Recently, a common process signal was reported by several PTAs. In this paper, we investigate the constraints on the BBH population model(s) by current PTA observations and further study the detections of both the GWB and individual BBHs by current/future PTAs. We find that the MBH--host galaxy scaling relation, an important ingredient of the BBH population model, is required to either evolve significantly with redshift or have a normalization $\sim0.86-1.1$ dex higher than the empirical ones, if the GWB is the same as the common process signal. For both cases, the estimated detection probability for individual BBHs is too small for a positive detection by current PTAs. By involving either the constrained scaling relations or those empirical ones into the BBH population models, we estimate that the GWB may be detected with a signal-to-noise ratio $\gtrsim3$ by the PTAs based on the Five hundred meter Aperture Spherical radio Telescope (CPTA) and the Square Kilometer Array (SKAPTA) after $\sim2-3$ (or $\sim6-11$) years' observation, if it is the same as (or an order of magnitude lower than) the common process signal. The detection time of individual BBHs by CPTA and SKAPTA is close to that of the GWB detection. We show that the BBH population model can be strongly constrained by the number and property distributions of BBHs to be detected by future PTAs. △ Less

Submitted 19 June, 2023; originally announced June 2023.

Comments: 25 pages, 9 figures, submitted to ApJ on April 10, 2023, and revised after reviewer's report

The correlations of stellar tidal disruption rates with properties of massive black holes and their host galaxies

Authors: Yunfeng Chen, Qingjuan Yu, Youjun Lu

Abstract: Stars can be either disrupted as tidal disruption events (TDEs) or swallowed as a whole by massive black holes (MBHs) at galactic centers when they approach sufficiently close to these MBHs. In this work, we investigate the correlations of such stellar consumption rates with both the MBH mass $M_{\rm bh}$ and the inner slope of the host galaxy mass density distribution $α$. We introduce a simplifi… ▽ More Stars can be either disrupted as tidal disruption events (TDEs) or swallowed as a whole by massive black holes (MBHs) at galactic centers when they approach sufficiently close to these MBHs. In this work, we investigate the correlations of such stellar consumption rates with both the MBH mass $M_{\rm bh}$ and the inner slope of the host galaxy mass density distribution $α$. We introduce a simplified analytical power-law model with a power-law stellar mass density distribution surrounding MBHs and separate the contributions of two-body relaxation and stellar orbital precession for the stellar orbital angular momentum evolution in nonspherical galaxy potentials. The stellar consumption rates derived from this simplified model can be well consistent with the numerical results obtained with a more realistic treatment of stellar distributions and dynamics around MBHs, providing an efficient way to estimate TDE rates. The origin of the correlations of stellar consumption rates with $M_{\rm bh}$ and $α$ are explained by the dependence of this analytical model on those MBH/host galaxy properties and by the separation of the stellar angular momentum evolution mechanisms. We propose that the strong positive correlation between the rates of stellar consumption due to two-body relaxation and $α$ provides one interpretation for the overrepresentation of TDEs found in some rare E+A/poststarburst galaxies. We find high TDE rates for giant stars, up to those for solar-type stars. The understanding of the origin of the correlations of the stellar consumption rates will be necessary for obtaining the demographics of MBHs and their host galaxies via TDEs. △ Less

Submitted 19 June, 2023; originally announced June 2023.

Comments: 17 pages, 6 figures, ApJ in press

HI Absorption in Low-power Radio AGNs Detected by FAST

Authors: Qingzheng Yu, Taotao Fang, Junfeng Wang, Jianfeng Wu

Abstract: We report the discovery of three HI absorbers toward low-power radio active galactic nuclei (AGNs) in a pilot HI absorption survey with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Compared to past studies, FAST observations have explored lower radio powers by $\sim$0.4 dex and detected these weakest absorbers at given redshifts. By comparing the gas properties and kinematics… ▽ More We report the discovery of three HI absorbers toward low-power radio active galactic nuclei (AGNs) in a pilot HI absorption survey with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Compared to past studies, FAST observations have explored lower radio powers by $\sim$0.4 dex and detected these weakest absorbers at given redshifts. By comparing the gas properties and kinematics of sources along radio powers, we aim to explore the interplay between AGN and the surrounding interstellar medium (ISM). Compared to brighter sources at similar redshifts, our observations suggest a slightly lower detection rate of HI absorption lines ($\sim$$11.5\%$) in low-power radio AGNs with $\text{log}(P_{\text{1.4 GHz}}/\text{W Hz}^{-1})=21.8-23.7$. The low-power sources with $\text{log}(P_{\text{1.4 GHz}}/\text{W Hz}^{-1})<23$ have a lower detection rate of $\sim$$6.7\%$. Due to the incompleteness of the sample, these detection rates may represent the lower limits. The selection of more extended sources and dilution by HI emission at lower redshifts may contribute to the lower detection rate of HI absorption lines. These detected absorbers present relatively narrow line widths and comparable column densities consistent with previous observations. One absorber has a symmetric profile with a large velocity offset, while the other two show asymmetric profiles that can be decomposed into multiple components, suggesting various possibilities of gas origins and kinematics. These HI absorbers may have connections with rotating disks, gas outflows, galactic gas clouds, gas fueling of the AGN, and jet-ISM interactions, which will be further investigated with the upcoming systematic survey and spatially resolved observations. △ Less

Submitted 3 August, 2023; v1 submitted 1 June, 2023; originally announced June 2023.

Comments: 12 pages, 5 figures, accepted for publication in ApJ

A dynamically discovered and characterized non-accreting neutron star -- M dwarf binary candidate

Authors: Tuan Yi, Wei-Min Gu, Zhi-Xiang Zhang, Ling-Lin Zheng, Mouyuan Sun, Junfeng Wang, Zhongrui Bai, Pei Wang, Jianfeng Wu, Yu Bai, Song Wang, Haotong Zhang, Yize Dong, Yong Shao, Xiang-Dong Li, Jia Zhang, Yang Huang, Fan Yang, Qingzheng Yu, Hui-Jun Mu, Jin-Bo Fu, Senyu Qi, Jing Guo, Xuan Fang, Chuanjie Zheng , et al. (4 additional authors not shown)

Abstract: Optical time-domain surveys can unveil and characterize exciting but less-explored non-accreting and/or non-beaming neutron stars (NS) in binaries. Here we report the discovery of such a NS candidate using the LAMOST spectroscopic survey. The candidate, designated LAMOST J112306.9+400736 (hereafter J1123), is in a single-lined spectroscopic binary containing an optically visible M star. The star's… ▽ More Optical time-domain surveys can unveil and characterize exciting but less-explored non-accreting and/or non-beaming neutron stars (NS) in binaries. Here we report the discovery of such a NS candidate using the LAMOST spectroscopic survey. The candidate, designated LAMOST J112306.9+400736 (hereafter J1123), is in a single-lined spectroscopic binary containing an optically visible M star. The star's large radial velocity variation and ellipsoidal variations indicate a relatively massive unseen companion. Utilizing follow-up spectroscopy from the Palomar 200-inch telescope and high-precision photometry from TESS, we measure a companion mass of $1.24_{-0.03}^{+0.03}~M_{\odot}$. Main-sequence stars with this mass are ruled out, leaving a NS or a massive white dwarf (WD). Although a massive WD cannot be ruled out, the lack of UV excess radiation from the companion supports the NS hypothesis. Deep radio observations with FAST yielded no detections of either pulsed or persistent emission. J1123 is not detected in numerous X-ray and gamma-ray surveys. These non-detections suggest that the NS candidate is not presently accreting and pulsing. Our work exemplifies the capability of discovering compact objects in non-accreting close binaries by synergizing the optical time-domain spectroscopy and high-cadence photometry. △ Less

Submitted 25 September, 2022; originally announced September 2022.

Comments: 53 pages, 15 figures, publication in Nature Astronomy

On Detecting Nearby Nano-Hertz Gravitational Wave Sources via Pulsar Timing Arrays

Authors: Xiao Guo, Youjun Lu, Qingjuan Yu

Abstract: Massive binary black holes (MBBHs) in nearby galactic centers, if any, may be nano-Hertz gravitational wave (GW) sources for pulsar timing arrays (PTAs) to detect. Normally the objective GWs for PTA experiments are approximated as plane waves because its sources are presumably located faraway. For nearby GW sources, however, this approximation may be inaccurate due to the curved GW wave front and… ▽ More Massive binary black holes (MBBHs) in nearby galactic centers, if any, may be nano-Hertz gravitational wave (GW) sources for pulsar timing arrays (PTAs) to detect. Normally the objective GWs for PTA experiments are approximated as plane waves because its sources are presumably located faraway. For nearby GW sources, however, this approximation may be inaccurate due to the curved GW wave front and the GW strength changes along the paths of PTA pulsar pulses. In this paper, we analyze the near-field effect in the PTA detection of nearby sources and find it is important if the source distance is less than a few tens Mpc, and ignoring this effect may lead to a significant signal-to-noise underestimation especially when the source distance is comparable to the pulsar distances. As examples, we assume a nano-Hertz MBBH source located at either the Galactic Center (GC) or the Large Magellanic Cloud (LMC) according to the observational constraints/hints on the MBBH parameter space, and estimate its detectability by current/future PTAs. We find that the GC MBBH may be detectable by the Square Kilometer Array (SKA) PTA. It is challenging for detecting the LMC MBBH; however, if a number ($N\gtrsim10$) of stable millisecond pulsars can be found in the LMC center, the MBBH may be detectable via a PTA formed by these pulsars. We further illustrate the near-field effects on the PTA detection of an isotropic GW background contributed mainly by nearby GW sources, and the resulting angular correlation is similar to the Hellings-Downs curve. △ Less

Submitted 13 October, 2022; v1 submitted 12 September, 2022; originally announced September 2022.

Comments: 24 pages, 11 figures, accepted for publication in The Astrophysical Journal

Searching for photon-ALPs mixing effects in AGN gamma-ray energy spectra

Authors: Qixin Yu, Dieter Horns

Abstract: High energy gamma-rays propagating in external magnetic fields may convert into axion-like particles (ALPs). In this case, the observed gamma-ray spectra are modified by the resulting energy-dependent conversion probability. In this study, we use the energy spectra of 20 extra-galactic gamma-ray sources recorded during 10 years of \textit{Fermi}-LAT observations. We define a test statistics based… ▽ More High energy gamma-rays propagating in external magnetic fields may convert into axion-like particles (ALPs). In this case, the observed gamma-ray spectra are modified by the resulting energy-dependent conversion probability. In this study, we use the energy spectra of 20 extra-galactic gamma-ray sources recorded during 10 years of \textit{Fermi}-LAT observations. We define a test statistics based upon the likelihood ratio to test the hypothesis for a spectral model without vs. a model with photon-ALPs coupling. The conversion probability is calculated for fixed values of the mass and two-photon coupling of the pseudo-scalar particle while the external magnetic field is characterized by the additional free parameters length scale $s$ and average field strength $B$. As a consistency check and in order to extend the analysis to include very high energy gamma-ray data, another test statistics is defined with the $χ^2$ method. We find for 18 of the 20 sources a favorable fit, particularly for Markarian~421 and NGC~1275 a significant improvement, with the hypothesis of photon-ALPs coupling in likelihood analysis. The test statistics of the sources are combined and the significance has been estimated $5.3~σ$ (test statistics summed in local maxima of all sources) and $6.0~σ$ (global maxima). The significance is estimated from dedicated simulations under the null hypotheses. The locally best-fitting values of $B$ and $s$ fall into the range that is expected for large scale magnetic fields present in relevant astrophysical environments. △ Less

Submitted 15 March, 2023; v1 submitted 29 July, 2022; originally announced August 2022.

Comments: 40pages, 45figures

Journal ref: JCAP 05 (2023) 029

On the HI Content of MaNGA Major Merger Pairs

Authors: Qingzheng Yu, Taotao Fang, Shuai Feng, Bo Zhang, C. Kevin Xu, Yunting Wang, Lei Hao

Abstract: The role of HI content in galaxy interactions is still under debate. To study the HI content of galaxy pairs at different merging stages, we compile a sample of 66 major-merger galaxy pairs and 433 control galaxies from the SDSS-IV MaNGA IFU survey. In this study, we adopt kinematic asymmetry as a new effective indicator to describe the merging stage of galaxy pairs. With archival data from the HI… ▽ More The role of HI content in galaxy interactions is still under debate. To study the HI content of galaxy pairs at different merging stages, we compile a sample of 66 major-merger galaxy pairs and 433 control galaxies from the SDSS-IV MaNGA IFU survey. In this study, we adopt kinematic asymmetry as a new effective indicator to describe the merging stage of galaxy pairs. With archival data from the HI-MaNGA survey and new observations from the Five-hundred-meter Aperture Spherical Radio Telescope (FAST), we investigate the differences in HI gas fraction ($f_{\text{HI}}$), star formation rate (SFR), and HI star formation efficiency ($\rm SFE_{\text{HI}}$) between the pair and control samples. Our results suggest that the HI gas fraction of major-merger pairs on average is marginally decreased by $\sim 15\%$ relative to isolated galaxies, implying mild HI depletion during galaxy interactions. Compared to isolated galaxies, pre-passage paired galaxies have similar $f_{\text{HI}}$, SFR and $\rm SFE_{\text{HI}}$, while pairs during pericentric passage have weakly decreased $f_{\text{HI}}$ ($-0.10\pm0.05$ dex), significantly enhanced SFR ($0.42\pm0.11$ dex) and $\rm SFE_{\text{HI}}$ ($0.48\pm0.12$ dex). When approaching the apocenter, paired galaxies show marginally decreased $f_{\text{HI}}$ ($-0.05\pm0.04$ dex), comparable SFR ($0.04\pm0.06$ dex) and $\rm SFE_{\text{HI}}$ ($0.08\pm0.08$ dex). We propose the marginally detected HI depletion may originate from the gas consumption in fuelling the enhanced $\rm H_2$ reservoir of galaxy pairs. In addition, new FAST observations also reveal an HI absorber ($N_{\text{HI}}\sim 4.7 \times 10^{21} \text{ cm}^{-2}$), which may suggest gas infalling and the triggering of AGN activity. △ Less

Submitted 13 June, 2022; originally announced June 2022.

Comments: 17 pages, 5 figures, accepted for publication in ApJ

Low Noise Frequency Domain Multiplexing of TES Bolometers using Sub-kelvin SQUIDs

Authors: Tucker Elleflot, Aritoki Suzuki, Kam Arnold, Chris Bebek, Robin H. Cantor, Kevin T. Crowley, John Groh, Tijmen de Haan, Amber Hornsby, John Joseph, Adrian T. Lee, Tiffany Liu, Joshua Montgomery, Megan Russell, Qingyang Yu

Abstract: Digital Frequency-Domain Multiplexing (DfMux) is a technique that uses MHz superconducting resonators and Superconducting Quantum Interference Device (SQUID) arrays to read out sets of Transition Edge Sensors. DfMux has been used by several Cosmic Microwave Background experiments, including most recently POLARBEAR-2 and SPT-3G with multiplexing factors as high as 68, and is the baseline readout te… ▽ More Digital Frequency-Domain Multiplexing (DfMux) is a technique that uses MHz superconducting resonators and Superconducting Quantum Interference Device (SQUID) arrays to read out sets of Transition Edge Sensors. DfMux has been used by several Cosmic Microwave Background experiments, including most recently POLARBEAR-2 and SPT-3G with multiplexing factors as high as 68, and is the baseline readout technology for the planned satellite mission LiteBIRD. Here, we present recent work focused on improving DfMux readout noise, reducing parasitic impedance, and improving sensor operation. We have achieved a substantial reduction in stray impedance by integrating the sensors, resonators, and SQUID array onto a single carrier board operated at 250 mK. This also drastically simplifies the packaging of the cryogenic components and leads to better-controlled crosstalk. We demonstrate a low readout noise level of 8.6 pA/Hz$^{-1/2}$, which was made possible by operating the SQUID array at a reduced temperature and with a low dynamic impedance. This is a factor of two improvement compared to the achieved readout noise level in currently operating Cosmic Microwave Background experiments using DfMux and represents a critical step toward maturation of the technology for the next generation of instruments. △ Less

Submitted 4 December, 2021; originally announced December 2021.

Investigating the existence of gravitomagnetic monopole in M87*

Authors: M. Ghasemi-Nodehi, Chandrachur Chakraborty, Qingjuan Yu, Youjun Lu

Abstract: We examine the possibility for the existence of gravitomagnetic monopole ($n_*$) in M87* by using the results obtained from its first Event Horizon Telescope image. By numerically deducing the shadow sizes in Kerr-Taub-NUT (KTN) spacetime, we show that the shadow size increases with increasing $|n_*|$ for a fixed Kerr parameter $|a_*|$ in case of the KTN black hole, whereas for a KTN naked singula… ▽ More We examine the possibility for the existence of gravitomagnetic monopole ($n_*$) in M87* by using the results obtained from its first Event Horizon Telescope image. By numerically deducing the shadow sizes in Kerr-Taub-NUT (KTN) spacetime, we show that the shadow size increases with increasing $|n_*|$ for a fixed Kerr parameter $|a_*|$ in case of the KTN black hole, whereas for a KTN naked singularity it increases with increasing $n_*$ for a fixed $a_* > 0$ if $n_* > -\cot 17^{\circ}$ . In general, the asymmetry of shadow shape increases if the central dark object in M87 is a KTN/Kerr naked singularity instead of a KTN/Kerr black hole. We find that a non-zero gravitomagnetic monopole is still compatible with the current EHT observations, in which case the upper limit of $n_*$ cannot be greater than $1.1$, i.e., $n_* \lesssim 1.1$ for the prograde rotation ($a_* > 0$), and the lower limit of $n_*$ cannot be less than $-1.1$, i.e., $ n_* \gtrsim -1.1$ for the retrograde rotation ($a_* < 0$). Moreover, if the circularity of the shadow can be measured on a precision of $\lesssim 1\%$, the Kerr and KTN naked singularities can be falsified for M87*. △ Less

Submitted 27 October, 2021; v1 submitted 30 September, 2021; originally announced September 2021.

Comments: 14 pages, 12 figures, accepted for publication in EPJC

Journal ref: Eur. Phys. J. C 81:939 (2021)

Gravitational Wave From Axion-like Particle Inflation

Authors: Wei Cheng, Tao Qian, Qing Yu, Hua Zhou, Rui-Yu Zhou

Abstract: In this paper, we investigate the Axion-like Particle inflation by applying the multi-nature inflation model, where the end of inflation is achieved through the phase transition (PT). The events of PT should not be less than $200$, which results in the free parameter $n\geq404$. Under the latest CMB restrictions, we found that the inflation energy is fixed at $10^{15} \rm{GeV}$. Then, we deeply di… ▽ More In this paper, we investigate the Axion-like Particle inflation by applying the multi-nature inflation model, where the end of inflation is achieved through the phase transition (PT). The events of PT should not be less than $200$, which results in the free parameter $n\geq404$. Under the latest CMB restrictions, we found that the inflation energy is fixed at $10^{15} \rm{GeV}$. Then, we deeply discussed the corresponding stochastic background of the primordial gravitational wave (GW) during inflation. We study the two kinds of $n$ cases, i.e., $n=404, 2000$. We observe that the magnitude of $n$ is negligible for the physical observations, such as $n_s$, $r$, $Λ$, and $Ω_{\rm{GW}}h^2$. In the low-frequency regions, the GW is dominated by the quantum fluctuations, and this GW can be detected by Decigo at $10^{-1}~\rm{Hz}$. However, GW generated by PT dominates the high-frequency regions, which is expected to be detected by future 3DSR detector. △ Less

Submitted 9 July, 2021; originally announced July 2021.

Comments: 7 pages, 5 figures

Seeking celestial Positronium with an OH-suppressed diffraction-limited spectrograph

Authors: J. Gordon Robertson, Simon Ellis, Qingshan Yu, Joss Bland-Hawthorn, Christopher Betters, Martin Roth, Sergio Leon-Saval

Abstract: Celestially, Positronium (Ps), has only been observed through gamma-ray emission produced by its annihilation. However, in its triplet state, a Ps atom has a mean lifetime long enough for electronic transitions to occur between quantum states. This produces a recombination spectrum observable in principle at near IR wavelengths, where angular resolution greatly exceeding that of the gamma-ray obse… ▽ More Celestially, Positronium (Ps), has only been observed through gamma-ray emission produced by its annihilation. However, in its triplet state, a Ps atom has a mean lifetime long enough for electronic transitions to occur between quantum states. This produces a recombination spectrum observable in principle at near IR wavelengths, where angular resolution greatly exceeding that of the gamma-ray observations is possible. However, the background in the NIR is dominated by extremely bright atmospheric hydroxyl (OH) emission lines. In this paper we present the design of a diffraction-limited spectroscopic system using novel photonic components - a photonic lantern, OH Fiber Bragg Grating filters, and a photonic TIGER 2-dimensional pseudo-slit - to observe the Ps Balmer alpha line at 1.3122 microns for the first time. △ Less

Submitted 18 June, 2021; originally announced June 2021.

Comments: 6 pages, 9 Figures, 2 Tables. Accepted to Applied Optics feature issue on Astrophotonics

Journal ref: Applied Optics Vol 60, No. 19, D122-D128, 2021

Extragalactic HI 21-cm absorption line observations with the Five-hundred-meter Aperture Spherical radio Telescope

Authors: Bo Zhang, Ming Zhu, Zhong-Zu Wu, Qing-Zheng Yu, Peng Jiang, You-Ling Yue, Meng-Lin Huang, Qiao-Li Hao

Abstract: We present a pilot study of extragalactic HI 21-cm absorption lines using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). We observed 5 continuum sources with HI absorption features firstly identified in the 40% data release of the Arecibo Legacy Fast Arecibo L-Band Feed Array (ALFA) Survey (ALFALFA), including two systems later detected by the Westerbork Synthesis Radio Telescop… ▽ More We present a pilot study of extragalactic HI 21-cm absorption lines using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). We observed 5 continuum sources with HI absorption features firstly identified in the 40% data release of the Arecibo Legacy Fast Arecibo L-Band Feed Array (ALFA) Survey (ALFALFA), including two systems later detected by the Westerbork Synthesis Radio Telescope (WSRT). Most of our observations were carried out during the FAST commissioning phase, and we have tested different observing modes, as well as data reduction methods, to produce the best spectra. Our observations successfully confirmed the existence of HI absorption lines in all these systems, including two sources that were marginally detected by ALFALFA. We fitted the HI profiles with single or double of Gaussian functions, and calculated the HI column densities of each source. The HI absorption profiles obtained by FAST show much higher spectral resolution and higher S/N ratio than the existing data in the literature, thus demonstrating the power of FAST in revealing detailed structures of HI absorption lines. Our pilot observations and tests have enabled us to develop a strategy to search for HI absorption sources using the data from the FAST extragalactic HI survey, which is one of the key projects undertaken at FAST. We expect that over 1,500 extragalactic HI absorbing systems could be detected with survey data, based on sensitivity level we achieved in pilot observations. △ Less

Submitted 11 March, 2021; originally announced March 2021.

Comments: 13 pages, 14 figures. Accepted by Monthly Notices of the Royal Astronomical Society

Black holes and the supermassive compact object at the Galactic center: multi-arts of thought and nature

Authors: Qingjuan Yu

Abstract: This is an invited commentary on the Nobel Prize in Physics 2020 which was awarded to Roger Penrose "for the discovery that black hole formation is a robust prediction of the general theory of relativity," and Reinhard Genzel and Andrea Ghez "for the discovery of a supermassive compact object at the center of our galaxy." This is an invited commentary on the Nobel Prize in Physics 2020 which was awarded to Roger Penrose "for the discovery that black hole formation is a robust prediction of the general theory of relativity," and Reinhard Genzel and Andrea Ghez "for the discovery of a supermassive compact object at the center of our galaxy." △ Less

Submitted 7 December, 2020; originally announced December 2020.

Comments: 5 pages, 1 figure. Published in The Innovation

Journal ref: 2020, The Innovation, 1, 100063

Cosmic distributions of stellar tidal disruptions by massive black holes at galactic centers

Authors: Yunfeng Chen, Qingjuan Yu, Youjun Lu

Abstract: Stars can be consumed (either tidally disrupted or swallowed whole) by massive black holes (MBHs) at galactic centers when they move into the vicinity of the MBHs. In this study, we investigate the rates of stellar consumption by central MBHs and their cosmic distributions, including the effects of triaxial galaxy shapes in enhancing the reservoir of low-angular-momentum stars and incorporating re… ▽ More Stars can be consumed (either tidally disrupted or swallowed whole) by massive black holes (MBHs) at galactic centers when they move into the vicinity of the MBHs. In this study, we investigate the rates of stellar consumption by central MBHs and their cosmic distributions, including the effects of triaxial galaxy shapes in enhancing the reservoir of low-angular-momentum stars and incorporating realistic galaxy distributions. We find that the enhancement in the stellar consumption rates due to triaxial galaxy shapes can be significant, by a factor of ~3 for MBH mass $M_{\rm BH}\sim10^5$-$10^6$Msun and up to more than one order of magnitude for $M_{\rm BH}\gtrsim10^8$Msun. Only for $M_{\rm BH}\lesssim10^7$Msun are the stellar consumption rates significantly higher in galaxies with steeper inner surface brightness profiles. The average (per galaxy) stellar consumption rates correlate with central MBH masses positively for $M_{\rm BH}\gtrsim10^7$Msun and negatively for $M_{\rm BH}\lesssim10^7$Msun. The volumetric stellar tidal disruption rates are ~$3\times10^{-5}$/yr/Mpc$^3$ for MBHs in the mass range of $10^5$-$10^8$Msun at z=0; and the volumetric stellar consumption rates by MBHs with higher masses are ~$10^{-6}$/yr/Mpc$^3$, which can be the stellar tidal disruption rate if the high-mass BHs are extremely spinning Kerr BHs or the rate of being swallowed if those BHs are Schwarzschild ones. The volumetric stellar consumption rates decrease with increasing redshift, and the decrease is relatively mild for $M_{\rm BH}\sim10^5$-$10^7$Msun and stronger for higher $M_{\rm BH}$. Most of the stellar tidal disruption events (TDEs) at z=0 occur in galaxies with mass $M_{\rm gal}\lesssim10^{11}$Msun, and about 1%-2% of the TDEs can occur in high-mass galaxies with $M_{\rm gal}\gtrsim10^{11}$Msun. △ Less

Submitted 7 December, 2020; originally announced December 2020.

Comments: 20 pages, 9 figures

Journal ref: 2020, ApJ, 900, 191

Reconciling hints on axion-like-particles from high-energy gamma rays with stellar bounds

Authors: Gautham Adamane Pallathadka, Francesca Calore, Pierluca Carenza, Maurizio Giannotti, Dieter Horns, Julian Kuhlmann, Jhilik Majumdar, Alessandro Mirizzi, Andreas Ringwald, Anton Sokolov, Franziska Stief, Qixin Yu

Abstract: It has been recently claimed by two different groups that the spectral modulation observed in gamma rays from Galactic pulsars and supernova remnants can be due to conversion of photons into ultra-light axion-like-particles (ALPs) in large-scale Galactic magnetic fields. While we show the required best-fit photon-ALP coupling, $g_{aγ} \sim 2 \times 10^{-10}$ GeV${}^{-1}$, to be consistent with con… ▽ More It has been recently claimed by two different groups that the spectral modulation observed in gamma rays from Galactic pulsars and supernova remnants can be due to conversion of photons into ultra-light axion-like-particles (ALPs) in large-scale Galactic magnetic fields. While we show the required best-fit photon-ALP coupling, $g_{aγ} \sim 2 \times 10^{-10}$ GeV${}^{-1}$, to be consistent with constraints from observations of photon-ALPs mixing in vacuum, this is in conflict with other bounds, specifically from the CAST solar axion limit, from the helium-burning lifetime in globular clusters, and from the non-observations of gamma rays in coincidence with SN 1987A. In order to reconcile these different results, we propose that environmental effects in matter would suppress the ALP production in dense astrophysical plasma, allowing to relax previous bounds and make them compatible with photon-ALP conversions in the low-density Galactic medium. If this explanation is correct, the claimed ALP signal would be on the reach of next-generations laboratory experiments such as ALPS II. △ Less

Submitted 10 November, 2021; v1 submitted 18 August, 2020; originally announced August 2020.

Comments: v2: major corrections to match the published version. 28 pages, 14 figures

Journal ref: JCAP 11 (2021) 036

Dynamical evolution of cosmic supermassive binary black holes and their gravitational wave radiation

Authors: Yunfeng Chen, Qingjuan Yu, Youjun Lu

Abstract: We investigate the evolution of supermassive binary black holes (BBHs) in galaxies with realistic property distributions and the gravitational-wave (GW) radiation from the cosmic population of these BBHs. We incorporate a comprehensive treatment of the dynamical interactions of the BBHs with their environments by including the effects of galaxy triaxial shapes and inner stellar distributions, and… ▽ More We investigate the evolution of supermassive binary black holes (BBHs) in galaxies with realistic property distributions and the gravitational-wave (GW) radiation from the cosmic population of these BBHs. We incorporate a comprehensive treatment of the dynamical interactions of the BBHs with their environments by including the effects of galaxy triaxial shapes and inner stellar distributions, and generate a large number of BBH evolution tracks. By combining these BBH evolution tracks, galaxy mass functions, galaxy merger rates, and supermassive black hole-host galaxy relations into our model, we obtain the statistical distributions of surviving BBHs, BBH coalescence rates, the strength of their GW radiation, and the stochastic GW background (GWB) contributed by the cosmic BBH population. About ~1%-3% (or ~10%) of supermassive BHs at nearby galactic centers are expected to be binaries with mass ratio >1/3 (or >1/100). The characteristic strain amplitude of the GWB at frequency 1/yr is estimated to be ~$2.0^{+1.4}_{-0.8}\times 10^{-16}$, and the upper bound of its results obtained with the different BH-host galaxy relations can be up to $5.4\times 10^{-16}$, which await testing by future experiments (e.g., the Square Kilometer Array, FAST, Next-Generation Very Large Array). The turnover frequency of the GWB spectrum is at ~0.25nHz. The uncertainties on the above estimates and prospects for detecting individual sources are also discussed. The application of the cosmic BBH population to the Laser Interferometer Space Antenna (LISA) band provides a lower limit to the detection rate of BBHs by LISA, ~0.9/yr. △ Less

Submitted 23 July, 2020; v1 submitted 21 May, 2020; originally announced May 2020.

Comments: 43 pages, 22 figures

Journal ref: 2020, ApJ, 897, 86

Cosmic muon flux measurement and tunnel overburden structure imaging

Authors: Ran Han, Qian Yu, Zhiwei Li, Jingtai Li, Yaping Cheng, Bin Liao, Lixiang Jiang, Sidao Ni, Tianfang Liu, Zheng Wang

Abstract: We present a cosmic ray muon tomographic experiment for measuring the muon flux and imaging the tunnel overburden structures in Changshu, China. The device used in this study is a tracking detector based on the plastic scintillator with SiPM technology, which can be conveniently operated in field works. The compact system with sensitive area of $6400 cm^2$ can measure the angular distribution of c… ▽ More We present a cosmic ray muon tomographic experiment for measuring the muon flux and imaging the tunnel overburden structures in Changshu, China. The device used in this study is a tracking detector based on the plastic scintillator with SiPM technology, which can be conveniently operated in field works. The compact system with sensitive area of $6400 cm^2$ can measure the angular distribution of cosmic muons. It's able to image the overburden density length from the surface of overburden to the detector along the muon tracks. The open sky muon flux measurement outside the tunnel has a good agreement with the modified Gassier Formula model. The distributions of muon flux at three positions inside the tunnel are very similar to that of open sky. Assuming the average density of overburden compact sandstone is $2.65 g/cm^3$, the overburden thickness can be obtained from the density length derived from the difference of muon flux inside and outside the tunnel. Moreover, for known penetrated lengths (i.e., topography of overburden), the density anomalies of the overburden can also been obtained. This study suggests a potential application for imaging and detecting subsurface structures in civil engineering, tunnels or caverns with the cosmic ray muon telescope. △ Less

Submitted 27 March, 2020; originally announced March 2020.

Journal ref: ournal of Instrumentation, Volume 15, June 2020

A PRESTO-based Parallel Pulsar Search Pipeline Used for FAST Drift Scan Data

Authors: Qiuyu Yu, Zhichen Pan, Lei Qian, Shen Wang, Youling Yue, Menglin Huang, Qiaoli Hao, Shanping You, Bo Peng, Yan Zhu, Lei Zhang, Zhijie Liu

Abstract: We developed a pulsar search pipeline based on PRESTO (PulsaR Exploration and Search Toolkit). This pipeline simply runs dedispersion, FFT (Fast Fourier Transformation), and acceleration search in process-level parallel to shorten the processing time. With two parallel strategies, the pipeline can highly shorten the processing time in both the normal searches or acceleration searches. This pipelin… ▽ More We developed a pulsar search pipeline based on PRESTO (PulsaR Exploration and Search Toolkit). This pipeline simply runs dedispersion, FFT (Fast Fourier Transformation), and acceleration search in process-level parallel to shorten the processing time. With two parallel strategies, the pipeline can highly shorten the processing time in both the normal searches or acceleration searches. This pipeline was first tested with PMPS (Parkes Multibeam Pulsar Survery) data and discovered two new faint pulsars. Then, it was successfully used in processing the FAST (Five-hundred-meter Aperture Spherical radio Telescope) drift scan data with tens of new pulsar discoveries up to now. The pipeline is only CPU-based and can be easily and quickly deployed in computing nodes for testing purposes or data processes. △ Less

Submitted 29 December, 2019; originally announced December 2019.

Comments: 10 pages, 2 figures, accepted by Research in Astronomy and Astrophysics (RAA)

The Origin of the Stellar Mass-Stellar Metallicity Relation in the Milky Way Satellites and beyond

Authors: Moran Xia, Qingjuan Yu

Abstract: Observations and semianalytical galaxy formation and evolution models (SAMs) have suggested the existence of a stellar mass-stellar metallicity relation (MZR), which is shown to be universal for different types of galaxies over a large range of stellar masses ($M_*\sim 10^3$-$10^{11}M_\odot$) and dark matter (DM) halo masses ($M_{\rm halo}\sim 10^9$-$10^{15}h^{-1}M_\odot$). In this work, we constr… ▽ More Observations and semianalytical galaxy formation and evolution models (SAMs) have suggested the existence of a stellar mass-stellar metallicity relation (MZR), which is shown to be universal for different types of galaxies over a large range of stellar masses ($M_*\sim 10^3$-$10^{11}M_\odot$) and dark matter (DM) halo masses ($M_{\rm halo}\sim 10^9$-$10^{15}h^{-1}M_\odot$). In this work, we construct a chemical evolution model to investigate the origin of the MZR, including both the effects of gas inflows and outflows in galaxies. We solve the MZR from the chemical evolution model, by assuming that the cold gas mass ($M_{\rm cold}$) and the stellar feedback efficiency ($β$) follow some power-law scaling relationships with $M_*$ during the growth of a galaxy, i.e., $M_{\rm cold}\propto M_*^{α_{\rm gs}}$ and $β\propto M_*^{α_{β{\rm s}}}$. We use the SAM to obtain these power-law scaling relations, which appear to be roughly universal over a large range of stellar masses for both satellites and central galaxies within a large range of halo masses. The range of the MZRs produced by our models is in a narrow space, which provides support to the universality of the MZRs. The formation of the MZR is a result caused jointly by that the cold gas fraction decreases with increasing $M_*$ and by that the stellar feedback efficiency decreases with increasing $M_*$ in the galaxy growth, and the exponent in the MZR is around $-α_{β{\rm s}}$ or $1-α_{\rm gs}$. The MZR represents an ''average'' evolutional track for the stellar metallicity of a galaxy. The comparison of our model with some previous models for the origin of MZRs is also discussed. △ Less

Submitted 30 May, 2019; originally announced May 2019.

Comments: 22 pages, 9 figures, accepted for publication in the Astrophysical Journal

Is the stellar mass-stellar metallicity relation universal in the Milky Way satellites and beyond?

Authors: Moran Xia, Qingjuan Yu

Abstract: Observations reveal a universal stellar mass-stellar metallicity relation (MZR) existing in Local Group dwarfs of different types, $Z_*\propto M_*^α$ with $α=0.30\pm0.02$. In this work, we investigate the "universality" of the MZRs for both satellites and central galaxies in a large number of different host dark matter halos covering a large mass range of $10^9$-$10^{15}h^{-1}M_\odot$, by using a… ▽ More Observations reveal a universal stellar mass-stellar metallicity relation (MZR) existing in Local Group dwarfs of different types, $Z_*\propto M_*^α$ with $α=0.30\pm0.02$. In this work, we investigate the "universality" of the MZRs for both satellites and central galaxies in a large number of different host dark matter halos covering a large mass range of $10^9$-$10^{15}h^{-1}M_\odot$, by using a semianalytical galaxy formation and evolution model. We obtain the following results. (1) The exponents ($α$) for the MZRs of the satellites in halos with the same mass as the Milky Way halo but different individual assembly histories are mostly $\sim$0.2-0.4, i.e., having a scatter of $\sim 0.2$; and the scatter of $α$ increases with decreasing halo masses. (2) The MZR relations are changed little by the variation of halo masses and the classification between central galaxies and satellites, if many halos with the same mass are stacked together. (3) A double power law exists in the MZR relations for both central galaxies and stacked satellites, with $α\sim$0.2-0.4 at $10^3M_\odot< M_*<10^{8} M_\odot$ and a relatively higher $α\sim0.5$ at $10^8M_\odot<M_*<10^{11}M_\odot$. (4) The high-mass satellites ($M_*>10^8M_\odot$) existing mostly in high-mass halos can lead to an apparent increase of $α$ (from $\sim0.2$ to $\sim0.4$) with increasing host halo masses shown in the single power law fitting results of stacked satellites. The universality of the MZR suggests the common physical processes in stellar formation and chemical evolution of galaxies can be unified over a large range of galaxy masses and halo masses. △ Less

Submitted 14 March, 2019; originally announced March 2019.

Comments: 13 pages, 4 figures, accepted for publication in the Astrophysical Journal

Simulating the galaxy cluster "El Gordo": gas motion, kinetic Sunyaev-Zel'dovich signal, and X-ray line features

Authors: Congyao Zhang, Qingjuan Yu, Youjun Lu

Abstract: The massive galaxy cluster "El Gordo" (ACT-CL J0102--4915) is a rare merging system with a high collision speed suggested by multi-wavelength observations and the theoretical modeling. Zhang et al. (2015) propose two types of mergers, a nearly head-on merger and an off-axis merger with a large impact parameter, to reproduce most of the observational features of the cluster, by using numerical simu… ▽ More The massive galaxy cluster "El Gordo" (ACT-CL J0102--4915) is a rare merging system with a high collision speed suggested by multi-wavelength observations and the theoretical modeling. Zhang et al. (2015) propose two types of mergers, a nearly head-on merger and an off-axis merger with a large impact parameter, to reproduce most of the observational features of the cluster, by using numerical simulations. The different merger configurations of the two models result in different gas motion in the simulated clusters. In this paper, we predict the kinetic Sunyaev-Zel'dovich (kSZ) effect, the relativistic correction of the thermal Sunyaev-Zel'dovich (tSZ) effect, and the X-ray spectrum of this cluster, based on the two proposed models. We find that (1) the amplitudes of the kSZ effect resulting from the two models are both on the order of $ΔT/T\sim10^{-5}$; but their morphologies are different, which trace the different line-of-sight velocity distributions of the systems; (2) the relativistic correction of the tSZ effect around $240 {\rm\,GHz}$ can be possibly used to constrain the temperature of the hot electrons heated by the shocks; and (3) the shift between the X-ray spectral lines emitted from different regions of the cluster can be significantly different in the two models. The shift and the line broadening can be up to $\sim 25{\rm\,eV}$ and $50{\rm\,eV}$, respectively. We expect that future observations of the kSZ effect and the X-ray spectral lines (e.g., by ALMA, XARM) will provide a strong constraint on the gas motion and the merger configuration of ACT-CL J0102--4915. △ Less

Submitted 22 November, 2017; originally announced November 2017.

Comments: 10 pages, 6 figures, submitted to ApJ

Did the 2000 November 8 solar flare accelerate protons to >=40 GeV?

Authors: R. G. Wang, L. K. Ding, Y. Q. Ma, X. H. Ma, Q. Q. Zhu, C. G. Yang, H. H. Kuang, Z. Q. Yu, Z. G. Yao, Y. P. Xu

Abstract: It has been reported that a 5.7sigma directional muon excess coincident with the 2000 July 14 solar flare was registered by the L3 precision muon spectrometer [Ruiguang Wang, Astroparticle Phys., 31(2009) 149]. Using a same analysis method and similar criteria of event selection, we have analyzed the L3 precision muon spectrometer data during November 2000. The result shows that a 4.7sigma muon ex… ▽ More It has been reported that a 5.7sigma directional muon excess coincident with the 2000 July 14 solar flare was registered by the L3 precision muon spectrometer [Ruiguang Wang, Astroparticle Phys., 31(2009) 149]. Using a same analysis method and similar criteria of event selection, we have analyzed the L3 precision muon spectrometer data during November 2000. The result shows that a 4.7sigma muon excess appeared at a time coincident with the solar flare of 8 of November 2000. This muon excess corresponds to above 40 GeV primary protons which came from a sky cell of solid angle 0.048 sr. The probability of being a background fluctuation is estimated to be about 0.1%. It has been convinced that solar protons could be accelerated to tens of GeV in those Class X solar flares which usually arose solar cosmic ray ground level enhancement (GLE) events. However, whether a Class M solar flare like the non-GLE event of 8 November 2000 may also accelerate solar protons to such high energies? It is interesting and noteworthy. △ Less

Submitted 8 August, 2017; originally announced August 2017.

X-ray Eclipses of Active Galactic Nuclei

Authors: Fupeng Zhang, Qingjuan Yu, Youjun Lu

Abstract: X-ray variation is a ubiquitous feature of active galactic nuclei (AGNs), however, its origin is not well understood. In this paper, we show that the X-ray flux variations in some AGNs, and correspondingly the power spectral densities (PSDs) of the variations, may be interpreted as being caused by absorptions of eclipsing clouds or clumps in the broad line region (BLR) and the dusty torus. By perf… ▽ More X-ray variation is a ubiquitous feature of active galactic nuclei (AGNs), however, its origin is not well understood. In this paper, we show that the X-ray flux variations in some AGNs, and correspondingly the power spectral densities (PSDs) of the variations, may be interpreted as being caused by absorptions of eclipsing clouds or clumps in the broad line region (BLR) and the dusty torus. By performing Monte-Carlo simulations for a number of plausible cloud models, we systematically investigate the statistics of the X-ray variations resulting from the cloud eclipsing and the PSDs of the variations. For these models, we show that the number of eclipsing events can be significant and the absorption column densities due to those eclipsing clouds can be in the range from 10^{21} to 10^{24} cm^{-2}, leading to significant X-ray variations. We find that the PSDs obtained from the mock observations for the X-ray flux and the absorption column density resulting from these models can be described by a broken double power law, similar to those directly measured from observations of some AGNs. The shape of the PSDs depend strongly on the kinematic structures and the intrinsic properties of the clouds in AGNs. We demonstrate that the X-ray eclipsing model can naturally lead to a strong correlation between the break frequencies (and correspondingly the break timescales) of the PSDs and the masses of the massive black holes (MBHs) in the model AGNs, which can be well consistent with the one obtained from observations. Future studies of the PSDs of the AGN X-ray (and possibly also the optical-UV) flux and column density variations may provide a powerful tool to constrain the structure of the BLR and the torus and to estimate the MBH masses in AGNs. △ Less

Submitted 14 August, 2017; v1 submitted 20 July, 2017; originally announced July 2017.

Comments: 25 pages, 10 figures

Journal ref: 2017, ApJ, 845, 88

Prospects for constraining the spin of the massive black hole at the Galactic center via the relativistic motion of a surrounding star

Authors: Qingjuan Yu, Fupeng Zhang, Youjun Lu

Abstract: In this paper, we investigate the quality of constraining the spin of the massive black hole (MBH) at the Galactic center (GC) by using full general relativistic simulations of the motion of a surrounding star. We obtain the dependence mapping of the spin-induced signals on any spin direction of the MBH for given example stars, which indicates the feasibility to test whether the spin direction is… ▽ More In this paper, we investigate the quality of constraining the spin of the massive black hole (MBH) at the Galactic center (GC) by using full general relativistic simulations of the motion of a surrounding star. We obtain the dependence mapping of the spin-induced signals on any spin direction of the MBH for given example stars, which indicates the feasibility to test whether the spin direction is the same as the normal of the young stellar disk located at the GC, and, further to provide insights into the assembly history of the MBH. We demonstrate the quality of constraining the MBH spin that may be achieved, given any set of the astrometric and the redshift precisions of observational facilities. We find that in the ranges of the astrometric and the velocity precisions with 1--30$μ$as and 0.1--10 km/s, an improvement in astrometric precision would be more effective at improving the quality of constraining the spin than an improvement in velocity precision. We obtain the parameter space of the semimajor axis and the eccentricity for the orbit of the target star that a high-precision constraint on the GC MBH spin can be obtained via the motion of the star. Our results show that the spin of the GC MBH can be constrained with 1-sigma error <~0.1 or even >~0.02 by monitoring the orbital motion of a star, if existing as expected, with semimajor axis <~300AU and eccentricity >~0.95 over a period shorter than a decade through future facilities. △ Less

Submitted 6 September, 2016; v1 submitted 24 June, 2016; originally announced June 2016.

Comments: 10 pages, 7 figures, updated to the journal version

Journal ref: 2016, ApJ, 827, 114

A Baryonic Effect on the Merger Timescale of Galaxy Clusters

Authors: Congyao Zhang, Qingjuan Yu, Youjun Lu

Abstract: Accurate estimation of the merger timescale of galaxy clusters is important to understand the cluster merger process and further the formation and evolution of the large-scale structure of the universe. In this paper, we explore a baryonic effect on the merger timescale of galaxy clusters by using hydrodynamical simulations. We find that the baryons play an important role in accelerating the merge… ▽ More Accurate estimation of the merger timescale of galaxy clusters is important to understand the cluster merger process and further the formation and evolution of the large-scale structure of the universe. In this paper, we explore a baryonic effect on the merger timescale of galaxy clusters by using hydrodynamical simulations. We find that the baryons play an important role in accelerating the merger process. The merger timescale decreases with increasing the gas fraction of galaxy clusters. For example, the merger timescale is shortened by a factor of up to 3 for merging clusters with gas fractions 0.15, compared with the timescale obtained with zero gas fractions. The baryonic effect is significant for a wide range of merger parameters and especially more significant for nearly head-on mergers and high merging velocities. The baryonic effect on the merger timescale of galaxy clusters is expected to have impacts on the structure formation in the universe, such as the cluster mass function and massive substructures in galaxy clusters, and a bias of "no-gas" may exist in the results obtained from the dark matter-only cosmological simulations. △ Less

Submitted 8 February, 2016; originally announced February 2016.

Comments: 9 pages, 6 figures, accepted for publication in the Astrophysical Journal

Simulating the galaxy cluster "El Gordo" and identifying the merger configuration

Authors: Congyao Zhang, Qingjuan Yu, Youjun Lu

Abstract: The observational features of the massive galaxy cluster "El Gordo" (ACT-CL J0102-4915), such as the X-ray emission, the Sunyaev-Zel'dovich (SZ) effect, and the surface mass density distribution, indicate that they are caused by an exceptional ongoing high-speed collision of two galaxy clusters, similar to the well-known Bullet Cluster. We perform a series of hydrodynamical simulations to investig… ▽ More The observational features of the massive galaxy cluster "El Gordo" (ACT-CL J0102-4915), such as the X-ray emission, the Sunyaev-Zel'dovich (SZ) effect, and the surface mass density distribution, indicate that they are caused by an exceptional ongoing high-speed collision of two galaxy clusters, similar to the well-known Bullet Cluster. We perform a series of hydrodynamical simulations to investigate the merging scenario and identify the initial conditions for the collision in ACT-CL J0102-4915. By surveying the parameter space of the various physical quantities that describe the two colliding clusters, including their total mass (M), mass ratio (ξ), gas fractions (f_b), initial relative velocity (V), and impact parameter (P), we find out an off-axis merger with P~800h_{70}^{-1}kpc, V~2500km/s, M~3x10^{15}Msun, and ξ=3.6 that can lead to most of the main observational features of ACT-CL J0102-4915. Those features include the morphology of the X-ray emission with a remarkable wake-like substructure trailing after the secondary cluster, the X-ray luminosity and the temperature distributions, and also the SZ temperature decrement. The initial relative velocity required for the merger is extremely high and rare compared to that inferred from currently available Lambda cold dark matter (LCDM) cosmological simulations, which raises a potential challenge to the LCDM model, in addition to the case of the Bullet Cluster. △ Less

Submitted 9 November, 2015; originally announced November 2015.

Comments: 30 pages, 12 figures

Journal ref: 2015, ApJ, 813, 129

On Testing the Kerr Metric of the Massive Black Hole in the Galactic Center via Stellar Orbital Motion: Full General Relativistic Treatment

Authors: Fupeng Zhang, Youjun Lu, Qingjuan Yu

Abstract: The S-stars in the Galactic center (GC) are anticipated to provide unique dynamical constraint on the spin of the GC massive black hole (MBH). In this paper, we develop a fast full general relativistic method to simultaneously constrain the MBH mass, spin, and spin direction by considering both the motion of a star and the propagation of photons from the star to a distant observer. Assuming some e… ▽ More The S-stars in the Galactic center (GC) are anticipated to provide unique dynamical constraint on the spin of the GC massive black hole (MBH). In this paper, we develop a fast full general relativistic method to simultaneously constrain the MBH mass, spin, and spin direction by considering both the motion of a star and the propagation of photons from the star to a distant observer. Assuming some example stars, we demonstrate that the spin-induced effects on the projected trajectory and redshift curve of a star depend on both the value and the direction of the spin. The maximum effects over a full orbit can differ by a factor upto more than one order of magnitude for cases with significantly different spin directions. Adopting the Markov Chain Monte Carlo fitting technique, we illustrate that the spin of the GC MBH is likely to be well constrained by using the motion of S0-2/S2 over a period of ~45yr if it is close to one and the astrometric and spectroscopic precisions (sigma_p,sigma_Z) can be as high as (10muas, 1km/s). In the mean time, the distance from the sun to the GC and the MBH mass can also be constrained to an unprecedented accuracy (0.01%-0.1%). If there exists a star with semimajor axis significantly smaller than that of S0-2/S2 and eccentricity larger than that of S0-2/S2, the MBH spin can be constrained with high accuracy over a period of <~10yr for (sigma_p,sigma_Z) ~ (10muas,1km/s), even if the spin is only moderately large (>~0.2). △ Less

Submitted 25 August, 2015; originally announced August 2015.

Comments: 29 pages, 17 figures

Journal ref: The Astrophysical Journal, 2015, 809, 127

A probable Milli-Parsec Supermassive Binary Black Hole in the Nearest Quasar Mrk 231

Authors: Chang-Shuo Yan, Youjun Lu, Xinyu Dai, Qingjuan Yu

Abstract: Supermassive binary black holes (BBHs) are unavoidable products of galaxy mergers and are expected to exist in the cores of many quasars. Great effort has been made during the past several decades to search for BBHs among quasars; however, observational evidence for BBHs remains elusive and ambiguous, which is difficult to reconcile with theoretical expectations. In this paper, we show that the di… ▽ More Supermassive binary black holes (BBHs) are unavoidable products of galaxy mergers and are expected to exist in the cores of many quasars. Great effort has been made during the past several decades to search for BBHs among quasars; however, observational evidence for BBHs remains elusive and ambiguous, which is difficult to reconcile with theoretical expectations. In this paper, we show that the distinct optical-to-UV spectrum of Mrk 231 can be well interpreted as emission from accretion flows onto a BBH, with a semimajor axis of ~590AU and an orbital period of ~1.2 year. The flat optical and UV continua are mainly emitted from a circumbinary disk and a mini-disk around the secondary black hole (BH), respectively; and the observed sharp drop off and flux deficit at wavelength lambda ~ 4000-2500 Angstrom is due to a gap (or hole) opened by the secondary BH migrating within the circumbinary disk. If confirmed by future observations, this BBH will provide a unique laboratory to study the interplay between BBHs and accretion flows onto them. Our result also demonstrates a new method to find sub-parsec scale BBHs by searching for deficits in the optical-to-UV continuum among the spectra of quasars. △ Less

Submitted 8 September, 2015; v1 submitted 25 August, 2015; originally announced August 2015.

Comments: a typo in equation (2) and also in equation (5) of the Appendix is fixed; 9 pages, 7 figures

Journal ref: The Astrophysical Journal, 2015, 809, 117

High-Resolution Spectroscopic Study of Extremely Metal-Poor Star Candidates from the SkyMapper Survey

Authors: Heather. R. Jacobson, Stefan Keller, Anna Frebel, Andrew R. Casey, Martin Asplund, Michael S. Bessell, Gary S. Da Costa, Karin Lind, Anna F. Marino, John E. Norris, Jose M. Pena, Brian P. Schmidt, Patrick Tisserand, Jennifer M. Walsh, David Yong, Qinsi Yu

Abstract: The SkyMapper Southern Sky Survey is carrying out a search for the most metal-poor stars in the Galaxy. It identifies candidates by way of its unique filter set that allows for estimation of stellar atmospheric parameters. The set includes a narrow filter centered on the Ca II K 3933A line, enabling a robust estimate of stellar metallicity. Promising candidates are then confirmed with spectroscopy… ▽ More The SkyMapper Southern Sky Survey is carrying out a search for the most metal-poor stars in the Galaxy. It identifies candidates by way of its unique filter set that allows for estimation of stellar atmospheric parameters. The set includes a narrow filter centered on the Ca II K 3933A line, enabling a robust estimate of stellar metallicity. Promising candidates are then confirmed with spectroscopy. We present the analysis of Magellan-MIKE high-resolution spectroscopy of 122 metal-poor stars found by SkyMapper in the first two years of commissioning observations. 41 stars have [Fe/H] <= -3.0. Nine have [Fe/H] <= -3.5, with three at [Fe/H] ~ -4. A 1D LTE abundance analysis of the elements Li, C, Na, Mg, Al, Si, Ca, Sc, Ti, Cr, Mn, Co, Ni, Zn, Sr, Ba and Eu shows these stars have [X/Fe] ratios typical of other halo stars. One star with low [X/Fe] values appears to be "Fe-enhanced," while another star has an extremely large [Sr/Ba] ratio: >2. Only one other star is known to have a comparable value. Seven stars are "CEMP-no" stars ([C/Fe] > 0.7, [Ba/Fe] < 0). 21 stars exhibit mild r-process element enhancements (0.3 <=[Eu/Fe] < 1.0), while four stars have [Eu/Fe] >= 1.0. These results demonstrate the ability to identify extremely metal-poor stars from SkyMapper photometry, pointing to increased sample sizes and a better characterization of the metal-poor tail of the halo metallicity distribution function in the future. △ Less

Submitted 13 July, 2015; v1 submitted 13 April, 2015; originally announced April 2015.

Comments: Minor corrections to text, missing data added to Tables 3 and 4; updated to match published version. Complete tables included in source

Journal ref: 2015 ApJ, 807, 171

Offsets between the X-ray and the Sunyaev-Zel'dovich-effect peaks in merging galaxy clusters and their cosmological implications

Authors: Congyao Zhang, Qingjuan Yu, Youjun Lu

Abstract: Observations reveal that the peaks of the X-ray map and the Sunyaev-Zel'dovich (SZ) effect map of some galaxy clusters are offset from each other. In this paper, we perform a set of hydrodynamical simulations of mergers of two galaxy clusters to investigate the spatial offset between the maxima of the X-ray and the SZ surface brightness of the merging clusters. We find that significantly large SZ-… ▽ More Observations reveal that the peaks of the X-ray map and the Sunyaev-Zel'dovich (SZ) effect map of some galaxy clusters are offset from each other. In this paper, we perform a set of hydrodynamical simulations of mergers of two galaxy clusters to investigate the spatial offset between the maxima of the X-ray and the SZ surface brightness of the merging clusters. We find that significantly large SZ-X-ray offsets (>100kpc) can be produced during the major mergers of galaxy clusters. The significantly large offsets are mainly caused by a `jump effect' occurred between the primary and secondary pericentric passages of the two merging clusters, during which the X-ray peak may jump to the densest gas region located near the center of the small cluster, but the SZ peak remains near the center of the large one. Our simulations show that merging systems with higher masses and larger initial relative velocities may result in larger offset sizes and longer offset time durations; and only nearly head-on mergers are likely to produce significantly large offsets. We further investigate the statistical distribution of the SZ-X-ray offset sizes and find that (1) the number distribution of the offset sizes is bimodal with one peak located at low offsets ~0 and the other at large offsets ~350-450kpc/h, but the objects with intermediate offsets are scarce; and (2) the probabilities of the clusters in the mass range higher than 2x10^{14}Msun/h that have offsets larger than 20, 50, 200, 300, and 500kpc/h are 34.0%, 11.1%, 8.0%, 6.5%, and 2.0% respectively at z=0.7. The probability is sensitive to the underlying pairwise velocity distribution and the merger rate of clusters. Future observations on the offsets for a large number of clusters may put strong constraints on the cosmic velocity fields on the cluster scale and the cluster merger rate. (Abridged) △ Less

Submitted 16 June, 2014; originally announced June 2014.

Comments: 25 pages, 15 figures

Probing baryonic processes and gastrophysics in the formation of the Milky Way dwarf satellites: I. metallicity distribution properties

Authors: Jun Hou, Qingjuan Yu, Youjun Lu

Abstract: In this paper, we study the chemical properties of the stars in the dwarf satellites around the MW-like host galaxies, and explore the possible effects of several baryonic processes, including supernova (SN) feedback, the reionization of the universe and H$_2$ cooling, on them and how current and future observations may put some constraints on these processes. We use a semi-analytical model to gen… ▽ More In this paper, we study the chemical properties of the stars in the dwarf satellites around the MW-like host galaxies, and explore the possible effects of several baryonic processes, including supernova (SN) feedback, the reionization of the universe and H$_2$ cooling, on them and how current and future observations may put some constraints on these processes. We use a semi-analytical model to generate MW-like galaxies, for which a fiducial model can reproduce the luminosity function and the stellar metallicity--stellar mass correlation of the MW dwarfs. Using the simulated MW-like galaxies, we focus on investigating three metallicity properties of their dwarfs: the stellar metallicity--stellar mass correlation of the dwarf population, and the metal-poor and metal-rich tails of the stellar metallicity distribution in individual dwarfs. We find that (1) the slope of the stellar metallicity--stellar mass correlation is sensitive to the SN feedback strength and the reionization epoch; (2) the extension of the metal-rich tails is mainly sensitive to the SN feedback strength; (3) the extension of the metal-poor tails is mainly sensitive to the reionization epoch; (4) none of the three chemical properties are sensitive to the H$_2$ cooling process; and (5) comparison of our model results with the current observational slope of the stellar metallicity--stellar mass relation suggests that the local universe is reionized earlier than the cosmic average and local sources may have a significant contribution to the reionization in the local region, and an intermediate to strong SN feedback strength is preferred. Future observations of metal-rich and metal-poor tails of stellar metallicity distributions will put further constraints on the SN feedback and the reionization processes. △ Less

Submitted 16 June, 2014; originally announced June 2014.

Comments: 22 pages, 16 figures, accepted for publication in the Astrophysical Journal

On the Existence of Pulsars in the Vicinity of the Massive Black Hole in the Galactic Center

Authors: Fupeng Zhang, Youjun Lu, Qingjuan Yu

Abstract: Pulsars, if existing and detectable in the immediate vicinity of the massive black hole (MBH) in the Galactic center (GC), may be used as a superb tool to probe both the environment and the metric of the central MBH. The recent discovery of a magnetized pulsar in the GC suggests that many more pulsars should exist near the MBH. In this paper, we estimate the number and the orbital distribution of… ▽ More Pulsars, if existing and detectable in the immediate vicinity of the massive black hole (MBH) in the Galactic center (GC), may be used as a superb tool to probe both the environment and the metric of the central MBH. The recent discovery of a magnetized pulsar in the GC suggests that many more pulsars should exist near the MBH. In this paper, we estimate the number and the orbital distribution of pulsars in the vicinity of the MBH in the GC by assuming that the pulsar progenitors, similar to the GC S-stars, were captured to orbits tightly bound to the MBH through the tidal breakup of stellar binaries. We use the current observations on both the GC S-stars and the hypervelocity stars to calibrate the injection rate(s) of and the dynamical model(s) for the stellar binaries. By including the relaxation processes, supernova kicks, and gravitational wave radiation in our simulations, we estimate that ~97-190 (9-14) pulsars may presently orbit the central MBH with semimajor axes <=4000AU (<=1000AU), which is compatible with the current observational constraints on the number of the GC pulsars. The semimajor axis and the pericenter distance of the pulsar closest to the central MBH are probably in the range of ~120-460AU and ~2-230AU, respectively. Future telescopes, such as the SKA, may be able to detect a significant number of pulsars with semimajor axis smaller than a few thousand AU in the GC. Long-term monitoring of these pulsars would be helpful in constraining both the environment and the metric of the central MBH. Our preferred model also results in about ten hyperfast pulsars with velocity >~1500km/s moving away from the Milky Way. △ Less

Submitted 11 February, 2014; originally announced February 2014.

Comments: 9 pages, 3 figures, accepted for publication in the Astrophysical Journal

Microlensing of Sub-parsec Massive Binary Black Holes in Lensed QSOs: Light Curves and Size-Wavelength Relation

Authors: Chang-Shuo Yan, Youjun Lu, Qingjuan Yu, Shude Mao, Joachim Wambsganss

Abstract: Sub-parsec binary massive black holes (BBHs) are long anticipated to exist in many QSOs but remain observationally elusive. In this paper, we propose a novel method to probe sub-parsec BBHs through microlensing of lensed QSOs. If a QSO hosts a sub-parsec BBH in its center, it is expected that the BBH is surrounded by a circum-binary disk, each component of the BBH is surrounded by a small accretio… ▽ More Sub-parsec binary massive black holes (BBHs) are long anticipated to exist in many QSOs but remain observationally elusive. In this paper, we propose a novel method to probe sub-parsec BBHs through microlensing of lensed QSOs. If a QSO hosts a sub-parsec BBH in its center, it is expected that the BBH is surrounded by a circum-binary disk, each component of the BBH is surrounded by a small accretion disk, and a gap is opened by the secondary component in between the circum-binary disk and the two small disks. Assuming such a BBH structure, we generate mock microlensing light curves for some QSO systems that host BBHs with typical physical parameters. We show that microlensing light curves of a BBH QSO system at the infrared-optical-UV bands can be significantly different from those of corresponding QSO system with a single massive black hole (MBH), mainly because of the existence of the gap and the rotation of the BBH (and its associated small disks) around the center of mass. We estimate the half-light radii of the emission region at different wavelengths from mock light curves and find that the obtained half-light radius vs. wavelength relations of BBH QSO systems can be much flatter than those of single MBH QSO systems at a wavelength range determined by the BBH parameters, such as the total mass, mass ratio, separation, accretion rates, etc. The difference is primarily due to the existence of the gap. Such unique features on the light curves and half-light radius-wavelength relations of BBH QSO systems can be used to select and probe sub-parsec BBHs in a large number of lensed QSOs to be discovered by current and future surveys, including the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS), the Large Synoptic Survey telescope (LSST) and Euclid. △ Less

Submitted 11 February, 2014; originally announced February 2014.

Comments: 18 pages, 17 figures, accepted for publication in the Astrophysical Journal

Orbital orientation evolution of massive binary black holes at the centres of non-spherical galaxies

Authors: Xiaobo Cui, Qingjuan Yu

Abstract: At the centre of a spherical and kinematically isotropic galaxy, the orientation of a massive binary black hole (BBH) orbit (i.e., the direction of the BBH orbital angular momentum) undergoes a random walk. If the stars in a spherical system have a non-zero total angular momentum, the BBH orbital orientation evolves towards aligning with the total stellar angular momentum direction. In this paper,… ▽ More At the centre of a spherical and kinematically isotropic galaxy, the orientation of a massive binary black hole (BBH) orbit (i.e., the direction of the BBH orbital angular momentum) undergoes a random walk. If the stars in a spherical system have a non-zero total angular momentum, the BBH orbital orientation evolves towards aligning with the total stellar angular momentum direction. In this paper, we show that a triaxial galaxy has an alignment-erasing effect, that is, the alignment of the BBH orientations towards the galaxy rotation axis can be decreased significantly or erased. We also show that in a non-rotating axisymmetric galaxy, the BBH orbital orientation evolves towards the axisymmetric axis and precesses about it in a retrograde direction. Our results provide a step towards understanding the spin orientations of the final merged BH (and hence probable orientation of any jet produced) within its host galaxy, and may help to constrain the recoiling velocity of the merged BH arose from gravitational wave radiation as well. △ Less

Submitted 8 October, 2013; originally announced October 2013.

Comments: 16 pages, 9 figures, MNRAS accepted

The Aquarius Co-Moving Group is Not a Disrupted Classical Globular Cluster

Authors: Andrew R. Casey, Stefan Keller, Alan Alves-Brito, Anna Frebel, Gary Da Costa, Amanda Karakas, David Yong, Kevin Schlaufman, Heather R. Jacobson, Qinsi Yu, Cherie Fishlock

Abstract: We present a detailed analysis of high-resolution, high S/N spectra for 5 Aquarius stream stars observed with the MIKE spectrograph on the Magellan Clay telescope. Our sample represents one third of the 15 known members in the stream. We find the stream is not mono-metallic: the metallicity ranges from [Fe/H] = -0.63 to -1.58. No anti-correlation in Na-O abundances is present, and we find a strong… ▽ More We present a detailed analysis of high-resolution, high S/N spectra for 5 Aquarius stream stars observed with the MIKE spectrograph on the Magellan Clay telescope. Our sample represents one third of the 15 known members in the stream. We find the stream is not mono-metallic: the metallicity ranges from [Fe/H] = -0.63 to -1.58. No anti-correlation in Na-O abundances is present, and we find a strong positive Mg-Al relationship, similar to that observed in the thick disk. We find no evidence that the stream is a result of a disrupted classical globular cluster, contrary to a previously published claim. High [(Na, Ni, alpha)/Fe] and low [Ba/Y] abundance ratios in the stream suggests it is not a tidal tail from a disrupted dwarf galaxy, either. The stream is chemically indistinguishable from Milky Way field stars with the exception of one candidate, C222531-145437. From its position, velocity, and detailed chemical abundances, C222531-145437 is likely a star that was tidally disrupted from omega-Centauri. We propose the Aquarius stream is Galactic in origin, and could be the result from a disk-satellite perturbation in the Milky Way thick disk on the order of a few Gyr ago: derived orbits, UVW velocities, and angular momenta of the Aquarius members offer qualitative support for our hypothesis. Assuming C222531-145437 is a tidally disrupted member of omega-Centauri, this system is the most likely disk perturber. In the absence of compelling chemical and/or dynamical evidence that the Aquarius stream is the tidal tail of a disrupted satellite, we advocate the "Aquarius group" as a more appropriate description. Like the Canis Major over-density, as well as the Hercules and Monoceros groups, the Aquarius group joins the list of kinematically-identified substructures that are not actually accreted material: they are simply part of the rich complexity of the Milky Way structure. △ Less

Submitted 23 May, 2014; v1 submitted 13 September, 2013; originally announced September 2013.

Comments: Accepted to MNRAS. Updated to journal version

Deriving Stellar Effective Temperatures of Metal-Poor Stars with the Excitation Potential Method

Authors: Anna Frebel, Andrew R. Casey, Heather R. Jacobson, Qinsi Yu

Abstract: It is well established that stellar effective temperatures determined from photometry and spectroscopy yield systematically different results. We describe a new, simple method to correct spectroscopically derived temperatures ("excitation temperatures") of metal-poor stars based on a literature sample with -3.3<[Fe/H]<-2.5. Excitation temperatures were determined from FeI line abundances in high-r… ▽ More It is well established that stellar effective temperatures determined from photometry and spectroscopy yield systematically different results. We describe a new, simple method to correct spectroscopically derived temperatures ("excitation temperatures") of metal-poor stars based on a literature sample with -3.3<[Fe/H]<-2.5. Excitation temperatures were determined from FeI line abundances in high-resolution optical spectra in the wavelength range of ~3700 to ~7000A, although shorter wavelength ranges, up to 4750 to 6800A, can also be employed, and compared with photometric literature temperatures. Our adjustment scheme increases the temperatures up to several hundred degrees for cool red giants, while leaving the near-main-sequence stars mostly unchanged. Hence, it brings the excitation temperatures in good agreement with photometrically derived values. The modified temperature also influences other stellar parameters, as the FeI-FeII ionization balance is simultaneously used to determine the surface gravity, while also forcing no abundance trend on the absorption line strengths to obtain the microturbulent velocity. As a result of increasing the temperature, the often too low gravities and too high microturbulent velocities in red giants become higher and lower, respectively. Our adjustment scheme thus continues to build on the advantage of deriving temperatures from spectroscopy alone, independent of reddening, while at the same time producing stellar chemical abundances that are more straightforwardly comparable to studies based on photometrically derived temperatures. Hence, our method may prove beneficial for comparing different studies in the literature as well as the many high-resolution stellar spectroscopic surveys that are or will be carried out in the next few years. △ Less

Submitted 8 April, 2013; originally announced April 2013.

Comments: 12 pages, emulateapj, accepted for publication in ApJ

The cosmic evolution of massive black holes and galaxy spheroids: global constraints at redshift z<~1.2

Authors: Xiaoxia Zhang, Youjun Lu, Qingjuan Yu

Abstract: We study the observational constraints on the cosmic evolution of the relationships between the massive black hole (MBH) mass (M_bh) and the stellar mass (M^*_sph; or velocity dispersion σ) of the host galaxy/spheroid. Assuming that the M_bh-M^*_sph (or M_bh-σ) relation evolves with redshift as \propto (1+z)^Γ, the MBH mass density can be obtained from either the observationally determined galaxy… ▽ More We study the observational constraints on the cosmic evolution of the relationships between the massive black hole (MBH) mass (M_bh) and the stellar mass (M^*_sph; or velocity dispersion σ) of the host galaxy/spheroid. Assuming that the M_bh-M^*_sph (or M_bh-σ) relation evolves with redshift as \propto (1+z)^Γ, the MBH mass density can be obtained from either the observationally determined galaxy stellar mass functions or velocity dispersion distribution functions over redshift z~0-1.2 for any given Γ. The MBH mass density at different redshifts can also be inferred from the luminosity function of QSOs/AGNs provided known radiative efficiency ε. By matching the MBH density inferred from galaxies to that obtained from QSOs/AGNs, we find that Γ=0.64^{+0.27}_{-0.29} for the M_bh-M^*_sph relation and Γ=-0.21^{+0.28}_{-0.33} for the M_bh-σrelation, and ε=0.11^{+0.04}_{-0.03}. Our results suggest that the MBH mass growth precedes the bulge mass growth but the galaxy velocity dispersion does not increase with the mass growth of the bulge after the quench of nuclear activity, which is roughly consistent with the two-phase galaxy formation scenario proposed by Oser et al. (2012) in which a galaxy roughly double its masses after z=1 due to accretion and minor mergers while its velocity dispersion drops slightly. △ Less

Submitted 15 October, 2012; originally announced October 2012.

Comments: 7 pages, 2 Figures, ApJ in press

S-stars in the Galactic center and hypervelocity stars in the Galactic halo: two faces of the tidal breakup of stellar binaries by the central massive black hole?

Authors: Fupeng Zhang, Youjun Lu, Qingjuan Yu

Abstract: In this paper, we investigate the link between the hypervelocity stars (HVSs) discovered in the Galactic halo and the S-stars moving in the Galactic center (GC), under the hypothesis that they are both the products of the tidal breakup of the same population of stellar binaries by the central massive black hole (MBH). By adopting several hypothetical models for binaries to be injected into the vic… ▽ More In this paper, we investigate the link between the hypervelocity stars (HVSs) discovered in the Galactic halo and the S-stars moving in the Galactic center (GC), under the hypothesis that they are both the products of the tidal breakup of the same population of stellar binaries by the central massive black hole (MBH). By adopting several hypothetical models for binaries to be injected into the vicinity of the MBH and doing numerical simulations, we realize the tidal breakup processes of the binaries and their follow-up evolution. We find that many statistical properties of the detected HVSs and S-stars can be reproduced under some binary injecting models, and their number ratio can be reproduced if the stellar initial mass function is top-heavy (e.g., with slope ~-1.6). The total number of the captured companions is ~50 that have masses in the range ~3-7Msun and semimajor axes <~4000 AU and survive to the present within their main-sequence lifetime. The innermost one is expected to have a semimajor axis ~300-1500 AU and a pericenter distance ~10-200 AU, with a significant probability of being closer to the MBH than S2. Future detection of such a closer star would offer an important test to general relativity. The majority of the surviving ejected companions of the S-stars are expected to be located at Galactocentric distances <~20 kpc, and have heliocentric radial velocities ~-500-1500 km/s and proper motions up to ~5-20 mas/yr. Future detection of these HVSs may provide evidence for the tidal-breakup formation mechanism of the S-stars. △ Less

Submitted 2 October, 2013; v1 submitted 5 October, 2012; originally announced October 2012.

Comments: updated to match the published version, 18 pages, 7 figures

The Rossiter-McLaughlin effect for exomoons or binary planets

Authors: Quntao Zhuang, Xun Gao, Qingjuan Yu

Abstract: In this paper we study possible signatures of binary planets or exomoons on the Rossiter-McLaughlin (R-M) effect. Our analyses show that the R-M effect for a binary planet or exomoon during its complete transit phase can be divided into two parts. The first is the conventional one similar to the R-M effect from the transit of a single planet, of which the mass and the projected area are the combin… ▽ More In this paper we study possible signatures of binary planets or exomoons on the Rossiter-McLaughlin (R-M) effect. Our analyses show that the R-M effect for a binary planet or exomoon during its complete transit phase can be divided into two parts. The first is the conventional one similar to the R-M effect from the transit of a single planet, of which the mass and the projected area are the combinations of the binary components; and the second is caused by the orbital rotation of the binary components, which may add a sine- or linear-mode deviation to the stellar radial velocity curve. We find that the latter effect can be up to several or several ten m/s. By doing numerical simulations as well as analytical analyses, we illustrate that the distribution and dispersion of the latter effects obtained from multiple transit events can be used to constrain the dynamical configuration of the binary planet, such as, how the inner orbit of the binary planet is inclined to its orbit rotating around the central star. We find that the signatures caused by the orbital rotation of the binary components are more likely to be revealed if the two components of binary planet have different masses and mass densities, especially if the heavy one has a high mass density and the light one has a low density. Similar signature on the R-M effect may also be revealed in a hierarchical triple star system containing a dark compact binary and a tertiary star. △ Less

Submitted 30 July, 2012; originally announced July 2012.

Comments: 30 pages, 9 figures

The low frequency of dual AGNs versus the high merger rate of galaxies: A phenomenological model

Authors: Qingjuan Yu, Youjun Lu, Roya Mohayaee, Jacques Colin

Abstract: Dual AGNs are natural byproducts of hierarchical mergers of galaxies in the LambdaCDM cosmogony. Recent observations have shown that only a small fraction (~ 0.1%-1%) of AGNs at redshift z<~ 0.3 are dual with kpc-scale separations, which is rather low compared to the high merger rate of galaxies. Here we construct a phenomenological model to estimate the number density of dual AGNs and its evoluti… ▽ More Dual AGNs are natural byproducts of hierarchical mergers of galaxies in the LambdaCDM cosmogony. Recent observations have shown that only a small fraction (~ 0.1%-1%) of AGNs at redshift z<~ 0.3 are dual with kpc-scale separations, which is rather low compared to the high merger rate of galaxies. Here we construct a phenomenological model to estimate the number density of dual AGNs and its evolution according to the observationally-estimated major merger rates of galaxies and various scaling relations on the properties of galaxies and their central massive black holes. We show that our model reproduces the observed frequency and separation distribution of dual AGNs provided that significant nuclear activities are triggered only in gas-rich progenitor galaxies with central massive black holes and only when the nuclei of these galaxies are roughly within the half-light radii of their companion galaxies. Under these constraints, the observed low dual AGN frequency is consistent with the relatively high merger rate of galaxies and supports the hypothesis that major mergers lead to AGN/QSO activities. We also predict that the number of kpc-scale dual AGNs decreases with increasing redshift and only about 0.02%--0.06% of AGNs are dual AGNs with double-peaked narrow line features at redshifts of z 0.5-1.2. Future observations of high-redshift dual AGNs would provide a solid test for this prediction. △ Less

Submitted 10 May, 2011; originally announced May 2011.

Comments: 8 pages, 3 figures

The spatial and velocity distributions of hypervelocity stars

Authors: Fupeng Zhang, Youjun Lu, Qingjuan Yu

Abstract: Hypervelocity stars (HVSs) found in the Galactic halo are probably the dynamical products of interactions between (binary) stars and the massive black hole(s) (MBH) in the Galactic center (GC). It has been shown that the detected HVSs are spatially consistent with being located on two thin disks (Lu et al.), one of which has the same orientation as the clockwise-rotating stellar disk in the GC. He… ▽ More Hypervelocity stars (HVSs) found in the Galactic halo are probably the dynamical products of interactions between (binary) stars and the massive black hole(s) (MBH) in the Galactic center (GC). It has been shown that the detected HVSs are spatially consistent with being located on two thin disks (Lu et al.), one of which has the same orientation as the clockwise-rotating stellar disk in the GC. Here we perform a large number of three-body experiments of the interactions between the MBH and binary stars bound to it, and find that the probability of ejecting HVSs is substantially enhanced by multiple encounters between the MBH and binary stars at distances substantially larger than the initial tidal breakup radii. Assuming that the HVS progenitors are originated from the two disks, the inclination distribution of the HVSs relative to the disk planes can be reproduced by either the mechanism of tidal breakup of binary stars or the mechanism of ejecting HVSs by a hypothetical binary black hole (BBH) in the GC. However, an isotropical origination of HVS progenitors is inconsistent with the observed inclination distribution. Assuming that the HVSs were ejected out by the tidal breakup mechanism, its velocity distribution can be reproduced if their progenitors diffuse onto low angular momentum orbits slowly and most of the progenitors were broken up at relatively large distances due to multiple encounters. Assuming that the HVSs were ejected out by a BBH within the allowed parameter space in the GC, our simulations produce relatively flatter velocity spectra compared to the observed ones; however, the BBH mechanism cannot be statistically ruled out, yet. Future surveys of HVSs and better statistics of their spatial and velocity distributions should enable to distinguish the ejection mechanisms of HVSs and shed new light on the dynamical environment of the MBH.(abridged) △ Less

Submitted 7 May, 2011; originally announced May 2011.

Comments: 19 pages, 16 figures

Journal ref: Astrophys.J. 722, 1744-1761 (2010)

On probing the properties of QSOs through their proximity effects on the intergalactic medium

Authors: Youjun Lu, Qingjuan Yu

Abstract: The proximity effect (PE) of QSOs is believed to be useful in constraining the QSO lifetime. Observations on the PE so far, however, give apparently contradictory results -- some are consistent with a long QSO lifetime (>~ a few 10^7 yr), but others appear to be only consistent with a short QSO lifetime <~ 10^6 yr. In this paper, we show that this apparent contradiction may be solved by simultaneo… ▽ More The proximity effect (PE) of QSOs is believed to be useful in constraining the QSO lifetime. Observations on the PE so far, however, give apparently contradictory results -- some are consistent with a long QSO lifetime (>~ a few 10^7 yr), but others appear to be only consistent with a short QSO lifetime <~ 10^6 yr. In this paper, we show that this apparent contradiction may be solved by simultaneously taking into account both the effect due to the density enhancement in the QSO near zones and that due to the obscuration of the tori associated with the QSOs, using a large number of Monte-Carlo generated synthetic Lyman alpha forest spectra. We demonstrate that the QSO properties and environment can be constrained simultaneously by the transverse PE and the line of sight PE of bright type 1 QSOs together. The current available measurements on the PEs of type 1 QSOs suggest that (1) the density is significantly enhanced in the vicinity of the QSOs; (2) the QSO lifetime is consistent with being as large as a few 10^7 yr and a substantially shorter lifetime (e.g., <~10^6 yr) is not required; and (3) the half opening angle of the tori associated with QSOs is ~60 deg, consistent with some other independent estimates. Our simulations also show that the TPE of type 2 QSOs can be significantly different from that of type 1 QSOs, which may be useful to put further constraints on the QSO properties and the QSO environment. △ Less

Submitted 7 May, 2011; originally announced May 2011.

Comments: 14 pages, 11 figures, to appear in ApJ

Investigation of excited states in $^{18}$Ne via resonant elastic scattering of $^{17}$F+p and its astrophysical implication in the stellar reaction of $^{14}$O($α$,$p$)$^{17}$F

Authors: J. Hu, J. J. He, S. W. Xu, Z. Q. Chen, X. Y. Zhang, J. S. Wang, X. Q. Yu, L. Li, L. Y. Zhang, Y. Y. Yang, P. Ma, X. H. Zhang, Z. G. Hu, Z. Y. Guo, X. Xu, X. H. Yuan, W. Lu, Y. H. Yu, Y. D. Zang, S. W. Tang, R. P. Ye, J. D. Chen, S. L. Jin, C. M. Du, S. T. Wang , et al. (13 additional authors not shown)

Abstract: Properties of proton resonances in $^{18}$Ne have been investigated efficiently by utilizing a technique of proton resonant elastic scattering with a $^{17}$F radioactive ion (RI) beam and a thick proton target. A 4.22~MeV/nucleon $^{17}$F RI beam was produced via a projectile-fragmentation reaction, and subsequently separated by a Radioactive Ion Beam Line in Lanzhou ({\tt RIBLL}). Energy spectra… ▽ More Properties of proton resonances in $^{18}$Ne have been investigated efficiently by utilizing a technique of proton resonant elastic scattering with a $^{17}$F radioactive ion (RI) beam and a thick proton target. A 4.22~MeV/nucleon $^{17}$F RI beam was produced via a projectile-fragmentation reaction, and subsequently separated by a Radioactive Ion Beam Line in Lanzhou ({\tt RIBLL}). Energy spectra of the recoiled protons were measured by two sets of $Δ$E-E silicon telescope at center-of-mass scattering angles of $θ_{c.m.}$$\approx$175${^\circ}$$\pm$5${^\circ}$, $θ_{c.m.}$$\approx$152${^\circ}$$\pm$8${^\circ}$, respectively. Several proton resonances in $^{18}$Ne were observed, and their resonant parameters have been determined by an $R$-matrix analysis of the differential cross sections in combination with the previous results. The resonant parameters are related to the reaction-rate calculation of the stellar $^{14}$O($α$,$p$)$^{17}$F reaction, which was thought to be the breakout reaction from the hot CNO cycles into the $rp$-process in x-ray bursters. Here, $J^π$=(3$^-$, 2$^-$) are tentatively assigned to the 6.15-MeV state which was thought the key 1$^-$ state previously. In addition, a doublet structure at 7.05 MeV are tentatively identified, and its contribution to the resonant reaction rate of $^{14}$O($α$,$p$)$^{17}$F could be enhanced by at least factors of about 4$\sim$6 in comparison with the previous estimation involving only a singlet. The present calculated resonant rates are much larger than those previous values, and it may imply that this breakout reaction could play a crucial role under x-ray bursters conditions. △ Less

Submitted 7 September, 2010; originally announced September 2010.

Comments: 8 figures, 2 tables

Reexamination of reaction rates for a key stellar reaction of $^{14}$O($α$,p)$^{17}$F

Authors: J. J. He, H. W. Wang, J. Hu, L. Li, L. Y. Zhang, M. L. Liu, S. W. Xu, X. Q. Yu

Abstract: The reaction rates of the key stellar reaction of $^{14}$O($α$,p)$^{17}$F have been reexamined. The previous conclusion, the 6.15-MeV state ($J^π$=1$^-$) dominating this reaction rate, has been overthrown by a careful reanalysis of the previous experimental data [J. Gómez del Campo {\it et al.}, Phys. Rev. Lett. {\bf 86}, 43 (2001)]. According to the present $R$-matrix analysis, the previous 1… ▽ More The reaction rates of the key stellar reaction of $^{14}$O($α$,p)$^{17}$F have been reexamined. The previous conclusion, the 6.15-MeV state ($J^π$=1$^-$) dominating this reaction rate, has been overthrown by a careful reanalysis of the previous experimental data [J. Gómez del Campo {\it et al.}, Phys. Rev. Lett. {\bf 86}, 43 (2001)]. According to the present $R$-matrix analysis, the previous 1$^-$ assignment for the 6.15-MeV state is definitely wrong. Most probably, the 6.286-MeV state is the 1$^-$ state and the 6.15-MeV state is a 3$^-$ one, and hence the resonance at $E_x$=6.286 MeV ($J^π$=1$^-$) dominates the reaction rates in the temperature region of astrophysical interests. The newly calculated reaction rates for the $^{14}$O($α$,p)$^{17}$F reaction are quite different from the previous ones, for instance, it's only about 1/6 of the previous value around 0.4 GK, while it's about 2.4 times larger than the previous value around 2 GK. The astrophysical implications have been discussed based on the present conclusions. △ Less

Submitted 12 January, 2010; originally announced January 2010.

Comments: 4 pages, 2 figures, 1 table

On the spatial distribution and the origin of hypervelocity stars

Authors: Youjun Lu, Fupeng Zhang, Qingjuan Yu

Abstract: Hypervelocity stars (HVSs) escaping away from the Galactic halo are dynamical products of interactions of stars with the massive black hole(s) (MBH) in the Galactic Center (GC). They are mainly B-type stars with their progenitors unknown. OB stars are also populated in the GC, with many being hosted in a clockwise-rotating young stellar (CWS) disk within half a parsec from the MBH and their form… ▽ More Hypervelocity stars (HVSs) escaping away from the Galactic halo are dynamical products of interactions of stars with the massive black hole(s) (MBH) in the Galactic Center (GC). They are mainly B-type stars with their progenitors unknown. OB stars are also populated in the GC, with many being hosted in a clockwise-rotating young stellar (CWS) disk within half a parsec from the MBH and their formation remaining puzzles. In this paper, we demonstrate that HVSs can well memorize the injecting directions of their progenitors using both analytical arguments and numerical simulations, i.e., the ejecting direction of an HVS is almost anti-parallel to the injecting direction of its progenitor. Therefore, the spatial distribution of HVSs maps the spatial distribution of the parent population of their progenitors directly. We also find that almost all the discovered HVSs are spatially consistent with being located on two thin disk planes. The orientation of one plane is consistent with that of the (inner) CWS disk, which suggests that most of the HVSs originate from the CWS disk or a previously existed disk-like stellar structure with an orientation similar to it. The rest of HVSs may be correlated with the plane of the northern arm of the mini-spiral in the GC or the plane defined by the outer warped part of the CWS disk. Our results not only support the GC origin of HVSs but also imply that the central disk (or the disk structure with a similar orientation) should persist or be frequently rejuvenated over the past 200 Myr, which adds a new challenge to the stellar disk formation and provides insights to the longstanding problem of gas fueling into massive black holes. △ Less

Submitted 19 October, 2009; originally announced October 2009.

Comments: 7 pages, 5 figures

Journal ref: Astrophys.J.709:1356-1361,2010

Toward precise constraints on growth of massive black holes

Authors: Qingjuan Yu, Youjun Lu

Abstract: Growth of massive black holes (MBHs) in galactic centers comes mainly from gas accretion during their QSO/AGN phases. In this paper we apply an extended Soltan argument, connecting the local MBH mass function with the time-integral of the QSO luminosity function, to the demography of MBHs and QSOs from recent optical and X-ray surveys, and obtain robust constraints on the luminosity evolution (o… ▽ More Growth of massive black holes (MBHs) in galactic centers comes mainly from gas accretion during their QSO/AGN phases. In this paper we apply an extended Soltan argument, connecting the local MBH mass function with the time-integral of the QSO luminosity function, to the demography of MBHs and QSOs from recent optical and X-ray surveys, and obtain robust constraints on the luminosity evolution (or mass growth history) of individual QSOs (or MBHs). We find that the luminosity evolution probably involves two phases: an initial exponentially increasing phase set by the Eddington limit and a following phase in which the luminosity declines with time as a power law (with a slope of -1.2--1.3) set by a self-similar long-term evolution of disk accretion. Neither an evolution involving only the increasing phase with a single Eddington ratio nor an exponentially declining pattern in the second phase is likely. The period of a QSO radiating at a luminosity higher than 10% of its peak value is about (2-3)x10^8 yr, during which the MBH obtains ~80% of its mass. The mass-to-energy conversion efficiency is $0.16\pm0.04 ^{+0.05}_{-0}$, with the latter error accounting for the maximum uncertainty due to Compton-thick AGNs. The expected Eddington ratios in QSOs from the constrained luminosity evolution cluster around a single value close to 0.5-1 for high-luminosity QSOs and extend to a wide range of lower values for low-luminosity ones. The Eddington ratios for high luminosity QSOs appear to conflict with those estimated from observations (~0.25) by using some virial mass estimators for MBHs in QSOs unless the estimators systematically over-estimate MBH masses by a factor of 2-4. We also infer the fraction of optically obscured QSOs ~60-80%. Further applications of the luminosity evolution of individual QSOs are also discussed. △ Less

Submitted 27 August, 2008; originally announced August 2008.

Comments: 25 pages, 13 figures, ApJ in press

Journal ref: Astrophys.J.689:732-754,2008

Hypervelocity binary stars: smoking gun of massive binary black holes

Authors: Youjun Lu, Qingjuan Yu, D. N. C. Lin

Abstract: The hypervelocity stars recently found in the Galactic halo are expelled from the Galactic center through interactions between binary stars and the central massive black hole or between single stars and a hypothetical massive binary black hole. In this paper, we demonstrate that binary stars can be ejected out of the Galactic center with velocities up to 10^3 km/s, while preserving their integri… ▽ More The hypervelocity stars recently found in the Galactic halo are expelled from the Galactic center through interactions between binary stars and the central massive black hole or between single stars and a hypothetical massive binary black hole. In this paper, we demonstrate that binary stars can be ejected out of the Galactic center with velocities up to 10^3 km/s, while preserving their integrity, through interactions with a massive binary black hole. Binary stars are unlikely to attain such high velocities via scattering by a single massive black hole or through any other mechanisms. Based on the above theoretical prediction, we propose a search for binary systems among the hypervelocity stars. Discovery of hypervelocity binary stars, even one, is a definitive evidence of the existence of a massive binary black hole in the Galactic center. △ Less

Submitted 21 July, 2007; v1 submitted 12 July, 2007; originally announced July 2007.

Comments: 5 pages, 3 figures, shortened version, ApJL in press

Journal ref: Astrophys.J.666:L89-L92,2007