分享材料方向的滑稽计算

AxiEJohn

Uus/pMeC6H4-/キ 发表于 2025-12-27 14:16
根据https://doi.org/10.1002/slct.202500161这篇理论计算文章的摘要，一个非整数计量比的混合氧化物钙钛矿 ...

ChemistrySelect - 2025 - Likhariya - Density Functional Theory DFT Study of La.pdf (3.76 MB, 下载次数 Times of downloads: 13)

2025-12-27 15:43 上传 Uploaded

点击下载
Click to download

LPOwenChan

Uus/pMeC6H4-/キ 发表于 2025-12-27 14:16
根据https://doi.org/10.1002/slct.202500161这篇理论计算文章的摘要，一个非整数计量比的混合氧化物钙钛矿 ...

粗略地看了下，承包了我今年的笑点：
    1. 把钙钛矿氧化物LaMn0.4Fe0.6O3作为孤立体系来计算，实际计算中是按照Table4中的LaMn3Fe5O12来计算的，不明觉厉；
    2. 用了B3LYP和PBE泛函来计算，并且还强调是Restricted方法，计算结果对于过渡金属离子来说估计是毫无意义的；
    3. Table5中直接将HOMO-LUMO Gap当作是Band Gap，算出来了-0.012eV这种金属体系才会有的数值，但是引用的实验文献却说是“宽禁带半导体”；
    4. Figure2和Table5中的La的Mulliken电荷为+3.6，La会有这种”价态“吗？
    5. 计算了这个”分子“的偶极矩与由气相频率导出的热容和熵，这真的有任何意义吗？
    6. 煞有介事地给出了这个计算级别下的IR，作者真的知道自己在干什么吗...?
其他一些相当”有意思“的细节：
    1. 经典GaussView直出；

    2. 似乎也是GaussView直出；

    3. 绝对是Excel的柱状图绘图直出，连色彩序列都是默认的。

其他估计还有很多笑点，但是实在不想浪费时间去细看了...
结论1: 写文章的人和审文章的人至少有一方不是人类；
结论2: 原来我不太相信AI能取代一部分人的工作，现在我相信了；
结论3: 仔细看完这篇文章能不笑的人应该可以确诊抑郁症了。
以上。
ChemistrySelect - 2025 - Likhariya - Density Functional Theory DFT Study of La.pdf (3.76 MB, 下载次数 Times of downloads: 19)

2025-12-27 17:05 上传 Uploaded

点击下载
Click to download

不想飞的猫头鹰

https://doi.org/10.1002/adma.202301398 作者声称用GROMACS进行了COMPASS力场的计算
https://doi.org/10.1002/adma.202512960 作者声称用VASP进行了B3LYP def2-TZVP的计算

lindaiyu

真的难崩啊，发的都还是化学优质期刊呢

Observer-Aya

Uus/pMeC6H4-/キ 发表于 2025-11-28 11:39
逛咱这论坛灌水区的人终归有限，建议像http://bbs.keinsci.com/thread-52582-1-1.html提到的那样，写好英文 ...

老师请问做上千K温度的AIMD有什么问题吗，我一度打算用CP2K做个一千多个原子的973K的AIMD来着，后来因计算量太大放弃了

SharkYYX2025

Observer-Aya 发表于 2026-3-14 14:23
老师请问做上千K温度的AIMD有什么问题吗，我一度打算用CP2K做个一千多个原子的973K的AIMD来着，后来因计 ...

好吓人的计算量，拿GFN1-xTB没准能行，DFT做不动吧
高温AIMD不是说一定不行，我个人理解是看物理意义。需要有理由用高温，不能无端设定。这个文章做的东西我也不懂，不好评价。如果他的物理模型不是设定在高温，或者对应实际应用一般是常温，那这个说不通。
比如我在我的研究中，为了研究两个有机小分子的热稳定性，拿ORCA做了高温AIMD，看共价键裂解。这种通常符合物理意义。或者CP2K研究某些有机物高温热解

Uus/pMeC6H4-/キ

Observer-Aya 发表于 2026-3-14 14:23
老师请问做上千K温度的AIMD有什么问题吗，我一度打算用CP2K做个一千多个原子的973K的AIMD来着，后来因计 ...

熔盐或者熔融金属之类的体系热浴控温做一两千K温度的AIMD当然还是符合实际的，另外质谱（参见http://bbs.keinsci.com/thread-20479-1-1.html提到的Grimme文章）之类高能量、非热平衡的体系也可以有更高的温度。（但是像http://bbs.keinsci.com/thread-53811-1-1.html最后的表格那样强行把常温量化计算所得反应速率外推到近万K的行为就比较可疑了，这种等离子体似的极端温度下连电子是否还与原子核相互作用、Born-Oppenheimer近似是否成立都不好说的。）

我说的那篇文章是理论计算了假想的“某某原子掺杂的某某纳米体系吸附某某小分子”体系，打算用AIMD说明体系不容易受热自发分解的，此时还跑出一两千K的轨迹且没有展示完整结构截图就非常离谱了。或者也很有可能那几条曲线是别的物理量但图例被误标成了温度……

Observer-Aya

SharkYYX2025 发表于 2026-3-14 19:24
好吓人的计算量，拿GFN1-xTB没准能行，DFT做不动吧
高温AIMD不是说一定不行，我个人理解是看物理意义。 ...

我们是做熔盐腐蚀的实验组，我本来打算用castep算个一百多个原子体系来研究腐蚀扩散行为的，但合作的指导老师告诉我至少得上千原子而且cp2k最多效率提三倍根本没法算，算了也不一定有想要的结果，最后就放弃了

wzkchem5

Observer-Aya 发表于 2026-3-15 15:12
我们是做熔盐腐蚀的实验组，我本来打算用castep算个一百多个原子体系来研究腐蚀扩散行为的，但合作的指导 ...

这种情况还有一个问题，就是实验上在秒级以上，甚至天级的时间尺度上观察到的腐蚀，计算在ps时间尺度上很可能完全观察不到，甚至腐蚀一个原子都看不到。又不好盲目大幅提高模拟温度来加速腐蚀，那样主要的腐蚀机理可能就变了

cyh时代

之前有人找我代算，就要这种等值面图，不明白有什么意义，这种图用VMD或者Chemcraft不好嘛，真的把话说得让人听不懂，让人听不明白，就以为自己很牛，很神，其实对于懂计算的，如同小丑

SharkYYX2025

cyh时代 发表于 2026-3-23 00:04
之前有人找我代算，就要这种等值面图，不明白有什么意义，这种图用VMD或者Chemcraft不好嘛，真的把话说得让 ...

哈哈哈哈，他意图故弄玄虚，实际上成了内行人的乐子

J.Liu

这么看来，PBE算弱相互作用不加D3矫正的都算比较强的了

hjlyyc0999

一般实验的审稿人也不太看计算，就这样过去啦

SharkYYX2025

一楼需要审核，并且我发现更新一次会少附件在这里继续更新，不知行不行

—————————新连载—————————
1. 神奇结合能竟是G>0
期刊 Advanced Functional Materials  通讯Mao Liang   doi.org/10.1002/adfm.75159
题目Multifunctional Pyrene-Carbazole Integration as Flexible-Skeletonized Hole Contact for Stable Inverted Perovskite Solar Cells

• B3LYP 没加色散校正。这算是小问题。他写计算偶极矩，没写计算弱相互作用。但这里对于量子化学计算就这一句话描述，默认弱相互作用也是这个计算级别。
• ESP在图S20写成EPS，偶极矩单位debye拼错（正文讨论也拼错了），ESP色彩轴看不清刻度上下限
• 图1c, IGMH没放颜色轴也没声明数值。相互作用能为啥用吉布斯自由能变（不用电子能量E差值），而且大于0。但正文讨论写的是负的和正的，也不知道到底正还是负。binding Gibbs free energy是什么新说法？
  对于这种实验为主的文章，咱尚且不要求什么泛函基组精度够不够，或者BSSE
• 图1d，动力学模拟不给图例，看不出这想表达啥