应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 用DL_POLY模拟了硅酸盐熔体后，如何计算硅氧四面体Qn结构？
91/1返回列表
查看: 1814  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

killbast

铜虫 (初入文坛)

[求助] 用DL_POLY模拟了硅酸盐熔体后，如何计算硅氧四面体Qn结构？

模拟了CaO,Al2O3，SiO2的体系，然后想计算其中硅氧四面体类型分布，即计算与所有硅配位的氧原子是否桥氧，
回复此楼

» 收录本帖的淘帖专辑推荐

动力学模拟

» 猜你喜欢
1楼 2014-04-28 15:33:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

我是砖头

银虫 (小有名气)
楼主您好；
我最近也在看硅酸盐熔体模拟文献，我能问您一个问题吗，就是我们模拟之前的初始构象怎么建？新人，求助
赞一下 回复此楼
2楼2015-09-14 15:12:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

15216105346

金虫 (著名写手)
这是统计程序，matlab的！
供你参考一下！

clear
Ob=0;Onb=0;Of=0;
A=xlsread('CS1');
[m,n]=size(A);
Si=0;Ca=0;O=0;
for i=2:m
    if A(i,4)==1
        Si=Si+1;
    elseif A(i,4)==2
        Ca=Ca+1;
    else
        O=O+1;
    end
end

for i=Si+Ca+2:Si+Ca+O+1 %O的个数
    for j=2:Si+Ca+1 %Si和Ca的个数
        B(i-Si-Ca-1,j-1)=sqrt((A(i,1)-A(j,1))^2+(A(i,2)-A(j,2))^2+(A(i,3)-A(j,3))^2);
    end
    [P,Q]=sort(B(i-Si-Ca-1,1:Si+Ca));
    if Q(1)<Si+2&&Q(2)<Si+2
        Ob=Ob+1;
    end
    if Q(1)<Si+2&&Q(2)>Si+1
        Onb=Onb+1;
    end
    if Q(1)>Si+1&&Q(2)<Si+2
        Onb=Onb+1;
    end
    if Q(1)>Si+1&&Q(2)>Si+1
        Of=Of+1;
    end
end

Ob=Ob/O*100
Onb=Onb/O*100
Of=Of/O*100

发自小木虫IOS客户端
赞一下 回复此楼
3楼2015-09-14 23:37:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

我是砖头

银虫 (小有名气)
引用回帖:
3楼: Originally posted by 15216105346 at 2015-09-14 23:37:39
这是统计程序，matlab的！
供你参考一下！
clear
Ob=0;Onb=0;Of=0;
A=xlsread('CS1');
=size(A);
Si=0;Ca=0;O=0;
for i=2:m
    if A(i,4)==1
        Si=Si+1;
    elseif A(i,4)==2
        Ca=Ca+1;
    else ...

你好，我想问一下，用matlab可以统计出桥氧与非桥氧的百分数吗？
赞一下 回复此楼
4楼2016-04-03 13:23:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

15216105346

金虫 (著名写手)
引用回帖:
4楼: Originally posted by 我是砖头 at 2016-04-03 13:23:13
你好，我想问一下，用matlab可以统计出桥氧与非桥氧的百分数吗？...

MATLAB是可以做到的！我做过这个工作的

发自小木虫IOS客户端
赞一下 回复此楼
5楼2016-04-03 13:47:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

我是砖头

银虫 (小有名气)
引用回帖:
5楼: Originally posted by 15216105346 at 2016-04-03 13:47:49
MATLAB是可以做到的！我做过这个工作的
...

我看了下这个编程，还有个问题想请教你，这个读入的CS1文件数据是什么数据啊，还有这个m、n、P、Q分别代表的是什么，谢谢
赞一下 回复此楼
6楼2016-04-07 14:57:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

15216105346

金虫 (著名写手)
引用回帖:
6楼: Originally posted by 我是砖头 at 2016-04-07 14:57:13
我看了下这个编程，还有个问题想请教你，这个读入的CS1文件数据是什么数据啊，还有这个m、n、P、Q分别代表的是什么，谢谢...

m n是矩阵的行列数，PQ是矩阵！看看matlab 教材吧

发自小木虫IOS客户端
赞一下 回复此楼
7楼2016-04-07 17:59:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

回形针的生活

新虫 (初入文坛)
楼主您好，我也想统计下硅酸盐的Qn结构，不知您现在是否找到方法，希望能得到您的指教
赞一下 回复此楼
8楼2016-04-27 21:37:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

回形针的生活

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by 15216105346 at 2015-09-14 23:37:39
这是统计程序，matlab的！
供你参考一下！
clear
Ob=0;Onb=0;Of=0;
A=xlsread('CS1');
=size(A);
Si=0;Ca=0;O=0;
for i=2:m
    if A(i,4)==1
        Si=Si+1;
    elseif A(i,4)==2
        Ca=Ca+1;
    else ...

您好，我现在想统计出Si-O-Si、Si-O-Al，Al-O-Al桥氧各自所占的比例，不知道您有没有这方面的程序，谢谢
赞一下 回复此楼
9楼2016-04-27 21:45:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 killbast 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿华中农业071010，总分320求调剂 +5 困困困困坤坤 2026-03-20 6/300 2026-03-22 17:41 by hxsm
[考研] 寻找调剂 +4 倔强芒? 2026-03-21 4/200 2026-03-22 16:14 by 木托莫露露
[考研] 319求调剂 +4 小力气珂珂 2026-03-20 4/200 2026-03-22 15:53 by ColorlessPI
[考博] 招收博士1-2人 +3 QGZDSYS 2026-03-18 4/200 2026-03-22 10:25 by QGZDSYS
[基金申请] 山东省面上项目限额评审 +4 石瑞0426 2026-03-19 4/200 2026-03-22 08:50 by Wei_ren
[考研] 085600材料与化工306 +4 z1z2z3879 2026-03-21 4/200 2026-03-21 23:44 by ms629
[考研] 278求调剂 +9 烟火先于春 2026-03-17 9/450 2026-03-21 17:47 by 学员8dgXkO
[考研] 307求调剂 +3 余意卿 2026-03-18 3/150 2026-03-21 17:31 by ColorlessPI
[考研] 22 350 本科985求调剂，求老登收留 +3 李轶男003 2026-03-20 3/150 2026-03-21 13:28 by 搏击518
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂！ +3 不想理你 2026-03-16 3/150 2026-03-21 03:59 by JourneyLucky
[考研] 初始318分求调剂（有工作经验） +3 1911236844 2026-03-17 3/150 2026-03-21 02:33 by JourneyLucky
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 一志愿重庆大学085700资源与环境专硕，总分308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-21 00:39 by JourneyLucky
[考研] 311求调剂 +5 冬十三 2026-03-18 5/250 2026-03-21 00:16 by JourneyLucky
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 材料与化工专硕调剂 +7 heming3743 2026-03-16 7/350 2026-03-20 19:31 by zhukairuo
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考研] 【同济软件】软件（085405）考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 0854可跨调剂，一作一项核心论文五项专利，省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 11408 一志愿西电，277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
信息提示
关闭
请填处理意见
关闭