BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Dhagash Mehta

Dhagash Mehta appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
hep-thcond-mat.stat-mechhep-latmath.AGmath-phmath.MPmath.NAmath.OChep-phmath.DSnlin.CDq-fin.CPSystems and ControlComputational Engineering, Finance, and Sciencecond-mat.dis-nncond-mat.softMachine LearningNumerical Analysisphysics.chem-phq-fin.STSymbolic ComputationApplicationsastro-ph.COCell Behaviorcond-mat.mtrl-scimath.CVmath.PRMathematical SoftwareMolecular NetworksNeural and Evolutionary ComputingNeurons and Cognitionnucl-thphysics.atm-clusphysics.comp-phQuantitative Methods

Catalog footprint

What is connected

58works
35topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

58 published item(s)

preprint2026arXiv

Deep Reinforcement Learning for Optimum Order Execution: Mitigating Risk and Maximizing Returns

Optimal Order Execution is a well-established problem in finance that pertains to the flawless execution of a trade (buy or sell) for a given volume within a specified time frame. This problem revolves around optimizing returns while minimizing risk, yet recent research predominantly focuses on addressing one aspect of this challenge. In this paper, we introduce an innovative approach to Optimal Order Execution within the US market, leveraging Deep Reinforcement Learning (DRL) to effectively address this optimization problem holistically. Our study assesses the performance of our model in comparison to two widely employed execution strategies: Volume Weighted Average Price (VWAP) and Time Weighted Average Price (TWAP). Our experimental findings clearly demonstrate that our DRL-based approach outperforms both VWAP and TWAP in terms of return on investment and risk management. The model's ability to adapt dynamically to market conditions, even during periods of market stress, underscores its promise as a robust solution.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Learning Embedded Representation of the Stock Correlation Matrix using Graph Machine Learning

Understanding non-linear relationships among financial instruments has various applications in investment processes ranging from risk management, portfolio construction and trading strategies. Here, we focus on interconnectedness among stocks based on their correlation matrix which we represent as a network with the nodes representing individual stocks and the weighted links between pairs of nodes representing the corresponding pair-wise correlation coefficients. The traditional network science techniques, which are extensively utilized in financial literature, require handcrafted features such as centrality measures to understand such correlation networks. However, manually enlisting all such handcrafted features may quickly turn out to be a daunting task. Instead, we propose a new approach for studying nuances and relationships within the correlation network in an algorithmic way using a graph machine learning algorithm called Node2Vec. In particular, the algorithm compresses the network into a lower dimensional continuous space, called an embedding, where pairs of nodes that are identified as similar by the algorithm are placed closer to each other. By using log returns of S&P 500 stock data, we show that our proposed algorithm can learn such an embedding from its correlation network. We define various domain specific quantitative (and objective) and qualitative metrics that are inspired by metrics used in the field of Natural Language Processing (NLP) to evaluate the embeddings in order to identify the optimal one. Further, we discuss various applications of the embeddings in investment management.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Learning Mutual Fund Categorization using Natural Language Processing

Categorization of mutual funds or Exchange-Traded-funds (ETFs) have long served the financial analysts to perform peer analysis for various purposes starting from competitor analysis, to quantifying portfolio diversification. The categorization methodology usually relies on fund composition data in the structured format extracted from the Form N-1A. Here, we initiate a study to learn the categorization system directly from the unstructured data as depicted in the forms using natural language processing (NLP). Positing as a multi-class classification problem with the input data being only the investment strategy description as reported in the form and the target variable being the Lipper Global categories, and using various NLP models, we show that the categorization system can indeed be learned with high accuracy. We discuss implications and applications of our findings as well as limitations of existing pre-trained architectures in applying them to learn fund categorization.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Machine learning the real discriminant locus

Parameterized systems of polynomial equations arise in many applications in science and engineering with the real solutions describing, for example, equilibria of a dynamical system, linkages satisfying design constraints, and scene reconstruction in computer vision. Since different parameter values can have a different number of real solutions, the parameter space is decomposed into regions whose boundary forms the real discriminant locus. This article views locating the real discriminant locus as a supervised classification problem in machine learning where the goal is to determine classification boundaries over the parameter space, with the classes being the number of real solutions. For multidimensional parameter spaces, this article presents a novel sampling method which carefully samples the parameter space. At each sample point, homotopy continuation is used to obtain the number of real solutions to the corresponding polynomial system. Machine learning techniques including nearest neighbor and deep learning are used to efficiently approximate the real discriminant locus. One application of having learned the real discriminant locus is to develop a real homotopy method that only tracks the real solution paths unlike traditional methods which track all~complex~solution~paths. Examples show that the proposed approach can efficiently approximate complicated solution boundaries such as those arising from the equilibria of the Kuramoto model.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Machine Learning Fund Categorizations

Given the surge in popularity of mutual funds (including exchange-traded funds (ETFs)) as a diversified financial investment, a vast variety of mutual funds from various investment management firms and diversification strategies have become available in the market. Identifying similar mutual funds among such a wide landscape of mutual funds has become more important than ever because of many applications ranging from sales and marketing to portfolio replication, portfolio diversification and tax loss harvesting. The current best method is data-vendor provided categorization which usually relies on curation by human experts with the help of available data. In this work, we establish that an industry wide well-regarded categorization system is learnable using machine learning and largely reproducible, and in turn constructing a truly data-driven categorization. We discuss the intellectual challenges in learning this man-made system, our results and their implications.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Decomposing the parameter space of biological networks via a numerical discriminant approach

Many systems in biology, physics and engineering can be described by systems of ordinary differential equation containing many parameters. When studying the dynamic behavior of these large, nonlinear systems, it is useful to identify and characterize the steady-state solutions as the model parameters vary, a technically challenging problem in a high-dimensional parameter landscape. Rather than simply determining the number and stability of steady-states at distinct points in parameter space, we decompose the parameter space into finitely many regions, the steady-state solutions being consistent within each distinct region. From a computational algebraic viewpoint, the boundary of these regions is contained in the discriminant locus. We develop global and local numerical algorithms for constructing the discriminant locus and classifying the parameter landscape. We showcase our numerical approaches by applying them to molecular and cell-network models.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Investigating the Maximum Number of Real Solutions to the Power Flow Equations: Analysis of Lossless Four-Bus Systems

The power flow equations model the steady-state relationship between the power injections and voltage phasors in an electric power system. By separating the real and imaginary components of the voltage phasors, the power flow equations can be formulated as a system of quadratic polynomials. Only the real solutions to these polynomial equations are physically meaningful. This paper focuses on the maximum number of real solutions to the power flow equations. An upper bound on the number of real power flow solutions commonly used in the literature is the maximum number of complex solutions. There exist two- and three-bus systems for which all complex solutions are real. It is an open question whether this is also the case for larger systems. This paper investigates four-bus systems using techniques from numerical algebraic geometry and conjectures a negative answer to this question. In particular, this paper studies lossless, four-bus systems composed of PV buses connected by lines with arbitrary susceptances. Computing the Galois group, which is degenerate, enables conversion of the problem of counting the number of real solutions to the power flow equations into counting the number of positive roots of a univariate sextic polynomial. From this analysis, it is conjectured that the system has at most 16 real solutions, which is strictly less than the maximum number of complex solutions, namely 20. We also provide explicit parameter values where this system has 16 real solutions so that the conjectured upper bound is achievable.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Kinetic Transition Networks for the Thomson Problem and Smale's 7th Problem

The Thomson Problem, arrangement of identical charges on the surface of a sphere, has found many applications in physics, chemistry and biology. Here we show that the energy landscape of the Thomson Problem for $N$ particles with $N=132, 135, 138, 141, 144, 147$ and $150$ is single funnelled, characteristic of a structure-seeking organisation where the global minimum is easily accessible. Algorithmically constructing starting points close to the global minimum of such a potential with spherical constraints is one of Smale's 18 unsolved problems in mathematics for the 21st century because it is important in the solution of univariate and bivariate random polynomial equations. By analysing the kinetic transition networks, we show that a randomly chosen minimum is in fact always `close' to the global minimum in terms of the number of transition states that separate them, a characteristic of small world networks.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Collection of Challenging Optimization Problems in Science, Engineering and Economics

Function optimization and finding simultaneous solutions of a system of nonlinear equations (SNE) are two closely related and important optimization problems. However, unlike in the case of function optimization in which one is required to find the global minimum and sometimes local minima, a database of challenging SNEs where one is required to find stationary points (extrama and saddle points) is not readily available. In this article, we initiate building such a database of important SNE (which also includes related function optimization problems), arising from Science, Engineering and Economics. After providing a short review of the most commonly used mathematical and computational approaches to find solutions of such systems, we provide a preliminary list of challenging problems by writing the Mathematical formulation down, briefly explaning the origin and importance of the problem and giving a short account on the currently known results, for each of the problems. We anticipate that this database will not only help benchmarking novel numerical methods for solving SNEs and function optimization problems but also will help advancing the corresponding research areas.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Algebraic Geometrization of the Kuramoto Model: Equilibria and Stability Analysis

Finding equilibria of the finite size Kuramoto model amounts to solving a nonlinear system of equations, which is an important yet challenging problem. We translate this into an algebraic geometry problem and use numerical methods to find all of the equilibria for various choices of coupling constants K, natural frequencies, and on different graphs. We note that for even modest sizes (N ~ 10-20), the number of equilibria is already more than 100,000. We analyze the stability of each computed equilibrium as well as the configuration of angles. Our exploration of the equilibrium landscape leads to unexpected and possibly surprising results including non-monotonicity in the number of equilibria, a predictable pattern in the indices of equilibria, counter-examples to popular conjectures, multi-stable equilibrium landscapes, scenarios with only unstable equilibria, and multiple distinct extrema in the stable equilibrium distribution as a function of the number of cycles in the graph.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An index-resolved fixed-point homotopy and potential energy landscapes

Stationary points (SPs) of the potential energy landscapes can be classified by their Morse index, i.e., the number of negative eigenvalues of the Hessian evaluated at the SPs. In understanding chemical clusters through their potential energy landscapes, only SPs of a particular Morse index are needed. We propose a modification of the "fixed-point homotopy" method which can be used to directly target stationary points of a specified Morse index. We demonstrate the effectiveness of our approach by applying it to the Lennard-Jones clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Global Structure of Curves from Generalized Unitarity Cut of Three-loop Diagrams

This paper studies the global structure of algebraic curves defined by generalized unitarity cut of four-dimensional three-loop diagrams with eleven propagators. The global structure is a topological invariant that is characterized by the geometric genus of the algebraic curve. We use the Riemann-Hurwitz formula to compute the geometric genus of algebraic curves with the help of techniques involving convex hull polytopes and numerical algebraic geometry. Some interesting properties of genus for arbitrary loop orders are also explored where computing the genus serves as an initial step for integral or integrand reduction of three-loop amplitudes via an algebraic geometric approach.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the Network Topology Dependent Solution Count of the Algebraic Load Flow Equations

A large amount of research activity in power systems areas has focused on developing computational methods to solve load flow equations where a key question is the maximum number of isolated solutions.Though several concrete upper bounds exist, recent studies have hinted that much sharper upper bounds that depend the topology of underlying power networks may exist. This paper establishes such a topology dependent solution bound which is actually the best possible bound in the sense that it is always attainable. We also develop a geometric construction called adjacency polytope which accurately captures the topology of the underlying power network and is immensely useful in the computation of the solution bound. Finally we highlight the significant implications of the development of such solution bound in solving load flow equations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Parallel degree computation for solution space of binomial systems with an application to the master space of $\mathcal{N}=1$ gauge theories

The problem of solving a system of polynomial equations is one of the most fundamental problems in applied mathematics. Among them, the problem of solving a system of binomial equations form a important subclass for which specialized techniques exist. For both theoretic and applied purposes, the degree of the solution set of a system of binomial equations often plays an important role in understanding the geometric structure of the solution set. Its computation, however, is computationally intensive. This paper proposes a specialized parallel algorithm for computing the degree on GPUs that takes advantage of the massively parallel nature of GPU devices. The preliminary implementation shows remarkable efficiency and scalability when compared to the closest CPU-based counterpart. Applied to the "master space problem of $\mathcal{N}=1$ gauge theories" the GPU-based implementation achieves nearly 30 fold speedup over its CPU-only counterpart enabling the discovery of previously unknown results. Equally important to note is the far superior scalability: with merely 3 GPU devices on a single workstation, the GPU-based implementation shows better performance, on certain problems, than a small cluster totaling 100 CPU cores.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Recent Advances in Computational Methods for the Power Flow Equations

The power flow equations are at the core of most of the computations for designing and operating electric power systems. The power flow equations are a system of multivariate nonlinear equations which relate the power injections and voltages in a power system. A plethora of methods have been devised to solve these equations, starting from Newton-based methods to homotopy continuation and other optimization-based methods. While many of these methods often efficiently find a high-voltage, stable solution due to its large basin of attraction, most of the methods struggle to find low-voltage solutions which play significant role in certain stability-related computations. While we do not claim to have exhausted the existing literature on all related methods, this tutorial paper introduces some of the recent advances in methods for solving power flow equations to the wider power systems community as well as bringing attention from the computational mathematics and optimization communities to the power systems problems. After briefly reviewing some of the traditional computational methods used to solve the power flow equations, we focus on three emerging methods: the numerical polynomial homotopy continuation method, Groebner basis techniques, and moment/sum-of-squares relaxations using semidefinite programming. In passing, we also emphasize the importance of an upper bound on the number of solutions of the power flow equations and review the current status of research in this direction.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Statistics of Stationary Points of Random Finite Polynomial Potentials

The stationary points (SPs) of the potential energy landscapes (PELs) of multivariate random potentials (RPs) have found many applications in many areas of Physics, Chemistry and Mathematical Biology. However, there are few reliable methods available which can find all the SPs accurately. Hence, one has to rely on indirect methods such as Random Matrix theory. With a combination of the numerical polynomial homotopy continuation method and a certification method, we obtain all the certified SPs of the most general polynomial RP for each sample chosen from the Gaussian distribution with mean 0 and variance 1. While obtaining many novel results for the finite size case of the RP, we also discuss the implications of our results on mathematics of random systems and string theory landscapes.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Toward Topologically Based Upper Bounds on the Number of Power Flow Solutions

The power flow equations, which relate power injections and voltage phasors, are at the heart of many electric power system computations. While Newton-based methods typically find the "high-voltage" solution to the power flow equations, which is of primary interest, there are potentially many "low-voltage" solutions that are useful for certain analyses. This paper addresses the number of solutions to the power flow equations. There exist upper bounds on the number of power flow solutions; however, there is only limited work regarding bounds that are functions of network topology. This paper empirically explores the relationship between the network topology, as characterized by the maximal cliques, and the number of power flow solutions. To facilitate this analysis, we use a numerical polynomial homotopy continuation approach that is guaranteed to find all complex solutions to the power flow equations. The number of solutions obtained from this approach upper bounds the number of real solutions. Testing with many small networks informs the development of upper bounds that are functions of the network topology. Initial results include empirically derived expressions for the maximum number of solutions for certain classes of network topologies.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An Inversion-Relaxation Approach for Sampling Stationary Points of Spin Model Hamiltonians

Sampling the stationary points of a complicated potential energy landscape is a challenging problem. Here we introduce a sampling method based on relaxation from stationary points of the highest index of the Hessian matrix. We illustrate how this approach can find all the stationary points for potentials or Hamiltonians bounded from above, which includes a large class of important spin models, and we show that it is far more efficient than previous methods. For potentials unbounded from above, the relaxation part of the method is still efficient in finding minima and transition states, which are usually the primary focus of attention for atomistic systems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Certification and the Potential Energy Landscape

Typically, there is no guarantee that a numerical approximation obtained using standard nonlinear equation solvers is indeed an actual solution, meaning that it lies in the quadratic convergence basin. Instead, it may lie only in the linear convergence basin, or even in a chaotic region, and hence not converge to the corresponding stationary point when further optimization is attempted. In some cases, these non-solutions could be misleading. Proving that a numerical approximation will quadratically converge to a stationary point is termed \textit{certification}. In this report, we provide details of how Smale's $α$-theory can be used to certify numerically obtained stationary points of a potential energy landscape, providing a \textit{mathematical proof} that the numerical approximation does indeed correspond to an actual stationary point, independent of the precision employed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Certified counting of roots of random univariate polynomials

A challenging problem in computational mathematics is to compute roots of a high-degree univariate random polynomial. We combine an efficient multiprecision implementation for solving high-degree random polynomials with two certification methods, namely Smale's $α$-theory and one based on Gerschgorin's theorem, for showing that a given numerical approximation is in the quadratic convergence region of Newton's method of some exact solution. With this combination, we can certifiably count the number of real roots of random polynomials. We quantify the difference between the two certification procedures and list the salient features of both of them. After benchmarking on random polynomials where the coefficients are drawn from the Gaussian distribution, we obtain novel experimental results for the Cauchy distribution case.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Communication: Newton homotopies for sampling stationary points of potential energy landscapes

One of the most challenging and frequently arising problems in many areas of science is to find solutions of a system of multivariate nonlinear equations. There are several numerical methods that can find many (or all if the system is small enough) solutions but they each exhibit characteristic problems. Moreover, traditional methods can break down if the system contains singular solutions. Here, we propose an efficient implementation of Newton homotopies, which can sample a large number of the stationary points of complicated many-body potentials. We demonstrate how the procedure works by applying it to the nearest-neighbor $ϕ^4$ model and atomic clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Energy Landscape of the Finite-Size Mean-field 2-Spin Spherical Model and Topology Trivialization

Motivated by the recently observed phenomenon of topology trivialization of potential energy landscapes (PELs) for several statistical mechanics models, we perform a numerical study of the finite size $2$-spin spherical model using both numerical polynomial homotopy continuation and a reformulation via non-hermitian matrices. The continuation approach computes all of the complex stationary points of this model while the matrix approach computes the real stationary points. Using these methods, we compute the average number of stationary points while changing the topology of the PEL as well as the variance. Histograms of these stationary points are presented along with an analysis regarding the complex stationary points. This work connects topology trivialization to two different branches of mathematics: algebraic geometry and catastrophe theory, which is fertile ground for further interdisciplinary research.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Enumerating Copies in the First Gribov Region on the Lattice in up to four Dimensions

The covariant gauges are known to suffer from the Gribov problem: even after fixing a gauge non-perturbatively, there may still exist residual copies which are physically equivalent to each other, called Gribov copies. While the influence of Gribov copies in the relevant quantities such as gluon propagators has been heavily debated in recent studies, the significance of the role they play in the Faddeev--Popov procedure is hardly doubted. We concentrate on Gribov copies in the first Gribov region, i.e., the space of Gribov copies at which the Faddeev--Popov operator is strictly positive (semi)definite. We investigate compact U($1$) as the prototypical model of the more complicated standard model group SU($N_{c}$). With our Graphical Processing Unit (GPU) implementation of the relaxation method we collect up to a few million Gribov copies per orbit. We show that the numbers of Gribov copies even in the first Gribov region increase exponentially in two, three and four dimensions. Furthermore, we provide strong indication that the number of Gribov copies is gauge orbit dependent.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Experiments on the zeros of harmonic polynomials using certified counting

Motivated by Wilmshurst's conjecture, we investigate the zeros of harmonic polynomials. We utilize a certified counting approach which is a combination of two methods from numerical algebraic geometry: numerical polynomial homotopy continuation to compute a numerical approximation of each zero and Smale's alpha-theory to certify the results. Using this approach, we provide new examples of harmonic polynomials having the most extreme number of zeros known so far; we also study the mean and variance of the number of zeros of random harmonic polynomials.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Exploring the Impact of Wind Penetration on Power System Equilibrium Using a Numerical Continuation Approach

In this paper we investigate how the equilibrium characteristics of conventional power systems may change with an increase in wind penetration. We first derive a differential-algebraic model of a power system network consisting of synchronous generators, loads and a wind power plant modeled by a wind turbine and a doubly-fed induction generator (DFIG). The models of these three components are coupled via nonlinear power flow equations. In contrast to the traditional approach for solving the power flows via iterative methods that often lead to only local solutions, we apply a recently developed parameter-homotopy based numerical continuation algorithm to compute all possible solutions. The method solves the power flow equations over multiple values of the wind penetration level with far less computational effort instead of solving them at each value individually. We observe that depending on the penetration limit and the setpoint value for the magnitude of the wind bus voltage, the system may exhibit several undesired or even unstable equilibria. We illustrate these results through a detailed simulation of a 5-machine power system model with wind injection, and highlight how the solutions may be helpful for small-signal stability assessment.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gauge-fixing on the Lattice via Orbifolding

When fixing a covariant gauge, most popularly the Landau gauge, on the lattice one encounters the Neuberger 0/0 problem which prevents one from formulating a Becchi--Rouet--Stora--Tyutin symmetry on the lattice. Following the interpretation of this problem in terms of Witten-type topological field theory and using the recently developed Morse theory for orbifolds, we propose a modification of the lattice Landau gauge via orbifolding of the gauge-fixing group manifold and show that this modification circumvents the orbit-dependence issue and hence can be a viable candidate for evading the Neuberger problem. Using algebraic geometry, we also show that though the previously proposed modification of the lattice Landau gauge via stereographic projection relies on delicate departure from the standard Morse theory due to the non-compactness of the underlying manifold, the corresponding gauge-fixing partition function turns out to be orbit independent for all the orbits except in a region of measure zero.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Numerical Polynomial Homotopy Continuation Method to Locate All The Power Flow Solutions

The manuscript addresses the problem of finding all solutions of power flow equations or other similar nonlinear system of algebraic equations. This problem arises naturally in a number of power systems contexts, most importantly in the context of direct methods for transient stability analysis and voltage stability assessment. We introduce a novel form of homotopy continuation method called the numerical polynomial homotopy continuation (NPHC) method that is mathematically guaranteed to find all the solutions without ever encountering a bifurcation. The method is based on embedding the real form of power flow equation in complex space, and tracking the generally unphysical solutions with complex values of real and imaginary parts of the voltage. The solutions converge to physical real form in the end of the homotopy. The so-called $γ$-trick mathematically rigorously ensures that all the paths are well-behaved along the paths, so unlike other continuation approaches, no special handling of bifurcations is necessary. The method is \textit{embarrassingly parallelizable} and can be applied to reasonably large sized systems. We demonstrate the technique by analysis of several standard test cases up to the 14-bus system size. Finally, we discuss possible strategies for scaling the method to large size systems, and propose several applications for transient stability analysis and voltage stability assessment.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On exact minimization of Higgs potentials

Minimizing the Higgs potential is an essential task in any model involving Higgs bosons. Exact minimization methods proposed in the literature are based on the polynomial form of the potential. These methods will in general no longer work if loop contributions to the potential are taken into account. We present a method to keep the tree level global minimum unchanged in passing to the effective potential. We illustrate the method for the case of the Minimal Supersymmetric Model (MSSM).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Potential Energy Landscape of the Two-Dimensional XY Model: Higher-Index Stationary Points

The application of numerical techniques to the study of energy landscapes of large systems relies on sufficient sampling of the stationary points. Since the number of stationary points is believed to grow exponentially with system size, we can only sample a small fraction. We investigate the interplay between this restricted sample size and the physical features of the potential energy landscape for the two-dimensional $XY$ model in the absence of disorder with up to $N=100$ spins. Using an eigenvector-following technique, we numerically compute stationary points with a given Hessian index $I$ for all possible values of $I$. We investigate the number of stationary points, their energy and index distributions, and other related quantities, with particular focus on the scaling with $N$. The results are used to test a number of conjectures and approximate analytic results for the general properties of energy landscapes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Certifying the Potential Energy Landscape

It is highly desirable for a numerical approximation of a stationary point for a potential energy landscape to lie in the quadratic convergence basin of that stationary point. However, it is possible that an approximation may lie only in the linear convergence basin, or even in a chaotic region, and hence not converge to the actual stationary point when further optimization is attempted. Proving that a numerical approximation will quadratically converge to the associated stationary point is termed certifying the numerical approximation. We employ Smale's α-theory to stationary points, providing a certification that serves as a mathematical proof that the numerical approximation does indeed correspond to an actual stationary point, independent of the precision employed. As a practical example, employing recently developed certification algorithms, we show how the α-theory can be used to certify all the known minima and transition states of Lennard-Jones LJ$_{N}$ atomic clusters for N = 7, ...,14.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Energy Landscape of the Finite-Size Mean-field 3-Spin Spherical Model

We study the 3-spin spherical model with mean-field interactions and Gaussian random couplings. For moderate system sizes of up to 20 spins, we obtain all stationary points of the energy landscape by means of the numerical polynomial homotopy continuation method. On the basis of these stationary points, we analyze the complexity and other quantities related to the glass transition of the model and compare these finite-system quantities to their exact counterparts in the thermodynamic limit.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Enumerating Gribov copies on the lattice

In the modern formulation of lattice gauge-fixing, the gauge fixing condition is written in terms of the minima or stationary points (collectively called solutions) of a gauge-fixing functional. Due to the non-linearity of this functional, it usually has many solutions called Gribov copies. The dependence of the number of Gribov copies, n[U] on the different gauge orbits plays an important role in constructing the Faddeev-Popov procedure and hence in realising the BRST symmetry on the lattice. Here, we initiate a study of counting n[U] for different orbits using three complimentary methods: 1. analytical results in lower dimensions, and some lower bounds on n[U] in higher dimensions, 2. the numerical polynomial homotopy continuation method, which numerically finds all Gribov copies for a given orbit for small lattices, and 3. numerical minimisation ("brute force"), which finds many distinct Gribov copies, but not necessarily all. Because n for the coset SU(N_c)/U(1) of an SU(N_c) theory is orbit-independent, we concentrate on the residual compact U(1) case in this article and establish that n is orbit-dependent for the minimal lattice Landau gauge and orbit-independent for the absolute lattice Landau gauge. We also observe that contrary to a previous claim, n is not exponentially suppressed for the recently proposed stereographic lattice Landau gauge compared to the naive gauge in more than one dimension.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Exploring the Potential Energy Landscape Over a Large Parameter-Space

Solving large polynomial systems with coefficient parameters are ubiquitous and constitute an important class of problems. We demonstrate the computational power of two methods--a symbolic one called the Comprehensive Gröbner basis and a numerical one called the cheater's homotopy-applied to studying both potential energy landscapes and a variety of questions arising from geometry and phenomenology. Particular attention is paid to an example in flux compactification where important physical quantities such as the gravitino and moduli masses and the string coupling can be efficiently extracted.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Numerical Analyses on Moduli Space of Vacua

We propose a new computational method to understand the vacuum moduli space of (supersymmetric) field theories. By combining numerical algebraic geometry (NAG) and elimination theory, we develop a powerful, efficient, and parallelizable algorithm to extract important information such as the dimension, branch structure, Hilbert series and subsequent operator counting, as well as variation according to coupling constants and mass parameters. We illustrate this method on a host of examples from gauge theory, string theory, and algebraic geometry.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Phase transitions and gluodynamics in 2-colour matter at high density

We investigate 2-colour QCD with 2 flavours of Wilson fermion at nonzero temperature T and quark chemical potential mu, with a pion mass of 700 MeV (m_pi/m_rho=0.8). From temperature scans at fixed mu we find that the critical temperature for the superfluid to normal transition depends only very weakly on mu above the onset chemical potential, while the deconfinement crossover temperature is clearly decreasing with mu. We find indications of a region of superfluid but deconfined matter at high mu and intermediate T. The static quark potential determined from the Wilson loop is found to exhibit a 'string tension' that increases at large mu in the 'deconfined' region. The electric (longitudinal) gluon propagator in Landau gauge becomes strongly screened with increasing temperature and chemical potential. The magnetic (transverse) gluon shows little sensitivity to temperature, and exhibits a mild enhancement at intermediate mu before becoming suppressed at large mu.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Potential Energy Landscapes for the 2D XY Model: Minima, Transition States and Pathways

We describe a numerical study of the potential energy landscape for the two-dimensional XY model (with no disorder), considering up to 100 spins and CPU and GPU implementations of local optimization, focusing on minima and saddles of index one (transition states). We examine both periodic and anti-periodic boundary conditions, and show that the number of stationary points located increases exponentially with increasing lattice size. The corresponding disconnectivity graphs exhibit funneled landscapes; the global minima are readily located because they exhibit relatively large basins of attraction compared to the higher energy minima as the lattice size increases.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Spectral functions of charmonium from 2 flavour anisotropic lattice data

The spectral functions of QCD can give us insight into properties of hadrons, and they are useful in probing the QCD vacuum. I will discuss the correlators and spectral functions of charmonium in high temperature two flavour QCD. The spectral functions have been obtained using the Maximum Entropy Method from anisotropic lattice data using the conserved vector current. This work has been done as part of the FASTSUM collaboration. We find that the spectral functions for zero momentum are stable. At non-zero momentum the spectral functions are less stable but still produce resonance and transport peaks. This work is part of our programme to calculate the heavy quark diffusion constant.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Tumbling through a landscape: Evidence of instabilities in high-dimensional moduli spaces

We argue that a generic instability afflicts vacua that arise in theories whose moduli space has large dimension. Specifically, by studying theories with multiple scalar fields we provide numerical evidence that for a generic local minimum of the potential the usual semiclassical bubble nucleation rate, Gamma = A e^{-B}, increases rapidly as function of the number of fields in the theory. As a consequence, the fraction of vacua with tunneling rates low enough to maintain metastability appears to fall exponentially as a function of the moduli space dimension. We discuss possible implications for the landscape of string theory. Notably, if our results prove applicable to string theory, the landscape of metastable vacua may not contain sufficient diversity to offer a natural explanation of dark energy.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Energy landscape analysis of the two-dimensional nearest-neighbor ϕ^4 model

The stationary points of the potential energy function of the ϕ^4 model on a two-dimensional square lattice with nearest-neighbor interactions are studied by means of two numerical methods: a numerical homotopy continuation method and a globally-convergent Newton-Raphson method. We analyze the properties of the stationary points, in particular with respect to a number of quantities that have been conjectured to display signatures of the thermodynamic phase transition of the model. Although no such signatures are found for the nearest-neighbor ϕ^4 model, our study illustrates the strengths and weaknesses of the numerical methods employed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Exploring the energy landscape of XY models

We investigate the energy landscape of two- and three-dimensional XY models with nearest-neighbor interactions by analytically constructing several classes of stationary points of the Hamiltonian. These classes are analyzed, in particular with respect to possible signatures of the thermodynamic phase transitions of the models. We find that, even after explicitly breaking the global O(2) symmetry of the XY spins, an exponentially large class of stationary points are singular and occur in continuous one-parameter families. This property may complicate the use of theoretical tools developed for the investigation of phase transitions based on stationary points of the energy landscape, and we discuss strategies to avoid these difficulties.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Finding all flux vacua in an explicit example

We explicitly construct all supersymmetric flux vacua of a particular Calabi-Yau compactification of type IIB string theory for a small number of flux carrying cycles and a given D3-brane tadpole. The analysis is performed in the large complex structure region by using the polynomial homotopy continuation method, which allows to find all stationary points of the polynomial equations that characterize the supersymmetric vacuum solutions. The number of vacua as a function of the D3 tadpole is in agreement with statistical studies in the literature. We calculate the available tuning of the cosmological constant from fluxes and extrapolate to scenarios with a larger number of flux carrying cycles. We also verify the range of scales for the moduli and gravitino masses recently found for a single explicit flux choice giving a Kähler uplifted de Sitter vacuum in the same construction.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Investigating the sign problem for two-dimensional $\mathcal{N}=(2,2)$ and $\mathcal{N}=(8,8)$ lattice super Yang--Mills theories

Recently there has been some controversy in the literature concerning the existence of a fermion sign problem in the $\mathcal{N}=(2,2)$ supersymmetric Yang--Mills (SYM) theories on the lattice. In this work, we address this issue by conducting Monte Carlo simulations not only for $\mathcal{N}=(2,2)$ but also for $\mathcal{N}=(8,8)$ SYM in two dimensions for the U(N) theories with N=2, using the new ideas derived from topological twisting followed by geometric discretization. Our results from simulations provide the evidence that these theories do {\it not} suffer from a sign problem as the continuum limit is approached. These results thus boost confidence that these new lattice formulations can be used successfully to explore the nonperturbative aspects of the four-dimensional $\mathcal{N}=4$ SYM theory.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Minimizing Higgs Potentials via Numerical Polynomial Homotopy Continuation

The study of models with extended Higgs sectors requires to minimize the corresponding Higgs potentials, which is in general very difficult. Here, we apply a recently developed method, called numerical polynomial homotopy continuation (NPHC), which guarantees to find all the stationary points of the Higgs potentials with polynomial-like nonlinearity. The detection of all stationary points reveals the structure of the potential with maxima, metastable minima, saddle points besides the global minimum. We apply the NPHC method to the most general Higgs potential having two complex Higgs-boson doublets and up to five real Higgs-boson singlets. Moreover the method is applicable to even more involved potentials. Hence the NPHC method allows to go far beyond the limits of the Gröbner basis approach.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Non-abelian gauged NJL models on the lattice

We use Monte Carlo simulation to probe the phase structure of a SU(2) gauge theory containing $N_f$ Dirac fermion flavors transforming in the fundamental representation of the group and interacting through an additional four fermion term. Pairs of physical flavors are implemented using the two tastes present in a reduced staggered fermion formulation of the theory. The resultant lattice theory is invariant under a set of shift symmetries which correspond to a discrete subgroup of the continuum chiral-flavor symmetry. The pseudoreal character of the representation guarantees that the theory has no sign problem. For the case of $N_f=4$ we observe a crossover in the behavior of the chiral condensate for strong four fermi coupling associated with the generation of a dynamical mass for the fermions. At weak gauge coupling this crossover is consistent with the usual continuous phase transition seen in the pure (ungauged) NJL model. However, if the gauge coupling is strong enough to cause confinement we observe a much more rapid crossover in the chiral condensate consistent with a first order phase transition

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Numerical Algebraic Geometry: A New Perspective on String and Gauge Theories

The interplay rich between algebraic geometry and string and gauge theories has recently been immensely aided by advances in computational algebra. However, these symbolic (Gröbner) methods are severely limited by algorithmic issues such as exponential space complexity and being highly sequential. In this paper, we introduce a novel paradigm of numerical algebraic geometry which in a plethora of situations overcomes these short-comings. Its so-called 'embarrassing parallelizability' allows us to solve many problems and extract physical information which elude the symbolic methods. We describe the method and then use it to solve various problems arising from physics which could not be otherwise solved.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Phase Structure of Lattice N=4 Super Yang-Mills

We make a first study of the phase diagram of four-dimensional N=4 super Yang-Mills theory regulated on a space-time lattice. The lattice formulation we employ is both gauge invariant and retains at all lattice spacings one exactly preserved supersymmetry charge. Our numerical results are consistent with the existence of a single deconfined phase at all observed values of the bare coupling.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Phase structure of QC2D at high temperature and density

We study two-color QCD with two flavors of Wilson fermion as a function of quark chemical potential mu and temperature T. We find evidence of a superfluid phase at intermediate mu and low T where the quark number density and diquark condensate are both very well described by a Fermi sphere of nearly-free quarks disrupted by a BCS condensate. This gives way to a region of deconfined quark matter at higher T and mu, with the deconfinement temperature decreasing only very slowly with increasing chemical potential. We find that heavy quarkonium bound states persist in the S-wave channels at all T and mu, with an energy reflecting the phase structure. P-wave states appear not to survive in the quarkyonic region.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Finding All the Stationary Points of a Potential Energy Landscape via Numerical Polynomial Homotopy Continuation Method

The stationary points (SPs) of a potential energy landscape play a crucial role in understanding many of the physical or chemical properties of a given system. Unless they are found analytically, there is, however, no efficient method to obtain 'all' the SPs of a given potential. We introduce a novel method, called the numerical polynomial homotopy continuation (NPHC) method, which numerically finds all the SPs, and is 'embarrassingly parallelizable'. The method requires the non-linearity of the potential to be polynomial-like, which is the case for almost all of the potentials arising in physical and chemical systems. We also certify the numerically obtained SPs so that they are independent of the numerical tolerance used during the computation. It is then straightforward to separate out the local and global minima. As a first application, we take the XY model with power-law interaction which is shown to have a polynomial-like non-linearity and apply the method.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Numerical Polynomial Homotopy Continuation Method and String Vacua

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 article, 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 \textit{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 \textit{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 and compare the results with the available results from the symbolic methods. We then show that all the vacua of a more complicated model of M theory compactified on the coset $\frac{SU(3)\times U(1)}{U(1)\times U(1)}$, which has an SU(3) structure, can be obtained by the NPHC method using a desktop machine in just about one 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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the sign problem in 2D lattice super Yang--Mills

In recent years a new class of supersymmetric lattice theories have been proposed which retain one or more exact supersymmetries for non-zero lattice spacing. Recently there has been some controversy in the literature concerning whether these theories suffer from a sign problem. In this paper we address this issue by conducting simulations of the N=(2, 2) and N=(8, 8) supersymmetric Yang--Mills theories in two dimensions for the U(N) theories with N=2,3,4, using the new twisted lattice formulations. Our results provide evidence that these theories do not suffer from a sign problem in the continuum limit. These results thus boost confidence that the new lattice formulations can be used successfully to explore non-perturbative aspects of four-dimensional N=4 supersymmetric Yang--Mills theory.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Phase transitions detached from stationary points of the energy landscape

The stationary points of the potential energy function V are studied for the ϕ^4 model on a two-dimensional square lattice with nearest-neighbor interactions. On the basis of analytical and numerical results, we explore the relation of stationary points to the occurrence of thermodynamic phase transitions. We find that the phase transition potential energy of the ϕ^4 model does in general not coincide with the potential energy of any of the stationary points of V. This disproves earlier, allegedly rigorous, claims in the literature on necessary conditions for the existence of phase transitions. Moreover, we find evidence that the indices of stationary points scale extensively with the system size, and therefore the index density can be used to characterize features of the energy landscape in the infinite-system limit. We conclude that the finite-system stationary points provide one possible mechanism of how a phase transition can arise, but not the only one.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Supersymmetric gauge theories on the lattice: Pfaffian phases and the Neuberger 0/0 problem

Recently a class of supersymmetric gauge theories have been successfully implemented on the lattice. However, there has been an ongoing debate on whether lattice versions of some of these theories suffer from a sign problem, with independent simulations for the ${\cal N} = (2, 2)$ supersymmetric Yang-Mills theories in two dimensions yielding seemingly contradictory results. Here, we address this issue from an interesting theoretical point of view. We conjecture that the sign problem observed in some of the simulations is related to the so called Neuberger 0/0 problem, which arises in ordinary non-supersymmetric lattice gauge theories, and prevents the realization of Becchi-Rouet-Stora-Tyutin symmetry on the lattice. After discussing why we expect a sign problem in certain classes of supersymmetric lattice gauge theories far from the continuum limit, we argue that these problems can be evaded by use of a non-compact parametrization of the gauge link fields.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Towards Modeling HIV Long Term Behavior

The precise mechanism that causes HIV infection to progress to AIDS is still unknown. This paper presents a mathematical model which is able to predict the entire trajectory of the HIV/AIDS dynamics, then a possible explanation for this progression is examined. A dynamical analysis of this model reveals a set of parameters which may produce two real equilibria in the model. One equilibrium is stable and represents those individuals who have been living with HIV for at least 7 to 9 years, and do not develop AIDS. The other one is unstable and represents those patients who developed AIDS in an average period of 10 years. However, further work is needed since the proposed model is sensitive to parameter variations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Stability and Bifurcation Analysis of Coupled Fitzhugh-Nagumo Oscillators

Neurons are the central biological objects in understanding how the brain works. The famous Hodgkin-Huxley model, which describes how action potentials of a neuron are initiated and propagated, consists of four coupled nonlinear differential equations. Because these equations are difficult to deal with, there also exist several simplified models, of which many exhibit polynomial-like non-linearity. Examples of such models are the Fitzhugh-Nagumo (FHN) model, the Hindmarsh-Rose (HR) model, the Morris-Lecar (ML) model and the Izhikevich model. In this work, we first prescribe the biologically relevant parameter ranges for the FHN model and subsequently study the dynamical behaviour of coupled neurons on small networks of two or three nodes. To do this, we use a computational real algebraic geometry method called the Discriminant Variety (DV) method to perform the stability and bifurcation analysis of these small networks. A time series analysis of the FHN model can be found elsewhere in related work[15].

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Stationary point analysis of the one-dimensional lattice Landau gauge fixing functional, aka random phase XY Hamiltonian

We study the stationary points of what is known as the lattice Landau gauge fixing functional in one-dimensional compact U(1) lattice gauge theory, or as the Hamiltonian of the one-dimensional random phase XY model in statistical physics. An analytic solution of all stationary points is derived for lattices with an odd number of lattice sites and periodic boundary conditions. In the context of lattice gauge theory, these stationary points and their indices are used to compute the gauge fixing partition function, making reference in particular to the Neuberger problem. Interpreted as stationary points of the one-dimensional XY Hamiltonian, the solutions and their Hessian determinants allow us to evaluate a criterion which makes predictions on the existence of phase transitions and the corresponding critical energies in the thermodynamic limit.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Lattice Landau Gauge and Algebraic Geometry

Finding the global minimum of a multivariate function efficiently is a fundamental yet difficult problem in many branches of theoretical physics and chemistry. However, we observe that there are many physical systems for which the extremizing equations have polynomial-like non-linearity. This allows the use of Algebraic Geometry techniques to solve these equations completely. The global minimum can then straightforwardly be found by the second derivative test. As a warm-up example, here we study lattice Landau gauge for compact U(1) and propose two methods to solve the corresponding gauge-fixing equations. In a first step, we obtain all Gribov copies on one and two dimensional lattices. For simple 3x3 systems their number can already be of the order of thousands. We anticipate that the computational and numerical algebraic geometry methods employed have far-reaching implications beyond the simple but illustrating examples discussed here.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Modified Lattice Landau Gauge

We propose a modified lattice Landau gauge based on stereographically projecting the link variables on the circle S^1 -> R for compact U(1) or the 3-sphere S^3 -> R^3 for SU(2) before imposing the Landau gauge condition. This can reduce the number of Gribov copies exponentially and solves the Gribov problem in compact U(1) where it is a lattice artifact. Applied to the maximal Abelian subgroup this might be just enough to avoid the perfect cancellation amongst the Gribov copies in a lattice BRST formulation for SU(N), and thus to avoid the Neuberger 0/0 problem. The continuum limit of the Landau gauge remains unchanged.

0 citations0 reviewsNo code linkedOpen access