毛细管内高分子溶液的浓度分布

采用连续自洽场理论分析了毛细管中发生凝胶化之前的聚合物溶液浓度分布的影响因素及其规律. 结果表明, 体系尺寸有限时, 改变聚合物链段、溶剂与壁面的相互作用参数之差, 聚合物溶液浓度分布会发生贫化/吸附转变; 临界作用参数与聚合物链长的倒数呈线性关系, 且拟合常数与体系尺寸、聚合物溶液平均体积分数有关; 聚合物分子量分布为多分散时, 分子量较低的组分更容易接近容器壁面, 分子量较大的组分则相反. 总之, 增加聚合物溶液浓度、链长, 选择优良溶剂, 减小体系尺寸等都会使浓度分布更加均匀.     

表面活性剂与高分子链混合体系的模拟

计算机模拟了高分子链对表面活性剂胶束形成过程的影响，以及高分子链构象性质随胶束化过程的变化.结果表明,当高分子链与表面活性剂之间的相互作用强度超过临界值后，高分子链的存在有利于表面活性剂胶束的形成.临界聚集浓度（CAC）与临界胶束浓度（CMC）的比值CAC/CMC随高分子链长的增大和相互吸引作用的增强而减小.在CAC之前，高分子链与表面活性剂分子只有动态的聚集；但在CAC之后，表面活性剂胶束随表面活性剂浓度X的增加而增大，并静态地吸附在高分子链上，形成表面活性剂/高分子聚集体.随着表面活性剂分子的加入，高分子链的均方末端距和平均非球形因子先保持恒定；从X略小于CAC开始， 和快速减小，至极小值后又逐渐增大.模拟结果支持高分子链包裹在胶束表面的实验模型.     

平行壁面间平衡聚合物吸附行为的自洽场分析

采用自洽场理论研究了平衡聚合物在平行壁面间的吸附规律.结果表明,平衡聚合物的吸附行为可根据壁面吸附强度划分为弱吸附和强吸附两个区间.在强吸附区间,平衡聚合物可以形成明确的吸附层,增大吸附强度能够引起平均分子链长的陡增,但不会改变链长分布的指数形式.平衡聚合物的长链分子在强吸附条件下较短链分子更靠近吸附壁面,在弱吸附条件下,则更接近体系中心.通过计算壁面压强发现:在弱吸附作用下,吸附壁面始终受到分子链的推斥作用;在强吸附作用下,分子链对壁面的作用随壁面间距的加大由推斥转变为吸引.平衡聚合物的溶解性和排空效应以及壁面吸附作用之间存在竞争关系,加大排除体积作用参数会引起分子链密度分布的均一化和壁面承受压强的降低.     

受限状态聚合物熔体的分子动力学模拟

用粗粒化的分子动力学(MD)模拟方法从分子层次研究了受限于粗糙壁内的聚合物熔体的动力学性质. 结果表明, 对于链长较短的受限聚合物熔体体系, 随着膜厚的增加, 体系内部高分子链的松弛时间逐渐减少; 然而对于链长较长的受限体系, 聚合物链的松弛时间随着膜厚的增加先减少后增加. 推测这种由于链长的变化所引起的动力学性质的差异源自受限熔体内聚合物链聚集状态的改变, 并且通过考察交叠参数对这种改变进行了分析. 结果表明, 在膜厚增加的过程中, 决定受限状态高分子长链松弛机理的因素逐渐从受限效应转变成为链间的缠结效应.     

规则星形高分子之间的相互作用势的自洽场研究

用自洽场理论方法计算了星形高分子之间的相互作用势, 并与标度理论的推测以及基于粒子的分子动力学模拟方法的结果进行了比较. 自洽场理论计算得到的星形高分子之间的二体相互作用势与标度理论推测的结果基本一致: 当两个星形高分子核心之间的距离小于其对应的等效软球的半径R时, 二体相互作用势U与星形高分子核心之间的距离d及官能度f之间的关系为U(d)∝-f 3/2 ln(d/R). 当星形高分子核心间距离较小, 相互作用势与星形高分子核心间距离的对数呈线性关系时, 三体相互作用为吸引作用, 但其强度约等于二体作用的10%. 因此, 以二体作用为基础的基于粒子的模拟方法是基本合理的.     

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

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

高分子链在微通道剪切流中的迁移机制及浓度的影响

提出了剪切流中高分子链在微通道内的迁移机制.该机制采用珠-簧链模型表示高分子链,高分子链受剪切作用而被拉伸,相邻珠子之间的流体力学相互作用产生了对称的扰动流场,由于在通道壁面附近对称的流场被破坏,壁面与高分子链间的流体力学相互作用使高分子远离壁面,在强受限时,这种壁面诱导的流体力学相互作用会被屏蔽掉.利用耗散粒子动力学数值模拟了高分子链在微通道压力流中的迁移行为.数值模拟结果表明,在受限较弱时,高分子链向远离壁面的方向迁移,并随着流场增强,远离壁面的趋势越强;在受限较强时,高分子链不会发生远离壁面的行为.实验研究了长链高分子λ-DNA在壁面附近的迁移行为,实验结果及模拟结果与迁移机制预测的结果相吻合,验证了迁移机制的正确性.高分子链浓度会影响高分子链的迁移行为,当高分子链浓度较大时,高分子链在通道宽度方向不会发生迁移现象,意味着随着浓度的增大,壁面与高分子链间的流体力学相互作用会逐渐被屏蔽.     

两亲性嵌段共聚物在选择性溶液中吸附行为的理论研究

通过自洽场理论方法研究了在圆形孔道以及水平基板受限情况下两亲性嵌段共聚物在选择性稀溶液中的自组装形态,考察受限管壁、管壁极性以及基板间距对聚合物吸附行为的影响.中性圆形孔道中,随着溶剂对疏水链的排斥作用增强和对亲水段的亲水作用增大,聚合物以蘑菇形状吸附于管壁上.疏水链、亲水链与溶剂的相互作用差异性越大,越容易产生吸附现象.在亲水性孔道中,吸附现象有所缓解,随着管壁对亲水段吸引作用的增强,聚合物在管壁附近形成规则排列的球状胶束.当疏水链、亲水链与溶剂的相互作用相差非常大时,即使管壁对亲水嵌段的吸附作用增强,也不能消除聚合物吸附现象.在中性水平基板中,随着基板间距的增加,两亲性嵌段共聚物在基板附近依次出现球状胶束-蘑菇状胶束-单层球状胶束-高分子刷-双层球状胶束-对称高分子刷的自组装形态.     

聚合物链长及表面吸附强度对三维受限下嵌段共聚物微观相分离的影响

利用动态密度泛函(Dynamic density functional theory, DDFT)方法研究了三维受限下嵌段共聚物的微观相分离, 讨论了共聚物链长和表面吸附强度对微观相形成与取向的影响. 体系中随机分布的等径微球提供三维限制结构, 体积分数为0.6. 增加微球的半径和体积分数, 能够使其从破坏微相规整结构的纳米掺杂过渡到提供三维限制结构. 调整嵌段共聚物与微球表面的相互作用对微相形成与取向有重要影响.     

弱电场驱动下高分子链在周期管道内输运的模拟研究

采用Langevin动力学方法模拟研究了弱电场驱动下高分子链在无限长周期管道中的输运过程. 管道由长度相等的α和β两部分周期排列而成, 其中高分子链与α管道间存在相互吸引作用, 而与β管道间存在纯排斥作用. 模拟结果表明, 高分子链在输运过程中存在明显的受限阶段, 其逃离受限的方式与管道宽度有关且满足不同的规律. 对于窄管道, 高分子链在输运过程中呈直线伸展构型且运动具有“蛇爬行”特征. 高分子链逃离受限过程伴随着整条链的运动, 从而导致迁移率随高分子链长呈周期变化, 而且在迁移率极值位置, 高分子链投影长度与管道半周期之间存在简单的整数倍关系. 对于宽管道, 高分子链在输运过程中出现弯折构型且运动具有“蠕虫运动”特征. 当链长比较长时, 高分子链可通过链前端部分的伸长逃离受限, 从而导致迁移率与高分子链长度无关. 模拟结果可能有助于利用周期管道对不同长度的高分子链进行分离及可控输运.     

SEGMENT DENSITY PROFILES OF COMPACT POLYMER CHAINS CONFINED BETWEEN TWO PARALLEL PLANE WALLS

We use the pruned-enriched Rosenbluth method to investigate systematically the segment density profiles of compact polymer chains confined between two parallel plane walls.The non-adsorption case of adsorption interaction energyε=0 and the weak adsorption case ofε=-1 are considered for the compact polymer chains with different chain lengths N and different separation distances between two walls D.Several special entropy effects on the confined compact polymer chains,such as a damped oscillation in the segment density profile for the large separation distance D,are observed and discussed for different separation distances D in the non-adsorption case.In the weak adsorption case,investigations on the segment density profiles indicate that the competition between the entropy and adsorption effects results in an obvious depletion layer.Moreover,the scaling laws of the damped oscillation period T_d and the depletion layer width L_d are obtained for the confined compact chains.Most of these results are obtained for the first time so far as we know,which are expected to understand the properties of the confined compact polymer chains more completely.     

Phase behavior of a blend of polymer-tethered nanoparticles with diblock copolymers

Using the self-consistent field theory (SCFT), we investigate the phase behavior of a mixture of diblock copolymers and nanoparticles with monodisperse polymer chains tethered to their surfaces. We assume the size of the nanoparticles to be much smaller than that of the attached polymer chains and therefore model the particles with their grafted polymer "shell" as star polymers. The polymer chains attached to the particles are of the same species as one of the blocks of the symmetric diblock copolymer. Of primary interest is how to tune the shell of the particle by changing both the length and number of tethered polymers in order to achieve higher loading of nanoparticles within an ordered structure without macrophase separation occurring. We find that the phase behavior of the system is very sensitive to the size of the particle including its tethered shell. The region of microphase separation is increased upon decreasing the star polymer size, which may be achieved by shortening and/or removing tethered polymer chains. To explore the possible structures in these systems we employ SCFT simulations that provide insight into the arrangement of the different species in these complex composites.     

Knotting of random ring polymers in confined spaces

Stochastic simulations are used to characterize the knotting distributions of random ring polymers confined in spheres of various radii. The approach is based on the use of multiple Markov chains and reweighting techniques, combined with effective strategies for simplifying the geometrical complexity of ring conformations without altering their knot type. By these means we extend previous studies and characterize in detail how the probability to form a given prime or composite knot behaves in terms of the number of ring segments N and confining radius R. For 50 < or =N < or =450 we show that the probability of forming a composite knot rises significantly with the confinement, while the occurrence probability of prime knots are, in general, nonmonotonic functions of 1R. The dependence of other geometrical indicators, such as writhe and chirality, in terms of R and N is also characterized. It is found that the writhe distribution broadens as the confining sphere narrows.     

Structure and Dynamics of Confined Polymer Melts from Attractive Interaction to Repulsive Interaction Between Polymer and Smooth Wall

We studied the static and dynamic properties of unentangled polymer chains which have a variable strength of interaction with the confining smooth walls by means of the lattice Monte Carlo simulation based on the bond-fluctuation model, that is, investigated the wall-polymer interactions which systematically vary from attraction to repulsion. A critical value of attractive potential(ε_wc) is found to be -0.6k_BT, and only below it can the adsorption layer of monomers be formed near the wall. At the critical point of attraction ε_wc, attractive interaction counterba- lances the wall-polymer excluded volume effect, which minimizes the confinement effects on both chain dimension and mobility. Influences on both chain dimension and mobility increase with the increasing of either attraction or repulsion imposed by the walls. Despite of the nature and strength of the wall-polymer interaction, with the decrease of film thickness, configurations more parallelly aligned and flattened are adopted by confined chains, and a systematic trend of deceleration is found. Variations of chain dynamics with both film thickness and wall-polymer interaction can be well explained by the corresponding changes in the confinement of the nearest-neighboring particles that surround the chains. Besides, the thickness of the interfacial layer inside polymer films, where chains adopt a flattened “pancake” shape, is about two times the bulk radius of gyration and independent of the wall-polymer interaction.     

Depletion-induced surface alignment of asymmetric diblock copolymer in selective solvents

Phase separation of asymmetric diblock copolymer near surfaces in selective solvents is theoretically investigated by using the real-space version of self-consistent field theory (SCFT). Several morphologies are predicted and the phase diagram is constructed by varying the distance between two parallel hard surfaces (or the film thickness) W and the block copolymer concentration f(P). Morphologies of the diblock copolymer in dilute solution are found to change significantly with different film thicknesses. In confined systems, stable morphologies found in the bulk solution become unstable due to the loss of polymer conformation entropy. The vesicle phase region contracts when the repulsive interaction between the blocks is strong (strong segregation regime). The mixture of vesicles, rodlike and spherelike micelles and the mixture of vesicles and sphere-like micelles disappear in contrast to the weakly segregating regime. The walls strongly affect the phase separation of block copolymer in selective solvent, and the depletion layer near the surface contributes much to the micelle formation of the block copolymer. Interestingly, the self-assembled morphologies stay near the walls with the distance on the order of the radius of gyration of the block copolymer. The oscillation of the polymer distribution near the walls allows the surface phase separation to be observed due to the strong repulsion between the blocks A and B.     

Effect of temperature on the structure of adsorption layers formed by adsorption of polymers in the transition region from dilute to semidilute solutions

The measurements of adsorption from solutions of polymers (poly(butyl methacrylate) and polystyrene) and their mixtures at various temperatures in the wide concentration region from dilute to semidilute solution have been made. The adsorption isotherms and fraction of bound segments confirm the existence of the transition concentration region near the critical concentration of the chain overlapping when the change of adsorption mechanisms take place. The effect of temperature on adsorption from the mixtures is different for both the polymers depending on the concentration regime. For PBMA, which adsorbs from the mixtures preferentially, the adsorption isotherms are of the same shape as for adsorption from binary solutions and are characterized by the presence of a small forepart in the transition region and by increasing adsorption with temperature. For less adsorbing PS by transition from dilute to semidilute regime the inversion of the temperature effect on adsorption is observed. These results are confirmed by the estimations of the parameters of preferential adsorption at various temperatures. The data on the fraction of bound segments for dilute regime corresponds to the extended conformation of chains at the surface. The transition to semidilute regime leads to the diminishing of the fraction of bound segments as a result of simultaneous adsorption of macromolecular aggregates. Values of adsorption layer thickness have been calculated for various solution regimes and concentration. The dependence of the adsorption layer thickness on the temperature and on the solution regime at which adsorption occurs was established.     

Solvent size effect on the depletion layer of a polymer solution near an interface

By varying polymer concentration phi0p and Flory-Huggins parameter chi, the effect of solvent size on the depletion interaction between polymer coils and a hard wall was investigated by the real-space version of self-consistent field theory (SCFT). The depletion profiles and depletion thickness indicated that the depletion effect is strong in less good solvent with large molecular volume. Through the analysis of the respective free energies of polymer coils and solvent molecules, we found that the increment in the translation entropy of the solvent is the key to strengthening the depletion interaction. On the basis of the SCFT results, we define a solvent with volume about one to six times that of the polymer segment as a "middle-sized solvent". The density oscillations previously studied by Van der Gucht et al. and Maassen et al. were also observed in our simulation, and the addition of middle-sized solvent will magnify the amplitude of the oscillations. The solvent-size-dependent depletion interaction may be an explanation for the reduced entanglement and promoted crystallization behavior of polymer coils prepared from the solution with middle-sized solvent.     

Adsorption of comb copolymers on weakly attractive solid surfaces

In this work continuum and lattice Monte Carlo simulation methods are used to study the adsorption of linear and comb polymers on flat surfaces. Selected polymer segments, located at the tips of the side chains in comb polymers or equally spaced along the linear polymers, are attracted to each other and to the surface via square-well potentials. The rest of the polymer segments are modeled as tangent hard spheres in the continuum model and as self-avoiding random walks in the lattice model. Results are presented in terms of segment-density profiles, distribution functions, and radii of gyration of the adsorbed polymers. At infinite dilution the presence of short side chains promotes the adsorption of polymers favoring both a decrease in the depletion-layer thickness and a spreading of the polymer molecule on the surface. The presence of long side chains favors the adsorption of polymers on the surface, but does not permit the spreading of the polymers. At finite concentration linear polymers and comb polymers with long side chains readily adsorb on the solid surface, while comb polymers with short side chains are unlikely to adsorb. The simple models of comb copolymers with short side chains used here show properties similar to those of associating polymers and of globular proteins in aqueous solutions, and can be used as a first approximation to investigate the mechanism of adsorption of proteins onto hydrophobic surfaces.