亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時(shí)熱門(mén)版塊排行榜    

查看: 3559  |  回復(fù): 24

bjmaowei25

銅蟲(chóng) (小有名氣)

[求助] 關(guān)于vasp收斂和優(yōu)化

請(qǐng)問(wèn)一下，關(guān)于vasp收斂，用能量來(lái)衡量，到什么程度可以算是收斂呢？
另外， 如果剛開(kāi)始精度不高，初始的晶體結(jié)構(gòu)不好，需要的時(shí)間長(zhǎng)，收斂后再優(yōu)化，將CONTCAR
復(fù)制到POSCAR在原目錄下運(yùn)行，其它的文件如WAVECAR 需要?jiǎng)h除嗎？
回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

第一性原理

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
追求真理
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖
回帖置頂 ( 共有11個(gè) )

bjmaowei25

銅蟲(chóng) (小有名氣)
引用回帖:
樓: Originally posted by cavediger at 2012-02-25 10:38:38:
1.收斂精度要看你的要求和計(jì)算的性質(zhì)，不知道時(shí)就先用默認(rèn)值
2.優(yōu)化收斂與否就是看是否達(dá)到了你的要求，一般可以取少數(shù)k點(diǎn)，待收斂后，增加k點(diǎn)，cp CONTCAR POSCAR，再進(jìn)一步優(yōu)化，此時(shí)不需要?jiǎng)h除波函數(shù)文件。
...

謝謝你的回復(fù)。我現(xiàn)在初步設(shè)計(jì)了一個(gè)包含80個(gè)原子的系統(tǒng)。為了收斂，我降低了精度，雖然運(yùn)行了幾個(gè)小時(shí)走了60步，自動(dòng)停下，但是不知到為什么看不到CONTCAR文件是空文件呢？這樣我沒(méi)法再優(yōu)化哦？
能問(wèn)下什么原因嗎？

附上關(guān)鍵的INCAR KPOINTS  OSZICAR  文件
INCAR
SYSTEM   = Y2O3
ENCUT    = 400
ISTART   = 0
ICHARG   = 2
ISMEAR   = 1
PREC     = M
EDIFF    = 1E-3
LREAL    = Auto
IALGO    = 48
POTIM    = 0.2


KPOINTS
Automatic generation
0
Monkhorst-Pack
4  4  4
0  0  0

OUTCAR
RMM:   7    -0.658806285728E+03   -0.86957E+02   -0.51125E+02  2380   0.463E+01
RMM:   8    -0.711753703723E+03   -0.52947E+02   -0.34621E+02  2380   0.264E+01
RMM:   9    -0.765155909509E+03   -0.53402E+02   -0.43352E+02  8456   0.207E+01
RMM:  10    -0.770473179417E+03   -0.53173E+01   -0.68668E+01  9489   0.616E+00
RMM:  11    -0.771646952082E+03   -0.11738E+01   -0.11839E+01  9404   0.334E+00
RMM:  12    -0.772403475153E+03   -0.75652E+00   -0.75905E+00  9512   0.232E+00    0.602E+01
RMM:  13    -0.795468219229E+03   -0.23065E+02   -0.55753E+02  7672   0.229E+01    0.126E+02
RMM:  14    -0.784297414099E+03    0.11171E+02   -0.44606E+02  6445   0.372E+01    0.165E+02
RMM:  15    -0.925576628452E+03   -0.14128E+03   -0.10698E+02  7465   0.174E+01    0.157E+02
RMM:  16    -0.772467558930E+03    0.15311E+03   -0.10494E+02  6573   0.202E+01    0.933E+01
RMM:  17    -0.784112451468E+03   -0.11645E+02   -0.12001E+01  7248   0.651E+00    0.890E+01
RMM:  18    -0.762522804367E+03    0.21590E+02   -0.19893E+01  6871   0.800E+00    0.766E+01
RMM:  19    -0.755086985068E+03    0.74358E+01   -0.18866E+01  6869   0.816E+00    0.111E+02
RMM:  20    -0.716594261115E+03    0.38493E+02   -0.11472E+01  6914   0.618E+00    0.704E+01
RMM:  21    -0.725383541584E+03   -0.87893E+01   -0.67816E+00  7051   0.484E+00    0.372E+01
RMM:  22    -0.746875517125E+03   -0.21492E+02   -0.12390E+01  6711   0.653E+00    0.333E+01
RMM:  23    -0.743159964378E+03    0.37156E+01   -0.34325E+01  6321   0.131E+01    0.800E+01
RMM:  24    -0.740404697184E+03    0.27553E+01   -0.77952E+00  6800   0.589E+00    0.725E+01
RMM:  25    -0.732389229235E+03    0.80155E+01   -0.63120E+00  6359   0.597E+00    0.686E+01
RMM:  26    -0.728137506081E+03    0.42517E+01   -0.65576E+00  6337   0.592E+00    0.306E+01
RMM:  27    -0.717403372303E+03    0.10734E+02   -0.45245E+00  6112   0.499E+00    0.484E+01
RMM:  28    -0.722015309987E+03   -0.46119E+01   -0.54068E+00  6104   0.525E+00    0.406E+01
RMM:  29    -0.713191034048E+03    0.88243E+01   -0.46201E+00  5953   0.501E+00    0.586E+01
RMM:  30    -0.717174822875E+03   -0.39838E+01   -0.43570E+00  5988   0.462E+00    0.136E+01
RMM:  31    -0.721963584515E+03   -0.47888E+01   -0.62275E+00  5929   0.573E+00    0.918E+01
RMM:  32    -0.715230010864E+03    0.67336E+01   -0.54929E+00  5953   0.508E+00    0.427E+01
RMM:  33    -0.718179287344E+03   -0.29493E+01   -0.54127E+00  5884   0.525E+00    0.162E+01
RMM:  34    -0.724784336523E+03   -0.66050E+01   -0.55183E+00  5898   0.521E+00    0.304E+01
RMM:  35    -0.725978235871E+03   -0.11939E+01   -0.66878E+00  5923   0.584E+00    0.378E+01
RMM:  36    -0.727131711019E+03   -0.11535E+01   -0.65556E+00  5927   0.566E+00    0.267E+01
RMM:  37    -0.722455076618E+03    0.46766E+01   -0.65257E+00  5882   0.579E+00    0.626E+01
RMM:  38    -0.724982520383E+03   -0.25274E+01   -0.59768E+00  5919   0.532E+00    0.224E+01
RMM:  39    -0.720469536465E+03    0.45130E+01   -0.64423E+00  5878   0.578E+00    0.210E+01
RMM:  40    -0.721702122896E+03   -0.12326E+01   -0.56120E+00  5897   0.524E+00    0.378E+01
RMM:  41    -0.724593909156E+03   -0.28918E+01   -0.61926E+00  5879   0.569E+00    0.221E+01
RMM:  42    -0.725702167740E+03   -0.11083E+01   -0.62888E+00  5908   0.548E+00    0.552E+01
RMM:  43    -0.726196661007E+03   -0.49449E+00   -0.63342E+00  5904   0.568E+00    0.137E+01
RMM:  44    -0.725656611274E+03    0.54005E+00   -0.64076E+00  5906   0.552E+00    0.444E+01
RMM:  45    -0.725124954625E+03    0.53166E+00   -0.64324E+00  5870   0.571E+00    0.207E+01
RMM:  46    -0.726934804249E+03   -0.18098E+01   -0.66043E+00  5910   0.570E+00    0.163E+01
RMM:  47    -0.721457790593E+03    0.54770E+01   -0.65760E+00  5866   0.581E+00    0.348E+01
RMM:  48    -0.724662371458E+03   -0.32046E+01   -0.61471E+00  5890   0.545E+00    0.213E+01
RMM:  49    -0.720103519073E+03    0.45589E+01   -0.66294E+00  5872   0.579E+00    0.483E+01
RMM:  50    -0.722503219999E+03   -0.23997E+01   -0.57810E+00  5861   0.521E+00    0.245E+01
RMM:  51    -0.720274022729E+03    0.22292E+01   -0.65966E+00  5884   0.577E+00    0.249E+01
RMM:  52    -0.721287189816E+03   -0.10132E+01   -0.58589E+00  5878   0.527E+00    0.148E+01
RMM:  53    -0.723427925172E+03   -0.21407E+01   -0.64708E+00  5865   0.576E+00    0.156E+01
RMM:  54    -0.725591465900E+03   -0.21635E+01   -0.63635E+00  5894   0.549E+00    0.420E+01
RMM:  55    -0.725480868979E+03    0.11060E+00   -0.64124E+00  5871   0.575E+00    0.193E+01
RMM:  56    -0.725426845415E+03    0.54024E-01   -0.64454E+00  5897   0.553E+00    0.460E+01
RMM:  57    -0.725127788882E+03    0.29906E+00   -0.63177E+00  5863   0.569E+00    0.151E+01
RMM:  58    -0.726513720555E+03   -0.13859E+01   -0.63990E+00  5895   0.554E+00    0.382E+01
RMM:  59    -0.724530447404E+03    0.19833E+01   -0.64467E+00  5870   0.576E+00    0.207E+01
RMM:  60    -0.727517073333E+03   -0.29866E+01   -0.64768E+00  5873   0.560E+00
   1 F= -.72751707E+03 E0= -.72736524E+03  d E =-.455502E+00
追求真理
4樓2012-02-25 13:32:28
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

bjmaowei25

銅蟲(chóng) (小有名氣)
引用回帖:
樓: Originally posted by liliangfang at 2012-02-25 10:55:00:
1樓說(shuō)的很明白了，建議多看看說(shuō)明書(shū)，有中文的，比如侯柱峰老師講的就很好

謝謝版主的回復(fù)。我現(xiàn)在初步設(shè)計(jì)了一個(gè)包含80個(gè)原子的系統(tǒng)。為了收斂，我降低了精度，雖然運(yùn)行了幾個(gè)小時(shí)走了60步，自動(dòng)停下，但是不知到為什么看不到CONTCAR文件是空文件呢？這樣我沒(méi)法再優(yōu)化哦？
能問(wèn)版主這是什么原因嗎？、是不是INCAR參數(shù)設(shè)置的問(wèn)題？

附上關(guān)鍵的INCAR KPOINTS  OSZICAR  文件
INCAR
SYSTEM   = Y2O3
ENCUT    = 400
ISTART   = 0
ICHARG   = 2
ISMEAR   = 1
PREC     = M
EDIFF    = 1E-3
LREAL    = Auto
IALGO    = 48
POTIM    = 0.2


KPOINTS
Automatic generation
0
Monkhorst-Pack
4  4  4
0  0  0

OUTCAR
RMM:   7    -0.658806285728E+03   -0.86957E+02   -0.51125E+02  2380   0.463E+01
RMM:   8    -0.711753703723E+03   -0.52947E+02   -0.34621E+02  2380   0.264E+01
RMM:   9    -0.765155909509E+03   -0.53402E+02   -0.43352E+02  8456   0.207E+01
RMM:  10    -0.770473179417E+03   -0.53173E+01   -0.68668E+01  9489   0.616E+00
RMM:  11    -0.771646952082E+03   -0.11738E+01   -0.11839E+01  9404   0.334E+00
RMM:  12    -0.772403475153E+03   -0.75652E+00   -0.75905E+00  9512   0.232E+00    0.602E+01
RMM:  13    -0.795468219229E+03   -0.23065E+02   -0.55753E+02  7672   0.229E+01    0.126E+02
RMM:  14    -0.784297414099E+03    0.11171E+02   -0.44606E+02  6445   0.372E+01    0.165E+02
RMM:  15    -0.925576628452E+03   -0.14128E+03   -0.10698E+02  7465   0.174E+01    0.157E+02
RMM:  16    -0.772467558930E+03    0.15311E+03   -0.10494E+02  6573   0.202E+01    0.933E+01
RMM:  17    -0.784112451468E+03   -0.11645E+02   -0.12001E+01  7248   0.651E+00    0.890E+01
RMM:  18    -0.762522804367E+03    0.21590E+02   -0.19893E+01  6871   0.800E+00    0.766E+01
RMM:  19    -0.755086985068E+03    0.74358E+01   -0.18866E+01  6869   0.816E+00    0.111E+02
RMM:  20    -0.716594261115E+03    0.38493E+02   -0.11472E+01  6914   0.618E+00    0.704E+01
RMM:  21    -0.725383541584E+03   -0.87893E+01   -0.67816E+00  7051   0.484E+00    0.372E+01
RMM:  22    -0.746875517125E+03   -0.21492E+02   -0.12390E+01  6711   0.653E+00    0.333E+01
RMM:  23    -0.743159964378E+03    0.37156E+01   -0.34325E+01  6321   0.131E+01    0.800E+01
RMM:  24    -0.740404697184E+03    0.27553E+01   -0.77952E+00  6800   0.589E+00    0.725E+01
RMM:  25    -0.732389229235E+03    0.80155E+01   -0.63120E+00  6359   0.597E+00    0.686E+01
RMM:  26    -0.728137506081E+03    0.42517E+01   -0.65576E+00  6337   0.592E+00    0.306E+01
RMM:  27    -0.717403372303E+03    0.10734E+02   -0.45245E+00  6112   0.499E+00    0.484E+01
RMM:  28    -0.722015309987E+03   -0.46119E+01   -0.54068E+00  6104   0.525E+00    0.406E+01
RMM:  29    -0.713191034048E+03    0.88243E+01   -0.46201E+00  5953   0.501E+00    0.586E+01
RMM:  30    -0.717174822875E+03   -0.39838E+01   -0.43570E+00  5988   0.462E+00    0.136E+01
RMM:  31    -0.721963584515E+03   -0.47888E+01   -0.62275E+00  5929   0.573E+00    0.918E+01
RMM:  32    -0.715230010864E+03    0.67336E+01   -0.54929E+00  5953   0.508E+00    0.427E+01
RMM:  33    -0.718179287344E+03   -0.29493E+01   -0.54127E+00  5884   0.525E+00    0.162E+01
RMM:  34    -0.724784336523E+03   -0.66050E+01   -0.55183E+00  5898   0.521E+00    0.304E+01
RMM:  35    -0.725978235871E+03   -0.11939E+01   -0.66878E+00  5923   0.584E+00    0.378E+01
RMM:  36    -0.727131711019E+03   -0.11535E+01   -0.65556E+00  5927   0.566E+00    0.267E+01
RMM:  37    -0.722455076618E+03    0.46766E+01   -0.65257E+00  5882   0.579E+00    0.626E+01
RMM:  38    -0.724982520383E+03   -0.25274E+01   -0.59768E+00  5919   0.532E+00    0.224E+01
RMM:  39    -0.720469536465E+03    0.45130E+01   -0.64423E+00  5878   0.578E+00    0.210E+01
RMM:  40    -0.721702122896E+03   -0.12326E+01   -0.56120E+00  5897   0.524E+00    0.378E+01
RMM:  41    -0.724593909156E+03   -0.28918E+01   -0.61926E+00  5879   0.569E+00    0.221E+01
RMM:  42    -0.725702167740E+03   -0.11083E+01   -0.62888E+00  5908   0.548E+00    0.552E+01
RMM:  43    -0.726196661007E+03   -0.49449E+00   -0.63342E+00  5904   0.568E+00    0.137E+01
RMM:  44    -0.725656611274E+03    0.54005E+00   -0.64076E+00  5906   0.552E+00    0.444E+01
RMM:  45    -0.725124954625E+03    0.53166E+00   -0.64324E+00  5870   0.571E+00    0.207E+01
RMM:  46    -0.726934804249E+03   -0.18098E+01   -0.66043E+00  5910   0.570E+00    0.163E+01
RMM:  47    -0.721457790593E+03    0.54770E+01   -0.65760E+00  5866   0.581E+00    0.348E+01
RMM:  48    -0.724662371458E+03   -0.32046E+01   -0.61471E+00  5890   0.545E+00    0.213E+01
RMM:  49    -0.720103519073E+03    0.45589E+01   -0.66294E+00  5872   0.579E+00    0.483E+01
RMM:  50    -0.722503219999E+03   -0.23997E+01   -0.57810E+00  5861   0.521E+00    0.245E+01
RMM:  51    -0.720274022729E+03    0.22292E+01   -0.65966E+00  5884   0.577E+00    0.249E+01
RMM:  52    -0.721287189816E+03   -0.10132E+01   -0.58589E+00  5878   0.527E+00    0.148E+01
RMM:  53    -0.723427925172E+03   -0.21407E+01   -0.64708E+00  5865   0.576E+00    0.156E+01
RMM:  54    -0.725591465900E+03   -0.21635E+01   -0.63635E+00  5894   0.549E+00    0.420E+01
RMM:  55    -0.725480868979E+03    0.11060E+00   -0.64124E+00  5871   0.575E+00    0.193E+01
RMM:  56    -0.725426845415E+03    0.54024E-01   -0.64454E+00  5897   0.553E+00    0.460E+01
RMM:  57    -0.725127788882E+03    0.29906E+00   -0.63177E+00  5863   0.569E+00    0.151E+01
RMM:  58    -0.726513720555E+03   -0.13859E+01   -0.63990E+00  5895   0.554E+00    0.382E+01
RMM:  59    -0.724530447404E+03    0.19833E+01   -0.64467E+00  5870   0.576E+00    0.207E+01
RMM:  60    -0.727517073333E+03   -0.29866E+01   -0.64768E+00  5873   0.560E+00
   1 F= -.72751707E+03 E0= -.72736524E+03  d E =-.455502E+00
追求真理
5樓2012-02-25 13:34:00
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

bjmaowei25

銅蟲(chóng) (小有名氣)
引用回帖:
樓: Originally posted by liliangfang at 2012-02-25 14:04:35:
初步建粒子位置的時(shí)候盡可能的準(zhǔn)確些，然后再優(yōu)化，ISIF的值也要設(shè)置，論壇里面有別人寫(xiě)好的優(yōu)化步驟，你可以參考，另外80個(gè)粒子很多了，建議先優(yōu)化小一點(diǎn)的結(jié)構(gòu)

謝謝版主的回復(fù)，我建立的粒子坐標(biāo)，就是一個(gè)單元晶胞里面包含的的原子數(shù)目，而且用VESTA
打開(kāi)POSCAR 文件，就可以得到單元晶胞，是不是沒(méi)有設(shè)置ISIF值，就沒(méi)有CONTCAR文件呢，
因?yàn)樯鲜龀绦蜻\(yùn)行結(jié)果 CONTCAR 是空文件。

如果按你所指，優(yōu)化小一點(diǎn)的結(jié)構(gòu)，比如20個(gè)原子，1/4單元晶胞，或者說(shuō)是原胞吧，可以設(shè)置精度高后收斂，
可以得到的總能free  energy   TOTEN .  對(duì)于整個(gè)單元晶胞 的優(yōu)化是不是，寫(xiě)出被優(yōu)化的小結(jié)構(gòu)的對(duì)稱結(jié)構(gòu)，
最后算出總能是原胞的體積比倍數(shù)嗎？

謝謝！
追求真理
7樓2012-02-25 14:25:49
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

bjmaowei25

銅蟲(chóng) (小有名氣)
引用回帖:
樓: Originally posted by cavediger at 2012-02-25 14:50:47:
如果你的體系就一晶胞，那就不能再小了，否則有違體相材料結(jié)構(gòu)；
POTIM再取小點(diǎn)，譬如0.15，再試一下，看看收斂情況；
把K mesh取成2x2x2，優(yōu)化，收斂后調(diào)為4x4x4，cp CONTCAR POSCAR，再接著優(yōu)化，比較得到的 ...

謝謝版 你的回復(fù)，按照你的意識(shí)，右必要在設(shè)置ISIF=3, IBRION =1進(jìn)入INCAR中去驗(yàn)證收斂嗎？

另外i，關(guān)于原子坐標(biāo)的選取，對(duì)于MxOy的化合物，在一個(gè)單元晶胞內(nèi)選取的原子個(gè)數(shù)比M/O必須等于x/y?
追求真理
11樓2012-02-25 16:12:51
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

bjmaowei25

銅蟲(chóng) (小有名氣)
引用回帖:
樓: Originally posted by liliangfang at 2012-02-25 15:38:12:
那到不是，CONTCAR沒(méi)計(jì)算完一步就會(huì)更新一次，估計(jì)沒(méi)有可能是沒(méi)有計(jì)算完

是不是，每次計(jì)算得到E和dE就會(huì)更新一次CONTCAR？
現(xiàn)在每計(jì)算完一次CHG CHGCAR    CONTCAR  都是空文件！
難道POSCAR？
追求真理
12樓2012-02-25 16:57:48
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

bjmaowei25

銅蟲(chóng) (小有名氣)
引用回帖:
樓: Originally posted by cavediger at 2012-02-25 14:50:47:
如果你的體系就一晶胞，那就不能再小了，否則有違體相材料結(jié)構(gòu)；
POTIM再取小點(diǎn)，譬如0.15，再試一下，看看收斂情況；
把K mesh取成2x2x2，優(yōu)化，收斂后調(diào)為4x4x4，cp CONTCAR POSCAR，再接著優(yōu)化，比較得到的 ...

謝謝你的回復(fù)，根據(jù)你的建議 K mesh取成2x2x2， 并且加進(jìn)了ISIF IBRION  POTIM等參數(shù)。
很幸運(yùn)的得到了 CONTCAR 坐標(biāo)和POSCAR一樣。但是運(yùn)行到60步后能量不收斂！
檢查了POSCAR 沒(méi)有問(wèn)題啊，因?yàn)槲矣檬謱?xiě)的坐標(biāo)輸入VESTA能打開(kāi)正確的晶體結(jié)構(gòu)。
我貼出INCAR  KPOINTS  POSCAR 和OSZICAR 文件，能幫助檢查一下嗎，怎么改可以收斂嗎？

謝謝了！
INCAR
SYSTEM   = M2O3
ENCUT    = 300
ISTART   = 0
ICHARG   = 2
ISMEAR   = 1
PREC     = M
EDIFF    = 1E-3
ISIF     = 3
IBRION   = 1
LREAL    = Auto
IALGO    = 48
POTIM    = 0.15



KPOINTS
Automatic generation
0
Monkhorst-Pack
2  2  2
0  0  0

POSCAR
M2O3                           
   10.5500000000000     
     1.0000000000000000    0.0000000000000000    0.0000000000000000
     0.0000000000000000    1.0000000000000000    0.0000000000000000
     0.0000000000000000    0.0000000000000000    1.0000000000000000
  32  48
Direct
  0.2500000000000000  0.0000000000000000  0.0000000000000000
  0.0000000000000000  0.2500000000000000  0.0000000000000000
  0.5000000000000000  0.2500000000000000  0.0000000000000000
  0.2500000000000000  0.5000000000000000  0.0000000000000000
  0.0000000000000000  0.0000000000000000  0.2500000000000000
  0.5000000000000000  0.0000000000000000  0.2500000000000000
  0.0000000000000000  0.5000000000000000  0.2500000000000000
  0.2500000000000000  0.2500000000000000  0.2500000000000000
  0.5000000000000000  0.5000000000000000  0.2500000000000000
  0.2500000000000000  0.0000000000000000  0.5000000000000000
  0.0000000000000000  0.2500000000000000  0.5000000000000000
  0.5000000000000000  0.2500000000000000  0.5000000000000000
  0.2500000000000000  0.5000000000000000  0.5000000000000000
  0.7500000000000000  0.0000000000000000  0.0000000000000000
  0.0000000000000000  0.7500000000000000  0.0000000000000000
  0.5000000000000000  0.7500000000000000  0.0000000000000000
  0.7500000000000000  0.5000000000000000  0.0000000000000000
  0.0000000000000000  0.0000000000000000  0.7500000000000000
  0.5000000000000000  0.0000000000000000  0.7500000000000000
  0.0000000000000000  0.5000000000000000  0.7500000000000000
  0.7500000000000000  0.7500000000000000  0.7500000000000000
  0.5000000000000000  0.5000000000000000  0.7500000000000000
  0.7500000000000000  0.0000000000000000  0.5000000000000000
  0.0000000000000000  0.7500000000000000  0.5000000000000000
  0.5000000000000000  0.7500000000000000  0.5000000000000000
  0.7500000000000000  0.5000000000000000  0.5000000000000000
  0.2500000000000000  0.7500000000000000  0.2500000000000000
  0.7500000000000000  0.2500000000000000  0.7500000000000000
  0.7500000000000000  0.7500000000000000  0.2500000000000000
  0.2500000000000000  0.2500000000000000  0.7500000000000000
  0.7500000000000000  0.2500000000000000  0.2500000000000000
  0.2500000000000000  0.7500000000000000  0.7500000000000000
  0.1250000000000000  0.1250000000000000  0.1250000000000000
  0.3750000000000000  0.1250000000000000  0.1250000000000000
  0.1250000000000000  0.3750000000000000  0.1250000000000000
  0.3750000000000000  0.1250000000000000  0.3750000000000000
  0.1250000000000000  0.3750000000000000  0.3750000000000000
  0.3750000000000000  0.3750000000000000  0.3750000000000000
  0.8750000000000000  0.8750000000000000  0.8750000000000000
  0.6250000000000000  0.8750000000000000  0.8750000000000000
  0.8750000000000000  0.6250000000000000  0.8750000000000000
  0.6250000000000000  0.8750000000000000  0.6250000000000000
  0.8750000000000000  0.6250000000000000  0.6250000000000000
  0.6250000000000000  0.6250000000000000  0.6250000000000000
  0.1250000000000000  0.8750000000000000  0.1250000000000000
  0.3750000000000000  0.8750000000000000  0.1250000000000000
  0.3750000000000000  0.6250000000000000  0.1250000000000000
  0.1250000000000000  0.8750000000000000  0.3750000000000000
  0.1250000000000000  0.6250000000000000  0.3750000000000000
  0.3750000000000000  0.6250000000000000  0.3750000000000000
  0.8750000000000000  0.1250000000000000  0.8750000000000000
  0.6250000000000000  0.1250000000000000  0.8750000000000000
  0.6250000000000000  0.3750000000000000  0.8750000000000000
  0.8750000000000000  0.1250000000000000  0.6250000000000000
  0.8750000000000000  0.3750000000000000  0.6250000000000000
  0.6250000000000000  0.3750000000000000  0.6250000000000000
  0.8750000000000000  0.8750000000000000  0.1250000000000000
  0.8750000000000000  0.6250000000000000  0.1250000000000000
  0.6250000000000000  0.6250000000000000  0.1250000000000000
  0.8750000000000000  0.8750000000000000  0.3750000000000000
  0.6250000000000000  0.8750000000000000  0.3750000000000000
  0.6250000000000000  0.6250000000000000  0.3750000000000000
  0.1250000000000000  0.1250000000000000  0.8750000000000000
  0.1250000000000000  0.3750000000000000  0.8750000000000000
  0.3750000000000000  0.3750000000000000  0.8750000000000000
  0.1250000000000000  0.1250000000000000  0.6250000000000000
  0.3750000000000000  0.1250000000000000  0.6250000000000000
  0.3750000000000000  0.3750000000000000  0.6250000000000000
  0.6250000000000000  0.1250000000000000  0.1250000000000000
  0.8750000000000000  0.3750000000000000  0.1250000000000000
  0.6250000000000000  0.3750000000000000  0.1250000000000000
  0.8750000000000000  0.1250000000000000  0.3750000000000000
  0.6250000000000000  0.1250000000000000  0.3750000000000000
  0.8750000000000000  0.3750000000000000  0.3750000000000000
  0.3750000000000000  0.8750000000000000  0.8750000000000000
  0.1250000000000000  0.6250000000000000  0.8750000000000000
  0.3750000000000000  0.6250000000000000  0.8750000000000000
  0.1250000000000000  0.8750000000000000  0.6250000000000000
  0.3750000000000000  0.8750000000000000  0.6250000000000000
  0.1250000000000000  0.6250000000000000  0.6250000000000000
追求真理
17樓2012-02-27 00:22:56
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

bjmaowei25

銅蟲(chóng) (小有名氣)
引用回帖:
樓: Originally posted by chgzhong at 2012-02-25 18:16:00:
默認(rèn)的電子步是60步，在你運(yùn)行60步仍然沒(méi)有收斂的情況下，自動(dòng)停下。此時(shí)離子位置并沒(méi)有動(dòng)過(guò)，當(dāng)然COUNTCAR中沒(méi)有內(nèi)容。
再看看你原子坐標(biāo)是否有錯(cuò)。

謝謝你的回復(fù)，我修改了 K ISIF IBRION  POTIM等參數(shù)。
很幸運(yùn)的得到了 CONTCAR 坐標(biāo)和POSCAR一樣。但是運(yùn)行到60步后能量不收斂！
檢查了POSCAR 沒(méi)有問(wèn)題啊，因?yàn)槲矣檬謱?xiě)的坐標(biāo)輸入VESTA能打開(kāi)正確的晶體結(jié)構(gòu)。
我貼出INCAR  KPOINTS  POSCAR 和OSZICAR 文件，能幫助檢查一下嗎，怎么改可以收斂嗎？

謝謝了！

INCAR
SYSTEM   = M2O3
ENCUT    = 300
ISTART   = 0
ICHARG   = 2
ISMEAR   = 1
PREC     = M
EDIFF    = 1E-3
ISIF     = 3
IBRION   = 1
LREAL    = Auto
IALGO    = 48
POTIM    = 0.15



KPOINTS
Automatic generation
0
Monkhorst-Pack
2  2  2
0  0  0

POSCAR
M2O3                           
   10.5500000000000     
     1.0000000000000000    0.0000000000000000    0.0000000000000000
     0.0000000000000000    1.0000000000000000    0.0000000000000000
     0.0000000000000000    0.0000000000000000    1.0000000000000000
  32  48
Direct
  0.2500000000000000  0.0000000000000000  0.0000000000000000
  0.0000000000000000  0.2500000000000000  0.0000000000000000
  0.5000000000000000  0.2500000000000000  0.0000000000000000
  0.2500000000000000  0.5000000000000000  0.0000000000000000
  0.0000000000000000  0.0000000000000000  0.2500000000000000
  0.5000000000000000  0.0000000000000000  0.2500000000000000
  0.0000000000000000  0.5000000000000000  0.2500000000000000
  0.2500000000000000  0.2500000000000000  0.2500000000000000
  0.5000000000000000  0.5000000000000000  0.2500000000000000
  0.2500000000000000  0.0000000000000000  0.5000000000000000
  0.0000000000000000  0.2500000000000000  0.5000000000000000
  0.5000000000000000  0.2500000000000000  0.5000000000000000
  0.2500000000000000  0.5000000000000000  0.5000000000000000
  0.7500000000000000  0.0000000000000000  0.0000000000000000
  0.0000000000000000  0.7500000000000000  0.0000000000000000
  0.5000000000000000  0.7500000000000000  0.0000000000000000
  0.7500000000000000  0.5000000000000000  0.0000000000000000
  0.0000000000000000  0.0000000000000000  0.7500000000000000
  0.5000000000000000  0.0000000000000000  0.7500000000000000
  0.0000000000000000  0.5000000000000000  0.7500000000000000
  0.7500000000000000  0.7500000000000000  0.7500000000000000
  0.5000000000000000  0.5000000000000000  0.7500000000000000
  0.7500000000000000  0.0000000000000000  0.5000000000000000
  0.0000000000000000  0.7500000000000000  0.5000000000000000
  0.5000000000000000  0.7500000000000000  0.5000000000000000
  0.7500000000000000  0.5000000000000000  0.5000000000000000
  0.2500000000000000  0.7500000000000000  0.2500000000000000
  0.7500000000000000  0.2500000000000000  0.7500000000000000
  0.7500000000000000  0.7500000000000000  0.2500000000000000
  0.2500000000000000  0.2500000000000000  0.7500000000000000
  0.7500000000000000  0.2500000000000000  0.2500000000000000
  0.2500000000000000  0.7500000000000000  0.7500000000000000
  0.1250000000000000  0.1250000000000000  0.1250000000000000
  0.3750000000000000  0.1250000000000000  0.1250000000000000
  0.1250000000000000  0.3750000000000000  0.1250000000000000
  0.3750000000000000  0.1250000000000000  0.3750000000000000
  0.1250000000000000  0.3750000000000000  0.3750000000000000
  0.3750000000000000  0.3750000000000000  0.3750000000000000
  0.8750000000000000  0.8750000000000000  0.8750000000000000
  0.6250000000000000  0.8750000000000000  0.8750000000000000
  0.8750000000000000  0.6250000000000000  0.8750000000000000
  0.6250000000000000  0.8750000000000000  0.6250000000000000
  0.8750000000000000  0.6250000000000000  0.6250000000000000
  0.6250000000000000  0.6250000000000000  0.6250000000000000
  0.1250000000000000  0.8750000000000000  0.1250000000000000
  0.3750000000000000  0.8750000000000000  0.1250000000000000
  0.3750000000000000  0.6250000000000000  0.1250000000000000
  0.1250000000000000  0.8750000000000000  0.3750000000000000
  0.1250000000000000  0.6250000000000000  0.3750000000000000
  0.3750000000000000  0.6250000000000000  0.3750000000000000
  0.8750000000000000  0.1250000000000000  0.8750000000000000
  0.6250000000000000  0.1250000000000000  0.8750000000000000
  0.6250000000000000  0.3750000000000000  0.8750000000000000
  0.8750000000000000  0.1250000000000000  0.6250000000000000
  0.8750000000000000  0.3750000000000000  0.6250000000000000
  0.6250000000000000  0.3750000000000000  0.6250000000000000
  0.8750000000000000  0.8750000000000000  0.1250000000000000
  0.8750000000000000  0.6250000000000000  0.1250000000000000
  0.6250000000000000  0.6250000000000000  0.1250000000000000
  0.8750000000000000  0.8750000000000000  0.3750000000000000
  0.6250000000000000  0.8750000000000000  0.3750000000000000
  0.6250000000000000  0.6250000000000000  0.3750000000000000
  0.1250000000000000  0.1250000000000000  0.8750000000000000
  0.1250000000000000  0.3750000000000000  0.8750000000000000
  0.3750000000000000  0.3750000000000000  0.8750000000000000
  0.1250000000000000  0.1250000000000000  0.6250000000000000
  0.3750000000000000  0.1250000000000000  0.6250000000000000
  0.3750000000000000  0.3750000000000000  0.6250000000000000
  0.6250000000000000  0.1250000000000000  0.1250000000000000
  0.8750000000000000  0.3750000000000000  0.1250000000000000
  0.6250000000000000  0.3750000000000000  0.1250000000000000
  0.8750000000000000  0.1250000000000000  0.3750000000000000
  0.6250000000000000  0.1250000000000000  0.3750000000000000
  0.8750000000000000  0.3750000000000000  0.3750000000000000
  0.3750000000000000  0.8750000000000000  0.8750000000000000
  0.1250000000000000  0.6250000000000000  0.8750000000000000
  0.3750000000000000  0.6250000000000000  0.8750000000000000
  0.1250000000000000  0.8750000000000000  0.6250000000000000
  0.3750000000000000  0.8750000000000000  0.6250000000000000
  0.1250000000000000  0.6250000000000000  0.6250000000000000
追求真理
18樓2012-02-27 00:26:14
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

bjmaowei25

銅蟲(chóng) (小有名氣)
引用回帖:
樓: Originally posted by liliangfang at 2012-02-25 19:22:22:
你檢查一下你的POSCAR看看結(jié)構(gòu)是否合理

謝謝版主的回復(fù)，我修改 K meshISIF IBRION  POTIM等參數(shù)。
很幸運(yùn)的得到了 CONTCAR 坐標(biāo)和POSCAR一樣。但是運(yùn)行到60步后能量不收斂！
檢查了POSCAR 沒(méi)有問(wèn)題啊，因?yàn)槲矣檬謱?xiě)的坐標(biāo)輸入VESTA能打開(kāi)正確的晶體結(jié)構(gòu)。
我貼出INCAR  KPOINTS  POSCAR 和OSZICAR 文件，能幫助檢查POSCAR一下嗎,看是否正確，
怎樣改其它參數(shù)才可以收斂呢？

謝謝了！

INCAR
SYSTEM   = M2O3
ENCUT    = 300
ISTART   = 0
ICHARG   = 2
ISMEAR   = 1
PREC     = M
EDIFF    = 1E-3
ISIF     = 3
IBRION   = 1
LREAL    = Auto
IALGO    = 48
POTIM    = 0.15



KPOINTS
Automatic generation
0
Monkhorst-Pack
2  2  2
0  0  0

POSCAR
M2O3                           
   10.5500000000000     
     1.0000000000000000    0.0000000000000000    0.0000000000000000
     0.0000000000000000    1.0000000000000000    0.0000000000000000
     0.0000000000000000    0.0000000000000000    1.0000000000000000
  32  48
Direct
  0.2500000000000000  0.0000000000000000  0.0000000000000000
  0.0000000000000000  0.2500000000000000  0.0000000000000000
  0.5000000000000000  0.2500000000000000  0.0000000000000000
  0.2500000000000000  0.5000000000000000  0.0000000000000000
  0.0000000000000000  0.0000000000000000  0.2500000000000000
  0.5000000000000000  0.0000000000000000  0.2500000000000000
  0.0000000000000000  0.5000000000000000  0.2500000000000000
  0.2500000000000000  0.2500000000000000  0.2500000000000000
  0.5000000000000000  0.5000000000000000  0.2500000000000000
  0.2500000000000000  0.0000000000000000  0.5000000000000000
  0.0000000000000000  0.2500000000000000  0.5000000000000000
  0.5000000000000000  0.2500000000000000  0.5000000000000000
  0.2500000000000000  0.5000000000000000  0.5000000000000000
  0.7500000000000000  0.0000000000000000  0.0000000000000000
  0.0000000000000000  0.7500000000000000  0.0000000000000000
  0.5000000000000000  0.7500000000000000  0.0000000000000000
  0.7500000000000000  0.5000000000000000  0.0000000000000000
  0.0000000000000000  0.0000000000000000  0.7500000000000000
  0.5000000000000000  0.0000000000000000  0.7500000000000000
  0.0000000000000000  0.5000000000000000  0.7500000000000000
  0.7500000000000000  0.7500000000000000  0.7500000000000000
  0.5000000000000000  0.5000000000000000  0.7500000000000000
  0.7500000000000000  0.0000000000000000  0.5000000000000000
  0.0000000000000000  0.7500000000000000  0.5000000000000000
  0.5000000000000000  0.7500000000000000  0.5000000000000000
  0.7500000000000000  0.5000000000000000  0.5000000000000000
  0.2500000000000000  0.7500000000000000  0.2500000000000000
  0.7500000000000000  0.2500000000000000  0.7500000000000000
  0.7500000000000000  0.7500000000000000  0.2500000000000000
  0.2500000000000000  0.2500000000000000  0.7500000000000000
  0.7500000000000000  0.2500000000000000  0.2500000000000000
  0.2500000000000000  0.7500000000000000  0.7500000000000000
  0.1250000000000000  0.1250000000000000  0.1250000000000000
  0.3750000000000000  0.1250000000000000  0.1250000000000000
  0.1250000000000000  0.3750000000000000  0.1250000000000000
  0.3750000000000000  0.1250000000000000  0.3750000000000000
  0.1250000000000000  0.3750000000000000  0.3750000000000000
  0.3750000000000000  0.3750000000000000  0.3750000000000000
  0.8750000000000000  0.8750000000000000  0.8750000000000000
  0.6250000000000000  0.8750000000000000  0.8750000000000000
  0.8750000000000000  0.6250000000000000  0.8750000000000000
  0.6250000000000000  0.8750000000000000  0.6250000000000000
  0.8750000000000000  0.6250000000000000  0.6250000000000000
  0.6250000000000000  0.6250000000000000  0.6250000000000000
  0.1250000000000000  0.8750000000000000  0.1250000000000000
  0.3750000000000000  0.8750000000000000  0.1250000000000000
  0.3750000000000000  0.6250000000000000  0.1250000000000000
  0.1250000000000000  0.8750000000000000  0.3750000000000000
  0.1250000000000000  0.6250000000000000  0.3750000000000000
  0.3750000000000000  0.6250000000000000  0.3750000000000000
  0.8750000000000000  0.1250000000000000  0.8750000000000000
  0.6250000000000000  0.1250000000000000  0.8750000000000000
  0.6250000000000000  0.3750000000000000  0.8750000000000000
  0.8750000000000000  0.1250000000000000  0.6250000000000000
  0.8750000000000000  0.3750000000000000  0.6250000000000000
  0.6250000000000000  0.3750000000000000  0.6250000000000000
  0.8750000000000000  0.8750000000000000  0.1250000000000000
  0.8750000000000000  0.6250000000000000  0.1250000000000000
  0.6250000000000000  0.6250000000000000  0.1250000000000000
  0.8750000000000000  0.8750000000000000  0.3750000000000000
  0.6250000000000000  0.8750000000000000  0.3750000000000000
  0.6250000000000000  0.6250000000000000  0.3750000000000000
  0.1250000000000000  0.1250000000000000  0.8750000000000000
  0.1250000000000000  0.3750000000000000  0.8750000000000000
  0.3750000000000000  0.3750000000000000  0.8750000000000000
  0.1250000000000000  0.1250000000000000  0.6250000000000000
  0.3750000000000000  0.1250000000000000  0.6250000000000000
  0.3750000000000000  0.3750000000000000  0.6250000000000000
  0.6250000000000000  0.1250000000000000  0.1250000000000000
  0.8750000000000000  0.3750000000000000  0.1250000000000000
  0.6250000000000000  0.3750000000000000  0.1250000000000000
  0.8750000000000000  0.1250000000000000  0.3750000000000000
  0.6250000000000000  0.1250000000000000  0.3750000000000000
  0.8750000000000000  0.3750000000000000  0.3750000000000000
  0.3750000000000000  0.8750000000000000  0.8750000000000000
  0.1250000000000000  0.6250000000000000  0.8750000000000000
  0.3750000000000000  0.6250000000000000  0.8750000000000000
  0.1250000000000000  0.8750000000000000  0.6250000000000000
  0.3750000000000000  0.8750000000000000  0.6250000000000000
  0.1250000000000000  0.6250000000000000  0.6250000000000000
追求真理
19樓2012-02-27 00:28:48
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

bjmaowei25

銅蟲(chóng) (小有名氣)
引用回帖:
樓: Originally posted by guodonglin at 2012-02-26 19:29:31:
你的CONTCAR時(shí)空的，是正確的。你使用了默認(rèn)設(shè)置NELM=60，所以跑了60步就停了。但是你的EDIFF= 1E-3，可是看你的輸出結(jié)果，還沒(méi)有達(dá)到電子步的收斂精度，你還沒(méi)有完成一個(gè)離子步，也沒(méi)有提示你優(yōu)化成功了，所以CO ...

謝謝你的回復(fù)，我修改了 K mesh ISIF IBRION  POTIM等參數(shù)。
很幸運(yùn)的得到了 CONTCAR 坐標(biāo)和POSCAR一樣。但是運(yùn)行到60步后能量還不不收斂！
檢查了POSCAR 沒(méi)有問(wèn)題啊，因?yàn)槲矣檬謱?xiě)的坐標(biāo)輸入VESTA能打開(kāi)正確的晶體結(jié)構(gòu)。
我貼出INCAR  KPOINTS  POSCAR 和OSZICAR 文件，能幫助檢查一下嗎，怎么改可以收斂嗎？

謝謝了！

INCAR
SYSTEM   = M2O3
ENCUT    = 300
ISTART   = 0
ICHARG   = 2
ISMEAR   = 1
PREC     = M
EDIFF    = 1E-3
ISIF     = 3
IBRION   = 1
LREAL    = Auto
IALGO    = 48
POTIM    = 0.15



KPOINTS
Automatic generation
0
Monkhorst-Pack
2  2  2
0  0  0

POSCAR
M2O3                           
   10.5500000000000     
     1.0000000000000000    0.0000000000000000    0.0000000000000000
     0.0000000000000000    1.0000000000000000    0.0000000000000000
     0.0000000000000000    0.0000000000000000    1.0000000000000000
  32  48
Direct
  0.2500000000000000  0.0000000000000000  0.0000000000000000
  0.0000000000000000  0.2500000000000000  0.0000000000000000
  0.5000000000000000  0.2500000000000000  0.0000000000000000
  0.2500000000000000  0.5000000000000000  0.0000000000000000
  0.0000000000000000  0.0000000000000000  0.2500000000000000
  0.5000000000000000  0.0000000000000000  0.2500000000000000
  0.0000000000000000  0.5000000000000000  0.2500000000000000
  0.2500000000000000  0.2500000000000000  0.2500000000000000
  0.5000000000000000  0.5000000000000000  0.2500000000000000
  0.2500000000000000  0.0000000000000000  0.5000000000000000
  0.0000000000000000  0.2500000000000000  0.5000000000000000
  0.5000000000000000  0.2500000000000000  0.5000000000000000
  0.2500000000000000  0.5000000000000000  0.5000000000000000
  0.7500000000000000  0.0000000000000000  0.0000000000000000
  0.0000000000000000  0.7500000000000000  0.0000000000000000
  0.5000000000000000  0.7500000000000000  0.0000000000000000
  0.7500000000000000  0.5000000000000000  0.0000000000000000
  0.0000000000000000  0.0000000000000000  0.7500000000000000
  0.5000000000000000  0.0000000000000000  0.7500000000000000
  0.0000000000000000  0.5000000000000000  0.7500000000000000
  0.7500000000000000  0.7500000000000000  0.7500000000000000
  0.5000000000000000  0.5000000000000000  0.7500000000000000
  0.7500000000000000  0.0000000000000000  0.5000000000000000
  0.0000000000000000  0.7500000000000000  0.5000000000000000
  0.5000000000000000  0.7500000000000000  0.5000000000000000
  0.7500000000000000  0.5000000000000000  0.5000000000000000
  0.2500000000000000  0.7500000000000000  0.2500000000000000
  0.7500000000000000  0.2500000000000000  0.7500000000000000
  0.7500000000000000  0.7500000000000000  0.2500000000000000
  0.2500000000000000  0.2500000000000000  0.7500000000000000
  0.7500000000000000  0.2500000000000000  0.2500000000000000
  0.2500000000000000  0.7500000000000000  0.7500000000000000
  0.1250000000000000  0.1250000000000000  0.1250000000000000
  0.3750000000000000  0.1250000000000000  0.1250000000000000
  0.1250000000000000  0.3750000000000000  0.1250000000000000
  0.3750000000000000  0.1250000000000000  0.3750000000000000
  0.1250000000000000  0.3750000000000000  0.3750000000000000
  0.3750000000000000  0.3750000000000000  0.3750000000000000
  0.8750000000000000  0.8750000000000000  0.8750000000000000
  0.6250000000000000  0.8750000000000000  0.8750000000000000
  0.8750000000000000  0.6250000000000000  0.8750000000000000
  0.6250000000000000  0.8750000000000000  0.6250000000000000
  0.8750000000000000  0.6250000000000000  0.6250000000000000
  0.6250000000000000  0.6250000000000000  0.6250000000000000
  0.1250000000000000  0.8750000000000000  0.1250000000000000
  0.3750000000000000  0.8750000000000000  0.1250000000000000
  0.3750000000000000  0.6250000000000000  0.1250000000000000
  0.1250000000000000  0.8750000000000000  0.3750000000000000
  0.1250000000000000  0.6250000000000000  0.3750000000000000
  0.3750000000000000  0.6250000000000000  0.3750000000000000
  0.8750000000000000  0.1250000000000000  0.8750000000000000
  0.6250000000000000  0.1250000000000000  0.8750000000000000
  0.6250000000000000  0.3750000000000000  0.8750000000000000
  0.8750000000000000  0.1250000000000000  0.6250000000000000
  0.8750000000000000  0.3750000000000000  0.6250000000000000
  0.6250000000000000  0.3750000000000000  0.6250000000000000
  0.8750000000000000  0.8750000000000000  0.1250000000000000
  0.8750000000000000  0.6250000000000000  0.1250000000000000
  0.6250000000000000  0.6250000000000000  0.1250000000000000
  0.8750000000000000  0.8750000000000000  0.3750000000000000
  0.6250000000000000  0.8750000000000000  0.3750000000000000
  0.6250000000000000  0.6250000000000000  0.3750000000000000
  0.1250000000000000  0.1250000000000000  0.8750000000000000
  0.1250000000000000  0.3750000000000000  0.8750000000000000
  0.3750000000000000  0.3750000000000000  0.8750000000000000
  0.1250000000000000  0.1250000000000000  0.6250000000000000
  0.3750000000000000  0.1250000000000000  0.6250000000000000
  0.3750000000000000  0.3750000000000000  0.6250000000000000
  0.6250000000000000  0.1250000000000000  0.1250000000000000
  0.8750000000000000  0.3750000000000000  0.1250000000000000
  0.6250000000000000  0.3750000000000000  0.1250000000000000
  0.8750000000000000  0.1250000000000000  0.3750000000000000
  0.6250000000000000  0.1250000000000000  0.3750000000000000
  0.8750000000000000  0.3750000000000000  0.3750000000000000
  0.3750000000000000  0.8750000000000000  0.8750000000000000
  0.1250000000000000  0.6250000000000000  0.8750000000000000
  0.3750000000000000  0.6250000000000000  0.8750000000000000
  0.1250000000000000  0.8750000000000000  0.6250000000000000
  0.3750000000000000  0.8750000000000000  0.6250000000000000
  0.1250000000000000  0.6250000000000000  0.6250000000000000
追求真理
20樓2012-02-27 00:31:12
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

bjmaowei25

銅蟲(chóng) (小有名氣)
引用回帖:
樓: Originally posted by liliangfang at 2012-02-27 08:43:25:
降低收斂精度，多跑幾步看看

謝謝版主的回復(fù)，你的 意識(shí)是將EDIFF 再增大嗎，現(xiàn)在為1x10-3 改為1x10-2 or 1x10-1?

另外有必要加入EDIFFG  參數(shù)嗎，離子做離子弛豫，初期設(shè)置為-0.03～-0.01？
追求真理
22樓2012-02-27 13:48:30
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

bjmaowei25

銅蟲(chóng) (小有名氣)
引用回帖:
樓: Originally posted by liliangfang at 2012-02-25 10:55:00:
1樓說(shuō)的很明白了，建議多看看說(shuō)明書(shū)，有中文的，比如侯柱峰老師講的就很好

版主有個(gè)小問(wèn)題，請(qǐng)教你一下。
我把系統(tǒng)改為10個(gè)原子�？纯此氖諗啃圆�(yōu)化。
設(shè)置運(yùn)行的NSW=1000, 可運(yùn)行60步后自動(dòng)停止，并告訴精度到達(dá)要求，但明顯dE 很大沒(méi)收斂，求點(diǎn)評(píng)，謝謝！


INCAR

SYSTEM   = M2O3
ENCUT    = 500
ISTART   = 0
ICHARG   = 2
ISMEAR   = 1
PREC     = M
EDIFF    = 1E-3
NSW      = 1000
ISIF     = 2
IBRION   = 2
IALGO    = 48
ALGO     = very fast
POTIM    = 0.1
LREAL    = Auto


OSZICAR
N       E                     dE             d eps       ncg     rms          rms(c)
RMM:   1     0.923112380146E+03    0.92311E+03   -0.50333E+04    74   0.357E+03
RMM:   2     0.175369114913E+03   -0.74774E+03   -0.10738E+04    74   0.394E+02
RMM:   3     0.174267355228E+03   -0.11018E+01   -0.15431E+03    74   0.342E+02
RMM:   4    -0.542524666320E+01   -0.17969E+03   -0.86607E+02    74   0.123E+02
RMM:   5    -0.373367180864E+02   -0.31911E+02   -0.29624E+02    74   0.151E+02
RMM:   6    -0.677436920725E+02   -0.30407E+02   -0.17804E+02    74   0.658E+01
RMM:   7    -0.800862495786E+02   -0.12343E+02   -0.74106E+01    74   0.524E+01
RMM:   8    -0.882359395166E+02   -0.81497E+01   -0.51421E+01    74   0.285E+01
RMM:   9    -0.955797020611E+02   -0.73438E+01   -0.58969E+01   248   0.225E+01
RMM:  10    -0.963124520811E+02   -0.73275E+00   -0.92937E+00   283   0.642E+00
RMM:  11    -0.964586535670E+02   -0.14620E+00   -0.14783E+00   277   0.327E+00
RMM:  12    -0.965482989444E+02   -0.89645E-01   -0.90834E-01   280   0.220E+00    0.216E+01
RMM:  13    -0.990880302596E+02   -0.25397E+01   -0.72834E+01   212   0.238E+01    0.466E+01
RMM:  14    -0.988030396415E+02    0.28499E+00   -0.71966E+01   192   0.442E+01    0.294E+01
RMM:  15    -0.901814394107E+02    0.86216E+01   -0.84060E+00   205   0.176E+01    0.184E+01
RMM:  16    -0.885813824920E+02    0.16001E+01   -0.98976E-01   235   0.470E+00    0.134E+01
RMM:  17    -0.906172378817E+02   -0.20359E+01   -0.11675E+00   233   0.458E+00    0.146E+01
RMM:  18    -0.916309571854E+02   -0.10137E+01   -0.14801E+00   211   0.617E+00    0.991E+00
RMM:  19    -0.889660004866E+02    0.26650E+01   -0.57077E+00   192   0.128E+01    0.581E+00
RMM:  20    -0.887204715274E+02    0.24553E+00   -0.11635E+00   211   0.564E+00    0.551E+00
RMM:  21    -0.864728355587E+02    0.22476E+01   -0.50053E-01   194   0.416E+00    0.479E+00
RMM:  22    -0.881407375740E+02   -0.16679E+01   -0.49187E-01   199   0.354E+00    0.315E+00
RMM:  23    -0.870375191558E+02    0.11032E+01   -0.34042E-01   159   0.352E+00    0.548E+00
RMM:  24    -0.900241807071E+02   -0.29867E+01   -0.71858E-01   153   0.443E+00    0.335E+00
RMM:  25    -0.871196437618E+02    0.29045E+01   -0.48867E-01   173   0.412E+00    0.515E+00
RMM:  26    -0.904592699743E+02   -0.33396E+01   -0.74835E-01   160   0.450E+00    0.380E+00
RMM:  27    -0.872593727356E+02    0.31999E+01   -0.52872E-01   159   0.449E+00    0.491E+00
RMM:  28    -0.908610841605E+02   -0.36017E+01   -0.67250E-01   155   0.461E+00    0.423E+00
RMM:  29    -0.876984696886E+02    0.31626E+01   -0.63891E-01   152   0.492E+00    0.464E+00
RMM:  30    -0.909369033777E+02   -0.32384E+01   -0.69348E-01   147   0.468E+00    0.460E+00
RMM:  31    -0.881459421380E+02    0.27910E+01   -0.85276E-01   150   0.536E+00    0.423E+00
RMM:  32    -0.906263470891E+02   -0.24804E+01   -0.68250E-01   148   0.472E+00    0.483E+00
RMM:  33    -0.884082531914E+02    0.22181E+01   -0.88297E-01   149   0.547E+00    0.392E+00
RMM:  34    -0.903375932469E+02   -0.19293E+01   -0.66028E-01   145   0.466E+00    0.494E+00
RMM:  35    -0.885481926372E+02    0.17894E+01   -0.89977E-01   148   0.547E+00    0.373E+00
RMM:  36    -0.901607324062E+02   -0.16125E+01   -0.67139E-01   145   0.464E+00    0.500E+00
RMM:  37    -0.886275118245E+02    0.15332E+01   -0.94080E-01   148   0.547E+00    0.362E+00
RMM:  38    -0.900544748007E+02   -0.14270E+01   -0.69192E-01   145   0.464E+00    0.503E+00
RMM:  39    -0.886695168003E+02    0.13850E+01   -0.94613E-01   149   0.549E+00    0.355E+00
RMM:  40    -0.899878387757E+02   -0.13183E+01   -0.68633E-01   146   0.463E+00    0.505E+00
RMM:  41    -0.886901460457E+02    0.12977E+01   -0.94424E-01   145   0.547E+00    0.350E+00
RMM:  42    -0.899453742520E+02   -0.12552E+01   -0.68720E-01   145   0.464E+00    0.507E+00
RMM:  43    -0.887034990119E+02    0.12419E+01   -0.95275E-01   145   0.545E+00    0.346E+00
RMM:  44    -0.899179872721E+02   -0.12145E+01   -0.69601E-01   147   0.464E+00    0.510E+00
RMM:  45    -0.887090983630E+02    0.12089E+01   -0.94088E-01   142   0.547E+00    0.342E+00
RMM:  46    -0.898958996580E+02   -0.11868E+01   -0.68101E-01   147   0.462E+00    0.512E+00
RMM:  47    -0.887119752044E+02    0.11839E+01   -0.93817E-01   142   0.543E+00    0.337E+00
RMM:  48    -0.898810921160E+02   -0.11691E+01   -0.68506E-01   149   0.464E+00    0.515E+00
RMM:  49    -0.887149032964E+02    0.11662E+01   -0.94807E-01   141   0.545E+00    0.333E+00
RMM:  50    -0.898688771465E+02   -0.11540E+01   -0.69777E-01   146   0.462E+00    0.518E+00
RMM:  51    -0.887154505611E+02    0.11534E+01   -0.94228E-01   141   0.544E+00    0.329E+00
RMM:  52    -0.898563194566E+02   -0.11409E+01   -0.68267E-01   151   0.464E+00    0.521E+00
RMM:  53    -0.887138792528E+02    0.11424E+01   -0.93836E-01   140   0.544E+00    0.324E+00
RMM:  54    -0.898450068695E+02   -0.11311E+01   -0.68500E-01   146   0.462E+00    0.524E+00
RMM:  55    -0.887150853364E+02    0.11299E+01   -0.95461E-01   140   0.543E+00    0.319E+00
RMM:  56    -0.898350463963E+02   -0.11200E+01   -0.69105E-01   148   0.462E+00    0.527E+00
RMM:  57    -0.887128917073E+02    0.11222E+01   -0.94546E-01   141   0.545E+00    0.314E+00
RMM:  58    -0.898213736552E+02   -0.11085E+01   -0.66544E-01   147   0.462E+00    0.529E+00
RMM:  59    -0.887064903571E+02    0.11149E+01   -0.93655E-01   142   0.542E+00    0.312E+00
RMM:  60    -0.898133215169E+02   -0.11068E+01   -0.66742E-01   148   0.459E+00
   1 F= -.89813322E+02 E0= -.89816821E+02  d E =-.898133E+02
curvature:   0.00 expect dE= 0.000E+00 dE for cont linesearch  0.000E+00
trial: gam= 0.00000 g(F)=  0.144E-58 g(S)=  0.000E+00 ort = 0.000E+00 (trialstep = 0.100E+01)
search vector abs. value=  0.100E-09
reached required accuracy - stopping structural energy minimisation
writing wavefunctions
追求真理
25樓2012-02-28 19:21:27
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖
回帖支持 ( 顯示支持度最高的前 50 名 )

cavediger

金蟲(chóng) (正式寫(xiě)手)
liliangfang(金幣+2): 多謝指教 2012-02-25 10:47:56
1.收斂精度要看你的要求和計(jì)算的性質(zhì)，不知道時(shí)就先用默認(rèn)值
2.優(yōu)化收斂與否就是看是否達(dá)到了你的要求，一般可以取少數(shù)k點(diǎn)，待收斂后，增加k點(diǎn)，cp CONTCAR POSCAR，再進(jìn)一步優(yōu)化，此時(shí)不需要?jiǎng)h除波函數(shù)文件。
3.初始結(jié)構(gòu)盡量給準(zhǔn)確點(diǎn)兒，這樣會(huì)節(jié)省不少時(shí)間的
搞不明白時(shí)，嘗試一下就懂了
我們都老得太快，卻聰明得太遲。當(dāng)我一路跌跌撞撞登上山巔，驀然發(fā)現(xiàn)，本來(lái)就有一條筆直的路，直達(dá)山頂....
2樓2012-02-25 10:38:38
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖
普通回帖

liliangfang

榮譽(yù)版主 (著名寫(xiě)手)
uuv2010(金幣+1): 多謝提示 2012-02-25 17:23:42
1樓說(shuō)的很明白了，建議多看看說(shuō)明書(shū)，有中文的，比如侯柱峰老師講的就很好
3樓2012-02-25 10:55:00
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

liliangfang

榮譽(yù)版主 (著名寫(xiě)手)
uuv2010(金幣+1): 多謝提示 2012-02-25 17:23:53
初步建粒子位置的時(shí)候盡可能的準(zhǔn)確些，然后再優(yōu)化，ISIF的值也要設(shè)置，論壇里面有別人寫(xiě)好的優(yōu)化步驟，你可以參考，另外80個(gè)粒子很多了，建議先優(yōu)化小一點(diǎn)的結(jié)構(gòu)
6樓2012-02-25 14:04:35
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

cavediger

金蟲(chóng) (正式寫(xiě)手)
liliangfang(金幣+1): 多謝交流 2012-02-26 08:53:37
引用回帖:
4樓: Originally posted by bjmaowei25 at 2012-02-25 13:32:28:
謝謝你的回復(fù)。我現(xiàn)在初步設(shè)計(jì)了一個(gè)包含80個(gè)原子的系統(tǒng)。為了收斂，我降低了精度，雖然運(yùn)行了幾個(gè)小時(shí)走了60步，自動(dòng)停下，但是不知到為什么看不到CONTCAR文件是空文件呢？這樣我沒(méi)法再優(yōu)化哦？
能問(wèn)下什么原因 ...

你這不是優(yōu)化，就一單點(diǎn)，而且電子步都沒(méi)收斂，好好研究一下IBRION.
另外，檢查一下你的OUTCAR，看看是否有些原子間的距離不合理
我們都老得太快，卻聰明得太遲。當(dāng)我一路跌跌撞撞登上山巔，驀然發(fā)現(xiàn)，本來(lái)就有一條筆直的路，直達(dá)山頂....
8樓2012-02-25 14:39:41
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

cavediger

金蟲(chóng) (正式寫(xiě)手)
liliangfang(金幣+1): 多謝提示 2012-02-26 08:53:55
引用回帖:
7樓: Originally posted by bjmaowei25 at 2012-02-25 14:25:49:
謝謝版主的回復(fù)，我建立的粒子坐標(biāo)，就是一個(gè)單元晶胞里面包含的的原子數(shù)目，而且用VESTA
打開(kāi)POSCAR 文件，就可以得到單元晶胞，是不是沒(méi)有設(shè)置ISIF值，就沒(méi)有CONTCAR文件呢，
因?yàn)樯鲜龀绦蜻\(yùn)行結(jié)果 CONTCAR ...

如果你的體系就一晶胞，那就不能再小了，否則有違體相材料結(jié)構(gòu)；
POTIM再取小點(diǎn)，譬如0.15，再試一下，看看收斂情況；
把K mesh取成2x2x2，優(yōu)化，收斂后調(diào)為4x4x4，cp CONTCAR POSCAR，再接著優(yōu)化，比較得到的能量，如果兩者差異較大，可再將K mesh調(diào)為5x5x5，在比較，如果沒(méi)明顯差異就可以了，如果差異不能接受，繼續(xù)增加k點(diǎn)......
我們都老得太快，卻聰明得太遲。當(dāng)我一路跌跌撞撞登上山巔，驀然發(fā)現(xiàn)，本來(lái)就有一條筆直的路，直達(dá)山頂....
9樓2012-02-25 14:50:47
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖

liliangfang

榮譽(yù)版主 (著名寫(xiě)手)
uuv2010(金幣+1): 多謝提示 2012-02-25 17:24:01
引用回帖:
7樓: Originally posted by bjmaowei25 at 2012-02-25 14:25:49:
謝謝版主的回復(fù)，我建立的粒子坐標(biāo)，就是一個(gè)單元晶胞里面包含的的原子數(shù)目，而且用VESTA
打開(kāi)POSCAR 文件，就可以得到單元晶胞，是不是沒(méi)有設(shè)置ISIF值，就沒(méi)有CONTCAR文件呢，
因?yàn)樯鲜龀绦蜻\(yùn)行結(jié)果 CONTCAR ...

那到不是，CONTCAR沒(méi)計(jì)算完一步就會(huì)更新一次，估計(jì)沒(méi)有可能是沒(méi)有計(jì)算完
10樓2012-02-25 15:38:12
已閱   回復(fù)此樓   關(guān)注TA 給TA發(fā)消息 送TA紅花 TA的回帖
相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 bjmaowei25 的主題更新
普通表情 龍 兔 虎 貓
最具人氣熱帖推薦 [查看全部] 作者 回/看 最后發(fā)表
[考研] 268求調(diào)劑 +7 好運(yùn)連綿不絕 2026-03-12 8/400 2026-03-17 20:28 by xilongliang
[考研] 301求調(diào)劑 +5 yy要上岸呀 2026-03-17 5/250 2026-03-17 20:20 by peike
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +3 生物工程調(diào)劑 2026-03-16 8/400 2026-03-17 19:03 by Wangjingyue
[考研] 271求調(diào)劑 +12 生如夏花… 2026-03-11 14/700 2026-03-17 10:56 by lovewei0727
[考研] 285化工學(xué)碩求調(diào)劑（081700） +9 柴郡貓_ 2026-03-12 9/450 2026-03-17 10:18 by Sammy2
[論文投稿] 有沒(méi)有大佬發(fā)小論文能帶我個(gè)二作 +3 增銳漏人 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 333求調(diào)劑 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
[基金申請(qǐng)] 今年的國(guó)基金是打分制嗎？ 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京萊茵潤(rùn)色
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 努力奮斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 265求調(diào)劑 +4 威化餅07 2026-03-12 4/200 2026-03-14 17:23 by userper
[考研] 復(fù)試調(diào)劑 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 學(xué)碩285求調(diào)劑 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[考研] 一志愿中科院，化學(xué)方向，295求調(diào)劑 +4 一氧二氮 2026-03-11 4/200 2026-03-13 22:35 by JourneyLucky
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9 2026調(diào)劑俠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 求調(diào)劑 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 307求調(diào)劑 +5 超級(jí)伊昂大王 2026-03-12 5/250 2026-03-13 15:56 by 棒棒球手
[考研] 求調(diào)劑 +3 程雨杭 2026-03-12 3/150 2026-03-13 15:06 by JourneyLucky
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導(dǎo)師收留 +3 土木小天才 2026-03-12 3/150 2026-03-13 15:01 by JourneyLucky
[考研] 290求調(diào)劑 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
信息提示
請(qǐng)?zhí)钐幚硪庖?jiàn)
91九色国产在线视频| 天天操天天干天天舔天天| 亚洲国产精品一区51动漫| 伊人免费观看视频一| 中文字幕日韩首页欧美在线激情| 色999日韩偷自拍拍免费| 天天操,天天射,天天爽| 亚洲码av一区二区三区| 青青草原在线播放日韩| 熟女一区二区三区综合| 伊人综合在线视频免费观看| 久久久久久免费观看av| 99亚偷拍自图区亚洲| 亚洲一区二区精品在线播放| 最近中文字幕免费视频一| 精产国品一二三产品区别97| 日本一本午夜在线播放| 91激情四射婷婷综合| 国产高清在线观看av| 91系列视频在线播放| 新香蕉视频香蕉视频2| 老司机伊人99久久精品| 亚洲欧美综合另类最新| 丰满放荡熟妇在线播放| 亚洲色大WWW永久网站| 亚洲中文字幕最新地址| 91精品91久久久久| 爱搞视频在线观看视频91| 92麻豆一区二区三区| 男生用大肌巴操美女骚穴| 色噜噜噜噜色噜噜色合久一| 日本四十路人妻熟女| 日韩人妻中文字幕区| 欧美精品激情在线不卡| 亚洲精品国产99999| 中文字幕福利视频第四页| 91精品视频在线观看视频| 2018中文字字幕人妻| 欧美久久一区二区伊人| 黄色片黄色片黄色片黄色片黄色 | 亚洲成人三级黄色片| 五月天色婷婷狠狠爱| 国产精品网站亚洲发布| 欧美 日韩 精品 中文| 日本少妇熟女乱码一区二区| 中字幕人妻熟女人妻a62v网| xxoo福利视频导航| av丝袜免费在线观看| 黄色av网址在线播放| 天天干夜夜操夜夜骑| 在线视频国产精品欧美| 又粗又长又硬又黄又爽| 国内销魂老女人老泬| 久久人人爽人人爽人人av东京热| 外国美女舔男人坤坤| 欧美精品一区二区三区观看 | 亚洲黄色免费在线观看网站| 欧美日韩久久丝袜在线| 婷婷综合缴情亚洲五月伊人| 老牛影视在线一区二区三区| 99久久国产精品免费热| 黄色大片一级老太太操逼| 黑人侵犯人妻森泽佳奈| 91九色人妻在线播放| 免费成人av麻豆| 亚洲成人三级黄色片| 国产 亚洲 欧美 自拍| 天海翼亚洲一区在线观看| 男人资源站中文字幕| 中国特黄色性生活片| 美女黄色啊啊啊啊视频| 黑人大巨屌操美女逼| 欧美日本亚欧在线观看| 免费成人av麻豆| 国产精美视频精品视频精品| 51vv精品视频在线观看| 亚洲国产综合久久精品| 亚洲国产精品一区51动漫| 免费24小时人妻视频| 一区二区三区午夜福利在线| 久99久视频免费观看中文字幕| 老司机免费视频福利0| 99精品久久精品一区二区| 中文字幕麻绳捆绑的人妻| 久久久精品人妻无码专区不卡| 中文字幕在线免费观看成人| 99在线视频精品观看高| 美女网站福利在线观看| 国产精品亚洲精品亚洲| 天天操天天干天天谢| 久久人人爽人人爽人人av东京热 | 特级aaaaa黄色片| 中文字幕亚洲乱码精品无限| 日本少妇熟女乱码一区二区| 亚洲永远av在线播放| 久久亚洲国产成人精品麻豆| 操人妻人妻天天爽天天偷| 在线视频自拍第三页| 天天干夜夜操91视频网站| 欧美成人久久久桃色aa| 黑吊操欧美极品美女| 99国产精品久久99久久久| 日韩av熟妇在线观看| 交换的一天中文字幕在线视频| 国长拍拍视频免费孕妇| 91佛爷视频在线观看| 国产福利三级在线观看| 欧美男女一区二区三区| 婷婷色九月综合激情丁香| 国产不卡免费在线观看| 911精产国品一二三产区区| 亚洲综合在线视频在线播放| 国产欧美福利在线观看| 涩涩黄片在线免费观看| 日本国产亚洲欧美色综合| —区二区三区女厕偷拍| 快使劲弄我视频在线播放| 人妻熟女 亚洲 一页二页| 51vv精品视频在线观看| 2020国产成人精品视频| 日韩成人精品久久久免费看| 久久综合狠狠综合久久综| 大乳人妻一区二区三区| 夜夜骚av一二三区| 欧美色区国产日韩亚洲区| 天堂网免费在线电影| 亚洲成人五月婷婷久久综合| 91精品一区一区三区| 精品一区二区三区免费毛片W| 日本黄页在线观看视频| 亚洲色图日韩在线视频观看| 亚洲成人中文无码在线| 一区二区三区五区六区| 久久国产精品久精国产爱| 96在线观看免费播放| 一二区二区不卡视频| 不卡一二三区别视频| 欧美日韩综合精品无人区| 午夜久久久久久av五月| 欧美男女一区二区三区| 不卡视频在线 欧美日韩| 大陆中文字幕视频在线| 农村大炕有肉大屁股熟妇| 69视频在线精品国自产拍| 国产激情免费在线视频| 午夜国产成人精品视频观看| 亚洲综合一区二区三区四区| 91青青青国产免费高清| 四虎精品久久免费最新| 国产av啊啊啊啊啊啊啊| 久久久亚洲综合国产精品| avtt中文字幕手机版| 亚洲va999天堂va| 欧美精品乱码99久久蜜桃免费| 亚洲国产日韩欧美一区二区三区,| 国产经典精品欧美日韩| 国产精品网站亚洲发布| 欧美成人屋影院在线视频观看 | 91大神福利视频网| 国产女人18毛片水真多精选| 日日躁夜夜躁狠狠操| 夜夜躁婷婷av蜜桃妖| 精品日本少妇久久久| aaaa级少妇高潮在线观看| 亚洲国产综合久久精品| 无码人妻丰满熟妇区五路| 日本少妇人妻凌辱在线| 人人妻人人爽人人爽欧美一区| 亚洲欧美另类丝袜另类自拍| 久久99精品热在线观看| 亚洲另类欧美综合久久| 午夜92福利1000| 97人妻在线视频自拍| 91九色人妻在线播放| 日韩加勒比精品在线看| 91精品在线视频免费视频| av资源中文字幕在线观看| 欧美亚洲愉拍一区二区三区| 欧美极品少妇高潮喷水| 亚洲天堂av最新在线| 亚洲另类激情视频在线看| 亚洲国产电影的一区| 亚洲综合另类欧美久久| 可以免费观看日韩av| 松本菜奈实最新av在线| 日本不卡 中文字幕| 精品日本少妇久久久| 国产视频1区2区3区| 国产午夜羞羞一区二区三区| 91九色91在线视频| 亚洲成人av在线一区二区| 人人妻人人爽人人爽欧美一区| 天天干天天操天天日天天日| 中文字幕 一区二区在线观看| 国产人妻熟女ⅹxx丝袜| 搞乱在线在线观看视频| 国产精品剧情在线亚洲| 国产福利一区二区三区在线观看| 最新久久这里只有精品| 中文字幕 人妻 熟女| 91精品国产欧美在线| 天天综合久久无人区| 手机看片福利一区二区三区四区| 高清国产美女a一级毛片| 黑人黄色免费一级av| 二十四小时日本高清在线观看| 亚洲综合成人精品成人精品| a级片特黄免费看| 开心五月综合激情婷婷| 中文字幕 人妻 熟女| 深夜福利免费观看在线看| 亚洲激情噜噜噜久久久| 欧美aaaa性bbbbaaaa| 日本小视频一区二区| 美女福利网站在线播放| 国产精品无码无卡免费观| 中日韩又粗又硬又大精品| 国产精品久久久久精品三级18| 中文字幕麻绳捆绑的人妻| www国产亚洲精品久久久| 51vv精品视频在线观看| 久久99精品热在线观看| 亚洲最强的25个城市| 亚洲国产精品久久久久久无码| 色哟哟亚洲乱码国产乱码精品精| 国产黑色丝袜 在线日韩欧美| 18禁网站在线点击观看| 91人妻人人爽色啊啊啊| 人妻激情偷乱一区二区三区av| 999精品视频免费在线观看| aa福利影视在线观看| 丰满少妇高潮喷水视频| 美女欧美视频在线观看免费| 深夜福利免费观看在线看| 国产福利小视频在线观看网站 | 麻豆午夜激情在线观看| 亚洲综合一区二区三区四区| 日本熟妇乱妇熟色视频| 国产精品美女免费视频观看| 午夜精品久久秘?18免费观看| 亚洲同性同志一二三专区| 一区二区三区婷婷中文字幕| 高清欧美色欧美综合网站 | 狂操鸡巴小骚逼视频免费观看| 精品人妻 色中文熟女 oo| 一区二区三区高清视频3| 午夜精品小视频在线播放 | 亚洲熟妇丰满多毛xxxx网站| av大尺度一区二区三区| 日本熟妇乱妇熟色视频| 日本亚洲精品视频在线观看| 色噜噜噜噜色噜噜色合久一| 午夜五十路久久福利| 天堂在线中文字幕av| 欧美啪啪一区二区三区| 97cao在线视频| 日本东京热最新中文字幕| 人妻少妇的va视频| 欧美精品一区二区三区观看| 国产伦理二区三区在干嘛呢| 顶级欧美色妇xxxx| 天堂网成人av电影| 熟妇人妻av无码中文字幕| 91porny九色视频偷拍| 黄很色很在线免费视频网站| 亚洲欧美精品日韩偷拍| 亚洲第一区av中文字幕| 美女张开腿给男人桶爽的软件| 性感美女人妻久久久| 麻豆午夜激情在线观看| 国内精品一区二区2021在线| 午夜精品小视频在线播放| 羞羞漫画无限免费观看秋蝉| 日韩激情亚洲国产欧美另类激情| 亚洲图片另类综合小说| 亚洲制服丝袜在线看| 91超碰九色porny| 精产国品一二三77777| 91精品视频在线观看视频| 欧美一区二区播放视频| 欧美最新一区二区三区| 岳母的诱惑电影在线观看| 伊人免费观看视频一| 欧美成人红桃视频在线观看| 亚洲av网站一区二区三区| 国产福利三级在线观看| 麻豆出品视频在线观看| 日韩美精品成人一区二区三区四区| 中文字幕麻绳捆绑的人妻| 国产av嗯嗯啊啊av| 骚穴被阴茎插免费视频| 99热99这里免费的精品| 91污污在线观看视频| 可以免费观看日韩av| 亚洲国产日韩精品在线| 久久久久久久久久久久久国产| 青青青青午夜手机国产视频| 亚洲精品国产99999| 最近最新欧美日韩精品| 亚洲乱码av一区二区蜜桃av| 国产自拍偷拍视频在线免费观看| 天天躁狠狠躁狠狠躁性色| 最新激情中文字幕视频| 欧美久久一区二区伊人| 亚洲无人区乱码中文字幕一区| 最新日韩中文字幕啪啪啪| 欧美日本在线免费视频| 久99久视频免费观看中文字幕| 色欲AV蜜桃一区二区三| 狠狠操深爱婷婷综合一区| 美女张开腿给男人桶爽的软件| 黑吊操欧美极品美女| 亚洲午夜国产末满十八岁勿进网站 | 内地精品毛片在线观看| 最近在线中文字幕免费| 成熟了的熟妇毛茸茸| 亚洲少妇视频在线观看| 日日夜夜免费视频精品| 亚洲|久久久久久一二三区丝袜| 人妻激情偷乱一区二区三区av| 啊不行啊操逼好爽大鸡吧视频| 亚洲男人天堂最新网址大全| 亚洲自拍偷拍一区二区中文字幕| —区二区三区女厕偷拍| 网友自拍第一页99热| 大成色亚洲一二三区| 久久热在线免费观看| 农村大炕有肉大屁股熟妇| 亚洲第一区av中文字幕| 亚洲乱熟女一区二区三区影片| 99热在线只有的精品| 人妻超清中文字幕在线乱码| 天天日 天天舔 天天射| 久久综合狠狠综合久久综| 亚洲在线观看中文字幕av| 18禁男女啪啪啪无遮挡| 日韩A级毛片免费视频| 久久久久久久精品乱码| 成人资源中文在线观看| 中文字幕亚洲乱码精品无限| 亚洲黄色免费在线观看网站| 国产美女主播av在线| 午夜福利国产精品久久久久| 日韩欧美黄色免费网站| 亚洲gay视频在线观看| 国产精品中文字幕丝袜| 亚洲国产电影的一区| 亚洲理论在线a中文字幕97| 亚洲制服丝袜美腿在线| 99 re国产精品| 夜夜躁av麻豆男| 深夜福利免费观看在线看| 顶级欧美色妇xxxx| 最新久久这里只有精品| 日本有码精品一区二区三区| 九热精品视频在线观看| 一区二区三区国产精华液区别大吗| 60路70路日本熟妇| 日本免费人爱做视频在线观看不卡| 美女把腿张开给男的捅| 自拍偷自拍亚洲精品10p| 天天摸天天舔天天操天天日| 2020国产成人精品视频| 91中文字幕视频网站| 先锋人妻啪啪中文字幕| 成人黄色录像在线观看| 操人妻人妻天天爽天天偷| 99久久碰碰人妻国产| 日韩国产欧美久久一区| 成人av在线视频免费| 亚洲一区亚洲二区成人福利| 青娱乐不卡视频在线| 日本少妇丰满大bbb的小乳沟| 国产成人在线观看视频播放| 欧美性受黑人猛交裸体视频| 欧美日韩福利视频网| 亚洲国产精品一区二区第二页| av在线免费在线观看| 99色在线观看免费观看| 美女露阴道让男人捅| 亚洲人精品午夜射精日韩| 亚洲女人自熨在线视频| av一区二区三区四区五区在线| 人妻熟女 亚洲 一页二页| 日本在线免费观看国产精品| 亚洲人人爽人人澡起碰av| 午夜福利国产精品久久久久| 亚洲熟女乱一区二区精品成人| 91人妻人人做人人爽高清| 熟妇精品午夜久久久久| 免费观看在线中文字幕视频| 91精品综合久久久久久五月天| 美女福利视频一区二区三区四区| 欧美成人屋影院在线视频观看| 中日韩又粗又硬又大精品| 特级aaaaa黄色片| 青娱乐不卡视频在线| 18禁网站在线点击观看| 午夜在线观看一级毛| 熟女阿高潮合集一区二区| 漂亮人妻口爆久久精品| 老司机伊人99久久精品| 在线看的免费网站黄| 亚洲黑人欧美二区三区| 人妻激情综合久久久久蜜桃| 国产青青青青草免费在线视频| 日韩激情文学在线视频 | 久久人妻人人草人人爽| 免费绝清毛片a在线播放| 2020精品视频在线| 日本高清在线观看不卡视频| 日本福利片在线播放| 欧美肥妇久久久久久| 日本免费人爱做视频在线观看不卡| 天天早上头和脸出汗是怎么办| 黑人3p日本女优中出| 成人超碰一区二区三区| 天天碰天天摸天天搞| 亚洲国产综合久久精品| 一级做性色a爱片久久片| 亚洲美女a级黄色在线播放| 亚洲熟女一区二区三区250p| 免费中文三级在线观看| 亚洲精品综合欧美精品综合| 亚洲综合熟女乱中文| 亚洲精品激情视频在线观看| 日本东京热视频欧美视频| 成人av在线视频免费| 中文字幕日本一二三区| 在线观看中文字幕精品av| 亚欧洲乱码视频一二三区| 人妻中文字幕亚洲在线| 亚洲成a人片777777张柏芝| 久久99精品热在线观看| 亚洲精品中文字幕手机在线免费看| 亚av一二三在线观看| 日本一区二区三区的资源| 91大神在线免费观看视频| 蜜乳av一区二区三区免费观看| 日韩人妻中文字幕区| 河北全程露脸对白自拍| 国产成人91色精品免费看片| 38av一区二区三区| 国际日韩日韩日韩日韩日韩| 五月婷婷伊人久久中文字幕| 欧美一区二区播放视频| 一区二区三区 国产日韩欧美| 1区3区4区产品乱入视频| 午夜宅男电影av网站| 成人超碰一区二区三区| 国产精品美女免费视频观看| 国产男女无套?免费网站下载 | 日本人妻少妇xxxxxxx| 日本小视频一区二区| 日韩三级精品电影久久久久| 久久久久久高清一区| 欧洲精品在线免费观看| 国产精品成人免费电影| 97人妻av人人澡人人爽| 人妻激情偷乱一区二区三区av| 在线视频国产精品欧美| 亚洲国产精品一区51动漫| 在线看日韩av不卡| 亚洲色视频在线播放网站| 蜜桃tv一区二区三区| 日韩加勒比精品在线看| av在线播放观看h| 最新激情中文字幕视频| 日韩人妻中文字幕区| 国产农村乱子伦精精品视频| 日本欧美高清在线观看视频| 69xx精品久久久久| 天天在线播放日韩av| 一区二区三区四区视频精品免费| 69视频在线精品国自产拍| 不卡在线一区二区三区| 国产精品久久久久久成人久| 亚洲综合首页综合在线观看| 伦理在线观看未删减中文字幕| 99久久久久久久久久久久久| 狠狠干狠狠操免费视频| 美女网站视频久久精品| av无限看熟女人妻另类av| 漂亮人妻口爆久久精品| 亚洲精品一区二区gif| 熟妇人妻丰满久久久久久久| 一区二区三区婷婷中文字幕| 亚洲自拍偷拍一区二区中文字幕 | 亚洲av日韩久久网站| 在宿舍强奷两个清纯校花| 亚洲中文字幕在线av| 日本有码精品一区二区三区| 大尺度久久久久久久| 蜜桃臀av在线一区二区| 亚洲av激情综合网| 天天躁狠狠躁狠狠躁性色| 亚洲欧美小说中文字幕| 五月婷婷伊人久久中文字幕| 成人午夜高清福利视频| 美利坚合众国av天堂| 91九色91在线视频| 亚洲春色av中文字幕| 中文字幕福利视频第四页| 男人的天堂av中文字幕| 精品人妻在线激情视频| 人妻系列级片在线观看视频| 青娱乐免费最新视频| 东北老女人熟女啪啪视频| 不卡高清一区二区三区| 免费绝清毛片a在线播放| 最新久久这里只有精品| 婷婷色九月综合激情丁香| 亚洲成人自拍av在线| 亚洲免费在线不卡视频| 久久久国产精品免费视频网| 美女露阴道让男人捅| 国产精品无码无卡免费观| 亚洲最强的25个城市| 久久久久久久精品乱码| 青青草原在线播放日韩| 亚洲欧美国产人成在线| 天天日天天亲天天操| 亚洲欧美韩国日本一区二区| 狠狠操av一区二区三区| 91人妻人人爽色啊啊啊| 4438x亚洲最大的成人| 欧美大鸡吧男操女啊啊啊视频| 日韩女同与成人用品电影免费看| av天堂hezyo| 欧美区一区二区三视频| 日韩av熟妇在线观看| 漂亮人妻口爆久久精品| 福利小视频免费在线| 91精品国产综合99| 黄色网络中文字幕日本| 可在线免费观看av| 精品一区二区三区喷水内射高潮 | 视频在线 一区二区| 亚洲中文字幕无线乱码人妻精品 | 日韩最近中文在线观看| 欧美老熟妇xxoo老妇| 男人电影天堂在线观看| 午夜久久久久欠久久久久| 不卡一二三区别视频| 亚洲天堂色综合久久| av资源中文字幕在线观看| julia人妻av一区二区三区| 中文人妻av一区二区三区| 在线 制服 中文字幕 日韩| 国产探花自拍亚洲av| 在线看日韩av不卡| 国产自拍偷拍视频在线免费观看| 9999久久久久老熟妇二区| 最近最新欧美日韩精品| 亚洲一区二区精品在线播放| 久久av色噜噜ai换脸| 亚洲春色av中文字幕| 男女爱爱好爽视频免费看| 免费在线小视频你懂的| 亚洲国产精品 久久久| 台湾18禁久久久久久久激情视频| www一区二区91| 五月婷婷伊人久久中文字幕| 91人妻人人爽色啊啊啊| 一级做性色a爱片久久片| 欧美日韩综合精品无人区| 亚洲国产精品久久久久久无码| 欧美极品少妇高潮喷水| 九十九步都是爱最后一步是尊严| 最新久久这里只有精品| 91九色91在线视频| 九十九步都是爱最后一步是尊严 | 亚洲欧洲无码一区2区无码| 狂操鸡巴小骚逼视频免费观看| 亚洲激情视频在线观看免费| 国产精品无码无卡免费观| 国产91九色视频在线观看| 男女爱爱好爽视频免费看| 日本清纯中文字幕版| 青青操天堂在线观看视频| 黄色av 在线观看| 日韩无码国产一区二区| 欧美一级特黄大片在线| 亚洲成人,国产精品| 日韩人妻中文字幕区| 亚洲 偷拍 自拍 欧美| 97视频人人爱麻豆| 在线成人教育平台排名| 911美女片黄在线观看| 福利在线国产小视频| 中文字字幕在线精品乱码| 午夜精品一区二区三区不卡顿| 国产成人在线观看hd| 国产成人在线观看视频播放| 亚洲精品色图1234| 每日更新日韩欧美在线| 超级黄肉动漫在线观看 | 久久内射天天玩天天懂色| 国产精品久久久久久成人久| 岛国av成人午夜高清| 豆豆专区操逼性视频在线| 国产成人情侣av在线| 亚洲色大WWW永久网站| 91精品视频在线观看视频| 欧美日韩久久丝袜在线| 亚洲欧美另类校园春色| 西野翔人妻中文字幕中字在| 欧美日韩福利视频网| 女生抠逼自慰啊啊啊啊啊啊啊下载| 午夜福利在线不卡视频| 97精品视频,全部免费| 大片a免费观看在线视频观看| 一区二区三区国产在线成人av| 日本人妻少妇xxxxxxx| 亚洲一区二区三区国产精品电影| 另类欧美激情校园春色| 午夜国产成人精品视频观看| 国产人妻777人伦精品hd超碰| 亚洲熟女乱色一区二区三区视频| 国产男人的天堂一区| 69久久夜色精品国产69乱电影 | 蜜乳视频一区二区三区| 日韩最近中文在线观看| 涩涩黄片在线免费观看| 新香蕉视频香蕉视频2| 夏目彩春av在线看| 麻豆国产91制片厂| 欧美视频亚洲视频在线| 国产男女无套?免费网站下载| 色老头一区二区三区四区五区| 国产精品无码无卡免费观| 青青青免费手机视频在线观看| 自拍偷拍 国产激情| 亚洲综合另类欧美久久| 久久精品四虎夜夜拍拍拍| 18福利视频在线观看| 久草视频在线视频在线视频| 桃色成人开心激情网| 两个人在一起靠逼啊啊啊| 色屁屁一区二区三区在线观看| 操人妻人妻天天爽天天偷| 中日韩又粗又硬又大精品| 女生抠逼自慰啊啊啊啊啊啊啊下载| 日韩激情文学在线视频| 精品视频一区二区三区◇| 伊人精品久久一区二区| aa福利影视在线观看| 日本人妻熟妇丰满成熟HD系列| 精品美女洗澡一区二区| 亚洲乱码国产乱码精品精视频 | 日本一区二区三区区别| 91色乱一区二区三区| 精品人妻 色中文熟女 oo| 九九九九九久久久国产| 五月婷婷伊人久久中文字幕| 日韩精品视频一区二区三区在线| 猫咪亚洲中文在线中文字幕| 国产视频成人一区二区| 一级毛片特级毛片免费的| 一区二区三区四区影片| 久久久国产精品免费视频网| 天天干夜夜撸天天操| 天天爽天天操天天插| 久久久久国产精品二区| 国产91免费在线观看| 黄色av网址在线播放| 色丁香久久激情综合网| 亚洲女人自熨在线视频| 国产成人在线观看视频播放| 欧美黑人性猛交小矮人| 漂亮人妻口爆久久精品| 婷婷色九月综合激情丁香| 东京热日韩av在线| 1区3区4区产品乱入视频| 亚洲综合在线视频在线播放| 中文字幕 人妻 熟女| 欧美极品少妇高潮喷水| 亚洲 偷拍 自拍 欧美| 2021国产剧情麻豆| 亚洲国产中文字幕在线看| 国产自拍偷拍视频在线免费观看| 最新激情中文字幕视频| 国产高清自拍偷拍在线| 欧美情色av在线观看| 亚洲人妻系列在线视频| 豆豆专区操逼性视频在线 | 青青在线免费手机播放视频| 午夜福利午夜福利影院| 中文字幕亚洲乱码精品无限| 操死你美女在线视频| 日本一区二区三区区别| 上床啪啪啪免费视频| 一区二区三区国产精华液区别大吗 | 一区二区在线观看视频网站| 不卡在线一区二区三区| 日韩加勒比精品在线看| 亭亭五月天在线观看| 精品国模一区二区三区欧美| 免费看超污视频在线观看| 亚洲人妻系列在线视频| 人妻系列在线免费视频| 99免费观看在线视频| 最新福利二区三区视频| 亚洲国产精品自产拍在线观看| 视频在线+欧美十亚洲曰本| 国产一区两区三区福利小视频| 免费高清av一区二区| 亚洲人成小说网站色| 国产美女高潮精品视频| 日韩人妻一区二区三区在线观看| 91porny九色视频偷拍| 18在线观看免费观看| 美女露阴道让男人捅| 欧美亚洲愉拍一区二区三区| 午夜精品小视频在线播放| 五月婷婷伊人久久中文字幕| 大屁股熟女一区二区视频| 亚洲精品久久久人妻| 亚洲欧美不卡专业视频| 在线免费观看视频18| 国产漂亮白嫩美女在线图片| 黑人巨大精品一区二区在线| 成年男女免费视频网站无毒| 久久久久久久岛国免费观看| 福利视频导航在线观看| 国产高清自拍偷拍在线| 亚洲成人五月婷婷久久综合| 真人一进一出抽搐大尺度视频| 黑人和日本人av一区二区| 亚洲 自拍 激情 另类| 夜夜操夜夜爱夜夜摸| 97精品视频,全部免费| 国产做A爱免费视频在线观看| 七色福利视频在线观看| 69国产在线视频网站| 99久9在线视频播放| 99 re国产精品| 久久99国产中文丝袜| 波多野结衣在线一区别| 中文字幕福利视频在线一区| 91香蕉国产亚洲一二三区| 国产av精品一区二区三区久久| 亚洲欧美国产人成在线| 亚洲一区二区偷拍女厕所| 18福利视频在线观看| 日韩欧美黄色免费网站| 午夜8050免费小说| 天天爽天天操天天插| 自拍偷拍亚洲综合第一页| 欧美aaaa性bbbbaaaa| 黄色片免费网站在线| 久久久国产精品免费视频网| 福利一二三在线视频观看| 亚洲欧美激情久久久| 二十四小时日本高清在线观看| 亚洲成人av在线一区二区| 69xx精品久久久久| 一区二区三区四区影片| 麻豆国产91制片厂| 神马午夜久久电影网| 最近最新最好看的中文字幕| 天堂一区二区三区在线等| 亚欧洲乱码视频一二三区| 精品免费一区二区三区四区视频| av在线男人的天堂亚洲| 日本高清有码在线视频| avjpm亚洲伊人久久| 两个奶被揉得又硬又翘怎么回事| 亚洲一区视频中文字幕在线播放 | 97香蕉久久国产超碰| 在线观看视频免费一区二区三区| 四虎国产精品国产精品国产精品| 欧美视频免费观看777| 日韩人妻精品久久久久| 河北全程露脸对白自拍| 福利美女视频在线观看| 4日日夜夜精品视频免费| 大香焦一道本一区二区三区| 日韩欧美一区二区三区免费看| 操死你美女在线视频| 黑人巨大精品一区二区在线| 天天干天天色综合久久| 亚洲人妻系列在线视频| 人妻激情偷乱一区二区三区av| 人妻熟女 亚洲 一页二页| 一级毛片特级毛片免费的| 亚成区一区二区人妻熟女| 一区二区三区四区影片| av在线观看视频免费| 中文字幕在线免费观看成人| 九九热视频1这里只有精品| 久草视频在线视频在线视频| 色视频免费观看网址| 青青草原在线播放日韩| 一区二区三区高清视频3| 亚洲成人五月婷婷久久综合| 欧美日韩高清片在线观看| 国产av高清二区三区| 亚洲综合天堂av网站在线观看| 中文字幕欧美一区二区视频| 久久人人爽人人爽人人av东京热 | 中文字幕亚洲乱码精品无限| 夜夜爽夜夜操夜夜爱| 欧美人与动欧交视频| 老司机免费视频福利0| 亚洲欧美另类丝袜另类自拍| 真人一进一出抽搐大尺度视频| 新香蕉视频香蕉视频2| 欧美成人久久久桃色aa| 人妻系列级片在线观看视频| 欧美日本国产一区二区| 天天干夜夜操夜夜骑| 蜜桃臀少妇白色紧身裤细高跟| 色噜噜噜噜色噜噜色合久一| 九九九九九久久久国产| 日韩人妻精品久久久久| 天天干天天色综合久久| jizzjizz国产精品传媒| 日本高清久久人人爽| 精品国产av虐杀两警花 | 大奶熟妇激情操逼逼| 久久久久久a女人处女| 国产精品蝌蚪自拍视频| 欧美大鸡吧男操女啊啊啊视频| 国产igao激情在线视频入口| 日本一区二区三区区别| 性感美女人妻久久久| 可在线免费观看av| 在线中文字幕人妻av| 91九色人妻在线播放| 91美女在线观看视频| 婷婷色综合五月天视频| 69国产在线视频网站| 久久99国产中文丝袜| 人妻少妇精品二三区| 在线观看网站伊人网| 日韩成人免费观看电影| 快进来插我的逼嗯啊视频| 大香蕉尹人在线最新| 一区二区三区国产在线成人av| 国产剧情av在线免费观看| 日韩免费黄色片在线观看| 老色鬼精品视频在线观看播放| 不卡一区二区视频在线| 中文字幕 一区二区在线观看| 亚洲在线观看中文字幕av| 国产精品剧情av在线播放| 黄色大片一级老太太操逼| 白白色在线免费视频发布视频| 国产经典精品欧美日韩| 3344永久在线观看视频下载| 十八禁黄色免费污污污亚洲| 日本高清 中文字幕| 久久99热精品免费观看视| 一区二区在线观看视频网站| 美女妩媚午夜诱惑网站| 天天躁狠狠躁狠狠躁性色| 快色视频在线观看免费| 欧美一区二区三区视频看| 天堂一区二区三区在线等| 91精品久久久久久久久99蜜臀| 得得爱在线视频观看| 亚洲乱熟女一区二区三区影片| 欧美成人短视频在线播放| 国产av在线免费视频| 亚洲va999天堂va| 亚洲精品国产99999| 亚洲熟妇在线视频观看| 偷拍熟女大胆免费视频| 91大神在线免费观看视频| 亚洲高清免费在线观看视频| 欧美一区二区三区爽爽| 天堂在线中文字幕av| 欧美精品一区二区三区观看| 欧美久久一区二区伊人| 放荡人妻极品少妇全集| av激情四射五月婷婷| 插鸡视频免费网站在线播放| 啪啪啪网站免费看视频| 亚洲中文字幕无线乱码人妻精品| 天堂网免费在线电影| 搞乱在线在线观看视频| 1区3区4区产品乱入视频| 亚洲色图日韩在线视频观看| 国产熟妇色xxⅹ交白浆视频| 亚洲一区二区精品三区视频| a级片特黄免费看| 欧美插插插插插插| 波多野结衣在线一区别| 91进入蜜桃臀在线播放| 福利小视频免费在线| 97精品人妻免费视频| 日本在线免费观看国产精品| 中文字幕亚洲乱码精品无限| 日韩成人在线电影首页| 国内自拍第一区二区三区| 日韩免费黄色片在线观看| 狠狠干狠狠操免费视频| 欧美国产精品久久久免费| 18禁网站在线点击观看| 久久久久久久久久久久久国产| av在线中文字幕在线| 国产高清自拍偷拍在线| 精品不卡一区二区三区| 亚洲最强的25个城市| 在线观看中文字幕视频成人| 亚洲欧美国产一本综合首页| 亚洲欧美国产一本综合首页| 中文字幕人妻一区二区视频系列| 国产精品性感美女视频| 偷拍欧美日韩另类图片| 国产91黑丝小视频在线观看| 丰满放荡熟妇在线播放 | 福利一二三在线视频观看| 精品久久久久久久久久久久久| 麻豆午夜激情在线观看| 天天操天天舔天天射天天日天天干| 97超碰人人爽人人做| 91精品麻豆91夜夜骚| 91亚洲精品久久蜜桃| 天天插天天透天天爽| 青青操91美女国产| 欧美vs亚洲vs日韩| 午夜福利国产精品久久久久| 人妻色综合aaaaaa网| 91久久精品美女高潮喷水白浆| 亚洲综合首页综合在线观看| 狠狠操深爱婷婷综合一区| 午夜福利国产精品久久久久 | 欧美丝袜亚洲国产日韩| 99亚偷拍自图区亚洲| 裸露视频免费在线观看| 亚洲色视频在线播放网站| 18禁网站在线点击观看| 无人区一码二码三码区别在哪| 欧美色视频网址大全| 三区美女视频在线观看| 青青青在线视频观看97| 午夜8050免费小说| 亚欧洲乱码视频一二三区| 中国精品人妻一区二区| 黄版视频在线免费观看| 五月婷婷激情视频网| 中文字幕在线观看亚洲情色| 日本老熟老熟妇七十路| 午夜福利在线不卡视频| 天堂在线中文字幕av| 国产精品剧情av在线播放| 丝袜美女诱惑佐佐三上| 亚洲av日韩久久网站| 91久久久精品成人国产| 两个人在一起靠逼啊啊啊| 日本一区二区三区的资源| 福利视频导航在线观看| 夜夜爽夜夜操夜夜爱| avjpm亚洲伊人久久| 精品美女洗澡一区二区| 三级欧美日韩一区二区三区| 亚洲熟女少妇中文字幕系列| 深夜福利免费观看在线看| 亚洲a区在线免费观看| 久久国产半精品99精品国产| 亚洲成人中文无码在线| 中文字幕在线字幕乱码怎么设置| av人摸人人人澡人人超碰小说| 亚洲综合首页综合在线观看| 丰满人妻熟女aⅴ一区| 国产探花自拍亚洲av| 91人妻人人爽色啊啊啊| 亚洲欧美成人激情在线| 青青免费观看视频| av丝袜免费在线观看| 伊人网在线免费观看| 五月婷婷激情视频网| 18禁网站在线点击观看| 国产资源网站在线播放| 亚洲国产精品青青草| 天天干天天操天天日天天日| 天天做天天日天天搞| 亚洲精品9999蜜桃| 天天插天天干天天狠| 最近最新欧美日韩精品| avgo成人短视频| 超碰在线免费观看视频97| 在线免费观看视频18| 亚洲gay视频在线观看| 女人高潮潮呻吟喷水网站| 日本特级黄片免费观看| 77亚洲视频在线观看| 美国十次了亚洲天堂网国产| 加勒比不卡在线视频| 日本黄页在线观看视频| 亚洲成人偷拍自拍在线| 黑川堇人妻88av| 亚洲精品9999蜜桃| 干逼又爽又黄又免费的视频| 日韩A级毛片免费视频| 亚洲欧美成人午夜一区二区| 人妻色综合aaaaaa网| 亚洲中文字幕最新地址| 日韩女同与成人用品电影免费看| 懂色av之国产精品| 欧美日本亚欧在线观看| 午夜精品久久久久久久精品乱码| 欧美在线视频不卡一区| 天天在线播放日韩av| 黄很色很在线免费视频网站| 加勒比不卡在线视频| 日本一区二区高清av中文| 男人的天堂av中文字幕| 97精品国产91久久久| 国长拍拍视频免费孕妇| 女人扒开逼让男人操| 日本少妇熟女乱码一区二区 | 欧美视频免费观看777| 美女av色播在线播放| 男女插鸡巴视频软件| 91porny九色视频偷拍| 免费成人av麻豆| 亚洲制服丝袜美腿在线| 国产黄色主播网址大全在线播放 | 午夜夫妻性生活视频| 成人精品影视一区二区| 亚洲一区二区中文字幕久久| 一区二区三区高清视频3| 国产激情免费在线视频| 日本男女免费福利视频| ass亚洲熟女ass| 成人超碰一区二区三区| 大陆中文字幕视频在线| 99久9在线视频播放| 99久久久久久亚洲精品免费| 欧美日本在线免费视频| 日本老熟妇av老熟妇| 亚洲欧美精品海量播放| 九九九九九久久久国产| 在线能看视频你懂的| 91污污在线观看视频| 不卡在线一区二区三区| 欧美强奸视频在线观看| 中文字幕人妻精品精品| 成人十欧美亚洲综合在线| 性色蜜桃臀x88av天美传媒| 2020国产成人精品视频| 成人精品影视一区二区| 99精品久久一区二区| 午夜精品秘一区二区三区| 91进入蜜桃臀在线播放| 青青青青青爽视频在线| 精品久久久久久久久久久久久| 男人的天堂aⅴ在线| 精品美女洗澡一区二区| 快使劲弄我视频在线播放 | 亚洲 偷拍 自拍 欧美| 国产精品蝌蚪自拍视频| 大香蕉在线欧美在线视频| 亚洲黑人欧美二区三区| 日本男女免费福利视频| 五月天男人的天堂中文字幕| 岳母的诱惑电影在线观看| 欧美一级aaaaaaa片| 97精品人妻免费视频| 美女福利网站在线播放| 国产视频成人自拍蝌蚪视频 | 欧美日韩黄片免费在线观看| 天天操天天射天天操天天日| 91青青青国产免费高清| 日本高清在线观看不卡视频| 亚洲精品激情视频在线观看| 亚洲激情噜噜噜久久久| 最新日韩中文字幕免费在线观看| 在线能看视频你懂的| 欧美视频亚洲视频在线| 欧美啪啪一区二区三区| 国产激情免费在线视频| 黑人3p日本女优中出| 午夜精品秘一区二区三区| 91超精品碰国产在线观看| 黑人大吊大战亚洲女人。| 久久久人妻免费视频| 免费观看在线中文字幕视频| 青青青在线视频免费播放| 国产主播诱惑毛片av| 51vv精品视频在线观看| 亚洲人妻系列在线视频| 视频免费在线观看网站| 熟女俱乐部jukujoclub| 视频免费在线观看网站| 免费中文三级在线观看| 大成色亚洲一二三区| 久久久久久高清一区| 久久综合狠狠综合久久综| av在线男人的天堂亚洲| 国产高清在线观看av| 黄色片免费网站在线| 国产精品久久久99| a级片特黄免费看| 成人精品影视一区二区| 精品美女洗澡一区二区| 2026天天操天天干| 五月天男人的天堂中文字幕| 天天操天天舔天天做| 国产清纯一区二区在线观看| 九九热精品视频在线播放| 国际日韩日韩日韩日韩日韩| 18在线观看免费观看| 亭亭五月天在线观看| 国产美女主播av在线| 成人超碰一区二区三区| 公侵犯人妻中文字幕巨| 大乳丰满人妻中文字幕韩国hd| 美女露阴道让男人捅| 182tv精品免费在线观看| 青青青在线视频观看97| 欧美成人屋影院在线视频观看| 亚洲av网站一区二区三区| 日韩在线 中文字幕| 欧美黄色性视频网站| 大香蕉伊人97在线| 99精品久久精品一区二区| 亚洲av网站一区二区三区| 国产成人在线观看视频播放| 猫咪亚洲中文在线中文字幕| 日韩加勒比精品在线看| 欧美日本在线免费视频| 911精产国品一二三产区区| 最新日韩中文字幕免费在线观看| 日本老熟老熟妇七十路| 河北全程露脸对白自拍| 日本不卡视频一二三区| 97人妻人人揉人人躁人人夜夜爽 | 国产黄色主播网址大全在线播放| 女人扒开逼让男人操| aaaa级少妇高潮在线观看 | 亚洲中文字幕在线视频观看二区| 男生用大肌巴操美女骚穴| 97超碰人人爽人人做| 女女抠逼白虎白丝袜| 91污污在线观看视频| 日韩三级黄色大片在线观看| 91精品资源在线观看| 亚洲国产精品 久久久| 天天干天天日天天弄| 天天操天天搞天天操| 亚洲免费在线不卡视频| 男人用大鸡巴狂操女人肉穴| 4438x亚洲最大的成人| 日本一区二区三区的资源| 青青在线免费手机播放视频| 92麻豆一区二区三区| 亚洲成人三级黄色片| 天天搞天天操天天干| 成人av中文字幕在线看| 中文字幕在线字幕乱码怎么设置| 亚洲综合色一区二区三区| 日韩成人在线电影首页| 综合久久伊人久久88| av里面的动作是真进去吗| 亚洲黄色成人一级片| 日本亚洲精品视频在线观看| 天天操天天搞天天操| 国产黄色主播网址大全在线播放| 成人免费视频现网站99在线观看| aaaa级少妇高潮在线观看| 十八禁黄色免费污污污亚洲| 中文字幕熟女人妻一区| 不卡一二三区别视频| 中文字幕福利视频在线一区| 天天爱天天日天天爽| 亚洲成人,国产精品| 在线成人教育平台排名| 日韩三级黄色大片在线观看| 99re这里是国产精品首页| www,日韩av,com| 中文字幕日韩人妻在线三区| 久久久久久a女人处女| 青娱乐这里只有精品| 午夜亚洲国产精品中字| 美国伦理片午夜理论片| 91九色国产在线视频| 一区二区三区不卡免费视频网站| 黄色大片一级老太太操逼| 1级黄色片在线观看| 日本国产亚洲欧美色综合| 中文字幕免费啪啪啪| 亚洲va999天堂va| 大尺度久久久久久久| 在线有码人妻自拍视频| 一区二区三区四区久久久久韩日| 亚洲精品激情视频在线观看| 丰满少妇人妻一区二区三区蜜桃| 国产一区二区三区四区精| 神马午夜久久电影网| 91中文字幕视频网站| 69精品人妻久久久久久久久久久| 夜夜操天天干夜夜操| 黄色片黄色片黄色片黄色片黄色| 中文字幕丰满子伦无码专区| 欧美成人短视频在线播放| 狠狠操狠狠操狠狠插| 亚洲色图日韩在线视频观看| 国产福利三级在线观看| 久久久久性感美女偷拍视频| 天天操天天舔天天爽| 老牛影视在线一区二区三区| 青青国产95免看视频| 欧美一区二区三区爽爽| 日韩激情文学在线视频| 日韩欧美国产一区二区在线观看| 99久久精品视频16| 亚洲韩精品一区二区三区 | 嗯~嗯~啊啊啊~高潮了软件| 免费观看在线中文字幕视频| 免费24小时人妻视频| 日韩一区二区在线播放观看| 97人妻在线视频自拍| 国产精品亚洲精品亚洲| 在线成人教育平台排名| 欧美日韩综合精品无人区| 欧洲成熟女人色惰片| 女生裸体视频免费网站| 全球高清中文字幕av| 99久久国产精品免费热| 宅男噜噜噜66国产在线观看| 69久久夜色精品国产69乱电影 | 亚洲激情噜噜噜久久久| 九九六视频,这里只有精品| 国产高清自拍偷拍在线| 91精品国产人妻麻豆| 97人妻av人人澡人人爽| 国产激情视频在线观看的| 性高潮视频在线观看日韩| 精品精品精品精品精品污污污污 | 国长拍拍视频免费孕妇| 美女激情久久久久久久| 亚洲成人三级黄色片| 亚洲成人中文无码在线| 中文字幕人妻一区色偷偷久久 | 一区二区三区国产精华液区别大吗| 自拍偷拍亚洲综合第一页| 青娱乐免费最新视频| 天天色 天天操 天天好逼| 天天在线播放日韩av| 四虎精品久久免费最新| 亚洲中文字幕无线乱码人妻精品| 天天色天天射天天日天天干| 内地精品毛片在线观看| 天天干夜夜操91视频网站| 日本欧美亚洲国产啊啊啊| 在宿舍强奷两个清纯校花| 青青在线视频看看| 丰满人妻熟女aⅴ一区| 亚洲一区二区精品三区视频| 国产白丝一区二区三区av| 欧美日韩高清片在线观看| 最新国产精品久久精品app| 夜夜操夜夜爱夜夜摸| 午夜国产精品免费视频| aaaa级少妇高潮在线观看 | 69精品人妻久久久久久久久久久| 亚洲 偷拍 自拍 欧美| 91精品资源在线观看| 久久久久久免费观看av| 亚洲熟女一区二区三区250p| 欧美久久一区二区伊人| 港台美女明星av天堂| 91福利高清在线播放| 91porny九色视频偷拍| 人妻激情综合久久久久蜜桃| 天天搞天天操天天干| 高潮喷水在线视频观看| 97人妻av人人澡人人爽| 国产高清自拍偷拍在线| 五十岁熟妇高潮喷水| 成年人免费福利在线| 妈妈的朋友2中文字幕在线| 亚洲国产精品自拍偷拍视频在线| 黑人和日本人av一区二区| 青青青在线视频免费播放| 国产成人综合久久婷婷| 天天搞天天操天天干| 天天干夜夜操夜夜骑| 久久中文字幕av一区二区| 公侵犯人妻中文字幕巨| 不用付费特黄特色亚洲特级黄色片| 欧美一区二区三区视频看| 最新日韩中文字幕免费在线观看| 日本少妇人妻中文在线| 国产精品久久久久久成人久| 精品国产污污污污免费观看| 国产成人情侣激情视频| 国产黄色主播网址大全在线播放| 99精品视频在线在线观看| 国产在线小视频一区二区| 黄版视频在线免费观看| 亚洲综合色一区二区三区| lutu玩弄人妻短视频| 天天操天天舔天天爽| 老熟女xxxⅹhd老熟女性| 亚洲中文字幕在线av| 公侵犯人妻中文字幕巨| 久久人妻人人草人人爽| 日产国产欧美精品另类| 东京热日本一区二区三区| 婷婷色九月综合激情丁香| 久久久久九九九九九12| 99久久人人爽亚洲精品美女| 韩国毛片w妈妈的朋友7| 日韩三级精品电影久久久久| 亚洲综合天堂av网站在线观看| 亚洲黄色成人一级片| 欧美视频免费观看777| 欧美大胆a级视频秒播| 国产av嗯嗯啊啊av| 亚洲第一中文字幕成人| 国模伊人久久精品一区二区三区| av天堂新资源在线| 在线观看免费啪啪啪| 久久久西西gogo日本美女人体| 亚洲女人自熨在线视频| 日韩av水蜜桃一区二区三区| 九色porny91国产| 午夜精品一区二区三区不卡顿| 啊~插得好快别揉我胸了视频 | 免费绝清毛片a在线播放| 日韩在线 中文字幕| 国产精品亚洲精品亚洲| 欧美大胆a级视频秒播| 夫亡人妻被强干中文字幕| 精品欧美黑人一区二区三区| 色999日韩偷自拍拍免费| 91精品久久久久久久99蜜月| 午夜国产成人精品视频观看| 最近在线中文字幕免费| 国产青青青青草免费在线视频| 午夜福利午夜福利影院| 午夜精品秘一区二区三区| 亚洲国产精品自产拍在线观看| 极品内射老女人操逼视频| 欧美黑人1区2区3区| 日韩免费黄色片在线观看| 全彩漫画口工18禁| 日本韩国福利在线播放| 深夜福利免费观看在线看| 国产经典精品欧美日韩| 男插女视频大全免费| 欧美啪啪一区二区三区| 77亚洲视频在线观看| xxnxx国产美女| 激情久久在线免费观看视频| 女人扒开逼让男人操| 亚洲午夜国产末满十八岁勿进网站| 夜夜操夜夜爱夜夜摸| 玖辛奈18禁同人污本子| 国产精品中文字幕丝袜| 午夜精品久久秘?18免费观看| 亚洲av在线免费播放| 精品人妻 色中文熟女 oo| 日本久久久久久黄色| 68福利精品在线视频| 99 re国产精品| 成人18禁高潮片免费日本| 亚洲天堂色综合久久| 91精品夜夜夜一区二区蜜桃| 亚洲综合首页综合在线观看| 亚洲综合天堂av网站在线观看| 可以免费观看日韩av| 欧美成人短视频在线播放| 中文字幕日本一二三区| 神马不卡视频在线视频| 亚洲欧美精品海量播放| 上床啪啪啪免费视频| 国产一区二区手机在线观看| 亚洲激情噜噜噜久久久| 国内自拍第一区二区三区| 熟妇高潮久久久久久久| 手机看电影一区二区三区| 免费在线观看亚洲福利| 国产人妻熟女ⅹxx丝袜| 午夜精品视频免费观看| 美女扒开逼逼给你看| 最新日韩中文字幕啪啪啪| 自拍偷拍亚洲综合第一页| 夜夜人人干人人爱人人操| 亚洲一级熟妇丰满的女人| 日本美女爱爱视频网站| 欧美aaaa性bbbbaaaa| 精品免费一区二区三区四区视频| 1区3区4区产品乱入视频| 狠狠操深爱婷婷综合一区| 精品人妻人人做人人爽| 久久久久久免费观看av| 人妻超清中文字幕在线乱码| 国产女主播在线观看一区| 亚洲黄色免费在线观看网站| 青青青青青爽视频在线| 黄色av日韩在线观看| 可在线免费观看av| 色屁屁一区二区三区在线观看| 一级做性色a爱片久久片| 亚洲免费午夜污福利| 爱搞视频在线观看视频91| 福利一二三在线视频观看| 5d蜜桃臀女无痕裸感| 91福利高清在线播放| 天天看片天天摸天天操| 色噜噜噜噜色噜噜色合久一| 午夜久久人妻一级内射av网址 | 天天操天天射天天操天天日| 五月婷婷伊人久久中文字幕| 精品视频一区二区三区◇| 日本人妻熟妇丰满成熟HD系列| 男生和女生羞羞91在线看| 天堂av国产av伦理av| 亚洲在线免费观看18| 欧美成人红桃视频在线观看| 蜜臀一区二区日韩美女少妇视频 | 无人区一码二码三码区别在哪| 黑人侵犯人妻森泽佳奈| tushy一区二区三区视频| 亚洲成人中文无码在线| 一区二区三区五区六区| 欧美情色av在线观看| 日日躁夜夜躁狠狠操| 亚洲综合另类欧美久久| 人妻熟女 亚洲 一页二页| 91国产精品乱码久久久久久| 久久久久高潮白浆久久| 99福利一区二区视频| 欧美色区国产日韩亚洲区| 午夜亚洲国产精品中字| tushy一区二区三区视频| 五月天天堂视频在线| 亚洲成人自拍av在线| 中文字幕中文字幕在线中…一区| 欧美情色av在线观看| 国际精品熟女一区二区| 欧美区日本区国产区| 在线视频自拍第三页| 成年人免费福利在线| avjpm亚洲伊人久久| 欧美极品少妇高潮喷水| 在线免费观看欧美小视频| 午夜久久人妻一级内射av网址| 国产人妻777人伦精品hd超碰| 中文字幕在线观看亚洲情色| 自拍偷拍 国产激情| 午夜福利午夜福利影院| 瑟瑟干视频在线观看| 国产午夜在线播放视频| www一区二区91| 中文乱码字幕人妻熟女人妻| 国产av高清二区三区| 亚洲自拍偷拍av在线| 91亚洲精品久久蜜桃| 性感人妻 中文字幕| 干逼又爽又黄又免费的视频| 天天操,天天射,天天爽| 久久综合狠狠综合久久综| 两个人在一起靠逼啊啊啊| 久久久久久久岛国免费观看| 性感美女人妻久久久| 鸡巴在里面福利视频在线观看| 大尺度久久久久久久| 91精产国品一二三产区区别网站| 国产视频1区2区3区| 黑人3p日本女优中出| 在宿舍强奷两个清纯校花| 亚洲国产日韩a在线欧美| 久久久西西gogo日本美女人体| 天天在线播放日韩av| 韩国一级片最火爆中文字幕| 少妇熟女天堂网av| 人妻系列级片在线观看视频| av激情四射五月婷婷| av 一区二区三区 熟女| 婷婷一区二区三区五月丁| 色就色综合偷拍区欧美在线| 九一精品人妻一区二区三区| 午夜夫妻性生活视频| 男人的天堂av中文字幕|