QuantumFDTD - A computational framework for the relativistic Schrödinger equation

We extend the publicly available quantumfdtd code. It was originally intended for solving the time-independent three-dimensional Schrödinger equation via the finite-difference time-domain (FDTD) method and for extracting the ground, first, and second excited states. We (a) include the case of the relativistic Schrödinger equation and (b) add two optimized FFT-based kinetic energy terms for the non-relativistic case. All the three new kinetic terms are computed using Fast Fourier Transform (FFT).We release the resulting code as version 3 of quantumfdtd. Finally, the code now supports arbitrary external filebased potentials and the option to project out distinct parity eigenstates from the solutions. Our goal is quark models used for phenomenological descriptions of QCD bound states, described by the three-dimensional Schrödinger equation. However, we target any field where solving either the non-relativistic or the relativistic three-dimensional Schrödinger equation is required

QuantumFDTD - A computational framework for the relativistic Schrödinger equation

We extend the publicly available quantumfdtd code. It was originally intended for solving the time-independent three-dimensional Schrödinger equation via the finite-difference time-domain (FDTD) method and for extracting the ground, first, and second excited states. We (a) include the case of the relativistic Schrödinger equation and (b) add two optimized FFT-based kinetic energy terms for the non-relativistic case. All the three new kinetic terms are computed using Fast Fourier Transform (FFT).We release the resulting code as version 3 of quantumfdtd. Finally, the code now supports arbitrary external filebased potentials and the option to project out distinct parity eigenstates from the solutions. Our goal is quark models used for phenomenological descriptions of QCD bound states, described by the three-dimensional Schrödinger equation. However, we target any field where solving either the non-relativistic or the relativistic three-dimensional Schrödinger equation is required

Static Energy in (2+1+12+1+1)-Flavor Lattice QCD: Scale Setting and Charm Effects

We present results for the static energy in ($2+1+1$)-flavor QCD over a wide range of lattice spacings and several quark masses, including the physical quark mass, with ensembles of lattice-gauge-field configurations made available by the MILC Collaboration. We obtain results for the static energy out to distances of nearly $1$~fm, allowing us to perform a simultaneous determination of the scales $r_{1}$ and $r_{0}$ as well as the string tension, $\sigma$. For the smallest three lattice spacings we also determine the scale $r_{2}$. Our results for $r_{0}/r_{1}$ and $r_{0}\sqrt{\sigma}$ agree with published ($2+1$)-flavor results. However, our result for $r_{1}/r_{2}$ differs significantly from the value obtained in the ($2+1$)-flavor case, which is most likely due to the effect of the charm quark. We also report results for $r_{0}$, $r_{1}$, and $r_{2}$ in~fm, with the former two being slightly lower than published ($2+1$)-flavor results. We study in detail the effect of the charm quark on the static energy by comparing our results on the finest two lattices with the previously published ($2+1$)-flavor QCD results at similar lattice spacing. We find that for $r > 0.2$~fm our results on the static energy agree with the ($2+1$)-flavor result, implying the decoupling of the charm quark for these distances. For smaller distances, on the other hand, we find that the effect of the dynamical charm quark is noticeable. The lattice results agree well with the two-loop perturbative expression of the static energy incorporating finite charm mass effects. This is the first time that the decoupling of the charm quark is observed and quantitatively analyzed on lattice data of the static energy.Comment: 50 pages, 37 figur