Temporal Hypergraphs | Understanding, Applications, and Research Directions

• Temporal hypergraphs are extensions of temporal networks that represent interactions among arbitrary sets of nodes over time, with varying membership, attributes, activation, duration, and recurrence. Memberships and interactions can be represented through timestamped hyperedges, interval-attributed hyperedges, and time-dependent adjacency tensors.
• Examples include distinguishing group chats in digital communication platforms, traffic patterns in road networks, and identifying shared medical attributes in patient-illness records. Hypergraph systems excel in representing interactions in complex networks efficiently, thus proving highly effective in higher-order tasks like distributing reinforcement learning and hypergraph logic reasoning.
• Applications of temporal hypergraphs have been tried and tested across dynamic scenarios which require decision-making based on group-level interactions, such as traffic control systems and autonomous navigation of vehicles.

Temporal hypergraphs generalize temporal networks by representing interactions among arbitrary sets of nodes whose membership, attributes, activation, duration, and recurrence may change over time. A temporal hypergraph may be represented as a sequence of timestamped hyperedges, $\mathcal E=\{(e_i,t_i)\}$, as a sequence of temporal snapshots $\{(V,\mathcal H(t))\}_{t=1}^{T}$, or as time-dependent adjacency tensors. Unlike pairwise temporal graphs, it preserves the distinction between one joint group event and several independent dyadic interactions. This distinction supports higher-order link prediction, motif analysis, temporal logic reasoning, consensus dynamics, topological persistence, controllability, generative modeling, and domain-specific applications including communication, collaboration, mobility, traffic control, biological interaction, and cybersecurity (Liu et al., 2021, Gallo et al., 2023, Dong et al., 2024).

1. Formal representation and temporal semantics

A static undirected hypergraph is $H=(V,E)$, where $V=\{v_1,\ldots,v_n\}$ is the vertex set and each hyperedge $e\in E$ is a subset of $V$. A hyperedge may contain any number of vertices; ordinary graphs are the special case in which every hyperedge has cardinality two. The incidence matrix $\mathbb I\in\{0,1\}^{|V|\times |E|}$ is defined by

$$i(v_j,e_k)= \begin{cases} 1,&v_j\in e_k,\ 0,&\text{otherwise}. \end{cases}$$

The neighborhood $\mathcal N(v_j)$ contains nodes occurring in at least one common hyperedge with $v_j$. In a temporal setting, the hypergraph may be denoted $\{(V,\mathcal H(t))\}_{t=1}^{T}$0, with both membership and attributes potentially evolving over time (Fischer et al., 2020).

Several temporal formalisms coexist:

• Timestamped hyperedges: a temporal hypergraph is an ordered sequence $\{(V,\mathcal H(t))\}_{t=1}^{T}$1, where $\{(V,\mathcal H(t))\}_{t=1}^{T}$2. Repeated occurrences of the same node set are allowed (Liu et al., 2021).
• Snapshot sequences: the system is represented as $\{(V,\mathcal H(t))\}_{t=1}^{T}$3, where the active hyperedge set changes at discrete times (Gallo et al., 2023).
• Interval-attributed hyperedges: each hyperedge $\{(V,\mathcal H(t))\}_{t=1}^{T}$4 has one or more active intervals $\{(V,\mathcal H(t))\}_{t=1}^{T}$5. This representation supports duration, overlap, and containment relations between events (Myers et al., 2023).
• Directed, predicate-labeled events: an event may be written $\{(V,\mathcal H(t))\}_{t=1}^{T}$6, where $\{(V,\mathcal H(t))\}_{t=1}^{T}$7 and $\{(V,\mathcal H(t))\}_{t=1}^{T}$8 are input and output entity sets and $\{(V,\mathcal H(t))\}_{t=1}^{T}$9 is an interval (Yang et al., 2022).
• Time-dependent tensors: a $H=(V,E)$0-way interaction is represented by an adjacency tensor $H=(V,E)$1, producing a time-varying polynomial dynamical system (Dong et al., 2024).

Temporal semantics therefore need not be restricted to an instantaneous binary event. They may encode exact ordering, inter-event time, persistence, recurrence, duration, or interval relations. The choice affects both inference and interpretation: a timestamped-event model emphasizes event order, whereas interval-based and snapshot-based models represent duration and simultaneous activity.

Temporal hypergraphs also differ in whether time is attached to hyperedges, vertices, or incidences. The topological-analysis framework permits temporal attributes on all three, although its experiments use edge-attributed temporal hypergraphs. Other work focuses on event streams in which a node set is active at a single timestamp, or on snapshots inferred from timestamped bipartite events (Myers et al., 2023, Kirkley, 2023).

2. Higher-order structure and representations

The defining property of a hypergraph is that group membership is represented directly rather than reconstructed from pairwise projections. A question tagged with several topics, a bill with several sponsors, a patient using multiple drugs, a discussion thread involving several users, or a paper with multiple authors can each be represented by one hyperedge (Liu et al., 2021). Clique expansion replaces a group of $H=(V,E)$2 nodes with a complete graph and therefore cannot distinguish a genuine joint event from independent pairwise interactions.

This loss is consequential in temporal settings. A three-node closure may mean that all three nodes appeared in one hyperedge, whereas a triangle may consist of two pairwise interactions through separate hyperedges. Similarly, a group may grow, fragment, or recur as one higher-order object. Pairwise projection preserves co-membership but loses the event-level identity, group size, and joint participation semantics (Liu et al., 2021, Zhu et al., 20 Jun 2025).

Several representations are used to retain or approximate higher-order structure:

Representation	Principal object	Temporal information
Incidence matrix	Node–hyperedge membership	Membership changes or predicted incidence
Timestamped event stream	Ordered pairs $H=(V,E)$3	Event order and recurrence
Interval hypergraph	Hyperedges with active intervals	Duration, overlap, and containment
Adjacency tensor	$H=(V,E)$4-way interaction tensor	Time-dependent interaction strength
Snapshot sequence	Static hypergraph per window	Window-level structural evolution
Star expansion	Nodes plus hyperedge nodes	Typed incidence paths across snapshots

The incidence representation is central to prediction systems. In Hyper-Matrix, rows represent vertices, columns represent hyperedges, and cells represent observed or predicted membership strength. Temporal prediction is formulated as

$H=(V,E)$5

so the target is the appearance, disappearance, or strength of node–hyperedge memberships rather than merely the formation of pairwise links (Fischer et al., 2020).

For higher-order temporal dynamics, tensor representations retain interaction order explicitly. For a $H=(V,E)$6-uniform hypergraph, a symmetric adjacency tensor $H=(V,E)$7 has entries

$H=(V,E)$8

Nonuniform systems use a collection $H=(V,E)$9. Time dependence is represented by $V=\{v_1,\ldots,v_n\}$0, allowing hyperedges to appear, disappear, or change interaction strength (Dong et al., 2024).

Snapshot construction introduces a temporal-resolution problem. Fixed windows can obscure short-lived events or merge unrelated interactions. A minimum-description-length method instead infers the number, duration, and boundaries of snapshots directly from timestamped bipartite events. It minimizes a description length consisting of snapshot widths, event counts, marginal counts, and residual incidence structure, with dynamic programming recovering the optimal binning (Kirkley, 2023). Hyper-core analysis likewise uses a resolution $V=\{v_1,\ldots,v_n\}$1, but independently decomposes each static snapshot and follows the resulting core structures through time (Mancastroppa et al., 2024).

3. Temporal patterns, motifs, and memory

Temporal hypergraphs support structural descriptions that combine overlap patterns with temporal order. TH-motifs define 96 canonical patterns involving three connected temporal hyperedges within a time window $V=\{v_1,\ldots,v_n\}$2. Their signatures are determined by the occupancy of the seven Venn regions generated by three hyperedges, together with the chronological order of their occurrences (Lee et al., 2021).

The 96 motifs comprise:

• Triple-inducing motifs: three distinct underlying static hyperedges;
• Pair-inducing motifs: exactly two temporal hyperedges induce the same static hyperedge;
• Single-inducing motif: all three temporal hyperedges are repeated occurrences of one static hyperedge.

The normalized motif vector

$V=\{v_1,\ldots,v_n\}$3

provides a fixed-length descriptor independent of the number of nodes, hyperedges, or observation duration. Characteristic profiles compare motif frequencies with randomized temporal hypergraphs. Across 11 real-world temporal hypergraphs, motif profiles were more similar within domains than across domains, and temporal motif features improved hyperedge prediction by up to $V=\{v_1,\ldots,v_n\}$4 relative to static h-motif features (Lee et al., 2021).

THyMe+ counts these motifs exactly through a sliding temporal window and a compressed projected graph $V=\{v_1,\ldots,v_n\}$5. Each node of $V=\{v_1,\ldots,v_n\}$6 represents a distinct static hyperedge, while $V=\{v_1,\ldots,v_n\}$7 stores active timestamps of repeated occurrences. Combinatorial procedures count triple-, pair-, and single-inducing motifs without explicitly enumerating every duplicate-induced temporal triple. The reported maximum speedup is $V=\{v_1,\ldots,v_n\}$8 over a static-enumeration baseline and $V=\{v_1,\ldots,v_n\}$9 over direct temporal enumeration (Lee et al., 2021).

Temporal organization also appears as memory across event order and group size. Higher-order correlation functions $e\in E$0 quantify temporal dependence among interactions of the same order, while $e\in E$1 measures directed dependence between different group sizes. The cross-order gap

$e\in E$2

captures asymmetry in temporal association. Empirical face-to-face, office, hospital, and university data show long-range intra-order correlations, order-dependent correlation ranges, banded cross-order dependence, and nonzero directional gaps. Larger groups generally have shorter correlation ranges, while large groups may also fragment or emerge through coordinated changes involving several smaller groups (Gallo et al., 2023).

Relational hyperevent models provide a likelihood-based framework for testing such mechanisms. For a candidate hyperedge $e\in E$3, the conditional probability is

$e\in E$4

Statistics may encode individual activity, repetition of pairs or triples, exact group repetition, triadic closure, subset assortativity, and homophily. RHEM condition on observed hyperedge size; they model participant composition rather than why an event has a particular size. Fitted models can serve as tailored null distributions that preserve selected mechanisms while testing additional higher-order configurations (Lerner et al., 2 Jun 2025).

Temporal locality is also evident in hypergraph ego-networks. Star, radial, and contracted ego-networks exhibit substantial overlap between consecutive simplices, while alter-network interactions are concentrated in time. A classifier trained to distinguish correct from shuffled orders achieved accuracies of $e\in E$5, $e\in E$6, and $e\in E$7 for coauthorship star, radial, and contracted ego-networks, respectively. A hill-climbing procedure reconstructed hidden orders with pairwise ordering accuracies of $e\in E$8, $e\in E$9, and $V$0 for those constructions (Comrie et al., 2021).

4. Dynamics and temporal processes

Temporal ordering can affect dynamics even when time-aggregated structure appears unchanged. For linear average consensus on temporal hypergraphs,

$V$1

the average state remains invariant under the stated connectivity assumptions. Temporal switching changes convergence speed but not the final average consensus. For suitable matched systems, the hierarchy is

$V$2

The result is established for the specified random-switching framework and matrices satisfying $V$3, rather than as a universal theorem for every temporal hypergraph (Neuhäuser et al., 2021).

The distinction changes for nonlinear consensus. A three-way interaction model uses a state-dependent factor $V$4, with the simulations taking $V$5 and $V$6. The resulting weights depend on the evolving state, so temporal order changes the operators applied later in the process. In balanced two-cluster simulations, the group whose local-majority hyperedges become active first exerts a first-mover advantage. This can shift the final consensus even when the numbers of opposing hyperedges are equal. The phenomenon is simulation-supported rather than a general theorem for arbitrary nonlinear temporal hypergraphs (Neuhäuser et al., 2021).

Temporal hypergraphs also support nonlinear contagion models in which simultaneous exposure has a group-level effect. In the EATH evaluation, a susceptible node in hyperedge $V$7 becomes infected with probability

$V$8

where $V$9 is the number of infected members and $\mathbb I\in\{0,1\}^{|V|\times |E|}$0. EATH-generated surrogates reproduced empirical final epidemic size, reproduction number, and temporal incidence more closely than a memoryless ablation. The ablated model produced faster spreading, a higher first peak, and greater total epidemic impact (Mancastroppa et al., 1 Jul 2025).

Consensus and contagion therefore illustrate two different roles of temporal hypergraph structure. In linear systems, higher-order structure may be reducible to weighted pairwise interactions for the dynamics under consideration. In nonlinear systems, state-dependent group effects preserve information that pairwise projections cannot represent.

Topological data analysis offers another dynamical perspective. Each temporal snapshot can be converted into an abstract simplicial complex, with a hyperedge generating a simplex and all of its faces. Because snapshots are not generally nested, zigzag persistence interleaves consecutive complexes with their unions:

$\mathbb I\in\{0,1\}^{|V|\times |E|}$1

Persistence intervals track the birth and disappearance of connected components, loops, and higher-dimensional holes. In Reddit data, a principal $\mathbb I\in\{0,1\}^{|V|\times |E|}$2 component persisted through most of the observation period, while $\mathbb I\in\{0,1\}^{|V|\times |E|}$3 features indicated changes in conversational centralization. In DARPA OpTC data, malicious activity was associated with more transient $\mathbb I\in\{0,1\}^{|V|\times |E|}$4 features and emerging $\mathbb I\in\{0,1\}^{|V|\times |E|}$5 features than benign activity (Myers et al., 2023).

Hyper-core analysis supplies a complementary structural-dynamical description. The $\mathbb I\in\{0,1\}^{|V|\times |E|}$6-hyper-core is the maximal subhypergraph in which each node participates in at least $\mathbb I\in\{0,1\}^{|V|\times |E|}$7 distinct hyperedges of size at least $\mathbb I\in\{0,1\}^{|V|\times |E|}$8. Temporal trajectories of relative hypercoreness

$\mathbb I\in\{0,1\}^{|V|\times |E|}$9

identify movement into and out of higher-order structural cores. Filling-profile similarity, central-core Jaccard similarity, and hypercoreness-rank correlation quantify global, mesoscopic, and microscopic stability. Scientific collaborations exhibited stable macroscopic profiles but substantial turnover in central-core membership, whereas hospital interactions showed stronger role-constrained stability (Mancastroppa et al., 2024).

5. Learning, inference, and prediction

Temporal hypergraph learning includes prediction of memberships, higher-order pattern types, event times, logical consequences, and future pairwise links induced by latent group structure.

Incidence prediction and interactive refinement

Hyper-Matrix combines a geometric deep-learning model with a multilevel matrix visualization. The incidence matrix is factorized as

$$i(v_j,e_k)= \begin{cases} 1,&v_j\in e_k,\ 0,&\text{otherwise}. \end{cases}$$0

and the predictor uses an explicit-characteristic Laplacian $$i(v_j,e_k)= \begin{cases} 1,&v_j\in e_k,\ 0,&\text{otherwise}. \end{cases}$$1 as auxiliary information:

$$i(v_j,e_k)= \begin{cases} 1,&v_j\in e_k,\ 0,&\text{otherwise}. \end{cases}$$2

The system supports interactive relevance feedback. An expert may update a user–topic connection strength between $$i(v_j,e_k)= \begin{cases} 1,&v_j\in e_k,\ 0,&\text{otherwise}. \end{cases}$$3 and $$i(v_j,e_k)= \begin{cases} 1,&v_j\in e_k,\ 0,&\text{otherwise}. \end{cases}$$4, after which the model recomputes effects on neighboring vertices and displays a diverging change scale. The provenance history stores interactions, model outputs, and fixed random seeds. In the law-enforcement case study, the model achieved ROC AUC $$i(v_j,e_k)= \begin{cases} 1,&v_j\in e_k,\ 0,&\text{otherwise}. \end{cases}$$5 and recall $$i(v_j,e_k)= \begin{cases} 1,&v_j\in e_k,\ 0,&\text{otherwise}. \end{cases}$$6 on internet-forum data (Fischer et al., 2020).

Higher-order pattern prediction

HIT predicts the first expansion of an interaction $$i(v_j,e_k)= \begin{cases} 1,&v_j\in e_k,\ 0,&\text{otherwise}. \end{cases}$$7 to a third node $$i(v_j,e_k)= \begin{cases} 1,&v_j\in e_k,\ 0,&\text{otherwise}. \end{cases}$$8. Its four mutually exclusive outcomes are Edge, Wedge, Triangle, and Closure. It uses temporal random walks, asymmetric distance encoding, RNN-based temporal encoding, learned Fourier time features, and task-specific decoders for pattern type, occurrence time, and structural explanation.

The asymmetric encoding is invariant to swapping $$i(v_j,e_k)= \begin{cases} 1,&v_j\in e_k,\ 0,&\text{otherwise}. \end{cases}$$9 and $\mathcal N(v_j)$0, but not to exchanging them with $\mathcal N(v_j)$1. This reflects the role distinction between the initial interacting pair and the newly added node. HIT achieved approximately a $\mathcal N(v_j)$2 average AUC gain over the compared baselines for pattern classification on five real-world temporal hypergraphs and produced uniformly more accurate time estimates (Liu et al., 2021).

Temporal inductive logic reasoning

TILR extends inductive logic programming to directed, predicate-labeled temporal hyperedges

$\mathcal N(v_j)$3

Its rules combine relational variables with Allen interval relations such as BEFORE, DURING, MEETS, OVERLAPS, and EQUAL. Because a many-to-one hyperedge is executable only when all required inputs are available, TILR uses the multi-start random B-walk and B-connectivity rather than ordinary one-node random walks. Path-consistency propagation then constrains temporal relations across sampled paths.

On YouCook2-HG and nuScenes-HG, the full parameterized TILR model achieved MRR $\mathcal N(v_j)$4 and $\mathcal N(v_j)$5, respectively, compared with $\mathcal N(v_j)$6 and $\mathcal N(v_j)$7 for its strongest reported clique-expanded comparisons. The results support retaining n-ary participation, event duration, and interval constraints for procedural and scene reasoning (Yang et al., 2022).

Neural models with inferred or heterogeneous hyperedges

HTGN infers latent temporal hyperedges from pairwise temporal events. In homogeneous graphs it identifies temporally localized maximal cliques; in bipartite graphs, nodes sharing a recent bipartite neighbor form a hyperedge whose membership changes over time. Hyperedge memories are updated with recurrent units and merged when smaller structures are absorbed into larger cliques. The model reports GPU-memory reductions of approximately $\mathcal N(v_j)$8–$\mathcal N(v_j)$9 relative to selected temporal graph baselines (Liu et al., 21 May 2025).

HTHGN targets heterogeneous temporal hypergraphs. It constructs $v_j$0-hop or $v_j$1-ring hyperedges, uses $v_j$2-uniform sampling to control hyperedge size, applies heterogeneous node–hyperedge attention, adds temporal attention across snapshots, and uses contrastive learning to preserve low-order structure. On DBLP, AMiner, and Yelp, the reported full-model AUC values for ordinary link prediction were $v_j$3, $v_j$4, and $v_j$5, respectively; for new-link prediction they were $v_j$6, $v_j$7, and $v_j$8 (Liu et al., 18 Jun 2025).

Traffic-signal control provides an application in which hyperedges are learned dynamically from spatio-temporal observations. HG-DRL constructs spatial hyperedges from current intersection embeddings and temporal hyperedges connecting current intersections to selected nodes from the previous time step. The resulting weighted incidence-like relation is processed by attention within the critic of multi-agent soft actor–critic. On the reported CityFlow scenarios, HG-DRL achieved lower average travel time than GCN-SAC, including reductions of $v_j$9, $\{(V,\mathcal H(t))\}_{t=1}^{T}$00, $\{(V,\mathcal H(t))\}_{t=1}^{T}$01, and $\{(V,\mathcal H(t))\}_{t=1}^{T}$02 on four synthetic settings (Wang et al., 2024).

6. Statistical models, generation, and structural inference

Temporal hypergraph models range from analytically tractable Markov processes to empirically calibrated surrogate generators.

The AR(1) dynamic hypergraph model assigns each possible hyperedge $\{(V,\mathcal H(t))\}_{t=1}^{T}$03 an independent binary process $\{(V,\mathcal H(t))\}_{t=1}^{T}$04, with appearance probability $\{(V,\mathcal H(t))\}_{t=1}^{T}$05 and disappearance probability $\{(V,\mathcal H(t))\}_{t=1}^{T}$06:

$\{(V,\mathcal H(t))\}_{t=1}^{T}$07

Under time homogeneity, the stationary presence probability is

$\{(V,\mathcal H(t))\}_{t=1}^{T}$08

and temporal autocorrelation decays geometrically. Transition-frequency estimators have uniform error bounds of order

$\{(V,\mathcal H(t))\}_{t=1}^{T}$09

under the stated regularity conditions. The model also defines a transition-based Laplacian using appearance probabilities and persistence probabilities separately, enabling spectral community recovery and likelihood-based change-point estimation (Zhu et al., 20 Jun 2025).

EATH is a two-stage surrogate model. Node activity switches between high- and low-activity states, while group formation depends on instantaneous activity, node-specific group-size propensities, long-term pair memory, and short-term group continuation. Groups may remain unchanged or gain or lose one node per time step. The model reproduces node and hyperedge burstiness, event trains, group-size distributions, projected-network structure, and higher-order contagion statistics across eight face-to-face datasets. Its principal limitations are frozen long-term memory, restricted group evolution, and the absence of explicit community structure (Mancastroppa et al., 1 Jul 2025).

A related node-driven model shows how Markovian activity states can generate non-Poissonian event sequences. For a hyperedge of size $\{(V,\mathcal H(t))\}_{t=1}^{T}$10, the event rate depends on the number $\{(V,\mathcal H(t))\}_{t=1}^{T}$11 of high-activity nodes. Under the AND rule, all members must be high to obtain the high rate; under the LIN rule, the rate interpolates linearly between low and high activity. Slow latent-state switching produces mixtures of geometric interevent-time components and slowly decaying autocorrelation, despite Markovian node dynamics (Jo et al., 9 Apr 2026).

The hypergraph stochastic block model in the AR(1) framework makes appearance and disappearance probabilities functions of latent community labels. Its transition-based spectral method exploits both $\{(V,\mathcal H(t))\}_{t=1}^{T}$12, constructed from appearance probabilities, and $\{(V,\mathcal H(t))\}_{t=1}^{T}$13, constructed from persistence probabilities. The primary-school application grouped students approximately by interaction patterns, while the Enron application identified August 2001 as a structural break in the temporal hypergraph (Zhu et al., 20 Jun 2025).

Controllability and observability provide a system-theoretic formulation for time-dependent tensor dynamics. The controlled system is

$\{(V,\mathcal H(t))\}_{t=1}^{T}$14

and the observed system is $\{(V,\mathcal H(t))\}_{t=1}^{T}$15. Recursive tensor-based matrices $\{(V,\mathcal H(t))\}_{t=1}^{T}$16 and $\{(V,\mathcal H(t))\}_{t=1}^{T}$17 generate rank conditions for local weak controllability and observability. The terms $\{(V,\mathcal H(t))\}_{t=1}^{T}$18 and $\{(V,\mathcal H(t))\}_{t=1}^{T}$19 represent directions created or destroyed by temporal variation itself. Greedy driver-node and sensor-node procedures are proposed, but they are heuristic rather than guaranteed globally optimal (Dong et al., 2024).

7. Applications, limitations, and research directions

Temporal hypergraphs have been applied to communication, collaboration, social proximity, online discussion, email, scientific publication, biological and ecological interaction, autonomous driving, traffic control, and cybersecurity. The domain-specific interpretation of a hyperedge varies: it may denote a group event, a shared attribute, a thread, a paper, a common destination, a co-active traffic pattern, or a typed relational event.

Several methodological limitations recur across the literature:

• Temporal resolution: fixed windows, inferred bins, timestamp ties, and coarse time units can alter motif counts, cores, correlations, and learned representations.
• Projection bias: clique or pairwise expansion can create artificial relations and erase joint-event semantics.
• Hyperedge inference: maximal cliques and neighborhood-based constructions are structural proxies and may produce false positives, missing groups, or overly large hyperedges.
• Scalability: hyperedge enumeration, dense incidence structures, large simplicial closures, and temporal attention can become computational bottlenecks.
• Limited event semantics: many models assume undirected, binary, instantaneous, or fixed-node hyperedges, whereas real systems may involve direction, roles, weights, durations, entry and exit, or simultaneous events.
• Interpretability: discriminative structural features, attention weights, and model feedback do not automatically constitute causal explanations.
• Model dependence: RHEM null distributions, AR(1) processes, EATH surrogates, and higher-order correlation models preserve selected mechanisms but leave other mechanisms unconstrained.
• Evaluation scope: several studies rely on formative evaluations, synthetic systems, selected case studies, or datasets constructed through preprocessing choices.

Current research directions include direct learning from observed temporal hyperedge events, uncertainty-aware hyperedge inference, interval-valued and duration-aware models, heterogeneous and directed hypergraphs, adaptive temporal resolution, higher-order motifs involving more than three events, scalable streaming algorithms, temporal hyper-core decompositions, multiscale memory kernels, non-Markovian and cross-order generative processes, and principled integration of topology with dynamical processes.

The principal methodological conclusion is that temporal and higher-order structure should be modeled jointly. Temporal ordering determines which interactions can influence subsequent events; hypergraph structure determines whether participation is pairwise or genuinely collective; nonlinear dynamics can amplify both effects; and statistical, topological, and learning-based methods provide complementary descriptions of the resulting evolution.