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

Ralf Everaers

Ralf Everaers 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 PhysicsBiomoleculesphysics.chem-ph

Catalog footprint

What is connected

15works
5topics
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

15 published item(s)

preprint2020arXiv

Single-molecule stretching experiments of flexible (wormlike) chain molecules in different ensembles: Theory and a potential application of finite chain length effects to nick-counting in DNA

We propose a formalism for deriving force-elongation and elongation-force relations for flexible chain molecules from analytical expressions for their radial distribution function, which provides insight into the factors controlling the asymptotic behavior and finite chain length corrections. In particular, we apply this formalism to our previously developed interpolation formula for the wormlike chain end-to-end distance distribution. The resulting expression for the asymptotic limit of infinite chain length is of similar quality as the numerical evaluation of Marko's and Siggia's variational theory and considerably more precise than their interpolation formula. A comparison to numerical data suggests, that our analytical expressions for the finite-chain length corrections are of similar quality. As an application of our results we discuss the possibility of inferring the changing number of nicks in a double-stranded DNA molecule in single-molecule stretching experiments from the accompanying changes in the effective chain length.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Characteristic time and length scales in melts of Kremer-Grest bead-spring polymers with wormlike bending stiffness

The Kremer-Grest (KG) model is a standard for studying generic polymer properties. Here we have equilibrated KG melts up to and beyond $200$ entanglements per chain for varying chain stiffness. We present methods for estimating the Kuhn length corrected for incompressibility effects, for estimating the entanglement length corrected for chain stiffness, for estimating bead frictions and Kuhn times taking into account entanglement effects. These are the key parameters for enabling quantitative, accurate, and parameter free comparisons between theory, experiment and simulations of KG polymer models with varying stiffness. We demonstrate this for the mean-square monomer displacements in moderately to highly entangled melts as well as for the shear relaxation modulus for unentangled melts, which are found to be in excellent agreement with the predictions from standard theories of polymer dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Kremer-Grest models for commodity polymer melts: Linking theory, experiment and simulation at the Kuhn scale

The Kremer-Grest (KG) polymer model is a standard model for studying generic polymer properties in Molecular Dynamics simulations. It owes its popularity to its simplicity and computational efficiency, rather than its ability to represent specific polymers species and conditions. Here we show, that by tuning the chain stiffness it is possible to adapt the KG model to model melts of real polymers. In particular, we provide mapping relations from KG to SI units for a wide range of commodity polymers. The connection between the experimental and the KG melts is made at the Kuhn scale, i.e. at the crossover from chemistry-specific small scale to the universal large scale behavior. We expect Kuhn scale-mapped KG models to faithfully represent universal properties dominated by the large scale conformational statistics and dynamics of flexible polymers. In particular, we observe very good agreement between entanglement moduli of our KG models and the experimental moduli of the target polymers.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Computer simulations of melts of randomly branching polymers

Randomly branching polymers with {\em annealed} connectivity are model systems for ring polymers and chromosomes. In this context, the branched structure represents transient folding induced by topological constraints. Here we present computer simulations of melts of annealed randomly branching polymers of $3 \le N \le 1800$ segments in $d=2$ and $d=3$ dimensions. In all cases, we perform a detailed analysis of the observed tree connectivities and spatial conformations. Our results are in excellent agreement with an asymptotic scaling of the average tree size of $R \sim N^{1/d}$, suggesting that the trees behave as compact, {\it territorial} fractals. The observed swelling relative to the size of ideal trees, $R\sim N^{1/4}$, demonstrates that excluded volume interactions are only partially screened in melts of annealed trees. Overall, our results are in good qualitative agreement with the predictions of Flory theory. In particular, we find that the trees swell by the combination of modified branching and path stretching. However, the former effect is subdominant and difficult to detect in $d=3$ dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Computer simulations of randomly branching polymers: Annealed vs. quenched branching structures

We present computer simulations of three systems of randomly branching polymers in d=3 dimensions: ideal trees and self-avoiding trees with annealed and quenched connectivities. In all cases, we performed a detailed analysis of trees connectivities, spatial conformations and statistical properties of linear paths on trees, and compare the results to the corresponding predictions of Flory theory. We confirm that, overall, the theory predicts correctly that trees with quenched ideal connectivity exhibit {\it less} overall swelling in good solvent than corresponding trees with annealed connectivity even though they are more strongly stretched on the path level. At the same time, we emphasize the inadequacy of the Flory theory in predicting the behaviour of other, and equally relevant, observables like contact probabilities between tree nodes. We show, then, that contact probabilities can be aptly characterized by introducing a novel critical exponent, $θ_{path}$, which accounts for how they decay as a function of the node-to-node path distance on the tree.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Multiscale approach to equilibrating model polymer melts

We present an effective and simple multiscale method for equilibrating Kremer Grest model polymer melts of varying stiffness. In our approach, we progressively equilibrate the melt structure above the tube scale, inside the tube and finally at the monomeric scale. We make use of models designed to be computationally effective at each scale. Density fluctuations in the melt structure above the tube scale are minimized through a Monte Carlo simulated annealing of a lattice polymer model. Subsequently the melt structure below the tube scale is equilibrated via the Rouse dynamics of a force-capped Kremer-Grest model that allows chains to partially interpenetrate. Finally the Kremer-Grest force field is introduced to freeze the topological state and enforce correct monomer packing. We generate $15$ melts of $500$ chains of $10.000$ beads for varying chain stiffness as well as a number of melts with $1.000$ chains of $15.000$ monomers. To validate the equilibration process we study the time evolution of bulk, collective and single-chain observables at the monomeric, mesoscopic and macroscopic length scales. Extension of the present method to longer, branched or polydisperse chains and/or larger system sizes is straight forward.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Inferring coarse-grain histone-DNA interaction potentials from high-resolution structures of the nucleosome

The histone-DNA interaction in the nucleosome is a fundamental mechanism of genomic compaction and regulation, which remains largely unkown despite a growing structural knowledge of the complex. Here, we propose a framework for the extraction of a nanoscale histone-DNA force-field from a collection of high-resolution structures, which may be adapted to a larger class of protein-DNA complexes. We apply the procedure on a large crystallographic database extended by snapshots from molecular dynamics simulations. The comparison of the structural models first shows that, at the sites of histone-DNA contact, the DNA base-pairs are locally shifted outwards, consistent with locally repulsive forces exerted by the histones. In a second step, we show that the various force profiles of the analyzed structures derive locally from a unique, sequence-independent, quadratic repulsive force field, while the sequence preferences are entirely due to the internal DNA mechanics. We thus obtain the first knowledge-derived nanoscale potential for the histone-DNA interaction in the nucleosome. The conformations obtained by relaxation of nucleosomal DNA with high-affinity sequences in this potential accurately reproduce experimental values of binding preferences. We finally address the more generic binding mechanisms relevant to the 80% genomic sequences incorporated in nucleosomes, by computing the conformation of nucleosomal DNA with sequence-averaged properties. This conformation is found to differ from those found in crystals, and the analysis suggests that repulsive histone forces are related to a local stretch tension in nucleosomal DNA, mostly between successive contact points. This tension could play a role in the stability of the complex.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Ring polymers in the melt state: the physics of crumpling

The conformational statistics of ring polymers in melts or dense solutions is strongly affected by their quenched microscopic topological state. The effect is particularly strong for non-concatenated unknotted rings, which are known to crumple and segregate and which have been implicated as models for the generic behavior of interphase chromosomes. Here we use a computationally efficient multi-scale approach to show that melts of rings of total contour length $L_r$ can be {\it quantitatively} mapped onto melts of {\em interacting} lattice trees with gyration radii $\langle R_g^2(L_r) \rangle \propto L_r^{2ν}$ and $ν=0.32\pm0.01$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Statistical Mechanics of Hamiltonian Adaptive Resolution Simulations

The Adaptive Resolution Scheme (AdResS) is a hybrid scheme that allows one to treat a molecular system with different levels of resolution depending on the location of the molecules. The construction of a Hamiltonian based on the this idea (H-AdResS) allows one to formulate the usual tools of ensembles and statistical mechanics. We present a number of exact and approximate results that provide a statistical mechanics foundation for this simulation method. We also present simulation results that illustrate the theory.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Monte Carlo adaptive resolution simulation of multicomponent molecular liquids

Complex soft matter systems can be efficiently studied with the help of adaptive resolution simulation methods, concurrently employing two levels of resolution in different regions of the simulation domain. The non-matching properties of high- and low-resolution models, however, lead to thermodynamic imbalances between the system's subdomains. Such inhomogeneities can be healed by appropriate compensation forces, whose calculation requires nontrivial iterative procedures. In this work we employ the recently developed Hamiltonian Adaptive Resolution Simulation method to perform Monte Carlo simulations of a binary mixture, and propose an efficient scheme, based on Kirkwood Thermodynamic Integration, to regulate the thermodynamic balance of multi-component systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Hamiltonian adaptive resolution simulation for molecular liquids

Adaptive resolution schemes allow the simulation of a molecular fluid treating simultaneously different subregions of the system at different levels of resolution. In this work we present a new scheme formulated in terms of a global Hamiltonian. Within this approach equilibrium states corresponding to well defined statistical ensembles can be generated making use of all standard Molecular Dynamics or Monte Carlo methods. Models at different resolutions can thus be coupled, and thermodynamic equilibrium can be modulated keeping each region at desired pressure or density without disrupting the Hamiltonian framework.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Bubble statistics and positioning in superhelically stressed DNA

We present a general framework to study the thermodynamic denaturation of double-stranded DNA under superhelical stress. We report calculations of position- and size-dependent opening probabilities for bubbles along the sequence. Our results are obtained from transfer-matrix solutions of the Zimm-Bragg model for unconstrained DNA and of a self-consistent linearization of the Benham model for superhelical DNA. The numerical efficiency of our method allows for the analysis of entire genomes and of random sequences of corresponding length ($10^6-10^9$ base pairs). We show that, at physiological conditions, opening in superhelical DNA is strongly cooperative with average bubble sizes of $10^2-10^3$ base pairs (bp), and orders of magnitude higher than in unconstrained DNA. In heterogeneous sequences, the average degree of base-pair opening is self-averaging, while bubble localization and statistics are dominated by sequence disorder. Compared to random sequences with identical GC-content, genomic DNA has a significantly increased probability to open large bubbles under superhelical stress. These bubbles are frequently located directly upstream of transcription start sites.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Topological versus rheological entanglement length in primitive path analysis protocols

Primitive path analysis algorithms are now routinely employed to analyze entanglements in computer simulations of polymeric systems, but different analysis protocols result in different estimates of the entanglement length, N_e. Here we argue that standard PPA measures the rheological entanglement length, typically employed by tube models and relevant to quantitative comparisons with experiment, while codes like Z or CReTA also determine the topological entanglement length. For loosely entangled systems, a simple analogy between between phantom networks and the mesh of entangled primitive paths suggests a factor of two between the two numbers. This result is in excellent agreement with reported values for poly-ethylene, poly-butadiene and bead-spring polymer melts.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Stress Relaxation in Entangled Polymer Melts

We present an extensive set of simulation results for the stress relaxation in equilibrium and step-strained bead-spring polymer melts. The data allow us to explore the chain dynamics and the shear relaxation modulus, $G(t)$, into the plateau regime for chains with $Z=40$ entanglements and into the terminal relaxation regime for $Z=10$. Using the known (Rouse) mobility of unentangled chains and the melt entanglement length determined via the primitive path analysis of the microscopic topological state of our systems, we have performed parameter -free tests of several different tube models. We find excellent agreement for the Likhtman-McLeish theory using the double reptation approximation for constraint release, if we remove the contribution of high-frequency modes to contour length fluctuations of the primitive chain.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

DNA: From rigid base-pairs to semiflexible polymers

The sequence-dependent elasticity of double-helical DNA on a nm length scale can be captured by the rigid base-pair model, whose strains are the relative position and orientation of adjacent base-pairs. Corresponding elastic potentials have been obtained from all-atom MD simulation and from high-resolution structural data. On the scale of a hundred nm, DNA is successfully described by a continuous worm-like chain model with homogeneous elastic properties characterized by a set of four elastic constants, which have been directly measured in single-molecule experiments. We present here a theory that links these experiments on different scales, by systematically coarse-graining the rigid base-pair model for random sequence DNA to an effective worm-like chain description. The average helical geometry of the molecule is exactly taken into account in our approach. We find that the available microscopic parameters sets predict qualitatively similar mesoscopic parameters. The thermal bending and twisting persistence lengths computed from MD data are 42 and 48 nm, respectively. The static persistence lengths are generally much higher, in agreement with cyclization experiments. All microscopic parameter sets predict negative twist-stretch coupling. The variability and anisotropy of bending stiffness in short random chains lead to non-Gaussian bend angle distributions, but become unimportant after two helical turns.

0 citations0 reviewsNo code linkedOpen access