Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userLattice-switch Monte Carlo for binary hard-sphere crystals
187
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We show how to generalize the Lattice Switch Monte Carlo method to calculate the phase diagram of a binary system. A global coordinate transformation is combined with a modification of particle diameters, enabling the multi-component system in question to be explored and directly compared to a suitable reference state in a single Monte Carlo simulation. We use the method to evaluate the free energies of binary hard sphere crystals. Calculations at moderate size ratios, alpha=0.58 and alpha=0.73, are in agreement with previous results, and confirm AB2 and AB13 as stable structures. We also find that the AB(CsCl) structure is not entropically stable at the size ratio and volume at which it has been reported experimentally, and therefore that those observations cannot be explained by packing effects alone.
rate research
Read More
We present a Monte Carlo method for the direct evaluation of the difference between the free energies of two crystal structures. The method is built on a lattice-switch transformation that maps a configuration of one structure onto a candidate configuration of the other by `switching one set of lattice vectors for the other, while keeping the displacements with respect to the lattice sites constant. The sampling of the displacement configurations is biased, multicanonically, to favor paths leading to `gateway arrangements for which the Monte Carlo switch to the candidate configuration will be accepted. The configurations of both structures can then be efficiently sampled in a single process, and the difference between their free energies evaluated from their measured probabilities. We explore and exploit the method in the context of extensive studies of systems of hard spheres. We show that the efficiency of the method is controlled by the extent to which the switch conserves correlated microstructure. We also show how, microscopically, the procedure works: the system finds gateway arrangements which fulfill the sampling bias intelligently. We establish, with high precision, the differences between the free energies of the two close packed structures (fcc and hcp) in both the constant density and the constant pressure ensembles.
We present a method for the direct evaluation of the difference between the free energies of two crystalline structures, of different symmetry. The method rests on a Monte Carlo procedure which allows one to sample along a path, through atomic-displacement-space, leading from one structure to the other by way of an intervening transformation that switches one set of lattice vectors for another. The configurations of both structures can thus be sampled within a single Monte Carlo process, and the difference between their free energies evaluated directly from the ratio of the measured probabilities of each. The method is used to determine the difference between the free energies of the fcc and hcp crystalline phases of a system of hard spheres.
We describe a Monte Carlo procedure which allows sampling of the disjoint configuration spaces associated with crystalline and fluid phases, within a single simulation. The method utilises biased sampling techniques to enhance the probabilities of gateway states (in each phase) which are such that a global switch (to the other phase) can be implemented. Equilibrium freezing-point parameters can be determined directly; statistical uncertainties prescribed transparently; and finite-size effects quantified systematically. The method is potentially quite general; we apply it to the freezing of hard spheres.
We construct asymptotic arguments for the relative efficiency of rejection-free Monte Carlo (MC) methods compared to the standard MC method. We find that the efficiency is proportional to $exp{({const} beta)}$ in the Ising, $sqrt{beta}$ in the classical XY, and $beta$ in the classical Heisenberg spin systems with inverse temperature $beta$, regardless of the dimension. The efficiency in hard particle systems is also obtained, and found to be proportional to $(rhoc -rho)^{-d}$ with the closest packing density $rhoc$, density $rho$, and dimension $d$ of the systems. We construct and implement a rejection-free Monte Carlo method for the hard-disk system. The RFMC has a greater computational efficiency at high densities, and the density dependence of the efficiency is as predicted by our arguments.
The Cooperative Motion Algorithm is an efficient lattice method to simulate dense polymer systems and is often used with two different criteria to generate a Markov chain in the configuration space. While the first method is the well-established Metropolis algorithm, the other one is an heuristic algorithm which needs justification. As an introductory step towards justification for the 3D lattice polymers, we study a simple system which is the binary equimolar uid on a 2D triangular lattice. Since all lattice sites are occupied only selected type of motions are considered, such the vacancy movements, swapping neighboring lattice sites (Kawasaki dynamics) and cooperative loops. We compare both methods, calculating the energy as well as heat capacity as a function of temperature. The critical temperature, which was determined using the Binder cumulant, was the same for all methods with the simulation accuracy and in agreement with the exact critical temperature for the Ising model on the 2D triangular lattice. In order to achieve reliable results at low temperatures we employ the parallel tempering algorithm which enables simultaneous simulations of replicas of the system in a wide range of temperatures.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
