Boundary Conditions for Sturm-Liouville Equation with Transition Regions and Barriers or Wells

Alfred Wünsche

Institut für Physik, Humboldt-Universität, Nonclassical Radiation, Berlin, Germany.

DOI: 10.4236/apm.2021.114018                

Abstract

By means of expansions of rapidly in infinity decreasing functions in delta functions and their derivatives, we derive generalized boundary conditions of the Sturm-Liouville equation for transitions and barriers or wells between two asymptotic potentials for which the solutions are supposed as known. We call such expansions “moment series” because the coefficients are determined by moments of the function. An infinite system of boundary conditions is obtained and it is shown how by truncation it can be reduced to approximations of a different order (explicitly made up to third order). Reflection and refraction problems are considered with such approximations and also discrete bound states possible in nonsymmetric and symmetric potential wells are dealt with. This is applicable for large wavelengths compared with characteristic lengths of potential changes. In Appendices we represent the corresponding foundations of Generalized functions and apply them to barriers and wells and to transition functions. The Sturm-Liouville equation is not only interesting because some important second-order differential equations can be reduced to it but also because it is easier to demonstrates some details of the derivations for this one-dimensional equation than for the full three-dimensional vectorial equations of electrodynamics of media. The article continues a paper that was made long ago.

Keywords

Schrödinger Equation, Drude Approximation, Transition Layer, Potential Barrier, Potential Well, Reflection, Refraction, Moment Series, Generalized Functions, Delta Function and Its Derivatives, Discrete or Bound Eigenstates

1. Introduction

It is well-known that differential equations in physics with discontinuous material properties can be solved by applying boundary conditions on the interface between the different media for the adaption of the primary independent solutions for the media. When the transitions between the different media are sufficiently abrupt, but still smooth, one can often solve such problems approximately by generalized boundary conditions. We will demonstrate this here for one of the simplest second-order differential equations in physics, the Sturm-Liouville equation. This equation is of great importance because it is almost identical to the one-dimensional Schrödinger equation and some special problems of wave propagation in optics can be also reduced to this equation.

Generalized boundary conditions in optics taking into account transition layers between two different media were first derived by Drude [1] [2]. Three additional

parameters (of the kind of moments) $l={\displaystyle {\int }_1^2}\text{d}z$, $p={\displaystyle {\int }_1^2}\text{d}z\epsilon$, $q={\displaystyle {\int }_1^2}\frac{\text{d}z}{\epsilon }$ ( $\epsilon$

permittivity) play a role, where the integration goes over the boundaries of “thin” transition layers with thickness l which, however, is difficult to define exactly. In reflection and refraction formulae they often go at once into numerator and denominator and partially compensate the incorrectness connected with this approach. They were applied, for example, in the monographs about thin films of Mayer [3], Heavens [4], Rozenberg [5], and Vašiček [6].

A relatively new procedure to derive generalized boundary conditions is to involve generalized functions and, in particular, the step function and all its derivatives which are the delta function with its derivatives. This method is appropriate in classical electrodynamics for the derivation of generalized boundary conditions for media with spatial dispersion [7] (including, e.g., optical activity) and for correct derivations in presence of transition layers. A preliminary step to such derivations of boundary conditions was made in [8] [9]. The method was also developed in [10] [11] in application to the optics of media with transition layers, and the reflection was calculated in the first-order approximation of the generalized boundary conditions. In our paper [12] (see also Acknowledgement) we developed this method for light propagation in media with spatial dispersion (e.g., optical activity or gyrotropy) and transition layers and tried to represent the general procedure of continuing it to higher-order approximations. The main role played expansions of rapidly in infinity decreasing functions into series of delta functions and their derivatives which we call moment series because their coefficients are moments of the expanded function. The further development in the field of reflection and refraction not considered here can be seen, in particular, in the monograph of Lekner [13]. It is fairly comprehensive and represents many different aspects of reflection and refraction in optics and acoustics. In chap. 3 it uses integral invariants which are partially analogous to the parameters which we call moments of a function and in chap. 18 it considers chiral media including boundary conditions (a special case of spatial dispersion [7] which in the new edition of vol. 8 of Landau and Lifshits [14], chap. XII and also in our paper [12] appear under this last name with one of its special cases as natural optical activity).

In the present article we try to demonstrate the method of derivation of generalized boundary conditions in detail at a simple one-dimensional equation and as such we have chosen the Sturm-Liouville equation (e.g., special monographs [15] [16] [17] and many books of functional analysis) which is nearly the one-dimensional Schrödinger equation and to which many other equations can be reduced (e.g., article [18] in [19] ). We need for this purpose exact solutions for asymptotic potentials of the one-dimensional Schrödinger equation (e.g., [20] [21] [22] [23] and others). As said it is easier to make and to explain for this equation the derivations of generalized boundary conditions between two asymptotic potentials up to high-order approximations than for the full Maxwell equations with transition layers between two media and with spatial dispersion in the constitutive equations.

In Section 2 we discuss the restrictions for the applicability of the approach to potential problems for the one-dimensional Schrödinger equation. From Section 3 on we develop systematically the approach of the derivation of boundary conditions for the Sturm-Liouville equation between two asymptotic potentials with solutions supposed as known, in principle. In Section 4 we show how we can obtain finite approximations of the infinite system of boundary conditions and explain the special new situation when we go to the third-order approximation. From Section 5 on we apply this to the special case of two different but asymptotically constant potentials and calculate reflection and refraction problems in Section 6 and discuss in Sections 7 and 8 the appearance of discrete eigenstates which are in certain sense equivalent to total reflection and to surface waves for the three-dimensional Maxwell equations. In Section 9 we give a matrix form of the boundary condition which is well appropriate for a layered system where the output of one layer may serve as input of the next layer. The Appendix A gives a detailed derivation of the moment series in delta functions and their derivatives with moments of the expanded functions as coefficients and some related formulae for Generalized functions. Appendix B brings special examples for such expansions which are interesting for potential barriers and wells and for transition between potentials.

2. One-Dimensional Schrödinger Equation for Potential Functions and the de-Broglie Wavelength

The one-dimensional time-dependent Schrödinger equation for a particle with mass m in a state $|\Psi \left(t\right)\rangle$ in abstract representation possesses the form

$\text{i}\hslash \frac{\partial }{\partial t}|\Psi \left(t\right)\rangle =H|\Psi \left(t\right)\rangle \mathrm{,}\left(H\equiv T\left(P\right)+U\left(Q\right)\right)\mathrm{,}$ (2.1)

where H is the Hamilton operator which is the sum of the operator of kinetic energy $T(\; P\; )$

$T\left(P\right)=\frac{P^2}{2m}\mathrm{,}$ (2.2)

and of potential energy $U\left(Q\right)$. With P and Q we denote a pair of canonical coordinates which in quantum mechanics obey the commutation relations

$\left[Q\mathrm{,}P\right]\equiv QP-PQ=\text{i}\hslash I\mathrm{,}$ (2.3)

where I is the identity operator in representation (Hilbert) space. In representation of the (non-normalizable) eigenstates $|q\rangle$ of the operator Q according to

$\begin{array}{l}Q|q\rangle =q|q\rangle \mathrm{,}\langle q|Q=q\langle q|\mathrm{,}P|q\rangle =\text{i}\hslash \frac{\partial }{\partial q}|q\rangle \mathrm{,}\\ \langle q|P=-\text{i}\hslash \frac{\partial }{\partial q}\langle q|\mathrm{,}\langle q^{\prime }|q\rangle =\delta \left(q-q^{\prime }\right)\mathrm{,}\end{array}$ (2.4)

we find from (2.1) the time-dependent Schrödinger equation in q-representation¹

$\text{i}\hslash \frac{\partial }{\partial t}\Psi \left(q\mathrm{,}t\right)=\left\{-\frac{{\hslash }^2}{2m}\frac{{\partial }^2}{\partial q^2}+U\left(q\right)\right\}\Psi \left(q\mathrm{,}t\right)\mathrm{,}\Psi \left(q\mathrm{,}t\right)\equiv \langle q|\Psi \left(t\right)\rangle \mathrm{.}$ (2.5)

With the supposition

$\Psi \left(q\mathrm{,}t\right)=\exp \left(-\text{i}\frac{Et}{\hslash }\right){\psi }_0\left(q\mathrm{,}E\right)\mathrm{,}$ (2.6)

where E is the energy of the particle one obtains from (2.1) the time-independent Schrödinger equation

$\left\{-\frac{{\hslash }^2}{2m}\frac{{\partial }^2}{\partial q^2}+U\left(q\right)\right\}{\psi }_0\left(q\mathrm{,}E\right)=E{\psi }_0\left(q\mathrm{,}E\right)\mathrm{.}$ (2.7)

It will be taken as our starting equation for the derivation of boundary conditions for the approximate solution of (one-dimensional) potential problems.

First, we make an estimation of the order of magnitudes with which we have to do in calculations with the Schrödinger Equation (2.7) if q is a spatial coordinate. Starting from Einstein’s equivalence of mass m and energy $E=m{c}^2$ de Broglie [24] conjectured that particles due to the dualism of waves and particles should possess also a wavelength which in analogy to Planck’s and Einstein’s relations for the energy E of particles and light is given by the momentum $p$ together with the dispersion equation and that according to hypothesis it is determined by (c light velocity, $m_0$ rest mass, $E_0$ rest energy)

$E=\hslash \omega ,p=\hslash k,E^2-c^2{p}^2={\left(m_0{c}^2\right)}^2,E_0=m_0{c}^2.$ (2.8)

Herein $\omega$ is the (circular) frequency, $k$ the wave vector and $\hslash$ Planck’s action quantum (divided by $2\pi$ )². From classical theory of a free particle it is

known and from relativistic theory for small velocities $\frac{v^2}{c^2}=\frac{c^2{p}^2}{E_0^2}\le 1$ (or $E\approx E_0$ ) it follows that its kinetic energy is

$E_{\text{kin}}\equiv E-E_0\approx \frac{p^2}{2m_0}=\frac{{\hslash }^2{k}^2}{2m_0}\mathrm{,}$ (2.9)

and from this dispersion equation resolved to $\left|p\right|=\hslash \left|k\right|=\hslash \frac{2\pi }{\lambda }=\sqrt{2m_0{E}_{\text{kin}}}$ one obtains for the wavelength $\lambda$ of the moving free particles with velocity $v\equiv \left|v\right|$ and momentum $p=m_{0}v$ and for frequency $\nu$

$\lambda \approx \frac{2\pi \hslash }{\sqrt{2m_0{E}_{\text{kin}}}}=\left(\frac{h}{m_{0}c}\right)\sqrt{\frac{m_0{c}^2}{2E_{\text{kin}}}}=\frac{h}{m_{0}v}\mathrm{,}\nu \approx \frac{m_0{c}^2}{2\pi \hslash }\mathrm{.}$ (2.10)

From this dispersion relations follows a phase velocity $\frac{\omega }{\left|k\right|}$ larger than c (is not a signal velocity)

$\frac{\omega }{\left|k\right|}\equiv \lambda \nu \approx c\sqrt{\frac{E_0}{2E_{\text{kin}}}}\gg c\mathrm{,}$ (2.11)

and from (2.8) the group velocity $v\equiv \frac{\partial \omega }{\partial k}$ for the free particle waves equal to particle velocity $v$

$\frac{\partial \omega }{\partial k}=\frac{\partial E}{\partial p}=\frac{c^{2}p}{E}=c\frac{ck}{\omega }\approx c\left(\frac{h}{m_{0}c}\right)\frac{k}{2\pi }=c\frac{v}{c}\frac{k}{\left|k\right|}\ll c\frac{k}{\left|k\right|}.$ (2.12)

The quantity $\frac{h}{m_{0}c}$ is the Compton length of the particle and $m_0{c}^2$ the rest

energy of the particle. For an electron with mass $m_{\text{e}}$ the following values are well known from experiment

$\frac{h}{m_{\text{e}}c}\approx 2.43\times {10}^{-10}\text{cm}\mathrm{,}{m}_{\text{e}}{c}^2\approx 0.511\text{MeV}=0.511\times {10}^6\text{eV}\mathrm{.}$ (2.13)

Therefore, the de-Broglie wavelength for a free electron with kinetic energy $E_{\text{kin}}$ is approximately

${\lambda }_{\text{e}}=\frac{1}{\sqrt{2}}\left(\frac{h}{m_{\text{e}}c}\right)\sqrt{\frac{m_{\text{e}}{c}^2}{E_{\text{kin}}}}\approx 1.23\times {10}^{-7}\sqrt{\frac{1\text{eV}}{E_{\text{kin}}}}\text{cm}\mathrm{.}$ (2.14)

For comparison, the first Bohr radius is $a_0=\frac{{\hslash }^2}{m_{\text{e}}{e}^2}\approx 5.29\times {10}^{-9}\text{cm}$ (e elementary charge), and the Rydberg constant $R_{\text{e}}$ (ionization potential) for infinitely high mass of the proton is $R_{\text{e}}=\frac{1}{2}{\alpha }^2{m}_{\text{e}}{c}^2\approx 13.6\text{eV}\mathrm{<mo>\; ;\; </mo>}\left(\alpha \approx 1/137\right)$. For energies of the electron up to the order of a few eV the wavelengths are large in comparison to the first Bohr radius. Furthermore, 1 eV is equivalent to energy ${\kappa }_{\text{B}}T$ corresponding to temperature T of 11,606 K ( ${\kappa }_{\text{B}}$ Boltzmann constant).

The approximate boundary conditions which we derive in the following for the Schrödinger Equation (2.7) are applicable if the wavelengths ${\lambda }_{\text{e}}$ (more generally, $\lambda$ ) are large in comparison to the order of lengths where the potential function $U\left(q\right)$ changes essentially. Oppositely, if the wave length of the particle is small in comparison to the order of lengths where the changes of the potential function play a role then another well-known approximation, the WKB-method³ or quasi-classical approximation, can be applied (e.g., [20] ). For high kinetic energies $E_{\text{kin}}$ of the order of the rest energy of the electron of about 0.5 MeV and higher the (non-relativistic) Schrödinger equation is no more relevant.

We usually consider in this paper Equation (2.7) in an infinite interval. If $U\left(q\right)$ in (2.7) is larger than $U\left(\pm \infty \right)$ then one speaks about a potential barrier and if it is smaller than $U\left(\pm \infty \right)$ then about a potential well. Both and more general cases which are not to classify in such simple way are included into the following calculations.

With the substitutions which preserve the coordinate as a quantity with dimension of a length

$q\equiv z\mathrm{,}U\left(q\right)\equiv \frac{{\hslash }^2}{2m}u\left(z\right)\mathrm{,}E\equiv \frac{{\hslash }^2}{2m}\lambda \mathrm{,}{\psi }_0\left(q\mathrm{,}E\right)\equiv \psi \left(z\mathrm{,}\lambda \right)\mathrm{,}$ (2.15)

Equation (2.7) can be transformed to the more neutrally looking form ( $\lambda$ now spectral parameter)

$\left\{-\frac{{\partial }^2}{\partial z^2}+u\left(z\right)-\lambda \right\}\psi \left(z\mathrm{,}\lambda \right)=0.$ (2.16)

This is the standard form of the Sturm-Liouville equation to which also some other differential equations of physics can be reduced as mentioned in the Introduction. An important example is the wave propagation in one-dimensional inhomogeneous media which is described by Landau and Lifshits [14] (chap. X, section 88). The case of polarization of the electric field $E\left(z\mathrm{<mo>\; ;\; </mo>}{k}_{\perp }\mathrm{,}\omega \right)$ perpendicular to the incidence plane is more directly connected with Equation (2.16)

(substitutions: $u\left(z\right)\to -\frac{{\omega }^2}{c^2}\epsilon \left(z\mathrm{<mo>\; ;\; </mo>}\omega \right)\mathrm{<mo>\; ;\; </mo>}\lambda \to -k_{\perp }^2$ ) than the case of polarization in

the incidence plane which requires more subtile substitutions. One has to remove then the arising first-order derivatives in the equation [14] [15] [16]. Central-symmetric potentials in two and three dimensions with transition regions can be dealt with as quasi-one-dimensional cases of the Schrödinger equation and can be transformed to a Sturm-Liouville equation with additional terms to the potential function (effective potentials) and possible applications to scattering problems. We cannot deal with this here.

If ${\psi }_1\left(z\mathrm{,}\lambda \right)$ and ${\psi }_2\left(z\mathrm{,}\lambda \right)$ are two solutions of the Sturm-Liouville equation to the same potential $u\left(z\right)$ and the same parameter $\lambda$ then in known way one derives the conservation theorem

$\frac{\partial }{\partial z}\left({\psi }_1\left(z\mathrm{,}\lambda \right)\frac{\partial }{\partial z}{\psi }_2\left(z\mathrm{,}\lambda \right)-{\psi }_2\left(z\mathrm{,}\lambda \right)\frac{\partial }{\partial z}{\psi }_1\left(z\mathrm{,}\lambda \right)\right)=0.$ (2.17)

The same is true if the second solution is the complex-conjugate solution of the first solution for real parameter $\lambda$ that means (we add then a factor “I” to get a real-valued conservation quantity)

$\text{i}\frac{\partial }{\partial z}\left(\psi \left(z\mathrm{,}\lambda \right)\frac{\partial }{\partial z}{\psi }^{\mathrm{<mo>\; *\; </mo>}}\left(z\mathrm{,}\lambda \right)-{\psi }^{\mathrm{<mo>\; *\; </mo>}}\left(z\mathrm{,}\lambda \right)\frac{\partial }{\partial z}\psi \left(z\mathrm{,}\lambda \right)\right)=\mathrm{0,}\left(\lambda ={\lambda }^{\mathrm{<mo>\; *\; </mo>}}\right)\mathrm{.}$ (2.18)

Other general properties we find, for example, in [15] [16] [17] [18].

3. Sturm-Liouville Equation for Transition Potentials between Two Asymptotic Potentials and Its Solution by Boundary Conditions

We derive in this Section at an example generalized boundary conditions using the theory of Generalized functions (e.g., [25] [26] [27] [28], see also Appendix A). As starting point of the further treatment we take the following ordinary second-order differential Equation (2.16)

$\left(L-\lambda \right)\psi \left(z\mathrm{,}\lambda \right)=\mathrm{0,}L\equiv -\frac{{\partial }^2}{\partial z^2}+u\left(z\right)\mathrm{,}$ (3.1)

which in this form is commonly called the Sturm-Liouville equation. The operator L is the Sturm-Liouville operator and $\lambda$ its eigenvalues. The eigensolutions $\psi =\psi \left(z\mathrm{,}\lambda \right)$ depend on $\lambda$ as parameter.

We assume that the real-valued potential $u\left(z\right)$ is defined on the entire real z-axis ( $-\infty <z<+\infty$ ). The complex-valued solutions $\psi \left(z\mathrm{,}\lambda \right)$ are considered in dependence on the spectral parameter $\lambda$ which can take on arbitrary complex values. As a linear second-order differential equation, the Sturm-Liouville equation possesses two linearly independent solutions for each complex value of the spectral parameter $\lambda$. As mentioned the Sturm-Liouville Equation (3.1) can also be seen as the eigenvalue problem of the Sturm-Liouville operator L. Then the real eigenvalues $\lambda$ of the discrete and of the continuous spectrum may be separated from the totality of complex values of the spectral parameter $\lambda$ by appropriate normalization conditions of the solutions. The one-dimensional Schrödinger equation in energy representation is identical with (3.1) after simple linear substitutions of the variables (2.15), where the energy becomes proportional to $\lambda$. It is by no means a disadvantage to assume that $u\left(z\right)$ and $\psi \left(z\mathrm{,}\lambda \right)$ are Generalized functions and that the derivatives in (2.1) are also Generalized functions because this includes or generalizes all classical cases and the point-wise definition of functions is perhaps an unnecessary (and may be unnatural) concept in physics.

For the purpose of the illustration only the asymptotic potentials $u_{\mp }\left(z\right)$ for $z\to \pm \infty$ are chosen as constant potentials $u_{-}\left(z\right)=u\left(-\infty \right),u_{+}\left(z\right)=u\left(+\infty \right)$. The first partial figure shows the genuine potential and below this its approximation by a step function. The upper partial figure on the right-hand side shows in hatched way the error which we make by the approximation and below this the difference of potential and its approximation. The derivations show that the parts to the left and to the right of the chosen boundary are, in general, separately to deal with.

We assume now that the potential $u\left(z\right)$ can be approximated for $z\to -\infty$ by $u_{-}\left(z\right)$ and for $z\to +\infty$ by $u_{+}\left(z\right)$ and that the totality of solutions of the following two Sturm-Liouville equations

$0=\left\{-\frac{{\partial }^2}{\partial z^2}+u_{\mp }\left(z\right)-\lambda \right\}{\psi }_{\mp }\left(z\mathrm{,}\lambda \right)\mathrm{,}$ (3.2)

is known. The asymptotic potentials $u_{-}\left(z\right)$ and $u_{+}\left(z\right)$ have to be of such kind that the differences $u\left(z\right)-u_{-}\left(z\right)$ and $u\left(z\right)-u_{+}\left(z\right)$ become rapidly decreasing functions for $z\to +\infty$ or $z\to -\infty$, respectively. The rapid decrease of functions in infinity means that they decrease in infinity more rapidly than any power function ${\left|z\right|}^n$ with arbitrary negative numbers n or, more generally expressed, that they can be defined as linear continuous functionals over the space of smooth and moderately increasing functions, especially, the polynomials. The effective lengths of the transition regions should be small compared with the characteristic lengths of changes in the solutions for the asymptotic potentials $u_{+}\left(z\right)$ and $u_{-}\left(z\right)$ (e.g., wavelengths). Usually, $u_{+}\left(z\right)$ and $u_{-}\left(z\right)$ are chosen as such potentials for which we know the totality of solutions in form of elementary or well-known higher transcendental functions as, for instance, for constant potential (Exponential or Trigonometric functions), linear potential (Airy functions) and quadratic potential (Hermite functions or more generally Parabolic-cylinder functions).

An exact boundary between the two regions of $z\to -\infty$ and $z\to +\infty$ is not defined and for approximations of the derived boundary conditions a “good” choice of the reference point z₀ is necessary. A “good” choice in our case is $z_0=0$ and is drawn in Figure 1. The partial Figures illustrate the influence of the reference point z₀ onto the hatched correction functions.

There are also some other well-known special potentials with known exact solutions of the Sturm-Liouville equation as, for instance, the reflection-less Bargmann potentials with their importance in soliton theory of the Korteweg-de Vries equation (the famous Hyperbolic Secans square as the simplest nontrivial one). The solutions for the most common potentials can be found, e.g., in Vol. III (Quantum mechanics) of Landau and Lifshitz [20] (chap. 3) and furthermore in [21] (chap. 12.3) and in [22]. The derivation of relations between the asymptotic solutions ${\psi }_{+}\left(z\mathrm{,}\lambda \right)$ and ${\psi }_{-}\left(z\mathrm{,}\lambda \right)$ belongs to reflection and refraction or scattering problems of the Sturm-Liouville equation. By means of the generalized boundary conditions which we will derive now, one can obtain approximate solutions of the direct scattering problem of Equation (3.1).

The almost smooth transition of the potential $u\left(z\right)$ from the asymptotic potential $u_{-}\left(z\right)$ for $z\to -\infty$ to the asymptotic potential $u_{+}\left(z\right)$ for $z\to +\infty$ has to be concentrated mainly on the neighborhood of a value z₀ of the variablez (Figure 1) which is not exactly determined and can be varied in certain way. We have then to choose a suitable point

$z=z_0,$ (3.3)

as the boundary or the reference point of the expansions in the proposed method.

Figure 1. Transition between two constant potentials at minus and plus infinity.

Figure 2. Transition between two constant potentials at infinity with variation of the reference point z₀ of the boundary.

A small variation of the reference point z₀ should not destroy the whole approximation (Figure 2). The special case of a potential barrier (Figure 3) or a potential well (by taking instead negative values of that for the barrier) (Figure 4) we consider more detailed in Sections 6 and 7.

Next we represent the solution $\psi \left(z\mathrm{,}\lambda \right)$ of Equation (2.16) after separation of

Figure 3. Potential barrier over constant potential at infinity with reference point at z₀.

Figure 4. Potential well under constant potential at infinity with reference point at z₀.

the asymptotic solutions ${\psi }_{\mp }\left(z\mathrm{,}\lambda \right)$ to the asymptotic potentials $u_{\mp }\left(z\right)$ by means of Heaviside’s step function $\theta \left(z\right)$ in the form ( $\theta \left(z\right)=0,\left(z<0\right),\theta \left(z\right)=1,\left(z>0\right)$ )

$\psi \left(z,\lambda \right)={\psi }_{-}\left(z,\lambda \right)\theta \left(z_0-z\right)+{\psi }_{+}\left(z,\lambda \right)\theta \left(z-z_0\right)+{\psi }^{\prime }\left(z_0;z,\lambda \right),$ (3.4)

where ${\psi }^{\prime }\left(z_0\mathrm{<mo>\; ;\; </mo>}z\mathrm{,}\lambda \right)$ depending in addition on the reference point z₀ is a small deviation from the sum of the two asymptotic solutions in the transition region and turns out to be a function of z rapidly decreasing for $z\to \pm \infty$ and mainly concentrated on the neighborhood of the point z₀. Below we make a moment-series expansion of ${\psi }^{\prime }\left(z_0\mathrm{<mo>\; ;\; </mo>}z\mathrm{,}\lambda \right)$ into a delta and its derivatives (see Appendix A).

Inserting (3.4) into the Sturm-Liouville Equation (2.16) the following equation results

$0=\left\{-\frac{{\partial }^2}{\partial z^2}+u\left(z\right)-\lambda \right\}\left({\psi }_{-}\left(z,\lambda \right)\theta \left(z_0-z\right)+{\psi }_{+}\left(z,\lambda \right)\theta \left(z-z_0\right)+{\psi }^{\prime }\left(z_0;z,\lambda \right)\right),$ (3.5)

which we have to solve. For this purpose we apply Generalized functions and an expansion of rapidly in infinity decreasing functions into delta function and its derivatives called moment series which we explain in Appendix A.

Using the Leibniz rule for the n-fold differentiation of a product of functions $f\left(z\right)g(\; z\; )$

$\frac{{\partial }^n}{\partial z^n}f\left(z\right)g\left(z\right)={\displaystyle \underset{k=0}{\overset{n}{\sum }}}\frac{n!}{\left(n-k\right)!k!}{f}^{\left(n-k\right)}\left(z\right)g^{\left(k\right)}\left(z\right),$ (3.6)

we can make the following identical transformation of a part of the right-hand side of (3.5) which brings the Heaviside step functions in front of the derivatives

$\begin{array}{l}-\frac{{\partial }^2}{\partial z^2}\left({\psi }_{-}\left(z,\lambda \right)\theta \left(z_0-z\right)+{\psi }_{+}\left(z,\lambda \right)\theta \left(z-z_0\right)\right)\\ =-\theta \left(z_0-z\right)\frac{{\partial }^2}{\partial z^2}{\psi }_{-}\left(z,\lambda \right)-\theta \left(z-z_0\right)\frac{{\partial }^2}{\partial z^2}{\psi }_{+}\left(z,\lambda \right)\\ +2\delta \left(z-z_0\right)\left(\frac{\partial }{\partial z}{\psi }_{-}\left(z,\lambda \right)-\frac{\partial }{\partial z}{\psi }_{+}\left(z,\lambda \right)\right)\\ +{\delta }^{\left(1\right)}\left(z-z_0\right)\left({\psi }_{-}\left(z,\lambda \right)-{\psi }_{+}\left(z,\lambda \right)\right)\end{array}$

$\begin{array}{l}=-\theta \left(z_0-z\right)\frac{{\partial }^2}{\partial z^2}{\psi }_{-}\left(z,\lambda \right)-\theta \left(z-z_0\right)\frac{{\partial }^2}{\partial z^2}{\psi }_{+}\left(z,\lambda \right)\\ +\delta \left(z-z_0\right)\left({\psi }_{-}^{\left(1\right)}\left(z_0,\lambda \right)-{\psi }_{+}^{\left(1\right)}\left(z_0,\lambda \right)\right)\\ +{\delta }^{\left(1\right)}\left(z-z_0\right)\left({\psi }_{-}\left(z_0,\lambda \right)-{\psi }_{+}\left(z_0,\lambda \right)\right).\end{array}$ (3.7)

We have applied the identity (A.5) to make the coefficients of the delta function and its first derivative independent on variable z. Equation (3.5) using this identity can be written

$\begin{array}{c}0=\theta \left(z_0-z\right)\left\{-\frac{{\partial }^2}{\partial z^2}+u\left(z\right)-\lambda \right\}{\psi }_{-}\left(z,\lambda \right)\\ +\theta \left(z-z_0\right)\left\{-\frac{{\partial }^2}{\partial z^2}+u\left(z\right)-\lambda \right\}{\psi }_{+}\left(z,\lambda \right)\end{array}$

$\begin{array}{c}+\delta \left(z-z_0\right)\left({\psi }_{-}^{\left(1\right)}\left(z_0,\lambda \right)-{\psi }_{+}^{\left(1\right)}\left(z_0,\lambda \right)\right)\\ +{\delta }^{\left(1\right)}\left(z-z_0\right)\left({\psi }_{-}\left(z_0,\lambda \right)-{\psi }_{+}\left(z_0,\lambda \right)\right)\\ +\left\{-\frac{{\partial }^2}{\partial z^2}+u\left(z\right)-\lambda \right\}{\psi }^{\prime }\left(z_0;z,\lambda \right).\end{array}$ (3.8)

With the asymptotic Equations (3.2) we may transform this equation to

$\begin{array}{c}0=+\delta \left(z-z_0\right)\left({\psi }_{-}^{\left(1\right)}\left(z_0\mathrm{,}\lambda \right)-{\psi }_{+}^{\left(1\right)}\left(z_0\mathrm{,}\lambda \right)\right)\\ +{\delta }^{\left(1\right)}\left(z-z_0\right)\left({\psi }_{-}\left(z_0\mathrm{,}\lambda \right)-{\psi }_{+}\left(z_0\mathrm{,}\lambda \right)\right)\\ +\theta \left(z_0-z\right)\left(u\left(z\right)-u_{-}\left(z\right)\right){\psi }_{-}\left(z\mathrm{,}\lambda \right)\\ +\theta \left(z-z_0\right)\left(u\left(z\right)-u_{+}\left(z\right)\right){\psi }_{+}\left(z\mathrm{,}\lambda \right)\\ +\left\{-\frac{{\partial }^2}{\partial z^2}+u\left(z\right)-\lambda \right\}{\psi }^{\prime }\left(z_0\mathrm{<mo>\; ;\; </mo>}z\mathrm{,}\lambda \right)\mathrm{.}\end{array}$ (3.9)

Our aim is to represent this equation as a superposition of the delta function and of its derivatives ${\delta }^{\left(n\right)}\left(z-z_0\right),\left(n=0,1,2,\cdots \right)$ with respect to variable z where the coefficients should not depend on variable z. For the sum term with the function ${\psi }^{\prime }\left(z_0\mathrm{<mo>\; ;\; </mo>}z\mathrm{,}\lambda \right)$ we make this below and consider now the two sum terms in the second line of (3.9).

The deviation $\Delta u\left(z_0\mathrm{<mo>\; ;\; </mo>}z\right)$ of the potential $u\left(z\right)$ in the transition region from the two asymptotic potentials $u_{-}\left(z\right)$ and $u_{+}\left(z\right)$ depends on the chosen boundary z₀ and will be defined by

$\begin{array}{c}\Delta u\left(z_0\mathrm{<mo>\; ;\; </mo>}z\right)\equiv u\left(z\right)-u_{-}\left(z\right)\theta \left(z_0-z\right)-u_{+}\left(z\right)\theta \left(z-z_0\right)\\ =\left(u\left(z\right)-u_{-}\left(z\right)\right)\theta \left(z_0-z\right)+\left(u\left(z\right)-u_{+}\left(z\right)\right)\theta \left(z-z_0\right)\mathrm{.}\end{array}$ (3.10)

However, since it does not find entrance in Equation (3.9) as a whole but is combined with the solutions ${\psi }_{-}\left(z\mathrm{,}\lambda \right)$ and ${\psi }_{+}\left(z\mathrm{,}\lambda \right)$ we have to deal with the two parts $u\left(z\right)-u_{-}\left(z\right)$ for $z\le z_0$ and $u\left(z\right)-u_{+}\left(z\right)$ for $z\ge z_0$ separately.

Since $u\left(z\right)-u_{-}\left(z\right)$ and $u\left(z\right)-u_{+}\left(z\right)$ by presupposition are concentrated in the neighborhood of the chosen boundary $z=z_0$ and are rapidly vanishing for $z\to \mp \infty$ one may make moment-series expansions of the form (see Equations (A.23) and (A.24) in Appendix A)

$\begin{array}{l}\left(u\left(z\right)-u_{-}\left(z\right)\right)\theta \left(z_0-z\right)={\displaystyle \underset{m=0}{\overset{\infty }{\sum }}}{u}_{-,m}\left(z_0\right){\delta }^{\left(m\right)}\left(z-z_0\right),\\ \left(u\left(z\right)-u_{+}\left(z\right)\right)\theta \left(z-z_0\right)={\displaystyle \underset{m=0}{\overset{\infty }{\sum }}}{u}_{+,m}\left(z_0\right){\delta }^{\left(m\right)}\left(z-z_0\right),\end{array}$ (3.11)

with coefficients $u_{-\mathrm{,}m}\left(z_0\right)$ and $u_{+\mathrm{,}m}\left(z_0\right)$ which are left-hand and right-hand moments of the hatched error parts of the potential which do not depend on variable z (but on reference point z₀) and are given by

$\begin{array}{l}{u}_{-,m}\left(z_0\right)\equiv \frac{{\left(-1\right)}^m}{m!}{\displaystyle {\int }_{-\infty }^{z_0}}\text{d}z\left(u\left(z\right)-u_{-}\left(z\right)\right){\left(z-z_0\right)}^m\\ =\frac{1}{m!}{\displaystyle {\int }_{-\infty }^{z_0}}\text{d}z\left(u\left(z\right)-u_{-}\left(z\right)\right)\underset{\ge 0}{\underset{︸}{{\left(z_0-z\right)}^m}},\end{array}$

$u_{+,m}\left(z_0\right)\equiv \frac{{\left(-1\right)}^m}{m!}{\displaystyle {\int }_{z_0}^{+\infty }}\text{d}z\left(u\left(z\right)-u_{+}\left(z\right)\right)\underset{\ge 0}{\underset{︸}{{\left(z-z_0\right)}^m}}.$ (3.12)

These are integrals over the hatched parts in last partial figure in Figure 2 multiplied by ${\left(z-z_0\right)}^m$.

We have now to multiply the expressions in (3.11) with ${\psi }_{-}\left(z\mathrm{,}\lambda \right)$ and ${\psi }_{+}\left(z\mathrm{,}\lambda \right)$, respectively, and using the multiplication formula (A.5) after reordering of the arising double sums we obtain

$\begin{array}{l}\left(u\left(z\right)-u_{-}\left(z\right)\right)\theta \left(z_0-z\right){\psi }_{-}\left(z,\lambda \right)+\left(u\left(z\right)-u_{+}\left(z\right)\right)\theta \left(z-z_0\right){\psi }_{+}\left(z,\lambda \right)\\ ={\displaystyle \underset{m=0}{\overset{\infty }{\sum }}}\left(u_{-,m}\left(z_0\right){\delta }^{\left(m\right)}\left(z-z_0\right){\psi }_{-}\left(z,\lambda \right)+u_{+,m}\left(z_0\right){\delta }^{\left(m\right)}\left(z-z_0\right){\psi }_{+}\left(z,\lambda \right)\right)\\ ={\displaystyle \underset{m=0}{\overset{\infty }{\sum }}}{\displaystyle \underset{k=0}{\overset{m}{\sum }}}\frac{{\left(-1\right)}^{k}m!}{k!\left(m-k\right)!}\left(u_{-,m}\left(z_0,\lambda \right){\psi }_{-}^{\left(k\right)}\left(z_0,\lambda \right)+u_{+,m}\left(z_0,\lambda \right){\psi }_{+}^{\left(k\right)}\left(z_0,\lambda \right)\right){\delta }^{\left(m-k\right)}\left(z-z_0\right)\\ ={\displaystyle \underset{n=0}{\overset{\infty }{\sum }}}{\displaystyle \underset{k=0}{\overset{\infty }{\sum }}}\frac{{\left(-1\right)}^k\left(n+k\right)!}{n!k!}\left(u_{-,n+k}\left(z_0\right){\psi }_{-}^{\left(k\right)}\left(z_0,\lambda \right)+u_{+,n+k}\left(z_0\right){\psi }_{+}^{\left(k\right)}\left(z_0,\lambda \right)\right){\delta }^{\left(n\right)}\left(z-z_0\right),\end{array}$ (3.13)

The hatched error function to the approximation does not depend on the reference point z₀ but the moments depend on it. However, the first moment ${\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}z\left(u\left(z\right)-u\left(\infty \right)\right)$ is independent of the chosen reference point z₀. Similar is the case of a potential well between two asymptotic equal potentials by changing the sign of $u\left(z\right)$.

The hatched error function to the approximation does not depend on the reference point z₀ but the moments depend on it. However, the first moment ${\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}z\left(u\left(z\right)-u\left(\infty \right)\right)$ is independent of the chosen reference point z₀.

There remains to represent the additional term ${\psi }^{\prime }\left(z_0\mathrm{<mo>\; ;\; </mo>}z\mathrm{,}\lambda \right)$ to the asymptotic solution in (3.5) which is concentrated onto the neighborhood of $z=z_0$ as a superposition of the delta function and of their derivatives that means the expansion into a moment series according to (see Appendix A)

${\psi }^{\prime }\left(z_0;z,\lambda \right)={\displaystyle \underset{m=0}{\overset{\infty }{\sum }}}{{\psi }^{\prime }}_m\left(z_0,\lambda \right){\delta }^{\left(m\right)}\left(z-z_0\right),$ (3.14)

The hatched error function to the approximation does not depend on the reference point z₀ but the moments dependon it. However, the first moment ${\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}z\left(u\left(z\right)-u\left(\infty \right)\right)$ is independent of the chosen reference point z₀.

With the coefficients ${{\psi }^{\prime }}_m\left(z_0\mathrm{,}\lambda \right)$ which do not depend on variable z but on the reference point z₀

${{\psi }^{\prime }}_m\left(z_0\mathrm{,}\lambda \right)\equiv \frac{{\left(-1\right)}^m}{m\mathrm{<mo>\; !\; </mo>}}{\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}z{\psi }^{\prime }\left(z_0\mathrm{<mo>\; ;\; </mo>}z\mathrm{,}\lambda \right){\left(z-z_0\right)}^m\mathrm{.}$ (3.15)

Now we have to apply the operator $L-\lambda$ of the Sturm-Liouville equation to the moment-series expansion of the “small” additional part ${\psi }^{\prime }\left(z_0\mathrm{<mo>\; ;\; </mo>}z\mathrm{,}\lambda \right)$ of the solution. Using the multiplication theorem (A.5) and the expansion (3.14) we find

$-\frac{{\partial }^2}{\partial z^2}{\psi }^{\prime }\left(z_0;z,\lambda \right)=-{\displaystyle \underset{n=2}{\overset{\infty }{\sum }}}{{\psi }^{\prime }}_{n-2}\left(z_0,\lambda \right){\delta }^{\left(n\right)}\left(z-z_0\right),$

$u\left(z\right){\psi }^{\prime }\left(z_0;z,\lambda \right)={\displaystyle \underset{n=0}{\overset{\infty }{\sum }}}\left({\displaystyle \underset{k=0}{\overset{\infty }{\sum }}}{\left(-1\right)}^k\frac{\left(n+k\right)!}{n!k!}{u}^{\left(k\right)}\left(z_0\right){{\psi }^{\prime }}_{n+k}\left(z_0,\lambda \right)\right){\delta }^{\left(n\right)}\left(z-z_0\right),$

$-\lambda {\psi }^{\prime }\left(z_0;z,\lambda \right)=-\lambda {\displaystyle \underset{n=0}{\overset{\infty }{\sum }}}{{\psi }^{\prime }}_n\left(z_0,\lambda \right){\delta }^{\left(n\right)}\left(z-z_0\right).$ (3.16)

Thus we have prepared the representation of all parts for insertion into Equation (3.9).

Inserting now (3.13) and (3.16) into Equation (3.9) and collecting all terms proportional to $\delta \left(z-z_0\right)$ and to its derivatives ${\delta }^{\left(n\right)}\left(z-z_0\right),\left(n=1,2,\cdots \right)$ we find

$\begin{array}{l}0=\{{\psi }_{-}^{\left(1\right)}\left(z_0,\lambda \right)-{\psi }_{+}^{\left(1\right)}\left(z_0,\lambda \right)+{\displaystyle \underset{k=0}{\overset{\infty }{\sum }}}{\left(-1\right)}^k(u_{-,k}\left(z_0\right){\psi }_{-}^{\left(k\right)}\left(z_0,\lambda \right)\\ \underset{}{\overset{}{}}+u_{+,k}\left(z_0\right){\psi }_{+}^{\left(k\right)}\left(z_0,\lambda \right)+u^{\left(k\right)}\left(z_0\right){{\psi }^{\prime }}_k\left(z_0,\lambda \right))-\lambda {{\psi }^{\prime }}_0\left(z_0,\lambda \right)\}\delta \left(z-z_0\right)\\ +\{{\psi }_{-}\left(z_0,\lambda \right)-{\psi }_{+}\left(z_0,\lambda \right)+{\displaystyle \underset{k=0}{\overset{\infty }{\sum }}}{\left(-1\right)}^k\frac{\left(1+k\right)!}{1!k!}(u_{-,1+k}\left(z_0\right){\psi }_{-}^{\left(k\right)}\left(z_0,\lambda \right)\\ \underset{}{\overset{}{}}+u_{+,1+k}\left(z_0\right){\psi }_{+}^{\left(k\right)}\left(z_0,\lambda \right)+u^{\left(k\right)}\left(z_0,\lambda \right){{\psi }^{\prime }}_{1+k}\left(z_0,\lambda \right))-\lambda {{\psi }^{\prime }}_1\left(z_0,\lambda \right)\}{\delta }^{\left(1\right)}\left(z-z_0\right)\\ +{\displaystyle \underset{n=2}{\overset{\infty }{\sum }}}\{{\displaystyle \underset{k=0}{\overset{\infty }{\sum }}}{\left(-1\right)}^k\frac{\left(n+k\right)!}{n!k!}(u_{-,n+k}\left(z_0\right){\psi }_{-}^{\left(k\right)}\left(z_0,\lambda \right)+u_{+,n+k}\left(z_0\right){\psi }_{+}^{\left(k\right)}\left(z_0,\lambda \right)\\ \underset{}{\overset{}{}}+u^{\left(k\right)}\left(z_0\right){{\psi }^{\prime }}_{n+k}\left(z_0,\lambda \right))-\left({{\psi }^{\prime }}_{n-2}\left(z_0,\lambda \right)+\lambda {{\psi }^{\prime }}_n\left(z_0,\lambda \right)\right)\}{\delta }^{\left(n\right)}\left(z-z_0\right).\end{array}$ (3.17)

All coefficients in front of the derivatives of the delta functions in this equation are now independent on the variable z and according to (3.12) and (3.13) they have to vanish since the left-hand side of (3.17) vanishes. This leads to the intended boundary conditions.

Summarizing, the whole solution of the Sturm-Liouville equation was taken in the general form

$\psi \left(z,\lambda \right)={\psi }_{-}\left(z,\lambda \right)\theta \left(z_0-z\right)+{\psi }_{+}\left(z,\lambda \right)\theta \left(z-z_0\right)+{\displaystyle \underset{n=0}{\overset{\infty }{\sum }}}{{\psi }^{\prime }}_n\left(z_0,\lambda \right){\delta }^{\left(n\right)}\left(z-z_0\right),$ (3.18)

and after application of the Sturm-Liouville operator to (3.18) it has to vanish and by separation of the asymptotic solutions (3.2) there remains a part which could be represented as a superposition of the delta function and its derivatives with constant coefficients (3.17) which do not depend on the variable z but only on the chosen boundary at $z=z_0$ and the parameter $\lambda$ and all the coefficients have to vanish. In this way we arrived at the intended boundary conditions. Since this system of boundary conditions is infinite we have to look for approximations. This is done in next Section.

4. The Infinite System of Boundary Conditions and Their Lowest Three Approximations

The boundary conditions which are as explained the vanishing of the (z-independent) coefficients in front of the delta function and of its derivatives in (3.17) are the following:

1) Vanishing of coefficient in front of $\delta \left(z-z_0\right)$

$\begin{array}{l}0={\psi }_{-}^{\left(1\right)}\left(z_0,\lambda \right)-{\psi }_{+}^{\left(1\right)}\left(z_0,\lambda \right)+{\displaystyle \underset{k=0}{\overset{\infty }{\sum }}}{\left(-1\right)}^k(u_{-,k}\left(z_0\right){\psi }_{-}^{\left(k\right)}\left(z_0,\lambda \right)\\ +u_{+,k}\left(z_0\right){\psi }_{+}^{\left(k\right)}\left(z_0,\lambda \right)+u^{\left(k\right)}\left(z_0\right){{\psi }^{\prime }}_k\left(z_0,\lambda \right))-\lambda {{\psi }^{\prime }}_0\left(z_0,\lambda \right),\end{array}$ (4.1)

2) Vanishing of coefficient in front of ${\delta }^{\left(1\right)}\left(z-z_0\right)$

$\begin{array}{l}0={\psi }_{-}\left(z_0,\lambda \right)-{\psi }_{+}\left(z_0,\lambda \right)+{\displaystyle \underset{k=0}{\overset{\infty }{\sum }}}{\left(-1\right)}^k\left(1+k\right)(u_{-,1+k}\left(z_0\right){\psi }_{-}^{\left(k\right)}\left(z_0,\lambda \right)\\ +u_{+,1+k}\left(z_0\right){\psi }_{+}^{\left(k\right)}\left(z_0,\lambda \right)+u^{\left(k\right)}\left(z_0\right){{\psi }^{\prime }}_{1+k}\left(z_0,\lambda \right))-\lambda {{\psi }^{\prime }}_1\left(z_0,\lambda \right),\end{array}$ (4.2)

3) Vanishing of coefficients in front of ${\delta }^{\left(n\right)}\left(z-z_0\right)\mathrm{,}\left(n=\mathrm{2,3,4,}\cdots \right)$

$\begin{array}{l}0={\displaystyle \underset{k=0}{\overset{\infty }{\sum }}}{\left(-1\right)}^k\frac{\left(n+k\right)!}{n!k!}(u_{-,n+k}\left(z_0\right){\psi }_{-}^{\left(k\right)}\left(z_0,\lambda \right)+u_{+,n+k}\left(z_0\right){\psi }_{+}^{\left(k\right)}\left(z_0,\lambda \right)\\ +u^{\left(k\right)}\left(z_0\right){{\psi }^{\prime }}_{n+k}\left(z_0,\lambda \right))-\left({{\psi }^{\prime }}_{n-2}\left(z_0,\lambda \right)+\lambda {{\psi }^{\prime }}_n\left(z_0,\lambda \right)\right).\end{array}$ (4.3)

Apart from the given asymptotic potentials $u_{\mp }\left(z\right)$ the properties of the complete potential function $u\left(z\right)$ are contained in the moments $u_{\mp \mathrm{,}m}\left(z_0\right)$ (lower indices $m=0,1,2,\cdots$ ) defined in (3.12) and in the k-th derivatives $u^{\left(k\right)}\left(z_0\right)$ of $u\left(z\right)$ at reference point $z=z_0$ (upper indices $\left(k=0,1,2,\cdots \right)$ ) and their definitions are given in last Section, where functions with lower indices mean moments and with upper indices derivatives of the corresponding functions.

All approximations of the boundary condition should be contained in (4.1), (4.2) and (4.3) but truncations for approximations have to be made in non-contradictory way. A difficulty is that in each of the conditions an infinite number of potential properties (summations over k) is involved and therefore that the truncation is not a simple problem. The properties of the transition region between the two asymptotic parts is described after the choice of a reference point z₀ by the moments $u_{\mp \mathrm{,}n}\left(z_0\right)$ defined by (3.12) and in connection with (3.11). As N-th approximation we will denote the case where all moments $n\ge N$ are neglected that means

$u_{\mp ,n}\left(z_0\right)=0,\left(n\ge N\right).$ (4.4)

It seems that an appropriate strategy to obtain finite approximation of boundary condition is to truncate the expansion of the solution (3.18) from a finite $n\ge N$ on as follows

${\psi }_n\left(z_0,\lambda \right)=0,\left(n\ge N\right).$ (4.5)

Then the approximate solution possesses the form

$\psi \left(z,\lambda \right)={\psi }_{-}\left(z,\lambda \right)\theta \left(z_0-z\right)+{\psi }_{+}\left(z,\lambda \right)\theta \left(z-z_0\right)+{\displaystyle \underset{n=0}{\overset{N-1}{\sum }}}{{\psi }^{\prime }}_n\left(z_0,\lambda \right){\delta }^{\left(n\right)}\left(z-z_0\right).$ (4.6)

We call both conditions (4.4) and (4.5) taken together the N-th approximation of the boundary conditions. The boundary conditions (4.3) are then identically satisfied in N-th approximation for $n\ge N+2$ and have to be not taken into account for $n\ge N$. Thus we have $N+2$ boundary conditions for $N+4$ quantities ${\psi }_{\mp }^{\left(1\right)}\left(z_0\mathrm{,}\lambda \right)$ and ${\psi }_{\mp }\left(z_0\mathrm{,}\lambda \right)$ (4 quantities) plus ${{\psi }^{\prime }}_0\left(z_0\mathrm{,}\lambda \right)\mathrm{,}\cdots \mathrm{,}{{\psi }^{\prime }}_{N-1}\left(z_0\mathrm{,}\lambda \right)$ (N quantities) in the N-th approximation. For $N\to \infty$ all terms in the boundary conditions are taken into account. In a sufficiently high approximation one may neglect additional terms if they prove to be negligible. We now consider explicitly the three lowest approximations to the basic boundary conditions (or zeroth approximation).

4.1. Zeroth-Order or Basic Approximation of the Boundary Conditions

In the zeroth-order (or basic) approximation we consider the solution in the approximate form

$\psi \left(z,\lambda \right)={\psi }_{-}\left(z,\lambda \right)\theta \left(z_0-z\right)+{\psi }_{+}\left(z,\lambda \right)\theta \left(z-z_0\right),$ (4.7)

and neglect all moments, i.e. $u_{\mp ,n}\left(z_0\right)=0,\left(n\ge 0\right)$ for the transition region. This leads according to (4.1) and (4.2) to the basic boundary conditions

$0={\psi }_{-}^{\left(1\right)}\left(z_0,\lambda \right)-{\psi }_{+}^{\left(1\right)}\left(z_0,\lambda \right),$

$0={\psi }_{-}\left(z_0,\lambda \right)-{\psi }_{+}\left(z_0,\lambda \right).$ (4.8)

These are the well-known conditions of continuity of the asymptotic solutions ${\psi }_{\mp }\left(z\mathrm{,}\lambda \right)$ and of their first derivatives ${\psi }_{\mp }^{\left(1\right)}\left(z\mathrm{,}\lambda \right)$ at the chosen boundary $z=z_0$. They correspond to the approximation of the genuine potential $u\left(z\right)$ by the abrupt transition of $u_{-}\left(z\right)$ into the potential $u_{+}\left(z\right)$ at the reference point $z=z_0$ (see Figure 1 and Figure 2).

The Sturm-Liouville equation provides for arbitrary fixed parameter $\lambda$ and asymptotic potentials $u_{-}\left(z\right)$ and $u_{+}\left(z\right)$ solutions ${\psi }_{-}\left(z\mathrm{,}\lambda \right)$ and ${\psi }_{+}\left(z\mathrm{,}\lambda \right)$, respectively, with the superposition of two possible basic solutions with 2 free amplitudes in each case and this total number of 4 amplitudes is related by the 2 boundary conditions to 2 free amplitudes (e.g., incident waves from each side of the boundary) and to 2 determinable amplitudes from them (e.g., for reflected or refracted waves in dependence from which side one sees this).

We now consider finite approximations to the boundary conditions (4.8) taking into account a transition region between the two asymptotic potentials with given properties. They have to be non-contradictory and each higher approximation should contain the preceding approximation by neglect of terms. In this way we can assume that in an infinite process which, clearly, is impossible to make every term in the boundary conditions will be taken into account.

4.2. First-Order Approximation of the Boundary Conditions

In the first-order approximation of the boundary conditions we consider the solution in the approximate form

$\psi \left(z,\lambda \right)={\psi }_{-}\left(z,\lambda \right)\theta \left(z_0-z\right)+{\psi }_{+}\left(z,\lambda \right)\theta \left(z-z_0\right)+{{\psi }^{\prime }}_0\left(z_0,\lambda \right)\delta \left(z-z_0\right),$ (4.9)

and in addition we take into account the moments $u_{\mp \mathrm{,0}}\left(z_0\right)$ and neglect all higher moments, i.e. $u_{\mp ,n}\left(z_0\right)=0,\left(n\ge 1\right)$.

Then we find from (4.1), (4.2) and (4.3) for $n=2$ the following boundary conditions

$\begin{array}{l}0={\psi }_{-}^{\left(1\right)}\left(z_0,\lambda \right)-{\psi }_{+}^{\left(1\right)}\left(z_0,\lambda \right)+(u_{-\mathrm{,0}}\left(z_0\right){\psi }_{-}\left(z_0,\lambda \right)\\ +u_{+\mathrm{,0}}\left(z_0\right){\psi }_{+}\left(z_0,\lambda \right)+\left(u\left(z_0\right)-\lambda \right){{\psi }^{\prime }}_0\left(z_0,\lambda \right))\end{array}$

$0={\psi }_{-}\left(z_0\mathrm{,}\lambda \right)-{\psi }_{+}\left(z_0\mathrm{,}\lambda \right)\mathrm{,}$

$0=-{{\psi }^{\prime }}_0\left(z_0\mathrm{,}\lambda \right)\mathrm{.}$ (4.10)

Due to ${\psi }_0\left(z_0\mathrm{,}\lambda \right)=0$ one obtains by its elimination from the first equation and using the equality ${\psi }_{-}\left(z_0\mathrm{,}\lambda \right)={\psi }_{+}\left(z_0\mathrm{,}\lambda \right)$ effectively the following two boundary conditions in first approximation

$0={\psi }_{-}^{\left(1\right)}\left(z_0\mathrm{,}\lambda \right)-{\psi }_{+}^{\left(1\right)}\left(z_0\mathrm{,}\lambda \right)+\frac{u_0\left(z_0\right)}{2}\left({\psi }_{-}\left(z_0\mathrm{,}\lambda \right)+{\psi }_{+}\left(z_0\mathrm{,}\lambda \right)\right)\mathrm{,}$

$0={\psi }_{-}\left(z_0\mathrm{,}\lambda \right)-{\psi }_{+}\left(z_0\mathrm{,}\lambda \right)\mathrm{.}$ (4.11)

In considered approximation the moments $u_{-\mathrm{,0}}\left(z_0\right)$ and $u_{+\mathrm{,0}}\left(z_0\right)$ are involved only as their sum $u_{-\mathrm{,0}}\left(z_0\right)+u_{+\mathrm{,0}}\left(z_0\right)$ (see also Figure 1 and Figure 2)

$\begin{array}{c}{u}_0\left(z_0\right)=u_{-\mathrm{,0}}\left(z_0\right)+u_{+\mathrm{,0}}\left(z_0\right)\\ ={\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}z\left(u\left(z\right)-u_{-}\left(z\right)\theta \left(z_0-z\right)-u_{+}\left(z\right)\theta \left(z-z_0\right)\right)\\ \equiv {\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}z\Delta u\left(z_0\mathrm{<mo>\; ;\; </mo>}z\right)\mathrm{.}\end{array}$ (4.12)

A quantity such as $u_{-\mathrm{,0}}\left(z_0\right)+u_{+\mathrm{,0}}\left(z_0\right)$ and its dependence on the choice of z₀ in optics and its connection to higher such moments was discussed by Lekner in his versatile monograph [13] in chap. 3 (in particular p. 81, Figure 3.1 and Eq. (3.27)) and they were called integral invariants but its origin from moment-series expansion of function (see Equations (A.17) and (A.23) in Appendix A) is not mentioned there. Such quantities played already a great role in our paper [12]. It became clear there that the functions $\left(u_{-}\left(z\right)-u\left(z\right)\right)\theta \left(z_0-z\right)$ and $\left(u_{+}\left(z\right)-u\left(z\right)\right)\theta \left(z-z_0\right)$ (correspondingly $\epsilon \left(z\mathrm{<mo>\; ;\; </mo>}\omega \right)$ in optics) and their moments have to be considered separately and that they usually are involved in higher-order approximations separately (see next subsections). In the old papers of Drude [1] [2] (see also [3] [5] [6] ) moments of the permeability play also a role but the asymptotic parts are not split there and the theory of Generalized functions they did not have to disposal.

If the asymptotic potentials $u_{-}\left(z\right)$ and $u_{+}\left(z\right)$ are smooth (existence of all its derivatives) then by Taylor-series expansions of the step functions $\theta \left(z_0-z\right)$ and $\theta \left(z-z_0\right)$ in powers of z₀ and using formulae (A.5) and (A.2) one obtains the following expansion of $u_0\left(z_0\right)$ in powers of z₀

$\begin{array}{c}{u}_0\left(z_0\right)=u_0\left(0\right)-{\displaystyle \underset{k=0}{\overset{\infty }{\sum }}}\left(u_{-}^{\left(k\right)}\left(0\right)-u_{+}^{\left(k\right)}\left(0\right)\right)\frac{z_0^{k+1}}{\left(k+1\right)!}\\ =u_0\left(0\right)-\left(u_{-}\left(0\right)-u_{+}\left(0\right)\right)z_0-\left(u_{-}^{\left(1\right)}\left(0\right)-u_{+}^{\left(1\right)}\left(0\right)\right)\frac{z_0^2}{2}-\cdots .\end{array}$ (4.13)

A favorable choice of the reference point z₀ is often the condition $u_0\left(z_0\right)=0$, approximately

$0\approx u_0\left(0\right)-\left(u_{-}\left(0\right)-u_{+}\left(0\right)\right)z_0\mathrm{,}\Rightarrow z_0\approx \frac{u_0\left(0\right)}{u_{-}\left(0\right)-u_{+}\left(0\right)}\mathrm{,}$ (4.14)

if $u_{-}\left(0\right)-u_{+}\left(0\right)$ is non-vanishing. Then $u_0\left(z_0\right)\equiv u_{-\mathrm{,0}}\left(z_0\right)+u_{+\mathrm{,0}}\left(z_0\right)$ is nearly vanishing and this choice is unique but, generally, it can be not unique due to nonlinearity of (4.13) in z₀.

The only effect in this approximation in comparison to zeroth approximation is then an influence on the phase relations between the asymptotic solutions on both sides. In case of not monotonically increasing or decreasing transitions the choice of vanishing moment $u_{-\mathrm{,0}}\left(z_0\right)+u_{+\mathrm{,0}}\left(z_0\right)$ is not in every case possible and unique, in particular, in case of a potential barrier or well between two equal asymptotic potentials we cannot suppose a monotonic increase or decrease of transition between the two potentials and this case has to be dealt with separately (see Section 5).

The unique choice of a boundary z₀ by the condition (4.14) can be also applied under the discussed suppositions in the higher-order approximations with favorable influence on the smallness of the involved higher moments and we do not discuss this again in the next considered approximations.

4.3. Second-Order Approximation of the Boundary Conditions

In the second-order approximation of the boundary conditions we use for the solution the following proposition

$\begin{array}{l}\psi \left(z\mathrm{,}\lambda \right)={\psi }_{-}\left(z\mathrm{,}\lambda \right)\theta \left(z_0-z\right)+{\psi }_{+}\left(z\mathrm{,}\lambda \right)\theta \left(z-z_0\right)\\ +{{\psi }^{\prime }}_0\left(z_0\mathrm{,}\lambda \right)\delta \left(z-z_0\right)+{{\psi }^{\prime }}_1\left(z_0\mathrm{,}\lambda \right){\delta }^{\left(1\right)}\left(z-z_0\right)\mathrm{,}\end{array}$ (4.15)

and take into account as non-vanishing parameters the moments $u_{\mp \mathrm{,0}}\left(z_0\right)$ and $u_{\mp \mathrm{,1}}\left(z_0\right)$ and neglect all others, i.e. $u_{\mp ,n}\left(z_0\right)=0,\left(n\ge 2\right)$. From (4.1) and (4.2) and from (4.3) for $n=2,3$ we find then

$\begin{array}{l}0={\psi }_{-}^{\left(1\right)}\left(z_0\mathrm{,}\lambda \right)-{\psi }_{+}^{\left(1\right)}\left(z_0\mathrm{,}\lambda \right)\\ +\left(u_{-\mathrm{,0}}\left(z_0\right){\psi }_{-}\left(z_0\mathrm{,}\lambda \right)+u_{+\mathrm{,0}}\left(z_0\right){\psi }_{+}\left(z_0\mathrm{,}\lambda \right)+\left(u\left(z_0\right)-\lambda \right){{\psi }^{\prime }}_0\left(z_0\mathrm{,}\lambda \right)\right)\\ -\left(u_{-\mathrm{,1}}\left(z_0\right){\psi }_{-}^{\left(1\right)}\left(z_0\mathrm{,}\lambda \right)+u_{+\mathrm{,1}}\left(z_0\right){\psi }_{+}^{\left(1\right)}\left(z_0\mathrm{,}\lambda \right)+u^{\left(1\right)}\left(z_0\right){{\psi }^{\prime }}_1\left(z_0\mathrm{,}\lambda \right)\right)\mathrm{,}\end{array}$

$\begin{array}{l}0={\psi }_{-}\left(z_0\mathrm{,}\lambda \right)-{\psi }_{+}\left(z_0\mathrm{,}\lambda \right)+(u_{-\mathrm{,1}}\left(z_0\right){\psi }_{-}\left(z_0\mathrm{,}\lambda \right)\\ +u_{+\mathrm{,1}}\left(z_0\right){\psi }_{+}\left(z_0\mathrm{,}\lambda \right)+\left(u\left(z_0\right)-\lambda \right){{\psi }^{\prime }}_1\left(z_0\mathrm{,}\lambda \right))\mathrm{,}\end{array}$

$0=-{{\psi }^{\prime }}_0\left(z_0\mathrm{,}\lambda \right)\mathrm{,}$

$0=-{{\psi }^{\prime }}_1\left(z_0\mathrm{,}\lambda \right)\mathrm{.}$ (4.16)

Using the last two conditions of vanishing of ${{\psi }^{\prime }}_0\left(z_0\mathrm{,}\lambda \right)$ and ${{\psi }^{\prime }}_1\left(z_0\mathrm{,}\lambda \right)$ the boundary conditions in the second-order approximation can be reduced to the two conditions

$\begin{array}{l}0=\left(1-u_{-\mathrm{,1}}\left(z_0\right)\right){\psi }_{-}^{\left(1\right)}\left(z_0\mathrm{,}\lambda \right)-\left(1+u_{+\mathrm{,1}}\left(z_0\right)\right){\psi }_{+}^{\left(1\right)}\left(z_0\mathrm{,}\lambda \right)\\ +\left(u_{-\mathrm{,0}}\left(z_0\right){\psi }_{-}\left(z_0\mathrm{,}\lambda \right)+u_{+\mathrm{,0}}\left(z_0\right){\psi }_{+}\left(z_0\mathrm{,}\lambda \right)\right)\mathrm{,}\end{array}$

$0=\left(1+u_{-\mathrm{,1}}\left(z_0\right)\right){\psi }_{-}\left(z_0\mathrm{,}\lambda \right)-\left(1-u_{+\mathrm{,1}}\left(z_0\right)\right){\psi }_{+}\left(z_0\mathrm{,}\lambda \right)\mathrm{,}$ (4.17)

and a favorable choice of the boundary z₀ is mostly possible by the condition (4.14).

4.4. Third-Order Approximation of the Boundary Conditions

In the third-order approximation of the boundary conditions we use the following proposition for the solution

$\begin{array}{c}\psi \left(z,\lambda \right)={\psi }_{-}\left(z,\lambda \right)\theta \left(z_0-z\right)+{\psi }_{+}\left(z,\lambda \right)\theta \left(z-z_0\right)+{{\psi }^{\prime }}_0\left(z_0,\lambda \right)\delta \left(z-z_0\right)\\ +{{\psi }^{\prime }}_1\left(z_0,\lambda \right){\delta }^{\left(1\right)}\left(z-z_0\right)+{{\psi }^{\prime }}_2\left(z_0,\lambda \right){\delta }^{\left(2\right)}\left(z-z_0\right),\end{array}$ (4.18)

and take into account as non-vanishing parameters the moments $u_{\mp \mathrm{,0}}\left(z_0\right)\mathrm{,}{u}_{\mp \mathrm{,1}}\left(z_0\right)$ and $u_{\mp \mathrm{,2}}\left(z_0\right)$. In this approximation follows from (4.1) and (4.2) and from (4.3) for $n=2,3,4$

$\begin{array}{l}0={\psi }_{-}^{\left(1\right)}\left(z_0\mathrm{,}\lambda \right)-{\psi }_{+}^{\left(1\right)}\left(z_0\mathrm{,}\lambda \right)\\ +\left(u_{-\mathrm{,0}}\left(z_0\right){\psi }_{-}\left(z_0\mathrm{,}\lambda \right)+u_{+\mathrm{,0}}\left(z_0\right){\psi }_{+}\left(z_0\mathrm{,}\lambda \right)+\left(u\left(z_0\right)-\lambda \right){{\psi }^{\prime }}_0\left(z_0\mathrm{,}\lambda \right)\right)\\ -\left(u_{-\mathrm{,1}}\left(z_0\right){\psi }_{-}^{\left(1\right)}\left(z_0\mathrm{,}\lambda \right)+u_{+\mathrm{,1}}\left(z_0\right){\psi }_{+}^{\left(1\right)}\left(z_0\mathrm{,}\lambda \right)+u^{\left(1\right)}\left(z_0\right){{\psi }^{\prime }}_1\left(z_0\mathrm{,}\lambda \right)\right)\\ +\left(u_{-\mathrm{,2}}\left(z_0\right){\psi }_{-}^{\left(2\right)}\left(z_0\mathrm{,}\lambda \right)+u_{+\mathrm{,2}}\left(z_0\right){\psi }_{+}^{\left(2\right)}\left(z_0\mathrm{,}\lambda \right)+u^{\left(2\right)}\left(z_0\right){{\psi }^{\prime }}_2\left(z_0\mathrm{,}\lambda \right)\right)\mathrm{,}\end{array}$

$\begin{array}{l}0={\psi }_{-}\left(z_0\mathrm{,}\lambda \right)-{\psi }_{+}\left(z_0\mathrm{,}\lambda \right)\\ +\left(u_{-\mathrm{,1}}\left(z_0\right){\psi }_{-}\left(z_0\mathrm{,}\lambda \right)+u_{+\mathrm{,1}}\left(z_0\right){\psi }_{+}\left(z_0\mathrm{,}\lambda \right)+\left(u\left(z_0\right)-\lambda \right){{\psi }^{\prime }}_1\left(z_0\mathrm{,}\lambda \right)\right)\\ -2\left(u_{-\mathrm{,2}}\left(z_0\right){\psi }_{-}^{\left(1\right)}\left(z_0\mathrm{,}\lambda \right)+u_{+\mathrm{,2}}\left(z_0\right){\psi }_{+}^{\mathrm{<mo\; stretchy="false">\; (\; </mo>1\; <mo\; stretchy="false">\; )\; </mo>}}\left(z_0\mathrm{,}\lambda \right)+u^{\left(1\right)}\left(z_0\right){{\psi }^{\prime }}_2\left(z_0\mathrm{,}\lambda \right)\right)\mathrm{,}\end{array}$

$0=\left(u_{-\mathrm{,2}}\left(z_0\right){\psi }_{-}\left(z_0\mathrm{,}\lambda \right)+u_{+\mathrm{,2}}\left(z_0\right){\psi }_{+}\left(z_0\mathrm{,}\lambda \right)+\left(u\left(z_0\right)-\lambda \right){{\psi }^{\prime }}_2\left(z_0\mathrm{,}\lambda \right)\right)-{{\psi }^{\prime }}_0\left(z_0\mathrm{,}\lambda \right)\mathrm{,}$

$0=-{{\psi }^{\prime }}_1\left(z_0,\lambda \right),$

$0=-{{\psi }^{\prime }}_2\left(z_0,\lambda \right).$ (4.19)

In these boundary conditions we have for the first time two new peculiarities. The first is the occurrence of the second-order derivatives ${\psi }_{\mp }^{\left(2\right)}\left(z_0\mathrm{,}\lambda \right)$ of the asymptotic solutions at the boundary z₀. Due to the asymptotic Equations (3.2) one may substitute them according to

${\psi }_{\mp }^{\left(2\right)}\left(z_0\mathrm{,}\lambda \right)=\left(u_{\mp }\left(z_0\right)-\lambda \right){\psi }_{\mp }\left(z_0\mathrm{,}\lambda \right)\mathrm{.}$ (4.20)

The second is the non-vanishing of the coefficients ${{\psi }^{\prime }}_0\left(z_0\mathrm{,}\lambda \right)$ for the first time. Taking into account the vanishing of ${\psi }_2\left(z_0\mathrm{,}\lambda \right)$ we can express them by

${{\psi }^{\prime }}_0\left(z_0,\lambda \right)=u_{-,2}\left(z_0\right){\psi }_{-}\left(z_0,\lambda \right)+u_{+,2}\left(z_0\right){\psi }_{+}\left(z_0,\lambda \right).$ (4.21)

Inserting this into (4.19) we obtain the following two effective boundary conditions

$\begin{array}{l}0=\left(1-u_{-,1}\left(z_0\right)\right){\psi }_{-}^{\left(1\right)}\left(z_0,\lambda \right)-\left(1+u_{+,1}\left(z_0\right)\right){\psi }_{+}^{\left(1\right)}\left(z_0,\lambda \right)\\ +\left(u_{-,0}\left(z_0\right)+\left(u\left(z_0\right)+u_{-}\left(z_0\right)-2\lambda \right)u_{-,2}\left(z_0\right)\right){\psi }_{-}\left(z_0,\lambda \right)\\ +\left(u_{+,0}\left(z_0\right)+\left(u\left(z_0\right)+u_{+}\left(z_0\right)-2\lambda \right)u_{+,2}\left(z_0\right)\right){\psi }_{+}\left(z_0,\lambda \right),\end{array}$

$\begin{array}{l}0=\left(1+u_{-,1}\left(z_0\right)\right){\psi }_{-}\left(z_0,\lambda \right)-\left(1-u_{+,1}\left(z_0\right)\right){\psi }_{+}\left(z_0,\lambda \right)\\ -2\left(u_{-,2}\left(z_0\right){\psi }_{-}^{\left(1\right)}\left(z_0,\lambda \right)+u_{+,2}\left(z_0\right){\psi }_{+}^{\left(1\right)}\left(z_0,\lambda \right)\right).\end{array}$ (4.22)

These are two boundary conditions for the asymptotic solutions ${\psi }_{\mp }\left(z\mathrm{,}\lambda \right)$ taken for $z=z_0$ at the boundary, i.e., for ${\psi }_{\mp }\left(z_0\mathrm{,}\lambda \right)$ and their first derivatives ${\psi }_{\mp }^{\left(1\right)}\left(z_0\mathrm{,}\lambda \right)$ at the boundary with the moments $u_{\mp \mathrm{,0}}\left(z_0\right)\mathrm{,}{u}_{\mp \mathrm{,1}}\left(z_0\right)\mathrm{,}{u}_{\mp \mathrm{,2}}\left(z_0\right)$. The quantities $u\left(z_0\right)\mathrm{,}{u}_{\mp }\left(z_0\right)$ involved have to be considered as known or as calculable parameters from the potential $u\left(z\right)$. The mentioned peculiarities with the additional Equation (4.21) for ${{\psi }^{\prime }}_0\left(z_0\mathrm{,}\lambda \right)$ are the proper reason that we explicitly wrote down the already complicated boundary conditions (4.22).

Higher-order approximations of the boundary condition than the considered third-order approximation become much more difficult since in these cases one has to insert also higher- than second-order derivatives of the asymptotic potentials ${\psi }_{\mp }\left(z\mathrm{,}\lambda \right)$ at the boundary $z=z_0$. Using (3.2) they can be expressed by zeroth- and first-order derivatives of ${\psi }_{\mp }\left(z\mathrm{,}\lambda \right)$ plus derivatives of $u_{\mp }\left(z\right)$ that, in principle, is makable but becomes very complicated.

In many cases it is possible to make further simplifications of the boundary conditions if we make suppositions about the symmetry properties of the deviations $\Delta u\left(z_0\mathrm{<mo>\; ;\; </mo>}z\right)$ of the potentials in the transition regions defined in (3.10) in connection with assumptions of monotonic increasing or decreasing.

5. Reflection and Refraction at Transition Layers between Two Asymptotically Constant Potentials in Second-Order Approximation

For the illustration of the derived boundary conditions and as important application we consider now problems of reflection and refraction for incident waves at a transition layer. Our notations for the amplitudes of the corresponding waves are schematically shown in Figure 5.

One of the simplest asymptotic potentials are potentials which go rapidly to constant values $u_{-}\left(z\right)\to u_{-}$ as $z\to -\infty$ and to $u_{+}\left(z\right)\to u_{+}$ as $z\to +\infty$ and, therefore, where the deviations $\Delta u\left(z_0\mathrm{<mo>\; ;\; </mo>}z\right)$ from these potentials are rapidly vanishing for $z\to \pm \infty$ in the inner region. The lengths in the inner region where the potential changes essentially must be small in comparison to characteristic lengths of the problem (e.g., wavelengths). Similar problems are sometimes calculated for concrete transition potentials, in case of potential barriers or potential wells for rectangular potentials. The WKB-approximation where the lengths of changes of the potentials are large in comparison to the characteristic lengths of the problem is the opposite case to our calculations (e.g., [20] ).

Thus we consider now the following asymptotic potentials $u_{-}\left(z\right)$ and $u_{+}(\; z\; )$

$u_{-}\left(z\to -\infty \right)\equiv u_{-}\mathrm{,}{u}_{+}\left(z\to -\infty \right)\equiv u_{+}\mathrm{.}$ (5.1)

Only solutions with the same value of the parameter $\lambda$ are coupled at the boundary. The general asymptotic solutions of Equation (3.2) for the constant

The upper indices “i” and “r” mean “incident” and “reflected or refracted”, correspondingly. The lower indices “−” and “+” are chosen to see in equations which part on the coordinate axis is meant.

Figure 5. Notations of amplitudes of incident and reflected or refracted waves.

asymptotic potentials $u_{\mp }$ for each value of the parameter $\lambda$ are superpositions of two independent basic solutions and we set

$\begin{array}{l}{\psi }_{-}\left(z\mathrm{,}\lambda \right)=A_{-}^{\text{i}}{\text{e}}^{\text{i}{k}_{-}z}+A_{-}^{\text{r}}{\text{e}}^{-\text{i}{k}_{-}z}\mathrm{,}\\ {\psi }_{+}\left(z\mathrm{,}\lambda \right)=A_{+}^{\text{i}}{\text{e}}^{-\text{i}{k}_{+}z}+A_{+}^{\text{r}}{\text{e}}^{\text{i}{k}_{+}z}\mathrm{,}\end{array}$ (5.2)

with the dispersion relations (in principle, also true for complex values $\lambda$ )

$\begin{array}{l}{k}_{-}=\sqrt{\lambda -u_{-}},k_{+}=\sqrt{\lambda -u_{+}},\\ \iff \lambda =k_{-}^2+u_{-}=k_{+}^2+u_{+},u_{-}-u_{+}=k_{+}^2-k_{-}^2.\end{array}$ (5.3)

The complex amplitudes of the two possible waves on each side of the boundary $z=z_0$ are denoted by $A_{-}^{\text{i}\mathrm{,}\text{r}}$ and $A_{+}^{\text{i}\mathrm{,}\text{r}}$ (upper index “i” stands for “incident” and upper index “r” for “reflected” or “refracted” in dependence on the side from which the incident and reflected or refracted waves are considered but only for real $\lambda$ and then chosen $k_{-}>0$ and $k_{+}>0$ these notations are uniquely justified). From (5.2) follows for the derivatives of the solutions

$\begin{array}{l}{\psi }_{-}^{\left(1\right)}\left(z\mathrm{,}\lambda \right)=+\text{i}{k}_{-}\left(A_{-}^{\text{i}}{\text{e}}^{\text{i}{k}_{-}z}-A_{-}^{\text{r}}{\text{e}}^{-\text{i}{k}_{-}z}\right)\mathrm{,}\\ {\psi }_{+}^{\left(1\right)}\left(z\mathrm{,}\lambda \right)=-\text{i}{k}_{+}\left(A_{+}^{\text{i}}{\text{e}}^{-\text{i}{k}_{+}z}-A_{+}^{\text{r}}{\text{e}}^{\text{i}{k}_{+}z}\right)\mathrm{.}\end{array}$ (5.4)

For $\lambda \ge u_{\mp }$ both solutions $k_{\mp }$ possess real values and denote then wavelengths and for $\lambda \le u_{\mp }$ both solutions $k_{\mp }$ possess imaginary values and denote decay constants. For $u_{-}<\lambda <u_{+}$ or $u_{-}>\lambda >u_{+}$ only one of the parameters $k_{\mp }$ is real and the other $k_{\pm }$ is imaginary, respectively.

The derived boundary conditions at $z=z_0$ in arbitrary approximation can be written in form of a matrix equation, for example, in the following way

${\text{S}}^{\text{i}}\left(\begin{array}{c}{A}_{-}^{\text{i}}{\text{e}}^{+\text{i}{k}_{-}{z}_0}\\ A_{+}^{\text{i}}{\text{e}}^{-\text{i}{k}_{+}{z}_0}\end{array}\right)={\text{S}}^{\text{r}}\left(\begin{array}{c}{A}_{-}^{\text{r}}{\text{e}}^{-\text{i}{k}_{-}{z}_0}\\ A_{+}^{\text{r}}{\text{e}}^{+\text{i}{k}_{+}{z}_0}\end{array}\right)\mathrm{,}$ (5.5)

where ${\text{S}}^{\text{i}\mathrm{,}\text{r}}$ are matrices depending on the considered approximations. These matrices ${\text{S}}^{\text{i}\mathrm{,}\text{r}}$ (or ${\text{S}}_{\mp }$ ) become more complicated with each higher approximation of the boundary conditions and the jump from the second approximation (4.17) to the third approximation (4.22) is already very large and leads to complicated expression.

In the second-order approximation (4.17) we find for the considered asymptotically constant potentials the following matrices in case of (5.9)

$\begin{array}{l}{\text{S}}^{\text{i}}=\left(\begin{array}{cc}\left(1-u_{-\mathrm{,1}}\left(z_0\right)\right)k_{-}-\text{i}{u}_{-\mathrm{,0}}\left(z_0\right)& \left(1+u_{+\mathrm{,1}}\left(z_0\right)\right)k_{+}-\text{i}{u}_{+\mathrm{,0}}\left(z_0\right)\\ 1+u_{-\mathrm{,1}}\left(z_0\right)& -1+u_{+\mathrm{,1}}\left(z_0\right)\end{array}\right)\mathrm{,}\\ {\text{S}}^{\text{r}}=\left(\begin{array}{cc}\left(1-u_{-\mathrm{,1}}\left(z_0\right)\right)k_{-}+\text{i}{u}_{-\mathrm{,0}}\left(z_0\right)& \left(1+u_{+\mathrm{,1}}\left(z_0\right)\right)k_{+}+\text{i}{u}_{+\mathrm{,0}}\left(z_0\right)\\ -1-u_{-\mathrm{,1}}\left(z_0\right)& 1-u_{+\mathrm{,1}}\left(z_0\right)\end{array}\right)\mathrm{.}\end{array}$ (5.6)

Due to two-dimensional Hamilton-Cayley identity of an arbitrary two-dimensional matrix $\text{S}$

$S^2-\langle S\rangle S+\left[S\right]I=\mathrm{0,}\left[S\right]\equiv \frac{1}{2}\left({\langle S\rangle }^2-\langle S^2\rangle \right)\mathrm{,}$ (5.7)

$S\equiv \left(\begin{array}{cc}{S}_{11}\mathrm{,}& S_{12}\\ S_{21}\mathrm{,}& S_{22}\end{array}\right)\mathrm{,}\bar{S}\equiv \left(\begin{array}{rr}\hfill S_{22}\mathrm{,}& \hfill -S_{12}\\ \hfill -S_{21}\mathrm{,}& \hfill S_{11}\end{array}\right)\mathrm{,}\langle S\rangle \equiv S_{11}+S_{22}\mathrm{,}\left[S\right]\equiv S_{11}{S}_{22}-S_{12}{S}_{21}\mathrm{.}$

However, already in 3D-case the advantage is mostly important and we prefer them everywhere.

where $\langle S\rangle$ is the trace and $\left[S\right]$ the determinant of $\text{S}$ and $I$ denotes the identity matrix one finds the inverse matrix $S^{-1}$ in coordinate-invariant way by⁴

$\bar{S}\equiv \langle S\rangle I-S\mathrm{,}\Rightarrow \langle \bar{S}\rangle =\langle S\rangle \mathrm{,}\left[\bar{S}\right]=\left[S\right]\mathrm{,}\bar{S}S=\left[S\right]I\mathrm{,}{S}^{-1}=\frac{\bar{S}}{\left[S\right]}\mathrm{,}$ (5.8)

where $\bar{S}$ by definition denotes the associated matrix to the matrix $\text{S}$.

By multiplication of (5.6) with the inverse matrix ${\left(S^{\text{r}}\right)}^{-1}$ to $S^{\text{r}}$ one may resolve the equations to the reflected or refracted wave amplitudes in dependence on the amplitudes of the incident waves.

$\left(\begin{array}{c}{A}_{-}^{\text{r}}{\text{e}}^{-\text{i}{k}_{-}{z}_0}\\ A_{+}^{\text{r}}{\text{e}}^{+\text{i}{k}_{+}{z}_0}\end{array}\right)=\frac{\overline{{\text{S}}^{\text{r}}}{\text{S}}^{\text{i}}}{\left[{\text{S}}^{\text{r}}\right]}\left(\begin{array}{c}{A}_{-}^{\text{i}}{\text{e}}^{+\text{i}{k}_{-}{z}_0}\\ A_{+}^{\text{i}}{\text{e}}^{-\text{i}{k}_{+}{z}_0}\end{array}\right)\equiv \text{R}\left(\begin{array}{c}{A}_{-}^{\text{i}}{\text{e}}^{+\text{i}{k}_{-}{z}_0}\\ A_{+}^{\text{i}}{\text{e}}^{-\text{i}{k}_{+}{z}_0}\end{array}\right)=\left(\begin{array}{cc}{R}_{--}& R_{-+}\\ R_{+-}& R_{++}\end{array}\right)\left(\begin{array}{c}{A}_{-}^{\text{i}}{\text{e}}^{+\text{i}{k}_{-}{z}_0}\\ A_{+}^{\text{i}}{\text{e}}^{-\text{i}{k}_{+}{z}_0}\end{array}\right)\mathrm{.}$ (5.9)

We find for the numerators of the matrix elements of the reflection-refraction matrix $\text{R}$

$\begin{array}{l}\left[{\text{S}}^{\text{r}}\right]R_{--}=\left(1+u_{-\mathrm{,1}}\left(z_0\right)u_{+\mathrm{,1}}\left(z_0\right)\right)\left(k_{-}-k_{+}\right)-\left(u_{-\mathrm{,1}}\left(z_0\right)+u_{+\mathrm{,1}}\left(z_0\right)\right)\left(k_{-}+k_{+}\right)\\ -\text{i}\left(u_{-\mathrm{,0}}\left(z_0\right)+u_{+\mathrm{,0}}\left(z_0\right)-u_{-\mathrm{,0}}\left(z_0\right)u_{+\mathrm{,1}}\left(z_0\right)+u_{+\mathrm{,0}}\left(z_0\right)u_{-\mathrm{,1}}\left(z_0\right)\right)\mathrm{,}\end{array}$

$\left[{\text{S}}^{\text{r}}\right]R_{-+}=2\left(1-u_{+\mathrm{,1}}^2\left(z_0\right)\right)k_{+}\mathrm{,}$

$\left[{\text{S}}^{\text{r}}\right]R_{+-}=2\left(1-u_{-\mathrm{,1}}^2\left(z_0\right)\right)k_{-}\mathrm{,}$

$\begin{array}{l}\left[{\text{S}}^{\text{r}}\right]R_{++}=-\left(1+u_{-\mathrm{,1}}\left(z_0\right)u_{+\mathrm{,1}}\left(z_0\right)\right)\left(k_{-}-k_{+}\right)+\left(u_{-\mathrm{,1}}\left(z_0\right)+u_{+\mathrm{,1}}\left(z_0\right)\right)\left(k_{-}+k_{+}\right)\\ -\text{i}\left(u_{-\mathrm{,0}}\left(z_0\right)+u_{+\mathrm{,0}}\left(z_0\right)-u_{-\mathrm{,0}}\left(z_0\right)u_{+\mathrm{,1}}\left(z_0\right)+u_{+\mathrm{,0}}\left(z_0\right)u_{-\mathrm{,1}}\left(z_0\right)\right)\mathrm{,}\end{array}$ (5.10)

and for the determinant of the matrix ${\text{S}}^{\text{r}}$

$\begin{array}{c}\left[{\text{S}}^{\text{r}}\right]=\left(1+u_{-\mathrm{,1}}\left(z_0\right)u_{+\mathrm{,1}}\left(z_0\right)\right)\left(k_{-}+k_{+}\right)-\left(u_{-\mathrm{,1}}\left(z_0\right)+u_{+\mathrm{,1}}\left(z_0\right)\right)\left(k_{-}-k_{+}\right)\\ +\text{i}\left(u_{-\mathrm{,0}}\left(z_0\right)+u_{+\mathrm{,0}}\left(z_0\right)-u_{-\mathrm{,0}}\left(z_0\right)u_{+\mathrm{,1}}\left(z_0\right)+u_{+\mathrm{,0}}\left(z_0\right)u_{-\mathrm{,1}}\left(z_0\right)\right)\mathrm{.}\end{array}$ (5.11)

Already into this second-order approximation the moments $u_{-\mathrm{,}\nu }\left(z_0\right)$ and $u_{+\mathrm{,}\nu }\left(z_0\right)$ are separately involved and not only as their sums $u_{-\mathrm{,}\nu }\left(z_0\right)+u_{+\mathrm{,}\nu }\left(z_0\right)\mathrm{,}\left(\nu =\mathrm{0,1}\right)$. The formulae (5.9) together with (5.10) and (5.11) describe the already complicated second-order approximation of reflection and refraction at a transition layer between two asymptotically constant potentials $u_{-}$ and $u_{+}$ which are contained in the expressions (5.3) for $k_{-}$ and $k_{+}$.

6. First-Order Approximation of Reflection and Refraction Formulae and Their Invariance with Small Variation of Boundary Choice

In first-order approximation we neglect the moments $u_{\mp \mathrm{,1}}\left(z_0\right)$ and obtain from (5.10) and (5.11) the following formulae for the amplitudes of the reflected and refracted waves in dependence on the amplitudes of the incident waves

$\left(\begin{array}{c}{A}_{-}^{\text{r}}{\text{e}}^{-\text{i}{k}_{-}{z}_0}\\ A_{+}^{\text{r}}{\text{e}}^{+\text{i}{k}_{+}{z}_0}\end{array}\right)=\left(\begin{array}{cc}\frac{k_{-}-k_{+}-\text{i}{u}_0\left(z_0\right)}{k_{-}+k_{+}+\text{i}{u}_0\left(z_0\right)}& \frac{2k_{+}}{k_{-}+k_{+}+\text{i}{u}_0\left(z_0\right)}\\ \frac{2k_{-}}{k_{-}+k_{+}+\text{i}{u}_0\left(z_0\right)}& \frac{k_{+}-k_{-}-\text{i}{u}_0\left(z_0\right)}{k_{-}+k_{+}+\text{i}{u}_0\left(z_0\right)}\end{array}\right)\left(\begin{array}{c}{A}_{-}^{\text{i}}{\text{e}}^{+\text{i}{k}_{-}{z}_0}\\ A_{+}^{\text{i}}{\text{e}}^{-\text{i}{k}_{+}{z}_0}\end{array}\right)\mathrm{,}$ (6.1)

with the inversion

$\left(\begin{array}{c}{A}_{-}^{\text{i}}{\text{e}}^{+\text{i}{k}_{-}{z}_0}\\ A_{+}^{\text{i}}{\text{e}}^{-\text{i}{k}_{+}{z}_0}\end{array}\right)=\left(\begin{array}{cc}\frac{k_{-}-k_{+}+\text{i}{u}_0\left(z_0\right)}{k_{-}+k_{+}-\text{i}{u}_0\left(z_0\right)}& \frac{2k_{+}}{k_{-}+k_{+}-\text{i}{u}_0\left(z_0\right)}\\ \frac{2k_{-}}{k_{-}+k_{+}-\text{i}{u}_0\left(z_0\right)}& \frac{k_{+}-k_{-}+\text{i}{u}_0\left(z_0\right)}{k_{-}+k_{+}-\text{i}{u}_0\left(z_0\right)}\end{array}\right)\left(\begin{array}{c}{A}_{-}^{\text{r}}{\text{e}}^{-\text{i}{k}_{-}{z}_0}\\ A_{+}^{\text{r}}{\text{e}}^{+\text{i}{k}_{+}{z}_0}\end{array}\right)\mathrm{.}$ (6.2)

The moment $u_0\left(z_0\right)$ is the sum of the moments $u_{-\mathrm{,0}}\left(z_0\right)$ and $u_{+\mathrm{,0}}\left(z_0\right)$ and corresponds to the integral over the deviation of the genuine potential from the asymptotic potentials (see Figure 1 and Figure 2)

$\begin{array}{c}{u}_0\left(z_0\right)\equiv u_{-\mathrm{,0}}\left(z_0\right)+u_{+\mathrm{,0}}\left(z_0\right)\\ \mathrm{<mo>\; =\; </mo>}{\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}z\left(u\left(z\right)-u_{-}\theta \left(z_0-z\right)-u_{+}\theta \left(z-z_0\right)\right)\mathrm{.}\end{array}$ (6.3)

Due to constant asymptotic potentials its Taylor series in powers of z₀ possesses only two terms (see also (4.13))

$u_0\left(z_0\right)=u_0\left(0\right)-\left(u_{-}-u_{+}\right)z_0=u_0\left(0\right)-\left(k_{+}^2-k_{-}^2\right)z_0\mathrm{,}$ (6.4)

where (5.3) is taken into account. For $u_{-}-u_{+}\ne 0$ (genuine transition between asymptotically different potentials) the reference point z₀ of the expansions (the “boundary”) can be uniquely chosen in such way that the moment $u_0\left(z_0\right)$ vanishes leading to

$u_0\left(z_0\right)=\mathrm{0,}\Rightarrow z_0=\frac{u_0\left(0\right)}{u_{-}-u_{+}}=\frac{u_0\left(0\right)}{k_{+}^2-k_{-}^2}\mathrm{.}$ (6.5)

This is likely the “best” choice of the point z₀ since in this case the matrix of reflection and refraction takes on the usual form without a transition layer and the presence of the last is then involved in the additional phase factors at the amplitudes in (6.1) in comparison to the usual formulae in case of step-wise (abrupt) transition at $z_0=0$. Formally, with substitutions they are even identical with the relations between the amplitudes of the electric field of incident, reflected and refracted waves for oblique incidence when the electric field is polarized perpendicular to the incidence plane (s-waves; see [14], chap. X, section 86, Eqs. (86.4)).

The matrix Equation (6.1) can be also represented in the form

$\left(\begin{array}{c}{A}_{-}^{\text{r}}\\ A_{+}^{\text{r}}\end{array}\right)=\left(\begin{array}{ll}\frac{k_{-}-k_{+}-\text{i}{u}_0\left(z_0\right)}{k_{-}+k_{+}+\text{i}{u}_0\left(z_0\right)}{\text{e}}^{+\text{i}2k_{-}{z}_0}\hfill & \frac{2k_{+}}{k_{-}+k_{+}+\text{i}{u}_0\left(z_0\right)}{\text{e}}^{+\text{i}\left(k_{-}-k_{+}\right)z_0}\hfill \\ \frac{2k_{-}}{k_{-}+k_{+}+\text{i}{u}_0\left(z_0\right)}{\text{e}}^{+\text{i}\left(k_{-}-k_{+}\right)z_0}\hfill & \frac{k_{+}-k_{-}-\text{i}{u}_0\left(z_0\right)}{k_{-}+k_{+}+\text{i}{u}_0\left(z_0\right)}{\text{e}}^{-\text{i}2k_{+}{z}_0}\hfill \end{array}\right)\left(\begin{array}{c}{A}_{-}^{\text{i}}\\ A_{+}^{\text{i}}\end{array}\right)\mathrm{.}$ (6.6)

We show now that the matrix in (6.6) which determines the relation between the amplitudes of the incident and reflected waves is nearly unsensible against small variations of the boundary z₀. For example, for the first matrix element in (6.6) we have

$\begin{array}{l}\frac{k_{-}-k_{+}-\text{i}{u}_0\left(z_0\right)}{k_{-}+k_{+}+\text{i}{u}_0\left(z_0\right)}{\text{e}}^{+\text{i}2k_{-}{z}_0}=\frac{\left\{k_{-}-k_{+}-\text{i}{u}_0\left(z_0\right)\right\}{\text{e}}^{+\text{i}\left(k_{-}+k_{+}\right)z_0}}{\left\{k_{-}+k_{+}+\text{i}{u}_0\left(z_0\right)\right\}{\text{e}}^{-\text{i}\left(k_{-}-k_{+}{z}_0\right)}}\\ \approx \frac{\left\{k_{-}-k_{+}-\text{i}\left(u_0\left(0\right)-\left(u_{-}-u_{+}\right)z_0\right)\right\}\left(1+\text{i}\left(k_{-}+k_{+}\right)z_0\right)}{\left\{k_{-}+k_{+}+\text{i}\left(u_0\left(0\right)-\left(u_{-}-u_{+}\right)z_0\right)\right\}\left(1-\text{i}\left(k_{-}-k_{+}\right)z_0\right)}\\ \approx \frac{k_{-}-k_{+}-\text{i}{u}_0\left(0\right)+\text{i}\left(u_{-}-u_{+}+\left(k_{-}^2-k_{+}^2\right)\right)z_0}{k_{-}+k_{+}+\text{i}{u}_0\left(0\right)-\text{i}\left(u_{-}-u_{+}+\left(k_{-}^2-k_{+}^2\right)\right)z_0}=\frac{k_{-}-k_{+}-\text{i}{u}_0\left(0\right)}{k_{-}+k_{+}+\text{i}{u}_0\left(0\right)}\mathrm{,}\end{array}$ (6.7)

due to (5.3) and where we neglected sum terms proportional to $z_0^2$ in numerator and denominator as “small” quantities. In full analogy we may treat the other matrix elements where in the non-diagonal elements we have only to transform the denominators. Thus if $z=z_0$ is already a “good” approximation for the boundary we may simplify (6.6) to

$\left(\begin{array}{c}{A}_{-}^{\text{r}}\\ A_{+}^{\text{r}}\end{array}\right)=\left(\begin{array}{cc}\frac{k_{-}-k_{+}-\text{i}{u}_0\left(0\right)}{k_{-}+k_{+}+\text{i}{u}_0\left(0\right)}& \frac{2k_{+}}{k_{-}+k_{+}+\text{i}{u}_0\left(0\right)}\\ \frac{2k_{-}}{k_{-}+k_{+}+\text{i}{u}_0\left(0\right)}& \frac{k_{+}-k_{-}-\text{i}{u}_0\left(0\right)}{k_{-}+k_{+}+\text{i}{u}_0\left(0\right)}\end{array}\right)\left(\begin{array}{c}{A}_{-}^{\text{i}}\\ A_{+}^{\text{i}}\end{array}\right)\mathrm{.}$ (6.8)

In view of the complexity of the second-order approximation of the reflection and refraction formulae (5.9), (5.10) and (5.11) we did not check them in analogous way with respect to their approximate invariance under small variations of the reference point z₀.

7. Semi-Bound and Bound States for Asymptotically of Constant Potentials

The formulae of the preceding Sections 5 and 6 are true, in principle, for complex values $\lambda$ in connection with (5.2) and (5.3).

For real values $\lambda$ in cases of $\lambda <u_{-}$ or $\lambda <u_{+}$ we find according to (5.3) imaginary values of $k_{-}$ or $k_{+}$ and we set in these special cases, respectively

$k_{-}\to \text{i}{\kappa }_{-}\mathrm{,}{k}_{+}\to \text{i}{\kappa }_{+}\mathrm{,}$ (7.1)

corresponding to exponentially decreasing or increasing factors for $z\to \mp \infty$. From these factors only ${\text{e}}^{+{\kappa }_{-}z}$ for $z\to -\infty$ and ${\text{e}}^{-{\kappa }_{+}z}$ for $z\to +\infty$ may correspond to finite sources in “infinity”. This implies that in addition to the considered general case we have to distinguish interesting special cases with exponential vanishing of the solutions in infinity. If we concretely assume now $u_{-}\le u_{+}$, we have first the case of semi-bound states

$u_{-}\le \lambda \le u_{+}\mathrm{,}$

$k_{-}\equiv \sqrt{\lambda -u_{-}}\mathrm{,}{\kappa }_{+}\equiv -\text{i}{k}_{+}=\sqrt{u_{+}-\lambda }\ge \mathrm{0,}$ (7.2)

with the form of the solutions (we set $A_{+}^{\text{i}}=\mathrm{0,}{A}_{+}^{\text{r}}\equiv A_{+}$ )

${\psi }_{-}\left(\lambda \mathrm{,}z\right)=A_{-}^{\text{i}}{\text{e}}^{+\text{i}{k}_{-}z}+A_{-}^{\text{r}}{\text{e}}^{-\text{i}{k}_{-}z}\mathrm{,}{\psi }_{+}\left(z\mathrm{,}\lambda \right)=A_{+}{\text{e}}^{-{\kappa }_{+}z}\mathrm{.}$ (7.3)

The second interesting case concerns bound states

$\lambda \le u_{-}\le u_{+}\mathrm{,}$

${\kappa }_{-}\equiv -\text{i}{k}_{-}=\sqrt{u_{-}-\lambda }\ge \mathrm{0,}{\kappa }_{+}\equiv -\text{i}{k}_{+}=\sqrt{u_{+}-\lambda }\ge \mathrm{0,}$ (7.4)

with the form of the solutions (we set $A_{-}^{\text{i}}=A_{+}^{\text{i}}=0,A_{-}^{\text{r}}\equiv A_{-}\mathrm{,}{A}_{+}^{\text{r}}\equiv A_{+}$ )

${\psi }_{-}\left(\lambda \mathrm{,}z\right)=A_{-}{\text{e}}^{+{\kappa }_{-}z}\mathrm{,}{\psi }_{+}\left(\lambda \mathrm{,}z\right)=A_{+}{\text{e}}^{-{\kappa }_{+}z}\mathrm{.}$ (7.5)

In the semi-bound case in the first-order approximation of the boundary conditions we find from (6.1) for the amplitudes of the reflected wave $A_{-}^{\text{r}}$ on the negative side and the exponentially decreasing wave $A_{+}$ on the positive side in dependence on the amplitude of the incident wave $A_{-}^{\text{i}}$ on the negative side and in dependence on the reference point z₀

$\begin{array}{l}{A}_{-}^{\text{r}}{\text{e}}^{-\text{i}{k}_{-}{z}_0}=\frac{k_{-}-\text{i}{\kappa }_{+}-\text{i}{u}_0\left(z_0\right)}{k_{-}+\text{i}{\kappa }_{+}+\text{i}{u}_0\left(z_0\right)}{A}_{-}^{\text{i}}{\text{e}}^{\text{i}{k}_{-}{z}_0}\mathrm{,}\\ A_{+}{\text{e}}^{-{\kappa }_{+}{z}_0}=\frac{2k_{-}}{k_{-}+\text{i}{\kappa }_{+}+\text{i}{u}_0\left(z_0\right)}{A}_{-}^{\text{i}}{\text{e}}^{\text{i}{k}_{-}{z}_0}\mathrm{,}\end{array}$ (7.6)

with

$k_{-}\equiv \sqrt{\lambda -u_{-}}\ge \mathrm{0,}{\kappa }_{+}\equiv \sqrt{u_{+}-\lambda }\ge 0.$ (7.7)

Semi-bound states as seen from (7.2) exist as solutions of the Sturm-Liouville equation for asymptotically constant potentials for a whole band of values $\lambda$ in the range $u_{-}\le \lambda \le u_{+}$ and therefore are not discrete states. They are non-degenerate that means for each possible value $\lambda$ in the given range exists one state. In some aspects they are analogous to total reflection in optics.

We come now to bound states which are characterized by the condition (7.4). The dispersion equation can be obtained from vanishing of the denominators in the matrix elements of in (6.1) but the amplitude relation are difficult to find directly from them. Therefore we go back for the case $\lambda \le u_{-}\le u_{+}$ directly to the boundary relations in first-order approximation with the proposition for the solution

$\begin{array}{l}{\psi }_{-}\left(z\mathrm{,}\lambda \right)=A_{-}{\text{e}}^{+{\kappa }_{-}z}\mathrm{,}{\psi }_{-}^{\left(1\right)}\left(z\mathrm{,}\lambda \right)=+{\kappa }_{-}{A}_{-}{\text{e}}^{+{\kappa }_{-}z}\mathrm{,}\\ {\psi }_{+}\left(z\mathrm{,}\lambda \right)=A_{+}{\text{e}}^{-{\kappa }_{+}z}\mathrm{,}{\psi }_{+}^{\left(1\right)}\left(z\mathrm{,}\lambda \right)=-{\kappa }_{-}{A}_{+}{\text{e}}^{-{\kappa }_{+}z}\mathrm{,}\end{array}$ (7.8)

with the relations (7.4)

${\kappa }_{+}=\sqrt{u_{+}-\lambda }\ge {\kappa }_{-}=\sqrt{u_{-}-\lambda }\ge 0.$ (7.9)

Inserting (7.8) into the boundary conditions (4.11) we find written in matrix form

$0=\left(\begin{array}{cc}{\kappa }_{-}+\frac{u_0\left(z_0\right)}{2}& {\kappa }_{+}+\frac{u_0\left(z_0\right)}{2}\\ -1& 1\end{array}\right)\left(\begin{array}{c}{A}_{-}{\text{e}}^{+{\kappa }_{-}{z}_0}\\ A_{+}{\text{e}}^{-{\kappa }_{+}{z}_0}\end{array}\right)\mathrm{.}$ (7.10)

Condition for the existence of non-trivial solutions is the vanishing of the determinant of the matrix in (7.10), i.e.

$0={\kappa }_{-}+{\kappa }_{+}+u_0\left(z_0\right)\mathrm{,}\Rightarrow u_0\left(z_0\right)=-{\kappa }_{-}-{\kappa }_{+}\le \mathrm{0,}$ (7.11)

and for the amplitude relations using (7.11) we obtain

$A_{-}{\text{e}}^{+{\kappa }_{-}{z}_0}=A_{+}{\text{e}}^{-{\kappa }_{+}{z}_0}\mathrm{.}$ (7.12)

We now resolve the dispersion Equation (7.11) with respect to $\lambda$ using (7.9). Squaring ${\kappa }_{+}+{\kappa }_{-}$ from (7.11) follows

$u_0^2\left(z_0\right)-\left({\kappa }_{+}^2+{\kappa }_{-}^2\right)=2{\kappa }_{+}{\kappa }_{-}\mathrm{.}$ (7.13)

Squaring both sides of this equation again we arrive at an equation which using (7.9) is already rational in $\lambda$ and can be written

$\begin{array}{c}0=u_0^4\left(z_0\right)-2\left({\kappa }_{+}^2+{\kappa }_{-}^2\right)u_0^2\left(z_0\right)+{\left({\kappa }_{+}^2-{\kappa }_{-}^2\right)}^2\\ =u_0^4\left(z_0\right)-2\left(u_{+}+u_{-}-2\lambda \right)u_0^2\left(z_0\right)+{\left(u_{+}-u_{-}\right)}^2\mathrm{.}\end{array}$ (7.14)

From this we find in first-order approximation of the boundary conditions only one solution for $\lambda$

$\lambda =-\frac{u_0^4\left(z_0\right)-2\left(u_{+}+u_{-}\right)u_0^2\left(z_0\right)+{\left(u_{+}-u_{-}\right)}^2}{4u_0^2\left(z_0\right)}\mathrm{.}$ (7.15)

This solution is below the lower of the values of $u_{-}$ and $u_{+}$, in considered case (7.4) lower than $u_{-}$ that is also seen from the following representation of (7.15)

$\lambda -u_{-}=-{\left(\frac{u_0^2\left(z_0\right)-\left(u_{+}-u_{-}\right)}{2u_0\left(z_0\right)}\right)}^2\le 0.$ (7.16)

According to (7.9) inserting (7.15) for $\lambda$ and choosing the signs of the arising square roots in the way that corresponding to (7.11) the moment $u_0\left(z_0\right)$ defined in (6.3) is non-positive this leads to

${\kappa }_{+}=-\frac{u_0^2\left(z_0\right)+\left(u_{+}-u_{-}\right)}{2u_0\left(z_0\right)}\ge {\kappa }_{-}=-\frac{u_0^2\left(z_0\right)-\left(u_{+}-u_{-}\right)}{2u_0\left(z_0\right)}\ge \mathrm{0,}$ (7.17)

and this unambiguously is nearer to a potential well than to a potential barrier (next Section).

The bound state (excluding here its limiting case ${\kappa }_{+}\ge {\kappa }_{-}=0$ )

$\psi \left(z\mathrm{,}\lambda \right)=A_{-}{\text{e}}^{{\kappa }_{-}z}\theta \left(z_0-z\right)+A_{+}{\text{e}}^{-{\kappa }_{+}z}\theta \left(z-z_0\right)\mathrm{,}$ (7.18)

is normalizable. As necessary and sufficient conditions for its presence follows from (7.11) and (7.17)

$u_0\left(z_0\right)<0,u_0^2\left(z_0\right)>u_{+}-u_{-}>0.$ (7.19)

The state (7.18) is concentrated to the neighborhood of the considered boundary z₀. In this sense it is analogous to surface waves at a boundary between two media in electrodynamics only with the difference that in last case there is also a real propagation constant parallel to the surface.

8. Asymptotically Equal Potentials and Bound States within a Potential Well

We consider now the special case of equal asymptotically constant potentials $u_{-}=u_{+}$

$u_{-}=u_{+}\equiv u_{\infty },$ (8.1)

from which follows that the first moment $u_0\left(z_0\right)\equiv u_{-}\left(z_0\right)+u_{+}\left(z_0\right)$ of the deviation of the real potential from the asymptotic potentials becomes independent on the reference point z₀, therefore

$u_0\left(z_0\right)\equiv {\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}z\left(u\left(z\right)-u_{\infty }\right)=u_0\left(0\right)\mathrm{.}$ (8.2)

For $u_0\left(0\right)>0$ one may speak about a potential barrier and for $u_0\left(0\right)<0$ about a potential well but one may use these notions also in a stronger sense of $u\left(z\right)-u_{\infty }>0$ and of $u\left(z\right)-u_{\infty }<0$ only for all z.

All formulae derived in the last two Sections for the general case $u_{-}\ne u_{+}$ remain true with the specialization (8.1) and only some formulae where $u_{-}-u_{+}$ stands in the denominator (such as (6.5) for a good choice of the reference point z₀) are no more applicable. For a “good” choice of the reference point z₀ one may choose instead that the first moment $u_1\left(z_0\right)$ defined now

$u_1\left(z_0\right)\equiv u_{-,1}\left(z_0\right)+u_{+,1}\left(z_0\right)=-{\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}z\left(u\left(z\right)-u_{\infty }\right)\left(z-z_0\right)=u_1\left(0\right)+z_0{u}_0\left(0\right),$ (8.3)

vanishes, from which follows for such choice

$0=u_1\left(z_0\right)=u_1\left(0\right)+z_0{u}_0\left(0\right),\Rightarrow z_0=-\frac{u_1\left(0\right)}{u_0\left(0\right)}.$ (8.4)

The normalizable bound states which vanish in infinity are now represented by (see (7.8))

$\begin{array}{l}{\psi }_{-}\left(z\mathrm{,}\lambda \right)=A_{-}{\text{e}}^{+\kappa z}\mathrm{,}{\psi }_{+}\left(z\mathrm{,}\lambda \right)=A_{+}{\text{e}}^{-\kappa z}\mathrm{,}\\ {\psi }_{-}^{\left(1\right)}\left(z\mathrm{,}\lambda \right)=\kappa A_{-}{\text{e}}^{+\kappa z}\mathrm{,}{\psi }_{+}^{\left(1\right)}\left(z\mathrm{,}\lambda \right)=-\kappa A_{+}{\text{e}}^{-\kappa z}\mathrm{.}\end{array}$ (8.5)

For the only discrete state in considered first-order approximation of the boundary conditions follows then from (7.15)

$\lambda -u_{\infty }=-\frac{u_0^2\left(0\right)}{4}<0,\kappa =\sqrt{u_{\infty }-\lambda }=-\frac{u_0\left(0\right)}{2}>0.$ (8.6)

and the amplitude relations are

$A_{-}{\text{e}}^{+\kappa z_0}=A_{+}{\text{e}}^{-\kappa z_0}\mathrm{.}$ (8.7)

The parameter $\lambda$ of the discrete state is in every case below the asymptotically constant potential $u_{\infty }$. All potential wells with the same area $u_0\left(0\right)$ below the asymptotic value $u_{\infty }$ (see also FigureA1 in Appendix B) are in first approximation equivalent

We now calculate the equation for normalizable discrete eigenstates in the third-order approximation of the boundary conditions. With the form of the solutions (8.5) follows from the boundary conditions in third-order approximation (4.22) two homogeneous equations which can be written in the following matrix form

${\text{S}}^{\left(2\right)}\left(\begin{array}{c}{A}_{-}{\text{e}}^{+\kappa z_0}\\ A_{+}{\text{e}}^{-\kappa z_0}\end{array}\right)=\left(\begin{array}{c}0\\ 0\end{array}\right)\mathrm{,}$ (8.8)

with the matrix ( $\lambda =-{\kappa }^2$ )

$\begin{array}{l}{\text{S}}^{\left(3\right)}\equiv (\begin{array}{c}\left(1-u_{-\mathrm{,1}}\left(z_0\right)\right)\kappa +u_{-\mathrm{,0}}\left(z_0\right)+\left(u\left(z_0\right)-2\lambda \right)u_{-\mathrm{,2}}\left(z_0\right)\\ 1+u_{-\mathrm{,1}}\left(z_0\right)-2u_{-\mathrm{,2}}\left(z_0\right)\kappa \end{array}\\ \begin{array}{c}\left(1+u_{+\mathrm{,1}}\left(z_0\right)\right)\kappa +u_{+\mathrm{,0}}\left(z_0\right)+\left(u\left(z_0\right)-2\lambda \right)u_{+\mathrm{,2}}\left(z_0\right)\\ -1+u_{+\mathrm{,1}}\left(z_0\right)+2u_{+\mathrm{,2}}\left(z_0\right)\kappa \end{array})\mathrm{.}\end{array}$ (8.9)

The condition of solvability of (8.8) is the vanishing of the determinant of ${\text{S}}^{\left(3\right)}$ which leads to the following equation of third degree for $\kappa$

$\begin{array}{l}0=8u_{-\mathrm{,2}}\left(z_0\right)u_{+\mathrm{,2}}\left(z_0\right){\kappa }^3-4\left\{u_{-\mathrm{,1}}\left(z_0\right)u_{+\mathrm{,2}}\left(z_0\right)-u_{+\mathrm{,1}}\left(z_0\right)u_{-\mathrm{,2}}\left(z_0\right)\right\}{\kappa }^2\\ -2\{1+u_{-\mathrm{,1}}\left(z_0\right)u_{+\mathrm{,1}}\left(z_0\right)-u_{-\mathrm{,0}}\left(z_0\right)u_{+\mathrm{,2}}\left(z_0\right)-u_{+\mathrm{,0}}\left(z_0\right)u_{-\mathrm{,2}}\left(z_0\right)\\ -2u\left(z_0\right)u_{-\mathrm{,2}}\left(z_0\right)u_{+\mathrm{,2}}\left(z_0\right)\}\kappa -\{u_0\left(z_0\right)-u_{-\mathrm{,0}}\left(z_0\right)u_{+\mathrm{,1}}\left(z_0\right)\\ +u_{+\mathrm{,0}}\left(z_0\right)u_{-\mathrm{,1}}\left(z_0\right)+u\left(z_0\right)\left(u_2\left(z_0\right)+u_{-\mathrm{,1}}\left(z_0\right)u_{+\mathrm{,2}}\left(z_0\right)-u_{+\mathrm{,1}}\left(z_0\right)u_{-\mathrm{,2}}\left(z_0\right)\right)\}\mathrm{,}\end{array}$ (8.10)

and for the amplitude relations we find from (8.8)

$\left(1+u_{-\mathrm{,1}}\left(z_0\right)-2u_{-\mathrm{,2}}\left(z_0\right)\kappa \right)A_{-}{\text{e}}^{\kappa z_0}=\left(1-u_{+\mathrm{,1}}\left(z_0\right)-2u_{+\mathrm{,2}}\left(z_0\right)\kappa \right)A_{+}{\text{e}}^{-\kappa z_0}\mathrm{.}$ (8.11)

Our notations specialized to the case (8.1) are ( $m=0,1,2$ )

$u_m\left(z_0\right)\equiv u_{-\mathrm{,}m}\left(z_0\right)+u_{+\mathrm{,}m}\left(z_0\right)=\frac{{\left(-1\right)}^m}{m\mathrm{<mo>\; !\; </mo>}}{\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}z\left(u\left(z\right)-u_{\infty }\right){\left(z-z_0\right)}^m\mathrm{.}$ (8.12)

As a third-degree equation with real coefficients it may possess three real solutions or one real plus two complex-conjugate solutions where due to the smallness of the coefficients in front of ${\kappa }^3$ and of ${\kappa }^2$ the real solution (8.6) is only slightly modified and the new solutions are small from their modulus in comparison to (8.6).

9. Favorable Form of Boundary Conditions for Layer Systems

In case of a layer system with each single layer of the considered form it is favorable to represent the boundary conditions in a form where the amplitudes to each side of the layer stand on different sides that means

${\text{S}}_{-}\left(\begin{array}{c}{A}_{-}^{\text{i}}{\text{e}}^{+\text{i}{k}_{-}{z}_0}\\ A_{-}^{\text{r}}{\text{e}}^{-\text{i}{k}_{-}{z}_0}\end{array}\right)={\text{S}}_{+}\left(\begin{array}{c}{A}_{+}^{\text{r}}{\text{e}}^{+\text{i}{k}_{+}{z}_0}\\ A_{+}^{\text{i}}{\text{e}}^{-\text{i}{k}_{+}{z}_0}\end{array}\right)\mathrm{.}$ (9.1)

Then the output of one side of the transition layer may serve as input for the next transition layer which have to be in the asymptotically range because in other case we must add new boundary relations for the patch position of the two transition layers.

The resolution of (9.1) to the right-hand side of the amplitudes is

$\left(\begin{array}{c}{A}_{+}^{\text{r}}{\text{e}}^{+\text{i}{k}_{+}{z}_0}\\ A_{+}^{\text{i}}{\text{e}}^{-\text{i}{k}_{+}{z}_0}\end{array}\right)=\frac{\overline{{\text{S}}_{+}}{\text{S}}_{-}}{\left[{\text{S}}_{+}\right]}\left(\begin{array}{c}{A}_{-}^{\text{i}}{\text{e}}^{+\text{i}{k}_{-}{z}_0}\\ A_{-}^{\text{r}}{\text{e}}^{-\text{i}{k}_{-}{z}_0}\end{array}\right)\equiv \text{T}\left(\begin{array}{c}{A}_{-}^{\text{i}}{\text{e}}^{+\text{i}{k}_{-}{z}_0}\\ A_{-}^{\text{r}}{\text{e}}^{-\text{i}{k}_{-}{z}_0}\end{array}\right)\equiv \left(\begin{array}{cc}{T}^{\text{ri}}& T^{\text{rr}}\\ T^{\text{ii}}& T^{\text{ir}}\end{array}\right)\left(\begin{array}{c}{A}_{-}^{\text{i}}{\text{e}}^{+\text{i}{k}_{-}{z}_0}\\ A_{-}^{\text{r}}{\text{e}}^{-\text{i}{k}_{-}{z}_0}\end{array}\right)\mathrm{.}$ (9.2)

In addition, it would be favorable to substitute in case of many layers when their connection is in the asymptotic ranges the lower indices “−” and “+” by specific notations of the involved layers to both sides of the transitions. The

matrices $\text{T}\equiv \frac{\overline{{\text{S}}_{+}}{\text{S}}_{-}}{\left[{\text{S}}_{+}\right]}$ we may call transmission matrices.

In second-order approximation we obtain for the matrices ${\text{S}}_{+}$ and ${\text{S}}_{-}$ under the same propositions (5.2) as before

$\begin{array}{l}{\text{S}}_{-}=\left(\begin{array}{cc}{k}_{-}\left(1-u_{-\mathrm{,1}}\left(z_0\right)\right)-\text{i}{u}_{-\mathrm{,0}}\left(z_0\right)& -k_{-}\left(1-u_{-\mathrm{,1}}\left(z_0\right)\right)-\text{i}{u}_{-\mathrm{,0}}\left(z_0\right)\\ 1+u_{-\mathrm{,1}}\left(z_0\right)& 1+u_{-\mathrm{,1}}\left(z_0\right)\end{array}\right)\mathrm{,}\\ {\text{S}}_{+}=\left(\begin{array}{cc}{k}_{+}\left(1+u_{+\mathrm{,1}}\left(z_0\right)\right)+\text{i}{u}_{+\mathrm{,0}}\left(z_0\right)& -k_{+}\left(1+u_{+\mathrm{,1}}\left(z_0\right)\right)+\text{i}{u}_{+\mathrm{,0}}\left(z_0\right)\\ 1-u_{+\mathrm{,1}}\left(z_0\right)& 1-u_{+\mathrm{,1}}\left(z_0\right)\end{array}\right)\mathrm{.}\end{array}$ (9.3)

The determinant $\left[{\text{S}}_{+}\right]$ of the matrix ${\text{S}}_{+}$ is

$\left[{\text{S}}_{+}\right]=2\left(1-u_{+\mathrm{,1}}\left(z_0\right)\right)\left(1+u_{+\mathrm{,1}}\left(z_0\right)\right)k_{+}\mathrm{.}$ (9.4)

Except trivial cases this determinant and the determinant $\left[{\text{S}}_{-}\right]$ of the matrix ${\text{S}}_{-}$ too cannot vanish and solutions which vanish only on one side of the transition layer do not exist. For the components of the transmission matrices $\text{T}$ we calculate from (9.2)

$\begin{array}{l}\left[{\text{S}}_{+}\right]T^{\text{ri}}=\left(1+u_{+\mathrm{,1}}\left(z_0\right)u_{-\mathrm{,1}}\left(z_0\right)\right)\left(k_{+}+k_{-}\right)+u_1\left(z_0\right)\left(k_{+}-k_{-}\right)\\ -\text{i}\left(u_0\left(z_0\right)+u_{+\mathrm{,0}}\left(z_0\right)u_{-\mathrm{,1}}\left(z_0\right)-u_{-\mathrm{,0}}\left(z_0\right)u_{+\mathrm{,1}}\left(z_0\right)\right)\mathrm{,}\end{array}$

$\begin{array}{l}\left[{\text{S}}_{+}\right]T^{\text{rr}}=\left(1+u_{+\mathrm{,1}}\left(z_0\right)u_{-\mathrm{,1}}\left(z_0\right)\right)\left(k_{+}-k_{-}\right)+u_1\left(z_0\right)\left(k_{+}+k_{-}\right)\\ -\text{i}\left(u_0\left(z_0\right)+u_{+\mathrm{,0}}\left(z_0\right)u_{-\mathrm{,1}}\left(z_0\right)-u_{-\mathrm{,0}}\left(z_0\right)u_{+\mathrm{,1}}\left(z_0\right)\right)\mathrm{,}\end{array}$

$\begin{array}{l}\left[{\text{S}}_{+}\right]T^{\text{ii}}=\left(1+u_{+\mathrm{,1}}\left(z_0\right)u_{-\mathrm{,1}}\left(z_0\right)\right)\left(k_{+}-k_{-}\right)+u_1\left(z_0\right)\left(k_{+}+k_{-}\right)\\ +\text{i}\left(u_0\left(z_0\right)+u_{+\mathrm{,0}}\left(z_0\right)u_{-\mathrm{,1}}\left(z_0\right)-u_{-\mathrm{,0}}\left(z_0\right)u_{+\mathrm{,1}}\left(z_0\right)\right)\mathrm{,}\end{array}$

$\begin{array}{l}\left[{\text{S}}_{+}\right]T^{\text{ir}}=\left(1+u_{+\mathrm{,1}}\left(z_0\right)u_{-\mathrm{,1}}\left(z_0\right)\right)\left(k_{+}+k_{-}\right)+u_1\left(z_0\right)\left(k_{+}-k_{-}\right)\\ +\text{i}\left(u_0\left(z_0\right)+u_{+\mathrm{,0}}\left(z_0\right)u_{-\mathrm{,1}}\left(z_0\right)-u_{-\mathrm{,0}}\left(z_0\right)u_{+\mathrm{,1}}\left(z_0\right)\right)\mathrm{.}\end{array}$ (9.5)

In the first-order approximation we neglect the moments $u_{+\mathrm{,1}}\left(z_0\right)$ and $u_{-\mathrm{,1}}\left(z_0\right)$ that we do not write down separately.

10. Conclusions

In this article, we have chosen the important Sturm-Liouville equation as one of the simplest ordinary second-order differential equations to demonstrate important aspects of the derivation of generalized boundary conditions by means of Generalized functions or, more especially, of step functions, delta functions and their derivatives. A countable infinite system of boundary conditions has been obtained and its first few low-order approximations are derived for this equation in case of sufficiently smooth potentials. It was necessary to choose two asymptotic “volume” potentials with known solutions of the corresponding Sturm-Liouville equation and an imaginary’ boundary surface in the transition region between these asymptotic potentials. The information about the deviations of the real potential from the chosen asymptotic potentials on both sides of the “imaginary” boundary surface is involved in the boundary conditions in form of the moments of these deviations. The condition for the applicability of low-order approximations of the boundary conditions requires that the “effective” thickness of the transition region must be small in comparison to the characteristic lengths in alterations of the “volume” solutions in the transition region (wavelengths in the case of asymptotic constant potentials, for instance). Although this is intuitively clear, it is convenient to discuss this more precisely not in general abstract form but in connection with concrete potentials and problems. We applied it to reflection and refraction problems for asymptotically constant potentials.

The general problems of the derivation of boundary conditions by means of the theory of Generalized functions could be solved for the Sturm-Liouville equation. It seems that there are no principal difficulties to derive boundary conditions in an analogous way for other linear scalar or vectorial differential or integro-differential equations as was already done for electrodynamics. However, we do not want to conceal that in the case of nonlinear equations there are some new principal difficulties with the presented method.

Acknowledgements

Lately, I am able to express my great gratitude to Prof. G. Vojta (from Editorial Board of Ann. d. Physik in former GDR) who advocated the publication of my paper [12] against the rejection without a written report by the Chief Editor of Ann. d. Physik Prof. G. Richter from GDR who at that same time was also a Member of our Central Institute of Optics and Spectroscopy (ZOS). The arguments of G. Richter were unjustified and unacceptable (see my home page)⁵. Furthermore, in particular, I thank Dr. Th. Jung from “Zentralinstitut für Elektronenphysik” der “Akademie der Wissenschaften der DDR” who managed to give a colloquium lecture in mentioned Institute and also Prof. V.M. Agranovich from Moscow where during a one-month visit I could give a seminar lecture in his “Institute of Spectroscopy” nearby Moscow to the problems dealt with in [12]. In my former “Zentralinstitut für Optik und Spektroskopy” I never gave a lecture about these problems to avoid further trouble with some persons.

Appendix A. Multiplication of Delta Function and of Its Derivatives with Smooth Functions and Moment Series of Generalized Functions

We derive in this Appendix some identities for generalized functions which we use in the main text. Generalized functions $f\left(x\right)$ of a function space K' are defined as linear continuous functionals over basis functions $\phi \left(x\right)$ of a space K of (sufficiently well-behaved, e.g., smooth, i.e., continuous and infinitely continuously differentiable) functions which we may write as generalization of the integral on the right-hand side in the relation [25] [26] [27] [28]

$\left(f\mathrm{,}\phi \right)\equiv {\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}xf\left(x\right)\phi \left(x\right)\mathrm{,}\underset{k\to \infty }{lim}\left(f_k\mathrm{,}\phi \right)=\left(f\mathrm{,}\phi \right)\mathrm{.}$ (A.1)

The delta function with real displacement parameter $x_0$ and their n-th derivatives ${\delta }^{\left(n\right)}\left(x-x_0\right)$ according to the rule of partial integration are defined by

${\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}x{\delta }^{\left(n\right)}\left(x-x_0\right)\phi \left(x\right)={\left(-1\right)}^n{\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}x\delta \left(x-x_0\right){\phi }^{\left(n\right)}\left(x\right)={\left(-1\right)}^n{\phi }^{\left(n\right)}\left(x_0\right)\mathrm{,}$ (A.2)

with the important special cases for $\phi \left(x\right)=1$

${\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}x{\delta }^{\left(n\right)}\left(x-x_0\right)={\delta }_{n\mathrm{,0}}=\{\begin{array}{ll}\mathrm{1,}\hfill & n=0\hfill \\ \mathrm{0,}\hfill & n\ne 0\hfill \end{array}\mathrm{,}$ (A.3)

showing that all “integrals” over genuine (i.e., $n\ge 1$ ) derivatives of the delta function taken over the real axis vanish. With possible (e.g., smooth) multiplication functions $\alpha \left(x\right)$ using the Leibniz rule for differentiation of a product one finds

$\begin{array}{l}{\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}x\alpha \left(x\right){\delta }^{\left(n\right)}\left(x-x_0\right)\phi \left(x\right)\\ ={\left(-1\right)}^n{\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}x\delta \left(x-x_0\right){\left(\alpha \left(x\right)\phi \left(x\right)\right)}^{\left(n\right)}\\ ={\left(-1\right)}^n{\displaystyle \underset{k=0}{\overset{n}{\sum }}}\frac{n!}{k!\left(n-k\right)!}{\alpha }^{\left(k\right)}\left(x_0\right){\phi }^{\left(n-k\right)}\left(x_0\right)\\ ={\displaystyle \underset{k=0}{\overset{n}{\sum }}}\frac{{\left(-1\right)}^{k}n!}{k!\left(n-k\right)!}{\alpha }^{\left(k\right)}\left(x_0\right){\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}x{\delta }^{\left(n-k\right)}\left(x-x_0\right)\phi \left(x\right),\end{array}$ (A.4)

and since this is correct for arbitrary basis functions $\phi (\; x\; )$

$\alpha \left(x\right){\delta }^{\left(n\right)}\left(x-x_0\right)={\displaystyle \underset{k=0}{\overset{n}{\sum }}}\frac{{\left(-1\right)}^{k}n!}{k!\left(n-k\right)!}{\alpha }^{\left(k\right)}\left(x_0\right){\delta }^{\left(n-k\right)}\left(x-x_0\right).$ (A.5)

On the right-hand side of this identity we have only constant coefficients in front of the derivatives of the delta function $\delta \left(x-x_0\right)$.

We give two examples for (A.5). The first example with $\alpha \left(x\right)=x^m$ is

$\begin{array}{c}{x}^m{\delta }^{\left(n\right)}\left(x-x_0\right)={\displaystyle \underset{k=0}{\overset{\left\{m,n\right\}}{\sum }}}\frac{{\left(-1\right)}^{k}m!n!}{k!\left(m-k\right)!\left(n-k\right)!}{x}_0^{m-k}{\delta }^{\left(n-k\right)}\left(x-x_0\right)\\ ={\text{L}}_{m\mathrm{,}n}\left(x_0\mathrm{,}\frac{\partial }{\partial x}\right)\delta \left(x-x_0\right)\mathrm{,}\end{array}$ (A.6)

with the special case $x_0=0$

$x^m{\delta }^{\left(n\right)}\left(x\right)=\frac{{\left(-1\right)}^{m}n\mathrm{<mo>\; !\; </mo>}}{\left(n-m\right)\mathrm{<mo>\; !\; </mo>}}{\delta }^{\left(n-m\right)}\left(x\right)={\text{L}}_{m\mathrm{,}n}\left(\mathrm{0,}\frac{\partial }{\partial x}\right)\delta \left(x\right)\mathrm{.}$ (A.7)

In (A.6) ${\text{L}}_{m\mathrm{,}n}\left(x\mathrm{,}y\right)$ denotes the Laguerre 2D polynomials defined by

$\begin{array}{c}{\text{L}}_{m\mathrm{,}n}\left(x\mathrm{,}y\right)\equiv {\displaystyle \underset{k=0}{\overset{\left\{m\mathrm{,}n\right\}}{\sum }}}\frac{{\left(-1\right)}^{k}m\mathrm{<mo>\; !\; </mo>}n\mathrm{<mo>\; !\; </mo>}}{k\mathrm{<mo>\; !\; </mo>}\left(m-k\right)\mathrm{<mo>\; !\; </mo>}\left(n-k\right)\mathrm{<mo>\; !\; </mo>}}{x}^{m-k}{y}^{n-k}\\ ={\left(-1\right)}^{n}n!x^{m-n}{\text{L}}_n^{m-n}\left(xy\right)={\left(-1\right)}^{m}m!y^{n-m}{\text{L}}_m^{n-m}\left(xy\right)\mathrm{,}\end{array}$ (A.8)

where ${\text{L}}_n^{\nu }\left(z\right)$ denotes the Generalized Laguerre polynomials or Laguerre-Sonin polynomials in standard notation.

A second example for illustration with $\alpha \left(x\right)=\exp \left(-\frac{x^2}{a^2}\right)$ is

$\begin{array}{l}\exp \left(-\frac{x^2}{a^2}\right){\delta }^{\left(n\right)}\left(x-x_0\right)\\ =\exp \left(-\frac{x_0^2}{a^2}\right){\displaystyle \underset{k=0}{\overset{n}{\sum }}}\frac{n\mathrm{<mo>\; !\; </mo>}}{k\mathrm{<mo>\; !\; </mo>}\left(n-k\right)\mathrm{<mo>\; !\; </mo>}}\frac{1}{a^k}{\text{H}}_k\left(\frac{x_0}{a}\right){\delta }^{\left(n-k\right)}\left(x-x_0\right)\\ =\exp \left(-\frac{x_0^2}{a^2}\right)\frac{1}{a^n}{\text{H}}_n\left(\frac{x_0}{a}+\frac{a}{2}\frac{\partial }{\partial x}\right)\delta \left(x-x_0\right)\mathrm{,}\end{array}$ (A.9)

with the special case $x_0=0$

$\exp \left(-\frac{x^2}{a^2}\right){\delta }^{\left(n\right)}\left(x\right)={\displaystyle \underset{l=0}{\overset{\left[\frac{n}{2}\right]}{\sum }}}\frac{{\left(-1\right)}^{l}n\mathrm{<mo>\; !\; </mo>}}{l!\left(n-2l\right)!}\frac{1}{a^{2l}}{\delta }^{\left(n-2l\right)}\left(x\right)=\frac{1}{a^n}{\text{H}}_n\left(\frac{a}{2}\frac{\partial }{\partial x}\right)\delta \left(x\right)\mathrm{,}$ (A.10)

where ${\text{H}}_n\left(z\right)$ denotes the Hermite polynomials

$\begin{array}{c}{\text{H}}_n\left(z\right)\equiv {\left(-1\right)}^n\exp \left(z^2\right)\frac{{\partial }^n}{\partial z^n}\exp \left(-z^2\right)=\exp \left(-\frac{1}{4}\frac{{\partial }^2}{\partial z^2}\right){\left(2z\right)}^n\\ ={\displaystyle \underset{k=0}{\overset{\left[\frac{n}{2}\right]}{\sum }}}\frac{{\left(-1\right)}^{k}n\mathrm{<mo>\; !\; </mo>}}{k\mathrm{<mo>\; !\; </mo>}\left(n-2k\right)\mathrm{<mo>\; !\; </mo>}}{\left(2z\right)}^{n-2k}\mathrm{,}\end{array}$ (A.11)

in different representations which can be used as alternative definitions.

An equation of the form

$0={\displaystyle \underset{n=0}{\overset{\infty }{\sum }}}{c}_n\left(x_0\right){\delta }^{\left(n\right)}\left(x-x_0\right)=c_0\left(x_0\right)\delta \left(x-x_0\right)+c_1\left(x_0\right){\delta }^{\left(1\right)}\left(x-x_0\right)+\cdots \mathrm{,}$ (A.12)

with coefficients $c_n\left(x_0\right)$ which are independent of the variable x is only satisfied if all coefficients $c_n\left(x_0\right)$ vanish. This follows by multiplication of this equation with the functions ${\left(x-x_0\right)}^m\mathrm{,}\left(m=\mathrm{0,1,2,}\cdots \right)$ using (A.7) and succeeding integration over the whole x-axis using (A.3) according to

$\begin{array}{c}0={\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}x{\displaystyle \underset{n=0}{\overset{\infty }{\sum }}}{c}_n\left(x_0\right){\left(x-x_0\right)}^m{\delta }^{\left(n\right)}\left(x-x_0\right)\\ ={\displaystyle \underset{n=0}{\overset{\infty }{\sum }}}{c}_n\left(x_0\right)\frac{{\left(-1\right)}^{m}n!}{\left(n-m\right)!}{\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}x{\delta }^{\left(n-m\right)}\left(x-x_0\right)\\ ={\left(-1\right)}^{m}m!c_m\left(x_0\right),\Rightarrow c_m\left(x_0\right)=0,\left(m=0,1,2,\cdots \right).\end{array}$ (A.13)

It means that the systems of functions ${\left(x-x_0\right)}^m$ and generalized functions ${\delta }^n\left(x-x_0\right)$ with arbitrary parameter x₀ are bi-orthogonal. For the derivation of boundary conditions for an equation one has to split the asymptotic parts and has to transform the remaining parts of the equation to the principal form (A.12) and find then the boundary condition from the vanishing of the coefficients $c_m\mathrm{,}\left(m=\mathrm{0,1,2,}\cdots \right)$. From (A.13) follows with variation of the parameter x₀

$c_n\left(x_0\right)={\displaystyle \underset{k=0}{\overset{n}{\sum }}}\frac{x_0^k}{k!}{c}_{n-k}\left(0\right).$ (A.14)

This means that with variation of the parameter x₀ the whole system of boundary conditions transforms into each other. However, by truncation of the boundary conditions with finite numbers $n=N$ this full equivalence becomes an approximate equivalence and a “best” choice of the parameter x₀ becomes important.

We now generalize the preceding considerations to expansions of generalized functions which we call moment series because the coefficients of these expansions are determined by the moments of the function. For this purpose we first define the Fourier transform $\tilde{f}\left(k\right)$ of $f\left(x\right)$ which together with its inversion is given by

$\tilde{f}\left(u\right)={\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}yf\left(y\right)\exp \left(-\text{i}yu\right)\mathrm{,}\iff f\left(x\right)=\frac{1}{2\pi }{\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}u\tilde{f}\left(u\right)\exp \left(\text{i}ux\right)\mathrm{.}$ (A.15)

Then we represent the function $f\left(x\right)$ as convolution (notation “*”) of this function with the delta function and apply the displacement operator of the argument of a function as follows

$\begin{array}{c}f\left(x\right)=f\left(x\right)\ast \delta \left(x\right)={\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}yf\left(y\right)\delta \left(x-y\right)\\ ={\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}yf\left(y\right)\exp \left(-y\frac{\partial }{\partial x}\right)\delta \left(x\right)=\tilde{f}\left(-\text{i}\frac{\partial }{\partial x}\right)\delta \left(x\right)\mathrm{,}\end{array}$ (A.16)

where we applied in last step the definition (A.15) of the Fourier transform. Making now expansions of the function $\delta \left(x-y\right)$ in a Taylor series with respect to powers of y and $\tilde{f}\left(k\right)$ in a Taylor series at $k=0$ we find from (A.16)

$f\left(x\right)={\displaystyle \underset{n=0}{\overset{\infty }{\sum }}}{f}_n\left(0\right){\delta }^{\left(n\right)}\left(x\right)={\displaystyle \underset{n=0}{\overset{\infty }{\sum }}}\frac{{\tilde{f}}^{\left(n\right)}\left(0\right)}{n\mathrm{<mo>\; !\; </mo>}}{\left(-\text{i}\frac{\partial }{\partial x}\right)}^n\delta \left(x\right)\mathrm{,}$ (A.17)

with the moments $f_n\left(0\right)$ of the function $f\left(x\right)$ defined by

$f_n\left(0\right)\equiv \frac{{\left(-1\right)}^n}{n\mathrm{<mo>\; !\; </mo>}}{\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}yf\left(y\right)y^n=\frac{{\left(-\text{i}\right)}^n}{n\mathrm{<mo>\; !\; </mo>}}{\tilde{f}}^{\left(n\right)}\left(0\right)\mathrm{.}$ (A.18)

This shows that the moments of $f\left(x\right)$ are related to the derivatives of the Fourier transforms $\tilde{f}\left(u\right)$ at $x=0$. Conversely, the derivatives of the function $f\left(x\right)$ at $x_0=0$ are related to the moments of its Fourier transform $\tilde{f}\left(u\right)$ by

$f^{\left(n\right)}\left(0\right)=\frac{{\left(-\text{i}\right)}^n}{2\pi }n!{\tilde{f}}_n\left(0\right)\mathrm{,}$ (A.19)

with ${\tilde{f}}_n\left(0\right)$ defined by

${\tilde{f}}_n\left(0\right)\equiv \frac{{\left(-1\right)}^n}{n\mathrm{<mo>\; !\; </mo>}}{\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}u\tilde{f}\left(u\right)u^n=2\pi \frac{{\text{i}}^n}{n\mathrm{<mo>\; !\; </mo>}}{f}^{\left(n\right)}\left(0\right)\mathrm{.}$ (A.20)

The relations (A.18) and (A.19) express a certain duality of the moments of a function and of the derivatives of its Fourier transforms at the coordinate origin. Upper indices in the form $f^{\left(n\right)}\left(x\right)$ mean in every case the n-th derivative of $f\left(x\right)$ and lower indices in the form $f_n\left(x_0\right)$ moments of a function to reference point x₀⁷.

The convergence of moment series is weak convergence in the sense of the convergence of the sequence of partial sums as linear functionals over the space of basis functions such as written in (A.1). The linear continuous functionals over basis functions $\phi \left(x\right)$ can be transformed by Fourier transformation (A.15) in linear continuous functionals over the space of Fourier transforms to the basis functions of the space $\mathcal{D}$ usually denoted by $\mathcal{Z}$ as follows

$\begin{array}{c}\left(f\mathrm{,}\phi \right)\equiv {\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}xf\left(x\right)\phi \left(x\right)\\ ={\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}x\frac{1}{2\pi }{\displaystyle {\int }_{-\infty }^{+\infty }}\tilde{f}\left(u\right){\text{e}}^{\text{i}ux}\phi \left(x\right)\\ ={\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}u\tilde{f}\left(u\right)\frac{1}{2\pi }{\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}x\phi \left(x\right){\text{e}}^{\text{i}xu}\\ =\frac{1}{2\pi }{\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}u\tilde{f}\left(u\right)\tilde{\phi }\left(-u\right)\equiv \frac{1}{2\pi }\left(\tilde{f}\mathrm{,}{\tilde{\phi }}^{\mathrm{<mo>\; *\; </mo>}}\right)\mathrm{,}\end{array}$ (A.21)

where we changed the order of integrations. The obtained relation is the Parseval theorem.

The expansion (A.17) can be generalized to an arbitrary reference point x₀ instead of $x_0=0$ and instead of (A.16) we have

$\begin{array}{c}f\left(x\right)={\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}yf\left(y\right)\delta \left(x-x_0+x_0-y\right)\\ ={\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}yf\left(y\right)\exp \left(\left(x_0-y\right)\frac{\partial }{\partial x}\right)\delta \left(x-x_0\right)\\ =\exp \left(x_0\frac{\partial }{\partial x}\right)\tilde{f}\left(-\text{i}\frac{\partial }{\partial x}\right)\delta \left(x-x_0\right)\mathrm{.}\end{array}$ (A.22)

The operators $\exp \left(x_0\frac{\partial }{\partial x}\right)$ and $\tilde{f}\left(-\text{i}\frac{\partial }{\partial x}\right)$ commute and using this it is easy to return to the form (A.16) which also follows as special case setting $x_0=0$. From (A.22) follows

$f\left(x\right)={\displaystyle \underset{n=0}{\overset{\infty }{\sum }}}{f}_n\left(x_0\right){\delta }^{\left(n\right)}\left(x-x_0\right)\mathrm{,}$ (A.23)

with the moments $f_n\left(x_0\right)$ determined by

$\begin{array}{c}{f}_n\left(x_0\right)\equiv \frac{{\left(-1\right)}^n}{n\mathrm{<mo>\; !\; </mo>}}{\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}yf\left(y\right){\left(y-x_0\right)}^n={\displaystyle \underset{k=0}{\overset{n}{\sum }}}\frac{x_0^k}{k\mathrm{<mo>\; !\; </mo>}}{f}_{n-k}\left(0\right)\\ ={\displaystyle \underset{k=0}{\overset{n}{\sum }}}\frac{{\left(-\text{i}\right)}^{n-k}}{k\mathrm{<mo>\; !\; </mo>}\left(n-k\right)\mathrm{<mo>\; !\; </mo>}}{x}_0^k{\tilde{f}}^{\left(n-k\right)}\left(0\right)=\frac{1}{n\mathrm{<mo>\; !\; </mo>}}{\left\{{\left(x_0-\text{i}\frac{\partial }{\partial u}\right)}^n\tilde{f}\left(u\right)\right\}}_{u=0}\mathrm{,}\end{array}$ (A.24)

where the relation of the moments at reference point $x=x_0$ to the moments at reference point $x=0$ is given. The zeroth moment $f_0\left(x_0\right)$ does not depend on the reference point x₀

$f_0\left(x_0\right)=f_0\left(0\right)=\tilde{f}\left(0\right)={\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}yf\left(y\right),$ (A.25)

and, therefore, is “absolute”. Furthermore, this formula shows how the moments to the reference points x₀ of a function can be expressed in simple way by the Fourier transform of this function.

For example, for a “normalized” Gaussian function together with its Fourier transform given by

$f\left(x\right)=\frac{1}{\sqrt{\pi a^2}}\exp \left(-\frac{x^2}{a^2}\right),\iff \tilde{f}\left(u\right)=\exp \left(-\frac{a^2{u}^2}{4}\right),$ (A.26)

we find the following moment series

$\frac{1}{\sqrt{\pi a^2}}\exp \left(-\frac{x^2}{a^2}\right)=\exp \left(\frac{a^2}{4}\frac{{\partial }^2}{\partial x^2}\right)\delta \left(x\right)={\displaystyle \underset{m=0}{\overset{\infty }{\sum }}}\frac{a^{2m}}{m!2^{2m}}{\delta }^{\left(2m\right)}\left(x\right),$ (A.27)

with the moments $f_n\left(0\right)$ determined by

$f_{2m}\left(0\right)=\frac{a^{2m}}{m!2^{2m}},f_{2m+1}\left(0\right)=0,\left(m=0,1,2,\cdots \right).$ (A.28)

Due to symmetry of the Gaussian function the odd moments are vanishing with respect to the reference point $x_0=0$ but not with respect to an arbitrary reference point x₀. Further examples for moment series of symmetric functions one may find in Appendix B.

As special cases we consider now the moment series for functions which are non-vanishing only to the left or to right of a reference point $x=x_0$.

For functions $f_{-}\left(x\right)$ which are non-vanishing only for variable $x<x_0$

$f_{-}\left(x\right)\equiv f\left(x\right)\theta \left(x_0-x\right)=\{\begin{array}{ll}f\left(x\right),\hfill & x<x_0\hfill \\ 0,\hfill & x>x_0\hfill \end{array},$ (A.29)

we find the moment series to reference point $x=x_0$

$f_{-}\left(x\right)={\displaystyle \underset{n=0}{\overset{\infty }{\sum }}}{f}_{-,n}\left(x_0\right){\delta }^{\left(n\right)}\left(x-x_0\right),f_{-}\left(x_0+y\right)={\displaystyle \underset{n=0}{\overset{\infty }{\sum }}}{f}_{-,n}\left(x_0\right){\delta }^{\left(n\right)}\left(y\right),$ (A.30)

with the moments

$\begin{array}{c}{f}_{-,n}\left(x_0\right)\equiv \frac{{\left(-1\right)}^n}{n!}{\displaystyle {\int }_{-\infty }^0}\text{d}yf\left(x_0+y\right)y^n=\frac{1}{n!}{\displaystyle {\int }_{-\infty }^{x_0}}\text{d}xf\left(x\right)\underset{\ge 0}{\underset{︸}{{\left(x_0-x\right)}^n}}\\ ={\displaystyle {\int }_{-\infty }^{x_0}}\text{d}{x}_n{\displaystyle {\int }_{-\infty }^{x_n}}\text{d}{x}_{n-1}\cdots {\displaystyle {\int }_{-\infty }^{x_2}}\text{d}{x}_1{\displaystyle {\int }_{-\infty }^{x_1}}\text{d}xf\left(x\right),\end{array}$ (A.31)

given together with an equivalent multiple integral.

For functions $f_{+}\left(x\right)$ which are non-vanishing only for variable $x>x_0$

$f_{+}\left(x\right)\equiv f\left(x\right)\theta \left(x-x_0\right)=\{\begin{array}{ll}f\left(x\right),\hfill & x>x_0\hfill \\ 0,\hfill & x<x_0\hfill \end{array},$ (A.32)

we find the moment series to reference point $x=x_0$

$f_{+}\left(x\right)={\displaystyle \underset{n=0}{\overset{\infty }{\sum }}}{f}_{+,n}\left(x_0\right){\delta }^{\left(n\right)}\left(x-x_0\right),f_{+}\left(x_0+y\right)={\displaystyle \underset{n=0}{\overset{\infty }{\sum }}}{f}_{+,n}\left(x_0\right){\delta }^{\left(n\right)}\left(y\right),$ (A.33)

with the moments

$\begin{array}{c}{f}_{+,n}\left(x_0\right)\equiv \frac{{\left(-1\right)}^n}{n!}{\displaystyle {\int }_0^{+\infty }}\text{d}yf\left(x_0+y\right)y^n=\frac{{\left(-1\right)}^n}{n!}{\displaystyle {\int }_{x_0}^{+\infty }}\text{d}xf\left(x\right)\underset{\ge 0}{\underset{︸}{{\left(x-x_0\right)}^n}}\\ ={\left(-1\right)}^n{\displaystyle {\int }_{x_0}^{+\infty }}\text{d}{x}_n{\displaystyle {\int }_{x_n}^{+\infty }}\text{d}{x}_{n-1}\cdots {\displaystyle {\int }_{x_2}^{+\infty }}\text{d}{x}_1{\displaystyle {\int }_{x_1}^{+\infty }}\text{d}xf\left(x\right),\end{array}$ (A.34)

given again together with an equivalent multiple integral.

The sum of the moments $f_{-,n}\left(x_0\right)$ and $f_{+,n}\left(x_0\right)$ is

$f\left(x\right)=f_{-}\left(x\right)+f_{+}\left(x\right),f_{-,n}\left(x_0\right)+f_{+,n}\left(x_0\right)=f_n\left(x_0\right),$ (A.35)

to the moment series (A.23) with the moments determined by (A.24).

It should be mentioned here that the identities in the formulae (A.31) and (A.34) are special cases of the more general identity

$\begin{array}{c}{f}^{\left(-\left(n+1\right)\right)}\left(a;y\right)\equiv {\displaystyle {\int }_a^y}\text{d}{x}_n{\displaystyle {\int }_a^{x_n}}\text{d}{x}_{n-1}\cdots {\displaystyle {\int }_a^{x_2}}\text{d}{x}_1{\displaystyle {\int }_a^{x_1}}\text{d}xf\left(x\right)\\ =\frac{1}{n!}{\displaystyle {\int }_a^y}\text{d}xf\left(x\right){\left(y-x\right)}^n,\end{array}$ (A.36)

with arbitrary parameter a that can be proved by complete induction with reordering of the arising double integral and accomplishing then the inner integral. With the choice $a=0$ and $n\to \mu -1$ as an arbitrary real number it is often taken as starting point for the introduction of “ $\mu$ -fold” fractional integration by the convolution integral on the right-hand side which in our written case corresponds to possible $\left(n+1\right)$ -fold integration of the function $f\left(x\right)$ (e.g., [29], Chap. XIII).

The moment series (A.23) can be made also for functions $f\left(x\right)$ which are not smooth and even not continuous and, roughly speaking, they are well defined if only all moments exist. The approximate reconstruction of the function $f\left(x\right)$ from the moment series can be obtained from the limits of the linear functionals $\left(f_k\mathrm{,}\phi \right)$ when $f_k$ denotes the truncated moment series up to the term with the k-th moment (see (A.1)). As basis functions $\phi \left(x\right)$ it seems to be sufficient to take, for example, all displaced normalized Gaussian functions

$\frac{1}{\sqrt{\pi {\epsilon }^2}}\exp \left(-\frac{{\left(x-x^{\prime }\right)}^2}{{\epsilon }^2}\right)$ with the limit $\epsilon \to 0$. However also the Fourier transforms $\tilde{f}\left(u\right)$ are appropriate if their Taylor series may be summed up and if then the inverse Fourier transformation $f\left(x\right)$ to $\tilde{f}\left(u\right)$ is calculable.

Appendix B. Examples of Potential Barriers and Potential Wells and of Transition Functions with Their Characteristics by Generalized Functions

In this Appendix we consider for illustration some examples of functions $f\left(x\right)$ which are appropriate as potential barriers or potential wells together with their Fourier transforms $\tilde{f}\left(u\right)$ (notations as in (A.15)) and their moment series and, furthermore, functions which are appropriate to simulate transitions between two potentials for which the Sturm-Liouville equation can be exactly solved (e.g., constant potentials).

In FigureA1 we present six from seven explicitly given symmetric functions $f_{\nu }\left(x\right)\mathrm{,}\left(\nu =\mathrm{1,}\cdots \mathrm{,7}\right)$ which are appropriate as potential barriers and, negatively taken, as potential wells

(1): $f_1\left(x\right)=\theta \left(1-\frac{2\left|x\right|}{a}\right)\mathrm{,}{\tilde{f}}_1\left(u\right)=\frac{2}{u}\sin \left(\frac{a}{2}u\right)\mathrm{,}$

(2): $f_2\left(x\right)=\sqrt{1-{\left(\frac{\pi x}{2a}\right)}^2}\theta \left(1-\frac{\pi \left|x\right|}{2a}\right)\mathrm{,}{\tilde{f}}_2\left(u\right)=\frac{\pi }{u}{\text{J}}_1\left(\frac{2a}{\pi }u\right)\mathrm{,}$

(3): $f_3\left(x\right)=\left(1-\frac{\left|x\right|}{a}\right)\theta \left(1-\frac{\left|x\right|}{a}\right)\mathrm{,}{\tilde{f}}_3\left(u\right)=a{\left(\frac{\sin \left(\frac{a}{2}u\right)}{\frac{a}{2}u}\right)}^2\mathrm{.}$

(4): $f_4\left(x\right)=\exp \left(-\frac{2\left|x\right|}{a}\right)\mathrm{,}{\tilde{f}}_4\left(u\right)=\frac{a}{1+{\left(\frac{a}{2}u\right)}^2}\mathrm{,}$

Figure A1. Examples for symmetric potential barriers and potential wells. The potential functions (5) and (6) and also (6) and (7) are very similar in some parts within the shown interval and it is difficult to distinguish them in the figure and therefore (7) is not drawn. The potential wells belong to the negatively taken functions $f\left(x\right)$.

(5): $f_5\left(x\right)=\exp \left(-\frac{\pi x^2}{a^2}\right)\mathrm{,}{\tilde{f}}_5\left(u\right)=a\exp \left(-\frac{a^2}{4\pi }{u}^2\right)\mathrm{,}$

(6): $f_6\left(x\right)=\frac{1}{{\text{ch}}^2\left(\frac{2x}{a}\right)}\mathrm{,}{\tilde{f}}_6\left(u\right)=a\frac{\pi au}{4\text{sh}\left(\frac{\pi a}{4}u\right)}\mathrm{,}$

(7): $f_7\left(x\right)=\frac{1}{\text{ch}\left(\frac{\pi x}{a}\right)}\mathrm{,}{\tilde{f}}_7\left(u\right)=\frac{a}{\text{ch}\left(\frac{a}{2}u\right)}\mathrm{.}$ (B.1)

Their first three initial terms of the moment series are

(1): $f_1\left(x\right)=a\left\{\delta \left(x\right)+\frac{a^2}{24}{\delta }^{\left(2\right)}\left(x\right)+\frac{a^4}{1920}{\delta }^{\left(4\right)}\left(x\right)+\cdots \right\}\mathrm{,}$

(2): $f_2\left(x\right)=a\left\{\delta \left(x\right)+\frac{a^2}{2{\pi }^2}{\delta }^{\left(2\right)}\left(x\right)+\frac{a^4}{12{\pi }^4}{\delta }^{\left(4\right)}\left(x\right)+\cdots \right\}\mathrm{,}$

(3): $f_3\left(x\right)=a\left\{\delta \left(x\right)+\frac{a^2}{12}{\delta }^{\left(2\right)}\left(x\right)+\frac{a^4}{360}{\delta }^{\left(4\right)}\left(x\right)+\cdots \right\}\mathrm{,}$

(4): $f_4\left(x\right)=a\left\{\delta \left(x\right)+\frac{a^2}{4}{\delta }^{\left(2\right)}\left(x\right)+\frac{a^4}{16}{\delta }^{\left(4\right)}\left(x\right)+\frac{a^6}{64}{\delta }^{\left(6\right)}\left(x\right)+\cdots \right\}\mathrm{,}$

(5): $f_5\left(x\right)=a\left\{\delta \left(x\right)+\frac{a^2}{4\pi }{\delta }^{\left(2\right)}\left(x\right)+\frac{a^4}{32{\pi }^2}{\delta }^{\left(4\right)}\left(x\right)+\cdots \right\}\mathrm{,}$

(6): $f_6\left(x\right)=a\left\{\delta \left(x\right)+\frac{{\pi }^2{a}^2}{96}{\delta }^{\left(2\right)}\left(x\right)+\frac{7{\pi }^4{a}^4}{92160}{\delta }^{\left(4\right)}\left(x\right)+\cdots \right\}\mathrm{,}$

(7): $f_7\left(x\right)=a\left\{\delta \left(x\right)+\frac{a^2}{8}{\delta }^{\left(2\right)}\left(x\right)+\frac{5a^4}{384}{\delta }^{\left(4\right)}\left(x\right)+\cdots \right\}\mathrm{.}$ (B.2)

All considered examples of functions $f_{\nu }\left(x\right)$ possess analytic Fourier transforms ${\tilde{f}}_{\nu }\left(u\right)$ and, therefore, well-defined moment series.

We now consider transition functions. Examples for transition functions can be constructed from the potential functions $f_{\nu }\left(x\right)$ in (B.2) by integration in the following way

$F_{\nu }\left(x\right)=\frac{1}{a}{\displaystyle {\int }_0^x}\text{d}y{f}_{\nu }\left(y\right)\mathrm{,}{F}_{\nu }^{\left(1\right)}\left(x\right)=\frac{f_{\nu }\left(x\right)}{a}\mathrm{,}{F}_{\nu }\left(\mp \infty \right)=\mp \frac{1}{2}\mathrm{.}$ (B.3)

Since we have chosen ${\displaystyle {\int }_{-\infty }^{+\infty }}\text{d}xf\left(x\right)=a$ as the zeroth moment of the function $f\left(x\right)$ we obtain by (B.3) a transition step of the amount $F\left(+\infty \right)-F\left(-\infty \right)=1$ that means a unit step and an increase in the center $x=0$ of $F^{\left(1\right)}\left(0\right)=\frac{1}{a}$. For the functions in (B.2) these are the functions

(1): $F_1\left(x\right)=\{\begin{array}{ll}-\frac{1}{2},\hfill & x\le -\frac{a}{2}\hfill \\ \frac{x}{a},\hfill & -\frac{a}{2}<x<+\frac{a}{2}\hfill \\ +\frac{1}{2},\hfill & +\frac{a}{2}<x\hfill \end{array},$

(2): $F_2\left(x\right)=\{\begin{array}{ll}-\frac{1}{2},\hfill & x\le -\frac{2a}{\pi }\hfill \\ \frac{x}{2a}\sqrt{1-{\left(\frac{\pi x}{2a}\right)}^2}+\frac{1}{\pi }\arcsin \left(\frac{\pi x}{2a}\right),\hfill & -\frac{2a}{\pi }<x<+\frac{2a}{\pi }\hfill \\ +\frac{1}{2},\hfill & +\frac{2a}{\pi }<x\hfill \end{array},$

(3): $F_3\left(x\right)=\{\begin{array}{ll}-\frac{1}{2},\hfill & x\le -a\hfill \\ \frac{x}{a}\left(1-\frac{\left|x\right|}{2a}\right),\hfill & -a<x<+a\hfill \\ +\frac{1}{2},\hfill & +a<x\hfill \end{array},$

(4): $F_4\left(x\right)=\frac{1}{2}\left(1-\exp \left(-\frac{2\left|x\right|}{a}\right)\right)\frac{x}{\left|x\right|}\mathrm{,}$

(5): $F_5\left(x\right)=\frac{1}{2}\text{erf}\left(\frac{\pi x}{2a}\right)\mathrm{,}$

(6): $F_6\left(x\right)=\frac{1}{2}\text{th}\left(\frac{2x}{a}\right)\mathrm{,}$

Figure A2. Examples for antisymmetric transition functions. The transition step is here normalized to unity $F\left(+\infty \right)-F\left(-\infty \right)=1$ and the increase in the center of the transition is equal to $F^{\left(1\right)}\left(0\right)=\frac{1}{a}$. Some of the transition functions, for example, (5) and (6) are very similar in some parts and are hardly to distinguish in this figure and (7) is also very near to (5) and (6) and therefore not drawn.

(7): $F_7\left(x\right)=\frac{2}{\pi }\text{arctg}\left(\text{th}\left(\frac{\pi x}{2a}\right)\right)\mathrm{,}$ (B.4)

shown in FigureA2. In the first approximation of the boundary conditions all these transition functions are equivalent and only in higher approximations they provide differences in the results.

NOTES

¹Due to non-commutativity of the operators $\left(Q\mathrm{,}P\right)$ a Hamilton operator H by writing it $H\left(Q\mathrm{,}P\right)$ is not uniquely defined. Fortunately in history of quantum mechanics was that in its simplest form $H=T+U$ two parts which depend on P and Q alone were additively separated and the ordering problem of these operators was not stringent.

²According to the practice of that time de Broglie used as notations for Planck’s action quantum $h\equiv 2\pi \hslash$ and instead of the circular frequency $\omega$ the usual frequency $\nu \equiv \frac{\omega }{2\pi }$ that means $E=h\nu =\hslash \omega$.

³WKB stands for “Wentzel, Kramers, Brillouin”.

⁴Coordinate-invariant calculations in 2D-case usually do not bring a great advantage in comparison to coordinate calculations due to the relatively simple formulae

⁵A later paper with application of the derived boundary condition from my published paper [12] onto reflection and refraction at the boundary between isotropic and optically uniaxial media was again rejected by G. Richter with the same unacceptable arguments and is not published.

⁶I overtook my exemplar (a translation into Russian from Publishing House ‘ONTI Leningrad’, 1935) from my former colleague Randolf Fischer after his death from cancer in 1988 in the age of about 53. The history of this exemplar before this I do not know. Inside I find handwritten “Iz knig B. Popova” (from the books of B. Popov). I do not know whether R. Fischer bought this book in the former Sovyet Union as second-hand book or did he directly overtake it from a Russian colleague. Likely, it is my oldest book in Russian.

⁷In comparison to my paper [12] from 1980 I changed here slightly the principles of index notations.

Conflicts of Interest

The author declares no conflicts of interest regarding the publication of this paper.

References

[1]	Drude, P. (1889) Ueber Oberflächenschichten. Wiedemann’s Annalen der Physik und Chemie, 36, I. Theil 532-560, II. Theil 865-897.
[2]	Drude, P. (1912) Lehrbuch der Optik [Dritte Erweiterte Auflage]. Hirzel, Leipzig.6
[3]	Mayer, H. (1950) Physik dünner Schichten. Wissenschaftliche Verlags-Gesellschaft, Stuttgart.
[4]	Heavens, O.S. (1955) Optical Properties of Thin Films. Butterworths, London.
[5]	Rozenberg, G.B. (1958) Optika tonkosloinykh pokrytiy [Optics of Thin-Layered Coverings]. Gosfizmatizdat, Moskva.
[6]	Vašiček, A. (1969) Optics of Thin Films. North-Holland, Amsterdam.
[7]	Agranovich, V.M. and Ginzburg, V.L. (1979) Kristallooptika s uchotom prostranstvennoi dispersii i teoriya eksitonov [Crystal Optics with Spatial Dispersion]. Nauka, Moskva.
[8]	Bedeaux, D. and Vlieger, J. (1973) A Phenomenological Theory of the Dielectric Properties of Thin Films. Physica, 67, 55-73. https://doi.org/10.1016/0031-8914(73)90022-0
[9]	Bedeaux, D. and Vlieger, J. (1974) A Statistical Theory of the Dielectric Properties of Thin Island Films: I. The Surface Material Coefficients. Physica, 73, 287-311. https://doi.org/10.1016/0031-8914(74)90002-0
[10]	Timashev, S.F. and Krykin, A.M. (1976) Optical and Electrooptical Properties of the Interfaces of Media. Physica Status Solidi (b), 76, 67-76. https://doi.org/10.1002/pssb.2220760106
[11]	Timashev, S.F. and Krykin, A.M. (1979) Appendix to Osnovy ellipsometrii [Foundations of Ellipsometry]. Nauka, Novosibirsk.
[12]	Wünsche, A. (1980) Grenzbedingungen der Elektrodynamik für Medien mit räumlicher Dispersion und Übergangsschichten. Annalen der Physik, 37, 121-142. https://doi.org/10.1002/andp.19804920205
[13]	Lekner, J. (2016) Theory of Reflection. Springer, Heidelberg.
[14]	Landau, L.D. and Lifshits, Ye.M. (1982) Elektrodinamika sploshnykh sred [Electrodynamics of Continuous Media]. Nauka, Moskva.
[15]	Levitan, B.M. and Sargsyan, I.S. (1970) Vvedenye v spektralnuyu teoriyu [Introduction to Spectral Theory]. Nauka, Moskva.
[16]	Levitan, B.M. and Sargsyan, I.S. (1988) Operatory Shturma-Liuvilly [Sturm-Liouville Operators). Nauka, Moskva.
[17]	Marchenko, V.A. (1977) Operatory Shturma-Liuvilla i yikh prilozhenya [Sturm-Liouville Operators and Their Application]. Naukova dumka, Kiev.
[18]	Birman, M.S. and Faddeyev, L.D. (1972) Operatory kvantovoy mekhaniki [Operators of Quantum Mechanics]. Nauka, Moskva .
[19]	Krein, S.G. (1972) Funktsional’nyi analiz [Functional Analysis]. 2nd Edition, Nauka, Moskva.
[20]	Landau, L.D. and Lifshits, Ye.M. (1974) Kvantovaya mekhanika [Quantum Mechanics]. Nauka, Moskva.
[21]	Morse, P.M. and Feshbach, H. (1953) Methods of Theoretical Physics, Part II. McGraw Hill, New York.
[22]	Flügge, S. (1999) Practical Quantum Mechanics. Springer-Verlag, Berlin.
[23]	Brekhovskikh, L.M. (1973) Volny v sloistyhk sredakh [Waves in Layered Media]. Nauka, Moskva.
[24]	de Broglie, L. (1925) Recherches sur la théorie des Quanta. Annales de Physique, 10, 22-128. https://doi.org/10.1051/anphys/192510030022
[25]	Gel’fand, I.M. and Shilov, G.Ye. (1959) Obobstchonnyje funktsii, I [Generalized Functions]. Vol. 1, Fizmatgiz, Moskva.
[26]	Schwartz, L. (1961) Méthods Mathematiques pour les Science Physiques. Hermann, Paris.
[27]	Vladimirov, V.S. (1971) Uravnyeniya matematicheskoi fiziki [Equations of Mathematical Physics]. Nauka, Moskva.
[28]	Vilyenkin, N.Ya. (1972) Obobstchonnyje funktsii [Generalized Functions]. Nauka, Moskva.
[29]	Bateman, H. and Erdélyi, A. (1954) Tables of Integral Transforms, Vol. II. McGraw-Hill, New York.