硅酸盐学报

碳酸钙在受限空间内成核生长过程的分子动力学模拟

PDF下载

王攀 ¹ , 路兴海 ¹ , 王慕涵 ¹ , 张悦 ¹ , 许洪建 ² , 侯东帅 ¹

硅酸盐学报 | 研究论文 2024,52(2): 579-591

收起

硅酸盐学报 |研究论文 2024 , 52 (2) : 579 -591

碳酸钙在受限空间内成核生长过程的分子动力学模拟

全屏

王攀¹, 路兴海¹, 王慕涵¹, 张悦¹, 许洪建², 侯东帅¹

作者信息

^1．青岛理工大学土木工程学院，山东青岛 266033

^2．青岛第一市政工程有限公司，山东青岛 266033

通讯作者:

侯东帅（1986—），男，博士，教授。

作者简介:

王攀（1986—），男，博士，教授。

Molecular Dynamics Simulation of Nucleation and Growth of Calcium Carbonate in Confined Space

Pan WANG¹, Xinghai LU¹, Muhan WANG¹, Yue ZHANG¹, Hongjian XU², Dongshuai HOU¹

Affiliations

^1.School of Civil Engineering, Qingdao University of Technology, Qingdao 266033, Shandong, China

^2.Qingdao First Municipal Engineering Co., Ltd, Qingdao 266033, Shandong, China

Corresponding Author:

HOU Dongshuai (1986-), male, Ph.D., Professor. E-mail: dshou@outlook.com

WANG Pan (1986-), male, Ph.D., Professor. E-mail: wangpan@qut.edu.cn

出版时间: 2024-02-20 doi: 10.14062/j.issn.0454-5648.20230477

文章导航

摘要

收起

碳酸钙的形成会一定程度上改变混凝土的内部结构，进而影响混凝土的宏观性能，然而碳酸钙在混凝土内部成核生长机理尚未明确。为了探究碳酸钙在混凝土内部的形成过程和成核机制，本工作采用分子动力学模拟方法，研究了碳酸钙在不同受限空间中的形成过程和成核机理。结果表明：碳酸钙的成核生长主要依靠离子之间结合形成离子键；不同基体对碳酸钙团簇的吸附效果不同，这是因为碳酸钙与基体界面主要是钙氧离子键相连接，而水化硅酸钙基体界面处具有更多游离钙离子和羟基，能够为碳酸钙的形成提供更多的成核位点；同时，不同的受限空间孔径的大小对碳酸钙的形成速度也有较大影响。通过分析碳酸钙成核过程发现，碳酸钙团簇的形成不遵循经典成核理论，而是与预成核机制相符。本研究有助于进一步理解碳酸钙在混凝土中的形成过程，为混凝土中碳酸钙成核过程提供理论指导。

关键词

碳酸钙 / 成核生长 / 混凝土 / 受限空间 / 分子动力学模拟

Abstract

收起

Introduction

Carbon dioxide emissions generated during the production of cement-based materials are one of the important sources of global greenhouse gas emissions. Various technological innovations and improvement measures are introduced to reduce carbon emissions in cement industry. Among these, CO₂ can be fixed in concrete to form a stable calcium carbonate by solidifying CO₂ or accelerating carbonization, thereby reducing the release of CO₂ into the atmosphere. This carbonization process can effectively reduce the carbon dioxide produced by concrete and improve the comprehensive performance of concrete to a certain extent. However, the carbonization process is complicated due to the diverse components of concrete materials. The existing studies indicate that the components such as calcium silicate hydrate, ettringite, and calcium hydroxide in the concrete can participate in the carbonization reaction, and different components have different effects on the formation process of calcium carbonate. Therefore, this paper investigated the nucleation and growth process of calcium carbonate in confined space of different components by a molecular dynamics simulation method to reveal the formation mechanism of calcium carbonate in confined space. In addition, the structural characteristics of each component in concrete and the mechanism of their interaction with carbon dioxide were also analyzed.

Methods

The adsorption models of calcium carbonate in C-S-H, AFt, CH and SiO₂ confined space were established, respectively, to investigate the cluster formation process at the interface between calcium carbonate and different components of concrete. For C-S-H, the model includes C-S-H matrix and aqueous solution containing calcium carbonate. The C-S-H substrate with a Ca/Si ratio of 1.7 was used as a substrate model, and the supercell was cut along the crystal plane parallel to (001) to form the surface of the C-S-H substrate. The overall model size is 43.20 Å×45.04 Å×80.00 Å, and the pore size of the substrate is 4 nm. A total of 60 CO₃^2-, 60 Ca²⁺ and 2 000 water molecules were randomly placed in the pores, and the ion concentration in the solution was 1.4 mol/L. The other three matrices have the same size as the C-S-H matrix, and the concentration of Ca²⁺ and CO₃^2- in the pores is the same as that of the above C-S-H system. Also, four matrix models with confined space of 3 nm and 5 nm were established, respectively, to explore the influence of pore size on calcium carbonate. Before the formal simulation, the conjugate gradient method was used to minimize the energy, and then the isothermal isobaric ensemble (NPT) was used in the simulation process. The temperature was set to 300 K, the running time was 10 ns, and the time step was set to 1 fs. The clay force field (ClayFF) was used to simulate the C-S-H, AFt, CH, and SiO₂ matrix, respectively, and the carbonate solution used a separate force field parameter. The application of the mixed force field to different components in the model could more accurately simulate the dynamic characteristics and interactions between the components.

Results and discussion

In the case of bonding and local structure, the C-S-H, AFt and CH matrix form a stable and strong bonding structure between the interface region and the ions in the solution. Especially in the C-S-H system, Ca_c and CO₃^2- in the pores form more bonds with the ions at the interface of the matrix and have a stronger interaction. In contrast, the adsorption capacity of the ions in the pores is obviously weak, and only a small amount of ions are adsorbed due to the lack of calcium ions on the surface of the SiO₂ matrix.

The nucleation mechanism of calcium carbonate clusters is more consistent with the pre-nucleation theory. The whole simulation process can be divided into three stages, i.e., rapid agglomeration, cluster growth and cluster densification. The stronger interaction between the ions at the interface of the C-S-H matrix and the free ions in the solution lead to the aggregation of calcium carbonate clusters at the interface, providing more nucleation sites for the further nucleation and growth of calcium carbonate. Therefore, calcium carbonate has the optimum nucleation and agglomeration effect in the C-S-H matrix. Also, the pore size of the confined space has a great influence on the formation process of calcium carbonate. When the pore size is 3 nm, the surface of the C-S-H, CH and AFt matrix has a stronger attraction to the free ions in the solution, which is conducive to promoting the rapid nucleation and growth of calcium carbonate. When the pore size increases to 5 nm, the nucleation and growth process of calcium carbonate clusters slows down due to the weakened effect of the matrix interface on the ions in the solution.

Conclusions

The ions at the interface of the matrix had a great influence on the formation of calcium carbonate clusters. For different components of concrete, calcium carbonate and matrix surface ions attracted each other to form ionic bonds, which were then adsorbed at the interface for nucleation and growth. It was easier to bond with carbonate and calcium ions in the pores because C-S-H had more free hydroxyl groups and calcium ions at the interface, providing more nucleation sites for the agglomeration and growth of calcium carbonate at the interface. Also, the attraction of the substrate surface to the free ions in the solution was stronger when the pore size was 3 nm promoting the nucleation and growth of calcium carbonate. The nucleation mechanism of calcium carbonate clusters in the matrix pores did not follow the classical nucleation theory, but calcium ions and carbonate ions combined in the form of ionic bonds to form smaller clusters and then aggregated to grow and nucleate. In this process, it was unnecessary to overcome a large nucleation energy barrier, which was consistent with the pre-nucleation mechanism.

Key words

calcium carbonate / nucleation and growth / concrete / confined space / molecular dynamics simulation

引用本文

王攀, 路兴海, 王慕涵, 张悦, 许洪建, 侯东帅. 碳酸钙在受限空间内成核生长过程的分子动力学模拟. 硅酸盐学报, 2024 , 52 (2) : 579 -591 . DOI: 10.14062/j.issn.0454-5648.20230477

Pan WANG, Xinghai LU, Muhan WANG, Yue ZHANG, Hongjian XU, Dongshuai HOU. Molecular Dynamics Simulation of Nucleation and Growth of Calcium Carbonate in Confined Space[J]. Journal of the Chinese Ceramic Society, 2024 , 52 (2) : 579 -591 . DOI: 10.14062/j.issn.0454-5648.20230477

基金

收起

国家重点研发计划项目(2022YFE0133800)
国家自然科学基金(U2006224)
山东省自然科学基金(ZR2022YQ55)

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

2024年第52卷第2期

PDF下载

引用本文

BibTeX

文章信息

doi: 10.14062/j.issn.0454-5648.20230477

接收时间：2023-07-06
首发时间：2026-10-22
出版时间：2024-02-20

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2023-07-06
修回日期：2023-08-11

基金

国家重点研发计划项目(2022YFE0133800)

国家自然科学基金(U2006224)

山东省自然科学基金(ZR2022YQ55)

作者信息

^1．青岛理工大学土木工程学院，山东青岛 266033

^2．青岛第一市政工程有限公司，山东青岛 266033

通讯作者:

侯东帅（1986—），男，博士，教授。

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/gsyxb/CN/10.14062/j.issn.0454-5648.20230477

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT