Paper deep dive
Learning to Trace Seiberg Dualities
Jonathan J. Heckman, Shani Meynet, Alessandro Mininno, Gary Shiu
Intelligence
Status: succeeded | Model: Gemma-4-26B-A4B | Prompt: intel-v1 | Confidence: 92%
Last extracted: 8/3/2026, 10:28:59 AM
Summary
This paper investigates the computational complexity of establishing Seiberg dualities between supersymmetric quiver gauge theories using machine learning. The authors propose using Graph Neural Networks (GNNs), specifically Distance GNNs (DGNN) and Adviser GNNs (AGNN), to guide heuristic pathfinding algorithms (A* and Beam Search) in tracing mutation sequences. Results indicate that for quivers with ~10 nodes, these hybrid AI approaches outperform deterministic algorithms, with a 'Hybrid LCA' pathfinder achieving the best efficiency and accuracy. The work positions this problem as a benchmark for AI in theoretical physics.
Entities (10)
Relation Signals (6)
Distance GNN → estimates → Mutation Distance
confidence 95% · The first NN we trained is a GNN that estimates the distance, i.e., the number of minimal mutations, between a pair of theories.
Seiberg Duality → modelsas → Quiver Mutation
confidence 95% · Mathematically, this involves establishing mutations of quivers
Adviser GNN → guides → Pathfinder Algorithm
confidence 92% · We find it most effective to use the NN as a policy guide for heuristic-driven pathfinders that search over possible ways of dualizing a quiver gauge theory.
A* search → combinedwith → GNN
confidence 90% · This approach of mixing a GNN with a heuristic search algorithm has been proven to be quite efficient
Hybrid LCA Pathfinder → outperforms → Deterministic Algorithms
confidence 90% · This pathfinder has the best overall performance of all the pathfinders we considered if we insist on having a 100% success rate.
Transformer → usedin → Graph Neural Network
confidence 88% · different network architectures consisting of transformers and multi-layer perceptrons tend to outperform deterministic algorithms.
Cypher Suggestions (0)
No Cypher suggestions yet.
Abstract
Abstract:Dualities play an important role in establishing both microscopic and emergent phenomena in a wide range of physical systems. In practice, though, it can often be computationally challenging to establish when two systems are dual, even when all of the "rules of the game" are well-known. Said differently, when confronted with two systems, how can one efficiently establish that they are in fact dual? In this paper we use machine learning methods to address this question for Seiberg dualities of supersymmetric quiver gauge theories. Mathematically, this involves establishing mutations of quivers, which is in turn a variation on the theme of "learning to unknot". On the one hand, this leads us to a practical tool for establishing the computational complexity of different dualities. On the other hand, it also allows us to study how different network architectures learn how to trace Seiberg dualities. We find that for quivers with a modest number of quiver nodes (of order $10$), different network architectures consisting of transformers and multi-layer perceptrons tend to outperform deterministic algorithms. Supplementing the network by well-established pathfinder algorithms (essentially "Google Maps for quivers") leads to an additional improvement in the efficiency and accuracy of the search strategy. We anticipate that this class of questions can serve as a useful benchmark for frontier AI models applied to theoretical physics.
Tags
Links
- Source: https://arxiv.org/abs/2607.28628v1
- Canonical: https://arxiv.org/abs/2607.28628v1
Trouble viewing inline? Open PDF directly →
Full Text
144,679 characters extracted from source content.
Expand or collapse full text
Learning to Trace Seiberg Dualities Jonathan J. Heckman, 1,2 Shani Nadir Meynet, 1 Alessandro Mininno, 3 and Gary Shiu 3 1 Department of Physics and Astronomy, University of Pennsylvania, 209 South 33rd Street, Philadelphia, PA 19104, USA 2 Department of Mathematics, University of Pennsylvania, 209 South 33rd Street, Philadelphia, PA 19104, USA 3 Department of Physics, University of Wisconsin–Madison, 1150 University Avenue, Madison, WI 53706, USA jheckman@sas.upenn.edu, smeynet@sas.upenn.edu, mininno@physics.wisc.edu, shiu@physics.wisc.edu Abstract Dualities play an important role in establishing both microscopic and emergent phenomena in a wide range of physical systems. In practice, though, it can often be computationally challenging to establish when two systems are dual, even when all of the “rules of the game” are well-known. Said differently, when confronted with two systems, how can one efficiently establish that they are in fact dual? In this paper we use machine learning methods to address this question for Seiberg dualities of supersymmetric quiver gauge theories. Mathematically, this involves establishing mutations of quivers, which is in turn a variation on the theme of “learning to unknot”. On the one hand, this leads us to a practical tool for establishing the computational complexity of different dualities. On the other hand, it also allows us to study how different network architectures learn how to trace Seiberg dualities. We find that for quivers with a modest number of quiver nodes (of order 10), different network architectures consisting of transformers and multi-layer perceptrons tend to outperform deterministic algorithms. Supplementing the network by well-established pathfinder algorithms (essentially “Google Maps for quivers”) leads to an additional improvement in the efficiency and accuracy of the search strategy. We anticipate that this class of questions can serve as a useful benchmark for frontier AI models applied to theoretical physics. arXiv:2607.28628v1 [hep-th] 30 Jul 2026 Contents 1 Introduction2 2 Tracing Seiberg Dualities6 2.1 Mutating Quiver Gauge Theories . . . . . . . . . . . . . . . . . . . . . . . .7 2.2 Mutation Trees . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .9 2.3 Seed and Validation Theories . . . . . . . . . . . . . . . . . . . . . . . . . .11 3 Dataset Generation12 3.1 Breadth-First Search over Duality Trees . . . . . . . . . . . . . . . . . . . .12 4 Graph Neural Networks for Seiberg Dualities15 4.1 Distance Graph Neural Network Architecture . . . . . . . . . . . . . . . . .15 4.2 Adviser Graph Neural Network Architecture . . . . . . . . . . . . . . . . . .20 5 Pathfinder for Seiberg Dualities25 5.1 Brute-Force Search Baseline and Branching Factor . . . . . . . . . . . . . . .25 5.2 Pathfinders for Seiberg Dualities . . . . . . . . . . . . . . . . . . . . . . . . .27 5.3 Complexity Comparison . . . . . . . . . . . . . . . . . . . . . . . . . . . . .32 5.4 Pathfinders Evaluation and Efficiency Metrics . . . . . . . . . . . . . . . . .33 6 Results35 6.1 In-Distribution Theories Analysis . . . . . . . . . . . . . . . . . . . . . . . .35 6.2 Out-of-Distribution Theories Analysis . . . . . . . . . . . . . . . . . . . . . .46 7 Failure Modes51 8 Conclusions57 A Seed Theories for Dataset Generation59 A.1 In-Distribution Theories . . . . . . . . . . . . . . . . . . . . . . . . . . . . .59 A.2 Out-of-Distribution Theories . . . . . . . . . . . . . . . . . . . . . . . . . . .62 B Overview of Graph Neural Networks64 B.1 Message Passing Paradigm . . . . . . . . . . . . . . . . . . . . . . . . . . . .64 B.2 Graph Transformers and Global Context . . . . . . . . . . . . . . . . . . . .64 C Overview of Search Algorithms65 C.1 A ∗ Search Algorithms . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .65 C.2 Beam Search Algorithms . . . . . . . . . . . . . . . . . . . . . . . . . . . . .66 1 1 Introduction Dualities provide important complementary perspectives in a wide range of physical systems. Celebrated examples include the Kramers–Wannier duality of the 2D Ising model, the electric- magnetic duality of Maxwell theory, and many modern examples ranging from interacting quantum field theories to string / M- / F-theory to holography. 1 Indeed, dualities provide a window into the structure of many physical systems at strong coupling. More precisely, a “duality” can often connect more than two presentations of a theory. Indeed, in many well-known examples, successive dualities lead to a wild proliferation of different possible descriptions. Turning the question around, if one is confronted with two theories, how does one establish that the pair is connected to one another via a chain of dualities, and if so, how many simple duality operations does it take to connect the pair? In this paper, we study this question for Seiberg dualities [1] applied to quiver gauge theories, i.e., mutations of quivers in the sense of [2, 3]. In the original examples of supersym- metric QCD, this duality relates an “electric” non-abelian gauge theory with a “magnetic” dual non-abelian gauge theory. This can be extended to a wide range of supersymmetric gauge theories with multiple simple gauge group factors, where the basic operation of duality involves acting on one gauge group factor at a time. A celebrated example of precisely this sort is the duality cascade of reference [4], in which there is a pair of gauge groups that share matter in bifundamental representations. This example also prominently features in many string theory applications since it characterizes fractional branes near the warped throat spec- ified by the conifold [4–6]. Physically, the flow of gauge couplings in these examples dictates how duality proceeds (i.e., which gauge group is running to strong coupling). In many cases, however, the structure of dualities is far more involved, and there can often be more than one prescribed answer depending on the parameters of the quiver. For example, [7] showed that the far more generic situation is that there is an infinite family of dual theories. The specific trajectory also exhibits a high level of sensitivity to the initial conditions of gauge couplings, resulting in a huge chaotic range of possible duality sequences. From this perspective, there is already a practical question: Given two supersymmetric quiver gauge theories, are they in fact related via a sequence of Seiberg dualities? If so, what are the most efficient/optimal simple duality operations available to realize these transfor- mations? In physical terms, this helps to establish the level of computational complexity in reaching dual phases of quiver gauge theories. From a holographic standpoint, this amounts to understanding how far into the bulk a given flow proceeds down a throat. 2 1 The literature is vast, so rather than providing a partial and biased list of references, we trust the reader to search the original literature for themselves. 2 It is also of interest in the context of statistical inference in the sense of [8–10]: what is the level of distinguishability between different string compactifications? 2 Treating this question in purely analytic terms is a daunting task, and so we shall instead seek a more practical, artificial intelligence/machine learning (AI/ML)-based approach. 3 As far as we are aware, previous AI/ML approaches to Seiberg dualities have focused on the special case where the quivers are of finite mutation type, i.e., they return to their original form after a finite number of dualities [13]. 4 As already mentioned, however, the far more generic case is the unwieldy scenario in which each successive duality dramatically increases the combinatorial possibilities. To focus on the essential points, in this paper we suppress some of the physical data of a quiver gauge theory. Along these lines, we shall specify a “simplified quiver” by a directed graph, with each node of the graph decorated with a choice of a non-negative integer (i.e., the integer N of an SU(N ) gauge theory). The quivers are dictated by the condition that the gauge theory is free of anomalies; namely, there is a weighted balancing condition on the arrows of the graph. A Seiberg duality on a node corresponds to a transformation of the gauge theory that alters the dualized node and changes the local connectivity of arrows on the graph (see Section 2 for details). A sequence of dualities can thus be viewed as specifying an alphabet, with words formed by concatenating all of these mutations. The underlying mathematical structure thus shares some similarities with the process of determining whether a given knot is related to an unknot [15]. Much as in [15], our aim here will be to compare a few AI methodologies to understand which neural network (N) architectures are most effective at learning how to trace Seiberg dualities. 5 Our purpose here is twofold: On the one hand, we use these different architectures for the practical goal of establishing the existence of dualities and obtaining an estimate of the overall computational complexity of duality transformations. On the other hand, there is also the question of performance, namely how efficient and accurate different approaches are in establishing a path to dualization. Indeed, while we focus our analysis on Seiberg dualities in quiver gauge theories, the same sort of combinatorial methods shows up in a broad class of other computational problems, ranging, e.g., from the mathematics of cluster algebras [2, 3] to scattering amplitudes [17]. As such, we anticipate that the analysis we provide will have far broader applications. Summary of Methodology and Results Let us now turn to a brief summary of our methodology. Additional details are presented later on. As in most machine learning tasks, we first collect/generate a suitable batch of training 3 One way to narrow the search is to also look for simple numerical invariants of the quiver which could be different, e.g., the defect group associated with a quiver, or its matrix Hilbert series. See e.g., [11, 12] for some recent examples of invariants along these lines. 4 See also [14] for the use of AI to study quivers with a finite number of mutations. 5 See also [16] for work on the study of how NNs learn the structure of arithmetic and group laws. 3 data. With this in place, we then train different neural network architectures (detailed below). We then consider the post-training phase, comparing their relative performance against some pre-specified performance metrics. In our case, the quivers used to generate the training data consist of D-branes probing a toric Calabi–Yau (CY) threefold singularity. This leads to a rather rich class of quiver gauge theories that have been extensively studied in the literature (see e.g., [18–20] for reviews). Starting from these “canonical presentations” of a quiver, we proceed to generate our training data by Seiberg dualizing these quivers multiple times. We primarily focus on quivers with 13 or fewer nodes, with the training set built from up to 12 Seiberg dualities. We note that after dualizing, one often reaches a non-toric phase (i.e., after flops of the geometry). In practice, generating the full dataset was accomplished with rather modest computational resources on personal laptops and the UPenn cluster over two to three days, with 256 GB of RAM across 8 CPUs. We then used this data to train different neural network architectures. To keep a broad perspective, we considered a few different choices of graph neural networks (GNNs) with transformer layers [21]. The first N we trained is a GNN that estimates the distance, i.e., the number of minimal mutations, between a pair of theories. We call this network a Distance GNN (DGNN) because, conceptually, it compares the two graphs and attempts to estimate their distance/difference. The accumulated difference is used to estimate how far apart the two quivers are. We show a schematic representation of the DGNN architecture in Figure 6a, while the detailed description of each layer is given in Section 4.1. The second N we trained is another GNN that aims to estimate, in a pair of theories, which vertices are most likely to be dualized to transition from one theory to the other. We call this network an Adviser GNN (AGNN), where the GNN functions as a policy adviser, i.e., it takes the input quiver and provides a probability that will be used by the pathfinder algorithm to make the choice of which node to dualize at each step. We show a schematic representation of the AGNN architecture that we will describe in Section 4.2, in Figure 6b. We find it most effective to use the N as a policy guide for heuristic-driven pathfinders that search over possible ways of dualizing a quiver gauge theory. We considered two kinds of pathfinders: one is a bidirectional A ∗ search pathfinder [22], and the other is a beam search pathfinder. For both, we used various definitions of cost and heuristic functions 6 which guides the search for a path connecting a pair of quivers. The output of the DGNN will be used as the heuristic function, while the AGNN acts as the cost function. A schematic representation 6 Search algorithms seek a path that minimizes specific criteria. These are formalized using a cost function that computes the accumulated penalty of the route taken so far and a heuristic function that estimates the remaining effort to reach the goal. In this sense, a beam search pathfinder is an A ∗ search without a heuristic function and with a reduced search frontier. 4 of the A ∗ search is shown in Figure 13a, with pseudocode in Algorithm 2, while a beam search [23] is shown in Figure 13b and described in Algorithm 3. This approach of mixing a GNN with a heuristic search algorithm has been proven to be quite efficient in the data science literature (see e.g. [24]), but as far as we are aware, the present work constitutes the first application to a problem in theoretical physics. Moreover, we designed a physics-informed pathfinder, which we call the Lowest Common Ancestor (LCA) pathfinder, that is also an A ∗ search pathfinder, with the policy determined by attempting to find the common quiver between a pair of theories with minimal ranks. This LCA pathfinder is inspired by the way (human) physicists approach the problem of finding the path between two theories related by Seiberg dualities, and it is similar to how we generated the database to train our NNs. The last pathfinder we studied is essentially a hybrid approach that unifies all the previous policies into a single pathfinder, guided both by the output of the NNs and the LCA policies. The policies for the pathfinders we designed are summarized in Table 1. We tested the performance of these pathfinders by comparing them first with a Breadth- First Search (BFS) algorithm. In this way, we defined the efficiency ratio as the average number of nodes expected to be explored by a BFS baseline over the number of nodes explored by the heuristic pathfinder. We also kept track of their success rate, namely the percentage of paths found before exhausting the maximum number of steps. Each pathfinder has its pros and cons based on maximizing either the efficiency ratios or their success rates. We analyzed the performance over the dataset used to train the NNs, i.e., the in-distribution (ID) dataset, and over out-of-distribution (OOD) theories. See Table 2a for a summary of these results. We then compared the performance of the pathfinders among themselves, using the LCA pathfinder as the baseline, as summarized in Table 2b. 7 In that table, one can also find the performance of the Hybrid LCA pathfinder, which is the final pathfinder we considered, with both the Hybrid and LCA policies added to guide the A ∗ search. This pathfinder has the best overall performance of all the pathfinders we considered if we insist on having a 100% success rate. The NNs (with our best checkpoints) and the pathfinders are available at our GitHub repository. In order to study the breaking point of the N-guided pathfinders, we also intro- 7 A careful comparison of Tables 2a and 2b reveals an apparent paradox when comparing the ERs of the OOD theories: the average EERs of the LCA relative to the BFS are superior to those of the Hybrid model. However, when evaluating the ER normalized by the LCA baseline, the Hybrid model outperforms the LCA. The reason is that the BFS baseline grows exponentially with the distances, whereas the LCA is efficient in finding the paths. On the other hand, the LCA pathfinder becomes trapped when the distance is not too large, whereas the Hybrid pathfinder is more efficient. We will discuss these performance tests further in Sections 6.1.4 and 6.1.6. 5 duce the concept of “complexity” defined as C = D log 10 K, where D is the distance between two quivers and K is the number of nodes of the quivers. The NNs for which we provide the checkpoints have been trained up to complexity ∼ 11.5, as we considered theories with up to 13 nodes and distances of up to 12 mutations. In the conclusion of the analysis of the breaking point of the N-guided pathfinders, we estimate that the Hybrid or Hybrid LCA pathfinders will perform better than pure LCA up to a complexity ∼ 1.5/2× of the largest complexity seen during the training of the NNs. Structure of the Paper The rest of this paper is organized as follows. In Section 2, we discuss some combinatorial aspects of Seiberg duality and pose the general question: What is the computational complex- ity of establishing whether two quiver gauge theories are Seiberg dual? With this in place, in Section 3, we provide further details on our method for generating training data for our neural networks. Section 4 discusses in greater detail the different GNN architectures, and Section 5 discusses pathfinder algorithms and their hybrid implementation with GNNs. In Section 6, we summarize the results from the various tests. Moreover, in Section 7, we analyze the breaking point of the pathfinders. Finally, we present our conclusions and potential future directions in Section 8. Additional details of GNNs are reviewed in Appendix B; while details on the different pathfinder algorithms used in the main body are summarized in Appendix C. The seed theories used for our training dataset are given in Appendix A. Note Added: As this paper was being finalized, reference [25] appeared, which also con- siders the task of identifying Seiberg dual quiver gauge theories. Our approach is centered on understanding the physics, computational complexity, and inference underlying such duality operations, whereas [25] focuses on evaluating AI models, specifically testing their ability to use a verifier to debug incorrect duality claims. 2 Tracing Seiberg Dualities In this section, we discuss in more detail the central physics goal of the present work: deter- mining how many “simple” Seiberg duality operations are necessary to connect two quiver gauge theories. To frame the discussion to follow, we briefly review the relevant physics and some motivating questions. 6 2.1 Mutating Quiver Gauge Theories We shall be interested in 4D N = 1 quiver gauge theories with SU-type gauge groups. 8 In physical terms, we introduce a collection of 4DN = 1 vector multiplets for the SU(N i ) gauge groups, which we denote by circles. Flavor symmetry factors can also be included, being formally viewed as gauge symmetry nodes at zero coupling. 9 In the quiver gauge theory literature, these are typically presented as squares. In what follows, we focus on quivers without such flavor symmetry factors. The “matter” of the theory is dictated by a choice of adjacency matrix A. While in top- down constructions A arises from a signed intersection pairing matrix, for our purposes it will suffice to consider a presentation in which all entries are non-negative i.e., A ij ∈ Z ≥0 . This entry specifies the number of 4DN = 1 chiral multiplets in the bifundamental representation (N i , N j ) of SU(N i )× SU(N j ). We draw this as a directed arrow that emanates from node i to node j. We also allow for arrows i → j as well as j → i, which in the physics literature are often referred to as “vector-like pairs”. We note that the adjacency matrix is different from its anti-symmetrization B = A− A T . This quiver data amounts to specifying a directed graph along with the decoration provided by rank assignments N i i∈Nodes (the rank of each SU(N i ) being N i − 1). The full specification of the theory is supplemented by additional data such as gauge-invariant superpotential interaction terms, gauge couplings, and other (possibly irrelevant) interaction terms of the QFT. While this additional physical data is quite important in fully specifying a theory, in what follows we shall primarily focus on the defining data of the directed graph (i.e., the adjacency matrix A) and the rank assignments. We denote this data as (A, N) in the obvious notation. This should be viewed as necessary but not sufficient data in fully specifying the 4D QFT. Now, for 4D quiver gauge theories, there are a few basic consistency conditions we need to maintain. First of all, the absence of SU(N i ) 3 gauge anomalies means that the rank as- signments are partially constrained. Consistency requires that for each node i, the following condition holds: 10 No SU(N i ) 3 Anomalies⇒ X j A ij N j − X j N j A ji = 0.(2.1) We do not impose any further constraints aside from the requirement that any duality op- eration produces another theory with rank assignments N i ≥ 0. In special cases such as 8 We will not be dealing with global properties of the gauge group, so we will conflate the usage of gauge group and gauge algebra since it affects nothing in our discussion. 9 One does not need to enforce 4D anomaly cancellation constraints on such nodes. 10 Recall that in our conventions, all entries of A ij are non-negative. Additionally, we are specifically ex- cluding quivers with flavor symmetry factors. These can be included by extending the sum over nodes, but one does not need to impose an anomaly cancellation condition on the non-abelian flavor symmetry itself. 7 superconformal field theories, one often imposes additional constraints on the N i . We do not impose such conditions here. 11 Given a 4D N = 1 quiver gauge theory, we can consider a Seiberg dual presentation of the IR theory. This is an IR equivalence among theories, whose prototypical example is given by Supersymmetric Quantum Chromodynamics (SQCD), i.e., SU(N c ) with N f flavors, which is dual to SU(N f − N c ) with N f flavors and a set of meson fields interacting via a cubic superpotential. This operation naturally extends to quiver gauge theories, where one considers the dualization of a specified node. We refer to this as a simple duality operation. Chaining together multiple such simple duality operations then relates many different quiver gauge theories. We refer to dualization of a quiver at node k by D k (A, N). Generalizing the SQCD example above, the mutation maps N k to N k f − N k , where N k f is the number of flavors at node k, given by the sum of the ranks of all nodes j connecting to k, weighted by the number of incoming arrows A jk . On the other hand, all the ranks of the other nodes in the quiver remain unchanged. In formulae, the dual theory has a rank vector N ′ i = K X j=1 N j A jk − N k if i = k, N i if i̸= k . (2.2) This operation is well-defined so long as the new rank N ′ k remains non-negative. 12 Note also that duality preserves the condition of anomaly cancellation, i.e., we have an equal number of incoming and outgoing arrows (weighted by the ranks of the neighboring nodes). 13 The mutation of the adjacency matrix modifies the matter content. The effects of dualizing node k include the reversal of the chiralities of the fundamental and antifundamental matter fields coupled to k, and the generation of composite meson fields corresponding to paths of length two that pass through k. The mesonic fields add new arrows between adjacent nodes i and j. We schematically show the action of Seiberg duality on node k in Figure 1. In addition to the changes to the directed graph and rank vectors, there is an accompa- 11 It is worth noting that even the condition of anomaly cancellation is a special feature of dealing with a 4D QFT. For example, in closely related 1D quiver quantum mechanics theories, no such condition needs to be imposed. 12 When the rank is too low compared to the number of flavors, a non-perturbative superpotential can be generated (see e.g, [26, 27]). This can often be diagnosed simply by dualizing to a theory with formally negative ranks, clearly a pathological situation. For our purposes, it suffices to work at the level of the directed graph and its ranks assignments. 13 We will also consider theories that are not necessarily anomaly free, but we will not use them as a training set for the neural networks; instead, they will only be used for our out-of-distribution tests. As already mentioned, such quivers also arise in physical systems such as quiver quantum mechanics. We will still impose that the resulting theory after dualization is given by a connected graph. 8 N k N j N i A ik A kj P l N l A lk − N k N j N i A ki A jk A ik A kj Figure 1: Schematic representation of Seiberg duality at the level of quiver theories. The labels on the arrows indicate the multiplicities of the fields in the various bifundamental representations. nying change to the superpotential, i.e., we add a new cubic interaction term that naturally generalizes the SQCD case. As already mentioned, we neglect such superpotential data in what follows, focusing exclusively on the data (A, N). The only way in which the physics of a superpotential enters is that any vector-like pairs generated by Seiberg duality/quiver mu- tation are automatically deleted. This naturally happens because vector-like pairs generated in this way also introduce a mass term that eliminates these contributions anyway. We note that this does not exclude the possibility of vector-like pairs in the original seed theory. Under this assumption, the Fomin–Zelevinsky mutation rule [2, 3] now becomes: A ′ ij = ( A ji if i = k or j = k, max(0, ∆ ij ) otherwise, (2.3) where the net flow of arrows ∆ ij from node i to node j is computed by taking the existing arrows A ij , adding the new mesonic arrows A ik A kj , and subtracting any oppositely oriented arrows: ∆ ij = (A ij + A ik A kj )− (A ji + A jk A ki ).(2.4) 2.2 Mutation Trees Having specified the defining data for mutating a quiver gauge theory, we can now proceed to introduce a sequence of dualities, i.e., we pick a node of the quiver gauge theory Q = (A, N), mutate some node k to reach D k Q, and then mutate again. Starting from a seed theory, we thus obtain a sequence of dualities, reaching a new quiver: Q (n) = D j n ◦·◦ D j 1 Q.(2.5) By inspection, each quiver obtained in this way amounts to a sequence of letters forming a mutation word. There can, in principle, be non-trivial relations between these words (i.e., a braid algebra), and in this sense, the structure of mutations shares many formal similarities with the problem of unknotting [15]. In a physical quiver gauge theory where we also include the data of the gauge couplings and superpotential couplings, the trajectory of an RG flow can often dictate a preferred se- quence of dualities. That being said, there is also a great deal of sensitivity to these boundary 9 conditions, i.e., the RG flow trajectory is chaotic [7]. From this perspective, one can, in prin- ciple, entertain any sequence of mutations. The set of quiver gauge theories that can mutate into each other, therefore, forms a treelike structure. Indeed, at each stage, we can opt to mutate any of the nodes. 14 Taking into account relations in the braid algebra, the number of distinct quivers reached in this way can be smaller. For example, finite-type quivers are precisely those that always return to the same starting point after a finite number of moves. On the other hand, the far more generic cas is that one can mutate endlessly. A simple ex- ample of this sort is the quiver obtained from N D3-branes probing the orbifold C 3 /Z (1,1,1) 3 (see [7]), as schematically depicted in Figure 2. We also show the example of dualizing one of the nodes of C 3 /Z (1,1,1) 3 in Figure 3. Q (0) Q (1) D 1 Q (2) Q (2) D 2 D 3 Q (3) Q (3) D 1 D 3 Q (3) Q (3) D 1 D 2 Q (1) D 2 Q (2) Q (2) D 3 D 1 Q (3) Q (3) D 2 D 1 Q (3) Q (3) D 2 D 3 Q (1) D 3 Q (2) Q (2) D 1 D 2 Q (3) Q (3) D 3 D 2 Q (3) Q (3) D 3 D 1 Figure 2: Seiberg Duality Tree for C 3 /Z 3 inspired by [32, Figure 2]. Each node in the tree is a gauge theory Q (n) with n being the number of non-trivial mutations D j done from Q (0) . Generally, the action of D j leads to different gauge theories. However, C 3 /Z 3 has a permutation symmetry such that each node of the same color corresponds to the same theory up to relabeling of the gauge nodes. 14 In the context of quiver gauge theories realized by D-branes probing a CY singularity, we comment that the braiding relations involve left- and right-mutation on an exceptional collection of sheaves [28–31]. 10 1 11 3 3 3 2 11 33 6 Figure 3: Example of dualization of C 3 /Z (1,1,1) 3 with respect to the red node (top quiver node). In the context of string based constructions, RG flows/duality moves correspond to motion along a warped throat. Moving further into the bulk of a string compactification thus has a direct bearing on how energy scales impact the computational complexity of a given collection of quivers. 2.3 Seed and Validation Theories A priori, there are an infinite number of possible quiver gauge theories one could introduce. To keep things concrete, we primarily focus on theories where the underlying quiver is dual to a well-known string theory construction, e.g., a D3-brane probing a toric Calabi–Yau singu- larity. At a small number of gauge group nodes, the entire list of such quiver gauge theories is completely known, and there are well-known algorithms for extracting more intricate quivers (see, e.g., [20] for an overview). In practice, then, we start with well-known toric Calabi–Yau geometries where the corresponding quiver gauge theory consists of up to 12 gauge group nodes. From this starting point, we then generate our training dataset by mutating many times. We comment that after mutation, the resulting geometry probed by the D3-brane may be in a non-toric phase. To test the efficacy of our method, we compare against a number of different kinds of quiv- ers. First of all, we consider validation on Seiberg-dualized quivers derived from D3-branes probing toric and non-toric Calabi–Yau singularities not included in the original training set. Additionally, we also consider more abstract quivers which are not known to be realized by any D3-brane probing a Calabi–Yau singularity. This case is especially interesting to consider because it directly points to another implicit question: when does a quiver gauge theory arise from a D-brane probing a singular geometry? Having stated the physics motivation and rules, we now build a machine to check for Seiberg dualities. 11 3 Dataset Generation In this section we describe the theories we have used to generate the database for the training of the NNs we considered in our work. We also explain the method used to generate the database via BFS over the duality tree of the theories in Appendix A. 3.1 Breadth-First Search over Duality Trees We now explain how we generated the pairs of theories that we used to train and test the NNs and the pathfinders. We considered the theories described in Appendix A, which we will refer to as seed theories from this moment on. Let us call Q the set of quivers that can be generated by applying Seiberg duality to any of the seed theories Q (0) ∈ Q. Every quiver Q ∈ Q is defined by its rank vector N and adjacency matrix A. We then define D k : Q → Q∪∅ to be the map of a Seiberg duality transformation at node k. The operator maps a quiver to its dual, D k Q (n) = Q (n+1) , according to the rules defined in Section 2. In the following, if mutating a node makes the rank of the dualized node non-positive, or if the quiver breaks into subquivers, we prune that branch; the operator returns ∅. We also introduce the concept of depth d of the duality tree, which is given by the number of non-trivial applications of D k to Q (0) . In order to generate the pairs, we start from Q (0) at depth d = 0. We then make a BFS exploration, layer by layer, withS d denoting the set of quivers Q discovered by mutating Q (0) d times. Moreover, we keep track of all visited quivers to avoid repetition in the database, and we call V d the set of quivers found up to depth d. 15 Starting from S 0 =Q (0) , V 0 =Q (0) , we proceed up to depth d, via the recursive definition: S d = Q (d) = D k Q (d−1) Q (d−1) ∈S d−1 , k ∈1,...,N, Q (d) ̸=∅, Q (d) /∈V d−1 ,(3.1) V d =V d−1 ∪S d .(3.2) A schematic representation of the generation is shown in Figure 4 and the pseudocode is in Algorithm 1. Once we have generated the full duality tree for the quiver Q (0) up to d max , we need to construct the pairs (Q A ,Q B ) of theories divided by their distance d(Q A ,Q B ), i.e., the minimal number of simple Seiberg duality operations that are necessary to go from Q A to 15 Introducing V d allows us to avoid artificially inflating the depth by applying D k twice consecutively to the same node, which returns the original quiver Q. This ensures only “non-trivial” applications of D k are considered. 12 Q (0) Q (1) a Q (1) b ∅ Q (2) a Q (2) b Q (2) c S 0 (Depth d = 0) S 1 (Depth d = 1) S 2 (Depth d = 2) D 1 D 2 D 3 D 2 D 3 D 1 D 3 D 1 → drop Figure 4: Schematic representation of the BFS generated by the mutation D k . Q B . This can be obtained by considering the graph obtained by placing each mutant of the quiver on the nodes, with the edges representing the mutations relating them; we will denote this graph as Q ∗ = (V ∗ ,E ∗ ). We define the duality distance between (Q A ,Q B ) as the length of the shortest sequence of edges connecting them in this graph: d(Q A ,Q B ) = minp : p number of edges between Q A and Q B in Q ∗ .(3.3) This distance is illustrated in Figure 5. The distance is computed during the BFS generation. Every discovered state Q i ∈V ∗ is associated with a specific mutation path sequence P i = (k i,1 ,k i,2 ,...,k i,d i ) originating from Q (0) . Since all the quivers Q∈Q are generated up to depth d max , the distance between any two quivers Q A and Q B is obtained from their path vectors. We define the LCA of the two theories as the state corresponding to the longest shared prefix between their path sequences: P LCA = P A ∩ P B .(3.4) The shortest path connecting Q A to Q B is constructed by taking the union of their path se- quences and removing the shared intersection (the path from the root to the LCA). Expressed as a formula, this is given by: d(Q A ,Q B ) =|P A | +|P B |− 2|P LCA |.(3.5) In generating the dataset, we keep track of the distance as well as the specific path we used to compute it. However, as is clear from Figure 5, this path is, in general, not unique. Moreover, if the quiver has permutation symmetries, the distance (3.5) becomes an upper 13 bound. As we will explain in Section 4.1.1, we trained the NNs using (3.5) as the target distance to be learned, without identifying quivers related by symmetries. This allows us to validate the training using the mean absolute error (MAE) between the predicted distance and (3.5). However, as we will see later, the DGNN loses accuracy at large distances, predicting shorter distances compared to the database distance (3.5). Once the permutations are taken into account, the accuracy of the network improves, meaning that the model was already accounting for these symmetries. Q A Q B k i k j k m Figure 5: Schematic representation of the duality distance d(Q A ,Q B ) = 3, defined as the shortest path (red) corresponding to the dualizations D k m ◦ D k j ◦ D k i Q A . Algorithm 1 BFS over Duality Trees 1: Input: Initial quiver Q (0) , maximum depth d max 2: Output: Set of visited quivers V 3: S 0 ←Q (0) 4: V ←Q (0) 5: for d← 1 to d max do 6: S d ←∅ 7: for Q (d−1) ∈S d−1 do 8:for k ← 1 to K do▷ Iterate over all nodes 9:Q (d) ← D k Q (d−1) 10:if Q (d) ̸=∅ and Q (d) /∈V then▷ Check validity and uniqueness 11:S d ←S d ∪Q (d) 12:V ←V ∪Q (d) 13:end if 14:end for 15: end for 16: if S d =∅ then 17:break▷ Terminate if no new states are generated 18: end if 19: end for 20: return V 14 4 Graph Neural Networks for Seiberg Dualities In this section, we describe the specific Graph Neural Networks (GNNs) we designed for this work. See Appendix B for a review of terminology as well as a general introduction to the Graph Neural Network paradigm. We primarily consider two types of GNN architec- tures, though they will serve two different purposes. The first one, called the Distance GNN (DGNN), is trained to estimate the distance d(Q A ,Q B ) between two theories Q A and Q B related by Seiberg duality. The second GNN, called the Adviser GNN (AGNN), predicts the first node to mutate in the sequence of Seiberg dualities that brings Q A to Q B . Both our architectures share the building blocks of standard MPNN [33] and Graph Transformers (similar to GraphGPS models, where GPS stands for General, Powerful, and Scalable [34]), but we adapt them to evaluate duality operations rather than classification tasks. Specifically, they use a message-passing tokenizer followed by a Transformer encoder. The tokenizer captures the structure of each node in the quiver, e.g., the rank of its gauge group or the number of edges connecting to it, while the global attention mechanism allows the network to process the totality of the graph. The differences lie in how the architectures handle the outputs. Graph Q A x A,i , w A,ij GNN Tokenizer (3 Layers) Transformer (2 Layers) Embeddings h A,i ∈ R H Global Mean z A Graph Q B x B,i , w B,ij GNN Tokenizer (3 Layers) Transformer (2 Layers) Embeddings h B,i ∈ R H Global Mean z B Node-wise |h A,i − h B,i | ↓ Mean (d local ) Concat [z A ⊕ z B ⊕ d local ] MLP Head (ReLU → Drop) Softplus Output ˆ d Shared WeightsShared Weights (a) DGNN Architecture. Graph Q A y A,i , LPE A Hybrid GPS Encoder Dense Nodes h A,i Graph Q B y B,i , LPE B Hybrid GPS Encoder Dense Nodes h B,i Node-wise Diff ∆ deep i = h B,i − h A,i ∆ raw i and ∆A i Concat Z i [h A ⊕ h B ⊕ ∆] MLP Classifier (Node-Level) Action Masking M Softmax P (D i |Q A , Q B ) Shared Weights (b) AGNN Architecture. Figure 6: Schematic representations of DGNN and AGNN architectures. 4.1 Distance Graph Neural Network Architecture We now explain how the DGNN architecture is constructed. Its purpose is to estimate the distance d(Q A ,Q B ) between two quivers Q A and Q B related by mutations. The DGNN ar- 15 chitecture is depicted in Figure 6a, showing how it processes a pair (Q A ,Q B ) through weight- sharing encoding branches to map them into a latent space. Processing the latent space leads to the prediction of a continuous distance ˆ d∈ R + . First of all, every quiver Q is represented by its ranks and adjacency matrix (N, A). These inputs are element-wise normalized as: x i = log(1 + N i ), w ij = log(1 + A ij ).(4.1) The normalized inputs are then given to an encoder consisting of a 3-layer GNN tokenizer. The primary role of this module is to construct locally-aware “tokens”, so that starting from the first layer encoding the information about the ranks of the nodes x i , the N can process more and more information about the edges connecting each node. Let h (ℓ) i ∈ R H represent the hidden feature vector of node i at layer ℓ, 16 where the initial input is the 1-dimensional normalized rank h (0) i ≡ x i ∈ R 1 . The convolution layers are parametrized by two distinct learnable weight matrices each: W (ℓ) self for self-loops and W (ℓ) msg for neighbor messages. 17 At each step, the unactivated update rule aggregates information from the immediate neighborhood: m (ℓ+1) i = W (ℓ) self h (ℓ) i + X j∈N (i) W (ℓ) msg h (ℓ) j · w ji .(4.2) The weights W (ℓ) self and W (ℓ) msg are shared between Q A and Q B . These vectors m (ℓ+1) i are normal- ized via LayerNorm and passed through a LeakyReLU activation function with a negative slope α = 0.2. 18 The fully activated node representation is therefore: ̃ h (ℓ+1) i = LeakyReLU LayerNorm m (ℓ+1) i .(4.3) Crucially, residual skip connections are enforced for all layers ℓ > 0 such that: h (ℓ+1) i = h (ℓ) i + ̃ h (ℓ+1) i .(4.4) Standard convolutional GNNs replace the node state at each step. This architecture differs by using residual connections to accumulate the outputs, adding the new features ̃ h (ℓ+1) i directly to the prior state h (ℓ) i . We accumulate these features to ensure that the network retains some of the initial information of the quiver extracted in previous layers. Consequently, the final token representation h (3) i is a sum of the information gathered by all the convolutional layers, and it will be the token passed through the transformer we describe next. 16 In our case, H = 64. 17 For the first layer ℓ = 0, the weight matrices project from the 1-dimensional input to the hidden dimension H, i.e., W (0) self , W (0) msg ∈ R H×1 . For subsequent layers ℓ > 0, the matrices are square, W (ℓ) self , W (ℓ) msg ∈ R H×H . 18 The precise definitions of LayerNorm and LeakyReLU are discussed in Appendix B.1. 16 While the GNN tokenizer has captured the local information of each node and its sur- roundings up to ℓ = 3, when comparing (Q A ,Q B ), we need a way to understand the global changes in the quiver that distinguish them after dualization. To achieve this, we pad h (3) i into X ∈ R K×H and process it using a 2-layer Transformer encoder via the Scaled Dot-Product Self-Attention mechanism [35]. We use Multi-Head Attention (MHA) with h = 4 independent heads and a projection dimension d k = H/h = 16, followed by a position-wise Feed-Forward Network (FFN). The progression of the node state X through a single encoder layer uses standard residual connections and LayerNorm: ̃ X = LayerNorm(X + MHA(X)), X out = LayerNorm( ̃ X + FFN( ̃ X)). (4.5) After passing through the 2 distinct Transformer layers, the model yields h A,i and h B,i for the paired quivers Q A and Q B , respectively. From these vectors, we define a node-wise difference vector d local ∈ R H : d local = 1 K K X i=1 |h A,i − h B,i | ,(4.6) which explicitly records the local discrepancies caused by the sequence of Seiberg dualities. Finally, we define the average of h A,i and h B,i over the nodes of the quiver: z A = 1 K K X i=1 h A,i , z B = 1 K K X i=1 h B,i ,(4.7) which provide global representations of each theory to the distance regressor. The averag- ing operation smooths out the individual node features, giving the macroscopic difference d global = |z A − z B |. For our purposes, we need the combined information of d global and d local . Indeed, unlike standard GNNs that typically pool node features to process the graph data (e.g., for classification purposes), the DGNN computes the absolute difference between corre- sponding nodes of two graphs, and it needs to retain finer information. If this discrepancy were evaluated solely after global pooling (yielding only d global ), localized node variations would be smoothed out. Incorporating the pre-pooling difference d local ensures that the regressor retains fine-grained structures that global averaging alone would discard. After the various encoding layers, the continuous distance is processed through a Multi- Layer Perceptron (MLP) head, with [z A ⊕z B ⊕d local ]∈ R 3H as input. The purpose of this part of the N is to reduce the 3H-dimensional vector down to a single scalar. To achieve this, we use standard linear projection with a ReLU activation function and a Dropout regularization layer (p = 0.2) to mitigate overfitting. The output of this projection is passed through a Softplus function, i.e. ˆ d = Softplus(s) = log(1 + e s ),(4.8) leading to the predicted distance ˆ d. 17 4.1.1 Training Evaluation 050100150200250 2 4 6 8 Epoch Loss Train Loss Val Loss (a) Training and Validation MSE Loss 050100150200250 1 1.5 2 Epoch MAE Val MAE (b) Validation MAE Figure 7: Training logs for the DGNN model The network parameters are optimized by minimizing the Mean Squared Error (MSE) between the prediction ˆ d and the path distance d computed at the generation of the database. For a given batch of size B, the loss function is defined as: L MSE = 1 B B X k=1 ( ˆ d k − d k ) 2 .(4.9) While MSE is utilized for gradient descent to penalize outlier predictions, the accuracy of the model is tracked using the Mean Absolute Error (MAE): MAE = 1 B B X k=1 | ˆ d k − d k |.(4.10) In the context of distance regression, MAE provides a direct measure of the average step deviation. The DGNN was trained for approximately 250 epochs to estimate the distance between quivers. Figure 7 displays the training trajectory. Checkpoints for the final evaluation were selected based on the epoch minimizing the MAE. Specifically, the checkpoint used for the pathfinder was reached at epoch 247 with a final MSE loss of 1.137. There can be multiple reasons why the MAE could not go lower. One possibility is that whenever the distance between two quivers is larger than the number of nodes, the graph difference taken by the DGNN of the hidden layers of the N becomes less reliable, leading to an inability to estimate larger distances than the number of nodes of the quivers. This is 18 corroborated by looking at the MAE by the number of nodes in Figure 8, where it is clear that the N struggled to predict larger distances for smaller quivers while being extremely accurate for larger quivers. 123456789101112 0 2 4 6 8 10 Distance MAE (a) MAE by Distance 345678910111213 0 2 4 6 8 6 . 547 3 . 785 2 . 680 2 . 232 1 . 551 1 . 476 1 . 132 1 . 280 1 . 155 1 . 229 1 . 224 Nodes (b) MAE by Nodes −12−10−8−6−4−2024 0 5k 10k 15k 20k 25k Error Pairs Median: −0.954 (c) Overall Error Distribution ( ˆ d− d) Figure 8: Benchmark results for DGNN. Because we apply the DGNN as the heuristic function of an A ∗ search, as we will ex- plain in Section 5.2, we have decided to perform further benchmark tests on the kind of heuristic function that the DGNN can provide. One test is to check if the heuristic function is monotonic, so that we reduce the possibility that the A ∗ search requires processing the same node multiple times. However, as we show in Figure 9, the DGNN output exhibits non-monotonicity. More precisely, we note that the DGNN is non-monotonic approximately 32.07% of the time. Consequently, during a pure A ∗ search, the priority queue can be delayed by repeatedly exploring uninformative branches. We will see that, despite the high success of the A ∗ search guided by the DGNN for both in distribution (ID) and out of distribution (OOD) theories, it is precisely the absence of monotonicity in the DGNN that causes the algorithm’s failure. 19 24681012 10 20 30 40 50 60 Distance Monotonicity Violation (a) Monotonicity Violations by Distance 0123456 0 50k 100k 150k 200k 250k 300k ∆h Pairs Violation: 32.07% (b) Jump Distribution Figure 9: Evaluation of DGNN monotonicity. The non-monotonicity of the heuristic results in excessive node expansions and timeouts. 4.2 Adviser Graph Neural Network Architecture We now proceed to describe the AGNN architecture. The objective of the AGNN is to predict the most probable node that should be mutated first in the sequence of Seiberg dualities from quiver Q A to quiver Q B . Unlike the DGNN, which predicts a distance, the AGNN acts as a localized policy function, yielding a probability distribution over the K nodes of Q A . The AGNN architecture is described in Figure 6b. First of all, much as in the DGNN case, we want to encode the information of the rank and adjacency matrix of the quiver via a suitable “tokenization”. The initial feature vector is 5-dimensional: y i = log(1 + N i ), log(1 + A in i ), log(1 + A out i ), log(1 + N in i ), log(1 + N out i ) (4.11) where A in i = K X j=1 A ji , A out i = K X j=1 A ij , N in i = X j A ji N j , N out i = X j A ij N j .(4.12) To allow the network to distinguish nodes, we evaluate the Laplacian Positional Encodings (LPE). 19 Let us denote as LPE i ∈ R d LPE the Laplacian eigenvectors. These LPE i vectors and the features x i are projected into the hidden layer of dimension H, summed, and normalized 19 Let A be the adjacency matrix of the directed quiver. We define the symmetric adjacency matrix of the underlying undirected graph as A sym = A + A T . The degree d i of node i is, then, the total number of edges connected to the node, i.e. d i = P j (A sym ) ij . The corresponding graph Laplacian is L =D− A sym , where D is the diagonal degree matrix with entries D i = d i and D i̸=j = 0. 20 to form the initial embeddings h (0) A,i and h (0) B,i : h (0) i = LayerNorm (W node y i + W LPE LPE i ) ,(4.13) where W node ∈ R H×5 and W LPE ∈ R H×d LPE are learnable projection matrices, shared by both quivers Q A and Q B . These initial embeddings are then passed through a Hybrid GPS Encoder, identical in structure to the combined GNN Tokenizer and Transformer used in the DGNN, from which we extract h A,i and h B,i , as before. However, despite what one might expect, we find that these more refined inputs empirically under-perform when used in the training of the DGNN. Indeed, the introduction of the LPE i vectors as an input degraded training performance, most likely because it made it harder to extract macroscopic properties of the graph, resulting in higher error rates when predicting distances. Similarly to the DGNN case, we also compute node-wise differences: ∆ deep i = h B,i − h A,i ,∆ raw i = h (0) B,i − h (0) A,i ,(4.14) as well as the element-wise absolute difference of the adjacency matrices: ∆A i = X j |A B,ij − A A,ij |.(4.15) In this way, we keep track of all possible local differences between the two quivers. The dense nodes and these local differences are concatenated into Z i ∈ R 4H+6 for each node: Z i = h h A,i ⊕ h B,i ⊕ ∆ deep i ⊕ ∆ raw i ⊕ y A,i ⊕ ∆A i i .(4.16) The concatenated vector Z i is then processed by an MLP classification head, which projects the features down to a single scalar value for each node i. This scalar, called logit, represents the network’s preliminary preference for mutating node i. However, not all nodes can, in general, be dualized since, for instance, mutating certain nodes might yield negative gauge ranks or disconnect the quivers. We generate a binary action mask M prior to evaluation, which takes the value M i = 1 if the mutation at node i is permissible, andM i = 0 otherwise. In practice, the network checks if N ′ i = P j A ji N j −N i is positive. Before converting the logits into probabilities, we apply this mask by manually overriding the scores of all invalid nodes to −∞. Finally, we apply a Softmax function across the N nodes. Because the exponential of−∞ is exactly zero, this masking guarantees that the network assigns a 0% probability to any mutation that is not possible, forcing the remaining valid nodes to absorb the entirety of the probability distribution. The final policy is thus given by: P (D i |Q A ,Q B ) = Softmax( ̃y) i ,where ̃y i = ( y i if M i = 1, −∞ if M i = 0, (4.17) 21 where y i denotes the logit output of the MLP. Once this valid probability distribution is obtained, we sort the probabilities P (D i |Q A ,Q B ) in descending order and extract the top-k nodes. 4.2.1 Training Evaluation 01020304050 0.8 1 1.2 Epoch Loss Train Loss Val Loss (a) Training and Validation Cross-Entropy Loss 01020304050 40 60 80 Epoch Accuracy (%) Top-1Top-2Top-3 (b) Next-Step Action Accuracy Figure 10: Training logs for the Autoregressive model. The AGNN is optimized by minimizing the Cross Entropy (CE) loss between the predicted probability distribution and the exact first-step mutation extracted from the pre-computed shortest path sequence. For a given mini-batch of size B, let t q ∈1,...,K denote the index of the node that must be dualized in the q-th quiver Q A,q . The loss function penalizes the network based on the predicted log-probability assigned to these actions and is defined as: L CE =− 1 B B X q=1 logP (D t q |Q A,q ,Q B,q ).(4.18) During training, we check the prediction accuracy of the top-k nodes, which monitor whether the node we predict to be the one to dualize, i.e., t q , appears within the highest probability choices predicted by the network. However, it is crucial to note that this metric does not perfectly reflect the true success rate of the AGNN. Due to quiver symmetries and the non- uniqueness of shortest paths in the duality tree, multiple nodes may represent equally optimal first-step mutations. The Cross Entropy loss penalizes the network for deviating from t q in our dataset, even if the AGNN can, in principle, propose a perfectly valid alternative node. Consequently, the top-k accuracy serves as a lower bound on the true success rate of the policy. 22 The AGNN was trained to predict the next-step Seiberg duality operation, with the N outputting the probabilities over the nodes. Figure 10 shows the cross-entropy loss and top-1 next-step accuracy across the 49 training epochs, at which point the training was stopped to prevent overfitting the data, since the training loss kept decreasing while the validation loss started increasing. However, the accuracy converges rapidly, with the optimal validation loss achieved at epoch 39. Due to the symmetry of the quivers, the top-1 accuracy does not exceed 70%, meaning that for a better policy, we need to rely on more than one probability. In fact, top-3 accuracy is, on average, around 96%, resulting in the correct node being dualized most of the time. This can be seen from Figure 11. We note that the accuracy actually increases with distance, resulting in a better oracle whenever two quivers are far apart, while the accuracy drops for smaller distances. However, unlike the DGNN, there is not much difference in accuracy for the different theories used for training. The AGNN will also be utilized to guide a pathfinding search. Initially, we evaluate it as a unidirectional pathfinder, where the trajectory is constructed by iteratively selecting the actions with the top-k probabilities predicted by the AGNN. To benchmark the reliability of the AGNN in this guiding role, we analyze the policy margin, defined as the difference between the probability assigned to the correct node (P correct ) and the highest probability assigned to an incorrect node (P best incorrect ). 20 A negative policy margin indicates an inversion, meaning the network prefers an incorrect operation over the correct one. As shown in Figure 12, the model exhibits such inversions 20% to 40% of the time. However, despite this inversion rate, the average policy margin remains positive across all target distances; this occurs because the network is highly confident when it is correct (yielding large positive margins), easily outweighing the margins of the incorrect predictions, which tend to be only slightly negative. From the discussion above, we can conclude that, considering the three most probable nodes as branching, the AGNN is correct over 90% of the time. This insight is critical: while the AGNN is not an infallible oracle and cannot navigate the duality tree on its own, it suppresses the vast majority of incorrect branches. This confidence motivates the last kind of pathfinder that will use the AGNN policy to restrict the branching factors, allowing the DGNN to evaluate only the most promising paths. Furthermore, analyzing the tail of the margin distribution in Figure 12 reveals that for a substantial subset of states (approximately 5.7× 10 4 theories), the network exhibits an extreme policy margin of P correct −P best incorrect ≈ 1.0. This signifies that the AGNN model is highly polarized: it appears to be uncertain between two options for only a small subset of nodes. Large regions of the duality tree are instead considered trivial and are solved with high 20 Because multiple valid paths can connect two theories, we define an ‘incorrect’ node as any node that does not belong to the specific shortest path recorded during our dataset generation. 23 0.110.660.580.820.840.870.910.860.971.001.001.00 0.230.770.810.790.760.780.750.800.820.880.910.94 0.400.860.820.740.690.680.680.690.700.770.800.83 0.590.880.790.710.640.620.590.630.620.680.670.78 0.760.880.830.770.700.660.600.630.650.680.690.68 0.830.830.780.710.630.620.580.600.570.620.600.67 0.530.790.710.740.730.670.650.660.640.630.630.61 0.390.680.610.670.620.610.570.560.540.570.550.61 0.320.660.600.700.610.580.590.570.560.590.540.61 0.350.570.530.610.580.560.540.520.510.530.510.54 0.360.570.560.650.620.580.600.540.550.540.550.49 123456789101112 3 4 5 6 7 8 9 10 11 12 13 Distance Nodes 0 0.2 0.4 0.6 0.8 1 (a) Top-1 Accuracy by Nodes and Distance 24681012 0 0.2 0.4 0.6 0.8 1 Distance Accuracy Top-1Top-2Top-3 (b) Top-k Accuracy Decay over Distance 24681012 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 Distance Cross-Entropy Loss (c) Mean Cross-Entropy Loss by Distance Figure 11: Benchmark results for the AGNN. While next-step accuracy degrades at larger duality distances, the validity rate remains high, indicating that the network has implicitly learned the rules of Seiberg dualities. 24 24681012 −1 −0.75 −0.5 −0.25 0 0.25 0.5 0.75 1 Distance P correct − P best incorrect Top-1Top-2Top-3 (a) Policy Margin Decay by Distance −1−0.500.51 0 5k 10k 15k 20k 25k 30k P correct − P best incorrect Pairs Top-1Top-2Top-3Inversion (b) Top-k Policy Margin Distributions Figure 12: Policy Margin analysis. While Top-1 inversions cause some paths to fail, the high Top-3 margins enable pruning of a large portion of the duality tree. confidence. From an algorithmic perspective, this is beneficial, as a guided search experiences a zero effective branching factor in these states. However, this same polarization means that when the network does make an error (inversions where the margin is < 0), it also often does so with high confidence. This extreme confidence in both correct and incorrect directions further highlights why a pure greedy search with AGNN is brittle, necessitating the safety net of the DGNN to backtrack out of confident dead-ends. 5 Pathfinder for Seiberg Dualities In this section, we discuss the other component of our inference model, namely the guided pathfinders that are used to determine the chain of dualities that connect two theories Q A and Q B . The pathfinders are of two kinds: a Bidirectional A ∗ Search Pathfinder and a Beam Search Pathfinder, as generally reviewed in Appendix C. For the latter kind, the beam is guided by the output of the AGNN model. On the other hand, we propose different kinds of cost and heuristic functions for the A ∗ Search Pathfinder, as summarized in Table 1. 5.1 Brute-Force Search Baseline and Branching Factor To benchmark the efficiency of guided search, we first formalize the complexity of unguided exploration. In standard graph search algorithms, the branching factor b characterizes the number of direct children a node possesses. In the present framework, the branching factor of a specific quiver Q = (N, A) corresponds to the number of valid Seiberg duality transformations it admits. Let K be the total number of nodes in the quiver. A priori, one might attempt 25 to dualize any of the K nodes, implying a maximum branching factor of K. However, as established in Section 3, a mutation D k Q is admissible if and only if it satisfies two conditions: 1. Rank Positivity : The mutated gauge rank must remain positive. 2. Graph Connectivity : The resulting adjacency matrix A ′ still describes a connected quiver. Therefore, the valid branching factor b Q for any state Q is defined as the number of valid mutations: b Q =|k ∈1,...,K| D k Q̸=∅|≤ K .(5.1) For a bidirectional search between Q A and Q B , b A and b B designate the valid branching factors: • b A ≡ b Q A : The number of valid dualizations originating from Q A . • b B ≡ b Q B : The number of valid dualizations originating from Q B . To establish a baseline for the expected number of mutations performed by an unguided bidirectional search (N BFS ), we observe the following fact: since Seiberg duality is an involu- tion, i.e. D k ◦ D k Q = Q, 21 the branching factor can be estimated to be, at depth i ≥ 1, at most b avg. (b avg. − 1) i−1 . This means that starting from Q, we reach distance d with N dir (d) dualities, which can be expressed as 22 N dir (d) = 1 + b Q d X i=1 (b Q − 1) i−1 = 1 + b Q (b Q − 1) d − 1 (b Q − 1)− 1 .(5.2) If we start a bidirectional search between two quivers at distance d, to guarantee finding an intersection, the unguided pathfinder that explores must reach depth d/2 from each side. To be precise, one of the two sides will reach⌊d/2⌋ depth, and the other⌈d/2⌉. The expectation is the sum of these two expansions: N BFS ≈ N dir (⌊d/2⌋) + N dir (⌈d/2⌉) = 2 + b avg. b avg. − 2 (b avg. − 1) ⌊d/2⌋ + (b avg. − 1) ⌈d/2⌉ − 2 . (5.3) This formulation implies asymmetric growth. When d changes from an even to an odd dis- tance, one tree must expand by an entire extra depth level, causing a spike in the total number of explored nodes. Conversely, when d goes from odd to even, the second tree catches up to the same depth, resulting in a smaller multiplicative increase. The exponential growth of N BFS therefore oscillates depending on the parity of d, an artifact that will be visible in the efficiency metrics defined in Section 6. 21 Strictly speaking, this requires a left- and right-mutation on the exceptional collection of sheaves rather than two mutations of the same sort (see e.g., [36]). Such distinctions will not matter in the present analysis but would be interesting to study further. 22 If b Q = 2, N dir (d) = 1 + 2d. 26 5.2 Pathfinders for Seiberg Dualities We now proceed to present the N-informed A ∗ pathfinder we developed in detail. A sum- mary of the algorithms is given in Appendix C, while their evaluations, costs, and objective functions are given in Table 1. Since quivers are defined up to node permutations, at each step of the pathfinder, we compare a canonicalized version of the quivers via a Weisfeiler–Lehman (WL) hashing of the graphs. 23 Once a pair is identified, the search terminates, and the duality sequence is reconstructed. Q A Q B Q n ′ g fwd g bwd h fwd (Q n ) h bwd (Q m ) V fwd V bwd (a) Bidirectional A ∗ search pathfinder. Q A Top 1 Top 2 Pruned Q B (b) Beam search pathfinder Figure 13: Schematic representations of the two kinds of pathfinders we consider in this work. In contrast to the unguided expansion N BFS , the total computational cost of the guided pathfinder is the sum of the unique quivers explored: N exp =|V fwd | +|V bwd |.(5.4) 23 This avoids missing the pair because of a “wrong” ordering of the nodes. 27 If the pathfinder is unidirectional, as will be the case for the AGNN pathfinder, thenV bwd =∅. If the search exceeds a predefined threshold of steps without finding an intersection, execution halts. This outcome indicates either that the two theories are not related by Seiberg dualities, or that their separation exceeds the maximum search depth. Pathfinder SearchCost g(.)Heuristic h(.) DGNNBidir. A ∗ 1DGNN(Q n ,Q target ) AGNNB = 3 Beam− logP (D i |Q n ,Q target )0 HybridBidir. A ∗ 1− λ AR logP (D i |Q n ,Q target )λ DGNN DGNN(Q n ,Q target ) LCABidir. A ∗ c(D i )0 Hybrid LCA Bidir. A ∗ (1− λ detcost ) + λ detcost c det (D k ) −λ AR logP (D k |Q n ,Q target ) λ DGNN DGNN(Q n ,Q target ) +λ LCA P rank(Q n ) P rank(Q root ) Table 1: Summary of Pathfinders and their cost and heuristic functions. BFS Pathfinder: As a fundamental baseline for evaluating our neural network-guided models, we implement an unguided Breadth-First Search (BFS) pathfinder. This algorithm operates by systematically exploring the Seiberg duality graph level by level, generating all admissible mutations at depth t before evaluating any transition at depth t+1. Structurally, it represents an unguided A ∗ search where the heuristic function is identically zero (h(Q)≡ 0) and the cumulative path cost simply tracks the traversal depth (g(Q n+1 ) = g(Q n ) + 1), assigning an equal step cost to every edge. Consequently, the BFS pathfinder is topologically agnostic, treating every valid mutation as equally promising without any predictive physical intuition. To optimize execution and match the structural framework of the DGNN and Hybrid models, the BFS pathfinder is implemented as a bidirectional search originating simultaneously from both the root theory Q A and the target theory Q B . DGNN Pathfinder: It consists of a bidirectional A ∗ search algorithm with a heuristic function h(·) given by the output of the DGNN. This means that for any intermediate quiver states Q n and Q m originating from mutating Q A and Q B , respectively, the evaluation func- tions for the forward and backward frontiers are defined as: f fwd (Q n ) = g fwd (Q n ) + DGNN(Q n ,Q B ) f bwd (Q m ) = g bwd (Q m ) + DGNN(Q m ,Q A ), (5.5) with g fwd (Q A ) = g bwd (Q B ) = 0. The dequeued element at each step is determined by which evaluation function is smaller. This represents the first of a set of pathfinders in which the N prediction provides the policy to guide the pathfinder. As we will discuss in Section 6, despite being quite efficient and successful, the main problem with using the DGNN as a heuristic 28 is its non-monotonicity; i.e., the distance between the two quivers is not monotonically de- creasing. In fact, sometimes the heuristic delays the search queue, eventually degrading to a BFS algorithm. AGNN Pathfinder: Unlike the distance-guided approach, this implementation utilizes the probability distributions output by the Adviser GNN (AGNN) to parameterize the cost function g(Q) of a beam search with beam width B. 24 Given an initial theory Q A and a target theory Q B , the algorithm explores the space of valid Seiberg dualities starting from Q A . Following our empirical evaluation of the network’s classification accuracy in Section 4.2, we set the fixed beam width to B = 3. At each expansion step, the algorithm processes all quivers currently residing in the beam. For a given state Q n ∈ B t , the AGNN predicts the mutation probability distribution P (D k |Q n ,Q B ) across all admissible nodes. The algorithm selects the B actions with the high- est predicted probabilities and generates the corresponding mutant quivers Q n+1 = D k Q n . To formulate the trajectory evaluation as a cost minimization problem, each valid succes- sor state Q n+1 is assigned a cumulative trajectory cost g(Q n+1 ), defined as the accumulated negative log-probability of the mutation sequence: g(Q n+1 ) = g(Q n )− logP (D k |Q n ,Q B ),(5.6) where the root theory is initialized at g(Q A ) = 0. Since logP ≤ 0, highly probable mutations contribute minimally to the accumulated cost, while unlikely ones heavily increase the cost. To prevent cyclic exploration, the algorithm checks the best candidates against the WL hashes of all previously evaluated quivers. If the candidate’s structure has already been visited, the candidate is discarded. 25 Once possible loop branches have been pruned, the algorithm collects all surviving candidates generated across the current beam, sorts them in ascending order based on their cumulative trajectory cost g(Q n+1 ), and retains only the top B quivers to form the new beam B t+1 . The search procedure terminates successfully the moment a generated candidate is con- firmed to be isomorphic to the target theory Q B . Conversely, the search returns a failure status if the exploration reaches a predefined maximum depth or if the beam becomes empty due to dead-end mutations. 26 24 Because beam search permanently discards alternative branches without a backtracking mechanism, the probability that a forward beam originating from Q A and a backward beam from Q B intersect in the vast intermediate mutation space is vanishingly small. Therefore, unlike the A ∗ pathfinder, beam search is strictly implemented as a unidirectional search originating from Q A . 25 Recall that the beam search algorithm simultaneously tracks up to B active trajectories at each step. In this sense, it can be conceptualized as an A ∗ search with a truncated frontier, trading the guarantee of completeness for memory efficiency. 26 A dead-end occurs when all active trajectories in the beam are pruned. This scenario arises if the AGNN 29 Hybrid Pathfinder: To overcome the limitations of the DGNN A ∗ search and the AGNN beam search, we developed a Hybrid Pathfinder. This approach combines the distance esti- mates of the DGNN with the branch probabilities of the AGNN. Specifically, the AGNN acts as a soft pruning mechanism, discouraging exploration along improbable branches and pre- venting the algorithm from stalling in flat or non-monotonic regions of the DGNN landscape. Retaining the A ∗ priority queue provides a mechanism to backtrack if the AGNN predictions diverge. The Hybrid Pathfinder inherits the bidirectional exploration structure of the DGNN A ∗ search but modifies the step-cost accumulation using the policy output of the AGNN. For any valid mutation Q n+1 = D k (Q n ), the evaluation function f (Q n+1 ) is defined again as: f (Q n+1 ) = g(Q n+1 ) + λ DGNN DGNN(Q n+1 ),(5.7) where λ DGNN is a scaling hyperparameter balancing the heuristic weight, and the cumulative trajectory cost g(Q n+1 ) updates recursively via: g(Q n+1 ) = g(Q n ) + 1− λ AR logP (D k | Q n ,Q target ).(5.8) Here, the base step cost of 1 ensures admissibility by guaranteeing strictly increasing path costs, while the parameter λ AR ≥ 0 controls the penalty imposed on actions deviating from the optimal AGNN policy. When λ AR = 0, the expansion reverts to the standard DGNN A ∗ search; conversely, as λ AR → ∞ with λ DGNN = 0, the trajectory cost dominates, recovering the behavior of the AGNN pathfinder. LCA Pathfinder: This algorithm works as an A ∗ search without a heuristic function, where at each iteration, the pathfinder expands the active frontier in the direction that has the minimum path cost determined by hard-coded rules, not learned from training. At every state Q n , the cost of mutating it with D k is based on the change in the rank of the dualized node k. The cost function g(D k ) is determined by three weights, defined as c dec : c eq : c inc ∼ 1 : 10 : 100, such that the cost of every step is g(D k ) = c dec if rank(k) decreases, c eq if rank(k) is unchanged, c inc if rank(k) increases. (5.9) Because we penalize the increase of rank for both frontiers, the search does not attempt to find the direct path between two theories. Instead, the cost function acts as a “gravity well” that funnels both search trees towards the model that has the lowest rank. policy repeatedly favors mutations that lead to previously visited states. Because the algorithm systematically filters out historical structures to prevent infinite loops, an overconfident policy trapped in a cycle will cause all candidate branches to be discarded, leaving B t+1 =∅. 30 The motivation behind the LCA is to construct an algorithm that mimics the process that (human) physicists follow when guessing if two theories are related by Seiberg dualities: one tries to obtain a “simpler theory” from both quivers, reducing ranks and arrows between nodes, hoping to reach a common ancestor that generates the pair of quivers. This heuristic serves as a more robust baseline than standard BFS for benchmarking our neural pathfinders. By explicitly incorporating the structural “complexity” of the quiver, ranks, and total number of arrows, this algorithm effectively prioritizes structural simplification. This contrasts with a standard BFS, which is topologically agnostic and blindly treats every valid mutation as equally promising, regardless of its impact on the theory’s complexity. Hybrid LCA Pathfinder: Finally, a natural possibility is the integration of the LCA and Hybrid pathfinders. The idea is to further help the Hybrid pathfinder overcome the non- monotonicity and inversion problems: if the guidance of the NNs is unsure, the LCA can move the search towards the minimal quiver. The LCA cost function enters into modifying both the cost and heuristic functions of the Hybrid pathfinder. To this end, we modify the total edge cost for the priority queue, which accumulates as: g(Q n+1 ) = g(Q n ) + (1− λ det cost ) + λ detcost c det (D k )− λ AR logP (D k |Q n ,Q target ),(5.10) where c det (D k ) is defined similarly to (5.9), i.e. c det (D k ) = c det dec if rank(k) decreases, c det eq if rank(k) is unchanged, c det inc if rank(k) increases. (5.11) Similarly, the heuristic incorporates a penalty proportional to the total rank of the quiver to encourage moving along quivers with smaller ranks. The modified function is defined as: h(Q n ) = λ DGNN DGNN(Q n ,Q target ) + λ LCA P rank(Q n ) P rank(Q root ) ,(5.12) where Q root denotes the initial graph Q A or Q B depending on the active frontier. By adjusting the weights λ det cost , λ AR , λ DGNN , and λ LCA , this pathfinder averages be- tween a purely N-guided framework and the LCA pathfinder. However, this integration can cause conflicting search signals between the two pathfinders. The Hybrid pathfinder uses the guidance of DGNN and AGNN, which have learned to find the shortest path between two theories. This means that such a path might necessitate moving along trajectories where the quivers increase their ranks, completely avoiding the LCA. On the other hand, the LCA pathfinder works like an effective potential, moving precisely to find the theories with lower 31 ranks. If the LCA is integrated into the Hybrid with the same hierarchy of weights defined for the LCA pathfinder, the resulting algorithm would simply reduce to the LCA pathfinder without any guidance given by the NNs. To balance the two guiding policies of the pathfinder, we introduce hyperparameters designed to maximize the efficiency ratio while maintaining a 100% success rate. The opti- mizer converged to c det dec = 0.3, c deteq = 2.7, and c detinc = 3.1, alongside a global weight λ det cost = 1.3. Substituting these into the definition of c base yields effective base step costs of c basedec = 0.09, c baseeq = 3.21, and c baseinc = 3.73 for the priority queue. This calibration imposes a penalty for rank-increasing mutations that is of the same order of magnitude as the neural network coefficients while contributing negligibly to the cost of rank-reducing mutations. Consequently, the LCA penalty does not dominate the Hybrid cost function. Instead, it biases search toward smaller-rank theories when network confidence is low while allowing high-probability shortcuts to proceed uninterrupted. 5.3 Complexity Comparison We, now, compare the computational complexity of the N-guided pathfinders against the unguided bidirectional BFS baseline established in Section 5.1. In an unguided expansion, the number of explored states scales exponentially as O b d/2 avg. , where b avg. is the average valid branching factor and d is the minimum mutation distance between the theories. In the ideal limit of a monotonic heuristic, an A ∗ search expands only the states along the optimal trajectory, achieving linear complexityO(d). However, as shown in Section 4.1.1, the distance estimates generated by the DGNN exhibit non-monotonic fluctuations across mutation sequences, violating the formal conditions for algorithmic consistency. Because A ∗ maintains an exhaustive frontier of unexplored branches, a locally misleading heuristic can attract the search into non-optimal trajectories. Consequently, in the worst-case scenario, the computational cost not only degrades to the exponential scaling of the BFS baseline but can even exceed it, exhausting the step budget before recovering the true path. The same worst-case bound applies to the LCA baseline, which operates as an A ∗ search guided purely by a deterministic topological cost function. Because the LCA heuristic strictly penalizes increases in gauge ranks and arrow counts, the algorithm can easily become trapped in a local minimum of topological “complexity”. In this case, the search might stagnate in minimal-rank configurations, eventually exhausting its computational budget. Conversely, the AGNN beam search strictly confines the expansion tree to a maximum of B candidate quivers per level. The number of evaluated states is strictly bounded by O(B·b avg. ·d), thereby guaranteeing predictable linear scaling with respect to distance at the cost of a decrease in N accuracy. Indeed, this massive reduction in memory and execution 32 time sacrifices algorithmic completeness: as discussed in Section 4.2.1, if an incorrect AGNN prediction prunes the true path from the active beam, the algorithm permanently loses the ability to backtrack and fails to find a solution. The Hybrid Pathfinder overcomes the limitations of both by integrating localized pruning policies with a backtracking memory. First, by applying a Top-k filter based on the AGNN policy distributions, the effective branching factor is bounded by b eff. ≤ B ≪ b avg. . Second, embedding the negative log-likelihood penalty− logP (D k |Q n ,Q target ) into the step-cost g(Q) counteracts the non-monotonicity of the DGNN heuristic. Furthermore, because P ≤ 1, the modified step cost is bounded from below by g(Q n+1 )− g(Q n )≥ 1. This positivity guarantees that, even under conditions where both NNs produce inaccurate predictions, the Hybrid Pathfinder will degrade to a BFS-like exploration. Finally, the Hybrid LCA variant preserves these theoretical complexity bounds while introducing a bias that favors rank-reducing mutations when the N-guidance is uncertain. 5.4 Pathfinders Evaluation and Efficiency Metrics To quantify the performance and efficiency of the pathfinders, we compare them to baseline algorithms, namely the BFS search and the LCA pathfinder. Although the AGNN pathfinder performs a unidirectional search, we nevertheless decided to compare it with the bidirectional baselines to facilitate comparisons with the other models. Success Rate (SR): We define, for each pathfinder, the success rate as the fraction of test pairs for which the model successfully finds a path, i.e., N succ. , out of the total number of test pairs evaluated, i.e., N tot. , obtaining SR = N succ. N tot. .(5.13) Efficiency Ratio (ER): This indicates how successful the guidance coming from the cost and heuristic functions is at pruning the search space compared to a given baseline search. An ER > 1 indicates an active reduction in search complexity. We define two different baselines: the BFS search, representing an uninformed baseline, and the LCA search, serving as a physics-informed baseline. We thus compute the ERs against either the BFS explored nodes N BFS or the LCA ones N LCA . The ratio is evaluated against the baseline expectation: ER bidir = N baseline N exp ,(5.14) where N exp is the total number of nodes explored by search algorithms as in (5.4), and N baseline corresponds to N BFS, bidir or N LCA depending on the comparison. 27 Note that when N baseline = 27 For the LCA pathfinder, the comparison is obviously done over the BFS baseline. 33 N BFS , the ratio inherits the oscillations originating from the asymmetric bidirectional BFS expansion. These spikes at odd distances are an artifact of the unguided baseline and do not reflect the underlying behavior of the pathfinders, as confirmed by the absence of such oscillations when evaluating the efficiency against the LCA pathfinder. Guidance Efficiency (GE): This metric quantifies the precision of the neural network policy by measuring the search relative to the theoretical minimum number of steps required to discover a path of distance d. Because it is defined as a ratio of observed evaluations to the ideal minimum, it is strictly bounded by GE≥ 1, where values approaching 1 indicate a near-perfect trajectory with minimal branching redundancy. An optimal search would never expand branches outside the shortest path. Accounting for the immediate neighbors evaluated to confirm each step along the trajectory, the theoretical minimum number of explored states is b avg. at the root level and b avg. − 1 for the remaining d− 1 transition steps. The bidirectional guidance efficiency is thus defined as: GE bidir = N exp b avg. + (d− 1)(b avg. − 1) .(5.15) For the AGNN pathfinder, however, the search is performed by only keeping a beam width B. Under perfect guidance, the AGNN evaluates only the single root quiver at depth 0 and at most min(B,b avg. ) candidate branches at each subsequent depth level. The unidirectional guidance efficiency is therefore defined as GE unidir = N exp 1 + (d− 1) min(B,b avg. ) .(5.16) Effective Efficiency Ratio (EER): Because the evaluated models exhibit varying SR across different duality distances and quiver complexities, and since the underlying ERs are formulated differently for A ∗ and beam search architectures, a unified comparison requires a composite metric. We, therefore, define the Effective Efficiency Ratio as the direct product of efficiency and accuracy: EER = ER× SR.(5.17) This formulation naturally balances the trade-off between search speed and algorithmic com- pleteness, strictly penalizing pathfinders that fail to converge by assigning zero effective efficiency to unresolved pairs. The benchmarking evaluations described above are initially conducted on an in-distribution (ID) test set drawn from the theory families detailed in Appendix A.1. This assesses pathfinder performance when the neural networks operate within familiar topological domains encoun- tered during training. Furthermore, because the search algorithms actively transform theory structures through sequential Seiberg dualities, it is important to evaluate their robustness 34 against unfamiliar architectures. Consequently, we extend our benchmarks to the out-of- distribution (OOD) theories cataloged in Appendix A.2 to test the generalization capabilities of the learned policies. 6 Results In this section, we evaluate the computational efficiency and routing performance of the GNN pathfinders introduced in Section 5. We begin our analysis by benchmarking the models on the in-distribution (ID) test set detailed in Appendix A.1, evaluating their search capabilities within familiar domains. In Section 6.2, we investigate their generalization properties when applied to the out-of-distribution (OOD) architectures listed in Appendix A.2. To ensure statistical consistency across all evaluations, each experimental benchmark is performed on a randomized sample of 500 quiver pairs, systematically categorized by node count and Seiberg duality distance. 6.1 In-Distribution Theories Analysis 6.1.1 DGNN Pathfinder Performance We first examine the baseline performance of the A ∗ search, guided purely by the DGNN dis- tance heuristic, as formulated in Section 5.2. As illustrated in Figure 14, the DGNN pathfinder demonstrates exceptional reliability across the ID evaluation set, achieving an overall SR of 99.94%. The algorithm successfully resolves the vast majority of duality trajectories, exhibit- ing only minor performance degradation when traversing large duality distances. As illustrated in Figure 15a, the search efficiency of the DGNN pathfinder remains stable and does not degrade significantly at large distances. However, the guidance of the DGNN model in Figure 15b reduces when considering pairs of quivers separated by larger distances. This means that the pathfinder is slowly degrading to an unguided pathfinder. Despite this, the EER does not show signs of decreasing with distance; hence, it still outperforms a pure BFS algorithm. However, if we had the possibility to improve the guidance of the N, the performance of the pathfinder could, in principle, improve. One of the reasons why the DGNN pathfinder becomes unguided at large distances is due to the non-monotonicity of the heuristic function given by the output of the DGNN. One way to compensate for this is to try to modify the cost function so that the pathfinder can prioritize certain directions when exploring the possible nodes to dualize. The AGNN pathfinder, discussed below, addresses this issue by modifying the cost function. 35 (a) Performance against BFS Baseline DatasetPathfinderSR (%)EREER ID dataset DGNN99.94430.25429.99 AGNN76.12928.59706.84 Hybrid100.00790.08790.07 LCA100.00763.66763.66 OOD dataset in Figure 37 DGNN99.51134.38133.72 AGNN50.91305.15155.35 Hybrid99.84221.99221.63 LCA99.97176.13176.08 OOD dataset in Figure 38 DGNN99.34714.42709.74 AGNN78.641175.18924.21 Hybrid99.92881.67880.96 LCA100.00568.25568.25 (b) Performance against LCA Baseline DatasetPathfinderSR (%)EREER ID dataset DGNN99.940.830.83 AGNN76.122.041.55 Hybrid100.001.101.10 Hybrid LCA100.001.141.14 OOD dataset in Figure 37 DGNN99.510.540.53 AGNN50.911.030.52 Hybrid99.840.830.83 Hybrid LCA100.000.980.98 OOD dataset in Figure 38 DGNN99.340.830.83 AGNN78.641.911.50 Hybrid99.921.161.16 Hybrid LCA100.001.181.18 Table 2: Summary of pathfinder performances evaluated as the average across all distances and node counts in each dataset. The EER is the effective efficiency ratio obtained by multi- plying the success rate by the corresponding efficiency ratio as in (5.17). Table 2a shows the comparison of the pathfinders against the BFS baseline, while Table 2b shows the comparison of the pathfinders against the LCA baseline. 36 1.001.001.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.001.000.99 1.001.001.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.001.000.95 1.001.001.001.001.001.001.001.001.001.001.000.99 1.001.001.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.00 123456789101112 3 4 5 6 7 8 9 10 11 12 13 Distance Nodes 0 0.2 0.4 0.6 0.8 1 Figure 14: Success Heatmap 24681012 10 0 10 1 10 2 10 3 10 4 Distance EER K=3K=4 K=5K=6 K=7K=8 K=9K=10 K=11K=12 K=13BFS (a) Efficiency Ratio 24681012 0 50 100 150 Distance GE K=3K=4 K=5K=6 K=7K=8 K=9K=10 K=11K=12 K=13Opt (b) Guidance Efficiency Figure 15: Global performance of the DGNN pathfinder tested on theories in the testing set of DGNN. 37 6.1.2 AGNN Pathfinder Performance We next evaluate the unidirectional beam search guided by the AGNN policy. The main difference of the beam search compared to the A ∗ search is the limited frontier that, at each step of the path, is maintained. However, limiting the frontier, if the policy is not perfect (as anticipated in Section 4.2.1), may result in the pathfinder being unable to find a path at all. For this reason, the AGNN pathfinder will have a success rate lower than that of the DGNN A ∗ pathfinder. As a benchmark, we evaluate the AGNN pathfinder across the same test set used in the previous section. As illustrated in Figure 16, the overall SR for the AGNN pathfinder converges to approximately 76.12%. Nevertheless, the AGNN pathfinder achieves high search efficiency. 1.001.000.660.550.590.460.430.380.350.340.300.33 0.900.900.880.840.770.710.630.600.510.470.380.39 0.920.920.930.890.860.800.710.660.570.510.380.35 0.960.950.960.940.890.810.700.640.530.480.340.35 0.960.960.950.950.910.860.770.720.600.580.430.38 0.970.960.970.940.910.850.770.690.520.510.350.33 0.930.910.870.850.830.790.780.740.700.680.600.56 0.950.960.920.890.870.860.790.710.690.650.440.38 0.980.960.950.930.950.930.940.920.940.87 0.960.930.900.860.860.820.800.730.690.65 0.990.970.970.940.950.930.930.920.910.87 123456789101112 3 4 5 6 7 8 9 10 11 12 13 Distance Nodes 0 0.2 0.4 0.6 0.8 1 Figure 16: Success Heatmap The AGNN pathfinder is the only pathfinder in this work that is unidirectional. However, for consistency in the analysis with the other pathfinders, we still compare it with the bidi- rectional BFS baseline. Additionally, the beam search we performed set the beam width to B = 3. This is the reason why, in Figure 17a, for quivers with K = 3 and small distances, the performance of the AGNN pathfinder is worse than that of the BFS baseline since, by con- struction, the AGNN pathfinder explores all the nodes of the quiver at each step. However, for larger quivers, the AGNN pathfinder is an order of magnitude better than the bidirec- tional BFS baseline, despite its reduced SR. Furthermore, the efficiency of the pathfinder grows with the distance, even for small quivers. This sustained efficiency is directly explained by the GE in Figure 17b, which remains 38 stable across all distances. The AGNN pathfinder resolves the large-distance breakdown seen in DGNN pathfinder for large distances because the pathfinder has a more definite criterion for choosing the direction of the path. However, this gain in efficiency comes at the expense of success rate 24681012 10 0 10 1 10 2 10 3 10 4 Distance EER K=3K=4 K=5K=6 K=7K=8 K=9K=10 K=11K=12 K=13BFS (a) Effective Efficiency Ratio 24681012 1 2 3 Distance GE K=3K=4 K=5K=6 K=7K=8 K=9K=10 K=11K=12 K=13Opt (b) Guidance Efficiency Figure 17: Global performance of the AGNN pathfinder evaluated on theories in the testing set of AGNN. Combining the AGNN and DGNN policies into a single pathfinder aims to recover the success rate while maintaining high GE. 28 6.1.3 Hybrid Pathfinder Performance We now evaluate the Hybrid pathfinder, designed to integrate the complementary guidance policies of the AGNN and DGNN to address their limitations. The evaluation is once again performed on the same test set. The most immediate result is the recovery of a 100% SR across all tested duality distances, as shown in Figure 18. Figure 19a confirms that the ER outperforms the unguided BFS baseline across all dis- tances, becoming more efficient also for small quivers. We also see that the GE in Figure 19b is closer, in order of magnitude, to the one observed for the AGNN pathfinder, starting to diverge only at large distances. The reason for this behavior is that the Hybrid pathfinder keeps the whole frontier of nodes to dualize, as does the DGNN pathfinder, but each node 28 Recall that, by definition in (5.15), the GE is better the closer it is to unity. 39 1.001.001.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.00 123456789101112 3 4 5 6 7 8 9 10 11 12 13 Distance Nodes 0 0.2 0.4 0.6 0.8 1 Figure 18: Success Heatmap has a cost function given by the probability distribution predicted by the AGNN model. 29 For ID theories, the Hybrid pathfinder resolves the guidance inefficiency of DGNN and the success-rate drop of the AGNN pathfinder, giving the so-far most efficient N-guided pathfinder with a 100% SR in finding a path (if the path exists). 6.1.4 LCA Pathfinder Performance In this section, we discuss the performance of the LCA pathfinder over the BFS baseline. Despite the fact that this pathfinder does not have a N-guided policy, we test it on the same sample of theories that we used for the previous pathfinders to be able to compare its performance. As shown in Figure 20, the LCA pathfinder reaches an SR near 100%, so we can use the ER as a ratio of efficiency when we compare it with the BFS baseline. We see from Figure 21a that the LCA pathfinder reaches levels of performance comparable to the Hybrid pathfinder in Figure 19a. However, the way in which each pathfinder is designed is different. The LCA pathfinder finds a common ancestor with lower ranks, while the Hybrid pathfinder has two N-guided policies that try to reduce the distance between the quivers and predict the nodes to be dualized at each step. It is true that the NNs have been trained to minimize the differences between the quivers; however, the distances predicted or the probability distributions are generated on the resulting graphs processed by the Transformers. 29 In Figure 19b, we can see that the GE can become smaller than 1. This happens when, given a starting quiver Q A , multiple nodes exist that can be dualized to reach Q B , while (5.15) assumes that there is only one path connecting each pair. 40 24681012 10 0 10 1 10 2 10 3 10 4 Distance EER K=3K=4 K=5K=6 K=7K=8 K=9K=10 K=11K=12 K=13BFS (a) Efficiency Ratio 24681012 1 1.5 2 Distance GE K=3K=4 K=5K=6 K=7K=8 K=9K=10 K=11K=12 K=13Opt (b) Guidance Efficiency Figure 19: Global performance of the Hybrid pathfinder tested on theories in the testing set of DGNN and AGNN. 1.001.001.001.001.001.00 1.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.001.001.00 1.00 1.001.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.001.001.00 1.00 1.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.00 1.00 1.001.001.001.001.001.001.001.001.00 123456789101112 3 4 5 6 7 8 9 10 11 12 13 Distance Nodes 0 0.2 0.4 0.6 0.8 1 Figure 20: Success Heatmap 41 The fact that both pathfinders perform similarly is, in some sense, surprising because it seems that the NNs have learned some pattern to guide the path towards a common quiver. This raises the question of how much the Hybrid pathfinder improves upon the LCA baseline. 24681012 10 0 10 1 10 2 10 3 10 4 Distance EER K=3K=4 K=5K=6 K=7K=8 K=9K=10 K=11K=12 K=13BFS (a) Efficiency Ratio 24681012 0.8 1 1.2 1.4 1.6 Distance GE K=3K=4 K=5K=6 K=7K=8 K=9K=10 K=11K=12 K=13Opt (b) Guidance Efficiency Figure 21: The global performance of the LCA pathfinder tested on the ID dataset. 6.1.5 Comparing Pathfinders with LCA Baseline Although the BFS baseline is the only objective comparison to judge the efficiency of a pathfinder, since it represents the blind search one would perform when given two graphs, the LCA pathfinder represents the algorithm that a physicist would write for the task of finding the chain of dualities connecting two theories; this should be, from now on, the baseline against which to compare our N-guided pathfinders. In summary, DGNN achieves a high SR with degrading GE, whereas AGNN yields a lower SR with superior EER and GE. The Hybrid pathfinder reaches a 100% SR, good performance, and stable GE. In Figure 22, we show how those same EERs compare when evaluated against the LCA pathfinder as the baseline. The first observation in Figure 22 is that the scale of the EER is now linear instead of logarithmic. Unlike the BFS comparisons, no model outperforms the LCA pathfinder by orders of magnitude; only the AGNN pathfinder in Figure 22b has a∼ 10× speedup for small distances. The Hybrid pathfinder achieves a global average speedup of 1.09× over the LCA baseline, but its performance depends on the quiver size and the distance: 42 24681012 0.5 1 Distance EER (a) DGNN Pathfinder 24681012 0 2 4 6 8 Distance EER (b) AGNN Pathfinder 24681012 0.6 0.8 1 1.2 1.4 Distance EER K=3K=4 K=5K=6 K=7K=8 K=9K=10 K=11K=12 K=13LCA (c) Hybrid Pathfinder Figure 22: EER for the three pathfinders using LCA baseline in the ID dataset. 43 • For small quivers (K < 5), the LCA pathfinder is more efficient, while Figure 23b shows that the Hybrid pathfinder outperforms LCA on larger graphs (K ≥ 5), with its average EER scaling from 1.05× up to 1.18× for K ≥ 10. • Figure 23c shows an immediate advantage for the Hybrid model across distances d ∈ [2, 8], peaking at ∼ 1.21× for d ∈ [4, 5]. The EER drops below unity only for small quivers at very high distances, d > 10. This is corroborated by the absolute node savings (∆N exp = N LCA − N Hyb ) in Figure 23a, confirming that N-guidance becomes important as soon as the size of the quiver increases. 24681012 −20 0 20 Distance ∆ N exp K=3K=4 K=5K=6 K=7K=8 K=9K=10 K=11K=12 K=13 (a) Difference of the nodes explored 4681012 0.9 1 1.1 1.2 Nodes EER EER Avg (1.09×) (b) Average EER vs nodes 24681012 0.9 1 1.1 1.2 Distance EER EER Avg (1.09×) (c) Average EER vs distance Figure 23: Nodes explored and average EER per nodes and distances of the Hybrid pathfinder compared with the LCA pathfinder in the ID dataset We must determine whether the Hybrid model’s advantage over LCA stems from finding shorter paths or simply exploring fewer nodes. It could then be that the better performance (within limits) of the Hybrid pathfinder over the LCA that we are seeing in Figure 23 is only due to the pathfinder being able to find shortcuts instead of finding the path to the 44 common ancestor. For this reason, in Figure 24, we computed the percentage of times that the Hybrid pathfinder found a shorter or longer path than the LCA pathfinder. 30 On average, in the ID dataset, the Hybrid pathfinder found a shorter path only ∼ 3.4% of the time and found a longer path (due to the confusion of the N policies) ∼ 12.6% of the time. This means that ∼ 84% of the time, the Hybrid pathfinder found a path of the same length as the LCA pathfinder, and its speedup is purely due to the exploration of fewer nodes during each step, as shown in Figure 23a. Therefore, despite the different logic of guidance between the pathfinders, the Hybrid pathfinder behaves similarly to LCA 84% of the time, while it is also able to find shortcuts when possible. 31 345678910111213 0 5 10 15 20 25 30 Nodes Percentage (%) Shorter Longer Figure 24: Total percentage of path deviations (Shorter and Longer) found in ID dataset by Hybrid pathfinder compared to the LCA pathfinder. 6.1.6 Hybrid LCA Pathfinder Performance The purpose of the Hybrid LCA pathfinder is to find a balance between the physics-informed policy defined by the LCA and the N-guided policies of the Hybrid pathfinder. In this way, similarly to when we discussed the Hybrid pathfinder, we can try to improve the performance of the pathfinders: when the Hybrid pathfinder is struggling, i.e., at larger distances or smaller 30 By the length of a path, we mean the number of mutations connecting the pair of theories. 31 The presence of longer paths can be due to symmetries of the quivers that allow for multiple paths connecting the pair or the non-perfect guidance of the N policies. Nevertheless, it is interesting to note that despite the path being longer, the number of nodes explored by the Hybrid pathfinder is generally smaller than those explored by the LCA pathfinder. 45 quivers, the LCA pathfinder can dominate the contribution to the heuristic. On the other hand, when the LCA pathfinder explores too many nodes in order to find the common ancestor, the N can select a shorter and better path to improve overall efficiency. Figure 25 shows that the Hybrid LCA pathfinder improves upon, with the parameters discussed in the previous section, is an improvement compared to the corresponding performance of the Hybrid (DGNN–AGNN only) pathfinder, which we can see in Figures 22c, 30c and 31c. 24681012 1 1.2 Distance EER K=3K=4 K=5K=6 K=7K=8 K=9K=10 K=11K=12 K=13LCA Figure 25: EER for the Hybrid LCA pathfinder on the ID dataset By inspection of Figure 26, the curves now show a positive trend in average EER for both nodes and distances, surpassing the pure LCA pathfinder as soon as the pairs of quivers are large and sufficiently separated. Despite the average EER being only slightly improved compared to what is shown in Figure 23, the trends seem to have a slower decay, guaranteeing that it is better than pure LCA for longer distances than the Hybrid pathfinder was. Finally, testing reveals that there is no significant competition between the Hybrid and LCA policies in the search for a path, leading to the preference for finding a path over finding the shortest path. We see in Figure 27 that there is no significant increase in longer paths over the pure Hybrid pathfinder statistics shown in Figure 24. We believe that this is one of the signals that the optimization of the parameters has been successful and that the two policies are assisting each other in finding the path (rather than interfering with each other). 6.2 Out-of-Distribution Theories Analysis To probe the generalization capabilities of our NNs, we test them on the out-of-distribution (OOD) quiver families detailed in Appendix A.2, comparing our neural pathfinders directly against the BFS and LCA baselines. This OOD dataset encompasses both finite- and infinite- mutation type quivers, i.e., quivers with finite or infinitely many mutants, as well as quivers that violate gauge anomaly cancellation conditions. Although such anomalous configurations lack a physical realization as 4D N = 1 supersymmetric gauge theories, mutations are fun- damentally well-defined combinatorial operations on directed graphs. Consequently, these 46 2 4681012 0 10 20 30 Distance ∆ N exp K=3K=4 K=5K=6 K=7K=8 K=9K=10 K=11K=12 K=13 (a) Difference of the nodes explored 4681012 1 1.1 1.2 Nodes EER EER Avg (1.13×) (b) Average EER per nodes 24681012 1 1.05 1.1 1.15 1.2 Distance EER Avg (1.13×) (c) Average EER per distances Figure 26: Nodes explored and average EER per nodes and distances of the Hybrid LCA pathfinder compared with the LCA pathfinder in the ID dataset 47 345678910111213 0 5 10 15 20 25 30 Nodes Percentage (%) Shorter Longer Figure 27: Comparison of paths found in ID dataset by Hybrid LCA pathfinder when com- pared to the LCA pathfinder. “unphysical” examples provide a rigorous mathematical stress-test to evaluate the topologi- cal robustness of our pathfinder algorithms beyond familiar physical constraints. 32 6.2.1 Success Rates on OOD Theories We now analyze the performance of the various N-guided pathfinders when tested on OOD theories. The heatmaps for the SRs for the theories in Figure 37 are shown in Figure 28. We see that the DGNN pathfinder has an SR over 97% in Figure 28a, proving the ability of DGNN to guide the pathfinder regardless of whether that theory was part of the training set. This, however, does not apply to the AGNN pathfinder, which once again struggles to keep the SR up for large distances. For special sorts of OOD theories, we see in Figure 28b that the AGNN pathfinder fails completely due to the special nature of theories such as Q1, Q2, Q3, Q13, and Q15 in Figure 37 that either do not admit any non-anomalous realizations or admit only a finite number of mutations. The struggle of the AGNN model to guide the pathfinder is also reflected at the level of the Hybrid pathfinder, which failed to find the path for quivers with a large number of nodes, despite the DGNN pathfinder managing to 32 Although these anomalous theories are inconsistent when interpreted as 4D gauge theories, they are consistent as quiver quantum mechanics theory, and appear, for example, as the BPS quivers for various 4D N = 2 QFTs. In that context, mutations correspond to moving from one BPS chamber to another (see e.g., [37]). 48 do so alone. The success rate of 100% is, however, achieved when the Hybrid pathfinder is augmented with the LCA policy, as shown in Figure 28d. 1.001.001.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.001.001.00 1.001.00 1.001.001.00 1.001.001.001.001.001.001.001.001.000.970.960.90 1.001.001.001.001.00 1.001.001.001.001.001.00 1.001.001.001.001.001.00 1.001.001.001.001.001.00 1.001.001.001.001.001.00 1.001.001.001.001.001.00 1.001.001.001.001.001.00 1.001.001.001.001.001.00 1.001.001.001.001.001.00 123456789101112 3 4 5 6 7 8 9 10 11 12 13 Distance Nodes 0 0.2 0.4 0.6 0.8 1 (a) DGNN Pathfinder 1.001.000.670.610.480.470.400.410.320.320.280.29 0.830.820.670.610.540.490.390.350.310.280.190.19 1.000.00 0.500.000.00 0.830.790.750.660.640.620.630.620.490.390.320.28 0.500.000.000.000.00 0.170.000.000.000.000.00 0.670.000.000.000.000.00 0.330.000.000.000.000.00 0.330.000.000.000.000.00 0.830.000.000.000.000.00 0.670.000.000.000.000.00 0.500.000.000.000.000.00 0.670.000.000.000.000.00 123456789101112 3 4 5 6 7 8 9 10 11 12 13 Distance Nodes 0 0.2 0.4 0.6 0.8 1 (b) AGNN Pathfinder 1.001.001.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.001.000.99 1.001.00 1.001.001.00 1.001.001.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.00 1.001.001.001.000.000.00 1.001.001.001.000.000.00 1.001.001.001.000.000.00 1.001.001.001.000.000.00 1.001.001.001.000.000.00 1.001.001.001.000.000.00 1.001.001.001.000.000.00 1.001.001.001.000.000.00 123456789101112 3 4 5 6 7 8 9 10 11 12 13 Distance Nodes 0 0.2 0.4 0.6 0.8 1 (c) Hybrid Pathfinder 1.001.001.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.001.001.00 1.001.00 1.001.001.00 1.001.001.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.00 1.001.001.001.001.001.00 1.001.001.001.001.001.00 1.001.001.001.001.001.00 1.001.001.001.001.001.00 1.001.001.001.001.001.00 1.001.001.001.001.001.00 1.001.001.001.001.001.00 1.001.001.001.001.001.00 123456789101112 3 4 5 6 7 8 9 10 11 12 13 Distance Nodes 0 0.2 0.4 0.6 0.8 1 (d) Hybrid LCA Pathfinder Figure 28: Success Heatmap for Theories in Figure 37 (OOD dataset). The blank white regions are a consequence of these theories being of finite mutation type, i.e., there is an upper bound on the distance between any pair of theories. The OOD theories in Figure 38 admit infinitely many mutations, so they are closer to the theories considered during the training, although they were never seen by the NNs during the training. However, we see that the pathfinders behave similarly to the analysis done for each pathfinder in Section 6.1: the DGNN pathfinder has a SR of over 90%, while the AGNN pathfinder struggles to find a path at large distances, as we see in Figures 29a and 29b. The Hybrid pathfinder has almost perfect SR due to the AGNN policy failing at large distances, but the perfect score is obtained by the Hybric LCA pathfinder, as we show in Figures 29c and 29d. 49 1.001.001.001.001.001.001.001.001.000.980.920.77 1.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.000.99 1.001.001.001.001.001.001.001.001.001.00 123456789101112 7 8 9 10 11 Distance Nodes 0 0.2 0.4 0.6 0.8 1 (a) DGNN Pathfinder 0.850.890.820.810.790.710.620.540.410.330.260.23 0.990.990.990.980.960.910.830.770.660.59 0.990.990.960.940.910.900.840.800.750.76 0.990.980.950.910.780.730.650.570.490.47 1.000.990.970.950.920.870.810.800.670.65 123456789101112 7 8 9 10 11 Distance Nodes 0 0.2 0.4 0.6 0.8 1 (b) AGNN Pathfinder 1.001.001.001.001.001.001.001.001.001.000.990.97 1.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.00 123456789101112 7 8 9 10 11 Distance Nodes 0 0.2 0.4 0.6 0.8 1 (c) Hybrid Pathfinder 1.001.001.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.00 1.001.001.001.001.001.001.001.001.001.00 123456789101112 7 8 9 10 11 Distance Nodes 0 0.2 0.4 0.6 0.8 1 (d) Hybrid LCA Pathfinder Figure 29: Success Heatmap for Theories in Figure 38 (OOD dataset) 50 6.2.2 Efficiency on OOD Theories The efficiencies of the N-guided pathfinders on the OOD datasets, when compared to the LCA baseline, show interesting results. We can see in Figure 30 that when the number of mutations of the quivers is finite or the theories are anomalous, all pathfinders perform worse than the LCA pathfinder in general ground. However, the worse performance of the Hybrid LCA is around the LCA performance. We believe that this is a sign that if the pairs of theories are not known to admit infinitely many mutations, the use of the Hybrid LCA pathfinder could still be more efficient, since it will only be outperformed by the LCA when the theories are anomalous. Moreover, beating the LCA pathfinder on the datasets we generated, as mentioned already in the previous sections, is the most difficult test for our pathfinders because the LCA pathfinder follows the same criterion (but opposite) that we used to generate the dataset in the first place, as explained in Section 3.1. Therefore, the datasets, in particular when the number of allowed mutations is finite, are biased towards pairs that share a “lowest common ancestor”, by construction. If, instead, we consider theories that admit an infinite number of mutations, such as those in Figure 38, the N-guided pathfinders behave similarly to how they performed on the ID datasets. We see in Figure 31 that, in particular, the Hybrid LCA pathfinder has an LCA baseline that is usually below the efficiency of the N-guided pathfinders. The overall performance is also consistent with the results found in the ID dataset, with a peak efficiency around 1.6× faster than the LCA pathfinder, showing symptoms of degradation for large distances. Overall, however, we can conclude that the NNs guidance is not specific to the theories that were in the dataset; instead, it is based solely on the number of nodes and distances that were observed during the training of the NNs. For this reason, in Section 7, we study the limitations of the pathfinders, i.e., where they start to exhibit failure modes, to understand how long we can achieve better performance than the LCA pathfinder given a certain training of the NNs. 7 Failure Modes So far, we have tested the performance of the different pathfinders against each other, analyz- ing the pros and cons of each strategy and trying to maximize efficiency while maintaining the highest possible success rate. In general, we see that all pathfinders, even those that are N-guided, do not distinguish much between ID and OOD theories, maintaining simi- lar efficiencies when compared to a BFS or LCA baseline. From the discussion in Section 6, the Hybrid LCA pathfinder appears to be the most consistent, outperforming pure LCA every time when the theories admit an infinite number of mutations. However, even for these theories, the EER starts to decrease when small quivers are very far apart in terms of the number of mutations. The purpose of this section is to understand what the common failure 51 24681012 0 0.5 1 Distance EER (a) DGNN Pathfinder 24681012 1 2 3 4 Distance EER (b) AGNN Pathfinder 24681012 0.4 0.6 0.8 1 Distance EER (c) Hybrid Pathfinder 24681012 0.7 0.8 0.9 1 Distance EER (d) Hybrid LCA Pathfinder LCAK=3K=4K=5K=6 K=7K=8K=9K=10K=11 K=12K=13K=14K=15K=16 Figure 30: EER for the pathfinders using LCA baseline in the OOD dataset of Figure 37. 52 24681012 0.5 1 Distance EER (a) DGNN Pathfinder 24681012 2 4 Distance EER (b) AGNN Pathfinder 24681012 0.6 0.8 1 1.2 1.4 Distance EER (c) Hybrid Pathfinder 24681012 0.6 0.8 1 1.2 1.4 Distance EER (d) Hybrid LCA Pathfinder LCAK=7K=8K=9K=10K=11 Figure 31: EER for the pathfinders using LCA baseline in the OOD dataset of Figure 38. 53 modes and limitations (i.e., “the breaking point”) of the pathfinders we propose are, given the training of the NNs over theories with at most K train nodes and at a distance D train from each other. In particular, we want to find the turning point when the efficiency ratio decreases and the maximum complexity after which the pathfinders degrade worse than either BFS or LCA baselines. For this analysis, we considered only theories with an infinite number of mutations, as found in the ID dataset and Figure 38 of the OOD dataset. We introduce two new parameters: C = d true log 10 K ,W = N exp d true + 1 EW = N exp d true + 1 1 SR .(7.1) The first is the logarithm of the complexity K d true for a pair of quivers with K nodes at a distance d true from each other, as saved from the database using BFS. The distance d true is not the distance found by the pathfinders; rather, it is the reference distance we save when generating the database, as explained in Section 3.1. The parameter W (resp. EW) is the “wandering” (resp. “effective wandering”), which measures how many nodes the model actually explored compared to the minimum required if it made zero mistakes. Throughout this work, and with the checkpoints available on the GitHub repository, the NNs have been trained up to distance 12 for theories in the ID dataset with up to 10 nodes and up to distance 10 for theories in the ID dataset with up to 13 nodes. The complexities associated with these extremal values of training are C 1 = 12 log 10 10 = 12, C 2 = 10 log 10 13 = 11.1394.(7.2) Therefore, despite being very close, we consider C 1 = C train as the reference complexity associated to the pathfinders for the moment. We can then plot in Figure 32 how the SR, EER (with respect to LCA baseline), and EW differ across complexity. Interestingly, we see in Figure 32a that the model, which started degrading even at small complexities, is the AGNN pathfinder. On the other hand, Figure 32b presents it as the best model. This is consistent with the discussion that the efficiency of the AGNN is unmatched, but the lack of backtracking makes it unreliable for large complexities. However, the most interesting plot is Figure 32c, where we see that the pathfinder that “wandered” more is actually the DGNN, while all other models remained consistent, as shown in Figure 32d. 33 In fact, the EER and EW of all models except for DGNN are better than the LCA baseline. The models show no signs of degradation even at the largest complexities used for training, i.e., C train . In order to find the breaking point at which the NNs’ guidance causes the pathfinder to degrade and become worse than an unguided LCA, we decided to artificially generate a 33 By definition, W and EW are “better” the closer they are to 1. 54 024681012 0 0.5 1 C SR (a) Success Rate 024681012 10 0 10 1 10 2 10 3 10 4 C EER (b) Effective Efficiency Ratio 024681012 10 1 10 2 C EW (c) Effective Wandering 024681012 10 0.5 10 1 C EW (d) Effective Wandering without DGNN AGNNDGNN Hybrid LCA Hybrid LCA C train Figure 32: Success Rate, Effective Efficiency Ratio (with respect to LCA baseline) and Effec- tive Wandering for theories when considering NNs trained up to distance 12 and 13 nodes. 55 situation where we could push the pathfinders to large complexities while training the N on smaller ones. We therefore re-trained the NNs to a maximum complexity C train = 6 log 10 9≃ 5.73 and tested the pathfinders up to a complexity C = 13.37, obtained by considering theories in the ID dataset with K = 13 at distance 12. 34 The results are shown in Figure 33. 024681012 0 0.5 1 C SR (a) Success Rate 024681012 0 1 2 C EER (b) Effective Efficiency Ratio 024681012 10 0 10 1 10 2 10 3 10 4 10 5 C EW (c) Effective Wandering 024681012 10 0 10 1 10 2 C EW (d) Effective Wandering without DGNN AGNNDGNN Hybrid LCA Hybrid LCA C train Figure 33: Success Rate, Effective Efficiency Ratio (with respect to LCA baseline) and Effec- tive Wandering for theories when considering the NNs trained up to distance 6 on theories with up to K = 9. This time, we clearly see that the pathfinders start to degrade when the C train threshold is passed. In fact, by focusing on the EER and EW for each model, we can predict the complexity value after which the LCA pathfinder will perform better than the N-guided 34 We considered theories in the ID dataset with up to K = 14 and up to distance 10, i.e., with complexity C = 11.46. 56 pathfinders: C peak (SR peak )C pred (SR pred )C pred /C train EEREWEEREWEEREW DGNN4.00 (1.00)4.00 (1.00)4.00 (1.00)5.90 (1.00)0.701.03 AGNN6.48 (0.78)10.79 (0.33)48.02 (0.00)18.61 (0.10)8.393.25 Hybrid6.00 (1.00)8.00 (1.00)8.27 (1.00)8.92 (1.00)1.441.56 Hybrid LCA 6.00 (1.00)8.00 (1.00)11.25 (1.00)10.93 (1.00)1.971.91 (7.3) We see that the DGNN pathfinder is the one degrading even before reaching C train , while the AGNN pathfinder appears to be the most efficient for larger complexities, but the predicted SR is rapidly dropping to zero. The most interesting pathfinders are the Hybrid and Hybrid LCA pathfinders, which are able to perform better than LCA for complexities 1.46 to 1.97× larger than the training complexity, respectively. Also, note that the breaking point of the Hybrid and Hybrid LCA occurs before their success rates drop below 100%, meaning that they will be reliable for the entire duration they perform better than LCA. We believe this to be a good estimate of how long the N-guided pathfinders can outperform the LCA pathfinder. 8 Conclusions Dualities between quantum field theories provide important access to deep non-perturbative phenomena. It is therefore natural to ask: given two QFTs, could it be that they are actually dual? If so, what is the sequence of simple duality operations that connects the two theories? In practice, this is a challenging question and is amenable to machine learning techniques. Given the encouraging success of Transformers in tackling complex problems in string theory [38–41] and conformal field theories [42], it is natural to ask whether this architecture can be similarly incorporated into the duality pathfinder we aim to develop. In this paper, we have constructed a dataset of dual QFTs with a seed theory given by a D3-brane probing a toric singularity. We used this dataset to train different graph neural networks. Combining this with pathfinder algorithms, we have studied the computational complexity of connecting Seiberg dual theories and established that a hybrid approach tends to be both more efficient and more accurate. In the remainder of this section, we discuss some potential avenues for future investigation. In this work, we have investigated the computational complexity of dualities between quiver gauge theories. From a holographic standpoint, such sequences of dualities prominently appear in warped throats of a string compactification, where motion down the throat triggers a specific sequence of duality moves. It would be interesting to investigate whether this notion 57 of computational complexity matches more precisely with other bulk proposals of complexity (see, e.g., [43–45]). 35 An ambitious aim for the future would be to classify all physically distinct quiver gauge theories. To accomplish this, one would also need to have a sharp understanding of how many quiver gauge theories are in fact dual (i.e., to identify the equivalence classes). We anticipate using this work to tackle this classification problem. We have primarily focused on 4D quiver gauge theories and their Seiberg duals, but we anticipate that the same techniques apply more broadly. It would be interesting to study the computational complexity of dualities in other well-motivated settings, including dualities of string, M-, or F-theory compactifications, as well as quantum field theories in other spacetime dimensions. At a practical level, it would be interesting to increase both the size of the training dataset and the complexity of the network architectures used to establish Seiberg dual pairs. In particular, studying possible scaling laws in the degree of inference as a function of dataset and network size would be instructive. Given recent advancements in large language models (LLMs), it is natural to ask whether the custom GNN architecture used here is really necessary. Said differently, can one simply provide two dualized quivers and ask an LLM to find the most efficient and accurate path between these theories? This provides a natural class of benchmark problems for studying the capabilities of state of the art inference models, but at the same time, it is difficult to provide a completely “fair comparison” since it can be challenging to properly compare the computational resources of an LLM with the specific (quite modest) computational budget used in the present work. 36 Acknowledgements We thank F. Carta, M. Danese, S. De, C. Lawrie, H. Lee, and D.S. Park for helpful com- ments and discussions. The work of JJH and SNM is supported by DOE (HEP) Award DE-SC0013528, BSF grant 2022100, and a University Research Foundation grant at the University of Pennsylvania. AM and GS are supported in part by DOE (HEP) Awards DE- SC0017647 and DE-SC0023719. This work was performed in part while JJH, AM, and GS were at the Aspen Center for Physics, which is supported by National Science Foundation grant PHY-2210452. JJH and GS especially thank the organizers of the Aspen workshop “Theoretical Physics for Artificial Intelligence” for providing a stimulating atmosphere that directly led to the start of this project. SNM and AM thank the Simons Center for Geometry 35 It would also be interesting to study various generalizations of complexity distance measures of the sort introduced in [46]. 36 Can one truly put a price tag on finding all Seiberg dual quiver gauge theories? We think not. 58 and Physics for the hospitality during the “23rd Simons Physics Summer Workshop: Theory, Experiment and the Emerging New Physics”. The code produced for this work was written with the help of Gemini Pro 3.1. A Seed Theories for Dataset Generation In this appendix, we describe in more detail the seed theories used for dataset generation. We focus primarily on D3-branes probing toric CY singularities. See [20] and references therein for additional discussion of the construction and study of these theories. A.1 In-Distribution Theories We trained our NNs on five families of 4d N = 1 SCFTs that arise from D3-branes probing toric CY threefolds: 1. The 4dN = 1 SCFTs on D3-branes probing C 3 /Z (1,1,−2) n , with n being odd and orbifold action (1, 1,−2) on the C 3 coordinates. These are orbifolds of N = 4 SU(N ) SYM, usually characterized by three chiral fields X,Y ,Z in the adjoint of SU(N ), with a superpotential term W = Tr (X[Y,Z]) ,(A.1) where the trace is over color indices. The action of the orbifold leads to a theory with n gauge nodes SU(N i ), with chiral fields in the bifundamental representation of SU(N i )× SU(N j ). Using the same notation of the parent theory, they can be grouped in X i|i−2 mod n , Y i|i+1 mod n , Z i|i+1 mod n ,(A.2) while the superpotential becomes W C 3 /Z n = n X i=1 Tr X i|i−2 Y i−2|i−1 Z i−1|i − X i|i−2 Z i−2|i−1 Y i−1|i .(A.3) Figure 34 shows representative C 3 /Z n quivers. For N training, we used odd orbifolds with 3≤ n≤ 13. 2. The 4d N = 1 SCFTs on D3-branes probing C(dP 1 ), where by C(.) we mean the complex cone over the dP 1 surface. This theory is characterized by four gauge groups and bifundamental fields X 14 , Y 14 , X 21 , X 24 , X 31 , X 32 , Y 32 , X 43 , Y 43 , Z 43 ,(A.4) subjected to the superpotential W C(dP 1 ) =X 14 Y 43 Y 32 X 21 − Y 14 Y 43 X 32 X 21 + X 24 Z 43 X 32 − X 24 X 43 Y 32 + X 31 Y 14 X 43 − X 31 X 14 Z 43 . (A.5) 59 3 3 3 (a) C 3 /Z 3 2 2 2 2 2 (b) C 3 /Z 5 2 2 2 2 2 2 2 (c) C 3 /Z 7 Figure 34: Examples of C 3 /Z n quivers theories. The corresponding superpotential can be derived from (A.3). The quiver for this theory is shown in Figure 35a. 3. The 4d N = 1 SCFTs on D3-branes probing C(dP 2 ). This theory is characterized by five gauge groups and bifundamental fields X 15 , Y 15 , X 21 , X 25 , X 31 , X 32 , X 42 , X 43 , X 53 , X 54 , Y 54 ,(A.6) subjected to the superpotential W C(dP 2 ) =X 15 X 53 X 31 − X 25 X 53 X 32 + X 25 X 54 X 42 − X 31 Y 15 X 54 X 43 − X 15 Y 54 X 42 X 21 + X 21 Y 15 Y 54 X 43 X 32 . (A.7) The quiver for this theory is shown in Figure 35b. 4. The 4d N = 1 SCFTs on D3-branes probing C(dP 3 ). This theory is characterized by six gauge groups and bifundamental fields X 15 , X 16 , X 21 , X 26 , X 31 , X 32 , X 42 , X 43 , X 53 , X 54 , X 64 , X 65 ,(A.8) subjected to the superpotential W C(dP 3 ) =X 15 X 53 X 31 − X 15 X 54 X 42 X 21 + X 26 X 64 X 42 − X 16 X 64 X 43 X 31 − X 26 X 65 X 53 X 32 + X 16 X 65 X 54 X 43 X 32 X 21 . (A.9) The quiver for this theory is shown in Figure 35c. 5. The 4d N = 1 SCFTs on D3-branes probing the C(Y p,q ) singularities. This family of theories is characterized by 2p gauge groups and bifundamental fields that can be categorized into SU(2) doublets U α , V α , and singlets Z, Y . With gauge node indices understood modulo 2p, they can be grouped as U α 2k−1|2k (k = 1...p), V α 2k|2k+1 (k = 1...q), Z 2k|2k+1 (k = q+1...p),(A.10) 60 2 3 2 (a) C(dP 1 ) 2 2 (b) C(dP 2 )(c) C(dP 3 ) Figure 35: Quivers of C(dP n ) theories. The corresponding superpotentials are in Eqs. (A.5), (A.7) and (A.9). Y 2k+1|2k−1 , Y 2k+2|2k (k = 1...q), Y 2k+2|2k−1 (k = q + 1...p),(A.11) where α = 1, 2 is the doublet index. These fields enter to a superpotential as W C(Y p,q ) = q X k=1 ε αβ U α 2k−1|2k V β 2k|2k+1 Y 2k+1|2k−1 − V α 2k|2k+1 U β 2k+1|2k+2 Y 2k+2|2k + p X k=q+1 (−1) k−q−1 ε αβ U α 2k−1|2k Z 2k|2k+1 U β 2k+1|2k+2 Y 2k+2|2k−1 . (A.12) Some examples of C(Y p,q ) quivers we considered in our database are shown in Figure 36. There are some known identities for the C(Y p,q ) theories, such as C(Y 2,1 )≃ C(dP 1 ), or C(Y p,0 )≃C/Z p , where C is the conifold theory. Moreover, C(Y p,p )≃ C 3 /Z 2p . Since we want to avoid 2-loops in the quiver, we exclude the latter theories, and C(Y 1,0 )≃C. For the training of the NNs we considered 2≤ p≤ 6 and 1≤ q ≤ 6 with q ̸= p, so that the total number of nodes in the quivers was up to 12. 2 2 2 2 (a) C(Y 3,1 ) 2 2 22 2 (b) C(Y 3,2 ) 2 2 2 22 2 2 (c) C(Y 4,3 ) Figure 36: Examples of C(Y p,q ) quivers theories. The corresponding superpotential can be derived from (A.12). 61 A.2 Out-of-Distribution Theories To evaluate the neural network’s generalization to unseen theories, we tested it on the fam- ily of quivers Q1–Q15 from [13] (Figure 37), using non-anomalous rank assignments when possible. We also evaluated quivers from D3-branes probing non-toric CY threefolds [47, 48], shown in Figure 38. · n− 2 (a) Q1 · n− 4 (b) Q2(c) Q3 2 2 2 2 (d) Q4 or C(F 0 )≃ C(Y 2,0 ) 3 3 2 (e) Q5 2 4 32 (f) Q6 3 3 3 3 (g) Q7 2 2 5 33 (h) Q8 3 33 (i) Q9 or C 3 /Z 3 . 2 4 22 (j) Q10 3 2 (k) Q11(l) Q12 2 22 (m) Q13 2 2 2 2 2 2 2 (n) Q14(o) Q15 Figure 37: Quiver families considered in [13]. We comment that some of these quivers are intrinsically anomalous when inerpreted as 4D N = 1 QFTs (i.e., inconsistent), but are perfectly consistent when viewed as the BPS quivers of 4D N = 2 QFTs. As such, they provide an interesting testing ground for studying mutations of quivers. 62 2 5 2 2 2 2 2 (a) C(dP 4 ) 22 22 (b) C(dP 5 ) (c) C(dP 6 ) 22 4 (d) C(dP 7 ) 33 3 (e) C 3 /∆ 27 Figure 38: Further quiver theories used to test the N. Nodes with a rank other than one are specified whenever it is necessary to make the theory non-anomalous. 63 B Overview of Graph Neural Networks Graph Neural Networks (GNNs) process data natively structured as graphs by assigning features to each node and iteratively updating these representations over successive layers. This refinement is regulated by a localized neighborhood aggregation mechanism known as message passing [33]. In the following, we provide a brief review of this technique. B.1 Message Passing Paradigm In a standard Message Passing Neural Network (MPNN) [33], the network updates the hidden representation h i of node i at layer ℓ by aggregating messages from its neighbors j ∈ N (i). The unactivated update rule commonly takes the form: m (ℓ+1) i = W (ℓ) self h (ℓ) i + X j∈N (i) W (ℓ) msg h (ℓ) j ,(B.1) where W self and W msg are learnable weight matrices. After aggregating the messages, the network normalizes the features. A frequent choice is Layer Normalization [49], which stan- dardizes a vector m i ∈ R H using its mean μ i and variance σ 2 i , and shifts it with learnable parameters γ,β: LayerNorm(m i ) = γ⊙ m i − μ i p σ 2 i + ε + β ,(B.2) where ⊙ is the Hadamard product, i.e., element-wise multiplication. The normalized vector is then passed through an activation function. Since Layer Normal- ization can shift feature values to be negative, the Leaky Rectified Linear Unit (LeakyReLU) [50] is often used to prevent “dead” neurons: LeakyReLU(x) = ( xif x≥ 0, αx if x < 0. (B.3) The parameter α is a small positive constant. Stacking several message passing layers allows a node to gather information from nodes multiple hops away in the graph. B.2 Graph Transformers and Global Context While standard MPNNs effectively capture local graph topology, their ability to gather long- range information is limited. In particular, the phenomenon of oversmoothing [51], where node representations become indistinguishable after too many layers, remains one of the primary architectural bottlenecks of deep MPNNs. To solve this, architectures like GraphGPS [34] and Graph Transformers combine local message passing with global attention mechanisms [21]. 64 In a Graph Transformer layer, the sequence of node embeddings X ∈ R K×H is projected into Queries (Q), Keys (K), and Values (V ) using learnable linear transformations: Q = XW Q , K = XW K , V = XW V .(B.4) The attention weights determine how much focus each node places on every other node in the graph: Attention(Q,K,V ) = softmax QK T √ d k V .(B.5) To evaluate different types of relationships simultaneously, the network splits the embed- dings into multiple independent heads using Multi-Head Attention (MHA). The outputs are concatenated and projected back to the original dimension: MHA(X) = Concat(head 1 ,..., head h )W O .(B.6) A single Transformer encoder layer typically wraps this operation with position-wise Feed- Forward Networks (FFN) and residual connections: ̃ X = LayerNorm(X + MHA(X)), X out = LayerNorm( ̃ X + FFN( ̃ X)). (B.7) C Overview of Search Algorithms In this Appendix, we review some generalities about A ∗ Search and Beam Search pathfinders. C.1 A ∗ Search Algorithms A ∗ is an informed search algorithm optimized for weighted graphs: starting from a specific node, the algorithm seeks a path to a given final node by minimizing a specified cost criterion. The algorithm identifies the trajectory by maintaining a tree of paths originating from the starting point, extending them until the target node is reached. In contrast to uninformed search strategies like Dijkstra’s algorithm, A ∗ determines which path to extend based on an evaluation function f (n) that it minimizes. This function is defined as f (n) = g(n) + h(n),(C.1) where n designates the next node on the path, g(n) denotes the accumulated cost of the path from the starting node to n, and h(n) represents the heuristic function estimating the cost of the cheapest remaining path from n to the goal. Typically, the A ∗ implementation relies on a priority queue known as the frontier. At each iteration, the node with the lowest f (·) value is selected (dequeued) from the frontier, compared against the target, and, if it differs, 65 all its neighbors are added to the queue. Crucially, these new states supplement rather than overwrite the existing frontier. Preserving previously explored branches is essential to allow the algorithm to backtrack if a promising path reaches a dead-end. This cycle of selecting the lowest-cost nodes and expanding the frontier repeats continuously until the dequeued node coincides with the target. We now see how bidirectional A ∗ search (illustrated schematically in Figure 13a and detailed in Algorithm 2) is implemented to explore the space of quiver mutations in Section 5.1. The search expands two frontiers simultaneously: a forward frontier (Q fwd ) originating from Q A and a backward frontier (Q bwd ) originating from Q B . For any intermediate state Q n reached from Q A or Q m reached from Q B , the forward and backward evaluation functions are defined as: f fwd (Q n ) = g fwd (Q n ) + h fwd (Q n ), f bwd (Q m ) = g bwd (Q m ) + h bwd (Q m ),(C.2) where g(·) and h(·) depend on the specific pathfinder considered in Section 5.2. At each iteration, the algorithm prioritizes the expansion of whichever frontier currently holds the minimum estimated path cost min(f ). When a quiver is selected for expansion, its valid Seiberg mutants are generated and compared against the states already visited by the opposing frontier. C.2 Beam Search Algorithms While A ∗ employs a selection strategy, prioritizing the state with the lowest f -score at each iteration, it preserves completeness by maintaining the entire frontier of unexplored branches. This memory allows the algorithm to backtrack from dead-ends, guaranteeing that a valid path will be found in a finite number of steps if one exists. In contrast, we now introduce a Beam Search pathfinder, which is a purely greedy pathfinder that retains only a subset of candidates at each step, discarding all alternative branches. While this approach reduces memory overhead and execution time, the absence of a backtracking mechanism makes it vulnerable to local minima: a single misleading evaluation by the heuristic can permanently trap the search in a dead-end. Formally, for a K-node quiver, we fix a “beam width” B ≤ K and define the “beam” B t as the set of candidate quivers retained in memory at depth t, with |B t | ≤ B. At each exploration step, the algorithm generates the pool of all valid Seiberg mutants originating from the current beam: C t+1 = [ Q∈B t D k Q| D k Q̸=∅ .(C.3) 66 Algorithm 2 Bidirectional A ∗ Search for Seiberg Dualities 1: Input: Initial quiver Q A , target quiver Q B 2: Output: Shortest path between Q A and Q B , or Failure 3: Initialize Q fwd and Q bwd as empty priority queues 4: V fwd ←Q A : 0, V bwd ←Q B : 0▷ Track states and accumulated costs 5: Push (g fwd (Q A ) + h fwd (Q A ), Q A ) to Q fwd 6: Push (g bwd (Q B ) + h bwd (Q B ), Q B ) to Q bwd 7: while Q fwd is not empty and Q bwd is not empty do 8: if min(f fwd )≤ min(f bwd ) then 9:Pop Q n from Q fwd 10:for k ← 1 to K do 11:Q n ′ ← D k Q n 12:if Q n ′ ̸=∅ then 13:if ∃ Q m ∈V bwd s.t. Q n ′ ∼ = Q m then 14:return ReconstructPath(V fwd ,V bwd , Q n ′ , Q m ) ▷ Intersection found up to isomorphism 15:end if 16:g new ←V fwd [Q n ] + step cost 17:if Q n ′ /∈V fwd or g new <V fwd [Q n ′ ] then 18:V fwd [Q n ′ ]← g new 19:Push (g new + h fwd (Q n ′ ), Q n ′ ) to Q fwd 20:end if 21:end if 22:end for 23: else 24:// Perform expansion for the backward frontier (Q bwd ) 25: end if 26: end while 27: return Failure 67 Rather than maintaining an exhaustive frontier, the new beam B t+1 is formed by selecting the B quivers from C t+1 that minimize a guiding cost function g(Q): B t+1 =Q∈C t+1 ||Q ′ ∈C t+1 \Q| g(Q ′ ) < g(Q)| < B ,(C.4) where ties at the boundary |B t+1 | = B are broken deterministically by graph lexicographical ordering. Morally, the Beam Search is therefore an A ∗ Search with f (n) = g(n) and h(n) = 0, and a limited frontier. In Section 5.2, g(·) is estimated using the output of the AGNN model. During expansion, every generated candidate Q ′ is checked against the target Q B . A schematic representation of the beam search is shown in Figure 13b, with the corresponding pseudocode provided in Algorithm 3. Algorithm 3 Beam Search Algorithm 1: Input: Initial quiver Q A , target quiver Q B , beam width B, max depth d max 2: Output: Path to target state Q B , or Failure 3: B 0 ←Q A 4: V ←Q A 5: for t← 0 to d max − 1 do 6: C t+1 ←∅▷ Candidate set 7: for Q∈B t do 8:for k ← 1 to K do 9:Q ′ ← D k Q 10:if Q ′ ∼ = Q B then 11:return Path(Q A →·→ Q ′ )▷ Success up to isomorphism 12:end if 13:if Q ′ ̸=∅ and Q ′ /∈V then 14:Compute cost g(Q ′ )← g(Q) + cost function(Q ′ ) 15:C t+1 ←C t+1 ∪(Q ′ , g(Q ′ )) 16:end if 17:end for 18: end for 19: if C t+1 =∅ then 20:return Failure▷ dead-end 21: end if 22: B t+1 ← arg min-B s∈C t+1 g(s)▷ Retain top B candidates 23: V ←V ∪B t+1 24: end for 25: return Failure 68 References [1] N. Seiberg, “Electric - Magnetic Duality in Supersymmetric Non-Abelian Gauge Theories,” Nucl. Phys. B 435 (1995) 129–146, arXiv:hep-th/9411149. [2] S. Fomin and A. Zelevinsky, “Cluster Algebras I: Foundations,” arXiv:math/0104151. [3] S. Fomin and A. Zelevinsky, “Cluster Algebras I: Finite Type Classification,” Invent. Math. 154 no. 1, (2003) 63–121, arXiv:math/0208229. [4] I. R. Klebanov and M. J. Strassler, “Supergravity and a Confining Gauge Theory: Duality Cascades and χSB Resolution of Naked Singularities,” JHEP 08 (2000) 052, arXiv:hep-th/0007191. [5] I. R. Klebanov and E. Witten, “Superconformal Field Theory on Three-Branes at a Calabi-Yau Singularity,” Nucl. Phys. B 536 (1998) 199–218, arXiv:hep-th/9807080. [6] I. R. Klebanov and A. A. Tseytlin, “Gravity Duals of Supersymmetric SU (N )× SU (N + M ) Gauge Theories,” Nucl. Phys. B 578 (2000) 123–138, arXiv:hep-th/0002159. [7] S. Franco, Y.-H. He, C. Herzog, and J. Walcher, “Chaotic Duality in String Theory,” Phys. Rev. D 70 (2004) 046006, arXiv:hep-th/0402120. [8] J. J. Heckman, “Statistical Inference and String Theory,” Int. J. Mod. Phys. A 30 no. 26, (2015) 1550160, arXiv:1305.3621 [hep-th]. [9] V. Balasubramanian, J. J. Heckman, and A. Maloney, “Relative Entropy and Proximity of Quantum Field Theories,” JHEP 05 (2015) 104, arXiv:1410.6809 [hep-th]. [10] R. Fowler and J. J. Heckman, “Misanthropic Entropy and Renormalization as a Communication Channel,” Int. J. Mod. Phys. A 37 no. 16, (2022) 2250109, arXiv:2108.02772 [hep-th]. [11] M. Del Zotto, J. J. Heckman, S. N. Meynet, R. Moscrop, and H. Y. Zhang, “Higher Symmetries of 5D Orbifold SCFTs,” Phys. Rev. D 106 no. 4, (2022) 046010, arXiv:2201.08372 [hep-th]. [12] V. Chakrabhavi, M. Cvetiˇc, J. J. Heckman, and S. Meynet, “Quiver Approach to Symmetry Theories,” arXiv:2605.30354 [hep-th]. [13] J. Bao, S. Franco, Y.-H. He, E. Hirst, G. Musiker, and Y. Xiao, “Quiver Mutations, Seiberg Duality and Machine Learning,” Phys. Rev. D 102 no. 8, (2020) 086013, arXiv:2006.10783 [hep-th]. 69 [14] F. Carta, A. Gauntlett, F. Griffin, and Y.-H. He, “BPS spectroscopy with reinforcement learning,” Phys. Lett. B 868 (2025) 139646, arXiv:2501.14863 [hep-th]. [15] S. Gukov, J. Halverson, F. Ruehle, and P. Su lkowski, “Learning to Unknot,” Mach. Learn. Sci. Tech. 2 no. 2, (2021) 025035, arXiv:2010.16263 [math.GT]. [16] A. Gromov, “Grokking modular arithmetic,” arXiv:2301.02679 [cs.LG]. [17] N. Arkani-Hamed, J. L. Bourjaily, F. Cachazo, A. B. Goncharov, A. Postnikov, and J. Trnka, Grassmannian Geometry of Scattering Amplitudes. Cambridge University Press, 4, 2016. arXiv:1212.5605 [hep-th]. [18] Y.-H. He, “Lectures on D-branes, gauge theories and Calabi-Yau singularities,” in 1st Hangzhou-Beijing International Summer School. 8, 2004. arXiv:hep-th/0408142. [19] P. S. Aspinwall, “D-branes on Calabi-Yau manifolds,” in Theoretical Advanced Study Institute in Elementary Particle Physics (TASI 2003): Recent Trends in String Theory, p. 1–152. 3, 2004. arXiv:hep-th/0403166. [20] S. Franco, Y.-H. He, C. Sun, and Y. Xiao, “A Comprehensive Survey of Brane Tilings,” Int. J. Mod. Phys. A 32 no. 23n24, (2017) 1750142, arXiv:1702.03958 [hep-th]. [21] A. Vaswani, N. Shazeer, N. Parmar, J. Uszkoreit, L. Jones, A. N. Gomez, L. Kaiser, and I. Polosukhin, “Attention Is All You Need,” arXiv e-prints (June, 2017) arXiv:1706.03762, arXiv:1706.03762 [cs.CL]. [22] P. E. Hart, N. J. Nilsson, and B. Raphael, “A Formal Basis for the Heuristic Determination of Minimum Cost Paths,” IEEE Transactions on Systems Science and Cybernetics 4 no. 2, (1968) 100–107. [23] B. T. Lowerre, The Harpy speech recognition system. PhD thesis, Carnegie Mellon University, Pennsylvania, Apr., 1976. [24] J. Wang, H. Liu, S. Di, Z. Wang, J. Wang, L. Chen, and X. Zhou, “Proficient Graph Neural Network Design by Accumulating Knowledge on Large Language Models,” arXiv e-prints (Aug., 2024) arXiv:2408.06717, arXiv:2408.06717 [stat.ML]. [25] X. Yu, “DualityCert: Verifier-Gated Language-Model Repair of Broken Duality Claims in Quantum Field Theory,” arXiv:2607.23614 [cs.CR]. [26] I. Affleck, M. Dine, and N. Seiberg, “Dynamical Supersymmetry Breaking in Supersymmetric QCD,” Nucl. Phys. B 241 (1984) 493–534. [27] K. A. Intriligator and N. Seiberg, “The Runaway quiver,” JHEP 02 (2006) 031, arXiv:hep-th/0512347. 70 [28] D. Berenstein and M. R. Douglas, “Seiberg Duality for Quiver Gauge Theories,” arXiv:hep-th/0207027. [29] C. P. Herzog, “Exceptional Collections and del Pezzo Gauge Theories,” JHEP 04 (2004) 069, arXiv:hep-th/0310262. [30] C. P. Herzog, “Seiberg Duality is an Exceptional Mutation,” JHEP 08 (2004) 064, arXiv:hep-th/0405118. [31] P. S. Aspinwall and I. V. Melnikov, “D-Branes on Vanishing del Pezzo Surfaces,” JHEP 12 (2004) 042, arXiv:hep-th/0405134. [32] S. Franco, A. Hanany, Y.-H. He, and P. Kazakopoulos, “Duality walls, duality trees and fractional branes,” arXiv:hep-th/0306092. [33] J. Gilmer, S. S. Schoenholz, P. F. Riley, O. Vinyals, and G. E. Dahl, “Neural Message Passing for Quantum Chemistry,” arXiv e-prints (Apr., 2017) arXiv:1704.01212, arXiv:1704.01212 [cs.LG]. [34] L. Ramp ́aˇsek, M. Galkin, V. P. Dwivedi, A. T. Luu, G. Wolf, and D. Beaini, “Recipe for a General, Powerful, Scalable Graph Transformer,” arXiv e-prints (May, 2022) arXiv:2205.12454, arXiv:2205.12454 [cs.LG]. [35] A. Vaswani, N. Shazeer, N. Parmar, J. Uszkoreit, L. Jones, A. N. Gomez, L. Kaiser, and I. Polosukhin, “Attention Is All You Need,” in Advances in Neural Information Processing Systems, vol. 30, p. 5998–6008. Curran Associates, Inc., 2017. [36] J. J. Heckman and C. Vafa, “Crystal Melting and Black Holes,” JHEP 09 (2007) 011, arXiv:hep-th/0610005. [37] M. Alim, S. Cecotti, C. Cordova, S. Espahbodi, A. Rastogi, and C. Vafa, “BPS Quivers and Spectra of Complete N=2 Quantum Field Theories,” Commun. Math. Phys. 323 (2013) 1185–1227, arXiv:1109.4941 [hep-th]. [38] J. H. T. Yip, C. Arnal, F. Charton, and G. Shiu, “Transforming Calabi-Yau Constructions: Generating New Calabi-Yau Manifolds with Transformers,” arXiv:2507.03732 [hep-th]. [39] M. Walden and M. Larfors, “Sampling string vacua using generative models,” Mach. Learn. Sci. Tech. 7 no. 1, (2026) 015018, arXiv:2509.16029 [hep-th]. [40] C. Arnal, J. H. T. Yip, F. Charton, and G. Shiu, “Generating Special Triangulations with Transformers,” in . 6, 2026. arXiv:2606.26660 [hep-th]. [41] J. H. T. Yip, A. Mininno, and G. Shiu, “Exploring Line Bundle Standard Models with Transformers,” arXiv:2607.00078 [hep-th]. 71 [42] H. Cao, G. Merz, K. Cranmer, and G. Shiu, “Reconstructing conformal field theoretical compositions with Transformers,” arXiv:2605.01072 [hep-th]. [43] A. R. Brown, D. A. Roberts, L. Susskind, B. Swingle, and Y. Zhao, “Holographic Complexity Equals Bulk Action?,” Phys. Rev. Lett. 116 no. 19, (2016) 191301, arXiv:1509.07876 [hep-th]. [44] D. Carmi, R. C. Myers, and P. Rath, “Comments on Holographic Complexity,” JHEP 03 (2017) 118, arXiv:1612.00433 [hep-th]. [45] V. Balasubramanian, M. Decross, A. Kar, and O. Parrikar, “Quantum Complexity of Time Evolution with Chaotic Hamiltonians,” JHEP 01 (2020) 134, arXiv:1905.05765 [hep-th]. [46] J. J. Heckman, R. J. Hicks, and C. Murdia, “Generalized Complexity Distances and Non-Invertible Symmetries,” arXiv:2604.14275 [hep-th]. [47] B. A. Burrington, J. T. Liu, and L. A. Pando Zayas, “Finite Heisenberg groups from nonAbelian orbifold quiver gauge theories,” Nucl. Phys. B 794 (2008) 324–347, arXiv:hep-th/0701028. [48] G. Beaujard, J. Manschot, and B. Pioline, “Vafa–Witten Invariants from Exceptional Collections,” Commun. Math. Phys. 385 no. 1, (2021) 101–226, arXiv:2004.14466 [hep-th]. [49] J. Lei Ba, J. R. Kiros, and G. E. Hinton, “Layer Normalization,” arXiv e-prints (July, 2016) arXiv:1607.06450, arXiv:1607.06450 [stat.ML]. [50] A. L. Maas, “Rectifier Nonlinearities Improve Neural Network Acoustic Models,” in . 2013. [51] X. Wu, A. Ajorlou, Z. Wu, and A. Jadbabaie, “Demystifying Oversmoothing in Attention-Based Graph Neural Networks,” arXiv e-prints (May, 2023) arXiv:2305.16102, arXiv:2305.16102 [cs.LG]. 72