# Influence of Rice Husk and Wood Biomass Properties on the Manufacture of Filaments for Fused Deposition Modeling

Marie-Joo Le Guen, Stefan Hill, Dawn Smith, Beatrix Theobald, Evamaria Gaugler, Abdellatif Barakat, Claire Mayer-Laigle

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Marie-Joo Le Guen, Stefan Hill, Dawn Smith, Beatrix Theobald, Evamaria Gaugler, et al. Influence of Rice Husk and Wood Biomass Properties on the Manufacture of Filaments for Fused Deposition Modeling. Frontiers in Chemistry, 2019, 7, 10.3389/fchem.2019.00735. [https://doi.org/10.3389/fchem.2019.00735](https://doi.org/10.3389/fchem.2019.00735)

HAL Id: hal-02342165

https://hal.science/hal-02342165

Submitted on 5 Nov 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

# Influence of Rice Husk and Wood Biomass Properties on the Manufacture of Filaments for Fused Deposition Modeling

Marie-Joo Le Guen¹*, Stefan Hill¹, Dawn Smith¹, Beatrix Theobald¹, Evamaria Gaugler¹, Abdellatif Barakat² and Claire Mayer-Laigle²

¹Scion, Rotorua, New Zealand, ²IATE, Univ Montpellier, CIRAD, INRA, Montpellier SupAgro, Montpellier, France

Additive manufacturing or 3D printing has the potential to displace some of the current manufacturing techniques and is particularly attractive if local renewable waste resources can be used. In this study, rice husk, and wood powders were compounded in polylactic acid (PLA) by twin screw extrusion to produce filaments for fused-deposition modeling 3D printing. The biomasses were characterized in terms of physical features (e.g., particle size, density) and chemical compositions (e.g., solid state nuclear magnetic resonance, ash content). The two biomasses were found to have a different impact on the rheological behavior of the compounds and the extrusion process overall stability. When comparing the complex viscosity of neat PLA to the biomass/PLA compounds, the integration of wood powder increased the complex viscosity of the compound, whereas the integration of rice husk powder decreased it. This significant difference in rheological behavior was attributed to the higher specific surface area (and chemical reactivity) of the rice husk particles and the presence of silica in rice husks compared to the wood powder. Color variations were also observed. Despite the biomass filler and rheological behavior differences, the mechanical properties of the 3D printed samples were similar and predominantly affected by the printing direction.

**Keywords:** 3D printing, rice husk, wood, biomass, composites, extrusion

## INTRODUCTION

The popularity of incorporating lignocellulosic biomasses in thermoplastic polymers has increased significantly due to their sustainability advantage, low density, thermal, and noise insulation capacity and stiffening performances (Le Guen and Newman, 2007; Müssig, 2010; Gallos et al., 2017; Lammi et al., 2018). The properties of the composites are related to (a) the amount of biomass incorporated and (b) the interactions between the polymeric matrix and the lignocellulosic particles in terms of physical (e.g., surface area, aspect ratio) and chemical interactions (e.g., interface compatibilisation). Amongst the lignocellulosic feedstocks available, rice husk is of particular interest due to its current abundance (i.e., estimated between 128 and 148 million tons of unutilized waste) (Giddel and Jivan, 2007; Pode, 2016), and chemical make-up. It is composed of around 80–85 wt.% lignocellulosic materials (e.g., cellulose, hemicellulose, lignin), 15–20 wt.% of amorphous silica and 1.1 to 2.5 wt.% of proteins (Juliano et al., 1987; Vadiveloo et al., 2009; Ummah et al., 2015).

## OPEN ACCESS

**Edited by:**
Caroline Remond,
Université de Reims
Champagne-Ardenne, France

**Reviewed by:**
Antti Ilmari Koponen,
VTT Technical Research Centre of
Finland Ltd., Finland
Rupinder Singh,
Guru Nanak Dev Engineering
College, India

*Correspondence:
Marie-Joo Le Guen
mariejoo.leguen@scionresearch.com

**Specialty section:**
This article was submitted to Chemical and Process Engineering, a section of the journal Frontiers in Chemistry

**Received:** 15 June 2019
**Accepted:** 14 October 2019
**Published:** 31 October 2019

**Citation:** Le Guen M-J, Hill S, Smith D, Theobald B, Gaugler E, Barakat A and Mayer-Laigle C (2019) Influence of Rice Husk and Wood Biomass Properties on the Manufacture of Filaments for Fused Deposition Modeling. Front. Chem. 7:735. doi: 10.3389/fchem.2019.00735
Previous research studies have reported its potential applications in composite structures (Kumar et al., 2012), as an energy feedstock (Chungsangunsit, 2009) or amorphous silica feedstock (Yalçin and Sevinç, 2001). In parallel, the incorporation of nano-silica in polylactic acid (PLA) was also reported to improve the thermal stability and mechanical properties of composites (Lv et al., 2016).

3D printing (3DP) or additive manufacturing is an emerging technology that enables the creation of innovative designs and the combination of materials that were previously impossible or impractical to make. The incorporation of fillers in 3DP has demonstrated that the filler's physical features could be transferred to the 3D printed object (Montalvo Navarrete et al., 2018; Zeidler et al., 2018). These attractive traits have been demonstrated with wood (Mirko et al., 2016; Tao et al., 2017), clays (Revelo and Colorado, 2018), metals (Gibson et al., 2018), and seashells (Graichen et al., 2017; Singamneni et al., 2018a; Zeidler et al., 2018).

The combination of (a) the inherent properties of polymer-filled compounds and (b) the control of the 3D printing process directionality, leads to new functionalities to printed objects (Liles et al., 2018; Wang et al., 2018). Examples of such functionalities for lignocellulosic feedstocks include hydromorphic properties (Le Duigou et al., 2016; Sydney Gladman et al., 2016), electrical conductivity (Shao et al., 2018) and vibration damping (Zeidler et al., 2018).

In this present research, the influence of rice husk powder integrated into a PLA matrix is investigated for fused deposition modeling (FDM) 3DP applications, and compared to wood flour PLA 3D prints. The impact of the two biomasses during the extrusion compounding process, and on the mechanical properties of 3D printed rectangular parallelepipeds are observed and quantified.

# MATERIALS AND METHODS

## Biomasses and Polymer

Carmargue rice (Oryza) husks were supplied from Silo des Tourelles- Comptoir Agricole du Languedoc (Aigues Mortes, France). The husks were first sieved with a screen of 0.8 mm using a vibratory sieve shaker 400 (Ritec, France) to remove soil dust and stones. The husks were subsequently milled in a vibratory ball mill model DM-1 (Sweco®, USA) consisting of a 36 L abrasion-resistant elastomer grinding chamber filled with 25 kg of ceramics cylinders (diameter and length of 13.5 mm) and 25 kg of ceramics balls (diameter of 13.5 mm). The chamber's motion was controlled by a vibrating mechanism composed of high-tensile steel springs. One kilogram of rice husk was placed in the chamber and milled for 10 h to obtain fine rice husk powder. The ground powder was dried in an oven at 60°C for 12 h and stored in aluminum foil vacuum packaging bags prior to utilization.

Wood flour (Pinus radiata D. Don) was supplied by a plywood mill (Carter Holt Harvey, Kinleith Mill, Tokoroa, New Zealand). Only the particles below 125 µm were used after separation using a 43305 F motorized sieve shaker (Humboldt, USA) and 125 µm

mesh sieve (Endercotts, UK). The wood flour powder was dried at 100°C for at least 24 h prior to extrusion.

PLA 3001D was obtained from Nature works (USA) and dried overnight prior to extrusion. The PLA and PLA compounds were dried for at least 4 h at 55°C prior to extrusion.

## Processing

### Compounding

The filaments were made in two sequential steps: (a) compounding the polymer and biomasses at 10 wt.% into pellets and (b) extrusion of the filament for FDM. Both steps were carried out on a 26-mm scientific twin screw, co-rotating extruder LTE26-40 (Labtech Engineering Ltd, Thailand), which had a 40 L/D ratio and was fitted with a vacuum crammer. During Step (a), materials were hand-mixed in plastic bags and fed into the volumetric feeder. The speed of the feeder was set to 20 rpm. The extruder screw speed was set to 200 rpm. The temperature was set to 190°C at the first 5 zones and 200°C at the last 5 zones. Pellets were cut to around 2 mm in length. During Step (b), the settings of the twin-screw extruder were modified (Table 1) and the pelletizer replaced by a winding unit set to produce a 1.75 mm diameter filament.

### 3D Printing

The objects were printed on a M2 fused deposition modeling printer (Makergear, USA) fitted with a 0.75 mm die nozzle. The platform temperature was set at 70°C and the nozzle temperature at 210°C. Rectangular beams (60 × 10 × 2.5 mm) for flexural testing were printed using slicing software (Simplify 3D®, USA) based on 0.2 mm layers thickness and no strand overlap. The 3D printing direction was carried out at either 0° or 90° of the beam length.

## Characterization of the Biomasses

### Physical Characterization

#### Particle size distribution

The particle size distribution of the risk husk powder and wood flour powder were measured by a laser diffraction particle size analyzer Hydro 2000S (Malvern Instruments Ltd., UK) with distilled water/ethylic alcohol solution (50:50 in volume) to allow good powder dispersion whilst avoiding swelling. The results are expressed using the Mie method with a refraction index of 1.53 corresponding to those of sawdust particles (Malvern, 2007). The results are generated in volume, based on the assumption that the particles are spherical. Five replicates were carried out for each biomass and the mean average distribution is reported.

<table><caption>TABLE 1 | Extrusion settings for Step (b) to manufacture the FDM filaments.</caption><tr><td>Temperature (°C)</td><td>160, 170, 170, 170, 180, 180, 190, 190, 200, 200</td></tr><tr><td>Screw speed (rpm)</td><td>200</td></tr><tr><td>Feed speed (rpm)</td><td>8</td></tr><tr><td>Speed of haul-off unit (rpm)</td><td>8-13</td></tr></table>
### Scanning electron microscopy (SEM)

Rice husk and wood powders were imaged with a scanning electron microscope Phenom Pro-X (Phenom, France) using a charge reduction sample holder and an acceleration voltage of 5 kV (wood powder) and 10 kV (rice husk powder), respectively.

### Tapped density and apparent density

The apparent and tapped density was measured in duplicate using a standard tapped density analyzer (Autotap, Quantachrome instrument®, USA). The experiment consisted of filling a graduated test tube with 200 mL of powder and submitting it to 1,500 taps. The apparent and tapped densities were calculated as the ratio of the mass of the powder to the initial and final volume respectively.

## Chemical Characterization

### pH measurements

One gram of biomass powder was stirred into 15 mL of distilled water at room temperature. Four pH measurements over 4 h were recorded using a ph meter (Mettler Toledo, USA). The four measurements were then averaged.

### Solid-state nuclear magnetic resonance (NMR) spectroscopy

Solid state CP-MAS (Cross Polarization-Magic Angle Spinning) ¹³C Nuclear Magnetic Resonance spectra were acquired at 50.3 MHz on a Bruker Avance III 200 spectrometer (Bruker BioSpin, Germany) fitted with a 4 mm solid state MAS probe (Bruker BioSpin, Switzerland). The samples were packed into a 4 mm rotor and spun at 5 kHz. Each 3 µs, a 90° proton preparation pulse was followed by a 1 ms contact time and a 30 ms acquisition time, with proton decoupling, and a 1.5 s recycle delay. All spectra were calibrated so that the cellulose Iβ interior C4 peak was assigned a value of 89.3 ppm relative to tetramethylsilane (Newman et al., 1993).

### Ash content

Around 2 g of powder was first dried in an oven at 130°C during 90 min and weighed to determine the dry mass (m<sub>s</sub>). The sample was subsequently heated at 900°C for 2 h and re-weighted (m<sub>r</sub>) after cooling in a desiccator to room temperature. The ash content was calculated according to:

$$ Ash = 100 \cdot m_r/m_s \quad (1) $$

Where:

Ash is in percentage,

m<sub>s</sub> is the sample dry mass in grams,

m<sub>r</sub> is the sample residual mass in grams.

Experiments were carried out in triplicate and averaged.

### Chemical composition

Carbohydrates and lignin content of the biomasses were measured after concentrated acid hydrolysis. The lignin content in samples was determined by the Klason method. Briefly, 100 mg of dried samples were treated with 72% v/v H₂SO₄ at ambient temperature for 2 h. The solutions were diluted with water to 12% v/v H₂SO₄ and autoclaved at 100°C for 3 h. The hydrolysates were filtered (10 µm) and the Klason lignin content

was determined as the weight of the residue after drying at 105°C for 24 h. The monomeric sugars glucose, xylose, and arabinose were determined using high-pressure liquid chromatography (HPLC). A Waters system, using a BioRad HPX-87H column at 40°C and 0.3 mL/min was used for analysis. All analyses were performed in triplicate.

## Characterization of the Biomass Compounds

### Rheology

Complex viscosity of the neat PLA and the compounds was measured on a AR 2000 rheometer (TA instruments, USA). The experiments were recorded in triplicate in frequency sweeps from 100 to 1 Hz at 200°C and 1% strain. The linear visco-elastic region was determined at 1 Hz.

### Molecular Weight

Gel Permeation Chromatography (GPC) experiments were carried out on a using a styrene divinylbenzene (SDVB) copolymer gel column set consisting of 2 analytical columns (SDV Lux Lin M; 5 µm; 300 × 8 mm; PSS Polymer Standards Service GmbH, Germany) and a precolumn (SDV analytical precolumn, 5 µm; 50 × 8 mm; PSS Polymer Standards Service GmbH, Germany). The flow rate was set to 1 mL/min, the system was conditioned to 30°C and an injection volume of 100 µL was used. Samples were prepared in triplicate by dissolving 6 to 8 mg of sample in chloroform to a concentration of 2.0 mg/mL overnight. The solutions were filtered through a 0.45 µm PTFE filter. The instrument (Knauer, Germany) was calibrated using 12 polystyrene standards (ReadyCal Polystyrene standards; PSS Polymer Standards Service GmbH, Germany) ranging from 370 to 2,520,000 g/mol. A conventional calibration curve was established by measuring the elution volume of the standards and plotting it against the logarithm of the molar mass. The calibration data was fitted with a fourth order polynomial function (r² = 1.00). Number average molecular weight (Mn), weight average molecular weight (Mw) and polydispersity index (PDI) were calculated from the detector response using the WinGPC UniChrom software (PSS Polymer Standards Service GmbH, Germany).

### Dynamic and Quasi-Static Mechanical Test

Specimens for dynamic mechanical thermal analysis (DMTA) and mechanical analysis in flexural mode were tested in 3-point bending mode on a RSA-G2 DMTA (TA instruments, USA) using a 40-mm span. Specimens were conditioned for 48 h at 23°C (±2°C) and 50% (± 5%) relative humidity.

Dynamic testing was carried out from 23 to 80°C at a rate of 5°C/min. Strain oscillations of 0.1% at a frequency of at 1 Hz were applied. The storage modulus (E') and tan δ were recorded throughout the test. Each formulation and printing direction was tested in triplicates.

Flexural testing was carried out on an Instron 5566 universal testing machine (Instron, USA) fitted with a 10 kN load cell at a rate of 2 mm/min. Flexural modulus, strength and strain were recorded throughout the test. At least 5 samples for each formulation were tested.
# RESULTS

## Characterization of the Biomass Particles

The wood powder particles were larger than those of the rice husks with a median size of 209.9 µm for the wood powder and 28.1 µm for rice husk powder (Table 2). Both powders exhibit a widespread distribution with a similar SPAN around 3.3. SPAN expresses the width of the distribution (D90-D10/D50). Rice husk powder contains a considerable proportion of fine particles below 20 µm whereas wood powder has a high content of coarse particle (above 400 µm) as observed in the shape of the particle size distribution graph (Figure 1). As a result, the specific surface area calculated from the particle size distribution is approximately 9 times higher for rice husk powder than for wood powder. However, the particle size distributions must be interpreted with caution as laser diffraction model calculations assume spherical homogeneous particles (Bohren and Huffman, 1983).

When looking at the biomass micro-morphology by SEM, large fibrous structures, >125 µm in length, are visible in the wood flour in contrast with low aspect ratio particles of the rice husk particles which supports the laser diffraction measurements (Figure 2). The presence of such particle sizes despite the sieving step is explained by their high aspect

ratio that allowed them to pass through the mesh. The size and morphological differences contribute to explain the lower apparent density of the wood biomass through a lesser packing efficiency.

In terms of chemical compositions determined by HPLC, the two biomasses are consistent with literature for lignocellulosic feedstocks (Table 3) (Müssig, 2010). The main differences between powders are (a) higher lignin content for the wood powder and (b) higher ash content for the rice husk powder (i.e., 15.7 wt.% compared to 0.16 wt.% for wood). The higher ash content in rice husks is explained by the presence of silica and contributes to the five-fold difference recorded in tapped density (Table 3) (Hamdan et al., 1997). The lack of a 100 % mass balance was attributed to the presence of water soluble extractives. Those were found to be present in higher quantities in wood powder than in rice husk powder.

The chemical compositions of the powders as determined by ¹³C solid-state NMR were consistent with lignocellulosic biomass. The major peaks observed in the wood powder were attributed to cellulose (C1-C6), lignin methoxyl (56 ppm) and aromatic/phenolic groups, with minor contributions from hemicelluloses (Figure 3A) (Gil and Neto, 1999). The rice husk powder showed similar signals with the addition of peaks attributed to amide carbonyls (~175 ppm) and signals predominantly from the Cβ/Cγ amino acid side chains (Cγ 15-25 ppm and Cβ 25-35 ppm) overlapping with signals from acetylated hemicelluloses (Figure 3B) (Bruker BioSpin, 2018). The amino acid signals were attributed to the known protein content present in rice husk (Juliano et al., 1987; Vadiveloo et al., 2009).

<table><caption>**TABLE 2** | Physical characterization of the wood and rice husk powders including particle size distribution indicators, specific surface area and densities (D10: 10th percentile, D50: Median Size, D90: 90th percentile).</caption><thead><tr><th></th><th>Wood powder</th><th>Rice husk powder</th></tr></thead><tbody><tr><td>D10 (µm)</td><td>48.4</td><td>3.5</td></tr><tr><td>D50 (µm)</td><td>209.9</td><td>28.1</td></tr><tr><td>D90 (µm)</td><td>748.1</td><td>95.2</td></tr><tr><td>Specific surface area</td><td>0.074</td><td>0.65</td></tr><tr><td>SPAN ((d90-d10)/d50)</td><td>3.26</td><td>3.33</td></tr><tr><td>Apparent density<br>(×10<sup>-2</sup> kg.L<sup>-1</sup>)</td><td>8.2</td><td>45.2</td></tr><tr><td>Tapped density<br>(×10<sup>-2</sup> kg.L<sup>-1</sup>)</td><td>10.0</td><td>56.9</td></tr></tbody></table>

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /><figcaption>**FIGURE 2** | SEM images of wood powder (a) and rice husk powder (b).</figcaption></figure>

<table><caption>**TABLE 3** | Chemical composition in mass percentage of the biomass.</caption><thead><tr><th></th><th>Wood powder</th><th>Rice husk powder</th></tr></thead><tbody><tr><td>Hemicellulose</td><td>17.1 ± 0.8</td><td>21.8 ± 0.7</td></tr><tr><td>Cellulose</td><td>44.3 ± 1.0</td><td>40.6 ± 1.0</td></tr><tr><td>Lignin</td><td>28.0 ± 0.8</td><td>19.4 ± 1.5</td></tr><tr><td>Ash (minerals)</td><td>0.16 ± 0.1</td><td>15.7 ± 0.7</td></tr><tr><td>Total</td><td>89.6 ± 2.7</td><td>97.5 ± 3.9</td></tr></tbody></table>
<figure><img src="image_6.png" /><figcaption>**FIGURE 3** | Solid-state NMR of the wood powder (A) and rice husk powder (B). P, protein; H, hemicelluloses; L, lignin; C, cellulose. Top left spectrum is a zoom in of the 10–40 ppm range.</figcaption></figure>

## Extrusion Process Stability

During the compounding step, the die pressure was significantly different between the wood and the rice husk powders. It ranged from 68 to 73 bars and 22 to 24 bars for the wood and rice husk compounds, respectively. Low die pressure, substantial die surge and color change were observed during rice husk compounding. The processing differences were initially attributed to the particle variation in apparent densities and sizes, however, it did not account for the variation in colors. It is known that compounding lignocellulose biomass in polymer can oxidize the carbohydrates contained in the biomass and lead to a darkening of the compounds depending on the processing conditions (i.e., high shear, high temperature) (Muniyasamy, 2013; Gallos et al., 2017). Yet the variation in color between the raw feedstocks and the compounds was more intense for the PLA rice compound than those of the PLA wood, and accentuated in the FDM filament (Figure 4). Ball milling of lignocellulosic biomass during an extended period (i.e., >120 min) is known to increase surface area, decrease cellulose crystallinity and decrease the carbohydrate's degree of polymerization that consequently increases the particle surface reactivity (Vaidya et al., 2016; Gao et al., 2017). This reactivity could explain the difference of coloration caused by chemical reactions triggered by the rice husk compared to the wood powder.

During the FDM filament manufacture, the FDM rice husk filament exhibited a lower melt strength and material surged at the die causing unevenness in the filament shape. These processing effects made the rice husk filament more difficult to manufacture.

<figure><img src="image_7.png" /><figcaption>**FIGURE 4** | Photographs of the wood flour (top) and rice husk (bottom) FDM filaments.</figcaption></figure>

## Influence of the Biomass on the PLA

In terms of rheology, the measurement of the complex viscosities of PLA and the PLA/biomass compounds as a function of frequency indicated that (a) the PLA-wood compound has a higher viscosity than neat PLA and that (b) both are exhibiting shear thinning behaviors (Figure 5). The findings are consistent with literature on biomass-filled polymers and the general effect of biomass in polymer melt (Li and Wolcott, 2004; Bettini et al., 2013; Gallos et al., 2017; Chun et al., 2018). Conversely,
the rice compound was half the complex viscosity of neat PLA and indicated a poor reproducibility (Figure 5). It was hypothesized that the rice husk particles acted as a solid lubricant in the polymer matrix. However, the lack of reproducibility also suggested additional effects such as polymer hydrolysis of the backbone.

When comparing the PLA Mn and Mw of the biomass compounds, the GPC analysis indicated no significant difference between the two fillers (Table 4). More specifically, the GPC analysis gave a Mw and Mn of 98,020 g/mol and 52,473 g/mol for the neat PLA. After extrusion with biomass powders, significantly lower average molecular weights were recorded, i.e., Mw = 66,790 g/mol and Mn = 33,780 g/mol for PLA-wood and Mw = 59,853 g/mol and Mn = 31,463 for PLA-rice (Table 4, Figure 6). This finding is consistent with literature, as PLA is known to hydrolyse in the presence of moisture or acid groups in the biomass during extrusion (Le Guen et al., 2017; Thumm et al., 2018). This translated to an average molecular weight decrease of between 32 and 39% for the PLA compounded with wood and rice husk powders. No significant difference was observed between fillers. The reduction was also observable

in the molecular weight distribution (MWD) profiles (Table 4). Similarly, no significant changes in the polydispersity indices (PDI) of pure PLA compared to PLA with biomasses were observed (Table 4).

When comparing the effect of the fillers onto PLA, the GPC analysis indicated that both resulted in similar level of PLA hydrolysis, hence, the difference rheological behavior could not be attributed to a difference in Mn/Mw.

## Mechanical Performances of the 3D Prints at 0° and 90°

The infill of the print direction was carried out at 0° and 90° of the length of the rectangular beam (Figure 7). The printing direction influenced the mechanical properties of all formulations, showing a decrease of flexural modulus and strength for the 90° printed beams (Figure 8). The anisotropy induced difference in properties between the lengthwise and widthwise printing direction that ranged between 28 and 41% variation in flexural strength, 4 to 15% in strain and 27 to 33% in modulus. Similar mechanical properties results were reported for 3D printed PLA wood rectangular beam where the printing direction (i.e., 0° and 90° of the longitudinal axis of the beam) influenced the tensile stiffness and strength by 20 and 35%, respectively (Le Duigou et al., 2016). They suggested that the mechanical properties differences between the lengthwise and widthwise printing directions were due to the lack of interlayer interactions.

No significant difference was found between the wood and the rice husk fillers in terms of mechanical property except for the flexural stiffness in the lengthwise direction where the PLA-wood sample's modulus is ~25% higher than the modulus of the PLA-rice samples. The decrease in strength compared to neat PLA and the lack of difference between the PLA-wood and PLA-rice compounds, indicated that the mechanical properties are effectively driven by the PLA inter-strands coalescence as observed by others (Le Duigou et al., 2016). This coalescence

<figure><img src="image_8.png" /><figcaption>**FIGURE 5** | Complex viscosity as a function of frequency for: PLA (Square), PLA-wood (Circle) and PLA-rice (Triangle). Triplicate measurements were carried out and all results are presented on the graph.</figcaption></figure>

<table><caption>**TABLE 4** | Molecular weight values and polydispersity indices of the PLA neat and compounded with biomasses.</caption><thead><tr><th></th><th>Mn<br>(g/mol)</th><th>St. Dev.<br>Mn (g/mol)</th><th>Mw<br>(g/mol)</th><th>St. Dev.<br>Mw (g/mol)</th><th>PDI<br>(-)</th><th>St. Dev.<br>PDI (-)</th></tr></thead><tbody><tr><td>PLA</td><td>52,473</td><td>5,245</td><td>98,020</td><td>6,294</td><td>1.9</td><td>0.24</td></tr><tr><td>PLA-wood</td><td>33,780</td><td>2,735</td><td>66,790</td><td>2,776</td><td>2.0</td><td>0.08</td></tr><tr><td>PLA-rice</td><td>31,463</td><td>3,526</td><td>59,853</td><td>5,963</td><td>1.9</td><td>0.03</td></tr></tbody></table>

n = 3.

<figure><img src="image_9.png" /><figcaption>**FIGURE 6** | Molecular weight distribution profiles of neat PLA (solid line), PLA-wood (dotted line) and PLA-rice (dashed line) compounds.</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>**FIGURE 7** | Typical 3D printed rectangular beams for flexural testing. From bottom to top PLA 0°, PLA 90°, PLA-rice 0°, PLA-rice 90°, PLA-wood 0° and PLA-wood 90°.</figcaption></figure>

is expected to decrease with the introduction of a solid filler that hinders the strand cohesion and creates weak points within the printed object. Singamenni et al., reported a gradual loss of inter-strand and inter-layer coalescence surface, measured by optical microscopy, with increasing amount of wood flour from 5 to 20 wt.% in polybutyrate-adipate-terephthalate-polymer matrix (PBAT) (Singamenni et al., 2018b). They attributed the decrease in flexural properties to the lower strands coalescence (98% coalescence for 5wt.% filler and 64% for 20wt.%). Their findings in terms of flexural strength are consistent with those presented in this study, where the addition of either filler decreases the flexural strength of the printed samples. However, their results relative to the flexural modulus differs. They reported a decrease in flexural modulus at 10 wt.% wood fraction whereas the present results indicate an increase of flexural modulus for the PLA-wood and decrease for the PLA-rice.

The difference was attributed to the polymer matrix as PBAT is more elastomeric than PLA but it raises the possibility of synergetic/competing effects between the type of filler and polymer, leading to different mechanical behavior of the 3D printed object. Similar findings were reported by Badouard et al. who compared the reinforcing ability of flax fibers in different biopolymer matrices (Badouard et al., 2019).

<figure><img src="image_11.png" /><figcaption>**FIGURE 8** | Mechanical properties comparison of the 3D printed beams PLA, (Top) Flexural modulus (Middle) Flexural strength and (Bottom) Flexural strain (Error bars indicate the ± 95 % confidence interval).</figcaption></figure>

The dynamic thermo-mechanical behavior of the samples tested from 23 to 80°C indicated a single transition, typical of the glass transition temperature of PLA samples (Table 5).

The stiffness variation due to the printing direction anisotropy agreed with the quasi-mechanical tests for the neat PLA and filled samples. Storage modulus (E') values indicated a decrease of 28.7, 25.1, and 19.7 % between 0° and 90° printing direction for respectively, the unfilled PLA, PLA-wood and PLA-rice samples, respectively (Table 5). The printing direction had no thermo-mechanical effect (e.g., Tg variation, thermal stress relaxation).
<table><caption>**TABLE 5** | Storage modulus values at 30 and 75°C, and glass transition temperature (Tg) for PLA, PLA-wood, and PLA-rice samples printed at 0 and 90° of the longitudinal direction.</caption><thead><tr><th></th><th>PLA<br>0°</th><th>PLA<br>90°</th><th>PLA-wood<br>0°</th><th>PLA-wood<br>90°</th><th>PLA-rice<br>0°</th><th>PLA-rice<br>90°</th></tr></thead><tbody><tr><td>E' at 30°C (GPa)</td><td>2.66 (0.22)</td><td>1.90 (0.15)</td><td>2.31 (0.07)</td><td>1.73 (0.16)</td><td>2.50 (0.15)</td><td>2.01 (0.43)</td></tr><tr><td>Tg (°C)</td><td>65.0 (0.35)</td><td>65.5 (0.35)</td><td>64.7 (1.7)</td><td>65.4 (1.8)</td><td>65.3 (0.9)</td><td>64.8 (0.4)</td></tr><tr><td>E' at 75°C (MPa)</td><td>5.6 (0.5)</td><td>3.9 (1.6)</td><td>12.8 (2.7)</td><td>7.8 (1.1)</td><td>8.6 (1.7)</td><td>13.4 (1.2)</td></tr></tbody></table>

The number in brackets is the standard deviation for n = 3.

However, the filler's contribution to the stiffness could be observed past the glass transition temperature where the storage modulus of the PLA-wood and PLA-rice samples was higher than the unfilled PLA samples. The retention of storage modulus indicated opposite trends between the PLA-wood samples 0° and 90°, and the PLA-rice samples 0° and 90°. Indeed, past the Tg, the wood filled PLA samples displayed a 39% difference in favor of the longitudinal printing direction. Conversely, the rice-husk filled samples displayed an unexpected 56% difference in favor of the transversal printing direction. The former finding is explained by an alignment of wood particles in the printing direction caused by shear forces in the nozzle (Sydney Gladman et al., 2016). The latter finding is more difficult to explain. It was here tentatively attributed to the smaller size of the rice husk particles and their better distribution within the PLA. The rice husk particles are around 10 times smaller than the wood ones (Figure 2, Table 2). Considering the low aspect ratio of the rice husk particles, it is conceivable that no or little alignment for the rice husk particles occurred. The difference in behavior between the rice husk and the wood was tentatively attributed to the finer particles which created a more uniform and stiffening network. This characteristic was only observed in the widthwise direction as the predominant factor (i.e., the alignment of the polymer chains and fillers) was potentially not present to mask it.

## DISCUSSION

Overall, the mechanical behavior of the 3D printed beams at room temperature was governed by the printing direction and the polymer physical properties rather than the type of fillers being used. However, when comparing the wood powder to the rice husk powder as a filler for 3D printing in PLA, particular observations were made on the interaction between the biomasses and PLA.

## Process Stability and Thinning Behavior

Whilst both biomass fillers caused a reduction in the PLA molecular chain length, the extrusion process was more unstable during rice husk powder compounding. GPC analysis indicated that the degree of hydrolysis during processing was slightly higher for the rice husk powder than the wood powder, however the difference was not significant and could not alone justify the rheological differences. The complex viscosity analysis indicated a dramatic thinning effect of the rice husk particles on the PLA (~2.5 times lower at 1 Hz). Thinning/plasticization additives is fairly common in polymer processing however, it is much less common

in the case of solids in the 10 to 100's micron size range. Coppola et al. found that the addition of 3 and 5 vol.%, powdered hemp shives in PLA reduced the complex viscosity of the compound compared to neat polymer (Coppola et al., 2018). This behavior was attributed to chemical degradation of the pectin contained in hemp shives. Conversely, Cipriano et al. reported the same effect on PLA caused by magnesium silicate particles (Cipriano et al., 2014). They attributed this lubricant effect to a physical interaction between the particles and the polymer chains.

In this study, both chemical and physical phenomena are valid. Rice husk is composed of 15.2 wt.% amorphous silicate which contributes to its higher density compared to wood or other lignocellulosics (e.g., flax) and could contribute to the rheological thinning behavior as reported Cipriano et al. (2014). However, the concurrent presence of organic materials (e.g., carbohydrates) in the particles also points to a chemical triggered mechanism. The difference in chemical reactivity between wood and rice powders was attributed to the particles size and the ball milling. Ball milling of biomass is known to break down the carbohydrates which increases their reactivity (e.g., better enzymatic adsorption) (Vaidya et al., 2016; Gao et al., 2017). The possibility of reacting these carbohydrates at temperatures above 150°C in an extruder is highly probable leading to the faster production of oxidation products such as alcohols, acids or other volatile products (Formela et al., 2018).

## Reactive Extrusion of Lignocellulosic Biomass

During compounding of biomass in PLA, hydrolysis of the polymer chains is known to occur (Le Guen et al., 2017; Thumm et al., 2018). A previous study comparing the mechanical behavior of PLA reinforced with protein-based filler and wood flour reported that under accelerated weathering, the protein-based filler hydrolysed the PLA matrix at least twice as fast as a wood flour filler due to the amino-acids produced from protein hydrolysis (Le Guen et al., 2017). Rice husks contain more proteins than wood (Hamdan et al., 1997; Vadiveloo et al., 2009), however, the GPC analysis indicated comparable molecular weights for the PLA compounded with either wood or rice husk powders. This finding ruled out a potential difference solely based on acid hydrolysis of the PLA chains.

In the case of the rice husk compound, the darkening color could be attributed to two chemical reactions:

- Maillard reaction between the amino acids and the carbohydrates (e.g., xylan), that has been previously reported
to occur during processing of the rice husks at 180°C (Vegas et al., 2008),

- Carbohydrate oxidation (i.e. caramelization) that is known to form volatile organic molecules (Jiang et al., 2008).

In both cases, the reactions release volatile compounds and water that can explain the instability of the rice husk compounding process and the plasticization of the PLA matrix. Whilst the oxidation reactions also occur with wood powder, the higher reactivity of the rice husk powder compared to the wood powder was attributed to the ball-milling pre-treatment.

## CONCLUSIONS

Wood and rice husk powders were blended by twin-screw extrusion with PLA to be used for 3D printing applications. It was found that the mechanical properties of the 3D printed samples at room temperature were mainly driven by the properties of PLA and anisotropy was only influenced by the 3D printing directionality. In addition, the mechanical properties were lowered by the integration of either filler which was attributed to a lowering of inter-strand cohesion in the 3D printed samples. Past the glass transition temperature of the PLA, the fillers enabled the retention of stiffness measured by the storage modulus.

In terms of visual assessment and process stability however, differences were observed and attributed to chemical reactions triggered by the (a) silica content of rice husk and/or (b) chemical reactions including Maillard reaction and carbohydrate oxidations. These chemical reactions were potentially enhanced by a particle size difference (i.e., rice husk powder was smaller and denser than the wood powder) and triggered during the temperature and shearing conditions of the extrusion process. Little polymer research has been conducted on the potential of using reactive biomass filler particles based on mechano-chemistry (i.e., ball milling). These

unique findings indicate that biomass can trigger chemical reactions that change the polymer rheological behaviors during processing dramatically and reinforces the need for careful attention when matching a biopolymer matrix with a biomass filler.

The ball-milled rice husk powder's effect on thinning PLA, if controlled during the compounding process, is of major interest in the development of reactive sustainable plastic additives.

## DATA AVAILABILITY STATEMENT

All datasets generated for this study are included in the manuscript/supplementary files.

## AUTHOR CONTRIBUTIONS

M-JL and CM-L contributed to the design of the research. All authors contributed to its implementation and the analysis of the results as per field of expertise and wrote the manuscript.

## ACKNOWLEDGMENTS

The authors thank the French GDR Symbiose, and the New Zealand Ministry of Business, Innovation, and Employment funding under High Value Manufacturing and Services (HVMS) Enabling Technologies investment contract C04X1205.

The authors are gratefully acknowledging the contribution of Mr. Ross Anderson and Mr. Bryce McQuillan for supplying and preparing the wood powder. The authors also want to acknowledge the team of the plant-based platform of the JRU IATE, particularly Charlène Fabre for preparing the rice husk powder. In addition, the authors want to express their gratitude to Dr. Mathias Sorieul, Mr. Armin Thumm, and Mr. Elouan Guillou for their helpful advice.

## REFERENCES

Badouard, C., Traon, F., Denoual, C., Mayer-Laigle, C., Paës, G., and Bourmaud, A. (2019). Exploring mechanical properties of fully compostable flax reinforced composite filaments for 3D printing applications. *Indus. Crops Products* 135, 246–250. doi: 10.1016/j.indcrop.2019.04.049

Bettini, S. H. P., de Miranda Josefovich, M. P. P., Muñoz, P. A. R., Lotti, C., and Mattoso, L. H. C. (2013). Effect of lubricant on mechanical and rheological properties of compatibilized PP/sawdust composites. *Carbohyd. Polymers* 94, 800–806. doi: 10.1016/j.carpol.2013.01.080

Bohren, C. F. and Huffman, D. R. (1983). *Absorption and Scattering of Light by Small Particles*. New York, NY: Wiley.

Bruker BioSpin (2018). *TopSpin 4.0.5 NMR Guide [Computer Software]*. Bruker BioSpin.

Chun, K. S., Husseinsyah, S., and Yeng, C. M. (2018). Torque rheological properties of agro-waste-based polypropylene composites: effect of filler content and green coupling agent. *J. Vinyl Addit. Technol.* 24, E30–E35. doi: 10.1002/vnl.21561

Chungsangunsit, T. (2009). Emission assessment of rice husk combustion for power production. *Int. J. Mech. Aerospace, Indus. Mechatr. Manufact. Eng.* 3. Available online at: https://pdfs.semanticscholar.org/7454/b2f04cd9b56de336abe6ac904f100519f123.pdf

Cipriano, T. F. S., da Silva, A. L. N., da Silva, A. H. M., de Sousa, A. M. F., da Silva, G. M., and Rocha, M. G. (2014). Thermal, rheological and morphological properties of poly (Lactic Acid) (PLA) and talc composites. *Polímeros* 24:2. doi: 10.17265/2161-6221/2013.11.002

Coppola, B., Garofalo, E., Maio, L. D., Scarfato, P., and Incarnato, L. (2018). Investigation on the use of PLA/hemp composites for the fused deposition modelling (FDM) 3D printing. *AIP Conference Proc.* 1981:020086. doi: 10.1063/1.5045948

Formela, K., Zedler, Ł., Hejna, A., and Tercjak, A. (2018). Reactive extrusion of bio-based polymer blends and composites – Current trends and future developments. *Exp. Polymer Lett.* 12, 24–57. doi: 10.3144/expresspolymlett.2018.4

Gallos, A., Paës, G., Allais, F., and Beaugrand, J. (2017). Lignocellulosic fibers: a critical review of the extrusion process for enhancement of the properties of natural fiber composites. *RSC Adv.* 7, 34638–34654. doi: 10.1039/C7RA05240E

Gao, C., Xiao, W., Ji, G., Zhang, Y., Cao, Y., and Han, L. (2017). Regularity and mechanism of wheat straw properties change in ball milling process at cellular scale. *Bioresour. Technol.* 241, 214–219. doi: 10.1016/j.biortech.2017.04.115

Gibson, M. A., Mykulowycz, N. M., Shim, J., Fontana, R., Schmitt, P., Roberts, A., et al. (2018). 3D printing metals like thermoplastics: fused
filament fabrication of metallic glasses. *Mater. Today* 21, 697–702. doi: 10.1016/j.mattod.2018.07.001

Giddel, M. R., and Jivan, A. P. (2007). "Waste to wealth, potential of rice husk in india a literature review," in *International Conference on Cleaner Technologies and Environmental Management PEC, Pondicherry*, 586–590.

Gil, A. M., and Neto, C. P. (1999). Solid-state NMR studies of wood and other lignocellulosic materials. *Ann. Rep. NMR Spectrosc.* 37, 75–117. doi: 10.1016/S0066-4103(08)60014-9

Graichen, F. H. M., Grigsby, W. J., Hill, S. J., Raymond, L. G., Sanglard, M., Smith, D. A., et al. (2017). Yes, we can make money out of lignin and other bio-based resources. *Indus. Crops Products* 106, 74–85. doi: 10.1016/j.indcrop.2016.10.036

Hamdan, H., Muhid, M. N. M., Endud, S., Listiorini, E., and Ramli, Z. (1997). 29Si MAS NMR, XRD and FESEM studies of rice husk silica for the synthesis of zeolites. *J. Non-Crystalline Solids* 211, 126–131. doi: 10.1016/S0022-3093(96)00611-4

Jiang, B., Liu, Y., Bhandari, B., and Zhou, W. (2008). Impact of caramelization on the glass transition temperature of several caramelized sugars. Part I: chemical analyses. *J. Agric. Food Chem.* 56, 5138–5147. doi: 10.1021/jf703791e

Juliano, B. O., Maniñat, C. C., Pascual, G. C. (1987). Properties of fraction of rice hull. *Phytochemistry* 26, 3261–3263. doi: 10.1016/S0031-9422(00)82483-8

Kumar, A., Mohanta, K., Kumar, D., and Parkash, O. (2012). Properties and industrial applications of rice husk: a review. *Int. J. Emerg. Technol. Adv. Eng.* 2, 86–90.

Lammi, S., Le Moigne, N., Djenane, D., Gontard, N., and Angellier-Coussy, H. (2018). Dry fractionation of olive pomace for the development of food packaging biocomposites. *Indus. Crops Products* 120, 250–261. doi: 10.1016/j.indcrop.2018.04.052

Le Duigou, A., Castro, M., Bevan, R., and Martin, N. (2016). 3D printing of wood fibre biocomposites: from mechanical to actuation functionality. *Mater. Design* 96, 106–114. doi: 10.1016/j.matdes.2016.02.018

Le Guen, M.-J., Thoury-Monbrun, V., Castellano Roldán, J. M., and Hill, S. J. (2017). Assessing the potential of farm dairy effluent as a filler in novel PLA biocomposites. *J. Polymers Environ.* 25, 419–426. doi: 10.1007/s10924-016-0824-1

Le Guen, M. J., and Newman, R. H. (2007). Pulped Phormium tenax leaf fibres as reinforcement for epoxy composites. Composites Part A. *Appl. Sci. Manufact.* 38, 2109–2115. doi: 10.1016/j.compositesa.2007.07.001

Li, T. Q., and Wolcott, M. P. (2004). Rheology of HDPE-wood composites. I. Steady state shear and extensional flow. Composites Part A. *Appl. Sci. Manufact.* 35, 303–311. doi: 10.1016/j.compositesa.2003.09.009

Liles, K. P., Greene, A. F., Danielson, M. K., Colley, N. D., Wellen, A., Fisher, J. M., Barnes, J. C. Photoredox-Based Actuation of an Artificial Molecular Muscle. *Macromol. Rapid Commun.* (2018) 39:1700781. doi: 10.1002/marc.201700781

Lv, H., Song, S., Sun, S., Ren, L., and Zhang, H. (2016). Enhanced properties of poly(lactic acid) with silica nanoparticles. *Polymers Adv. Technol.* 27, 1156–1163. doi: 10.1002/pat.3777

Malvern (2007). *Sample Dispersion and Refractive Index Guide - Mastersizer* 2000. Malvern.

Mirko, K., Sernek, M., and Kitek Kuzman, M. (2016). Use of wood powder and adhesive as a mixture for 3D printing. *Eur. J. Wood Products* 74:123. doi: 10.1007/s00107-015-0987-9

Montalvo Navarrete, J. I., Hidalgo-Salazar, M. A., Escobar Nunez, E., and Rojas Arciniegas, A. J. (2018). Thermal and mechanical behavior of biocomposites using additive manufacturing. *Int. J. Interact. Design Manufact.* 12, 449–458. doi: 10.1007/s12008-017-0411-2

Muniyasamy, S. (2013). Biodegradability and compostability of lignocellulosic based composite Materials. *J. Renew. Mater.* 1:253.

Müssig, J. (2010). *Industrial Applications of Natural Fibres: Structure, Properties and Technical Applications*. John Wiley & Sons, Ltd. doi: 10.1002/9780470660324

Newman, R. H., Hemmingson, J. A., and Suckling, I. D. (1993). Carbon13 nuclear magnetic resonance studies of kraft pulping. *Holzforschung*, doi: 10.1515/hfsg.1993.47.3.234

Pode, R. (2016). Potential applications of rice husk ash waste from rice husk biomass power plant. *Renew. Sustain. Energy Rev.* 53, 1468–1485. doi: 10.1016/j.rser.,2015.09.051

Revelo, C. F., and Colorado, H. A. (2018). 3D printing of kaolinite clay ceramics using the Direct Ink Writing (DIW) technique. *Ceram. Int.* 44, 5673–5682. doi: 10.1016/j.ceramint.2017.12.219

Shao, Y., Guizani, C., Grosseau, P., Chaussy, D., and Beneventi, D. (2018). Use of lignocellulosic materials and 3D printing for the development of structured monolithic carbon materials. *Compos. Part B Eng.* 149, 206–215. doi: 10.1016/j.compositesb.2018.05.035

Singamneni, S., Behera, M. P., Le Guen, M., and Zeidler, H. (2018a). Mechanism of bonding in seashell powder based ceramic composites used for binder-jet 3D printing. *Bioceram. Dev. Appl.* 108. Available online at: https://www.omicsonline.org/open-access/mechanism-of-bonding-in-seashell-powder-based-ceramic-compositesused-for-binderjet-3d-printing-2090-5025-1000108.pdf

Singamneni, S., Smith, D., LeGuen, M.-J., and Truong, D. (2018b). Extrusion 3D printing of polybutyrate-adipate-terephthalate-polymer composites in the pellet form. *Polymers* 10:922. doi: 10.3390/polym10080922

Sydney Gladman, A., Matsumoto, E. A., Nuzzo, R. G., Mahadevan, L., and Lewis, J. A. (2016). Biomimetic 4D printing. *Nat. Mater.* 15:413. doi: 10.1038/nmat4544

Tao, Y., Wang, H., Li, Z., Li, P., and Shi, S. Q. (2017). Development and application of wood flour-filled polylactic acid composite filament for 3D printing. *Materials* 10:339. doi: 10.3390/ma10040339

Thumm, A., Even, D., Gini, P.-Y., and Sorieul, M. (2018). Processing and properties of MDF fibre-reinforced biopolyesters with chain extender additives. *Int. J. Polymer Sci.* 2018:9. doi: 10.1155/2018/9601753

Ummah, H., Suriamihardja, D., Selintung, M., and Wahab, A. W. (2015). Analysis of chemical composition of rice husk used as absorber plates sea water into clean water. *ARPN J. Eng. Appl. Sci.* 10, 6046–6050.

Vadiveloo, J., Nurfariza, B., and Fadel, J. G. (2009). Nutritional improvement of rice husks. *Anim. Feed Sci. Technol.* 151, 299–305. doi: 10.1016/j.anifeedsci.2009.03.002

Vaidya, A. A., Donaldson, L. A., Newman, R. H., Suckling, I. D., Campion, S. H., Lloyd, J. A., et al. (2016). Micromorphological changes and mechanism associated with wet ball milling of Pinus radiata substrate and consequences for saccharification at low enzyme loading. *Bioresour. Technol.* 214, 132–137. doi: 10.1016/j.biortech.2016.04.084

Vegas, R., Kabel, M., Schols, H. A., Alonso, J. L., and Parajó, J. C. (2008). Hydrothermal processing of rice husks: effects of severity on product distribution. *J. Chem. Technol. Biotechnol.* 83, 965–972. doi: 10.1002/jctb.1896

Wang, Q., Sun, J., Yao, Q., Ji, C., Liu, J., and Zhu, Q. (2018). 3D printing with cellulose materials. *Cellulose* 25, 4275–4301. doi: 10.1007/s10570-018-1888-y

Yalçin, N., and Sevinç, V. (2001). Studies on silica obtained from rice husk. *Ceram. Int.* 27, 219–224. doi: 10.1016/S0272-8842(00)00068-7

Zeidler, H., Klemm, D., Böttger-Hiller, F., Fritsch, S., Le Guen, M. J., and Singamneni, S. (2018). 3D printing of biodegradable parts using renewable biobased materials. *Proc. Manufact.* 21, 117–124. doi: 10.1016/j.promfg.2018.02.101

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2019 Le Guen, Hill, Smith, Theobald, Gaugler, Barakat and Mayer-Laigle. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.