Quantum-chemical research of endohedral yttrium metallofullerenes Текст научной статьи по специальности «Физика»

УДК 544.137

Quantum-chemical Research of Endohedral Yttrium Metallofullerenes

Anastasia S. Kholtobina* Darya I. Tsyplenkova^ Aleksandr A. Kuzubov* Maxim A. Visotin§ Aleksandr S. Fedorov^

Kirensky Institute of Physics Akademgorodok, 50/38, Krasnoyarsk, 660036 Siberian Federal University Svobodny, 79, Krasnoyarsk, 660041

Russia

Received 01.12.2016, received in revised form 07.03.2017, accepted 26.05.2017 The structural and electronic properties of single molecule Y@C82, their join couple and crystal structure of Y@C82 were investigated by DFT-GGA approach. The calculations show that Y@Ce2 form stable crystal structures wich may have ferroelectric properties, so they can be applied in electronics as a ferroelectric memory.

Keywords: endohedral metallofullerenes, Y@Ce2, ferroelectric properties, density functional theory. DOI: 10.17516/1997-1397-2017-10-4-422-428.

Introduction

Endohedral fullerenes are an interesting class of fullerenes because electron transfer from encaged metal atom to carbon cage has been known to occur and this oftentimes alters the electronic and magnetic properties of the fullerenes [1]. Particularly, endohedral yttrium-fullerenes have been obtained by several research groups [2,3]. In 1995, Takata and co-workers performed synchrotron X-ray studies on a powder of Y@C82 to confirm the endohedral nature of EMFs for the first time [3]. However, it is still not clear whether the sample contained a pure Y@C82 isomer or if it was a mixture of two or more regioisomers. Nowadays, the structures of new EMFs can be routinely predicted from the first principles (knowing only a formula of the molecule) with the high reliability rivaling that of singlecrystal X-ray diffraction studies [4]. Potential applications of endodedral fullerenes were also predicted on the basis of their peculiar electronic, physical and chemical properties, including superconductors [5], metallofullerene lasers [5] ferroelectric materials [6,7], nanomemory devices [5] quantum computers, etc. [8].

* akholtobina93@mail.ru

tasteralin@gmail.com

talexxkuzubov@gmail.com

§ visotin.maxim@gmail.com ^alex99@ iph.krasn.ru © Siberian Federal University. All rights reserved

The structural and electronic properties of endohedral yttrium-fullerenes were investigated using software package OpenMX (Open source package for Material eXplorer), designed for nano-scale material simulations based on density functional theories (DFT) norm-conserving pseudopotentials, and pseudo-atomic localized basis functions. Simulations of single complex of Y@C82, structure, consisting of two complexes of Y@C82, and crystal of Y@C82 (face centered cubic cell with a=11.457A) [9] were performed. Potential barriers of yttrium transitions were then calculated using NEB method. To search a minimum energy path (MEP) in geometrical phase space connecting two stable structures, a nudged elastic band (NEB) method based on is supported in OpenMX. The Monkhorst-Pack [10] k-point Brillouin sampling was used. The k-point grid contained 1 x 1 x 1 points along a, b and c directions, respectively

2. Results and discussion

2.1. Yttrium positions in fullerene cage

Since a fullerene molecule has pillow-like shape, an yttrium atom occupies several unequal states into fullerene cavity. All these states, namely high, lowl and low2, that mean the top part of molecule, left and right on the bottom part of molecule, respectively, were investigated (see Fig. 1B). It was found, both low states have the lowest bond energies (-5.154 and 5.152 eV,

Fig. 1. (A) Direction of the dipole moment for the state 'low', (B) transitions and corresponding energy differenses (eV) of Y atom jump between states 'high'-'low', 'low'-'high', 'low'-'low'; (C) directions of the applied electric field: t - along 'high'-'low' states transition, n - along the dipole moment direction

respectively), which were calculated by using formula 1, and it indicates their equivalence, while bond energy for state high is -4.788 eV. Moreover, the dipole moment was calculated for systems high and low (2.638 and 2.859 D, respectively).

1. Methods

1.282-■>

•4— 1.282

C

Ebond E complex Eful EY

(1)

where Ecompiex is total energy of endohedral yttrium-fullerene Y@C82, Eful is total energy of fullerene, EY is energy of yttrium atom.

e kT

P = ^ Aji (2)

Ee kT i

where AEi is difference between total energies of endohedral yttrium-fullerene Y@C82 with distinct states of yttrium atom. The population probabilities, calculated by using formula (2), are 2.6 x 10-7 and ~ 0.5 for high and low, respectively. Since the population probability high is short, the polarizability of the system will be defined by the dipole moment of low state.

Furthermore, the yttrium atom transition barriers between all potential states, namely from high to low, from low to high, from lowl to low2, were studied and values of barriers are 0.230, 0.596 and 1.282 eV, respectively. Consequently, it is more possible, that yttrium atom transition will be performed from state high to low along 111 direction in crystal.

2.2. Structure, consisting of two complexes of Y@C82

Simulation of endohedral yttrium-fullerenes with all potential states of yttrium atom, signed as high1-high2 (h1h2), low1right-high2 (l1rh2), low1right-low2left (l1r_l2l), low1left-low2right (l1l_l2r), low1left-low2left (l1l_l2l), were performed for the structure, consisting of two complexes of Y@C82, which are arranged as in crystal (see Fig. 2). Among all researched structures the lowest bond energy was found for l1r_l2l geometry (-5.186 eV). Fig. 2 indicates the barriers of all potential transitions of yttrium atom in these structures.

The bond energies, the bond energies of molecule couple (the analogue of crystal lattice energy for structures, consisting of two complexes of Y@C82) (see Tab. 1) show that the lowest metal bond energy belongs to the geometry with low states of yttrium atom as in case of single molecule. Obviously, the yttrium atom charge insignificantly changes, so only displacement of yttrium atom influences on the dipole moment.

Table 1. The metal bond energies, the bond energies of molecule couple and yttrium atom charges for structures consisting of two complexes of Y@C82

Structure Ebond Me (eV) Ebond 2mol (eV) Charge of atom Y (a. u.)

h1h2 -4.800 -0.034 0.428

l1r h2 -5.007 -0.083 0.307

l1r l2l -5.186 -0.076 0.305

l1l l2r -5.155 -0.012 0.307

l1l l2l -5.158 -0.017 0.304

2.3. Crystal structure of Y@C82

Moreover, the calculations, connected with crystal structure of endohedral yttrium-fullerenes, were executed by analogy with the previous one. In particular, different states of yttrium atom (high, low1, and low2) were investigated and, it was found, that the most advantageous yttrium atom state is low (see Tab. 2). In addition, the crystal lattice energies were calculated for all

Fig. 2. The potential barriers (eV) of yttrium atom transitions for structures consisting of two Y@C82 molecules

crystal structures using formula 3.

Ecrystal _lat Ecrystal nEcomple

(3)

where Ecrystai is total energy of crystal model, n is number of molecules Y@C82 in crystal model, Ecomplex is total energy of Y@C82

In addition, a correlation was obtained between the dipole moment and the electric field and the correlation between the yttrium atom charge and the electric field. Furthermore, the electric field was applied along two directions (see Fig. 1C), the first of them runs along the transition from state high to one low, and the second one runs along the dipole moment direction. It was found, the first and the second directions correspond to the complex and linear dependences (see Fig. 3). With respect to the yttrium atom charge, electric field correlation in the case of electric field applied along the dipole moment direction, charge stops increasing and remains constant in contradistinction to case of electric field applied along the transition from state high to low.

The research of transition barriers from state high to low (0.661 eV) and from lowl to low2 (1.229) has shown that they insignificantly differ from analog transition barriers in sin-

gle molecules (see Fig. 1). It is interesting to note that in this case synchronous displacement of yttrium atoms descents in all endohedral yttrium-fullerenes. Also, the values of the transition barriers are comparable to those of common ferroelectric materials (e.g. 0.44 eV in potassium dihydrogen phosphate [11] ), thus, such fullerene crystal structures can be considered as a promic-ing material for ferroelectric applications. Furthermore, the crystal model of Y@C82 comprises four endohedral yttrium-fullerenes, which are arranged as in the face-centered lattice. The displacement of one yttrium atom from state high to low was investigated in this model (see Fig. 4). The yttrium atom transition barrier between potential states, namely from A to C, was studied and value of barrier are 1.127 eV.

Table 2. The metal bond energies, the crystal lattice energy and yttrium atom charges for crystals of Y@C82

Structure Ebond Me (eV) Ecryst (eV) Charge of atom Y (a. u.)

high -5.240 -0.452 0.424

lowl -5.610 -0.458 0.295

low2 -5.579 -0.425 0.299

Fig. 3. The dipole moment electric field (a) and Y atom charge electric field dependencies (b) for different directions of applied field: t is along 'high'-'low' states transition, n is along the dipole moment direction

Conclusion

Endohedral yttrium-fullerenes Y@C82 simulations result in formation of stable structures like Y@C82, both in the form of a single molecule and a crystal. The most favorable positions of Y atom inside the carbon cage are similar in the each type of structures, confirming the common trend in all researched models. Furthermore, the transitions barriers also indicate that Y atom tends to occupy the same positions. It is interesting to note that ferroelectric properties were found in crystal structure. The complex and linear dependences were obtained for electric field applied along the transition from state high to low and for electric field applied along the

Fig. 4. The geometries of crystal model Y@C82 with different yttrium atom positions in fullerene

cage

dipole moment direction, respectively. Researching of the yttrium atom charge - electric field correlation indicated that in the case of electric field applied along the dipole moment direction charge stops increasing and remains constant. Endohedral yttrium-fullerenes Y@C82 can be applied in radiotronics and automatics and the development of piezoelectric devices, condensers and temperature gages.

This work was supported by the Russian Foundation for Basic Research (no. 16-43-242148), Krasnoyarsk Regional Science Foundation (agreement 35/16 dated 18.11.2016) and by the president of Russia Scientific School program NSh 7559.2016.2.

References

[1] S.Yung, C.-R.Wang, Word Scientific, 2014.

[2] M.M.Ross, H.H.Nelson, J.H.Callahan, S.W.McElvany, J. Phys. Chem., 96(1992), 5231.

[3] S.W.McElvany, J. Phys. Chem., 96(1992), 4935.

[4] M.Takata, B.Umeda, E.Nishibori, M.Sakata, Y.Saitot, M.Ohno, H.Shinohara, Nature, 377(1995), 46-49.

[5] L.Bao, C.Pan, Z.Slanina, F.Uhlik, T.Akasaka, X.Lu, Angew. Chem. Int. Ed., 55(2016), 1-6.

[6] D.S.Bethune, R.D.Johnson, J.R.Salem, M.S.Devries, C.S.Yannoni, Nature, 366(1993), 123.

[7] J.J.Cioslowski, Am. Chem. Soc., 113(1992), 4139.

[8] A.A.Popov, S.Yang, L.Dunsch, Chemical reviews, 113(2013), 989-6113.

[9] H.Suematsu, Y.Murakami, H.Kawata, Y.Fujii, N.Hamaya, O.Shimomura, K.Kikuchi, Y.Achiba, I.Ikemoto, Mat. Res. Soc. Symp. Proc., 349(1994), 01994.

[10] H.J.Monkhorst, J.D.Pack, Phys. Rev. B, 13(1976), 5188.

[11] B.A.Strukov, I.V.Shaidshtein, T.V.Pavlovskaya, S.V.Grabovskii, Y.Uesu, M.Fukunaga, L.Carman, Ferroelectrics, 267(2002), 329.

Квантово-химическое исследование эндоэдральных иттриевых металлофуллеренов

Анастасия С. Холтобина Дарья И. Цыпленкова Александр А. Кузубов Максим А. Высотин Александр С. Федоров

Институт физики им. Л. В. Киренского СО РАН Академгородок, 50/38, Красноярск, 660036

Сибирский федеральный университет Свободный, 79, Красноярск, 660041

Россия

Структурные и электронные свойства одиночного комплекса Y@Cs2, структуры, состоящей из двух комплексов Y@C82, а также данного комплекса в периодических условиях были исследованы в рамках DFT-GGA-подхода. Моделирование эндоэдральных металлофуллеренов Y@C82 показало возможность формирования из них стабильных кристаллических структур, которые обладают сегнетоэлектрическими свойствами, что может быть применено в электронике и сегнетоэлек-трической памяти.

Ключевые слова: эндоэдральные металлофуллерены Y@C82, сегнетоэлектрические свойства, метод функционала плотности.