Diversity of potential hydrogen bonds in cellulose I revealed by molecular dynamics simulation