Ligand Exchange Processes on Solvated Beryllium Cations. IV [Be(H2O)2(imidazole‐based Chelate)]1)

The water exchange reaction of [Be(H₂O)₂(1H‐imidazole‐4,5‐dicarboxylate)] and [Be(H₂O)₂(1H‐imidazol‐3‐ium‐4,5‐dicarboxylate)]⁺ in water was studied by DFT calculations (RB3LYP/6‐311+G**) and identified as an associative interchange mechanism. The activation barriers for [Be(H₂O)₂(1H‐imidazole‐4,5‐dicarboxylate)] (16.6 kcal/mol) and [Be(H₂O)₂(1H‐imidazol‐3‐ium‐4,5‐dicarboxylate)]⁺ (13.8 kcal/mol) are similar to the barrier for [Be(H₂O)₄)]²⁺ and independent of the overall charge. NICS calculations show no indication that the aromaticity of the imidazole ring is affected during the water exchange process.     

Ligand Exchange Processes on Solvated Beryllium Cations. II [Be(solvent)(12‐Crown‐4)]2+

The solvation and solvent exchange mechanism of [Be(12‐crown‐4)]²⁺ in water and ammonia was studied by DFT calculations (RB3LYP/6‐311+G**). In solution, five‐fold coordinated Be²⁺ species of quadratic pyramidal [Be(H₂O)(12‐crown‐4)]²⁺ and [Be(NH₃)(12‐crown‐4)]²⁺ exist. The water and ammonia exchange reactions follow an associative interchange mechanism, similar to that found for the pure solvent complexes [Be(H₂O)₄]²⁺ and [Be(NH₃)₄]²⁺. The activation barriers are clearly smaller than for the pure solvent complexes, viz. [Be(H₂O)(12‐crown‐4)]²⁺: 6.0 kcal/mol and [Be(NH₃)(12‐crown‐4)]²⁺: 15.3 kcal/mol.     

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     

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     

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.     

Ligand‐Exchange Processes on Solvated Zinc Cations: Water Exchange on [Zn(H2O)4(L)]2+⋅2 H2O (L=Heterocyclic Ligand)

The water‐exchange mechanisms of [Zn(H₂O)₄(L)]²⁺?2 H₂O (L=imidazole, pyrazole, 1,2,4‐triazole, pyridine, 4‐cyanopyridine, 4‐aminopyridine, 2‐azaphosphole, 2‐azafuran, 2‐azathiophene, and 2‐azaselenophene) have been investigated by DFT calculations (RB3LYP/6‐311+G**). The results support limiting associative reaction pathways that involve the formation of six‐coordinate intermediates [Zn(H₂O)₅(L)]²⁺?H₂O. The basicity of the coordinated heterocyclic ligands shows a good correlation with the activation barriers, structural parameters, and stability of the transition and intermediate states.     

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.     

Ligand Exchange Processes on Solvated Beryllium Cations V. Water Exchange on [Be(X)(H2O)3]+ (X: H–,F–, Cl–, Br–, OH–, CN–, NCNCN–)‡

The water exchange reaction of Be^II complexes in the series [Be(X)(H₂O)₃]⁺ (X = H^–, F^–, Cl^–, Br^–, OH^–, CN^–, NCNCN^–) was studied by DFT calculations (RB3LYP/6‐311+G**) and identified as an associative interchange mechanism. The influence of X on the activation energies was examined and found to be largely negligible, thus making them all act as spectator ligands. The energies for addition of a fourth water molecule, representing the second coordination sphere, were approximately half of that found for similar dicationic complexes and close to that found for monocationic species like [Li(H₂O)₄]⁺.     

Internal conrotation and disrotation in H2BCH2BH2 and diborylmethane 1,3 H exchange

The paths of correlated internal disrotation (barrier less than 0.4 kcal/mol) and conrotation (barrier around 1.9 kcal/mol) of the two BH₂ groups in H₂BCH₂BH₂ have been computed employing ab initio [MP2(full)/6–31G**] and density functional theory (Becke3LYP/6–311+G**) methods. Two B(SINGLE BOND)C(DOTTED BOND)B(p) hyperconjugative interactions stabilize the C_s symmetric H₂BCH₂BH₂ isomer ( 1 ). The B(SINGLE BOND)C(DOTTED BOND)B(p) hyperconjugative stabilization, evaluated by homodesmotic reactions and using the orbital deletion procedure (which “deactivates” the “vacant” born p orbital), is less than 6 kcal/mol in diborylmethane. The B(SINGLE BOND)C(DOTTED BOND)B(p) stabilization is shown to be remarkably large in C₄B₆H₁₀ (Td). At MP2(fu)/6–31G**, disproportionation into 1 and methane is only 5.6 kcal/mol exothermic. The 1,3 H exchange in diborylmethane is an asynchronous process and proceeds via a doubly bridged cyclic intermediate with 9.3 kcal/mol barrier. Structures with “planar tetracoordinate” carbon are stabilized considerably by BH₂ substituents, but they are still high in energy. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 1792–1803, 1997     

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.     

Bonding and correlation analysis of various Si2CO isomers on the potential energy surface

At various levels of theory, singlet and triplet potential energy surfaces (PESs) of Si₂CO, which has been studied using matrix isolation infrared spectroscopy, are investigated in detail. A total of 30 isomers and 38 interconversion transition states are obtained at the B3LYP/6‐311G(d) level. At the higher CCSD(T)/6‐311+G(2d)//QCISD/6‐311G(2d)+ZPVE level, the global minimum ¹1 (0.0 kcal/mol) corresponds to a three‐membered ring singlet O‐cCSiSi (¹A′). On the singlet PES, the species ¹2 (0.2 kcal/mol) is a bent SiCSiO structure with a ¹A′ electronic state, followed by a three‐membered ring isomer Si‐cCSiO (¹A′) ¹3 (23.1 kcal/mol) and a linear SiCOSi ¹4 (¹Σ⁺) (38.6 kcal/mol). The isomers ¹1, ¹2, ¹3 , and ¹4 possess not only high thermodynamic stabilities, but also high kinetic stabilities. On the triplet PES, two isomers ³1 (³B₂) (18.8 kcal/mol) and ³7 (³A″) (23.3 kcal/mol) also have high thermodynamic and kinetic stabilities. The bonding natures of the relevant species are analyzed. The similarities and differences between C₃O, C₃S, SiC₂O, and SiC₂S are discussed. The present results are also expected to be useful for understanding the initial growing step of the CO‐doped Si vaporization processes. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Gas-phase electron diffraction and quantum chemical study of the structure of the 4-nitrobenzene sulfonyl chloride molecule

A combined gas-phase electron diffraction and quantum chemical (B3LYP/6-311+G**, B3LYP/cc-pVTZ, MP2/6-31G*, and MP2/cc-pVTZ) study of the structure of the 4-nitrobenzene sulfonyl chloride molecule is performed. It is found that at a temperature of 391(3) K only one conformer with C _s symmetry is present in the gas phase. The following experimental values of structural parameters are obtained: r _h1(C-H)_av = 1.086(6) Å, r _h1(C-C)_av = 1.395(3) Å, r _h1(C1-S) = 1.773(4) Å, r _h1(S=O) = 1.423(3) Å, r _h1(S-Cl) = 2.048(4) Å, r _h1(N-O) = 1.224(3) Å, r _h1(N-C4) = 1.477(3) Å, ∠(C1-S=O) = 109.0(4)°, ∠(Cl-S-O) = 106.7(2)°, ∠C1-S-Cl = 100.2(13)°, ∠O=S=O = 122.9(11)°, ∠O=N=O = 123.6(5)°. The C2-C1-S-Cl torsion angle that characterizes the position of the S-Cl bond relative to the benzene ring plane is 89(4)°. The NO₂ group lies in the benzene ring plane. Internal rotation barriers calculated by B3LYP/6-311+G** and MP2/6-31G* methods are: V ₁ = 4.7 kcal/mol and 5.3 kcal/mol for the sulfonyl chloride group; V ₂ = 4.9 kcal/mol and 6.0 kcal/mol for the nitro group.     

Ligand Exchange Processes on Solvated Zinc Cations – DFT Analysis of Hydrogen Cyanide Exchange on [Zn(HCN)6]2+

The mechanisms for hydrogen cyanide exchange on [Zn(HCN)₆]²⁺ were studied using density functional theory (B3LYP/6‐311+G**), and showed that the limiting dissociative (D) pathway is more favorable than the associative interchange (I_a) mechanism. The activation barrier for the dissociative mechanism (7 kcal · mol^–1) is clearly lower than for the interchange mechanism (15.9 kcal · mol^–1).     

A computational mechanistic study of the deamination reaction of melamine

A detailed computational study of the deamination reaction of melamine by OH^–, n H₂O/OH^–, n H₂O (where n = 1, 2, 3), and protonated melamine with H₂O, has been carried out using density functional theory and ab initio calculations. All structures were optimized at M06/6‐31G(d) level of theory, as well as with the B3LYP functional with each of the basis sets: 6‐31G(d), 6‐31 + G(d), 6‐31G(2df,p), and 6‐311++G(3df,3pd). B3LYP, M06, and ω B97XD calculations with 6‐31 + G(d,p) have also been performed. All structures were optimized at B3LYP/6‐31 + G(d,p) level of theory for deamination simulations in an aqueous medium, using both the polarizable continuum solvation model and the solvation model based on solute electron density. Composite method calculations have been conducted at G4MP2 and CBS‐QB3. Fifteen different mechanistic pathways were explored. Most pathways consisted of two key steps: formation of a tetrahedral intermediate and in the final step, an intermediate that dissociates to products via a 1,3‐proton shift. The lowest overall activation energy, 111 kJ mol^?1 at G4MP2, was obtained for the deamination of melamine with 3H₂O/OH^?.     

Water exchange on beryllium complexes: part VIII – influence of neutral electron pair donors

In this article, water exchange reactions on [Be(L)(H₂O)₃]²⁺ (L?=?NH_{3? x} (CH₃) _x , PH_{3? x} (CH₃) _x , AsH_{3? x} (CH₃) _x , OH_{2? x} (CH₃) _x , SH_{2? x} (CH₃) _x , SeH_{2? x} (CH₃) _x , pyridine, 4-fluoropyridine, 4-bromopyridine, 4-chloropyridine, 4-hydroxypyridine, 4-thiolopyridine, 4-selenidopyridine, 4-nitrilopyridine, 1,4-diazine, 1,3,5-triazine, HCN, acetonitrile, and benzonitrile) are examined, utilizing the B3LYP//6-311?+?G** density functional for geometry optimizations, and B3LYP//6-311?+?G** both with and without the CPCM solvent model as well as MP2(full)//6-311?+?G** for subsequent single-point energy calculations. In all examined cases, the results prove that these complexes show associative interchange mechanisms for water exchange. With the exception of the NH _x (CH₃)_{3? x} series of ligands, activation energy barriers vary little, making these ligands mostly spectator ligands. Geometrical parameters vary mainly with the ligand size.     

Quantum chemical study of the coordination of glycolic acid conformers and their conjugate bases to [Ca(OH2)n] (n=0–4) ions

A detailed exploration of the configurational and conformational space of glycolic acid and their conjugate bases has been carried out with the aid of first principles quantum chemical techniques at the B3LYP/6-311+G(d,p) and CCSD(T)/6-31G(d,p) levels of theory. The most stable configuration among the eight possible glycolic acid conformers corresponds to the E-s-cis, s-trans configuration, while the highest energy E-s-trans, s-cis conformer was found at 10.88 and 12.17 kcal mol⁻¹ higher in energy at the B3LYP/6-311+G(d,p) and CCSD(T)/6-31G(d,p) levels of theory, respectively. Upon dissociation of glycolic acid the s-cis(syn), and s-trans(anti) configurations of the glycolate anion can be formed. The anti conformer was found to be less stable than the syn one by 14.20 and 16.87 kcal mol⁻¹ at the B3LYP/6-311+G(d,p) and CCSD(T)/6-31G(d,p)) levels of theory, respectively. The computed B3LYP/6-311+G(d,p) proton affinity of the syn conformer for the protonation process affording the more stable E-s-cis, s-trans conformer, in vacuum was found to be 325.35 kcal mol⁻¹ (ΔG⁰ value). From a methodological point of view, our results confirm the reliability of the integrated computational tool formed by the B3LYP density functional model. This model has subsequently been used to investigate the interaction of Ca²⁺ ions with the glycolic acid conformers and their conjugate bases in vacuum and in the presence of extra water ligands. For the complexes of glycolic acid conformers the η²–O,O–(COOH) coordination, that is the structure that arises from the coordination of the Ca²⁺ to the carboxylic group, is the global minimum of the PES, while the η²–O(OH),O–(COOH) coordination is a local minimum found at only 1.0 and 1.3 kcal mol⁻¹ higher in energy at the B3LYP/6-311+G(d,p) and CCSD(T)/6-31G(d,p) levels of theory, respectively. Moreover, the two isomers exhibit nearly the same binding affinities, which are predicted to be 89 and 85 kcal mol⁻¹ at the B3LYP/6-311+G(d,p) and CCSD(T)/6-31G(d,p) levels of theory, respectively. The same holds also true for the complexes of the glycolate anion. The η²–O,O–(COO⁻) coordination involving the syn conformer of the glycolato ligand, is the global minimum, while the η²–O(OH),O–(COO⁻) one lies at 1.5 and 5.6 kcal mol⁻¹ higher in energy at the B3LYP/6-311+G(d,p) and CCSD(T)/6-31G(d,p) levels of theory, respectively. The other conformer with an η²–O,O–(COO⁻) coordination involving the anti conformer of the glycolato ligand, is less stable by only 0.2 kcal mol⁻¹ at both levels of theory. Noteworthy is the trend seen for the incremental binding energy due to the successive addition of water molecules to [HOCH₂C(O)O]Ca²⁺ species; the computed values are 30.4, 26.8, 22.9 and 16.2 kcal mol⁻¹ at the B3LYP/6-311+G(d,p) level of theory for the mono-, di-, tri- and tetraaqua complexes, respectively. This trend arising from the repulsion of the dipoles between the water ligands and from unfavorable many body interactions is in accordance with those anticipated from electrostatic considerations. The Ca(II)-water interaction weakens with increasing coordination of the metal. Obviously, it is the electrostatic nature of the Ca(II)-water interactions that accounts well for the computed coordination geometries of the cationic (aqua)(glycolato)calcium complexes. Calculated structures, relative stability and bonding properties of the conformers and their complexes with [Ca(OH₂)_n]²⁺ (n=0–4) ions are discussed with respect to computed electronic and spectroscopic properties, such as charge density distribution, harmonic vibrational frequencies and NMR chemical shifts.     

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     

[Be3(μ3‐O)3(MeCN)6{Be(MeCN)3}3](I)6 – ein Berylliumkomplex mit Cyclo‐(Be3O3)‐Kern und sein Hydrolyseprodukt [Be(H2O)4](I)2·2MeCN

Single crystals of [Be₃(μ₃‐O)₃(MeCN)₆{Be(MeCN)₃}₃](I)₆·4CH₃CN ( 1 ·4CH₃CN) were obtained in low yield by the reaction of beryllium powder with iodine in acetonitrile suspension, which probably result from traces of beryllium oxide containing the applied beryllium metal. The compound 1 ·4CH₃CN forms moisture sensitive, colourless crystal needles, which were characterized by IR spectroscopy and X‐ray diffraction (Space group Pnma, Z = 4, lattice dimensions at 100(2) K: a = 2317.4(1), b = 2491.4(1), c = 1190.6(1) pm, R₁ = 0.0315). The hexaiodide complex cation 1 ⁶⁺consists of a cyclo‐Be₃O₃ core with slightly distorted chair conformation, stabilized by coordination of two acetonitrile ligands at each of the beryllium atoms and by a {Be(CH₃CN)₃}²⁺ cation at each of the oxygen atoms. This unique coordination behaviour results in coplanar OBe₃ units with short Be–O distances of 155.0 pm and 153.6 pm on average of bond lengths within the cyclo‐Be₃O₃ unit and of the peripheric BeO bonds, respectively. Exposure of compound 1 ·4CH₃CN to moist air leads to small orange crystal plates of [Be(H₂O)₄]I₂·2CH₃CN ( 3 ·2CH₃CN). According to the crystal structure determination (Space group C2/c, Z = 4, lattice dimensions at 100(2) K: a = 1220.7(1), b = 735.0(1), c = 1608.5(1) pm, β = 97.97(1)°, R₁ = 0.0394), all hydrogen atoms of the dication [Be(H₂O)₄]²⁺ are involved to form O–H ··· N and O–H ··· I hydrogen bonds with the acetonitrile molecules and the iodide ions, respectively. Quantum chemical calculations (B3LYP/6‐311+G**) at the model [Be₃(μ₃‐O)₃(HCN)₆{Be(HCN)₃}₃]⁶⁺ show that chair and boat conformation are stable and that the distorted chair conformation is stabilized by packing effects.     

Structural and energetic exploration of a boron‐rich sulfide cluster B6S

Boron and mixed‐boron clusters have received considerable attention because of their wide applications and their essential roles in advancing chemical bonding models. Bearing the bright prospects as building blocks to form novel polymeric materials, the sulfur‐rich boron sulfides have been greatly studied. However, the knowledge of the boron‐rich boron sulfides is much rare. In this article, we report an extensive theoretical study on the structural, energetic, and stability features of a hitherto unknown septa‐atomic cluster B₆S at the CCSD(T)/6‐311+G(2df)//B3LYP/6‐311+G(d) level. The local minimum isomers were obtained through our recently developed program “grid‐based comprehensive isomeric search algorithm.” The results show that the planar knife‐like isomer B₅(?BS) 01 (0.0 kcal/mol) containing the ?BS moiety is the lowest energy, followed by the quasi‐planar belt‐like isomer B₆(>S) 02 (6.7 kcal/mol) and the pyramid‐like isomer B₆(>S) 03 (8.4 kcal/mol). Notably, the three singlet isomers all have good kinetic stability on the basis of the potential energy surface analysis. The B/S‐centered wheel‐like isomers are unfavorable in thermodynamics and kinetics. The triplet state structures generally can not compete with the singlet ones. The results are compared to the analogous and isoelectronic cluster B₆O. Our work is expected to provide useful information for understanding the structures and stability of boron sulfides. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012