PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of materialsLink to Publisher's site
 
Materials (Basel). 2016 September; 9(9): 726.
Published online 2016 August 25. doi:  10.3390/ma9090726
PMCID: PMC5457092

Mechanical and Electronic Properties of XC6 and XC12

Martin O. Steinhauser, Academic Editor

Abstract

A series of carbon-based superconductors XC6 with high Tc were reported recently. In this paper, based on the first-principles calculations, we studied the mechanical properties of these structures, and further explored the XC12 phases, where the X atoms are from elemental hydrogen to calcium, except noble gas atoms. The mechanically- and dynamically-stable structures include HC6, NC6, and SC6 in XC6 phases, and BC12, CC12, PC12, SC12, ClC12, and KC12 in XC12 phases. The doping leads to a weakening in mechanical properties and an increase in the elastic anisotropy. C6 has the lowest elastic anisotropy, and the anisotropy increases with the atomic number of doping atoms for both XC6 and XC12. Furthermore, the acoustic velocities, Debye temperatures, and the electronic properties are also studied.

Keywords: first-principles calculations, carbides, mechanical properties, elastic anisotropy

1. Introduction

Elemental carbon exhibits a rich diversity of structures and properties, due to its flexible bond hybridization. A large number of stable or metastable phases of the pure carbon, including the most commonly known, graphite and diamond, and other various carbon allotropes [1,2,3,4] (such as lonsdaleite, fullerene, and graphene, etc.), and diversified carbides [5,6,7,8,9,10,11], have been studied in experiments and theoretical calculations. Graphite, which is the most stable phase at low pressure, has a sp2-hybridized framework and is ultrasoft semimetallic, whereas diamond, stable at high pressure, is superhard, insulating with a sp3 network. Recently, a novel one-dimensional metastable allotrope of carbon with a finite length was first synthesized by Pan et al. [1], called Carbyne. It has a sp-hybridized network and shows a strong purple-blue fluorescence. The successful synthesis of Carbyne is a great promotion for the further analysis on properties and applications. The 2D material MXenes as a promising electrode material, which is early transition metal carbides and carbon nitrides, is reported [11], owing to its metallic conductivity and hydrophilic nature. These properties of different carbides are appealing. To find superhard superconductors, researches designed some carbide superconductors, such as boron carbides and XC6 structure with cubic symmetry. The diamond-like BxCy system, which is superhard and superconductive, has also attracted much interest [5,6,7,8,9,10]. The best simulated structure of the synthesized d-BC3 (Pmma-b phase) has a Vickers hardness of 64.8 GPa, showing a superhard nature, and its Tc reaches 4.9–8.8 K [5]. The P-4m2 polymorph of d-BC7 with a low energy also has a high Vickers hardness of 75.2 GPa [8]. Furthermore, Wang et al. [9] explored more potential superhard structures of boron carbide, uncovering the stability is mainly contributed by the elemental boron at low pressure, and by the carbon at high pressure. The novel metastable carbon structure C6 bcc is predicted with a cubic symmetry [12]. It is an indirect band gap semiconductor with 2.5 eV, calculated by the local density approximation. Recently, doped with simple metals, Lu et al. [13] studied a series of sodalite-based carbon structures, similar to the boron-doped diamond. Although they found these structures are all metastable, some of these structures show a superconductivity, e.g., the critical temperature of NaC6 is 116 K. In this paper, we mainly study the mechanical properties of these eleven XC6 phases (HC6, LiC6, NC6, OC6, FC6, NaC6, AlC6, SiC6, PC6, SC6, and ClC6) which is of dynamical stability and, for comparison, C6 is also calculated. In addition, the XC12 structures are systematically explored, in which the X atom is from H to Ca, except He, Ne, and Ar. The doping-induced changes in elastic constant, modulus, the anisotropy of elasticity and acoustic velocity, Debye temperature, and the electronic structures are also studied.

2. Results and Discussion

As shown in Figure 1a, the structure of XC6 is obtained by doping the X atom into the C6 bcc structure at (0, 0, 0). It is of Im-3m symmetry (No. 229), consisting of two formula units (f.u.) per unit cell. Each C atom has four nearest neighbors with the bond angle of 90° or 120°. The XC6 structure has four C4 rings and eight C6 rings. In Table 1, the calculated lattice parameter a of C6 has a good agreement with the available result [12], and is smaller than that of the XC6 structures. By removing the corner atoms and only leaving the center X atom, the XC12 structure is obtained (Figure 1b). All of the XC12 phases are smaller than the corresponding XC6 phases, but larger than the C6 phase in the lattice parameter.

Figure 1
Unit cell of XC6 (a) and XC12 (b). The black and blue spheres represent C and X atoms, respectively.
Table 1
Calculated lattice parameter a, elastic constants Cij (GPa), mechanical stability, bulk modulus B (GPa), shear modulus G (GPa), Young’s modulus E (GPa), Poisson’s ratio ν, and B/G ratio.

The formation enthalpies of XC6 in [13] and XC12 structures are calculated reference to diamond and the most stable X phase at ambient pressure. The equations are given by ΔHXC6=(HXC6HX6HC)/7, and ΔHXC12=(HXC12HX12HC)/13, and the calculated results are shown in Figure 2. The positive values indicate these phases are metastable. The two curves of the formation enthalpy follow a similar trend, where the F-doped carbides have the lowest ΔH, and the PC6 and CC12 have the largest ΔH in XC6 and XC12, respectively. Compared to other doped elements of the second and the third periods in the XC6 and XC12, fluorine (F) possesses the largest electronegativity difference relative to C, leading to a stronger interaction between F and C atoms; thus, FC6 and FC12 phases are more stable.

Figure 2
Formation enthalpy of XC6 and XC12.

The calculated elastic constants and moduli are listed in Table 1. The generalized Born’s mechanical stability criteria of cubic phase are given by [15]: C11>0,C44>0, C11>|C12|, and (C11+2C12)>0. In Table 1, the C6 and HC6, NC6, and SC6 have the mechanical stability, and they are also dynamically stable [13]. The XC12 has ten mechanically stable phases, but only six of these phases have the dynamical stability (BC12, CC12, PC12, SC12, ClC12, and KC12) due to the absence of the imaginary frequency in the whole Brillouin zone (see Figure 3 and Figure 4). The S is the only element that is capable to make not only XC6, but also XC12, stable.

Figure 3
Phonon spectra of dynamically stable phases (a) BC12; (b) CC12; (c) PC12; (d) SC12; (e) ClC12; and (f) KC12.
Figure 4
Phonon spectra of dynamically unstable phases (a) LiC12; (b) BeC12; (c) MgC12; and (d) AlC12.

By Voigt-Reuss-Hill approximations [16,17,18], the bulk modulus B and shear modulus G can be obtained, and the Young’s modulus E and Poisson’s ratio ν are defined as [19,20] E=9BG/(3B+G) and ν=(3B2G)/[2(3B+G)]. HC6 has the largest bulk modulus of 346 GPa, showing the best ability to resist the compression. The shear modulus is often used to qualitatively predict the hardness, and Young’s modulus E is defined as the ratio between stress and strain to measure the stiffness of a solid material. In Table 1, C6 is the largest in shear modulus and Young’s modulus, which means that doping leads to a weakening in mechanical properties. The Poisson’s ratio exhibits the plasticity; usually, the larger the value, the better the plasticity. According to Pugh [21], C6, HC6, BC12, CC12, and PC12 are brittle materials (B/G < 1.75), while NC6, SC6, SC12, ClC12, and KC12 are ductile materials (B/G > 1.75). This conforms the calculated results of Poisson’s ratio.

The elastic anisotropy is important for the analysis on the mechanical property and, thus, the universal elastic anisotropy index (AU), Zener anisotropy index (A), and the percentage anisotropy in compressibility and shear are calculated. For the cubic phase, the universal elastic anisotropy index [22] is defined as: AU=5GV/GR+BV/BR6, the nonzero value suggests an anisotropy characteristic. Furthermore, it is known that C44 represents the resistance to deformation with respect to a shear stress applied across the (100) plane in the [010] direction, and (C11C12)/2 represents the resistance to shear deformation by a shear stress applied across the (110) plane in the [11¯0] direction. For an isotropic crystal, the two shear resistances turn to identical. Therefore, Zener [23] introduced A=2C44/(C11C12) to quantify the extension of anisotropy. The value of 1.0 represents the isotropy, and any deviation from 1.0 indicates the degree of the shear anisotropy. The percentage anisotropy in compressibility and shear are given by: AB=(BVBR)/(BV+BR) and AG=(GVGR)/(GV+GR) [24]. The AB is always 0.0 for a cubic phase. As shown in Table 2, C6 has the lowest anisotropy. The universal elastic anisotropy index and the percentage anisotropy in shear is increasing with the atomic number of doped element for both XC6 and XC12, and the anisotropy which obtains from the shear anisotropic factor is also increasing, except SC6 and KC12. Furthermore, owing to the percentage anisotropy in shear of C6, BC12, and CC12 being slight, they are almost isotropic.

Table 2
Universal elastic anisotropy index (AU), Zener anisotropy index (A), and percentage anisotropy in shear (AG).

The elastic anisotropies are calculated with the elastics anisotropy measures (ElAM) code [25,26] which makes the representations of non-isotropic materials easy and visual. For the cubic phase, the representation in xy, xz, and yz planes are identical, as a result, only the xy plane is presented. The 2D figures of the differences in each direction of Poisson’s ratio are shown in Figure 5. The maximum value curves and minimum positive value curves of C6 and XC6 stable phases are illustrated in Figure 5a,b, and those of XC12 stable phases are shown in Figure 5c,d. Particularly, the SC12 and ClC12 have the negative minimum Poisson’s ratio. It is seen that all of the structures are anisotropic and C6 has the lowest anisotropy, suggesting the doping increase the elastic anisotropy. The largest value of maximum curve is in the same direction of the lowest value of minimum positive value curve for each structure. Furthermore, for XC12 phases, the anisotropy of Poisson’s ratio is increasing with the atomic number. The negative minimum Poisson’s ratio of SC12 and ClC12 indicate these two phases have auxeticity [27], and ClC12 is more prominent than SC12.

Figure 5
2D representations of Poisson’s ratio. (a) Maximum of C6 and XC6 stable phases; (b) minimum positive of C6 and XC6 stable phases; (c) maximum of XC12 stable phases; and (d) minimum positive and minimum negative of XC12 stable phases; particularly, ...

The directional dependence of the Young’s modulus [28] are demonstrated in Figure 6 and Figure 7. The distance from the origin of system of coordinate to the surface equals the Young’s modulus in this direction, and thus any departure from the sphere indicates the anisotropy. As shown, all of the phases are anisotropic, and the anisotropy of Young’s modulus is increasing with the doping atomic number. For the S-doped phases, which have stable XC6 and XC12 structures, the maximum (minimum) values of SC6 and SC12 are 650 (291) and 371 (175) GPa, respectively. The Emax/Emin ratio of SC6 (2.23) is slightly larger than that of SC12 (2.12), indicating the SC6 is more anisotropic.

Figure 6
Directional dependence of the Young’s modulus of C6 (a); HC6 (b); NC6 (c); and SC6 (d).
Figure 7
Directional dependence of the Young’s modulus of BC12 (a); CC12 (b); PC12 (c); SC12 (d); ClC12 (e); and KC12 (f).

The acoustic velocity is a fundamental parameter to measure the chemical bonding characteristics, and it is determined by the symmetry of the crystal and propagation direction. Brugger [29] provided an efficient procedure to calculate the phase velocities of pure transverse and longitudinal modes from the single crystal elastic constants. The cubic structure only has three directions [001], [110], and [111] for the pure transverse and longitudinal modes and other directions are for the qusi-transverse and qusi-longitudinal waves. The acoustic velocities of a cubic phase in the principal directions are [30]:

  • for [100], vl=C11/ρ, [010]vt1=[001]vt2=C44/ρ,
  • for [110], vl=(C11+C12+2C44)/2ρ, [11¯0]vt1=(C11C12)/2ρ, [001]vt2=C44/ρ,
  • for [111], vl=(C11+2C12+4C44)/3ρ, [112¯]vt1=vt2=(C11C12+C44)/3ρ.

where ρ is the density of the structure, vl is the longitudinal acoustic velocity, and vt1 and vt2 refer the first transverse mode and the second transverse mode, respectively. It should be noted that there is a misprint for equation of [11¯0]vt1 in [30]. Here, the correct expression is given. Based on the elastic constants, the anisotropic properties of acoustic velocities indicate the elastic anisotropy in these crystals. As a fundamental physical parameter which correlates with many physical properties of solids, the Debye temperature can be obtained from the average acoustic velocity: ΘD=hkB[3n4π(NAρM)]1/3vm, where h and kB are the Planck and Boltzmann constants, respectively; NA is Avogadro’s number; n is the total number of atoms in the formula unit; M is the mean molecular weight, and ρ is the density. The average acoustic velocity is vm=[(2/vtm3+1/vlm3)/3]1/3, where vlm=(B+4G/3)/ρ is the average longitudinal acoustic velocity, and vtm=G/ρ is the average transverse acoustic velocity.

All of the calculated acoustic velocities and Debye temperatures of diamond and stable XC6 and XC12 phases are shown in Table 3. Diamond is larger than C6 and doped structures in anisotropic and average acoustic velocity. The densities are increasing and the average acoustic velocities are decreasing with the atomic number, except NC6, which has a much smaller shear modulus. Compared to C6, the doping results in a decrease in the average acoustic velocity and Debye temperature. For the element S, which makes both XC6 and XC12 phases stable, the average acoustic velocity of SC6 decreases by 38.65% than C6, and that of SC12 by 35.96%. Furthermore, it can be found that the Debye temperature is decreasing with the atomic number, except SC6. The ΘD characterizes the strength of the covalent bond in solids, so the strength of the covalent bond is lower for the phase which has the larger atomic number of doping atom.

Table 3
Density (g/cm3), anisotropic acoustic velocities (m/s) and average acoustic velocity (m/s).

Figure 8 shows the electronic band structure and density of state (DOS) of XC12 stable phases. The dash line represents the Fermi level (EF). The electronic properties of XC6 have been studied in [13]. For XC12, all of the band structures cross the Fermi level in the Brillouin zone, showing the metallic nature. The conduction band and valence band are mainly characterized by the contributions of C-p states, whereas the DOS near the Fermi level originated from the p orbital electrons of the doped element, except the ClC12 and KC12.

Figure 8Figure 8Figure 8
Electronic band structure and density of state of BC12 (a); CC12 (b); PC12 (c); SC12 (d); ClC12 (e); and KC12 (f).

3. Computational Methods

The calculations are performed with the first-principles calculations. The structural optimizations are using the density functional theory (DFT) [31,32] with the generalized gradient approximation (GGA), which is parameterized by Perdew, Burke, and Ernzerrof (PBE) [33]. The Broyden-Fletcher-Goldfarb-Shanno (BFGS) minimization scheme [34] was used in the geometry optimization, and the total energy convergence tests are within 1 meV/atom. When the total energy is 5.0 × 10−6 eV/atom, the maximum ionic Hellmann-Feynman force is 0.01 eV/Å, the maximum stress is 0.02 GPa and the maximum ionic displacement is 5.0 × 10−4 Å, the structural relaxation will stop. The energy cutoff is 400 eV, and the K-points separation is 0.02 Å−1 in the Brillouin zone.

4. Conclusions

By using the first-principles calculations, the analyses on the mechanical properties of XC6 and the further exploration of XC12 structures are given. The formation enthalpies of dynamically stable XC6 phases and all of the XC12 structures, and the elastic constants, are calculated. There are ten structures which have the mechanical and dynamical stability (C6, HC6, NC6, SC6, BC12, CC12, PC12, SC12, ClC12, and KC12). The elastic modulus and anisotropy of the ten structures are studied and, in these structures, C6 has the lowest elastic anisotropy and the anisotropy increases with the atomic number. The doping leads to the weakening in mechanical properties and the increase in the elastic anisotropy. In addition, Debye temperatures and the anisotropy of acoustic velocities are also studied. The electronic properties studies show the metallic characteristic for XC6 and XC12 phases.

Acknowledgments

This work was financially supported by the Natural Science Foundation of China (No. 11204007), Natural Science Basic Research plan in Shaanxi Province of China (grant No.: 2016JM1026, 20161016), and Education Committee Natural Science Foundation in Shaanxi Province of China (grant No.: 16JK1049).

Author Contributions

Author Contributions

Qun Wei and Meiguang Zhang designed the project; Quan Zhang and Qun Wei performed the calculations, Qun Wei and Quan Zhang prepared the manuscript, Meiguang Zhang revised the paper, all authors discussed the results and commented on the manuscript.

Conflicts of Interest

Conflicts of Interest

The authors declare no conflict of interest.

References

1. Pan B., Xiao J., Li J., Liu P., Wang C., Yang G. Carbyne with finite length: The one-dimensional sp carbon. Sci. Adv. 2015;1 doi: 10.1126/sciadv.1500857. [PMC free article] [PubMed] [Cross Ref]
2. Boulfelfel S.E., Oganov A.R., Leoni S. Understanding the nature of “superhard graphite” Sci. Rep. 2012;2:471. doi: 10.1038/srep00471. [PMC free article] [PubMed] [Cross Ref]
3. Wei Q., Zhang M., Yan H., Lin Z., Zhu X. Structural, electronic and mechanical properties of Imma-carbon. EPL. 2014;107:27007. doi: 10.1209/0295-5075/107/27007. [Cross Ref]
4. Falcao E.H.L., Wudl F. Carbon allotropes: Beyond graphite and diamond. J. Chem. Technol. Biotechnol. 2007;82:524–531. doi: 10.1002/jctb.1693. [Cross Ref]
5. Liu H., Li Q., Zhu L., Ma Y. Superhard and superconductive polymorphs of diamond-like BC3. Phys. Lett. A. 2011;375:771–774. doi: 10.1016/j.physleta.2010.12.034. [Cross Ref]
6. Li Q., Wang H., Tian Y., Xia Y., Cui T., He J., Ma Y., Zou G. Superhard and superconducting structures of BC5. J. Appl. Phys. 2010;108:023507. doi: 10.1063/1.3452374. [Cross Ref]
7. Xu L., Zhao Z., Wang Q., Wang L.M., Xu B., He J., Tian Y. Prediction of a superconductive superhard material: Diamond-like BC7. J. Appl. Phys. 2011;110:013501. doi: 10.1063/1.3601349. [Cross Ref]
8. Liu H., Li Q., Zhu L., Ma Y. Superhard polymorphs of diamond-like BC7. Solid State Commun. 2011;151:716–719. doi: 10.1016/j.ssc.2011.02.013. [Cross Ref]
9. Wang D.Y., Yan Q., Wang B., Wang Y.X., Yang J., Yang G. Predicted boron-carbide compounds: A first-principles study. J. Chem. Phys. 2014;140:224704. doi: 10.1063/1.4882071. [PubMed] [Cross Ref]
10. Zheng B., Zhang M., Chang S. Mechanical properties and atomic explanation of plastic deformation for diamond-like BC2. Materials. 2016;9:514. doi: 10.3390/ma9070514. [Cross Ref]
11. Naguib M., Mochalin V.N., Barsoum M.W., Gogotsi Y. 25th anniversary article: MXenes: A new family of two-dimensional materials. Adv. Mater. 2014;26:992–1005. [PubMed]
12. Ribeiro F.J., Tangney P., Louie S.G., Cohen M.L. Hypothetical hard structures of carbon with cubic symmetry. Phys. Rev. B. 2006;74 doi: 10.1103/PhysRevB.74.172101. [Cross Ref]
13. Lu S., Liu H., Naumov I.I., Meng S., Li Y., Tse J.S., Yang B., Hemley R.J. Superconductivity in dense carbon-based materials. Phys. Rev. B. 2016;93 doi: 10.1103/PhysRevB.93.104509. [Cross Ref]
14. Fan Q., Wei Q., Chai C., Yan H., Zhang M., Lin Z., Zhang Z., Zhang J., Zhang D. Structural, mechanical, and electronic properties of P3m1-BCN. J. Phys. Chem. Solids. 2015;79:89–96. doi: 10.1016/j.jpcs.2014.12.008. [Cross Ref]
15. Wu Z., Zhao E., Xiang H., Hao X., Liu X., Meng J. Crystal structures and elastic properties of superhard IrN2 and IrN3 from first principles. Phys. Rev. B. 2007;76:054115. doi: 10.1103/PhysRevB.76.054115. [Cross Ref]
16. Voigt W. Lehrburch der Kristallphysik. Teubner; Leipzig, Germany: 1928.
17. Reuss A. Berechnung der Fließgrenze von Mischkristallen auf Grund der Plastizitätsbedingung für Einkristalle. J. Appl. Math. Mech. 1929;9:49–58. doi: 10.1002/zamm.19290090104. [Cross Ref]
18. Hill R. The Elastic behaviour of a crystalline aggregate. Phys. Soc. Lond. Sect. A. 1952;65:349–354. doi: 10.1088/0370-1298/65/5/307. [Cross Ref]
19. Wei Q., Zhang M., Guo L., Yan H., Zhu X., Lin Z., Guo P. Ab initio studies of novel carbon nitride phase C2N2(CH2) Chem. Phys. 2013;415:36–43. doi: 10.1016/j.chemphys.2013.02.015. [Cross Ref]
20. Wei Q., Zhang Q., Zhang M. Crystal structures and mechanical properties of Ca2C at high pressure. Materials. 2016;9:570. doi: 10.3390/ma9070570. [Cross Ref]
21. Pugh S.F. Relations between the elastic moduli and the plastic properties of polycrystalline pure metals. Lond. Edinb. Dublin Philos. Mag. J. Sci. 1954;45:823–843. doi: 10.1080/14786440808520496. [Cross Ref]
22. Ranganathan S.I., Ostoja-Starzewski M. Universal elastic anisotropy index. Phys. Rev. Lett. 2008;101:055504. doi: 10.1103/PhysRevLett.101.055504. [PubMed] [Cross Ref]
23. Zener C. Elasticity and Anelasticity of Metals. University of Chicago; Chicago, IL, USA: 1948.
24. Ravindran P., Fast L., Korzhavyi P.A., Johansson B., Wills J., Eriksson O. Density functional theory for calculation of elastic properties of orthorhombic crystals: Application to TiSi2. J. Appl. Phys. 1998;84:4891–4904. doi: 10.1063/1.368733. [Cross Ref]
25. Marmiera 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. doi: 10.1016/j.cpc.2010.08.033. [Cross Ref]
26. Zhang Q., Wei Q., Yan H., Fan Q., Zhu X., Zhang J., Zhang D. Mechanical and electronic properties of P42/mnm silicon carbides. Z. Naturforsch. A. 2016;71:387–396. doi: 10.1515/zna-2015-0539. [Cross Ref]
27. Baughman R.H., Shacklette J.M., Zakhidov A.A., Stafström S. Negative Poisson’s ratios as a common feature of cubic metals. Nature. 1998;392:362–365. doi: 10.1038/32842. [Cross Ref]
28. Wei Q., Zhang M., Yan H., Li R., Zhu X., Lin Z., Yao R. A new superhard phase of C3N2 polymorphs. Z. Naturforsch. A. 2015;70:1001. doi: 10.1515/zna-2015-0270. [Cross Ref]
29. Brugger K. Determination of third-order elastic coefficients in crystals. J. Appl. Phys. 1965;36:768. doi: 10.1063/1.1714216. [Cross Ref]
30. Sun L., Gao Y., Xiao B., Li Y., Wang G. Anisotropic elastic and thermal properties of titanium borides by first-principles calculations. J. Alloys Compd. 2013;579:457–467. doi: 10.1016/j.jallcom.2013.06.119. [Cross Ref]
31. Hohenberg P., Kohn W. Inhomogeneous electron gas. Phys. Rev. 1964;136:B864–B871. doi: 10.1103/PhysRev.136.B864. [Cross Ref]
32. Kohn W., Sham L.J. Self-consistent equations including exchange and correlation effects. Phys. Rev. 1965;140:A1133–A1138. doi: 10.1103/PhysRev.140.A1133. [Cross Ref]
33. Perdew J.P., Burke K., Ernzerhof M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996;77:3865–3868. doi: 10.1103/PhysRevLett.77.3865. [PubMed] [Cross Ref]
34. 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. doi: 10.1006/jcph.1996.5612. [Cross Ref]

Articles from Materials are provided here courtesy of Multidisciplinary Digital Publishing Institute (MDPI)