Abstract
In this work, the effects of cellulose nanofibrils (CNFs) on workability, hydration reaction, microstructure, early age shrinkage, fracture properties, flexural strength, and compressive strength of cement paste were investigated. Six batches with variable CNF concentrations with the same water-to-cement (w/c) ratio (0.35) were tested. Flow table test showed a reduction in the workability as CNF dosage increased. Isothermal calorimetry (IC) tests showed that after 3 days, degree of hydration (DOH) improved up to 8% because of the addition of CNFs. Thermogravimetric analysis (TGA) tests at 7 and 28 days showed no significant changes in DOH for all pastes. After 7 days, mixture with 0.15% CNF resulted in up to 31% improvement in compressive strength. For 0.09% CNF addition, cement paste showed 26% increase in compressive strength after 28 days. Tests revealed that adding a small quantity of CNF (0.06%) along with entraining 0.05 extra water reduces autogenous shrinkage by 49% at a cement paste with w/c = 0.30. For interpreting the results, a tunnels, reservoirs, and bridges (TR&B) model is proposed. This model suggests that, as proposed by others, CNFs can modify microstructure by providing tunnels for transporting water to unhydrated cement grain. Because of their hydrophilicity, CNFs retain water and work as reservoirs (internal curing), which explains the improvement in properties at low w/c ratios. Significant increases in fracture energy (up to 60%) and flexural strength (up to 116%) suggest that CNFs are an effective toughening mechanism, acting as bridges that increase the energy required for crack propagation.
Nanotechnology has shown, and continues to show, tremendous promise for enhancing the short-term and long-term properties of cement-based composites. While investigations of nanotechnology applications to cement and concrete have been proceeding for some time, recent successes illustrate some profound effects that nanomaterials can have on concrete properties ( 1 – 4 ).
A class of nanomaterial that holds promise in concrete technology is that based on cellulose. Plant-based fibers have long played a role in concrete for their role in reinforcing as well as their role in early age properties ( 5 – 8 ). Cellulosic materials have the advantages of being abundant, low toxicity, and relatively low cost.
Cellulosic nanomaterials have gradually been making their way into concrete technology ( 9 ). The effect of cellulose microcrystalline particles (MCC) on the properties of cement-based composites was studied by Gomez Hoyos et al., who found adding MCC increases the degree of hydration (DOH) ( 10 ). Onuaguluchi et al. worked on the effects of cellulose nanofiber gel suspension on the properties of cement pastes, showing that DOH, flexural strength, and energy absorption improved in pastes with cellulose nanofiber ( 11 ). Cao et al. found that the addition of cellulose nanocrystals (CNCs) enhances the performance of cement paste ( 12 ). They showed that CNCs can improve DOH and the flexural properties of cement pastes. They further showed that adding CNC can increase the flexural strength of cement paste by up to 50% ( 13 ). Other researchers worked on the effects of CNCs on the microstructure of cement paste ( 14 , 15 ). Fu et al. used nine various CNCs for improving the hydration and flexural properties of Portland cement pastes ( 16 ). Hisseine et al. performed some investigations on cellulose filaments in cement-based composites. They found that adding cellulose filaments enhances the mechanical performance of cement systems and improves DOH ( 17 – 19 ).
In the work described here, attention is focused on the potential role of cellulose nanofibrils (CNFs). CNFs are hydrophilic nanoparticles (typically, 20–50 nm in width and less than 0.2 mm in length). They can be extracted from various plants, trees, and recycled paper stock. CNFs are usually branched or forked. They are new nanoscale material given their high specific surface (31–33 m2/g), high aspect ratio, low density (1.0 g/cc slurry), low toxicity, and low cost ($1.25/lb. = $2.67/kg slurry) that enable functionalization ( 20 ). Previous work showed that the addition of small quantities of CNFs and micro cellulose fibers improved fracture toughness of high-performance concrete by up to 50% ( 21 ). Haddad Kolour et al. showed that adding a small quantity of CNF can reduce free shrinkage by 13% and improve compressive strength by up to 28% ( 22 , 23 ).
Despite these preliminary successes, important questions remain about the specific role that CNFs play in cement hydration and the development of microstructure. Therefore, the objective of this work was to examine the changes in hydration brought on by small dosages of CNF to batches of cement paste, and to tie those changes to observed early age shrinkage and strength properties. To realize this objective, an array of laboratory techniques for measuring both reaction kinetics and resulting properties were employed. In this work, the focus was on low water-to-cement (w/c) ratio systems, as preliminary work indicated that these systems realize the greatest benefit from CNF additions ( 22 ).
The CNF materials used in this study were the product of the University of Maine Process Development Center. CNFs were extracted mechanically (not chemically/thermally) from northern bleached softwood Kraft pulp (spruce/fir). As-received CNF materials were in a white slurry form. The concentration of solids was 3.0% (3 g of nanofibrils and 97 g of water in 100 g of suspension), and the density of aqueous gel was 1.0 g/cm3 ( 20 ). Because of their high surface area and surface hydroxyl (OH–) groups, CNFs are hydrophilic and ready for surface reactivity and interactions ( 10 ).
Methods
Specimen Preparation
The paste specimens were used to study workability, hydration kinetics, early age shrinkage, compressive strength, flexural strength, and fracture properties ( 23 ). Ordinary Portland cement (OPC) Type I/II (commercial grade) that complies with ASTM C150/C150M-17, Standard Specification for Portland Cement, was used in this study ( 24 ). The OPC contained 20.1% SiO2, 62.0% CaO, 3.6% Al2O3, 3.0% Fe2O3, 3.4% SO3, 3.3% MgO, 0.4% Na2O, and 1.0% K2O with a Blaine fineness of 378 m2/kg. The cement pastes were mixed with a conventional 8-qt rotary kitchen mixer. The following procedure was used for mixing the pastes: (1) the CNF suspension was mixed with water in the kitchen mixer for 180 s at high speed (homogenization/dispersion); (2) the solution from step 1 was combined with cement powder and mixed at low speed for 60 s; (3) the mix was allowed to rest for 15 s; a spatula was used to scrape the wall and bottom of the mixing bowl; (4) mixing for 240 s at high speed. Sufficient periods of mixing were used to improve CNF dispersion ( 6 ).
The CNF concentrations of each batch were calculated based on weight with respect to cement. Cement pastes were prepared at six different CNF concentrations. Test matrix for CNF-reinforced cement paste is shown in Table 1. Water demand for each batch was modified by subtracting the amount of water brought by the CNF-water suspension (3% solid nanofibrils and 97% water).
Test Matrix for Cellulose Nanofibril (CNF)-Reinforced Cement Pastes
An experimental program was developed to measure the effects of CNF dosages on the following properties: workability, hydration kinetics, early age shrinkage, fracture properties, flexural strength, and compressive strength. The specific laboratory tests are detailed below.
Workability
For each batch, and immediately after finishing the mixing procedure, a flow table test following ASTM C1437-15, Standard Test Method for Flow of Hydraulic Cement Mortar, was conducted ( 25 ). The result in percentage is the flow of that batch.
Zeta Potential
In colloidal chemistry, the zeta potential (ζ) is the potential between the fixed layer of liquid adjacent to the solid phase and the liquid constituting the bulk liquid phase ( 26 ). This test quantifies the strength of the electrostatic attraction or repulsion between particles. For studying the affinity between particles in the fresh cement paste, the zeta potentials of the cement and CNF particles were measured ( 12 ). The CNF and cement particles were, respectively, diluted in deionized (DI) water or synthetic pore solution to a concentration of about 100 mg/l ( 27 ). Synthetic pore solution (0.35 M KOH + 0.05 M NaOH in DI water) was the same solution suggested by Rajabipour et al. ( 28 ). Since the zeta potential is dependent on pH values, the tests were performed in two different pH environments at the age of 1 h: a neutral environment with a pH of 7, and the simulated pore solution with a pH of 12.7 ( 12 , 27 ).
Cement Hydration
Isothermal Calorimetry
Heat flow rate and cumulative heat release were monitored in this test in an attempt to study the hydration reaction and the DOH of the cement pastes. Immediately after mixing, 19–24 g samples of the paste were transferred to a 20 ml glass ampoule, 1 in. (25.4 mm) in diameter and 2¼ in. (57.15 mm) in height. Measurement has been performed in a sealed isothermal calorimetry (IC) chamber (maintained at 23 ± 0.1°C). Before starting the data collection, the isothermal condition was held for 45 min to reach equilibration and then steady heat measurement was run for almost 3 days. Two replicate tests were performed for each of the six different CNF dosages.
Thermogravimetric Analysis
The thermogravimetric analysis (TGA) was performed to find the DOH of CNF-reinforced cement pastes at two different ages: 7 and 28 days. At later ages, the heat release rate from IC is so small and the measuring error could be too large. So TGA was run to obtain the DOH at larger ages ( 12 ). The TGA tests were also done for the individual materials of CNF and cement. After mixing the paste, two 1 by 1 in. (25.4 by 25.4 mm) plastic cylinder molds were used to prepare two specimens for each batch (one specimen for each age). Molds were covered by plastic films and kept inside zipped bags (sealed condition) for 24 h. Then, specimens were demolded, covered by plastic films, and cured inside other zipped bags (sealed condition) for 7 and 28 days. Two replicate tests were performed for each of the six different CNF dosages.
At the ages of testing, the paste specimens were removed from the zipped bags. Immediately after removal, hydration of the cement paste was stopped by solvent exchange technique. Cement specimens were submerged in a relatively large amount of solvent for 30–40 min. Water present in the pores of the cement paste is diluted and removed by solvent. In this research, ethanol was used as solvent ( 29 ).
After stopping the hydration, specimens were stored in a sealed condition, ready for TGA. At the day of testing, they were ground into powders with mortar and pestle while evaporation was minimum. Approximately 35 mg of powder was transferred into the TGA chamber for running the test. First, the specimen was held for 5 min to reach equilibration temperature in the chamber. At the second step, the specimen was heated to 1,000°C at a rate of 10°C/min. This step is done to remove all chemically bound water (CBW).
In TGA tests, the total mass of the chemically bound water (CBW) in the hardened cement pastes is calculated. CBW is the weight loss between 105°C and 1,000°C ( 30 ). At temperatures below 105°C, some relatively small amounts of CBW could be missed from some hydration products (particularly C-S-H, ettringite, and monosulfate) ( 29 ). As a simplifying assumption, this part of CBW has been ignored in this study. The decrease in mass observed between 400°C and 500°C shows the decomposition of Ca(OH)2 ( 29 ). Almost all of the weight loss above 600°C is related to CO2 loss because of different forms of CaCO3 decomposition (decarbonation). This part of weight loss, because of CO2 evaporation, was not considered as CBW.
CBW divided by the final weight of the material shows the mass of CBW per unit gram of unhydrated cement. The assumption for the maximum amount of CBW per unit gram of fully hydrated cement was 0.23 g ( 12 ). Then, DOH is the mass of CBW per unit gram of unhydrated cement divided by 0.23 g.
Compressive Strength
After mixing the paste, six 2 by 4 in. (50.8 by 101.6 mm) cylinder molds were used to mold specimens for each batch to obtain the compressive strength at two different ages: 7 and 28 days (three specimens for each age). Molds were covered by plastic films and kept inside zipped bags (sealed condition) for 24 h. Then, specimens were demolded, covered by plastic films and cured inside other zipped bags (sealed condition) for 7 and 28 days. ASTM C39/C39M-18 (AASHTO T22), Standard Test Method for Compressive Strength of Cylindrical Concrete Specimens, was used for breaking the specimens and measuring the compressive strengths ( 31 ). Results are the average value of the compressive strengths of three specimens for each batch at that particular age.
Early Age Shrinkage
The positive effects of CNF on both hydration reaction (3 days) and compressive strength showed that the CNF is acting as an internal curing agent in a way similar to a super-absorbent polymer (SAP) ( 32 – 35 ). That is, some of the mix water is being captured and held by the CNF during mixing but is released later on as hydration progresses. Such a mechanism—water entraining—would explain both the DOH and strength improvements, as it would play a role in minimizing self-desiccation and autogenous shrinkage.
Autogenous deformation was measured using a dilatometer equipped with an automatic data-logger and a digital length gauge. The test was performed using the corrugated tube based on ASTM C1698-09, Standard Test Method for Autogenous Strain of Cement Paste and Mortar ( 36 ). The length change of corrugated polyethylene molds, approximately 29 mm in diameter and 420 mm in length, were monitored in this test. Corrugated molds should be sealed with two end plugs and filled with fresh paste. A tamping rod was used for consolidating the fresh paste along the length of the tube. The corrugated shape of the tubes makes them much stiffer in the radial direction than in the longitudinal direction. This higher stiffness in radial direction translates almost all of the volumetric deformation of the paste to linear deformation of the mold ( 6 ).
For examining the hypothesis of water entraining by CNF, pastes with and without CNF were compared in a series of autogenous shrinkage tests. The basic w/c ratio is 0.30 in all mixes. Mixes with CNF, additionally contain entrained water. An amount of entrained water equal to (w/c)e = 0.05 was added to mixes with CNF in an attempt to decrease self-desiccation and autogenous shrinkage based on a method proposed by Jensen et al. ( 32 , 33 ). The changes in lengths of the two samples were recorded for each mix. Each test was conducted for 7 days in a thermostatically controlled room set at 23°C. Every 30 min, for the whole duration of the test, the lengths of the samples were continuously monitored. The time of the final set is zero point for all of the readings. Before that, the major deformation is deformation of the polyethylene corrugated tube. Also, it is known that before the final set, the paste is still plastic and the reading is actually related to chemical shrinkage ( 6 ). The time of final set for each mix was determined by the Vicat needle test ( 37 ).
Notched Beam
Notched beam tests were done after 28 days for all cement pastes (four specimens for each batch in Table 1), using a simplified hybrid fracture test that combines the measurements required for the two-parameter fracture model of Jenq and Shah and the measurements required for concrete fracture energy (RILEM 1985) ( 38 ). Molds were covered by plastic films and kept inside zipped bags (sealed condition) for 24 h. Then, the specimens were demolded, covered by plastic films and paper towel, and cured inside other zipped bags (sealed condition) for 28 days. In this simplified method, midspan displacement was not measured, so Gf was redefined as the area under the complete load-crack mouth opening displacement (CMOD) curve. By this method, it was possible to measure the modulus of elasticity (E in GPa) and the material fracture parameters, critical crack length (ac in mm), and fracture energy Gf (area under Load-CMOD curve in N.mm). The span, depth, and width of the beams were 180 mm, 40 mm, and 30 mm, respectively. The beams had an initial notch of 13 mm (1/3 of depth), and the notch was created by a diamond wet saw. The beams were tested in an Instron closed-loop testing machine. The three-point bending test of notched beams was done in CMOD control mode, and the rate of loading was 0.015 mm/min. A clip gage, mounted on the notch, controls the CMOD opening. The test consisted of loading the beam until a crack happens. Then the specimen was unloaded and reloaded so that the crack grew steadily ( 21 , 38 ).
Flexural Strength (Three-Point Bending)
Flexural strength tests were done after 28 days for all cement pastes (four specimens for each batch in Table 1). The beam sizes and curing procedure were the same as described above for the notched beam tests. By this test, it was possible to measure the flexural strength (modulus of rupture) of all specimens. Three-point bending tests were performed in an Instron closed-loop testing machine. The test consisted of loading the beam until breaking the specimen. The test was displacement (deflection) control and the rate of loading was 0.5 mm/min.
Results
Workability
Table 2 shows the effect of CNF on the flow table results for different cement batches with various amounts of CNF. The results show that increasing CNF content leads to reduced workability. It can be seen that every 0.05% of additional CNF incorporated in the cement paste decreases the flow value by almost 15%. Results from the flow table tests suggest that water in CNF slurry is not immediately available to mix with the cement, as it is being held back by the CNF. Indeed, the loss of workability suggests that some of the added mix water is held back by the CNF. Possible CNFs agglomeration because of poor dispersion (especially at high dosages) likely also plays an important role in decreasing workability.
Flow Table Results (Water/Cement = 35%)
Zeta Potential
The zeta potentials for the CNF and cement particles at these two different pH environments are shown in Table 3. Results show that the absolute value of the zeta potential for CNF fibers is much higher than that for cement particles. This means that cement particles have much more potential for agglomeration in comparison with CNF fibers. At pH = 12.7, synthetic pore solution, the affinities between the particles have following order:
The Zeta Potentials of Cellulose Nanofibril (CNF) and Cement Particles
Note: DI = deionized.
The zeta potential results show that the affinity between cement and CNF is stronger than that between cement particles. It means that CNF fibers tend to adhere onto cement particles rather than to agglomerate themselves. For CNFs-cement adhesion, well-dispersion of both cement particles and CNF fibers is a necessity. Relatively low zeta potentials of cement (+0.4 mV) and CNF (–9.5 mV) at synthetic pore solution show that there is a potential of agglomeration within the fresh cement paste for both CNFs and cement particles.
Cement Hydration
Isothermal Calorimetry
Results of the cumulative heat of hydration for the first 3 days (72 h) of the six cement pastes with different CNF dosages are shown in Figure 1. It can be observed that the cumulative heat showed increase in all CNF-modified pastes. This increment continues until the end of the test (3 days). All CNF dosages showed higher DOH than the plain paste. After 3 days, the CNF dosages of 0.09% and 0.15% resulted in 8% improvement in DOH, while for the pastes with 0.015%, 0.03%, and 0.06% CNF, changes in DOH are not significant. The increased DOH could result from internal curing effects and short circuit diffusion because of the CNF presence ( 12 ). Figure 2 displays hydration reaction heat flow curves for the first 40 h of the six mixtures described in Table 1.

Cumulative heat of cellulose nanofibril (CNF)-reinforced cement pastes for the first 3 days.

Heat flow curves of the cellulose nanofibril (CNF)-reinforced cement pastes for the first 40 h.
Thermogravimetric Analysis
Following calculations described in the TGA methods section, after 7 days, the DOH for the control specimen (without CNF) was 61. For cement pastes with 0.015%, 0.03%, 0.06%, 0.09%, and 0.15% of CNF, DOHs were 60, 60, 60, 61, and 60, respectively. At 28 days, the DOH for control specimen was 64. For cement pastes with 0.015%, 0.03%, 0.06%, 0.09%, and 0.15% of CNF, DOHs were 64, 64, 64, 64, and 67, respectively. These results show that incorporating CNF does not change DOH significantly at the ages of 7 and 28 days. Figure 3 shows TGA results at the age of 28 days.

28-day thermogravimetric analysis (TGA) results.
Compressive Strength
Results of compressive strength tests for cement pastes in Table 1 can be seen in Figure 4 for the age 7 and 28 days. After 7 days, the CNF dosages of 0.015%, 0.03%, 0.06%, 0.09%, and 0.15% resulted 20%, 3%, 18%, 28%, and 31% improvement in compressive strength, respectively. For the same CNF contents and after 28 days, compressive strengths are 2%, 23%, 12%, 26%, and 16% higher in comparison with the plain system.

Cement pastes compressive strength results.
Early Age Shrinkage
Figure 5 shows the results of the autogenous deformation of cement paste specimens with and without CNF. It can be observed that adding 0.06% and 0.09% of CNF reduces autogenous shrinkage up to 49% and 26%, respectively. Reduction in internal curing effects of CNF at 0.09% could be because of potential CNFs agglomeration at higher dosage.

Autogenous shrinkage versus time for pastes.
It is interesting to consider the expansion right after set in the mixtures with CNF. It happens because of water released by CNF, because of a drop in internal relative humidity inside the sealed (closed) cement paste. This phenomenon has already been observed in the cement system using SAP as an internal curing agent ( 33 ).
Compressive strengths of these mixtures have been depicted in Table 4. It shows that after 7 days, the CNF dosages of 0.06% and 0.09% resulted in 5% and 15% improvement in compressive strength, respectively. After 28 days, results for 0.06% addition of CNF are almost the same as for the plain system. At the same time, for the specimen with 0.09% CNF, compressive strength is 12% higher.
Compressive Strength Results for Autogenous Shrinkage Mixtures
Notched Beam
Table 5 shows the measured modulus of elasticity (E), the critical crack length (ac), and fracture energy (Gf). Small changes in critical crack length show that CNFs are not very important in preventing crack initiation because they may not always be available at the crack tip. After 28 days, the CNF dosages of 0.015%, 0.03%, 0.06%, 0.09%, and 0.15% resulted in 0%, 17%, 11%, 24%, and 60% improvement in fracture energy, respectively. Significant changes in fracture energy suggest that CNFs are effective in preventing crack growth and the CNFs play an important role when the cracks try to propagate within the notched beam specimen ( 21 ).
Notched Beam Tests Results
Flexural Strength (Three-Point Bending)
Results for pastes in Table 1 can be seen in Table 6. After 28 days, the CNF dosages of 0.015%, 0.03%, 0.06%, 0.09%, and 0.15% showed 44%, 39%, 55%, 72%, and 116% improvement, respectively. Significant changes in flexural strength suggest that CNFs are effective in improving the modulus of rupture in cement pastes.
Flexural Strength Tests Results
Discussion
Because of their high surface area and surface hydroxyl (OH–) groups, CNFs are ready for surface reactivity and interactions ( 10 ). Water which contains hydrogen have the potential to make hydrogen bonds with surface hydroxyl (OH–) groups of CNFs. Because of these potential hydrogen bonds, CNFs are hydrophilic and can capture some water (depending on CNF dispersion/agglomeration condition) at early ages and release it at later ages.
The improvement in 7- and 28-day compressive strength observed in all mixes can be explained once again by an effect of internal curing, which suggests CNFs retain water at early ages and release it later. Tests revealed that, like SAPs, incorporating CNFs helps to entrain extra water and reduces autogenous shrinkage. At the same time, compressive strength was enhanced in those specimens. However, looking at Figure 4 shows that this improvement is not systematic with respect to CNF dosage. The possible interpretation of these tests is that flaw sizes in the cement pastes could be increased by the potential CNF agglomerating. At low dosages, the influence of this potential agglomeration is small, but it likely plays a more important role at high dosages.
The CNF system appear to work with a mechanism that is here named as the “tunnels, reservoirs, and bridges” (TR&B) mechanism ( 23 ). The tunnels mechanism described here is based on Cao et al.’s hypothesis ( 12 ). Zeta potential results showed that CNFs tend to adhere onto cement particles. Porous space around these adhered CNFs will work as “tunnels.” It means that, similar to CNCs, CNFs provide tunnels (channels) for transporting water to unhydrated cement grain which in turn help to produce more hydration products and enhance total hydration reaction ( 12 ).
CNFs are hydrophilic and can retain some water (depending on the CNF dispersion/agglomeration situation) and, like SAPs, work as “reservoirs.” This stored water is released at later ages and works as a supplementary source of water. Effect of this supplementary water can be seen in reducing autogenous shrinkage in Figure 5.
Significant changes (60% increase) in fracture energy (Table 5) suggest that CNFs are effective in preventing crack growth, and the CNFs work as “bridges” and play an important role when the cracks try to propagate within the notched beam specimen. Figure 6 shows the schematic of the TR&B model based on the above discussions.

Schematic of tunnels, reservoirs, and bridges (TR&B) model.
Conclusion
This project examined how the addition of CNFs modifies the performance of cement paste and mortar. Flow table tests on cement pastes showed that increasing CNF dosages decrease workability. Based on these results, it is possible to work on using CNFs as viscosity modifying agents, with or without other admixtures (e.g., in self-compacting concretes). IC tests disclosed that incorporating CNF improved DOH after 3 days. TGA tests at the ages of 7 and 28 days showed no significant changes in DOH because of adding CNF. After 7 and 28 days, CNF-reinforced specimens showed higher compressive strength. A new TR&B model has been proposed for interpreting the results.
Based on CNF dosage, and dispersion efficiency, the effects of CNF on the cement-based systems can vary. High dosages of CNF increase the likelihood of agglomeration and covering the surface of cement particles and, as a result, slowing down the hydration reaction at early ages ( 12 ). Good dispersion optimizes the CNF dosage, while poor dispersion results agglomeration and dysfunction. Agglomerated CNFs work as flaws in the system. At low dosages, this effect is negligible, but at high dosages these agglomerations (flaws) can deteriorate the microstructure and mechanical properties of cement paste. CNF source, cement type, mixing, and curing procedures are the other important factors that can affect the test results.
Understanding the status of water during mixing and hydration allows better control of both the effects on workability and the effects on microstructure. As a future research, studying the water kinetics of CNFs during the mixing procedure and cement hydration is highly recommended.
Footnotes
Acknowledgements
The authors gratefully acknowledge the Public-Private Partnership for Nanotechnology (P3Nano) and the USDA Forest Service for their support of this work. Also, the authors would like to acknowledge the contributions of both the University of Maine Process Development Center, who provided the CNF for this study, and Parivash Takasi, who assisted with laboratory work.
Author Contributions
The authors confirm contribution to the paper as follows: study conception and design: H. Haddad Kolour, W. Ashraf, E. N. Landis; data collection: H. Haddad Kolour; analysis and interpretation of results: H. Haddad Kolour, W. Ashraf, E. N. Landis; draft manuscript preparation: H. Haddad Kolour, W. Ashraf, E. N. Landis. All authors reviewed the results and approved the final version of the manuscript.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
