Abstract: The Pierre Auger Observatory is making significant contributions towards understanding the nature and origin of ultra-high energy cosmic rays. One of its main challenges is the monitoring of the atmosphere, both in terms of its state variables and its optical properties. The aim of this work is to analyse aerosol optical depth tau(a)(z) values measured from 2004 to 2012 at the observatory, which is located in a remote and relatively unstudied area of Pampa Amarilla, Argentina. The aerosol optical depth is in average quite low – annual mean tau(a)(3.5 km) similar to 0.04 – and shows a seasonal trend with a winter minimum – tau(a)(3.5 km) – 0.03 -, and a summer maximum – tau(a)(3.5 km) similar to 0.06 -, and an unexpected increase from August to September tau(a)(35 km) similar to 0.055. We computed backward trajectories for the years 2005 to 2012 to interpret the air mass origin. Winter nights with low aerosol concentrations show air masses originating from the Pacific Ocean. Average concentrations are affected by continental sources (wind-blown dust and urban pollution), whilst the peak observed in September and October could be linked to biomass burning in the northern part of Argentina or air pollution coming from surrounding urban areas.

Abstract: We provide three Fortran programs which evaluate the QCD analytic (holomorphic) couplings A(v)(Q(2)) for complex or real squared momenta Q(2). These couplings are holomorphic analogs of the powers a(Q(2))(v) of the underlying perturbative QCD (pQCD) coupling a(Q(2)) equivalent to alpha(s)(Q(2))/pi, in three analytic QCD models (anQCD): Fractional Analytic Perturbation Theory (FAPT), Two-delta analytic QCD (2 delta anQCD), and Massive Perturbation Theory (MPT). The index v can be noninteger. The provided programs do basically the same job as the Mathematica package anQCD.m published by us previously (Ayala and Cvetic, 2015), but are now written in Fortran. Program summary Program title: AanQCDext Catalogue identifier: AEYKv10 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEYICv1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 12105 No. of bytes in distributed program, including test data, etc.: 98822 Distribution format: tar.gz Programming language: Fortran. Computer: Any work-station or PC where Fortran 95/200312008 (gfortran) is running. Operating system: Operating system Linux (Ubuntu and Scientific Linux), Windows (in all cases using gfortran). Classification: 11.1, 11.5. Nature of problem: Calculation of values of the running analytic couplings A(v)(Q(2); N-f) for general complex squared momenta Q(2) equivalent to -q(2), in three analytic QCD models, where A(v)(Q(2); N-f) is the analytic (holomorphic) analog of the power (alpha(s)(Q(2); N-f)/pi)(v). Here, A(v)(Q(2); N-f) is a holomorphic function in the Q(2) complex plane, with the exception of the negative semiaxis (-infinity, -M-thr(2)], reflecting the analyticity properties of the spacelike renormalization invariant quantities D(Q(2)) in QCD. In contrast, the perturbative QCD power (alpha(s)(Q(2); N-f)/pi)(v) has singularities even outside the negative semiaxis (Landau ghosts). The three considered models are: Analytic Perturbation theory (APT); Two-delta analytic QCD (2 delta anQCD); Massive Perturbation Theory (MPT). We refer to Ref. [1] for more details and literature. Solution method: The Fortran programs for FAPT and 2 delta anQCD models contain routines and functions needed to perform two-dimensional numerical integrations involving the spectral function, in order to evaluate A(v)(Q(2)) couplings. In MPT model, one-dimensional numerical integration involving A(1)(Q(2)) is sufficient to evaluate any A(v)(Q(2)) coupling. Restrictions: For unphysical choices of the input parameters the results are meaningless. When Q(2) is close to the cut region of the couplings (Q(2) real negative), the calculations can take more time and can have less precision. Running time: For evaluation of a set of about 10 related couplings, the times vary in the range t similar to 10(1)-10(2) s. MPT requires less time, t similar to 1-10(1) s. References: [1] C. Ayala and G. Cvetic, anQCD: a Mathematica package for calculations in general analytic QCD models, Comput. Phys. Commun. 190 (2015) 182.

Abstract: We present a pedagogical overview of the nonperturbative mechanism that endows gluons with a dynamical mass. This analysis is performed based on pure Yang-Mills theories in the Landau gauge, within the theoretical framework that emerges from the combination of the pinch technique with the background field method. In particular, we concentrate on the Schwinger-Dyson equation satisfied by the gluon propagator and examine the necessary conditions for obtaining finite solutions within the infrared region. The role of seagull diagrams receives particular attention, as do the identities that enforce the cancellation of all potential quadratic divergences. We stress the necessity of introducing nonperturbative massless poles in the fully dressed vertices of the theory in order to trigger the Schwinger mechanism, and explain in detail the instrumental role of these poles in maintaining the Becchi-Rouet-Stora-Tyutin symmetry at every step of the mass-generating procedure. The dynamical equation governing the evolution of the gluon mass is derived, and its solutions are determined numerically following implementation of a set of simplifying assumptions. The obtained mass function is positive definite, and exhibits a power law running that is consistent with general arguments based on the operator product expansion in the ultraviolet region. A possible connection between confinement and the presence of an inflection point in the gluon propagator is briefly discussed.