| 查看: 374 | 回复: 0 | ||
joyce79928铜虫 (小有名气)
|
[求助]
如何做出力對於偏壓的圖利用transeista
|
|
在'Current-induced energy barrier suppression for electromigration from first principles'这篇文献中 里面有做出力对于偏压的图 这应该要如何做出呢? 在transeista的output档里会有计算出来的力作用于每一个原子上的说明 siesta: Atomic forces (eV/Ang): siesta: 1 0.119023 -0.068722 -6.937691 siesta: 2 -0.119064 -0.068720 -6.937699 siesta: 3 -0.000016 -0.066163 -6.850128 siesta: 4 -0.000017 -0.007928 -6.901799 siesta: 5 -0.057327 0.033093 -6.850128 siesta: 6 0.057290 0.033096 -6.850128 siesta: 7 0.006852 0.003977 -6.901803 siesta: 8 -0.000017 0.137475 -6.937697 siesta: 9 -0.006889 0.003976 -6.901804 siesta: 10 -0.000009 0.065853 -23.295369 siesta: 11 -0.097958 -0.056551 -22.815127 siesta: 12 0.097953 -0.056544 -22.815131 siesta: 13 -0.107622 0.062138 -23.309638 siesta: 14 0.000003 0.113097 -22.815117 siesta: 15 0.107606 0.062125 -23.309633 siesta: 16 -0.057044 -0.032906 -23.295355 siesta: 17 0.057017 -0.032900 -23.295352 siesta: 18 0.000003 -0.124246 -23.309636 siesta: 19 0.000023 0.000008 -45.383366 siesta: 20 0.813486 0.507253 -45.206940 siesta: 21 -0.813424 0.507268 -45.206959 siesta: 22 0.846037 0.450818 -45.206932 siesta: 23 0.000021 0.000019 -43.763899 siesta: 24 -0.845994 0.450820 -45.206978 siesta: 25 0.032612 -0.958114 -45.206951 siesta: 26 -0.032557 -0.958123 -45.206971 siesta: 27 -0.000015 0.000012 -44.835979 siesta: 28 0.898675 -0.518807 -76.236054 siesta: 29 -0.898588 -0.518817 -76.236049 siesta: 30 0.000055 -1.171231 -78.496302 siesta: 31 0.000114 0.308294 -76.701739 siesta: 32 -1.014221 0.585563 -78.496266 siesta: 33 1.014356 0.585575 -78.496269 siesta: 34 -0.266946 -0.154181 -76.701738 siesta: 35 0.000059 1.037629 -76.236051 siesta: 36 0.267111 -0.154146 -76.701732 siesta: 37 -0.000041 -0.468888 38.222911 siesta: 38 -1.617373 -0.933794 34.644977 siesta: 39 1.617246 -0.933805 34.644977 siesta: 40 -1.020184 0.588848 37.198775 siesta: 41 -0.000061 1.867487 34.644648 siesta: 42 1.020085 0.588933 37.198796 siesta: 43 0.406184 0.234590 38.222877 siesta: 44 -0.406286 0.234599 38.222869 siesta: 45 -0.000044 -1.177913 37.198790 siesta: 46 0.001396 0.001810 -0.008797 siesta: 47 -0.005711 -0.001807 -0.001017 siesta: 48 0.007792 -0.001778 0.004357 siesta: 49 -0.003466 0.001849 0.005785 siesta: 50 0.000001 -0.154646 -38.078172 siesta: 51 -1.668779 -0.963666 -34.684690 siesta: 52 1.668936 -0.963644 -34.684797 siesta: 53 -1.160005 0.669567 -37.942066 siesta: 54 0.000068 1.927018 -34.684697 siesta: 55 1.160012 0.669676 -37.942044 siesta: 56 0.133920 0.077309 -38.078129 siesta: 57 -0.133663 0.077205 -38.077870 siesta: 58 -0.000089 -1.339047 -37.942073 siesta: 59 -0.000127 0.000288 76.677987 siesta: 60 -0.865981 -0.738084 76.693409 siesta: 61 0.865619 -0.738128 76.693619 siesta: 62 -1.072092 -0.380625 76.693510 siesta: 63 -0.000101 -0.000008 79.863168 siesta: 64 1.071756 -0.380643 76.693596 siesta: 65 -0.206486 1.118754 76.693514 siesta: 66 0.206192 1.118928 76.693464 siesta: 67 -0.000138 0.000019 75.519516 siesta: 68 -0.732727 0.423028 41.825032 siesta: 69 0.733775 0.422942 41.825102 siesta: 70 0.000498 1.948940 40.103647 siesta: 71 0.000462 -0.440445 42.048351 siesta: 72 1.688620 -0.974934 40.103842 siesta: 73 -1.687592 -0.974941 40.103677 siesta: 74 0.381708 0.219787 42.048352 siesta: 75 0.000564 -0.846919 41.825163 siesta: 76 -0.380688 0.219778 42.048316 siesta: 77 -0.000029 -0.114614 25.032750 siesta: 78 -0.377372 -0.217859 24.253679 siesta: 79 0.377310 -0.217852 24.253659 siesta: 80 -0.233848 0.135022 24.932312 siesta: 81 -0.000025 0.435752 24.253679 siesta: 82 0.233774 0.135013 24.932314 siesta: 83 0.099246 0.057323 25.032754 siesta: 84 -0.099306 0.057328 25.032747 siesta: 85 -0.000034 -0.269956 24.932300 siesta: 86 0.000106 -0.000061 8.979956 siesta: 87 -0.038638 -0.150365 8.943876 siesta: 88 0.038851 -0.150368 8.943879 siesta: 89 -0.149431 0.041536 8.943887 siesta: 90 0.000100 -0.000055 8.731317 siesta: 91 0.149648 0.041533 8.943881 siesta: 92 -0.110689 0.108644 8.943882 siesta: 93 0.110897 0.108646 8.943870 siesta: 94 0.000102 -0.000060 9.023879 siesta: ---------------------------------------- siesta: Tot 0.004457 -0.002760 -4.563117 但是在input档里面有一个对力的收敛的设定 MD.MaxForceTol=0.04 eV/Ang(预设值) 所以最后我必定会得到一个介于-0.04~0.04 得力作用在原子上 但是从一开始提到的文献的图可以看到，力在不同偏压下作用在原子上的利大于0.04甚至到达0.5 这是怎么办到的? 还是我应该要设置MD.MaxForceTol=0.5 eV/Ang? 拜托各位大侠求助阿... 最近我放上来的帖子都没有人回应... 很无助阿... |
回复此楼» 猜你喜欢
300求调剂
已经有9人回复
-
291求调剂
已经有7人回复
-
求调剂推荐
已经有7人回复
-
289 分105500药学专硕求调剂(找B区学校)
已经有4人回复
-
0854求调剂
已经有13人回复
-
初试324 中药学 一志愿天中医 求调剂
已经有3人回复
-
药学求调剂
已经有14人回复
-
327求调剂
已经有27人回复
-
急需调剂
已经有5人回复
-
271求调剂
已经有33人回复
