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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 【求助】castep计算弹性常数问题
51/1返回列表
查看: 5185  |  回复: 14
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

truemantwo

木虫 (正式写手)

[交流] 【求助】castep计算弹性常数问题

如题
用CASTEP模块计算材料的弹性常数时，必须对模型进行“impose symmetry”吗？如果没有施加对称性，对结果会有影响吗？我自己对比了一下两种操作的结果，大致差不多，只是有一点差别，但是“impose symmetry”后速度快很多，理论上，每个体系对应有自己的晶系和空间群，如果是自己手动建的模应该要施加对称性，计算结果才能体现体系的性质吧？欢迎大家讨论
回复此楼

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

仿真建模与计算 晶体模型计算分析

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2011-01-13 13:21:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sars518

木虫 (著名写手)
★ ★
zhang668(金币+2): 多谢交流 2011-03-05 13:09:53
castep 用的是 应力应变关系来计算cij

加入对称性后，速度更快，因为有些cij是对称的，这样减少了计算量

用P1的话，算起来不是一般的慢啊

此外，castep 计算cij，有的体系受cutoff和kpoint的影响非常大，所以计算前需要好好测试下
赞一下(2人) 回复此楼
7楼2011-01-18 18:34:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

cash_ms

木虫 (正式写手)
★ ★ ★
truemantwo(金币+2):O(∩_∩)O谢谢 2011-01-14 09:09:39
zzy870720z(金币+3):谢谢讨论交流 2011-01-14 09:39:37
引用回帖:
Originally posted by truemantwo at 2011-01-13 13:21:22:
如题
用CASTEP模块计算材料的弹性常数时，必须对模型进行“impose symmetry”吗？如果没有施加对称性，对结果会有影响吗？我自己对比了一下两种操作的结果，大致差不多，只是有一点差别，但是“impose symmetry” ...

对结构施加对称性，不是必要的，但是是需要的。
施加对称性之后，计算的时候因为某些对称性的结构的某些方向上Cij是相同的，或者可以由其他方向的Cij计算出，即不是独立的Cij，所以在计算的时候，会不去计算这一部分，所以计算的速度快了很多。
即使不施加对称性，但是这个结构的这种对称性还是存在的，所以算出来的结果差不多，某些方向的Cij应该是0值的，不施加对称性的话，会近似是0，可以不予考虑。
拙见，欢迎讨论。
赞一下(4人) 回复此楼
2楼2011-01-14 08:18:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

truemantwo

木虫 (正式写手)
这样啊，谢谢你的回答。
其他朋友的看法也一样吗？
赞一下 回复此楼
3楼2011-01-14 09:10:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yefeixjtu

金虫 (小有名气)
★
gzqdyouxia(金币+1):谢谢交流 2011-01-14 13:20:17
truemantwo(金币+1):谢谢 2011-01-14 13:21:12
不知道用的是应力应变关系计算，还是用应变能应变关系计算。如果是前者跟前述虫友差不多；后者的话，自己要对晶体矩阵施加应变模式，让体系具有一定的变形，这样也就改变了原有体系的空间群，可以用P1来算，但计算量比较大。
赞一下(2人) 回复此楼
4楼2011-01-14 13:01:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026国自然函评费到账 +15 羊腰板 2026-08-21 16/800 2026-08-23 10:45 by process2012
[基金申请] 让我中一个面上吧！ +12 大萍1987 2026-08-20 14/700 2026-08-23 10:30 by wrm
[基金申请] 今天放榜吗？ +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 07:41 by OEbVnUOu01ol
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[论文投稿] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 04:04 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:55 by OEbVnUOu01ol
[论文投稿] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 03:52 by OEbVnUOu01ol
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:43 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 03:40 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[基金申请] filecode，4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 今天基金会出结果吗？20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
关闭
请填处理意见
关闭