Small fullerenes with BN belts: a density functional theory investigation

Density functional calculations have been carried out on a series of BCN hybrid fullerenes with certain substitution patterns in comparison with their parent compounds Cn (n = 30, 32, 36, 38, 40, 44, 48, 50, 52). The substitutional structures, energy gaps between the highest occupied molecular orbital and the lowest unoccupied molecular orbital, ionization potentials, electron affinities, as well as molecular electrostatic potentials have been systematically investigated. The following important points of BCN hybrid fullerenes are stressed: The present studied fullerenes, comprising tubular "belt" and polar "cap", could be divided into three types of structure; each has different indexes of tubular structure and terminal caps. The properties of BCN hybrid fullerenes depend on the type of "tubular belt + polar cap" structures, especially, the HOMO and LUMO characters and MEPs of BCN fullerene are strongly governed by their structure types.     

Structure, spectroscopy, and reactivity properties of porphyrin pincers: a conceptual density functional theory and time-dependent density functional theory study

Porphyrin and pincer complexes are both important categories of compounds in biological and catalytic systems. The idea to combine them is computationally investigated in this work. By employment of density functional theory (DFT), conceptual DFT, and time-dependent DFT approaches, structure, spectroscopy, and reactivity properties of porphyrin pincers are systematically studied for a selection of divalent metal ions. We found that the porphyrin pincers are structurally and spectroscopically different from their precursors and are more reactive in electrophilic and nucleophilic reactions. A few quantitative linear/exponential relationships have been discovered between bonding interactions, charge distributions, and DFT chemical reactivity indices. These results are implicative in chemical modification of hemoproteins and understanding chemical reactivity in heme-containing and other biologically important complexes and cofactors.     

Structure,epr parameters,and reactivity of organic free radicals from a density functional approach

Summary A recently introduced density functional incorporating gradient corrections and some Hartree-Fock exchange has been used to study the structures, properties, and reactivity of representative organic free radicals. A general theoretical model has been introduced, in which standardized grid, functional, and orbital basis set are used to compute geometrical parameters, vibrational frequencies, and one-electron properties. The results are compared with available experimental data from diatomic to polyatomic radicals. All the geometric and electronic parameters compare favourably with available experimental data and with the results of refined post Hartree-Fock computations. Also the thermodynamics and kinetics of a representative unimolecular reaction (isomerization of formaldehyde radical cation) are well reproduced. These findings together with the very favourable scaling of the computations with the number of electrons suggest that the density functional approach is a promising theoretical tool for the study of relationships between structure and properties of large free radicals.     

Broken-symmetry density functional theory investigation on bis-nitronyl nitroxide diradicals: influence of length and aromaticity of couplers

A series of nitronyl nitroxide (NN) diradicals with linear conjugated couplers and another series with aromatic couplers have been investigated by the broken-symmetry (BS) DFT approach. The overlap integral between the magnetically active orbitals in the BS state has been explicitly computed and used for the evaluation of the magnetic exchange coupling constant (J). The calculated J values are in very good agreement with the observed values in the literature. The magnitude of J depends on the length of the coupler as well as the conformation of the radical units. The aromaticity of the spacer decreases the strength of the exchange coupling constant. The SOMO-SOMO energy splitting analysis, where SOMO stands for the singly occupied molecular orbital, and the calculation of electron paramagnetic resonance (EPR) parameters have also been carried out. The computed hyperfine coupling constants support the intramolecular magnetic interactions. The nature of magnetic exchange coupling constant can also be predicted from the shape of the SOMOs as well as the spin alternation rule in the unrestricted Hartree-Fock (UHF) treatment. It is found that pi-conjugation along with the spin-polarization plays the major role in controlling the magnitude and sign of the coupling constant.     

Electronic structure and reactivity of cobalt oxide dimers and their hexacarbonyl complexes: a density functional study

The dimers of cobalt oxide (CoO)(2) with cyclic and open bent structure are studied with the B1LYP density functional; the ordering of states is validated by the CCSD(T) method. The D(2h)-symmetry rhombic dioxide Co(2)O(2) with antiferromagnetically ordered electrons on cobalt centers is the global minimum. The cyclic peroxide Co(2)(O(2)) with side-on-bonded dioxygen in (7)B(2) ground state is separated from the global minimum by an energy gap of 3.15 eV. The dioxide is highly reactive as indicated by the high value of proton affinity and chemical reactivity indices. The four-member ring structures are more stable than those with three-member ring or chain configuration. The thermodynamic stability toward dissociation to CoO increases upon carbonylation, whereas proton affinity and reactivity with release of molecular oxygen also increase. The global minimum of Co(2)O(2)(CO)(6) corresponds to a triplet state (3)A" with oxygen atoms shifted above the molecular plane of the rhombic dioxide Co(2)O(2). The SOMO-LUMO gap in the ground-state carbonylated dioxide is wider, compared to the same gap in the bare dicobalt dioxide. The peroxo-isomer Co(2)(O(2))(CO)(6) retains the planar Co(2)(O(2)) ring and is only stable in a high-spin state (7)A". The carbonylated clusters have increased reactivity in both redox and nucleophilic reactions, as a result of the increased electron density in the Co(2)O(2)-ring area.     

Size dependences of the diradical character and the second hyperpolarizabilities in dicyclopenta-fused acenes: relationships with their aromaticity/antiaromaticity

Using long-range corrected density functional theory, the relationships between the electronic, magnetic, and nonlinear optical properties are drawn for two families of organic compounds, the dicyclopenta-fused acenes (DPAs) and the polyacenes (PAs), containing up to N = 12 fused rings. First, the longitudinal second hyperpolarizability (γ) of singlet DPAs is significantly enhanced with increasing system size, in comparison to PAs. This behavior is associated with an increase in the longitudinal spin polarization between the terminal five-membered rings of DPAs and is consistent with previous studies where γ is maximized for intermediate diradical character. The size dependence of the diradical character is also found to cause a hump in the γ/N evolution for singlet DPAs around N = 8. In fact, in the case of singlet PAs, the diradical characters y(0) and y(1), the various magnetic properties and the γ/N values vary monotonically with N, whereas for singlet DPAs, the shielding, the magnetizability, and the γ/N values exhibit extrema near N = 8 due to the appearance of transversal spin polarization in the middle six-membered rings in addition to the longitudinal spin polarization between the terminal five-membered rings. Moreover, it is shown that for singlet DPAs the longitudinal spin polarization (characterized by y(0)) is associated with the antiaromaticity (N ≤ 3) and the slight- or non-aromaticity (N ≥ 4) of the terminal five-membered rings, whereas the appearance of transversal spin-polarization (characterized by y(1)) is associated with the decrease in the aromaticity in the inner six-membered rings as shown for large PAs. Therefore, the exceptional behaviors in singlet DPAs for small N (N < 9) are caused by the increase in diradical character y(0) correlated with the anti-aromaticity or the slight-/non-aromaticity of terminal rings and the corresponding emergence of a global aromatic character. Such a relationship between the aromaticity/antiaromaticity and the diradical character is useful for designing real open-shell NLO molecules through the control of their diradical characters.     

Stability, reactivity, and aromaticity of compounds of a multivalent superatom

In this article, we analyze the stability, reactivity, and possible aromatic behavior of two recently reported clusters (Reveles, J. U.; Khanna, S. N.; Roach, P. J.; Castleman, A. W., Jr. Proc. Natl. Acad. Sci. 2006, 103, 18405), viz., Al(7)C(-) and Al(7)O(-) in the light of the principles of the maximum hardness and minimum electrophilicity as well as the nucleus-independent chemical shift values. Stability of these clusters in the context of addition/removal of an electron or an Al atom is now clearly understood.     

Electronic and magnetic properties of substituted BN sheets: a density functional theory study

Using density functional calculations, we investigate the geometries, electronic structures and magnetic properties of hexagonal BN sheets with 3d transition metal (TM) and nonmetal atoms embedded in three types of vacancies: V(B), V(N), and V(B+N). We show that some embedded configurations, except TM atoms in V(N) vacancy, are stable in BN sheets and yield interesting phenomena. For instance, the band gaps and magnetic moments of BN sheets can be tuned depending on the embedded dopant species and vacancy type. In particular, embedment such as Cr in V(B+N), Co in V(B), and Ni in V(B) leads to half-metallic BN sheets interesting for spin filter applications. From the investigation of Mn-chain (C(Mn)) embedments, a regular 1D structure can be formed in BN sheets as an electron waveguide, a metal nanometre wire with a single atom thickness.     

Nitrobenzoxadiazoles and related heterocycles: a relationship between aromaticity, superelectrophilicity and pericyclic reactivity

A study of the dual electrophilic and pericyclic reactivity of 4,6-dinitrobenzofurazan (DNBZ, 2), 4,6-dinitro-2,1,3-benzothiadiazole (DNBS, 3), 4,6-dinitro-2,1,3-benzoselenadiazole (DNBSe, 4) is reported. Kinetic and thermodynamic measurements of the ease of covalent hydration of 2-4 to give the corresponding hydroxy sigma-adducts C-2-C-4 have been carried out over a large pH range in aqueous solution. Analysis of the data has allowed a determination of the rate constants k1(H2O) pertaining to the susceptibility of 2-4 to water attack as well as the pKa values for the sigma-complexation processes. With pKa values ranging from 3.92 for DNBZ to 6.34 for DNBSe to 7.86 for DNBS, the electrophilic character of the three heteroaromatics is much closer to that of the superelectrophilic reference, i.e. 4,6-dinitrobenzofuroxan (DNBF, 1; pKa = 3.75), than that of the standard Meisenheimer electrophile 1,3,5-trinitrobenzene (TNB, pKa = 13.43). Most importantly, water is found to be an efficient nucleophile which contributes strongly to the formation of the adducts C-2 and C-4. This confirms a previous observation that a pKa value of ca. 8 is a primary requirement for having H2O competing effectively as a nucleophile with OH- in the formation of hydroxy sigma-adducts. On the other hand, 2-4 are found to exhibit dienophilic and/or heterodienic behaviour on treatment with isoprene, 2,3-dimethylbutadiene, cyclopentadiene or cyclohexadiene, affording Diels-Alder mono- or di-adducts which have all been structurally characterized. A major finding is that the order of Diels-Alder reactivity follows clearly the order of electrophilicity, pointing to a direct relationship between superelectrophilic and pericyclic reactivity. This relationship is discussed.     

Cobalt analogues of Roussin's red salt esters: a density functional study

Density functional theory calculations on the Co(2)(NO)(4)(SR)(2) compounds (R = CH(3), CF(3) and C(4)H(9)) predict butterfly and open isomeric structures with and without a direct Co-Co bond. The open Co(2)(NO)(4)(SR)(2) structures are favored over the butterfly isomers, in terms of relative energy. Furthermore the open structures are predicted to have approximately twice as large HOMO-LUMO gaps than the butterfly Co(2)(NO)(4)(SR)(2) isomers. For the related Co(2)(CO)(6)(SR)(2) species, competing open and butterfly structures with similar HOMO-LUMO gaps were predicted. This could explain why the Co(2)(NO)(4)(μ-SR)(2) compounds have already been synthesized and why no genuine Co(2)(CO)(6)(SR)(2) derivatives have yet been reported.     

Pentapnictogen heterocyclic monoanions: Structure,stability, and aromaticity

Inorganic planar ring-shape molecules with 4n + 2 π electrons are always the focus of experimental synthesis and theoretical research due to their potential aromaticity and stability. In this work, the whole series of five-membered heterocycle monoanions X _nY_5-n⁻ (X, Y = group 15 elements; n = 1-4) were thoroughly investigated by means of density functional theory calculations. They all have large formation energies and HOMO-LUMO gap energies, suggesting the potential thermodynamic and kinetic stability. Their aromaticities are comparable to that of typical aromatic hydrocarbons. Their thermal stabilities were firmly established by the ab initio molecular dynamics simulations. As most of them are predicted for the first time, their various spectra were simulated for experimental characterization. Furthermore, we demonstrate that these five-membered cyclic anions can be employed as η⁵-ligand to construct novel all-inorganic metallocenes, which may serve as the building blocks of low-dimensional nanomaterials.     

Electronic structures, intramolecular interactions, and aromaticity of substituted 1-(2-iminoethylidene) silan amine: a density functional study

The intramolecular hydrogen bond, molecular structure, π electrons delocalization, and vibrational frequencies in 1-(2-iminoethylidene) silan amine and its derivatives have been investigated by means of density functional method with 6-311++G** basis set, in gas phase, water, and carbon tetrachloride solutions. The obtained results showed that the hydrogen bond strength is mainly governed by resonance variations inside the chelate ring induced by the substituent groups. Furthermore, the topological properties of the electron density distributions for N–H···N intramolecular hydrogen bond were analyzed in terms of the Bader's theory of atoms in molecules. On the other hand, the aromaticity of the ring formed is measured using several well-established indices of aromaticity such as nucleus-independent chemical shift, harmonic oscillator models of the aromaticity, para-delocalization index, average two-center indices, aromatic fluctuation index, and π-fluctuation aromatic index. Natural population analysis data, the electron density and Laplacian properties, as well as γ(NH) and ν(NH) were further used for estimation of the hydrogen bonding interactions and the forces driving their formation.     

Structure of inhomogeneous polymer solutions: a density functional approach

The structure of polymer solutions confined between surfaces is studied using a density functional theory where the polymer molecules have been modeled as a pearl necklace of freely jointed hard spheres and the solvent as hard spheres. The present theory uses the concept of universality of the free energy density functional to obtain the first-order direct correlation function of the nonuniform system from that of the corresponding uniform system, calculated through the Verlet-modified type bridge function. The uniform bulk fluid direct correlation function required as input has been calculated from the reference interaction site model integral equation theory using the Percus-Yevick closure relation. The calculated results on the density profiles of the polymer as well as the solvent are shown to compare well with computer simulation results.     

Approximation of the molecular electrostatic potential in a gaussian density functional method

The recently developed Asymptotic Density Model (ADM) [6, 9] is here implemented in the density functional framework using the program deMon-KS [13]. While the original implementation divided the atoms into a core shell and a valence shell, the present version allows for an arbitrary number of shells making it therefore more flexible and, as shown with benzene, potentially more accurate. Moreover, since this method is derived through Poisson's equation, an expression for the electronic charge density is also obtained. However, the present discussion will restrict itself to the electrostatic potential. Finally, even though this method requires parametrization, it is shown that the parameters obtained for homonuclear diatomic species, and used as is in molecular calculations, yield satisfactory results. Indeed, the ADM reproduces almost all basic features of the MEP for all molecules presented here, (water, ammonia, ethylene, acetylene, hydrogen cyanide, carbon monoxide, benzene, nitrous acid). Received: 5 July 1996 / Accepted: 12 November, 1996     

Dispersion, static correlation, and delocalisation errors in density functional theory: an electrostatic theorem perspective

Dispersion, static correlation, and delocalisation errors in density functional theory are considered from the unconventional perspective of the force on a nucleus in a stretched diatomic molecule. The electrostatic theorem of Feynman is used to relate errors in the forces to errors in the electron density distortions, which in turn are related to erroneous terms in the Kohn-Sham equations. For H(2), the exact dispersion force arises from a subtle density distortion; the static correlation error leads to an overestimated force due to an exaggerated distortion. For H(2)(+), the exact force arises from a delicate balance between attractive and repulsive components; the delocalisation error leads to an underestimated force due to an underestimated distortion. The net force in H(2)(+) can become repulsive, giving the characteristic barrier in the potential energy curve. Increasing the fraction of long-range exact orbital exchange increases the distortion, reducing delocalisation error but increasing static correlation error.     

Structure and ionization energies of some analogues of iron-only hydrogenases studied by density functional theory methods

Density functional theory calculations using both the B3LYP and BP86 functional in conjunction with a medium and large size basis set have been used to predict the structures and ionization energies of 12 models of iron-only hydrogenases. Although the structural predictions do not allow a clear discrimination between the different computational models, these models do yield significantly different adiabatic and vertical ionization energies. The closest agreement with experiment is given by the BP86 functional and the large all-electron basis. At this level of theory the adiabatic ionization energies are very close to experiment, but the vertical values are uniformly too small, leading to an underestimation of the reorganization energies. The calculations also suggest that measured ionization energies may help in identifying both the bridge-head group and whether CO bridging takes place upon ionization.