WJCMP  Vol.6 No.1 , February 2016
The Boguslawski Melting Model
Abstract: The anharmonic vibrator, whose expression of potential energy contains second and third powers of coordinates, is treated on the basis of dynamical procedure, which presents the state of motion by means of mean position and mean amplitude of vibration. The divergent statistical integral comes here not into consideration. The free energy is represented through mean atomic displacement and developed in power series, retaining fourth degree. The graphs show that at certain temperature, the minimum in free energy disappears, and the atom escapes from the potential pit. A simple atomic model that represents this phenomenon is proposed and the influence of model dimension and pressure on melting temperature will be presented.

Received 15 December 2015; accepted 21 February 2016; published 24 February 2016

1. Introduction

A simple atomic model, which describes the phenomenon of melting, was elaborated by Boguslawski [1] [2] . According to this model, the atom is exposed to the action of two forces―attraction force, directed to the equilibrium position, and repulsive force, directed away from it. The total potential energy of an atom as a function of abscissa x of the atom with origin in the position of equilibrium (minimum of energy) is represented by the expression


Here is a constant, which is absent in the original paper, and c, a are elastic and anharmonic coefficients. Figure 1 is the graph of polynom (1).

At low temperature the atom vibrates in the region of minimum. When the temperature is high enough, the atom approaches the top of the energy curve and escapes from the potential pit, which signifies the beginning of melting phenomenon. For determining of melting temperature this author tries to calculate the statistical integral,

Figure 1. Potential energy of anharmonic vibrator.

but, unfortunately, it diverges in infinity for represented law of force. Boguslawski assumes that repulsive force is zero in space outside the top of potential curve, but attractive force increases as second power of abscissa, then the calculations is possible to carry out.

The behavior of potential pit outside the pit could have no meaning, if the atom was restricted in his movement to the pit only, but the statistics does not impose such a restriction.

The statistics permits the existents of a time moment, when the atom makes a jump over the potential barrier. In a regular crystal the movements of atoms are ordered and form in simple case a plane wave, so are mutually determined, and the determination lies in the laws of motion.

Braunbeck [3] solved the dynamical equation of motion of an atom in a potential field, described by sinusoidal function of coordinate, which represents the potential field in the crystal lattice. He established two types of motion―in the pit at low temperatures and above the top at high, which are in relation with solid and liquid states. So the melting phenomenon was demonstrated, but mathematical difficulties are considerable. It is worthwhile to use the form of dynamical equations that is able to calculate the mean thermal displacement of an atom and mean amplitude of vibration without solving time-dependent dynamical equations [4] [5] . The problem contains in solving algebraic equations, which may be carried out in some cases.

It is not excluded that the atoms in the lattice are fable bounded and behave almost independently. Then they can jump one by one over the potential barrier, which is an occasionally event, and the theory of transition must be developed in other direction.

Briefly, the described way of calculation is the following. Abscissa x is represented in the form

where s is the mean value of x and u is the mean value of amplitude of vibration, which is determined by equation

The almost periodical function of time satisfies the conditions and. It will be assumed that, because the vibration function may be represented as a sum of trigonometric functions with incommensurable frequencies. Further, mean value of potential energy may be calculated, and the variation principle for variables s, u leads to the equations of state, where the temperature is a designation for mean kinetics energy


The expression for free energy F has the form


The terms that depend on temperature only, are omitted, because they are not represented in the condition of minimum.

Boguslawski used his model of vibrator also for construction of the theory of pyroelectricity, and this work was mentioned in the article [6] . But for consideration of melting, it is useful to apply the way, elaborated by Ginsburg for considering Segnettelectricity [7] [8] . According to this way one must express the free energy through polarization P, which is proportional to s, develop in power series, and construct graphs for different temperatures, which represent obviously the picture of phase transition.

2. The Picture of Melting Phenomenon

If the variable, that describes the state, is volume, then the conditions for melting are formulated in the work [9] . Volum and atom displacement (or polarizations) are similar, because they happen from place changing of atoms. As is shown in [9] , the energy, arises from the change of volume, may be represented like for Boguslawski vibrator.

The free energy of Boguslawski vibrator is presented through expression


One of two equations of thermal equilibrium is written as


So, the amplitude is


Substituting in (4), we obtain after some calculation


Differentiating with respect to s and equating the derivative to zero, we receive the second equation of thermal equilibrium and a connection between temperature and displacement




which gives an expression of temperature through displacement and with (6) also temperature through amplitude. The dependence is represented on Figure 2, where physical meaning has the part of the graph from origin to maximum.

The dependence T(u2) of the form


is represented through the curve with one maximum, which characterize the limiting values of variables. Further increase of the temperature leads to cease of the existence of state and to escape of the atom from the pit. Following expressions for coordinates at escaping are obtained

Figure 2. Dependence of temperature from displacement.


Here the index m marks the maximum of the curve.

So, the relation of amplitude to the displacement at melting point does not depend from the parameters of vibrator. Also, the coefficient of thermal elongation multiplied with temperature of melting does not contain parameters of vibrator. On may receive


These formulas can be compared with Grüneisen assertions [10] [11] . It may be remarked, that is a thermal elongation of interatomic distance, but not an interatomic distance itself, which cannot be calculated from the model. Because of acting of anharmonic forces, which diminish the restoring force, the quantity is greater, then, namely.

The temperature of minimum on Figure 2 has negative sign, but the same magnitude, as in expression (10).

Definite value of amplitude at melting point in relation to interatomic distance was offered by Lindemann as condition of melting [12] , and this point of view was elaborated in the work of Gilvarry [13] . For Boguslawski vibrator the parameters of movement are obtained from the equations of state.

Further, the expansion will be used


The expression for free energy takes the form of an expansion


We may think that the first term includes the second term and under first term understand both. Graphs for low and high temperatures are represented on Figure 3 and Figure 4.

The abscissa of minimum (thermal equilibrium) is determined through the equation


Figure 3. at low temperature.

Figure 4. at high temperature.

The temperature can be expressed as function of displacement and represented as power series.

By neglecting higher then second power of displacement this equation coincides with (9), so is valid for that part of Figure 2, which is placed between origin and maximum.

For low temperature one may put


The first term is zero and the second is obtained through differentiating of (14). So, coordinates of minimum are given by means of expressions, where are retained terms quadratic in temperature


When the temperature rises, the linear and cubic terms in (13) decreases as negative, the quadratic term decreases as positive, and the curve on the Figure 3 strives to be smooth. When both extremes coincide in a bending point with horizontal tangent, the corresponding temperature will be determined by means of (14) and condition


Excluding from (14) and (17) the displacement, we receive the equation for the temperature



we obtain an equation


So, and temperature and displacement have in the bend point B value


which are not far from the values (10). The free energy at this point has the value


These values may be regarded as parameters of melting. The corresponding graph is shown on Figure 5.

Here this phenomenon of melting is understood as disappearance of stable crystalline state without consideration of the liquid, but liquid state is also state with own properties, so that can hinder the movement of atom. The atom will vibrate during some time about equilibrium position, so that the model of anharmonic vibrator can be used for liquid state. Then the expression (16) for free energy with other values of coefficients (denoted by asterisk) may be an approximation to the free energy of liquid state


The condition of phase equilibrium consists in equality of expressions for free energies of both states

Figure 5. for transition temperature.


The temperature of equilibrium is equal


Nominator and denominator must be either positive or both negative. The denominator contains differences of barriers heights reciprocal values. The second case is probable, because in solid phase the depth of minimum is more significant and the height of barrier greater than in the liquid one. When the vibrations are considered as harmonic, the expression like (24) was obtained in the work [14] .

In the work of Born [15] , devoted to the melting phenomenon, condition of melting is formulated as disappearance of shearing resistance thanks to the action of fluidity. But this opinion is not precise, because fluidity exists, when mechanical action is of long duration. By short duration (firing in the liquid) fluidity is not manifested. When duration of mechanical action is very long, the solid corps also reveals fluidity, such behavior show a glacier. The question about rigidity and fluidity was elaborated in the work of Frenkel [16] .

Brilloin also considers the question of temperature influence on elastic properties of solids, and disappearance of elastically resistance at some temperature is according his opinion a melting phenomenon [17] .

The theory of Grüneisen was précised in our work [5] and about melting may be formulated following assertions. In Figure 6 of the mentioned work a family of isotherms is drawn, that are curves with maximum. The height of maximum (value of maximal applied stretching force) diminishes with increasing temperature till such temperature is achieved, that the top of the curve reaches the abscissa axis. The slop of the tangent to isotherm in the point of intersection with abscissa axis increases to the value, then elasticity disappears. Probably the melting arises? But such opinion is imprudent, it can be only supposed, that the state is now a non equilibrium one and cannot be described by formulated equations. The moving atoms form a new structure, which comes to thermal equilibrium with primary and after a time, which can be much greater, then the period of atomic vibration, the fluid will be developed. In this time rearrangements of atoms between different sites can be realized.

Figure 6. Four atomic model of vibrator.

3. Atomic Model of Anharmonic Vibrator

The Boguslawski model is formal, because the parameters of the model are introduced as quantities determined from experiment. Later Frenkel [18] [19] developed a three atomic model of melting (three identical atoms placed on a straight line), where the middle atom interacts by means of atomic forces with others, that are immobile, and moves between them. It turns out, that in the middle point of interval between right and left atoms is placed a minimum of potential energy, if the atoms are near one to another, and maximum, if they are far. So, by changing of the length of the model phase transition may occur. But in this model the middle atom moves only between the others and cannot escape from the pit. It is worthwhile to replace the right or left atom through two identical atoms placed on ordinate from both sides of abscissa axis on equal distances, as is shown on Figure 6, where are marked the dimensions. Such model from other reasons was considered in the work [20] . Here is assumed, that the interactions, that are considered, are stipulated through the forces of repulsion, which is admissible, if the distances between atoms are small. The potential energy of the middle atom has the form


It will be supposed, that h < l.

Further expansions will be used

The potential energy of moving atom may be represented by the expression


It will be assumed, that the action of right atom in the region of origin is small and may be omitted, so the linear term falls out and following expression may be accepted


The graph is represented on Figure 7.

The deduction is not rigorous, but physically we return to Boguslawski vibrator, where the ordinate axis goes through the maximum of potential energy. The elastic coefficient c is now negative and a positive, so the direction of temperature growing on Figure 2 is now reversed. Physical meaning has now the part of the graph, where displacement diminishes with rising of temperature, beginning from the value c/a (minimum of energy). The general expression for transition temperature rests the same, as before. Using expression (20) for transition temperature, we receive


This formula is not suited for numerical calculations, because the ratio l/h is presented in high power, so small variations of the fraction leads to great changes of temperature,. But it is possible to investigate the influence of length by constant ratio l/h on transition temperature. In experiments the dimension is not considered as independent variable, but such is the stretching or, usually, compressing force. This influence can be calculated, using a simple model. When on two interacting with atomic forces atoms, separated through distance l (molecule), are acting stretching forces f and −f. the potential energy is

Figure 7. Potential curve in four atomic models.

The equilibrium distance in the absence of stretching forces is

When the force has maximal value, corresponding to destruction of molecule (two extreme coincide in a bend point with horizontal tangent), the limiting distance has the value

Corresponding limiting force is

So, by increasing of stretching force from zero value to maximal one the distance between atoms increases from to. Supposing linear dependence between distance and external force, we may put


So, stretching leads to decrease and compressing to increase of the melting temperature according linear law. The latter phenomenon was observed in experiments [21] . More elaborate investigation shows, that connection between pressure and temperature for melting phenomenon is represented by means of non linear function [22] .

In article [23] was pointed out, that the phenomenon of melting occurs, when the moving atom reaches the bend point of potential energy. The bend point in potential energy curve actually leads to peculiarities in behavior of the atom [20] , but gives no possibility to escape from the potential pit.

The consideration of melting phenomenon with the help of a cubic cell model was undertaken in the work [24] . There the melting was treated as escaping from the cell, and as independent variable was regarded the length of the cell edge without introducing of the stretching force. This investigation will be continued.

So, the simple Boguslawski model of anharmonic vibrator can be used to describe in general features the phenomenon of melting on the base of dynamical principles.

4. Concluding Remarks

Boguslawski in many of his works elaborated the theory of solid corps, based on the expression of potential energy of anharmonic vibrator with terms of second and third powers of coordinate. The statistical integral with such expression of energy diverges in infinity, and one cannot be sure that the way of calculation, employed by this author, leads to correct results. So, the received expression for the temperature of transition to the liquid state seems some uncertain. Calculations on the base of dynamical principles lead to utterly definite expression for the temperature of transition. This is so, because dynamical calculations treat the movement of atom in the pit as anharmonic vibrations with finite amplitude and without shift of atom in the infinity. But the calculation is not so simple, as in original paper. An atomic model, which creates a potential of Boguslawski vibrator is considered, and the parameters of Boguslawski expression of potential energy are presented in their dependence from parameters of atomic interaction, which include interatomic distances. So the influence of external pressure on the phenomenon of melting can be included in dynamical calculations.

In what relations to the experimental data are the results of calculation? After eliminating of infinite quantity, the Boguslawski formula for melting temperature could be compared with experimental data for some metals. It turned out that his formula differs from our only unimportant―the coefficient in denominator is marked 2η2, where η is determined from experimental data and approximately has for all metals the same value 1.89, so the coefficient in denominator is 7.14 (in our formula 10.4). This means that a tenth of atoms escaped from regular places when melting begins. Boguslawski gives also values for the ratio c/a, calculated from experimental data, so it is possible to calculate Sm and compare with experimental determined interatomic distances, which are given in the table.

Cite this paper: Kozlovskiy, V. (2016) The Boguslawski Melting Model. World Journal of Condensed Matter Physics, 6, 45-55. doi: 10.4236/wjcmp.2016.61007.

[1]   Boguslawski, S.A. (1915) About the Mechanism of Melting Process. Transactions of Russian Physics—Chemical Society, Physics, 47, 540-545.

[2]   Boguslawski, S.A. (1961) About the Mechanism of Melting Process. Selected Works in Physics, Edition with Notes of Prof. V.K. Sementchenko. State Edition for Physics-Mathematical Literature. Moscow, 99-104.

[3]   Braunbeck, W. (1926) Grating Dynamical Theory of Melting Phenomenon. Physical Magazine, 38, 549-572.

[4]   Kozlovskiy, V.Kh. (1991) Dynamical Equations for Time-Averaged Coordinates in Thermal Equilibrium. Soviet Physics Journal, 34, 353-359.

[5]   Kozlovskiy, V.Kh. (2012) To the Question of Validity Grüneisen Solid State Equation. World Journal of Condensed Matter Physics, 2, 219-227.

[6]   Born, M. and Brody, E. (1921) About the Vibration of Mechanical System with Finite Amplitude and Their Quantization. Physical Magazine, 6, 140-152.

[7]   Ginsburg, V.L. (1945) About Dielectric Properties of Segnettelectrics and Barium Titanat. Journal of Experimental and Theoretical Physics, 15, 739-749.

[8]   Ginsburg, V.L. (1949) The Theory of Segnettelectric Phenomena. Successes in Physical Sciences, 38, 490-525.

[9]   Herzfeld, K.F. and Mayer, M.G. (1934) On the Theory of Fusion. Physical Review, 46, 995-1001.

[10]   Grüneisen, E. (1911) To the Theory of One Atomic Solid Corps. Bulletin of German Physical Society, 13, .836-841.

[11]   Grüneisen, E. (1912) Theory of One Atomic Elements Solid State. Annals of Physics, 39, 257-306.

[12]   Lindemann, F.A. (1910) About the Calculation of Own Molecular Frequencies. Physical Papers, 11, 609-612

[13]   Gilvarry, J.J. (1956) The Lindemann and Grüneisen Laws. Physical Review, 102, 308-316.

[14]   Kozlovskiy, V.Kh. (2003) Vibration-Diffusional Model of Melting. Transactions of VI International Conference of Crystals: Growth, Properties, Real Structure, Application, Alexandrov, 8-12 September 2003, 573-594.

[15]   Born, M. (1939) Thermodynamics of Crystals and Melting. Journal of Chemical Physics, 7, 591-603.

[16]   Frenkel, J.I. (1936) On the Liquid State and the Theory of Fusion. Physical Magazine, 47, 58-65.

[17]   Brillouin, L. (1938) On Thermal Dependence of Elasticity in Solids. Physical Review, 54, 916-917.

[18]   Frenkel, J.I. (1935) The Liquid State and the Theory of Fusion. Physico-Chemical Works of URSS, 3, 632-648.

[19]   Frenkel, J.I. (1958) The Liquid State and the Theory of Fusion. Collected Papers, 2, 269-279.

[20]   Kozlovskiy, V.K. (2010) Model of Crystal Behavior in the Vicinity of the Bend Point of the Potential Curve for an Atom. Russian Physics Journal, 53, 1158-1166.

[21]   Bridgeman, P.W. (1931) The Physics of High Pressure. Chap. VII, George Bell & Sons, London.

[22]   Simon, F.E. (1953) The Melting of Iron at high Pressure. Nature, 172, 746-747.

[23]   Rashevsky, N.V. (1926) To the Theorie of Melting Heats. Physical Magazine, 40, 214-219.

[24]   Kozlovskiy, V.K. (1999) Dynamical Examination of Mechanical Behavior According to Three Atomic and Melting According to Nine Atomic Models of a Crystal. Transactions of IV International Conference Сrystals: Growth, Properties, Real Structure, Application, 2, 466-476.