---
title: 'FACT: Fibril Analysis for Cellulose Technology'
url: https://www.emergentmind.com/topics/fibril-analysis-for-cellulose-technology-fact
type: topic
---

# FACT: Fibril Analysis for Cellulose Technology

Fibril Analysis for Cellulose Technology (FACT) denotes a set of integrated analytical workflows for cellulose-based materials in which modality-specific measurements are converted into quantitative descriptors of fibril morphology, orientation, porosity, mechanics, interfacial state, or process response. In the reported literature, FACT is used for high-resolution synchrotron phase-contrast microtomography and 3D pore morphometrics in flax fibres, for machine-learning segmentation and morphological thinning in negative contrast SEM images of cellulose nanofibers, for AFM–CLSM–SEM mapping of humidity-dependent local mechanics in single cellulose fibres, and for multiscale orientation analysis in neutron tomography and flow-stop birefringence experiments [2303.18127, 2509.06618, 2012.10207, 2407.06728, 1801.07558]. The term therefore refers not to a single instrument or a single algorithm, but to a methodological framework in which fibril-resolved measurements are standardized, spatially registered, and interpreted against processing history and macroscopic performance.

## 1. Conceptual scope and analytical logic

FACT is unified by a recurring sequence: acquisition of structurally informative data, segmentation or reconstruction of cellulose-relevant features, extraction of physically interpretable metrics, and correlation of those metrics with mechanics, transport, or processing behavior. In the flax-fibre study, the framework is explicitly described as leveraging high-resolution synchrotron phase-contrast microtomography, advanced 3D segmentation, and quantitative morphometrics to link cell-wall ultrastructure, porosity, microfibrillar angle perturbation, and mechanical performance [2303.18127]. In the CNF image-analysis study, FACT is defined through binary segmentation, topology-preserving thinning, and whole-network width measurement from NegC-SEM images [2509.06618]. Auernhammer et al. formulate a FACT workflow in which maps of $E(x)$, $F_a(x)$, $E_{diss}(x)$, $S(x)$, and $\theta(x)$ are correlated to segment a single fibre into regions of distinct behavior under controlled relative humidity [2012.10207].

This common logic makes FACT inherently multiscale. The reported workflows encompass sub-elementary features at $\approx 0.44$ nm and $\approx 0.88$ nm in TEMPO-oxidized nanofibers, pore layers with median radial width $0.957\,\mu$m in flax-fibre kink-bands, and centimetre-scale foam specimens imaged by neutron tomography [2010.15686, 2303.18127, 2407.06728]. Taken together, these reports suggest that FACT is best understood as a quantitative bridge between cellulose ultrastructure and engineering observables, rather than as a narrowly defined characterization protocol.

## 2. Structural imaging, reconstruction, and segmentation

A central branch of FACT is high-resolution imaging of fibrillar architecture and its defects. In flax fibres, measurements were performed on individual fibres at the ANATOMIX beamline using a polychromatic X-ray beam centered at $\sim 12$ keV, with a $2048\times 2048$-pixel detector, effective pixel size $0.325\,\mu$m, voxel size $0.325\,\mu\text{m}^3$, a $0.65\times 0.65$ mm$^2$ field of view, and total scan time per $2$ mm fibre segment of $\approx 4$ min. The raw projections were reconstructed with PyHST2 using filtered back-projection and a Paganin phase retrieval with $\delta/\beta$ kernel length $=8\,\mu$m, followed by a non-local means filter, global thresholding of voids, manual separation of the lumen from defect-induced pores in Avizo 2021.1, and “Separate Objects” analysis for individual-pore morphometrics [2303.18127]. This pipeline isolates ultrastructural void organization inside kink-bands rather than treating the fibre as a homogeneous solid.

In image-based CNF metrology, FACT requires a binary image separating foreground from background and implements two pixel-classification approaches: Weka for small datasets and rapid per-image training, and a modified U-Net CNN for larger, heterogeneous datasets. The U-Net reduces parameters from $31$ M to $7.7$ M by halving feature-map counts in each layer; the input workflow mirror-pads each $2048\times 2048$ image to $2736\times 2736$, extracts overlapping $700\times 700$ patches with $342$ px overlap, applies rotations by $0^\circ$, $90^\circ$, $180^\circ$, $270^\circ$ and mirror at $0^\circ$, and splits data into $62\%$ train, $33\%$ validation, and $5\%$ test, totaling $\sim 2\,560$ patches from $8$ images. Inference takes $\sim 5$ s per $2048\times 2048$ image on NVIDIA RTX A5000, and the full pipeline with a pre-trained U-Net is reported as less than $5$ min per image [2509.06618]. The key structural operation after segmentation is iterative thinning,
$$
S_{n+1}=T(S_n),
$$
continued until convergence to a one-pixel-wide skeleton.

A further extension of FACT is multiscale 3D orientation imaging by multi-directional dark-field neutron tomography. Here a single-absorption grating with period $p=350\,\mu$m creates a 2D spatial intensity modulation, the sample is placed $\sim 25$ cm upstream of the detector, and the real-space correlation length is set by
$$
\xi=\lambda L_s/p,
$$
with $\xi\approx 200$ nm under the stated geometry. Three in-plane scattering-sensitivity directions, $0^\circ$, $60^\circ$, and $120^\circ$, are acquired over $72$ projections across $360^\circ$, each with $15$ min exposure, and reconstructed with the SIRT algorithm in the ASTRA toolbox [2407.06728]. Unlike X-ray or electron workflows, this route is explicitly positioned as non-destructive and suitable for fragile hierarchical biomaterials.

## 3. Quantitative descriptors: porosity, width, and orientation

FACT is defined as much by its metrics as by its instrumentation. In the flax-fibre implementation, local porosity in a sub-volume is
$$
\phi=\frac{V_{pores}}{V_{total}},
$$
with $V_{pores}=\sum_i V_p(i)$, and the local orientation proxy for cellulose microfibrils is expressed through the inclination of pore-layer long axes relative to the lumen axis, with first two statistical moments
$$
\mu_\theta=\frac{1}{N}\sum_{j=1}^N \theta_j,\qquad
\sigma_\theta=\sqrt{\frac{1}{N}\sum_{j=1}^N (\theta_j-\mu_\theta)^2 }.
$$
Using this framework, $12$ kink-bands were studied over a cumulative fibre length of $2.66$ mm; pore-layer thickness had median $0.957\,\mu$m with range $0.447$–$2.235\,\mu$m; inter-layer spacing had median $0.412\,\mu$m with range $0.325$–$0.612\,\mu$m; and $517$ individual pores were extracted, with volumes from $0.04$ to $289\,\mu\text{m}^3$ and median $V_p\approx 0.34\,\mu\text{m}^3$. Local $\phi$ rises from $<1\%$ in intact regions to peaks approaching $6\%$ in kink-band zones, while pore inclination angles lie between $22^\circ$ and $51^\circ$, with $\mu_\theta\approx 35^\circ$ and $\sigma_\theta\approx 8^\circ$ over $N=13$ measurements [2303.18127]. The linear increase of pore-layer radius versus distance from the lumen with slope $\approx 1$ is taken to demonstrate a concentric “onion-peel” arrangement at successive interfaces of the S2(G) growth layers.

In whole-network CNF morphology analysis, the defining FACT observable is width derived from the Euclidean distance transform. If $K$ is the final skeleton and $B$ the unfiltered binary foreground, then
$$
D(u,v)=\min_{(x,y)\in \partial B}\sqrt{(u-x)^2+(v-y)^2},\qquad
w(u,v)=2D(u,v).
$$
The width set $\{w_i\}_{i=1}^N$ is summarized by
$$
\bar w=\frac{1}{N}\sum_{i=1}^N w_i,\qquad
\sigma_w=\sqrt{\frac{1}{N-1}\sum_{i=1}^N (w_i-\bar w)^2 }.
$$
Validation on idealized rectangular branches with true widths $30$, $62$, $120$, and $240$ px yielded FACT peak means $[30.2,62.3,119.8,240.5]\pm 1$ px, corresponding to error $<1\%$. On low-branching CNFs, five images at $5.43$ nm/px gave FACT $\bar w=22.1\pm 8.9$ nm versus manual $\bar w=22.7\pm 8.2$ nm, and one image at $1.79$ nm/px gave FACT $16.9\pm 5.2$ nm versus manual $17.1\pm 2.5$ nm. For high-branching CNFs, FACT produced $\bar w=510\pm 532$ nm versus manual $\bar w=320\pm 310$ nm, with the reported skew toward larger widths attributed to the fact that longer, thicker fibrils contribute more skeleton pixels [2509.06618].

Orientation analysis is also formalized in tensorial or order-parameter terms. In neutron tomography, the second-moment orientation tensor is
$$
\mathbf{T}=\langle \mathbf{u}\otimes \mathbf{u}\rangle,
$$
and Herman’s orientation factor in a reference direction $\mathbf{e}$ is
$$
S=\left\langle P_2(\cos\gamma)\right\rangle.
$$
The dark-field implementation does not reconstruct a full continuous ODF, but instead infers local anisotropy through an eccentricity $e$ computed from orthogonal dark-field channels. Cross-validation by SAXS gave Herman factors at $d\approx 150$ nm of $S\approx -0.395$ for the CNC shell, $S\approx -0.162$ for the CNF shell, and $S\approx +0.864$ for the CNF-uit core [2407.06728]. A plausible implication is that FACT descriptors can remain comparable across modalities only when the underlying orientation proxy is stated explicitly.

## 4. Mechanical and hygro-mechanical mapping

A second major branch of FACT concerns local mechanics under controlled environmental state. Auernhammer et al. describe a single cotton-linter fibril suspended over a $1$ mm trench and clamped at both ends inside a humidity-controlled AFM chamber at RH $=2$, $40$, $75$, and $90\%$, with equilibration for $\approx 45$ min after each RH step. A spherical SiO$_2$ colloidal probe of diameter $D=50\,\mu$m is attached to a V-shaped cantilever with $k\approx 2$–$5$ N/m, and static force–distance curves are acquired at points spaced by $\Delta x=5\,\mu$m. Local stress and strain are computed as
$$
\sigma(x)=F(x)/A(x),\qquad \epsilon(x)=\Delta L(x)/L,
$$
and the Young’s modulus map is extracted from the linear regime of $\sigma$–$\epsilon$ curves. Combined with CLSM-based swelling, $S(x,\text{RH})=[d(x,\text{RH})-d_0(x)]/d_0(x)$, and SEM-derived fibril orientation $\theta(x)$, the workflow identifies “wet spots” characterized by local maxima in $S(x,40\%)$ and minima in $E(x,40\%)$ and $D(x)$, accompanied by peaks in $F_a(x)$ and $E_{diss}(x)$ [2012.10207]. The reported quantitative ranges are sharply humidity dependent: $E(x,\text{RH})$ changes from $25$–$95$ GPa at RH $=2\%$ to $4$–$18$ GPa at RH $=90\%$; adhesion rises from $\sim 50$ nN at $2\%$ RH to $\sim 700$–$800$ nN at $90\%$ RH; dissipated energy increases from $1\times 10^{-17}$ J to $\approx 1\times 10^{-16}$ J; and swelling at $40\%$ RH reaches up to $+70\%$ in “loose” ROIs and $+30\%$ in “tight” ROIs.

Humidity sensitivity is corroborated in nanofibrillated cellulose films. Simão et al. report dry, local water-free Brillouin Light Scattering values $\nu_{LA}=3427\pm 22$ m s$^{-1}$ and $\nu_{TA}=1783\pm 19$ m s$^{-1}$, corresponding to $C_{11}=17.61\pm 0.23$ GPa, $C_{12}=8.07\pm 0.29$ GPa, bulk modulus $B=11.0\pm 0.5$ GPa, Poisson ratio $\mu=0.31\pm 0.01$, and $E_{BLS}=12.5\pm 0.4$ GPa. Under QNM-AFM, the same material shows $E_{QNM}=16.4$ GPa at $0\%$ RH and $6$ GPa at $100\%$ RH, while adhesion changes from $8$ nN to $\approx 0$ nN, with repeatability of $E$ and $F_{ad}$ better than $5\%$ over three RH cycles [1504.00230]. The stated modulus change, $\Delta E/E_0\approx 65\%$, emphasizes that RH is not a secondary test variable but a primary state variable in FACT mechanics.

At the cell-wall scale, Amando de Barros et al. combine micropillar compression and DIC on Norway spruce tracheids to resolve MFA-dependent stiffness and failure. Pillars with height $l_0\approx 5\,\mu$m and diameter $2$–$3\,\mu$m are tested at MFA $\approx 0^\circ$, $20^\circ$, $70^\circ$, and $90^\circ$ under displacement-controlled compression to $\delta_{meas}=1\,\mu$m at nominal strain rate $10^{-3}$ s$^{-1}$. DIC-based measurements at $2$ kV imaging yield $E=42\pm 3$ GPa and $\sigma_y=0.18\pm 0.04$ GPa for MFA $=0^\circ$, $E=26\pm 5$ GPa and $\sigma_y=0.14\pm 0.03$ GPa at $20^\circ$, $E=11\pm 1$ GPa and $\sigma_y=0.10\pm 0.02$ GPa at $70^\circ$, and $E=7\pm 1$ GPa and $\sigma_y=0.09\pm 0.02$ GPa at $90^\circ$. Low-MFA pillars show fibril-aligned kink bands, whereas high-MFA pillars show shear-related catastrophic failure. Continuous $5$ kV exposure causes extensive surface shrinkage and more than $50\%$ drop in $E$ and $\sigma_y$ versus “no-beam” measurements [2506.11177]. This result directly connects imaging protocol to apparent cell-wall mechanics.

## 5. Chemical disassembly and interfacial interactions

FACT also includes chemically driven fibril-state control. In TEMPO-mediated oxidation of sugarcane bagasse cellulose pulp, the catalyst system consists of $0.1$ mmol TEMPO and $1.0$ mmol NaBr per gram of cellulose, with three NaClO loadings: SC-5 at $5$ mmol NaClO/g cellulose, SC-25 at $25$ mmol NaClO/g, and SC-50 at $50$ mmol NaClO/g. The reaction is maintained at $\text{pH}=10.0\pm 0.2$, temperature $\sim 25^\circ$C, and total reaction time $\sim 2$ h, then quenched with ethanol and washed until conductivity plateau. The degree of oxidation rises from $0.40$ mmol COO$^-$/g for SC-5 to $1.10$ mmol/g for SC-25 and $1.40$ mmol/g for SC-50, with zeta potentials $\zeta\approx -15$ mV and $\zeta\approx -65$ mV for the more oxidized states. AFM analysis of $\sim 300$ individual nanofibrils shows aggregated bundles $>10$ nm in SC-5, but individualized nanofibers of average length $243$–$370$ nm with a major width peak at $\approx 4.0$ nm in SC-25 and SC-50, plus sub-elementary populations at $\approx 0.44$ nm and $\approx 0.88$ nm assigned to single-chain and double-chain nanofibers [2010.15686]. First-principles calculations quantify the accompanying energetic changes: in pristine CNFs, $E^b_{IS}=-1.351$ eV/unit-chain and $E^b_{IC}=-0.676$ eV/unit-chain, whereas at $50\%$ carboxylation and $q=0.5e$, $E^b_{IS}$ weakens from $-0.60$ to $-0.51$ eV and $E^b_{IC}$ from $-0.48$ to $-0.23$ eV. The stronger relative weakening of interchain O–H$\cdots$O hydrogen-bond interactions is presented as the mechanistic basis for liberation of single and double chains.

Interfacial adsorption assays based on cellulose nanocrystals extend FACT to colloidal formulation science. Cotton linters hydrolyzed with $65\%$ (w/v) H$_2$SO$_4$ at $63^\circ$C for $30$ min yield lath-shaped particles with length $L=180\pm 30$ nm, width $w=17\pm 4$ nm, thickness $h\simeq 7$ nm, and $\zeta_{CNC}=-38$ mV at $0.1$ wt\% and pH $4.5$. The cationic surfactant TEQ has $\zeta_{TEQ}=+65$ mV, self-assembles into unilamellar vesicles of $100$–$500$ nm diameter at $c\leq 0.2$ wt\%, multivesicular structures of $300$ nm–$1\,\mu$m at $c\simeq 1$ wt\%, and large bilayer stacks above $c\geq 8$ wt\%. Continuous Variation plots of Rayleigh ratio $R(X)$ and hydrodynamic diameter $D_H(X)$ show maxima at $X_{max}\simeq 0.3$, whereas electrophoretic mobility crosses zero near $X\simeq 1$, indicating charge-driven neutralization. The binding constant $K_b$ extracted from a Langmuir-type isotherm is stated to fall typically in the $10^5$–$10^6$ M$^{-1}$ range for strong electrostatic adsorption [1703.10668]. Here FACT functions as a sensitive surrogate assay for deposition on cotton.

Ion-specific perturbation of cellulose provides a further chemical axis. Replica-exchange MD of cellotetraose in NaCl solution and NPT-LD simulation of a cellulose I$_\beta$ fibril show preferred Na$^+$ contacts at O$_2$, O$_3$, and O$_6$. For the solvated tetramer at $298$ K, the reported probabilities are $P(\text{Na--O}_2)=0.32$, $P(\text{Na--O}_3)=0.37$, and $P(\text{Na--O}_6)=0.30$, with O$_2$–O$_3$–Na$^+$ bridging at $P\approx 0.22$. For the fibril surface at $298$ K, $P(\text{Na--O}_2)=0.24$, $P(\text{Na--O}_3)=0.33$, and $P(\text{Na--O}_6)=0.41$, while core hydroxymethyl populations shift from tg/gt/gg $=92.8/4.6/2.6\%$ in pure water to $59.9/25.2/14.9\%$ with NaCl [1306.1818]. The paper interprets this as disruption of native intrachain hydrogen bonds and promotion of alternative intersheet interactions.

## 6. Process analytics, limitations, and technological significance

FACT is not limited to static imaging or ex situ assays; it also includes process-relevant orientation dynamics. In the flow-stop POM method for semi-dilute CNF dispersions, a thin channel is placed between crossed polarizers at $+45^\circ$ and $-45^\circ$ relative to the flow axis, illuminated by a $660$ nm laser over a $\sim 10$ mm spot, and imaged at $100$ fps while fast three-way valves arrest the flow quasi-instantly. The birefringence-derived signal obeys
$$
I_{POM}(z,t)\approx I_0\Bigl(\frac{2\pi d}{\lambda}\Bigr)^2(\Delta n(z,t))^2,
$$
and, under a purely Brownian dilute-rod model,
$$
I_{POM}(z,t)=I_{POM,0}(z)\exp(-12D_r t).
$$
Piecewise fitting reveals two rotary-diffusion regimes, with $D_{r,\text{fast}}\approx 1$–$1.3\ \text{rad}^2\text{s}^{-1}$ over $0$–$0.1$ s and $D_{r,\text{slow}}\approx 0.05\ \text{rad}^2\text{s}^{-1}$ over $2$–$4$ s. The fast process depends on prior deformation history, and de-alignment is faster in shear-dominated flow than in pure extensional flow [1801.07558]. The same report explicitly proposes FACT quality criteria such as maximum $D_{r,\text{fast}}$ at design shear below $1.0\ \text{rad}^2\text{s}^{-1}$, $D_{r,\text{fast}}/D_{r,\text{slow}}<10$, and steady-state $S_{\phi,0}$ curves within $\pm 10\%$ of a validated reference dispersion.

The methodological breadth of FACT also creates recurring limitations. In ML-based NegC-SEM analysis, reliable performance requires high-contrast images and at least $5$ px across the narrowest fibril; segmentation errors in low-contrast or highly overlapped regions produce extraneous skeleton spurs, junction-point handling depends on SST and SSF tuning by trial-error, and width resolution is limited to $2\times$ the pixel size [2509.06618]. In cell-wall micropillar compression, uncontrolled electron-beam exposure degrades pillar integrity and can explain data scatter and mechanical underestimation in earlier studies, leading to the recommendation that continuous SEM imaging be limited to $\leq 2$ kV with beam currents $\leq 50$ pA and short acquisition windows [2506.11177]. In neutron dark-field tomography, only three scattering directions are probed, attenuation can bias dark-field unless corrected, and the spatial resolution is limited by detector pixel size, so complementary SAXS or SEM remains necessary below the $\sim 20$ nm regime [2407.06728]. In humidity-sensitive mechanics, moisture strongly modulates $E$ and adhesion, so RH must be standardized or at minimum reported alongside mechanical values [1504.00230].

These limitations do not diminish the integrative value of FACT; rather, they define the conditions under which its descriptors are interpretable. In flax-fibre defects, localized porosity spikes and MFA misalignment are linked to stress concentrators and weakened fibre-matrix adhesion in composites [2303.18127]. In single-fibre AFM–CLSM–SEM mapping, the output is a multi-parametric fibre “fingerprint” that feeds improved finite-element or network models [2012.10207]. In oxidation-controlled disassembly, measured $E^b_{IS}$ and $E^b_{IC}$ are presented as inputs for process-simulation modules and for targeting diameter distributions under given ionic strength and oxidation conditions [2010.15686]. Collectively, the reported literature suggests that FACT is becoming a general framework for converting fibril-resolved measurements into design rules for composites, membranes, paper networks, foams, and formulation-relevant cellulose interfaces.

Source: https://www.emergentmind.com/topics/fibril-analysis-for-cellulose-technology-fact