Two-body Dirac equation approach to the deuteron

The two-body Dirac(Breit) equation with potentials associated to one-boson-exchanges with cutoff masses is solved for the deuteron and its observables calculated. The 16-component wave-function for the Jπ = 1+ state contains four independent radial functions which satisfy a system of four coupled differential equations of first order. This system is numerically integrated, from infinity towards the origin, by fixing the value of the deuteron binding energy and imposing appropriate boundary conditions at infinity. For the exchange potential of the pion, a mixture of direct plus derivative couplings to the nucleon is considered. We varied the pion-nucleon coupling constant, and the best results of our calculations agree with the lower values recently determined for this constant.