Long-term multi-wavelength variability and extreme spectral properties of the TeV-detected blazar 1ES 0033+595

Context. The multi-wavelength (MWL) properties of some TeV-detected blazars are still poorly understood. By studying the MWL properties of the BL Lacertae source 1ES 0033+595, we make further progress in our understanding of the instable physical processes responsible for particle acceleration to ultra-relativitic energies and the observed emission properties. Aims. This paper presents the results of a detailed X-ray timing and spectral analysis of the source 1ES 0033+595, allowing us to draw conclusions about the physical mechanisms responsible for particle acceleration and the generation of X-ray emission. We also examined the long-term MWL behaviour of the source and interband cross-correlations. Methods. Our study focused on the observations performed with the X-Ray Telescope on board the Neil Gehrels Swift Observatory (Swift-XRT) in the period 2005–2022. The MWL data sets were obtained with several space and ground-based instruments. We analysed the MWL flaring behaviour during the entire period, as well as in particular subperiods selected on the basis of the X-ray flaring activity of the target; constructed histograms characterising the distributions of the different spectral parameters and MWL fluxes; and checked cross-correlations between these parameters and/or MWL fluxes. The obtained results are compared with the predictions of the various theoretical studies and simulations presented by the different authors to date. Results. The source showed extreme spectral properties with dominance of high spectral curvature, frequent occurrence of very and extremely hard photon indices, and the presence of a synchrotron spectral energy distribution (SED) peak in the hard X-ray range (sometimes at energies beyond 10 keV). These properties reveal the importance of first-order Fermi acceleration with very low initial particle energy distribution, along with the co-existence of stochastic acceleration and hadronic processes. The source was characterised by very uneven and erratic flaring activity in diverse epochs: the period of strong flares (2013–2016) was preceded by a moderate variability and followed by a gradual long-term decrease in MWL flaring activity. We detected a number of instances of intraday 0.3–10 keV variability, which were sometimes observed within a few hundred seconds and explained by the interaction between the relativistic shock front and jet inhomogeneities with strong magnetic fields. The X-ray and γ-ray fluxes showed a lognormal distribution, which hints at the imprint of accretion disc instabilities on the blazar jet