Methods For Differential Equations Research Articles

Fourier calculus from intersection theory

Building on recent advances in studying the cohomological properties of Feynman integrals, we apply intersection theory to the computation of Fourier integrals. We discuss applications pertinent to gravitational bremsstrahlung and deep inelastic scattering in the saturation regime. After identifying the bases of master integrals, the latter are evaluated by means of the differential equation method. Finally, new results with exact dependence on the spacetime dimension D are presented. Published by the American Physical Society 2024

Galerkin-finite difference method for fractional parabolic partial differential equations

The fractional form of the classical diffusion equation embodies the super-diffusive and sub-diffusive characteristics of any flow, depending on the fractional order. This study aims to approximate the solution of parabolic partial differential equations of fractional order in time and space. For this, firstly, we briefly discussed on some existing methods to solve Partial Differential Equations (PDEs) and fractional Differential Equations (DEs), then introduce a combined technique such as the Galerkin weighted residual method for the space fractional term with modified Bernoulli polynomials as basis functions, and the finite difference approximation for the time fractional term, respectively. The mathematical formulation of the proposed method is explained elaborately. Then we describe the order of convergence for the time fractional term only, as the convergence of the Galerkin method is obvious. We impose this technique on the fractional Black-Scholes model subsequently. Finally, we experimented with our proposed technique on some numerical problems. All the results are depicted in both tabular and 3D visualizations as well. We compare our results with the available methods in the literature, and our accuracy is considerable. To summarize: •The paper introduces an approach that integrates the Galerkin weighted residual method with modified Bernoulli polynomials to handle space fractional terms, alongside employing a finite difference approximation for time fractional terms.•The convergence analysis is focused.•The technique is implemented on the fractional Black-Scholes model and other numerical problems, with outcomes depicted through tables and 3D visualizations.

Continuous-discrete extended Kalman filtering based on the neural ordinary differential equations method

In bearings-only tracking (BOT), the accuracy of state estimation significantly diminishes when confronted with an unknown target motion model. Therefore, this study proposed a novel continuous-discrete filtering algorithm, the neural ordinary differential equations-based continuous-discrete extended Kalman filtering (NODE-CD-EKF). The approach utilizes the neural ordinary differential equation (NODE) to build the target's system model, which is then used in the time update to precisely predict the target state at the next time point within the continuous-discrete Kalman framework. The NODE concisely and accurately describes the motion model in the form of ordinary differential equations in a data-driven way. This approach addresses the limitations of traditional Kalman filters, which suffer from degraded performance due to a lack of a precise mathematical model of target states. Our method is evaluated and compared to existing multiple-model and other data-driven algorithms in simulations of BOT and another scenario using the Accumulated Root Mean Square Error (ARMSE) of position. The ARMSE of our method is only 16.30 % of the multiple-model methods at least. Overall, this approach outperforms other multiple-model and data-driven algorithms and possesses the potential for solving other nonlinear tasks in real-world applications.

Finite Elements with Switch Detection for direct optimal control of nonsmooth systems

This paper introduces Finite Elements with Switch Detection (FESD), a numerical discretization method for nonsmooth differential equations. We consider the Filippov convexification of these systems and a transformation into dynamic complementarity systems introduced by Stewart (Numer Math 58(1):299–328, 1990). FESD is based on solving nonlinear complementarity problems and can automatically detect nonsmooth events in time. If standard time-stepping Runge–Kutta (RK) methods are naively applied to a nonsmooth ODE, the accuracy is at best of order one. In FESD, we let the integrator step size be a degree of freedom. Additional complementarity conditions, which we call cross complementarities, enable exact switch detection, hence FESD can recover the high order accuracy that the RK methods enjoy for smooth ODE. Additional conditions called step equilibration allow the step size to change only when switches occur and thus avoid spurious degrees of freedom. Convergence results for the FESD method are derived, local uniqueness of the solution and convergence of numerical sensitivities are proven. The efficacy of FESD is demonstrated in several simulation and optimal control examples. In an optimal control problem benchmark with FESD, we achieve up to five orders of magnitude more accurate solutions than a standard time-stepping approach for the same computational time.

Analysis of Natural Daftardar–Jafari Method for Fractional Delay Differential Equations

Analysis of Natural Daftardar–Jafari Method for Fractional Delay Differential Equations

On minimum vertex bisection of random d-regular graphs

Minimum vertex bisection is a graph partitioning problem in which the aim is to find a partition of the vertices into two equal parts that minimizes the number of vertices in one partition set that have a neighbor in the other set. In this work we are interested in providing asymptotically almost surely upper bounds on the minimum vertex bisection of random d-regular graphs, for constant values of d. Our approach is based on analyzing a greedy algorithm by using the differential equation method. In this way, we obtain the first known non-trivial upper bounds for the vertex bisection number in random regular graphs. The numerical approximations of these theoretical bounds are compared with the emprical ones, and with the lower bounds from Kolesnik and Wormald (2014) [30].

An adaptive certified space-time reduced basis method for nonsmooth parabolic partial differential equations

In this paper, a nonsmooth semilinear parabolic partial differential equation (PDE) is considered. For a reduced basis (RB) approach, a space-time formulation is used to develop a certified a-posteriori error estimator. This error estimator is adopted to the presence of the discrete empirical interpolation method (DEIM) as approximation technique for the nonsmoothness. The separability of the estimated error into an RB and a DEIM part then guides the development of an adaptive RB-DEIM algorithm, combining both offline phases into one. Numerical experiments show the capabilities of this novel approach in comparison with classical RB and RB-DEIM approaches.

Near mass-shell double boxes

Two-loop multi-leg form factors in off-shell kinematics require knowledge of planar and nonplanar double box Feynman diagrams with massless internal propagators. These are complicated functions of Mandelstam variables and external particle virtualities. The latter serve as regulators of infrared divergences, thus making these observables finite in four space-time dimensions. In this paper, we use the method of canonical differential equations for the calculation of (non)planar double box integrals in the near mass-shell kinematical regime, i.e., where virtualities of external particles are much smaller than the Mandelstam variables involved. We deduce a basis of master integrals with uniform transcendental weight based on the analysis of leading singularities employing the Baikov representation as well as an array of complementary techniques. We dub the former asymptotically canonical since it is valid in the near mass-shell limit of interest. We iteratively solve resulting differential equations up to weight four in terms of multiple polylogarithms.

Higher-loop integrated negative geometries in ABJM

In the three-dimensional \\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathcal{N}$$\\end{document} = 6 Chern-Simons matter (ABJM) theory, the integrand for the logarithm of the scattering amplitude admits a decomposition in terms of negative geometries, which implies that all the infrared divergences concentrate in the last loop integration. We compute the infrared-finite functions that arise from performing a three-loop integration over the four-loop integrand for the logarithm of the four-point amplitude, for which we use the method of differential equations. Our results provide a direct computation of the four-loop cusp anomalous dimension of the theory, in agreement with the current all-loop integrability-based proposal. We find an apparent simplicity in the leading singularities of the integrated results, provided one works in the frame in which the unintegrated loop variable goes to infinity. Finally, our results suggest an alternating sign pattern for the integrated negative geometries in the Euclidean region.

A Collocation Method for Nonlinear Stochastic Differential Equations Driven by Fractional Brownian Motion and its Application to Mathematical Finance

A Collocation Method for Nonlinear Stochastic Differential Equations Driven by Fractional Brownian Motion and its Application to Mathematical Finance

Alleviated Shifted Gegenbauer Spectral Method for Ordinary and Fractional Differential Equations

This article introduces two numerical methods to address boundary value problems associated with secondorder and fractional differential equations. These methods employ two parameters related to shifted Gegenbauer polynomials as their basis functions. The process involves establishing a differentiation operational matrix for the shifted Gegenbauer polynomials. Subsequently, the initial/boundary value problems for ordinary and fractional differential equations are transformed into a system of equations through the Galerkin, collocation, and tau methods. The convergence analysis is ensured by leveraging theorems pertaining to the shifted Gegenbauer polynomials. To validate the accuracy of the approach, numerous numerical examples are presented.

Two-loop non-planar four-point topology with massive internal loop

We study a set of two-loop non-planar master integrals needed for the NNLO QCD corrections to diphoton and dijet production at hadron colliders. The top-sector topology contains an internal massive fermion loop and is known to contain elliptic curves. Leveraging the method of differential equations, we provide a comprehensive discussion for deriving an ϵ-factorized differential equation related to the most intricate sector within the Feynman integral family. Despite the dependence on multiple scales and the presence of two elliptic sectors, we demonstrate how to leverage the properties of their maximal cuts and the factorization of the Picard-Fuchs operator to deal with the complexity of the analytic computation. In particular, we construct a transformation matrix that brings the differential equations into a format enabling the convenient expression of analytic results in terms of Chen’s iterated integrals.

On the convergence of the Galerkin method for random fractional differential equations

On the convergence of the Galerkin method for random fractional differential equations

Tamed Stochastic Runge-Kutta-Chebyshev Methods for Stochastic Differential Equations with Non-Globally Lipschitz Coefficients

Tamed Stochastic Runge-Kutta-Chebyshev Methods for Stochastic Differential Equations with Non-Globally Lipschitz Coefficients

Convergence order of one point large deviations rate functions for backward Euler method of stochastic delay differential equations with small noise

The present work focuses on the convergence order of one point large deviations rate functions for the backward Euler method of stochastic delay differential equations (SDDEs) with small noise. The drift and diffusion coefficients of SDDEs are allowed to grow super-linearly with respect to both the state variables and the delay variables. It is shown that the backward Euler method satisfies the one point large deviations principle (LDP) with good rate function. Further, the local uniform convergence orders of the one point large deviations rate function are derived by means of the equivalent characterizations for the original and numerical rate functions. These theoretical results are finally supported by a series of numerical experiments.

Rotationally symmetric transverse magnetic vector wave propagation for nonlinear optics

In this paper the theory and simulation results are presented for 3D cylindrical rotationally symmetric spatial soliton propagation in a nonlinear medium using a modified finite-difference time-domain general vector auxiliary differential equation method for transverse magnetic polarization. The theory of 3D rotationally symmetric spatial solitons is discussed, and compared with two (1 + 1)D, termed “2D” for this paper, hyperbolic secant spatial solitons, with a phase difference of π (antiphase). The simulated behavior of the 3D rotationally symmetric soliton was compared with the interaction of the two antiphase 2D solitons for different source hyperbolic secant separation distances. Lastly, we offer some possible explanations for the simulated soliton behavior.

A Streamline Upwind Petrov-Galerkin Reduced Order Method for Advection-Dominated Partial Differential Equations Under Optimal Control

Abstract In this paper we will consider distributed Linear-Quadratic Optimal Control Problems dealing with Advection-Diffusion PDEs for high values of the Péclet number. In this situation, computational instabilities occur, both for steady and unsteady cases. A Streamline Upwind Petrov–Galerkin technique is used in the optimality system to overcome these unpleasant effects. We will apply a finite element method discretization in an optimize-then-discretize approach. Concerning the parabolic case, a stabilized space-time framework will be considered and stabilization will also occur in both bilinear forms involving time derivatives. Then we will build Reduced Order Models on this discretization procedure and two possible settings can be analyzed: whether or not stabilization is needed in the online phase, too. In order to build the reduced bases for state, control, and adjoint variables we will consider a Proper Orthogonal Decomposition algorithm in a partitioned approach. It is the first time that Reduced Order Models are applied to stabilized parabolic problems in this setting. The discussion is supported by computational experiments, where relative errors between the FEM and ROM solutions are studied together with the respective computational times.

Runge-Kutta Method Based on the Linear Combination of Arithmetic Mean, Geometric Mean, and Centroidal Mean: A Numerical Solution for the Hybrid Fuzzy Differential Equation

Objectives: To find the absolute error of the Hybrid fuzzy differential equation with triangular fuzzy number since its membership function has a triangular form. Methods: The third order runge-kuttta method based on the linear combination of Arithmetic mean, Geometric mean and Centroidal mean is introduced to solve Hybrid fuzzy differential equation. Seikkala’s derivatives are considered, and Numerical example is given to demonstrate the effectiveness of the proposed method. Findings: The comparative study have been done between the proposed method and the existing third order runge-kutta method based on Arithmetic mean. The proposed method gives better error tolerance than the third order runge-kutta method based on Arithmetic mean. Novelty: In this study a new formula has been developed by combining three means Arithmetic Mean, Geometric Mean and Centroidal Mean using Khattri’s formula and it is solved by the third order runge-kutta method for the first order hybrid fuzzy differential equation. The error tolerance has been calculated for the proposed method and the Arithmetic mean and it is seen that the error tolerance is better for the proposed method. Keywords: Hybrid fuzzy differential equations, Triangular fuzzy number, Seikkala’s derivative, Third order runge-kutta method, Initial value problem

The improvement of the truncated Euler-Maruyama method for non-Lipschitz stochastic differential equations

The improvement of the truncated Euler-Maruyama method for non-Lipschitz stochastic differential equations

A new method for calculating the soft anomalous dimension matrix for massive particle scattering

The general structure of infrared divergences in the scattering of massive particles is captured by the soft anomalous dimension matrix. The latter can be computed from a correlation function of multiple Wilson lines. The state-of-the-art two-loop result has a tantalizingly simple structure that is not manifest in the calculations. We argue that the complexity in intermediate steps of the known calculations comes from spurious, regulator-dependent terms. Based on this insight we propose a different infrared regulator that is associated to only one of the Wilson lines. We demonstrate that this streamlines obtaining the two-loop result: computing the required Feynman integrals via the differential equations method, only multiple polylogarithmic functions appear (to all orders in the dimensional regulator), as opposed to elliptic polylogarithms. We show that the new method is promising for higher-loop applications by computing a three-loop diagram of genuine complexity, and provide the answer in terms of multiple polylogarithms. The relatively simple symbol alphabet we obtain may be of interest for bootstrap approaches.