| 查看: 474 | 回复: 1 | |||
[交流]
晶格常数 不能输出，贴出输入的.py文件。大家看看哪里写错了
|
| 任务是：优化sic，以下是输入的.py文件内容，可以运行输出结果，但是，输出内容中，没有晶格常数 |
回复此楼» 猜你喜欢
327求调剂
已经有28人回复
-
一志愿中科大材料与化工，353分还有调剂学校吗
已经有6人回复
-
0854求调剂
已经有17人回复
-
售SCI一区文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
药学求调剂
已经有15人回复
-
初试324 中药学 一志愿天中医 求调剂
已经有4人回复
-
279学硕食品专业求调剂院校
已经有29人回复
-
294求调剂
已经有10人回复
-
恳请有学校收留
已经有6人回复
-
求调剂推荐
已经有8人回复
» 抢金币啦！回帖就可以得到:
-
物理学调剂
+1/88
-
教程+代算 GROMACS模拟蛋白质-小分子、药物自组装、高分子弛豫、分子转运等
+1/88
-
物理学0702 调剂
+1/88
-
厦门华侨大学就读的泰国排名第一世界排名221位朱拉隆功大学学位正式招生了
+1/88
-
湖北大学材料学院 矿物与固废资源化利用环境能源新材料研究团队 第二轮调剂 少量名额
+1/86
-
湖北文理学院，材料，化学，环境，机械，电子，工程等专业均可报名调剂！
+3/86
-
新加坡国立大学化学系刘小钢教授课题组招收博士生3-4名
+1/85
-
中国林科院博士招生/读博
+2/80
-
上海应用技术大学刘敏课题组招生---欢迎材料，物理和化学相关背景
+1/41
-
平顶山学院 化工材料专业研究生报考调剂名额充足
+1/40
-
2026年盐城工学院化学工程与技术专业（学硕）招收08所有专业和0703方向考生
+1/39
-
湖南理工大学化学化工学院的化学和化学工程与技术专业调剂
+1/10
-
重庆三峡科技大学物理电子学、场波专业接收调剂，考数二可调，系统开至4月10日
+1/9
-
接收资源与环境专业硕士调剂
+1/8
-
材料与化工专业还有大量指标，欢迎大家调剂，调剂系统已经开放
+1/7
-
与其发展宠物经济，不如改变年轻人的婚恋观
+1/5
-
延安大学化学与化工学院接收化学专业调剂生，11号下午复试（工科、理学均可调剂）
+5/5
-
2026年北京服装学院- 材料学院还有调剂名额，可招收化学材料化工相关专业考研调剂
+1/3
-
中国海洋大学工程学院科研博士后招聘启事（数学/能源化工方向）
+1/2
-
西安工程大学环境与化学工程学院研究生调剂（理学，农学可调，专硕可调学硕）
+1/2
|
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
