| 查看: 468 | 回复: 1 | |||
[交流]
晶格常数 不能输出，贴出输入的.py文件。大家看看哪里写错了
|
| 任务是：优化sic，以下是输入的.py文件内容，可以运行输出结果，但是，输出内容中，没有晶格常数 |
回复此楼» 猜你喜欢
材料调剂
已经有7人回复
-
化学工程085602 305分求调剂
已经有10人回复
-
289求调剂
已经有15人回复
-
291 求调剂
已经有7人回复
-
274求调剂
已经有14人回复
-
北京林业大学硕导招生广告
已经有3人回复
-
309求调剂
已经有5人回复
-
292求调剂
已经有8人回复
-
求调剂
已经有4人回复
-
一志愿 西北大学 总分282 英语一62 求调剂
已经有3人回复
» 抢金币啦！回帖就可以得到:
-
山东建筑大学 （市政与环境工程学院）硕士研究生招生啦 静电纺丝纤维与盐水分离应用
+1/90
-
Postdoc position in Hangzhou Normal University （陈健翔组博后招聘）
+1/83
-
华南理工-广中医联合招收生物医学与材料交叉学科博士生
+1/81
-
福建农林大学材料工程学院高分子材料课题组招生
+1/79
-
南方油校——资源与环境专硕招生，化学环境生物类专业均可，就业好，有专业壁垒
+1/48
-
哈尔滨工业大学（深圳）-何自开教授团队诚招化学工程与技术专业博士研究生
+1/31
-
欢迎08或部分07开头的有缘人报考
+1/31
-
天津理工大学国家杰青王铁团队招聘
+1/30
-
上海中医药大学创新中药研究院 招收审核制博士生一名
+1/19
-
下周见一个合作者，应该问些什么
+1/16
-
【博士招生】广东工业大学国家优青课题组招收2026年环境方向博士生
+1/16
-
北京某研究院结构生物学相关专业接收调剂研究生
+1/10
-
西南某211招27届博士
+1/9
-
shu湘潭大学化工学院招收化工、能源动力、工程热物理相关研究生调剂
+1/8
-
南京林业大学化学工程学院柏惺峰课题组诚聘科研助理、水杉英才E类/F类青年教师或博后
+1/7
-
宁波东方理工大学余鹏课题组2026年上海交大联培博士生、博士后、科研助理招聘
+1/6
-
一志愿中国海洋大学095137，初试成绩312，不挑学校和地区，有学上就行！
+1/6
-
计算机(11408)考研361分求调剂
+1/5
-
2026调剂硕士生招生—北京双一流学科高校环境科学与工程专业
+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
