The one dimensional harmonic oscillator or square well, for example, for which the energy -dependent Green function G0(x,x′;E) is known, have been taken for many years as solvable models for semi-conductor quantum wells [1]. Frequently delta function potentials are placed at various points to simulate defects or impurities. In the case of a single impurity potential, V(x) = λδ(x−a), the Green function for the composite system is known to be [2]

In this note a corresponding formula is derived for the case V(x) = λδ(x−a)+μδ(x−b) On the basis of the analogy of the algebraic structure the result is extended to N-impurities and for a standard interpretation of the Dirac delta function, to higher dimension.

We first note that the same argument can be used for the time dependent-, as well as the energy-dependent Green functions, so we shall omit the third argument and write simply G(x, x′).

Beginning with the Dyson equation, noting that G₀(x, y) = G₀(y, x)

where the integration extends over the system domain, one has the set of equations

The linear Equations (4) and (5) are easily solved for G(a, x′) and G(b, x′):

with

By inserting (6) and (7) into (3) we obtain the desired expression

By setting μ to 0 (9) reduces to (1), proving this expression as well. The most salient feature of (9) is the denominator D whose zeros form the exact spectrum of the composite system. For example, when a and b coincide, D reduces to 1−(λ+μ)G(a, a) and (9) reduces to (1) with λ replaced by the amplitude λ+μ. I.e., the two impurities combine to form one with combined amplitude. This generalizes the result of Fasssari and Rinaldi [3], for two identical defects symmetrically placed with respect to the center of a harmonic oscillator. An expression similar to (9) has been derived recently by Horing (private communication) for the case of a quantum dot in a magnetic field.

Two further points can be made. Nothing in the derivation of (9) restricts it to the line. If we accept the standard definition δ(x→)=Πj=1dδ(xj), then (9), and its consequences, are valid for d-dimensional quantum systems. This has been proven function-theoretically for the three dimensional quantum dot with two symmetrically placed identical impurities by Albeverio et al. [4].

A second observation is that D is simply the Cramer determinant for the pair of simultaneous linear Equations (4) and (5). In the case of impurity potential

there will be N such equations and the determinant is easily evaluated. The general result is

If a quantum system having Green function G₀(x, y) is decorated with N delta function impurities λ_jδ(x−a_j), j = 1, 2, ⋯N, then the new energy levels are the roots of

where A_lm = δ_lm−λ_lG₀(a_l, a_m).

Thus,

which reduces to the N = 1 and N = 2 cases appropriately and shows that any two coinciding impurities coalesce as indicated above.

Note that if all the λs and a's coincide then

Equation (11) might offer a new approach to Kronig-Penney-type systems for periodic or random unit cells.

Finally, it should be pointed out that the work in this note is paralleled in the theory of quantum graphs introduced by Linus Pauling about 1930 to describe electrons in molecules which has developed into a sophisticated and important branch of quantum physics [5]. For relations of this discipline to the present work see the papers by Andrade et al. [6] and Andrade and Severini [7].

The author confirms being the sole contributor of this work and has approved it for publication.