Context. It has become evident that one-zone synchrotron self-Compton models are not always adequate for very-high-energy (VHE) gamma-ray emitting blazars. While two-component models are performing better, they are difficult to constrain due to the large number of free parameters. Aims. In this work, we make a first attempt to take into account the observational constraints from Very Long Baseline Interferometry (VLBI) data, long-term light curves (radio, optical, and X-rays) and optical polarisation to limit the parameter space for a two-component model and test if it can still reproduce the observed spectral energy distribution (SED) of the blazars. Methods. We selected five TeV BL Lac objects based on the availability of VHE gamma-ray and optical polarisation data. We collected constraints for the jet parameters from VLBI observations. We evaluated the contributions of the two components to the optical flux by means of decomposition of long-term radio and optical light curves as well as modeling of the optical polarisation variability of the objects. We selected eight epochs for these five objects, based on the variability observed at VHE gamma rays, for which we constructed the SEDs that we then modeled with a two-component model. Results. We found parameter sets which can reproduce the broadband SED of the sources in the framework of two-component models considering all available observational constraints from VLBI observations. Moreover, the constraints obtained from the long-term behavior of the sources in the lower energy bands could be used to determine the region where the emission in each band originates. Finally, we attempted to use optical polarisation data to shed new light on the behavior of the two components in the optical band. Our observationally constrained two zone model allows explanation of the entire SED from radio to VHE with two co-located emission regions.

We report on detailed spectral analysis of emission from X-ray pulsar 4U 1901+03 using the data obtained by the Insight-HXMT and NuSTAR observatories during the 2019 outburst of the source. Thanks to the extensive coverage of the outburst by Insight-HXMT, we were able to investigate spectral evolution of the source as a function of flux, and compare these results to the previous reports with focus on the properties of a putative absorption feature at around 10 keV. In particular, we demonstrate that the broadband X-ray continuum of 4U 1901+03 can be well described with a two-component continuum model without absorption line at 10 keV, which casts doubt on its interpretation as a cyclotron line. The high quality of the data has also allowed us to perform both phase-averaged and phase-resolved spectral analysis as function of luminosity. Finally, we performed a detailed investigation of another absorption feature in the spectrum of the source around 30 keV recently reported in the NuSTAR data. In line with the previous findings, we show that this feature appears to be significantly detected, albeit only in the pulse-phase resolved spectra at relatively high luminosities. The non-detection of the line in the averaged spectrum is possibly associated with the observed strong variations of its parameters with the pulse phase.

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.

We report on the analysis of the spin evolution of a slow X-ray pulsar GX 301$-$2 along the orbit using long-term monitoring by Fermi/GBM. Based on the observationally confirmed accretion scenario and an analytical model for the accretion of angular momentum we demonstrate that in this system, the neutron star spins retrogradely, that is, in a direction opposite to the orbital motion. This first-of-a-kind discovery of such a system proves the principal possibility of retrograde rotation in accreting systems with suitable accretion torque, and might have profound consequences for our understanding of the spin evolution of X-ray pulsars, estimates of their initial spin periods, and the ultimate result of their evolution.

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.

We present the X-ray analysis of dwarf nova EY Cyg using the 45 ks \textit{XMM-Newton} observatory archival data obtained in quiescence. We find orbital modulations in X-rays. We simultaneously fitted EPIC pn, MOS1 and MOS2 data using a model for interstellar medium absorption (\textit{tbabs}) and a multi-temperature plasma emission model with a power-law distribution of temperatures (CEVMKL) as expected from low accretion rate quiescent dwarf novae. The \textit{XMM-Newton} EPIC spectra of the source yields a maximum temperature $kT_{\rm max}$ $\sim$ 14.9$^{+3.3}_{-2.2}$ keV with an unabsorbed X-ray flux and luminosity of (1.8--2.0) $\times$ 10$^{-12}$ ergs$^{-1}$ cm$^{-1}$ and (8.7--9.7) $\times$ 10$^{31}$ ergs$^{-1}$, respectively, in the energy range 0.1 to 50 keV. There is 3--4 sigma excess at energies below 0.5 keV, we model the excess using \textit{MEKAL}, POWERLAW and BBODY models and favor the model \textit{MEKAL} which is physical. According to previous studies, the secondary in this system is thought to be a K-type star which may radiate in the soft X-ray region. The fit with an additive \textit{MEKAL} model gives a temperature of $kT$ $\sim$ 0.1 keV with an unabsorbed X-ray flux and luminosity of (2.7--8.8) $\times$ 10$^{-14}$ ergs$^{-1}$ cm$^{-1}$ and (1.3--4.2) $\times$ 10$^{30}$ ergs$^{-1}$, respectively, for the companion star. Based on the results from the timing and spectral analysis, we highly suggest that the secondary of EY Cyg is a K-type star.

We report on \textit{NuSTAR} observations of the well-known wind-accreting X-ray pulsar \source\ during a strong spin-up episode that took place in January-March 2019. A high luminosity of the source in a most recent observation allowed us to detect a positive correlation of the cyclotron line energy with luminosity. Beyond that, only minor differences in spectral and temporal properties of the source during the spin-up, presumably associated with the formation of a transient accretion disk, and the normal wind-fed state could be detected. We finally discuss conditions for the formation of the disk and possible reasons for lack of any appreciable variations in most of the observed source properties induced by the change of the accretion mechanism, and conclude that the bulk of the observed X-ray emission is still likely powered by direct accretion from the wind.

We present results of investigation of the poorly studied X-ray pulsar Swift J1816.7--1613 during its transition from the type I outburst to the quiescent state. Our studies are based on the data obtained from X-ray observatories \textit{Swift}, \textit{NuSTAR} and \textit{Chandra} alongside with the latest IR data from UKIDSS/GPS and \textit{Spitzer}/GLIMPSE surveys. The aim of the work is to determine parameters of the system: the strength of the neutron star magnetic field and the distance to the source, which are required for the interpretation of the source behaviour in the framework of physically motivated models. No cyclotron absorption line was detected in the broad-band energy spectrum. However, the timing analysis hints at the typical for the X-ray pulsars magnetic field from a few $\times 10^{11}$ to a few $\times 10^{12}$ G. We also estimated type of the IR-companion as a B0-2e star located at distance of 7--13~kpc.