No Arabic abstract
A lattice model is presented for the simulation of dynamics in polymeric systems. Each polymer is represented as a chain of monomers, residing on a sequence of nearest-neighbor sites of a face-centered-cubic lattice. The polymers are self- and mutually avoiding walks: no lattice site is visited by more than one polymer, nor revisited by the same polymer after leaving it. The dynamics occurs through single-monomer displacements over one lattice spacing. To demonstrate the high computational efficiency of the model, we simulate a dense binary polymer mixture with repelling nearest-neighbor interactions between the two types of polymers, and observe the phase separation over a long period of time. The simulations consist of a total of 46,080 polymers, 100 monomers each, on a lattice with 13,824,000 sites, and an interaction strength of 0.1 kT. In the final two decades of time, the domain-growth is found to be d(t) ~ t^1/3, as expected for a so-called Model B system.
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 $rho(z)$ in the $z$-direction perpendicular to the surface for rings and linear chains are practically identical, $rho_R(2 N_L, z)=rho_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.
Surface segregation of the low-molecular weight component in a polymeric mixture leads to degradation of industrial formulations. We report a simultaneous phase separation and surface migration phenomena in oligomer-polymer and oligomer-gel systems following a temperature quench. We compute equilibrium and time varying migrant density profiles and wetting layer thickness using coarse grained molecular dynamics and mesoscale hydrodynamics simulations to demonstrate that surface migration in oligomer-gel systems is significantly reduced due to network elasticity. Further, phase separation processes are significantly slowed in gels, modifying the Lifshitz-Slyozov-Wagner (LSW) law $ell(tau) sim tau^{1/3}$. Our work allows for rational design of polymer/gel-oligomer mixtures with predictable surface segregation characteristics.
One of the most promising applications in nanoscience is the design of new materials to improve water permeability and selectivity of nanoporous membranes. Understanding the molecular architecture behind these fascinating structures and how it impacts the water flow is an intricate but necessary task. We studied here, the water flux through multi-layered nanoporous molybdenum disulfide (MLNMoS$_2$) membranes with different nanopore sizes and length. Molecular dynamics simulations show that the permeability do not increase with the inverse of the membrane thickness, violating the classical hydrodynamic behavior. The data also reveals that the water dynamics is slower than that observed in frictionless carbon nanotubes and multi-layer graphene membranes, which we explain in terms of an anchor mechanism observed in between layers. We show that the membrane permeability is critically dependent on the nanopore architecture, bringing important insights into the manufacture of new desalination membranes.
We study the phenomenon of migration of the small molecular weight component of a binary polymer mixture to the free surface using mean field and self-consistent field theories. By proposing a free energy functional that incorporates polymer-matrix elasticity explicitly, we compute the migrant volume fraction and show that it decreases significantly as the sample rigidity is increased. Estimated values of the bulk modulus suggest that the effect should be observable experimentally for rubber-like materials. This provides a simple way of controlling surface migration in polymer mixtures and can play an important role in industrial formulations, where surface migration often leads to decreased product functionality.
The phase diagram of star polymer solutions in a good solvent is obtained over a wide range of densities and arm numbers by Monte Carlo simulations. The effective interaction between the stars is modeled by an ultrasoft pair potential which is logarithmic in the core-core distance. Among the stable phases are a fluid as well as body-centered cubic, face-centered cubic, body-centered orthogonal, and diamond crystals. In a limited range of arm numbers, reentrant melting and reentrant freezing transitions occur for increasing density.