Theoretical study on interactions between thiophene/dibenzothiophene/cyclohexane/toluene and 1-methyl-3-octylimidazolium tetrafluoroborate

It is critically important to understand the interactions between thiophene/dibenzothiophene/cyclohexane/toluene and 1-methyl-3-octylimidazolium tetrafluoroborate ([C8MIM]⁺[BF₄]^?) due to desulfurization by ionic liquids. In this work, the structures of thiophene, dibenzothiophene, cyclohexane, toluene, [C8MIM]⁺[BF₄]^?, [C8MIM]⁺[BF₄]^?-thiophene, [C8MIM]⁺[BF₄]^?-dibenzothiophene, [C8MIM]⁺[BF₄]^?-cyclohexane, and [C8MIM]⁺[BF₄]^?-toluene were optimized systematically at the GGA/PW91/DNP level, and the most stable geometries were performed by NBO and AIM analyses. It was found that [BF₄]^? anion tends to locate near C2–H2 and four hydrogen bonds between [C8MIM]⁺ and [BF₄]^? form in [C8MIM]⁺[BF₄]^?. There exist hydrogen bonds and C–H···π interactions between [C8MIM]⁺[BF₄]^? and thiophene/cyclohexane/toluene, while the hydrogen bonding interactions, π···π and C–H···π interactions occur between [C8MIM]⁺[BF₄]^? and dibenzothiophene confirmed by NBO and AIM analyses. The interaction energies between [C8MIM]⁺[BF₄]^? and thiophene, dibenzothiophene, cyclohexane, toluene are 18.83, 20.93, 6.83, 12.99 kcal/mol, showing the preferential adsorption of dibenzothiophene and thiophene by ionic liquid, in agreement with the experimental results of dibenzothiophene and thiophene extraction by [C8MIM]⁺[BF₄]^?.     

二苯并噻吩及其氧化物与离子液体相互作用的理论研究

采用密度泛函理论方法比较了DBT/DBTO₂和[BMIM]⁺[PF₆]^-/[BMIM]⁺[BF₄]^-的相互作用。对最稳定的[BMIM]⁺[PF₆]^-、[BMIM]⁺[PF₆]^--DBT、[BMIM]⁺[PF₆]^--DBTO₂、[BMIM]⁺[BF₄]^-、[BMIM]⁺[BF₄]^--DBT、[BMIM]⁺[BF₄]^--DBTO₂进行了NBO和AIM分析。结果表明,DBT和[BMIM]⁺[PF₆]^-/[BMIM]⁺[BF₄]^-中的咪唑环彼此相互平行,NBO和AIM分析表明它们之间发生了π-π相互作用。H1'和H9'形成的F…H氢键有利于π-π堆积作用的形成。DBTO₂倾向于趋近C2-H2和甲基基团形成O…H相互作用;DBTO₂优先吸附在[BMIM]⁺[PF₆]^-/[BMIM]⁺[BF₄]^-。在模拟油中,[BMIM]⁺[PF₆]^-和[BMIM]⁺[BF₄]^-离子液体对DBTO₂的萃取能力大于DBT,其原因是可能是DBTO₂具有较大的极性和O…H与F…H的氢键作用。     

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.     

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.     

Adsorption of nitrogen dioxide on C30B15N15 heterofullerene: AIM and NBO study via DFT

Adsorption of nitrogen dioxide on the exterior surface of C₃₀B₁₅N₁₅ through B atom of adsorbent is studied using the DFT-B3LYP/6-31G^* level. Independent to the orientation, the NO₂ molecule presents strong chemisorption with large charge transfer from C₃₀B₁₅N₁₅ to adsorbed gas. NBO (natural electron configuration of adsorbing B atom), AIM analysis and spin density distributions indicate that the chemisorption of NO₂ on the C₃₀B₁₅N₁₅ in all studied complexes are covalent in nature.     

β‐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     

Physical and electrochemical properties of 1-alkyl-3-methylimidazolium tetrafluoroborate for electrolyte

Three kinds of ionic liquids, 1-alkyl-3-methylimidazolium tetrafluoroborate (n=2–4), were prepared and fundamental properties of ionic liquids and those mixed with industrially used organic solvents (PC, GBL and AN) were investigated compared to solid salts, TEMABF₄. It was found that degree of ionization of the ionic liquids were almost same as that of TEMABF₄ from the conductivity measurement in diluted system of PC. The ionic liquids and the organic solvents intermingle with each other. Some enhancement in conductivity was observed compared to TEMABF₄.     

DFT study of the monomers and dimers of 2-pyrrolidone: equilibrium structures, vibrational, orbital, topological, and NBO analysis of hydrogen-bonded interactions

A computational study of the monomers and hydrogen-bonded dimers of 2-pyrrolidone was executed at different DFT levels and basis sets. The above dimeric complexes were treated theoretically to elucidate the nature of the intermolecular hydrogen bonds, geometry, thermodynamic parameters, interaction energies, and charge transfer. The processes of dimer formation from monomers and concerted reactions of double proton transfer were considered. The evolution of geometry, vibrational frequencies, charge distribution, and AIM properties in going from monomers to dimers was systematically followed. The solvent effects upon dimer formation were investigated in terms of the self-consistent reaction field (SCRF Onsager model). For the monomers and three dimers, vibrational frequencies were calculated and the changes in frequencies of the vibrations most sensitive to complexation were discussed. The orbital interactions were shown to lengthen the X-H (X = N, O) bond and lower its vibrational frequency (a red shift). To better understand the nature of the corresponding intermolecular interactions, we performed natural bond orbital (NBO) analysis. Topological analysis of electron density at bond critical points (BCP) was executed for complex molecules using the Bader's atoms in molecules (AIM) theory. The interaction energies were calculated, and the basis set superposition errors (BSSE) were estimated systematically. Satisfactory correlations between the structural parameters, interaction energies, and electron density characteristics at BCP were found.     

Theoretical analysis of bonding properties in 1,1-dicyanocyclopentane and 1,1-dichlorocyclopentane by applying the AIM and NBO approaches

Using Bader’s quantum-topological theory of atoms in molecules (AIM) and Weinhold’s Natural Bond Orbital (NBO) analysis we could rationalize the impact of the geminal substitution by C≡N and Cl on the geometry and electronic structure of the cyclopentane ring in 1,1-dicyanocyclopentane (DCCP) and 1,1-dichlorocyclopentane (DClCP). Among the crucial results we obtained are: 1. The topological quantities, particularly the bond ellipticity of 0.035 for the C–CN bond, indicate that this bond possesses higher bond order than a single bond. This conclusion is clearly supported by the NBO results. 2. The AIM theory as well as the NBO analysis confirm uniformly the non-linearity of the C–C≡N moiety. 3. Regardless the quantum mechanical method that has been employed for a variety of nitriles the sign of the Laplacian of the charge density, ∇²ρ(r), of the C≡N group changes its sign from negative to positive upon moving from the triple zeta to the double zeta basis set. A possible explanation for this striking behavior has been provided. By invoking the AIM and NBO approaches the different endocyclic C–C bond lengths in DCCP and DClCP as a consequence of the geminal substitution could be explained. Also the variations of these bond lengths upon moving from the more stable C _s to the energetically less favorable C ₂ conformer of both compounds could be rationalized. For the purpose of comparison and verification of some findings of this work, we also carried out AIM and NBO calculations on various related cyclic and non-cyclic compounds.     

DFT studies and AIM analysis of AlN-polycycles

The structure and properties of AlN-polycycles were studied by DFT (density functional theory) method. The results of calculations were obtained at B3LYP/6-311G(d, p) level on model species. Topological parameters such as electron density, its Laplacian, kinetic electron energy density, potential electron energy density, and total electron energy density at the ring critical points (RCP) from Bader’s ‘Atoms in molecules’ (AIM) theory were analyzed in detail. These results indicate a good correlation between ρ(3, +1), G(r), H(r), and V(r) averaged values and hardness of AlN-polycycles. The aromaticity of all molecules has been studied by nucleus-independent chemical shift. There is a linear correlation between ΣNICS(0.0)_molecule values and polarizability.     

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.     

DFT,AIM, and NBO study of the interaction of simple and sulfur-doped graphenes with molecular halogens,CH3OH,CH3SH,H2O,and H2S

Graphene is an important material in adsorption processes because of its high surface. In this work, the interactions between graphene (G), S-doped graphene (SG), and 2S-doped graphene (2SG) with eight small molecules including molecular halogens, CH₃OH, CH₃SH, H₂O, and H₂S were studied using density functional theory calculations. The adsorption energies showed that the SG was the best adsorbent, fluorine was the best adsorbate, and all molecular halogens were adsorbed on graphenes better than the other molecules. Most adsorption processes in the gas phase were exothermic with small positive ΔG _ads. Moreover, the solvent effect on the adsorption process was examined and all ΔH _ads and ΔG _ads values for adsorption processes tended to be more negative in all solvents. Therefore, most adsorption processes in the solvents were thermodynamically favorable. The second order perturbation energies obtained from NBO calculations confirmed that the interactions between molecular halogens and our molecules had more strength than those of other molecules. The Laplacian of ρ values obtained from the AIM calculations indicated that the type of interaction in all our complexes was one of closed shell interaction. The MO results and DOS plots also revealed that sulfur doping could increase the conductivity of graphene and this conductivity was enhanced more when they interacted with molecular halogens.     

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     

Control of a photoswitching chelator by metal ions: DFT,NBO, and QTAIM analysis

Ability of aroylhydrazones to change conformation upon interaction with light makes them promising candidates for molecular switches. Isomerization can be controlled through complexation with selected metal ions which bind with different affinity. N′‐[1‐(2‐hydroxyphenyl)ethyliden]iso‐nicotinoylhydrazide (HAPI) is an example of a dual‐wavelenght photoswitching molecule, whose complexation with metal ions was recently experimentally investigated (Franks et al. J. Inorg. Chem. 2014, 53, 1397). In this contribution, complexes between HAPI and K⁺, Ca²⁺, Mn²⁺, Fe²⁺, Fe³⁺, Cu⁺, Cu²⁺, and Zn²⁺ ions were investigated using Density Functional Theory, Natural Bond Order analysis, and Quantum Theory of Atoms in Molecules. The most important parameters that determine complex stability are found to be ion radius and charge transferred from ligands to the ion: smaller ion radii and larger CT values characterize formation of more stable complexes. Our results explain experimentally observed effect of different metal ions on photoisomerization through determination of metal ion affinity (MIA): photoisomerization is inhibited if MIA exceeds 100 kcal/mol; for MIA between 50 and 100 kcal/mol excess of metal ions prevents isomerization, whereas in case of MIA below 50 kcal/mol metal ions have no influence on light–HAPI interaction. © 2015 Wiley Periodicals, Inc.     

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.     

Surface tension of 1-alkyl-3-methylimidazolium based ionic liquids with trifluoromethanesulfonate and tetrafluoroborate anion

The amount of available accurate experimental data on the surface tension of ionic liquids is still limited; in many cases the data are rare or even absent. In the present study, air-liquid interfacial tension data were determined experimentally for five 1-C_n-3-methylimidazolium based ionic liquids (n = 2, 4, and 6), three with trifluoromethanesulfonate and two with tetrafluoroborate anion, at atmospheric pressure in the temperature range from 268 to 356 K. The resultant surface tension data are average values of the measurements repeated many times at each set point temperature. The accuracy of the results, was confirmed by employing the Wilhelmy plate and the du Noüy ring methods in parallel, using the Krüss K100MK2 tensiometer. For the Wilhelmy plate data the combined standard uncertainty is estimated to be about 0.05 mN m⁻¹. The data obtained by du Noüy method show about up to seven times greater scatter than those obtained by the Wilhelmy plate method. To the 50 up to now published surface tension values for the five studied ionic liquids the present study adds further 175 data points. In contrast to that of n-alkanes, the surface tension of 1-alkyl-3-methylimidazolium based ionic liquids decreases and their surface entropy increases with the cation alkyl chain length.     

Study of noncovalent interactions of end-caped sulfur-doped carbon nanotubes using DFT,QTAIM, NBO and NCI calculations

In the current study, the interactions of carbon nanotube and sulfur-doped carbon nanotubes (SCNTs) with methanol, methanethiol, water and dihydrogen sulfide at on-body and dead-end positions of nanotubes have been studied. Interaction energies in the gas and solvent (via PCM model) were calculated using density functional theory calculations. Atomic charges, interaction energies, electron densities and their Laplacians at bond critical points have been calculated. Moreover, noncovalent interaction isosurfaces have been visualized using NCI index calculations. Interactions in gaseous phase were more favorable than those in solvent phase, and among considered solvents (benzene, chloroform and cyclohexane), cyclohexane showed the most preferred interactions. In addition, oxygen-bearing molecules (methanol and water) showed more favorable interactions compared with sulfur-bearing ones. NBO analyses revealed the stronger donor–acceptor interactions with methanol and methanethiol. QTAIM calculation results indicated the reasonable electron densities at BCPs, and the Laplacians of electron densities showed ionic-like (closed shell) interactions. Moreover, isosurfaces of these interactions were also studied to depict the interaction surfaces, and DOS plots for SCNTs were obtained to define their HOMO–LUMO levels and electric conduction properties. The increasing of the global softness and decreasing of total hardness was resulted by sulfur doping of nanotubes, which causes the heterodoped nanotubes to become less electrophilic species.