A markov chain Monte Carlo method for estimating population mixing using Y-chromosome markers: Mixing of the Han people in China

We present a new approach for estimating mixing between populations based on non-recombining markers, specifically Y-chromosome microsatellites. A Markov chain Monte Carlo (MCMC) Bayesian statistical approach is used to calculate the posterior probability

Bayesian inference of population expansions in domestic bovines

< p > The past population dynamics of four domestic and one wild species of bovine were estimated using Bayesian skyline plots, a coalescent Markov chain Monte Carlo method that does not require an assumed parametric model of demographic history. Four dom

Effects of chain flexibility on polymer conformation in dilute solution studied by lattice Monte Carlo simulation

Effects of chain flexibility on the conformation of homopolymers in good solvents have been investigated by Monte Carlo simulation. Bond angle constraint coupled with persistence length of polymer chains has been introduced in the modified eight-site bond fluctuation simulation model. The study about the effects of chain flexibility on polymer sizes reveals that the orientation of polymer chains under confinement is driven by the loss of conformation entropy. The conformation of polymer chains undergoing a gradual change from spherical iso-diametric ellipsoid to rodlike iso-diametric ellipsoid with the decrease of polymer chain flexibility in a wide region has been clearly illustrated from several aspects. Furthermore, a comparison of the freely jointed chain (FJC) model and the wormlike chain (WLC) model has also been made to describe the polymer sizes in terms of chain flexibility and quasi-quantitative boundary toward the suitability of the models.

Monte Carlo simulation of miscibility of polymer blends with repulsive interactions: effect of chain structure

The effects of the chain structure and the intramolecular interaction energy of an A/B copolymer on the miscibility of the binary blends of the copolymer and homopolymer C have been studied by means of a Monte Carlo simulation. In the system, the interactions between segments A, B and C are more repulsive than those between themselves. In order to study the effect of the chain structure of the A/B copolymer on the miscibility, the alternating, random and block copolymers were introduced in the simulations, respectively. The simulation results show that the miscibility of the binary blends strongly depends on the intramolecular interaction energy ((ε) over bar (AB)) between segments A and B within the A/B copolymers. The higher the repulsive interaction energy, the more miscible the A/B copolymer and homopolymer C are. For the diblock copolymer/homopolymer blends, they tend to form micro phase domains. However, the phase domains become so small that the blend can be considered as a homogeneous phase for the alternating copolymer/ homopolymer blends. Furthermore, the investigation of the average end-to-end distance ((h) over bar) in different systems indicates that the copolymer chains tend to coil with the decrease Of (ε) over bar (AB) whereas the (h) over bar of the homopolymer chains depends on the chain structure of the copolymers.

高分子溶液中物理凝胶化的Monte Carlo模拟和小角中子散射研究

高分子凝胶广泛地存在于自然界以及日常生活中，按其形成作用力不同分为化学凝胶和物理凝胶两大类。由于高分子物理凝胶具有凝胶化的可逆性及其对环境条件强烈的响应性，因此，在近半个世纪的研究与应用中受到极大的关注。高分子溶液中的物理凝胶因其结构及形成机制复杂，在实验方面，除了散射技术及流变技术能够有效地揭示它的部分信息外，其它的实验手段很难用于这个领域的研究；在理论方面，化学凝胶的理论已经比较成熟，而物理凝胶的粘弹性质以及凝胶化是一个远离平衡态的松弛过程，除了一些特征的标度指数外，人们还没有得到适用于高分子物理凝胶的普适规律。当前，由于计算机模拟理论及模拟方法的发展，使得计算机模拟成为除了实验和理论研究方法之外的第三个重要的研究方法。但是，由于物理凝胶化行为的复杂性，用实验和理论获得的信息很难较好地描述凝胶化过程，而计算机模拟的高度透明性及反映信息的完整性，有助于理解这一复杂过程中所涉及的物理本质。因此，利用计算机模拟结合实验及理论方法深入研究高分子物理凝胶的形成机制、结构与性能关系已成为目前最有效的手段之一。 本论文主要运用Monte Carlo模拟方法，并结合小角中子散射（Small-Angle Neutron Scattering, SANS）和流变（Rheology）等实验手段从多个角度探讨了以下几类典型的高分子溶液物理凝胶化行为。 1. 温度对遥爪型三嵌段共聚物在选择性溶剂中的自组装及凝胶化行为影响的研究：采用二维简单方格子Monte Carlo模拟方法，结合逾渗（Percolation）理论，建立了溶胶－凝胶转变相图在统计热力学中的确定方法；甄别了具有特征构象的链，讨论了链及胶束的聚集，明晰了相互作用（体现为约化温度）、构象转变、聚集与凝胶化的一致的关联关系；提出了构象转变模型，进而明确了此体系的凝胶化过程，在微观尺度上表现为桥型链和环型链之间的竞争。 2. 模拟模型改进及其应用到持续长度对稀溶液中高分子链构象影响的研究：考虑到原始八位置键涨落模型效率低，实现复杂且不能应用到复杂的高分子体系，对该模型进行了改进，使其实现简单、效率高，并拓宽了该模型的应用范围。然后，以刚性对均聚物构象的影响为例，发现随着刚性增加，均聚物构象从球形椭球到棒状椭球的转变，并对比了自由连接链（Free Joint Chain, FJC）模型和蠕虫链（Wormlike Chain, WLC）模型在不同刚性范围内对高分子链末端距预测的偏差，首次给出了这两个经典模型的半定量的适用边界。 3. 溶剂尺寸对遥爪型三嵌段共聚物在选择性溶剂中的自组装及凝胶化行为影响的研究：用改进后的八位置键涨落Monte Carlo模型，研究了遥爪型三嵌段共聚物在选择性溶剂条件下的聚集和凝胶化对溶剂尺寸的依赖性，发现溶剂尺寸效应对凝胶化的作用是非单调的。由一个均聚物体系的对比模拟证明这种作用主要是由熵驱动的，并给出了中分子溶剂的半定量定义。在均聚物和嵌段共聚物溶液中，不同尺寸的溶剂分子可以使溶液由于高分子聚集不同而具有不同的微结构，并影响高分子链构象和溶液的性质。从多个角度研究了三嵌段共聚物在不同尺寸溶剂的溶液中所遵循的三种不同的凝胶化机理。 4. 聚氧化乙烯－氧化丙稀－氧化乙烯三嵌段共聚物（poly(ethylene oxide)-poly (propylene oxide)-poly-(ethylene oxide), PEO-PPO-PEO）重水溶液凝胶化的小角中子散射（SANS）和Monte Carlo研究：结合Pluronic F127（EO65PO99EO65）/D2O三嵌段共聚物溶液的特征，对照SANS数据，用改进后的八位置键涨落模型成功地从模拟中获得了F127/D2O的溶胶－凝胶转变相图。详细地考察了体系的微观结构，提出此类高分子溶液中形成的物理凝胶包含高分子逾渗网络的生成，以及被束缚溶剂（Bound Solvent）必须超过离散组分体系逾渗的临界体积分数的机理。着重研究了一定浓度的F127水溶液随温度升高引起的溶胶－凝胶转变以及凝胶－溶胶转变的Reentrant相行为，发现体系在低温区域的溶胶－凝胶转变遵循相同的机理，而在中等温度和较高温度以及不同浓度区域中的凝胶－溶胶转变遵循不同的机理。 5. 极性基团饱和度和溶剂条件对两亲性聚合物在溶液中的聚集行为和凝胶化影响的研究：用改进后的八位置键涨落模型，针对两亲性聚合物在不同溶剂条件的溶液建立了粗粒化模型，以两亲性聚合物中极性基团的饱和度，溶剂条件和高分子浓度为变量，考察了其对链构象、聚集及其凝胶化的影响。 6. 多糖水溶液凝胶化的流变和小角中子散射研究：用流变和SANS考察了两个多糖水溶液中物理凝胶化过程，针对由氢键主导的水基凝胶体系的典型特征进行了讨论，从分子链构象，聚集体结构及其关联以及流变特征等方面对聚强电解质角叉胶（Carrageenan）水溶液和聚弱电解质明胶（Pectin）水溶液进行了详细的讨论。考察了不同多糖的种类（聚合物链的电荷密度），盐的种类和浓度，溶液温度等对凝胶化和凝胶结构的影响，分析了不同多糖溶液的凝胶化机理。

Rigidity effect on phase behavior of symmetric ABA triblock copolymers: A Monte Carlo simulation

The phase behavior of symmetric ABA triblock copolymers containing a semiflexible midblock is studied by lattice Monte Carlo simulation. As the midblock evolves from a fully flexible state to a semiflexible state in terms of increase in its persistence length, different phase behaviors are observed while cooling the system from an infinite high temperature to a temperature below T-ODT (order-disorder transition temperature). Within the midblock flexibility range we studied (l(p)/N-c <= 0.105), a lamellar structure is formed at equilibrium state as the situation for fully flexible chains. The fraction of bridge chain is evaluated for the lamellar structures. We find that the increase in midblock rigidity indeed results in the increase in bridge chain fraction within the range from 44.9% to 51.8%.

Monte Carlo simulation on symmetric ABA/AB copolymer blends in confined thin films

The effects of blend composition on morphology, order-disorder transition (ODT), and chain conformation of symmetric ABA/AB copolymer blends confined between two neutral hard walls have been investigated by lattice Monte Carlo simulation. Only lamellar structure is observed in all the simulation morphologies under thermodynamic equilibrium state, which is supported by theoretical prediction. When the composition of AB diblock copolymer (phi) increases, both lamellar spacing and the corresponding ODT temperature increase, which can be attributed to the variation of conformation distribution of the diblock and the triblock copolymer chains. In addition, both diblock and triblock copolymer, chains with bridge conformation extend dramatically in the direction parallel to the surface when the system is in ordered state. Finally, the copolymer chain conformation depends strongly on both the blend composition and the incompatibility parameter chi N.

Monte Carlo simulation of the compatibility of graft copolymer compatibilized two incompatible homopolymer blends: Effect of graft structure

Compatibility of graft copolymer compatibilized two incompatible homopolymer A and B blends was simulated by using Monte Carlo method in a two-dimensional lattice model. The copolymers with various graft structures were introduced in order to study the effect of graft structure on the compatibility. Simulation results showed that incorporation of both A-g-B (A was backbone) and B-g-A (B was backbone) copolymers could much improve the compatibility of the blends. However, A-g-B copolymer was more effective to compatibilize the blend if homopolymer A formed dispersed phase. Furthermore, simulation results indicated that A-g-B copolymers tended to locate at the interface and anchor two immiscible components when the side chain is relatively long. However, most of A-g-B copolymers were likely to be dispersed into the dispersed homopolymer A phase domains if the side chains were relatively short. On the other hand, B-g-A copolymers tended to be dispersed into the matrix formed by homopolymer B. Moreover, it was found that more and more B-g-A copolymers were likely to form thin layers at the phase interface with decreasing the length of side chain.

Grafting of maleic anhydride onto linear polyethylene: A Monte Carlo study

Monte Carlo simulation was used to study the graft of maleic anhydride (MAH) onto linear polyethylene (PE-g-MAH) initiated by dicumyl peroxide (DCP). Simulation results revealed that major MAH monomers attached onto PE chains as branched graft at higher MAH content. However, at extremely low MAH content, the fraction of bridged graft was very close to that of branched graft. This conclusion was somewhat different from the conventional viewpoint, namely, the fraction of bridged graft was always much lower than that of branched graft under any condition. Moreover, the results indicated that the grafting degree increased almost linearly to MAH and DCP concentrations. On the other hand, it was found that the amount of grafted MAH dropped sharply with increasing the length of grafted MAH, indicating that MAH monomers were mainly attached onto the PE chain as single MAH groups or very short oligomers. With respect to the crosslink of PE, the results showed that the fraction of PE-(MAH)(n)-PE crosslink structure increased continuously, and hence the fraction of PE-PE crosslink decreased with increasing MAH concentration.

A Monte Carlo simulation for the micellization of ABC 3-miktoarm star terpolymers in a selective solvent

We have used Monte Carlo simulation to study the micellization of ABC 3-miktoarm star terpolymers in a selective solvent (good to A segment, bad to B and C segments). The simulation results reveal that the self-assembled morphology is determined by the block length, molecular architecture, terpolymer concentration and insolubility of insoluble block in the solvent. In dilute solution, symmetric terpolymers (N-B = N-C = 30) tend to aggregate into a novel wormlike pearl-necklace structure linked by an alternating arrangement of B and C spheres, whereas the asymmetric terpolymers (NB = 10, NC = 50) are likely to aggregate into spherical or cylindrical micelles (formed by C blocks) connected with some small B spheres, when the concentration of terpolymer is relatively low (chain number is 100). However, when the concentration of terpolymer is relatively high (chain number is 250), the symmetric terpolymers tend to aggregate into a netlike structure linked by an alternation of B and C spheres, whereas the asymmetric terpolymers are likely to aggregate into wormlike micelles (formed by C blocks) connected with some of small spheres (formed by B blocks). Moreover, when the insolubility of insoluble block in the solvent is weak, the insoluble blocks aggregate into some incompact micelles.

Monte Carlo simulation of the grafting of maleic anhydride onto polypropylene at higher temperature

The graft of maleic anhydride (MAH) onto isotactic polypropylene (iPP) initiated by dicumyl peroxide (DCP) at 190 degreesC was studied by means of the Monte Carlo method. The ceiling temperature theory, i.e., no possibility for the homopolymerization of MA-H to occur at higher temperatures, was used in this study. The simulation results show that most MAH monomers were grafted onto the radical chain ends arising from beta scission at a lower MAH concentration, whereas the amount of MAH monomers attached to the tertiary carbons was much larger than that grafted onto the radical chain ends at a higher MAH concentration for various DCP concentrations. This conclusion gives a good interpretation for the disagreement on the grafting sites along a PP chain. Moreover, it was found that the grafting degree increased considerably up to a peak value; thereafter, it decreased continuously with increasing MA-H concentration. The peak shifted in the lower MAH concentration direction and became lower and lower with increasing DCP concentration. When the DCP concentration was below 0.1 wt %, the peak was hardly observed. Those results are in good agreement with the experiments.

Monte Carlo simulation of pulsed-laser polymerization in emulsion

The pulsed-laser polymerization in emulsions has been simulated by the Monte Carlo method. Our simulation shows that the best measure of the propagation rate coefficients K-p is the peak maximum of molecular weight distribution for microemulsions when the droplets are small. However, the inflection point at the low-molecular-weight side of the peaks provides the best measure of K-p of bigger droplets. (C) 2000 Elsevier Science Ltd. All rights reserved.

Monte Carlo simulation for the modification of polymer via grafting

Monte Carlo method has been applied to investigate the kinetic of grafting reaction in free radical copolymerization. The simulation is quits in agreement with that of theoretical and experimental results. It proves that the Monte Carlo simulation is an effective method for investigating the grafting reaction of free radical copolymerization. (C) 2000 Elsevier Science Ltd. All rights reserved.

Monte Carlo simulation of phase behavior of polymer blends with special interactions

Full Paper: The phase, behavior of A-B-random copolymer/C-homopolymer, blends with special interaction was studied by a. Monte, Carlo simulation in two dimensions. The interaction between I segment A and segment C was repulsive, whereas it was attractive between segment B and segment C. The simulation results showed that the blend became two large co-continuous phase domains at lower segment-B component compositions, indicating that the blend showed spinodal decomposition. With an increase of the segment-B component, the miscibility between the copolymer,and the polymer was gradually improved up to being miscible. In addition, it was found that segment B tended to move to the surface of the copolymer phase in the case of a lower component of segment B. On the other hand, if was observed that the average, end-to-end distances ((h) over bar) for both copolymer and polymer changed slowly with increasing segment-B component of the copolymer up to 40%, thereafter they increased considerably with increasing segment B component. Moreover, it was found that the (h) over bar of the copolymer was obviously shorter than that of the homopolymer for the segment-B composition, region from 0% to 80%. Finally, a, phase diagram showing I phase and - II phase regions under the condition of constant-temperature is presented.

Monte-Carlo simulation of surface crack growth rate for offshore structural steel E36-Z35

The Monte- Carlo method is used to simulate the surface fatigue crack growth rate for offshore structural steel E36-Z35, and to determine the distributions and relevance of the parameters in the Paris equation. By this method, the time and cost of fatigue crack propagation testing can be reduced. The application of the method is demonstrated by use of four sets of fatigue crack propagation data for offshore structural steel E36-Z35. A comparison of the test data with the theoretical prediction for surface crack growth rate shows the application of the simulation method to the fatigue crack propagation tests is successful.