اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدEfficient simulation of non-crossing fibers and chains in a hydrodynamic solvent
50
0
0.0
(
0
)
تأليف
J.T. Padding
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
An efficient simulation method is presented for Brownian fiber suspensions, which includes both uncrossability of the fibers and hydrodynamic interactions between the fibers mediated by a mesoscopic solvent. To conserve hydrodynamics, collisions between the fibers are treated such that momentum and energy are conserved locally. The choice of simulation parameters is rationalised on the basis of dimensionless numbers expressing the relative strength of different physical processes. The method is applied to suspensions of semiflexible fibers with a contour length equal to the persistence length, and a mesh size to contour length ratio ranging from 0.055 to 0.32. For such fibers the effects of hydrodynamic interactions are observable, but relatively small. The non-crossing constraint, on the other hand, is very important and leads to hindered displacements of the fibers, with an effective tube diameter in agreement with recent theoretical predictions. The simulation technique opens the way to study the effect of viscous effects and hydrodynamic interactions in microrheology experiments where the response of an actively driven probe bead in a fiber suspension is measured.
قيم البحث
اقرأ أيضاً
The mobility of polymer chains in perfect polyethylene (PE) crystal was calculated as a function of temperature and chain length through Molecular dynamics (MD) in united atom approximation. The results demonstrate that the chain mobility drastically
increases in the vicinity of the phase transition from the orthorhombic to quasi-hexagonal phase. In the quasi-hexagonal phase, the chain mobility is almost independent on temperature and inversely proportional to the chain length.
Using a recently developed bead-spring model for semiflexible polymers that takes into account their natural extensibility, we report an efficient algorithm to simulate the dynamics for polymers like double-stranded DNA (dsDNA) in the absence of hydr
odynamic interactions. The dsDNA is modelled with one bead-spring element per basepair, and the polymer dynamics is described by the Langevin equation. The key to efficiency is that we describe the equations of motion for the polymer in terms of the amplitudes of the polymers fluctuation modes, as opposed to the use of the physical positions of the beads. We show that, within an accuracy tolerance level of $5%$ of several key observables, the model allows for single Langevin time steps of $approx1.6$, 8, 16 and 16 ps for a dsDNA model-chain consisting of 64, 128, 256 and 512 basepairs (i.e., chains of 0.55, 1.11, 2.24 and 4.48 persistence lengths) respectively. Correspondingly, in one hour, a standard desktop computer can simulate 0.23, 0.56, 0.56 and 0.26 ms of these dsDNA chains respectively. We compare our results to those obtained from other methods, in particular, the (inextensible discretised) WLC model. Importantly, we demonstrate that at the same level of discretisation, i.e., when each discretisation element is one basepair long, our algorithm gains about 5-6 orders of magnitude in the size of time steps over the inextensible WLC model. Further, we show that our model can be mapped one-on-one to a discretised version of the extensible WLC model; implying that the speed-up we achieve in our model must hold equally well for the latter. We also demonstrate the use of the method by simulating efficiently the tumbling behaviour of a dsDNA segment in a shear flow.
Electrostatic interactions between dielectric objects are complex and of a many-body nature, owing to induced surface bound charge. We present a collection of techniques to simulate dynamical dielectric objects. We calculate the surface bound charge
from a matrix equation using the Generalized Minimal Residue method (GMRES). Empirically, we find that GMRES converges very quickly. Indeed, our detailed analysis suggests that the relevant matrix has a very compact spectrum for all non-degenerate dielectric geometries. Each GMRES iteration can be evaluated using a fast Ewald solver with cost that scales linearly or near-linearly in the number of surface charge elements. We analyze several previously proposed methods for calculating the bound charge, and show that our approach compares favorably.
Multiblock copolymer chains in implicit nonselective solvents are studied by Monte Carlo method which employs a parallel tempering algorithm. Chains consisting of 120 $A$ and 120 $B$ monomers, arranged in three distinct microarchitectures: $(10-10)_{
12}$, $(6-6)_{20}$, and $(3-3)_{40}$, collapse to globular states upon cooling, as expected. By varying both the reduced temperature $T^*$ and compatibility between monomers $omega$, numerous intra-globular structures are obtained: diclusters (handshake, spiral, torus with a core, etc.), triclusters, and $n$-clusters with $n>3$ (lamellar and other), which are reminiscent of the block copolymer nanophases for spherically confined geometries. Phase diagrams for various chains in the $(T^*, omega)$-space are mapped. The structure factor $S(k)$, for a selected microarchitecture and $omega$, is calculated. Since $S(k)$ can be measured in scattering experiments, it can be used to relate simulation results to an experiment. Self-assembly in those systems is interpreted in term of competition between minimization of the interfacial area separating different types of monomers and minimization of contacts between chain and solvent. Finally, the relevance of this model to the protein folding is addressed.
A variety of electrostatic phenomena, including the structure of electric double layers and the aggregation of charged colloids and proteins, are affected by nonuniform electric permittivity. These effects are frequently ignored in analytical and com
putational studies, and particularly difficult to handle in situations where multiple dielectric contrasts are present, such as in colloids that are heterogeneous in permittivity. We present an extension to the Iterative Dielectric Solver developed by Barros and Luijten [Phys. Rev. Lett. 113, 017801 (2014)] that makes it possible to accurately compute the polarization of anisotropic particles with multiple dielectric contrasts. This efficient boundary-element method-based approach is applicable to geometries that are not amenable to other solvers, opening the possibility of studying collective phenomena of dielectrically anisotropic particles. We provide insight into the underlying physical reasons for this efficiency.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
