Search

GRB221009A is the brightest gamma-ray burst ever detected. To probe the very-high-energy (VHE, $>$\!100 GeV) emission, the High Energy Stereoscopic System (H.E.S.S.) began observations 53 hours after the triggering event, when the brightness of the moonlight no longer precluded observations. We derive differential and integral upper limits using H.E.S.S. data from the third, fourth, and ninth nights after the initial GRB detection, after applying atmospheric corrections. The combined observations yield an integral energy flux upper limit of $\Phi_\mathrm{UL}^{95\%} = 9.7 \times 10^{-12}~\mathrm{erg\,cm^{-2}\,s^{-1}}$ above $E_\mathrm{thr} = 650$ GeV. The constraints derived from the H.E.S.S. observations complement the available multiwavelength data. The radio to X-ray data are consistent with synchrotron emission from a single electron population, with the peak in the SED occurring above the X-ray band. Compared to the VHE-bright GRB190829A, the upper limits for GRB221009A imply a smaller gamma-ray to X-ray flux ratio in the afterglow. Even in the absence of a detection, the H.E.S.S. upper limits thus contribute to the multiwavelength picture of GRB221009A, effectively ruling out an IC dominated scenario., Comment: 10 pages, 4 figures. Accepted for publication in APJL. Corresponding authors: J. Damascene Mbarubucyeye, H. Ashkar, S. J. Zhu, B. Reville, F. Sch\"ussler

In this paper, we present the results of the second census of pulsars in the decametre wave range at the Ukrainian T-shaped Radio telescope, second modification (UTR-2). Over the past 10 yr, the worldwide number of discovered nearby pulsars has doubled, which means that there is some urgency to search for a low-frequency radio emission from these newly discovered sources. To increase the sensitivity of this census, the integration time was doubled compared with the first census of 2010–2013. As a result, the decametre radio emission of 20 pulsars was detected, and their flux densities and the shape of their pulses were obtained. The dispersion measures for 10 pulsars and the rotation periods for eight pulsars were refined. For several pulsars, the scattering time constant and FWHM were estimated in the decametre wave range. The upper limits of flux densities for 102 pulsars that have not yet been detected were also estimated.

Energetic particle precipitation associated with pulsating aurora (PsA) can reach down to lower mesospheric altitudes and deplete ozone. It is well documented that pulsating aurora is a common phenomenon during substorm recovery phases. This indicates that using magnetic indices to model the chemistry induced by PsA electrons could underestimate the energy deposition in the atmosphere. Integrating satellite measurements of precipitating electrons in models is considered to be an alternative way to account for such an underestimation. One way to do this is to test and validate the existing ion chemistry models using integrated measurements from satellite and ground-based observations. By using satellite measurements, an average or typical spectrum of PsA electrons can be constructed and used as an input in models to study the effects of the energetic electrons in the atmosphere. In this study, we compare electron densities from the EISCAT (European Incoherent Scatter scientific radar system) radars with auroral ion chemistry and the energetics model by using pulsating aurora spectra derived from the Polar Operational Environmental Satellite (POES) as an energy input for the model. We found a good agreement between the model and EISCAT electron densities in the region dominated by patchy pulsating aurora. However, the magnitude of the observed electron densities suggests a significant difference in the flux of precipitating electrons for different pulsating aurora types (structures) observed.

Solar eclipses provide an excellent opportunity to study the effects of a sudden localized change in photoionization flux in the Earth's ionosphere and its consequent repercussion in the Geomagnetic field. We have focused on a subset of the data available from the North American 2017 eclipse in order to study VTEC measurements from GNSS data and geomagnetic field estimations from INTERMAGNET observatories near the eclipse path. Our simultaneous analysis of both datasets allowed us to quantify the ionosphere and magnetic field reaction to the eclipse event with which allowed us to compare how differently these take place in time. We found that studying the behaviour of VTEC differences with respect to reference values provides better insight of the actual eclipse effect and were able to characterize the dependence of parameters such as time delay of maximum depletion and recovery phase. We were also able to test models that link the ionospheric variations in a quantitative manner. Total electron content depletion measured from GNSS were fed into an approximation of Ashour-Chapman model at the locations of geomagnetic observatories and its predictions match the behaviour of magnetic field components in time and magnitude strikingly accurately., Comment: Accepted for publication in Advances in Space Research

The astrophysical journal / 2 944(2), L34 (2023). doi:10.3847/2041-8213/acb6d7, We present a multi-messenger model for the prompt emission from GRB 221009A within the internal shock scenario. We consider the time-dependent evolution of the outflow with its impact on the observed light curve from multiple collisions, as well as the self-consistent generation of the electromagnetic spectrum in synchrotron and inverse Compton-dominated scenarios. Our lepto-hadronic model includes UHE protons potentially accelerated in the outflow, and their feedback on spectral energy distribution and on the neutrino emission. We find that we can roughly reproduce the observed light curves with an engine with varying ejection velocity of ultrarelativistic material, which has an intermediate quiescent period of about 200 s and a variability timescale of ∼1 s. We consider baryonic loadings of 3 and 30 that are compatible with the hypothesis that the highest-energetic LHAASO photons might come from UHECR interactions with the extragalactic background light, and the paradigm that energetic GRBs may power the UHECR flux. For these values and the high dissipation radii considered, we find consistency with the nonobservation of neutrinos and no significant signatures on the electromagnetic spectrum. Inverse Compton-dominated scenarios from the prompt emission are demonstrated to lead to about an order of magnitude higher fluxes in the HE range; this enhancement is testable via its spectral impact in the Fermi-GBM and LAT ranges., Published by Institute of Physics Publ., London

A deep survey of the Large Magellanic Cloud at ~0.1-100TeV photon energies with the Cherenkov Telescope Array is planned. We assess the detection prospects based on a model for the emission of the galaxy, comprising the four known TeV emitters, mock populations of sources, and interstellar emission on galactic scales. We also assess the detectability of 30 Doradus and SN 1987A, and the constraints that can be derived on the nature of dark matter. The survey will allow for fine spectral studies of N157B, N132D, LMC P3, and 30 Doradus C, and half a dozen other sources should be revealed, mainly pulsar-powered objects. The remnant from SN 1987A could be detected if it produces cosmic-ray nuclei with a flat power-law spectrum at high energies, or with a steeper index 2.3-2.4 pending a flux increase by a factor >3-4 over ~2015-2035. Large-scale interstellar emission remains mostly out of reach of the survey if its >10GeV spectrum has a soft photon index ~2.7, but degree-scale 0.1-10TeV pion-decay emission could be detected if the cosmic-ray spectrum hardens above >100GeV. The 30 Doradus star-forming region is detectable if acceleration efficiency is on the order of 1-10% of the mechanical luminosity and diffusion is suppressed by two orders of magnitude within, Accepted for publication in MNRAS. Corresponding authors: Pierrick Martin, Maria Isabel Bernardos Martin, Fabio Iocco

We present a search for transient radio sources on timescales of 2 to 9 yr at 150 MHz. This search is conducted by comparing the first Alternative Data Release of the TIFR GMRT Sky Survey (TGSS ADR1) and the second data release of the LOFAR Two-metre Sky Survey (LoTSS DR2). The overlapping survey area covers 5570 $\rm{deg}^2$ on the sky, or 14 per cent of the total sky. We introduce a method to compare the source catalogues that involves a pair match of sources, a flux density cutoff to meet the survey completeness limit and a newly developed compactness criterion. This method is used to identify both transient candidates in the TGSS source catalogue that have no counterpart in the LoTSS catalogue and transient candidates in LoTSS without a counterpart in TGSS. We find that imaging artefacts and uncertainties and variations in the flux density scales complicate the transient search. Our method to search for transients by comparing two different surveys, while taking into account imaging artefacts around bright sources and misaligned flux scales between surveys, is universally applicable to future radio transient searches. No transient sources were identified, but we are able to place an upper limit on the transient surface density of $, 14 pages, 11 figures

Direct collapse black holes~(BH) are promising candidates for producing massive $z\gtrsim 6$ quasars, but their formation requires fine-tuned conditions. In this work, we use cosmological zoom simulations to study systematically the impact of requiring: 1) low gas angular momentum, and 2) a minimum incident Lyman-Werner~(LW) flux in order to form BH seeds. We probe the formation of seeds (with initial masses of $M_{\rm seed} \sim 10^4$ - $10^6 M_{\odot}/h)$ in halos with a total mass $> 3000\times M_{\mathrm{seed}}$ and a dense, metal poor gas mass $> 5\times M_{\mathrm{seed}}$. We find that the seed-forming halos have a prior history of star formation and metal enrichment, but contain pockets of dense, metal poor gas. When seeding is further restricted to halos with low gas spins, the number of seeds formed is suppressed by factors of $\sim6$ compared to the baseline model, regardless of the seed mass. Seed formation is much more strongly impacted if the dense, metal poor gas is required to have a critical LW flux ($J_{\mathrm{crit}}$). Even for $J_{\mathrm{crit}}$ values as low as $50J_{21}$, no $8\times10^{5}M_{\odot}/h$ seeds are formed. While lower mass ($1.25\times10^{4},1\times10^{5} M_{\odot}/h$) seeds do form, they are strongly suppressed~(by factors of $\sim10-100$) compared to the baseline model at gas mass resolutions of $\sim10^4~M_{\odot}/h$ (with even stronger suppression at higher resolutions). As a result, BH merger rates are also similarly suppressed. Since early BH growth is dominated by mergers in our models, no seeds are able to grow to the supermassive regime~($\gtrsim10^6 M_{\odot}/h$) by $z=7$. Our results hint that producing the bulk of the $z\gtrsim6$ supermassive BH population may require alternate seeding scenarios that do not depend on the LW flux, early BH growth dominated by rapid or super-Eddington accretion, or a combination of these possibilities., 21 pages, 13 figures

We use the Sherwood-Relics suite of hybrid hydrodynamical and radiative transfer simulations to model the effect of inhomogeneous reionisation on the 1D power spectrum of the \Lya forest transmitted flux at redshifts $4.2\leq z \leq 5$. Relative to models that assume a homogeneous UV background, reionisation suppresses the power spectrum at small scales, $k \sim 0.1\rm\,km^{-1}\,s$ by $\sim 10$ per cent because of spatial variations in the thermal broadening kernel and the divergent peculiar velocity field associated with over-pressurised intergalactic gas. On larger scales, $k4$, $0.006 \rm \,km^{-1}\,s\leq k \leq 0.2 \rm\, km^{-1}\,s$, we find inhomogeneous reionisation does not introduce any significant bias in thermal parameter recovery for the current measurement uncertainties of $\sim 10$ per cent. However, for $5$ per cent uncertainties, $\sim 1\sigma$ shifts between the estimated and true parameters occur., Comment: Accepted by MNRAS (21 November 2021), 20 pages, 11 figures

Sagittarius A* exhibits regular variability in its multiwavelength emission, including daily X-ray flares and roughly continuous near-infrared (NIR) flickering. The origin of this variability is still ambiguous since both inverse Compton and synchrotron emission are possible radiative mechanisms. The underlying particle distributions are also not well constrained, particularly the non-thermal contribution. In this work, we employ the GPU-accelerated general relativistic magnetohydrodynamics (GRMHD) code H-AMR perform a study of flare flux distributions, including the effect of particle acceleration for the first time in high-resolution 3D simulations of Sgr A*. For the particle acceleration, we use the general relativistic ray-tracing (GRRT) code BHOSS to perform the radiative transfer, assuming a hybrid thermal+non-thermal electron energy distribution. We extract ~60 h lightcurves in the sub-millimetre, NIR and X-ray wavebands, and compare the power spectra and the cumulative flux distributions of the lightcurves to statistical descriptions for Sgr A* flares. Our results indicate that non-thermal populations of electrons arising from turbulence-driven reconnection in weakly magnetised accretion flows lead to moderate NIR and X-ray flares and reasonably describe the X-ray flux distribution while fulfilling multiwavelength flux constraints. These models exhibit high rms% amplitudes, >~150% both in the NIR and the X-rays, with changes in the accretion rate driving the 230~GHz flux variability, in agreement with Sgr A* observations., Accepted by MNRAS

Neutron stars in low mass X-ray binaries are hypothesised to emit continuous gravitational waves that may be detectable by ground-based observatories. The torque balance model predicts that a higher accretion rate produces larger-amplitude gravitational waves, hence low mass X-ray binaries with high X-ray flux are promising targets for gravitational wave searches. The detection of X-ray pulsations would identify the spin frequency of these neutron stars, and thereby improve the sensitivity of continuous gravitational-wave searches by reducing the volume of the search parameter space. We perform a semi-coherent search for pulsations in the two low mass X-ray binaries Scorpius X-1 and Cygnus X-2 using X-ray data from the \textit{ Rossi X-ray Timing Explorer} Proportional Counter Array. We find no clear evidence for pulsations, and obtain upper limits (at $90\%$ confidence) on the fractional pulse amplitude, with the most stringent being $0.034\%$ for Scorpius X-1 and $0.23\%$ for Cygnus X-2. These upper limits improve upon those of Vaughan et al. (1994) by factors of $\sim 8.2$ and $\sim 1.6$ respectively., 10 pages, 11 figures

Solar activity modulates cosmic ray (CR) flux with different magnitudes on different time scales. Several investigations from locations with different, mainly low, cutoff rigidity utilizing different solar activity parameters have been carried out to characterize their influence on the observed CR flux. In this study, the effect of seven solar, interplanetary, and geophysical parameters on the secondary CR muons has been investigated and established using correlation analyses. Daily mean values of the CR data for the period between 2002 and 2012 were obtained from the King Abdulaziz City for Science and Technology (KACST) muon detector, Riyadh, central Saudi Arabia (Rc = 14.4 GV). Considered parameters are sunspot number, radial component of the interplanetary magnetic field, Kp index, solar radio emission flux at 10.7 cm, Dst index, solar wind speed, and solar wind density. Same analyses have been carried out using data from the Oulu neutron monitor (Rc = 0.8 GV) for comparison purposes. Obtained results showed that the intensities of the secondary CRs from both stations are significantly correlated with the Dst index and plasma density, and anti-correlated with the rest of the variables. These results are in agreement with those obtained by several investigators. Magnitude and strength of the correlations between the considered variables and CRs were higher and stronger for the CR neutrons than the CR muons. These findings may be due to the greater influence of seasonal variations on the CR muons than on the CR neutrons. Additionally, the muon detector is sensitive (due to its high cutoff rigidity) to higher-energy CR particles, which are less affected by solar disturbances. Time-lag cross-correlation analyses between the monthly mean CR values from both stations and the monthly mean values of the considered variables were conducted. Comparable results with the previous investigations were found. Multivariable models using the seven parameters were developed to predict the CR variations for both sites. For the data from the KACST detector, the model was able to predict the measured data with a correlation coefficient of 0.48 and a standard deviation of 0.60%. On the other hand, the developed model for the Oulu neutron monitor has a correlations coefficient of 0.77 and a standard deviation of 3.7%.

The Diffuse Ionized Gas (DIG) contributes to the nebular emission of galaxies, resulting in emission line flux ratios that can be significantly different from those produced by HII regions. Comparing the emission of [SII]6717,31 between pointed observations of HII regions in nearby galaxies and integrated spectra of more distant galaxies, it has been recently claimed that the DIG can also deeply affect the emission of bright, star-forming galaxies, and that a large correction must be applied to observed line ratios to recover the genuine contribution from HII regions. Here we show instead that the effect of DIG on the integrated spectra of star-forming galaxies is lower than assumed in previous work. Indeed, aperture effects on the spectroscopy of nearby HII regions are largely responsible for the observed difference: when spectra of local HII regions are extracted using large enough apertures while still avoiding the DIG, the observed line ratios are the same as in more distant galaxies. This result is highly relevant for the use of strong-line methods to measure metallicity., Comment: MNRAS, in press

In this work, we investigate upper limits on the global escape fraction of ionizing photons ($f_{\rm esc/global}^{\rm abs}$) from a sample of galaxies probed for Lyman-continuum (LyC) emission characterized as non-LyC and LyC leakers. We present a sample of nine clean non-contaminated (by low-redshift interlopers, CCD problems, and internal reflections of the instrument) galaxies that do not show significant (>3σ) LyC flux in the range 880

BL Lacertae (BL Lac) is categorized as TeV blazar and considered as a possible source of astrophysical neutrinos. In 2020, the brightest X-ray flare ever detected from it. A detailed study can answer many puzzling questions related to multiband emissions and fast-flux variability often seen in this kind of source. We have performed the temporal and spectral analysis of the brightest flare. The variability is characterized by the fractional variability amplitude and the variability time. We found that the source has crossed all its previous limits of flux and reached the maximum ever seen from it in optical and X-rays. It is highly variable in X-rays with fractional variability above 100$\%$ (1.8397$\pm$0.0181) and the fastest variability time of 11.28 hours within a day. The broadband light curves correlation with X-ray suggest a time lag of one day. A broadband SED modeling is pursued to understand the possible physical mechanisms responsible for broadband emission. Modeling requires two emission regions located at two different sites to explain the low and high flux states. A significant spectral change is observed in the optical-UV and X-ray spectrum during the high state, which eventually leads to shifts in the location of the synchrotron peak towards higher energy, suggesting an emergence of a new HBL component., Comment: 13 pages, 11 figures, 3 tables, Accepted in MNRAS

Modern-Era Retrospective Analysis for Research and Applications, Version-2 (MERRA-2) total column ozone (TCO) data from 1980 to 2018 and ground based TCO data were used to investigate its responses to periodic variations concerning stratospheric dynamics such as tropopause heights and Quasi-Biennial Oscillations (QBO), and solar activity signatures of Magnesium II core-to-wing ratio (Mg II) and F10.7 cm solar radio flux (F10.7) observed over Mountain (Mt.) Waliguan (36.17°N, 100.53°E), China. Cross-correlation and wavelet analyses were used for carrying out the investigation. The MgII and F10.7 show a periodicity of around an 11-year solar cycle (SC) while the TCO trend is though periodic and it does not exhibit one to one correspondence with the 11-year SC. The highest value in TCO is 355.7 DU noted in February (1982) whereas, the lowest value of TCO is 256.02 DU in October (1988) during the entire study period. The average value in TCO is 291.17 ± 0.94 DU with an amplitude of variation of 99.68 DU (∼ 28%). The trend of TCO depicts a persistent distinct seasonal pattern with maxima in winter/spring (December through April) and minima in summer/autumn (August through November). The long-term variability shows a notable decreasing trend in TCO by ∼ 30%, 13%, and 5% per decade for SCs 21-23, where SC 24 represents an increasing trend of ∼ 1% per decade. Hence, this illustrates the significant recovery of TCO over Mt. Waliguan, China which is as a result of the reduction of the ozone depletion agents arising as a positive outcome of the implementation of the Montreal regulatory policy. From wavelet analysis, the F10.7, MgII, tropopause height, QBO and TCO show a significant dominant periodicity of 10-12 months. It is concluded that the variations in mean TCO are strongly correlated with variations of the F10.7 cm and Mg-II solar indices, and the mean TCO responds rapidly to solar activity variations.

Several studies have been showing that cosmic ray intensity variations observed at the Earth surface can be a useful tool to analyze and predict space weather and climate phenomena. However, before reach the ground level, the cosmic ray particles interact with different atmospheric atoms. This results in the generation of different secondary particles that are differently affected by atmospheric effects. For example, data from neutron monitors, which are widely used to study cosmic ray flux changes related to space phenomena, are significantly affected by atmospheric pressure changes. Thus, for a more detailed study of these instruments’ data, it is necessary to know very well how works the pressure affects secondary neutrons. In this way, we seek to understand how this effect changes according to the cosmic rays’ energy. Thus, in this work, we compare how the barometric coefficient ( β ), which quantifies the pressure effect on the cosmic ray intensity, changes with the geomagnetic cutoff rigidity ( R C ), which quantifies the minimal rigidity/energy that a primary particle needs to reach a given location in a chosen direction. To do that, we used data recorded between January 2008 to December 2016 by several neutron monitors around the globe with R C from 0.1 to 8.3 GV. We experimentally obtained β in two ways: (I) the typical one, which does not explicitly consider isotropic variations unrelated to the pressure effect; and (II) one that considers such variations by using coupling coefficients and a reference neutron monitor station. Like other works, we found that the pressure effect seems to present a logarithmic relation with the cutoff rigidity when we estimate β through the typical method used. On the other hand, we found that the pressure effect linearly decreases with the cutoff rigidity when considering an adjustment in this typical method to eliminate external isotropic variations. This linear relation appears to not change accordingly to the reference monitor chosen. These results enhance our knowledge of how the pressure effect globally acts on the cosmic ray intensity observed at ground level and, consequently, improves future analysis about global changes in the cosmic ray flux related to space weather phenomena.

The astrophysical journal 938(1), 38 (2022). doi:10.3847/1538-4357/ac8de4, The majority of astrophysical neutrinos have undetermined origins. The IceCube Neutrino Observatory has observed astrophysical neutrinos but has not yet identified their sources. Blazars are promising source candidates, but previous searches for neutrino emission from populations of blazars detected in ≳GeV gamma rays have not observed any significant neutrino excess. Recent findings in multimessenger astronomy indicate that high-energy photons, coproduced with high-energy neutrinos, are likely to be absorbed and reemitted at lower energies. Thus, lower-energy photons may be better indicators of TeV–PeV neutrino production. This paper presents the first time-integrated stacking search for astrophysical neutrino emission from MeV-detected blazars in the first Fermi Large Area Telescope low energy (1FLE) catalog using ten years of IceCube muon–neutrino data. The results of this analysis are found to be consistent with a background-only hypothesis. Assuming an E$^{−2}$ neutrino spectrum and proportionality between the blazars MeV gamma-ray fluxes and TeV–PeV neutrino flux, the upper limit on the 1FLE blazar energy-scaled neutrino flux is determined to be 1.64 × 10$^{−12}$ TeV cm$^{−2}$ s$^{−1}$ at 90% confidence level. This upper limit is approximately 1% of IceCube’s diffuse muon–neutrino flux measurement., Published by Univ., Chicago, Ill. [u.a.]

The origin and evolution of the 1/f power law observed in the energy spectrum of solar coronal and solar wind fluctuations at scales of around an hour is not entirely understood. Several existing theories aim at explaining it, involving both linear and nonlinear mechanisms. An often overlooked property of the solar corona and solar wind is their highly inhomogeneous nature. In this paper we investigate the linear evolution of pure Alfv\'en and surface Alfv\'en waves propagating through a plasma which is inhomogeneous across the magnetic field. The inhomogeneity is given by density, which we model to be two-dimensional colored noise, with power spectral slopes ranging from -2 to -1. Alfv\'en waves propagate independently on individual magnetic field lines, and eventually get completely out of phase through the process of phase mixing, leading to unrealistic spectra. When the coupling between the inhomogeneous background and the propagating waves is fully accounted for, transverse waves such as surface Alfv\'en waves (also referred to as kink or Alfv\'enic) appear, showing collective wave behavior of neighboring magnetic field lines with different Alfv\'en speeds. We show that the linear cascade of surface Alfv\'en wave energy, induced by phase mixing and resonant absorption, leads to a perpendicular wave energy spectrum which tends to the perpendicular power spectrum of the background density. Based on our model, we propose that a perpendicular density power spectrum of 1/f in the solar corona can induce, through linear processes, the 1/f spectrum of the fluctuations that is observed at the largest scales., Comment: Accepted in ApJ

Soft Ultra-Luminous X-ray (ULXs) sources are a subclass of the ULXs that can switch from a supersoft spectral state, where most of the luminosity is emitted below 1 keV, to a soft spectral state with significant emission above 1 keV. In a few systems, dips have been observed. The mechanism behind this state transition and the dips nature are still debated. To investigate these issues, we obtained a long XMM-Newton monitoring campaign of a member of this class, NGC 247 ULX-1. We computed the hardness-intensity diagram for the whole dataset and identified two different branches: the normal branch and the dipping branch, which we study with four and three hardness-intensity resolved spectra, respectively. All seven spectra are well described by two thermal components: a colder ($kT_{\rm bb}$ $\sim$ 0.1-0.2 keV) black-body, interpreted as emission from the photo-sphere of a radiatively-driven wind, and a hotter ($kT_{\rm disk}$ $\sim$ 0.6 keV) multicolour disk black-body, likely due to reprocessing of radiation emitted from the innermost regions. In addition, a complex pattern of emission and absorption lines has been taken into account based on previous high-resolution spectroscopic results. We studied the evolution of spectral parameters and the flux of the two thermal components along the two branches and discuss two scenarios possibly connecting the state transition and the dipping phenomenon. One is based on geometrical occultation of the emitting regions, the other invokes the onset of a propeller effect., 13 pages, accepted by MNRAS on 2021/08/16

Blazar variability can be described as flaring activities on a wide range of time-scales over a baseline flux level. It is important to detect and distinguish baseline flux changes from long-term flare variations, since the former may reflect state transitions caused by a secular change in bulk properties such as the bulk Lorentz factor or the viewing angle. We report such a transition observed in the 11 yr Fermi γ-ray light curve of the blazar 3C 66A, where the baseline flux of ∼1.8 × 10−7 ph cm−2 s−1 persisted for three years and then changed over a month time-scale to ∼0.8 × 10−7 ph cm−2 s−1 and remained in that level for the next eight years. Moreover, there is evidence for a similar shift in baseline flux in the optical band. Broad-band spectral energy distribution modelling for the two activity states reveals that the baseline flux change is consistent with an overall decrease in the Doppler factor, which in turn implies an increase in the jet inclination angle by ∼1° or a decrease in the bulk Lorentz factor by ∼25 per cent. We discuss the implication of such a variation occurring on a month time-scale.

We present a detailed spectral and timing analysis of Cygnus X-1 with multi-epoch observations, during $2016$ to $2019$, by SXT and LAXPC on-board AstroSat. We model the spectra in broad energy range of $0.5\!-\!70.0\,\rm{keV}$ to study the evolution of spectral properties while Cygnus X-1 transited from hard state to an extreme soft state via intermediate states in 2017. Simultaneous timing features are also examined by modelling the power density spectra in $3.0\!-\!50.0\,\rm{keV}$ . We find that during high-soft state observations, made by AstroSat on Oct $24,\,2017$ (MJD $58050$), the energy spectrum of the source exhibits an inner disk temperature (kT$\rm_{in}$) of $0.46\!\pm\!0.01\,\rm{keV}$ , a very steep photon index ($\Gamma$) of $3.15\!\pm\!0.03$ along with a fractional disk flux contribution of $\sim\!45\%$. The power density spectrum in the range of $0.006\!-\!50.0\,\rm{Hz}$ is also very steep with a power-law index of $1.12\!\pm\!0.04$ along with a high RMS value of $\sim\!25\%$. Comparing the spectral softness of high-soft state with those of previously reported, we confirm that {\it AstroSat} observed Cygnus X-1 in the `softest' state. The lowest MAXI spectral hardness ratio of $\sim\!0.229$ corroborates the softest nature of the source. Moreover, we estimate the spin of the black hole by continuum-fitting method, which indicates that Cygnus X-1 is a maximally rotating `hole'. Further, Monte Carlo (MC) simulations are performed to estimate the uncertainty in spin parameter, which is constrained as a$_{\ast}>0.9981$ with $3\sigma$ confidence interval. Finally, we discuss the implications of our findings., Comment: 13 pages, 8 figures

We present a Suzaku X-ray study for the neutron star low mass X-ray binary Aquila X-1. The observations of the system carried out during three times in 2011 from October 18 to October 24, where the repeated X-ray transient Aquila X-1 was in the rising phase of the outburst. During this period the source was successfully detected in its different states, namely: the soft, transition (hard to soft) and hard state. Our results show that the broadband spectra detected in Oct.18 using XIS and HXD-PIN within the range from 0.8 to 50 KeV, have an absorbed flux of 3.7× 10 − 9 erg s−1 cm2 and luminosity of 1.3× 1037 erg s−1 (assuming a distance of 5.2 Kpc), although the luminosity is relatively high, the observed spectra confirm that the source in the hard state. The spectrum obtained on October 21 was detected over a 0.8 − 30 keV and displays a clear transitions from the hard state to the soft state at an absorbed flux of 1.07× 10 − 8 erg s−1 cm2 and luminosity of 3.45× 1037 ergs s−1. While the 24 October spectrum detected over a 0.8-30 keV, show a typical case of the soft state at an absorbed flux of 2.3× 10 − 8 erg s−1 cm2 and luminosity of 7.537× 1037 erg s−1. The continuum of the different spectra states successfully described by a two components model: a multi-color blackbody emission which originating from the accretion disk plus a Comptonized blackbody component originating from the boundary layer and/or neutron star surface. As the results, all model parameters showed change continuously from the ideal in the hard state to the ideal in the soft state. Prominently, the inner disk radius decreased from 41 km to 14 km, the effect of Comptonization changed from strong in the hard state to weak in the soft state by electron temperature which decreased from 21 keV to 3 keV. Inferred parameters sight that the coronal flow changed from optically thick in the hard state to optically thin flow in the hard state and the flow gradually shrink towards the soft state.

We present a new release of the RELTRANS model to fit the complex cross-spectrum of accreting black holes as a function of energy. The model accounts for continuum lags and reverberation lags self-consistently in order to consider the widest possible range of X-ray variability timescales. We introduce a more self-consistent treatment of the reverberation lags, accounting for how the time variations of the illuminating flux change the ionisation level of the accretion disc. This process varies the shape of the reflection spectrum in time causing an additional source of lags besides the light crossing delay. We also consider electron densities in the accretion disc up to $10^{20}$ cm$^{-3}$, which are found in most of the stellar mass black holes and in some AGN. These high densities increase the amplitude of the reverberation lags below $1$ keV since the reflection flux enhances in the same energy range. In addition, we investigate the properties of hard lags produced by variations in the power-law index of the continuum spectrum, which can be interpreted as due to roughly $3\%$ variability in the corona's optical depth and temperature. As a test case, we simultaneously fit the lag energy spectra in a wide range of Fourier frequency for the black hole candidate MAXI J1820+070 observed with NICER. The best fit shows how the reverberation lags contribute even at the longer timescales where the hard lags are important. This proves the importance of modelling these two lags together and self-consistently in order to constrain the parameters of the system., Accepted for publication in MNRAS

We present an analysis of all the archival high resolution spectra of the Narrow-line Seyfert 1 galaxy Mrk~335 obtained with Reflection Grating Spectrometer (RGS) on board \textit{XMM-Newton}. The spectra show rich emission and absorption features in low and intermediate flux intervals. We model the emission lines with the \textsc{pion\_xs} grid and try to find any possible correlation between the properties of the emitting gas and the source flux. Current data does not allow detailed trace of the response of the line emitting gas to the X-ray flux of Mrk~335, but the flux of the X-ray lines is significantly less variable than the X-ray continuum. We also find that the warm absorber's properties are not correlated with the flux variability. From the latest \textit{XMM-Newton} observation in 2019 December, we find that the photoionized emission and distant reflection components have not responded to the flux drop of Mrk~335 from 2018 July. However, the possible existence of partial covering absorber in the 2018--2019 low state of Mrk~335 makes it difficult to constrain the scale of the emitting gas using this lack of response., 11 pages, 6 figures

Simulations of the ionospheric response to solar flux changes driven by the 27 d solar rotation have been performed using the global 3-D Coupled Thermosphere Ionosphere Plasmasphere electrodynamics (CTIPe) physics-based numerical model. Using the F10.7 index as a proxy for solar extreme ultraviolet (EUV) variations in the model, the ionospheric delay at the solar rotation period is well reproduced and amounts to about 1 d, which is consistent with satellite and in situ measurements. From mechanistic CTIPe studies with reduced and increased eddy diffusion, we conclude that the eddy diffusion is an important factor that influences the delay of the ionospheric total electron content (TEC). We observed that the peak response time of the atomic oxygen to molecular nitrogen ratio to the solar EUV flux changes quickly during the increased eddy diffusion compared with weaker eddy diffusion. These results suggest that an increase in the eddy diffusion leads to faster transport processes and an increased loss rate, resulting in a decrease in the ionospheric time delay. Furthermore, we found that an increase in solar activity leads to an enhanced ionospheric delay. At low latitudes, the influence of solar activity is stronger because EUV radiation drives ionization processes that lead to compositional changes. Therefore, the combined effect of eddy diffusion and solar activity leads to a longer delay in the low-latitude and midlatitude region.

The use of continuum emission radio galaxies as cosmological tracers of the large-scale structure will soon move into a new phase. Upcoming surveys from the Australian Square Kilometre Array Pathfinder (ASKAP), MeerKAT, and the Square Kilometre Array project (SKA) will survey the entire available sky down to an ~100uJy flux limit, increasing the number of detected extra-galactic radio sources by several orders of magnitude. External data and machine learning algorithms will also enable some low resolution radial selection (photometric redshift binning) of the sample, increasing the cosmological utility of the sample observed. In this paper, we discuss the flux limit required to detect enough galaxies to decrease the shot noise term in the error to be 10% of the total. We show how future surveys of this type will be limited by available technology. The confusion generated by the intrinsic sizes of galaxies may have the consequence that surveys of this type eventually reach a hard flux limit of ~100nJy, as is predicted by the current modelling of AGN sizes by simulations such as the Tiered Radio Extragalactic Continuum Simulation (T-RECS). Finally, when considering the multi-tracer approach, where galaxies are split by type to measure some bias ratio, we find that there are not enough AGN present to achieve a reasonable level of shot noise for this kind of measurement., 10 pages, 8 figures, submitted to the journal

We present a systematic search for gamma-ray emission from supernovae (SNe) in the Fermi Large Area Telescope (LAT) Pass 8 data. The sample of targets consists of 55,880 candidates from the Open Supernova Catalog. We searched for gamma rays from SNe by means of a variable-size sliding-time-window analysis. Our results confirm the presence of transient gamma-ray emission from the sources of non-AGN classes, including transitional pulsars, solar flares, gamma-ray bursts, novae, and the Crab Nebula, which are projected near some of these SN's positions, and also strengthen support to the variable signal in the direction of SN iPTF14hls. The analysis is successful in finding both short (e.g. solar flares) and long (e.g. transitional pulsars) high flux states. Our search reveals two new gamma-ray transient signals occurred in 2019 in the directions of optical transients that are SN candidates, AT2019bvr and AT2018iwp, with their flux increases within 6 months after the dates of SN's discoveries. These signals are bright and their variability is at a higher statistical level than that of iPTF14hls. An exploration of archival multi-wavelength observations towards their positions is necessary to establish their association with SNe or other classes of sources. Our analysis, in addition, shows a bright transient gamma-ray signal at low Galactic latitudes in the direction of PSR J0205+6449. In addition, we report the results of an all-sky search for gamma-ray transient sources. This provided two additional candidates to gamma-ray transient sources., Published in MNRAS. Updated to address typos

We present results of phase- and time-resolved study of iron spectral features in the emission of the Be/X-ray transient pulsar V0332+53, during its type II outburst in 2004 using archival RXTE/PCA data. Coherent pulsations of both fluorescent iron line at {\approx}6.4 keV and neutral iron K-edge at {\approx}7.1 keV have been detected throughout the entire outburst. The pulsating iron K-edge is reported for the first time for this object. Near the peak of the outburst, the 3-12 keV pulse profile shows two deep, F_max / F_min ~ 2, and narrow dips of nearly identical shape, separated by exactly {\Delta}{\phi}=0.5 in phase. The dip spectra are nearly identical to each other and very similar in shape to the spectra outside the dips. The iron K-edge peaks at the phase intervals corresponding to the dips, although its optical depth {\tau}_K ~ 0.05 is by far insufficient to explain the dips. The iron line shows pulsations with a complex pulse profile without any obvious correlation with the total flux or optical depth of the K-edge. Accounting for the component associated with reprocessing of the pulsar emission by the surface of the donor star and circumstellar material, we find a very high pulsation amplitude of the iron line flux, F_max / F_min ~ 10. We demonstrate that these properties of V0332+53 can not be explained by contemporary emission models for accreting X-ray pulsars and speculate about the origin of the observed iron spectral features., Comment: 16 pages, 6 figures. Monthly Notices of the Royal Astronomical Society Main Journal Accepted

This paper examines long-term (more than four solar cycles) temporal and spatial fluctuations in the solar rotation by investigating radio-emission escapes from various layers of the solar atmosphere during the years 1967–2010. The flux modulation approach can also be used to investigate variations in solar rotation, which is a contentious topic in solar physics. This study makes use of a time-series of radio flux data at various frequencies (245–15 400 MHz) obtained at Sagamore Hill Solar Radio Observatory in Massachusetts, USA, and at other observatories from 1967 to 2010. The periodicity present in the temporal variation of the time-series is estimated through a Lomb–Scargle periodogram. The rotation period estimated for five radio emissions (606, 1415, and 2695 MHz from the corona, and 4995 and 8800 MHz from the transition region) through a statistical approach shows continuous temporal and spatial variations throughout the years. The smoothed rotation period shows the presence of periodic ∼22-yr and ∼11-yr components. The 22-yr component could be linked to the reversal of the solar magnetic field (Hale) cycle, while the 11-yr component is most likely related to the sunspot (Schwabe) cycle. In addition to these two components, random components are also prominently present in the analysed data. The cross-correlation between the sunspot number and the rotation period obtained shows a strong correlation with the 11-yr Schwabe and 22-yr Hale cycle. The corona rotates faster or slower than the transition region in different epochs. The alternation of the faster rotation speed between the corona and transition region also follows the 22-yr cycle.

(Abridged) Several phenomena in astrophysics generate light curves with time delays. Among these are reverberation mapping, and lensed quasars. In some systems, the measurement of the time-delay is complicated by the fact that the delayed components are unresolved and that the light curves are generated from a red-noise process. We derive the likelihood function of the observations given a model of either a combination of time-delayed light curves or a single light curve. This likelihood function is different from the auto-correlation function. We demonstrate that given a single-band light curve that is a combination of two (or more) time-shifted copies of an original light curve, generated from a red-noise probability distribution, we can test if the total-flux light curve is a composition of time-delayed copies or, alternatively, is consistent with being the original light curve. Furthermorew, in some realistic cases, it is possible to measure the time delays and flux ratios between these unresolved components even when the flux ratio is about 1/10. This method is useful for identifying lensed quasars and simultaneously measuring their time delays, and for estimating the reverberation time scales of active galactic nuclei. In a companion paper, we derive a method that uses the center-of-light position (e.g., of a lensed quasar) along with the combined flux. This allow us to identify lensed quasars and supernovae and measure their time delays, with higher fidelity compared to the flux-only method. The astrometry + flux method, however, is not suitable for quasar reverberation mapping. We also comment on the commonly used method of fitting a power-law model to a power spectrum, and present the proper likelihood function for such a fit. We test the new method on simulations and provide Python and MATLAB implementations., Submitted to MNRAS, 16 pp, 17 figures

We present a kinematic study of ionized extraplanar gas in two low-inclination late-type galaxies (NGC 3982 and NGC 4152) using integral field spectroscopy data from the DiskMass H α sample. We first isolate the extraplanar gas emission by masking the H α flux from the regularly rotating disc. The extraplanar gas emission is then modelled in the 3D position–velocity domain using a parametric model described by three structural and four kinematic parameters. Best-fitting values for the model are determined via a Bayesian MCMC approach. The reliability and accuracy of our modelling method are carefully determined via tests using mock data. We detect ionized extraplanar gas in both galaxies, with scale heights $0.83^{+0.27}_{-0.40}\, \mathrm{kpc}$ (NGC 3982) and $1.87^{+0.43}_{-0.56}\, \mathrm{kpc}$ (NGC 4152) and flux fraction between the extraplanar gas and the regularly rotating gas within the disc of 27 and 15 per cent, respectively, consistent with previous determinations in other systems. We find lagging rotation of the ionized extraplanar gas in both galaxies, with vertical rotational gradients $-22.24^{+6.60}_{-13.13} \, \mathrm{km\, s^{-1}\, kpc^{-1}}$ and $-11.18^{+3.49}_{-4.06}\, \mathrm{km\, s^{-1}\, kpc^{-1}}$, respectively, and weak evidence for vertical and radial inflow in both galaxies. The above results are similar to the kinematics of the neutral extraplanar gas found in several galaxies, though this is the first time that 3D kinematic modelling of ionized extraplanar gas has been carried out. Our results are broadly consistent with a galactic fountain origin combined with gas accretion. However, a dynamical model is required to better understand the formation of ionized extraplanar gas.

The relationship between mass and radius (M-R relation) is the key for inferring the planetary compositions and thus valuable for the studies of formation and migration models. However, the M-R relation alone is not enough for planetary characterization due to the dependence of it on other confounding variables. This paper provides a non-trivial extension of the M-R relation by including the incident flux as an additional variable. By using Bayesian hierarchical modeling (BHM) that leverages the flexibility of finite mixture models, a probabilistic mass-radius-flux relationship (M-R-F relation) is obtained based on a sample of 319 exoplanets. We find that the flux has nonnegligible impact on the M-R relation, while such impact is strongest for hot-Jupiters. On the population level, the planets with higher level of flux tend to be denser, and high flux could trigger significant mass loss for plants with radii larger than $13R_{\oplus}$. As a result, failing to account for the flux in mass prediction would cause systematic over or under-estimation. With the recent advent of computing power, although a lot of complex statistical models can be fitted using Monte Carlo methods, it has largely remain illusive how to validate these complex models when the data are observed with large measurement errors. We present two novel methods to examine model assumptions, which can be used not only for the models we present in this paper but can also be adapted for other statistical models., Accepted for publication in MNRAS

Results of 111-MHz monitoring observations carried out on the Big Scanning Antenna of the Pushchino Radio Astronomy Observatory during September 1-28, 2015 are presented. Fifty-four pulsating sources were detected at declinations $-9^o < δ< +42^o$. Forty-seven of these are known pulsars, five are new sources, and two are previously discovered transients. Estimates of the peak flux densities and dispersion measures are presented or all these sources., published in Astronomy Report, translated by Yandex translator with correction of scientific lexis, 8 pages, 2 figures, 3 tables

Context. The Rosetta spacecraft followed comet 67P/Churyumov-Gerasimenko (67P) for more than two years at a slow walking pace (~1 m s−1) within 1500 km from the nucleus. During one of the radial movements of the spacecraft in the early phase of the mission, the radial distribution of the plasma density could be estimated, and the ionospheric density was found to be inversely proportional to the cometocentric distance r from the nucleus (a 1/r distribution). Aims. This study aims to further characterise the radial distribution of plasma around 67P throughout the mission and to expand on the initial results. We also aim to investigate how a 1/r distribution would be observed during aflyby with a fast (~10’s km s−1) spacecraft, such as the upcoming Comet Interceptor mission, when there is also an asymmetry introduced to the outgassing over the comet surface. Methods. To determine the radial distribution of the plasma, we used data from the Langmuir probe and Mutual Impedance instruments from the Rosetta Plasma Consortium during six intervals throughout the mission, for which the motion of Rosetta was approximately radial with respect to the comet. We then simulated what distribution a fast flyby mission would actually observe during its passage through a coma when there is a 1/r plasma density distribution as well as a sinusoidal variation with a phase angle (and then a sawtooth variation) multiplied to the outgassing rate. Results. The plasma density around comet 67P is found to roughly follow a 1/r dependence, although significant deviations occur in some intervals. If we normalise all data to a common outgassing rate (or heliocentric distance) and combine the intervals to a radial range of 10–1500 km, we find a 1/r1.19 average distribution. The simulated observed density from a fast spacecraft flying through a coma with a 1/r distribution and an asymmetric outgassing can, in fact, appear anywhere in the range from a 1/r distribution to a 1/r2 distribution, or even slightly outside of this range. Conclusions. The plasma density is distributed in such a way that it approximately decreases in a manner that is inversely proportional to the cometocentric distance. This is to be expected from the photoionisation of a collision-less, expanding neutral gas at a constant ionisation rate and expansion speed. The deviation from a pure 1/r distribution is in many cases caused by asymmetric outgassing over the surface, additional ionisation sources being present, electric fields accelerating plasma, and changing upstream solar wind conditions. A fast flyby mission can observe a radial distribution that deviates significantly from a 1/r trend if the outgassing is not symmetric over the surface. The altitude profile that will be observed depends very much on the level of outgassing asymmetry, the flyby velocity, the comet rotation rate, and the rotation phase. It is therefore essential to include data from both the inbound and outbound legs, as well as to compare plasma density to neutral density to get a more complete understanding of the radial distribution of the plasma.

After updating the status of the measurements of the cosmic neutrino flux by the IceCube experiment, we summarize the observations of the first identified source of cosmic rays and speculate on the connection between the two observations.

Astronomy and astrophysics 662, A65 (2022). doi:10.1051/0004-6361/202243096, Observations with imaging atmospheric Cherenkov telescopes (IACTs) have enhanced our knowledge of nearby supernova (SN) remnants with ages younger than 500 yr by establishing Cassiopeia A and the remnant of Tycho’s SN as very-high-energy (VHE) γ-ray sources. The remnant of Kepler’s SN, which is the product of the most recent naked-eye SN in our Galaxy, is comparable in age to the other two, but is significantly more distant. If the γ-ray luminosities of the remnants of Tycho’s and Kepler’s SNe are similar, then the latter is expected to be one of the faintest γ-ray sources within reach of the current generation IACT arrays. Here we report evidence at a statistical level of 4.6σ for a VHE signal from the remnant of Kepler’s SN based on deep observations by the High Energy Stereoscopic System (H.E.S.S.) with an exposure of 152 h. The measured integral flux above an energy of 226 GeV is ∼0.3% of the flux of the Crab Nebula. The spectral energy distribution (SED) reveals a γ-ray emitting component connecting the VHE emission observed with H.E.S.S. to the emission observed at GeV energies with Fermi-LAT. The overall SED is similar to that of the remnant of Tycho’s SN, possibly indicating the same nonthermal emission processes acting in both these young remnants of thermonuclear SNe., Published by EDP Sciences, Les Ulis

Paleogeography is the study of the changing surface of Earth through time. Driven by plate tectonics, the configuration of the continents and ocean basins has been in constant flux. Plate tectonics pushes the land surface upward or pulls it apart, causing its collapse. All the while, the unrelenting forces of climate and weather slowly reduce mountains to sand and mud and redistribute these sediments to the sea. This article reviews the changing paleogeography of the past 750 million years. It describes the broad patterns of Phanerozoic paleogeography as well as many of the specific paleogeographic events that have shaped the modern continents and ocean basins. The focus is on the changing latitudinal distribution of the continents, fluctuations in sea level, the opening and closing of oceanic seaways, mountain building, and how these paleogeographic changes have affected global climate, ocean circulation, and the evolution of life. This review presents an atlas of 114 paleogeographic maps that illustrate how Earth's surface has evolved during the past 750 million years. During that time interval, Earth has witnessed the formation and breakup of two supercontinents: Pannotia and Pangea. The continents have been transformed from low-lying flooded platforms to high-standing land areas crisscrossed by the scars of past continental collisions. Oceans have opened and closed, and then opened again in a seemingly never-ending cycle. ▪ The changing configuration of the continents and ocean basins during the past 750 million years is illustrated in 114 paleogeographic maps. ▪ These maps describe how the surface of Earth has been continually modified by mountain building and erosion. ▪ The changing paleogeography has affected global climate, ocean circulation, and the evolution of life. ▪ The data and methods used to produce the maps are described in detail.

The properties of axisymmetric magnetohydrodynamic wave modes are studied for the model of coronal magnetic tube, consisting of a central cord with homogeneous axial magnetic field, surrounded by an annulus with twisted magnetic field. This model mimics the coronal loops with radially localized magnetic twist. The derived dispersion relation is solved numerically. A number of limiting cases are examined analytically. The two families of axisymmetric modes are found to exist in the model. The first one includes an infinite number of fast-sausage modes (FSMs), modified by the twist and the second one is a set of modes with frequencies, lying in a narrow band, closed to Alfvén frequency of a twisted layer ωA0. The fundamental FSM (the mode of the lowest radial order) exists as a trapped mode for the entire range of axial wavenumbers. Its phase speed is always below the Alfvén speed of a magnetically twisted layer. This mode has a weak dispersion in the range of long and intermediate wavelengths. The higher radial order FSMs were found to be less sensitive to the presence of the magnetic twist. Both the fundamental FSM in the case of a weak magnetic twist and a set of symmetric modes for arbitrary twist have very similar frequencies, which are nearly equal to ωA0. This property implies that aforementioned wave modes behave like Alfvén torsional modes in the twisted annulus and like the radial modes in the untwisted core and environment, moreover they are strongly coupled. The main results of the study are discussed in the framework of their applications to coronal seismology.

We present $^{56}$Ni mass estimates for 110 normal Type II supernovae (SNe II), computed here from their luminosity in the radioactive tail. This sample consists of SNe from the literature, with at least three photometric measurements in a single optical band within 95-320 d since explosion. To convert apparent magnitudes to bolometric ones, we compute bolometric corrections (BCs) using 15 SNe in our sample having optical and near-IR photometry, along with three sets of SN II atmosphere models to account for the unobserved flux. We find that the $I$- and $i$-band are best suited to estimate luminosities through the BC technique. The $^{56}$Ni mass distribution of our SN sample has a minimum and maximum of 0.005 and 0.177 M$_{\odot}$, respectively, and a selection-bias-corrected average of $0.037\pm0.005$ M$_{\odot}$. Using the latter value together with iron isotope ratios of two sets of core-collapse (CC) nucleosynthesis models, we calculate a mean iron yield of $0.040\pm0.005$ M$_{\odot}$ for normal SNe II. Combining this result with recent mean $^{56}$Ni mass measurements for other CC SN subtypes, we estimate a mean iron yield $$36 per cent. We also find that the empirical relation between $^{56}$Ni mass and steepness parameter ($S$) is poorly suited to measure the $^{56}$Ni mass of normal SNe II. Instead, we present a correlation between $^{56}$Ni mass, $S$, and absolute magnitude at 50 d since explosion. The latter allows to measure $^{56}$Ni masses of normal SNe II with a precision around 30 per cent., Comment: 33 pages, 20 figures, 6 figures in appendix, accepted for publication to MNRAS

One of the mainstreams in modern X-ray astronomy is research into extragalactic X-ray sources on the basis of the data acquired at the X-ray Multi-Mirror Newton (XMM-Newton) space observatory. According to observations, X-rays coming from galaxies are mainly radiated from their central regions, i.e., active galactic nuclei and groups of X-ray sources in the galactic disks. In this paper, we consider the cross-correlation between the 4XMM-DR9 catalog and the Hyper-Linked Extragalactic Databases and Archives (HyperLeda) of galaxies. The 4XMM-DR9 catalog is a large, up-to-date catalog of observations, which contains 550 124 unique sources and covers 2.85% of the sky, while the HyperLeda database comprises 1.5 million galaxies. Our analysis resulted in a sample of more than 5000 X-ray galaxies, most of which are active galactic nuclei of low luminosity. From this sample, we selected galaxies whose the X-ray flux exceeds F = 10–20 J/cm2s. The sources of this kind are of particular interest since it is easier to construct an informative spectrum for them. The identified and classified catalog of 1172 manually verified galaxies—the X-ray galaxy catalog named Xgal—was created. In the Xgal catalog, most galaxies have an active X-ray nucleus; Seyfert galaxies predominate among them at short distances, while quasars are prevalent at large distances. We revealed 169 galaxies that exhibit extended nuclei with a visible surface brightness distribution and 173 galaxies with more than one X-ray source. Based on the Xgal catalog, we created a catalog of elongated X-ray galaxies (the optical angular sizes of which are a > 60″) that have X-ray sources outside the nucleus. Both catalogs are freely accessible. The Xgal catalog may serve to construct the spectra of objects of a certain class in different ranges, to develop or improve the theory of their emission, and to survey bright and extended quasars. Moreover, the entire cross-sample may be used to study active galactic nuclei with low luminosity and a large-scale structure of the universe in the X-ray range.

We present results from \textit{Chandra} X-ray observations and 325 MHz Giant Metrewave Radio Telescope (GMRT) observations of the massive and X-ray luminous cluster of galaxies Abell S1063. We report the detection of large-scale \lq\lq excess brightness\rq\rq\ in the residual \textit{Chandra} X-ray surface brightness map, which extends at least 2.7 Mpc towards the north-east from the center of the cluster. We also present a high fidelity X-ray flux and temperature map using \textit{Chandra} archival data of 122 ksec, which shows the disturbed morphology in the cluster. The residual flux map shows the first observational confirmation of the merging axis proposed by earlier simulation by \citet{Gomez2012AJ....144...79G}. The average temperature within $R_{500}$ is $11.7 \pm 0.56$ keV, which makes AS1063 one of the hottest clusters in the nearby Universe. The integrated radio flux density at 325 MHz is found to be $62.0\pm6.3$ mJy. The integrated spectrum of the radio halo follows a power-law with a spectral index $\alpha=-1.43\pm 0.13$. The radio halo is found to be significantly under-luminous, which favored for both the hadronic as well as the turbulent re-acceleration mechanism for its origin., Comment: Accepted for publication in MNRAS (12 pages, 8 figures, 2 tables)

Emission-line variability studies have the potential to provide key information about the structures of the near-nuclear outflow regions of AGN. Here we present a VLT/Xshooter spectrum of the nucleus of the ULIRG F01004-2237 that was taken in August 2018, about 8 yr after a major flare in its integrated optical emission. Compared with our WHT/ISIS spectrum from September 2015, the broad, red wings of the emission lines most closely associated with the flaring event, including HeII(4686), NIII(4100,4640) and HeI(5876), have substantially declined in flux. In contrast, the broad, blue wings that dominate the [OIII], [NeIII], [NeV] and [OI] forbidden lines have increased in flux by a factor 1.4 - 4.4 (depending on the line). Moreover, the [FeVII](6087) line is detected in the new spectrum for the first time. We interpret these results in terms of a light echo from the outflowing warm gas: the direct emission from the flaring event is continuing to fade, but due to light travel time effects we are only now observing the impact of the flare on the emission from the extended outflow region. Unless the outflow is confined to a small range of angles close to our line of sight, these observations imply that the outflow must be relatively compact (r, 12 pages, 8 figures, accepted for publication in MNRAS

We use the Gemini NIFS instrument to map the H$_2 2.1218\mu$m and Br$\gamma$ flux distributions in the inner 0.04-2 kpc of a sample of 36 nearby active galaxies ($0.001\lesssim z\lesssim0.056$) at spatial resolutions from 4 to 250 pc. We find extended emission in 34 galaxies. In $\sim$55% of them, the emission in both lines is most extended along the galaxy major axis, while in the other 45% the extent follows a distinct orientation. The emission of H$_2$ is less concentrated than that of Br$\gamma$, presenting a radius that contains half of the flux 60% greater, on average. The H$_2$ emission is driven by thermal processes - X-ray heating and shocks - at most locations for all galaxies, where $0.46$ (seen in 40% of the galaxies), shocks are the main H$_2$ excitation mechanism, while in regions with H$_2$/Br$\gamma, Comment: 21 pages, accepted for publication in MNRAS

We report here results on the analysis of correlated flux variations between the optical and GeV $\gamma$-ray bands in three bright BL Lac objects, namely AO\, 0235+164, OJ 287 and PKS 2155$-$304. This was based on the analysis of about 10 years of data from the {\it Fermi} Gamma-ray Space Telescope covering the period between 08 August 2008 to 08 August 2018 along with optical data covering the same period. For all the sources, during the flares analysed in this work, the optical and $\gamma$-ray flux variations are found to be closely correlated. From broad band spectral energy distribution modelling of different epochs in these sources using the one zone leptonic emission model, we found that the optical-UV emission is dominated by synchrotron emission from the jet. The $\gamma$-ray emission in the low synchrotron peaked sources AO\, 0235+164 and OJ 287 are found to be well fit with external Compton (EC) component, while, the $\gamma$-ray emission in the high synchrotron peaked source PKS 2155$-$304 is well fit with synchrotron self Compton component. Further we note that the $\gamma$-ray emission during the high flux state of AO 0235+164 (epochs A and B) requires seed photons from both the dusty torus and broad line region, while the $\gamma$-ray emission in OJ 287 and during epochs C and D of AO\,0235+164 can be modelled by EC scattering of infra-red photons from the torus., Comment: 17 pages, 26 figures, Accepted for publication in MNRAS (In the previous edition, a minor correction was made.)

PSR J1420-6048 is a young gamma-ray pulsar with recurrent glitches. Utilizing long-term monitoring data obtained from the Fermi Gamma-ray Space Telescope, we found that PSR J1420-6048 has shown gamma-ray flux variation and we also detected four glitches between 2008 and 2019. Two of the glitches are previously unknown, and their gamma-ray spectrum also shows variability between each glitch. Since the results might be contaminated by background sources, we discuss whether the observed changes in flux and spectra were caused by artificial misallocations of photons from a nearby pulsar wind nebula (HESS J1420-607) and a pulsar (PSR J1418-6058), or a change of the emission geometry from the target pulsar itself. We examine the correlation of the flux changes and the alternating pulse structure to investigate whether the emission geometry in the outer magnetosphere was changing. By assuming the observational features were not totally resulted from the background environment, we compare our results with similar phenomena observed in other gamma-ray pulsars and propose that a strong crust crack can cause timing anomaly of a neutron star, which can affect the particle accelerations or pair creation regions resulting in the changes of emission behaviors., Comment: 11 pages, 5 figures and 1 table; accepted by MNRAS

Model of a partial current-carrying torus loop anchored to the photosphere is analyzed. Conditions of the catastrophic loss of equilibrium are considered and corresponding value of the critical decay index of external magnetic field is found. Taking into account line-tying conditions leads to non-monotonous dependence of the critical decay index on the height of the apex and length of the flux rope (its endpoints separation). For relatively short flux ropes, the critical decay index is significantly lower than unity, which is in contrast to widespread models with the typical critical decay index above unity. The steep decrease of the critical index with height at low heights is due to the sharp increase of the curvature of the flux-rope axis that transforms from a nearly straight line to a crescent., Comment: 11 pages, 4 figures. Accepted for publication in MNRAS