Exchanged Metal-Hydrogen Anagostic Bonds and Resonance of Dithiocarbamate and Thioureide Mesomers**

The pressure-induced transformation of plane-square complex nickel(II) bis(N,N-diethyldithiocarbamate) between its soft dithiocarbamate (form I) and thioureide (form II) mesomeres is coupled to the interchange of anagostic Ni⋅⋅⋅H−C interactions from methylene to the methyl group, respectively. At 1.23 GPa, the clearly visible giant anomalous compressibility of the crystal reveals a potential-energy difference of 5.4 kJ mol⁻¹ between the two complex forms. The structural and spectroscopic results, which are supported by quantum-mechanical calculations, connect this solid-state phase transition with the mesomeric transition, and this is accompanied by the conformational transformation of anagostic Ni⋅⋅⋅H−C rearrangement and formation of the charge-assisted S⁻⋅⋅⋅H−C bond under pressure.     

Implementation of molecular symmetry in valence bond calculation

A novel strategy for the construction of many-electron symmetry-adapted wave function is proposed for ab initio valence bond (VB) calculations and is implemented for valence bond self-consistent filed (VBSCF) and breathing orbital valence bond (BOVB) methods with various orbital optimization algorithms. Symmetry-adapted VB functions are constructed by the projection operator of symmetry group. The many-electron symmetry-adapted wave function is expressed in terms of symmetry-adapted VB functions, and thus the VB calculations can be performed with the molecular symmetry restriction. Test results show that molecular symmetry reduces the computational cost of both the iteration numbers and CPU time. Furthermore, excited states with specific symmetry can be conveniently obtained in VB calculations by using symmetry-adapted VB functions.     

Microwave spectra of the Xe-N2 van der Waals complex: a comparison of experiment and theory

Rotational transitions for the Xe-N2 complex were measured in the frequency region from 4 to 18 GHz using a pulsed-nozzle Fourier-transform microwave spectrometer. Twelve (four) a-type transitions were recorded for the 132Xe-14N2 and 129Xe-14N2 (131Xe-15N)) isotopomers. In addition, the nuclear quadrupole hyperfine structures due to the presence of the 14N (nuclear-spin quantum number I=1) and 131Xe (I=32) nuclei were detected and analyzed. Two ab initio potential-energy surfaces were calculated at the coupled-cluster level of theory with single, double, and pertubatively included triple excitations. Dunning's augmented correlation-consistent polarized valence triple-zeta basis set was used for the nitrogen atoms. For the first surface, a well-tempered basis set with additional polarization functions was used for the Xe atom; for the second surface, a newly developed augmented correlation-consistent polarized valence quintuple-zeta basis set employing small-core relativistic pseudopotentials was used for the Xe atom. The basis sets were supplemented with bond functions for the van der Waals bond. The counterpoise correction was applied to reduce the basis-set superposition error. The resulting two surfaces both have a single minimum at a T-shaped geometry, with well depths of 122.4 and 119.3 cm(-1), respectively. Bound-state energies supported by the potential-energy surface were determined. The quality of the ab initio potential-energy surfaces was evaluated by comparison of the experimental transition frequencies and rotational and centrifugal distortion constants with those derived from the bound-state energies. A scaled potential-energy surface was obtained which has excellent agreement with the experimental data.     

Combined valence bond-molecular mechanics potential-energy surface and direct dynamics study of rate constants and kinetic isotope effects for the H + C2H6 reaction

This article presents a multifaceted study of the reaction H+C(2)H(6)-->H(2)+C(2)H(5) and three of its deuterium-substituted isotopologs. First we present high-level electronic structure calculations by the W1, G3SX, MCG3-MPWB, CBS-APNO, and MC-QCISD/3 methods that lead to a best estimate of the barrier height of 11.8+/-0.5 kcal/mol. Then we obtain a specific reaction parameter for the MPW density functional in order that it reproduces the best estimate of the barrier height; this yields the MPW54 functional. The MPW54 functional, as well as the MPW60 functional that was previously parametrized for the H+CH(4) reaction, is used with canonical variational theory with small-curvature tunneling to calculate the rate constants for all four ethane reactions from 200 to 2000 K. The final MPW54 calculations are based on curvilinear-coordinate generalized-normal-mode analysis along the reaction path, and they include scaled frequencies and an anharmonic C-C bond torsion. They agree with experiment within 31% for 467-826 K except for a 38% deviation at 748 K; the results for the isotopologs are predictions since these rate constants have never been measured. The kinetic isotope effects (KIEs) are analyzed to reveal the contributions from subsets of vibrational partition functions and from tunneling, which conspire to yield a nonmonotonic temperature dependence for one of the KIEs. The stationary points and reaction-path potential of the MPW54 potential-energy surface are then used to parametrize a new kind of analytical potential-energy surface that combines a semiempirical valence bond formalism for the reactive part of the molecule with a standard molecular mechanics force field for the rest; this may be considered to be either an extension of molecular mechanics to treat a reactive potential-energy surface or a new kind of combined quantum-mechanical/molecular mechanical (QM/MM) method in which the QM part is semiempirical valence bond theory; that is, the new potential-energy surface is a combined valence bond molecular mechanics (CVBMM) surface. Rate constants calculated with the CVBMM surface agree with the MPW54 rate constants within 12% for 534-2000 K and within 23% for 200-491 K. The full CVBMM potential-energy surface is now available for use in variety of dynamics calculations, and it provides a prototype for developing CVBMM potential-energy surfaces for other reactions.     

H atom elimination from the pi sigma* state in the photodissociation of phenol

Photodissociation of phenol at 248 nm was studied using multimass ion imaging techniques. Photofragment translational energy distribution of H atom elimination was measured. The results demonstrate that H atom elimination occurs on the pi sigma(*) excited state which has repulsive potential-energy functions with respect to the stretching of OH bond. It supports the recent ab initio calculation.     

Ab initio calculations of the Ar–ethane intermolecular potential energy surface using bond function basis sets

The intermolecular potential energy surface (PES) of argon with ethane has been studied by ab initio calculations at the levels of second‐order Møller–Plesset perturbation (MP2) theory and coupled‐cluster theory with single, double, and noniterative triple configurations (CCSD(T)) using a series of augmented correlation‐consistent basis sets. Two sets of bond functions, bf1 (3s3p2d) and bf2 (6s6p4d2f), have been added to the basis sets to show a dramatic and systematic improvement in the convergence of the entire PES. The PES of Ar–ethane is characterized by a global minimum at a near T‐shaped configuration with a well depth of 0.611 kcal mol^?1, a second minimum at a collinear configuration with a well depth of 0.456 kcal mol^?1, and a saddle point connecting the two minima. It is shown that an augmented correlation‐consistent basis set with a set of bond functions, either bf1 or bf2, can effectively produce results equivalent to the next larger augmented correlation‐consistent basis set, that is, aug‐cc‐pVDZ‐bf1 ≈ aug‐cc‐pVTZ, aug‐cc‐pVTZ‐bf1 ≈ aug‐cc‐pVQZ. Very importantly, the use of bond functions improves the PES globally, resulting accurate potential anisotropy. Finally, MP2 method is inadequate for accurate calculations, because it gives a potentially overestimated well depth and, more seriously, a poor potential anisotropy. © 2012 Wiley Periodicals, Inc.     

键能的分子轨道理论研究 1: 理论公式

从LCAO-MO出发, 给出了一个计算键能的近似方法, 即EAB(i)-∑∑CaiSabCbiεi为第i个占据分子轨道(MO)中的一对电子对A-B键键能的贡献。对所有分子轨道求和即为该键的键能: EAB=∑EAB(i)。按该方法, 不仅可以计算各种不同分子中每两个相键连原子间的键能, 还可以从MO及AO角度分析每一具体键, 如σ, π, δ键的键能以及各AO对键能的贡献。该方法虽有别于求键焓和平衡离解能De, 但计算结果和De的实验值甚相符合。通过对键能的分析研究, 能较好地揭示原子间的相互作用关系及化学键的强弱, 从而可进一步探讨化学反应活性, 反应速率等化学性质。     

The strength and length of donor-acceptor bonds in molecular complexes

Numerous published data on the structure and thermodynamics of formation of molecular complexes are analyzed. The enthalpies of complexation (−ΔH) are related to the characteristic parameter Δr = [r _DA−a ₁(r _D+r _A)], where r _DA is the donor-acceptor bond length determined by microwave spectroscopy and X-ray analysis, r _D and r _A are the tabulated values of the homopolar covalent radii of the heteroatoms that form the donor-acceptor bond, and a ₁ is an empirical coefficient equal to 0.901±0.007. The relation between −ΔH and Δr values has the form −ΔH = a ₂/Δr (a ₂ = 21.6±1.6 kJ Å mol⁻¹), with a mean relative error of approximation of about 15% and a correlation coefficient of 0.97. As the strength of the complex increases, the donor-acceptor bond length approaches the sum of the heteropolar covalent radii of the atoms involved in the bond (Δr tends to zero). At Δr ≫ 1, the strength of complexes is determined by weak van der Waals interactions between the complex components and the −ΔH values tend to zero. Dedicated to the memory of E. N. Guryanova (1911–2004), a prominent scientist who formulated the basic principles elaborated by her disciples in this review. Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 10, pp. 1869–1878, October, 2007.     

A molecular orbital explanation for the B?N bond shortening in H3BNH3 on going from the gaseous to the solid state

A simple molecular orbital model has been applied to explanation of the B? N bond shortening in H₃BNH₃ on going from the gaseous to the solid state. In this model, the shortening is attributed to the bond order increase that is caused by the fact that each atom in the crystal experiences different external electrostatic potential to each other and thus the orbital energy level of each atom is changed. To illustrate this model, Effective Fragment Potential (EFP) method has been applied to the system consisting of a H₃BNH₃ molecule and 30 dipole moments whose magnitudes are determined by Lorentz's local field theory. This EFP computation has brought significant B? N bond shortening (1.668 → 1.623 Å), which is about 50% of the actual shortening. The factor of the remaining discrepancy has been analyzed by Morokuma decomposition under EFP and localized orbital analysis. These analyses have revealed that the remaining discrepancy is almost compensated by incorporating the dihydrogen bonds (B? H···H? N) that are formed by the orbital interaction between the bonding orbital of the B? H and the antibonding orbital of the N? H. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Conformational analysis of substituted polyacetylenes

A conformational analysis of isolated chains of polymethylacetylene (PMA), polypentylacetylene (PPA), and poly(t-butylacetylene) (PTA) was carried out taking into account interactions between nonbonded atoms and torsional potentials. It was found that the trans configuration of all three polymers is more stable than the cis configuration, the difference in potential energy between the trans and the cis isomers however being very small for PTA, leading to the possibility of observing a trans^→ cis isomerization in some solvents. The calculations show that the substituted polyacetylenes are not found in a planar conformation, the larger deviations from planarity being found with the bulkier substituents: PTA > PPA > PMA. A correlation could be established between the UV absorption limit of the samples and the minimum torsional angle of the potential-energy functions. This relation predicts that the absorption limit is shifted to long wavelengths on increasing the planarity of the molecule. Moreover, UV spectra could be calculated from the potential-energy functions, and it is shown that the potential-energy functions of other substituted polyacetylenes can be calculated from their experimental UV spectra.     

Domain‐Averaged Exchange‐Correlation Energies as a Physical Underpinning for Chemical Graphs

A novel solution to the problem of assigning a molecular graph to a collection of nuclei (i.e. how to draw a molecular structure) is presented. Molecules are universally understood as a set of nuclei linked by bonds, but establishing which nuclei are bonded and which are not is still an empirical matter. Our approach borrows techniques from quantum chemical topology, which showed for the first time the construction of chemical graphs from wave functions, shifting the focus on energetics. This new focus resolves issues surrounding previous topological analyses, in which domain‐averaged exchange‐correlation energies (V_xc), quantities defined in real space between each possible atom pair, hold the key. Exponential decay of V_xc in non‐metallic systems as the intercenter distance increases guarantees a well‐defined hierarchy for all possible V_xc values in a molecule. Herein, we show that extracting the set of atom pairs that display the largest V_xc values in the hierarchy is equivalent to retrieving the molecular graph itself. Notably, domain‐averaged exchange‐correlation energies are transferable, and they can be used to calculate bond strengths. Fine‐grained details resulted to be related to simple stereoelectronic effects. These ideas are demonstrated in a set of simple pilot molecules.     

XMVB: a program for ab initio nonorthogonal valence bond computations

An ab initio nonorthogonal valence bond program, called XMVB, is described in this article. The XMVB package uses Heitler-London-Slater-Pauling (HLSP) functions as state functions, and calculations can be performed with either all independent state functions for a molecule or preferably a few selected important state functions. Both our proposed paired-permanent-determinant approach and conventional Slater determinant expansion algorithm are implemented for the evaluation of the Hamiltonian and overlap matrix elements among VB functions. XMVB contains the capabilities of valence bond self-consistent field (VBSCF), breathing orbital valence bond (BOVB), and valence bond configuration interaction (VBCI) computations. The VB orbitals, used to construct VB functions, can be defined flexibly in the calculations depending on particular applications and focused problems, and they may be strictly localized, delocalized, or bonded-distorted (semidelocalized). The parallel version of XMVB based on MPI (Message Passing Interface) is also available.     

Intramolecular Natural Energy Decomposition Analysis: Applications to the Rational Design of Foldamers

We describe an intramolecular version of the natural energy decomposition analysis (NEDA), with the aim of evaluating interactions between molecular fragments across covalent bonds. The electronic energy in intramolecular natural energy decomposition analysis (INEDA) is divided into electrical, core, and charge transfer components. The INEDA method describes the fragments using the nonfragmented electronic density, and, therefore, there are no limitations in how to choose the boundary orbital. We used INEDA to evaluate the interaction energies that give origin to barriers of rotation around C_amide C_aromatic (C_am C_ar) and N_amide C_aromtaic (N_am C_ar) bonds in arylamide‐foldamer building blocks. We found that differences of barrier height between models with different ortho‐aryl substituents stem from charge transfer and core interactions. In three‐center hydrogen‐bond (H‐bond) models with an NH proton donor H‐bound to two electronegative ortho‐aryl substituents, the interaction energy of the three‐center system is larger than in either of the two‐center H‐bond subsystem alone, indicating an increase of overall rigidity. The combination of INEDA and NEDA allows the evaluation of intermolecular and intramolecular interactions using a consistent theoretical framework. © 2018 Wiley Periodicals, Inc.     

键能的分子轨道理论研究2: 键能与Mulliken布居对键强度的 判断的比较

用键能E~A~B和Mulliken布居对化学键强度的判别进行了分析比较。结果表明,键能判据比Mulliken布居判据所得结论更符合实际情况。作为衡量原子间化学键强度的尺度,不仅应考虑原子轨道间的布居因素,还应考虑分子轨道(或原子轨道)的能量因素。     

取代基对N-H…O=C氢键三聚体中氢键强度的影响

使用MP2方法研究了氢键三聚体中N—H…O=C氢键强度, 探讨了氢键受体分子中不同取代基对N—H…O=C氢键强度的影响. 研究表明, 不同取代基对氢键三聚体中N—H…O=C氢键强度的影响是不同的: 取代基为供电子基团, 氢键键长r(H…O)缩短, 氢键强度增强; 取代基为吸电子基团, 氢键键长r(H…O)伸长, 氢键强度减弱. 自然键轨道(NBO)分析表明, N—H…O=C氢键强度越强, 氢键中氢原子的正电荷越多, 氧原子的负电荷越多, 质子供体和受体分子间的电荷转移越多. 供电子基团使N—H…O=C氢键中氧原子的孤对电子n(O)对N—H的反键轨道滓*(N—H)的二阶相互作用稳定化能增加, 吸电子基团使这种二阶相互作用稳定化能减小. 取代基对与其相近的N—H…O=C氢键影响更大.     

Synthesis, spectroscopic and structural properties of CF3SO2OCCl3

Trichloromethyl trifluoromethanesulphonate, CF₃SO₂OCCl₃, was prepared by quantitative reaction between Ag(CF₃SO₂O) and BrCCl₃. The conformational and structural properties of the gaseous molecule were studied by vibrational spectroscopy (IR (gas, liquid), Raman (liquid) and quantum chemical calculations (DFT and ab initio methods)).Theoretical and experimental vibrational results evidenced the presence of a single conformer with C₁ symmetry. This result is in agreement with the adopted geometry of covalent sulphonates. The conformational preference was studied using the total energy scheme and natural bond orbital partition scheme. Additionally, the total potential-energy has been deconvoluted using six fold decomposition in terms of a Fourier-type expansion.     

Ab initio calculation of the NH(3sigma-)-NH(3sigma-) interaction potentials in the quintet, triplet, and singlet states

We present the ab initio potential-energy surfaces of the NH-NH complex that correlate with two NH molecules in their 3sigma- electronic ground state. Three distinct potential-energy surfaces, split by exchange interactions, correspond to the coupling of the S(A) = 1 and S(B) = 1 electronic spins of the monomers to dimer states with S = 0, 1, and 2. Exploratory calculations on the quintet (S = 2), triplet (S = 1), and singlet (S = 0) states and their exchange splittings were performed with the valence bond self-consistent-field method that explicitly accounts for the nonorthogonality of the orbitals on different monomers. The potential surface of the quintet state, which can be described by a single Slater determinant reference function, was calculated at the coupled cluster level with single and double excitations and noniterative treatment of the triples. The triplet and singlet states require multiconfiguration reference wave functions and the exchange splittings between the three potential surfaces were calculated with the complete active space self-consistent-field method supplemented with perturbative configuration interaction calculations of second and third orders. Full potential-energy surfaces were computed as a function of the four intermolecular Jacobi coordinates, with an aug-cc-pVTZ basis on the N and H atoms and bond functions at the midpoint of the intermolecular vector R. An analytical representation of these potentials was given by expanding their dependence on the molecular orientations in coupled spherical harmonics, and representing the dependence of the expansion coefficients on the intermolecular distance R by the reproducing kernel Hilbert space method. The quintet surface has a van der Waals minimum of depth D(e) = 675 cm(-1) at R(e) = 6.6a0 for a linear geometry with the two NH electric dipoles aligned. The singlet and triplet surfaces show similar, slightly deeper, van der Waals wells, but when R is decreased the weakly bound NH dimer with S = 0 and S = 1 converts into the chemically bound N2H2 diimide (also called diazene) molecule with only a small energy barrier to overcome.     

卤键在超分子化学中的应用进展

卤键是指作用在卤原子(路易斯酸)和具有孤对电子的原子或π电子体系(路易斯碱)之间的新型弱相互作用,其在超分子多维自组装和分子识别(如超分子催化、超分子选择拆分、超分子传感)等领域有着广泛的应用。本文介绍了卤键的类别、特性、功能及在超分子化学结构与功能领域中的应用。