Open Access
An asymptotic expansion of the hyberbolic umbilic catastrophe integral
(Springer, 2022) Ferreira González, Chelo; López García, José Luis; Pérez Sinusía, Ester; Estadística, Informática y Matemáticas; Estatistika, Informatika eta Matematika; Institute for Advanced Materials and Mathematics - INAMAT2; Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa
We obtain an asymptotic expansion of the hyperbolic umbilic catastrophe integral Ψ(H) (x,y,z) := ∫∞−∞∫∞−∞exp(i(s3+t3+zst +yt+xs))ds dt for large values of |x| and bounded values of |y| and |z|. The expansion is given in terms of Airy functions and inverse powers of x. There is only one Stokes ray at argx=π . We use the modified saddle point method introduced in (López et al. J Math Anal Appl 354(1):347–359, 2009). The accuracy and the asymptotic character of the approximations are illustrated with numerical experiments.
Open Access
Analytical ray transfer matrix for the crystalline lens
(OSA Publishing, 2022) Navarro, Rafael; Lockett-Ruiz, Verónica; López García, José Luis; Estadística, Informática y Matemáticas; Estatistika, Informatika eta Matematika
We present the formulation of a paraxial ray transfer or ABCD matrix for onion-type GRIN lenses. In GRIN lenses, each iso-indicial surface (IIS) can be considered a refracting optical surface. If each IIS is a shell or layer, the ABCD matrix of a GRIN lens is computed by multiplying a typically high number of translation and refraction matrices corresponding to the K layers inside the lens. Using a differential approximation for the layer thickness, this matrix product becomes a sum. The elements A, B, C, and D of the approximated GRIN ray transfer matrix can be calculated by integrating the elements of a single-layer matrix. This ABCD matrix differs from a homogeneous lens matrix in only one integration term in element C, corresponding to the GRIN contribution to the lens power. Thus the total GRIN lens power is the sum of the homogeneous lens power and the GRIN contribution, which offers a compact and simple expression for the ABDC matrix. We then apply this formulation to the crystalline lens and implement both numerical and analytical integration procedures to obtain the GRIN lens power. The analytical approximation provides an accurate solution in terms of Gaussian hypergeometric functions. Last, we compare our numerical and analytical procedures with published ABCD matrix methods in the literature, and analyze the effect of the iso-indicial surface’s conic constant (Q) and inner curvature gradient (G) on the lens power for different lens models.
Open Access
Convergent and asymptotic expansions of solutions of differential equations with a large parameter: Olver cases II and III
(Rocky Mountain Mathematics Consortium, 2015) Ferreira González, Chelo; López García, José Luis; Pérez Sinusía, Ester; Ingeniería Matemática e Informática; Matematika eta Informatika Ingeniaritza
This paper continues the investigation initiated in [Lopez, 2013]. We consider the asymptotic method designed by F. Olver [Olver, 1974] for linear differential equations of the second order containing a large (asymptotic) parameter . We consider here the second and third cases studied by Olver: differential equations with a turning point (second case) or a singular point (third case). It is well-known that his method gives the Poincar´e-type asymptotic expansion of two independent solutions of the equation in inverse powers of . In this paper we add initial conditions to the differential equation and consider the corresponding initial value problem. By using the Green function of an auxiliary problem, we transform the initial value problem into a Volterra integral equation of the second kind. Then, using a fixed point theorem, we construct a sequence of functions that converges to the unique solution of the problem. This sequence has also the property of being an asymptotic expansion for large (not of Poincar´e-type) of the solution of the problem. Moreover, we show
Open Access
Convergent expansions of the Bessel functions in terms of elementary functions
(Springer US, 2018) López García, José Luis; Pagola Martínez, Pedro Jesús; Matematika eta Informatika Ingeniaritza; Institute for Advanced Materials and Mathematics - INAMAT2; Ingeniería Matemática e Informática; Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa
We consider the Bessel functions Jν (z) and Yν (z) for ν > −1/2 and z ≥ 0. We derive a convergent expansion of Jν (z) in terms of the derivatives of (sin z)/z, and a convergent expansion of Yν (z) in terms of derivatives of (1−cos z)/z, derivatives of (1 − e−z)/z and (2ν, z). Both expansions hold uniformly in z in any fixed horizontal strip and are accompanied by error bounds. The accuracy of the approximations is illustrated with some numerical experiments.
Open Access
An extension of the multiple Erdélyi-Kober operator and representations of the generalized hypergeometric functions
(De Gruyter, 2018) Karp, D. B.; López García, José Luis; Ingeniería Matemática e Informática; Matematika eta Informatika Ingeniaritza; Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa
In this paper we investigate the extension of the multiple Erd elyi-Kober fractional integral operator of Kiryakova to arbitrary complex values of parameters by the way of regularization. The regularization involves derivatives of the function in question and the integration with respect to a kernel expressed in terms of special case of Meijer's G function. An action of the regularized multiple Erd elyi-Kober operator on some simple kernels leads to decomposition formulas for the generalized hypergeometric functions. In the ultimate section, we de ne an alternative regularization better suited for representing the Bessel type generalized hypergeometric function p􀀀1Fp. A particular case of this regularization is then used to identify some new facts about the positivity and reality of zeros of this function.
Open Access
Convergent and asymptotic methods for second-order difference equations with a large parameter
(Springer, 2018) Ferreira González, Chelo; López García, José Luis; Pérez Sinusía, Ester; Matematika eta Informatika Ingeniaritza; Institute for Advanced Materials and Mathematics - INAMAT2; Ingeniería Matemática e Informática; Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa
We consider the second-order linear difference equation y(n+2)−2ay(n+1)−Λ2y(n)=g(n)y(n)+f(n)y(n+1) , where Λ is a large complex parameter, a≥0 and g and f are sequences of complex numbers. Two methods are proposed to find the asymptotic behavior for large |Λ|of the solutions of this equation: (i) an iterative method based on a fixed point method and (ii) a discrete version of Olver’s method for second-order linear differential equations. Both methods provide an asymptotic expansion of every solution of this equation. The expansion given by the first method is also convergent and may be applied to nonlinear problems. Bounds for the remainders are also given. We illustrate the accuracy of both methods for the modified Bessel functions and the associated Legendre functions of the first kind.
Open Access
Series representations of the Volterra function and the Fransén–Robinson constant
(Elsevier, 2021) López García, José Luis; Pagola Martínez, Pedro Jesús; Palacios Herrero, Pablo; Estatistika, Informatika eta Matematika; Institute for Advanced Materials and Mathematics - INAMAT2; Estadística, Informática y Matemáticas; Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa
The Volterra function μ(t,β,α) was introduced by Vito Volterra in 1916 as the solution to certain integral equations with a logarithmic kernel. Despite the large number of applications of the Volterra function, the only known analytic representations of this function are given in terms of integrals. In this paper we derive several convergent expansion of μ(t,β,α) in terms of incomplete gamma functions. These expansions may be used to implement numerical evaluation techniques for this function. As a particular application, we derive a numerical series representation of the Fransén–Robinson constant F := µ(1, 1, 0) = R ∞ 0 1 Γ(x) dx. Some numerical examples illustrate the accuracy of the approximations
Open Access
Effect of high-energy ball-milling on the magnetostructural properties of a Ni45Co5Mn35Sn15 alloy
(Elsevier, 2021) López García, José Luis; Sánchez-Alarcos Gómez, Vicente; Recarte Callado, Vicente; Rodríguez Velamazán, José Alberto; Unzueta, Iraultza; García, José Ángel; Plazaola, Fernando; La Roca, Paulo Matías; Pérez de Landazábal Berganzo, José Ignacio; Zientziak; Institute for Advanced Materials and Mathematics - INAMAT2; Ciencias; Gobierno de Navarra / Nafarroako Gobernua, PC017-018 AMELEC
The effect of high-energy ball-milling on the magnetostructural properties of a Ni45Co5Mn35Sn15 alloy in austenitic phase at room temperature has been analyzed by neutron and high-resolution X-ray diffraction. The ball milling promotes a mechanically-induced martensitic transformation as well as the appearance of amorphous-like non-transforming regions, following a double stage; for short milling times (below 30 min), a strong size reduction and martensite induction occur. On the opposite, for longer times, the increase of strains predominates and consequently a larger amount of non-transforming regions appears. The effect of the microstructural defects brought by milling (as dislocations) on both the enthalpy change at the martensitic transformation and the high field magnetization of the austenite has been quantitatively estimated and correlated to the internal strains. Contrary to what occurs in ternary Ni-Mn-Sn alloys, the mechanically-induced defects do not change the ferromagnetic coupling between Mn atoms, but just cause a net reduction on the magnetic moments.
Open Access
Uniform convergent expansions of integral transforms
(American Mathematical Society, 2021) López García, José Luis; Palacios Herrero, Pablo; Pagola Martínez, Pedro Jesús; Estadística, Informática y Matemáticas; Estatistika, Informatika eta Matematika; Institute for Advanced Materials and Mathematics - INAMAT2; Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa
Several convergent expansions are available for most of the special functions of the mathematical physics, as well as some asymptotic expansions [NIST Handbook of Mathematical Functions, 2010]. Usually, both type of expansions are given in terms of elementary functions; the convergent expansions provide a good approximation for small values of a certain variable, whereas the asymptotic expansions provide a good approximation for large values of that variable. Also, quite often, those expansions are not uniform: the convergent expansions fail for large values of the variable and the asymptotic expansions fail for small values. In recent papers [Bujanda & all, 2018-2019] we have designed new expansions of certain special functions, given in terms of elementary functions, that are uniform in certain variables, providing good approximations of those special functions in large regions of the variables, in particular for large and small values of the variables. The technique used in [Bujanda & all, 2018-2019] is based on a suitable integral representation of the special function. In this paper we face the problem of designing a general theory of uniform approximations of special functions based on their integral representations. Then, we consider the following integral transform of a function g(t) with kernel h(t, z), F(z) := 1 0 h(t, z)g(t)dt. We require for the function h(t, z) to be uniformly bounded for z ∈D⊂ C by a function H(t) integrable in t ∈ [0, 1], and for the function g(t) to be analytic in an open region Ω that contains the open interval (0, 1). Then, we derive expansions of F(z) in terms of the moments of the function h, M[h(·, z), n] := 1 0 h(t, z)tndt, that are uniformly convergent for z ∈ D. The convergence of the expansion is of exponential order O(a−n), a > 1, when [0, 1] ∈ Ω and of power order O(n−b), b > 0, when [0, 1] ∈/ Ω. Most of the special functions F(z) having an integral representation may be cast in this form, possibly after an appropriate change of the integration variable. Then, special interest has the case when the moments M[h(·, z), n] are elementary functions of z, because in that case the uniformly convergent expansion derived for F(z) is given in terms of elementary functions. We illustrate the theory with several examples of special functions different from those considered in [Bujanda & all, 2018-2019].
Open Access
A simplification of the stationary phase method: application to the Anger and Weber functions
(Kent State University, 2017) López García, José Luis; Matematika eta Informatika Ingeniaritza; Institute for Advanced Materials and Mathematics - INAMAT2; Ingeniería Matemática e Informática; Universidad Pública de Navarra / Nafarroako Unibertsitate Publikoa
The main difficulty in the practical use of the stationary phase method in asymptotic expansions of integrals is originated by a change of variables. The coefficients of the asymptotic expansion are the coefficients of the Taylor expansion of a certain function implicitly defined by that change of variables. In general, this function is not explicitly known, and then the computation of those coefficients is cumbersome. Using the factorization of the exponential factor used in previous works of [Tricomi, 1950], [Erdélyi and Wyman, 1963], and [Dingle, 1973], we obtain a variant of the method that avoids that change of variables and simplifies the computations. On the one hand, the calculation of the coefficients of the asymptotic expansion is remarkably simpler and explicit. On the other hand, the asymptotic sequence is as simple as in the standard stationary phase method: inverse powers of the asymptotic variable. New asymptotic expansions of the Anger and Weber functions Jλx(x) and Eλx(x) for large positive x and real parameter λ 6= 0 are given as an illustration.