Influence of NO2 attachment on the nuclear magnetic shielding tensors of N and B nuclei in C30B15N15 heterofullerene: a DFT study

Adsorption of nitrogen dioxide in three different configurations on the exterior surface of C₃₀B₁₅N₁₅ is studied using density functional theory calculations. To this end, we optimized the structures of raw C₃₀B₁₅N₁₅ and nine NO₂–C₃₀B₁₅N₁₅ complexes at the B3LYP/6-31G^* level of theory and then calculated chemical shielding (CS) tensors at the GIAO-B3LYP/6-311G^** level for the optimized structures. The calculated chemical shielding isotropy (CSI), chemical shielding anisotropy (CSA), and orientation of CS tensors (Euler angles) reveal that the adsorption configurations (nitro, trans-nitrite, and cis-nitrite) have different effects on the electronic structure of C₃₀B₁₅N₁₅. Natural atomic charges based on natural population analysis (NPA) were used to justify the changes in CSI values after gas sorption.     

Theoretical study of interactions between 1-alkyl-3-methyimidazolium tetrafluoroborate and dibenzothiophene: DFT,NBO, and AIM analysis

Density functional theory is employed to study the interaction energies between dibenzothiophene (DBT) and 1-alkyl-3-methylimidazolium tetrafluoroborate ([C _n mim]⁺[BF₄]^?). The structures of DBT, 1-ethyl-3-methylimidazolium tetrafluoroborate ([C₂mim]⁺[BF₄]^?), 1-butyl-3-methylimidazolium tetrafluoroborate ([C4mim]+[BF₄]?), 1-hexyl-3-methylimidazolium tetrafluoroborate ([C₆mim]⁺[BF₄]^?), 1-octyl-3-methylimidazolium tetrafluoroborate ([C₈mim]⁺[BF₄]^?), [C₂mim]⁺[BF₄]^?–DBT, [C₄mim]⁺[BF₄]^?–DBT, [C₆mim]⁺[BF₄]^?–DBT and [C₈mim]⁺[BF₄]^?–DBT systems are optimized systematically at the B3LYP/6-31G(d,p) level, and the most stable geometries are obtained by NBO and AIM analyses. The results indicate that DBT and imidazolium rings of ionic liquids are parallel to each other. It is found that the [BF₄]^? anion prefers to be located close to a C1–H9 proton ring in the vicinity of the imidazolium ring and the most stable gas-phase structure of [C _n mim]⁺[BF₄]^? has four hydrogen bonds between [C _n mim]+ and [BF₄]?. There are hydrogen bonding interactions, π–π and C–H–π interactions between [C₈mim]⁺[BF₄]^? and DBT, which is confirmed by NBO and AIM analyses. The calculated interaction energies for the studied ionic liquids can be used to interpret a better extracting ability of [C₈mim]⁺[BF₄]^? to remove DBT, due to stronger interactions between [C₈mim]⁺[BF₄]^? and DBT, in agreement with the experimental results of dibenzothiophene extraction by [C _n mim]⁺[BF₄]^?.     

Investigation of inclusion complex of metformin into selective cyclic peptides as novel drug delivery system: Structure,electronic properties,AIM, and NBO study via DFT

In this study, the density functional theory computational method is used to investigate the encapsulation process of metformin into three types of the cyclic peptides composed of eight serine (CP1), eight glycine (CP2), and four serine‐glycine (CP3) cyclic peptides as a new model in the process of drug delivery in the gas phase. The obtained results using the B3LYP/6‐31++G (d,p) method indicate that the complexes formed are energetically favored. Furthermore, results reveal that the drug encapsulation process is typically chemisorption. The natural bonding orbital analysis shows that the intermolecular interaction of the C2 complex (metformin/CP2) is stronger than the C1 (Metformin/CP1) and C3 (Metformin/CP3) complexes due to greater total charge transfer energy, and the C1 complex is found to be the most favored complex. The theory of atoms in molecule (AIM) method is used to analyze the nature of interactions in different molecular systems. The results show the investigated cyclic peptides as effective carriers of metformin in the nanomedicine field.     

Prediction of heterofullerene stabilities: a combined DFT and chemometric study of C56Pt2, C57Pt2 and C81Pt2

A systematic search of the regioisomers of the heterofullerenes, C57Pt2 and C56Pt2, has been carried out by means of density functional calculations to find the most stable structures. Both heterofullerenes incorporate two metal atoms into the fullerene surface. In the case of C57Pt2, one platinum atom substitutes one carbon atom of C60 and the other platinum atom replaces a C--C bond, whereas in C56Pt2 each platinum atom replaces one C--C bond. Several geometric factors were studied, three of which have particularly important effects on the relative stabilities of the regioisomers: the Pt--Pt separation, the number of C--C bonds remaining after substitution, and the type of C--C bond that is substituted. All these factors indicate that the deformation of the carbon framework is a general factor that governs the relative stabilities of the regioisomers. Because a high number of factors affect the stability of the heterofullerenes we also used chemometric techniques in this study. Partial least-squares (PLS) regression was used to establish the structure-energy relationships of C57Pt2 and C56Pt2 heterofullerenes. The understanding gained of the factors that affect the relative isomers stabilities has allowed us to predict the stabilities of larger disubstituted carbon cages, for example, C81Pt2 heterofullerene.     

Intramolecular hydrogen bonds in the sulfonamide derivatives of oxamide, dithiooxamide, and biuret. FT-IR and DFT study, AIM and NBO analysis

The hydrogen bonding in [(1-arylsulfonylamino-2,2,2-trichloro)ethyl]biuret 1, [(1-arylsulfonylamino-2,2,2-trichloro)ethyl]oxamide 2, and [(1-arylsulfonylamino-2,2,2-trichloro)ethyl]dithiooxamide 3, the sulfonamide derivatives of biuret 4, oxamide 5, and dithiooxamide 6, has been studied by molecular spectroscopy and DFT theoretical calculations including frequency calculations, at the B3LYP/6-311+G(d,p) level of theory. The analysis of the CO?HN and CS?HN intramolecular hydrogen bonds closing the five- and six-membered rings employing the atoms-in-molecules (AIM) method using the MP2(full)/6-311++G(d,p) wave functions has shown that their stability is increased in comparison to the original molecules and is much higher in the thiocarbonyl compounds. The results of the AIM and the NBO analysis of donor-acceptor interactions are in good agreement with each other and with the experimental FT-IR spectroscopy data.     

Molecular and electronic structure of several heterofullerene BNC58 and B2N2C56 oligomers and [B2N2C56] n macromolecule

Molecular and electronic structure of heterofullerene BNC₅₈ (C _s) and B₂N₂C₅₆ (C _2h) monomers, B₂N₂C₁₁₆ and B₄N₄C₁₁₂ dimers, and B₆N₆C₁₆₈ trimer (the last three molecules withC _2h symmetry) was simulated by the MNDO method. Clusters BNC₅₈ and B₂N₂C₅₆ are formed by replacement of carbon atoms participating in one or two of the most distant oppositely lying (6,6)-type C−C bonds in fullerene C₆₀ by B and N atoms. In one of the two studied isomers of the B₂N₂C₁₁₆ dimer, the monomers are linked by the four-membered carbon cycle, while the heteroatoms form the most distant oppositely lying bonds of the dimer. In the other isomer of the B₂N₂C₁₁₆ dimer, as well as in the B₄N₄C₁₁₂ dimer and B₆N₆C₁₆₈ trimer, the monomers are linked by four-membered B₂N₂ cycles with alternation of the atoms. For all the systems studied, the optimum geometric parameters, heats of formation, ionization potentials, and atomic charges were calculated. Dimerization energies of heterofullerenes BNC₅₈ and B₂N₂C₅₆ lie in the range from 33 to 49 kcal mol⁻¹. It was found that the B₂N₂C₁₁₆ dimer, in which the monomers are linked by the four-mernbered carbon cycle, is the most stable system. In the case of B₂N₂C₅₆ trimerization, the energy gain (compared to the triple monomer energy) is about twice as large as the dimerization energy. Molecular structure of the quasi-linear [B₂N₂C₅₆) _n macromolecule was simulated, and extended Hückel calculations of its energy band structure by the crystal orbital method were performed. It was found that the electron energy spectrum is of semiconducting type (the band gap is equal to 1.27 eV. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 3, pp. 431–435, March, 1999.     

DFT,AIM, and NBO analyses of 1-methyl-2-thioxoimidazolidin-4-one tautomers and their complexes with iodine

Thioimidazoline derivatives can be used to treat hyperthyroidism due to their ability to make complexes with iodine. In this research designed to find new structures with the same ability, 1-methyl-2-thioxoimidazolidin-4-one (MTIO) and the structures of MTIO tautomers (5 tautomers), their isomers (total 9 isomers) and their complexes with iodine are optimized using the B3LYP method with two different basis sets to obtain their molecular parameters, relative energies, and vibrational frequencies. The relative energies show that in all tautomers and complexes, ketone and thione forms are more stable than enol and thienol forms, and also Z isomers are more stable than E isomers. Moreover, the NBO calculation is carried out for tautomers and complexes to obtain atomic charges and acceptor-donor interactions. These results confirm the ability of MTIO tautomers to form complexes and show that the planar complexes have more effective interaction than the perpendicular complexes. The essence and important complexation properties are also calculated and confirmed using the AIM analysis.     

Structural and solvent effects on the 13C and 15N NMR chemical shifts of indoloquinoline alkaloids: experimental and DFT study

Indoloquinoline alkaloids represent an important class of antimalarial, antibacterial and antiviral compounds. They have been shown to bind to DNA via intercalation preferentially at GC-rich sequences containing nonalternating CC sites. The stability of complexes formed with biological macromolecules depends on noncovalent binding. In the present study, the ability of indoloquinolines to form intermolecular interactions with solvents was investigated by using NMR spectroscopy and density functional theory (DFT) (B3LYP/6-31G**) calculations. NMR data measured for indoloquinoline bases and the corresponding hydrochlorides are discussed in relation to the structure. DFT calculations of shielding constants in vacuo and in solution allowed the investigation of the influence of the environment on the NMR parameters. Calculations incorporating solvent effects indicated significant changes in the anisotropy of the electron distribution, reflected in the span of the chemical shielding tensor (Omega = sigma11 - sigma33). Solvent effects on the span of the 13C and 15N shielding tensor depended on the type of atom and the data indicated a significant influence of solute-solvent interactions.     

Adsorption of the nitrosamine and thionitrosamine molecules as carcinogen compounds on the BN and B3Al N nanotubes: A DFT study

DFT calculations were performed to investigation of the influence of doping three atoms of aluminum on the electronic properties of the (4,0) zigzag boron nitride nanotube (BNNT). Also, adsorption properties of nitrosamine (NA) and thionitrosamine (TNA) molecules as carcinogen agents onto BN and B_Al3N nanotubes were studied. The results show that the B_3AlN nanotube is the most energetically favorable candidates for adsorption of these molecules. Also, B(B_3Al)NNT/TNA complexes are more stable than B(B_3Al)NNT/NA complexes. The HOMO–LUMO gap, electronic chemical potential (μ), hardness (?), softness (S), the maximum amount of electronic charge (ΔN_max) and electrophilicity index (ω) for monomers and complexes in the gas and polar solvent phases were calculated. The results show that the conductivity and reactivity of BNNT increase by doping Al atoms instead of B atoms. Also, the interaction of NA and TNA molecules with BN and B_Al3N nanotubes results in significant changes in the electronic properties of nanotubes. Based on the natural bond orbital (NBO) analysis, in all complexes charge transfer occurs from NA and TNA molecules to nanotubes. Theory of atoms in molecules (AIM) was applied to characterize the nature of interactions in nanotubes. It is predicted that, BN and B_3AlN nanotubes can be used to as sensor for detection of NA and TNA molecules.     

Reinvestigation of intramolecular hydrogen bond in malonaldehyde derivatives: An ab initio,AIM and NBO study

The RAHB systems in malonaldehyde and its derivatives at MP2/ 6‐311++G(d,p) level of theory were studied and their intramolecular hydrogen bond energies by using the related rotamers method was obtained. The topological properties of electron density distribution in O? H···O intramolecular hydrogen bond have been analyzed in term of quantum theory of atoms in molecules (QTAIM). Correlations between the H‐bond strength and topological parameters are probed. The results of QTAIM clearly showed that the linear correlation between the electron density distribution at HB critical point and RAHB ring critical point with the corresponding hydrogen bond energies was obtained. Moreover, it was found a linear correlation between the electronic potential energy density, V(r_cp), and hydrogen bond energy which can be used as a simple equation for evaluation of HB energy in complex RAHB systems. Finally, the similar linear treatment between the geometrical parameters, such as O···O or O? H distance, and Lp(O)→σ*_OH charge transfer energy with the intramolecular hydrogen bond energy is observed. © 2010 Wiley Periodicals, Inc., Int J Quantum Chem, 2011     

A Theoretical Study of NBO,NICS, and 14N NQR Parameters of Adenine Tautomers in the Gas Phase via DFT

The natural bond orbital (NBO) analysis, nucleus independent chemical shift (NICS), and ¹⁴N NQR parameters of the most stable tautomers of adenine in the gas phase were predicted using density functional theory method. The NBO analysis revealed that the resonance interaction between lone pair of the nitrogen atom and empty non‐Lewis NBO increases with increasing the p character of the nitrogen lone pair. The present investigation indicated the π clouds in both the considered heterocyclic rings containing six electrons, and these tautomers has the aromatic character. The NICS study utilizing the gauge‐invariant atomic orbital method showed that there are diatropic currents in the heterocyclic rings of the tautomers, so we determined the order of overall aromaticity of these tautomers. The results of NQR parameter calculations showed three parameters are effective on nuclear quadrupole coupling constant; the p character value of lone pair electrons of nitrogens, and the related occupancies and whenever, the lone pair electrons of nitrogens participate in the formation of chemical bond and/or π system of the ring, the q_zz and consequently its χ decreases.     

β‐Aminoacrolein: An ab initio,AIM and NBO study

The molecular structure and the intramolecular hydrogen bonding of β‐aminoacrolein and its simple derivatives were investigated at the MP2 and B3LYP levels of theory using the standard 6‐311++G(d, p) basis set. The “atoms in molecules” or AIM theory of Bader which is based on topological properties of the electron density (ρ), was used. Additionally, an analysis of the critical points was performed to study the nature hydrogen bonding in these systems. Natural bond orbital (NBO) analysis was also carried out for to better comprehend the nature of the intramolecular interactions in β‐aminoacrolein and its derivatives. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

DFT and NBO theoretical study of protonation of tri-tert-butoxysilanethiol and its anion

Density functional theory (DFT) calculations of geometry optimization at the B3LYP/6-31+G(d) level have been made for the tri-tert-butoxysilanethiolato anion, i.e. (^tBuO)₃SiS⁻ [coded A], and its O- and S-protonated derivatives: (^tBuO)₃SiSH (B), (^tBuO⁺H)(^tBuO)₂SiS⁻ (C), (^tBuO)₃SiSH₂⁺ (D) and (^tBuO⁺H)(^tBuO)₂SiSH (E). Energy balances for the protonation of silanethiol B in vacuum and the aqueous phase (using the polarized continuum IPCM model) were considered. Based on theoretical calculations, it was found that protonation on the oxygen atom causes Si–S bond shortening, while protonation on the sulfur atoms results in elongation of the bond. Changeability of the Si–S and Si–O bond lengths was analyzed in terms of the natural bond orbital (NBO) method. Principal electronic delocalizations caused by anomeric effects of the sulfur and oxygen lone pairs interaction with the anti-bonding Si–O orbitals were confirmed by the correlation between the shortening of the Si–S bond length and the elongation of the Si–O bond within the structures of the species A–E.     

Adsorption and diffusion on a phosphorene monolayer: a DFT study

A computational study of the adsorption and diffusion behavior of alkali and alkaline earth metal atoms on a phosphorene monolayer is reported. Our calculations were performed within the framework of density functional theory using the Perdew–Burke–Ernzerhof functional and projector augmented wave potentials, as derived from the generalized gradient approximation. Our binding energy calculations for various potential adsorption sites showed that the site located above the center of a triangle formed by three surface phosphorus atoms is the most attractive to all adatoms. In addition, simulation of the diffusion of adatoms across the surface of the phosphorene monolayer showed that the diffusion is anisotropic, with K having the lowest diffusion barrier (0.02 eV along the zigzag pathway). To the best of our knowledge, this is the lowest diffusion barrier of any metal adatom on a single layer of phosphorene. While phosphorene exhibited significantly better adatom adsorption and diffusion than graphene, it also showed a reduced storage capacity compared to graphene, most probably due to the structural distortion induced by the oversaturated phosphorene surface. This finding strongly suggests that a phosphorene–graphene hybrid system could be employed as a promising high-capacity ion anode.     

Boosting BeONT Reactivity with HCN by Calcium and Magnesium Doping: A DFT Investigation of Electronic Structure,AIM, NMR,NQR and NBO Analysis

Adsorption of the HCN molecule is very important in environment and industrial applications. The BeONT may be good candidate for HCN capture because of large surface. Unfortunately, BeONT shows limited HCN detection. Therefore, we investigate the possibility of HCN adsorption on Ca and Mg-doped BeONT by density functional theory calculations. It was found that HCN adsorption on doped nanotube has relatively higher adsorption energy as compared with the perfect one. Furthermore, there exists a strong adsorption between HCN molecule and doped nanotubes, which exhibited more active interaction and larger net charge transfer than that of pristine nanotube. As well as, calculated geometrical parameters and electronic properties for studied systems indicate that the Ca-doped BeONT and Mg-doped BeONT present high sensitivity to HCN, compared with the pristine BeONT. Theoretical results reveal that the adsorption of the HCN on the doped nanotube is influenced on the electronic conductance of the doped-BeONT. Therefore, Ca and Mg-doped nanotube can be considered as promising sensor for detecting HCN molecule. According to NBO analysis, electron flow is spontaneous from doped nanotube to HCN molecule.     

The effect of boron nitride nanotubes size on the HArF interaction by NBO and AIM analysis

Three molecules HArF@BNNT(5,0), HArF@BNNT(6,0), and HArF@BNNT(7,0) have been formed by HArF encapsulated in boron nitride nanotubes (BNNTs) with different sizes. Due to the interaction between the HArF and the BNNTs, the H? Ar bond lengths are in a decrease trend, while the Ar? F bond lengths are in an increase trend compared with those of HArF. To investigate the nature of the interaction between H and Ar and the interaction between Ar and F, the quantum theory of “atoms in molecules” was carried out. The Laplacian (?²ρ_b) values of H? Ar suggest that the covalent interaction plays a key role in the H‐Ar interaction. For Ar? F, the results indicate that the Ar‐F interaction has a dominant noncovalent character. Moreover, the results obtained from the ratio of the kinetic‐energy density (G_b) and the potential‐energy density (V_b) (?G_b/V_b) and the total energy density (H_b) are in good agreement with that of ?²ρ_b values. In addition, the results of natural bond orbital charge and electron density difference between the HArF and BNNTs show that less electrons transfer away from the HArF to BNNTs with the gradual increase in the diameters of the BNNTs. © 2014 Wiley Periodicals, Inc.