Full‐chain dynamics of entangled linear and star polymers

The full‐chain dynamics and the linear viscoelastic properties of monodisperse, entangled linear and star polymers are simulated consistently via an equilibrium stochastic algorithm, based on a recently proposed full‐chain reptation theory 1 that is able to treat self‐consistently mechanisms of chain reptation, chain‐length fluctuations, and constraint release. In particular, it is the first time that the full‐chain simulation for star polymers is performed without subjecting to the great simplifications usually made. To facilitate the study on linear viscoelasticity, we employ a constraint release mechanism that resembles the idea of tube dilation, in contrast to the one used earlier in simulating flows, where constraint release was performed in a fashion similar to double reptation. Predictions of the simulation are compared qualitatively and quantitatively with experiments, and excellent agreement is found for all investigated properties, which include the scaling laws for the zero‐shear‐rate viscosity and the steady‐state compliance as well as the stress relaxation and dynamic moduli, for both polymer systems. The simulation for linear polymers indicates that the full‐chain reptation theory considered is able to predict very well the rheology of monodisperse linear polymers under both linear viscoelastic and flow conditions. The simulation for star polymers, on the other hand, strongly implies that double reptation alone is insufficient, and other unexplored mechanisms that may further enhance stress relaxation of the tube segments near the star center seem crucial, in explaining the linear viscoelasticity of star polymers. © 2000 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys 38: 248–261, 2000     

Thermodynamic Study of the Molten Salt Binary System KHSO4−NaHSO4

The partial molar enthalpies of mixing of NaHSO₄ and KHSO₄ have been measured at 528 K by dropping samples of pure compounds into molten mixtures of NaHSO₄ and KHSO₄ in Calvet calorimeter. From these values the molar enthalpy of mixing has been deduced. The same method has been used for the determination of the heat capacity of the two pure compounds in the solid and liquid states. The phase diagram of this system has been confirmed by conductometric and thermal analysis methods. By an optimization method the excess entropy of the liquid mixtures was also calculated. This revised version was published online in August 2006 with corrections to the Cover Date.     

Interdiffusion of polymers across interfaces

Neutron Reflection (NR) and Dynamic Secondary Ion Mass Spectroscopy (DSIMS) experiments were conducted on symmetrically deuterated polystyrene triblock bilayers (HDH/DHD) which directly probed the interdiffusion dynamics of the chains during welding. The HDH chains had their centers deuterated 50%, the DHD chains had their ends deuterated (25% at each end) such that each chain contained approximately 50% D. During welding, anisotropic motion of the chains produces a time-dependent oscillation (ripple) in the H and D concentration at the interface, which bears the characteristic signature of the polymer dynamics. These oscillations were compared with those predicted by Rouse, polymer mode coupling (PMC), and reptation dynamics. The following conclusions can be made from this study. (a) During the interdiffusion of high molecular weight HDH/DHD pairs, higher mobility of the chain ends caused a concentration oscillation which increased to a maximum amplitude, and eventually vanished at times, t > τ_D. The amplitude, or excess enrichment found, was appreciably more than that predicted by Rouse and PMC simulations, and was only slightly less than that predicted from reptation simulations. (b) The oscillations were completely missing in the 30 and 50K HDH/DHD polymers, which are only weakly entangled. The lack of oscillations for the 30 and 50K pairs may be due to a combination of surface roughness and fluctuations of order 30 Å. (c) It was found that the position of the maximum in this ripple stayed at the interface during its growth. This is also consistent with reptation and has not been explained by other theories. (d) All dynamics models for linear polymers produce ripples, many of which are qualitatively similar to that predicted for reptation. However, each ripple bears the fingerprint of the dynamics in terms of its time-dependent shape, position, and magnitude, and the models are clearly distinguishable. Our results, in summary, support reptation as a candidate mechanism of interdiffusion at polymer(SINGLEBOND) polymer interfaces and its uniqueness is being further pursued. © 1996 John Wiley & Sons, Inc.     

柱阵列中高分子的爬行过程

利用Monte Carlo算法模拟了嵌入柱阵列中高分子的爬行过程.模拟结果显示窄间距的柱阵列可以有效地抑制嵌于其中的高分子链的长度涨落,并使得高分子链的弛豫时间符合爬行理论给出的标度律.通过对比爬行理论与模拟得到的链端关联函数可知,柱阵列中的爬行过程与高分子熔体中的爬行过程存在着微妙的差别.在柱阵列中高分子的链端无法完全消灭原有的管道,从而导致了链端关联函数被爬行理论低估.本文给出的二维自回避格点模型则可以很好地描述柱阵列中高分子的爬行过程.     

分子印迹手性分离过程的热力学研究

以丙烯酰胺为功能单体, 以二甲基丙烯酸乙二醇酯为交联剂, 在模板分子N-叔丁氧羰酰-L-色氨酸(N-Boc-L-Trp)和N-叔丁氧羰酰-L-酪氨酸(N-Boc-L-Tyr)的存在下, 分别采用光引发聚合和热引发聚合制备了N-Boc-L-Trp和N-Boc-L-Tyr的分子印迹聚合物(MIPs), 进行分子印迹手性分离过程的热力学研究. 测定了分离过程的熵变、焓变和自由能变化. 结果显示, 在流动相中添加异丙醇或甲醇等强氢键竞争性溶剂时, 熵变对分离起到了主要作用, 而且分离过程中的溶剂化对分离的影响也非常大. 分子印迹聚合物对印迹分子和非印迹分子进行分子识别的主要作用是印迹聚合物与印迹分子匹配的三维空间结构.     

Welding dynamics in an atomistic model of an amorphous polymer blend with polymer–polymer interface

We consider an atomistic model of thermal welding at the polymer-polymer interface of a polyetherimide/polycarbonate blend, motivated by applications to 3D manufacturing in space. We follow diffusion of semiflexible chains at the interface and analyze strengthening of the samples as a function of the welding time t_w by simulating the strain–stress and shear viscosity curves. The time scales for initial wetting, and for fast and slow diffusion, are revealed. It is shown that each component of the polymer blend has its own characteristic time of slow diffusion at the interface. Analysis of strain–stress demonstrates saturation of the Young's modulus at t_w = 240 ns, while the tensile strength continues to increase. The shear viscosity is found to have a very weak dependence on the welding time for t_w > 60 ns. It is shown that both strain–stress and shear viscosity curves agree with experimental data.     

Thermodynamic properties of novolac-type phenolic resin blended with poly(ethylene oxide)

The thermodynamic properties of novolac type phenolic resin blended with poly(ethylene oxide) (PEO) were investigated by the Painter–Coleman association model (PCAM). Equilibrium constants and enthalpy corresponding to the interaction between phenolic and poly(ethylene oxide) were calculated from the Fourier transform infrared spectroscopy of low molecular weight analogues in dilute solutions. The association parameters of the model compounds are transferred to the corresponding polymers, to predict the Gibbs free energy, phase behavior, and the degree of hydrogen bonding in the polymer blend. The heat capacity (C_P) and the excess heat capacity (ΔC_P) are used to verify the validity of PCAM model on predicting the thermodynamics properties of phenolic/PEO blend. It is found that the hydrogen bonding interaction dominates at moderate temperatures, which is outweighed by the dispersion force at higher temperature or high PEO compositions. © 1998 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys 36: 1647–1655, 1998     

无盐聚电解质溶液的分子热力学模型

聚电解质在水处理和胶体领域应用广泛 ,在生物系统中 ,蛋白质和 DNA等均是聚电解质 .文献中已发表了不少聚电解质溶液的理论和模型 ,如基于柱型胞腔模型的 Poison- Boltzmann( PB)方程理论[1]及其改进型 ( MPB) [2 ] 、Manning反离子凝聚理论 [3 ] 、积分方程理论 [4] 及适合于工程应用的半经验模型 [5 ]等 ,但离工程实际应用的要求还相差较远 .最近 ,姜建文等 [6]建立的分子热力学模型 ,其预测结果与 MD模拟数据非常吻合 .本文在此基础上进一步考虑聚离子与反离子间的缔合作用 ,可以在中高浓度范围内很好地关联实验数据 .1　模　　型…     

精馏节能过程非平衡热力学分析——模型方程的建立

通过对精馏过程进行热力学分析,建立了通用性较强的非平衡热力学数学模型.该模型可对精馏过程中用能情况进行非平衡热力学分析和熵均分分析,也可用于评价塔板性能.模型建立的基础只适用于线性非平衡热力学区,即近平衡区,有较大的局限性.     

Predicting the proton conductivity of perfluorosulfonic acid membrane via combining statistical thermodynamics and molecular dynamics simulation

The electrochemical properties of a perfluorosulfonic acid (PFSA) membrane are estimated using a combination of molecular dynamics simulation and statistical thermodynamic model. We obtain all parameters in an ionic conductivity model from an atomistic simulation and remove all adjusted model parameters. From a microscopic point of view, the hydrated PFSA membrane shows micro‐phase segregation which separated into hydrophilic and hydrophobic phases. Our present work originates with this phenomenon and we treat this phase segregation as if it is a continuous phase for each of which the proton (H⁺) is transported inside the PFSA membrane/solvent (water and alcohols) mixture. The chemical potential for a given system is estimated using a molecular simulation technique to predict the van der Waals interaction energy between the polymer and solvent. In addition, the self diffusion coefficients are calculated from the molecular dynamics simulation. We study various polymer/solvent compositions to understand the concentration dependence of self diffusion coefficient. Our self diffusion coefficients and also the predicted final ionic conductivity agree well with previously reported experimental data. © 2011 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 49: 1455–1463, 2011     

卤代硅烯异构化为卤代甲硅烷基硅烯的热力学与动力学分析

利用密度泛函理论(DFT) 中的B3LYP方法,全参数优化了卤代硅烯HXSi=SiXH（X=F、Cl、Br、I）异构化反应的反应物、产物及过渡态的几何构型,计算出了它们的振动频率、零点振动能(ZPVE)和总能量,并对它们进行了振动分析,以确定过渡态的真实性。又计算了反应的热力学函数变化,平衡常数及速率常数,比较了不同卤素对反应的影响。热力学与动力学计算结果表明,该异构化反应过程是一个放热的、且在低温下可自发进行的反应,但对于溴代硅烯和碘代硅烯而言,当温度达到1000 K时,反应开始转化为非自发反应。     

Structure,Dynamics and Properties of Materials with Polymers Having Complex Architectures

Various non-linear highly branched polymers such as multiarm stars, block copolymer micelles and bottlebrush-like polymers have been studied in order to analyze their intramolecular structure and effects of spatial ordering resulting from their specific macromolecular architecture. These polymers constitute a class of compact macromolecules which, due to the high intramolecular density, interact strongly, excluding each other in space. Investigations of the structure and dynamics in such systems, using various experimental methods as well as computer simulations, have been performed. Small angle X-ray scattering is used to characterize the structure and mechanical spectroscopy is used for the detection of the dynamic behavior of the systems. Simulations have been performed using the cooperative motion algorithm with lattice polymers equivalent to the considered macromolecules. Results of both experiments and the simulation seem to support the concept of slow structural cooperative rearrangements controlling the flow in such systems. The effects are analogous to those related to flow in low molecular liquids but take place on another size scale. The new slow relaxation processes creates new super soft-states which are characterized by shear modulus plateau lower than 10⁴ Pa.     

从统计热力学角度理解重要的热力学概念

基于玻尔兹曼统计法,简要总结独立粒子的非定位系统的统计热力学基础,并应用一些实例来分析统计热力学和热力学方法之间的区别和联系,帮助学生从统计热力学角度加深对一些重要的热力学概念(如热容、吉布斯自由能、焓和化学势)的理解,建立微观物理图像,以更好地学习两部分热力学知识,并能够灵活运用于化学问题的求解中。     

A Unified Thermodynamic Model of Flow-induced Crystallization of Polymer

We propose a unified thermodynamic model of flow-induced crystallization of polymer (uFIC),which incorporates not only the conformational entropy reduction but also the contributions of flow-induced chain orientation,the interaction of ordered segments,and the free energy of crystal nucleus and crystal morphology.Specifically,it clarifies the determining parameters of the critical crystal nucleus size,and is able to account for the acceleration of nucleation,the emergence of precursor,different crystal morphologies and structures induced by flow.Based on the nucleation barrier under flow,we analyze at which condition precursor may occur and how flow affects the competition among different crystal forms such as orthorhombic and hexagonal phases of polyethylene.According to the uFIC model,the different crystal morphologies and structures in the flow-temperature space have been clarified,which give a good agreement with experiments of FIC.     

Frequency and molecular-mass-dependent nonexponential proton NMR relaxation in polymer melts

A theoretical treatment of the nonexponential relaxation behavior of the different proton nuclear magnetic resonance (NMR) relaxation processes in polymer melts is presented. Formulas are derived for a three-component model given by two versions and a homogeneous distribution of correlation times. The theoretical results were tested with measurements of T₁, T_2e, and T₂ as functions of frequency and molecular mass in linear fractionated polyethylene samples. While the T₁ relaxation always yields exponential magnetization decays, the T_2e and T₂ measurements show biexponential relaxation behavior. From the calculations it was found that the correlation time of the local motion is independent of the molecular mass, whereas the correlation time of the slowest motional process increases with M^2.8_w for the three-component model and with M^2.2_w for the distribution of correlation times, respectively. © 1992 John Wiley & Sons, Inc.     

Analysis of Thermodynamic Parameters of Glass Forming Polymeric Melts

The temperature dependence of the Gibbs free energy difference (ΔG), enthalpy difference (ΔH) and entropy difference (ΔS) between the undercooled meltand the corresponding equilibrium solid has been analysed for glass forming polymeric materials by calculating ΔG, ΔH and ΔS within the framework of the hole theory of liquids. The study is made for nine samples of glass forming polymeric melts; polypropylene oxide (PPO), polyamid-6 (PA-6), polytetramethylene oxide (PTMO), polyethylene oxide (PEO), polystyrene (PS), polypropylene (PP), polyethylene (PE), polyethylene terephthalate (PET) and polybutadiene (PB) and three simple organic liquids: tri-α-naphthyl benzene (tri-α-NB), o-terphenyl (o-ter) and phenyl salicylate (salol) in the entire temperature range T _m (melting temperature) to T _g (glass transition temperature). The ideal glass transition temperature (T _K) and the residual entropy (ΔS _R) of these samples have also been studied due to their important role in the study of the glass forming ability of materials. This revised version was published online in August 2006 with corrections to the Cover Date.     

Disentanglement and reptation during dissolution of rubbery polymers

The dissolution mechanism of rubbery polymers was analyzed by dividing the penetrant concentration field into three regimes that delineate three distinctly different transport processes. The solvent penetration into the rubbery polymer was assumed to be Fickian. The mode of mobility of the polymer chains was shown to undergo a change at a critical penetrant concentration expressed as a change in the diffusion coefficient of the polymer. It was assumed that beyond the critical penetrant concentration, reptation was the dominant mode of diffusion. Molecular arguments were invoked to derive expressions for the radius of gyration, the plateau modulus, and the reptation time, thus leading to an expression for the reptation diffusivity. The disentanglement rate was defined as the ratio between the radius of gyration of the polymer and the reptation time. Transport in the second penetrant concentration regime was modeled to occur in a diffusion boundary layer adjacent to the polymer-solvent interface, where a Smoluchowski type diffusion equation was obtained. The model equations were numerically solved using a fully implicit finite difference technique. The results of the simulation were analyzed to ascertain the effect of the polymer molecular weight and its diffusivity on the dissolution process. The results show that the dissolution can be either disentanglement or diffusion controlled depending on the polymer molecular weight and the thickness of the diffusion boundary layer. © 1996 John Wiley & Sons, Inc.     

A steered molecular dynamics simulation on the elastic behavior of adsorbed star polymer chains

<正>Elastic behavior of 4-branched star polymer chain with different chain length N adsorbed on attractive surface is investigated using steered molecular dynamics(SMD) simulation method based on the united-atom(UA) model for branched alkanes.The simulation is realized by pulling up the chain via a linear spring with a constant velocity v = 0.005 nm/ps.At the beginning,the chain lies extensionally on adsorbed surface and suffers continuous deformations during the tensile process.Statistical parameters as mean-square radii of gyration S~2_(xy),S~2_z,shape factor δ,describing the conformational changes,sectional density den which gives the states of the chain,and average surface attractive energy U_a,average total energy U,average force f probed by the spring,which characterize the thermodynamic properties, are calculated in the stimulant process.Remarkably,distinguishing from the case in linear chains that there only exists one long plateau in the curve of f,the force plateau in our study for star chains is multiple,denoting different steps of desorption,and this agrees well with the experimental results in essence.We find during the tensile process,there are three characteristic distances Z_c,Z_t and Z_0 from the attractive surface,and these values vary with N.When Z=Z_c,the chain is stripped from the surface,but due to the form of wall-monomer interaction,the surface retains weak influence on the chain till Z = Z_c.From Z=Z_t,parameters U_a,U and f respectively reach a stable value,while the shape and the size of the chain still need adjustments after Z_t till Z_0 to reach their equilibrium states.Specifically,for short chain of N= 41,Z_t and Z_0 are incorporated.These results may help us to deepen the knowledge about the elastic behavior of adsorbed star polymer chains.