Related papers: High-precision minmax solution of the two-center D…
We present a precise fully relativistic numerical solution of the two-center Coulomb problem. The special case of unit nuclear charges is relevant for the accurate description of the ${\rm H}_2^+$ molecular ion and its isotopologues,…
The Dirac equation for H$_2^+$ is solved numerically by expansion in a basis set of two-center exponential functions, using different kinetic balance schemes. Very high precision (27-32 digits) is achieved, either with the dual kinetic…
We introduce \texttt{featom}, an open source code that implements a high-order finite element solver for the radial Schr\"odinger, Dirac, and Kohn-Sham equations. The formulation accommodates various mesh types, such as uniform or…
The electronic spectra of light one-electron quasi-molecular compounds H-H$^+$, He$^+$-He$^2+$ and He$^+$-H$^+$ are analyzed. To this end, the two-center Dirac equation is solved by the dual-kinetically balanced finite-basis-set method for…
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 rigorous two-center approach based on the dual-kinetically balanced finite-basis-set expansion is applied to one-electron, heteronuclear diatomic Bi-Au, U-Pb, and Cf-U quasimolecules. The obtained $1\sigma$ ground-state energies are…
The solutions of elliptic problems with a Dirac measure in right-hand side are not H1 and therefore the convergence of the finite element solutions is suboptimal. Graded meshes are standard remedy to recover quasi-optimality, namely…
A finite element method using B-splines is presented and compared with a conventional finite element method of Lagrangian type. The efficiency of both methods has been investigated at the example of a coupled non-linear system of Dirac…
In this chapter we focus first on the theoretical methods and relevant computational approaches to calculate the electronic structure of atoms, molecules, and clusters containing heavy elements for which relativistic effects become…
Obtaining high-precision guaranteed lower eigenvalue bounds remains difficult, even though the standard high-order conforming finite element (FEM) easily yields extremely sharp upper bounds. Recently developed rigorous approaches using such…
The Dirac equation for H$_2^+$ is solved numerically using an iterative method proposed by Kutzelnigg [Z. Phys. 11, 15 (1989]. The four-component wavefunction is expanded in a newly introduced kinetically balanced exponential basis set. The…
Using the standard finite element method (FEM) to solve general partial differential equations, the round-off error is found to be proportional to $N^{\beta_{\rm R}}$, with $N$ the number of degrees of freedom (DoFs) and $\beta_{\rm R}$ a…
Particle breakage due to collisional interactions plays a vital role in the development of several phenomena in science and engineering. The nonlinear collisional breakage equations (NCBEs) are a significant set of equations in this…
The adaptive nonconforming Morley finite element method (FEM) approximates a regular solution to the von K\'{a}rm\'{a}n equations with optimal convergence rates for sufficiently fine triangulations and small bulk parameter in the D\"orfler…
We present and analyze a linearized finite element method (FEM) for the dynamical incompressible magnetohydrodynamics (MHD) equations. The finite element approximation is based on mixed conforming elements, where Taylor--Hood type elements…
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…
In this paper, we present a finite element method (FEM) framework enhanced by an operator-adapted wavelet decomposition algorithm designed for the efficient analysis of multiscale electromagnetic problems. Usual adaptive FEM approaches,…
Thanks to a finite element method, we solve numerically parabolic partial differential equations on complex domains by avoiding the mesh generation, using a regular background mesh, not fitting the domain and its real boundary exactly. Our…
The Morley finite element method (FEM) is attractive for semilinear problems with the biharmonic operator as a leading term in the stream function vorticity formulation of 2D Navier-Stokes problem and in the von K\'{a}rm\'{a}n equations.…
In this paper, we provide a procedure to solve the eigen solutions of Dirac equation with complicated potential approximately. At first, we solve the eigen solutions of a linear Dirac equation with complete eigen system, which approximately…