Continuous melting through a hexatic phase in confined bilayer water

Liquid water is not only of obvious importance but also extremely intriguing, displaying many anomalies that still challenge our understanding of such an a priori simple system. The same is true when looking at nanoconfined water: The liquid between constituents in a cell is confined to such dimensions, and there is already evidence that such water can behave very differently from its bulk counterpart. A striking finding has been reported from computer simulations for two-dimensionally confined water: The liquid displays continuous or discontinuous melting depending on its density. In order to understand this behavior, we have analyzed the melting exhibited by a bilayer of nanoconfined water by means of molecular dynamics simulations. At high density we observe the continuous melting to be related to the phase change of the oxygens only, with the hydrogens remaining liquidlike throughout. Moreover, we find an intermediate hexatic phase for the oxygens between the liquid and a triangular solid ice phase, following the Kosterlitz-Thouless-Halperin-Nelson-Young theory for two-dimensional melting. The liquid itself tends to maintain the local structure of the triangular ice, with its two layers being strongly correlated yet with very slow exchange of matter. The decoupling in the behavior of the oxygens and hydrogens gives rise to a regime in which the complexity of water seems to disappear, resulting in what resembles a simple monoatomic liquid. This intrinsic tendency of our simulated water may be useful for understanding novel behaviors in other confined and interfacial water systems

Enhanced Configurational Entropy in High-Density Nanoconfined Bilayer Ice

A novel kind of crystal order in high-density nanoconfined bilayer ice is proposed from molecular dynamics and density-functional theory simulations. A first-order transition is observed between a low-temperature proton-ordered solid and a high-temperature proton-disordered solid. The latter is shown to possess crystalline order for the oxygen positions, arranged on a close-packed triangular lattice with AA stacking. Uniquely among the ice phases, the triangular bilayer is characterized by two levels of disorder (for the bonding network and for the protons) which results in a configurational entropy twice that of bulk ice.This work was partly funded by Grants No. FIS2012- 37549-C05 from the Spanish Ministry of Science, and Exp. 97/14 (Wet Nanoscopy) from the Programa Red Guipuzcoana de Ciencia, Tecnología e Innovación, Diputación Foral de Gipuzkoa. We thank José M. Soler and M.-V. Fernández-Serra for useful discussions. The calculations were performed on the arina HPC cluster (Universidad del País Vasco/Euskal Herriko Unibertsitatea, Spain). SGIker (UPV/EHU, MICINN, GV/EJ, ERDF and ESF) support is gratefully acknowledged

Knock-on damage in bilayer graphene: Indications for a catalytic pathway

We study by high-resolution transmission electron microscopy the structural response of bilayer graphene to electron irradiation with energies below the knock-on damage threshold of graphene. We observe that one type of divacancy, which we refer to as the butterfly defect, is formed for radiation energies and doses for which no vacancies are formed in clean monolayer graphene. By using first principles calculations based on density-functional theory, we analyze two possible causes related with the presence of a second layer that could explain the observed phenomenon: an increase of the defect stability or a catalytic effect during its creation. For the former, the obtained formation energies of the defect in monolayer and bilayer systems show that the change in stability is negligible. For the latter, ab initio molecular dynamics simulations indicate that the threshold energy for direct expulsion does not decrease in bilayer graphene as compared with monolayer graphene, and we demonstrate the possibility of creating divacancies through catalyzed intermediate states below this threshold energy. The estimated cross section agrees with what is observed experimentally. Therefore, we show the possibility of a catalytic pathway for creating vacancies under electron radiation below the expulsion threshold energy. © 2013 American Physical Society