蛭石晶胞同构置换后，用forcite模块clayff力场分配力场和电荷，不呈电中性

MatSrookie

原文描述：

                              

文章中化学式M1.3125（Si5.875Al2.125）(Al0.8125Mg5.1875)O20(OH)4，根据后文M为补偿阳离子，Al3+替换Si4+,Al3+替换Mg2+,推测最初的化学式为Mg3 Si4 O10 (O H)2。无法找到Tambach等人的蛭石结构cif文件。用下图中cif文件代替：

将其导入MS软件后，点自动成键，将Mg-O八面体手动相连、加氢，建立4*2超胞

[ cell form: H64 O384 Mg96 Si128 ]，通过SubstitutionalDisorder脚本用13个Al替换Mg（13.54%，+13），34个Al替换Si（26.5625%，-34）, 最终得到cellform为Si94Al47 Mg83 O384 H64的结构，整理为[(Si94 Al34)(Al13Mg83)O320(OH)64]，即(Si5.8752 Al2.125)(Al0.8125Mg5.1875)O20(OH)4，为了平衡同构置换取代引起的负电荷，加入21个K+，用clayff力场分配力场和电荷，发现整体电荷不为0。该怎么解决？

张同学快做实验

楼主方便分享一下建好的蛭石XSD文件分享一下吗，我最近也想做蛭石的相关计算，但是Mg-O一直不成键，谢谢~

yaneeee

请问您解决了吗

slxc920113

MS里面的力场电荷是写死的，本来就不能保证电荷中性，AuToFF（https://cloud.hzwtech.com/web/pe ... /inorganic-crystals）将ClayFF力场的电荷改成了BCI计算规则，在固定阳离子电荷的前提下，可以自动调整O的电荷来保证体系的电中性，前提是用户自己上传的结构是正确的，比如Al或者Mg置换Si之后需要补充阳离子来保证电荷平衡。