Why different water models predict different structures under 2D confinement

James Dix, Leo Lue, Paola Carbone

    Research output: Contribution to journalArticlepeer-review


    Experiments of nanoconfined water between graphene sheets at high pressure suggest that it forms a square ice structure (Algara‐Siller et al., Nature, 2015, 519, 443). Molecular dynamics (MD) simulations have been used to attempt to recreate this structure, but there have been discrepancies in the structure formed by the confined water depending on the simulation set‐up that was employed and particularly on the choice of water model. Here, using classical molecular dynamics simulations, we have systematically investigated the effect that three different water models (SPC/E, TIP4P/2005 and TIP5P) have on the structure of water confined between two rigid graphene sheets with a 0.9 nm separation. We show that the TIP4P/2005 and the TIP5P water models form a hexagonal AA‐stacked structure, whereas the SPC/E model forms a rhombic AB‐stacked structure. Our work demonstrates that the formation of these structures is driven by differences in the strength of hydrogen bonds predicted by the three water models, and that the nature of the graphene/water interaction only mildly affects the phase diagram. Considering the available experimental data and first‐principle simulations we conclude that, among the models tested, the TIP4P/2005 and TIP5P force fields are for now the most reliable when simulating water under confinement.
    Original languageEnglish
    Pages (from-to)2051-2059
    Number of pages9
    JournalJournal of Computational Chemistry
    Issue number25
    Early online date2018
    Publication statusPublished - 30 Sept 2018


    • SPC/E
    • TIP4P/2005
    • TIP5P
    • confinement
    • graphene
    • hydrogen bond
    • molecular dynamics
    • water


    Dive into the research topics of 'Why different water models predict different structures under 2D confinement'. Together they form a unique fingerprint.

    Cite this