| 查看: 556 | 回复: 1 | |||
[交流]
晶格常数 不能输出，贴出输入的.py文件。大家看看哪里写错了
|
| 任务是：优化sic，以下是输入的.py文件内容，可以运行输出结果，但是，输出内容中，没有晶格常数 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有7人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有7人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有10人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
» 抢金币啦！回帖就可以得到:
-
坐标苏州，诚聘女友
+1/278
-
大连理工大学化工学院郑跃楠副教授计划招收2027级博士研究生
+1/192
-
大连理工大学机械工程学院招收微纳及原子制造方向博士生
+1/188
-
云南大学丁军桥团队2027年博士招生
+1/95
-
海南大学徐月山老师招收2027年第一批博士研究生1名，欢迎联系咨询。
+1/47
-
中国海洋大学 物理学博士点 固体发光方向 招收申请-考核制博士生1–2名
+2/44
-
中国海洋大学 物理学博士点 凝聚态物理-固体发光方向 招收申请-考核制博士生1–2名
+2/44
-
哈尔滨工业大学招收博士研究生（欢迎环境、生物、市政、农业、化学等专业）长期有效
+1/39
-
北京化工大学化学学院优青招收硕博士生
+1/18
-
27申博自荐，储能电池方向
+1/16
-
深圳理工大学核酸药物递送课题组博士后招聘启事
+1/16
-
南科大陈鹏课题组招收2027级博士生
+1/12
-
华南理工大学黄佳豪课题组招聘-铁电压电聚合物、电介质聚合物、智能柔性器件
+1/7
-
您身边的科研非标定制专家--赛思睿
+1/4
-
【2027博士招生】福州大学姚建年院士课题组招收博士研究生2名
+1/3
-
27申博 寻膜方向博导
+1/3
-
英国杜伦大学Dr Keith Andrews课题组(化学催化)，招收2027年CSC公派博士学生。
+1/2
-
中国药科大学陈达威课题组诚聘师资博士后(2023级博后留校率65%)
+1/2
-
【博士招生】湖南科技大学 陈建 教授 课题组招收博士研究生
+1/2
-
请问手中有可以提取3微米，6微米的球形颗粒，提取设备，这个技术有市场吗？
+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
