Tzu-Jen
Lin
*a and
Cheng-Chau
Chiu
b
aDepartment of Chemical Engineering, Chung Yuan Christian University, 200 Chung Pei Road, Chung Li District, Taoyuan City, 32023, Taiwan. E-mail: tzujenlin999@gmail.com
bInstitute of Atomic and Molecular Sciences, Academia Sinica, No. 1, Roosevelt Rd., Sec. 4, Taipei, 10617, Taiwan
First published on 27th November 2017
Since the first detection of pre-nucleation clusters during the formation of calcium phosphate minerals, determining such clusters’ compositions and structures has become crucial for understanding the early-stage nucleation of these minerals in solutions. In previous experimental studies, the composition and sizes of pre-nucleation clusters have been calculated, but their structural information has been difficult to determine because they are very small (<1 nm). In this study, we examined the structures and infrared spectra of small- and medium-sized calcium phosphate clusters using ab initio calculations combined with implicit solvation models. Adding solvent effects increased the possibility of the existence of alternative configurations of calcium phosphate clusters other than their compact configurations. The calcium atoms had a tendency to be located outside of the clusters to coordinate with water molecules in the aqueous environment. The computed infrared spectra of extended small calcium phosphate clusters captured some of the features measured in the in situ infrared spectra, which supports the network structures proposed by large-scale molecular dynamics studies and X-ray adsorption near-edge spectra. The relative stabilities of medium-sized Ca9(PO4)6 clusters with respect to the stability of Posner's cluster in water were also reviewed. We found that in water, alternative structures with low symmetry or large dipole moments had lower energies than Posner's cluster.
In recent experiments, the appearance of stable pre-nucleation clusters in calcium phosphate systems has been observed. Wang10 observed stable nanosized clusters in solution without any organic template through atomic force microscopy (AFM), and Sommerdijk11 observed clusters with an average diameter of 0.87 ± 0.2 nm through high-resolution cryo-TEM. The cluster sizes in both experimental observations suggested that the observed clusters were Posner's clusters12 (Ca9(PO4)6), even though the structure details of these clusters could not be determined. In recent experimental studies, pre-nucleation calcium phosphate clusters in solution have been investigated. Using in situ cryo-TEM, DeYoreo and Sommerdijk13 detected calcium triphosphate clusters, [Ca2(HPO4)3]2−, and suggested that they were the building blocks for the formation of amorphous calcium phosphate. Qun Zhang14 used synchrotron X-ray adsorption near-edge spectroscopy (XANES) at the calcium K-edge spectrum to probe calcium phosphate clusters in the early stages of crystallization, and suggested that the most abundant component of the clusters in solution was Ca(η2-PO4)2L2 (L = H2O or η1-PO4). Qun Zhang proposed that these clusters would subsequently form network structures. Tian-Lan Zhang15 used in situ XANES, ex situ X-ray diffraction, continuous pH measurements, and coordination chemistry to propose an improved Posner's cluster, Ca9(PO4)6(H2O)30, in which the eight calcium atoms are coordinated with water molecules except that at the center of the cluster.
Determining the structure of pre-nucleation clusters in solution is highly challenging because their sizes are nearly or less than 1 nm. Theoretical calculations provide alternative methods to overcome this challenge. Treboux16–18 conducted pioneer studies on calcium phosphate clusters through theoretical calculations. He systematically examined calcium phosphate clusters such as Ca3(PO4)2 and Ca9(PO4)6 through ab initio calculations, and showed that the clusters preferred compact arrangements with D3h and S6 symmetry, respectively. The structure of the Ca9(PO4)6 cluster with S6 symmetry was close to that of Posner's cluster. However, Treboux's studies were performed in vacuo and did not consider solvent effects. Recent theoretical studies have begun to include solvent effects. de Leeuw19 used ab initio molecular dynamics (AIMD) combined with umbrella sampling to determine the stability of calcium phosphate clusters in an aqueous environment. The AIMD study demonstrated that the structure of the pre-nucleation calcium phosphate cluster was Ca(η2-HPO42−)2(η1-HPO42−)(H2O)2, in which the calcium ion was sevenfold-coordinated. Furthermore, de Leeuw20 employed classical molecular dynamics to study the aggregation of calcium and phosphate ions in water, and the results showed that these ions preferred to aggregate into a structure similar to Posner's cluster.
With such diverse information available in the literature, it would not be trivial to investigate the stable structures of calcium phosphate clusters in an aqueous environment by using affordable simulation methods. In this study, we selected ab initio methods rather than classical molecular dynamics to evaluate the structure of calcium phosphate clusters because the performance and accuracy of the force field parameters of mineral ions is still uncertain.5,21,22 In addition, force field parameters usually have to be tuned for specific chemical and physical conditions, which could lead to case-dependent parameters. AIMD can include solvent effects, but is computationally demanding, thus leading to limited system sizes and sampling times. Therefore, we employed compromise methods, ab initio calculations combined with implicit solvation models, for our study. We compared optimized structures of small calcium phosphate clusters in vacuo and in water, and then re-examined the relative stability of medium-sized Ca9(PO4)6 clusters with respect to the stability of Posner's cluster in vacuo and in water. Furthermore, we computed the infrared spectra of calcium phosphate clusters. Previous theoretical studies have focused on the infrared spectra of hydroxyapatite or carbonate-substituted hydroxyapatite crystals.23,24 Boldeskul obtained the infrared spectra of calcium phosphate clusters, but the structure information of the calcium phosphate clusters in that study was elusive.25 In this study, we attempted to correlate the computed infrared spectra of calcium phosphate clusters with their in situ Fourier transform infrared (FTIR) spectra.13
The optimized structures and relative energies of calcium phosphate clusters were evaluated by density functional theory (DFT).26,27 Calcium phosphate clusters are highly polar and contain big phosphate anions. The meta GGA hybrid M06-2X functional shows good performance not only for main group thermochemistry, but also for anionic systems.28,29 6-31+G(d) basis sets were implemented for geometry optimizations, and this was followed by single-point calculations conducted using the same functional with 6-311+G(d) basis sets. Solvent effects were included by employing an SMD solvation model.30 The SMD solvation model is an extension of the SM8 (The 8 generation of universal Solvation Model) solvation model which was parameterized by a large training set and developed by Cramer and Truhlar.31 The “D” in SMD stands for “density”, which means the solvation model is based on the polarized continuous quantum mechanical charge density of the solute. The performance of the SMD solvation model has been tested for a broad range of systems.32,33 The infrared spectra of calcium phosphate clusters in water were on the basis of harmonic approximation, and were computed by using M062X/6-31+G(d) with the SMD solvation model. All DFT calculations were conducted by using Gaussian 09, Revision E.01.34 In this study, no explicit water molecules were considered because the number and placements of water molecules would have made the simulations unmanageable and beyond the capability of the computational resources in our lab.
Environment | [Ca(PO4)]1− | [Ca(PO4)2]4− | [Ca2(PO4)]1+ | |||
---|---|---|---|---|---|---|
Vacuo | Water | Vacuo | Water | Vacuo | Water | |
Ca–O distance (Å) | 2.251 | 2.418 | 2.509 | 2.443 |
2.400
2.307 |
2.446 |
Calcium coordinated
P–O bond lengths (Å) |
1.621 | 1.580 | 1.588 | 1.578 |
1.523
1.634 |
1.565 |
Calcium uncoordinated
P–O bond lengths (Å) |
1.499 | 1.555 | 1.533 | 1.558 | N.A. | N.A. |
O–P–O angles (°) |
100.0
117.8 |
109.5 |
102.9
115.5 |
109.5 | 93–128 | 109.5 |
Ca–O coordination number | 3 | 2 | 6 | 4 | 3, 3 | 2, 2 |
We then investigated the optimized structures of [Ca2(PO4)2]2− [Ca2(PO4)3]5− and Ca3(PO4)2 clusters in water. We constructed two configurations for each cluster. One was a compact configuration, and the other was an extended configuration (see Fig. 2). The compact configuration was constructed by maximizing Ca–O coordination, thus the calcium atoms were placed between phosphate ions. The extended configuration was constructed to verify the experimental conjecture that pre-nucleation clusters would aggregate into network structures or liquid-like ionic polymers.5,14Table 2 shows that the Ca–O coordination numbers in the compact configurations of [Ca2(PO4)2]2− [Ca2(PO4)3]5− and Ca3(PO4)2 clusters are higher than those of the extended configurations. In the compact configuration of [Ca2(PO4)2]2− and Ca3(PO4)2 clusters, the Ca–O distances were converged to 2.466 Å and 2.493 Å, respectively. This indicated that the strengths of Ca–O interactions in these two clusters were consistent. By contrast, the variation of Ca–O distances in the extended configurations of [Ca2(PO4)2]2− and Ca3(PO4)2 clusters indicated the different strengths of Ca–O interactions in these two clusters. In the [Ca2(PO4)3]5− cluster, different strengths of Ca–O interactions were in both compact and extended configuration based on the Ca–O distances ranging from 2.439 Å to 2.497 Å. The calcium-uncoordinated P–O bond lengths of [Ca2(PO4)2]2−, [Ca2(PO4)3]5−, and Ca3(PO4)2 clusters in water, no matter whether in the compact or extended configuration, were 1.53–1.55 Å, which were much shorter than the P–O bond lengths of a free phosphate ion in water. Based on DFT calculations combined with the implicit solvation model, we discovered that the compact configurations of [Ca2(PO4)2]2−, [Ca2(PO4)3]5−, and Ca3(PO4)2 clusters in water all had lower energies than their extended counterparts. The close-packed configurations of [Ca2(PO4)2]2− and [Ca2(PO4)3]5− clusters had 15 kcal mol−1 lower energies than their extended configurations, and the close-packed configuration of the Ca3(PO4)2 cluster had 25 kcal mol−1 lower energies than its extended configurations. This was mainly because of the higher Ca–O coordination number of the compact configurations than those of the extended configurations. The extended configurations gave much more space for coordination with water molecules in an aqueous environment than the close-packed configurations, but they also had less Ca–O bond interactions inside the clusters at the same time. As the calcium phosphate clusters grew, the Ca–O coordination number in the cluster seemed to be a dominant factor for creating stable structures. Although, the extended configurations had higher energies than the compact configurations, it was still difficult to deny their existence in water. We later found that the extended configurations captured some features of the measured in situ FTIR spectra of the clusters. The relative energies of the extended and compact configurations of [Ca2(PO4)2]2− [Ca2(PO4)3]5− and Ca3(PO4)2 clusters may need to be determined by AIMD where the clusters are solvating with numbers of explicit water molecules.
Configuration | [Ca2(PO4)2]2− | [Ca2(PO4)3]5− | Ca3(PO4)2 | |||
---|---|---|---|---|---|---|
Compact | Extended | Compact | Extended | Compact | Extended | |
Ca–O distance | 2.466 |
2.439
2.444 2.467 |
2.445–2.497 | 2.439–2.474 | 2.493 |
2.441
2.462 |
Ca-coordinated
Ca–O bond |
1.567
1.586 |
1.565
1.579 |
1.565–1.585 | 1.565–1.587 | 1.575 | 1.565 |
Ca-uncoordinated Ca–O bond | 1.546 | 1.558 |
1.547
1.559 |
1.546
1.557 |
1.535 | N.A. |
Ca–O coordination number | 4, 4 | 2, 4 | 4, 6 | 4, 4 | 4, 4, 4 | 2, 4, 2 |
Relative energy | 0 | 16.1 | 0 | 13.5 | 0 | 24.6 |
Fig. 3 The computed infrared spectra of (a) [Ca(PO4)]1−, (b) [Ca(PO4)2]4−, and (c) [Ca2(PO4)]1+ clusters in water. The blue curves represent the computed spectrum of a free phosphate ion in water. |
The computed infrared spectra of the extended and compact configurations of [Ca2(PO4)2]2−, [Ca2(PO4)3]5−, and [Ca3(PO4)2] clusters are presented in Fig. 4. The profile of infrared spectra of the compact configurations exhibited sharp peaks, and the profile of the extended configurations showed a distribution character especially in the extended [Ca2(PO4)2]2− cluster. The calcium-coordinated P–O stretching frequencies were mostly located between 900 cm−1 and 1025 cm−1. The calcium uncoordinated P–O stretching frequencies were located at 1050 cm−1, and the stretching frequency shifts higher at 1085 cm−1 in the case of the compact Ca3(PO4)2 cluster. As mentioned in the previous section, the Ca–O interactions in the compact configurations are more consistent than those of the extended configurations based on their Ca–O distances. This lead to the sharp P–O stretching peaks in the infrared spectra of the compact configurations. By contrast, the P–O bonds experienced different strengths of Ca–O interactions in the extended configurations, leading to the distributed P–O stretching peaks in the infrared spectra. The in situ FTIR spectrum at an early stage of calcium phosphate nucleation is given in Fig. 4(d).13 The P–O stretching mode corresponds to a broad distribution between 950 cm−1 and 1150 cm−1. Based on the computed infrared spectra, P–O bonds of the phosphate ions in the calcium phosphate clusters may feel different strengths of Ca–O interactions at the early stage of nucleation, and the configuration of calcium phosphate clusters was extended rather than compact. We already knew that from an energy perspective, the extended configurations of [Ca2(PO4)2]2−, [Ca2(PO4)3]5− and Ca3(PO4)2 clusters are not as stable as their compact counterparts. Nevertheless, by comparing the computed and in situ infrared spectra, we cannot rule out the possibility of their existence in water. As mentioned in the previous section, the relative stability of these compact and extended configurations should be further examined by using AIMD simulations.
Fig. 4 The computed infrared spectra of (a) [Ca2(PO4)2]2−, (b) [Ca2(PO4)3]5−, and (c) Ca3(PO4)2 clusters in water. (d) The sketched in situ FTIR spectrum based on ref. 13. |
We then conducted geometry optimizations for the 39 Ca9(PO4)6 clusters in combination with the SMD solvation model, and the relative energies and the change of dipole moments in vacuo and in water are presented in Fig. 5(c) and (d), respectively. Surprisingly, only five structures had relative energies higher than 20 kcal mol−1 with respect to the S6 isomer. Solvent effects lowered the relative energies of the low symmetry structures, and these structures generally exhibited a considerable dipole moment change in water compared with in vacuo. Moreover, solvent effects greatly reduced the relative energies of Ca9(PO4)6 clusters starting in high-symmetry configurations. Most of these optimized structures did not maintain their starting symmetry after geometry optimizations in water. Td, D3d-1, and CS-14 clusters were the most stable clusters with respect to the S6 cluster in water, and the relative energies were −6.7, −8.7, and −9.8 kcal mol−1, respectively. Stott relaxed calcium phosphate clusters which were extracted from hydroxyapatite crystals without symmetry restrictions, and discovered that the structure has lower energy than Posner's cluster by about 1.3 eV using HF and B3LYP in vacuo.36 We examined the structure by the theoretical method used in this study, and the relative energy with respect to the S6 isomer was only −2.6 kcal mol−1 in water.
T d and D3d-1 structures were notable, and their structures and structural information in vacuo and in water are presented in Fig. 6 and Table 3. The relative energies of the Td and D3d-1 clusters were 387 and 162 kcal mol−1 with respect to the S6 isomer, respectively, without considering solvent effects, but −6.7 and −8.7 kcal mol−1, respectively, when solvent effects were included. The optimized structures of Td and D3d-1 had zero dipole moments in vacuo, but had dipole moments of 10 and 31 D, respectively, in water. The considerable dipole moment of the D3d-1 cluster in water could be attributed to its optimized structure in water. Fig. 6 displays fully relaxed Td, D3d-1, S6, and CS-14 clusters in vacuo and in water, and the optimized structure of Td and D3d-1 clusters in water were very different from their optimized structures in vacuo. In vacuo, the optimized structure of Td and D3d-1 clusters maintained their highly symmetric structures, but their structures became relaxed in water shown in Fig. 6. The histogram of Ca–O distance in Td and D3d-1 clusters in vacuo and in water also showed an obvious demonstration of the structural character of these two clusters (see Fig. 7). The histograms of Ca–O distance in Td and D3d-1 clusters showed a sharp distribution in vacuo rather than a broad distribution in water. This indicated that considerable structural changes lead to large changes in dipole moment of these two structures in vacuo and in water. By contrast, the changes of the Ca–O distance histogram in S6 and CS-14 clusters between in vacuo and in water were not as obvious as the change in Td and D3d-1 clusters. The small changes in the structures of S6 and CS-14 clusters lead to the small differences of dipole moments in vacuo and in water. In the case of D3d-1, the relaxed structure in water resembles a pyramid. Two of the six phosphate ions are displaced to the lower part of the cluster. Therefore, the distribution of phosphate ions in the upper and lower parts of the D3d-1 cluster was 2-to-4. The displaced phosphate ions create an imbalanced charge distribution in the cluster, leading to an enormous dipole moment. This polar molecule was, however, stable in an aqueous environment, which is consistent with the like-dissolves-like principle.
Cluster | Relative energy (kcal mol−1) | Dipole moment (Debye) | Ca–Ca distances of hexahedron (Å) | |||
---|---|---|---|---|---|---|
In vacuo | Water | In vacuo | Water | In vacuo | Water | |
T d | 387 | −6.7 | 0 | 10 | 4.28 | 4.42 ± 0.45 |
D 3d-1 | 162 | −8.7 | 0 | 31.3 | 3.58, 4.13 | 4.43 ± 0.45 |
S 6 | 0 | 0 | 0 | 0 | 3.84, 4.40 | 4.33 ± 0.02 |
C S-14 | −4.7 | −9.8 | 1.9 | 6 | 4.05 ± 0.35 | 4.30 ± 0.32 |
Fig. 7 Ca–O distance distributions in vacuo and in water of (a) Td, (b) D3d-1, (c) S6, and (d) CS-14 clusters. The range of P–O bond lengths was from 2 to 9 Å with an increment of 0.5 Å. |
T d and D3d-1 clusters exhibited large Ca–Ca side length variations of the calcium hexahedron being in vacuo and in water. The Ca–Ca side lengths of the calcium hexahedron in the D3d-1 cluster in vacuo were only 4.13 Å and 3.58 Å listed in Table 3. In water, these lengths were increased to 4.43 ± 0.45 Å, with the largest increase being to 5.2 Å. For the Td cluster, the Ca–Ca side lengths were 4.28 Å in vacuo and 3.87–5.32 Å in water. The increase in the Ca–Ca side lengths from vacuum to water was also observed in S6 and CS-14 clusters, with increases of approximately 0.3–0.7 Å, which were not as apparent as those in the D3d-1 and Td clusters. The increase in the Ca–Ca side lengths of the calcium hexahedron in Ca9(PO4)6 clusters in water indicated that the calcium atoms had a tendency to be located at the surface of the clusters. This deduction is consistent with the simulation results for the small calcium phosphate clusters described in previous sections.
The computed infrared spectra of D3d-1, S6, and CS-14 clusters are shown in Fig. 8. The computed infrared spectra of D3d-1 and CS-14 clusters have different features than the spectrum of the S6 cluster. The S6 cluster exhibited three sharp peaks at 965 cm−1, 995 cm−1, and 1055 cm−1 in the infrared spectrum, and the D3d-1 and CS-14 clusters exhibited broad distributions between 960 cm−1 and 1085 cm−1 in the spectra. The sharp distribution of P–O bond lengths (see Fig. 9) of the S6 cluster in water indicated that the P–O bonds experienced consistent Ca–O interactions in the ordered cluster. Therefore, the computed infrared spectrum of the S6 cluster in water exhibited sharp peaks. By contrast, the P–O bond length had broad distributions in the less ordered D3d-1 and CS-14 clusters in water, which indicated that the P–O bonds experienced different strengths of Ca–O interactions. This lead to broad distributions of infrared spectra in these two clusters. The broad distribution in the in situ FTIR spectrum at an early stage of nucleation gradually transformed into a sharp peak at 1030 cm−1 as nucleation continued (Fig. 8(d)). This indicated that the some calcium phosphate clusters transformed from disordered to ordered structures during the nucleation process. Nevertheless, the measured spectrum also exhibited broad distributions at frequencies higher than 1030 cm−1 as nucleation continued, and not all calcium phosphate clusters aggregate into an ordered structure such as the S6 cluster. The broad distribution at frequencies higher than 1030 cm−1 indicated that the calcium phosphate clusters may aggregate into a less ordered structure such as D3d-1 and CS-14 clusters. The incoherence between the computed and measured IR spectra was the frequencies over 1050 cm−1 in D3d-1, S6, and CS-14 clusters. This was because the outermost P–O bonds of phosphate ions are always coordinated with ions in solutions and protonated in experiments, which was not considered in our simulation models. The coordination of ions or protonation of P–O bonds lowers the P–O stretching frequencies as mentioned in previous sections.
Fig. 8 The simulated vibration spectrum of (a) S6, (b) D3d-1, and (c) CS-14 Ca9(PO4)6 clusters in water. (d) The sketched in situ FTIR spectrum based on ref. 13. |
Fig. 9 P–O bond length distributions in water of (a) S6, (b) D3d-1, and (c) CS-14 clusters. The range of P–O bond lengths was from 1.5 to 1.6 Å with an increment of 0.01 Å. |
This study indicated that there could be many candidate structures with similar or even lower energies than Posner's cluster in an aqueous environment. This study was the initial search for calcium phosphate clusters in an aqueous environment. Ab initio methods combined with implicit solvation models appeared to provide an efficient approach to determining the structural features of mineral clusters in an aqueous environment, although the calculations of accurate relative energies may require the solvation of mineral clusters by explicit water molecules in theoretical calculations. Combining implicit solvation models and advanced algorithms such as ab initio random structure searching (AIRSS)37 would be worthwhile for conducting a detailed search for pre-nucleation calcium phosphate clusters in an aqueous environment in the future.
This journal is © the Owner Societies 2018 |