Solar Neutrino Masses and Mixing from Bilinear R-Parity Broken Supersymmetry: Analytical versus Numerical Results

We give an analytical calculation of solar neutrino masses and mixing at one-loop order within bilinear R-parity breaking supersymmetry, and compare our results to the exact numerical calculation. Our method is based on a systematic perturbative expansion of R-parity violating vertices to leading order. We find in general quite good agreement between approximate and full numerical calculation, but the approximate expressions are much simpler to implement. Our formalism works especially well for the case of the large mixing angle MSW solution (LMA-MSW), now strongly favoured by the recent KamLAND reactor neutrino data.Comment: 34 pages, 14 ps figs, some clarifying comments adde

Numerical solutions of the one-dimensional nucleon-meson cascade equations

Numerical integration of meson-nucleon cascade equations for accelerator shielding calculation

Numerical Calculation of Bessel Functions

A new computational procedure is offered to provide simple, accurate and flexible methods for using modern computers to give numerical evaluations of the various Bessel functions. The Trapezoidal Rule, applied to suitable integral representations, may become the method of choice for evaluation of the many Special Functions of mathematical physics.Comment: 10 page

Numerical calculation of Bessel, Hankel and Airy functions

The numerical evaluation of an individual Bessel or Hankel function of large order and large argument is a notoriously problematic issue in physics. Recurrence relations are inefficient when an individual function of high order and argument is to be evaluated. The coefficients in the well-known uniform asymptotic expansions have a complex mathematical structure which involves Airy functions. For Bessel and Hankel functions, we present an adapted algorithm which relies on a combination of three methods: (i) numerical evaluation of Debye polynomials, (ii) calculation of Airy functions with special emphasis on their Stokes lines, and (iii) resummation of the entire uniform asymptotic expansion of the Bessel and Hankel functions by nonlinear sequence transformations. In general, for an evaluation of a special function, we advocate the use of nonlinear sequence transformations in order to bridge the gap between the asymptotic expansion for large argument and the Taylor expansion for small argument ("principle of asymptotic overlap"). This general principle needs to be strongly adapted to the current case, taking into account the complex phase of the argument. Combining the indicated techniques, we observe that it possible to extend the range of applicability of existing algorithms. Numerical examples and reference values are given.Comment: 18 pages; 7 figures; RevTe

Associated production of charged Higgs bosons and top quarks with POWHEG

The associated production of charged Higgs bosons and top quarks at hadron colliders is an important discovery channel to establish the existence of a non-minimal Higgs sector. Here, we present details of a next-to-leading order (NLO) calculation of this process using the Catani-Seymour dipole formalism and describe its implementation in POWHEG, which allows to match NLO calculations to parton showers. Numerical predictions are presented using the PYTHIA parton shower and are compared to those obtained previously at fixed order, to a leading order calculation matched to the PYTHIA parton shower, and to a different NLO calculation matched to the HERWIG parton shower with MC@NLO. We also present numerical predictions and theoretical uncertainties for various Two Higgs Doublet Models at the Tevatron and LHC.Comment: 36 page

Loop-after-loop contribution to the second-order Lamb shift in hydrogenlike low-Z atoms

We present a numerical evaluation of the loop-after-loop contribution to the second-order self-energy for the ground state of hydrogenlike atoms with low nuclear charge numbers Z. The calculation is carried out in the Fried-Yennie gauge and without an expansion in Z \alpha. Our calculation confirms the results of Mallampalli and Sapirstein and disagrees with the calculation by Goidenko and coworkers. A discrepancy between different calculations is investigated. An accurate fitting of the numerical results provides a detailed comparison with analytic calculations based on an expansion in the parameter Z \alpha. We confirm the analytic results of order \alpha^2 (Z\alpha)^5 but disagree with Karshenboim's calculation of the \alpha^2 (Z \alpha)^6 \ln^3(Z \alpha)^{-2} contribution.Comment: RevTex, 19 pages, 4 figure

Numerical calculation of transonic boattail flow

A viscid-inviscid interaction procedure for the calculation of subsonic and transonic flow over a boattail was developed. This method couples a finite-difference inviscid analysis with an integral boundary-layer technique. Results indicate that the effect of the boundary layer is as important as an accurate inviscid method for this type of flow. Theoretical results from the solution of the full transonic-potential equation, including boundary layer effects, agree well with the experimental pressure distribution for a boattail. Use of the small disturbance transonic potential equation yielded results that did not agree well with the experimental results even when boundary-layer effects were included in the calculations

Transverse distinguishability of entangled photons with arbitrarily shaped spatial near- and far-field distributions

Entangled photons generated by spontaneous parametric down conversion inside a nonlinear crystal exhibit a complex spatial photon count distribution. A quantitative description of this distribution helps with the interpretation of experiments that depend on this structure. We developed a theoretical model and an accompanying numerical calculation that includes the effects of phase matching and the crystal properties to describe a wide range of spatial effects in two-photon experiments. The numerical calculation was tested against selected analytical approximations. We furthermore performed a double-slit experiment where we measured the visibility V and the distinguishability D and obtained $D^2 + V^2 = 1.43$. The numerical model accurately predicts these experimental results.Comment: 10 pages, 10 figure