1-烷基-4-氨基-1,2,4-三唑硝酸盐含能离子液体的蒸气压及其摩尔汽化焓的理论研究

合成了1-烷基-4-氨基-1,2,4-三唑硝酸盐含能离子液体([RATZ]NO₃),并通过核磁和红外进行了结构表征;采用Gaussian09/B3LYP/6-311+G(d,p)密度泛函理论,计算了[RATZ]NO₃的离子间相互作用能及摩尔体积;在298 K-323 K温度范围内,测定了不同配比[RATZ]NO₃-EtOH混合溶液的饱和蒸气压,其中乙醇摩尔分数分别为0.984、0.996、0.999以及1.000。系统研究了[RATZ]NO₃-EtOH混合溶液的饱和蒸气压、[RATZ]NO₃的蒸气压及摩尔汽化焓与温度、离子间相互作用能以及结构之间的关系。结果表明：[RATZ]NO₃-EtOH混合溶液的饱和蒸气压随着温度的升高、离子间作用能的减小以及阳离子体积的增大而增大,其沸点比纯溶剂高,且在298 K-323 K温度范围内[RATZ]NO₃的平衡蒸气压均低于250 mPa,因此说明含能离子液体具有不挥发性,蒸气压极低,并通过理论计算得到的离子间相互作用能及体积,解释了[RATZ]NO₃的摩尔汽化焓随烷基链增长而降低的原因。     

醚基功能化离子液体[MOEMIm]Cl和[EOEMIm]Cl热力学性质

通过核磁共振氢谱,核磁共振碳谱,元素分析和热重分析对醚基功能化的离子液体[MOEMIm]Cl和[EOEMIm]Cl进行了表征。在温度范围T=288.15–328.15 K内,测定了离子液体[MOEMIm]Cl和[EOEMIm]Cl的密度(ρ)、表面张力(γ)和折光率(nD)。根据这些实验数据,讨论并计算了离子液体[MOEMIm]Cl和[EOEMIm]Cl的体积性质。计算出离子液体[MOEMIm]Cl和[EOEMIm]Cl的摩尔表面吉布斯自由能(gs)、摩尔表面熵(s)、摩尔表面焓(h)、摩尔极化度(Rm)和摩尔极化率(αp),h均近似为一个常数说明这两种离子液体从内部到表面的过程是一个等库仑过程,同时这两种离子液体的Rm和αp均与温度无关。本文还用摩尔表面Gibbs自由能改进Lorentz-Lorenz方程并预测离子液体表面张力,预测值与实验值高度相关。     

硝酸1-烷基-4-氨基-1,2,4-三唑鎓含能离子液体的蒸气压及其摩尔汽化焓的理论研究

合成了硝酸1-烷基-4-氨基-1,2,4-三唑鎓含能离子液体([RATZ]NO3),并通过核磁和红外进行了结构表征;采用Gaussian09/B3LYP/6-311+G(d,p)密度泛函理论计算了[RATZ]NO3的离子间相互作用能及摩尔体积;在298K～323K温度范围内,测定了不同配比[RATZ]NO3-Et OH混合溶液的饱和蒸气压。研究了[RATZ]NO3-Et OH混合溶液的饱和蒸气压、[RATZ]NO3的蒸气压及摩尔汽化焓与温度、离子间相互作用能以及结构之间的关系。结果表明,[RATZ]NO3-Et OH混合溶液的饱和蒸气压随温度的升高、离子间作用能的减小以及阳离子体积的增大而增大,其沸点比纯溶剂高,且在298K～323K温度范围内[RATZ]NO3的平衡蒸气压均低于250mPa,说明含能离子液体具有不挥发性;通过理论计算得到离子间相互作用能及体积,解释了[RATZ]NO3的摩尔汽化焓随烷基链增长而降低的原因。     

HTPB固体推进剂增塑剂选取分子模拟研究

固体推进剂中增塑剂要求同粘合剂体系相容性良好,并提高体系的低温性能.本文采用分子动力学(MD)方法,首先计算了端羟基聚丁二烯(HTPB)粘合剂及增塑剂癸二酸二辛酯(DOS)、己二酸二辛酯(DOA)、壬酸异癸酯(TOA)、邻苯二甲酸二丁酯(DBP)和邻苯二甲酸二辛酯(DOP)的溶度参数,以此从相容性角度选取推进剂增塑剂;计算数值基本吻合实验值,表明常用的增塑剂从相容性都能满足要求.其次模拟获取了HTPB及HTPB/增塑剂混合体系的比体积-温度关系得到了体系的玻璃化转变温度(Tg),揭示增塑剂对HTPB体系低温性能的影响.结果显示:(1)HTPB的Tg模拟值为202K,基本吻合实验值196K.(2)HTPB/DOS混合体系中,当增塑剂DOS的质量含量从12%、22%、29%到36%(摩尔含量分别为50%、66%、75%和90%)增加时,体系的Tg线性降低;TOA和DOP增塑的粘合剂体系(摩尔含量为75%)Tg也降低,而增塑剂DOA和DBP对体系的Tg影响不大.因此,基于相容性及提高粘合剂低温性能考虑,DOS、DOA和DOP作为HTPB的增塑剂优于TOA和DBP.     

Nd(NO_3)_3-HNO_3-H_2O体系的过量吉布斯自由能和Nd(NO_3)_3-HNO_3-H_2O-HDEHP-AMSCO体系水相中各组分活度系数的计算

我们简化了混合电解质溶液过量吉布斯自由能的Pitzer公式。用此简化公式并结合我们以往的工作计算了Nd(NO_3)_3-HNO_3-H_2O体系的过量吉布斯自由能和Nd(NO_3)_3-HNO_3-H_2O-HDEHP-AMSCO体系水相中各组分的活度系数,计算值与实验值均很接近。此工作为镧系金属萃取分离进行热力学分析提供了方便。     

季铵盐型双子表面活性剂与十八醇的混合单分子膜

研究了双子表面活性剂12-2-16和12-2-12分别与十八醇(C18H37OH)在空气-水界面上混合单分子膜的π-A等温线. 在相分离表面压以下, 比较了不同表面压下和不同混合比单分子膜的混合表面过剩自由能ΔGMexo, 分析了双子表面活性剂与脂肪醇在空气-水界面上混合膜中的相容性. 结果表明, 12-2-16与C18H37OH在所有混合摩尔比下随着表面压增高, 自由能增大. 12-2-12与C18H37OH混合膜体系的相容性取决于两者的混合比, ΔGMexo随所加入C18H37OH摩尔分数的增加逐渐增大, 从异种分子间净的吸引作用转变到相互排斥作用体系, 转变点为C18H37OH加入量的摩尔分数0.65. 当混合为热力学自发过程时, 增大表面压将有利于混合; 而对相互排斥体系, 增加表面压将使体系内异种分子之间的相互排斥作用更大.     

热解条件对LiFePO4/C表面自由能的影响

以葡萄糖为碳前驱体对水热合成的磷酸铁锂粉末进行碳包覆.测定磷酸铁锂粉末与三种探测液的接触角,并据此通过Young方程采用Wu方法计算了粉末的表面自由能.研究了表面自由能与电导率、放电比容量的关系及对工业化涂布性能的影响.结果表明,在电导率相同的前提下,粉末表面自由能中色散分量与极性分量的比值(γd/γp)对低倍率放电性能没有影响,但对高倍率放电性能影响较大.提高热解温度和延长热解时间可使γd/γp值增加,这有利于粉体与聚偏氟乙烯(PVdF)的粘合和高倍率放电容量的提高.     

磷酯酰乙醇胺与硬脂酸和十八醇二元混合体系单层膜的热力学特性和AFM观测

利用表面压力-平均分子面积(π-A)曲线的关系, 分别研究了在水/空气界面上形成的磷酯酰乙醇胺(PE)与硬脂酸(SA)和十八醇(OD)二元混合体系的热力学特性. 根据对弹性模量(C_S^-1)、过量分子面积(A_ex)以及表面过量吉布斯自由能(ΔG_ex)等热力学参数的数据计算, 定量分析了混合单层膜分子之间的相互作用. 实验结果表明, PE/SA和PE/OD单层膜两种组分之间摩尔比对其热力学参数有影响. 热力学分析表明, 在X_SA=0.2, 0.8和X_OD=0.8处, PE/SA和PE/OD二元混合体系的热力学参数[过量分子面积(A_ex)和表面过量吉布斯自由能(ΔGex)]相对理想状态均表现为负偏差作用, 说明分子之间的作用为引力作用. 相反, PE/OD二元体系在X_OD=0.2, 0.4, 0.6处, 表现为正偏差作用, 这说明分子之间的作用为斥力作用. AFM观测为PE/SA和PE/OD单层膜热力学特性的理论分析提供了有力的支持.     

293.15K时咪唑醋酸盐-乙醇二元体系的体积性质及分子间相互作用

用比重瓶法测定了293.15 K时1-甲基咪唑醋酸盐([Mim]Ac)/1,3-二甲基咪唑醋酸盐([Mmim]Ac)/1-乙基-3-甲基咪唑醋酸盐([Emim]Ac)-乙醇(EtOH)二元体系在全组成范围内的密度.计算出[Mim]Ac/[Mmim]Ac/[Emim]Ac和EtOH的表观摩尔体积和体系的超额摩尔体积.用三参数多项式关联拟合了表观摩尔体积与摩尔分数的关系,外推出组分的极限偏摩尔体积和摩尔体积.[Mim]Ac/[Mmim]Ac/[Emim]Ac和EtOH的摩尔体积的外推值与实验值分别在±0.07和±0.04 cm3/mol范围内相一致.计算出了[Mim]Ac/[Mmim]Ac/[Emim]Ac和EtOH分别在无限稀溶液中的溶剂化系数.用Redlich-Kister方程关联拟合了超额摩尔体积与摩尔分数的关系.分别根据极限偏摩尔体积、摩尔体积与极限偏摩尔体积的差值、溶剂化系数和超额摩尔体积对照讨论了分子间相互作用的强弱.结果显示,在[Mim]Ac/[Mmim]Ac/[Emim]Ac的浓度无限稀溶液中,[Mim]Ac/[Mmim]Ac/[Emim]Ac-EtOH分子对间相互作用的强弱顺序为[Mim]Ac-EtOH[Mmim]Ac-EtOH[Emim]Ac-EtOH;在EtOH的浓度无限稀溶液中,以及体系中[Mim]Ac/[Mmim]Ac/[Emim]Ac的摩尔分数在0.15~0.95间时,[Mim]Ac/[Mmim]Ac/[Emim]Ac-EtOH分子对间相互作用的强弱顺序都为[Emim]Ac-EtOH[Mmim]Ac-EtOH[Mim]Ac-EtOH.     

化学平衡的本质——浅谈混合熵对化学平衡的贡献

化学平衡是物理化学的重要内容之一。化学反应的等温方程式,给出了平衡常数与标准摩尔反应Gibbs自由能的关系。然而,仅从化学反应等温方程式还很难理解混合熵对化学平衡贡献的本质。本文建立了理想气体体系Gibbs自由能与反应进度的关系,给出了三个理想气体反应体系Gibbs自由能与反应进度的具体关系,分析了混合熵对化学平衡的本质贡献。     

Thermodynamics of mixing of polymers

A method of calculating free energy of mixing of two polymers has been proposed. From a study of the sorption of organic solvents vapour by different polymers and their mixtures, Gibbs free energy of mixing for polymer-solvent systems are calculated. The free energy of mixing of two polymers has been calculated by the method proposed. Heats of mixing of polymers have been measured with a Calvet type microcalorimeter. Entropy of mixing has been calculated. Stabilities of polymer-polymer systems have been discussed.     

"Exact" surface free energies of iron surfaces using a modified embedded atom method potential and lambda integration

Previously a new universal lambda-integration path and associated methodology was developed for the calculation of "exact" surface and interfacial free energies of solids. Such a method is in principle applicable to any intermolecular potential function, including those based on ab initio methods, but in previous work the method was only tested using a relatively simple embedded atom method iron potential. In this present work we apply the new methodology to the more sophisticated and more accurate modified embedded atom method (MEAM) iron potential, where application of other free- energy methods would be extremely difficult due to the complex many-body nature of the potential. We demonstrate that the new technique simplifies the process of obtaining "exact" surface free energies by calculating the complete set of these properties for the low index surface faces of bcc and fcc solid iron structures. By combining these data with further calculations of liquid surface tensions we obtain the first complete set of exact surface free energies for the solid and liquid phases of a realistic MEAM model system. We compare these predictions to various experimental and theoretical results.     

Exploring reaction pathways with transition path and umbrella sampling: application to methyl maltoside

The transition path sampling (TPS) method is a powerful approach to study chemical reactions or transitional properties on complex potential energy landscapes. One of the main advantages of the method over potential of mean force methods is that reaction rates can be directly accessed without knowledge of the exact reaction coordinate. We have investigated the complementary nature of these two differing approaches, comparing transition path sampling with the weighted histogram analysis method to study a conformational change in a small model system. In this case study, the transition paths for a transition between two rotational conformers of a model disaccharide molecule, methyl beta-D-maltoside, were compared with a free energy surface constrained by the two commonly used glycosidic (phi,psi) torsional angles. The TPS method revealed a reaction channel that was not apparent from the potential of mean force method, and the suitability of phi and psi as reaction coordinates to describe the isomerization in vacuo was confirmed by examination of the transition path ensemble. Using both transition state theory and transition path sampling methods, the transition rate was estimated. We have estimated a characteristic time between transitions of approximately 160 ns for this rare isomerization event between the two conformations of the carbohydrate. We conclude that transition path sampling can extract subtle information about the dynamics not apparent from the potential of mean force method. However, in calculating the reaction rate, the transition path sampling method required 27.5 times the computational effort than was needed by the potential of mean force method.     

基于扩散势能和自由体积的自扩散系数模型

A new model for self-diffusion coe±cients was proposed based on both the concepts of molecular free volume and activation energy. The unknown parameters of this model were clearly defined and compared with the Chapman-Enskog model. At the same time a new method for calculating activation energy was devised and applied to the new model. In addition, the free volume was defined by implementing the generic van der Waals equation of state, the radial distribution function of which was obtained by using the MorsaliGoharshadi empirical formula. Under the same conditions, the new model was better than the original free volume model.     

Chemical potential perturbation: a method to predict chemical potentials in periodic molecular simulations

A new method, called chemical potential perturbation (CPP), has been developed to predict the chemical potential as a function of density in periodic molecular simulations. The CPP method applies a spatially varying external force field to the simulation, causing the density to depend upon position in the simulation cell. Following equilibration the homogeneous (uniform or bulk) chemical potential as a function of density can be determined relative to some reference state after correcting for the effects of the inhomogeneity of the system. We compare three different methods of approximating this correction. The first method uses the van der Waals density gradient theory to approximate the inhomogeneous Helmholtz free energy density. The second method uses the local pressure tensor to approximate the homogeneous pressure. The third method uses the Triezenberg-Zwanzig definition of surface tension to approximate the inhomogeneous free energy density. If desired, the homogeneous pressure and Helmholtz free energy can also be predicted by the new method, as well as binodal and spinodal densities of a two-phase fluid region. The CPP method is tested using a Lennard-Jones (LJ) fluid at vapor, liquid, two-phase, and supercritical conditions. Satisfactory agreement is found between the CPP method and an LJ equation of state. The efficiency of the CPP method is compared to that for Widom's method under the tested conditions. In particular, the new method works well for dense fluids where Widom's method starts to fail.     

Application of a volume-translated Peng-Robinson equation of state on vapor-liquid equilibrium calculations

A volume-translated Peng-Robinson (VTPR) equation of state (EOS) is developed in this study. Besides the two parameters in the original Peng-Robinson equation of state, a volume correction term is employed in the VTPR EOS. In this equation, the temperature dependence of the EOS energy parameter was regressed by an improved expression which yields better correlation of pure-fluid vapor pressures. The volume correction parameter is also correlated as a function of the reduced temperature. The VTPR EOS includes two optimally fitted parameters for each pure fluid. These parameters are reported for over 100 nonpolar and polar components. The VTPR EOS shows satisfactory results in calculating the vapor pressures and both the saturated vapor and liquid molar volumes. In comparison with other commonly used cubic EOS, the VTPR EOS presents better results, especially for the saturated liquid molar volumes of polar systems. VLE calculations on fluid mixtures were also studied in this work. Traditional van der Waals one-fluid mixing rules and other mixing models using excess free energy equations were employed in the new EOS. The VTPR EOS is comparable to other EOS in VLE calculations with various mixing rules, but yields better predictions on the molar volumes of liquid mixtures.     

Prediction of thermodynamic properties of polymer solutions by a group-contribution method

The Flory–Huggins lattice-theory expression for solvent activity in a polymer-solution is commonly used to calculate the thermodynamic interaction parameter χ with the aid of experimental data from vapor pressure osmometry. This expression assumes that χ is independent of composition. However, experimental data for a variety of polymer-solvent mixtures indicate that χ exhibits an appreciable concentration dependence. A group contribution method, UNIFAC (UNIQUAC Functional-Group Activity Coefficients) incorporating the free-volume correction of Oishi and Prausnitz is used to predict the dependence of χ on solvent concentration. Agreement with previously reported experimental data is within 15%. Calculated values of χ obtained from the Flory–Huggins expression for solvent activity and from the corresponding Gibbs free energy of mixing (which does not assume that χ is independent of composition) are compared. Calculations based on the Gibbs free energy of mixing predict a somewhat larger value of χ relative to those based on solvent activity. The specific Gibbs free energy of mixing for polystyrene-solvent mixtures is calculated using the UNIFAC model, and is found to represent qualitatively the phase equilibrium behavior. Quantitative discrepancies are observed, however, for the polystyrene-acetone system in light of the actual experimental solubility reported by Suh and Clark (20). Most of the thermodynamic predictions for polymer-solvent systems investigated herein are correlated qualitatively with the relative mismatch between solubility parameters of both components.     

Temperature and molecular weight dependence of interfacial tension between immiscible polymer pairs by the square gradient theory combined with the Flory–Orwoll–Vrij equation-of-state theory

Interfacial tension between immiscible polymer pairs was predicted by using a square gradient theory in conjunction with the Flory–Orwoll–Vrij equation-of-state expression for the free energy of mixing. The contact interaction parameter was determined by fitting the equation-of-state theory to experimental cloud points taken from the literature, and the square gradient coefficient was estimated from the relation derived from a scattering function. The modified square gradient theory could successfully predict both the magnitude and temperature dependence of interfacial tension between polystyrene and poly(methyl methacrylate), although no adjustable parameters were used in calculating interfacial tension. The molecular weight dependence of interfacial tension was also successfully predicted. The contribution of free volume on interfacial tension is analyzed for two systems: polystyrene/poly(methyl methacrylate) and polystyrene/poly(dimethyl siloxane) blends. © 1998 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys 36: 2683–2689, 1998     

Activating Mixing in Preparation and Modification of Polymeric Composites

The concept of activating mixing as a procedure for preparing and modification of polymeric composites was developed. A new classification and generalized parameters of mixing, characterizing conditions of the energy supply to polymeric composition through the device working elements, were proposed. The dependence for determining the treatment modes providing controllable development of chemical transformations in a polymer system, which is invariant with respect to treatment modes and working installation, was obtained.