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.

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 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.

Ab initio molecular dynamics simulation on the formation process of He@C60 synthesized by explosion

The applications of endohedral non-metallic fullerenes are limited by their low production rate. Recently, an explosive method developed in our group shows promise to prepare He@C₆₀ at fairly high yield, but the mechanism of He inserting into C₆₀ cage at explosive conditions was not clear. Here, ab initio molecular dynamics analysis has been used to simulate the collision between C₆₀ molecules at high-temperature and high-pressure induced by explosion. The results show that defects formed on the fullerene cage by collidsion can effectively decrease the reaction barrier for the insertion of He into C₆₀, and the self-healing capability of the defects was also observed.

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.

Interactions of calcium with the external surfaces of fullerenes and endofullerenes doped with radioactive sodium iodide

We report first-principles calculations carried out to analyze the adsorption of calcium on the outer surface of the fullerene C₆₀, yielding [C₆₀?+?mCa]. Geometric optimization (GO) and molecular dynamics (MD) simulation were performed using the plane-wave pseudopotential method within the framework of density functional theory (DFT) and time-dependent DFT (TD-DFT) to investigate the configurations, the associated energies in the ground state, and the stabilities of fullerenes and endofullerenes doped with radioactive sodium iodide when they interact with calcium atoms on the outer fullerene surface (i.e., [nNa¹³¹I@C₆₀?+?mCa]). The reason for investigating these calcium-functionalized (endo)fullerene systems was to gauge their potential stability when used as vectors to deliver radioiodine to cancerous tissue in the human body. In the simulations, we found that the geometric limit on the number of calcium atoms that can be physisorbed on the outer surface of an empty fullerene while maintaining its structural stability is 28 calcium atoms, which also takes into account the proportional expansion of the fullerene as the number of absorbed calcium atoms increases. However, the stability of a fullerene system during calcium adsorption also strongly depends on whether any atoms or molecules are being encapsulated by the fullerene, as these encapsulated atoms/molecules can also interact with the fullerene and influence its stability. A Mulliken electronegativity analysis revealed that, when atoms inside and/or outside the fullerene donate charge (electrons) to the fullerene, the fullerene expands. The excess charge on the carbon atoms of the fullerene weakens some of the carbon–carbon bonds, potentially causing them to break, in which case the fullerene loses its ability to encapsulate molecules and releases them.

Open image in new window

Graphical Abstract DFT simulation of a endo fullerene doped with radioactive sodium iodide interacting with 28 calcium atoms in a geometric arrangement

     

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 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).

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.

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.

Relation between electron exchange and nuclear vibrations in the H 2 + +H 2 + → H 3 + +p exothermic ion-molecular reaction

The effective cross section for the H ₂ ⁺ +H ₂ ⁺ → H ₃ ⁺ +p reaction in the energy range 5.7–11.5 eV is measured by the split beam method. The maximum of the cross section at an energy of ～8 eV is related to the production of the H ₄ ⁺⁺ compound system. The reaction threshold W _thr≈5 eV provides evidence in favor of the classical model of the H ₂ ⁺ ion with the charge fixed on one of the nuclei throughout the collision event.     

Novel open-framework architecture in strontium(II) phosphonate

Hydrothermal reaction of strontium nitrate with N,N′-piperazine-bis(methylenephosphonic acid) (H₄L) afforded a novel three-dimensional strontium(II) diphosphonate, namely, [Sr(H₂L)(H₂O)₂] · 3H₂O (1). The structure of compound 1 can be viewed as a pillar-layered architecture in which 1D chains of [Sr(PO₃)] along a-axes are further crosslinked by the piperazine bismethylene spacers along b- and c-axes.     

Synthesis of trialkyl[60]fullerene C60(CH2SiMe3)3H and its potassium and rhodium(I) complexes

The reaction of [60]fullerene with Me₃SiCH₂MgCl in 1,2-Cl₂C₆H₄/THF (1/1) under dry air afforded a bis-alkyl adduct, C₆₀(CH₂SiMe₃)₂ (1), in 54% yield. Treatment of 1 with Me₃SiCH₂MgCl in THF under argon then afforded a trialkyl[60]fullerene, C₆₀(CH₂SiMe₃)₃H (2), in 37% yield. Further treatment of 2 with KO^tBu gave a potassium salt, [K(thf)_n][C₆₀(CH₂SiMe₃)₃] (3), which was then converted to a C_s-symmetric Rh(I) complex, Rh[η⁵-C₆₀(CH₂SiMe₃)₃](1,5-cyclooctadiene) (4), in 91% yield.     

s-Block metallabenzene: aromaticity and hydrogen adsorption

Second group metal dimers can replace the carbon atom in benzene to form metallabenzene (C₅H₆M₂) compounds. These complexes possess some aromatic character and promising hydrogen adsorption properties. In this study, we investigated the aromatic character of these compounds using aromaticity indices and molecular orbital analysis. To determine the nature of interactions between hydrogen and the metallic center, variation-perturbational decomposition of interaction energy was applied together with ETS-NOCV analysis. The results obtained suggest that the aromatic character comes from three π orbitals located mainly on the C₅H₅ ^? fragment. The high hydrogen adsorption energy (up to 6.5 kcal mol^?1) results from two types of interaction. In C₅H₆Be₂, adsorption is controlled by interactions between the empty metal orbital and the σ orbital of the hydrogen molecule (Kubas interaction) together with corresponding back-donation interactions. Other C₅H₆M₂ compounds adsorb H₂ due to Kubas interactions enhanced by H₂–π interactions.

Open image in new window

Graphical Abstract First π orbital in C₅H₆Be₂

     

New routes toward metallated methanofullerene terdentate bisaminoaryl ligands

The preparation of a series of new ligands, comprising a NCN ligand (C₆H₂(CH₂NMe₂)₂-2,6-X-4, X=H, Br, I) substituted on the para position with a [1,2]-methanofullerene group is described. Metallation of these ligands is achieved via an oxidative addition reaction of nickel(0) or palladium(0) complexes. The first step of this reaction comprises the formation of a π-complex between the metal center and the fullerene moiety of the ligand, which upon heating in a subsequent step is followed by oxidative addition of the NCN-halide bond to the metal center. The nickel(II) complex C₆₀C(Me)NCN·NiBr (9) was isolated in pure form. Low temperature NMR spectroscopy and theoretical calculations showed that the NCN moiety in 9 rotates about the (Me)CNCN bond at room temperature. UV-Vis spectroscopic measurements of 9 suggest that binding a metal center via a covalent metal-carbon bond to the NCN moiety does not lead to a substantial change in the electronic structure of the fullerene moiety. The cyclo voltammograms of 9 showed the absence of a Ni(II)/Ni(III) oxidation/reduction process, while the second and fourth C₆₀ reduction potentials were anodically shifted.     

Electron correlation effects and density analysis of the first-order hyperpolarizability of neutral guanine tautomers

Dipole moments (μ), charge distributions, and static electronic first-order hyperpolarizabilities (β _μ ) of the two lowest-energy keto tautomers of guanine (7H and 9H) were determined in the gas phase using Hartree–Fock, Møller–Plesset perturbation theory (MP2 and MP4), and DFT (PBE1PBE, B97-1, B3LYP, CAM-B3LYP) methods with Dunning’s correlation-consistent aug-cc-pVDZ and d-aug-cc-pVDZ basis sets. The most stable isomer 7H exhibits a μ value smaller than that of the 9H form by a factor of ca. 3.5. The β _μ value of the 9H tautomer is strongly dependent on the computational method employed, as it dramatically influences the β _μ (9H)/β _μ (7H) ratio, which at the highest correlated MP4/aug-cc-pVDZ level is predicted to be ca. 5. The Coulomb-attenuating hybrid exchange-correlation CAM-B3LYP method is superior to the conventional PBE1PBE, B3LYP, and B97-1 functionals in predicting the β _μ values. Differences between the largest diagonal hyperpolarizability components were clarified through hyperpolarizability density analyses. Dipole moment and first-order hyperpolarizability are molecular properties that are potentially useful for distinguishing the 7H from the 9H tautomer.

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.

CO2 adsorption on polar surfaces of ZnO

Physical and chemical adsorption of CO₂ on ZnO surfaces were studied by means of two different implementations of periodic density functional theory. Adsorption energies were computed and compared to values in the literature. In particular, it was found that the calculated equilibrium structure and internuclear distances are in agreement with previous work. CO₂ adsorption was analyzed by inspection of the density of states and electron localization function. Valence bands, band gap and final states of adsorbed CO₂ were investigated and the effect of atomic displacements analyzed. The partial density of states (PDOS) of chemical adsorption of CO₂ on the ZnO(0001) surface show that the p orbitals of CO₂ were mixed with the ZnO valence band state appearing at the top of the valence band and in regions of low-energy conduction band.

On the identification of ionic species of neutral halogen dimers, monomers and pincer type palladacycles in solution by electrospray mass and tandem mass spectrometry

Electrospray (ESI) mass spectra analysis of acetonitrile solutions of a series of neutral chloro dimers, pincer type, and monomeric palladacycles has enabled the detection of several of their derived ionic species. The monometallic cationic complexes Pd[κ¹-C,κ¹-N,κ¹-S-C(CH₃S-2-C₆H₄)C(Cl)CH₂N(CH₃)₂]⁺ (1a) and [Pd[κ¹-C,κ¹-N,κ¹-S-C(CH₃S-2-C₆H₄)C(Cl)CH₂N(CH₃)₂](CH₃CN)]⁺ (1b) and the bimetallic cationic complex [κ¹-C,κ¹-N,κ¹-S-C(CH₃S-2-C₆H₄)C(Cl)CH₂N(CH₃)₂]Pd-Cl-Pd[κ¹-C,κ¹-N,κ¹-S-C(CH₃S-2-C₆H₄)C(Cl)CH₂N(CH₃)₂]⁺ (1c) were detected from an acetonitrile solution of the pincer palladacycles Pd[κ¹-C,κ¹-N,κ¹-S-C(CH₃S-2-C₆H₄)C(Cl)CH₂N(CH₃)₂](Cl) 1. For the dimeric compounds {Pd[κ¹-C,κ¹-N-C(Y-2-C₆H₄)C(Cl)CH₂N(CH₃)₂](μ-Cl)}₂ (2, Y=H and 3, CF₃), highly electronically unsaturated palladacycles [Pd[κ¹-C,κ¹-N-C(Y-2-C₆H₄)C(Cl)CH₂N(CH₃)₂]⁺ (2d, 3d) and their mono and di-acetonitrile adducts, namely, [Pd[κ¹-C,κ¹-N-C(Y-2-C₆H₄)C(Cl)CH₂N(CH₃)₂](CH₃CN)]⁺ (2e, 3e) and [Pd[κ¹-C,κ¹-N-C(Y-2-C₆H₄)C(Cl)CH₂N(CH₃)₂](CH₃CN)₂]⁺ (2f and 3f) were detected together with the bimetallic complex [Pd[κ¹-C,κ¹-N-C(Y-2-C₆H₄)C(Cl)CH₂N(CH₃)₂]-Cl-Pd[κ¹-C,κ¹-N-C(Y-2-C₆H₄)C(Cl)CH₂N](CH₃)₂]⁺ (2a, 3a) and its acetonitrile adducts [κ¹-C,κ¹-N-C(Y-2-C₆H₄)C(Cl)CH₂N(CH₃)₂](CH₃CN)Pd-Cl-Pd[ κ¹-C,κ¹-N-C(Y-2-C₆H₄)C(Cl)CH₂N(CH₃)₂]⁺ (2b, 3b) and [κ¹-C,κ¹-N-C(Y-2-C₆H₄)C(Cl)CH₂N(CH₃)₂](CH₃CN)Pd-Cl-Pd[κ¹-C, κ¹-N-C(Y-2-C₆H₄)C(Cl)CH₂N(CH₃)₂(CH₃CN)]⁺ (2c, 3c). The dimeric palladacycle {Pd[κ¹-C,κ¹-N-C(CH₃O-2-C₆H₄)C(Cl)CH₂N(CH₃)₂](μ-Cl)}₂ (4) is unique as it behaves as a pincer type compound with the OCH₃ substituent acting as an intramolecular coordinating group which prevents acetonitrile full coordination, thus forming the cationic complexes [(C₆H₄(o-CH₃O)CC(Cl)CH₂N(CH₃)₂-κO,κC,κN)Pd]⁺ (4b), [(C₆H₄(o-CH₃O)CC(Cl)CH₂N(CH₃)₂- κO,κC,κN)Pd(CH₃CN)]⁺ (4c) and [(C₆H₄ (o-MeO)CC(Cl)CH₂N(CH₃)₂-κO, κC,κN)Pd-Cl-Pd(C₆H₄(o-CH₃O)CC(Cl)CH₂N(CH₃)₂-κO,κC,κN)]⁺ (4a). ESI-MS spectra analysis of acetonitrile solutions of the monomeric palladacycles Pd[κ¹-C,κ¹-N-C(Y-2-C₆H₄)C(Cl)CH₂N(CH₃)₂](Cl)(Py) (5, Y=H and 6, Y=CF₃) allows the detection of some of the same species observed in the spectra of the dimeric palladacycles, i.e., monometallic cationic 2d-3d, 2e-3e and {Pd[κ¹-C,κ¹-N-C(Y-2-C₆H₄)C(Cl)CH₂N(CH₃)₂](Py)}⁺ (5a, 6a) and {Pd[κ¹-C,κ¹-N-C(Y-2-C₆H₄)C(Cl)CH₂N(CH₃)₂](CH₃CN)(Py)}⁺ (5b, 6b) and the bimetallic 2a, 3a, 2b, 3b, 2c and 3c. In all cationic complexes detected by ESI-MS, the cyclometallated moiety was intact indicating the high stability of the four or six electron anionic chelate ligands. The anionic (chloride) or neutral (pyridine) ligands are, however, easily replaced by the acetonitrile solvent.