Source author record

J. Segura

J. Segura appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2022arXiv

Complete asymptotic expansions for the relativistic Fermi-Dirac integral

Fermi-Dirac integrals appear in problems in nuclear astrophysics, solid state physics or in the fundamental theory of semiconductor modeling, among others areas of application. In this paper, we give new and complete asymptotic expansions for the relativistic Fermi-Dirac integral. These expansions could be useful to obtain a correct qualitative understanding of Fermi systems. The performance of the expansions is illustrated with numerical examples.

preprint2020arXiv

Asymptotic computation of classical orthogonal polynomials

The classical orthogonal polynomials (Hermite, Laguerre and Jacobi) are involved in a vast number of applications in physics and engineering. When large degrees $n$ are needed, the use of recursion to compute the polynomials is not a good strategy for computation and a more efficient approach, such as the use of asymptotic expansions,is recommended. In this paper, we give an overview of the asymptotic expansions considered in [8] for computing Laguerre polynomials $L^{(α)}_n(x)$ for bounded values of the parameter $α$. Additionally, we show examples of the computational performance of an asymptotic expansion for $L^{(α)}_n(x)$ valid for large values of $α$ and $n$. This expansion was used in [6] as starting point for obtaining asymptotic approximations to the zeros. Finally, we analyze the expansions considered in [9], [10] and [11] to compute the Jacobi polynomials for large degrees $n$.

preprint2020arXiv

Asymptotic inversion of the binomial and negative binomial cumulative distribution functions

The computation and inversion of the binomial and negative binomial cumulative distribution functions play a key role in many applications. In this paper, we explain how methods used for the central beta distribution function (described in [2]) can be used to obtain asymptotic representations of these functions, and also for their inversion. The performance of the asymptotic inversion methods is illustrated with numerical examples.

preprint2020arXiv

Fast and reliable high accuracy computation of Gauss--Jacobi quadrature

Iterative methods with certified convergence for the computation of Gauss--Jacobi quadratures are described. The methods do not require a priori estimations of the nodes to guarantee its fourth-order convergence. They are shown to be generally faster than previous methods and without practical restrictions on the range of the parameters. The evaluation of the nodes and weights of the quadrature is exclusively based on convergent processes which, together with the fourth order convergence of the fixed point method for computing the nodes, makes this an ideal approach for high accuracy computations, so much so that computations of quadrature rules with even millions of nodes and thousands of digits are possible in a typical laptop.

preprint2020arXiv

Simplified error bounds for turning point expansions

Recently, the present authors derived new asymptotic expansions for linear differential equations having a simple turning point. These involve Airy functions and slowly varying coefficient functions, and were simpler than previous approximations, in particular being computable to a high degree of accuracy. Here we present explicit error bounds for these expansions which only involve elementary functions, and thereby provide a simplification of the bounds associated with the classical expansions of F. W. J. Olver.

preprint2016arXiv

A new type of sharp bounds for ratios of modified Bessel functions

The bounds for the ratios of first and second kind modified Bessel functions of consecutive orders are important quantities appearing in a large number of scientific applications. We obtain new bounds which are accurate in a large region of parameters and which are shaper than previous bounds. The new bounds are obtained by a qualitative analysis of the Riccati equation satisfied by these ratios. A procedure is considered in which the bounds obtained from the analysis of the Riccati equation are used to define a new function satisfying a new Riccati equation which yields sharper bounds. Similar ideas can be applied to other functions.

preprint2016arXiv

Computation of the incomplete gamma function for negative values of the argument

An algorithm for computing the incomplete gamma function $γ^*(a,z)$ for real values of the parameter $a$ and negative real values of the argument $z$ is presented. The algorithm combines the use of series expansions, Poincaré-type expansions, uniform asymptotic expansions and recurrence relations, depending on the parameter region. A relative accuracy $\sim 10^{-13}$ in the parameter region $(a,z) \in [-500,\,500] \times [-500,\,0)$ can be obtained when computing the function $γ^*(a,z)$ with the Fortran 90 module IncgamNEG implementing the algorithm.

preprint2015arXiv

Computing the Kummer function U(a,b,z) for small values of the arguments

We describe methods for computing the Kummer function $U(a,b,z)$ for small values of $z$, with special attention to small values of $b$. For these values of $b$ the connection formula that represents $U(a,b,z)$ as a linear combination of two ${}_1F_1$-functions needs a limiting procedure. We use the power series of the ${}_1F_1$-functions and consider the terms for which this limiting procedure is needed. We give recursion relations for higher terms in the expansion, and we consider the derivative $U^\prime(a,b,z)$ as well. We also discuss the performance for small $\vert z\vert$ of an asymptotic approximation of the Kummer function in terms of modified Bessel functions.

preprint2015arXiv

GammaCHI: a package for the inversion and computation of the gamma and chi-square cumulative distribution functions (central and noncentral)

A Fortran 90 module (GammaCHI) for computing and inverting the gamma and chi-square cumulative distribution functions (central and noncentral) is presented. The main novelty of this package are the reliable and accurate inversion routines for the noncentral cumulative distribution functions. Additionally, the package also provides routines for computing the gamma function, the error function and other functions related to the gamma function. The module includes the routines cdfgamC, invcdfgamC, cdfgamNC, invcdfgamNC, errorfunction, inverfc, gamma, loggam, gamstar and quotgamm for the computation of the central gamma distribution function (and its complementary function), the inversion of the central gamma distribution function, the computation of the noncentral gamma distribution function (and its complementary function), the inversion of the noncentral gamma distribution function, the computation of the error function and its complementary function, the inversion of the complementary error function, the computation of: the gamma function, the logarithm of the gamma function, the regulated gamma function and the ratio of two gamma functions, respectively.

preprint2014arXiv

Computation of a numerically satisfactory pair of solutions of the differential equation for conical functions of non-negative integer orders

We consider the problem of computing satisfactory pairs of solutions of the differential equation for Legendre functions of non-negative integer order $μ$ and degree $-\frac12+iτ$, where $τ$ is a non-negative real parameter. Solutions of this equation are the conical functions ${\rm{P}}^μ_{-\frac12+iτ}(x)$ and ${Q}^μ_{-\frac12+iτ}(x)$, $x>-1$. An algorithm for computing a numerically satisfactory pair of solutions is already available when $-1<x<1$ (see \cite{gil:2009:con}, \cite{gil:2012:cpc}).In this paper, we present a stable computational scheme for a real valued numerically satisfactory companion of the function ${\rm{P}}^μ_{-\frac12+iτ}(x)$ for $x>1$, the function $\Re\left\{e^{-iπμ} {Q}^μ_{-\frac{1}{2}+iτ}(x) \right\}$. The proposed algorithm allows the computation of the function on a large parameter domain without requiring the use of extended precision arithmetic.

preprint2014arXiv

Monotonicity properties and bounds for the chi-square and gamma distributions

The generalized Marcum functions $Q_μ(x,y)$ and $P_μ(x,y)$ have as particular cases the non-central $χ^2$ and gamma cumulative distributions, which become central distributions (incomplete gamma function ratios) when the non-centrality parameter $x$ is set to zero. We analyze monotonicity and convexity properties for the generalized Marcum functions and for ratios of Marcum functions of consecutive parameters (differing in one unity) and we obtain upper and lower bounds for the Marcum functions. These bounds are proven to be sharper than previous estimations for a wide range of the parameters. Additionally we show how to build convergent sequences of upper and lower bounds. The particularization to incomplete gamma functions, together with some additional bounds obtained for this particular case, lead to combined bounds which improve previously exiting inequalities.

preprint2014arXiv

Recent software developments for special functions in the Santander-Amsterdam project

We give an overview of published algorithms by our group and of current activities and future plans. In particular, we give details on methods for computing special functions and discuss in detail two current lines of research. Firstly, we describe the recent developments for the computation of central and non-central chi-square cumulative distributions (also called Marcum Q-functions), and we present a new quadrature method for computing them. Secondly, we describe the fourth-order methods for computing zeros of special functions recently developed, and we provide an explicit example for the computation of complex zeros of Bessel functions. We end with an overview of published software by our group for computing special functions.

preprint2014arXiv

The asymptotic and numerical inversion of the Marcum $Q-$function

The generalized Marcum functions appear in problems of technical and scientific areas such as, for example, radar detection and communications. In mathematical statistics and probability theory these functions are called the noncentral gamma or the noncentral chi-squared cumulative distribution functions. In this paper we describe a new asymptotic method for inverting the generalized Marcum $Q-$function and for the complementary Marcum $P-$function. Also, we show how monotonicity and convexity properties of these functions can be used to find initial values for reliable Newton or secant methods to invert the function. We present details of numerical computations that show the reliability of the asymptotic approximations.

preprint2013arXiv

Computation of the Marcum Q-function

Methods and an algorithm for computing the generalized Marcum $Q-$function ($Q_μ(x,y)$) and the complementary function ($P_μ(x,y)$) are described. These functions appear in problems of different technical and scientific areas such as, for example, radar detection and communications, statistics and probability theory, where they are called the non-central chi-square or the non central gamma cumulative distribution functions. The algorithm for computing the Marcum functions combines different methods of evaluation in different regions: series expansions, integral representations, asymptotic expansions, and use of three-term homogeneous recurrence relations. A relative accuracy close to $10^{-12}$ can be obtained in the parameter region $(x,y,μ) \in [0,\,A]\times [0,\,A]\times [1,\,A]$, $A=200$, while for larger parameters the accuracy decreases (close to $10^{-11}$ for $A=1000$ and close to $5\times 10^{-11}$ for $A=10000$).

preprint2007arXiv

A Fast Concurrent Power-Thermal Model for Sub-100nm Digital ICs

As technology scales down, the static power is expected to become a significant fraction of the total power. The exponential dependence of static power with the operating temperature makes the thermal profile estimation of high-performance ICs a key issue to compute the total power dissipated in next-generations. In this paper we present accurate and compact analytical models to estimate the static power dissipation and the temperature of operation of CMOS gates. The models are the fundamentals of a performance estimation tool in which numerical procedures are avoided for any computation to set a faster estimation and optimization. The models developed are compared to measurements and SPICE simulations for a 0.12mm technology showing excellent results.

preprint2007arXiv

Smart Temperature Sensor for Thermal Testing of Cell-Based ICs

In this paper we present a simple and efficient built-in temperature sensor for thermal monitoring of standard-cell based VLSI circuits. The proposed smart temperature sensor uses a ring-oscillator composed of complex gates instead of inverters to optimize their linearity. Simulation results from a 0.18$μ$m CMOS technology show that the non-linearity error of the sensor can be reduced when an adequate set of standard logic gates is selected.

preprint1997arXiv

Neutrino Oscillation and Magnetic Moment from $ν- e^{-}$ Elastic Scattering

We discuss how the measurement of the $\barν_{e}-e^{-}$ elastic cross section at reactor energies can be used to extract new information on the neutrino oscillation parameters. We also consider the magnetic moment contribution and show how both effects tend to cancel each other when the total cross section is measured; to achieve the separation of each of the effects, experiments capable of measuring angular and energy distributions with respect to the outgoing electron become necessary. We study how their different energy and angle dependence enables such a separation; then, the sensitivity to magnetic moments, masses and mixings is discussed. We also show how these experiments can be sensitive to the magnetic moment of $τ$ neutrinos via $\barν_{e}\leftrightarrow \barν_μ$ oscillation.

preprint1994arXiv

A Novel Kind of Neutrino Oscillation Experiment

A novel method to look for neutrino oscillations is proposed based on the elastic scattering process $\barν_{i} e^{-}\rightarrow \barν_{i} e^{-}$, taking advantage of the dynamical zero present in the differential cross section for $\barν_{e} e^{-}\rightarrow \barν_{e} e^{-}$. An effective tunable experiment between the "appearance" and "disappearance" limits is made possible. Prospects to exclude the allowed region for atmospheric neutrino oscillations are given.

preprint1994arXiv

Neutrino Magnetic Moment and the Process $νe \rightarrow νe γ$

The contribution of a neutrino magnetic moment $μ_ν$ to the cross section of the process $νe^{-}\rightarrow νe^{-} γ$ has been calculated and compared with the Standard Electroweak one. The radiative process allows to reach low enough values of $Q^2$ without the need to operate at very small energies of recoil electrons. Regions in the phase space which are more favourable to set bounds on $μ_ν$ are suggested.

preprint1994arXiv

Neutrino Mass and Magnetic Moment from Neutrino-Electron Scattering

We study both the elastic ($νe \rightarrow νe$) and the radiative process ($νe \rightarrow νe γ$) and discuss how these processes can shed light on some current topics in neutrino physics such as a neutrino magnetic moment and neutrino oscillations. The radiative process allows to reach low values of $Q^2$ without the need to operate at very small energies of recoil electrons, a favourable scenario to search for a neutrino magnetic moment. The elastic cross section contains a dynamical zero at $E_ν=m/(4 sin^2θ_{W})$ and forward electrons for the electron antineutrino channel, which is reachable at reactor facilities and accessible after the convolution with the antineutrino spectrum. The implication for lepton flavour changing transitions in that energy region searched for in neutrino oscillation experiments, which combine disappearance and appearance rates, is discussed.

preprint1993arXiv

Dynamical zeros in neutrino-electron elastic scattering at leading order

We show the existence of dynamical zeros in the helicity amplitudes for neutrino-electron elastic scattering at lowest order in the standard theory. In particular, the $λ=1/2$ non-flip electron helicity amplitude in the electron antineutrino process vanishes for an incident neutrino energy $E_ν=m_{e}/(4sin^{2}θ_{W})$ and forward electrons (maximum recoil energy). The rest of helicity amplitudes show kinematical zeros in this configuration and therefore the cross section vanishes. Prospects to search for neutrino magnetic moment are discussed.