Reactive molecular dynamics (MD) is used to simulate the equilibrium process of water confined between two fully hydroxylatedα-quartz (001) surfaces with separation distances from 7 to 20 ?. Effect of different patterns of interfacial hydrogen bonds on the structure and dynamics of confined water is investigated. Density profiles, radial distribution functions, number of interfacial hydrogen bonds, and mean square displacements are calculated. The α-quartz (001) surface is cut from an α-quartz crystal at a certain depth to construct a surface with geminal silanols after being fully hydroxylated. The silanol groups on the surface are treated in two different ways in the MD simulations. One of the silanol groups are treated as to be fixed, and the other one is treated as no constraint for the movement of surface silanols. Our results show that different patterns of hydrogen bonds are formed at the interface between SiO2 surface and water. For the fixed silanol surface there is one type of strong hydrogen bonds interacting between the oxygen atoms of water and the hydrogen atoms of surface silanols, leading to the dipole moment of water molecules pointing out from the surface. For the movable silanol surface there are two types of strong hydrogen bonds formed at the interface. One is between the oxygen atoms of water and the hydrogen atoms of surface silanols, and the other is between the oxygen atoms of surface silanols and the hydrogen atoms of water. The number of hydrogen bonds of the first type is much less than those of the second type, leading to the dipole moment of water molecules pointing to the surface. Moreover, the total number of interfacial hydrogen bonds formed on the fixed silanol surfaces is larger than that on the movable silanol surfaces. The density profiles of the confined water indicate the formation of a strong layering of water in the vicinity of the fixed silanol surface, and the water layer is also more ordered with an ice-like structure, as compared with a dense water layer with a liquid-like structure in the case of movable silanol surfaces. Thus the mean square displacements of confined water show that, as compared with interfacial hydrogen bonds formed on the fixed silanol surfaces, the weaker and the lesser interfacial hydrogen bonds formed on the movable silanol surfaces may be responsible for more intense movement of confined water between the movable silanol surfaces. Our simulation suggests that the different pattern of interfacial hydrogen bonds could signifiantly affect the structure and dynamic behaviors of the confined water between two fully hydroxylated silica surfaces.
受限水 界面氢键 二氧化硅晶体 分子动力学模拟