The proposed nonstandard diffraction theory is constructed directly from the Maxwell equations for the crystalline medium in the X-ray wavelength range. Analysis of Maxwell?s equations for dynamic diffraction is possible using the method of multiple scales which is modified to the vector character of the problem. In this case, the small parameter of the expansion is the Fourier component of the polarizability of the crystal. The second-order wave equation is analyzed without any assumptions about the possibility of the interaction between the refracted and scattered waves which automatically leads to the dynamic character of the scattering. The unified consideration of different geometrical schemes of diffraction including grazing geometry is possible. This is due to the construction of a unified wave field in the crystal and obtaining the field amplitudes according to the boundary conditions. The proposed theory allows generalization to the case of an imperfect crystal. Thus, a unified approach to account for deformations and other crystal structure disturbances in all diffraction schemes is implemented. The determination of a unified wave field without separation of the refracted and scattered waves is of the greatest importance in the analysis of secondary processes.
- X-ray diffraction
- dynamic theory
- imperfect crystal
- perturbation theory
- the method of multiple scales
- extinction length
- boundary conditions
- reflection coefficient
The Darwin theory  is based on the Bragg model for the crystal considered as a family of parallel crystal planes. The X-ray wave reflection is considered as a result of successive transmission and multiple reflections from planes. In this case determination of the diffracted wave amplitude is reduced to the solution of recurrent relations between the amplitudes of transmitted and scattered waves in passing through the specified atomic plane. In essence, the Darwin theory represents a direct extrapolation of the optical task of propagation of light in a layer continuum to the case of wavelengths of the X-ray range.
The Evald-Laue theory  was the next stage in the development of theoretical approximations about the character of X-ray wave propagation in the crystal under dynamic scattering. The model approximations of the X-ray crystal interaction were formulated in the framework of this theory, meaning that the X-ray wavelength is comparable to the interatomic distances.
Therefore, the standard continual approximation for electrodynamics of continua proves to be unacceptable, and the scattering from individual charges should be taken into consideration. As is known, taking this into account results in the formalism of 3D periodic dielectric permeability ε(
In spite of a series of unconditional advantages in the interpretation and theoretical prediction of experimental results on dynamic X-ray scattering in crystals, both the Evald-Laue theory and, substantially, the Darwin theory have the principal limitation that they describe the dynamical diffraction in perfect crystals only.
The necessity of taking into accounts the different deviations from ideal periodicity in the crystal and, first of all, deformations resulted in the creation of the generalized theory developed by Takagi and Taupin . This is based on the approximation of a wave field in the form of superposition of the transmitted and diffracted waves with slowly varied amplitudes depending on the coordinates that leads to the Takagi-Taupin equations to the system of differential equations relative to the field amplitudes. This formalism gives the possibility to describe the dynamical diffraction in a distorted crystal since the distortions of ideal periodicity can be taken into account in the explicit form in the wave field approximation. Correspondingly, the Takagi-Taupin equations become the system of differential equations with variable coefficients.
It is important that the Takagi-Taupin system is shortened and the coordinate second derivatives of field amplitudes are neglected in it. On the one hand, this significantly facilitates the theoretical consideration and makes observable the solution of a series of diffraction problems in standard diffraction geometries when this simplification proves to be justified. On the other hand, the Takagi-Taupin equations become inapplicable under the conditions, for example, of grazing diffraction geometry; then, it is necessary to solve the third- or even fourth-order differential equations .
The principal disadvantage of the procedure of equation shortening is related to the impossibility to state correctly the boundary conditions at the crystal-vacuum interface for field amplitudes. Instead of the known classical continuity conditions for tangential components of electric and magnetic fields, the boundary conditions of the type of setting of the normal components of field amplitudes on the crystal surface become vaguely clear but not in line with the Maxwell equations. Of cause the solutions of these boundary problems prove to be applicable only for rather large (more correctly significantly exceeding the angle of full external reflection) angles of radiation incidence and yield.
As a result the theory faces the difficulties related to the necessity to solve the third- or fourth-order equations which become virtually overwhelming for the case of the crystal with lattice deformation when the diffraction schemes of the type of sliding diffraction are considered.
At the same time, the Maxwell equations are the first-order equations or the second-order ones in the case of a single wave equation when passing, for example, to an electric field. Namely, this structure of equations agrees with the mentioned classical boundary conditions. This means that the requirement of taking correctly into account the boundary conditions in any theoretical diffraction scheme leads virtually unambiguously to the known structure of wave equation following from the Maxwell equation.
Thus, to overcome the aforementioned difficulties, it is necessary to create the theory directly based on the Maxwell equations using model approximations of crystal polarizability in the X-ray wavelength range.
2. A covariant dynamical theory of X-ray scattering in perfect crystals
2.1. Physical diffraction model and basic equation
The model approximations for crystal and X-rays propagating in it are reduced to the following.
The plane monochromatic wave falls from vacuum onto the crystal. The crystal-vacuum interface is considered as a geometrical one so that the classical boundary conditions for optics prove to be applicable. The uniform wave field in the crystal is described by the fundamental Maxwell equations supplemented by the constitutive equation
The Maxwell equations, as the equations of electromagnetic waves in dielectric, when dispersion is absent, can be written in the form (designations here and below are standard):
These equations are supplemented by the constitutive relations:
As usual, we will consider that the time dependence of E and H is harmonic:
We have the system
For economy we use here and below factor 2π into k, 2πk → k. We come from system (6) to the basic equation in the standard way:
In the present section, we will consider the case of a perfect crystal. This will allow us to compare the conclusions of the proposed formalism to the known theoretical results.
We choose the crystal model χ(
Below, we will bring 2π into
To take correctly into account the contribution of different terms into Eq. (9) when using the methods of perturbation theory, it is necessary to bring it to the dimensionless form. This procedure assumes the choice of a certain characteristic spatial scale of the problem. Apparently, this is the parameter determining the reciprocal lattice in our problem, namely, the modulus of the reciprocal lattice vector:
Eq. (10) cannot be exactly solved. Correspondingly, it is necessary to use some method of approximate solution. There are two main ways to analyze Eq. (10). In the first case, taking into account that χ(
The approximate solution is searched in the space beyond the scattering crystal at large distances from it in the form of the first term of a series of the Born expansion. The applicability of this approach is limited by the smallness of the scattering cross section as compared to the geometrical area of the crystal section. As is known this approach results in the kinematical theory .
In the second case, two variants of the dynamical theory are developed from Eq. (10). For the first variant, the solution of Eq. (10) is sought in the form of Bloch wave represented in the form of an infinite series of plane waves with wave vectors corresponding to the refracted wave and diffracted waves in the crystal. This Bloch wave is interpreted as a multiwave solution of the dynamical theory. As a result, the use of expansion of χ(
Since it is not possible to solve an infinite system, one has to be limited as a rule by two equations, that is, by the two-wave approximation.
For the second variant of the theory, the solution of Eq. (10) is presented in the form of plane waves with the slowly varying amplitude. As a result of the two-wave approximation, we obtain the Takagi-Taupin equations which can be interpreted as the recurrent Darwin relations written in the differential form .
We propose here to use a new approach to the analysis of Eq. (10). The physical justification of the proposed method indicates the choice of the perturbation parameter and is as follows. The propagation of an X-ray wave in the crystal unaccompanied by the appearance of diffracted beams is adequately described by the uniform wave equation with ε = 1 + χ0. This corresponds to the propagation of an X-ray wave in the crystal as a continuum with the refraction factor with respect to the usual laws of optics. This situation can be considered typical. On the contrary the appearance of diffracted beams requires that the definite geometrical conditions be fulfilled for the wave vectors and reciprocal lattice vector. Apparently, χ
Thus, in spite of the fact that all quantities χ0, χ
2.2. Direct expansion and geometric diffraction conditions
We use the simplest perturbation method, direct expansion, in χ
We will restrict ourselves to the first-order expansion. We have
The zeroth approximation corresponds to zero power of the perturbation parameter χ
We select the solution of this vector wave equation in the form of superposition of two plane waves:
The reasons for this selection are the following. The zeroth approximation corresponds to the propagation of two plane transverse waves in opposite directions in the continuum with χ = const. This is a singular analog of the total field in the crystal for the case of an empty lattice. The propagation directions and the wave amplitudes remain indeterminate and are specified below by the boundary conditions at the vacuum-crystal interface.
The first approximation is obtained when all terms in Eq. (11) proportional to the first power of χ
We obtained the inhomogeneous wave equation. According to the perturbation theory, it is necessary to find its particular solution. Since the continuum is uniform with an accuracy of the reciprocal lattice vector
takes the form
Eq. (13) is obtained taking into account the condition
from which follows
Now, the solution in the first order of the perturbation theory can be written:
Finally, the direct expansion with an accuracy of takes the form
It follows from Eq. (15) that in addition to the direct (
However, this position radically changes when any denominator in Eq. (15) approaches zero. In this case
, and we cannot consider a small correction to
This condition is well known: it is the Laue condition for X-rays, and therefore there is no need to detail its physical sense. Note only that all geometric constructions following from the Laue condition appear in this case as a natural consequence of validity violation of the direct field expansion in the parameter χ
Thus, the wave field structure principally changes for certain
Thus, it is necessary to modify the direct expansion near the
There are different methods to modify the direct expansion. All of them are directed to solve one problem: to obtain a so-called uniformly acceptable expansion near the values of parameters interesting for us. The method of multiple scales is most favorable for our investigation . However, method modification is necessary having in mind the vector character of the problem.
2.3. The method of multiple scales
The main idea of the method for the considered problem is the following. The wave field singularities appear for different spatial scales determined by a small parameter of expansion of χ
Thus, we will search for the approximate solution of Eq. (10) for the most interesting case in the area of the Bragg maximum when the Laue condition
Hereinafter, we will restrict ourselves to the first order of expansion; correspondingly, we will consider two spatial scales
Then, for rot
that is, operator rot is linear relative to the carried-out substitution. The index of the operator signifies the space in which it operates. Using this property we obtain
As was indicated above, the interaction of the field with continuum has a parametric character. This means that along with the field expansion it is also necessary to expand the wave vector
We substitute all expansions into Eq. (10):
In this case, χ(
Namely, the initial approximation (unperturbed state) and subsequent ones are obtained when the constants of powers of the perturbation parameter χ
Let us demonstrate this procedure. As before the zeroth approximation, (the unperturbed equation) has the form of the standard vector wave equation for transverse waves propagating in continuum:
However, in contrast to the direct expansion, the operator rot acts here only on the single spatial scale
This structure of the wave field supposes strict fulfillment of the diffraction Laue condition
The condition of wave transversity (
The following first-order approximation to leads to the inhomogeneous equation:
Operator (gradient) acts here on r1. Eq. (21) is obtained taking into account and the additional condition . This means the following. The quantities are considered as the perturbed amplitudes of the corresponding plane waves. The two first terms in the right side of Eq. (21) generate the secular components in expansion which is seen from the abovementioned particular solution. In other words they provide parametric resonance in the system. Consequently, it is necessary to eliminate these terms in order to obtain the uniform approximation near the Laue condition. Then, we obtain the following system of vector equations:
In contrast to the usual scalar system, additional limitations are necessary to solve Eq. (22). The issue is that
The particular solution of the inhomogeneous equation
in the resonance case takes the form
It is seen from here that the infinite increase in the wave amplitude is associated with the wave vector
in the direction coinciding with
Thus, multiplying scalarly the first equation of the system by
The obtained system (23) is virtually the dispersion relation written in the differential form for the transmitted and diffracted waves in the two-wave approximation. As follows from the presented conclusion, the possibility to obtain this system is dictated by the choice of the zeroth approximation. Namely, the wave vectors of transmitted and diffracted waves must have the identical component that leads to the scattering interpretation as a result of reflection from the atomic plane.
We pass from the differential form of system (23) to the algebraic one; for this purpose we make the following substitution
This system should be considered as the condition of field limitation in the specified direction of field propagation that is related to the structure parameters of the crystal and to the diffraction geometry. As for the physical meaning of the considered problem, this limitation should be in the direction of the normal inside crystal.
Here, are the direction cosines of the corresponding wave vectors, = .
The nontrivial solution of Eq. (25) requires the corresponding determinant to vanish:
The solution of quadratic equation relative to
Then, solving, for example, the first equation of system (25) relative to
Finally, going back to the initial variables, we represent the wave field in the crystal in the following form:
Here, the constants
The equation for the wave field is simplified for the case of semi-infinite crystal when the wave reflection from a lower crystal face is eliminated. In this case the choice of sign in
Finally, the wave field in the crystal takes the form
where the constant
Eq. (32) describes the uniform wave field in a perfect crystal near the Bragg maximum.
To compare the obtained expression to the known results and to the experiment, it is necessary to go to the angular variable which connects with the deviation from the exact Bragg angle Δθ. In this case, the parameter
should be expressed by Δθ. As was mentioned the parametric resonance (diffraction Laue condition
) is determined only by the component
Here, we introduce the standard (with an accuracy of refraction) angular variable:
As is known, two principal schemes are considered in the diffraction theory, by Laue ( ) and by Bragg ( ). In the case of Bragg diffraction, the wave field structure will qualitatively differ depending on the considered angular range.
In particular in the range
the waves will exponentially attenuate along the normal surface into the crystal. In terms of the qualitative theory of differential equations, the solution will be unstable. The known interpretation [10, 11] leads to the conclusion of the expulsion of the wave from the crystal and the formation of a diffraction maximum. Thus, the indicated condition separates the stable solutions (oscillating type) from unstable ones (exponential type), that is, provides the equations of transition curves of the parametric plane ( ) [10, 11].
The width of the unstable region, the region of exponential wave attenuation, is determined by the expression
or proceeding to the angular variable
This is the known expression (with an accuracy of refraction) for the angular width of the Bragg table for the case of the semi-infinite perfect non-absorbing crystal. The extinction length is determined as a decrement of wave attenuation in the point Bragg position:
We now summarize the intermediate stage. The use of the generalized method of many scales allowed us to obtain the system of basic equations describing the behavior of wave field near the Bragg maximum in the two-wave approximation. This system is a direct analog of the dispersion relations of the Ewald-Laue theory and the Takagi-Taupin system of the generalized dynamic theory. The substantial difference of the developed variant of the theory consists in the cancelation of the shortening procedure of equations by neglecting the second derivatives. The principal moment here is the expansion in χ
Comparison of the obtained results to that known from the Takagi-Taupin theory shows the complete correspondence both in the qualitative interpretation of types of the solutions obtained in the different angular ranges in the case of Bragg diffraction and in the analytical expressions for the width of the Bragg maximum and the extinction length. This correspondence indicates that in spite of the formal representation of χ(
However, the value of the theory developed here shows itself to a large degree when the boundary conditions are taken into account that takes the explicit expressions for the reflection coefficient. Therefore, we consider now the boundary conditions and determine the principal difference between our approach and known variants of the dynamical diffraction theory.
2.4. Boundary conditions and the amplitude reflection coefficient
According to Eq. (32), the expression above obtained for the wave field in the crystal depends on the constant
We find out now the differences appearing when the boundary conditions are strictly taken into account in our theory.
The reflection coefficient is determined as a ratio of the averaged values of normal components of the pointing vector for the diffracted and incident waves:
As is known, the boundary conditions require continuity of the tangential components of electric and magnetic fields which is a consequence of the uniformity of the problem along the surface. The boundary problem disintegrates into successive steps related to the determination of . In this case the elementary problem is considered at each stage, namely, the determination of the relation between the amplitudes of incident, transmitted, and secularly reflected waves. The solution of this problem leads to the known Fresnel formulas:
The obtained relations allow us to determine not only the diffracted wave but also the specularly reflected wave, which principally distinguishes our approach from the formalism of the Takagi-Taupin equations.
Eqs. (38)–(40) solve the problem of the determination of field amplitudes under the conditions of sliding noncoplanar diffraction when the incident and diffracted waves are near the critical angles of total external reflection (TER). They are analogous to the relations obtained in  where this problem was solved by the fourth-order dispersion equation.
The correspondence with the Takagi-Taupin theory must be undoubtedly fulfilled for the case of large angles of incidence and yield of the diffraction wave. Really in this case, the amplitude of specular wave tends to be zero, and the reflection coefficient takes the form
This is the known expression for the coefficient of reflection from a perfect half-infinite crystal, i.e. the Bragg table. In addition to this in the case of extremely asymmetric diffraction when the diffraction wave leaving the crystal is almost parallel to the surface, the amplitude is modulated by the factors taking into account the refraction of transmitted and diffracted waves at the crystal-vacuum interface and the diffraction wave interaction related to the vector
The proposed covariant (it may be named nonstandard) theory allows the generalization for the case of the crystal with lattice deformations. Therefore, a uniform approach to the account of deformations and other distortions in all diffraction schemes is realized.
3. Generalization of a covariant dynamic diffraction theory to the case of deformed crystal
It is known that Takagi-Taupin equations were obtained using the model concept of the character of lattice distortions which makes it possible to directly take into account the displacement of atomic planes in the three-dimensional periodic function of crystal polarizability. This concept allows for describing lattice displacements and strain using the methods of classical theory of elasticity formulated within the continuum approximation. To satisfy the condition of the dynamic character of scattering in the Takagi-Taupin theory, the lattice distortion is assumed to be rather weak; correspondingly, the strain is small. The character of variation in strain is implicitly taken into account only when the wave field is chosen in the form of a Bloch function with slowly varying amplitudes; thus, the question of applicability of this concept remains open.
There is another limitation of the Takagi-Taupin theory which is related to the correct statement of boundary conditions. Mathematically, the Takagi-Taupin equations form a first-order differential system with respect to the scalar amplitudes of transmitted and diffracted waves. The procedure of determining these amplitudes at the crystal-vacuum interface does not correspond to the classical boundary conditions.
This discrepancy is due to the fact that the Takagi-Taupin equations are obtained disregarding the second derivatives of the field amplitudes with respect to coordinates. Thus, the boundary conditions impose fundamental limitations on the applicability of the Takagi-Taupin equations (e.g., when analyzing extremely asymmetric diffraction schemes).
In this section we generalized the covariant theory of dynamic diffraction which was presented in the previous section, to a crystal with a distorted lattice. This approach makes it possible to formulate the limitation on the character of variation in strain for the applicability of the Takagi-Taupin equations. In addition, the equations obtained can be applied (as in the case of an ideal crystal) for arbitrary diffraction schemes.
3.1. Wave field in the absence of diffraction
The polarizability χ(
It can be seen that the general structure of χ(
As in the case of an ideal crystal, Eq. (43) should be reduced to a dimensionless form. To this end we will use the length
The approximate solution to Eq. (44) in the Takagi-Taupin theory is known to be sought after in the form of plane waves with slowly varying amplitudes. Finally, the Takagi-Taupin equations can be derived from Eq. (44) with allowance for the two-wave approximation .
The covariant dynamic theory in the case of an ideal crystal is due to the fact that the Fourier component of polarizability χ
Recall that in the case of an ideal crystal the direct expansion of the solution to Eq. (44) in the parameter χ
In this case, the scattered wave amplitude increases unlimitedly. This is due to the fact that the diffraction wave can be found as a particular solution to the inhomogeneous wave equation:
which according to Eq. (13) has the form
For a deformed crystal, the plane wave is replaced with
At the same time from the physical point of view, the displacement field changes significantly at distances much larger than the lattice parameter. This limitation is substantiated in particular by the fact that we describe lattice distortions within the continuum approximation. In this case the approximate solution to Eq. (46) can be obtained similarly to Eq. (47):
Correspondingly, the direct expansion of the solution to Eq. (44) for a deformed crystal is similar to that for an ideal crystal and is determined by the wave superposition in the form
The most general limitation on the strain in the crystal is imposed by the requirement for the smallness of the strain tensor
where ε0 is the strain amplitude in the structure. For sufficiently regular displacement fields, this requirement is reduced to the condition ε0 ≪ 1. It follows from Eq. (50) that in the case of a deformed crystal when the aforementioned conditions are satisfied the applicability of direct expansion remains limited because of the Laue resonant condition (45) where the amplitudes of waves excited in the crystal increase unlimitedly. Thus, the method for obtaining an approximate solution must also be modified in the case of a crystal with a distorted lattice. The choice of the modification technique depends on the character of the displacement field in the crystal and obviously cannot provide a universal solution for all physically possible cases. We will consider the most widespread situation where the displacement field changes at distances comparable with the extinction length. In this case the multiscale method which is the basis of the covariant theory of diffraction in an ideal crystal can directly be extended to a deformed structure.
3.2. Derivation of the main equations for strained crystal
The strain field in a crystal may have various forms depending on the nature of lattice distortions. These forms can mathematically be represented by setting different structural parameters (e.g., the thicknesses of epitaxial layers and transition regions between layers, sizes of lattice-strain regions caused by various defects, etc.). The values of these parameters are determined by not only the strain amplitudes ε0
Thus, the parameter χ
Thus, a consideration of different types of lattice distortions generally calls for taking into account different characteristic spatial regions; this approach completely corresponds to the main concept of the multiscale method . Obviously, different modifications of the method are required depending on the specific structure of Eq. (52). We will consider the simplest case of one scale which leads in a particular case to Takagi-Taupin equations.
Let us consider the atomic plane displacement occurring at some effective layer thickness
where is the displacement model and ε0 is the strain amplitude.
According to the multiscale method , the main equation (Eq. (44)) is analyzed near the Bragg maximum on different spatial scales (determined by χ
Then, according to Section 2, the field, the operator rot, and the wave vector
In contrast to the case of an ideal crystal, χ(
The parameter 1/
We will follow the scheme for solving Eq. (55) that was reported in Section 1. As for an ideal crystal, the initial approximation can be written in the form as in Eq. (18). Accordingly, the operator rot acts on only one spatial scale
The first-order approximation with respect to χ
Here, the gradient
acts in the space
As can easily be found from Eq. (47), two first terms in the right-hand part of Eq. (60) generate the divergence of a particular solution. The following vector system can be obtained by excluding these terms from Eq. (60):
According to Section 2, to satisfy the field boundedness condition, the projections of Eq. (61) on the corresponding unit vectors
According to Section 2, the parameter can be expressed in terms of the deviation from the exact Bragg angle Δθ as follows:
Eq. (62) describes the changes in the wave amplitudes on the scale
Accordingly, Eq. (62) is reduced to the form
In a particular case of coplanar diffraction in the
Under the assumption that the displacement field in Eq. (64) depends on only
3.3. Reflection coefficient
System (64) has a more general character because it was obtained without additional limitations on the Takagi-Taupin equations which is related to the rejection of the second derivatives of the field amplitudes with respect to coordinates. It is known that the boundary conditions at the crystal-vacuum interface cannot be correctly taken into account due to these limitations. Finally, the Takagi-Taupin equations cannot be applied to extremely asymmetric diffraction schemes.
Let us show how the consideration of the boundary conditions yields explicit expressions for the amplitudes of diffracted and specularly reflected waves for arbitrary angles of incidence. For simplicity we will consider diffraction in a semi-infinite crystal in the case of σ polarization where the vectors
where the constant
Finally, the wave field in the semi-infinite crystal near the Bragg angle of incidence can be described by an expression that formally corresponds to an ideal crystal:
The procedure of solving the boundary problem which can be reduced to a successive establishment of relations between wave amplitudes having common tangential components of the electric and magnetic fields at the crystal-vacuum interface remains the same. Therefore, one can use the expressions for amplitudes obtained in Section 2. As a result we arrive at formulas that are similar to the Fresnel formulas:
The following designations are introduced here:
Eqs. (68)–(70) can be used to find the field amplitudes for arbitrary angles of incidence including those in the vicinity of the critical total reflection angles. In particular the formula for reflectance can be found from Eq. (70) as follows:
In a particular case of large angles of incidence, Eq. (71) can be simplified, and Eq. (70) yields the following expression for the reflectance (which can be derived from the Takagi-Taupin equations):
The variant of the dynamic X-ray diffraction presented in the present work is based on direct analysis of Maxwell equations for the definite model representations of the field-medium interaction taking into account the lattice presence which agree as a whole with the Ewald-Laue theory. This analysis proves to be available when the method of many scales adapted to the vector character of the problem is used. In this case the magnitude is the parameter of expansion that corresponds in full to the physical character of the problem. This correspondence is reflected in the mathematical structure of the analyzed field equation in the crystal under conditions of dynamic scattering.
The expressions obtained for the main field characteristics in the Bragg maximum region following from the qualitative singularities of the field propagation correspond to the known results of the dynamical theory. However, the correct use of the boundary conditions leads to an expression for the reflection coefficient that substantially differs from the classical one for the case of extremely asymmetric diffraction schemes. In addition the presented approach provides the amplitude of specularly reflected wave under conditions of dynamic diffraction, which cannot be apparently obtained in the framework of traditional approaches.
In the present work, we do not state the problem to analyze the features of dynamic scattering in the sliding diffraction geometry.
In conclusion, we note the most important in our opinion differences and advantages of the approach developed in the present work.
The second-order wave equation analyzed without any additional assumptions of the possibility of the interaction of refracted and scattered waves automatically results in dynamical scattering character; in this case the kinematical scattering can be considered to a certain extent as an artificial process having limited applicability. The diffraction Laue conditions appear as a result of natural limitations of the direct expansion of the solution in the resonance case.
In the framework of the developed theory, the total consideration of different geometrical diffraction schemes including sliding geometry and other surface variants proves to be possible. In this case the order of the dispersion equations does not change. This situation is related to the effective decomposition of the problem into the construction of the uniform wave field in the crystal and the determination of field amplitudes according to the boundary conditions.
Determination of the wave field as a whole without decomposition into refracted and scattered waves is the advantage of the theory. It is clear that this feature of the theory is most important for analysis of secondary diffraction processes.
We have generalized the covariant theory of dynamic X-ray diffraction to the case of a crystal with lattice deformation. In this case the displacement field is specified, starting from model representations used in Takagi-Taupin dynamic theory. In our case (in contrast to the formalism of the Takagi-Taupin equations), lattice distortions have been taken into account on various spatial scales that were different from the scale of the lattice period.
The displacement field was also a slowly varying function of coordinates. If the displacement field is considered on one spatial scale on the order of the extinction length, then the particular case of fundamental equations for the field amplitudes is obtained as a result.
In precisely this case, we have the same result as that of the Takagi-Taupin equations. By doing so we have shown possible restrictions on the applicability of the Takagi-Taupin equations to describing dynamic diffraction in crystals using various deformation models.
At the same time, the presented theory offers an opportunity for successively taking into account displacement fields of various types implemented on different spatial scales (that are larger or significantly smaller than the extinction length).
The possibility of the correct application of boundary conditions including cases of extremely asymmetric diffraction schemes in covariant theory for ideal crystals is also wholly retained for crystals with lattice distortions. Such a situation is due to the fact that the solution of the diffraction problem proper is not related to the boundary conditions; in particular the order of fundamental equations of the theory remains the same for arbitrary diffraction geometry.
James RW. The Optical Principles of the Diffraction of X Rays. Vol. II. London: G. Bell and Sons Ltd; 1962. p. 664
Pinsker ZG. X-Ray Crystal Optics. Moscow: Nauka; 1982. p. 392. [in Russian]
Andreeva MA, Kuzmin RN. Mössbauer and X-ray surface optics. Moscow: Obthcenatc. Akad. Znaniii; 1996. p. 128. [in Russian]
Dyshekov AA. Non-Standard Theory of X-Ray Scattering. Aktualnye Voprosy Sovremennogo Estestvoznaniya. 2009. No. 7. pp. 3-20. [in Russian]
Dyshekov AA. Spatial Scaling in the Theory of X-Ray Scattering. Non-Standard Dynamic Theory. Metallofiz. Noveishie Tekhnol. 2010. Vol. 32. No 1. pp. 13-21. [in Russian]
Dyshekov AA. A nonstandard dynamical theory of x-ray scattering in crystals. Journal of Surface Investigation: X-ray, Synchrotron and Neutron Techniques. 2010; 4(6):956-964
Dyshekov AA, Khapachev YuP. Novye Analiticheskie Podkhody k Zadacham Rentgenodifraktsionnoi Kristallooptiki. Nalchik. KBSU. 2010. p. 75. [in Russian]
Landau LD, Pitaevskii LP, Lifshitz EM. Electrodynamics of Continuous Media. 2nd ed. Oxford: Butterworth-Heinemann; 1984. p. 460
Nayfeh Ali H. Perturbation Methods. New York: The Wiley Classics Library; 2000. p. 437
Khapachev YP. The theory of dynamical X-ray diffraction on a superlattice. Physica Status Solidi (b). 1983; 120(1):155-163
Molodkin VB, Shpak AP, Dyshekov AA, Khapachev YuP. Dynamic scattering of X-ray and synchrotron radiation in superlattices. X-ray diffraction crystal optics of superlattices. Kiev. Akademperiodika. 2004. p. 120. [in Russian]
Bushuev VA, Oreshko AP. Specular X-ray reflection from a crystal coated with an amorphous film under the conditions for strongly asymmetric noncoplanar diffraction. Physics of the Solid State. 2001; 43(5):941-948
Dyshekov AA. Generalization of a nonstandard approach in dynamic diffraction theory to the case of deformed crystal. Journal of Surface Investigation: X-Ray, Synchrotron and Neutron Techniques. 2013; 7(1):56-61
Dyshekov AA. Generalization of the nonstandard approach in the dynamic theory of diffraction for deformed crystals. Crystallography Reports. 2013; 58(7):984-989