Related papers: Resolving the spurious-state problem in Dirac equa…
To solve the Dirac equation with the finite difference method, one has to face up to the spurious-state problem due to the fermion doubling problem when using the conventional central difference formula to calculate the first-order…
A finite difference scheme for the numerical treatment of the (3+1)D Dirac equation is presented. Its staggered-grid intertwined discretization treats space and time coordinates on equal footing, thereby avoiding the notorious fermion…
A numerical scheme utilizing a grid which is staggered in both space and time is proposed for the numerical solution of the (2+1)D Dirac equation in presence of an external electromagnetic potential. It preserves the linear dispersion…
We study a multigrid method for solving large linear systems of equations with tensor product structure. Such systems are obtained from stochastic finite element discretization of stochastic partial differential equations such as the…
A new approach to finite basis sets for the Dirac equation is developed. It solves the problem of spurious states and, as a result, improves the convergence properties of basis set calculations. The efficiency of the method is demonstrated…
We present a simple method to solve spherical harmonics moment systems, such as the the time-dependent $P_N$ and $SP_N$ equations, of radiative transfer. The method, which works for arbitrary moment order $N$, makes use of the specific…
A simple and reliable finite difference approach is presented for solution of the Dirac equation eigenproblem for states confined in rotationally symmetric systems. The method sets the boundary condition for the spinor wave function…
Multigrid methods were invented for the solution of discretized partial differential equations in ordered systems. The slowness of traditional algorithms is overcome by updates on various length scales. In this article we discuss…
We develop an algebraic multigrid method for solving the non-Hermitian Wilson discretization of the 2-dimensional Dirac equation. The proposed approach uses a bootstrap setup algorithm based on a multigrid eigensolver. It computes test…
In this paper we propose a novel arbitrary high order accurate semi-implicit space-time discontinuous Galerkin method for the solution of the two dimensional incompressible Navier-Stokes equations on staggered unstructured triangular…
In this paper we propose a novel arbitrary high order accurate semi-implicit space-time DG method for the solution of the three-dimensional incompressible Navier-Stokes equations on staggered unstructured curved tetrahedral meshes. As…
Critical slowing down in Krylov methods for the Dirac operator presents a major obstacle to further advances in lattice field theory as it approaches the continuum solution. Here we formulate a multi-grid algorithm for the Kogut-Susskind…
An efficient method, preconditioned conjugate gradient method with a filtering function (PCG-F), is proposed for solving iteratively the Dirac equation in 3D lattice space for nuclear systems. The filtering function is adopted to avoid the…
A challenging difficulty in solving the radial Dirac eigenvalue problem numerically is the presence of spurious (unphysical) eigenvalues among the correct ones that are neither related to mathematical interpretations nor to physical…
The Dirac-Bergmann algorithm is a recipe for converting a theory with a singular Lagrangian into a constrained Hamiltonian system. Constrained Hamiltonian systems include gauge theories -- general relativity, electromagnetism, Yang Mills,…
This work focuses on developing high-order energy-stable schemes for wave-dominated problems in closed domains using staggered finite-difference summation-by-parts (SBP FD) operators. We extend the previously presented uniform staggered…
Staggered grid finite difference scheme is widely used for the first order elastic wave equation, which constitutes the basis for least-squares reverse time migration and full waveform inversion. It is of great importance to improve the…
A finite difference scheme is presented for the Dirac equation in (1+1)D. It can handle space- and time-dependent mass and potential terms and utilizes exact discrete transparent boundary conditions (DTBCs). Based on a space- and…
The gauge symmetries of a general dynamical system can be systematically obtained following either a Hamiltonean or a Lagrangean approach. In the former case, these symmetries are generated, according to Dirac's conjecture, by the first…
In this paper we develop a staggered discontinuous Galerkin method for the Stokes and Darcy-Forchheimer problems coupled with the \Red{Beavers-Joseph-Saffman} conditions. The method is defined by imposing staggered continuity for all the…