给博主留言

如果你有相关的技术问题,或对博主有任何建议、想法,欢迎在此留言。

不便留言的,可通过邮箱联系博主:xianbao.d#gmail.com(将#替换成@)

说明

1. 讨论技术问题,请直接留言,不要用邮件,谢谢!

2. 由于博主工作较忙,难以做到及时回复,还望理解。

添加新评论 (无需注册,可直接评论)

已有 146 条评论

  1. 云飞扬

    段老师好,您的这个帖子http://www.52souji.net/how-to-convert-hexagonal-cell-to-cubic-cell.html,图片不能显示了,能更新下吗,谢谢

  2. 王雅兰

    老师您好,感谢您的分享,我想请问一下老师有没有了解用lammps做油水分离这一块,主要是利用水与基板之间的亲疏水性来实现油水分离,期待您的回复,感谢!

  3. wyf

    老师 您好!非常感谢您整理出来的资料,受益颇多。目前我在用lammps做模拟的过程中总是会在跑一段时间之后出现outof range-cant compute pppm,尝试修改了时间步长以及neighbor bin的值,可还是出错。希望老师可以提供一些解决的思路,谢谢老师!

  4. LSW

    我是初学者,感谢您提供的资源。谢谢,希望能与您多多交流。不断进步。

  5. 中科院MD

    图片看不了啦,请及时更下下哈,新手留~~辛苦

  6. candy

    宝哥,MD初学者,最近遇到两个大问题。

    我用atomsk建立多晶合金,多晶里只能包含两种元素,随机比例替换的。按照你的思路是建立一个晶胞再进行扩大,这样的前提是合金里面一般有多种元素,而且好像是随机比例添加进去金属基体里面的,这样的话,如果只知道元素比例,该如何建模? 我做纳米切削,想加上tip4p水分子模型模拟切削液,但是出现各种各样的问题,你能不能写一个切削加水的博文出来给我们这群菜鸟解疑?
  7. 科研路上人

    博主您好,我是做lammps切削的,最近仿真出现了一些问题,如果您在百忙之中可以看一下万分感谢,问题如下:
    1.我使用nvt进行300k下驰豫,仿真时出现温度由低到高升致300k,而非初始在300K附近大幅震荡之后平稳,请问这种情况正确么?
    2.在上述nvt系综下进行非固定层单精铁驰豫,出现固定层与非固定层结合处少量断裂,之后愈合的现象,请问造成这种现象的可能的原因?
    3.如果可以,我想将我的仿真图发给您,这样您可以更好的判断。
    万分感谢

  8. 雕刻时光

    段老师,你好!感谢您分享这么多宝贵的资源!想请教你一些问题:
    我现在在尝试做Zr-X二元合金2NN MEAM势函数的优化,想优化里面的屏蔽参数(可能会优化4到6个),我很详细的阅读了你的博士论文,因为2NN MEAM势函数中有个Zrx.meam文件,手动去改变Zrx.meam文件中屏蔽参数的值来达到优化是不可行的,我想请教的是,如何使用MATLAB、pso优化算法来实现最终参数的优化,由于课题组之前没有过相关的经验与探索,所以不了解具体的操作细节,可否请老师解答一下我的疑惑,不胜感激!

    1. 关于势参数优化的过程,在我的博士论文里有写。PSO的核心代码网络上可以下到,我博客里也有(见下链接)。你自己需要些的主要就是流程控制代码。
      http://www.52souji.net/matlab-code-of-particle-swarm-optimization.html

  9. 博主,你好!在铜晶体弥散分布三氧化二铝(随便一个含量),请问如何建模且有没有放入lammps计算的脚本

  10. 逸秋雨

    博主有没有将POSCAR文件转换为lmp的python脚本呀?

  11. EMK

    SQ076478684020180825002 的

  12. gerrywong

    博主您好,我在用lammps计算一个161616区域内的原子的rdf曲线,之前计算的是整体的rdf,现在我想计算其中某个原子的周边原子对该原子的rdf曲线,所有原子都是一样的,我尝试修改compute rdf命令,但是一直报错,是不是我的想法有问题

  13. 大工的硕士

    老师,你好,我最近在做有关激光辅热加工石英玻璃的仿真,有两个问题请教一下:
    1:我想设置一个随时间变化的热源,我考虑用region 加上move,你看可行吗?
    2.:有关于fix heat 这个命令,我想用它模拟激光强度,您看可以吗?

  14. huajia

    老师,您好,我刚开始接触这个软件,现在想用这个软件计算Cu的不同晶界类型的晶界能。可是现在我不知道要怎么构建不同的晶界类型(对称晶界、非对称晶界、扭转晶界),甚至也不清楚要用什么来算晶界能?请老师帮帮我。(ps:想用的是双晶模型及Born-Von Karman 3D periodic boundary conditions)

  15. 2NN MEAM势函数如何在短程部分用ZBL修正,我看很多文献都用这个势函数与ZBL函数耦合,我想用这种势函数来做辐照级联,但是短程部分我一直弄不好,想请问下大神
    谢谢

    1. 六六

      我也想了解这类问题。还请老师 指教

  16. 老师您好:请教您一个问题,就是我要计算一个不规则形状的气泡的体积大小,能用 compute voronoi/atom command 这个命令,然后再对它们进行求和吗?万分感谢!!

    1. 抱歉,没有处理过这类问题,请在网络上搜索。

  17. 小唐

    博主,能不能用matlab写一个键角分布函数的统计的程序,大致是通过循环找到离Si原子最近的两个O原子坐标,这样计算这三个原子之间的距离a,b,c,再通过余弦定理求得夹角,对lammps中的数据做处理!!

    1. 这个还是请自行解决吧。

  18. 小唐

    希望博主可以用matlab写一个键角分布函数的统计的小程序

  19. 一点点进步

    博主好,新手入门,想重复别人论文里的孪晶例子,这里建模使用lammps建模还是其他软件更方便?我在官网找到了一个晶界生成的例子,是不是用这个例子就可以?另外,如何在模型中建立椭圆形孔,我只在手册里找到了圆柱形孔的建模方法,求教,非常感谢!

    1. 可以通过编程实现复杂模型的构建。

  20. 奋斗

    我在进行多晶钛铝合金模型替换原子时,发现可以用Nb替换一定比例的Ti原子,但是没法用Nb替换一定比例的Al原子,这是什么呢

    1. 抱歉,你的描述不够准确,无法帮助到你。