Search

We introduce a new version v1.1 of NeatIBP. In this new version, a Kira interface is included. It allows the user to reduce the integration-by-parts (IBP) identity systems generated by NeatIBP using Kira in a highly automated way. This new version also implements the so-called spanning cuts method. It helps to reduce the total computational complexity of IBP reduction for certain hard problems. Another important feature of this new version is an algorithm to simplify the solution module of the syzygy equations hinted by the idea of maximal cuts., Comment: 57pages, 5 tables, 7 figures

We give illustrative examples of how the computer algebra system OSCAR can support research in commutative algebra and algebraic geometry. We start with a thorough introduction to Groebner basis techniques, with particular emphasis on the computation of syzygies, then apply these techniques to deal with ideal and ring theoretic concepts such as primary decomposition and normalization, and finally use them for geometric case studies which concern curves and surfaces, both from a local and global point of view., Comment: 50 pages, 10 figures

Our recent extension of Arnold's classification includes all singularities of corank up to two equivalent to a germ with a non-degenerate Newton boundary, thus broadening the classification's scope significantly by a class which is unbounded with respect to modality and Milnor number. This method is based on proving that all right-equivalence classes within a mu-constant stratum can be represented by a single normal form derived from a regular basis of a suitably selected special fiber. While both Arnold's and our preceding work on normal forms addresses the determination of a normal form family containing the given germ, this paper takes the next natural step: We present an algorithm for computing for a given germ the values of the moduli parameters in its normal form family, that is, a normal form equation in its stable equivalence class. This algorithm will be crucial for understanding the moduli stacks of such singularities. The implementation of this algorithm, along with the foundational classification techniques, is implemented in the library arnold.lib for the computer algebra system Singular., Comment: 28 pages

Computations over the rational numbers often encounter the problem of intermediate coefficient growth. A solution to this is provided by modular methods, which apply the algorithm under consideration modulo a number of primes and then lift the results to the rationals. We present a novel, massively parallel framework for modular computations with polynomial data, which is able to cover a broad spectrum of applications in commutative algebra and algebraic geometry. We demonstrate the framework's effectiveness in Groebner basis computations over the rationals and algorithmic methods from birational geometry. In particular, we develop algorithms to compute images and domains of rational maps, as well as determining invertibility and computing inverses. Our implementation is based on the Singular/GPI-Space framework, which uses the computer algebra system Singular as computational backend, while coordination and communication of parallel computations is handled by the workflow management system GPI-Space, which relies on Petri nets as its mathematical modeling language. Convenient installation is realized through the package manager Spack. Relying on Petri nets, our approach provides automated parallelization and balancing of the load between computation, lifting, stabilization testing, and potential verification. We use error tolerant rational reconstruction to ensure termination as long as for a fixed computation there exist only finitely many bad primes. Via stabilization testing, our approach automatically finds with high probablity a minimal set of primes required for the successful reconstruction. We present timings to illustrate the potential for a game changing improvement of performance over previous modular and non-modular methods. In particular, we illustrate that the approach scales very well with the number of processor cores used for the computation., Comment: 56 pages, 10 figures, 6 tables

We present an enhanced algorithm for exploring mirror symmetry for elliptic curves through the correspondence of algebraic and tropical geometry, focusing on Gromov-Witten invariants of elliptic curves and, in particular, Hurwitz numbers. We present a new highly efficient algorithm for computing generating series for these numbers. We have implemented the algorithm both using Singular and OSCAR. The implementations outperform by far the current method provided in Singular. The OSCAR implementation, benefiting in particular from just-in-time compilation, again by far outperforms the implementation of the new algorithm in Singular. This advancement in computing the Gromov-Witten invariants facilitates a study of number theoretic and geometric properties of the generating series, including quasi-modularity and homogeneity., Comment: 20 pages, 5 figures, 2 tables. arXiv admin note: text overlap with arXiv:1309.5893

In this work, we present the package {\sc NeatIBP}, which automatically generates small-size integration-by-parts (IBP) identities for Feynman integrals. Based on the syzygy and module intersection techniques, the generated IBP identities' propagator degree is controlled and thus the size of the system of IBP identities is shorter than that generated by the standard Laporta algorithm. This package is powered by the computer algebra systems {\sc Mathematica} and {\sc Singular}, and the library {\sc SpaSM}. It is parallelized on the level of Feynman integral sectors. The generated small-size IBP identities can subsequently be used for either finite field reduction or analytic reduction. We demonstrate the capabilities of this package on several multi-loop IBP examples., Comment: update to published version

We uncover a connection between two seemingly separate subjects in integrable models: the representation theory of the affine Temperley-Lieb algebra, and the algebraic structure of solutions to the Bethe equations of the XXZ spin chain. We study the solution of Bethe equations analytically by computational algebraic geometry, and find that the solution space encodes rich information about the representation theory of Temperley-Lieb algebra. Using these connections, we compute the partition function of the completely-packed loop model and of the closely related random-cluster Potts model, on medium-size lattices with toroidal boundary conditions, by two quite different methods. We consider the partial thermodynamic limit of infinitely long tori and analyze the corresponding condensation curves of the zeros of the partition functions. Two components of these curves are obtained analytically in the full thermodynamic limit., Comment: 66 pages

Multivariate partial fractioning is a powerful tool for simplifying rational function coefficients in scattering amplitude computations. Since current research problems lead to large sets of complicated rational functions, performance of the partial fractioning as well as size of the obtained expressions are a prime concern. We develop a large scale parallel framework for multivariate partial fractioning, which implements and combines an improved version of Leinartas' algorithm and the {\sc MultivariateApart} algorithm. Our approach relies only on open source software. It combines parallelism over the different rational function coefficients with parallelism for individual expressions. The implementation is based on the \textsc{Singular}/\textsc{GPI-Space framework} for massively parallel computer algebra, which formulates parallel algorithms in terms of Petri nets. The modular nature of this approach allows for easy incorporation of future algorithmic developments into our package. We demonstrate the performance of our framework by simplifying expressions arising from current multiloop scattering amplitude problems., Comment: Change the title and the paper format, for the journal Computer Physics Communications

In his groundbreaking work on classification of singularities with regard to right and stable equivalence of germs, Arnold has listed normal forms for all isolated hypersurface singularities over the complex numbers with either modality less than or equal to two or Milnor number less than or equal to 16. Moreover, he has described an algorithmic classifier, which determines the type of a given such singularity. In the present paper, we extend Arnold's work to a large class of singularities which is unbounded with regard to modality and Milnor number. We develop an algorithmic classifier, which determines a normal form for any singularity with corank less than or equal to two which is equivalent to a germ with non-degenerate Newton boundary in the sense of Kouchnirenko. In order to realize the classifier, we prove a normal form theorem: Suppose K is a mu-constant stratum of the jet space which contains a germ with a non-degenerate Newton boundary. We first observe that all germs in K are equivalent to some germ with the same fixed non-degenerate Newton boundary. We then prove that all right-equivalence classes of germs in K can be covered by a single normal form obtained from a regular basis of an appropriately chosen special fiber. All algorithms are implemented in the library arnold.lib for the computer algebra system Singular., Comment: 32 pages

In this manuscript, which is to appear in the proceedings of the conference "MathemAmplitude 2019" in Padova, Italy, we provide an overview of the module intersection method for the the integration-by-parts (IBP) reduction of multi-loop Feynman integrals. The module intersection method, based on computational algebraic geometry, is a highly efficient way of getting IBP relations without double propagator or with a bound on the highest propagator degree. In this manner, trimmed IBP systems which are much shorter than the traditional ones can be obtained. We apply the modern, Petri net based, workflow management system GPI-Space in combination with the computer algebra system Singular to solve the trimmed IBP system via interpolation and efficient parallelization. We show, in particular, how to use the new plugin feature of GPI-Space to manage a global state of the computation and to efficiently handle mutable data. Moreover, a Mathematica interface to generate IBPs with restricted propagator degree, which is based on module intersection, is presented in this review., Comment: 19 pages, 5 figures. To appear in the proceedings of "MathemAmplitudes 2019: Intersection Theory & Feynman Integrals", held in Padova 18-20 December 2019

We present an efficient method to shorten the analytic integration-by-parts (IBP) reduction coefficients of multi-loop Feynman integrals. For our approach, we develop an improved version of Leinartas' multivariate partial fraction algorithm, and provide a modern implementation based on the computer algebra system Singular. Furthermore, We observe that for an integral basis with uniform transcendental (UT) weights, the denominators of IBP reduction coefficients with respect to the UT basis are either symbol letters or polynomials purely in the spacetime dimension $D$. With a UT basis, the partial fraction algorithm is more efficient both with respect to its performance and the size reduction. We show that in complicated examples with existence of a UT basis, the IBP reduction coefficients size can be reduced by a factor of as large as $\sim 100$. We observe that our algorithm also works well for settings without a UT basis., Comment: minor changes, typos corrected

In this article, we present a massively parallel framework for computing tropicalizations of algebraic varieties which can make use of finite symmetries. We compute the tropical Grassmannian TGr$_0(3,8)$, and show that it refines the $15$-dimensional skeleton of the Dressian Dr$(3,8)$ with the exception of $23$ special cones for which we construct explicit obstructions to the realizability of their tropical linear spaces. Moreover, we propose algorithms for identifying maximal-dimensional tropical cones which belong to the positive tropicalization. These algorithms exploit symmetries of the tropical variety even though the positive tropicalization need not be symmetric. We compute the maximal-dimensional cones of the positive Grassmannian TGr$^+(3,8)$ and compare them to the cluster complex of the classical Grassmannian Gr$(3,8)$., Comment: 32 pages, 9 figures

The behavior of a certain random growth process is analyzed on arbitrary regular and non-regular graphs. Our argument is based on the Expander Mixing Lemma, which entails that the results are strongest for Ramanujan graphs, which asymptotically maximize the spectral gap. Further, we consider Erd\H{o}s--R\'enyi random graphs and compare our theoretical results with computational experiments on flip graphs of point configurations. The latter is relevant for enumerating triangulations., Comment: 22 pages, 7 figures, 1 table. This version makes several changes based on feedback of the first version. This includes a change to the title and a new section with results on Erd\H{o}s--R\'{e}nyi random graphs

We introduce an algebro-geometrically motived integration-by-parts (IBP) reduction method for multi-loop and multi-scale Feynman integrals, using a framework for massively parallel computations in computer algebra. This framework combines the computer algebra system Singular with the workflow management system GPI-Space, which is being developed at the Fraunhofer Institute for Industrial Mathematics (ITWM). In our approach, the IBP relations are first trimmed by modern algebraic geometry tools and then solved by sparse linear algebra and our new interpolation methods. These steps are efficiently automatized and automatically parallelized by modeling the algorithm in GPI-Space using the language of Petri-nets. We demonstrate the potential of our method at the nontrivial example of reducing two-loop five-point nonplanar double-pentagon integrals. We also use GPI-Space to convert the basis of IBP reductions, and discuss the possible simplification of IBP coefficients in a uniformly transcendental basis., Comment: minor corrections; We encourage researchers in high energy community to send us IBP reduction problems (mailto: alessandro.georgoudis@physics.uu.se)

We study generating series of Gromov-Witten invariants of $E\times\mathbb{P}^1$ and their tropical counterparts. Using tropical degeneration and floor diagram techniques, we can express the generating series as sums of Feynman integrals, where each summand corresponds to a certain type of graph which we call a pearl chain. The individual summands are --- just as in the case of mirror symmetry of elliptic curves, where the generating series of Hurwitz numbers equals a sum of Feynman integrals --- complex analytic path integrals involving a product of propagators (equal to the Weierstrass-$\wp$-function plus an Eisenstein series). We also use pearl chains to study generating functions of counts of tropical curves in $E_{\mathbb{T}}\times\mathbb{P}^1_\mathbb{T}$ of so-called leaky degree., Comment: 26 pages, 12 figures. We thank an anonymous referee for pointing out a mistake which we fixed in this version

In this note for the joint meeting of DMV and GDM we illustrate with examples the role of computer algebra in university mathematics education. We discuss its potential in teaching algebra, but also computer algebra as a subject in its own right, its value in the context of practical programming projects and its role as a research topic in student papers., Comment: German, 5 pages

The design and implementation of parallel algorithms is a fundamental task in computer algebra. Combining the computer algebra system Singular and the workflow management system GPI-Space, we have developed an infrastructure for massively parallel computations in commutative algebra and algebraic geometry. In this note, we give an overview on the current capabilities of this framework by looking into three sample applications: determining smoothness of algebraic varieties, computing GIT-fans in geometric invariant theory, and determining tropicalizations. These applications employ algorithmic methods originating from commutative algebra, sheaf structures on manifolds, local geometry, convex geometry, group theory, and combinatorics, illustrating the potential of the framework in further problems in computer algebra., Comment: 9 pages, 6 figures

We prove a tropical mirror symmetry theorem for descendant Gromov-Witten invariants of the elliptic curve, generalizing the tropical mirror symmetry theorem for Hurwitz numbers of the elliptic curve, Theorem 2.20 in [B\"ohm J., Bringmann K., Buchholz A., Markwig H., J. Reine Angew. Math. 732 (2017), 211-246, arXiv:1309.5893]. For the case of the elliptic curve, the tropical version of mirror symmetry holds on a fine level and easily implies the equality of the generating series of descendant Gromov-Witten invariants of the elliptic curve to Feynman integrals. To prove tropical mirror symmetry for elliptic curves, we investigate the bijection between graph covers and sets of monomials contributing to a coefficient in a Feynman integral. We also soup up the traditional approach in mathematical physics to mirror symmetry for the elliptic curve, involving operators on a Fock space, to give a proof of tropical mirror symmetry for Hurwitz numbers of the elliptic curve. In this way, we shed light on the intimate relation between the operator approach on a bosonic Fock space and the tropical approach.

Introducing parallelism and exploring its use is still a fundamental challenge for the computer algebra community. In high performance numerical simulation, on the other hand, transparent environments for distributed computing which follow the principle of separating coordination and computation have been a success story for many years. In this paper, we explore the potential of using this principle in the context of computer algebra. More precisely, we combine two well-established systems: The mathematics we are interested in is implemented in the computer algebra system Singular, whose focus is on polynomial computations, while the coordination is left to the workflow management system GPI-Space, which relies on Petri nets as its mathematical modeling language, and has been successfully used for coordinating the parallel execution (autoparallelization) of academic codes as well as for commercial software in application areas such as seismic data processing. The result of our efforts is a major step towards a framework for massively parallel computations in the application areas of Singular, specifically in commutative algebra and algebraic geometry. As a first test case for this framework, we have modeled and implemented a hybrid smoothness test for algebraic varieties which combines ideas from Hironaka's celebrated desingularization proof with the classical Jacobian criterion. Applying our implementation to two examples originating from current research in algebraic geometry, one of which cannot be handled by other means, we illustrate the behavior of the smoothness test within our framework, and investigate how the computations scale up to 256 cores., Comment: 39 pages, 5 figures, 2 tables

We present the powerful module-intersection integration-by-parts (IBP) method, suitable for multi-loop and multi-scale Feynman integral reduction. Utilizing modern computational algebraic geometry techniques, this new method successfully trims traditional IBP systems dramatically to much simpler integral-relation systems on unitarity cuts. We demonstrate the power of this method by explicitly carrying out the complete analytic reduction of two-loop five-point non-planar hexagon-box integrals, with degree-four numerators, to a basis of $73$ master integrals., Comment: minor corrections, updated references

Integration-by-parts identities between loop integrals arise from the vanishing integration of total derivatives in dimensional regularization. Generic choices of total derivatives in the Baikov or parametric representations lead to identities which involve dimension shifts. These dimension shifts can be avoided by imposing a certain constraint on the total derivatives. The solutions of this constraint turn out to be a specific type of syzygies which correspond to logarithmic vector fields along the Gram determinant formed of the independent external and loop momenta. We present an explicit generating set of solutions in Baikov representation, valid for any number of loops and external momenta, obtained from the Laplace expansion of the Gram determinant. We provide a rigorous mathematical proof that this set of solutions is complete. This proof relates the logarithmic vector fields in question to ideals of submaximal minors of the Gram matrix and makes use of classical resolutions of such ideals., Comment: 12 pages, two-column format, 3 figures. Added example with internal masses; minor corrections; journal version

Computations over the rational numbers often suffer from intermediate coefficient swell. One solution to this problem is to apply the given algorithm modulo a number of primes and then lift the modular results to the rationals. This method is guaranteed to work if we use a sufficiently large set of good primes. In many applications, however, there is no efficient way of excluding bad primes. In this note, we describe a technique for rational reconstruction which will nevertheless return the correct result, provided the number of good primes in the selected set of primes is large enough. We give a number of illustrating examples which are implemented using the computer algebra system Singular and the programming language Julia. We discuss applications of our technique in computational algebraic geometry., Comment: 8 pages, 1 figure, 1 table

This note is based on the plenary talk given by the second author at MACIS 2015, the Sixth International Conference on Mathematical Aspects of Computer and Information Sciences. Motivated by some of the work done within the Priority Programme SPP 1489 of the German Research Council DFG, we discuss a number of current challenges in the development of Open Source computer algebra systems. The main focus is on algebraic geometry and the system Singular., Comment: 18 pages, 13 figures

In (Arnold, 1985), V.I. Arnold has obtained normal forms and has developed a classifier for, in particular, all isolated hypersurface singularities over the complex numbers up to modality 2. Building on a series of 105 theorems, this classifier determines the type of the given singularity. However, for positive modality, this does not fix the right equivalence class of the singularity, since the values of the moduli parameters are not specified. In this paper, we present a simple classification algorithm for isolated hypersurface singularities of corank and modality up to two. For a singularity given by a polynomial over the rationals, the algorithm determines its right equivalence class by specifying a polynomial representative in Arnold's list of normal forms., Comment: 19 pages, 5 figures, minor revisions

We propose an algorithm to compute the GIT-fan for torus actions on affine varieties with symmetries. The algorithm combines computational techniques from commutative algebra, convex geometry and group theory. We have implemented our algorithm in the Singular library gitfan.lib. Using our implementation, we compute the Mori chamber decomposition of the cone of movable divisors of $\bar{M}_{0,6}$., Comment: 16 pages, 4 figures, 2 tables, minor revisions

Based on an idea in Hironaka's proof of resolution of singularities, we present an algorithmic smoothness test for algebraic varieties. The test is inherently parallel and does not involve the calculation of codimension-sized minors of the Jacobian matrix of the variety. We also describe a hybrid method which combines the new method with the Jacobian criterion, thus making use of the strengths of both approaches. We have implemented all algorithms in the computer algebra system Singular, and compare the different approaches with respect to timings and memory usage. The test examples originate from questions in algebraic geometry, where the use of the Jacobian criterion is impractical due to the number and size of the minors involved., Comment: 12 pages, 4 algorithms

We present a classification algorithm for isolated hypersurface singularities of corank 2 and modality 1 over the real numbers. For a singularity given by a polynomial over the rationals, the algorithm determines its right equivalence class by specifying all representatives in Arnold's list of normal forms (Arnold et al. 1985) belonging to this class, and the corresponding values of the moduli parameter. We discuss how to computationally realize the individual steps of the algorithm for all singularities in consideration, and give explicit examples. The algorithm is implemented in the Singular library realclassify.lib., Comment: 31 pages, 5 figures, 1 table, improvements in the algorithms

3D printing and other layer manufacturing processes are challenged by dimensional accuracy. Several techniques are used to validate and calibrate dimensional accuracy through the complete building envelope. The validation process involves the growing and measuring of a shape with known parameters. The measured result is compared with the intended digital model. Processes with the risk of deformation after time or post processing may find this technique beneficial. We propose to use objects from algebraic geometry as test shapes. A cubic surface is given as the zero set of a 3rd degree polynomial with 3 variables. A class of cubics in real 3D space contains exactly 27 real lines. We provide a library for the computer algebra system Singular which, from 6 given points in the plane, constructs a cubic and the lines on it. A surface shape derived from a cubic offers simplicity to the dimensional comparison process, in that the straight lines and many other features can be analytically determined and easily measured using non-digital equipment. For example, the surface contains so-called Eckardt points, in each of which three of the lines intersect, and also other intersection points of pairs of lines. Distances between these intersection points can easily be measured, since the points are connected by straight lines. At all intersection points of lines, angles can be verified. Hence, many features distributed over the build volume are known analytically, and can be used for the validation process. Due to the thin shape geometry the material required to produce an algebraic surface is minimal. This paper is the first in a series that proposes the process chain to first define a cubic with a configuration of lines in a given print volume and then to develop the point cloud for the final manufacturing. Simple measuring techniques are recommended., Comment: 8 pages, 1 figure, 1 table

We present a new algorithm for computing integral bases in algebraic function fields of one variable, or equivalently for constructing the normalization of a plane curve. Our basic strategy makes use of the concepts of localization and completion, together with the Chinese remainder theorem, to reduce the problem to the task of finding integral bases for the branches of each singularity of the curve. To solve the latter task, in turn, we work with suitably truncated Puiseux expansions. In contrast to van Hoeij's algorithm, which also relies on Puiseux expansions (but pursues a different strategy), we use Hensel's lemma as a key ingredient. This allows us at some steps of the algorithm to compute factors corresponding to conjugacy classes of Puiseux expansions, without actually computing the individual expansions. In this way, we make substantially less use of the Newton-Puiseux algorithm. In addition, our algorithm is inherently parallel. As a result, it outperforms in most cases any other algorithm known to us by far. Typical applications are the computation of adjoint ideals and, based on this, the computation of Riemann-Roch spaces and the parametrization of rational curves., Comment: 47 pages; revised structure

We present new algorithms for computing adjoint ideals of curves and thus, in the planar case, adjoint curves. With regard to terminology, we follow Gorenstein who states the adjoint condition in terms of conductors. Our main algorithm yields the Gorenstein adjoint ideal G of a given curve as the intersection of what we call local Gorenstein adjoint ideals. Since the respective local computations do not depend on each other, our approach is inherently parallel. Over the rationals, further parallelization is achieved by a modular version of the algorithm which first computes a number of the characteristic p counterparts of G and then lifts these to characteristic zero. As a key ingredient, we establish an efficient criterion to verify the correctness of the lift. Well-known applications are the computation of Riemann-Roch spaces, the construction of points in moduli spaces, and the parametrization of rational curves. We have implemented different variants of our algorithms together with Mnuk's approach in the computer algebra system Singular and give timings to compare the performance of the algorithms., Comment: 32 pages