计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: madhatter
打印 Print 上一主题 Last thread 下一主题 Next thread

[综合交流] 计算锌离子和葡萄糖结合能初始模型设置问题

[复制链接 Copy URL]

63

帖子

0

威望

723

eV
积分
786

Level 4 (黑子)

16#
 楼主 Author| 发表于 Post on 2022-1-20 22:55:02 | 只看该作者 Only view this author
sobereva 发表于 2022-1-18 12:49
AIMD本来就属于MD的一种,当前语境你应当把MD称为基于经典力场的MD

1 如果是基于DFT做的AIMD,只要计 ...

在了解aimd,使用cp2k中用的xtb是gnf1,
1.这个精度对于模拟锌离子和中性分子和氯离子的动力学是否可以。
2.常规的aimd使用的是什么方法和基组,简单搜了几篇电解液文章只写dftmd的字样。
3.我做电解液的aimd,目的是看锌离子配位情况,800原子用xtb还是再小一些的体系用b3lyp之类的。打扰sob老师了

6万

帖子

99

威望

5万

eV
积分
120160

管理员

公社社长

17#
发表于 Post on 2022-1-20 23:18:52 | 只看该作者 Only view this author
madhatter 发表于 2022-1-20 22:55
在了解aimd,使用cp2k中用的xtb是gnf1,
1.这个精度对于模拟锌离子和中性分子和氯离子的动力学是否可以 ...

那叫GFN1-xTB

1 大概率可以,但貌似没人测过,不和DFT对比谁也不好打包票
2 文章里不可能没写明,正文里没写的话SI里也会写,SI里真若没写就发邮件问作者。具体用什么合适看用什么程序
3 xtb不怎么适合周期性。CP2K里的GFN1-xTB可以用跑。杂化泛函在第一性原理程序里极度昂贵,甭指望杂化泛函在CP2K里跑得动AIMD
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

63

帖子

0

威望

723

eV
积分
786

Level 4 (黑子)

18#
 楼主 Author| 发表于 Post on 2022-1-21 21:45:46 | 只看该作者 Only view this author
本帖最后由 madhatter 于 2022-1-21 22:03 编辑
sobereva 发表于 2022-1-20 23:18
那叫GFN1-xTB

1 大概率可以,但貌似没人测过,不和DFT对比谁也不好打包票

我用超算跑800原子的aimd用的gfn1-xtb,1fs跑了10000步,300k,输入文件建立方式是模仿sob老师的216个水的那个例子改的,用vmd处理溶液盒子周期性后,观察到,一共放了五个锌离子,基本倾向与氯离子配位,其中有两个锌离子都形成了zncl4的结构,其他也基本是zncl2eg的结构,现象基本符合预期,不过sob老师说xtb不太适合周期性,想问一下需要考虑变换dft方法算吗,这算是第一次跑出理想的结构有点开心?。谢谢sob老师!

6万

帖子

99

威望

5万

eV
积分
120160

管理员

公社社长

19#
发表于 Post on 2022-1-21 23:18:49 | 只看该作者 Only view this author
madhatter 发表于 2022-1-21 21:45
我用超算跑800原子的aimd用的gfn1-xtb,1fs跑了10000步,300k,输入文件建立方式是模仿sob老师的216个水 ...

分清楚xtb和GFN-xTB,完全不是同一个层面上的东西
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

63

帖子

0

威望

723

eV
积分
786

Level 4 (黑子)

20#
 楼主 Author| 发表于 Post on 2022-1-24 23:04:11 | 只看该作者 Only view this author
sobereva 发表于 2022-1-21 23:18
分清楚xtb和GFN-xTB,完全不是同一个层面上的东西

sob老师晚上好,我用cp2k 里gfn1-  xtb,跑了一个包括锌离子 乙二醇  胆碱阳离子  氯离子的800原子体系,300k,跑了15ps,用rdf可以得到zn周围cl的配位为2。1.我从能量和温度来看在10ps就趋于稳定了,结构初步复合我预期,但是纠结是否还有继续跑的必要,继续跑结构会不会更理想,还是换pbe方法再跑。

2.虽然文献里实验得到的体系锌离子团簇唯一是zncl4 2-,我体系模拟的rdf全出的配位是zncl2,这个方法算出的配位,与真实配位是否一样。

3.我想算氯离子和尿素的氢键作用   氯离子和乙二醇的氢键作用做比较,一般有什么做法,结合能计算?还是用社长的软件可视化氢键互相比较?
4.做锌电极带电表面的电解液动力学,使用cp2k是不是比较方便,因为不太涉及反应,只涉及到锌金属表面的分子排布,当然有氯离子和中性分子,这种也是aimd去做吗,qm/mm?不太了解,打扰sob老师了。

6万

帖子

99

威望

5万

eV
积分
120160

管理员

公社社长

21#
发表于 Post on 2022-1-26 01:45:56 | 只看该作者 Only view this author
madhatter 发表于 2022-1-24 23:04
sob老师晚上好,我用cp2k 里gfn1-  xtb,跑了一个包括锌离子 乙二醇  胆碱阳离子  氯离子的800原子体系, ...

1 肯定得换DFT继续跑,否则缺乏说服力。就算真的只用GFN-xTB来跑,也得从优化的单体结构、配体解离能等角度证明GFN-xTB和DFT(或其它靠谱的方法)算出来的结果相近。

2 DFT跑完了再说

3 看你想比较什么。比较结构、结合能,就把Cl放氢对面然后优化结构、计算结合能之类的
可视化作用区域可以用Multiwfn做IRI或IGMH分析。

4 CP2K跑DFT做AIMD研究这个没问题。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

63

帖子

0

威望

723

eV
积分
786

Level 4 (黑子)

22#
 楼主 Author| 发表于 Post on 2022-2-9 19:01:33 | 只看该作者 Only view this author
本帖最后由 madhatter 于 2022-2-9 21:41 编辑
sobereva 发表于 2022-1-26 01:45
1 肯定得换DFT继续跑,否则缺乏说服力。就算真的只用GFN-xTB来跑,也得从优化的单体结构、配体解离能等角 ...

sob老师过年好,我用dft跑完了两种体系,结果很符合预期,现在打算做锌离子与尿素中氧的径向分布函数,因为锌离子放得少,所以打算做稳定阶段最后5ps的轨迹里锌离子与尿素中氧的径向分布函数。vmd中单纯一帧里用name O and index xxx to xxx可以得到不同结构里的氧和锌的rdf,但是对于5000帧我就不知道怎么去处理了。我将两种结构做了静电势分析,zncl4明显更负,我打算说zncl4得电子能力较弱,锌离子不容易被还原。不知道合不合理。
其次是sob老师有推荐关于固液界面建模的方案吗,再搜这部分内容没有太详细的教程。打算做锌电极表面做电解液的aimd。

图片1.png (5.54 KB, 下载次数 Times of downloads: 48)

图片1.png

图片2.png (5.43 KB, 下载次数 Times of downloads: 64)

图片2.png

6万

帖子

99

威望

5万

eV
积分
120160

管理员

公社社长

23#
发表于 Post on 2022-2-10 10:11:33 | 只看该作者 Only view this author
madhatter 发表于 2022-2-9 19:01
sob老师过年好,我用dft跑完了两种体系,结果很符合预期,现在打算做锌离子与尿素中氧的径向分布函数,因 ...

Extensions - Analysis - Radial pair distribution function直接就能用来对轨迹计算rdf

没什么问题

建模没什么费劲的。金属界面弄好了,用packmol之类插入液体分子到它上面,然后再做控压MD消掉真空区就完了。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

63

帖子

0

威望

723

eV
积分
786

Level 4 (黑子)

24#
 楼主 Author| 发表于 Post on 2022-2-12 23:23:03 | 只看该作者 Only view this author
本帖最后由 madhatter 于 2022-2-12 23:27 编辑
sobereva 发表于 2022-2-10 10:11
Extensions - Analysis - Radial pair distribution function直接就能用来对轨迹计算rdf

没什么问题

谢谢sob老师,1.vmd做rdf的原子选择如何弄, 我用的type zn       name O and index xxx to xxx(xyz格式取最后一帧的尿素编号),这样算下来应该是一帧里锌和尿素里氧的rdf还是统计了所有帧的锌与尿素中的氧,我想得到的是5000帧里锌与尿素中的氧的rdf平均(体系中还有别的分子含有氧)。
2.其次目前建模了固液界面,建立cp2k的输入文件,显示的这个模型是xyz周期性吗?不太理解这种固液界面建模原理。我用ms的bulid layer建的模,先固定了固体金属,然后优化(对于离子液体这种溶剂,需要用溶剂连续那个选项吗?),优化后加了电场去跑md,结果是基本不收敛,每一步都是达到200上限,想问一下如何解决。

3.最近在看二代cpmd,想问一下关键词输入是Multiwfn中控温选项里的LANGEVIN吗?还是有别的关键词使用二代的cpmd,这种方法可以用在过渡金属上吗,例如锌晶面的电解液动力学。 打扰sob老师了
&MD

ENSEMBLE LANGEVIN

&LANGEVIN

GAMMA 0.001 # gamma L

NOISY_GAMMA 1E-4 #gamma D

&END LANGEVIN

6万

帖子

99

威望

5万

eV
积分
120160

管理员

公社社长

25#
发表于 Post on 2022-2-13 03:17:04 | 只看该作者 Only view this author
madhatter 发表于 2022-2-12 23:23
谢谢sob老师,1.vmd做rdf的原子选择如何弄, 我用的type zn       name O and index xxx to xxx(xyz格式 ...

1 VMD自带的rdf插件里直接就能设考虑的帧号范围

2 语义不明。是否当周期性考虑看你是否恰当设置了盒子,以及输入文件里怎么设的关键词
把问题充分描述清楚,“溶剂连续”是什么鬼,如果是连续溶剂模型(隐式溶剂模型),干嘛要开它?加电场要做甚?

3 没事甭鼓捣非主流的东西,比起直接跑BOMD又没什么非常显著的好处
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

63

帖子

0

威望

723

eV
积分
786

Level 4 (黑子)

26#
 楼主 Author| 发表于 Post on 2022-2-25 21:59:16 | 只看该作者 Only view this author
本帖最后由 madhatter 于 2022-2-25 22:08 编辑
sobereva 发表于 2022-2-13 03:17
1 VMD自带的rdf插件里直接就能设考虑的帧号范围

2 语义不明。是否当周期性考虑看你是否恰当设置了盒子 ...

sob老师打扰了,我做了界面aimd,首先建模了四层锌金属原子和电解液的界面,我把四层锌原子都固定,对电解液进行了热裕,真空层10A,根据Multiwfn做了输入文件,电场强度用了0.0002a.u.我先进行了结构优化,然后在进行的md(无法加K点,一加K点就报错说这种方法不能用K点)
1.我想做电池负极界面的模拟,这个电场强度方向怎么分辨,我用正的方向可以进行md,如果用负的则很难收敛,建模是否有误(我另一种个想法是做两个锌金属层夹着电解液 从经典md文章找到类似的)
2.使用OT选型收敛速度更快些,OT和另一个应该怎么选
3.目前用64核节点跑该体系290原子,12h跑了500步,想问一下有别的方法优化这个速度吗?aimd需要跑10ps以上吗?(我想讨论的仅仅是电极表面受电场作用下分子分布和移动,例如氯离子的吸附等问题,不计划引入反应,因此电场给的不高)
4.在思考如果讨论电场作用下的界面分子排布是否更大体系的经典md更适合,再看gmx的教程,关于拓扑文件一直不太懂,如果的建模是多种类分子原子,是不是我需要生成每种分子原子的拓扑文件然后放在一起去合成拓扑文件,还是有类似软件可以读取模型结合力场生成拓扑文件。
最后推荐一篇文章,最近浙江工业的science,https://www.science.org/doi/10.1126/science.abn1818用的vasp做的界面电解液aimd,使用促进阴离子分解用于形成sei膜,实现长循环电池。他的研究是在0-450fs内完成了阴离子分解再与锂结合的反应,一共跑了12ps。很好奇,一直到12ps是为了什么?

微信图片_20220225215259.png (108.74 KB, 下载次数 Times of downloads: 16)

微信图片_20220225215259.png

微信图片_20220225215411.png (186.73 KB, 下载次数 Times of downloads: 19)

微信图片_20220225215411.png

science1.png (121.4 KB, 下载次数 Times of downloads: 20)

文献电解液分解

文献电解液分解

science2.png (203.01 KB, 下载次数 Times of downloads: 19)

文献模拟全程

文献模拟全程

layer2.inp

22.8 KB, 下载次数 Times of downloads: 2

6万

帖子

99

威望

5万

eV
积分
120160

管理员

公社社长

27#
发表于 Post on 2022-2-25 22:05:12 | 只看该作者 Only view this author
madhatter 发表于 2022-2-25 21:59
sob老师打扰了,我做了界面aimd,首先建模了四层锌金属原子和电解液的界面,我把四层锌原子都固定,对电 ...

2 “另一个”语义不明

3 Z方向没必要考虑PBC,这样真空层还能减小不少来节约时间。其它地方没什么油水,虽然用GFN-xTB便宜得多但精度和合理性上有风险

4 看有无类似问题的研究文章,如果有用经典力场跑的,也可以用GROMACS效仿着跑。各种离子和分子的拓扑信息都得有才能跑。用单原子离子有现成的,或者基于文献里的参数可以轻易手写,分子或多原子离子的可以用Sobtop(http://sobereva.com/soft/Sobtop)产生拓扑文件。
北京科音分子动力学与GROMACS培训班对动力学模拟的原理和GROMACS有十分系统的讲解,对各类体系给了海量例子,也有金属的模拟部分:http://www.keinsci.com/workshop/KGMX_content.html,学过后自然就能搞懂怎么实现模拟
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

63

帖子

0

威望

723

eV
积分
786

Level 4 (黑子)

28#
 楼主 Author| 发表于 Post on 2022-2-25 22:14:49 | 只看该作者 Only view this author
sobereva 发表于 2022-2-25 22:05
2 “另一个”语义不明

3 Z方向没必要考虑PBC,这样真空层还能减小不少来节约时间。其它地方没什么油水 ...

谢谢sob老师,最近感觉aimd还是很有用的,就是比较吃资源,我打算经典和aimd都试一下。

63

帖子

0

威望

723

eV
积分
786

Level 4 (黑子)

29#
 楼主 Author| 发表于 Post on 2022-2-28 22:31:03 | 只看该作者 Only view this author
本帖最后由 madhatter 于 2022-2-28 22:37 编辑
sobereva 发表于 2022-2-25 22:05
2 “另一个”语义不明

3 Z方向没必要考虑PBC,这样真空层还能减小不少来节约时间。其它地方没什么油水 ...

想问一下结合能计算的问题
目前两个结构,一个是[zncl3eg]-是a结构,一个是zncl2egurea是二结构,eg是乙二醇,urea是尿素
1我想说明尿素与锌结合能更高,是E总-Ezn-2Ecl-Eurea-Eeg和E总-Ezn-E3cl-Eeg相比较?还是a结构解离一个cl和b结构解离一个尿素的能量相比,有点懵。
2我通过aimd算出来的rdf显示尿素要比氯离锌离子更近,是不是可以说尿素的结合能更高
3.最近研究想验证一个想法,多种原子的氢键作用是否可以固定氯原子,例如经常电化学反应过程中,氯离子在金属表面吸附,我想讨论讨论不同组分与氯离子产生多种氢键,使它固定在溶液‘’笼子‘’中,吸附在金属表面数量变少。初步想法是用orca跑电场下的aimd,不同电场大小下,氯离子氢键的断裂,去比较断裂发生的电场大小?或者计算不同构形和氯离子的氢键作用大小?不知道有别的类似的研究吗?
4 因为电化学过程出现了双电层超吸附现象,电化学还原过程出现成核过电位,即抑制了还原,所以我通过对基体和基体加吸附体系中各种分子氯离子等的计算,去比较他们之间的禁带宽度,讨论导电性,如果导电性变差是否可以证明这种抑制作用,还是有其他的方法去讨论。
打扰sob老师了

63

帖子

0

威望

723

eV
积分
786

Level 4 (黑子)

30#
 楼主 Author| 发表于 Post on 2022-3-3 23:12:53 | 只看该作者 Only view this author
sobereva 发表于 2022-2-25 22:05
2 “另一个”语义不明

3 Z方向没必要考虑PBC,这样真空层还能减小不少来节约时间。其它地方没什么油水 ...

1.最近再跑金属界面电解液的aimd,先优化了结构,然后在电场作用下跑,我想问的是,aimd从开始跑例如1-1000fs阶段出现的变化,可以认为是可能发生的现象吗?例如我的体系500步以为氯离子去和尿素分子形成氢键,并稳定结构到5000fs,在5000fs后氯离子脱离然后吸附在锌金属表面,我想讨论这个形成氢键的过程,我认为是电化学反应中可能存在的现象,可以这么讨论吗?经典md似乎需要讨论长时间md平衡后的结构,而aimd从md开始在电场下就在不断运动,我理解后者的每一步似乎都是可能存在的。不知道这么理解是否正确。2.我跑了溶液盒子体系的aimd,跑了20ps,打算用msd算锌离子的扩散系数,我看其他文献横坐标都是ns或几千ps级别的,我的横坐标只有8ps,是否aimd可以做短时间的msd,如果要加时间,我看有关文献做的40ps,是否可行。

mmexport1646320340891.jpg (300.61 KB, 下载次数 Times of downloads: 13)

mmexport1646320340891.jpg

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-16 16:57 , Processed in 0.252362 second(s), 22 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list