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

查看: 3561  |  回復(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è)置，論壇里面有別人寫好的優(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)化是不是，寫出被優(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)槲矣檬謱懙淖鴺?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)槲矣檬謱懙淖鴺?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)槲矣檬謱懙淖鴺?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)槲矣檬謱懙淖鴺?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) (正式寫手)
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ù)版主 (著名寫手)
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ù)版主 (著名寫手)
uuv2010(金幣+1): 多謝提示 2012-02-25 17:23:53
初步建粒子位置的時(shí)候盡可能的準(zhǔn)確些，然后再優(yōu)化，ISIF的值也要設(shè)置，論壇里面有別人寫好的優(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) (正式寫手)
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) (正式寫手)
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ù)版主 (著名寫手)
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ā)表
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +3 脫穎而出 2026-03-16 6/300 2026-03-17 21:19 by z1z2z3879
[考研] 材料與化工求調(diào)劑 +6 為學(xué)666 2026-03-16 6/300 2026-03-17 20:15 by peike
[考研] 326求調(diào)劑 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 本人考085602 化學(xué)工程 專碩 +16 不知道叫什么！ 2026-03-15 18/900 2026-03-17 17:05 by ruiyingmiao
[考研] 梁成偉老師課題組歡迎你的加入 +8 一鴨鴨喲 2026-03-14 10/500 2026-03-17 15:07 by 一鴨鴨喲
[考研] 285化工學(xué)碩求調(diào)劑（081700） +9 柴郡貓_ 2026-03-12 9/450 2026-03-17 10:18 by Sammy2
[考研] 環(huán)境工程調(diào)劑 +6 大可digkids 2026-03-16 6/300 2026-03-16 17:16 by barlinike
[考研] 285求調(diào)劑 +6 ytter 2026-03-12 6/300 2026-03-16 15:05 by njzyff
[考研] 一志愿華中師范071000，325求調(diào)劑 +6 RuitingC 2026-03-12 6/300 2026-03-16 14:50 by 可淡不可忘
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +5 ZY-05 2026-03-13 6/300 2026-03-14 10:51 by Jy?
[考研] 0703化學(xué)調(diào)劑 +4 快樂(lè)的香蕉 2026-03-11 4/200 2026-03-13 22:41 by JourneyLucky
[考研] 336求調(diào)劑 +6 Iuruoh 2026-03-11 6/300 2026-03-13 22:06 by JourneyLucky
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9 2026調(diào)劑俠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 301求調(diào)劑 +6 Liyouyumairs 2026-03-11 6/300 2026-03-13 20:11 by JourneyLucky
[考研] 材料專碩350 求調(diào)劑 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
[考研] 328化工專碩求調(diào)劑 +4 。，。，。，。i 2026-03-12 4/200 2026-03-13 14:44 by JourneyLucky
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3 T123 tt 2026-03-12 3/150 2026-03-13 10:49 by houyaoxu
[考研] 290求調(diào)劑 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
[考研] 化工學(xué)碩306求調(diào)劑 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考研] 321求調(diào)劑（食品/專碩） +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
信息提示
請(qǐng)?zhí)钐幚硪庖?jiàn)
港台美女明星av天堂| 久久久久性感美女偷拍视频| yy4080黄色片| 国产漂亮白嫩美女在线图片| 亚洲熟女少妇中文字幕系列| av资源中文字幕在线观看| 丝袜美女诱惑佐佐三上| 手机视频在线观看一区| 日本熟妇乱妇熟色视频| avgo成人短视频| 九色91操最新在线观看网址| 亚洲综合一区二区三区四区| 天天干夜夜爽狠狠操| 麻豆国产91制片厂| 亚洲成人偷拍自拍在线| 国产一区二区手机在线观看| 99精品久久一区二区| 青青青在线视频免费播放| 91精品国产91久久久久久密臀| 亚洲日本欧美韩国另类综合| 人人妻人人狠人人爽| 午夜精品小视频在线播放| 欧美成人区一区二区三| 亚洲综合另类欧美久久| 91精品夜夜夜一区二区| 日本高清在线观看不卡视频| 亚洲欧美另类校园春色| 亚洲欧美精品海量播放| 精产国品一二三产品区别97| 中文字幕丰满子伦无码专区| 女生裸体视频免费网站| 国产成人在线观看hd| 不卡高清一区二区三区| 日本在线免费观看国产精品| 五十岁熟女高潮喷水| 蜜臀一区二区日韩美女少妇视频| 日韩加勒比精品在线看| 99女福利女女视频在线播放| 亚洲综合色一区二区三区| 美利坚合众国av天堂| 91亚洲精品久久蜜桃| 亚洲图片另类综合小说| 国产成人情侣激情视频| 91人妻人人爽色啊啊啊| 人妻在线中文视频视频| 亚洲美女午夜激情视频在线观看| 伊人精品久久一区二区| 两个人在一起靠逼啊啊啊| 特级aaaaa黄色片| 国产大桥未久一区二区| 日本一区二区高清av中文| 国产av高清二区三区| 国产成人91色精品免费看片| 九一精品人妻一区二区三区| 操死你美女在线视频| 92午夜免费福利视频www| 大片a免费观看在线视频观看| 午夜宅男电影av网站| 成熟了的熟妇毛茸茸| 日本老熟妇av老熟妇| 天天干天天色综合久久| 美国十次了亚洲天堂网国产| 99久久免费播放在线观看视频| 亚洲一区在线视频观看地址| 操人妻人妻天天爽天天偷| 夏目彩春av在线看| 欧美日韩久久丝袜在线| 久久国产精品久精国产爱 | 欧美日韩成人高清中文网| 伊人网在线观看 视频一区| 美女把逼扒开让男人桶| 日本免费人爱做视频在线观看不卡| 色欲天天媓色媓香视频综合网| 麻豆白洁少妇在线播放| 中国精品人妻一区二区| 亚洲激情视频在线观看免费| 性色蜜桃臀x88av天美传媒| 国产成人91色精品免费看片| 国产激情一区二区视频| 男人的天堂aⅴ在线| 国内销魂老女人老泬| 人人妻人人爽人人爽欧美一区| 青青青青青爽视频在线| 亚洲经典av中文字幕| 国产igao激情在线视频入口| 91九色人妻在线播放| 国产高清视频www夜色资源| 国产自拍偷拍视频在线免费观看| 中文字幕在线观看亚洲情色| 久久久久夜色国产精品电影| 一区二区三区五区六区| 国产午夜在线播放视频| 老司机免费视频福利0| 在宿舍强奷两个清纯校花| 亚洲同性同志一二三专区| 玖玖资源站在线观看亚洲| 狠狠操狠狠操狠狠插| 天天爽天天操天天插| 69av精品国产探花| 亚洲春色av中文字幕| 亭亭五月天在线观看| jizzjizz国产精品传媒| 欧洲精品在线免费观看| 在线观看中文字幕少妇av| 亭亭五月天在线观看| 亚洲色视频在线播放网站| 91激情四射婷婷综合| 亚洲黑人欧美二区三区| aa福利影视在线观看| 在线人成视频免费观看尤物| 青娱乐免费最新视频| 亚洲AV无码久久精品国产一区老| 夜夜人人干人人爱人人操| 欧美不卡一二三区精品| 日韩人妻一区二区三区在线观看| 99久久碰碰人妻国产| 自拍偷自拍亚洲精品10p| 鸡巴在里面福利视频在线观看| 欧美男女一区二区三区| 男人的天堂aⅴ在线| 中文字幕在线观看av观看| 欧美成人久久久桃色aa| 日本少妇精品免费视频| 日本丰满熟妇浓密多毛| 另类欧美激情校园春色| 天天操天天舔天天爽| 最近日韩免费在线观看| 91 精品视频在线看| 99亚偷拍自图区亚洲| 日本老熟老熟妇七十路| 亚洲欧美另类丝袜另类自拍| 亚洲熟妇在线视频观看| 亚洲欧洲一区二区三区在线| 国产乱码有码一区二区三区| 国产精品久久久久精品三级18| 日本午夜福利免费在线播放| 亚洲码av一区二区三区| 午夜宅男电影av网站| 蜜桃臀少妇白色紧身裤细高跟| 91亚洲精品久久蜜桃| 国产精品蝌蚪自拍视频| 天天插天天操天天射天天干| 松本菜奈实最新av在线| 欧美vs亚洲vs日韩| 亚洲精品9999蜜桃| 国产女主播在线观看一区| 成人资源中文在线观看| 国产成人综合久久婷婷| 国产精品剧情av在线播放| 五月激情婷婷四射基地| 手机看片福利一区二区三区四区| 60路70路日本熟妇| 亚洲国产美女主播在线观看| 久久热在线免费观看| 快色视频在线观看免费| 色视频免费观看网址| 自拍偷自拍亚洲精品10p| 欧美久久蜜臀蜜桃资源吧| 午夜国产精品免费视频| 天天摸天天干夜夜操| 亚洲资源在线免费观看| 亚洲国产日韩精品在线| 日本东京热最新中文字幕| 精品国模一区二区三区欧美| 天天插天天透天天爽| 天天透天天舔天天操| 国产激情视频在线观看的| 老熟妇一区二区三区v∧88| 中文字幕在线免费观看成人| 交换的一天中文字幕在线视频| 天天操天天舔天天爽| 92午夜免费福利视频www| 大秀成年人国产精品视频| 中文在线字幕免费观看日韩视频| 欧洲精品在线免费观看| 亚洲国产精品一区51动漫| 中文人妻av一区二区三区| 丝袜美腿日韩av一区| 上床啪啪啪免费视频| 天海翼亚洲一区在线观看| 亚洲美女黄色福利视频网站大全| 九九六视频,这里只有精品| 国产三级自拍视频在线观看网站| 午夜精品小视频在线播放| 成熟了的熟妇毛茸茸| 中文字幕熟女乱一区二区| 亚洲午夜精品视频节目| 漂亮人妻口爆久久精品| 欧美精品激情在线不卡| 四虎国产精品国产精品国产精品| 天天早上头和脸出汗是怎么办| 熟女人妻精品视频一区| 国产精品igao为爱寻找激情| 亚洲成人,国产精品| 嗯~嗯~啊啊啊~高潮了软件| 蜜桃tv一区二区三区| caopeng97在线观看视频| 女人的天堂av在线网| 一区二区三区资源视频| 91麻豆精品国产在线| 一区二区在线观看视频观看| 激情久久在线免费观看视频| 日韩人妻中文字幕二区 | 天天想要天天操天天干| 神马不卡视频在线视频| 福利美女视频在线观看| 91精品久久久久久久99蜜月| 日本高清激情乱一区二区三区| 老熟女 露脸 嗷嗷叫| 福利一二三在线视频观看| 欧美精品999不卡| 成年人免费黄色av| 免费观看在线中文字幕视频| 日韩成人精品久久久免费看| 天堂av在线最新地址| 国产精品久久人人添| 2020国产成人精品视频| 色视频在线播放免费观看| 亚洲人妻系列在线视频| 国产主播诱惑毛片av| 中国特黄色性生活片| 大乳丰满人妻中文字幕韩国hd| aaaa级少妇高潮在线观看| 国产视频成人一区二区| xxoo福利视频导航| 人妻系列中文字幕大乳丰满人妻| 国产精品午夜无码AV体验区| 亚洲成人偷拍自拍在线| 国产伦理二区三区在干嘛呢| 亚洲成人av在线一区二区| 日本香港韩国三级黄色| 在宿舍强奷两个清纯校花| 国产午夜在线播放视频| 国产精品中文字幕丝袜| 18福利视频在线观看| 国产视频成人自拍蝌蚪视频| 天天日天天玩天天摸| 东京热日韩av影片| 青青青在线视频观看97| 国产欧美福利在线观看| 国产探花自拍亚洲av| 精品国模一区二区三区欧美| 在线观看免费啪啪啪| aaaa级少妇高潮在线观看| 欧美亚洲国产一区二区| 日本东京热视频欧美视频| 夜夜爽夜夜操夜夜爱| 日本人妻少妇xxxxxxx| 自拍偷拍 国产激情| av资源中文字幕在线观看| 亚洲一区视频中文字幕在线播放 | 亚洲综合天堂av网站在线观看| 天天操天天搞天天操| 精品国产久久久久午夜精品av| 360偷拍蜜桃臀69式| 日韩国产欧美一区二区三区粉嫩 | 男人的天堂aⅴ在线| 九九六视频,这里只有精品| 黄片视频免费观看视频| 老司机伊人99久久精品| 亚洲|久久久久久一二三区丝袜 | 久久午夜免费鲁丝片| 亚洲美女午夜激情视频在线观看| 日本少妇精品免费视频| 成人av在线视频免费| 九热精品视频在线观看| 欧美精品乱码99久久蜜桃免费| 亚洲同性同志一二三专区| 成人精品影视一区二区| 操人妻人妻天天爽天天偷| 五十岁熟妇高潮喷水| 性感美女人妻久久久| 五月天男人的天堂中文字幕 | 三级欧美日韩一区二区三区| 天天看片天天摸天天操| 中文字幕 中文字幕 亚洲| 日韩免费黄色片在线观看| 啊~插得好快别揉我胸了视频| 亚洲三级综合在线观看| 久久久亚洲熟女一区二区| 欧美精品熟妇免费在线| 每日更新日韩欧美在线| 亚洲国产日韩欧美一区二区三区,| 核xp工厂精品久久亚洲| 一区二区三区国产在线成人av| 青青操久久综合激情| 蜜桃臀少妇白色紧身裤细高跟| 一区二区欧美 国产日韩| 日本特级黄片免费观看| 成人免费视频现网站99在线观看| 91九色人妻在线播放| 日韩av水蜜桃一区二区三区| 中文字幕在线观看亚洲情色| 精品日本少妇久久久| 全彩漫画口工18禁| 国产在线观看av一区| 国产igao激情在线视频入口| 免费啪啪啪网站在线观看| www一区二区91| 人妻人妻在线视频网站| 美女网站福利在线观看| 国产熟女五十路一区二区三区| 国产女人18毛片水真多精选| 538欧美在线观看一区二区三区| 欧美人与动欧交视频| 制服丝袜中文字幕熟女人妻| 情趣视频在线观看91| 日本人妻熟妇丰满成熟HD系列| 国产福利三级在线观看| 91九色尤物无套内射| 东京热男人的天堂视频| 亚成区一区二区人妻熟女| 裸露视频免费在线观看| 亚洲18片综合国产av| 成年男女免费视频网站无毒| 人妻激情综合久久久久蜜桃| avtt中文字幕手机版| 亚洲午夜国产末满十八岁勿进网站| 顶级欧美色妇xxxx| 91中文字幕视频网站| 亚洲熟女一区二区三区250p| 欧美日韩久久丝袜在线| 新香蕉视频香蕉视频2| 国产91九色视频在线观看| 亚洲综合成人精品成人精品| 伊人网在线欧美日韩在线| 欧美一级aaaaaaa片| 97香蕉久久国产超碰| 99国产精品国产精品毛片19| 婷婷一区二区三区五月丁| 亚洲无码专区中文字幕专区| 天天插天天干天天狠| 亚洲av毛片一区二区三区网| 中文字幕久久久国产| 亚洲精品国品乱码久久久久| 91麻豆精品国产在线| 91精品资源在线观看| 亚洲综合首页综合在线观看 | 深夜福利免费观看在线看| 天天综合久久无人区| 国语精品视频自产自拍| 琪琪日本福利伦理视频| 国产最新av在线免费观看| 欧美日本在线免费视频| 91精品麻豆91夜夜骚| 日本老女人日比视频| 人妻系列级片在线观看视频| av在线免费在线观看| 最新日韩av电影在线播放| 一区二区三区国产在线成人av| 熟妇人妻av无码中文字幕| 日本一区二区三区的资源| xxxx69在线观看视频| 天堂网免费在线电影| 欧美大胆a级视频秒播| 国产 少妇 一区二区| 亚洲欧美日韩电影一区| 熟妇人妻av无码中文字幕| 国产高清在线观看av| 一区二区三区四区久久久久韩日| 婷婷色九月综合激情丁香| 顶级欧美色妇4khd| 黑鸡巴肏少妇逼视频| 69精品人妻久久久久久久久久久| 午夜精品一区二区三区不卡顿| 亚洲国产日韩a在线欧美| 女生抠逼自慰啊啊啊啊啊啊啊下载 | 午夜3p福利视频合集| 国产精品免费看一区二区三区| 亚洲 自拍 激情 另类| 国产一区两区三区福利小视频| 最近在线中文字幕免费| 在线成人教育平台排名| 快进来插我的逼嗯啊视频| 亚洲在线观看中文字幕av| 黄色大片一级老太太操逼| 日韩av电影中文在线免费观看| 国产欧美福利在线观看| 成年人免费福利在线| 69av精品国产探花| 在线成人教育平台排名| 国产在线小视频一区二区| 激情九月天在线视频| 免费绝清毛片a在线播放| 成人人妻h在线观看| 亚洲蜜桃久久久久久| 亚洲中文字幕最新地址| 成人精品影视一区二区| 裸日本资源在线午夜| 成人午夜麻豆大胆视频| 久久久久久久久久久久久国产| 亚洲国产日韩欧美一区二区三区,| 丝袜美腿日韩av一区| 亚洲黄色成人一级片| 精品国产av虐杀两警花| 日韩少妇免费在线播放| 国际精品熟女一区二区| 午夜精品视频免费观看| 中文字幕 首页 人妻| 国产一区两区三区福利小视频 | 欧美日韩亚洲国产视频二区| 好看的日本中文字幕在线观看二区| 亚洲免费在线不卡视频| 亚洲日本欧美韩国另类综合| 4438x亚洲最大的成人| 日本少妇熟女乱码一区二区| 蜜乳视频一区二区三区| alisontyler和黑人| 午夜呻吟亚洲精品中文字幕在上面| 日本成人福利电影网| 午夜精品久久久久久久久久蜜桃| 伊人免费观看视频一| 小妹妹爱大棒棒免费观看视频 | 97视频人人爱麻豆| 97人妻av人人澡人人爽| 黄色片黄色片黄色片黄色片黄色| 亚洲一区在线视频观看地址| 99久久99九九九99九| 中文字幕日韩首页欧美在线激情 | 日韩三级黄色大片在线观看| ysl蜜桃色7425| 女同性恋av在线播放| 台湾18禁久久久久久久激情视频| 中文字幕在线免费观看成人| 91进入蜜桃臀在线播放| 91久久精品美女高潮喷水白浆| 亚洲第一成年偷拍视频| 亚洲成人激情在线综合| 国产91免费在线观看| 最新福利二区三区视频| 一区二区三区不卡免费视频网站 | 福利视频免费在线播放| 国产中文亚洲熟女日韩| 中文字幕国产一区在线视频| 美女av色播在线播放| 农村大炕有肉大屁股熟妇| 亚洲熟女少妇中文字幕系列| 午夜福利片无码10000| 青青操久久综合激情| 久久久久久a女人处女| 大秀成年人国产精品视频| 最新日韩av电影在线播放| 国产精品久久久99| 日韩男女视频网站在线观看| 天天早上头和脸出汗是怎么办| 99久久精品视频16| avjpm亚洲伊人久久| 国产成人深夜福利短视频99| 亚洲综合成人精品成人精品| 全国熟妇精品一区二区免费视频| 免费24小时人妻视频| av网页免费在线观看| 91精品91久久久久| 成年男女免费视频网站无毒| 日韩精品视频一区二区三区在线| 天堂av国产av伦理av| 精品人妻在线激情视频| 91 精品视频在线看| 国内精品一区二区2021在线| 亚洲综合天堂av网站在线观看| 麻豆国产91制片厂| 黄在线看片免费人成视频| 一区二区三区午夜福利在线| 黄片视频免费观看视频| 亚洲人妻系列在线视频| 国产女主播在线观看一区| 二十四小时日本高清在线观看| 男女爱爱好爽视频免费看| 人妻系列在线免费视频| 免费高清av一区二区| 天天操天天舔天天爽| 丰满少妇高潮喷水视频| 麻豆白洁少妇在线播放| 夜色福利视频免费观看| 男女啪啪啪啪91av日韩| 欧美成人一二三在线网| 2021国产剧情麻豆| 一看就是假奶的av| 欧美第一激情综合网欧美激情| 裸日本资源在线午夜| 男人的天堂aⅴ在线| 欧美日韩黄片免费在线观看| 黄片操操操操操操c| 九九视频在线观看全部| 黑人大巨屌操美女逼| 国产成人av在线你懂得| 99re这里是国产精品首页| 二十四小时日本高清在线观看 | 一区二区三区观看在线| 免费啪啪啪网站在线观看| 二十四小时日本高清在线观看| 亚洲理论在线a中文字幕97| 亚洲国产综合久久精品| 深夜福利免费观看在线看| 中文字幕熟女乱一区二区| 青青青在线视频免费播放| 天堂在线中文字幕av| 蜜乳av一区二区三区免费观看| 色哟哟亚洲乱码国产乱码精品精| 成人精品影视一区二区| 久久99国产中文丝袜| 国产高清自拍偷拍在线| 99色在线观看免费观看| 女女抠逼白虎白丝袜| 女同大尺度视频网站在线观看| 夜夜躁av麻豆男| 精品欧美乱码久久久| 亚洲中文字幕在线av| 亚洲AV无码久久精品国产一区老| 黄片视频免费观看视频| 欧美男男在线观看视频网站| 亚洲男人天堂最新网址大全| 亚洲精品综合欧美精品综合| 91 精品视频在线看| 国产一区二区三区四区精| 精产国品一二三产品区别91| 国色天香一二三期区别大象| 成人免费视频现网站99在线观看 | 午夜亚洲国产精品中字| jandara在线观看| 2021国产剧情麻豆| 中文字幕熟女人妻一区| 亚洲同性同志一二三专区| 在线观看黄页网站视频网站| 天天摸天天舔天天操天天日| 日韩免费黄色片在线观看| 国产自拍偷拍视频在线免费观看 | 五十岁熟女高潮喷水| 欧美成人短视频在线播放| 天天干夜夜操夜夜骑| 成年人黄色日本视频| 欧美精品999不卡| 精品国产久久久久午夜精品av| 在线中文字幕人妻av| 十八禁黄色免费污污污亚洲| 中日韩又粗又硬又大精品| 美女网站视频久久精品| 熟女一区二区视频在线 | 久久久久久高清一区| 操操操操操操操操操网| 国语精品视频自产自拍| 中文字幕精品人妻久久久久| 国产毛片特级Av片| 欧美vs亚洲vs日韩| 久久精品国产亚洲av热软件| 日本成人福利电影网| 核xp工厂精品久久亚洲| 国产91黑丝小视频在线观看| 亚洲av在线免费播放| 美国伦理片午夜理论片| 国产自拍偷拍在线精品| 97人妻在线视频自拍| 外国美女舔男人坤坤| 日韩黄色在线观看网站上 | 亚洲激情视频在线观看免费| 免费在线观看黄色小网站| 久久99精品热在线观看| 亚洲国产电影的一区| 亚洲av中文无码网站| 色哟哟亚洲乱码国产乱码精品精| 亚洲熟女乱一区二区精品成人| 国产精品网站的黄色| 亚洲一区二区三区国产精品电影| avjpm亚洲伊人久久| 中文字幕福利视频在线一区| 亚洲av激情综合网| 男人的天堂av中文字幕| avtt中文字幕手机版| 欧美精品激情在线不卡| av在线免费在线观看| 欧美视频亚洲视频在线| 在线 制服 中文字幕 日韩| 欧美男女一区二区三区| 熟女阿高潮合集一区二区| 国产精品igao为爱寻找激情| 奇米网首页神马久久| 免费24小时人妻视频| 538欧美在线观看一区二区三区 | 久久久视频在线播放| 五十岁熟女高潮喷水| 最新日韩中文字幕啪啪啪| 全国熟妇精品一区二区免费视频 | 国产大桥未久一区二区| 最新久久这里只有精品| 日韩激情亚洲国产欧美另类激情| 99久久人人爽亚洲精品美女| 国产青青青青草免费在线视频| 亚洲高清免费在线观看视频| 熟女一区二区视频在线| 亚洲免费午夜污福利| 午夜久久久久久av五月| 欧美日韩亚洲tv不卡久久| av人摸人人人澡人人超碰小说| av网页免费在线观看| 日本高清 中文字幕| 制服丝袜 中文字幕 日韩| 日本成人福利电影网| 在线免费观看视频18| 日本欧美国产在线一区| 在线观看中文字幕精品av| 欧美日本在线免费视频| 亚洲国产精品一区51动漫| 天天干夜夜操夜夜骑| 18在线观看免费观看| 精品不卡一区二区三区| 日本午夜福利免费在线播放| 亚洲妹妹我爱你在线观看| 日本亚洲午夜福利一区二区三区| 成人午夜av电影网| av中文字幕国产精品| 亚洲永远av在线播放| 黄很色很在线免费视频网站| 伊人免费观看视频一| 国产成人情侣av在线| 国产 亚洲 欧美 自拍| 精品人妻在线激情视频| 最新日韩av电影在线播放| 亚洲精品激情视频在线观看| 午夜一区二区三区视频在线观看| 精品欧美乱码久久久| 人妻系列级片在线观看视频| 最近最新欧美日韩精品| 国产精品久久久久久成人久| 高潮喷水一区二区三区| 免费24小时人妻视频| 欧美男女一区二区三区| 色狠狠色综合久久久绯色| 伊人久久综合国产精品| 夜夜爽夜夜操夜夜爱| 五十岁熟女高潮喷水| 天堂网免费在线电影| 五月婷婷激情视频网| 男人电影天堂在线观看| 中文字幕人妻一区色偷偷久久| 午夜免费福利老司机| 亚洲美女午夜激情视频在线观看| 亚洲妹妹我爱你在线观看| 青青在线免费手机播放视频| 中文字幕在线观看av观看| 日本黄色一级电影网址| 天天干夜夜撸天天操| 亚洲av日韩久久网站| jiee日本美女视频网站| 男人用大鸡巴狂操女人肉穴| 大尺度av毛片在线网址| 麻豆国产91制片厂| 91久久久久久最新网站| 日本高清在线观看不卡视频| 可在线免费观看av| 亚av一二三在线观看| 不卡一区二区视频在线| 亚洲一区二区三区国产精品电影| www,日韩av,com| 美女网站视频久久精品| 亚洲欧美日韩中文在线观看| 二十四小时日本高清在线观看| 亚洲一区在线视频观看地址| 91精品91久久久久| 国产人妻777人伦精品hd超碰| 国产一区二区三区四区精| 中文字幕丰满子伦无码专区| 成人资源中文在线观看| 亚洲色视频在线播放网站| 老鸭窝在线毛片观看免费播放| 久久视频 在线播放| 久久久久久久久久久久久国产| 久草视频在线看免费| 女同性恋av在线播放| xxnxx国产美女| 天天色 天天操 天天好逼| 自拍偷拍 国产激情| 老色鬼精品视频在线观看播放| 狠狠操av一区二区三区| 最新久久这里只有精品| yy4080黄色片| 宅男噜噜噜66国产在线观看| 亚洲av中文免费在线| 亚洲男人的天堂最新网址| iga肾三级算严重吗| 午夜美女福利视频在线| 激情久久在线免费观看视频| 在线免费观看欧美小视频 | 美女一区二区四区六区八区| 97精品久久久久久无码人妻| 美女欧美视频在线观看免费| 亚洲成人激情在线综合| 国产高清视频www夜色资源| 亚洲综合熟女乱中文| 黑川堇人妻88av| 天天操天天日天天插天天舔| 最近最新最好看的中文字幕| 国语对白性爱三级片免费看| 在线视频自拍第三页| 亚洲午夜精品一级毛片app| 久久国产半精品99精品国产| 欧美成人性生活视频播放| 日韩A级毛片免费视频| 日韩三级黄色大片在线观看| 9久re热视频在线精品| 91大神在线免费观看视频| 亚洲国产日韩a在线欧美| 亚洲综合成人精品成人精品| 91 精品视频在线看| 麻豆出品视频在线观看| 中文字幕欧美一区二区视频| 中国特黄色性生活片| 青娱乐免费最新视频| 夏目彩春av在线看| 中文字幕亚洲乱码精品无限| 91日本精产品一区二区三区| 亚洲无人区乱码中文字幕一区| 一二区二区不卡视频| 成年人黄色日本视频| 免费在线观看视频啪啪| 91性高湖久久久久久久久久| 欧美vs亚洲vs日韩| 亚洲精品久久久人妻| 久久午夜免费鲁丝片| 亚洲唯美激情综合四射| 亚洲欧美不卡专业视频| 国产中文亚洲熟女日韩| 婷婷一区二区三区五月丁| 东京热男人的天堂视频| 又粗又长又硬又黄又爽| 人妻少妇精品二三区| 天天早上头和脸出汗是怎么办| 天天操天天舔天天射天天日天天干| 亚洲国产综合久久精品| 久久精品久久久久观看99水蜜桃| 人妻系列中文字幕大乳丰满人妻 | 色视频在线播放免费观看| 插鸡视频免费网站在线播放| 日本美女爱爱视频网站| 91超碰国产在线观看| 亚洲一区二区精品在线播放| 色欲AV亚洲AV无码精品| 美女精品久久久久久久久| 人人妻人人澡人人爽97| 国产av嗯嗯啊啊av| 欧美一区日韩二区三区四区| 丝袜美女诱惑佐佐三上| 亚洲精品国产99999| 91精品久久久久久久久99蜜臀| 60路70路日本熟妇| 9999久久久久老熟妇二区| 91精品资源在线观看| 国产精品性感美女视频| 中文字幕熟女乱一区二区| 91九色国产在线视频| 天天干夜夜爽狠狠操| 日本福利视频网站导航| 亚洲永远av在线播放| 亚洲av三级电影在线观看| 国产精品蝌蚪自拍视频| 岛国av成人午夜高清| av激情四射五月婷婷| 国产资源网站在线播放| 久草视频在线看免费| 日本熟女0930视频| 看女人大BB群伦交| 激情久久在线免费观看视频| 亚洲色图日韩在线视频观看| 国产不卡免费在线观看| 伊人网在线免费观看| 国产igao激情在线视频入口| 亚洲美女露隐私av一区二区精品| 欧美日韩成人高清中文网| 老司机伊人99久久精品| 人妻超清中文字幕在线乱码| 熟女俱乐部jukujoclub| 漂亮人妻口爆久久精品| 亚洲熟女一区二区六区| 丰满少妇高潮喷水视频| 亭亭五月天在线观看| 天海翼亚洲一区在线观看| 少妇熟女天堂网av| 顶级欧美色妇xxxx| 国产精品久久久久精品三级18| 国产高清视频www夜色资源| 岳母的诱惑电影在线观看| 欧美亚洲愉拍一区二区三区| 日本黄页在线观看视频| 亚洲男人的天堂最新网址| 欧洲成熟女人色惰片| 国产精品蝌蚪自拍视频| 福利美女视频在线观看| 国产毛片特级Av片| 老司机免费视频福利0| 黄在线看片免费人成视频| 日本欧美视频在线免费| 黄片操操操操操操c| a级黄片免费观看| 亚洲唯美激情综合四射| 十八禁黄色免费污污污亚洲| 最近最新最好看的中文字幕| 色老头一区二区三区四区五区| 亚州av嫩草av极品在线观看| 天天搞天天操天天干| 欧美成人少妇人妻精品| 亚洲熟女人妻自拍在线视频| 国产伦理二区三区在干嘛呢| 操死你美女在线视频| 99精品视频在线在线观看| 91九色国产在线视频| 日本高清久久人人爽| 国产精品无码无卡免费观| av一区二区三区蜜桃| 久久久亚洲熟女一区二区| 亚洲国产精品自产拍在线观看| 91精品久久久久久久99蜜月| 亚洲欧美激情国产综合久久久| 一区二区三区国产精华液区别大吗| 中文字幕日韩首页欧美在线激情| 欧美日本亚欧在线观看| 50熟妇一区二区三区| 亚洲在线免费观看18| 亚洲无人区乱码中文字幕一区| 久久久久性感美女偷拍视频| 3344永久在线观看视频下载| 国产成人在线观看hd| 亚洲精品1卡2卡3卡| iga肾三级算严重吗| 美女激情久久久久久久| 丰满人妻被猛烈进入中文字幕| 两个奶被揉得又硬又翘怎么回事 | 天天操天天干天天谢| 国产精品 亚洲欧美 自拍偷拍| av 资源在线播放| 伊人久久综合国产精品| 男女啪啪啪啪91av日韩| 熟女一区二区视频在线| 老司国产精品视频免费观看| 国产三级自拍视频在线观看网站| 成人黄色录像在线观看| 福利在线国产小视频| 欧美精品999不卡| 日韩一级视频一区二区三区| 日韩激情亚洲国产欧美另类激情| 日韩A级毛片免费视频| 国产精品久久人人添| 91九色人妻在线播放| 亚洲熟女一区二区六区| 日韩成人在线电影首页| 欧美日韩综合精品无人区 | 中文字幕 一区二区在线观看| 高清国产美女a一级毛片| 亚洲精品乱码久久久久app| 操烂你的骚逼天天欧美| 天天插天天透天天爽| 亚洲av手机免费在线| 欧美黄色一区二区三区视频| 国产美女高潮精品视频| 大奶熟妇激情操逼逼| 亚洲第一区av中文字幕| 亚洲综合首页综合在线观看| 在线中文字幕人妻av| 亚洲成人动漫av在线| 女人的天堂 av在线| 日本福利视频网站导航| 亚洲一区二区精品在线播放| 免费中文字幕a级激情| 日本一道中文字幕99| 狠狠操狠狠操狠狠插| 中文人妻av一区二区三区| 黄色片黄色片黄色片黄色片黄色| 天天想要天天操天天干| 午夜五十路久久福利| 中文字幕久久久国产| 亚洲av 综合av| 国产精品国产三级在线高清观看| 一区二区三区五区六区| 青青操天堂在线观看视频| 蜜臀一区二区日韩美女少妇视频| 操人妻人妻天天爽天天偷| 秋霞成人午夜鲁丝一区二区三区| 午夜久久人妻一级内射av网址 | 最近中文字幕免费视频一| 99精品视频在线在线观看| 精品视频一区二区三区◇| 日韩在线 中文字幕| 2020精品视频在线| 黑人和日本人av一区二区| 极品少妇高潮喷水日出白浆| 成年男女免费视频网站无毒| 日本老女人日比视频| 新香蕉视频香蕉视频2| 天堂网免费在线电影| 小妹妹爱大棒棒免费观看视频| 日本黄色一级电影网址| 中文字幕一区二区三区久久久| 成人18禁高潮片免费日本| 天天天天天天天天干夜夜| 亚洲欧美日韩中文在线观看| 夜夜骚av一二三区| 久久99久久99久久97的人| 激情久久在线免费观看视频| 精品美女洗澡一区二区| 伦理在线观看未删减中文字幕| 欧美日韩成人高清中文网| 国产视频成人自拍蝌蚪视频 | 日韩黄色在线观看网站上| 日韩久久不卡免费视频| 久99久视频免费观看中文字幕| 久久99精品久久久久久三级| 亚洲国产美女主播在线观看| 东京热男人的天堂视频| 天天透天天舔天天操| 亚洲一区二区在线激情| 国产av高清二区三区| 亚洲熟女人妻自拍在线视频 | 日本在线免费观看国产精品| 欧美一级日韩一级亚洲一级va| 91九色91在线视频| 福利视频免费在线播放| 大屁股熟女一区二区视频| 国产成人在线观看视频播放| 美女欧美视频在线观看免费| 第一福利视频在线观看| 国产自拍偷拍在线精品| 国产女人18毛片水真多精选| 加勒比不卡在线视频| 欧美视频亚洲视频在线| 国产在线观看av一区| 在线看日韩av不卡| 亚洲av网站一区二区三区| 大成色亚洲一二三区| 亭亭五月天在线观看| 天堂在线中文字幕av| v天堂国产精品久久| 顶级欧美色妇4khd| 日本午夜福利免费在线播放| 琪琪日本福利伦理视频| 黄色片免费网站在线| 69精品人妻久久久久久久久久久| av在线播放观看h| 精品一区二区三区喷水内射高潮| 91香蕉国产亚洲一二三区| 中文字幕在线观看av观看| 97精品国产91久久久| 成年人黄色日本视频| 91超碰国产在线观看| 精品一区二区三区免费毛片W| 人妻少妇视频系列视频在线| 天天操天天日天天插天天舔| 日本一区二区三区区别| 69视频在线精品国自产拍| 五月婷婷激情视频网| 久久久精品人妻无码专区不卡 | 亚洲成人欧洲成人在线| 在线免费观看欧美小视频| 亚洲永远av在线播放| 欧美日本国产一区二区| 亚洲精品综合欧美精品综合| 九一精品人妻一区二区三区| 亚洲妹妹我爱你在线观看| 久久热在线免费观看| 可以免费观看日韩av| 熟女国内精品一区二区三区| 亚洲国产日韩a在线欧美| 青青青国产精品视频| 中字幕人妻熟女人妻a62v网| 北野中文字幕一区二区| 欧美成人性生活视频播放| 人人妻人人澡人人爽97| 不卡视频在线 欧美日韩| 北野中文字幕一区二区| 日韩av熟妇在线观看| 日本国产亚洲欧美色综合| 亚洲在线观看中文字幕av| 精品高潮呻吟久久av| 裸日本资源在线午夜| 每日更新日韩欧美在线| 午夜精品久久久久久久精品乱码| 蜜乳av一区二区三区免费观看| 欧美精品999不卡| 国产黄色主播网址大全在线播放| 松本菜奈实最新av在线| 360偷拍蜜桃臀69式| 日韩人妻中文字幕二区| 插鸡视频免费网站在线播放| 全彩漫画口工18禁| 麻豆国产精品777777在| 鸡巴在里面福利视频在线观看| 丰满少妇人妻一区二区三区蜜桃| 久久亚洲国产成人精品麻豆| 无人区一码二码三码区别在哪| 国产原创一区二区三区在线播放| 猫咪亚洲中文在线中文字幕| 91超精品碰国产在线观看| 夜夜骚av一二三区| 91精品国产欧美在线| 天天操天天干天天舔天天| 大秀成年人国产精品视频 | 手机看片福利一区二区三区四区 | 蜜臀一区二区日韩美女少妇视频| 69精品互换人妻4p| 欧美最新一区二区三区| 亚洲国产日韩欧美一区二区三区,| 日韩成人精品久久久免费看| 亚洲一区二区三区四区入口| 国产又粗又长又大视频| 中文字幕欧美一区二区视频| 99精品久久精品一区二区| 中文字幕在线字幕乱码怎么设置| 亚洲a级视频在线播放| 久久99国产中文丝袜| 日本免费人爱做视频在线观看不卡| 亚洲欧美综合另类最新| 中文字幕人妻一区二区视频系列 | 免费观看在线中文字幕视频| 在线免费观看a视频免费| 欧美大胆a级视频秒播| 手机看片福利一区二区三区四区| 亚洲熟女少妇中文字幕系列| 亚洲国产精品青青草| 九九九九九久久久国产| 超peng视频在线免费播放97| 妈妈的朋友中字在线免费观看| 亚洲宅男噜噜噜66在线观看| 欧美久久一区二区伊人| 午夜精品久久久久久久精品乱码| 伊人网在线免费观看| 天天日天天亲天天操| 久久国产半精品99精品国产| 成人免费视频现网站99在线观看| 免费看超污视频在线观看| 成人大片男人的天堂| 九色91操最新在线观看网址| 日韩精品欧美一区二区| 中文字幕在线观看av观看| 高潮喷水在线视频观看| 欧美日韩综合精品无人区| 午夜国产成人精品视频观看| 亚洲人成大片在线观看| 天天操天天干加勒比久久| 情趣视频在线观看91| 欧美最新一区二区三区| 亚洲美女色www色| 国产91九色视频在线观看| 在线观看网站伊人网| 欧美成人屋影院在线视频观看| 青青操久久综合激情| 亚洲熟女一区二区六区| 9999久久久久老熟妇二区| 女生裸体视频免费网站| 天天操天天舔天天射天天日天天干| 玖玖资源站在线观看亚洲| 丝袜美腿日韩av一区| 顶级欧美色妇4khd| 小妹妹爱大棒棒免费观看视频| 中文字幕 中文字幕 亚洲| 日本特级黄片免费观看| 婷婷一区二区三区五月丁| 国产做A爱免费视频在线观看| 日韩精品视频一区二区三区在线| 91精品国产欧美在线| 中文字幕av人妻一区二区三区| 亚洲制服丝袜美腿在线| 国语精品视频自产自拍| 欧美最新一区二区三区| 国产夫妻视频在线观看免费| 五月天天堂视频在线| 东京热日韩av在线| 欧美一级aaaaaaa片| 4438全国成人免费视频| 九一精品人妻一区二区三区| 中文字幕av特黄毛片| 天堂av在线最新地址| 亚洲美女黄色福利视频网站大全| 熟女国内精品一区二区三区 | 2020精品视频在线| 久久精品国产亚洲av清纯| 日本福利视频网站导航| 乱子伦国产一区二区三区| 麻豆白洁少妇在线播放| 亚洲天堂av最新在线| 亚洲熟女一区二区三区250p| 午夜夫妻性生活视频| 精品国产人伦一区二区三区| 亚洲国产美女主播在线观看| 国产精品美女免费视频观看| 亚av一二三在线观看| 日韩在线 中文字幕| 99久久精品视频16| 熟女一区二区三区综合| 中文字幕人妻一区色偷偷久久 | ysl蜜桃色7425| 天天操天天舔天天做| www国产亚洲精品久久久| 亚洲av激情综合网| 人妻熟女 亚洲 一页二页| 快使劲弄我视频在线播放| 精品免费一区二区三区四区视频| 天天干天天日天天弄| 久久99嫩草99久久精品| 99精品久久一区二区| tobu8日本高清| 美利坚合众国av天堂| 一级毛片特级毛片免费的| 亚洲熟女在线免费观看| 国产肥胖熟女又色又爽免费视频| 天天操天天射天天操天天日| 亚洲男人天堂最新网址大全| 欧美肥妇久久久久久| 猫咪亚洲中文在线中文字幕| 99免费观看在线视频| 制服丝袜 中文字幕 日韩| 亚洲国产日韩欧美一区二区三区,| 日韩成人精品久久久免费看| 国产igao激情在线视频入口| 中文字幕精品人妻久久久久| 91福利高清在线播放| 欧美视频亚洲视频在线| av一区二区三区四区五区在线| 亚洲欧美精品日韩偷拍| ysl蜜桃色7425| 色丁香久久激情综合网| 亚洲同性同志一二三专区| 亚洲人成大片在线观看| 伊人网在线观看 视频一区| 4日日夜夜精品视频免费| 亚洲成人,国产精品| 人妻激情综合久久久久蜜桃| 欧美在线视频不卡一区| 黄色av网址在线播放| 欧美久久蜜臀蜜桃资源吧| 亚洲图片另类综合小说| 国产女人18毛片水真多精选| 伦理在线观看未删减中文字幕| 黄片操操操操操操c| 东北老女人熟女啪啪视频| 欧美视频免费观看777| 91亚洲精品久久蜜桃| 美女把腿张开给男的捅| 亚洲另类欧美综合久久| 女女抠逼白虎白丝袜| 丰满人妻被猛烈进入中文字幕| 免费在线观看视频啪啪| 国产91九色视频在线观看| 开心五月综合激情婷婷| 538欧美在线观看一区二区三区| 操烂你的骚逼天天欧美| 亚洲色图日韩在线视频观看| 二十四小时日本高清在线观看| 国产女主播在线观看一区| 琪琪日本福利伦理视频| 在线有码人妻自拍视频| 国产激情一区二区视频| 久久久亚洲综合国产精品| 国产一级一国产一级毛片| 日本人妻少妇xxxxxxx| 99国产精品国产精品毛片19| 美女张开腿给男人桶爽的软件| 欧美精品激情在线不卡| 亚州av嫩草av极品在线观看| 久久99久久99久久97的人| 日韩av电影中文在线免费观看| 操烂你的骚逼天天欧美| 青娱乐这里只有精品| 亚洲国产精品自产拍在线观看| 一区二区三区资源视频| 黑人巨大精品一区二区在线 | 国产黑色丝袜 在线日韩欧美| 韩日一级人添人人澡人人妻精品| 五月天色婷婷狠狠爱| 青青青青青爽视频在线| 日韩激情文学在线视频| 999久久久人妻精品一区| 在线观看中文字幕少妇av| av在线男人的天堂亚洲| 日本高清在线观看不卡视频| 18岁禁一二三区免费体验| 夜色17s精品人妻熟女av| 亚洲欧美另类丝袜另类自拍| 天天插天天干天天狠| 天天摸天天舔天天操天天日| 日韩最近中文在线观看| 国产av啊啊啊啊啊啊啊| 大秀成年人国产精品视频 | 久久一级片三上悠亚| 欧美巨大另类极品video| 亚洲欧美国产人成在线| 亚洲国产日韩欧美一区二区三区, 精久久久久久久久久久久 | 97人妻人人揉人人躁人人夜夜爽| av福利免费体验观看| 最近日韩免费在线观看| 日本少妇三级交换做爰做| 黑川堇人妻88av| 中文字幕免费啪啪啪| 国产精品午夜无码AV体验区| 视频免费在线观看网站| 女生裸体视频免费网站| 91九色人妻在线播放| 日韩一级视频一区二区三区| 东京热日韩av影片| 91精品国产人妻麻豆| 成人免费电影二区三区| 欧美色区国产日韩亚洲区| 黄色大片一级老太太操逼| 国产大桥未久一区二区| 男人电影天堂在线观看| 91进入蜜桃臀在线播放| 上床啪啪啪免费视频| 乱子伦国产一区二区三区| 乌克兰美女操逼高清内射视频| 国产女主播在线观看一区| 国产黑色丝袜 在线日韩欧美| 亚洲 偷拍 自拍 欧美| 午夜久久久久久av五月| 午夜久久久久久av五月| 妈妈的朋友2中文字幕在线| 东北老女人熟女啪啪视频| 亚洲午夜熟女在线观看| 日韩人妻一区二区三区在线观看| 亚洲国产日韩精品在线| 亚洲另类激情视频在线看| 国产 少妇 一区二区| 男生和女生羞羞91在线看| 国产精品久久久久精品三级18| 最新日韩中文字幕啪啪啪| 在线观看黄页网站视频网站| 天天操天天日天天碰| 国产精品性感美女视频| 日本高清激情乱一区二区三区| 91色乱一区二区三区| 裸露视频免费在线观看| 伊人精品久久一区二区| 日韩免费黄色片在线观看| 在线观看网站伊人网| 日韩久久九九精品视频| 日韩一级欧美一级片| 美女黄色啊啊啊啊视频| 一区二区三区四区影片| 亚洲熟女人妻自拍在线视频| 户外露出视频在线观看| 手机视频在线观看一区| 都市激情校园春色 亚洲| 午夜五十路久久福利| 国产精品成人免费电影| 裸日本资源在线午夜| 污网址在线观看视频| 欧美色区国产日韩亚洲区| 黄在线看片免费人成视频| 一区二区三区免费版在线| 久久一级片三上悠亚| 久久内射天天玩天天懂色| 91精品国产成人久久久久久| 97精品视频,全部免费| 久久99精品热在线观看| 色老头一区二区三区四区五区| 天天爽天天操天天插| 日本在线免费观看国产精品| 午夜一区二区三区视频在线观看| 国产 少妇 一区二区| 69xx精品久久久久| 午夜国产一区二区三区| 开心五月综合激情婷婷| 久久99久久99久久97的人| 9999久久久久老熟妇二区| 婷婷色九月综合激情丁香| 一区二区三区不卡免费视频网站| 综合激情网,激情五月| 91精品综合久久久久久五月天| 日本成年视频在线免费观看| 中文字幕麻绳捆绑的人妻| 日本老熟妇av老熟妇| 中文人妻av一区二区三区| 亚洲唯美激情综合四射| 亚洲色大WWW永久网站| 亚洲av综合av一去二区三区| 欧美成人短视频在线播放| 天天干天天操天天要| 国产午夜羞羞一区二区三区| 亚洲乱熟女一区二区三区影片| 国产精品亚洲精品亚洲| 白白色在线免费视频发布视频| 川上优所有中文字幕在线| 快使劲弄我视频在线播放| 男人电影天堂在线观看| 亚洲成人 国产精品| 亚洲欧美一级特黄大片| 午夜情色一区二区三区| 欧美性受黑人猛交裸体视频| 亚洲最强的25个城市| 欧美在线观看视频欧美| 爱搞视频在线观看视频91| 日本人妻少妇xxxxxxx| 高清国产美女a一级毛片| 在线 制服 中文字幕 日韩| 黑人侵犯人妻森泽佳奈| 九九热在线精品播放| 高潮喷水一区二区三区| 日本老女人日比视频| 在线观看免费啪啪啪| 91色哟哟视频在线观看| 99re这里是国产精品首页| yellow在线亚洲精品一区| 九热精品视频在线观看| 精品精品精品精品精品污污污污| 欧美一级日韩一级亚洲一级va| 97视频人人爱麻豆| 国产美女主播av在线| 亚洲中文字幕在线视频观看二区| 成人午夜麻豆大胆视频| 川上优所有中文字幕在线| 精产国品一二三产品区别97| 亚洲成人偷拍自拍在线| 中文字幕av特黄毛片| 欧美在线观看一区二区不卡| 丰满少妇高潮喷水视频| 91porny九色视频偷拍| 成人免费视频现网站99在线观看| 亚洲码av一区二区三区| 麻豆白洁少妇在线播放| 色哟哟亚洲乱码国产乱码精品精 | 日本不卡视频一二三区| 狠狠操av一区二区三区| 亚洲蜜桃久久久久久| 黄色片免费网站在线| 日韩av电影中文在线免费观看 | 国产精品国产三级在线高清观看| 亚洲无人区乱码中文字幕一区| 97人妻在线视频自拍| 国产农村乱子伦精精品视频| 中文字幕人妻一区二区视频系列| 欧美成人区一区二区三| 九十九步都是爱最后一步是尊严 | 最近最新最好看的中文字幕 | 天天插天天透天天爽| 亚洲无码专区中文字幕专区| 人妻被强av系列一区二区| 99999久久久精品| 99福利一区二区视频| —区二区三区女厕偷拍| 中文字幕日韩首页欧美在线激情| 538欧美在线观看一区二区三区 | 国产精品剧情在线亚洲| 91 精品视频在线看| 欧美一区二区三区视频看| 麻豆国产91制片厂| 乌克兰美女操逼高清内射视频| 伊人情人成综合视频| 97人妻av人人澡人人爽| 中国精品人妻一区二区| 国产av剧变态维修工虐杀美女| 男女69视频在线观看免费| 亚洲一区视频中文字幕在线播放| 9662av在线视频| 午夜亚洲国产精品中字| 青青青在线观看国产| 91久久久久久最新网站| 国产igao激情在线视频入口| 蜜桃tv一区二区三区| 国产精品亚洲精品亚洲| 首页欧美日韩中文字幕| 91精品国产综合99| 亚洲欧美成人激情在线| 久久热在线免费观看| 日本熟女0930视频| 人妻被强av系列一区二区| 2021国产剧情麻豆| 午夜宅男电影av网站| 久久av色噜噜ai换脸| 国产大桥未久一区二区| 久久99久久99久久97的人| 快使劲弄我视频在线播放| 中文字幕丰满子伦无码专区| 亚洲成人动漫av在线| 不用付费特黄特色亚洲特级黄色片| 久久精品国产亚洲av热软件| 一级毛片特级毛片免费的| 日日躁夜夜躁狠狠操| 天天插天天操天天射天天干| 国产精美视频精品视频精品 | 夫亡人妻被强干中文字幕| 日本韩国福利在线播放| 午夜在线观看一级毛| 午夜宅男电影av网站| 熟女俱乐部jukujoclub| 欧美日韩亚洲tv不卡久久| 免费在线小视频你懂的| 午夜一区二区三区视频在线观看| 国产剧情av在线免费观看| 中文字幕一区二区人妻视频| 国产成人综合久久婷婷| 在线有码人妻自拍视频| 亚洲国产精品一区51动漫| 亚洲一区二区精品三区视频| 一区二区三区内射美女| 亚洲成人,国产精品| 欧美一区二区三区爽爽| 欧美人与动欧交视频| 青青青青午夜手机国产视频| 日本一区二区三区调教性奴视频| 久久久精品人妻无码专区不卡| 久久99嫩草99久久精品| 久99久视频免费观看中文字幕| 亚成区一区二区人妻熟女| 美女张开腿给男人桶爽的软件| 天天操天天干加勒比久久| 日韩久久九九精品视频| 69精品互换人妻4p| 亚洲另类激情视频在线看| 呻吟求饶的人妻中文字幕| 97香蕉久久国产超碰| 精品美女洗澡一区二区| 亚洲综合一区二区三区四区| 人妻激情综合久久久久蜜桃| 天天插天天干天天狠| 亚洲一区二区三区无码在线| 天天爱天天日天天爽| 绿巨人浩克在线视频观看| 亚洲一区二区三区国产精品电影| 国产精品网站亚洲发布| 夜夜躁婷婷av蜜桃妖| 人妻激情偷乱一区二区三区av| 大乳丰满人妻中文字幕韩国hd| 婷婷六月天在线视频| 韩国毛片w妈妈的朋友7| 黑人大吊大战亚洲女人。| 国产91九色视频在线观看| 午夜呻吟亚洲精品中文字幕在上面| 久久av色噜噜ai换脸| 欧美vs亚洲vs日韩| 松本菜奈实最新av在线| 天天色 天天操 天天好逼| 亚洲一区亚洲二区成人福利| 开心五月综合激情婷婷| 天天干天天操天天日天天日| 亚洲av综合av一去二区三区| 50熟妇一区二区三区| 老熟女xxxⅹhd老熟女性| 夜夜骚av一二三区| 3344永久在线观看视频下载| 中文字幕日韩人妻在线三区| 国产精品视频网站污污污| 99久久国产精品免费消防器材| 成人人妻h在线观看| 人妻被强av系列一区二区| 国际精品熟女一区二区| 夜色17s精品人妻熟女av| 男女爱爱好爽视频免费看| 午夜精品久久久久久久精品乱码| 38av一区二区三区| 日本福利网站一区二区| 天天综合久久无人区| 96在线观看免费播放| 50熟妇一区二区三区| 欧美日韩亚洲国产视频二区| 欧美男女一区二区三区| 国产夫妻视频在线观看免费| 最近在线中文字幕免费| av无限看熟女人妻另类av| 91精品一区一区三区| 欧美人与动欧交视频| 男女真人做带声音视频图片| av丝袜免费在线观看| 天天操天天干加勒比久久| 91精品视频在线观看视频| 99女福利女女视频在线播放| 农村大炕有肉大屁股熟妇| 精品不卡一区二区三区| 91人妻人人做人人爽高清| 午夜精品一区二区三区不卡顿 | 中文在线字幕免费观看日韩视频| 天天看片天天摸天天操| avgo成人短视频| 日本一道中文字幕99| 成人做爰av在线观看网站| 人妻熟女 亚洲 一页二页| 日本熟女0930视频| 在线免费观看视频18| 久久久久久久久久久久久国产| 成人免费视频现网站99在线观看| av毛片在线观看网址| 午夜亚洲国产精品中字| 有码一区二区三区四区五区| 99热99这里免费的精品| 9662av在线视频| 天天日天天干天天日天天干天天| 亚洲春色av中文字幕|