To accelerate the solution of large eigenvalue problems arising from many-body calculations in nuclear physics on distributed-memory parallel systems equipped with general-purpose Graphic Processing Units (GPUs), we modified a previously developed hybrid MPI/OpenMP implementation of an eigensolver written in FORTRAN 90 by using an OpenACC directives based programming model. Such an approach requires making minimal changes to the original code and enables a smooth migration of large-scale nuclear structure simulations from a distributed-memory many-core CPU system to a distributed GPU system. However, in order to make the OpenACC based eigensolver run efficiently on GPUs, we need to take into account the architectural differences between a many-core CPU and a GPU device. Consequently, the optimal way to insert OpenACC directives may be different from the original way of inserting OpenMP directives. We point out these differences in the implementation of sparse matrix-matrix multiplications (SpMM), which constitutes the main cost of the eigensolver, as well as other differences in the preconditioning step and dense linear algebra operations. We compare the performance of the OpenACC based implementation executed on multiple GPUs with the performance on distributed-memory many-core CPUs, and demonstrate significant speedup achieved on GPUs compared to the on-node performance of a many-core CPU. We also show that the overall performance improvement of the eigensolver on multiple GPUs is more modest due to the communication overhead among different MPI ranks.
We study the semileptonic decay of the $B_c (0^-)$ to charmonia through the bottom-to-charm-quark electroweak current in the framework of basis light-front quantization. Explicitly, we calculate the weak transition form factors for processes of $B_c$ decaying into $eta_c$ or $J/psi$ based on the corresponding initial and final valence light-front wave functions both obtained from the basis light-front quantization. We also present the corresponding differential decay width and branching ratios, as well as the branching ratios for decays into the $eta_c$ and the $psi$. We observe unphysical frame dependence of the calculated form factors, which is attributed to only including the valence light-front wave functions. Based on the analysis of current component, we propose a preferred set of frames to calculate these form factors.
We investigate the scattering of a quark on a heavy nucleus at high energies using the time-dependent basis light-front quantization (tBLFQ) formalism, which is the first application of the tBLFQ formalism in QCD. We present the real-time evolution of the quark wave function in a strong classical color field of the relativistic nucleus, described as the Color Glass Condensate. The quark and the nucleus color field are simulated in the QCD SU(3) color space. We calculate the total and the differential cross sections, and the quark distribution in coordinate and color spaces using the tBLFQ approach. We recover the eikonal cross sections in the eikonal limit. We find that the differential cross section from the tBLFQ simulation is in agreement with a perturbative calculation at large $p_perp$, and it deviates from the perturbative calculation at small $p_perp$ due to higher-order contributions. In particular, we relax the eikonal limit by letting the quark carry realistic finite longitudinal momenta. We study the sub-eikonal effect on the quark through the transverse coordinate distribution of the quark with different longitudinal momentum, and we find the sub-eikonal effect to be sizable. Our results can significantly reduce the theoretical uncertainties in small $p_perp$ region which has important implications to the phenomenology of the hadron-nucleus and deep inelastic scattering at high energies.
The light-front wave functions of hadrons allow us to calculate a wide range of physical observables; however, the wave functions themselves cannot be measured. We discuss recent results for quarkonia obtained in basis light-front quantization using an effective Hamiltonian with a confining model in both the transverse and longitudinal directions and with explicit one-gluon exchange. In particular, we focus on the numerical convergence of the basis expansion, as well as the asymptotic behavior of the light-front wave functions. We also illustrate that, for mesons with unequal quark masses, the maxima of the light-front wave functions depend in a non-trivial way on the valence quark-mass difference.
We solve the Minkowski-space Schwinger-Dyson equation (SDE) for the fermion propagator in quantum electrodynamics (QED) with massive photons. Specifically, we work in the quenched approximation within the rainbow-ladder truncation. Loop-divergences are regularized by the Pauli-Villars regularization. With moderately strong fermion-photon coupling, we find that the analytic structure of the fermion propagator consists of an on-shell pole and branch-cuts located in the timelike region. Such structures are consistent with the direct solution of the fermion propagator as functions of the complex momentum. Our method paves the way towards the calculation of the Minkowski-space Bethe-Salpeter amplitude using dressed fermion propagator.
We calculate the transition form factor between vector and pseudoscalar quarkonia in both the timelike and the spacelike region using light-front dynamics. We investigate the frame dependence of the form factors for heavy quarkonia with light-front wavefunctions calculated from the valence Fock sector. This dependence could serve as a measure for the Lorentz symmetry violation arising from the Fock-space truncation. We suggest using frames with minimal longitudinal momentum transfer for calculations in the valence Fock sector, namely the Drell-Yan frame for the space-like region and a specific longitudinal frame for the timelike region; at $q^2=0$ these two frames give the same result. We also use the transition form factor in the timelike region to investigate the electromagnetic Dalitz decay $psi_Ato psi_B l^+l^-$ ($l = e,mu$) and predict the effective mass spectrum of the lepton pair.
We study exclusive charmonium production in diffractive deep inelastic scattering and ultra-peripheral heavy-ion collisions within the dipole picture. The mass spectrum and light-front wavefunctions of charmonium are obtained from the basis light-front quantization approach, using the one-gluon exchange interaction plus a confining potential inspired by light-front holography. We apply these light-front wavefunctions to exclusive charmonium production. The resulting cross sections are in reasonable agreement with electron-proton collision data at HERA and ultra-peripheral nucleus collision measurements at RHIC and LHC. The charmonium cross-section has model dependence on the dipole model. We observe that the cross-section ratio of excited states to the ground state has a weaker dependence than the cross-section itself. We suggest that measurements of excited states of heavy quarkonium production in future electron-ion collision experiments will impose rigorous constraints on heavy quarkonium light-front wavefunctions, thus improving our understanding of meson structure, which eventually will help us develop a precise description of the gluon distribution function in the small-$x$ regime.
We solve for properties of $^6$Li in the ab initio No-Core Full Configuration approach and we separately solve for its ground state and $J^{pi}=2_{2}^{+}$ resonance with the Gamow Shell Model in the Berggren basis. We employ both the JISP16 and chiral NNLO$_{opt}$ realistic nucleon-nucleon interactions and investigate the ground state energy, excitation energies, point proton root-mean-square radius and a suite of electroweak observables. We also extend and test methods to extrapolate the ground state energy, point proton root-mean-square radius, and electric quadrupole moment. We attain improved estimates of these observables in the No-Core Full Configuration approach by using basis spaces up through N$_{max}$=18 that enable more definitive comparisons with experiment. Using the Density Matrix Renormalization Group approach with the JISP16 interaction, we find that we can significantly improve the convergence of the Gamow Shell Model treatment of the $^6$Li ground state and $J^{pi}=2_{2}^{+}$ resonance by adopting a natural orbital single-particle basis.
Basis Light-front Quantization has been proposed as a nonperturbative framework for solving quantum field theory. We apply this approach to Quantum Electrodynamics and explicitly solve for the light-front wave function of a physical electron. Based on the resulting light-front wave function, we evaluate the electron anomalous magnetic moment. Nonperturbative mass renormalization is performed. Upon extrapolation to the infinite basis limit our numerical results agree with the Schwinger result obtained in perturbation theory to an accuracy of 0.06%.
We introduce a nonperturbative, first-principles approach to time-dependent problems in quantum field theory. In this approach, the time-evolution of quantum field configurations is calculated in real time and at the amplitude level. This method is particularly suitable for treating systems interacting with a time-dependent background field. As a test problem, we apply this approach to QED and study electron acceleration and the associated photon emission in a time- and space-dependent electromagnetic background field.
