The two loop crossed ladder vertex diagram with two massive exchanges

A numerical routine for the crossed vertex diagram with a massive-particle loop

We present analytical formulas for the calculation of the two-center two-electron integrals in the basis of Slater geminals and products of Slater orbitals. Our derivation starts with establishing a inhomogeneous fourth-order ordinary differential equation that is obeyed by the master integral, the simplest integral with inverse powers of all interparticle distances. To solve this equation it was necessary to introduce a new family of special functions which are defined through their series expansions around regular singular points of the differential equation. To increase the power of the interparticle distances under the sign of the integral we developed a family of open-ended recursion relations. A handful of special cases of the integrals is also analysed with some remarks on simplifications that occur. Additionally, we present some numerical examples of the master integral that validate the usefulness and correctness of the key equations derived in this paper. In particular, we compare our results with the calculations based on the series expansion of the exp(-\gamma r12) term in the master integral.

In this work, we propose a method for the derivation of prolate spheroidal wave functions (PSWFs) and Slepian functions on continuous and disjoint intervals on the real number line. The proposed method uses Fourier series to obtain a closed-form approximation for Slepian functions on the real line. With this closed-form expression, Slepian functions can be evaluated at arbitrary points in the region of interest with high accuracy. The conventional method uses properties of the Slepian concentration problem to evaluate PSWFs on finite number of points in an interval. The conventional method is computationally expensive and does not allow for easy storage. By approximating an interval containing regions of interest as periodic, we express the Slepian concentration problem as a finite dimensional problem using the Fourier series domain. Solutions to the Slepian concentration problem in this form are Fourier series coefficients corresponding to the Slepian functions. Reconstruction in Fourier series basis, scaling and subsequent truncation provides the closed-form expression for the Slepian problem. Upon comparison with PSWFs obtained by the conventional method, we find negligible difference.

In this paper, a new numerical method, the half boundary method (HBM) is studied. HBM reduced the order of the basic equation by substitution, and then calculated the relationship between the adjacent two nodes on the geometric model by integral solution, and recurred to the relationship between any node on the model and the boundary point. HBM is a small-storage, high-precision, fast and efficient numerical calculation method, with few unknowns and non-iterative innovation nodes, which can meet the requirements of accurate acquisition of real-time temperature field in the case of less computational resources. HBM calculation schemes of one and two dimensions in Cartesian coordinate system are derived in detail. The result of HBM calculation is compared with the result of FVM calculation in the corresponding example. It is found that when the number of nodes is small, the calculation accuracy of HBM is much higher than that of FVM. When the number of nodes is large, both HBM and FVM have high accuracy, but HBM calculation time is shorter. HBM calculation results is provided by MATLAB programs.

A bstract In this paper, we present analytical results for the two-loop QCD corrections to the production of two partons or a photon and a parton in hadronic collisions, mediated by loops of massive quarks. These amplitudes involve Feynman integrals defined on an elliptic curve. We compute them by generalizing our recent results for the production of two photons to include additional crossings of the corresponding master integrals, which we compute in terms of the same basis of independent iterated integrals. We discuss the analytical properties of the amplitudes, highlighting the cancellations of a large number of elliptic differential forms in their finite remainders. Finally, we elaborate on a strategy for their numerical evaluation based on generalized series expansions at singular points of the physical amplitude, through the introduction of suitable sets of variables that allow us to resolve all singularities.

We compute the top quark contribution to the two-loop amplitude for on-shell Z boson pair production in gluon fusion, gg → ZZ. Exact dependence on the top quark mass is retained. For each phase space point the integral reduction is performed numerically and the master integrals are evaluated using the auxiliary mass flow method, allowing fast computation of the amplitude with very high precision.

We calculate master integrals for the two-loop QED corrections to e+e− → 2γ in terms of generalized power series with respect to electron mass. The coefficients of this series are expressed via Goncharov’s polylogarithms. Our approach exploits a number of modern multiloop methods: IBP reduction, differential equations for master integrals, Frobenius method, reduction to ϵ-form, and DRA method.

The work addresses studying the possibilities of a numerical-analytical variant of the boundary elements method (BEM) in determining the internal force factors and stresses at singular points when bending thin isotropic plates. The simplest type of singularity has been investigated, a point of application of external concentrated forces and moments. The importance of a given problem is due to the fact that at these points the internal force factors tend towards infinity and it is not possible to determine the size using elementary methods. At the same time, these singular points are the significant stress concentrators (both tangential and normal), which is why calculating the limits to which the internal forces and moments tend is essential to analyze the strength of plate structures. In order to describe an external load, it is proposed to apply the Dirac delta function of two variables. The models of external loads are presented. A given proposal makes it possible to accurately calculate the limits to which the transverse forces, as well as bending and torsional moments, tend at singular points of thin plates. We simulated plate bending using the variational Kantorovich-Vlasov method, which is fully compatible with the models of external load. The internal force factors at the singular points of plates were determined while solving the boundary value problems, formed based on the algorithm of BEM. The MATLAB environment was used for programming and computation. Results of the calculations are characterized by high accuracy and reliability, in particular the errors in determining the deflections of plates at singular points do not exceed 2.0 % and the errors for bending moments are not above 3.0 %. Recommendations have been given to solving different types of boundary problems on bending the plates with singular points based on the proposed approach. It has been established that an accurate model of the external load in the form of concentrated forces and moments fundamentally enables determining the internal forces and moments at the singular points of thin plates applying an algorithm of the variational Kantorovich-Vlasov method. Up to now, there are no data on the importance of internal forces and moments at the singular points of plates. It is also shown that when calculating the internal forces and moments of plates, it is inappropriate to apply a single term from a series of the Kantorovich-Vlasov method; the errors amount to significant magnitudes of the order of 43‒44 %

We present a new method for visualizing implicit real algebraic curves inside a bounding box in the $2$-D or $3$-D ambient space based on numerical continuation and critical point methods. The underlying techniques work also for tracing space curve in higher-dimensional space. Since the topology of a curve near a singular point of it is not numerically stable, we trace only the curve outside neighborhoods of singular points and replace each neighborhood simply by a point, which produces a polygonal approximation that is $\epsilon$-close to the curve. Such an approximation is more stable for defining the numerical connectedness of the complement of the projection of the curve in $\mathbb{R}^2$, which is important for applications such as solving bi-parametric polynomial systems. The algorithm starts by computing three types of key points of the curve, namely the intersection of the curve with small spheres centered at singular points, regular critical points of every connected components of the curve, as well as intersection points of the curve with the given bounding box. It then traces the curve starting with and in the order of the above three types of points. This basic scheme is further enhanced by several optimizations, such as grouping singular points in natural clusters, tracing the curve by a try-and-resume strategy and handling "pseudo singular points". The effectiveness of the algorithm is illustrated by numerous examples. This manuscript extends our preliminary results that appeared in CASC 2018.

Analysis of traffic flow signals plays an important role in traffic prediction and management. As an intrinsic property, the singular point of a traffic flow signal labels a new nonsteady status. Therefore, detecting the singular point is an effective approach to determine the moment of traffic flow prediction. In this paper, an improved wavelet transform is proposed to detect singular points of real-time traffic flow signals. The number of detected singular points is output via the heuristic selection of multiple scales. Then, a weighted similarity measurement of historical traffic flow signals is utilized to predict the next singular point. The position of the next singular point decides the duration of prediction adaptively. The detected and predicted singular points are applied to dynamically update the unequal interval prediction of traffic flow. Furthermore, a Vasicek model is used to predict the traffic flow by minimizing the sum of the relative mean standard error (RMSE) between the traffic flow increment in the predicted interval and the sampled increments of previous intervals. A decomposition method is used to solve the unequal matrix problem. Based on the scenario and traffic flow imported from the real-world map, the simulation results show that the proposed algorithm outperforms existing approaches with high prediction accuracy and much lower computing cost.

It is a very demanding task to design a reliable fingerprint matching approach with high accuracy and speed. An algorithm based on the novel structure combining the singular point with its neighborhood minutiae proposed in this paper can solve this problem efficiently. The structure introduced in this paper has two novel ideas as follows: First, we give an efficient singular points detection method by the inter-relationship between the singular points and the minutiae around them. It can reject the spurious singular points detected by the Poincare index. Second, an improvement of minutiae pairing strategy is introduced, which can sharply decrease the number of candidate minutiae pairs. Experiment results show that these improvements can highly speed up the matching with a preferable accuracy. This algorithm can be used in the one-to-many matching of the on-line fingerprint identification system.KeywordsFingerprintMatchingSingular PointMinutiae Pairing

The Stokes phenomenon in the confluence of the hypergeometric equation using Riccati equation

Optimized computation of the voigt and complex probability functions

We consider a regular singular Sturm‐Liouville operator equation image on the line segment (0,1]. We impose certain boundary conditions such that we obtain a semi‐bounded self‐adjoint operator. It is known (cf. Theorem 1.1 below) that the ζ‐function of this operator equation image has a meromorphic continuation to the whole complex plane with 0 being a regular point. Then, according to [RS] the ζ ‐ regularized determinant of L is defined by equation image In this paper we are going to express this determinant in terms of the solutions of the homogeneous differential equation Ly = 0 generalizing earlier work of S. Levit and U. Smilansky [LS], T. Dreyfus and H. Dym [DD], and D. Burghelea, L. Friedlander and T. Kappeler [BFK1, BFK2). More precisely we prove the formula equation image Here ϕ ψ is a certain fundamental system of solutions for the homogeneous equation Ly = 0, W(ϕ ψ), denotes their Wronski determinant, and v0, v1 are numbers related to the characteristic roots of the regular singular points 0, 1.

In the paper we deal with the Heun functions --- solutions of the Heun\nequation, which is the most general Fuchsian equation of second order with four\nregular singular points. Despite the increasing interest to the equation and\nnumerous applications of the functions in a wide variety of physical problems,\nit is only Maple amidst known software packages which is able to evaluate the\nHeun functions numerically. But the Maple routine is known to be imperfect:\neven at regular points it may return infinities or end up with no result.\nImproving the situation is difficult because the code is not publicly\navailable. The purpose of the work is to suggest and develop alternative\nalgorithms for numerical evaluation of the Heun functions. A procedure based on\npower series expansions and analytic continuation is suggested which allows us\nto avoid numerical integration of the differential equation and to achieve\nreasonable efficiency and accuracy. Results of numerical tests are given.\n