scispace - formally typeset
Search or ask a question

Showing papers in "arXiv: High Energy Astrophysical Phenomena in 2021"


Posted ContentDOI
TL;DR: In this paper, the authors reported the detection of 1652 independent bursts with a peak burst rate of 122~hr^{-1}, in 59.5 hours spanning 47 days, and the burst energy distribution is bimodal, and well characterized by a combination of a log-normal function and a generalized Cauchy function.
Abstract: The event rate, energy distribution, and time-domain behaviour of repeating fast radio bursts (FRBs) contains essential information regarding their physical nature and central engine, which are as yet unknown. As the first precisely-localized source, FRB 121102 has been extensively observed and shows non-Poisson clustering of bursts over time and a power-law energy distribution. However, the extent of the energy distribution towards the fainter end was not known. Here we report the detection of 1652 independent bursts with a peak burst rate of 122~hr^{-1}, in 59.5 hours spanning 47 days. A peak in the isotropic equivalent energy distribution is found to be ~4.8 x 10^{37} erg at 1.25~GHz, below which the detection of bursts is suppressed. The burst energy distribution is bimodal, and well characterized by a combination of a log-normal function and a generalized Cauchy function. The large number of bursts in hour-long spans allow sensitive periodicity searches between 1 ms and 1000 s. The non-detection of any periodicity or quasi-periodicity poses challenges for models involving a single rotating compact object. The high burst rate also implies that FRBs must be generated with a high radiative efficiency, disfavoring emission mechanisms with large energy requirements or contrived triggering conditions.

91 citations


Journal ArticleDOI
TL;DR: In this article, the authors report the observation of gravitational waves from two compact binary coalescences in LIGO's and Virgo's third observing run with properties consistent with neutron star-black hole (NSBH) binaries.
Abstract: We report the observation of gravitational waves from two compact binary coalescences in LIGO's and Virgo's third observing run with properties consistent with neutron star-black hole (NSBH) binaries. The two events are named GW200105_162426 and GW200115_042309, abbreviated as GW200105 and GW200115; the first was observed by LIGO Livingston and Virgo, and the second by all three LIGO-Virgo detectors. The source of GW200105 has component masses $8.9^{+1.2}_{-1.5}\,M_\odot$ and $1.9^{+0.3}_{-0.2}\,M_\odot$, whereas the source of GW200115 has component masses $5.7^{+1.8}_{-2.1}\,M_\odot$ and $1.5^{+0.7}_{-0.3}\,M_\odot$ (all measurements quoted at the 90% credible level). The probability that the secondary's mass is below the maximal mass of a neutron star is 89%-96% and 87%-98%, respectively, for GW200105 and GW200115, with the ranges arising from different astrophysical assumptions. The source luminosity distances are $280^{+110}_{-110}$ Mpc and $300^{+150}_{-100}$ Mpc, respectively. The magnitude of the primary spin of GW200105 is less than 0.23 at the 90% credible level, and its orientation is unconstrained. For GW200115, the primary spin has a negative spin projection onto the orbital angular momentum at 88% probability. We are unable to constrain spin or tidal deformation of the secondary component for either event. We infer a NSBH merger rate density of $45^{+75}_{-33}\,\mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$ when assuming GW200105 and GW200115 are representative of the NSBH population, or $130^{+112}_{-69}\,\mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$ under the assumption of a broader distribution of component masses.

67 citations


Journal ArticleDOI
TL;DR: The ART-XC telescope as mentioned in this paper is a hard X-ray instrument with grazing incidence imaging optics on board the Spektr-Roentgen-Gamma (SRG) observatory, which was successfully launched into orbit around the second Lagrangian point (L2) with a Proton rocket from the Baikonur cosmodrome on 13 July 2019.
Abstract: ART-XC (Astronomical Roentgen Telescope - X-ray Concentrator) is the hard X-ray instrument with grazing incidence imaging optics on board the Spektr-Roentgen-Gamma (SRG) observatory. The SRG observatory is the flagship astrophysical mission of the Russian Federal Space Program, which was successively launched into orbit around the second Lagrangian point (L2) of the Earth-Sun system with a Proton rocket from the Baikonur cosmodrome on 13 July 2019. The ART-XC telescope will provide the first ever true imaging all-sky survey performed with grazing incidence optics in the 4-30 keV energy band and will obtain the deepest and sharpest map of the sky in the energy range of 4-12 keV. Observations performed during the early calibration and performance verification phase as well as during the on-going all-sky survey that started on 12 Dec. 2019 have demonstrated that the in-flight characteristics of the ART-XC telescope are very close to expectations based on the results of ground calibrations. Upon completion of its 4-year all-sky survey, ART-XC is expected to detect ~5000 sources (~3000 active galactic nuclei, including heavily obscured ones, several hundred clusters of galaxies, ~1000 cataclysmic variables and other Galactic sources), and to provide a high-quality map of the Galactic background emission in the 4-12 keV energy band. ART-XC is also well suited for discovering transient X-ray sources. In this paper, we describe the telescope, results of its ground calibrations, major aspects of the mission, the in-flight performance of ART-XC and first scientific results.

62 citations


Posted Content
TL;DR: In this article, an equatorial plasmoid-unstable current sheet formed in a transient, non-axisymmetric, low-density magnetosphere within the inner few Schwarzschild radii.
Abstract: Magnetic reconnection can power bright and rapid flares originating from the inner magnetosphere of accreting black holes. We conduct extremely high resolution ($5376\times2304\times2304$ cells) general-relativistic magnetohydrodynamics simulations, capturing plasmoid-mediated reconnection in a 3D magnetically arrested disk for the first time. We show that an equatorial, plasmoid-unstable current sheet forms in a transient, non-axisymmetric, low-density magnetosphere within the inner few Schwarzschild radii. Magnetic flux bundles escape from the event horizon through reconnection at the universal plasmoid-mediated rate in this current sheet. The reconnection feeds on the highly-magnetized plasma in the jets and heats the plasma that ends up trapped in flux bundles to temperatures proportional to the jet's magnetization. The escaped flux bundles can complete a full orbit as low-density hot spots, consistent with Sgr A$^{*}$ observations by the GRAVITY interferometer. Reconnection near the horizon produces sufficiently energetic plasma to explain flares from accreting black holes, such as the TeV emission observed from M87. The drop in mass accretion rate during the flare, and the resulting low-density magnetosphere make it easier for very high energy photons produced by reconnection-accelerated particles to escape. The extreme resolution results in a converged plasmoid-mediated reconnection rate that directly determines the timescales and properties of the flare.

58 citations


Journal ArticleDOI
TL;DR: In this paper, general relativistic magnetohydrodynamics (MHD) simulations have greatly improved our understanding of jet production and its dependence on the black hole spin, showing that both disks and jets are sensitive to the magnetic flux threading the accretion flow as well as possible misalignment between the angular momentum of the fluid and the spin.
Abstract: There is a broad consensus that accretion onto supermassive black holes and consequent jet formation power the observed emission from active galactic nuclei (AGNs). However, there has been less agreement about how jets form in accretion flows, their possible relationship to black hole spin, and how they interact with the surrounding medium. There have also been theoretical concerns about instabilities in standard accretion disk models and lingering discrepancies with observational constraints. Despite seemingly successful applications to X-ray binaries, the standard accretion disk model faces a growing list of observational constraints that challenge its application to AGNs. Theoretical exploration of these questions has become increasingly reliant on numerical simulations owing to the dynamic nature of these flows and the complex interplay between hydrodynamics, magnetic fields, radiation transfer, and curved spacetime. We conclude the following: The advent of general relativistic magnetohydrodynamics (MHD) simulations has greatly improved our understanding of jet production and its dependence on black hole spin. Simulation results show both disks and jets are sensitive to the magnetic flux threading the accretion flow as well as possible misalignment between the angular momentum of the accretion flow and the black hole spin. Radiation MHD simulations are providing new insights into the stability of luminous accretion flows and highlighting the potential importance of radiation viscosity, UV opacity from atoms, and spiral density waves in AGNs.

57 citations


Journal ArticleDOI
TL;DR: In this paper, the relativistic Shapiro time delay of the radio pulsar J0740+6620 has been improved to 1.14+0.17 + 0.15 kpc by using a combination of the 100m Green Bank Telescope and the Canadian Hydrogen Intensity Mapping Experiment telescope.
Abstract: We report results from continued timing observations of PSR J0740+6620, a high-mass, 2.8-ms radio pulsar in orbit with a likely ultra-cool white dwarf companion. Our data set consists of combined pulse arrival-time measurements made with the 100-m Green Bank Telescope and the Canadian Hydrogen Intensity Mapping Experiment telescope. We explore the significance of timing-based phenomena arising from general-relativistic dynamics and variations in pulse dispersion. When using various statistical methods, we find that combining $\sim 1.5$ years of additional, high-cadence timing data with previous measurements confirms and improves upon previous estimates of relativistic effects within the PSR J0740+6620 system, with the pulsar mass $m_{\rm p} = 2.08^{+0.07}_{-0.07}$ M$_\odot$ (68.3\% credibility) determined by the relativistic Shapiro time delay. For the first time, we measure secular variation in the orbital period and argue that this effect arises from apparent acceleration due to significant transverse motion. After incorporating contributions from Galactic differential rotation and off-plane acceleration in the Galactic potential, we obtain a model-dependent distance of $d = 1.14^{+0.17}_{-0.15}$ kpc (68.3\% credibility). This improved distance confirms the ultra-cool nature of the white dwarf companion determined from recent optical observations. We discuss the prospects for future observations with next-generation facilities, which will likely improve the precision on $m_{\rm p}$ for J0740+6620 by an order of magnitude within the next few years.

54 citations


Journal ArticleDOI
TL;DR: Quasi-Periodic Eruptions (QPEs) are extreme high-amplitude bursts of X-ray radiation recurring every few hours and originating near the central supermassive black holes in galactic nuclei as discussed by the authors.
Abstract: Quasi-Periodic Eruptions (QPEs) are extreme high-amplitude bursts of X-ray radiation recurring every few hours and originating near the central supermassive black holes in galactic nuclei. It is currently unknown what triggers these events, how long they last and how they are connected to the physical properties of the inner accretion flows. Previously, only two such sources were known, found either serendipitously or in archival data, with emission lines in their optical spectra classifying their nuclei as hosting an actively accreting supermassive black hole. Here we present the detection of QPEs in two further galaxies, obtained with a blind and systematic search over half of the X-ray sky. The optical spectra of these galaxies show no signature of black hole activity, indicating that a pre-existing accretion flow typical of active nuclei is not required to trigger these events. Indeed, the periods, amplitudes and profiles of the newly discovered QPEs are inconsistent with current models that invoke radiation-pressure driven accretion disk instabilities. Instead, QPEs might be driven by an orbiting compact object. Furthermore, their observed properties require the mass of the secondary object to be much smaller than the main body and future X-ray observations may constrain possible changes in the period due to orbital evolution. This scenario could make QPEs a viable candidate for the electromagnetic counterparts of the so-called extreme mass ratio inspirals, with considerable implications for multi-messenger astrophysics and cosmology.

51 citations


Posted Content
TL;DR: In this article, it was shown that the muon discrepancy starts at the TeV scale, which suggests that this change in hadron composition is observable at the Large Hadron Collider.
Abstract: High-energy cosmic rays are observed indirectly by detecting the extensive air showers initiated in Earth's atmosphere. The interpretation of these observations relies on accurate models of air shower physics, which is a challenge and an opportunity to test QCD under extreme conditions. Air showers are hadronic cascades, which eventually decay into muons. The muon number is a key observable to infer the mass composition of cosmic rays. Air shower simulations with state-of-the-art QCD models show a significant muon deficit with respect to measurements; this is called the Muon Puzzle. The origin of this discrepancy has been traced to the composition of secondary particles in hadronic interactions. The muon discrepancy starts at the TeV scale, which suggests that this change in hadron composition is observable at the Large Hadron Collider. An effect that can potentially explain the puzzle has been discovered at the LHC, but needs to be confirmed for forward produced hadrons with LHCb, and with future data on oxygen beams.

48 citations


Journal ArticleDOI
TL;DR: The GPPS survey as discussed by the authors was designed to discover pulsars within the Galactic plane from the FAST L-band 19-beam receiver, and the results of the survey have been published.
Abstract: Discovery of pulsars is one of the main goals for large radio telescopes. The Five-hundred-meter Aperture Spherical radio Telescope (FAST), that incorporates an L-band 19-beam receiver with a system temperature of about 20~K, is the most sensitive radio telescope utilized for discovering pulsars. We designed the {\it snapshot} observation mode for a FAST key science project, the Galactic Plane Pulsar Snapshot (GPPS) survey, in which every four nearby pointings can observe {\it a cover} of a sky patch of 0.1575 square degrees through beam-switching of the L-band 19-beam receiver. The integration time for each pointing is 300 seconds so that the GPPS observations for a cover can be made in 21 minutes. The goal of the GPPS survey is to discover pulsars within the Galactic latitude of $\pm10^{\circ}$ from the Galactic plane, and the highest priority is given to the inner Galaxy within $\pm5^{\circ}$. Up to now, the GPPS survey has discovered 201 pulsars, including currently the faintest pulsars which cannot be detected by other telescopes, pulsars with extremely high dispersion measures (DMs) which challenge the currently widely used models for the Galactic electron density distribution, pulsars coincident with supernova remnants, 40 millisecond pulsars, 16 binary pulsars, some nulling and mode-changing pulsars and rotating radio transients (RRATs). The follow-up observations for confirmation of new pulsars have polarization-signals recorded for polarization profiles of the pulsars. Re-detection of previously known pulsars in the survey data also leads to significant improvements in parameters for 64 pulsars. The GPPS survey discoveries are published and will be updated at this http URL .

48 citations


RepositoryDOI
R. Alves Batista, Amin, G. Barenboim, Nicola Bartolo, Daniel Baumann, A. Bauswein, Emilio Bellini, D. Benisty, Gianfranco Bertone, P. Blasi, C.G. Böhmer, Ž Bošnjak, T. Bringmann, Clare Burrage, M. Bustamante, J. Calderón Bustillo, C.T. Byrnes, F. Calore, R. Catena, D.G. Cerdeño, S.S. Cerri, Marco Chianese, Katy Clough, Alex Cole, P. Coloma, Adam Coogan, L. Covi, D. Cutting, A.C. Davis, C. De Rham, A. Di Matteo, G. Domènech, M. Drewes, T. Dietrich, T.D.P. Edwards, I. Esteban, R. Erdem, C. Evoli, Matteo Fasiello, S.M. Feeney, R.Z. Ferreira, A. Fialkov, N. Fornengo, S. Gabici, T. Galatyuk, Daniele Gaggero, D. Grasso, C. Guépin, J. Harz, M. Herrero-Valea, T. Hinderer, N.B. Hogg, D.C. Hooper, D. Iocco, J. Isern, Konstantin Karchev, Bradley J. Kavanagh, M. Korsmeier, K. Kotera, Kazuya Koyama, B. Krishnan, Julien Lesgourgues, J. Levi Said, Lucas Lombriser, C.S. Lorenz, S. Manconi, M. Mapelli, A. Marcowith, S.B. Markoff, D.J. Marsh, M. Martinelli, C.J.A.P. Martinsolami, P. Millington, P. Moesta, K. Nippel, V. Niro, E. O'Connor, F. Oikonomou, C.F. Paganini, G. Pagliaroli, Paolo Pani, C. Pfrommer, Silvia Pascoli, L. Pinol, L. Pizzuti, R.A. Porto, A. Pound, F. Quevedo, G.G. Raffelt, Alvise Raccanelli, E. Ramirez-Ruiz, M. Raveri, S. Renaux-Petel, Angelo Ricciardone, A. Rida Khalifeh, Antonio Riotto, R. Roiban, J. Rubio, M. Sahlén, Nashwan Sabti, L. Sagunski, N. Šarčević, K. Schmitz, Pedro Schwaller, T. Schwetz, A. Sedrakian, E. Sellentin, A. Serenelli, P.D. Serpico, E.I. Sfakianakis, S. Shalgar, A. Silvestri, I. Tamborra, Konstantinos Tanidis, D. Teresi, A.A. Tokareva, L. Tolos, S. Trojanowski, R. Trotta, C. Uhlemann, F.R. Urban, Filippo Vernizzi, A. van Vliet, F.L. Villante, A. Vincent, J. Vink, E. Vitagliano, Christophe Weniger, A. Wickenbrock, W. Winter, S. Zell, M. Zeng 
TL;DR: The European Consortium for Astroparticle Theory (EuCAPT) white paper as mentioned in this paper explores upcoming theoretical opportunities and challenges for our field of research with particular emphasis on the possible synergies among different subfields, and the prospects for solving the most fundamental open questions with multi-messenger observations.
Abstract: Astroparticle physics is undergoing a profound transformation, due to a series of extraordinary new results, such as the discovery of high-energy cosmic neutrinos with IceCube, the direct detection of gravitational waves with LIGO and Virgo, and many others. This white paper is the result of a collaborative effort that involved hundreds of theoretical astroparticle physicists and cosmologists, under the coordination of the European Consortium for Astroparticle Theory (EuCAPT). Addressed to the whole astroparticle physics community, it explores upcoming theoretical opportunities and challenges for our field of research, with particular emphasis on the possible synergies among different subfields, and the prospects for solving the most fundamental open questions with multi-messenger observations.

46 citations


Journal ArticleDOI
TL;DR: In this article, the authors used spectrophotometry of PSR J1810+1744 to show a flat, but asymmetric light-curve maximum and a deep, narrow minimum, which indicates strong gravity darkening near the L_1 point, along with a heated pole and surface winds.
Abstract: Keck-telescope spectrophotometry of the companion of PSR J1810+1744 shows a flat, but asymmetric light-curve maximum and a deep, narrow minimum. The maximum indicates strong gravity darkening near the L_1 point, along with a heated pole and surface winds. The minimum indicates a low underlying temperature and substantial limb darkening. The gravity darkening is a consequence of extreme pulsar heating and the near-filling of the Roche lobe. Light-curve modeling gives a binary inclination i=65.7+/-0.4deg. With the Keck-measured radial-velocity amplitude K_c=462.3+/-2.2km/s, this gives an accurate neutron star mass M_NS=2.13+/-0.04M_o, with important implications for the dense-matter equation of state. A classic direct-heating model, ignoring the L_1 gravitational darkening, would predict an unphysical M_NS>3M_o. A few other ``spider" pulsar binaries have similar large heating and fill factor; thus, they should be checked for such effects.

Posted Content
TL;DR: In this article, the authors presented a catalog of 536 fast radio bursts (FRBs) detected by the Canadian Hydrogen Intensity Mapping Experiment Fast Radio Burst (CHIME/FRB) Project between 400 and 800 MHz from 2018 July 25 to 2019 July 1, including 62 bursts from 18 previously reported repeating sources.
Abstract: We present a catalog of 536 fast radio bursts (FRBs) detected by the Canadian Hydrogen Intensity Mapping Experiment Fast Radio Burst (CHIME/FRB) Project between 400 and 800 MHz from 2018 July 25 to 2019 July 1, including 62 bursts from 18 previously reported repeating sources. The catalog represents the first large sample, including bursts from repeaters and non-repeaters, observed in a single survey with uniform selection effects. This facilitates comparative and absolute studies of the FRB population. We show that repeaters and apparent non-repeaters have sky locations and dispersion measures (DMs) that are consistent with being drawn from the same distribution. However, bursts from repeating sources differ from apparent non-repeaters in intrinsic temporal width and spectral bandwidth. Through injection of simulated events into our detection pipeline, we perform an absolute calibration of selection effects to account for systematic biases. We find evidence for a population of FRBs - comprising a large fraction of the overall population - with a scattering time at 600 MHz in excess of 10 ms, of which only a small fraction are observed by CHIME/FRB. We infer a power-law index for the cumulative fluence distribution of $\alpha=-1.40\pm0.11(\textrm{stat.})^{+0.06}_{-0.09}(\textrm{sys.})$, consistent with the $-3/2$ expectation for a non-evolving population in Euclidean space. We find $\alpha$ is steeper for high-DM events and shallower for low-DM events, which is what would be expected when DM is correlated with distance. We infer a sky rate of $[820\pm60(\textrm{stat.})^{+220}_{-200}({\textrm{sys.}})]/\textrm{sky}/\textrm{day}$ above a fluence of 5 Jy ms at 600 MHz, with scattering time at $600$ MHz under 10 ms, and DM above 100 pc cm$^{-3}$.

Posted ContentDOI
TL;DR: In this paper, the authors discuss the impact of laboratory experiments and astrophysical observation of neutron stars (NSs) on its equation of state (EoS): (a) The new measurement of neutron skin thickness of $\rm ^{208} \! Pb$, $R_{\rm skin}^{208} = 0.29 \pm 0.07$ fm by the PREX-II experiment.
Abstract: In this paper, we discuss the impact of the following laboratory experiments and astrophysical observation of neutron stars (NSs) on its equation of state (EoS): (a) The new measurement of neutron skin thickness of $\rm ^{208} \! Pb$, $R_{\rm skin}^{208} = 0.29 \pm 0.07$ fm by the PREX-II experiment. (b) The mass measurement of PSR J0740+6620 has been slightly revised down by including additional $\sim 1.5$ years of pulsar timing data. As well as the radius measurement of PSR J0740+6620 by joint NICER/XMM-Newton collaboration which has a similar size to PSR J0030+0451. We combine these information using Bayesian statistics along with the previous LIGO/Virgo and NICER observations of NS using a hybrid nuclear+piecewise polytrope (PP) EoS parameterization. Our findings are as follows: (a). Adding PREX-II result yields the value of empirical parameter $L = 69^{+21}_{-19}$ MeV, $R_{\rm skin}^{208} = 0.20_{-0.05}^{+0.05}$ fm, and radius of a $1.4 M_{\odot}$ ($R_{1.4}) = 12.75_{-0.54}^{+ 0.42}$ km at $1 \sigma$ confidence interval (CI). We find these inferred values are mostly dominated by the combined astrophysical observations as the measurement uncertainty in $R_{\rm skin}^{208}$ by PREX-II is much broader. Also, a better measurement of $R_{\rm skin}^{208}$ might have a little effect on the radius of low mass NSs, but for the high masses there will be almost no effect. (b) After adding The revised mass and radius measurement of PSR J0740+6620, we find the inferred radius of NSs is slightly pushed towards the larger values and the uncertainty on the radius of a $2.08 M_{\odot}$ NS is moderately improved.

Journal ArticleDOI
TL;DR: In this article, the authors proposed a new scenario in which fast radio bursts are powered by short-lived relativistic outflows from accreting black holes or neutron stars, which propagate into the cavity of the pre-existing (quiescent'') jet.
Abstract: The discovery of periodicity in the arrival times of the fast radio bursts (FRBs) poses a challenge to the oft-studied magnetar scenarios. However, models that postulate that FRBs result from magnetized shocks or magnetic reconnection in a relativistic outflow are not specific to magnetar engines; instead, they require only the impulsive injection of relativistic energy into a dense magnetized medium. Motivated thus, we outline a new scenario in which FRBs are powered by short-lived relativistic outflows (``flares'') from accreting black holes or neutron stars, which propagate into the cavity of the pre-existing (``quiescent'') jet. In order to reproduce FRB luminosities and rates, we are driven to consider binaries of stellar-mass compact objects undergoing super-Eddington mass-transfer, similar to ultraluminous X-ray (ULX) sources. Indeed, the host galaxies of FRBs, and their spatial offsets within their hosts, show broad similarities with ULXs. Periodicity on timescales of days to years could be attributed to precession (e.g., Lens-Thirring) of the polar accretion funnel, along which the FRB emission is geometrically and relativistically beamed, which sweeps across the observer line of sight. Accounting for the most luminous FRBs via accretion power may require a population of binaries undergoing brief-lived phases of unstable (dynamical-timescale) mass-transfer. This will lead to secular evolution in the properties of some repeating FRBs on timescales of months to years, followed by a transient optical/IR counterpart akin to a luminous red nova, or a more luminous accretion-powered optical/X-ray transient. We encourage targeted FRB searches of known ULX sources.

Journal ArticleDOI
TL;DR: In this article, the gamma-ray emission from the supernova remnant G106.3+2.7 above 10 teraelectronvolts (TeV) has been observed and it was shown that the very high-energy gamma ray emission above 10 TeV is well correlated with a molecular cloud rather than the pulsar PSR J2229+6114.
Abstract: Cosmic rays (protons and other atomic nuclei) are believed to gain energies of petaelectronvolts (PeV) and beyond at astrophysical particle accelerators called 'PeVatrons' inside our Galaxy. Although a characteristic feature of a PeVatron is expected to be a hard gamma-ray energy spectrum that extends beyond 100 teraelectronvolts (TeV) without a cutoff, none of the currently known sources exhibits such a spectrum due to the low maximum energy of accelerated cosmic rays or insufficient detector sensitivity around 100 TeV. Here we report the observation of gamma-ray emission from the supernova remnant G106.3+2.7 above 10 TeV. This work provides flux data points up to and above 100 TeV and indicates that the very-high-energy gamma-ray emission above 10 TeV is well correlated with a molecular cloud rather than the pulsar PSR J2229+6114. Regarding the gamma-ray emission mechanism of G106.3+2.7, this morphological feature appears to favor a hadronic origin via the {\pi}0 decay caused by accelerated relativistic protons over a leptonic one via the inverse-Compton scattering by relativistic electrons. Furthermore, we point out that an X-ray flux upper limit on the synchrotron spectrum would provide important information to firmly establish the hadronic scenario as the mechanism of particle acceleration at the source.

Posted Content
TL;DR: In this paper, a large ensemble of randomly generated model-independent equations of state (EoSs) and the corresponding rotating stellar structures without the use of quasi-universal relations is analyzed.
Abstract: Recent rapid progress in multimessenger observations of neutron stars (NSs) offers great potential to constrain the properties of strongly interacting matter under the most extreme conditions. In order to fully exploit the current observational inputs and to study the impact of future observations, we analyze a large ensemble of randomly generated model-independent equations of state (EoSs) and the corresponding rotating stellar structures without the use of quasi-universal relations. We discuss the compatibility and the impact on the EoS of various hypotheses and measurements, including those involving the merger product in GW170817, the binary merger components in GW190814, and radius measurements of PSR J0740+6620. We obtain an upper limit for the dimensionless spin of a rigidly rotating NS, |chi| < 0.84, and find that the conservative hypothesis that the remnant in GW170817 ultimately collapsed to a black hole strongly constrains the EoS and the maximal mass of NSs, implying M_TOV < 2.53M_sol (or M_TOV < 2.19M_sol if we assume that a hypermassive NS was created). Furthermore, we find that radius measurements of massive NSs stars can dramatically limit the EoS, with strict lower limits around 12 km and above offering particularly efficient constraints. Finally, we find that radius values for a two-solar mass NS around 11-13 km would be completely compatible with the presence of quark matter inside massive NSs.

Journal ArticleDOI
TL;DR: The Cherenkov Telescope Array (CTA) provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics as mentioned in this paper, which can be explored as part of the Key Science Projects of CTA through simulated observations of active galactic nuclei and their relativistic jets.
Abstract: The Cherenkov Telescope Array (CTA), the new-generation ground-based observatory for gamma-ray astronomy, provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics. We study some of the salient areas of gamma-ray cosmology that can be explored as part of the Key Science Projects of CTA, through simulated observations of active galactic nuclei (AGN) and of their relativistic jets. Observations of AGN with CTA will enable a measurement of gamma-ray absorption on the extragalactic background light with a statistical uncertainty below 15% up to a redshift z = 2 and to constrain or detect gamma-ray halos up to intergalactic-magnetic-field strengths of at least 0.3 pG. Extragalactic observations with CTA also show promising potential to probe physics beyond the Standard Model. The best limits on Lorentz invariance violation from gamma-ray astronomy will be improved by a factor of at least two to three. CTA will also probe the parameter space in which axion-like particles could constitute a significant fraction, if not all, of dark matter. We conclude on the synergies between CTA and other upcoming facilities that will foster the growth of gamma-ray cosmology.

Book ChapterDOI
TL;DR: In this article, the main formation channels of binary black hole (BBH) from stellar evolution and dynamics are discussed, and the importance of reconstructing the cosmic evolution of BBHs is discussed.
Abstract: These are exciting times for binary black hole (BBH) research. LIGO and Virgo detections are progressively drawing a spectacular fresco of BBH masses, spins and merger rates. In this review, we discuss the main formation channels of BBHs from stellar evolution and dynamics. Uncertainties on massive star evolution (e.g., stellar winds, rotation, overshooting and nuclear reaction rates), core-collapse supernovae and pair instability still hamper our comprehension of the mass spectrum and spin distribution of black holes (BHs), but substantial progress has been done in the field over the last few years. On top of this, the efficiency of mass transfer in a binary system and the physics of common envelope substantially affect the final BBH demography. Dynamical processes in dense stellar systems can trigger the formation of BHs in the mass gap and intermediate-mass BHs via hierarchical BH mergers and via multiple stellar collisions. Finally, we discuss the importance of reconstructing the cosmic evolution of BBHs.

Journal ArticleDOI
TL;DR: In this paper, the authors analyzed 5 epochs of NICER data of the black hole X-ray binary MAXI J1820+070 during the bright hard-to-soft state transition in its 2018 outburst with both reflection spectroscopy and Fourier-resolved timing analysis.
Abstract: We analyze 5 epochs of NICER data of the black hole X-ray binary MAXI J1820+070 during the bright hard-to-soft state transition in its 2018 outburst with both reflection spectroscopy and Fourier-resolved timing analysis. We confirm the previous discovery of reverberation lags in the hard state, and find that the frequency range where the (soft) reverberation lag dominates decreases with the reverberation lag amplitude increasing during the transition, suggesting an increasing X-ray emitting region, possibly due to an expanding corona. By jointly fitting the lag-energy spectra in a number of broad frequency ranges with the reverberation model reltrans, we find the increase in reverberation lag is best described by an increase in the X-ray coronal height. This result, along with the finding that the corona contracts in the hard state, suggests a close relationship between spatial extent of the X-ray corona and the radio jet. We find the corona expansion (as probed by reverberation) precedes a radio flare by ~5 days, which may suggest that the hard-to-soft transition is marked by the corona expanding vertically and launching a jet knot that propagates along the jet stream at relativistic velocities.

Journal ArticleDOI
TL;DR: In this article, the authors estimate the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740$+$6620, conditional on pulse profile modeling of Neutron Star Interior Composition Explorer X-ray Timing Instrument (NICER XTI) event data.
Abstract: We report on Bayesian estimation of the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740$+$6620, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray Timing Instrument (NICER XTI) event data. We condition on informative pulsar mass, distance, and orbital inclination priors derived from the joint NANOGrav and CHIME/Pulsar wideband radio timing measurements of arXiv:2104.00880. We use XMM European Photon Imaging Camera spectroscopic event data to inform our X-ray likelihood function. The prior support of the pulsar radius is truncated at 16 km to ensure coverage of current dense matter models. We assume conservative priors on instrument calibration uncertainty. We constrain the equatorial radius and mass of PSR J0740$+$6620 to be $12.39_{-0.98}^{+1.30}$ km and $2.072_{-0.066}^{+0.067}$ M$_{\odot}$ respectively, each reported as the posterior credible interval bounded by the 16% and 84% quantiles, conditional on surface hot regions that are non-overlapping spherical caps of fully-ionized hydrogen atmosphere with uniform effective temperature; a posteriori, the temperature is $\log_{10}(T$ [K]$)=5.99_{-0.06}^{+0.05}$ for each hot region. All software for the X-ray modeling framework is open-source and all data, model, and sample information is publicly available, including analysis notebooks and model modules in the Python language. Our marginal likelihood function of mass and equatorial radius is proportional to the marginal joint posterior density of those parameters (within the prior support) and can thus be computed from the posterior samples.

Posted Content
TL;DR: In this paper, the authors reported a radius measurement based on fits of rotating hot spot patterns to Neutron Star Interior Composition Explorer (NICER) and X-ray Multi-Mirror (XMM-Newton) observations.
Abstract: PSR J0740$+$6620 has a gravitational mass of $2.08\pm 0.07~M_\odot$, which is the highest reliably determined mass of any neutron star. As a result, a measurement of its radius will provide unique insight into the properties of neutron star core matter at high densities. Here we report a radius measurement based on fits of rotating hot spot patterns to Neutron Star Interior Composition Explorer (NICER) and X-ray Multi-Mirror (XMM-Newton) X-ray observations. We find that the equatorial circumferential radius of PSR J0740$+$6620 is $13.7^{+2.6}_{-1.5}$ km (68%). We apply our measurement, combined with the previous NICER mass and radius measurement of PSR J0030$+$0451, the masses of two other $\sim 2~M_\odot$ pulsars, and the tidal deformability constraints from two gravitational wave events, to three different frameworks for equation of state modeling, and find consistent results at $\sim 1.5-3$ times nuclear saturation density. For a given framework, when all measurements are included the radius of a $1.4~M_\odot$ neutron star is known to $\pm 4$% (68% credibility) and the radius of a $2.08~M_\odot$ neutron star is known to $\pm 5$%. The full radius range that spans the $\pm 1\sigma$ credible intervals of all the radius estimates in the three frameworks is $12.45\pm 0.65$ km for a $1.4~M_\odot$ neutron star and $12.35\pm 0.75$ km for a $2.08~M_\odot$ neutron star.

Posted Content
TL;DR: In this paper, the authors examine the possible future evolution of Melnick 34, the most massive known binary star system, and show that various assumptions regarding stellar and binary physics predict a wide variety of outcomes: from a close BH-BH binary (which would lead to a potentially detectable coalescence), through a wide BHbH binary, or a Thorne-Zytkow object, to a complete disruption of both objects by pair-instability supernovae.
Abstract: The LIGO/Virgo gravitational--wave observatories have detected 50 BH-BH coalescences. This sample is large enough to have allowed several recent studies to draw conclusions about the branching ratios between isolated binaries versus dense stellar clusters as the origin of double BHs. It has also led to the exciting suggestion that the population is highly likely to contain primordial black holes. Here we demonstrate that such conclusions cannot yet be robust, because of the large current uncertainties in several key aspects of binary stellar evolution. These include the development and survival of a common envelope, the mass and angular momentum loss during binary interactions, mixing in stellar interiors, pair-instability mass loss and supernova outbursts. Using standard tools such as the population synthesis codes StarTrack and COMPAS and the detailed stellar evolution code MESA, we examine as a case study the possible future evolution of Melnick 34, the most massive known binary star system. We show that, despite its well-known orbital architecture, various assumptions regarding stellar and binary physics predict a wide variety of outcomes: from a close BH-BH binary (which would lead to a potentially detectable coalescence), through a wide BH-BH binary (which might be seen in microlensing observations), or a Thorne-Zytkow object, to a complete disruption of both objects by pair-instability supernovae. Thus since the future of massive binaries is inherently uncertain, sound predictions about the properties of BH-BH systems are highly challenging at this time. Consequently, drawing conclusions about the formation channels for the LIGO/Virgo BH-BH merger population is premature.

Posted Content
TL;DR: In this paper, the authors employ large-scale 3D particle-in-cell simulations of reconnection in magnetically-dominated pair plasmas to study the energization physics of high-energy particles.
Abstract: Magnetic reconnection is invoked as one of the primary mechanisms to produce energetic particles. We employ large-scale three-dimensional (3D) particle-in-cell simulations of reconnection in magnetically-dominated ($\sigma=10$) pair plasmas to study the energization physics of high-energy particles. We identify a novel acceleration mechanism that only operates in 3D. For weak guide fields, 3D plasmoids / flux ropes extend along the $z$ direction of the electric current for a length comparable to their cross-sectional radius. Unlike in 2D simulations, where particles are buried in plasmoids, in 3D we find that a fraction of particles with $\gamma\gtrsim 3\sigma$ can escape from plasmoids by moving along $z$, and so they can experience the large-scale fields in the upstream region. These "free" particles preferentially move in $z$ along Speiser-like orbits sampling both sides of the layer, and are accelerated linearly in time -- their Lorentz factor scales as $\gamma\propto t$, in contrast to $\gamma\propto \sqrt{t}$ in 2D. The energy gain rate approaches $\sim eE_{\rm rec}c$, where $E_{\rm rec}\simeq 0.1 B_0$ is the reconnection electric field and $B_0$ the upstream magnetic field. The spectrum of free particles is hard, $dN_{\rm free}/d\gamma\propto \gamma^{-1.5}$, contains $\sim 20\%$ of the dissipated magnetic energy independently of domain size, and extends up to a cutoff energy scaling linearly with box size. Our results demonstrate that relativistic reconnection in GRB and AGN jets may be a promising mechanism for generating ultra-high-energy cosmic rays.

Journal ArticleDOI
TL;DR: Observations of the γ-ray burst GRB 200415A have a rapid onset, very fast time variability, flat spectra and substantial sub-millisecond spectral evolution, which match well with those expected for a giant flare from an extragalactic magnetar.
Abstract: Magnetars are slowly-rotating neutron stars with extremely strong magnetic fields ($10^{13-15}$ G), episodically emitting $\sim100$ ms long X-ray bursts with energies of $\sim10^{40-41}$ erg. Rarely, they produce extremely bright, energetic giant flares that begin with a short ($\sim0.2$ s), intense flash, followed by fainter, longer lasting emission modulated by the magnetar spin period (typically 2-12 s), thus confirming their origin. Over the last 40 years, only three such flares have been observed in our local group; they all suffered from instrumental saturation due to their extreme intensity. It has been proposed that extra-galactic giant flares likely constitute a subset of short gamma-ray bursts, noting that the sensitivity of current instrumentation prevents us from detecting the pulsating tail, while the initial bright flash is readily observable out to distances $\sim 10-20$ Mpc. Here, we report X- and gamma-ray observations of GRB 200415A, which exhibits a rapid onset, very fast time variability, flat spectra and significant sub-millisecond spectral evolution. These attributes match well with those expected for a giant flare from an extra-galactic magnetar, noting that GRB 200415A is directionally associated with the galaxy NGC 253 ($\sim$3.5 Mpc away). The detection of $\sim3$ MeV photons provides definitive evidence for relativistic motion of the emitting plasma. The observed rapid spectral evolution can naturally be generated by radiation emanating from such rapidly-moving gas in a rotating magnetar.

Journal ArticleDOI
TL;DR: In this article, the authors studied X-ray spectra from the outburst rise of the accreting black-hole binary MAXI J1820+070 and found that having the disk inclinations within those of either the binary or the jet imply significant changes of the disk inner radius during the luminous part of the hard spectral state, with that radius changing from $>$100 to $sim$10 gravitational radii.
Abstract: We study X-ray spectra from the outburst rise of the accreting black-hole binary MAXI J1820+070. We find that models having the disk inclinations within those of either the binary or the jet imply significant changes of the accretion disk inner radius during the luminous part of the hard spectral state, with that radius changing from $>$100 to $\sim$10 gravitational radii. The main trend is a decrease with the decreasing spectral hardness. Our analysis requires the accretion flow to be structured, with at least two components with different spectral slopes. The harder component dominates the bolometric luminosity and produces strong, narrow, X-ray reflection features. The softer component is responsible for the underlying broader reflection features. The data are compatible with the harder component having a large scale height, located downstream the disk truncation radius, and reflecting mostly from remote parts of the disk. The softer component forms a corona above the disk up to some transition radius. Our findings can explain the changes of the characteristic variability time scales, found in other works, as being driven by the changes of the disk characteristic radii.

Journal ArticleDOI
TL;DR: In this paper, the effects of a dark matter spike on the neighborhood of Sgr A*, the black hole (BH) in the center of the Milky Way galaxy were investigated.
Abstract: In this paper, we investigate the effects of a dark matter (DM) spike on the neighborhood of Sgr A*, the black hole (BH) in the center of the Milky Way galaxy. Our main goal is to investigate whether current and future astronomical observations of Sgr A* could detect the presence of such a DM spike. At first, we construct the spacetime metric around a static and spherically symmetric BH with a DM spike, and later this solution is generalized for a rotating BH using the Newman-Janis-Azreg Ainou algorithm. For the static BH metric, we use the data of the S2 star orbiting the Sgr A* to determine and analyze the constraints on the two free parameters characterizing the density and the innermost boundary of the DM halo surrounding the BH. Furthermore, by making use of the available observational data for the DM spike density $\rho_\text{sp}$ and the DM spike radius $R_\text{sp}$ in the Milky Way galaxy, we consider a geometrically-thick accretion disk model around the Sgr A* BH and demonstrate that the effect of DM distribution on the shadow radius and the image of the BH is considerably weak for realistic DM densities, becoming significant only when the DM density is of the order $\rho_\text{sp} \sim (10^{-19}-10^{-20})$ g/cm$^3$ near the BH. We further analyze the possibility of observing this effect with radio interferometry, simulating observations with an EHT--like array, and find that it is unlikely to be detectable in the near future.

Journal ArticleDOI
Bei You1, Y. L. Tuo2, C. K. Li1, Wei Wang1, Shuang-Nan Zhang2, Shu Zhang2, Ming-Yu Ge2, Chong Luo1, Bifang Liu2, Weimin Yuan2, Zi-Gao Dai3, J. C. Liu2, J. C. Liu1, Erlin Qiao2, Chichuan Jin2, Zhu Liu2, Bozena Czerny4, Qingwen Wu5, Qingcui Bu6, Qingcui Bu2, C. Cai2, X. L. Cao2, Zhi Chang2, Gang Chen2, Li Chen7, Tian-Xiang Chen2, Y. P. Chen8, Yong Chen2, Yu-Peng Chen2, Wei Cui8, Wei-Wei Cui2, Jing-Kang Deng8, Yongwei Dong2, Yuan-Yuan Du2, Min-Xue Fu8, Guan-Hua Gao2, He Gao2, M. Gao2, Yu-Dong Gu2, Ju Guan2, Cheng-Cheng Guo2, Da-Wei Han2, Yue Huang2, Jia Huo2, Shu-Mei Jia2, Luhua Jiang2, W. C. Jiang2, Jing Jin2, YongJie Jin1, L. D. Kong2, Bing Li2, Cheng-Kui Li2, Gang Li2, Mao-Shun Li2, Ti-Pei Li2, Ti-Pei Li8, Wei Li2, Xian Li2, Xiaobo Li2, Xufang Li2, Yanguo Li2, Zheng-Wei Li2, Xiaohua Liang2, Jin-Yuan Liao2, C. X. Liu2, GuoQing Liu8, Hongwei Liu2, Xiao-Jing Liu2, Y. N. Liu8, Bo Lu2, Fangjun Lu2, Xue-Feng Lu2, Qi Luo2, T. Luo2, Xiang Ma2, Bin Meng2, Y. Nang2, Jianyin Nie2, Ge Ou2, Jin-Lu Qu2, N. Sai2, R. C. Shang8, Li-Ming Song2, X. Y. Song2, Liang Sun2, Ying Tan2, Lian Tao2, Chen Wang2, G. F. Wang2, Juan Wang2, Ling-Jun Wang2, Wenshuai Wang2, Yu-Sa Wang2, XiangYang Wen2, B. B. Wu2, Bobing Wu2, Mei Wu2, Guang-Cheng Xiao2, Shuo Xiao2, Shaolin Xiong2, Yupeng Xu2, J. W. Yang2, Sheng Yang2, Yanji Yang2, Qi-Bin Yi9, Qi-Bin Yi2, Qian-Qing Yin2, Yuan You2, Aimei Zhang2, Chengmo Zhang2, Fan Zhang2, Hong-Mei Zhang9, Juan Zhang2, Tong Zhang2, Wan-Chang Zhang2, Wei Zhang2, Wen-Zhao Zhang7, Yi Zhang2, Yifei Zhang2, Yongjie Zhang2, Yue Zhang2, Zhao Zhang8, Zi-Liang Zhang2, Hai-Sheng Zhao2, Xiao-Fan Zhao2, Shi-Jie Zheng2, D. K. Zhou2, Jian-Feng Zhou8, Yu-Xuan Zhu2, Yu-Xuan Zhu10, Yue Zhu2 
TL;DR: In this paper, a jet-like corona model was proposed, in which the corona can be understood as a standing shock where the material flowing through. But the model was not applied to the case of MAXI J1820+070, which is the ratio of the coronal intensity illuminating the disk to the intensity reaching the observer.
Abstract: A black hole X-ray binary produces hard X-ray radiation from its corona and disk when the accreting matter heats up. During an outburst, the disk and corona co-evolves with each other. However, such an evolution is still unclear in both its geometry and dynamics. Here we report the unusual decrease of the reflection fraction in MAXI J1820+070, which is the ratio of the coronal intensity illuminating the disk to the coronal intensity reaching the observer, as the corona is observed to contrast during the decay phase. We postulate a jet-like corona model, in which the corona can be understood as a standing shock where the material flowing through. In this dynamical scenario, the decrease of the reflection fraction is a signature of the corona's bulk velocity. Our findings suggest that as the corona is observed to get closer to the black hole, the coronal material might be outflowing faster.

Journal ArticleDOI
TL;DR: In this paper, the authors presented a catalog of central engine properties, i.e., black hole mass (mbh) and accretion luminosity (ld), for a sample of 1077 blazars detected with the Fermi~Large Area Telescope.
Abstract: We present a catalog of central engine properties, i.e., black hole mass (mbh) and accretion luminosity (ld), for a sample of 1077 blazars detected with the Fermi~Large Area Telescope. This includes broad emission line systems and blazars whose optical spectra lack emission lines but dominated by the absorption features arising from the host galaxy. The average mbh for the sample is $\langle \log~M_{{\rm BH,all}~msun} \rangle=8.60$ and there are evidences suggesting the association of more massive black holes with absorption line systems. Our results indicate a bi-modality of ld~in Eddington units (ld/ledd) with broad line objects tend to have a higher accretion rate (ld/ledd$>$0.01). We have found that ld/ledd and Compton dominance (CD, the ratio of the inverse Compton to synchrotron peak luminosities) are positively correlated at $>$5$\sigma$ confidence level, suggesting that the latter can be used to reveal the state of accretion in blazars. Based on this result, we propose a CD based classification scheme. Sources with CD$>$1 can be classified as High-Compton Dominated or HCD blazars, whereas, that with CD$\lesssim$1 are Low-Compton Dominated (LCD) objects. This scheme is analogous to that based on the mass accretion rate proposed in the literature, however, it overcomes the limitation imposed by the difficulty in measuring ld and mbh for objects with quasi-featureless spectra. We conclude that the overall physical properties of Fermi blazars are likely to be controlled by the accretion rate in Eddington units. The catalog is made public at this http URL and Zenodo.

Posted Content
TL;DR: In this paper, the authors presented a localisation of the repeating source FRB 20201124A, and its association with a host galaxy (SDSS J050803.48+260338.0, z=0.098) and persistent radio source.
Abstract: The physical properties of fast radio burst (FRB) host galaxies provide important clues towards the nature of FRB sources. The 16 FRB hosts identified thus far span three orders of magnitude in mass and specific star-formation rate, implicating a ubiquitously occurring progenitor object. FRBs localised with ~arcsecond accuracy also enable effective searches for associated multi-wavelength and multi-timescale counterparts, such as the persistent radio source associated with FRB 20121102A. Here we present a localisation of the repeating source FRB 20201124A, and its association with a host galaxy (SDSS J050803.48+260338.0, z=0.098) and persistent radio source. The galaxy is massive ($\sim3\times10^{10} M_{\odot}$), star-forming (few solar masses per year), and dusty. Very Large Array and Very Long Baseline Array observations of the persistent radio source measure a luminosity of $1.2\times10^{29}$ erg s$^{-1}$ Hz$^{-1}$, and show that is extended on scales $\gtrsim50$ mas. We associate this radio emission with the ongoing star-formation activity in SDSS J050803.48+260338.0. Deeper, more detailed observations are required to better utilise the milliarcsecond-scale localisation of FRB 20201124A reported from the European VLBI Network, and determine the origin of the large dispersion measure ($150-220$ pc cm$^{-3}$) contributed by the host. SDSS J050803.48+260338.0 is an order of magnitude more massive than any galaxy or stellar system previously associated with a repeating FRB source, but is comparable to the hosts of so far non-repeating FRBs, further building the link between the two apparent populations.

Posted Content
TL;DR: In this article, the authors used a long GRB sample with X-ray plateau phases produced by electromagnetic dipole emissions from central new-born magnetars and found a tight correlation between the plateau luminosity and the end time of the plateau in Xray afterglows out to the redshift.
Abstract: Gamma-ray bursts (GRBs) are the most luminous explosions and can be detectable out to the edge of Universe. It has long been thought they can extend the Hubble diagram to very high redshifts. Several correlations between temporal or spectral properties and GRB luminosities have been proposed to make GRBs cosmological tools. However, those correlations cannot be properly standardized. In this paper, we select a long GRB sample with X-ray plateau phases produced by electromagnetic dipole emissions from central new-born magnetars. A tight correlation is found between the plateau luminosity and the end time of the plateau in X-ray afterglows out to the redshift $z=5.91$. We standardize these long GRBs X-ray light curves to a universal behavior by this correlation for the first time, with a luminosity dispersion of 0.5 dex. The derived distance-redshift relation of GRBs is in agreement with the standard $\Lambda$CDM model both at low and high redshifts. The evidence of accelerating universe from this GRB sample is $3\sigma$, which is the highest statistical significance from GRBs to date.