Abstract Classical molecular dynamics simulations are based on solving Newton’s equations of motion. Using a small timestep, numerical integrators such as Verlet generate trajectories of particles as solutions to Newton’s equations. We introduce operators derived using recurrent neural networks that accurately solve Newton’s equations utilizing sequences of past trajectory data, and produce energy-conserving dynamics of particles using timesteps up to 4000 times larger compared to the Verlet timestep. We demonstrate significant speedup in many example problems including 3D systems of up to 16 particles.
more »
« less
SOLVING DIFFERENCE EQUATIONS IN SEQUENCES: UNIVERSALITY AND UNDECIDABILITY
We study solutions of difference equations in the rings of sequences and, more generally, solutions of equations with a monoid action in the ring of sequences indexed by the monoid. This framework includes, for example, difference equations on grids (for example, standard difference schemes) and difference equations in functions on words. On the universality side, we prove a version of strong Nullstellensatz for such difference equations under the assumption that the cardinality of the ground field is greater than the cardinality of the monoid and construct an example showing that this assumption cannot be omitted. On the undecidability side, we show that the following problems are undecidable:
more »
« less
- PAR ID:
- 10169850
- Date Published:
- Journal Name:
- Forum of Mathematics, Sigma
- Volume:
- 8
- ISSN:
- 2050-5094
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
This paper develops a finite approximation approach to find a non-smooth solution of an integral equation of the second kind. The equation solutions with non-smooth kernel having a non-smooth solution have never been studied before. Such equations arise frequently when modeling stochastic systems. We construct a Banach space of (right-continuous) distribution functions and reformulate the problem into an operator equation. We provide general necessary and sufficient conditions that allow us to show convergence of the approximation approach developed in this paper. We then provide two specific choices of approximation sequences and show that the properties of these sequences are sufficient to generate approximate equation solutions that converge to the true solution assuming solution uniqueness and some additional mild regularity conditions. Our analysis is performed under the supremum norm, allowing wider applicability of our results. Worst-case error bounds are also available from solving a linear program. We demonstrate the viability and computational performance of our approach by constructing three examples. The solution of the first example can be constructed manually but demonstrates the correctness and convergence of our approach. The second application example involves stationary distribution equations of a stochastic model and demonstrates the dramatic improvement our approach provides over the use of computer simulation. The third example solves a problem involving an everywhere nondifferentiable function for which no closed-form solution is available.more » « less
-
Recent work has reemphasized the importance of cardinality estimates for query optimization. While new techniques have continuously improved in accuracy over time, they still generally allow for under-estimates which often lead optimizers to make overly optimistic decisions. This can be very costly for expensive queries. An alternative approach to estimation is cardinality bounding, also called pessimistic cardinality estimation, where the cardinality estimator provides guaranteed upper bounds of the true cardinality. By never underestimating, this approach allows the optimizer to avoid potentially inefficient plans. However, existing pessimistic cardinality estimators are not yet practical: they use very limited statistics on the data, and cannot handle predicates. In this paper, we introduce SafeBound, the first practical system for generating cardinality bounds. SafeBound builds on a recent theoretical work that uses degree sequences on join attributes to compute cardinality bounds, extends this framework with predicates, introduces a practical compression method for the degree sequences, and implements an efficient inference algorithm. Across four workloads, SafeBound achieves up to 80% lower end-to-end runtimes than PostgreSQL, and is on par or better than state of the art ML-based estimators and pessimistic cardinality estimators, by improving the runtime of the expensive queries. It also saves up to 500x in query planning time, and uses up to 6.8x less space compared to state of the art cardinality estimation methods.more » « less
-
We consider the KZ differential equations over C in the case, when the hypergeometric solutions are one-dimensional integrals. We also consider the same differential equations over a finite field F_p. We study the polynomial solutions of these differential equations over F_p, constructed in a previous work joint with V. Schechtman and called the F_p-hypergeometric solutions. The dimension of the space of F_p-hypergeometric solutions depends on the prime number p. We say that the KZ equations have ample reduction for a prime p, if the dimension of the space of F_p-hypergeometric solutions is maximal possible, that is, equal to the dimension of the space of solutions of the corresponding KZ equations over C. Under the assumption of ample reduction, we prove a determinant formula for the matrix of coordinates of basis F_p-hypergeometric solutions. The formula is analogous to the corresponding formula for the determinant of the matrix of coordinates of basis complex hypergeometric solutions, in which binomials (z_i−z_j)^{M_i+M_j} are replaced with (z_i−z_j)^{Mi+Mj−p} and the Euler gamma function Γ(x) is replaced with a suitable F_p-analog defined on F_pmore » « less
-
Abstract The goal of this paper is to show that valuation theory and Hopf theory are compatible on the class of generalized permutahedra. We prove that the Hopf structure $$\textbf {GP}^+$$ on these polyhedra descends, modulo the inclusion-exclusion relations, to an indicator Hopf monoid $$\mathbb {I}(\textbf {GP}^+)$$ of generalized permutahedra that is isomorphic to the Hopf monoid of weighted ordered set partitions. This quotient Hopf monoid $$\mathbb {I}(\textbf {GP}^+)$$ is cofree. It is the terminal object in the category of Hopf monoids with polynomial characters; this partially explains the ubiquity of generalized permutahedra in the theory of Hopf monoids. This Hopf theoretic framework offers a simple, unified explanation for many new and old valuations on generalized permutahedra and their subfamilies. Examples include, for matroids: the Chern–Schwartz–MacPherson cycles, Eur’s volume polynomial, the Kazhdan–Lusztig polynomial, the motivic zeta function, and the Derksen–Fink invariant; for posets: the order polynomial, Poincaré polynomial, and poset Tutte polynomial; for generalized permutahedra: the universal Tutte character and the corresponding class in the Chow ring of the permutahedral variety. We obtain several algebraic and combinatorial corollaries; for example, the existence of the valuative character group of $$\textbf {GP}^+$$ and the indecomposability of a nestohedron into smaller nestohedra.more » « less