Complexation of Phenol and Thiophenol by Amine N‐Oxides: Isothermal Titration Calorimetry and ab Initio Calculations

To develop a new solvent‐impregnated resin (SIR) system for removal of phenols from water, the complex formation of dimethyldodecylamine N‐oxide (DMDAO), trioctylamine N‐oxide (TOAO), and tris(2‐ethylhexyl)amine N‐oxide (TEHAO) with phenol (PhOH) and thiophenol (PhSH) is studied. To this end we use isothermal titration calorimetry (ITC) and quantum chemical modeling (on B3LYP/6‐311G(d,p)‐optimized geometries: B3LYP/6‐311+G(d,p), B3LYP/6‐311++G(2d,2p), MP2/6‐311+G(d,p), and spin component scaled (SCS) MP2/6‐311+G(d,p); M06‐2X/6‐311+G(d,p)//M06‐2X/6‐311G(d,p), MP2 with an extrapolation to the complete basis set limit (MP2/CBS), as well as CBS‐Q). The complexes are analyzed in terms of structural (e.g., bond lengths) and electronic elements (e.g., charges). Furthermore, complexation and solvent effects (in benzene, toluene, and mesitylene) are investigated by ITC measurements, yielding binding constants K, enthalpies ΔH⁰, Gibbs fre energies ΔG⁰, and entropies ΔS⁰ of complex formation, and stoichiometry N. The ITC measurements revealed strong 1:1 complex formation between both DMDAO–PhOH and TOAO–PhOH. The binding constant (K=1.7–5.7×10⁴ M ^?1) drops markedly when water‐saturated toluene was used (K=5.8×10³ M ^?1), and π–π interaction with the solvent is shown to be relevant. Quantum mechanical modeling confirms formation of stable 1:1 complexes with linear hydrogen bonds that weaken on attachment of electron‐withdrawing groups to the amine N‐oxide moiety. Modeling also showed that complexes with PhSH are much weaker than those with PhOH, and in fact too weak for ITC determination. CBS‐Q incorrectly predicts equal or even higher binding enthalpies for PhSH than for PhOH, which invalidates it as a benchmark for other calculations. Data from the straightforward SCS‐MP2 method without counterpoise correction show very good agreement with the MP2/CBS values.     

C6H5—H…X分子间氢键的理论计算

利用密度泛函理论B3LYP方法, 在6-311+G(3df,2p)水平上对C6H5—H…X型分子间氢键进行了几何构型优化、氢键相互作用能、电子密度分布等计算. 其中C6H6为质子供体, HCOH、H2O、NH3、CH2NH和HCN为质子受体. 从电荷布居分析、自然键轨道等角度详细地讨论了C6H5—H…X 体系中, 共轭π键、O和N的不同键型结构对氢键形成的影响以及孤电子对与C—H 反键轨道之间的相互作用(n→σ*)等.     

Neutral homoaromaticity in some heterocyclic systems

[reaction: see text] Neutral homoaromaticity has been evaluated in heterocyclic systems related to the bicyclo[3.2.1]octane skeleton with replacement of CH(2) at C-2 in bicyclo[3.2.1]octa-3,6-diene with X = BH, AlH, Be, Mg, O, S, PH, NH (12); replacement of CH at C-3 in bicyclo[3.2.1]octa-3,6-dien-2-yl anion with PH, S, NH, O (13); and replacement at C-2 and C-3 with N and O (14). Stabilization energies (SE) are evaluated using density functional theory and homodesmotic equations at the B3LYP/6-311+G(3df,2p)//B3LYP/6-31G(d) level for series 12-14. Stabilization energies are compared with diamagnetic susceptibility exaltations, Lambda, CSGT-B3LYP/6-31G(d)//B3LYP/6-31G(d), and nucleus-independent chemical shifts (NICS), GIAO-B3LYP/6-311+G(2d,p)//B3LYP/6-31G(d). Analysis of frontier orbitals and geometries, B3LYP/6-31G(d)//B3LYP/6-31G(d), and proton affinities of 2-azabicyclo[3.2.1]octa-3,6-diene, pyrrole, and related model systems, B3LYP/6-311+G(2d,2p)//B3LYP/6-31G(d), provide complementary evidence supporting the division of the substrates evaluated into antihomoaromatic (12, X = BH, AlH, and Be), nonhomoaromatic (12, X = O, S, NH, PH), and homoaromatic (13, X = S, PH, NH, O and 14 X = ON), with 12 (X = Mg) appearing as transitional between anti- and nonhomoaromatic.     

BeH2与HX(X=F,Cl, Br,I)形成的二氢键复合物的结构特征与本质

对BeH2与HX(X=F, Cl, Br, I)形成的二氢键复合物的结构特征及本质进行了探讨. 在MP2/6-311++G(3d,3p)水平优化、频率验证, 得到复合物的分子结构, 用分子间距离及电子密度拓扑理论确认BeH2与卤化氢已形成了二氢键型复合物. 在MP2/6-311++G(3d, 3p)水平下进行基函数重叠误差(BSSE)校正后的结合能在-14.468 kJ·mol-1到-5.464 kJ·mol-1之间.用对称匹配微扰理论(SAPT)对复合物的结合能进行分解, 结果表明, BeH2…HX二氢键复合物中静电能对总吸引能的贡献都是最主要的, 但交换排斥能、诱导能、色散能对总结合能的贡献也很重要. 从BeH2…HF到BeH2…HI, 诱导能对总吸引能的贡献从37.8%逐渐减小到24.0%. 而色散能对总吸引能的贡献从BeH2…HF体系中的16.0%逐渐增加到BeH2…HI体系中的33.8%.     

SiH4与HX(X=F,Cl,Br,I)形成的二氢键复合物的结构特征及本质

本文就SiH4与HX形成的二氢键复合物的结构特征及本质进行了探讨。在MP2/6-311++G(3d,3p)水平优化、频率验证得到复合物的分子结构,通过分子间距离及电子密度等值线图,我们确认SiH4与卤化氢已形成了二氢键复合物。MP2/6-311++G(3d,3p)水平下进行BSSE校正后的结合能为2.703－4.439 KJ/mol。用对称匹配微绕理论(SAPT)对结合能进行分解,分解结果显示,SiH4匟X(X=F,Cl,Br,I)二氢键复合物中静电能对总吸引能的贡献小于28％,并且相对稳定,这就是说SiH4匟X二氢键复合物的本质并非静电作用,而是静电能、诱导能、色散能、交换能对总结合能的贡献都非常重要。     

An evaluation of various computational methods for the treatment of organoselenium compounds

A reliable computational method for the prediction of organoselenium geometries and bond dissociation energies (BDEs) has been determined on the basis of the performance of density functional theory (DFT: B3LYP and B3PW91) and ab initio molecular orbital procedures (Hartree-Fock (HF)) in conjunction with various Pople basis sets including (but not limited to) the 6-31G(d), 6-31G(d,p), 6-311G(d), 6-311G(d,p), 6-311G(2df,p), and 6-311G(3df,3pd) sets. Predicted geometries and BDEs are compared with available experimental data and quadratic configuration interaction including single and double substitutions (QCISD) results. The B3PW91/6-311G(2df,p) level of theory is recommended for the prediction of the geometries and energetics of organoselenium compounds.     

Structural and Energetic Effects in the Molecular Recognition of Acetylated Amino Acids by 18-Crown-6

Absolute 18-crown-6 (18C6) binding affinities of four protonated acetylated amino acids (AcAAs) are determined using guided ion beam tandem mass spectrometry techniques. The AcAAs examined in this work include: N-terminal acetylated lysine (N_???CAcLys), histidine (N_???CAcHis), and arginine (N_???CAcArg) as well as side chain acetylated lysine (N_???CAcLys). The kinetic-energy-dependent cross sections for collision-induced dissociation (CID) of the (AcAA)H⁺(18C6) complexes are analyzed using an empirical threshold law to extract absolute 0 and 298?K (AcAA)H⁺?18C6 bond dissociation energies (BDEs) after accounting for the effects of multiple collisions, kinetic and internal energy distributions of the reactants, and unimolecular dissociation lifetimes. Theoretical electronic structure calculations are performed to determine stable geometries and energetics for neutral and protonated 18C6 and the AcAAs as well as the proton bound complexes of these species, (AcAA)H⁺(18C6), at the B3LYP/6-311+G(2d,2p)//B3LYP/6-31?G* and M06/6-311+G(2d,2p)//B3LYP/6-31G* levels of theory. For all four (AcAA)H⁺(18C6) complexes, loss of neutral 18C6 corresponds to the most favorable dissociation pathway. At elevated energies, products arising from sequential dissociation of the primary CID product, H⁺(AcAA), are also observed. Protonated N_???CAcLys exhibits a greater 18C6 binding affinity than other protonated N_???CAcAAs, suggesting that the side chains of Lys residues are the preferred binding sites for 18C6 complexation to peptides and proteins. N_???CAcLys exhibits a greater 18C6 binding affinity than N_???CAcLys, suggesting that binding of 18C6 to the side chain of Lys residues is more favorable than to the N-terminal amino group of Lys.     

脯氨酸的构象及性质

用X3LYP法在6-311++G(d, p)和6-311++G(3df, 3pd)基组水平上对脯氨酸15种构象的几何结构、相对能量、电子结构、红外光谱、偶极矩、极化率等性质进行了研究, 并与PBE1PBE/6-311++G(d, p)的结果和文献相比较, 从而得到: (1) 的脯氨酸的15种构象中能量最低的有4种, 不同构象中存在着强弱不同的5种氢键, 其中以N…H—O氢键最强, 并存在特殊的C—H…O=C氢键. 两种方法计算的几何结构数据相近, 均与实验值吻合; (2) 在构象相对能差计算方面, X3LYP具有明显的优势, 用中等基组就可以得到与高水平从头算法和大基组相同的结果, 而PBE1PBE法计算的相对能值则相差较大; (3) 脯氨酸不同构象中偶极矩最大和极化率最小的是最稳定的构象1和2, 两种方法计算的结果一致.     

传统氢键(X—H…Y)还是X—H…π相互作用？ ——无机苯与卤化氢相互作用本质的量子化学研究

在MP2/6-311＋＋G**水平对无机苯(B3N3H6)与卤化氢HX (X＝F, Cl, Br, I)相互作用体系进行了系统研究. 结果表明在B3N3H6-HX (X＝F, Cl, Br, I)体系的平衡几何结构中, HX的H原子倾向于指向B3N3H6环上的N原子, 且从HF到HI相互作用强度依次减弱. 与苯-卤化氢体系比较, 除与HF相互作用B3N3H6较C6H6强外, 其余体系B3N3H6均较C6H6弱(结合能数值相差4 kJ/mol左右). 对称匹配微扰理论(SAPT)能量分解结果说明静电、诱导和色散力对描述B3N3H6-卤化氢体系的相互作用都很重要, 从HF到HI静电能占总吸引作用能的百分比逐渐减少, 色散能占总吸引作用能的百分比逐渐增加, 这种变化趋势与苯-卤化氢体系比较类似, 表明B3N3H6与卤化氢的相互作用随着卤素原子序数的递增, 传统氢键作用趋势减弱, X—H…π相互作用趋势增强.     

Conformational energies for 2-substituted butanes

The conformational free energies for some 2-substituted butanes where X = F, Cl, CN, and CCH were calculated using G3-B3, CBS-QB3, and CCSD(T)/6-311++G(2d,p) as well as other theoretical levels. The above methods gave consistent results with free energies relative to the trans conformers as follows: X = CCH, g+ = 0.77 +/- 0.05 kcal/mol. g- = 0.88 +/- 0.05 kcal/mol; X = CN, g+ = 0.85 +/- 0.05 kcal/mol, g- = 0.75 +/- 0.05 kcal/mol; X = Cl, g+ = 0.70 +/- 0.05 kcal/ml, g- = 0.80 +/- 0.05 kcal/mol; and X = F, g+ = 0.53 +/- 0.05 kcal/mol, g- = 0.83 +/- 0.05 kcal/mol. The conformational free energies also were estimated using the observed liquid phase IR spectra and intensities calculated using B3LYP/6-311++G** and MP2/6-311++G**. The rotational free energy profiles for all of the compounds were estimated at the G3-B3 level.     

Experimental and theoretical studies of sodium cation interactions with the acidic amino acids and their amide derivatives

The binding of Na+ to aspartic acid (Asp), glutamic acid (Glu), asparagine (Asn), and glutamine (Gln) is examined in detail by studying the collision-induced dissociation (CID) of the four sodiated amino acid complexes with Xe using a guided ion beam tandem mass spectrometer (GIBMS). Analysis of the energy-dependent CID cross sections provides 0 K sodium cation affinities for the complexes after accounting for unimolecular decay rates, internal energy of the reactant ions, and multiple ion-molecule collisions. Quantum chemical calculations for a number of geometric conformations of each Na+(L) complex are determined at the B3LYP/6-311+G(d,p) level with single-point energies calculated at MP2(full), B3LYP, and B3P86 levels using a 6-311+G(2d,2p) basis set. This coordinated examination of both experimental work and quantum chemical calculations allows the energetic contributions of individual functionalities as well as steric influences of relative chain lengths to be thoroughly explored. Na+ binding affinities for the amide complexes are systematically stronger than those for the acid complexes by 14 +/- 1 kJ/mol, which is attributed to an inductive effect of the OH group in the carboxylic acid side chain. Additionally, the Na+ binding affinity for the longer-chain amino acids (Glx) is enhanced by 4 +/- 1 kJ/mol compared to the shorter-chain Asx because steric effects are reduced.     

Diels-Alder Reactivity of a Chiral Anthracene Template with Symmetrical and Unsymmetrical Dienophiles: A DFT Study

In this work, we used Density Functional Theory calculations to assess the factors that control the reactivity of a chiral anthracene template with three sets of dienophiles including maleic anhydrides, maleimides and acetoxy lactones in the context of Diels-Alder cycloadditions. The results obtained here (at the M06-2X/6-311++G(d,p) level of theory) suggest that the activation energies for maleic anhydrides and acetoxy lactones are dependent on the nature of the substituent in the dienophile. Among all studied substituents, only −CN reduces the energy barrier of the cycloaddition. For maleimides, the activation energies are independent of the heteroatom of the dienophile and the R group attached to it. The analysis of frontier molecular orbitals, charge transfer and the activation strain model (at the M06-2X/TZVP level based on M06-2X/6-311++G(d,p) geometries) suggest that the activation energies in maleic anhydrides are mainly controlled by the amount of charge transfer from the diene to the dienophile during cycloaddition. For maleimides, there is a dual control of interaction and strain energies on the activation energies, whereas for the acetoxy lactones the activation energies seem to be controlled by the degree of template distortion at the transition state. Finally, calculations show that considering a catalyst on the studied cycloadditions changes the reaction mechanism from concerted to stepwise and proceed with much lower activation energies.     

Transition state structure and energetics of the N(2)O + X (X = Cl,Br) reactions

The structural and vibrational properties of the transition state of the N(2)O + X (X = Cl,Br) reactions have been characterized by ab initio methods using density functional theory. We have employed Becke's hybrid functional (B3LYP), and transition state optimizations were performed with 6-31G(d), 6-311G(2d,2p), 6-311+G(3d,2p), and 6-311+G(3df,2p) basis sets. For the chlorine atom reaction the coupled-cluster method (CCSD(T)) with 6-31G(d) basis set was also used. All calculations resulted in transition state structures with a planar cis arrangement of atoms for both reactions. The geometrical parameters of transition states at B3LYP are very similar, and the reaction coordinates involve mainly the breaking of the N-O bond. At CCSD(T)/6-31G(d) level a contribution of the O-Cl forming bond is also observed in the reaction coordinate. In addition, several highly accurate ab initio composite methods of Gaussian-n (G1, G2, G3), their variations (G2(MP2), G3//B3LYP), and complete basis set (CBS-Q, CBS-Q//B3LYP) series of models were applied to compute reaction energetics. All model chemistries predict exothermic reactions. The G3 and G2 methods result in the smallest deviations from experiment, 1.8 and 0 kcal mol(-1), for the enthalpies of reaction for N(2)O reaction with chlorine and bromine, respectively. The G3//B3LYP and G1 methods perform best among the composite methods in predicting energies of the transition state, with a deviation of 1.9 and 3.0 kcal mol(-1), respectively, in the activation energies for the above processes. However, the B3LYP/6-311+G(3df,2p) method gives smaller deviations of 0.4 and -1.0 kcal mol(-1), respectively. The performance of the methodologies applied in predicting transition state energies was analyzed.     

Hydrogen bonding interaction in sarcosine–water complex using ab initio and DFT method

The hydrogen bonding interaction in the Sarcosine (N‐methylglycine)–water complex is studied using ab initio, MP2, and density functional theory (DFT/B3LYP). For this complex, binding energies, dipole–dipole interactions, chemical hardness, and chemical potential have been calculated. Three different basis sets, viz. 6‐311+G, 6‐311++G, and 6‐311++G*, have been used to optimize the geometries by all three methods. The basis set superposition errors are also calculated, and the corrected binding energies are reported for this complex. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Cooperative enhancement of water binding to crownophane by multiple hydrogen bonds: analysis by high level ab initio calculations.

The intermolecular interaction energy of the model system of the water-crownophane complex was analyzed. The water molecule has four hydrogen bonds, with the two hydrogen-donating phenolic hydroxy groups and two hydrogen-accepting oxygen atoms of the poly-oxyethylene chain of the crownophane in the complex. The MP2/6-311G(2d,2p) level calculations of the model system of the complex (hydrogen donating unit + hydrogen accepting unit + water) indicate that the binding energy of the water is 21.85 kcal/mol and that the hydrogen bond cooperativity increases the binding energy as much as 3.67 kcal/mol. The calculated interaction energies depend on the basis set, while the basis set dependence of the cooperative increment is negligible. Most of the cooperative increment is covered by the HF level calculation, which suggests that the major source of the hydrogen bond cooperativity in this system has its origin in induction. The BLYP/6-311G** and PW91/6-311G** level interaction energies of the model system are close to the MP2/6-311G** interaction energies, which suggests that the DFT calculations with these functionals are useful methods to evaluated the interactions of hydrogen bonded systems.     

H2CO与双卤分子间卤键的电子密度分析

运用B3LYP和MP2方法在6-311++G(d,p)基组水平上, 对H2CO-XY(XY=F2、Cl2、Br2、ClF、BrF、BrCl)卤键体系进行构型全优化, 得到了O…X—Y型卤键复合物. 结果表明, MP2/6-311++G(d,p)计算结果与实验值较吻合. 并在MP2水平下计算了分子间的相互作用能, 用完全均衡校正CP(counterpoise procedure)方法对基函数重叠误差(BSSE)进行了校正. 利用电子密度拓扑分析方法对卤键复合物的电子密度拓扑性质进行了分析研究.     

甲醛与乙腈相互作用的理论研究

在HF/6-311G(d,p)、 MP2/6-311G(d,p)和B3LYP/6-311G(d,p)水平上,对H2CO和CH3CN以及设计的4种结构H2CO…CH3CN复合物等进行几何全优化和振动频率计算,排除振动频率为负值的非局域极小点结构,并对稳定的环状构型复合物结合能进行基组重叠误差校正和零点振动能校正.分子间相互作用的能量分解分析显示,静电能在H2CO...CH3CN相互作用能量中占主导地位,电荷转移能居第二位.     

DFT computational studies on rotation barriers,tautomerism, intramolecular hydrogen bond,and solvent effects in 8‐hydroxyquinoline

Intramolecular hydrogen binding interactions in 8‐hydroxyquinoline, both in its zwitterionic tautomer and in the rotamer without the intramolecular hydrogen bond (IHB), have been computed using the B3LYP and MPW1K density functionals. The rotation of the O? H bond and intramolecular proton transfer reactions were studied theoretically. The following theory levels have been applied: B3LYP/6‐31G(d,p), B3LYP/6‐311++G(d,p), MPW1K/6‐311++G(d,p), and MPW1K/6‐311++G(2d,3p)//MPW1K/6‐311++G(d,p). Natural bond orbital (NBO) analysis has also been carried out. The effect of medium (benzene, chloroform, tetrahydrofuran, 1,2‐dichloroethane, acetone, water) was simulated using the self‐consistent reaction field (SCRF) method within the framework of the polarizable continuum model (PCM), at the MPW1K/6‐311++G(d,p) level. The evolution of geometry, relative energies, heights of rotation (around the O? H bond) and tautomerization barriers, IHB energies, and ΔG(solv) have been systematically investigated. The results obtained have shown the failure to neglect some changes of the above characteristics in polar media with respect to the gaseous phase. The series of stability of the forms under study in the gaseous phase remains the same in solution. Thus, in spite of the important role of the solvent electrostatic effects, the intrinsic stability of those species overcomes the solvent effects. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Theoretical investigation of the weakly dihydrogen bonded dimers H2 − nXnAlHHArF and H2 − nXnAlHHKrF (n = 0–2; X = F, Cl)

An ab initio computational study of the properties of the dihydrogen-bonded complexes of H_{2 − n}X_nAlH (n = 0–2; X = F, Cl) with the rare gas (Rg) compounds HArF and HKrF was carried out at the MP2(full)/6-311++G(2d,2p) level of theory. For all the studied complexes, we found a low zero-point corrected binding energies. Large red shifts of the H–Rg vibrational stretching frequency in both complexes were predicted. Electrostatic interactions between the individual monomers are also predicted.