Recent advances in higher order quasi-Monte Carlo methods

In this article we review some of recent results on higher order quasi-Monte Carlo (HoQMC) methods. After a seminal work by Dick (2007, 2008) who originally introduced the concept of HoQMC, there have been significant theoretical progresses on HoQMC in terms of discrepancy as well as multivariate numerical integration. Moreover, several successful and promising applications of HoQMC to partial differential equations with random coefficients and Bayesian estimation/inversion problems have been reported recently. In this article we start with standard quasi-Monte Carlo methods based on digital nets and sequences in the sense of Niederreiter, and then move onto their higher order version due to Dick. The Walsh analysis of smooth functions plays a crucial role in developing the theory of HoQMC, and the aim of this article is to provide a unified picture on how the Walsh analysis enables recent developments of HoQMC both for discrepancy and numerical integration

Hybrid discrete (H TN) approximations to the equation of radiative transfer

The linear kinetic transport equations are ubiquitous in many application areas, including as a model for neutron transport in nuclear reactors and the propagation of electromagnetic radiation in astrophysics. The main computational challenge in solving the linear transport equations is that solutions live in a high-dimensional phase space that must be sufficiently resolved for accurate simulations. The three standard computational techniques for solving the linear transport equations are the (1) implicit Monte Carlo, (2) discrete ordinate(S$_N$), and (3) spherical harmonic(P$_N$) methods. Monte Carlo methods are stochastic methods for solving time-dependent nonlinear radiative transfer problems. In a traditional Monte Carlo method when photons are absorbed, they are reemitted in a distribution which is uniform over the entire spatial cell where the temperature is assumed constant, resulting in loss of information. In implicit Monte Carlo(IMC) methods, photons are reemitted from the place where they were actually absorbed, which improves the accuracy. Overall, IMC method improves stability, flexibility, and computational efficiency \cite{fleck}. The S$_N$ method solves the transport equation using a quadrature rule to reconstruct the energy density. This method suffers from so-called ray effect , which are due to the approximation of the double integral over a unit sphere by a finite number of discrete angular directions \cite{chai}. The P$_N$ approximation is based on expanding the part of the solution that depends on velocity direction (i.e., two angular variables) into spherical harmonics. A big challenge with the P$_N$ approach is that the spherical harmonics expansion does not prevent the formation of negative particle concentrations. The idea behind my research is to develop on an alternative formulation of P$_N$ approximations that hybridizes aspects of both P$_N$ and S$_N$. Although the basic scheme does not guarantee positivity of the solution, the new formulation allows for the introduction of local limiters that can be used to enforce positivity

Mathematical Properties of a New Levin-Type Sequence Transformation Introduced by \v{C}\'{\i}\v{z}ek, Zamastil, and Sk\'{a}la. I. Algebraic Theory

\v{C}\'{\i}\v{z}ek, Zamastil, and Sk\'{a}la [J. Math. Phys. \textbf{44}, 962 - 968 (2003)] introduced in connection with the summation of the divergent perturbation expansion of the hydrogen atom in an external magnetic field a new sequence transformation which uses as input data not only the elements of a sequence $\{s_n \}_{n=0}^{\infty}$ of partial sums, but also explicit estimates $\{\omega_n \}_{n=0}^{\infty}$ for the truncation errors. The explicit incorporation of the information contained in the truncation error estimates makes this and related transformations potentially much more powerful than for instance Pad\'{e} approximants. Special cases of the new transformation are sequence transformations introduced by Levin [Int. J. Comput. Math. B \textbf{3}, 371 - 388 (1973)] and Weniger [Comput. Phys. Rep. \textbf{10}, 189 - 371 (1989), Sections 7 -9; Numer. Algor. \textbf{3}, 477 - 486 (1992)] and also a variant of Richardson extrapolation [Phil. Trans. Roy. Soc. London A \textbf{226}, 299 - 349 (1927)]. The algebraic theory of these transformations - explicit expressions, recurrence formulas, explicit expressions in the case of special remainder estimates, and asymptotic order estimates satisfied by rational approximants to power series - is formulated in terms of hitherto unknown mathematical properties of the new transformation introduced by \v{C}\'{\i}\v{z}ek, Zamastil, and Sk\'{a}la. This leads to a considerable formal simplification and unification.Comment: 41 + ii pages, LaTeX2e, 0 figures. Submitted to Journal of Mathematical Physic

A numerical adaptation of SAW identities from the honeycomb to other 2D lattices

Recently, Duminil-Copin and Smirnov proved a long-standing conjecture by Nienhuis that the connective constant of self-avoiding walks on the honeycomb lattice is $\sqrt{2+\sqrt{2}}.$ A key identity used in that proof depends on the existence of a parafermionic observable for self-avoiding walks on the honeycomb lattice. Despite the absence of a corresponding observable for SAW on the square and triangular lattices, we show that in the limit of large lattices, some of the consequences observed on the honeycomb lattice persist on other lattices. This permits the accurate estimation, though not an exact evaluation, of certain critical amplitudes, as well as critical points, for these lattices. For the honeycomb lattice an exact amplitude for loops is proved.Comment: 21 pages, 7 figures. Changes in v2: Improved numerical analysis, giving greater precision. Explanation of why we observe what we do. Extra reference

Quarkonium mass splittings in three-flavor lattice QCD

We report on calculations of the charmonium and bottomonium spectrum in lattice QCD. We use ensembles of gauge fields with three flavors of sea quarks, simulated with the asqtad improved action for staggered fermions. For the heavy quarks we employ the Fermilab interpretation of the clover action for Wilson fermions. These calculations provide a test of lattice QCD, including the theory of discretization errors for heavy quarks. We provide, therefore, a careful discussion of the results in light of the heavy-quark effective Lagrangian. By and large, we find that the computed results are in agreement with experiment, once parametric and discretization errors are taken into account.Comment: 21 pages, 17 figure

Calculation of Hydrogenic Bethe Logarithms for Rydberg States

We describe the calculation of hydrogenic (one-loop) Bethe logarithms for all states with principal quantum numbers n <= 200. While, in principle, the calculation of the Bethe logarithm is a rather easy computational problem involving only the nonrelativistic (Schroedinger) theory of the hydrogen atom, certain calculational difficulties affect highly excited states, and in particular states for which the principal quantum number is much larger than the orbital angular momentum quantum number. Two evaluation methods are contrasted. One of these is based on the calculation of the principal value of a specific integral over a virtual photon energy. The other method relies directly on the spectral representation of the Schroedinger-Coulomb propagator. Selected numerical results are presented. The full set of values is available at quant-ph/0504002.Comment: 10 pages, RevTe

Mass Spectra and Decay Constants of Heavy-Light Mesons: A Case Study of QCD Sum Rules and Quark Model

In this paper we visited mass spectra and decay constants of pseudoscalar and vector heavy-light mesons ($B$, $B_s$, $D$ and $D_S$) in the framework of QCD sum rule and quark model. The harmonic oscillator wave function was used in quark model while a simple interpolating current was used in QCD sum rule calculation. We obtained good results in accordance with the available experimental data and theoretical studies.Comment: 11 pages, 2 tables, 8 figure