Molecular dynamics simulation based on the multi-component molecular orbital method: Application to H5O2+,D5O2+,andT5O2+

► Molecular dynamics method with nuclear quantum effect was developed. ► Multi-component molecular orbital method was used as ab initio MO calculation. ► Developed method applied to basic hydrogen bonding system, H5O2+, and isotopomers. ► O⋯O vibrational stretching reflected to the distribution of protonic wavefunctions. ► H/D/T isotope effect was also analyzed.