Lignin and cellulose are found ubiquitously in plant cell walls, making them two of the most abundant biopolymers on Earth. In spite of their ubiquity, comparatively little is known quantitatively about their interactions within native biomass, and by extension how their interaction may play a role in biomaterials. The possible binding faces for lignin on crystalline cellulose I-ss and the heterogeneous chemistries found in native lignin typically preclude the isolation of specific interactions for quantitative comparison in the laboratory. In this study, we use molecular simulation to examine individual combinations of binding face and lignin chemistry to quantify lignin-cellulose interactions, including contributions from both specific hydrogen bonds and nonspecific van der Waals interactions. For all lignin chemistries simulated, the hydrophobic 200 face is the preferred for cellulose interactions. For hydrophilic crystalline faces, lignin is found to bind most effectively to the 110 face. The affinity between lignin and cellulose increases as additional lignin methoxy groups are added. These groups increase the lignin-cellulose contact area, increasing the binding affinity. These trends are extended to larger lignin polymers, where the relationships between lignin polymer size and binding affinity are quantified through biased simulation. In this case, the binding surface area between the lignin polymer and the cellulose surface dictates the interaction strength.