硅酸盐学报

机器学习优化Clausius-Mossotti介电常数计算方程

PDF下载

秦锦成 ¹ , 马名生 ¹ , 张发强 ¹ , 李永祥 ² , 刘志甫 ¹

硅酸盐学报 | 研究论文 2025,53(9): 2643-2650

收起

硅酸盐学报 |研究论文 2025 , 53 (9) : 2643 -2650

机器学习优化Clausius-Mossotti介电常数计算方程

全屏

秦锦成¹, 马名生¹, 张发强¹, 李永祥², 刘志甫¹

作者信息

^1．中国科学院上海硅酸盐研究所，关键陶瓷材料全国重点实验室，上海 201899

^2．皇家墨尔本理工大学工程学院，澳大利亚墨尔本 3000

通讯作者:

刘志甫（1975—），男，博士，研究员；

马名生（1986—），男，博士，研究员。

作者简介:

秦锦成（1995—），男，博士。

Machine Learning-Based Optimization of Clausius-Mossotti Equation for Dielectric Constant Calculation

Jincheng QIN¹, Mingsheng MA¹, Faqiang ZHANG¹, Yongxiang LI², Zhifu LIU¹

Affiliations

^1.State Key Laboratory of High Performance Ceramics and Superfine Microstructure, Shanghai Institute of Ceramics, Chinese Academy of Sciences, Shanghai 201899, China

^2.School of Engineering, RMIT University, Melbourne 3000, Australia

Corresponding Author:

LIU Zhifu (1975-), male, Ph.D., Professor. E-mail: liuzf@mail.sic.ac.cn;

MA Mingsheng (1986-), male, Ph.D., Professor. mamingsheng@mail.sic.ac.cn

QIN Jincheng (1995-), male, Ph.D. E-mail: qinjincheng@mail.sic.ac.cn

出版时间: 2025-09-20 doi: 10.14062/j.issn.0454-5648.20250039

文章导航

摘要

收起

介电常数的机器学习预测模型准确性虽然较经典Clausius-Mossotti方程更高，但模型解释性仍需提升。本工作提出了领域知识和机器学习双向嵌入的方法，在C-M方程中引入机器学习修正项以提升方程准确性。以机器学习模型所识别出的介电常数关键参量作为修正项的自变量，结合符号回归和遗传算法生成修正项的数学表达式。基于Pareto前沿选出了平衡精准度和复杂度的结果，并采用偏导数揭示了“低值补偿、高值抑制”的极化纠偏机制。对于采用Shannon极化率数据的介电常数计算，修正后的C-M方程R²提升了42.30%，RMSE下降了14.72%。本工作初步探索了“双向嵌入”方法在兼顾模型准确性和解释性方面的潜力，并尝试为材料设计提供新的思考角度。

关键词

介电常数 / Clausius-Mossotti方程 / 机器学习 / 符号回归 / 微波介质陶瓷

Abstract

收起

Introduction

The Clausius-Mossotti (C-M) equation is one of the most fundamental theoretical models to calculate the dielectric constant of materials. As was proposed over a century ago, it has been widely applied in dielectric material research due to its relatively simple mathematical form and clear physical interpretation. The equation, derived based on the Lorentz effective field approximation and the assumption of constant polarization, establishes a relationship among the macroscopic dielectric constant and the microscopic polarizability and unit cell volume. However, with the growing understanding of dielectric materials, it becomes clear that dielectric constant is affected by a range of additional factors, such as bond length, packing density and symmetry. These factors are not captured by the conventional C-M equation, indicating the limitations of the equation. Recently, machine learning (ML) methods show a promise in improving the prediction accuracy of dielectric properties. The ML models can outperform classical equations in terms of accuracy, but they often lack interpretability, which hinders their widespread use in material design and optimization. To address these issues, this study proposed a novel bidirectional embedding approach with domain knowledge and machine learning. The ML-based correction term was integrated into the classical C-M equation, enhancing its predictive accuracy, while retaining physical interpretability. This work could explore a potential of combining physical formulas with ML methods to create more robust and explainable models for material design.

Methods

The data of single-phase microwave dielectric ceramics were collected from published literatures and the Materials Project database. The dataset included materials whose dielectric constants were measured by the Hakki-Coleman method in a frequencies range from 5 to 18 GHz. A machine learning model identified molecular dielectric polarizability per volume (ppv), average bond length (blm) and unit cell volume (va) as the most significant contributors to the dielectric constant in a previous study. To enhance the accuracy of the C-M equation, correction terms were introduced based on the three key features. Symbolic regression, combined with genetic algorithms, was used to derive mathematical expressions for the correction terms. The symbolic regression technique utilized evolutionary algorithms to generate and evolve potential expressions. The structures and parameters of expressions were optimized to minimize prediction errors. Hyperparameters were optimized by a grid search method to identify the best-performing models. The process was carried out iteratively. A total of 16,000 candidate mathematical expressions were generated and compared. Finally, the Pareto front analysis was used to select the optimal expression that balances accuracy and complexity.

Results and Discussion

The revised C-M equation results in significant improvements in its predictive performance. For the Shannon polarizability dataset, the R² value of the revised C-M equation is increased by 42.29%, and the RMSE value is decreased by 14.72%. The correction term compensates for the inaccuracies of the original equation, particularly in cases where the molecular dielectric polarizability per volume of Shannon value (ppv_s) is either overestimated or underestimated. For low-polarizability materials with ppv_s values less than 0.2087, the Shannon database tends to underestimate the dielectric constant, and the correction term compensates for this error. Conversely, for high-polarizability materials, the dielectric constant is often overestimated, and the correction term reduces the overestimation. The analysis of the relationship between the correction term (Δppv_s) and features reveals that the correction term is most sensitive to ppv_s, followed by blm, with va having the least impact. This finding indicates that the dielectric polarizability plays a dominant role in determining the accuracy of the revised C-M equation. The influence of the features on the correction term is further explored through partial derivatives. The first-order partial derivatives show that the correction term's contribution from ppv_s is much greater than that from blm or va. The second-order partial derivatives reveal a non-linear relationship between the correction term and ppv_s and blm, while the relationship with va is linear. The results indicate that the compensatory effect of the correction term gradually diminishes, and the suppression effect becomes more pronounced as ppv_s increases. This behavior is consistent with the correction mechanism of "elevating the underrated, reducing the overrated".

Conclusions

This study introduced a bidirectional embedding approach that could integrate machine learning with domain knowledge to enhance both the accuracy and interpretability of the Clausius-Mossotti equation. The revised equation achieved a higher prediction accuracy with a 42.29% increase in R² and a 14.72% decrease in RMSE via incorporating the correction term derived from symbolic regression. The revised equation retained the physical meaning of the original one. The research could highlight a potential of bidirectional embedding approach to create models that could balance both accuracy and interpretability, offering a promising perspective on optimizing classical equations and material design.

Key words

dielectric constant / clausius-Mossotti equation / machine learning / symbolic regression / microwave dielectric ceramics

引用本文

秦锦成, 马名生, 张发强, 李永祥, 刘志甫. 机器学习优化Clausius-Mossotti介电常数计算方程. 硅酸盐学报, 2025 , 53 (9) : 2643 -2650 . DOI: 10.14062/j.issn.0454-5648.20250039

Jincheng QIN, Mingsheng MA, Faqiang ZHANG, Yongxiang LI, Zhifu LIU. Machine Learning-Based Optimization of Clausius-Mossotti Equation for Dielectric Constant Calculation[J]. Journal of the Chinese Ceramic Society, 2025 , 53 (9) : 2643 -2650 . DOI: 10.14062/j.issn.0454-5648.20250039

基金

收起

国家自然科学基金联合基金(U24A2052)
国家资助博士后研究人员计划C档(GZC20232825)
上海市东方英才计划(QNKJ2024026)

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

2025年第53卷第9期

PDF下载

引用本文

BibTeX

文章信息

doi: 10.14062/j.issn.0454-5648.20250039

接收时间：2025-01-15
首发时间：2026-10-21
出版时间：2025-09-20

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2025-01-15
修回日期：2025-02-12

基金

国家自然科学基金联合基金(U24A2052)

国家资助博士后研究人员计划C档(GZC20232825)

上海市东方英才计划(QNKJ2024026)

作者信息

^1．中国科学院上海硅酸盐研究所，关键陶瓷材料全国重点实验室，上海 201899

^2．皇家墨尔本理工大学工程学院，澳大利亚墨尔本 3000

通讯作者:

刘志甫（1975—），男，博士，研究员；

马名生（1986—），男，博士，研究员。

参考文献

分享链接

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

分享至

全文二维码

扫描看全文

引用本文

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