In this paper, we use the theory of fractional powers of linear operators to construct a general (analytic) representation theory for the square-root energy operator of relativistic quantum theory, which is valid for all values of the spin. We focus on the spin 1/2 case, considering a few simple yet solvable and physically interesting cases, in order to understand how to interpret the operator. Our general representation is uniquely determined by the Green's function for the corresponding Schrödinger equation. We find that, in general, the operator has a representation as a nonlocal composite of (at least) three singularities. In the standard interpretation, the particle component has two negative parts and one (hard core) positive part, while the antiparticle component has two positive parts and one (hard core) negative part. This effect is confined within a Compton wavelength such that, at the point of singularity, they cancel each other providing a finite result. Furthermore, the operator looks like the identity outside a few Compton wavelengths (cutoff). To our knowledge, this is the first example of a physically relevant operator with these properties.When the magnetic field is constant, we obtain an additional singularity, which could be interpreted as particle absorption and emission. The physical picture that emerges is that, in addition to the confined singularities and the additional attractive (repulsive) term, the effective mass of the composite acquires an oscillatory behavior.We also derive an alternate relationship between the Dirac equation (with minimal coupling) and the square-root equation that is much closer than the one obtained via the FoldyWouthuysen method, in that there is no change in the wave function. This is accomplished by considering the scalar potential to be a part of the mass. This approach leads to a new KleinGordon equation and a new square-root equation, both of which have the same eigenvalues and eigenfunctions as the Dirac equation. Finally, we develop a perturbation theory that allows us to extend the range of our theory to include suitable spacetime-dependent potentials.2