应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 请教有关DMOL3的设置问题
111/2返回列表12下一页
查看: 3722  |  回复: 10
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 1ST强帖 ，点击这里进行查看

yuhuashi1

铜虫 (初入文坛)

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

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

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

真诚请教有经验的虫友们，先谢谢大家了！
回复此楼

» 收录本帖的淘帖专辑推荐

Materials Modeling Graphene DFT加电场

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2011-07-03 16:01:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)

优秀版主
★ ★
mazuju028(金币+2): 谢谢交流 2011-07-03 22:47:30
1. 默认就是不优化晶格常数以及夹角的
2. 不知道
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

yuhuashi1
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的回帖

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
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

yuhuashi1
4楼2011-07-03 22:36:08
 已阅   回复此楼   关注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的回帖

zjuer

至尊木虫 (著名写手)
★ ★
franch(金币+2): 谢谢回帖交流 2011-07-04 18:08:47
引用回帖:
Originally posted by yuhuashi1 at 2011-07-03 21:35:13:
非常感谢！我计算时只要选中optimize cell，并且不加对称性，最后得到的晶格常数都发生了变化。不知道哪里出了问题

选中optimize cell就是对晶格参数和原子位置同时进行优化，不选则只优化原子位置。
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

yuhuashi1
6楼2011-07-04 11:30:59
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuhuashi1

铜虫 (初入文坛)
送鲜花一朵
引用回帖:
Originally posted by zjuer at 2011-07-04 11:30:59:
选中optimize cell就是对晶格参数和原子位置同时进行优化，不选则只优化原子位置。

恩，是的！之前不选optimize cell计算了几次都失败了，我再调整下参数试试！
赞一下 回复此楼
7楼2011-07-04 11:47:54
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zjuer

至尊木虫 (著名写手)
★ ★
zzy870720z(金币+2): 谢谢指教 2011-07-06 14:10:35
第一个问题，最后您建议固定xyz三边的四个原子间的角度，很受启发。可我现在的模型是一个单层薄片，在Z方向（垂直薄片方向）加的真空层，该怎么保证Z轴垂直，还是有点迷惑。


只有一层的时候哪怕优化过程中C轴与A、B轴的夹角发生变化不再垂直，这个只是看起来不好看，但是对问题本身似乎是没影响的，因为你保证相邻层之间足够远，没有相互作用，所以C轴没什么意义，你可以计算完了在写论文作图的时候再把它做成垂直的。
赞一下(1人) 回复此楼
8楼2011-07-05 10:08:13
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

WJ-MS

新虫 (小有名气)
引用回帖:
1899033楼: 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 ...

想请教一下，modify-constraint 里面不是有一系列的约束选项么，能通过那个来固定吗？
赞一下 回复此楼
加油！
9楼2012-05-07 10:31:32
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaosong2011

银虫 (初入文坛)
我请教前辈们一个问题：我在计算场发射特性时，施加一个外加电场，单位是.a.u，那么如果我加（0， 0 ，0.37）的电场，就是相当于在Z轴施加了1eV/A的电场强度吗？
棘手的问题，恳请知情的前辈指教
赞一下 回复此楼
友爱互助，用爱相处。
10楼2012-06-18 10:24:03
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yuhuashi1 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭