Structure and dynamics of confined water in naphthalene-diimide based molecular crystals

Under confinement, water changes its structure and dynamics, displaying new properties with respect to the bulk. We studied water confined in naphthalene diimide (NDI)-based molecular crystals via classical molecular dynamics simulations. We examined NDIs functionalized with either hydrophilic linear triethylene glycol side chains (NDI-TEG) or hydrophobic n-hexyl chains (NDI-C6), increasing the hydration (i.e., amount of water molecules within each crystal) up to a 1:3 NDI:water molar ratio. Static and dynamical properties of confined water are compared to those of room-temperature and supercooled (200 K) bulk water, which serve as references for liquid and glassy states, respectively. The impact of confinement is analyzed through structural order parameters, time-dependent correlation functions, and hydrogen-bond (HB) analyses. At low hydration (approximate to 1:1 NDI:water), confined water assumes a structural order that is far from the tetrahedral one, showing different spatial organizations when inserted within hydrophilic NDI-TEG or hydrophobic NDI-C6 crystals. By increasing the hydration level (approximate to 1:3 NDI:water), the structure of confined water clusters shifts toward a more bulk-like (liquid) arrangement, while the chemical nature of the NDI side chains plays a marginal role. Notably, fast (sub-ps) and slow (hundreds of ns) water dynamics are not much influenced by the hydrophilicity/hydrophobicity of the side chains but rather by confinement effects. The analysis of the HB network autocorrelation functions highlights how finite-size effects and restricted connectivity are the main factors controlling HB dynamics. Our study paves the way toward an atomistic understanding of both structural and dynamical functions of water confined in molecular crystals, opening new paths for the rationalization of transport phenomena in the emerging class of organic mixed ionic-electronic conductors.