1. Introduction
According to recent study from 2023, chemical industries contributes approximately 5 % to global CO
2 emissions [
1], securing its position as the third-largest emitter due to reliance on fossil fuels [
2]. To address this issue, reinvention through innovations in raw materials and reaction engineering is needed to generate environmentally friendly products based on green chemistry principles.
One such approach is the holistic biorefineries model, where the entire lignocellulosic biomass (LCB) from waste is valorized into various high value bioproducts using interesting eco-friendly processes. Nowadays, almost 70 % of the LCB, consisting of cellulose and hemicelluloses polysaccharides [
3], is efficiently valued into chemicals [
4], biofuels [
5], pulp [
6], fibers and nanofibers [
7], while the remaining 20-30 % composed of lignin biopolymer [
8], mainly serves as an energy source [
9]. Indeed, due to its important binding with polysaccharides and complex chemical structure, lignin has often been seen as a critical barrier in biorefinery model and therefore categorized as a waste [
4].
Lignin, the predominant aromatic biopolymer in nature [
10], involves an amorphous and irregular shape mainly composed of three monomeric units, named p-hydroxyphenyl (H), guaiacyl (G), and syringyl (S) interconnected among others [
11]. Considering the annual global production of biomass waste estimated at 140 Gt [
12], lignin is an aromatic, renewable, biodegradable and abundant material [
13], making it a central resource in mitigating carbon emissions associated with traditional chemical processes.
Given the intricate complexity and heterogeneous nature of lignin structure, pretreatment technologies play a pivotal role in facilitating its applications, as it is the case with organosolv pulping method [
14]. In a pressurized heated and agitated reactor, organosolv pretreatment involves LCB breakdown and lignin solubilization in organic solvents such as EtOH providing environmental benefits by using non-toxic chemicals, minimizing emissions, reducing equipment corrosion and ensuring low energy consumption through the recovery of high-volatility alcohol via simple distillation [
15]. Unlike the kraft and sulfite processes of the conventional paper industry, which increase lignin dispersity [
10], organosolv technology enables the extraction of sulfur-free and pure lignin with a poor condensed structure that closely match to the native one, opening up possibilities for improved valuations [
16,
17]. Furthermore, organosolv also improves cellulose isolation with due to organic solvents positioning this technology as a promising asset for the biorefinery model [
18].
Positioned as the potential driving force behind the next industrial revolution, bio-based nanoscale technology may be an alternative approach to overcome lignin heterogeneity while offering a wide array of benefits tailored to meet evolving consumer needs across various sectors including energy, transportation, agriculture, food, materials, electronics, and medicine [
19]. Recently, manipulating lignin at the nanoscale level has been considered as a significant process for a sustainable biorefinery model [
20]. Interestingly, reduce lignin size enhances both morphological and chemical uniformity, yielding advancements in its structural integrity and compositional consistency [
21]. Due to their heightened surface area to volume ratio, LNPs exhibit distinctive adjustable and multifunctional properties, including improved antibacterial, anti-oxidant and UV protection properties such as greater homogeneity [
22,
23]. The growing interest around LNPs has encouraged significant advancements across fields including energy, environment, materials, medicine, cosmetics and feed [
24,
25], owing their advantages of biocompatibility, non-toxicity, medium-term biodegradability, environmental resistance, and improved properties [
26,
27]. The demonstrated low-cytotoxicity of LNPs has been crucial for generating industrial concern and ensuring the viability of lignin-based biorefinery models by opening up new opportunities in markets involving human contact [
28].
LNPs can be synthesized using diverse methodologies, including assisted or self-assembly formation, with supercritical CO2 treatment [
29], mechanical or ultrasonication homogenization [
30,
31,
32], aerosol processing [
33], electrospinning [
34], ice segregation [
35], cross-linking/polymerization [
36], solvent exchange and pH shifting [
37,
38,
39,
40]. The most common and promising technique for leverage green chemistry approaches to industrial applications is based on the self-assembly of LNPs through the antisolvent precipitation [
26]. This method also referred to nanoprecipitation and ouzo effect [
41], relies on lignin dissolution into water-miscible solvents, (e. g., ethanol, acetone, THF, DMSO among others) shifted with an excess amount of an antisolvent (e.g., water) [
42]. The LNPs formation is facilitated through orderly rearrangement and aggregation of lignin macroparticles (LMPs), harnessing lignin amphiphilic nature and both hydrophilic and hydrophobic interactions with solvents, driven respectively by aromatic structures and internal H-bonding [
43].
Despite the numerous advantages offered by LNPs, their production and valorization remain significant industrial challenges due to the complexity of controlling multiple parameters in a large-scale process [
17]. Indeed, in a previous study [
44], Girard et al. reinforced earlier findings [
16,
43,
45,
46] that highlighted the significance of the inherent structures of lignin oligomers, impacted by both biomass nature (hardwood, softwood, herbaceous material) and chemical extraction methodology, as a key factor affecting the final LNPs properties (size, morphology, stability) and therefore the associated industrial valuations. It was demonstrated that the lignin self-assembly process was partly influenced by lignin molecular weight and its amphiphilic nature, intricately related to the content of phenolic and aliphatic units, the building unit type (H,G,S) and non-covalent interactions such as π–π interactions [
47]. To this end, hardwoods materials were identified as optimal feedstocks for LNPs production, leading to enhanced nucleation characterized by formation of smaller and spherical particles compared to softwoods and herbaceous species [
44]. Similarly, within the context of a specific biomass specie, the organosolv extraction method produced improved lignin structure without sulfur [
44], leading to enhanced aggregation and greater LNPs properties compared to other extraction methods such as kraft process [
48].
Then, other studies also demonstrated that the self-assembly process is also partly driven by precipitation method itself, which depends on parameters such as antisolvent flow rate [
37,
39,
49,
50,
51], initial lignin concentration [
43,
52,
53], solvent temperature [
52,
54], solvent type [
40,
49,
52,
55], solvent ratio [
43,
49,
52,
56], stirring speed [
51,
52], and pH value [
37,
40,
51]. However, a significant issue arises from the fact that the majority of previous studies have analyzed factors in isolation over kraft lignin, which fails to adequately address process variability and serve as a basis for potential industrialization.
Hence, for the first time, in this work the efforts are built upon the previous findings to devise a controlled, environmentally friendly method using only water and ethanol for producing tailored and predictable LNPs from organosolv isolation. This approach may easily be transposed at industrial scale to faster the development of value-added applications. To achieve this objective, LMPs were first extracted from beech by-products using 10 liters semi-industrial organosolv reactor, generating sufficient amounts of material for parametric studies. Then, a parametric study over LNPs production from beech lignin using only water and ethanol has been performed, in order to define the relevant range of values in which each parameter can vary. It has been followed by a fractional factorial design plan, including parameters studied in the parametric study: (1) initial lignin concentration (1 – 50 g/L), (2) solvent flow rate (0.5 – 500 mL/min), (3) antisolvent composition (100/0 – 50/50 H2O/EtOH v/v), (4) antisolvent volume (1/2 – 1/20 solvent/antisolvent v/v), and (5) antisolvent stirring speed (150 – 1200 rpm). The extracted LMPs underwent detailed characterization utilizing nuclear magnetic resonance spectroscopy (NMR) and size exclusion chromatography (SEC). The size distribution of LNPs was analyzed using dynamic light scattering (DLS) and transmission electron microscopy (TEM). This approach enabled both to rank over the impacts of the LNP manufacturing parameters and process parameters, as well as to predict LNPs sizes, thanks to a model which rely on manufacturing parameters. These results can be used to design an environmentally friendly manufacturing process of lignin and LNPs at production scale.
4. Conclusions
This study has extensively explored the antisolvent precipitation method for producing LNPs, which is essential for the future of biorefineries and large-scale lignin valorization. For several years, research has focused on transitioning to the nanoscale to enhance the intrinsic properties of polymers. However, in the case of lignin, current applications are primarily hindered by the highly condensed and sulfur-rich chemical structure of kraft lignin, as well as by existing nanoscale reduction processes that are energy-consuming, lack control and repeatability, and rely on toxic solvents such as THF and DMSO.
Here, we generate a pure lignin from beech biomass residues using a semi-industrial organosolv process, which is chemically ideal for fabrication of nanoparticles as shown by Girard et al.[
44]. By studying 5 different parameters of the antisolvent precipitation method with ethanol and water for the first time, this work demonstrates that the antisolvent composition and solvent/antisolvent ratio are the most crucial factors for controlling the nucleation phenomenon according to the classical nucleation theory. Starting from an average size of 577 nm, the LNPs were reduced to 134 nm by shifting the antisolvent composition from 50/50 v/v water/ethanol to 100/0 v/v water/ethanol. Additionally, increasing the solvent/antisolvent ratio from 1/2 to 1/20 reduced the average LNPs size by 315 nm. The proposed optimization of these different parameters in
Section 3.3 led to an improvement of PDI, ζ-potential and particle morphology, resulting in tailored environmentally friendly suspensions with desired properties. Additionally, by conducting these experiments within an experimental design framework, the former hypothesis regarding the most important factors was demonstrated and an accurate predictive model for LNPs size properties was developed. With a linear correlation coefficient of R² = 0.982 and a model validity of 97.9 %, the anti-solvent precipitation model presents essential information to implement the process towards production scale.
This study describes a globally controlled top-down process using biomass residues to produce tailored lignin nanoparticles with enhanced properties and reduced environmental impact. The process allows for optimization and prediction of LNP characteristics for industrial applications paving the way for future developments.
Figure 1.
(A) Summary of the different experimental parameters used for LNPs synthesis. xi represents the different variables (values in bold brackets) with x1 (initial lignin concentration, g/L), x2 (solvent flow rate, ml/min), x3 (antisolvent composition, water/EtOH, v/v), x4 (antisolvent ratio, solvent/antisolvent, v/v), and x5 (antisolvent stirring speed, rpm). (B) Summary of the different parameters used for the FFD model. xj represents the same variables as A) but with different values (in bold brackets) to maintain the self-assembly mechanism during precipitation. The specific values for each variable for the screening design construction, along with the corresponding results, are presented in Figure S5.
Figure 1.
(A) Summary of the different experimental parameters used for LNPs synthesis. xi represents the different variables (values in bold brackets) with x1 (initial lignin concentration, g/L), x2 (solvent flow rate, ml/min), x3 (antisolvent composition, water/EtOH, v/v), x4 (antisolvent ratio, solvent/antisolvent, v/v), and x5 (antisolvent stirring speed, rpm). (B) Summary of the different parameters used for the FFD model. xj represents the same variables as A) but with different values (in bold brackets) to maintain the self-assembly mechanism during precipitation. The specific values for each variable for the screening design construction, along with the corresponding results, are presented in Figure S5.
Figure 2.
Graphs from DLS show the effect of the different parameters on LNPs size and polydispersity index (PDI) with a) x1 (initial lignin concentration, g/L), b) x2 (solvent flow rate, ml/min), c) x3 (antisolvent composition, water/EtOH, v/v), d) x4 (antisolvent ratio, solvent/antisolvent, v/v), e), and x5 (antisolvent stirring speed, rpm). The precise particle size distribution from DLS is given in Figure S5. Results of the different parameters on LNPs zeta potential is given in Figure S6.
Figure 2.
Graphs from DLS show the effect of the different parameters on LNPs size and polydispersity index (PDI) with a) x1 (initial lignin concentration, g/L), b) x2 (solvent flow rate, ml/min), c) x3 (antisolvent composition, water/EtOH, v/v), d) x4 (antisolvent ratio, solvent/antisolvent, v/v), e), and x5 (antisolvent stirring speed, rpm). The precise particle size distribution from DLS is given in Figure S5. Results of the different parameters on LNPs zeta potential is given in Figure S6.
Figure 4.
Correlation between the predictive model and 26 experiments from the experimental design with the associated response (LNPs size average, nm). The model considers the following factors and their alias: x1 (initial lignin concentration, g/L), x2 (solvent flow rate, ml/min), x3 (antisolvent composition, water/EtOH, v/v), x4 (antisolvent ratio, solvent/antisolvent, v/v), and x5 (antisolvent stirring speed, rpm).
Figure 4.
Correlation between the predictive model and 26 experiments from the experimental design with the associated response (LNPs size average, nm). The model considers the following factors and their alias: x1 (initial lignin concentration, g/L), x2 (solvent flow rate, ml/min), x3 (antisolvent composition, water/EtOH, v/v), x4 (antisolvent ratio, solvent/antisolvent, v/v), and x5 (antisolvent stirring speed, rpm).
Table 1.
1 - Chemical composition of the raw biomass used in the study (results were taken from Girard et al. work [
44]). 2 - Main data for lignin isolation from 10 liters organosolv process. Extraction repetitions were performed from same dry biomass with equivalent granulometry. The initial particle diameter was 8 mm, and no acid catalysis was used. Cellulose, hemicellulose and lignin content were based on the solid residue’s analysis. Lignin macroparticles (LMPs) extraction yields (wtc %) are based on raw lignin biomass content. Each value is the mean of three repetitions.
Table 1.
1 - Chemical composition of the raw biomass used in the study (results were taken from Girard et al. work [
44]). 2 - Main data for lignin isolation from 10 liters organosolv process. Extraction repetitions were performed from same dry biomass with equivalent granulometry. The initial particle diameter was 8 mm, and no acid catalysis was used. Cellulose, hemicellulose and lignin content were based on the solid residue’s analysis. Lignin macroparticles (LMPs) extraction yields (wtc %) are based on raw lignin biomass content. Each value is the mean of three repetitions.
1. Raw Material (%) |
Cellulose |
Hemicellulose |
Lignin |
Extractives |
Ashes |
Total |
47.8 ± 1.5 |
22.5 ± 0.9 |
23.7 ± 0.2 |
2.7 ± 0.3 |
0.7 ± 0 |
97.4 ± 2.9 |
2. Organosolv solid residue (%) |
Cellulose |
Hemicellulose |
Lignin |
Mass loss yield (wt %) |
LMPs purity (%) |
LMPs isolation yields (wtc %) |
51.5 ± 1.8 |
13.0 ± 0.7 |
13.0 ± 0.4 |
46.0 ± 0.1 |
93.8 ± 0.3 |
70.8 ± 0.2 |
Table 2.
The experimental design summary shows the different factors, runs and the associated response (LNPs size average, nm). The factors are x1 (initial lignin concentration, g/L), x2 (solvent flow rate, ml/min), x3 (antisolvent composition, water/EtOH, v/v), x4 (antisolvent ratio, solvent/antisolvent, v/v), and x5 (antisolvent stirring speed, rpm).
Table 2.
The experimental design summary shows the different factors, runs and the associated response (LNPs size average, nm). The factors are x1 (initial lignin concentration, g/L), x2 (solvent flow rate, ml/min), x3 (antisolvent composition, water/EtOH, v/v), x4 (antisolvent ratio, solvent/antisolvent, v/v), and x5 (antisolvent stirring speed, rpm).
Run |
Design Factors (25-1) |
Design response |
x1 |
x2 |
x3 |
x4 |
x5 |
LNPs size average |
1 |
10 |
100 |
80 |
5 |
150 |
323 ± 6 |
2 |
20 |
2 |
100 |
5 |
1000 |
112 ± 3 |
3 |
10 |
100 |
80 |
20 |
1000 |
109 ± 2 |
4 |
20 |
2 |
100 |
20 |
150 |
114 ± 4 |
5 |
20 |
100 |
80 |
5 |
1000 |
266 ± 2 |
6 |
10 |
2 |
100 |
5 |
150 |
128 ± 3 |
7 |
20 |
100 |
80 |
20 |
150 |
173 ± 3 |
8 |
10 |
2 |
100 |
20 |
1000 |
50 ± 2 |
9 |
20 |
2 |
80 |
5 |
150 |
340 ± 6 |
10 |
10 |
100 |
100 |
5 |
1000 |
87 ± 3 |
11 |
20 |
2 |
80 |
20 |
1000 |
127 ± 3 |
12 |
10 |
100 |
100 |
20 |
150 |
62 ± 2 |
13 |
10 |
2 |
80 |
5 |
1000 |
237 ± 5 |
14 |
20 |
100 |
100 |
5 |
150 |
134 ± 3 |
15 |
10 |
2 |
80 |
20 |
150 |
153 ± 3 |
16 |
20 |
100 |
100 |
20 |
1000 |
73 ± 2 |
17 |
15 |
51 |
90 |
12.5 |
575 |
132 ± 2 |