Conformations and electronic structures of axially coordinated fullerene-porphyrin-fullerene triad (C60CHCOO)2-Sn(IV) porphyrin

The conformational (cis and trans) stability and electronic structures of (C(60)CHCOO)(2)-Sn(IV) porphyrin, recently synthesized as a novel fullerene-porphyrin-fullerene triad linked by metal axial coordination, have been studied by ab initio calculations. The cis conformer was found to be slightly more stable than the trans by 1.38 kcal/mol in the neutral compound. Upon the addition of an electron to the triad, the relative stability of the cis conformer was found to be higher (3.29 kcal/mol) than that in the neutral one. From the investigation of frontier molecular orbitals, for the cis conformer, it was found that the electrons are localized in HOMO of the porphyrin, while the electrons are localized in LUMO of the syn-fullerene. For the trans conformer, it was found that the electrons are localized in HOMO of the porphyrin, while the electrons are localized in LUMO of one of the two fullerene moieties, and the electrons are localized in LUMO2 of the other fullerene moiety, but the LUMO and LUMO2 have the same orbital energy. Thus, the PET may take place unidirectionally in the cis conformer from the porphyrin to the syn-fullerene, while it is bidirectional from the porphyrin to both of the fullerene moieties.     

金属酞菁配合物的从头计算研究

利用从头计算方法研究了6种金属酞菁MPc(M=Zn2+、Cu2 +、Ni2+、Co2+、Fe2+、Mn2+)。得到了它们的基态能量,基态自旋多重度,分子轨道组成与能级,电荷分布与键序。其中,自旋多重度的计算结果与实验相符。中心离子d轨道参与HOMO、LUMO构成的程度可以解释各MPc光敏活性不同的实验现象。     

Zeolite-Y entrapped bivalent transition metal complexes as hybrid nanocatalysts: density functional theory investigation and catalytic aspects

The intriguing research toward the exploitation of zeolite-Y-based hybrid nanocatalysts for catalytic oxidation reactions has been growing significantly. In the present investigation, we describe the synthesis of zeolite-Y entrapped transition metal complexes of the general formulae [M(SFCH)·xH₂O]-Y (where, M = Mn, Fe, Co, Ni (x = 3) and Cu (x = 1)); H₂SFCH = (E)-N′-(2-hydroxybenzylidene)furan-2-carbohydrazide]. These nanocatalysts have been characterized by various physicochemical techniques. Density functional theory calculations are performed to address the relaxed geometry, bond angle, bond length, dihedral angle, highest occupied molecular orbital (HOMO)–lowest unoccupied molecular orbital (LUMO) energy gap, and electronic density of states of H₂SFCH ligand and their neat transition metal complexes. The observed HOMO–LUMO gap and the Fermi energy is higher for Cu(II) complexes, which demonstrates the better catalytic activity of this nanocatalyst. The catalytic activity was performed in liquid-phase oxidation of cyclohexane using hydrogen peroxide as oxidant to give cyclohexanone (CyONE) and cyclohexanol (CyOL). Among them, [Cu(SFCH)·H₂O]-Y catalyst has the highest selectivity toward CyONE (84.5%).     

The role of chelating ligands and central metals in the oxygen reduction reaction activity: a DFT study

A detailed investigation of oxygen reduction reaction (ORR) catalyzed by various metal chelates has been performed by DFT study. The results indicate that the ORR activity is determined by both of the central metal ions and chelating ligands, among which the former play a key role. For the same ligand, the central metal ions Fe, Co, or Mn give higher ORR activity, while the others almost have no catalytic activity, which is due to the fact that the O₂ and oxygen containing species are either excessively adsorbed (on central Cr) or difficult to be adsorbed on the active sites (for central Zn, Cu, or Ni). Furthermore, the ORR activity for Fe chelates is slightly increased with the increase of ligand field strength, while for other metal chelates there seems to be no clear trends between ligand field strength and ORR activity. The origin of the ORR activity for the studied metal chelates is mainly attributed to the appropriate energy gaps between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO).     

Molecular details from computational reaction dynamics for the cellobiohydrolase I glycosylation reaction

Glycosylation of cellobiose hydrolase I (CBHI), is a key step in the processing and degradation of cellulose. Here the pathways and barriers of the reaction are explored using the free energy from adaptive reaction coordinate forces (FEARCF) reaction dynamics method coupled with SCC-DFTB/MM. In many respects CBHI follows the expected general GH7 family mechanism that involves the Glu-X-Asp-X-X-Glu motif. However, critical electronic and conformational details, previously not known, were discovered through our computations. The central feature that ensures the success of the glycosylation reaction are the Glu212 nucleophile's hydrogen bond to the hydroxyl on C2, of the glucose in the -1 position of the cellulosic strand. This Glu212 function restricts the C2 hydroxyl in such a way as to favor the formation of the (4)E ring pucker of the -1 position glucose. A frontier molecular orbital analysis of the structures along the reaction surface proves the existence of an oxocarbenium ion, which has both transition state and intermediate character. The transition state structure is able to descend down the glycosylation pathway through the critical combination of Asp214 (HOMO), ring oxygen (LUMO), and Glu212 (HOMO), anomeric carbon (LUMO) interactions. Using the fully converged FEARCF SCC-DFTB/MM reaction surface, we find a barrier of 17.48 kcal/mol separating bound cellulose chain from the glycosylated CBHI. Taking recrossing into account gives k(cat) = 0.415 s(-1) for cellobiohydrolase glycosylation.     

Elaborated spectral,modeling, QSAR,docking, thermal,antimicrobial and anticancer activity studies for new nanosized metal ion complexes derived from sulfamerazine azodye

New azodye ligand (H₂L) and its relative Cr(III)-, Mn(II)-, Fe(III)-, Co(II)-, Ni(II)-, Cu(II)-, Zn(II)- and Cd(II)-nanosized complexes were prepared. A new synthesized compounds were characterized using spectral (mass, IR, UV–Vis, XRD, and ESR) and analytical (elemental, molar conductance, thermal and magnetic moment measurements) tools. Infrared spectra showed that the ligand behaves as a monobasic bidentate, coordinating with central atoms through carbonyl oxygen and α-hydroxyl group. The geometrical structures of Cr(III) and Fe(III) complexes were found to be in octahedral configuration, whereas Mn(II), Co(II), Ni(II), Cu(II), Zn(II) and Cd(II) complexes have tetrahedral forms. XRD patterns reflect an amorphous appearance of all investigated complexes. TEM images showed nanosized particles and identical distribution over the complex surface. Molecular modeling for the drug ligand and its metal ion complexes were performed using Gaussian09 program to assert on their structural formulae. Some essential parameters were extracted using HOMO and LUMO energies. AutoDock tools 4.2 was used to simulate the interaction process with infected cell proteins to expect the experimental pathway. The inhibition activity of drug ligand and its metal ion complexes was evaluated towards different types of bacteria and fungi through in vitro antimicrobial activities. The antitumor activities of all compounds are straightened towards human liver carcinoma (HEPG2) cell lines. Fe(III) and Co(II) complexes exhibited IC₅₀ of 2.90 and 4.23 µg mL^?1, respectively, which means they are more potent anticancer drug than the standard (doxorubicin, IC₅₀ = 4.73 µg mL^?1). Therefore, the two complexes may consider promising anticancer drugs.

     

Fe<Subscript><Emphasis Type="Italic">n</Emphasis></Subscript> Clusters (<Emphasis Type="Italic">n</Emphasis> = 2–7) Interaction with Furan Ring: DFT Studies over Iron Surface Suitability for Furan Adsorption

The interaction of different iron clusters (Fe₂, Fe₃, Fe₄, Fe₅, Fe₆ and Fe₇) with furan compound was studied by density functional theory (DFT) and the results show that the compound possess suitable structural and electronic parameters for the metal adhesion. After analyzing the binding energy and molecular orbital studies, it is found that there is a bond between furan ring and metal. In the molecular orbital, since the HOMO localizes over furan ring, especially over C2=C3–C4=C5, the furan acts as donator (HOMO) and metal performs as acceptor (LUMO) in the interaction of furan with iron surface, transferring the high charge density mainly from the delocalization region of furan ring to the metal (L(σ) → Fe).     

Molecular and vibrational structure of tetroxo d0 metal complexes in their excited states. a study based on time-dependent density functional calculations and Franck-Condon theory

We have applied time dependent density functional theory to study excited state structures of the tetroxo d(0) transition metal complexes MnO(4)(-), TcO(4)(-), RuO(4), and OsO(4). The excited state geometry optimization was based on a newly implemented scheme [Seth et al. Theor. Chem. Acc. 2011, 129, 331]. The first excited state has a C(3v) geometry for all investigated complexes and is due to a "charge transfer" transition from the oxygen based HOMO to the metal based LUMO. The second excited state can uniformly be characterized by "charge transfer" from the oxygen HOMO-1 to the metal LUMO with a D(2d) geometry for TcO(4)(-), RuO(4), and OsO(4) and two C(2v) geometries for MnO(4)(-). It is finally found that the third excited state of MnO(4)(-) representing the HOMO to metal based LUMO+1 orbital transition has a D(2d) geometry. On the basis of the calculated excited state structures and vibrational modes, the Franck-Condon method was used to simulate the vibronic structure of the absorption spectra for the tetroxo d(0) transition metal complexes. The Franck-Condon scheme seems to reproduce the salient features of the experimental spectra as well as the simulated vibronic structure for MnO(4)(-) generated from an alternative scheme [Neugebauer J. J. Phys. Chem. A 2005, 109, 1168] that does not apply the Franck-Condon approximation.     

Exohedral functionalization of C60 by [4+2] cycloaddition of multiple anthracenes

By using density functional theory, we have investigated [4+2] cycloaddition of one to four anthracene (ANT) molecules to a C₆₀ fullerene which has been already studied by experimentalists. It was found that the reaction is regioselective and the ANT molecule prefers to be adsorbed atop a C–C bond of the C₆₀ which is shared between two hexagons with reaction energy of ?25.2 kcal/mol (for one ANT). The HOMO of the ANT interacts with the LUMO of the C₆₀ via a cycloaddition reaction. Also five regioisomeric bis-adducts of ANT/C₆₀ complexes were compared from stand point of stability. Increasing the number of attached ANTs, the reaction energy becomes less negative. The HOMO–LUMO energy gap of C₆₀ is slightly changed and the potential barrier of the field electron emission from its surface may be reduced upon the reaction.     

Monomeric MnIII/II and FeIII/II complexes with terminal hydroxo and oxo ligands: probing reactivity via O-H bond dissociation energies

Non-heme manganese and iron complexes with terminal hydroxo or oxo ligands are proposed to mediate the transfer of hydrogen atoms in metalloproteins. To investigate this process in synthetic systems, the monomeric complexes [M(III/II)H(3)1(OH)](-/2-) and [M(III)H(3)1(O)](2-) have been prepared, where M(III/II) = Mn and Fe and [H(3)1](3-) is the tripodal ligand, tris[(N'-tert-butylureaylato)-N-ethyl)]aminato. These complexes have similar primary and secondary coordination spheres, which are enforced by [H(3)1](3-). The homolytic bond dissociation energies (BDEs(O-H)) for the M(III/II)-OH complexes were determined, using experimentally obtained values for the pK(a)(M-OH) and E(1/2) measured in DMSO. This thermodynamic analysis gave BDEs(O-H) of 77(4) kcal/mol for [Mn(II)H(3)1(O-H)](2-) and 66(4) kcal/mol for [Fe(II)H(3)1(O-H)](2-). For the M(III)-OH complexes, [Mn(III)H(3)1(OH)]- and [Fe(III)H(3)1(OH)]-, BDEs(O-H) of 110(4) and 115(4) kcal/mol were obtained. These BDEs(O-H) were verified with reactivity studies with substrates having known X-H bond energies (X = C, N, O). For instance, [Fe(II)H(3)1(OH)](2-) reacts with a TEMPO radical to afford [Fe(III)H(3)1(O)](2-) and TEMPO-H in isolated yields of 60 and 75%, respectively. Consistent with the BDE(O-H) values for [Mn(II)H(3)1(OH)](2-), TEMPO does not react with this complex, yet TEMPO-H (BDE(O-H) = 70 kcal/mol) reacts with [Mn(III)H(3)1(O)](2-), forming TEMPO and [Mn(II)H(3)1(OH)](2-). [Mn(III)H(3)1(O)](2-) and [Fe(III)H(3)1(O)](2-) react with other organic substrates containing C-H bonds less than 80 kcal/mol, including 9,10-dihydroanthracene and 1,4-cyclohexadiene to produce [M(II)H(3)1(OH)](2-) and the appropriate dehydrogenated product in yields of greater than 80%. Treating [Mn(III)H(3)1(O)](2-) and [Fe(III)H(3)1(O)](2-) with phenolic compounds does not yield the product expected from hydrogen atom transfer but rather the protonated complexes, [Mn(III)H(3)1(OH)]- and [Fe(III)H(3)1(OH)]-, which is ascribed to the highly basic nature of the terminal oxo group.     

The effect and influence of cis-ligands on the electronic and oxidizing properties of nonheme oxoiron biomimetics. A density functional study

Density functional theory studies on the nature of the cis effect and cis influence of ligands on oxoiron nonheme complexes have been performed. A detailed analysis of the electronic and oxidizing properties of [Fe(IV)O(TPA)L](+) with L = F(-), Cl(-), and Br(-) and TPA = tris-(2-pyridylmethyl)amine are presented and compared with [Fe(IV)O(TPA)NCCH(3)](2+). The calculations show that the electronic cis effect is determined by favorable orbital overlap between first-row elements with the metal, which are missing between the metal and second- and third-row elements. As a consequence, the metal 3d block is split into a one-below-two set of orbitals with L = Cl(-) and Br(-), and the HOMO/LUMO energy gap is widened with respect to the system with L = F(-). However, this larger HOMO/LUMO gap does not lead to large differences in electron affinities of the complexes. Moreover, a quantum mechanical analysis of the binding of the ligand shows that it is built up from a large electric field effect of the ligand on the oxoiron species and a much smaller quantum mechanical effect due to orbital overlap. These contributions are of similar strength for the three tested halogen cis ligands and result in similar reactivity patterns with substrates. The calculations show that [Fe(IV)O(TPA)L](+) with L = F(-), Cl(-), and Br(-) have closely lying triplet and quintet spin states, but only the quintet spin state is reactive with substrates. Therefore, the efficiency of the oxidant will be determined by the triplet-quintet spin state crossing of the reaction. The reaction of styrene with a doubly charged reactant, that is, [Fe(V)O(TPA)L](2+) with L = F(-), Cl(-), and Br(-) or [Fe(V)O(TPA)NCCH(3)](3+), leads to an initial electron transfer from the substrate to the metal followed by a highly exothermic epoxidation mechanism. These reactivity differences are mainly determined by the overall charge of the system rather than the nature of the cis ligand.     

Comparative study of Li, Na, and K adsorptions on graphite by using ab initio method

A comprehensive study has been conducted to compare the adsorptions of alkali metals (including Li, Na, and K) on the basal plane of graphite by using molecular orbital theory calculations. All three metal atoms prefer to be adsorbed on the "middle hollow site" above a hexagonal aromatic ring. A novel phenomenon was observed, that is, Na, instead of Li or K, is the weakest among the three types of metal atoms in adsorption. The reason is that the SOMO (single occupied molecular orbital) of the Na atom is exactly at the middle point between the HOMO and the LUMO of the graphite layer in energy level. As a result, the SOMO of Na cannot form a stable interaction with either the HOMO or the LUMO of the graphite. On the other hand, the SOMO of Li and K can form a relatively stable interaction with either the HOMO or the LUMO of graphite. Why Li has a relatively stronger adsorption than K on graphite has also been interpreted on the basis of their molecular-orbital energy levels.     

Interaction of metal porphyrins with fullerene C60: a new insight

The electronic structure and bonding in the noncovalent, supramolecular complexes of fullerene C60 with a series of first-row transition metal porphines MP (M=Fe, Co, Ni, Cu, Zn) have been re-examined with DFT methods. A dispersion correction was made for the C60-MP binding energy through an empirical method (J. Comput. Chem. 2004, 25, 1463). Several density functionals and two types of basis sets were employed in the calculations. Our calculated results are rather different from those obtained in a recent paper (J. Phys. Chem. A 2005, 109, 3704). The ground state of C60.FeP is predicted to be high spin (S=2); the low-spin (S=0), closed-shell state is even higher in energy than the intermediate-spin (S=1) state. With only one electron in the Co-dz2 orbital, the calculated Co-C60 distance is in fact rather short, about 0.1 A longer than the Fe-C60 distance in high-spin C60.FeP. Double occupation of an M-dz2 orbital in MP prevents close association of any axial ligand, and so the Ni-C60, Cu-C60, and Zn-C60 distances are much longer than the Co-C60 one. The evaluated MP-C60 binding energies (Ebind) are 0.8 eV (18.5 kcal/mol) for M=Fe/Co and 0.5 eV (11.5 kcal/mol) for M=Ni/Cu/Zn (Ebind is about 0.2 eV larger in the case of C60-MTPP). They are believed to be reliable and accurate based on our dispersion-corrected DFT calculations that included the counterpoise (CP) correction. The effects of the C60 contact on the redox properties of MP were also examined.     

Nitrogen atom transfer from iron(IV) nitrido complexes: a dual-nature transition state for atom transfer

The mechanism of nitrogen atom transfer from four-coordinate tris(carbene)borate iron(IV) nitrido complexes to phosphines and phosphites has been investigated. In the absence of limiting steric effects, the rate of nitrogen atom transfer to phosphines increases with decreasing phosphine σ-basicity. This trend has been quantified by a Hammett study with para-substituted triarylphosphines, and is contrary to the expectations of an electrophilic nitrido ligand. On the basis of electronic structure calculations, a dual-nature transition state for nitrogen atom transfer is proposed, in which a key interaction involves the transfer of electron density from the nitrido highest occupied molecular orbital (HOMO) to the phosphine lowest unoccupied molecular orbital (LUMO). Compared to analogous atom transfer reactions from a 5d metal, these results show how the electronic plasticity of a 3d metal results in rapid atom transfer from pseudotetrahedral late metal complexes.     

Quantum chemical simulation of 1,3-dipolar cycloaddition of nitrones to alkenes and their (η6-arene)(tricarbonyl)-chromium complexes

Charge distribution and frontier orbital energies of styrene, C,N-diphenylnitrone, and their (arene)-(tricarbonyl)chromium complexes were calculated by quantum chemical methods. The difference in the HOMO and LUMO energies of the chromium complexes was found to be smaller than in the free ligands, and the reactions with (arene)(tricarbonyl)chromium complexes were characterized by higher rate and selectivity.     

Structural and spectroscopic characterization of an electrophilic iron nitrido complex

We have isolated and structurally characterized a terminal iron nitrido complex supported by a bulky tris(carbene)borate ligand. The electronic structure of this complex reveals that the a1 LUMO (formerly Fe(dz2)) is strongly stabilized by reduced antibonding interactions with the carbene sigma-donor ligands and configurational mixing (hybridization) with higher lying Fe 4s and 4p atomic orbitals. This unusual bonding interaction results in a low-lying Fe nitrido acceptor orbital (LUMO) that possesses electrophilic character. Reaction with PPh3 results in nitrogen atom transfer to the phosphine, supporting a reaction mechanism involving nucleophilic attack of the triphenylphosphine HOMO at the electrophilic LUMO of the iron nitrido complex.     

Electron smearing in DFT calculations: A case study of doxorubicin interaction with single‐walled carbon nanotubes

To address the choice of an appropriate value of electron smearing to facilitate self‐consistent field (SCF) convergence, we studied the interaction of doxorubicin with short armchair and zigzag single‐walled carbon nanotube models with closed caps, at the PWC/DNP level of density functional theory. By gradually reducing the electron smearing value from a large and most commonly used one of 0.005 Ha to zero (Fermi occupation), we monitored the changes in close contacts between the interacting species, total energy of the molecular system, highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) energy and isosurfaces, HOMO‐LUMO gap energy, and plots of electrostatic potential. It became evident that the commonly used smearing values of ≥0.001 Ha can alter the results significantly (for example, by one order of magnitude for HOMO–LUMO gap energy). We suggest the setting of electron smearing value at 0.0001 Ha, which does not imply too high computation cost and can guarantee the results close to the ones obtained with Fermi occupation. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Trifluorosulfane ligand as an analogue of the nitrosyl ligand: highly exothermic fluorine transfer reactions from sulfur to metal in the chemistry of SF3 metal carbonyls of the first row transition metals

The variety of known very stable PF(3) metal derivatives analogous to metal carbonyls suggests the synthesis of SF(3) metal derivatives analogous to metal nitrosyls. However, the only known SF(3) metal complex is the structurally uncharacterized (Et(3)P)(2)Ir(CO)(Cl)(F)(SF(3)) synthesized by Cockman, Ebsworth, and Holloway in 1987 and suggested by electron counting to have a one-electron donor SF(3) group rather than a three-electron donor SF(3) group. In this connection, the possibility of synthesizing SF(3) metal derivatives analogous to metal nitrosyls has been investigated using density functional theory. The [M]SF(3) derivatives with [M] = V(CO)(5), Mn(CO)(4), Co(CO)(3), Ir(CO)(3), (C(5)H(5))Cr(CO)(2), (C(5)H(5))Fe(CO), and (C(5)H(5))Ni analogous to known metal nitrosyl derivatives are all predicted to be thermodynamically disfavored with respect to the corresponding [M](SF(2))(F) derivatives by energies ranging from 19.5 kcal/mol for Mn(SF(3))(CO)(4) to 5.4 kcal/mol for Co(SF(3))(CO)(3). By contrast, the isoelectronic [M]PF(3) derivatives with [M] = Cr(CO)(5), Fe(CO)(4), Ni(CO)(3), (C(5)H(5))Mn(CO)(2), (C(5)H(5))Co(CO), and (C(5)H(5))Cu are all very strongly thermodynamically favored with respect to the corresponding [M](PF(2))(F) derivatives by energies ranging from 64.3 kcal/mol for Cr(PF(3))(CO)(5) to 31.6 kcal/mol for (C(5)H(5))Co(PF(3))(CO). The known six-coordinate (Et(3)P)(2)Ir(CO)(Cl)(F)(SF(3)) is also predicted to be stable relative to the seven-coordinate (Et(3)P)(2)Ir(CO)(Cl)(F)(2)(SF(2)). Most of the metal SF(3) complexes found in this work are singlet structures containing three-electron donor SF(3) ligands with tetrahedral sulfur coordination. However, two examples of triplet spin state metal SF(3) complexes, namely, the lowest energy (C(5)H(5))Fe(SF(3))(CO) structure and a higher energy Co(SF(3))(CO)(3) structure, are found containing one-electron donor SF(3) ligands with pseudo square pyramidal sulfur coordination with a stereochemically active lone electron pair.     

The effect of conjugated spacer on novel carbazole derivatives for dye‐sensitized solar cells: Density functional theory/time‐dependent density functional theory study

The ground‐state structure and frontier molecular orbital of D‐π‐A organic dyes, CFT1A, CFT2A, and CFT1PA were theoretically investigated using density functional theory (DFT) on B3LYP functional with 6‐31G(d,p) basis set. The vertical excitation energies and absorption spectra were obtained using time‐dependent DFT (TD‐DFT). The adsorptions of these dyes on TiO₂ anatase (101) were carried out by using a 38[TiO₂] cluster model using Perdew–Burke–Ernzerhof functional with the double numerical basis set with polarization (DNP). The results showed that the introduction of thiophene–thiophene unit (T–T) as conjugated spacer in CFT2A could affect the performance of intramolecular charge transfer significantly due to the inter‐ring torsion of T–T being decreased compared with phenylene–phenylene (P–P) spacer of CFP2A in the researhcers' previous report. It was also found that increasing the number of π‐conjugated unit gradually enhanced charge separation between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) of these dyes, leading to a high‐efficiency photocurrent generation. The HOMO–LUMO energy gaps were calculated to be 2.51, 2.37, and 2.50 eV for CFT1A, CFT2A, and CFT1PA respectively. Moreover, the calculated adsorption energies of these dyes on TiO₂ cluster were ～14 kcal/mol, implying that these dyes strongly bind to TiO₂ surface. Furthermore, the electronic HOMO and LUMO shapes of all dye–TiO₂ complexes exhibited injection mechanism of electron via intermolecular charge‐transfer transition. © 2012 Wiley Periodicals, Inc.     

Structural, electronic, and optical properties of 9-heterofluorenes: a quantum chemical study

Density-functional theory studies were applied to investigate the structural, electronic, and optical properties of 9-heterofluorenes achieved by substituting the carbon at 9 position of fluorene with silicon, germanium, nitrogen, phosphor, oxygen, sulfur, selenium, or boron. These heterofluorenes and their oligomers up to pentamers are highly aromatic and electrooptically active. The alkyl and aryl substituents of the heteroatom have limited influence, but the oxidation of the atom has significant influence on their molecular structures and properties. The highest occupied molecular orbital (HOMO)-lowest occupied molecular orbital (LUMO) interaction theory was successfully applied to analyze the energy levels and the frontier wave functions of these heterofluorenes. Most heterofluorenes belong to type B of interaction with low-lying LUMO and have the second kind of wave function. Carbazole and selenafluorene have type C of interaction with high-lying HOMO and the third kind of wave function. Types C and D of heterofluorenes, such as carbazole, oxygafluorene, sulfurafluorene, and selenafluorene also have high triplet state energies. The extrapolated HOMO and LUMO for polyheterofluorenes indicate that polyselenonafluorene has the lowest LUMO; polycarbazole has the highest HOMO; polyselenafluorene has the highest bandgap (E(g)); and polyborafluorene has the lowest E(g). Heterofluorenes and their oligomers and polymers are of great experimental interests, especially those having extraordinary properties revealed in this study.