1. Introduction
A wide variety of synthetic materials, including metals, polymers, and ceramics, are currently being investigated with regard to clinical applications as bone substitutes [
1,
2,
3,
4]. Of particular interest are sintered calcium phosphates, especially hydroxyapatite ceramics (Ca
10(PO
4)
6(OH)
2; HAp) [
5,
6], β-tricalcium phosphate ceramics (β-Ca
3(PO
4)
2; β-TCP) [
5,
7] and biphasic calcium phosphate ceramic composites (BCP) [
8,
9,
10]. This interest is primarily related to the fact that bone is largely composed of calcium phosphate apatite minerals. HAp implants are more crystalline than biological bone and so these implants exhibit greater chemical stability in vivo and tend not to degrade after implantation [
11,
12]. In the case that resorption does not occur in such implants, bone deformation [
13] and the risk of fracture around the replacement bone may increase in the long term [
11]. In contrast to HAp ceramics, β-TCP ceramic implants are bioabsorbable and easily replaced by new autogenous bone [
14,
15]. Because the solubility of these materials is close to that of living bone, these ceramics are also not soluble under physiological conditions (meaning a pH of 7.4 and temperature of 37 °C) [
16] and instead are typically absorbed by osteoclasts. The associated mechanism is believed to involve local acidification leading to the dissolution of the sintered β-TCP [
14,
17]. In addition to this osteoconductive capacity, sintered β-TCP has been shown to exhibit significant osteoinductive characteristics [
18]. For these reasons, sintered β-TCP is one of the most attractive bone replacement materials. In particular, porous β-TCP can be used as a scaffold material to promote the repair of bone tissue because it has suitable pore sizes to promote the growth of bone tissue within the porous material [
19,
20]. β-TCP also tends to contain macropores (100-600 µm) together with micropores (0.1-10 µm) and connections between these pores allow nutrients, blood, cells and biological tissue to enter, thereby promoting autogenous ossification. Many previous studies have demonstrated these effects [
21,
22,
23,
24,
25,
26].
Porous bodies can be produced with relative ease using foaming methods, in which a foaming agent is added and ambient air is entrained by mechanical agitation. In these techniques, the foam is typically stabilized by adding a surfactant based on the Gibbs-Marangoni effect, in which surfactant molecules compensate for surface tension gradients [
27]. Prior research by the authors has demonstrated that porous materials prepared by foaming methods exhibit biases in pore size and pore distribution. This bias occurs due to various factors involved in the formation of bubbles, such as flow due to gravity, surface viscosity, the contribution of capillary attraction at the plateau boundary to membrane liquid drainage and gas diffusion associated with the coalescence and redistribution of bubbles [
28]. An especially important factor is capillary attraction at the plateau boundary at which bubbles are adjacent to one another. The rate of flow to the plateau boundary has been shown to be correlated with the viscosity of the adsorption layer near the surface. It has also been reported that a higher viscosity will increase the thickness of the film between bubbles such that the bubble lifetime is increased [
29]. The Young-Laplace equation [
30] describing gas diffusion states that, in the case that two bubbles of different sizes come into contact with one another, the air from the smaller bubble (which has a higher internal pressure) will diffuse through the bubble film towards the larger bubble. As a consequence, smaller bubbles are absorbed by larger bubbles, causing the bubbles to merge and resulting in the disappearance of some bubbles.
The authors previously attempted to fabricate porous materials by adding a nonionic surfactant to β-TCP slurries as a foaming agent together with cellulose nanofibers (CNFs) as a thickening agent and bubble stabilizer. These additives were employed in conjunction with a physical foaming method in which pores were formed by mechanical stirring. This prior work showed that, in trials using a conventional foaming method, pores were uniformly distributed throughout the specimens and porous bodies with high mechanical strength of the slurry hardened body could be obtained. It was also found that the gravity-induced outflow of water into the lower layers and the coalescence of bubbles prevented the pores from concentrating at the top of the sample and prevented the pores from coarsening [
31]. However, using this method, it was challenging to properly disperse the CNFs and these fibers tended to agglomerate during drying in association with heating. This agglomeration resulted in sintered porous bodies having irregular pore shapes rather than spherical pores. In later work, porous materials were generated using CNFs in which some of the -OH groups on the CNFs were acetylated [
32] or replaced by phosphite groups [
33]. The resulting porous materials were found to have spherical pores. Furthermore, we prepared porous β-TCP by a similar method, instead of CNFs with the thickening effect of carboxymethylcellulose ammonium salt, and found that porous β-TCP could be prepared with suitable pore size, porosity and mechanical strength for bone filling material, without any difference in pore size distribution according to the sample location [
34].
In the present study, foamed β-TCP slurries were rapidly frozen and solidified, after which dried porous materials were obtained from the frozen specimens, from which water was removed by freeze-drying. This was an alternative to slower heat-based drying. The freeze-drying method used in this work was superior in that the water in the slurry was rapidly frozen at reduced pressure, after which the material was returned to room temperature such that the water evaporated without going through the liquid phase. This process allowed the samples to be dried without destroying the pore structure of the frozen slurry. The effects of varying the fabrication conditions were investigated with a focus on quick freezing. Porous specimens were also fabricated using nonionic surfactants with different hydrophilic/lipophilic balance (HLB) values [
35] to control the pore size, and the physical properties of the resulting porous materials were evaluated.
2. Experimental
2.1. Preparation of Raw Powder
In this procedure, 11.32 g of calcium carbonate (CaCO3, purity 99.0%, Fujifilm Wako Pure Chemicals, Osaka, Japan) and 39.51 g of calcium hydrogen phosphate dihydrate (CaHPO4⋅2H2O, purity 99.5%, Fujifilm Wako Pure Chemicals, Osaka, Japan) were mixed in a Ca/P molar ratio of 1.50. This mixture was subsequently combined with 450 ml pure water heated to 80 °C, and then ground using zirconia balls (Nikkato Corporation, Osaka, Japan; YSZ balls, 400 g) in a zirconia pot (As One Corporation, Osaka, Japan) having dimensions of 100 mm (width) × 130 mm (height) for 24 h. Each mixed sample was then spread in a plastic container and dried at 70°C for 24 h. Following this, each dried specimen was ground using an agate mortar and then calcined in an electric furnace by heating to 750 °C at 3 °C/min with a hold time of 10 h. The calcined powders were characterized by X-ray diffraction (XRD, MiniFlex 600: Rigaku Corporation, Tokyo, Japan), Fourier transform infrared (FT-IR) spectroscopy (FT-IR-4200: JASCO Corporation, Tokyo, Japan) and scanning electron microscopy (SEM, VE-7800: KEYENCE Corporation, Tokyo, Japan).
2.2. Preparation of Porous Materials by Freeze-Drying
In each trial, 30 g of the calcined powder was mixed with 30 g of a 3.0 wt% dispersion of CNFs in water (fiber width: 0.069 µm, fiber length: 9.66 µm: Daio Paper Corporation, Shikoku, Japan) together with 30 ml of ammonium polyacrylate solution 70-110 (molecular weight: approximately 10,000, Fujifilm Wako Pure Chemicals, Osaka, Japan), diluted to 5 vol% and dispersed in cold water using a hand blender (HB-502WJ: Cuisinart, Tokyo, Japan) in an ultrasonication device (As One, Tokyo, Japan) filled with ice water for 5 min. Subsequently, a 4 ml quantity of a polyoxyethylene alkyl ether nonionic surfactant (BT-7, Nikko Chemicals Corporation, Tokyo, Japan) was added as a foaming agent. Foaming was then carried out by mechanical agitation under the same conditions as used to produce the original dispersions. The resulting foam was tapped at a height of approximately 20 mm and poured into a polypropylene mold (25 mm in length × 25 mm in width × 25 mm in height). Each specimen was subsequently frozen by being suspended above the surface of liquid nitrogen in a stainless-steel dewar (at which location the ambient temperature was −180 °C), as shown in
Figure 1. The sample was then pre-frozen by direct immersion in the liquid nitrogen (at a temperature of −196 °C) and then stored in a −40 °C freezer. Following this, each pre-frozen sample was freeze-dried using an FDU-1200 apparatus (EYELA, Tokyo, Japan) at −40 °C and 20 Pa. Each of the dried specimens was then sintered in an electric furnace by heating at a rate of 5 °C/min to 180 °C and holding at that temperature for 4 h followed by a 4 h hold at 300 °C, a 4 h hold at 400 °C and a 4 h hold at 1100 °C.
The time spans over which each specimen was held above the liquid nitrogen surface and the immersion times in the liquid nitrogen are given in
Table 1. Note that each sample is referred to herein using the notation x-y, where x is the hold time above the liquid nitrogen [s] and y is the immersion time in the liquid nitrogen [s]. To allow a comparison of drying methods, samples were also prepared by heat-drying at 40 °C for 72 h using an NDO-450ND apparatus (EYELA, Tokyo, Japan) after which the same sintering program was applied.
2.3. Preparation of Porous Materials Using Nonionic Surfactants with Different HLB Values
In these procedures, 30 ml of a 5 wt% ammonium polyacrylate solution, 30 g of a 3.0 wt% dispersion of CNFs in water and 30 g of the calcined β-TCP powder were mixed for 5 min using a hand mixer in an ultrasonication device filled with ice water. Following this, 4 ml of a polyoxyethylene alkyl ether surfactant (either BT-5, BT-7, BT-9 or BT-12, Nikko Chemicals Corporation, Tokyo, Japan) was added as a foaming agent and the specimen was foamed for 5 min. After tapping, the mixture was poured into a polypropylene mold and first frozen by holding for 30 s above the surface of liquid nitrogen then pre-frozen by direct immersion for 20 s in the liquid nitrogen. Samples were subsequently stored in a freezer. Following this, each frozen specimen was freeze-dried at −40 °C and 20 Pa. The freezing conditions were 30-20 using the notation described above. The dried samples were sintered using the method described in the previous section. The HLB values of the foaming agents used in this study were 10.5, 12.0, 13.5 and 14.5 for the BT-5, BT-7, BT-9 and BT-12, respectively [
35].
2.4. Evaluation of Porous Materials
Some of the porous materials were crushed and the resulting powders were characterized by XRD and FT-IR spectroscopy. The samples were also cut into sections along the vertical direction and the pore structures of these sections were evaluated by SEM. These sections comprised the top, center and bottom parts with respect to the height direction of each sample. From the SEM images, the largest diameters of 300 randomly selected pores were determined for each specimen using the ImageJ image processing software. In addition, 30 pores were randomly selected from SEM images of the top, center and bottom sections of each specimen and the areas and perimeters of these pores were ascertained. The corresponding circularity value for each pore (having a value between 0 and 1) was then calculated as
The porosity and bulk density of each material were determined by the Archimedes method following the procedure in the JIS R 1634 standard. Compressive strength tests based on the JIS R 1680 standard were carried out using a universal testing machine (Little Senstar, Tokyo Testing Machinery, Tokyo, Japan). Pore size distributions and pore volumes were obtained using the mercury injection method, employing an Autopore IV 9520 instrument (Shimadzu Corporation, Kyoto, Japan).
3. Results and Discussion
Figure 2 presents XRD patterns obtained for sintered samples prepared with various hold and immersion times during pre-freezing. The diffraction peaks generated by all samples can be attributed to β-TCP (ICDD No. 055-0898). No peaks related to by-products such as hydroxyapatite appeared, indicating that each of the samples comprised a single β-TCP phase.
Figure 3 provides the FT-IR spectra of the same samples as shown in
Figure 2. These spectra exhibit peaks ascribed to the bending vibration of PO
4 groups in β-TCP at 430 cm
-1 (ν
2) and 560-600 cm
-1 (ν
4) and the stretching vibration of PO
4 groups at 960 cm
-1 (ν
1), 1020 cm
-1 (ν
3) and 1120 cm
-1 (ν
3). The absence of peaks related to the potential by-products hydroxyapatite (that is, -OH peaks) and calcium pyrophosphate (meaning P-O-P peaks) in these spectra provides further evidence that each material was pure β-TCP [
36].
Figure 4 presents photographic images of the samples. Here,
Figure 4a shows the material immersed in liquid nitrogen for 30 s and the sample is seen from above. In this case, many cracks were formed during the freezing process and it was difficult to maintain the shape of the specimen. The image in
Figure 4b shows a sample immersed in liquid nitrogen but not held above the liquid nitrogen surface, as viewed from the side. Internal bubbles evidently exploded at the top of the sample and the original shape of the specimen was changed. These effects are attributed to the sudden temperature change caused by immersion in the liquid nitrogen. The image in
Figure 4c shows a sample immersed in the liquid nitrogen for 20 s or longer with a hold time of more than 30 s above the liquid nitrogen surface, seen from the side. Compared with the specimen in
Figure 4b, the shape of the sample in
Figure 4c was hardly distorted, indicating that the sample morphology could be maintained by using hold and immersion times of 20 s or longer.
Figure 5 shows cross-sectional SEM images of a sample immersed in the liquid nitrogen for 10 s with no hold time above the liquid nitrogen and a sample held above the liquid nitrogen for 240 s but with no immersion. The 240-0 sample exhibits spherical pores along with many crack-like pores that distort the spherical pore shapes. Holding this sample for 240 s above the liquid nitrogen caused it to freeze more slowly due to the lower thermal conductivity of the gaseous nitrogen. This, in turn, promoted the growth of ice crystals in the sample, and the sublimation of this ice during subsequent freeze-drying is thought to have produced the crack-like pores. The SEM images of the sample immersed in the liquid nitrogen for 10 s confirm the formation of spherical pores in the top and bottom parts of the sample whereas the center part contained numerous cracks similar to those in the 240-0 specimen. These results suggest that a 10 s immersion in the liquid nitrogen was not sufficient to freeze the specimen such that subsequent cooling froze the center part.
Figure 6 presents cross-sectional SEM images of the sample immersed in the liquid nitrogen for 20 s and the samples produced with heat-drying. The latter material exhibits many distorted sphere-like pores with sizes of 200-300 µm, whereas the former sample shows spherical pores with sizes of 50-200 µm together with an interconnected pore structure. These findings indicate that the preparation of materials with spherical pores is possible based on pre-freezing with liquid nitrogen followed by freeze-drying. The sample immersed in liquid nitrogen for 20 s was evidently sufficiently frozen as no cracked pores were observed. Hence, ice crystals were not grown in this material as had occurred in the 240-0 sample.
Table 2 summarizes the pore circularity values calculated from the SEM images of the 30-20 sample and the heat-dried material. Note that the pore circularity increases as the values approach 1.0. The values for the samples made using heat-drying differed slightly, with an average of 0.85. In contrast, the 30-20 sample produced using freeze-drying showed no difference in circularity between different sites, with an average value of 0.97. These data confirm that improved circularity was obtained relative to the heat-dried sample.
Figure 7 provides pore diameter histograms obtained from analysis of SEM images of the 0-20, 10-20, 20-20, 30-20 and Heat-drying samples. Here, the horizontal axis represents the pore diameter while the vertical axis indicates the frequency. These data confirm that the pore diameters for all samples were distributed over a wide range from 50 to 200 µm. The mean and median pore sizes for the samples immersed in liquid nitrogen were all approximately 100 µm while the values for the Heat-drying sample were drastically increased. This discrepancy is attributed to the fact that cracked pores in the specimens disrupted the pore shapes and also indicates an effect whereby pores smaller than 50 µm were very difficult to observe. It is also apparent that the pore sizes for the foams did not change when the freezing conditions were changed, presumably because the foams were frozen instantaneously.
The formation of pores within these materials using the different drying methods is illustrated in
Figure 8. Drying with hot air required a longer time span such that bubbles had more time to undergo coalescence. During sintering, the CNFs were burnt off and produced pores and so the samples generated using heat-drying were more likely to have distorted sphere-like pores. In contrast, freeze-drying after instantaneous pre-freezing prevented the coalescence of bubbles, resulting in the formation of many pores with smaller diameters. The freezing of the CNFs while they were dispersed in the slurry prevented their agglomeration, providing more spherical pores [
31].
Figure 9 summarizes the open porosity, closed porosity and bulk density values for each sample as ascertained using the Archimedes method. The open porosity for each of the freeze-dried specimens was approximately 70%, indicating that highly porous materials were prepared. The closed porosity values for all samples were less than 1% of the total porosities, demonstrating that these bodies should be more permeable than those prepared using heat-drying. Freeze-drying evidently increased pore connectivity because the CNFs were dispersed throughout each sample and generated pores as they were removed during sintering [
31].
The compressive strength test results for all samples are shown in
Figure 10. The compressive strength of each freeze-dried product exceeds 2 MPa, indicating that sufficient mechanical strength has been achieved to withstand handling. There was a correlation between the porosity of each specimen and its compressive strength, such that more porous materials had lower compressive strength. The 240-0 sample had the lowest compressive strength, which was assumed to be due to the crack-like pores observed in the SEM images. The specimens made using heat-drying had distorted sphere-like pores and therefore exhibited reduced compressive strength despite their low porosity values. Each of the materials containing only spherical pores was found to exhibit both high porosity and good compressive strength, indicating that spherical pores promoted both mechanical strength and porosity.
Figure 11 shows the results of pore distribution measurements by the mercury injection method for samples immersed in liquid nitrogen for 30 s, immersed for 20 s and freeze-dried or heat-dried [
31]. The heat-dried sample was found to contain micropores with a maximum pore size of 0.46 µm and a wide range of macropores with sizes of 50-200 µm. The freeze-dried samples exhibited micropores with a maximum pore size of 1.24 µm and macropores with a maximum pore size of 25.8 µm in the 10-80 µm range. These data indicate a transition of the micropores to larger sizes and of the macropores to smaller sizes compared with the heat-dried samples. The micropores likely became larger due to the absence of shrinkage in the case of the freeze-dried materials, as the pore shapes changed only minimally when the specimens were dried while frozen. The macropores are thought to have shifted to smaller sizes because bubbles did not coalesce during drying. These data also suggest that the pores in these materials were connected to the exterior surfaces.
In additional trials, the HLB value [
34] for the nonionic surfactants used as foaming agents were varied. The HLB value is an indicator of the hydrophobicity/hydrophilicity of each surfactant, with a higher value reflecting a more hydrophilic compound. A nonionic surfactant having a low HLB value typically shows reduced foaming properties and generates smaller pores. Cross-sectional SEM images of porous β-TCP bodies prepared using nonionic surfactants with different HLB values are shown in
Figure 12. Pores less than 50 µm in size were predominant in the sample made with the BT-5 (having an HLB of 10.5) while larger pores were obtained when using the BT-7 (HLB 12.0), BT-9 (HLB 13.5) and BT-12 (HLB 14.5). Materials with more macro-sized pores having larger pore diameters were obtained by using foaming agents with higher HLB values.
Figure 13 provides pore size histograms for each sample as obtained by analysis of the SEM images. The average pore size and median diameter in the BT-5 sample were 88 and 73 μm, respectively, while the values obtained with the BT-7 were 100 and 93 μm, respectively, those with the BT-9 were 140 and 128 μm, respectively, those with the BT-9 were 140 and 128 μm, respectively, and those with the BT-12 were 136 and 120 μm, respectively. These results show that the addition of foaming agents with high HLB values increased both the average pore size and the median diameter. In particular, the use of the BT-9 or BT-12 came close to producing saturation and gave macro-sized pore diameters of more than 180 µm.
The bulk densities and porosities of the porous materials prepared by adding foaming agents with different HLB values are summarized in
Figure 14. The bulk density and porosity values obtained with the BT-5 and BT-7 were approximately 0.9 g/cm
3, along with open porosity values of 70%. However, the BT-9 and BT-12 provided values of approximately 0.6 g/cm
3 with open porosities of 80%. In both cases, the closed porosity of the sample was almost 0%.
The compressive strengths of the various materials are provided in
Figure 15. The compressive strengths of the specimens evidently decreased with the use of foaming agents having higher HLB values. The BT-9 and BT-12, both of which provided greater porosity, produced materials with strengths on the order of 1.9 MPa. Based on the above results, it appears that a foaming agent with an HLB value intermediate between those of BT-7 and BT-9 would be ideal and could produce materials simultaneously exhibiting both high porosity and good compressive strength.
Figure 1.
Schematic illustration showing pre-freezing process.
Figure 1.
Schematic illustration showing pre-freezing process.
Figure 2.
XRD patterns for porous samples fabricated in this work.
Figure 2.
XRD patterns for porous samples fabricated in this work.
Figure 3.
FT-IR spectra of porous samples fabricated in this work.
Figure 3.
FT-IR spectra of porous samples fabricated in this work.
Figure 4.
Photographic images of (a) 30-30 (top view), (b) 0-30 (side view) and (c) 30-20 (side view) specimens.
Figure 4.
Photographic images of (a) 30-30 (top view), (b) 0-30 (side view) and (c) 30-20 (side view) specimens.
Figure 5.
SEM images of cross-sections of (a) 0-10, (b) 20-10, (c) 30-10 and (d) 240-0 specimens, showing (i) top, (ii) center and (iii) bottom sections.
Figure 5.
SEM images of cross-sections of (a) 0-10, (b) 20-10, (c) 30-10 and (d) 240-0 specimens, showing (i) top, (ii) center and (iii) bottom sections.
Figure 6.
SEM images of cross-sections of (a) 10-20, (b) 20-20, (c) 30-20 and (d) heat-dried specimens, showing (i) top, (ii) center and (iii) bottom sections.
Figure 6.
SEM images of cross-sections of (a) 10-20, (b) 20-20, (c) 30-20 and (d) heat-dried specimens, showing (i) top, (ii) center and (iii) bottom sections.
Figure 7.
Pore size histograms for various porous β-TCP specimens.
Figure 7.
Pore size histograms for various porous β-TCP specimens.
Figure 8.
Schematic illustration showing pore formation in foams using different drying methods.
Figure 8.
Schematic illustration showing pore formation in foams using different drying methods.
Figure 9.
Bulk density and porosity values for porous β-TCP specimens as determined using Archimedes method. The directions of the arrows indicate the relationship between the axes in the figure. The numerical values for porosity indicate closed porosity in the upper part and open porosity in the lower part.
Figure 9.
Bulk density and porosity values for porous β-TCP specimens as determined using Archimedes method. The directions of the arrows indicate the relationship between the axes in the figure. The numerical values for porosity indicate closed porosity in the upper part and open porosity in the lower part.
Figure 10.
Compressive strength for porous β-TCP specimens fabricated in this work.
Figure 10.
Compressive strength for porous β-TCP specimens fabricated in this work.
Figure 11.
Pore distributions for porous β-TCP specimens fabricated by different drying methods as determined by mercury injection method.
Figure 11.
Pore distributions for porous β-TCP specimens fabricated by different drying methods as determined by mercury injection method.
Figure 12.
SEM images of cross-sections of porous β-TCP specimens fabricated using nonionic surfactants (a) BT-5, (b) BT-7, (c) BT-9 and (d) BT-12, showing (i) top, (ii) center and (iii) bottom sections.
Figure 12.
SEM images of cross-sections of porous β-TCP specimens fabricated using nonionic surfactants (a) BT-5, (b) BT-7, (c) BT-9 and (d) BT-12, showing (i) top, (ii) center and (iii) bottom sections.
Figure 13.
Pore size histograms for porous β-TCP specimens fabricated using different nonionic surfactants.
Figure 13.
Pore size histograms for porous β-TCP specimens fabricated using different nonionic surfactants.
Figure 14.
Bulk density and porosity values for porous β-TCP specimens fabricated using different nonionic surfactants. The directions of the arrows indicate the relationship between the axes in the figure. The numerical values for porosity indicate closed porosity in the upper part and open porosity in the lower part.
Figure 14.
Bulk density and porosity values for porous β-TCP specimens fabricated using different nonionic surfactants. The directions of the arrows indicate the relationship between the axes in the figure. The numerical values for porosity indicate closed porosity in the upper part and open porosity in the lower part.
Figure 15.
Compressive strength of porous β-TCP specimens fabricated using different nonionic surfactants.
Figure 15.
Compressive strength of porous β-TCP specimens fabricated using different nonionic surfactants.
Table 1.
Hold times for samples above liquid nitrogen surface and immersion times in liquid nitrogen.
Table 1.
Hold times for samples above liquid nitrogen surface and immersion times in liquid nitrogen.
Table 2.
Pore circularity values for porous β-TCP specimens obtained using different drying methods.
Table 2.
Pore circularity values for porous β-TCP specimens obtained using different drying methods.