Theoretical predictions of the structure,gas‐phase acidity,and aromaticity of tetrathiosquaric acid

Results of ab initio self‐consistent‐field and density functional theory calculations of the gas‐phase structure, acidity (free energy of deprotonation, ΔG⁰), and aromaticity of tetrathiosquaric acid (3,4‐dithiohydroxy‐3‐cyclobutene‐1,2‐dithione, H₂C₄S₄) are reported. The global minimum found on the potential energy surface of tetrathiosquaric acid presents a planar conformation. The ZZ isomer was found to have the lowest energy among the three planar conformers and the ZZ and ZE isomers are very close in energy. The optimized geometric parameters exhibit a bond length equalization relative to reference compounds, cyclobutanedithione, and cyclobutenedithiol. The computed aromatic stabilization energy by homodesmotic reaction is −18.4 (MP2(fu)/6‐311+G**//RHF/6‐311+G**) and −15.1 kcal/mol (B3LYP//6‐311+G**// B3LYP/6‐311+G**). The aromaticity of tetrathiosquaric acid is indicated by the calculated diamagnetic susceptibility exaltation (Λ) −11.77 (CSGT(IGAIM)‐RHF/6‐311+G**// RHF/6‐311+G**) and −18.08 (CSGT(IGAIM)‐B3LYP/6‐311+G**// B3LYP/6‐311+G**). Thus, tetrathiosquaric acid fulfils the geometric, energetic and magnetic criteria of aromaticity. The most reliable theoretical gas‐phase acidities are $\Delta G^{0}_{1(298\mathrm{K})}=303.7$ and $\Delta G^{0}_{2(298\mathrm{K})}=394.1$ kcal/mol. Hence, tetrathiosquaric acid is a stronger acid than squaric acid (3,4‐dihydroxy‐3‐ cyclobutene‐1,2‐dione, H₂C₄O₄). Comparisons of the computed results of tetrathiosquaric acid with squaric acid have also been made. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 78: 443–449, 2000     

Carbon,nitrogen, and oxygen hypercoordination in half-sandwich and sandwich structures

The structures and stabilities of a number of neutral and charged half-sandwich (pyramidal) and sandwich compounds, which obey the “electron octet” rule and contain hypercoordinate carbon, nitrogen, and oxygen atoms, were studied by ab initio MP2(full)/6-311+G** and density functional B3LYP/6-311+G** methods. Introduction of lithium counterions or bridging hydrogen atoms can provide an additional stabilization of non-classical systems with hypercoordinate centers.__________Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 3, pp. 527–540, March, 2005.     

AB initio study of thermodynamic stability and structure of cage molecules B4N4H8 and Be4O4H8

Stable molecular structures of heterocubane systems B₄N₄H₈ 2 and Be₄O₄H₈, isoelectronic to the cubane molecule, are investigated by ab initio (RHF/6-31G**, MP2(full)/6-31C**, and MP2(full)/6-311+ + G**) methods and are shown to be highly thetmodynamically stable. Decomposition of structure 2 into two 1,3,2,4-diazadiboroethidine molecules 6 or four iminobomne NBNH molecules 11 is an endothermal process taking 10.1 (RHF/6-31G**), 39.6 (MP2(full)/6-31G**) kcaUmole and 140.6 (RHF/6-31G**), 161.4 (MP2(full)/6-31G**) kcal/mole, respectively. Decomposition of structure 3 into two 1,3,2,4-dioxydi-beryllothidine molecules 12 or four molecules 13 is also an endothermal reaction taking 22.1 (RHF/6-31G**), 39.8 (MP2(full)/6-31G**) kcal/mole and 127.1 (RHF/6-31G**), 155.2 (MP2(full)/631G**) kcal/mole, respectively. The geometrical characteristics of simple molecules BeH₂ 15, Be₂ 16 and 17, Be₂H₂ 18, Be₂H₄ 19, BeO 20, and Be₂O₂ 21 are calculated. Translated from Zhumal Struktumoi Khim ii, Vol. 41, No. 1, pp. 3-13, January–February, 2000     

Effect of the counterion on the steric and electronic structure of pyrylium cation

The effect of counterions on the stabilization of aromatic pyrylium cation, a heterocyclic analog of benzene, was examined by ab initio [MP2(full)/6-311++G**] method and density functional theory [B3LYP/6-311++G(d,p)(6d,10f)]. Structures of the charge-transfer molecular complexes formed by interaction of pyrylium cation with counterion are predicted. The principal role of the counterion in the stabilization of the pyrylium cation aromatic system is demonstrated. The stabilization results from both electron density transfer and covalent bonding.     

Theoretical predictions of the structure, gas-phase acidity and aromaticity of 1,2-diseleno-3,4-dithiosquaric acid

Results of ab initio self-consistent-field (SCF) and density functional theory (DFT) calculations of the gas-phase structure, acidity (free energy of deprotonation, ΔG^o), and aromaticity of 1,2-diseleno-3,4-dithiosquaric acid (3,4-dithiohydroxy-3-cyclobutene-1,2-diselenone, H₂C₄Se₂S₂) are reported. The global minimum found on the potential energy surface of 1,2-diseleno-3,4-dithiosquaric acid presents a planar conformation. The ZZ isomer was found to have the lowest energy among the three planar conformers and the ZZ and ZE isomers are very close in energy. The optimized geometric parameters exhibit a bond length equalization relative to reference compounds, cyclobutanediselenone, and cyclobutenedithiol. The computed aromatic stabilization energy (ASE) by homodesmotic reaction (Eq 1) is −20.1 kcal/mol (MP2(fu)/6-311+G** //RHF/6-311+G**) and −14.9 kcal/mol (B3LYP//6-311+G**//B3LYP/6-311+G**). The aromaticity of 1,2-diseleno-3,4-dithiosquaric acid is indicated by the calculated diamagnetic susceptibility exaltation (Λ) −17.91 (CSGT(IGAIM)-RHF/6-311+G**//RHF/6-311+G**) and −31.01 (CSGT(IGAIM)-B3LYP/6-311+G**//B3LYP/6-311+G**). Thus, 1,2-diseleno-3,4-dithiosquaric acid fulfils the geometric, energetic and magnetic criteria of aromaticity. The calculated theoretical gas-phase acidity is ΔG^o _1(298K)=302.7 kcal/mol and ΔG^o _2(298K)=388.4 kcal/mol. Hence, 1,2-diseleno-3,4-dithiosquaric acid is a stronger acid than squaric acid(3,4-dihydroxy-3-cyclobutene-1,2-dione, H₂C₄O₄). Received: 11 April 2000 / Accepted: 7 July 2000 / Published online: 27 September 2000     

Predicting a New Δ-Proton Sponge-Base of 4,12-Dihydrogen-4,8,12-triazatriangulene through Proton Affinity,Aromatic Stabilization Energy,and Aromatic Magnetism

Herein, we report designing a new Δ (delta-shaped) proton sponge base of 4,12-dihydrogen-4,8,12-triazatriangulene (compound 1 ) and calculating its proton affinity (PA), aromatic stabilization, natural bond orbital (NBO), electron density ρ(r), Laplacian of electron density ∇²ρ(r), (2D-3D) multidimensional off-nucleus magnetic shielding (σ_zz(r) and σ_iso(r)), and scanning nucleus-independent chemical shift (NICS_zz and NICS). Density functional theory (DFT) at B3LYP/6-311+G(d,p), ωB97XD/6-311+G(d,p), and PW91/def2TZVP were used to compute the magnetic shielding variables. In addition, relevant bases like pyridine, quinoline, and acridine were also studied and compared. The protonation of compound 1 yields a highly symmetric carbocation of three Hückel benzenic rings. Comparing our findings of the studied molecules showed that compound 1 precedes others in PA, aromatic isomerization stabilization energy, and basicity. Therefore, the basicity may be enhanced when a conjugate acid gains higher aromatic features than its unprotonated base. Both multidimensional σ_zz(r) and σ_iso(r) off-nucleus magnetic shieldings outperformed electron-based techniques and can visually monitor changes in aromaticity that occur by protonation. The B3LYP/6-311+G(d,p), ωB97XD/6-311+G(d,p), and PW91/def2TZVP levels showed no significant differences in detailing isochemical shielding surfaces.     

Tetra- and hexaatomic cyclic clusters of main-group elements (XY)2 and (XY)3: an ab initio and density functional study

The electronic and molecular structure of planar (cyclic and linear) tetra- and hexaatomic clusters (XY) _n (XY = CC, BN, BeO, LiF; n = 2, 3) was studied using the ab initio CCD(full)/6-311+G** method and density functional approach (B3LYP/6-311+G**). The stability of cyclic clusters C₆, B₃N₃, and Be₃O₃ with D_3h symmetry is mainly determined by the aromaticity of their -electron systems.     

Can Substituted Cyclopentadiene Become Aromatic or Antiaromatic?

Cyclopentadiene derivatives with electronegative (F, Cl) or electropositive (H(3)Si, Me(3)Si) bis-5,5-substituents were studied at the B3LYP/6-311G* level of theory. It was found that there is no special stabilization or destabilization for any of the derivatives; the energetic effects that were previously attributed to aromatic stabilization or antiaromatic destabilization are the result of interactions in the reference systems. A nucleus-independent chemical shift (NICS) scan study at the HF-GIAO/6-311+G* theoretical level of these and similar derivatives suggest that they all show different magnitudes of diamagnetic ring current. None of the derivatives shows a paramagnetic ring current. Thus, cyclopentadienes are neither aromatic nor antiaromatic. It is also concluded that a diamagnetic ring current is perhaps necessary but certainly not a sufficient condition for aromaticity. The NICS scan procedure describes the type of ring current in the system, whereas a single isotropic NICS value (i.e., NICS(1)) may wrongly assign the type of ring current. It is shown that neither NICS(1) nor the NICS scan procedure can be used as a single aromaticity criterion.     

Structure and stability of closo-hexaboranes and their heteroanalogs

The molecular and electronic structures of closo-hexaboranes B₆H₆ ^2–, B₆H₇ ^–, and B₆H₈ and closo-heterohexaboranes XYB₄H₄ (X = Y = CH, N; X = BH, Y = CH^–, N^–, NH, O) were studed by the ab initio (MP2(full)/6-311+G**) and density functional (B3LYP/6-311+G**) methods. The bridging H atoms in closo-hexaboranes B₆H₇ ^– and B₆H₈ can undergo facile low-barrier migrations around the boron cage (the barrier heights are about 10—15 kcal mol^–1). All heteroboranes having octahedron-like structures with hypercoordinated N and O atoms are rather stable and can be the subject of synthetic research efforts.     

A new 3D-aromatic organoboron species on the basis of CB6 unit: Two states of carbon hypercoordination and structural isomerism of non-classical forms

Calculations using the DFT B3LYP/6-311+G(df,2p) and wb97XD/6-311+G(df,2p) methods predict the stability of a new three-dimensional aromatic organoboron species based on the CB₆ unit. Derived systems containing carbon atoms in two hypercoordination states (planar hexacoordination and inverted-umbrella heptacoordination) can be used as building blocks for constructing various stable condensed systems with many hypercoordinated carbon centers.     

硝酸甲与硝酸乙酯的结构、振动频率和热力学性质的密度函数 理论(DFT)研究

用密度函数理论(DFT)的BLYP和B3LYP方法,取6-31G,6-31G^*,6-31G^*^*,6-311G,6-311G^*和6-311G^*^*六种基组,对硝酸甲酯和硝酸乙酯的几何构型和红外振动频率进行了计算研究.结果表明,B3LYP方法在采用极化基组(6-31G^*,6-31G^*^*,6-311G^*和6-311G^*^*)时计算得到的结果均较好,适用于硝酸酯类化合物的研究.而BLYP方法无论采用何种基组均不适用;运用校正后的B3LYP/6-31G^*频率(校正因子0.975)计算得到的热力学性质(C^o~p,H^o和S^o)与实验结果较吻合。     

Polar substituent effects in the bicyclo[1.1.1]pentane ring system: acidities of 3-substituted bicyclo[1.1.1]pentane-1-carboxylic acids

Experimental gas-phase acidities are reported for a series of 3-substituted (X) bicyclo [1.1.1]pent-1-yl carboxylic acids (1, Y = COOH). A comparison with available calculated data (MP2/6-311++G**// B3LYP/6-311+G**) reveals good agreement. The relative substituent effects are shown to be adequately described by a much lower level of theory (B3LYP/6-31+G*). Various correlations are presented which clearly point to polar field effects as being the origin of the relative acidities.     

儿茶素结构和振动光谱的密度泛函理论研究

本文利用密度泛函理论方法B3LYP在6-311 G(d,p)水平上研究了儿茶素及其异构体表儿茶素分子的几何构型,计算结果与实验所得结构参数一致.计算了儿茶素分子平衡构型下的力常数,使用Wilson的GF矩阵方法计算了振动基频以及相应的势能分布,据此结合理论计算的光谱强度,对儿茶素分子的振动基频进行了完善合理的理论归属.     

Stabilization of octacoordinate carbon center in metal-containing derivatives of orthocarbonic acid

The electronic and spatial structures of alkali metal compounds CO₃M₂, CO₃M₃ ⁺, and CO₄M₄ (M = Li, Na, K) were investigated by the ab initio (MP2(full)/6-311+G**) and density functional (B3LYP/6-311+G**) methods. The calculated energies of formation decrease in the order E ^Li > E ^Na > E ^K for all structural types, being determined by steric and orbital interactions. Stable structures with octacoordinate carbon are formed in the case of CO₄M₄ salts. Dedicated to Academicians A. L. Buchachenko and N. S. Zefirov on the occasions of their 70th birthdays. Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 9, pp. 1929–1938, September, 2005.     

Ab initio and density functional theory studies of the structure,gas-phase acidity and aromaticity of tetraselenosquaric acid

Results of ab initio self-consistent-field (SCF) and density functional theory (DFT) calculations of the gas-phase structure,acidity (free energy of deprotonation,G0) and aro-maticity of tetraselenosquaric acid (3,4-diselenyl-3-cy-dobutene-1,2-diselenone,H2C4Se4) are reported.The global minimum found on the potential energy surface of tetraselenosquaric acid presents a planar conformation.The ZZ iso-mer was found to have the lowest energy among the three planar conformers and the ZZ and ZE isomers are very dose in energy.The optimized geometric parameters exhibit a bond length equalization relative to reference compounds,cyclobu-tanediselenone,and cydobutenediselenol.The computed aromatic stabilization energy (ASE) by homodesmotic reaction is -77.4 (MP2(fu)/6 - 311 G //RHF/6 - 311 G) and - 54.8 kJ/mol (B3LYP/6 - 311 G //B3LYP/6 -311 G).The aromaticity of tetraselenosquaric add is indicated by the calculated diamagnetic susceptibility exaltation (A) - 19.13 (CSGT(IGAEM) - RHF/6 - 311 G// RHF/6-     

Double carbon-halogen bond in the compounds H2CX+, H2CCX+, and H2BCX (X = F, Cl)

The structure of formaldehyde and ketene analogs H₂CX and H₂CCX (X = O, F⁺, Ne²⁺, S, Cl⁺, Ar²⁺), and also of boron-containing compounds H₂BCX (X = F, Cl), was studied by ab initio [CCD(full)/6-311+G**] and DFT (B3LYP/6-311+G**) calculations. In all the halogen-containing species except H₂BCF, a double carbon-halogen bond is formed.Translated from Zhurnal Obshchei Khimii, Vol. 74, No. 10, 2004, pp. 1649–1654.Original Russian Text Copyright © 2004 by Minyaev, Gribanova.This revised version was published online in April 2005 with a corrected cover date.     

Theoretical analysis of the internal rotation in aminoborane and borylphosphine

Using a recently proposed orbital deletion procedure and the block-localized wavefunction method, the rotational barriers in H₂BNH₂ and H₂BPH₂ are analyzed in terms of conjugation, hyperconjugation, steric effect and pyramidalization. With the zero-point energy corrections, the π-binding strengths in the planar H₂BNH₂ and H₂BPH₂ are both around 20 kcal/mol at the HF level using the 6-311+G** basis set. With the deactivation of the π atomic orbitals on the boron atom and the evolution from a planar structure to a 90°-twisted structure, the steric repulsion between the B‐H and the N‐H or P‐H is relieved and moreover, the negative hyperconjugation from the lone electron pair or pairs on the nitrogen or phosphorus atoms to the antibonding orbital ^χ* _B H ₂ of the BH₂ group stabilizes the twisted structure by 7.4(8.8) or 4.0(5.0) kcal/mol at the HF/6-31G*(6-311+G**) level. However, the repulsive interaction between the lone pair(s) and the two BH σ bonds is so prominent that the overall steric effect contributes 20.3(22.9) and 19.3(19.8) kcal/mol to the rotational barriers in H₂BNH₂ and H₂BPH₂ with the 6-31G*(6-311+G**) basis set. The present techniques and analyses may also give some clues to justify the parameterization in the empirical molecular mechanics methods. Received: 17 April 1998 / Accepted: 17 September 1998 / Published online: 1 February 1999     

Geometry,stability, and isomerization of BnN2 (n = 1−6) isomers

We perform a systematic study on the geometry, stability, nature of bonding, and potential energy surface of low‐lying isomers of planar and cyclic B_nN₂ (n = 1?6) at the CCSD(T)/6‐311+G(d)//B3LYP/6‐311+G(d) level. B_nN₂ (n = 2?4) clusters are structurally similar to pure boron clusters. The evolution of the binding energy per atom, incremental binding energy, and second‐order difference of total energy with the size of B_nN₂ reveals that the lowest energy isomer of B₃N₂ has high stability. B₅N₂ and B₆N₂ possess π‐aromaticity according to Hückel (4n + 2) rule. The aromaticity of some isomers of B₄N₂ and B₆N₂ is examined based on their valence molecular orbitals. At the CCSD(T)/6‐311+G(d)//B3LYP/6‐311+G(d) level, several B₂N₂, B₃N₂, B₄N₂, and B₅N₂ isomers are predicted to be stable both thermodynamically and kinetically, and detectable in future experiments. © 2013 Wiley Periodicals, Inc.     

Theoretical study of electroreductive intramolecular coupling of nonconjugated olefinic and aromatic ketones

The electroreductive cyclization of hept-6-en-2-one, octa-7-en-2-one, and 5-phenylpentan-2-one was investigated by ab initio (UHF/6-311++G**) and density functional (UB3LYP/6-311++G**) calculation methods. The high regio- and stereoselectivities previously reported for olefinic ketones (exo-trans) and 5-arylpentan-2-ones (endo-trans) were reconfirmed. These experimental results well agree with the computational outcome for the transition states in the intramolecular coupling of ketyl radicals generated by one-electron transfer to the ketones.     

HCN exchange on [Cu(HCN)4]+: a quantum chemical investigation

Density functional (B3LYP, B3PW91, X3LYP, BP86, PBEPBE, PW91PW91, and M06) and ab initio (MP2, MP4sdq, CCSD, and CCSD(T)) calculations with extended basis sets (6-311+G**, TZVP, LANL2DZ+p, and SDD+p, the latter including extra polarization and diffuse functions) indicate that HCN exchange on [Cu(HCN)₄]⁺ proceeds via an associative interchange (I_a) mechanism and a D_3h transition structure {[Cu(HCN)₅]⁺}^?. The activation barrier, relative to the model complex [Cu(HCN)₄]⁺·HCN, varies modestly, depending on the computational level. Typical values are 8.0?kcal?M^?1 (B3LYP/6-311+G**), 6.0?kcal?M^?1 (M06/6-311+G**), and 4.8?kcal?M^?1 (CCSD(T)/6-311+G**//MP2(full)/6-311+G**). Inclusion of an implicit solvent model (B3LYP(CPCM)/6-311+G**) leads to an activation barrier of 5.8?kcal?mol^?1. Comparison of the HCN exchange mechanisms on [Li(HCN)₄]⁺ (limiting associative, A) and [Cu(HCN)₄]⁺ (associative interchange, I_a) reveals that π back donation in the equatorial Cu–N bonds in the transition state determines the mechanism.