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

Kurt Binder

Kurt Binder 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-mechBiological Physicscond-mat.mtrl-sciphysics.comp-phcond-mat.dis-nnphysics.chem-ph

Catalog footprint

What is connected

46works
7topics
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

46 published item(s)

preprint2021arXiv

Wetting transitions of polymer solutions: Effects of chain length and chain stiffness

Wetting and drying phenomena are studied for flexible and semiflexible polymer solutions via coarse-grained molecular dynamics simulations and density functional theory calculations. The study is based on the use of Young's equation for the contact angle, determining all relevant surface tensions from the anisotropy of the pressure tensor. The solvent quality (or effective temperature, equivalently) is varied systematically, while all other interactions remain unaltered. For flexible polymers, the wetting transition temperature $T_{\rm w}$ increases monotonically with chain length $N$, while the contact angle at temperatures far below $T_{\rm w}$ is independent of $N$. For semiflexible polymer solutions, $T_{\rm w}$ varies non-monotonically with the persistence length: Initially, $T_{\rm w}$ increases with increasing chain stiffness and reaches a maximum, but then a sudden drop of $T_{\rm w}$ is observed, which is associated with the isotropic-nematic transition of the system.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Entropic Unmixing in Nematic Blends of Semiflexible Polymers

Binary mixtures of semiflexible polymers with the same chain length but different persistence lengths separate into two coexisting different nematic phases when the osmotic pressure of the lyotropic solution is varied. Molecular Dynamics simulations and Density Functional Theory predict phase diagrams either with a triple point, where the isotropic phase coexists with two nematic phases, or a critical point of unmixing within the nematic mixture. The difference in locally preferred bond angles between the constituents drives this unmixing without any attractive interactions between monomers.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

When does Wenzel's extension of Young's equation for the contact angle of droplets apply? A density functional study

he contact angle of a liquid droplet on a surface under partial wetting conditions differs for a nanoscopically rough or periodically corrugated surface from its value for a perfectly flat surface. Wenzel's relation attributes this difference simply to the geometric magnification of the surface area (by a factor $r_{\rm w}$), but the validity of this idea is controversial. We elucidate this problem by model calculations for a sinusoidal corrugation of the form $z_{\rm wall}(y) = Δ\cos(2πy/λ)$ , for a potential of short range $σ_{\rm w}$ acting from the wall on the fluid particles. When the vapor phase is an ideal gas, the change of the wall-vapor surface tension can be computed exactly, and corrections to Wenzel's equation are typically of order $σ_{\rm w}Δ/λ^2$. For fixed $r_{\rm w}$ and fixed $σ_{\rm w}$ the approach to Wenzel's result with increasing $λ$ may be nonmonotonic and this limit often is only reached for $λ/σ_{\rm w}>30$. For a non-additive binary mixture, density functional theory is used to work out the density profiles of both coexisting phases both for planar and corrugated walls, as well as the corresponding surface tensions. Again, deviations from Wenzel's results of similar magnitude as in the above ideal gas case are predicted. Finally, a crudely simplified description based on the interface Hamiltonian concept is used to interpret corresponding simulation results along similar lines. Wenzel's approach is found to generally hold when $λ/σ_{\rm w}\gg 1$, $Δ/λ<1$, and conditions avoiding proximity of wetting or filling transitions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Anomalous fluctuations of nematic order in solutions of semiflexible polymers

The nematic ordering in semiflexible polymers with contour length $L$ exceeding their persistence length $\ell_p$ is described by a confinement of the polymers in a cylinder of radius $r_{eff}$ much larger than the radius $r_ρ$, expected from the respective concentration of the solution. Large scale Molecular Dynamics simulations combined with Density Functional Theory are used to locate the Isotropic-Nematic ($I-N$)-transition and to validate this cylindrical confinement. Anomalous fluctuations, due to chain deflections from neighboring chains in the nematic phase are proposed. Considering deflections as collective excitations in the nematically ordered phase of semiflexible polymers elucidates the origins of shortcomings in the description of the $I-N$ transition by existing theories.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Semiflexible Polymers Under Good Solvent Conditions Interacting With Repulsive Walls

Solutions of semiflexible polymers confined by repulsive planar walls are studied by density functional theory and Molecular Dynamics simulations, to clarify the competition between the chain alignment favored by the wall and the depletion caused by the monomer-wall repulsion. A coarse-grained bead-spring model with bond bending potential is studied, varying both the contour length and the persistence length of the polymers, as well as the monomer concentration in the solution (good solvent conditions are assumed throughout, and solvent molecules are not included explicitly). The profiles of monomer density and pressure tensor components near the wall are studied, and the surface tension of the solution is obtained. While the surface tension slightly decreases with chain length for flexible polymers, it clearly increases with chain length for stiff polymers. Thus, at fixed density and fixed chain length the surface tension also increases with increasing persistence length. Chain ends always are enriched near the wall, but this effect is much larger for stiff polymers than for flexible ones. Also the profiles of the mean square gyration radius components near the wall and the nematic order parameter are studied to clarify the conditions where wall-induced nematic order occurs.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Crystal nuclei in melts: A Monte Carlo simulation of a model for attractive colloids

As a model for a suspension of hard-sphere like colloidal particles where small nonadsorbing dissolved polymers create a depletion attraction, we introduce an effective colloid-colloid potential closely related to the Asakura-Oosawa model but that does not have any discontinuities. In simulations, this model straightforwardly allows the calculation of the pressure from the Virial formula, and the phase transition in the bulk from the liquid to crystalline solid can be accurately located from a study where a stable coexistence of a crystalline slab with a surrounding liquid phase occurs. For this model, crystalline nuclei surrounded by fluid are studied both by identifying the crystal-fluid interface on the particle level (using suitable bond orientational order parameters to distinguish the phases) and by "thermodynamic" means. I.e., the latter method amounts to compute the enhancement of chemical potential and pressure relative to their coexistence values. We show that the chemical potential can be obtained from simulating thick films, where one wall with a rather long range repulsion is present, since near this wall the Widom particle insertion method works, exploiting the fact that the chemical potential in the system is homogeneous. Finally, the surface excess free energy of the nucleus is obtained, for a wide range of nuclei volumes. From this method, it is established that classical nucleation theory works, showing that for the present model the anisotropy of the interface excess free energy of crystals and their resulting nonspherical shape has only a very small effect on the barrier.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Finite size effects on liquid-solid phase coexistence and the estimation of crystal nucleation barriers

A fluid in equilibrium in a finite volume $V$ with particle number $N$ at a density $ρ= N/V$ exceeding the onset density $ρ_f $ of freezing may exhibit phase coexistence between a crystalline nucleus and surrounding fluid. Using a method suitable for the estimation of the chemical potential of dense fluids we obtain the excess free energy due to the surface of the crystalline nucleus. There is neither a need to precisely locate the interface nor to compute the (anisotropic) interfacial tension. As a test case, a soft version of the Asakura-Oosawa model for colloid polymer-mixtures is treated. While our analysis is appropriate for crystal nuclei of arbitrary shape, we find the nucleation barrier to be compatible with a spherical shape, and consistent with classical nucleation theory.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Anisotropic interfacial tension, contact angles, and line tensions: A graphics-processing-unit-based Monte Carlo study of the Ising model

As a generic example for crystals where the crystal-fluid interface tension depends on the orientation of the interface relative to the crystal lattice axes, the nearest neighbor Ising model on the simple cubic lattice is studied over a wide temperature range, both above and below the roughening transition temperature. Using a thin film geometry $L_x \times L_y \times L_z$ with periodic boundary conditions along the z-axis and two free $L_x \times L_y$ surfaces at which opposing surface fields $\pm H_{1}$ act, under conditions of partial wetting, a single planar interface inclined under a contact angle $θ< π/2$ relative to the yz-plane is stabilized. In the y-direction, a generalization of the antiperiodic boundary condition is used that maintains the translational invariance in y-direction despite the inhomogeneity of the magnetization distribution in this system. This geometry allows a simultaneous study of the angle-dependent interface tension, the contact angle, and the line tension (which depends on the contact angle, and on temperature). All these quantities are extracted from suitable thermodynamic integration procedures. In order to keep finite size effects as well as statistical errors small enough, rather large lattice sizes (of the order of 46 million sites) are found necessary, availability of very efficient code implementation of graphics processing units (GPUs) was crucial for the feasibility of this study.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Conformations, Transverse Fluctuations and Crossover Dynamics of a Semi-Flexible Chain in Two Dimensions

We present a unified scaling description for the dynamics of monomers of a semiflexible chain under good solvent condition in the free draining limit. We consider both the cases where the contour length $L$ is comparable to the persistence length $\ell_p$ and the case $L\gg \ell_p$. Our theory captures the early time monomer dynamics of a stiff chain characterized by $t^{3/4}$ dependence for the mean square displacement(MSD) of the monomers, but predicts a first crossover to the Rouse regime of $t^{2ν/{1+2ν}}$ for $τ_1 \sim \ell_p^3$, and a second crossover to the purely diffusive dynamics for the entire chain at $τ_2 \sim L^{5/2}$. We confirm the predictions of this scaling description by studying monomer dynamics of dilute solution of semi-flexible chains under good solvent conditions obtained from our Brownian dynamics (BD) simulation studies for a large choice of chain lengths with number of monomers per chain N = 16 - 2048 and persistence length $\ell_p = 1 - 500$ Lennard-Jones (LJ) units. These BD simulation results further confirm the absence of Gaussian regime for a 2d swollen chain from the slope of the plot of $\langle R_N^2 \rangle/2L \ell_p \sim L/\ell_p$ which around $L/\ell_p \sim 1$ changes suddenly from $\left(L/\ell_p \right) \rightarrow \left(L/\ell_p \right)^{0.5} $, also manifested in the power law decay for the bond autocorrelation function disproving the validity of the WLC in 2d. We further observe that the normalized transverse fluctuations of the semiflexible chains for different stiffness $\sqrt{\langle l_{\bot}^2\rangle}/L$ as a function of renormalized contour length $L/\ell_p$ collapse on the same master plot and exhibits power law scaling $\sqrt{\langle l_{\bot}^2\rangle}/L \sim (L/\ell_p)^η$ at extreme limits, where $η= 0.5$ for extremely stiff chains ($L/\ell_p \gg 1$), and $η= -0.25$ for fully flexible chains.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Determination of the origin and magnitude of logarithmic finite-size effects on interfacial tension: Role of interfacial fluctuations and domain breathing

The ensemble-switch method for computing wall excess free energies of condensed matter is extended to estimate the interface free energies between coexisting phases very accurately. By this method, system geometries with linear dimensions $L$ parallel and $L_z$ perpendicular to the interface with various boundary conditions in the canonical or grandcanonical ensemble can be studied. Using two- and three-dimensional Ising models, the nature of the occurring logarithmic finite size corrections is studied. It is found crucial to include interfacial fluctuations due to "domain breathing".

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Logarithmic Finite-Size Effects on Interfacial Free Energies: Phenomenological Theory and Monte Carlo Studies

The computation of interfacial free energies between coexisting phases (e.g.~saturated vapor and liquid) by computer simulation methods is still a challenging problem due to the difficulty of an atomistic identification of an interface, and due to interfacial fluctuations on all length scales. The approach to estimate the interfacial tension from the free energy excess of a system with interfaces relative to corresponding single-phase systems does not suffer from the first problem but still suffers from the latter. Considering $d$-dimensional systems with interfacial area $L^{d-1}$ and linear dimension $L_z$ in the direction perpendicular to the interface, it is argued that the interfacial fluctuations cause logarithmic finite-size effects of order $\ln (L) / L^{d-1}$ and order $\ln (L_z)/L ^{d-1}$, in addition to regular corrections (with leading order $\text{const}/L^{d-1}$). A phenomenological theory predicts that the prefactors of the logarithmic terms are universal (but depend on the applied boundary conditions and the considered statistical ensemble). The physical origin of these corrections are the translational entropy of the interface as a whole, "domain breathing" (coupling of interfacial fluctuations to the bulk order parameter fluctuations of the coexisting domains), and capillary waves. Using a new variant of the ensemble switch method, interfacial tensions are found from Monte Carlo simulations of $d=2$ and $d=3$ Ising models and a Lennard Jones fluid. The simulation results are fully consistent with the theoretical predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Phase Behavior of Active Swimmers in Depletants: Molecular Dynamics and Integral Equation Theory

We study the structure and phase behavior of a binary mixture where one of the components is self-propelling in nature. The inter-particle interactions in the system were taken from the Asakura-Oosawa model, for colloid-polymer mixtures, for which the phase diagram is known. In the current model version the colloid particles were made active using the Vicsek model for self-propelling particles. The resultant active system was studied by molecular dynamics methods and integral equation theory. Both methods produce results consistent with each other and demonstrate that the Vicsek model based activity facilitates phase separation, thus broadening the coexistence region.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Transitions of tethered chain molecules under tension

An applied tension force changes the equilibrium conformations of a polymer chain tethered to a planar substrate and thus affects the adsorption transition as well as the coil-globule and crystallization transitions. Conversely, solvent quality and surface attraction are reflected in equilibrium force-extension curves that can be measured in experiments. To investigate these effects theoretically, we study tethered chains under tension with Wang-Landau simulations of a bond-fluctuation lattice model. Applying our model to pulling experiments on biological molecules we obtain a good description of experimental data in the intermediate force range, where universal features dominate and finite size effects are small. For tethered chains in poor solvent, we observe the predicted two-phase coexistence at transitions from the globule to stretched conformations and also discover direct transitions from crystalline to stretched conformations. A phase portrait for finite chains constructed by evaluating the density of states for a broad range of solvent conditions and tensions shows how increasing tension leads to a disappearance of the globular phase. For chains in good solvents tethered to hard and attractive surfaces we find the predicted scaling with the chain length in the low-force regime and show that our results are well described by an analytical, independent-bond approximation for the bond-fluctuation model for the highest tensions. Finally, for a hard or slightly attractive surface the stretching of a tethered chain is a conformational change that does not correspond to a phase transition. However, when the surface attraction is sufficient to adsorb a chain it will undergo a desorption transition at a critical value of the applied force. Our results for force-induced desorption show the transition to be discontinuous with partially desorbed conformations in the coexistence region.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Unconventional ordering behavior of semi-flexible polymers in dense brushes under compression

Using a coarse-grained bead-spring model for semi-flexible macromolecules forming a polymer brush, structure and dynamics of the polymers is investigated, varying chain stiffness and grafting density. The anchoring condition for the grafted chains is chosen such that their first bonds are oriented along the normal to the substrate plane. Compression of such a semi-flexible brush by a planar piston is observed to be a two-stage process: for small compressions the chains contract by "buckling" deformation whereas for larger compression the chains exhibit a collective (almost uniform) bending deformation. Thus, the stiff polymer brush undergoes a 2-nd order phase transition of collective bond reorientation. The pressure, required to keep the stiff brush at a given degree of compression, is thereby significantly smaller than for an otherwise identical brush made of entirely flexible polymer chains! While both the brush height and the chain linear dimension in the z-direction perpendicular to the substrate increase monotonically with increasing chain stiffness, lateral (xy) chain linear dimensions exhibit a maximum at intermediate chain stiffness. Increasing the grafting density leads to a strong decrease of these lateral dimensions, compatible with an exponential decay. Also the recovery kinetics after removal of the compressing piston is studied, and found to follow a power-law / exponential decay with time. A simple mean-field theoretical consideration, accounting for the buckling/bending behavior of semi-flexible polymer brushes under compression, is suggested.

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
preprint2013arXiv

Monte Carlo Tests of Nucleation Concepts in the Lattice Gas Model

The conventional theory of homogeneous and heterogeneous nucleation in a supersaturated vapor is tested by Monte Carlo simulations of the lattice gas (Ising) model with nearest-neighbor attractive interactions on the simple cubic lattice. The theory considers the nucleation process as a slow (quasi-static) cluster (droplet) growth over a free energy barrier $ΔF^*$, constructed in terms of a balance of surface and bulk term of a "critical droplet" of radius $R^*$, implying that the rates of droplet growth and shrinking essentially balance each other for droplet radius $R=R^*$. For heterogeneous nucleation at surfaces, the barrier is reduced by a factor depending on the contact angle. Using the definition of "physical" clusters based on the Fortuin-Kasteleyn mapping, the time-dependence of the cluster size distribution is studied for "quenching experiments" in the kinetic Ising model, and the cluster size $\ell ^*$ where the cluster growth rate changes sign is estimated. These studies of nucleation kinetics are compared to studies where the relation between cluster size and supersaturation is estimated from equilibrium simulations of phase coexistence between droplet and vapor in the canonical ensemble. The chemical potential is estimated from a lattice version of the Widom particle insertion method. For large droplets it is shown that the "physical clusters" have a volume consistent with the estimates from the lever rule. "Geometrical clusters" (defined such that each site belonging to the cluster is occupied and has at least one occupied neighbor site) yield valid results only for temperatures less than 60% of the critical temperature, where the cluster shape is non-spherical. We show how the chemical potential can be used to numerically estimate $ΔF^*$ also for non-spherical cluster shapes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Non mean-field behaviour of critical wetting transition for short-range forces

Critical wetting transition for short-range forces in three dimensions ($d=3$) is reinvestigated by means of Monte Carlo simulation. Using anisotropic finite size scaling approach, as well as approaches that do not rely on finite size scaling, we show that the critical wetting transition shows clear deviation from mean-field behaviour. We estimate that the effective critical exponent $ν_{\|}^{\textrm{eff}}=1.76\pm 0.08$ for $J/kT=0.35$ and $ν_{\|}^{\textrm{eff}}=1.85\pm 0.07$ for $J/kT=0.25$. These values are in accord with predictions of Parry {\it et al.} [Phys. Rev. Lett. {\bf 100}, 136105 (2008)]. We also point out that the anisotropic finite size scaling approach in $d=3$ requires additional modification in order to reach full consistency of simulational results.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Universal monomer dynamics of a two dimensional semi-flexible chain

We present a unified scaling theory for the dynamics of monomers for dilute solutions of semiflexible polymers under good solvent conditions in the free draining limit. Our theory encompasses the well-known regimes of mean square displacements (MSDs) of stiff chains growing like t^{3/4} with time due to bending motions, and the Rouse-like regime t^{2 ν/ (1+ 2ν)} where νis the Flory exponent describing the radius R of a swollen flexible coil. We identify how the prefactors of these laws scale with the persistence length l_p, and show that a crossover from stiff to flexible behavior occurs at a MSD of order l^2_p (at a time proportional to l^3_p). A second crossover (to diffusive motion) occurs when the MSD is of order R^2. Large scale Molecular Dynamics simulations of a bead-spring model with a bond bending potential (allowing to vary l_p from 1 to 200 Lennard-Jones units) provide compelling evidence for the theory, in D=2 dimensions where ν=3/4. Our results should be valuable for understanding the dynamics of DNA (and other semiflexible biopolymers) adsorbed on substrates.

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
preprint2012arXiv

Simulation of fluid-solid coexistence in finite volumes: A method to study the properties of wall-attached crystalline nuclei

The Asakura-Oosawa model for colloid-polymer mixtures is studied by Monte Carlo simulations at densities inside the two-phase coexistence region of fluid and solid. Choosing a geometry where the system is confined between two flat walls, and a wall-colloid potential that leads to incomplete wetting of the crystal at the wall, conditions can be created where a single nanoscopic wall-attached crystalline cluster coexists with fluid in the remainder of the simulation box. Following related ideas that have been useful to study heterogeneous nucleation of liquid droplets at the vapor-liquid coexistence, we estimate the contact angles from observations of the crystalline clusters in thermal equilibrium. We find fair agreement with a prediction based on Young's equation, using estimates of interface and wall tension from the study of flat surfaces. It is shown that the pressure versus density curve of the finite system exhibits a loop, but the pressure maximum signifies the "droplet evaporation-condensation" transition and thus has nothing in common with a van der Waals-like loop. Preparing systems where the packing fraction is deep inside the two-phase coexistence region, the system spontaneously forms a "slab state", with two wall-attached crystalline domains separated by (flat) interfaces from liquid in full equilibrium with the crystal in between; analysis of such states allows a precise estimation of the bulk equilibrium properties at phase coexistence.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Transitions between imperfectly ordered crystalline structures: A phase switch Monte Carlo study

A model for two-dimensional colloids confined laterally by "structured boundaries" (i.e., ones that impose a periodicity along the slit) is studied by Monte Carlo simulations. When the distance D between the confining walls is reduced at constant particle number from an initial value D_0, for which a crystalline structure commensurate with the imposed periodicity fits, to smaller values, a succession of phase transitions to imperfectly ordered structures occur. These structures have a reduced number of rows parallel to the boundaries (from n to n-1 to n-2 etc.) and are accompanied by an almost periodic strain pattern, due to "soliton staircases" along the boundaries. Since standard simulation studies of such transitions are hampered by huge hysteresis effects, we apply the phase switch Monte Carlo method to estimate the free energy difference between the structures as a function of the misfit between D and D_0, thereby locating where the transitions occurs in equilibrium. For comparison, we also obtain this free energy difference from a thermodynamic integration method: the results agree, but the effort required to obtain the same accuracy as provided by phase switch Monte Carlo would be at least three orders of magnitude larger. We also show for a situation where several "candidate structures" exist for a phase, that phase switch Monte Carlo can clearly distinguish the metastable structures from the stable one. Finally, applying the method in the conjugate statistical ensemble (where the normal pressure conjugate to D is taken as an independent control variable) we show that the standard equivalence between the conjugate ensembles of statistical mechanics is violated.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Wetting Transition in the Two-Dimensional Blume-Capel Model: A Monte Carlo study

The wetting transition of the Blume-Capel model is studied by a finite-size scaling analysis of $L \times M$ lattices where competing boundary fields $\pm H_1$ act on the first row or last row of the $L$ rows in the strip, respectively. We show that using the appropriate anisotropic version of finite size scaling, critical wetting in $d=2$ is equivalent to a "bulk" critical phenomenon with exponents $α=-1$, $β=0$, and $γ=3$. These concepts are also verified for the Ising model. For the Blume-Capel model it is found that the field strength $H_{1c} (T)$ where critical wetting occurs goes to zero when the bulk second-order transition is approached, while $H_{1c}(T)$ stays nonzero in the region where in the bulk a first-order transition from the ordered phase, with nonzero spontaneous magnetization, to the disordered phase occurs. Interfaces between coexisting phases then show interfacial enrichment of a layer of the disordered phase which exhibits in the second order case a finite thickness only. A tentative discussion of the scaling behavior of the wetting phase diagram near the tricritical point also is given.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Anomalous Structure and Scaling of Ring Polymer Brushes

A comparative simulation study of polymer brushes formed by grafting at a planar surface either flexible linear polymers (chain length $N_L$) or (non-catenated) ring polymers (chain length $N_R=2 N_L$) is presented. Two distinct off-lattice models are studied, one by Monte Carlo methods, the other by Molecular Dynamics, using a fast implementation on graphics processing units (GPUs). It is shown that the monomer density profiles $ρ(z)$ in the $z$-direction perpendicular to the surface for rings and linear chains are practically identical, $ρ_R(2 N_L, z)=ρ_L(N_L, z)$. The same applies to the pressure, exerted on a piston at hight z, as well. While the gyration radii components of rings and chains in $z$-direction coincide, too, and increase linearly with $N_L$, the transverse components differ, even with respect to their scaling properties: $R_{gxy}^{(L)} \propto N_L^{1/2}$, $R_{gxy}^{(R)} \propto N_L^{0.4}$. These properties are interpreted in terms of the statistical properties known for ring polymers in dense melts.

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

Hard sphere fluids confined between soft repulsive walls: A comparative study using Monte Carlo and density functional methods

Hard-sphere fluids confined between parallel plates a distance $D$ apart are studied for a wide range of packing fractions, including also the onset of crystallization, applying Monte Carlo simulation techniques and density functional theory. The walls repel the hard spheres (of diameter $σ$) with a Weeks-Chandler-Andersen (WCA) potential $V_{WCA}(z) = 4 ε[(σ_w/z)^{12}-(σ_w/z)^6 + 1/4]$, with range $σ_w = σ/2$. We vary the strength $ε$ over a wide range and the case of simple hard walls is also treated for comparison. By the variation of $ε$ one can change both the surface excess packing fraction and the wall-fluid $(γ_{wf})$ and wall-crystal $(γ_{wc})$ surface free energies. Several different methods to extract $γ_{wf}$ and $γ_{wc}$ from Monte Carlo (MC) simulations are implemented, and their accuracy and efficiency is comparatively discussed. The density functional theory (DFT) using Fundamental Measure functionals is found to be quantitatively accurate over a wide range of packing fractions; small deviations between DFT and MC near the fluid to crystal transition need to be studied further. Our results on density profiles near soft walls could be useful to interpret corresponding experiments with suitable colloidal dispersions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Mesophase formation in two-component cylindrical bottle-brush polymers

When two types of side chains (A,B) are densely grafted to a (stiff) backbone and the resulting bottle-brush polymer is in a solution under poor solvent conditions, an incompatibility between A and B leads to microphase separation in the resulting cylindrical brush. The possible types of ordering are reminiscent of the ordering of block copolymers in cylindrical confinement. Starting from this analogy, Leibler's theory of microphase separation in block copolymer melts is generalized to derive a description of the system in the weak segregation limit. Also molecular dynamics simulation results of a corresponding coarse-grained bead-spring model are presented. Using side chain lengths up to N = 50 effective monomers, the ratio of the Lennard-Jones energy parameter between unlike monomers $(ε_{AB})$ and monomers of the same kind $(ε_{AA} = ε_{BB})$ is varied. Various correlation functions are analyzed to study the conditions when (local) Janus cylinder-type ordering and when (local) microphase separation in the direction along the cylinder axis occurs. Both the analytical theory and the simulations give evidence for short range order due to a tendency towards microphase separation in the axial direction, with a wavelength proportional to the side chain gyration radius, irrespective of temperature and grafting density, for a wide range of these parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Methods to Compute Pressure and Wall Tension in Fluids containing Hard Particles

Colloidal systems are often modelled as fluids of hard particles (possibly with an additional soft attraction, e.g. caused by polymers also contained in the suspension). in simulations of such systems, the virial theorem cannot be straightforwardly applied to obtain the components of the pressure tensor. In systems confined by walls, it is hence also not straightforward to extract the excess energy due to the wall (the "wall tension") from the pressure tensor anisotropy. A comparative evaluation of several methods to circumvent this problem is presented, using as examples fluids of hard spheres and the Asakura-Oosawa model of colloid-polymer mixtures with a size ratio $q=0.15$ (for which the effect of the polymers can be integrated out to yield an effective attractive potential between the colloids). Factors limiting the accuracy of the various methods are carefully discussed, and controlling these factors very good mutual agreement between the various methods is found.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Monte Carlo Methods for Estimating Interfacial Free Energies and Line Tensions

Excess contributions to the free energy due to interfaces occur for many problems encountered in the statistical physics of condensed matter when coexistence between different phases is possible (e.g. wetting phenomena, nucleation, crystal growth, etc.). This article reviews two methods to estimate both interfacial free energies and line tensions by Monte Carlo simulations of simple models, (e.g. the Ising model, a symmetrical binary Lennard-Jones fluid exhibiting a miscibility gap, and a simple Lennard-Jones fluid). One method is based on thermodynamic integration. This method is useful to study flat and inclined interfaces for Ising lattices, allowing also the estimation of line tensions of three-phase contact lines, when the interfaces meet walls (where "surface fields" may act). A generalization to off-lattice systems is described as well. The second method is based on the sampling of the order parameter distribution of the system throughout the two-phase coexistence region of the model. Both the interface free energies of flat interfaces and of (spherical or cylindrical) droplets (or bubbles) can be estimated, including also systems with walls, where sphere-cap shaped wall-attached droplets occur. The curvature-dependence of the interfacial free energy is discussed, and estimates for the line tensions are compared to results from the thermodynamic integration method. Basic limitations of all these methods are critically discussed, and an outlook on other approaches is given.

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

Capillary Condensation in Cylindrical Pores: Monte Carlo Study of the Interplay of Surface and Finite Size Effects

When a fluid that undergoes a vapor to liquid transition in the bulk is confined to a long cylindrical pore, the phase transition is shifted (mostly due to surface effects at the walls of the pore) and rounded (due to finite size effects). The nature of the phase coexistence at the transition depends on the length of the pore: For very long pores the system is axially homogeneous at low temperatures. At the chemical potential where the transition takes place fluctuations occur between vapor-like and liquid-like states of the cylinder as a whole. At somewhat higher temperatures (but still far below bulk criticality) the system at phase coexistence is in an axially inhomogeneous multi-domain state, where long cylindrical liquid-like and vapor-like domains alternate. Using Monte Carlo simulations for the Ising/lattice gas model and the Asakura-Oosawa model of colloid-polymer mixtures the transition between these two different scenarios is characterized. It is shown that the density distribution changes gradually from a double-peak structure to a triple-peak shape, and the correlation length in axial direction (measuring the equilibrium domain length) becomes much smaller than the cylinder length. The (rounded) transition to the disordered phase of the fluid occurs when the axial correlation length has decreased to a value comparable to the cylinder diameter. It is also suggested that adsorption hysteresis vanishes when the transition from the simple domain state to the multi-domain state of the cylindrical pore occurs. We predict that the difference between the pore critical temperature and the hysteresis critical emperature should increase logarithmically with the length of the pore.

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

Does Young's equation hold on the nanoscale? A Monte Carlo test for the binary Lennard-Jones fluid

When a phase-separated binary ($A+B$) mixture is exposed to a wall, that preferentially attracts one of the components, interfaces between A-rich and B-rich domains in general meet the wall making a contact angle $θ$. Young's equation describes this angle in terms of a balance between the $A-B$ interfacial tension $γ_{AB}$ and the surface tensions $γ_{wA}$, $γ_{wB}$ between, respectively, the $A$- and $B$-rich phases and the wall, $γ_{AB} \cos θ=γ_{wA}-γ_{wB}$. By Monte Carlo simulations of bridges, formed by one of the components in a binary Lennard-Jones liquid, connecting the two walls of a nanoscopic slit pore, $θ$ is estimated from the inclination of the interfaces, as a function of the wall-fluid interaction strength. The information on the surface tensions $γ_{wA}$, $γ_{wB}$ are obtained independently from a new thermodynamic integration method, while $γ_{AB}$ is found from the finite-size scaling analysis of the concentration distribution function. We show that Young's equation describes the contact angles of the actual nanoscale interfaces for this model rather accurately and location of the (first order) wetting transition is estimated.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Monte Carlo Test of the Classical Theory for Heterogeneous Nucleation Barriers

Flat walls facilitate the condensation of a supersaturated vapor: Classical theory of heterogeneous nucleation predicts that the free energy barrier $ΔF_{\rm het}^*$ which needs to be overcome for the formation of sphere-cap shaped nucleation seeds is smaller than the barrier $ΔF^*_{\rm hom}$ for spherical droplets in the bulk by a factor $0<f(θ)<1$, which only depends on the contact angle $θ$. In this letter we compute both $ΔF^*_{\rm hom}$ and $ΔF^*_{\rm het}$ from Monte Carlo simulations and test the theory for the lattice gas model (for which $θ$ can be readily controlled). Even though the theory is only based on macroscopic arguments, it is shown to hold for experimentally relevant nanoscopic nucleation seeds ($20\leqΔF^*_{\rm hom}/k_BT\leq 200)$ if (independently estimated) line tension effects are considered.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Out of Equilibrium Characteristics of a Forced Translocating Chain through a Nanopore

Polymer translocation through a nano-pore in a thin membrane is studied using a coarse-grained bead-spring model and Langevin dynamics simulation with a particular emphasis to explore out of equilibrium characteristics of the translocating chain. We analyze the out of equilibrium chain conformations both at the $cis$ and the $trans$ side separately either as a function of the time during the translocation process or as as function of the monomer index $m$ inside the pore. A detailed picture of translocation emerges by monitoring the center of mass of the translocating chain, longitudinal and transverse components of the gyration radii and the end to end vector. We observe that polymer configurations at the $cis$ side are distinctly different from those at the $trans$ side. During the translocation, and immediately afterwards, the chain is clearly out of equilibrium, as different parts of the chain are characterized by a series of effective Flory exponents. We further notice that immediately after the translocation the last set of beads that have just translocated take a relatively compact structure compared to the first set of beads that translocated earlier, and the chain immediately after translocation is described by an effective Flory exponent $0.45 \pm 0.01$. The analysis of these results is further strengthened by looking at the conformations of chain segments of equal length as they cross from the $cis$ to the $trans$ side, We discuss implications of these results to the theoretical estimates and numerical simulation studies of the translocation exponent reported by various groups.

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
preprint2010arXiv

Rounding of Phase Transitions in Cylindrical Pores

Phase transitions of systems confined in long cylindrical pores (capillary condensation, wetting, crystallization, etc.) are intrinsically not sharply defined but rounded. The finite size of the cross section causes destruction of long range order along the pore axis by spontaneous nucleation of domain walls. This rounding is analyzed for two models (Ising/lattice gas and Asakura-Oosawa model for colloid-polymer mixtures) by Monte Carlo simulations and interpreted by a phenomenological theory. We show that characteristic differences between the behavior of pores of finite length and infinitely long pores occur. In pores of finite length a rounded transition occurs first, from phase coexistence between two states towards a multi-domain configuration. A second transition to the axially homogeneous phase follows near pore criticality.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Spherical Polymer Brushes Under Good Solvent Conditions: Molecular Dynamics Results Compared to Density Functional Theory

A coarse grained model for flexible polymers end-grafted to repulsive spherical nanoparticles is studied for various chain lengths and grafting densities under good solvent conditions, by Molecular Dynamics methods and density functional theory. With increasing chain length the monomer density profile exhibits a crossover to the star polymer limit. The distribution of polymer ends and the linear dimensions of individual polymer chains are obtained, while the inhomogeneous stretching of the chains is characterized by the local persistence lengths. The results on the structure factor of both single chain and full spherical brush as well as the range of applicability of the different theoretical tools are presented. Eventually an outlook on experiments is given.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Two-state protein-like folding of a homopolymer chain

Many small proteins fold via a first-order "all-or-none" transition directly from an expanded coil to a compact native state. Here we study an analogous direct freezing transition from an expanded coil to a compact crystallite for a simple flexible homopolymer. Wang-Landau sampling is used to construct the 1D density of states for square-well chains of length 128. Analysis within both the micro-canonical and canonical ensembles shows that, for a chain with sufficiently short-range interactions, the usual polymer collapse transition is preempted by a direct freezing or "folding" transition. A 2D free-energy landscape, built via subsequent multi-canonical sampling, reveals a dominant folding pathway over a single free-energy barrier. This barrier separates a high entropy ensemble of unfolded states from a low entropy set of crystallite states and the transition proceeds via the formation of a transition-state folding nucleus. Despite the non-unique homopolymer ground state, the thermodynamics of this direct freezing transition are identical to the thermodynamics of two-state protein folding. The model chain satisfies the van't Hoff calorimetric criterion for two-state folding and an Arrhenius analysis of the folding/unfolding free energy barrier yields a Chevron plot characteristic of small proteins.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Effective interactions and melting of a one dimensional defect lattice within a two-dimensional confined colloidal solid

We report Monte Carlo studies of a two-dimensional soft colloidal crystal confined in a strip geometry by parallel walls. The wall-particle interaction has corrugations along the length of the strip. Compressing the crystal by decreasing the distance between the walls induces a structural transition characterized by the sudden appearance of a one-dimensional array of extended defects each of which span several lattice parameters, a "soliton staircase". We obtain the effective interaction between these defects. A Lindemann criterion shows that the reduction of dimensionality causes a finite periodic chain of these defects to readily melt as the temperature is raised. We discuss possible experimental realizations and speculate on potential applications.

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
preprint2009arXiv

Soliton Staircases and Standing Strain Waves in Confined Colloidal Crystals

We show by computer simulation of a two-dimensional crystal confined by corrugated walls that confinement can be used to impose a controllable mesoscopic superstructure of predominantly mechanical elastic character. Due to an interplay of the particle density of the system and the width D of the confining channel, "soliton staircases" can be created along both parallel confining boundaries, that give rise to standing strain waves in the entire crystal. The periodicity of these waves is of the same order as D. This mechanism should be useful for structure formation in the self-assembly of various nanoscopic materials.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Scaling exponents of Forced Polymer Translocation through a nano-pore

We investigate several scaling properties of a translocating homopolymer through a thin pore driven by an external field present inside the pore only using Langevin Dynamics (LD) simulation in three dimension (3D). Specifically motivated by several recent theoretical and numerical studies that are apparently at odds with each other, we determine the chain length dependence of the scaling exponents of the average translocation time, the average velocity of the center of mass, $<v_{CM}>$, the effective radius of gyration during the translocation process, and the scaling exponent of the translocation coordinate ($s$-coordinate) as a function of the translocation time. We further discuss the possibility that in the case of driven translocation the finite pore size and its geometry could be responsible that the veclocity scaling exponent is less than unity and discuss the dependence of the scaling exponents on the pore geometry for the range of $N$ studied here.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Static and dynamic properties of the interface between a polymer brush and a melt of identical chains

Molecular dynamics simulations of a short-chain polymer melt between two brush-covered surfaces under shear have been performed. The end-grafted polymers which constitute the brush have the same chemical properties as the free chains in the melt and provide a soft deformable substrate. Polymer chains are described by a coarse-grained bead-spring model with Lennard-Jones interactions between the beads and a FENE potential between nearest neighbors along the backbone of the chains. The grafting density of the brush layer offers a way of controlling the behavior of the surface without altering the molecular interactions. We perform equilibrium and non-equilibrium Molecular Dynamics simulations at constant temperature and volume using the Dissipative Particle Dynamics thermostat. The equilibrium density profiles and the behavior under shear are studied as well as the interdigitation of the melt into the brush, the orientation on different length scales (bond vectors, radius of gyration, and end-to-end vector) of free and grafted chains, and velocity profiles. The viscosity and slippage at the interface are calculated as functions of grafting density and shear velocity.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Finite-size scaling above the upper critical dimension revisited: The case of the five-dimensional Ising model

Monte Carlo results for the moments <M^k> of the magnetization distribution of the nearest-neighbor Ising ferromagnet in a L^d geometry, where L (4 \leq L \leq 22) is the linear dimension of a hypercubic lattice with periodic boundary conditions in d=5 dimensions, are analyzed in the critical region and compared to a recent theory of Chen and Dohm (CD) [X.S. Chen and V. Dohm, Int. J. Mod. Phys. C (1998)]. We show that this finite-size scaling theory (formulated in terms of two scaling variables) can account for the longstanding discrepancies between Monte Carlo results and the so-called ``lowest-mode'' theory, which uses a single scaling variable tL^{d/2} where t=T/T_c-1 is the temperature distance from the critical temperature, only to a very limited extent. While the CD theory gives a somewhat improved description of corrections to the ``lowest-mode'' results (to which the CD theory can easily be reduced in the limit t \to 0, L \to \infty, tL^{d/2} fixed) for the fourth-order cumulant, discrepancies are found for the susceptibility (L^d <M^2>). Reasons for these problems are briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Molecular Dynamics Computer Simulation of the Dynamics of Supercooled Silica

We present the results of a large scale computer simulation of supercooled silica. We find that at high temperatures the diffusion constants show a non-Arrhenius temperature dependence whereas at low temperature this dependence is also compatible with an Arrhenius law. We demonstrate that at low temperatures the intermediate scattering function shows a two-step relaxation behavior and that it obeys the time temperature superposition principle. We also discuss the wave-vector dependence of the nonergodicity parameter and the time and temperature dependence of the non-Gaussian parameter.

0 citations0 reviewsNo code linkedOpen access