Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userMonte Carlo simulation of electrostatic interactions in inhomogeneous dielectric media: Correct sampling for the local lattice simulation algorithm
219
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present a lattice Monte Carlo algorithm based on the one originally proposed by Maggs and Rossetto for simulating electrostatic interactions in inhomogeneous dielectric media. The original algorithm is known to produce attractive interactions between particles of the same dielectric constant in the medium of different dielectric constant. We demonstrate that such interactions are spurious, caused by incorrectly biased statistical weight arising from particle motion during the Monte Carlo moves. We propose a simple parallel tempering algorithm that corrects this unphysical bias. The efficacy of our algorithm is tested on a simple binary mixture and on an uncharged polymer in a solvent, and applied to salt-doped polymer solutions.
rate research
Read More
Comptonization is the process in which photon spectrum changes due to multiple Compton scatterings in the electronic plasma. It plays an important role in the spectral formation of astrophysical X-ray and gamma-ray sources. There are several intrinsic limitations for the analytical method in dealing with the Comptonization problem and Monte Carlo simulation is one of the few alternatives. We describe an efficient Monte Carlo method that can solve the Comptonization problem in a fully relativistic way. We expanded the method so that it is capable of simulating Comptonization in the media where electron density and temperature varies discontinuously from one region to the other and in the isothermal media where density varies continuously along photon paths. The algorithms are presented in detail to facilitate computer code implementation. We also present a few examples of its application to the astrophysical research.
Electrostatic correlations and variable permittivity of electrolytes are essential for exploring many chemical and physical properties of interfaces in aqueous solutions. We propose a continuum electrostatic model for the treatment of these effects in the framework of the self-consistent field theory. The model incorporates a space-or field-dependent dielectric permittivity and an excluded ion-size effect for the correlation energy. This results in a self-energy modified Poisson-Nernst-Planck or Poisson-Boltzmann equation together with state equations for the self energy and the dielectric function. We show that the ionic size is of significant importance in predicting a finite self energy for an ion in an inhomogeneous medium. Asymptotic approximation is proposed for the solution of a generalized Debye-Huckel equation, which has been shown to capture the ionic correlation and dielectric self energy. Through simulating ionic distribution surrounding a macroion, the modified self-consistent field model is shown to agree with particle-based Monte Carlo simulations. Numerical results for symmetric and asymmetric electrolytes demonstrate that the model is able to predict the charge inversion at high correlation regime in the presence of multivalent interfacial ions which is beyond the mean-field theory, and also show strong effect to double layer structure due to the space- or field-dependent dielectric permittivity.
A coarse-grained simulation model eliminates microscopic degrees of freedom and represents a polymer by a simplified structure. A priori, two classes of coarse-grained models may be distinguished: those which are designed for a specific polymer and reflect the underlying atomistic details to some extent, and those which retain only the most basic features of a polymer chain (chain connectivity, short-range excluded-volume interactions, etc.). In this review we mainly focus on the second class of generic polymer models, while the first class of specific coarse-grained models is only touched upon briefly.
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 interaction 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.
Monte Carlo simulations are widely used in many areas including particle accelerators. In this lecture, after a short introduction and reviewing of some statistical backgrounds, we will discuss methods such as direct inversion, rejection method, and Markov chain Monte Carlo to sample a probability distribution function, and methods for variance reduction to evaluate numerical integrals using the Monte Carlo simulation. We will also briefly introduce the quasi-Monte Carlo sampling at the end of this lecture.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
