Related papers: Widom insertion method in simulations with Ewald s…
The electrostatic potential and energy of point charges in a real crystal, in the presence of thermal vibrations, is obtained as a special case of the Fourier method. Incorporating the role of thermal vibrations in electrostatic energy…
An integral equation method is presented for the 1D steady-state Poisson-Nernst-Planck equations modeling ion transport through membrane channels. The differential equations are recast as integral equations using Green's 3rd identity…
Using multiple ion beam analysis measurements, or techniques, combined with self-consistent data processing, generally allows extracting more (or more accurate) information from the measurements than processing separately data from single…
The previously applied independent atom model (IAM) for highly charged ion-molecule collisions which implemented the suppression of multiple ionization and capture on the basis of geometric overlaps of cross-sectional areas representing…
Quantum simulation - the use of one quantum system to simulate a less controllable one - may provide an understanding of the many quantum systems which cannot be modeled using classical computers. Impressive progress on control and…
The constant potential molecular dynamics simulation method proposed by Siepmann and Sprik and reformulated later by Reed (SR-CPM) has been widely employed to investigate the metallic electrolyte/electrode interfaces, especially for…
The ion distribution around charged colloids in solution has been investigated intensely during the last decade. However, few theoretical approaches have included the influence of variation in the dielectric permittivity within the system,…
The hydration free energies of ions exhibit an approximately quadratic dependence on the ionic charge, as predicted by the Born model. We analyze this behavior using second-order perturbation theory. This provides effective methods to…
We consider diffusion of particles on a lattice in the so-called dynamical mean-field regime (memory effects are neglected). Interactions are local, that is, only among particles at the same lattice site. It is shown that a statistical…
Artifacts arise in the simulations of electrolytes using periodic boundary conditions (PBC). We show the origin of these artifacts are the periodic image charges and the constraint of charge neutrality inside the simulation box, both of…
Molecular dynamics simulations at a constant electric potential are an essential tool to study electrochemical processes, providing microscopic information on the structural, thermodynamic, and dynamical properties. Despite the numerous…
We present a method to facilitate Monte Carlo simulations in the grand canonical ensemble given a target mean particle number. The method imposes a fictitious dynamics on the chemical potential, to be run concurrently with the Monte Carlo…
We present a simple modification of the direct-forcing immersed boundary method (IBM) proposed by Uhlmann [J. Comput. Phys, 2005] in order to enable it to be applied to particulate flows with solid-to-fluid density ratios around unity. The…
Ewald summation and physically equivalent methods such as particle-mesh Ewald, kubic-harmonic expansions, or Lekner sums are commonly used to calculate long-range electrostatic interactions in computer simulations of polar and charged…
Two-dimensional ionic liquids with single site anion and cation-neutral dimer are studied by computer simulations and integral equation techniques, with the aim of characterizing differences with single site anion-cation mixtures, and also…
When evaluating the electrostatic potential, periodic boundary conditions in one, two or three of the spatial dimensions are often needed for different applications. The triply periodic Ewald summation formula is classical, and Ewald…
We study a model of dyons for SU(2) Yang-Mills theory at finite temperature T < T_c, in particular its ability to generate a confining force between a static quark antiquark pair. The interaction between dyons corresponds to a long-range…
A simple and efficient method for quantum Monte Carlo simulation is presented, based on discretization of the action in the path integral, and a Gaussian averaging of the potential, which works well e.g. with the Coulomb potential.
We utilize the domain integral equation formulation to simulate two-dimensional transverse electric scattering in a homogeneous medium and a summation of modulated Gaussian functions to approximate the dual Gabor window. Then we apply Ewald…
Within the reweighting approach, one has the freedom to choose the Monte Carlo action so that it provides a good overlap with the finite-\mu measure but remains simple to simulate. We explore several choices of action in the regime of small…