## Abstract

The chapter illustrates how simple quantum mechanics can sometimes provide quite precise description of nanophysics phenomena. From this perspective, both exact and approximate solutions for the bound-state energy of an electron in a square well are exposed. These results are used to improve the calculation of quantum size effects (QSEs) in ultrathin metallic films, obtained by several authors with simpler models of quantum wells. We show that, for a small (less than 5) number of monolayers, the differences between the predictions of these simpler models, and our approach, are important. Methods to improve the accuracy in the evaluation of various quantum size effects are shortly discussed. Using quantum mechanical-electromagnetic analogies, our results can be used in the study of light propagation in dielectric wave guides.

### Keywords

- ultrathin metallic films
- quantum wells
- finite square well
- quantum size effects
- heterojunctions

## 1. Introduction

If the dimension of a physical system is reduced, on one or several directions, up to the nanometric range, electron confinement generates states specific to quantum wells, quantum dots, or to other nanostructures, studied by new branches of science, nanophysics and nanotechnology, extremely interesting from both applicative and fundamental perspectives. A fascinating aspect of nanophysics is that it can be sometimes understood using elementary, one-particle quantum mechanics; for instance, many phenomena specific to quantum wells can be treated using the simple model of a particle in a rectangular potential.

A class of physical systems which can be studied in the frame of this model is a particular kind of quantum wells—the ultrathin metallic films. In the last decades, they were object for active investigation, mainly due to the quantum size effects (QSEs), like the oscillatory behavior of the film stability [1], of the lattice deformation [2], of the work function [3], etc., in dependence of the number of atomic monolayers. The QSEs, predicted in the pioneering papers of Sandomirskii [4] and Schulte [5], are important for both practical and theoretical reasons. The ultrathin metallic films have a special relevance for ferromagnetic materials, as they are responsible for the giant magnetoresistivity of the Fe/Cr antiferromagnetic lattice [6]. Also, the possibility of obtaining ultrathin metallic films, having a specific number of monolayers, allows the experimentalist to tune the work function, controlling the chemistry of the metallic surface [3]. All these effects can be satisfactorily explained with a quite simple physics, whose basic ingredient is the different quantization imposed to electrons moving on longitudinal and transversal directions. Namely, the electrons moving parallel to the surface of the metallic film are quantized by cyclic conditions; the result is that the wave vectors are quasi-continuous. The electrons moving perpendicular to the film are considered as confined in a rectangular well, so they are quantized according to the theory of quantum wells; the result is that the spectrum is discrete.

How simple can the model of the well be, in order to provide a quantitative understanding of the physics of ultrathin metallic films? In spite of its simplicity, even the model of the infinite rectangular well gives sometimes good results, for instance, for the calculation of lattice deformation [2] or of Fermi energy [7]. These successes can be explained by the fact that, if the number of monolayers is not very small,

In this chapter, we shall present exact or approximate analytic results for the energy levels of a finite square well and show how they can improve the simple theoretical models which give a quantitative understanding of the behavior of ultrathin metallic films, especially the QSEs. Its structure is the following: in the second section, we shall discuss the quantum problem of the finite square well, mainly in order to put the eigenvalue equations in an appropriate form. The next one is a short review of the various attempts of solving these transcendental eigenvalue equations. The fourth section describes a simple algebraic approximation of the solution of the eigenvalue equations—the parabolic approximation—mentioning also similar but more precise approaches. In the next one, we put the eigenvalue equations in differential form and obtain the exact solution as a series expansion. The sixth section is devoted to the applications in the quantum statistical physics of the ultrathin metallic films of the analytic results obtained for the bound-state energy in a finite square well. By analyzing the predictions of the three models frequently used in the physics of ultrathin metallic films (infinite, semi-infinite, and finite square well) for the Fermi wave vector, we show the key role played by the finitude of the well, in the evaluation of QSEs. In the last section, we describe how our results can improve the current theory of this class of metallic films.

## 2. The bound states of a particle in a finite rectangular well

Until the mid-1980 of the previous century, the finite rectangular well was just an elementary problem of quantum mechanics, with applications in finding the energy levels of the quasi-free electrons on long molecules [8] or of the Ramsauer-Townsend effect [9]. The progress of solid-state physics, which finally led to the fabrication of quantum wells [10], quantum dots, or ultrathin metallic films [11, 12] and to the observation of QSEs associated with them, transformed these simple systems from problems of elementary quantum mechanics into theoretical models of devices of great practical interest.

We shall study now the movement of a particle in a finite rectangular well. There are, in principle, two ways of defining the potential of the well, choosing the origin of the energy

where

As the potential is invariant at spatial inversion,

where the quantities

For the particle outside the well, it is

The even solutions are

and the odd ones are

The continuity of the derivative in

and for odd states

Defining the dimensionless parameter

sometimes called potential strength, which actually characterizes both the particle

Also, the energy is

Using well-known trigonometric identities, Eqs. (11) and (12) take the form:

The sign must be chosen in agreement with Eqs. (11) and (12), so to satisfy the conditions

In other words, to solve the eigenvalue, Eq. (14) means to find the functions

This is, of course, a difficult task. If we write Eq. (15) in a slightly different form

to solve Eq. (15) means to invert the function

The functions

We shall write in a more explicit form Eq. (15), taking into account both the sign of the

where

and so on.

Similarly, the extremum points of the function

where

and so on. Each of Eqs. (18)–(20) and (22)–(23) has a unique solution,

According to Eq. (13), the energy eigenvalues are

If the particle moves not in potential

then the energy levels will be given by

According to the parity of

As already mentioned, the advantage of using the potential (1) is that the energy of a particle “inside the well,” so in a bound state, is negative, corresponding to the most usual convention of quantum mechanics. However, the form (26) of the potential has the advantage that its levels approach, in the limit of a very deep well, the levels of the infinite well. Indeed, for

and Eq. (27) gives the expression of the wave vector corresponding to the

## 3. Solving the eigenvalue equation of the finite well: a historical perspective

The eigenvalue equations for the wave vectors (18)–(20) and (22)–(24) are transcendental equations, and their solutions cannot be written as a finite combination of elementary functions. More than this, till now, they cannot be expressed neither in terms of the special functions of the mathematical physics. There are a large number of papers devoted to this subject, in the last 60 years.

The first one, due to Pitkanen [14], writes the eigenvalue Eqs. (8) and (9) in the simpler form (18)–(20) and (22)–(24), providing an interesting visualization of the solutions. The second one, due to Cantrell [15] (who does not cite [14], producing a delay in the circulation of this paper), also proposes the replacement of Eqs. (8) and (9) with (18)–(20) and (22)–(24)—in fact, a repetition of Pitkanen’s contribution—and notices that the eigenvalue equation for odd states is also the eigenvalue equation for a particle moving in a semi-infinite well, i.e., in a potential given by

Graphical solutions are proposed by Guest [16], who made visible the similarities between the bound-state energies in a finite well and the modes of a metallic wave guide ([17]; fig. (8.14)); actually, both the electrodynamic and quantum mechanical problems are equivalent forms of the same Sturm-Liouville problem [18]. Aronstein and Stroud [19] wrote the eigenvalue equation as

This elegant form had been already given in the first edition of Landau’s textbook of quantum mechanics, in the late 1940s of the twentieth century (for the English version of a more recent edition, see [20]) but remained unknown to Western physicists—a minor but significant consequence of the poor circulation of scientific information during the Cold War.

A completely different approach was proposed by Siewert [21], who obtained an exact solution in an integral form; unfortunately, it is very complicated and of limited practical use. Recently, Siewert’s solutions were discussed in the context of generalized Lambert functions [22], a subject under intense investigation.

Among the papers which provide approximate analytical solutions of the eigenvalue Eqs. (18)–(20) and (22)–(24), the most popular one, authored by Barker et al. [23], is essentially a low-order algebraic approximation of

as the magnitude of the domain outside the well, where the wave function can penetrate significantly, decreasing however exponentially. This concept is similar to the concept of skin depth in electromagnetism [17] or to the concept of viscous penetration depth in fluids [25],

In the context of various approximations, it is worth to mention the “algebraization” of trigonometric functions, proposed by de Alcantara Bonfim and Griffiths [26], which transforms the transcendental equations for the eigenvalues of the finite well in approximate, tractable, algebraic equations. For instance, we can use the approximations:

where the pair of constants can be chosen as

## 4. The parabolic approximation

To solve the eigenvalue equations, or—more generally—Eq. (16), with

Clearly, only the monotonic functions can be inverted; for instance, in our case, the function

Solving this equation for

and making the change

Following exactly the same steps, similar expressions can be obtained for

The method cannot be applied, evidently, for

The explicit expressions of the parabolic approximation for the functions

It is possible to improve the parabolic approximation in two ways:

(1) To express the numerical coefficients in formulas similar to Eq. (36) using analytic approximations for the roots of the equations

(2) To approximate the bumps of the functions

For an algebraic approximation of

proposed in [30].

The finite square well is a good starting point for similar quantum mechanical problems, i.e., the asymmetric well (when the walls of the well, see Figure 1, have different heights), the semi-infinite well (when one of the walls is infinite), or more realistic cases, when the walls are rounded (see [11, 12]). These potentials can model a semiconductor heterojunction (a thin semiconductor slice sandwiched between two different, larger semiconductors), a metallic film deposited on a semiconductor (in vacuum), and so on.

## 5. The differential form of transcendental equations

We shall indicate now an approach for solving the eigenvalue Eqs. (18)–(20) and (22)–(24) providing an exact solution, written as a series expansion. We shall first illustrate this method with the function

Taking the derivative with respect to

we get

Using Eq. (40) and taking into account that we are in the second quadrant

we obtain the differential form of the equation for

with the initial condition:

Putting

replacing

with the initial condition:

With Eq. (46), we can obtain the derivatives of any order of

The parameters

and so on. For the explicit expression of

Let us also remark that, in spite of the fact that the equivalence of Sturm-Liouville problems for electromagnetic fields and for wave functions was noticed many years ago, the results obtained for the finite rectangular well remain unused by the researchers studying wave propagation in wave guides or in other simple geometries. Reciprocally, the very detailed solutions of the equations for the normal modes of electromagnetic waves (see, for instance, the references [90, 92] in [31]) were apparently overlooked by researchers working in quantum mechanics.

## 6. Applications to the statistical physics of ultrathin metallic films

With few exceptions, the physics of ultrathin metallic films can be satisfactorily explained using different types of infinite well for the potential of electrons moving normally to the film plane. The model of the infinite well can be improved, for instance, by the phase accumulation theory [11, 12, 32], quite popular among the scientist working in surface physics. The theory satisfactorily explains the quantum scale effects (QSEs) appearing in such systems and predicted theoretically in the pioneering papers of Sandomirskii [4] and Schulte [5].

If, for thin films, such theoretical models can be successfully applied, for ultrathin films, with only few (typically, less than 5) monolayers, obtained experimentally in the last two decades, the approximation of the infinite well is inadequate. This is why in such cases we have to use the exact solutions for the bound-state energy of the finite well or, at least, their analytic approximations. In order to make clear the differences between the predictions of the two models—the first one is based on the infinite well, and the second one is based on the finite well—we shall evaluate some QSE for an ultrathin metallic film for three potentials: infinite, semi-infinite, and finite wells.

### 6.1. The infinite well model for the quantum well in an ultrathin metallic film

Let us consider a rectangular metallic films, with edges

and the electron energy is

where we put

The differences between the values taken by the integers

For ultrathin films, the discrete spectrum of

Let us consider a numeric example. For a metallic film with two atomic monolayers, the typical values are

We shall compute the number of occupied electronic states and the Fermi wave vector of the ultrathin film. The total number of subbands, which cut the Fermi sphere is

where

As the occupied states belonging to the subband of index

occupied states in the subband

where the factor of 2 is due to the electron spin. Putting

and introducing the number density of electrons

giving the dependence of the Fermi wave vector on the thickness

Introducing Eq. (61) in Eq. (57), we get

The last two equations define the QSEs on the Fermi wave vector; they can be considered as the starting point of all other similar QSEs of various physical quantities characterizing the ultrathin film.

Choosing

where

### 6.2. The semi-infinite well model for the quantum well in an ultrathin metallic film

As already mentioned (see Eq. (45) and the remark just below Eq. (27)), the relation between the solutions of the eigenvalue Eq. (46), namely, the functions

and the bound states of the semi-infinite well are described by the odd states of a finite well with the same length. In other words

It is convenient to define

So, the inverse strength of the quantum well, similar to Eq. (10), can be defined as

According to Eq. (68), the wave vector depends on both

with

An important difference which occurs at semi-infinite wells, compared to the infinite wells, is that it keeps a finite number of bound states. Consequently, the energy spectrum of the electron gas of the metallic film contains a finite number of subbands, in dependence of the value of

and exactly one state

This corresponds, usually, to a film with one or two monolayers. We have two states in the well,

This corresponds, usually, to a film with up to four monolayers, etc. These conditions are purely mathematical, i.e., consequences of the specific form of the eigenvalue equations.

Now, we shall impose physical conditions, due to the

Let us presume that the electron gas contains exactly

or

Therefore, instead of Eq. (62), we have

The term corresponding to the r.h.s. of the inequality (62) is missing in this case, as the number of roots (solutions) is completely determined by the condition imposed to

Replacing the electron number density with a typical value

This restriction on

#### 6.2.1. The finite well model for the quantum well in an ultrathin metallic film

The situation is quite similar to the previous one—the semi-infinite well. However, in this case, there is at least a solution for each value of

and so on. In Eqs. (76) and (77), the replacement

These solutions, or their analytic approximations (for instance, the cubic one), can be used directly in the models already proposed for the infinite well [33], in order to obtain the electron density, the surface free energy, the surface dipolar moment, or other similar quantities, in the more realistic case of a finite rectangular well.

## 7. Conclusions

This chapter illustrates how solutions of a simple quantum mechanical problem can be used for the description of certain interesting phenomena of nanophysics. Specifically, we referred to the exact solutions of the eigenvalue equations for the eigenenergy of the bound states of a particle in a rectangular well. If the physical problem is elementary, and the wave functions are simply written in terms of elementary functions, the equations for the eigenvalues of energy (or of the wave vector) are transcendental—and highly nontrivial. We obtain both exact solutions (series expansions) of these transcendental equations and approximate ones—with various degrees of complexity and accuracy. The value of the Fermi wave vector of the electrons in the metallic film, calculated for the finite well model, differs drastically from those calculated with the infinite well one.

Our results for the one-electron wave functions of the finite barrier model can be used as Kohn-Sham state in the self-consistent calculations of surface energy [34], for more accurate calculations of the stability of the films [1] and of other QSEs [33]. They can be also used as zero-order approximations for more realistic potentials, e.g., with rounded walls or undulate bottom—in a Rayleigh-Schroedinger or Dalgarno-Lewis perturbation theory [35].

Using the analogy between the movement of electrons in time-independent potentials and propagation of electromagnetic waves in dielectrics or metallic wave guides [18], mathematically, they are identical Sturm-Liouville problems; our results can be extended to several problems of electromagnetism and optics. This analogy can be easily developed for planar dielectric waveguides, namely, for “step-index” dielectrics, consisting of a slab of higher refractive index (core), sandwiched between two half spaces of lower refractive index (cladding). In such a situation, the quantum counterpart of the dielectric guide is a square well. This issue is discussed in detail by Casey and Panish in the context of heterostructure lasers [36]. It is easy to notice that the eigenvalue equations for transverse electric and magnetic modes, (2.4–45, 54, 60, 66) in [36], are essentially identical with our Eqs. (8) and (9).

## Acknowledgments

The author acknowledges the financial support of the IFIN-HH–ANCSI project PN 16 42 0101/2016 and of the IFIN-HH–JINR grant 04-4-1121-2015/2017.