计算化学公社

 找回密码 Forget password
 注册 Register
Views: 12782|回复 Reply: 12
打印 Print 上一主题 Last thread 下一主题 Next thread

[CP2K] 如何做单双轴应变加载呢

[复制链接 Copy URL]

148

帖子

0

威望

887

eV
积分
1035

Level 4 (黑子)

各位大佬,请问对于一个斜着六方晶格怎么实现单双轴的应变加载呢,例如如图所示的模型,单轴只加Armchair方向5%的应变,是不是把x方向格子长度a乘以1.05*cos30°呀。双轴加载5%的应变是不是同时再把斜着zigzag方向晶格长度b乘以1.05就行了。然后把改变长度的晶格参数固定来opt cell呀。活着有啥软件需要什么特殊计算嘛,欢迎大佬们指点!先手动@两位大佬 @卡开发发  @kyuu

202107041727303250..png (79.47 KB, 下载次数 Times of downloads: 52)

202107041727303250..png

3621

帖子

3

威望

1万

eV
积分
18429

Level 6 (一方通行)

第一原理惨品小作坊

2#
发表于 Post on 2021-7-4 20:16:53 | 只看该作者 Only view this author
抱歉,力学性质计算和cp2k的使用都不是我所擅长的,可以看看其他人有什么好的点子
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。
本周忙

148

帖子

0

威望

887

eV
积分
1035

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2021-7-4 21:13:11 | 只看该作者 Only view this author
卡开发发 发表于 2021-7-4 20:16
抱歉,力学性质计算和cp2k的使用都不是我所擅长的,可以看看其他人有什么好的点子。

嗯嗯,多谢大大捧场

353

帖子

0

威望

2029

eV
积分
2382

Level 5 (御坂)

真 掘墓者

4#
发表于 Post on 2021-7-8 19:19:42 | 只看该作者 Only view this author
两种方法 可以固定应变或固定应力
小体系固定应变就好了 大体系后者
然后看弹性形变变化

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
ycl0426 + 1

查看全部评分 View all ratings

圣诞刨坟忙

148

帖子

0

威望

887

eV
积分
1035

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2021-7-8 21:17:39 | 只看该作者 Only view this author
kyuu 发表于 2021-7-8 19:19
两种方法 可以固定应变或固定应力
小体系固定应变就好了 大体系后者
然后看弹性形变变化

狗爷,久仰久仰!您看咱这上面的应变加载没毛病吧;另外为啥大体系倾向于固定应力呢,cp2k里也能施加应力的嘛

353

帖子

0

威望

2029

eV
积分
2382

Level 5 (御坂)

真 掘墓者

6#
发表于 Post on 2021-7-18 16:27:07 | 只看该作者 Only view this author
Aridea 发表于 2021-7-8 21:17
狗爷,久仰久仰!您看咱这上面的应变加载没毛病吧;另外为啥大体系倾向于固定应力呢,cp2k里也能施加应力 ...

因为前者是定容体系, 大体系施加应变或应力或引起体积变化,所以做固定应力做定压计算,在两者皆可应用的范围内固定应变和应力结果应当是相同的
圣诞刨坟忙

148

帖子

0

威望

887

eV
积分
1035

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2021-7-18 19:42:12 | 只看该作者 Only view this author
kyuu 发表于 2021-7-18 16:27
因为前者是定容体系, 大体系施加应变或应力或引起体积变化,所以做固定应力做定压计算,在两者皆可应用 ...

狗爷,这个体积变化应该靠体积模量和泊松比来确定吧,对于二维材料(一二层原子这种)的体积变化是不是不用考虑呀;

353

帖子

0

威望

2029

eV
积分
2382

Level 5 (御坂)

真 掘墓者

8#
发表于 Post on 2021-7-21 15:05:36 | 只看该作者 Only view this author
Aridea 发表于 2021-7-18 19:42
狗爷,这个体积变化应该靠体积模量和泊松比来确定吧,对于二维材料(一二层原子这种)的体积变化是不是不 ...

体积模量的英文叫bulk modulus
and so.....
圣诞刨坟忙

148

帖子

0

威望

887

eV
积分
1035

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2021-7-21 15:07:48 | 只看该作者 Only view this author
kyuu 发表于 2021-7-21 15:05
体积模量的英文叫bulk modulus
and so.....

狗爷,您这也太忙了,听冰姐说如果我学QE,狗爷可以手把手教学,我现在怎么QE都安装不上呢

353

帖子

0

威望

2029

eV
积分
2382

Level 5 (御坂)

真 掘墓者

10#
发表于 Post on 2021-7-23 14:51:52 | 只看该作者 Only view this author
Aridea 发表于 2021-7-21 15:07
狗爷,您这也太忙了,听冰姐说如果我学QE,狗爷可以手把手教学,我现在怎么QE都安装不上呢

BBJ那么nb 让他手把手教你安装
圣诞刨坟忙

148

帖子

0

威望

887

eV
积分
1035

Level 4 (黑子)

11#
 楼主 Author| 发表于 Post on 2021-7-23 16:24:58 | 只看该作者 Only view this author
kyuu 发表于 2021-7-23 14:51
BBJ那么nb 让他手把手教你安装

哈哈,术业有专攻嘛,冰姐让我yum安装的

353

帖子

0

威望

2029

eV
积分
2382

Level 5 (御坂)

真 掘墓者

12#
发表于 Post on 2021-7-24 02:07:50 | 只看该作者 Only view this author
Aridea 发表于 2021-7-23 16:24
哈哈,术业有专攻嘛,冰姐让我yum安装的

足够用了
圣诞刨坟忙

1

帖子

0

威望

27

eV
积分
28

Level 2 能力者

13#
发表于 Post on 2024-7-30 18:17:53 | 只看该作者 Only view this author
本帖最后由 匹诺曹 于 2024-7-30 22:09 编辑

应该构建一个矩形晶胞,让矩形晶胞两个边对应扶手椅方向和之字形方向来施加应变。你那样按照cos30度算的施加应变的方法还是改变的Armchair方向,因为a方向和b方向根据对称性是一样的。

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-24 06:19 , Processed in 0.175624 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list