Next Article in Journal
Dedifferentiated Human Articular Chondrocytes Redifferentiate to a Cartilage-Like Tissue Phenotype in a Poly(ε-Caprolactone)/Self-Assembling Peptide Composite Scaffold
Next Article in Special Issue
A Stress-Induced Martensitic Transformation in Aged Ti49Ni51 Alloy after High-Velocity Impact
Previous Article in Journal
Oil-in-Water Self-Assembled Synthesis of Ag@AgCl Nano-Particles on Flower-like Bi2O2CO3 with Enhanced Visible-Light-Driven Photocatalytic Activity
Previous Article in Special Issue
The Lateral Compressive Buckling Performance of Aluminum Honeycomb Panels for Long-Span Hollow Core Roofs
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Reinvestigation of a Superhard Tetragonal sp3 Carbon Allotrope

Faculty of Information Engineering and Automation, Kunming University of Science and Technology, Kunming 650051, China
*
Author to whom correspondence should be addressed.
Materials 2016, 9(6), 484; https://doi.org/10.3390/ma9060484
Submission received: 5 May 2016 / Revised: 8 June 2016 / Accepted: 13 June 2016 / Published: 17 June 2016
(This article belongs to the Special Issue Computational Multiscale Modeling and Simulation in Materials Science)

Abstract

:
I 4 ¯ –carbon was first proposed by Zhang et al., this paper will report regarding this phase of carbon. The present paper reports the structural and elastic properties of the three-dimensional carbon allotrope I 4 ¯ –carbon using first-principles density functional theory. The related enthalpy, elastic constants, and phonon spectra confirm that the newly-predicted I 4 ¯ –carbon is thermodynamically, mechanically, and dynamically stable. The calculated mechanical properties indicate that I 4 ¯ –carbon has a larger bulk modulus (393 GPa), shear modulus (421 GPa), Young’s modulus (931 GPa), and hardness (55.5 GPa), all of which are all slightly larger than those of c-BN. The present results indicate that I 4 ¯ –carbon is a superhard material and an indirect-band-gap semiconductor. Moreover, I 4 ¯ –carbon shows a smaller elastic anisotropy in its linear bulk modulus, shear anisotropic factors, universal anisotropic index, and Young’s modulus.
PACS:
61.50.-f; 62.20.-x; 63.20.dk; 71.20.-b

1. Introduction

The group 14 elements, such as carbon, silicon, and germanium, have attracted much interest and have been extensively studied [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20]. Carbon is found on the Earth mainly in the form of graphite and diamond. The quest for carbon materials with desired properties is of great interest in both fundamental science and advanced technology. One of the most famous carbon materials is graphene. Recently, some scholars of great research found several carbon allotropes with low-energy metastable structures, such as monoclinic M-carbon [1,2], F-carbon [3], orthorhombic W-carbon [4], Z-carbon [5], H-carbon and S-carbon [6], C-carbon [7], Imma-carbon [8,9], M585-carbon [10], T12-carbon [11], C2/m-16 carbon [12], P2221-carbon [13], and Cco-carbon [14]. We easily found that these carbon allotropes, with low-energy metastable structures with sp3 hybridization, could possibly explain the superhard property of materials. Hardness is an important property that determines many of the technological applications of materials. The mechanical properties of some carbon allotropes (such as T-carbon [21] and Y-carbon [22]) are not excellent, as discussed above, but can be improved through modulation [23,24]. At the moment, there are a multitude of metastable carbon allotropes that have been predicted using the quantum-chemical methods of calculation, e.g., evolutionary metadynamics technique [25,26,27,28], and Universal Structure Prediction: Evolutionary Xtallography (USPEX) [29,30,31,32], or post-graphite superhard phase synthesized by cold-compressing graphite [33,34,35,36,37,38,39]. Compressed graphite appears as an initial material for transition to different superhard phases. An ocean of possible carbon modifications can be obtained by combining four-, five-, six-, seven-, even eight-membered carbon rings. The discussion above of carbon allotropes basically adopts this structure. Wang et al. [40] found a dynamically stable and energetically favourable carbon allotrope, i.e., J-carbon. It has a wider band gap (5.89 eV within LDA-HSE06) than that of diamond and a larger bulk modulus (395 GPa). Currently, Z-carbon [33,34] is the most stable and hardest material predicted by the theoretical methods compared with other theoretical structures, but it has not exceeded diamond’s bulk modulus and hardness. The inclusion in the Z-carbon of additional diamond blocks generates another superhard carbon allotrope family, investigated in [41]. Li et al. [42] predicted a new cubic carbon allotrope, namely, C96 carbon. Unfortunately, its mechanical properties are not excellent and it cannot be used as a potential superhard material. Recently, another sp3 carbon allotrope, i.e., diamond nanothread, was synthesized experimentally; Fitzgibbons et al. [43] and Zhan et al. [44] then reported regarding its mechanical properties. These nanothreads show extraordinary properties, such as strength and stiffness, higher than those of sp2 carbon nanotubes.
Moreover, I 4 ¯ –carbon was first predicted in [45]. However, the physical properties, i.e., the structural, mechanical, and electronic properties, of I 4 ¯ –carbon were not studied in [45], nor in other studies. The physical properties of a novel superhard carbon allotrope with tetragonal I 4 ¯ symmetry (16 atoms/cell) will be detailed in this paper using first-principles calculations.

2. Materials and Methods

The calculations were performed using the density functional theory (DFT) [46,47] by using the Cambridge Serial Total Energy Package (CASTEP) code [48]. The electron–ionic core interaction was represented using the ultrasoft pseudopotentials [49]. The structural optimizations were conducted using the Broyden–Fletcher–Goldfarb–Shanno (BFGS) minimization [50]. The calculations were performed using the local density approximation (LDA) [51,52] and generalized gradient approximation (GGA) in the form of the Perdew–Burke–Ernzerhof (PBE) [53] exchange correlation potential. The total energy convergence tests showed that convergence to within 1 m∙eV/atom was achieved using the above calculation parameters. The electron and core interactions were included by using the ultrasoft pseudopotentials method, with a plane-wave energy cutoff energy of 400 eV, where 2s22p2 are treated as valence electrons for C. Highly dense k-point [54] sampling with a grid spacing of less than 2π × 0.025 Å−1 (7 × 7 × 7 for I 4 ¯ –carbon) in the Brillouin zone was used. The self-consistent convergence of the total energy was 5 × 10−6 eV/atom; the maximum force on the atom was 0.01 eV/Å; the maximum ionic displacement was within 5 × 10−4 Å; the maximum stress was within 0.02 GPa.

3. Results and Discussion

3.1. Structural Properties

There is a new tetragonal carbon phase, namely, I 4 ¯ –carbon, which belongs to the I 4 ¯ space group. The crystal structure of I 4 ¯ –carbon is shown in Figure 1. Different atomic positions are denoted by different colours of spheres. There are five inequivalent carbon atoms in its conventional cell, located at C1: 8g (−0.1219, −0.2709, −0.6225), C2: 8g (0.4044, −0.1809, 0.0492), C3: 8g (0.4812, 0.2514, −0.3772), C4: 2a (0.5000, 0.5000, 0.50000) and C5: 2d (0.0000, 0.5000, 0.7500). The optimized lattice parameters within the GGA and LDA level of I 4 ¯ –carbon at ambient pressure are listed in Table 1. As listed, a = 5.5628 (5.5017) Å and c = 5.5082 (5.4471) Å are within the GGA (LDA) level and c is slightly smaller than a. These lattice parameters are in excellent agreement with the theoretical values predicted by Zhang et al. [45]. The results within GGA are closer to the theoretical values predicted by Zhang et al. [45]; thus, the following discussions use the results of the GGA level. I 4 ¯ –carbon has seven bond lengths, which are also listed in Table 1. Between C1 and C2, C1 and C3, there are two different bond lengths, i.e., 1.5586 and 1.5224 Å, 1.5555 and 1.5207 Å, respectively, within the GGA level. The average bond length of I 4 ¯ –carbon is 1.5510 Å, while it is 1.567 Å and 1.551 Å for C2/m–16 carbon and M carbon, respectively. They are all slightly greater than that of diamond (1.535 Å). The hardness of I 4 ¯ –carbon, C2/m–16 carbon, M carbon, c-BN and diamond is calculated by using Lyakhov and Oganov’s model [55]. The hardness of I 4 ¯ –carbon is 55.5 GPa, which is slightly smaller than that of C2/m–16 carbon (59.5 GPa) [12] and M carbon (66.6 GPa), approximately half of that of diamond (89.7 GPa), and slightly larger than that of c-BN (49.9 GPa); the hardness of diamond is very close to the result of Ref [55] (91.0 GPa) using the Lyakhov and Oganov’s model. The hardness of I 4 ¯ –carbon is 83.0 GPa in Ref [45], the main reason for this situation is that the empirical formula may over or underestimate the value of the material’s hardness. Most researchers agree that “superhard” materials are those with Hv exceeding 40 GPa. Although there are slight differences between the results of the empirical models above, all of them greatly exceed 40 GPa, indicating that I 4 ¯ –carbon is a superhard material.

3.2. Mechanical Properties

For the tetragonal I 4 ¯ –carbon, seven independent elastic constants Cij were determined within the GGA level from the stress of the strained structure with a finite strain. The elastic constants of I 4 ¯ –carbon under different pressures are listed in Table 2. From Table 2, one can find that the mechanical stability of I 4 ¯ –carbon satisfies Born’s criterion for a tetragonal crystal [56]: Cii > 0, i = 1, 3, 4, 6; C11C12 > 0; C11 + C33 − 2C13 > 0 and 2(C11 + C12) + C33 + 4C13 > 0, indicating that I 4 ¯ –carbon is mechanically stable at 0–100 GPa. Moreover, almost all elastic constants increase with increasing pressure. C11 increases the fastest (dC11/dP = 4.57), followed by C33 (dC33/dP = 4.24); C44 is the slowest (dC44/dP = 0.63). C16 decreases with increasing pressure (dC16/dP = −0.35). Furthermore, it is important to explore the dynamic stability and thermodynamic stability for further experimental synthesis. The phonon spectra and the enthalpies as a function of pressure are shown in Figure 2a–d. Clearly, no imaginary mode is shown in the phonon spectra. This strongly confirms that the I 4 ¯ –carbon with the I 4 ¯ structure is dynamically stable; thus, it is a new metastable phase in the carbon family. The enthalpy of different carbon allotropes is quantified in terms of the following formation enthalpies formula: ΔH = Hcarbon allotropes/n1 − Hdiamond/n2, where n denotes the number of atoms in a conventional cell for carbon allotropes or graphite. As shown in Figure 2c,d, the calculated formation enthalpies of I 4 ¯ –carbon are slightly smaller than those of M-carbon below 32.25 GPa. Bulk modulus B and shear modulus G are calculated by using the Voigt–Reuss–Hill approximation. It is known that the Voigt bound is obtained by using the average polycrystalline modulus, based on an assumption of uniform strain throughout a polycrystalline, and is the upper limit of the actual effective modulus, while the Reuss bound is obtained by assuming a uniform stress and is the lower limit of the actual effective modulus. The arithmetic average of Voigt and Reuss bounds is referred to as the Voigt–Reuss–Hill approximations. Young’s modulus E and Poisson’s ratio v are obtained using the following equations: E = 9BG/(3B + G), v = (3B − 2G)/(6B + 2G), respectively. The calculated results of elastic modulus for I 4 ¯ –carbon are also shown in Table 2. The bulk modulus, shear modulus, Poisson’s ratio, and Young’s modulus all increase with increasing pressure. The bulk modulus increases the fastest, while shear modulus increases the slowest, i.e., dB/dP = 3.30 and dG/dP = 0.86, respectively. Pugh [57] proposed the ratio of bulk to shear modulus (B/G) and Lewandowski [58] proposed Poisson’s ratio as indications of a ductile versus brittle characteristic. B/G > 1.75 and v > 0.26 for a solid material represents ductility, while B/G < 1.75 and v < 0.26 usually indicate brittleness. From 0 to 100 GPa, I 4 ¯ –carbon is brittle; namely, the brittleness of I 4 ¯ –carbon decreases with increasing pressure. Compared with other carbon allotropes (P2221–carbon: 0.872 [13]; C2/m–16 carbon: 0.867 [12]; Imma–carbon: 0.858 [12]; diamond: 0.831 [12]), I 4 ¯ –carbon (0.933) has the least brittleness at ambient pressure.
Incompressibility can reflect the size of the bulk modulus from the other side. The V/V0 of I 4 ¯ –carbon, C2/m–16 carbon, M–carbon, P2221–carbon, c–BN, and diamond as a function of pressure are shown in Figure 3a. From Figure 3a, it is clear that c–BN shows a weaker incompressibility than that of the other carbon allotropes, while I 4 ¯ –carbon exhibits the weakest incompressibility among the carbon allotropes. The order of incompressibility for carbon allotropes and c-BN is diamond > P2221–carbon > M–carbon > C2/m–16 carbon > I 4 ¯ –carbon > c–BN. The bulk moduli of carbon allotropes and c-BN are also in this order: diamond (439 GPa) > P2221–carbon (409 GPa) > M–carbon (399 GPa) > C2/m–16 carbon (398 GPa) > I 4 ¯ –carbon (393 GPa) > c–BN (370 GPa). Regarding the incompressibility of I 4 ¯ –carbon in detailed investigations, the lattice parameters a/a0 and c/c0 as functions of pressure are illustrated in Figure 3b. Regarding the incompressibility of lattice parameters a/a0 and c/c0, they are very close. The incompressibility of lattice parameter c is slightly larger than that of a. By fitting the calculated data using the polynomial fitting methods, the following relationships at 0 GPa and 0 K are obtained:
a/a0 = 1.96545 × 10−6P2 − 7.85386 × 10−4P + 0.99970
c/c0 = 1.93918 × 10−6P2 − 7.79164 × 10−4P + 0.99965
where the unit of pressure is GPa. By fitting using the polynomial, it is clear that the incompressibility of lattice parameter c is slightly better than that of a.
The defined Debye temperature parameter, which can be interpreted as the temperature at which the highest-frequency mode (and hence every mode) is excited, is of great importance, allowing us to predict the heat capacity at any temperature in an atomic solid. The formula for evaluating the Debye temperature is given by [59]:
Θ D = h k B [ 3 n 4 π ( N A ρ M ) ] 1 3 v m
where h denotes Planck’s constant, kB denotes Boltzmann’s constant, n denotes the number of atoms per formula unit, NA denotes Avogadro’s number, M denotes molar mass, and ρ denotes density. The average sound velocity vm is described by the following formula:
v m = [ 1 3 ( 2 v t 3 + 1 v l 3 ) ] 1 3
The transverse and longitudinal elastic wave velocity, vt and vl, can be obtained using Navier’s equation, as follows [60]:
v t = G ρ   and   v l = 3 B + 4 G 3 ρ
where B and G are the bulk modulus and shear modulus, respectively. Based on Equation (3), the Debye temperature of I 4 ¯ –carbon ΘD = 2024 K. The Debye temperature of I 4 ¯ –carbon as a function of pressure is plotted in Figure 4a. It is found that the Debye temperature increases with increasing pressure. Furthermore, the greater the pressure, the more slowly the Debye temperature increases. The calculated results of the density, Debye temperature, transverse and longitudinal elastic wave velocity vt and vl, and the mean sound velocity vm at different pressures are listed in Table 3. The density and vl increase with increasing pressure. With the increase of pressure, the transverse longitudinal elastic wave velocity vt shows no monotonic increase or decrease. However, the mean sound velocity vm increases with increasing pressure from 0 to 80 GPa and decreases with increasing pressure from 80 to 100 GPa. The main reason is that sometimes vt increases with increasing pressure and sometimes decreases, i.e., there is no law.

3.3. Anisotropic Properties

It is well known that the anisotropy of elasticity is an important implication in engineering science and crystal physics [61], hence, it is worthwhile to investigate the elastic anisotropy of materials. The shear anisotropic factors provide a measure of Ba, Bb, and Bc, i.e., the bulk modulus along the a, b and c axes, respectively, which can be calculated by using the following equations:
B a = a d P d a = Λ 1 + α + β ,
B b = b d P d b = B a α ,
B c = c d P d c = B a β ,
Λ = C 11 + 2 C 12 α + C 22 α 2 + C 33 β 2 + 2 C 13 β ( 1 + α ) ,
α = ( C 11 C 12 ) ( C 22 C 12 ) ,
β = ( C 22 C 12 ) ( C 11 C 13 ) ( C 11 C 12 ) ( C 13 C 12 ) ( C 22 C 12 ) ( C 33 C 13 )
The calculated Ba, Bb, and Bc at different pressures are shown in Figure 5a. Due to the lattice constant a being equal to b, the bulk modulus along the a-axis is equal to that along the b axis, in other words, Ba = Bb. It is clear that Ba, Bb, and Bc increase with increasing pressure and that Ba, Bb (87.69%) increase more than Bc (76.50%). The calculated directional bulk modulus suggests that it is the largest along the a-axis and the smallest along the c-axis, indicating that the compressibility along the c-axis is the smallest, while along the a-axis, is the largest. This is in accordance with the relationships between the ratios a/a0, c/c0 and pressure, as shown in Figure 3b. Therefore, the anisotropy of the linear bulk modulus should also be considered. The anisotropy of the bulk modulus along the a-axis and c-axis with respect to the b-axis can be estimated by:
A B a = B a B b ,   and   A B c = B c B b
The anisotropy factors of the bulk modulus along the a-axis and c-axis for I 4 ¯ –carbon at T = 0 K as a function of pressure are shown in Figure 5b. Note that a value of 1.0 indicates elastic isotropy and any departure from 1.0 represents elastic anisotropy. The anisotropy of the bulk modulus along the a-axis and c-axis with respect to the b-axis shows that ABa is elastic isotropy and ABc is elastic anisotropy.
After discussing the anisotropy of the bulk modulus, we now discuss the anisotropy of the shear modulus. The shear anisotropic factors provide a measure of the degree of anisotropy in the bonding between atoms in different planes. The shear anisotropic factor for the (1 0 0) shear planes between [0 1 1] and [0 1 0] directions and the (0 1 0) shear planes between [1 0 1] and [0 0 1] directions is [62]:
A 1 = 4 C 44 C 11 + C 33 2 C 13 ,
For the (0 1 0) shear planes between [1 0 1] and [0 0 1] directions, it is [60]:
A 2 = 4 C 55 C 22 + C 33 2 C 23 ,
For the (0 0 1) shear planes between [1 1 0] and [0 1 0] directions, it is [60]:
A 3 = 4 C 66 C 11 + C 22 2 C 12
For an isotropic crystal, the factors A1, A2, and A3 are 1.0, while any value smaller or larger than 1.0 is a measure of the elastic anisotropy possessed by the materials. The anisotropy factors of I 4 ¯ –carbon at T = 0 K as a function of pressure are shown in Figure 5c. Due to the symmetry of the crystal structure, C11 = C22, C44 = C55 and C13 = C23; thus, A1 = A2. From Figure 5c, A1 first increases and then decreases with increasing pressure, while A3 increases monotonically with increasing pressure. I 4 ¯ –carbon shows elastic anisotropy in the bulk modulus and shear modulus. The universal anisotropic index (AU = 5GV/GR + BV/BR − 6) combines the shear modulus and bulk modulus to exhibit the anisotropy of the material. The AU of I 4 ¯ –carbon increases with pressure (see Figure 5d), indicating that the anisotropy of I 4 ¯ –carbon increases with increasing pressure. The AU of I 4 ¯ –carbon is 0.0369, which is smaller than that of C2/m–carbon (0.0766) and P2221–carbon (0.0526).
Young’s modulus is not always the same in all orientations of a material; it will change depending on the direction of the force vector. Engineers can use this directional phenomenon to their advantage in creating structures. The directional dependence of the anisotropy in Young’s modulus is calculated by using the Elastic Anisotropy Measures (ELAM) [63,64] code. The 2D representation of Young’s modulus for I 4 ¯ –carbon in the (001), (010), (100), (111), ( 00 1 ¯ ), and (101) planes is shown in Figure 6a–f, respectively. The black, red, and blue solid curves represent Young’s modulus for I 4 ¯ –carbon in the (001), (010), (100), (111), ( 00 1 ¯ ), and (101) planes at 0, 50 and 100 GPa, respectively. It is obvious that the anisotropy of Young’s modulus increases with increasing pressure in the (001), (010), (100), (111), ( 00 1 ¯ ), and (101) planes. To investigate the anisotropy of Young’s modulus in detail, we calculate the maximal and minimal values of Young’s modulus for I 4 ¯ –carbon in the (001), (010), (100), (111), ( 00 1 ¯ ), and (101) planes and for all possible directions, together with the ratio Emax/Emin, which are listed in Table 4. The conclusion from Table 4 is consistent with the conclusion we obtained from Figure 6, i.e., that the anisotropy of Young’s modulus increases with increasing pressure in the (001), (010), (100), (111), ( 00 1 ¯ ), and (101) planes. Regarding the distribution of Young's modulus in the (001) plane and ( 00 1 ¯ ) plane, the maximal and minimal values are the same, while in the (010) plane and (100) plane, they are also the same. The (001) plane and ( 00 1 ¯ ) plane exhibit the largest anisotropy of Young’s modulus for I 4 ¯ –carbon, while the (101) plane exhibits the smallest anisotropy of Young’s modulus. This situation also appears under high pressure. All of the degrees of special plane anisotropy are smaller than the overall performance of the material anisotropy because the minimal value of the entire Young’s modulus was not obtained. The anisotropy of Young’s modulus of I 4 ¯ –carbon (Emax/Emin = 1.133) is slightly smaller than that of C2/m-16 carbon (Emax/Emin = 1.275) and P2221-carbon (Emax/Emin = 1.173). Thus, I 4 ¯ –carbon shows a smaller elastic anisotropy in its linear bulk modulus, shear anisotropic factors, universal anisotropic index, and Young’s modulus.

3.4. Electronic Properties

The band structures of I 4 ¯ –carbon are shown in Figure 7a,b at 0 and 100 GPa, respectively. Semiconductor materials have direct or indirect band gaps. In direct band gap materials, the maximum of the valence band (VBM) and minimum of the conduction band (CBM) occur at the same value of the wave vector, while in indirect band gap materials, they occur for different values of the wave vector. It is obvious that I 4 ¯ –carbon is an indirect-band-gap and wide-band-gap semiconductor material. The VBM of I 4 ¯ –carbon is located at the M point (0.500, 0.500, 0.000), while the CBM is located from point (0.000, 0.000, 0.3667) through G (0.000, 0.000, 0.000) to the Z point (0.000, 0.000, 0.500), whether at ambient pressure or under high pressure. From Figure 7a,b, it is shown that the band gap increases with increasing pressure; more details are shown in Figure 4b; the band gap of I 4 ¯ –carbon increases from 5.19 eV to 5.59 eV. The greater the pressure, the less the band gap increases. To determine the reason why, we calculate the energy values of the CBM and VBM under different pressures. The calculated results are illustrated in Figure 7c. The energy of VBM for I 4 ¯ –carbon is 7.96 eV at 0 GPa and 10.58 eV at 100 GPa; ΔEVBM = 2.62 eV. The energy of CBM for I 4 ¯ –carbon is 13.15 eV at 0 GPa and 16.17 eV at 100 GPa; ΔECBM = 3.02 eV. ΔECBMΔEVBM = 0.40 eV, which is the same as with the band gap increase. In other words, as pressure increases, the conduction band and valence band energy will increase, but the energy of the conduction band will increase more than that of the valence band. The partial density of states (PDOS) of five inequivalent carbon atoms in I 4 ¯ –carbon at 0 and 100 GPa are shown in Figure 8. Figure 8a,c,e,g,i show the partial density of states of five inequivalent carbon atoms at 0 GPa; Figure 8b,d,f,h,j show the partial density of states of five inequivalent carbon atoms at 100 GPa. No matter the position of the carbon atoms, the partial density of states in the conduction band part is, with increasing pressure intensity, near low energy, and the valence band deviates from the Fermi level. The PDOS in the lower energy part is contributed mainly by the C-2s orbitals, while in the other energy part, is contributed mainly by the C-2p orbitals.

4. Conclusions

In conclusion, the low-energy metastable structures of I 4 ¯ –carbon have been systematically investigated based on the density functional theory. The related enthalpy, elastic constants and phonon spectra confirm that the newly-predicted I 4 ¯ –carbon is thermodynamically, mechanically, and dynamically stable. As a potential superhard material for engineering applications, the mechanical properties of I 4 ¯ –carbon indicate that it is a superhard material. I 4 ¯ –carbon shows a smaller anisotropy in its linear bulk modulus, shear anisotropic factors, universal anisotropic index, and Young’s modulus. The band structure shows that I 4 ¯ –carbon is an indirect-band-gap and wide-band-gap semiconductor material. Finally, the calculations of its mechanical properties reveal that I 4 ¯ –carbon possesses a high bulk and shear modulus as well as a low Poisson’s ratio and B/G ratio (<1.75). Moreover, I 4 ¯ –carbon has a larger Debye temperature (ΘD = 2024 K). Due to their wide-band-gap and higher bulk moduli, hardness, they are attractive for semiconductor device applications and superhard material with potential technological and industrial applications.

Acknowledgments

This work was supported by the Fund for Talents of Yunnan Province, China (Grant No. KKSY201403006) and the National Natural Science Foundation of China (No. 61564005).

Author Contributions

Mengjiang Xing and Binhua Li designed the project; Mengjiang Xing, Binhua Li and Zhengtao Yu performed the calculations; Mengjiang Xing, Zhengtao Yu and Qi Chen determined the results; Mengjiang Xing and Binhua Li wrote the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Oganov, A.R.; Glass, C.W. Crystal structure prediction using ab initio evolutionary techniques: Principles and applications. J. Chem. Phys. 2006, 124, 244704. [Google Scholar] [CrossRef] [PubMed]
  2. Li, Q.; Ma, Y.M.; Oganov, A.R.; Wang, H.B.; Wang, H.; Xu, Y.; Cui, T.; Mao, H.K.; Zou, G.T. Superhard Monoclinic Polymorph of Carbon. Phys. Rev. Lett. 2009, 102, 175506. [Google Scholar] [CrossRef] [PubMed]
  3. Tian, F.; Dong, X.; Zhao, Z.S.; He, J.L.; Wang, H.T. Superhard F-carbon predicted by ab initio particle-swarm optimization methodology. J. Phys. Condens. Matter 2012, 24, 165504. [Google Scholar] [CrossRef] [PubMed]
  4. Wang, J.T.; Chen, C.; Kawazoe, Y. Low-Temperature Phase Transformation from Graphite to sp3 Orthorhombic Carbon. Phys. Rev. Lett. 2011, 106, 075501. [Google Scholar] [CrossRef] [PubMed]
  5. Li, Z.P.; Gao, F.M.; Xu, Z.M. Strength, hardness, and lattice vibrations of Z-carbon and W-carbon: First-principles calculations. Phys. Rev. B 2012, 85, 144115. [Google Scholar] [CrossRef]
  6. He, C.Y.; Sun, L.Z.; Zhang, C.X.; Peng, X.Y.; Zhang, K.W.; Zhong, J.X. New superhard carbon phases between graphite and diamond. Solid State Commun. 2012, 152, 1560–1563. [Google Scholar] [CrossRef]
  7. Li, D.; Bao, K.; Tian, F.B.; Zeng, Z.W.; He, Z.; Liu, B.B.; Cui, T. Lowest enthalpy polymorph of cold-compressed graphite phase. Phys. Chem. Chem. Phys. 2012, 14, 4347–4350. [Google Scholar] [CrossRef] [PubMed]
  8. Wei, Q.; Zhang, M.G.; Yan, H.Y.; Lin, Z.Z.; Zhu, X.M. Structural, electronic and mechanical properties of Imma-carbon. EPL 2014, 107, 27007. [Google Scholar] [CrossRef]
  9. Liu, Y.M.; Lu, M.C.; Zhang, M. First-principles study of a novel superhard sp3 carbon allotrope. Phys. Lett. A 2014, 378, 3326. [Google Scholar] [CrossRef]
  10. He, C.Y.; Zhong, J.X. M585, a low energy superhard monoclinic carbon phase. Solid State Commun. 2014, 181, 24–27. [Google Scholar] [CrossRef]
  11. Zhao, Z.S.; Tian, F.; Dong, X.; Li, Q.; Wang, Q.Q.; Wang, H.; Zhong, X.; Xu, B.; Yu, D.L.; He, J.L.; et al. Tetragonal Allotrope of Group 14 Elements. J. Am. Chem. Soc. 2012, 134, 12362–12365. [Google Scholar] [CrossRef] [PubMed]
  12. Xing, M.J.; Li, B.H.; Yu, Z.T.; Chen, Q. C2/m-carbon: Structural, mechanical, and electronic properties. J. Mater. Sci. 2015, 50, 7104–7114. [Google Scholar] [CrossRef]
  13. Xing, M.J.; Li, B.H.; Yu, Z.T.; Chen, Q. Structural, Elastic, and Electronic Properties of a New Phase of Carbon. Commun. Theor. Phys. 2015, 64, 237–243. [Google Scholar] [CrossRef]
  14. Zhao, Z.S.; Xu, B.; Zhou, X.F.; Wang, L.M.; Wen, B.; He, J.L.; Liu, Z.Y.; Wang, H.T.; Tian, Y.J. Novel Superhard Carbon: C-Centered Orthorhombic C8. Phys. Rev. Lett. 2011, 107, 215502. [Google Scholar] [CrossRef] [PubMed]
  15. Guo, Y.G.; Wang, Q.; Kawazoe, Y.; Jena, P. A New Silicon Phase with Direct Band Gap and Novel Optoelectronic Properties. Sci. Rep. 2015, 5, 14342. [Google Scholar] [CrossRef] [PubMed]
  16. Fan, Q.Y.; Chai, C.C.; Wei, Q.; Yan, H.Y.; Zhao, Y.B.; Yang, Y.T.; Yu, X.H.; Liu, Y.; Xing, M.J.; Zhang, J.Q.; et al. Novel silicon allotropes: Stability, mechanical, and electronic properties. J. Appl. Phys. 2015, 118, 185704. [Google Scholar] [CrossRef]
  17. Lee, I.H.; Lee, J.Y.; Oh, Y.J.; Kim, S.; Chang, K.J. Computational search for direct band gap silicon crystals. Phys. Rev. B 2014, 90, 115209. [Google Scholar] [CrossRef]
  18. Fan, Q.Y.; Chai, C.C.; Wei, Q.; Yang, Y.T.; Yang, Q.; Chen, P.Y.; Xing, M.J.; Zhang, J.Q.; Yao, R.H. Prediction of novel phase of silicon and Si–Ge alloys. J. Solid State Chem. 2016, 233, 471–483. [Google Scholar] [CrossRef]
  19. De, A.; Pryor, C.E. Electronic structure and optical properties of Si, Ge and diamond in the lonsdaleite phase. J. Phys. Condens. Matter 2014, 26, 045801. [Google Scholar] [CrossRef] [PubMed]
  20. Wang, Q.Q.; Xu, B.; Sun, J.; Liu, H.Y.; Zhao, Z.S.; Yu, D.L.; Fan, C.Z.; He, J.L. Direct Band Gap Silicon Allotropes. J. Am. Chem. Soc. 2014, 136, 9826–9829. [Google Scholar] [CrossRef] [PubMed]
  21. Sheng, X.L.; Yan, Q.B.; Ye, F.; Zheng, Q.R.; Su, G. T-Carbon: A Novel Carbon Allotrope. Phys. Rev. Lett. 2011, 106, 155703. [Google Scholar] [CrossRef] [PubMed]
  22. Jo, J.Y.; Kim, B.G. Carbon allotropes with triple bond predicted by first-principle calculation: Triple bond modified diamond and T-carbon. Phys. Rev. B 2012, 86, 075151. [Google Scholar] [CrossRef]
  23. Srinivasu, K.; Ghosh, S.K. Electronic Structure, Optical Properties, and Hydrogen Adsorption Characteristics of Supercubane-Based Three-Dimensional Porous Carbon. J. Phys. Chem. C 2012, 116, 25015. [Google Scholar] [CrossRef]
  24. Li, D.; Tian, F.B.; Duan, D.F.; Zhao, Z.L.; Liu, Y.X.; Chu, B.H.; Sha, X.J.; Wang, L.; Liu, B.B.; Cui, T. Modulated T carbon-like carbon allotropes: An ab initio study. RSC Adv. 2014, 4, 17364–17369. [Google Scholar] [CrossRef]
  25. Martoňák, R.; Oganov, A.R.; Glass, C.W. Crystal structure prediction and simulations of structural transformations: Metadynamics and evolutionary algorithms. Phase Transit. 2007, 80, 277–298. [Google Scholar] [CrossRef]
  26. Zhu, Q.; Oganov, A.R.; Lyakhov, A.O. Evolutionary metadynamics: A novel method to predict crystal structures. Crystengcomm 2012, 14, 3596–3601. [Google Scholar] [CrossRef]
  27. Boulfelfel, S.E.; Zhu, Q.; Oganov, A.R. Novel sp3 forms of carbon predicted by evolutionary metadynamics and analysis of their synthesizability using transition path sampling. J. Superhard Mater. 2012, 34, 350–359. [Google Scholar] [CrossRef]
  28. Finkelstein, G.J.; Dera, P.K.; Jahn, S.; Oganov, A.R.; Holl, C.M.; Meng, Y.; Duffy, T.S. Phase transitions and equation of state of forsterite to 90 GPa from single-crystal X-ray diffraction and molecular modelling. Am. Mineral. 2014, 99, 35–43. [Google Scholar] [CrossRef]
  29. Oganov, A.R.; Glass, C.W. Evolutionary crystal structure prediction as a tool in materials design. J. Phys. Condens. Matter 2008, 20, 064210. [Google Scholar] [CrossRef] [PubMed]
  30. Oganov, A.R.; Ma, Y.M.; Lyakhov, A.O.; Valle, M.; Gatti, C. Evolutionary Crystal Structure Prediction as a Method for the Discovery of Minerals and Materials. Theor. Comput. Methods Miner. Phys. Geophys. Appl. Rev. Miner. Geochem. 2010, 71, 271–298. [Google Scholar] [CrossRef]
  31. Oganov, A.R.; Lyakhov, A.O.; Valle, M. How Evolutionary Crystal Structure Prediction Works—And Why. Accounts Chem. Res. 2011, 44, 227–237. [Google Scholar] [CrossRef] [PubMed]
  32. Lyakhov, A.O.; Oganov, A.R.; Stokes, H.T.; Zhu, Q. New developments in evolutionary structure prediction algorithm USPEX. Comput. Phys. Commun. 2013, 184, 1172–1182. [Google Scholar] [CrossRef]
  33. Selli, D.; Baburin, I.A.; Martonák, R.; Leoni, S. Superhard sp3 carbon allotropes with odd and even ring topologies. Phys. Rev. B 2011, 84, 161411. [Google Scholar] [CrossRef]
  34. Amsler, M.; Flores-Livas, J.A.; Lehtovaara, L.; Balima, F.; Ghasemi, S.A.; Machon, D.; Pailhès, S.; Willand, A.; Caliste, D.; Botti, S.; et al. Crystal Structure of Cold Compressed Graphite. Phys. Rev. Lett. 2012, 108, 065501. [Google Scholar] [CrossRef] [PubMed]
  35. Wang, J.T.; Chen, C.; Kawazoe, Y. Orthorhombic carbon allotrope of compressed graphite: Ab initio calculations. Phys. Rev. B 2012, 85, 033410. [Google Scholar] [CrossRef]
  36. He, C.Y.; Sun, L.Z.; Zhang, C.X.; Zhong, J.X. Two viable three-dimensional carbon semiconductors with an entirely sp2 configuration. Phys. Chem. Chem. Phys. 2013, 15, 680–684. [Google Scholar] [CrossRef] [PubMed]
  37. He, C.Y.; Sun, L.Z.; Zhang, C.X.; Peng, X.Y.; Zhang, K.W.; Zhong, J.X. Four superhard carbon allotropes: A first-principles study. Phys. Chem. Chem. Phys. 2012, 14, 8410–8414. [Google Scholar] [CrossRef] [PubMed]
  38. Kvashnina, Y.A.; Kvashnin, A.G.; Sorokin, P.B. Investigation of new superhard carbon allotropes with promising electronic properties. J. Appl. Phys. 2013, 114, 183708. [Google Scholar] [CrossRef]
  39. Xie, H.X.; Yin, F.X.; Yu, T. Mechanism for direct graphite-to-diamond phase transition. Sci. Rep. 2014, 4, 5930. [Google Scholar] [CrossRef] [PubMed]
  40. Wang, J.T.; Chen, C.F.; Kawazoe, Y. Phase conversion from graphite toward a simple monoclinic sp3-carbon allotrope. J. Chem. Phys. 2012, 137, 024502. [Google Scholar] [CrossRef] [PubMed]
  41. Zhou, R.; Zeng, X.C. Polymorphic Phases of sp3-hybridized carbon under cold compression. J. Am. Chem. Soc. 2012, 134, 7530–7538. [Google Scholar] [CrossRef] [PubMed]
  42. Li, D.; Tian, F.B.; Chu, B.H.; Duan, D.F.; Wei, S.L.; Lv, Y.Z.; Zhang, H.D.; Wang, L.; Lu, N.; Liu, B.B.; et al. Cubic C96: A novel carbon allotrope with a porous nanocube network. J. Mater. Chem. A 2015, 3, 10448–10452. [Google Scholar] [CrossRef]
  43. Fitzgibbons, T.C.; Guthrie, M.; Xu, E.S.; Crespi, V.H.; Davidowski, S.K.; Cody, G.D.; Alem, N.; Badding, J.V. Benzene-derived carbon nanothreads. Nat. Mater. 2015, 14, 43–47. [Google Scholar] [CrossRef] [PubMed]
  44. Zhan, H.F.; Zhang, G.; Tan, V.B.C.; Cheng, Y.; Bell, J.M.; Zhang, Y.W.; Gu, Y.T. From Brittle to Ductile: A Structure Dependent Ductility of Diamond Nanothread. Nanoscale 2016, 8, 11177–11184. [Google Scholar] [CrossRef] [PubMed]
  45. Zhang, X.X.; Wang, Y.C.; Lv, J.; Zhu, C.Y.; Li, Q.; Zhang, M.; Li, Q.; Ma, Y.M. First-principles structural design of superhard materials. J. Chem. Phys. 2013, 138, 114101. [Google Scholar] [CrossRef] [PubMed]
  46. Hohenberg, P.; Kohn, W. Inhomogeneous electron gas. Phys. Rev. 1964, 136, B864. [Google Scholar] [CrossRef]
  47. Kohn, W.; Sham, L.J. Self-consistent equations including exchange and correlation effects. Phys. Rev. 1965, 140, A1133. [Google Scholar] [CrossRef]
  48. Clark, S.J.; Segall, M.D.; Pickard, C.J.; Hasnip, P.J.; Probert, M.I.J.; Refson, K.; Payne, M.C. First principles methods using CASTEP. Z. Kristallogr. 2005, 220, 567–570. [Google Scholar] [CrossRef] [Green Version]
  49. Vanderbilt, D. Soft self-consistent pseudopotentials in a generalized eigenvalue formalism. Phys. Rev. B 1990, 41, 7892R. [Google Scholar] [CrossRef]
  50. Pfrommer, B.G.; Côté, M.; Louie, S.G.; Cohen, M.L. Relaxation of crystals with the quasi-newton method. J. Comput. Phys. 1997, 131, 233–240. [Google Scholar] [CrossRef]
  51. Ceperley, D.M.; Alder, B.J. Ground state of the electron gas by a stochastic method. Phys. Rev. Lett. 1980, 45, 566. [Google Scholar] [CrossRef]
  52. Perdew, J.P.; Zunger, A. Self-interaction correction to density-functional approximations for many-electron systems. Phys. Rev. B 1981, 23, 5048–5079. [Google Scholar] [CrossRef]
  53. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef] [PubMed]
  54. Monkhorst, H.J.; Pack, J.D. Special points for Brillouin-zone integrations. Phys. Rev. B 1976, 13, 5188–5192. [Google Scholar] [CrossRef]
  55. Lyakhov, A.O.; Oganov, A.R. Evolutionary search for superhard materials: Methodology and applications to forms of carbon and TiO2. Phys. Rev. B 2011, 84, 092103. [Google Scholar] [CrossRef]
  56. Wu, Z.J.; Zhao, E.J.; Xiang, H.P.; Hao, X.F.; Liu, X.J.; Meng, J. Crystal structures and elastic properties of superhard IrN2 and IrN3 from first principles. Phys. Rev. B 2007, 76, 054115. [Google Scholar] [CrossRef]
  57. Pugh, S.F. XCII. Relations between the elastic moduli and the plastic properties of polycrystalline pure metals. Lond. Edinb. Dublin Philos. Mag. J. Sci. Ser. 7 1954, 45, 823–843. [Google Scholar] [CrossRef]
  58. Lewandowski, J.J.; Wang, W.H.; Greer, A.L. Intrinsic plasticity or brittleness of metallic glasses. Philos. Mag. Lett. 2005, 85, 77–87. [Google Scholar] [CrossRef]
  59. Anderson, O.L. A simplified method for calculating the debye temperature from elastic constants. J. Phys. Chem. Solids 1963, 24, 909–917. [Google Scholar] [CrossRef]
  60. Panda, K.B.; Ravi, K.S. Determination of elastic constants of titanium diboride (TiB2) from first principles using FLAPW implementation of the density functional theory. Comput. Mater. Sci. 2006, 35, 134–150. [Google Scholar] [CrossRef]
  61. Fan, Q.Y.; Wei, Q.; Yan, H.Y.; Zhang, M.G.; Zhang, Z.X.; Zhang, J.Q.; Zhang, D.Y. Elastic and electronic properties of Pbca-BN: First-principles calculations. Comput. Mater. Sci. 2014, 85, 80–87. [Google Scholar] [CrossRef]
  62. Connetable, D.; Thomas, O. First-principles study of the structural, electronic, vibrational, and elastic properties of orthorhombic NiSi. Phys. Rev. B 2009, 79, 094101. [Google Scholar] [CrossRef]
  63. Marmier, A.; Lethbridge, Z.A.D.; Walton, R.I.; Smith, C.W.; Parker, S.C.; Evans, K.E. ElAM: A computer program for the analysis and representation of anisotropic elastic properties. Comput. Phys. Commun. 2010, 181, 2102–2115. [Google Scholar] [CrossRef] [Green Version]
  64. Fan, Q.Y.; Wei, Q.; Chai, C.C.; Yan, H.Y.; Zhang, M.G.; Lin, Z.Z.; Zhang, Z.X.; Zhang, J.Q.; Zhang, D.Y. Structural, mechanical, and electronic properties of P3m1-BCN. J. Phys. Chem. Solids 2015, 79, 89–96. [Google Scholar] [CrossRef]
Figure 1. The crystal structure of I 4 ¯ –carbon.
Figure 1. The crystal structure of I 4 ¯ –carbon.
Materials 09 00484 g001
Figure 2. Phonon spectra for I 4 ¯ –carbon at 0 GPa (a) and 100 GPa (b); enthalpies of I 4 ¯ –carbon and other carbon allotropes relative to graphite as a function of pressure (c,d).
Figure 2. Phonon spectra for I 4 ¯ –carbon at 0 GPa (a) and 100 GPa (b); enthalpies of I 4 ¯ –carbon and other carbon allotropes relative to graphite as a function of pressure (c,d).
Materials 09 00484 g002
Figure 3. The primitive cell volume V/V0 as a function of pressure for I 4 ¯ –carbon, P2221-carbon, C2/m-16 carbon, M-carbon, c-BN and diamond (a); lattice constants a/a0, c/c0 compression of I 4 ¯ –carbon as a function of pressure (b).
Figure 3. The primitive cell volume V/V0 as a function of pressure for I 4 ¯ –carbon, P2221-carbon, C2/m-16 carbon, M-carbon, c-BN and diamond (a); lattice constants a/a0, c/c0 compression of I 4 ¯ –carbon as a function of pressure (b).
Materials 09 00484 g003
Figure 4. Debye temperature (a) and band gap (b) of I 4 ¯ –carbon as a function of pressure.
Figure 4. Debye temperature (a) and band gap (b) of I 4 ¯ –carbon as a function of pressure.
Materials 09 00484 g004
Figure 5. Linear bulk modulus Ba, Bb, and Bc of I 4 ¯ –carbon at 0 K as a function of pressure (a); anisotropy factors of I 4 ¯ –carbon at 0 K as a function of pressure, (b) ABb and ABc; (c) A1, A2 and A3; (d) AU.
Figure 5. Linear bulk modulus Ba, Bb, and Bc of I 4 ¯ –carbon at 0 K as a function of pressure (a); anisotropy factors of I 4 ¯ –carbon at 0 K as a function of pressure, (b) ABb and ABc; (c) A1, A2 and A3; (d) AU.
Materials 09 00484 g005
Figure 6. 2D representation of Young’s modulus in the (001), (010), (100), (111), ( 00 1 ¯ ), and (101) planes for I 4 ¯ –carbon, shown in (af), respectively.
Figure 6. 2D representation of Young’s modulus in the (001), (010), (100), (111), ( 00 1 ¯ ), and (101) planes for I 4 ¯ –carbon, shown in (af), respectively.
Materials 09 00484 g006
Figure 7. Electronic band structure for the I 4 ¯ –carbon at ambient pressure (a) and 100 GPa (b); the calculated energy values of CBM and VBM under different pressures (c).
Figure 7. Electronic band structure for the I 4 ¯ –carbon at ambient pressure (a) and 100 GPa (b); the calculated energy values of CBM and VBM under different pressures (c).
Materials 09 00484 g007
Figure 8. Partial density of state for I 4 ¯ –carbon; (a,c,e,g,i) are at ambient pressure, while (b,d,f,h,j) are at 100 GPa.
Figure 8. Partial density of state for I 4 ¯ –carbon; (a,c,e,g,i) are at ambient pressure, while (b,d,f,h,j) are at 100 GPa.
Materials 09 00484 g008aMaterials 09 00484 g008b
Table 1. The calculated lattice parameters (in Å) and bond lengths (in Å) for I 4 ¯ –carbon.
Table 1. The calculated lattice parameters (in Å) and bond lengths (in Å) for I 4 ¯ –carbon.
MethodsacdC1-C2dC1-C3dC2-C3dC2-C5dC3-C4
GGA5.56285.50821.55861.55551.56941.58711.5430
1.52241.5207
LDA5.50175.44711.54351.54081.55311.56831.5253
1.50361.5039
Table 2. The calculated elastic constants (in GPa) and elastic modulus (in GPa) of I 4 ¯ –carbon at different pressures.
Table 2. The calculated elastic constants (in GPa) and elastic modulus (in GPa) of I 4 ¯ –carbon at different pressures.
PC11C12C13C16C33C44C66BG
0926133101−281017398472393421
101015175142−351028408493441434
201034189153−401081422511460446
301079220182−431132428528495455
401128248208−461181436544529466
501174274235−491229442559563475
601220304261−531275449573596484
701263332287−561320454586629491
801306360303−591364458598661499
901347387340−621403462608692504
1001383419366−631441461614723507
Table 3. The calculated density (ρ, in g/cm3), the longitudinal, transverse, and mean elastic wave velocity (vs, vp, vm, in m/s), and the Debye temperature (ΘD, in K) for I 4 ¯ –carbon.
Table 3. The calculated density (ρ, in g/cm3), the longitudinal, transverse, and mean elastic wave velocity (vs, vp, vm, in m/s), and the Debye temperature (ΘD, in K) for I 4 ¯ –carbon.
PρvlvtvmΘD
03.276417,06911,33612,3982024
103.356417,43411,37312,4672051
203.433217,52911,40212,5032073
303.505917,73011,39512,5132089
403.575017,93911,41512,5492109
503.64121812411417125662124
603.70451830311427125892141
703.76541846511424125982154
803.82431861711417126022166
903.88141875111400125942176
1003.93651885211345125472178
Table 4. The maximum and minimum values of Young’s modulus (in GPa) in different planes and for pressures for I 4 ¯ –carbon.
Table 4. The maximum and minimum values of Young’s modulus (in GPa) in different planes and for pressures for I 4 ¯ –carbon.
PlanePMaximumMinimumEmax/Emin
(001)09998871.126
50125010541.186
100142711681.222
(010)09988891.123
50115310531.095
100129211501.123
(100)09988891.123
50115310531.095
100129211501.123
(111)09928881.117
50123310451.180
100140211391.231
( 00 1 ¯ )09998871.126
50125010541.186
100142711681.222
(101)09498851.072
50116010501.105
100129811471.132
All09998821.133
50125010421.200
100142711381.254

Share and Cite

MDPI and ACS Style

Xing, M.; Li, B.; Yu, Z.; Chen, Q. A Reinvestigation of a Superhard Tetragonal sp3 Carbon Allotrope. Materials 2016, 9, 484. https://doi.org/10.3390/ma9060484

AMA Style

Xing M, Li B, Yu Z, Chen Q. A Reinvestigation of a Superhard Tetragonal sp3 Carbon Allotrope. Materials. 2016; 9(6):484. https://doi.org/10.3390/ma9060484

Chicago/Turabian Style

Xing, Mengjiang, Binhua Li, Zhengtao Yu, and Qi Chen. 2016. "A Reinvestigation of a Superhard Tetragonal sp3 Carbon Allotrope" Materials 9, no. 6: 484. https://doi.org/10.3390/ma9060484

APA Style

Xing, M., Li, B., Yu, Z., & Chen, Q. (2016). A Reinvestigation of a Superhard Tetragonal sp3 Carbon Allotrope. Materials, 9(6), 484. https://doi.org/10.3390/ma9060484

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop