^{1}

^{*}

^{1}

Theoretical solid-state physicists formulate their models usually in the form of a Hamiltonian. In quantum mechanics, the Hamilton operator (Hamiltonian) is of fundamental importance in most formulations of quantum theory. Mentioned operator corresponds to the total energy of the system and its spectrum determines the set of possible outcomes when one measures the total energy. Interpretation of results obtained by the applying of models based on the Hamiltonian indicates very specific mechanisms of some observed phenomena that are not fully consistent with the experience. Such approach may occasionally lead to surprises when obtained results are confronted with expectations. The aim of this work is to find Hamilton operator of acoustic phonons in inhomogeneous solids. The transport of energy in the vibrating crystal consisting of ions whose properties differ over long distances is described in the work. We modeled crystal lattice by 1D “inhomogeneous” ionic chain vibrating by acoustic frequencies and found the Hamiltonian of such system in the second quantization. The influence of long-distance inhomogeneities on the acoustic phonons quantum states can be discussed on basis of our results.

Many authors have been studied the influence of phonons on various processes in solids. Phonons affect certain aspects of the solid structure such as electronic density distribution [_{c} superconductors [

Procedure mentioned above is well applicable in fully homogeneous structures where phonon-mediated energy transport can be successfully introduced.

With the understanding of the importance of acoustic excitations in solids the procedure of canonical quantization attracts a significant interest. During past years some effort has been devoted to the understanding of the influence of quantization on the vibrational properties of heterostructures. Great interest has been oriented mainly on effects of acoustic-phonon quantization in restricted geometries [

In this paper we present a description of the long-distance structural inhomogenities influence on acoustic phonon quantum-mechanical states.

Main ideas of the work are outlined in the next sections which are organized as follows: We present motivation of the work in Section 2. Section 3 describes the dynamics of ions in the 1D model of “inhomogeneous” crystal lattice. Section 4 illustrates continuum limit for the model presented in Section 3. Section 5 examines canonical coordinates of the investigated system. Application of second quantization formalism to inhomogeneous vibrating solid is presented in Section 6. Section 7 is devoted to short discussion of our result. Finally we give a brief conclusion in Section 8.

If there are a weak interatomic forces in the ideal crystal, then the Hamiltonian of vibrating crystal can be written in the one-phonon approximation:

where is the phonon creation (annihilation) operator for the w-th phonon mode. In the harmonic approximation, a crystal can be described in terms of noninteracting phonons. The different phonon modes are independent and do not interact with each other. However, this simple picture is only an approximation, and in the case of large atom displacements, we should take into account the higher-order terms obtained by expanding of the potential energy. Keeping the cubic term, we can describe the interaction between phonon modes and thus we explicitly consider the phonon-phonon interaction between two modes [

The concept of phonon mentioned above remains valid when the anharmonic contribution in crystal vibration is small compared with the harmonic. In this case, the quasi-particle approximation can be made as well. Then some of the phenomena can be explained on the basis of the idea of interaction of phonons with other elementary excitation in solid-state matter (excitons, magnons, plasmons, etc.).

The question is, what is the effect of long-distance spatial inhomogeneities of structure in solids on acoustic phonon states and how it is necessary to modify results of “phonon based” quantum mechanical models of some phenomena in the case of inhomogeneous continuum. In the modern literature, numerous investigations were, devoted to the chracter of phonon mediated interactions in solids. However, only a little work had been done about the quantum-mechanical states of phonons in inhomogeneous solids. Thus, it is very difficult to evaluate theoretical and practical importance of consideration of structure inhomogeneities in the frame of second quantization method. Next we propose the basis for the solution of this problem.

Next, the 1D ion chain model is studied (

Consider N ions of unequal masses and charges with nonlinear nearest-neighbor interactions, described by a potential energy U. The position of the n-th ion will be denoted by x_{n}, where the ions are numbered in order of increasing value of their axial positions, so that implies that (see

Ions are bound relatively strongly in their equilibrium positions in the crystal lattice. The following analysis is based on the classical description of the dynamics of ions in the 1D crystal lattice in harmonic approximation.

We consider one-dimensional crystal represented by a linear chain of N ions (

If displacements of ions from the equilibrium are small, the potential energy U can be expanded to the Taylor series. Next harmonic approximation can be applied, i.e. only the second degree polynomial function can be considered. Moreover is it possible to take account the fact that potential energy U has a minimum for each equilibrium positions. Therefore, the potential energy U in the harmonic approximation is given by:

where:

and. (3.3)

Then the force acting on the n-th ion in chain can be written as:

In case if displacements of all ions from equilibrium are the same, the force (3.4) acting on any ion in chain is zero because it corresponds to the displacement of whole crystal:

This means that the coefficients are not linearly independent. Next condition results from the relation (3.5):

Assuming inhomogeneous case of 1D chain there are different masses of individual ions (see

and applying the tight binding approximation the equation of motion for the n-th ion in the chain is given by:

In addition, as can be seen from the Equation (3.6):

the Equation (3.8) can be written in the next form:

System of Equations (3.10) for allows to describe the dynamics of ions in the frame of the proposed 1D chain model. Solution of this system is the set of functions, which describes time dependence of the displacements of ions from equilibrium in the chain.

We consider the continuum limit for model of solid presented above. The term “continuum limit” usually relates to discrete models where it means “becoming less discrete”. To make this let us consider a transitions:

where d is the distance of neighboring ions in the chain acting oscillations.

Linear density and elastic coefficient change in the inhomogeneous case along the chain. They depend on coordinate x in the 1D model (see

and next we suggest to replace discrete variables in the Equation (4.2) by continuous functions on the basis of (4.1), i.e.:

After substituting of relations (4.3) into Equation (4.2) we get:

We believe that the Equation (4.4) describes the transport of acoustic waves in inhomogeneous continuum in 1D case.

On the basis of the result (4.4) it can be concluded that transport of acoustic waves in 3D case can be described by next equation:

Lagrangian formulation of the oscillating system mechanics is based on the formal introduction of canonical variables to Equation (4.5). Firstly, it is necessary to find Lagrangian L of the oscillating system considered above. Suppose that after the substituting of the Lagrangian L into the Euler-Lagrange equation:

this equation changes into the form (4.5). Consider the Lagrangian in an analogous form as in the homogeneous case:

where we integrate over the entire sample volume. One can easily check that by the substituting of the relation (5.2) into Equation (5.1) this equation takes the form (4.5).

Displacement from equilibrium state is a one of complex conjugated variables-canonical coordinate. Complex conjugated momentum can be found using the next relationship:

By substituting (5.2) into (5.3) the conjugated momentum can be written in the form:

Then the Hamilton function of oscillating inhomogeneous continuum is determined by formula:

Canonical variables of the inhomogeneous oscillating continuum can be written as:

We propose the canonical coordinates of the investigated system in an analogous form as in the homogeneous case. The question is, what is the impact of spatial dependence of mass density and elastic coefficient on the energy states of the system.

For the quantization of Equation (4.5) it is necessary to define operators corresponding to physical quantities. In the considered case it is necessary to replace the canonical variables by their operators:

Mentioned operators satisfy well known commutation relations:

Subsequently, the coefficients in the relations (5.6) and (5.7) become operators defined by the next way:

Formulas for the annihilation and creation operators and can be found by inversion of relations (6.5) and (6.6). The quantization of the elastic vibrations is made through the commutation relations:

We assume that operators and obey mentioned commutation relations even in the inhomoneneous case. Hamiltonian of the vibrating system can be obtained easily by the replacing of the canonical variables in formula (5.5) by the operators (6.5) and (6.6), i.e.:

It is necessary to substitute relations (6.5) and (6.6) into the Equation (6.8) and determine the Hamiltonian by means of operators and. If we raise to the second power the formula (6.6) and next we substitute it to the first integral in (6.8), we get:

i.e.:

Because it holds:

;

where a = 1, 2; b = 1, 2, the Equation (6.9) goes into the next simple form:

We take into account the fact that:

It is necessary to find the gradient of the operator (6.5) for the calculation of the second integral in the expression (6.8). We must take into account the spatial dependence of in this calculation. The gradient of operator (6.5) may be determined as follows:

In order to find the operator we must raise to the second power the Equation (6.12) and next substitute the obtained result into the second integral in Equation (6.8).

After the substituting of (6.12) to the second integral in the (6.8) it takes the form:

where:

and.

Operator in (6.13) is determined as follows:

where.

Functions (for n = 1, 2, 3, 4) are related to spatial inhomogeneities of investigated continuum:

Mentioned functions can be expanded into the Fourier series:

where:

and then coefficients determined by formula (6.14) take the form:

Next the coefficients in (6.18) can be determined by means of integrals (6.17). It is necessary to substitute functions in the form (6.15) to mentioned integrals. For example:

By the introduction of the:

the coefficients take the form:

Then we consider:

;

where a = 1, 2; b = 1, 2 and the relation (6.13) can be simplified:

Finally after the substituting of both results (6.10) and (6.27) into the Equation (6.8) the Hamilton operator of the inhomogeneous continuum vibrating by the acoustic frequencies can be written as:

where:

Operator must be hermitian (self adjoint) and therefore can be modified to the following form:

where coefficients are:

We believe that the second term in the Hamiltonian (6.28) determined by Equation (6.30) can make the basis for quantum-mechanical analysis of acoustic phonons scattering process in solid structure caused by long-distance spatial inhomogeneities.

In this paper we presented a simple approach to the quantization of energy of acoustic waves in inhomogeneous continuum. A fully canonical quantization of mechanical energy of vibrating media was introduced in the presence of inhomogeneities of structure represented by the spatial dependencies of mass density and elastic coefficient. Subsequently, the Hamilton operator of the vibrating inhomogeneous continuum was determined by standard technique of second quantization. The contribution of the work lies in the fact that our result makes it possible to estimate the impact of long-distance inhomogeneities on the energy spectrum of acoustic phonons.

Hamiltonian of the inhomogeneous continuum vibrating by acoustic frequencies can be written in the form (6.28) where is well-known Hamiltonian of homogeneous vibrating system identical to (2.1). Operator is generated by inhomogeneities. In that case the coefficients can be interpreted as phonon scattering form factors that reflect the structure of an extended energy distribution. Mentioned coefficients determined by integrals (6.31) represent non-locality of phonons scattering centre and in the case of homogenous continuum disappear.

Summing up the above-said, one can conclude that nonlocal changes of mechanical properties affect the mechanical energy of acoustic vibration of the continuum. This may affect the formation of acoustic excitation and cause anomalies in related phenomena. Possible influence of long-distance structural inhomogenities on acoustic phonon spectrum is expected. For example a heat capacity anomalies in some materials observed at low temperatures by many authors [15,16] could be partly caused by their structural transition, which is associated with origin of structural inhomogeneities. In our opinion the influence of these structural inhomogeneities on the low-T behavior of heat capacity of such system could be discussed on the basis of our results.

Finally, we would like to point out that the presented contribution is only intended to preliminarily obtain some aspects of second quantization in inhomogeneous cases which could be useful within other theoretical models.