Calculation of the proton affinities of polyfluorobenzenes and the activation energies for 1,2-hydrogen shifts in their carbocations

Isolated polyfluorobenzene (PFB) molecules and their protonated forms are investigated by the AM1 method with full geometry optimization. The proton affinities of PFB are estimated for different protonated positions. The proton affinity of PFB averaged over all isomers is shown to decrease monotonically as the number of fluorine atoms in the molecule increases. The relative populations of different isomers of arenonium ions (AI) formed by PFB protonation are determined. From the calculated data, the value of ⁺ for the F atom in theipso-position is estimated as 1.00. The activation energies of the 1,2-hydrogen shifts in AI are calculated. The dependences of the proton affinity and the activation energies of 1,2-hydrogen shifts on the number of halogen atoms are found to have distinct characters for PFB and polychlorobenzenes. The physical reasons for these difference are discussed.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 11, pp. 1878–1882, November, 1993.     

Electronic structure of polychlorophenols and their carbocations

Polychlorophenols (PCP) with different numbers of chlorine atoms and their protonated forms have been calculated by the AM1 method with full optimization of geometry. The proton affinity of PCP with various types of coordination of the proton to the molecule has been estimated. The calculations show that the proton affinity averaged over isomers decreases monotonically as the number of chlorine atoms increases. Based on the calculations of the carbocations, the ⁺ constant of the OH-substituent at theipso-position is equal to 0.69. There is fair agreement between the relative energies of the isomers of the carbocations determined from the constants of the substituents and calculated by the AM1 method. These energies can be recommended for qualitative estimation of the proton affinity in definite positions of aromatic molecules with Cl- and OH-substituents.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 7, pp. 1198–1201, July, 1994.     

Theoretical study on the structures and properties of the isomers of Nn(CH)8−nH8 (n = 0–7)

With replacement of N atoms by CH groups in the most stable chain isomer of N8H8, 34 possible isomers of N_n(CH)_8−nH₈ (n = 0–7) have been designed and optimized at the B3LYP/6-311++G** level of theory. The natural bond orbital (NBO) and atoms in molecules (AIM) analysis are carried out to study the bonding nature and relative stabilities of these conformers. G3MP2 method is applied to calculate energies and heats of formation. The results indicate that the hyperconjugation effect from lone pairs of nitrogen atoms to germinal C–N bonds is the major factor which caused the change of the C–N bond length. With the more replacement of nitrogen atoms by CH groups, the heats of formation of the isomers of N_n(CH)_8−nH₈ (n = 0–7) decrease gradually, but the energies increase linearly.     

硼碳团簇BnC2 (n＝1～6)的理论研究

应用密度泛函理论在B3LYP/6-311＋G(d)水平上研究了硼碳团簇B_nC₂ (n＝1～6)的几何结构、生长机制和相对稳定性. 计算结果表明, 对于n＝2～6的簇, 平面多环状构型为最稳定的结构, 其中C原子分布于环的顶点、有尽可能多的三配位硼原子和尽可能多的B—C键. 碳原子作为杂原子倾向掺杂于团簇的顶点位置, 它的掺杂不改变硼团簇的主体结构. 与平面多环状结构相比, 随着簇尺寸的增大, 三维结构和线性链结构更不稳定. 在低能线性结构中, C原子位于链两侧的第二个位置. 计算的碎片分裂能、递增键能以及HOMO-LUMO能隙表明, B₄C₂为幻数簇.     

An ab initio study of the structures and enthalpies of the hydrogen-bonded complexes of the acids H2O,H2S,HCN, and HCl with the anions OH−, SH−, CN−, and Cl−

Ab initio calculations at second-order Møller-Plesset perturbation theory with the 6-31 + G(d,p) basis set have been performed to determine the equilibrium structures and energies of a series of negative-ion hydrogen-bonded complexes with H₂O, H₂S, HCN, and HCl as proton donors and OH^–, SH^–, CN^–, and Cl^– as proton acceptors. The computed stabilization enthalpies of these complexes are in agreement to within the experimental error of 1 kcal mol^–1 with the gas-phase hydrogen bond enthalpies, except for HOHOH^–, in which case the difference is 1.8 kcal mol^–1. The structures of these complexes exhibit linear hydrogen bonds and directed lone pairs of electrons except for complexes with H₂O as the proton donor, in which cases the hydrogen bonds deviate slightly from linearity. All of the complexes have equilibrium structures in which the hydrogen-bonded proton is nonsymmetrically bound, although the symmetric structures of HOHOH^– and ClHCl^– are only slightly less bound than the equilibrium structures. MP2/6-31 + G(d,p) hydrogen bond energies calculated at optimized MP2/B-31 + G(d,p) and at optimized HF/6-31G(d) geometries are similar. Using HF/6-31G(d) frequencies to evaluate zero-point and thermal vibrational energies does not introduce significant error into the computed hydrogen bond enthalpies of these complexes provided that the hydrogen-bonded proton is definitely nonsymmetrically bound at both Hartree-Fock and MP2.     

Electronic spectra of the linear magnesium-containing carbon chains MgC2nH (n = 1–5): A CASPT2 study

Complete active space self-consistent-field (CASSCF) approach has been used for the geometry optimization of the X²Σ⁺ and A²Π electronic states for the linear magnesium-containing carbon chains MgC_2nH (n = 1–5). Multireference second-order perturbation theory (CASPT2) has been used to calculate the vertical excitation energies from the ground to selected seven excited states, as well as the potential energy curves of two ²Σ⁺ and two ²Π electronic states. The studies indicate that the vertical excitation energies of the A²Π ← X²Σ⁺ transition for MgC_2nH (n = 1–5) are 2.837, 2.793, 2.767, 2.714, and 2.669 eV, respectively, showing remarkable linear size dependence. Compared with the previous TD-DFT and RCCSD(T) results, our estimates for MgC_2nH (n = 1–3) are in the best agreement with the available observed data of 2.83, 2.78, and 2.74 eV, respectively. In addition, the dissociation energies in MgC_2nH (n = 1–5) are also been evaluated.     

Force constant calculations of the LAM-1 modes ofn-alkanes and their substituted α-Monohalo- and α,ω-dihalo- species in their extended zig-zag planar configurations as found in urea inclusion compounds

Force constants were determined for the C₈, C₁₀, C₁₂ and C₁₄ series ofn-alkanes C _n H_{2n + 2} using an approximate SVFF calculation and observed LAM = 1 wave numbers. In this calculation the hydrogen atoms were neglected and only the carbon backbone chain and terminal atoms were considered; this was valid since only low-frequency vibrations were under consideration. Using force constant transfer, the wavenumbers of the LAM = 1 accordion modes for the analogous -C_n H_2n + 1 X and ,-C_n H_2nX₂ species, where X = C1, Br or I were calculated. For -chloroalkanes and ,-dichloroalkanes, them = 1 accordion modes are calculated to be in the 220–130 cm^–1 and 200–120 cm^–1 regions, respectively. For the bromo- and iodo-analogues them = 1 accordion modes are calculated to be in the 200–100 cm^–1, 150–90 cm^–1 and in the 170–100, 135–80 cm^–1 regions, respectively.     

Modeling the Proton Transport in Orthoperiodic and Orthotelluric Acids and Their Salts

Orthoperiodic and orthotelluric acids, their salts MIO₆H₄ (M = Li, Rb, Cs) and CsH₅TeO₆, and dimers of the salt · acid type are calculated within density functional theory B3LYP and basis set LanL2DZ complemented by the polarizationd,p-functions. According to calculations, the salt · acid dimerization is energetically favorable for compounds MIO₆H₄ · H₅IO₆ (M = Rb, Cs) and CsIO₆H₄ · H₆TeO₆. The dimerization energy is equal to 138–146 kJ mol^–1. With relatively small activation energies equal to 4 kJ mol^–1 (M = Li) and 11 kJ mol^–1 (M = Rb, Cs), possible is rotation of octahedron IO₆ relative to the M atom in monomers of salt molecules. The proton transfer along an octahedron occurs with activation energies of 63–84 kJ mol^–1. The activation energy for the proton transfer between neighboring octahedrons of the type salt · acid acid · salt equals 8–17 kJ mol^–1. Quantum-chemical calculations nicely conform to x-ray diffraction and electrochemical data.     

Theoretical studies of uracil–(H2O)n (n = 1–7) clusters by ab initio and ABEEMσπ/MM fluctuating charge model

Uracil–(H₂O)_n (n = 1–7) clusters were systemically investigated by ab initio methods and the newly constructed ABEEMσπ/MM fluctuating charge model. Water molecules have been gradually placed in an average plane containing uracil. The geometries of 38 uracil–water complexes were obtained using B3LYP/6-311++G** level optimizations, and the energies were determined at the MP2/6-311++G** level with BSSE corrections. The ABEEMσπ/MM potential model gives reasonable properties of these clusters when comparing with the present ab initio data. For interaction energies, the root mean square deviation is 0.96 kcal/mol, and the linear coefficient reaches 0.997. Furthermore, the ABEEMσπ charges changed when H₂O interacted with the uracil molecule, especially at the sites where the hydrogen bond form. These results show that the ABEEMσπ/MM model is fine giving the overall characteristic hydration properties of uracil–water systems in good agreement with the high-level ab initio calculations.     

Quantum-Chemical Evaluation of the Protolytic Properties of Thiophenols

The enthalpies and free energies of proton affinity (PA) were calculated by the PM3 method for 27 thiophenoxide anions. For thiophenols, linear correlation has been found between the PAs and the values of pK _a measured by different authors in aqueous acetone, aqueous ethanol, and methanol media. The dependences found permit one to predict a priori the protolytic properties of thiophenols. The differential effects of aqueous ethanol media with different ethanol contents on thiophenol acidity were evaluated based on quantum-chemical data.Original Russian Text Copyright © 2004 by A. N. Pankratov and A. V. Shalabai__________Translated from Zhurnal Strukturnoi Khimii, Vol. 45, No. 5, pp. 800–806, September–October, 2004.     

Structure of carbon fluorochlorides CF n Cl m (n,m ≤ 3) and their singly charged negative ions

The electronic and geometrical structures of carbon fluorochlorides with low coordination numbers (n 3) and their singly charged anions are calculated using the functional density method. The results of the calculations are used to evaluate the electron affinities (EA) of the neutral compounds and the first ionization potentials of the anions as well as the energies of fragmentation through different decay channels of both series. The adiabatic EA of carbon fluorochloride CF _k Cl _3–k is shown to be determined mainly by the presence of a CX₂ unit in these compounds. There are no monotonic changes in stability of either the neutral compounds withn = 3 or the anions withn = 2 or 3 upon successive substitution of one halogen by another.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 6, pp. 1044–1049, June, 1993.     

Theoretical studies of the benzene oxide—oxepin valence tautomerism

We have calculated the geometry and energy of the valence tautomers benzene oxide and oxepin using the semiempirical AM1 model and the 6–31G and 6–31G^* basis sets utilizing full geometry optimization. In the oxide the folding angle, the angle between the epoxide ring and the adjacent plane containing four carbon atoms, is about 106°. The carbon skeleton is almost planar, the folding angle, the angle between the two four-carbon atom planes being about 175°. In contrast, oxepin is found to have a marked boat-shaped structure with the corresponding and angles about 137° and 159°, respectively. The AM1, 6–31G, and 6–31G^* calculations give –11.4, –10.8, and –2.9 kcal mol^–1 for the energy change that accompanies the valence tautomerism, oxide-oxepin, compared to an experimental value of about +0.3 kcal mol^–1. Single point calculations of the energies at the 6–31 G^* geometry using Møller-Plesset perturbation theory to second order (MP2/6–31 G^*) and third order (MP3/6–31G^*) give E _T =+3.3 and +0.8 kcal mol^–1. The values for the energy change in the transfer of epoxide oxygen from ethylene oxide to benzene using AM1, 6–31G, and 6–31G^* are in good agreement, viz., +31.1, +34.5, and +33.6 kcal mol^–1, respectively. A large positive energy change is to be expected in view of the loss of benzene aromaticity.     

DFT study on the intermolecular interactions between Aun (n = 2–4) and thymine

Density functional B3LYP method with 6-31++G** basis set is applied to optimize the geometries of the luteolin, water and luteolin–(H₂O)_n complexes. The vibrational frequencies are also studied at the same level to analyze these complexes. We obtained four steady luteolin–H₂O, nine steady luteolin–(H₂O)₂ and ten steady luteolin–(H₂O)₃, respectively. Theories of atoms in molecules (AIM) and natural bond orbital (NBO) are used to investigate the hydrogen bonds involved in all the systems. The interaction energies of all the complexes corrected by basis set superposition error, are within −13.7 to −82.5 kJ/mol. The strong hydrogen bonding mainly contribute to the interaction energies, Natural bond orbital analysis is performed to reveal the origin of the interaction. All calculations also indicate that there are strong hydrogen bonding interactions in luteolin–(H₂O)_n complexes. The OH stretching modes of complexes are red-shifted relative to those of the monomer.     

Die Kristallstruktur der Tieftemperaturmodifikation von Ag8GeS6

The low temperature form of Ag₈GeS₆ (synthetic argyrodite) is orthorhombic, space group Pna2₁, witha=15.149 (1),b=7.476 (2),c=10.589 (1), andZ=4. The crystal structure has been determined from 2-circle diffractometer data by means of direct methods and refined toR=0.081 for 3431 intensities. The structure consists of slightly distorted isolated GeS₄ tetrahedra (mean Ge–S bond length 2.212 Å) and two further sulphur atoms without bonds to the germanium atoms. The GeS₄ tetrahedra and the S atoms are connected by the Ag atoms to form a three-dimensional framework. Three types of Ag coordination by S atoms are encountered. Three Ag positions have a strongly distorted tetrahedral environment, four Ag positions an approximately planar threefold coordination, while one Ag atom is almost linearly coordinated by two S atoms. The Ag–S distances are 2.56–2.94 Å, 2.49–2.76 Å, and 2.42–2.44 Å, resp. All Ag atoms have at least one near Ag neighbour between 2.93 and 3.11 Å.

---

---

Mit 3 Abbildungen     

Structures and bonding patterns of nanoannular carbon clusters (C4–C20) through AIM analyses

Structures, binding energies, harmonic frequencies, dipole moments, HOMO–LUMO energy gaps and particularly atoms in molecules (AIM) analyses of some nanoannular carbon clusters (C₄–C₂₀) are investigated at B3LYP/6-31+G(d) level of theory. No correlation is found by plotting the calculated binding energies as a functional number of carbon atoms of carbon clusters. The calculated binding energies sharply increase from C₄ to C₁₀ while slowly from C₁₀ to C₂₀. The binding energies of C_4n+2 clusters including C₆, C₁₀, C₁₄, and C₁₈ have a clear increase when compared with others indicating their aromatic characters which is confirmed by results of HOMO–LUMO energy gaps and geometrical parameters. AIM analyses show that most of our carbon clusters are topologically normal (non-conflict) with stable structures. Nevertheless, the topological networks of small antiaromatic rings, C₄ and C₈, at their equilibrium geometries may change via molecular vibrations. The existence of straight bond paths in 3D molecular graphs of carbon clusters with n > 10 implies that ring strains are decreased as the ring sizes grow. Except for C₄ and C₈, the ellipticity values for the remaining carbon clusters are small indicating that the C–C bond is conserved in these clusters. Dipole moments of even-numbered structures are negligible, whereas odd-numbered ones have μ values of 0.09−0.73 D.     

Formation of mixed aryl–, alkyl–lithium aggregates in the heteroatom assisted lithiation of α,α′-dialkyl substituted 1,3-bis[(dimethylamino)methyl]benzene

The heteroatom assisted lithiation of 1,3-bis[1-(dimethylamino)ethyl]benzene with n-BuLi afforded 2,6-bis[1-(dimethylamino)ethyl]phenyllithium. An X-ray crystal structure determination revealed a dimeric aggregate in which the four benzylic chiral centers are identical, pointing to stereoselective crystallization. In contrast, reaction of 1,3-bis[1-(dimethylamino)propyl]benzene with n-BuLi afforded a dimeric aggregate comprising the parent lithiated compound and n-BuLi in a 1:1 molar ratio. The four Li atoms and the four bridging carbon atoms are arranged in a unique ladder-type C–Li₂–C₂′–Li₂–C framework.     

Theoretical study on the structures and properties of LiCn, , and (n = 1–10) clusters

The lithium-doped carbon clusters LiC_n, , and (n = 1–10) have been investigated systemically with density functional theory (DFT) method at the B3LYP/6-311+G* level. According to the total energies of different kinds of isomers, the LiC_n, , and (n = 1–10) clusters have Li-terminated linear ground states structures, except for LiC₂, LiC₃, , and (n = 4–6). The incremental binding energies are evaluated to elucidate the stabilities of the clusters with different numbers of carbon atoms for neutral molecules, cations, and anions, respectively. Clear even–odd alternation effects are observed for the stability of the cationic clusters and anionic clusters, while for neutral LiC_n clusters the alternation effect is less pronounced. Similarly, the ionization potentials and electron affinities of LiC_n also express an obvious parity alternation. In addition, the most favorable dissociation channels are acquired according to the fragmentation energies accompanying various pathways.     

Characteristic Features of the Rhodium Stereochemistry in the Structure of Carbonyls and Organometallic Compounds

Characteristics of the Voronoi–Dirichlet polyhedra (VDP) were used to carry out crystal chemical analysis of 160 compounds containing 308 RhC _n or RhC _n Rh _m coordination polyhedra. The volume of the rhodium VDP virtually does not depend on the coordination number (C.N.), which varies from 4 to 12, in spite of the pronounced variation of the Rh–C (1.72–2.83 Å) and Rh–Rh (2.55–2.97 Å) bond lengths. It is shown that the VDP parameters allow one to estimate quantitatively the main features of rhodium stereochemistry irrespective of the nature (Rh or C) and the number (n varies from 3 to 12 and m varies from 0 to 6) of atoms in the first coordination sphere.     

Effect of the potential of an external electron donor on C<Subscript>2</Subscript>H<Subscript>2</Subscript> reduction catalyzed by the nitrogenase active center (FeMoco) isolated from the enzyme

The effect of potential value and chemical properties of an external electron donor on C₂H₂ reduction catalyzed by nitrogenase active center (cluster [(₆-N)Fe₇MoS₉·homocitrate] FeMoco isolated from the enzyme) has been investigated in the presence of proton donors of different acidity. The temperature—reaction rate dependences of these reactions have been studied. It has been shown that the rate-limiting steps of the reactions differ depending on the proton donor used. When thiophenol or water are used as proton donors, and electrochemical step — the electron transfer from cathode to adsorbed catalytic cluster — has been found to be a rate-limiting one. The effective activation energy of ethane formation as a product of four-electron C₂H₂ reduction is found to be 1.5 times lower than that of ethylene, namely, 13 kcal mol^–1. When stronger acid, pentafluorothiophenol, is used as a proton donor, the chemical step of intramolecular rearrangement of the catalyst—substrate complex taking place in solution becomes a rate-limiting one. The effective activation energies of both ethylene and ethane become equal to 32 kcal mol^–1.Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 8, pp. 1583–1591, August, 2004.     

A density functional theory study on boundary of “superreduced” transition metal carbonyl anions [M(CO)n](M=Cr, n=5, 4, 3, z=2, 4, 6; M=Mn, n=5, 4, 3, z=1, 3, 5; M=Fe, n=4, 3, 2, z=2, 4, 6; M=Co, n=4, 3, 2, z=1, 3, 5)

A nonlocal density functional theory (DFT) method has been applied to the calculations on optimized geometry, Mulliken atomic net charges and interatomic Mulliken bond orders as well as total bonding energies (E) in the binary transition metal carbonyl anions with different reduced states [M(CO)_n]^z− (M=Cr, n=5, 4, 3, z=2, 4, 6; M=Mn, n=5, 4, 3, z=1, 3, 5; M=Fe, n=4, 3, 2, z=2, 4, 6; M=Co, n=4, 3, 2, z=1, 3, 5). For comparison of relative stability, a relative stabilization energy D is defined as D=E([M(CO)_n]^z−)−nE(CO). The calculated C–O distances are lengthened monotonously with the increase of the anionic charge, but the M–C distances are significantly lengthened only in the higher reduced states. The relative stabilization energy calculated is a considerable negative value in the lower reduced states, but a larger positive value in the higher reduced states. The DFT calculations show that with the increase of the anionic charge, the Mulliken net charges on the M, C, and O atoms all increase, however, an excess of the anionic charge is mainly located at the central metal atom. The calculated C–O Mulliken bond orders decrease consistently with the increase of the anionic charge, but the M–C bond orders exhibit an irregular behavior. However, the total bond orders calculated clearly explain the higher reduced states to be considerably unstable. From analysis of the calculated results, it is deduced that the stability of the binary transition metal carbonyl anions [M(CO)_n]^z− studied are associated with the coordination number n and the anionic charge z, further, it is possible for the anions studied to be stable if n≥z, conversely, it is impossible when n<z.