24小时热门版块排行榜    

查看: 1161  |  回复: 0

赵奇一

新虫 (小有名气)

[求助] QE的一个例子不是太理解,求解释

在研究例子中对这两个过程的区别不是太清楚,求解释一下:
# self-consistent calculation
cat > al001.rx.in << EOF
&CONTROL
  calculation = "relax",
  pseudo_dir  = "$PSEUDO_DIR",
  outdir      = "$TMP_DIR",
/
&SYSTEM
  ibrav       = 6,
  celldm(1)   = 5.3033D0,
  celldm(3)   = 8.D0,
  nat         = 7,
  ntyp        = 1,
  ecutwfc     = 12.D0,
  occupations = "smearing",
  smearing    = "methfessel-paxton",
  degauss     = 0.05D0,
/
&ELECTRONS
  conv_thr    = 1.D-6,
  mixing_beta = 0.3D0,
/
&IONS
  bfgs_ndim         = 3,
/
ATOMIC_SPECIES
Al  1.0  Al.pz-vbc.UPF
ATOMIC_POSITIONS
Al   0.5000000      0.5000000     -2.121320
Al   0.0000000      0.0000000     -1.414213
Al   0.5000000      0.5000000     -0.707107
Al   0.0000000      0.0000000      0.000000
Al   0.5000000      0.5000000      0.707107
Al   0.0000000      0.0000000      1.414213
Al   0.5000000      0.5000000      2.121320
K_POINTS
3
0.125  0.125  0.0  1.0
0.125  0.375  0.0  2.0
0.375  0.375  0.0  1.0
关于上边,说明上是说:make a geometry relaxation for a Al (001) slab performing a series of  self-consistent calculations and computing the forces on atoms。
下边这个计算,说明上是这么描述的:make a geometry relaxation for the same Al (001) slab used in step 2  performing a series of self-consistent calculations, computing the   forces on atoms and evolving the atomic positions according to Newton   equation.
# self-consistent calculation
rm -f e eal ave p avec tv
cat > al001.mm.in << EOF
&CONTROL
  calculation = "relax",
  dt          = 30.D0,
  pseudo_dir  = "$PSEUDO_DIR",
  outdir      = "$TMP_DIR",
/
&SYSTEM
  ibrav       = 6,
  celldm(1)   = 5.3033D0,
  celldm(3)   = 8.D0,
  nat         = 7,
  ntyp        = 1,
  ecutwfc     = 12.D0,
  occupations = "smearing",
  smearing    = "methfessel-paxton",
  degauss     = 0.05D0,
/
&ELECTRONS
  conv_thr    = 1.D-7,
  mixing_beta = 0.3D0,
/
&IONS
  ion_dynamics      = "damp",
  pot_extrapolation = "second_order",
  wfc_extrapolation = "second_order",
/
ATOMIC_SPECIES
Al  1.D0  Al.pz-vbc.UPF
ATOMIC_POSITIONS
Al   0.5000000      0.5000000     -2.121320
Al   0.0000000      0.0000000     -1.414213
Al   0.5000000      0.5000000     -0.707107
Al   0.0000000      0.0000000      0.000000
Al   0.5000000      0.5000000      0.707107
Al   0.0000000      0.0000000      1.414213
Al   0.5000000      0.5000000      2.121320
K_POINTS
3
0.125  0.125  0.0  1.0
0.125  0.375  0.0  2.0
0.375  0.375  0.0  1.0
对于两个的区别不是太能理解,希望大神能尽可能详细的解答一下...
回复此楼
天使可以飞起来是因为他们把自己看得很轻
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 赵奇一 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 19:13 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 18:53 by hZiFeudZyoGR
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 18:51 by hZiFeudZyoGR
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-30 18:29 by hZiFeudZyoGR
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 gy1nBQXYQJqL 2026-08-29 9/450 2026-08-30 17:56 by hZiFeudZyoGR
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 17:14 by hZiFeudZyoGR
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +8 ASdOkHsho7FD 2026-08-28 12/600 2026-08-30 17:11 by hZiFeudZyoGR
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 14:05 by GkCVyODlx4mR
[基金申请] 有没有仍没收到信息的 +6 德尚中行 2026-08-27 7/350 2026-08-30 12:53 by Bxndjjd
[基金申请] 国社科又开始会评了,不知道这次命运如何 +3 雨打竹帘 2026-08-30 4/200 2026-08-30 12:35 by 白大多亮
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 11:51 by l0VvVHGBGRLv
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 为什么到现在没收到通知? +4 tannykie 2026-08-29 4/200 2026-08-29 17:14 by lmz0216
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
请填处理意见