Modeling non-Poissonian temporal hypergraphs by Markovian node dynamics
Pith reviewed 2026-05-10 18:18 UTC · model grok-4.3
The pith
Markovian node activity switches produce non-Poissonian event timing in temporal hypergraphs
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Despite the Markovian dynamics of node states, the event processes on hyperedges emerge as mixtures of Poissonian components. This mixture arises because the effective rate for each hyperedge fluctuates according to the instantaneous number of high-activity nodes it contains, producing interevent time distributions with heavier tails and autocorrelation functions that decay more slowly than those of homogeneous Poisson processes. The strength of these deviations increases with hyperedge size, in agreement with empirical temporal hypergraph data.
What carries the argument
Two-state Markovian node dynamics combined with hyperedge event probability that depends on the number of high-state nodes in the hyperedge.
If this is right
- The inter-event time distribution for a hyperedge is a weighted sum of exponential distributions, with weights given by the binomial probabilities of having k high nodes.
- Autocorrelation functions of event sequences decay as a sum of exponentials whose rates are set by the node transition rates.
- For larger hyperedges the effective event process approaches a Poisson process with average rate because the fraction of high nodes concentrates around its mean.
- Node-level inter-event times remain shorter-tailed than the induced hyperedge-level times under the same dynamics.
Where Pith is reading between the lines
- Testing the model on hypergraphs with heterogeneous hyperedge sizes could reveal whether size-dependent burstiness holds across domains.
- If node states are not fully independent the mixture may produce even longer tails, suggesting a direction for model refinement.
- The same mechanism might apply to modeling temporal networks with higher-order interactions beyond hypergraphs.
- One could use the analytical expressions to fit node transition rates directly from observed group event data.
Load-bearing premise
The event generation rate for each hyperedge depends only on the current count of high-activity nodes within it and not on any history or external factors.
What would settle it
Collect inter-event time histograms for hyperedges of varying sizes from empirical temporal hypergraph data and check whether the tail heaviness or coefficient of variation scales with size exactly as predicted by the mixture-of-Poissons formula.
Figures
read the original abstract
Temporal hypergraphs capture time-resolved group interactions among nodes. Empirical data support that time-stamped group interactions show bursty event sequences and non-trivial temporal correlations. In the present study, we introduce node-driven temporal hypergraph models in which each node stochastically alternates between low- and high-activity states, and a hyperedge produces time-stamped events with a probability that depends on the number of high-state nodes in the hyperedge. For two event-generation rules, we analytically derive interevent time distributions and autocorrelation functions of event sequences, both for hyperedges and nodes. Despite Markovian node-state dynamics, the induced event processes become mixtures of Poissonian, short-tailed components, resulting in longer-tailed interevent time distributions and slowly decaying autocorrelation. The theory further shows the dependence of these features on the size of hyperedge, which largely agrees with various empirical data. We expect our models to provide a simple, interpretable framework for connecting individual-level activity fluctuations to the timing patterns observed in real group interactions.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper introduces node-driven models for temporal hypergraphs in which each node follows an independent continuous-time Markov chain alternating between low- and high-activity states. A hyperedge generates events according to one of two deterministic rules that depend only on the instantaneous count of high-activity nodes inside it. For both rules the authors derive closed-form expressions for the inter-event time distribution and the autocorrelation function of the event sequence, both at the hyperedge and at the node level. The resulting aggregate processes are mixtures of Poissonian components whose tails are longer than exponential and whose autocorrelations decay slowly; both features strengthen with hyperedge size, in qualitative agreement with several empirical data sets.
Significance. If the derivations are correct, the work supplies a minimal, fully analytical mechanism that converts Markovian individual dynamics into the bursty, long-range-correlated event sequences observed in real group interactions. The explicit dependence on hyperedge cardinality furnishes testable predictions and connects the model to the broader literature on how microscopic fluctuations generate macroscopic temporal heterogeneity.
major comments (2)
- [§3.2] §3.2, Eq. (8)–(11): the master equation for the count process N(t) is written, but the explicit transition rates q_{k→k±1} induced by the independent node chains are not displayed; without them the claim that the event process is exactly a mixture of Poisson processes cannot be verified by the reader.
- [§5.3] §5.3, Fig. 4: the comparison of model predictions with empirical inter-event time distributions is visual only; no quantitative goodness-of-fit measure (KS statistic, likelihood ratio, or parameter uncertainty) is reported, weakening the statement that the size dependence “largely agrees” with data.
minor comments (4)
- The two event-generation rules are introduced in §2.2 but never given short mnemonic names; consistent labels (e.g., “linear rule” and “threshold rule”) would improve readability throughout the derivations and figures.
- [§4] Notation for the stationary probability π_k of k active nodes inside a hyperedge of size m is introduced only in §3.1; repeating the definition once in the results section would help readers who skip the methods.
- [§1] Several sentences in the introduction cite “temporal hypergraphs” without distinguishing between the static hypergraph and its temporal realization; a brief clarifying sentence would avoid ambiguity.
- The supplementary material is referenced for “full derivations,” yet the main text does not indicate which equations are proved there versus derived in place; an explicit pointer at the end of each derivation would be helpful.
Simulated Author's Rebuttal
We thank the referee for the careful reading, positive assessment, and constructive suggestions. We address each major comment below and will revise the manuscript accordingly.
read point-by-point responses
-
Referee: [§3.2] §3.2, Eq. (8)–(11): the master equation for the count process N(t) is written, but the explicit transition rates q_{k→k±1} induced by the independent node chains are not displayed; without them the claim that the event process is exactly a mixture of Poisson processes cannot be verified by the reader.
Authors: We agree that the explicit transition rates q_{k→k±1} should be displayed to allow direct verification of the mixture-of-Poisson property. In the revised manuscript we will derive and insert these rates from the independent node-level continuous-time Markov chains, showing how they enter the master equation for the count process N(t) and confirming that the resulting event process is indeed a mixture of Poisson processes whose weights depend on hyperedge size. revision: yes
-
Referee: [§5.3] §5.3, Fig. 4: the comparison of model predictions with empirical inter-event time distributions is visual only; no quantitative goodness-of-fit measure (KS statistic, likelihood ratio, or parameter uncertainty) is reported, weakening the statement that the size dependence “largely agrees” with data.
Authors: The referee correctly notes that a purely visual comparison limits the strength of the agreement claim. In the revision we will augment Fig. 4 and the accompanying text with quantitative measures: Kolmogorov-Smirnov statistics comparing model and empirical inter-event time distributions for each hyperedge size, together with the fitted parameter values and their uncertainties obtained from the maximum-likelihood procedure. These additions will provide an objective basis for the reported size dependence. revision: yes
Circularity Check
Derivation is self-contained from model definitions
full rationale
The paper starts from explicit definitions of independent continuous-time Markov chains for each node's activity state (low/high) and two deterministic event-rate rules that depend only on the instantaneous count of high-activity nodes within a hyperedge. Inter-event time distributions and autocorrelation functions are then obtained by solving the resulting higher-level CTMC for the count process and applying standard renewal theory or embedded Markov equations. These steps use only the model primitives and binomial stationary occupancy; no fitted parameters are renamed as predictions, no self-citations supply load-bearing uniqueness theorems, and no ansatz is smuggled in. The size dependence of tail weight and autocorrelation decay follows directly from the hyperedge-size parameter in the binomial count distribution. The derivation therefore contains no circular reduction.
Axiom & Free-Parameter Ledger
free parameters (2)
- node state transition rates
- event generation probability parameters
axioms (2)
- domain assumption Node dynamics are independent Markov chains
- domain assumption Event probability depends only on count of high-state nodes
Reference graph
Works this paper leans on
- [1]
-
[2]
Holme, The European Physical Journal B88, 1 (2015)
P. Holme, The European Physical Journal B88, 1 (2015)
work page 2015
-
[3]
P. Holme and J. Saram¨ aki, eds.,Temporal Network The- ory, Computational Social Sciences (Springer Interna- tional Publishing, Cham, 2019). 10
work page 2019
-
[4]
N. Masuda and R. Lambiotte,A Guide to Temporal Net- works, 2nd ed., Series on Complexity Science No. vol. 6 (World Scientific, Singapore, 2020)
work page 2020
-
[5]
L. Kovanen, M. Karsai, K. Kaski, J. Kert´ esz, and J. Saram¨ aki, Journal of Statistical Mechanics: Theory and Experiment2011, P11005 (2011)
work page 2011
-
[6]
A. Paranjape, A. R. Benson, and J. Leskovec, inProceed- ings of the Tenth ACM International Conference on Web Search and Data Mining, WSDM ’17 (ACM, Cambridge, United Kingdom, 2017) pp. 601–610
work page 2017
-
[7]
P. Liu, V. Guarrasi, and A. E. Sarıy¨ uce, IEEE Trans- actions on Knowledge and Data Engineering35, 945 (2023)
work page 2023
- [8]
- [9]
- [10]
-
[11]
M. Starnini, A. Baronchelli, A. Barrat, and R. Pastor- Satorras, Physical Review E85, 056115 (2012)
work page 2012
- [12]
- [13]
- [14]
- [15]
-
[16]
P. Van Mieghem and R. Van De Bovenkamp, Physical Review Letters110, 108701 (2013)
work page 2013
-
[17]
M. Rosvall, A. V. Esquivel, A. Lancichinetti, J. D. West, and R. Lambiotte, Nature Communications5, 4630 (2014)
work page 2014
-
[18]
N. Masuda and P. Holme, inTemporal Network Epidemi- ology, edited by N. Masuda and P. Holme (Springer Sin- gapore, Singapore, 2017) pp. 1–16
work page 2017
-
[19]
A. Li, S. P. Cornelius, Y. Y. Liu, L. Wang, and A. L. Barab´ asi, Science358, 1042 (2017)
work page 2017
-
[20]
R. Lambiotte, M. Rosvall, and I. Scholtes, Nature Physics 15, 313 (2019)
work page 2019
-
[21]
F. Battiston, G. Cencetti, I. Iacopini, V. Latora, M. Lu- cas, A. Patania, J.-G. Young, and G. Petri, Physics Re- ports874, 1 (2020)
work page 2020
-
[22]
F. Battiston, E. Amico, A. Barrat, G. Bianconi, G. Ferraz De Arruda, B. Franceschiello, I. Iacopini, S. K´ efi, V. La- tora, Y. Moreno, M. M. Murray, T. P. Peixoto, F. Vac- carino, and G. Petri, Nature Physics17, 1093 (2021)
work page 2021
-
[23]
Bianconi,Higher-Order Networks: An Introduction to Simplicial Complexes, Cambridge Elements
G. Bianconi,Higher-Order Networks: An Introduction to Simplicial Complexes, Cambridge Elements. Elements in the Structure and Dynamics of Complex Networks (Cam- bridge University Press, Cambridge, 2021)
work page 2021
- [24]
-
[25]
C. Bick, E. Gross, H. A. Harrington, and M. T. Schaub, SIAM Review65, 686 (2023)
work page 2023
-
[26]
G. Ferraz De Arruda, A. Aleta, and Y. Moreno, Nature Reviews Physics6, 468 (2024)
work page 2024
-
[27]
F. Battiston, C. Bick, M. Lucas, A. P. Mill´ an, P. S. Skardal, and Y. Zhang, Nature Reviews Physics8, 146 (2026)
work page 2026
-
[28]
G. Cencetti, F. Battiston, B. Lepri, and M. Karsai, Sci- entific Reports11, 7028 (2021)
work page 2021
- [29]
-
[30]
I. Iacopini, M. Karsai, and A. Barrat, Nature Communi- cations15, 7391 (2024)
work page 2024
-
[31]
M. Mancastroppa, G. Cencetti, and A. Barrat, Physical Review E112, 054305 (2025)
work page 2025
- [32]
-
[33]
A. V´ azquez, J. G. Oliveira, Z. Dezs¨ o, K.-I. Goh, I. Kon- dor, and A.-L. Barab´ asi, Physical Review E73, 036127 (2006)
work page 2006
-
[34]
M. Karsai, H.-H. Jo, and K. Kaski,Bursty Human Dy- namics(Springer International Publishing, Cham, 2018)
work page 2018
- [35]
-
[36]
H. (Bart) Peters, A. Ceria, and H. Wang, PLOS ONE 20, e0323753 (2025)
work page 2025
-
[37]
S. Chowdhary, A. Kumar, G. Cencetti, I. Iacopini, and F. Battiston, Journal of Physics: Complexity2, 035019 (2021)
work page 2021
-
[38]
G. St-Onge, H. Sun, A. Allard, L. H´ ebert-Dufresne, and G. Bianconi, Physical Review Letters127, 158301 (2021)
work page 2021
-
[39]
L. Neuh¨ auser, R. Lambiotte, and M. T. Schaub, Physical Review E104, 064305 (2021)
work page 2021
-
[40]
X. Wang, L. Zhou, A. McAvoy, Z. Tian, and A. Li, Proceedings of the National Academy of Sciences123, e2516380123 (2026)
work page 2026
- [41]
-
[42]
L. Di Gaetano, F. Battiston, and M. Starnini, Physical Review Letters132, 037401 (2024)
work page 2024
-
[43]
Z. Han, L. Liu, X. Wang, Y. Hao, H. Zheng, S. Tang, and Z. Zheng, Chaos: An Interdisciplinary Journal of Nonlinear Science34, 023137 (2024)
work page 2024
-
[44]
M. Mancastroppa, I. Iacopini, G. Petri, and A. Barrat, EPJ Data Science13, 50 (2024)
work page 2024
-
[45]
J. Lerner, M.-G. Hˆ ancean, and M. Perc, Journal of Com- plex Networks13, cnaf054 (2025)
work page 2025
- [46]
-
[47]
C. L. Vestergaard, M. G´ enois, and A. Barrat, Physical Review E90, 042805 (2014)
work page 2014
- [48]
-
[49]
T. Hiraoka, N. Masuda, A. Li, and H.-H. Jo, Physical Review Research2, 023073 (2020)
work page 2020
-
[50]
K.-I. Goh, B. Kahng, and D. Kim, Physical Review Let- ters87, 278701 (2001)
work page 2001
-
[51]
G. Caldarelli, A. Capocci, P. De Los Rios, and M. A. Mu˜ noz, Physical Review Letters89, 258702 (2002)
work page 2002
-
[52]
F. Chung and L. Lu, Proceedings of the National Academy of Sciences99, 15879 (2002)
work page 2002
- [53]
- [54]
-
[55]
R. D. Malmgren, D. B. Stouffer, A. E. Motter, and L. A. N. Amaral, Proceedings of the National Academy of Sciences105, 18153 (2008)
work page 2008
-
[56]
R. D. Malmgren, D. B. Stouffer, A. S. L. O. Campanharo, and L. A. Amaral, Science325, 1696 (2009)
work page 2009
- [57]
- [59]
-
[60]
S. M. Ross,A First Course in Probability, 9th ed. (Pear- son, Boston, 2014)
work page 2014
-
[61]
Yannaros, Annals of the Institute of Statistical Math- ematics46, 641 (1994)
N. Yannaros, Annals of the Institute of Statistical Math- ematics46, 641 (1994)
work page 1994
-
[63]
N. Masuda and C. L. Vestergaard,Gillespie Algorithms for Stochastic Multiagent Dynamics in Populations and Networks(Cambridge University Press, Cambridge, UK, 2022)
work page 2022
-
[64]
K.-I. Goh and A.-L. Barab´ asi, EPL (Europhysics Letters) 81, 48002 (2008)
work page 2008
- [65]
-
[66]
H.-H. Jo, T. Birhanu, and N. Masuda, Chaos: An In- terdisciplinary Journal of Nonlinear Science34, 083110 (2024)
work page 2024
-
[67]
A. R. Benson, Austin R. Benson datasets, https://www.cs.cornell.edu/˜arb/data/ (2025)
work page 2025
-
[68]
A. R. Benson, R. Abebe, M. T. Schaub, A. Jadbabaie, and J. Kleinberg, Proceedings of the National Academy of Sciences115, E11221 (2018)
work page 2018
-
[69]
E. Fonseca dos Reis and N. Masuda, Physical Review Research4, 013050 (2022)
work page 2022
-
[70]
Z. Chen, S. Vijayan, R. Barbieri, M. A. Wilson, and E. N. Brown, Neural Computation21, 1797 (2009)
work page 2009
- [71]
-
[72]
X. Zhao, H. Yu, S. Li, S. Liu, J. Zhang, and X. Cao, Phys- ica A: Statistical Mechanics and its Applications606, 128073 (2022)
work page 2022
-
[73]
L. Zino and A. Rizzo, in2024 IEEE 63rd Conference on Decision and Control (CDC)(IEEE, Milan, Italy, 2024) pp. 3912–3917
work page 2024
- [74]
-
[75]
Y. Xu, J. Wang, C. Xia, and Z. Wang, Science China Information Sciences66, 222208 (2023)
work page 2023
- [76]
- [77]
-
[78]
Modeling non-Poissonian temporal hypergraphs by Markovian node dynamics
R. Liu, M. Ogura, E. Fonseca dos Reis, and N. Ma- suda, European Journal of Applied Mathematics35, 430 (2024). Supplementary Information on “Modeling non-Poissonian temporal hypergraphs by Markovian node dynamics” Hang-Hyun Jo 1,∗ and Naoki Masuda 2, 3, 4,† 1Department of Physics, The Catholic University of Korea, Bucheon, Republic of Korea 2Gilbert S. Om...
work page 2024
-
[79]
IET DISTRIBUTIONS IN CONTINUOUS TIME Here we analyze the continuous-time variant of the AND and LIN models, deriving IET distributions for hyperedges and nodes. We largely reuse the notations introduced in Section III in the main text, with the caveat thatr ℓh,r hℓ, λh, andλ ℓ are rates, not probabilities, in this section. 1.1. IET distribution for hypere...
-
[80]
(7), (8), (10), (11), and (12), we can compactly writeP e(τ) in Eq
IET DISTRIBUTION FOR HYPEREDGES AS MA TRIX PRODUCTS By combining Eqs. (7), (8), (10), (11), and (12), we can compactly writeP e(τ) in Eq. (6) as Pe(τ) = ⃗Λ⊤ e,f(WeLe)τ−1 We⃗Λe,i,(S12) where ⃗Λe,f ≡ λe(m,0) λe(m,1) ... λe(m, m) ,(S13) Le ≡ ¯λe(m,0) 0· · ·0 0 ¯λe(m,1)· · ·0 ... ... ... 0 0· · · ¯λe(m, m) ,(S14) 3 ⃗Λe,i ≡ 1 Ωe ...
-
[81]
DERIV A TION OF THE IET DISTRIBUTION FOR A NODE IN THE GENERAL CASE The IET distribution for the node,P v(τ), can be written as Pv(τ) = Pr[v(τ) = 1, v(τ−1) = 0, . . . , v(1) = 0|v(0) = 1] = X ⃗ n0,...,⃗ nτ Pr[v(τ) = 1|⃗ µh(τ) =⃗ nτ] "τ−1Y t=1 Pr[⃗ µh(t+ 1) =⃗ nt+1|⃗ µh(t) =⃗ nt] Pr[v(t) = 0|⃗ µh(t) =⃗ nt] # ×Pr[⃗ µh(1) =⃗ n1|⃗ µh(0) =⃗ n0] Pr[⃗ µh(0) =⃗ n...
-
[82]
DERIV A TION OF THE IET DISTRIBUTION FOR A NODE IN A MINIMAL NONTRIVIAL CASE 4.1. AND rule Under the AND rule, we consider a simple nontrivial case in which the focal node belongs to two hyperedges of the same sizem, i.e.,k= 2 andm 1 =m 2 =m. It implies that⃗ µ= (1, m−1, m−1). Usingλ e given by Eq. (1), we write λv,AND in Eq. (22) as λv,AND(⃗ µ, ⃗ n) = 1−...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.