Related papers: Solving the Dirac equation on a GPU for strong-fie…
Immersed boundary methods (IBMs) facilitate the simulation of flows around stationary, moving, and deforming bodies on Cartesian grids. However, extending these simulations to the large grid sizes required for realistic flow problems…
The Dirac equation is an important model in relativistic quantum mechanics. In the semi-classical regime $\epsilon\ll1$, even a spatially spectrally accurate time splitting method \cite{HuJi:05} requires the mesh size to be $O(\epsilon)$,…
Methods for solving hyperbolic systems typically depend on unknown ordering (e.g., Gauss-Seidel, or sweep/wavefront/marching methods) to achieve good convergence. For many discretisations, mesh types or decompositions these methods do not…
An efficient solver for the three dimensional free-space Poisson equation is presented. The underlying numerical method is based on finite Fourier series approximation. While the error of all involved approximations can be fully controlled,…
Einstein-Dirac equations for two spinor fields are considered. It is shown that in this case one can obtain self-consistent equations set for these gravitating spinors. The key idea to obtain Einstein-Dirac equations is to use special…
Graphics Processing Units (GPUs) are high performance co-processors originally intended to improve the use and quality of computer graphics applications. Once, researchers and practitioners noticed the potential of using GPU for general…
The Dirac equation is solved for triangular and hexagonal graphene quantum dots for different boundary conditions in the presence of a perpendicular magnetic field. We analyze the influence of the dot size and its geometry on their energy…
The Dirac equation is a cornerstone of modern particle physics, which integrates special relativity and quantum mechanics into a consistent framework, yielding the prediction of electron and its antiparticle counterpart, positron. The Dirac…
We study massive Dirac fields in the background of the near-horizon limit of the extremal Myers-Perry black hole in five dimensions. We consider the case in which both angular momenta have equal magnitude. The resulting Dirac equation can…
Efficient methods for generating samples of wave packet trajectories are needed to build machine learning models for quantum dynamics. However, simulating such data by direct integration of the time-dependent Schrodinger equation can be…
We present a massively parallel, 3D phase-field simulation framework for modeling ferro-electric materials based scalable logic devices. We self-consistently solve the time-dependent Ginzburg Landau (TDGL) equation for ferroelectric…
The numerical solution of the Kadanoff-Baym nonlinear integro-differential equations, which yields the non-equilibrium Green's functions (NEGFs) of quantum many-body systems, poses significant computational challenges due to its high…
We present a new open source software for the integration of the radial Dirac equation developed specifically with muonic atoms in mind. The software, called mudirac, is written in C++ and can be used to predict frequencies and…
Solving quantifier-free non-linear real arithmetic (NRA) problems is a computationally hard task. To tackle this problem, prior work proposed a promising approach based on gradient descent. In this work, we extend their ideas and combine…
In this article we study the existence of solutions for the Dirac systems \begin{equation}\label{e:0.1} \left\{ \begin{array}{c} Pu=\frac{\partial H}{\partial v}(x,u,v) \quad\hbox{on} \ M, Pv=\frac{\partial H}{\partial u}(x,u,v)…
In this paper, we study the following nonlinear Dirac equations \begin{align*} \begin{cases} -i\sum\limits_{k=1}^3\alpha_k\partial_k u+m\beta u=f(x,|u|)u+\omega u, \displaystyle \int_{\mathbb{R}^3} |u|^2dx=a^2, \end{cases} \end{align*}…
In this paper, we investigate GPU based parallel triangular solvers systematically. The parallel triangular solvers are fundamental to incomplete LU factorization family preconditioners and algebraic multigrid solvers. We develop a new…
We present an algorithm specifically tailored for solving kinetic equations onto GPUs. The efficiency of the algorithm is demonstrated by solving the one-dimensional shock wave structure problem and a two-dimensional low Mach number driven…
Solving discretized versions of the Dirac equation represents a large share of execution time in lattice Quantum Chromodynamics (QCD) simulations. Many high-performance computing (HPC) clusters use graphics processing units (GPUs) to offer…
We present our implementation of the RHMC algorithm for staggered fermions on Graphics Processing Units using the NVIDIA CUDA programming language. While previous studies exclusively deal with the Dirac matrix inversion problem, our code…