Abstract
We consider electrons in the presence of interfaces with different effective electron mass, and electromagnetic fields in the presence of a highpermittivity interface in bulk material. The equations of motion for these dimensionally hybrid systems yield analytic expressions for Green’s functions and electromagnetic potentials that interpolate between the twodimensional logarithmic potential at short distance, and the threedimensional r^{−1} potential at large distance. This also yields results for electron densities of states which interpolate between the wellknown twodimensional and threedimensional formulas. The transition length scales for interfaces of thickness L are found to be of order Lm/2m_{*} for an interface in which electrons move with effective mass m_{*}, and for a dielectric thin film with permittivity in a bulk of permittivity . We can easily test the merits of the formalism by comparing the calculated electromagnetic potential with the infinite series solutions from image charges. This confirms that the dimensionally hybrid models are excellent approximations for distances r ≳ L/2.
Keywords:
Density of states; Coulomb and exchange interactions in nanostructures; Dielectric thin filmsIntroduction
When we suppress motion of particles in certain directions through confining potentials, e.g. in quantum wells or quantum wires, we often model the residual low energy excitations in the system through lowdimensional quantum mechanical systems. Prominent examples of this concern layered heterostructures, and one instance where the number d of spatial dimensions enters in a manner which is of direct relevance to technology is in the density of states. In the standard parabolic band approximation, this takes the form (with two helicity or spin states)
These are densities of states per ddimensional volume and per unit of energy. The corresponding dependence of the relation between the Fermi energy and the density n of electrons on d is
Variants of these equations (including summation over subbands) are often used for d = 2 or d = 1 to estimate carrier densities in quasi twodimensional systems or nanowires, and the density of states plays a crucial role in all transport and optical properties of materials. Indeed, the obvious relevance for electrical conductivity properties in micro and nanotechnology implies that densities of states for d = 1, 2, or 3 are now commonly discussed in engineering textbooks, but there is another reason why I anticipate that variants of Eq. (1) will become ever more prominent in the technical literature. Densities also play a huge role in data storage, but with us still relying on binary logic switching between two stable states (spin up or down, charge or no charge, conductivity or no conductivity), data storage densities are limited by the physical densities of the systems which provide the dual states. We could (and likely will) drive information technology and integration much further if we can find ways to utilize more than just two states of a physical system to store and process information. Then, data storage densities should become proportional to energy integrals of local densities of states. Equation (1) for d = 1 or d = 2 is certainly applicable for particles which have low energies compared to the confinement energy of a nanowire or a quantum well, but how can we effectively model particles which are weakly confined to a nanowire or quantum well, or which are otherwise affected by the presence of a lowdimensional substructure? In these cases, we can devise dimensionally hybrid models [1,2] which yield e.g. densities of states which interpolate between d = 2 and d = 3 [3,4]. This construction will be reviewed in Sect. 2. Based on the experience gained with dimensionally hybrid Hamiltonians for massive particles, we can also construct interdimensional Hamiltonians for photons which should be applicable to photons in the presence of highpermittivity thin films or interfaces. These models can also be solved in terms of infinite series expansions using image charges, and the merits of this approach can easily be tested. The case of highpermittivity thin films and testing the theory against image charge solutions will be discussed in Sect. 3.
Dimensionally Hybrid Hamiltonians and Green’s Functions for Massive Particles in the Presence of Thin Films or Interfaces
We use the connection between Green’s functions and the density of states to generalize Eq. (1) for massive particles in the presence of a thin film or interface.
The energydependent Green’s function for a Hamiltonian H with spectrum E_{n} and eigenstates n, ν〉 is
Here, ν is a degeneracy index and the notation implies that continuous components in the indices (n, ν) are integrated. The first equation simply states the relation between the resolvent of the Hamiltonian and the Green’s function G(E) which is normalized as lim_{m→0,E→0}G(E)_{d=3} = (4πr)^{−1}.
The zeroenergy Green’s function G(0) determines e.g. 2particle correlation functions and electromagnetic interaction potentials, and the energydependent Green’s function G(E) determines e.g. scattering amplitudes for particles of energy E. Application for resistivity calculations is therefore another technologically relevant application of Green’s functions. However, in the present section we are interested in this function because it also determines the local density of states in a system with Hamiltonian H through the relation
Here, we explicitly included a factor 2 for the number of spin or helicity states, because the summation over degeneracy indices in (3,4) usually only involves orbital indices.
For our present investigation, the distinctive feature of the interface is that the particles move in it with an effective mass m_{*}, while their mass in the surrounding bulk is m. We use coordinates parallel to a plane interface, which is located at z = z_{0}. Bold vector notation is used for quantities parallel to the interface, e.g. and .
We assume that the interface has a thickness L. If the wavenumber component orthogonal to the interface is small compared to the inverse width, k_{⊥}L ≪ 1, i.e. if the de Broglie wavelength and the incidence angle satisfy λ ≫ 2πLcosϑ, we can approximate the kinetic energy of the particles through a second quantized Hamiltonian
where μ = m_{*}/L. The corresponding first quantized Hamiltonian is
The interesting aspect of the Hamiltonians (5,6) is the linear superposition of twodimensional and threedimensional kinetic terms. The formalism presented here could and will certainly be extended to include also kinetic terms which are linear in derivatives, in particular in the interface term. This would be motivated either by a Rashba term arising from perpendicular fields penetrating the interface [511] of from the dispersion relation in Graphene [1215]. However, for the present investigation we will use a parabolic band approximation in the bulk and in the interface.
The energydependent Green’s function describes scattering effects in the presence of the interface but also applies to scattering off perturbations which are not located on the interface. In an axially symmetric mixed representation
the first order approximation to scattering of an orthogonally incoming plane wave off an impurity potential
corresponds to
Green’s functions for surfaces or interfaces are commonly parametrized in an axially symmetric mixed representation like . In braket notation, this corresponds for the free Green’s function G_{0}(E), which is also translation invariant in z direction, to
We will briefly recall the explicit form of the free Green’s function G_{0}(E) in the axially symmetric mixed parametrization for later comparison. The equation
yields
To study how this is modified in the presence of the interface, we observe that the Hamiltonians (5) or (6) yield a Schrödinger equation
The corresponding equation for the Green’s function or 2point correlation function is
The solution of this equation is described in the Appendix. In particular, we find the representation (see Eq. (27))
where the definition ℓ≡m/2μ = Lm/2m_{*} was used. The ℓindependent terms in (10) correspond to the free Green’s function G_{0}(E) (8).
The interface at z_{0} breaks translational invariance in z direction, and we have with Eq. (7)
We will use the result (10) to calculate the density of states in the interface. Substitution yields
and after evaluation of the integral
This is a more complicated result than the density (1) for d = 2 or d = 3. However, it reduces to either the twodimensional or threedimensional density of states in the appropriate limits, see Fig. 1. For large energies, i.e. if the states only probe length scales smaller than the transition length scale ℓ, we find the twodimensional density of states properly rescaled by a dimensional factor to reflect that it is a density of states per threedimensional volume,
Figure 1. The red line is the twodimensional limit (12). The blue line is the threedimensional density of states. The it black line is the interdimensional density of states (11) for ℓ = 50 nm
For small energies, i.e. if the states probe length scales larger than ℓ, we find the threedimensional density of states
This limiting behavior for interpolation between two and three dimensions is consistent with what is also observed for the zeroenergy Green’s function in the interface, see equations (21–22) below.
Equation (11) also implies interpolating behavior for the relation between electron density and Fermi energy on the interface. The full relation is
This approximates twodimensional behavior for ,
and threedimensional behavior for ,
It is intuitively understandable that the presence of a layer reduces the available density of states for given energy, or equivalently increases the Fermi energy for a given density of electrons. The presence of a layer generically implies boundary or matching conditions which reduce the number of available states at a given energy.
A condition for relevance of the interdimensional behavior is a large transition scale compared to the layer thickness, ℓ ≫ L, see also Fig. 2. In terms of effective particle mass, this means
i.e. the energy band in the interface should be more strongly curved than in the bulk matrix for the transition to twodimensional behavior to be observable.
Figure 2. The upper dotted (blue) line is the threedimensional Green’s function (4πr)^{−1} in units of ℓ^{−1}, the continuous line is the Green’s function (19) in units of ℓ^{−1}, and the lower dotted (red) line is the twodimensional logarithmic Green’s function ℓ·G = − (γ + ln(r/2ℓ))/(4π)
Electric Fields in the Presence of HighPermittivity Thin Films or Interfaces
The zeroenergy Green’s function determines electrostatic and exchange interactions through the electrostatic potential . Here, q is an electric charge in a dielectric material of permittivity . The zeroenergy Green’s function in d spatial dimensions is given by
We cannot infer from the previous section that the zero energy limit of the interdimensional Green’s function calculated there also yields a dimensionally hybrid potential, because we were dealing with solutions of Schrödinger’s equation instead of the Gauss law. However, we can rederive the zero energy limit of that Green’s function from the Gauss law for electromagnetic fields in the presence of a highpermittivity interface.
Suppose we have charge carriers of charge q and mass m in the presence of an interface with permittivity and permeability μ_{*}, We continue to denote vectors parallel to the interface in bold face notation, , , etc.
If the photon wavelengths and incidence angles satisfy the condition λ ≫ 2πLcosϑ, we can approximate the system with an action
Variation with respect to the electrostatic potential, , yields the Gauss law in the form
and the continuity condition E_{z}(z_{0} − 0) = E_{z}(z_{0} + 0).
We solve Eq. (16) in Coulomb gauge,
where the Green’s function has to satisfy
This equation is the zero energy limit of Eq. (9) with the substitution
We can therefore read off the solution from the results of the previous section with E = 0 and now .
Equation (10) yields in particular
with . Fourier transformation yields
The zeroenergy Green’s function in the interface is given in terms of a Struve function and a Neumann function^{1},
This yields logarithmic behavior of interaction potentials at small distances r ≪ ℓ and 1/r behavior for large separation r ≫ ℓ of charges in highpermittivity thin films,
see also Fig. 2.
For the comparison with image charges, we set z_{0} = 0 and recall that the solution for the potential of a charge q at , z = 0 proceeds through the ansatz
and symmetric continuation to z < −L/2.
This yields electric fields
and the junction conditions at z = L/2 yield for n ≥ 0 from the continuity of E_{r},
and from the continuity of D_{z},
These conditions can be solved through
In particular, the potential at z = 0 is
We have
and therefore for
The solution from image charges is in very good agreement with the analytic model for distances r ≳ L/2, where both the image charge solution and the analytic model show strong deviations from the bulk r^{−1} behavior. This is illustrated in Fig. 3 by plotting the reduced electrostatic potential for a charge q, in the interface.
Figure 3. Different reduced electrostatic potentials are plotted for . The upper dotted (green) line is the threedimensional reduced potential L/(4πr). The central dotted (blue) line is the reduced potential following from the image charge solution (22). The solid (black) line is the potential from the analytic model (19). The lower dotted (red) line is the reduced logarithmic potential. The reduced potentials from our analytic model and from image charges are indistinguishable for , see also Fig. 4
It is also instructive to plot the relative deviation between the dimensionally hybrid potential which follows from (20) and the potential (23) from image charges.
Figure 4 shows that for r ≳ L/2, the dimensionally hybrid model is a very good approximation to the potential from image charges with accuracy better than 10^{−2} if . For , the accuracy is still better than 4 × 10^{−2}.
Figure 4. The relative deviation between the dimensionally hybrid potential from (19) and the potential (22) from image charges for
Summary
An analysis of models for particles in the presence of a low effective mass interface, and for electromagnetic fields in the presence of a highpermittivity thin film, yields dimensionally hybrid densities of states (11) and electrostatic potentials (17,20) which interpolate between twodimensional behavior and threedimensional behavior. The analytic model for the electromagnetic fields is in very good agreement with the infinite series solution already for small distance scales r ≳ L/2, where the potential strongly deviates from the standard bulk r^{−1} potential. At distance scales smaller than L/2, r^{−1}, behavior seems to dominate again for the electrostatic potential, in agreement with expectations that for distances which are small compared to the lateral extension of a dielectric slab, bulk behavior should be restored. However, note that neither the interdimensional analytic model nor the solution from image charges is trustworthy for very small distances, because both models rely on a continuum approximation through the use of effective permittivities, but the continuum approximation should break down at subnanometer scales.
The most important finding is that interfaces and thin films of width L should exhibit transitions between twodimensional and threedimensional distance laws for physical quantities at length scales of order Lm/2m_{*} or , respectively. Interfaces with strong band curvature or high permittivity should provide good samples for experimental study of the transition between twodimensional and threedimensional behavior.
Appendix: Solution of Eq. 9
Substitution of the fourier transform
into Eq. 9 yields
This yields with (7) the condition
Fournier transformation with respect to z yields
This result implies that has the form
with the yet to be determined satisfying
For the treatment of the integrals, we should be consistent with the calculation of the free retarded Green's function (8),
This yields
And therefore
where the definition ℓ⊝m/2μ = Lm/2m*. Fourier transformation of Eq (26) with respect to k_{⊥} yields finally
The Green's function with only k space variables is found from the Fournier transform of Eq. 25,
and the ensuing equations
This yields
It is easily verified that Fournier transformation yields again the result (26).
Footnotes
^{1}Our notations for special functions follow the conventions of Abramowitz and Stegun [16].
Acknowledgements
This research was supported by NSERC Canada.
Open Access
This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.
References

Int. J. Theor. Phys.. 2003, 42:569. Publisher Full Text

Nanoscale Res. Lett.. 2008, 3:140.
Bibcode number [2008NRL.....3..140D]
Publisher Full Text 
Phys. E. 2008, 40:524. Publisher Full Text

Phys. E. 2008, 40:2973.
COI number [1:CAS:528:DC%2BD1cXntVOns7o%3D]
Publisher Full Text 
JETP Lett.. 1984, 39:78.
Bibcode number [1984JETPL..39...78B]

J. Phys. C. 1984, 17:6039.
Bibcode number [1984JPhC...17.6039B]
Publisher Full Text 
Cappelluti E, Grimaldi C, Marsiglio F:
Phys. Rev. Lett.. 2007, 98:167002.
COI number [1:STN:280:DC%2BD2s3ptVyqsA%3D%3D]; Bibcode number [2007PhRvL..98p7002C]
PubMed Abstract  Publisher Full Text 
Cappelluti E, Grimaldi C, Marsiglio F:
Phys. Rev. B. 2007, 76:085334.
Bibcode number [2007PhRvB..76h5334C]
Publisher Full Text 
Srisongmuang B, Pairor P, Berciu M:
Phys. Rev. B. 2008, 78:155317.
Bibcode number [2008PhRvB..78o5317S]
Publisher Full Text 
Phys. E. 2008, 40:1729. Publisher Full Text

Nanoscale Res. Lett.. 2009, 4:178.
COI number [1:CAS:528:DC%2BD1MXht1Kntrw%3D]; Bibcode number [2009NRL.....4..178L]
Publisher Full Text 
Phys. Rev. Lett.. 1984, 53:2449.
Bibcode number [1984PhRvL..53.2449S]
Publisher Full Text 
Apalkov V, Wang XF, Chakraborty T:
Mod. Phys. B. 2007, 21:1165.
COI number [1:CAS:528:DC%2BD2sXltl2jtb0%3D]
Publisher Full Text 
Phys. Rev. Lett. 2008, 100:256405.
Bibcode number [2008PhRvL.100y6405C]
PubMed Abstract  Publisher Full Text 
Nanoscale Res. Lett.. 2010, 5:169.
Bibcode number [2010NRL.....5..169L]
Publisher Full Text 
Abramowitz M., Stegun I.A. (Eds): Handbook of Mathematical Functions 9th printing. Dover Publications, New York; 1970.