Subscriber access provided by BUFFALO STATE
Article
Measuring molecular strain in rewetted and neverdried eucalypt wood with Raman spectroscopy Fei Guo, and Clemens M. Altaner Biomacromolecules, Just Accepted Manuscript • DOI: 10.1021/acs.biomac.9b00808 • Publication Date (Web): 17 Jul 2019 Downloaded from pubs.acs.org on July 21, 2019
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biomacromolecules
1
Measuring molecular strain in rewetted and never-
2
dried eucalypt wood with Raman spectroscopy
3
Fei Guo1, 2, Clemens M. Altaner2*
4
AUTHOR ADDRESS
5
1 College
6
China
7
2 New
8
New Zealand
9
KEYWORDS: band shift, E. regnans, E. quadrangulata, growth stress, molecular deformation,
10
of Furnishings and Industrial Design, Nanjing Forestry University, 210037, Nanjing,
Zealand School of Forestry, University of Canterbury, Private Bag 4800, Christchurch,
tensile strain, PLS.
11
1 Environment ACS Paragon Plus
Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
12
ABSTRACT
13
To measure growth strain in wood using Raman spectroscopy, we investigated the Raman spectra
14
of rewetted (water-saturated) E. regnans and green E. quadrangulata wood during tensile and
15
four-point bending tests. Partial least squares models to predict tensile strain were built from the
16
Raman spectra. The best model could predict tensile strain with a root mean square error of 427.5
17
>?: Apart from the widely reported band shift at 1095 cm@ upon mechanical strain, spectral
18
changes at 1420, 1120, 895, 456 cm@ were identified. The assignments of these bands were
19
discussed in relation to the molecular deformation of cellulose. The band shift rates during tensile
20
tests were @ :,. and @ : $ cm@ /% for rewetted E. regnans and green E. quadrangulata wood,
21
respectively. We successfully detected the release of molecular growth strain in green eucalyptus
22
wood with Raman spectroscopy by observing band shifts of the 1095 cm@ signal. Further, there
23
was a moderate correlation (r = 0.48) between growth-strain measured with strain gauges and the
24
1095 cm@ band position. The precision of the prediction of growth strain using Raman
25
spectroscopy was negatively affected by variation attributed to the inhomogeneity of wood on the
26
millimeter scale and instrumental instability.
27 28
1 INTRODUCTION
29
Raman and mid-infrared (MIR) spectroscopy are complementary techniques measuring
30
vibrational frequencies of chemical bonds 1. The former is based on the Raman scattering
31
phenomenon, while the latter detects the absorption of electromagnetic radiation. With
32
technological advances, Raman spectroscopy is gaining popularity in various applications
33
including material identification 2, forensic and pharmaceutical analysis 3 as well as chemical
2 Environment ACS Paragon Plus
Page 2 of 37
Page 3 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biomacromolecules
34
distribution imaging of plant cell walls 4, 5. Portable Raman spectrometers are commercially
35
available 6, opening up the possibility of measuring wood properties in the field rapidly and non-
36
destructively. Compared to infrared spectroscopy, Raman spectroscopy is insensitive to water,
37
which interferes with the structural information of plant material 7 as it is typically water
38
saturated when alive. Challenges of using Raman spectroscopy on lignocellulosic materials are
39
the laser-excited florescence from lignin, which causes an intense background masking the
40
Raman signal, and the high energy of the laser, which leads to thermal decomposition of the
41
sample. Using 1064 nm near infrared (NIR) lasers dramatically reduces florescence from lignin 8
42
and water immersion sampling lessens thermal degradation 9.
43
Mechanical stress induces both macroscopic and molecular deformation of wood and cellulose.
44
Axial elongation and lateral contraction of the cellulose crystal lattice in wood was observed
45
upon tensile stress by X-ray diffraction 10, 11. Deformation of crystalline cellulose will cause
46
changes in vibrational frequencies of chemical bonds, which were observed as band shifts in
47
MIR 12, 13, near infrared (NIR) 14 and Raman spectra 15. IR spectroscopy detected the deformation
48
of the COC glycosidic bonds at 1160 cm@ and intramolecular O3H O5 groups at 3348 cm@
49
NIR band shifts at 6286 ± 5 cm@ and 6470 ± 10 cm@ indicated intramolecular hydrogen bonds
50
of crystalline cellulose elongated under tensile strain 14.
51
Raman spectroscopy has been a powerful tool to study the molecular deformation of polymer
52
fibers. 16 reported that the Raman frequencies assigned to the backbone structure of
53
polydiacetylene decreased with the elongation of fibers. Following this work, Raman
54
spectroscopy was used to study the molecular deformation of regenerated cellulose 17. The
55
Raman band shift at around 1095 cm@ , assigned to the stretching of the cellulose ring structure,
3 Environment ACS Paragon Plus
13.
Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
56
has also been reported for wood, flax and hemp fibers 18-20. However, the effect of mechanical
57
strain on other regions in the Raman spectra is not well understood.
58
Growth stresses in trees, accumulated during cell wall maturation 21, cause problems like end-
59
splitting, heart checking and collapse and are the major cause restricting plantation-grown
60
eucalypts from in solid wood processing 22. Growth strain, similar to mechanical strain, causes
61
deformation of wood and cellulose on the molecular level 23. Methods to measure growth strain
62
in trees, the Nicholson technique 24, the CIRAD-Forêt method 25, log splitting 26 and the strain
63
gauge method 27, are destructive or time consuming. Rapid and non-destructive evaluation of
64
growth stress is needed to segregate logs with low growth strain for high-value solid wood uses.
65
Raman spectroscopy could potentially measure growth strain quickly and non-destructively by
66
accurate determination of vibrational frequencies. Previous studies on the deformation
67
mechanism of cellulosic materials focused on small areas of single fibers using a Raman
68
microscope, requiring elaborate sample preparation. However, Raman spectrometers capable of
69
direct surface sampling are available, enabling the collection of spectra without the limitation of
70
sample size. This is also advantageous as mechanical properties of wood have been known to
71
depend on sample size 28. Thin wood samples were shown to have different load-displacement
72
curves compared to bulk wood 29 and the modulus of elasticity of microtomed wood sections
73
were lower than that of normal sized specimens 30.
74
In this study, thick wood strips were mechanically stretched to confirm the effect of strain on
75
Raman spectra, previously only reported for microscopic samples, differing in their mechanical
76
behavior. Band shift rates were calculated to quantify molecular strain in rewetted and green
77
wood. Following on from this, the release of growth strain was monitored by Raman
4 Environment ACS Paragon Plus
Page 4 of 37
Page 5 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biomacromolecules
78
spectroscopy on green, never-dried stems before finally evaluating the possibility of measuring
79
growth strain in green stems using Raman spectroscopy.
80 81
2 MATERIALS AND METHODS
82
2.1 Tensile tests of rewetted and green wood monitored by Raman spectroscopy
83
An air-dry Eucalyptus regnans board with an oven-dry density of 0.51 g/cm3 was sourced from
84
Canterbury, New Zealand. The microfibril angle (MFA) of the E. regnans wood was 7°
85
(standard deviation of 2.9°), as previously measured by X-ray diffraction. Ten wood strips with a
86
thickness of 0.5 to 1 mm (radial) were cut with a circular saw from the sapwood of the air-dry E.
87
regnans board and sanded before being cut into the shape shown in Figure 1. Wood sample used
88
for tensile tests while collecting Raman spectra. The grain (longitudinal direction) runs along the
89
longest dimension of the sample. Dimensions are labelled in millimeters.
90
. The samples were re-saturated through immersion in water for at least 7 days and were referred
91
to as “rewetted” wood in this study. A 2-year old E. quadrangulata tree with an oven-dry density
92
of 0.68 g/cm3 was harvested from a site in Woodville, New Zealand. The tree was kept in water
93
during transport to the lab. Eight wood strips were prepared from the never-dried green tree in
94
accordance to the procedure described above for the E. regnans samples. E. quadrangulata
95
samples were put in a sealed plastic bag and stored in a freezer until tested.
5 Environment ACS Paragon Plus
Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
96 97
Figure 1. Wood sample used for tensile tests while collecting Raman spectra. The grain
98
(longitudinal direction) runs along the longest dimension of the sample. Dimensions are labelled
99
in millimeters.
100
A purpose-built tensile test rig 14 was placed under a dispersive Raman spectrometer (Sierra 2.0,
101
Metrohm Raman, USA) equipped with a 1064 nm NIR laser and an InGaAs array detector. The
102
wood samples were mechanically stretched while simultaneously collecting Raman spectra at
103
twelve strain levels, increasing progressively from “Strain 0” (approx. 0 >?< to “Strain 11”
104
(maximum strain of approx. 4500 >?< with increments of approximately 410 >?) before returning
105
to the relaxed (RL) state. Applied macroscopic strain was measured by strain gauges (FLA-5-11,
106
Tokyo Sokki Kenkyujo, Japan) glued (Loctite 454, Australia) onto the samples and connected to
107
a TC-31K strain-meter (Tokyo Sokki Kenkyujo, Japan). Strain was defined as the amount of
108
deformation per unit length and 1 >? equals to 1×10@. strain. The wood surface was kept wet
109
during spectra collection to restrict thermal decomposition 31 and improve spectral quality 9. Two
110
Raman spectra were taken at each strain level in the region from 2300 to 200 cm@ with a laser
111
power of 500 mW and an integration time of 30 to 35 s. In total 26 spectra were acquired for
112
each sample. The spectrometer featured orbital raster scan (ORS), a technique to reduce sample
113
damage by moving the laser over a measuring area of 1 to 2 mm in diameter. The polarization
114
direction of the incident light was parallel with the longitudinal direction of the wood samples.
6 Environment ACS Paragon Plus
Page 6 of 37
Page 7 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biomacromolecules
ACS Paragon Plus Environment
Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
ACS Paragon Plus Environment
Page 8 of 37
Page 9 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biomacromolecules
151
2.4 Data analysis
152
In order to visualize the band shift caused by mechanical strain at around 1095 cm@ upon tensile
153
stress, Raman spectra were processed and plotted in R 33. Raman spectra were smoothed using
154
the Savitzky-Golay filter with a window length of 9 in the “signal” package 34 before being
155
normalized using the standard normal variate (SNV) method 35 in the range of 1180 to 950 cm@ .
156
For E. regnans wood upon tensile tests, PLS models were built to predict the applied mechanical
157
strain based on Raman spectra in the region between 2300 and 200 cm@ . The 258 spectra from
158
10 E. regnans samples were split randomly into a training group with 180 spectra and a
159
validation group with 78 spectra using the “sample” function in R. Five PLS models were
160
developed using leave-one-out (LOO) cross-validation with the “pls” package 36 based on
161
different spectra pre-processing, including baseline correction, SNV normalization, second
162
derivatives and spectra truncation (1180 ~ 950 cm@ ). Baseline correction was conducted with a
163
modified polynomial algorithm in the “baseline” package 37. Coefficients of determination (R2)
164
and root mean squared errors (RMSE) of the PLS models for the training and the validation
165
group were extracted.
166
Spectral changes caused by mechanical strain were also analyzed by the regression coefficients
167
plot of the PLS model and ‘slope’ spectra. ‘Slope’ spectra were calculated to show the rate of
168
change in spectral intensity at each wavenumber induced by mechanical strain as described
169
previously 14, 38. It was calculated from baseline corrected and normalized spectra. The average
170
‘slope’ spectrum was calculated as the average of 10 ‘slope’ spectra of 10 rewetted E. regnans
171
samples. To quantify the band shift at around 1095 cm@ , a Cauchy-Lorentzian function was used
9 Environment ACS Paragon Plus
Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
172
to fit the spectra in the range from 1103 to 1084 cm@ using the nonlinear least squares method in
173
R. Band positions were extracted from the fitted model.
174 175
3 RESULTS AND DISCUSSIONS
176
3.1 Raman spectra of rewetted wood under tensile strain
177
3.1.1 Visualization of Raman band shifts
178
The signal intensity of Raman bands mainly depends on laser power and its frequency,
179
polarizability and concentration of the chemical group, as well as the sample orientation 39. For
180
biological materials, fluorescence and various physical effects such as surface roughness, particle
181
size and sample inhomogeneity have noticeable influence on the Raman scattering intensity as
182
well 40, 41. Due to these interfering factors, the intensities of the Raman spectra of the rewetted
183
wood fluctuated and it was difficult to identify the small spectral changes caused by mechanical
184
strain from the original spectra (Figure 4A). The fluctuations in the spectra could be removed by
185
SNV normalization in range from 1180 to 950 cm@ , allowing to determine the effect of
186
mechanical tensile strain on the position of the band centered at around 1095 cm@ (Figure 4A
187
inset). The 1095 cm@ band was assigned to CC and CO stretching of the cellulose ring structure
188
17, 42,
189
cellulose backbone. This band shift had been reported for regenerated cellulose as well as single
190
wood and other natural fibers using Raman microscopy 18, 19, 43-45. Our results showed that this
191
tensile strain induced Raman band shift can also be observed in relatively thick (>0.5 mm) E.
192
regnans wood strips, with mechanical behavior similar to solid wood 30. A strong linear
indicating the elongation of cellulose chains and weakening of the covalent bonds in the
10 Environment ACS Paragon Plus
Page 10 of 37
Page 11 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biomacromolecules
193
relationship (r = 0.997) was found between the peak position and the applied strain (Figure 4B).
194
As strain increased from @$ >? to 4649 >?) the band position shifted from 1095.35 to 1093.57
195
cm@ . The band shift rate for this sample was @ :-4 cm@ /%, which was greater compared to the
196
values ;@ :
197
to 3.2.
to @ :. cm@ /%) reported for single wood fibers 19. For a detailed discussion refer
198 199
Figure 4 Unmodified Raman spectra at 12 tensile strain levels of one rewetted E. regnans wood
200
sample (A) and the corresponding band shifts of the 1095 cm@ band (B). Inset (A) shows the
201
band at 1095 cm@ after SNV normalization in the region of 1180 to 950 cm@ . “Strain 0” to
202
“Strain 3” were represented with solid lines, “Strain 4” to “Strain 7” with dashed lines and
203
“Strain 8” to “Strain 11” with dotted lines, each in the color sequence blue, red, pink and green.
204
The band at 1095 cm@ shifted 1.78 cm@ to lower frequencies with a strain increase of 4702 >?:
205
3.1.2 Predicting strain based on Raman spectra using PLS modelling
206
Table 1 shows the effects of spectral processing on the PLS models to predict mechanical strain
207
in wood. Baseline correction aimed to remove the effects of lignin fluorescence, while
208
normalization attempted to compensate multiplicative effects such as laser fluctuation and
11 Environment ACS Paragon Plus
Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 12 of 37
209
uncertainties in focusing 40. Second derivative is another way to reduce baseline effects and
210
truncation to the region from 1180 to 950 cm@ eliminated many variables unaffected by
211
mechanical strain. The root mean square error of cross-validation (RMSECV) ranged from 354.7
212
to 493.1 >?) while the coefficient of determination of cross-validation (R2CV) varied between
213
0.89 and 0.94 (Table 1). Using the models with the validation data resulted in comparable
214
statistics, indicating no over-fitting. Baseline correction followed by normalization resulted in
215
the best performing PLS model (RMSEP = 397.2 >?< using 4 components. The predicted and
216
measured strain for the training group for this PLS model was shown in Figure 5A. The
217
distribution of residuals suggested an unbiased model (Figure 5B). The models based on
218
truncated spectra in the region of 1180 to 950 cm@ and unmodified spectra achieved slightly
219
reduced model statistics with RMSEP of 426.7 and 443.6 >?) respectively. However, more
220
components were required for the unmodified spectra (Table 1). PLS modelling coupled with
221
reflection NIR spectroscopy was successfully used to predict bending (R2 V 0.90) and tensile
222
load (R2 = 0.96) on small air-dry wood samples 46, 47. As NIR is sensitive to moisture, this
223
technique is of limited use for wood with high moisture content. Non-destructive reflection NIR
224
measurements of mechanical strain in green wood were found to be less accurate (R2 = 0.36),
225
due to the high moisture content interfering with structural information and a low signal to noise
226
ratio 7.
227
Table 1 Metrics of PLS models utilizing unmodified and processed Raman spectra to predict
228
mechanical strain in rewetted E. regnans wood. The spectral region was 2300 to 200 cm@ , except
229
for the truncated spectra, for which ranged from 1180 to 950 cm@ .
Pre-processing method
Training group Number of 2 components R CV RMSECV
12 Environment ACS Paragon Plus
Validation group R2P RMSEP
Page 13 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biomacromolecules
Raw data Baseline correction Baseline correction + SNV Second derivative Truncation + SNV
0.91 0.91 0.94 0.89 0.93
440.5 432.8 354.7 493.1 385.9
6 4 4 3 3
0.92 0.92 0.94 0.90 0.93
443.6 456.1 397.2 510.5 426.7
230 231
SNV: standard normal variate (SNV) normalization; CV: Leave-one-out cross-validation;
232
RMSECV: root mean squared error of cross-validation using the training data; RMSEP: root mean
233
squared error of prediction for the validation data; R2: coefficient of determination.
234 235
Figure 5 Measured versus predicted tensile strain of rewetted E. regnans wood (A) and distribution
236
of the residuals (B) for the PLS model based on baseline corrected and SNV normalized Raman
237
spectra of the training data.
238
3.1.3 Spectral changes analyzed by regression coefficients and the average ‘slope’ spectrum
239
Both, regression coefficients of the PLS model and ‘slope’ spectra reflect spectral changes
240
caused by mechanical strain, but were calculated with different algorithms. The regression
241
coefficients of the PLS model represent the change in the response variable (i.e. tensile strain)
13 Environment ACS Paragon Plus
Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
242
induced by one unit change in the predictor variables (i.e. the Raman intensities at each
243
wavenumber). As shown in Figure 6, the PLS regression coefficients and the average ‘slope’
244
spectrum showed almost identical features, with spectral changes in the regions of 1500 to 850
245
cm@ and 500 to 280 cm@ . Split peaks at 1420, 1120, 1095, 895, 456 cm@ with negative peaks
246
immediately followed by positive peaks indicated band shifts to lower frequencies 38, 48.
247
The most dominant Raman band shift induced by mechanical strain at around 1095 cm@ (Figure
248
4 and Figure 6) was reported for various (ligno)cellulosic materials 18-20. The band at 1127 cm@
249
was also previously reported to shift to lower wavenumbers for wood fibers under tensile strain
250
19.
251
deformation of the cellulose skeleton 5. Additionally, in regenerated cellulose fibers the band at
252
895 cm@ was also reported to shift upon stress, but with a lower sensitivity to strain than the
253
band at 1095 cm@
254
HCC at C6 of cellulose 9, 42, indicating that HCO and HCC at C6 were deformed upon tensile
255
strain. NIR spectroscopy observed the deformation of 2OH, involved in intramolecular hydrogen
256
bonds to O6, upon stretching 14, likely affecting the other covalent bonds connected to this group.
257
The region from 1430 to 1350 cm@ was assigned to the bending vibrations of COH groups 42.
258
Therefore, the Raman band shift at 1420 cm@ suggested the molecular deformation of COH
259
groups, which was consistent with MIR and NIR evidence of tensile strain induced deformation
260
of OH groups involved in hydrogen bonds 14, 49-51. Consistent with our observations, Raman band
261
shifts at 1414, 1095 and 895 cm@ were reported for high-modulus cellulose fibers under tensile
262
stress 52 and additionally at around 459 cm@ in single ramie fibers during tensile tests 53. The
263
region between 250 and 550 cm@ was mainly attributed to the skeletal bending motions of
264
cellulose, including CCC, COC, OCC and OCO bonds 42, 54.
Both bands were assigned to the stretching of cellulose CC and CO bonds, confirming the
17.
The band at 895 cm@ was tentatively assigned to the bending of HCO and
14 Environment ACS Paragon Plus
Page 14 of 37
Page 15 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biomacromolecules
265 266
Figure 6 Regression coefficients (grey) of the PLS model based on baseline corrected and
267
normalized Raman spectra and the average ‘slope’ spectrum (black) for 10 rewetted E. regnans
268
wood samples upon stretching in the region from 2300 to 200 cm@ . The dashed line shows a
269
representative unmodified Raman spectrum of rewetted E. regnans wood.
270
3.2 Quantification of Raman band shifts caused by stretching
271
Apart from PLS modelling, another approach to predict strain is to use linear models based on
272
the positions of the Raman band at around 1095 cm@ . Raman band positions at around 1095
273
cm@ were determined by fitting a Cauchy-Lorentzian function to the signal. The dependence of
274
the Raman band position on tensile strain was quantified for rewetted E. regnans (Figure 7A)
275
and green E. quadrangulata (Figure 7B) wood (Figure 7). As tensile strain increased from the
276
relaxed state to the maximal applied strain level, the band at 1095 cm@ shifted to lower
277
wavenumbers for both, previously dried water-saturated and never-dried green wood, indicating
278
elongation of the cellulose ring structure. With the release of the applied tensile stress, the
15 Environment ACS Paragon Plus
Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
279
Raman band returned to its original position. No significant differences were found between the
280
band positions before (“Strain 0”) and after the tensile tests (RL) (p = 0.31 and 0.27 for rewetted
281
and green wood, respectively). Similar behavior has been reported for single wood fibers 19 but
282
no similar studies have been done for bulk wood. This study showed that the Raman band at
283
around 1095 cm@ was also sensitive to strain for wood samples with a thickness of 0.5 to 1 mm.
284
Thin wood samples, only a few cells thick, behave differently from bulk wood during tensile
285
tests due to unrestricted cellular and supramolecular deformations caused by removing the
286
supporting structure of neighboring tissue 29. With eucalyptus fibers having diameters of 10 to 20
287
µm 55, this effect was minimized here, as the wood strips used in this study were comprised of
288
more than 25 cells in thickness. For even larger wood samples MOE was shown to increase with
289
sample thickness 28, 30.
290
16 Environment ACS Paragon Plus
Page 16 of 37
Page 17 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biomacromolecules
291
Figure 7 Dependence of the position of the Raman band at around 1095 cm@ on tensile strain for
292
rewetted E. regnans (A) wood (n = 10) and green E. quadrangulata (B) wood (n = 8). Raman
293
spectra were collected as mechanical strain increased stepwise from “Strain 0” (approx. 0 >?< to
294
“Strain 11” (approx. 4500 >?< with a step of approx. 410 >? and then returned to the relaxed (RL)
295
state.
296
The mean Raman band position in the relaxed state for rewetted E. regnans and green E.
297
quadrangulata was 1095.41 (±0.15) and 1095.26 (±0.17) cm@ , respectively (Table 2). The
298
difference in band position between E. regnans and E. quadrangulata was not statistically
299
significant (p = 0.08). The observed variation could be caused by instrumentation like laser
300
fluctuation or variation in chemical composition and wood anatomy features such as microfibril
301
angle and grain angle 56. Additionally, the drying history was reported to affect the cellulose
302
structure 57, 58.
303
The band shift rate with respect to strain for rewetted E. regnans and green E. quadrangulata
304
during the tensile tests was @ :,. and @ : $ cm@ /% (Table 2). The significantly lower (p ?:
cm@ /% for air-dried hemp fibers 20 and @ :
Species
Sample condition
Sample size
E. regnans
rewetted
10
E. quadrangulata
green
8
to @ :. cm@ /% for air-dried single wood
Band position (cm@ ) 1095.41 (0.15) 1095.26 (0.17)
Band shift rate (cm@ /%)
R2
@ :,. (0.48)
0.95 (0.02)
@ : $ (0.19)
0.89 (0.05)
327 328
3.3 Monitoring the release of growth stress using Raman spectroscopy
329
Wood formed in the upper side of leaning hardwood trees, referred to as tension wood, is usually
330
associated with high tensile stress 62. Correspondingly, growth strain ranged from 1294 to 1446
18 Environment ACS Paragon Plus
Page 19 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biomacromolecules
331
>? in the tension sides and from @ -3 to 64 >? in the opposite sides of the e investigated E.
332
bosistoana stems (Figure 8). Raman spectra were taken from these stems before and after
333
releasing growth stress. For the tension sides, the band position shifted to higher wavenumbers (p
334
< 0.016) after cutting the grooves, suggesting the release of tensile strain (Figure 8). In contrast,
335
no significant differences were found for the opposite sides (p > 0.140), attributed to the
336
relatively low stress levels in opposite wood and the detection limit of the Raman method.
337 338
Figure 8 Raman band positions of three leaning E. bosistoana trees before (open boxes) and after
339
(filled boxes) releasing growth stress. Positive strain values stand for tensile strain (contraction)
340
and negative values for compressive strain (expansion). “Side T” represents tension wood and
341
“Side O” opposite wood.
342
The variation in the Raman band positions was greater between the samples than within the
343
samples. The reasons are not fully understood, but could be related to the variation of chemical
344
composition 63 and anatomical structures 55 at the millimeter scale. The variation between
345
samples is a challenge for the prediction of growth strain based on Raman band positions.
346
3.4 Correlation of growth strain to the position of the 1095 cm@ Raman band
19 Environment ACS Paragon Plus
Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
347
Methods to measure growth strain in trees, including the Nicholson technique 24, the CIRAD-
348
Forêt method 25 and the strain gauge method, had been reviewed by 64. The splitting test
349
proposed by 26 estimates growth strain of the stem based on the opening at the ends after splitting
350
the sample along the pith. This estimation will be affected by the three-dimensional distribution
351
of growth strain throughout the whole sample. In contrast, the other methods are a localized
352
measurement of surface growth strain. Although linear correlations were found between the
353
methods 26, 27, it is still challenging to build quantitative relationships between the different
354
methods 22.
355
As shown in Figure 9A, a strong linear relationship (r = 0.78) was found between the average
356
growth strain measured by strain gauges and the EGS using the splitting method. But the growth
357
strain measured by strain gauges was significantly lower compared to the EGS using the splitting
358
method. Additionally, the growth strain measured with the strain gauges at each side of the same
359
stem differed by up to more than four times (data not shown), confirming asymmetric SGS
360
around the stem 24.
361
The positions of the Raman band at around 1095 cm@ were moderately correlated (r = 0.48) with
362
growth strain measured using strain gauges (Figure 9B), but only a lower correlation (r = 0.27)
363
was found between the Raman band positions and the EGS (data not shown). This agreed with
364
the understanding that the Raman measurements reflected the localized molecular strain in the
365
stem and that this strain is variable throughout the stem 24.
366
The intercept of the linear fit was 1095.24 cm@ , consistent with the band position in the relaxed
367
state for green E. quadrangulata wood strips (Table 2). However, the slope of the linear fit was
368
@-: cm@ /%, more than twice the band shift rate of mechanically deformed samples (Table 2).
20 Environment ACS Paragon Plus
Page 20 of 37
Page 21 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biomacromolecules
369
This suggested that the growth strain measured by strain gauges were underestimated, while the
370
band shift rate based on EGS calculated from the splitting test was @ :$ cm@ /%, similar to
371
those obtained by mechanical strain (Table 2). Similar discrepancies were reported between the
372
CIRAD-Forêt method and strain gauge measurements (Yoshida & Okuyama, 2002). The groove
373
depth and the distance to the strain gauge can affect the measured growth strain 27. The
374
surrounding wood tissue will impede the release of growth stresses 21. Additionally, according to
375
an X-ray diffraction study (Toba et al., 2013), molecular growth stress remained in a wood block
376
after being removed from the living stem. It is likely that not all growth stresses were released in
377
this experiment.
378 379
Figure 9 Relationship between growth strain measured by strain gauges and the estimated growth
380
strain (EGS) by the splitting method (A) and the dependence of the Raman band at around 1095
21 Environment ACS Paragon Plus
Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
381
cm@ on growth strain measured by strain gauges (B). Each band position was calculated as the
382
average of 9 Raman spectra collected on the wood surface at the location of the strain gauge.
383 384
While this study has shown that molecular strain causing surface growth strain in stems can be
385
detected with Raman spectroscopy, technological advances are needed for Raman spectroscopy
386
to become commercially feasible for non-destructive segregation of eucalyptus logs for solid
387
wood processing. The relatively low correlation (r = 0.48) in the present study implied that
388
estimations of growth strain by Raman spectroscopy will result in falsely rejected usable logs.
389
This might be acceptable for a pulp chip business, aiming to increase revenue by identifying top
390
logs suitable for solid wood production, but is likely to be uneconomic for high value solid wood
391
plantations, which will lose revenue by falsely rejecting quality logs. Technological challenges
392
like measurement time and the focusing procedure of the dispersive system renders industrial
393
applications inefficient. Furthermore, assuming a band shift rate of @ cm@ /%, 1000 >? only
394
induces a band shift of 0.3 cm@ , making it necessary to address factors like instrumental
395
uncertainties, sample surface conditions, fluorescence, as well as variations in wood properties
396
and growth strain 24, 25, 65. Better prediction of growth strain can be expected with a Raman
397
instrument that can cover a larger surface area and averaging out local variations in chemical
398
composition and growth strain.
399 400 401
4 CONCLUSION
22 Environment ACS Paragon Plus
Page 22 of 37
Page 23 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biomacromolecules
402
A partial least squares model based on pre-processed spectra was able to predict mechanical
403
strain with a root mean square error of 397.2 >? and R2 of 0.94. The major spectral changes
404
caused by mechanical strain were in the regions of 1500 to 850 cm@ and 500 to 280 cm@ . Apart
405
from the predominant band shift at 1095 cm@ , detailed analysis suggested that mechanical strain
406
also caused Raman band shifts at 1420, 1120, 895, 456 cm@ . The spectral changes could be
407
interpreted as the molecular deformation of the cellulose skeleton, HCO and HCC groups at C6
408
as well as COH structures. Further, studies with green wood indicated that the Raman band shift
409
rate with respect to strain was lower in never-dried green wood than previously dried rewetted
410
wood. The release of growth strain was successfully detected by Raman spectroscopy based on
411
the band positions at around 1095 cm@ . A moderate correlation was found between the band
412
positions at around 1095 cm@ and the growth strain measured by strain gauges, indicating that
413
the technology can potentially be used to segregate logs for solid wood production in certain
414
circumstances. However, improvements in Raman instrumentation, allowing more precise
415
measurements and spatial averaging, are needed.
416 417
AUTHOR INFORMATION
418
Corresponding Author
419
*New Zealand School of Forestry, University of Canterbury, Private Bag 4800, Christchurch,
420
New Zealand.
421
Tel: +64 33695949. Email:
[email protected].
422
Author Contributions
23 Environment ACS Paragon Plus
Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
423
The manuscript was written through contributions of both authors. FG conducted the
424
experiments and drafted the manuscript. CA provided guidance to the experimental design and
425
manuscript revision. All authors have given approval to the final version of the manuscript.
426
ACKNOWLEDGMENT
427
This work was financially supported by the Ministry of Business, Innovation and Employment,
428
New Zealand, through the Specialty Wood Products Partnership (FFRX1501). Nigel Pink and
429
Gert Hendriks (University of Canterbury) helped with the preparation of the wood samples.
430 431
REFERENCES
432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458
1. Larkin, P., Infrared and raman spectroscopy: principles and spectral interpretation. Elsevier: 2011. 2. Amer, M., Raman spectroscopy for soft matter applications. John Wiley & Sons: 2009. 3. Vankeirsbilck, T.; Vercauteren, A.; Baeyens, W.; Van der Weken, G.; Verpoort, F.; Vergote, G.; Remon, J. P., Applications of Raman spectroscopy in pharmaceutical analysis. TrAC, Trends Anal. Chem. 2002, 21 (12), 869-877. 4. Richter, S.; Müssig, J.; Gierlinger, N., Functional plant cell wall design revealed by the Raman imaging approach. Planta 2011, 233 (4), 763-772. 5. Gierlinger, N.; Schwanninger, M., Chemical imaging of poplar wood cell walls by confocal Raman microscopy. Plant Physiol. 2006, 140 (4), 1246-1254. 6. Carron, K.; Cox, R., Qualitative analysis and the answer box: A perspective on portable Raman spectroscopy. Anal. Chem. 2010, 82 (9), 3419-3425. 7. Guo, F. Molecular deformation of wood and cellulose studied by near infrared and Raman spectroscopy. University of Canterbury, Christchurch, New Zealand, 2019. 8. Agarwal, U. P., 1064 nm FT-Raman spectroscopy for investigations of plant cell walls and other biomass materials. Front. Plant Sci. 2014, 5 (490). 9. Agarwal, U. P., An overview of Raman spectroscopy as applied to lignocellulosic materials. In Advances in lignocellulosics characterization, Argyropoulos, D. S., Ed. TAPPI Press: 1999; pp 201-225. 10. Peura, M.; Kölln, K.; Grotkopp, I.; Saranpää, P.; Müller, M.; Serimaa, R., The effect of axial strain on crystalline cellulose in Norway spruce. Wood Sci. Technol. 2007, 41 (7), 565-583. 11. Zabler, S.; Paris, O.; Burgert, I.; Fratzl, P., Moisture changes in the plant cell wall force cellulose crystallites to deform. J. Struct. Biol. 2010, 171 (2), 133-141. 12. Åkerholm, M.; Salmén, L., Interactions between wood polymers studied by dynamic FTIR spectroscopy. Polymer 2001, 42 (3), 963-969. 13. Salmén, L.; Bergström, E., Cellulose structural arrangement in relation to spectral changes in tensile loading FTIR. Cellulose 2009, 16 (6), 975-982.
24 Environment ACS Paragon Plus
Page 24 of 37
Page 25 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504
Biomacromolecules
14. Guo, F.; Altaner, C. M., Molecular deformation of wood and cellulose studied by near infrared spectroscopy. Carbohydr. Polym. 2018, 197, 1-8. 15. Eichhorn, S. J.; Young, R. J., The Young's modulus of a microcrystalline cellulose. Cellulose 2001, 8 (3), 197-207. 16. Mitra, V. K.; Risen Jr, W. M.; Baughman, R. H., A laser Raman study of the stress dependence of vibrational frequencies of a monocrystalline polydiacetylene. J. Chem. Phys. 1977, 66 (6), 2731-2736. 17. Hamad, W. Y.; Eichhorn, S., Deformation micromechanics of regenerated cellulose fibers using Raman spectroscopy. J. Eng. Mater. Technol. 1997, 119 (3), 309-313. 18. Eichhorn, S. J.; Sirichaisit, J.; Young, R. J., Deformation mechanisms in cellulose fibres, paper and wood. J. Mater. Sci. 2001, 36 (13), 3129-3135. 19. Gierlinger, N.; Schwanninger, M.; Reinecke, A.; Burgert, I., Molecular changes during tensile deformation of single wood fibers followed by Raman microscopy. Biomacromolecules 2006, 7 (7), 2077-2081. 20. Peetla, P.; Schenzel, K. C.; Diepenbrock, W., Determination of mechanical strength properties of hemp fibers using near-infrared Fourier transform Raman microspectroscopy. Appl. Spectrosc. 2006, 60 (6), 682-691. 21. Alméras, T.; Clair, B., Critical review on the mechanisms of maturation stress generation in trees. J. R. Soc. Interface 2016, 13 (122), 20160550. 22. Yang, J. L.; Waugh, G., Growth stress, its measurement and effects. Aust. For. 2001, 64 (2), 127-135. 23. Clair, B.; Alméras, T.; Pilate, G.; Jullien, D.; Sugiyama, J.; Riekel, C., Maturation stress generation in poplar tension wood studied by synchrotron radiation microdiffraction. Plant Physiol. 2011, 155 (1), 562-570. 24. Nicholson, J. E., A rapid method for estimating longitudinal growth stresses in logs. Wood Sci. Technol. 1971, 5 (1), 40-48. 25. Gerard, J.; Bailleres, H.; Fournier, M.; Thibaut, B., Wood quality in plantation Eucalyptus - a study of variation in three reference properties. Bois For. Trop. 1995, 245, 101110. 26. Chauhan, S. S.; Entwistle, K., Measurement of surface growth stress in Eucalyptus nitens Maiden by splitting a log along its axis. Holzforschung 2010, 64 (2), 267-272. 27. Yoshida, M.; Okuyama, T., Techniques for measuring growth stress on the xylem surface using strain and dial gauges. Holzforschung 2002, 56 (5), 461. 28. Buchelt, B.; Pfriem, A., Influence of wood specimen thickness on its mechanical properties by tensile testing: solid wood versus veneer. Holzforschung 2011, 65 (2), 249-252. 29. Navi, P.; Rastogi, P. K.; Gresse, V.; Tolou, A., Micromechanics of wood subjected to axial tension. Wood Sci. Technol. 1995, 29 (6), 411-429. 30. Yu, Y.; Jiang, Z.-h.; Tian, G.-l., Size effect on longitudinal MOE of microtomed wood sections and relevant theoretical explanation. For. Stud. China 2009, 11 (4), 243-248. 31. Conti, C.; Botteon, A.; Bertasa, M.; Colombo, C.; Realini, M.; Sali, D., Portable sequentially shifted excitation Raman spectroscopy as an innovative tool for in situ chemical interrogation of painted surfaces. Analyst 2016, 141 (15), 4599-4607. 32. Davies, N. T.; Apiolaza, L. A.; Sharma, M., Heritability of growth strain in Eucalyptus bosistoana: a Bayesian approach with left-censored data. N. Z. J. For. Sci. 2017, 47, 5. 33. R Core Team, A language and environment for statistical computing. In R foundation for statistical computing, Vienna, Austria, 2017.
25 Environment ACS Paragon Plus
Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549
34. Signal developers, Signal: Signal processing. 2013. 35. Stevens, A.; Ramirez-Lopez, L., An introduction to the prospectr package. In R package version 0.1.3, 2013. 36. Mevik, B.-H.; Wehrens, R.; Liland, K. H., Pls: Partial least squares and principal component regression. In R package version 2.6-0, 2016. 37. Liland, K. H.; Mevik, B.-H., Baseline: Baseline xorrection of spectra. In R package version 1.2-1, 2015. 38. Altaner, C. M.; Horikawa, Y.; Sugiyama, J.; Jarvis, M. C., Cellulose (c investigated by IR-spectroscopy at low temperatures. Cellulose 2014, 21 (5), 3171-3179. 39. Smith, E.; Dent, G., Modern Raman spectroscopy: a practical approach. John Wiley & Sons, Ltd: 2005. 40. Liland, K. H.; Kohler, A.; Afseth, N. K., Model-based pre-processing in Raman spectroscopy of biological samples. J. Raman Spectrosc. 2016, 47 (6), 643-650. 41. Huang, J.; Romero-Torres, S.; Moshgbar, M., Practical considerations in data pretreatment for NIR and Raman spectroscopy. American Pharmaceutical Review 2010, 13 (6), 116. 42. Wiley, J. H.; Atalla, R. H., Band assignments in the raman spectra of celluloses. Carbohydr. Res. 1987, 160 (Supplement C), 113-129. 43. Hsieh, Y.-C.; Yano, H.; Nogi, M.; Eichhorn, S. J., An estimation of the Young’s modulus of bacterial cellulose filaments. Cellulose 2008, 15 (4), 507-513. 44. Šturcová, A.; Davies, G. R.; Eichhorn, S. J., Elastic modulus and stress-transfer properties of tunicate cellulose whiskers. Biomacromolecules 2005, 6 (2), 1055-1061. 45. Eichhorn, S. J.; Young, R. J.; Davies, G. R., Modeling crystal and molecular deformation in regenerated cellulose fibers. Biomacromolecules 2005, 6 (1), 507-513. 46. André, N.; Labbé, N.; Rials, T. G.; Kelley, S. S., Assessment of wood load condition by Near Infrared (NIR) spectroscopy. J. Mater. Sci. 2006, 41 (7), 1879-1886. 47. Sandak, J.; Sandak, A.; Pauliny, D.; Krasnoshlyk, V.; Hagman, O., Near infrared spectroscopy as a tool for estimation of mechanical stresses in wood. Adv. Mater. Res. (DurntenZurich, Switz.) 2013, 778, 448-453. 48. Åkerholm, M.; Hinterstoisser, B.; Salmén, L., Characterization of the crystalline structure of cellulose using static and dynamic FT-IR spectroscopy. Carbohydr. Res. 2004, 339 (3), 569-578. 49. Hinterstoisser, B.; Åkerholm, M.; Salmén, L., Load distribution in native cellulose. Biomacromolecules 2003, 4 (5), 1232-1237. 50. Altaner, C. M.; Thomas, L. H.; Fernandes, A. N.; Jarvis, M. C., How cellulose stretches: synergism between covalent and hydrogen bonding. Biomacromolecules 2014, 15 (3), 791-798. 51. Kong, K.; Eichhorn, S. J., The influence of hydrogen bonding on the deformation micromechanics of cellulose fibers. J. Macromol. Sci., Part B: Phys. 2005, 44 (6), 1123-1136. 52. Eichhorn, S. J.; Young, R. J.; Davies, R. J.; Riekel, C., Characterisation of the microstructure and deformation of high modulus cellulose fibres. Polymer 2003, 44 (19), 59015908. 53. Gierlinger, N., New insights into plant cell walls by vibrational microspectroscopy. Appl. Spectrosc. Rev. 2018, 53 (7), 517-551. 54. Österberg, M.; Schmidt, U.; Jääskeläinen, A.-S., Combining confocal Raman spectroscopy and atomic force microscopy to study wood extractives on cellulose surfaces. Colloids and Surfaces A: Physicochemical and Engineering Aspects 2006, 291 (1), 197-201.
26 Environment ACS Paragon Plus
Page 26 of 37
Page 27 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577
Biomacromolecules
55. Wilkes, J., Variations in wood anatomy within species of Eucalyptus. IAWA J 1988, 9 (1), 13-23. 56. Agarwal, U. P., Raman imaging to investigate ultrastructure and composition of plant cell walls: distribution of lignin and cellulose in black spruce wood (Picea mariana). Planta 2006, 224 (5), 1141. 57. Hill, S. J.; Kirby, N. M.; Mudie, S. T.; Hawley, A. M.; Ingham, B.; Franich, R. A.; Newman, R. H., Effect of drying and rewetting of wood on cellulose molecular packing. Holzforschung 2010, 64 (4), 421-427. 58. Toba, K.; Yamamoto, H.; Yoshida, M., Mechanical interaction between cellulose microfibrils and matrix substances in wood cell walls induced by repeated wet-and-dry treatment. Cellulose 2012, 19 (4), 1405-1412. 59. Agarwal, U. P.; Ralph, S. A.; Reiner, R. S.; Baez, C., Probing crystallinity of neverdried wood cellulose with Raman spectroscopy. Cellulose 2016, 23 (1), 125-144. 60. Keckes, J.; Burgert, I.; Fruhmann, K.; Muller, M.; Kolln, K.; Hamilton, M.; Burghammer, M.; Roth, S. V.; Stanzl-Tschegg, S.; Fratzl, P., Cell-wall recovery after irreversible deformation of wood. Nat. Mater. 2003, 2 (12), 810-813. 61. Green, D. W.; Winandy, J. E.; Kretschmann, D. E., Mechanical properties of wood. In Wood handbook: wood as an engineering material, Ross, R. J., Ed. Department of Agriculture, Forest Service, Forest Products Laboratory: Madison, Wisconsin, 2010. 62. Gardiner, B.; Barnett, J.; Saranpää, P.; Gril, J., The biology of reaction wood. Springer: 2014. 63. Downes, G. M.; Touza, M.; Harwood, C.; Wentzel-Vietheer, M., NIR detection of nonrecoverable collapse in sawn boards of Eucalyptus globulus. Eur. J. Wood Wood Prod. 2014, 72 (5), 563-570. 64. Yang, J. L.; Baillères, H.; Okuyama, T.; Muneri, A.; Downes, G., Measurement methods for longitudinal surface strain in trees: a review. Aust. For. 2005, 68 (1), 34-43. 65. Yang, J. L.; Fife, D.; Matheson, A. C., Growth strain in three provenances of plantationgrown Eucalyptus globulus Labill. Aust. For. 2001, 64 (4), 248-256.
578
27 Environment ACS Paragon Plus
Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
159x80mm (600 x 600 DPI)
ACS Paragon Plus Environment
Page 28 of 37
Page 29 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biomacromolecules
Figure 1. Wood sample used for tensile tests while collecting Raman spectra. The grain (longitudinal direction) runs along the longest dimension of the sample. Dimensions are labelled in millimeters. 83x37mm (600 x 600 DPI)
ACS Paragon Plus Environment
Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Figure 2. Schematic drawing of the experimental setup for monitoring the release of growth strain by Raman spectroscopy on green stems. The dashed red lines indicate the grooves cut to release growth stress. T: tension wood side; O: opposite wood side. 85x37mm (600 x 600 DPI)
ACS Paragon Plus Environment
Page 30 of 37
Page 31 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biomacromolecules
Figure 3. Illustration of evaluating growth strain of green stems using Raman spectroscopy. In total, 18 spectra were collected from each sample. The dashed red line indicates the cut for the splitting test. 85x34mm (600 x 600 DPI)
ACS Paragon Plus Environment
Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Figure 4 Unmodified Raman spectra at 12 tensile strain levels of one rewetted E. regnans wood sample (A) and the corresponding band shifts of the 1095 cm−1 band (B). Inset (A) shows the band at 1095 cm−1 after SNV normalisation in the region of 1180 to 950 cm−1. “Strain 0” to “Strain 3” were represented with solid lines, “Strain 4” to “Strain 7” with dashed lines and “Strain 8” to “Strain 11” with dotted lines, each in the color sequence blue, red, pink and green. The band at 1095 cm−1 shifted 1.78 cm−1 to lower frequencies with a strain increase of 4702 µε. 169x67mm (600 x 600 DPI)
ACS Paragon Plus Environment
Page 32 of 37
Page 33 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biomacromolecules
Figure 5 Measured versus predicted tensile strain of rewetted E. regnans wood (A) and distribution of the residuals (B) for the PLS model based on baseline corrected and SNV normalized Raman spectra of the training data 219x92mm (600 x 600 DPI)
ACS Paragon Plus Environment
Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Figure 6 Regression coefficients (grey) of the PLS model based on baseline corrected and normalized Raman spectra and the average ‘slope’ spectrum (black) for 10 strained rewetted E. regnans wood samples in the region from 2300 to 200 cm−1. The dashed line shows a representative unmodified Raman spectrum of rewetted E. regnans wood. 85x85mm (600 x 600 DPI)
ACS Paragon Plus Environment
Page 34 of 37
Page 35 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biomacromolecules
Figure 7 Dependence of the position of the Raman band at around 1095 cm−1 on tensile strain for rewetted E. regnans (A) wood (n = 10) and green E. quadrangulata (B) wood (n = 8). Raman spectra were collected as mechanical strain increased stepwise from “Strain 0” (approx. 0 µε) to “Strain 11” (approx. 4500 µε) with a step of approx. 410 µε and then returned to the relaxed (RL) state. 85x103mm (600 x 600 DPI)
ACS Paragon Plus Environment
Biomacromolecules 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Figure 8 Raman band positions of three leaning E. bosistoana trees before (open boxes) and after (filled boxes) releasing growth stress. Positive strain values stand for tensile strain (contraction) and negative values for compressive strain (expansion). “Side T” represents tension wood and “Side O” opposite wood. 169x66mm (600 x 600 DPI)
ACS Paragon Plus Environment
Page 36 of 37
Page 37 of 37 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biomacromolecules
Figure 9 Relationship between growth strain measured by strain gauges and the estimated growth strain (EGS) by the splitting method (A) and the dependence of the Raman band at around 1095 cm−1 on growth strain measured by strain gauges (B). Each band position was calculated as the average of 9 Raman spectra collected on the wood surface at the location of the strain gauge. 85x103mm (600 x 600 DPI)
ACS Paragon Plus Environment