First-principles vdW-DF investigation on the interaction between the oxazepam molecule and C60 fullerene

The interaction between oxazepam and C₆₀ fullerene was explored using first-principles vdW-DF calculations. It was found that oxazepam binds weakly to the fullerene cage via its carbonyl group. The binding of oxazepam to C₆₀ is affected drastically by nonlocal dispersion interactions, while vdW forces affect the corresponding geometries only a little. Furthermore, aqueous solution affects the geometries of the oxazepam approaching to fullerene slightly, while oxazepam binds slightly farther away from the nanocage. The results presented provide evidence for the applicability of the vdW-DF method and serve as a practical benchmark for the investigation of host–guest interactions in biological systems.

Dynamic motion of La atom inside the C74 (D 3h) cage: a relativistic DFT study

The interaction between lanthanum atom (La) and C₇₄ (D _3h) was investigated by all-electron relativistic density function theory (DFT). With the aid of the representative patch of C₇₄ (D _3h), we studied the interaction between C₇₄ (D _3h) and La and obtained the interaction potential. Optimized structures show that there are three equivalent stable isomers, with La located about 1.7 Å off center. There is one transition state between every two stable isomers. According to the minimum energy pathway, the possible movement trajectory of La atoms in the C₇₄ (D _3h) cage was explored. The calculated energy barrier for La atoms moving from the stable isomer to the transition state is 18.4 kcal mol^?1. In addition, the dynamic NMR spectra of La@C₇₄ according to the trajectory was calculated.

Synthesis and biological evaluation of cationic fullerene quinazolinone conjugates and their binding mode with modeled Mycobacterium tuberculosis hypoxanthine-guanine phosphoribosyltransferase enzyme

The present work reports a series of novel cationic fullerene derivatives bearing a substituted-quinazolinone moiety as a side arm. Fullerene-quinazolinone conjugates synthesized using the 1,3-dipolar cycloaddition reaction of C₆₀ with azomethine ylides generated from the corresponding Schiff bases of substituted quinazolinone were characterized by elemental analysis, FT-IR, ¹H NMR, ¹³C NMR and ESI-MS and screened for their antibacterial activity against Mycobacterium tuberculosis (H ₃₇ Rv strain). All the compounds exhibited significant activity with the most effective having MIC in the range of 1.562–3.125 μg/mL. Compound 9f exhibited good biological activity compared to standard drugs. We developed a computational strategy based on the modeled M. tuberculosis hypoxanthine-guanine phosphoribosyltransferase (HGPRT) using homology modeling techniques and studied its binding pattern with synthesized fullerene derivatives. We then explored the surface geometry of the protein to place the cage adjacent to the active site while optimizing its quinazolinone side arm to establish H bonding with active site residues.

Theoretical study on aluminum carbide endohedral fullerene-Al4C@C80

The possibility of a new endohedral fullerene with a trapped aluminum carbide cluster, Al₄C @C₈₀-I _h , was theoretical investigated. The geometries and electronic properties of it were investigated using density functional theory methods. The Al₄C unit formally transfers six electrons to the C₈₀ cage which induces stabilization of Al₄C@C₈₀. A favorable binding energy, relatively large HOMO-LUMO gap, electron affinities and ionization potentials suggested the Al₄C@C₈₀ is rather stable. The analysis of vertical ionization potential and vertical electron affinity indicate Al₄C@C₈₀ is a good electron acceptor.

The effect of C-vacancy on hydrogen storage and characterization of H2 modes on Ti functionalized C60 fullerene A first principles study

Density functional theory calculations were performed to examine the effect of a C vacancy on the physisorption of H₂ onto Ti-functionalized C₆₀ fullerene when H₂ is oriented along the x-, y-, and z-axes of the fullerene. The effect of the C vacancy on the physisorption modes of H₂ was investigated as a function of H₂ binding energy within the energy window (?0.2 to ?0.6 eV) targeted by the Department of Energy (DOE), and as functions of a variety of other physicochemical properties. The results indicate that the preferential orientations of H₂ in the defect-free (i.e., no C vacancy) C₆₀TiH₂ complex are along the x- and y-axes of C₆₀ (with adsorption energies of ?0.23 and ?0.21 eV, respectively), making these orientations the most suitable ones for hydrogen storage, in contrast to the results obtained for defect-containing fullerenes. The defect-containing (i.e., containing a C vacancy) C₅₉TiH₂ complex do not exhibit adsorption energies within the targeted energy range. Charge transfer occurs from Ti 3d to C 2p of the fullerene. The binding of H₂ is dominated by the pairwise support–metal interaction energy E(i)^Cn...Ti, and the role of the fullerene is not restricted to supporting the metal. The C vacancy enhances the adsorption energy of Ti, in contrast to that of H₂. A significant reduction in the energy gap of the pristine C₆₀ fullerene is observed when TiH₂ is adsorbed by it. While the C _n fullerene readily participates in nucleophilic processes, the adjacent TiH₂ fragment is available for electrophilic processes.

Direct ab initio study on the rate constants of radical C2(A3Πu) + C3H8 reaction

The mechanism and kinetics of the radical ³C₂?+?C₃H₈ reaction have been investigated theoretically by direct ab initio kinetics over a wide temperature range. The potential energy surfaces have been constructed at the CCSD(T)/B3//UMP2/B1 levels of theory. The electron transfer was also analyzed by quasi–restricted orbital (QRO) in detail. It was shown that all these channels proceed exclusively via hydrogen abstraction. The overall ICVT/SCT rate constants are in agreement with the available experimental results. The prediction shows that the secondary hydrogen of C₃H₈ abstraction by ³C₂ radical is the major pathway at low temperatures (below 700 K), while as the temperature increases, the primary hydrogen of C₃H₈ abstraction becomes more important and more favorable. A negative temperature dependence of the rate constants for the reaction of ³C₂?+?C₃H₈ was observed. The three–(k ₃) and four–parameter (k ₄) rate-temperature expressions were also provided within 243–2000 K to facilitate future experimental studies.

B30H8, B39H9 2−, B42H10, B48H10, and B72H12: polycyclic aromatic snub hydroboron clusters analogous to polycyclic aromatic hydrocarbons

Calculations performed at the ab initio level using the recently reported planar concentric π-aromatic B₁₈H₆ ²⁺(1) [Chen Q et al. (2011) Phys Chem Chem Phys 13:20620] as a building block suggest the possible existence of a new class of B_3n H _m polycyclic aromatic hydroboron (PAHB) clusters—B₃₀H₈(2), B₃₉H₉ ^2?(3), B₄₂H₁₀(4/5), B₄₈H₁₀(6), and B₇₂H₁₂(7)—which appear to be the inorganic analogs of the corresponding C _n H _m polycyclic aromatic hydrocarbon (PAHC) molecules naphthalene C₁₀H₈, phenalenyl anion C₁₃H₉ ^?, phenanthrene/anthracene C₁₄H₁₀, pyrene C₁₆H₁₀, and coronene C₂₄H₁₂, respectively, in a universal atomic ratio of B:C?=?3:1. Detailed canonical molecular orbital (CMO), adaptive natural density partitioning (AdNDP), and electron localization function (ELF) analyses indicate that, as they are hydrogenated fragments of a boron snub sheet [Zope RR, Baruah T (2010) Chem Phys Lett 501:193], these PAHB clusters are aromatic in nature, and exhibit the formation of islands of both σ- and π-aromaticity. The predicted ionization potentials of PAHB neutrals and electron detachment energies of small PAHB monoanions should permit them to be characterized experimentally in the future. The results obtained in this work expand the domain of planar boron-based clusters to a region well beyond B₂₀, and experimental syntheses of these snub B_3n H _m clusters through partial hydrogenation of the corresponding bare B_3n may open up a new area of boron chemistry parallel to that of PAHCs in carbon chemistry.

A B3LYP and MP2(full) theoretical investigation into the strength of the C-NO2 bond upon the formation of the molecule-cation interaction between Na+ and the nitro group of nitrotriazole or its methyl derivatives

The changes of bond dissociation energy (BDE) in the C–NO₂ bond and nitro group charge upon the formation of the molecule-cation interaction between Na⁺ and the nitro group of 14 kinds of nitrotriazoles or methyl derivatives were investigated using the B3LYP and MP2(full) methods with the 6-311++G**, 6-311++G(2df,2p) and aug-cc-pVTZ basis sets. The strength of the C–NO₂ bond was enhanced in comparison with that in the isolated nitrotriazole molecule upon the formation of molecule-cation interaction. The increment of the C–NO₂ bond dissociation energy (ΔBDE) correlated well with the molecule-cation interaction energy. Electron density shifts analysis showed that the electron density shifted toward the C-NO₂ bond upon complex formation, leading to the strengthened C-NO₂ bond and the possibly reduced explosive sensitivity.

A B3LYP and MP2(full) theoretical investigation into the strength of the C–NO2 bond upon the formation of the intermolecular hydrogen-bonding interaction between HF and the nitro group of nitrotriazole or its methyl derivatives

The changes of bond dissociation energy (BDE) in the C–NO₂ bond and nitro group charge upon the formation of the intermolecular hydrogen-bonding interaction between HF and the nitro group of 14 kinds of nitrotriazoles or methyl derivatives were investigated using the B3LYP and MP2(full) methods with the 6-311++G**, 6-311++G(2df,2p) and aug-cc-pVTZ basis sets. The strength of the C–NO₂ bond was enhanced and the charge of nitro group turned more negative in complex in comparison with those in isolated nitrotriazole molecule. The increment of the C–NO₂ bond dissociation energies correlated well with the intermolecular H-bonding interaction energies. Electron density shifts analyses showed that the electron density shifted toward the C–NO₂ bond upon complex formation, leading to the strengthened C–NO₂ bond and the possibly reduced explosive sensitivity.

Scytonin,a novel cyanobacterial photoprotective pigment: calculations of Raman spectroscopic biosignatures

Experimental studies have shown that chemical vapor decomposition method by using C₂H₂/C₂H₄ as carbon source could dramatically decrease the defects in prepared CNT. However, the inherent mechanism with regards to reduction of defects is quite unclear. In the present paper, density functional theory is used to study the healing process of CNT with divacancy defect by C₂H₂/C₂H₄ molecule. The healing processes undergo three evolution steps: (i) the chemisorption of the first C₂H₂/C₂H₄ molecule on defective CNT; (ii) the insertion of C atoms from C₂H₂/C₂H₄ molecule into defective CNT; (iii) the removal of the H atoms on CNT, forming perfect CNT. The estimated adsorption energy barrier of C₂H₂/C₂H₄ molecules on defective CNT is within the range from 1.10 to 1.63 eV, and the eventual formation of CNT is strongly exothermic (4.40/4.54 eV in (8, 0) CNT). In light of the unique conditions of CNT synthesis, i.e., high temperature in a closed container, such healing processes could most likely take place. Therefore, we propose that during CNT synthesis procedures, both C₂H₂ and C₂H₄ could act as a carbon source and the defect healer.

A comparative theoretical investigation into the strength of the trigger-bond in the Na+, Mg2+ and HF complexes involving the nitro group of R–NO2 (R = –CH3, –NH2 and –OCH3) or the C = C bond of (E)-O2N–CH = CH–NO2

A comparative theoretical investigation into the change in strength of the trigger-bond upon formation of the Na⁺, Mg²⁺ and HF complexes involving the nitro group of RNO₂ (R?=? –CH₃, –NH₂, –OCH₃) or the C?=?C bond of (E)-O₂N–CH?=?CH–NO₂ was carried out using the B3LYP and MP2(full) methods with the 6-311++G**, 6-311++G(2df,2p) and aug-cc-pVTZ basis sets. Except for the Mg²⁺?π system with (E)-O₂N–CH?=?CH–NO₂ (i.e., C₂H₂N₂O₄?Mg²⁺), the strength of the trigger-bond X–NO₂ (X?=?C, N or O) was enhanced upon complex formation. Furthermore, the increment of bond dissociation energy of the X–NO₂ bond in the Na⁺ complex was far greater than that in the corresponding HF system. Thus, the explosive sensitivity in the former might be lower than that in the latter. For C₂H₂N₂O₄?Mg²⁺, the explosive sensitivity might also be reduced. Therefore, it is possible that introducing cations into the structure of explosives might be more efficacious at reducing explosive sensitivity than the formation of an intermolecular hydrogen-bonded complex. AIM, NBO and electron density shifts analyses showed that the electron density shifted toward the X–NO₂ bond upon complex formation, leading to a strengthened X–NO₂ bond and possibly reduced explosive sensitivity.

Stability and isomerization of complexes formed by metal ions and cytosine isomers in aqueous phase

We present a systematic study of the stability of the formation of complexes produced by four metal ions (M^+/2+) and 14 cytosine isomers (C_n). This work predicts theoretically that predominant product complexes are associated with higher-energy C₄M^+/2+ and C₅M^+/2+ rather than the most stable C₁M^+/2+. The prediction resolves successfully several experimental facts puzzling two research groups. Meanwhile, in-depth studies further reveal that direct isomerization of C₁?C₄ is almost impossible, and also that the isomerization induced by either metalation or hydration, or by a combination of the two unfavorable. It is the single water molecule locating between the H1(?N1) and O2 of the cytosine that plays the dual roles of being a bridge and an activator that consequently improves the isomerization greatly. Moreover, the cooperation of divalent metal ion and such a monohydration actually leads to an energy-free C₁←C₄ isomerization in the gas phase. Henceforth, we are able to propose schemes inhibiting the free C₁←C₄ isomerization, based purely on extended hydration at the divalent metal ion.

Relationship between structure and entropy contributions in an anthraquinone mercapto derivative

The structural and thermodynamic properties of an anthraquinone derivative were studied by means of quantum-chemical calculations. Conformational analysis using ab initio and density functional theory methods revealed 14 low-energy conformers. In order to discuss similarities and differences in entropy of the conformers, the rotational and vibrational contributions to entropy were correlated with changes in conformer structure. The component of the moment of inertia perpendicular to the molecular plane gives significant input to ΔS _rot , whereas the largest contributions to the ΔS _vib have vibrations associated with the τ _S1C20 coordinate.

DFT study of the mechanism of the reaction of aminoguanidine with methylglyoxal

Density functional theory calculations were used in the theoretical investigation of the adsorption properties of sumanene towards molecules considered as common air pollutants: CO, CO₂ and NH₃. The insignificant perturbation of sumanene after adsorption and the adsorption energies obtained indicate a physisorption mechanism. It was shown that, contrary to carbon nanotubes, sumanene is able to adsorb CO molecules, and that adsorption of CO₂ by sumanene is stronger than adsorption of CO₂ by C₆₀. To better understand the adsorption characteristics of sumanene, density of states and natural bond order analyses were performed, which showed that chemical interactions exist and that these are more important mostly on the convex side. Better adsorption properties were obtained for the concave side as adsorption is dictated by physisorption mechanisms due to the specific bowl-shaped geometry of sumanene, because of which more negative charge is located precisely on the concave side. Molecular electrostatic potential surfaces were also used in order to better locate the adsorption sites and gain additional details about adsorption.

Theoretical description of halogen bonding – an insight based on the natural orbitals for chemical valence combined with the extended-transition-state method (ETS-NOCV)

In the present study we have characterized the halogen bonding in selected molecules H₃N–ICF₃ (1-NH ₃ ), (PH₃)₂C–ICF₃ (1-CPH ₃ ), C₃H₇Br–(IN₂H₂C₃)₂C₆H₄ (2-Br), H₂–(IN₂H₂C₃)₂C₆H₄ (2-H ₂ ) and Cl–(IC₆F₅)₂C₇H₁₀N₂O₅ (3-Cl), containing from one halogen bond (1-NH ₃ , 1-CPH ₃ ) up to four connections in 3-Cl (the two Cl–HN and two Cl–I), based on recently proposed ETS-NOCV analysis. It was found based on the NOCV-deformation density components that the halogen bonding C–X ^… B (X-halogen atom, B-Lewis base), contains a large degree of covalent contribution (the charge transfer to X ^… B inter-atomic region) supported further by the electron donation from base atom B to the empty σ*(C–X) orbital. Such charge transfers can be of similar importance compared to the electrostatic stabilization. Further, the covalent part of halogen bonding is due to the presence of σ-hole at outer part of halogen atom (X). ETS-NOCV approach allowed to visualize formation of the σ-hole at iodine atom of CF₃I molecule. It has also been demonstrated that strongly electrophilic halogen bond donor, [C₆H₄(C₃H₂N₂I)₂][OTf]₂, can activate chemically inert isopropyl bromide (2-Br) moiety via formation of Br–I bonding and bind the hydrogen molecule (2-H ₂ ). Finally, ETS-NOCV analysis performed for 3-Cl leads to the conclusion that, in terms of the orbital-interaction component, the strength of halogen (Cl–I) bond is roughly three times more important than the hydrogen bonding (Cl–HN).

DFT study on the reactions of ClO–/BrO– with RCl (R = CH3, C2H5, and C3H7) in gas phase

Gas-phase reactions of ClO^–/BrO^– with RCl (R = CH₃, C₂H₅, and C₃H₇) have been investigated in detail using the popular DFT functional BHandHLYP/aug-cc-pVDZ level of theory. As a result, our findings strongly suggest that the type of reaction is firstly initiated by a typical S_N2 fashion. Subsequently, two competitive substitution steps, named as S_N2-induced substitution and S_N2-induced elimination, respectively, would proceed before the initial S_N2 product ion-dipole complex separates, in which the former exhibits less reactivity than the latter. Those are consistent with relevant experimental results. Moreover, we have also explored reactivity difference for the title reactions in term of some factors derived from methyl group, p-π electronic conjugation, ionization energy (IE), as well as molecular orbital (MO) analysis.

The average local ionization energy as a tool for identifying reactive sites on defect-containing model graphene systems

In a continuing effort to further explore the use of the average local ionization energy $ \overline{\mathrm{I}}\left( \mathbf{r} \right) $ as a computational tool, we have investigated how well $ \overline{\mathrm{I}}\left( \mathbf{r} \right) $ computed on molecular surfaces serves as a predictive tool for identifying the sites of the more reactive electrons in several nonplanar defect-containing model graphene systems, each containing one or more pentagons. They include corannulene (C₂₀H₁₀), two inverse Stone-Thrower-Wales defect-containing structures C₂₆H₁₂ and C₄₂H₁₆, and a nanotube cap model C₂₂H₆, whose end is formed by three fused pentagons. Coronene (C₂₄H₁₂) has been included as a reference planar defect-free graphene model. We have optimized the structures of these systems as well as several monohydrogenated derivatives at the B3PW91/6-31G* level, and have computed their $ \overline{\mathrm{I}}\left( \mathbf{r} \right) $ on molecular surfaces corresponding to the 0.001 au, 0.003 au and 0.005 au contours of the electronic density. We find that (1) the convex sides of the interior carbons of the nonplanar models are more reactive than the concave sides, and (2) the magnitudes of the lowest $ \overline{\mathrm{I}}\left( \mathbf{r} \right) $ surface minima (the $ {{\overline{\mathrm{I}}}_{{\mathrm{S}\text{,}\min }}} $ ) correlate well with the interaction energies for hydrogenation at these sites. These $ {{\overline{\mathrm{I}}}_{{\mathrm{S}\text{,}\min }}} $ values decrease in magnitude as the nonplanarity of the site increases, consistent with earlier studies. A practical benefit of the use of $ \overline{\mathrm{I}}\left( \mathbf{r} \right) $ is that a single calculation suffices to characterize the numerous sites on a large molecular system, such as graphene and defect-containing graphene models.

Theoretical study on cooperative effects between X⋯N and X⋯Carbene halogen bonds (X = F,Cl,Br and I)

Quantum chemical calculations are performed to study the interplay between halogen?nitrogen and halogen?carbene interactions in NCX?NCX?CH₂ complexes, where X?=?F, Cl, Br and I. Molecular geometries and interaction energies of dyads and triads are investigated at the MP2/aug-cc-pVTZ level of theory. It is found that the X?N and X?C_carbene interaction energies in the triads are larger than those in the dyads, indicating that both the halogen bonding interactions are enhanced. The estimated values of cooperative energy E _coop are all negative with much larger E _coop in absolute value for the systems including iodine. The nature of halogen bond interactions of the complexes is analyzed using parameters derived from the quantum theory atoms in molecules methodology and energy decomposition analysis.

Microsolvation and hydration enthalpies of CaC2O4(H2O) n (n = 0-16) and C2O4 2-(H2O) n (n = 0-14): an ab initio study

We studied hydrated calcium oxalate and its ions at the restricted Hartree–Fock RHF/6-31G* level of theory. Performing a configurational search seems to improve the fit of the HF/6-31G* level to experimental data. The first solvation shell of calcium oxalate contains 13 water molecules, while the first solvation shell of oxalate ion is formed by 14 water molecules. The first solvation shell of Ca(II) is formed by six water molecules, while the second shell contains five. At 298.15 K, we estimate the asymptotic limits (infinite dilution) of the total standard enthalpies of hydration for Ca(II), oxalate ion and calcium oxalate as ?480.78, –302.78 and –312.73 kcal mol^?1, resp. The dissociation of hydrated calcium oxalate is an endothermic process with an asymptotic limit of +470.84 kcal mol^?1.

A B3LYP and MP2(full) theoretical investigation into the cooperativity effect between dihydrogen-bonding and H–M∙∙∙π (M = Li, Na, K) interactions among HF, MH with the π-electron donor C2H2, C2H4 or C6H6

The DFT-B3LYP/6-311++G(3df,2p) and MP2(full)/6-311++G(3df,2p) calculations were carried out on the binary complex formed by HM (M?=?Li, Na, K) and HF or the π-electron donor (C₂H₂, C₂H₄, C₆H₆), as well as the ternary system FH???HM???C₂H₂/C₂H₄/C₆H₆. The cooperativity effect between the dihydrogen-bonding and H–M???π interactions was investigated. The result shows that the equilibrium distances R _H???H and R _M???π in the ternary complex decrease and both the H???H and H–M???π interactions are strengthened when compared to the corresponding binary complex. The cooperativity effect of the dihydrogen bond on the H–M???π interaction is more pronounced than that of the M???π bond on the H???H interaction. Furthermore, the values of cooperativity effect follow the order of FH???HNa???π?>?FH???HLi???π?>?FH???HK???π and FH???HM???C₆H₆?>?FH???HM???C₂H₄?>?FH???HM???C₂H₂. The nature of the cooperativity effect was revealed by the analyses of the charge of the hydrogen atoms in H???H moiety, atom in molecule (AIM) and electron density shifts methods.