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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 晶格常数 不能输出，贴出输入的.py文件。大家看看哪里写错了
21/1返回列表
查看: 470  |  回复: 1
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

tiantian2921

新虫 (小有名气)

[交流] 晶格常数 不能输出，贴出输入的.py文件。大家看看哪里写错了

任务是：优化sic，以下是输入的.py文件内容，可以运行输出结果，但是，输出内容中，没有晶格常数
回复此楼

» 猜你喜欢

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

查看全部散金贴
1楼 2011-05-05 19:20:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiantian2921

新虫 (小有名气)
引用回帖:
Originally posted by tiantian2921 at 2011-05-05 19:20:28:
任务是：优化sic，以下是输入的.py文件内容，可以运行输出结果，但是，输出内容中，没有晶格常数

from ATK.KohnSham import *
# Specify bulk configuration with 4 silicon atoms and 4 carbon stoms
# as basis in a Hex Bravais lattice
SiC_Hex = BulkConfiguration(
   bravais_lattice = Hexagonal(a=3.078*Angstrom,c=10.046*Angstrom),        
    elements = 4*[Carbon]+4*[Silicon],
    fractional_coordinates = [
        [0.00, 0.00, 0.188],
        [0.00, 0.00, 0.688],
        [0.667, 0.333, 0.439],
        [-0.667, -0.333, 0.939],
        [0.00, 0.00, 0.00],
        [0.00, 0.00, 0.50],
        [0.667, 0.333, 0.251],
        [-0.667, -0.333, 0.751]
        ]
        )
# Display Cartesian coordinates for SiC crystal
#from printBulkConfig import printBulkConfig
#printBulkConfig(SiC_Hex)

# Store crystal structure in a VNL file for later use
vnl_file=VNLFile("SiC.vnl"
vnl_file.addToSample(SiC_Hex,"SiC_Hex"

#optimaize the  crystal

import ATK                          
verbosity_level=ATK.verbosityLevel()
ATK.setVerbosityLevel(2)            
                                                                                                                 
optimized_configuration = calculateOptimizedAtomicGeometry(                                                      
    atomic_configuration =SiC_Hex,                                                                                   
    method=KohnShamMethod(                                                                                   
        [basisSetParameters(type=SingleZeta,element=Carbon),basisSetParameters(type=SingleZeta,element=Silicon)],
        exchange_correlation_type=GGA.PBE,                                                                       
        electron_density_parameters=electronDensityParameters(mesh_cutoff=250.0*Rydberg),                        
        eigenstate_occupation_parameters=eigenstateOccupationParameters(                                         
        temperature=300*Kelvin),                                                                                 
        brillouin_zone_integration_parameters=brillouinZoneIntegrationParameters(monkhorst_pack_parameters=[6,6,2])                                    
        ),                                                                                                      
    optimization_parameters=geometricOptimizationParameters(                                                                  
        force_tolerance=0.001*eV/Ang,                                                                           
        max_steps = 300,                                                                                         
        optimizer = Optimizer.SteepestDescent,
        time_step = 0.5*femtosecond)                                                                           
        )

#print the lattice constant
                                                   
lattice=optimized_configuration                                                     
print "a = %g Bohr" % (lattice.getA().inUnitsOf(Bohr))
print "c = %g Bohr" % (lattice.getC().inUnitsOf(Bohr))
赞一下 回复此楼
2楼2011-05-05 19:20:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tiantian2921 的主题更新
21/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿 西北大学 总分282 英语一62 求调剂 +3 18419759900 2026-03-25 3/150 2026-03-25 23:20 by peike
[考研] 0856求调剂 +5 zhn03 2026-03-25 6/300 2026-03-25 22:58 by peike
[考研] 求调剂 +3 QiMing7 2026-03-25 3/150 2026-03-25 21:13 by 给你你注意休息
[考研] 材料与化工 322求调剂 +6 然11 2026-03-19 6/300 2026-03-25 18:37 by haxia
[考研] 招08考数学 +8 laoshidan 2026-03-20 17/850 2026-03-25 17:52 by 一个红太阳
[考研] 293求调剂 +7 加一一九 2026-03-24 7/350 2026-03-25 12:02 by userper
[考研] 0854电子信息求调剂 324 +4 Promise-jyl 2026-03-23 4/200 2026-03-25 11:36 by Sugarlight
[考研] 一志愿武理085500机械专业总分300求调剂 +3 an10101 2026-03-24 7/350 2026-03-25 00:00 by 山鬼0-
[考研] 0856材料专硕353求调剂 +6 NIFFFfff 2026-03-20 6/300 2026-03-24 21:02 by hello七七
[考研] 调剂 +4 13853210211 2026-03-24 4/200 2026-03-24 19:44 by ms629
[考研] 求调剂 +6 研研，接电话 2026-03-24 7/350 2026-03-24 17:01 by barlinike
[基金申请] 请教下大家 2026年国家基金申请是双盲审吗？ +3 lishucheng1 2026-03-22 5/250 2026-03-24 08:22 by gltch
[考研] 一志愿武理材料工程348求调剂 +6 ￣^￣゜汗 2026-03-19 9/450 2026-03-23 19:53 by pswait
[考研] 275求调剂 +6 shansx 2026-03-22 8/400 2026-03-22 15:27 by barlinike
[基金申请] 山东省面上项目限额评审 +4 石瑞0426 2026-03-19 4/200 2026-03-22 08:50 by Wei_ren
[考研] 一志愿苏州大学材料求调剂，总分315（英一） +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
[考研] 求调剂一志愿南京航空航天大学289分 +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
[考研] 一志愿西安交通大学 学硕 354求调剂211或者双一流 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
信息提示
关闭
请填处理意见
关闭