| 查看: 469 | 回复: 1 | |||
[交流]
晶格常数 不能输出，贴出输入的.py文件。大家看看哪里写错了
|
| 任务是：优化sic，以下是输入的.py文件内容，可以运行输出结果，但是，输出内容中，没有晶格常数 |
回复此楼» 猜你喜欢
材料调剂
已经有7人回复
-
化学工程085602 305分求调剂
已经有10人回复
-
289求调剂
已经有15人回复
-
291 求调剂
已经有7人回复
-
274求调剂
已经有14人回复
-
北京林业大学硕导招生广告
已经有3人回复
-
309求调剂
已经有5人回复
-
292求调剂
已经有8人回复
-
求调剂
已经有4人回复
-
一志愿 西北大学 总分282 英语一62 求调剂
已经有3人回复
» 抢金币啦！回帖就可以得到:
-
中国散裂中子源诚聘“AI材料计算”方向博士后/职工
+1/415
-
西南交通大学医学院招收生物医学工程和材料与化工专业硕士研究生
+1/85
-
太原科技大学材料科学与工程学院招收金属材料类博士一名
+1/85
-
招收调剂生，动物医学专业
+1/84
-
中国石油大学（北京）国家级大人才团队博士招生1名：化学、材料、石油工程：油田化学
+1/38
-
哈尔滨工业大学（深圳）-何自开教授团队诚招化学工程与技术专业博士研究生
+1/31
-
同济大学环境学院 肖倩特聘研究员课题组 招聘硕士/博士（长期有效）
+1/17
-
双一流大学湘潭大学“化工过程模拟与强化”国家地方联合工程研究中心招收各类博士生
+1/17
-
【实战型】【生物医药】2026青岛大学招博士生 含少数民族骨干计划2名！
+1/17
-
上海工程技术大学纺织服装学院辛斌杰院长、郑元生教授团队招收学硕/专硕研究生
+1/9
-
欢迎来英国拉夫堡大学交流-结晶，蛋白质，AI，水处理等
+1/8
-
各种IgG抗体表达系统
+1/6
-
山东济南某省属重点高校-环境生态类专业有调剂名额
+1/6
-
化学工程与技术 国家双一流学科 211高校 国家重点实验室 博士研究生
+1/5
-
一志愿西南交大机械专硕343求调剂
+1/4
-
上海师范大学有机化学专业研究生招生
+1/3
-
26年招收07 08学术型/专业型硕士研究生
+1/3
-
山东重点高校招收08开头工科调剂考生
+1/3
-
Top-88悉尼科技大学数据科学/AI 招收2027年入学 校奖 博士生1到2名(国际和本地学生）
+1/2
-
中国人民大学物理学系蔡鹏课题组招收优秀硕士博士研究生
+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
