اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدEfficient and accurate simulation of dynamic dielectric objects
256
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
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.
قيم البحث
اقرأ أيضاً
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.
We have investigated thermodynamic and dynamic properties as well as the dielectric constant of water-metha-nol model mixtures in the entire range of composition by using constant pressure molecular dynamics simulations at ambient conditions. The SPC
/E and TIP4P/Ew water models are used in combination with the OPLS united atom modelling for methanol. Changes of the average number of hydrogen bonds between particles of different species and of the fractions of differently bonded molecules are put in correspondence with the behavior of excess mixing volume and enthalpy, of self-diffusion coefficients and rotational relaxation times. From the detailed analyses of the results obtained in this work, we conclude that an improvement of the description of an ample set of properties of water-methanol mixtures can possibly be reached, if a more sophisticated, carefully parameterized, e.g., all atom, model for methanol is used. Moreover, exploration of parametrization of the methanol force field, with simultaneous application of different combination rules for methanol-water cross interactions, is required.
We have employed molecular dynamics simulations based on the TIP4P/2005 water model to investigate the local structural, dynamical, and dielectric properties of the two recently reported body-centered-cubic and face-centered-cubic plastic crystal pha
ses of water. Our results reveal significant differences in the local orientational structure and rotational dynamics of water molecules for the two polymorphs. The probability distributions of trigonal and tetrahedral order parameters exhibit a multi-modal structure, implying the existence of significant local orientational heterogeneities, particularly in the face-centered-cubic phase. The calculated hydrogen bond statistics and dynamics provide further indications of the existence of a strongly heterogeneous and rapidly interconverting local orientational structural network in both polymorphs. We have observed a hindered molecular rotation, much more pronounced in the body-centered-cubic phase, which is reflected by the decay of the fourth-order Legendre reorientational correlation functions and angular Van Hove functions. Molecular rotation, however, is additionally hindered in the high-pressure liquid compared to the plastic crystal phase. The results obtained also reveal significant differences in the dielectric properties of the polymorphs due to the different dipolar orientational correlation characterizing each phase.
Because of consuming energy to drive their motion, systems of active colloids are intrinsically out of equilibrium. In the past decade, a variety of intriguing dynamic patterns have been observed in systems of active colloids, and they offer a new pl
atform for studying non-equilibrium physics, in which computer simulation and analytical theory have played an important role. Here we review the recent progress in understanding the dynamic assembly of active colloids by using numerical and analytical tools. We review the progress in understanding the motility induced phase separation in the past decade, followed by the discussion on the effect of shape anisotropy and hydrodynamics on the dynamic assembly of active colloids.
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.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
