Search

The stellar stripping of satellites in cluster haloes is understood to play an important role in the production of intracluster light. Increasingly, cosmological simulations have been utilised to investigate its origin and assembly. However, such simulations typically model individual galaxies at relatively coarse resolutions, raising concerns about their accuracy. Although there is a growing literature on the importance of numerical resolution for the accurate recovery of the mass loss rates of dark matter (DM) haloes, there has been no comparable investigation into the numerical resolution required to accurately recover stellar mass loss rates in galaxy clusters. Using N-body simulations of satellite galaxies orbiting in a cluster halo represented by a static external potential, we conduct a set of convergence tests in order to explore the role of numerical resolution and force softening length on stellar stripping efficiency. We consider a number of orbital configurations, satellite masses and satellite morphologies. We find that stellar mass resolution is of minor importance relative to DM resolution. Resolving the central regions of satellite DM halos is critical to accurately recover stellar mass loss rates. Poorly resolved DM haloes develop cored inner profiles and, if this core is of comparable size to the stellar component of the satellite galaxy, this leads to significant over-stripping. To prevent this, relatively high DM mass resolutions of around $m_{\rm DM}\sim10^{6}$ M$_{\odot}$, better than those achieved by many contemporary cosmological simulations, are necessary., Comment: Accepted for publication in MNRAS. 19 pages, 13 figures

According to CMB measurements, baryonic matter constitutes about $5\%$ of the mass-energy density of the universe. A significant population of these baryons, for a long time referred to as `missing', resides in a low density, warm-hot intergalactic medium (WHIM) outside galaxy clusters, tracing the ``cosmic web'', a network of large scale dark matter filaments. Various studies have detected this inter-cluster gas, both by stacking and by observing individual filaments in compact, massive systems. In this paper, we study short filaments (< 10 Mpc) connecting massive clusters ($M_{500} \approx 3\times 10^{14} M_{\odot}$) detected by the Atacama Cosmology Telescope (ACT) using the scattering of CMB light off the ionised gas, a phenomenon known as the thermal Sunyaev-Zeldovich (tSZ) effect. The first part of this work is a search for suitable candidates for high resolution follow-up tSZ observations. We identify four cluster pairs with an intercluster signal above the noise floor (S/N $>$ 2), including two with a tentative $>2\sigma$ statistical significance for an intercluster bridge from the ACT data alone. In the second part of this work, starting from the same cluster sample, we directly stack on ${\sim}100$ cluster pairs and observe an excess SZ signal between the stacked clusters of $y=(7.2^{+2.3}_{-2.5})\times 10^{-7}$ with a significance of $3.3\sigma$. It is the first tSZ measurement of hot gas between clusters in this range of masses at moderate redshift ($\langle z\rangle\approx 0.5$). We compare this to the signal from simulated cluster pairs with similar redshifts and separations in the THE300 and MAGNETICUM Pathfinder cosmological simulations and find broad consistency. Additionally, we show that our measurement is consistent with scaling relations between filament parameters and mass of the embedded halos identified in simulations., Comment: 37 pages, 17 images

Massive stars end their life as core-collapse supernovae, amongst which some extremes are Type Ic broad-lined supernovae associated with long-duration gamma-ray bursts (LGRBs) having powerful relativistic jets. Their less-extreme brethren make unsuccessful jets that are choked inside the stars, appearing as X-ray flashes or low-luminosity GRBs. On the other hand, there exists a population of extragalactic fast X-ray transients (EFXTs) with timescales ranging from seconds to thousands of seconds, whose origins remain obscure. Known sources that contribute to the observed EFXT population include the softer analogs of LGRBs, shock breakouts of supernovae, or unsuccessful jets. Here, we report the discovery of the bright X-ray transient EP240414a detected by the Einstein Probe (EP), which is associated with the Type Ic supernova SN 2024gsa at a redshift of 0.401. The X-ray emission evolution is characterised by a very soft energy spectrum peaking at < 1.3 keV, which makes it distinct from known LGRBs, X-ray flashes, or low-luminosity GRBs. Follow-up observations at optical and radio bands revealed the existence of a weak relativistic jet that interacts with an extended shell surrounding the progenitor star. Located on the outskirts of a massive galaxy, this event reveals a new population of explosions of Wolf-Rayet stars characterised by a less powerful engine that drives a successful but weak jet, possibly owing to a progenitor star with a smaller core angular momentum than in traditional LGRB progenitors., Comment: 43 pages, 9 figures, 4 tables, submitted. Comments are welcome

The unique characteristics of the brightest group galaxies (BGGs) link the evolutionary continuum between galaxies like the Milky Way and more massive BCGs in dense clusters. This study investigates the stellar properties of BGGs over cosmic time (z = 0.08-1.30), extending our previous work (Gozaliasl et al. 2016, 2018; Paper I and Paper II). We analyze data of 246 BGGs from our X-ray galaxy group catalog in the COSMOS field, examining stellar age, mass, star formation rate (SFR), specific SFR (sSFR), and halo mass. Comparisons are made with Millennium and Magneticum simulations. We explore the variation of stellar properties with the projected offset from the X-ray peak or host halo center. Using a mock galaxy catalog, we evaluated the accuracy of SED-derived stellar ages, finding a mean absolute error of about one Gyr. Observed BGG age distributions show a bias towards younger ages compared to semi-analytical models and the Magneticum simulation. Our analysis of stellar age versus mass reveals trends with a positive slope, suggesting complex evolutionary pathways across redshifts. We observe a negative correlation between stellar age and SFR across all redshift ranges. Using a cosmic-time-dependent main sequence framework, we identify star-forming BGGs, finding that about 20% of BGGs in the local universe exhibit star-forming characteristics, increasing to 50% at $z=1.0$. Our findings support an inside-out formation scenario for BGGs, where older stellar populations are near the X-ray peak and younger populations at larger offsets indicate ongoing star formation. The distribution of stellar ages for lower-mass BGGs ($10^{10-11} M_\odot$) deviates from constant ages predicted by models, highlighting current models' limitations in capturing galaxies' complex star formation histories., Comment: 24 Pages, 45 figures

Long gamma-ray bursts (GRBs) are believed to originate from core collapse of massive stars. High-redshift GRBs can probe the star formation and reionization history of the early universe, but their detection remains rare. Here we report the detection of a GRB triggered in the 0.5--4 keV band by the Wide-field X-ray Telescope (WXT) on board the Einstein Probe (EP) mission, designated as EP240315a, whose bright peak was also detected by the Swift Burst Alert Telescope and Konus-Wind through off-line analyses. At a redshift of $z=4.859$, EP240315a showed a much longer and more complicated light curve in the soft X-ray band than in gamma-rays. Benefiting from a large field-of-view ($\sim$3600 deg$^2$) and a high sensitivity, EP-WXT captured the earlier engine activation and extended late engine activity through a continuous detection. With a peak X-ray flux at the faint end of previously known high-$z$ GRBs, the detection of EP240315a demonstrates the great potential for EP to study the early universe via GRBs., Comment: 41 pages, 8 figures, 7 tables

The abundance of galaxy clusters with mass and redshift is a well-known cosmological probe. The cluster mass is a key parameter for studies that aim to constrain cosmological parameters using galaxy clusters, making it critical to understand and properly account for the errors in its estimates. Subsequently, it becomes important to correctly calibrate scaling relations between observables like the integrated Compton parameter and the mass of the cluster. The NIKA2 Sunyaev-Zeldovich Large program (LPSZ) enables one to map the intracluster medium profiles in the mm-wavelength band with great details (resolution of $11 \ \mathrm{\&}\ 17^{\prime \prime}$ at $1.2 \ \mathrm{\&}\ 2 $ mm, respectively) and hence, to estimate the cluster hydrostatic mass more precisely than previous SZ observations. However, there are certain systematic effects which can only be accounted for with the use of simulations. For this purpose, we employ THE THREE HUNDRED simulations which have been modelled with a range of physics modules to simulate galaxy clusters. The so-called twin samples are constructed by picking synthetic clusters of galaxies with properties close to the observational targets of the LPSZ. In particular, we use the Compton parameter maps and projected total mass maps of these twin samples along 29 different lines of sight. We investigate the scatter that projection induces on the total masses. Eventually, we consider the statistical values along different lines of sight to construct a kind of 3D scaling law between the integrated Compton parameter, total mass, and overdensity of the galaxy clusters to determine the overdensity that is least impacted by the projection effect., Comment: to appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

Cluster number counts will be a key cosmological probe in the next decade thanks to the Euclid satellite mission. For this purpose, cluster detection algorithm performance, which are sensitive to the spatial distribution of the cluster galaxy members and their luminosity function, need to be accurately characterized. Using The Three Hundred hydrodynamical and dark matter only simulations we study a complete sample of massive clusters beyond 7 (5) $\times$ 10$^{14}$ M$_{\odot}$ at redshift 0 (1) on a $(1.48 \ \mathrm{Gpc})^3$ volume. We find that the mass resolution of the current hydrodynamical simulations (1.5 $\times$ 10$^9$ M$_{\odot}$) is not enough to characterize the luminosity function of the sample in the perspective of Euclid data. Nevertheless, these simulations are still useful to characterize the spatial distribution of the cluster substructures assuming a common relative mass threshold for the different flavours and resolutions. By comparing with the dark matter only version of these simulations, we demonstrate that baryonic physics preserves significantly low mass subhalos (galaxies) as have also been observed in previous studies with less statistics. Furthermore, by comparing the hydro simulations with higher resolution dark matter only simulations of the same objects and taking the same limit in subhalo mass we find significantly more cuspy galaxy density profiles towards the center of the clusters, where the low mass substructures would tend to concentrate. We conclude that using dark matter only simulation may lead to some biases on the spatial distribution and density of galaxy cluster members. Based on the preliminary analysis of few high resolution hydro simulations we conclude that a mass resolution of 1.8 $\times$ 10$^8$ h$^{-1}$ M$_{\odot}$ will be needed for The Three Hundred simulations to approach the expected magnitude limits for the Euclid survey.

The determination of the mass of galaxy clusters from observations is subject to systematic uncertainties. Beyond the errors due to instrumental and observational systematic effects, in this work we investigate the bias introduced by modelling assumptions. In particular, we consider the reconstruction of the mass of galaxy clusters from convergence maps employing spherical mass density models. We made use of The Three Hundred simulations, selecting clusters in the same redshift and mass range as the NIKA2 Sunyaev-Zel'dovich Large Programme sample: $3 \leq M_{500}/ 10^{14} \mathrm{M}_{\odot} \leq 10$ and $0.5 \leq z \leq 0.9$. We studied different modelling and intrinsic uncertainties that should be accounted for when using the single cluster mass estimates for scaling relations. We confirm that the orientation of clusters and the radial ranges considered for the fit have an important impact on the mass bias. The effect of the projection adds uncertainties to the order of $10\%$ to $16\%$ to the mass estimates. We also find that the scatter from cluster to cluster in the mass bias when using spherical mass models is less than $9\%$ of the true mass of the clusters.

GRB 221009A is the brightest Gamma-Ray Burst (GRB) detected in more than 50 years of study. In this paper, we present observations in the X-ray and optical domains after the GRB obtained by the GRANDMA Collaboration (which includes observations from more than 30 professional and amateur telescopes) and the Insight-HXMT Collaboration. We study the optical afterglow with empirical fitting from GRANDMA+HXMT data, augmented with data from the literature up to 60 days. We then model numerically, using a Bayesian approach, the GRANDMA and HXMT-LE afterglow observations, that we augment with Swift-XRT and additional optical/NIR observations reported in the literature. We find that the GRB afterglow, extinguished by a large dust column, is most likely behind a combination of a large Milky-Way dust column combined with moderate low-metallicity dust in the host galaxy. Using the GRANDMA+HXMT-LE+XRT dataset, we find that the simplest model, where the observed afterglow is produced by synchrotron radiation at the forward external shock during the deceleration of a top-hat relativistic jet by a uniform medium, fits the multi-wavelength observations only moderately well, with a tension between the observed temporal and spectral evolution. This tension is confirmed when using the extended dataset. We find that the consideration of a jet structure (Gaussian or power-law), the inclusion of synchrotron self-Compton emission, or the presence of an underlying supernova do not improve the predictions, showing that the modelling of GRB22109A will require going beyond the most standard GRB afterglow model. Placed in the global context of GRB optical afterglows, we find the afterglow of GRB 221009A is luminous but not extraordinarily so, highlighting that some aspects of this GRB do not deviate from the global known sample despite its extreme energetics and the peculiar afterglow evolution., Comment: Accepted to ApJL for the special issue, 37 pages, 23 pages main text, 6 tables, 13 figures

The photometric catalogue of galaxy clusters extracted from ESA Euclid data is expected to be very competitive for cosmological studies. Using state-of-the-art hydrodynamical simulations, we present systematic analyses simulating the expected weak lensing profiles from clusters in a variety of dynamic states and at wide range of redshifts. In order to derive cluster masses, we use a model consistent with the implementation within the Euclid Consortium of the dedicated processing function and find that, when jointly modelling mass and the concentration parameter of the Navarro-Frenk-White halo profile, the weak lensing masses tend to be, on average, biased low by 5-10% with respect to the true mass, up to z=0.5. Using a fixed value for the concentration $c_{200} = 3$, the mass bias is diminished below 5%, up to z=0.7, along with its relative uncertainty. Simulating the weak lensing signal by projecting along the directions of the axes of the moment of inertia tensor ellipsoid, we find that orientation matters: when clusters are oriented along the major axis, the lensing signal is boosted, and the recovered weak lensing mass is correspondingly overestimated. Typically, the weak lensing mass bias of individual clusters is modulated by the weak lensing signal-to-noise ratio, related to the redshift evolution of the number of galaxies used for weak lensing measurements: the negative mass bias tends to be larger toward higher redshifts. However, when we use a fixed value of the concentration parameter, the redshift evolution trend is reduced. These results provide a solid basis for the weak-lensing mass calibration required by the cosmological application of future cluster surveys from Euclid and Rubin., Comment: Accepted for publication in Astronomy & Astrophysics

Soft x-ray emissions induced by solar wind ions that collide with neutral material in the solar system have been detected around planets, and were proposed as a remote probe for the solar wind interaction with the Martian exosphere. A multi-fluid three-dimensional magneto-hydrodynamic model is adopted to derive the global distributions of solar wind particles. Spherically symmetric exospheric H, H$_2$, He, O, and CO$_2$ density profiles and a sophisticated hybrid model that includes charge-exchange and proton/neutral excitation processes are used to study the low triplet line ratio $G=\frac{i+f}{r}$ (0.77$\pm$0.58) of O VII and total x-ray luminosity around Mars. We further calculate the emission factor $\alpha$-value with different neutrals over a wide ion abundance and velocity ranges. Our results are in good agreement with those of previous reports. The evolution of the charge stage of solar wind ions shows that sequential recombination due to charge-exchange can be negligible at the interaction region. This only appears below the altitude of 400~km. The anonymous low disk $G$ ratio can be easily explained by the collisional quenching effect at neutral densities higher than 10$^{11}$cm$^{-3}$. However, the quenching contribution is small in Mars' exosphere and only appears below 400~km. Charge-exchange with H$_2$ and N$_2$ is still the most likely reason for this low $G$-ratio. X-ray emissivity maps in collisions with different neutrals differ from each other. A clear bow shock in the collision with all the neutrals is in accordance with previous reports. The resulting total x-ray luminosity of 6.55~MW shows a better agreement with the XMM-Newton observation of 12.8$\pm$1.4~MW than that of previous predictions., Comment: 17 pages, 13 figures; Accpted by ApJ

We develop a machine learning algorithm to infer the 3D cumulative radial profiles of total and gas mass in galaxy clusters from thermal Sunyaev-Zel'dovich effect maps. We generate around 73,000 mock images along various lines of sight using 2,522 simulated clusters from the \thethreehundred{} project at redshift $z< 0.12$ and train a model that combines an autoencoder and a random forest. Without making any prior assumptions about the hydrostatic equilibrium of the clusters, the model is capable of reconstructing the total mass profile as well as the gas mass profile, which is responsible for the SZ effect. We show that the recovered profiles are unbiased with a scatter of about $10\%$, slightly increasing towards the core and the outskirts of the cluster. We selected clusters in the mass range of $10^{13.5} \leq M_{200} /(\hMsun) \leq 10^{15.5}$, spanning different dynamical states, from relaxed to disturbed halos. We verify that both the accuracy and precision of this method show a slight dependence on the dynamical state, but not on the cluster mass. To further verify the consistency of our model, we fit the inferred total mass profiles with an NFW model and contrast the concentration values with those of the true profiles. We note that the inferred profiles are unbiased for higher concentration values, reproducing a trustworthy mass-concentration relation. The comparison with a widely used mass estimation technique, such as hydrostatic equilibrium, demonstrates that our method recovers the total mass that is not biased by non-thermal motions of the gas., Comment: Accepted for publication on Monthly Notices of the Royal Astronomical Society, 8 pages, 10 figures

We present a spectral study of the black hole candidate MAXI J1348-630 during its 2019 outburst, based on monitoring observations with Insight-HXMT and Swift. Throughout the outburst, the spectra are well fitted with power-law plus disk-blackbody components. In the soft-intermediate and soft states, we observed the canonical relation L ~ T_in^4 between disk luminosity L and peak colour temperature T_in, with a constant inner radius R_in (traditionally identified with the innermost stable circular orbit). At other stages of the outburst cycle, the behaviour is more unusual, inconsistent with the canonical outburst evolution of black hole transients. In particular, during the hard rise, the apparent inner radius is smaller than in the soft state (and increasing), and the peak colour temperature is higher (and decreasing). This anomalous behaviour is found even when we model the spectra with self-consistent Comptonization models, which take into account the up-scattering of photons from the disk component into the power-law component. To explain both those anomalous trends at the same time, we suggest that the hardening factor for the inner disk emission was larger than the canonical value of ~1.7 at the beginning of the outburst. A more physical trend of radii and temperature evolution requires a hardening factor evolving from ~3.5 at the beginning of the hard state to ~1.7 in the hard intermediate state. This could be evidence that the inner disk was in the process of condensing from the hot, optically thin medium and had not yet reached a sufficiently high optical depth for its emission spectrum to be described by the standard optically-thick disk solution., Comment: 20 pages, 10 figures, 4 tables, accepted by ApJ

Based on Insight-HXMT data, we report on the pulse fraction evolution during the 2017-2018 outburst of the newly discovered first Galactic ultraluminous X-ray source (ULX) Swift J0243.6+6124. The pulse fractions of 19 observation pairs selected in the rising and fading phases with similar luminosity are investigated. The results show a general trend of the pulse fraction increasing with luminosity and energy at super-critical luminosity. However, the relative strength of the pulsation between each pair evolves strongly with luminosity. The pulse fraction in the rising phase is larger at luminosity below $7.71\times10^{38}$~erg~s$^{-1}$, but smaller at above. A transition luminosity is found to be energy independent. Such a phenomena is firstly confirmed by Insight-HXMT observations and we speculate it may have relation with the radiation pressure dominated accretion disk.

We search for the signature of cosmological shocks in stacked gas pressure profiles of galaxy clusters using data from the South Pole Telescope (SPT). Specifically, we stack the latest Compton-y maps from the 2500 deg^2 SPT-SZ survey on the locations of clusters identified in that same dataset. The sample contains 516 clusters with mean mass = 1e14.9 msol and redshift = 0.55. We analyze in parallel a set of zoom-in hydrodynamical simulations from The Three Hundred project. The SPT-SZ data show two features: (i) a pressure deficit at R/R200m = $1.08 \pm 0.09$, measured at $3.1\sigma$ significance and not observed in the simulations, and; (ii) a sharp decrease in pressure at R/R200m = $4.58 \pm 1.24$ at $2.0\sigma$ significance. The pressure deficit is qualitatively consistent with a shock-induced thermal non-equilibrium between electrons and ions, and the second feature is consistent with accretion shocks seen in previous studies. We split the cluster sample by redshift and mass, and find both features exist in all cases. There are also no significant differences in features along and across the cluster major axis, whose orientation roughly points towards filamentary structure. As a consistency test, we also analyze clusters from the Planck and Atacama Cosmology Telescope Polarimeter surveys and find quantitatively similar features in the pressure profiles. Finally, we compare the accretion shock radius (Rsh_acc) with existing measurements of the splashback radius (Rsp) for SPT-SZ and constrain the lower limit of the ratio, Rsh_acc/Rsp > $2.16 \pm 0.59$., Comment: [v1]: 8 Figures, 16 Pages in Main text. [v2]: Added text to discussion. Version accepted in MNRAS

The simulation database of THE THREE HUNDRED Project has been used to pick synthetic clusters of galaxies with properties close to the observational targets of the NIKA2 camera Sunyaev-Zeldovich (SZ) Large Program. Cross-matching of cluster parameters such as mass and redshift of the cluster in the two databases has been implemented to generate the so-called twin samples for the Large Program. This SZ Large Program is observing a selection of galaxy clusters at intermediate and high redshift $\left( 0.5 < z < 0.9 \right)$, covering one order of magnitude in mass. These are SZ-selected clusters from the Planck and Atacama Cosmology Telescope catalogs, wherein the selection is based on their integrated Compton parameter values, $Y_{500}$: the value of the parameter within the characteristics radius $R_{500}$. THE THREE HUNDRED hydrodynamical simulations provide us with hundreds of clusters satisfying these redshift, mass, and $Y_{500}$ requirements. In addition to the standard post-processing analysis of the simulation, mock observational maps are available mimicking X-ray, optical, gravitational lensing, radio, and SZ observations of galaxy clusters. The primary goal of employing the twin samples is to compare different cluster mass proxies from synthetic X-ray, SZ effect and optical maps (via the velocity dispersion of member galaxies and lensing $\kappa$-maps) of the clusters. Eventually, scaling laws between different mass proxies and the cluster mass will be cross-correlated to reduce the scatter on the inferred mass and the mass bias will be related to various physical parameters., Comment: To appear in the Proceedings of the International Conference entitled "mm Universe @ NIKA2", Rome (Italy), June 2021, EPJ Web of conferences

We investigate how large-scale cosmic filaments impact the quenching of galaxies within one virial radius of 324 simulated clusters from The Three Hundred project. We track cosmic filaments with the versatile, observation-friendly program DisPerSE and identify halos hosting galaxies with VELOCIRaptor. We find that cluster galaxies close to filaments tend to be more star-forming, bluer, and contain more cold gas than their counterparts further away from filaments. This effect is recovered at all stellar masses. This is in stark contrast with galaxies residing outside of clusters, where galaxies close to filaments show clear signs of density related pre-processing. We first show that the density contrast of filaments is reduced inside the intra-cluster medium. Moreover, examination of flows around and into cluster galaxies shows that the gas flows in intra-cluster filaments are colder and tend to stream along with galaxies in their midst, partially shielding them from strangulation by the hot, dense intra-cluster medium. This also preserves accretion onto satellites and limit ram pressure., Comment: 20 pages, 21 figures, submitted to MNRAS

The High Energy X-ray telescope (HE) on-board the Hard X-ray Modulation Telescope (\textit{Insight}-HXMT) can serve as a wide Field of View (FOV) gamma-ray monitor with high time resolution ($\mu$s) and large effective area (up to thousands cm$^2$). We developed a pipeline to search for Gamma-Ray Bursts (GRBs), using the traditional signal-to-noise ratio (SNR) method for blind search and the coherent search method for targeted search. By taking into account the location and spectrum of the burst and the detector response, the targeted coherent search is more powerful to unveil weak and sub-threshold bursts, especially those in temporal coincidence with Gravitational Wave (GW) events. Based on the original method in literature, we further improved the coherent search to filter out false triggers caused by spikes in light curves, which are commonly seen in gamma-ray instruments (e.g. \textit{Fermi}/GBM, \textit{POLAR}). We show that our improved targeted coherent search method could eliminate almost all false triggers caused by spikes. Based on the first two years of \textit{Insight}-HXMT/HE data, our targeted search recovered 40 GRBs, which were detected by either \textit{Swift}/BAT or \textit{Fermi}/GBM but too weak to be found in our blind search. With this coherent search pipeline, the GRB detection sensitivity of \textit{Insight}-HXMT/HE is increased to about 1.5E-08 erg/cm$^2$ (200 keV$-$3 MeV). We also used this targeted coherent method to search \textit{Insight}-HXMT/HE data for electromagnetic (EM) counterparts of LIGO-Virgo GW events (including O2 and O3a runs). However, we did not find any significant burst associated with GW events., Comment: 12 pages, 14 figures, 5 tables; accepted for publication in MNRAS

We report on the observation of the accreting pulsar GRO J1008-57 performed by Insight-HXMT at the peak of the source's 2017 outburst. Pulsations are detected with a spin period of 93.283(1) s. The pulse profile shows double peaks at soft X-rays, and only one peak above 20 keV. The spectrum is well described by the phenomenological models of X-ray pulsars. A cyclotron resonant scattering feature is detected with very high statistical significance at a centroid energy of $E_{\rm cyc}=90.32_{-0.28}^{+0.32}$ keV, for the reference continuum and line models, HIGHECUT and GABS respectively. Detection is very robust with respect to different continuum models. The line energy is significantly higher than what suggested from previous observations, which provided very marginal evidence for the line. This establishes a new record for the centroid energy of a fundamental cyclotron resonant scattering feature observed in accreting pulsars. We also discuss the accretion regime of the source during the Insight-HXMT observation., Comment: 8 pages, 3 figures, accepted for publication in ApJL

Fast radio bursts (FRBs) are short pulses observed in radio band from cosmological distances. One class of models invoke soft gamma-ray repeaters (SGRs), or magnetars, as the sources of FRBs. Some radio pulses have been observed from some magnetars, however, no FRB-like events had been detected in association any magnetar burst, including one giant flare. Recently, a pair of FRB-like bursts (FRB 200428 hereafter) separated by milliseconds (ms) were detected from the general direction of the Galactic magnetar SGR J1935+2154. Here we report the detection of a non-thermal X-ray burst in the 1-250 keV energy band with the Insight-HXMT satellite, which we identify as emitted from SGR J1935+2154. The burst showed two hard peaks with a separation of 34 ms, broadly consistent with that of the two bursts in FRB 200428. The delay time between the double radio and X-ray peaks is about 8.57 s, fully consistent with the dispersion delay of FRB 200428. We thus identify the non-thermal X-ray burst is associated with FRB 200428 whose high energy counterpart is the two hard peaks in X-ray. Our results suggest that the non-thermal X-ray burst and FRB 200428 share the same physical origin in an explosive event from SGR J1935+2154., Comment: 24 pages, 9 figures, 6 tables; initial submission to a journal on May 9th, 2020. Significant changes include updated localization and detailed spectral evolution of the X-ray burst, and better determination of the two narrow X-ray peaks corresponding to the two radio pulses. Conclusions are strengthened. Nature Astronomy online on Feb. 18, 2021

We use the In data collected during the 2019 outburst from X-ray pulsar 4U 1901+03 to complement the orbital parameters reported by Fermi/GBM. Using the Insight-HXMT, we examine the correlation between the derivative of the intrinsic spin frequency and bolometric flux based on accretion torque models. It was found that the pulse profiles significantly evolve during the outburst. The existence of two types of the profile's pattern discovered in the Insight-HXMT data indicates that this source experienced transition between a super-critical and a sub-critical accretion regime during its 2019 outburst. Based on the evolution of the pulse profiles and the torque model, we derive the distance to 4U 1901+03 as 12.4+-0.2 kpc., Comment: 8 pages, 5 figures, accepted by JHEAP

We study the evolution of the temporal properties of MAXI 1820+070 during the 2018 outburst in its hard state from MJD 58190 to 58289 with Insight-HXMT in a broad energy band 1-150 keV. We find different behaviors of the hardness ratio, the fractional rms and time lag before and after MJD 58257, suggesting a transition occurred at around this point. The observed time lags between the soft photons in the 1-5 keV band and the hard photons in higher energy bands, up to 150 keV, are frequency-dependent: the time lags in the low-frequency range, 2-10 mHz, are both soft and hard lags with a timescale of dozens of seconds but without a clear trend along the outburst; the time lags in the high-frequency range, 1-10 Hz, are only hard lags with a timescale of tens of milliseconds; first increase until around MJD 58257 and decrease after this date. The high-frequency time lags are significantly correlated to the photon index derived from the fit to the quasi-simultaneous NICER spectrum in the 1-10 keV band. This result is qualitatively consistent with a model in which the high-frequency time lags are produced by Comptonization in a jet., Comment: Accepted for publication in ApJ

Glitches correspond to sudden jumps of rotation frequency ($\nu$) and its derivative ($\dot{\nu}$) of pulsars, the origin of which remains not well understood yet, partly because the jump processes of most glitches are not well time-resolved. There are three large glitches of the Crab pulsar, detected in 1989, 1996 and 2017, which were found to have delayed spin-up processes before the normal recovery processes. Here we report two additional glitches of the Crab pulsar occurred in 2004 and 2011 for which we discovered delayed spin up processes, and present refined parameters of the largest glitch occurred in 2017. The initial rising time of the glitch is determined as $<0.48$ hour. We also carried out a statistical study of these five glitches with observed spin-up processes. The two glitches occurred in 2004 and 2011 have delayed spin-up time scales ($\tau_{1}$) of $1.7\pm0.8$\,days and $1.6\pm0.4$\,days, respectively. We find that the $\Delta{\nu}$ vs. $|\Delta{\dot\nu}|$ relation of these five glitches is similar to those with no detected delayed spin-up process, indicating that they are similar to the others in nature except that they have larger amplitudes. For these five glitches, the amplitudes of the delayed spin-up process ($|\Delta{\nu}_{\rm d1}|$) and recovery process ($\Delta{\nu}_{\rm d2}$), their time scales ($\tau_{1}$, $\tau_{2}$), and permanent changes in spin frequency ($\Delta{\nu}_{\rm p}$) and total frequency step ($\Delta{\nu}_{\rm g}$) have positive correlations. From these correlations, we suggest that the delayed spin-up processes are common for all glitches, but are too short and thus difficult to be detected for most glitches., Comment: 25 pages, 8 figures

No robust detection of prompt electromagnetic counterparts to fast radio bursts (FRBs) has yet been obtained, in spite of several multi-wavelength searches carried out so far. Specifically, X/gamma-ray counterparts are predicted by some models. We planned on searching for prompt gamma-ray counterparts in the Insight-Hard X-ray Modulation Telescope (Insight-HXMT) data, taking advantage of the unique combination of large effective area in the keV-MeV energy range and of sub-ms time resolution. We selected 39 FRBs that were promptly visible from the High-Energy (HE) instrument aboard Insight-HXMT. After calculating the expected arrival times at the location of the spacecraft, we searched for a significant excess in both individual and cumulative time profiles over a wide range of time resolutions, from several seconds down to sub-ms scales. Using the dispersion measures in excess of the Galactic terms, we estimated the upper limits on the redshifts. No convincing signal was found and for each FRB we constrained the gamma-ray isotropic-equivalent luminosity and the released energy as a function of emission timescale. For the nearest FRB source, the periodic repeater FRB180916.J0158+65, we find $L_{\gamma,iso}<5.5\times 10^{47}$ erg/s over 1 s, whereas $L_{\gamma,iso}<10^{49}-10^{51}$ erg/s for the bulk of FRBs. The same values scale up by a factor of ~100 for a ms-long emission. Even on a timescale comparable with that of the radio pulse itself no keV-MeV emission is observed. A systematic association with either long or short GRBs is ruled out with high confidence, except for subluminous events, as is the case for core-collapse of massive stars (long) or binary neutron star mergers (short) viewed off axis. Only giant flares from extra-galactic magnetars at least ten times more energetic than Galactic siblings are ruled out for the nearest FRB., Comment: 15 pages, 3 figures, 6 tables, accepted by A&A

We report on our analysis of the 2019 outburst of the X-ray accreting pulsar 4U 1901+03 observed with Insight-HXMT and NICER. Both spectra and pulse profiles evolve significantly in the decaying phase of the outburst. Dozens of flares are observed throughout the outburst. They are more frequent and brighter at the outburst peak. We find that the flares, which have a duration from tens to hundreds of seconds, are generally brighter than the persistent emission by a factor of $\sim$ 1.5. The pulse profile shape during the flares can be significantly different than that of the persistent emission. In particular, a phase shift is clearly observed in many cases. We interpret these findings as direct evidence of changes of the pulsed beam pattern, due to transitions between the sub- and super-critical accretion regimes on a short time scale. We also observe that at comparable luminosities the flares' pulse profiles are rather similar to those of the persistent emission. This indicates that the accretion on the polar cap of the neutron star is mainly determined by the luminosity, i.e., the mass accretion rate., Comment: 11 pages, 8 figures, accepted for publication in MNRAS

A new black hole X-ray binary (BHXRB) MAXI J1535-571 was discovered by MAXI during its outburst in 2017. Using observations taken by the first Chinese X-ray satellite, the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT), we perform a joint spectral analysis (2-150 keV) in both energy and time domains. The energy spectra provide the essential input for probing the intrinsic Quasi-Periodic Oscillation (QPO) fractional rms spectra (FRS). Our results show that during the intermediate state, the energy spectra are in general consistent with those reported by Swift/XRT and NuSTAR. However, the QPO FRS become harder and the FRS residuals may suggest the presence of either an additional power-law component in the energy spectrum or a turn-over in the intrinsic QPO FRS at high energies.

We study the 2018 outburst of Aql X-1 via the monitor of all sky X-ray image (MAXI) data. We show that the outburst starting in February 2018 is a member of short-low class in the frame of outburst duration and the peak count rate although the outburst morphology is slightly different from the other fast-rise-exponential-decay (FRED) type outbursts with a milder rising stage. We study the partial accretion in the weak propeller stage of Aql X-1 via the MAXI data of the 2018 outburst. We report on the spectral analysis of 3 observations of Aquila X-1 obtained by Insight - hard X-ray modulation telescope (Insight-HXMT) during the late decay stage of the 2018 outburst. We discuss that the data taken by Insight-HXMT is just after the transition to the weak propeller stage. Our analysis shows the necessity of a comptonization component to take into account the existence of an electron cloud resulting photons partly up-scattered., Comment: 8 pages, 4 figures, accepted for publication in JHEAp

GRO J1008-57, as a Be/X-ray transient pulsar, is considered to have the highest magnetic field in known neutron star X-ray binary systems. Observational data of the X-ray outbursts in GRO J1008-57 from 2017 to 2020 were collected by the Insight-HXMT satellite. In this work, the spin period of the neutron star in GRO J1008-57 was determined to be about 93.28 seconds in August 2017, 93.22 seconds in February 2018, 93.25 seconds in June 2019 and 93.14 seconds in June 2020. GRO J1008-57 evolved in the spin-up process with a mean rate of $-(2.10\pm 0.05)\times$10$^{-4}$ s/d from 2009 -- 2018, and turned into a spin down process with a rate of $(6.7\pm 0.6)\times$10$^{-5}$ s/d from Feb 2018 to June 2019. During the type II outburst of 2020, GRO J1008-57 had the spin-up torque again. During the torque reversals, the pulse profiles and continuum X-ray spectra did not change significantly, and the cyclotron resonant scattering feature around 80 keV was only detected during the outbursts in 2017 and 2020. Based on the observed mean spin-up rate, we estimated the inner accretion disk radius in GRO J1008-57 (about 1 - 2 times of the Alfv\'{e}n radius) by comparing different accretion torque models of magnetic neutron stars. During the spin-down process, the magnetic torque should dominate over the matter accreting inflow torque, and we constrained the surface dipole magnetic field $B\geq 6\times 10^{12}$ G for the neutron star in GRO J1008-57, which is consistent with the magnetic field strength obtained by cyclotron line centroid energy., Comment: 10 pages, 5 figures, and 3 tables, the Journal of High Energy Astrophysics in press

In this paper, we presented a detailed timing analysis of a prominent outburst of 4U 0115+63 detected by \textit{Insight}-HXMT in 2017 August. The spin period of the neutron star was determined to be $3.61398\pm 0.00002$ s at MJD 57978. We measured the period variability and extract the orbital elements of the binary system. The angle of periastron evolved with a rate of $0.048\pm0.003$ $yr^{-1}$. The light curves are folded to sketch the pulse profiles in different energy ranges. A multi-peak structure in 1-10 keV is clearly illustrated. We introduced wavelet analysis into our data analysis procedures to study QPO signals and perform a detailed wavelet analysis in many different energy ranges. Through the wavelet spectra, we report the discovery of a QPO at the frequency $\sim 10$ mHz. In addition, the X-ray light curves showed multiple QPOs in the period of $\sim 16-32 $ s and $\sim 67- 200 $ s. We found that the $\sim100$ s QPO was significant in most of the observations and energies. There exist positive relations between X-ray luminosity and their Q-factors and S-factors, while the QPO periods have no correlation with X-ray luminosity. In wavelet phase maps, we found that the pulse phase of $\sim 67- 200 $ s QPO drifting frequently while the $\sim 16-32 $ s QPO scarcely drifting. The dissipation of oscillations from high energy to low energy was also observed. These features of QPOs in 4U 0115+63 provide new challenge to our understanding of their physical origins., Comment: 14 pages, 9 figures, and 6 tables. This work has been submitted to MNRAS after the referee's report

Carbon nanotubes tend to collapse when their diameters exceed a certain threshold, or when a sufficiently large external pressure is applied on their walls. Their radial stability of tubes has been studied in each of these cases, however a general theory able to predict collapse is still lacking. Here, we propose a simple model predicting stability limits as a function of the tube diameter, the number of walls and the pressure. The model is supported by atomistic simulations, experiments, and is used to plot collapse phase diagrams. We have identified the most stable carbon nanotube, which can support a maximum pressure of 18 GPa before collapsing. The latter was identified as a multiwall tube with an internal tube diameter of 12nm and 30 walls. This maximum pressure is lowered depending on the internal tube diameter and the number of walls. We then identify a tube diameter domain in which the radial mechanical stability can be treated as equivalent to macroscopic tubes, known to be described by the canonical L\'evy-Carrier law. This multiscale behavior is shown to be in good agreement with experiments based on O-ring gaskets collapse, proposed as a simple macroscopic parallel to nanotubes in this domain., Comment: 32 pages, 9 figures

We report on the Insight-HXMT observations of the new black hole X-ray binary MAXI J1820+070 during its 2018 outburst. Detailed spectral analysis via the continuum fitting method shows an evolution of the inferred spin during its high soft sate. Moreover, the hardness ratio, the non-thermal luminosity and the reflection fraction also undergo an evolution, exactly coincident to the period when the inferred spin transition takes place. The unphysical evolution of the spin is attributed to the evolution of the inner disc, which is caused by the collapse of a hot corona due to condensation mechanism or may be related to the deceleration of a jet-like corona. The studies of the inner disc radius and the relation between the disc luminosity and the inner disc radius suggest that, only at a particular epoch, did the inner edge of the disc reach the innermost stable circular orbit and the spin measurement is reliable. We then constrain the spin of MAXI J1820+070 to be a*=0.2^{+0.2}_{-0.3}. Such a slowly spinning black hole possessing a strong jet suggests that its jet activity is driven mainly by the accretion disc rather than by the black hole spin., Comment: 14 pages, 13 figures, 5 tables, accepted for publication in MNRAS

We investigate the absorption and emission features in observations of GX 301-2 detected with Insight-HXMT/LE in 2017-2019. At different orbital phases, we found prominent Fe Kalpha, Kbeta and Ni Kalpha lines, as well as Compton shoulders and Fe K-shell absorption edges. These features are due to the X-ray reprocessing caused by the interaction between the radiation from the source and surrounding accretion material. According to the ratio of iron lines Kalpha and Kbeta, we infer the accretion material is in a low ionisation state. We find an orbital-dependent local absorption column density, which has a large value and strong variability around the periastron. We explain its variability as a result of inhomogeneities of the accretion environment and/or instabilities of accretion processes. In addition, the variable local column density is correlated with the equivalent width of the iron Kalpha lines throughout the orbit, which suggests that the accretion material near the neutron star is spherically distributed., Comment: 10 pages, 5 figures, 2 tables, accepted for publication in MNRAS

We present the analysis of the brightest flare that was recorded in the \emph{Insight}-HMXT data set, in a broad energy range (2$-$200 keV) from the microquasar GRS~1915+105 during an unusual low-luminosity state. This flare was detected by \emph{Insight}-HXMT among a series of flares during 2 June 2019 UTC 16:37:06 to 20:11:36, with a 2-200 keV luminosity of 3.4$-$7.27$\times10^{38}$ erg s$^{-1}$. Basing on the broad-band spectral analysis, we find that the flare spectrum shows different behaviors during bright and faint epochs. The spectrum of the flare can be fitted with a model dominated by a power-law component. Additional components show up in the bright epoch with a hard tail and in the faint epoch with an absorption line $\sim$ 6.78 keV. The reflection component of the latter is consistent with an inner disk radius $\sim$ 5 times larger than that of the former. These results on the giant flare during the "unusual" low-luminosity state of GRS~1915+105 may suggest that the source experiences a possible fast transition from a jet-dominated state to a wind-dominated state. We speculate that the evolving accretion disk and the large-scale magnetic field may play important roles in this peculiar huge flare.

We present a detailed timing study of the brightest persistent X-ray source Sco X-1 using the data collected by the Hard X-ray Modulation Telescope ($Insight$-HXMT) from July 2017 to August 2018. A complete $Z$-track hardness-intensity diagram (HID) is obtained. The normal branch oscillations (NBOs) at $\sim$ 6 Hz in the lower part of the normal branch (NB) and the flare branch oscillations (FBOs) at $\sim$ 16 Hz in the beginning part of the flaring branch (FB) are found in observations with the Low Energy X-ray Telescope (LE) and the Medium Energy X-ray Telescope (ME) of $Insight$-HXMT, while the horizontal branch oscillations (HBOs) at $\sim$ 40 Hz and the kilohertz quasi-periodic oscillations (kHz QPOs) at $\sim$ 800 Hz are found simultaneously up to 60 keV for the first time on the horizontal branch (HB) by the High Energy X-ray Telescope (HE) and ME. We find that for all types of the observed QPOs, the centroid frequencies are independent of energy, while the root mean square (rms) increases with energy; the centroid frequencies of both the HBOs and kHz QPOs increase along the $Z$-track from the top to the bottom of the HB; and the NBOs show soft phase lags increasing with energy. A continuous QPO transition from the FB to NB in $\sim$ 200 s are also detected. Our results indicate that the non-thermal emission is the origin of all types of QPOs, the innermost region of the accretion disk is non-thermal in nature, and the corona is nonhomogeneous geometrically.

The evidences for the influence of thermonuclear (type-I) X-ray bursts upon the surrounding environments in neutron star low-mass X-ray binaries (LMXB) were detected previously via spectral and timing analyses. Benefitting from a broad energy coverage of Insight-HXMT, we analyze one photospheric radius expansion (PRE) burst, and find an emission excess at soft X-rays. Our spectral analysis shows that, such an excess is not likely relevant to the disk reflection induced by the burst emission and can be attributed to an enhanced pre-burst/persistent emission. We find that the burst and enhanced persistent emissions sum up to exceed Eddington luminosity by $\sim$ 40 percentages. We speculate that the enhanced emission is from a region beyond the PRE radius, or through the Comptonization of the corona., Comment: accepted by JHEA(Journal of High Energy Astrophysics)

Corona cooling was detected previously from stacking a series of short type-I bursts occurred during the low/had state of atoll outburst. Type-I bursts are hence regarded as sharp probe to our better understanding on the basic property of the corona. The launch of the first Chinese X-ray satellite Insight-HXMT has large detection area at hard X-rays which provide almost unique chance to move further in this research field. We report the first detection of the corona cooling by Insight-HXMT from single short type-I burst showing up during {\bf flare} of 4U 1636-536. This type-I X-ray burst has a duration of $\sim$13 seconds and hard X-ray shortage is detected with significance 6.2~$\sigma$ in 40-70 keV. A cross-correlation analysis between the lightcurves of soft and hard X-ray band, shows that the corona shortage lag the burst emission by 1.6 $\pm$1.2~s. These results are consistent with those derived previously from stacking a large amount of bursts detected by RXTE/PCA within a series of {\bf flares} of 4U 1636-536. Moreover, the broad bandwidth of Insight-HXMT allows as well for the first time to infer the burst influence upon the continuum spectrum via performing the spectral fitting of the burst, which ends up with the finding that hard X-ray shortage appears at around 40 keV in the continuum spectrum. These results suggest that the evolution of the corona along with the outburst{\bf /flare} of NS XRB may be traced via looking into a series of embedded type-I bursts by using Insight-HXMT., Comment: published in 2018, ApJL,864, L30

We present a study of timing properties of the accreting pulsar 2S 1417-624 observed during its 2018 outburst, based on Swift/BAT, Fermi/GBM, Insight-HXMT and NICER observations. We report a dramatic change of the pulse profiles with luminosity. The morphology of the profile in the range 0.2-10.0keV switches from double to triple peaks at $\sim2.5$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$ and from triple to quadruple peaks at $\sim7$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$. The profile at high energies (25-100keV) shows significant evolutions as well. We explain this phenomenon according to existing theoretical models. We argue that the first change is related to the transition from the sub to the super-critical accretion regime, while the second to the transition of the accretion disc from the gas-dominated to the radiation pressure-dominated state. Considering the spin-up as well due to the accretion torque, this interpretation allows to estimate the magnetic field self-consistently at $\sim7\times 10^{12}$G., Comment: 7 pages, 4 figures, 1 tables, accepted for publication in MNRAS

The long-term evolution of the centroid energy of the CRSF in Her X-1 is still a mystery. We report a new measurement from a campaign between {\sl Insight}-HXMT and {\sl NuSTAR} performed in February 2018. Generally, the two satellites show well consistent results of timing and spectral properties. The joint spectral analysis confirms that the previously observed long decay phase has ended, and that the line energy instead keeps constant around 37.5 keV after flux correction.

We report on analysis of observations of the bright transient X-ray pulsar \src obtained during its 2017-2018 giant outburst with Insight-HXMT, \emph{NuSTAR}, and \textit{Swift} observatories. We focus on the discovery of a sharp state transition of the timing and spectral properties of the source at super-Eddington accretion rates, which we associate with the transition of the accretion disk to a radiation pressure dominated (RPD) state, the first ever directly observed for magnetized neutron star. This transition occurs at slightly higher luminosity compared to already reported transition of the source from sub- to super-critical accretion regime associate with onset of an accretion column. We argue that this scenario can only be realized for comparatively weakly magnetized neutron star, not dissimilar to other ultra-luminous X-ray pulsars (ULPs), which accrete at similar rates. Further evidence for this conclusion is provided by the non-detection of the transition to the propeller state in quiescence which strongly implies compact magnetosphere and thus rules out magnetar-like fields., Comment: Submitted to MNRAS