Crystallisation Competition Between Cubic And Hexagonal Ice Structures: Molecular-Dynamics Insight

Under ambient pressure, two different ice polytypes of ice, cubic (I-c ) and hexagonal (I-h ) can be crystallised. Here, we elucidate the mechanism and kinetics of the phase transition of cubic ice with no stacking disorder to a perfect hexagonal ice at various temperatures, via molecular simulation. The study was performed through non-isothermal and isothermal molecular-dynamics simulations; these conditions were enforced on an equilibrated I-c /LDA/I-h system. Different results have been drawn from analysis in terms of local density, annealing kinetics, and phase transition. The evolution from cubic to hexagonal ice becomes prevalent at higher temperatures during the performed non-isothermal studies; to this end, a model-free kinetic study was proposed to address this phenomenon's kinetic properties over a characteristic temperature range. Ice crystallisation is also studied in solid/liquid systems, where cubic or hexagonal ice was in contact with liquid water at sub-zero temperatures. It was revealed that water molecules have a propensity to follow and crystalise in the existing ice structure rather than forming a new interface.