The doublet C4 peaks at * 85 and * 89 ppm in solid-state 13 C NMR spectra of native cellulose have been attributed to signals of C4 atoms on the surface (solvent-exposed) and in the interior of microfibrils, designated as sC4 and iC4, respectively. The relative intensity ratios of sC4 and iC4 observed in NMR spectra of cellulose have been used to estimate the degree of crystallinity of cellulose and the number of glucan chains in cellulose microfibrils. However, the molecular structures of cellulose responsible for the specific surface and interior C4 peaks have not been positively confirmed. Using density functional theory (DFT) methods and structures produced from classical molecular dynamics simulations, we investigated how the following four factors affect 13 C NMR chemical shifts in cellulose: conformations of exocyclic groups at C6 (tg, gt and gg), H 2 O molecules H-bonded on the surface of the microfibril, glycosidic bond angles (U, W) and the distances between H4 and HO3 atoms. We focus on changes in the d 13 C4 value because it is the most significant observable for the same C atom within the cellulose structure. DFT results indicate that different conformations of the exocyclic groups at C6 have the greatest influence on d 13 C4 peak separation, while the other three factors have secondary effects that increase the spread of the calculated C4 interior and surface peaks.