Electronic structure,chemical bonding,and oxidation numbers of first‐row transition metals in [MePIm2] complexes and their cations

The qualitative structures of the upper one‐electron energy levels of imidazole‐coordinated first‐row transition metal porphyrin [MePIm₂] complexes established in the present study have shown that the second oxidation number of the first‐row transition metals in the neutral complexes do not change in their cations and double cations. It was found that occupied orbitals of the density functional theory method obtained with B3LYP functional are not correctly ordered. Therefore, they cannot be used in investigations of the orbital structure of the upper molecular orbitals. A qualitative analysis of density functional theory method wave functions in terms of Mulliken and natural charges of atoms, together with an analysis of electrostatic potentials of the neutral [MePIm₂] complex, its single and double cations, demonstrates that the highest occupied orbitals of these complexes are mainly formed by atomic orbitals of the porphyrin ring atoms. Therefore, transition metal atoms are not active in chemical reactions with these complexes unless the 3d electrons of transition metal atoms are excited, for example by light. A mechanism of an electron transfer reaction that occurs between a heme cytochrome and Fe‐oxide mineral surface is discussed in the light of the obtained results. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Why do the second row transition metal atoms prefer 5s14dm+1 to 5s24dm?

Numerical Hartree-Fock (NHF) calculations have been performed for 332 ground and low-lying excited states of the fifth period atoms Rb through Xe, with our special interest in the states arising from the 5s ²4d ^m , 5s ¹4d ^m +1, and 5s ⁰ 4d ^m +2 configurations of the second row transition metal atoms. Among various properties, orbital energies and mean values ofr of the outermost orbitals of each symmetry are presented as well as total energies. It is discussed in some detail why the second row transition metal atoms have a tendency to prefers ¹ d ^m +1 as the ground configuration in contrast to the preferreds ² d ^m configuration in the first row transition metal atoms. Our systematic NHF computations reported in this and the previous papers conclude that the Hartree-Fock method correctly predicts the experimental ground state of the atoms He through Xe with the sole exception for Zr.     

Excited state properties,fluorescence energies,and lifetimes of a poly(fluorene‐phenylene), based on TD‐DFT investigation

The structural and electronic properties of fluorene‐phenylene copolymer (FP)_n, n = 1–4 were studied by means of quantum chemical calculations based on density functional theory (DFT) and time dependent density functional theory (TD‐DFT) using B3LYP functional. Geometry optimizations of these oligomers were performed for the ground state and the lowest singlet excited state. It was found that (FP)_n is nonplanar in its ground state while the electronic excitations lead to planarity in its S₁ state. Absorption and fluorescence energies were calculated using TD‐B3LYP/SVP and TD‐B3LYP/SVP+ methods. Vertical excitation energies and fluorescence energies were obtained by extrapolating these values to infinite chain length, resulting in extrapolated values for vertical excitation energy of 2.89 and 2.87 eV, respectively. The S₁ ← S₀ electronic excitation is characterized as a highest occupied molecular orbital to lowest unoccupied molecular orbital transition and is distinguishing in terms of oscillator strength. Fluorescence energies of (FP)_n calculated from TD‐B3LYP/SVP and TD‐B3LYP/SVP+ methods are 2.27 and 2.26 eV, respectively. Radiative lifetimes are predicted to be 0.55 and 0.51 ns for TD‐B3LYP/SVP and TD‐B3LYP/SVP+ calculations, respectively. These fundamental information are valuable data in designing and making of promising materials for LED materials. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Combined experimental and computational modeling studies on 3,5‐dimethyl‐pyrazole‐1‐carbodithioic acid benzyl ester

The title compound, 3,5‐Dimethyl‐pyrazole‐1‐carbodithioic acid benzyl ester, has been synthesized and structurally characterized by X‐ray single crystal diffraction, elemental analysis, IR spectra, and UV‐Vis spectrum. The crystal belongs to orthorhombic, space group P212121, with a = 5.3829(15), b = 11.193(3), c = 21.824(6) Å, V = 1315.0(6) ų, and Z = 4. The molecules are connected via intermolecular C–H···N hydrogen bonds into 1D infinite chains. The crystal structure is consolidated by the intramolecular C–H···S hydrogen bonds. Furthermore, Density functional theory (DFT) calculations of the structure, stabilities, orbital energies, composition characteristics of some frontier molecular orbitals and Mulliken charge distributions of the title compound were performed by means of Gaussian 03W package and taking B3LYP/6‐31G(d) basis set. The time‐dependent DFT (TD‐DFT) calculations have been employed to calculate the electronic spectrum of the title compound, and the UV‐Vis spectra has been discussed on this basis. The results show that DFT method at B3LYP/6‐31G(d) level can well reproduce the structure of the title compound. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Quantum study on the electrocyclic reactions of CH2CHCHCHX (XH,F, Cl)

The microscopic mechanisms of the electrocyclic reactions for cis‐1,3‐butadiene and its monofluoro‐, monochloroderivatives have been studied by density functional theory (DFT), using the B3LYP method and 6‐311++G** basis sets. We optimized the geometric configurations of reactants, transition states, and products; verified all the probable transition states through vibrational analysis; and calculated the relative single‐point energies at the QCISD(T)/6‐311++G**//B3LYP/6‐311++G**. The results show that the monofluoro‐, monochloroderivatives of cis‐1,3‐butadiene both have two conformers; the reactant favors the electrocyclic reaction when one outboard hydrogen atom of the CH₂ groups is substituted by the fluorine or chlorine atom. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

NO migration from N‐methyl‐N‐nitrosobenzene‐sulfonamide to 3,6‐dibromocarbazole: Concerted or stepwise reaction path?

The NO migration from N‐methyl‐N‐nitrosobenzene‐sulfonamide to 3,6‐dibromocarbazole was proposed in a recent literature to follow a stepwise reaction path. However, the present density functional theory calculations at the MP2/6–31G(d,p)//B3LYP/6–31G(d,p) level show that this reaction exclusively proceeds via a concerted mechanism involving a four‐membered ring transition state. The calculated barrier is in good agreement with the experimental finding. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Theoretical Study on Gas—phase Pyrolytic Reactions of N—Ethyl,Isopropyl and N—t—Butyl Substituted 2—Aminopyrazine

Density functional theory (DFT) and ab initio methods were used to study gas‐phase pyrolytic reaction mechanisms of iV‐ethyl, N‐isopropyl and N‐t‐butyl substituted 2‐aminopyrazine at B3LYP/6–31G* and MP2/6–31G*, respectively. Single‐point energies of all optimized molecular geometries were calculated at B3LYP/6–311 + G(2d,p) level. Results show that the pyrolytic reactions were carried out through a unimolecular first‐order mechanism which were caused by the migration of atom H(17) via a six‐member ring transition state. The activation energies which were verified by vibrational analysis and correlated with zero‐point energies along the reaction channel at B3LYP/6–311 + G(2d,p) level were 252.02 kJ. mo^?1 (N‐ethyl substituted), 235.92 kJ‐mol^?1 (N‐t‐isopropyl substituted) and 234.27 kJ‐mol^?1 (N‐t‐butyl substituted), respectively. The results were in good agreement with available experimental data.     

Origin of the chemical shift in X‐ray absorption near‐edge spectroscopy at the Mn K‐edge in manganese oxide compounds

The absorption edge in Mn K‐edge X‐ray absorption spectra of manganese oxide compounds shows a shift of several electronvolts in going from MnO through LaMnO₃ to CaMnO₃. On the other hand, in X‐ray photoelectron spectra much smaller shifts are observed. To identify the mechanisms that cause the observed chemical shifts, 1s ionization as well as 1s → “4p” transition energies have been determined by electronic structure calculations on embedded Mn ions and embedded MnO₆ clusters. Systematic variation of the cluster geometry and the cluster embedding showed that the chemical shifts are predominantly determined by two effects: the changes in the Mn 3d occupation and the changes in the Madelung potential. The large chemical shift in the 1s → 4p transition energies between different materials occurs because the two effects do not compensate each other. The chemical shifts obtained for the embedded MnO₆ clusters agree reasonably with the experimental shifts. The small sensitivity to the material observed for the Mn 1s ionization energies is explained by the near cancellation of the effects of the Madelung potential and the 3d occupation of the Mn ion. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

Long‐range corrected functionals satisfy Koopmans' theorem: Calculation of correlation and relaxation energies

In this article, we show that the long‐range‐corrected (LC) density functionals LC‐BOP and LCgau‐BOP reproduce frontier orbital energies and highest‐occupied molecular orbital (HOMO)—lowest‐unoccupied molecular orbital (LUMO) gaps better than other density functionals. The negative of HOMO and LUMO energies are compared with the vertical ionization potentials (IPs) and electron affinities, respectively, using CCSD(T)/6‐311++G(3df,3pd) for 113 molecules, and we found LC functionals to satisfy Koopmans' theorem. We also report that the frontier orbital energies and the HOMO‐LUMO gaps of LC‐BOP and LCgau‐BOP are better than those of recently proposed ωM05‐D (Lin et al., J. Chem. Phys. 2012, 136 , 154109). We express the exact IP in terms of orbital relaxation, and correlation energies and hence calculate the relaxation and correlation energies for the same set of molecules. It is found that the LC functionals, in general, includes more relaxation effect than Hartree–Fock and more correlation effect than the other density functionals without LC scheme. Finally, we scan μ parameter in LC scheme from 0.1 to 0.6 bohr^?1 for the above test set molecules with LC‐BOP functional and found our parameter value, 0.47 bohr^?1, is usefully applicable to our tested systems. © 2013 Wiley Periodicals, Inc.     

Quantum‐chemical,IR, NMR,and X‐ray diffraction studies on 2‐(4‐chlorophenyl)‐1‐methyl‐ 1H‐benzo[d]imidazole

The title molecule, 2‐(4‐chlorophenyl)‐1‐methyl‐1H‐benzo[d]imidazole (C₁₄H₁₁ClN₂), was prepared and characterized by ¹H NMR, ¹³C NMR, IR, and single‐crystal X‐ray diffraction. The molecular geometry, vibrational frequencies, and gauge including atomic orbital (GIAO) ¹H and ¹³C NMR chemical shift values of the title compound in the ground state have been calculated by using the Hartree‐Fock (HF) and density functional theory (DFT/B3LYP) method with 6‐31G(d) basis sets, and compared with the experimental data. The calculated results show that the optimized geometries can well reproduce the crystal structural parameters, and the theoretical vibrational frequencies and GIAO ¹H and ¹³C NMR chemical shifts show good agreement with experimental values. The energetic behavior of the title compound in solvent media has been examined using B3LYP method with the 6‐31G(d) basis set by applying the Onsager and the polarizable continuum model (PCM). Besides, molecular electrostatic potential (MEP), frontier molecular orbitals (FMO) analysis, and nonlinear optical (NLO) properties of the title compound were investigated by theoretical calculations. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Performance of DFT hybrid functionals in the theoretical treatment of H‐bonds: Analysis term‐by‐term

The performance of DFT to reproduce noncovalent interactions like H‐bonds was assessed. Three DFT hybrid functionals (B3LYP, BHandHLYP from Jaguar5.5, and BHandHLYP from G03) were used to calculate the interaction energies and H‐bond distances of several host–guest complexes theoretically designed. Two reference systems (whose experimental data were available) were also calculated at the same levels of theory. In all the cases B3LYP and BHandHLYP from G03 gave rather the same results of interaction energy and distances, whereas the functional BHandHLYP from Jaguar5.5 overestimated the interaction energies and underestimated the H‐bond distances. Considering the construction of each functional, the terms responsible for such differences are the gradient correction to the exchange functional (Becke88) and the correlation functional (VWN, LYP) and not the Hartree–Fock contribution. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Theoretical study on the molecular structures of toluene,para‐fluorotoluene,para‐chlorotoluene,and 4‐methylpyridine and their sixfold internal rotational barriers

Two kinds of sixfold internal rotational configurations of toluene, para‐fluorotoluene, para‐chlorotoluene, and 4‐methylpyridine were calculated using Hartree–Fock (HF), second‐order Møller–Plesset (MP2), and Beck's three parameter hybrid functional using the LYP correlation functional (B3LYP) theory methods with various high‐level basis sets. Structures and energies were compared for different configurations. Calculations indicate that the orthogonal configuration has a local minimum while the planar configuration is a transition structure. Furthermore, geometries of the orthogonal and the planar configurations are quite similar, except for a methyl CH bond. Sixfold internal rotational barriers were calculated from the energy difference of two different configurations. For the calculated results, HF methods underestimated the rotational barriers, but MP2 calculations overestimated them. However, the density functional theory (DFT) method is a reliable method since the calculated internal rotational barriers are similar to the experimental ones. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 77: 772–778, 2000     

反式-1,4,5,8-四硝基-1,4,5,8-四氮杂萘烷异构体热解机理与稳定性的理论研究

运用密度泛函理论和半经验分子轨道方法,对一系列高能杂环硝胺—反式-1,4,5,8-四硝基-1,4,5,8-四氮杂萘烷异构体的热解机理和稳定性进行了系统地计算研究。在B3LYP/6-31G**和PM3水平上,分别计算了标题物的化学键离解能(BDE)和热解反应活化能(Ea),并根据BDE和Ea数值考察了硝胺取代基对化合物稳定性和热解机理的影响;同时,还详细考察了BDE与Ea、化学键重叠布居数、前线轨道能级以及能隙之间的相关性。结果表明,由BDE、Ea和静态电子结构参数推断的标题物热稳定性和热解机理的结论基本是一致的,N-NO2键均裂是标题物的热解引发步骤,间位取代异构体较对位取代异构体稳定,而邻位取代的异构体稳定性最差。     

Theoretical investigation of solvent effects on tautomeric equilibrium of 2‐diazo‐4,6‐dinitrophenol

The density functional theory has been used to study the tautomeric equilibrium of 2‐diazo‐4,6‐dinitrophenol(DDNP) in the gas phase and in 14 solvents at the B3LYP/6‐31G* level. The solvent effects on the tautomeric equilibria were investigated by the self‐consistent reaction field theory (SCRF) based on conductor polarized continuum model (CPCM) in apolar and polar solvents and by the hybrid continuum‐discrete model in protic solvent, respectively. Solvent effects on the computed molecular properties, such as molecular geometries, dipole moments, E_LUMO, E_HOMO, total energies for DDNP tautomers and transition state, tautomerization energies and solvation energies have been found to be evident. The tautomeric equilibrium of DDNP is solvent‐dependent to a certain extent. The tautomer I (cyclic azoxy form) is preferred in the gas phase, while in nonpolar solvents tautomer I and II (quinold form) exist in comparable amounts, and in highly polar solvents, the tautomeric equilibrium is shifted in favor of the more polar tautomer II . © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Experimental and theoretical investigation of the molecular and electronic structure of N′‐benzylidene‐N‐[4‐(3‐methyl‐3‐phenyl‐cyclobutyl)‐thiazol‐2‐yl]‐chloro‐acetic acid hydrazide

The title compound, N′‐benzylidene‐N‐[4‐(3‐methyl‐3‐phenyl‐cyclobutyl)‐thiazol‐2‐yl]‐chloro‐acetic acid hydrazide, has been synthesized and characterized by elemental analysis, IR, ¹H and ¹³C NMR, and X‐ray single crystal diffraction. The compound crystallizes in the orthorhombic space group P 21 21 21 with a = 5.8671 (3) Å, b = 17.7182 (9) Å, and c = 20.6373 (8) Å. Moreover, the molecular geometry from X‐ray experiment, the molecular geometry, vibrational frequencies, and gauge‐including atomic orbital ¹H and ¹³C chemical shift values of the title compound in the ground state have been calculated by using the Hartree–Fock and density functional methods (B3LYP) with 6‐31G(d) and 6‐31G(d,p) basis sets. The results of the optimized molecular structure are exhibited and compared with the experimental X‐ray diffraction. Besides, molecular electrostatic potential, Frontier molecular orbitals, and thermodynamic properties of the title compound were determined at B3LYP/6‐31G(d) levels of theory. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Nature of closed‐ and open‐shell interactions between noble metals and rare gas atoms

Interactions between noble metals and rare gases have become an interesting topic over the last few years. In this work, a computational study of the open‐shell (d¹⁰s¹) and closed‐shell (d¹⁰s and d¹⁰s²) noble metals (M = Cu, Ag, and Au) with three heaviest rare gas atoms (Rg = Kr, Xe, and Rn) has been performed. Potential energy curves based on ab initio [MP2, MP4, QCISD, and CCSD(T)] and DFT functionals (M06‐2X and CAM‐B3LYP) were obtained for ionic and neutral AuXe complexes. Dissociation energies indicate that neutral metals have the lowest and cationic metals have the highest affinities for interaction with rare gas atoms. For the same metals, there is a continuous increase in dissociation energies (D_e) from Kr to Rn. The nature of bonding and the trend of D_e and equilibrium bond lengths (R_e) have been interpreted by means of quantum theory of atoms in molecules, natural bond orbital, and energy decomposition analysis. © 2013 Wiley Periodicals, Inc.     

A density functional theory study on diazotization and nitration of 3,5‐diamino‐1,2,4‐triazole

The reaction mechanism, thermodynamic and kinetic properties for diazotization and nitration of 3,5‐diamino‐1,2,4‐triazole were studied by a density functional theory. The geometries of the reactants, transition states, and intermediates were optimized at the B3LYP/6‐31G (d, p) level. Vibrational analysis was carried out to confirm the transition state structures, and the intrinsic reaction coordinate (IRC) method was used to explore the minimum energy path. The single‐point energies of all stagnation points were further calculated at the B3LYP (MP2)/6‐311+G (2d, p) level. The statistical thermodynamic method and Eyring transition state theory with Wigner correction were used to study the thermodynamic and kinetic characters of all reactions within 0–25°C. Two reaction channels are computed, including the diazotization and nitration of 3‐NH₂ or 5‐NH₂, and there are six steps in each channel. The reaction rate in each step is increased with temperature. The last step in each channel is the slowest step. The first, second, and fifth steps are exothermic reactions, and are favored at lower temperature in the thermodynamics. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011