Atomic modelling of crystal/complex fluid/crystal contacts-Part I. The genetic iterative multi-species (GIMS) approach and case of kaolinite/brine/kaolinite
Résumé
Predicting the impact of underground engineering on the environment requires the knowledge of natural media at different scales. In particular, understanding basic phenomena controlling the properties of rocks in the presence of complex fluids necessitates a detailed atomic description of the solid/fluid/solid contacts, the subject of Part I of the present study. First, building the solid interspace between two different crystals in a non-periodic situation is achieved using the ab initio and molecular mechanics code GenMol (TM). A description of the fluid confined within the interspace is then derived from the original genetic iterative multi-species (GIMS) algorithm implemented in the same code. This approach consists of equilibrating chemical potentials, cycle after cycle and species after species, between the confined fluid and the free natural fluid. An elementary iteration for a species k consists of different steps incrementing the number N-k of particles k, with other numbers N-k remaining constant. At each step, an optimum fluid composition is obtained by a genetic process distributing the fluid particles on a grid by stochastic shots, followed in fine by a refining process. The effectiveness of the GIMS approach is demonstrated in the case study of a fluid confined between two (0 0 1) kaolinite faces, with apertures h varying between 4 and 10 A. connected to a 9-species external solution [H2O, Cl-, Na+ CO2(aq), NaCl(aq), Ca2+, Mg2+, HCO3-, H3O+] where concentrations are ranging from 55 to 10(-4) mol/L. The results show a drastic variation in the solute/solvent and cations/ions ratios in the confined fluid when aperture h is lowered to less than 1 nm. These results obtained with a very rapid convergence of the iterative algorithm combined with a very competitive genetic optimizer are validated with high precision on a free solution. This description of contacts between crystals is original and unattainable by standard crystal interface approaches. It opens the way to understanding and/or predicting phenomena between crystals in a complex natural environment, such as adhesion or repulsion investigated in the Part II of this study.