nav emailalert searchbtn searchbox tablepage yinyongbenwen piczone journalimg journalInfo journalinfonormal searchdiv searchzone qikanlogo popupnotification paper paperNew
2024, 02, v.52 579-591
碳酸钙在受限空间内成核生长过程的分子动力学模拟
基金项目(Foundation): 国家重点研发计划项目(2022YFE0133800); 国家自然科学基金(U2006224); 山东省自然科学基金(ZR2022YQ55)
邮箱(Email): dshou@outlook.com;
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, CO2 can be fixed in concrete to form a stable calcium carbonate by solidifying CO2 or accelerating carbonization, thereby reducing the release of CO2 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 Si O2 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 CO32–, 60 Ca2+ 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 Ca2+ and CO32– 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(Clay FF) was used to simulate the C-S-H, AFt, CH, and Si O2 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, Cac and CO32– 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 Si O2 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.

参考文献

[1] SOUSA V, BOGAS J A, REAL S, et al. Industrial production of recycled cement:Energy consumption and carbon dioxide emission estimation[J]. Environ Sci Pollut Res, 2023, 30(4):8778–8789.

[2] MAKUL N. Advanced smart concrete-A review of current progress,benefits and challenges[J]. J Clean Prod, 2020, 274:122899.

[3] NIE S, ZHOU J, YANG F, et al. Analysis of theoretical carbon dioxide emissions from cement production:methodology and application[J]. J Clean Prod, 2022, 334:130270.

[4] LIU J H, WANG Y, LI Y Q, et al. Carbonated concrete brick capturing carbon dioxide from cement kiln exhaust gas[J]. Case Stud Constr Mater, 2022, 17:e01474.

[5] BARCELO L, KLINE J, WALENTA G, et al. Cement and carbon emissions[J]. Mater Struct, 2014, 47(6):1055–1065.

[6] WU M, ZHANG Y S, JI Y S, et al. Reducing environmental impacts and carbon emissions:Study of effects of superfine cement particles on blended cement containing high volume mineral admixtures[J]. J Clean Prod, 2018, 196:358–369.

[7] MADDALENA R, ROBERTS J J, HAMILTON A. Can Portland cement be replaced by low-carbon alternative materials? A study on the thermal properties and carbon emissions of innovative cements[J]. J Clean Prod, 2018, 186:933–942.

[8] YONG L, YU R, SHUI Z H, et al. Development of an environmental Ultra-High Performance Concrete(UHPC)incorporating carbonated recycled coarse aggregate[J]. Constr Build Mater, 2023, 362:129657.

[9] LIU J, ZHANG W Z, JIN H S, et al. Exploring the carbon capture and sequestration performance of biochar-artificial aggregate using a new method[J]. Scie Total Environ, 2023, 859:160423.

[10] GALAN I, ANDRADE C, MORA P, et al. Sequestration of CO2 by concrete carbonation[J]. Environ Sci Technol, 2010, 44(8):3181–3186.

[11] VOGLER N, DRABETZKI P, LINDEMANN M, et al. Description of the concrete carbonation process with adjusted depth-resolved thermogravimetric analysis[J]. J Therm Anal Calorim, 2022, 147(11):6167–6180.

[12] BORGES P H R, COSTA J O, MILESTONE N B, et al. Carbonation of CH and C–S–H in composite cement pastes containing high amounts of BFS[J]. Cem Concr Res, 2010, 40(2):284–292.

[13] MORANDEAU A, THIéRY M, DANGLA P. Investigation of the carbonation mechanism of CH and C-S-H in terms of kinetics,microstructure changes and moisture properties[J]. Cem Concr Res,2014, 56:153–170.

[14] CHANG J, FANG Y F. Quantitative analysis of accelerated carbonation products of the synthetic calcium silicate hydrate(C-S-H)by QXRD and TG/MS[J]. J Therm Anal Calorim, 2015, 119(1):57–62.

[15] WANG D, FANG Y F, ZHANG Y Y, et al. Changes in mineral composition, growth of calcite crystal, and promotion of physico-chemical properties induced by carbonation of β-C2S[J]. J CO2Util, 2019, 34:149–162.

[16] MARTíNEZ-RAMíREZ S, FERNáNDEZ-CARRASCO L. Carbonation of ternary cement systems[J]. Constr Build Mater, 2012, 27(1):313–318.

[17] LANGE L C, HILLS C D, POOLE A B. Preliminary investigation into the effects of carbonation on cement-solidified hazardous wastes[J].Environ Sci Technol, 1996, 30(1):25–30.

[18] CHANG J, LI Y, CAO M L, et al. Influence of magnesium hydroxide content and fineness on the carbonation of calcium hydroxide[J].Constr Build Mater, 2014, 55:82–88.

[19] SULAPHA P, WONG S F, WEE T H, et al. Carbonation of concrete containing mineral admixtures[J]. J Mater Civ Eng, 2003, 15(2):134–143.

[20] ZHAN B J, POON C S, LIU Q, et al. Experimental study on CO2curing for enhancement of recycled aggregate properties[J]. Constr Build Mater, 2014, 67:3–7.

[21] FANG Y F, CHANG J. Microstructure changes of waste hydrated cement paste induced by accelerated carbonation[J]. Constr Build Mater, 2015, 76:360–365.

[22] FENG H T, LI X, XING Y H, et al. Adsorption of CO32–/HCO3–on a quartz surface:Cluster formation, pH effects, and mechanistic aspects[J]. Phys Chem Chem Phys, 2023, 25(11):7951–7964.

[23] YIN B, XU T Y, HOU D S, et al. Superhydrophobic anticorrosive coating for concrete through in situ bionic induction and gradient mineralization[J]. Constr Build Mater, 2020, 257:119510.

[24] DEMICHELIS R, RAITERI P, GALE J D, et al. Stable prenucleation mineral clusters are liquid-like ionic polymers[J]. Nat Commun, 2011,2(1):1–8.

[25] HOU D S, MA H Y, YU Z, et al. Calcium silicate hydrate from dry to saturated state:Structure, dynamics and mechanical properties[J]. Acta Mater, 2014, 67:81–94.

[26] PELLENQ R J M, KUSHIMA A, SHAHSAVARI R, et al. A realistic molecular model of cement hydrates[J]. Proc Natl Acad Sci USA, 2009,106(38):16102–16107.

[27] QIN L, MAO X T, CUI Y F, et al. New insights into the early stage nucleation of calcium carbonate gels by reactive molecular dynamics simulations[J]. J Chem Phys, 2022, 157(23):9.

[28] CYGAN R T, LIANG J J, KALINICHEV A G. Molecular models of hydroxide, oxyhydroxide, and clay phases and the development of a general force field[J]. J Phys Chem B, 2004, 108(4):1255–1266.

[29] HOU D S, YANG Q R, WANG P, et al. Unraveling disadhesion mechanism of epoxy/CSH interface under aggressive conditions[J].Cem Concr Res, 2021, 146:106489.

[30] HOU D S, JIA Y T, YU J, et al. Transport properties of sulfate and chloride ions confined between calcium silicate hydrate surfaces:A molecular dynamics study[J]. J Phys Chem C, 2018, 122(49):28021–28032.

[31] RAITERI P, DEMICHELIS R, GALE J D. Thermodynamically consistent force field for molecular dynamics simulations of alkaline-earth carbonates and their aqueous speciation[J]. J Phys Chem C, 2015, 119(43):24447–24458.

[32] HEBERLING F, KLA?I?T, RAITERI P, et al. Structure and surface complexation at the calcite(104)–water interface[J]. Environ Sci Technol, 2021, 55(18):12403–12413.

[33] SCHUITEMAKER A, AUFORT J, KOZIARA K B, et al. Simulating the binding of key organic functional groups to aqueous calcium carbonate species[J]. Phys Chem Chem Phys, 2021, 23(48):27253–27265.

[34] REISCHL B, RAITERI P, GALE J D, et al. Atomistic simulation of atomic force microscopy imaging of hydration layers on calcite,dolomite, and magnesite surfaces[J]. J Phys Chem C, 2019, 123(24):14985–14992.

[35] HOSSEINI E, ZAKERTABRIZI M, HABIBNEJAD KORAYEM A, et al. Orbital overlapping through induction bonding overcomes the intrinsic delamination of 3D-printed cementitious binders[J]. ACS Nano, 2020, 14(8):9466–9477.

[36] DI TOMMASO D, DE LEEUW N H. The onset of calcium carbonate nucleation:A density functional theory molecular dynamics and hybrid microsolvation/continuum study[J]. J Phys Chem B, 2008, 112(23):6965–6975.

[37] AVARO J, MOON E M, ROSE J, et al. Calcium coordination environment in precursor species to calcium carbonate mineral formation[J]. Geochim Cosmochim Acta, 2019, 259:344–357.

[38] LOPEZ-BERGANZA J A, DIAO Y J, PAMIDIGHANTAM S, et al. Ab initio studies of calcium carbonate hydration[J]. J Phys Chem A, 2015,119(47):11591–11600.

[39] KELLERMEIER M, PICKER A, KEMPTER A, et al. A straightforward treatment of activity in aqueous CaCO3 solutions and the consequences for nucleation theory[J]. Adv Mater, 2014, 26(5):752–757.

[40] WANG B B, XIAO Y, XU Z M. Variation in properties of pre-nucleation calcium carbonate clusters induced by aggregation:A molecular dynamics study[J]. Crystals, 2021, 11(2):102.

[41] GEBAUER D, KELLERMEIER M, GALE J D, et al. Pre-nucleation clusters as solute precursors in crystallisation[J]. Chem Soc Rev, 2014,43(7):2348–2371.

[42] SHEN X Y, et al. New insights into the non-classical nucleation of C-S-H[J]. Cem Concr Res, 2023, 168:12.

[43] HOU D S, LI T, HAN Q H, et al. Insight on the sodium and chloride ions adsorption mechanism on the ettringite crystal:Structure,dynamics and interfacial interaction[J]. Comput Mater Sci, 2018, 153:479–492.

[44] QIAO G, HOU D S, WANG P, et al. Insights on failure modes of calcium-silicate-hydrate interface strengthened by polyacrylamides:Structure, dynamic and mechanical properties[J]. Constr Build Mater,2021, 278:122406.

[45] WANG P, YANG Q R, WANG M H, et al. Theoretical investigation of epoxy detachment from C-S-H interface under aggressive environment[J]. Constr Build Mater, 2020, 264:120232.

[46] SUN M, GENG G Q, XIN D B, et al. Molecular quantification of the decelerated dissolution of tri-calcium silicate(C3S)due to surface adsorption[J]. Cem Concr Res, 2022, 152:106682.

基本信息:

DOI:10.14062/j.issn.0454-5648.20230477

中图分类号:TQ132.32

引用信息:

[1]王攀,路兴海,王慕涵,等.碳酸钙在受限空间内成核生长过程的分子动力学模拟[J].硅酸盐学报,2024,52(02):579-591.DOI:10.14062/j.issn.0454-5648.20230477.

基金信息:

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

发布时间:

2023-12-02

出版时间:

2023-12-02

网络发布时间:

2023-12-02

检 索 高级检索

引用

GB/T 7714-2015 格式引文
MLA格式引文
APA格式引文