Reliability of energy landscape analysis of resting-state functional MRI data.

Energy landscape analysis is a data-driven method to analyse multidimensional time series, including functional magnetic resonance imaging (fMRI) data. It has been shown to be a useful characterization of fMRI data in health and disease. It fits an Ising model to the data and captures the dynamics of the data as movement of a noisy ball constrained on the energy landscape derived from the estimated Ising model. In the present study, we examine test-retest reliability of the energy landscape analysis. To this end, we construct a permutation test that assesses whether or not indices characterizing the energy landscape are more consistent across different sets of scanning sessions from the same participant (i.e. within-participant reliability) than across different sets of sessions from different participants (i.e. between-participant reliability). We show that the energy landscape analysis has significantly higher within-participant than between-participant test-retest reliability with respect to four commonly used indices. We also show that a variational Bayesian method, which enables us to estimate energy landscapes tailored to each participant, displays comparable test-retest reliability to that using the conventional likelihood maximization method. The proposed methodology paves the way to perform individual-level energy landscape analysis for given data sets with a statistically controlled reliability.

State-transition dynamics of resting-state functional magnetic resonance imaging data: model comparison and test-to-retest analysis.

Electroencephalogram (EEG) microstate analysis entails finding dynamics of quasi-stable and generally recurrent discrete states in multichannel EEG time series data and relating properties of the estimated state-transition dynamics to observables such as cognition and behavior. While microstate analysis has been widely employed to analyze EEG data, its use remains less prevalent in functional magnetic resonance imaging (fMRI) data, largely due to the slower timescale of such data. In the present study, we extend various data clustering methods used in EEG microstate analysis to resting-state fMRI data from healthy humans to extract their state-transition dynamics. We show that the quality of clustering is on par with that for various microstate analyses of EEG data. We then develop a method for examining test-retest reliability of the discrete-state transition dynamics between fMRI sessions and show that the within-participant test-retest reliability is higher than between-participant test-retest reliability for different indices of state-transition dynamics, different networks, and different data sets. This result suggests that state-transition dynamics analysis of fMRI data could discriminate between different individuals and is a promising tool for performing fingerprinting analysis of individuals.

Dimension reduction in higher-order contagious phenomena.

We investigate epidemic spreading in a deterministic susceptible-infected-susceptible model on uncorrelated heterogeneous networks with higher-order interactions. We provide a recipe for the construction of one-dimensional reduced model (resilience function) of the N-dimensional susceptible-infected-susceptible dynamics in the presence of higher-order interactions. Utilizing this reduction process, we are able to capture the microscopic and macroscopic behavior of infectious networks. We find that the microscopic state of nodes (fraction of stable healthy individual of each node) inversely scales with their degree, and it becomes diminished due to the presence of higher-order interactions. In this case, we analytically obtain that the macroscopic state of the system (fraction of infectious or healthy population) undergoes abrupt transition. Additionally, we quantify the network's resilience, i.e., how the topological changes affect the stable infected population. Finally, we provide an alternative framework of dimension reduction based on the spectral analysis of the network, which can identify the critical onset of the disease in the presence or absence of higher-order interactions. Both reduction methods can be extended for a large class of dynamical models.

Amplification of explosive width in complex networks.

We present an adaptive coupling strategy to induce hysteresis/explosive synchronization in complex networks of phase oscillators (Sakaguchi-Kuramoto model). The coupling strategy ensures explosive synchronization with significant explosive width enhancement. Results show the robustness of the strategy, and the strategy can diminish (by inducing enhanced hysteresis loop) the contrarian impact of phase frustration in the network, irrespective of the network structure or frequency distributions. Additionally, we design a set of frequency for the oscillators, which eventually ensure complete in-phase synchronization behavior among these oscillators (with enhanced explosive width) in the case of adaptive-coupling scheme. Based on a mean-field analysis, we develop a semi-analytical formalism, which can accurately predict the backward transition of the synchronization order parameter.

Reliability of energy landscape analysis of resting-state functional MRI data.

Energy landscape analysis is a data-driven method to analyze multidimensional time series, including functional magnetic resonance imaging (fMRI) data. It has been shown to be a useful characterization of fMRI data in health and disease. It fits an Ising model to the data and captures the dynamics of the data as movement of a noisy ball constrained on the energy landscape derived from the estimated Ising model. In the present study, we examine test-retest reliability of the energy landscape analysis. To this end, we construct a permutation test that assesses whether or not indices characterizing the energy landscape are more consistent across different sets of scanning sessions from the same participant (i.e., within-participant reliability) than across different sets of sessions from different participants (i.e., between-participant reliability). We show that the energy landscape analysis has significantly higher within-participant than between-participant test-retest reliability with respect to four commonly used indices. We also show that a variational Bayesian method, which enables us to estimate energy landscapes tailored to each participant, displays comparable test-retest reliability to that using the conventional likelihood maximization method. The proposed methodology paves the way to perform individual-level energy landscape analysis for given data sets with a statistically controlled reliability.

Perfect synchronization in complex networks with higher-order interactions.

Achieving perfect synchronization in a complex network, specially in the presence of higher-order interactions (HOIs) at a targeted point in the parameter space, is an interesting, yet challenging task. Here we present a theoretical framework to achieve the same under the paradigm of the Sakaguchi-Kuramoto (SK) model. We analytically derive a frequency set to achieve perfect synchrony at some desired point in a complex network of SK oscillators with higher-order interactions. Considering the SK model with HOIs on top of the scale-free, random, and small world networks, we perform extensive numerical simulations to verify the proposed theory. Numerical simulations show that the analytically derived frequency set not only provides stable perfect synchronization in the network at a desired point but also proves to be very effective in achieving a high level of synchronization around it compared to the other choices of frequency sets. The stability and the robustness of the perfect synchronization state of the system are determined using the low-dimensional reduction of the network and by introducing a Gaussian noise around the derived frequency set, respectively.

Impact of phase lag on synchronization in frustrated Kuramoto model with higher-order interactions.

The study of first order transition (explosive synchronization) in an ensemble (network) of coupled oscillators has been the topic of paramount interest among the researchers for more than one decade. Several frameworks have been proposed to induce explosive synchronization in a network and it has been reported that phase frustration in a network usually suppresses first order transition in the presence of pairwise interactions among the oscillators. However, on the contrary, by considering networks of phase frustrated coupled oscillators in the presence of higher-order interactions (up to 2-simplexes) we show here, under certain conditions, phase frustration can promote explosive synchronization in a network. A low-dimensional model of the network in the thermodynamic limit is derived using the Ott-Antonsen ansatz to explain this surprising result. Analytical treatment of the low-dimensional model, including bifurcation analysis, explains the apparent counter intuitive result quite clearly.

Transition to synchrony in degree-frequency correlated Sakaguchi-Kuramoto model.

We investigate transition to synchrony in degree-frequency correlated Sakaguchi-Kuramoto (SK) model on complex networks both analytically and numerically. We analytically derive self-consistent equations for group angular velocity and order parameter for the model in the thermodynamic limit. Using the self-consistent equations we investigate transition to synchronization in SK model on uncorrelated scale-free (SF) and Erdos-Rényi (ER) networks in detail. Depending on the degree distribution exponent (γ) of SF networks and phase-frustration parameter, the population undergoes from first-order transition [explosive synchronization (ES)] to second-order transition and vice versa. In ER networks transition is always second order irrespective of the values of the phase-lag parameter. We observe that the critical coupling strength for the onset of synchronization is decreased by phase-frustration parameter in case of SF network where as in ER network, the phase-frustration delays the onset of synchronization. Extensive numerical simulations using SF and ER networks are performed to validate the analytical results. An analytical expression of critical coupling strength for the onset of synchronization is also derived from the self-consistent equations considering the vanishing order parameter limit.