We present the results of the \textit{Fermi}-LAT 10-years-long light curves (LCs) modeling of selected blazars: six flat spectrum radio quasars (FSRQs) and five BL Lacertae (BL Lacs), examined in 7-, 10-, and 14-day binning. The LCs and power spectral densities (PSDs) were investigated with various methods: Fourier transform, Lomb-Scargle periodogram (LSP), wavelet scalogram, autoregressive moving average (ARMA) process, continuous-time ARMA (CARMA), Hurst exponent ($H$), and the $\mathcal{A}-\mathcal{T}$ plane. First, with extensive simulations we showed that parametric modeling returns unreliable parameters, with a high dispersion for different realizations of the same stochastic model. Hence any such analysis should be supported with Monte Carlo simulations. For our blazar sample, we find that the power law indices $\beta$ calculated from the Fourier and LSP modeling mostly fall in the range $1\lesssim\beta\lesssim 2$. Using the wavelet scalograms, we confirm a quasi-periodic oscillation (QPO) in PKS~2155$-$304 at a $3\sigma$ significance level, but do not detect any QPOs in other objects. The ARMA fits reached higher orders for 7-day binned LCs and lower orders for 10- and 14-days binned LCs for the majority of blazars, suggesting there might exist a characteristic timescale for the perturbations in the jet and/or accretion disk to die out. ARMA and CARMA modeling revealed breaks in their PSDs at timescales of few hundred days. The estimation of $H$ was performed with several methods. We find that most blazars exhibit $H > 0.5$, indicating long-term memory. Finally, the FSRQ and BL Lac subclasses are clearly separated in the $\mathcal{A}-\mathcal{T}$ plane.

In order to explain the origin of skewness in the gamma-ray burst (GRB) duration distribution, a statistical model, based on the cosmological time dilation and taking into account the GRB formation rate and instrumental effects, is examined. It is concluded that the considered effects are able to account for only a small fraction of the observed skewness. Therefore, its origin needs to be searched for in the processes governing the progenitors of GRBs.

It was recently shown that the $T_{90}-H_{32}$ distributions of gamma-ray bursts from CGRO/BATSE and Fermi/GBM are well described by a mixture of only two skewed components, making the presumed third, intermediate class unnecesary. The Swift/BAT, Konus-Wind, RHESSI and Suzaku/WAM data sets are found to be consistent with a two-class description as well.

We present an optical variability study of 44 newly identified blazar candidates behind the Magellanic Clouds, including 27 flat spectrum radio quasars (FSRQs) and 17 BL Lacertae objects (BL Lacs). All objects in the sample possess high photometric accuracy and irregularly sampled optical light curves (LCs) in I filter from the long-term monitoring conducted by the Optical Gravitational Lensing Experiment. We investigated the variability properties to look for blazar-like characteristics and to analyze the long-term behaviour. We analyzed the LCs with the Lomb-Scargle periodogram to construct power spectral densities (PSDs), found breaks for several objects, and linked them with accretion disk properties. In this way we constrained the black hole (BH) masses of 18 FSRQs to lie within the range $8.18\leq\log (M_{\rm BH}/M_\odot)\leq 10.84$, assuming a wide range of possible BH spins. By estimating the bolometric luminosities, we applied the fundamental plane of active galactic nuclei variability as an independent estimate, resulting in $8.4\leq\log (M_{\rm BH}/M_\odot)\leq 9.6$, with a mean error of 0.3. Many of the objects have very steep PSDs, with high frequency spectral index in the range $3-7$. An alternative attempt to classify the LCs was made using the Hurst exponent, $H$, and the $\mathcal{A}-\mathcal{T}$ plane. Two FSRQs and four BL Lacs yielded $H>0.5$, indicating presence of long-term memory in the underlying process governing the variability. Additionally, two FSRQs with exceptional PSDs, stand out also in the $\mathcal{A}-\mathcal{T}$ plane.

We present results of a variability study in the optical band of 44 newly identified blazar candidates behind the Magellanic Clouds. Our sample contains 27 flat spectrum radio quasars (FSRQs) and 17 BL Lacertae objects (BL Lacs). However, only nine of them are considered as secure blazar candidates, while the classification of the remaining 35 objects is still uncertain. All studied blazar candidates possess infrequently sampled optical light curves (LCs) in I filter provided by the Optical Gravitational Lensing Experiment group. The LCs were analysed with the Lomb-Scargle periodogram, the Hurst exponent $H$, and the $\mathcal{A}-\mathcal{T}$ plane, to look for blazar-like characteristic features and to study the long-term behaviour of the optical fluxes. The power law (PL) indices of the Lomb-Scargle power spectral density (PSD) of the FSRQ blazar candidates mostly lie in the range (1,2). In case of the BL Lacs they are located in the range (1,1.8). The PL PSD is indicative of a self-affine stochastic process characterised by $H$, underlying the observed variability. We find that the majority of analysed objects have $H\leq 0.5$, indicating short-term memory, whereas four BL Lacs and two FSRQs have $H>0.5$, implying long-term memory. 41 blazar candidates are located in the $\mathcal{A}-\mathcal{T}$ plane in the region available to PL plus Poisson noise processes. Interestingly, one FSRQ is located marginally below this region, while two FSRQs lie above the line $\mathcal{T}=2/3$, i.e. they are even more noisy than white noise. The BL Lac candidates are characterised by higher $\mathcal{A}$ values than FSRQs, i.e. $0.71\pm 0.06$ and $0.29\pm 0.05$, respectively.

The number of classes of gamma-ray bursts (GRBs), besides the well-established short and long ones, remains a debatable issue. It was already shown, however, that when invoking skewed distributions, the $\log T_{90}$ and $\log T_{90} - \log H_{32}$ spaces are adequately modeled with mixtures of only two such components, implying two GRB types. Herein, a comprehensive multivariate analysis of several multi-dimensional parameter spaces is conducted for the BATSE sample of GRBs, with the usage of skewed distributions. It is found that the number of extracted components varies between the examined parameter sets, and ranges from 2 to 4, with higher-dimensional spaces allowing for more classes. A Monte Carlo testing implies that these additional components are likely to be artifacts owing to the finiteness of the data and be a result of examining a particular realization of the data as a random sample, resulting in spurious identifications.

Gamma-ray bursts (GRBs) were confirmed to be of extragalactic origin due to their isotropic angular distribution, combined with the fact that they exhibited an intensity distribution that deviated strongly from the $-3/2$ power law. This finding was later confirmed with the first redshift, equal to at least $z=0.835$, measured for GRB970508. Despite this result, the data from $CGRO$/BATSE and $Swift$/BAT indicate that long GRBs are indeed distributed isotropically, but the distribution of short GRBs is anisotropic. $Fermi$/GBM has detected 1669 GRBs up to date, and their sky distribution is examined in this paper. A number of statistical tests is applied: nearest neighbour analysis, fractal dimension, dipole and quadrupole moments of the distribution function decomposed into spherical harmonics, binomial test, and the two point angular correlation function. Monte Carlo benchmark testing of each test is performed in order to evaluate its reliability. It is found that short GRBs are distributed anisotropically on the sky, and long ones have an isotropic distribution. The probability that these results are not a chance occurence is equal to at least 99.98\% and 30.68\% for short and long GRBs, respectively. The cosmological context of this finding and its relation to large-scale structures is discussed.

The two widely accepted classes of gamma-ray bursts (GRBs), short and long, are with confidence ascribed to mergers of compact objects and collapse of massive stars, respectively. A third, intermediate/soft class, remains putative. Its existence was claimed based on univariate and bivariate analyses of GRB observables modeled with Gaussian distributions. This, however, may not be the appropriate approach, as it has been already shown that the univariate distributions of durations are better described by mixtures of two skewed components rather than three Gaussian ones. This paper investigates whether data in the duration--hardness ratio plane is better modeled by mixtures of skewed bivariate distributions than by normal ones. The archival data set of the Compton Gamma-Ray Observatory/BATSE and Fermi/GBM data from the most recent catalogue release are examined. The preferred model is chosen based on two information criteria, Akaike ($AIC$) and Bayesian ($BIC$). It is found that the best description is given by a two-component mixture of skewed Student-$t$ distributions, which outperforms any other model considered. This implies that the distribution of the studied parameters is intrinsically skewed, introducing spurious Gaussian components, and hence the third class is unlikely to be a real phenomenon. Its existence, based on statistical inference, is therefore rejected as unnecessary to explain the observations.

The mechanism responsible for the prompt emission of gamma-ray bursts (GRBs) is still a debated issue. The prompt phase-related GRB correlations can allow to discriminate among the most plausible theoretical models explaining this emission. We present an overview of the observational two-parameter correlations, their physical interpretations, their use as redshift estimators and possibly as cosmological tools. The nowadays challenge is to make GRBs, the farthest stellar-scaled objects observed (up to redshift $z=9.4$), standard candles through well established and robust correlations. However, GRBs spanning several orders of magnitude in their energetics are far from being standard candles. We describe the advances in the prompt correlation research in the past decades, with particular focus paid to the discoveries in the last 20 years.

Aims. This paper investigates the chaotic rotation of an oblate satellite in the context of chaos control. Methods. A model of planar oscillations, described with the Beletskii equation, was investigated. The Hamiltonian formalism was utilized to employ a control method for suppressing chaos. Results. An additive control term, which is an order of magnitude smaller than the potential, is constructed. This allows not only for significantly diminished diffusion of the trajectory in the phase space, but turns the purely chaotic motion into strictly periodic motion.

This paper presents bifurcation and generalized bifurcation diagrams for a rotational model of an oblate satellite. Special attention is paid to parameter values describing one of Saturn's moons, Hyperion. For various oblateness the largest Lyapunov Characteristic Exponent (LCE) is plotted. The largest LCE in the initial condition as well as in the mixed parameter-initial condition space exhibits a fractal structure, for which the fractal dimension was calculated. It results from the bifurcation diagrams of which most of the parameter values for preselected initial conditions lead to chaotic rotation. The First Recurrence Time (FRT) diagram provides an explanation of the birth of chaos and the existence of quasi-periodic windows occuring in the bifurcation diagrams.

The gravitational influence of a second satellite on the rotation of an oblate moon is numerically examined. A simplified model, assuming the axis of rotation perpendicular to the (Keplerian) orbit plane, is derived. The differences between the two models, i.e. in the absence and presence of the second satellite, are investigated via bifurcation diagrams and by evolving compact sets of initial conditions in the phase space. It turns out that the presence of another satellite causes some trajectories, that were regular in its absence, to become chaotic. Moreover, the highly structured picture revealed by the bifurcation diagrams in dependence on the eccentricity of the oblate body's orbit is destroyed when the gravitational influence is included, and the periodicities and critical curves are destroyed as well. For demonstrative purposes, focus is laid on parameters of the Saturn-Titan-Hyperion system, and on oblate satellites on low-eccentric orbits, i.e. $e\approx 0.005$.

Two classes of GRBs have been confidently identified thus far and are prescribed to different physical scenarios -- NS-NS or NS-BH mergers, and collapse of massive stars, for short and long GRBs, respectively. A third, intermediate in duration class, was suggested to be present in previous catalogs, such as BATSE and $Swift$, based on statistical tests regarding a mixture of two or three log-normal distributions of $T_{90}$. However, this might possibly not be an adequate model. This paper investigates whether the distributions of $\log T_{90}$ from BATSE, $Swift$, and $Fermi$ are described better by a mixture of skewed distributions rather than standard Gaussians. Mixtures of standard normal, skew-normal, sinh-arcsinh and alpha-skew-normal distributions are fitted using a maximum likelihood method. The preferred model is chosen based on the Akaike information criterion. It is found that mixtures of two skew-normal or two sinh-arcsinh distributions are more likely to describe the observed duration distribution of $Fermi$ than a mixture of three standard Gaussians, and that mixtures of two sinh-arcsinh or two skew-normal distributions are models competing with the conventional three-Gaussian in the case of BATSE and $Swift$. Based on statistical reasoning, existence of a third (intermediate) class of GRBs in $Fermi$ data is rejected, and it is shown that other phenomenological models may describe the observed $Fermi$, BATSE, and $Swift$ duration distributions at least as well as a mixture of standard normal distributions.

The duration distribution of 408 GRBs with measured both duration $T_{90}$ and redshift $z$ is examined. Mixtures of a number of distributions (standard normal, skew-normal, sinh-arcsinh, and alpha-skew-normal) are fitted to the observed and intrinsic durations using the maximum log-likelihood method. The best fit is chosen via the Akaike information critetion. The aim of this work is to assess the presence of the presumed intermediate GRB class, and to provide a phenomenological model more appropriate than the common mixture of standard Gaussians. While $\log T^{obs}_{90}$ are well described by a truly trimodal fit, after moving to the rest frame the statistically most significant fit is unimodal. To trace the source of this discrepancy, 334 GRBs observed only by $Swift$/BAT are examined in the same way. In the observer frame, this results in a number of statistically plausible descriptions, being uni- and bimodal, and with the number of components ranging from one to three. After moving to the rest frame, no unambiguous conclusions may be put forward. It is concluded that the size of the sample is not big enough to infer reliably GRB properties based on a univariate statistical reasoning only.

The duration distribution of 947 GRBs observed by $Swift$/BAT, as well as its subsample of 347 events with measured redshift, allowing to examine the durations in both the observer and rest frames, are examined. Using a maximum log-likelihood method, mixtures of two and three standard Gaussians are fitted to each sample, and the adequate model is chosen based on the value of the difference in the log-likelihoods, Akaike information criterion and Bayesian information criterion. It is found that a two-Gaussian is a better description than a three-Gaussian, and that the presumed intermediate-duration class is unlikely to be present in the $Swift$ duration data.

Two classes of gamma-ray bursts (GRBs), short and long, have been determined without any doubts, and are usually ascribed to different progenitors, yet these classes overlap for a variety of descriptive parameters. A subsample of 46 long and 22 short $Fermi$ GRBs with estimated Hurst Exponents (HEs), complemented by minimum variability time-scales (MVTS) and durations ($T_{90}$) is used to perform a supervised Machine Learning (ML) and Monte Carlo (MC) simulation using a Support Vector Machine (SVM) algorithm. It is found that while $T_{90}$ itself performs very well in distinguishing short and long GRBs, the overall success ratio is higher when the training set is complemented by MVTS and HE. These results may allow to introduce a new (non-linear) parameter that might provide less ambiguous classification of GRBs.

Two classes of GRBs have been identified thus far without doubt and are prescribed to different physical scenarios -- NS-NS or NS-BH mergers, and collapse of massive stars, for short and long GRBs, respectively. The existence of two distinct populations was inferred through a bimodal distribution of the observed durations $T_{90}$, and the commonly applied $2\,{\rm s}$ limit between short and long GRBs was obtained by fitting a parabola between the two peaks in binned data from BATSE 1B. Herein, by means of a maximum likelihood (ML) method a mixture of two Gaussians is fitted to the datasets from BATSE, $Swift$, $BeppoSAX$, and $Fermi$ in search for a local minimum that might serve as a new, more proper, limit for the two GRB classes. It is found that $Swift$ and $BeppoSAX$ distributions are unimodal, hence no local minimum is present, $Fermi$ is consistent with the conventional limit, whereas BATSE gives the limit significantly longer (equal to $3.38\pm 0.27\,{\rm s}$) than $2\,{\rm s}$. These new values change the fractions of short and long GRBs in the samples examined, and imply that the observed $T_{90}$ durations are detector dependent, hence no universal limiting value may be applied to all satellites due to their different instrument specifications. Because of this, and due to the strong overlap of the two-Gaussian components, the straightforward association of short GRBs to mergers and long ones to collapsars is ambiguous.

Two classes of GRBs, short and long, have been determined without any doubts, and are usually prescribed to different physical scenarios. A third class, intermediate in $T_{90}$ durations, has been reported to be present in the datasets of BATSE, Swift, RHESSI and possibly BeppoSAX. The latest release of $>1500$ GRBs observed by Fermi gives an opportunity to further investigate the duration distribution. The aim of this paper is to investigate whether a third class is present in the $\log T_{90}$ distribution, or is it described by a bimodal distribution. A standard $\chi^2$ fitting of a mixture of Gaussians is applied to 25 histograms with different binnings. Different binnings give various values of the fitting parameters, as well as the shape of the fitted curve. Among five statistically significant fits none is trimodal. Locations of the Gaussian components are in agreement with previous works. However, a trimodal distribution, understood in the sense of having three separated peaks, is not found for any binning. It is concluded that the duration distribution in Fermi data is well described by a mixture of three log-normal distributions, but it is intrinsically bimodal, hence no third class is present in the $T_{90}$ data of Fermi. It is suggested that the log-normal fit may not be an adequate model.

Hyperion is a satellite of Saturn that was predicted to remain in a chaotic rotational state. This was confirmed to some extent by Voyager 2 and Cassini series of images and some ground-based photometric observations. The aim of this aticle is to explore conditions for potential observations to meet in order to estimate a maximal Lyapunov Exponent (mLE), which being positive is an indicator of chaos and allows to characterise it quantitatively. Lightcurves existing in literature as well as numerical simulations are examined using standard tools of theory of chaos. It is found that existing datasets are too short and undersampled to detect a positive mLE, although its presence is not rejected. Analysis of simulated lightcurves leads to an assertion that observations from one site should be performed over a year-long period to detect a positive mLE, if present, in a reliable way. Another approach would be to use 2---3 telescopes spread over the world to have observations distributed more uniformly. This may be achieved without disrupting other observational projects being conducted. The necessity of time-series to be stationary is highly stressed.

The case of Hyperion has been studied excesively due to the fact it is the largest known celestial body of a highly aspherical shape. It also has a low mass density and remains in a 4:3 orbital resonance with Titan. Its lightcurve, obtained through photometric observations by (Klavetter 1989a,b), was initialy used to show that Hyperion's rotation exhibits no periodicity. Herein it is analyzed in the means of time series analysis. The Hurst Exponent was estimated to be H=0.87, indicating a persistent behaviour. The largest Lyapunov Exponent $\lambda_{max}$ unfortunately could not be given a reliable estimate because of the shortness of the dataset, consisting 38 observational points. These results are compared with numerical simulations, which gave a value H=0.88 for the chaotic zone of the phase space. The Lyapunov time $T_{Lyap}=1/\lambda_{max}$ is about 30 days, which is roughly 1.5 times greater than the orbital period. By conducting observations over a longer period an insight in the dynamical features of the present rotational state is possible.