An investigation of the fluid-holding cavities in a lignocellulose-based bamboo matrix via a combined X-ray microtomography and proton time-domain NMR approach

Bamboo (Dendrocalamus giganteus) is a functionally graded material with well-organized hierarchical structures. Its micrometer-sized vascular bundles and parenchymatic living cells allow an efficient upward flow of water and nutrients, endowing the organism with remarkably fast growth. As demonstrated recently, the hollow microstructure channels can be explored as a natural template for microfluidics applications in chemical synthesis, analytical detection, solar steam generation, and electrochemical devices. Thus, the knowledge of the kinetics of the imbibition and spatial distribution of fluid through the microcavities of the bamboo vegetal tissue became of interest. Here, we employed a combination of X-ray microtomography (µCT) and proton time-domain NMR (TD-NMR) to identify, measure, and investigate empty volumes embedded in the bamboo’s tissue as experienced by different organic and inorganic fluids, namely dimethyl sulfoxide (DMSO) and distilled water (H2O). Results have shown that the extensive communication between the voids (parenchymal cells and vascular channels) does not prevent the individuation of fluid reservoirs with different behaviors, as evidenced by the H2O desorption profile. Bamboo parenchyma was the tissue that retained more residual H2O after desorption. The difference in wettability of bamboo tissues by DMSO and H2O was detected.


Introduction
The bamboo (Dendrocalamus giganteus) plant tissue matrix is chemically and morphologically complex. It is considered an intelligent functionally graded material with oriented, hierarchical, and shaped biological tissues that developed over long-term evolution and selection (Nogata and Takahashi 1995). At the molecular level, it consists mainly of cellulose, a linear polysaccharide, highly organized due to hydrogen bonds between the abundant hydroxyl moieties of β-(1-4)-d-glucopyranose monomers; hemicellulose, a branched polysaccharide of low crystallinity consisting of monosaccharide units of 5 to 6 carbons (mainly glucose, xylose, galactose, and mannose) and lignin, a heterogeneous phenolic macromolecule (Youssefian and Rahbar 2015). Both cellulose and lignin are responsible for the high mechanical properties displayed by the plant and are found in large quantities (55% and 25%, respectively) in the bamboo matrix (Ji et al. 2013;Youssefian and Rahbar 2015).
At the micrometric level, the bamboo tissues concentrate cavities where liquids can flow and be stored, being greatly responsible for its fast growth. In the longitudinal direction of the bamboo's culm, straight and parallel channels with diameters varying from a few tens to a few hundred micrometers compose the phloem, protoxylem, and metaxylem vascular systems. These vascular bundles are mechanically protected by the sclerenchymatous tissue, mainly constituted of highly oriented, crystalline, and hierarchical ligno-carbohydrate complex (LCC) (Hao et al. 2018). Also present in the culm, polyhedral-shaped micrometric cells, with highly lignified walls interconnected via micro-and nanostructures called pits, constitute the parenchyma tissue (Xu et al. 2011;Ji et al. 2013;Peng et al. 2013;Youssefian and Rahbar 2015;Ginoble Pandoli et al. 2019).
Functional gradients, heterogeneities, and highly anisotropic morphology of biological materials have recently attracted much attention in investigations concerning novel sustainable and high-performance materials (Liu et al. 2017). Bamboo, less than wood, has been explored for new challenging applications, adding mechanical (Jin et al. 2015a, b;Li et al. 2017a), optical (Wang et al. 2022), biological (Pandoli et al. 2016) and chemical properties (Palma et al. 2020;De Sá et al. 2022). In the literature reports concerning bamboo specimens, most chemical modifications were conducted on the external surface of the material. However, few reports have been published exploring the internal structure of bamboo as biotemplates for microfluidic chemical platforms (Kuan et al. 2015;De Sa et al. 2019;Li et al. 2022), steam generation (Sheng et al. 2020;Chen et al. 2021), and electrochemical devices . These works addressed different ways of taking advantage of the internal structure of bamboo tissues, describing the filling and chemical functionalization of the parenchyma and/or vascular bundles with transition metal salts or nanoparticles, such as copper and silver (Ginoble Pandoli et al. 2022Pandoli et al. , 2023. Recently, the anisotropic behavior of water vapor sorption into the hierarchical structure of bamboo was demonstrated, revealing different sorption directions, longitudinal (L) and transversal (T), for a saturated salt solution and water vapor (Chen et al. 2021). Thus, the need for a deeper understanding of how vascular channels and cell tissues interact with different fluids mainly used in organic chemistry has become evident.
X-ray microtomography (µCT) is a non-destructive technique widely used for the characterization of the microstructure of plants and has been used for bamboo as well (Palombini et al. 2020). Time-domain nuclear magnetic resonance (TD-NMR) is an efficient, fast, and non-destructive way to probe fluid behavior in confined environments (Steele et al. 2016;Singer et al. 2018). It is widely used in the study of oil reservoir rock matrices (Muhammad et al. 2014;Plastino et al. 2017;Rudszuck et al. 2019), water trapped in hydrogels (Rodrigues et al. 2016(Rodrigues et al. , 2017Ratzsch et al. 2017) and to study the molecular dynamics of liquids during the hardening of cementitious materials (Valckenborg et al. 2001;Jansen et al. 2018).
TD-NMR has seen some uses in the study of liquids confined in vegetal tissue, including quantification of oil in seeds (Prestes et al. 2007;Colnago et al. 2011Colnago et al. , 2014Li and Ma 2021;. Zhang et al. (2016) used TD-NMR methods to characterize hierarchical pore structures in cotton fibers by measuring the different relaxation behaviors of the liquids imbibed in the cellulose matrix. The use of single-sided TD-NMR equipment is a very practical approach to study liquids in plant tissues (Menegazzo et al. 2020(Menegazzo et al. , 2021. Mohebbi et al. (2019) used the technique to quantify the amount of water in open foams and fibrous materials like paper. Li et al. (2017b) determined the rate of humidity desorption in wood (Populus platyphylla and Pinus sylvestris) samples via TD-NMR. Mikac et al. (2021) used magnetic resonance imaging (MRI) and TD-NMR methods to investigate the behavior of moisture content in Fagus sylvatica branches. Water uptake and moisture content in plant tissues were also investigated by Gezici-Koç et al. (2017) and Longo et al. (2023). Recently, Li and Ma (2022) proposed the investigation of the cell wall water of wood in the hygroscopic range using 2D TD-NMR, successfully differentiating the water mobility in cellulosic material. Some other examples of 2D TD-NMR experiments applied to the investigation of water molecules behavior in wood samples were reported by Cox et al. (2010), Bonnet et al. (2017) and Stagno et al. (2021).
The interest in investigating the behavior of fluids in the cavities of harvested and treated bamboo emerges when its use as a renewable platform for carrying out physicochemical processes becomes evident. However, morphological heterogeneities in the cavities that retain and transport fluids inside the living plant can be found in each internode. Therefore, the structure-properties relationship in this context is highly relevant for the current and future applications of this class of renewable porous material.
In this work, we analyzed the behavior of organic (dimethyl sulfoxide) and inorganic (distilled water) fluids in the lignocellulosic matrix of bamboo (D. giganteus) through proton TD-NMR experiments, with the support of qualitative and quantitative structural data extracted from µCT analysis.

Bamboo sampling
Bamboo specimens consisted of cylinders of 1 cm in height and 0.6 cm in diameter, extracted from internode number 3 closer to the ground. The sampling process is illustrated in Fig. 1. Considering the radial direction from the outer to the inner wall of the internode, we cut specimens at the external (close to the bark), intermediate and internal parts. The bamboo samples will be referred to as External for the outermost fraction, Middle for the intermediate fraction, and Internal for the innermost fraction. To eliminate organic residues, samples were washed under magnetic stirring in distilled water at 60 °C for 2 h. The samples were dried in an oven at 100 °C until constant mass.

Probe fluids
Initially, we selected a wide range of organic and inorganic fluids to serve as a relaxation probe in the experiments. However, due to issues such as high volatility (methanol, ethanol, chloroform, tetrahydrofuran), which could prevent some of the lengthier experiments, we finally chose two solvents with high boiling points. Distilled water (H 2 O) (resistivity of 18 MΩ, filtered through a 0.22 μm mesh membrane, Milli-Q Gradient System A10, Millipore, USA) and dimethyl sulfoxide (DMSO) (ACS Reagent, Sigma-Aldrich, St. Louis, USA) were selected due to their low volatility, varying polarity values and different capacity for establishing physicochemical interactions with the lignocellulosic matrix of the bamboo plant.

Probe fluids saturation of the bamboo sample fractions
We conducted the saturation of the bamboo samples (external, middle, and internal) with each probe fluid using a vacuum chamber (Fig. 2). First, the samples were oven-dried at 100 °C until constant mass. Then, they were submerged in a given fluid volume (usually 5 mL of DMSO or H 2 O) inside  a Schlenk flask connected to a vacuum line. At this point, the samples remained afloat in the liquid. The system remained under reduced pressure (760 mmHg) for 2 h. When the vacuum was broken, the samples sank into the fluid, suggesting extensive filling of their structural cavities. To avoid interfering with the water content inside the bamboo sample's tissues we only surface-dried them with absorbent paper, weighed them again, and taken to the proton time-domain nuclear magnetic resonance equipment for analysis. The samples were taken to an oven at 100 °C for drying to constant mass before changing the type of probe fluid. All gravimetric measurements were conducted with a precision scale, with a resolution of 0.0001 g. All weighings were conducted in triplicate and reported as a mean with standard deviation values.

X-ray microtomography (µCT)
The 3D images of dry bamboo samples were obtained using a Zeiss Xradia 510 Versa microtomograph to visualize and quantify the 3D structure. The experimental parameters for all samples can be found in Table 1.
The images were processed using the CTan software (version 1.18.8.0, Bruker, 2018) and Avizo© software (version 9.0.1, FEI, 2018). Image pre-processing using a non-local means filter was used to improve the signal-to-noise ratio. Smaller regions of interest were delimited around the sample's image to remove the background and reduce computational requirements. Watershed segmentation (Ther-moFisher Scientific 2018) as used to automatically separate the structures of interest (sclerenchyma, parenchyma, and vascular channels). Segmented images with each structure, i.e., parenchyma and vascular channels were then used for quantification.
The entire volume was measured, and the mean diameter was calculated for the whole 2D sections of the samples.
. 1 H time-domain nuclear magnetic resonance (TD-NMR) All measurements were performed on a MARAN Ultra (Oxford Instruments®) low-resolution NMR spectrometer operating at 0.54 T (23 MHz, 1 H Larmor frequency) equipped with an 18 mm diameter probe at 30 ± 1 °C. The 90° (p90) RF pulse was calibrated to 7.5 μs, and the 180° (p180) pulse was calibrated to 15 μs. Coil dead time and filter stabilization time were set to 6.0 and 3.0 μs, respectively.
We employed the CPMG (Carr-Purcell-Meiboom-Gill) pulse sequence with 2048 points, with a time between echoes (2τ) of 600 μs. For each sample, 64 scans were accumulated, with a recycle delay of 5 s. Total acquisition time was around 7 min. We carefully chose this interval between echoes to exclude the bulk of the relaxation signal from the lignocellulosic matrix's protons.
For measurements of the fluids' molecular dynamics in non-dynamic conditions, the samples saturated with each fluid were inserted in the TD-NMR tubes, and CPMG curves were acquired and stored for later processing. For the fluid desorption experiment, a given sample was filled with H 2 O with the aid of the vacuum chamber, then superficially dried with a paper towel, and then left inside an open TD-NMR tube at 30 °C. The equipment was then set to acquire and store a CPMG curve immediately after the sample insertion and then at 2 h intervals for five days (120 h).

TD-NMR data processing
Data processing was conducted with the aid of commercial software. For the Inverse Laplace Transform (ILT), we employed WinDXP® Toolbox 1.8.1.0 (Resonance Instruments®), with a smoothing parameter set at 2 × 10 -4 and tolerance at 10 -7 . A maximum of 50 iterations was allowed. None of the 2048 points obtained from the CPMG relaxation curve was excluded during the ILT processing, and all distribution curves were constructed with 256 points between 10 μs and 10 s. All curve fitting procedures were made in Origin 8.5 (OriginLab®). The Solid-Echo pulse sequence (SEPS) was used to evaluate the influence of the presence of the fluid inside the polymer matrix of the inner walls of the bamboo sample. The decay signal was acquired with 3072 points spaced by 1 µs, with a receiver gain of 40 dB and a recycle delay of 2 s. The obtained FID signal was adjusted according to the following function: where, A(t) is the total amplitude of the FID signal, F r is attributed to highly rigid molecular fractions and is written as a stretched exponential function to better fit the Gaussian contribution arising from residual dipolar coupling, F sr is associated with molecular fractions with intermediate molecular mobility and F m is the term related to highly mobile molecular fractions. T 2* refers to a transverse relaxation time, a constant related to the loss of phase coherence in the total magnetization of a given sample, provoked by inhomogeneities in the magnetic field and distinct physico-chemical environments in the sample. Finally, k is an offset parameter. Physical interpretations of the fractions and their respective relaxation times are given in the discussion section.
From the decomposition of the signal into the respective components, the average residual dipole coupling <M 2 > was extracted from the adjustment of the component Fr with the function below: Furthermore, the theoretical background for attributions of the behavior of liquid molecules in confined spaces such as pores, voids, and channels is given below. The ILT can be applied to the NMR decay signal to detect and discriminate the relaxation behavior of water protons in different families of pore sizes. The transverse relaxation rate of a liquid within porous materials is governed by the rates of free solvent relaxation (T 2F ), surface relaxation (T 2S ), (1) and diffusion resulting from molecular diffusion (T 2D ) (Zhang et al. 2016): 1/T 2D is the relaxation time caused by the diffusion of solvent molecules in the presence of internal field gradients being equal to (τγG) 2 D/3. The time between echoes is τ, γ is the gyromagnetic ratio, G is the fieldstrength gradient and D is the molecular diffusion coefficient. The internal gradient is dependent on the strength of the applied magnetic field. For PSD analysis where the time between echoes is short enough (2 ms or smaller) in low-field spectrometers, this term of the equation also becomes greatly reduced, allowing the simplification of the equation. Thus, considering that T 2F >> T 2S and a system within a fast diffusion limit, the equation can be approximated to (Al-Mahrooqi et al. 2006;Capitani et al. 2012;Meyer et al. 2018;Stagno and Capuani 2022): where S is the pore surface area, V is the pore volume, and ρ is the relaxivity constant of the confined fluid. For a given ρ, as S V is proportional to 1 R for a sphere, one sees that the transverse relaxation rate is directly proportional to the pore size (Zhang et al. 2016).
The surface contour maps used to report the behavior of the transversal relaxation domain curves of the water during its desorption at room temperature in the bamboo samples were generated by the Origin software. From the amplitudes of the 256 points generated by the ILT on the transverse relaxation data of the samples varying as a function of the desorption time, the software builds a matrix and assigns colors to the different values found in said matrix. Thus, it is possible to follow the evolution of the behavior expressed by the data over time and investigate dynamic processes in a visual and relatively simple way.
The values displayed in Table 4 were obtained using the value of the areas of the transverse relaxation curves generated by CPMG. For that, only the curves at the beginning of the desorption process were used (Time 0), corresponding to the maximum saturation value by fluids of those samples. Thus, for a given reservoir (free fluid, parenchyma, or vascular channels) the area under the corresponding curve was divided by the total area.

Results and discussion
The highly anisotropic 3D architecture of the bamboo vegetal tissues was investigated by µCT, as shown in Fig. 3. The variation and spatial distribution of the different biological tissues were analyzed through segmentation and post-processing µCT images. Sclerenchyma tissue with its lignin-rich fibers (green, Fig. 3b) presents a variation in morphology and density along the radial direction from the outer to the inner wall. Accordingly, the diameters of longitudinally-oriented vascular channels (blue, Fig. 3d) imbibed in the sclerenchyma tissue varied from 121 ± 27 to 192 ± 32 µm. Here it is helpful to note that we call the plant's phloem, xylem, and metaxylem channels collectively "vascular bundle" and "vascular channels." The parenchyma (yellow, Fig. 3c) is distributed throughout the volume of the sample. The volume occupied by parenchyma tissue cells and vascular channels in bamboo samples from each fraction is described in Table 2. Together, the vascular channels and the living cells of the parenchyma represent The cells of the parenchyma cortex and the fibrous bundles of the sclerenchyma, where the vascular channels are found, increase in size from the external to the internal fraction. This change directly affects the total fluid's mass that a given fraction of comparable dimensions can hold when fully saturated, as described in Table 3.
The discrepancy found in the masses of H 2 O and DMSO in the saturated samples can be attributed to the different wettability capacities of the molecular structure of bamboo tissues composed of cellulose, hemicellulose, and lignin. The filling fluids used in this study differ in molecular mass, polarity, and the number of hydrogen bonds they can establish, with H 2 O having a greater physicochemical affinity for the bamboo plant structure (Reichardt and Welton 2010).
The observation of cross sections obtained from µCT images (Fig. 4) of the three fractions indicates distinct volumes occupied by the plant's structures. For instance, in the internal fraction, one can clearly see the predominance of the parenchyma tissues. In fact, the abundance of the parenchyma tissue cavities in the internal fraction becomes evident when we observe the darkening of the yellow color of the generated CT image, suggestive of a less dense region. Although it is true that density plays a role in the different shades and colors on the images, we cannot ascribe the differences observed only to density gradients. Some of the variations in the color have been affected by the image processing software, even though we were trying to keep the same   color palette for all images. The cell wall of the parenchyma cortex is populated by pit holes, cellwall openings hundreds of nanometers in diameter, responsible for exchanging materials in the living plant. In the case of internal and middle fractions, the parenchyma cells present more surface area in contact with the vascular bundle channels offering extra possibilities of fluid exchange via pit holes (Pandoli et al. 2019). This can be seen in Fig. 4, where the channels (blue) are less embedded in the sclerenchyma tissues (green). In the external fraction, however, the architecture of the different tissues reduces the possibility of the facilitated fluid exchange found in the innermost fractions. In addition, the external fraction exhibits more massive fibrous bundles, unlike the innermost fractions of bamboo. TD-NMR is very useful in probing the molecular mobility of a liquid phase embedded in a solid matrix (Meyer et al. 2018;Stagno and Capuani 2022). Through interferences in the relaxation processes, such as those described earlier, we can extract very accurate data about the environment in which the molecules in a liquid are cofined. In Fig. 5a, b, we have the transverse relaxation time domain curves, obtained by the CPMG pulse sequence, of H 2 O and DMSO molecules in four different environments. The free water/DMSO curves correspond to the behavior expected from a fluid when it is free from spatial confinement. That is, most of its molecules can undergo free rotation, and have higher mobility and free diffusion processes, implying bigger T 2 values (in the order of seconds). A fluid confined in a solid matrix suffers the effects of interaction between its molecules and the surrounding walls. In this scenario, surface relaxation, as described in Eq. 4, plays a relevant part and lowers the molecules mobility when compared to the bulk solvent (Zhang et al. 2016). Therefore, the smaller the pore size, the lower the molecular mobility of the solvent. (Al-Mahrooqi et al. 2006;Capitani et al. 2012). These effects interfere with the magnetic relaxation of the fluid molecules, widening the bases of the domain curves and significantly reducing the relaxation time values, as seen in Fig. 5.
The area under a relaxation domain curve can be related to the number of protons that relax in that time interval. We observed that the water molecules confined in bamboo cavities have at least two large relaxation domains: one at shorter times (tens of milliseconds) and the other at longer times (around 100 ms). Relaxation time values are inversely proportional to the confinement experienced by the fluid molecules (Zhang et al. 2016). Therefore, we can attribute the shortest relaxation times to the water molecules occupying the parenchyma cavities (left side) and the longest times to the water molecules occupying the vascular channels (right side). The same rationale can be applied to DMSO molecules. There is a clear difference in the populations of fluid molecules in each of the bamboo fractions. In the internal fraction, the largest population is found inside the vascular channels; in the middle fraction, the population relative to the molecules of the fluid in the parenchyma starts to increase and, finally, in the external fraction, the amount of fluid contained in the parenchyma is slightly larger than that inside the vascular channels. Therefore, the importance of choosing which region of the plant becomes evident when technological applications of the material are envisioned. The results show that each fraction has a distinctive way of distributing fluid throughout its volume (Table 4).
From the distributions of transverse relaxation times, two relaxation domains are observed, both for the water-saturated bamboo and the DMSOsaturated bamboo. The domain with the lowest T 2 value (approximately 2 × 10 -3 s) is related to the solvent present in more confined spaces of the bamboo matrix, such as the parenchyma and sclerenchyma and the internal channels that connect the different regions of the bamboo. In addition, there is a fraction of the solvent on the surface of the inner walls of the vascular channels, where the effect of surface relaxation is more pronounced, which promotes the reduction of the T 2 value compared to the free solvent (bulk) molecules.
The domain with the highest T 2 value (approximately 0.1 s) is related to the solvent present in the vascular channels, more specifically in their central region, away from the internal walls. In this region's cross-section, the solvent molecules can be found in an unimpeded environment, translating to a relaxation regimen closer to their bulk state.
Similar trends are observed when we compare the internal, intermediate, and external bamboo fractions distribution for the water-and DMSO-logged samples. The domain with the lowest T 2 value shifts to the shorter relaxation times for the outer fraction of the bamboo. However, the domain with the highest T 2 value shifts to the right, increasing its value for the outermost fraction of the bamboo. The shortening of the relaxation times corroborates the data obtained from microtomography. As noted, parenchyma and sclerenchyma have the smallest pore sizes, which become even smaller as we advance from the inside to the outside of the bamboo. This higher confinement of the fluid translates into a reduction in the value of T 2 .
The increase in the T 2 value for the domain around 0.1 s can be explained by the more significant number of bundles of vascular channels present in the outermost fraction, as highlighted in Fig. 3. Although the diameters of the bundles are smaller in the outermost region of the bamboo, the channels have three-dimensional cylindrical structures, and each cylinder has a central region in which the fluids' molecules can relax their magnetization with less encumbrance; that is, they do not experience the effect of surface relaxation of the channel walls. Since the outermost fraction has more vascular channels, the amount of water in the bulk state is larger than in the innermost fraction of the bamboo. This, in turn, translates into an increase in the T 2 value for the second domain shown in Fig. 5.
The difference in distribution profiles between water and DMSO should be highlighted. While the first solvent has narrower and more well-defined domains, the transverse relaxation distribution for DMSO shows a broader profile. From the presented data, we can infer that, while the bamboo tissues' water-saturation capabilities are extensive, as expected, the DMSO molecules can permeate with ease through the plant matrix. This information may be relevant when choosing the best solvent in a microfluidics process (Machado et al. 2014).
The Solid-echo pulse sequence (SEPS) was employed to investigate the effects of the saturation of water or DMSO (Fig. 6) molecules on the structure of the bamboo matrix. Solvents of different polarities act on the lignocellulosic chains. These observations seek the best performance of bamboo in microfluidics operations. The FID measurements obtained by SEPS reveal three well-defined regions in the water-saturated (Fig. 6, on the left side) or DMSO-saturated bamboo structure (Fig. 6 on the right side). The rigid fraction (F r ) is attributed to the sclerenchyma (crystalline cellulose), the semi-rigid fraction (F sr ) is associated with the parenchyma, and the high mobility fraction (F m ) refers to the inner walls of the vascular bundles. These vessels contain the solvent molecules relaxing closest to the bulk state, with high T 2* values.
Relaxation signals for dry bamboo were obtained to evaluate the effect of the presence of water or DMSO on the three (F r , F sr, and F m ) relaxation regions of the bamboo. In the absence of solvent molecules, the bamboo polymeric matrix displays rigidity, as observed in the fast decay regimen of the NMR signal (Fig. 6, middle column). The signal was fitted by a Gaussian function with a transverse relaxation time of about 10 µs and a residual dipole coupling (M 2 ) close to 9000 rad 2 ms −2 . As expected, the presence of liquid molecules with chemical affinity with the natural polymer matrix increases the molecular mobility of the lignocellulosic structures constituting the samples, leading to the observed increase in the T 2* values. Similar behavior was observed in the NMR relaxation of wood samples saturated with different fluids (Zhang et al. 2016).
For the samples saturated with water, a rising trend in the F r contribution to the total magnetization in each sample is observed when we move from inner to outer fractions of the sample (35-55% of F r ). As seen in the µCT results (Figs. 3, 4), the external fraction contains larger quantities of sclerenchyma, which displays very little permeability. The percentage of F sr Fig. 6 Obtained SEPS NMR signal intensities fitted to Eq. 1 (signal components) and Eq. 2 (M 2 values) for the samples saturated with distilled H 2 O (left) and DMSO (right). The middle column refers to the dry bamboo samples remains higher than F m since the parenchyma volume is larger than the vascular bundles in the samples, as seen in Table 2.
Comparing the values of M 2 for the dry (M 2average = 9241 rad 2 ms −2 ) and water-saturated bamboo samples (M 2average = 8159 rad 2 ms −2 ) for each of the sampling fractions, the internal fraction shows a reduction of 16% in this parameter, the highest variation in all three sampling fractions. This reduction reflects a small volume of sclerenchyma tissues (Table 4), implying a higher water-soaking ability. Furthermore, the analysis of the M 2 values informs that solvation of the sclerenchyma with water and DMSO molecules has happened. We tentatively attribute this to the micro and nanometric pits interconnecting vascular bundles and sclerenchyma tissues, allowing for solvent flow between these regions.
The bamboo samples saturated with DMSO also display three well-defined regions, as described above. The F r values (from 47 to 65%) are higher when compared to the results obtained with distilled water (from 35 to 55%). This feature, associated with the higher values of M 2 , suggests that DMSO has lower permeability in the sclerenchyma than water. We attribute the rise in the F r percentage from 47 (internal) to 65% (external) to the higher volume of sclerenchyma tissue in the external fraction.
Another important aspect is the T 2* values for the region of high mobility (F m ) of the bamboo samples saturated with DMSO (T 2*average = 1563 ms), which are higher than those obtained with the water (T 2*average = 1089 ms). This suggests that DMSO molecules tend to occupy, preferably, the vascular bundle. Therefore, it is observed that ultrapure water, as expected, has a greater capacity to diffuse through the bamboo's internal structures, mainly due to polarity and molecular size and the ability to perform hydrogen bonds. Both solvents remain primarily in the parenchyma and vascular channels, but water also has a small capacity to permeate the sclerenchyma.
The evolution of the magnetic relaxation domain curves in TD-NMR as a function of time can be represented through contour surface maps (Rodrigues et al. 2018). The intensities and relaxation time values of the domains at each moment represent the dynamics of the water evaporating inside the bamboo cavities as a function of time. The water evaporation from bamboo samples was monitored in real-time by TD-NMR to describe the behavior and distribution of  Figure 7 shows the results obtained.
Initially, 120 h of process monitoring was recorded, but from 72 h onwards, the remaining amount of water becomes very small, and the signalto-noise ratio greatly decreases. This compromises the quality of the Inverse Laplace Transform that generates the domain curves, which significantly reduces the reliability of the results.
A qualitative analysis of each contour surface map shows that the behavior of confined water during evaporation is different for each bamboo fraction. There are two large populations of water molecules represented by the two trails on the map: one at longer relaxation times relative to water in vascular channels and one at shorter relaxation times relative to the environment of the parenchyma cortex. The coalescence of the confined water molecules' transverse relaxation domains progressively increases from the outer to the inner fraction, again indicating the existence of varying degrees of connection between the two cavities in the tissue. We propose that this behavior becomes apparent in TD-NMR experiments from a certain threshold of water saturation, such as those used in this work.
Previous experiments with wood samples containing varying amounts of water did not report the present findings, as they operated inside a framework of water vapor saturation or natural humidity present in their samples (Mikac et al. 2021;Cocusse et al. 2022;Li and Ma 2022). The extent of communication between the parenchyma cavities and the vascular channels became evident during the desorption experiment. The population of water molecules in vascular channels disappears after 48 h in the three fractions, albeit at different evaporation rates for each fraction. The water present in the smaller cavities of the parenchyma does not disappear completely, although the amplitude of these domains decreased considerably compared to the beginning of the evaporation process. It has been demonstrated that a residual amount of water is retained in lignocellulosic materials due to hydrogen interactions with the cellulose chains in plant tissue (Zhang et al. 2016;Li et al. 2017c).
Both the coalescence of the relaxation domains of water molecules and the time until the amplitude of these domains is considerably reduced by evaporation substantiate the results of µCT. Previous studies have shown that the interconnected pore structure constitutes a highly specialized fluid transportation system (Liu et al. 2021). The efficiency of this system is evident when we examine the behavior of the confined water volume in the samples of the internal and middle fractions. In highly confined environments, the trapped water molecules have their evaporation rate delayed at temperatures close to room temperature, like those used in our experiments, which may explain the remaining amplitude of the observed relaxation domains even though it has long test times. We demonstrated that the efficiency of fluid exchange between the two types of microcavities (vascular bundles and parenchyma cells) becomes a relevant factor affecting the fluid's evaporation rate.

Conclusions
µCT and TD-NMR showed great potential for analysis of the internal structure of bamboo, as well as the distribution of solvents of different polarities and their effects on the molecular mobility of the fully solvent-saturated bamboo structure. The µCT results supported the reasoning used in the interpretation of the TD-NMR data. Most of the interpretations of the behavior of the molecular dynamics of the fluids confined in the bamboo structure were only possible with the aid of this technique's qualitative and quantitative information.
The CPMG sequence, through the values of T 2 and the area under the distribution of domains, can map the distribution of solvents in the parenchyma and vascular channels, being sensitive to the morphology changes as a function of the bamboo cut region (internal, intermediate, or external), corroborating the data obtained by microtomography. Ultrapure water and DMSO showed similar distribution profiles. Monitoring water desorption indicates that it takes 40-48 h for almost all of the solvent to leave the bamboo matrix, leaving only a small fraction confined in the parenchyma.
The Solid Echo sequence, through the T 2* values and the percentages of rigid, semi-rigid, and mobile fractions, was sensitive to the mobility variations of the lignocellulosic chain segments in the presence of ultrapure water or DMSO. The results indicate that water has a greater permeation capacity and affects the rigid fraction of bamboo more, while DMSO settles more easily in vascular channels.
We believe that the data found in this work may be fundamental in choosing the most suitable solvent for a particular type of bamboo during a microfluidic process, as well as its effects on the structure of the segments of the bamboo chain.
In general, the results indicate dependence between the physicochemical characteristics of the fluids and lignocellulose-based microfluidic reactors with hierarchical architecture, such as bamboo. Studies with other fluids and ionic liquids are being conducted by our group to understand better the molecular behavior in confined environments made up of natural polymers. Data availability All data generated or analyzed during this study are included in this published article.

Conflict of interest
There are no conflicts of interest/competing interests.

Ethical approval
No results of studies involving humans or animals are reported.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.