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 Kremer

Kurt Kremer 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-phcond-mat.mtrl-scicond-mat.str-elphysics.comp-phSubcellular Processes

Catalog footprint

What is connected

30works
9topics
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

30 published item(s)

preprint2026arXiv

Roadmap for Condensates in Cell Biology

Biomolecular condensates govern essential cellular processes yet elude description by traditional equilibrium models. This roadmap, distilled from structured discussions at a workshop and reflecting the consensus of its participants, clarifies key concepts for researchers, funding bodies, and journals. After unifying terminology that often separates disciplines, we outline the core physics of condensate formation, review their biological roles, and identify outstanding challenges in nonequilibrium theory, multiscale simulation, and quantitative in-cell measurements. We close with a forward-looking outlook to guide coordinated efforts toward predictive, experimentally anchored understanding and control of biomolecular condensates.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

FAIR data enabling new horizons for materials research

The prosperity and lifestyle of our society are very much governed by achievements in condensed matter physics, chemistry and materials science, because new products for sectors such as energy, the environment, health, mobility and information technology (IT) rely largely on improved or even new materials. Examples include solid-state lighting, touchscreens, batteries, implants, drug delivery and many more. The enormous amount of research data produced every day in these fields represents a gold mine of the twenty-first century. This gold mine is, however, of little value if these data are not comprehensively characterized and made available. How can we refine this feedstock; that is, turn data into knowledge and value? For this, a FAIR (findable, accessible, interoperable and reusable) data infrastructure is a must. Only then can data be readily shared and explored using data analytics and artificial intelligence (AI) methods. Making data 'findable and AI ready' (a forward-looking interpretation of the acronym) will change the way in which science is carried out today. In this Perspective, we discuss how we can prepare to make this happen for the field of materials science.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Proline isomerization regulates the phase behavior of elastin-like polypeptides in water

Responsiveness of polypeptides and polymers in aqueous solution plays an important role in biomedical applications and in designing advanced functional materials. Elastin-like polypeptides (ELPs) are a well-known class of synthetic intrinsically disordered proteins (IDPs), which exhibit a lower critical solution temperature (LCST) in pure water. The LCST transition can be further tuned by proline isomerization. Here, we study and compare the influence of cis/trans proline isomerization on the collapse of single ELPs in aqueous solution. Our results reveal that cis isomers play an important role in tuning the phase behavior of ELPs by hindering peptide-water hydrogen bonding while promoting intramolecular interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Free-energy landscape of polymer-crystal polymorphism

Polymorphism rationalizes how processing can control the final structure of a material. The rugged free-energy landscape and exceedingly slow kinetics in the solid state have so far hampered computational investigations. We report for the first time the free-energy landscape of a polymorphic crystalline polymer, syndiotactic polystyrene. Coarse-grained metadynamics simulations allow us to efficiently sample the landscape at large. The free-energy difference between the two main polymorphs, $α$ and $β$, is further investigated by quantum-chemical calculations. The two methods are in line with experimental observations: they predict $β$ as the more stable polymorph at standard conditions. Critically, the free-energy landscape suggests how the $α$ polymorph may lead to experimentally observed kinetic traps. The combination of multiscale modeling, enhanced sampling, and quantum-chemical calculations offers an appealing strategy to uncover complex free-energy landscapes with polymorphic behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Ligand-protein interactions in lysozyme investigated through a dual-resolution model

A fully atomistic modelling of biological macromolecules at relevant length- and time-scales is often cumbersome or not even desirable, both in terms of computational effort required and it a posteriori analysis. This difficulty can be overcome with the use of multi-resolution models, in which different regions of the same system are concurrently described at different levels of detail. In enzymes, computationally expensive atomistic detail is crucial in the modelling of the active site in order to capture e.g. the chemically subtle process of ligand binding. In contrast, important yet more collective properties of the remainder of the protein can be reproduced with a coarser description. In the present work, we demonstrate the effectiveness of this approach through the calculation of the binding free energy of hen egg white lysozyme (HEWL) with the inhibitor di-N-acetylchitotriose. Particular attention is posed to the impact of the mapping, i.e. the selection of atomistic and coarse-grained residues, on the binding free energy. It is shown that, in spite of small variations of the binding free energy with respect to the active site resolution, the separate contributions coming from different energetic terms (such as electrostatic and van der Waals interactions) manifest a stronger dependence on the mapping, thus pointing to the existence of an optimal level of intermediate resolution.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Why Do Elastin-Like Polypeptides Possibly Have Different Solvation Behaviors in Water-Ethanol and Water-Urea Mixtures?

The solvent quality determines the collapsed or the expanded state of a polymer. For example, a polymer dissolved in a poor solvent collapses, whereas in a good solvent it opens up. While this standard understanding is generally valid, there are examples when a polymer collapses even in a mixture of two good solvents. This phenomenon, commonly known as co-non-solvency, is usually associated with smart polymers. Moreover, recent experiments have shown that the elastin-like polypeptides (ELPs) show co-non-solvency behavior in aqueous-ethanol mixtures. In this study, we investigate the phase behavior of ELPs in aqueous binary mixtures using molecular dynamics simulations of all-atom and complementary explicit solvent generic models. The model is parameterized by mapping the solvation free energy obtained from the all-atom simulations onto the generic interaction parameters. For this purpose, we derive segment based generic parameters for four different peptides, namely proline (P), valine (V), glycine (G) and alanine (A). Here we compare the conformational behavior of two ELP sequences, namely VPGGG and VPGVG, in aqueous-ethanol and -urea mixtures. Consistent with recent experiments, we find that ELPs show co-non-solvency in aqueous-ethanol mixtures. Ethanol molecules have preferential binding with all ELP residues and thus driving the coil-to-globule transition. On the contrary, ELP conformations show weak variation in aqueous-urea mixtures. Our simulations suggest that the glycine residues dictate the overall behavior of ELPs in aqueous-urea, where urea molecules have a rather weak preferential binding with glycine, i.e., less than kT. While the validation of the latter findings will require more detailed experimental investigation, the results presented here may provide a new twist to the present understanding of cosolvent interactions with peptides and proteins.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Open-Boundary Hamiltonian adaptive resolution. From grand canonical to non-equilibrium molecular dynamics simulations

We propose an open-boundary molecular dynamics method in which an atomistic system is in contact with an infinite particle reservoir at constant temperature, volume and chemical potential. In practice, following the Hamiltonian adaptive resolution strategy, the system is partitioned into a domain of interest and a reservoir of non-interacting, ideal gas, particles. An external potential, applied only in the interfacial region, balances the excess chemical potential of the system. To ensure that the size of the reservoir is infinite, we introduce a particle insertion/deletion algorithm to control the density in the ideal gas region. We show that it is possible to study non-equilibrium phenomena with this open-boundary molecular dynamics method. To this aim, we consider a prototypical confined liquid under the influence of an external constant density gradient. The resulting pressure-driven flow across the atomistic system exhibits a velocity profile consistent with the corresponding solution of the Navier-Stokes equation. In contrast to available computational methods in which external forces drive the system far from equilibrium, this approach conserves momentum and closely resembles experimental conditions. The presented method can be used to study various direct and indirect out-of-equilibrium conditions in complex molecular systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A multi-resolution model to capture both global fluctuations of an enzyme and molecular recognition in the ligand-binding site

In multi-resolution simulations, different system components are simultaneously modelled at different levels of resolution, these being smoothly coupled together. In the case of enzyme systems, computationally expensive atomistic detail is needed in the active site to capture the chemistry of substrate binding. Global properties of the rest of the protein also play an essential role, determining the structure and fluctuations of the binding site; however, these can be modelled on a coarser level. Similarly, in the most computationally efficient scheme only the solvent hydrating the active site requires atomistic detail. We present a methodology to couple atomistic and coarse-grained protein models, while solvating the atomistic part of the protein in atomistic water. This allows a free choice of which protein and solvent degrees of freedom to include atomistically, without loss of accuracy in the atomistic description. This multi-resolution methodology can successfully model stable ligand binding, and we further confirm its validity via an exploration of system properties relevant to enzymatic function. In addition to a computational speedup, such an approach can allow the identification of the essential degrees of freedom playing a role in a given process, potentially yielding new insights into biomolecular function.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Consistent interpretation of molecular simulation kinetics using Markov state models biased with external information

Molecular simulations can provide microscopic insight into the physical and chemical driving forces of complex molecular processes. Despite continued advancement of simulation methodology, model errors may lead to inconsistencies between simulated and reference (e.g., from experiments or higher-level simulations) observables. To bound the microscopic information generated by computer simulations within reference measurements, we propose a method that reweights the microscopic transitions of the system to improve consistency with a set of coarse kinetic observables. The method employs the well-developed Markov state modeling framework to efficiently link microscopic dynamics with long-time scale constraints, thereby consistently addressing a wide range of time scales. To emphasize the robustness of the method, we consider two distinct coarse-grained models with significant kinetic inconsistencies. When applied to the simulated conformational dynamics of small peptides, the reweighting procedure systematically improves the time scale separation of the slowest processes. Additionally, constraining the forward and backward rates between metastable states leads to slight improvement of their relative stabilities and, thus, refined equilibrium properties of the resulting model. Finally, we find that difficulties in simultaneously describing both the simulated data and the provided constraints can help identify specific limitations of the underlying simulation approach.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Equilibration of High Molecular-Weight Polymer Melts: A Hierarchical Strategy

A strategy is developed for generating equilibrated high molecular-weight polymer melts described with microscopic detail by sequentially backmapping coarse-grained (CG) configurations. The microscopic test model is generic but retains features like hard excluded volume interactions and realistic melt densities. The microscopic representation is mapped onto a model of soft spheres with fluctuating size, where each sphere represents a microscopic subchain with $N_{\rm b}$ monomers. By varying $N_{\rm b}$ a hierarchy of CG representations at different resolutions is obtained. Within this hierarchy, CG configurations equilibrated with Monte Carlo at low resolution are sequentially fine-grained into CG melts described with higher resolution. A Molecular Dynamics scheme is employed to slowly introduce the microscopic details into the latter. All backmapping steps involve only local polymer relaxation thus the computational efficiency of the scheme is independent of molecular weight, being just proportional to system size. To demonstrate the robustness of the approach, microscopic configurations containing up to $n=1000$ chains with polymerization degrees $N=2000$ are generated and equilibration is confirmed by monitoring key structural and conformational properties. The extension to much longer chains or branched polymers is straightforward.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

One size fits all: equilibrating chemically different polymer liquids through universal long-wavelength description

Mesoscale behavior of polymers is frequently described by universal laws. This physical property motivates us to propose a new modeling concept, grouping polymers into classes with a common long-wavelength representation. In the same class samples of different materials can be generated from this representation, encoded in a single library system. We focus on homopolymer melts, grouped according to the invariant degree of polymerization. They are described with a bead-spring model, varying chain stiffness and density to mimic chemical diversity. In a renormalization group-like fashion library samples provide a universal blob-based description, hierarchically backmapped to create configurations of other class-members. Thus large systems with experimentally-relevant invariant degree of polymerizations (so far accessible only on very coarse-grained level) can be microscopically described. Equilibration is verified comparing conformations and melt structure with smaller scale conventional simulations.

0 citations0 reviewsNo code linkedOpen access
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
preprint2016arXiv

Why does high pressure destroy co-non-solvency of PNIPAm in aqueous methanol?

It is well known that poly(N-isopropylacrylamide) (PNIPAm) exhibits an interesting, yet puzzling, phenomenon of co-non-solvency. Co-non-solvency occurs when two competing good solvents for PNIPAm, such as water and alcohol, are mixed together. As a result, the same PNIPAm collapses within intermediate mixing ratios. This complex conformational transition is driven by preferential binding of methanol with PNIPAm. Interestingly, co-non-solvency can be destroyed when applying high hydrostatic pressures. In this work, using a large scale molecular dynamics simulation employing high pressures, we propose a microscopic picture behind the suppression of the co-non-solvency phenomenon. Based on thermodynamic and structural analysis, our results suggest that the preferential binding of methanol with PNIPAm gets partially lost at high pressures, making the background fluid reasonably homogeneous for the polymer. This is consistent with the hypothesis that the co-non-solvency phenomenon is driven by preferential binding and is not based on depletion effects.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Advantages and challenges in coupling an ideal gas to atomistic models in adaptive resolution simulations

In adaptive resolution simulations, molecular fluids are modeled employing different levels of resolution in different subregions of the system. When traveling from one region to the other, particles change their resolution on the fly. One of the main advantages of such approaches is the computational efficiency gained in the coarse-grained region. In this respect the best coarse-grained system to employ in the low resolution region would be the ideal gas, making intermolecular force calculations in the coarse-grained subdomain redundant. In this case, however, a smooth coupling is challenging due to the high energetic imbalance between typical liquids and a system of non-interacting particles. In the present work, we investigate this approach, using as a test case the most biologically relevant fluid, water. We demonstrate that a successful coupling of water to the ideal gas can be achieved with current adaptive resolution methods, and discuss the issues that remain to be addressed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Co-non-solvency: Mean-field polymer theory does not describe polymer collapse transition in a mixture of two competing good solvents

Smart polymers are a modern class of polymeric materials that often exhibit unpredictable behavior in mixtures of solvents. One such phenomenon is co-non-solvency. Co-non-solvency occurs when two (perfectly) miscible and competing good solvents, for a given polymer, are mixed together. As a result, the same polymer collapses into a compact globule within intermediate mixing ratios. More interestingly, polymer collapses when the solvent quality remains good and even gets increasingly better by the addition of the better cosolvent. This is a puzzling phenomenon that is driven by strong local concentration fluctuations. Because of the discrete particle based nature of the interactions, Flory-Huggins type mean field arguments become unsuitable. In this work, we extend the analysis of the co-non-solvency effect presented earlier [Nature Communications 5, 4882 (2014)]. We explain why co-non-solvency is a generic phenomenon that can be understood by the thermodynamic treatment of the competitive displacement of (co)solvent components. This competition can result in a polymer collapse upon improvement of the solvent quality. Specific chemical details are not required to understand these complex conformational transitions. Therefore, a broad range of polymers are expected to exhibit similar reentrant coil-globule-coil transitions in competing good solvents.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

From classical to quantum and back: Hamiltonian coupling of classical and Path Integral models of atoms

In computer simulations, quantum delocalization of atomic nuclei can be modeled making use of the Path Integral (PI) formulation of quantum statistical mechanics. This approach, however, comes with a large computational cost. By restricting the PI modeling to a small region of space, this cost can be significantly reduced. In the present work we derive a Hamiltonian formulation for a bottom-up, theoretically solid simulation protocol that allows molecules to change their resolution from quantum-mechanical to classical and vice versa on the fly, while freely diffusing across the system. This approach renders possible simulations of quantum systems at constant chemical potential. The validity of the proposed scheme is demonstrated by means of simulations of low temperature parahydrogen. Potential future applications include simulations of biomolecules, membranes, and interfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Computational studies of biomembrane systems: Theoretical considerations, simulation models, and applications

This chapter summarizes several approaches combining theory, simulation and experiment that aim for a better understanding of phenomena in lipid bilayers and membrane protein systems, covering topics such as lipid rafts, membrane mediated interactions, attraction between transmembrane proteins, and aggregation in biomembranes leading to large superstructures such as the light harvesting complex of green plants. After a general overview of theoretical considerations and continuum theory of lipid membranes we introduce different options for simulations of biomembrane systems, addressing questions such as: What can be learned from generic models? When is it expedient to go beyond them? And what are the merits and challenges for systematic coarse graining and quasi-atomistic coarse grained models that ensure a certain chemical specificity?

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

From a melt of rings to chromosome territories: The role of topological constraints in genome folding

We review pro and contra of the hypothesis that generic polymer properties of topological constraints are behind many aspects of chromatin folding in eukaryotic cells. For that purpose, we review, first, recent theoretical and computational findings in polymer physics related to concentrated, topologically-simple (unknotted and unlinked) chains or a system of chains. Second, we review recent experimental discoveries related to genome folding. Understanding in these fields is far from complete, but we show how looking at them in parallel sheds new light on both.

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

Adaptive resolution molecular dynamics simulation through coupling to an internal particle reservoir

For simulation studies of (macro) molecular liquids it would be of significant interest to be able to adjust or increase the level of resolution within one region of space, while allowing for the free exchange of molecules between open regions of different resolution or representation. We generalize the adaptive resolution idea and suggest an interpretation in terms of an effective generalized grand canonical approach. The method is applied to liquid water at ambient conditions.

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

Molecular Dynamics Simulation Study of Nonconcatenated Ring Polymers in a Melt: I. Statics

Molecular dynamics simulations were conducted to investigate the structural properties of melts of nonconcatenated ring polymers and compared to melts of linear polymers. The longest rings were composed of N=1600 monomers per chain which corresponds to roughly 57 entanglement lengths for comparable linear polymers. For the rings, the radius of gyration squared was found to scale as N to the 4/5 power for an intermediate regime and N to the 2/3 power for the larger rings indicating an overall conformation of a crumpled globule. However, almost all beads of the rings are "surface beads" interacting with beads of other rings, a result also in agreement with a primitive path analysis performed in the following paper (DOI: 10.1063/1.3587138). Details of the internal conformational properties of the ring and linear polymers as well as their packing are analyzed and compared to current theoretical models.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Molecular Dynamics Simulation Study of Nonconcatenated Ring Polymers in a Melt: II. Dynamics

Molecular dynamics simulations were conducted to investigate the dynamic properties of melts of nonconcatenated ring polymers and compared to melts of linear polymers. The longest rings were composed of N=1600 monomers per chain which corresponds to roughly 57 entanglement lengths for comparable linear polymers. The ring melts were found to diffuse faster than their linear counterparts, with the self-diffusion coefficient for both architectures scaling as approximately N to the -2.4 power for large N. The mean-square displacement of the center-of-mass of the rings follows a sub-diffusive behavior for times and distances beyond the mean-square gyration radius, neither compatible with the Rouse nor the reptation model. The rings relax stress much faster than linear polymers and the zero-shear viscosity was found to vary as approximately N to the 1.4 power which is much weaker than the N to the 3.4 power of linear chains, not matching any commonly known model for polymer dynamics when compared to the observed mean-square displacements. These findings are discussed in view of the conformational properties of the rings presented in the preceding paper (DOI: 10.1063/1.3587137).

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Rheology of Ring Polymer Melts: From Linear Contaminants to Ring/Linear Blends

Ring polymers remain a major challenge to our current understanding of polymer dynamics. Experimental results are difficult to interpret because of the uncertainty in the purity and dispersity of the sample. Using both equilibrium and non-equilibrium molecular dynamics simulations we have systematically investigated the structure, dynamics and rheology of perfectly controlled ring/linear polymer blends with chains of such length and flexibility that the number of entanglements is up to about 14 per chain, which is comparable to experimental systems examined in the literature. The smallest concentration at which linear contaminants increase the zero-shear viscosity of a ring polymer melt of these chain lengths by 10% is approximately one-fifth of their overlap concentration. When the two architectures are present in equal amounts the viscosity of the blend is approximately twice as large as that of the pure linear melt. At this concentration the diffusion coefficient of the rings is found to decrease dramatically, while the static and dynamic properties of the linear polymers are mostly unaffected. Our results are supported by a primitive path analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Coupling different levels of resolution in molecular simulations

Simulation schemes that allow to change molecular representation in a subvolume of the simulation box while preserving the equilibrium with the surrounding introduce conceptual problems of thermodynamic consistency. In this work we present a general scheme based on thermodynamic arguments which ensures thermodynamic equilibrium among the molecules of different representation. The robustness of the algorithm is tested for two examples, namely an adaptive resolution simulation, atomistic/coarse-grained, for a liquid of tetrahedral molecules and an adaptive resolution simulation of a binary mixture of tetrahedral molecules and spherical solutes.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Modeling Diffusive Dynamics in Adaptive Resolution Simulation of Liquid Water

We present a dual-resolution molecular dynamics (MD) simulation of liquid water employing a recently introduced Adaptive Resolution Scheme (AdResS). The spatially adaptive molecular resolution procedure allows for changing from a coarse-grained to an all-atom representation and vice-versa on-the-fly. In order to find the most appropriate coarse-grained water model to be employed with AdResS we first study the accuracy of different coarse-grained water models in reproducing the structural properties of the all-atom system. Typically, coarse-grained molecular models have a higher diffusion constant than the corresponding all-atom models due to the reduction in degrees of freedom (DOFs) upon coarse-graining that eliminates the fluctuating forces associated with those integrated-out molecular DOFs. Here, we introduce the methodology to obtain the same diffusional dynamics across different resolutions. We show that this approach leads to the correct description of essential thermodynamic, structural and dynamical properties of liquid water at ambient conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Corrections to Scaling in the Hydrodynamic Properties of Dilute Polymer Solutions

We discuss the hydrodynamic radius $R_H$ of polymer chains in good solvent, and show that the leading order correction to the asymptotic law $R_H \propto N^ν$ ($N$ degree of polymerization, $ν\approx 0.59$) is an ``analytic'' term of order $N^{-(1 - ν)}$, which is directly related to the discretization of the chain into a finite number of beads. This result is further corroborated by exact calculations for Gaussian chains, and extensive numerical simulations of different models of good--solvent chains, where we find a value of $1.591 \pm 0.007$ for the asymptotic universal ratio $R_G / R_H$, $R_G$ being the chain's gyration radius. For $Θ$ chains the data apparently extrapolate to $R_G / R_H \approx 1.44$, which is different from the Gaussian value 1.5045, but in accordance with previous simulations. We also show that the experimentally observed deviations of the initial decay rate in dynamic light scattering from the asymptotic Benmouna--Akcasu value can partly be understood by similar arguments.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Structure of Polyelectrolytes in Poor Solvent

We present simulations on charged polymers in poor solvent. First we investigate in detail the dilute concentration range with and without imposed extension constraints. The resulting necklace polymer conformations are analyzed in detail. We find strong fluctuations in the number of pearls and their sizes leading only to small signatures in the form factor and the force-extension relation. The scaling of the peak in the structure factor with the monomer density shows a pertinent different behavior from good solvent chains.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Strong electrostatic interactions in spherical colloidal systems

We investigate spherical macroions in the strong Coulomb coupling regime within the primitive model in salt-free environment. We first show that the ground state of an isolated colloid is naturally overcharged by simple electrostatic arguments illustrated by the Gillespie rule. We furthermore demonstrate that in the strong Coulomb coupling this mechanism leads to ionized states and thus to long range attractions between like-charged spheres. We use molecular dynamics simulations to study in detail the counterion distribution for one and two highly charged colloids for the ground state as well as for finite temperatures. We compare our results in terms of a simple version of a Wigner crystal theory and find excellent qualitative and quantitative agreement.

0 citations0 reviewsNo code linkedOpen access