硅酸盐学报

纳米纤维素改性水泥基材料的分子动力学模拟

PDF下载

林子韵 ¹ , 刘巧玲 ¹^,² , 孙晓宇 ¹ , 张一晨 ¹ , 任勇 ¹ , 苏元鑫 ¹ , 王桢 ¹^,²

硅酸盐学报 | 研究论文 2026,54(10): 3628-3637

收起

硅酸盐学报 |研究论文 2026 , 54 (10) : 3628 -3637

纳米纤维素改性水泥基材料的分子动力学模拟

全屏

林子韵¹, 刘巧玲¹^,², 孙晓宇¹, 张一晨¹, 任勇¹, 苏元鑫¹, 王桢¹^,²

作者信息

^1.山东建筑大学土木工程学院，济南　250101

^2.山东建筑大学建筑结构加固改进与地下空间工程教育部重点实验室，济南　250101

通讯作者:

刘巧玲（1977-），女，教授

王桢（1991-），女，讲师。

作者简介:

林子韵（2001-），女，硕士研究生。

Molecular Dynamics Simulation of Nanocellulose-Modified Cement-Based Materials

Ziyun LIN¹, Qiaoling LIU¹^,², Xiaoyu SUN¹, Yichen ZHANG¹, Yong REN¹, Yuanxin SU¹, Zhen WANG¹^,²

Affiliations

^1.School of Civil Engineering, Shandong Jianzhu University, Jinan 250101, China

^2.Key Laboratory of Building Structural Retrofitting and Underground Space Engineering, Ministry of Education, Shandong Jianzhu University, Jinan 250101, China

Corresponding Author:

LIU Qiaoling (1977-), female, Ph.D., Professor. E-mail: lql263@163.com

WANG Zhen (1991-), female, Lecturer. zhenwang25@126.com

LIN Ziyun (2001-), female, Master candidate. E-mail: 308550692@qq.com

出版时间: 2026-07-29 doi: 10.14062/j.issn.0454-5648.20250755

文章导航

摘要

收起

纳米纤维素对水泥基材料性能的优化常因界面作用机制不明确而受限。采用分子动力学模拟方法，研究了水化硅酸钙（C-S-H）与羧基化改性前后的纳米纤维素之间的界面结构、动力学行为、能量特性与力学性能。通过将纳米纤维素嵌入C-S-H基体中，构建了具有不同表面官能团的复合模型。模拟结果表明，羧基化纳米纤维素中的-COOH基团和羟基化纳米纤维素中的-OH基团，为C-S-H提供了大量氢键位点。界面结合能大小与官能团类型密切相关，羧基化纳米纤维素的界面结合能明显高于羟基化纳米纤维素。此外，羧基化纤维素中更多的非桥氧位点可与C-S-H中的钙离子形成更稳定的配位结构，并构建更密的氢键网络。本工作从微观层面揭示了纳米纤维素对C-S-H的增强机制，为设计高性能的绿色水泥基复合材料提供了理论支持与选材依据。

关键词

纳米纤维素 / 分子动力学模拟 / 水泥基材料 / 水化硅酸钙凝胶

Abstract

收起

Introduction

The application of nanomaterials in cement-based composites has attracted considerable research interest in recent years, because of their ability to enhance mechanical properties such as compressive, flexural, and fracture strength. Among these materials, nanocellulose, a bio-derived nanomaterial obtained from abundant renewable resources and characterized by excellent properties, has considerable potential for use in cement composites. Its surface functional groups can act as nucleation sites for calcium silicate hydrate (C-S-H) gel, leading to the formation of denser hydration products and refined pore structures. Furthermore, nanocellulose inhibits microcrack propagation through nano-filling and bridging effects, thereby synergistically improving the mechanical performance of cement-based materials. Despite these advances, the interfacial bonding behavior between nanocellulose and C-S-H gel at the nanoscale is not well understood because of limitations in conventional macroscopic experiments and microscopic characterization techniques. The reinforcing effect of nanocellulose has been proposed to arise primarily from a dual micromechanism of "short-circuit diffusion and bridging crack resistance". However, direct atomistic evidence supporting this hypothesis is lacking. In particular, the influence of nanocellulose surface functionalization (e.g., carboxylation) on interfacial interactions has not been systematically investigated. In this study, molecular dynamics (MD) simulations were used to investigate the interfacial structure, interaction energy, hydrogen bonding evolution, and failure behaviors of nanocellulose/C-S-H composites. A novel embedded composite model was constructed using a Monte Carlo-based strategy. The effects of hydroxylated and carboxylated nanocellulose on interfacial binding and mechanical performance were systematically compared.

Methods

The C-S-H gel model with a chemical composition of (CaO)_1.65(SiO₂)(H₂O)_1.75 was constructed following the realistic disordered structure proposed by Pellenq et al., which matches the experimental density (~2.6 g/cm³) and Ca/Si ratio (1.6-1.8). The modeling process started from the 11 Å tobermorite crystal structure, followed by random removal of bridging silicate tetrahedra, water saturation at 298 K, and subsequent energy minimization and NPT dynamic relaxation (298 K, 120 ps, 1 fs timestep) using the COMPASS force field. C-S-H models with different Ca/Si ratios (0.75, 1.31, 1.48, and 1.67) were also constructed by controlling the removal of silicate tetrahedra. Amorphous cellulose models (both hydroxylated and carboxylated) were built using cellobiose as the repeating unit. A three-dimensional periodic amorphous cellulose model (43.0 Å × 43.0 Å × 22.0 Å) was constructed by randomly packing five cellulose chains using the Grand Canonical Monte Carlo (GCMC) method, followed by NVT relaxation (298 K, 100 ps, 1 fs timestep). The final density of 1.32 g/cm³ agreed well with experimental values. The C-S-H/cellulose composite model was constructed by inserting the amorphous cellulose box into the interlayer of C-S-H, with a 5 Å vacuum layer to avoid initial atomic overlap. The system was first energy-minimized using the smart algorithm (energy tolerance 1.0×10^-4 kcal/mol, force tolerance 0.001 kcal/mol/Å), followed by NVT equilibration (298 K, 100 ps, 1 fs timestep). Uniaxial tensile simulations were performed along X, Y, and Z directions to obtain stress-strain curves, Young's modulus and tensile strength. Interfacial interaction energy was calculated as E_inter = E_{_total}-(E_{_C-S-H} + E_{_CNF}). Hydrogen bonding, radial distribution functions (RDF) and mean squared displacement (MSD) were also analyzed.

Results and discussion

With increasing Ca/Si ratio from 0.75 to 1.67, the degree of polymerization of silicate chains decreased systematically, as evidenced by the decrease in Q₂ content from 100% to 20%. The C-S-H model with Ca/Si = 1.67 yielded a density of 2.44 g/cm³ and Q_n distribution (Q₀ = 11.6%, Q₁ = 65.12%, Q₂ = 20%), which closely matched the real molecular structure of C-S-H gel. Uniaxial tensile simulations revealed that Young's modulus decreased with increasing Ca/Si ratio in all three directions (X, Y, Z). For Ca/Si = 1.67, Young's modulus values ranged between 10-30 GPa, consistent with nanoindentation experimental results. The Y-direction exhibited the highest modulus due to the dominant contribution of silicate chains, while the X and Z directions were governed by the fractured hydrogen bond network. The interfacial interaction energies for both hydroxylated and carboxylated nanocellulose with C-S-H were negative, indicating spontaneous adsorption and thermodynamically stable interfaces. The absolute value of interfacial interaction energy for carboxylated nanocellulose/C-S-H (-1.70×10² kJ/mol) was significantly larger than that for hydroxylated nanocellulose/C-S-H (-0.73×10² kJ/mol), demonstrating stronger spontaneous adsorption and tighter interfacial bonding after carboxylation. The interaction energy mainly consisted of van der Waals and electrostatic contributions, with electrostatic energy accounting for over 90% of the total adsorption energy. Hydrogen bond analysis showed that in the C-S-H/carboxylated nanocellulose composite model, 56.1 hydrogen bonds were formed in the interfacial region, corresponding to a hydrogen bond density of 3.12 bonds/nm². In contrast, the hydroxylated nanocellulose/C-S-H interface gave only 25.7 hydrogen bonds with a density of 1.43 bonds/nm². The significantly higher hydrogen bond density after carboxylation strongly enhanced interfacial binding. RDF analysis revealed characteristic peaks at 0.97 Å and 2.05 Å for oxygen and hydrogen atoms in C-S-H, with atomic distances smaller than the hydrogen bond cutoff distance of 2.5 Å. The carboxylated system exhibited more pronounced peaks, indicating stronger bonding ability of carboxyl oxygen with interfacial hydrogen. MSD analysis showed that atoms at the interface of the carboxylated nanocellulose/C-S-H composite had lower MSD values than those in the hydroxylated system, indicating stronger atomic constraints and lower motional freedom. This enhanced interfacial stability enabled more effective stress transfer from the matrix to the reinforcing phase.

Conclusions

In this work, molecular dynamics simulations were used to investigate the interfacial interaction and reinforcement mechanism between nanocellulose and C-S-H. The results showed that the C-S-H model with a Ca/Si ratio of 1.67 constructed by the Monte Carlo method yielded a density of 2.44 g/cm³ and a Q_n distribution (Q₀ = 11.6%, Q₁ = 65.12%, Q₂ = 20%) that closely matches actual C-S-H gel. As the Ca/Si ratio increases from 0.75 to 1.67, the degree of polymerization of silicate chains decreases, and Young's modulus in all three directions decreases accordingly. Spontaneous adsorption occurs at the nanocellulose/C-S-H interface, where functional groups on the nanocellulose surface form hydrogen bonds to promote interfacial bonding. The number of hydrogen bonds formed at the interface of carboxylated nanocellulose/C-S-H (56.1) was significantly greater than that of hydroxylated nanocellulose/C-S-H (25.7). Both types of nanocellulose reinforce C-S-H, but the interfacial interaction energy of carboxylated nanocellulose/C-S-H (approximately -1.70×10² kJ/mol) is higher than that of hydroxylated nanocellulose, with electrostatic energy contributing over 90% of the total adsorption energy. These findings demonstrate that carboxylation modification significantly enhances interfacial bonding, providing an atomistic-scale theoretical basis for designing high-performance cement-based composites.

Key words

nanocellulose / molecular dynamics simulations / cement-based materials / calcium silicate hydrate gel

引用本文

林子韵, 刘巧玲, 孙晓宇, 张一晨, 任勇, 苏元鑫, 王桢. 纳米纤维素改性水泥基材料的分子动力学模拟. 硅酸盐学报, 2026 , 54 (10) : 3628 -3637 . DOI: 10.14062/j.issn.0454-5648.20250755

Ziyun LIN, Qiaoling LIU, Xiaoyu SUN, Yichen ZHANG, Yong REN, Yuanxin SU, Zhen WANG. Molecular Dynamics Simulation of Nanocellulose-Modified Cement-Based Materials[J]. Journal of the Chinese Ceramic Society, 2026 , 54 (10) : 3628 -3637 . DOI: 10.14062/j.issn.0454-5648.20250755

基金

收起

国家自然科学基金面上项目(52078282)

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

[1]

黄春龙, 王栋民, 白亚飞. 纳米微晶纤维素应用于水泥基材料研究综述[J]. 材料导报, 2018, 32(z1): 462-465.

HUANG Chunlong, WANG Dongmin, BAI Yafei. Mater Rep, 2018, 32(21): 462-465.

[2]

TRACHE D, TARCHOUN A F, ABDELAZIZ A, et al. A comprehensive review on processing, characteristics, and applications of cellulose nanofibrils/graphene hybrid-based nanocomposites: Toward a synergy between two-star nanomaterials[J]. Int J Biol Macromol, 2024, 268: 131633.

[3]

AL-ASKARY A S J, KOPECSKÓ K, OKTAY D. Cellulose nanocrystals and their utilization in cement-based composites: A comprehensive review[J]. Case Stud Constr Mater, 2025, 23: e05550.

[4]

HISSEINE O A, WILSON W, SORELLI L, et al. Nanocellulose for improved concrete performance: A macro-to-micro investigation for disclosing the effects of cellulose filaments on strength of cement systems[J]. Constr Build Mater, 2019, 206: 84-96.

[5]

ALREKABI S, CUNDY A B, LAMPROPOULOS A, et al. Mechanical performance of novel cement-based composites prepared with nano-fibres, and hybrid nano- and micro-fibres[J]. Compos Struct, 2017, 178: 145-156.

[6]

GÓMEZ HOYOS C, CRISTIA E, VÁZQUEZ A. Effect of cellulose microcrystalline particles on properties of cement based composites[J]. Mater Des, 2013, 51: 810-818.

[7]

SUN X X, WU Q L, LEE S, et al. Cellulose nanofibers as a modifier for rheology, curing and mechanical performance of oil well cement[J]. Sci Rep, 2016, 6: 31654.

[8]

AKHLAGHI M A, BAGHERPOUR R, KALHORI H. Application of bacterial nanocellulose fibers as reinforcement in cement composites[J]. Constr Build Mater, 2020, 241: 118061.

[9]

BHALERAO N, WAYAL A S, PATIL P G, et al. A review on effect of nano cellulose on concrete[J]. Int J Civ Struct Eng Res, 2015, 3(1): 251-254.

[10]

聂苗苗. 纳米纤维素改性水泥基材料力学性能研究[D]. 济南: 山东建筑大学, 2023.

NIE Miaomiao. Research on the mechanical properties of cement-based materials modified with nanofibrillated cellulose[D]. Jinan: Shandong Jianzhu University, 2023.

[11]

邓方茜, 唐可, 曹聪, 等. 纳米纤维素改性的再生混凝土骨料-基体界面微结构表征[J]. 水泥, 2024(11): 8-13.

DENG Fangqian, TANG Ke, CAO Cong, et al. Cement, 2024(11): 8-13.

[12]

黄卓琳, 田冬梅, 尹健, 等. 纳米纤维素对UHPC收缩性能的影响及机理[J]. 铁道科学与工程学报, 2025, 22(12): 5389-5397.

HUANG Zhuolin, TIAN Dongmei, YIN Jian, et al. J Railway Sci Eng, 2025, 22(12): 5389-5397.

[13]

CAO Y Z, ZAVATERRI P, YOUNGBLOOD J, et al. The influence of cellulose nanocrystal additions on the performance of cement paste[J]. Cem Concr Compos, 2015, 56: 73-83.

[14]

张静, 刘巧玲, 张姣, 等. 纳米纤维素改性水泥砂浆收缩变形性能研究[J]. 混凝土, 2025(3): 172-176.

ZHANG Jing, LIU Qiaoling, ZHANG Jiao, et al. Concrete, 2025(3): 172-176.

[15]

樊其昌, 王展鹏, 孟雪, 等. 纤维素纳米晶须水泥基复合材料试验及分子动力学模拟[J]. 复合材料学报, 2023, 40(9): 5315-5329.

FAN Qichang, WANG Zhanpeng, MENG Xue, et al. Acta Mater Compos Sin, 2023, 40(9): 5315-5329.

[16]

BASQUIROTO DE SOUZA F, SAGOE-CRENTSIL K, DUAN W H. A century of research on calcium silicate hydrate (C-S-H): Leaping from structural characterization to nanoengineering[J]. J Am Ceram Soc, 2022, 105(5): 3081-3099.

[17]

MASOUMI S, EBRAHIMI D, VALIPOUR H, et al. Nanolayered attributes of calcium-silicate-hydrate gels[J]. J Am Ceram Soc, 2020, 103(1): 541-557.

[18]

BAUCHY M, QOMI M J A, BICHARA C, et al. Rigidity transition in materials: Hardness is driven by weak atomic constraints[J]. Phys Rev Lett, 2015, 114(12): 125502.

[19]

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

[20]

CONG X D, JAMES KIRKPATRICK R. 17O and ²⁹Si MAS NMR study of β-C₂S hydration and the structure of calcium-silicate hydrates[J]. Cem Concr Res, 1993, 23(5): 1065-1077.

[21]

NISHIYAMA Y, LANGAN P, CHANZY H. Crystal structure and hydrogen-bonding system in cellulose Iβ from synchrotron X-ray and neutron fiber diffraction[J]. J Am Chem Soc, 2002, 124(31): 9074-9082.

[22]

IOELOVICH M. Preparation, characterization and application of amorphized cellulose: A review[J]. Polymers, 2021, 13(24): 4313.

[23]

MAZEAU K, HEUX L. Molecular dynamics simulations of bulk native crystalline and amorphous structures of cellulose[J]. J Phys Chem B, 2003, 107(10): 2394-2403.

[24]

CHEN W, LICKFIELD G C, YANG C Q. Molecular modeling of cellulose in amorphous state. Part I: Model building and plastic deformation study[J]. Polymer, 2004, 45(3): 1063-1071.

[25]

CHEN J J, THOMAS J J, TAYLOR H F W, et al. Solubility and structure of calcium silicate hydrate[J]. Cem Concr Res, 2004, 34(9): 1499-1519.

[26]

RAMASAMY J, AMANULLAH M. Nanocellulose for oil and gas field drilling and cementing applications[J]. J Petrol Sci Eng, 2020, 184: 106292.

[27]

SANCHEZ F, ZHANG L. Molecular dynamics modeling of the interface between surface functionalized graphitic structures and calcium-silicate-hydrate: Interaction energies, structure, and dynamics[J]. J Colloid Interface Sci, 2008, 323(2): 349-358.

[28]

ABDOLHOSSEINI QOMI M J, KRAKOWIAK K J, BAUCHY M, et al. Combinatorial molecular optimization of cement hydrates[J]. Nat Commun, 2014, 5: 4960.

[29]

HOU D S, ZHAO T J, MA H Y, et al. Reactive molecular simulation on water confined in the nanopores of the calcium silicate hydrate gel: Structure, reactivity, and mechanical properties[J]. J Phys Chem C, 2015, 119(3): 1346-1358.

[30]

HOU D S, ZHAO T J, JIN Z Q, et al. Molecular simulation of calcium silicate composites: Structure, dynamics, and mechanical properties[J]. J Am Ceram Soc, 2015, 98(3): 758-769.

2026年第54卷第10期

PDF下载

引用本文

BibTeX

文章信息

doi: 10.14062/j.issn.0454-5648.20250755

接收时间：2025-10-16
首发时间：2026-10-22
出版时间：2026-07-29

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2025-10-16
修回日期：2025-12-21

基金

国家自然科学基金面上项目(52078282)

作者信息

^1.山东建筑大学土木工程学院，济南　250101

^2.山东建筑大学建筑结构加固改进与地下空间工程教育部重点实验室，济南　250101

通讯作者:

刘巧玲（1977-），女，教授

王桢（1991-），女，讲师。

参考文献

分享链接

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

分享至

全文二维码

扫描看全文

引用本文

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