An embedding scheme is developed for the Dirac Hamiltonian H. Dividing space
into regions I and II separated by surface S, an expression is derived for the
expectation value of H which makes explicit reference to a trial function
defined in I alone, with all details of region II replaced by an effective
potential acting on S and which is related to the Green function of region II.
Stationary solutions provide approximations to the eigenstates of H within I.
The Green function for the embedded Hamiltonian is equal to the Green function
for the entire system in region I. Application of the method is illustrated for
the problem of a hydrogen atom in a spherical cavity and an Au(001)/Ag/Au(001)
sandwich structure using basis sets that satisfy kinetic balance.Comment: 16 pages, 5 figure