Lithium speciation in aqueous fluids at high P and T studied by ab initio molecular dynamics and consequences for Li-isotope fractionation between minerals and fluids

Ab initio molecular dynamics simulations are performed to study the speciation changes in lithium bearing aqueous fluids at high temperature (T = 1000 K) and high pressures, P, between about 0.3 and 6.0 GPa. The simulations show a linear increase in Li coordination with fluid density, from 3.2 to about five in the considered pressure range. Towards low densities, associated LiF complexes are becoming increasingly stable, which is quantified by evaluating the dynamic behavior of the respective species. In the high-density region, HF complexes are observed. The differences in speciation may be related to structural changes of the solvent under compression. At a fluid density of 1.2 g/cm3, kinks in the pressure dependences of the oxygen–oxygen nearest neighbor distance and the oxygen–oxygen coordination are observed, which indicates a change in compaction mechanism. Assuming that the Li coordination difference between crystal and fluid is a major determinant for the isotopic fractionation between minerals and fluids, we expect only a small pressure dependence of the Li isotopic fractionation between Li bearing fluids and minerals. Our simulation results are consistent with experimental data that show reverse fractionation of 7Li between fluid and mineral, when Li is in tetrahedral instead of octahedral coordination in the crystal.