Imdad Hussainab,
Adnan Ali Khanac,
Iftikhar Ahmad*ab,
Rashid Ahmadac and
Saif Ullahab
aDepartment of Physics, University of Malakand, Chakdara, Pakistan. E-mail: ahma5532@gmail.com
bCenter for Computational Materials Science, University of Malakand, Chakdara, Pakistan
cDepartment of Chemistry, University of Malakand, Chakdara, Pakistan
First published on 18th March 2024
In this article we report novel composite materials of bucky ball (C60 fullerene) and III-nitrides (BN, AlN, GaN, InN). The experimental feasibility of the novel composite materials is confirmed through negative binding energies and molecular dynamics simulations performed at 500 K. The structural properties of the novel composites are explored through density functional theory. An unusual phenomenon of surface bowing is observed in the 2D structure of the III-nitride monolayers due to the C60 sticking. This surface bowing systematically increases as one proceeds from BN → AlN → GaN → InN. The electron density difference and Hirshfeld charge density analysis show smaller charge transfer during the complexation, which is probably due to weak van der Waal's forces. The presence of van der Waal's forces is also confirmed by the Atom in Molecule analysis, Reduced Density Gradient Technique and Non-covalent Interaction analysis. This work provides a foundation for further theoretical and experimental studies of the novel phenomenon of systematic bowing in the 2D structure of III-nitride monolayers.
This reduction in size is possible by either reducing the size of a material in one direction (2D materials), or two directions (nanotubes etc.) or three directions (nanodots). The reduced dimensions of materials, with enhanced properties caused by quantum confinement effects, have introduced a new class of nanomaterials. Nanomaterials have attracted much attention in both science and technology for their potential high-tech applications in modern technologies.
The wide range of direct band gap nature of III nitrides (BN, AlN, GaN and InN) make them extensively applicable in the optoelectronic industry. The band gap of III nitrides can also be engineered for a specific device application using various techniques like doping, pressure, strain etc. Rare-earth doped III nitrides are commonly used as phosphor materials in optical displays.1,2
Since the discovery of graphene in 2005,3 two dimensional (2D) monolayers have attracted much attention for their unusual properties, due to the quantum confinement effect caused by the reduced dimensionality. The search for suitable non-toxic, environment-friendly, cheaper and long-lasting novel materials for various technological applications is the growing need of the day. Two dimensional (2D) monolayers, 2D heterostructures4–6 and 2D Janus materials are recently at the forefront of this search for their potential applications.5–10 Among 2D materials, 2D III-nitrides (BN, AlN, GaN, InN) are wide band gap, sp2 hybridized semiconductors. The 2D III-nitrides were synthesized by different experimental techniques like Chemical-Vapor Deposition (CVD), Physical-Vapor Transport (PVT), Metal–Organic Vapor Phase Epitaxy (MOVPE), Molecular-Beam Epitaxy (MBE) and Metal–Organic Chemical Vapor Deposition (MOCVD) methods.11–13 The 2D III-nitrides semiconductors have numerous applications like field emitter, biosensor, ideal layer for fabricating flexible devices, as a substrate in graphene electronic devices, sensors, anode material for light emitting diode, field effect transistors, piezoelectric materials, solar cells and high-speed electronics.
Fullerene (C60) is a carbon allotrope having closed cage like structure, discovered in 1985 by Kroto et al. has opened a new era of nanoscience.14 It consists of sp2 hybridized carbon atoms and is the only soluble material among all elemental carbon allotropes.15 The C60 is also called Buckminster fullerene is the mostly studied fullerene due to its unique 0D cage like structure, charge acceptor character (holding up to six electrons), high temperature, nanometer dimension and large polarizability.16 The properties of fullerene can be tuned by interaction with metals and semiconductors. The unique cage like closed 0D structure of fullerenes enables it to be used as important building blocks for the construction of supramolecular assemblies and micro/nano hybrid materials, having numerous technological applications like fast phototransistor, hydrogen storage, photocatalyst, thermoelectric devices and UV-visible photodetectors.17–21 Fullerene can also be regarded as an electron deficient-polyalkene, so, it is chemically reactive.22 Hence, to extend the functionalities of Fullerene, derivation is feasible. Many efforts have been made to adsorb fullerene on graphene and its derivative for to modify its properties.23–25 The adsorption of C20, C34, C42 and C60 fullerenes on graphene was studied and a sizeable band gap opening was observed.26 The adsorption of C20 and C60 fullerenes on the nitrogenated-holy graphene (NHG) sheet for investigating band gap engineering was performed by DFT method. A decrease in the band gap energy of NHG sheet by fullerene adsorption was observed. Furthermore, the calculations revealed that C20 adsorption on NHG sheet to p-type semiconductor and C60 fullerene adsorption changes NHG sheet to n-type semiconductor.
In this paper we report novel composite materials of C60 fullerene and 2D-III nitrides (h-BN, h-AlN, h-GaN, h-InN monolayer). The experimental feasibility of the novel materials is confirmed by using the binding energies and molecular dynamics simulations at 500 K. Density functional theory calculations are used to explore the structural properties of these compounds. The fascinating, aspect of these studies was the bowing effect in the III nitride monolayers. The bowing effect was increased as we proceeded from BN to InN. The Atom in Molecule, Electron Density Difference, Reduced Density Gradient and Non-covalent Interaction analysis are performed to understand the bonding nature (forces) among the complexes and stability of the complexes. This study is a baseline for further experimental and theoretical studies on the surface bowing effect and physical properties for the potential applications of the novel composite materials.
Eb = (EC60–XN) − (EC60 + EX–N) | (1) |
The optimized geometries of the composites (C60@h-BN, C60@h-AlN, C60@h-GaN and C60@h-InN) material is illustrated in Fig. 1, while the optimized structural parameters like lattice parameters, bond lengths are collected in Table 2. The composite materials are designed as the C60 fullerene is placed over the top of the III-nitrides monolayers and fully geometry relaxation is performed. Initially the fullerene molecule is kept onto the top of monolayers by keeping it at about 2.00 Å distance. After geometry relaxation it is found that the fullerene system is physiosorbed over the surface of nitrides monolayers by interacted through van der Waal's forces having an average bond distance up to 3.00 Å between the C atoms of C60 fullerene and X–N (X = B, Al, Ga, In) atoms of monolayers. The bond lengths X–N increases from 0.03 to 0.07 Å shown by Δ in Table 2 due to the fullerene adsorptions on these monolayers. Moreover, the geometry relaxation shows some interesting results about the planar structure of III-nitrides. As we moved from C60@h-BN composite to C60@h-BN, C60@GaN and AlN@C60 composite the structure distortion enhanced. The planar geometry of monolayers (except h-BN monolayer) changes from planar to bowl shape. This variation is more pronounced in the order C60@h-InN > C60@h-GaN > C60@h-AlN > C60@h-BN (Fig. 1). This result shows that the C60@h-BN composite is more stable compared to other studied composites. This distortion may be due to the enhanced repulsive forces between the pi-conjugated electrons of fullerene system and lone-pair electrons of N atoms of III-nitrides monolayer. This can be explained as there is opposite forces operational between these system like the positive B atoms of III-nitrides and pi-electron of fullerene atoms and the electropositive C atoms of fullerene and lone-pair of N atoms of III-nitrides. Thus, there is a competition between the attractive and repulsive forces between the two materials. Finally, the repulsive forces dominant over the attractive forces and both systems repels each other and they are stabilized at certain distance and formed stable composite material. The sequence of distortion is in the order C60@h-InN > C60@h-GaN > C60@h-AlN > C60@h-BN because lattice parameter i.e. X–N (X = B, Al, Ga, In) distance is larger as we moved down the group. This resulted dominant vibrations of X–N bond moving from B–N to In–N system. Thus, the rigidity in h-BN system is more compared to other III-nitrides which has the ability to these vibrations without any structural distortion and evince the most stable composite structure.
Fig. 1 The optimized geometry of C60@h-BN (a), C60@h-AlN (b), C60@h-GaN (c) and C60@h-InN (d) composite. |
Composite | X–C | X–N | C–N | a | Δ | Eb | V |
---|---|---|---|---|---|---|---|
C60@BN | 2.9 | 1.48 | 3.07 | 12.78 | 0.03 | −0.94 | 2903.19 |
C60@AlN | 3.17 | 1.82 | 3.18 | 15.64 | 0.02 | −1.8 | 4241.9 |
C60@GaN | 3.24 | 1.87 | 3.29 | 15.94 | 0.07 | −2.17 | 4445.8 |
C60@InN | 3.0 | 2.11 | 2.92 | 17.75 | 0.06 | −2.20 | 5713.49 |
The binding energy for all these composites (calculated by eqn (1)) are negative (−0.94 to −2.20 eV), confirm the thermodynamic stability and feasibility of experimental synthesis. The binding energy increases from −0.45 eV which represent the chemisorption state but actually the composites are in physiosorbed state. This increasement in the binding energies is due to the dominant distortion in these systems. This order of binding energy is C60@h-InN (−2.20 eV) > C60@h-GaN (−2.17 eV) > C60@h-AlN (−1.80 eV) > C60@h-BN (−0.94 eV) as shown in Table 2.
The molecular dynamic (MD) simulations of these composites were performed to check the thermal stability of a material at 500 K for 2.0 ps equilibrium at a time period of 1 fs. The potential energy plots and structures variation after MD simulation are displayed in Fig. 2. Due to out of planes vibrations modes at 500 K temperature, there are dangling bonds in these structures and a small variations in potential energies specially for C60–GaN and C60–InN which are consistent with the reported result.37 It is clear from the Fig. 2 the potential energy varies nearby a specific magnitude and no breakdown of any bond in any material is observed so all these materials are thermodynamically stable at 500 K, which are consistent with the energetic analysis.
Fig. 2 Geometries of C60–XN (X = B, Al, Ga, In) composites after MD simulations and corresponding potential energy fluctuation plots. |
The charge transfer was determined by Hirshfeld charge transfer analysis. Here it is clear from Table 3 in case of C60–BN the charge transfer take place from fullerene to BN, same is the case of C60–AlN. In these composites the fullerene is a charge donor molecule and behaves as n-type dopant. While in case of C60–GaN and C60–InN charge is transferred from monolayers sheets to fullerene, so fullerene is charge acceptor molecule and behave as p-type dopants. The charge transfer between the 2D III-nitrides monolayers and fullerene molecule is small. The large bond distances and small charge transfer identify the weak interactions between 2D III-nitrides sheets and fullerene molecule. Hence fullerene is weakly physiosorbed on these monolayers.
Materials | ρ | ∇2ρ | H | QC | QX | QN | Gravimetric density W (%) |
---|---|---|---|---|---|---|---|
a QC, QX, QN represent charge on carbon atom, X atom and N atom. | |||||||
C60–BN | 0.001 | 0.029 | 0.0011 | 0.0088 | 0.1901 | −0.2000 | 54 |
C60–AlN | 0.01 | 0.027 | 0.008 | 0.01 | 0.4243 | −0.4346 | 41 |
C60–GaN | 0.012 | 0.026 | 0.005 | −0.001 | 0.3950 | −0.3940 | 26 |
C60–lnN | 0.012 | 0.025 | 0.004 | −0.016 | 0.4500 | −0.4364 | 18 |
We used Bader quantum atom in molecule theory of atoms in molecules38 to understand the bond strength and nature formed between the fullerene and III-nitrides in composite. For studying the nature of bonding in molecules we need to exploit charge density or electron density (ρ). The AIM analysis of the fullerene binding onto 2D III-nitrides is shown in Fig. 4. In characteristics of bonding, the Laplacian (∇2ρ) plays an important role. In case if ∇2ρ < 0 clarify charge concentrations towards interaction line. This results in lowering potential energy creating attractive forces. On the other hand, ∇2ρ > 0 indicates the repulsive forces act on the nuclei. The electron energy density (H) is another parameter also helpful in identifying the nature of bonding. According to Bader and Essen the interactions are covalent if values of ∇2ρ and H are negative and electron density at BCP > 10−1 a.u while values of electron density less than 10−1 a.u and positive ∇2ρ and H attribute to Vander Waal interactions. The calculated values of electron density (ρ), Laplacian (∇2ρ) and total electron energy density (H) are shown in Table 3. The calculated values of electron density (ρ) for the C60–XN (X = B, Al, Ga and In) are 0.001 a.u, 0.01 a.u, 0.012 a.u and 0.012 a.u are less than 10−1 a.u and the values of Laplacian and total electron energy density for these materials are all positive. The calculated values of electron density for these composites are less than 10−1 a.u and positive ∇2ρ and H attribute to van der Waal's interactions. The calculated bond lengths X–N of these materials shown in Table 1 are in the range of 1.00 to 2.00 Å. The bond distance i.e. X–C and C–N are nearly in the range of 3.00 Å. It is also clear from our calculations of bond length and AIM analysis that the bond strength decreases as we move from C60–BN to C60–InN. The same situations appear in AIMD calculations performed at 500 K, where the wiggles in the structure increase from h-BN to h-InN composites.
Fig. 4 Atom in molecule analysis plots of (a) BN-C60, (b) AlN-C60, (c) GaN-C60, and (d) InN-C60 composites. The number represents the bond critical points. |
The maximum storage capacity of fullerene over the 2D III-nitrides is estimated from gravimetric density39 which was calculated from the following equation
(2) |
According to Table 3 the storage capacity of fullerenes are 54, 41, 26 and 18 for BN, AlN, GaN and InN. Next, the reduced density gradient (RDG) and non-covalent interaction (NCI) analysis are investigated to further understand the nature of interactions and structural distortion in the studied composite. The 2D-RDG and 3D-NCI plots are illustrated in Fig. 5, the red, green and dark blue spikes in the RDG plot represent the repulsive, van der Waal's and strong electrostatic (covalent or ionic) forces in the complex system. The RDG clearly shows that there are green spikes appear at x-axis at the range from −0.005 au to −0.02 au, confirm the presence of van der Waal's forces between the C60 fullerene and III-nitrides materials. This is justified by NCI analysis, the appearance of 3D green spot between the C60 fullerene and III-nitrides evidenced the presence of van der Waal's forces. While the appearance of red spots over the 3D green spot in the composite materials indicates the existence of repulsive forces. The red-spots over green spot is more dominant by moving from C60@h-BN composite to C60@h-InN. While this red spot teardown in case of C60@h-GaN and C60@h-InN composite which justify the more prominent repulsive forces between these systems, which caused larger distortion and instability. These results are in good agreement with the geometrical and AIM analysis.
Fig. 5 Reduced density gradient and non-covalent interaction plots for the studied composite materials. |
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4ra01310g |
This journal is © The Royal Society of Chemistry 2024 |