Structures and stability of N7+ and N7− clusters

Ab initio molecular orbital theory and density functional theory have been used to study nine isomers of N₇ ionic clusters with low spin at the HF/6-31G*, MP2/6-31G*, B3LYP/6-31G*, and B3LYP/6-311(+)G* levels of theory. All stationary points are examined with harmonic vibrational frequency analyses. Four N₇ ⁺ isomers and five N₇ ⁻ isomers are determined to be local minima or very close to the minima on their potential-energy hypersurfaces, respectively. For N₇ ⁺ and N₇ ⁻, the energetically low lying isomers are open-chain structures (C _{2 v} and C _{2 v} or C₂). The results are very similar to those of other known odd-number nitrogen ions, such as N₅ ⁺, N₉ ⁺, and N₉ ⁻, for which the open-chain structures are also the global minima. This research suggests that the N₇ ionic clusters are likely to be stable and to be potential high-energy-density materials if they could be synthesized. Received: 16 July 2001 / Accepted: 8 October 2001 / Published online: 21 January 2002     

Calculations of thermodynamic stability of CpCoC<Subscript>2</Subscript>B<Subscript>9</Subscript>H<Subscript>11</Subscript> cobaltacarborane isomers

Quantum-chemical calculations of the geometry and energies of nine possible isomers of 12-vertex cobaltacarborane CpCoC₂B₉H₁₁ (1) were carried out by the DFT method (PBEPBE/DGDZVP/DGA1). Thermodynamic stability of the isomers increases with increasing distance between the carbon atoms in the cage and is virtually independent of the position of the CpCo vertex. The relative stabilities of the 1,2,3-(17.57 kcal mol⁻¹), 1,2,4-(3.72 kcal mol⁻¹), and 1,2,9-isomers of 1 (0 kcal mol⁻¹) are similar to the corresponding values for the ortho (17.61 kcal mol⁻¹), meta (3.21 kcal mol⁻¹), and para isomers (0 kcal mol⁻¹) of carborane C₂B₁₀H₁₂. The results of the present study confirm a close similarity of the CpCo and BH fragments in metallacarborane chemistry. Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 7, pp. 1557–1559, July, 2005.     

Structures and stability of N9, N9 − and N9 + clusters

 Ab initio molecular orbital calculations for N₉, N⁻ ₉ and N⁺ ₉ isomers were carried out at the HF/ 6-31G*, B3PW91/6-31G*, B3LYP/6-31G* and MP2/ 6-31G* levels of theory. Stable equilibrium geometric structures were determined by harmonic vibrational frequency analyses at the HF/6-31G*, B3PW91/6-31G* and B3LYP/6-31G* levels of theory. The most stable free-radical N₉ cluster is structure 1 with C _{2 v} symmetry and that of anion N⁻ ₉ is structure 3 with C _s symmetry. Only one stable structure of the N⁺ ₉ cation with C _{2 v} symmetry was predicted. Their potential application as high-energy-density materials has been examined. Received: 15 June 1999 / Accepted: 11 October 1999 / Published online: 14 March 2000     

Anisotropy in the physical properties of UNiGa under high pressure

The linear thermal expansion, compressibility and magnetostriction of UNiGa have been measured under high pressure. Huge anisotropic behaviors are observed in these physical quantities. The linear thermal expansion coefficients are α _a ∼ 16·10⁻⁶ K⁻¹ along thea-axis anda _c ∼5·10⁻⁶ K⁻¹ along thec-axis, and the linear compressibilities at room temperature are κ _a ∼ 3.6·10⁻³ GPa⁻¹ and κ _c ∼ 1.7·10⁻³ GPa⁻¹ alonga-axis andc-axes, respectively. UNiGa orders antiferromagnetically belowT _N=39 K and shows a metamagnetic transition at 4.2 K in magnetic fieldB _C=1 T. It is found thatT _N decreases andB _C increases with increasing pressure.     

Heteroorganic betaines

The structures of 6,6-dimethyl-6-silafulvene C₅H₄SiMe₂ (3), its donor-acceptor complex with ammonia. C₅H₄SiMe₂·NH₃, dimethylfulvene, a number of cyclopentadienylides, methylenetrimethylphosphorane (6), and silicon-containing organophosphorus betaine⁻C₅H₄SiMe₂CH₂PMe₃ ⁺ (13), the product of nucleophilic addition of6 to3, were calculated using the density functional approach. For compound13, the potential energy minimum corresponds to the conformation withgauche-arrangement of the cyclopentadienyl anionie and trimethylphosphonium cationic centers and a C−Si−C−P dihedral angle of 30.5°, which is due to the Coulomb attraction between these centers. According to calculations, betaine13 is rather stable toward decomposition into3 and6 (ΔH ^o=42 kcal mol⁻¹, ΔG ^Δ=30 kcal mol⁻¹). The main channel of thermal decomposition of compound13 involves an intramolecular nucleophilic substitution, which proceeds with elimination of trimethylphosphine and results in 1,1-dimethyl-1-silaspiro[2,4]hepta-4,6-diene, which then undergoes a ready and irreversible isomerization into 6,6-dimethyl-6-silabicyclo[3.2.0]hepta-1,3-diene owing to the [1.5]-sigmatropic shift of the C−Si bond. For Part 4, see Ref. 1. Published inIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 11, pp. 1850–1857, November, 2000.     

A theoretical investigation on structures of tripodal thiourea derivatives and their anion recognition

The structural geometries of three tripodal thiourea receptors, i.e. 1,3,5-triethyl-2,4,6-tris[(N′-methylthioureido)methyl]benzene (1), tris[N′-methyl-N-(2-aminoethyl)thiourea]methane (2), tris[N′-methyl-N-(2-aminoethyl)thiourea]amine (3), and their complexes with F⁻, Cl⁻, Br⁻, I⁻, NO₃ ⁻, CO₃ ²⁻, SO₄ ²⁻, HSO₄ ⁻, PO₄ ³⁻, HPO₄ ²⁻ and H₂PO₄ ⁻ were obtained using the density functional theory calculations. Electronic and thermodynamic properties of anion binding complexes of the receptors 1, 2 and 3 were investigated. Recognition abilities of all the receptors in terms of selectivity coefficients are reported. Intermolecular interactions in all the studied complexes occurring via multi-point hydrogen bonding were found. The receptors 1, 2 and 3 were found to be excellent selectivity for phosphate ion and their binding free energy for the phosphate ion are −292.57, −291.77 and −295.01 kcal/mol, respectively.     

Isomerism of metal complexes with the boron cluster anions B10H 102? and B12H 122?

A new type of inner-sphere isomerism in coordination chemistry observed in metal complexes with boron cluster anions B _n H _n ²⁻ (n = 10, 12) is discussed. Specific structure of the closo-borohydride anions gives rise to edge and facial isomers, among which mirror isomers have been found. Examples of edge and facial isomers of metal complexes with the B₁₀H₁₀²⁻ and B₁₂H₁₂²⁻, anions obtained experimentally and studied by X-ray diffraction are considered. Analysis of their structure revealed a new type of isomerism, which was called “positional.”     

Liquid-crystalline CuII and PdII complexes with nonmesogenic ferrocene-containing β-aminovinyl ketone

New liquid-crystalline heteropolynuclear complexes L₂M (M=Cu²⁺ (2a), Pd²⁺ (2b)) were synthesized by the reactions of C₅H₅FeC₅H₄−C₆H₄NH−C₂H₂−(CO)−C₆H₄OC₁₂H₂₅ (1, LH) with copper(ii) and palladium(ii) acetates. Compound2b was found to possess monotropic nematic and smectic phases;2a exhibits the monotropic nematic phase and a phenomenon of “double melting”. The compositions and structures of compounds1 and2a,b were established by elemental analysis,¹H and¹³C NMR, ESR, and IR spectroscopy. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 2, pp. 381–383, February, 1999.     

Carbometallic derivatives of hydroborons C2h and D5d: A list of substitution isomers

The problem of determining the number and type of X-substituted (where X represents certain substituents) carbometallic derivatives of hydroborons [Co(C₂B₉H₁₁)₂]⁻ C _2h and [Co(C₅B₆H₁₁)₂]⁻ D _5d (from apexes) was solved on the basis of G. Pólya’s theorem by means of combinatory analysis. The formulae of symmetry Z and generating functions of the number of achiral substitution isomers were determined. The family distributions of isomers (depending on the type and number of substituents) and the sites of possible substitution depending on the number m were determined. Mono and di-X-substituted isomers of [Co(C₅B₆H₁₁)₂]⁻ C _2h , as well as mono and di-X-substituted isomers of [Co(C₅B₆H₁₁)₂]⁻ D _5d , were identified. The procedure for plotting an additive model of calculating the properties of isomers of apical substitution of [Co(C₅B₆H₁₁)₂]⁻ D _5d was described on the basis of splitting of polygonal numbers (K ₃, K _TE, K ₅, and so on) of Pascal triangle, upon use of which there is no randomization in the choice of parameters of calculation model. The additive schemes containing 7 and 25 parameters for the calculation of properties of X-substituted carbometallic derivative of hydroboron [Co(C₅B₆H₁₁)₂]⁻ D _5d to various approximations were obtained.     

Synthesis, structure, and properties of a binuclear Fe(III) complex with N-(1-propanol)-N,N-bis(3-tert-butyl-5-methyl-2-hydroxybenxyl)amine

Mannich reaction of 2-Amino propanol, 2-tert-butyl-4-methylphenol, and formaldehyde in the ratio of 1:2:2 provides a new compound, N-(1-propanol)-N,N-bis(3-tert-butyl-5-methyl-2-hydroxybenxyl)amine (H₃L), which has been characterized by X-ray crystallography and elemental analysis. In the presence of Et₃N, the reaction of H₃L and FeCl₃·6H₂O gives a dinuclear Fe(III) complex [Fe₂L₂] 1, which has been characterized by X-ray crystallography, magnetic measurement, and cyclic voltammetry. The value of μ_eff at room temperature (5.97 μ_B) is much less than the expected spin-only value (8.37 μ_B) of two high spin (hs) Fe³⁺ (S = 5/2) ions [μ = g[∑ZS(S + 1)]^1/2], indicating there are strong coupling interactions between Fe³⁺ ions. The magnetic behavior of 1 denotes the occurrence of intramolecular antiferromagnetic interactions (J = −13.35 cm⁻¹ ). CV of 1 reveals two reversible waves at 0.433 and 1.227 V versus AgCl/Ag, which can be ascribed to the successive redox coupling of Fe^IIFe^II/Fe^IIIFe^II and Fe^IIIFe^II/Fe^IIIFe^III, respectively.     

Simulation of molecular and electronic structure of polyhydrogenated (n,0)-tubulenes and their analogs intercalated with lithium

Molecular and electronic structure of four polyhydrogenated (n,0)-tubulenes, namely, [−C₂₄H₄−] _m (1), two isomers of composition [−C₂₈H₄−] _m (2 and3), and [−C₃₂H₄−] _m (4) withn benzene rings in the cross section (n=6, 7, 7, and 8, respectively), was simulated atm>1 (m is the number of repeating fragemnts). It was assumed that hydrogen atoms are attached to all carbon atoms lying on the two most distant elements of the cylinders of the corresponding tubulenes. The energy band structures of macromolecules1–4 and their Li-intercalated analogs [−C₂₄H₄Li−] _m (5) [−C₂₈H₄Li−] _m (two isomers,6 and7), and [−C₃₂H₄Li−] _m (8), containing one Li atom per repeating unit at each center, were obtained in the EHT approximation by the crystal orbital method. Geometric parameters of repeating units of structures1–8 were found after MNDO/PM3 optimization of the energies of hydrocarbon molecules C₇₂H₂₄, C₈₄H₂₆ (two geometric isomers), and C₉₆H₂₈, containing three repeating units of corresponding tubulenes1–4 each. The conductivity types of polyhydrogenated tubulenes1–4 are the same as those of their precursors, (6,0)-, (7,0)-, and (8,0)-tubulenes. Dispersion curves of systems5–8 are much the same as those of macromolecules1–4; however, electron energy spectra of5–8 possess metallic conductivity type and the positions of Fermi levels for these systems are higher than for compounds1–4. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 11, pp. 2061–2067, November, 1999.     

Kinetics and mechanism of the base hydrolysis of cis-eq-[Cr(NCS)(S-pdtra)]−, trans-eq-[Cr(NCS)(edtrp)]− and trans-eq-[Cr(NCS)(R-pdtrp)]− complexes. Strong rate enhancement for NCS− ligand release with an increase of [OH−] only for the trans-equatorial isomers

The [Cr(NCS)(edtrp)]⁻, [Cr(NCS)(R-pdtrp)]⁻ and [Cr(NCS)(S-pdtra)]⁻ complexes, that are derivatives of the trans-equatorial isomers of [Cr(edtrp)(H₂O)]° and [Cr(R-pdtrp)(H₂O)]° and the cis-equatorial isomer of [Cr(S-pdtra)-(H₂O)]° (edtrp = ethylenediamine-N,N,N′-tripropionate, R-pdtrp = R-propane-1,2-diamine-N,N,N′-tripropionate, S-pdtra = S-propane-1,2-diamine-N,N,N′-triacetate) undergo aquation in alkaline media with a strong dependence of the rate on [OH⁻] for the trans-equatorial isomers and a very weak dependence for the cis-equatorial isomer. The thiocyanate ligand release follows a stereoretentive course for all reactants. Based on kinetic data the reaction mechanism has been discussed. Rate differences between the isomers are interpreted in terms of an interchange via a conjugate base (I c.b.) mechanism, assuming an equilibrium between the cis-equatorial-Cr^III-S-pdtra complexes with penta- and tetradentate coordination of the edta-like ligand. This revised version was published online in June 2006 with corrections to the Cover Date.     

Isomerism in OBe3F3 + cation: anab initio study

Ab initio MP2/6-31G^*//HF/6-31G^*+ZPE(HF/6-31G^*) calculations of the potential energy surface in the vicinity of stationary points and the pathways of intramolecular rearrangements between low-lying structures of the OBe₃F₃ ⁺ cation detected in the mass spectra of μ₄-Be₄O(CF₃COO)₆ were carried out. Ten stable isomers with di- and tricoordinate oxygen atoms were localized. The relative energies of six structures lie in the range 0–8 kcal mol⁻¹ and those of the remaining four structures lie in the range 20–40 kcal mol⁻¹. Two most favorable isomers, aC _2v isomer with a dicoordinate oxygen atom, planar six-membered cycle, and one terminal fluorine atom and a pyramidalC _3v isomer with a tricoordinate oxygen atom and three bridging fluorine atoms, are almost degenerate in energy. The barriers to rearrangements with the breaking of one fluorine bridge are no higher than 4 kcal mol⁻¹, except for the pyramidalC _3v isomer (∼16 kcal mol⁻¹). On the contrary, rearrangements with the breaking of the O−Be bond occur with overcoming of a high energy barrier (∼24 kcal mol⁻¹). A planarD _3h isomer with a tricoordinate oxygen atom and linear O−Be−H fragments was found to be the most favorable for the OBe₃H₃ ⁺ cation, a hydride analog of the OBe₃F₃ ⁺ ion; the energies of the remaining five isomers are more than 25 kcal mol⁻¹ higher. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 3, pp. 420–430, March, 1999.     

Boat form contributions in crowded piperidines: Theoretical and experimental study

The preferred conformations of N-nitroso-t(3)-alkyl-r(2),c(6)-bis(2′-furyl)-piperidin-4-ones 1–3 [alkyl = CH₃, C₂H₅ and CH(CH₃)₂] and N-nitroso-t(3),t(5)-dimethyl-r(2),c(6)-bis(2′-furyl)piperidin-4-one 4 in solutions were assigned by means of ¹H and ¹³C NMR studies. The results derived from NMR spectra indicate the presence of an equilibrium mixture of boat conformation B ₁ and alternate chair conformation CA for the E isomers of 1–3 and Z isomers of 2–3. For the Z isomer of 1 boat form B ₂ is predicted to be the major conformer. The N-nitroso-3,5-dimethyl derivative 4 exists in the boat form B ₁ only. Conformational analysis performed through semiempirical molecular orbital calculations also supports the conformations for 3–4. The presence of one conformer in the equilibrium can be predicted to a reasonable accuracy by theoretical studies in 1–2. The effects due to N-nitrosation on ¹H and ¹³C chemical shifts are also interpreted in terms of these conformations. The conformation of isopropyl group at C(3) was also predicted by spectral and theoretical studies.     

Theoretical study of the structure and stability of cage-substituted octahedral boranes,alanes, and gallanes

The equilibrium geometric parameters and energetic and spectroscopic characteristics of low-lying isomers for a series of successively cage-substituted octahedral closo-borane, alane, and gallane dianions Al_{6 − i} B _i H₆²⁻ and Ga_{6 − i} Al _i H₆²⁻ (i = 0−6) and mono- and disubstituted anions and neutral molecules of boranes, XB₅H₆⁻ and XYB₄H₆ and alanes, XAl₅H₆⁻ and XYAl₄H₆, with the same or different cage heteroatoms of Group IV elements (X and Y = C, Si, Ge, Ti) were calculated at the B3LYP level of the density functional theory using the 6-31G* and 6-311+G** basis sets. Differences in structure and stability between borane and alane clusters of like composition are revealed. The mutual influence of the X and Y heteroatoms in the trans and cis isomers and its manifestations in the behavior of molecular characteristics are considered for the disubstituted octahedral clusters XYB₄H₆ and XYAl₄H₆.     

Supramolecular self-assembly of two Cu(II) complexes with 1,2,4-triazole derivatives: syntheses,crystal structures and magnetic properties

Two two-dimensional coordination complexes, {[Cu₄(BTM)₆(OPA)₄] · 4DMF · 3H₂O} _n (1) and {[Cu(BDTM)(OH)](ClO₄) · 2H₂O} _n (2) (BTM = bis(1,2,4-triazol-1-yl)methane, BDTM = bis(3,5-dimethyl-1,2,4-triazol-1-yl)methane, OPA²⁻ = ortho-phthalic dianion, DMF = N,N-dimethylformamide), were synthesized and structurally characterized. Each Cu(II) ion locates in a distorted square pyramidal geometry in 1, in which OPA²⁻ ligands bridge Cu²⁺ ions along a axis to form a magnetic transmission chain and BTM ligands act as flexible spacers to construct the two-dimensional layer structure. In 2, each Cu²⁺ ion adopts tetra-coordination geometry to two hydroxyl groups and two triazolyl nitrogen atoms from two different BDTM ligands. Two hydroxyl groups bridge two Cu²⁺ ions to form a rhombic diamond, and four BDTM ligands connect four diamonds to form a 36-membered macrocyclic structure with large channels along a axis. Magnetic properties revealed that both OPA²⁻ and OH⁻ mediate anti-ferromagnetic interactions between Cu²⁺ ions with J = − 0.06(3) and −301.9(2) cm⁻¹ for 1 and 2, respectively. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Matrix IR-spectroscopic and quantum-chemical study of difluorostannylene complexation with dinitrogen

Complexes of difluorostannylene with dinitrogen of composition 1∶1 and 1∶2 were stabilized in Ar matrix (12 K) and characterized by IR spectra. The bands at 588, 565, and 583, 557 cm⁻¹, respectively, were assigned to these complexes. Potential energy surfaces of the systems SnF₂+N₂ and SnF₂+2N₂ were studied by theab initio MP2/3-21G(d2)//HF/3-21G(d2) method using the basis set including polarization functions at Sn, F, and N atoms. Equilibrium structures of the complexes haveC _s andC _2v symmetry and correspond to coordination of lone electron pairs of nitrogen molecules with vacant p-AO of the carbenic center. The calculated complexation energies are equal to 4.6 and 8.9 kcal mol⁻¹, respectively. Based on results of quantum-chemical calculations an interpretation of the IR spectra of the complexes was given and it was shown that cycloaddition of SnF₂ to a triple N≡N bond with formation ofcyclo-SnF₂N₂ is energetically unfavorable. The absorption band belonging to SiF₄·N₂ complex in Ar matrix was detected and assigned. Dedicated to the memory of Academician M. E. Vol'pin timed to his 75th birthday. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 6, pp. 1087–1093, June, 1998.     

Double H-bridged and single H-bridged diboryl radicals

The structures of B₂H₅·, B₂H₅CO·, and B₂H₅N₂· radicals are investigated using the 6–31G^* basis set. Both double H-bridged and single H-bridged isomers are found to be local minima on the potential energy surface. The effects of electron correlation are taken into account using single point MP4/6–31G^* calculations and, for the diboryl radicals, complete MP3/6–31G^* optimizations. In all cases the single H-bridged isomers are found to be more stable than the corresponding double H-bridged isomers.The transition state for the double H-bridged to single H-bridged B₂H₅· isomerization reaction is calculated to be 2.54 kcal mol^–1 above the double H-bridged radical at the MP4SDTQ/6-31G^*//UHF/ 6–31 G^* level when corrected for zero point energy. Barrier tunneling increased the reaction rate by a factor of 2.5–3.0, strongly suggesting the system is fluxional at this temperature.The addition of CO and N₂ to the diboryl radicals leads to relocation of the unpaired electron and rehybridization of the C and N atoms adjacent to the boron atoms. The isomers of B₂H₅CO· and B₂H₅N₂· are different and should be distinguishable experimentally. While the CO moiety is bound to the diboryl radicals isomers by over 19 kcal mol^–1, no binding energy is evident for N₂.     

Theoretical refinements of theendo and theexo structures of tetraborane(8) carbonyl

The molecular structures of theendo (1a) andexo (1b) isomers of B₄H₈CO have been optimized at the ab initio MP2(Full)/6-31G^* level of theory. The agreement of the computed geometrical parameters with the recently published electron-diffraction (GED) data is very good, even though a number of geometrical constraints were applied in the experimental determination. The IGLO (individual gauge for localized orbitals)¹¹B NMR chemical shifts, calculated at the II//MP2/6-31G^* level, are also in accord with experiment. The formation of1a and1b by association of B₄H₈ and CO is computed to be exothermic by 22.8 and 22.2 kcal/mol, respectively, at the MP2(Full)/6-31G^*//MP2(Full)/6-31G^* + ZPE(6-31G^*) level of theory. The Lewis acid strength of B₄H₈ toward CO is comparable to that of BH₃.     

Effects of substituents and n-donors on the energies of Si←N,Si←O,and Si=C bonds in hypervalent silenes

The effect of the nature of substituents at sp²-hybridized silicon atom in the R₂Si=CH₂ (R = SiH₃, H, Me, OH, Cl, F) molecules on the structure and energy characteristics of complexes of these molecules with ammonia, trimethylamine, and tetrahydrofuran was studied by the ab initio (MP4/6-311G(d)//MP2/6-31G(d)+ZPE) method. As the electronegativity, χ, of the substituent R increases, the coordination bond energies, D(Si← N(O)), increase from 4.7 to 25.9 kcal mol⁻¹ for the complexes of R₂Si=CH₂ with NH₃, from 10.6 to 37.1 kcal mol⁻¹ for the complexes with Me₃N, and from 5.0 to 22.2 kcal mol⁻¹ for the complexes with THF. The n-donor ability changes as follows: THF ≤ NH₃ < Me₃N. The calculated barrier to hindered internal rotation about the silicon—carbon double bond was used as a measure of the Si=C π-bond energy. As χ increases, the rotational barriers decrease from 18.9 to 5.2 kcal mol⁻¹ for the complexes with NH₃ and from 16.9 to 5.7 kcal mol⁻¹ for the complexes with Me₃N. The lowering of rotational barriers occurs in parallel to the decrease in D _π(Si=C) we have established earlier for free silenes. On the average, the D _π(Si=C) energy decreases by ∼25 kcal mol⁻¹ for NH₃· R₂Si=CH₂ and Me₃N·R₂Si=CH₂. The D(Si←N) values for the R₂Si=CH₂· 2Me₃N complexes are 11.4 (R = H) and 24.3 kcal mol⁻¹ (R = F). sp²-Hybridized silicon atom can form transannular coordination bonds in 1,1-bis[N-(dimethylamino)acetimidato]silene (6). The open form (I) of molecule 6 is 35.1 and 43.5 kcal mol⁻¹ less stable than the cyclic (II, one transannular Si←N bond) and bicyclic (III, two transannular Si←N bonds) forms of this molecule, respectively. The D(Si←N) energy for structure III was estimated at 21.8 kcal mol⁻¹. Dedicated to Academician N. S. Zefirov on the occasion of his 70th birthday. Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 9, pp. 1952–1961, September, 2005.