Contour Integral Research Articles

Emergent potentials and non-perturbative open topological strings

We show that integrating out M2 branes ending on M5 branes inside Calabi-Yau manifolds captures non-perturbative open topological string physics. The integrating out is performed using a contour integral in complexified Schwinger proper time. For the resolved conifold, this contour can be extended to include the zero pole, which we argue captures the ultraviolet completion of the integrating out and yields the tree-level polynomial terms in the free energy. This is a manifestation of the Emergence Proposal, and provides further evidence for it. Unlike the case of closed strings, where the emergent terms are kinetic terms in the action, for these open strings it is tree-level potential terms which are emergent. This provides a first quantitative example of the proposal that classical tree-level potentials in string theory emerge from integrating out co-dimension one states.

Analytical Solution for the Problem of Point Location in Arbitrary Planar Domains

This paper presents a general analytical solution for the problem of locating points in planar regions with an arbitrary geometry at the boundary. The proposed methodology overcomes the traditional solutions used for polygonal regions. The method originated from the explicit evaluation of the contour integral using the Residue and Cauchy theorems, which then evolved toward a technique very similar to the winding number and, finally, simplified into a variant of ray-crossing approach slightly more informed and more universal than the classic approach, which had been used for decades. The very close relation of both techniques also emerges during the derivation of the solution. The resulting algorithm becomes simpler and potentially faster than the current state of the art for point locations in arbitrary polygons because it uses fewer operations. For polygonal regions, it is also applicable without further processing for special cases of degeneracy, and it is possible to use in fully integer arithmetic; it can also be vectorized for parallel computation. The major novelty, however, is the extension of the technique to virtually any shape or segment delimiting a planar domain, be it linear, a circular arc, or a higher order curve.

Application of J-integral to adhesive contact under general plane loading for rolling resistance

In the present work, a mechanical model for two-dimensional non-slipping adhesive contact between dissimilar elastic solids under general loading, namely, normal forces, tangential forces and moments is proposed. The general solutions are obtained analytically with the stresses at the contact edges exhibiting oscillatory singularity, similar to those at a bimaterial interface crack. The well-known J-integral under the current context is analyzed. Its application under the selected integration contour readily gives the relationship between the stress intensity factors and energy release rates at the contact edges. With the results rolling adhesion between two solids with parabolic profiles is considered further. The applied moment can be directly determined by the difference in energy release rates at the trailing and leading edges and hence the rolling resistance even for adhesive contact with cohesive zones. These results provide the foundation for understanding some tribological phenomena associated with adhesion.

A fourth order Runge-Kutta type of exponential time differencing and triangular spectral element method for two dimensional nonlinear Maxwell's equations

In this paper, we study a numerical scheme to solve the nonlinear Maxwell's equations. The discrete scheme is based on the triangular spectral element method (TSEM) in space and the exponential time differencing fourth-order Runge-Kutta (ETDRK4) method in time. TSEM has the advantages of spectral accuracy and geometric flexibility. The ETD method involves exact integration of the linear part of the governing equation followed by an approximation of an integral involving the nonlinear terms. The RK4 scheme is introduced for the time integration of the nonlinear terms. The stability region of the ETDRK4 method is depicted. Moreover, the contour integral in the complex plan is utilized and improved to compute the matrix function required by the implementation of ETDRK4. The numerical results demonstrate that our proposed method is of exponential convergence with the order of basis function in space and fourth order accuracy in time.

The semi-classical saddles in three-dimensional gravity via holography and mini-superspace approach

We determine the complex geometries dual to the semi-classical saddles in three-dimensional gravity with positive or negative cosmological constant. We examine the semi-classical saddles in Liouville field theory and interpret them in terms of gravity theory. For this, we describe the gravity theory by Chern-Simons theory and classify the possible saddles based on the homotopy group argument. We further realize the semi-classical saddles using the mini-superspace model of quantum gravity and explicitly determine the integral contour. In the case of positive cosmological constant, we recovered the geometry used for no-boundary proposal of Hartle and Hawking. In the case of negative cosmological constant, the geometry can be identified with Euclidean anti-de Sitter space attached with imaginary radius spheres. The geometry should be unphysical and several arguments on this issue are provided. Partial results were already presented in our earlier letter, and more detailed derivations and explanations on the results are given along with additional results. In particular, we reproduce the classical Liouville action from the Chern-Simons formulation of dual gravity theory.

Preconditioned flow as a solution to the hierarchical growth problem in the generalized Lefschetz thimble method

The generalized Lefschetz thimble method is a promising approach that attempts to solve the sign problem in Monte Carlo methods by deforming the integration contour using the flow equation. Here we point out a general problem that occurs due to the property of the flow equation, which extends a region on the original contour exponentially to a region on the deformed contour. Since the growth rate for each eigenmode is governed by the singular values of the Hessian of the action, a huge hierarchy in the singular value spectrum, which typically appears for large systems, leads to various technical problems in numerical simulations. We solve this hierarchical growth problem by preconditioning the flow so that the growth rate becomes identical for every eigenmode. As an example, we show that the preconditioned flow enables us to investigate the real-time quantum evolution of an anharmonic oscillator with the system size that can hardly be achieved by using the original flow.

The error bounds of Gaussian quadratures for one rational modification of Chebyshev measures

Abstract For an analytic integrand, the error term in the Gaussian quadrature can be represented as a contour integral, where the contour is commonly taken to be an ellipse. Thus, finding its upper bound can be reduced to finding the maximum of the modulus of the kernel on the ellipse. The location of this maximum was investigated in many special cases, particularly, for the Gaussian quadrature with respect to the Chebyshev measures modified by a quadratic divisor (known as the Bernstein–Szeg̋ measures). Here, for the Gaussian quadratures with respect to the Chebyshev measures modified by a linear over linear rational factor, we examine the kernel and describe sufficient conditions for the maximum to occur on the real axis. Furthermore, an assessment of the kernel is made in each case, since in some cases the true maximum is hard to reach. Hence, we derive the error bounds for these quadrature formulas. The results are illustrated by the numerical examples. An alternative approach for estimating the error of the Gaussian quadrature with respect to the same measure can be found in [Djukić, D. L., Djukić, R. M. M., Reichel, L. & Spalević, M. M. (2023, Weighted averaged Gaussian quadrature rules for modified Chebyshev measure. Appl. Numer. Math., ISSN 0168-9274)].

2D numerical modeling of crack propagation using SFEMD method of a LEHI material

Numerical methods today play a useful and important role in solving various problems related to fracture mechanics including modeling crack propagation, fretting fatigue and cohesion... Furthermore, these methods have been widely used to solve problems in linear and nonlinear fracture mechanics in cases of elastic, plastic fracture problems. The evaluation of stress intensity factor in 2D and 3D geometries, thus these techniques widely used for non-standard crack configurations. The objective of this work is study the effects of the crack length and the number of structural elements on the two crack parameters such as the stress intensity factor KI and KII and the contour integral (J). As well as presenting a numerical modeling of crack propagation, for a LEHIM (Linear Elastic Homogeneous Isotopic Material) with mechanical characteristics by the method SFEMD (Stretching Finite Element Method Developed), the model chosen with quadratic elements with 4 nodes (CPE4). However, this method is based on the effect of the number of contours and the number of elements around the crack tip on the variation of the stress intensity factors. In addition, the crack propagation criterion (MCSC) was used, a computer program was created in FORTRAN language to develop and evaluate the stress intensity factors and the contour integral (J). Several examples of crack lengths a = 0.7, 1.4, 2.1, 2.8 and 3.5 mm were used. Additionally, the number of items has been changed several times. The stress intensity factors of modes I and II and direction angles (α) are calculated to solve the problem using the ABAQUS finite element code. The results obtained by the SFEMD method and the results of the analytical method are very close.

Simplex slicing: An asymptotically-sharp lower bound

We show that for the regular n-simplex, the 1-codimensional central slice that's parallel to a facet will achieve the minimum area (up to a 1−o(1) factor) among all 1-codimensional central slices, thus improving the previous best known lower bound (Brzezinski 2013) by a factor of 23e≈1.27. In addition to the standard technique of interpreting geometric problems as problems about probability distributions and standard Fourier-analytic techniques, we rely on a new idea, mainly changing the contour of integration of a meromorphic function.

A Modification to the Novel Toy Model of the Riemann Zeta Function Roots Equation

It is true that my previous series of papers about the Riemann Hypothesis has raised lots of discussions or the concerns among the mathematics society. Indeed, some professionals may think that the Riemann Zeta function at s = 1 or ∑_(n=1)^∞▒1/n is actually divergent and may tend to an infinity, how it can be acted as a denominator in a proof or something like 1/∞ ? This author’s idea is there may be no need to evaluate such kind of summation as one may consider the sum just as an (imaginary) constant (but NOT the complex number). Or we may only consider such summation as a kind of improper integrals with a black-boxed answer. In other words, the true final value of the summation ∑_(n=1)^∞▒1/n is NOT an concerning serious issue. The function ∑_(n=1)^∞▒1/n^s is in practice a smooth function and equals to log N or the logarithmic function for s = 1 or have the infinite differentiable properties except n = 0, that is why the commercial software Maple Soft can compute the Zeta function’s associated Taylor series successfully. In fact, from the above series together with the harmonic properties, one may calculate the artifical Zeta function’s root model equation such as the 0.5 +/- β*cot(ln(x))/(αx+1)n (will be discussed in my next paper of the series). Indeed, this is NOT the best model. Actually, the true Riemann Zeta non-trivial Root model equation is: {z = x + yI | z∈0.5±(y_1⩽y⩽y_2 )I} where y1= (±cot(k))/ln(k) , y2 = (±tan(k))/ln(k) with 0<k<2π. Certainly, this author cannot avoid mistakes. Some of the defects in the previous paper such as one in the linguistic proof of the RH will be stated and amended. Lastly, this author wants to remind that the present series of paper is first to establish a model equation for zeta function’s roots and the Plank’s like constant. This author then applies the telescopic and logarithmic methods for a 3 cases proof of the RH. Also, we have a Matlab code for the evaluation of RH complex contour integral. Finally, the RH problem will be transformed into a linguistic one and hence we may have solved the Riemann Hypothesis issue completely through a 4 cases of truth table collaborations with the logical inference [7] & [15] to the causality (structural) organization of the sentences [16] or just named as “Logical & Organized Context”.

KPZ fluctuations in finite volume

These lecture notes, adapted from the habilitation thesis of the author, survey in a first part various exact results obtained in the past few decades about KPZ fluctuations in one dimension, with a special focus on finite volume effects describing the relaxation to its stationary state of a finite system starting from a given initial condition. The second part is more specifically devoted to an approach allowing to express in a simple way the statistics of the current in the totally asymmetric simple exclusion process in terms of a contour integral on a compact Riemann surface, whose infinite genus limit leads to KPZ fluctuations in finite volume.

Fermionic sign problem minimization by constant path integral contour shifts

Fermionic sign problem minimization by constant path integral contour shifts

A Contour Integral and Complex Power Series Proof of Stirling Formula

A Contour Integral and Complex Power Series Proof of Stirling Formula

On the integral solution of hyperbolic Kepler’s equation

In a recent paper of Philcox, Goodman and Slepian, the solution of the elliptic Kepler’s equation is given as a quotient of two contour integrals along a Jordan curve that contains in its interior the unique real solution of the elliptic Kepler’s equation and does not include other complex zeroes. In this paper, we show that a similar explicit integral solution can be given for the hyperbolic Kepler’s equation. With this purpose, we carry out a study of the complex zeros of the hyperbolic Kepler’s equation in order to define suitable Jordan contours in the integrals. Even more, we show that appropriate elliptic Jordan contours can be defined for such integrals, which reduces the computing time. Moreover, using the ideas behind the fast Fourier transform (FFT) algorithm, these integrals can be approximated by the composite trapezoidal rule which gives an algorithm with spectral convergence as a function of the number of nodes. The results of some numerical experiments are presented to show that this implementation is a reliable and very accurate algorithm for solving the hyperbolic Kepler’s equation.

Version 2 — RPExpand: Software for Riesz projection expansion of resonance phenomena

The software RPExpand, implemented in MATLAB®, provides modal expansions of physical observables in resonant systems, based on contour integrals. Version 2.0 improves the workflow and adds new features. In addition to the poles corresponding to eigenmodes of the systems, zeros of various quantities can be computed. Their derivatives are evaluated if the required information is supplied at the integration points. Furthermore, expansions based on quasinormal modes are supported.

Non-polynomial q-Askey Scheme: Integral Representations, Eigenfunction Properties, and Polynomial Limits

AbstractWe construct a non-polynomial generalization of the q-Askey scheme. Whereas the elements of the q-Askey scheme are given by q-hypergeometric series, the elements of the non-polynomial scheme are given by contour integrals, whose integrands are built from Ruijsenaars’ hyperbolic gamma function. Alternatively, the integrands can be expressed in terms of Faddeev’s quantum dilogarithm, Woronowicz’s quantum exponential, or Kurokawa’s double sine function. We present the basic properties of all the elements of the scheme, including their integral representations, joint eigenfunction properties, and polynomial limits.

On the Sums over Inverse Powers of Zeros of the Hurwitz Zeta Function and Some Related Properties of These Zeros

Recently, we have applied the generalized Littlewood theorem concerning contour integrals of the logarithm of the analytical function to find the sums over inverse powers of zeros for the incomplete gamma and Riemann zeta functions, polygamma functions, and elliptical functions. Here, the same theorem is applied to study such sums for the zeros of the Hurwitz zeta function ζ(s,z), including the sum over the inverse first power of its appropriately defined non-trivial zeros. We also study some related properties of the Hurwitz zeta function zeros. In particular, we show that, for any natural N and small real ε, when z tends to n = 0, −1, −2… we can find at least N zeros of ζ(s,z) in the ε neighborhood of 0 for sufficiently small |z+n|, as well as one simple zero tending to 1, etc.

Numerical Analysis of Residual Stresses Effect in Multi-materials

Background: In this study, the effect of mechanical and physical properties of the metal and temperature is highlighted. The purpose of this study is to analyze the effect of these stresses on multi-materials. This work aimed to conduct a three-dimensional numerical study by the finite element method on the levels and distributions of the stresses in the Al2O3/NI/HAYNES multi-materials. These stresses of thermal and mechanical origin are generally detrimental to the service life of multi-material. Methods: The use of numerical resolution by finite element method is the most suitable for complex mechanical problems. It allows a more in-depth analysis of all points of the structure. In this study, a fundamental tool was constructed to resolve the mechanical behavior of materials subjected to complex solicitations. Therefore, the ABAQUS calculation code version 6.14 was used to analyze the residual stresses. Results: By interpreting the results in terms of stress variation, we identified the areas at risk; in particular, the nature of a joint effect plays a decisive role in the assembly of the right material in its mechanical resistance. This nature is defined in terms of stiffness (Young's modulus) and differential expansion (coefficient of thermal expansion). Conclusion: The presence of strong residual stresses can constitute a risk of damage to several materials. The difference between the coefficients of thermal expansion of the two materials (metal and ceramic) linked together induces, in these two constituents, normal internal stresses. This difference determines the level and distribution of these constraints. Moreover, the sign of this difference determines the state of the normal stresses.

An isoparametric inclusion model for determining the thermo-elastic fields produced by varying Eigen-temperature gradients

The existence of inclusions or inhomogeneities may disrupt heat flow, resulting in increased stresses and temperature fluctuations at the material interface. Analyzing the thermo-elastic fields surrounding such microstructures is crucial for unveiling the intricate failure mechanisms within materials. Based on the method of Green's function and contour integral, this work presents a complete set of explicit analytical solutions for thermal and elastic fields of an arbitrary polygonal inclusion subjected to linearly varying eigen-temperature gradients or thermal eigenstrains. In contrast to the previous studies on potentials showing analogy with anti-plane elasticity, a mathematical analogy of Green's functions between steady-state heat conduction and thermal inclusion for plane elasticity is accordingly established. By employing the linear interpolation of the nodal eigen-temperature gradients or eigenstrains, the work further extends the traditional isoparametric element idea to the thermos-elastic inclusions, and proposes a novel isoparametric inclusion model to solve arbitrarily shaped inclusions with any distributed eigen-temperature gradients. Unlike the finite element method, which requires meshing for both the matrix and inclusion, the present model exhibits efficiency, flexibility, and versatility with the mesh discretization only conducted inside the inclusion domain. The inclusion based isoparametric method may have potential applications in mechanical engineering and material science for analyzing the thermo-structural behavior of various microstructures.

Lifting of two-mode states in the D1-D5 CFT

We consider D1-D5-P states in the untwisted sector of the D1-D5 orbifold CFT where one copy of the seed CFT has been excited by a pair of oscillators, each being either bosonic or fermionic. While such states are BPS at the orbifold point, they will in general ‘lift’ as the theory is deformed towards general values of the couplings. We compute the expectation value of this lift at second order in the deformation parameter for the above mentioned states. We write this lift in terms of a fixed number of nested contour integrals on a given integrand; this integrand depends on the mode numbers of the oscillators in the state. We evaluate these integrals to obtain the explicit value of the lift for various subfamilies of states. At large mode numbers one observes a smooth increase of the lift with the dimension of the state h; this increase appears to follow a∼h\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \ extrm{a}\\sim \\sqrt{h} $$\\end{document} behavior similar to that found analytically in earlier computations for other classes of states.