ترغب بنشر مسار تعليمي؟ اضغط هنا

An Accelerated Multiboson Algorithm for Coulomb Gases with Dynamical Dielectric Effects

43   0   0.0 ( 0 )
 نشر من قبل Robert Sedgewick
 تاريخ النشر 2006
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

A recent reformulation [1] of the problem of Coulomb gases in the presence of a dynamical dielectric medium showed that finite temperature simulations of such systems can be accomplished on the basis of completely local Hamiltonians on a spatial lattice by including additional bosonic fields. For large systems, the Monte Carlo algorithm proposed in Ref. [1] becomes inefficient due to a low acceptance rate for particle moves in a fixed background multiboson field. We show here how this problem can be circumvented by use of a coupled particle-multiboson update procedure that improves acceptance rates on large lattices by orders of magnitude. The method is tested on a one-component plasma with neutral dielectric particles for a variety of system sizes.



قيم البحث

اقرأ أيضاً

We discuss the application of the local lattice technique of Maggs and Rossetto to problems that involve the motion of objects with different dielectric constants than the background. In these systems the simulation method produces a spurious interac tion force which causes the particles to move in an unphysical manner. We show that this term can be removed using a variant of a method known from high-energy physics simulations, the multiboson method, and demonstrate the effectiveness of this corrective method on a system of neutral particles. We then apply our method to a one-component plasma to show the effect of the spurious interaction term on a charged system.
The net charge of solvated entities, ranging from polyelectrolytes and biomolecules to charged nanoparticles and membranes, depends on the local dissociation equilibrium of individual ionizable groups. Incorporation of this phenomenon, emph{charge re gulation}, in theoretical and computational models requires dynamic, configuration-dependent recalculation of surface charges and is therefore typically approximated by assuming constant net charge on particles. Various computational methods exist that address this. We present an alternative, particularly efficient charge regulation Monte Carlo method (CR-MC), which explicitly models the redistribution of individual charges and accurately samples the correct grand-canonical charge distribution. In addition, we provide an open-source implementation in the LAMMPS molecular dynamics (MD) simulation package, resulting in a hybrid MD/CR-MC simulation method. This implementation is designed to handle a wide range of implicit-solvent systems that model discreet ionizable groups or surface sites. The computational cost of the method scales linearly with the number of ionizable groups, thereby allowing accurate simulations of systems containing thousands of individual ionizable sites. By matter of illustration, we use the CR-MC method to quantify the effects of charge regulation on the nature of the polyelectrolyte coil--globule transition and on the effective interaction between oppositely charged nanoparticles.
114 - E. Rico , M. Dalmonte , P. Zoller 2018
An ab initio calculation of nuclear physics from Quantum Chromodynamics (QCD), the fundamental SU(3) gauge theory of the strong interaction, remains an outstanding challenge. Here, we discuss the emergence of key elements of nuclear physics using an SO(3) lattice gauge theory as a toy model for QCD. We show that this model is accessible to state-of-the-art quantum simulation experiments with ultracold atoms in an optical lattice. First, we demonstrate that our model shares characteristic many-body features with QCD, such as the spontaneous breakdown of chiral symmetry, its restoration at finite baryon density, as well as the existence of few-body bound states. Then we show that in the one-dimensional case, the dynamics in the gauge invariant sector can be encoded as a spin S=3/2 Heisenberg model, i.e., as quantum magnetism, which has a natural realization with bosonic mixtures in optical lattices, and thus sheds light on the connection between non-Abelian gauge theories and quantum magnetism.
430 - Yuwen Chen 2020
Derivative-free optimization (DFO) has recently gained a lot of momentum in machine learning, spawning interest in the community to design faster methods for problems where gradients are not accessible. While some attention has been given to the conc ept of acceleration in the DFO literature, existing stochastic algorithms for objective functions with a finite-sum structure have not been shown theoretically to achieve an accelerated rate of convergence. Algorithms that use acceleration in such a setting are prone to instabilities, making it difficult to reach convergence. In this work, we exploit the finite-sum structure of the objective in order to design a variance-reduced DFO algorithm that provably yields acceleration. We prove rates of convergence for both smooth convex and strongly-convex finite-sum objective functions. Finally, we validate our theoretical results empirically on several tasks and datasets.
Inspired by the allure of additive fabrication, we pose the problem of origami design from a new perspective: how can we grow a folded surface in three dimensions from a seed so that it is guaranteed to be isometric to the plane? We solve this proble m in two steps: by first identifying the geometric conditions for the compatible completion of two separate folds into a single developable four-fold vertex, and then showing how this foundation allows us to grow a geometrically compatible front at the boundary of a given folded seed. This yields a complete marching, or additive, algorithm for the inverse design of the complete space of developable quad origami patterns that can be folded from flat sheets. We illustrate the flexibility of our approach by growing ordered, disordered, straight and curved folded origami and fitting surfaces of given curvature with folded approximants. Overall, our simple shift in perspective from a global search to a local rule has the potential to transform origami-based meta-structure design.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا