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 have explored the question of whether the renormalization group running of the neutrino mixing parameters in the Minimal Supersymmetric Standard Model is detectable with ultra-high energy neutrinos from active galactic nuclei (AGN). We use as observables the ratios of neutrino fluxes produced at the AGN, focusing on four different neutrino production models: (Phi(0)(v epsilon+(v) over bar epsilon) : Phi(0)(v mu+(v) over bar mu) : Phi(0)(v tau+(v) over bar tau)) = (1 : 2 : 0), (0 : 1 : 0), (1 : 0 : 0), and (1 : 1 : 0). The prospects for observing deviations experimentally are taken into consideration, and we find out that it is necessary to impose a cut-off on the transferred momentum of Q(2) >= 10(7) GeV(2). However, this condition, together with the expected low value of the diffuse AGN neutrino flux, yields a negligible event rate at a km-scale. Cerenkov detector such as IceCube.

Abstract: During single clock inflation, hidden fields (i.e. fields coupled to the inflaton only gravitationally) in their adiabatic vacua can ordinarily only affect observables through virtual effects. After renormalizing background quantities (fixed by observations at some pivot scale), all that remains are logarithmic runnings in correlation functions that are both Planck and slow roll suppressed. In this paper we show how a large number of hidden fields can partially compensate this suppression and generate a potentially observable running in the tensor two point function, consistently inferable courtesy of a large N resummation. We detour to address certain subtleties regarding loop corrections during inflation, extending the analysis of [1]. Our main result is that one can extract bounds on the hidden field content of the universe from bounds on violations of the consistency relation between the tensor spectral index and the tensor to scalar ratio, were primordial tensors ever detected. Such bounds are more competitive than the naive bound inferred from requiring inflation to occur below the strong coupling scale of gravity if deviations from the consistency relation can be bounded to within the sub-percent level. We discuss how one can meaningfully constrain the parameter space of various phenomenological scenarios and constructions that address naturalness with a large number of species (such as N-naturalness') with CMB observations up to cosmic variance limits, and possibly future 21cm and gravitational wave observations.

Abstract: The nature of the three narrow hidden-charm pentaquark P-c states, i.e., P-c (4312), P-c (4440) and P-c (4457), is under intense discussion since their discovery from the updated analysis of the process Lambda(0)(b) -> I ) J/psi pK(-) by LHCb. In this work we extend our previous coupled-channel approach [Phys. Rev. Lett. 124, 072001 (2020)], in which the Pc states are treated as Sigma(()(c)*()) (D) over bar (()*()) molecules, by including the Lambda(c)(D) over bar (()*()) and eta(c)p as explicit inelastic channels in addition to the J/psi p, as required by unitarity and heavy quark spin symmetry (HQSS), respectively. Since inelastic parameters are very badly constrained by the current data, three calculation schemes are considered: (a) scheme I with pure contact interactions between the elastic, i.e., Sigma(()(c)*()) (D) over bar (()*()), and inelastic channels and without the Lambda(c)(D) over bar (()*()) interactions, (b) scheme II, where the one-pion exchange (OPE) is added to scheme I, and (c) scheme III, where the Lambda(c)(D) over bar (()*()) interactions are included in addition. It is shown that to obtain cutoff independent results, OPE in the multichannel system is to be supplemented with S-wave-to-D-wave mixing contact terms. As a result, in line with our previous analysis, we demonstrate that the experimental data for the J/psi p invariant mass distribution are consistent with the interpretation of the P-c(4312) and P-c(4440)/P-c(4457) as Sigma(c)(D) over bar and Sigma(c)(D) over bar* hadronic molecules, respectively, and that the data show clear evidence for a new narrow state, P-c(4380), identified as a Sigma(c)*(D) over bar molecule, which should exist as a consequence of HQSS. While two statistically equally good solutions are found in scheme I, only one of these solutions with the quantum numbers of the P-c (4440) and P-c (4457) being J(P) = 3/2(-) and 1/2(-), respectively, survives the requirement of regulator independence once the OPE is included. Moreover, we predict the line shapes in the elastic and inelastic channels and demonstrate that those related to the P-c (4440) and the P-c (4457) in the Sigma(()(c)*())<(D)over ( )anf eta(c)p mass distributions from Lambda(0)(b) ->( )Sigma(()(c)*()) (D) over barK(-) and Lambda(0)(b) -> eta(c)pK(-) will shed light on the quantum numbers of those states, once the data are available. We also investigate possible pentaquark signals in the Lambda(c)(D) over bar (()*()) final states.

Abstract: We extend the adiabatic regularization method by introducing an arbitrary mass scale μin the construction of the subtraction terms. This allows us to obtain, in a very robust way, the running of the coupling constants by demanding mu-invariance of the effective semiclassical (Maxwell-Einstein) equations. In particular, we get the running of the electric charge of perturbative quantum electrodynamics. Furthermore, the method brings about a renormalization of the cosmological constant and the Newtonian gravitational constant. The running obtained for these dimensionful coupling constants has new relevant (non-logarithmic) contributions, not predicted by dimensional regularization.