Theoretical investigation on geometries and aromaticity of heterocyclic platinabenzenes

The electronic structure and properties of the heterocyclic platinabenzenes isomers have been investigated using the hybrid density functional B3LYP theory. Basic measures of aromatic character were derived from the structure and nucleus-independent chemical shift (NICS). The energetic criterion suggests that the ortho-isomer enjoy conspicuous stabilization where heteroatom is P or As. But the meta-isomer is most stable isomer, where heteroatom is N. The NICS values calculated at several points above the ring center to gave the result consistent with that cased on the relative energy. The atoms in molecules analysis indicates a correlation between NICS(1.0) and the electron density of ring critical point in all species. The natural bond orbital analysis has been used to study the bond characterizations in all species.     

Theoretical investigations on electronics structure and chemical bonding on iridathiabenzene and iridaoxabenzene

The electronic structure and properties of the iridathiabenzene and iridaoxabenzene isomers have been investigated using the hybrid density functional mpw1pw91 theory. The energetic aspect shows that trans-ortho-isomer is the most stable isomer. This is compatible with principles of minimum energy and minimum polarizability. Molecular orbital analysis shows a linear correlation between hardness and anisotropic polarizability values for Iridathiabenzene and iridaoxabenzene isomers. The structural and natural bond analysis (NBO) results illustrate electronic delocalization in these rings. Also, the study of non linear optical properties of these molecules indicate a good correlation between β_tot and E(HOMO) for iridathiabenzene. The results from natural bond orbital (NBO) analysis have provided insights into Ir—ligand, P—H_apical and P—H_basal bonding.     

Ab initio calculation of carbon clusters. II. Relative stabilities of fullerene and nonfullerene C24

Chemical stabilities of six low-energy isomers of C24 derived from global-minimum search are investigated. The six isomers include one classical fullerene (isomer 1) whose cage is composed of only five- and six-membered rings (56-MRs), three nonclassical fullerene structures whose cages contain at least one four-membered ring (4-MR), one plate, and one monocyclic ring. Chemical and electronic properties of the six C24 isomers are calculated based on a density-functional theory method (hybrid PBE1PBE functional and cc-pVTZ basis set). The properties include the nucleus-independent chemical shifts (NICS), singlet-triplet splitting, electron affinity, ionization potential, and gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital (HOMO-LUMO) gap. The calculation suggests that the neutral isomer 2, a nonclassical fullerene with two 4-MRs, may be more chemically stable than the classical fullerene (isomer 1). Analyses of molecular orbital NICS show that the incorporations of 4-MRs into the cage considerably reduce paratropic contributions from HOMO, HOMO-1, and HOMO-2, which are mainly responsible for the sign change in NICS from positive for isomer 1 (42) to negative (-19) for isomer 2, although C24 clusters satisfy neither 4N+2 nor 2(N+1)2 aromaticity rule. Anion photoelectron spectra of four cage isomers, one plate, one monocyclic ring, and one tadpole isomer, as well as three bicyclic ring isomers are calculated. The simulated photoelectron spectra of mono- and bicyclic rings (with C1 symmetry) appear to match the measured HOMO-LUMO gap (between the first and second band in the experimental spectra) [S. Yang et al., Chem. Phys. Lett. 144, 431 (1988)]. Nevertheless, the nonclassical fullerene isomers 3 and 4 apparently also match the measured vertical detachment energy (2.90 eV) reasonably well. These results suggest possible coexistence of nonclassical fullerene isomers with the mono- and bicyclic ring isomers of C24(-) under the experimental conditions.     

Theoretical studies on the structures,properties, and aromaticities of fluorinated arsabenzenes

The electronic structures and properties of the fluorinated arsabenzenes series have been investigated using the basis set 6–311+G(d,p) and hybrid density functional theory. The basic measures of aromatic character derived from molecular orbitals and magnetic criteria (anisotropic susceptibilities and nucleus-independent chemical shifts) are considered. The energy criteria suggest that the F3, F36, H36, and H3 isomers are the most stable isomers in the mono-, di-, tri-, and tetrafluorinated species, respectively. Analysis of χ_aniso and the HOMO-LUMO gaps showed that these were not compatible with NICS data. The NICS values show that aromaticity is greater in the fluorinated derivatives.     

A quantum chemistry study of ruthenabenzene complexes

The electronic structure and properties of the ruthenabenzenes and substituted ruthenabenzenes have been explored using the hybrid density functional B3LYP theory. Systematic studies on the substituent effect in para-substituted ruthenabenzenes complexes have been studied. The following substituents were taken into consideration: H, NO₂, CN, CHO, COOH, F, CH₃, OH, and NH₂. Basic measures of aromatic character were derived from the structure and nucleus-independent chemical shift (NICS). The NICS calculations indicate a correlation between NICS(1.5) and the hardness in all species. The atoms in molecule analysis indicates a correlation between r(Ru-C) bonds and the electron density of bond critical point in all species.     

Theoretical study of structure, bonding, and aromaticity of borazyne and B-substituted borazynes

The density functional theory (DFT) is used to study the geometries, and electronic structures of triplet and singlet of borazyne and B-substituted of borazyne. The aromaticity of these systems is analyzed in the light of nucleus-independent chemical shift (NICS), average of two-center indices (ATI). These methods show increasing of aromaticity in deactivating groups. The relation between electron density in ring critical point (RCP) and NICS(1.0) is observed. The most important interaction in these molecules has been investigated by natural bonding orbital method (NBO).     

A quantum chemical study of three isomers of N20

Ab initio molecular orbital (MP2) and density functional theory (B3LYP) calculations using different basis sets have been employed to study the structures, energetics and vibrational frequencies of the large homonuclear polynitrogen compound, N₂₀. In the present study, three distinct forms were found to represent local minima on the potential energy surface. They are the fullerene-type cage form of I_h symmetry, a corannulene-like bowl form of C_5v symmetry, and a ring isomer with D₅ symmetry of which the cage form turns out to be the highest energy form. Both the bowl and ring forms are calculated to be more stable than the cage form by about 200 kcal/mol. The molecular properties calculated for these isomers may serve as valuable predictions for future experimental searches for new high energy density materials (HEDM).     

Theoretical study on the aromaticity of the bimetallic clusters X2M2 (X=Si, Ge, M=Al, Ga)

Four neutral bimetallic clusters X₂M₂ (X=Si, Ge, M=Al, Ga) are investigated using density functional theory (DFT) and post-HF methods. The calculated results show that each of four X₂M₂ species has two energetically close stable isomers with rhombic structure (D_2h symmetry) and trapezoidal structure (C_2v symmetry) respectively. For the Ge₂Al₂ species the rhombic (D_2h) isomer is the ground state, whereas for other three species Ge₂Ga₂, Si₂Al₂, and Si₂Ga₂, the trapezoidal (C_2v) isomers are the ground states. The calculated magnetic susceptibility anisotropy (χ_anis) and nucleus-independent chemical shift (NICS) indicate that a strong diatropic ring current exists in the two heterocyclic planar isomers, suggesting they are highly aromatic. A detailed molecular orbital analysis further reveals that both heterocyclic isomers possess multiple aromaticity derived from one delocalized π MOs and two delocalized σ MOs.     

The maximum hardness and minimum polarizability principles in accordance with the Bent rule

Energies, hardness values and polarizabilities for different isomers of SF₄, SF₄O molecules and a family of PCl_xF_5−x, (x=1,2,3,4) compounds are calculated at the density functional level of theory (B3LYP) using different basis sets (6-31G**, 6-311G**, 6-31++G and 6-311++G). For all molecules and in all cases the Bent rule is confirmed with the stability orders, which is obtained from the calculated energies. It is also found that for each molecule the isomer in which the more electronegative atom occupies the axial position has maximum hardness; and therefore according to the Maximum Hardness Principle (MHP), it must be the most stable isomer. This result is consistent with the Bent rule, but for some cases in which there are more than two isomers, the obtained trends for the stability are not in accordance with this rule. On the other hand, since the calculation of dipole polarizability is very sensitive to the quality of the selected basis set, only for those calculations in which diffuse basis sets are used the obtained trends of stability from the Minimum Polarizability Principle (MPP) are the same as those which are predicted from the Bent rule. It seems that when there are more than two isomers for a molecule both the MHP and MPP can only predict the most stable isomer and cannot predict the correct trend for the stability of the isomers.     

The substitution effect on heavy versions of cyclobutadiene

The structural and electronic characteristics of the M₄R₄ with M = C, Si, Ge and R = H, F, Cl, OH, MgH, BH₂, NH₂ Li are studied by ab initio calculations based on density functional theory using the hybrid B3LYP functional and 6‐311++G* basis set. The results of natural bond orbital analysis (NBO), the analysis of structural parameters (bond lengths and bond angle, dihedral angle), atoms in molecules based theory (AIM) topological parameters such as the characteristics of bond critical points (electron densities and their Laplacians) leads to the conclusion that interamolecular interaction due with OH, MgH₂, and NH₂ substitutions are different from other derivatives. The cyclobutadiene ring and its heavy versions display very different geometrical structures. The cations and anions of C₄R₄, Si₄R₄, and Ge₄R₄ were studied to determine the effect of substitution on structure and stability. It has been shown that substitutions with ability of forming hydrogen bonds, bound to the silicon and germanium atoms are able to change the properties of the four‐member rings more than the other substitutions. Comprehensive analyses of calculated NICS (nucleus‐independent chemical shifts), polarizability, electrophilicity, and electron density are done. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

A Theoretical Study of NBO,NICS, and 14N NQR Parameters of Adenine Tautomers in the Gas Phase via DFT

The natural bond orbital (NBO) analysis, nucleus independent chemical shift (NICS), and ¹⁴N NQR parameters of the most stable tautomers of adenine in the gas phase were predicted using density functional theory method. The NBO analysis revealed that the resonance interaction between lone pair of the nitrogen atom and empty non‐Lewis NBO increases with increasing the p character of the nitrogen lone pair. The present investigation indicated the π clouds in both the considered heterocyclic rings containing six electrons, and these tautomers has the aromatic character. The NICS study utilizing the gauge‐invariant atomic orbital method showed that there are diatropic currents in the heterocyclic rings of the tautomers, so we determined the order of overall aromaticity of these tautomers. The results of NQR parameter calculations showed three parameters are effective on nuclear quadrupole coupling constant; the p character value of lone pair electrons of nitrogens, and the related occupancies and whenever, the lone pair electrons of nitrogens participate in the formation of chemical bond and/or π system of the ring, the q_zz and consequently its χ decreases.     

A Quantum Chemical Study of N14 Cluster

Ab initio molecular orbital theory and density functional theory have been employed to study N₁₄ cluster with low spin at the HF/6-31G*, B3LYP/6-31G*, B3PW91/6-31G*, BP86/6-31G*, and BHLYP/6-31G* levels of theory. Twelve isomers were studied, including one previously investigated cage molecule. The most stable isomer of N₁₄ is a C _2h -symmetric molecule that contains two separated five-membered nitrogen rings connected by a —N=N—N=N— bridge. The second, third, and fifth most stable isomers each have one five-membered nitrogen ring. The theoretical results suggest that the five-membered nitrogen ring gives rise to a particularly stable structural unit, and the more side chains that the five-membered nitrogen ring links with, the less stable the structure will become.     

Spectroscopic studies and molecular orbital analysis on platinanaphthalenes and ring-fused B-N platinanaphthalenes

The electronic structure and properties of the platinanaphthalenes, and ring-fused B-N platinanaphthalenes isomers have been explored using the hybrid density functional mpw1pw91 theory. The energetic aspect shows that stability of Ptb and Pta isomers are isoenergetic in platinanaphthalenes. On the other hand, BNa1 isomer is the most stable isomer of ring-fused B-N platinanaphthalenes. This is compatible with principles of minimum energy and minimum polarizability in ring-fused B-N platinanaphthalenes. Molecular orbital analysis shows increasing of hardness in ring-fused B-N platinanaphthalenes isomers. Also, electronic spectra analysis indicates that, in all the molecules HOMO-1 → LUMO transition makes the major contribution in most intense electronic transition.     

Solvent impact on the planarity and aromaticity of free and monohydrated zinc phthalocyanine: a theoretical study

A theoretical investigation on the planarity of molecular structure of zinc phthalocyanine (ZnPc) and its aromaticity has been performed using B3LYP and M06-2X density functionals combined with selected Pople-type basis sets. The effect of the applied calculation method on the optimized structure of ZnPc and ZnPc???H₂O, both in the gas phase and in the polar solvent, was analyzed. To quantify the aromaticity of the ZnPc and ZnPc???H₂O complexes, both the geometric and magnetic criteria, i.e., Harmonic Oscillator Model of Aromaticity (HOMA) index and the nucleus-independent chemical shift (NICS) values at the centers or 1 Å above the centers of structural subunits, were calculated. The energies of highest energy occupied molecular orbital (HOMO) and lowest energy unoccupied molecular orbital (LUMO) and energy gaps were also estimated. The results show that the free ZnPC molecule is flat in the gas phase and nonplanar in the polar environments (DMSO and water). ZnPC???H₂O is nonpolar in the gas phase and polar solvent which is in agreement with recent X-ray reports. Both HOMA and NICS indexes indicate the presence of highly aromatic macrocycle and benzene rings while these parameters for pyrrolic ring are significantly smaller than in free pyrrole. The presence of polar solvents practically does not change aromaticity of the ring subunits of the studied compounds.     

Theoretical Study of the Solvent Effect on the Electronic and Vibrational Properties of [CpFe(CO)<Subscript>2</Subscript>(NCS)] and [CpFe(CO)<Subscript>2</Subscript>(SCN)] Linkage Isomers

The present study illustrates the stability of [CpFe(CO)₂(NCS)] and [CpFe(CO)₂(SCN)] linkage isomers by the use of MPW1PW91 quantum method in the gas and solution phases. Our results reveal that the [CpFe(CO)₂(NCS)] isomer is more stable than the [CpFe(CO)₂(SCN)] isomer. Based on the polarizable continuum model, the effect of the solvent polarity on the stability, structural parameters, frontier orbital energies, and vibrational modes of carbonyl ligands (ν_CO) of these linkage complexes is explored. The molecular orbital analysis suggests that the major contributions to HOMO and LUMO arise from the ambidentate ligand and Fe in two isomers, respectively. In addition, the bonding interaction between the CpFe(CO)₂ fragment and the ambidentate ligand is studied by means of the energy decomposition analysis. The back-bonding effect in Fe–CO bonds is revealed in the calculation of the quadrupole polarization of the carbon atom by the QTAIM analysis. The character of Fe–N and Fe–S bonds in these complexes is analyzed by the natural bond orbital analysis.     

A quantum chemical study of Cr(CO)3(B3N3H6 ? n F n ) (n = 1–3) complexes

The electronic structure and properties of Cr(CO)₃(B₃N₃H_{6 ? n} F _n ) (n = 1?C3) complexes have been explored using hybrid density functional B3LYP theory. Calculations indicate B-fluorinated isomers are more stable, and less polarizable, than N-fluorinated isomers. The aromatic natures of the borazine rings have been analyzed by nucleus independent chemical shift (NICS). The atoms in molecules (AIM) analysis indicates that Cr-C and Cr-N bonds distance is well correlated with the electron density of critical point (??_cp) in all species.     

Comparative semiempirical and DFT study of styrylnaphthalenes and styrylquinolines and their photocyclization products

Semiempirical molecular orbital (PM3, PM6, and RM1) and density functional theory (DFT) (B3LYP/6‐31G*) studies are carried out for 1‐ and 2‐styrylnaphthalenes and their aza‐derivatives—2‐ and 4‐styrylquinolines. Relative stabilities of three isomeric forms: E‐ and Z‐isomers and the closed‐ring dihydrocyclophotoproduct (derivative of dihydrophenanthrene) are calculated. Compared to PM3, PM6 and especially RM1 understate heats of formation; in some cases, PM6 and RM1 even place Z‐isomer in energy below E‐isomer. PM3 rather close to DFT predicts heats of isomerization reaction, whereas PM6 and especially RM1 underestimate these values. Semiempirical methods in comparison with DFT markedly underestimate heats of cyclization reaction; however, reproduce trends in relative stabilities of different isomers in dependence on the structure of styrylnaphthalenes and styrylquinolines. Qualitative correlation is found between calculated relative stabilities of the closed‐ring forms (heats of cyclization reaction) and experimental data: cyclized products with low heats of cyclization are observed in steady‐state photolysis and those with high heats of cyclization are not. In the latter case, the closed‐ring compounds, if formed in the excited state, due to thermal instability decompose rapidly with ring opening in the ground state that prevents their observation. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Diarylethene Self‐Assembled Monolayers: Cocrystallization and Mixing‐Induced Cooperativity Highlighted by Scanning Tunneling Microscopy at the Liquid/Solid Interface

Stimulus control over 2D multicomponent molecular ordering on surfaces is a key technique for realizing advanced materials with stimuli‐responsive surface properties. The formation of 2D molecular ordering along with photoisomerization was monitored by scanning tunneling microscopy at the octanoic acid/highly oriented pyrolytic graphite interface for a synthesized amide‐containing diarylethene, which underwent photoisomerization between the open‐ and closed‐ring isomers and also a side‐reaction to give the annulated isomer. The nucleation (K_n) and elongation (K_e) equilibrium constants were determined by analysis of the concentration dependence of the surface coverage by using a cooperative model at the liquid/solid interface. It was found that the annulated isomer has a very large equilibrium constant, which explains the predominantly observed ordering of the annulated isomer. It was also found that the presence of the closed‐ring isomer induces cooperativity into the formation of molecular ordering composed of the open‐ring isomer. A quantitative analysis of the formation of ordering by using the cooperative model has provided a new view of the formation of 2D multicomponent molecular ordering.     

DFT studies and AIM analysis of AlN-polycycles

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