A Network Topology Dependent Upper Bound on the Number of Equilibria
of the Kuramoto Model

arXiv:1603.05905v1 [math.OC] 18 Mar 2016

Tianran Chen1, Dhagash Mehta2 and Matthew Niemerg3
Abstract— We begin with formulating the stationary equations of the Kuramoto model as a system of
polynomial equations in a novel way. Then, based on
an algebraic geometric root count, we give a prescription of computing an upper bound on the number of
equilibria of the Kuramoto model for the most general
case, i.e., defined on an arbitrary graph and having
generic values of natural frequencies and inhomogeneous couplings. We demonstrate with computational
experiments utilizing the numerical polynomial homotopy continuation method that our bound is tight for
the number of complex equilibria for the Kuramoto
model in the most general case. We then discuss the
implications or our results in relation to finding all
the real equilibria of the Kuramoto model.

I. Introduction
The Kuramoto model [1]–[4] is a prototypical model for
studying many phenomena including the synchronization
power systems, neural networks, chemical oscillators,
particle coordinations, rhythmic applause, and so on.
The general Kuramoto model is given by the diﬀerential
equation:
N
1 X
dθi
Ki,j sin(θi − θj ), for i = 1, ..., N, (1)
= ωi −
dt
N j=1

where N is the number of oscillators, Ki,j is the coupling strength between the i-th and j-th oscillators.
The matrix K = [Ki,j ] may also be viewed as the
adjacency matrix for the underlying weighted graph.
Ω = (ω1 , . . . , ωN ) contains the natural frequencies of
the N oscillators. Of crucial importance in studying
the phase space of this system of ordinary diﬀerential
equations are the equilibria which are values of θ1 , . . . , θN
i
for which dθ
dt = 0 for all i = 1, . . . , N . The stability
analysis of the equilibria can reveal the behavior of the
dynamical system near the equilibria. One can also study
synchronization phenomenon of the Kuramoto model
with the knowledge of the equilibria of the model [5].
Interpreting as the special case [6] of the power ﬂow
equations the equilibria of the Kuramoto model provide
1 Department of Mathematics, Michigan State University, East
Lansing MI, USA. chentia1@msu.edu
2 D. Mehta is with Department of Applied and Computational Mathematics and Statistics, University of Notre Dame,
IN, USA, and was supported by NSF-ECCS award ID 1509036.

dmehta@nd.edu
3 IBM T.J. Watson Research Center and was supported in part
by the Fields Institute for Research in Mathematical Sciences.
mniemer@us.ibm.com.

crucial information for planning and designing power
grids.
Since the equilibrium conditions are invariant under
transformations of the form θi → θi +α for all θi with any
ﬁxed α ∈ (−π, π], they possess inﬁnitely many solutions.
To remove this degree of freedom, we ﬁx θN = 0 and
remove the N -th equation dθdtN = 0. Thus, we are left
with N − 1 nonlinear equations in N − 1 angles.
The problem of counting the number of equilibria of
the ﬁnite N Kuramoto model has been tackled for diﬀerent graphs such as complete graphs [7], nearest-neighbor
coupling on one-dimensional lattice graphs [7]–[11] and
with long-range interactions [12], two-dimensional lattice graphs [8], [13]–[17], and three-dimensional lattice
graphs [17]–[19]; and the homogeneous frequencies case
(also known as the power ﬂow equations for the lossless
power system) [20], in diﬀerent disguises.
For N = 3, the Kuramoto model on a complete
graph of 3 nodes was shown to have at most 6 equilibria
[21]–[23]. Several necessary and suﬃcient conditions for
equilibria to exist and other conjectures on the types of
equilibria for diﬀerent graphs are provided in [24]–[32]
though some of which were disproved by counterexamples in [5].
It was also shown in Ref. [7] that even after ﬁxing
the trivial zero mode by setting θN = 0, there existed
inﬁnitely many equilibria (also called incoherent manifolds [33]). Surprisingly, such inﬁnite families of equilibria
were also shown to exist in the one-, two- and threedimensional lattice model with nearest neighbour interaction in [8], [15], [34], where it was demonstrated that
the number of inﬁnite families of equilibria exponentially
grows with N .
As a signiﬁcant leap forward from the recent studies
of the Kuramoto model and the closely related load
ﬂow equations (for power networks) from algebraic view
points by the authors [5], [35], in the present contribution, we establish an upper bound on the number of equilibria for (1) using a novel polynomial formulation of the
and the theory of Bernshtein-Kushnirenko-Khovanskii
(BKK) bound [36]–[38]. The advantage of this bound
over existing bounds is that it takes into consideration
of the sparsity of the connections in the underlying
network. We present strong numerical evidence that for
many Kuramoto models, the proposed bound is best
possible bound for the number of complex equilibria. In
§ II-A, we describe a novel polynomial formulation of
the equilibrium equations of the Kuramoto model. After

reviewing existing bounds on the number of equilibria in
§ II-B, we propose an upper bound on the number of
equilibria based on the theory of BKK bound in § II-C.
We compute this bound for the Kuramoto models deﬁned
on diﬀerent graphs. Employing the numerical polynomial
homotopy continuation method, we provide a strong
numerical evidence to demonstrate that this bound to be
tight in these cases. In IV, we discuss the implications of
our results and conclude.
II. Upper Bound on the Number of Equilibria
In this section we discuss a novel polynomial formulation of the Kuramoto model and provide an upper bound
based on the theory of BKK bound.
A. Polynomial Formulations
In studying nonlinear systems of equations like the
equilibrium equations of (1), it is a common practice
to ﬁrst transform them into algebraic equations which
would allow the use of powerful tools from algebraic
geometry. Previously (e.g. , [5], [8], [16], [20], [39]), the
equilibrium conditions of Eqs. (1) were transformed into
a system of polynomial equations by ﬁrst using the
identities sin(θi − θj ) = sin θi cos θj − sin θj cos θi and
then the substitution si := sin θi and ci := cos θi for
all i = 1, . . . , N − 1, as well as adding the constraint
equations s2i +c2i −1 = 0, for all i = 1, . . . , N −1, to ﬁnally
obtain the following system of polynomial equations:
ωi −

N
1 X
Ki,j (si cj − sj ci ) = 0
N j=1

(2)

s2i + c2i − 1 = 0,
for i = 1, . . . , N − 1.
In the present work, we adopt a diﬀerent transformation: using the trigonometric identity,
sin(θi − θj ) =
√
1
I(θi −θj )
−I(θi −θj )
(e
−e
)
where
I
:=
−1
is
the imaginary
2I
unit, the equilibrium equations of (1) becomes
ωi −

N
X
Ki,j
j=1

IN

(eIθi e−Iθj − e−Iθi eIθj ) = 0.

(3)

To formulate the above system of equations as an algebraic system, we let
xi := eIθi

and

yi := e−Iθi ,

(4)

for all i = 1, . . . , N − 1. With this substitution (3)
becomes an enlarged system of 2(N − 1) equations in
2(N − 1) variables
N
X
Ki,j
j=1

IN

(xi yj − xj yi ) = ωi
xi yi = 1

for i = 1, . . . , N − 1
for i = 1, . . . , N − 1.

(5)

For example, with N = 3, the system becomes
K1,2
(x1 y2 − x2 y1 ) +
3I
K2,1
(x2 y1 − x1 y2 ) +
3I

K1,3
(x1 − y1 ) − ω1 = 0
3I
K2,3
(x2 − y2 ) − ω2 = 0
3I
x1 y1 − 1 = 0

(6)

x2 y2 − 1 = 0.

It can be readily veriﬁed that the equilibria of (1) (with
the translation symmetry removed) are in one-to-one
correspondence with the special solutions of the above
system (5) that satisfy the additional restriction that
|xi | = |yi | = 1 for i = 1, . . . , N − 1.
B. Existing Bounds
Via the transformations given in (2) or (5), the problem of counting equilibria of the Kuramoto model is
turned into a root counting problem for systems of
polynomial equations. A well known bound on the root
count of a system of polynomial equations is the Bézout
bound which is simply the product of the degrees of all
the equations. In the example shown in (6), since each of
the four equations is quadratic (degree 2), the highest
possible number of isolated solutions as given by the
Bézout bound is therefore 24 = 16. In general, the Bézout
bound for the system (5) is 22(N −1) .
The Bézout bound is a basic result in intersection
theory [40]. Based on this theory, a more reﬁned root
count can be derived for the equilibrium equations for
(1). This was ﬁrst studied in Ref. [20] in which
it was

−1)
shown to be the binomial coeﬃcient 2(N
,
which
will
N −1
be referred to as the binomial bound.1
C. BKK Bound
In this subsection, we brieﬂy review the Bernstein–
Khovanskii–Kushnirenko (BKK) bound [36]–[38] for the
number of isolated nonzero complex solutions which is
a reﬁnement of the Bézout bound that takes into consideration the monomials that appear in the polynomial
system: Given a polynomial, each of its terms give rise to
an exponent vector. For instance, for the term x3 y 2 z 1 , the
exponent vector is simply the vector whose entries are
the exponents of x, y and z, respectively, i.e., (3, 2, 1).
The choice of this ordering is inconsequential as long
as it is kept the same for each equation. The set of
all exponent vectors derived from the nonzero terms of
an polynomial equation is called the support of that
equation. For example, if we arrange the variables in
the order of (x1 , y1 , x2 , y2 ), then the supports of the four
equations in (6) are
S1
S2
S3
S4

= {(1, 0, 0, 1), (0, 1, 1, 0), (1, 0, 0, 0), (0, 1, 0, 0), (0, 0, 0, 0)}
= {(0, 1, 1, 0), (1, 0, 0, 1), (0, 0, 1, 0), (0, 0, 0, 1), (0, 0, 0, 0)}
= {(1, 1, 0, 0), (0, 0, 0, 0)}
= {(0, 0, 1, 1), (0, 0, 0, 0)}.

1 In Ref. [20], the system (1) was termed as a special case of
the power ﬂow equations. Later on, the same binomial bound was
proposed for the complete power ﬂow equations [41] (see [42] for an
alternative proof for this bound).

A convex set is a set of points in which the line segment
connecting any pair of points in the set also lie in that
set. The convex hull of a set is the minimal convex set
containing that set. For a polynomial, the convex hull
of its support is known as the Newton polytope of that
polynomial. In the study of convex polytopes, the mixed
volume of several polytopes is an important concept. In
the simplest case, the mixed volume of two line segments
on the plane is precisely the area of the parallelogram
spanned by translations of these two line segments. The
general deﬁnition is included in the Appendix. Bernstein’s Theorem [36] provides an upper bound for the
number solutions to a system of polynomial equations in
terms of mixed volume:
Theorem 1 (Bernshtein [36]): Given a system of n
polynomial equations in n variables, the number of isolated complex solutions for which no variable is zero
is bounded above by the mixed volume of the Newton
polytopes of the equations.
This upper bound is known as the BKK bound, after a circle of closely related works by Bernstein[36],
Kushnirenko[37], and Khovanskii[38]. Recall that via the
change of variables given in (4), each (real) equilibrium of
(1) corresponds to a unique nonzero complex solution of
(5). Therefore the BKK bound given above provides an
upper bound to the number of isolated (real) equilibria.
A crucial fact is that the BKK bound is generically
exact: when the coeﬃcients in the polynomial system
are chosen at random, with probability one the number of
isolated complex solutions for which no variable is zero is
exactly the BKK bound. In the polynomial formulation
of the Kuramoto model given in (5), if certain relations
are imposed on the coeﬃcients (e.g., the coeﬃcients of
x1 y2 and x2 y1 in (6) must be the same) the generic exactness still holds true with a mild additional condition:
Theorem 2: If there exists a choice of Ki,j ’s and ωi ’s
for which the number of nonzero complex solutions of (5)
is the BKK bound, then for almost all choices of complex
Ki,j ’s and ωi , the number of nonzero complex solutions
of (5) will be the BKK bound.
In other words, among the systems (5) for all possible
choices of the Ki,j ’s and ωi ’s, if the BKK bound is
attainable then it must also be generically exact. In
the following, we shall compute the BKK bound for the
polynomial system (5) induced by a number of graphs.
Then, the attainability and hence the generic exactness
of the BKK bound in each case is veriﬁed by solving the
system (5) for some speciﬁc chosen set of Ki,j ’s and ωi ’s.
III. Results
In this section, we provide results for the Kuramoto
model on diﬀerent graphs. We compare the BKK bound
for both formulations (2) and (5) with the help of
the computational packages HOM4PS-3.0 [43], [44] and
Bertini [45], [46]. We also compare the bounds with the
corresponding binomial bound.

To analyze how tight the above mentioned bounds
are we provide results from our numerical experiments
using the NPHC method. The NPHC method guarantees
that one will obtain all isolated complex solutions for a
system of polynomial equations by following the strategy
as brieﬂy explained below[47], [48]: to solve a system of
polynomial equations, one starts with an upper bound
on the number of complex solutions of the system. Then,
another system is created such that the system has exactly the same number of complex solutions as the upper
bound, and it is easy to solve. Finally, each solution
of this new system is evolved over a single parameter
towards the system to be solved. The so-called gammatrick is employed to ensure to achive all the isolated
complex solutions of the system.
A. BKK for the Complete Graphs
The complete graph is a graph in which each node
is connected to all other nodes. Hence, this is the most
densely connected graph. We compare the Bézout bound,
binomial bound, the BKK bound for the traditional
polynomial formulation (2), the BKK bound for the
new formulation (5), and the actual number of complex
solutions (with 100 random samples for each graph) in
Table I.
B. BKK for the Path Graphs
Now we turn our attention to one of the most sparsely
connected graphs — path graph. In a path graph, the
i-th node for 1 < i < N − 1 is connected to two of its
neighbors: the (i − 1)-th node and the (i + 1)-th node
forming a path. The results are given in Table II.
C. BKK for the Ring Graphs
A ring graph is a graph in which each node is connected
only to its nearest neighbors, same as the path graph,
and the last node is connected to the ﬁrst node making
a closed loop. The results are in Table III.
D. BKK for Erdős-Rényi Random Graphs
The Erdős-Rényi model for random graphs was initially studied in [49]–[51]. There are two formulations of
the model, but we will use the convention given originally
by Gilbert, that is the existence of each edge is i.i.d. with
probability p. Such graphs are denoted as G(N, p) where
N is the number of vertices of the graph and p is the
probability that any edge exists.
We construct a random graph G(N, p) by initially
creating an N × N empty adjacency matrix, A. Then,
for each Ai,j where i > j, we test for the existence of
edge Ei,j with the given probability p and update Ai,j
and Aj,i accordingly.
Connectedness of a graph is necessary to have nontrivial equilibria for the Kuramoto model. Since the
probability that a random graph G(N, p) is connected
increases if p > pc = ln(N )/N , we let p vary between pc
and 1. Testing for connectedness is done by computing

No. of Nodes
Bézout
Binomial
BKK for (2)
BKK for (5)
Generic root count

3
16
6
8
6
6

4
64
20
40
20
20

5
256
70
192
70
70

6
1024
252
864
252
252

7
4096
924
3712
924
924

8
16384
3432
15488
3432
3432

9
65536
12870
63488
12870
12870

10
262144
48620
257536
48620
48620

11
1048576
184756
1038336
184756
184756

12
4194304
705432
4171776
705432
705432

13
16777216
2704156
16728064
2704156
2704156

14
67108864
10400600
67002368
10400600
10400600

15
268435456
40116600
268206080
40116600
40116600

TABLE I
Comparision among different bounds on the number of equilibria and the actual number of complex solutions for
generic parameter-values for the Kuramoto model on the complete graph.

Nodes
Bezout
Binomial
BKK for (2)
BKK for (5)
Generic root count

3
16
6
8
4
4

4
64
20
24
8
8

5
256
70
80
16
16

6
1024
252
256
32
32

7
4096
924
832
64
64

8
16384
3432
2688
128
128

9
65536
12870
8704
256
256

10
262144
48620
28160
512
512

11
1048576
184756
91136
1024
1024

12
4194304
705432
294912
2048
2048

13
16777216
2704156
954368
4096
4096

14
67108864
10400600
3088384
8192
8192

15
268435456
40116600
9994240
16384
16384

16
1073741824
155117520
32342016
32768
32768

TABLE II
Comparison among different bounds on the number of equilibria and the actual number of complex solutions for
generic parameter-values for the Kuramoto model on the path graph.

Nodes
Bézout
Binomial
BKK for (2)
BKK for (5)
Generic root count

3
16
6
8
6
6

4
64
20
24
12
12

5
256
70
80
30
30

6
1024
252
256
60
60

7
4096
924
832
140
140

8
16384
3432
2688
280
280

9
65536
12870
8704
630
630

10
262144
48620
28160
1260
1260

11
1048576
184756
91136
2772
2772

12
4194304
705432
294912
5544
5544

13
16777216
2704156
954368
12012
12012

14
67108864
10400600
3088384
24024
24024

15
268435456
40116600
9994240
51480
51480

16
1073741824
155117520
32342016
102960
102960

TABLE III
Comparision among different bounds on the number of equilibria and the actual number of complex solutions for
generic parameter-values for the Kuramoto model on the ring graph.

the size of the largest component of G(N, p) and testing
for equality with N .
For both formulations of the Kuramoto model, the
Bézout bound is always 22(N −1) . For all the random
Erdős-Rényi connected graphs with N = 5 and 6 that
we tested, we found that the complex solution count of
the associated polynomial system under the Kuramoto
model using the formulation of (5) was the same as
the BKK bound. For N = 5 and 6, we calculated the
number of complex roots for 10 random instances of the
Kuramoto model using (2) and (5). Tables IV and V
compare the BKK root count of polynomial systems arising from the Kuramoto model associated to 10 instances
of random Erdős-Rényi connected graphs using the two
formulations.
IV. Discussion and Conclusion
In this article, we have translated the stationary equations of the Kuramoto model to the polynomial form in
a diﬀerent way than the traditional polynomial formulation to utilize complex algebraic geometry results. We
have considered a general form of the Kuramoto model,
i. e. , with arbitrary natural frequencies and on graphs
with arbitrary weights. We then provide a prescription to
compute a tight upper bound, called the BKK bound, on
the number of equilibria of the model for a given graph
topology, i.e., the bound captures the graph topology
accurately.
We ﬁrst reproduced
the previously available upper

−1)
bound, 2(N
,
by
computing
the BKK bound for
N −1
the complete graph with arbitrary (and inhomogeneous)

coupling strengths and natural frequencies: the binomial
bound is the tightest possible bound for the Kuramoto
model on the complete graph. The previous upper bound,
however, does not take into account the sparsity of the
graphs unlike ours. For other sparser graphs such as path
graphs, ring graphs and Erdős-Rényi random graphs,
the BKK bound for the new polynomial formulation is
signiﬁcantly lower than the binomial bound as well as the
BKK bound for the traditional polynomial formulation.
Furthermore, with extensive numerical experiments
with the NPHC method, which guarantees to ﬁnd all
isolated complex solutions of systems of polynomial equations, we ﬁnd that for generic natural frequencies and
weights the number of complex solutions is always equal
to the BKK bound of the new polynomial formulation.
Hence, this is the best possible upper bound on the
number of complex solutions for generic values of natural
frequencies and inhomogeneous couplings.
Our prescription is also constructive in solving the
stationary equations of the model in that we not only
provide the tight-most upper bound but also the optimum way of constructing a starting system for the NPHC
method. Because our numerical experiments suggest that
the BKK bound is always equal to the number of complex
solutions for generic parameter values, no path diverges
in the NPHC computation.
This tight upper bound not only helps the NPHC
method to solve the correpsonding systems more eﬃciently, but it also provides a concrete stopping criterion
to other stochastic nonlinear equation solving method.
Our results also give rise to the following interesting

Sample no.
Bézout bound
Binomial
BKK for (2)
BKK for (5)
Complex Count

1
256
70
48
24
24

2
256
70
144
32
32

3
256
70
112
24
24

4
256
70
80
24
24

5
256
70
16
16
16

6
256
70
112
24
24

7
256
70
80
16
16

8
256
70
48
24
24

9
256
70
80
32
32

10
256
70
48
16
16

TABLE IV
Comparision among different bounds on the number of equilibria and the actual number of complex solutions for
generic parameter-values for the Kuramoto model on the random graph of 5 nodes for 10 samples.

Sample no.
Bézout bound
Binomial
BKK for (2)
BKK for (5)
Complex Count

1
1024
252
128
48
48

2
1024
252
320
48
48

3
1024
252
288
48
48

4
1024
252
544
88
88

5
1024
252
480
104
104

6
1024
252
544
144
144

7
1024
252
224
48
48

8
1024
252
160
64
64

9
1024
252
160
32
32

10
1024
252
608
140
140

TABLE V
Comparision among different bounds on the number of equilibria and the actual number of complex solutions for
generic parameter-values for the Kuramoto model on the random graph of 6 nodes for 10 samples.

and mathematically challenging questions: can there be
some values of Ki,j ’s and ωi ’s so that the number of
real equilibria attains the BKK bound for the given
graph? If not, what is the maximum number of real
solutions a particular graph can have? Though there are
a few algebraic geometry based methods[52] have been
developed to answer such questions, they can not yet deal
with systems beyond simple ones. Moreover, it will also
be interesting to investigate how the number of complex
and real equilibria behave for non-generic values of the
couplings and frequencies. We plan to study all these
issues in future.
In the future, we plan to explore the advantages and
disadvantages of both polynomial formulations in terms
of the numerical stability, ease to implement, and so on.
Systems with sine and cosine of angles and diﬀerence of
angles are not unique to the Kuramoto model, and hence
such a comparative analysis can have far reaching implications. We also point out a remarkable parallel between
upper bounds for the equilibria of the Kuramoto model
and those for the equilibria of the complete power ﬂow
equations [53]–[55]. We plan to quantitatively establish
this similarity in the future.
Appendix
Given n convex polytopes Q1 , . . . , Qn ⊂ Rn and positive real numbers λ1 , . . . , λn Minkowski’s Theorem [48]
states that the n-dimensional volume of the Minkowski
sum λ1 Q1 + · · · + λn Qn , deﬁned as
{ λ1 q1 + · · · + λn qn | qi ∈ Qi for i = 1, . . . , n }
is a homogeneous polynomial of degree n in the variables
λ1 , . . . , λn . The coeﬃcient associated with the monomial
λ1 · · · λn in this polynomial is known as the mixed volume
of the polytopes Q1 , . . . , Qn .

References
[1] Yoshiki Kuramoto, “Self-entrainment of a population of coupled non-linear oscillators,” in International symposium on
mathematical problems in theoretical physics. Springer, 1975,
pp. 420–422.
[2] Juan A Acebrón, Luis L Bonilla, Conrad J Pérez Vicente, Félix
Ritort, and Renato Spigler, “The Kuramoto model: A simple
paradigm for synchronization phenomena,” Reviews of modern
physics, vol. 77, no. 1, pp. 137, 2005.
[3] Steven H Strogatz, “From Kuramoto to Crawford: exploring
the onset of synchronization in populations of coupled oscillators,” Physica D: Nonlinear Phenomena, vol. 143, no. 1, pp.
1–20, 2000.
[4] F. Dörﬂer and F. Bullo, “Synchronization in complex oscillator
networks: A survey,” Automatica, vol. 50, no. 6, pp. 1539–1564,
June 2014.
[5] Dhagash Mehta, Noah S. Daleo, Florian Dörﬂer, and
Jonathan D. Hauenstein, “Algebraic geometrization of the
kuramoto model: Equilibria and stability analysis,” Chaos,
vol. 25, no. 5, pp. –, 2015.
[6] Florian Dorﬂer and Francesco Bullo, “Synchronization and
transient stability in power networks and non-uniform Kuramoto oscillators,” SIAM Journal on Control and Optimization, vol. 50, no. 3, pp. 1616–1642, 2010.
[7] L. Casetti, M. Pettini, and E. G. D. Cohen, “Phase transitions
and topology changes in conﬁguration space,” Journal of
Statistical Physics, vol. 111, pp. 1091–1123(33), June 2003.
[8] Dhagash Mehta, “Lattice vs. Continuum: Landau Gauge
Fixing and ’t Hooft-Polyakov Monopoles,” Ph.D. Thesis, The
Uni. of Adelaide, Australasian Digital Theses Program, 2009.
[9] Dhagash Mehta and Michael Kastner, “Stationary point
analysis of the one-dimensional lattice Landau gauge ﬁxing
functional, aka random phase XY Hamiltonian,” Annals
Phys., vol. 326, pp. 1425–1440, 2011.
[10] Lorenz von Smekal, Dhagash Mehta, Andre Sternbeck, and
Anthony G. Williams, “Modiﬁed Lattice Landau Gauge,”
PoS, vol. LAT2007, pp. 382, 2007.
[11] Lorenz von Smekal, Alexander Jorkowski, Dhagash Mehta,
and Andre Sternbeck, “Lattice Landau gauge via Stereographic Projection,” PoS, vol. CONFINEMENT8, pp. 048,
2008.
[12] Michael Kastner, “Stationary-point approach to the phase
transition of the classical xy chain with power-law interactions,” Physical Review E, vol. 83, no. 3, pp. 031114, 2011.
[13] Lapo Casetti, Cesare Nardini, and Rachele Nerattini, “Microcanonical relation between continuous and discrete spin

[14]

[15]
[16]
[17]
[18]

[19]
[20]

[21]
[22]
[23]
[24]
[25]

[26]
[27]

[28]

[29]
[30]

[31]
[32]

[33]
[34]

models,” Physical review letters, vol. 106, no. 5, pp. 057208,
2011.
Cesare Nardini, Rachele Nerattini, and Lapo Casetti, “Density
of states of continuous and discrete spin models: a case study,”
Journal of Statistical Mechanics: Theory and Experiment, vol.
2012, no. 02, pp. P02007, 2012.
Rachele Nerattini, Michael Kastner, Dhagash Mehta, and
Lapo Casetti, “Exploring the energy landscape of XY models,”
Phys.Rev., vol. E87, no. 3, pp. 032140, 2013.
Dhagash Mehta, Andre Sternbeck, Lorenz von Smekal, and
Anthony G Williams, “Lattice Landau Gauge and Algebraic
Geometry,” PoS, vol. QCD-TNT09, pp. 025, 2009.
Ciaran Hughes, Dhagash Mehta, and Jon-Ivar Skullerud,
“Enumerating Gribov copies on the lattice,” Annals Phys.,
vol. 331, pp. 188–215, 2013.
Dhagash Mehta, Ciaran Hughes, Mario Schröck, and David J.
Wales, “Potential Energy Landscapes for the 2D XY Model:
Minima, Transition States and Pathways,” J. Chem. Phys.,
vol. 139, pp. 194503, 2013.
Dhagash Mehta and Mario Schröck, “Enumerating copies in
the ﬁrst gribov region on the lattice in up to four dimensions,”
2014.
John Baillieul and Christopher I Byrnes, “Geometric critical
point analysis of lossless power system models,” Circuits and
Systems, IEEE Transactions on, vol. 29, no. 11, pp. 724–737,
1982.
Carlos J Tavora and Otto JM Smith, “Equilibrium analysis
of power systems,” Power Apparatus and Systems, IEEE
Transactions on, , no. 3, pp. 1131–1137, 1972.
J Baillieul, “The critical point analysis of electric power
systems,” in The 23rd IEEE Conference on Decision and
Control, 1984, number 23, pp. 154–159.
A Klos and J Wojcicka, “Physical aspects of the nonuniqueness of load ﬂow solutions,” International Journal of Electrical
Power & Energy Systems, vol. 13, no. 5, pp. 268–276, 1991.
A. Araposthatis, S. Sastry, and P. Varaiya, “Analysis of powerﬂow equation,” International Journal of Electrical Power &
Energy Systems, vol. 3, no. 3, pp. 115–126, 1981.
D. Aeyels and J. A. Rogge, “Existence of partial entrainment
and stability of phase locking behavior of coupled oscillators,”
Progress on Theoretical Physics, vol. 112, no. 6, pp. 921–942,
2004.
R. E. Mirollo and S. H. Strogatz, “The spectrum of the locked
state for the Kuramoto model of coupled oscillators,” Physica
D: Nonlinear Phenomena, vol. 205, no. 1-4, pp. 249–266, 2005.
Mark Verwoerd and Oliver Mason, “On computing the critical
coupling coeﬃcient for the Kuramoto model on a complete bipartite graph,” SIAM Journal on Applied Dynamical Systems,
vol. 8, no. 1, pp. 417–453, 2009.
Mark Verwoerd and Oliver Mason, “Global phase-locking in
ﬁnite populations of phase-coupled oscillators,” SIAM Journal
on Applied Dynamical Systems, vol. 7, no. 1, pp. 134–160,
2008.
J Ochab and PF Góra, “Synchronization of coupled oscillators
in a local one-dimensional Kuramoto model,” arXiv preprint
arXiv:0909.0043, 2009.
F. Dörﬂer, M. Chertkov, and F. Bullo, “Synchronization in
complex oscillator networks and smart grids,” Proceedings of
the National Academy of Sciences, vol. 110, no. 6, pp. 2005–
2010, 2013.
Florian Dörﬂer and Francesco Bullo, “On the critical coupling
for kuramoto oscillators,” SIAM Journal on Applied Dynamical Systems, vol. 10, no. 3, pp. 1070–1099, 2011.
Richard Taylor, “There is no non-zero stable ﬁxed point
for dense networks in the homogeneous Kuramoto model,”
Journal of Physics A: Mathematical and Theoretical, vol. 45,
no. 5, pp. 055102, 2012.
Steven H Strogatz and Renato E Mirollo, “Stability of incoherence in a population of coupled oscillators,” Journal of
Statistical Physics, vol. 63, no. 3-4, pp. 613–635, 1991.
Dhagash Mehta and Michael Kastner, “Stationary point
analysis of the one-dimensional lattice landau gauge ﬁxing
functional, aka random phase xy hamiltonian,” Annals of
Physics, vol. In Press, pp. –, 2011.

[35] Tianran Chen and Dhagash Mehta, “On the Network Topology Dependent Solution Count of the Algebraic Load Flow
Equations,” arXiv:1512.04987 [cs, math], Dec. 2015.
[36] D. N. Bernstein, “The number of roots of a system of
equations,” Funkts. Anal. Pril., vol. 9, pp. 1–4, 1975.
[37] A. G. Kushnirenko, “Newton polytopes and the bezout
theorem,” Funkts. Anal. Pril., vol. 10, no. 3, 1976.
[38] A. G. Khovanski, “Newton polyhedra and the genus of
complete intersections,” Funkts. Anal. Pril., vol. 12, no. 1,
pp. 51–61, 1978.
[39] D. Mehta, “Finding All the Stationary Points of a Potential
Energy Landscape via Numerical Polynomial Homotopy Continuation Method,” Phys.Rev., vol. E84, pp. 025702, 2011.
[40] William Fulton, Intersection Theory, Springer New York, Jan.
1998.
[41] Tien-Yien Li, Tim Sauer, and James A Yorke, “Numerical
solution of a class of deﬁcient polynomial systems,” SIAM
journal on numerical analysis, vol. 24, no. 2, pp. 435–451,
1987.
[42] Jakub Marecek, Timothy McCoy, and Martin Mevissen,
“Power ﬂow as an algebraic system,”
arXiv preprint
arXiv:1412.8054, 2014.
[43] T L Lee, T Y Li, and C H Tsai, “Hom4ps-2.0, a software
package for solving polynomial systems by the polyhedral
homotopy continuation method,” Computing, vol. 83, pp. 109–
133, 2008.
[44] T. Chen, Tsung-Lin Lee, and Tien-Yien Li, “Hom4ps-3: A
parallel numerical solver for systems of polynomial equations
based on polyhedral homotopy continuation methods,” in
Mathematical Software–ICMS 2014, pp. 183–190. Springer,
2014.
[45] D. J. Bates, J. D. Hauenstein, A. J. Sommese, and C. W.
Wampler, ,” Available at www.nd.edu/∼sommese/bertini.
[46] D.J. Bates, J.D. Hauenstein, A.J. Sommese, and C.W.
Wampler,
Numerically solving polynomial systems with
Bertini, vol. 25, SIAM, 2013.
[47] T Y Li, “Solving polynomial systems by the homotopy
continuation method,” Handbook of numerical analysis, vol.
XI, pp. 209–304, 2003.
[48] A.J. Sommese and C.W. Wampler, The Numerical Solution of
Systems of Polynomials Arising in Engineering and Science,
World Scientiﬁc Publishing, Hackensack, NJ, 2005.
[49] Paul Erdős and Alfréd Rényi, “On random graphs,” Publicationes Mathematicae, vol. 6, no. 290-297, pp. 5, 1959.
[50] Paul Erdős and Alfréd Rényi, “On the strength of connectedness of a random graph,” Acta Mathematica Hungarica, vol.
12, no. 1-2, pp. 261–267, 1961.
[51] Edgar N Gilbert, “Random graphs,” The Annals of Mathematical Statistics, pp. 1141–1144, 1959.
[52] Zachary A Griﬃn and Jonathan D Hauenstein, “Real solutions
to systems of polynomial equations and parameter continuation,” Advances in Geometry, vol. 15, no. 2, pp. 173–187, 2015.
[53] Tianran Chen and Dhagash Mehta, “On the network topology
dependent solution count of the algebraic load ﬂow equations,”
arXiv preprint arXiv:1512.04987, 2015.
[54] Dhagash Mehta, Daniel K Molzahn, and Konstantin Turitsyn,
“Recent advances in computational methods for the power
ﬂow equations,” arXiv preprint arXiv:1510.00073, 2015.
[55] Daniel K Molzahn, Dhagash Mehta, and Matthew Niemerg,
“Toward topologically based upper bounds on the number of
power ﬂow solutions,” arXiv preprint arXiv:1509.09227, 2015.

