Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userGrand-canonical Monte-Carlo simulation of solutions of salt mixtures: theory and implementation
68
0
0.0
(
0
)
Authors
Toan T. Nguyen
Ask ChatGPT about the research
No Arabic abstract
A Grand-canonical Monte-Carlo simulation method extended to simulate a mixture of salts is presented. Due to charge neutrality requirement of electrolyte solutions, ions must be added to or removed from the system in groups. This leads to some complications compared to regular Grand Canonical simulation. Here, a recipe for simulation of electrolyte solution of salt mixture is presented. It is then implemented to simulate solution of 1:1, 2:1 and 2:2 salts or their mixtures at different concentrations using the primitive ion model. The osmotic pressures of the electrolyte solutions are calculated and shown to depend linearly on the salt concentrations within the concentration range simulated. We also show that at the same concentration of divalent anions, the presence of divalent cations make it easier to insert monovalent cations into the system. This can explain some quantitative differences observed in experiments of the MgCl$_2$ salt mixture and MgSO$_4$ salt mixture.
rate research
Read More
We introduce a simple but efficient method for grand-canonical twist averaging in quantum Monte Carlo calculations. By evaluating the thermodynamic grand potential instead of the ground state total energy, we greatly reduce the sampling errors caused by twist-dependent fluctuations in the particle number. We apply this method to the electron gas and to metallic lithium, aluminum, and solid atomic hydrogen. We show that, even when using a small number of twists, grand-canonical twist averaging of the grand potential produces better estimates of ground state energies than the widely used canonical twist-averaging approach.
The large surface density changes associated with the (100) noble metals surface hex-reconstruction suggest the use of non-particle conserving simulation methods. We present an example of a surface Grand Canonical Monte Carlo applied to the transformation of a square non reconstructed surface to the hexagonally covered low temperature stable Au(100). On the other hand, classical Molecular Dynamics allows to investigate microscopic details of the reconstruction dynamics, and we show, as an example, retraction of a step and its interplay with the surface reconstruction/deconstruction mechanism.
We extensively investigated the critical behavior of mixtures of colloids and polymers via the two-component Asakura-Oosawa model and its reduction to a one-component colloidal fluid using accurate theoretical and simulation techniques. In particular the theoretical approach, hierarchical reference theory [Adv. Phys. 44, 211 (1995)], incorporates realistically the effects of long-range fluctuations on phase separation giving exponents which differ strongly from their mean-field values, and are in good agreement with those of the three-dimensional Ising model. Computer simulations combined with finite-size scaling analysis confirm the Ising universality and the accuracy of the theory, although some discrepancy in the location of the critical point between one-component and full-mixture description remains. To assess the limit of the pair-interaction description, we compare one-component and two-component results.
We show how canonical ensemble expectation values can be extracted from quantum Monte Carlo simulations in the grand canonical ensemble. In order to obtain results for all particle sectors, a modest number of grand canonical simulations must be performed, each at a different chemical potential. From the canonical ensemble results, grand canonical expectation values can be extracted as a continuous function of the chemical potential. Results are presented from the application of this method to the two-dimensional Hubbard model.
We have studied the electrostatic screening effect of NaCl solutions on the interactions between anionic lipid bilayers in the fluid lamellar phase using a Poisson-Boltzmann based mean-field approach with constant charge and constant potential limiting charge regulation boundary conditions. The full DLVO potential, including the electrostatic, hydration and van der Waals interactions, was coupled to thermal bending fluctuations of the membranes via a variational Gaussian Ansatz. This allowed us to analyze the coupling between the osmotic pressure and the fluctuation amplitudes and compare them both simultaneously with the measured dependence on the bilayer separation, determined by the small-angle X-ray scattering experiments. High-structural resolution analysis of the scattering data revealed no significant changes of membrane structure as a function of salt concentration. Parsimonious description of our results is consistent with the constant charge limit of the general charge regulation phenomenology, with fully dissociated lipid charge groups, together with a four-fold reduction of the membranes bending rigidity upon increasing NaCl concentration.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
