
Asymptotic approximations of Jacobi polynomials are given for large values of
the $\beta$parameter and of their zeros. The expansions are given in terms of
Laguerre polynomials and of their zeros. The levels of accuracy of the
approximations are verified by numerical examples.

Asymptotic approximations to the zeros of Hermite and Laguerre polynomials
are given, together with methods for obtaining the coefficients in the
expansions. These approximations can be used as a standalone method of
computation of Gaussian quadratures for high enough degrees, with Gaussian
weights computed from asymptotic approximations for the orthogonal polynomials.
We provide numerical evidence showing that for degrees greater than $100$ the
asymptotic methods are enough for a double precision accuracy computation
($15$$16$ digits) of the nodes and weights of the GaussHermite and
GaussLaguerre quadratures.

Uniform asymptotic expansions involving exponential and Airy functions are
obtained for Laguerre polynomials $L_{n}^{(\alpha)}(x)$, as well as
complementary confluent hypergeometric functions. The expansions are valid for
$n$ large and $\alpha$ small or large, uniformly for unbounded real and complex
values of $x$. The new expansions extend the range of computability of
$L_n^{(\alpha)}(x)$ compared to previous expansions, in particular with respect
to higher terms and large values of $\alpha$. Numerical evidence of their
accuracy for real and complex values of $x$ is provided.

Conical functions appear in a large number of applications in physics and
engineering. In this paper we describe an extension of our module CONICAL for
the computation of conical functions. Specifically, the module includes now a
routine for computing the function ${{\rm R}}^{m}_{\frac{1}{2}+i\tau}(x)$, a
realvalued numerically satisfactory companion of the function ${\rm
P}^m_{\tfrac12+i\tau}(x)$ for $x>1$. In this way, a natural basis for solving
Dirichlet problems bounded by conical domains is provided.

Linear second order differential equations having a large real parameter and
turning point in the complex plane are considered. Classical asymptotic
expansions for solutions involve the Airy function and its derivative, along
with two infinite series, the coefficients of which are usually difficult to
compute. By considering the series as asymptotic expansions for two explicitly
defined analytic functions, Cauchy's integral formula is employed to compute
the coefficient functions to high order of accuracy. The method employs a
certain exponential form of LiouvilleGreen expansions for solutions of the
differential equation, as well as for the Airy function. We illustrate the use
of the method with the high accuracy computation of Airytype expansions of
Bessel functions of complex argument.

Nearly 70 Ground Level Enhancements (GLEs) of cosmic rays have been recorded
by the worldwide neutron monitor network since the 1950s depicting a big
variety of energy spectra of solar energetic particles (SEP). Here we studied a
statistical relation between the eventintegrated intensity of GLEs (calculated
as countrate relative excess, averaged over all available polar neutron
monitors, and expressed in percenthours) and the hardness of the solar
particle energy spectra. For each event the integral omnidirectional
eventintegrated fluences of particles with energy above 30 MeV ($F_{30}$) and
above 200 MeV ($F_{200}$) were computed using the reconstructed spectra, and
the ratio between the two fluences was considered as a simple index of the
event's hardness. We also provided a justification of the spectrum estimate in
the form of the Bandfunction, using direct PAMELA data for GLE 71
(17May2012). We found that, while there is no clear relation between the
intensity and the hardness for weak events, all strong events with the
intensity greater 100 \%*hr are characterized by a very hard spectrum. This
implies that a hard spectrum can be securely assumed for all extreme GLE
events, e.g., those studied using cosmogenic isotope data in the past.

An efficient algorithm and a Fortran 90 module (LaguerrePol) for computing
Laguerre polynomials $L^{(\alpha)}_n(z)$ are presented. The standard threeterm
recurrence relation satisfied by the polynomials and different types of
asymptotic expansions valid for $n$ large and $\alpha$ small, are used
depending on the parameter region.
Based on tests of contiguous relations in the parameter $\alpha$ and the
degree $n$ satisfied by the polynomials, we claim that a relative accuracy
close or better than $10^{12}$ can be obtained using the module LaguerrePol
for computing the functions $L^{(\alpha)}_n(z)$ in the parameter range $z \ge
0$, $1 < \alpha \le 5$, $n \ge 0$.

An algorithm for computing the incomplete gamma function $\gamma^*(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\'etype 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 $\gamma^*(a,z)$ with the Fortran 90
module IncgamNEG implementing the algorithm.

Accurate and efficient algorithms for the inversion of the cumulative central
beta distribution are described. The algorithms are based on the combination of
a fourthorder fixed point method with good nonlocal convergence properties
(the SchwarzianNewton method), asymptotic inversion methods and sharp bounds
in the tails of the distribution function.

The tagged quasifree $np \to np\pi^+\pi^$ reaction has been studied
experimentally with the High Acceptance DiElectron Spectrometer (HADES) at GSI
at a deuteron incident beam energy of 1.25 GeV/nucleon ($\sqrt s \sim$ 2.42
GeV/c for the quasifree collision). For the first time, differential
distributions for $\pi^{+}\pi^{}$ production in $np$ collisions have been
collected in the region corresponding to the large transverse momenta of the
secondary particles. The invariant mass and angular distributions for the
$np\rightarrow np\pi^{+}\pi^{}$ reaction are compared with different models.
This comparison confirms the dominance of the $t$channel with $\Delta\Delta$
contribution. It also validates the changes previously introduced in the
Valencia model to describe twopion production data in other isospin channels,
although some deviations are observed, especially for the $\pi^{+}\pi^{}$
invariant mass spectrum. The extracted total cross section is also in much
better agreement with this model. Our new measurement puts useful constraints
for the existence of the conjectured dibaryon resonance at mass M$\sim$ 2.38
GeV and with width $\Gamma\sim$ 70 MeV.

We present the stochastic model of the galactic cosmic ray (GCR) particles
transport in the heliosphere. Based on the solution of the Parker transport
equation we developed models of the shorttime variation of the GCR intensity,
i.e. the Forbush decrease (Fd) and the 27day variation of the GCR intensity.
Parker transport equation being the FokkerPlanck type equation delineates
nonstationary transport of charged particles in the turbulent medium. The
presented approach of the numerical solution is grounded on solving of the set
of equivalent stochastic differential equations (SDEs). We demonstrate the
method of deriving from Parker transport equation the corresponding SDEs in the
heliocentric spherical coordinate system for the backward approach. Features
indicative the preeminence of the backward approach over the forward is
stressed. We compare the outcomes of the stochastic model of the Fd and 27day
variation of the GCR intensity with our former models established by the finite
difference method. Both models are in an agreement with the experimental data.

We present the newly developed stochastic model of the galactic cosmic ray
(GCR) particles transport in the heliosphere. Mathematically Parker transport
equation (PTE) describing nonstationary transport of charged particles in the
turbulent medium is the FokkerPlanck type. It is the second order parabolic
timedependent 4dimensional (3 spatial coordinates and particles
energy/rigidity) partial differential equation. It is worth to mention that, if
we assume the stationary case it remains as the 3D parabolic type problem with
respect to the particles rigidity R. If we fix the energy it still remains as
the 3D parabolic type problem with respect to time. The proposed method of
numerical solution is based on the solution of the system of stochastic
differential equations (SDEs) being equivalent to the Parker's transport
equation. We present the method of deriving from PTE the equivalent SDEs in the
heliocentric spherical coordinate system for the backward approach. The
obtained stochastic model of the Forbush decrease of the GCR intensity is in an
agreement with the experimental data. The advantages and disadvantages of the
forward and the backward solution of the PTE are discussed.

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.

A Fortran 90 module (GammaCHI) for computing and inverting the gamma and
chisquare 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.

Ionization and scintillation produced by nuclear recoils in gaseous xenon at
approximately 14 bar have been simultaneously observed in an electroluminescent
time projection chamber. Neutrons from radioisotope $\alpha$Be neutron sources
were used to induce xenon nuclear recoils, and the observed recoil spectra were
compared to a detailed Monte Carlo employing estimated ionization and
scintillation yields for nuclear recoils. The ability to discriminate between
electronic and nuclear recoils using the ratio of ionization to primary
scintillation is demonstrated. These results encourage further investigation on
the use of xenon in the gas phase as a detector medium in dark matter direct
detection experiments.

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 chisquared 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.

We study the distribution of zeros of general solutions of the Airy and
Bessel equations in the complex plane. Our results characterize the patterns
followed by the zeros for any solution, in such a way that if one zero is known
it is possible to determine the location of the rest of zeros.

We consider the problem of computing satisfactory pairs of solutions of the
differential equation for Legendre functions of nonnegative integer order
$\mu$ and degree $\frac12+i\tau$, where $\tau$ is a nonnegative real
parameter. Solutions of this equation are the conical functions
${\rm{P}}^{\mu}_{\frac12+i\tau}(x)$ and ${Q}^{\mu}_{\frac12+i\tau}(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}}^{\mu}_{\frac12+i\tau}(x)$ for $x>1$, the function
$\Re\left\{e^{i\pi \mu} {{Q}}^{\mu}_{\frac{1}{2}+i\tau}(x) \right\}$. The
proposed algorithm allows the computation of the function on a large parameter
domain without requiring the use of extended precision arithmetic.

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 noncentral chisquare cumulative distributions (also called Marcum
Qfunctions), and we present a new quadrature method for computing them.
Secondly, we describe the fourthorder 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.

A technical description of NEXTMM and its commissioning and first
performance is reported. Having an active volume of ~35 cm drift $\times$ 28 cm
diameter, it constitutes the largest Micromegasread TPC operated in Xenon ever
constructed, made by a sectorial arrangement of the 4 largest single wafers
manufactured with the Microbulk technique to date. It is equipped with a
suitably pixelized readout and with a sufficiently large sensitive volume (~23
l) so as to contain long (~20 cm) electron tracks. First results obtained at 1
bar for Xenon and trimethylamine (Xe(2 %)TMA) mixture are presented. The TPC
can accurately reconstruct extended background tracks. An encouraging
fullwidth halfmaximum of 11.6 % was obtained for ~29 keV gammas without
resorting to any data postprocessing.

NEXTMM is a generalpurpose high pressure (10 bar, $\sim25$ l active volume)
Xenonbased TPC, read out in charge mode with an 8 cm $\times$8 cmsegmented
700 cm$^2$ plane (1152 ch) of the latest microbulkMicromegas technology. It
has been recently commissioned at University of Zaragoza as part of the R&D of
the NEXT $0\nu\beta\beta$ experiment, although the experiment's first stage is
currently being built based on a SiPM/PMTreadout concept relying on
electroluminescence. Around 2 million events were collected during the last
months, stemming from the low energy $\gamma$rays emitted by a $^{241}$Am
source when interacting with the Xenon gas ($\epsilon$ = 26, 30, 59.5 keV). The
localized nature of such events above atmospheric pressure, the long drift
times, as well as the possibility to determine their production time from the
associated $\alpha$ particle in coincidence, allow the extraction of primordial
properties of the TPC filling gas, namely the drift velocity, diffusion and
attachment coefficients. In this work we focus on the little explored
combination of Xe and trimethylamine (TMA) for which, in particular, such
properties are largely unknown. This gas mixture offers potential advantages
over pure Xenon when aimed at Rare Event Searches, mainly due to its Penning
characteristics, wavelength shifting properties and reduced diffusion, and it
is being actively investigated by our collaboration. The chamber is currently
operated at 2.7 bar, as an intermediate step towards the envisaged 10 bar. We
report here its performance as well as a first implementation of the
calibration procedures that have allowed the extension of the previously
reported energy resolution to the whole readout plane (10.6%FWHM@30keV).

Methods and an algorithm for computing the generalized Marcum $Q$function
($Q_{\mu}(x,y)$) and the complementary function ($P_{\mu}(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 noncentral chisquare 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 threeterm homogeneous
recurrence relations. A relative accuracy close to $10^{12}$ can be obtained
in the parameter region $(x,y,\mu) \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$).

NEXT is an experiment dedicated to neutrinoless double beta decay searches in
xenon. The detector is a TPC, holding 100 kg of highpressure xenon enriched in
the $^{136}$Xe isotope. It is under construction in the Laboratorio
Subterr\'aneo de Canfranc in Spain, and it will begin operations in 2015. The
NEXT detector concept provides an energy resolution better than 1% FWHM and a
topological signal that can be used to reduce the background. Furthermore, the
NEXT technology can be extrapolated to a 1ton scale experiment.

NEXTDEMO is a highpressure xenon gas TPC which acts as a technological
testbed and demonstrator for the NEXT100 neutrinoless double beta decay
experiment. In its current configuration the apparatus fully implements the
NEXT100 design concept. This is an asymmetric TPC, with an energy plane made
of photomultipliers and a tracking plane made of silicon photomultipliers
(SiPM) coated with TPB. The detector in this new configuration has been used to
reconstruct the characteristic signature of electrons in dense gas.
Demonstrating the ability to identify the MIP and "blob" regions. Moreover, the
SiPM tracking plane allows for the definition of a large fiducial region in
which an excellent energy resolution of 1.82% FWHM at 511 keV has been measured
(a value which extrapolates to 0.83% at the xenon Qbetabeta).

Highpressure xenon gas is an attractive detection medium for a variety of
applications in fundamental and applied physics. In this paper we study the
ionization and scintillation detection properties of xenon gas at 10 bar
pressure. For this purpose, we use a source of alpha particles in the NEXTDEMO
time projection chamber, the large scale prototype of the NEXT100 neutrinoless
double beta decay experiment, in three different drift electric field
configurations. We measure the ionization electron drift velocity and
longitudinal diffusion, and compare our results to expectations based on
available electron scattering cross sections on pure xenon. In addition, two
types of measurements addressing the connection between the ionization and
scintillation yields are performed. On the one hand we observe, for the first
time in xenon gas, large eventbyevent correlated fluctuations between the
ionization and scintillation signals, similar to that already observed in
liquid xenon. On the other hand, we study the field dependence of the average
scintillation and ionization yields. Both types of measurements may shed light
on the mechanism of electronion recombination in xenon gas for highlyionizing
particles. Finally, by comparing the response of alpha particles and electrons
in NEXTDEMO, we find no evidence for quenching of the primary scintillation
light produced by alpha particles in the xenon gas.