CC BY
30
OPTO-IT
OPTICAL PROPERTIES OF LOWEST-ENERGY CARBON ALLOTROPES FROM FIRST-PRINCIPLES CALCULATIONS
V.A. Saleev1, A.V. Shipilova1 1 Samara National Research University, Samara, Russia
Abstract
We study optical properties of lowest-energy carbon allotropes in the infrared, visible and ultraviolet spectral ranges in the general gradient approximation of the density functional theory. In our calculations we use an all-electron approach as well as a pseudo-potential approximation. In the infrared range, complex dielectric functions, infrared and Raman spectra have been calculated using a CRYSTAL14 program. Electronic properties and energy-dependent dielectric functions in the visible and ultraviolet spectral ranges are calculated using a VASP program. We describe with good accuracy the experimentally known optical properties of a cubic diamond crystal. Using the obtained set of relevant calculation parameters, we predict the optical constants, dielectric functions and Raman spectra of the lowest-energy hypothetical carbon allotropes and lonsdaleite.
Keywords: optical properties, Raman spectrum, first-principles calculations, density functional theory, crystal structure, carbon allotropes.
Citation: Saleev VA, Shipilova AV. Optical properties of lowest-energy carbon allotropes from first-principles calculations. Computer Optics 2017; 41(4): 476-483. DOI: 10.18287/2412-61792017-41-4-476-483.
Acknowledgments: The work was partially funded by the Ministry of Education and Science of the Russian Federation under the Competitiveness Enhancement Program of Samara University for 20132020, project 3.5093.2017/8.9.
Introduction
Diamond and different forms of carbon materials are the subject of the intensive theoretical and experimental study [1, 2]. It is well known that carbon can form sp2 and sp3 hybridized bonds which are realized under ambient conditions in cubic diamond, graphite, fullerene and graphene [3]. At high temperature and pressure graphite can be converted to cubic diamond or to the lonsdaleite (2H hexagonal diamond) [4, 5]. For a long time, hexagonal diamond has been produced artificially by static and shock wave compression of well-crystallized graphites [6]. Recently it was shown that hexagonal diamond can be also obtained from cubic diamond [7]. It was found that graphite cold-compression leads to a creation of the new sp3-bonded stable forms of carbon allotropes [8]. Nowadays, there are about two hundreds of different predicted hypothetical sp3 carbon allotropes, which are collected in the SACADA database [9]. The small difference of energy of several these allotropes relatively to the diamond, 0.01 - 0.10 eV per atom, raises the idea of a possibility to find these allotropes in the mixed carbon phases. The experimental search of such carbon allotropes should be based on some new physical information about formation of these new phases. The most probable signals can be connected with optical properties of the materials, such as Raman spectra, optical coefficients in the regions of infrared (IR), visible and ultraviolet (UV) spectra, as well as the energy dependence of the absorption and refractive indices. Here we calculate different optical properties of the six lowest-energy sp3 carbon allotropes: cubic diamond, lonsdaleite [4], 4H-diamond [10], SiC12 [11], C28 [12] and mtn [13]. The main goal of our study is to find the quantitative level of difference between optical properties of cubic diamond and relevant carbon allotropes, which can be measured experi-
mentally. We start with a prediction of Raman shift spectra and IR active mode spectra, than we will study energy-dependent complex dielectric functions and their derivative in the ranges of visible and ultraviolet light. Our calculations are based on density-functional-theory (DFT) methods [14, 15] as it is implemented in CRYSTAL14 [16] and VASP [17] program packages. The calculations for the IR range have been done using CRYSTAL14 program, which uses the all-electron approach with atomic orbital basis sets. The relevant calculations in the visible and ultraviolet ranges of light have been done using VASP program package which uses plane-wave basis sets and the pseudo-potential approach.
Computational methods and details
Raman and IR spectra of crystal structure are defined by the set of harmonic phonon frequencies at the r point which can be obtained from the diagonalization of the Hessian matrix of the second derivatives with respect to atomic displacements [18]:
Hai,bj =
1
4MMb
d2 E du du.
J
where uai and ubj are displacements of atoms a and b in the reference cell along the i-th and j-th Cartesian directions, respectively.
The Raman intensity of the Stokes line of a phonon mode Qp, active due to the ap component of the polariza-bility tensor a, is expressed as follows:
ip ~(da„ /dQp )2. (1)
The scheme of calculation, recently implemented in the CRYSTAL14 program [16], explores second-order Coupled-
Perturbed Hartree-Fock/Kohn-Sham (CPHF/KS) equations [19]. The Raman spectrum is then computed by considering the transverse optical (TO) modes and by adopting a pseudoVoigt functional form: a linear combination of a Lorentzian and a Gaussian curve with full width at half maximum of 8 cm-1. Raman intensities are normalized so that the largest value is conventionally set to 100 a.u.
To calculate IR spectra we should know the complex dielectric tensor £,,(v) which is computed for each inequivalent polarization direction on the basis of a classical Drude-Lorentz model:
e„ (n) = Ep,„ + X fPA / (np -v2 - ivgp), (2)
p
where ii indicates the polarization direction, eopta is the optical dielectric tensor, vv, f and gp are the TO frequency, oscillator strength and damping factor for the p-th vibration mode, respectively. The real and imaginary parts of e,i(v) are computed and the maxima of this function correspond to the TO frequencies. The optical or high-frequency dielectric tensor is computed in a quasi-free electron approximation via coupled perturbed Hartree-Fock (Kone-Sham) method [20]. The refractive (n) and absorption (k) indices are computed as real and imaginary
parts of the complex refractive index n*(v) = Ve , also for each inequivalent polarization direction.
In the projector augmented plane wave method realized in VASP the frequency-dependent dielectric functions is obtained in the random phase approximation, where the imaginary part of the frequency-dependent dielectric tensor is written as
e® (w) = (4p2e2/W)lim(1/q2) X 2wk x
q®0 c vk
q cv (3)
xd(eck - evk - W) (uck+eaa I Uvk) (uck+e,q I Uvk) *
where the indices c and v refer to conduction and valence band states respectively in the sum over the empty states, ec,vk are the corresponding eigenenergies, W is the volume of a primitive cell, k-point weights wk are defined such that they sum to 1, ea,b are the unit vectors for the three Cartesian directions, and uck is the cell-periodic part of the orbitals at point k. The real part of the dielectric tensor e^ (w) is obtained by the Kramers-Kronig transformation
2 p
m/ N 2 f- £(„2R (w')w'
eab(w) = 1+-J dw,
Jo w'2 -w2 + ih
(4)
where P denotes the principal value. By cubic symmetry, the following relation is satisfied for diamond and mtn
(5)
e(1'2) = e(1'2) = e(1'2),
e^ = 0, a^ß,
so the real and imaginary parts of the complex dielectric constant e = E1+/E2 can be determined by e1 2 = s .j . In
the case of the other, anisotropic, structures we use the average values:
Ej,2 =E1; = (1/3)(eXX:2) +E^2) +E^).
(6)
Computational setup, CRYSTAL14
We use the two schemes of DFT calculations, both are implemented in the CRYSTAL14 program [16]. The
first one incorporates the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional [21] at general gradient approximation (GGA) of DFT and POB-TZVP all-electron basis set [22], and the second one - the B3LYP hybrid functional [23] and adopted Pople's 6-21G Gaussian all-electron basis set [24, 25]. The level of accuracy of calculating the energies of Coulomb and Hartree-Fock exchanges is controlled by a set of TOLINTEG parameters, which were chosen as {8, 8, 8, 8, 18}. The convergence threshold on energy for the self-consistent-field (SCF) calculations is 10-7 Hartree for structural optimization and 10-8 Hartree for vibration frequency calculations. The number of basis vectors in the irreducible Brillouin zone is given by the shrink parameter IS = 8 for structural optimization and IS = 16 for vibration frequency calculations. The relaxation of cell parameters and atomic positions to equilibrium values was carried out until the lattice stress became less than 0.02 GPa.
Computational setup, VASP
The computations of the complex dielectric functions in the visible and ultraviolet regions were performed by the VASP package [17] at the level of GGA and with two functionals: the above-mentioned PBE [21] and the screened Heyd, Scuseria, and Ernzerhof (HSE06) hybrid functional [26], since the latter is well-suited to reproduce the electronic properties including the band gaps of a number of elemental and binary insulators and semiconductors. For this purpose, at first we performed the relaxation of the discussed structures with a cutoff energy of 700 eV for the plane-wave basis set, until the energy was converged up 10-7 eV per unit cell and the residual stress was less than 10-5 eV per atom. The integration over the Brillouin zone was performed using Monkhorst-Pack grids with the number of k-points along each direction of the reciprocal cell equal to an integer divisor of 25A over the length of corresponding lattice vector. The chosen set of input parameters gives the reliable results for the main physical properties of diamond as it was shown in our previous work [12].
Structural properties
The full geometry optimization of studied structures have been performed using the quasi-Newton algorithm in CRYSTAL14 program. The relaxed lattice constants and atomic Wyckoff positions are presented in the Tables 6 and 7. The differences between lattice constants obtained in two schemes of calculations are about 1 %. While the energy differences per atom relatively diamond are equal, see Tables 1 and 2. To check a mechanical stability of hypothetical allotropes under study, we calculated the elastic constants and phonon spectra which are collected in the Fig. 4. Our results for structural properties agree with previous DFT calculations [10 - 13].
Raman and IR spectra
Raman identification of new carbon allotropes needs a precise calculation of vibration spectra. We performed these calculations applying quasi-harmonic approximation at the r point as it was done for Raman identification of
lonsdaleite in Refs. [27, 28] and for Raman spectroscopy of nanocrystalline diamond in Ref. [29]. We have found for cubic diamond that the position of single Raman peak, corresponding to the first-order scattering of F2g symmetry, is 1297 cm-1 (PBE, POB-TZVP ) and 1332 cm-1 (B3LYP, 6-21G). The last one is in a best agreement with experimental value. Such a way, we perform calculations of Raman and IR spectra here and after using B3LYP exchange-correlation functional and modified Pople's 6 - 21 basis set [25]. The presented in Fig. 1 polycrystalline (powder) Raman spectra were computed by averaging over the all possible orientations of the crystallites.
Table 1. Energy difference per atom relatively diamond, density and bulk modulus of the allotrope (PBE, POB-TZVP)
Table 2. Energy difference per atom relatively diamond, density and bulk modulus of the allotrope (B3LYP, 6-21G)
For the lonsdaleite we predict three Raman active modes: E:g (1336), A:g (1312) and E2g (1209). Up to now, the pure monocrystalline lonsdaleite have not been found or synthesized and there is a well-known problem of correct diagnostics of lonsdaleite phase within cubic dia-mond-lonsdaleite intergrowths. Possible attendance of nanocrystalline cubic diamond can be a reason of incorrect interpretation of Raman spectra in experimental study. It was found that the Raman modes dynamics under laser heating allows splitting of cubic diamond and lons-daleite modes and recognize diamond polymorphs [28]. The measurements in the region of 1300 - 1340 cm-1 show the presence of two Raman active modes, A:g (1319) and F1g (1322). The control measurement for co-sized cubic diamond particles give downshift of a Raman active mode F2g (1332) up to 1326 cm-1. Taking into account this downshift effect, we can declare the agreement of our prediction with the experimental data for lonsdaleite [28] at least for two high-frequency modes. We also obtain a rough agreement with recent measurements of Ref. [29], in which the most intensive band at 1292 - 1303 cm-1 and at 1219-1244 cm-1 are interpreted as contributions from A1g and E2g vibration modes of lonsdaleite phase in Pop-igai impact rock. Early DFT calculations, performed in the local density approximation (LDA), gave the following results: E:g(1312), A:g(1305), E2g(1193) [30] and E1g(1338), A1g (1280), E2g(1221) [31].
In Fig. 1a - e, we present our results for polycrystal-line (powder) Raman spectra with intensity plotted in arbitrary units. Of course, it has the dependence on direction for monocrystalline structures, but it is small and we discuss it later for real refractive indices only. The single peak of cubic diamond corresponding F2g(1332) mode is shown in all figures for comparison.
The number of active Raman modes or peaks depends on symmetry of crystal lattice and on the number of non-equivalent atoms in the asymmetric unit. As we see, quantum-mechanical approach predicts a very specific Raman spectrum for every allotrope, which can be considered as "finger prints" of these structures in the experimental search. The improvement of experimental methods to the level of accuracy of theoretical calculations will be crucial for search and recognition of diamond polymorphs.
The high-frequency refractive indices nii for different direction (ii = xx, yy, zz), accordingly Voigt notation, are collected in the Table 3. The calculation with PBE functional and POB-TZVP basis set gives answer which is very close to the experimental value for cubic diamond (~2.40) instead of the calculation with B3LYP functional used successfully for Raman spectrum calculation. The refractive asymmetry factor can be defined as follows Ap = 2(nn- nj / (nii+np). The values of Ap for lonsdaleite, C28, sic12 and 4H-diamond are about 2.5 - 3.5 %. This is enough large values for experimental test and recognition. The list of calculated IR active modes which initiate absorption of light in IR region are presented in Table 4. The ideal cubic diamond and lonsdaleite crystals do not absorb IR light. It means that experimental viewing of absorption for diamond-like structures may be used for recognition of new sp3 carbon allotropes.
Table 3. High-frequency refractive indices along the different directions
Basis set Structure POB-TZVP, PBE nxx, nyy, nzz Pople-6-21G, B3LYP nxx, nyy, nzz
Diamond 2.40, 2.40, 2.40 2.34, 2.34, 2.34
lonsdaleite 2.36, 2.36, 2.44 2.30, 2.30, 2.36
Mtn 2.11, 2.11, 2.11 2.10, 2.10, 2.10
SiC12 2.38, 2.38, 2.42 2.32, 2.32, 2.36
C28 2.35, 2.38, 2.43 2.29, 2.32, 2.36
4H-diamond 2.37, 2.37, 2.43 2.32, 2.32, 2.34
Table 4. IR active modes
^^^^^Basis set Structure Pople-6-21G, B3LYP
SiC12 Eu(425), Eu(578), A2u(728), A2u(1075), Eu(1271)A 2u(1330)
4H-diamond E1u(1248), A2u(1316)
C28 B2u(443), B2u(539), B3u(547), Bm(586), B1u(719), B3u(742), B2u(786), B2u(862), B3u(878), B2u(962), B3u(1000), B2u(1002), B1u(1071), B3u(1074), B3u(1172), B3u(1203), B2u(1211), B1u(1220), B2u(1239), B1u(1260), B2u(1278), B3u(1279), B3u(1298), B3u(1352), B2u(1356), B2u(1400)
mtn F1u(735), F1u(871), F1u(891), F1u(1059), F1u(1144), F1u(1223),
AE/atom, P, B,
Structure eV g/cm3 GPa
diamond 0.00 3.51 443
4H-diamond 0.01 3.61 445
SiC12 0.01 3.53 444
lonsdaleite 0.03 3.51 445
C28 0.07 3.35 427
mtn 0.11 3.06 383
AE/atom, eV p, g/cm3 B, GPa
diamond 0.00 3.43 428
4H-diamond 0.01 3.53 429
SiC12 0.01 3.43 428
lonsdaleite 0.03 3.43 429
C28 0.07 3.37 413
mtn 0.08 2.98 371
1200 1300
Wavenumber, cm'1
1200
Wavenumber, cm'1
1300
Wavenumber, cm'1
Fig. 1. Raman shift spectra in arbitrary units for carbon allotropes: a) lonsdaleite, b) mtn, c) SiC12, d) C28, e) 4H-diamond. The
peak at the 1332 cm1 in all panels corresponds cubic diamond
Optical and UV spectra
Optical properties are immediately connected with the electronic band structure of the crystal. To calculate the interband and intraband electron optical transitions one needs to solve the band structure of the system and know its band gap in different symmetry points of reciprocal space. Because of this reason, we calculated the electronic band structure for all the allotropes under study (see Fig. 5) and collected to the table our results for indirect and direct (in r-point) band gaps for all the considered structures within PBE and HSE06 functionals (see Table 5). It is well-known that the former fails to reproduce the optical band gaps while the latter provides the reliable results. Our result of 5.38 eV band gap for diamond successfully reproduces the well-known experimental value of 5.47 eV [1].
Then, we calculated the complex dielectric function for diamond, lonsdaleite and considered lowest-energy allotropes in the random phase approximation as it was explained in the Sec. "Computational methods and Details". In the Figs. 2 and 3 we present our predictions for diamond obtained with PBE exchange-correlation functional (bold solid line) and HSE06 functional (dashed line) together with experimental data (solid line) from the work [1] and PBE based predictions for lonsdaleite (dash-dotted line). We find a good agreement between theory and experiment within PBE functional calculations in the position and highness of the peaks, while with the HSE06 functional calculations, although providing correct optical band gaps, leads to the shift of the peaks approximately to 1 eV towards the high frequencies. Also, our results for cubic diamond and lonsdaleite are consistent with the earlier calculations in the work [32].
b) hv(eV)
Fig. 2. Real (a) and imaginary (b) dielectric functions: solid line -
experiment for cubic diamond, bold solid line - predictions for cubic diamond in PBE, dashed line - predictions for cubic diamond in HSE, dash-dotted line - lonsdaleite in PBE
The predictions for allotropes C28, mtn, sic are presented in the Fig. 3 in comparison with diamond.The most of new allotropes (SiC12, H-Carbon, C28) demon-
strate the very similar behavior of Si,2(v) to diamond or loncdaleite in the optical and UV region, so they do not lose the desirable optical parameters of diamond. In contrast, the mtn shows 2 times shorter and the smoothest peaks and has the smallest band gap in comparison to other allotropes, but the position of the peaks remains the same. Such a way, a presence of mtn in the polycrystal-line sample can contribute a visual opacity.
The effect of the anisotropy of the optical properties for non-cubic structures is found to be small, so we do not show the dependence of the dielectric tensor from the crystal directions in the figures, but only the average values as defined by formula (2.6).
Table 5. Band gaps for carbon allotropes
PBE, indirect HSE06, indirect PBE, direct HSE06, direct
diamond 4.67 5.38 5.64 7.04
lonsdaleite 3.34 4.91 4.96 6.37
mtn 3.76 5.09 3.76 5.03
Sic21 4.44 5.64 5.25 6.65
C28 4.77 5.96 4.77 6.06
4H-diamond 4.53 5.73 5.29 6.68
Conclusion
In summary, the first-principles quantum-mechanical calculations have been performed to obtain the frequency spectra and optical properties of lowest-energy sp3 carbon allotropes, such as: cubic diamond, lonsdaleite, SiC12, 4H-diamond, C28 and mtn. We have obtained Raman and IR spectra for all discussed allotropes and study their peculiarities. The electronic structure as well as the linear photon energy-dependent complex dielectric functions and related optical properties were computed.
Table 6. Lattice parameters for carbon allotropes calculated with POB-TZVP basis set and PBE exchange-correlation functional
Our investigations are beneficial to the experimental search and to the practical applications of these hypothetic carbon allotropes in IR, visible and UV photonics [33, 34].
Fig. 3. Real (a) and imaginary (b) dielectric functions in predicted for carbon allotropes in PBE: solid line - cubic diamond, bold solid line - C28, dashed line -mtn, dash-dotted line - SiC12
Structure N a, À b, À c, À Atomic positions
diamond 227 3.569 - - (0.125, 0.125, 0.125)
lonsdaleite 194 2.507 - 4.169 (0.333, 0.667, 0.062)
mtn 227 9.625 - - (0.880, 0.067, 0.683), (0.125, 0.125, 0.625) (0.033, 0.033, 0.716)
SiC12 166 2.511 - 24.822 (0.667, 0.333, 0.740), (0.667, 0.333, 0.572) (0.667, 0.333, 0.511), (0.000, 0.000, 0.656)
C28 55 7.298 7.728 2.535 (0.398, 0.039, 0.0), (0.376, 0.151, 0.5) (0.207, 0.270, 0.5), (0.239, 0.384, 0.0) (0.015, 0.196, 0.5), (0.460, 0.408, 0.0)
4H-diamond 195 2.534 - 8.355 (0.000, 0.000, 0.093), (0.667, 0.333, 0.155)
Table 7. Lattice parameters for carbon allotropes calculated with Pople's 6-21 basis set and B3LYP exchange-correlation
functional
Structure N a, À b, À c, À Atomic positions
diamond 227 3.594 - - (0.125, 0.125, 0.125)
lonsdaleite 194 2.527 - 4.205 (0.333, 0.667, 0.063)
mtn 227 9.688 - - (0.880, 0.068, 0.683), (0.125, 0.125, 0.625) (0.033, 0.033, 0.716)
SiC12 166 2.534 - 25.067 (0.333, 0.667, 0.073), (0.333, 0.667, 0.916) (0.333, 0.667, 0.844), (0.667, 0.333, 0.999)
C28 55 7.301 7.735 2.538 (0.398, 0.039, 0.0), (0.375, 0.151, 0.5) (0.268, 0.270, 0.5), (0.239, 0.384, 0.0) (0.015, 0.197, 0.5), (0.460, 0.408, 0.0)
4H-diamond 195 2.513 - 8.266 (0.000, 0.000, 0.093), (0.667, 0.333, 0.156)
Band structure
Density of states
Density of states
a)
d)
Band structure
Density of states
Band structure
Density of states
e)
Band structure
Density of states
Band structure
Density of states
c)
Fig. 4. Phonon band structures and density of states of allotropes: References
[1] Adachi S. Optical constants of crystalline and amorphous semiconductors. New York: Springer Science & Business Media; 1999. ISBN: 978-0-7923-8567-7.
[2] Yang N, ed. Novel aspects of diamond: From growth to applications. Cham, Heidelberg, New York, Dordrecht, London: Springer International Publishing Switzerland; 2015. ISBN: 978-3-319-09833-3.
[3] Zaitsev AM. Optical properties of diamond. Berlin, Heidelberg: Springer-Verlag; 2001. ISBN: 978-3-540-66582-3.
[4] Bundy FP, Kasper JS. Hexagonal diamond-A new form of carbon. J Chem Phys 1967; 46(9): 3437-3446. DOI: 10.1063/1.1841236.
[5] Frondel C, Marvin UB. Lonsdaleite, a hexagonal polymorph of diamond. Nature 1967; 214: 587-589. DOI: 10.1038/214587a0.
[6] Hirai H, Kenichi K. Modified phases of diamond formed under shock compression and rapid quenching. Science 1991; 253(5021): 772-774. DOI: 10.1126/science.253.5021.772.
[7] Hongliang H, Sekine T, Kobayashi T. Direct transformation of cubic diamond to hexagonal diamond. Appl Phys Lett 2002; 81: 610. DOI: 10.1063/1.1495078.
[8] Mao WL, Mao HK, Eng PJ, Trainor TP, Newville M, Kao CC, Heinz DL, Shu J, Meng Y, Hemley RJ. Bonding
ß
a) diamond, b) SiC12, c) mtn, d) lonsdaleite, e) C28, f) 4H-diamond
changes in compressed superhard graphite. Science 2003; 302(5644): 425-427. DOI: 10.1126/science.1089713.
[9] Hoffmann R, Kabanov AA, Golov AA, Proserpio DM. Homo Citans and carbon allotropes: for an ethics of citation. Angew Chem Int Ed 2016; 55(37): 10962-10976. DOI: 10.1002/anie.201600655.
[10] Hu M, Huang Q, Zhao Z, Xu B, Yu D, He J. Superhard and high-strength yne-diamond semimetals. Diamond and Related Materials 2014; 46: 15-20. DOI: 10.1016/j.dia-mond.2014.04.005.
[11] Wang JT, Chen C, Kawazoe Y. Mechanism for direct conversion of graphite to diamond. Phys Rev B 2011; 84(1): 012102. DOI: 10.1103/PhysRevB.84.012102.
[12] Baburin IA, Proserpio DM, Saleev VA, Shipilova AV. From zeolite nets to sp3 carbon allotropes: A topology-based multiscale theoretical study. Physical Chemistry Chemical Physics 2015; 17(2): 1332-1338. DOI: 10.1039/C4CP04569F.
[13] Nesper R, Vogel K, Blöchl PE. Hypothetical carbon modifications derived from zeolite frameworks. Angewandte Chemie International Edition 1993; 32(5): 701-703. DOI: 10.1002/anie. 1993 07011.
[14] Hohenberg P, Kohn W. Inhomogeneous electron gas. Phys Rev 1964; 136: B864-B871. DOI: 10.1103/PhysRev.136.B864.
arb. units
Fig. 5. Electronic band structures and density of states of allotropes: a) diamond, b)SiC12, c) mtn, d)lonsdaleite, e) C28, f) 4H-diamond
[15] Kohn W, Sham LJ. Self-consistent equations including exchange and correlation effects. Phys Rev 1965; 140: A1133-A1138. DOI: 10.1103/PhysRev.140.A1133.
[16] Dovesi R, et al. CRYSTAL14: A program for the ab initio investigation of crystalline solids. Int J Quantum Chem 2014: 114(19): 1287-1317. DOI: 10.1002/qua.24658.
[17] Kresse G, Furthmuller J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys Rev B 1996; 54(16): 11169-11186. DOI: 10.1103/PhysRevB.54.11169.
[18] Pascale F, Zicovich-Wilson CM, Lopez F, Civalleri B, Orlando R, Dovesi R. The calculation of the vibration frequencies of crystalline compounds and its implementation in the CRYSTAL code. J Comput Chem 2004; 25(6): 888897. DOI: 10.1002/jcc.20019.
[19] Maschio L, Kirtman B, Rerat M, Orlando R, Dovesi R. Ab initio analytical Raman intensities for periodic systems through a coupled perturbed Hartree-Fock/Kohn-Sham method in an atomic orbital basis. I. Theory. J Chem Phys 2013; 139(16): 164101. DOI: 10.1063/1.4824442.
[20] Ferrero M, Rerat M, Orlando R, Dovesi R. The calculation of static polarizabilities in 1-3D periodic compounds. The implementation in the CRYSTAL code. J Comput Chem 2008; 29(9): 1450-1459. DOI: 10.1002/jcc.20905.
[21] Perdew JP, Burke K, Ernzerhof M. Generalized gradient approximation made simple. Phys Rev Lett 1996; 77(18): 3865-3868. DOI: 10.1103/PhysRevLett.77.3865.
[22] Peintinger MF, Oliveira DV, Bredow T. Consistent gaussian basis sets of triple-zeta valence with polarization quality for solid-state calculations. J Comp Chem 2013; 34(6): 451-459. DOI: 10.1002/jcc.23153.
[23] Becke AD. Density-functional thermochemistry. III. The role of exact exchange. J Chem Phys 1993; 98(7): 56485652. DOI: 10.1063/1.464913.
[24] Gordon MS, Binkley JS, Pople JA, Pietro WJ, Hehre WJ. Self-consistent molecular orbital methods. 22. Small split-valence basis sets for second-row elements. J Am Chem Soc 1982; 104(10): 2797-2803. DOI: DOI: 10.1021/ja00374a017.
[25] Baima J, Zelferino A, Olivero P, Erba A, Dovesi R. Raman spectroscopic features of the neutral vacancy in diamond from ab ignition quantum-mechanical calculations. Phys Chem Chem Phys 2016; 18(3): 1961-1968. DOI: 10.1039/C5CP06672G.
[26] Heyd J, Scuseria GE, Ernzerhof M. Erratum: Hybrid func-tionals on a screened Coulomb potential. J Chem Phys 2006; 124(21): 219906. DOI: 10.1063/1.2204597.
[27] Isaenko S, Shumilova T. Thermostimulated Raman spectrum dynamics of lonsdaleite. EGU General Assembly 2012; 14: 608.
[28] Goryainov SV, Likhacheva AY, Rashchenko SV, Shubin AS, Afanas'ev VP, Pokhilenko NP. Raman identification of lonsdaleite in Popigai impactites. J Raman Spectrosc 2014; 45(4): 305-313. DOI: 10.1002/jrs.4457.
[29] Filik J, Harvey JN, Allan NL, May PW, Dahl JEP, Liu S, Carlson RMK. Raman spectroscopy of nanocrystalline diamond: An ab initio approach. Phys Rev B 2006; 74(3): 035423. DOI: 10.1103/PhysRevB.74.035423.
[30] Wu BR, Xu J. Total energy calculations of the lattice properties of cubic and hexagonal diamond. Phys Rev B 1998; 57(21): 13355-13359. DOI: 10.1103/PhysRevB.57.13355.
[31] Denisov VN, Mavrin BN, et. al. First-principles, UV Raman, X-ray diffraction and TEM study of the structure and lattice dynamics of the diamond-lonsdaleite system. Diamond & Relared Matherials 2011; 20(7): 951-953. DOI: 10.1016/j.diamond.2011.05.013.
[32] Wang Z, Zhang RJ, Zheng YX, et al. Electronic and optical properties of novel carbon allotropes. Carbon 2016; 101: 77-85. DOI: 10.1016/j.carbon.2016.01.078.
[33] Kaminskii AA, Ral'chenko VG, Yoneda H, Bol'shakov AP, Inyushkin AV. Stimulated Raman scattering-active isotopi-cally pure 12C and 13C diamond crystals: A milestone in the development of diamond photonics. JETP Letters 2016; 104(5): 347-352. DOI: 10.1134/S0021364016170082.
[34] Salvatori S, Girolami M, Oliva P, Conte G, Bolshakov A, Ralchenko V, Konov V. Diamond device architectures for UV laser monitoring. Laser Phys 2016; 26: 084005. DOI: 10.1088/1054-660X/26/8/084005.
Authors' information
Vladimir Anatolievich Saleev (b. 1959) graduated from Kuibyshev State University (presently, Samara National Research University) in 1981, majoring in Theoretical Physics. Currently he works as Professor at the Samara National Research University. Research interests are quantum mechanics, quantum field theory, theoretical materials science and high energy physics. E-mail: saleev@samsu.ru .
Alexandra Victorovna Shipilova (b. 1987) graduated from Samara State University (presently, Samara National Research University) in 2008, Department of Physics. Currently she works as Associated Professor at the Samara National Research University. Research interests are quantum mechanics, quantum field theory, theoretical materials science and high energy physics. E-mail: alexshipilova@samsu.ru .
Code of State Categories Scientific and Technical Information (in Russian - GRNTI)): 29.31.21, 29.31.27, 29.19.24.
Received May 10, 2017. The final version - June 23, 2017.
