Characterization of the polycaprolactone melt crystallization: complementary optical microscopy, DSC, and AFM studies.
The phenomenon of crystallization of polymers has been studied since the start of the macromolecular science. Despite the great efforts spent by the scientific community on this topic, the formation of crystalline structures between highentan-gled polymer chains is still debated.
A complete investigation of the crystallization process requires the use of several analytical techniques, since a single technique does not supply sufficient understanding of the complex multiscale transformation. Heat transition involved in the process is typically measured by differential scanning calorimeter (DSC-DTA), structures formation is generally monitored by microscopy (optical, AFM, TEM, or SEM), and order alteration is observed by X-ray diffraction and infrared spectroscopy, whereas changes in flow properties are commonly measured with rheological tests [1, 2]. However, an accurate comparison of the quantities measured on individual instruments is not straightforward, since not only the scale but also the samples conditions are always different.
The application of atomic force microscopy (AFM) to polymer system is often devoted to the observation of a solidified sample after quenching the crystal structure, as a static subsequent step. Most dynamic hot stage studies of morphological development in situ have been performed using optical microscopy, with its limited resolution.
Recently, the possibility to follow real-time dynamics of processes at nanoscale, even if limited on time scale, is opened to AFM. Since the initial in situ observations of the polymer crystallization reported by Pearce and Vancso [3-5] and by Schultz and Miles  several additional real-time AFM studies were reported [7-14]. Some of these studies were performed by controlling the temperature of the sample [10,11] and sometimes the temperature of the AFM tip [13,14].
Pearce and Vancso [3, 4] first reported the consistency of in situ AFM results with conventional optical microscopy results. They found that growth rates of individual lamellae agree with results obtained at lower magnification by optical microscopy. Similar growth rates were observed at the growth front by Hobbs et al.  whereas observations at the scale of the lamella revealed for the first time that dominant lamellae do not grow at a constant rate. Lamellae grow at sporadic or constant rates depending on the presence or not of surrounding lamellae competing at the growth front 10,11,14].
In situ AFM studies gave the opportunity to observe developing morphologies at a lamellar scale, from the generation of the primary stable nucleus  to the development of dominant lamellae, followed by the birth of branching lamellae that splay, leading to the spherulitic morphology [4, 6, 8] with, finally, impingement of adjacent spherulites .
These observations highlighted the roughness of the growth front, seen as being smooth at the larger scale of optical microscopy .
In particular for the investigation of material transformations that are associated with changes in rheological and thermal properties, such as melting and crystallization, an accurate comparison between these techniques is crucial. In the literature, a few attempts for combining morphology, calorimetry, and rheology can be found .
Martins et al.  developed a shear differential thermal analyzer that allowed for applying controlled shear pulses during isothermal or nonisothermal solidification of polymers. Hereto, a capillary rheometer and a differential thermal analyzer were combined into a single device. The instrument has two capillaries, one for the investigated sample and one for a reference material, which are placed inside two separate identical ovens.
Nagatake et al.  modified a commercial rheometer by adding two thermocouples on the base plate, which allowed the measurement of differences between sample and base plate temperatures.
In this paper, the investigation of the morphology development of melt-grown crystal lamellae of polycaprolactone by real-time hot stage AFM is reported. The morphological insights are also investigated by using polarized optical microscopy (POM), differential scanning calorimetry, and rheology tests.
2.1. Material. The polymer investigated in this work is polycaprolactone ([M.sub.n] 70000-90000 by GPC, [M.sub.w]/[M.sub.n] < 2, density 1.145 g/mL at 25[degrees]C), supplied by Sigma Aldrich.
2.2. Differential Scanning Calorimetry. A differential scanning calorimeter DSC 822 from Mettler Toledo Inc. was used for determination and measurement of the calorimetric behavior of PCL. The calibration of the temperature was done with the extrapolated onset temperature of the phase transition of indium. The samples with a weight of about 10 mg were put into an aluminum pan. All the experiments were conducted in nitrogen (flow rate 50mL/min), in order to prevent oxidative degradation.
2.3. Optical Microscopy. A PCL thin film, of about 100 [micro]m, is sandwiched between two cover glasses and this is observed in light transmission using an optical Olympus BX51 Microscope equipped with digital camera. To apply the same thermal protocol to the glass/polymer system a special calorimeter (Linkam DSC600, calibrated with indium, cooled with liquid nitrogen, and purged with nitrogen ) was adopted for the experiments. The morphology and the isothermal growth of the crystalline structures were observed using the microscope coupled to the Linkam DSC600 during the calorimetric tests.
2.4. Atomic Force Microscopy. In order to characterize the local morphology of the samples a Veeco Dimension 3100 AFM with Nanoscope III controller was adopted. All the sample images were kept using the tapping mode in order to minimize the interaction with the sample surface. Two tests were conducted with the same thermal history, but analyzing different surfaces: one test explored a square surface having a side of 10 [micro]m and the other one a square surface having a side of 30 [micro]m. In both experiments, the scan rate was 2 Hz and 256 sample lines. The Linkam hot stage THMS600 is used for conditioning the sample during the in situ AFM characterization. The equilibrium melting temperature of PCL is relatively low, about 60[degrees] C, so that it was possible to perform the AFM acquisition without altering the linearity of the AFM piezo.
2.5. Rheological Tests. Rheological experiments were per formed using a stress controlled rotational rheometer (HAAKE RheoScope) equipped with 25 mm parallel plates and using a gap thickness of 0.1mm. To follow both the cooling process and the isothermal crystallization, the evolution of the mechanical properties has been monitored by an oscillatory test at the constant stress, [sigma] = 100 Pa, and constant frequency, [omega] = 0.1 rad/sec. The choice of these values is crucial in order to reduce, as much as possible, the disturbance during the crystallization and ensure that the whole crystallization process occurs in the range of sensitivity of the equipment. To verify whether these parameters could affect the experimental results, tests with stress values ranging from 1 to 1000 Pa have been performed.
Optical observations were made directly in the rheometer by a CCD camera (HAAKE RheoScope) equipped with a 5x magnification objective.
2.6. Methods. The same thermal history has been applied to all the experiments. In particular, an annealing treatment has been conducted at T = 120[degrees]C for 15 min. The cooling step from the annealing temperature ([T.sub.a]) to the crystallization temperature ([T.sub.c]) has been performed in two steps. The first step from [T.sub.a] to the melting temperature [T.sub.m] = 60[degrees]C was performed at a rate compatible with the limitations of the different experimental techniques (-10[degrees]C [min.sup.-1] for DSC, AFM, and POM; -3[degrees]C [min.sup.-1] for rheological tests) whereas the second step from [T.sub.m] to [T.sub.c] was performed at a rate of -0.5[degrees]C [min.sup.-1]. The temperature profile has been verified by an independent thermocouple placed inside the molten sample. The temperature sensors of the rheometer, the calorimeter, and the hot stage were carefully calibrated to minimize temperature differences among the instruments.
3. Results and Discussion
3.1. DSC Experiments. A preliminary calorimetric test was performed in order to determine the relevant temperatures for the crystallization phenomena of PCL. The sample was heated to 120[degrees]C at 30 C/min, kept at this temperature for 15 minutes, and then cooled down to 60[degrees]C at -10 C/min and from 60 to 25[degrees]C at -0.5 C/min. The result of this test is reported in Figure 1.
The calorimetric plot evidences an exothermic peak due to the crystallization of PCL, which starts at 44[degrees] C, has a maximum at 39[degrees]C, and ends at 34[degrees]C.
A series of isothermal crystallization tests were thus carried out at temperatures between 50[degrees]C and 43[degrees]C. The corresponding exothermic peaks are reported in Figure 2 as a function of the time of the isothermal step. As expected, by increasing the crystallization temperature the peaks decrease in height and spread along the time axis. The area delimited by the peaks, characteristic of the heat released during the crystallization, is almost constant for all the experiments.
The crystallization heat, obtained by the measurement of the area under the exothermic peak, can be transformed into the relative degree of crystallinity ([[alpha].sub.DSC](t)) by the division of the heat that develops at each crystallization time t ([DELTA][H.sub.t]) by the total area under the exothermic peak, that is, the total heat ([DELTA][H.sub.0]) generated up to the complete crystallization:
[[alpha].sub.DSC] = (t) = [[integral].sup.t.sub.0] (dH/dt) dt/[[integral].sup.[infinity].sub.0] (dH/dt)dt = [DELTA][H.sub.t]/[DELTA][H.sub.0]. (1)
Figure 3 displays the development of [[alpha].sub.DSC] as a function of the isothermal time at several crystallization temperatures. A shift of the curves toward lower times and an increment of the slope of their linear portion can be seen as the crystallization temperature becomes lower. This implies that a higher undercooling degree leads to a higher crystallization rate.
An approximation of the crystallization rate can be made by the overall crystallization rate (1/[t.sub.0.5]), where [t.sub.0.5] (half the crystallization time) is the time at which [[alpha].sub.DSC] (t) approaches the value of 50%. This parameter is proportional to both the primary nucleation rate and the crystal/spherulite growth . The induction time ([t.sub.0]), that is, the time at which [[alpha].sub.DSC](t) starts to increase, is also clear from Figure 3. As expected, [t.sub.0] increases exponentially with the crystallization temperature.
3.2. POM Experiments. The morphology developing during the crystallization of PCL was investigated by polarized optical microscopy (POM). The micrographs were acquired at suitable time intervals with crossed polarizer during the same isothermal tests performed during the DSC tests.
The development of crystalline structures at 45[degrees]C is shown in the optical micrographs reported in Figures 4, 5, and 6. On observing the morphology evolution by the polarized optical microscope, it was possible to measure the evolution of birefringence. Birefringence was then normalized according to the following equation:
[[alpha].sub.POM] (t)= [DELTA]n (t) - [DELTA][n.sub.0]/[DELTA][n.sub.[infinity]] - [DELTA][n.sub.0], (2)
in which [DELTA]n is the measured birefringence at a given time, [DELTA][n.sub.0] is the birefringence at the start of the isothermal test (which is nearly zero), and [DELTA][n.sub.[infinity]] is the birefringence at the end of crystallization (when a plateau is reached). Being birefringence essentially proportional to crystallinity, [[alpha].sub.POM] (i) defined by (2) can be used to monitor the evolution of crystallization in thin samples.
In Figure 7, the evolution of [[alpha].sub.POM] is plotted as a function of time for a series of isothermal tests similar to those presented in Figure 2. As the crystallization temperature increases, the sigmoidal curves are shifted parallel to each other on a linear-log scale. At the highest temperature used here, the time to reach the birefringence plateau has been increased by an order of magnitude.
It is important to point out that evolution intensity expressed by (2) does not provide a direct measure of the degree of crystallinity. Polymer crystallizing systems consist of many birefringent spots, arranged in the light path either side by side or in series. Ziabicki  argued that single-crystal formula (2) should be applied only to very dilute systems and/or very thin samples, in which probability of the appearance of more than one crystal in the light path is negligible. On the other hand, Binsbergen  developed a model of stacks of plates with small optical retardation, claiming applicability of the single-crystal formula to nondilute systems. More in general the depolarization ratio is not proportional to crystallinity alone but to crystallinity multiplied by a function of average plate dimensions . However, assuming directly observable quantity aPOM as a crystallization characteristic per se, information about structure formation (induction time) and structural changes can be obtained.
3.3. Rheological Experiments. In Figures 8(a) and 8(b) the evolution of the storage modulus (G'), the loss modulus (G") and the complex viscosity ([eta]*) versus time during both the cooling (after the annealing process at T = 120[degrees]C) and the isothermal crystallization is reported together with the corresponding temperature profiles for two different quiescent crystallization tests.
The experiment shows that both moduli (G' and G") are susceptible to structural changes in the polymer. As an example, G' after a first increase during the cooling step, which can be ascribed to the effect of temperature, shows a plateau at T = [T.sub.c], followed by a second increase representative of the crystallization process. Finally, a second plateau is found when crystallinity reaches a maximum value.
From the rheological experiment, the induction time for the isothermal crystallization can be detected at the time when the storage modulus abruptly increases. The transformed crystallized fraction [[alpha].sub.RHEO] (i) has been estimated from the time dependent G'(t) values according to the following equation :
[[alpha].sub.RHEO] (t) = G' (t) - [G'.sub.0]/[G'.sub.[infinity]] (t) - [G'.sub.0], (3)
where [G'.sub.0], G'(t), and [G'.sub.[infinity]] are the melt stiffness values at time 0, t, and infinity, respectively.
The results for the structure development during the isothermal quiescent crystallization have been obtained at different crystallization temperatures. In Figures 9(a) and 9(b) the evolution of [[alpha].sub.RHEO] to the isothermal crystallization at 45[degrees]C and 47[degrees]C, respectively, is reported.
The time scale starts at the end of the cooling step for both the rheological and the morphological tests. The times read from the figures are, then, directly related to the evolution of the crystallization in isothermal condition.
The morphological investigation shows that spherulites develop and grow until the impingement occurs. As expected, the number of nuclei increases as the isothermal crystallization temperature decreases. What is interesting is the comparison between the structure development and the G' trend (Figures 9(a) and 9(b)). Indeed, in all the cases it has been observed that spherulites appear and grow when the storage modulus is still constant at the first plateau value (representative of G' for the undercooled melt at [T.sub.c]). This result is in agreement with the results reported by Acierno et al. , Bove and Nobile , and Pogodina et al. . The increase in G' during time after the impingement is still related to the crystallization process itself.
3.4. AFM Experiments. The same thermal protocol adopted before was replicated during AFM measurements. In that case, however, the recording of the data starts in parallel with the second cooling step (T < 60[degrees]C). In contrast to that observed with the calorimetric experiments, the first signal of structures formation was observed during the second cooling step, at temperatures just below 50[degrees]C. For that, it was impossible to carry out isothermal AFM tests at temperature lower than 50[degrees] C.
The evolution of the samples surface morphology was recoded in terms of distance from a reference position (height) and the characteristic oscillation (amplitude and phase) of the AFM cantilever during the observation times. The three acquired signals are strictly related to the structure formation on the sample surface during the crystallization phenomena. However, in our case, the amplitude signal was found the most sensible to the change in surface morphology.
For the present analysis, two different tests were considered: one with an observed square area with a side of 10 [micro]m and another with wider area having a side of 30 [micro]m. As shown in Figure 10, some evidence of crystalline structures appears on the surface after about 1500 s in isothermal conditions.
The crystalline structure development was followed by AFM while growing and it was found that the growth front propagates by the advancement of primary arms, with the regions between the arms filled in by secondary arms growing behind the primary front .
During the isothermal step new crystal aggregates were not formed but just the development of the previously formed structures was evident. This observation is coherent with the hypothesis of nucleation dependence on temperature and thus during the isothermal step no new stable nuclei are formed.
The development of new crystalline structures on the sample surface was followed by means of the evolution of the average value of the amplitude signal, [bar.A](t), during the whole test.
In particular, by assuming as initial value of the averaged amplitude ([[bar.A].sub.0]) the value that [bar.A](t) assumes at starting crystallization time, the variations of the averaged amplitude for each time ([[bar.A].sub.(t)] - [[bar.A].sub.0]) divided by the total variation, recorded during the whole test, can be transformed into a variation index ([[alpha].sub.AFM](t)):
[[alpha].sub.AFM] (t) = [[integral].sup.t.sub.0] ([bar.A] (t) - [[bar.A].sub.0])dt/[[integral].sup.[infinity].sub.0]( [bar.A](t) - [[bar.A].sub.0])dt. (4)
In Figure 11 is reported, for both tests conditions, the evolution of [[alpha].sub.AFM](t) as function of the isothermal crystallization time.
At early stage of the crystallization process, a small part of the investigated surface is covered by crystalline structure, so the variations of the averaged amplitude signal are very low. As the crystallization proceeds, the variations of the averaged amplitude become larger until a large part of the surface is covered by the crystalline structure and therefore the variations of the averaged amplitude are constant.
As expected, even on the same sample, the morphology appears very different by changing the surface analyzed. This is clearly shown in Figure 12, where the same sample shows a different texture by changing the scanned area.
All the techniques adopted in this work were analysed in such a way to obtain sigmoidal curves representative of the evolution of the crystallization. From all the curves, it was possible to identify an induction time as the time at which the measured signal deviates from the value that it presents at the start of the isothermal step when the sample is still molten, without any crystallinity. A comparison between the results obtained with the different techniques presented is shown in Figure 13.
Results clearly reveal that the induction time found on the entire range of crystallization is dependent on the technique used. As expected, the observation scale plays an important role in its determination. Rheology appears to be the least sensitive technique, probably because in order to measure an increase in rheological parameters, an interaction among crystals is necessary. On the contrary, AFM tests that are able to resolve the surface sample at dimensions comparable to the primary nuclei can detect crystallization at initial stages.
The morphological development taking place during the crystallization of polycaprolactone (PCL) was studied using several techniques.
Results show that the induction time is dependent on the technique used. As expected, the observation scale plays an important role in its determination. In particular, despite its limitations, the AFM that is able to resolve the surface sample at dimensions comparable to the primary nuclei is able to detect crystallization at initial stages. On the contrary, rheology appears to be the least sensitive technique, probably because in order to measure an increase in rheological parameters an interaction among crystals is necessary. DSC and POM techniques supply closer results. It suggests that the definition of induction time in the polymer crystallization is a vague concept that, in any case, requires the definition of the technique used for its characterization.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
 R. Pantani, I. Coccorullo, V. Speranza, and G. Titomanlio, "Modeling of morphology evolution in the injection molding process of thermoplastic polymers," Progress in Polymer Science, vol. 30, no. 12, pp. 1185-1222, 2005.
 R. Pantani, V. Speranza, A. Sorrentino, and G. Titomanlio, "Molecular orientation and strain in injection moulding of thermoplastics," Macromolecular Symposia, vol. 185, no. 1, pp. 293-307, 2002.
 R. Pearce and G. J. Vancso, "Imaging of melting and crystallization of poly (ethylene oxide) in real-time by hot-stage atomic force microscopy," Macromolecules, vol. 30, no. 19, pp. 5843-5848, 1997.
 R. Pearce and G. J. Vancso, "Real-time imaging of melting and crystallization in poly (ethylene oxide) by atomic force microscopy," Polymer, vol. 39, no. 5, pp. 1237-1242, 1998.
 R. Pearce and G. J. Vancso, "Observations of crystallization and melting in poly (ethylene oxide)/poly (methyl methacrylate) blends by hot-stage atomic-force microscopy," Journal of Polymer Science B, vol. 36, no. 14, pp. 2643-2651, 1998.
 J. M. Schultz and M. J. Miles, "AFM study of morphological development during the melt-crystallization of poly (ethylene oxide)," Journal of Polymer Science B, vol. 36, no. 13, pp. 2311-2325, 1998.
 T. J. McMaster, J. K. Hobbs, P. J. Barham, and M. J. Miles, "AFM study of in situ real time polymer crystallization and spherulite structure," Probe Microscopy, vol. 1, no. 1, pp. 43-56, 1997.
 L. Li, C.-M. Chan, K. L. Yeung, J.-X. Li, K.-M. Ng, and Y. Lei, "Direct observation of growth of lamellae and spherulites of a semicrystalline polymer by AFM" Macromolecules, vol. 34, no. 2, pp. 316-325, 2001.
 J. K. Hobbs, T. J. McMaster, M. J. Miles, and P. J. Barham, "Direct observations of the growth of spherulites of poly (hydroxybutyrate-co-valerate) using atomic force microscopy," Polymer, vol. 39, no. 12, pp. 2437-2446, 1998.
 L. G. M. Beekmans, R. Valle, and G. J. Vancso, "Nucleation and growth of poly (e-caprolactone) on poly (tetrafluoroethylene) by in-situ AFM," Macromolecules, vol. 35, no. 25, pp. 9383-9390, 2002.
 W. Zhou, S. Z. D. Cheng, S. Putthanarat et al., "Crystallization, melting and morphology of syndiotactic polypropylene fractions. 4. In situ lamellar single crystal growth and melting in different sectors," Macromolecules, vol. 33, no. 18, pp. 6861-6868, 2000.
 Y. Jiang, D.-D. Yan, X. Gao et al., "Lamellar branching of poly (bisphenol A-co-decane) spherulites at different temperatures studied by high-temperature AFM," Macromolecules, vol. 36, no. 10, pp. 3652-3655, 2003.
 D. A. Ivanov, Z. Amalou, and S. N. Magonov, "Real-time evolution of the lamellar organization of poly (ethylene terephthalate) during crystallization from the melt: high-temperature atomic force microscopy study," Macromolecules, vol. 34, no. 26, pp. 8944-8952, 2001.
 Y. Kikkawa, H. Abe, M. Fujita, T. Iwata, Y. Inoue, and Y. Doi, "Crystal growth in poly (l-lactide) thin film revealed by in situ atomic force microscopy," Macromolecular Chemistry and Physics, vol. 204, no. 15, pp. 1822-1831, 2003.
 R. Pantani, V Speranza, and G. Titomanlio, "Relevance of crystallisation kinetics in the simulation of the injection molding process," International Polymer Processing, vol. 16, no. 1, pp. 61-71, 2001.
 A. Martins, W. Zhang, M. Brito, U. Infante, M. Romero, and F. O. Soares, "Isothermal and nonisothermal crystallization of polymers: analysis with a shear differential thermal analyzer," Review of Scientific Instruments, vol. 76, no. 10, Article ID105105, 2005.
 W. Nagatake, T. Takahashi, Y. Masubuchi, J.-I. Takimoto, and K. Koyama, "Development of shear flow thermal rheometer for direct measurement of crystallization fraction of polymer melts under shear deformation," Polymer, vol. 41, no. 2, pp. 523-531, 2000.
 F. de Santis, A. R. Vietri, and R. Pantani, "Morphology-evolution during polymer crystallization simultaneous calorimetric and optical measurements," Macromolecular Symposia, vol. 234, pp. 7-12, 2006.
 A. Sorrentino, R. Pantani, and G. Titomanlio, "Kinetics of melting and characterization of the thermodynamic and kinetic properties of syndiotactic polystyrene," Journal of Polymer Science B, vol. 45, no. 2, pp. 196-207, 2007
 A. Ziabicki, "Applicability of light transmission technique to studies of crystalline polymers texture," Kolloid-Zeitschrift & Zeitschrift fur Polymere, vol. 219, no. 1, pp. 2-12, 1967
 F. L. Binsbergen, "Depolarization of linearly polarized light by a polycrystalline specimen," Journal of Macromolecular Science B, vol. 4, pp. 837-851, 1970.
 A. Ziabicki and B. Misztal-Faraj, "Applicability of light depolarization technique to crystallization studies," Polymer, vol. 46, no. 7, pp. 2395-2403, 2005.
 Y. P. Khanna, "Rheological mechanism and overview of nucleated crystallization kinetics," Macromolecules, vol. 26, no. 14, pp. 3639-3643, 1993.
 S. Acierno, E. di Maio, S. Iannace, and N. Grizzuti, "Structure development during crystallization of polycaprolactone," Rheologica Acta, vol. 45, no. 4, pp. 387-392, 2006.
 L. Bove and M. R. Nobile, "Shear-induced crystallization of isotactic poly (1-butene)," Macromolecular Symposia, vol. 185, no. 1, pp. 135-147, 2002.
 N. V. Pogodina, V. P. Lavrenko, S. Srinivas, and H. H. Winter, "Rheology and structure of isotactic polypropylene near the gelpoint: quiescent and shear-induced crystallization," Polymer, vol. 42, no. 21, pp. 9031-9043, 2001.
 L. G. M. Beekmans and G. J. Vancso, "Real-time crystallization study of poly ([epsilon]-caprolactone) by hot-stage atomic force microscopy," Polymer, vol. 41, no. 25, pp. 8975-8981, 2000.
V. Speranza, (1) A. Sorrentino, (2) F. De Santis, (1) and R. Pantani (1)
(1) Department of Industrial Engineering (DIIN), University of Salerno, Via Giovanni Paolo II, 84084 Fisciano, Italy
(2) Institute for Composite and Biomedical Materials (IMCB), National Research Council (CNR), Piazzale Enrico Fermi 1, 80055 Portici, Italy
Correspondence should be addressed to A. Sorrentino; firstname.lastname@example.org
Received 28 August 2013; Accepted 8 October 2013; Published 9 January 2014
Academic Editors: V. Amigo, P. Reis, and S. Wu
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Speranza, V.; Sorrentino, A.; De Santis, F.; Pantani, R.|
|Publication:||The Scientific World Journal|
|Date:||Jan 1, 2014|
|Previous Article:||Gaining surgical access for repositioning the inferior alveolar neurovascular bundle.|
|Next Article:||Effect of catalytic pyrolysis conditions using pulse current heating method on pyrolysis products of wood biomass.|