Dipartimento di Ingegneria Meccanica e Aerospaziale, Sapienza Università di Roma, Rome, Italy.
J Chem Phys. 2023 May 28;158(20). doi: 10.1063/5.0147647.
Understanding intrusion and extrusion in nanoporous materials is a challenging multiscale problem of utmost importance for applications ranging from energy storage and dissipation to water desalination and hydrophobic gating in ion channels. Including atomistic details in simulations is required to predict the overall behavior of such systems because the statics and dynamics of these processes depend sensitively on microscopic features of the pore, such as the surface hydrophobicity, geometry, and charge distribution, and on the composition of the liquid. On the other hand, the transitions between the filled (intruded) and empty (extruded) states are rare events that often require long simulation times, which are difficult to achieve with standard atomistic simulations. In this work, we explored the intrusion and extrusion processes using a multiscale approach in which the atomistic details of the system, extracted from molecular dynamics simulations, informed a simple Langevin model of water intrusion/extrusion in the pore. We then used the Langevin simulations to compute the transition times at different pressures, validating our coarse-grained model by comparing it with nonequilibrium molecular dynamics simulations. The proposed approach reproduces experimentally relevant features such as the time and temperature dependence of the intrusion/extrusion cycles, as well as specific details about the shape of the cycle. This approach also drastically increases the timescales that can be simulated, reducing the gap between simulations and experiments and showing promise for more complex systems.
理解纳米多孔材料中的侵入和挤出是一个具有挑战性的多尺度问题,对于从能量存储和耗散到水淡化和离子通道中的疏水性门控等应用至关重要。为了预测这些系统的整体行为,需要在模拟中包含原子细节,因为这些过程的静态和动态特性敏感地依赖于孔的微观特征,如表面疏水性、几何形状和电荷分布,以及液体的组成。另一方面,填充(侵入)和空(挤出)状态之间的转变是罕见事件,通常需要很长的模拟时间,这在标准的原子模拟中很难实现。在这项工作中,我们使用多尺度方法探索了侵入和挤出过程,该方法从分子动力学模拟中提取系统的原子细节,为水在孔中的侵入/挤出的简单朗之万模型提供信息。然后,我们使用朗之万模拟在不同压力下计算跃迁时间,通过将其与非平衡分子动力学模拟进行比较来验证我们的粗粒化模型。所提出的方法再现了实验相关的特征,例如侵入/挤出循环的时间和温度依赖性,以及循环形状的具体细节。这种方法还极大地增加了可以模拟的时间尺度,缩小了模拟和实验之间的差距,并为更复杂的系统展示了前景。