端基附壁模型聚合物环形链的构象统计理论

环形链定义为两个端基均附壁的线形聚合物链所得的环 .采用精确计数和蒙特卡罗模拟方法 ,研究了自回避行走 (SAW)模型表示的环形链的构象 ,求得构象数和均方回转半径随链长的变化 .对于二维和三维SAW环形链 ,精确计数的最大链长分别为N =2 9和 1 9.用标度理论处理了数值结果 ,所得的标度指数和其他参数与理论预示值进行了比较 .模拟结果表明 ,SAW环形链限制壁平行方向的尺寸大于垂直方向的尺寸 ,与SAW尾形链尺寸的变化正好相反 .     

端基附壁高分子链的构象统计理论——Ⅱ.SAW尾形链

研究了SAW尾形链,采用精确计数和Monte Carlo模拟方法求得d维(d=2,3)SAW尾形链的构象数C_1~d和均方末端距[h_1~d(N)]~2及其分量随链长N的变化关系.发现它们与自由SAW链一样都服从标度律,从模拟数据拟合求出了这些SAW尾形链的临界指数.还表明短链SAW在壁的法向与NRW尾形链一样有所伸展,均方末端距的法向分量几乎是平行分量的2倍;但随N→∞,链自回避效应对壁的作用有所屏蔽,这与NRW尾形链有明显不同.     

嵌段高分子尾形链构象性质的Monte Carlo研究

基于简立方格点模型对AB两嵌段高分子尾形链的构象性质及其链节的空间分布进行了MonteCarlo模拟.结果表明,链的尺寸、形状和链节的空间分布等统计性质和B链节与平面壁之间的相互吸引能有关.随着B链节与壁之间的吸引能的增加,链的尺寸和形状均呈现出先下降后升高的变化趋势,而且B链节的比例越大,这种变化越明显.     

金刚石格点上自避行走尾形链的构象统计理论

研究了金刚石格点上自避随机行走（ＳＡＷ）尾形链，采用精确计数和ＭｏｎｔｅＣａｒｌｏ模拟方法求得该ＳＡＷ尾形链的构象数Ｃ（Ｄ）１（Ｎ）和均方末端距［ｈ（Ｄ）１（Ｎ）］２及其分量随链长Ｎ的变化关系．发现它们与自由ＳＡＷ链一样都服从标度律，从这些量的计算机实验数据拟合求出了金刚石格点上ＳＡＷ尾形链的临界指数和格点指数．计算结果还表明短链ＳＡＷ在壁的法向与ＮＲＷ尾形链一样有所伸展，均方末端距的法向分量几乎是平行分量的２倍；但随Ｎ→∞，链自回避效应对壁的作用有所屏蔽．这些都与简立方格子模型上得到的结果一致．     

聚乙烯型尾形链构象统计的Monte Carlo研究

以聚乙烯型尾形链为例 ,根据链的实际键长、键角和Θ条件下每一链段处于不同旋转异构态的条件概率 ,用MonteCarlo模拟方法生成样本链分子 ,计算了尾形链的均方末端矩及其分量、构象数 ,并与完全计算法及理论推导值作了比较 ,讨论了尾形链构象数和沿边界方向及法线方向均方末端矩分量随链段数n及键角Φ的变化 .     

近壁模型高分子链的计算机模拟

采用简立方格点上的普通无规行走 (NRW)和自避行走 (SAW )为模型链 ,研究了不可穿透的刚性壁附近高分子的构象性质 ,得出近壁链的构象熵和尺寸随一端点与壁之间距离z0 的变化 .所用的计算机模拟技术包括精确计数和MonteCarlo模拟 .数值结果表明 ,近壁链的构象熵降低遵循一种简单规律 ,即当链长愈长或一端愈接近于壁时 ,与长度相同的自由链相比 ,链因壁限制所致构象熵的降低愈大 .当链十分接近于壁时 ,其均方末端距〈R2 〉(尤其是与壁垂直的分量〈R2z〉)大于自由链的相应值 ;随z0 增大 ,〈R2 〉及〈R2z〉开始减小 ,通过某一极小值 ,然后上升 ;当z0 →∞时 ,趋于自由链的极限值 .换言之 ,近壁模型链的线团经历一个收缩再逐渐扩张的过程     

磁性高分子链的磁性质和构象性质的Monte Carlo模拟

采用退火 (Annealing)MonteCarlo方法 ,从高温到低温顺序模拟了简立方格点上考虑最近邻Ising相互作用的磁性高分子链在不同温度的磁性质和构象性质 .磁性高分子链在低温下存在自发磁矩 ,无限长链的临界温度Tc=1 77± 0 0 5J kB.在临界温度附近 ,高分子链经历了从伸展的无规线团到紧缩球体的塌缩相变 .对链的尺寸、形状、近邻数及能量的分析表明 ,高分子链的构象性质从温度Tc=1 77开始发生较明显的变化 ,这表明高分子Ising链的相变是Ising相互作用和链节运动协同作用的结果 .     

端基附壁高分子链的构象统计理论——Ⅲ.NRW环形链

当线形高分子链两端被固体壁吸附后,形成了“环形链”。采用半空间格点中的普通随机行走(NRW)模型讨论了环形链构象的统计理论。在NRW尾形链构象分布函数的基础上,容易推导出NRW环形链构象数的公式。当链长N→∞时,环形链构象数与自由链构象数之比随N~(-3/2)而变化。用递归方程也导出了同样的结果。环形链的均方末端距(?)与均方键长l~2之比为2N/3,与长度相同的自由链比较,有一定收缩。所导出的基本关系得到了精确计数和Monte Carlo模拟的支持。     

环形高分子单链塌缩的熔球中间态及其热力学解释

采用动态蒙特卡洛分子模拟研究了环形高分子单链在不良溶剂中发生塌缩转变时可逆地出现具有核-壳结构特征的熔球中间态,发现该结构特征与相同链长的线形单链基本相同,表明其只与链的长短有关,而与链端基的特殊效应无关.本工作将这一现象与单链单晶在其平衡熔点附近出现的类似现象相互关联,采用表面预溶模型来解释单链塌缩出现熔球中间态的热力学机理.分子量越低,熔球越小,表面预溶现象就越显著,塌缩转变随热力学条件变化就越缓慢.实际的高分子体系由于链内拓扑缠结,在表面未必能充分释放片段链,达不到理论预期的平衡态.表面预溶使得相分离临界点或晶体熔点附近在界面厚度方向上存在链单元能量状态不连续分布,这在微观分子水平上与临界界面连续浓度梯度的传统理论处理不一致,为我们深入理解高分子流体界面的微观结构带来帮助.     

高分子链构象性质和溶解性的Monte Carlo模拟

运用Monte Carlo方法对线型高分子链格点模型的构象进行了模拟,研究了单链体系的构象的尺寸（采用均方末端距,均方回转半径来表征）、形状（采用平均非球形因子来表征）和溶解性随溶剂与链段间相互作用能的变化情况.结果表明,、和随着εPS的增大而减小,具有相同的变化规律;随着Δε增大,溶剂变得越来越不利于溶解,高分子链的形状蜷曲为椭球形,高分子线团相互穿插交叠.     

Conformational statistics of polymer chain terminally attached to wall (II)

The SAW tail chains were studied. The permitted conformational number and the mean square end-to-end distance as a function of the chain length N for such a model tail chain were obtained by computer simulations, including the exact enumeration and Monte Carlo method. These two basic quantities obeyed the relations deduced from the scaling law. The critical exponents and the lattice indexes were given by fitting the data of the computer experiments. It has been shown that there is a certain extension in the size of the SAW tail chains as well as the NRW tail chains in the direction normal to the wall. The normal component of the mean square end-to-end distance is almost twice as large as the parallel component of the short chain SAW. However, as N →∞, the effect of the wall on the chain conformation becomes a little weak because of the self-avoiding behavior for the model. That is quite different from the case of the NRW tail chain. Project supported by the National Natural Science Foundation of China     

Conformational statistics of polymer chain terminally attached to wall (III) —NRW model loop chain

When the two end groups of a linear polymer chain are absorbed on a solid surface, the polymer chain forms the “loop” conformation. Investigation has been made on the conformational statistics of a model loop chain by the normal random walk (NRW) on a lattice confined in the half-infinite space. Based on the conformational distribution function of the NRW model tail chain, it is easy to deduce an analytical formula expressing the conformational number of the model loop chain. It was found that the ratio of the conformational number of the model loop chain to that of the free chain varies with the power functionN ^-2/3 when the chain lengthN→οο. The same result -was obtained by means of the recursion equation. The ratio of the mean square end-to-end distanceh ² for the model loop chain to its mean square bond lengthl ² is 2N/3. Compared with the free chain with the same lengthN, the mean square end-to-end distance of the model loop chain contracts to a certain extent. The basic relationships deduced were supported by the exact enumeration and Monte Carlo simulations. Project supported by the National Natural Science Foundation of China.     

The directed self-assembly of amphiphilic diblock copolymers in selective solvents

The directed self-assembly of diblock copolymers in solvents is studied systematically using a simulated annealing method. Effects of the shape, scale, and adsorption capacity of the induced surface on the morphology of the aggregates are examined. A variety of morphologies are predicted. By increasing the scale of induced surface, the micellar shape transforms from cylinder to sheet with a tail and finally to thin sheet without tail. The shape of induced surface determines the sheet’s shape, such as rounded and square. Configurations of hydrophobic blocks and interfacial energies are investigated by calculating the mean square end-to-end distances and the contact numbers between hydrophobic monomer and other species, respectively.     

模型高分子尾型链的构象统计理论

基于金刚石格点上受限于半空间的随机行走，借概率论中的反射原理推导出了模型高分子尾型链的构象分布．得到了链长为Ｎ的模型尾型链的允许构象数Ｃ１（Ｎ）和末端距分布函数．在四选择模型中，发现Ｃ１（Ｎ）／４Ｎ与Ｎ－１／２成正比．模型尾型链在与壁正交方向上的均方末端距分布与平行方向比较扩张至２倍，而后者与自由链对应的分量相等．这些解析结果得到了精确计数和ＭｏｎｔｅＣａｒｌｏ模拟计算机“实验”的支持．     

Molecular Weight Effect on Latex Film Formation from Polystyrene Particles: A Photon Transmission Study

The photon transmission method was used to probe the time evolution of film formation from latex particles. Two different latex films were prepared from high molecular weight (HM) and low molecular weight (LM) polystyrene particles at room temperature and were annealed at various temperatures in 2.5-min time intervals above the glass transition. The increase in the transmitted photon intensity (Itr) is attributed to the increase in "crossing density" at the junction surface. The Prager-Tirrell model was employed to interpret the increase in crossing density at the junction surface. The back and forth activation energies were measured for HM and LM films and found to be around 59 and 87 kcal/mol for a reptating polymer chain across the junction surface. Monte Carlo simulations were performed for photon transmission through a rectangular lattice. The number of transmitted photons (Ntr) was calculated as a function of the mean free path of photons. It was observed that Ntr, similar to Itr, increases as the square of the mean free path of photons is increased. Copyright 1999 Academic Press.     

拉力诱导下的紧密高分子链相变行为

采用相互作用自回避行走(interacting self-avoiding walks,ISAWS)模型研究了一端固定的紧密高分子链在拉伸过程中的低温相变行为,观察到在拉伸过程中当温度T<0.1时平均拉力会出现一个震荡,随着温度的升高这种震荡现象又渐渐消失,这是由于紧密高分子链在低温时类似于β折叠的"冻结构象"被拉开而引起的.比较吸附条件下和无吸附作用下平均拉力、自由能以及相变行为的差别,发现在吸附条件下在拉伸的初始阶段为了克服表面吸附的相互作用,拉力会出现一个峰.吸附作用也使得外界作用到高分子链上的实际有效拉力减小,造成崩塌相态(collapsed phase)区域面积减少.另外发现在吸附条件下平均拉力还受温度变化的影响.在拉伸的初期由于单体间存在体积排除效应,平均拉力是随着温度的升高而降低,随着拉伸的深入当末端距到达一定长度时平均拉力是随着温度的升高而增加.并同Kumar等人在不考虑吸附作用下拉伸紧密高分子链得到的结果进行了比较.这些研究对于进一步研究外力诱导下吸附紧密高分子的相变有一定的参考价值.     

COMPUTER SIMULATION OF A SINGLE POLYMER CHAIN IN DIFFERENT SOLVENTS

In the present paper, the behavior of a single polymer chain under various solvent conditions was modeled by self-avoiding walks (SAW) with nearest neighbors attraction Δε on a simple cubic lattice. Determination of the θ-condition wasbased on the numerical results of the mean square radius of gyration and end-to-end distance. It was found that at the θtemperatue Δε/kT equals -0.27. The exponents a in the Mark-Houwink equation with different interaction parameters areconsistent with the results of experiments: under θ-condition, a = 0.5, and for a good solvent α= 0.74-0.84, respectively.