Self-similar factor approximants for evolution equations and boundary-value problems

The method of self-similar factor approximants is shown to be very convenient for solving different evolution equations and boundary-value problems typical of physical applications. The method is general and simple, being a straightforward two-step procedure. First, the solution to an equation is represented as an asymptotic series in powers of a variable. Second, the series are summed by means of the self-similar factor approximants. The obtained expressions provide highly accurate approximate solutions to the considered equations. In some cases, it is even possible to reconstruct exact solutions for the whole region of variables, starting from asymptotic series for small variables. This can become possible even when the solution is a transcendental function. The method is shown to be more simple and accurate than different variants of perturbation theory with respect to small parameters, being applicable even when these parameters are large. The generality and accuracy of the method are illustrated by a number of evolution equations as well as boundary value problems.     

Optimisation of self-similar factor approximants

The problem is analysed for extrapolating power series, derived for an asymptotically small variable, to the region of finite values of this variable. The consideration is based on the self-similar approximation theory. A new method is suggested for defining the odd self-similar factor approximants by employing an optimisation procedure. The method is illustrated by several examples having the mathematical structure typical of the problems in statistical and chemical physics. It is shown that the suggested method provides a good accuracy even when the number of terms in the perturbative power series is small.     

Calculation of critical exponents by self-similar factor approximants

The method of self-similar factor approximants is applied to calculating the critical exponents of the O(N)-symmetric ϕ⁴ theory and of the Ising glass. It is demonstrated that this method, being much simpler than other known techniques of series summation in calculating the critical exponents, at the same time, yields the results that are in very good agreement with those of other rather complicated numerical methods. The principal advantage of the method of self-similar factor approximants is the combination of its extraordinary simplicity and high accuracy.     

A method for summation of perturbation series

A method for summation of perturbation series is developed. It consists of a proper rearrangement of the power series through an appropriate redefinition of the perturbation parameter. Well-known divergent power series appearing in a φ⁴-scalar field theory or in the treatment of the linear confining potential model and the Stark effect in hydrogen are used as illustrative examples. A comparison between present results and accurate, numerical ones shows that our series converge even faster than sequences of Padé or Borel-Padé approximants.     

Solution of singular two-point boundary value problems using differential transformation method

In this Letter, we implemented relatively new, exact series method of solution known as the differential transform method for solving singular two-point boundary value problems. Several illustrative examples are given to demonstrate the effectiveness of the present method.     

Exp-function method for N-soliton solutions of nonlinear evolution equations in mathematical physics

In this Letter, the Exp-function method is generalized to construct N-soliton solutions of a KdV equation with variable coefficients. As a result, 1-soliton, 2-soliton and 3-soliton solutions are obtained, from which the uniform formula of N-soliton solutions is derived. It is shown that the Exp-function method may provide us with a straightforward and effective mathematical tool for generating N-soliton solutions of nonlinear evolution equations in mathematical physics.     

Exact solutions of semiclassical non-characteristic Cauchy problems for the sine-Gordon equation

The use of the sine-Gordon equation as a model of magnetic flux propagation in Josephson junctions motivates studying the initial-value problem for this equation in the semiclassical limit in which the dispersion parameter ε tends to zero. Assuming natural initial data having the profile of a moving −2π kink at time zero, we analytically calculate the scattering data of this completely integrable Cauchy problem for all ε>0 sufficiently small, and further we invert the scattering transform to calculate the solution for a sequence of arbitrarily small ε. This sequence of exact solutions is analogous to that of the well-known N-soliton (or higher-order soliton) solutions of the focusing nonlinear Schrödinger equation. We then use plots obtained from a careful numerical implementation of the inverse-scattering algorithm for reflectionless potentials to study the asymptotic behavior of solutions in the semiclassical limit. In the limit ε↓0 one observes the appearance of nonlinear caustics, i.e. curves in space-time that are independent of ε but vary with the initial data and that separate regions in which the solution is expected to have different numbers of nonlinear phases.In the appendices, we give a self-contained account of the Cauchy problem from the perspectives of both inverse scattering and classical analysis (Picard iteration). Specifically, Appendix A contains a complete formulation of the inverse-scattering method for generic L¹-Sobolev initial data, and Appendix B establishes the well-posedness for L^p-Sobolev initial data (which in particular completely justifies the inverse-scattering analysis in Appendix A).     

On the evaluation of a certain class of Feynman diagrams in x-space: Sunrise-type topologies at any loop order

We review recently developed new powerful techniques to compute a class of Feynman diagrams at any loop order, known as sunrise-type diagrams. These sunrise-type topologies have many important applications in many different fields of physics and we believe it to be timely to discuss their evaluation from a unified point of view. The method is based on the analysis of the diagrams directly in configuration space which, in the case of the sunrise-type diagrams and diagrams related to them, leads to enormous simplifications as compared to the traditional evaluation of loops in momentum space. We present explicit formulae for their analytical evaluation for arbitrary mass configurations and arbitrary dimensions at any loop order. We discuss several limiting cases in their kinematical regimes which are e.g. relevant for applications in HQET and NRQCD. We completely solve the problem of renormalization using simple formulae for the counterterms within dimensional regularization. An important application is the computation of the multi-particle phase space in D-dimensional space-time which we discuss. We present some examples of their numerical evaluation in the general case of D-dimensional space-time as well as in integer dimensions D = D₀ for different values of dimensions including the most important practical cases D₀ = 2, 3, 4. Substantial simplifications occur for odd integer space-time dimensions where the final results can be expressed in closed form through elementary functions. We discuss the use of recurrence relations naturally emerging in configuration space for the calculation of special series of integrals of the sunrise topology. We finally report on results for the computation of an extension of the basic sunrise topology, namely the spectacle topology and the topology where an irreducible loop is added.     

Perturbation theory and Padé approximants in realistic large-matrix models of the nuclear effective interaction

Realistic extended shell-model calculations are used to construct exact effective Hamiltonians, the perturbation series for the effective Hamiltonian to any order, and the [N+1, N] Padé approximants to the series. It is found that the Padé approximants give reliable results even when the series diverge, but that for both convergent and divergent series reasonably accurate results can be obtained only in fifth, or even seventh order. In addition, the poles of the low-order Padé approximants are not always reliable indicators of singularities of the perturbation series. The perturbation series and Padé approximants for the Q-box (energy-dependent effective Hamiltonian) are no more accurate in low orders than those for the usual effective Hamiltonian. Explicit formulas for the matrix Padé approximants are given in an appendix.     

A quantum exactly solvable non-linear oscillator with quasi-harmonic behaviour

The quantum version of a non-linear oscillator, previously analyzed at the classical level, is studied. This is a problem of quantization of a system with position-dependent mass of the form m = (1 + λx²)⁻¹ and with a λ-dependent non-polynomial rational potential. This λ-dependent system can be considered as a deformation of the harmonic oscillator in the sense that for λ → 0 all the characteristics of the linear oscillator are recovered. First, the λ-dependent Schrödinger equation is exactly solved as a Sturm-Liouville problem, and the λ-dependent eigenenergies and eigenfunctions are obtained for both λ > 0 and λ < 0. The λ-dependent wave functions appear as related with a family of orthogonal polynomials that can be considered as λ-deformations of the standard Hermite polynomials. In the second part, the λ-dependent Schrödinger equation is solved by using the Schrödinger factorization method, the theory of intertwined Hamiltonians, and the property of shape invariance as an approach. Finally, the new family of orthogonal polynomials is studied. We prove the existence of a λ-dependent Rodrigues formula, a generating function and λ-dependent recursion relations between polynomials of different orders.     

The inverse problem in the case of bound states

We investigate the inverse problem for bound states in the D = 3 dimensional space. The potential is assumed to be local and spherically symmetric. The present method is based on relationships connecting the moments of the ground state density to the lowest energy of each state of angular momentum ?. The reconstruction of the density ρ(r) from its moments is achieved by means of the series expansion of its Fourier transform F(q). The large q-behavior is described by Padé approximants. The accuracy of the solution depends on the number of known moments. The uniqueness is achieved if this number is infinite. In practice, however, an accuracy better than 1% is obtained with a set of about 15 levels.The method is tested on a simple example, and applied to three different spectra.     

Differential transform method for solving linear and non-linear systems of partial differential equations

In this Letter, we propose a reliable algorithm to develop exact and approximate solutions for the linear and non-linear systems of partial differential equations. The approach rest mainly on two-dimensional differential transform method which is one of the approximate methods. The method can easily be applied to many linear and non-linear problems and is capable of reducing the size of computational work. Exact solutions can also be achieved by the known forms of the series solutions. Several illustrative examples are given to demonstrate the effectiveness of the present method.     

Reduction of systems of first-order differential equations via Λ-symmetries

The notion of λ-symmetries, originally introduced by C. Muriel and J.L. Romero, is extended to the case of systems of first-order ODEs (and of dynamical systems in particular). It is shown that the existence of a symmetry of this type produces a reduction of the differential equations, restricting the presence of the variables involved in the problem. The results are compared with the case of standard (i.e. exact) Lie-point symmetries and are also illustrated by some examples.     

Generalized differential transform method to differential-difference equation

In this Letter, we generalize the differential transform method to solve differential-difference equation for the first time. Two simple but typical examples are applied to illustrate the validity and the great potential of the generalized differential transform method in solving differential-difference equation. A Padé technique is also introduced and combined with GDTM in aim of extending the convergence area of presented series solutions. Comparisons are made between the results of the proposed method and exact solutions. Then we apply the differential transform method to the discrete KdV equation and the discrete mKdV equation, and successfully obtain solitary wave solutions. The results reveal that the proposed method is very effective and simple. We should point out that generalized differential transform method is also easy to be applied to other nonlinear differential-difference equation.     

The q-gamma and (q,q)-polygamma functions of Tsallis statistics

An axiomatic definition is given for the q-gamma function of Tsallis (non-extensive) statistical physics, the continuous analogue of the q-factorial of Suyari [H. Suyari, Physica A 368 (1) (2006) 63], and the q-analogue of the gamma function of Euler and Gauss. A working definition in closed form, based on the Hurwitz and Riemann zeta functions (including their analytic continuations), is shown to satisfy this definition. Several relations involving the q-gamma and other functions are obtained. The (q,q)-polygamma functions , defined by successive derivatives of , where ln_qa=(1−q)⁻¹(a^1−q−1),a>0 is the q-logarithmic function, are also reported. The new functions are used to calculate the inferred probabilities and multipliers for Tsallis systems with finite numbers of particles N?∞.     

Delay-dependent exponential stability for neural networks with discrete and distributed time-varying delays

This Letter studies the exponential stability for a class of neural networks (NNs) with both discrete and distributed time-varying delays. Under weaker assumptions on the activation functions, by defining a more general type of Lyapunov functionals and developing a new convex combination technique, new less conservative and less complex stability criteria are established to guarantee the global exponential stability of the discussed NNs. The obtained conditions are dependent on both discrete and distributed delays, are expressed in terms of linear matrix inequalities (LMIs), and contain fewer decision variables. Numerical examples are given to illustrate the effectiveness and the less conservatism of the proposed conditions.     

Approximate Homotopy Direct Reduction Method: Infinite Series Reductions to Perturbed mKdV Equations

An approximate homotopy direct reduction method is proposed and applied to two perturbed modified Korteweg- de Vries （mKdV） equations with fourth-order dispersion and second-order dissipation. The similarity reduction equations are derived to arbitrary orders. The method is valid not only for single soliton solutions but also for the Painlevd Ⅱ waves and periodic waves expressed by Jacobi elliptic functions for both fourth-order dispersion and second-order dissipation. The method is also valid for strong perturbations.     

The nature of the Extended Analog Computer

This paper defines the relationship between Rubel’s Extended Analog Computer (EAC) model and Indiana University’s implementation of it, a fundamental distinction that is necessary to understand the design, implementation, operation, and applications of the EAC. Beginning with the detailed structure of the current EAC, its method of computing is illustrated with diagrams and three applications. Next, the Δ-digraph, a labeled directed graph, is introduced to show how unconventional and conventional computers relate nature, mathematics and computer architecture. The Δ-digraph defines the paradigms of analogy and algorithm, illustrates how applications for the EAC are analogies developed by choosing the semantics for a machine configuration, and suggests how partial differential equations might be compiled to EAC configurations (a difficult problem compared to compiling a digital computer program, and one that is still not solved). The Δ-digraph also distinguishes between the EAC’s explicit functions (whose operation is defined in the structure of its components) and implicit functions (whose operation is inherent in the properties of matter and energy). Brief case studies show how industry and academia are adopting the EAC. The paper closes with topics and open questions for future research.     

Accurate summation of the perturbation series for periodic eigenvalue problems

We show that algebraic approximants prove suitable for the summation of the perturbation series for the eigenvalues of periodic problems. Appropriate algebraic approximants constructed from the perturbation series for a given eigenvalue provide information about other eigenvalues connected with the chosen one by branch points in the complex plane. Such approximants also give those branch points with remarkable accuracy. We choose Mathieu's equation as illustrative example. Received 6 December 2000