The subdivision of forces, acting on the atoms in crystals, on harmonic and anharmonic, where the latter presents the second approximation in the expansion of crystal lattice potential energy in series relative to the atomic displacements, was undertaken in works of Debye  and Born  . Detailed examination of anharmonic forces action in dielectric crystals was carried out in the works of Boguslawski    , Frenkel   Braunbek   , and Brillouin   . In the works of Boguslawski and Frenkel the anharmonic energy arises in the exponent of statistical distribution, which is expanded in power under the premise, that this energy is small and only first power can be taken into account. In some phenomena (thermal expansion) it is so, but for considering of gliding such approximation is insufficient. Also for constructing of adiabatic curve the behavior of anharmonic vibrator must be considered more precise. Braunbek integrates the dynamical equation of motion under the action of anharmonic force, and Brilluin calculates thermodynamically the influence of temperature on the elastic moduls and occurrence the phenomena of melting, when the modul turns to zero.
In the Born works   is asserted, that all distinctions of real crystals properties from that of ideal (harmonic approximation) are caused by anharmonic forces. An approximate formula for calculation of anharmonic forces actions is proposed. The development of such calculations is presented in works   . The kinetic theory of strength was developed in the work  . In the work  the Boguslawski model of anharmonic vibrator was applied to consideration of electric discharge of a dielectric crystal, when an electric field is applied. But when the force is of other nature, the calculations will be almost the same; only the comprehension of the phenomenon will be changed and considered as plasticity or destruction. The phenomenon of plasticity and destruction was considered in some works; here the model of anharmonic vibrator is used for describing these phenomena. The calculations are carried out according to dynamical methodic (further development of calculations that use the virial theorem), which has in solid corps an advantage thanks to the lesser number of assumptions, which can be seen from comparing of work  with the works of Boguslawski. It can be also marked, that in Boguslawski calculations the positive and negative crystal directions are not equivalent, what is essential for arising of thermal expansion and pyroelectricity, but for description of transition to gliding motion the both directions of axis may be equivalent. In presented work the potential energy of the atom contains second and fourth powers of coordinate. For calculation of vibrator behavior dynamical procedure is used   . Another example of the action of anharmonic forces is adiabatic stretching or compression, which causes change of temperature. Works in this direction raised long ago     , but in crystals these effects remained feebly expressed. Modern investigations use high pressures, created also by shock waves, generated in nuclear explosions  . In such conditions the compression is considerable, and structures conversions in crystals can be observed. Approximately the process can be considered as adiabatic, therefore to investigation of adiabatic compression of condensed corps some works are devoted  -  . The equation of adiabatic process for Boguslawski vibrator is deduced in this paper.
2. Dynamical Consideration of Transition to Glide
A symmetrical in relation to both directions of coordinate axis of molecular forces anharmonic vibrator, that potential energy as a function of coordinate is
(f-applied force), will be considered. The graph of , (f is absent) is presented on Figure 1.
The coordinates of maxima have the values
The free energy as a function of mean displacement s and root mean square amplitude u is presented in the form 
Condition of minimum in relation to leads to equations
If the thermal movement is absent , the dependence of the applied force from the displacement is represented through the equation
The extreme points have the coordinates
The stability part of the curve is placed between extreme points. When the applied force is absent, consequently s = 0, the dependence of temperature from amplitude represents the expression
The coordinates of maximum have the values
Figure 1. Graph of vibrator potential energy.
The states are stabile in the range from zero to abscissa of maximum amplitude. By acting of the force the expression of amplitude through displacement with retaining quadratic terms has according (5) the form
It can be introduced the marking
Equations (4), (5) receive the form
The dependence of force from the displacement for a number of temperatures is presented on Figure 2
The coordinates of force maximum have the values
For the maximum of the amplitude will be found the expression
which increases with rising of temperature and at temperature (9) reaches maximal value.
By rising of temperature aspires to zero, but, as can be seen from (12), the coefficient of expansion turns to infinity, so the solution of Equations (4), (5) must be found exact. The exact expression of amplitude is
Figure 2. Dependence of force from displacement for a number of temperatures.
If we put here for the temperature the value (9), corresponding , the radical will be imaginary for values of displacement different from zero, so this value will be the single solution with (great values of displacement, which can change this conclusion, as is evident from Figure 2, are excluded). Correct limiting value of amplitude is given trough the first term of (12), when we put . Because the calculation is approximate, the conditions of stability are not formulated. Usually, the point of extremity divides the curve in stability and not stability parts. When the force overcomes the maximal value (14), the vibrator begins to glide.
3. The Equation of Adiabatic Process
We deduce the equation of adiabatic curve, using the results of Boguslawski vibrator calculations   . The inner energy as a sum of kinetic and potential energies or their mean values
for Boguslawski vibrator with stretching force f has the form
His entropy is presented as follows
Entropy is constant, so one may put
and regard A as a positive constant, depending from the quantity of entropy. The temperature as a function of amplitude is introduced in (18), and the condition of minimum energy leads to equations
The conditions of stability are the following
From (22) the amplitude will be obtained
and from (20), (22) the temperature
From (21) and (26) follows the expression for the applied force
This dependence is represented in Figure 3
From the condition of maximum
follows the values of corresponding coordinates
The conditions of stability (23), (24) are satisfied in this point, in condition (25) the value (26) of amplitude will be substituted, then this inequality transforms to following
So, the point of maximum is the limiting point of stability, which is realized for smaller values of displacement. Excluding displacement from (27), (29), we receive connection between force and temperature
The graph of this connection is represented in Figure 4.
The coordinates of the curve top are
Figure 3. Dependence of force from displacement.
Figure 4. Dependence between force and temperature.
The stability states are placed to right side of the top, because diminishing of displacement corresponds to increase of temperature. So, the Equation (33) represents the adiabatic curve.
4. Concluding Remarks
Considering of anharmonic vibrator simple models gave possibility to receive the expression for the force that corresponds to the limit of elastic behavior and beginning of gliding. It is calculated the dependence of crystal temperature from the applied force in the condition of thermoisolation. More complicated models can give possibility to examine some other properties of the object. In paper  is considered the model of vibrator with the potential energy, which includes the coordinate in second, third and fourth powers. The graph of energy shows two differences in height maxima. It is evident, that overcoming of the low maximum leads to the lost of elastic properties, and overcoming the high to destruction. The calculation of these phenomena can be carried out, when the term in third power is small and can be neglected in zero approximation. In the first approximation when the applied force is absent the phenomenon of pyroelectricity can be taken into account. It can be of interest the calculating of gliding phenomena considering real structure of crystals with real forces of interactions between atoms, as was undertaken in  , but with use of dynamical calculations. This is a heavy task and it is worthwhile to divide them in steps that can be considered one after another. The behavior of vibrators with more complicated dependence of potential energy from deflection can also be considered.