A Monte Carlo method for grand canonical and grand isoshear ensemble simulations has been used to characterize the free energy, energy, and entropy of clay mineral swelling. The Monte Carlo approach was found to be more efficient at simulating water content fluctuations in the highly constrained clay environment than a previously developed molecular dynamics method. Swelling thermodynamics calculated for Cs-, Na-, and Sr-montmorillonite clays indicate a strong dependence of swelling on the interlayer ion identity, in agreement with various experimental measurements. The Sr clay swells most readily, and both the Na and Sr clays prefer expanded states (two-layer hydrate or greater) when in contact with bulk water. In contrast, swelling is inhibited in the Cs clay. Differences in swelling behavior are traced directly to the tendency of the different ions to hydrate. The swelling free energies are decomposed into their energetic and entropic components, revealing an overall energetic driving force for the swelling phenomena. Entropic effects provide a smaller, mediating role in the swelling processes. The results provide a unique molecular perspective on experimentally well-characterized swelling phenomena.
We study the problem of electron-ion temperature equilibration in plasmas. We consider pure H at various densities and temperatures, and Ar-doped H at temperatures high enough so that the Ar is fully ionized. Two theoretical approaches are used: classical molecular dynamics (MD) with statistical 2-body potentials, and a generalized Lenard-Balescu (GLB) theory capable of treating multi-component weakly-coupled plasmas. The GLB is used in two modes: 1) with the quantum dielectric response in the random-phase approximation (RPA) together with the pure Coulomb interaction, and 2) with the classical (h −→ 0) dielectric response (both with and without local-field corrections) together with the statistical potentials. We find that the MD results are described very well by classical GLB including the statistical potentials and without local-field corrections (RPA only); worse agreement is found when static local-field effects are included, in contradiction to the classical pure-Coulomb case with like charges. The results of the various approaches are all in excellent agreement with pure-Coulomb quantum GLB when the temperature is high enough.In addition, we show that classical calculations with statistical potentials derived from the exact quantum 2-body density matrix produce results in far better agreement with pure-Coulomb quantum GLB than classical calculations performed with older existing statistical potentials.
A new method for the determination of clay swelling thermodynamics from computer simulation is discussed and evaluated. This method allows for the determination of temperature, pressure, and water chemical potential dependence of clay swelling from simulations at a single thermodynamic state point. The temperature dependence and pressure dependence of clay swelling are shown to be directly related to the composite system entropy and volume change, respectively, that accompany swelling. Expressions for the chemical potential dependence of clay swelling are used to determine constant pressure layer spacing and adsorption isotherms, quantities that are well suited for comparison with experimental measurements. This method is evaluated through grand isoshear ensemble simulations of Na-montmorillonite, a prototypical swelling clay. Approximations associated with all expressions are discussed with explicit calculations used to demonstrate their regimes of validity.
We compute electrical and thermal conductivities of hydrogen plasmas in the nondegenerate regime using Kohn-Sham density functional theory (DFT) and an application of the Kubo-Greenwood response formula, and demonstrate that for thermal conductivity, the mean-field treatment of the electron-electron (e-e) interaction therein is insufficient to reproduce the weak-coupling limit obtained by plasma kinetic theories. An explicit e-e scattering correction to the DFT is posited by appealing to Matthiessen's Rule and the results of our computations of conductivities with the quantum Lenard-Balescu (QLB) equation. Further motivation of our correction is provided by an argument arising from the Zubarev quantum kinetic theory approach. Significant emphasis is placed on our efforts to produce properly converged results for plasma transport using Kohn-Sham DFT, so that an accurate assessment of the importance and efficacy of our e-e scattering corrections to the thermal conductivity can be made.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.