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

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

受限于平行界面之间高分子链力学行为的研究

采用PERM(pruned-enriched-rosenbluth method)算法,和简立方格点上的自避行走模型,研究了受限于平行界面之间的星型高分子链和线型高分子链的力学行为,并将两种结果进行了比较.模拟结果表明,如果界面对高分子链存在吸附作用,界面之间距离D的增大需要外力的作用;如果界面对链没有吸附作用,则D的增大是一个自发的过程,无需外力的作用.随着界面间距D的增大,通过计算星型链和线型链的均方回转半径〈S2〉和形状因子〈δ*〉,研究了链尺寸和形状的变化情况.另外,为了更细微地了解星型链和线型链的结构及力学行为,还研究了受限高分子链在平行界面之间的轨链(train)、环链(bridge)和尾链(tail).     

支链对高分子链导热性能影响的分子动力学模拟

采用非平衡分子动力学(NEMD)方法模拟了立构规整的乙丙交替共聚物(alternating isotactic PEP)和聚乙烯(PE)分子链的热传导过程,研究了支链对高分子链导热率的影响,并通过分析均方回转半径、径向分布函数以及均方位移与导热率的关系,进一步探讨了高分子链中热输运的微观机理.通过比较发现,主链上含支链的高分子链导热率较低;均方回转半径显示,高分子链的构象越稳定,导热率越高;径向分布函数显示,主链上碳原子分布越紧密,导热率越低;均方位移分析结果表明,主链上的支链使高分子链中的原子运动加剧,从而导致导热率降低.     

平行板间高分子链构象熵的动力学Monte Carlo模拟

用动力学Monte Carlo方法模拟了受限于两平行板之间的高分子链,并用扫描法计算了链的构象熵S,研究了构象熵相对于自由链的减小量(S0-S)与平行板间距D和高分子链长n的关系.结果证实了de Gennes的自由能标度关系,并给出了标度关系适用的范围.当D非常小时,高分子链受到强烈限制,S0-S与n成正比,表明单链节受到平行板的平均排斥作用力与链长无关.随着D增大,平行板对构象熵的影响越来越弱,单链节受到平行板的平均排斥作用力随链长的增长而增大.当D比较大时,平行板对构象熵的影响近似可以忽略,高分子链构象熵与自由空间中的结果一致.     

溶剂性质对星形聚电解质构象的影响

采用分子动力学方法模拟仿真溶剂对星形聚电解质构象的影响.选择的溶剂是从良溶剂到不良溶剂.研究了溶剂性质对单体的密度曲线、均方末端距、均方旋转半径、均方位移和系统渗透压的影响,结果表明星形聚电解质的构象随着溶剂参数β的增加而逐渐收缩,另外随着溶剂参数β的不断变大,系统中反离子的扩散能力不断下降,反离子凝聚能力增强,系统的...     

良溶剂中接枝型聚两性电解质单链构象的布朗动力学模拟

采用布朗动力学研究了在良溶剂中荷电平衡的接枝聚两性电解质(GPA)的单链构象转变行为,讨论了主链链长、支链数及电荷密度对GPA分子链构象转变的影响.研究发现,随着静电相互作用的增强,GPA分子链构象转变过程由线团、主链与支链间的折叠、链段塌缩和电荷配对形成偶极子与四极子等4个阶段构成.与线型聚两性电解质不同,GPA存在的额外支链间空间排斥与静电排斥作用随着分子结构的变化而改变,并影响构象转变行为.在强静电相互作用下,良溶剂中的GPA链由于溶剂化作用会再伸展,以保证偶极子完全配对成四极子.减小主链长度或电荷密度或增加支链数目都会增大体系的排斥力和主链的刚性,阻滞分子链的塌缩,并使得分子链再伸展的幅度增大.     

聚乙烯醇水溶液反复冷冻过程中的溶剂化效应

采用示差扫描量热技术(DSC)对聚乙烯醇(PVA)水溶液反复冰冻过程中的溶剂化效应进行研究.引入水化数的概念来表征溶剂化效应的大小.结果表明不同浓度区间的PVA水溶液其在反复冰冻过程中溶剂化效应显著不同,主要归因于高分子链分子内和分子间缠结程度对溶剂分子"参与"溶剂化的程度和方式的不同.作者把极稀高分子溶液的研究结果拓展到高分子稀溶液或亚浓溶液区间,阐述了高分子溶液中高分子链的物理图像.冷冻次数的增加导致链间缠结增加,部分溶剂则被包裹在由链间缠结点所形成的网圈内成为分子链的一部分.溶液溶剂化程度的变化受到包裹溶剂与高分子链脱溶剂化的综合影响.     

不同螺旋结构类蛋白质分子管道传输动力学研究

采用动态蒙特卡洛模拟(Dynamical Monte Carlo,DMC)对不同螺旋结构的类蛋白质分子穿越管道的动力学过程进行了研究.计算了不同螺旋能参数下类蛋白质链的传输时间,传输成功率,管道内单体数目,形状因子<δ>,均方位移g等.结果发现螺旋能较小的链传输时间较短,变化不明显,而螺旋能较大的链随着螺旋能的增强传输时间迅速增大.计算表明相同类型的分子链其传输成功率随着电场的增大而减小.而对于不同螺旋结构的链,螺旋能参数εh=0的链传输成功率最大,达到80%左右,略小于自避链;而εh=-0.5的类蛋白质链传输成功率突然减小.随着螺旋能的增强,传输成功率又开始增大.另外的计算表明自避链的传输过程比类蛋白质链平稳.<δ>在传输过程中的变化都存在了最大值,说明链在传输时有一个最松散的构象.对于均方位移的计算得到εh=0的链均方位移变化最快,说明传输过程中扩散最快.同时发现εh=-3.0的链比εh=-1.5的链均方位移扩散快.形状因子和均方位移的趋势均能很好的解释不同螺旋参数下的类蛋白质链不同的传输性质.这些研究对生物大分子的输运行为有一定的指导意义.     

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

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

高分子链形状与尺寸关联的Monte Carlo模拟

运用MonteCarlo方法对线型高分子链格点模型的构型进行了模拟,研究了构型的尺寸(采用平方末端距R²,平方回转半径S²来表征)和形状(由非球形因子A表征)之间的关联.对任何长度的高分子链,其关联系数CA,R²和CA,S²均为正值,表明高分子链的形状与尺寸之间存在正关联,即尺寸小的构型其非球形因子A一般也小,反之尺寸大的构型其非球形因子A一般也大.关联系数CA,R²和CA,S²均随链长的增大而减小,近似地与链长的倒数(n^-1)成正比.研究还表明,关联系数的极限值(链长n很大时)与格点的类型无关,与链样本产生的方式也无关,但与链是否考虑排斥体积有关,考虑了排斥体积后,关联系数增大.     

三维吸附紧密高分子单链的力学性质研究

采用PERM(pruned-enriched Rosenbluth method)算法,研究了吸附在界面附近的紧密高分子链力学行为.发现当界面的吸附能比较大时,紧密高分子链从紧贴于吸附界面到逐渐远离的过程中,其外形会经历4种典型的变化.同时紧密高分子链的尺寸大小如/N、xy/N、z/N,形状参数<δ*>,热力学性质如每个键的平均自由能A/N,平均相互作用能/N等,甚至所受外力的大小都会同时做出相应的变化,其出现变化的位置也一致.特别是随着紧密高分子链离开吸附界面的过程中,作用于高分子链上的外力明显出现几个力学平台,这与实验得到的结果完全一致.同时还研究了弱吸附能的情况,在这种情况下实验是很难进行的.     

Synthesis, helicity, and chromism of optically active poly(phenylacetylene)s carrying different amino acid moieties and pendant terminal groups

Functional phenylacetylene derivatives containing l-alanine and l-leucine moieties with chiral menthyl and achiral n-octyl terminal groups {HC[triple bond]C-C6H4-p-CONHCH(R)CO2R': R = CH3, R'= (-)-(1R,2S,5R)-menthyl [1(-)]; R = CH2CH(CH2)3, R' = (-)-(1R,2S,5R)-menthyl [2(-)]; R'= CH2CH(CH2)3, R' = (+)-(1S,2R,5S)-menthyl [2(+)]; R'= CH2CH(CH2)3, R' = (CH2)7CH3 (2o)} are synthesized. Polymerizations of the acetylene monomers are effected by organorhodium catalysts, giving corresponding polymers P1(-), P2(-), P2(+), and P2o of high molecular weights (Mw up to 1.2 x 10(6)) in high yields (up to 89%). The polymers are thermally stable (Td >or= 300 degrees C) and soluble in common organic solvents. The polymer structures are characterized by IR, NMR, UV, and CD spectroscopies. Intense CD signals are observed in the visible spectral region, indicating that the polymer chains are taking a helical conformation with an excess of preferred handedness. The backbone conjugation and chain helicity of the polymers can be tuned by changing their molecular structures [(a)chiral pendant groups] and by applying external stimuli (solvent and pH). Addition of trifluoroacetic acid to the polymer solutions decreases their molar ellipticities and enhances their backbone conjugations, inducing a halochromism with a continuous and reversible color change (yellow <==> red).     

疏水缔合聚丙烯酰胺在盐水中的自组装行为

使用静态光散射、动态光散射以及原子力显微镜(AFM)研究了疏水缔合聚丙烯酰胺(HAPAM)在盐水溶液中的自组装行为.研究了聚合物分子在不同盐浓度中的表观重均分子量(Mw,a),均方根回转半径(),流体力学半径(),第二维里系数(A2)的变化,并根据/的比值得出聚合物分子的聚集形态.实验结果表明:随着盐浓度的增强,聚合物溶液的分子链由舒展变得卷曲.     

Elastic behavior of adsorbed polymer chains

Elastic behaviors of single polymer chains adsorbed on the attractive surface are first investigated using Monte Carlo simulation method based on the bond fluctuation model. We investigate the chain size and shape of adsorbed chains, such as mean-square radius of gyration S2, mean-square bond length b2, shape factors sf(i) and delta*, and the orientation of chain segments P2, to illuminate how the shape of polymer chains changes during the process of tensile elongation. There are some special behaviors of the chain size and shape at the beginning of elongation, especially for strong attraction interaction. For example, mean fraction of adsorbed segments decreases abruptly in the region of small elongation ratio and then decreases slowly with increasing elongation ratio. In fact, the chain size and shape also changes abruptly for small elongation ratio with strong attraction interaction. Some thermodynamics properties are also investigated here. Average Helmholtz free energy increases fast for elongation ratio lambda<1.15, especially with strong attraction, and increases slowly for lambda>1.15. Similar behaviors are obtained for average energy per bond. Elastic force (f ) and energy contribution to force (f(U)) are also studied, and we find that elastic force decreases abruptly for lambda<1.15, and there is a minimum of elastic force for strong attraction interaction, then increases very slowly with increasing elongation ratio. However, there are different behaviors for weak attraction interaction. For energy contribution to force (f(U)), there is a maximum value for strong attraction interaction in the region of lambda<1.15. Some comparisons with the atomic force microscopy experiments are also made. These investigations may provide some insights into the elastic behaviors of adsorbed polymer chains.     

Effect of solvent quality on the conformations of a model comb polymer

The effect of solvent quality on the equilibrium structure of a densely branched comb polymer is investigated based on the structure factor analyses by off-lattice Monte Carlo simulations. First, theta temperature (theta(infinity)) must be determined to identify the solvent condition. We locate the characteristic temperature theta(A)(N) at which the second virial coefficient vanishes and the transition temperature theta(R)(N) at which radius of gyration R(g) of the chain varies most rapidly with temperature, i.e., d(2)R(g)/dT(2)|(theta(R)) = 0. N represents the total number of monomers of a comb. As N --> infinity, theta(A) and theta(R) coincide to a point that is identified as the true theta temperature (theta(infinity)). The structure factors of the main chain, the side chain, and the whole polymer are calculated, respectively. It is found that at T = theta(infinity), the structural factors S(qR(g)) for the overall comb polymers match quite well with those of their Gaussian counterparts. When T< theta(infinity), the overall comb polymer assumes collapsed conformations, similar to a homogeneous sphere. However, the structure factor of the side chain indicates that it always remains in an expanded state regardless of the solvent condition. It is attributed to the strong interactions between side chains. The same effect leads to enhanced rigidity of the main chain in comparison to the linear chain, as clearly observed from the rescaled Kratky plot.     

Contributions of short-range and excluded-volume interactions to unperturbed polymer chain dimensions

A Monte Carlo study is made of the mean-square radius of gyration for the freely rotating chain with such fictitious excluded-volume interactions that the Lennard-Jones 6-12 potentials at the Theta temperature act only between the fourth-through (3+Delta)th-neighbor beads (Delta > or = 1) along the chain. The behavior of the asymptotic value (/n)infinity of the ratio /n as a function of the number n of bonds in the chain in the limit of n --> infinity is examined as a function of Delta. It is shown that the approach of (/n)infinity to its value for the real unperturbed chain with Delta = infinity is so slow that the interactions between even up to about 100th-neighbor beads should be taken into account in order to reproduce nearly its dimension. The result implies that the unperturbed polymer chain dimension as experimentally observed at the Theta temperature depends not only on short-range interactions but also to a considerable extent on the long-range excluded-volume interactions, and that the asymptotic value Cinfinity of the characteristic ratio Cn for the rotational isomeric state model in the limit of n --> infinity, which is determined only by the very local conformational energy, cannot be directly compared with the corresponding experimental value.     

Molecular Dimension and Interaction Parameters of Polyacrylamide in Water‐N,N‐Dimethylformamide Mixtures

The intrinsic viscosities [η] of polyacrylamide (PAM) having different average molecular weights are measured in various mixtures of water (good solvent) and N, N dimethyl formamide (DMF, nonsolvent) at different temperatures. The observed results show a significant variation of cosolvency as a function of solvent composition (?_DMF). The nature of curves in [η] vs. ?_DMF plot at different temperatures indicates the existence of two antagonistic effects. The unperturbed dimensions (K_θ) of the polymer are determined by a number of methods, which agree well with each other. The temperature coefficient of unperturbed dimension (K′), molecular extension factors (α_n), characteristic ratio (C_α) and chain rigidity (σ) are evaluated and the effects of temperature, solvent composition are discussed. The volume related parameter V_E and shape factor ν were also computed, which shows the shape of polymer molecules to be more or less spherical in solution.     

Synthesis and esterolytic reactions of linear and cross-linked asymmetric imidazole-containing polymers

Asymmetric linear and cross-linked imidazole-containing polymers were prepared from 1-[2(S)-methylbutyl]-4-vinylimidazole and 1-[2(S)-methylbutyl]-5-vinylimidazole. The esterolytic reactions of these linear and cross-linked asymmetric polymers with the enantiomeric substrates (R)- and (S)-4-(3-methylpentadecanoyl)-3-nitrobenzoic acid, (R)- and (S)-S, were studied by measuring the pseudo-first order kinetics of the solvolysis of these enantiomeric substrates in the presence of these asymmetric polymers. The linear homopolymers and copolymers of 1-[2(S)-methylbutyl]-4- and 5-vinylimidazole showed hydrophobic and electrostatic effects in the solvolysis of the enantiomeric substrates with these linear asymmetric polymers. Cross-links were introduced into these asymmetric polymers to increase the rigidity and reduce the number of conformations available to the polymer. The reduced conformational mobility was expected to enhance any enantioselectivity in the solvolysis of the enantiomeric substrates with these polymers. Using these cross-linked polymers, hydrophobic interactions were observed in the solvolysis of a series of substrates with increasing alkyl chain length. Also, on changing the polarity of the solvent, a bell-shaped rate profile was observed in the solvolysis of the long chain substrate S. This effect was attributed to a combination of the coiling of the cross-linked polymer chains and hydrophobic interactions with the substrate on changing solvent polarity. Even with the increased rigidity of these cross-linked polymers, no significant enantioselectivity in the solvolysis of the enantiomeric substrates was observed. So, neither the linear nor the cross-linked asymmetric polyvinylimidazoles showed enantioselectivity in the solvolysis of these enantiomeric substrates. In this case, the hydrophobic interactions and the reduced conformational mobility through crosslinking were not strong enough to bring about enantioselectivity in the solvolysis of these enantiomeric substrates.     

关于聚乙烯分子链均方回转半径的进一步讨论

本文考虑到氢原子的影响, 用矩阵代数的方法重新推导了聚乙烯分子链的均方回转半径。在我们得到的公式中若忽略氢原子的质量时, 就简化到Flory的公式。具体数值计算得到聚乙烯分子链均方回转半径对分子量的依赖关系是:〈S^2〉^1^/^2=0.44M^1^/^2比较θ-溶剂的实验数据: ⅹ〈S^2〉^1^/^2~w=0.45±0.08M^1^/^2~w是符合得很好的。