| 查看: 539 | 回复: 1 | |||
[交流]
晶格常数 不能输出，贴出输入的.py文件。大家看看哪里写错了
|
| 任务是：优化sic，以下是输入的.py文件内容，可以运行输出结果，但是，输出内容中，没有晶格常数 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有6人回复
-
面上没中，邀请各位路过的虫友分析一下分数
已经有7人回复
-
哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生
已经有3人回复
-
找导师
已经有8人回复
-
国社科又开始会评了，不知道这次命运如何
已经有14人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有13人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有14人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有17人回复
-
学科评审组评审是指会评吗？
已经有5人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
|
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