CC BY
66
Физика
ELABORATION OF ATOMIC MODEL FOR AB INITIO CALCULATION OF THE FERRITE / CEMENTITE INTERFACE1
A.V. Verkhovykh2, K.Yu. Okishev3, A.A. Mirzoev4, D.A. Mirzaev5
The work describes constructing a supercell for ab initio modelling of ferrite / cementite interface in pearlite microsctructure of iron-carbon alloys. Baga-ryatskii type orientation relationship between the two phases was accepted, with habit plane (101)c || (112)a. Optimization of ab initio modelling parameters was also performed using WIEN2k package.
Keywords: ferrite, cementite, interface, Bagaryatskii relationship, ab initio modelling, supercell.
Introduction
It is well known that solubility of carbon in the low-temperature bcc polymorphous modification of iron is very low, so that almost all carbon in commercial iron alloys (steels) usually containing 0,5 to 5 at. % C is present in the form of iron carbide Fe3C (cementite) co-existing with almost pure alpha iron (ferrite). In slowly cooled steels cementite is a part of a structural constituent called pearlite, a mixture of thin lamellae of ferrite and cementite [1]. Average interlamellar spacing of pearlite formed isothermally from prior fcc gamma solid solution (austenite) depends on transformation temperature, ranging from 0,5...0,3 ^m when supercooling below the eutectoid temperature Ai = 1000 K is AT = 30...50 K to ~ 0,1 ^m and even smaller if austenite is supercooled to about 830 K (AT ~ 170 K). Interlamellar spacing of pearlite significantly affects hardness and strength of steels [2].
Ferrite / cementite interfaces also play a role in trapping dissolved atoms of hydrogen. Hydrogen is one of important components in iron and steels. It has the smallest atomic radius among interstitial impurities (0,46 A according to [3], which is confirmed by ab initio calculations [4]) and high diffusivity but low solubility in iron. From technical point of view hydrogen is a harmful impurity causing embrittlement, cracking etc. [5-7]. One of the ways of suppressing harmful effects of hydrogen is creating structural traps, i.e. positions where hydrogen atoms are bound and incapable of moving on the lattice [6-8]. One type of such traps are interphase interfaces, including a / Fe3C [6, 9].
Beside direct experimental studies, an important tool for quantitative analysis of hydrogen trapping is ab initio modelling. Previously the authors carried out such modelling of hydrogen interaction with substitutional impurity atoms [10] and vacancies [11, 12]. Ab initio modelling of hydrogen interaction with a / Fe3C interfaces requires determination of optimal modelling parameters and constructing a supercell describing the interface. This is the purpose of the present paper.
Atomic model
It is well known that cementite has an orthorhombic lattice of Pnma space group. Unit cell with room-temperature parameters a = 4,524; b = 5,088; c = 6,741 A [13] contains 12 iron and 4 carbon atoms. Iron atoms occupy two crystallographically nonequivalent positions [14]: G (“general”):
1 The work was supported by the Ministry of Education and Science of the Russian Federation and by the Russian Foundation of Basic Research (grants no. 13-02-00048 and 14-03-00618).
2 Verkhovykh Anastasia Vladimirovna is Post-graduate Student, General and Theoretical Physics Department, South Ural State University. E-mail: ursaeva@physics.susu.ac.ru
3 Okishev Konstantin Yur’evich is Dr. Sc. (Physics and Mathematics), Professor, Physical Metallurgy and Solid State Physics Department, South Ural State University.
E-mail: okishev@physmet.susu.ac.ru
4 Mirzoev Aleksandr Aminulaevich is Dr.Sc. (Physics and Mathematics), Professor, General and Theoretical Physics Department, South Ural State University.
E-mail: mirzoev@physics.susu.ac.ru
5 Mirzaev Dzhalal Aminulovich is Dr. Sc. (Physics and Mathematics), Professor, Physical Metallurgy and Solid State Physics Department, South Ural State University.
E-mail: mirzayev@physmet.susu.ac.ru
± [[ х, у, z ]]; ± [[ 1/2 - х, 1/2 + y, z ]]; ± [[ x, у, 1/2 - z ]]; ± [[ 1/2 - х, 1/2 + у, 1/2 - z ]] (1а)
and S (“special”):
± [[ u, v, 1/4 ]]; ± [[ 1/2 - и, 1/2 + v, 1/4 ]]. (16)
Carbon atoms located in the centres of “normal” prismatic interstices [15] have coordinates
± [[ r, s, 1/4 ]]; ± [[ 1/2 - r, 1/2 + s, 1/4 ]]. (1c)
According to [16], parameters in eqs. (1) are x = 0,333; у = 0,175; z = 0,065; и = -0,167; v = 0,040; r=0,43; s = -0,13.
There are two main types of orientation relationships (OR) between ferrite and cementite in pearlite, namely the Bagaryatskii (Isaichev) and Pitsch OR. For pearlite formed at small supercoolings fractions of colonies having each of the two OR’s are approximately equal, while with decreasing transformation temperature the probability of Bagaryatskii OR increases [17]. We will then concentrate our analysis on the Bagaryatskii type of OR. It supposes that the following planes are parallel in both phases:
(100)c || (0T1)a; (010)c || (1TT)a; (001)c || (211)a, (2)
and the habit plane of lamellae is (101)c || (112)a [17-19].
Fig. 1 shows cementite lattice projected on the (010) plane. Planes (101) in the iron sublattice of
cementite form 12-layer stacking with average interplanar spacing of т2^/a 2 + c2 = 0,677 A (fig. 1).
There are two types of these planes containing G and S atoms. It is most probable that the (101)c plane at the coherent ferrite / cementite interface will contain S type atoms because atomic pattern in S planes is more similar to that of (112)a planes [20].
For adequate description of macroscopic amounts of phases by ab initio modelling, periodic boundary conditions are to be set: the supercell in which electronic structure is calculated is supposed to be translated infinitely along all three axes. In order to provide integrity of the interface under these boundary conditions the interface has to be parallel to one of the faces of the supercell. In our case the coordinate system is to be changed in such a way that the habit plane (101)c || (112)a become parallel to one of coordinate planes. The supercell should consist of two parts, cementite and ferrite, each containing sufficient number of atomic planes to consider interaction of neighbouring interfaces insignificant.
Thus a reasonable way to construct the cementite part of the supercell is to change the three base vectors [a, 0, 0], [0, b, 0], [0, 0, c] to [-a, 0, c], [0, b, 0], [a, 0, c], see fig. 2. The unit cell of cementite built up on new base vectors has the volume twice greater than the initial one and contains 32 atoms. The new cell is monoclinic; its parameters are ai = 8,118 A; bi = 5,088 A; ci = 8,118 A; Pi = 2 atan(a/c) = 67,73°. It is also convenient to shift the origin of coordinates of the new cell along the initial х axis by (1/4 - 0,333) = -0,083 a in order to place the first S plane at z1 = 0 (see fig. 2).
The distance between the first and the last S planes in this cell (planes 0 and 13 in fig. 1) is
1,083 c1 = 8,792 A. The distance from the last cementite S plane located at z1 = 1,083 to the first ferrite plane is to be chosen so that the average distance between iron atoms in these two planes be as close as possible to the interatomic distance in a iron aa^bj2 = 2,482 A, where aa = 2,866 A is a iron lattice parameter. This condition is fulfilled if the distance between the two planes is 0,820 A (i.e. 0,109 c1) and the initial atom in the first ferrite plane is located at x1 = 0,368 и y1 = 1/4 being shifted from the S atom [[0,167 a1; 0,040 b1; 1,083 c1]] of the last cementite plane by [0,201 a1; 0,210 b1; 0,109 c1]. Distances between iron atoms in the two planes are then 2,38...2,69 A.
Interplanar distance for subsequent ferrite planes is to be aa/ V6 = 1,170 A. Iron atoms in each next (112)a plane are shifted from the previous one by a vector 6aa [5 IT]; in new coordinates it corresponds
to [-4 a1 - d cos в , 1 b1, (d/ sin P^q] = [0,195 a1; 1 b1; 0,156 c1].
Initial atom of the last ferrite plane is to be shifted from the initial S atom of the next (translated along z1) supercell by [-0,201 a1; -0,210 b1; -0,109 c1]. Since this S atom has coordinates x1 = 0,250 and y1 = 0,540, then the number of ferrite planes in the supercell is to be taken such as to make the coordinates of the atom of the last plane close to x1 = 0,049 and y1 = 0,330. To provide this, seven ferrite planes are to be inserted, which yields the coordinates of the initial atom of the seventh plane x1 = 0,040 and y1 = 0,250. If the distance between neighbouring ferrite and cementite planes is again taken 0,886 A, then the distances between iron atoms of the two planes are 2,43...2,65 A.
With the added seven ferrite planes the Oxzx size of the supercell increases from c = 8,118 A to c2 = 18,150 A, i.e. 2,236 times. The supercell now contains 64 atoms: 56 iron (28 in both ferrite and cementite parts) and 8 carbon (in the cementite part). It is shown in fig. 3.
13 12 11 10 9 8 7 6 5 4 3 2 1 0
Fig. 1. Cementite lattice projected along [010]. Traces of (101) planes are shown and numbered. Iron atoms in S positions are drawn brighter; smaller circles are carbon atoms in the centres of “normal” prismatic interstices
Fig. 2. Cementite lattice projected along [010] with base vectors of the initial orthorhombic and having double volume monoclinic unit cells. Crosses indicate atoms included in the cementite part of the supercell
Fig. 3. Supercell containing both cementite and ferrite parts
Choice of ab initio modelling parameters
The above described supercell was used as initial structure for ab initio computer modelling. All calculations were performed by WIEN2k code [21] implementing linearized augmented plane wave (LAPW) method. Generalized gradient approximation (GGA) was chosen for calculations.
Supercells with monoclinic angles significantly different from 90° (fi\ = 67,73°) are inconvenient for ab initio modelling in the WIEN2k. So the above described supercell was rebuilt by the program: translation periods changed to a3 = 8,118; b3 = 5,088; c3 = 16,842 A with fi3 = 94,221° and a3 = y3 = 90°.
Then optimal modelling parameters were defined. In the LAPW method the unit cell is divided into non-overlapping atomic spheres (muffin-tin spheres) having radius Rmt, and the interstitial region. It is assumed that wave functions within muffin-tin spheres are combinations of spherical harmonics, and the charge density is thus spherically symmetric. Wave functions in the interstitial region are expanded into plane waves with the Kmax vector. Its value is one of the principal parameters that affect calculation accuracy because it determines the number of basis functions (matrix size). Another important parameter is the number of ^-points. This value characterizes the discrete mesh used for numerical integration over the Brillouin zone. The higher it is, the more points are taken into account in the integration, the greater the accuracy, and the greater computer costs. The purpose of choosing optimal parameters is to provide high calculation accuracy (0,01 eV ~ 1 mRy) at acceptable consumption of computer resources.
First optimization was performed for the number of k-points, since this parameter is independent of
the others. Calculation parameters were the following: R^ = 2,0 a.u. [22]; Rnt = 1,4 a.u.; Kmax = 5,0
[22]; cut-off energy Ecut = 7,0 Ry. As can be seen from fig. 4, when the number of k-points exceeds 64, the total energy E changes not greater than by 1 mRy, so 64 k-points may be accepted as the optimal value.
Optimal value of Rmt for iron atoms is R^ = 2,0 a.u. [22]. Dependence of the system energy calculated using previously optimized parameters on Rt is shown in fig. 5. As can be seen, at R^ > 1,45 a.u.
the total energy changes by less than 0,5 mRy, so the value R^ = 1,45 a.u. was taken as optimal. Parameter Rmt Kmax controlling convergence is then equal to 7,25.
Thus, the optimal parameters are following: number of k-points Nk = 64; parameter controlling the
convergence Rmt Kmax = 7,25; radii of atomic spheres R^ = 2,0 a.u. and = 1,45 a.u.
Number of ^-points Rmt, a. u.
Fig. 4. Dependence of the total energy of the system Fig. 5. Dependence of the total energy of the system
on the number of ^-points on the muffin-tin sphere radius of carbon
Resume
The supercell for ab initio modelling of the ferrite / cementite interface with Bagaryatskii orientation relationship and (101)c || (112)a habit plane has been built up and optimal values of parameters necessary for ab initio modelling have been determined. After additional volume optimization this supercell can be used for calculation of ferrite / cementite interfacial energy, interaction of the interface with impurities (e.g., hydrogen) and other lattice defects etc.
References
1. Smirnov M.A., Schastlivtsev V.M., Zhuravlev L.G. Osnovy termicheskoy obrabotki stali (Basics of Heat Treatment of Steel). Ekaterinburg, Ural Branch RAS Publ., 1999. 496 p. (in Russ.).
2. Pickering F.B. Physical Metallurgy and the Design of Steels (Materials Science Series). Applied Science Publishers Ltd., 1978. 275 p.
3. Shtremel’ M.A. Prochnost’ splavov. Ch. II: Deformatsiya (Strength of Alloys. Pt. II: Deformation). Moscow, MISiS Publ., 1997. 527 p. (in Russ.).
4. Rakitin M.S., Mirzoev A.A., Mirzayev D.A. Bulletin of the South Ural State University. Series “Metallurgy”. 2010. no. 13(189). Issue 14. pp. 67-71.
5. Gel’d P.V., Ryabov R.A. Vodorod v metallakh i splavakh (Hydrogen in Metals and Alloys). Moscow, Metallurgiya Publ., 1974. 272 p. (in Russ.).
6. Hirth J.P. Effects of Hydrogen on the Properties of Iron and Steel. Metallurgical Transactions. 1980. Vol. 11A, no. 6. pp. 861-890. DOI: 10.1007/BF02654700.
7. Kolachev B.A. Vodorodnaya khrupkost’ metallov (Hydrogen Brittleness of Metals). Moscow, Metallurgiya Publ., 1985. 218 p. (in Russ.).
8. Bernstein I.M., Pressouyre G.M. The Role of Traps in the Microstructural Control of Hydrogen Embrittlement of Steels. Hydrogen Degradation of Ferrous Alloys (ed. by R.A. Oriani, J.P. Hirth and M. Smialowski). Park Ridge, N.J.: Noyes Publications, 1985, pp. 641-685.
9. Takai K., Nozue A. Effect of Isothermal Transformation Temperature and Cold-Drawing on Hydrogen Occlusion Behavior in Eutectoid Steel. J. Japan Inst. Metals, 2000. Vol. 64, no. 8. pp. 669-676.
10. Mirzaev D.A., Mirzoev A.A., Okishev K.Yu., Rakitin M.S. Theory of Hydrogen Solubility in Binary Iron Alloys Based on Ab Initio Calculation Results. Molecular Physics. 2012. Vol. 110, no. 1112. pp. 1299-1304. DOI: 10.1080/00268976.2011.645895.
11. Ursaeva A.V., Mirzoev A.A., Ruzanova G.E. Bulletin of the South Ural State University. Series “Metallurgy”, 2011. Issue 17. no. 36 (253), pp. 51-55. (in Russ.).
12. Mirzaev D.A., Mirzoev A.A., Okishev K.Yu., Verkhovykh A.V. Hydrogen-Vacancy Interaction in BCC Iron: Ab Initio Calculations and Thermodynamics. Molecular Physics. 2014. DOI: 10.1080/00268976.2013.861087.
13. Andrews K.W., Dyson D.J., Keown S.R. Interpretation of Electron Diffraction Patterns. London, 1968.188 p.
14. Fasiska E.J., Jeffrey G.A. On the Cementite Structure. Acta Cryst. 1965. Vol. 19, pt. 3. pp. 463471. DOI: 10.1107/S0365110X65003602.
15. Schastlivtsev V.M., Yakovleva I.L., Mirzaev D.A., Okishev K.Yu. On the Possible Positions of Carbon Atoms in the Cementite Lattice. The Physics of Metals and Metallography. 2003. Vol. 96, no. 3. pp.313-320.
16. Andrews K.W. The Structure of Cementite and Its Relation to Ferrite. Acta Metallurgica, 1963. Vol. 11, no. 8. pp. 939-946. DOI: 10.1016/0001-6160(63)90063-4.
17. Schastlivtsev V.M., Mirzaev D.A., Yakovleva I.L. Struktura termicheski obrabotannoy stali (The Structure of Heat Treated Steel). Moscow, Metallurgiya, 1994. 288 p. (in Russ.).
18. Tushinskiy L.I., Bataev A.A., Tikhomirova L.B. Struktura perlita i konstruktivnaya prochnost’ stali (Pearlite Structure and the Structural Strength of Steel). Novosibirsk, VO “Nauka” Publ., 1993. 280 p. (in Russ.).
19. Kurdyumov G.V., Utevskiy L.M., Entin R.I. Prevrashcheniya v zheleze i stali (Transformations in Iron and Steel). Moscow, Nauka Publ., 1977. 240 p. (in Russ.).
20. Okishev K.Yu., Mirzaev D.A., Verkhovykh A.V. Bulletin of the South Ural State University. Series “Metallurgy”. 2013. Vol. 13, no. 2, pp. 95-102. (in Russ.).
21. Schwarz K., Blaha P. Solid State Calculations Using WIEN2k. Computational Materials
Science. 2003. Vol. 28, no. 2. pp. 259-273. DOI: 10.1016/S0927-0256(03)00112-5.
22. Ursaeva A.V., Ruzanova G.E., Mirzoev A.A. Vybor optimal'nykh parametrov dlya postroeniya
maksimal'no tochnoy modeli OTsK-zheleza (Selection of Optimal Parameters for Formation the Most Accurate Model of Bcc Iron). Bulletin of the South Ural State University. Series “Mathematics. Me-
chanics. Physics”. 2010. No. 9 (185). Issue 2. pp. 97-101. (in Russ.).
Received 3 March 2014
Bulletin of the South Ural State University Series “Mathematics. Mechanics. Physics” ________________2014, vol. 6, no. 2, pp. 49-55
УДК 538.915 + 669.112.227.322
ПОСТРОЕНИЕ АТОМИСТИЧЕСКОЙ МОДЕЛИ ФАЗОВОЙ ГРАНИЦЫ ПЕРЛИТ/ЦЕМЕНТИТ ДЛЯ AB INITIO МОДЕЛИРОВАНИЯ
А.В. Верховых1, К.Ю. Окишев2, А.А. Мирзоев3, Д.А. Мирзаев4
Работа описывает построение атомной структуры поверхности раздела феррита / цементит в системе Fe-C предназначенную для проведения первопринципных расчетов параметров фазовой границы методом суперячейки для железоуглеродистых сплавов. Для согласования соотношений между решетками двух фаз был использован были использованы ориентационные соотношения типа Багарятского, при которых плоскости (101)с (112)а параллельны между собой. Проведена оптимизация параметров первопринципного моделирования для точного воспроизведения параметров решетки феррита с использованием пакета WIEN2k.
Ключевые слова: феррит, цементит, поверхность, соотношение Багарятского, неэмпирическое моделирование, суперячейка.
Литература
1. Смирнов, M.A. Основы термической обработки стали / M.A. Смирнов, В.М. Счастливцев, Л.Г. Журавлев. - Екатеринбург, УрО РАН, 1999. - 496 с.
2. Pickering, F.B. Physical Metallurgy and the Design of Steels (Materials Science Series) / F.B. Pickering. - Applied Science Publishers Ltd., 1978. - 275 p.
3. Штремель M.A. Прочность сплавов. В 2 ч. Ч. 2: Деформация / M.A. Штремель. - М.: МИ-СиС, 1997. - 527 с.
4. Ракитин, М.С. Изменение электронной структуры альфа-железа, содержащего внедренные атомы водорода / М.С. Ракитин, А.А. Мирзоев, Д.А. Мирзаев // Вестник Южно-Уральского государственного университета. Серия «Mеталлургия». - 2010. - Вып. 14. - № 13(189). - С. 67-71.
5. Гельд, П.В. Водород в металлах и сплавах / П.В. Гельд, РА. Рябов. - М.: Mеталлургия, 1974. - 272 с.
6. Hirth, J.P. Effects of Hydrogen on the Properties of Iron and Steel. Metallurgical Transactions. -1980. - Vol. 11A, no. 6. - P. 861-890. doi: 10.1007/BF02654700.
7. Кюлачев, БА. Водородная хрупкость металлов / БА. Кюлачев. - M.: Mеталлургия, 1985. -218 с.
8. Bernstein, I.M. The Role of Traps in the Microstructural Control of Hydrogen Embrittlement of Steels. Hydrogen Degradation of Ferrous Alloys / I.M. Bernstein, G.M. Pressouyre. - Park Ridge, N.J.: Noyes Publications, 1985. - P. 641-685.
9. Takai, K. Effect of Isothermal Transformation Temperature and Cold-Drawing on Hydrogen Occlusion Behavior in Eutectoid Steel / K. Takai, A. Nozue. - J. Japan Inst. Metals. - 2000. - Vol. 64, no.8. - P. 669-676.
10. Theory of Hydrogen Solubility in Binary Iron Alloys Based on Ab Initio Calculation Results / D.A. Mirzaev, A.A. Mirzoev, K.Yu. Okishev, M.S. Rakitin. Molecular Physics. 2012. Vol. 110, no. 11-
12. - P. 1299-1304. doi: 10.1080/00268976.2011.645895.
11. Урсаева, А.В. Ab initio моделирование комплексов водород-вакансия (VHn) в ОЦК-железе / А.В. Урсаева, А.А. Мирзоев, Г.Е. Рузанова // Bестник Южно-Уральского государственного университета. Серия «Металлургия». - 2011. - Вып. 17. - № 36 (253). - С. 51-55.
1 Верховых Анастасия Владимировна - аспирант, кафедра общей и теоретической физики, Южно-Уральский государственный университет.
E-mail: ursaeva@physics.susu.ac.ru
2 Окишев Константин Юрьевич - профессор, доктор физико-математических наук, кафедра физического металловедения и физики твердого тела, Южно-Уральский государственный университет.
E-mail: okishev@physmet.susu.ac .ru
3 Мирзоев Александр Аминулаевич - профессор, доктор физико-математических наук, кафедра общей и теоретической физики, Южно-Уральский государственный университет.
E-mail: mirzoev@physics.susu.ac.ru
4 Мирзаев Джалал Аминулович - профессор, доктор физико-математических наук, кафедра физического металловедения и физики твёрдого тела, Южно-Уральский государственный университет.
E-mail: mirzayev@physmet.susu.ac.ru
12. Hydrogen-Vacancy Interaction in bcc Iron: Ab Initio Calculations and Thermodynamics / D.A. Mirzaev, A.A. Mirzoev, K.Yu. Okishev, A.V. Verkhovykh // Molecular Physics: сб. науч. тр., 2014. doi: 10.1080/00268976.2013.861087.
13. Andrews, K.W. Interpretation of Electron Diffraction Patterns / K.W. Andrews, D.J. Dyson,
S.R. Keown. London, 1968. - 188 p.
14. Fasiska, E.J. On the Cementite Structure / E.J. Fasiska, G.A. Jeffrey // Acta Cryst. - 1965. -Vol. 19, pt. 3. - P. 463-471. doi: 10.1107/S0365110X65003602.
15. On the Possible Positions of Carbon Atoms in the Cementite Lattice / V.M. Schastlivtsev, I.L. Yakovleva, D.A. Mirzaev, K.Yu. Okishev // The Physics of Metals and Metallography. - 2003. -Vol. 96, no. 3. - P. 313-320.
16. Andrews, K.W. The Structure of Cementite and Its Relation to Ferrite / K.W. Andrews. - Acta Metallurgica. - 1963. - Vol. 11, no. 8. - P. 939-946. doi: 10.1016/0001-6160(63)90063-4.
17. Счастливцев, В.М. Структура термически обработанной стали / В.М. Счастливцев, Д.А. Mирзаев, И.Л. Яковлева. - M.: Металлургия, 1994. - 288 с.
18. Турчинский, Л.И. Структура перлита и конструктивная прочность стали / Л.И. Турчинский, A.A. Батаев, Л.Б. ^хомирова. - Новосибирск, ВО «Наука», 1993. - 280 с.
19. ^рдюмов, Г.В. Превращения в железе и стали / Г.В. Kурдюмов, Л.М. Утевский, Р.И. Энтин. - M.: Наука, 1977. - 240 с.
20. Oкишев, ^Ю. Возможность захвата водорода на границе феррит/цементит в сталях / ^Ю. Oкишев, Д.А. Mирзаев, A^. Верховых // Вестник Южно-Уральского государственного университета. Серия «Металлургия». - 2013. - Т. 13, № 2. - С. 95-102.
21. Schwarz, K. Solid State Calculations Using WIEN2k / K. Schwarz, P. Blaha // Computational Materials Science. - 2003. - Vol. 28, no. 2. - P. 259-273. doi: 10.1016/S0927-0256(03)00112-5.
22. Урсаева, А.В. Выбор оптимальных параметров для построения максимально точной модели ОЦК-железа / А.В. Урсаева, Г.Е. Рузанова, А.А. Мирзоев // Bестник Южно-Уральского государственного университета. Серия «Математика. Механика. Физика». - 2010. - Вып. 2. - № 9 (185). - С. 97-101.
Поступила в редакцию 3 марта 2014 г.
