#### Abstract

Finding vacua for the four-dimensional effective theories for supergravity which descend from flux compactifications and analyzing them according to their stability is one of the central problems in string phenomenology. Except for some simple toy models, it is, however, difficult to find all the vacua analytically. Recently developed algorithmic methods based on symbolic computer algebra can be of great help in the more realistic models. However, they suffer from serious algorithmic complexities and are limited to small system sizes. In this paper, we review a numerical method called the numerical polynomial homotopy continuation (NPHC) method, first used in the areas of lattice field theories, which by construction finds *all* of the vacua of a given potential that is known to have only isolated solutions. The NPHC method is known to suffer from no major algorithmic complexities and is *embarrassingly parallelizable*, and hence its applicability goes way beyond the existing symbolic methods. We first solve a simple toy model as a warm-up example to demonstrate the NPHC method at work. We then show that all the vacua of a more complicated model of a compactified M theory model, which has an structure,
can be obtained by using a desktop machine in just about an hour, a feat which was reported to be prohibitively difficult by the existing symbolic methods. Finally, we compare the various technicalities between the two methods.

#### 1. Introduction

A lot of current research in string phenomenology is focused on developing methods to find and analyze vacua of four-dimensional effective theories for supergravity descended from flux compactifications. Stated in explicit terms, one is interested in finding all the vacua (usually, isolated stationary points) of the scalar potential of such a theory. In particular, given a Kähler potential and a superpotential , for uncharged moduli fields, the scalar potential is given by where is the Kähler derivative and is the inverse of . Once the vacua are found, one can then classify them by either using the eigenvalues of the Hessian matrix of or by introducing further constraints such as .

Finding all the stationary points of a given potential , amounts to solving the stationary equations, that is, solving the system of equations consisting of the first derivatives of , with respect to all the fields, equated to zero. The stationary equations for arising in the string phenomenological models are usually nonlinear. In the perturbative limit, usually has a polynomial form. This is an important observation since we can then use the algebraic geometry concepts and methods to extract a lot of information about . Solving systems of nonlinear equations is usually a highly nontrivial task. However, if the system of stationary equations has polynomial-like nonlinearity, then the symbolic methods based on the Gröbner basis technique can be used to solve the system [1]. These symbolic methods ensure that all the stationary points are obtained when the computation finishes. Roughly speaking, for a given system of multivariate polynomial equations, a set of which is called an ideal, the so-called Buchberger Algorithm (BA) or its refined variants can compute a new system of equations, called a Gröbner basis [2]. For the systems known to have only isolated solutions, called 0-dimensional ideals, a Gröbner basis always has at least one univariate equation and the subsequent equations consist of increasing number of variables, that is, it is in a *triangular form *(Note that this is only true for a few specific types of monomial orderings. For other monomial ordering, the new system of equations may not have a triangular form.) The solutions of a Gröbner basis are always the same as the original system, but the former is easier to solve due to its triangular form as the univariate equation can be solved either analytically or numerically quite straightforwardly. Then by back-substituting the solutions in the subsequent equations and continually solving them we can find all the solutions of the system. (Using the Gröbner basis methods, one can also deal with systems which have at least one free variable, called positive dimensional ideals. However, in this paper we only focus on the 0-dimensional ideals.) It should be noted that the BA reduces to Gaussian elimination in the case of linear equations, that is, it is a generalization of the latter. Similarly it is also a generalization of the Euclidean algorithm for the computation of the Greatest Common Divisors of a univariate polynomial. Recently, more efficient variants of the BA have been developed to obtain a Gröbner basis, for example, F4 [3], F5 [4], and Involution Algorithms [5]. Symbolic computation packages such as Mathematica, Maple, and Reduce, have built-in commands to calculate a Gröbner basis. Singular [6], COCOA [7], and MacCaulay2 [8] are specialized packages for Gröbner basis and Computational Algebraic Geometry, available as freeware. MAGMA [9] is also such a specialized package available commercially.

In [10–12], it was shown that one does not need to solve the system using the Gröbner basis techniques, in the usual sense, in order to extract some of the important information such as the dimensionality of the ideal and the number of real roots in the system. But one can indirectly obtain this information by computing the so-called primary decomposition of the ideal (still using the Gröbner basis technique internally). This was a remarkable success as it allowed one to work on nontrivial models and extract a lot of information using a regular desktop machine only. The authors of these papers also made a very helpful computational package, called Stringvacua [11], publicly available. Stringvacua is a Mathematica interface to Singular and has string phenomenology-specific utilities which makes the package quite useful to the users.

However, even with such tricks, there are a few problems with the symbolic methods: the BA is known to suffer from exponential space complexity, that is, the memory (Random Access Memory) required by the machine blows up exponentially with the number of variables, equations, terms in each polynomial, and so forth. So even for small sized systems, one may not be able to compute a Gröbner basis, nor the related objects such as primary decomposition of the ideal. It is also usually less efficient for systems with irrational coefficients. Another drawback is that the BA is highly sequential, that is, very difficult to efficiently parallelize.

Below we explain a novel numerical method, called the numerical polynomial homotopy continuation (NPHC) method, which overcomes all the shortcomings of the Gröbner basis methods. The method was first introduced in particle physics and condensed matter theory areas in [13–15], where all the stationary points of a multivariate function called the lattice Landau gauge fixing functional [16–19] were found using the NPHC method. Below, we begin by describing the NPHC method for the univariate case before generalizing it to the multivariate case. We then consider a toy model that is used in the Stringvacua manual and also a compactified M theory model. Finding all the vacua using the symbolic methods for both these models is already known to be prohibitively difficult. We briefly describe the models and explain how the corresponding stationary equations can be viewed as having polynomial form. With the help of the NPHC method, we find all the isolated vacua for the model and give a technical comparison between both the symbolic and numerical methods. After mentioning a few other important aspects of the NPHC method in the Frequently Asked Questions section, we conclude the paper.

#### 2. The Numerical Polynomial Homotopy Continuation Method

Here, we explain the numerical polynomial homotopy continuation method. Let us begin by exemplifying the method for the univariate case.

Firstly, we know that for a single-variable equation, , with coefficients and the variable both defined over , the number of solutions is exactly if , counting multiplicities. This powerful result comes from the Fundamental Theorem of Algebra. To get all roots of such single-variable polynomials, there exist many numerical methods such as the *companion matrix* trick for low-degree polynomials and the divide-and-conquer techniques for high-degree polynomials. Here we present the Numerical Polynomial Homotopy Continuation (NPHC) by first describing it for the univariate case which can then be extended to the multivariate case in a straightforward manner. We follow [20, 21] throughout this section unless specified otherwise.

The strategy behind the NPHC method is as follows: first write down the equation or system of equations to be solved in a more general parametric form, solve this system at a point in parameter space where its solutions can be easily found, and finally track these solutions from this point in parameter space to the point in parameter space corresponding to the original system/problem. This approach can be applied to many types of equations (e.g., nonalgebraic equations) which exhibit a continuous dependence of the solutions on the parameters, but there exist many difficulties in making this method a primary candidate method to solve a set of nonalgebraic equations. However, for reasons that will be clear below, this method works exceptionally well for polynomial equations.

To clarify how the method works, we first take a univariate polynomial, say , pretending that we do not know its solutions (i.e., ). We then begin by defining the more general parametric family
where is a parameter. For , we have and at we recover our original problem. The problem of getting all solutions of the original problem now reduces to tracking solutions of from where we know the solutions, that is, , to . The choice of in (2.1), called the *start system*, should be clear now: this system has the same number of solutions as the original problem and is easy to solve. For multivariate systems, a clever choice of a start system is essential in reducing the computation, and the discussion about this issue will follow soon. Here, we briefly mention the numerical methods used in path-tracking from to . One of the ways to track the paths is to solve the differential equation that is satisfied along all solution paths, say for the th solution path,
This equation is called the Davidenko differential equation. Inserting (2.1) in this equation, we have
We can solve this initial value problem numerically (again, pretending that an exact solution is hitherto unknown) with the initial conditions as and . The other approach is to use Euler's predictor and Newton's corrector methods. This approach works well too. We do not intend to discuss the actual path tracker algorithm used in practice, but it is important to mention that in these path tracker algorithms, almost all apparent difficulties have been resolved, such as tracking singular solutions, multiple roots, solutions at infinity, and so forth. It is also important to mention here that in the actual path tracker algorithms the homotopy is randomly complexified to avoid singularities, that is, taking
where with chosen randomly.

It is shown that for a generic value of the complex the paths are well behaved for , that is, for the whole path except the endpoint. This makes sure that there is no singularity or bifurcation along the paths. This is a remarkable trick, called the -trick, since this is the reason why we can claim that the NPHC method is guaranteed to find *all* solutions. Note that , for example, is not a generic value.

There are several sophisticated numerical packages well equipped with path trackers such as Bertini [22], PHCpack [23], PHoM [24], HOMPACK [25], and HOM4PS2 [26, 27]. They all are available freely from their respective research groups.

In the above example, the PHCpack with its default settings gives the solutions Thus, it gives the expected two solutions of the system with a very high numerical precision.

##### 2.1. Multivariate Polynomial Homotopy Continuation

We can now generalize the NPHC method to find all the solutions of a system of multivariate polynomial equations, say , where and , that is *known to have isolated solutions* (i.e., a 0-dimensional ideal). To do so, we first need to have some knowledge about the expected number of solutions of the system. There is a classical result, called the *Classical Bèzout Theorem*, that asserts that for a system of polynomial equations in variables the maximum number of solutions in is , where is the degree of the th polynomial. This bound, called the Classical Bèzout Bound (CBB), is exact for generic values (i.e., roughly speaking, nonzero random values) of coefficients. The *genericity* is well defined, and the interested reader is referred to [21] for details.

Based on the CBB, we can construct a *homotopy*, or a set of problems, similar to the aforementioned one-dimensional case, as
where is a system of polynomial equations, with the following properties (1)The solutions of are known or can be easily obtained. is called the *start system,* and the solutions are called the *start solutions*. (2)The number of solutions of is equal to the CBB for . (3)The solution set of for consists of a finite number of smooth paths, called homotopy paths, each parametrized by . (4)Every isolated solution of can be reached by some path originating at a solution of .

We can then track all of the paths corresponding to each solution of from to and reach . By implementing an efficient path tracker algorithm, we can get all the isolated solutions of a system of multivariate polynomials just as in the univariate case.

The homotopy constructed using the CBB is called the *Total Degree Homotopy*. The start system can be taken, for example, as
where is the degree of the polynomial of the original system . Equation (2.7) can be easily solved and its total number of solutions (the start solutions) is , all of which are nonsingular. The Total Degree Homotopy is a very effective and popular homotopy whose variants are used in the actual path trackers.

For the multivariate case, a solution is a set of numerical values of the variables which satisfies each of the equations within a given tolerance, ( in our set up). Since the variables are allowed to take complex values, all the solutions come with real and imaginary parts. A solution is a real solution if the imaginary part of each of the variables is less than or equal to a given tolerance, ( is a suitable choice for the equations we will be dealing with in the next section, below which the number of real solutions does not change). All of these solutions can be further refined to an *arbitrary precision* limited by the machine precision.

The obvious question at this stage would be if the number of real solutions depends on . To resolve this issue, we use a recently developed algorithm called alphaCertified which is based on the so-called Smale's -theory [28]. This algorithm certifies the real nonsingular solutions of polynomial systems using both exact rational arithmetic and arbitrary precision floating point arithmetic. This is a remarkable step, because using alphaCertified we can prove that a solution classified as a real solution is actually a real solution independent of , and hence these solutions are as good as the *exact solutions*.

#### 3. A Toy Model

Here, we apply the NPHC method to a toy model from the examples given in the Stringvacua package. The Kähler potential for this model is given as and the superpotential is given as Here, and are parameters. Note that the field comes along with its complex conjugate. So even though they can be treated as different variables by merely relabeling them, they are not actually independent variables. To avoid this problem, we can write them in terms of real and imaginary parts, that is, with , and are real. Finally, we get the potential as which has 2 variables. To find the stationary points of , we need to solve the system of equations consisting of the first-order derivatives of , with respect to both variables and , equated to zero, that is,

We also note that the stationary equations in this example involve denominators. Since we are not interested in the solutions for which the denominators are zero, we clear them out by multiplying them with the numerators appropriately.

Using the symbolic methods, this task is known to be difficult for general numerical (i.e., floating points) values of parameters and , with the computation continuing indefinitely [1, 12].

Firstly, we used the Stringvacua package to compute the dimension of the ideal which turned out to be 0 for generic values of and , that is, the system of equations has only isolated solutions. Note that to actually find the solutions of the system, we have to put some numerical values for and . The Gröbner basis techniques, as mentioned above, work much better for the cases where parameters are rational. We first use the same values, and , as used in the Stringvacua manual. Then, we use the command “NumRoots” which computes the number of real roots of the system, that is, 7 in this case, in less than a minute on a desktop machine.

Let us now turn our attention to solving this system using the NPHC method. Firstly, the CBB for this system is 182. We used both Bertini and HOM4PS2 to track all these paths. Both took around one minute to solve this system: there are 86 complex (including real) finite solutions, out of which 36 solutions are real. Out of the 36 real solutions, six of them are distinct solutions (multiplicity one) and the only other distinct solution which comes with multiplicity 30. Thus, there are 7 distinct solutions as expected from the Stringvacua's “NumRoots” command. However, we should mention that the Stringvacua package does not give any information about the multiplicity of the solutions, as seen in this example, whereas the NPHC method gives all the solutions with its multiplicities making the method already useful for this simple example. Not only that, but the NPHC also gives the infinite solutions (which are the solutions on the projective space but not on the affine space): the running example has 2 infinite solutions both coming with multiplicity 48. Thus, the total number of solutions in this case, , is indeed the same as the CBB.

Note that in these equations all the denominators were multiples of . The condition that none of the denominators is zero can be imposed algebraically by adding a constraint equation as with being an additional variable. Thus there are now 3 equations in 3 variables. Note that in the Stringvacua package the denominators are thrown away by multiplying each equation appropriately, but the additional equation is not included in the final ideal. In the package, one can of course use the “Saturation” command in order to ensure that this equation is properly taken into account.

We can again solve the above system 3 equations in 3 variables using the Bertini and HOM4PS2. The CBB of this new system is 364. In the end, there are 56 finite complex solutions out of which there are six real solutions, all with multiplicity 1. There are no infinite solutions in this case. This should be expected since the only multiple real solution in the previous system was when the denominator was zero. After adding the constraint equation, we have got rid of this solution and hence left with the rest of the six distinct solutions. Finally, the real solutions (throwing the very small imaginary parts out) are

Since we have all the real solutions, we can now compute the Hessian of at these solutions and separate out the physically interesting vacua. Since the purpose of this paper is to introduce the NPHC method only, we refrain from discussing the interesting physics of these solutions here. A detailed analysis of these solutions and the solutions of other systems will be published elsewhere. For now we discuss how the two methods, the symbolic algebra methods and the NPHC, compare with each other.

#### 4. A Model of Compactified M Theory

Here, we take an example of M theory compactified on the coset from [29] which is also considered in [12]. The coset has structure. The corresponding Kähler and superpotential are

Here, we use , for , and . Then the potential is

We need to solve the stationary equations, that is, the derivatives of with respect to , and equated to zero. We also need to add an additional equation to ensure that none of the denominators of the stationary equations are zero. Thus, in total there are 9 equations in 9 variables. This system only has isolated solutions, (In [12], this system is reported to have positive dimensional components in its solution space. However, the denominator equation was not included in the analysis there. Once we include the denominator equation in the system, the combined system has no positive dimensional components. Hence, there is no discrepancy here.) The equations are quite complicated, and we avoid writing all of them down here. This system of equations is not only *prohibitively difficult* to be solved completely but also not tractable even using the primary decomposition techniques (except that some information about the solutions may be obtained if one further restricts the system such as taking ) [1, 12]. In short, it is not possible to handle this system in its full glory using the available symbolic methods.

Now, let us move to the NPHC method. Firstly, since there are four equations of degree 3, another four equations of degree 4, and one equation of degree 5, the CBB is 103680. This system is actually quite straightforward to solve using the NPHC method. The HOM4PS2 package, for example, solves the full system in around 1 hour on a regular desktop machine: there are 516 total solutions for this system, out of which there are only 12 real solutions. The solutions in the order are

It is easy to recognize that some of the numbers in the above list of solutions are rational numbers, for example, . We can now easily compute the eigenvalues of the Hessian of the potential and other related quantities of all these solutions and hence classify the vacua in terms of physics. However, again we refrain from discussing the interesting physics of these solutions here. The full analysis will be published elsewhere.

#### 5. Comparison between Gröbner Basis Techniques and the NPHC Method

Here, we compare the two different methods. Firstly, the Gröbner basis techniques solve the system symbolically. This is immensely significant since one then has a *proof* for the results and/or the results in closed form. There is caveat here however: if the univariate equation in a Gröbner basis is of degree 5 or higher, then the Abel-Ruffini theorem prevents us from solving it exactly, in general, at least in terms of the radicals of its coefficients (this does not mean that the univariate equation cannot be solved exactly at all). In such a situation, one may end up solving this equation numerically, and hence the above-mentioned feature of the symbolic method no longer applies. The NPHC method is a numerical method. That said, the method by construction gives *all *of the isolated solutions for the system known to have only isolated solutions, up to a numerical precision. The solutions then can be refined to within an arbitrary precision up to the machine precision by the Newton's corrector method or otherwise. Moreover, using the alphaCertified method, we can certify if the real nonsingular solutions obtained by the above packages are actually the real nonsingular solutions of the system independent of the numerical precision used during the computation. Hence, though the solutions cannot be obtained in a closed form using the NPHC method, the solutions are as good as exact solutions for all the practical purposes.

We should emphasize here that using the methods presented in [10–12] one can learn quite a lot about a system without having to necessarily obtain its solutions. In particular, one can use the so-called primary decomposition of the ideal (though making use of the Gröbner basis technique only) to obtain information such as the dimensionality of the solution space and number of isolated real roots. This is indeed a clever way to resolve the above-mentioned issue up to a certain level. However, here, the next difficulty comes in the form of algorithmic complexity. The BA is known to suffer from exponential space complexity, which roughly means that the memory (Random Access Memory) required by the machine blows up exponentially with increasing number of polynomials, variables, monomials, and/or degree of the polynomials involved in the system. Hence, even the computation for the primary decomposition may not finish for large sized systems, whereas the NPHC method is strikingly different from the Gröbner basis techniques in that the algorithm for the former suffers from no known major complexities. Hence one can in principle find all solutions of bigger systems.

The BA is a highly sequential algorithm, that is, each step in the algorithm requires knowledge of the previous one. Thus, although recently there are certain parts of the BA which have been parallelized, in general, it is extremely difficult to parallelize the algorithm. On the other hand, in the NPHC method, the path tracking is *embarrassingly parallelizable*, because each start solution can be tracked completely independently of the others. This feature along with the rapid progress towards the improvements of the algorithms makes the NPHC well suited for a large class of physical problems arising not only in string phenomenology but in condensed matter theory, lattice QCD, and so forth.

The BA is mainly defined for systems with rational coefficients, while in real-life applications, the systems may have real coefficients. The NPHC method being a purely numerical method by default incorporates floating point coefficients as well.

In conclusion, both the Gröbner basis techniques and the NPHC have advantages and disadvantages. However, for practical purposes, the NPHC method is a far more efficient and promising method for realistic systems.

#### 6. Frequently Asked Questions

In this section, we collect the frequently asked questions and their answers. (1)What does the NPHC method tell us about systems which do not have any solutions?As mentioned above, the NPHC method by construction (in conjunction with the -trick) gives all real and complex solutions of a system of multivariate polynomial equations that is known to have only isolated solutions. Hence, we are always sure that we have got all the solutions numerically. This statement is true for all cases such as when the system has no complex solutions and/or no real solutions or no solution at all. One possible issue, as mentioned above, regards the classification of the real solutions independent of the tolerance used. This can be resolved by using, for example, the alphaCertified algorithm which certifies when a solution is a real.(2)For many practical problems, only real solutions are required. Thus, when implementing the NPHC, a huge amount of computational effort is wasted in getting the other types of solutions. Would it not be helpful to track only the real solutions?It would be much more useful if there was a way of getting only real solutions. However, for a number of technical reasons nicely discussed in [21], a path tracker does not know in advance if a given start solution will end up being a real solution of the original system. Moreover, one can wonder if a root count exists only for the real solutions of a system. This would involve obtaining a corresponding fundamental theorem of algebra on the real space for the multivariate case. This, however, has yet to be achieved. Hence, the best way for now is to track all complex (including real paths) solutions and then filter out the real solutions.(3)Can the NPHC method be used as a global or local minimization method?Absolutely, most of the conventional methods used to minimize a function are based on the Newton-Raphson method, where a start solution is guessed and is then refined by successive iterations in the direction of the minima. By performing this algorithm several times on the functions, one can obtain many minima of the given potential. Recently, more efficient methods such as the basin-hopping method are available for local minimization [30]. However, we are never sure if we have got all the minima from any of these methods. For the global minimization, we may use more efficient methods such as Simulated Annealing and Genetic Algorithm. However, these methods are known to fail for larger systems since it can easily get trapped at a local minimum when trying to find the global minimum. Thus, in addition to the usual error from the numerical precision of the machine, there can be an error of an *unknown* order (i.e., we do not know if the found one is the global minimum!). But if the function has a polynomial-like nonlinearity, in theory the NPHC method can give all the minima since it obtains all the stationary points. So it solves the local minimization problem. Moreover, it is then easy to identify *the* global minimum out of the minima, and hence we are sure that the found one is actually the global minimum.(4)As in the example system in this paper, for many systems the number of actual solutions may be well below the CBB. Is there any remedy for this issue?The main reason why the number of actual solutions is less than the CBB for many systems is that the CBB does not take the sparsity (i.e., very few monomials in each polynomial in the system) of the system into account. There is indeed a tighter upper bound on the number of complex solutions, called the Bernstein-Khovanskii-Kushnirenko (BKK) count [20, 21], which takes this sparsity into account and thus in most cases is much lower than the CBB. In many cases, it is in fact equal to the number of solutions. The BKK bound can thus save a lot of computation time since the number of paths to be tracked is less than the CBB. The details on the BKK count relating to string phenomenology problems will be published elsewhere.(5)Are there any alternative/supplementary numerical methods?There are not many methods to find the stationary points of a multivariate function around, compared to the number of methods to find minima. One of the methods that can find stationary points is the Gradient-minimization method which finds all the minima of an auxiliary function whose minima are the stationary points of provided we further restrict to be zero [31]. One can find many minima of using some conventional minimization method such as the Conjugate Gradient method or the Simulated Annealing method. However, it is known that as the system size increases, the number of minima of that are not the minima of , that is, , increases rapidly, making the method inefficient [32]. Another method is the Newton-Raphson method (and its sophisticated variants) [30, 32–34]. There, an initial guess is refined iteratively to a given precision. It should be emphasized, however, that no matter how many different random initial guesses are fed into the algorithm, we can never be sure to get all the solutions in the end, unlike the NPHC method. However, these two methods can be supplementary methods for bigger systems to get an idea on what to expect there.(6)This paper mainly deals with the potentials having polynomial-like nonlinearity which may be usual in the perturbation limit. What about the fully nonperturbative potentials?The most interesting application for this method would be in the nonperturbative regime, certainly. This question can be stated in different words: is it possible to translate the stationary equations for the nonperturbative potential (i.e., the potentials which have logarithm and exponential terms), and if so, how? Once we can translate the equations in the polynomial form, we can again use the NPHC method as before. The answer is already available in [10]. In this work, Gray et al. have already prescribed how to translate the corresponding equations arising in the nonperturbative regime which usually involve logarithms and/or exponentials, by using dummy variables. After that, we can solve the system using algebraic geometry methods, such as the Gröbner basis, or for more complicated cases, the NPHC method presented in this paper. Once we have all the solutions, we can extract the solutions in terms of the original variables which were logarithms and/or exponentials of the fields. This trick makes all the algebraic geometry methods, not only the NPHC method, applicable to finding the vacua of the potentials in the nonperturbative regime.(7)This method assumes that one knows that the system under consideration has only isolated solutions. But, in general, one may not know if a given system has only isolated solutions or it contains some positive dimensional components. In that case, do not we need to rely on the Gröbner basis techniques only, at least to check the dimension of the system?Firstly, in many systems, once we add the constraint equation (i.e., the denominators are never zero), they usually turn out to have only isolated solutions. Thus, there are way too many interesting systems in string phenomenology which only have isolated solutions. Of course, there may be many more systems which would have positive dimensional solution components. To solve such systems, there is a recently developed generalization of the numerical homotopy continuation method, called the Numerical Algebraic Geometry method. This method finds out each of the positive dimensional solution components with its dimensionality. This method is also *embarrassingly parallelizable* and hence goes far beyond the reach of the Gröbner basis methods. The details of this method are much more involved and beyond the scope of the present paper. But, in short, to find out the dimensionality of the system we do not necessarily need to rely on the Gröbner basis methods. The details of this method with applications will be published elsewhere.

#### 7. Summary

In this paper, we have reviewed a novel method, called the numerical polynomial homotopy continuation (NPHC) method, which can find all the string vacua of a given potential. It does not suffer from any major algorithmic complexities compared to the existing symbolic algebra methods based on the Gröbner basis techniques, which are known to suffer from exponential space complexity. Moreover, the NPHC method is *embarrassingly parallelizable, *making it a very efficient alternative to the existing symbolic algebra methods*. *As an example, we studied a toy model and, using the NPHC method, found all the vacua within less than a minute using a regular desktop machine. Note that this system with the irrational coefficients is already a difficult task using the Gröbner basis techniques. In addition to that, using the NPHC method, with just about an hour of computation on a regular desktop machine, we found all vacua of an M theory model compactified on the coset , which has an structure. This system was reported to be a prohibitively difficult problem using the symbolic method. Thus, we have already shown how efficiently the NPHC method can solve the problems that are yet far beyond the reach of the traditional symbolic methods. We also emphasize that using the procedure prescribed in [10] to translate the stationary equations arising in the nonperturbative regime, by replacing logarithm and exponential terms of the field variables by dummy variables, into the polynomial form, we can use the NPHC method to find the vacua for the nonperturbative potentials as well. It is this application of the method which makes it quite promising. With the help of the NPHC method it is thus hoped that we can go far beyond the reach of the existing methods and study realistic models very efficiently.

#### Acknowledgments

D. Mehata was supported by the US Department of Energy grant under contract no. DE-FG02-85ER40237 and Science Foundation Ireland grant 08/RFP/PHY1462. D. Mehata would like to thank James Gray, Yang-Hui He, and Anton Ilderton for encouraging and helping him throughout in this work.