Investigating Lignin Aggregation and Interactions with Solvents during γ-Valerolactone (GVL) Pretreatment: A Combined Small Angle Neutron Scattering and Molecular Simulations Study

The strong tendency of lignin to aggregate in solution, coupled with limited understanding of how its molecular structure governs this behavior, hinders its effective utilization in biorefineries. Here, we investigated the solution behavior of lignin extracted from poplar using γ-valerolactone/water (GVL/H2O, 9:1 wt/wt) through combined small-angle neutron scattering (SANS) and molecular dynamics (MD) simulations. Lignin samples obtained at 100 °C (L100) and 120 °C (L120) differed in β–O–4 content, hydroxyl distribution, and S/G ratio, enabling direct assessment of how molecular composition governs solvation and aggregation. SANS showed that L120 formed rigid and elongated cylindrical aggregates at 25 °C that transitioned to more flexible spheroidal structures by 50 °C and remained stable up to 80 °C, whereas L100 adopted globular aggregates that progressively collapsed with increasing temperature. MD simulations reinforced these observations by showing that S-rich (L120-like) oligomers had larger radii of gyration, stronger solvent coordination driven by methoxy groups, and fewer lignin–lignin contacts. In contrast, G-rich (L100-like) oligomers displayed persistent aggregation and lower solubility. Collectively, these results reveal that increased aromatic methoxylation enhances lignin–solvent interactions and suppresses self-association, whereas reduced methoxylation and higher β–O–4 content promote persistent aggregation into colloid-like structures with restricted solvent penetration into the aggregate interior.