Variational theory for a single polyelectrolyte chain

Variational methods are applied to a single polyelectrolyte chain. The polymer is modeled as a Gaussian chain with screened electrostatic repulsion between all monomers. As a variational Hamiltonian, the most general Gaussian kernel, including the possibility of a classical or mean polymer path, is employed. The resulting self-consistent equations are systematically solved both for large and small monomer-monomer separations along the chain. In the absence of screening, the polymer is stretched on average. It is described by a straight classical path with Gaussian fluctuations around it. If the electrostatic repulsion is screened, the polymer is isotropically swollen for large separations, and for small separations the polymer correlation function is calculated as an analytic expansion in terms of the monomer-monomer separation along the chain. The electrostatic persistence length and the electrostatic blobsize are inferred from the crossover between distinct scaling ranges. We perform a global analysis of the scaling behavior as a function of the screening length and electrostatic interaction strength , where is the Bjerrum length and A is the distance of charges along the polymer chain. We find three different scaling regimes. i) A Gaussian-persistent regime with Gaussian behavior at small, persistent behavior at intermediate, and isotropically swollen behavior at large length scales. This regime occurs for weakly charged polymers and only for intermediate values of the screening length. The electrostatic persistence length is defined as the crossover length between the persistent and the asymptotically swollen behavior and is given by and thus disagrees with previous (restricted) variational treatments which predict a linear dependence on the screening length .ii) A Gaussian regime with Gaussian behavior at small and isotropically swollen behavior at large length scales. This regime occurs for weakly charged polymers and/or strong screening, and the electrostatic repulsion between monomers only leads to subfluent corrections to Gaussian scaling at small separations. The concept of a persistence length is without meaning in this regime. iii) A persistent regime , where the chain resembles a stretched rod on intermediate and small scales. Here the persistence length is given by the original Odijk prediction, , if the overstretching of the chain is avoided. We also investigate the effects of a finite polymer length and of an additional excluded-volume interaction, which modify the resultant scaling behavior. Applications to experiments and computer simulations are discussed. Received 24 December 1997     

Grafted polymers with annealed excluded volume: A model for surfactant association in brushes

We present an analytical self-consistent-field (SCF) theory for a neutral polymer brush (a layer of long polymer chains end-grafted to a surface) with annealed excluded volume interactions between the monomer units. This model mimics the reversible adsorption of solute molecules or aggregates, such as small globular proteins or surfactant micelles, on the grafted chains. The equilibrium structural properties of the brush (the brush thickness, the monomer density profile, the distribution of the end segments of the grafted chains) as well as the overall adsorbed amount and the adsorbate density profile are analyzed as a function of the grafting density, the excluded volume parameters and the chemical potential (the concentration) of the adsorbate in the solution. We demonstrate that, when the grafting density is varied, the overall adsorbed amount always exhibits a maximum, whereas the root-mean-square brush thickness either increases monotonically or passes through a (local) minimum. At high grafting densities the chains are loaded by adsorbed aggregates preferentially in the distal region of the brush, whereas in the region proximal to the grafting surface depletion of aggregates occurs and the polymer brush retains an unperturbed structure. Depending on the relative strength of the excluded volume interactions between unloaded and loaded monomers both the degree of loading of the chains and the polymer density profile are either continuous or they exhibit a discontinuity as a function of the distance from the grafting surface. In the latter case intrinsic phase separation occurs in the brush: the dense phase consists of unloaded and weakly extended chains and occupies the region proximal to the surface, whereas a more dilute phase consisting of highly loaded and strongly extended chains forms the periphery of the brush. Received 26 November 1998 and Received in final form 2 April 1999     

Polymer adsorption on heterogeneous surfaces

The adsorption of a single ideal polymer chain on energetically heterogeneous and rough surfaces is investigated using a variational procedure introduced by Garel and Orland (Phys. Rev. B 55, 226 (1997)). The mean polymer size is calculated perpendicular and parallel to the surface and is compared to the Gaussian conformation and to the results for polymers at flat and energetically homogeneous surfaces. The disorder-induced enhancement of adsorption is confirmed and is shown to be much more significant for a heterogeneous interaction strength than for spatial roughness. This difference also applies to the localization transition, where the polymer size becomes independent of the chain length. The localization criterion can be quantified, depending on an effective interaction strength and the length of the polymer chain. Received: 29 October 1997 / Revised: 19 December 1997 / Accepted: 6 March 1998     

考虑扩散和吸附作用的聚合物驱替过程渗流数值模拟

采用显式求解饱和度、隐式求解浓度的思路,对考虑扩散和吸附作用的聚合物驱替过程渗流模型进行了数值求解.饱和度方程求解应用了显式全变差递减(TVD)法;浓度方程求解过程中,空间项离散采用Crank Nicolson差分格式,时间项变量进行拟线性处理,保证了计算的稳定性.通过与解析解对比,验证了该方法的有效性.计算实例分析表明,扩散使聚合物在溶液中稀释,导致浓度传播分散;吸附使聚合物损耗,导致浓度传播滞后.同时,计算结果直观反映了聚合物驱重要的"油墙"形成机理.在段塞注入情况下,原油富集区在出口端的突破是介于聚合物浓度前缘突破和聚合物浓度峰值突破之间.     

Field theoretical approach to inhomogeneous ionic systems: thermodynamic consistency with the contact theorem,Gibbs adsorption and surface tension

Using a statistical field approach we investigate the structure of an electrolyte solution in contact with a neutral impenetrable wall. The Hamiltonian contains the Coulomb interaction and the ideal entropy. At the level of the quadratic approximation, the Hamiltonian yields the Debye-Hückel theory in the bulk. Analytic expressions of the charge-charge and potential-potential inhomogeneous correlation functions are obtained. Exact asymptotic results for point ion charge correlation functions are obtained and the profile for the fluctuation of the electric potential is calculated. We also consider the term beyond the quadratic expansion of the ideal entropy in the Hamiltonian. With this term a higher order coupling between charge density and number density produces a non-trivial profile for the total ion density. This density profile is consistent with the contact theorem and the related surface tension calculated from the Gibbs adsorption isotherm.     

Local entropic effects of polymers grafted to soft interfaces

In this paper, we study the equilibrium properties of polymer chains end-tethered to a fluid membrane. The loss of conformational entropy of the polymer results in an inhomogeneous pressure field that we calculate for Gaussian chains. We estimate the effects of excluded volume through a relation between pressure and concentration. Under the polymer pressure, a soft surface will deform. We calculate the deformation profile for a fluid membrane and show that close to the grafting point, this profile assumes a cone-like shape, independently of the boundary conditions. Interactions between different polymers are also mediated by the membrane deformation. This pair-additive potential is attractive for chains grafted on the same side of the membrane and repulsive otherwise. Received 20 April 2000     

Thermodynamics of polystyrene solutions. I. Interaction parameters and polymer order

The Maron theory was used to deduce from mostly osmotic pressure and light-scattering data the interaction parameters as a function of concentration, temperature, and polymer molecular weight for solutions of polystyrene in seven different solvents. These interaction parameters were then utilized to calculate the activity coefficients of the solvent in more concentrated solutions for comparison with experimental values deduced from vapor pressure measurements. When order was absent in the polymers (all cases but two), the agreement between calculated and observed activity coefficients was found to be very good. In the two cases where order was present, the activity coefficients observed in the concentrated solutions were higher than those calculated. These deviations were employed to find the degree of polymer order both in solution and in the original polymers.     

The state space of short-range Ising spin glasses: the density of states

The correlations between the segments of a semidilute polymer solution are found to induce correlations in the positions of small particles added to the solution. Small means a diameter much less than the polymer's correlation length. In the presence of polymer the particles behave as if they attracted each other. It is shown how the polymer's correlation length may be determined from a scattering experiment performed on the spheres. Received: 7 July 1997 / Received in final form: 12 November 1997 / Accepted: 19 November 1997     

Adsorption of polymers anchored to membranes

Polymers, which are attached to a membrane at one of their ends, exert an entropic pressure, which curves the membrane away from the polymers. It is shown that adsorption which arises from a short-ranged potential between the polymer and the membrane has a large influence on the curvature of the membrane, leading to a decrease of the entropically induced curvature. If one ignores the finite size of the anchor segment, the polymer-induced curvature does not change sign and vanishes in the limit of strong adsorption and a pure contact potential. If one includes the finite size of the anchor segment, the membrane bends towards the polymer for sufficiently strong adsorption. Received 20 October 2000     

快中子照相中的点扩展函数计算

 快中子照相中,基于反冲核原理探测快中子的有机闪烁体平板是普遍采用的快中子辐射转换体。模拟了D-T中子垂直入射BC400闪烁体平板,计算了不同厚度平板闪烁体的点扩展函数,对14.1 MeV快中子照相中闪烁体固有分辨率随厚度的变化进行了研究。计算结果表明,在不考虑二次中子与闪烁体作用及背景噪声等情况时,点扩展函数几乎不依赖于闪烁体厚度。同时,计算还表明在一定的分辨率范围内,由于荧光收集效率的限制,闪烁体厚度增加并不会改善图像对比度。     

快中子照相中的点扩展函数计算

快中子照相中,基于反冲核原理探测快中子的有机闪烁体平板是普遍采用的快中子辐射转换体。模拟了D-T中子垂直入射BC400闪烁体平板,计算了不同厚度平板闪烁体的点扩展函数,对14.1 MeV快中子照相中闪烁体固有分辨率随厚度的变化进行了研究。计算结果表明,在不考虑二次中子与闪烁体作用及背景噪声等情况时,点扩展函数几乎不依赖于闪烁体厚度。同时,计算还表明在一定的分辨率范围内,由于荧光收集效率的限制,闪烁体厚度增加并不会改善图像对比度。     

Behavior of a wetting phase near a solid boundary: vapor near a weakly attractive surface

Density profiles of a LJ vapor near a weakly attractive surface with long-range fluid wall potential was studied along the pore coexistence curve. There are two localized density maxima near the pore wall: the first one is caused by localization of the molecules in the minimum of the fluid-wall potential, and the second one reflects adsorption of molecules at the first layer at higher densities. In addition, a third, weak density maximum is observed close to the critical temperature due to the competition between the long-range attractive tail of the fluid-wall potential and the effect of missing neighbors. This maximum separates the region of a gradual density depletion toward the surface due to the missing neighbor effect and the adsorption region further from the surface, where the density gradually increases toward the surface due to the attractive fluid-wall potential. When approaching the bulk critical temperature, this maximum moves away from the surface due to the divergence of the bulk correlation length. Applicability of various equations to describe the vapor density profiles is examined. Excess adsorption of vapor at low temperatures turns into excess depletion at higher temperatures. The crossover temperature increases with increasing pore size and strengthening fluid-wall interaction. The problems of the theory of the surface critical behavior of Ising models in a case of a non vanishing surface field and its mapping on a fluid is discussed.     

Analysis of polymer adsorption onto colloidal particles

The structure of the layer formed after polymer adsorption onto a spherical particle is numerically studied by means of the application of the Single-Chain Mean-Field theory. We have determined several overall layer properties including the monomer volume fraction profiles, the layer thickness, adsorbances related to loops and to tails, as well as the variation of the crossover distance between loops and tails for different particle radii and fixed polymer length. When the radius of the sphere is small enough to affect the loop layer, one enters a single-adsorbed-chain regime, characterized by a critical sphere radius. In this regime, structural changes in the adsorbed layer arise. For such small sphere, the loop layer is confined to a region whose thickness is of the order of the radius of the adsorbing sphere, and two long tails dominate the outer layer and the adsorbance due to tails dominates that due to loops. An analysis of the structure of the outer tail layer for this small sphere case is also presented.     

Phase separation of polymer blends in solution near adsorbing surface

In this paper, we consider a mixture of two polymers A and B of different chemical nature, dissolved in a common good solvent, in contact with an interacting surface. We start from a mixture of two incompatible homopolymers A and B in the molten state, and assume that the surface adsorbs strongly one or both polymer species at high temperature. It is assumed that this is a strong adsorption, so that chains cannot desorb once they are linked to the surface. This constrains the system to a quenched composition on the surface. Once the adsorption process is finished, a quantity of a good solvent is added to get a semi-dilute solution. We assume that demixing transition in the presence of solvent occurs at lower temperature. The purpose is to discuss the influence of the quenched surface fluctuations on the critical properties of the mixture. Within the framework of the so-called blob model, we determine the exact shape of the composition profile as a function of the distance z to the surface, for any value of the relevant parameters, namely, the temperature T, the molecular weight M, the monomer concentration c and the surface composition x₀. Our analysis reveals a universal character of the composition profile for , where the characteristic size D is some known length depending on the relevant parameters of the problem, and not on temperature, and is the thermal correlation length. Near surface, for (a is the monomer size), the profile is no longer universal, and in particular, it is sensitive to the boundary condition. Far from the surface, that is , the profile tends exponentially to its bulk value. We show that the length Dapproaches its lowest value as the surface composition reaches its saturated value l. In this limit, we find that the profile shape is a characteristic of critical adsorption in simple binary fluid mixtures. Finally, this work must be regarded as a natural extension of a previous one, which was concerned with the same problem, but in the absence of solvent. Received 24 June 1999 and Received in final form 5 November 1999     

Free energy of semiflexible polymers and structure of interfaces

The free energy of semiflexible polymers is calculated as a functional of the compositional scalar order parameter and the orientational order parameter of second-rank tensor S_ij on the basis of a microscopic model of wormlike chains with variable segment lengths. We use a density functional theory and a gradient expansion to evaluate the entropic part of the free energy, which is given in a power series of .The interaction term of the free energy is derived with a random phase approximation. For the rigid rod limit, the nematic-isotropic transition point is given by , N and w being the degree of polymerization and the anisotropic interaction parameter, respectively, and the degree of ordering at the transition point is 0.33448. We also find that the contour length of polymer chains becomes larger in a nematic phase than in an isotropic phase. Interface profiles are obtained numerically for some typical cases. In the neighborhood of isotropic-isotropic interfaces, polymer chains tend to align parallel to the interface on the polymer-rich side and perpendicular on the poor side. When an isotropic region and a nematic region coexist, orientational order parallel to the interface is preferred in the nematic region. Received: 28 May 1998 / Revised: 12 August 1998 / Accepted: 8 September 1998     

熵驱动下柱状高分子刷螺旋结构的研究

自然界中广泛存在螺旋结构,在特定情形下熵能驱动高分子链形成螺旋结构.本文采用分子动力学方法研究柱状高分子刷吸附在无限长圆柱表面时的构象行为.发现其构象与嫁接支链条数、柱状高分子刷与圆柱表面之间的吸附能密切相关.在较弱的吸附能下,具有较多支链条数的柱状高分子刷能形成完整的螺旋结构,其本质就是熵驱动下形成的螺旋结构.该研究有助于加深对生物大分子螺旋结构的理解.     

外力驱动作用下高分子链在表面吸附性质的计算机模拟

采用退火法模拟研究受外力F驱动的高分子链在吸引表面的吸附特性.通过高分子链的平均表面接触数〈M〉与温度T之间的关系计算临界吸附温度T_c,并发现T_c随着F的增加而减小;进而通过高分子链的均方回转半径分析外力驱动作用对高分子链构象的影响,并从回转半径极小值或者垂直外力方向的y和z分量的变化交叉校验临界吸附点T_c.模拟计算了处于吸附状态的高分子链随着外力F的增加是否会发生吸附状态到脱附状态的相变以及发生相变所需施加的外力是否由温度所决定.模拟结果表明:两种不同温度下高分子链的吸附性质和构象性质受外力驱动作用而产生不同现象,在温度区间T*_cTT_c时会发生脱附现象,而在TT*_c时不会发生脱附现象.