Classifying Directional Trajectories Near Criticality in the Three-State Majority-Vote Model with Deep Belief Networks and Bidirectional GRUs
Authors: Mauricio A. Valle, Gonzalo A. Ruz
Organizations: Faculty of Engineering and Sciences Universidad Adolfo Ibáñez Millennium Nucleus for Social Data Science (SODAS) Santiago, Chile · Faculty of Engineering and Sciences Universidad Adolfo Ibáñez Millennium Nucleus for Social Data Science (SODAS) Millennium Nucleus in Data Science for Plant Resilience (PhytoLearning) Santiago, Chile
In this work, we investigate whether the latent representations learned by a Deep Belief Network (DBN) and a Bidirectional Gated Recurrent Unit (Bi-GRU) can discriminate among four dynamically distinct trajectory types in the three-state majority vote model (MV3): approach from disorder, approach from order, departure to disorder, and departure to order. The DBN, pre-trained in an unsupervised manner on static equilibrium samples via a Gaussian-Bernoulli Restricted Boltzmann Machine input layer and architecture 784→4096→225→81, encodes each lattice snapshot into an 81-dimensional latent vector. A t-SNE analysis of the DBN latent space reveals only partial separation of the four trajectory types, reflecting the fact that a model trained on static configurations cannot fully resolve directional temporal structure. A two-layer Bi-GRU classifier, trained on sequences of DBN-encoded snapshots of length T=50, achieves near-perfect separation of all four trajectory types in its hidden state space, as confirmed by t-SNE visualization on both training and test sets. Furthermore, a sliding-window application of the trained Bi-GRU to continuous MV3 dynamics demonstrates its ability to sense the system's current dynamical regime in real-time. These results establish a principled hierarchical architecture for detecting and classifying critical transitions in agent-based opinion dynamics models.
Figures & tables
Fig. 1: Methodology developed to analyze the criticality dynamics of MV3 systems.
Fig. 2: Examples of systems S(p) of L2=282=784 for different values of the parameter p , from more to less noise.
Fig. 3: Example of non-equilibrium phase transition dynamics. The graph shows the magnetization evolution for 100 trajectories in each mode, with T=50 time steps. The thick line represents the average.
Fig. 4: 2D t-SNE projection of the latent representations zt or L(ps→pe) last layer of the DBN. Left: The samples have been colored according to the trajectory mode to which they belong. Right: The samples have been colored according to the system’s magnetization value.
Fig. 5: 2D t-SNE projection of the latent representations htg or G(ps→pe) , the last layer of the 32 hidden units of the Bi-GRU, using T=50 time steps. Left: the resulting projection for the train set. Right: the resulting projection for the test set. These projections were performed using 960 samples from the train and test set respectively, with approximately 240 samples for each of the four trajectories.
Fig. 6: System dynamics and GRU predictions. (a) Evolution of the parameter p (orange dashed line) and the magnetization ∣M∣ (blue line). (b) Probabilities predicted by the Bi-GRU for the four phase transitions. The dotted vertical line indicates the onset of the quench.
Faculty of Engineering and Sciences Universidad Adolfo Ibáñez Millennium Nucleus for Social Data Science (SODAS) Santiago, Chile · Faculty of Engineering and Sciences Universidad Adolfo Ibáñez Millennium Nucleus for Social Data Science (SODAS) Millennium Nucleus in Data Science for Plant Resilience (PhytoLearning) Santiago, Chile
Department of Physics, National University of Singapore Singapore · Institute of High Performance Computing Agency for Science, Technology and Research (A*STAR) Singapore