| 查看: 787 | 回复: 0 | |||
yoghurt117木虫 (正式写手)
|
[求助]
盒子，盒子 ：（
|
|
最近看cpmd官网的教程，其中一个bulk系统的动力学计算，其输入文件是： &INFO md at T=300K with thermostats &END &CPMD MOLECULAR DYNAMICS RESTART WAVEFUNCTIONS COORDINATES VELOCITIES LATEST TRAJECTORY SAMPLE XYZ 10 RESCALE OLD VELOCITIES NOSE IONS 400.0d0 3000.0d0 NOSE ELECTRONS 0.0080d0 10000.0d0 MAXSTEP 5000 TIMESTEP 4.0 EMASS 400.0 MEMORY BIG &END &SYSTEM SYMMETRY 1 ANGSTROM CELL 9.95991 1.0 1.0 0.0 0.0 0.0 DUAL 5.0 CUTOFF 25.0 CHARGE 1 &END &DFT FUNCTIONAL BLYP GC-CUTOFF 5.0d-6 &END &ATOMS ISOTOPES 16.00 2.016 *O_VDB_BLYP.psp FORMATTED LMAX=P 33 0.7327 -4.2440 2.1886 -0.8347 -4.7227 4.3691 2.6518 -1.5714 -0.1309 3.2623 -4.4641 2.6462 -4.7361 4.8510 -2.0459 -2.3920 0.3963 3.9231 -0.5493 -4.5740 -2.0499 2.2782 -3.3931 -4.3667 -3.5142 -4.3400 2.4342 -3.8442 0.5428 -2.4979 2.6625 2.2735 2.5634 4.0301 2.4810 -3.1204 -1.3385 -0.8110 -0.2060 3.3646 -4.3723 -0.2574 0.7415 -0.6233 -1.8793 4.5585 3.6260 4.3279 1.4453 3.8230 -1.3437 0.7650 1.2950 4.2396 -2.5713 3.1043 -4.8839 -4.1207 3.2564 1.0086 0.4496 3.1128 1.4121 4.0998 0.0228 4.2341 -0.5407 -2.4217 -3.9645 -4.5374 -1.1430 -0.2440 -1.6032 1.7664 0.5732 -3.1993 -1.5826 -4.4357 0.5592 -1.4779 1.7951 -1.8413 2.8694 -2.0738 0.9314 1.6635 -3.0237 4.8851 0.6455 1.8791 -4.7775 -3.9950 -4.5931 -2.9360 -3.2201 -0.0728 1.8651 -1.3161 4.1881 *H_VDB_BLYP.psp FORMATTED LMAX=P 67 0.1804 -4.3091 3.0196 0.4583 -3.4352 1.6684 -0.8907 -3.8505 4.8551 -1.4470 -5.3849 4.8012 2.6813 -2.5699 -0.0836 3.5705 -1.2189 -0.3089 2.2641 -4.5114 2.6070 3.6143 -5.2567 3.1440 -5.2043 3.9908 -2.2480 -5.2257 5.3344 -1.3202 -2.8338 0.3816 3.0261 -2.6987 -0.3905 4.4587 0.2035 -5.0868 -1.6374 -1.2460 -5.2114 -2.3791 3.1464 -3.8096 -4.0969 1.5234 -3.8864 -3.9344 -4.2719 -3.7659 2.7445 -2.7608 -4.2804 3.0891 -3.5930 -0.0389 -3.2716 -3.0803 1.1482 -2.2743 3.3247 1.5747 2.2926 3.0265 2.7969 3.3338 4.7065 1.8430 -2.7522 3.1319 2.0416 -3.1370 -1.9556 -1.5794 -0.0367 -1.7355 0.0251 0.1726 3.1616 -4.3000 0.7191 2.7360 -5.0244 -0.6811 1.5443 -1.0463 -1.4591 -0.0408 -0.7072 -1.2621 5.5462 3.4712 4.3092 4.1734 3.2213 5.1573 1.3289 3.2796 -0.5124 1.3588 3.2301 -2.1443 1.2158 1.6650 3.4272 -0.2192 1.4615 4.1796 -2.3393 2.1786 -5.1828 -2.3958 3.1917 -3.9033 -4.0469 3.9980 1.6755 -4.8820 3.4416 0.3872 1.3081 2.8588 1.8577 0.2815 4.0896 1.5449 4.4722 0.2911 3.3456 4.8505 -0.1433 4.8736 -0.5906 -3.0520 -3.1897 -0.0056 -1.6175 -3.7060 -4.2505 -0.6992 -1.0929 -3.9724 -1.9523 -0.0834 -2.3879 2.3426 0.8019 -0.7686 2.1804 0.9365 -3.8725 -2.3207 -4.3898 -2.2782 -1.9717 -4.4453 1.3775 -1.3132 1.2443 -0.2523 -1.2143 1.2736 -1.1466 2.2371 -2.4168 -2.0815 2.6232 -1.1348 1.0048 0.8099 -2.5081 0.8075 1.4570 -3.9943 5.3467 1.4906 1.6095 4.9527 -0.0223 1.1379 -4.4389 -4.2859 -3.6983 -5.0180 -4.7996 -5.1361 -2.9873 -3.6787 0.8143 -2.9264 -3.9045 -0.8018 1.8677 -2.1810 4.7284 2.7693 -0.8629 4.0578 1.2248 -1.2816 3.3950 &END 我把原子和坐标部分做成了xyz格式的文件 100 O 0.7327 -4.2440 2.1886 O -0.8347 -4.7227 4.3691 O 2.6518 -1.5714 -0.1309 O 3.2623 -4.4641 2.6462 O -4.7361 4.8510 -2.0459 O -2.3920 0.3963 3.9231 O -0.5493 -4.5740 -2.0499 O 2.2782 -3.3931 -4.3667 O -3.5142 -4.3400 2.4342 O -3.8442 0.5428 -2.4979 O 2.6625 2.2735 2.5634 O 4.0301 2.4810 -3.1204 O -1.3385 -0.8110 -0.2060 O 3.3646 -4.3723 -0.2574 O 0.7415 -0.6233 -1.8793 O 4.5585 3.6260 4.3279 O 1.4453 3.8230 -1.3437 O 0.7650 1.2950 4.2396 O -2.5713 3.1043 -4.8839 O -4.1207 3.2564 1.0086 O 0.4496 3.1128 1.4121 O 4.0998 0.0228 4.2341 O -0.5407 -2.4217 -3.9645 O -4.5374 -1.1430 -0.2440 O -1.6032 1.7664 0.5732 O -3.1993 -1.5826 -4.4357 O 0.5592 -1.4779 1.7951 O -1.8413 2.8694 -2.0738 O 0.9314 1.6635 -3.0237 O 4.8851 0.6455 1.8791 O -4.7775 -3.9950 -4.5931 O -2.9360 -3.2201 -0.0728 O 1.8651 -1.3161 4.1881 H 0.1804 -4.3091 3.0196 H 0.4583 -3.4352 1.6684 H -0.8907 -3.8505 4.8551 H -1.4470 -5.3849 4.8012 H 2.6813 -2.5699 -0.0836 H 3.5705 -1.2189 -0.3089 H 2.2641 -4.5114 2.6070 H 3.6143 -5.2567 3.1440 H -5.2043 3.9908 -2.2480 H -5.2257 5.3344 -1.3202 H -2.8338 0.3816 3.0261 H -2.6987 -0.3905 4.4587 H 0.2035 -5.0868 -1.6374 H -1.2460 -5.2114 -2.3791 H 3.1464 -3.8096 -4.0969 H 1.5234 -3.8864 -3.9344 H -4.2719 -3.7659 2.7445 H -2.7608 -4.2804 3.0891 H -3.5930 -0.0389 -3.2716 H -3.0803 1.1482 -2.2743 H 3.3247 1.5747 2.2926 H 3.0265 2.7969 3.3338 H 4.7065 1.8430 -2.7522 H 3.1319 2.0416 -3.1370 H -1.9556 -1.5794 -0.0367 H -1.7355 0.0251 0.1726 H 3.1616 -4.3000 0.7191 H 2.7360 -5.0244 -0.6811 H 1.5443 -1.0463 -1.4591 H -0.0408 -0.7072 -1.2621 H 5.5462 3.4712 4.3092 H 4.1734 3.2213 5.1573 H 1.3289 3.2796 -0.5124 H 1.3588 3.2301 -2.1443 H 1.2158 1.6650 3.4272 H -0.2192 1.4615 4.1796 H -2.3393 2.1786 -5.1828 H -2.3958 3.1917 -3.9033 H -4.0469 3.9980 1.6755 H -4.8820 3.4416 0.3872 H 1.3081 2.8588 1.8577 H 0.2815 4.0896 1.5449 H 4.4722 0.2911 3.3456 H 4.8505 -0.1433 4.8736 H -0.5906 -3.0520 -3.1897 H -0.0056 -1.6175 -3.7060 H -4.2505 -0.6992 -1.0929 H -3.9724 -1.9523 -0.0834 H -2.3879 2.3426 0.8019 H -0.7686 2.1804 0.9365 H -3.8725 -2.3207 -4.3898 H -2.2782 -1.9717 -4.4453 H 1.3775 -1.3132 1.2443 H -0.2523 -1.2143 1.2736 H -1.1466 2.2371 -2.4168 H -2.0815 2.6232 -1.1348 H 1.0048 0.8099 -2.5081 H 0.8075 1.4570 -3.9943 H 5.3467 1.4906 1.6095 H 4.9527 -0.0223 1.1379 H -4.4389 -4.2859 -3.6983 H -5.0180 -4.7996 -5.1361 H -2.9873 -3.6787 0.8143 H -2.9264 -3.9045 -0.8018 H 1.8677 -2.1810 4.7284 H 2.7693 -0.8629 4.0578 H 1.2248 -1.2816 3.3950 然后转化成mol用gaussview量了一下，这个系统各个方向上平均10个多埃，最大的有16埃，不知道symmetry部分的 CELL 9.95991 1.0 1.0 0.0 0.0 0.0 里面的9.995991×9.995991×9.995991埃这样的格子怎么放下的，初学不太明白，请教一下，谢谢。 |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有56人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有9人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
关于Filecode分析方法
已经有13人回复
-
时间戳又变了8-15
已经有9人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有8人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有9人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有6人回复
找到一些相关的精华帖子，希望有用哦~
求问关于真空盒子的作用。
已经有12人回复
- 跳出盒子--一本不错的经济管理书籍推荐【转载】
已经有21人回复
- 求助一个带有真空层盒子的建立问题……
已经有5人回复
- 请教做MD计算时，蛋白质分子有一部分跑出水盒子
已经有9人回复
- 【求助】怎样才能使整个盒子充满水,用gromacs
已经有6人回复
- 【求助】溶剂盒子的密度？
已经有5人回复
- 【求助】dlpoly跑动力学生成的体系，分散到盒子的顶点上了，这是怎么一回事？应该怎么
已经有5人回复
- 【交流】将8个完全相同的球放到3个不同的盒子中
已经有20人回复
- 【求助】如何把刚出锅的牛肉快速、安全的分装到小盒子中？
已经有5人回复
- 【求助】水盒子问题
已经有17人回复
- 【求助】做MD时，溶剂盒子应该加多大算合适呢？
已经有6人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
