Unique stoichiometric representation for computational thermochemistry

Evaluation of the enthalpy of formation of species via quantum chemical methods, as well as the evaluation of their performance, is mainly based on single reaction schemes, i.e., reaction schemes that involve a minimal number of reference species where minimal means that, if a reference species is omitted, there is no way to write a balanced reaction scheme involving the remaining species. When the number of reference species exceeds the minimal number, the main problem of computational thermochemistry is inevitably becoming an optimization problem. In this communication we present an exact and unique solution of the optimization problem in computational thermochemistry along with a stoichiometric interpretation of the solution. Namely, we prove that the optimization problem may be identically solved by enumerating a finite and unique set of reactions referred to as group additivity (GA) response reactions (RERs).     

An exact stoichiometric representation of the resonance energy

The energetic measure of aromaticity usually referred to as resonance energy (RE) is shown to possess a remarkable stoichiometric interpretation. Namely, the RE may be uniquely partitioned into a linear sum of contributions associated with group additivity (GA) response reactions (RERs). This new result is a powerful tool for critical analyses of various energetic approaches to RE. In particular, the single reaction scheme approach that is routinely used to evaluate RE is shown to be a particular case of the general GA method.     

Thermal intramolecular electron transfer reactions of [N,N′-ethylene-bis(salicylideneaminato)](2,4-pentanedionato)cobalt(III) complexes in the solid state

The title complexes with the formula Co(salen)L where L is a series of 2,4-pentanedionates underwent thermally induced one-electron transfer reactions from L to Co(III). The reaction left behind a stoichiometric amount of the crystalline Co^II(salen) complex which took up oxygen in a molar ratio of Co:O₂ = 2:1. The kinetic analyses showed that the electron transfer reaction rate was apparently dominated by activation entropy rather than by activation enthalpy.     

Interaction between oppositely charged polyelectrolytes in aqueous solution

Oppositely charged polyelectrolytes interact in solution, forming polyelectrolyte complexes, which often appear as gel-like precipitates. This kind of complex formation was studied by means of calorimetric and rheological measurements. The enthalpy effects, though being fairly small, give some information about the binding strength of counterions to the macroion. We studied the system poly(p-styrene sulfonate)/poly(trimethylammonium-2-ethyl methacrylate) (PSS-PTMA), varying systematically the low molar mass counterions of PSS. In every case, the maximum of enthalpy was found around a 1:1 (mol:mol monomer units) composition of the complexes, with the shape of enthalpy versus composition-curve indicating a stoichiometric interaction. The maximum enthalpy decreased with increasing atomic mass of the counterion when the alkaline metal salts of PSS were used and no change was made on the side of the cationic polyelectrolyte. The salts of the alkaline earth metals gave a distinctly higher enthalpy. On the contrary, viscosity measurements showed a very broad minimum as a function of composition, indicating that the formation of non-stoichiometric complexes is also occurring. The conclusion of these observations is that the complex formation is stoichiometric with respect to the monomeric units, but not necessarily stoichiometric with respect to the entire macromolecules.     

Identification of general linear relationships between activation energies and enthalpy changes for dissociation reactions at surfaces

The activation energy to reaction is a key quantity that controls catalytic activity. Having used ab inito calculations to determine an extensive and broad ranging set of activation energies and enthalpy changes for surface-catalyzed reactions, we show that linear relationships exist between dissociation activation energies and enthalpy changes. Known in the literature as empirical Br?nsted-Evans-Polanyi (BEP) relationships, we identify and discuss the physical origin of their presence in heterogeneous catalysis. The key implication is that merely from knowledge of adsorption energies the barriers to catalytic elementary reaction steps can be estimated.     

Stoichiometric zirconium oxide cations as potential building blocks for cluster assembled catalysts

Employing guided-ion-beam mass spectrometry, we identified a series of positively charged stoichiometric zirconium oxide clusters that exhibit enhanced activity and selectivity for three oxidation reactions of widespread chemical importance. Density functional theory calculations reveal that these clusters all contain the same active site consisting of a radical oxygen center with an elongated zirconium-oxygen bond. Calculated energy profiles demonstrate that each oxidation reaction is highly favorable energetically and involves easily surmountable barriers. Furthermore, the active stoichiometric clusters may be regenerated by reacting oxygen-deficient clusters with a strong oxidizer. This indicates that these species may promote multiple cycles of oxidation reactions and, therefore, exhibit true catalytic behavior. The stoichiometric clusters, having structures that resemble specific sites in bulk zirconia, are promising candidates for potential incorporation into a cluster assembled catalyst material.     

反相高效液相色谱中溶质保留过程的热力学研究 Ⅱ. 焓变及其分量, 吸附及解吸附焓变

依据液相色谱中溶质计量置换保留模型, 从理论上分别推导出了在反相高效液相色谱(RP-HPLC)中lgI(与1mol溶质对固定相的亲合势有关的常数)和Z(对应于1mol溶质被固定相吸附时释放出置换剂的摩尔数)对绝对温度导数的线性关系式。发现非极性与极性小分子溶质的lgI和Z均具有热力学衡常数的特征。如在前报[1]报告的吉布斯自由能变一样, 也推导出了可将溶质保留过程中的总吸附焓变△H(Pa)分成两个独立的分量, 吸附焓变△(1,a)和解吸附焓变△H(Z,D)的表达式。可同时对在流动相中不同置换剂浓度条件下的△H(Pa)和△(Z,D)进行估算,并与实验值进行了比较, 偏差小于±10%。并用△H(1,a)和△H(Z,D)值对溶质在RP-HPLC的保留过程中热变化进行了解释。     

Estimation of kinetic parameters from thermochemical data

A method is described by which kinetic parameters may be calculated from the measured temperature changes caused by the heat produced during a chemical reaction. An isoperibol titration calorimeter with an ampoule-breaking facility is used to obtain the temperature data. The temperature changes resulting from the reaction between tri-isopropyl phosphite and sulphur (S(8)) are used as an example to demonstrate the effectiveness of the method. The temperature changes are used to calculate an enthalpy of reaction. From the enthalpy of reaction and intermediate heats, instantaneous concentrations of the reactants may be calculated.     

Retro-[3 + 2]-cycloaddition reactions in the decomposition of five-membered nitrogen-containing heterocycles

Nitrogen-containing heterocycles form the basis for a new generation of high-energy density materials, and they serve as model compounds for nitrogen-containing fuels, such as coal and biomass, and they form the backbone of ionic liquids. A novel retro-[3 + 2]-cycloaddition to a three-membered diene and a two-membered dienophile, analogous to a retro-Diels-Alder reaction, may constitute an important initial reaction step in the thermal decomposition of these heterocyclic compounds. We investigate the kinetics and thermodynamics of these reactions for the heterocycles pyrrole, pyrazole, imidazole, 1,2,3-triazole, 1,2,4-triazole, 1,2,5-triazole, 1,3,4,-triazole, 1H-tetrazole, and 2H-tetrazole, using theoretical computational chemistry. The retro-cycloadditions are shown to form one of the three-membered products: hydrazoic acid (NH=N=N), nitrilimine (NH=N=CH), carbodiimide (NH=C=NH), or ketenimine (NH=C=CH2) plus one of the two-membered products acetylene, hydrogen cyanide, or N2. Accurate enthalpies of formation are calculated for the reaction products using the high-level W1 computational protocol, providing the previously undetermined enthalpy values of 70.09, 88.75, 35.03, and 44.28 kcal mol(-1) for hydrazoic acid, nitrilimine, carbodiimide, and ketenimine, respectively. We apply a variable-order form of the Marcus equation to the dissociation reactions in correlating the enthalpy of reaction with the activation enthalpy. Typical molecular elimination reactions from the heterocycles proceed with an intrinsic activation enthalpy of 36.8 kcal mol(-1) and intrinsic activation free energy of 42.1 kcal mol(-1). However, dissociation reactions resulting in the formation of either NH=C=NH or NH=C=CH2 demonstrate intrinsic barriers ca. 30 kcal mol(-1) higher, as a result of a concerted intramolecular hydrogen shift. Rate constants calculated between 300 and 3000 K indicate that the proposed dissociation reactions should be important in the decomposition of tetrazole and 1,2,3-triazole. This is confirmed by comparison with available experimental data. Decomposition of 1,2,4-triazole to HCN + nitrilimine may also be important at high temperatures. From extrapolation of our Marcus equation relationship, we predict pentazole to decompose to N2 + NHNN with an activation enthalpy of 19.5 kcal mol-1 and a half-life of only 14 s at 298 K.     

The cure of elastomers by dicumyl peroxide as observed in differential scanning calorimetry

Differential scanning calorimetry (DSC) has been used to determine the enthalpy changes accompanying the thermal homolysis and subsequent radical reactions occuring in the dicumyl peroxide (DCP) cure of elastomers (NR, EPDM, SBR, BR, NBR, and EPM). The thermal degradation of DCP alone, dispersed on kaolin clay in a hydrocarbon solvent, was also studied.The degradation of DCP alone results in an exothermic enthalpy change of ～ 215 kJ mole^?1. In elastomer systems, the observed enthalpy exchange can be ordered BR >; SBR ; NBR (34% ACN) > NBR (27% ACN) > NR > EPDM ～ EPM.Low enthalpy values are associated with systems containing predominantly secondary hydrogen atoms (EPDM, EPM). The high enthalpy of cure for BR appears associated with the known ability of the polybutadienyl radicals to propagate a limited chain reaction. The relative extent of cumyloxy radical disproportionation to hydrogen abstraction was determined in each system and this has no apparent effect on the observed enthalpy change.The method of Borchardt and Daniels was used for data reduction and calculation of apparent activation energies. For DCP degradation alone the calculated activation energy is in good agreement with literature values. In elastomer systems, the calculated activation energies must be treated with caution because, as pointed out by Borchardt and Daniels, their method does not apply to solid state reactions.     

计算饱和醇异构体标准生成焓的新方法

基于极性叠加原理,在成功设计烷烃异构体和多氯代烷烃生成焓计算新方法的基础上,设计了一种计算多元醇异构体生成焓的新方法,并合理地假定任一异构体的原子化焓等于三种键(C-C、C-H和C-O-H键)的键能、极性叠加能项以及氢键能项的加和.用这一模型拟合24种原子化焓数据,得到了标准生成焓的估算公式.为了检验预测的精确性,又设计了一种预测方法,使用在排除被预测的化合物条件下回归得到的参数,预测该化合物的生成焓.按这种方法,预测了24种异构体的生成焓.通过该5参数预测的相对于实验值的各种误差(平均绝对误差、均方根误差和最大绝对误差)不仅比7参数的基团法预测的对应误差小得多,而且比相应实验数据的误差还要小.与键加和法比较,该方法的模型包含了极性叠加能和氢键能量,该两项代表了主要的非键相互作用能,表征了不同异构体的结构差异,并大大减少了参数.     

A genetic algorithm encoded with the structural information of amino acids and dipeptides for efficient conformational searches of oligopeptides

The genetic algorithm (GA) is an intelligent approach for finding minima in a highly dimensional parametric space. However, the success of GA searches for low energy conformations of biomolecules is rather limited so far. Herein an improved GA scheme is proposed for the conformational search of oligopeptides. A systematic analysis of the backbone dihedral angles of conformations of amino acids (AAs) and dipeptides is performed. The structural information is used to design a new encoding scheme to improve the efficiency of GA search. Local geometry optimizations based on the energy calculations by the density functional theory are employed to safeguard the quality and reliability of the GA structures. The GA scheme is applied to the conformational searches of Lys, Arg, Met‐Gly, Lys‐Gly, and Phe‐Gly‐Gly representative of AAs, dipeptides, and tripeptides with complicated side chains. Comparison with the best literature results shows that the new GA method is both highly efficient and reliable by providing the most complete set of the low energy conformations. Moreover, the computational cost of the GA method increases only moderately with the complexity of the molecule. The GA scheme is valuable for the study of the conformations and properties of oligopeptides. © 2016 Wiley Periodicals, Inc.     

Thermal behaviors of electrolytes in lithium-ion batteries determined by differential scanning calorimeter

Lithium-ion batteries have been widely used in daily electric appliance, but abusive accidents occurred from time to time. Lithium-ion batteries composed of various electrolytes (containing organic solvents, inorganic salts), binder, and electrode materials, which may be unstable under some abnormal conditions. The formulation or components of electrolytes played a crucial factor in affecting reactive hazards within the cell. In order to meet safety requirements of the lithium-ion batteries on electronic device, reseachers and scholars are continuing to do further studies on the important issues in relation to the lithium-ion batteries. This study aims to apply the differential scanning calorimeter for measuring the thermal or reactive hazards of the organic solvents related to the formulation of the electrolytes. Besides, thermal instabilities of lithiated graphite or deposited lithium with electrolytes were simulated by the reactions between metallic lithium (Li) and organic carbonates of linear and cyclic structures. Exothermic onset temperatures and enthalpy changes were measured and analyzed. Results showed that the thermal behaviors of these organic carbonates themselves or with lithium hexafluorophosphate liberated less enthalpy changes. However, violent exothermic reactions were detected between the linear or cyclic carbonates and lithium metal with larger enthalpy change larger than 1,000 J g^?1 of lithium. This can be observed by Li reacted with diethyl carbonate at surprisingly lower onset temperature about 46.6 °C.     

反相高效液相色谱中溶质保留过程的热力学研究Ⅲ: 熵变及其分量,吸附及解吸附熵变

从理论上推导出了用柱相比和不用柱相比时溶质在反相高效液相色谱保留过程中的总熵变△S~(~P~a~)的两个表达式, 后者亦为计量置换保留模型中线性参数与该过程热力学函数间的定量表达式。△S~(~P~a~)也能够被分成两个独立的分量,吸附熵变△S~(~Ⅰ~,~a~)和解吸附熵变△S~(~Z~,~D~)。依据液相色谱中溶质计量置换保留模型中容量因子与流动相中置换剂活度a~D间的定量关系, 对在不同a~D条件下熔质的总熵变、吸附及解吸附熵变进行了估算, 并用实验测定的总熵变△S~(~P~a~,~e~)与估算总熵变△S~(~P~a~,~c~)进行了比较, 获得了结果的一致性。本文对传统的"吸附伴随着熵减小, 相反, 解吸伴随着熵增大"的定性规律予以定量的说明。     

Modulation of Prins Cyclization by Vibrational Strong Coupling

Light‐molecule strong coupling has emerged within the last decade as a new method to control chemical reactions. A few years ago it was discovered that chemical reactivity could be altered by vibrational strong coupling (VSC). Only a limited number of reactions have been investigated under VSC to date, including solvolysis and deprotection reactions. Here the effect of VSC on a series of aldehydes and ketones undergoing Prins cyclization, an important synthetic step in pharmaceutical chemistry, is investigated. A decrease of the second‐order rate constant with VSC of the reactant carbonyl stretching groups is observed. We also observe an increased activation energy due to VSC, but proportional changes in activation enthalpy and entropy, suggesting no substantive change in reaction pathway. The addition of common cycloaddition reactions to the stable of VSC‐modified chemical reactions is another step towards establishing VSC as a genuine tool for synthetic chemistry.     

Modulation of Prins Cyclization by Vibrational Strong Coupling

Light-molecule strong coupling has emerged within the last decade as a new method to control chemical reactions. A few years ago it was discovered that chemical reactivity could be altered by vibrational strong coupling (VSC). Only a limited number of reactions have been investigated under VSC to date, including solvolysis and deprotection reactions. Here the effect of VSC on a series of aldehydes and ketones undergoing Prins cyclization, an important synthetic step in pharmaceutical chemistry, is investigated. A decrease of the second-order rate constant with VSC of the reactant carbonyl stretching groups is observed. We also observe an increased activation energy due to VSC, but proportional changes in activation enthalpy and entropy, suggesting no substantive change in reaction pathway. The addition of common cycloaddition reactions to the stable of VSC-modified chemical reactions is another step towards establishing VSC as a genuine tool for synthetic chemistry.     

反相高效液相色谱中溶质保留过程的热力学研究Ⅰ: 柱相比, 吉氏自由能变及其分量, 吸附和解吸附自由能变

从热力学观点给反相高效液相色谱(RP-HPLC)中的相比以新的定义, 提出了测定该相比的方法并准确测定了RP-HPLC中溶质计量置换过程中的吉氏自由能变△G~(Pα)。发现在RP-HPLC中相比对保留过程中溶质总自由能变的贡献值几乎可与溶质本身的保留值相当。依据液相色谱中溶质计量置换保留模型。将溶质的△G~(Pα)分成与溶质吸附及解吸有关的两项独立的自由能变分量。△G~(Ⅰ,α)和△G~(Z,D)。对不同流动相组成条件下的△G~(Pα)和△G~(Z,D)进行了估算, 并与实验结果进行了比较, 偏差一般小于±2%。     

Molar volume, excess enthalpy, and Prigogine-Defay ratio of some silicate glasses with different (P,T) histories

Structural relaxation in silicate glasses with different (p,T) histories was experimentally examined by differential scanning calorimetry and measurements of molar volume under ambient pressure. Temperature and pressure-dependent rates of changes in molar volume and generation of excess enthalpy were determined for sodium trisilicate, soda lime silicate, and sodium borosilicate (NBS) compositions. From the derived data, Prigogine-Defay ratios are calculated and discussed. Changes of excess enthalpy are governed mainly by changes in short-range structure, as is shown for NBS where boron coordination is highly sensitive to pressure. For all three glasses, it is shown how the relaxation functions that underlie volume, enthalpy, and structural relaxation decouple for changes in cooling rates and pressure of freezing, respectively. The magnitude of the divergence between enthalpy and volume may be related to differences in structural sensitivity to changes in the (p,V,T,t) space on different length scales. The findings suggest that the Prigogine-Defay ratio is related to the magnitude of the discussed decoupling effect.     

MM3 parametrization of four‐ and five‐coordinated rhenium complexes by a genetic algorithm—Which factors influence the optimization performance?

Genetic algorithms (GA) were used to solve one of the multidimensional problems in computational chemistry, the optimization of force field parameters. The correlation between the composition of the GA, its parameters (p(c), p(m)) and the quality of the results were investigated. The composition was studied for all combinations of a Simple GA/Steady State GA with a Roulette Wheel/Tournament Selector using different values each for crossover (0.5, 0.7, 0.9) and mutation rates (0.01, 0.02, 0.05, 0.10, 0.20). The results show that the performance is strongly dependent on the GA scheme, where the Simple GA/Tournament Selector yields the best results. Two new MM3 parameters were introduced for rhenium compounds with coordination number four (204) and coordination number five (205), the formal oxidation states of rhenium ranging from +V to +VII. A manifold of parameters (Re-C, N, O, S) was obtained by using a diverse set of CSD structures. The advantage of the GA vs. UFF calculations is shown by comparison of several examples. The GA optimized parameters were able to reproduce the geometrical data of the X-ray structures.     

Interfacial thermodynamics of gallic acid adsorption on a chargeable hydrophobic surface

The thermodynamics of adsorption of gallic acid (GA, 3,4,5-trihydroxylbenzoic acid) on the hanging mercury drop electrode (HMDE) surface was studied by temperature-dependent stripping voltammetry (TD-SV), at physiological pH 7.4. The thermodynamic parameters, e.g., Gibbs free energy, ΔG(ADS), enthalpy, ΔΗ(ADS) and entropy, ΔS(ADS), of adsorption have been determined at physiological temperatures 2-40 °C. Chemisorption of the radical species ≡[GA(OH)(2)(O(-))]* is the energetically important reaction. The thermodynamic data show a complex mechanism of adsorption of GA on the electrode surface, which is strongly dependent on temperature. At low-temperatures T<12 °C, adsorption is controlled by enthalpy, while at T>22 °C, adsorption is entropy driven. In the temperature range 12 °C and 22 °C, a combined enthalpy-entropy stabilization occurs. A mechanism is proposed which analyses the implication of thermodynamics to the interfacial adsorption of polyphenols with cell membranes under physiological conditions.