Monte Carlo Simulations of the Yukawa One-Component Plasma

The excess free energy f of the Yukawa one-component plasma is investigated by means of Monte Carlo simulations. These simulations are performed in the canonical ensemble within hyperspherical boundary conditions and f is computed for various values of the coupling parameter in the range 0.1100 and of the screening parameter * in the range 0.1*6.     

Partial auxeticity induced by nanoslits in the Yukawa crystal

Monte Carlo simulations in the isobaric–isothermal ensemble with variable shape of the periodic box were used to investigate Poisson's ratios of the hard‐core repulsive Yukawa crystals with periodically distributed nanoslits. Each nanoslit, oriented perpendicularly to the crystallographic direction [010], was filled by a monolayer of hard spheres. It is shown that presence of the nanoslits leads to negative Poisson's ratio, as low as –0.57(2), in the [100][001]‐direction. This direction is not auxetic in the pure Yukawa crystal, i.e. shows positive Poisson's ratio, equal to 0.43(1).     

Molecular dynamics study of two-and three-dimensional classical fluids using double Yukawa potential

We have carried out a molecular dynamics simulation of two- and three-dimensional double Yukawa fluids near the triple point. We have compared some of the static and dynamic correlation functions with those of Lennard—Jones, when parameters occurring in double Yukawa potential are chosen to fit Lennard-Jones potential. The results are in good agreement. However, when repulsive and attractive parameters occurring in double Yukawa potential are varied, we found distinct differences in static and dynamic correlation functions. We have also compared the two-dimensional correlation functions with those of three-dimensional to study the effect of dimensionality, near the triple point region.     

Bound states and critical behavior of the Yukawa potential

We investigate the bound states of the Yukawa potential V (r)=−λexp(−αr)/r, using different algorithms: solving the Schr?dinger equation numerically and our Monte Carlo Hamiltonian approach. There is a critical α = α_C, above which no bound state exists. We study the relation between α_C and λ for various angular momentum quantum number l, and find in atomic units, α_C(l) = λ[A ₁ exp(−l/B ₁) + A ₂ exp(−l/B ₂)], with A ₁ = 1.020(18), B ₁ = 0.443(14), A ₂ = 0.170(17), and B ₂ = 2.490(180).     

Lennard-Jones binary fluids: A comparative study between the molecular dynamics and Monte Carlo descriptions of their structural properties

Molecular dynamics and Monte Carlo techniques are employed for the study of binary Lennard-Jones fluids. Systematic comparisons between the predictions of both techniques are discussed, with particular emphasis on the dependency of the structural properties with respect to temperature and Lennard-Jones potential parameters.     

An efficientmethod to determine chemical potential of mixtures in the isothermal and isobaric bulk phase with kineticMonte Carlo simulation

A kinetic Monte Carlo (kMC) scheme in the NPT ensemble (constant number of molecules, pressure and temperature) has been developed to determine accurate chemical potentials for all components in a homogeneous mixture. The simulation requires two moves: (1) a displacement move and (2) a volume change move. In the former, the mobility rate of a selected molecule is determined by its interaction with all the other molecules in the system and is moved to a random position within the simulation box, according to the Rosenbluth algorithm, without any rejections (entropic sampling). The volume change move is decided by a comparison between either the instant pressure or the partial average pressure (with long-range correction) and the specified pressure and is carried out much less frequently than the displacement move. We applied this NPT scheme to a number of mixtures in both the gaseous and liquid phases, and show that the derived chemical potentials are accurate and reproducible. The method is recommended for obtaining chemical potentials in mixtures that are required as input in a grand canonical ensemble simulation.     

A three-component model on the structure of colloidal solution with size-asymmetric electrolytes

ABSTRACT

A three-component model on the structure of colloidal solution with size asymmetric electrolytes is attempted here using density functional theory and Monte Carlo simulation. The solvent is represented as an individual component along with that as a dielectric continuum. The theory uses a weighted density approximation for the hard-sphere contribution to the free energy, whereas the ionic contribution is evaluated through a perturbation expansion around the bulk density. The theory is found to reproduce the simulation data quite well for a wide range of parametric conditions. The present study reflects the importance of the presence of the solvent in determining the structural behaviour of spherical double layers.     

Effect of ionic size on the structure of spherical double layers: a Monte Carlo simulation and density functional theory study

The effect of ionic size on the diffuse layer characteristics of a spherical double layer is studied using Monte Carlo simulation and density functional theory within the restricted primitive model. The macroion is modelled as an impenetrable charged hard sphere carrying a uniform surface charge density, surrounded by the small ions represented as charged hard spheres and the solvent is taken as a dielectric continuum. The density functional theory uses a partially perturbative scheme, where the hard sphere contribution to the one particle correlation function is evaluated using weighted density approximation and the ionic interactions are calculated using a second-order functional Taylor expansion with respect to a bulk electrolyte. The Monte Carlo simulations have been performed in the canonical ensemble. The detailed comparison is made in terms of zeta potentials for a wide range of physical conditions including different ionic diameters. The zeta potentials show a maximum or a minimum with respect to the polyion surface charge density for a divalent counterion. The ionic distribution profiles show considerable variations with the concentration of the electrolyte, the valency of the ions constituting the electrolyte, and the ionic size. This model study shows clear manipulations of ionic size and charge correlations in dictating the overall structure of the diffuse layer.     

Monte Carlo simulation of overdamped motion in a double well potential under the influence of coloured noise

We present a Monte Carlo simulation algorithm for evaluating the stationary probability and the mean and the variance of first passage times in any dynamical system under the influence of additive coloured Gaussian and Marcovian noise (Ornstein-Uhlenbeck process). Our algorithm generates the Ornstein-Uhlenbeck process by a superposition of a finite number of random telegraph processes. We obtain our results from a direct evaluation of the trajectories. We apply our method to the overdamped motion of a particle in a double well potential. We compare our simulation results with various analytic approximations for the stationary probability and the mean first passage times.