#### Abstract

We present exact solutions of the Dirac equation with Yukawa potential in the presence of a Coulomb-like tensor potential. For this goal we expand the Yukawa form of the nuclear potential in its mesonic clouds by using Taylor extension to the power of seventh and bring out its simple form. In order to obtain the energy eigenvalue and the corresponding wave functions in closed forms for this potential (with great powers and inverse exponent), we use ansatz method. We also regard the effects of spin-spin, spin-isospin, and isospin-isospin interactions on the relativistic energy spectra of nucleon. By using the obtained results, we have calculated the deuteron mass. The results of our model show that the deuteron spectrum is very close to the ones obtained in experiments.

#### 1. Introduction

It is well known that the exact energy eigenvalues of the bound state play an important role in quantum mechanics. In particular, the Dirac equation which describes the motion of a spin-1/2 particle has been used in solving many problems of nuclear and high-energy physics. Recently, there has been an increase in searching for analytic solution of the Dirac equation [1–11]. Recently, tensor couplings have been used widely in the studies of nuclear properties [12–22], and they were introduced into the Dirac equation by substitution [16, 23], where is one of the particles, mass and *ω* refers to harmonic oscillator. In this work, we are going to solve the relativistic Dirac equation for Yukawa potential in fresh approaches. We first discuss the Dirac equation with a tensor potential for Yukawa potential. Then, we solve the resulting equations for deuteron nuclei and obtain the parameters for it. We are also able to number the mass of Deuteron by using obtained parameters, and the equation is in conformity with experimental measure.

#### 2. Yukawa Potential

The screened Coulomb potential, also known as the Yukawa potential in atomic physics and the Debye-Huckel potential in plasma physics, is of interest in many areas of physics. It was originally used to model strong nucleon-nucleon interactions due to meson exchange in nuclear physics by Yukawa [19, 20]. It is also used to represent a screened Coulomb potential due to the cloud of electronic charges around the nucleus in atomic physics or to account for the shielding by outer charges of the Coulomb field experienced by an atomic electron in hydrogen plasma. However, the Schrödinger equation for this potential cannot be solved exactly; hence, various numerical and pertubative methods have been devised to obtain the energy levels and related physical quantities [24]. The generic form of this potential is given by that is the potential strength, and

In this part, we expand the Yukawa potential in its meson clouds, , by using Taylor extension to the power of seventh: We can reduce this equation to Now, we have the general Yukawa equation, that , , , , , , , and are

#### 3. Exact Analytical Solution of the Dirac Equation for Yukawa Potential

We are interested in the particular solution of the Dirac equation with the Yukawa interaction potential (4). In this section, we anticipate the answer and then obtain the parameters for highest nuclei, deuteron. A deuteron consists of a neutron and a proton. Let us assume one nucleon is fixed and another nucleon is moving around the center of mass. At last, by using single particle model, we investigate the effects of one pion exchange for a single nucleon. The Dirac Hamiltonian is where is the scalar potential, is the vector potential, and is the tensor part. We have and . The Dirac equation can be solved exactly for the cases and and which are the usual Dirac matrices. We will indicate upper and lower components of Dirac spinor with and . Expressing the Dirac matrices in terms of the Pauli matrices: we obtain the following coupled equations, [25, 26]We assume that , , and are radial potentials. Using the identity With (7a) and (7b), we find the following second-order differential equations for and : Here stands for the spin operator and for orbital angular momentum operator. The Hamiltonian operator commutes with the total angular momentum operator which is given by and with the parity operator. Therefore, the eigenfunctions defined by Hamiltonian can be expressed as where denotes the spin spherical harmonics. The quantum number is related to and as follows: For spin spherical harmonics we have Equations (9a) and (9b) reduced to a set of coupled equations for the radial wave functions and :Using the following relations [13]: for particle of mass and energy, we have We assume the case [27–31] and consider , then the upper component in (13a) is as follows: As we know in deuteron, two nucleons have in D-state, so Therefore, we have a Schrödinger-like equation [31] By setting In order to solve (18), we suppose the following form for the wave function: Now, for the functions and , respectively, we make use of the following ansatz [27, 32, 33]:

For a particular grand angular quantum number , there are different solutions which are labeled by ( determines the number of the nodes of the wave function). By substitution of and into (21) and then taking the second-order derivative of the obtained equation, we can get
We consider the ground state, which is called the 0th node solution of the above differential equation. By equating (18) and (23) for , it can be found that
The upper component is as follows:
and another component of wave function is
By substituting (25) into (26), we obtained the following equation:
Therefore, the total wave function is
Now we have to obtain parameters , , , , , , , and , by using of (5):
We know that [34], by using of (24), we can get amount of parameters , , , and *:*
The bonding energy for deuteron is reported 2.224 MeV [35]. This is useful for calculating the numeric values of parameters in (30):
By substituting the amounts of (31) and , we obtained the numeric values of parameters in (30). The parameters are reported in Table 1, while the parameters of Yukawa potential are reported in Table 2.

#### 4. The Hyperfine Interaction

As the Baryons have spin and isospin, the deuteron neutron and proton can interact together. We introduce the hyperfine interaction potential. The complete potential is

In this work, we have added the hyperfine interaction potentials (, , and ) which yield properties very close to the experimental results. By regarding as the nonpertubative potential and the other terms in (32) as pertubative ones according to this explanation. The spin and isospin potential contains a -like term. We have modified it by a Gaussian function of the nucleons pair relative distance [36]
where is the spin operator of the th nucleon and is the relative nucleon pair coordinate*. * and are constants: and [34, 37]. We know that the deuteron consists of two nucleons. Where 1 indicate the first nucleon, and 2 indicate the second nucleon. We have
Furthermore, we add two hyperfine interaction terms to the Hamiltonian nucleon pairs similar to (33) [38]
where is the isospin operator of the th nucleon, and the constants are described as and [37, 39]. The another one is a spin-isospin interaction, given by similar (33) [38, 40]
The fitted parameters are and [35, 41], where is the perturbed wave function, and we write it as
The deuteron mass are given by two nucleons masses, and the eigenenergies of the Dirac equation , ( is a function of , , , , ) with the first-order energy correction from potential can be obtained by using the unperturbed wavefunction (33), (35), and (36). The total potential for the ground state as well as the other states can be written as [40, 42]
We first assume . The potentials can be extracted from (33), (35), and (36):
Now, we can calculate the deuteron mass according to the following formula:
By substituting , energy eigenvalue (), and the expectation values of in (40) we have

By comparing the experimental amount of deuteron mass () with our calculated mass for deuteron, we found that a good agreement has been obtained by our model [35].

#### 5. Conclusion

In this work, we offer a new approach for the solving of Yukawa potential. We expand this potential around of its mesonic cloud that gets a new form with great powers and inverse exponent. We calculate wave function of the Dirac equation for a new form of Yukawa potential, including a coulomb-like tensor potential. By considering the effects of pertubative interaction potential, we see that the theoretical and experimental masses are in complete agreement. These improvements in reproduction of deuteron mass is obtained by using a suitable form for confinement potential and exact analytical solution of the Dirac equation for our proposed potential. Finally, one can use this model for another nuclei, by relativistic or nonrelativistic equation and get the properties of them and the amount of “”, the strength of nuclear force.