CCES Unicamp

Uniaxial-deformation behavior of ice I h as described by the TIP4P/Ice and mW water models

Using molecular dynamics simulations, we assess the uniaxial deformation response of ice I h as described by two popular water models, namely, the all-atom TIP4P/Ice potential and the coarse-grained mW model. In particular, we investigate the response to both tensile and compressive uniaxial deformations along the [0001] and [0¯110] crystallographic directions for a series of different temperatures. We classify the respective failure mechanisms and assess their sensitivity to strain rate and cell size. While the TIP4P/Ice model fails by either brittle cleavage under tension at low temperatures or large-scale amorphization/melting, the mW potential behaves in a much more ductile manner, displaying numerous cases in which stress relief involves the nucleation and subsequent activity of lattice dislocations. Indeed, the fact that mW behaves in such a malleable manner even at strain rates that are substantially higher than those applied in typical experiments indicates that the mW description of ice I h is excessively ductile. One possible contribution to this enhanced malleability is the absence of explicit protons in the mW model, disregarding the fundamental asymmetry of the hydrogen bond that plays an important role in the nucleation and motion of lattice dislocations in ice I h .

J. Chem. Phys. 149 164711 (2018)

https://aip.scitation.org/doi/10.1063/1.5048517

Related posts

Insights on the mechanism of water-alcohol separation in multilayer graphene oxide membranes: Entropic versus enthalpic factors

escience

Analysis of three-dimensional hexagonal and cubic polycrystals using the boundary element method

cces cces

CCES researchers accelerate scientific cloud computing

cces cces