Related papers: Widom insertion method in simulations with Ewald s…
The calculation of chemical potential has traditionally been a challenge in atomistic simulations. One of the most used approaches is Widom's insertion method in which the chemical potential is calculated by periodically attempting to…
The numerical computation of chemical potential in dense, non-homogeneous fluids is a key problem in the study of confined fluids thermodynamics. To this day several methods have been proposed, however there is still need for a robust…
The Ewald summation technique is generalised to power-law 1/|r|^k potentials in three-, two- and one-dimensional geometries with explicit formulae for all the components of the sums. The cases of short-range, long-range and "marginal"…
The gradual insertion method for direct calculation of the chemical potential by molecular simulation is applied in the NpT ensemble to different quadrupolar two-centre Lennard-Jones fluids at high density state points. The results agree…
A modified 3D-Ewald summation is presented for accurately simulating the ion-dipole mixture under dielectric confinement. The method is based on the combination of image charges and image dipoles with the conventional Ewald summation and…
While the excess chemical potential is the key quantity in determining phase diagrams, its direct computation for high-density liquids of long polymer chains has posed a significant challenge. Computationally, the excess chemical potential…
We present a new approach to efficiently simulate electrolytes confined between infinite charged walls using a 3d Ewald summation method. The optimal performance is achieved by separating the electrostatic potential produced by the charged…
The evaluation of the electrostatic potential is fundamental to the study of condensed phase systems. We discuss the calculation of the relevant lattice summations by Ewald-type techniques. A model charge density is introduced, that cancels…
Computer simulations of model systems are widely used to explore striking phenomena in promising applications spanning from physics, chemistry, biology, to materials science and engineering. The long range electrostatic interactions between…
A crucial aspect in the simulation of electrochemical interfaces consists in treating the distribution of electronic charge of electrode materials that are put in contact with an electrolyte solution. Recently, it has been shown how a…
We derive analytic solutions for the potential and field in a one-dimensional system of masses or charges with periodic boundary conditions, in other words Ewald sums for one dimension. We also provide a set of tools for exploring the…
Using the specific model of a system of like charged ions confined between two planar like charged surfaces, we compare the predictions for the energy and density profile of four simulation methods available to treat the long range Coulomb…
We develop a computational method for modeling electrostatic interactions of arbitrarily-shaped, polarizable objects on colloidal length scales, including colloids/nanoparticles, polymers, and surfactants, dispersed in explicit ion…
To minimise systematic errors in Monte Carlo simulations of charged particles, long range electrostatic interactions have to be calculated accurately and efficiently. Standard approaches, such as Ewald summation or the naive application of…
An efficient real space method is derived for the evaluation of the Madelung's potential of ionic crystals. The proposed method is an extension of the Evjen's method. It takes advantage of a general analysis for the potential convergence in…
We present a simulation method that allows us to calculate the titration curves for systems undergoing protonation/deprotonation reactions -- such as charged colloidal suspensions with acidic/basic surface groups, polyelectrolytes,…
A physical model based on a Monte-Carlo approach is proposed to calculate the ionization dynam- ics of warm dense matters (WDM) within particle-in-cell simulations, and where the impact (col- lision) ionization (CI), electron-ion…
The charge neutralization of an ion beam by electron injection is investigated using a two-dimensional electrostatic particle-in-cell code. The simulation results show that electrostatic solitary waves (ESWs) can be robustly generated in…
Many popular methods for the calculation of chemical potentials rely on the insertion of test particles into the target system. In the case of liquids and liquid mixtures, this procedure increases in difficulty upon increasing density or…
We present a technique to efficiently compute long-range interactions in systems with periodic boundary conditions. We extend the well-known Ewald method by using a linear combination of screening Gaussian charge distributions instead of…