Hypercoordinate atoms of second-row elements in dodecahedrane endohedral complexes

The energy characteristics and geometric parameters of the dodecahedrane endohedral complexes X@C₂₀H₂₀ (X = C⁴⁻, N³⁻, O²⁻, F⁻, Ne) were studied by the density functional theory B3LYP method with the 6-311G(d,p), 6-311+G(d,p), and 6-311G(df,p)) basis sets. In all structures the central atoms X are characterized by a coordination number of 20. The energy of formation of the complexes decreases in the order X = C⁴⁻, N³⁻, O²⁻, F⁻, Ne. The coordination number of the central atom remains unchanged upon adding Li⁺ counterions to anionic systems. Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 5, pp. 824–830, May, 2007.     

氟代乙烯阳离子的理论研究

用B3LYP和MP2方法及6-31G(d, p)、6-31＋G(d, p)、6-311G(d, p)和6-311＋G(d, p) 基组，对六种氟代乙烯阳离子做了理论研究，优化了它们的基电子态的结构，计算了对应分子的垂直电离势（VIP）和绝热电离势(AIP).结果表明，与具有非平面结构的乙烯阳离子不同，六种氟代乙烯阳离子都只具有平面结构；与分子结构相比，离子结构的C－C键增长， C－F键缩短， CCF键角变小. 自然布居分析计算表明，这些离子的正电荷主要分布在与F原子相连的C原子和各H原子上. B3LYP/6-311＋G(d, p) 级别上计算的各分子的VIP和AIP值和实验值符合得很好. 使用含弥散基函数的基集可以明显提高这类分子的电离势的计算精度.     

Synthesis of Cationic R2P5+ Cages and Subsequent Chalcogenation Reactions

Cationic R₂P₅⁺ cage compounds ( 1 ⁺) have been synthesized by the stoichiometric reaction of R₂PCl, GaCl₃ and P₄. The reaction conditions depend on the substituent R. Alkyl‐substituted derivatives ( 1 a – 1 d [GaCl₄]) are best synthesized under solvent‐free conditions, whereas aryl‐substituted derivatives ( 1 e – 1 h [GaCl₄]) are formed in C₆H₅F. All compounds have been prepared on a multi‐gram scale in good to excellent yields and have been fully characterized with an emphasis on ³¹P NMR spectroscopy in solution and single‐crystal structure determination. Subsequent chalcogenation reactions of cations R₂P₅⁺ ( 1 a ⁺, 1 e ⁺) and trication Ph₆P₇³⁺ ( 3 ³⁺) with elemental sulfur (α‐S₈) or grey selenium (Se_grey) yielded a series of unique polyphosphorus–chalcogen cations ( 4 a ⁺, 4 e ⁺, 5 a ⁺, 6 ²⁺ and 7 ²⁺), possessing nortricyclane‐type molecular structures. An in‐depth study of the ³¹P{¹H} and ⁷⁷Se NMR spectroscopic parameters is presented, and correlations between the substitution pattern and the observed structural features have been investigated in detail.     

Ab initio molecular orbital analysis of acetaldehyde dimers.: Thermodynamic properties

The surface of the acetaldehyde dimer was studied by MP2/6-31+G^* ab initio calculations and some of the minima thus located were used in subsequent optimizations and single-point calculations at the B3-LYP/6-311+G(2d,p), MP2/6-311+G(2d,p), MP2/aug-cc-pDVZ and MP2/aug-cc-pTVZ levels. An overall of six minima in two distinct groups were detected, one group consisting of planar configurations and the other of more stable, non-planar (spatial) configurations of the two monomers. Some of the dimer thermodynamic properties were calculated from its harmonic frequencies and a comparison between the experimental values for the thermodynamic properties and the calculated values is discussed.     

Inductive effects in radicals calculated from DFT energies; substituted bicyclo[2.2.2]octan-1-yloxy radicals

Energies of a series of 4-substituted 1-oxybicyclo[2.2.2]octan-1-yloxy radicals with 18 various substituents were calculated within the framework of the DFT theory at the levels UB3LYP/6-311+G(d,p)//UB3LYP/6-311+G(d,p) and UB3LYP/6-311++G(2df,p)//UB3LYP/6-311+G(d,p) and compared with similar series of the parent alcohols, their deprotonated and protonated forms calculated at the levels B3LYP/6-311+G(d,p)//B3LYP/6-311+G(d,p) and B3LYP/6-311++G(2df,p)//B3LYP/6-311+G(d,p). The two levels are of the same performance and both are sufficient for molecules of this type according to comparison with scarce experimental gas-phase acidities and basicities. The substituent effects were analyzed in terms of isodesmic equations. In addition to strong dependence on the substituent inductive effect, a slight dependence on the electronegativity of the first atom of the substituent was proven in certain cases. In all aspects, there is no qualitative difference between the effects on radicals and on similar closed shell species. Radicals behave as slightly electron deficient; the substituent effect is weaker than that on the ions but stronger than on neutral molecules.     

Investigation on transition States of [Alanine + M2+] (M = Ca,Cu, and Zn) complexes: A quantum chemical study

The ground state geometries of [Alanine (Ala) + M²⁺] [M = Ca, Cu, and Zn) complexes were calculated in gas phase at B3LYP/6‐311++G(d,p) level of theory. Transition states (TSs) between different stable conformers of [Ala + M²⁺] complexes were also calculated. Among the different [Ala + M²⁺] complexes, the complex where metal cations coordinated to carboxylate group (? COO^?) is found to be energetically most favorable. To calculate TSs, the ground state structures of any two conformers of [Ala + M²⁺] complexes were used. The ground state energies of two stable conformers and their TS structures were used to calculate the activation energy. The reactivity of different conformers of [Ala + M²⁺] complexes have been discussed in terms of energy difference between their highest occupied molecular orbital and lowest unoccupied molecular orbital. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Aromaticity in Group 14 homologues of the cyclopropenylium cation

The nature of the bonding and the aromaticity of the heavy Group 14 homologues of cyclopropenylium cations E₃H₃⁺ and E₂H₂E′H⁺ (E, E′=C–Pb) have been investigated systematically at the BP86/TZ2P DFT level by using several methods. Aromatic stabilization energies (ASE) were evaluated from the values obtained from energy decomposition analysis (EDA) of charged acyclic reference molecules. The EDA‐ASE results compare well with the extra cyclic resonance energy (ECRE) values given by the block localized wavefunction (BLW) method. Although all compounds investigated are Hückel 4n+2 π electron species, their ASEs indicate that the inclusion of Group 14 elements heavier than carbon reduces the aromaticity; the parent C₃H₃⁺ ion and Si₂H₂CH⁺ are the most aromatic, and Pb₃H₃⁺ is the least so. The higher energies for the cyclopropenium analogues reported in 1995 employed an isodesmic scheme, and are reinterpreted by using the BLW method. The decrease in the strength of both the π cyclic conjugation and the aromaticity in the order C?Si>Ge>Sn>Pb agrees reasonably well with the trends given by the refined nucleus‐independent chemical shift NICS(0)_πzz index.     

Lewis Base Free Oxophosphonium Ions: Tunable,Trigonal‐Planar Lewis Acids

Oxophosphonium ions (R₂P=O)⁺ are fascinating chemical intermediates related to the well‐known acylium cations (RC=O)⁺, and comprise a tricoordinate phosphorus(V) center with a phosphorus–oxygen double bond. Here, we report the synthesis of two oxophosphonium ions stabilized by bulky imidazolin‐2‐imine and imidazolin‐2‐olefin substituents attached to phosphorus. The novel species were characterized by NMR spectroscopy and single‐crystal X‐ray diffraction analysis, and the bonding situation was probed by DFT calculations. Determination of the acceptor number and the fluoride ion affinity revealed that the choice of the substituents has a strong influence on the electrophilicity of the phosphorus center. Additionally, the formation of Lewis base adducts with pyridine derivatives and the reactivity with isopropyl alcohol was explored.     

The acidity of analogous ammonium cations: A description of the solvent effect through the attainment of hydration clusters using the AGOA methodology

The optimized geometries of ethylaminium (1), ethyldiaminium (2), n-butylaminium (3) and tri-ethylaminium (4) cations were calculated using the B3LYP/6-311++G(d,p) theoretical level. Subsequently, the hydration clusters of these structures were explored using the AGOA methodology. Analysing the molecular electrostatic potential of the solute and using the TIP4P model to orient the water molecules around the ammonium cations, systematic increases were observed in hydration energy, from 1 to 4. This result does not correlate with the experimental acidity of 1–4, suggesting that the AGOA methodology is not efficient in modelling cationic structures in an aqueous medium. The AGOA hydration clusters were, therefore, re-optimized (AGOA-OPT) using the B3LYP/6-311++G (d,p) theoretical level. This new AGOA-OPT procedure revealed the hydration energy of the 1, 2, 3 and 4 cations to be in accordance with the experimental acidity.     

Quantum chemical study of the interaction between selenium analog of methimazole as an anti-thyroid drug and metal cations

Energetic and structural properties of complexes formed from interaction between selenium analog of methimazole (MSeI) as an anti-thyroid drug and M^z+ (Li⁺, Na⁺, K⁺, Be²⁺, Mg²⁺ and Ca²⁺) cations have been investigated using B3LYP, M062X, PBE1PBE, and MP2 methods with 6-311++G(d,p) and 6-311++G(2d,2p) basis sets. Two planar and perpendicular complexes were predicted from interaction of MSeI and M^z+ cations. From the Gibbs free energy difference between the planar and perpendicular forms of MSeI–M^z+ complexes, it is found that the perpendicular forms are the predominant ones. In addition, the comparison of interaction energies shows that the order of energies increases in the following order: K⁺ < Na⁺ < Li⁺ < Ca²⁺ < Mg²⁺ < Be²⁺. The results of natural bond orbital analysis showed that the charge transfer occurs from MSeI to metal cations. The atom in molecule analysis shows that the charge density and its Laplacian at the Se–M^z+ bond critical point of the MSeI–M²⁺ complexes are greater than the MSeI–M¹⁺ ones. Also, it was revealed that the Se–M^z+ interactions in perpendicular complexes of alkali and alkaline metal cations are electrostatic and partially covalent in nature, respectively.     

Conjugation in radical cations of benzene, ethylene, and acetylene derivatives. Nonuniversality of the resonance parameters of group 14 organoelement substituents

The first vertical ionization potentialsl ₁ of molecules R_πX (X=Ph, H₂C=CH, and HC≡C) depend on the joint influence of the inductive, resonance, and polarizability effects of substituents X, which are characterized by parameters σ₁, σ_R+, and σ_α, respectively. The mechanism of conjugation in radical cations formed upon ionization of R_πX is changed as compared to neutral R_πX molecules, while the substituent X becomes polarized. The conjugation and polarizability effects are strenthened in the sequence Ph < H₂C=CH <HC≡C as R_π changes from Ph to H₂C=CH and HC≡C. The σ_R+ parameters of Si-, Ge-, and Sn-containing substituents X are dependent on the type of R_π but are connected by linear dependences in the series of benzene, ethylene, and acetylene derivatives. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 8, pp. 1481–1486, August, 1998.     

Isomerizations of CH3SO2: Quantum Chemistry and Topological Study

The structures, stabilities and the isomerization reactions of CH₃SO₂ isomers in a doublet electronic state have been studied at B3LYP/6‐311+ +G (d,p), MP2/6‐311++G (d,p) and CCSD(T)/6‐311++G (d,p) levels. The three different levels of calculation give the similar results: thirteen minimum isomers were located and they were connected by eleven transition states. Among the thirteen isomers, cis‐CH₃OSO, trans‐CH₃OSO and CH₃SO₂ are the most stable species, and they should be detected easily in experiment. This is well consistent with the experimental result. These isomers could isomerize to each other by chemical bond vibration, chemical bond rotation and atom migration. The non‐planar ring structure transition state (STS), which was found in this paper, extended the concept of ring STS to the non‐planar systems.     

Ionization potentials of non-metal monosulfides. Conjugation in radical cations containing Group IV elements

The first vertical ionization potentialsI(n_s) of 69 monosulfides XSY (X, Y=H, Hal, organic, or heteroorganic substituent) are related to the inductive σ_I resonance (σ _R ⁺ ) and polarizability (σ_α) constants of the substituents by dependences of theI(n_S)=a+bΣσ_I+bΣσ_R+bΣσ_α type. TheI(n_s) values are also affected by hyperconjugation which increases on going from XSH to XSY (Y≠H) compounds. The first calculations of the σ _R ⁺ parameters characterizing the conjugation of Si-, Ge-, Sn-, and Pb-containing substituents with the S^.+ radical cation center are reported. The reasons for weakening of resonance donor properties of heteroorganic substituents of the +M-type in the systems studied as compared to those of the same substituents in the corresponding aromatic radical cations are considered. Translated fromIzvestiya Akademii Nauk. Seriya Khmicheskaya, No. 1, pp. 25–31, January, 2000.     

A DFT study of hydride transfers to the carbonyl oxygen of DDQ

Hydride‐transfer reactions between benzylic substrates and 2,3‐dichloro‐5,6‐dicyano‐1,4‐benzoquinone (DDQ) were investigated by DFT (density functional theory) calculations. The lowest unoccupied molecular orbital of DDQ has the largest extension on two carbonyl oxygens, which comes from two‐step mixing of antisymmetric orbitals of fragment π MOs. Transition‐state (TS) geometries and activation energies of reactions of four benzylic substrates R²? CH₂‐para‐C₆H₄? R¹ (R¹, R² = H and/or OCH₃) with DDQ were calculated. M06‐2X/6‐311(+*)G* was found to be a practical computational method, giving energies and geometries similar to those of M06‐2X/6‐311++G(3df,2pd) and wB97xD/6‐311++G(3df,2pd). For toluene (R¹ = R² = H), an initiation‐propagation model was suggested, and the calculated kinetic isotope effect k(H)/k(D) = 5.0 with the tunnel correction at the propagating step is in good agreement with the experimental value 5.2. A reaction of para‐MeO? C₆H₄? CH₂(OMe) + DDQ + (H₂O)₁₄ → para‐MeO? C₆H₄? C(?O)H + HOMe + DDQH₂ + (H₂O)₁₃ was investigated by M06‐2X/6‐311(+*)G*. Four elementary processes were found and the hydride transfer (TS1) is the rate‐determining step. The hydride transfer was promoted by association with the water cluster. The size of the water cluster, (H₂O)_n, at TS1 was examined. Three models of n = 14, 20, and 26 were found to give similar activation energies. Metal‐free neutral hydride transfers from activated benzylic substrates to DDQ were proposed to be ready processes both kinetically and thermodynamically. © 2015 Wiley Periodicals, Inc.     

Ionization potentials of non-metal monosulfides. Conjugation in radical cations containing Group IV elements

The first vertical ionization potentialsI(n_s) of 69 monosulfides XSY (X, Y=H, Hal, organic, or heteroorganic substituent) are related to the inductive σ_I resonance (σ _R ⁺ ) and polarizability (σ_α) constants of the substituents by dependences of theI(n_S)=a+bΣσ_I+bΣσ_R+bΣσ_α type. TheI(n_s) values are also affected by hyperconjugation which increases on going from XSH to XSY (Y≠H) compounds. The first calculations of the σ _R ⁺ parameters characterizing the conjugation of Si-, Ge-, Sn-, and Pb-containing substituents with the S^.+ radical cation center are reported. The reasons for weakening of resonance donor properties of heteroorganic substituents of the +M-type in the systems studied as compared to those of the same substituents in the corresponding aromatic radical cations are considered. Translated fromIzvestiya Akademii Nauk. Seriya Khmicheskaya, No. 1, pp. 25–31, January, 2000.     

戊烯自由基阳离子的密度泛函理论研究

使用密度泛函理论B3LYP方法和6-31G(d,p)、6-31＋G(d,p)、6-311G(d,p)及6-311＋G(d,p)基组，分别对2-C5H10＋和1-C5H10＋的各种构象进行了几何构型优化，并用B3LYP/6-311G(d,p)进行了频率分析计算.计算预言1-C5H10＋具有非平面构型，与以往报导的从头算计算结论相反.在两个自由基阳离子的各种构象的B3LYP几何构型上，进行了B3LYP和UMP2(full)方法的超精细偶合常数计算，得到了比以往更好的结果.     

Are They Linear,Bent, or Cyclic? Quantum Chemical Investigation of the Heavier Group 14 and Group 15 Homologues of HCN and HNC

The singlet potential‐energy surface (PES) of the system involving the atoms H, X, and E (the (H, X, E) system) in which X=N–Bi and E=C–Pb has been explored at the CCSD(T)/TZVPP and BP86/TZ2P+ levels of theory. The nature of the X? E bonding has been analyzed with charge‐ and energy‐partitioning methods. The calculations show that the linear isomers of the nitrogen systems lin ‐HEN and lin ‐HNE are minima on the singlet PES. The carbon compound lin ‐HCN (HCN=hydrogen cyanide) is 14.9 kcal mol^?1 lower in energy than lin ‐HNC but the heavier group 14 homologues lin ‐HEN (E=Si–Pb) are between 64.8 and 71.5 kcal mol^?1 less stable than the lin ‐HNE isomers. The phosphorous system (H, P, E) exhibits significant differences concerning the geometry and stability of the equilibrium structures compared with the nitrogen system. The linear form lin ‐HEP of the former system is much more stable than lin ‐HPE . The molecule lin ‐HCP is the only minimum on the singlet PES. It is 78.5 kcal mol^?1 lower in energy than lin ‐HPC , which is a second‐order saddle point. The heavier homologues lin ‐HPE , in which E=Si–Pb, are also second‐order saddle points, whereas the bent ‐HPE structures are the global minima on the PES. They are between 10.3 (E=Si) and 36.5 kcal mol^?1 (E=Pb) lower in energy than lin ‐HEP . The bent ‐HPE structures possess rather acute bending angles H‐P‐E between 60.1 (E=Si) and 79.7° (E=Pb). The energy differences between the heavier group 15 isomers lin ‐HEX (X=P–Bi) and the bent structures bent ‐HXE become continuously smaller. The silicon species lin ‐HSiBi is even 3.1 kcal mol^?1 lower in energy than bent ‐HBiSi . The bending angle H‐X‐E becomes more acute when X becomes heavier. The drastic energy differences between the isomers of the system (H, X, E) are explained with three factors that determine the relative stabilities of the energy minima: 1) The different bond strength between the hydrogen bonds H? X and H? E. 2) The electronic excitation energy of the fragment HE from the X ²Π ground state to the ⁴Σ^? excited state, which is required to establish a E≡X triple bond in the molecules lin ‐HEX . 3) The strength of the intrinsic X? E interactions in the molecules. The trends of the geometries and relative energies of the linear, bent, and cyclic isomers are explained with an energy‐decomposition analysis that provides deep insight into the nature of the bonding situation.     

Effect of metal cations [Li+, Na+, K+, Be2+, Mg2+, and Ca2+] on the structure of 2‐(3′‐hydroxy‐2′‐pyridyl)benzoxazole: A theoretical investigation

Density functional theory calculations were performed at the B3LYP/6‐311++G(d,p) level to systematically explore the geometrical multiplicity and binding strength for the complexes formed by alkaline and alkaline earth metal cations, viz. Li⁺, Na⁺, K⁺, Be²⁺, Mg²⁺, and Ca²⁺ (Mⁿ⁺, hereinafter), with 2‐(3′‐hydroxy‐2′‐pyridyl)benzoxazole. A total of 60 initial structures were designed and optimized, of which 51 optimized structures were found, which could be divided into two different types: monodentate complexes and bidentate complexes. In the cation‐heteroatom complex, bidentate binding is generally stronger than monodentate binding, and of which the bidentate binding with five‐membered ring structure has the strongest interaction. Energy decomposition revealed that the total binding energies mainly come from electrostatic interaction for alkaline metal ion complexes and orbital interaction energy for alkaline earth metal ion complex. In addition, the electron localization function analysis show that only the Be? O and Be? N bond are covalent character, and others are ionic character. © 2012 Wiley Periodicals, Inc.     

氟代吡啶阳离子的理论研究

Cations of fluorinated pyridines（pentafluoropyridine,2,6-difluoropyridine,and 2-fluoropyridine）have been studied by using density functional B3LYP method in conjunction with 6-31G（d,p）,6-311G（d,p）,6-31+G（d,p）,and 6-311+G（d,p）basis sets. B3LYP geometry optimization and frequency analysis calculations indicate that the pentafluoropyridine cation,2,6-difluoropyridine cation,and 2-fluoropyridine cation have C2v,C2v,and Cs structures in the 2A2,2A2,and 2A" ground states,respectively. The calculated geometries of the cations and the parent molecules were compared. The natural population analysis calculations at the B3LYP level with different basis sets were performed on the three cations and the three parent molecules. The isotropic hyperfine coupling constants in the three cations（radicals）were calculated. The vertical and adiabatic ionization potential（VIP and AIP）values of the pentafluoropyridine,2,6-difluoropyridine,and 2-fluoropyridine molecules were calculated by using the B3LYP method,and the calculated VIP values are in excellent agreement with experiment.     

Predominantly sigma-bonded metal carbonyl cations (sigma-carbonyls): new synthetic,structural, and bonding concepts in metal carbonyl chemistry

Predominantly σ‐bonded metal carbonyl cations (σ‐carbonyls) are conveniently generated in the Lewis superacid SbF₅ or the conjugate Brønsted–Lewis superacid HF? SbF₅, primarily by solvolytic or reductive carbonylations. Thermally stable salts are formed with the fluoroantimonate(V ) ions [SbF₆]^? and [Sb₂F₁₁]^?. The salts are characterised by analytical, structural, spectroscopic and computational methods. Most homoleptic carbonyl cations have very regular geometries, comensurate with their d‐electron configurations: linear (d¹⁰), square planar (d⁸) or octahedral(d⁶). The cations with metals in oxidation states of +2 or +3 are termed “superelectrophilic”. Extended molecular structures form by significant interionic C? F contacts with electrophilic carbon as acceptor. To account for all experimental observations, a conceptually simple synergetic bonding model is proposed. An outlook at anticipated future developments based on very recent results is provided.