Search

Thanks to the rapidly increasing time-domain facilities, we are entering a golden era of research on gamma-ray bursts (GRBs). In this Letter, we report our observations of GRB 240529A with the Burst Optical Observer and Transient Exploring System, the 1.5-meter telescope at Observatorio Sierra Nevada, the 2.5-meter Wide Field Survey Telescope of China, the Large Binocular Telescope, and the Telescopio Nazionale Galileo. The prompt emission of GRB 240529A shows two comparable energetic episodes separated by a quiescence time of roughly 400 s. Combining all available data on the GRB Coordinates Network, we reveal the simultaneous apparent X-ray plateau and optical re-brightening around $10^3-10^4$ s after the burst. Rather than the energy injection from the magnetar as widely invoked for similar GRBs, the multi-wavelength emissions could be better explained as two shocks launched from the central engine separately. The optical peak time and our numerical modeling suggest that the initial bulk Lorentz factor of the later shock is roughly 50, which indicates that the later jet should be accretion-driven and have a higher mass loading than a typical one. The quiescence time between the two prompt emission episodes may be caused by the transition between different accretion states of a central magnetar or black hole, or the fall-back accretion process. A sample of similar bursts with multiple emission episodes in the prompt phase and sufficient follow-up could help to probe the underlying physics of GRB central engines., Comment: Resubmitted to ApJL after addressing the referee's comments; comments are welcome

Magnetars are young, highly magnetized neutron stars that are associated with magnetar short bursts (MSBs), magnetar giant flares (MGFs), and at least some fast radio bursts (FRBs). In this work, we consider a magnetar and a main sequence star in a binary system and analyze the properties of the electromagnetic signals generated by the interaction between the magnetar bursts and the companion star. During the preburst period, persistent radiation could be generated by the interaction between the $e^+e^-$-pair wind from the magnetar and the companion or its stellar wind. We find that for a newborn magnetar, the persistent preburst radiation from the strong magnetar wind can be dominant, and it is mainly at the optical and ultraviolet (UV) bands. For relatively old magnetars, the re-emission from a burst interacting with the companion is larger than the persistent preburst radiation and the luminosity of the companion itself. The transient re-emission produced by the heating process has a duration of $0.1 - 10^5 {\rm~s}$ at the optical, UV, and X-ray bands. Additionally, we find that if these phenomena occur in nearby galaxies within a few hundred kiloparsecs, they could be detected by current or future optical telescopes., Comment: 15 pages, 8 figures; Published in A&A

Gamma-ray bursts (GRBs) have long been considered potential sources of ultra-high-energy cosmic rays (UHECRs; with energy $\gtrsim 10^{18} {\rm~eV}$). In this work, we propose a novel model generating MeV emission lines in GRB, which can constrain the properties of heavy nuclei that potentially exist in GRB jets. Specifically, we find that relativistic hydrogen-like high-atomic-number ions originating from the $\beta$ decay of unstable nuclei and/or the recombination entrained in the GRB jet can generate narrow MeV emission lines through the de-excitation of excited-electrons. This model can successfully explain the MeV emission line observed in the most luminous GRB ever recorded, GRB~221009A, with suitable parameters including a Lorentz factor $\gamma \sim 820-1700$ and a total mass of heavy nuclei $M_{\rm tot} \sim 10^{23} - 10^{26}$~g. Especially, the emission line broadening can be reasonably attributed to both the expansion of the jet shell and the thermal motion of nuclei, naturally resulting in a narrow width ($\sigma_{\rm line} / E_{\rm line} \lesssim 0.2$) consistent with the observation. Furthermore, we predict that different GRBs can exhibit lines in different bands with various evolving behaviors, which might be confirmed with further observations. Finally, our model provides indirect evidence that GRBs may be one of the sources of UHECRs., Comment: 13 pages, 4 figures; Published in ApJL, https://doi.org/10.3847/2041-8213/ad4ce1

Rapidly spinning magnetars are potential candidates for the energy source of supernovae (SNe) and gamma-ray bursts and the most promising sources for continuous gravitational waves (GWs) detected by ground-based GW detectors. Continuous GWs can be radiated from magnetars due to magnetic-induced deformation or fluid oscillations, compatible with magnetic dipole (MD) radiation for spin-down energy. In this paper, we investigate the diverse light curves of magnetar-driven SNe in the scenario that the spin-down is dominated by GW radiation and/or MD radiation. By simulating the light curves of SNe and employing the Markov chain Monte Carlo method, we constrain the parameters of the magnetars and SN explosions and show that the signature of GW radiation may be indicated by the bolometric luminosity curves of SNe Ic-BL 2007ru and 2009bb. We find that the ellipticity of magnetars in the order of $10^{-3}$ can be induced by the magnetic field of $\sim 10^{16} \mbox{ G}$. If such continuous GWs associated with SNe can be detected in the future by the Advanced LIGO and Virgo detectors, this would be a smoking gun for a magnetar engine powering SNe., Comment: 13pages, 9 figures, Accepted for publication in ApJ

In 2019, Neutron star Interior Composition ExploreR (NICER) mission released its findings on the mass and radius of the isolated neutron star (INS) PSR J0030+0451, revealing a mass of approximately 1.4 solar masses ($M_{\odot}$) and a radius near 13 kilometers. However, the recent re-analysis by the NICER collaboration \citep{vinciguerra2024updated} suggests that the available data primarily yields a precise inference of the compactness for this source while the resulting mass and radius are strongly model-dependent and diverse (the 68.3\% credible regions just overlap slightly for the ST+PDT and PDT-U models). By integrating this compactness data with the equation of state (EoS) refined by our latest investigations, we have deduced the mass and radius for PSR J0030+0451, delivering estimates of $M=1.48^{+0.09}_{-0.10}~M_\odot$ and $R=12.38_{-0.70}^{+0.51}~{\rm km}$ for the compactness found in ST+PDT model, alongside $M=1.47^{+0.14}_{-0.20}~M_\odot$ and $R=12.37_{-0.69}^{+0.50}~{\rm km}$ for the compactness in PDT-U model. These two groups of results are well consistent with each other and the direct X-ray data inference within the ST+PDT model seems to be favored. Additionally, we have calculated the tidal deformability, moment of inertia, and gravitational binding energy for this NS. Furthermore, employing these refined EoS models, we have updated mass-radius estimates for three INSs with established gravitational redshifts., Comment: 10 pages,6 figures, Accepted for publication in ApJ

The measurements of the bulk properties of most isolated neutron stars (INSs) are challenging tasks. Tang et al. (2020) have developed a new method, based on the equation of state (EoS) of neutron star (NS) material constrained by the observational data, to infer the gravitational masses of a few INSs whose gravitational redshifts are available. However, in that work, the authors only considered the constraints on the EoS from nuclear experiments/theories and the gravitational wave data of GW170817; the possible phase transition has not been taken into account. In this work, we adopt three EoS models (including the one incorporates a first-order strong phase transition) that are constrained by the latest multimessenger NS data, including in particular the recent mass\textendash radius measurements of two NSs by Neutron Star Interior Composition Explorer, to update the estimation of the gravitational masses of RBS 1223, RX J0720.4-3125, and RX J1856.5-3754. In comparison to our previous approach, the new constraints are tighter, and the gravitational masses are larger by about $0.1M_\odot$. All the inferred gravitational masses are within the range of the NS masses measured in other ways. We have also calculated the radius, tidal-deformability, and moment of inertia of these sources. The inclusion of the first-order strong phase transition has little influence on modifying the results., Comment: 8pages,4figures

A stellar-mass black hole (BH) or a millisecond magnetar is believed to be born as the central engine of Gamma-ray bursts (GRBs). The presence of plateaus in the X-ray extended emission or afterglow of GRBs is widely accepted as an indicator of magnetar central engine, particularly those with a sharp decay (faster than $t^{-3}$), so-called internal plateau. However, an alternative model, by taking the evolution of the magnetic flux at the BH horizon into account, suggests that an internal plateau can also arise from a Blandford-Znajek (BZ) mechanism powered jet (hereafter referred to as the BZ jet). In this study, we propose that a precessional BZ jet would manifest a Quasi-Periodic Oscillation (QPO) signature on the internal plateau and the subsequent sharp decay. Such lightcurves cannot be readily explained by the activity of a short-lived, supermassive magnetar, thus favoring a Kerr BH as the central engine. The X-ray afterglow of GRB 050904, comprising nine flares, is characterized by a QPO-modulated plateau and sharp decay, which can be well reproduced by a precessional BZ jet model. Therefore, one potential clue for distinguishing between these two engines lies in whether QPO signature throughout the entire plateau and the subsequent sharp decay, as the magnetar scenario suggests a collapse at the end of the plateau., Comment: Accepted for publication in ApJ

Direct mechanical coupling is known to be critical for establishing synchronization among cilia. However, the actual role of the connections is still elusive - partly because controlled experiments in live samples are challenging. Here, we employ an artificial ciliary system to address this issue. Two cilia are formed by chains of self-propelling robots and anchored to a shared base so that they are purely mechanically-coupled. The system mimics biological ciliary beating but allows fine control over the beating dynamics. We find that the artificial cilia exhibit rich motion behaviors, depending on the mechanical coupling scheme. Particularly, their synchronous beating display two distinct modes - analogous to those observed in C. reinhardtii, the biciliated model organism for studying synchronization. Close examination suggests that the system evolves towards the most dissipative mode. Using this guideline in both simulations and experiments, we are able to direct the system into a desired state by altering the modes' respective dissipation. Our results have significant implications in understanding the synchronization of cilia.

Recently, an association of GW190425 and FRB 20190425A had been claimed and a highly magnetized neutron star (NS) remnant was speculated. Given the $\sim 2.5$-h delay of the occurrence of FRB 20190425A, a uniformly rotating supramassive magnetar is favored since the differential rotation would have been promptly terminated by the magnetic braking. The required maximum gravitational mass ($M_{\rm TOV}$) of the nonrotating NS is $\approx 2.77M_\odot$, which is strongly in tension with the relatively low $M_{\rm TOV}\approx 2.25M_\odot$ obtained in current equation of state (EOS) constraints incorporating perturbative quantum chromodynamics (pQCD) information. However, the current mass-radius and mass-tidal deformability measurements of NSs alone do not convincingly exclude the high $M_{\rm TOV}$ possibility. By performing EOS constraints with mock measurements, we find that with a $2\%$ determination for the radius of PSR J0740+6620-like NS it is possible to distinguish between the low and high $M_{\rm TOV}$ scenarios. We further explore the prospect to resolve the issue of the appropriate density to impose the pQCD constraints with future massive NS observations or determinations of $M_{\rm TOV}$ and/or $R_{\rm TOV}$. It turns out that measuring the radius of a PSR J0740+6620-like NS is insufficient to probe the EOSs around 5 nuclear saturation density, where the information from pQCD becomes relevant. The additional precise $M_{\rm TOV}$ measurements anyhow could provide insights into the EOS at such a density. Indeed, supposing the central engine of GRB 170817A is a black hole formed via the collapse of a supramassive NS, the resulting $M_{\rm TOV}\approx 2.2M_\odot$ considerably softens the EOS at the center of the most massive NS, which is in favor of imposing the pQCD constraint at density beyond the one achievable in the NSs., Comment: 12 pages, 5 figures, PRD published

The long gamma-ray burst GRB 191221B has abundant observations in X-ray, optical and radio bands. In the literature, the observed optical light curve of GRB 191221B displays a plateau around 0.1-day, which is rather peculiar in gamma-ray bursts. Here we performed detailed analysis of the observational data from Swift/UVOT, VLT and LCO, obtained the light curve of the multi-band afterglow of GRB 191221B. By examining optical, ultraviolet, X-ray, and radio data for this event, we demonstrate that an on-axis two-component jet model can explain the observations. Our analysis suggests that the narrow component has an initial Lorentz factor of 400 and a jet opening half-angle of $1.4^{\circ}$, while the wide component has an initial Lorentz factor of 25 and a jet opening half-angle of $2.8^{\circ}$. The narrow jet dominates the early decay, whereas the wider jet causes the optical plateau and dominates late decay. According to this model, the reason for the absence of the X-ray plateau is due to the steeper spectral index of the wide component, resulting in a less significant flux contribution from the wide jet in the X-ray bands than in the optical bands. Moreover, we have explained the inconsistency in the decay indices of the UVOT and Rc-band data around 2000 seconds using reverse shock emission., Comment: 10 pages, 7 figures, Accepted by the Monthly Notices of the Royal Astronomical Society

The mass ($M_{\rm TOV}$) and radius ($R_{\rm TOV}$) of the maximum-mass nonrotating neutron star (NS) play a crucial role in constraining the elusive equation of state of cold dense matter and in predicting the fate of remnants from binary neutron star (BNS) mergers. In this study, we introduce a novel method to deduce these parameters by examining the mergers of second-generation (2G) black holes (BHs) with NSs. These 2G BHs are assumed to originate from supramassive neutron stars (SMNSs) formed in BNS mergers. Since the properties of the remnant BHs arising from the collapse of SMNSs follow a universal relation governed by $M_{\rm TOV}$ and $R_{\rm TOV}$, we anticipate that by analyzing a series ($\sim 100$ detections) of mass and spin measurements of the 2G BHs using the third-generation ground-based gravitational-wave detectors, $M_{\rm TOV}$ and $R_{\rm TOV}$ can be determined with a precision of $\sim 0.01M_\odot$ and $\sim 0.6$ km, respectively., Comment: 11 pages, 5 figures, ApJ published

For bright transients such as Gamma-Ray Bursts (GRBs), the Ultra-Violet/Optical Telescope (UVOT) operates under event mode at early phases, which records incident positions and arrival time for each photon. The event file is able to be screened into many exposures to study the early light curve of GRBs with a high time resolution, including in particular the rapid brightening of the UV/Optical emission. Such a goal, however, is hampered for some extremely bright GRBs by the saturation in UVOT event images. For moderately saturated UVOT sources, in this work we develop the method proposed in Jin et al. (2023) to recover their photometries. The basic idea is to assume a stable point spread function (PSF) of UVOT images, for which the counts in the core region (i.e., an aperture of a radius of 5 arcsec) and the wing region (i.e., an annulus ranging from 15 arcsec to 25 arcsec) should be a constant and the intrinsic flux can be reliably inferred with data in the ring. We demonstrate that in a given band, a tight correlation does hold among the background-removed count rates in the core and the wing. With the new method, the bright limit of measuring range for UVOT V and B bands increases ~ 1.7 mag, while only ~ 0.7 mag for U band due to the lack of bright calibration sources. Systematic uncertainties are ~ 0.2 mag for V, B and U bands., Comment: 27 pages, 19 figures, 3 tables, accepted to ApJS

Turbulent flows are observed in low-Reynolds active fluids. They are intrinsically different from the classical inertial turbulence and behave distinctively in two- and three-dimensions. Understanding the behaviors of this new type of turbulence and their dependence on the system dimensionality is a fundamental challenge in non-equilibrium physics. We experimentally measure flow structures and energy spectra of bacterial turbulence between two large parallel plates spaced by different heights $H$. The turbulence exhibits three regimes as H increases, resulting from the competition of bacterial length, vortex size and H. This is marked by two critical heights ($H_0$ and $H_1$) and a $H^{0.5}$ scaling law of vortex size in the large-$H$ limit. Meanwhile, the spectra display distinct universal scaling laws in quasi-two-dimensional (2D) and three-dimensional (3D) regimes, independent of bacterial activity, length and $H$, whereas scaling exponents exhibit transitions in the crossover. To understand the scaling laws, we develop a hydrodynamic model using image systems to represent the effect of no-slip confining boundaries. This model predicts universal 1 and -4 scaling on large and small length scales, respectively, and -2 and -1 on intermediate length scales in 2D and 3D, respectively, which are consistent with the experimental results. Our study suggests a framework for investigating the effect of dimensionality on non-equilibrium self-organized systems.

Studies have proposed that there is evidence for cosmological coupling of black holes (BHs) with an index of $k\approx 3$; hence, BHs serve as the astrophysical source of dark energy. However, the data sample is limited for the redshifts of $\leq 2.5$. In recent years, the James Webb Space Telescope (JWST) has detected many high-redshift active galactic nuclei (AGNs) and quasars. Among the JWST NIRSpec-/NIRCam-resolved AGNs, three are determined to be in early-type host galaxies with a redshift of $z\sim 4.5--7$. However, their $M_{\star}$ and $M_{\rm BH}$ are in tension with the predicted cosmological coupling of black holes with $k = 3$ at a confidence level of $\sim 2\sigma$, which challenges the hypothesis that BHs serve as the origin of dark energy. Future work on high-redshift AGNs using the JWST will further assess such a hypothesis by identifying more early-type host galaxies in the higher mass range., Comment: 9 pages, 3 figures, 2 tables; Comments are welcome!

A simple prediction of the well-known unification model of active galactic nuclei is that a sample of sources should exhibit an anti-correlation between the solid angle of the dusty torus and of the ionization cone (as the sum of them shall equal 4$\pi$), which however has never been detected. In this work, we analyze the correlation between [OIII] 5007 narrow emission line equivalent width and $L_{\rm IR}(\lambda)/L_{\rm bol}$ for a large sample of luminous quasars. For the first time, we detect a clear intrinsic anti-correlation between them, which immediately verifies the torus/ionization-cone geometry in luminous quasars. More interestingly, the anti-correlation significantly weakens with increasing wavelength from $\sim$ 2 to 12 $\mu$m, and disappears at $\sim$ 12 $\mu$m. Simulations show a cool dust component (in addition to equatorial torus) with its strength positively correlating with the solid angle of the ionization cone is required to explain the observations. This shows that the polar dust seen in nearby active galaxies also exists in luminous quasars, with its contribution to total dust emission increasing with $\lambda$ (from $\sim$ 2 to 12 $\mu$m) and reaching between 39%-62% (model dependent) at rest frame 12 $\mu$m. Our findings provide a unique approach to map the otherwise spatially unresolvable inner structure of quasars., Comment: 10 pages, 8 figures, accepted for publication in MNRAS

Eukaryotes swim with coordinated flagellar (ciliary) beating and steer by fine-tuning the coordination. The model organism for studying flagellate motility, C. reinhardtii (CR), employs synchronous, breast-stroke-like flagellar beating to swim, and it modulates the beating amplitudes differentially to steer. This strategy hinges on both inherent flagellar asymmetries (e.g. different response to chemical messengers) and such asymmetries being effectively coordinated in the synchronous beating. In CR, the synchrony of beating is known to be supported by a mechanical connection between flagella, however, how flagellar asymmetries persist in the synchrony remains elusive. For example, it has been speculated for decades that one flagellum leads the beating, as its dynamic properties (i.e. frequency, waveform, etc.) appear to be copied by the other one. In this study, we combine experiments, computations, and modeling efforts to elucidate the roles played by each flagellum in synchronous beating. With a non-invasive technique to selectively load each flagellum, we show that the coordinated beating essentially responds to only load exerted on the cis flagellum; and that such asymmetry in response derives from a unilateral coupling between the two flagella. Our results highlight a distinct role for each flagellum in coordination and have implication for biflagellates tactic behaviors.

Hyperluminous optical-ultraviolet fares have been detected in gamma-ray bursts and the luminosity record was held by naked-eye event GRB 080319B. Such fares are widely attributed to internal shock or external reverse shock radiation. Here, with a new method developed to derive reliable photometry from saturated sources of Swift/UVOT, we carry out time-resolved analysis of the initial white-band 150s exposure of GRB 220101A, a burst at a redshift of 4.618, and report a rapidly evolving optical-ultraviolet fare with a high absolute AB magnitude of -39.4$\pm$0.2. In contrast to GRB 080319B, the temporal behaviour of this new fare does not trace the gamma-ray activity. Instead of either internal shocks or reverse shock, this extremely energetic optical-ultraviolet fare is most likely to originate from the refreshed shocks induced by the late-ejected extremely energetic material catching up with the earlier-launched decelerating outflow. This finding reveals the diverse origins of the extremely energetic optical-ultraviolet fares and demonstrates the necessity of high-time-resolution observations at early times., Comment: Published in Nature Astronomy, the definitive version is available at https://www.nature.com/articles/s41550-023-02005-w

The temporal behavior of the very dim optical afterglow of GRB 080503 is at odds with the regular forward shock afterglow model and a sole kilonova component responsible for optical emission has been speculated in some literature. Here we analyze the optical afterglow data available in archive and construct time-resolved spectra. The significant detection by Keck-I in {\it G/R} bands at $t\sim 3$ day, which has not been reported before, as well as the simultaneous Gemini-North {\it r} band measurement, are in favor of a power-law spectrum that is well consistent with the optical to X-ray spectrum measured at $t\sim 4.5$ day. However, for $t\leq 2$ day, the spectra are thermal-like and a straightforward interpretation is a kilonova emission from a neutron star merger, making it, possibly, the first detection of a very early kilonova signal at $t\sim 0.05$ day. A non-thermal nature of optical emission at late times ($t\sim 2$ day), anyhow, can not be ruled out because of the large uncertainty of the {\it g}-band data. We also propose to classify the neutron star merger induced optical transients, according to the temporal behaviors of the kilonova and the non-thermal afterglow emission, into four types. GRB 080503 would then represent the first observation of a sub-group of neutron star merger driven optical transients (i.e., Type IV) consisting of an early blue kilonova and an adjacent non-thermal afterglow radiation., Comment: Accepted to ApJ. 15 pages, 4 figures, 5 tables

The LIGO, Virgo, and KAGRA (LVK) collaboration has announced 90 coalescing binary black holes (BBHs) with $p_{\rm astro} > 50\%$ to date, however, the origin of their formation channels is still an open scientific question. Given various properties of BBHs (BH component masses and individual spins) inferred using the default priors by the LVK, independent groups have been trying to explain the formation of the BBHs with different formation channels. Of all formation scenarios, the chemically homogeneous evolution (CHE) channel has stood out with distinguishing features, namely, nearly-equal component masses and preferentially high individual spins aligned with the orbital angular momentum. We perform Bayesian inference on the BBH events officially reported in GWTC-3 with astrophysically-predicted priors representing different formation channels of the isolated binary evolution (CEE: common-envelope evolution channel; CHE; SMT: stable mass transfer). Given assumed models, we report strong evidence for GW190517\_055101 being most likely to have formed through the CHE channel. Assuming the BBH events in the subsample are all formed through one of the isolated binary evolution channels, we obtain the lower limits on the local merger rate density of these channels at $11.45 ~\mathrm{Gpc^{-3}~yr^{-1}}$ (CEE), $0.18 ~\mathrm{Gpc^{-3}~yr^{-1}}$ (CHE), and $0.63 ~\mathrm{Gpc^{-3}~yr^{-1}}$ (SMT) at $90\%$ credible level., Comment: 13 pages, 4 figures, 1 table

Hawking's area theorem is one of the fundamental laws of black holes (BHs), which has been tested at a confidence level of $\sim 95\%$ with gravitational wave (GW) observations by analysing the inspiral and ringdown portions of GW signals, independently. In this work, we propose to carry out the test in a new way with the hierarchical triple merger (i.e., two successive BH mergers occurred sequentially within the observation window of GW detectors), for which the properties of the progenitor BHs and the remnant BH of the first coalescence can be inferred from the inspiral portions of the two mergers. As revealed in our simulations, the BH area law can be well confirmed for some plausible hierarchical triple merger events detectable in LIGO/Virgo/KAGRA/LIGO-India's O4/O5 runs. Our proposed method provides significant facilitation for testing the area law and complements previous methods., Comment: 8 pages, 8 figures, 2 tables, published in MNRAS

The origins of coalescing binary black holes (BBHs) detected by the advanced LIGO/Virgo are still under debate, and clues may be present in the joint mass-spin distribution of these merger events. Here we construct phenomenological models containing two sub-populations to investigate the BBH population detected in gravitational wave observations. We find that our models can explain the GWTC-3 data rather well, and several constraints to our model are required by the data: first, the maximum mass for the component with a stellar-origin, $m_{\rm max}$, is $39.1^{+2.4}_{-2.7}M_{\odot}$ at 90\% credibility; second, about $15\%$ of the mergers happen in dynamical environments, in which $7-16\%$ of events are hierarchical mergers, and these BHs have an average spin magnitude significantly larger than the first-generation mergers, with ${\rm d}\mu_{\rm a} > 0.4 $ at $99\%$ credibility; third, the dynamical component BHs tend to pair with each other with larger total mass and higher mass ratio. An independent analysis focusing on spins is also carried out, and we find that the spin amplitude of component BHs can be divided into two groups according to a division mass $m_{\rm d} = 46.1^{+5.6}_{-5.1}M_{\odot}$. These constraints can be naturally explained by current formation channels, and our results suggest that some of the observed events were likely from AGN disks., Comment: 21 pages, 8 figures. Published in The Astrophysical Journal Letters

The X-ray plateau emission observed in many Long gamma-ray bursts (LGRBs) has been usually interpreted as the spin-down luminosity of a rapidly spinning, highly magnetized neutron star (millisecond magnetar). If this is true, then the magnetar may emit extended gravitational wave (GW) emission associated with the X-ray plateau due to non-axisymmetric deformation or various stellar oscillations. The advanced LIGO and Virgo detectors have searched for long-duration GW transients for several years, no evidence of GWs from any magnetar has been found until now. In this work, we attempt to search for signature of GW radiation in the electromagnetic observation of 30 LGRBs under the assumption of the magnetar model. We utilize the observations of the LGRB plateau to constrain the properties of the new-born magnetar, including the initial spin period $P_0$, diploe magnetic field strength $B_p$ and the ellipticity $\epsilon$. We find that there are some tight relations between magnetar parameters, e.g., $\epsilon \propto B_p^{1.29}$ and $B_p\propto P_0^{1.14}$. In addition, we derive the GW strain for magnetar sample via their spin-down processes, and find that the GWs from these objects may not be detectable by the aLIGO and ET detectors. For a rapidly spinning magnetar ($P\sim 1\mbox{ ms}$, $ B \sim10^{15}\mbox{ G}$), the detection horizon for advanced LIGO O5 detector is $\sim 180\mbox{ Mpc}$. The detection of such GW signal associated with the X-ray plateau would be a smoking gun that the central engine of GRB is a magnetar., Comment: 15 pages, 7 figures,accepted for publication in ApJ

Coalescing binary black hole (BBH) systems are likely formed via several channels, and it is challenging to understand their formation / evolutionary processes. Some features in the mass function of the primary components ($m_1$), such as the distinct Gaussian-like peak located at $\sim 34M_\odot$, have been previously found. In this work, we investigate the possible dependence of the mass ratio ($q=m_2/m_1$) distribution on the primary mass. We find a Bayesian odds ratio of 18.1 in favor of divergence in the mass ratio distributions between the low- and high-mass ranges over an invariable mass ratio distribution. The BBHs with $m_1\gtrsim29M_{\odot}$ have a stronger preference to be symmetric compared to those with $m_1\lesssim29M_{\odot}$ at a 97.6\% credible level. Additionally, we find mild evidence that the BBHs with $m_1$ located in the Gaussian-like peak have a mass ratio distribution different from that of other BBHs. Our findings may be in favor of some formation channels, such as the chemically homogeneous evolution and the dynamical assembly in globular clusters/nuclear star clusters, which are more likely to provide symmetric BBHs in the high-mass range., Comment: 11 pages, 6 figures, 2 tables, published In ApJL

The DArk Matter Particle Explorer (DAMPE) is well suitable for searching for monochromatic and sharp $\gamma$-ray structures in the GeV$-$TeV range thanks to its unprecedented high energy resolution. In this work, we search for $\gamma$-ray line structures using five years of DAMPE data. To improve the sensitivity, we develop two types of dedicated data sets (including the BgoOnly data which is the first time to be used in the data analysis for the calorimeter-based gamma-ray observatories) and adopt the signal-to-noise ratio optimized regions of interest (ROIs) for different DM density profiles. No line signals or candidates are found between 10 and 300 GeV in the Galaxy. The constraints on the velocity-averaged cross section for $\chi\chi \to \gamma\gamma$ and the decay lifetime for $\chi \to \gamma\nu$, both at 95% confidence level, have been calculated and the systematic uncertainties have been taken into account. Comparing to the previous Fermi-LAT results, though DAMPE has an acceptance smaller by a factor of $\sim 10$, similar constraints on the DM parameters are achieved and below 100 GeV the lower limits on the decay lifetime are even stronger by a factor of a few. Our results demonstrate the potential of high-energy-resolution observations on dark matter detection., Comment: 14 pages, 8 figures. Update the content to keep up with the published version

The origins of the coalescing binary black holes (BBHs) detected by the advanced LIGO/Virgo are still in debate and clues may present in the mass and effective spin ($\chi_{\rm eff}$) distributions of these merger events. Here we analyze the GWTC-1 and GWTC-2.1 BBH events, and find evidence for a primary mass dependent $\chi_{\rm eff}$ distribution. Below (above) the mass of $\sim 50M_\odot$, a nearly zero (relatively high) effective spin is typical. In the joint analysis of the primary mass and $\chi_{\rm eff}$ distributions, we identify a sub-population of BBHs characterized by both a bump-like primary mass distribution extending to $\sim 100M_\odot$ and a relatively high $\chi_{\rm eff}\sim 0.3$, consistent with the expectation for the $\geq 2$g merger events in the disks of active galactic nuclei (AGNs). Our results suggest that the merger rate in AGN disks is around $1.7~{\rm Gpc}^{-3}{\rm yr}^{-1}$., Comment: 7 pages, 3 figures

GRB 060505 was the first well-known nearby (at redshift 0.089) "hybrid" gamma-ray burst that has a duration longer than 2 seconds but without the association of a supernova down to very stringent limits. The prompt $\gamma-$ray flash lasting $\sim 4$ sec could consist of an intrinsic short burst and its tail emission, but the sizable temporal lag ($\sim 0.35$ sec) as well as the environment properties led to the widely-accepted classification of a long duration gamma-ray burst originated from the collapse of a massive star. Here for the $ first$ time we report the convincing evidence for a thermal-like optical radiation component in the spectral energy distribution of the early afterglow emission. In comparison to AT2017gfo, the thermal radiation is $\sim 2$ times brighter and the temperature is comparable at similar epochs. The optical decline is much quicker than that in X-rays, which is also at odds with the fireball afterglow model but quite natural for the presence of a blue kilonova. Our finding reveals a neutron star merger origin of the hybrid GRB 060505 and strongly supports the theoretical speculation that some binary neutron stars can merge ultra-quickly (within $\sim 1$ Myr) after their formation when the surrounding region is still highly star-forming and the metallicity remains low. Gravitational wave and electromagnetic jointed observations are expected to confirm such scenarios in the near future., Comment: 20 pages, 7 figures, 2 tables

We perform a hierarchical Bayesian inference to investigate the population properties of the coalescing compact binaries involving at least one neutron star (NS). With the current observation data, we can not rule out either of the Double Gaussian, Single Gaussian and Uniform NS mass distribution models, although the mass distribution of the Galactic NSs is slightly preferred by the gravitational wave (GW) observations. The mass distribution of black holes (BHs) in the neutron star-black hole (NSBH) population is found to be similar to that for the Galactic X-ray binaries. Additionally, the ratio of the merger rate densities between NSBHs and BNSs is estimated to be about 3 : 7. The spin properties of the binaries, though constrained relatively poor, play nontrivial role in reconstructing the mass distribution of NSs and BHs. We find that a perfectly aligned spin distribution can be ruled out, while a purely isotropic distribution of spin orientation is still allowed.

In a binary system, the gravitational potential of the primary black hole may play an important role in enhancing the fallback accretion onto the lighter compact object newly formed in the second supernova explosion. As a result, the final masses of the binary compact objects would be correlated, as suggested recently by Safarzadeh \& Wysocki. In this work, we analyze the mass distribution of four gravitational-wave events, which are characterized by both a small mass ratio and a low mass ($\leq 5M_\odot$) of the light component, and find tentative evidence for a mass correlation among the objects. To evaluate the feasibility of testing such a hypothesis with upcoming observations, we carry out simulations with a mock population and perform Bayesian hierarchical inference for the mass distribution. We find that with dozens of low mass ratio events, whether there exists correlation in the component mass distributions or not can be robustly tested and the correlation, if it exists, can be well determined., Comment: 8 pages, 3 figures, ApJ published

We propose and demonstrate a high performance four-wavelength erbium-doped fiber laser (EDFL), enabled by a figure-8 compound-ring-cavity (F8-CRC) filter for single-longitudinal-mode (SLM) selection and a polarization-managed four-channel filter (PM-FCF) for defining four lasing wavelengths. We introduce a novel methodology utilizing signal-flow graph combined with Mason's rule to analyze a CRC filter in general and apply it to obtain the important design parameters for the F8-CRC filter used in this paper. By combining the functions of the F8-CRC filter and the PM-FCF filter assisted by the enhanced polarization hole-burning and polarization dependent loss, we achieve the EDFL with fifteen lasing states, including four single-, six dual-, four tri- and one quad-wavelength lasing operations. In particular, all the four single-wavelength operations are in stable SLM oscillation, typically with a linewidth of <600 Hz, a RIN of <=-154.58 dB/Hz@>=3 MHz and an output power fluctuation of <=+/-3.45%. In addition, all the six dual-wavelength operations have very similar performances, with the performance parameters close to those of the single-wavelength lasing operations. Finally, we achieve the wavelength spacing tuning of the dual-wavelength operations for the photonic generation of tunable microwave signals, and successfully obtain a signal at 23.10 GHz as a demonstration., Comment: 21 pages, 14 figures