Modelling the structural evolution of ternary phosphate glasses from melts to solid amorphous materials

The local and medium-range structural properties of phosphate-based melts and glasses have been characterized by means of first principles (density functional theory) and classical (shell-model) molecular dynamics simulations. The structure of glasses with biomedically active molecular compositions, (P2O5)0.45(CaO)x(Na2O)0.55−x (x = 0.30, 0.35 and 0.40), have been generated using first principles molecular dynamics simulations for the full melt-and-quench procedure and the changes in the structural properties as the 3000 K melt is cooled down to room temperature have been compared extensively with those of the final glasses. The melts are characterized by a significant fraction of threefold (P3c) and fivefold (P5c) phosphorus atoms, but structural defects rapidly decrease during the cooling phase and for temperatures lower than 1800 K the system is free of under- and over-coordinated species. The analysis of the structures of the glasses at 300 K shows a prevalence of the metaphosphate Q2 and pyrophosphate Q1 species, whereas the number of Q3 units, which constitute the three-dimensional phosphate network, significantly decreases with the increase of calcium content in the glass. The radial and angular distribution functions indicate that higher calcium concentration in the glass leads to an increase of the rigidity of the phosphate tetrahedral network, which has been explained in terms of the calcium's higher field strength compared to that of sodium. The structural characterization of the melts and glasses obtained from first principles simulations was used to assess and validate a recently developed interatomic shell-model forcefield for phosphate-based materials. For all three compositions, our potential model is in good agreement with the first principles data. In the glass network, the forcefield provides a very good description of the split between the shorter distances of phosphorus to non-bonded oxygen and the longer distances of the phosphorus to bonded oxygen; the phosphorus–phosphorus medium-range distribution; and the coordination environment around the Na and Ca glass modifiers. Moreover, the distribution of the Qn species in the melts and glasses is in excellent agreement with the values extracted from the first principles simulations. In contrast, simulations using standard rigid ion potentials do not provide a satisfactory description of the local short-range structure of phosphate-based glasses and are therefore less suitable to model this class of multicomponent amorphous system

Modelling the effects of salt solutions on the hydration of calcium ions

Classical molecular dynamics simulations of several aqueous alkali halide salt solutions have been used to determine the effect of electrolytes on the structure of water and the hydration properties of calcium ions. Compared with the simulations of Ca2+ ions in pure liquid water, the frequency of water exchange in the first hydration shell of calcium, which is a fundamental process in controlling the reactivity of calcium(II) aqua-ions, is drastically reduced in the presence of other electrolytes in solution. The strong stabilization of the hydration shell of Ca2+ occurs not only when the halide anions are directly coordinated to calcium, but also when the alkali and halide ions are placed at or outside the second coordination shell of Ca2+, suggesting that the reactivity of the first solvation shell of the calcium ion can be influenced by the specific affinity of other ions in solution for the water molecules coordinated to Ca2+. Analysis of the hydrogen-bonded structure of water in the vicinity of the calcium ion shows that the average number of hydrogen bonds per water molecules, which is 1.8 in pure liquid water, decreases as the concentration of alkali–halide salts in solution increases, and that the temporal fluctuations of hydrogen bonds are significantly larger than those obtained for Ca2+ in pure liquid water. This effect has been explained in terms of the dynamics of reorganization of the O–H X (X = F, Cl and Br) hydrogen bond. This work shows the importance of solution composition in determining the hydrogen-bonding network and ligand-exchange dynamics around metal ions, both in solution and at the mineral–water interfaces, which in turn has implications for interactions occurring at the mineral–water interface, ultimately controlling the mobilization of ions in the environment as well as in industrial processes

Reconsidering Calcium Dehydration as the Rate-Determining Step in Calcium Mineral Growth

The dehydration of cations is generally accepted as the rate-limiting step in many processes. Molecular dynamics (MD) can be used to investigate the dynamics of water molecules around cations, and two different methods exist to obtain trajectory-based water dehydration frequencies. Here, these two different post-processing methods (direct method versus survival function) have been implemented to obtain calcium dehydration frequencies from a series of trajectories obtained using a range of accepted force fields. None of the method combinations reproduced the commonly accepted experimental water exchange frequency of 10–8.2 s–1. Instead, our results suggest much faster water dynamics, comparable with more accurate ab initio MD simulations and with experimental values obtained using neutron scattering techniques. We obtained the best agreement using the survival function method to characterize the water dynamics, and we show that different method combinations significantly affect the outcome. Our work strongly suggests that the fast water exchange kinetics around the calcium ions is not rate-limiting for reactions involving dissolved/solvated calcium. Our results further suggest that, for alkali and most of the earth alkali metals, mechanistic rate laws for growth, dissolution, and adsorption, which are based on the principle of rate-limiting cation dehydration, need careful reconsideration

Calcite surface structure and reactivity: molecular dynamics simulations and macroscopic surface modelling of the calcite-water interface

Calcite–water interactions are important not only in carbon sequestration and the global carbon cycle, but also in contaminant behaviour in calcite-bearing host rock and in many industrial applications. Here we quantify the effect of variations in surface structure on calcite surface reactivity. Firstly, we employ classical Molecular Dynamics simulations of calcite surfaces containing an etch pit and a growth terrace, to show that the local environment in water around structurally different surface sites is distinct. In addition to observing the expected formation of more calcium–water interactions and hydrogen-bonds at lower-coordinated sites, we also observed subtle differences in hydrogen bonding around acute versus obtuse edges and corners. We subsequently used this information to refine the protonation constants for the calcite surface sites, according to the Charge Distribution MUltiSite Ion Complexation (CD-MUSIC) approach. The subtle differences in hydrogen bonding translate into markedly different charging behaviour versus pH, in particular for acute versus obtuse corner sites. The results show quantitatively that calcite surface reactivity is directly related to surface topography. The information obtained in this study is not only crucial for the improvement of existing macroscopic surface models of the reactivity of calcite towards contaminants, but also improves our atomic-level understanding of mineral–water interactions

Molecular Dynamics Simulations of Hydroxyapatite Nanopores in Contact with Electrolyte Solutions: The Effect of Nanoconfinement and Solvated Ions on the Surface Reactivity and the Structural, Dynamical, and Vibrational Properties of Water

This research utilized Queen Mary’s MidPlus computational facilities, supported by QMUL Research-IT and funded by EPSRC grant EP/K000128/1. Via our membership of the UK’s HEC Materials Chemistry Consortium, which is funded by EPSRC (EP/L000202

Anisotropic diffusion of water molecules in hydroxyapatite nanopores

Funded by EPSRC Grant EP/K000128/1

Properties of water confined in hydroxyapatite nanopores as derived from molecular dynamics simulations

N.H. de Leeuw is grateful to ‘Université ParisEst Créteil’ (UPEC) for financial support received during the course of this research. T.T. Pham is grateful to the ‘Institut des sciences de l’ingénierie et des systèmes’ (INSIS) of the ‘Centre national de la recherche scientifique’ (CNRS) for financial support received during the course of this research. D. Di Tommaso would like to thank the Royal Society, UK, for the award of a Royal Society Industry Fellowship