Nutrient-Doped Hydroxyapatite: Structure, Synthesis and Properties

: Complex inorganic powders based on calcium phosphates have found a plethora of practical applications. Of particular interest are the CaO-P 2 O 5 system-based multi-component material powders and granules as the source of major-and micronutrients for the plants. The emerging strategy is to use nano fertilizers based on hydroxyapatite (HAP) for phosphorus and other nutrient delivery. The doping of micronutrients into HAP structure presents an interesting challenge in obtaining speciﬁc phase compositions of these calcium phosphates. Various techniques, including mechanochemical synthesis, have been employed to fabricate doped HAP. Mechanochemical synthesis is of particular interest in this review since it presents a relatively simple, scalable, and cost-effective method of calcium phosphate powder processing. The method involves the use of mechanical force to promote chemical reactions and create nanometric powders. This technique has been successfully applied to produce HAP nanoparticles alone, and HAP doped with other elements, such as zinc and magnesium. Nanofertilizers developed through mechanochemical synthesis can offer several advantages over conventional fertilizers. Their nanoscale size allows for rapid absorption and controlled release of nutrients, which leads to improved nutrient uptake efﬁciency by plants. Furthermore, the tailored properties of HAP-based nano fertilizers, such as controlled porosity and degradation levels, contribute to their effectiveness in providing plant nutrition.


Introduction
Synthesis of multicomponent inorganic powder based on CaO and P 2 O 5 systems is critical in agriculture.Phosphorus (P) fertilizers commonly used in agriculture suffer from low efficiency due to P immobilization in the form of aluminum (Al) and iron (Fe) oxides formed in the soil immediately after P release [1].Additionally, a substantial amount of P that is released from organic matter in the soil rapidly undergoes fixation and becomes trapped in insoluble inorganic compounds [2].This process significantly limits the availability of P to crops and results in deficiencies even in soils with high overall P content [3].In other words, the released P is unable to remain in a bioavailable form, which hinders its uptake by plants and impedes their growth and development.This highlights the importance of developing strategies that can effectively address P cycling and ensure an adequate supply of accessible P to support optimal crop growth.To address this challenge, there is a need to develop fertilizers that release P at a slower rate while increasing the availability of bioactive P in a controlled manner.One potential solution being explored is the use of nanoporous hydroxyapatite (HAP) as a fertilizer [4].Previous studies have suggested that HAP nanoparticles may release P more slowly than conventional fertilizers while still mobilizing P faster than the forms immobilized in the soil [5][6][7].
The utilization of HAP as a fertilizer holds significant promise for enhancing nutrient availability and improving agricultural productivity [8,9].By providing a slow-release mechanism for P, HAP can mitigate the issue of P immobilization and increase bioavailable P levels for crop uptake.Furthermore, the ability of HAP to serve as a carrier for various macronutrients and micronutrients opens up avenues for targeted and efficient multi-nutrient delivery to plants.Water soluble mono ammonium phosphate, diammonium phosphate, and triple superphosphate are widely utilized as P fertilizers in agriculture [10][11][12].These fertilizers are applied to the soil releasing water-soluble forms of P that are readily accessible to crops [13].Their use results in substantial P losses due to leaching and surface run-off, which present challenges to food security and environmental sustainability [14].The use of less soluble P forms, such as phosphate rocks and apatite, can mitigate these losses and also show P availability for plants [15,16].Exploring the potential of nanoscale inorganic powdered materials, such as HAP, offers a promising avenue, albeit necessitating thorough and systematic investigations of the complex multi-component system synthesis techniques and the resulting fundamental properties.
Previous studies have examined the use of HAP nanoparticles as a phosphorus (P) fertilizer, capitalizing on the distinctive characteristics of nanosized particles to enhance their mobility in soil and release P [1,17].The work of Liu and Lal in 2014 introduced the concept of utilizing nanosized solid P as a fertilizer while considering the trade-off between agricultural benefits and environmental risks.They designed a greenhouse trial focusing on soybean cultivation and demonstrated that the application of HAP as a phosphorus source significantly enhanced root and aerial biomass along with increased grain yield compared to control and conventional P fertilizers [18].Another study by Marchiol et al. explored the potential of HAP not only as a P supplier but also as a carrier for other elements or molecules [19].Their germination trial with tomato plants showcased the versatility of HAP in fulfilling these dual roles.In light of the growing interest in nano HAP as a promising fertilizer component, there is a pressing need within the industry for a straightforward and efficient method to produce high-quality HAP.
The conventional synthesis processes of HAP typically involve complex chemical reactions and energy-intensive procedures, which can be both time-consuming and costly.Moreover, the large-scale production must be aligned with the nutritional requirements of plants (essential micro and macronutrients present in the synthesized nanocomposite).In this regard, the mechanochemical approach can be a potential solution to address these scale challenges.Pohshna et al. examined the use of urea/HAP nanomaterials in rice cultivation and compared their performance with commercial bulk urea fertilizer in addition to a control treatment.They reported that urea/HAP nano fertilizer showed improved grain yield, nutrient use efficiency, and reduced nutrient leaching.Therefore, urea/HAP nano fertilizer has the potential to serve as a sustainable nitrogen and phosphorus supplement for rice crops by enhancing productivity [20].Pradhan et al. synthesized nanoporous HAP nanoparticles using a sol-gel technique and loaded them with the commonly used nitrogen fertilizer urea to create a urea/HAP composite.Their results showed sustained release of urea from urea/HAP in simulated environmental conditions.Surface migration, water retention, and leachate analysis demonstrated that both HAP and urea/HAP exhibited superior water retention properties compared to conventional fertilizers.They also exhibited stability and increased nitrogen use efficiency under high temperatures.Furthermore, urea/HAP exhibited promising results in promoting rice seed germination, which was attributed to enhanced α-amylase activity and starch content.However, that indicates the ability to fabricate engineered urea/HAP promise for enhancing crop growth and productivity while minimizing environmental impact [21].Elsayed et al. introduced HAP nanoparticles as a novel nano-fertilizer for rosemary plants.HAP nanoparticles, with an average size of approximately 24 nm and a mixed morphology of rods and spherical nanoparticles, were applied via foliar application at various concentrations.HAP showed significant results compared to traditional nitrogen-phosphorus-potassium fertilizers.The treatment by HAP nanoparticles exhibited significant improvements in growth parameters, essential oil yield, and physiological indices of the rosemary plants.The stem and leaf characteristics, such as diameter, the thickness of different tissues, and pith diameter, were developed and enhanced.Moreover, the quality of the essential oil obtained from the rosemary plant was notably higher with the use of the proposed nano fertilizer [5].
Mechanochemical synthesis refers to the production of materials through mechanical force-induced reactions, typically by milling or grinding solid reactants together [22].This approach offers several advantages over traditional chemical synthesis routes.It can fulfill the need for low temperatures and decrease the usage of many chemical reagents [23][24][25], resulting in a more environmentally friendly process with reduced energy consumption.In addition, mechanochemical methods enable to be controlled regarding the composition, size, and morphology of the resulting particles [25], hence, leading to improved product quality and performance.Further, the simplicity and scalability of mechanochemical techniques make them attractive for industrial-scale production [26,27].Moreover, one of the key benefits of mechanochemical production for HAP is the ability to achieve nanoscale particle sizes and consequently enhance the reactivity [28].Furthermore, the mechanical forces applied during the milling process promote the formation of finely dispersed HAP nanoparticles, which possess a larger surface area and higher chemical reactivity compared to their bulk counterparts.This increased reactivity is particularly advantageous for fertilizers as it facilitates the efficient release of phosphorus and enhances nutrient uptake by plants [29,30].
There are several studies have demonstrated the successful application of mechanochemical methods for HAP synthesis [31][32][33].For instance, high-energy ball milling has been employed to produce nanostructured HAP powders with tailored properties.By controlling the milling parameters, such as milling time, speed, and reactants, researchers have achieved HAP particles with desired characteristics, including size, crystallinity, and surface area [34][35][36][37][38].These finely tuned properties can positively influence the fertilizer's effectiveness and performance in agricultural applications.In this regard, there is a notable increment in the number of publications about using HAP as a fertilizer.It can be seen in Figure 1 that there has been a significant increase over the years, from one article in the year 2000, to rise gradually and reach around 50-51 articles by the year 2022.In addition, there has been a growing number of articles investigating nanohydroxyapatite as a fertilizer.Over the past decade, the number of scientific Js has seen a significant increase.However, it is noteworthy that research investigating the use of nano HAP as a fertilizer did not witness a substantial surge until after 2020.The limited attention to nano HAP as a fertilizer in the years leading up to 2020 might be attributed to several factors.One possibility is that researchers were primarily focused on exploring other aspects of nano HAP applications in fields, such as biomedicine.Consequently, the potential agricultural and agronomic benefits of nano HAP as a fertilizer may have been overshadowed by other promising applications during that period.However, the subsequent upturn in research interest post-2020 could indicate a growing recognition of the untapped potential and benefits of incorporating nano HAP in agricultural practices.
In summary, the mechanochemical synthesis of HAP can offer opportunities for incorporating additional elements into the nanostructure by doping or additional compounds as nanocomposites.This approach allows for the development of tailored HAP-based materials with enhanced nutrient-release properties or the ability to deliver specific additives or growth-promoting agents to plants.In summary, the mechanochemical synthesis of HAP can offer opportunities for incorporating additional elements into the nanostructure by doping or additional compounds as nanocomposites.This approach allows for the development of tailored HAPbased materials with enhanced nutrient-release properties or the ability to deliver specific additives or growth-promoting agents to plants.

Overview
Many calcium phosphates have been studied extensively, such as HAP, with the chemical formula Ca10(PO4)6(OH)2 and a Ca/P ratio of 1.667, as well as α-and β-tricalcium phosphate (TCP) with the chemical formula Ca3(PO4)2 and a Ca/P ratio of 1.500 [39][40][41].They are widely used in various fields, such as biomedical applications, specifically in the areas of dentistry, orthopedics, and bone tissue engineering [41,42].These materials serve various purposes, including their use as dental fillers due to their mechanical properties, coatings for titanium dental implants for their biocompatibility, fillers in biodegradable composites, and substitutes for bone reconstruction and regeneration [43][44][45].The major reason for their extensive usage in these applications is their chemical resemblance to the mineral component, which is naturally found in hard tissues [46,47].
HAP exhibits two distinct crystal structures: Hexagonal and monoclinic, which are represented in Figure 2a,b.The hexagonal form has lattice parameters of a = b = 9.432 Å, c = 6.881Å, and γ = 120°, while the monoclinic form has lattice parameters of a = 9.421 Å, b = 2a, c = 6.881Å, and γ = 120° [48].Therefore, it seems that one monoclinic unit cell contains two hexagonal cells.In other words, it is a pseudo-monoclinic structure.These crystal forms both adhere to a stoichiometric Ca/P ratio of 1.67 [49].It is widely recognized that the lattice structure of HAP possesses a remarkable ability to accommodate nearly half of the elements found in the periodic system [50][51][52][53][54].This unique characteristic allows for the substitution of various ions in its three main components: Ca 2+ , PO4 3− , and OH − groups [55,56].The high structural flexibility of hydroxyapatite enables the incorporation of both bivalent ions, such as Sr 2+ , Ba 2+ , Mg 2+ , and others, as well as monovalent cations, such as Na + and K + , which typically replace Ca 2+ ions [56,57].In terms of anionic substitutions, they can occur either in the phosphate or hydroxyl positions.For example, F − and Cl − ions

Overview
Many calcium phosphates have been studied extensively, such as HAP, with the chemical formula Ca 10 (PO 4 ) 6 (OH) 2 and a Ca/P ratio of 1.667, as well as αand β-tricalcium phosphate (TCP) with the chemical formula Ca 3 (PO 4 ) 2 and a Ca/P ratio of 1.500 [39][40][41].They are widely used in various fields, such as biomedical applications, specifically in the areas of dentistry, orthopedics, and bone tissue engineering [41,42].These materials serve various purposes, including their use as dental fillers due to their mechanical properties, coatings for titanium dental implants for their biocompatibility, fillers in biodegradable composites, and substitutes for bone reconstruction and regeneration [43][44][45].The major reason for their extensive usage in these applications is their chemical resemblance to the mineral component, which is naturally found in hard tissues [46,47].
HAP exhibits two distinct crystal structures: Hexagonal and monoclinic, which are represented in Figure 2a,b.The hexagonal form has lattice parameters of a = b = 9.432 Å, c = 6.881Å, and γ = 120 • , while the monoclinic form has lattice parameters of a = 9.421 Å, b = 2a, c = 6.881Å, and γ = 120 • [48].Therefore, it seems that one monoclinic unit cell contains two hexagonal cells.In other words, it is a pseudo-monoclinic structure.These crystal forms both adhere to a stoichiometric Ca/P ratio of 1.67 [49].It is widely recognized that the lattice structure of HAP possesses a remarkable ability to accommodate nearly half of the elements found in the periodic system [50][51][52][53][54].This unique characteristic allows for the substitution of various ions in its three main components: Ca 2+ , PO 4 3− , and OH − groups [55,56].The high structural flexibility of hydroxyapatite enables the incorporation of both bivalent ions, such as Sr 2+ , Ba 2+ , Mg 2+ , and others, as well as monovalent cations, such as Na + and K + , which typically replace Ca 2+ ions [56,57].In terms of anionic substitutions, they can occur either in the phosphate or hydroxyl positions.For example, F − and Cl − ions tend to occupy the OH − sites, while CO 3 2-can substitute for both PO 4 3− and OH − simultaneously or one by one [55].Furthermore, SiO 4 4− ions can replace PO 4 3− , and there are instances where hydroxyl group sites may become unoccupied [58].This inherent flexibility in the HAP lattice structure allows for the incorporation of a diverse range of ions, which provides significant opportunities for tailoring and changing its properties to meet specific application requirements.In addition, certain substitutional ions can distort the lattice by either expanding or contracting its size [59].This depends on the difference in radius between the substituted ions, further influencing the overall properties and behavior of HAP.
tend to occupy the OH − sites, while CO3 2-can substitute for both PO4 3− and OH − simultaneously or one by one [55].Furthermore, SiO4 4− ions can replace PO4 3− , and there are instances where hydroxyl group sites may become unoccupied [58].This inherent flexibility in the HAP lattice structure allows for the incorporation of a diverse range of ions, which provides significant opportunities for tailoring and changing its properties to meet specific application requirements.In addition, certain substitutional ions can distort the lattice by either expanding or contracting its size [59].This depends on the difference in radius between the substituted ions, further influencing the overall properties and behavior of HAP.

Structure
Naturally, HAP is one of the apatite family, which is calcium apatite represented by the chemical formula Ca5(PO4)3(OH) [60].However, it is commonly denoted as Ca10(PO4)6(OH)2 to indicate that the crystal unit cell consists of two molecular units.One distinguishing characteristic of HAP is the OH − ions existence that creates inner channels toward the c-axis [49,61].These channels exhibit various random orientations of OH-dipoles under different external conditions.This unique structural arrangement allows for the transfer of protons along the c-axis under specific circumstances, such as elevated temperatures or the application of an electric field with critical values [62].Further, HAP typically presents two primary crystalline structures, which are hexagonal, such as the P63/m space group, and monoclinic, such as the P21/b space group, as shown in Figure 3a,b [55].There is a third structure which is an amorphous structure with no existing symmetry or long-range order [63].However, the orientation of OH-dipoles in HAP leads to two distinct order variations.In the hexagonal disordered structure (P63/m), OH-dipoles are spatially randomly distributed throughout the crystal and encompass neighboring unit cells, as demonstrated in Figure 3.In contrast, the hexagonal ordered structure (P63) demonstrates a parallel orientation of OH-dipoles.The same concept applies to the monoclinic phase where the disordered structure (P21/b) features oppositely oriented OH-dipoles in adjacent OH-columns, while the ordered structure (P21) exhibits parallel OH-dipoles as shown in Figure 3b [62].These various HAP structures have been experimentally determined using X-ray diffraction techniques [64,65].Further, the hexagonal structure of HAP has 44 atoms per unit cell, while the monoclinic contains a double number of atoms, which is 88 atoms per unit cell [49].Both the hexagonal and monoclinic phases of HAP exhibit a periodic repetition of their crystallographic structures throughout the entire crystal.These

Structure
Naturally, HAP is one of the apatite family, which is calcium apatite represented by the chemical formula Ca 5 (PO 4 ) 3 (OH) [60].However, it is commonly denoted as Ca 10 (PO 4 ) 6 (OH) 2 to indicate that the crystal unit cell consists of two molecular units.One distinguishing characteristic of HAP is the OH − ions existence that creates inner channels toward the c-axis [49,61].These channels exhibit various random orientations of OH-dipoles under different external conditions.This unique structural arrangement allows for the transfer of protons along the c-axis under specific circumstances, such as elevated temperatures or the application of an electric field with critical values [62].Further, HAP typically presents two primary crystalline structures, which are hexagonal, such as the P6 3 /m space group, and monoclinic, such as the P2 1 /b space group, as shown in Figure 3a,b [55].There is a third structure which is an amorphous structure with no existing symmetry or long-range order [63].However, the orientation of OH-dipoles in HAP leads to two distinct order variations.In the hexagonal disordered structure (P6 3 /m), OH-dipoles are spatially randomly distributed throughout the crystal and encompass neighboring unit cells, as demonstrated in Figure 3.In contrast, the hexagonal ordered structure (P6 3 ) demonstrates a parallel orientation of OH-dipoles.The same concept applies to the monoclinic phase where the disordered structure (P2 1 /b) features oppositely oriented OH-dipoles in adjacent OH-columns, while the ordered structure (P2 1 ) exhibits parallel OH-dipoles as shown in Figure 3b [62].These various HAP structures have been experimentally determined using X-ray diffraction techniques [64,65].Further, the hexagonal structure of HAP has 44 atoms per unit cell, while the monoclinic contains a double number of atoms, which is 88 atoms per unit cell [49].Both the hexagonal and monoclinic phases of HAP exhibit a periodic repetition of their crystallographic structures throughout the entire crystal.These structures adhere to strict stoichiometry requirements.Their ideal stoichiometric ratio of Ca atoms to P atoms must be precisely 1.67 for the ideal chemical composition of the material.This strict adherence to stoichiometry is crucial for maintaining the desired properties and functionalities of HAP.It is common to observe deviations from the ideal crystallographic structures of HAP in actual mineral samples that existed in nature.These deviations mainly come from defects within the crystal lattice, such as vacancies, interstitials, and substitutions of multiple atoms.These defects can arise due to various factors, including impurities in the mineral composition, growth conditions, or external influences during crystal formation.These changes in the Ca/P ratio and the difference from the ideal structure can significantly impact the physical and chemical properties of HAP.They introduce structural irregularities and alter the arrangement of atoms within the crystal lattice.Therefore, studying and understanding these defects is crucial for exploiting the behavior and properties of HAP in practical applications.Moreover, there are two different positions for Ca 2+ ions in the HAP structure [66].( 1) Ca(1) or Ca(I) is surrounded by six oxygen atoms from four orthophosphate groups, while (2) Ca(2) or Ca(II) is surrounded by seven oxygen atoms [67].In Ca(2), six of the oxygen atoms come from four orthophosphate groups and one from hydroxyl ion, which is demonstrated in Figure 4.The crystallographic studies reported that Ca-O distance is longer in Ca(1) than Ca(2) [68].
structures adhere to strict stoichiometry requirements.Their ideal stoichiometri Ca atoms to P atoms must be precisely 1.67 for the ideal chemical composition o terial.This strict adherence to stoichiometry is crucial for maintaining the desired ties and functionalities of HAP.It is common to observe deviations from the idea lographic structures of HAP in actual mineral samples that existed in nature.Th ations mainly come from defects within the crystal lattice, such as vacancies, int and substitutions of multiple atoms.These defects can arise due to various factor ing impurities in the mineral composition, growth conditions, or external influen ing crystal formation.These changes in the Ca/P ratio and the difference from structure can significantly impact the physical and chemical properties of HAP.troduce structural irregularities and alter the arrangement of atoms within the cr tice.Therefore, studying and understanding these defects is crucial for exploitin havior and properties of HAP in practical applications.Moreover, there are two positions for Ca 2+ ions in the HAP structure [66].( 1) Ca(1) or Ca(I) is surrounde oxygen atoms from four orthophosphate groups, while (2) Ca(2) or Ca(II) is sur by seven oxygen atoms [67].In Ca(2), six of the oxygen atoms come from four or phate groups and one from hydroxyl ion, which is demonstrated in Figure 4. Th lographic studies reported that Ca-O distance is longer in Ca(1) than Ca(2) [68].HAP demonstrates limited ability to dissolve in water and remains insoluble in alkaline solutions.Nonetheless, in the presence of acids (low pH), the phosphate and hydroxide ions within the HAP structure can react with hydrogen ions, leading to the dissolution of HAP, according to (1) [69].HAP demonstrates limited ability to dissolve in water and remains insoluble in alkaline solutions.Nonetheless, in the presence of acids (low pH), the phosphate and hydroxide ions within the HAP structure can react with hydrogen ions, leading to the dissolution of HAP, according to (1) [69].
Studies showed that doping HAP with dopants, such as manganese, magnesium, and carbonate, has the effect of increasing solubility and, therefore, decreasing chemical stability [69,70].The preparation of nanosized HAP, however, might lead to higher reactivity and increased solubility in water.Moreover, previous studies revealed that HAP exhibits two distinct binding sites on its surface.(1) Ca-site and (2) P-site.The Ca-sites are composed of calcium ions and carry a positive charge, which leads to being more inclined to attract and adsorb negatively charged species, such as acidic molecules.On the other hand, the P-sites consist of phosphate groups that carry a negative charge and consequently interact preferentially and adsorb positive charges, such as basic molecules.This property of HAP allows for selective binding and interaction with molecules based on their acidic or basic nature.In addition, the selective binding phenomenon is a result of the distribution of charges and electrostatic forces in the HAP structure.Furthermore, HAP shows two different crystallographic planes with differences in their atomic arrangements.The a and b-planes predominantly expose positively charged calcium ions, while the c-planes expose negatively charged by the existence of phosphate and hydroxyl groups.This unique feature has prompted researchers to explore the potential for achieving selective adsorption behavior based on the orientation of the prepared HAP material.Therefore, the Ca/P ratio directly affects the adsorption behavior due to adjusting the ability to interact with negative charges to positive charges.The obtaining of this rare property of having both basic and acidic sites makes HAP a significantly important material that can be modified for the desired function [69].

Ionic Doping of HAP
Seventeen elements or nutrients are needed for plant growth and reproduction.While not all of the reviewed doped HAP materials were synthesized for nutrient transport and release purposes, the example shown in this section provides for the scope and potential of doped-HAP material synthesis for agricultural applications.

Carbonate Anion (CO 3 2− )
Carbonated HAP (CHAP) is particularly similar to naturally formed HAP that is found in the human body [71].CHAP exhibits more favorable formation during precipitation compared to the ideal stoichiometric HAP.The structure of CHAP is relatively complex due to its ability to be in various sites and consequently has three types: Type-A, where a CO 3 2− ion replaces two OH − ions, and Type-B, where a CO 3 2− ion substitutes PO 4 3− ion.In Type-AB CHAP, two CO 3 2− ions simultaneously replace one PO 4 3− and one OH − .These types can be determined experimentally using Fourier Transform infrared spectroscopy (FTIR) [58,[72][73][74].

Chlorine Anion (Cl − )
Fahami et al. produced ClHAP and examined the crystal structure by XRD, and reported an increase in crystallite size upon substitution by Cl -.In addition, an increase in the lattice volume was noticed [75].Naqshbandi et al. prepared ClHAP and studied the crystallinity behavior.They reported that the introduction of Cl − ions into the crystal structure of HAP led to the substitution of OH − ions.Therefore, a noticeable impact on the lattice parameters has been observed.The larger size of Cl − ions relative to OH − ions resulted in expanding the lattice parameters.Therefore, the incorporation of Cl − ions results in a slight increase in both the crystallite size and the unit cell volume [76].
Moreover, Hsu et al. prepared ClHAP/polydopamine as a coating for Ti64 alloy.They reported that chloride release increases with time in SBF [77].This indicates the ability to dope HAP with Cl − and being used as a fertilizer for providing phosphorus, calcium, and chlorine.

Nitrate Anion (NO 3 − )
Gafurov et al. mentioned that N can be substituted in HAP as NO 3 − and they can behave as CO 3 2− because they are isoelectronic ions [78].They reported that these impurities can be incorporated from the preparation by using calcium nitrate (Ca(NO 3 ) 2 ).Further, Biktagirove et al. reported that NO 3 − can locate in A-site and B-site instead of OH, and PO 4 3− respectively [79].These results open the ability to obtain HAP that might be used in N and P fertilizers at the same time without being in a composite.That means if urea is incorporated into the structure, the results might be effective in the agricultural field.

Sulfate Anion (SO 4 2− )
Alshemary et al. reported a decrease in crystal volume at small amounts of substitution, then increased upon increasing the amount of SO 4 2− due to its relatively large ionic radius.In addition, crystallite size decreased upon increasing the content of SO 4 2− [80].Another work by Alshemary et al. investigated sulfur-doped HAP (SHAP) in vitro towards osteoblast-like cells.They showed that ALP activity increased dramatically upon the existence of sulfur in the HAP structure.Furthermore, SHAP showed higher biocompatibility than unmodified HAP [81].Furthermore, Radha G et al. examined the increment of S substitution and its biological effect on MC3T3-E1 cells.They demonstrated that cell proliferation enhanced until 10SHAP, while at 15SHAP, the proliferation decreased after 1, 3, and 7 days of culture [82].

Borate Anion (BO 3 3− )
Boron is considered a micronutrient that is important for growth.In this regard, Atila et al. prepared nanocomposite of cellulose/B-HAP/gelatin, and they reported that the incorporation of B in HAP led to a decrease in the crystallinity of HAP.Further, mechanical properties were directly affected by this incorporation which led to an increase in Young's modulus from 75.26 to 94.9 MPa.Moreover, the thermal stability in the existence of B in the HAP lattice was improved [83].In addition, Tuncay et al. prepared B-HAP, and they mentioned the ability of B to enhance cell response and cell viability in vitro [84].Furthermore, Arslan et al. reported that doping by B in HAP can directly affect biocompatibility, which improved as a result of this doping [85].

Iron (Fe 2+ and Fe 3+ )
The inclusion of ferrous cation (Fe 2+ ) and ferric cation (Fe 3+ ) in hydroxyapatite affects the crystal parameters and local ion reorganization.The different electron configurations and ionic radii (70 pm for Fe 2+ and 60 pm for Fe 3+ ) contribute to their distinct interactions and positions within the crystal.A study was performed on an 88-atom monoclinic HAP supercell for incorporating ferrous and ferric cations in HAP.Makshakova et al. showed that cations can be introduced into HAP at both Ca(1) and Ca(2) positions.Further, Fe 2+ can occupy either position without requiring additional charge compensation and resulting in the formula Ca 9.5 Fe 0.5 2+ (PO 4 ) 6 (OH) 2 .On the other hand, the substitution of Ca 2+ by Fe 3+ can lead to (Ca 9.5 Fe 0.5 3+ (PO 4 ) 6 O 2 H 1.5 ), which is element perturbation.This possibly causes vacancy formation in the H + or Ca 2+ positions by releasing one Ca 2+ for every two substituted atoms (for neutralizing the charge).Further, the release of H + from OH − was also considered to maintain a neutral overall charge and to minimize the cell contraction caused by charge compensation.The presence of OH -within the anion channel, surrounded by calcium ions, can show a minimal impact on the overall packing and density of the crystal [86].In addition, they reported that Fe 3+ substitution is more reactive with the surrounding.

Magnesium (Mg 2+ )
Karunakaran et al. showed that doping HAP with Mg 2+ ions leads to less crystallinity and higher distortion in the lattice using XRD in addition to a growth in (a) crystal parameter [87].Another study by Laurencin et al. reported the lattice distortion by Mg incorporation but mentioned falls in the lattice parameter (a) due to the substitution [88].Further, Mg 2+ can substitute Ca (1) or Ca (2) sites and lead to a shortening distance between Mg and O [89].In addition, Goldberg et al. prepared Mg-HAP and reported that the substituting using Mg 2+ ions led to an increase in BET-specific surface area from 23.35 to 86.91 m 2 /g from 0 to 1% substitution.They also reported that cell volume decreased upon increasing the amount of Mg 2+ doping in addition to the adsorption capacity.Further, Mg-HAP led to an increase in the viable cells and proliferation of MG-63 cells after 24 h [68].Predoi et al. reported the antibacterial properties of Mg-HAP in addition to its biocompatibility as a biofilm [90].In another study by Predoi et al., they prepared a thin layer of Mg-HAP by dip coating procedure.It was reported that 25Mg-HAP has antibacterial activity against Staphylococcus aureus, Escherichia coli, and Candida albicans strains.In addition, the layer improved biocompatibility without affecting the cell behavior [91].Further, Sharma et al. prepared Mg-HAP/urea, and they reported the improvement of nanohybrid compared to the traditional fertilizers [92].

Zinc (Zn 2+ )
Sharma et al. doped HAP with zinc and integrated ZnHAP with urea to be used as a nano-fertilizer.The incorporation of zinc can be applied to the growth of plants as an essential element of this process.Further, nitrogen can be taken up from urea, while phosphorus and calcium essentially come from HAP.They reported that growth is highly improved compared to the control sample, which confirms the ability of doped-HAP with the micronutrients to be significantly effective [92].There are previous reports by Guerra-Lopez et al. that mentioned that Zn 2+ prefers the Ca(2) site [93].Furthermore, the same results are done by Chathurika et al.They reported that Zn-HAP/urea nano seed coating showed around 69% growth for Zea mays seeds [94].Abeywardana et al. prepared Zn-HAP/urea nanohybrids to offer a sustainable solution for addressing urea and phosphorus deficiencies in plants while providing zinc supplementation.They reported that these biodegradable nanoparticles effectively deliver essential nutrients during early plant development and resulting in improved growth, development, and nutrient distribution in maize.The plant height reached around 140 cm for the nanocomposite and 80 cm for the control sample.In addition, the number of leaves and root lengths reached around 20 and 70 cm using Zn-HAP/urea, while 14 and 30 cm for the control sample [95].

Nickel (Ni 2+ )
Kurinjinathan et al. synthesized pure and Ni-doped hydroxyapatite (HAP) nanoparticles using ultrasonication-assisted wet-chemical techniques.The synthesized nanoparticles were systematically characterized using various techniques.X-ray diffraction (XRD) analysis revealed that the concentration of Ni 2+ ions increased in the HAP structure while the crystallite sizes decreased.The average crystallite sizes for pure HAP and Ni-doped HAP with concentrations of 0.01 M, 0.05 M, and 0.1 M were found to be 29.09,16.57, 14.33, and 13.66 nm, respectively [96].Sebastimmal et al. employed the sol-gel method to synthesize Ni-doped HAP nanoparticles.SEM examination revealed that the particles exhibited a spherical or nearly spherical shape.Furthermore, Ni 2+ doped HAP nanoparticles displayed remarkable antibacterial activity against a range of bacteria, including S. flexneri, E. coli, P. aeruginosa, K. pneumoniae, and S. aureus, as demonstrated by the disc diffusion technique.Thus, the importance of Ni dopants in HAP stricture is denoted by the applicability of these nanoparticles for medical applications as antibacterial agents.Moreover, these nanoparticles showed high levels of antifungal, antioxidant, and anticancer activities, making them highly promising for biomedical applications [97].

Copper (Cu 2+ )
Bazin et al. doped HAP with Cu, and they reported the ability of Cu dopant to increase lattice parameters and volume [98].Unabia et al. prepared Cu-HAP coating, and they also mentioned the increase in lattice parameters and volume of HAP lattice upon increasing Cu content [99].Imrie et al. prepared Cu-HAP, and they reported the expansion of the unit cell due to the doping, where volume and lattice parameters increased [100].In addition, Stanic et al. synthesized Cu-HAP using the neutralization method, and they examined the released ions.They reported that HAP in phosphate buffer solution released extremely small amounts of Cu ions which are non-detectable [101].Hidalgo-Robatto et al. studied the incorporation of Cu in HAP up to 3.31% and their effect on pre-osteoblast cells.They reported that cells proliferation and had healthy growth after 7 and 21 days in addition to their ability to prevent bacterial formation [102].It should be mentioned that Cu is a plant micronutrient.

Potassium (K + )
Potassium is considered one of the essential macronutrients that are required for plant growth.It is a monovalent cation, however, it might replace Ca(II) ions partially.In this regard, to maintain neutrality, vacancies can be formed through the HAP structure, which might affect the chemical stability and the degradability of the formed crystals.Chahkandi et al. prepared K-HAP and revealed that particle size was, on average, 65-70 nm through XRD and TEM.Further, they reported the antibacterial activity of K-HAP in addition to the dye removal efficiency against Congo red (CR) [103].

Molybdenum (Mo 2+ )
The Molybdenum element is also considered one of the plant micronutrients.Abutalib et al. studied the incorporation of Mo in HAP lattice and reported that volume gradually decreased upon doping (from 0.05% to 5.00%) [104]. 2− to obtain Mg-HAP, CHAP, and Mg-CHAP.They found that the incorporation of the lattice structure by both elements can be more efficient energetically [89].On the other hand, Landi et al. synthesized SrMgHAP as a filler sample to be used in osteoporotic bone regeneration.They reported that the co-dopant substitution improved the cell response and bone metabolism, in particular at 10 mg/mL dose [105].Sharma et al. prepared ZnMg-HAP/urea nanohybrids for improving wheat growth and uptake of nitrogen.They found that doping HAP nanoparticles enables the delivery of both macro and micronutrients to plants and makes them a biocompatible system suitable for agricultural use.Coating these nanoparticles with half and quarter doses of urea results in improved crop yield, nutrient levels, and leaf nitrogen content while simultaneously reducing nitrogen loss.The controlled and gradual release of nitrogen from the urea-coated ZnMg-HAP nanoparticles ensures an adequate and sustained supply of nitrogen during the crucial early growth stage of wheat plants.Furthermore, these nanohybrids effectively prevent the leaching or volatilization of urea from the soil, thereby ensuring a consistent and long-lasting nutrient delivery system [92].However, the application of these nanohybrids showed no detrimental effects on wheat growth, which offers a promising alternative to conventional fertilizer methods in agriculture.Tang et al. prepared NiMo-HAP, and they reported that the incorporation of Ni and Mo showed a negative effect on the stability as a catalyst [106].Moreover, Alshemary et al. demonstrated the incorporation of Sr 2+ ions in SHAP, and they reported a high increase in biocompatibility and cytocompatibility of Sr-SHAP [81].

Synthesis of HAP
The synthesis of HAP involves several methods with distinct advantages and yielding HAP particles with specific properties.There are two main ways for obtaining HAP, which are (a) wet chemical methods, as described in Figure 5, and (b) solid-state reactions shown in Figures 6 and 7.
in agriculture.Tang et al. prepared NiMo-HAP, and they reported that the incorporation of Ni and Mo showed a negative effect on the stability as a catalyst [106].Moreover, Alshemary et al. demonstrated the incorporation of Sr 2+ ions in SHAP, and they reported a high increase in biocompatibility and cytocompatibility of Sr-SHAP [81].

Synthesis of HAP
The synthesis of HAP involves several methods with distinct advantages and yielding HAP particles with specific properties.There are two main ways for obtaining HAP, which are (a) wet chemical methods, as described in Figure 5, and (b) solid-state reactions shown in Figures 6 and 7.

Coprecipitation
The wet chemical synthesis allows producing HAP nanoparticles with various ranges in scale, different morphologies, and crystallinity.One of the most simple and basic wet chemical synthesis techniques is the coprecipitation method which is frequently used [67,107].The coprecipitation method involves the simultaneous precipitation of calcium and phosphate ions from a solution to form HAP nanoparticles.This is done by mixing calcium and phosphate precursors under controlled conditions, such as temperature, pH, and stirring rate [108].This leads to the nucleation and growth of HAP particles.The  The wet chemical synthesis allows producing HAP nanoparticles with various ranges in scale, different morphologies, and crystallinity.One of the most simple and basic wet chemical synthesis techniques is the coprecipitation method which is frequently used [67,107].The coprecipitation method involves the simultaneous precipitation of calcium and phosphate ions from a solution to form HAP nanoparticles.This is done by mixing calcium and phosphate precursors under controlled conditions, such as temperature, pH, and stirring rate [108].This leads to the nucleation and growth of HAP particles.The coprecipitation process can be performed at ambient temperature or at elevated temperatures which depends on the desired properties of the nanoparticles.One of the key advantages of the coprecipitation method is its simplicity and cost-effectiveness, as presented in Table 1.[109].This method has no complex requirements, such as complex equipment or specialized devices.Therefore, it is suitable for researchers to obtain samples and industries for large-scale production [110][111][112].In addition, the materials needed for the synthesis are readily available and relatively inexpensive, which makes it more cost-effective.Another advantage of the coprecipitation method is its ability to provide high yields of HAP nanoparticles within a relatively short time.Furthermore, the reaction kinetics can be easily controlled by adjusting the temperature, pH, and concentrations in addition to the started precursors [113].This availability of controlling various conditions can directly affect the properties of obtained HAP powder.By adjusting the concentrations of the calcium and phosphate precursors, the stoichiometry of the HAP can be finely tuned.Furthermore, controlling the reaction parameters, such as temperature (almost low temperature in this method) and pH, enables the modulation of particle size.This flexibility in composition and particle size control is crucial for tailoring the properties of HAP nanoparticles to suit specific applications.Another advantage from the point of view of ceramics production is the suboptimal crystallinity of the precipitated HAP nanoparticles in addition to the rapid precipitation process, often resulting in the formation of nanocrystalline and amorphous HAP.Therefore, lower crystallinity compared to other synthesis methods might affect the desired application [71].Needs special devices, nanoparticles tend to aggregate, contamination possibility from milling media, a wide range of particle sizes and shapes (inhomogeneous particulate) [114,116] Despite these advantages, the coprecipitation method has limitations and disadvantages.One of the main challenges is achieving a uniform particle size distribution and specific crystal morphology [117,118].The coprecipitated nanoparticles tend to agglomerate, which leads to poor dispersion and affects their performance in certain applications.This issue can be addressed through additional processing steps, such as sonication or surface modification techniques [119].These additional steps might be considered disadvantages.Moreover, the coprecipitation method can be sensitive to reaction conditions such as temperature, pH, and reaction time, which need to be carefully controlled to ensure the desired properties of the HAP nanoparticles.Any deviation in these parameters can lead to variations in particle size, composition, and crystallinity.Therefore, it affects the overall quality of the synthesized nanoparticles.Sciena et al. prepared HAP by the coprecipitation method and reported the improvement of P ion release [120].

Hydrothermal Method
The hydrothermal method is a widely used technique for the synthesis of HAP nanoparticles with controlled size, morphology, and crystallinity [121,122].This method involves the reaction of calcium and phosphate sources in a closed system under high pressure and temperature [122].The hydrothermal synthesis process promotes the transformation of amorphous calcium phosphate precursors into well-defined crystalline HAP nanoparticles [123].One of the main advantages of the hydrothermal method is the ability to produce highly crystalline HAP nanoparticles with a narrow size distribution [123,124].Well-defined crystal structures and high crystallinity can be obtained due to exposing a closed system to high-pressure and high-temperature conditions.In addition, controlling reaction time, temperature, and precursor concentrations can lead to tuning the particle size and crystallinity precisely.Further, the system under high temperature and pressure can directly affect the purity leaving the pure desired shape of HAP nanoparticles.The precise control over the crystalline nature of the nanoparticles contributes to their superior mechanical strength, stability, and biocompatibility, making them suitable for various applications.In addition, various morphologies can be obtained by adjusting the reaction conditions [125].On the other hand, hydrothermal methods need specialized equipment to stand in these high conditions.The usage of high-pressure vessels and autoclaves increases the complexity and cost of the synthesis process, which directly limits its accessibility.Moreover, the reaction time, which ranges from several hours to several days, is another drawback of this method [126].Moreover, the high-temperature and high-pressure conditions can promote particle agglomeration and lead to the formation of larger aggregates rather than individual nanoparticles.However, additional processing steps, such as son-ication or surface modification techniques, may be employed to disperse and stabilize the HAP nanoparticles.In this regard, Tang et al. prepared HAP to be used in fertilizer application, and the preparation included a hydrothermal process.They obtained nanorod particles, which are efficient as P-delivery where phosphorus use efficiency reached around 47% compared to the chemical phosphorus fertilizer, which reached around 23% [7].Solid-state reactions are employed for HAP synthesis in addition to the wet methods.These methods involve the direct heating or sintering of calcium and phosphate precursors at high temperatures.Usually, solid-state reactions involve high temperatures, which is a major requirement to activate the solid reactants for chemical interaction or crossing the activation barrier [127].One widely used technique in nanoparticle preparation, in general, is the solid-state thermal decomposition method [128][129][130].The precursors are heated to temperatures above 800 • C, and hence nucleation and growth occur to form crystalline HAP.The thermal decomposition method is relatively simple and cost-effective for small production but results in larger particle sizes and lower crystallinity compared to wet chemical synthesis methods.

Mechanochemical Synthesis
The mechanochemical method or mechanical activation is based on merging mechanical energy with solid-state reaction [27].The conventional process of solid-state reactions consumes significant time and energy, but it can be replaced by incorporating mechanical energy [131].Introducing mechanical energy into the process can achieve various benefits: (1) Increased reactivity, (2) overcoming the activation barrier, (3) improved diffusion, and (4) reduced reaction time.The technique can approach its purpose through grinding or milling by inducing mechanical stress to the solid reactants.This stress can disturb the surface-bound compounds leading to enhance the reaction.Further, grinding or milling leads directly to breaking down the particles into smaller ones.In more detail, when sample powders are crushed between colliding balls inside a ball mill pot, they undergo plastic deformation.That can lead to cold welding, fracturing, and rewelding.The impact force breaks crystallographic bonds and creates new surfaces on the particles.These new surfaces facilitate easy welding of particles together, resulting in an accelerated dissolution rate of the solid material [131].Continued mechanical deformation further reduces particle size and increases surface energy, causing various changes in the properties such as surface-tovolume ratio and, consequently, reactivity [131].In addition, chemical reactions typically require a certain amount of energy, which is commonly achieved by raising the temperature.However, the introduction of mechanical energy can contribute to the required energy and reduce the reliance on temperature and allow for lower temperatures to overcome the activation barrier.Furthermore, milling and grinding apply intense forces causing deformations, defects, and fractures.The decrease in particle size leads to more contact points between reactants which increases diffusivity and reactivity, then directly improving the reaction rate.The mechanochemical approach is a technically straightforward method used in the large-scale production of nanocrystalline HAP powder.When compared to the wet approach, it offers advantages in the aspect of simplicity and versatility [132].The mechanochemical process does not require parameters, such as reactant concentrations, temperature, aqueous solution, or long processing times.That means the final morphology and size can be affected by other factors.The mechanochemical process requires different parameters that can be easily controlled such as speed or frequency, duration of milling, and milling media.To apply mechanical energy during the synthesis of HAP, commonly used equipment includes ball-milling or planetary mills.The device offers the option of adjusting the speed and duration of milling.Then, by utilizing high-energy milling, the interaction between the calcium source and the phosphorus source is stimulated.This leads to improved reaction rates.Further, this milling process enhances diffusion processes and allows for the synthesis of nanocrystalline HAP powder at lower temperatures and shorter reaction times.Therefore, the high-energy milling approach is highly advantageous for large-scale production of HAP powder with desirable properties.Moreover, there is no need for external temperature while the process is working.The collision and friction between the balls can generate heat (from kinetic energy), and this heat can facilitate the reaction.Various factors, including the reactant types, milling medium, milling ball type, and diameter, duration of milling steps and interval pauses, powder-to-ball mass ratio, and rotational speed, are crucial processing variables in the mechanochemical method [133].In addition, it is reasonable to assume that the nanoscale characteristics of HAP obtained mechanochemically could influence the release of P ions in the soil.This is because nano scaling can increase the surface area and reactivity of HAP particles.That potentially affects the dissolution behavior and subsequent release of P ions.Studies have shown that nanosized materials can exhibit enhanced reactivity and solubility compared to their bulk counterparts, which may lead to altered release patterns of nutrient ions, such as P, in soil [134].Moreover, a study conducted by Silva et al. prepared urea-loaded HAP using the mechanochemical method.They investigated the release of N from the urea-loaded HAP and found that the release kinetics followed the Higuchi model indicating a diffusioncontrolled mechanism.However, it should be noted that in this study, they did not measure the release of P ions [135].

Economic Feasibility of the Mechanochemical HAP Synthesis Method
Han et al. highlighted the advantages of using mechanochemistry in sustainable large-scale processing as a green, solvent-free, scalable, and low-temperature process.In their study, they achieved significantly higher ammonia concentration compared to the traditional Haber-Bosch process.While the traditional method reached a concentration of only 25 vol%, the mechanochemical approach yielded an impressive 82.5 vol% ammonia concentration.The milling conditions for the mechanochemical method involved a temperature of around 45 • C and 1 bar pressure, whereas the Haber-Bosch process required much higher conditions with temperatures of around 450 • C and pressures of 200 bar.Despite the higher concentration achieved with mechanochemistry, the total energy consumption was remarkably lower, suggesting a potential for many practical applications.They evaluated the total energy consumption for the mechanochemical method at 4.5 × 10 12 J•m/t, while the thermal Haber-Bosch and electrochemical methods consumed 39 × 10 12 J•m/t and 228 × 10 12 J•m/t, respectively.It must be mentioned that these energy consumption values have relatively small yields compared to the mechanochemical method [136].AlShamaileh et al. emphasized the need to reduce fertilizer consumption, particularly nitrogen (N) fertilizers, to mitigate environmental pollution.They pointed to the mechanochemical method as a potential solution for reducing consumption and its associated environmental impact [137].Furthermore, Mohamed et al. discussed the benefits of the mechanochemical method for synthesizing cathode material containing lithium (Li), and iron (Fe), with a chemical formula of (Li 2 Fe)SO, in one step, leading to cost reduction [138].In another study conducted by Palaz et al., they highlighted the advantages of milling as a continuous process, offering high productive capacity and low energy consumption.They compared this method to traditional approaches for producing phosphate fertilizer, which typically require 450 kWh/t in a periodic process lasting for 20 days and involve the use of reagents such as sulfuric acid (H 2 SO 4 ).In contrast, the mechanochemical process was found to require 500 kWh/t, but only for a significantly shorter duration of 24 h, which aligns with the production process [139].The remarkable difference lies in the elimination of the need for reagents in the mechanochemical approach, making it a more efficient and eco-friendly alternative.Moreover, its continuous nature ensures a more streamlined and continuous production, enhancing overall productivity [139].Finally, these studies collectively demonstrate the potential of mechanochemistry as a sustainable and efficient approach for various applications, including fertilizer, such as HAP, production, and all contributing to reducing energy consumption and environmental impact.Silva et al. conducted a study focused on HAP synthesis using ball milling and investigated the effect of milling time.They found that crystallite size increased from 17.5 nm after 5 h of milling to 20.5 nm after 15 h [140].This suggests that prolonged milling leads to a gradual increase in crystallite size.Nasiri-Tabrizi et al. focused on fluorapatite (FAP) (Ca 10 (PO 4 ) 6 F 2 ) and observed changes in crystal size with milling time.They found that the crystal size of FAP increased to approximately 90 nm after 50 min of milling, reached around 55 nm after 100 min, and eventually reduced to 35 nm after 300 min of milling [141].The milling process also induced lattice strain, which increased from 0.2% to approximately 0.5% over 300 min.They noted that the milling process gradually eliminated presence of CaO, with complete disappearance observed after 300 min of milling.Further, the morphology of the particles changed from flake-like shapes to more rounded particles as the milling time increased.These observations indicate that the milling duration and conditions can have a significant impact on the final crystal size of the synthesized HAP or its derivatives.Ebrahimi-Kahrizsangi et al. focused on the synthesis of spheroidal nanoparticles of FAP.They achieved an average crystalline size of 25 nm and an average lattice strain of 0.434%.The milling process included the use of zirconia balls with a diameter of 20 mm, a low rotational speed of 600 rpm, and a high ball-to-powder weight ratio (BPR) of 20:1.To prevent excessive temperature rise during milling, they conducted the process in steps, with 45 min of milling followed by a 15-min rest period for a total duration of 60 h.[142].Fakharzadeh et al. reported the results of a study on HAP synthesis through ball milling and investigated the effects of milling time.They achieved a crystallite size of approximately 30 nm and a lattice strain of 0.5% after 3 h of milling.The milling conditions included a BPR of 15:1, a rotational speed of 600 rpm, and the use of steel balls with a diameter of 20 mm.The resulting HAP particles exhibited a spherical morphology and showed a tendency to agglomerate [143].Further, it can be seen that at low rotational speed, more time is needed to obtain small crystallites.Bulina et al. investigated the synthesis of HAP particles through ball milling.They reported that after 30 min of milling at a high rotational speed of 1800 rpm, they obtained HAP particles with a non-uniform size distribution and a crystallite size of approximately 24.2 nm [35].This study demonstrates the influence of milling parameters, such as rotational speed, on the particle size distribution and crystallite size of the synthesized HAP.It can be seen that 1800 rpm is relatively high, which indicates less time to obtain nanocrystals.Therefore, the use of high rotational speed allows for the efficient production of nanosized HAP particles within a short milling time.Chaikina et al. revealed a study focusing on the impact of milling time on the crystallite size of HAP.They demonstrated that longer milling times increased the crystallite size, ranging from approximately 20 to 32.7 nm.They conducted milling for durations of 3, 5, 10, 20, and 30 min, resulting in crystallite sizes of approximately 20, 20.6, 28.9, and 32.7 nm, respectively.The milling conditions included a high rotational speed of 1200 rpm and a high BPR of 20:1 [144].These values shed light on the relationship between milling time, BPR, and rotational speed with crystallite size.Therefore, high speed and high BPR lead to reaching nanosized particles in a short time, and when milling time increases, the crystal size increases.In another study by Bulina et al., they observed spherical HAP particles with sizes less than 50 nm after 30 min of milling at 1200 rpm, as revealed by TEM [145].The crystallite size determined by XRD was approximately 20 nm.The influence of milling parameters on particle morphology can be noticed.This study demonstrated the efficient synthesis of nanosized HAP particles within a short milling duration when the high rotational speed is incorporated.
Abdel-Aal et al. investigated the effect of milling time on the crystallite size of HAP.They found that the crystallite size decreased from around 13 to 6 nm after 3 and 5 h of milling, respectively.The experimental setup involved the use of a fine particle grinder instrument and a rotational speed of 2000 rpm [146].This study demonstrated the impact of milling duration on the crystallite size of HAP, which showed that longer milling times resulted in a reduction in crystallite size.Moreover, it can be observed that the absence of grinding balls and being replaced with another instrument affect the ability to grind directly.Iwasaki et al. employed intermittent milling to avoid temperature rise during the process.The net wet milling time was 1 h, and they successfully synthesized HAP/La 0.7 Sr 0.3 MnO composite.The HAP particles showed sizes ranging from 50 to 100 nm.The milling conditions included a rotational speed of 600 rpm, zirconia balls with a diameter of 5 mm, and a BPR of 18:1 [147].This study showcases the use of intermittent milling to control the synthesis of composite materials and achieve specific particle sizes.In addition, it can also be noticed that low rotational speed, low diameter, and high BPR increase the needed time to approach the nanoscale.
Galotta et al. investigated the synthesis of HAP particles through ball milling and showed that after 4 h of milling, then followed by drying at 150 • C, HAP particles had a crystallite size of approximately 21 nm.The particles exhibited a flaky morphology with nanoscale thickness.The experimental setup involved the use of zirconia balls with diameters of 6 mm, a rotational speed of 50 Hz (3000 rpm), and a BPR of 5:1 [37].This reveals the influence of milling parameters on the morphology and size of the synthesized HAP particles.They used low BPR, low diameter, and high rotational speed.Fahami et al. focused on the synthesis of HAP/Ti nanocomposites through ball milling, and they observed a decrease in the crystallite size from 20 to 13 nm after 20 h of milling, which was accompanied by an increase in lattice strain.Morphologically flaky and ellipsoidal particles were observed, indicating a notable reduction in size after 10 h of milling.In addition, the experimental conditions included zirconia balls with a diameter of 20 mm, a rotational speed of 600 rpm, and a BPR of 20:1 [148].It can be seen that the incorporation of TI led to an increase in the needed time for reducing the crystallite size, especially at a low rotational speed.Ferro et al. studied the synthesis of HAP from chicken eggshell and cuttlefish bone.The used balls are made of alumina with a diameter of 10 mm and BPR of 10:1 in addition to 400 and 600 rpm.They found that HAP structure was obtained after prolonged milling for approximately 10-12 h using high milling energy.XRD analysis revealed broad HAP peaks, indicating the nanosized of the HAP crystals [149].That indicates the precursors used to synthesize HAP have a potential effect on the milling time, especially at the low rotational speed.Coreno A et al. conducted a study on HAP synthesis by milling the reactant mixture every 2 h.All reactants transformed into carbon-substituted HAP after 6 h of milling with particle sizes ranging from 10-20 nm.They used hardened steel balls with a BPR of 20:1 throughout the study [150].Lala et al. investigated the doping effect of Zn, Mn, and Mg on HAP, which are prepared by ball milling.The preparation was done using a high planetary ball mill with a hardened chrome steel vial for 10 h in addition to pausing 10 min every 20 min.They reported that sintering the samples at 950 • C led to decomposing the undoped Mn, and Mg-HAP partially to β-TCP, while Zn-HAP decomposed partially to the tetracalcium phosphate phase [151].
Bulina et al. prepared magnesium-substituted HAP during a 20-min milling duration at a high rotational speed of 1200 rpm.They reported a crystallite size of approximately 18 nm, and the SEM images showed particles in spherical and rod shapes with sizes reaching 20 nm [152].Youness et al. systematically investigated the effect of milling time on B-CHAP, varying the milling time from 1 to 10 h with 1-h increments and including a 30-min pause every 60 min.The milling process was at a high rotational speed of 1500 rpm using alumina balls with a diameter of 10 mm and a BPR of 10:1.They observed an increase in crystal size from 13 to 30 nm, while the lattice strain decreased from 0.9506% to 0.4420% [153].The increase in the crystal size is due to using high BPR, high diameter, high speed, and long milling time.Hussain et al. studied preparing HAP from waste seashells by hydrothermal and ball milling under conditions of 40 min with a 10-min break every 15 min at 350 rpm and BPR 10:1.They showed that heat treatment at 900 • C led to the maximum crystallinity with near stoichiometric ratio with rod, spherical and needle shapes in addition to crystallite size around 36 nm [154].
Ravindranadh et al. synthesized calcium-lithium HAP (CLHAP) nanopowder using mechanochemical milling with a low rotational speed and 3 h of milling.They included pauses in the milling process, with a 15-min rest period every 60 min.The obtained CLHAP powder exhibited a crystallite size of 16 nm and a morphology resembling stone-like particles [155].Furthermore, Nasiri-Tabrizi et al. prepared HAP under various mechanochemical conditions such as milling time, milling atmosphere, milling media, and precursor compositions.They reported that they fixed BPR 20:1 of zirconia balls of 20 mm diameter and rotational speed of 600 rpm.They studied the crystallite size concerning milling time and they found that crystallite size decreased from around 40 to 30 nm after 40 and 80 h, respectively.These results are estimated when they started the reaction using CaHPO 4 +Ca(OH) 2 with air and argon atmosphere in addition to polyamide-6 (PA6) medium.There was no difference in the results estimated by changing the milling atmosphere.Moreover, they changed the precursor compounds to CaCO 3 +CaHPO 4 , and the results changed.The crystallite size of 40 h showed a smaller size in the air with a value of 32 nm.Furthermore, the argon medium showed 40 nm at the same duration of milling, and by increasing the milling time, the crystallite size decreased to 28 and 32 nm for the air and argon media, respectively.These results indicate the complex effect between precursors, milling media, and size.Furthermore, under the fixed previous conditions of speed, BPR, with 80 h of milling, changed atmosphere, precursors, and milling media.They found that in the polymeric media (PA6) starting with CaHPO 4 +Ca(OH) 2 in air, the morphology of nanorods, cluster-like structures, and agglomerations.In the argon medium, a coalescence of nanorods is observed, which are chain-like agglomerations.In the metallic medium (only argon atmosphere), the only observed are coalescence agglomeration nanoparticles.By starting with CaCO 3 +CaHPO 4 , nanospheres are observed that make obvious chain-like clusters, especially in the argon atmosphere of polymeric medium [156].This study showed that not only the milling factors are important to study, but also precursors, vial media, and even atmosphere are essentials.
A consistent trend can be observed in the previous studies investigating HAP synthesis using the mechanochemical method.It is evident that as the milling duration increases, there is a general decrease in the crystallite size of the material until a critical point is reached.Beyond this critical point, recrystallization occurs and results in an increase in the crystal size accompanied by a decrease in lattice strain.The duration of milling plays a significant role in this trend, but other factors related to the milling process also contribute.All the factors are translated into the material as energy which transfers into the sample.
Burgio et al. developed an equation to express the energy transferred into the sample from the planetary mechanical milling according to (2) [149,157]: The parameters involved in the (1) are as follows: ϕ is a numerical value dependent on vial filling, n b is the number of balls, m b refers to the mass of each ball, t is the milling time, w p refers to the angular rate of the mill, w v is the angular rate of the vial, r v indicates the vial radius, d b is the ball diameter, and r p indicates the distance between the center of the plate and the center of the vial.This means all these parameters can affect the transferred energy into the sample.However, there are main parameters that can be mainly focused and controlled, such as rotational speed, time, and the ratio between the weight of balls and powder.For instance, Chaikina et al. reported an unexpected increase in crystallite size after only 3 min of milling.This can be attributed to the high rotational speed of 1200 rpm and a BPR value of 20:1 employed in their study.The combination of these factors led to the rapid and efficient breakage of grains, facilitating the cold welding and recrystallization process even within a short milling time [144].Further, in the morphological aspect, the up-bottom technique of ball-milling breaks the grains, which tend to be in spherical and ellipsoidal shapes, unlike coprecipitation in which HAP particles tend to be in nanorod shape.

Effect of Ball Milling on Biological Properties
Angioni et al. investigated the biological response of HAP/bioglass named BGMS10 with 30 min of milling.They reported that this milling time was optimized with 800 • C of heat treatment at 70 MPa of pressure which led to covering SBF completely with HAP after 7 days of immersion.On the other hand, 120 min of milling led to lower bioactivity and lower mechanical properties [158].Sharifianjazi et al. showed that PHAP showed around 60% cell viability in 7 days and around 100% ALP activities [159].Lala et al. prepared CHAP by mechanochemical method at high energy for various milling times.They found that the best biocompatibility near the natural bone was in the sample after 2 h of milling.The high biocompatibility was estimated even at relatively high concentrations of 100 µg/mL, which was around 85% of cells were alive after 24 h of incubation.Furthermore, after 10 h of milling, the sample showed around 100% cell viability from the low concentrations to the relatively moderate of 60 µg/mL after 6, 12, and 24 h of incubation [160].Fernandez et al. prepared HAP and β-TCP through chemical methods and combined them with bioactive glass containing 45SiO 2 -24.5CaO-24.5Na 2 O-6P 2 O 5 .Further, the biphasic bioceramic was treated mechanically, and they reported the development of biological response by decreasing the time of implantation into tibias with no inflammation for 30 days post-implantation [161].Lala et al. prepared Mg-CHAP using a ball mill and showed that more than 95% of the cells were viable at low concentrations after 12, 24, and 48 h of incubation [162].

The Potential to Incorporate Dopants into HAP Structure Using Mechanochemical Synthesis
Lala et al. prepared Mg-doped-CHAP using a dry ball milling method in the existence of a hardened chrome steel vial for 10 h and 10 min pauses every 30 min with high energy planetary mill.They showed that crystallinity falls upon the increase of Mg ions.Furthermore, the lattice parameters were affected by Mg addition, and led to a decrease in the lattice parameters (a) and (c) [162].Shu et al. prepared CHAP using zirconia balls with a maximum diameter of 15 mm in addition to a speed of 1000 rpm.Then they obtained a nanofibrous structure based on CHAP and reported that HAP was completely synthesized after 1 h of milling.The nanofiber-based bell-milled HAP exhibited a diameter of 8 nm and a length in the range of 50-150 nm [163].Bulina et al. used water-cooled steel vials with a speed of 1200 rpm and BPR 1:20 for 30 min of milling.They estimated Si-HAP and revealed that substitution led to vacancy formation.In addition, the substitution increased the cell volume for the same milling conditions [164].Coreno et al. reported that a completed reaction is achieved for the preparation of CHAP after 6 h of milling under conditions of 20:1 BPR with high energy mixing and hardened steel balls.Further, the particle size was in nanosize in the range of 10-50 nm [150].In another work by Bulina et al., co-substituted HAP was prepared with the incorporation of Sr, and SiO 4 .They reported that SrSi-HAP prepared after 25 min of milling, which is estimated at 1300 rpm [165].Furthermore, Chaikina et al. prepared Zn-HAP, and the structure was obtained after 30 min of milling with a high speed of 1800 rpm and 20:1 BPR of steel balls [34].Moreover, Bulina et al. prepared Zn-HAP, and Cu-HAP at 1800 rpm after also 30 min in addition to Mg-HAP at 1200 rpm after 20 min [35,152].The successful synthesis of Zn, Cu, Mg, and CHAP holds great significance in the field of fertilizer applications.These materials play a crucial role in enhancing plant growth and productivity.The process of milling has proven to be highly effective in achieving the desired outcomes within a short duration, which makes it suitable for industrial-scale production.By optimizing the milling parameters, such as milling time and energy input, as shown in Figure 8, high-quality doped-HAP nanoparticles with desirable properties for efficient fertilizer utilization can be obtained.The ability to produce these materials through efficient milling techniques paves the way for their widespread adoption in the agricultural industry, leading to improved crop yields and sustainable agricultural practices.
proven to be highly effective in achieving the desired outcomes within a short duration, which makes it suitable for industrial-scale production.By optimizing the milling parameters, such as milling time and energy input, as shown in Figure 8, high-quality doped-HAP nanoparticles with desirable properties for efficient fertilizer utilization can be obtained.The ability to produce these materials through efficient milling techniques paves the way for their widespread adoption in the agricultural industry, leading to improved crop yields and sustainable agricultural practices.A schematic diagram for describing the formation of nanoparticles through short-time ball milling due to the existence of high frequency, high BPR, which leads to high input energy.These conditions should respect the lower and higher boundaries of the process.If the frequency or BPR is very high, the process stops, and no energy is transferred to the particles.Therefore, it can be done using high energy input, which is allowed to be transferred to the particles.

Summary and Future Perspectives
In this review, methods and examples of the doped-HAP synthesis and resulting properties are provided with an emphasis on the potential use of doped-HAP as a nutrient carrier.Several challenges have been identified in the preparation of HAP for nano fertilizers, especially using the mechanochemical method: One of the challenges is the limited research on doping HAP with various elements for agricultural purposes specifically.While elements, such as Zn, have been widely studied for doping HAP, other elements, such as molybdenum (Mo), have not been extensively explored.
"Nitrogen-doped" HAP using urea are of particular interest and recently enjoyed lots of publicity.However, it is unclear from the crystallographic point of view what these entities entail as urea molecules cannot simply substitute ions in HAP, and the resulting composite materials are likely inclusion compounds at best.The resulting crystal structure of these HAP-molecular urea complexes needs to be scrutinized in light of the reported enhanced nitrogen release properties.
The number of studies that specifically discussed the combination of mechanochemical synthesis and doping of HAP is relatively low.Therefore, the combination between mechanochemical doped HAP and agricultural applications is still to be developed.
Practical scalability, energy input, and overall sustainability of mechanochemically synthesized nutrient-doped-HAP compared to the conventional large-scale methods need to be developed.
Overall, more research efforts are required to explore the potential of different elements for doping HAP as a nano fertilizer by investigating the ionic release rate and use efficiency of the plants.This essentially requires the elements which are important for plant growth, including the macronutrients, such as K, Mg, S, and N.In addition to the micronutrients involving Fe, Mn, Cu, Zn, B, Mo, Ni, and Cl.This will allow for a broader range of nutrients to be incorporated into the HAP structure, catering to specific agricultural needs.Moreover, further research and innovation in doped HAP based on mechanochemical synthesis can lead to an investigation of new nano fertilizer based on HAP for large-scale production.In addition, exploring modifications to the process parameters, such as milling duration and speed, can enhance the efficiency and scalability of the HAP synthesis process.This will enable the production of high-quality nano fertilizers with controlled properties.

Figure 1 .
Figure 1.The number of publications over the years including (hydroxyapatite, and fertilizer) and (hydroxyapatite, fertilizer, and nano) in the keywords, title, or abstract from Scopus.

Figure 1 .
Figure 1.The number of publications over the years including (hydroxyapatite, and fertilizer) and (hydroxyapatite, fertilizer, and nano) in the keywords, title, or abstract from Scopus.

Figure 2 .
Figure 2. The crystal structure of (a) Hexagonal HAP, and (b) monoclinic HAP, where monoclinic has a c-axis longer than a and b axes.

Figure 2 .
Figure 2. The crystal structure of (a) Hexagonal HAP, and (b) monoclinic HAP, where monoclinic has a c-axis longer than a and b axes.

Figure 3 .Figure 3 .Figure 4 .
Figure 3. OH orientation in (a) hexagonal HAP with disordered OH dipoles, (b) in mono HAP with ordered dipoles of OH; where a and b represent the HAP unit cell axis.

Figure 4 .
Figure 4.The various calcium sites in HAP hexagonal crystal structure which are mainly two sites Ca(1) and Ca(2); where a and b represent the HAP unit cell axis.

3. 13 .
Co-Doping (Two or More Ions) Makshakova et al. studied the doping of HAP by Mg and/or CO 3

Figure 5 .
Figure 5.The schematic diagram for the general preparation pathway of the HAP coprecipitation method.Figure 5.The schematic diagram for the general preparation pathway of the HAP coprecipitation method.

Figure 5 . 12 Figure 6 .
Figure 5.The schematic diagram for the general preparation pathway of the HAP coprecipitation method.Figure 5.The schematic diagram for the general preparation pathway of the HAP coprecipitation method.Ceramics 2023, 6, FOR PEER REVIEW 12

Figure 6 .
Figure 6.The general preparation route of HAP using the hydrothermal method.Figure 6.The general preparation route of HAP using the hydrothermal method.

Figure 6 .
Figure 6.The general preparation route of HAP using the hydrothermal method.

Figure 7 .
Figure 7.The preparation steps of preparing HAP using the mechanochemical synthesis technique.

Figure 7 .
Figure 7.The preparation steps of preparing HAP using the mechanochemical synthesis technique.

Figure 8 .
Figure 8.A schematic diagram for describing the formation of nanoparticles through short-time ball milling due to the existence of high frequency, high BPR, which leads to high input energy.These conditions should respect the lower and higher boundaries of the process.If the frequency or BPR

Figure 8 .
Figure8.A schematic diagram for describing the formation of nanoparticles through short-time ball milling due to the existence of high frequency, high BPR, which leads to high input energy.These conditions should respect the lower and higher boundaries of the process.If the frequency or BPR is very high, the process stops, and no energy is transferred to the particles.Therefore, it can be done using high energy input, which is allowed to be transferred to the particles.

Table 1 .
The advantages and disadvantages of mentioned preparation methods.