凝固点降低法研究某些脂肪胺的自缔合作用

本工作利用计算机技术改进加热式凝固点测量装置的数据采集系统. 改进后装置的混合物凝固点测量再现性为±0.001 K; 溶剂偏摩尔过量自由能测量再现性为±0.1 J.mol~(-1). 用凝固点降低法测量正丙胺、异丙胺、正丁胺, 异丁胺、正戊胺、异戊胺、二乙胺、二异丙胺、三乙胺、三正丁胺分别与苯和对二甲苯构成二元混合物的溶剂活度系数. 并使用一个包含物理相互作用的三参数连续缔合模型拟合实验数据. 结果表明该模型能较好地刻画所研究的二元缔合混合物。     

Capture and dissociation in the complex-forming CH + H2 → CH2 + H, CH + H2 reactions

The rate coefficients for the capture process CH + H(2)→ CH(3) and the reactions CH + H(2)→ CH(2) + H (abstraction), CH + H(2) (exchange) have been calculated in the 200-800 K temperature range, using the quasiclassical trajectory (QCT) method and the most recent global potential energy surface. The reactions, which are of interest in combustion and in astrochemistry, proceed via the formation of long-lived CH(3) collision complexes, and the three H atoms become equivalent. QCT rate coefficients for capture are in quite good agreement with experiments. However, an important zero point energy (ZPE) leakage problem occurs in the QCT calculations for the abstraction, exchange and inelastic exit channels. To account for this issue, a pragmatic but accurate approach has been applied, leading to a good agreement with experimental abstraction rate coefficients. Exchange rate coefficients have also been calculated using this approach. Finally, calculations employing QCT capture/phase space theory (PST) models have been carried out, leading to similar values for the abstraction rate coefficients as the QCT and previous quantum mechanical capture/PST methods. This suggests that QCT capture/PST models are a good alternative to the QCT method for this and similar systems.     

Redox Titrations-II Location of inflection points of titration curves for homogeneous redox reactions

The relative positions of the inflection points and equivalence point of a homogeneous redox reaction have been studied by using the redox buffer capacity to derive an equation for the titration curve. The position of the inflection point corresponding to the maximum slope of the titration curve relative to the equivalence point depends on the electron transfers of the analyte and the titrant (the stoichiometric coefficients of the reaction equation) and on the difference between the formal potentials of the redox couples in a more complicated way than has been described previously.     

2-氨基-5巯基-1,3,4-噻二唑异构化反应机理的量子化学研究

用密度泛函理论(DFT)中的B3LYP方法,采取6-31+G* *基组对2-氨基-5巯基-1,3,4-噻二唑(简称AMT)的异构化反应机理进行了量子化学研究,全参数优化了异构化过程中反应物、产物的几何构型,找出了异构化途径中的过渡态,并通过振动分析加以确认,同时进行零点能校正.研究结果表明,异构化过程存在六种不同的异构化通道,有六个过渡态,相对而言, A→C之间的异构化反应最易发生,C是最稳定的异构化产物.     

直链烷烃结构型和凝聚型性质的递变规律研究

用图示法和非线性规划的方法对直链烷烃同系物最高成键分子轨道能级、最低成键分子轨道能级,电离电位、氧化半波电位、正常沸点、正常熔点、临界压力、临界温度、密度、折光率、表面张力和粘度等12种结构型性质和凝聚型性质的变化规律进行研究,结果发现直链烷烃结构型性质和凝聚型性质一般能遵守同系对数递变规律,各种结构型性能和凝聚型性能均与对数递变函数呈优良的相关性,相关系数均大于0.99．用同系对数递变规律对直链烷烃各种性能进行预测的结果表明,除少数同系列起始化合物的偏差较大外,大多数预测值与实验值非常吻合,实验散点几乎与对数递变函数曲线相重合．     

Recent computational developments with the self-consistent electron pairs method and application to the stability of glycine conformers

Significant improvements in calculational efficiency and capability with the self-consistent electron pairs (SCEP ) method has resulted from several new computational developments. A new procedure for constructing the important internal Coulomb and exchange operators has substantially reduced the preiteration time. A general scheme for utilizing molecular symmetry has been used to advantage in reducing the number of pair functions and external operators that must be found explicitly at each iteration. A projection operator tool has been implemented and found to be quite effective at minimizing the number of iterations required at some point on a potential energy surface when an SCEP wavefunction exists for some nearby point. These and other improvements in the program construction have yielded sizable reductions in time for some representative test cases, including water and a potential energy curve for formaldehyde. The new SCEP program also performs low-order perturbation theory treatments and coupled electron pair approximation (CEPA ) calculations using the same operator approach. The usefulness of the approach is demonstrated by very large scale calculations on the stability of the two interstellar glycine conformers. These calculations involve the variational treatment of 82,205 symmetry-adapted singly and doubly substituted configurations involving 225 internal electron pairs.     

BrO与CH3SH反应机理的量子化学及拓扑研究

利用密度泛函和电子密度拓扑分析方法对BrO与CH₃SH反应的微观机理进行了理论研究. 在B3LYP/6-311G (d, p)水平上对反应势能面上的各驻点进行几何构型的全优化; 振动分析和IRC计算证实了中间体和过渡态的真实性和相互连接关系; 计算得到了各反应通道的活化能, 并进行了零点能校正. 计算结果表明: 该反应存在7个反应通道, 其中生成CH₃S＋HOBr和CH₃SO＋HBr的通道为主要反应通道. 通过对反应过程中部分驻点的电子密度拓扑分析, 首次发现了接近平面的四元环状过渡态, 从而拓展了原来对环状结构过渡态定义的适用范围.     

Modeling a liquid crystal dynamics by atomistic simulation with an ab initio derived force field

Atomistic molecular dynamics (MD) simulations of 4-n-pentyl 4'-cyano-biphenyl (5CB) have been performed, adopting a specific ab initio derived force field. Two state points in the nematic phase and three in the isotropic phase, as determined in a previous work, have been considered. At each state point, at least 10 ns have been produced, allowing us to accurately calculate single-molecule properties. In the isotropic phase, the values of the translational diffusion coefficient, and even more so the activation energy for the process, agree well with experimental data. Qualitatively, also the dynamic anisotropy of the nematic phase is correctly accounted for. Rotational diffusion coefficients, which describe spinning and tumbling motions, fall well within the range of experimental values. The reorientational dynamics of our model 5CB covers diverse time regimes. The longest one is strongly temperature dependent and characterized by a relaxation time in accord with experimental dielectric relaxation data. Shear viscosity and Landau-de Gennes relaxation times, typically collective variables, reproduce the experimental results very well in the isotropic phase. In the nematic phase, despite a large statistical uncertainty due to the extremely slow relaxation of the correlation functions involved, our simulation yields the correct relative order of the three experimental Miesowicz viscosities.     

Molecular orbital anharmonic estimates for the infrared spectrum of CO2

The vibrational spectrum of CO2 up to second overtones has been calculated at four different ab initio levels using second order perturbation theory equations in a simplified manner, in which just a few cross-terms suitable for numerical estimation are considered in the Taylor series representing the potential energy and dipole moment functions. The series coefficients are obtained through polynomial regression of estimated single point energy and dipole values for a few distorted geometries along each normal coordinate. The effect of Fermi resonance on near-degenerate energy levels was also taken into account through the usual first order perturbation equations. MP2/6-31G(extended) frequency estimates have a root mean square error of just 32.14 cm(-1). This accuracy is achieved partly due to the underestimation of the harmonic frequencies, which compensates for the neglect of the cross-term chi(ij) anharmonic constants. The chi(ii) constants which depend on cubic and quartic energy coefficients are reasonably well estimated at all ab initio levels. The energy coefficient beta(sbb) responsible for the magnitude of the Fermi resonance is estimated with a maximum error of just 13%. Despite the inclusion of anharmonicities, errors for band intensities are still much larger than for the frequencies. Both electrical and mechanical anharmonicities may be equally important to the band intensity.     

Vibrations of H+(D+) in stoichiometric LiNbO3 single crystal

A first principles quantum mechanical calculation of the vibrational energy levels and transition frequencies associated with protons in stoichiometric LiNbO(3) single crystal has been carried out. The hydrogen contaminated crystal has been approximated by a model one obtains by translating a supercell, i.e., a cluster of LiNbO(3) unit cells containing a single H(+) and a Li(+) vacancy. Based on the supercell model an approximate Hamiltonian operator describing vibrations of the proton sublattice embedded in the host crystal has been derived. It is further simplified to a sum of uncoupled Hamiltonian operators corresponding to different wave vectors (ks) and each describing vibrations of a quasi-particle (quasi-proton). The three dimensional (3D) Hamiltonian operator of k=0 has been employed to calculate vibrational levels and transition frequencies. The potential energy surface (PES) entering this Hamiltonian operator has been calculated point wise on a large set of grid points by using density functional theory, and an analytical approximation to the PES has been constructed by non-parametric approximation. Then, the nuclear motion Schro?dinger equation has been solved by employing the method of discrete variable representation. It has been found that the (quasi-)H(+) vibrates in a strongly anharmonic PES. Its vibrations can be described approximately as a stretching, and two orthogonal bending vibrations. The theoretically calculated transition frequencies agree within 1% with those experimentally determined, and they have allowed the assignment of one of the hitherto unassigned bands as a combination of the stretching and the bending of lower fundamental frequency.     

Improved treatment for matrix elements of spin-dependent operators

In this paper a general method for the evaluation of the matrix elements of spin-dependent operators is proposed to improve the treatment primitively suggesteed by Cooper and Musher. This approach is largely based on the recent results which the present authors have achieved in the representation theory for the inner- and outer-product reduction of the symmetric group. It is shown that the so-called outer-product coupling coefficients (OPCC ) can be used to generalize the method for constructing the irreducible tensor operators of group S_n. Together with the use of inner-product coupling coefficients (IPCC ), an expression for the matrix elements of spin-dependent operators is presented as the product of a Racah coefficient for S_n and a reduced matrix element which can be expressed in terms of IPCC, OPCC , and the related integrals. The treatment for one- and two-electron spin-dependent operators is discussed in detail.     

Molecular dynamics investigations of ozone on an ab initio potential energy surface with the utilization of pattern-recognition neural network for accurate determination of product formation

The singlet-triplet transformation and molecular dissociation of ozone (O(3)) gas is investigated by performing quasi-classical molecular dynamics (MD) simulations on an ab initio potential energy surface (PES) with visible and near-infrared excitations. MP4(SDQ) level of theory with the 6-311g(2d,2p) basis set is executed for three different electronic spin states (singlet, triplet, and quintet). In order to simplify the potential energy function, an approximation is adopted by ignoring the spin-orbit coupling and allowing the molecule to switch favorably and instantaneously to the spin state that is more energetically stable (lowest in energy among the three spin states). This assumption has previously been utilized to study the SiO(2) system as reported by Agrawal et al. (J. Chem. Phys. 2006, 124 (13), 134306). The use of such assumption in this study probably makes the upper limits of computed rate coefficients the true rate coefficients. The global PES for ozone is constructed by fitting 5906 ab initio data points using a 60-neuron two-layer feed-forward neural network. The mean-absolute error and root-mean-squared error of this fit are 0.0446 eV (1.03 kcal/mol) and 0.0756 eV (1.74 kcal/mol), respectively, which reveal very good fitting accuracy. The parameter coefficients of the global PES are reported in this paper. In order to identify the spin state with high confidence, we propose the use of a pattern-recognition neural network, which is trained to predict the spin state of a given configuration (with a prediction accuracy being 95.6% on a set of testing data points). To enhance the prediction effectiveness, a buffer series of five points are validated to confirm the spin state during the MD process to gain better confidence. Quasi-classical MD simulations from 1.2 to 2.4 eV of total internal energy (including zero-point energy) result in rate coefficients of singlet-triplet transformation in the range of 0.027 ps(-1) to 1.21 ps(-1). Also, we find very low dissociation probability up to 2.4 eV of internal energy during the investigating period (5 ps), which suggests that dissociation does not occur directly from the singlet ground-state, but it involves the excited triplet-state as an intermediate step and requires more reaction time to occur.     

Variational properties of the discrete variable representation: Discrete variable representation via effective operators

A variational finite basis representation/discrete variable representation (FBR/DVR) Hamiltonian operator has been introduced. By calculating its matrix elements exactly one obtains, depending on the choice of the basis set, either a variational FBR or a variational DVR. The domain of grid points on which the FBR/DVR is variational has been shown to consist of the subsets of the set of grid points one obtains by diagonalizing commuting variational basis representations of the coordinate operators. The variational property implies that the optimal of the subsets of a fixed number of points, i.e., the subset which gives the possible highest accuracy eigenpairs, gives the DVR of the smallest trace. The symmetry properties of the variational FBR/DVR Hamiltonian operator are analyzed and methods to incorporate symmetry into FBR/DVR calculations are discussed. It is shown how the Fourier-basis FBR/DVR suitable to solving periodic systems arise within the theory presented. Numerical examples are given to illustrate the theoretical results. The use of variational effective Hamiltonian and coordinate operators has been instrumental in this study. They have been introduced in a novel way by exploiting quasi-Hermiticity.     

The Diels-Alder reaction of ethene and 1,3-butadiene: an extended multireference ab initio investigation.

The concerted and stepwise mechanisms of the Diels-Alder reaction between 1,3-butadiene and ethene have been investigated using highly correlated multireference methods (MRAQCC) and extended basis sets. Full MRAQCC geometry optimizations have been performed in all cases. The best estimate for the energy barrier of the Diels-Alder reaction is 22 kcalmol(-1). Anti- and gauche-out minima for the biradical structures and corresponding fragmentation saddle points have been determined. The biradical anti fragmentation saddle point is located 6.5 kcalmol(-1) above the concerted saddle point. The gauche-in structure does not correspond to a local minimum, but leads on geometry optimization directly to cyclohexene.     

Three-dimensional numerical integration for electronic structure calculations

Two three-dimensional numerical schemes are presented for molecular integrands such as matrix alements of one-electron operators occuring in the Fock operator and expectation values of one-electron operators describing molecular properties. The schemes are based on a judicious partitioning of space so that product-Gauss integration rules can be used in each region. Convergence with the number of integration points is such that very high accuracy (8–10 digits) may be obtained with obtained with a modest number of points. The use of point group symmetry to reduce the required number of points is discussed. Examples are given for overlap, nuclear potential, and electric field gradient integrals.     

Isopropylcyclopropane + OH gas phase reaction: a quantum chemistry + CVT/SCT approach

A theoretical study of the mechanism and kinetics of the OH hydrogen abstraction from isopropylcyclopropane (IPCP) is presented. Optimum geometries, frequencies and gradients have been computed at the BHandHLYP/6-311++G(d,p) level of theory for all stationary points, as well as for additional points along the minimum energy path (MEP). Energies have been improved by single-point calculations at the above geometries using CCSD(T)/6-311++G(d,p) to produce the potential energy surface. The rate coefficients are calculated for the temperature range 260-350 K by using canonical variational theory (CVT) with small-curvature tunneling (SCT) corrections. Our analysis suggests a stepwise mechanism involving the formation of a reactant complex in the entrance channel and a product complex in the exit channel, for all the modeled paths. The reactant complexes are examined in detail, because they exhibit alkene-like structure. The excellent agreement between the overall calculated and experimental rate coefficients at 298 K supports the reliability of the parameters obtained for the temperature dependence and branching ratios of the IPCP + OH reaction, proposed here for the fist time. The expression that best describes the studied reaction is k(overall) = 6.15 x 10(-13)e1747/RT cm3 x molecule(-1) x s(-1). The predicted activation energy is -0.89 kcal/mol.     

Thermodynamic studies of drug-alpha-cyclodextrin interactions in water at 298.15 K: promazine hydrochloride/chlorpromazine hydrochloride + alpha-cyclodextrin + H(2)O systems

Data on osmotic coefficients have been obtained for a binary aqueous solution of two drugs, namely, promazine hydrochloride (PZ) and chlorpromazine hydrochloride (CPZ) using a vapor pressure osmometer at 298.15 K. The observed critical micelle concentration (cmc) agrees excellently with the available literature data. The measurements are extended to aqueous ternary solutions containing fixed a concentration of alpha-cyclodextrin (alpha-CD) of 0.1 mol kg(-1) and varied concentrations (approximately 0.005-0.2 mol kg(-1)) of drugs at 298.15 K. It has been found that the cmc values increase by the addition of alpha-CD. The mean molal activity coefficients of the ions and the activity coefficient of alpha-CD in binary as well as ternary solutions were obtained, which have been further used to calculate the excess Gibbs free energies and transfer Gibbs free energies. The lowering of the activity coefficients of ions and of alpha-CD is attributed to the existence of host-guest (inclusion)-type complex equilibria. It is suggested that CPZ forms 2:1 and 1:1 complexed species with alpha-CD, while PZ forms only 1:1 complexed species. The salting constant (ks) values are determined at 298.15 K for promazine-alpha-CD and chlorpromazine-alpha-CD complexes, respectively, by following the method based on the application of the McMillan-Mayer theory of virial coefficients to transfer free energy data. It is noted that the presence of chlorine in the drug molecule imparts better complexing capacity, the effect of which gets attenuated as a result of hydrophobic interaction. The results are discussed from the point of view of associative equilibria before the cmc and complexed equilibria for binary and ternary solutions, respectively.     

Ab initio-based global double many-body expansion potential energy surface for the first 2A" electronic state of NO2

An ab initio-based global double many-body expansion potential energy surface is reported for the first electronic state of the NO(2)((2)A") manifold. Up to 1700 ab initio energies have been employed to map the full configuration space of the title molecule, including stationary points and asymptotic channels. The calculated grid energies have been scaled to account for the incompleteness of the basis set and truncation of the MRCI expansion and fitted analytically with a total root-mean-squared-deviation smaller than 1.0 kcal mol(-1). The lowest point of the potential energy surface corresponds to the (2)B(1) linear minimum, which is separated from the C(s) distorted minimum by a C(2v) barrier of ≈9.7 kcal mol(-1). As usual, the proposed form includes a realistic representation of long-range interactions. Preliminary work indicates its reliability for reaction dynamics.     

Boiling points for five binary systems of sulfolane with aromatic hydrocarbons at 101.33 kPa

Boiling points have been determined at 101.33 kPa for the binary mixtures of sulfolane+o-xylene, sulfolane+m-xylene, sulfolane+p-xylene, sulfolane+ethylbenzene and sulfolane+1,2,4-trimethylbenzene. Calculations of the non-ideality of the vapor phase were made with the second virial coefficients evaluated from the Hayden–O’Connell method. The binary parameters for five activity coefficient models (Margules, van Laar, Wilson, NRTL and UNIQUAC) have been fitted with the experimental boiling points measured in this work. A comparison of model performances has been carried out using the criterion of the average absolute deviations in boiling point. The activity coefficient of the component in the liquid phase is discussed based on the UNIFAC model with the consideration of the dipole–dipole interactions.