Born-Oppenheimer molecular dynamics simulations on structures of high-density and low-density water: a comparison of the SCAN meta-GGA and PBE GGA functionals.

Using Born-Oppenheimer ab initio molecular dynamics (BOAIMD) simulations, the high-density water (HDW) and low-density water (LDW) structures based on SCAN meta-GGA are compared with those based on PBE GGA. Compared with Car-Parrinello ab initio molecular dynamics (CPAIMD) simulations, BOAIMD simulations can produce more accurate results because no fictitious electron mass is introduced. At each state point, our simulations continue for 100 ps after the system reached equilibrium, which is the longest ab initio simulations of liquid water reported so far and can ensure an accurate statistical average.