| 查看: 476 | 回复: 1 | |||
[交流]
晶格常数 不能输出，贴出输入的.py文件。大家看看哪里写错了
|
| 任务是：优化sic，以下是输入的.py文件内容，可以运行输出结果，但是，输出内容中，没有晶格常数 |
回复此楼» 猜你喜欢
300求调剂
已经有10人回复
-
26年电池方向博士申请
已经有5人回复
-
急需调剂
已经有6人回复
-
RY：中国产出的科学垃圾论文，绝对数量和比例都世界第一
已经有18人回复
-
一志愿中科大材料与化工，353分还有调剂学校吗
已经有7人回复
-
291求调剂
已经有8人回复
-
271求调剂
已经有36人回复
-
0854调剂
已经有9人回复
-
322求调剂
已经有15人回复
-
085404 298分求调剂
已经有13人回复
» 抢金币啦！回帖就可以得到:
-
26年二次电池储能方向招博士-中山大学四青人才课题组
+1/99
-
欧洲“CO2捕集与加氢催化转化”博士招聘（捷克-德国-斯洛文尼亚三国项目）
+1/92
-
伦敦布鲁内尔大学诚招国家公派读博，联合培养博士生，物流供应链方向
+1/84
-
郑州大学物理学院倪佩楠教授课题组诚邀博士后、青年教师及副教授加盟
+1/43
-
化学、化学工程、生物与医药，微生物学、生物化学与分子生物学专业调剂
+1/38
-
盐城师范学院生物与医药（生物技术与工程领域）还有调剂名额
+1/36
-
北华大学 08代码 调剂虚位以待
+1/36
-
上海工作，94年的博士
+1/35
-
武汉纺织大学学硕调剂
+1/19
-
吉林大学电合成与组装方向
+1/18
-
中科院合肥物质研究院招聘水系锌电博士后（长期有效）诚聘海内外优秀博士后
+1/18
-
哈尔滨工业大学（深圳）赵怡潞课题组诚招博士后、博士生
+1/10
-
中南林院士团队 招 生物质能源与材料 博士生
+1/9
-
佛山大学动物科技学院2026年硕士研究生招生调剂公告（第三批）
+1/8
-
西安工业大学光电学院接收调剂
+1/7
-
大连大学-贵州省煤炭洁净利用重点实验室联合培养研究生 化学5人+环境工程7人
+1/7
-
莆田学院环境与生物工程学院资源与环境专硕（环境工程方向） 开放调剂招生
+1/7
-
26申博，ai+材料方向，求课题组推荐
+1/3
-
以色列理工-生物质塑料等催化转化或多相流体力学方向---全奖博士研究生和科研助理
+2/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
