Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DIRAC EQUATION”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Applications of Two Body Dirac Equations to Hadron and Positronium Spectroscopy

We review recent applications of the Two Body Dirac equations of constraint dynamics to meson spectroscopy and describe new extensions to three-body problems in their use in the study of baryon spectroscopy. We outline unique aspects of these equations for QED bound states that distinguish them among the various other approaches to the relativistic two body problem. Finally we discuss recent theoretical solutions of new peculiar bound states for positronium arising from the Two Body Dirac equations of constraint dynamics, assuming point particles for the electron and the positron.

Crater, H. W.↗
DOE OSTI2017

Approximate analytic solutions to coupled nonlinear Dirac equations

Here, we consider the coupled nonlinear Dirac equations (NLDEs) in 1+11+1 dimensions with scalar–scalar self-interactions g 1 2 /2($\bar{ψ}$ψ) 2 + g 2 2 /2($\bar{Φ}$Φ) 2 + g 2 3 ($\bar{ψ}$ψ)($\bar{Φ}$Φ) as well as vector–vector interactions g 1 2 /2($\bar{ψ}$γμψ)($\bar{ψ}$γμψ) + g 2 2 /2($\bar{Φ}$γμΦ)($\bar{Φ}$γμΦ) + g 2 3 ($\bar{ψ}$γμψ)($\bar{Φ}$γμΦ). Writing the two components of the assumed rest frame solution of the coupled NLDE equations in the form ψ=e –iω1t R 1 cosθ,R 1 sinθΦ=e –iω2t R 2 cosη,R 2 sinη, and assuming that θ(x),η(x) have the same functional form they had when g3 = 0, which is an approximation consistent with the conservation laws, we then find approximate analytic solutions for Ri(x) which are valid for small values of g 3 2 /g 2 2 and g 3 2 /g 1 2 . In the nonrelativistic limit we show that both of these coupled models go over to the same coupled nonlinear Schrödinger equation for which we obtain two exact pulse solutions vanishing at x → ±∞.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Error estimates of finite difference methods for the Dirac equation in the massless and nonrelativistic regime

We present four frequently used finite difference methods and establish the error bounds for the discretization of the Dirac equation in the massless and nonrelativistic regime, involving a small dimensionless parameter 0 < ε &NestedLessLess; 1 inversely proportional to the speed of light. In the massless and nonrelativistic regime, the solution exhibits rapid motion in space and is highly oscillatory in time. Specifically, the wavelength of the propagating waves in time is at O(ε), while in space, it is at O(1) with the wave speed at O(ε -1 ). We adopt one leap-frog, two semi-implicit, and one conservative Crank-Nicolson finite difference methods to numerically discretize the Dirac equation in one dimension and establish rigorously the error estimates which depend explicitly on the time step τ, mesh size h, and the small parameter ε. The error bounds indicate that, to obtain the “correct” numerical solution in the massless and nonrelativistic regime, i.e., 0 < ε &NestedLessLess; 1, all these finite difference methods share the same ε-scalability as time step τ = O(ε 3/2 ) and mesh size h = O(ε 1/2 ). A large number of numerical results are reported to verify the error estimates.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Quantum structure emerging from the transformation design of the Dirac equation

It is shown that a quantum structure can be created by a set of chosen constraint conditions that emerge from the transformation design of the Dirac equation in general relativity. As an explanation, the constraints that cause novel bound states with the quantization rule of a 2D Coulomb system are presented. The discussion in this paper provides a systematic way to look for constraints that generate a required quantization rule. -- Highlights: •We perform the transformation design of space and time for spin-1/2 matter waves. •A quantum rule could naturally emerge as constraints imposed on the Dirac equation itself. •New fermion states share the quantum spectrum of a 2D Coulomb system. •Transformation design uncovers a new exact solvable model. •A quantum spectrum can be created by a geometric structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Phase-integral method for the radial Dirac equation

A phase-integral (WKB) solution of the radial Dirac equation is calculated up to the third order of approximation, retaining perfect symmetry between the two components of the wave function and introducing no singularities except at the zeroth-order transition points. The potential is allowed to be of scalar, vector, or tensor type, or any combination of these. The connection problem is investigated in detail. Explicit formulas are given for single-turning-point phase shifts and single-well energy levels.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Spatial resolution of different discretizations over long-time for the Dirac equation with small potentials

In this report we compare the long-time error bounds and spatial resolution of finite difference methods with different spatial discretizations for the Dirac equation with small electromagnetic potentials characterized by $\varepsilon \in (0, 1]$ a dimensionless parameter. We begin with the simple and widely used finite difference time domain (FDTD) methods, and establish rigorous error bounds of them, which are valid up to the time at $O(1/\varepsilon)$. In the error estimates, we pay particular attention to how the errors depend explicitly on the mesh size $h$ and time step $\tau$ as well as the small parameter $\varepsilon$. Based on the results, in order to obtain "correct" numerical solutions up to the time at $O(1/\varepsilon)$, the $\varepsilon$-scalability (or meshing strategy requirement) of the FDTD methods should be taken as $h = O(\varepsilon^{1/2})$ and $\tau = O(\varepsilon^{1/2})$. To improve the spatial resolution capacity, we apply the Fourier spectral method to discretize the Dirac equation in space. Error bounds of the resulting finite difference Fourier pseudospectral (FDFP) methods show that they exhibit uniform spatial errors in the long-time regime, which are optimal in space as suggested by the Shannon's sampling theorem. Extensive numerical results are reported to confirm the error bounds and demonstrate that they are sharp.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Analytical continuation from bound to resonant states in the Dirac equation with quadrupole-deformed potentials

Here, resonances with pronounced single-particle characteristics are crucial for quantitative descriptions of exotic nuclei near and beyond the drip lines, and often impact halo formation and nucleon decay processes. Since the majority of nuclei are deformed, the interplay between deformation and orbital structure near threshold can lead to improved descriptions of exotic nuclei. Develop a method to study single-particle resonant orbital structure in the Dirac equation with a quadrupole-deformed Woods-Saxon potential. Determine the structure evolution of bound and resonant levels with deformation in this scheme, and examine the impact on halo formation in loosely bound systems, with a focus on the recent halo candidate nucleus 37 Mg. Analytical continuation of the coupling constant (ACCC) method is developed on the basis of the Dirac equation with a deformed Woods-Saxon potential. The scalar and vector terms in the deformed potential are determined by the energies of the valence neutron and nearby orbitals, which are extracted from a self-consistent relativistic Hartree-Bogoliubov (RHB) calculation with the PC-PK1 density functional. As a result, we compare the energies and widths of resonant orbitals in the recent halo nucleus candidate 37 Mg using the ACCC method based on the Dirac coupled-channel equations with those determined from the scattering phase shift (SPS) method. It is found that the results from the two methods agree well for narrow resonances, whereas the SPS method fails for broad resonances. Nilsson levels for bound and resonant orbitals from the ACCC method are calculated over a wide range of deformations and show some decisive hints of halo formation in 37 Mg. In conclusion, in our ACCC model for deformed potentials in the coupled-channel Dirac equations, the crossing of the configuration 1/2[321] and 5/2[312] orbitals at a deformation of approximately 0.5 enhances the probability to occupy the 1/2[321] orbital coming from 2 p3/2 thereby explaining the recent observation of a p-wave one-neutron halo configuration in 37 Mg. The resonant 1/2[301] configuration plays a crucial role in halo formation in the magnesium isotopes beyond A = 40 for a wide range of deformations larger than 0.2.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Speed-of-light pulses in the massless nonlinear Dirac equation with a potential

In this work, we consider the massless nonlinear Dirac (NLD) equation in 1 + 1 dimension with scalar-scalar self-interaction g 2 2 ( Ψ ¯ Ψ ) 2 in the presence of three external electromagnetic real potentials V ( x ) , a potential barrier, a constant potential, and a potential well. By solving numerically the NLD equation, we find different scenarios depending on initial conditions, namely, propagation of the initial pulse along one direction, splitting of the initial pulse into two pulses traveling in opposite directions, and focusing of two initial pulses followed by a splitting. For all considered cases, the final waves travel with the speed of light and are solutions of the massless linear Dirac equation. During these processes the charge and the energy are conserved, whereas the momentum is conserved when the solutions possess specific symmetries. For the case of the constant potential, we derive exact analytical solutions of the massless NLD equation that are also solutions of the massless linearized Dirac equation. Decay or growth of the initial pulse is also predicted from the evolution of the charge for the case of a non-zero imaginary part of the potential.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Analytical solutions of the Dirac equation under Hellmann–Frost–Musulin potential

The approximate analytical solutions of the Dirac equation with Hellmann–Frost–Musulin potential have been studied by using the generalized parametric Nikiforov–Uvarov (NU) method for arbitrary spin–orbit quantum number k under the spin and pseudospin symmetries. The Hellmann–Frost–Musulin potential is a superposition potential that consists of Yukawa potential, Coulomb potential, and Frost–Musulin potential. As a particular case, we found the energy levels of the non-relativistic limit of the spin symmetry. The energy equation of Yukawa potential, Coulomb potential, Hellmann potential and Frost–Musulin potential are obtained. Energy values are generated for some diatomic molecules.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

General solution of the Dirac equation for quasi-two-dimensional electrons

The general solution of the Dirac equation for quasi-two-dimensional electrons confined in an asymmetric quantum well, is found. The energy spectrum of such a system is exactly calculated using special unitary operator and is shown to depend on the electron spin polarization. This solution contains free parameters, whose variation continuously transforms one known particular solution into another. As an example, two different cases are considered in detail: electron in a deep and in a strongly asymmetric shallow quantum well. The effective mass renormalized by relativistic corrections and Bychkov–Rashba coefficients are analytically obtained for both cases. It is demonstrated that the general solution transforms to the particular solutions, found previously (Eremko et al., 2015) with the use of spin invariants. The general solution allows to establish conditions at which a specific (accompanied or non-accompanied by Rashba splitting) spin state can be realized. These results can prompt the ways to control the spin degree of freedom via the synthesis of spintronic heterostructures with the required properties.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Exact solutions of Dirac equation on a static curved space–time

In this work we study the Dirac equation in a curved static space–time, where the metric depends on two arbitrary functions f(r) and g(r), namely the line element is ds{sup 2}=e{sup 2f(r)}dt{sup 2}-e{sup 2g(r)}dr{sup 2}-r{sup 2}dθ{sup 2}-r{sup 2}sin{sup 2}θd&straightphi;{sup 2}. We exactly solve the angular equation for arbitrary f(r) and g(r). Then, we uncouple the radial component of the spinor for a metric model defined by f=g and find the exact solutions for three quantum systems: Hydrogen atom, Dirac–Morse oscillator and a linear radial potential by V(r)=ar.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Analytic Solutions to Coherent Control of the Dirac Equation

A simple framework for Dirac spinors is developed that parametrizes admissible quantum dynamics and also analytically constructs electromagnetic fields, obeying Maxwell’s equations, which yield a desired evolution. In particular, we show how to achieve dispersionless rotation and translation of wave packets. Additionally, this formalism can handle control interactions beyond electromagnetic. This work reveals unexpected flexibility of the Dirac equation for control applications, which may open new prospects for quantum technologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Absorbing layers for the Dirac equation

This work is devoted to the construction of perfectly matched layers (PML) for the Dirac equation, that not only arises in relativistic quantum mechanics but also in the dynamics of electrons in graphene or in topological insulators. While the resulting equations are stable at the continuous level, some care is necessary in order to obtain a stable scheme at the discrete level. This is related to the so-called fermion doubling problem. For this matter, we consider the numerical scheme introduced by Hammer et al. [19], and combine it with the discretized PML equations. We state some arguments for the stability of the resulting scheme, and perform simulations in two dimensions. The perfectly matched layers are shown to exhibit, in various configurations, superior absorption than the absorbing potential method and the so-called transport-like boundary conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Stationary solutions of the Dirac equation in the gravitational field of a charged black hole

A stationary solution of the Dirac equation in the metric of a Reissner-Nordstroem black hole has been found. Only one stationary regular state outside the black hole event horizon and only one stationary regular state below the Cauchy horizon are shown to exist. The normalization integral of the wave functions diverges on both horizons if the black hole is non-extremal. This means that the solution found can be only the asymptotic limit of a nonstationary solution. In contrast, in the case of an extremal black hole, the normalization integral is finite and the stationary regular solution is physically self-consistent. The existence of quantum levels below the Cauchy horizon can affect the final stage of Hawking black hole evaporation and opens up the fundamental possibility of investigating the internal structure of black holes using quantum tunneling between external and internal states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Solution of D dimensional Dirac equation for coulombic potential using NU method and its thermodynamics properties

The analytical solution of Ddimensional Dirac equation for Coulombic potential is investigated using Nikiforov-Uvarov method. The D dimensional relativistic energy spectra are obtained from relativistic energy eigenvalue equation by using Mat Lab software.The corresponding D dimensional radial wave functions are formulated in the form of generalized Jacobi and Laguerre Polynomials. In the non-relativistic limit, the relativistic energy equation reduces to the non-relativistic energy which will be applied to determine some thermodynamical properties of the system. The thermodynamical properties of the system are expressed in terms of error function and imaginary error function.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Speed-of-light pulses in a massless nonlinear Dirac equation

Here, in this work, we explore a massless nonlinear Dirac equation, i.e., a nonlinear Weyl equation. We study the dynamics of its pulse solutions and find that a localized one-hump initial condition splits into a localized two-hump pulse, while an associated phase structure emerges in suitable components of the spinor field. For times larger than a transient time t s this pulse moves with the speed of light, effectively featuring linear wave dynamics and maintaining its shape (both in two and three dimensions). We show that for the considered nonlinearity, this pulse represents an exact solution of the nonlinear equation. Finally, we briefly comment on the generalization of the results to a broader class of nonlinearities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Phenomenological description of relativistic-quark interaction on the basis of the Dirac equation with the Cornell potential

Relativistic-quark interaction is described phenomenologically on the basis of the Dirac equation with the Cornell potential. A general form of the initial equation involving the vector and scalar components of the Cornell potential is used for the case of an arbitrary relation between them. The Hamiltonian in the Foldy–Wouthuysen representation is derived in a general form with allowance for electromagnetic interaction. In contrast to earlier investigations, it is relativistic and exact for the zeroth- and first-order terms in the Planck constant and also for those second-order terms that describe contact interactions. General quantum-mechanical equations of motion for the momentum and spin are derived, and the classical limit of the Hamiltonian and for the equations of motion is found for the first time. A relation between the angular velocity of quark spin precession and the force acting on the quark is obtained. The energy of spin–orbit interaction is rather high (on the order of 100 MeV). Terms that describe spin–orbit and contact interactions have opposite signs for the vector and scalar components of the Cornell potential. The evolution of the quark helicity and the spin–spin interaction of the quarks are also calculated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗