Научная статья на тему 'Variations of structure modeling methods and Raman spectral characteristics for the iodine crystal'

Variations of structure modeling methods and Raman spectral characteristics for the iodine crystal Текст научной статьи по специальности «Химические науки»

CC BY
124
34
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
КРИСТАЛЛИЧЕСКАЯ СТРУКТУРА ЙОДА / ГАЛОГЕННЫЕ СВЯЗИ / ТОПОЛОГИЧЕСКИЙ АНАЛИЗ ЭЛЕКТРОННОЙ ПЛОТНОСТИ / РЕЛЯТИВИСТСКИЙ ЭФФЕКТ / СПЕКТРОСКОПИЯ КОМБИНАЦИОННОГО РАССЕЯНИЯ / IODINE CRYSTAL STRUCTURE / HALOGEN BOND / QUANTUM TOPOLOGICAL ANALYSIS OFELECTRON DENSITY / RELATIVISTIC EFFECT / RAMAN SCATTERING SPECTROSCOPY

Аннотация научной статьи по химическим наукам, автор научной работы — Yushina I.D., Bulatova L.M., Nasibullina S.E., Bartashevich E.V.

On the example of the iodine crystal structure we have selected the optimal basis sets, allowing reproduction of interatomic distances and Raman spectral characteristics by 3D periodic Kohn-Sham calculations. Advantages of two approaches, taking into account the relativistic effect, have been compared: Effective Core Pseudopotentials and the Gaussian type basis set, constructed on the basis of the Douglas-Kroll-Hess approach. It has been shown that the latter approach not only correctly reproduce the experimentally observed geometric parameters of the iodine interactions and characteristics of Raman spectra, but also it reveals the electron density accumulation and depletion in the area of outermost valence shell of an iodine atom. That is directly illustrated by the Laplacian of electron density function.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Variations of structure modeling methods and Raman spectral characteristics for the iodine crystal»

Физическая химия

VARIATIONS OF STRUCTURE MODELING METHODS AND RAMAN SPECTRAL CHARACTERISTICS FOR THE IODINE CRYSTAL

I.D. Yushina, South Ural State University, Chelyabinsk, Russian Federation, idu-xda@mail.ru L.M. Bulatova, South Ural State University, Chelyabinsk, Russian Federation, luciyabulat@yandex. ru

S.E. Nasibullina, South Ural State University, Chelyabinsk, Russian Federation, svetlaya.n.e@ yandex.ru

E. V. Bartashevich, South Ural State University, Chelyabinsk, Russian Federation, kbartash@yandex. ru

On the example of the iodine crystal structure we have selected the optimal basis sets, allowing reproduction of interatomic distances and Raman spectral characteristics by 3D periodic Kohn-Sham calculations. Advantages of two approaches, taking into account the relativistic effect, have been compared: Effective Core Pseudopotentials and the Gaussian type basis set, constructed on the basis of the Douglas-Kroll-Hess approach. It has been shown that the latter approach not only correctly reproduce the experimentally observed geometric parameters of the iodine interactions and characteristics of Raman spectra, but also it reveals the electron density accumulation and depletion in the area of outermost valence shell of an iodine atom. That is directly illustrated by the Laplacian of electron density function.

Keywords: iodine crystal structure, halogen bond, Quantum Topological Analysis of Electron Density, relativistic effect, Raman scattering spectroscopy.

Introduction

Structure-forming non-covalent interactions of halogens Hal...Hal in crystals and solid states in many ways determine thermodynamic, thermophysical, spectral and other properties [1, 2]. It is known that because of anisotropy of electrostatic potential on the Van der Waals surfaces of a molecule, the halogen atom, bound in the molecule, can form two types of non-covalent interactions [3, 4]. The Type I interactions are purely Van der Waals interactions in their nature; they are characterized by random positioning of two covalent bonds of halogens, belonging to different molecules. The Type II interactions are strongly directed. They are noted for mutual orientation of molecules, in which two covalent bonds of halogens are situated at the right angle to each other. In this case the area of electron density accumulation in one halogen atom is directed to the area of electron density depletion in the other atom. This area of electron density depletion is always formed on the extension of the covalent bond of a halogen atom; it is called o-hole [5]. In this area the nucleus is shielded by valence electrons to a lesser degree, and the area of generally positive values of electrostatic potential is formed [6, 7]. This type of non-covalent interactions is called halogen bonds [8, 9].

Halogen bond properties can be successfully studied from the perspective of QTAIMC - Quantum Theory of Atoms in Molecules and Crystals [10], because this theory is aimed at searching for binding interactions, including those which occur among non-covalent ones. Such interactions are characterized by bond paths: two lines, each point of which is different from other neighboring points of space by larger values of electron density. These lines connect atomic nuclei, which are separated by the general interatomic surface, where through the vector gradient of electron density equals zero. The QTAIMC approach includes the topological analysis of electron density p(r), obtained by quantum-chemical calculations, or within precise X-ray diffraction experiment. In its turn the topological analysis allows us to identify critical points of electron density. The electron density and its properties in the bond critical points are important characteristic values describing chemical bond properties in molecules and crystals.

Halogen bonds in a chlorine crystal were first described from the perspective of electron density distribution properties, obtained through high-resolution X-ray diffraction experiment [11]. The paper

Вестник ЮУрГУ. Серия «Химия». 2015. Т. 7, № 3. С. 25-32

presents the visual presentation of mutually consistent orientation of the exhaustion area on to the area of electron accumulation with the Laplacian of electron density V2p(r). The Laplacian of electron density characterizes three-dimensional curvature of its dropping with increase of distance from the atom nucleus. Depending on the ratio of radial positive curvature of electron density and the curvature in the orthogonal directions, the Laplacian will alternate in signs. Alternating minimums and maximums of V2p(r) correspond to areas of accumulation V2p(r)<0 and depletion V2p(r)>0 of the electron density around the nucleus, so they show atom electron shells. However, ranges of the Laplacian negative values, which should correspond to outermost electron shells, do not appear at atomic numbers Z > 29, as it is specified in the paper [12]. That is why, as a rule, other functions of electron density, for instance, one-electron potential, are used for the Br and I compounds for the sake of outermost electron shells localization [9, 13].

Another important feature at reproduction of some experimentally observed properties of halogens with high atomic numbers, such as I or At, is the importance of relativistic effect, this is particularly manifested when analyzing properties, connected with outermost electron shells [14, 15]. Together with structural information, obtained from the electron density distribution analysis in the iodine crystal, the special attention shall be paid to its vibrational properties. The most informative method in studying I-I covalent bonds in crystals and solids is Raman spectroscopy [16, 17]. A great number of works deal with the theoretical frequency rates, obtained for structures of isolated polyiodide anions or cation-anions systems [18, 19]. However, this approach does not allow us to take into consideration the influence of crystal environment, effects of intermolecular interactions in solids, that makes it difficult to compare with the experimental spectral data. Such differences can be partly due to significant disturbance of isolated structure geometry of a molecular complex or a cluster in comparison with the crystal structure. The dynamic approach to crystal lattice with consideration of atomic vibrations around equilibrium positions allows us to explain physicochemical crystal properties, connected with thermal effects, phase transitions, conduction properties. Theoretical calculations allow us to deduce vibrational spectra, to forecast crystal structure stability and to obtain thermodynamic properties, such as heat of formation and sublimation, entropy, and others [20].

The purpose of this paper is to select basis sets, allowing to provide the accurate modeling of cova-lent and halogen bonds and their characteristics in the iodine crystal, as well as to calculate wavenumb-ers for vibrations, that are active in Raman spectra. The method should provide data, suitable for the to-pological analysis of electron density in conditions of periodic quantum-chemical calculations, and should reproduce experimental data as accurately as possible. For this purpose several of basis sets for the iodine atom available in literature have been tested in this paper, and possibilities and advantages of the two approaches to relativistic effect consideration have been compared. These two approaches are: usage of Effective Core Pseudopotentials and the Gaussian-type basis set, constructed according to Douglas-Kroll-Hess approach [21, 22]. We make the comparison of calculated geometric characteristics of the iodine crystal structure and of electron density properties, calculated for covalent and halogen bonds, with the data obtained on the basis of the high resolution X-ray diffraction experiment [23]. Theoretical vibrational characteristics are compared to the experimental polarized Raman spectra for the single crystal of iodine [24].

Calculations

As a part of the study the periodic calculations of the wave function in the iodine crystal were made using the program CRYSTAL14, by the Kohn-Sham method (B3LYP) and various basis sets, shown in Table 1. Two groups of basis sets were tested. One group included Stuttgart fully-relativistic energy-consistent pseudopotentials ECP-mdf28 and ECP-mdf46 including 28 and 46 core electrons, respectively [25]. The pseudopotential ECP-mdf28 could be attributed to the group of potentials with the small core, and ECP-mdf46 had the member of the group with the large core. When using the pseudopotential ECP-mdf46 only outer-shell electrons 5 s2 h 5p5 were included into the valence part. The valence part in both cases was described by the three-time split basis sets of the VTZ type. The other group was represented by the DZVP basis set, including 14 shells [26] and its analogue DZVPmod including 11 shells [27]. In the DZVPmod basis set the sp-type hybrid shells were used to describe internal electron levels, so that up to 8 electrons that could be situated existed in that shell. Accounting of relativistic effect was implemented with the help of the DZPDKH basis set in terms of the Douglas-Kroll-Hess ap-

proach [21, 22]. The basis set was obtained through optimization of standard Gaussian function coefficients [28, 29], approximating equations, obtained from the electronic part of Dirac Hamiltonian [30].

Iodine crystal structure optimization was made for all atoms of the irreducible cell part with fixed cell parameters. Allowed atom coordinate variations were only those, which did not cause changes in the crystal symmetry. The Hessian matrix was calculated for the found optimal atom configuration at the Г point in the center of the Brillouin zone. Vibration frequencies were calculated at the Г point in the harmonic approximation. On the basis of the obtained data the total Raman intensities were calculated for the single crystal of iodine.

Results and Discussion

Mutual arrangement of molecules with halogen bonds in the iodine crystal and with the interatomic distances, stated in the paper [23], is shown in Fig. 1. It has been found that in the case of using fully-relativistic core pseudopotentials ECP-mdf28 and the valence part, described by the VTZ basis set, such geometric characteristics, as the covalent and halogen bond lengths, are reproduced most accurately: ARi.i = 0.04 A, ARl i = -0.03 A, where AR = R calculated - R experimental. According to the data, represented in Table 1, other basis sets also show satisfactory results when localizing equilibrium geometry in the crystal. In all considered cases the positive values of ARI-I have been observed, besides, the I...I halogen bond lengths have been underestimated in most of the cases.

Fig. 1. Mutual arrangement of molecules with halogen bonds in the iodine crystal

Table 1

Distances (A), the values of the electron density in the bond critical point (a.u.) in the crystal structure of iodine, optimized in different basis sets

Basis set I(1)-I(2) Covalent bond I(1)...I(3) Halogen bond I(1)...I(1) Van der Waals interactions

Rcalculation Р(гь) Rcalculation Р(гь) Rcalculation Р(гь)

ECP-mdf46 VTZ 2.876 0.045 3.709 0.011 - -

ECP-mdf28 VTZ 2.755 0.059 3.474 0.018 3.977 0.008

DZVPmod 2.791 0.062 3.473 0.018 3.940 0.009

DZVP 2.802 0.061 3.467 0.019 3.939 0.009

DZPDKH 2.812 0.060 3.414 0.022 4.060 0.007

Experimental [23] Rexp =2.718 p(rb)exp =0.050 Rexp =3.501 p(rb)exp =0.015 Rexp =3.980 p(rb)exp =0.009

The results of the theoretical topological analysis of electron density have been compared to the data from the paper [23]. There the adjusted values of electron density restored on the basis of the extended multipole modeling Hansen and Coppens approach [31] have been analyzed for I-I covalent bonds and the strongest I.I noncovalent interactions. The application of basis sets with core pseudopo-

Вестник ЮУрГУ. Серия «Химия». 2015. Т. 7, № 3. С. 11-18

tentials leads to underestimation (ECP-mdf46 VTZ) and overestimation (ECP-mdf28 VTZ) of electron density in the I-I covalent bond critical points. The DZPDKH basis set provides a significantly underestimated halogen bond length and, consequently, an overestimated electron density in the halogen bond critical points Ap(rb) = 0.007 atomic units. However, the observed range of Ap(rb) values both for cova-lent and halogen bonds may be considered to be reasonable.

The following characteristic features have been observed while analyzing the Laplacian of electron density distribution in the crystal plane, containing halogen bonds between diiodine molecules. When using the DZVP basis sets with core pseudopotentials, the covalent bond area and the space between neighbouring molecules are described in the similar way. In both cases the covalent bond areas have similar level of observed details. However, the Laplacian of electron density in Fig. 2a,b does not show the electron density accumulation in the outermost electron shell area of iodine. This fact evidences low informative value of the V2p(r) function when describing iodine ability to form halogen bonds. On the contrary, the outline map of the electron density Laplacian, obtained through the DZPDKH basis set and indicated in Figure 2c demonstrates the electron density accumulation in the iodine atom equatorial area and exhaustion formed on the continuation of the covalent bond. Consequently, the DZPDKH basis set appears to be the only one of the studied sets, which demonstrates the principle of halogen bond formation: orientation of the area of electron density accumulation of one atom up on the area of electron density exhaustion of the other. Figure 2c shows the I-I covalent bond area in greater details (contour lines of 0.002, 0.004 and 0.008 atomic units) and it also has the contour line of 0.002 atomic units in the halogen bond area unlike Fig. 2a and 2b.

a) b) c)

Fig. 2. The contour lines of the Laplacian of the electron density, accompanied by bond paths in iodine crystal in a variety of basis sets: in fully-relativistic pseudopotential basis set type ECP-28mdf VTZ (a), in the DZVP basis set (b)

and DZPDKH (c)

According to the data of experimental polarized spectra [24], two lines are observed for the iodine crystal. These lines correspond to the in-phase and out-of-phase valence symmetrical vibrations. In the first case, the change of lengths of neighboring molecules is in coordination, while in the second case, it is out of phase: the stretch of the covalent bond in one molecule corresponds to the contraction in the neighboring molecule. Symmetry of obtained in-phase Ag and out-of-phase B2g vibrations corresponds to the experimental data from the polarized spectra [24] (Table 2).

Table 2

The wavenumbers of valence vibrations in the crystal of iodine calculated using a variety of basis sets

Basis set и (Ag) I-I, cm"1 и (B2g) I-I, cm-1

ECP-mdf46 VTZ 166.5 177.5

ECP-mdf28 VTZ 196.5 202.2

DZVPmod 178.6 187.4

DZVP 180.7 188.5

DZPDKH 180.5 188.8

Experimental [25] 180 189

mod

100

DZVP

ЕСР mdf46-VTZ \

\ V

о

150 160 170 180 190 200 210

Wavenumber. cm

-l

Fig. 3. Calculated unpolarized Raman spectra of iodine crystal optimized in various basis sets

On the basis of obtained theoretical wavenumbers for stretching vibrations in the iodine crystal it is possible to select a group of preferable basis sets, reproducing the data of experimental polarized spectra within the error of 2 cm-1. They are DZVPmod, DZVP and DZPDKH. The application of the core pseudopotential with a big core ECP-mdf46 leads to substantial decrease of obtained wavenumbers, and application of the pseudopotential with a small core ECP-mdf28 leads to underestimation, where deviations from the experimental value in these cases are comparable.

Display of the calculated integrated Raman spectra (Fig. 3) allows us to visually assess the ratio of in-phase and out-of-phase vibrations and overlapping of the lines, corresponding to them. Maximum resolution of lines is attained through the use of the basis set ECP-mdf46 VTZ, the difference between stretching peaks equals 13 cm-1. The greatest overlapping of lines is observed in the basis set ECP-mdf28 VTZ, where the difference between stretching peaks equals 6 cm-1. Basis sets of the DZVP group demonstrate roughly the same difference between wavenumbers of in-phase and out-of-phase vibrations; however, in case of the DZPDKH basis the out-of-phase vibration is characterized by relatively lower intensity.

As can be seen from the above, basis sets with fully-relativistic core pseudopotentials slightly underestimate electron density in the critical point of the covalent bond. However, lack of overlapping in grade lines in the area linking to the core makes it vulnerable to apply basis sets with core pseudopotentials in the electron density topological analysis tasks.

Calculated wavenumbers of iodine vibrations in the crystal, which have been obtained when analyzing the basis sets, show that the application of core pseudopotentials ECP-mdf46 and ECP-mdf28 leads to significant deviations from experimentally observed values. Other basis sets reproduce experimentally observed wavenumbers within the accuracy up to ±2 cm"1, only underestimating the difference between in-phase and out-of-phase vibrations a little.

The Laplacian of electron density, obtained with the use of the DZPDKH basis set, demonstrates accumulation of electron density in the area of the outermost valence shell of iodine atoms. That is why the DZPDKH basis set, where the relativistic effect on the basis of the Douglas-Kroll-Hess approach is taken into consideration, is the most suitable for modeling and analyzing electron density topological properties in crystals of iodine containing compounds with halogen bonds.

The work was supported by the Russian Ministry for Education and Science GZ729, and was carried out on a supercomputer "TORNADO" SUSU.

Conclusion

Вестник ЮУрГУ. Серия «Химия». 2015. Т. 7, № 3. С. 11-18

References

1. Kupper F.C., Feiters M.C., Olofsson B., Kaiho T., Yanagida S., Zimmermann M.B., Carpenter L.J., Luther G.W., Lu Z., Jonsson M., Kloo L. [Commemorating Two Centuries of Iodine Research: An Interdisciplinary Overview of Current Research]. Angew. Chem., Int. Ed., 2011, vol. 50, pp. 1159811620.

2. Svensson P. H., Kloo L. [A Vibrational Spectroscopic, Structural and Quantum Chemical Study of the Triiodide Ion]. Journal of the Chemical Society, Dalton Trans, 2000, pp. 2449-2455.

3. Sakurai T., Sundaralingam M., Jeffrey G.A. [A Nuclear Quadrupole Resonance and X-ray Study of the Crystal Structure of 2,5-dichloroaniline]. Acta Crystallogr., 1963, vol. 16, pp. 354-363.

4. Desiraju G.R. [Crystal Engineering: the Design of Organic Solids]. J. Appl. Cryst., 1991, vol. 16, pp. 265.

5. Politzer P., Murray J.S., Clark T. [Halogen Bonding: an Electrostatically-driven Highly Directional Noncovalent Interaction]. Phys. Chem. Chem. Phys., 2010, vol. 12, pp. 7748-7757.

6. Politzer P., Riley K.E., Bulat F.A., Murray J.S. [Perspectives on Halogen Bonding and Other o-hole Interactions: Lex Parsimoniae (Occam's Razor)]. Comput. Theor. Chem., 2012, pp. 2-8.

7. Clark T., Hennemann M., Murray J.S., Politzer P. [Halogen Bonding: the o-Hole]. J. Mol. Model., 2007, vol. 13, pp. 291-296.

8. Desiraju G.R. Ho P.S., Kloo L., Legon A.C., Marquardt R., Metrangolo P., Politzer A. P., Res-nati G., Rissanen K. [Definition of the Halogen Bond (IUPAC Recommendations 2013)]. Pure Appl. Chem., 2013, vol. 85, pp. 1711-1713.

9. Bartashevich E.V., Yushina I.D., Stash A.I., Tsirelson V.G. [Halogen Bonding and Other Iodine Interactions in Crystals of Dihydrothiazolo(Oxazino)Quinolinium Oligoiodides from the Electron-Density Viewpoint]. Cryst Growth Des., 2014, vol. 14, pp. 5674-5684.

10. R.F.W. Bader, Atoms in Molecules. A Quantum Theory, Oxford University Press, New York, 1990, 532 p.

11. Tsirelson V.G., Zou P.F., Tang T.H., Bader R. [Topological Definition of Crystal Structure: Determination of the Bonded Interactions in Solid Molecular Chlorine]. Acta Crystallogr., 1995, vol. 51, pp.143-153.

12. Kohout M., Savin A., Preuss H. [Contribution to the Electron Distribution Analysis. I. Shell Structure of Atoms. J. Chem. Phys, 1991, vol. 95, pp. 1928-1929.

13. Hunter G. [Conditional Probability Amplitudes in Wave Mechanics]. Int. J. Quant. Chem., 1975, vol. 9, pp. 237-238.

14. Silvi B., Savin A., Causa M. [Classification of Chemical B based on Topological Analysis of Electron Localization Functions]. Nature, 1994, vol. 371, pp. 683-686.

15. Pilme J., Renault E., Ayed T., Montavon G., Galland N. J. [Introducing the ELF Topological Analysis in the Field of Quasirelativistic Quantum Calculations]. Chem. Theor. Comput., 2012, vol. 8, pp. 2985-2990.

16. Svensson P. H., Kloo L. [Synthesis, Structure, and Bonding in Polyiodide and Metal Iodide-Iodine Systems]. Chemical Reviews, 2003, vol. 103, pp. 1649-1684.

17. Deplano P., Ferraro J. R., Mercuri M. L., Trogu E. F. Structural and Raman Spectroscopic Studies as Complementary Tools in Elucidating the Nature of the Bonding in Polyiodides and in Donor-I2 Adducts. Coordination Chemistry Reviews, 1999, vol. 188, pp. 71-95.

18. Al-Hashimi N. A., Hussein Y. H. A. [Ab Initio Study on the Formation of Triiodide CT Complex from the Reaction of Iodine with 2,3-Diaminopyridine]. Spectrochimica Acta Part A, 2010, vol. 75, pp.198-202.

19. Otsuka M., Mori H., Kikuchi H., Takano K. [Density Functional Theory Calculations of Iodine Cluster Anions: Structures, Chemical Bonding Nature, and Vibrational Spectra]. Computational and Theoretical Chemistry, 2011, vol. 973, pp. 69-75.

20. Matta C.F., Boyd R.J. [The Quantum Theory of Atoms in Molecules. From Solid State to DNA and Drug Design]. Wiley-VCH Verlag GmbH & Co. KGaA, 2007, pp. 527.

21. Douglas M., Kroll N.M. [Quantum Electrodynamical Corrections to Fine-structure of Gelium]. Ann Phys, 1974, vol. 82, pp. 89-155.

22. Hess B.A. [Relativistic Electronic-structure Calculations Employing a Two-component No-pair Formalism with External-field Projection Operators]. Phys. Rev., 1986, vol. 33, pp. 3742-3748.

23. Bertolotti F., Tsirelson V. G. Shishkina A. V., Forni A., Gervasio G., Stash A. I. [The Intermolecular Bonding Features in Solid Iodine]. Crystal Growth & Design, 2014, pp. 1-20.

24. Congeduti A., Nardone M., Postorino P. [Polarized Raman Spectra of a Single Crystal of Iodine]. Chemical Physics, 2000, vol. 256, pp. 117-123.

25. Peterson K.A. et al. [On the Spectroscopic and Thermochemical Properties of ClO, BrO, IO, and Their Anions]. J. Phys. Chem, 2006. vol. 110, pp. 13877-13878.

26. Godbout N., Salahub D. R. et al. [Optimization of Gaussian-type Basis Sets for Local Spin Density Functional Calculations. Part I. Boron Through Neon, Optimization Technique and Validation]. Can. J. Chem, 1992, vol. 70, pp. 560-562.

27. http://www.tcm.phy.cam.ac.uk/~mdt26/basis_sets/I_basis.txt

28. Jorge F.E., Canal Neto A., Camiletti G.G., Machado S.F. [Contracted Gaussian Basis Sets for Douglas-Kroll-Hess Calculations: Estimating Scalar Relativistic Effects of Some Atomic and Molecular Properties]. J. Chem. Phys., 2009, vol. 130, pp. 064108.

29. Barros C.L., Jorge F.E., Canal Neto A., Campos M. [Gaussian Basis Set of Double Zeta Quality for Atoms Rb Through Xe: Application in Non-relativistic and Relativistic Calculations of Atomic and Molecular Properties]. Mol. Phys., 2010, vol. 108, pp. 1965-1972.

30. Reiher M. [Relativistic Douglas-Kroll-Hess theory]. Wiley Interdisciplinary Reviews: Computational Molecular Science, 2012, vol. 2, pp. 139-149.

31. Hansen N.K., Coppens P. [Testing Aspherical Atom Refinements on Small-molecule Data Sets]. Acta Cryst, 1978, vol. 34, pp. 909-921.

Received 23 May 2015

УДК 541.2, 544.18

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

ВАРИАЦИИ МЕТОДОВ МОДЕЛИРОВАНИЯ ЭЛЕКТРОННЫХ И СПЕКТРАЛЬНЫХ СВОЙСТВ В КРИСТАЛЛИЧЕСКОЙ СТРУКТУРЕ ЙОДА

И.Д. Юшина, Л.М. Булатова, С.Э. Насибуллина, Е.В. Барташевич

Южно-Уральский государственный университет, г. Челябинск

На примере кристаллической структуры йода произведен подбор оптимальных базисных наборов, позволяющих в условиях периодических расчетов воспроизводить межъядерные расстояния, распределение электронной плотности и спектральные свойства йодсодержащих молекулярных кристаллов. Сопоставлены преимущества двух подходов к учету эффектов релятивизма: эффективные остовные псевдопотенциалы и базисный набор гауссового типа, конструируемый на основе методологии Дугласа-Кролла-Гесса. Показано, что последний позволяет корректно воспроизвести не только экспериментально наблюдаемые геометрические параметры кристалла йода и характеристические колебания в спектрах комбинационного рассеяния, но и накопление электронной плотности в области внешней валентной оболочки атома йода, что наглядно иллюстрирует функция лапласиана электронной плотности.

Ключевые слова: кристаллическая структура йода, галогенные связи, топологический анализ электронной плотности, релятивистский эффект, спектроскопия комбинационного рассеяния.

Вестник ЮУрГУ. Серия «Химия». 2015. Т. 7, № 3. С. 11-18

Юшина Ирина Дмитриевна - аспирант, кафедра экологии и природопользования, ЮжноУральский государственный университет. 454080, г. Челябинск, пр. им. В.И. Ленина, 76. E-mail: idu-xda@mail.ru

Булатова Люция Марсельевна - студентка химического факультета, Южно-Уральский государственный университет. 454080, г. Челябинск, пр. им. В.И. Ленина, 76. E-mail: luciyabu-lat@yandex.ru

Насибуллина Светлана Эдуардовна - студентка химического факультета, ЮжноУральский государственный университет. 454080, г. Челябинск, пр. им. В.И. Ленина, 76. E-mail: svetlaya.n.e@yandex.ru

Барташевич Екатерина Владимировна - кандидат химических наук, доцент, кафедра органической химии, Южно-Уральский государственный университет. 454080, г. Челябинск, пр. им. В.И. Ленина, 76. E-mail: kbartash@yandex.ru

Поступила в редакцию 23 мая 2015 г.

ОБРАЗЕЦ ЦИТИРОВАНИЯ

FOR CITATION

Variations of structure modeling methods and raman spectral characteristics for the iodine crystal / I.D. Yushina, L.M. Bulatova, S.E. Nasibullina, E.V. Bartashevich // Вестник ЮУрГУ. Серия «Химия». - 2015. - Т. 7, № 3. -С. 25-32.

Yushina I.D., Bulatova L.M., Nasibullina S.E., Bartashevich E.V. Variations of Structure Modeling Methods and Raman Spectral Characteristics for the Iodine Crystal. Bulletin of the South Ural State University. Ser. Chemistry. 2015, vol. 7, no. 3, pp. 25-32.

i Надоели баннеры? Вы всегда можете отключить рекламу.