Symbolic Computation Research Articles

Symbolic equation solving via reinforcement learning

Machine-learning methods are rapidly being adopted in a wide variety of social, economic, and scientific contexts, yet they are notorious for struggling with exact mathematics. A typical example is computer algebra, which includes tasks like simplifying mathematical terms, calculating formal derivatives, or finding exact solutions of algebraic equations. Traditional software packages for these purposes are commonly based on a huge database of rules for how a specific operation (e.g., differentiation) transforms a certain term (e.g., sine function) into another one (e.g., cosine function). These rules have usually needed to be discovered and subsequently programmed by humans. Efforts to automate this process by machine-learning approaches are faced with challenges like the singular nature of solutions to mathematical problems, when approximations are unacceptable, as well as hallucination effects leading to flawed reasoning. We propose a novel deep-learning interface involving a reinforcement-learning agent that operates a symbolic stack calculator to explore mathematical relations. By construction, this system is capable of exact transformations and immune to hallucination. Using the paradigmatic example of solving linear equations in symbolic form, we demonstrate how our reinforcement-learning agent autonomously discovers elementary transformation rules and step-by-step solutions.

Automatic differential kinematics of serial manipulator robots through dual numbers

Dual Numbers are an extension of real numbers known for its capability of performing exact automatic differentiation of one-valued functions theoretically without error approximation. Also, Differential Kinematics of robots involves the computation of the Jacobian, which is a matrix of partial derivatives of the Forward Kinematic equations with respect to the robot’s joints. Thus, to perform the automatic calculation of the Jacobian matrix, this paper presents an extension of dual numbers composed of a scalar real part and a vector dual part, where the real part represents the angular value of the robot joint, and the dual part represents the direction of the corresponding partial derivative for each joint. The presented method was implemented in Matlab through Object Orientes Programming (OOP), and the results for calculating the Jacobian of the KUKA KR 500 robot model for 1000 random postures were subsequently compared in terms of execution time and Mean Squared Error (MSE) with other conventional methods: the geometric method, the symbolic method, and the finite difference method. The results showed a significant improvement in the computing time for calculating the Jacobian of the robotic model compared to the other methods, as well as a minimum MSE having as reference the numerical value of the symbolic calculations.

Study of Fifth-Order Nonlinear Caudrey–Dodd–Gibbon Model for Multiwave, M-Shaped Solitons, Lumps, Generalized Breathers, Manifold Periodic, Rogue Wave and Kink Wave Interactions

Numerous novel wave structures, such as lump, lump with one-kink, lump with two-kink soliton, generalized breather, Ma breather, Kuznetsov-Ma breather, Akhmediev breather, manifold periodic, and rogue wave solutions to [Formula: see text]-dimensional fifth-order nonlinear Caudrey–Dodd–Gibbon (FNCDG) model via symbolic computation are studied based on the ansatz function scheme. Studying the lump-type solution involves selecting the function as the generic quadratic polynomial function. The lump with one- and two-kink solutions is created by mixing a quadratic function with one or two exponential functions. By adding hyperbolic function with a quadratic function, the rogue wave solutions are manipulated. In addition, the multiwave, [Formula: see text]-shaped rational solitons and their interactions with kink waves are analytically evaluated. Furthermore, different 3D, 2D, and contour profiles are created in different dimensions by changing the parameter values.

Fully charmed tetraquarks from LHC to FCC: natural stability from fragmentation

We investigate the inclusive production of fully charmed tetraquarks, T4c(0++)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$T_{4c}(0^{++})$$\\end{document} or T4c(2++)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$T_{4c}(2^{++})$$\\end{document} radial excitations, in high-energy proton collisions. We build our study upon the collinear fragmentation of a single parton in a variable-flavor number scheme, suited to describe the tetraquark formation mechanism from moderate to large transverse-momentum regimes. To this extent, we derive a novel set of DGLAP-evolving collinear fragmentation functions, named TQ4Q1.0 determinations. They encode initial-scale inputs corresponding to both gluon and heavy-quark fragmentation channels, defined within the context of quark-potential and spin-physics inspired models, respectively. We work within the NLL/NLO+\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\hbox {NLO}^+$$\\end{document} hybrid factorization and make use of the JETHAD numeric interface along with the symJETHAD symbolic calculation plugin. With these tools, we provide predictions for high-energy observables sensitive to T4c\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$T_{4c}$$\\end{document} plus jet emissions at center-of-mass energies ranging from 14 TeV at the LHC to the 100 TeV nominal energy of the FCC.

Symmetries and symmetry-generated averages of elastic constants up to the sixth order of nonlinearity for all crystal classes, isotropy and transverse isotropy.

Algebraic expressions for averaging linear and nonlinear stiffness tensors from general anisotropy to different effective symmetries (11 Laue classes elastically representing all 32 crystal classes, and two non-crystalline symmetries: isotropic and cylindrical) have been derived by automatic symbolic computations of the arithmetic mean over the set of rotational transforms determining a given symmetry. This approach generalizes the Voigt average to nonlinear constants and desired approximate symmetries other than isotropic, which can be useful for a description of textured polycrystals and rocks preserving some symmetry aspects. Low-symmetry averages have been used to derive averages of higher symmetry to speed up computations. Relationships between the elastic constants of each symmetry have been deduced from their corresponding averages by resolving the rank-deficient system of linear equations. Isotropy has also been considered in terms of generalized Lamé constants. The results are published in the form of appendices in the supporting information for this article and have been deposited in the Mendeley database.

Innovative approache for the nonlinear atangana conformable Klein-Gordon equation unveiling traveling wave patterns

The aim of current work is to establish novel traveling wave solutions of the nonlinear Atangana conformable Klein - Gordon equation using a new extended direct algebraic technique. The Klein - Gordon equation is the relativistic state of the Schrödinger equation with a second - order time derivative and zero spin. Complex wave variable transformation is used to convert Atangana conformable nonlinear differential equation into an ordinary differential equation. Using the proposed technique based on Maple software structure, various types of solutions, such as, generalized trigonometric, generalized hyperbolic, and exponential functions, are established. When special parameteric values are considered for this method, solitary wave solutions can be obtained through other methods, such as the (G′G)-expansion method, the modified Kudryashov method, the sub-equation method, and so forth. A physical explanation is provided for the solutions under consideration to enhance comprehension of the physical phenomena resulting from the obtained solutions, provided that the physical parameters are set appropriately using 3D, 2D, and contour simulations. The results demonstrated that the new extended direct algebraic method provides a more potent mathematical tool for solving numerous more nonlinear partial differential equations with the aid of symbolic computation.

Kink Wave Phenomena in the Nonlinear Partial Differential Equation Representing the Transmission Line Model of Microtubules for Nanoionic Currents

This paper provides several new traveling wave solutions for a nonlinear partial differential equation (PDE) by applying symbolic computation and a new approach, the Riccati–Bernoulli sub-ODE method, in a computer algebra system. Herein, employing the Bäcklund transformation, we solve a nonlinear PDE associated with nanobiosciences and biophysics based on the transmission line model of microtubules for nanoionic currents. The equation introduced here in this form is suitable for critical nanoscience concerns like cell signaling and might continue to explain some of the basic cognitive functions in neurons. We employ advanced procedures to replicate the previously detected solitary waves. We offer our solutions in graphical forms, such as 3D and contour plots, using Mathematica. We can generalize the elementary method to other nonlinear equations in physics, requiring only a few steps.

Lattice points on polyominoes of inversion sequences

Inversion sequences of length n are positive integer sequences e 1 e 2 …e n such that 1 ≤ ei ≤ i for all 1 ≤ i ≤ n. These sequences are in bijection with the permutations of [n]. This paper focuses on the polyomino or barpgraph representation of the inversion sequences. More precisely, we study the distribution of lattice points on these polyominoes. We find the generating functions respect to the length, the number of interior vertices, corners, and vertices of a given degree. We also give simple explicit formulas for the total values of these parameters over all polyominoes of inversion sequences of length n. Throughout this work, we use symbolic computer algebra to facilitate the computations.

Numerical and Symbolic Analysis for Mathematical Problem-Solving with Maple

This study explores the versatile capabilities of Maple, a widely used mathematical software, in addressing a wide range of numerical and symbolic computations essential for scientific and engineering applications. The researchers investigated Maple's diverse suite of tools, including numerical integration, nonlinear equation solving, polynomial interpolation, symbolic integration, and various numerical methods. Through an in-depth literature review, illustrated case studies, and detailed performance evaluations, the paper demonstrates the effectiveness and accuracy of Maple's computational approaches in dealing with complex problems in various areas of applied mathematics. This study's findings underscored Maple's tremendous value as a reliable and comprehensive software package for researchers, scientists, and professionals involved in advanced mathematical analysis and scientific computing. Furthermore, the paper highlighted Maple's versatility in creating high-quality three-dimensional plots, crucial for visualizing and analyzing complex mathematical and scientific data. Using either sets or lists, the ability to display multiple surfaces in a single three-dimensional plot showcases Maple's power in data visualization and communicating complex ideas. By positioning Maple as a powerful platform for solving versatile mathematical problems, this study highlights the software's indispensable role in advancing scientific discoveries and engineering innovations.

Exploring Multiwave Solutions to the Integrable Combined pKP-BKP Equation in (3+1)-Dimensions

Under consideration in the current paper is a new combined Painleve integrable equation in (3+1)-dimensions, namely the potential Kadomtsev-Petviashvili equation incorporating the B-type Kadomtsev-Petviashvili equation (pKP-BKP equation). Maple symbolic calculations are made to present abundant multiwave solutions, which cover a breather-kink wave interacting with one-kink wave, and a breather-kink wave interacting with two-kink waves, as well as a breather-kink wave interacting with three-kink waves. Particularly, the dynamic and structural characteristics of some derived solutions are illustrated through some vivid 3D graphics.

Modern Tensor-Spinor Symbolic Algebra Algorithms and Computing Non-Closure Geometry & Holoraumy in 11D, 𝒩 = 1 Supergravity

Modern Tensor-Spinor Symbolic Algebra Algorithms and Computing Non-Closure Geometry & Holoraumy in 11D, 𝒩 = 1 Supergravity

Lessons on Datasets and Paradigms in Machine Learning for Symbolic Computation: A Case Study on CAD

Symbolic Computation algorithms and their implementation in computer algebra systems often contain choices which do not affect the correctness of the output but can significantly impact the resources required: such choices can benefit from having them made separately for each problem via a machine learning model. This study reports lessons on such use of machine learning in symbolic computation, in particular on the importance of analysing datasets prior to machine learning and on the different machine learning paradigms that may be utilised. We present results for a particular case study, the selection of variable ordering for cylindrical algebraic decomposition, but expect that the lessons learned are applicable to other decisions in symbolic computation. We utilise an existing dataset of examples derived from applications which was found to be imbalanced with respect to the variable ordering decision. We introduce an augmentation technique for polynomial systems problems that allows us to balance and further augment the dataset, improving the machine learning results by 28% and 38% on average, respectively. We then demonstrate how the existing machine learning methodology used for the problem—classification—might be recast into the regression paradigm. While this does not have a radical change on the performance, it does widen the scope in which the methodology can be applied to make choices.

Solitary wave type solutions of nonlinear improved mKdV equation by modified techniques

In this paper, the improved and modified version of the Sardar sub-equation method (IMSSEM) and the improved generalized Riccati equation mapping method (IGREMM) are manipulated effectively and generously to determine the exact solitary wave soliton solutions of the improved modified KdV (mKdV) equation. The purpose of this study is to provide novel exact solutions to the improved mKdV equation. Specifically, we utilized IMSSEM and IGREMM to study different solutions of the nonlinear improved mKdV equation, focusing on exponential, trigonometric, and trigonometric hyperbolic type solutions. Furthermore, the plotting of various solutions for direct viewing analysis is provided in two and three-dimensional graphs. The new strategies are straightforward, quick, and efficient and have many other advantages, whereas, they provide the most accurate and unique solution to many other types of nonlinear partial differential equations (NPDEs), which usually arise in engineering and applied sciences. It should be noted that these methodologies are novel mathematical instruments that have shown to be the most effective mathematical tools for solving higher-order nonlinear partial differential equations in mathematical physics. Symbolic computation was used to validate all of the solutions that were established. Thus, it is also hoped that these techniques will ultimately reduce the cumbersome workload involved during the process of solutions to complicated NLPDEs.

Discrete differential geometry-based model for nonlinear analysis of axisymmetric shells

In this paper, we propose a novel one-dimensional (1D) discrete differential geometry (DDG)-based numerical method for geometrically nonlinear mechanics analysis (e.g., buckling and snapping) of axisymmetric shell structures. Our numerical model leverages differential geometry principles to accurately capture the complex nonlinear deformation patterns exhibited by axisymmetric shells. By discretizing the axisymmetric shell into interconnected 1D elements along the meridional direction, the in-plane stretching and out-of-bending potentials are formulated based on the geometric principles of 1D nodes and edges under the Kirchhoff–Love hypothesis, and elastic force vector and associated Hessian matrix required by equations of motion are later derived based on symbolic calculation. Through extensive validation with available theoretical solutions and finite element method (FEM) simulations in literature, our model demonstrates high accuracy in predicting the nonlinear behavior of axisymmetric shells. Importantly, compared to the classical theoretical model and three-dimensional (3D) FEM simulation, our model is highly computationally efficient, making it suitable for large-scale real-time simulations of nonlinear problems of shell structures such as instability and snap-through phenomena. Moreover, our framework can easily incorporate complex loading conditions, e.g., boundary nonlinear contact and multi-physics actuation, which play an essential role in the use of engineering applications, such as soft robots and flexible devices. This study demonstrates that the simplicity and effectiveness of the 1D discrete differential geometry-based approach render it a powerful tool for engineers and researchers interested in nonlinear mechanics analysis of axisymmetric shells, with potential applications in various engineering fields.

Monomial-agnostic computation of vanishing ideals

Approximate basis computation of vanishing ideals has recently been studied extensively in computational algebra and data-driven applications such as machine learning. However, symbolic computation and the dependency on term order remain essential gaps between the two fields. In this study, we present the first monomial-agnostic basis computation, which works fully numerically with proper normalization and without term order. This is realized by gradient normalization, a newly proposed data-dependent normalization that normalizes a polynomial with the magnitude of gradients at given points. Its data-dependent nature brings various advantages: i) efficient resolution of the spurious vanishing problem, the scale-variance issue of approximately vanishing polynomials, without accessing coefficients of terms, ii) scaling-consistent basis computation, ensuring that input scaling does not lead to an essential change in the output, and iii) robustness against input perturbations, where the upper bound of error is determined only by the magnitude of the perturbations. Existing studies did not achieve any of these. As further applications of gradient information, we propose a monomial-agnostic basis reduction method and a regularization method to manage positive-dimensional ideals.

Dynamics of quasi-periodic, bifurcation, sensitivity and three-wave solutions for (n + 1)-dimensional generalized Kadomtsev-Petviashvili equation.

This study endeavors to examine the dynamics of the generalized Kadomtsev-Petviashvili (gKP) equation in (n + 1) dimensions. Based on the comprehensive three-wave methodology and the Hirota's bilinear technique, the gKP equation is meticulously examined. By means of symbolic computation, a number of three-wave solutions are derived. Applying the Lie symmetry approach to the governing equation enables the determination of symmetry reduction, which aids in the reduction of the dimensionality of the said equation. Using symmetry reduction, we obtain the second order differential equation. By means of applying symmetry reduction, the second order differential equation is derived. The second order differential equation undergoes Galilean transformation to obtain a system of first order differential equations. The present study presents an analysis of bifurcation and sensitivity for a given dynamical system. Additionally, when an external force impacts the underlying dynamic system, its behavior resembles quasi-periodic phenomena. The presence of quasi-periodic patterns are identified using chaos detecting tools. These findings represent a novel contribution to the studied equation and significantly advance our understanding of dynamics in nonlinear wave models.

Algorithms for representations of quiver Yangian algebras

In this note, we aim to review algorithms for constructing crystal representations of quiver Yangians in detail. Quiver Yangians are believed to describe an action of the BPS algebra on BPS states in systems of D-branes wrapping toric Calabi-Yau three-folds. Crystal modules of these algebras originate from molten crystal models for Donaldson-Thomas invariants of respective three-folds. Despite the fact that this subject was originally at the crossroads of algebraic geometry with effective supersymmetric field theories, equivariant toric action simplifies applied calculations drastically. So the sole pre-requisite for this algorithm’s implementation is linear algebra. It can be easily taught to a machine with the help of any symbolic calculation system. Moreover, these algorithms may be generalized to toroidal and elliptic algebras and exploited in various numerical experiments with those algebras. We illustrate the work of the algorithms in applications to simple cases of Ysl2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \ extrm{Y}\\left({\\mathfrak{sl}}_2\\right) $$\\end{document}, Ygl̂1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \ extrm{Y}\\left({\\hat{\\mathfrak{gl}}}_1\\right) $$\\end{document} and Ygl̂11\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \ extrm{Y}\\left({\\hat{\\mathfrak{gl}}}_{\\left.1\\right|1}\\right) $$\\end{document}.

Multiple rogue wave, double-periodic soliton and breather wave solutions for a generalized breaking soliton system in (3 + 1)-dimensions

We focused on solitonic phenomena in wave propagation which was extracted from a generalized breaking soliton system in (3 + 1)-dimensions. The model describes the interaction phenomena between Riemann wave and long wave via two space variable in nonlinear media. Abundant double-periodic soliton, breather wave and the multiple rogue wave solutions to a generalized breaking soliton system by the Hirota bilinear form and a mixture of exponentials and trigonometric functions are presented. Periodic-soliton, breather wave and periodic are studied with the usage of symbolic computation. In addition, the symbolic computation and the applied methods for governing model are investigated. Through three-dimensional graph, density graph, and two-dimensional design using Maple, the physical features of double-periodic soliton and breather wave solutions are explained all right. The findings demonstrate the investigated model’s broad variety of explicit solutions. All outcomes in this work are necessary to understand the physical meaning and behavior of the explored results and shed light on the significance of the investigation of several nonlinear wave phenomena in sciences and engineering.

Interactions among lump optical solitons for coupled nonlinear Schrödinger equation with variable coefficient via bilinear method

In this paper, a non-autonomous (3+1) dimensional coupled nonlinear Schrödinger equation (NLSE) with variable coefficients in optical fiber communication is analyzed. By means of bilinear technique and symbolic computations, new multi-soliton solutions of the coupled model in different trigonometric and lump functions are given. Then, in terms of perturbed waves, considering the steady state solution and the small perturbation on the three directions x, y, z and the time t, the soliton transmission are also considered. The behaviour of interaction among lump periodic soliton is studied and optical soliton solutions are reached. This study has certain significance for the analysis of other nonlinear dispersion systems and the application of optical physics. The results are presented through graphs generated by using Maple. The important feature of the proposed study is to show different behaviour of the soliton at each component. The behaviour of solitons, their interactions, and their transformations are all governed by the fundamental concept of energy conservation in all three examples. We demonstrate the efficiency of our suggested methodology for analyzing the NLSE equations using the numerical simulations and analytical tools, yielding fresh insights into their behaviour and solutions. Our findings help to develop mathematical tools for investigating nonlinear partial differential equation (NLPDEs) and provide new insights on the dynamics of NLSE equations, which have implications for many domains of physics and applied mathematics

New insight to the shallow-water studies on a (2+1)-dimensional generalized Broer-Kaup system

Of current interest, there have recently appeared a series of the shallow-water papers in Mod. Phys. Lett. B. Illuminated by those papers and to make the story more complete, for the nonlinear long waves in the shallow water, this paper is scheduled to investigate a (2+1)-dimensional generalized Broer-Kaup system. As for the wave height and wave horizontal velocity, we employ symbolic computation, with the view of constructing out two groups of the similarity reductions.