Next Article in Journal
The Use of Boron Fibers and Particles for Creating Functionally Graded Material Based on Ti64 Using the Laser Additive Manufacturing Method
Next Article in Special Issue
First-Principles Study on the Mechanical Properties of Gd-Doped BCZT Ceramics
Previous Article in Journal
Colloidal Synthesis and Ultraviolet Luminescence of Rb2AgI3 Nanocrystals
Previous Article in Special Issue
Si-Doped Nitrogenated Holey Graphene (C2N) as a Promising Gas Sensor for O-Containing Volatile Organic Compounds (VOCs) and Ammonia
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Crystals
Volume 13
Issue 7
10.3390/cryst13071111
crystals-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Crystal Chemistry and First Principles Studies of Novel Superhard Tetragonal C7, C5N2, and C3N4

by
Samir F. Matar
1 and
Vladimir L. Solozhenko
2,*
1
CMMS, Lebanese German University, Jounieh P.O. Box 206, Lebanon
2
LSPM–CNRS, Université Sorbonne Paris Nord, 93430 Villetaneuse, France
*
Author to whom correspondence should be addressed.
Crystals 2023, 13(7), 1111; https://doi.org/10.3390/cryst13071111
Submission received: 20 June 2023 / Revised: 11 July 2023 / Accepted: 12 July 2023 / Published: 17 July 2023
(This article belongs to the Special Issue First Principles Calculation for Crystalline Materials)
Browse Figures
Versions Notes

Abstract

:
Tetragonal C7, C5N2, and C3N4, characterized by mixed tetrahedral and trigonal atomic hybridizations, have been devised based on crystal chemistry rationale and structural optimization calculations within density functional theory (DFT). Substitution of C(sp2) and C(sp3) in C7 for nitrogen yields α-C5N2 and β-C5N2, respectively, both of which are superhard, cohesive, and stable mechanically (elastic properties) and dynamically (phonon band structures). tet-C3N4 with both nitrogen sites within the C7 structure was found to be cohesive and classified as ductile with a Vickers hardness of 65 GPa. Due to the delocalization of π electrons of the sp2-like hybridized atoms, metallic behavior characterizes all four phases.

1. Introduction

Since the 1950s, interest has focused on carbon allotropes, especially those that resemble diamond in hardness and other electronic properties, such as diamond at the nanoscale for applications in electronics [1]. Expanding to carbon’s neighbors, B and N, a superhard cubic boron nitride was synthesized in 1957 [2], which often replaces diamond in tools because it is more stable at high temperatures, although its hardness is less than that of diamond. As for binary compounds of the C-N system, ultrahard phases of C3N4 stoichiometry have been predicted already in the 1990s [3,4], with allegedly higher bulk moduli and hardness than those of diamond due to the short and strong C-N bonds. However, the existence of such phases has been strongly questioned, and even in the case of experimentally synthesized C-N phases, their stoichiometry and structures are still controversial due to the heterogeneity of the samples obtained. In particular, it should be noted that strong N–N repulsions induced by short N–N nonbonded distances do not favor the synthesis of the nitrogen-rich C-N phases, especially when CVD and PVD methods are used. The hypothesis of gaseous carbon dinitrogen release during vapor phase deposition led to the proposal of model systems such as CN2 [5,6], based on studies within the well-established quantum mechanics framework of density functional theory (DFT) [7,8].
Thus, the search for new C-N phases is mainly limited to the carbon-rich systems (i.e., N/C < 1), where many hypothetical binary compounds have been predicted in recent years, such as rhombohedral and cubic C17N4 [9], cubic C4N [9], monoclinic [10] and orthorhombic [11] C3N, tetragonal and orthorhombic C11N4 [12], cubic α-C3N2 and β-C3N2 [13], cubic C4N3 [14], orthorhombic CN [15], etc. Among the hypothetical nitrogen-rich C-N phases, only CN2 is worth mentioning, whose stability has been predicted over a wide range of pressures [16]. However, none of these many predicted structures have been synthesized, with the sole exception of the orthorhombic CN [17] and diamond-like C2N [18]. The progress and problems in predicting novel ultrahard binary and ternary phases of the B-C-N system from first principles in conjunction with available experimental data have been recently discussed in our review paper [19].
In the present paper, we further develop the problem of carbon nitrides, which remains of great interest, as follows from recent works [20,21]. In the framework of accurate DFT calculations, we have designed a new cohesive and stable (mechanically and dynamically) tetragonal heptacarbon (C7) allotrope, characterized by both tetrahedral (sp3-like) and trigonal (sp2-like) carbon sites, which was used as a template for selective substitutions of carbon for nitrogen, leading to two new carbon-rich nitrides of C5N2 stoichiometry, as well as to a new tetragonal C3N4 (Appendix A).

2. Computational Framework

The identification of the ground state structures corresponding to the energy minima and prediction of their mechanical and dynamical properties were carried out via DFT-based calculations using the Vienna Ab initio Simulation Package (VASP) code [22,23] and the projector augmented wave (PAW) method [23,24] for the atomic potentials. Exchange correlation (XC) effects were considered using the generalized gradient functional approximation (GGA) [25]. Relaxation of the atoms to the ground state structures was performed with the conjugate gradient algorithm according to Press et al. [26]. The Blöchl tetrahedron method [27] with corrections according to the Methfessel and Paxton scheme [28] was used for geometry optimization and energy calculations, respectively. Brillouin-zone (BZ) integrals were approximated using a special k-point sampling according to Monkhorst and Pack [29]. Structural parameters were optimized until atomic forces were below 0.02 eV/Å and all stress components were <0.003 eV/Å3. The calculations were converged at an energy cutoff of 400 eV for the plane-wave basis set in terms of the k-point integration in the reciprocal space from kx(6) × ky(6) × kz(4) up to kx(12) × ky(12) × kz(8) for the final convergence and relaxation to zero strains. In the post-processing of the ground-state electronic structures, the charge density projections were operated on the carbon atomic sites. The charge density was analyzed using Bader’s atoms in molecules (AIM) theory [30] to assess the trends of charge transfer. The elastic constants Cij were calculated to evaluate the mechanical stability and hardness as described in detail below in the Section 3.3. The dynamic stabilities were evaluated from the phonon band structures according to Togo et al. [31]. The electronic band structures were obtained using the all-electron DFT-based ASW method [32] and the GGA XC functional [25]. The VESTA (Visualization for Electronic and Structural Analysis) program [33] was used to visualize the crystal structures and charge densities.

3. Results and Discussion

3.1. Crystal Chemistry and Structure Characterization

The calculated structures at energy minima from fully geometry-optimized setups are detailed in Table 1 and the structures are shown in Figure 1. The body-centered tetragonal ‘neoglitter’ C6 [34] has been used to devise a novel carbon allotrope, tetragonal C7. The ‘neoglitter’ structure, shown in ball-and-stick and polyhedral representations, consists of C4 tetrahedra at the eight corners and body-center, connected by C-C pairs along the vertical direction. C6 is then expressed as (Ctet)2(Ctrg)4 as made up of tetrahedral and trigonal sites, readily available to be used as a template for CN2 carbon dinitride. However, such hypothesis was rejected as preliminary calculations resulted in a weakly cohesive compound (−0.62 eV/atom) compared to the magnitudes of the other dinitrides in Table 1. Furthermore, the CN2 dinitride was characterized by unstable dynamic behavior due to negative phonons.
The novel C7 consists of an elongated tetragonal box with tetrahedral carbon atoms occupying the eight corners as in ‘neoglitter’ and in-between trigonal carbons leading to a loss of the body-centered character with subsequent change of space group I-4m2 (No. 119). The fully optimized geometric structure is shown in Figure 1b, and the polyhedral representation on the right highlights the sp3-sp2 hybrid structure of C7 with tetrahedral and trigonal motifs. The first data column of Table 1 provides the numerical description of the structure inscribed as the following phases in the tetragonal P-4m2 (No. 115) space group. Four different carbon sites build the structure. The C4 tetrahedra have an angle of 106.5°, slightly different from the perfect 109.47° of C(sp3), and the corresponding distance is d(C-C) = 1.52 Å, also smaller than in diamond (1.54 Å). The other shorter distances correspond to the vertical C-C bonds. Subtracting the atomic C contribution of −6.6 eV from the total energy gives an atomic average cohesive energy of −2.09 eV/atom. This value is lower than the cohesive energy of diamond (−2.49 eV/atom).
Substitution of C for N was performed along with two working hypotheses to obtain the stoichiometry of C5N2:
  • The first hypothesis was to create a dinitride with N-N pairs as shown in Figure 1c, where the white spheres represent the N@C (2e) 0,0,z (Table 1, second data column). After geometric optimization, the final structure shows little deviation from pristine C7, as can be seen from the internal parameters, while the distances show a short d(N-N) = 1.36 Å, larger than the length of the N≡N bond (~1.09 Å) in the N2 molecule. The volume is smaller, resulting in a slightly higher density. The atom-averaged cohesive energy is significantly reduced compared to C7, but the carbonitride labeled α-C5N2 remains cohesive with Ecoh/atom = −1.50 eV/at.
  • The second hypothesis was to replace N@C (2g) ½, 0, z (Table 1, third data column), thus creating CN4 tetrahedra at the eight corners as schematically shown in Figure 1d. Therefore, d(C-N) = 1.49 Å is shorter than in C7 and α-C5N2. The remarkable result of such β-C5N2 is the higher (by 0.23 eV) cohesive atomic energy of this configuration, which is clearly more favorable than that of α-C5N2. Finally, the density is higher at 3.252 g/cm3. Therefore, the tetrahedral configuration prevails.
Finally, we found that the simultaneous N@C (2g) and N@C (2e) substitutions resulted in the known C3N4 nitrogen-rich stoichiometry. The structure thus constructed was then calculated for the sake of completeness. The results given in the last column of Table 1 are also presented in Figure 1e, with N shown as red spheres. The structure shows a small d(N-N) and the volume is the smallest, resulting in the highest density. The cohesive energy is smaller than those of both α-C5N2 and β-C5N2. C3N4, proposed by Teter and Hemley [4], calculated for comparison, gives a cohesive energy of −1.29 eV/atom, i.e., it is more cohesive than tet-C3N4 (Table 1).
Simulated X-ray diffraction patterns of three new phases are shown in Figure 2. It is easy to see that they all have the same topology, which is significantly different from diamond topology.

3.2. Trends of Charge Transfer

For Pauling negativities χN = 3.04 and χC = 2.55, charge transfers from C to N of different magnitudes are expected. From the high BZ integration, the charge density output (CHGCAR) was analyzed based on the AIM theory according to Bader [30]. The analyses of the carbon nitrides studied showed a small amount of charge transfer.
In α-C5N2, ΔQC = 0.676 and ΔQN = −1.69; in β-C5N2, ΔQC = 0.588 and ΔQN = −1.47; and in tet-C3N4, ΔQC = 2.22 and ΔQN = −1.67. Regarding C7, despite the different atomic positions and chemical behaviors from sp3/sp2, the charge transfers were found to be negligible, letting the allotrope be considered as a covalent chemical system. Both C5N2 polymorphs can be considered polar covalent, but the large amount of charge transfer in tet-C3N4 allows it to be characterized as ionocovalent, reinforced by the large number of electrons transferred from C to N. Mechanical and dynamic properties of the novel phases were further analyzed.

3.3. Mechanical Properties

3.3.1. Elastic Constants

The mechanical properties were obtained through determining the elastic constants from the strain–stress relationship. For the tetragonal crystal systems, the terms of mechanical stability are defined by elastic constants Cij as follows:
Cii (i = 1, 3, 4, 6) > 0; C11 > C12, C11 + C33 − 2C13 > 0; 2C11 + C33 + 2C12 + 4C13 > 0
The calculated Cij of C7 and new carbon nitrides are given in Table 2. All values are positive, and their combinations given above obey the rules of mechanical stability. The bulk (BV) and shear (GV) moduli were calculated via Voigt averaging the single-crystal elastic constants using ELATE software [35]. As a general trend, the bulk moduli are close to 390 GPa for β-C5N2 and tet-C3N4, proportional to their densities, and close to 360 GPa for C7 and α-C5N2, which have lower densities. The shear moduli of the new phases are close to 230 GPa, except for a slightly lower value of 205 GPa for tet-C3N4. There is a general trend of the Pugh modulus ratio (GV/BV) < 1 with the lowest value of 0.52 for tet-C3N4. This is indicative of ductile systems with low brittleness, probably due to the presence of trigonal atomic fragments in the structures of all phases.

3.3.2. Hardness

To estimate the Vickers hardness (HV) of the proposed new phases, two contemporary theoretical models were used. The thermodynamic (T) model [36], which is based on thermodynamic properties and crystal structure, generally shows perfect agreement with experimental data and is therefore recommended for the hardness evaluation of superhard phases [19]. The Lyakhov–Oganov (LO) model [37] considers the topology of the crystal structure, the strength of covalent bonding, the degree of ionicity, and directionality. As shown earlier [19], empirical models using elastic properties are not reliable in the case of superhard compounds of light elements, so we have not considered these models. Fracture toughness (KIc) was evaluated using the Mazhnik–Oganov model [38]. The results for the currently proposed phases compared to the properties of diamond and hypothetical cubic C3N4 [3] are summarized in Table 3 and Table 4.
The hardness of tetragonal C7 is the highest of all the new phases, but it is ~15% less than that of diamond. This is to be expected since the density of this phase is significantly lower (by a factor of 1.2). The hardness of the three new C-N phases is even lower and is ~65 GPa for all of them, closer to the hardness of cubic C3N4 [3].
The elastic moduli of all four new phases are close, with the Pugh modulus ratio (G/B) being less than 1 and varying from 0.52 for tet-C3N4 to 0.64 for C7.
These trends in mechanical properties are consistent with the fact that nitrogen plays a major role in the structure softening of the three C-N compounds, while the presence of trigonal carbon in all new phases further reinforces this tendency.

3.4. Dynamical Properties with the Phonons

To verify the dynamic stability of the new phases, the phonon band structures were calculated. The methodology consists of calculating the phonon modes by finite displacements of the atoms from the equilibrium position found through the preliminary optimization of the crystal geometry—another heavier method exists as density functional perturbation theory (DFPT). This leads the forces on all atoms constituting the crystal structure to raise off their equilibrium positions. Analysis of the forces associated with a systematic set of displacements yields a series of phonon frequencies. The phonon dispersion curves along the main lines of the tetragonal Brillouin zone are then obtained using the Python code “Phonopy” (see [31] for further details).
Figure 3 shows the phonon band structures for the new carbon allotrope and carbon nitride phases that develop along the main lines of the primitive tetragonal Brillouin zone.
Frequencies (ω) expressed in units of THz are along the vertical direction. In all four panels the frequencies are positive, providing the signature of dynamically stable systems. Between 0 and ~10 THz, the first three bands correspond to lattice rigid translations: two in-plane and one out-of-plane; they are labeled as acoustic modes. At higher energies and up to the highest ω, the remaining bands are the optical modes. Knowing that diamond has the highest band observed using Raman spectroscopy at ω~40 THz [41], all frequency maxima of the four phases in Figure 3a–d are above this magnitude, up to ω~49 THz for α-C5N2 (Figure 3b), assigned to N-N stretching. Then, while diamond has pure sp3 hybridization, the currently studied C7 and new carbon nitrides have mixed sp3/sp2 hybridization, the latter leading to high frequencies for very short distances (Table 1).

3.5. Electronic Band Structures

Figure 4 shows the electronic band structure plots obtained using the all-electron DFT-based augmented spherical method (ASW) [32] and the GGA exchange correlation DFT functional. In all four panels along the vertical direction of the energy, there is a continuous development of bands. Consequently, there is no separation (gap) between the valence band (VB) filled with electrons and the empty conduction band (CB), indicating metallic-like behavior. Then, zero energy is with respect to the Fermi level (EF). Such electronic behavior can be attributed to the delocalization induced by the π electrons of the sp2-like hybridized atoms present in all systems. Insulating behavior requires electron localization as found in diamond, which is known to be an electronic insulator with a wide band gap of ~5 eV between VB and CB. In fact, diamond is characterized only by tetrahedral C(sp3).

3.6. Density of States (DOS)

For a better assessment of the role of each site and chemical constituent in the electronic structure and, further, on the band structures, the density of states (DOS) was calculated and plotted. The energy reference is now along the x-axis. As in the band structures, zero energy is considered with respect to the Fermi level EF, visualized with a straight vertical red line. Along the y-axis, the DOSs are in units of 1/eV. Figure 5a, showing C7′s DOS, exhibits a continuous skyline of the different constituents’ valence states from the lowest VB energy (−22 eV–−15 eV), where s states predominate, up to −5 eV, where p states predominate. Such similar shapes of the site-resolved DOS are indicative of quantum mixing of the respective valence states describing the chemical bonding. From −5 eV up to EF, lower-intensity DOS of C1 is observed with almost zero intensity at EF, while the other C DOSs show large intensity at EF. This illustrates the different behavior of tetrahedral C1 versus trigonal C2.
Upon introduction of nitrogen (C5N2, Figure 5b,c), additional features appear with a broadening of VB by almost 5 eV due to the larger electronegativity of N compared to C. VB is divided into s-like VB (from −27 eV to −20 eV) dominated by N DOS, and the part with p-states (from −20 eV to EF), which also shows dominant N states. Less intense states are found at EF, mainly originating from C2, indicating s-like conductivity. Additionally, in Figure 5d, the C3N4 DOSs show increasingly larger nitrogen contribution throughout VB and above EF in CB.

4. Conclusions

The presence of mixed sp3/sp2 carbon hybridization in a novel tetragonal C7 allotrope allowed us to devise three new carbon nitrides along hypotheses related to the selective occupations of either atomic site. Two C5N2 phases were found to be mechanically and dynamically stable, with the tetrahedral CN4 phase being more cohesive than the phase with linear N-N motifs. Tetragonal C3N4 was similarly devised through allowing N to occupy both sp3 and sp2 carbon sites in C7, resulting in a ductile metallic phase. C7 was found to be cohesive and classified as ductile with a Vickers hardness of 65 GPa. The ductility was found to be enhanced with a decrease of the Pugh modulus ratio G/B along the series from 0.65 to 0.52 with a resultant decrease in hardness. Due to the delocalization of π electrons of the sp2-like hybridized atoms, metallic behavior characterizes all four phases. If synthesized, either in bulk via solid-state synthesis at high pressures and high temperatures, or through sputtering and other deposition techniques, the new phases are likely to improve our understanding of the long-standing problem of nitrogen chemistry in carbon nitrides.

Author Contributions

Conceptualization, S.F.M.; methodology, S.F.M. and V.L.S.; validation, S.F.M. and V.L.S.; formal analysis, S.F.M. and V.L.S.; investigation, S.F.M. and V.L.S.; resources, S.F.M.; data curation, S.F.M. and V.L.S.; writing—original draft preparation, S.F.M.; writing—review and editing, V.L.S.; visualization, S.F.M. and V.L.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available on request.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Crystal structures publications in the Cambridge Crystallographic Data Centre (CCDC):
C7: Deposition No 2255470. Refcode: JIJCUE. DOI: 10.5517/ccdc.csd.cc2fq039
α-C5N2: Deposition No 2252296. Refcode: DIMGOZ. DOI: 10.5517/ccdc.csd.cc2flpqg
β-C5N2: Deposition No 2252297. Refcode: DIMGUF. DOI: 10.5517/ccdc.csd.cc2flprh
tet-C3N4: Deposition No 2271012. Refcode: TINMIQ. DOI: 10.5517/ccdc.csd.cc2g75gc

References

  1. Qin, J.-X.; Yang, X.-G.; Lv, C.-F.; Li, Y.-Z.; Liu, K.-K.; Zang, J.-H.; Yang, X.; Dong, L.; Shan, C.-X. Nanodiamonds: Synthesis, properties, and applications in nanomedicine. Mater. Des. 2021, 210, 110091. [Google Scholar] [CrossRef]
  2. Wentorf, R.H. Cubic form of boron nitride. J. Chem. Phys. 1957, 26, 956. [Google Scholar] [CrossRef]
  3. Liu, A.; Wentzcovitch, R. Stability of carbon nitride solids. Phys. Rev. B 1994, 50, 10362–10365. [Google Scholar] [CrossRef] [PubMed]
  4. Teter, D.; Hemley, R. Low-compressibility carbon nitrides. Science 1996, 271, 53–55. [Google Scholar] [CrossRef]
  5. Weihrich, R.; Eyert, V.; Matar, S.F. Structure and electronic properties of new model dinitride systems: A density-functional study of CN2, SiN2, and GeN2. Chem. Phys. Lett. 2003, 373, 636–641. [Google Scholar] [CrossRef] [Green Version]
  6. Weihrich, R.; Matar, S.F.; Betranhandy, E.; Eyert, V. A model study for the breaking of N2 from CNx within DFT. Solid State Sci. 2003, 5, 701–703. [Google Scholar] [CrossRef]
  7. Hohenberg, P.; Kohn, W. Inhomogeneous electron gas. Phys. Rev. B 1964, 136, 864–871. [Google Scholar] [CrossRef] [Green Version]
  8. Kohn, W.; Sham, L. Self-consistent equations including exchange and correlation effects. Phys. Rev. A 1965, 140, 1133–1138. [Google Scholar] [CrossRef] [Green Version]
  9. Li, Z.; Luo, K.; Liu, B.; Sun, L.; Ying, P.; Liu, C.; Hu, W.; He, J. Superhard carbon-rich C-N compounds hidden in compression of the mixture of carbon black and tetracyanoethylene. Carbon 2021, 184, 846–854. [Google Scholar] [CrossRef]
  10. Li, X.; Xing, M. Novel carbon-rich nitride C3N: A superhard phase in monoclinic symmetry. Comp. Mater. Sci. 2019, 158, 170–177. [Google Scholar] [CrossRef]
  11. Hao, J.; Liu, H.; Lei, W.; Tang, X.; Lu, J.; Liu, D.; Li, Y. Prediction of a superhard carbon-rich C-N compound comparable to diamond. J. Phys. Chem. C 2015, 119, 28614–28619. [Google Scholar] [CrossRef]
  12. Ding, Y. Mechanical properties and hardness of new carbon-rich superhard C11N4 from first-principles investigations. Physica B 2012, 407, 2282–2288. [Google Scholar] [CrossRef]
  13. Tian, F.; Wang, J.; He, Z.; Ma, Y.; Wang, L.; Cui, T.; Chen, C.; Liu, B.; Zou, G. Superhard semiconducting C3N2 compounds predicted via first-principles calculations. Phys. Rev. B 2008, 78, 235431. [Google Scholar] [CrossRef]
  14. Qian, Y.; Wu, H. D-C4N3: A superhard ferromagnetic half-metal predicted by first-principles study. Phys. Lett. A 2022, 423, 127814. [Google Scholar] [CrossRef]
  15. Wang, X. Polymorphic phases of sp3-hybridized superhard CN. J. Chem. Phys. 2012, 138, 184506. [Google Scholar] [CrossRef]
  16. Dong, H.F.; Oganov, A.R.; Zhu, Q.; Qian, G.R. The phase diagram and hardness of carbon nitrides. Sci. Rep. 2015, 5, 9870. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Stavrou, E.; Lobanov, S.; Dong, H.; Oganov, A.R.; Prakapenka, V.B.; Konopkova, Z.; Goncharov, A.F. Synthesis of ultra-incompressible sp3-hybridized carbon nitride with 1:1 stoichiometry. Chem. Mater. 2016, 28, 6925–6933. [Google Scholar] [CrossRef]
  18. Komatsu, T. Shock synthesis and characterization of new diamond-like carbon nitrides. Phys. Chem. Chem. Phys. 2004, 6, 878–880. [Google Scholar] [CrossRef]
  19. Solozhenko, V.L.; Matar, S.F. Prediction of novel ultrahard phases in the B-C-N system from first principles: Progress and problems. Materials 2023, 16, 886. [Google Scholar] [CrossRef]
  20. Pickard, C.J.; Salamat, A.; Bojdys, M.J.; Needs, R.J.; McMillan, P.F. Carbon nitride frameworks and dense crystalline polymorphs. Phys. Rev. B 2016, 94, 094104. [Google Scholar] [CrossRef] [Green Version]
  21. Talukdar, M.; Deb, P. Recent progress in research on multifunctional graphitic carbon nitride: An emerging wonder material beyond catalyst. Carbon 2022, 192, 308–331. [Google Scholar] [CrossRef]
  22. Kresse, G.; Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 1996, 54, 11169. [Google Scholar] [CrossRef] [PubMed]
  23. Kresse, G.; Joubert, J. From ultrasoft pseudopotentials to the projector augmented wave. Phys. Rev. B 1999, 59, 1758–1775. [Google Scholar] [CrossRef]
  24. Blöchl, P.E. Projector augmented wave method. Phys. Rev. B 1994, 50, 17953–17979. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Perdew, J.; Burke, K.; Ernzerhof, M. The Generalized Gradient Approximation made simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef] [Green Version]
  26. Press, W.; Flannery, B.; Teukolsky, S.; Vetterling, W. Numerical Recipes, 2nd ed.; Cambridge University Press: New York, NY, USA, 1986. [Google Scholar]
  27. Blöchl, P.; Jepsen, O.; Anderson, O. Improved tetrahedron method for Brillouin-zone integrations. Phys. Rev. B 1994, 49, 16223–16233. [Google Scholar] [CrossRef]
  28. Methfessel, M.; Paxton, A. High-precision sampling for Brillouin-zone integration in metals. Phys. Rev. B 1989, 40, 3616–3621. [Google Scholar] [CrossRef] [Green Version]
  29. Monkhorst, H.; Pack, J. Special k-points for Brillouin Zone integration. Phys. Rev. B 1976, 13, 5188–5192. [Google Scholar] [CrossRef]
  30. Bader, R.F.W. A bond path:  a universal indicator of bonded interactions. J. Phys. Chem. A 1998, 102, 7314–7323. [Google Scholar] [CrossRef]
  31. Togo, A.; Tanaka, I. First principles phonon calculations in materials science. Scr. Mater. 2015, 108, 1–5. [Google Scholar] [CrossRef] [Green Version]
  32. Eyert, V. Basic notions and applications of the augmented spherical wave method. Int. J. Quantum Chem. 2000, 77, 1007–1031. [Google Scholar] [CrossRef]
  33. Momma, K.; Izumi, F. VESTA3 for three-dimensional visualization of crystal, volumetric and morphology data. J. Appl. Crystallogr. 2011, 44, 1272–1276. [Google Scholar] [CrossRef]
  34. Matar, S.F.; Solozhenko, V.L. Novel ultrahard sp2/sp3 hybrid carbon allotrope from crystal chemistry and first principles: Body centered tetragonal C6 (‘neoglitter’). Diam. Relat. Mater. 2023, 133, 109747. [Google Scholar] [CrossRef]
  35. Gaillac, R.; Pullumbi, P.; Coudert, F.-X. ELATE: An open-source online application for analysis and visualization of elastic tensors. J. Phys. Condens. Matter 2016, 28, 275201. [Google Scholar] [CrossRef] [Green Version]
  36. Mukhanov, V.A.; Kurakevych, O.O.; Solozhenko, V.L. The interrelation between hardness and compressibility of substances and their structure and thermodynamic properties. J. Superhard Mater. 2008, 30, 368–378. [Google Scholar] [CrossRef]
  37. 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] [Green Version]
  38. Mazhnik, E.; Oganov, A.R. A model of hardness and fracture toughness of solids. J. Appl. Phys. 2019, 126, 125109. [Google Scholar] [CrossRef]
  39. Bindzus, N.; Straasø, T.; Wahlberg, N.; Becker, J.; Bjerg, L.; Lock, N.; Dippel, A.-C.; Iversen, B.B. Experimental determination of core electron deformation in diamond. Acta Cryst. A 2014, 70, 39–48. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Brazhkin, V.V.; Solozhenko, V.L. Myths about new ultrahard phases: Why materials that are significantly superior to diamond in elastic moduli and hardness are impossible. J. Appl. Phys. 2019, 125, 130901. [Google Scholar] [CrossRef]
  41. Krishnan, R.S. Raman spectrum of diamond. Nature 1945, 155, 171. [Google Scholar] [CrossRef]
Crystals 13 01111 g001 550
Figure 1. Neoglitter’ C6 [34] (a) and tetragonal C7 (b) with ball-and-stick and polyhedral representations; α-C5N2 with linear N-N fragments (white spheres) (c); tet-C5N2 with nitrogen forming CN4 tetrahedra (yellow spheres) (d); tet-C3N4 with red (regrouping white + yellow) nitrogen spheres (e). Brown spheres represent carbon.
Figure 1. Neoglitter’ C6 [34] (a) and tetragonal C7 (b) with ball-and-stick and polyhedral representations; α-C5N2 with linear N-N fragments (white spheres) (c); tet-C5N2 with nitrogen forming CN4 tetrahedra (yellow spheres) (d); tet-C3N4 with red (regrouping white + yellow) nitrogen spheres (e). Brown spheres represent carbon.
Crystals 13 01111 g001
Crystals 13 01111 g002 550
Figure 2. Simulated X-ray diffraction patterns of new tetragonal phases.
Figure 2. Simulated X-ray diffraction patterns of new tetragonal phases.
Crystals 13 01111 g002
Crystals 13 01111 g003 550
Figure 3. Phonon band structures of tetragonal C7 (a), α-C5N2 (b), tet-C5N2 (c), and tet-C3N4 (d).
Figure 3. Phonon band structures of tetragonal C7 (a), α-C5N2 (b), tet-C5N2 (c), and tet-C3N4 (d).
Crystals 13 01111 g003
Crystals 13 01111 g004a 550Crystals 13 01111 g004b 550
Figure 4. Electronic band structures of tetragonal C7 (a), α-C5N2 (b), tet-C5N2 (c), and tet-C3N4 (d).
Figure 4. Electronic band structures of tetragonal C7 (a), α-C5N2 (b), tet-C5N2 (c), and tet-C3N4 (d).
Crystals 13 01111 g004aCrystals 13 01111 g004b
Crystals 13 01111 g005a 550Crystals 13 01111 g005b 550
Figure 5. Site-projected density of states (DOS) of tetragonal C7 (a), α-C5N2 (b), β-C5N2 (c), and tet-C3N4 (d).
Figure 5. Site-projected density of states (DOS) of tetragonal C7 (a), α-C5N2 (b), β-C5N2 (c), and tet-C3N4 (d).
Crystals 13 01111 g005aCrystals 13 01111 g005b
Table
Table 1. Crystal structure parameters of new tetragonal (P-4m2) phases.
Table 1. Crystal structure parameters of new tetragonal (P-4m2) phases.
C7α-C5N2β-C5N2tet-C3N4
a (Å)2.52902.52622.49292.4848
c (Å)7.37207.14287.22797.1433
C1 (1a) 0, 0, 0✓✓✓✓
C2 (2g) ½, 0, zz = 0.3039z = 0.1158z = 0.3050z = 0.3094
C(N) (2g) ½, 0, zz = 0.1146(C) z = 0.3057(N) z = 0.113(N) z = 0.1138
C’(N’) (2e) 0, 0, zz = 0.3991(N’) z = 0.4043(C) z = 0.3977(N’) z = 0.4053
Cell volume (Å3)47.1545.5844.9244.10
Density (g/cm3)2.9613.2083.2563.466
Bond length (Å)1.396
1.454
1.521		1.356 (C-C)
1.367 (N-N)
1.446 (C-N)
1.510 (C-C)		1.389 (C-N)
1.415 (C-C)
1.476 (C-C)
1.491 (C-N)		1.353 (N-N)
1.412 (C-N)
1.485 (C-N)
Etotal (eV)−60.88−57.13−58.70−54.23
Ecoh/atom (eV)−2.09−1.50−1.73−1.03
Table
Table 2. Elastic constants and bulk (BV) and shear (GV) moduli (in GPa) of new tetragonal phases. G/B is the Pugh modulus ratio.
Table 2. Elastic constants and bulk (BV) and shear (GV) moduli (in GPa) of new tetragonal phases. G/B is the Pugh modulus ratio.
C11C12C13C33C44C66BVGVGV/BV
C7771141221224781383652380.65
α-C5N27382092137543883622200.61
β-C5N283220145126480983942300.58
tet-C3N454215428210271272173942050.52
Table
Table 3. Vickers hardness (HV) and bulk moduli (B0) of carbon allotropes and carbon nitrides calculated in the framework of the thermodynamic model of hardness [36].
Table 3. Vickers hardness (HV) and bulk moduli (B0) of carbon allotropes and carbon nitrides calculated in the framework of the thermodynamic model of hardness [36].
Space Groupa = b (Å)c (Å)ρ (g/cm3)HV (GPa)B0 (GPa)
DiamondFd-3m3.56661 *3.51798445 †
C7 #115P-4m22.52907.37202.96182373
α-C5N2 #115P-4m22.52627.14283.20865383
β-C5N2 #115P-4m22.49297.22793.25666389
tet-C3N4#115P-4m22.48487.14333.46665385
c-C3N4#215 [3]P-43m3.4300 [3]3.78871421
* Ref. [39]; † Ref. [40].
Table
Table 4. Mechanical properties of carbon allotropes and carbon nitrides: Vickers hardness (HV), bulk modulus (B), shear modulus (G), Young’s modulus (E), Poisson’s ratio (ν), and fracture toughness (KIc).
Table 4. Mechanical properties of carbon allotropes and carbon nitrides: Vickers hardness (HV), bulk modulus (B), shear modulus (G), Young’s modulus (E), Poisson’s ratio (ν), and fracture toughness (KIc).
HVBGVE ‡v ‡KIc §
T *LO †B0 *BV
GPa MPa·m½
Diamond9890445 **530 **11380.0746.4
C7 #11582763733652385860.2324.7
α-C5N2 #11565643833622205490.2474.2
β-C5N2 #11566693893942305770.2564.9
tet-C3N4#11565673853942055240.2784.8
c-C3N4#215 [3]7173421422 ††397 ††9070.1427.2
* Thermodynamic model [36]; † Lyakhov–Oganov model [37]; ‡ E and ν values calculated from elastic constants using ELATE software [35]; § Mazhnik–Oganov model [38]; ** Ref. [40]; †† calculated from the literature data on elastic constants [3] using the Voigt approach.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Matar, S.F.; Solozhenko, V.L. Crystal Chemistry and First Principles Studies of Novel Superhard Tetragonal C7, C5N2, and C3N4. Crystals 2023, 13, 1111. https://doi.org/10.3390/cryst13071111

AMA Style

Matar SF, Solozhenko VL. Crystal Chemistry and First Principles Studies of Novel Superhard Tetragonal C7, C5N2, and C3N4. Crystals. 2023; 13(7):1111. https://doi.org/10.3390/cryst13071111

Chicago/Turabian Style

Matar, Samir F., and Vladimir L. Solozhenko. 2023. "Crystal Chemistry and First Principles Studies of Novel Superhard Tetragonal C7, C5N2, and C3N4" Crystals 13, no. 7: 1111. https://doi.org/10.3390/cryst13071111

APA Style

Matar, S. F., & Solozhenko, V. L. (2023). Crystal Chemistry and First Principles Studies of Novel Superhard Tetragonal C7, C5N2, and C3N4. Crystals, 13(7), 1111. https://doi.org/10.3390/cryst13071111

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