https://www.selleckchem.com/products/rgt-018.html itate the development of future force fields based on the pGM electrostatics for applications in biomolecular systems and processes where polarization plays crucial roles.We report a molecular dynamics simulation study of dense ice modeled by the reactive force field (ReaxFF) potential, focusing on the possibility of phase changes between crystalline and plastic phases as observed in earlier simulation studies with rigid water models. It is demonstrated that the present model system exhibits phase transitions, or crossovers, among ice VII and two plastic ices with face-centered cubic (fcc) and body-centered cubic (bcc) lattice structures. The phase diagram derived from the ReaxFF potential is different from those of the rigid water models in that the bcc plastic phase lies on the high-pressure side of ice VII and does the fcc plastic phase on the low-pressure side of ice VII. The phase boundary between the fcc and bcc plastic phases on the pressure, temperature plane extends to the high-temperature region from the triple point of ice VII, fcc plastic, and bcc plastic phases. Proton hopping, i.e., delocalization of a proton, along between two neighboring oxygen atoms in dense ice is observed for the ReaxFF potential but only at pressures and temperatures both much higher than those at which ice VII-plastic ice transitions are observed.We implement Epstein-Nesbet perturbative corrections in the third-order for the initiator approximation of the configuration space quantum Monte Carlo. An improved sampling algorithm is proposed to address the stochastic noise of the corrections. The stochastic error for the perturbative corrections is considerably larger than that for the reference energy, and it fails to provide reasonable results unless a very long imaginary time integration is performed. The new sampling algorithm accumulates rejected walkers from multiple independent steps to cover a larger portion of the secondary