Critical transitions occur in complex dynamical systems, when the system dynamics undergoes a regime shift. These can often occur with little change in the mean amplitude of system response prior to the actual time of transition. The recent dimming and brightening event in Betelgeuse occured as a sudden shift in the brightness and has been the subject of much debate. Internal changes or an external dust cloud have been suggested as reasons for this change in variability. We examine whether the dimming and brightening event of 2019-20 could be due to a critical transition in the pulsation dynamics of Betelgeuse, by studying the characteristics of the light curve prior to transition. We calculate the quantifiers hypothesised to rise prior to a critical transition for the light curve of Betelgeuse upto the dimming event of 2019-20. These include the autocorrelation at lag-1, variance and the spectral coefficient calculated from detrended fluctation analysis (DFA), apart from two measures that quantify the recurrence properties of the light curve. Significant rises are confirmed using the Mann-Kendall trend test. We see a significant increase in all quantifiers (p < 0.05) prior to the dimming event of 2019-20. This suggests that the event was a critical transition related to the underlying nonlinear dynamics of the star. Together with results that suggests minimal change in $T_{eff}$ and infra-red flux, a critical transition in the pulsation dynamics could be a possible reason for the unprecedented dimming of Betelgeuse. The rise in the studied quantifiers prior to the dimming event, supports this possibility.

Close binary stars are binary stars where the component stars are close enough such that they can exchange mass and/or energy. They are subdivided into semi-detached, overcontact or ellipsoidal binary stars. A challenging problem in the context of close binary stars, is their classification into these subclasses, based solely on their light curves. Conventionally, this is done by observing subtle features in the light curves like the depths of adjacent minima, which is tedious when dealing with large datasets. In this work we suggest the use of machine learning algorithms applied to measures of recurrence networks and nonlinear time series analysis to differentiate between classes of close binary stars. We show that overcontact binary stars occupy a region different from semi-detached and ellipsoidal binary stars in a plane of characteristic path length(CPL) and average clustering coefficient(CC), computed from their recurrence networks. We use standard clustering algorithms and report that the clusters formed corresponds to the standard classes with a high degree of accuracy.

Overcontact binary stars are systems of two stars where the component stars are in contact with each other. This implies that they share a common envelope of gas. In this work we seek signatures of nonlinearity and chaos in these stars by using time series analysis techniques. We use three main techniques, namely the correlation dimension,f (\alpha) spectrum and the bicoherence. The former two are calculated from the reconstructed dynamics, while the latter is calculated from the Fourier transforms of the time series of intensity variations(light curves) of these stars. Our dataset consists of data from 463 overcontact binary stars in the Kepler field of view [1]. Our analysis indicates nonlinearity and signatures of chaos in almost all the light curves. We also explore whether the underlying nonlinear properties of the stars are related to their physical properties like fill-out-factor, a measure of the extend of contact between the components of an overcontact binary system . We observe that significant correlations exist between the fill out factor and the nonlinear quantifiers. This correlation is more pronounced in specific subcategories constructed based on the mass ratios and effective temperatures of the binaries. The correlations observed can be indicative of variations in the nonlinear properties of the star as it ages. We believe that this study relating nonlinear and astrophysical properties of binary stars is the first of its kind and is an important starting point for such studies in other astrophysical objects displaying nonlinear dynamical behaviour.

Recurrence networks and the associated statistical measures have become important tools in the analysis of time series data. In this work, we test how effective the recurrence network measures are in analyzing real world data involving two main types of noise, white noise and colored noise. We use two prominent network measures as discriminating statistic for hypothesis testing using surrogate data for a specific null hypothesis that the data is derived from a linear stochastic process. We show that the characteristic path length is especially efficient as a discriminating measure with the conclusions reasonably accurate even with limited number of data points in the time series. We also highlight an additional advantage of the network approach in identifying the dimensionality of the system underlying the time series through a convergence measure derived from the probability distribution of the local clustering coefficients. As examples of real world data, we use the light curves from a prominent black hole system and show that a combined analysis using three primary network measures can provide vital information regarding the nature of temporal variability of light curves from different spectroscopic classes.

Deriving meaningful information from observational data is often restricted by many limiting factors, the most important of which is the presence of noise. In this work, we present the use of the bicoherence function to extract information about the underlying nonlinearity from noisy time series. We show that a system evolving in the presence of noise which has its dynamical state concealed from quantifiers like the power spectrum and correlation dimension D2, can be revealed using the bicoherence function. We define an index called main peak bicoherence function as the bicoherence associated with the maximal power spectral peak. We show that this index is extremely useful while dealing with quasi-periodic data as it can distinguish strange non chaos from quasi periodicity even with added noise. We demonstrate this in a real world scenario, by taking the bicoherence of variable stars showing period doubling and strange non-chaotic behavior. Our results indicate that bicoherence analysis can also bypass the method of surrogate analysis using Fourier phase randomization, used to differentiate linear stochastic processes from non linear ones, in conventional methods involving measures like D2.

Observational data, especially astrophysical data, is often limited by gaps in data that arises due to lack of observations for a variety of reasons. Such inadvertent gaps are usually smoothed over using interpolation techniques. However the smoothing techniques can introduce artificial effects, especially when non-linear analysis is undertaken. We investigate how gaps can affect the computed values of correlation dimension of the system, without using any interpolation. For this we introduce gaps artificially in synthetic data derived from standard chaotic systems, like the R{\"o}ssler and Lorenz, with frequency of occurrence and size of missing data drawn from two Gaussian distributions. Then we study the changes in correlation dimension with change in the distributions of position and size of gaps. We find that for a considerable range of mean gap frequency and size, the value of correlation dimension is not significantly affected, indicating that in such specific cases, the calculated values can still be reliable and acceptable. Thus our study introduces a method of checking the reliability of computed correlation dimension values by calculating the distribution of gaps with respect to its size and position. This is illustrated for the data from light curves of three variable stars, R Scuti, U Monocerotis and SU Tauri. We also demonstrate how a cubic spline interpolation can cause a time series of Gaussian noise with missing data to be misinterpreted as being chaotic in origin. This is demonstrated for the non chaotic light curve of variable star SS Cygni, which gives a saturated D$_{2}$ value, when interpolated using a cubic spline. In addition we also find that a careful choice of binning, in addition to reducing noise, can help in shifting the gap distribution to the reliable range for D$_2$ values.

GRS 1915+105 is a prominent black hole system exhibiting variability over a wide range of time scales and its observed light curves have been classified into 12 temporal states. Here we undertake a complete analysis of these light curves from all the states using various quantifiers from nonlinear time series analysis, such as, the correlation dimension (D_2), the correlation entropy (K_2), singular value decomposition (SVD) and the multifractal spectrum ($f(\alpha)$ spectrum). An important aspect of our analysis is that, for estimating these quantifiers, we use algorithmic schemes which we have proposed recently and tested successfully on synthetic as well as practical time series from various fields. Though the schemes are based on the conventional delay embedding technique, they are automated so that the above quantitative measures can be computed using conditions prescribed by the algorithm and without any intermediate subjective analysis. We show that nearly half of the 12 temporal states exhibit deviation from randomness and their complex temporal behavior could be approximated by a few (3 or 4) coupled ordinary nonlinear differential equations. These results could be important for a better understanding of the processes that generate the light curves and hence for modelling the temporal behavior of such complex systems. To our knowledge, this is the first complete analysis of an astrophysical object (let alone a black hole system) using various techniques from nonlinear dynamics.

We present an adaptation of the standard Grassberger-Proccacia (GP) algorithm for estimating the Correlation Dimension of a time series in a non subjective manner. The validity and accuracy of this approach is tested using different types of time series, such as, those from standard chaotic systems, pure white and colored noise and chaotic systems added with noise. The effectiveness of the scheme in analysing noisy time series, particularly those involving colored noise, is investigated. An interesting result we have obtained is that, for the same percentage of noise addition, data with colored noise is more distinguishable from the corresponding surrogates, than data with white noise. As examples for real life applications, analysis of data from an astrophysical X-ray object and human brain EEG, are presented.

Using non-linear time series analysis, along with surrogate data analysis, it is shown that the various types of long term variability exhibited by the black hole system GRS 1915+105, can be explained in terms of a deterministic non-linear system with some inherent stochastic noise. Evidence is provided for a non-linear limit cycle origin of one of the low frequency QPO detected in the source, while some other types of variability could be due to an underlying low dimensional chaotic system. These results imply that the partial differential equations which govern the magneto-hydrodynamic flow of the inner accretion disk, can be approximated by a small number ($\approx 3 -5$) of non-linear but {\it ordinary} differential equations. While this analysis does not reveal the exact nature of these approximate equations, they may be obtained in the future, after results of magneto-hydrodynamic simulation of realistic accretion disks become available.

A modified non-linear time series analysis technique, which computes the correlation dimension $D_2$, is used to analyze the X-ray light curves of the black hole system GRS 1915+105 in all twelve temporal classes. For four of these temporal classes $D_2 $ saturates to $\approx 4-5$ which indicates that the underlying dynamical mechanism is a low dimensional chaotic system. Of the other eight classes, three show stochastic behavior while five show deviation from randomness. The light curves for four classes which depict chaotic behavior have the smallest ratio of the expected Poisson noise to the variability ($ < 0.05$) while those for the three classes which depict stochastic behavior is the highest ($ > 0.2$). This suggests that the temporal behavior of the black hole system is governed by a low dimensional chaotic system, whose nature is detectable only when the Poisson fluctuations are much smaller than the variability.