24小时热门版块排行榜    

查看: 934  |  回复: 1

15216105346

金虫 (著名写手)

[求助] 怎样求总径向分布函数 已有1人参与

各位大侠!
液体的总径向分布函数和各离子对的偏径向分布函数是什么关系!
怎样求出总径向分布函数?
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abinitio

版主 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
其实木虫上有许多代码,搜索一下即可

比如:by @chuntao118
原贴
关于径向分布函数RDFs的代码(求指正)! - 分子模拟 - 小木虫论坛 - 学术科研第一站
http://muchong.com/bbs/viewthread.php?tid=8813127&fpage=1&target=blank
CODE:
我把写的C代码(源自于《the art of MD)附上,自己参考一下,
void EvalRdf ()
{
        VecR dr;
        real deltaR, normFac, rr;
        int j1, j2, n;
        
        if (countRdf == 0)
        {
                for (n = 0; n < sizeHistRdf; n ++)
                {
                        histRdf[n] = 0.;
                }
        }
        deltaR = rangeRdf / sizeHistRdf;

        for (j1 = 0; j1 < nMol - 1; j1 ++)
        {
                for (j2 = j1 + 1; j2 < nMol; j2 ++)
                {
                        VSub (dr, mol[j1].r, mol[j2].r);
                        VWrapAll (dr);
                        rr = VLenSq (dr);
                        if (rr < Sqr (rangeRdf))
                        {
                                n = sqrt (rr) / deltaR;
                                ++ histRdf[n];
                        }
                }
        }
        ++ countRdf;
        if (countRdf == limitRdf)
        {
                normFac = VProd (region) / (2. * M_PI * Cube (deltaR) * Sqr (nMol) * countRdf);
                for (n = 0; n < sizeHistRdf; n ++)
                {
                        histRdf[n] *= normFac / Sqr (n - 0.5);
                }
                PrintRdf (stdout);
                countRdf = 0;
        }
}
fortan的一些结构写法我不太明白。你自己对照改一下就行了,算法都一样。

2楼2015-10-29 08:34:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 15216105346 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0854电子信息求调剂 +3 α____ 2026-03-22 3/150 2026-03-22 21:28 by zhq0425
[考研] 328求调剂,英语六级551,有科研经历 +6 生物工程调剂 2026-03-17 10/500 2026-03-22 20:22 by edmund7
[考研] 287求调剂 +8 晨昏线与星海 2026-03-19 9/450 2026-03-22 17:01 by i_cooler
[考研] 308求调剂 +3 墨墨漠 2026-03-21 3/150 2026-03-22 16:54 by i_cooler
[考研] 311求调剂 +3 26研0 2026-03-20 3/150 2026-03-22 14:46 by ColorlessPI
[考研] 求调剂 +5 Zhangbod 2026-03-21 7/350 2026-03-22 13:13 by Zhangbod
[考研] 一志愿华中科技大学071000,求调剂 +4 沿岸有贝壳6 2026-03-21 4/200 2026-03-22 07:21 by ilovexiaobin
[考研] 286分人工智能专业请求调剂愿意跨考! +4 lemonzzn 2026-03-17 8/400 2026-03-21 22:49 by lemonzzn
[考研] 材料学硕301分求调剂 +7 Liyouyumairs 2026-03-21 7/350 2026-03-21 22:31 by peike
[考研] 求助 +5 梦里的无言 2026-03-21 6/300 2026-03-21 17:51 by 学员8dgXkO
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] 330求调剂0854 +3 assdll 2026-03-21 3/150 2026-03-21 13:01 by 搏击518
[考研] 332求调剂 +3 凤凰院丁真 2026-03-20 3/150 2026-03-21 10:27 by luoyongfeng
[考研] 一志愿西南交大,求调剂 +5 材化逐梦人 2026-03-18 5/250 2026-03-21 00:26 by JourneyLucky
[考研] 一志愿西安交通大学 学硕 354求调剂211或者双一流 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
[考研] 生物学调剂招人!!! +3 山海天岚 2026-03-17 4/200 2026-03-19 21:34 by 怎么释怀
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[论文投稿] 有没有大佬发小论文能带我个二作 +3 增锐漏人 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
信息提示
请填处理意见