BZPEERReading, trust and collaboration for researchDiscover
Workspaces
Network
Opportunities
Account
Researcher profile
Mariusz Tarnopolski contributes to research discovery and scholarly infrastructure.
Trust snapshot
Actions
Identity and collaboration
Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.
Log in to claimDirect collaboration
Claim this author entity first to unlock direct invitations.
Research graph
Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
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.
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.
Closed-form expressions, parametrized by the Hurst exponent $H$ and the length $n$ of a time series, are derived for paths of fractional Brownian motion (fBm) and fractional Gaussian noise (fGn) in the $\mathcal{A}-\mathcal{T}$ plane, composed of the fraction of turning points $\mathcal{T}$ and the Abbe value $\mathcal{A}$. The exact formula for $\mathcal{A}_{\rm fBm}$ is expressed via Riemann $ζ$ and Hurwitz $ζ$ functions. A very accurate approximation, yielding a simple exponential form, is obtained. Finite-size effects, introduced by the deviation of fGn's variance from unity, and asymptotic cases are discussed. Expressions for $\mathcal{T}$ for fBm, fGn, and differentiated fGn are also presented. The same methodology, valid for any Gaussian process, is applied to autoregressive moving average processes, for which regions of availability of the $\mathcal{A}-\mathcal{T}$ plane are derived and given in analytic form. Locations in the $\mathcal{A}-\mathcal{T}$ plane of some real-world examples as well as generated data are discussed for illustration.
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.
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.
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 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.
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.
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 $χ^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.
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.
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.
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.
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 $λ_{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/λ_{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.
The box counting dimension $d_C$ and the correlation dimension $d_G$ change with the number of numerically generated points forming the attractor. At a sufficiently large number of points the fractal dimension tends to a finite value. The obtained values are $d_C\approx 1.43$ and $d_G\approx 1.38$.