计算化学公社

 找回密码 Forget password
 注册 Register
Views: 992|回复 Reply: 1
打印 Print 上一主题 Last thread 下一主题 Next thread

[Lammps] Lammps如何通过variables 方法设置初始速度

[复制链接 Copy URL]

26

帖子

0

威望

1370

eV
积分
1396

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
想要修改lammps动力学过程中的初速度。读取文件里面有速度信息如下:
Velocities

     1    -4.47890057    -1.17385902    -0.94094528
     2     3.05323177    12.14933185    19.67673848
     3    -6.42438377     4.91325661    -1.91292388
     4    -5.11973815     5.33623365     2.04415182
     5     4.58948674     3.53064422    -9.57939583
     6   -13.45035851   -16.34786554    -6.46497131
     7    -2.62383340    19.62462569   -20.26197574
     8    -4.83349727    -0.15072031    11.32092111
     9    -0.49705471     4.59045988     8.84395494
    10   -12.18502954   -29.19157946    13.11627328
    11   -10.10542350     0.63133858    -1.99135360
    12    -3.61184068     6.42655097     9.69511583
    13    10.57703558     3.01256147     9.57409972
    14    -5.61177974     3.48757832   -12.87966363
    15     3.61537634    -5.96994269     1.92956211

这里想通过命令set atom 100 vx 0.0 vy 0.0 vz -1.0 的方式达到设置初速度的效果,想问下具体的文件格式和这三个参数vx, vy, vz的具体命令行设置

26

帖子

0

威望

1370

eV
积分
1396

Level 4 (黑子)

2#
 楼主 Author| 发表于 Post on 2023-8-24 10:32:35 | 只看该作者 Only view this author
Atoms

  1 molecule-tag atom-type q x y z nx ny nz  (nx,ny,nz are optional -
  ...                                    see "true flag" input command)
  ...               
  N molecule-tag atom-type q x y z nx ny nz  (N = # of atoms)

Velocities

  1 vx vy vz
  ...
  ...               
  N vx vy vz                        (N = # of atoms)
可以按照上面的方式写入conf_new.lmp文件,通过read_data方法读取。参考网页https://docs.lammps.org/2001/data_format.html

本版积分规则 Credits rule

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

GMT+8, 2025-8-14 14:17 , Processed in 0.137458 second(s), 20 queries , Gzip On.

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