Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNumerical Solution of the Time-Dependent Dirac Equation in Coordinate Space without Fermion-Doubling
103
0
0.0
(
0
)
Added by
Francois Fillion-Gourdeau
Publication date
2011
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
The validation and parallel implementation of a numerical method for the solution of the time-dependent Dirac equation is presented. This numerical method is based on a split operator scheme where the space-time dependence is computed in coordinate space using the method of characteristics. Thus, most of the steps in the splitting are calculated exactly, making for a very efficient and unconditionally stable method. We show that it is free from spurious solutions related to the fermion-doubling problem and that it can be parallelized very efficiently. We consider a few simple physical systems such as the time evolution of Gaussian wave packets and the Klein paradox. The numerical results obtained are compared to analytical formulas for the validation of the method.
rate research
Read More
We develop an approach to solving numerically the time-dependent Schrodinger equation when it includes source terms and time-dependent potentials. The approach is based on the generalized Crank-Nicolson method supplemented with an Euler-MacLaurin expansion for the time-integrated nonhomogeneous term. By comparing the numerical results with exact solutions of analytically solvable models, we find that the method leads to precision comparable to that of the generalized Crank-Nicolson method applied to homogeneous equations. Furthermore, the systematic increase in precision generally permits making estimates of the error.
Two numerical methods are used to evaluate the relativistic spectrum of the two-centre Coulomb problem (for the $H_{2}^{+}$ and $Th_{2}^{179+}$ diatomic molecules) in the fixed nuclei approximation by solving the single particle time-independent Dirac equation. The first one is based on a min-max principle and uses a two-spinor formulation as a starting point. The second one is the Rayleigh-Ritz variational method combined with kinematically balanced basis functions. Both methods use a B-spline basis function expansion. We show that accurate results can be obtained with both methods and that no spurious states appear in the discretization process.
A theoretical study of the intense-field multiphoton ionization of hydrogenlike systems is performed by solving the time-dependent Dirac equation within the dipole approximation. It is shown that the velocity-gauge results agree to the ones in length gauge only if the negative-energy states are included in the time propagation. On the other hand, for the considered laser parameters, no significant difference is found in length gauge if the negative-energy states are included or not. Within the adopted dipole approximation the main relativistic effect is the shift of the ionization potential. A simple scaling procedure is proposed to account for this effect.
We propose an improvement of the basis for the solution of the stationary two-centre Dirac equation in Cassini coordinates using the finite-basis-set method presented in Ref. [1]. For the calculations in Ref. [1], we constructed the basis for approximating the energy eigenfunctions by using smooth piecewise defined polynomials, called B-splines. In the present article, we report that an analysis of the employed representation of the Dirac matrices shows that the above approximation is not efficient using B-splines only. Therefore, we include basis functions which are defined using functions with step-like behaviour instead of B-splines. Thereby, we achieve a significant increase of accuracy of results as compared to Ref. [1].
We examine the performance of various time propagation schemes using a one-dimensional model of the hydrogen atom. In this model the exact Coulomb potential is replaced by a soft-core interaction. The model has been shown to be a reasonable representation of what occurs in the fully three-dimensional hydrogen atom. Our results show that while many numerically simple (low order) propagation schemes work, they often require quite small time-steps. Comparing them against more accurate methods, which may require more work per time-step but allow much larger time-steps, can be illuminating. We show that at least in this problem, the compute time for a number of the more accurate methods is actually less than lower order schemes. Finally, we make some remarks on what to expect in generalizing our findings to more than one dimension.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
