24小时热门版块排行榜    

查看: 360  |  回复: 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?

拜托各位大侠求助阿...
最近我放上来的帖子都没有人回应...
很无助阿...
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 joyce79928 的主题更新
信息提示
请填处理意见