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 propose a realization of the so-called bimodal/schizophrenic model proposed recently. We assume 54, the permutation group of four objects as flavor symmetry giving tri-bimaximal lepton mixing at leading order. In these models the second massive neutrino state is assumed quasi-Dirac and the remaining neutrinos are Majorana states. In the case of inverse mass hierarchy, the lower bound on the neutrinoless double beta decay parameter m(ee) is about two times that of the usual lower bound, within the range of sensitivity of the next generation of experiments.

Abstract: The framed standard model (FSM) is obtained from the standard model by incorporating, as field variables, the frame vectors (vielbeins) in internal symmetry space. It gives the standard Higgs boson and 3 generations of quarks and leptons as immediate consequences. It gives moreover a fermion mass matrix of the form: m = mT alpha alpha dagger, where alpha is a vector in generation space independent of the fermion species and rotating with changing scale, which has already been shown to lead, generically, to up-down mixing, neutrino oscillations and mass hierarchy. In this paper, pushing the FSM further, one first derives to 1-loop order the RGE for the rotation of alpha, and then applies it to fit mass and mixing data as a first test of the model. With 7 real adjustable parameters, 18 measured quantities are fitted, most (12) to within experimental error or to better than 0.5 percent, and the rest (6) not far off. (A summary of this fit can be found in Table 2 of this paper.) Two notable features, both generic to FSM, not just specific to the fit, are: (i) that a theta-angle of order unity in the instanton term in QCD would translate via rotation into a Kobayashi-Maskawa phase in the CKM matrix of about the observed magnitude (J similar to 10(-5)), (ii) that it would come out correctly that m(u) < m(d), despite the fact that m(t) >> m(b), m(c) >> m(s). Of the 18 quantities fitted, 12 are deemed independent in the usual formulation of the standard model. In fact, the fit gives a total of 17 independent parameters of the standard model, but 5 of these have not been measured by experiment.

Abstract: It is shown that if, from the starting point of a universal rank-one mass matrix long favored by phenomenologists, one adds the assumption that it rotates (changes its orientation in generation space) with changing scale, one can reproduce, in terms of only six real parameters, all the 16 mass ratios and mixing parameters of quarks and leptons. Of these 16 quantities so reproduced, 10 for which data exist for direct comparison (i.e. the CKM elements including the CP-violating phase, the angles theta(12), theta(13), theta(23) in nu-oscillation, and the masses m(c), m(mu), m(e)) agree well with experiment, mostly to within experimental errors; four others (m(s), m(u), m(d), m(nu 2)), the experimental values for which can only be inferred, agree reasonably well; while two others (m(nu 1), delta(CP) for leptons), not yet measured experimentally, remain as predictions. In addition, one gets as bonuses, estimates for (i) the right-handed neutrino mass m(nu R) and (ii) the strong CP angle theta inherent in QCD. One notes in particular that the output value for sin(2) 2 theta(13) from the fit agrees very well with recent experiments. By inputting the current experimental value with its error, one obtains further from the fit two new testable constraints: (i) that theta(23) must depart from its “maximal” value: sin(2) 2 theta(23) similar to 0.935 +/- 0.021, (ii) that the CP-violating (Dirac) phase in the PMNS would be smaller than in the CKM matrix: of order only vertical bar sin delta(CP)vertical bar <= 0.31 if not vanishing altogether.