Fan Yang,
Jinxuan Zhou,
Haiyang Huang and
Meng Li*
Guangxi Key Laboratory of Natural Polymer Chemistry and Physics, Nanning Normal University, Nanning 530001, PR China. E-mail: limeng_2016@126.com
First published on 4th August 2025
Phosphorus bioavailability in soils is largely governed by the nucleation and transformation of calcium phosphate (Ca-P) minerals, processes influenced by organic acids like citrate exuded by plant roots. Despite extensive research, the concentration-dependent effects of citrate on Ca-P crystallization remain unclear. This study explores how varying citrate concentrations regulate the nucleation, morphology, and phase transitions of Ca-P under weakly alkaline conditions. Real-time pH monitoring indicated that low citrate concentrations (≤1 μM) accelerate nucleation by promoting the transformation from amorphous calcium phosphate (ACP) to crystalline hydroxyapatite (HAP), whereas higher citrate levels (≥2 μM) inhibit nucleation, extending induction times markedly. Morphological analyses revealed distinct crystal shapes: flake-like structures at low citrate concentrations and dense spherical aggregates at higher levels. Raman spectroscopy and XPS characterization confirmed that citrate adsorbs onto Ca-P surfaces and modulates the electronic states of calcium and phosphorus. At low concentrations, citrate likely donates electron density to surface ions, facilitating nucleation, while at higher concentrations, electron-withdrawing carboxyl groups dominate, stabilizing intermediate amorphous phases and hindering crystallization. These findings elucidate the molecular mechanism behind citrate's dual regulatory role in ACP transformation, providing insights into phosphorus cycling in the rhizosphere and strategies for improving phosphate fertilizer efficiency in alkaline soils.
However, the efficiency of phosphorus fertilizer utilization remains low, typically ranging from 5% to 20%, due to factors such as surface runoff, soil erosion, and strong phosphorus fixation in soils.8–10 This inefficiency is particularly pronounced in alkaline calcareous soils, which typically exhibit pH values ranging from 7.5 to 8.5 due to the buffering capacity of carbonate minerals, where phosphorus fixation is intensified.11 In such soils, phosphorus undergoes an initial phase of adsorption and immobilization, followed by chemical fixation. The initially formed amorphous calcium phosphate (ACP) subsequently transforms into more crystalline phases such as calcium hydrogen phosphate dihydrate (CaHPO4·2H2O, DCPD) and octacalcium phosphate (Ca8H2(PO4)6·5H2O, OCP). Ultimately, these phases mature into hydroxyapatite (Ca10(PO4)6(OH)2, HAP), the most thermodynamically stable and least soluble form under neutral to alkaline conditions.12 This sequential crystallization significantly limits phosphorus bioavailability due to reduced solubility.
To mitigate the stress of phosphorus deficiency, plants have evolved adaptive strategies to enhance phosphorus mobilization in the rhizosphere. One prominent mechanism involves the secretion of low molecular weight organic acids, including citric, oxalic, tartaric, and malic acids,13,14 which can account for up to 20% of the carbon fixed by photosynthesis.15 Among these, citric acid is particularly effective in mobilizing phosphorus from soil minerals. Studies utilizing X-ray absorption near edge structure (XANES) spectroscopy reveal that soil phosphorus availability dramatically improves when citrate concentrations exceed approximately 10 mM.16 Although typical bulk soil solutions contain citrate concentrations below 100 μM,17 rhizosphere concentrations are substantially elevated, often by 10 to 40 times due to root exudation.18
The interaction between organic acids and calcium phosphate precipitates plays a critical role in regulating phosphorus bioavailability.19 Functional groups within the organic acids, especially carboxyl and hydroxyl groups, influence their capacity to solubilize calcium phosphate phases. For instance, organic acids with higher numbers of hydroxyl groups demonstrate greater solubility towards DCPD under alkaline conditions.20 Citrate, possessing one alcoholic hydroxyl and three carboxyl groups, is notably abundant in the plant rhizosphere. However, atomic force microscopy (AFM) studies on the role of citrate in phosphate dissolution have reported contradictory findings: citrate concentrations above 100 μM appear to promote dissolution, whereas intermediate concentrations (10–100 μM) may inhibit it.21 Additionally, citrate can facilitate Ca-P mineralization by reducing the nucleation energy barrier through interactions with interfacial Ca2+ ions or by displacing coordinated water molecules.22,23 Conversely, several studies have demonstrated that organic acids, including citrate, can inhibit calcium phosphate nucleation and crystallization.24 This suggests a complex, concentration-dependent regulatory mechanism that encompasses both promotive and inhibitory effects.
Given the paradoxical evidence regarding citrate's impact on calcium phosphate phase behavior—enhancing dissolution at higher concentrations while inhibiting nucleation or growth at lower levels—it is imperative to elucidate the underlying mechanisms governing these dual effects. We hypothesize that citrate modulates Ca-P crystallization in a concentration-dependent manner, exerting both promotive and inhibitory influences during different crystallization stages. To investigate this, we conducted real-time pH monitoring of solutions supersaturated with respect to hydroxyapatite (σHAP >0) under weakly alkaline conditions to track the kinetic progression of Ca-P phase transformation in the presence of varying citrate concentrations. Our approach enabled the capture and characterization of intermediate crystalline phases, revealing distinct morphological features—flaky structures under low citrate concentrations and spherical morphologies at higher concentrations. Notably, low citrate levels accelerated the phase transition from ACP to HAP, facilitating HAP nucleation. Although plant root exudation of citrate plays an important role in the rhizosphere, the current study focuses on model crystallization systems under controlled conditions. Our objective is to elucidate the fundamental role of citrate adsorption in the transformation of amorphous calcium phosphate, rather than directly modeling the complexity of soil or root environments. This study provides novel insights into the dualistic role of citrate in calcium phosphate nucleation and crystallization, advancing our understanding of phosphorus cycling in soils and optimizing fertilizer use efficiency in agricultural systems, particularly in calcareous soils where phosphorus fixation severely limits its accessibility.
![]() | ||
Fig. 1 (a) Effect of citrate concentration on calcium phosphate crystallization. (b) Statistical analysis of the induction time under different citrate concentrations. |
It indicates that citrate plays an important role in the regulation of calcium phosphate nucleation and morphology. To further confirm the effect of citrate on Ca-P morphology at different nucleation stages, TEM analyses were performed by dropping the reaction solution of three different stages (stage I: 10 minutes before rapid nucleation; stage II, pH = 6.9; stage III: reaction after 7 days) onto the copper network of the carbon-supported film. It showed that the spontaneous crystallization of Ca-P was initiated by a branched polymeric network with developed granules of ACP with dimensions of 30–50 nm (Fig. 3a). The ACP granules subsequently aggregated and developed into ribbons at stage II (Fig. 3b). Finally, the Ca-P crystals showed spherical clusters (Fig. S1a), in accordance with the spherical morphology observed by SEM (Fig. 2b), with lamellar crystals embedded after 7 days of reaction (Fig. 3c). Addition of citrate at the relatively low concentration of 0.5 μM showed a flocculent morphology of ACP (Fig. 3d), which might be the intermediate state of granules and ribbons morphology in Fig. 3a and b, considering the acceleration of Ca-P nucleation as shown in Fig. 1. These flocculent ACP transformed into needle-like microcrystals concurrently with the formation of ribbon crystals (Fig. 3e). The crystals finally transformed into short rod morphology after 7 days of reaction (Fig. 3f). With the increase of citrate to 50 μM, rapid crystallization proceeded with a few nodules with sizes of less than 30 nm (Fig. 3g). These nodules grew out to form a branched polymeric network (Fig. 3h) and subsequently aggregated into a spherical morphology (Fig. S1b), which was consistent with the morphology observed by SEM. The stage II crystals in the presence of 50 μM citrate did not show diffraction rings, which might be due to the relatively small crystal nuclei formed and not captured by selected area electron diffraction (SAED). The presence of both flaky crystals and short rod-like crystals was also noted (Fig. 3c, f and i) after 7 days of nucleation, which are identical to the crystal morphologies of OCP and HAP,32,33 respectively.
![]() | ||
Fig. 3 TEM images of calcium phosphate crystals formed under different concentrations of citrate. (a–c) 0 μM, (d–f) 0.5 μM, and (g–i) 50 μM citrate. |
![]() | ||
Fig. 4 HRTEM images of calcium phosphate crystals formed in the presence of different citric acid concentrations. Images correspond to (a–c) 0 μM, (d–f) 0.5 μM, and (g–i) 50 μM citrate. |
Ca-P crystallization undergoes a series of metastable mesophases, a phenomenon that is also observable during the crystallization of both calcium carbonate and calcium oxalate.36,37 Raman microscopy was further used to investigate the influence of citrate on the final crystallization products (Fig. 5). The analysis revealed that after 7 days of maturation, the calcium phosphate consisted of a mixture of OCP and HAP. Due to the structural similarities between OCP and HAP, many of their characteristic peaks overlap. Notably, common characteristic peaks at 577, 587, and 607 cm−1 correspond to v4 PO43− bending vibrations.38 Additionally, peaks at 960 cm−1 and 1043 cm−1 are attributed to the symmetric v1 PO43− stretching vibration and the asymmetric v3 PO43− stretching vibration, respectively. However, the v2 PO43− bending vibration peaks at 426 and 445 cm−1 is distinctive to that of HAP.37 In contrast, the v4 HPO42− bending vibration at 527 cm−1, the v1 HPO42− stretching vibrations at 878, 906, and 1000 cm−1, and the asymmetric v3 PO43− stretching vibration peak at 1073 cm−1 are characteristic peaks of OCP.39,40 Specifically, the peaks at 527 cm−1 and 878 cm−1 are attributed to the hydration layer of OCP, while the peak at 906 cm−1 is associated with the apatite layer of OCP.41 This differentiation of Raman peaks is essential for accurately characterizing the composition of calcium phosphate after citrate-influenced crystallization.
![]() | ||
Fig. 5 Raman spectra of calcium phosphate crystals matured for 7 days in the presence of different concentrations of citrate. |
According to the above Raman characteristic peak analysis, the presence of 0.5 μM citrate significantly enhances the distinct characteristic peaks associated with OCP and HAP compared to the absence of citrate. In contrast, when the citrate concentration is increased to 50 μM, only several major phosphate group peaks are prominent. This indicates that at a citrate concentration of 50 μM, calcium phosphate has not yet fully crystallized into the relatively stable OCP and HAP phases.
To elucidate the potential mechanisms by which citrate influences the crystallization process of HAP, X-ray photoelectron spectroscopy (XPS) was employed to analyze the chemical states of the constituent elements within the calcium phosphate mixtures. As depicted in Fig. 6a, the presence of Ca, P, and O elements was confirmed in all samples. A more detailed examination (Fig. 6b and c) reveals that in the absence of citrate, the binding energies at 347.34 eV and 133.18 eV correspond to the Ca 2p3/2 and P 2p3/2 peaks.42 The presence of 0.5 μM and 50 μM citrate resulted in noticeable shifts in the binding energies of both Ca 2p3/2 and P 2p3/2 peaks. Upon introduction of citrate at low concentration (0.5 μM), both signals exhibited substantial negative shifts: P 2p3/2 decreased by 0.31 eV to 132.87 eV, while Ca 2p3/2 shifted downward by 0.33 eV to 347.01 eV. This coordinated reduction in binding energies signifies enhanced electron density surrounding both Ca and P atoms at the mineral surface. Strikingly, this electronic effect reversed when citrate concentration was increased to 50 μM, where binding energies shifted positively: P 2p3/2 increased by 0.39 eV to 133.57 eV, and Ca 2p3/2 shifted upward by 0.08 eV to 347.42 eV. This inversion indicates significant electron depletion around these atomic centers at higher citrate loadings.43
![]() | ||
Fig. 6 XPS survey spectrum of (a) calcium phosphate matured for 7 days under different concentrations of citrate, and high-resolution XPS spectra of (b) P 2p and (c) Ca 2p. |
These concentration-dependent shifts in binding energies provide compelling evidence for distinct interaction mechanisms between citrate and the calcium phosphate surface. At low concentrations, weak coordination involving hydroxyl groups or possibly partial electron donation from carboxylate groups may locally increase electron density.44 Conversely, at higher concentrations, the collective electron-withdrawing effect of multiple adsorbed citrate molecules, particularly through the carboxylate groups, becomes dominant, leading to a decrease in electron density around the surface Ca and P atoms. This modification of the electronic environment of key structural elements on the calcium phosphate surface by adsorbed citrate is crucial. It directly influences surface charge, binding sites, and hydration layers, thereby fundamentally altering the nucleation and growth kinetics during HAP crystallization.45
The induction time (t) for the nucleation process is primarily governed by supersaturation (S) and interfacial energy (γ). As derived by Liu et al.,46 lower γ and higher S lead to shorter nucleation induction times. Meanwhile, under conditions of constant supersaturation, interfacial energy is often considered the primary variable influencing t; solely considering the influence of interfacial energy is clearly insufficient in the context of the present experiments.
Under the influence of a low citrate concentration (0.5 μM), the reduction in supersaturation is negligible, with only approximately 0.0375% of Ca2+ ions estimated to be complexed. A small amount of citrate adsorbed onto the calcium phosphate surface at this concentration is likely not sufficient to inhibit calcium phosphate crystallization. Instead, it is possible that this small amount of citrate, by complexing Ca2+ in solution (typically forming [Ca(Cit)]− and [Ca(HCit)]0, with logK values of −3.45 and −2.13 at 25 °C,47 respectively), can lead to localized increases in free Ca2+ concentration. This can potentially accelerate the initial nucleation process.
Conversely, at a high citrate concentration (50 μM), two main factors contribute to the observed effects. Firstly, a significant proportion of Ca2+ (estimated at 3.75%) is complexed by citrate, leading to a reduction in supersaturation.48 This reduction in supersaturation is reflected in the prolonged stabilization time of ACP, as the driving force for its transformation or further crystallization is diminished. Secondly, the high concentration of citrate in solution, coupled with the strong affinity of carboxylate groups (RCOO−) for Ca2+ ions,49 promotes significant adsorption of citrate onto the calcium phosphate surface.50 This adsorption forms a protective coating layer that effectively hinders both the initial nucleation of calcium phosphate and its subsequent crystalline growth.26
Besides, the electronic modification represents the molecular basis for citrate's observed effects on nucleation kinetics and crystallization pathways. At low concentrations (0.5–1 μM), limited citrate adsorption optimizes the electronic structure of nucleation sites, accelerating nucleation by 56% compared to controls. Conversely, at higher concentrations (2–50 μM), extensive citrate coverage depletes electron density around Ca and P atoms, inhibiting nucleation by extending induction time up to 356%. This concentration-dependent electronic modification manifests in distinct transformation pathways: ACP → OCP → HAP at low citrate concentrations versus ACP → DCPD → HAP at higher concentrations.
The electron density perturbation directly affects three critical aspects of calcium phosphate crystallization: (1) surface charge distribution, altering electrostatic interactions with ions and water molecules,51 (2) binding site geometry and reactivity, influencing ion incorporation rates,52,53 and (3) interfacial hydration layer structure, modifying water-mediated interactions during crystal growth.54 These effects collectively explain the observed morphological transitions from spherical aggregates to sheet-like structures as citrate concentration changes. The concentration gradient of citrate typically observed in soils—from bulk soil to rhizosphere—may create spatial heterogeneity in phosphorus availability through these electronic modification mechanisms. Near root surfaces where citrate concentrations are highest, the electron-depleted state of calcium phosphate surfaces would inhibit crystallization, potentially increasing phosphorus bioavailability.
This journal is © The Royal Society of Chemistry 2025 |