BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Hsiao-Ping Hsu

Hsiao-Ping Hsu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.softcond-mat.stat-mechphysics.chem-ph

Catalog footprint

What is connected

17works
3topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2016arXiv

Static and dynamic properties of large polymer melts in equilibrium

We present a detailed study of the static and dynamic behavior of long semiflexible polymer chains in a melt. Starting from previously obtained fully equilibrated high molecular weight polymer melts [{\it Zhang et al.} ACS Macro Lett. 3, 198 (2014)] we investigate their static and dynamic scaling behavior as predicted by theory. We find that for semiflexible chains in a melt, results of the mean square internal distance, the probability distributions of the end-to-end distance, and the chain structure factor are well described by theoretical predictions for ideal chains. We examine the motion of monomers and chains by molecular dynamics simulations using the ESPResSo++ package. The scaling predictions of the mean squared displacement of inner monomers, center of mass, and relations between them based on the Rouse and the reptation theory are verified, and related characteristic relaxation times are determined. Finally we give evidence that the entanglement length $N_{e,PPA}$ as determined by a primitive path analysis (PPA) predicts a plateau modulus, $G_N^0=\frac{4}{5}(ρk_BT/N_e)$, consistent with stresses obtained from the Green-Kubo relation. These comprehensively characterized equilibrium structures, which offer a good compromise between flexibility, small $N_e$, computational efficiency, and small deviations from ideality provide ideal starting states for future non-equilibrium studies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Lattice Monte Carlo Simulations of Polymer Melts

We use Monte Carlo simulations to study polymer melts consisting of fully flexible and moderately stiff chains in the bond fluctuation model at a volume fraction $0.5$. In order to reduce the local density fluctuations, we test a pre-packing process for the preparation of the initial configurations of the polymer melts, before the excluded volume interaction is switched on completely. This process leads to a significantly faster decrease of the number of overlapping monomers on the lattice. This is useful for simulating very large systems, where the statistical properties of the model with a marginally incomplete elimination of excluded volume violations are the same as those of the model with strictly excluded volume. We find that the internal mean square end-to-end distance for moderately stiff chains in a melt can be very well described by a freely rotating chain model with a precise estimate of the bond-bond orientational correlation between two successive bond vectors in equilibrium. The plot of the probability distributions of the reduced end-to-end distance of chains of different stiffness also shows that the data collapse is excellent and described very well by the Gaussian distribution for ideal chains. However, while our results confirm the systematic deviations between Gaussian statistics for the chain structure factor $S_c(q)$ [minimum in the Kratky-plot] found by Wittmer et al.~\{EPL {\bf 77} 56003 (2007).\} for fully flexible chains in a melt, we show that for the available chain length these deviations are no longer visible, when the chain stiffness is included. The mean square bond length and the compressibility estimated from collective structure factors depend slightly on the stiffness of the chains.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Monte Carlo Simulations of Lattice Models for Single Polymer Systems

Single linear polymer chains in dilute solutions under good solvent conditions are studied by Monte Carlo simulations with the pruned-enriched Rosenbluth method up to the chain length $N \sim {\cal O}(10^4)$. Based on the standard simple cubic lattice model (SCLM) with fixed bond length and the bond fluctuation model (BFM) with bond lengths in a range between $2$ and $\sqrt{10}$, we investigate the conformations of polymer chains described by self-avoiding walks (SAWs) on the simple cubic lattice, and by random walks (RWs) and non-reversible random walks (NRRWs) in the absence of excluded volume (EV) interactions. In addition to flexible chains, we also extend our study to semiflexible chains for different stiffness controlled by a bending potential. The persistence lengths of chains extracted from the orientational correlations are estimated for all cases. We show that chains based on the BFM are more flexible than those based on the SCLM for a fixed bending energy. The microscopic differences between these two lattice models are discussed and the theoretical predictions of scaling laws given in the literature are checked and verified. Our simulations clarify that a different mapping ratio between the coarse-grained models and the atomistically realistic description of polymers is required in a coarse-graining approach due to the different crossovers to the asymptotic behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Estimation of Persistence Lengths of Semiflexible Polymers: Insight from Simulations

The persistence length of macromolecules is one of their basic characteristics, describing their intrinsic local stiffness. However, it is difficult to extract this length from physical properties of the polymers, different recipes may give answers that disagree with each other. Monte Carlo simulations are used to elucidate this problem, giving a comparative discussion of two lattice models, the self-avoiding walk model extended by a bond bending energy, and bottle-brush polymers described by the bond fluctuation model. The conditions are discussed under which a description of such macromolecules by Kratky-Porod worm-like chains holds, and the question to what extent the persistence length depends on external conditions (such as solvent quality) is considered. The scattering function of semiflexible polymers is discussed in detail, a comparison to various analytic treatments is given, and an outlook to experimental work is presented.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Scattering function of semiflexible polymer chains under good solvent conditions

Using the pruned-enriched Rosenbluth Monte Carlo algorithm, the scattering functions of semiflexible macromolecules in dilute solution under good solvent conditions are estimated both in $d=2$ and $d=3$ dimensions, considering also the effect of stretching forces. Using self-avoiding walks of up to $N = 25600$ steps on the square and simple cubic lattices, variable chain stiffness is modeled by introducing an energy penalty $ε_b$ for chain bending; varying $q_b=\exp (- ε_b/k_BT)$ from $q_b=1$ (completely flexible chains) to $q_b = 0.005$, the persistence length can be varied over two orders of magnitude. For unstretched semiflexible chains we test the applicability of the Kratky-Porod worm-like chain model to describe the scattering function, and discuss methods for extracting persistence length estimates from scattering. While in $d=2$ the direct crossover from rod-like chains to self-avoiding walks invalidates the Kratky-Porod description, it holds in $d=3$ for stiff chains if the number of Kuhn segments $n_K$ does not exceed a limiting value $n^*_K$ (which depends on the persistence length). For stretched chains, the Pincus blob size enters as a further characteristic length scale. The anisotropy of the scattering is well described by the modified Debye function, if the actual observed chain extension $<X>$ (end-to-end distance in the direction of the force) as well as the corresponding longitudinal and transverse linear dimensions $<X^2> - <X>^2$, $<R_{g,\bot}^2>$ are used.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A fast Monte Carlo algorithm for studying bottle-brush polymers

Obtaining reliable estimates of the statistical properties of complex macromolecules by computer simulation is a task that requires high computational effort as well as the development of highly efficient simulation algorithms. We present here an algorithm combining local moves, the pivot algorithm, and an adjustable simulation lattice box for simulating dilute systems of bottle-brush polymers with a flexible backbone and flexible side chains under good solvent conditions. Applying this algorithm to the bond fluctuation model, very precise estimates of the mean square end-to-end distances and gyration radii of the backbone and side chains are obtained, and the conformational properties of such a complex macromolecule are studied. Varying the backbone length (from $N_b=67$ to $N_b=1027$), side chain length (from N=0 to N=24 or 48), the scaling predictions for the backbone behavior as well as the side chain behavior are checked. We are also able to give a direct comparison of the structure factor between experimental data and the simulation results.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A review of Monte Carlo simulations of polymers with PERM

In this review, we describe applications of the pruned-enriched Rosenbluth method (PERM), a sequential Monte Carlo algorithm with resampling, to various problems in polymer physics. PERM produces samples according to any given prescribed weight distribution, by growing configurations step by step with controlled bias, and correcting "bad" configurations by "population control". The latter is implemented, in contrast to other population based algorithms like e.g. genetic algorithms, by depth-first recursion which avoids storing all members of the population at the same time in computer memory. The problems we discuss all concern single polymers (with one exception), but under various conditions: Homopolymers in good solvents and at the $Θ$ point, semi-stiff polymers, polymers in confining geometries, stretched polymers undergoing a forced globule-linear transition, star polymers, bottle brushes, lattice animals as a model for randomly branched polymers, DNA melting, and finally -- as the only system at low temperatures, lattice heteropolymers as simple models for protein folding. PERM is for some of these problems the method of choice, but it can also fail. We discuss how to recognize when a result is reliable, and we discuss also some types of bias that can be crucial in guiding the growth into the right directions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Breakdown of the Kratky-Porod Wormlike Chain Model for Semiflexible Polymers in Two Dimensions

By large-scale Monte Carlo simulations of semiflexible polymers in $d=2$ dimensions the applicability of the Kratky-Porod model is tested. This model is widely used as "standard model" for describing conformations and force versus extension curves of stiff polymers. It is shown that semiflexible polymers in $d=2$ show a crossover from hard rods to self-avoiding walks, the intermediate Gaussian regime (implied by the Kratky-Porod model) is completely absent. Hence the latter can also describe force versus extension curves only if the contour length is only a few times larger than the persistence length. Consequences for experiments on biopolymers at interfaces are briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Computer simulation of bottle brush polymers with flexible backbone: Good solvent versus Theta solvent conditions

By Molecular Dynamics simulation of a coarse-grained bead-spring type model for a cylindrical molecular brush with a backbone chain of $N_b$ effective monomers to which with grafting density $σ$ side chains with $N$ effective monomers are tethered, several characteristic length scales are studied for variable solvent quality. Side chain lengths are in the range $5 \le N \le 40$, backbone chain lengths are in the range $50 \le N_b \le 200$, and we perform a comparison to results for the bond fluctuation model on the simple cubic lattice (for which much longer chains are accessible, $N_b \le 1027$, and which corresponds to an athermal, very good, solvent). We obtain linear dimensions of side chains and the backbone chain and discuss their $N$-dependence in terms of power laws and the associated effective exponents. We show that even at the Theta point the side chains are considerably stretched, their linear dimension depending on the solvent quality only weakly. Effective persistence lengths are extracted both from the orientational correlations and from the backbone end-to-end distance; it is shown that different measures of the persistence length (which would all agree for Gaussian chains) are not mutually consistent with each other, and depend distinctly both on $N_b$ and the solvent quality. A brief discussion of pertinent experiments is given.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

New Development of Monte Carlo Techniques for Studying Bottle-brush Polymers

Due to the complex characteristics of bottle-brush polymers, it became a challenge to develop an efficient algorithm for studying such macromolecules under various solvent conditions or some constraints in the space by using computer simulations. In the limit of a bottle-brush polymer with a rather stiff backbone (straight rigid backbone), we generalize the variant of the biased chain growth algorithm, the pruned-enriched Rosenbluth method, for simulating polymers with complex architecture, from star polymers to bottle-brush polymers, on the simple cubic lattice. With the high statistics of our Monte Carlo results, we check the theoretical predictions of side chain behavior and radial monomer density profile. For the comparison of the experimental data for bottle-brush polymers with a flexible backbone and flexible side chains, based on the bond fluctuation model we propose another fast Monte Carlo algorithm combining the local moves, the pivot move, and an adjustable simulation lattice box. By monitoring the autocorrelation functions of gyration radii for the side chains and for the backbone, we see that for fixed side chain length there is no change in the behavior of these two functions as the backbone length increases. Our extensive results cover the range which is accessible for the comparison to experimental data and for the checking of the theoretically predicted scaling laws.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Scaling behaviour of lattice animals at the upper critical dimension

We perform numerical simulations of the lattice-animal problem at the upper critical dimension d=8 on hypercubic lattices in order to investigate logarithmic corrections to scaling there. Our stochastic sampling method is based on the pruned-enriched Rosenbluth method (PERM), appropriate to linear polymers, and yields high statistics with animals comprised of up to 8000 sites. We estimate both the partition sums (number of different animals) and the radii of gyration. We re-verify the Parisi-Sourlas prediction for the leading exponents and compare the logarithmic-correction exponents to two partially differing sets of predictions from the literature. Finally, we propose, and test, a new Parisi-Sourlas-type scaling relation appropriate for the logarithmic-correction exponents.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Stretching Semiflexible Polymer Chains: Evidence for the Importance of Excluded Volume Effects from Monte Carlo Simulation

Semiflexible macromolecules in dilute solution under very good solvent conditions are modeled by self-avoiding walks on the simple cubic lattice ($d=3$ dimensions) and square lattice ($d=2$ dimensions), varying chain stiffness by an energy penalty $ε_b$ for chain bending. In the absence of excluded volume interactions, the persistence length $\ell_p$ of the polymers would then simply be $\ell_p=\ell_b(2d-2)^{-1}q_b^{-1}$ with $q_b= \exp(-ε_b/k_BT)$, the bond length $\ell_b$ being the lattice spacing, and $k_BT$ is the thermal energy. Using Monte Carlo simulations applying the pruned-enriched Rosenbluth method (PERM), both $q_b$ and the chain length $N$ are varied over a wide range $(0.005 \leq q_b \leq 1, \; N \leq 50000$), and also a stretching force $f$ is applied to one chain end (fixing the other end at the origin). In the absence of this force, in $d=2$ a single crossover from rod-like behavior (for contour lengths less than $\ell_p$) to swollen coils occurs, invalidating the Kratky-Porod model, while in $d=3$ a double crossover occurs, from rods to Gaussian coils (as implied by the Kratky-Porod model) and then to coils that are swollen due to the excluded volume interaction. If the stretching force is applied, excluded volume interactions matter for the force versus extension relation irrespective of chain stiffness in $d=2$, while theories based on the Kratky-Porod model are found to work in $d=3$ for stiff chains in an intermediate regime of chain extensions. While for $q_b \ll 1$ in this model a persistence length can be estimated from the initial decay of bond-orientational correlations, it is argued that this is not possible for more complex wormlike chains (e.g. bottle-brush polymers). Consequences for the proper interpretation of experiments are briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Conformational Studies of bottle-brush polymers absorbed on a flat solid surface

The adsorption of a bottle-brush polymer end-grafted with one chain end of its backbone to a flat substrate surface is studied by Monte Carlo simulation of a coarse-grained model, that previously has been characterized in the bulk, assuming a dilute solution under good solvent conditions. Applying the bond fluctuation model on the simple cubic lattice, we vary the backbone chain length $N_b$ from $N_b=67$ to $N_b = 259$ effective monomeric units, the side chain length $N$ from N=6 to N=48, and the grafting density $σ=1$, i.e., parameters that correspond well to the experimentally accessible range. When the adsorption energy strength $ε$ is varied, we find that the adsorption transition (which becomes well-defined in the limit $N_b \rightarrow \infty$, for arbitrary finite $N$) roughly occurs at the same value $ε_c$ as for ordinary linear chains (N=0), at least within our statistical errors. Mean square end-to-end distances and gyration radii of the side chains are obtained, as well as the monomer density profile in the direction perpendicular to the adsorbing surface. We show that for longer side chains the adsorption of bottle-brushes is a two step process, the decrease of the perpendicular linear dimension of side chains with adsorption energy strength can even be non-monotonic. Also the behavior of the static structure factor $S(q)$ is analyzed, evidence for a quasi-two-dimensional scaling is presented, and consequences for the interpretation of experiments are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Polymer chain stiffness versus excluded volume: A Monte Carlo study of the crossover towards the wormlike chain model

When the local intrinsic stiffness of a polymer chain varies over a wide range, one can observe both a crossover from rigid-rod-like behavior to (almost) Gaussian random coils and a further crossover towards self-avoiding walks in good solvents. Using the pruned-enriched Rosenbluth method (PERM) to study self-avoiding walks of up to $N_b=50000$ steps and variable flexibility, the applicability of the Kratky-Porod model is tested. Evidence for non-exponential decay of the bond-orientational correlations $<\cos θ(s) >$ for large distances $s$ along the chain contour is presented, irrespective of chain stiffness. For bottle-brush polymers on the other hand, where experimentally stiffness is varied via the length of side-chains, it is shown that these cylindrical brushes (with flexible backbones) are not described by the Kratky-Porod wormlike chain model, since their persistence length is (roughly) proportional to their cross-sectional radius, for all conditions of practical interest.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

How to Define Variation of Physical Properties Normal to an Undulating One-Dimensional Object

One-dimensional flexible objects are abundant in physics, from polymers to vortex lines to defect lines and many more. These objects structure their environment and it is natural to assume that the influence these objects exert on their environment depends on the distance from the line-object. But how should this be defined? We argue here that there is an intrinsic length scale along the undulating line that is a measure of its "stiffness" (i.e., orientational persistence), which yields a natural way of defining the variation of physical properties normal to the undulating line. We exemplify how this normal variation can be determined from a computer simulation for the case of a so-called bottle-brush polymer, where side chains are grafted onto a flexible backbone.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Polymers grafted to porous membranes

We study a single flexible chain molecule grafted to a membrane which has pores of size slightly larger than the monomer size. On both sides of the membrane there is the same solvent. When this solvent is good, i.e. when the polymer is described by a self avoiding walk, it can fairly easily penetrate the membrane, so that the average number of membrane crossings tends, for chain length $N\to\infty$, to a positive constant. The average numbers of monomers on either side of the membrane diverges in this limit, although their ratio becomes infinite. For a poor solvent, in contrast, the entire polymer is located, for large $N$, on one side of the membrane. For good and for theta solvents (ideal polymers) we find scaling laws, whose exponents can in the latter case be easily understood from the behaviour of random walks.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Universal scaling functions for bond percolation on planar random and square lattices with multiple percolating clusters

Percolation models with multiple percolating clusters have attracted much attention in recent years. Here we use Monte Carlo simulations to study bond percolation on $L_{1}\times L_{2}$ planar random lattices, duals of random lattices, and square lattices with free and periodic boundary conditions, in vertical and horizontal directions, respectively, and with various aspect ratio $L_{1}/L_{2}$. We calculate the probability for the appearance of $n$ percolating clusters, $W_{n},$ the percolating probabilities, $P$, the average fraction of lattice bonds (sites) in the percolating clusters, $<c^{b}>_{n}$ ($<c^{s}>_{n}$), and the probability distribution function for the fraction $c$ of lattice bonds (sites), in percolating clusters of subgraphs with $n$ percolating clusters, $f_{n}(c^{b})$ ($f_{n}(c^{s})$). Using a small number of nonuniversal metric factors, we find that $W_{n}$, $P$, $<c^{b}>_{n}$ ($<c^{s}>_{n}$), and $f_{n}(c^{b})$ ($f_{n}(c^{s})$) for random lattices, duals of random lattices, and square lattices have the same universal finite-size scaling functions. We also find that nonuniversal metric factors are independent of boundary conditions and aspect ratios.

0 citations0 reviewsNo code linkedOpen access