Characterizations of B and N heteroatoms as substitutional doping on structure,stability, and aromaticity of novel heterofullerenes evolved from the smallest fullerene cage C20: a density functional theory perspective

The structural stabilities and electronic properties of C₂₀ fullerene and some its incorporated boron and nitrogen derivatives are probed at B3LYP/AUG‐cc‐pVTZ level of theory. According to density functional theory results, the topology of inserted B or N heteroatoms in [20]‐fullerene perturbs strongly the stability, energy, geometry, charge, polarity, nucleus‐independent chemical shifts, aromaticity, and highest‐occupied molecular orbital and lowest‐unoccupied molecular orbital (HOMO–LUMO) gap of the resulting heterofullerenes. Vibrational frequency (υ_min) calculations show that except N₁₀C₁₀, all other B_bN_nC_{20‐(b + n)} heterofullerenes with b, and n = 0, 4, 5, 8, and 10 are true minima. The calculated band gaps (?E_HOMO–LUMO) of B₈C₁₂, and N₈C₁₂ (2.86 eV), show them the most stable heterofullerenes against electronic excitations. While 10 B substituting in equatorial position increase the conductivity of B₁₀C₁₀ through decreasing its band gaps, 10 N doping in equatorial position enhance stability of N₁₀C₁₀ against electronic excitations via increasing its band gaps. High natural bond orbital and Mulliken charge transfer on the surfaces of B atoms, especially B₅N₅C₁₀with five B–N bonds in the equatorial position, provokes further investigation on its possible application for hydrogen storage. Nucleus‐independent chemical shift values show that B₅N₅C₁₀ is the most aromatic species. The calculated heat of atomization per carbon (ΔH_at/C) of B₈C₁₂ shows it the most thermodynamic stable heterofullerenes of each. Copyright © 2016 John Wiley & Sons, Ltd.     

Adsorption of NH3 and NO2 molecules on C48B6N6 heterofullerene: A DFT study on electronic properties

Adsorption of NH₃ and NO₂ molecules on the external surface of C₄₈B₆N₆ heterofullerene is investigated using DFT method. Attachment of NH₃ and NO₂ on C₄₈B₆N₆ heterofullerenes are compared with the bare C₄₈B₆N₆ model optimized at the B3LYP/6-31G^? level. The high surface binding energies indicates that ammonia undergoes chemical adsorption and could be compatible with the long recovery time but C₄₈B₆N₆ should be good NO₂ sensors with quick response as well as short recovery time. Total (TDOS) and partial (PDOS) density of state calculations is also considered to elucidate the difference in the NH₃ and NO₂ gas detection mechanism of C₄₈B₆N₆. The overlap population density of state (OPDOS) indicated that the chemical adsorption is due to the overlap of atomic orbitals below the Fermi level. The calculated results suggest that the C₄₈B₆N₆ heterofullerene is a suitable sensor material for NO₂ and is an ideal material for elimination and filtering of ammonia.     

A comparative study on the Ge6C14 heterofullerene nanocages: a density functional survey

Ten C₁₄Ge₆ heterofullerene isomers of C₂₀ have been investigated by density functional theory (DFT) methods with Becke 3‐Parameter (Exchange), Lee, Yang and Parr (B3LYP) functional at the 6‐311 + G*, 6‐311++G** and AUG‐cc‐pVTZ levels. In contrast to identical bonds in the latter, contractions of C═C double bonds are encountered at the expense of longer C―Ge bonds in the former. Vibrational frequency analysis confirms that all of the nanocages are true minima. In contrast to the common belief, for obtaining highly silicon‐doped stable heterofullerenes, that the silicon dopants must be completely isolated from each other by means of strong C═C double bonds. Here, linking the germanium substitutions together is an applicable strategy for obtaining highly doped stable isolated heterofullerenes since it avoids weak heteroatom─heteroatom bonds. Therefore, none of the computed heterofullerenes collapses to open, to deform, or to segregate fullerenic cages. As to band gaps (ΔE_HOMO‐LUMO), and nucleus‐independent chemical shifts at cage centers (NICS (0)), C₁₄Ge₆‐2 immerges with the highest value. Hence, it is predicted to be the most stable against electronic excitation. It contains 2 Ge─Ge single bonds at the cap‐equatorial positions. On the other hand, as to zero‐point vibrational energy and heat of atomization (ΔH_at), C₁₄Ge₆‐8 appears with the lowest and highest value, respectively. It contains 6 alternating germanium atoms in the equatorial and cap positions. Thus, it is predicted to be the most thermodynamically stable. So, germanium substitution leads to a high charge distribution on the surfaces of all the isomers specially C₁₄Ge₆‐9 with +1.496 charged germanum atoms. C₁₄Ge₆ isomers seem to be a good candidate for the hydrogen storage material.     

Amantadine antiparkinsonian drug adsorption on the AlN and BN nanoclusters: A computational study

To find a sensor for Amantadine (AM) antiparkinsonian drug, we studied its interaction with Al₁₂N₁₂ and B₁₂N₁₂ nanoclusters by density functional theory calculations. The AM molecule attaches via its –NH₂ group to the Al or B atoms of Al₁₂N₁₂ or B₁₂N₁₂ with Gibbs free energy change about −31.5 or −26.1 kcal/mol. Increasing the AM concentration, the interaction becomes weaker due to steric effects. The AM adsorbs on the Al₁₂N₁₂ and B₁₂N₁₂ with two different mechanisms, including electrostatic and charge transfer, respectively. The AM significantly reduces the Al₁₂N₁₂ work function from 4.50 to 3.66 eV, increasing the electron field emission. Thus, the AlN cluster may be a work function type sensor. Upon the AM adsorption on the BN cage, the HOMO level is largely destabilized, reducing the E_g from 6.84 to 5.01 eV which largely increases the electrical conductivity. This indicates that the BN cluster may be a potential electronic sensor.     

The Hartree-Fock exchange effect on the CO adsorption by the boron nitride nanocage

We studied the effect of Hartree-Fock (HF) exchange percentage of a density functional on the adsorption properties and electronic sensitivity of the B₁₂N₁₂ nanocluster to CO molecule. It was found that by an increase in the %HF, the LUMO level is nearly constant while the HOMO level is strongly stabilized, expanding the HOMO-LUMO gap (E_g). Also, the volume of the all structures decreased and the sensitivity of the B₁₂N₁₂ is slightly increased to CO molecule. For the pristine B₁₂N₁₂ cluster, the B66 and B64 bonds are about 1.43 and 1.49 Å at 10% HF, and 1.23 and 1.26 Å at 100% HF, respectively. The HF exchange between 10–20% may predict an accurate E_g for the B₁₂N₁₂ system. We concluded that functionals with a large %HF such as M06-HF, and M06-2X may significantly overestimate the E_g, and bond strength. We obtained a parabolic relationship between the %HF and the adsorption energy of CO molecule on the B₁₂N₁₂ cluster. Also, an increase in the %HF predicts a larger charge transfer from the CO molecule to the cage.     

Computer simulations and analysis of structural and energetic features of crystalline cage energetic compound: 2, 4, 6, 8, 12‐pentanitro‐10‐(3, 5, 6‐trinitro (2‐pyridyl))‐2, 4, 6, 8, 12‐hexaazatetracyclo [5.5.0.03,11.05,9]dodecane

First principles molecular orbital and plane‐wave ab initio calculations have been used to investigate the structural and energetic properties of a new cage compound 2, 4, 6, 8, 12‐pentanitro‐10‐(3, 5, 6‐trinitro (2‐pyridyl))‐2, 4, 6, 8, 12‐hexaazatetracyclo [5.5.0.0^3,11.0^5,9]dodecane (PNTNPHATCD) in both the gas and solid phases. The molecular orbital calculations using the density functional theory methods at the B3LYP/6‐31G(d,p) level indicate that both the heat of formation and strain energy of PNTNPHATCD are larger than those of 2, 4, 6, 8, 10, 12‐hexanitro‐2, 4, 6, 8, 10, 12‐hexaazatetracyclo [5.5.0.0.0] dodecane (CL‐20). The infrared spectra and the thermodynamic property in gas phase were predicted and discussed. The calculated detonation characteristics of PNTNPHATCD estimated using the Kamlet–Jacobs equation equally matched with those of CL‐20. Bond‐breaking results on the basis of natural bond orbital analysis imply that C–C bond in cage skeleton, C–N bond in pyridine, and N–NO₂ bond in the side chain of cage may be the trigger bonds in the pyrolysis. The structural properties of PNTNPHATCD crystal have been studied by a plane‐wave density functional theory method in the framework of the generalized gradient approximation. The crystal packing predicted using the Condensed‐phase Optimized Molecular Potentials for Atomistic Simulation Studies (COMPASS) force fields belongs to the Pbca space group, with the lattice parameters a = 20.87 Å, b = 24.95 Å, c = 7.48 Å, and Z = 8, respectively. The results of the band gap and density of state suggest that the N–NO₂ bond in PNTNPHATCD may be the initial breaking bond in the pyrolysis step. As the temperature increases, the heat capacity, enthalpy, and entropy of PNTNPHATCD crystal all increase, whereas the free energy decreases. Considering that the cage compound has the better detonation performances and stability, it may be a superior high energy density compound. Copyright © 2015 John Wiley & Sons, Ltd.     

Electronic structure of crystal-forming fulborenes B<Subscript><Emphasis Type="Italic">n</Emphasis></Subscript>N<Subscript><Emphasis Type="Italic">n</Emphasis></Subscript>

A general approach is formulated to the design of crystal-forming fullerene-like clusters X _n Y _n from which zeolite-like covalent crystals based on IV-IV, III-V, and II-VI binary semiconductor compounds with diamond-like sp ³ bonds can be constructed and synthesized by means of copolymerization through faces. A number of the smallest sized crystal-forming boron nitride clusters are constructed, such as the B₁₂N₁₂, B₁₆N1₆, B₁₈N₁₈, B₂₄N₂₄, B₃₆N₃₆, and B ₆₀N₆₀ fulborenes. The optimized configurations, electronic structures, charge transfers, band gaps, total energies, cohesive energies, and electron density maps of the clusters are calculated using the spin-restricted Hartree-Fock method in the 6–31G basis set. Comparative calculations of the B₆₀N₆₀ fulborene with the use of the density functional theory method have demonstrated that the spin-restricted Hartree-Fock method in the 6–31G basis set is optimum from the standpoint of the accuracy and efficiency.     

Structural, elastic, and electronic properties of icosahedral boron subcarbides (B12C3, B13C2), subnitride B12N2, and suboxide B12O2 from data of SCC-DFTB calculations

The structural, elastic, and electronic properties of a series of icosahedral phases, such as boron subcarbides B₁₂C₃ and B₁₃C₂, subnitride B₁₂N₂, and suboxide B₁₂O₂, have been studied in the framework of the SCC-DFTB method. It has been found that the B₁₂C₂ and B₁₃C₂ phases manifest metal-like properties, while B₁₂C₃ and B₁₂O₂ are semiconductors. The estimates have shown that the insertion of 2p atoms (C, N, or O) into intericosahedral pores of elemental boron can cause both a decrease in its elastic modulus (an increase in the compressibility of B₁₂N₂) and a sharp increase in the modulus B (in subcarbides B₁₂C₃ and B₁₂BCC). On the other hand, the insertion of 2p atoms into α-B₁₂ will favor an increase in its hardness (suboxide B₁₂O₂ will have a maximum hardness).     

Theoretical study on C100 fullerenes and C96X4 (X=N,P, B,Si)

The geometrical structures and electronic properties of six fullerene isomers of C₁₀₀ were studied at the HF/6-31G^? and B3LYP/6-31G^? levels, respectively. The results of the fully optimized calculations show that three C₁₀₀ isomers 449:D₂, 425:C₁ and 442:C₂ are near isoenergetic isomers. The energies and properties of C₁₀₀ hexaanions were calculated. The C₁₀₀^6? (450:D₅) isomer is predicted to be the most stable isomer at the B3LYP/6-31G^? level, and the C₁₀₀^6? (449:D₂) isomer is 44.1 kcal/mol higher in energy. The heterofullerenes C₉₆X₄ (X=N, P, B, Si) formed from the initial C₁₀₀ (449:D₂) have also been investigated at the B3LYP/6-31G^? level. The HOMO–LUMO gaps and aromaticities show that the replacement of fullerene carbon atoms with four heteroatoms can enhance the electronic stabilization of C₁₀₀ (449:D₂).     

The electronic and structural properties of BN and BP nano-cages interacting with OCN: A DFT study

The adsorption of OCN⁻ (cyanato anion) on boron nitride (B₁₂N₁₂ and B₁₆N₁₆) and boron phosphide nano-cages (B₁₂P₁₂ and B₁₆P₁₆) in terms of energetic, geometric, and electronic properties are studied using density functional theory calculations. Our study results indicated that the first OCN⁻ strongly prefers to be adsorbed from its N atom upon B atoms of the nano-cages than the O atoms of OCN⁻. These findings have been rationalized using frontier molecular orbitals and total electron density plots. The energy gap of the B₁₂P₁₂ is significantly reduced upon the adsorption of OCN⁻ compared to B₁₂N₁₂, thus leading to the increase in electrical conductance of nano-cage.     

On the validity of formal electron counting rules in lithium silicides

The validity of widely used electron counting rules (i.e. (8-N) rule) is investigated in two of the four stable crystalline lithium silicides (Li₁₂Si₇ and Li₂₁Si₅) by means of semiempirical MO calculations of the INDO-type on the stage of localized orbitals which are derived from the canonical one-particle set by a unitary transformation. This representation allows for a transparent interpretation of the nature of the chemical bond in the employed systems. Finite cluster models are used in both cases. The system Li₁₂Si₇ is approximated by [Li₂₁Si₄]⁹⁺ units, the latter one by [Li₂₂Si₄]⁴⁺ fragments. It is demonstrated that the extension of the octet configuration is only pretended in these systems. The close cage of Li atoms allows for the formation of an additional bonding one-electron state formed by the 2s AO's at the metal centers. This cage orbital is highly independent from the Si subspace. It can be occupied by two electrons that are not considered in formal counting rules. The geometric condition supporting the formation of such a function is a close cage of Li atoms. This solid-state effect has a molecular counter part in hyperlithiated molecules (i.e. [CLi₄]²⁻, [CLi₅]⁻, CLi₆). It is shown that other crystalline lithium silicides are stabilized by Li 2s cage orbitals too.     

Structural and optical properties of the M@C59X cages (X=N,B and M=Li,Na)

Using B3LYP/6-31G* density functional level of theory, the structural and optical properties of the C₆₀ and M@C₅₉X cages have been investigated. Results indicate that the charge on C atoms and band gap of C₆₀ cage are changed dramatically with the substitution of one B or N atom at one of the C sites and the Li and Na atom encapsulations in the C₆₀ cage. The Mulliken analyses show that the charge is transferred completely between the alkali atoms and the C₅₉X cage. The substitutional and encapsulation doping (SED) reduce the optical gaps of the C₆₀ cage. Also, the oscillator strengths of the absorption peaks are dependent on dopant types.     

Piperidine‐appended imidazolium ionic liquids as task‐specific catalysts: computational study,synthesis, and multinuclear NMR

Imidazolium ionic liquids (IMILs) with a piperidine moiety appended via variable length methylene spacers (with n = 1–4) were studied computationally to assess their potential to act as internal base for N‐heterocyclic carbene (NHC) generation. Proton transfer energies computed by B3LYP/6‐311+G(2d,p) were least endothermic for the basic‐IL with n = 3, whose optimized structure showed the shortest C₂‐H‐‐‐‐N(piperidine) distance. Inclusion of counter anion (Cl or NTf₂) caused dramatic conformational changes to enable close contact between the acidic C₂‐H and the anions. To examine the prospect for internal C₂‐H‐‐‐‐N coordination, multinuclear NMR data (¹H, ¹⁵N, and ¹³C) were computed by gauge independent atomic orbitals–density functional theory (GIAO‐DFT) in the gas phase and in several solvents by the PCM method for comparison with the experimental NMR data for the basic ILs (with n = 2–4) synthesized in the laboratory. These studies indicate that interactions with solvent and counter ion are dominant forces that could disrupt internal C₂‐H‐‐‐‐N coordination/proton transfer, making carbene generation from these basic‐ILs unlikely without an added external base. Therefore, the piperidine‐appended IMILs appear suitable for application as dual solvent/base in organic/organometallic transformations that require the use of mild base, without the necessity to alkylate at C‐2 to prevent N‐heterocyclic carbene formation. Copyright © 2016 John Wiley & Sons, Ltd.     

Substituent effect on structure,stability, and aromaticity of novel BnNmC20–(n+m) heterofullerenes

We are focusing our calculations on the structural stabilities and electronic properties of 26 novel B_nN_mC_20–(n+m) heterofullerenes, with n, m = 1 ? 5, at B3LYP/6‐311++G** and B3LYP/AUG‐cc‐pVTZ levels of theory. Vibrational frequency calculations on C₂₀ and its analogues show that except B₂N₂C₁₆ (1) and B₂N₂C₁₆ (2), all other heterofullerenes are true minima. The heats of atomization energies, binding energy, band gaps (ΔE_HOMO‐LUMO), aromaticity, nucleus‐independent chemical shifts, thermodynamic stability, kinetic stability against electronic excitation, binding energy as a stability criterion of different configurations, geometrical parameters, conformational structures, conductivity, charge transfer, and possibility for hydrogen storage of these heterofullerenes strongly depend on their number of heteroatoms, topology, filling patterns, and locations as well as “B‐site and N‐site attachments.” B₅N₅C₁₀ contains 5 alternating boron and nitrogen atoms in the equatorial position. It is predicted to be thermodynamically and kinetically the most stable against electron excites. Thus, it is energetically favorable and its electronic properties as well as stabilities make it perhaps a good candidate for an experimental investigation and testing verification.     

Photo-induced site-specific nitridation of plasma-deposited B10C2Hx films: A new pathway toward post-deposition doping of semiconducting boron carbides

We show that dopant impurities can be introduced in a controlled, site-specific manner into pre-deposited semiconducting boron carbide films. B―N bond formation has been characterized by X-ray photoelectron spectroscopy for semiconducting B₁₀C₂H_x films exposed to vacuum ultraviolet photons in the presence of NH₃. Core level photoemission data indicate that B―NH₂ bonds are formed at B sites bonded to other boron atoms (B―B), and not at boron atoms adjacent to carbon atoms (B―C) or at carbon atom sites. Nitridation obeys diffusion-limited kinetics. These results indicate that dopant species can be introduced in a controlled, site-specific manner into pre-deposited boron carbide films, as opposed to currently required dopant incorporation during the deposition process.     

Ammonia adsorption on the C30B15N15 heterofullerene: DFT study of nuclear magnetic shielding and electric field gradient tensors of N and B nuclei

Ammonia adsorption on the external surface of C₃₀B₁₅N₁₅ heterofullerene was studied using density functional calculations. Three models of the ammonia-attached C₃₀B₁₅N₁₅ together with the perfect model were optimized at the B3LYP/6-31G^? level. The optimization process reveals that dramatic influences occurred for the geometrical structure of C₃₀B₁₅N₁₅ after ammonia adsorption; the B atom relaxes outwardly and consequently the heterofullerene distorts from the spherical form in the adsorption sites. The chemical shielding (CS) tensors and nuclear quadrupole coupling constants of B and N nuclei were calculated at the B3LYP/6-311G^** level. Our calculations reveal that the B atom is chemically bonded to NH₃ molecule. The B atom in the NH₃-attached form has the largest chemical shielding isotropic (CSI) value among the other boron nuclei. The C_Q parameters of B nuclei at the interaction sites are significantly decreased after ammonia adsorption.     

First-principles calculations of the structural, electronic and magnetic properties of BnN20−n (n = 6−18) clusters

The structures of B _n N_20 ? n    (n = 6?18), the clusters of boron nitride, are investigated by the density functional theory calculations. The structures of the obtained low-lying isomers can be described by the following six prototypes: single ring, double ring, three-ring, graphitic-like sheet, fullerene and others. B₁₀N₁₀ is demonstrated to be the most stable cluster against the nonstoichiometric ones. Nonzero magnetic moments, 1.999, 1.998, 2.000, 3.999 and 1.999μ _B respectively, are found in five B _n N_20?n (n = 6, 7, 11, 12, 13) clusters. Further analysis indicates that the magnetic moment of the B₆N₁₄ cluster is mainly originated from the N atoms, while those of others are from the B atoms. The magnetic moment are finally attributed to the interesting issues of the 2p electrons due to the breaking of local symmetries, the change of coordination number, charge distribution and orbital hybridization.     

Electronic structure,physico-chemical,linear and non linear optical properties analysis of coronene, 6B-, 6N-, 3B3N- substituted C24H12 using RHF,B3LYP and wB97XD methods

In this work we have investigated the effects of substituting carbon atoms with B, N and BN on the electronic structure, physico-chemical, linear and non linear optical properties of Coronene (C24H12) using HF and DFT methods. We have calculated total electronic energy E₀, zero point vibrational energy ZPVE, the enthalpy H, entropy S, molar heat capacity at constant volume C_v, ionization potential IP, electron affinity EA, hardness \(\kappa\), softness \(\vartheta\), electronegativity EN, dipole moment µ, average polarizability \(< \alpha >\), anisotropy \(\Delta \alpha\), the first molecular hyperpolarizability β_mol, second order hyperpolarizability \(\gamma_{av}\), HOMO–LUMO Energy gap E_gap, work function E_F, refractive index n, susceptibility χ, dielectric constant ε and molar refractivity MR of coronene (C24H12), the 6B-, 6N- and 3B3N- substitute-doped C24H12 C18B6H12 C18N6H12 and C18B3N3H12. The E_gap values of the molecules are between 0.91 and 2.36 eV. We observed that β_mol changes slightly when C24H12 is doped with either 6B or 6N even though their β_mol values are too small. However, by doping C24H12 with both 3B and 3N, creating a strong donor–acceptor system, a very large increase in µ and β_mol was found for C18B3N3H12. This study was done using RHF, B3LYP and wB97XD methods with the cc-pVDZ basis set. The studies have shown that doping decreases some of the above properties significantly while some increases significantly compared to pure coronene, suggesting that 6B-, 6N-, and 3B3N-doped Coronene as serious candidates for electronics, optoelectronics and photonic devices.     

Characterization of BCN films synthesized by radiofrequency plasma enhanced chemical vapor deposition

Boron carbonitride (BCN) films have been synthesized on Si(1 0 0) substrate by radio frequency plasma enhanced chemical vapor deposition using tris-(dimethylamino)borane (TDMAB) as a precursor. The deposition was performed at the different RF powers of 400-800 W, at the working pressure of 2×10⁻¹ Torr. The formation of the sp²-bonded BCN phase was confirmed by Fourier transform infrared spectroscopy. X-ray photoelectron spectroscopy measurements showed that B atoms were bonded to C and N atoms to form the BCN atomic hybrid configurations with the chemical compositions of B₅₂C₁₂N₃₆ (sample 1; prepared at the RF power of 400 W), B₅₂C₁₀N₃₈ (sample 2; at 500 W) and B₄₆C₁₈N₃₆ (sample 3; at 800 W), respectively. Near-edge X-ray absorption fine structure (NEXAFS) measurements indicated that B atoms were bonded not only to N atoms but also to C atoms to form various configurations of sp²-BCN atomic hybrids. The polarization dependence of NEXAFS suggested that the predominant hybrid configuration of sp²-BCN films oriented in the direction perpendicular to the Si substrate.