24小时热门版块排行榜    

查看: 4756  |  回复: 10
本帖产生 1 个 1ST强帖 ,点击这里进行查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

yuhuashi1

铜虫 (初入文坛)

[求助] 请教有关DMOL3的设置问题

1. 用dmol3优化晶胞,想使晶胞的 xyz 轴夹角始终不变该怎么设置呢?help 里面说dmol不支持lattice parameters的固定啊!之前发现“use symmetry”打上勾可以保持xyz夹角不变,现在想请教大家如果不使用对称性,还有其他办法实现整个胞的夹角不变吗?

2. 关于dmol3加电场的问题。一般建立模型后直接run是不是默认没有电场?为什么和我 save files 之后加electric_field 0 0 0 算出来的总能不一样呢?想不通啊!

真诚请教有经验的虫友们,先谢谢大家了!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zjuer

至尊木虫 (著名写手)

【答案】应助回帖

★ ★ ★
mazuju028(金币+3, 1ST强帖+1): 谢谢交流 2011-07-03 22:49:35
yuhuashi1(金币+10): 非常感谢您的耐心回复!受益匪浅!谢谢~O(∩_∩)O~ 2011-07-04 11:02:15
Opt_Constraint
Purpose
The lines following the Opt_Constraint keyword define constrained internal coordinates (bonds, angles, torsions) in molecules or solids. These internals are held fixed at the values of the initial structure.

Syntax
Opt_Constraint
constraint specifications ...
blank line

Description
For molecular structures the format is:

Opt_Constraint
i j k l





where the integers i, j, k, and l correspond to the numbers of the atoms in the .car file. The following designations are used:

i j 0 0: constrain the bond between atoms i and j
i j k 0: constrain the angle formed by atoms i, j, and k
i j k l: constrain the torsion formed by all 4 atoms.
For periodic structures the format is:

i ia ib ic j ja jb jc k ka kb kc l la lb lc





where i,j,k,l have the same meaning as for molecules, and ia ib ic, etc. indicate which unit cell contains the atom.

These constraints are added automatically when you run via the Materials Studio user interface. Use the Constraints option on the Modify menu to define constraints via the interfaces.

Example
In a molecule to freeze the bond between atoms 5 and 6, the bond between atoms 7 and 8, and the torsion 4-5-6-7 the input is:

Opt_constraint
5  6  0  0
7  8  0  0
4  5  6  7


In a periodic structure, to constrain the bond between atoms 1 and 2 in the same unit cell:

Opt_constraint
1 0 0 0 2 0 0 0 0 0 0 0 0 0 0 0


In a periodic structure, to constrain the bond between atom 3 in one unit cell and the periodic image of atom 4 one cell in the "a" direction:

Opt_constraint
3 0 0 0 4 1 0 0 0 0 0 0 0 0 0 0


To include a second constraint in the above example use:

Opt_constraint
3  0 0 0 4 1 0 0 0 0 0 0 0 0 0 0
2 -1 0 0 3 0 0 0 4 0 1 0 0 0 0 0


This constrains the angle, including the periodic image of atom 2 in the cell in the "-a" direction, atom number 3 in the reference cell, and the periodic image of atom number 4 in the "+b" direction.

将构XYZ三边的四个原子的相互角度固定可能实现你的目的。

2、加电场需要手动加关键词
electric_field
Purpose
The electric_field keyword allows you to specify the value of the external electric field used in the calculation. The field is described in terms of three cartesian components.

Syntax
electric_field value_in_X value_in_Y value_in_Z

command keyword  values  default  meaning  
electric_field  real real real  0.0 0.0 0.0  The x, y and z components of the electric field applied to the system.  

Example
This example is for electric field with strength 0.05 a.u. applied in the Z direction only:

electric_field  0.000 0.000 0.0500

» 本帖已获得的红花(最新10朵)

4楼2011-07-03 22:36:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

uuv2010

荣誉版主 (职业作家)

优秀版主

★ ★
mazuju028(金币+2): 谢谢交流 2011-07-03 22:47:30
1. 默认就是不优化晶格常数以及夹角的
2. 不知道

» 本帖已获得的红花(最新10朵)

2楼2011-07-03 20:02:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuhuashi1

铜虫 (初入文坛)

送鲜花一朵
引用回帖:
Originally posted by uuv2010 at 2011-07-03 20:02:32:
1. 默认就是不优化晶格常数以及夹角的
2. 不知道

非常感谢!我计算时只要选中optimize cell,并且不加对称性,最后得到的晶格常数都发生了变化。不知道哪里出了问题
3楼2011-07-03 21:35:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuhuashi1

铜虫 (初入文坛)

送鲜花一朵
引用回帖:
Originally posted by zjuer at 2011-07-03 22:36:08:
Opt_Constraint
Purpose
The lines following the Opt_Constraint keyword define constrained internal coordinates (bonds, angles, torsions) in molecules or solids. These internals are held fixed at t ...

非常感谢您的耐心回复!对我挺有帮助!可还有些地方不太清晰:
第一个问题,最后您建议固定xyz三边的四个原子间的角度,很受启发。可我现在的模型是一个单层薄片,在Z方向(垂直薄片方向)加的真空层,该怎么保证Z轴垂直,还是有点迷惑。
第二个问题,我也是按您说的步骤手动加的关键词,问题是当我加0 0 0电场时和直接优化不加电场所得的体系总能不一致,本来0电场和不加电场应该效果一样才对啊。
再次表示感谢了!
5楼2011-07-04 10:57:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 PbwVEMK5haaE 2026-08-25 4/200 2026-08-26 03:20 by cNXvBfCpiZOM
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 Qj98d49dHFb9 2026-08-25 4/200 2026-08-26 03:08 by cNXvBfCpiZOM
[基金申请] 今天务委会开完了,明天出结果吗 +18 angus9576 2026-08-25 22/1100 2026-08-26 00:41 by merchancy
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 19/950 2026-08-25 17:57 by chick875
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 2026国自然函评费到账 +20 羊腰板 2026-08-21 23/1150 2026-08-25 16:34 by zsna
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 只有每年这种时候来逛逛小木虫 +25 yaoyewhu2008 2026-08-20 27/1350 2026-08-25 08:01 by Equinoxhua
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
信息提示
请填处理意见