Piaggi Pablo M, Selloni Annabella, Panagiotopoulos Athanassios Z, Car Roberto, Debenedetti Pablo G
Department of Chemistry, Princeton University, Princeton, NJ 08544, USA.
Department of Chemical and Biological Engineering, Princeton University, Princeton, NJ 08544, USA.
Faraday Discuss. 2024 Feb 6;249(0):98-113. doi: 10.1039/d3fd00100h.
The formation of ice in the atmosphere affects precipitation and cloud properties, and plays a key role in the climate of our planet. Although ice can form directly from liquid water under deeply supercooled conditions, the presence of foreign particles can aid ice formation at much warmer temperatures. Over the past decade, experiments have highlighted the remarkable efficiency of feldspar minerals as ice nuclei compared to other particles present in the atmosphere. However, the exact mechanism of ice formation on feldspar surfaces has yet to be fully understood. Here, we develop a first-principles machine-learning model for the potential energy surface aimed at studying ice nucleation at microcline feldspar surfaces. The model is able to reproduce with high-fidelity the energies and forces derived from density-functional theory (DFT) based on the SCAN exchange and correlation functional. Our training set includes configurations of bulk supercooled water, hexagonal and cubic ice, microcline, and fully-hydroxylated feldspar surfaces exposed to a vacuum, liquid water, and ice. We apply the machine-learning force field to study different fully-hydroxylated terminations of the (100), (010), and (001) surfaces of microcline exposed to a vacuum. Our calculations suggest that terminations that do not minimize the number of broken bonds are preferred in a vacuum. We also study the structure of supercooled liquid water in contact with microcline surfaces, and find that water density correlations extend up to around 10 Å from the surfaces. Finally, we show that the force field maintains a high accuracy during the simulation of ice formation at microcline surfaces, even for large systems of around 30 000 atoms. Future work will be directed towards the calculation of nucleation free-energy barriers and rates using the force field developed herein, and understanding the role of different microcline surfaces in ice nucleation.
大气中冰的形成会影响降水和云的特性,并在我们星球的气候中发挥关键作用。尽管在深度过冷条件下冰可以直接从液态水形成,但外来颗粒的存在可以在温度高得多的情况下促进冰的形成。在过去十年中,实验突出了长石矿物作为冰核相对于大气中存在的其他颗粒的显著效率。然而,长石表面上冰形成的确切机制尚未完全理解。在这里,我们开发了一种针对势能面的第一性原理机器学习模型,旨在研究微斜长石表面的冰核形成。该模型能够以高保真度重现基于SCAN交换和关联泛函的密度泛函理论(DFT)得出的能量和力。我们的训练集包括大块过冷水、六方冰和立方冰、微斜长石以及暴露于真空、液态水和冰的完全羟基化长石表面的构型。我们应用机器学习力场来研究暴露于真空的微斜长石(100)、(010)和(001)表面的不同完全羟基化终端。我们的计算表明,在真空中,不使断键数量最小化的终端更受青睐。我们还研究了与微斜长石表面接触的过冷液态水的结构,发现水密度相关性从表面延伸到约10 Å左右。最后,我们表明,即使对于约30000个原子的大系统,该力场在微斜长石表面冰形成的模拟过程中也能保持高精度。未来的工作将致力于使用本文开发的力场计算成核自由能垒和速率,并理解不同微斜长石表面在冰核形成中的作用。