CC BY
37
Supercritical fluid synthesis and possible properties of "cubic graphite"
A. V. Pokropivny1'2, A.N. Enyashin3, A. S. Smolyar1'4, V. A. Kuts5, V. G. Gurin5, S. A. Antipov5, P. M. Silenko1, Yu. M. Solonin1
1Frantsevich Institute for Problems of Materials Science NASU, Kyiv, Ukraine 2G. V. Kurdyumov Institute for Metal Physics, NASU, Kyiv, Ukraine 3Institute of Solid State Chemistry UB RAS, Ekaterinburg, Russia 4M.P. Semenenko Institute of Geochemistry, Mineralogy and Ore Formation NASU, Kyiv, Ukraine 5 Institute of Magnetism NASU, Kyiv, Ukraine [email protected], [email protected]
PACS 81.10.Dn, 61.66.Bi, 61.48.+c, 61.50.Ah DOI 10.17586/2220-8054-2020-11-4-408-416
We report on supercritical fluid synthesis of an intermediate carbon phase - austite - at a pressure of 180 MPa and temperatures 500-700° C from soot as a precursor and supercritical carbon dioxide as a solvent. According to the results of electron and X-ray diffractions, spectral measurements and density-functional theory calculations, the observed carbon phase is proved to be cubic with a lattice parameter value of 8.96±0.05 A and a possible structural type as for KFI zeolite.
Keywords: crystal growth, carbon allotropes, electron microscopy, X-ray diffraction, DFT calculations. Received: 15 June 2020 Revised: 10 July 2020
1. Introduction
Supercritical fluid synthesis is a common method for the high-pressure high-temperature synthesis of different oxides, sulfides, nitrides, carbon phases as well as self-assembling of nanomaterials [1,2]. Different gases can be employed for creating pressure during supercritical fluid synthesis as an inert medium (He, Ar) or a solvent and a chemical reagent (H2, N2, H2O, CO2 etc.). The dissolving capacity of fluids is rapidly increased with increasing of the pressure, which is also important for supercritical fluid extraction. On the other hand, the velocity of crystallization in supercritical medium is increased up to 0.1-0.3 mm/h, which is at least five times greater than for hydrothermal methods. Moreover, crystals formed from the gas phases may have a high degree of perfection, in comparison with other methods.
In 1963 Aust and Drickamer claimed on synthesis of a new cubic carbon phase (henceforth, austite) by transformation from single crystal of graphite at 15 GPa and 180 K [3]. There was a sharp rise of resistance with pressure at this point, accompanied by drifting upward with time - a behavior typical of a sluggish first-order phase transition. The mixture of new phase with graphite was found to have the density of 2.35-2.4 g/cm3 with the float-sink tests. X-ray diffraction analysis of the partially transformed mixture showed that this phase could be cubic with the lattice parameter of 5.545 A, having 24 atoms per unit cell and a crystal density of ~ 2.8 g/cm3 [3].
In 1967 Bandy reported, that Drickamer rejected on the facts of the synthesis of austite in a private communication [4]. Therefore, the discovery of a new phase in carbon phase diagram was not accepted by the scientific community. It should be noted that authors [3] had not published official refutation of their work. Presumably the indisputable authority of Bandy misled other scientists, who could confirm or refute completely their extraordinary results. As a result, in 1974 the pattern of this phase (No 18-311) was deleted from the database. Nevertheless, at the moment several forms of diamond, nanotubes, fullerenes, a large number of carbon fullerites, and other allotropic modifications have been discovered on carbon phase diagram, hence, an observation of new carbon phase is not a surprise since the last half of twentieth century [5].
However, the austite-like phases have been repeatedly observed after Aust and Drickamer. Namely, Fedoseev et. al. reproduced an austite-like phase during the growth of diamond from the gas phase [6], by Smolyar et. al. with a supercritical fluid treatment of graphite [7], and, possibly, by Shterenberg et. al. from channel black with a heat treatment at 1250°C in the presence of Ni at 3.7 MPa [8]. The hardness of austite was found to be intermediate, namely 1-5 GPa [7]. Shumilova et. al. discovered austite-like phase in the nature among the products of graphite mineralization [9]. The mineral lignite from Dakota Star Mines [10], the existence of which was questionable, could also be hypothesized as this phase. The structure was interpreted for the first time as a cubic lattice of C24 fullerenes arranged as zeolite-like LTA lattice [11,12].
In the last two decades, other new phases of carbon and boron nitride attract outstanding interest of the same or even higher level, as fullerenes and nanotubes do. All of these new materials can be used in the architectures of future devices or, alternatively, as advanced materials with prominent absorption properties or superhardness. Among these new phases as the C8 phase [13], E-phase of boron nitride [14,15] and fulborenite B12N12 [16] are remarkable. Crystal structures of these phases were proposed quite recently, although they were synthesized in the last half of the twentieth century. The problem of structure determination and attribution of new forms of carbon and BN is ordinary due to their small mass yields after synthesis, unlike that, for example, nanotubes and fullerenes, and due to low atom weight of constituent elements. Theoretically, several hypothetical zeolite-like carbon structures (ATO, KAN, ATN, AFI etc.) are proposed for future identifications [17-26].
The goal of this work is to obtain a sample of austite phase ("cubic graphite") using the supercritical fluid synthesis from soot and to revise the structural data on the received samples. The density-functional theory methods are employed to supply the results of experiments.
2. Experimental part
2.1. Synthesis and characterization
Equipment for supercritical fluid synthesis represents the reactor with an inner holder for the samples, aggregates for compressions and cleaning of gas, heating system, and automatic control system. The equipment allows one to create the gas pressures up to 1 GPa for a short time period (0.5-1 h) and to maintain the pressure during the long time period (hours, days) [1-7].
The soot was used as a precursor and supercritical carbon dioxide as a solvent. It was saturated in the reactor at pressure 95 MPa for 25 hours at ambient temperature. After that, the temperature was cyclically increased and reduced in the intervals of 500-700° C during 4 hours at pressure 180 MPa. The reaction product was observed as polycrystalline plates and as separate crystals in the soot.
Transformed soot powders were studied using different physical methods. X-ray powder diffraction (XRPD) patterns were received on DRON-2 diffractometer in the range of 26 = 5°-90° at CuKa radiation. Electron diffraction (ED) patterns were obtained at room temperature with a Philips EM-400T transmission electron microscope, operating at acceleration voltage of 100 kV. The samples were prepared by crushing in a mortar and dispersal of the synthesized crystallites in ethanol, and subsequent deposition on a grid containing holes. The elemental composition of the plates of new phase was studied with the JEOL spectrometer.
2.2. Computational details
The calculations were performed within the framework of the density-functional theory (DFT) using the SIESTA 4.0 implementation [27, 28]. The exchange-correlation potential within the Generalized Gradient Approximation (GGA) with the Perdew-Burke-Ernzerhof parametrization was used. The core electrons were treated within the frozen core approximation, applying norm-conserving Troullier-Martins pseudopotentials. The valence electrons for C were taken as 2s22p2. The pseudopotential core radii were chosen, as suggested by Martins, and equal to 1.3 aB for sand p-states of C. In all calculations double-Z basis set was used. The k-point mesh was generated by the method of Monkhorst and Pack with a cutoff 10 A for the k-point sampling. The real-space grid used for the numeric integrations was set to correspond to the energy cutoff of 200 Ry. All calculations were performed using variable-cell and atomic position relaxations, with convergence criteria corresponding to the maximum residual stress of 0.1 GPa for each component of the stress tensor, and the maximum residual force component of 0.01 eV/A.
Apart from the geometry optimization using DFT method, the stability, mechanical and thermal properties were tested using geometry optimization and molecular-dynamics (MD) simulations within the framework of the density-functional tight-binding method (DFTB) [29] in r-point approximation as implemented in the deMon software [30]. The structures were MD annealed during 20 ps with the time step of 1.0 fs as canonical (NVT) ensembles at temperatures T = 300 and 600 K using global Berendsen thermostat with the time constant 100 fs.
The optimized geometry of all carbon allotropes was used to simulate theoretical XRD spectra for wavelength A = 1.5406 A (CuKa radiation). All structures were considered as periodic 3D crystals. The smearing of reflection profiles was approximated with correction for the isotropic atomic temperature factor [31].
3. Results and Discussion
3.1. Electron diffraction analysis
ED patterns from both untransformed soot and synthesized crystal are collected in Fig. 1. It should be noted that synthesized particles were found to be hard, in comparison with starting material, in good agreement with previous research of this phase [7]. The test samples consisted of crystalline microparticles of 10 to 100 ^m in size. The
observed crystals were solid, unlike carbon black, but crumble well in an agate mortar. Therefore, their hardness may be estimated as the one not exceeding 6.5-7.0 units on the Mohs scale. Microparticles of untransformed soot were amorphous or badly-crystalline with the characteristic parameter of 3.6 A, hence, larger than that for graphite. The patterns of the synthesized sample can be attributed without any restriction to a cubic phase with the characteristic lattice parameter 8.96±0.05 A.
(a) [001] * (b) tool] # 1 nm
(c) [012] -•* - i (d)
Fig. 1. Electron diffractograms in different projections (a-c) for the carbon samples prepared using supercritical fluid synthesis and in comparison to the diffractogram for the carbon precursor (d). Common scale bar is given in (b)
Qualitative analysis of the elemental composition using TEM spectrometer reveals only carbon, and, possibly, light elements in the samples under investigations. Using spectral quantitative analysis at the surface of the polycrystal plates, the composition was found to be 98-99 at% of carbon, with the 0.5-1 at% of oxygen and much less then 1 % of other elements. Thus, the phase is really a pure carbon, while oxygen is likely adsorbed on the surface, and the rest elements are impurities arising from the reactor and from the sample preparation for the spectrometer.
3.2. X-ray diffraction analysis
Results of XRPD analysis of the sample of powders are listed in column 1 of Table 1. For comparison, the theoretical XRD reflections of austite have been calculated with the JSV computer code for the structure proposed as carbon KFI zeolite [32] with the space group Im-3m (column 2 of Table 1). The hkl indices, dhkl-spacings (prior to brackets), relative intensities (in brackets) and recalculated lattice parameters (after brackets) are tabulated. Interpretation of reciprocal spaces clearly indicates the presence of a single cubic phase of 8.95±0.15 A, in excellent agreement with the results of electron diffraction. Simultaneously, the lattice parameter of the austite sample from [3] (column 3 of Table 1) can be found slightly greater, which may be related to the difference in sample preparation. Rietveld refinements are not possible because the powders of new phase are pyrogenetically bound to graphite or some other form of carbon. This is not the only problem associated with carbon or boron nitride phases.
Column 2 of Table 1 shows the results of model powder diffraction for zeolite-like KFI carbon with the lattice constant a = 8.95 A. The calculated diffraction patterns from the proposed KFI structure correlate well with those from our sample and from report of Aust et al. [3]. There are (200), (222) and (330) main calculated reflexes for the proposed structure. Experimentally, two of them are presented for our sample in good agreement with the calculations. Though, main experimental reflexes (110), (200), (220), (332) have only weak analogs of theoretical reflexes as calculated for zeolite-like KFI carbon. The synthesized phase should have a texture along (220) and (110) planes, since their two reflexes rank first in intensity. The (222) strong reflex is not observed by Aust et al. [3], but it is observed from our samples. And vice versa (422) reflex is observed by Aust et al. [3], but it is not observed from our sample.
Table 1. Theoretically derived as for carbon KFI zeolite (column 2), and experimental XRPD patterns (column 3 and 4) with corresponding lattice parameters (in A after brackets). Subscripts in column 1 denote the formerly ascribed indices for this phase [3]. X-ray diffraction pattern of the experimental samples contains additionally several peaks of different intensities, which correspond to other phases
hkl Theoretical KFI structure for a = 8.95 A Experimental in our work Experimental [3]
100 9.10 (5), a = 9.10
110 6.329 (6) 6.27 (30), a = 8.87
111 5.30 (5), a = 9.17
200 4.475 (100) 4.42 (15), a = 8.84
210(110) 3.95 (10), a = 8.83
211 3.654 (< 1) 3.61 (10), a = 8.84
220(111) 3.164(3) 3.17(100), a = 8.96 3.208 (10), a = 9.07
310(200) 2.830 (10) 2.80 (10), a = 8.85 2.770 (50), a = 8.76
222 2.584 (41) 2.55 (10), a = 8.83 2.467 (50), a = 8.90
321 2.392 (< 1) 2.365 (10), a = 8.85
400(211) 2.238 (3) 2.214(10), a = 8.86
330,411 2.110(25)
420 2.001 (5) 2.00 (20), a = 8.94
332 1.908 (5) 1.932 (30), a = 9.06 1.961 (10), a = 8.99
421(220) 1.932 (30), a = 8.85
422(221) 1.827 (4) 1.844(10), a = 9.03
510, 431 1.755 (< 1)
521 1.634 (1)
440(222) 1.582(1) 1.589 (5), a = 8.99 1.600(10), a = 9.05
433, 530 1.535 (3)
442, 600(321) 1.492 (< 1) 1.485 (10), a = 8.91
532 1.452 (6)
It should be noted here that a set of structural transformations influence on the powder diffraction patterns. The essential differences in peaks distribution exist not only for different method of genesis, but also from probe to probe of the same method. In particular, sample of [3] is more textured with a poor set of diffraction lines in comparison with our sample. Nevertheless, first three XRD peaks can be fairly attributed to the peaks of our sample.
3.3. DFT computational data
The preliminary experimental results for observed cubic carbon phase with the lattice parameter ~ 8.95 A suggest the model corresponding to zeolite-like KFI lattice, which impose the unit cell consisting of 96 C atoms - C96 (vide infra). Here, we employ several DFT approaches to infix the possible crystal structure of synthesized carbon phase as austite, confirming its dynamical and mechanical stabilities. Apart from the zeolite-like KFI C96, the diamond with face-centered cubic lattice was chosen as a reference system to study structural, electronic and elastic properties. In addition, similar calculations were also performed for other porous carbon - C24 fullerite with simple cubic lattice (or LTA carbon zeolite), where C24 fullerenes are connected by square-like edges [11,12]. While the properties of diamond are profoundly established from both experimental and theoretical views, the C24 fullerite was selected due to its structural relativity to C96 (Fig. 2). The convenient KFI structure of C96 can be represented as a primitive body-centered cubic lattice of C48 building blocks. The latter have the shape of truncated cuboctahedra and are assembled by means of hexagonal edges creating a characteristic motif of hexagonal prisms. Yet, simple cubic lattice of C24 fullerite may be also represented as composed of the same C48 building blocks linked via square-like edges, which form a primitive face-centered cubic lattice.
The calculated structural, electronic and mechanical properties of diamond have validated the calculational approaches used in this work and demonstrated their fair adequacy for the reproduction of basic properties of carbon
3
Fig. 2. Structures of subtle carbon allotropes studied in this work: the fragments of austite C96 with zeolite-like KFI lattice (bcc-C48) (1) and C24 with zeolite-like LTA lattice (fcc-C48) (2) after GGA-DFT geometry optimization; the fragments of austite C96 after DFTB molecular-dynamics simulations at 600 K (3)
allotropes (Table 2). The error in the estimation of lattice parameter does not exceed 0.3 %. In agreement with the experimental data, the cubic diamond is characterized as a wide band gap semiconductor with non-direct type of transition from r-point to (2/3) T-X point (Fig. 3). As in many DFT approaches, the value of the band gap obtained after DFT-GGA calculations is underestimated by ~ 1 eV. In turn, DFTB approach typically overestimates the band gap, while the features of the band structure are as well fairly reproduced. The values of the bulk modulus estimated either from total energy curve or calculated from independent elastic constants may be found in the fair agreement with the experimental value (Tables 2 and 3). In general, apart of the band gap estimations, both employed calculational methods yield equivalent results.
The calculations of C96 allotrope were performed using primitive body-centered cubic unit cell, which contains 48 atoms. The lattice of subtle C96 allotrope preserves bcc-structure after the geometry optimization (Fig. 2). Contrary to the diamond, the band structure calculation at DFT-GGA level reveals that C96 is a semiconductor with the smaller and direct band gap of about 2.6 eV (Table 2). However, the partial densities of states at the valence band of both carbon allotropes have dominant 2pC character.
A highly strained coordination of sp3-hybridized carbon atoms in C96 cannot be energy gainful and the total energy calculations evidence the lower stability of C96 carbon comparing to the diamond almost on ~ 1 eV per C-atom. The first hint at the dynamical stability of this lattice was given using molecular dynamics simulations of the supercell consisting of 16 C48 building blocks. A reasonably high stability of C96 allotrope can be revealed in the absence of external pressure and at temperature of 300 K, when no change in the structure can be obtained and the
Table 2. Main characteristics of C96 carbon allotrope (KFI-like lattice), C24 fullerite (LTA-like lattice) and diamond in conventional cubic lattice representation: number of atoms in the unit cell Z, relative formation energy (AE, in eV/atom), lattice parameter (a, in A), mass density (p, in g/cm3), bulk modulus (B, GPa) and band gap (Eg, in eV), as calculated at different DFT levels
System Method Z AE a P B Eg
fcc-diamond DFT GGA 8 0.0000 3.58 3.47 493 4.20
DFTB 0.0000 3.58 3.47 487 7.73
exp.[25] 3.57 3.52 442 5.46-5.6
sc-C96 or bcc-C48 DFT GGA 96 0.8382 9.29 2.39 285 2.61
DFTB 1.0659 9.35 2.34 276 7.15
sc-C24 or fcc-C48 DFT GGA 24 0.7501 5.93 2.30 290 2.42
DFTB 0.8892 5.95 2.27 280 7.15
Fig. 3. Band structure and densities of states (DOS) for cubic diamond (1), for austite C96 with zeolite-like KFI lattice (or bcc-C48 ) (2) and for C24 fullerite with zeolite-like LTA lattice (or fcc-C48 ) (3) calculated within DFT-GGA approach. 2pC and 2sC partial densities of states are painted in light gray and gray, respectively
atom motions are not accompanied with the break of chemical bonds. However, the structure of C96 undergoes a considerable destruction with the partial destruction of hexagonal prisms and occurrence of double C=C bonds at T = 600 K (Fig. 2).
Thereafter, the stability of the lattice was estimated in more details by means of the widely used condition of intrinsic stability, i.e. when the Gibbs free energy of the crystal is in a local minimum with respect to small structural deformations. In terms of elastic constants (Cj), this stability criterion (so-called mechanical stability) for cubic crystals requires that (Cn - C12) > 0, (Cn + 2Ci2 ) > 0, and C44 > 0. In the framework of less sophisticated DFTB method the values of three independent elastic constants (C11, C12 and C44) were evaluated and compared with those of diamond (Table 3). These constants satisfy the aforementioned conditions. Thus, we may assert that the proposed C96 allotrope should be intrinsically stable.
Calculated elastic constants allow to estimate several other important mechanical properties, which are listed along with available experimental values in Table 3: bulk modulus B = (Cn + 2C12)/3, compressibility fi = 1/B, shearmodulus G = (3C44+Cn —Ci2)/5, Pugh's indicator k = G/B, tetragonal shearmoduli G' = 1/2(Cn —Ci2), Young's modulus Y = 9BG/(3B + G), Cauchy pressure CP = C12 — C44, Poisson's ratio v = 1/2(1 — Y/3B), Zener anisotropy index, AZ = 2C44/(C11 — C12).
The Young's modulus Y of a material is defined as a ratio of linear stress to linear strain, which determines its stiffness. In our study, the Young's modulus of C96 lattice is found to be Y = 475 GPa, thus, the material will show a rather high stiffness.
Brittle/ductile behavior is one of many important mechanical characteristics of a material, which is closely related to its reversible compressive deformation and fracture ability. Often, the ductile/brittle nature of a material is related to its Cauchy pressure CP: if CP is negative, the material is expected to be brittle. Another most widely used materials' malleability measure is the Pugh's index k [34]. Namely, at k < 0.57, a material behaves in a ductile manner, and vice versa, at k > 0.57, it demonstrates brittleness. According to both indicators (Table 3), C96 carbon allotrope should behave as a brittle material.
Elastic anisotropy of crystals reflects a different bonding character in different directions and has an important implication since it correlates with the possibility to induce microcracks in materials. The elastic anisotropy can be estimated using the Zener's anisotropy index AZ [35]. For isotropic crystals, AZ = 1, and the deviations of AZ from unity define the extent of elastic anisotropy. From this point of view, C96 allotrope with AZ = 0.65 should be elastically quite isotropic.
The calculated elastic parameters allow estimate Vickers hardness of the material (By) in a few simplified ways. Generally, the hardness of a material is characterized experimentally by indentation and depends strongly on plastic (irreversible) deformation. However, a set of empirical relationships between Vickers hardness and elastic properties was proposed. Herein, the following semi-empirical correlations between Vickers hardness and shear, Young's modulus, and Pugh's index were applied [36,37]: By (1) = 0.1769G -2.899; By (2) = 0.0608Y; By (3) = 2(k2G)0'585 - 3. Present numerical estimations of Hy (Table 3) demonstrate that all three correlations give comparable values and suggest that the C96 allotrope should be a hard material.
As was previously mentioned, C96 or bcc-C48 lattice is related to the lattice of C24 fullerite or fcc-C48. Thus, the simultaneous fabrication of both these carbon allotropes might be expected. In addition, the stability, electronic and mechanical properties were estimated also for C24 fullerite at the same levels of theory (Tables 2 and 3, Fig. 3). Noteworthy, the latter phase characterized earlier by Pokropivny et al [11,12] is slightly more stable, than C96, and posses very similar trends in the mechanical stability of lattice and in the electronic properties. Obviously, the coproduction or the production of individual C24 fullerite under the same experimental conditions of supercritical fluid synthesis cannot be fully excluded.
Table 3. Mechanical properties of sp3 carbon allotropes as obtained after DFTB calculations: independent elastic constants (Cj, in GPa), bulk modulus (B, in GPa), compressibility (fi, in GPa-1), shear modulus (G, in GPa), Pugh's indicator (k), tetragonal shear moduli (G', in GPa), Young's modulus (Y, in GPa), Cauchy's pressure (CP, in GPa), Poisson's ratio (v), Zener's anisotropy index (AZ), and Vickers hardness (Bv, in GPa)
System C11 C12 G44 B £
fcc-diamond 1103.5(1079*) 178.8 (124*) 602.2 (578*) 487 (442*) 0.0021
sc-C9e (bcc-C48) 607.0 109.7 161.7 275 0.0036
SC-C24 (fcc-C48) 550.7 145.4 259.7 281 0.0036
System G k G' Y CP
fcc-diamond 546 1.12 462 (478*) 1192(1050*) -423.4
sc-C96 (bcc-C48) 196 0.71 249 475 -52.0
sc-C24 (fcc-C48) 237 0.84 203 555 -114.3
System V Az Hv (1) Hv (2) Hv (3)
fcc-diamond 0.09 (0.1*) 1.30 93.7 72.5 88.2
sc-C96 (bcc-C48) 0.21 0.65 31.8 28.9 30.8
sc-C24 (fcc-C48) 0.17 1.28 39.0 33.7 37.0
* experimental values from Ref. [33]
4. Summary
In 1963 Aust and Drickamer [3] found a new phase during transformation of graphite at high pressures. In 2004 a few members of our group offered the structure of this phase on the basis of available data. At that time, the phase was indexed in the cubic system with a period of ~ 5.55 A with the structure of LTA zeolite [11,12]. With the support from electron diffraction technique the lattice parameter has been refined and its revised value is found to be approximately ^9 A instead of ~ 5.55 A, according to previous experimental data [3] and theoretical [11,12] reports. The synthesized austite is suggested as carbon phase Cg6 with zeolite-like KFI lattice (space group Im - 3m) in accordance with analysis of electron and X-ray diffraction patterns. As independent method, DFT calculations confirm the mass density, dynamical and mechanical stability of the proposed structure.
In the present experimental work the high level of nucleation of new crystalline carbon phase from disordered material is for the first time observed without additional catalysts, with only supercritical carbon dioxide as a solvent. The presence of graphite and, possibly, chaoite and C8 sodalite phases [13] as concomitant products of synthesis indirectly confirms an assumption that new carbon phases are relatively stable and not rare. Indeed, more then 15 allotropes of carbon have been discovered in nature [38]. Apart from the Cg6 carbon allotrope with KFI-like structure a wealth of other carbon allotropes might be potentially fabricated under extreme conditions of fluid synthesis. Being fabricated under the conditions far from equilibrium, not all these phases can be hosted at the classical phase diagram of carbon, what can entangle their confident identification and characterization.
While the basic structural, electronic and mechanical properties have been established theoretically, the experimental characterization of newly fabricated carbon C96 phase still remains in its infancy. A set of new synthesis experiments, in particular, under different pressure-temperature and kinetic parameters are planned to increase the yield of cubic carbon and to reveal in details the formation mechanism. The latter would facilitate the understanding of phase transformations between numerous known and hypothetical carbon allotropes from experimental point-of-view, and further stimulate the search of new applications of these phases, e.g., in electronics.
Acknowledgments
We acknowledge M. Czank and W. Depmeier for their assistance with Electron diffraction method according to DFG visiting grant. A.E. acknowledges a support from ISSC UB RAS (theme AAAAA19-119031890025-9). Visualizations of atomic structures and MD trajectories were performed using VMD software (www.ks.uiuc.edu/Research/vmd/) [39].
References
[1] Holmes J.D., Lyons D.M., Ziegler K.J. Supercritical fluid synthesis of metal and semiconductor nanomaterials. Chemistry Eur. J., 2003, 23, P. 2144-2150.
[2] Smolyar A.S., Barcholenko V.A., Kuts V.A., Gurin V.G., Archipov A.P., Maloshtan S.A., Razvadovsky N.A., Svyato V.P., Klimovich A.P. UA patent No 98084606. Bulletin No 6, 2002.
[3] Aust R.B., Drickamer H.C. Carbon: a new crystalline phase. Science, 1963, 140, P. 817; Deleted pattern No 18-311.
[4] Bandy F.P., Kasper J.S. Hexagonal diamond - a new form of carbon. J. Chem. Phys., 1967, 46, P. 3437-3446.
[5] Pokropivny V.V., Ivanovsky A.L. New nanoforms of carbon and boron nitride. Russ. Chem. Rev., 2008, 77, P. 837-874.
[6] Fedoseev D.V., Deryagin B.V., Varnin V.P., VnukovS.P., TeremetskayaI.G., PolyanskayaN.D. Polymorphism in carbon-boron nitride systems. Dokl. Akad. NaukSSSR, 1976, 228, P. 371-374.
[7] Smolyar A.S., Sozin Yu.I., Barholenko V.A., Maloshtan S.N., Kuts V.A., Gurin V.G., Arhipov A.P., Gerasimov A.Yu., Razvadovskii N.A., Titenko A.N.Fluid synthesis of carbon phases. J. SuperhardMater., 2002, 2, P. 79.
[8] Shterenberg I., Bogdanova S. Influence of nickel on graphitization of carbon materials at high-pressures and temperatures. Inorg. Mater., 1979, 15, P. 632-636.
[9] Shumilova T.G., Kablis G.N., Pushkarev E.V. Typomorphic features of graphite mineralization of probable alternative high-pressure sources of diamond: cubic graphite. Dokl. Earth Sciences, 2002, 387, P. 958-961.
[10] Wild R., Schellenbaum R. X-ray diffraction studies of mineral matter in North Dakota Lignite. Proc. North Dakota Acad. Sci., 1951, 5, P. 40-42.
[11] Pokropivny V.V., Pokropivny A.V. Structure of "cubic graphite": simple cubic fullerite C24. Phys. Solid State, 2004, 46, P. 392-394.
[12] Pokropivny A.V. X-ray diffraction analysis of austite (cubic graphite), a novel carbon [C24]-LTA zeolite. Physics Low.-Dim. Struct., 2006, 2, P. 64-68.
[13] Pokropivny A.V., Volz S. C8 phase: supercubane, tetrahedral, BC-8 or carbon sodalite? Phys. Status Sol. B, 2012, 249, P. 1704-1708.
[14] Pokropivny A.V. Structure of the boron nitride E-phase: diamond lattice of Bi2Ni2 fullerenes. Diam. Relat. Mat., 2006, 15, P. 1492-1495.
[15] Pokropivny A.V., Volz S. Hybrid porous nanotube crystal networks for nanostructured device applications. J. Mater. Sci., 2013, 48, P. 29532960.
[16] Pokropivny V.V., Smolyar A.S., Pokropivny A.V. Fluid synthesis and structure of a new boron nitride polymorph-hyperdiamond fulborenite B12N12 (E phase). Phys. Solid State, 2007, 49, P. 591-598.
[17] Baburin I.A., Proserpio D.M., Saleev V.A., Shipilova A.V. From zeolite nets to sp3 carbon allotropes: a topology-based multiscale theoretical study. Phys. Chem. Chem. Phys., 2015, 17, P. 1332-1338.
[18] Hu M., Zhao Z., Tian F., Oganov A.R., Wang Q., Xiong M., Fan C., Wen B., He J., Yu D., Wang H.-T., Xu B., Tian Y. Compressed carbon nanotubes: a family of new multifunctional carbon allotropes. Sci. Rep., 2013, 3, 1331.
[19] Belenkov E.A., Brzhezinskaya M.M., Greshnyakov V.A. Novel carbon diamond-like phases LA5, LA7 and LA8. Diam. Relat. Mat., 2014, 50, P. 9-14.
[20] Belenkov E.A., Greshnyakov V.A. Diamond-like phases formed from fullerene-like clusters. Phys. Solid State, 2015, 57, P. 2331-2341.
[21] Belenkov E.A., Brzhezinskaya M.M., Greshnyakov V.A. Crystalline structure and properties of diamond-like materials. Nanosystems: Physics Chemistry Mathematics, 2017, 8, P. 127-136.
[22] Belenkov E.A., Greshnyakov V.A. Modeling of phase transitions of graphites to diamond-like phases. Phys. Solid State, 2018, 60, P. 12941302.
[23] Greshnyakov V.A., Belenkov E.A. Structures of diamond-like phases. J. Exp. Theor. Physics, 2011, 113, P. 86-95.
[24] Kvashnina Y.A., Kvashnin D.G., Kvashin A.G., Sorokin P.B. New allotropic forms of carbon based on D-60 and D-20 fullerenes with specific mechanical characteristics. JETP Lett., 2017,105, P. 419-425.
[25] Kvashnina Y.A., Kvashnin A.G., Popov M.Y., Kulnitskiy A.A., Perezhogin I.A., Tyukalova E.V., Chernozatonskii L.A., Sorokin P.B., Blank V.D. Toward ultra-incompressible carbon materials: computational simulation and experimental observation. J. Phys. Chem. Letters, 2015, 6, P. 2147-2152.
[26] Kvashnina Y.A., Kvashnin A.G., Sorokin P.B. Investigation of new superhard carbon allotropes with promising electronic properties. J. Appl. Phys., 2013,114, P. 183708.
[27] Ordejon P., Artacho E., Soler J.M., Self-consistent order-N denisity-functional calculations for very large systems. Phys. Rev. B, 1996, 53, P. R10441.
[28] Soler J.M., Artacho E., Gale J.D., Garcia A., Junquera J., Ordejon P., Sanchez-Portal D. The SIESTA method for ab initio order-N materials simulation. J. Phys. Condens. Matter, 2002, 14, P. 2745-2780.
[29] Porezag D., Frauenheim T., Köhler T., Seifert G., Kaschner R., Construction of tight-binding-like potentials on the basis of density-functional theory: application to carbon. Phys. Rev. B, 1995, 51, P. 12947.
[30] Koster A.M., Flores R., Geudtner G., Goursot A., Heine T., Patchkovskii S., Reveles J.U., Vela A., Salahub D., deMon, Version 1.2; NRC: Ottawa, 2004.
[31] Caglioti G., Paoletti A., Ricci F.P. Choice of collimators for a crystal spectrometer for neutron diffraction. Nucl. Instr. Meth., 1958, 3, P. 223228.
[32] Baerlocher Ch., McCusker L.B., Olson D.H. Atlas of zeolite framework types. Elsevier, Amsterdam, 2007, 404 p.
[33] http://www.ioffe.ru/SVA/NSM/Semicond/Diamond/
[34] Pugh S.F. Relations between the elastic moduli and plastic properties of polycrystalline pure metals. Phil. Mag., 1954, 45, P. 822-843.
[35] Zener C.M., Siegel S. Elasticity and anelasticity of metals. J. Phys. Chem., 1949, 53, P. 1468.
[36] Teter D.M. Computational alchemy: the search for new superhard materials. MRS Bull, 1998, 23, P. 22-27.
[37] Chen X.Q., Niu H., Li D., Li Y. Modeling hardness of polycrystalline materials and bulk metallic glasses. Intermetallics, 2011, 19, P. 12751281.
[38] Shumilova T.G. Mineralogy of natural carbon (in Russian). Ural Branch of Russian Academy of sciences, Ekaterinburg, 2003, 316 p.
[39] Humphrey W., Dalke A., Schulten K., VMD: visual molecular dynamics. J. Mol. Graphics, 1996, 14, P. 33-38.
