| 查看: 489 | 回复: 1 | |||
[交流]
晶格常数 不能输出，贴出输入的.py文件。大家看看哪里写错了
|
| 任务是：优化sic，以下是输入的.py文件内容，可以运行输出结果，但是，输出内容中，没有晶格常数 |
回复此楼» 猜你喜欢
青C资助名额大幅增加！
已经有13人回复
-
风把牡丹吹跑了
已经有7人回复
-
教学课件你会给同学吗
已经有8人回复
-
26应届毕业生考博求助
已经有3人回复
-
重磅！青年科学基金项目（C类）资助增幅预计超过50%
已经有7人回复
-
求助2，4-二氯-5-嘧啶甲醛的合成方法
已经有5人回复
-
本人最近太闲了，谁有问题可以提，每天会统一回复
已经有19人回复
-
申博自荐
已经有6人回复
-
这年头没有找到涵评专家，还有中面上的可能吗
已经有8人回复
-
材料类只有一篇综述能申博么
已经有4人回复
» 抢金币啦！回帖就可以得到:
-
南京大学化学学院有机学科徐骞课题组诚聘博士后及研究助理
+1/76
-
南京赫尔辛基大气与地球系统科学学院HALONEN ROOPE MIKKO SANTERI课题组招聘博士后
+1/61
-
华南师范大学（211）- 生物光子学研究院招收生物学、物理学、神经科学等专业博士
+1/31
-
谷歌学术引用量
+1/29
-
双一流高校-南京林业大学-化学工程学院-国家海外优青团队招2026级博士（5月15号截止）
+1/28
-
清华大学机械系吴爱萍教授课题组招收博士后
+1/9
-
双一流高校-南京林业大学-化学工程学院-国家海外优青团队招2026级博士（5月15号截止）
+1/5
-
南方科技大学物理系 CSLM Lab 招生招聘 （软物质/生物物理，Physics of AI）
+1/5
-
检测XPS/XRD/SEM/拉曼/ICPQQ:3392739119
+1/5
-
华侨大学（厦门）——急招化工/材料/化学/能动——26年9月份博士
+1/5
-
华南理工大学微电子学院招电子信息专业博士研究生（2026年秋季入学）
+1/5
-
ChemicalBook英文词条兼职招聘
+1/5
-
【招收2026年秋博士研究生】哈工大（深圳）医工学院任卫红老师（15号截止！！！）
+1/5
-
中山大学招收2027年秋季入学的博士研究生（计算机视觉/机器学习）
+1/5
-
安徽农业大学植物保护博士研究生招生（报名截止到 5月20日）-生物或者化学方向
+1/4
-
海南大学黄玲/徐从军团队2026年博士研究生
+1/4
-
双一流高校-南京林业大学-化学工程学院-国家海外优青团队招青年师资（正式事业编制）
+1/3
-
【26年9月入学；博士招募】海南大学玄萍教授课题组（人工智能与信息处理方向）
+1/3
-
试遍所有AI写作工具，最终锁定这个组合 Prose + DeepSeek V4 Pro
+1/1
-
ACS inorganic chemistry 拒稿重投，editorial review一般多久会下一步
+1/1
|
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
精华III: