路径积分方法计算聚合物溶液热力学性质

提出了把聚合物溶液的自由能分解为 3部分 ,平动部分、构象部分和热修正项部分 .其中后两部分具有相对于链段的可加性 ,在引入了“有效浓度”的概念后可以用路径积分方法计算 .当采用体积分数代替有效浓度时 ,就回到了李晓毅和赵得禄的溶液理论 .还利用本理论研究了不同分子量聚苯乙烯 环己烷体系的相分离曲线 ,同Flory Huggins(FH)理论相比 ,大大提高了同实验数据的符合程度     

高分子的良溶剂和不良溶剂的两种定义辩析

介绍了目前国际学术界存在的对高分子溶液体系的溶剂优劣性的两种判断标准 ,即相互作用参数χ=1 /2和χ=0 ,前者基于混合自由能 ,后者基于混合热。作者阐明了其倾向于采用后者的观点 ,而后者在目前国内所有高分子物理学教科书中均未被介绍。     

评估石杉碱甲杂合体衍生物作为AD症药物活性顺序的简单方法

为了合理化设计、指导具有高活性、低毒性的石杉碱甲杂合体(Huprine X)衍生物的合成, 建立了一个计算结合自由能的新方案, 以实现对该类衍生物活性排列顺序的预测. 该结合自由能由三部分组成: 抑制剂和靶酶的相互作用能; 活性位点残基在结合前后的构象绝对自由能的变化; 抑制剂从稳定构象变为活性构象的自由能增加值. 通过计算已合成的14个杂合体衍生物的结合自由能, 结果显示理论值和实验测定的生物活性值有很好的等级相关性, 其斯皮尔曼相关系数为0.85, 证明了该方法的可行性.     

修正的格子空间的密度泛函理论在狭缝中的应用

对描述单原子分子溶液在狭缝中的吸附现象的格子空间的密度泛函理论 (LDFT, lattice density functional theory)进行了修正, 在系统Helmholtz函数的推导中引入了平均场近似校正和Gibbs-Helmholtz方程. 对比Monte Carlo (MC)模拟结果, 发现LDFT理论对吸附分子在狭缝中的吸附浓度分布的预测与模拟数据有较大的偏差, 而修正模型的结果与模拟数据吻合较好 .随着体相浓度的变化,分子在狭缝中具有多级吸附行为, 具体表现为在特定体相浓度区, 对相同的体相浓度,狭缝中同时存在不同的分子浓度分布, 而在Gibbs等温线上可以明显看出多级吸附的性质. 对比修正前后的结果发现,两者均可以预测多级吸附行为, 但仍存在着较大的差异.     

高分子的良溶剂和不良溶剂的两种定义辨析

介绍了目前国际学术界存在的对高分子溶液体系的溶剂优劣性的两种判断标准，即相互作用参数ｘ＝１／２和ｘ＝０，前者基于混合自由能，后者基于混合热。作者阐明了其倾向于采用后者的观点，而后者在目前国内所有高分子物理学教科书中均未被介绍。     

糖与非电解质在甲酰胺中相互作用的Gibbs自由能参数

用气液色谱法测定了298.15K时10个非电解质溶质在不同浓度甲酰胺糖(葡萄糖,果糖,蔗糖)溶液中的保留参数.依据McMillanMayer理论,求得非电解质溶质与糖间的Gibbs自由能相互作用参数g_ij和g_ijj.采用OkamotoWood的基团加合法,求得各基团的Gibbs自由能相互作用参数G_x,y,讨论了非电解质溶质与糖间的相互作用以及溶剂性质的影响.     

分子动力学模拟研究两性离子结构对生物分子谷胱甘肽正常构象的维持

用分子动力学方法模拟了正常构象下谷胱甘肽分子(GSH)以及谷胱甘肽分子与两性离子材料表面和疏水性材料表面相互作用的过程. 均方根偏移、构象伸缩性和二面角分布结果显示, 两性离子材料表面能够较好地维持谷胱甘肽分子的正常构象.     

类蛋白质分子的二级结构对力学行为的影响

采用Kolinski等建立的类蛋白质分子的格点模型,研究了由典型的(HHPPHPP)x重复单元构成、含有α螺旋结构的类蛋白质分子链在拉伸过程中的构象性质和力学行为.发现不同强度的α螺旋相互作用会直接影响其拉伸过程.α螺旋相互作用强的类蛋白质分子链,具有更低的内能,更小的应力,在拉伸过程中更容易失去紧密接触对,同时也更容易被拉成“棒状”结构,但在整个拉伸过程中,α螺旋结构且能保持稳定;还发现类蛋白质分子的链长对拉伸也有影响,对较长的类蛋白质分子链,其内能更低,弹性力更小,自由能更大,紧密接触对的含量比例也更高,而“棒状化”程度较小.这些研究能够帮助我们加深对蛋白质分子的构象和弹性力学行为的理解.     

Gibbs吸附公式在固-气和固-液界面吸附中的应用

Gibbs吸附公式本适用于各种界面的吸附研究 ,但由于有固体参与形成的界面界面能测定困难 ,至今Gibbs公式多只用于气 液界面的研究。本文讨论了Gibbs吸附公式在固 气和固 液界面吸附研究中的应用 ,主要是介绍将Gibbs公式与二维气体状态方程结合导出吸附等温式 ,因吸附而引起的固 气和固 液界面能变化规律 ,以及吸附过程中分子状态变化的有关信息     

单嘧磺隆晶体-活性构象转换的分子动力学模拟

应用分子动力学模拟方法对单嘧磺隆在水、正辛醇和正辛烷3种不同溶剂中的构象行为、单嘧磺隆与3种溶剂之间的相互作用能及氢键相互作用进行了计算研究. 计算结果表明, 在3种不同的溶剂中, 单嘧磺隆的优势构象不同; 其构象转换过程, 特别是转换成活性构象的过程主要发生在水溶液中; 与溶剂分子间的相互作用是分子构象行为的决定因素; 单嘧磺隆的脲桥部分可以和含氢键接受体的溶剂形成氢键, 分子间与分子内氢键的竞争可能是从晶体构象转换成活性构象的主要驱动力.     

Morphology development of main‐chain liquid crystalline polymer fibers during solvent evaporation

A nonequilibrium thermodynamic approach has been developed for describing the emergence of fiber morphologies from a liquid crystalline polymer solution undergoing solvent evaporation, including fibrillar structures, concentric rings, and spiral structures. We utilized Matsuyama–Kato free energy for main‐chain liquid crystalline polymer (MCLCP) solutions, which is an extension of Maier–Saupe theory for nematic ordering and incorporates a chain‐stiffening, combined with Flory‐Huggins free energy of mixing. Temporal evolution of the concentration and nematic order parameters pertaining to the above free energy density of liquid crystalline polymer solution was simulated in the context of time‐dependent Ginzburg–Landau theory coupled with the solvent evaporation rate equation under the quasi‐steady state assumption. The emerged morphological patterns are discussed in relation to the phase diagram of the MCLCP solution and the rate of solvent evaporation. © 2007 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 45: 429–435, 2007     

Phase separation phenomena in solutions of poly(2,6-dimethyl-1,4-phenylene oxide). IV. Thermodynamic parameters for solutions in a series of homologous solvents: Toluene to hexylbenzene

Melting-point curves for solutions of poly(2,6-dimethyl-1,4-phenylene oxide) (PPO) in a series of homologous solvents (toluene to n-hexylbenzene) have been obtained from visual and differential scanning calorimetry measurements. The measured melting points were used to calculate thermodynamic interaction parameters. It was found that consistent values were obtained with the Flory–Hoffman melting-point depression equation, if the assumption was made that solvent molecules are incorporated in the crystal lattice. To this end, an adapted dependence of the enthalpy of melting per polymer unit on the cocrystallizing solvent was used. The values of the thermodynamic interaction parameters in the series of solvents and their dependence on polymer weight fraction are explained qualitatively with simplified versions of equation-of-state theory and solubility parameter theory.     

Effect of attractions on the structure of polymer solutions confined between surfaces: a density functional approach

A density functional theory is presented to study the effect of attractions on the structure of polymer solutions confined between surfaces. The polymer molecules have been modeled as a pearl necklace of freely jointed hard spheres and the solvent as hard spheres, both having Yukawa-type attractions and the mixture being confined between attractive Yukawa-type surfaces. The present theory treats the ideal gas free energy functional exactly and uses weighted density approximation for the hard chain and hard sphere contributions to the excess free energy functional. The attractive interactions are calculated using the direct correlation function obtained from the polymer reference interaction site model theory along with the mean spherical approximation closure. The theoretical predictions on the density profiles of the polymer and the solvent molecules are found to agree quite well with the Monte Carlo simulation results for varying densities, chain lengths, wall separations, and different sets of interaction potentials.     

基于链端体积分数的聚合物溶液理论

提出“链端体积分数”的概念以克服Flory-Huggins(FH)理论中的平均场假设.将聚合物多链体系看成一条“间隔连续链”,使一些热力学量写成路径积分的形式.在聚合物/溶剂体系中,保留了FH理论的无热熵假设,把总熵分为两部分,建立了一个新的聚合物溶液理论.计算了聚苯乙烯/环己烷体系的FH相互作用参数和相分离曲线,并与实验数据和FH理论作了比较,结果表明,我们的理论结果比FH的理论结果有了很大的改进.     

Phase equilibria and phase separation dynamics in a polymer composite containing a main-chain liquid crystalline polymer

Various topological phase diagrams of blends of main-chain liquid crystalline polymer (MCLCP) and flexible polymer have been established theoretically in the framework of Matsuyama–Kato theory by combining Flory–Huggins (FH) free energy for isotropic mixing, Maier–Saupe (MS) free energy for nematic ordering in the constituent MCLCP, and free energy pertaining to polymer chain-rigidity. As a scouting study, various phase diagrams of binary flexible polymer blends have been solved self-consistently that reveal a combined lower critical solution temperature (LCST) and upper critical solution temperature (UCST), including an hourglass phase diagram. The calculated phase diagrams exhibit liquidus and solidus lines along with a nematic–isotropic (NI) transition of the constituent MCLCP. Depending on the strengths of the FH interaction parameters and the anisotropic (nematic–nematic) interaction parameters, the self-consistent solution reveals an hourglass type phase diagram overlapping with the NI transition of the constituent MCLCP. Subsequently, thermodynamic parameters estimated from the phase diagrams hitherto established have been employed in the numerical computation to elucidate phase separation dynamics and morphology evolution accompanying thermal-quench induced phase separation of the MCLCP/polymer mixture. © 2006 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 3621-3630, 2006     

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.     

Swelling of polyacrylamide gels in polyacrylamide solutions

Swelling behavior of polyacrylamide (PAAm) and polyacrylamide-co-polyacrylic acid (PAAm-co-PAAc) gels was investigated in aqueous solutions of monodisperse PAAms with molecular weights (M_w) ranging from 1.5 × 10³ to 5 × 10⁶ g/mol. The volume of the gels decreases as the PAAm concentration in the external solution increases. This decrease becomes more pronounced as the molecular weight of PAAm increases. The classical Flory–Huggins (FH) theory correctly predicts the swelling behavior of nonionic PAAm gels in PAAm solutions. The polymer–polymer interaction parameter χ₂₃ was found to decrease as the molecular weight of PAAm increases. The swelling behavior of PAAm-co-PAAc gels in PAAm solutions deviates from the predictions of the FH theory. This is probably due to the change of the ionization degree of AAc units depending on the polymer concentration in the external solution. © 1998 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys 36: 1313–1320, 1998     

The Flory mixing free energy: the effects of solute topology and monomer size

This note presents a lattice‐based derivation of the Flory free energy that does not invoke assumptions on the chain character of the solute. This derivation highlights two crucial features of the Flory free energy: (i) The mixing entropy S_mix is insensitive to the solute structural details. It depends only on the number of lattice sites occupied by the solute. (ii) In marked distinction, the mixing energy does vary with the topology of the solute, the monomer size and the dimensionality. Our discussion suggests re‐examination of the applicability of the Flory free energy to non‐polymeric solutions.     

A molecular thermodynamic model for the swelling of thermo-sensitive hydrogels

A new molecular thermodynamic model for describing the swelling behavior of thermo-sensitive hydrogels was developed. The model consists of two terms. One is the contribution of the mixing of hydrogel network and water, which is dependent on the local polymer concentration and the interaction between polymer segment and solvent. A closed packed lattice model for polymer solution developed by Yang et al. was adopted for this term. The other is the elastic contribution derived from the network elasticity, which is dependent on the cross-linking degree of gel network. The elastic Gibbs energy model based on the Gaussian chain model developed by Flory was adopted. The model equation has two parameters. One is an energy parameter ? reflecting the interaction between water and gel network, the other is a size parameter V^* that represents the cross-linking degree of the hydrogel. When the energy parameter ? is expressed as a quadratic of inverse temperature, this model can describe the swelling equilibrium behavior of neutral thermo-sensitive hydrogels quite well. The influences of model parameters were discussed in details. The experimental swelling curves of two kinds of polyacrylamide-based gels were correlated and good agreement was obtained.     

Opposite trends in thermodynamic compatibility between copolyamide and chitosan in their binary blend

Gibbs energy, enthalpy, and entropy of mixing in binary blends of chitosan with ter‐copolyamide 6/66/610 at ambient conditions have been determined over the entire concentration range using thermodynamic cycle based on dissolution of individual polymers and their blends of different composition in a common solvent – formic acid. Experimental procedure included stepwise equilibrium vapor sorption of glacial formic acid on the cast films and isothermal microcalorimetry of dissolution of these films in liquid glacial formic acid at 25 °C. Formic acid appeared to be a very good solvent for individual polymers and their blends. Flory‐Huggins interaction parameter determined from sorption isotherms was negative and varied from ?2.56 to ?1.79 depending upon blend composition. The enthalpies of dissolution of individual polymers and their blends were strongly exothermic and varied from ?200 to ?40 Joule/g. Independent thermodynamic cycles for Gibbs free energy and enthalpy remarkably revealed similar trends in concentration dependence of different thermodynamic functions of mixing between chitosan and copolyamide. At high chitosan content, the binary blend is characterized by large and negative values of Gibbs free energy, enthalpy, and entropy of mixing that provide high polymer compatibility. On the contrary, at high copolyamide content the blends are incompatible and are characterized by positive values of enthalpy, entropy, and Gibbs free energy of mixing. Such complicated thermodynamic behavior is the result of the superposition of strong molecular interactions (H‐bonds) between polymers in the blend and isothermal fusion of copolyamide crystallites. Thermodynamic analysis correlates well with the data obtained by polarized microscopy. © 2007 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 45: 2603–2613, 2007