| 查看: 1315 | 回复: 2 | |||
| 【悬赏金币】回答本帖问题，作者大力2010将赠送您 10 个金币 | |||
[求助]
请教大师，vasp做分子动力学怎样设置INCAR，当前提示：very serious probelem！
|
|||
|
请教大师，vasp做分子动力学怎样设置incar，当前提示：very serious probelem！ 我想采用从头算分子动力学（aimd）模拟的熔融淬火和几何优化，采用正则系综（nvt）下 nose进行温度调节。将系统升温至 8000 k，平衡 1 ps， 时间步长为 1 fs，，然后以 1.6×104 k/ps 的冷却速率将系统温度从 8000 k 快速冷却到 1 k，随后对系统进行基于共轭梯度法的几何优化，令原子位置完全弛豫。截断能设置为 500 e v，并对布里渊区进行收敛计算。 global parameters istart = 1 (read existing wavefunction; if there) ispin = 1 (non-spin polarised dft) # icharg = 11 (non-self-consistent: gga/lda band structures) lreal = .false. (projection operators: automatic) encut = 550 (cut-off energy for plane wave basis set, in ev) prec = normal (precision level) lwave = .true. (write wavecar or not) lcharg = .true. (write chgcar or not) addgrid= .true. (increase grid; helps gga convergence) # lvtot = .true. (write total electrostatic potential into locpot or not) # lvhar = .true. (write ionic + hartree electrostatic potential into locpot or not) # nelect = (no. of electrons: charged cells; be careful) # lplane = .true. (real space distribution; supercells) # npar = 4 (max is no. nodes; don't set for hybrids) # nwrite = 2 (medium-level output) # kpar = 2 (divides k-grid into separate groups) # ngx = 500 (fft grid mesh density for nice charge/potential plots) # ngy = 500 (fft grid mesh density for nice charge/potential plots) # ngz = 500 (fft grid mesh density for nice charge/potential plots) electronic relaxation ismear = 0 sigma = 0.05 ediff = 1e-06 molecular dynamics ibrion = 0 (activate md) nsw = 1000 (max ionic steps) ediffg = -1e-02 (ionic convergence; ev/a) potim = 3 (timestep in fs) smass = 3 (md algorithm: -3-microcanonical ensemble; 0-canonical ensemble) ! tebeg = 100 (start temperature k) ! teend = 100 (final temperature k) ! mdalgo = 1 (andersen thermostat) ! isym = 0 (switch symmetry off) nwrite = 0 (for long md-runs use nwrite=0 or nwrite=1) |
回复此楼» 猜你喜欢
中国科学技术大学材料与化工281求调剂，有科研和获奖经历
已经有12人回复
-
申博
已经有4人回复
-
085701环境工程 求调剂
已经有7人回复
-
考研282分求调剂，接受跨专业
已经有5人回复
-
一志愿南京大学资源与环境专业283分求调剂
已经有3人回复
-
欢迎211本科同学，过A区国家线，A区非偏远一本，交叉学科课题组
已经有14人回复
-
【求调剂】293分环境工程求调剂材料/化工，服从调剂，抗压能力强！
已经有6人回复
-
347求调剂
已经有6人回复
-
274环境工程求调剂
已经有4人回复
-
环境工程专硕307 一志愿211 四六级已过 求调剂
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2022-10-06 17:11:51
3楼2022-11-04 18:53:40
6