Related papers: Fast convergence to equilibrium for long-chain pol…
The widely used double-bridging hybrid (DBH) method for equilibrating simulated entangled polymer melts [R. Auhl et al., J. Chem. Phys. v. 119, p. 12718, 2003] loses its effectiveness as chain stiffness increases into the semiflexible…
Several methods for preparing well equilibrated melts of long chains polymers are studied. We show that the standard method in which one starts with an ensemble of chains with the correct end-to-end distance arranged randomly in the…
Equilibration of polymer melts containing highly entangled long polymer chains in confinement or with free surfaces is a challenge for computer simulations. We approach this problem by first studying polymer melts based on the soft-sphere…
Dimerization and subsequent aggregation of polymers and biopolymers often occur under nonequilibrium conditions. When the initial state of the polymer is not collapsed or the final folded native state, the dynamics of dimerization can…
An off-lattice Monte Carlo algorithm for solutions of equilibrium polymers (EP) is proposed. At low and moderate densities this is shown to reproduce faithfully the (static) properties found recently for flexible linear EP using a lattice…
In this paper, a hybrid quasi-static atomistic simulation method at finite temperature is developed, which combines the advantages of MD for thermal equilibrium and atomic-scale finite element method (AFEM) for efficient equilibration. Some…
We discuss in detail a recently proposed hybrid particle-continuum scheme for complex fluids and evaluate it at the example of a confined homopolymer solution in slit geometry. The hybrid scheme treats polymer chains near the impenetrable…
Monte Carlo simulations of proteins are hindered by the system's high density which often makes local moves ineffective. Here we devise and test a set of long range moves that work well even when all sites of a lattice simulation are…
We report a novel hybrid method of simultaneous atomistic simulation of solids in critical regions (contacts surfaces, cracks areas, etc.), along with continuum modeling of other parts. The continuum is treated in terms of quasi-atoms of…
In this paper, we present a family of new mixed finite element methods for linear elasticity for both spatial dimensions $n=2,3$, which yields a conforming and strongly symmetric approximation for stress. Applying…
We generalize a recently investigated lattice model of semiflexible polymers formed under equilibrium polymerization in a solution and conduct a comprehensive investigation of its melting properties. The model is characterized by six…
We propose an efficient Monte Carlo algorithm for the off-lattice simulation of dense hard sphere polymer melts using cluster moves, called event chains, which allow for a rejection-free treatment of the excluded volume. Event chains also…
We present a computationally efficient multiscale method for preparing equilibrated, isotropic long chain model polymer melts. As an application we generate Kremer-Grest melts of $1000$ chains with $200$ entanglements and $25000$-$2000$…
The paper proposes a novel hybrid method for solving equilibrium problems and fixed point problems. By constructing specially cutting-halfspaces, in this algorithm, only an optimization program is solved at each iteration without the…
Mean-field techniques provide a rather accurate description of single-chain conformations in spatially inhomogeneous polymer systems containing interfaces or surfaces. Intermolecular correlations, however, are not described by the…
A previously-developed hybrid particle-continuum method [J. B. Bell, A. Garcia and S. A. Williams, SIAM Multiscale Modeling and Simulation, 6:1256-1280, 2008] is generalized to dense fluids and two and three dimensional flows. The scheme…
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…
It was recently demonstrated that a simple Monte Carlo (MC) algorithm involving the swap of particle pairs dramatically accelerates the equilibrium sampling of simulated supercooled liquids. We propose two numerical schemes integrating the…
Momentum methods for convex optimization often rely on precise choices of algorithmic parameters, based on knowledge of problem parameters, in order to achieve fast convergence, as well as to prevent oscillations that could severely…
The Langevin Equation for cooperative dynamics represents the dynamics of polymer melts with chains of increasing degree of polymerization, covering the full range of behavior from the unentangled to the entangled regime. This equation…